Genetic diagnosis of Mendelian disorders via RNA sequencing

Clinical sequencing: From raw data to diagnosis with lifetime value

High-throughput sequencing (HTS) has revolutionized genetics by enabling the detection of sequence variants at hitherto unprecedented large scale. Despite these advances, however, there are still remaining challenges in the complete coverage of targeted regions (genes, exome or genome) as well as in HTS data analysis and interpretation. Moreover, it is easy to get overwhelmed by the plethora of available methods and tools for HTS. Here, we review the step-by-step process from the generation of sequence data to molecular diagnosis of Mendelian diseases. Highlighting advantages and limitations, this review addresses the current state of (1) HTS technologies, considering targeted, whole-exome, and whole-genome sequencing on short- and long-read platforms; (2) read alignment, variant calling and interpretation; as well as (3) regulatory issues related to genetic counseling, reimbursement, and data storage.

Rewiring of the Human Mitochondrial Interactome during Neuronal Reprogramming Reveals Regulators of the Respirasome and Neurogenesis

Mitochondrial protein (MP) assemblies undergo alterations during neurogenesis, a complex process vital in brain homeostasis and disease. Yet which MP assemblies remodel during differentiation remains unclear. Here, using mass spectrometry-based co-fractionation profiles and phosphoproteomics, we generated mitochondrial interaction maps of human pluripotent embryonal carcinoma stem cells and differentiated neuronal-like cells, which presented as two discrete cell populations by single-cell RNA sequencing. The resulting networks, encompassing 6,442 high-quality associations among 600 MPs, revealed widespread changes in mitochondrial interactions and site-specific phosphorylation during neuronal differentiation. By leveraging the networks, we show the orphan C20orf24 as a respirasome assembly factor whose disruption markedly reduces respiratory chain activity in patients deficient in complex IV. We also find that a heme-containing neurotrophic factor, neuron-derived neurotrophic factor [NENF], couples with Parkinson disease-related proteins to promote neurotrophic activity. Our results provide insights into the dynamic reorganization of mitochondrial networks during neuronal differentiation and highlights mechanisms for MPs in respirasome, neuronal function, and mitochondrial diseases.

•

Rewiring of mitochondrial (mt) protein interaction network in distinct cell states

•

Dramatic changes in site-specific phosphorylation during neuronal differentiation

•

C20orf24 is a respirasome assembly factor depleted in patients deficient in CIV

•

NENF binding with DJ-1/PINK1 promotes neurotrophic activity and neuronal survival

P5CS expression study in a new family with ALDH18A1-associated hereditary spastic paraplegia SPG9

In 2015–2016, we and others reported ALDH18A1 mutations causing dominant (SPG9A) or recessive (SPG9B) spastic paraplegia. In vitro production of the ALDH18A1 product, Δ¹‐pyrroline‐5‐carboxylate synthetase (P5CS), appeared necessary for cracking SPG9 disease‐causing mechanisms. We now describe a baculovirus–insect cell system that yields mgs of pure human P5CS and that has proven highly valuable with two novel P5CS mutations reported here in new SPG9B patients. We conclude that both mutations are disease‐causing, that SPG9B associates with partial P5CS deficiency and that it is clinically more severe than SPG9A, as reflected in onset age, disability, cognitive status, growth, and dysmorphic traits.

Genomic testing in myeloid malignancy

Clinical genetic testing in the myeloid malignancies is undergoing a rapid transition from the era of cytogenetics and single-gene testing to an era dominated by next-generation sequencing (NGS). This transition promises to better reveal the genetic alterations underlying disease, but there are distinct risks and benefits associated with different NGS testing platforms. NGS offers the potential benefit of being able to survey alterations across a wider set of genes, but analytic and clinical challenges associated with incidental findings, germ line variation, turnaround time, and limits of detection must be addressed. Additionally, transcriptome-based testing may offer several distinct benefits beyond traditional DNA-based methods. In addition to testing at disease diagnosis, research indicates potential benefits of genetic testing both prior to disease onset and at remission. In this review, we discuss the transition from the era of cytogenetics and single-gene tests to the era of NGS panels and genome-wide sequencing-highlighting both the potential and drawbacks of these novel technologies.

Genome sequencing and implications for rare disorders

The practice of genomic medicine stands to revolutionize our approach to medical care, and to realize this goal will require discovery of the relationship between rare variation at each of the ~ 20,000 protein-coding genes and their consequent impact on individual health and expression of Mendelian disease. The step-wise evolution of broad-based, genome-wide cytogenetic and molecular genomic testing approaches (karyotyping, chromosomal microarray [CMA], exome sequencing [ES]) has driven much of the rare disease discovery to this point, with genome sequencing representing the newest member of this team. Each step has brought increased sensitivity to interrogate individual genomic variation in an unbiased method that does not require clinical prediction of the locus or loci involved. Notably, each step has also brought unique limitations in variant detection, for example, the low sensitivity of ES for detection of triploidy, and of CMA for detection of copy neutral structural variants. The utility of genome sequencing (GS) as a clinical molecular diagnostic test, and the increased sensitivity afforded by addition of long-read sequencing or other -omics technologies such as RNAseq or metabolomics, are not yet fully explored, though recent work supports improved sensitivity of variant detection, at least in a subset of cases. The utility of GS will also rely upon further elucidation of the complexities of genetic and allelic heterogeneity, multilocus rare variation, and the impact of rare and common variation at a locus, as well as advances in functional annotation of identified variants. Much discovery remains to be done before the potential utility of GS is fully appreciated.

Molecular diagnostics of Mendelian disorders via combined DNA and RNA sequencing

The diagnostic yield in rare disorders is currently less than 50% although sequencing technologies in use are able to detect the majority of possible variants in our genome. The diagnostic gap is in part due to limitations in prioritizing and interpreting identified variants. The integration of functional data, such as transcriptomics, is emerging as a powerful complementary tool in diagnostics. It is able to quantify aberrant splicing, validate nonsense-mediated mRNA decay for potential loss-of-function variants, identify mono-allelically expressed variants, and help prioritize variants not predicted to change the encoded protein. Moreover, RNA-sequencing has been validated as a tool for the discovery of pathogenic variants in novel Mendelian disease genes. As RNA sequencing provides complementary information to DNA sequencing and can easily be established in addition to DNA sequencing, it has great potential for implementation as a routine tool for improving molecular diagnosis.

