Actionable exomic incidental findings in 6503 participants: challenges of variant classification.

Recommendations for laboratories to report incidental findings from genomic tests have stimulated interest in such results. In order to investigate the criteria and processes for assigning the pathogenicity of specific variants and to estimate the frequency of such incidental findings in patients of European and African ancestry, we classified potentially actionable pathogenic single-nucleotide variants (SNVs) in all 4300 European- and 2203 African-ancestry participants sequenced by the NHLBI Exome Sequencing Project (ESP). We considered 112 gene-disease pairs selected by an expert panel as associated with medically actionable genetic disorders that may be undiagnosed in adults. The resulting classifications were compared to classifications from other clinical and research genetic testing laboratories, as well as with in silico pathogenicity scores. Among European-ancestry participants, 30 of 4300 (0.7%) had a pathogenic SNV and six (0.1%) had a disruptive variant that was expected to be pathogenic, whereas 52 (1.2%) had likely pathogenic SNVs. For African-ancestry participants, six of 2203 (0.3%) had a pathogenic SNV and six (0.3%) had an expected pathogenic disruptive variant, whereas 13 (0.6%) had likely pathogenic SNVs. Genomic Evolutionary Rate Profiling mammalian conservation score and the Combined Annotation Dependent Depletion summary score of conservation, substitution, regulation, and other evidence were compared across pathogenicity assignments and appear to have utility in variant classification. This work provides a refined estimate of the burden of adult onset, medically actionable incidental findings expected from exome sequencing, highlights challenges in variant classification, and demonstrates the need for a better curated variant interpretation knowledge base.

Sporadic autism exomes reveal a highly interconnected protein network of de novo mutations.

It is well established that autism spectrum disorders (ASD) have a strong genetic component; however, for at least 70% of cases, the underlying genetic cause is unknown. Under the hypothesis that de novo mutations underlie a substantial fraction of the risk for developing ASD in families with no previous history of ASD or related phenotypes--so-called sporadic or simplex families--we sequenced all coding regions of the genome (the exome) for parent-child trios exhibiting sporadic ASD, including 189 new trios and 20 that were previously reported. Additionally, we also sequenced the exomes of 50 unaffected siblings corresponding to these new (n = 31) and previously reported trios (n = 19), for a total of 677 individual exomes from 209 families. Here we show that de novo point mutations are overwhelmingly paternal in origin (4:1 bias) and positively correlated with paternal age, consistent with the modest increased risk for children of older fathers to develop ASD. Moreover, 39% (49 of 126) of the most severe or disruptive de novo mutations map to a highly interconnected ß-catenin/chromatin remodelling protein network ranked significantly for autism candidate genes. In proband exomes, recurrent protein-altering mutations were observed in two genes: CHD8 and NTNG1. Mutation screening of six candidate genes in 1,703 ASD probands identified additional de novo, protein-altering mutations in GRIN2B, LAMC3 and SCN1A. Combined with copy number variant (CNV) data, these results indicate extreme locus heterogeneity but also provide a target for future discovery, diagnostics and therapeutics.

Whole-exome sequencing identifies rare and low-frequency coding variants associated with LDL cholesterol.

Elevated low-density lipoprotein cholesterol (LDL-C) is a treatable, heritable risk factor for cardiovascular disease. Genome-wide association studies (GWASs) have identified 157 variants associated with lipid levels but are not well suited to assess the impact of rare and low-frequency variants. To determine whether rare or low-frequency coding variants are associated with LDL-C, we exome sequenced 2,005 individuals, including 554 individuals selected for extreme LDL-C (>98(th) or <2(nd) percentile). Follow-up analyses included sequencing of 1,302 additional individuals and genotype-based analysis of 52,221 individuals. We observed significant evidence of association between LDL-C and the burden of rare or low-frequency variants in PNPLA5, encoding a phospholipase-domain-containing protein, and both known and previously unidentified variants in PCSK9, LDLR and APOB, three known lipid-related genes. The effect sizes for the burden of rare variants for each associated gene were substantially higher than those observed for individual SNPs identified from GWASs. We replicated the PNPLA5 signal in an independent large-scale sequencing study of 2,084 individuals. In conclusion, this large whole-exome-sequencing study for LDL-C identified a gene not known to be implicated in LDL-C and provides unique insight into the design and analysis of similar experiments.

