Identification of low-frequency variants associated with gout and serum uric acid levels

We tested 16 million SNPs, identified through whole-genome sequencing of 457 Icelanders, for association with gout and serum uric acid levels. Genotypes were imputed into 41,675 chip-genotyped Icelanders and their relatives, for effective sample sizes of 968 individuals with gout and 15,506 individuals for whom serum uric acid measurements were available. We identified a low-frequency missense variant (c.1580C>G) in ALDH16A1 associated with gout (OR = 3.12, P = 1.5 × 10(-16), at-risk allele frequency = 0.019) and serum uric acid levels (effect = 0.36 s.d., P = 4.5 × 10(-21)). We confirmed the association with gout by performing Sanger sequencing on 6,017 Icelanders. The association with gout was stronger in males relative to females. We also found a second variant on chromosome 1 associated with gout (OR = 1.92, P = 0.046, at-risk allele frequency = 0.986) and serum uric acid levels (effect = 0.48 s.d., P = 4.5 × 10(-16)). This variant is close to a common variant previously associated with serum uric acid levels. This work illustrates how whole-genome sequencing data allow the detection of associations between low-frequency variants and complex traits.     

Genome-wide association study identifies a new melanoma susceptibility locus at 1q21.3

We performed a genome-wide association study of melanoma in a discovery cohort of 2,168 Australian individuals with melanoma and 4,387 control individuals. In this discovery phase, we confirm several previously characterized melanoma-associated loci at MC1R, ASIP and MTAP-CDKN2A. We selected variants at nine loci for replication in three independent case-control studies (Europe: 2,804 subjects with melanoma, 7,618 control subjects; United States 1: 1,804 subjects with melanoma, 1,026 control subjects; United States 2: 585 subjects with melanoma, 6,500 control subjects). The combined meta-analysis of all case-control studies identified a new susceptibility locus at 1q21.3 (rs7412746, P = 9.0 × 10(-11), OR in combined replication cohorts of 0.89 (95% CI 0.85-0.95)). We also show evidence suggesting that melanoma associates with 1q42.12 (rs3219090, P = 9.3 × 10(-8)). The associated variants at the 1q21.3 locus span a region with ten genes, and plausible candidate genes for melanoma susceptibility include ARNT and SETDB1. Variants at the 1q21.3 locus do not seem to be associated with human pigmentation or measures of nevus density.     

A transition zone complex regulates mammalian ciliogenesis and ciliary membrane composition

Mutations affecting ciliary components cause ciliopathies. As described here, we investigated Tectonic1 (Tctn1), a regulator of mouse Hedgehog signaling, and found that it is essential for ciliogenesis in some, but not all, tissues. Cell types that do not require Tctn1 for ciliogenesis require it to localize select membrane-associated proteins to the cilium, including Arl13b, AC3, Smoothened and Pkd2. Tctn1 forms a complex with multiple ciliopathy proteins associated with Meckel and Joubert syndromes, including Mks1, Tmem216, Tmem67, Cep290, B9d1, Tctn2 and Cc2d2a. Components of this complex co-localize at the transition zone, a region between the basal body and ciliary axoneme. Like Tctn1, loss of Tctn2, Tmem67 or Cc2d2a causes tissue-specific defects in ciliogenesis and ciliary membrane composition. Consistent with a shared function for complex components, we identified a mutation in TCTN1 that causes Joubert syndrome. Thus, a transition zone complex of Meckel and Joubert syndrome proteins regulates ciliary assembly and trafficking, suggesting that transition zone dysfunction is the cause of these ciliopathies.     

The genome of the extremophile crucifer Thellungiella parvula

Thellungiella parvula is related to Arabidopsis thaliana and is endemic to saline, resource-poor habitats, making it a model for the evolution of plant adaptation to extreme environments. Here we present the draft genome for this extremophile species. Exclusively by next generation sequencing, we obtained the de novo assembled genome in 1,496 gap-free contigs, closely approximating the estimated genome size of 140 Mb. We anchored these contigs to seven pseudo chromosomes without the use of maps. We show that short reads can be assembled to a near-complete chromosome level for a eukaryotic species lacking prior genetic information. The sequence identifies a number of tandem duplications that, by the nature of the duplicated genes, suggest a possible basis for T. parvula's extremophile lifestyle. Our results provide essential background for developing genomically influenced testable hypotheses for the evolution of environmental stress tolerance.     

