Resequencing of 200 human exomes identifies an excess of low-frequency non-synonymous coding variants

Targeted capture combined with massively parallel exome sequencing is a promising approach to identify genetic variants implicated in human traits. We report exome sequencing of 200 individuals from Denmark with targeted capture of 18,654 coding genes and sequence coverage of each individual exome at an average depth of 12-fold. On average, about 95% of the target regions were covered by at least one read. We identified 121,870 SNPs in the sample population, including 53,081 coding SNPs (cSNPs). Using a statistical method for SNP calling and an estimation of allelic frequencies based on our population data, we derived the allele frequency spectrum of cSNPs with a minor allele frequency greater than 0.02. We identified a 1.8-fold excess of deleterious, non-syonomyous cSNPs over synonymous cSNPs in the low-frequency range (minor allele frequencies between 2% and 5%). This excess was more pronounced for X-linked SNPs, suggesting that deleterious substitutions are primarily recessive.     

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.     

Mosaic overgrowth with fibroadipose hyperplasia is caused by somatic activating mutations in PIK3CA

The phosphatidylinositol 3-kinase (PI3K)-AKT signaling pathway is critical for cellular growth and metabolism. Correspondingly, loss of function of PTEN, a negative regulator of PI3K, or activating mutations in AKT1, AKT2 or AKT3 have been found in distinct disorders featuring overgrowth or hypoglycemia. We performed exome sequencing of DNA from unaffected and affected cells from an individual with an unclassified syndrome of congenital progressive segmental overgrowth of fibrous and adipose tissue and bone and identified the cancer-associated mutation encoding p.His1047Leu in PIK3CA, the gene that encodes the p110α catalytic subunit of PI3K, only in affected cells. Sequencing of PIK3CA in ten additional individuals with overlapping syndromes identified either the p.His1047Leu alteration or a second cancer-associated alteration, p.His1047Arg, in nine cases. Affected dermal fibroblasts showed enhanced basal and epidermal growth factor (EGF)-stimulated phosphatidylinositol 3,4,5-trisphosphate (PIP(3)) generation and concomitant activation of downstream signaling relative to their unaffected counterparts. Our findings characterize a distinct overgrowth syndrome, biochemically demonstrate activation of PI3K signaling and thereby identify a rational therapeutic target.     

De novo mutations in ATP1A3 cause alternating hemiplegia of childhood

Alternating hemiplegia of childhood (AHC) is a rare, severe neurodevelopmental syndrome characterized by recurrent hemiplegic episodes and distinct neurological manifestations. AHC is usually a sporadic disorder and has unknown etiology. We used exome sequencing of seven patients with AHC and their unaffected parents to identify de novo nonsynonymous mutations in ATP1A3 in all seven individuals. In a subsequent sequence analysis of ATP1A3 in 98 other patients with AHC, we found that ATP1A3 mutations were likely to be responsible for at least 74% of the cases; we also identified one inherited mutation in a case of familial AHC. Notably, most AHC cases are caused by one of seven recurrent ATP1A3 mutations, one of which was observed in 36 patients. Unlike ATP1A3 mutations that cause rapid-onset dystonia-parkinsonism, AHC-causing mutations in this gene caused consistent reductions in ATPase activity without affecting the level of protein expression. This work identifies de novo ATP1A3 mutations as the primary cause of AHC and offers insight into disease pathophysiology by expanding the spectrum of phenotypes associated with mutations in ATP1A3.     

De novo mutations in the actin genes ACTB and ACTG1 cause Baraitser-Winter syndrome

Brain malformations are individually rare but collectively common causes of developmental disabilities. Many forms of malformation occur sporadically and are associated with reduced reproductive fitness, pointing to a causative role for de novo mutations. Here, we report a study of Baraitser-Winter syndrome, a well-defined disorder characterized by distinct craniofacial features, ocular colobomata and neuronal migration defect. Using whole-exome sequencing of three proband-parent trios, we identified de novo missense changes in the cytoplasmic actin-encoding genes ACTB and ACTG1 in one and two probands, respectively. Sequencing of both genes in 15 additional affected individuals identified disease-causing mutations in all probands, including two recurrent de novo alterations (ACTB, encoding p.Arg196His, and ACTG1, encoding p.Ser155Phe). Our results confirm that trio-based exome sequencing is a powerful approach to discover genes causing sporadic developmental disorders, emphasize the overlapping roles of cytoplasmic actin proteins in development and suggest that Baraitser-Winter syndrome is the predominant phenotype associated with mutation of these two genes.     

De novo germline and postzygotic mutations in AKT3, PIK3R2 and PIK3CA cause a spectrum of related megalencephaly syndromes