Frequent PIK3CA Mutations in Colorectal and Endometrial Tumors With 2 or More Somatic Mutations in Mismatch Repair Genes.

BACKGROUND & AIMS: Some colorectal and endometrial tumors with microsatellite instability not attributable to MLH1 hypermethylation or germline mutations contain 2 or more somatic mutations in genes encoding mismatch repair (MMR) proteins. We sought to define the molecular phenotype of this newly recognized tumor subtype. METHODS: From 2 prospective studies of the efficacy of screening for Lynch syndrome, we identified patients with colorectal and endometrial tumors who had 2 or more somatic (but not germline) mutations in genes encoding MMR proteins (double somatic). We determined the frequencies of tumor mutations in PIK3CA, BRAF, KRAS, NRAS, and PTEN by targeted next-generation sequencing and used logistic-regression models to compare them with those from patients with Lynch syndrome, MLH1-hypermethylated, or microsatellite-stable tumors. We validated our findings using independent data sets from The Cancer Genome Atlas. RESULTS: Among colorectal cancer cases, we found that 14 of 21 (67%) patients with double somatic tumors also had PIK3CA mutations, compared with 4 of 18 (22%) tumors from patients with Lynch syndrome, 2 of 10 (20%) tumors with MLH1 hypermethylation, and 12 of 78 (15%) tumors with microsatellite stability (P < .0001 for patients with double somatic tumors vs other subgroups). Mutations in PIK3CA were detected in all 13 patients with double somatic endometrial cancers (P = .04 compared with other subgroups). We did not detect BRAF mutations in patients with double somatic colorectal tumors or Lynch syndrome. We found highly similar results in a validation cohort from The Cancer Genome Atlas (113 patients with colorectal tumors, 178 endometrial tumors); 100% of double somatic cases had a somatic mutation in PIK3CA (P < .0001 compared with other subgroups). CONCLUSIONS: Most patients with colorectal or endometrial tumors with 2 or more somatic (but not germline) mutations in MMR proteins also have mutations in PIK3CA; mutations in PIK3CA are detected at substantially higher frequencies in these double somatic tumors than in other microsatellite-instability subgroups. PIK3CA mutation status might be used to identify a specific group of colorectal tumors, and to select treatment or determine prognosis.

A new congenital disorder of glycosylation caused by a mutation in SSR4, the signal sequence receptor 4 protein of the TRAP complex.

Nearly 50 congenital disorders of glycosylation (CDG) are known, but many patients biochemically diagnosed with CDG do not have mutations in known genes. Here, we describe a 16-year-old male who was born with microcephaly, developed intellectual disability, gastroesophageal reflux and a seizure disorder. We identified a de novo variant in the X-linked SSR4 gene which encodes a protein of the heterotetrameric translocon-associated protein (TRAP) complex. The c.316delT causes a p.F106Sfs*53 in SSR4 and also reduces expression of other TRAP complex proteins. The glycosylation marker Glyc-ER-GFP was used to confirm the underglycosylation in fibroblasts from the patient. Overexpression of the wild-type SSR4 allele partially restores glycosylation of the marker and of the other members of the TRAP complex. This is the first evidence that the TRAP complex, which binds to the oligosaccharyltransferase complex, is directly involved in N-glycosylation.

WNT1 mutations in families affected by moderately severe and progressive recessive osteogenesis imperfecta.

