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.     

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.     

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.     

Exome sequencing of liver fluke-associated cholangiocarcinoma

Opisthorchis viverrini-related cholangiocarcinoma (CCA), a fatal bile duct cancer, is a major public health concern in areas endemic for this parasite. We report here whole-exome sequencing of eight O. viverrini-related tumors and matched normal tissue. We identified and validated 206 somatic mutations in 187 genes using Sanger sequencing and selected 15 genes for mutation prevalence screening in an additional 46 individuals with CCA (cases). In addition to the known cancer-related genes TP53 (mutated in 44.4% of cases), KRAS (16.7%) and SMAD4 (16.7%), we identified somatic mutations in 10 newly implicated genes in 14.8-3.7% of cases. These included inactivating mutations in MLL3 (in 14.8% of cases), ROBO2 (9.3%), RNF43 (9.3%) and PEG3 (5.6%), and activating mutations in the GNAS oncogene (9.3%). These genes have functions that can be broadly grouped into three biological classes: (i) deactivation of histone modifiers, (ii) activation of G protein signaling and (iii) loss of genome stability. This study provides insight into the mutational landscape contributing to O. viverrini-related CCA.     

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.