Pneumococcal genome sequencing tracks a vaccine escape variant formed through a multi-fragment recombination event

Streptococcus pneumoniae ('pneumococcus') causes an estimated 14.5 million cases of serious disease and 826,000 deaths annually in children under 5 years of age(1). The highly effective introduction of the PCV7 pneumococcal vaccine in 2000 in the United States(2,3) provided an unprecedented opportunity to investigate the response of an important pathogen to widespread, vaccine-induced selective pressure. Here, we use array-based sequencing of 62 isolates from a US national monitoring program to study five independent instances of vaccine escape recombination(4), showing the simultaneous transfer of multiple and often large (up to at least 44 kb) DNA fragments. We show that one such new strain quickly became established, spreading from east to west across the United States. These observations clarify the roles of recombination and selection in the population genomics of pneumococcus and provide proof of principle of the considerable value of combining genomic and epidemiological information in the surveillance and enhanced understanding of infectious diseases.     

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.     

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.     

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.     

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.     

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.