Identification of rare-disease genes using blood transcriptome sequencing and large control cohorts

It is estimated that 350 million individuals worldwide suffer from rare diseases, which are predominantly caused by mutation in a single gene1. The current molecular diagnostic rate is estimated at 50%, with whole-exome sequencing (WES) among the most successful approaches2-5. For patients in whom WES is uninformative, RNA sequencing (RNA-seq) has shown diagnostic utility in specific tissues and diseases6-8. This includes muscle biopsies from patients with undiagnosed rare muscle disorders6,9, and cultured fibroblasts from patients with mitochondrial disorders7. However, for many individuals, biopsies are not performed for clinical care, and tissues are difficult to access. We sought to assess the utility of RNA-seq from blood as a diagnostic tool for rare diseases of different pathophysiologies. We generated whole-blood RNA-seq from 94 individuals with undiagnosed rare diseases spanning 16 diverse disease categories. We developed a robust approach to compare data from these individuals with large sets of RNA-seq data for controls (n = 1,594 unrelated controls and n = 49 family members) and demonstrated the impacts of expression, splicing, gene and variant filtering strategies on disease gene identification. Across our cohort, we observed that RNA-seq yields a 7.5% diagnostic rate, and an additional 16.7% with improved candidate gene resolution.

Identification and Validation of Novel Reference Genes in Acute Lymphoblastic Leukemia for Droplet Digital PCR

Droplet digital PCR is the most robust method for absolute nucleic acid quantification. However, RNA is a very versatile molecule and its abundance is tissue-dependent. RNA quantification is dependent on a reference control to estimate the abundance. Additionally, in cancer, many cellular processes are deregulated which consequently affects the gene expression profiles. In this work, we performed microarray data mining of different childhood cancers and healthy controls. We selected four genes that showed no gene expression variations (PSMB6, PGGT1B, UBQLN2 and UQCR2) and four classical reference genes (ACTB, GAPDH, RPL4 and RPS18). Gene expression was validated in 40 acute lymphoblastic leukemia samples by means of droplet digital PCR. We observed that PSMB6, PGGT1B, UBQLN2 and UQCR2 were expressed ~100 times less than ACTB, GAPDH, RPL4 and RPS18. However, we observed excellent correlations among the new reference genes (p < 0.0001). We propose that PSMB6, PGGT1B, UBQLN2 and UQCR2 are housekeeping genes with low expression in childhood cancer.

Long-Read Sequencing Emerging in Medical Genetics

The wide implementation of next-generation sequencing (NGS) technologies has revolutionized the field of medical genetics. However, the short read lengths of currently used sequencing approaches pose a limitation for the identification of structural variants, sequencing repetitive regions, phasing of alleles and distinguishing highly homologous genomic regions. These limitations may significantly contribute to the diagnostic gap in patients with genetic disorders who have undergone standard NGS, like whole exome or even genome sequencing. Now, the emerging long-read sequencing (LRS) technologies may offer improvements in the characterization of genetic variation and regions that are difficult to assess with the prevailing NGS approaches. LRS has so far mainly been used to investigate genetic disorders with previously known or strongly suspected disease loci. While these targeted approaches already show the potential of LRS, it remains to be seen whether LRS technologies can soon enable true whole genome sequencing routinely. Ultimately, this could allow the de novo assembly of individual whole genomes used as a generic test for genetic disorders. In this article, we summarize the current LRS-based research on human genetic disorders and discuss the potential of these technologies to facilitate the next major advancements in medical genetics.

Biallelic mutations in PIGP cause developmental and epileptic encephalopathy

Developmental and epileptic encephalopathies are characterized by infantile seizures and psychomotor delay. Glycosylphosphatidylinositol biosynthesis defects, resulting in impaired tethering of various proteins to the cell surface, represent the underlying pathology in some patients. One of the genes involved, PIGP, has recently been associated with infantile seizures and developmental delay in two siblings. Here, we report the second family with a markedly overlapping phenotype due to a homozygous frameshift mutation (c.456delA;p.Glu153Asnfs*34) in PIGP. Flow cytometry of patient granulocytes confirmed reduced expression of glycosylphosphatidylinositol-anchored proteins as functional consequence. Our findings corroborate PIGP as a monogenic disease gene for developmental and epileptic encephalopathy.

A patient with homozygous nonsense variants in two Leigh syndrome disease genes: Distinguishing a dual diagnosis from a hypomorphic protein-truncating variant

Leigh syndrome is a mitochondrial disease caused by pathogenic variants in over 85 genes. Whole exome sequencing of a patient with Leigh-like syndrome identified homozygous protein-truncating variants in two genes associated with Leigh syndrome; a reported pathogenic variant in PDHX (NP_003468.2:p.(Arg446*)), and an uncharacterized variant in complex I (CI) assembly factor TIMMDC1 (NP_057673.2:p.(Arg225*)). The TIMMDC1 variant was predicted to truncate 61 amino acids at the C-terminus and functional studies demonstrated a hypomorphic impact of the variant on CI assembly. However, the mutant protein could still rescue CI assembly in TIMMDC1 knockout cells and the patient's clinical phenotype was not clearly distinct from that of other patients with the same PDHX defect. Our data suggest that the hypomorphic effect of the TIMMDC1 protein-truncating variant does not constitute a dual diagnosis in this individual. We recommend cautious assessment of variants in the C-terminus of TIMMDC1 and emphasize the need to consider the caveats detailed within the American College of Medical Genetics and Genomics (ACMG) criteria when assessing variants.

A toolkit for genetics providers in follow-up of patients with non-diagnostic exome sequencing

There are approximately 7,000 rare diseases affecting 25-30 million Americans, with 80% estimated to have a genetic basis. This presents a challenge for genetics practitioners to determine appropriate testing, make accurate diagnoses, and conduct up-to-date patient management. Exome sequencing (ES) is a comprehensive diagnostic approach, but only 25%-41% of the patients receive a molecular diagnosis. The remaining three-fifths to three-quarters of patients undergoing ES remain undiagnosed. The Stanford Center for Undiagnosed Diseases (CUD), a clinical site of the Undiagnosed Diseases Network, evaluates patients with undiagnosed and rare diseases using a combination of methods including ES. Frequently these patients have non-diagnostic ES results, but strategic follow-up techniques identify diagnoses in a subset. We present techniques used at the CUD that can be adopted by genetics providers in clinical follow-up of cases where ES is non-diagnostic. Solved case examples illustrate different types of non-diagnostic results and the additional techniques that led to a diagnosis. Frequent approaches include segregation analysis, data reanalysis, genome sequencing, additional variant identification, careful phenotype-disease correlation, confirmatory testing, and case matching. We also discuss prioritization of cases for additional analyses.

