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.     

A20 (TNFAIP3) deficiency in myeloid cells triggers erosive polyarthritis resembling rheumatoid arthritis

A20 (TNFAIP3) is a protein that is involved in the negative feedback regulation of NF-κB signaling in response to specific proinflammatory stimuli in different cell types and has been suggested as a susceptibility gene for rheumatoid arthritis. To define the contribution of A20 to rheumatoid arthritis pathology, we generated myeloid-specific A20-deficient mice and show that specific ablation of Tnfaip3 in myeloid cells results in spontaneous development of a severe destructive polyarthritis with many features of rheumatoid arthritis. Myeloid-A20-deficient mice have high levels of inflammatory cytokines in their serum, consistent with a sustained NF-κB activation and higher TNF production by macrophages. Destructive polyarthritis in myeloid A20 knockout mice was TLR4-MyD88 and IL-6 dependent but was TNF independent. Myeloid A20 deficiency also promoted osteoclastogenesis in mice. Together, these observations indicate a critical and cell-specific function for A20 in the etiology of rheumatoid arthritis, supporting the idea of developing A20 modulatory drugs as cell-targeted therapies.     

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.     

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.