Megalencephaly-capillary malformation (MCAP) and megalencephaly-polymicrogyria-polydactyly-hydrocephalus (MPPH) syndromes are sporadic overgrowth disorders associated with markedly enlarged brain size and other recognizable features. We performed exome sequencing in 3 families with MCAP or MPPH, and our initial observations were confirmed in exomes from 7 individuals with MCAP and 174 control individuals, as well as in 40 additional subjects with megalencephaly, using a combination of Sanger sequencing, restriction enzyme assays and targeted deep sequencing. We identified de novo germline or postzygotic mutations in three core components of the phosphatidylinositol 3-kinase (PI3K)-AKT pathway. These include 2 mutations in AKT3, 1 recurrent mutation in PIK3R2 in 11 unrelated families with MPPH and 15 mostly postzygotic mutations in PIK3CA in 23 individuals with MCAP and 1 with MPPH. Our data highlight the central role of PI3K-AKT signaling in vascular, limb and brain development and emphasize the power of massively parallel sequencing in a challenging context of phenotypic and genetic heterogeneity combined with postzygotic mosaicism.     

Exome sequencing identifies somatic mutations of DNA methyltransferase gene DNMT3A in acute monocytic leukemia

Abnormal epigenetic regulation has been implicated in oncogenesis. We report here the identification of somatic mutations by exome sequencing in acute monocytic leukemia, the M5 subtype of acute myeloid leukemia (AML-M5). We discovered mutations in DNMT3A (encoding DNA methyltransferase 3A) in 23 of 112 (20.5%) cases. The DNMT3A mutants showed reduced enzymatic activity or aberrant affinity to histone H3 in vitro. Notably, there were alterations of DNA methylation patterns and/or gene expression profiles (such as HOXB genes) in samples with DNMT3A mutations as compared with those without such changes. Leukemias with DNMT3A mutations constituted a group of poor prognosis with elderly disease onset and of promonocytic as well as monocytic predominance among AML-M5 individuals. Screening other leukemia subtypes showed Arg882 alterations in 13.6% of acute myelomonocytic leukemia (AML-M4) cases. Our work suggests a contribution of aberrant DNA methyltransferase activity to the pathogenesis of acute monocytic leukemia and provides a useful new biomarker for relevant cases.     

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 affecting components of the SWI/SNF complex cause Coffin-Siris syndrome

By exome sequencing, we found de novo SMARCB1 mutations in two of five individuals with typical Coffin-Siris syndrome (CSS), a rare autosomal dominant anomaly syndrome. As SMARCB1 encodes a subunit of the SWItch/Sucrose NonFermenting (SWI/SNF) complex, we screened 15 other genes encoding subunits of this complex in 23 individuals with CSS. Twenty affected individuals (87%) each had a germline mutation in one of six SWI/SNF subunit genes, including SMARCB1, SMARCA4, SMARCA2, SMARCE1, ARID1A and ARID1B.     

Mutations in NNT encoding nicotinamide nucleotide transhydrogenase cause familial glucocorticoid deficiency

Using targeted exome sequencing, we identified mutations in NNT, an antioxidant defense gene, in individuals with familial glucocorticoid deficiency. In mice with Nnt loss, higher levels of adrenocortical cell apoptosis and impaired glucocorticoid production were observed. NNT knockdown in a human adrenocortical cell line resulted in impaired redox potential and increased reactive oxygen species (ROS) levels. Our results suggest that NNT may have a role in ROS detoxification in human adrenal glands.     

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.     

ISPD loss-of-function mutations disrupt dystroglycan O-mannosylation and cause Walker-Warburg syndrome

Walker-Warburg syndrome (WWS) is clinically defined as congenital muscular dystrophy that is accompanied by a variety of brain and eye malformations. It represents the most severe clinical phenotype in a spectrum of diseases associated with abnormal post-translational processing of a-dystroglycan that share a defect in laminin-binding glycan synthesis1. Although mutations in six genes have been identified as causes of WWS, only half of all individuals with the disease can currently be diagnosed on this basis2. A cell fusion complementation assay in fibroblasts from undiagnosed individuals with WWS was used to identify five new complementation groups. Further evaluation of one group by linkage analysis and targeted sequencing identified recessive mutations in the ISPD gene (encoding isoprenoid synthase domain containing). The pathogenicity of the identified ISPD mutations was shown by complementation of fibroblasts with wild-type ISPD. Finally, we show that recessive mutations in ISPD abolish the initial step in laminin-binding glycan synthesis by disrupting dystroglycan O-mannosylation. This establishes a new mechanism for WWS pathophysiology.     

Mutations in MEGF10, a regulator of satellite cell myogenesis, cause early onset myopathy, areflexia, respiratory distress and dysphagia (EMARDD)

Infantile myopathies with diaphragmatic paralysis are genetically heterogeneous, and clinical symptoms do not assist in differentiating between them. We used phased haplotype analysis with subsequent targeted exome sequencing to identify MEGF10 mutations in a previously unidentified type of infantile myopathy with diaphragmatic weakness, areflexia, respiratory distress and dysphagia. MEGF10 is highly expressed in activated satellite cells and regulates their proliferation as well as their differentiation and fusion into multinucleated myofibers, which are greatly reduced in muscle from individuals with early onset myopathy, areflexia, respiratory distress and dysphagia.     

Mutations in genes encoding subunits of RNA polymerases I and III cause Treacher Collins syndrome