Osteogenesis imperfecta (OI) is a heritable disorder that ranges in severity from death in the perinatal period to an increased lifetime risk of fracture. Mutations in COL1A1 and COL1A2, which encode the chains of type I procollagen, result in dominant forms of OI, and mutations in several other genes result in recessive forms of OI. Here, we describe four recessive-OI-affected families in which we identified causative mutations in wingless-type MMTV integration site family 1 (WNT1). In family 1, we identified a homozygous missense mutation by exome sequencing. In family 2, we identified a homozygous nonsense mutation predicted to produce truncated WNT1. In family 3, we found a nonsense mutation and a single-nucleotide duplication on different alleles, and in family 4, we found a homozygous 14 bp deletion. The mutations in families 3 and 4 are predicted to result in nonsense-mediated mRNA decay and the absence of WNT1. WNT1 is a secreted signaling protein that binds the frizzled receptor (FZD) and the coreceptor low-density lipoprotein-receptor-related protein 5 (LRP5). Biallelic loss-of-function mutations in LRP5 result in recessive osteoporosis-pseudoglioma syndrome with low bone mass, whereas heterozygous gain-of-function mutations result in van Buchem disease with elevated bone density. Biallelic loss-of-function mutations in WNT1 result in a recessive clinical picture that includes bone fragility with a moderately severe and progressive presentation that is not easily distinguished from dominant OI type III.

Mosaicism of the UDP-galactose transporter SLC35A2 causes a congenital disorder of glycosylation.

Biochemical analysis and whole-exome sequencing identified mutations in the Golgi-localized UDP-galactose transporter SLC35A2 that define an undiagnosed X-linked congenital disorder of glycosylation (CDG) in three unrelated families. Each mutation reduced UDP-galactose transport, leading to galactose-deficient glycoproteins. Two affected males were somatic mosaics, suggesting that a wild-type SLC35A2 allele may be required for survival. In infancy, the commonly used biomarker transferrin showed abnormal glycosylation, but its appearance became normal later in childhood, without any corresponding clinical improvement. This may indicate selection against cells carrying the mutant allele. To detect other individuals with such mutations, we suggest transferrin testing in infancy. Here, we report somatic mosaicism in CDG, and our work stresses the importance of combining both genetic and biochemical diagnoses.

Whole-genome analysis reveals that mutations in inositol polyphosphate phosphatase-like 1 cause opsismodysplasia.

Opsismodysplasia is a rare, autosomal-recessive skeletal dysplasia characterized by short stature, characteristic facial features, and in some cases severe renal phosphate wasting. We used linkage analysis and whole-genome sequencing of a consanguineous trio to discover that mutations in inositol polyphosphate phosphatase-like 1 (INPPL1) cause opsismodysplasia with or without renal phosphate wasting. Evaluation of 12 families with opsismodysplasia revealed that INPPL1 mutations explain ~60% of cases overall, including both of the families in our cohort with more than one affected child and 50% of the simplex cases.

Actionable, pathogenic incidental findings in 1,000 participants' exomes.

The incorporation of genomics into medicine is stimulating interest on the return of incidental findings (IFs) from exome and genome sequencing. However, no large-scale study has yet estimated the number of expected actionable findings per individual; therefore, we classified actionable pathogenic single-nucleotide variants in 500 European- and 500 African-descent participants randomly selected from the National Heart, Lung, and Blood Institute Exome Sequencing Project. The 1,000 individuals were screened for variants in 114 genes selected by an expert panel for their association with medically actionable genetic conditions possibly undiagnosed in adults. Among the 1,000 participants, 585 instances of 239 unique variants were identified as disease causing in the Human Gene Mutation Database (HGMD). The primary literature supporting the variants' pathogenicity was reviewed. Of the identified IFs, only 16 unique autosomal-dominant variants in 17 individuals were assessed to be pathogenic or likely pathogenic, and one participant had two pathogenic variants for an autosomal-recessive disease. Furthermore, one pathogenic and four likely pathogenic variants not listed as disease causing in HGMD were identified. These data can provide an estimate of the frequency (∼3.4% for European descent and ∼1.2% for African descent) of the high-penetrance actionable pathogenic or likely pathogenic variants in adults. The 23 participants with pathogenic or likely pathogenic variants were disproportionately of European (17) versus African (6) descent. The process of classifying these variants underscores the need for a more comprehensive and diverse centralized resource to provide curated information on pathogenicity for clinical use to minimize health disparities in genomic medicine.

Exome sequencing identifies mutations in CCDC114 as a cause of primary ciliary dyskinesia.

