Search

BackgroundThe spectrum of mitochondrial disease is genetically and phenotypically diverse, resulting from pathogenic variants in over 400 genes, with aerobic energy metabolism defects as a common denominator. Such heterogeneity poses a significant challenge in making an accurate diagnosis, critical for precision medicine.MethodsIn an international collaboration initiated by the European Network for Mitochondrial Diseases (GENOMIT) we recruited 2,023 pediatric patients at 11 specialist referral centers between October 2010 and January 2021, accumulating exome sequencing and HPO-encoded phenotype data. An exome-wide search for variants in known and potential novel disease genes, complemented by functional studies, followed ACMG guidelines.Results1,109 cases (55%) received a molecular diagnosis, of which one fifth have potential disease-modifying treatments (236/1,109, 21%). Functional studies enabled diagnostic uplift from 36% to 55% and discovery of 62 novel disease genes. Pathogenic variants were identified within genes encoding mitochondrial proteins or RNAs in 801 cases (72%), while, given extensive phenotype overlap, the remainder involved proteins targeted to other cellular compartments. To delineate genotype-phenotype associations, our data was complemented with registry and literature data to develop “GENOMITexplorer”, an open access resource detailing patient- (n=3,940), gene- (n=427), and variant-level (n=1,492) associations (prokischlab.github.io/GENOMITexplorer/).ConclusionsReaching a molecular diagnosis was essential for implementation of precision medicine and clinical trial eligibility, underlining the need for genome-wide screening given inability to accurately define mitochondrial diseases clinically. Key to diagnostic success were functional studies, encouraging early acquisition of patient- derived tissues and routine integration of high-throughput functional data to improve patient care by uplifting diagnostic rate.

RNA sequencing (RNA-seq) has emerged as a powerful approach to discover disease-causing gene regulatory defects in individuals affected by genetically undiagnosed rare disorders. Pioneering studies have shown that RNA-seq could increase the diagnosis rates over DNA sequencing alone by 8–36%, depending on the disease entity and tissue probed. To accelerate adoption of RNA-seq by human genetics centers, detailed analysis protocols are now needed. We present a step-by-step protocol that details how to robustly detect aberrant expression levels, aberrant splicing and mono-allelic expression in RNA-seq data using dedicated statistical methods. We describe how to generate and assess quality control plots and interpret the analysis results. The protocol is based on the detection of RNA outliers pipeline (DROP), a modular computational workflow that integrates all the analysis steps, can leverage parallel computing infrastructures and generates browsable web page reports. The authors describe a modular pipeline to detect aberrant gene expression events (expression level, splicing and mono-allelic expression) from patient RNA sequencing data, which can complement DNA-based diagnosis by enhancing the functional interpretation of variants.

The regulation of mitochondrial proteome is unique in that its components have origins in both mitochondria and nucleus. With the development of OMICS technologies, emerging evidence indicates an interaction between mitochondria and nucleus based not only on the proteins but also on the non-coding RNAs (ncRNAs). It is now accepted that large parts of the non‐coding genome are transcribed into various ncRNA species. Although their characterization has been a hot topic in recent years, the function of the majority remains unknown. Recently, ncRNA species microRNA (miRNA) and long-non coding RNAs (lncRNA) have been gaining attention as direct or indirect modulators of the mitochondrial proteome homeostasis. These ncRNA can impact mitochondria indirectly by affecting transcripts encoding for mitochondrial proteins in the cytoplasm. Furthermore, reports of mitochondria-localized miRNAs, termed mitomiRs, and lncRNAs directly regulating mitochondrial gene expression suggest the import of RNA to mitochondria, but also transcription from the mitochondrial genome. Interestingly, ncRNAs have been also shown to hide small open reading frames (sORFs) encoding for small functional peptides termed micropeptides, with several examples reported with a role in mitochondria. In this review, we provide a literature overview on ncRNAs and micropeptides found to be associated with mitochondrial biology in the context of both health and disease. Although reported, small study overlap and rare replications by other groups make the presence, transport, and role of ncRNA in mitochondria an attractive, but still challenging subject. Finally, we touch the topic of their potential as prognosis markers and therapeutic targets.

