Lessons learned from additional research analyses of unsolved clinical exome cases.

BACKGROUND: Given the rarity of most single-gene Mendelian disorders, concerted efforts of data exchange between clinical and scientific communities are critical to optimize molecular diagnosis and novel disease gene discovery. METHODS: We designed and implemented protocols for the study of cases for which a plausible molecular diagnosis was not achieved in a clinical genomics diagnostic laboratory (i.e. unsolved clinical exomes). Such cases were recruited to a research laboratory for further analyses, in order to potentially: (1) accelerate novel disease gene discovery; (2) increase the molecular diagnostic yield of whole exome sequencing (WES); and (3) gain insight into the genetic mechanisms of disease. Pilot project data included 74 families, consisting mostly of parent-offspring trios. Analyses performed on a research basis employed both WES from additional family members and complementary bioinformatics approaches and protocols. RESULTS: Analysis of all possible modes of Mendelian inheritance, focusing on both single nucleotide variants (SNV) and copy number variant (CNV) alleles, yielded a likely contributory variant in 36% (27/74) of cases. If one includes candidate genes with variants identified within a single family, a potential contributory variant was identified in a total of ~51% (38/74) of cases enrolled in this pilot study. The molecular diagnosis was achieved in 30/63 trios (47.6%). Besides this, the analysis workflow yielded evidence for pathogenic variants in disease-associated genes in 4/6 singleton cases (66.6%), 1/1 multiplex family involving three affected siblings, and 3/4 (75%) quartet families. Both the analytical pipeline and the collaborative efforts between the diagnostic and research laboratories provided insights that allowed recent disease gene discoveries (PURA, TANGO2, EMC1, GNB5, ATAD3A, and MIPEP) and increased the number of novel genes, defined in this study as genes identified in more than one family (DHX30 and EBF3). CONCLUSION: An efficient genomics pipeline in which clinical sequencing in a diagnostic laboratory is followed by the detailed reanalysis of unsolved cases in a research environment, supplemented with WES data from additional family members, and subject to adjuvant bioinformatics analyses including relaxed variant filtering parameters in informatics pipelines, can enhance the molecular diagnostic yield and provide mechanistic insights into Mendelian disorders. Implementing these approaches requires collaborative clinical molecular diagnostic and research efforts.

MIPEP recessive variants cause a syndrome of left ventricular non-compaction, hypotonia, and infantile death.

BACKGROUND: Mitochondrial presequence proteases perform fundamental functions as they process about 70 % of all mitochondrial preproteins that are encoded in the nucleus and imported posttranslationally. The mitochondrial intermediate presequence protease MIP/Oct1, which carries out precursor processing, has not yet been established to have a role in human disease. METHODS: Whole exome sequencing was performed on four unrelated probands with left ventricular non-compaction (LVNC), developmental delay (DD), seizures, and severe hypotonia. Proposed pathogenic variants were confirmed by Sanger sequencing or array comparative genomic hybridization. Functional analysis of the identified MIP variants was performed using the model organism Saccharomyces cerevisiae as the protein and its functions are highly conserved from yeast to human. RESULTS: Biallelic single nucleotide variants (SNVs) or copy number variants (CNVs) in MIPEP, which encodes MIP, were present in all four probands, three of whom had infantile/childhood death. Two patients had compound heterozygous SNVs (p.L582R/p.L71Q and p.E602*/p.L306F) and one patient from a consanguineous family had a homozygous SNV (p.K343E). The fourth patient, identified through the GeneMatcher tool, a part of the Matchmaker Exchange Project, was found to have inherited a paternal SNV (p.H512D) and a maternal CNV (1.4-Mb deletion of 13q12.12) that includes MIPEP. All amino acids affected in the patients' missense variants are highly conserved from yeast to human and therefore S. cerevisiae was employed for functional analysis (for p.L71Q, p.L306F, and p.K343E). The mutations p.L339F (human p.L306F) and p.K376E (human p.K343E) resulted in a severe decrease of Oct1 protease activity and accumulation of non-processed Oct1 substrates and consequently impaired viability under respiratory growth conditions. The p.L83Q (human p.L71Q) failed to localize to the mitochondria. CONCLUSIONS: Our findings reveal for the first time the role of the mitochondrial intermediate peptidase in human disease. Loss of MIP function results in a syndrome which consists of LVNC, DD, seizures, hypotonia, and cataracts. Our approach highlights the power of data exchange and the importance of an interrelationship between clinical and research efforts for disease gene discovery.