Primary ciliary dyskinesia (PCD) is a genetically heterogeneous, autosomal-recessive disorder, characterized by oto-sino-pulmonary disease and situs abnormalities. PCD-causing mutations have been identified in 14 genes, but they collectively account for only ~60% of all PCD. To identify mutations that cause PCD, we performed exome sequencing on six unrelated probands with ciliary outer dynein arm (ODA) defects. Mutations in CCDC114, an ortholog of the Chlamydomonas reinhardtii motility gene DCC2, were identified in a family with two affected siblings. Sanger sequencing of 67 additional individuals with PCD with ODA defects from 58 families revealed CCDC114 mutations in 4 individuals in 3 families. All 6 individuals with CCDC114 mutations had characteristic oto-sino-pulmonary disease, but none had situs abnormalities. In the remaining 5 individuals with PCD who underwent exome sequencing, we identified mutations in two genes (DNAI2, DNAH5) known to cause PCD, including an Ashkenazi Jewish founder mutation in DNAI2. These results revealed that mutations in CCDC114 are a cause of ciliary dysmotility and PCD and further demonstrate the utility of exome sequencing to identify genetic causes in heterogeneous recessive disorders.

Mutations in KCTD1 cause scalp-ear-nipple syndrome.

Scalp-ear-nipple (SEN) syndrome is a rare, autosomal-dominant disorder characterized by cutis aplasia of the scalp; minor anomalies of the external ears, digits, and nails; and malformations of the breast. We used linkage analysis and exome sequencing of a multiplex family affected by SEN syndrome to identify potassium-channel tetramerization-domain-containing 1 (KCTD1) mutations that cause SEN syndrome. Evaluation of a total of ten families affected by SEN syndrome revealed KCTD1 missense mutations in each family tested. All of the mutations occurred in a KCTD1 region encoding a highly conserved bric-a-brac, tram track, and broad complex (BTB) domain that is required for transcriptional repressor activity. KCTD1 inhibits the transactivation of the transcription factor AP-2α (TFAP2A) via its BTB domain, and mutations in TFAP2A cause cutis aplasia in individuals with branchiooculofacial syndrome (BOFS), suggesting a potential overlap in the pathogenesis of SEN syndrome and BOFS. The identification of KCTD1 mutations in SEN syndrome reveals a role for this BTB-domain-containing transcriptional repressor during ectodermal development.

Mutations in SPAG1 cause primary ciliary dyskinesia associated with defective outer and inner dynein arms.

Primary ciliary dyskinesia (PCD) is a genetically heterogeneous, autosomal-recessive disorder, characterized by oto-sino-pulmonary disease and situs abnormalities. PCD-causing mutations have been identified in 20 genes, but collectively they account for only ∼65% of all PCDs. To identify mutations in additional genes that cause PCD, we performed exome sequencing on three unrelated probands with ciliary outer and inner dynein arm (ODA+IDA) defects. Mutations in SPAG1 were identified in one family with three affected siblings. Further screening of SPAG1 in 98 unrelated affected individuals (62 with ODA+IDA defects, 35 with ODA defects, 1 without available ciliary ultrastructure) revealed biallelic loss-of-function mutations in 11 additional individuals (including one sib-pair). All 14 affected individuals with SPAG1 mutations had a characteristic PCD phenotype, including 8 with situs abnormalities. Additionally, all individuals with mutations who had defined ciliary ultrastructure had ODA+IDA defects. SPAG1 was present in human airway epithelial cell lysates but was not present in isolated axonemes, and immunofluorescence staining showed an absence of ODA and IDA proteins in cilia from an affected individual, thus indicating that SPAG1 probably plays a role in the cytoplasmic assembly and/or trafficking of the axonemal dynein arms. Zebrafish morpholino studies of spag1 produced cilia-related phenotypes previously reported for PCD-causing mutations in genes encoding cytoplasmic proteins. Together, these results demonstrate that mutations in SPAG1 cause PCD with ciliary ODA+IDA defects and that exome sequencing is useful to identify genetic causes of heterogeneous recessive disorders.

The million mutation project: a new approach to genetics in Caenorhabditis elegans.