Insights into genetics, human biology and disease gleaned from family based genomic studies

Identifying genes and variants contributing to rare disease phenotypes and Mendelian conditions informs biology and medicine, yet potential phenotypic consequences for variation of >75% of the ~20,000 annotated genes in the human genome are lacking. Technical advances to assess rare variation genome-wide, particularly exome sequencing (ES), enabled establishment in the United States of the National Institutes of Health (NIH)-supported Centers for Mendelian Genomics (CMGs) and have facilitated collaborative studies resulting in novel "disease gene" discoveries. Pedigree-based genomic studies and rare variant analyses in families with suspected Mendelian conditions have led to the elucidation of hundreds of novel disease genes and highlighted the impact of de novo mutational events, somatic variation underlying nononcologic traits, incompletely penetrant alleles, phenotypes with high locus heterogeneity, and multilocus pathogenic variation. Herein, we highlight CMG collaborative discoveries that have contributed to understanding both rare and common diseases and discuss opportunities for future discovery in single-locus Mendelian disorder genomics. Phenotypic annotation of all human genes; development of bioinformatic tools and analytic methods; exploration of non-Mendelian modes of inheritance including reduced penetrance, multilocus variation, and oligogenic inheritance; construction of allelic series at a locus; enhanced data sharing worldwide; and integration with clinical genomics are explored. Realizing the full contribution of rare disease research to functional annotation of the human genome, and further illuminating human biology and health, will lay the foundation for the Precision Medicine Initiative.

Proteome Analysis of Human Neutrophil Granulocytes From Patients With Monogenic Disease Using Data-independent Acquisition

Neutrophil granulocytes are critical mediators of innate immunity and tissue regeneration. Rare diseases of neutrophil granulocytes may affect their differentiation and/or functions. However, there are very few validated diagnostic tests assessing the functions of neutrophil granulocytes in these diseases. Here, we set out to probe omics analysis as a novel diagnostic platform for patients with defective differentiation and function of neutrophil granulocytes. We analyzed highly purified neutrophil granulocytes from 68 healthy individuals and 16 patients with rare monogenic diseases. Cells were isolated from fresh venous blood (purity >99%) and used to create a spectral library covering almost 8000 proteins using strong cation exchange fractionation. Patient neutrophil samples were then analyzed by data-independent acquisition proteomics, quantifying 4154 proteins in each sample. Neutrophils with mutations in the neutrophil elastase gene ELANE showed large proteome changes that suggest these mutations may affect maturation of neutrophil granulocytes and initiate misfolded protein response and cellular stress mechanisms. In contrast, only few proteins changed in patients with leukocyte adhesion deficiency (LAD) and chronic granulomatous disease (CGD). Strikingly, neutrophil transcriptome analysis showed no correlation with its proteome. In case of two patients with undetermined genetic causes, proteome analysis guided the targeted genetic diagnostics and uncovered the underlying genomic mutations. Data-independent acquisition proteomics may help to define novel pathomechanisms in neutrophil diseases and provide a clinically useful diagnostic dimension.

Genomic Analysis in the Age of Human Genome Sequencing

Affordable genome sequencing technologies promise to revolutionize the field of human genetics by enabling comprehensive studies that interrogate all classes of genome variation, genome-wide, across the entire allele frequency spectrum. Ongoing projects worldwide are sequencing many thousands-and soon millions-of human genomes as part of various gene mapping studies, biobanking efforts, and clinical programs. However, while genome sequencing data production has become routine, genome analysis and interpretation remain challenging endeavors with many limitations and caveats. Here, we review the current state of technologies for genetic variant discovery, genotyping, and functional interpretation and discuss the prospects for future advances. We focus on germline variants discovered by whole-genome sequencing, genome-wide functional genomic approaches for predicting and measuring variant functional effects, and implications for studies of common and rare human disease.

National Rare Diseases Registry System of China and Related Cohort Studies: Vision and Roadmap

Rare diseases are major challenges in healthcare and medical research and are the basis of national development strategies in many countries. However, inadequate definition of rare diseases and lags in orphan drug development in China hinder rare disease research. In response, the first National Rare Diseases Registry System of China (NRDRS) was established, and various cohort studies have been launched since 2016. More than 20 top academic institutions in China are currently participating in this joint effort to carry out nationwide registration of rare diseases. The primary objectives are to establish standardization for the registration platform, build biobanks of genomic data, and create partnerships for data sharing and research collaboration. Innovative informatics technologies have been implemented to develop the NRDRS, including employment of ontological and knowledge bases to render standardization and support standard of care. Development of informatics analysis tools will facilitate accurate and more efficient diagnoses for rare diseases. Long-term research collaboration is encouraged to create additional national rare disease networks for research translation and to benefit patients with rare diseases. The NRDRS of China and related cohort studies are anticipated to enlighten rare disease research significantly in China.

Expanding the Boundaries of RNA Sequencing as a Diagnostic Tool for Rare Mendelian Disease

Gene-panel and whole-exome analyses are now standard methodologies for mutation detection in Mendelian disease. However, the diagnostic yield achieved is at best 50%, leaving the genetic basis for disease unsolved in many individuals. New approaches are thus needed to narrow the diagnostic gap. Whole-genome sequencing is one potential strategy, but it currently has variant-interpretation challenges, particularly for non-coding changes. In this study we focus on transcriptome analysis, specifically total RNA sequencing (RNA-seq), by using monogenetic neuromuscular disorders as proof of principle. We examined a cohort of 25 exome and/or panel "negative" cases and provided genetic resolution in 36% (9/25). Causative mutations were identified in coding and non-coding exons, as well as in intronic regions, and the mutational pathomechanisms included transcriptional repression, exon skipping, and intron inclusion. We address a key barrier of transcriptome-based diagnostics: the need for source material with disease-representative expression patterns. We establish that blood-based RNA-seq is not adequate for neuromuscular diagnostics, whereas myotubes generated by transdifferentiation from an individual's fibroblasts accurately reflect the muscle transcriptome and faithfully reveal disease-causing mutations. Our work confirms that RNA-seq can greatly improve diagnostic yield in genetically unresolved cases of Mendelian disease, defines strengths and challenges of the technology, and demonstrates the suitability of cell models for RNA-based diagnostics. Our data set the stage for development of RNA-seq as a powerful clinical diagnostic tool that can be applied to the large population of individuals with undiagnosed, rare diseases and provide a framework for establishing minimally invasive strategies for doing so.