Secondary findings and carrier test frequencies in a large multiethnic sample.

BACKGROUND: Besides its growing importance in clinical diagnostics and understanding the genetic basis of Mendelian and complex diseases, whole exome sequencing (WES) is a rich source of additional information of potential clinical utility for physicians, patients and their families. We analyzed the frequency and nature of single nucleotide variants (SNVs) considered secondary findings and recessive disease allele carrier status in the exomes of 8554 individuals from a large, randomly sampled cohort study and 2514 patients from a study of presumed Mendelian disease having undergone WES. METHODS: We used the same sequencing platform and data processing pipeline to analyze all samples and characterized the distributions of reported pathogenic (ClinVar, Human Gene Mutation Database (HGMD)) and predicted deleterious variants in the pre-specified American College of Medical Genetics and Genomics (ACMG) secondary findings and recessive disease genes in different ethnic groups. RESULTS: In the 56 ACMG secondary findings genes, the average number of predicted deleterious variants per individual was 0.74, and the mean number of ClinVar reported pathogenic variants was 0.06. We observed an average of 10 deleterious and 0.78 ClinVar reported pathogenic variants per individual in 1423 autosomal recessive disease genes. By repeatedly sampling pairs of exomes, 0.5 % of the randomly generated couples were at 25 % risk of having an affected offspring for an autosomal recessive disorder based on the ClinVar variants. CONCLUSIONS: By investigating reported pathogenic and novel, predicted deleterious variants we estimated the lower and upper limits of the population fraction for which exome sequencing may reveal additional medically relevant information. We suggest that the observed wide range for the lower and upper limits of these frequency numbers will be gradually reduced due to improvement in classification databases and prediction algorithms.

COPA mutations impair ER-Golgi transport and cause hereditary autoimmune-mediated lung disease and arthritis.

Unbiased genetic studies have uncovered surprising molecular mechanisms in human cellular immunity and autoimmunity. We performed whole-exome sequencing and targeted sequencing in five families with an apparent mendelian syndrome of autoimmunity characterized by high-titer autoantibodies, inflammatory arthritis and interstitial lung disease. We identified four unique deleterious variants in the COPA gene (encoding coatomer subunit α) affecting the same functional domain. Hypothesizing that mutant COPA leads to defective intracellular transport via coat protein complex I (COPI), we show that COPA variants impair binding to proteins targeted for retrograde Golgi-to-ER transport. Additionally, expression of mutant COPA results in ER stress and the upregulation of cytokines priming for a T helper type 17 (TH17) response. Patient-derived CD4(+) T cells also demonstrate significant skewing toward a TH17 phenotype that is implicated in autoimmunity. Our findings uncover an unexpected molecular link between a vesicular transport protein and a syndrome of autoimmunity manifested by lung and joint disease.

DVL1 frameshift mutations clustering in the penultimate exon cause autosomal-dominant Robinow syndrome.

Robinow syndrome is a genetically heterogeneous disorder characterized by mesomelic limb shortening, genital hypoplasia, and distinctive facial features and for which both autosomal-recessive and autosomal-dominant inheritance patterns have been described. Causative variants in the non-canonical signaling gene WNT5A underlie a subset of autosomal-dominant Robinow syndrome (DRS) cases, but most individuals with DRS remain without a molecular diagnosis. We performed whole-exome sequencing in four unrelated DRS-affected individuals without coding mutations in WNT5A and found heterozygous DVL1 exon 14 mutations in three of them. Targeted Sanger sequencing in additional subjects with DRS uncovered DVL1 exon 14 mutations in five individuals, including a pair of monozygotic twins. In total, six distinct frameshift mutations were found in eight subjects, and all were heterozygous truncating variants within the penultimate exon of DVL1. In five families in which samples from unaffected parents were available, the variants were demonstrated to represent de novo mutations. All variant alleles are predicted to result in a premature termination codon within the last exon, escape nonsense-mediated decay (NMD), and most likely generate a C-terminally truncated protein with a distinct -1 reading-frame terminus. Study of the transcripts extracted from affected subjects' leukocytes confirmed expression of both wild-type and variant alleles, supporting the hypothesis that mutant mRNA escapes NMD. Genomic variants identified in our study suggest that truncation of the C-terminal domain of DVL1, a protein hypothesized to have a downstream role in the Wnt-5a non-canonical pathway, is a common cause of DRS.