We have created a library of 2007 mutagenized Caenorhabditis elegans strains, each sequenced to a target depth of 15-fold coverage, to provide the research community with mutant alleles for each of the worm's more than 20,000 genes. The library contains over 800,000 unique single nucleotide variants (SNVs) with an average of eight nonsynonymous changes per gene and more than 16,000 insertion/deletion (indel) and copy number changes, providing an unprecedented genetic resource for this multicellular organism. To supplement this collection, we also sequenced 40 wild isolates, identifying more than 630,000 unique SNVs and 220,000 indels. Comparison of the two sets demonstrates that the mutant collection has a much richer array of both nonsense and missense mutations than the wild isolate set. We also find a wide range of rDNA and telomere repeat copy number in both sets. Scanning the mutant collection for molecular phenotypes reveals a nonsense suppressor as well as strains with higher levels of indels that harbor mutations in DNA repair genes and strains with abundant males associated with him mutations. All the strains are available through the Caenorhabditis Genetics Center and all the sequence changes have been deposited in WormBase and are available through an interactive website.

Improving performance of multigene panels for genomic analysis of cancer predisposition.

PURPOSE: Screening multiple genes for inherited cancer predisposition expands opportunities for cancer prevention; however, reports of variants of uncertain significance (VUS) may limit clinical usefulness. We used an expert-driven approach, exploiting all available information, to evaluate multigene panels for inherited cancer predisposition in a clinical series that included multiple cancer types and complex family histories. METHODS: For 1,462 sequential patients referred for testing by BROCA or ColoSeq multigene panels, genomic DNA was sequenced and variants were interpreted by multiple experts using International Agency for Research on Cancer guidelines and incorporating evolutionary conservation, known and predicted variant consequences, and personal and family cancer history. Diagnostic yield was evaluated for various presenting conditions and family-history profiles. RESULTS: Of 1,462 patients, 12% carried damaging mutations in established cancer genes. Diagnostic yield varied by clinical presentation. Actionable results were identified for 13% of breast and colorectal cancer patients and for 4% of cancer-free subjects, based on their family histories of cancer. Incidental findings explaining cancer in neither the patient nor the family were present in 1.7% of subjects. Less than 1% of patients carried VUS in BRCA1 or BRCA2. For all genes combined, initial reports contained VUS for 10.5% of patients, which declined to 7.5% of patients after reclassification based on additional information. CONCLUSIONS: Individualized interpretation of gene panels is a complex medical activity. Interpretation by multiple experts in the context of personal and family histories maximizes actionable results and minimizes reports of VUS.Genet Med 18 10, 974-981.

Whole-exome capture and sequencing identifies HEATR2 mutation as a cause of primary ciliary dyskinesia.

Motile cilia are essential components of the mucociliary escalator and are central to respiratory-tract host defenses. Abnormalities in these evolutionarily conserved organelles cause primary ciliary dyskinesia (PCD). Despite recent strides characterizing the ciliome and sensory ciliopathies through exploration of the phenotype-genotype associations in model organisms, the genetic bases of most cases of PCD remain elusive. We identified nine related subjects with PCD from geographically dispersed Amish communities and performed exome sequencing of two affected individuals and their unaffected parents. A single autosomal-recessive nonsynonymous missense mutation was identified in HEATR2, an uncharacterized gene that belongs to a family not previously associated with ciliary assembly or function. Airway epithelial cells isolated from PCD-affected individuals had markedly reduced HEATR2 levels, absent dynein arms, and loss of ciliary beating. MicroRNA-mediated silencing of the orthologous gene in Chlamydomonas reinhardtii resulted in absent outer dynein arms, reduced flagellar beat frequency, and decreased cell velocity. These findings were recapitulated by small hairpin RNA-mediated knockdown of HEATR2 in airway epithelial cells from unaffected donors. Moreover, immunohistochemistry studies in human airway epithelial cells showed that HEATR2 was localized to the cytoplasm and not in cilia, which suggests a role in either dynein arm transport or assembly. The identification of HEATR2 contributes to the growing number of genes associated with PCD identified in both individuals and model organisms and shows that exome sequencing in family studies facilitates the discovery of novel disease-causing gene mutations.