Translational Metabolism: A multidisciplinary approach towards precision diagnosis of inborn errors of metabolism in the omics era

The laboratory diagnosis of inborn errors of metabolism has been revolutionized in recent years, thanks to the amazing developments in the field of DNA sequencing including whole exome and whole genome sequencing (WES and WGS). Interpretation of the results coming from WES and/or WGS analysis is definitely not trivial especially since the biological relevance of many of the variants identified by WES and/or WGS, have not been tested experimentally and prediction programs like POLYPHEN-2 and SIFT are far from perfect. Correct interpretation of WES and/or WGS results can only be achieved by performing functional studies at multiple levels (different metabolomics platforms, enzymology, in vitro and in vivo flux analysis), often requires studies in model organisms like zebra fish, Caenorhabditis elegans, Saccharomyces cerevisiae, mutant mice and others, and also requires the input of many different disciplines to make this Translational Metabolism approach effective.

Quantification and discovery of sequence determinants of protein-per-mRNA amount in 29 human tissues

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 are still lacking. Here, we quantified PTR ratios for 11,575 proteins across 29 human tissues using matched transcriptomes and proteomes. We estimated by regression the contribution of known sequence determinants of protein synthesis and degradation in addition to 45 mRNA and 3 protein sequence motifs that we found by association testing. While PTR ratios span more than 2 orders of magnitude, our integrative model predicts PTR ratios at a median precision of 3.2-fold. A reporter assay provided functional support for two novel UTR motifs, and an immobilized mRNA affinity competition-binding assay identified motif-specific bound proteins for one motif. Moreover, our integrative model led to a new metric of codon optimality that captures the effects of codon frequency on protein synthesis and degradation. Altogether, this study shows that a large fraction of PTR ratio variation in human tissues can be predicted from sequence, and it identifies many new candidate post-transcriptional regulatory elements.

Homozygous frame shift variant in ATP7B exon 1 leads to bypass of nonsense-mediated mRNA decay and to a protein capable of copper export

Wilson disease (WD) is an autosomal recessive disease of copper excess due to pathogenic variants in the ATP7B gene coding for a copper-transporting ATPase. We present a 5-year-old girl with the homozygous frame shift variant NM_000053.3: c.19_20del in exon 1 of ATP7B (consecutive exon numbering with c.1 as first nucleotide of exon 1), detected by whole-exome sequencing as a secondary finding. The variant leads to a premature termination codon in exon 2. The girl exhibited no WD symptoms and no abnormalities in liver biopsy. ATP7B liver mRNA expression was comparable to healthy controls suggesting that nonsense-mediated mRNA decay (NMD) could be bypassed by the mechanism of translation reinitiation. To verify this hypothesis, a CMV-driven ATP7B minigene (pcDNA3) was equipped with the authentic ATP7B 5' untranslated region  and a truncated intron 2. We introduced c.19_20del by site-directed mutagenesis and overexpressed the constructs in HEK293T cells. We analyzed ATP7B expression by qRT-PCR, northern and western blot, and examined protein function by copper export capacity assays. Northern blot, qRT-PCR, and western blot revealed that c.19_20del ATP7B mRNA and protein is expressed in size and amount comparable to wild-type. Copper export capacity was also comparable to wild-type. Our results indicate that c.19_20del in ATP7B is able to bypass NMD by translation reinitiation, demonstrating that the classification of truncating variants as pathogenic without additional investigations should be done carefully.

Spatially clustering de novo variants in CYFIP2, encoding the cytoplasmic FMRP interacting protein 2, cause intellectual disability and seizures

CYFIP2, encoding the evolutionary highly conserved cytoplasmic FMRP interacting protein 2, has previously been proposed as a candidate gene for intellectual disability and autism because of its important role linking FMRP-dependent transcription regulation and actin polymerization via the WAVE regulatory complex (WRC). Recently, de novo variants affecting the amino acid p.Arg87 of CYFIP2 were reported in four individuals with epileptic encephalopathy. We here report 12 independent patients harboring a variety of de novo variants in CYFIP2 broadening the molecular and clinical spectrum of a novel CYFIP2-related neurodevelopmental disorder. Using trio whole-exome or -genome sequencing, we identified 12 independent patients carrying a total of eight distinct de novo variants in CYFIP2 with a shared phenotype of intellectual disability, seizures, and muscular hypotonia. We detected seven different missense variants, of which two occurred recurrently (p.(Arg87Cys) and p.(Ile664Met)), and a splice donor variant in the last intron for which we showed exon skipping in the transcript. The latter is expected to escape nonsense-mediated mRNA decay resulting in a truncated protein. Despite the large spacing in the primary structure, the variants spatially cluster in the tertiary structure and are all predicted to weaken the interaction with WAVE1 or NCKAP1 of the actin polymerization regulating WRC-complex. Preliminary genotype-phenotype correlation indicates a profound phenotype in p.Arg87 substitutions and a more variable phenotype in other alterations. This study evidenced a variety of de novo variants in CYFIP2 as a novel cause of mostly severe intellectual disability with seizures and muscular hypotonia.

Autosomal recessive congenital ichthyosis: Genomic landscape and phenotypic spectrum in a cohort of 125 consanguineous families

