Genetic diagnosis of Mendelian disorders via RNA sequencing
Laura S Kremer, Daniel M Bader, Christian Mertes, Robert Kopajtich, Garwin Pichler, Arcangela Iuso, Tobias B Haack, Elisabeth Graf, Thomas Schwarzmayr, Caterina Terrile, Eliška Koňaříková, Birgit Repp, Gabi Kastenmüller, Jerzy Adamski, Peter Lichtner, Christoph Leonhardt, Benoit Funalot, Alice Donati, Valeria Tiranti, Anne Lombes, Claude Jardel, Dieter Gläser, Robert W Taylor, Daniele Ghezzi, Johannes A Mayr, Agnes Rötig, Peter Freisinger, Felix Distelmaier, Tim M Strom, Thomas Meitinger, Julien Gagneur, Holger Prokisch, Laura S Kremer, Daniel M Bader, Christian Mertes, Robert Kopajtich, Garwin Pichler, Arcangela Iuso, Tobias B Haack, Elisabeth Graf, Thomas Schwarzmayr, Caterina Terrile, Eliška Koňaříková, Birgit Repp, Gabi Kastenmüller, Jerzy Adamski, Peter Lichtner, Christoph Leonhardt, Benoit Funalot, Alice Donati, Valeria Tiranti, Anne Lombes, Claude Jardel, Dieter Gläser, Robert W Taylor, Daniele Ghezzi, Johannes A Mayr, Agnes Rötig, Peter Freisinger, Felix Distelmaier, Tim M Strom, Thomas Meitinger, Julien Gagneur, Holger Prokisch
Abstract
Across a variety of Mendelian disorders, ∼50-75% of patients do not receive a genetic diagnosis by exome sequencing indicating disease-causing variants in non-coding regions. Although genome sequencing in principle reveals all genetic variants, their sizeable number and poorer annotation make prioritization challenging. Here, we demonstrate the power of transcriptome sequencing to molecularly diagnose 10% (5 of 48) of mitochondriopathy patients and identify candidate genes for the remainder. We find a median of one aberrantly expressed gene, five aberrant splicing events and six mono-allelically expressed rare variants in patient-derived fibroblasts and establish disease-causing roles for each kind. Private exons often arise from cryptic splice sites providing an important clue for variant prioritization. One such event is found in the complex I assembly factor TIMMDC1 establishing a novel disease-associated gene. In conclusion, our study expands the diagnostic tools for detecting non-exonic variants and provides examples of intronic loss-of-function variants with pathological relevance.
Conflict of interest statement
The authors declare no competing financial interests.
Figures
References
- Wortmann S. B., Koolen D. A., Smeitink J. A., van den Heuvel L. & Rodenburg R. J. Whole exome sequencing of suspected mitochondrial patients in clinical practice. J. Inherit. Metab. Dis. 38, 437–443 (2015).
- 1000 Genomes Project Consortium. et al.. A global reference for human genetic variation. Nature 526, 68–74 (2015).
- Sudmant P. H. et al.. An integrated map of structural variation in 2,504 human genomes. Nature 526, 75–81 (2015).
- Taylor J. C. et al.. Factors influencing success of clinical genome sequencing across a broad spectrum of disorders. Nat. Genet. 47, 717–726 (2015).
- Li X. et al.. The impact of rare variation on gene expression across tissues. bioRxiv doi:; DOI: 10.1101/074443 (2016).
- Zeng Y. et al.. Aberrant gene expression in humans. PLoS Genet. 11, 1–20 (2015).
- Guan J. et al.. Exploiting aberrant mRNA expression in autism for gene discovery and diagnosis. Hum. Genet. 135, 1–15 (2016).
- Zhao J. et al.. A burden of rare variants associated with extremes of gene expression in human peripheral blood. Am. J. Hum. Genet. 98, 299–309 (2016).
- Albers C. A. Compound inheritance of a low-frequency regulatory SNP and a rare null mutation in exon-junction complex subunit RBM8A causes TAR syndrome. Nat. Genet. 44, 435–439 S1-2 (2012).
- Reinius B. & Sandberg R. Random monoallelic expression of autosomal genes: stochastic transcription and allele-level regulation. Nat. Rev. Genet. 16, 653–664 (2015).
- Eckersley-Maslin M. A. & Spector D. L. Random monoallelic expression: regulating gene expression one allele at a time. Trends Genet. 30, 237–244 (2014).
- Tazi J., Bakkour N. & Stamm S. Alternative splicing and disease. Biochim. Biophys. Acta 1792, 14–26 2009.
- Scotti M. M. & Swanson M. S. RNA mis-splicing in disease. Nat. Rev. Genet. 17, 19–32 (2015).
- Singh R. K. & Cooper T. A. Pre-mRNA splicing in disease and therapeutics. Trends Mol. Med. 18, 472–482 (2012).
- Xiong H. Y. et al.. The human splicing code reveals new insights into the genetic determinants of disease. Science 347, 1254806–1254806 (2015).
- Muntoni F., Torelli S. & Ferlini A. Dystrophin and mutations: one gene, several proteins, multiple phenotypes. Lancet Neurol. 2, 731–740 (2003).
- Gonorazky H. et al.. RNAseq analysis for the diagnosis of muscular dystrophy. Ann. Clin. Transl. Neurol. 3, 55–60 (2016).
- Morel C. F. et al.. A LMNA splicing mutation in two sisters with severe dunnigan-type familial partial lipodystrophy type 2. J. Clin. Endocrinol. Metab. 91, 2689–2695 (2006).
- Qu Y. et al.. A rare variant (c.863G>T) in exon 7 of SMN1 disrupts mRNA splicing and is responsible for spinal muscular atrophy. Eur. J. Hum. Genet. 24, 864–870 (2016).
- Gorman G. S. et al.. Mitochondrial diseases. Nat. Rev. Dis. Primer 2, 16080 (2016).
- Elstner M., Andreoli C. & Ahting U. MitoP2: an integrative tool for the analysis of the mitochondrial proteome. Mol. Biotechnol. 40, 306–315 (2008).
- Mayr J. A. et al.. Spectrum of combined respiratory chain defects. J. Inherit. Metab. Dis. 38, 629–640 (2015).
- Haack T. B. et al.. ELAC2 mutations cause a mitochondrial RNA processing defect associated with hypertrophic cardiomyopathy. Am. J. Hum. Genet. 93, 211–223 (2013).
- Haack T. B. et al.. Exome sequencing identifies ACAD9 mutations as a cause of complex I deficiency. Nat. Genet. 42, 1131–1134 (2010).
- Li Y. I. et al.. RNA splicing is a primary link between genetic variation and disease. Science 352, 600–604 (2016).
- Zhang B. et al.. Proteogenomic characterization of human colon and rectal cancer. Nature 513, 382–387 (2014).
- Liu Y., Beyer A. & Aebersold R. On the dependency of cellular protein levels on mRNA abundance. Cell 165, 535–550 (2016).
- Lee K. K., Shimoji M., Hossain Q. S., Sunakawa H. & Aniya Y. Novel function of glutathione transferase in rat liver mitochondrial membrane: role for cytochrome c release from mitochondria. Toxicol. Appl. Pharmacol. 232, 109–118 (2008).
- Holzerova E. et al.. Human thioredoxin 2 deficiency impairs mitochondrial redox homeostasis and causes early-onset neurodegeneration. Brain 139, 346–354 (2016).
- Guarani V. et al.. TIMMDC1/C3orf1 functions as a membrane-embedded mitochondrial complex I assembly factor through association with the MCIA complex. Mol. Cell. Biol. 34, 847–861 (2014).
- Andrews B., Carroll J., Ding S., Fearnley I. M. & Walker J. E. Assembly factors for the membrane arm of human complex I. Proc. Natl Acad. Sci. USA 110, 18934–18939 (2013).
- Pervouchine D. D., Knowles D. G. & Guig R. Intron-centric estimation of alternative splicing from RNA-seq data. Bioinformatics 29, 273–274 (2013).
- Lek M. et al.. Analysis of protein-coding genetic variation in 60,706 humans. Nature 536, 285–291 (2016).
- Halperin T., Zheng B., Itzhaki H., Clarke A. K. & Adam Z. Plant mitochondria contain proteolytic and regulatory subunits of the ATP-dependent Clp protease. Plant Mol. Biol. 45, 461–468 (2001).
- Jenkinson E. M. et al.. Perrault syndrome is caused by recessive mutations in CLPP, encoding a mitochondrial ATP-dependent chambered protease. Am. J. Hum. Genet. 92, 605–613 (2013).
- Jenkinson E. M. et al.. Perrault syndrome: further evidence for genetic heterogeneity. J. Neurol. 259, 974–976 (2012).
- Szczepanowska K. et al.. CLPP coordinates mitoribosomal assembly through the regulation of ERAL1 levels. EMBO J. 35, 2566–2583 (2016).
- Piva F., Giulietti M., Burini A. B. & Principato G. SpliceAid 2: a database of human splicing factors expression data and RNA target motifs. Hum. Mutat. 33, 81–85 (2012).
- Dogan R. I., Getoor L., Wilbur W. J. & Mount S. M. SplicePort--an interactive splice-site analysis tool. Nucleic Acids Res. 35, W285–W291 (2007).
- Timmermans M. J. T. N. et al.. Why barcode? High-throughput multiplex sequencing of mitochondrial genomes for molecular systematics. Nucleic Acids Res. 38, e197–e197 (2010).
- Desmet F.-O. et al.. Human Splicing Finder: an online bioinformatics tool to predict splicing signals. Nucleic Acids Res. 37, e67 (2009).
- Yeo G., Hoon S., Venkatesh B. & Burge C. B. Variation in sequence and organization of splicing regulatory elements in vertebrate genes. Proc. Natl Acad. Sci. USA 101, 15700–15705 (2004).
- Burge C. & Karlin S. Prediction of complete gene structures in human genomic DNA. J. Mol. Biol. 268, 78–94 (1997).
- Kapustin Y. et al.. Cryptic splice sites and split genes. Nucleic Acids Res. 39, 5837–5844 (2011).
- Adams E. & Frank L. Metabolism of proline and the hydroxyprolines. Annu. Rev. Biochem. 49, 1005–1061 (1980).
- Baumgartner M. R. et al.. Hyperammonemia with reduced ornithine, citrulline, arginine and proline: a new inborn error caused by a mutation in the gene encoding delta(1)-pyrroline-5-carboxylate synthase. Hum. Mol. Genet. 9, 2853–2858 (2000).
- Fischer-Zirnsak B. et al.. Recurrent de novo mutations affecting residue Arg138 of pyrroline-5-carboxylate synthase cause a progeroid form of autosomal-dominant cutis laxa. Am. J. Hum. Genet. 97, 483–492 (2015).
- Coutelier M. et al.. Alteration of ornithine metabolism leads to dominant and recessive hereditary spastic paraplegia. Brain 138, 2191–2205 (2015).
- Sibley C. R., Blazquez L. & Ule J. Lessons from non-canonical splicing. Nat. Rev. Genet. 17, 407–421 (2016).
- Cummings B. B. et al.. Improving genetic diagnosis in Mendelian disease with transcriptome sequencing. Sci. Transl. Med 9, eaal5209 (2017).
- Gibson G. Human genetics. GTEx detects genetic effects. Science 348, 640–641 (2015).
- Vafai S. B. & Mootha V. K. Mitochondrial disorders as windows into an ancient organelle. Nature 491, 374–383 (2012).
- Gagneur J. et al.. Genotype-environment interactions reveal causal pathways that mediate genetic effects on phenotype. PLoS Genet. 9, e1003803 (2013).
- Mayr J. A. et al.. Lack of the mitochondrial protein acylglycerol kinase causes sengers syndrome. Am. J. Hum. Genet. 90, 314–320 (2012).
- Li H. & Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 25, 1754–1760 (2009).
- Li H. A statistical framework for SNP calling, mutation discovery, association mapping and population genetical parameter estimation from sequencing data. Bioinformatics 27, 2987–2993 (2011).
- Li H. et al.. The Sequence Alignment/Map format and SAMtools. Bioinformatics 25, 2078–2079 (2009).
- McLaren W. et al.. Deriving the consequences of genomic variants with the Ensembl API and SNP Effect Predictor. Bioinformatics 26, 2069–2070 (2010).
- Eilbeck K. et al.. The Sequence Ontology: a tool for the unification of genome annotations. Genome Biol. 6, R44 (2005).
- Dobin A. et al.. STAR: ultrafast universal RNA-seq aligner. Bioinformatics 29, 15–21 (2013).
- Hsu F. et al.. The UCSC known genes. Bioinformatics 22, 1036–1046 (2006).
- Lawrence M. et al.. Software for computing and annotating genomic ranges. PLoS Comput. Biol. 9, e1003118 (2013).
- Anders S. & Huber W. DESeq: differential expression analysis for sequence count data. Genome Biol. 11, R106 (2010).
- Love M. I., Huber W. & Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15, 550 (2014).
- Hochberg Y. A sharper bonferroni procedure for multiple tests of significance. Biometrika 75, 800–802 (1988).
- Li Y. I., Knowles D. A. & Pritchard J. K. LeafCutter: annotation-free quantification of RNA splicing. bioRxiv doi:; DOI: 10.1101/044107 (2016).
- Harrow J. et al.. GENCODE: the reference human genome annotation for The ENCODE Project. Genome Res. 22, 1760–1774 (2012).
- Van Haute L. et al.. Deficient methylation and formylation of mt-tRNA(Met) wobble cytosine in a patient carrying mutations in NSUN3. Nat. Commun. 7, 12039 (2016).
Source: PubMed