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.     

Mutations in SWI/SNF chromatin remodeling complex gene ARID1B cause Coffin-Siris syndrome

We identified de novo truncating mutations in ARID1B in three individuals with Coffin-Siris syndrome (CSS) by exome sequencing. Array-based copy-number variation (CNV) analysis in 2,000 individuals with intellectual disability revealed deletions encompassing ARID1B in 3 subjects with phenotypes partially overlapping that of CSS. Taken together with published data, these results indicate that haploinsufficiency of the ARID1B gene, which encodes an epigenetic modifier of chromatin structure, is an important cause of CSS and is potentially a common cause of intellectual disability and speech impairment.     

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.     

Estimating the proportion of variation in susceptibility to schizophrenia captured by common SNPs

Schizophrenia is a complex disorder caused by both genetic and environmental factors. Using 9,087 affected individuals, 12,171 controls and 915,354 imputed SNPs from the Schizophrenia Psychiatric Genome-Wide Association Study (GWAS) Consortium (PGC-SCZ), we estimate that 23% (s.e. = 1%) of variation in liability to schizophrenia is captured by SNPs. We show that a substantial proportion of this variation must be the result of common causal variants, that the variance explained by each chromosome is linearly related to its length (r = 0.89, P = 2.6 × 10(-8)), that the genetic basis of schizophrenia is the same in males and females, and that a disproportionate proportion of variation is attributable to a set of 2,725 genes expressed in the central nervous system (CNS; P = 7.6 × 10(-8)). These results are consistent with a polygenic genetic architecture and imply more individual SNP associations will be detected for this disease as sample size increases.     

Common variants near TARDBP and EGR2 are associated with susceptibility to Ewing sarcoma

Ewing sarcoma, a pediatric tumor characterized by EWSR1-ETS fusions, is predominantly observed in populations of European ancestry. We performed a genome-wide association study (GWAS) of 401 French individuals with Ewing sarcoma, 684 unaffected French individuals and 3,668 unaffected individuals of European descent and living in the United States. We identified candidate risk loci at 1p36.22, 10q21 and 15q15. We replicated these loci in two independent sets of cases and controls. Joint analysis identified associations with rs9430161 (P = 1.4 × 10(-20); odds ratio (OR) = 2.2) located 25 kb upstream of TARDBP, rs224278 (P = 4.0 × 10(-17); OR = 1.7) located 5 kb upstream of EGR2 and, to a lesser extent, rs4924410 at 15q15 (P = 6.6 × 10(-9); OR = 1.5). The major risk haplotypes were less prevalent in Africans, suggesting that these loci could contribute to geographical differences in Ewing sarcoma incidence. TARDBP shares structural similarities with EWSR1 and FUS, which encode RNA binding proteins, and EGR2 is a target gene of EWSR1-ETS. Variants at these loci were associated with expression levels of TARDBP, ADO (encoding cysteamine dioxygenase) and EGR2.     

A genome-wide association study in Han Chinese identifies new susceptibility loci for ankylosing spondylitis

To identify susceptibility loci for ankylosing spondylitis, we performed a two-stage genome-wide association study in Han Chinese. In the discovery stage, we analyzed 1,356,350 autosomal SNPs in 1,837 individuals with ankylosing spondylitis and 4,231 controls; in the validation stage, we analyzed 30 suggestive SNPs in an additional 2,100 affected individuals and 3,496 controls. We identified two new susceptibility loci between EDIL3 and HAPLN1 at 5q14.3 (rs4552569; P = 8.77 × 10(-10)) and within ANO6 at 12q12 (rs17095830; P = 1.63 × 10(-8)). We also confirmed previously reported associations in Europeans within the major histocompatibility complex (MHC) region (top SNP, rs13202464; P < 5 × 10(-324)) and at 2p15 (rs10865331; P = 1.98 × 10(-8)). We show that rs13202464 within the MHC region mainly represents the risk effect of HLA-B*27 variants (including HLA-B*2704, HLA-B*2705 and HLA-B*2715) in Chinese. The two newly discovered loci implicate genes related to bone formation and cartilage development, suggesting their potential involvement in the etiology of ankylosing spondylitis.     

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.