Autosomal recessive congenital ichthyosis (ARCI), a phenotypically heterogeneous group of non-syndromic Mendelian disorders of keratinization, is caused by mutations in as many as 13 distinct genes. We examined a cohort of 125 consanguineous families with ARCI for underlying genetic mutations. The patients' DNA was analyzed with a gene-targeted next generation sequencing panel comprising 38 ichthyosis associated genes. The interpretations of results of genomic data were assisted by genome-wide homozygosity mapping and transcriptome sequencing. Sequence data analysis identified biallelic mutations in 106 families out of a total of 125 (85%), most of them (102, 96.2%) being homozygous, reflecting consanguinity in these families. Among the 85 distinct mutations in 10 different genes, 45 (53%) were previously unreported. Phenotype-genotype correlations allowed assignment of specific genes in the majority of the families to a specific subtype of ARCI, lamellar ichthyosis (LI) versus congenital ichthyosiform erythroderma (CIE). Interestingly, mutations in several genes could give rise to an overlapping phenotype consistent with either LI or CIE. Also, this is the third report for SDR9C7 and SULT2B1, and fourth report for CERS3 mutations. Direct comparison of our results with previously published regional cohorts highlights the global mutation landscape of ARCI, however, population specific differences were noted.

Solving Mendelian Mysteries: The Non-coding Genome May Hold the Key

Despite revolutionary advances in sequencing approaches, many mendelian disorders have remained unexplained. In this issue of Cell, Aneichyk et al. combine genomic and cell-type-specific transcriptomic data to causally link a non-coding mutation in the ubiquitous TAF1 gene to X-linked dystonia-parkinsonism.

A Multiplexed Assay for Exon Recognition Reveals that an Unappreciated Fraction of Rare Genetic Variants Cause Large-Effect Splicing Disruptions

Mutations that lead to splicing defects can have severe consequences on gene function and cause disease. Here, we explore how human genetic variation affects exon recognition by developing a multiplexed functional assay of splicing using Sort-seq (MFASS). We assayed 27,733 variants in the Exome Aggregation Consortium (ExAC) within or adjacent to 2,198 human exons in the MFASS minigene reporter and found that 3.8% (1,050) of variants, most of which are extremely rare, led to large-effect splice-disrupting variants (SDVs). Importantly, we find that 83% of SDVs are located outside of canonical splice sites, are distributed evenly across distinct exonic and intronic regions, and are difficult to predict a priori. Our results indicate extant, rare genetic variants can have large functional effects on splicing at appreciable rates, even outside the context of disease, and MFASS enables their empirical assessment at scale.

Early-onset inflammatory bowel disease as a model disease to identify key regulators of immune homeostasis mechanisms

Rare, monogenetic diseases present unique models to dissect gene functions and biological pathways, concomitantly enhancing our understanding of the etiology of complex (and often more common) traits. Although inflammatory bowel disease (IBD) is a generally prototypic complex disease, it can also manifest in an early-onset, monogenic fashion, often following Mendelian modes of inheritance. Recent advances in genomic technologies have spurred the identification of genetic defects underlying rare, very early-onset IBD (VEO-IBD) as a disease subgroup driven by strong genetic influence, pinpointing key players in the delicate homeostasis of the immune system in the gut and illustrating the intimate relationships between bowel inflammation, systemic immune dysregulation, and primary immunodeficiency with increased susceptibility to infections. As for other human diseases, it is likely that adult-onset diseases may represent complex diseases integrating the effects of host genetic susceptibility and environmental triggers. Comparison of adult-onset IBD and VEO-IBD thus provides beautiful models to investigate the relationship between monogenic and multifactorial/polygenic diseases. This review discusses the present and novel findings regarding monogenic IBD as well as key questions and future directions of IBD research.

Mitochondrial Neurogastrointestinal Encephalomyopathy: Into the Fourth Decade, What We Have Learned So Far

Mitochondrial neurogastrointestinal encephalomyopathy (MNGIE) is an ultra-rare metabolic autosomal recessive disease, caused by mutations in the nuclear gene TYMP which encodes the enzyme thymidine phosphorylase. The resulting enzyme deficiency leads to a systemic accumulation of the deoxyribonucleosides thymidine and deoxyuridine, and ultimately mitochondrial failure due to a progressive acquisition of secondary mitochondrial DNA (mtDNA) mutations and mtDNA depletion. Clinically, MNGIE is characterized by gastrointestinal and neurological manifestations, including cachexia, gastrointestinal dysmotility, peripheral neuropathy, leukoencephalopathy, ophthalmoplegia and ptosis. The disease is progressively degenerative and leads to death at an average age of 37.6 years. As with the vast majority of rare diseases, patients with MNGIE face a number of unmet needs related to diagnostic delays, a lack of approved therapies, and non-specific clinical management. We provide here a comprehensive collation of the available knowledge of MNGIE since the disease was first described 42 years ago. This review includes symptomatology, diagnostic procedures and hurdles, in vitro and in vivo disease models that have enhanced our understanding of the disease pathology, and finally experimental therapeutic approaches under development. The ultimate aim of this review is to increase clinical awareness of MNGIE, thereby reducing diagnostic delay and improving patient access to putative treatments under investigation.

Diagnosing rare diseases after the exome

High-throughput sequencing has ushered in a diversity of approaches for identifying genetic variants and understanding genome structure and function. When applied to individuals with rare genetic diseases, these approaches have greatly accelerated gene discovery and patient diagnosis. Over the past decade, exome sequencing has emerged as a comprehensive and cost-effective approach to identify pathogenic variants in the protein-coding regions of the genome. However, for individuals in whom exome-sequencing fails to identify a pathogenic variant, we discuss recent advances that are helping to reduce the diagnostic gap.

Pathogenic copy number variants that affect gene expression contribute to genomic burden in cerebral palsy

Cerebral palsy (CP) is the most frequent movement disorder of childhood affecting 1 in 500 live births in developed countries. We previously identified likely pathogenic de novo or inherited single nucleotide variants (SNV) in 14% (14/98) of trios by exome sequencing and a further 5% (9/182) from evidence of outlier gene expression using RNA sequencing. Here, we detected copy number variants (CNV) from exomes of 186 unrelated individuals with CP (including our original 98 trios) using the CoNIFER algorithm. CNV were validated with Illumina 850 K SNP arrays and compared with RNA-Seq outlier gene expression analysis from lymphoblastoid cell lines (LCL). Gene expression was highly correlated with gene dosage effect. We resolved an additional 3.7% (7/186) of this cohort with pathogenic or likely pathogenic CNV while a further 7.7% (14/186) had CNV of uncertain significance. We identified recurrent genomic rearrangements previously associated with CP due to 2p25.3 deletion, 22q11.2 deletions and duplications and Xp monosomy. We also discovered a deletion of a single gene, PDCD6IP, and performed additional zebrafish model studies to support its single allele loss in CP aetiology. Combined SNV and CNV analysis revealed pathogenic and likely pathogenic variants in 22.7% of unselected individuals with CP.

