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.

Multiple acyl-CoA dehydrogenase deficiencies (MADDs) are a heterogeneous group of metabolic disorders with combined respiratory-chain deficiency and a neuromuscular phenotype. Despite recent advances in understanding the genetic basis of MADD, a number of cases remain unexplained. Here, we report clinically relevant variants in FLAD1, which encodes FAD synthase (FADS), as the cause of MADD and respiratory-chain dysfunction in nine individuals recruited from metabolic centers in six countries. In most individuals, we identified biallelic frameshift variants in the molybdopterin binding (MPTb) domain, located upstream of the FADS domain. Inasmuch as FADS is essential for cellular supply of FAD cofactors, the finding of biallelic frameshift variants was unexpected. Using RNA sequencing analysis combined with protein mass spectrometry, we discovered FLAD1 isoforms, which only encode the FADS domain. The existence of these isoforms might explain why affected individuals with biallelic FLAD1 frameshift variants still harbor substantial FADS activity. Another group of individuals with a milder phenotype responsive to riboflavin were shown to have single amino acid changes in the FADS domain. When produced in E. coli, these mutant FADS proteins resulted in impaired but detectable FADS activity; for one of the variant proteins, the addition of FAD significantly improved protein stability, arguing for a chaperone-like action similar to what has been reported in other riboflavin-responsive inborn errors of metabolism. In conclusion, our studies identify FLAD1 variants as a cause of potentially treatable inborn errors of metabolism manifesting with MADD and shed light on the mechanisms by which FADS ensures cellular FAD homeostasis.

The 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 the estimation of different bioenergetic measures by monitoring the oxygen consumption rates (OCR) of living cells in multi-well plates. However, studies of the statistical best practices for determining aggregated OCR measurements and comparisons have been lacking. Therefore, to understand how OCR behaves across different biological samples, wells, and plates, we performed mitochondrial stress tests in 126 96-well plates involving 203 fibroblast cell lines. 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 the 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-plate variations. This led to a significant reduction of the coefficient of variation across plates of basal respiration by 45% and of maximal respiration by 29%. Moreover, using positive and negative controls, we show that our statistical test outperforms the existing methods, which suffer from an excess of either false positives (within-plate methods), or false negatives (between-plate methods). Altogether, this study provides statistical good practices to support experimentalists in designing, analyzing, testing, and reporting the results of mitochondrial stress tests using this high throughput platform.