Mutations in COL27A1 cause Steel syndrome and suggest a founder mutation effect in the Puerto Rican population.

Osteochondrodysplasias represent a large group of developmental structural disorders that can be caused by mutations in a variety of genes responsible for chondrocyte development, differentiation, mineralization and early ossification. The application of whole-exome sequencing to disorders apparently segregating as Mendelian traits has proven to be an effective approach to disease gene identification for conditions with unknown molecular etiology. We identified a homozygous missense variant p.(Gly697Arg) in COL27A1, in a family with Steel syndrome and no consanguinity. Interestingly, the identified variant seems to have arisen as a founder mutation in the Puerto Rican population.

A drosophila genetic resource of mutants to study mechanisms underlying human genetic diseases.

Invertebrate model systems are powerful tools for studying human disease owing to their genetic tractability and ease of screening. We conducted a mosaic genetic screen of lethal mutations on the Drosophila X chromosome to identify genes required for the development, function, and maintenance of the nervous system. We identified 165 genes, most of whose function has not been studied in vivo. In parallel, we investigated rare variant alleles in 1,929 human exomes from families with unsolved Mendelian disease. Genes that are essential in flies and have multiple human homologs were found to be likely to be associated with human diseases. Merging the human data sets with the fly genes allowed us to identify disease-associated mutations in six families and to provide insights into microcephaly associated with brain dysgenesis. This bidirectional synergism between fly genetics and human genomics facilitates the functional annotation of evolutionarily conserved genes involved in human health.

De novo truncating mutations in AHDC1 in individuals with syndromic expressive language delay, hypotonia, and sleep apnea.

Clinical whole-exome sequencing (WES) for identification of mutations leading to Mendelian disease has been offered to the medical community since 2011. Clinically undiagnosed neurological disorders are the most frequent basis for test referral, and currently, approximately 25% of such cases are diagnosed at the molecular level. To date, there are approximately 4,000 "known" disease-associated loci, and many are associated with striking dysmorphic features, making genotype-phenotype correlations relatively straightforward. A significant fraction of cases, however, lack characteristic dysmorphism or clinical pathognomonic traits and are dependent upon molecular tests for definitive diagnoses. Further, many molecular diagnoses are guided by recent gene-disease association discoveries. Hence, there is a critical interplay between clinical testing and research leading to gene-disease association discovery. Here, we describe four probands, all of whom presented with hypotonia, intellectual disability, global developmental delay, and mildly dysmorphic facial features. Three of the four also had sleep apnea. Each was a simplex case without a remarkable family history. Using WES, we identified AHDC1 de novo truncating mutations that most likely cause this genetic syndrome.

Clinical utility of whole-exome sequencing in rare diseases: Galactosialidosis.

Rare genetic disorders can go undiagnosed for years as the entire spectrum of phenotypic variation is not well characterized given the reduced number of patients reported in the literature and the low frequency at which these occur. Moreover, the current paradigm for clinical diagnostics defines disease diagnosis by a specified spectrum of phenotypic findings; when such parameters are either missing, or other findings not usually observed are seen, the phenotype driven approach to diagnosis may result in a specific etiological diagnosis not even being considered within the differential diagnosis. The novel implementation of genomic sequencing approaches to investigate rare genetic disorders is allowing not only the discovery of new genes, but also the phenotypic expansion of known Mendelian genetic disorders. Here we report the detailed clinical assessment of a patient with a rare genetic disorder with undefined molecular diagnosis. We applied whole-exome sequencing to this patient and unaffected parents in order to identify the molecular cause of her disorder. We identified compound heterozygous mutations in the CTSA gene, responsible for causing galactosialidosis; the molecular diagnosis was further confirmed by biochemical studies. This report expands on the clinical spectrum of this rare lysosomal disorder and exemplifies how genomic approaches are further elucidating the characterization and understanding of genetic diseases.