OUTRIDER: A Statistical Method for Detecting Aberrantly Expressed Genes in RNA Sequencing Data

RNA sequencing (RNA-seq) is gaining popularity as a complementary assay to genome sequencing for precisely identifying the molecular causes of rare disorders. A powerful approach is to identify aberrant gene expression levels as potential pathogenic events. However, existing methods for detecting aberrant read counts in RNA-seq data either lack assessments of statistical significance, so that establishing cutoffs is arbitrary, or rely on subjective manual corrections for confounders. Here, we describe OUTRIDER (Outlier in RNA-Seq Finder), an algorithm developed to address these issues. The algorithm uses an autoencoder to model read-count expectations according to the gene covariation resulting from technical, environmental, or common genetic variations. Given these expectations, the RNA-seq read counts are assumed to follow a negative binomial distribution with a gene-specific dispersion. Outliers are then identified as read counts that significantly deviate from this distribution. The model is automatically fitted to achieve the best recall of artificially corrupted data. Precision-recall analyses using simulated outlier read counts demonstrated the importance of controlling for covariation and significance-based thresholds. OUTRIDER is open source and includes functions for filtering out genes not expressed in a dataset, for identifying outlier samples with too many aberrantly expressed genes, and for detecting aberrant gene expression on the basis of false-discovery-rate-adjusted p values. Overall, OUTRIDER provides an end-to-end solution for identifying aberrantly expressed genes and is suitable for use by rare-disease diagnostic platforms.

Expanded Insights Into Mechanisms of Gene Expression and Disease Related Disruptions

Definitive molecular diagnoses in disorders apparently due to genetic or genomic defects are still lacking in a significant number of investigated cases, despite use of studies designed to discover defects in the protein coding regions of the genome. Increasingly studies are being designed to search for defects in the non-protein coding genome, and for alterations in gene expression. Here we review new insights into genomic elements involved in control of gene expression, including methods to analyze chromatin that is accessible for transcription factor binding, enhancers, chromatin looping, transcription, RNA binding proteins, and alternative splicing. We review new studies on levels of genome organization, including the occurrence of transcriptional domains and their boundary elements. Information is presented on specific malformation syndromes that arise due to structural genomic changes that impact the non-protein coding genome and sometimes impact specific transcriptional domains. We also review convergence of genome-wide association with studies of gene expression, discoveries related to expression quantitative trait loci and splicing quantitative trait loci and the relevance of these to specific complex common diseases. Aspects of epigenetic mechanisms and clinical applications of analyses of methylation signatures are also discussed.

Perspectives on the Genomics of HSP Beyond Mendelian Inheritance

Hereditary Spastic Paraplegia is an extraordinarily heterogeneous disease caused by over 50 Mendelian genes. Recent applications of next-generation sequencing, large scale data analysis, and data sharing/matchmaking, have discovered a quickly expanding set of additional HSP genes. Since most recently discovered HSP genes are rare causes of the disease, there is a growing concern of a persisting diagnostic gap, estimated at 30-40%, and even higher for sporadic cases. This missing heritability may not be fully closed by classic Mendelian mutations in protein coding genes. Here we show strategies and published examples of broadening areas of attention for Mendelian and non-Mendelian causes of HSP. We suggest a more inclusive perspective on the potential final architecture of HSP genomics. Efforts to narrow the heritability gap will ultimately lead to more precise and comprehensive genetic diagnoses, which is the starting point for emerging, highly specific gene therapies.

RNA sequencing solved the most common but unrecognized NEB pathogenic variant in Japanese nemaline myopathy

PURPOSE

The diagnostic rate for Mendelian diseases by exome sequencing (ES) is typically 20-40%. The low rate is partly because ES misses deep-intronic or synonymous variants leading to aberrant splicing. In this study, we aimed to apply RNA sequencing (RNA-seq) to efficiently detect the aberrant splicings and their related variants.

METHODS

Aberrant splicing in biopsied muscles from six nemaline myopathy (NM) cases unresolved by ES were analyzed with RNA-seq. Variants related to detected aberrant splicing events were analyzed with Sanger sequencing. Detected variants were screened in NM patients unresolved by ES.

RESULTS

We identified a novel deep-intronic NEB pathogenic variant, c.1569+339A>G in one case, and another novel synonymous NEB pathogenic variant, c.24684G>C (p.Ser8228Ser) in three cases. The c.24684G>C variant was observed to be the most frequent among all NEB pathogenic variants in normal Japanese populations with a frequency of 1 in 178 (20 alleles in 3552 individuals), but was previously unrecognized. Expanded screening of the variant identified it in a further four previously unsolved nemaline myopathy cases.

CONCLUSION

These results indicated that RNA-seq may be able to solve a large proportion of previously undiagnosed muscle diseases.

Somatic alterations compromised molecular diagnosis of DOCK8 hyper-IgE syndrome caused by a novel intronic splice site mutation

In hyper-IgE syndromes (HIES), a group of primary immunodeficiencies clinically overlapping with atopic dermatitis, early diagnosis is crucial to initiate appropriate therapy and prevent irreversible complications. Identification of underlying gene defects such as in DOCK8 and STAT3 and corresponding molecular testing has improved diagnosis. Yet, in a child and her newborn sibling with HIES phenotype molecular diagnosis was misleading. Extensive analyses driven by the clinical phenotype identified an intronic homozygous DOCK8 variant c.4626 + 76 A > G creating a novel splice site as disease-causing. While the affected newborn carrying the homozygous variant had no expression of DOCK8 protein, in the index patient molecular diagnosis was compromised due to expression of altered and wildtype DOCK8 transcripts and DOCK8 protein as well as defective STAT3 signaling. Sanger sequencing of lymphocyte subsets revealed that somatic alterations and reversions revoked the predominance of the novel over the canonical splice site in the index patient explaining DOCK8 protein expression, whereas defective STAT3 responses in the index patient were explained by a T cell phenotype skewed towards central and effector memory T cells. Hence, somatic alterations and skewed immune cell phenotypes due to selective pressure may compromise molecular diagnosis and need to be considered with unexpected clinical and molecular findings.

CHD3 helicase domain mutations cause a neurodevelopmental syndrome with macrocephaly and impaired speech and language

