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.     

A mixed-model approach for genome-wide association studies of correlated traits in structured populations

Genome-wide association studies (GWAS) are a standard approach for studying the genetics of natural variation. A major concern in GWAS is the need to account for the complicated dependence structure of the data, both between loci as well as between individuals. Mixed models have emerged as a general and flexible approach for correcting for population structure in GWAS. Here, we extend this linear mixed-model approach to carry out GWAS of correlated phenotypes, deriving a fully parameterized multi-trait mixed model (MTMM) that considers both the within-trait and between-trait variance components simultaneously for multiple traits. We apply this to data from a human cohort for correlated blood lipid traits from the Northern Finland Birth Cohort 1966 and show greatly increased power to detect pleiotropic loci that affect more than one blood lipid trait. We also apply this approach to an Arabidopsis thaliana data set for flowering measurements in two different locations, identifying loci whose effect depends on the environment.     

Activation of the AXL kinase causes resistance to EGFR-targeted therapy in lung cancer

Human non-small cell lung cancers (NSCLCs) with activating mutations in EGFR frequently respond to treatment with EGFR-targeted tyrosine kinase inhibitors (TKIs), such as erlotinib, but responses are not durable, as tumors acquire resistance. Secondary mutations in EGFR (such as T790M) or upregulation of the MET kinase are found in over 50% of resistant tumors. Here, we report increased activation of AXL and evidence for epithelial-to-mesenchymal transition (EMT) in multiple in vitro and in vivo EGFR-mutant lung cancer models with acquired resistance to erlotinib in the absence of the EGFR p.Thr790Met alteration or MET activation. Genetic or pharmacological inhibition of AXL restored sensitivity to erlotinib in these tumor models. Increased expression of AXL and, in some cases, of its ligand GAS6 was found in EGFR-mutant lung cancers obtained from individuals with acquired resistance to TKIs. These data identify AXL as a promising therapeutic target whose inhibition could prevent or overcome acquired resistance to EGFR TKIs in individuals with EGFR-mutant lung cancer.     

An efficient multi-locus mixed-model approach for genome-wide association studies in structured populations

Population structure causes genome-wide linkage disequilibrium between unlinked loci, leading to statistical confounding in genome-wide association studies. Mixed models have been shown to handle the confounding effects of a diffuse background of large numbers of loci of small effect well, but they do not always account for loci of larger effect. Here we propose a multi-locus mixed model as a general method for mapping complex traits in structured populations. Simulations suggest that our method outperforms existing methods in terms of power as well as false discovery rate. We apply our method to human and Arabidopsis thaliana data, identifying new associations and evidence for allelic heterogeneity. We also show how a priori knowledge from an A. thaliana linkage mapping study can be integrated into our method using a Bayesian approach. Our implementation is computationally efficient, making the analysis of large data sets (n > 10,000) practicable.     

Mutations in NMNAT1 cause Leber congenital amaurosis and identify a new disease pathway for retinal degeneration

Leber congenital amaurosis (LCA) is a blinding retinal disease that presents within the first year after birth. Using exome sequencing, we identified mutations in the nicotinamide adenine dinucleotide (NAD) synthase gene NMNAT1 encoding nicotinamide mononucleotide adenylyltransferase 1 in eight families with LCA, including the family in which LCA was originally linked to the LCA9 locus. Notably, all individuals with NMNAT1 mutations also have macular colobomas, which are severe degenerative entities of the central retina (fovea) devoid of tissue and photoreceptors. Functional assays of the proteins encoded by the mutant alleles identified in our study showed that the mutations reduce the enzymatic activity of NMNAT1 in NAD biosynthesis and affect protein folding. Of note, recent characterization of the slow Wallerian degeneration (Wld(s)) mouse model, in which prolonged axonal survival after injury is observed, identified NMNAT1 as a neuroprotective protein when ectopically expressed. Our findings identify a new disease mechanism underlying LCA and provide the first link between endogenous NMNAT1 dysfunction and a human nervous system disorder.     

Fast and accurate genotype imputation in genome-wide association studies through pre-phasing

The 1000 Genomes Project and disease-specific sequencing efforts are producing large collections of haplotypes that can be used as reference panels for genotype imputation in genome-wide association studies (GWAS). However, imputing from large reference panels with existing methods imposes a high computational burden. We introduce a strategy called 'pre-phasing' that maintains the accuracy of leading methods while reducing computational costs. We first statistically estimate the haplotypes for each individual within the GWAS sample (pre-phasing) and then impute missing genotypes into these estimated haplotypes. This reduces the computational cost because (i) the GWAS samples must be phased only once, whereas standard methods would implicitly repeat phasing with each reference panel update, and (ii) it is much faster to match a phased GWAS haplotype to one reference haplotype than to match two unphased GWAS genotypes to a pair of reference haplotypes. We implemented our approach in the MaCH and IMPUTE2 frameworks, and we tested it on data sets from the Wellcome Trust Case Control Consortium 2 (WTCCC2), the Genetic Association Information Network (GAIN), the Women's Health Initiative (WHI) and the 1000 Genomes Project. This strategy will be particularly valuable for repeated imputation as reference panels evolve.     

Maize HapMap2 identifies extant variation from a genome in flux

Whereas breeders have exploited diversity in maize for yield improvements, there has been limited progress in using beneficial alleles in undomesticated varieties. Characterizing standing variation in this complex genome has been challenging, with only a small fraction of it described to date. Using a population genetics scoring model, we identified 55 million SNPs in 103 lines across pre-domestication and domesticated Zea mays varieties, including a representative from the sister genus Tripsacum. We find that structural variations are pervasive in the Z. mays genome and are enriched at loci associated with important traits. By investigating the drivers of genome size variation, we find that the larger Tripsacum genome can be explained by transposable element abundance rather than an allopolyploid origin. In contrast, intraspecies genome size variation seems to be controlled by chromosomal knob content. There is tremendous overlap in key gene content in maize and Tripsacum, suggesting that adaptations from Tripsacum (for example, perennialism and frost and drought tolerance) can likely be integrated into maize.     

Genome-wide copy number variation study associates metabotropic glutamate receptor gene networks with attention deficit hyperactivity disorder

Attention deficit hyperactivity disorder (ADHD) is a common, heritable neuropsychiatric disorder of unknown etiology. We performed a whole-genome copy number variation (CNV) study on 1,013 cases with ADHD and 4,105 healthy children of European ancestry using 550,000 SNPs. We evaluated statistically significant findings in multiple independent cohorts, with a total of 2,493 cases with ADHD and 9,222 controls of European ancestry, using matched platforms. CNVs affecting metabotropic glutamate receptor genes were enriched across all cohorts (P = 2.1 × 10(-9)). We saw GRM5 (encoding glutamate receptor, metabotropic 5) deletions in ten cases and one control (P = 1.36 × 10(-6)). We saw GRM7 deletions in six cases, and we saw GRM8 deletions in eight cases and no controls. GRM1 was duplicated in eight cases. We experimentally validated the observed variants using quantitative RT-PCR. A gene network analysis showed that genes interacting with the genes in the GRM family are enriched for CNVs in ～10% of the cases (P = 4.38 × 10(-10)) after correction for occurrence in the controls. We identified rare recurrent CNVs affecting glutamatergic neurotransmission genes that were overrepresented in multiple ADHD cohorts.