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.     

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.     

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.     

Complex interactions between genes controlling trafficking in primary cilia

Cilia-associated human genetic disorders are striking in the diversity of their abnormalities and their complex inheritance. Inactivation of the retrograde ciliary motor by mutations in DYNC2H1 causes skeletal dysplasias that have strongly variable expressivity. Here we define previously unknown genetic relationships between Dync2h1 and other genes required for ciliary trafficking. Mutations in mouse Dync2h1 disrupt cilia structure, block Sonic hedgehog signaling and cause midgestation lethality. Heterozygosity for Ift172, a gene required for anterograde ciliary trafficking, suppresses cilia phenotypes, Sonic hedgehog signaling defects and early lethality of Dync2h1 homozygotes. Ift122, like Dync2h1, is required for retrograde ciliary trafficking, but reduction of Ift122 gene dosage also suppresses the Dync2h1 phenotype. These genetic interactions illustrate the cell biology underlying ciliopathies and argue that mutations in intraflagellar transport genes cause their phenotypes because of their roles in cilia architecture rather than direct roles in signaling.     

miRNA regulation of Sdf1 chemokine signaling provides genetic robustness to germ cell migration

microRNAs (miRNAs) function as genetic rheostats to control gene output. Based on their role as modulators, it has been postulated that miRNAs canalize development and provide genetic robustness. Here, we uncover a previously unidentified regulatory layer of chemokine signaling by miRNAs that confers genetic robustness on primordial germ cell (PGC) migration. In zebrafish, PGCs are guided to the gonad by the ligand Sdf1a, which is regulated by the sequestration receptor Cxcr7b. We find that miR-430 regulates sdf1a and cxcr7 mRNAs. Using target protectors, we demonstrate that miR-430-mediated regulation of endogenous sdf1a (also known as cxcl12a) and cxcr7b (i) facilitates dynamic expression of sdf1a by clearing its mRNA from previous expression domains, (ii) modulates the levels of the decoy receptor Cxcr7b to avoid excessive depletion of Sdf1a and (iii) buffers against variation in gene dosage of chemokine signaling components to ensure accurate PGC migration. Our results indicate that losing miRNA-mediated regulation can expose otherwise buffered genetic lesions leading to developmental defects.