Chromatin remodeling is of crucial importance during brain development. Pathogenic alterations of several chromatin remodeling ATPases have been implicated in neurodevelopmental disorders. We describe an index case with a de novo missense mutation in CHD3, identified during whole genome sequencing of a cohort of children with rare speech disorders. To gain a comprehensive view of features associated with disruption of this gene, we use a genotype-driven approach, collecting and characterizing 35 individuals with de novo CHD3 mutations and overlapping phenotypes. Most mutations cluster within the ATPase/helicase domain of the encoded protein. Modeling their impact on the three-dimensional structure demonstrates disturbance of critical binding and interaction motifs. Experimental assays with six of the identified mutations show that a subset directly affects ATPase activity, and all but one yield alterations in chromatin remodeling. We implicate de novo CHD3 mutations in a syndrome characterized by intellectual disability, macrocephaly, and impaired speech and language.

Evidence of mild founder LMOD3 mutations causing nemaline myopathy 10 in Germany and Austria

OBJECTIVE

To expand the clinical and genetic spectrum of nemaline myopathy 10 by a series of Austrian and German patients with a milder disease course and missense mutations in LMOD3.

METHODS

We characterized the clinical features and the genetic status of 4 unrelated adolescent or adult patients with nemaline myopathy.

RESULTS

The 4 patients showed a relatively mild disease course. They all have survived into adulthood, 3 of 4 have remained ambulatory, and all showed marked facial weakness. Muscle biopsy specimens gave evidence of nemaline bodies. All patients were unrelated but originated from Austria (Tyrol and Upper Austria) and Southern Germany (Bavaria). All patients carried the missense variant c.1648C>T, p.(Leu550Phe) in the LMOD3 gene, either on both alleles or in trans with another missense variant (c.1004A>G, p.Gln335Arg). Both variants were not reported previously.

CONCLUSIONS

In 2014, a severe form of congenital nemaline myopathy caused by disrupting mutations in LMOD3 was identified and denoted as NEM10. Unlike the previously reported patients, who had a severe clinical picture with a substantial risk of early death, our patients showed a relatively mild disease course. As the missense variant c.1648C>T is located further downstream compared to all previously published LMOD3 mutations, it might be associated with higher protein expression compared to the reported loss-of-function mutations. The apparent clusters of 2 mild mutations in Germany and Austria in 4 unrelated families may be explained by a founder effect.

Bi-allelic Mutations in NDUFA6 Establish Its Role in Early-Onset Isolated Mitochondrial Complex I Deficiency

Isolated complex I deficiency is a common biochemical phenotype observed in pediatric mitochondrial disease and often arises as a consequence of pathogenic variants affecting one of the ∼65 genes encoding the complex I structural subunits or assembly factors. Such genetic heterogeneity means that application of next-generation sequencing technologies to undiagnosed cohorts has been a catalyst for genetic diagnosis and gene-disease associations. We describe the clinical and molecular genetic investigations of four unrelated children who presented with neuroradiological findings and/or elevated lactate levels, highly suggestive of an underlying mitochondrial diagnosis. Next-generation sequencing identified bi-allelic variants in NDUFA6, encoding a 15 kDa LYR-motif-containing complex I subunit that forms part of the Q-module. Functional investigations using subjects’ fibroblast cell lines demonstrated complex I assembly defects, which were characterized in detail by mass-spectrometry-based complexome profiling. This confirmed a marked reduction in incorporated NDUFA6 and a concomitant reduction in other Q-module subunits, including NDUFAB1, NDUFA7, and NDUFA12. Lentiviral transduction of subjects’ fibroblasts showed normalization of complex I. These data also support supercomplex formation, whereby the ∼830 kDa complex I intermediate (consisting of the P- and Q-modules) is in complex with assembled complex III and IV holoenzymes despite lacking the N-module. Interestingly, RNA-sequencing data provided evidence that the consensus RefSeq accession number does not correspond to the predominant transcript in clinically relevant tissues, prompting revision of the NDUFA6 RefSeq transcript and highlighting not only the importance of thorough variant interpretation but also the assessment of appropriate transcripts for analysis.

Mutations in known disease genes account for the majority of autosomal recessive retinal dystrophies

Retinal dystrophies (RDs) are hereditary blinding eye conditions that are highly variable in their clinical presentation. The remarkable genetic heterogeneity that characterizes RD was a major challenge in establishing the molecular diagnosis in these patients until the recent advent of next-generation sequencing. It remains unclear, however, what percentage of autosomal recessive RD remain undiagnosed when all established RD genes are sequenced. We enrolled 75 families in which RD segregates in an apparently autosomal recessive manner. We show that the yield of a multigene panel that contains known RD genes is 67.5%. The higher yield (82.3%) when whole exome sequencing was implemented instead was often due to hits in genes that were not included in the original design of the panel. We also show the value of homozygosity mapping even during the era of exome sequencing in uncovering cryptic mutations. In total, we describe 45 unique likely deleterious variants (of which 18 are novel including one deep intronic and one genomic deletion mutation). Our study suggests that the genetic heterogeneity of autosomal recessive RD is approaching saturation and that any new RD genes will probably account for only a minor role in the mutation burden.

Bi-allelic mutations in TRAPPC2L result in a neurodevelopmental disorder and have an impact on RAB11 in fibroblasts

Background

The combination of febrile illness-induced encephalopathy and rhabdomyolysis has thus far only been described in disorders that affect cellular energy status. In the absence of specific metabolic abnormalities, diagnosis can be challenging.

Objective

The objective of this study was to identify and characterise pathogenic variants in two individuals from unrelated families, both of whom presented clinically with a similar phenotype that included neurodevelopmental delay, febrile illness-induced encephalopathy and episodes of rhabdomyolysis, followed by developmental arrest, epilepsy and tetraplegia.

Methods

Whole exome sequencing was used to identify pathogenic variants in the two individuals. Biochemical and cell biological analyses were performed on fibroblasts from these individuals and a yeast two-hybrid analysis was used to assess protein-protein interactions.

Results