Dense genotyping identifies and localizes multiple common and rare variant association signals in celiac disease

Using variants from the 1000 Genomes Project pilot European CEU dataset and data from additional resequencing studies, we densely genotyped 183 non-HLA risk loci previously associated with immune-mediated diseases in 12,041 individuals with celiac disease (cases) and 12,228 controls. We identified 13 new celiac disease risk loci reaching genome-wide significance, bringing the number of known loci (including the HLA locus) to 40. We found multiple independent association signals at over one-third of these loci, a finding that is attributable to a combination of common, low-frequency and rare genetic variants. Compared to previously available data such as those from HapMap3, our dense genotyping in a large sample collection provided a higher resolution of the pattern of linkage disequilibrium and suggested localization of many signals to finer scale regions. In particular, 29 of the 54 fine-mapped signals seemed to be localized to single genes and, in some instances, to gene regulatory elements. Altogether, we define the complex genetic architecture of the risk regions of and refine the risk signals for celiac disease, providing the next step toward uncovering the causal mechanisms of the disease.     

A genetic interaction network of five genes for human polycystic kidney and liver diseases defines polycystin-1 as the central determinant of cyst formation

Autosomal dominant polycystic liver disease results from mutations in PRKCSH or SEC63. The respective gene products, glucosidase IIβ and SEC63p, function in protein translocation and quality control pathways in the endoplasmic reticulum. Here we show that glucosidase IIβ and Sec63p are required in mice for adequate expression of a functional complex of the polycystic kidney disease gene products, polycystin-1 and polycystin-2. We find that polycystin-1 is the rate-limiting component of this complex and that there is a dose-response relationship between cystic dilation and levels of functional polycystin-1 following mutation of Prkcsh or Sec63. Reduced expression of polycystin-1 also serves to sensitize the kidney to cyst formation resulting from mutations in Pkhd1, the recessive polycystic kidney disease gene. Finally, we show that proteasome inhibition increases steady-state levels of polycystin-1 in cells lacking glucosidase IIβ and that treatment with a proteasome inhibitor reduces cystic disease in orthologous gene models of human autosomal dominant polycystic liver disease.     

Genome-wide association study reveals three susceptibility loci for common migraine in the general population

Migraine is a common, heterogeneous and heritable neurological disorder. Its pathophysiology is incompletely understood, and its genetic influences at the population level are unknown. In a population-based genome-wide analysis including 5,122 migraineurs and 18,108 non-migraineurs, rs2651899 (1p36.32, PRDM16), rs10166942 (2q37.1, TRPM8) and rs11172113 (12q13.3, LRP1) were among the top seven associations (P < 5 × 10(-6)) with migraine. These SNPs were significant in a meta-analysis among three replication cohorts and met genome-wide significance in a meta-analysis combining the discovery and replication cohorts (rs2651899, odds ratio (OR) = 1.11, P = 3.8 × 10(-9); rs10166942, OR = 0.85, P = 5.5 × 10(-12); and rs11172113, OR = 0.90, P = 4.3 × 10(-9)). The associations at rs2651899 and rs10166942 were specific for migraine compared with non-migraine headache. None of the three SNP associations was preferential for migraine with aura or without aura, nor were any associations specific for migraine features. TRPM8 has been the focus of neuropathic pain models, whereas LRP1 modulates neuronal glutamate signaling, plausibly linking both genes to migraine pathophysiology.     

Variation in the DEPDC5 locus is associated with progression to hepatocellular carcinoma in chronic hepatitis C virus carriers

Chronic viral hepatitis is the most important risk factor for progression to hepatocellular carcinoma (HCC). To identify genetic risk factors for progression to HCC in individuals with chronic hepatitis C virus (HCV), we analyzed 467,538 SNPs in 212 Japanese individuals with chronic HCV with HCC and 765 individuals with chronic HCV without HCC. We identified one intronic SNP in the DEPDC5 locus on chromosome 22 associated with HCC risk and confirmed the association using an independent case-control population (710 cases and 1,625 controls). The association was highly significant when we analyzed the stages separately as well as together (rs1012068, P(combined) = 1.27 × 10(-13), odds ratio = 1.75). The significance level of the association further increased after adjustment for gender, age and platelet count (P = 1.35 × 10(-14), odds ratio = 1.96). Our findings suggest that common variants within the DEPDC5 locus affect susceptibility to HCC in Japanese individuals with chronic HCV infection.