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.