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.

44Across a large variety of Mendelian disorders, ~50-75% of patients do not receive a 45 genetic diagnosis by whole exome sequencing indicative of underlying disease-causing 46 variants in non-coding regions. In contrast, whole genome sequencing facilitates the 47 discovery of all genetic variants, but their sizeable number, coupled with a poor 48 understanding of the non-coding genome, makes their prioritization challenging. Here, we 49 demonstrate the power of transcriptome sequencing to provide a confirmed genetic 50 diagnosis for 10% (5 of 48) of undiagnosed mitochondrial disease patients and identify 51 strong candidate genes for patients remaining without diagnosis. We found a median of 1 52 aberrantly expressed gene, 5 aberrant splicing events, and 6 mono-allelically expressed 53 rare variants in patient-derived fibroblasts and established disease-causing roles for each 54 kind. Private exons often arose from sites that are weakly spliced in other individuals, 55providing an important clue for future variant prioritization. One such intronic exon-56 creating variant was found in three unrelated families in the complex I assembly factor 57 TIMMDC1, which we consequently established as a novel disease-associated gene. In 58 conclusion, our study expands the diagnostic tools for detecting non-exonic variants of 59Mendelian disorders and provides examples of intronic loss-of-function variants with 60 pathological relevance. 61Despite the revolutionizing impact of whole exome sequencing (WES) on the molecular 62 genetics of Mendelian disorders, ~50-75% of the patients do not receive a genetic diagnosis after 63 WES [1][2][3][4][5][6] . The disease-causing variants might be detected by WES but remain as variants of 64 unknown significance (VUS, Methods) or they are missed due to the inability to prioritize them. 65Many of these VUS are synonymous or non-coding variants that may affect RNA abundance or 66 isoform but cannot be prioritized due to the poor understanding of regulatory sequence to date 67 compared to coding sequence. Furthermore, WES covers only the 2% exonic regions of the 68 genome. Accordingly, it is mostly blind to regulatory variants in non-coding regions that could 69 affect RNA sequence and abundance. While the limitation of genome coverage is overcome by 70 whole genome sequencing (WGS), prioritization and interpretation of variants identified by 71 WGS is in turn limited by their amount [7][8][9] . 72With RNA sequencing (RNA-seq), limitations of the sole genetic information can be 73 complemented by directly probing variations in RNA abundance and in RNA sequence, 74 including allele-specific expression and splice isoforms. At least three extreme situations can be 75 directly interpreted to prioritize candidate disease-causing genes for a rare disorder. First, the 76 expression level of a gene can lie outside its physiological range. Genes with expression outside 77 their physical range can be identified as expression outliers, often using a stringent cutoff on 78 expression variat...

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.

Leber’s hereditary optic neuropathy (LHON) is the most frequent mitochondrial disease and was the first to be genetically defined by a point mutation in mitochondrial DNA (mtDNA). A molecular diagnosis is achieved in up to 95% of cases, the vast majority of which are accounted for by 3 mutations within mitochondrial complex I subunit–encoding genes in the mtDNA (mtLHON). Here, we resolve the enigma of LHON in the absence of pathogenic mtDNA mutations. We describe biallelic mutations in a nuclear encoded gene, DNAJC30 , in 33 unsolved patients from 29 families and establish an autosomal recessive mode of inheritance for LHON (arLHON), which to date has been a prime example of a maternally inherited disorder. Remarkably, all hallmarks of mtLHON were recapitulated, including incomplete penetrance, male predominance, and significant idebenone responsivity. Moreover, by tracking protein turnover in patient-derived cell lines and a DNAJC30 -knockout cellular model, we measured reduced turnover of specific complex I N-module subunits and a resultant impairment of complex I function. These results demonstrate that DNAJC30 is a chaperone protein needed for the efficient exchange of complex I subunits exposed to reactive oxygen species and integral to a mitochondrial complex I repair mechanism, thereby providing the first example to our knowledge of a disease resulting from impaired exchange of assembled respiratory chain subunits.

Overall survival of patients with BA in Switzerland compares favorably with current international standards, whereas results of the Kasai operation could be improved to reduce the need for LTs in infancy and early childhood.