Probands shared a homozygous TRAPPC2L variant (c.109G>T) resulting in a p.Asp37Tyr missense variant. TRAPPC2L is a component of transport protein particle (TRAPP), a group of multisubunit complexes that function in membrane traffic and autophagy. Studies in patient fibroblasts as well as in a yeast system showed that the p.Asp37Tyr protein was present but not functional and resulted in specific membrane trafficking delays. The human missense mutation and the analogous mutation in the yeast homologue Tca17 ablated the interaction between TRAPPC2L and TRAPPC10/Trs130, a component of the TRAPP II complex. Since TRAPP II activates the GTPase RAB11, we examined the activation state of this protein and found increased levels of the active RAB, correlating with changes in its cellular morphology.

Conclusions

Our study implicates a RAB11 pathway in the aetiology of the TRAPPC2L disorder and has implications for other TRAPP-related disorders with similar phenotypes.

TIMMDC1 Knockdown Inhibits Growth and Metastasis of Gastric Cancer Cells through Metabolic Inhibition and AKT/GSK3β/β-Catenin Signaling Pathway

TIMMDC1 (C3orf1), a predicted 4-pass membrane protein, which locates in the mitochondrial inner membrane, has been demonstrated to have association with multiple member of mitochondrial complex I assembly factors and core mitochondrial complex I subunits. The expression level of TIMMDC1 in highly-metastatic tumor cells is higher than that in lowly- metastatic tumor cells. However, the role of TIMMDC1 in human gastric cancer progression is unclear. In this study, human gastric cancer cells SGC-7901 and BGC-823 cells were used, and TIMMDC1 was knockdown with small interfering RNA. The data showed that TIMMDC1 knockdown caused inhibitory effects on the cell proliferation in vitro and tumor progression in vivo. Knockdown of TIMMDC1 significantly and exclusively reduced the activity of mitochondrial complex I but not complex II~ IV, and caused an obvious inhibition in mitochondrial respiration and ATP-linked oxygen consumption. Besides, the glycolysis pathway was also attenuated by TIMMDC1 knockdown, and the ATP content in the group of shTIMMDC1 cells was significantly lower than that in the shCont cells. The expression levels of phosphoylated AKT(Ser473) and GSK-3β (Ser9), as well as the downstream protein β-catenin and c-Myc were also markedly reduced in the group of shTIMMDC1 cells. Taken together, these findings suggest that TIMMDC1 may play an important role in human gastric cancer development, and its underlying mechanism is not only associated with mitochondrial complex I inhibition and reduced mitochondrial respiration, but is also associated with reduced glycolysis activity and the AKT/GSK3β/β-catenin signaling pathways.

Human diseases associated with defects in assembly of OXPHOS complexes

The structural biogenesis and functional proficiency of the multiheteromeric complexes forming the mitochondrial oxidative phosphorylation system (OXPHOS) require the concerted action of a number of chaperones and other assembly factors, most of which are specific for each complex. Mutations in a large number of these assembly factors are responsible for mitochondrial disorders, in most cases of infantile onset, typically characterized by biochemical defects of single specific complexes. In fact, pathogenic mutations in complex-specific assembly factors outnumber, in many cases, the repertoire of mutations found in structural subunits of specific complexes. The identification of patients with specific defects in assembly factors has provided an important contribution to the nosological characterization of mitochondrial disorders, and has also been a crucial means to identify a huge number of these proteins in humans, which play an essential role in mitochondrial bioenergetics. The wide use of next generation sequencing (NGS) has led to and will allow the identifcation of additional components of the assembly machinery of individual complexes, mutations of which are responsible for human disorders. The functional studies on patients' specimens, together with the creation and characterization of in vivo models, are fundamental to better understand the mechanisms of each of them. A new chapter in this field will be, in the near future, the discovery of mechanisms and actions underlying the formation of supercomplexes, molecular structures formed by the physical, and possibly functional, interaction of some of the individual respiratory complexes, particularly complex I (CI), III (CIII), and IV (CIV).

OCR-Stats: Robust estimation and statistical testing of mitochondrial respiration activities using Seahorse XF Analyzer

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.

Bi-allelic Mutations in Phe-tRNA Synthetase Associated with a Multi-system Pulmonary Disease Support Non-translational Function

The tRNA synthetases catalyze the first step of protein synthesis and have increasingly been studied for their nuclear and extra-cellular ex-translational activities. Human genetic conditions such as Charcot-Marie-Tooth have been attributed to dominant gain-of-function mutations in some tRNA synthetases. Unlike dominantly inherited gain-of-function mutations, recessive loss-of-function mutations can potentially elucidate ex-translational activities. We present here five individuals from four families with a multi-system disease associated with bi-allelic mutations in FARSB that encodes the beta chain of the alpha2beta2 phenylalanine-tRNA synthetase (FARS). Collectively, the mutant alleles encompass a 5'-splice junction non-coding variant (SJV) and six missense variants, one of which is shared by unrelated individuals. The clinical condition is characterized by interstitial lung disease, cerebral aneurysms and brain calcifications, and cirrhosis. For the SJV, we confirmed exon skipping leading to a frameshift associated with noncatalytic activity. While the bi-allelic combination of the SJV with a p.Arg305Gln missense mutation in two individuals led to severe disease, cells from neither the asymptomatic heterozygous carriers nor the compound heterozygous affected individual had any defect in protein synthesis. These results support a disease mechanism independent of tRNA synthetase activities in protein translation and suggest that this FARS activity is essential for normal function in multiple organs.

Mitochondrial medicine in the omics era

Mitochondria are dynamic bioenergetic organelles whose maintenance requires around 1500 proteins from two genomes. Mutations in either the mitochondrial or nuclear genome can disrupt a plethora of cellular metabolic and homoeostatic functions. Mitochondrial diseases represent one of the most common and severe groups of inherited genetic disorders, characterised by clinical, biochemical, and genetic heterogeneity, diagnostic odysseys, and absence of disease-modifying curative therapies. This Review aims to discuss recent advances in mitochondrial biology and medicine arising from widespread use of high-throughput omics technologies, and also includes a broad discussion of emerging therapies for mitochondrial disease. New insights into both bioenergetic and biosynthetic mitochondrial functionalities have expedited the genetic diagnosis of primary mitochondrial disorders, and identified novel mitochondrial pathomechanisms and new targets for therapeutic intervention. As we enter this new era of mitochondrial medicine, underpinned by global unbiased approaches and multifaceted investigation of mitochondrial function, omics technologies will continue to shed light on unresolved mitochondrial questions, paving the way for improved outcomes for patients with mitochondrial diseases.