Extremely low-coverage sequencing and imputation increases power for genome-wide association studies

Genome-wide association studies (GWAS) have proven to be a powerful method to identify common genetic variants contributing to susceptibility to common diseases. Here, we show that extremely low-coverage sequencing (0.1-0.5×) captures almost as much of the common (>5%) and low-frequency (1-5%) variation across the genome as SNP arrays. As an empirical demonstration, we show that genome-wide SNP genotypes can be inferred at a mean r(2) of 0.71 using off-target data (0.24× average coverage) in a whole-exome study of 909 samples. Using both simulated and real exome-sequencing data sets, we show that association statistics obtained using extremely low-coverage sequencing data attain similar P values at known associated variants as data from genotyping arrays, without an excess of false positives. Within the context of reductions in sample preparation and sequencing costs, funds invested in extremely low-coverage sequencing can yield several times the effective sample size of GWAS based on SNP array data and a commensurate increase in statistical power.     

Common variants at 12q14 and 12q24 are associated with hippocampal volume

Aging is associated with reductions in hippocampal volume that are accelerated by Alzheimer's disease and vascular risk factors. Our genome-wide association study (GWAS) of dementia-free persons (n = 9,232) identified 46 SNPs at four loci with P values of <4.0 × 10(-7). In two additional samples (n = 2,318), associations were replicated at 12q14 within MSRB3-WIF1 (discovery and replication; rs17178006; P = 5.3 × 10(-11)) and at 12q24 near HRK-FBXW8 (rs7294919; P = 2.9 × 10(-11)). Remaining associations included one SNP at 2q24 within DPP4 (rs6741949; P = 2.9 × 10(-7)) and nine SNPs at 9p33 within ASTN2 (rs7852872; P = 1.0 × 10(-7)); along with the chromosome 12 associations, these loci were also associated with hippocampal volume (P < 0.05) in a third younger, more heterogeneous sample (n = 7,794). The SNP in ASTN2 also showed suggestive association with decline in cognition in a largely independent sample (n = 1,563). These associations implicate genes related to apoptosis (HRK), development (WIF1), oxidative stress (MSR3B), ubiquitination (FBXW8) and neuronal migration (ASTN2), as well as enzymes targeted by new diabetes medications (DPP4), indicating new genetic influences on hippocampal size and possibly the risk of cognitive decline and dementia.     

Loss-of-function mutations in TGFB2 cause a syndromic presentation of thoracic aortic aneurysm

Loeys-Dietz syndrome (LDS) associates with a tissue signature for high transforming growth factor (TGF)-β signaling but is often caused by heterozygous mutations in genes encoding positive effectors of TGF-β signaling, including either subunit of the TGF-β receptor or SMAD3, thereby engendering controversy regarding the mechanism of disease. Here, we report heterozygous mutations or deletions in the gene encoding the TGF-β2 ligand for a phenotype within the LDS spectrum and show upregulation of TGF-β signaling in aortic tissue from affected individuals. Furthermore, haploinsufficient Tgfb2(+/-) mice have aortic root aneurysm and biochemical evidence of increased canonical and noncanonical TGF-β signaling. Mice that harbor both a mutant Marfan syndrome (MFS) allele (Fbn1(C1039G/+)) and Tgfb2 haploinsufficiency show increased TGF-β signaling and phenotypic worsening in association with normalization of TGF-β2 expression and high expression of TGF-β1. Taken together, these data support the hypothesis that compensatory autocrine and/or paracrine events contribute to the pathogenesis of TGF-β-mediated vasculopathies.     

PNPLA1 mutations cause autosomal recessive congenital ichthyosis in golden retriever dogs and humans

Ichthyoses comprise a heterogeneous group of genodermatoses characterized by abnormal desquamation over the whole body, for which the genetic causes of several human forms remain unknown. We used a spontaneous dog model in the golden retriever breed, which is affected by a lamellar ichthyosis resembling human autosomal recessive congenital ichthyoses (ARCI), to carry out a genome-wide association study. We identified a homozygous insertion-deletion (indel) mutation in PNPLA1 that leads to a premature stop codon in all affected golden retriever dogs. We subsequently found one missense and one nonsense mutation in the catalytic domain of human PNPLA1 in six individuals with ARCI from two families. Further experiments highlighted the importance of PNPLA1 in the formation of the epidermal lipid barrier. This study identifies a new gene involved in human ichthyoses and provides insights into the localization and function of this yet uncharacterized member of the PNPLA protein family.     

Conditional and joint multiple-SNP analysis of GWAS summary statistics identifies additional variants influencing complex traits

We present an approximate conditional and joint association analysis that can use summary-level statistics from a meta-analysis of genome-wide association studies (GWAS) and estimated linkage disequilibrium (LD) from a reference sample with individual-level genotype data. Using this method, we analyzed meta-analysis summary data from the GIANT Consortium for height and body mass index (BMI), with the LD structure estimated from genotype data in two independent cohorts. We identified 36 loci with multiple associated variants for height (38 leading and 49 additional SNPs, 87 in total) via a genome-wide SNP selection procedure. The 49 new SNPs explain approximately 1.3% of variance, nearly doubling the heritability explained at the 36 loci. We did not find any locus showing multiple associated SNPs for BMI. The method we present is computationally fast and is also applicable to case-control data, which we demonstrate in an example from meta-analysis of type 2 diabetes by the DIAGRAM Consortium.     

Biogeographic and taxonomic relationships among the mountain snails (Gastropoda: Oreohelicidae) of the central Great Basin

Described here are 4 species of mountain snails, Oreohelix , isolated on mountains in the central Great Basin of Nevada and Utah since the end of the Pleistocene. Forty-three mountains were searched during an 18-year period, resulting in 24 mountains found with no oreohelicids present. One population, Oreohelix loisae (19 mm to 23 mm in shell diameter), is described here as a new species related to, but geographically isolated from, the species Oreohelix nevadensis (17 mm to 22 mm diameter). Oreohelix loisae is present only in the Goshute Mountains while O. nevadensis is represented in 3 geographically adjacent ranges in the central Great Basin. These 2 species are possibly related to the Oreohelix haydeni group from the northern Wasatch Range. The subspecies Oreohelix strigosa depressa (15 mm to 21 mm diameter) is present on 11 ranges from western Utah west to east central Nevada. This subspecies is closely related to populations found today in the northern Wasatch Mountains of Utah. The smallest species in diameter (8 mm to 14 mm), Oreohelix hemphilli , is centered in the central Great Basin and found on 16 ranges often in sympatry with 1 or 2 of the larger conspecifics. Both qualitative and quantitative information on shell characters and soft anatomy is provided here for these 4 species. Shell characters, soft anatomy, geographical isolation, and statistical analysis suggest that 4 distinct species inhabit the central Great Basin today. Xeric and calciphilic species include O. hemphilli and O. loisae , while O. strigosa and O. nevadensis typically are associated with permanent water and both metamorphic and limestone mountains.     

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.