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.

Applying Ancestry and Sex Computation as a Quality Control Tool in Targeted Next-Generation Sequencing.

OBJECTIVES: To apply techniques for ancestry and sex computation from next-generation sequencing (NGS) data as an approach to confirm sample identity and detect sample processing errors. METHODS: We combined a principal component analysis method with k-nearest neighbors classification to compute the ancestry of patients undergoing NGS testing. By combining this calculation with X chromosome copy number data, we determined the sex and ancestry of patients for comparison with self-report. We also modeled the sensitivity of this technique in detecting sample processing errors. RESULTS: We applied this technique to 859 patient samples with reliable self-report data. Our k-nearest neighbors ancestry screen had an accuracy of 98.7% for patients reporting a single ancestry. Visual inspection of principal component plots was consistent with self-report in 99.6% of single-ancestry and mixed-ancestry patients. Our model demonstrates that approximately two-thirds of potential sample swaps could be detected in our patient population using this technique. CONCLUSIONS: Patient ancestry can be estimated from NGS data incidentally sequenced in targeted panels, enabling an inexpensive quality control method when coupled with patient self-report.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

"Mandibulofacial dysostosis with microcephaly" caused by EFTUD2 mutations: expanding the phenotype.

Heterozygous mutations in the EFTUD2 were identified in 12 individuals with a rare sporadic craniofacial condition termed Mandibulofacial dysostosis with microcephaly (MIM 610536). We present clinical and radiographic features of three additional patients with de novo heterozygous mutations in EFTUD2. Although clinical features overlap with findings of the original report (choanal atresia, cleft palate, maxillary and mandibular hypoplasia, and microtia), microcephaly was present in two of three patients and cognitive impairment was milder in those with head circumference proportional to height. Our cases expand the phenotypic spectrum to include epibulbar dermoids and zygomatic arch clefting. We suggest that craniofacial computed tomography studies to assess cleft of zygomatic arch may assist in making this diagnosis. We recommend consideration of EFTUD2 testing in individuals with features of oculo-auriculo-vertebral spectrum and bilateral microtia, or individuals with atypical CHARGE syndrome who do not have a CHD7 mutation, particularly those with a zygomatic arch cleft. The absence of microcephaly in one patient indicates that it is a highly variable phenotypic feature.

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.

Multiplex targeted sequencing identifies recurrently mutated genes in autism spectrum disorders.

Exome sequencing studies of autism spectrum disorders (ASDs) have identified many de novo mutations but few recurrently disrupted genes. We therefore developed a modified molecular inversion probe method enabling ultra-low-cost candidate gene resequencing in very large cohorts. To demonstrate the power of this approach, we captured and sequenced 44 candidate genes in 2446 ASD probands. We discovered 27 de novo events in 16 genes, 59% of which are predicted to truncate proteins or disrupt splicing. We estimate that recurrent disruptive mutations in six genes-CHD8, DYRK1A, GRIN2B, TBR1, PTEN, and TBL1XR1-may contribute to 1% of sporadic ASDs. Our data support associations between specific genes and reciprocal subphenotypes (CHD8-macrocephaly and DYRK1A-microcephaly) and replicate the importance of a ß-catenin-chromatin-remodeling network to ASD etiology.

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.