Targeted capture and massively parallel sequencing of 12 human exomes.

Genome-wide association studies suggest that common genetic variants explain only a modest fraction of heritable risk for common diseases, raising the question of whether rare variants account for a significant fraction of unexplained heritability. Although DNA sequencing costs have fallen markedly, they remain far from what is necessary for rare and novel variants to be routinely identified at a genome-wide scale in large cohorts. We have therefore sought to develop second-generation methods for targeted sequencing of all protein-coding regions ('exomes'), to reduce costs while enriching for discovery of highly penetrant variants. Here we report on the targeted capture and massively parallel sequencing of the exomes of 12 humans. These include eight HapMap individuals representing three populations, and four unrelated individuals with a rare dominantly inherited disorder, Freeman-Sheldon syndrome (FSS). We demonstrate the sensitive and specific identification of rare and common variants in over 300 megabases of coding sequence. Using FSS as a proof-of-concept, we show that candidate genes for Mendelian disorders can be identified by exome sequencing of a small number of unrelated, affected individuals. This strategy may be extendable to diseases with more complex genetics through larger sample sizes and appropriate weighting of non-synonymous variants by predicted functional impact.

Gorilla genome structural variation reveals evolutionary parallelisms with chimpanzee.

Structural variation has played an important role in the evolutionary restructuring of human and great ape genomes. Recent analyses have suggested that the genomes of chimpanzee and human have been particularly enriched for this form of genetic variation. Here, we set out to assess the extent of structural variation in the gorilla lineage by generating 10-fold genomic sequence coverage from a western lowland gorilla and integrating these data into a physical and cytogenetic framework of structural variation. We discovered and validated over 7665 structural changes within the gorilla lineage, including sequence resolution of inversions, deletions, duplications, and mobile element insertions. A comparison with human and other ape genomes shows that the gorilla genome has been subjected to the highest rate of segmental duplication. We show that both the gorilla and chimpanzee genomes have experienced independent yet convergent patterns of structural mutation that have not occurred in humans, including the formation of subtelomeric heterochromatic caps, the hyperexpansion of segmental duplications, and bursts of retroviral integrations. Our analysis suggests that the chimpanzee and gorilla genomes are structurally more derived than either orangutan or human genomes.

Microsatellite instability detection by next generation sequencing.

BACKGROUND: Microsatellite instability (MSI) is a useful phenotype in cancer diagnosis and prognosis. Nevertheless, methods to detect MSI status from next generation DNA sequencing (NGS) data are underdeveloped. METHODS: We developed an approach to detect the MSI phenotype using NGS (mSINGS). The method was used to evaluate mononucleotide microsatellite loci that were incidentally sequenced after targeted gene enrichment and could be applied to gene or exome capture panels designed for other purposes. For each microsatellite locus, the number of differently sized repeats in experimental samples were quantified and compared to a population of normal controls. Loci were considered unstable if the experimental number of repeats was statistically greater than in the control population. MSI status was determined by the fraction of unstable microsatellite loci. RESULTS: We examined data from 324 samples generated using targeted gene capture assays of 3 different sizes, ranging from a 0.85-Mb to a 44-Mb exome design and incorporating from 15 to 2957 microsatellite markers. When we compared mSING results to MSI-PCR as a gold standard for 108 cases, we found the approach to be both diagnostically sensitive (range of 96.4% to 100% across 3 panels) and specific (range of 97.2% to 100%) for determining MSI status. The fraction of unstable microsatellite markers calculated from sequencing data correlated with the number of unstable loci detected by conventional MSI-PCR testing. CONCLUSIONS: NGS data can enable highly accurate detection of MSI, even from limited capture designs. This novel approach offers several advantages over existing PCR-based methods.

Parallel, tag-directed assembly of locally derived short sequence reads.

We demonstrate subassembly, an in vitro library construction method that extends the utility of short-read sequencing platforms to applications requiring long, accurate reads. A long DNA fragment library is converted to a population of nested sublibraries, and a tag sequence directs grouping of short reads derived from the same long fragment, enabling localized assembly of long fragment sequences. Subassembly may facilitate accurate de novo genome assembly and metagenome sequencing.