ADP-ribosylation is a reversible posttranslational modification used to regulate protein function. ADP-ribosyltransferases transfer ADP-ribose from NAD(+) to the target protein, and ADP-ribosylhydrolases, such as ADPRHL2, reverse the reaction. We used exome sequencing to identify five different bi-allelic pathogenic ADPRHL2 variants in 12 individuals from 8 families affected by a neurodegenerative disorder manifesting in childhood or adolescence with key clinical features including developmental delay or regression, seizures, ataxia, and axonal (sensori-)motor neuropathy. ADPRHL2 was virtually absent in available affected individuals’ fibroblasts, and cell viability was reduced upon hydrogen peroxide exposure, although it was rescued by expression of wild-type ADPRHL2 mRNA as well as treatment with a PARP1 inhibitor. Our findings suggest impaired protein ribosylation as another pathway that, if disturbed, causes neurodegenerative diseases.

Recently, heterozygous de novo mutations in SCL1A2 have been reported to underlie severe early-onset epileptic encephalopathy. In one male presenting with epileptic seizures and visual impairment, we identified a novel homozygous splicing variant in SCL1A2 (c.1421 + 1G > C) by using exome sequencing. Functional studies on cDNA level confirmed a consecutive loss of function. Our findings suggest that not only de novo mutations but also biallelic variants in SLC1A2 can cause epilepsy and that there is an additional autosomal recessive mode of inheritance. These findings also contribute to the understanding of the genetic mechanism of autosomal dominant SLC1A2-related epileptic encephalopathy as they exclude haploinsufficiency as exclusive genetic mechanism.

Leber's hereditary optic neuropathy (LHON) is an inherited mitochondrial disease that usually leads to acute or subacute bilateral central vision loss. In 95% of cases, LHON is caused by one of three primary mutations of the mitochondrial DNA (mtDNA), m.11778G>A in the MT-ND4 gene, m.14484T>C in the MT-ND6 gene, or m.3460G>A in the MT-ND1 gene. Here we characterize clinically, genetically, and biochemically a LHON family with multiple patients harboring two of these primary LHON mutations, m.11778G>A homoplasmic and m.14484T>C heteroplasmic. The unusually low male-to-female ratio of affected family members is also seen among the other patients previously reported with two primary LHON mutations m.11778G>A and m.14484T>C. While the index patient had very late onset of symptoms at 75years and severe visual loss, her two daughters had both onset in childhood (6 and 9years), with moderate to mild visual loss. A higher degree of heteroplasmy of the m.14484T>C mutation was found to correlate with an earlier age at onset in this family. Ours is the first LHON family harboring two primary LHON mutations where functional studies were performed in several affected family members. A more pronounced bioenergetic defect was found to correlate with an earlier age at onset. The patient with the earliest age at onset had a more significant complex I dysfunction than all controls, including the LHON patient with only the m.11778G>A mutation, suggesting a synergistic effect of the two primary LHON mutations in this patient.

Aberrant splicing is a major cause of rare diseases, yet its prediction from genome sequence remains in most cases inconclusive. Recently, RNA sequencing has proven to be an effective complementary avenue to detect aberrant splicing. Here, we developed FRASER, an algorithm to detect aberrant splicing from RNA sequencing data. Unlike existing methods, FRASER captures not only alternative splicing but also intron retention events. This typically doubles the number of detected aberrant events and identified a pathogenic intron retention in MCOLN1. FRASER automatically controls for latent confounders, which are widespread and substantially affect sensitivity. Moreover, FRASER is based on a count distribution and multiple testing correction, reducing the number of calls by two orders of magnitude over commonly applied z score cutoffs, with a minor sensitivity loss. The application to rare disease diagnostics is demonstrated by reprioritizing a pathogenic aberrant exon truncation in TAZ from a published dataset. FRASER is easy to use and freely available.