We identified a deletion of a gene encoding a subunit of RNA polymerases I and III, POLR1D, in an individual with Treacher Collins syndrome (TCS). Subsequently, we detected 20 additional heterozygous mutations of POLR1D in 252 individuals with TCS. Furthermore, we discovered mutations in both alleles of POLR1C in three individuals with TCS. These findings identify two additional genes involved in TCS, confirm the genetic heterogeneity of TCS and support the hypothesis that TCS is a ribosomopathy.     

A de novo paradigm for mental retardation

The per-generation mutation rate in humans is high. De novo mutations may compensate for allele loss due to severely reduced fecundity in common neurodevelopmental and psychiatric diseases, explaining a major paradox in evolutionary genetic theory. Here we used a family based exome sequencing approach to test this de novo mutation hypothesis in ten individuals with unexplained mental retardation. We identified and validated unique non-synonymous de novo mutations in nine genes. Six of these, identified in six different individuals, are likely to be pathogenic based on gene function, evolutionary conservation and mutation impact. Our findings provide strong experimental support for a de novo paradigm for mental retardation. Together with de novo copy number variation, de novo point mutations of large effect could explain the majority of all mental retardation cases in the population.     

FAN1 mutations cause karyomegalic interstitial nephritis, linking chronic kidney failure to defective DNA damage repair

Chronic kidney disease (CKD) represents a major health burden. Its central feature of renal fibrosis is not well understood. By exome sequencing, we identified mutations in FAN1 as a cause of karyomegalic interstitial nephritis (KIN), a disorder that serves as a model for renal fibrosis. Renal histology in KIN is indistinguishable from that of nephronophthisis, except for the presence of karyomegaly. The FAN1 protein has nuclease activity and acts in DNA interstrand cross-link (ICL) repair within the Fanconi anemia DNA damage response (DDR) pathway. We show that cells from individuals with FAN1 mutations have sensitivity to the ICL-inducing agent mitomycin C but do not exhibit chromosome breakage or cell cycle arrest after diepoxybutane treatment, unlike cells from individuals with Fanconi anemia. We complemented ICL sensitivity with wild-type FAN1 but not with cDNA having mutations found in individuals with KIN. Depletion of fan1 in zebrafish caused increased DDR, apoptosis and kidney cysts. Our findings implicate susceptibility to environmental genotoxins and inadequate DNA repair as novel mechanisms contributing to renal fibrosis and CKD.     

Exome sequencing in sporadic autism spectrum disorders identifies severe de novo mutations

Evidence for the etiology of autism spectrum disorders (ASDs) has consistently pointed to a strong genetic component complicated by substantial locus heterogeneity. We sequenced the exomes of 20 individuals with sporadic ASD (cases) and their parents, reasoning that these families would be enriched for de novo mutations of major effect. We identified 21 de novo mutations, 11 of which were protein altering. Protein-altering mutations were significantly enriched for changes at highly conserved residues. We identified potentially causative de novo events in 4 out of 20 probands, particularly among more severely affected individuals, in FOXP1, GRIN2B, SCN1A and LAMC3. In the FOXP1 mutation carrier, we also observed a rare inherited CNTNAP2 missense variant, and we provide functional support for a multi-hit model for disease risk. Our results show that trio-based exome sequencing is a powerful approach for identifying new candidate genes for ASDs and suggest that de novo mutations may contribute substantially to the genetic etiology of ASDs.     

Exome sequencing of extreme phenotypes identifies DCTN4 as a modifier of chronic Pseudomonas aeruginosa infection in cystic fibrosis

Exome sequencing has become a powerful and effective strategy for the discovery of genes underlying Mendelian disorders. However, use of exome sequencing to identify variants associated with complex traits has been more challenging, partly because the sample sizes needed for adequate power may be very large. One strategy to increase efficiency is to sequence individuals who are at both ends of a phenotype distribution (those with extreme phenotypes). Because the frequency of alleles that contribute to the trait are enriched in one or both phenotype extremes, a modest sample size can potentially be used to identify novel candidate genes and/or alleles. As part of the National Heart, Lung, and Blood Institute (NHLBI) Exome Sequencing Project (ESP), we used an extreme phenotype study design to discover that variants in DCTN4, encoding a dynactin protein, are associated with time to first P. aeruginosa airway infection, chronic P. aeruginosa infection and mucoid P. aeruginosa in individuals with cystic fibrosis.     

Exome sequencing identifies recurrent somatic MAP2K1 and MAP2K2 mutations in melanoma

We performed exome sequencing to detect somatic mutations in protein-coding regions in seven melanoma cell lines and donor-matched germline cells. All melanoma samples had high numbers of somatic mutations, which showed the hallmark of UV-induced DNA repair. Such a hallmark was absent in tumor sample-specific mutations in two metastases derived from the same individual. Two melanomas with non-canonical BRAF mutations harbored gain-of-function MAP2K1 and MAP2K2 (MEK1 and MEK2, respectively) mutations, resulting in constitutive ERK phosphorylation and higher resistance to MEK inhibitors. Screening a larger cohort of individuals with melanoma revealed the presence of recurring somatic MAP2K1 and MAP2K2 mutations, which occurred at an overall frequency of 8%. Furthermore, missense and nonsense somatic mutations were frequently found in three candidate melanoma genes, FAT4, LRP1B and DSC1.