Whole-exome sequencing reveals GPIHBP1 mutations in infantile colitis with severe hypertriglyceridemia.

Severe congenital hypertriglyceridemia (HTG) is a rare disorder caused by mutations in genes affecting lipoprotein lipase (LPL) activity. Here we report a 5-week-old Hispanic girl with severe HTG (12,031 mg/dL, normal limit 150 mg/dL) who presented with the unusual combination of lower gastrointestinal bleeding and milky plasma. Initial colonoscopy was consistent with colitis, which resolved with reduction of triglycerides. After negative sequencing of the LPL gene, whole-exome sequencing revealed novel compound heterozygous mutations in GPIHBP1. Our study broadens the phenotype of GPIHBP1-associated HTG, reinforces the effectiveness of whole-exome sequencing in Mendelian diagnoses, and implicates triglycerides in gastrointestinal mucosal injury.

Heterozygous de novo and inherited mutations in the smooth muscle actin (ACTG2) gene underlie megacystis-microcolon-intestinal hypoperistalsis syndrome.

Megacystis-microcolon-intestinal hypoperistalsis syndrome (MMIHS) is a rare disorder of enteric smooth muscle function affecting the intestine and bladder. Patients with this severe phenotype are dependent on total parenteral nutrition and urinary catheterization. The cause of this syndrome has remained a mystery since Berdon's initial description in 1976. No genes have been clearly linked to MMIHS. We used whole-exome sequencing for gene discovery followed by targeted Sanger sequencing in a cohort of patients with MMIHS and intestinal pseudo-obstruction. We identified heterozygous ACTG2 missense variants in 15 unrelated subjects, ten being apparent de novo mutations. Ten unique variants were detected, of which six affected CpG dinucleotides and resulted in missense mutations at arginine residues, perhaps related to biased usage of CpG containing codons within actin genes. We also found some of the same heterozygous mutations that we observed as apparent de novo mutations in MMIHS segregating in families with intestinal pseudo-obstruction, suggesting that ACTG2 is responsible for a spectrum of smooth muscle disease. ACTG2 encodes γ2 enteric actin and is the first gene to be clearly associated with MMIHS, suggesting an important role for contractile proteins in enteric smooth muscle disease.

Mutations in VRK1 associated with complex motor and sensory axonal neuropathy plus microcephaly.

IMPORTANCE: Patients with rare diseases and complex clinical presentations represent a challenge for clinical diagnostics. Genomic approaches are allowing the identification of novel variants in genes for very rare disorders, enabling a molecular diagnosis. Genomics is also revealing a phenotypic expansion whereby the full spectrum of clinical expression conveyed by mutant alleles at a locus can be better appreciated. OBJECTIVE: To elucidate the molecular cause of a complex neuropathy phenotype in 3 patients by applying genomic sequencing strategies. DESIGN, SETTING, AND PARTICIPANTS: Three affected individuals from 2 unrelated families presented with a complex neuropathy phenotype characterized by axonal sensorimotor neuropathy and microcephaly. They were recruited into the Centers for Mendelian Genomics research program to identify the molecular cause of their phenotype. Whole-genome, targeted whole-exome sequencing, and high-resolution single-nucleotide polymorphism arrays were performed in genetics clinics of tertiary care pediatric hospitals and biomedical research institutions. MAIN OUTCOMES AND MEASURES: Whole-genome and whole-exome sequencing identified the variants responsible for the patients' clinical phenotype. RESULTS: We identified compound heterozygous alleles in 2 affected siblings from 1 family and a homozygous nonsense variant in the third unrelated patient in the vaccinia-related kinase 1 gene (VRK1). In the latter subject, we found a common haplotype on which the nonsense mutation occurred and that segregates in the Ashkenazi Jewish population. CONCLUSIONS AND RELEVANCE: We report the identification of disease-causing alleles in 3 children from 2 unrelated families with a previously uncharacterized complex axonal motor and sensory neuropathy accompanied by severe nonprogressive microcephaly and cerebral dysgenesis. Our data raise the question of whether VRK1 mutations disturb cell cycle progression and may result in apoptosis of cells in the nervous system. The application of unbiased genomic approaches allows the identification of potentially pathogenic mutations in unsuspected genes in highly genetically heterogeneous and uncharacterized neurological diseases.