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.     

Dominant missense mutations in ABCC9 cause Cantú syndrome

Cantú syndrome is characterized by congenital hypertrichosis, distinctive facial features, osteochondrodysplasia and cardiac defects. By using family-based exome sequencing, we identified a de novo mutation in ABCC9. Subsequently, we discovered novel dominant missense mutations in ABCC9 in 14 of the 16 individuals with Cantú syndrome examined. The ABCC9 protein is part of an ATP-dependent potassium (K(ATP)) channel that couples the metabolic state of a cell with its electrical activity. All mutations altered amino acids in or close to the transmembrane domains of ABCC9. Using electrophysiological measurements, we show that mutations in ABCC9 reduce the ATP-mediated potassium channel inhibition, resulting in channel opening. Moreover, similarities between the phenotype of individuals with Cantú syndrome and side effects from the K(ATP) channel agonist minoxidil indicate that the mutations in ABCC9 result in channel opening. Given the availability of ABCC9 antagonists, our findings may have direct implications for the treatment of individuals with Cantú syndrome.     

CCDC103 mutations cause primary ciliary dyskinesia by disrupting assembly of ciliary dynein arms

Cilia are essential for fertilization, respiratory clearance, cerebrospinal fluid circulation and establishing laterality. Cilia motility defects cause primary ciliary dyskinesia (PCD, MIM244400), a disorder affecting 1:15,000-30,000 births. Cilia motility requires the assembly of multisubunit dynein arms that drive ciliary bending. Despite progress in understanding the genetic basis of PCD, mutations remain to be identified for several PCD-linked loci. Here we show that the zebrafish cilia paralysis mutant schmalhans (smh(tn222)) encodes the coiled-coil domain containing 103 protein (Ccdc103), a foxj1a-regulated gene product. Screening 146 unrelated PCD families identified individuals in six families with reduced outer dynein arms who carried mutations in CCDC103. Dynein arm assembly in smh mutant zebrafish was rescued by wild-type but not mutant human CCDC103. Chlamydomonas Ccdc103/Pr46b functions as a tightly bound, axoneme-associated protein. These results identify Ccdc103 as a dynein arm attachment factor that causes primary ciliary dyskinesia when mutated.     

Mutations in the RNA exosome component gene EXOSC3 cause pontocerebellar hypoplasia and spinal motor neuron degeneration

RNA exosomes are multi-subunit complexes conserved throughout evolution and are emerging as the major cellular machinery for processing, surveillance and turnover of a diverse spectrum of coding and noncoding RNA substrates essential for viability. By exome sequencing, we discovered recessive mutations in EXOSC3 (encoding exosome component 3) in four siblings with infantile spinal motor neuron disease, cerebellar atrophy, progressive microcephaly and profound global developmental delay, consistent with pontocerebellar hypoplasia type 1 (PCH1; MIM 607596). We identified mutations in EXOSC3 in an additional 8 of 12 families with PCH1. Morpholino knockdown of exosc3 in zebrafish embryos caused embryonic maldevelopment, resulting in small brain size and poor motility, reminiscent of human clinical features, and these defects were largely rescued by co-injection with wild-type but not mutant exosc3 mRNA. These findings represent the first example of an RNA exosome core component gene that is responsible for a human disease and further implicate dysregulation of RNA processing in cerebellar and spinal motor neuron maldevelopment and degeneration.     

Mutations in ISPD cause Walker-Warburg syndrome and defective glycosylation of α-dystroglycan

Walker-Warburg syndrome (WWS) is an autosomal recessive multisystem disorder characterized by complex eye and brain abnormalities with congenital muscular dystrophy (CMD) and aberrant a-dystroglycan glycosylation. Here we report mutations in the ISPD gene (encoding isoprenoid synthase domain containing) as the second most common cause of WWS. Bacterial IspD is a nucleotidyl transferase belonging to a large glycosyltransferase family, but the role of the orthologous protein in chordates is obscure to date, as this phylum does not have the corresponding non-mevalonate isoprenoid biosynthesis pathway. Knockdown of ispd in zebrafish recapitulates the human WWS phenotype with hydrocephalus, reduced eye size, muscle degeneration and hypoglycosylated a-dystroglycan. These results implicate ISPD in a-dystroglycan glycosylation in maintaining sarcolemma integrity in vertebrates.     

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.     

Mouse ooplasm confers context-specific reprogramming capacity

Enucleated oocytes have the distinctive ability to reprogram somatic nuclei back to totipotency. Here, we investigate genome-scale DNA methylation patterns after nuclear transfer and compare them to the dynamics at fertilization. We identify specific targets for DNA demethylation after nuclear transfer, such as germline-associated promoters, as well as unique limitations that include certain repetitive element classes.     

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.     

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.