Inherited optic neuropathies include complex phenotypes, mostly driven by mitochondrial dysfunction. We report an optic atrophy spectrum disorder, including retinal macular dystrophy and kidney insufficiency leading to transplantation, associated with mitochondrial DNA (mtDNA) depletion without accumulation of multiple deletions. By whole-exome sequencing, we identified mutations affecting the mitochondrial single-strand binding protein (SSBP1) in 4 families with dominant and 1 with recessive inheritance. We show that SSBP1 mutations in patient-derived fibroblasts variably affect the amount of SSBP1 protein and alter multimer formation, but not the binding to ssDNA. SSBP1 mutations impaired mtDNA, nucleoids, and 7S-DNA amounts as well as mtDNA replication, affecting replisome machinery. The variable mtDNA depletion in cells was reflected in severity of mitochondrial dysfunction, including respiratory efficiency, OXPHOS subunits, and complex amount and assembly. mtDNA depletion and cytochrome c oxidase-negative cells were found ex vivo in biopsies of affected tissues, such as kidney and skeletal muscle. Reduced efficiency of mtDNA replication was also reproduced in vitro, confirming the pathogenic mechanism. Furthermore, ssbp1 suppression in zebrafish induced signs of nephropathy and reduced optic nerve size, the latter phenotype complemented by WT mRNA but not by SSBP1 mutant transcripts. This previously unrecognized disease of mtDNA maintenance implicates SSBP1 mutations as a cause of human pathology.

Accurate quantification of cellular and mitochondrial bioenergetic activity is of great interest in medicine and biology. Mitochondrial stress tests performed with Seahorse Bioscience XF Analyzers allow estimating different bioenergetic measures by monitoring oxygen consumption rates (OCR) of living cells in multi-well plates. However, studies of statistical best practices for determining OCR measurements and comparisons have been lacking so far. Therefore, we performed mitochondrial stress tests in 126 96-well plates involving 203 fibroblast cell lines to understand how OCR behaves across different biosamples, wells, and plates. We show that the noise of OCR is multiplicative, that outlier data points can concern individual measurements or all measurements of a well, and that the inter-plate variation is greater than intra-plate variation. Based on these insights, we developed a novel statistical method, OCR-Stats, that: i) robustly estimates OCR levels modeling multiplicative noise and automatically identifying outlier data points and outlier wells; and ii) performs statistical testing between samples, taking into account the different magnitudes of the between- and within-plates variations. This led to a significant reduction of the coefficient of variation across plates of basal respiration by 36% and of maximal respiration by 32%. Moreover, using positive and negative controls, we show that our statistical test outperforms existing methods, which either suffer from an excess of false positives (within-plates methods), or of false negatives (between-plates methods). Altogether, the aim of this study is to propose statistical good practices to support experimentalists in designing, analyzing, testing and reporting results of mitochondrial stress tests using this high throughput platform.

Despite their importance in determining protein abundance, a comprehensive catalogue of sequence features controlling protein-to-mRNA (PTR) ratios and a quantification of their effects is still lacking. Here we quantified PTR ratios for 11,575 proteins across 29 human tissues using matched transcriptomes and proteomes. We analyzed the contribution of known sequence determinants of protein synthesis and degradation and 15 novel mRNA and protein sequence motifs that we found by association testing. While the dynamic range of PTR ratios spans more than 2 orders of magnitude, our integrative model predicts PTR ratios at a median precision of 3.2-fold. A reporter assay provided significant functional support for two novel UTR motifs and a proteome-wide competition-binding assay identified motif-specific bound proteins for one motif. Moreover, our direct comparison of protein to RNA levels led to a new metrics of codon optimality. Altogether, this study shows that a large fraction of PTR ratio variance across genes can be predicted from sequence and identified many new candidate post-transcriptional regulatory elements in the human genome.