Functional genomics and small molecules in mitochondrial neurodevelopmental disorders

Mitochondria are critical for brain development and homeostasis. Therefore, pathogenic variation in the mitochondrial or nuclear genome which disrupts mitochondrial function frequently results in developmental disorders and neurodegeneration at the organismal level. Large-scale application of genome-wide technologies to individuals with mitochondrial diseases has dramatically accelerated identification of mitochondrial disease-gene associations in humans. Multi-omic and high-throughput studies involving transcriptomics, proteomics, metabolomics, and saturation genome editing are providing deeper insights into the functional consequence of mitochondrial genomic variation. Integration of deep phenotypic and genomic data through allelic series continues to uncover novel mitochondrial functions and permit mitochondrial gene function dissection on an unprecedented scale. Finally, mitochondrial disease-gene associations illuminate disease mechanisms and thereby direct therapeutic strategies involving small molecules and RNA-DNA therapeutics. This review summarizes progress in functional genomics and small molecule therapeutics in mitochondrial neurodevelopmental disorders.


Introduction
Human neurodevelopment occurs in an orderly fashion from neurulation onto adulthood.Patterns of developmental attainment are wellcharacterized and include clinical (e.g., motor, language, social, cognitive) and biological (e.g., anatomic, cellular, network, biochemical) milestones.The precise timing and order of milestone attainment varies from individual-to-individual; significant delays or disordered patterns of milestone attainment are the defining feature of neurodevelopmental disorders.The precise definition of neurodevelopmental disorders varies.The most expansive definitions encompass 11-23.9% of the population and include common disorders like attention deficient hyperactivity disorder (ADHD), dyslexia, and learning disabilities [1].The spectrum of neurodevelopmental disorders ranges from common conditions to rare disorders with more substantial impacts on motor, behavioral, and cognitive function like autism, intellectual disability (ID), and cerebral palsy.While neurodevelopmental disorders may result from environmental factors or acquired brain injuries, genetic factors play an important role across the entire developmental disorder spectrum.Human genetics has demonstrated common disorders like ADHD are primarily complex or polygenic traits [2], whereas rare, severe neurodevelopmental disorders like intellectual disability are often monogenic or Mendelian traits [3].Mendelian neurodevelopmental disorders are highly heterogeneous [4,5].Their discovery rate has dramatically accelerated with advances in genomics; this is attested by the near doubling of ID-associated genes in a seven-year span (2016-2023) from approximately 800 to over 1500 [5].These discoveries provide answers for patients and families with neurodevelopmental disorders and frequently conclude a long 'diagnostic odyssey.'They also illuminate the complex biological pathways underlying human brain development.
Mitochondrial dysfunction is often implicated as a common pathway in neurodevelopmental disorders [6,7].Mitochondria are the "powerhouse of the cell" which generate adenosine triphosphate (ATP) through aerobic respiration.As only liver and muscle exceed the brain's energy demands, any failure of mitochondrial energy production can have disastrous neurological consequences [8].Mitochondrial energy production is particularly important during the neurodevelopmental transition from glycolysis-dependent neuronal stem cells to oxidative phosphorylation (OXPHOS)-dependent neurons [6].Mitochondria also have essential functions outside energy production including calcium buffering (Fig. 1).They are dynamic organelles migrating throughout neurons and undergoing continuous cycles of fusion and fission to replenish mitochondrial pools [7].Consequently, pathogenic variation in genes with an integral role in mitochondrial migration, fission and fusion like MSTO1, MFN2, OPA1, YME1L1, DNM1L, MFF and FBXL4 often manifests clinically in neurodevelopmental disabilities [9,10].
Each mitochondrion contain its own 16.5 kilobase genome (mtDNA) with 37 genes encoding 13 proteins, 22 transfer RNAs (tRNA), and 2 ribosomal RNAs (rRNA) [11,12].Unlike the nuclear genome, mtDNA is exclusively inherited from an individual's mother.Each cell contains multiple mitochondria and therefore multiple copies of the mitochondrial genome allowing for homoplasmy and heteroplasmy.Homoplasmy refers to a state of uniformity in which all copies of mtDNA contain a particular genetic variation.Heteroplasmy occurs when only a fraction of the multiple copies of mtDNA contains a particular genetic variation.The degree of heteroplasmy in separate tissues can increase or decrease throughout life depending on selective pressure; consequently, mitochondrial dysfunction may worsen or improve throughout the individual lifespan.RBCs have a short life cycle, and therefore mtDNA variation heteroplasmy levels may decrease over time secondary to mtDNA turnover and lead to difficulty with detection and diagnosis based on blood alone.Thus, muscle or liver testing may be required to detect heteroplasmic mtDNA variations.Maternal mtDNA inheritance, heteroplasmy, and homoplasmy also result in mitochondrial disease intrafamilial variability.An asymptomatic mother with low level heteroplasmy for a pathogenic mtDNA variant may have children with highly variable disease severity depending on each child's degree of heteroplasmy.While mtDNA is essential for mitochondrial function, mtDNA alone is insufficient for mitochondrial homeostasis.As the mitochondrial proteome is estimated to contain over 1500 proteins while the mitochondrial genome contains only 37 genes, most mitochondrial proteins are encoded by the nuclear genome [13,14].These nuclear encoded mitochondrial genes are often essential for mitochondrial function and typical neurodevelopment; consequently, individuals with mitochondrial neurodevelopmental disorders are more likely to have pathogenic variation in nuclear encoded genes than genes within the mitochondrial genome.

Traditional approaches and cellular models of mitochondrial neurodevelopmental disorders
There are a multitude of well-established laboratory and imaging biomarkers of mitochondrial neurodevelopmental disorders including elevated plasma and cerebrospinal fluid lactate levels, organic aciduria, low plasma carnitine, elevated alanine, deep gray matter signal hyperintensity on magnetic resonance imaging (MRI), and lactate peak of magnetic resonance spectroscopy (MRS) [15].Serum levels of growth factors GDF15 and FGF21 are elevated in some individuals with mitochondrial disorders but can be elevated in other clinical settings such as inflammation and acute illness [16][17][18].Other more sensitive and specific assays include quantitative assessment of OXPHOS enzyme activity, blue-native electrophoresis (BN-PAGE), mitochondrial copy number in tissue for DNA depletion, and gene-specific assays [12,15].
Cellular models to characterize the functional impact of mitochondrial genetic variation ex vivo include 1) cybrids, cell lines generated from the fusion of patient-derived mitochondria with a cell line lacking functional mitochondria; 2) patient-derived primary cells or immortalized cell lines including skin fibroblasts and lymphoblastoid cell lines; 3) patient-derived induced pluripotent stem cells; and 4) patient-derived induced neurons.Animal models including rodents, pigs, zebrafish, roundworms, and fruit flies have also been extensively used to characterize the impact of mitochondrial gene deletion or even human genetic variation through genetic variant knock-in at the organismal level [19,20].
Induced pluripotent stem cells (iPSCs) are produced by reprogramming primary cells like fibroblasts and peripheral blood mononuclear cells (PBMCs) into an embryonic stem cell-like state via the Yamanaka factors [21].iPSCs maintain pluripotency in culture and can be expanded and differentiated into many cell types within the developing brain including neuronal progenitor cells (NPC), neurons, astrocytes, oligodendroglial cells, and microglia.Abnormalities identified in genetic neurodevelopmental disorders have been recapitulated in iPSC-derived NPCs and neurons from patients with neurodevelopmental disorders as well as iPSCs genetically engineered to contain pathogenic variants [22].Thus, iPSC-derived neurons are a well-established human cellular model for the study of mitochondrial neurodevelopmental disorders with the potential to be used in high-throughput functional genomics studies.
The production, maintenance, and differentiation of iPSCs is costly, technically challenging, and time-consuming.An alternative ex vivo model of human neurons with the potential to be less costly and timeconsuming than iPSCs are inducible neurons (iNeurons) [23,24].iNeurons are produced via the induction of neuronal transcription factors in fibroblasts or PBMCs.They closely resemble primary neurons morphologically and express neuronal markers.They also fire action potentials and form mature synapses.To date, iNeurons have not received much study in neurodevelopmental disorders but have shown promise in the study of neurodegenerative diseases.

Genomic sequencing
The cornerstone of functional genomics is genomic sequencing (Fig. 2).Advances in genomic technologies over the past twenty years have enabled high-throughput genomic investigations in cellular models, model organisms, primary cells, and biospecimens.Short-read sequencing (SRS), also known as large-scale massively parallel sequencing or next generation sequencing (NGS), has transformed genomics from its introduction (Fig. 2A) [25,26].SRS generates a tremendous volume of short (~50-150 base pairs, bp) DNA sequence fragments known as reads.Each nucleotide is sequenced through multiple reads; genomic sequencing may involve anywhere from 10 to 100s of reads per nucleotide with ultra-high read depths used to study somatic mosaicism.Each read must be computationally mapped to the reference genome and assembled to generate gene or genome-level sequence.At present, SRS is the most widely used genomic sequencing approach.With adequate read depth, SRS has excellent sensitivity for detection of single nucleotide variants (SNVs), indels, some structural variants, somatic mosaicism, and tandem repeat expansions [27][28][29][30][31].
Another major milestone in genomics was the development of capture methods allowing for targeted sequencing of genomic regions of interest Fig. 2. Genomic approaches to mitochondrial neurodevelopmental disorders.A) Short-read sequencing (SRS, also known as next generation sequencing, NGS) is the most common technology used in genomic sequencing.In SRS, sequencing data consists of small (50-150 bp) sequences called reads.Short-read data must be computationally assembled to generate genomic sequences.Gaps in genome coverage (gray rectangles) in SRS can limit resolution of complex genomic regions (CGR) and structural variants (SV).In contrast, long-read sequencing technologies generate kilobase to megabase sequencing reads.Long sequencing reads provide better coverage without gaps across the genome.They improve SV and CGR resolution and played an integral role in the Telomere-to-Telomere (T2T) project.B) Exome sequencing (ES) involves the selective capture and sequencing of the exome.The exome is composed of all protein-coding genes and makes up ~1 % of the genome.Capture regions (purple rectangles) often include flanking intronic regions since intronic genetic variants near splicing junctions often cause splicing dysfunction.Genome sequencing (GS) involves sequencing of the entire genome including introns, promoters, and intergenic regions including regulatory elements like enhancers and suppressors.C) The advantages and disadvantages of SRS, LRS, ES, and GS are shown.SVs ¼ structural variants; CGRs ¼ complex genomic regions; SNVs ¼ single nucleotide variants; indels ¼ small insertion and deletions.Figure created using Biorender.com.
by SRS [32,33].Capture methods enabled exome sequencing (ES), the selective sequencing of all exons within the genome (Fig. 2B).Exons are the protein-coding portions of genes which are transcribed as messenger RNA (mRNA).The exome makes up approximately 1 % of the whole genome.ES is highly sensitive for the detection of coding SNVs and indels.It can also detect large copy number variants (deletions and duplications) and unbalanced translocations.ES has proven highly efficient at identifying pathogenic coding variation in mitochondrial disease genes encoded within the nuclear genome.
In recent years, genome sequencing (GS) has been increasingly utilized in research and clinical settings (Fig. 2B).GS encompasses the entire nuclear genome as well as the mitochondrial genome.Consequently, GS is more costly but more comprehensive than ES.GS is superior to ES for the detection of structural variants; it also detects noncoding variants in promoters, deep intronic regions, and other regulatory elements like enhancers and suppressors.Noncoding variants are increasingly recognized as a cause of Mendelian neurodevelopmental disorders [34,35].
Another recent advance is long-read sequencing technologies (LRS) like single-molecule real-time (SMRT) and Nanopore sequencing.LRS can produce sequence reads in excess of 10 kilobases [27,36].The major advantage of LRS over SRS is a significant reduction in sequencing gaps (Fig. 2A-C).This improves read mapping and thereby helps resolve repetitive complex genomic regions (CGRs) like Alu elements, long interspersed nuclear elements (LINE), segmental duplications, and centromeres.It also improves the precise resolution of structural variants.LRS can also determine the phase of variants at a single locus if both variants are close enough to permit their capture within single reads [37].LRS played an integral role in the Telomere-to-Telomere (T2T) project generating the first complete haploid genome assembly [38] and the human pangenome reference, a collection of 47 phased, diploid genome assemblies from a genetically diverse cohort [39].In addition, LRS platforms can simultaneously profile chromatin accessibility and methylation states.The major disadvantages of LRS at present are high cost, requirement for high molecular weight DNA to maximize sequencing read length, and higher error rate (Fig. 2C).Technical advances like circular consensus SMRT sequencing have improved error rates and accuracy of LRS [40].The application of LRS to rare Mendelian mitochondrial neurodevelopmental disorders is ongoing [41,42].Bulk genomic sequencing of tissues, biospecimens, or cell cultures should also be distinguished from single cell genomic sequencing.The recent application of single cell genomic sequencing to mitochondrial disorders and has demonstrated distinctive patterns of purifying selection across blood and immune cell lineages [43,44].
The interpretation of genetic variation has been facilitated by the development of genomic databases.Since mitochondrial neurodevelopmental disorders are rare because of negative selection due to their severe phenotypic consequences, their underlying pathogenic genetic variants are generally rare (allele frequencies <1 %).Common genetic variants (allele frequencies >1 %) within mitochondrial disease genes can be concluded to represent benign polymorphisms using the same logic.Therefore, estimations of variant allele frequencies within human populations are essential for variant interpretation.One major resource of variant allele frequencies is the Genome Aggregation Database (gnomAD) [45,46].gnomAD was launched in 2014 to aggregate exome and genome sequencing data for the clinical and research human genetics communities.The latest version v4.0.0 contains ES data from 730,947 individuals and GS data from 76,215 individuals.The dataset can be downloaded and used for variant annotation and filtering; it is also accessible via a user-friendly web portal to quickly assess the allele frequency of any variant of interest.Pathogenic de novo variants in individuals with sporadic severe neurodevelopmental disorders should be absent or extremely rare in population databases; higher allele frequencies are incompatible with a severe de novo disease model.Allele frequencies of pathogenic variants in genes associated with severe recessive diseases can be higher than those associated with severe dominant disorders, but they should not be present in the homozygous state.As some genetic variants can be extremely rare in one population but more common in others, the lack of diversity in genomic datasets is an important limitation which must be considered [47].In addition to allele frequencies, gnomAD also contains helpful gene-level metrics including constraint against missense and predicted loss-of-function (pLoF) variation and regional missense variation constraint.Genes with the highest degree of pLoF constraint (probability of being loss-of-function intolerant [pLI] >0.9) are enriched for haploinsufficiency.Haploinsufficiency refers a requirement for two functional copies of a gene for normal health: loss of one copy of a haploinsufficient gene due to a loss-of-function allele results in a disease phenotype.High levels of regional missense variation constraint correspond to critical functional domains like helicase domains and protein binding sites.Constraint metrics can help prioritize and interpret genetic variation in known and novel disease genes.Genetic variants in highly constrained genes or gene regions are anticipated to be more deleterious than genetic variants in less constrained genes or gene regions.
Other useful genomic databases include ClinVar, Human Gene Mutation Database (HGMD), Leiden Open Variation Database (LOVD), and MITOMAP [48][49][50][51].These databases aggregate human genetic variation reported in the medical and scientific literature and genetic variation identified in clinical diagnostic laboratories.Gene variant entries in ClinVar, HGMD, and LOVD include clinical classificationspathogenic, likely pathogenic, variant of uncertain significance, likely benign, and benignaccording to American College of Medical Genetics and Genomics criteria [52].The MITOMAP database specifically compiles genetic variation in the mitochondrial genome and thus can be useful to discriminate between benign mitochondrial polymorphisms and pathogenic genetic variants.Thus, genomic databases like gnomAD, ClinVar, and MITOMAP can be quickly queried for any variant to determine allele frequency and clinical significance if known.
In silico software prediction programs also play an important role in human genetic variant interpretation and functional genomics.These tools leverage evolutionary conservation, variant allele frequencies, and protein structure/function to predict a genetic variant's impact on gene function.Ensemble methods like CADD and REVEL incorporating multiple in silico tools and/or trained on known pathogenic and benign genetic variants have demonstrate superior performance over individual tools [53,54].More recently, programs like Alpha Missense, ESM1b, PrimateAI-3D, EVE, MAVERICK, SpliceAI and Pangolin have utilized unsupervised deep learning models to provide genome-wide predictions of variant effect [55][56][57][58][59][60][61].These models consistently outperform older in silico predictive models.

Transcriptomics
Despite advances in our understanding of regulatory elements, transcription, and RNA splicing, it remains difficult to predict the functional consequence of genetic variation within noncoding regions of the genome.This is problematic as pathogenic noncoding variation can alter gene expression or RNA splicing and is increasingly found to underlie Mendelian neurodevelopmental disorders [34,62].For example, genetic variants within promoters, enhancers, or long noncoding RNAs may increase or decrease gene expression and thereby cause disease.Similarly, genetic variation in deep intronic regions may disrupt intronic splicing enhancers and silencers and therefore cause aberrant inclusion of introns or exon skipping.While deep learning models like SpliceAI and Pangolin have improved predictions of splicing abnormalities resulting from coding and non-coding single nucleotide variants and indels [59,60,63], these in silico predictive models continue to have significant limitations demonstrated by numerous examples of pathogenic splicing variants not predicted by these programs [64,65].
RNA sequencing (RNA-seq) provides a relatively unbiased survey of all transcripts through which aberrant gene expression, aberrant splicing, and allelic dropout/imbalance can be detected (Fig. 3) [12,[65][66][67].Gene expression outliers within samples can be identified by comparison with control samples (Fig. 3A); increased or decreased gene expression may point to a genetic variant disrupting a promoter or enhancer.Genetic variants which introduce premature termination codons (PTC) may reduce gene expression via transcript degradation due to the quality control mechanism nonsense-mediated decay (NMD).Similarly, splicing outliers including exon skipping and intron inclusion can be detected (Fig. 3B); these outliers may occur due to genetic variants which destroy or introduce novel splice sites or alter exonic and intronic splicing enhancers and silencers.Finally, RNA-seq can detect the preferential expression of a single allele, so-called allelic imbalance or dropout (Fig. 3C).Allelic imbalance may result from genetic variants in regulatory regions or genetic variants which introduce PTCs and trigger NMD.
RNA-seq may be performed in bulk using tissue specimens like muscle biopsies or cell cultures or at the single cell level.Single-cell transcriptomics has provided an unprecedented window into the heterogeneity and complexity within tissues and cells and is a very active area of investigation [68].Single-cell transcriptomics has been used to study the signaling pathways and cell types involved in mitochondrial disease pathogenesis [69,70].
Bulk RNA-seq using patient blood samples, fibroblasts and muscle tissue has been shown to boost diagnostic rates in rare Mendelian disorders.Characterization of gene expression outliers, allelic dropout/ imbalance, and aberrant splicing via RNA-seq can help identify pathogenic genetic variants missed during the analysis of ES and GS data (Fig. 3) [42,[65][66][67][71][72][73].RNA-seq can also functionally validate intronic variant pathogenicity through demonstration of splicing abnormalities.For example, transcriptomics identified the mitochondrial disease gene TIMMDC1 through recognition of a homozygous deep intronic variant introducing a pseudoexon and premature termination codon [72].Moreover, RNA-seq can identify the genes and pathways impacted by gene knockout/knockdown or pathogenic genetic variation and thereby provide insights into the biology underlying mitochondrial neurodevelopmental disorders.Thus, transcriptomics is an important component of functional genomics toolkit.
Transcriptomic data interpretation remains challenging for several reasons.First, each tissue has distinctive patterns of gene splicing and expression.Comparisons between tissues are therefore inherently limited.Important genes and transcripts in the brain may not be present in clinically accessible tissues like blood, fibroblasts, or skeletal muscle.Second, Fig. 3. Transcriptomic approaches to mitochondrial neurodevelopmental disorders.RNA sequencing (RNA-seq) can identify gene expression outliers, splicing outliers, and allelic imbalance/ dropout within cells, tissues, or other biospecimens.A) Normal gene expression within cells or tissues is established by performing RNA-seq on control samples (blue bars).Each bar represents gene expression levels within a sample.Gene expression within the sample (red bar) is found to be an outlier.The lower expression within the sample indicates that a genetic variant reducing gene expression may be present within the gene or gene regulatory elements.B) RNA splicing results in the removal of intronic sequences and the splicing of exons (colored rectangles).Abnormal splicing patterns (splicing outliers) may result from pathogenic splicing variants; these can be detected in comparison to control samples.The pathogenic splicing variant in the figure results in exon skipping in half of transcripts (yellow box).C) The relative expression of two alleles at a particular locus can be quantified by examining how many reads contain a single nucleotide variant or polymorphism versus the reference sequence.Pathogenic variation within a gene or gene regulatory element that reduce one allele's expression will result in allelic imbalance or dropout as shown in the figure.sex, age, genotype, and other health conditions influence transcriptomics, and age-and sex-matched healthy control samples and data are seldom available.Third, standardized approaches to transcriptomic data analysis remain in development.Finally, nearly all transcriptomic studies to date have utilized short-read RNA-seq.The short sequencing reads of SRS cannot capture the entire coding sequence of most genes; therefore, transcript analysis requires computationally reconstruction of each transcript and is believed to lack many isoforms [74].LRS can capture the entire coding sequence of transcripts and may therefore provide a more accurate transcriptomic survey, but its high-cost limits coverage and may limit sensitivity for detection of low frequency but biologically meaningful isoforms [74].Alternative approaches to study splicing dysfunction like high throughput splicing assays are in development [64].

Proteomics
Biological complexity increases at each step along the central dogma: while the human genome contains approximately 20,000 genes, the human transcriptome and proteomes are estimated to contain approximately 300,000 transcripts and 10 6 -10 9 protein isoforms, respectively [75].The large jump in complexity from transcriptome to proteome reflects post-translational modifications, e.g., glycosylation, phosphorylation, ubiquitination, and protein cleavage.Comprehensive characterization of the proteome can be achieved through mass spectroscopy (MS) (Fig. 4) [76,77].MS permits the quantification of posttranslational modifications and protein isoforms within biospecimens, cells, and tissues.The combination of liquid chromatography (LC) or high-performance liquid chromatography (HPLC) to separate complex mixtures of proteins with MS is particularly suitable for high-throughput proteomics.
There are three main MS proteomic strategies: top-down, bottomup, and shotgun approaches [77].The top-down approach involves fragmentation of full-length proteins within the mass spectrometer followed by quantification of fragment sizes.In the bottom-up approach, proteins are enzymatically or chemically digested before performing MS.Finally, in the shotgun approach, complex solutions like serum, urine, and cell or tissue lysates are fragmented followed by HPLC before MS.
MS has also been applied to single-cell proteomics [78].Like single-cell transcriptomics, single-cell proteomics provides an opportunity to dissect cellular heterogeneity in mitochondrial neurodevelopmental disorders.However, single-cell proteomics faces significant technical challenges [78].While transcripts can be amplified by polymerase chain reaction (PCR) to facilitate single-cell transcriptomics, an analogous process for proteins does not exist.Thus, sample preparation must be carefully optimized to limit protein loss.The throughput of current single-cell proteomic approaches is also limited but is enhanced by multiplexed analysis.Other proteomic approaches including protein pathway array, proximity extension assay, and single-molecule protein sequencing exist but are beyond the scope of this review [77,79].
High-throughput proteomics thus has the potential to identify biomarkers useful for mitochondrial neurodevelopmental disorders diagnosis and assessment of response to therapeutic interventions.It also has the potential to identify biomarkers to inform the development of novel therapeutic modalities.Proteomics can also aid in the clinical diagnostic interpretation of single nucleotide variants, indels, and structural variants in mitochondrial neurodevelopmental disorders by providing a functional readout of protein stability and post-translational modification state [12].An important limitation of MS is that most missense variants do not alter protein stability and abundance.Unfortunately, the integration of unbiased, high-throughput proteomics into mitochondrial neurodevelopmental disorders or any rare Mendelian disorders has been limited to date.

Metabolomics
Metabolomics is the systematic characterization of small molecule metabolites within tissues, cells, or biospecimens like plasma, urine, or cerebrospinal fluid [80,81].Multiple analytical platforms exist for metabolomic studies including MS and nuclear magnetic resonance (NMR) [82].MS is widely used in metabolomics: it has greater sensitivity than NMR, is less expensive, and has a broader detection range.However, it requires sample manipulation prior to analysis and results in sample destruction.MS is typically coupled with liquid or gas chromatography to reduce complexity and improve sample throughput and sensitivity.Advantages of NMR include high reproducibility, rapid analytic time, and no need for sample manipulation prior to analysis.Additionally, NMR preserves samples for downstream analysis.The major downside of NMR compared with MS is lower sensitivity.The detection of low abundance compounds with NMR is particularly challenging.At present, neither technology can characterize all possible metabolites; the choice of technology must be tailored to each application [82].
There are two approaches to metabolomics: targeted and untargeted metabolomics [83,84].Targeted metabolomics strategies involve the quantification of a defined set of metabolites, whereas untargeted metabolomics quantify all metabolites within a sample.Targeted metabolomics encompasses traditional clinical diagnostic metabolic studies like plasma amino acids and urine organic acids.Targeted approaches generally result in the absolute quantification of metabolite levels and have reduced false positives and artifacts than untargeted approaches.Thus, targeted metabolomics is ideal for testing hypotheses about specific metabolites or metabolic pathways.Targeted metabolomics is not ideal if the metabolic pathway of interest is not fully understood.Furthermore, targeted metabolomics is typically limited to a relatively small number of metabolites.Untargeted metabolomics can provide data about 100s-1000s of small molecule metabolites within a sample and therefore can be very helpful for non-hypothesis driven investigational research.Disadvantages of untargeted metabolomics include semi-quantitative measurements of metabolites and greater susceptibility to false positives and artifacts than targeted metabolomics.
Untargeted metabolomic approaches using liquid chromatographymass spectroscopy (LC-MS) have recently entered the clinical realm [80].In patients with inborn errors of metabolism (IEM), untargeted metabolomics has superior diagnostic ability compared to traditional metabolic testing [85].Untargeted metabolomics has been applied to multiple biospecimens including plasma, urine, and cerebrospinal fluid.As many metabolites are unstable or labile, careful sample preparation and storage is critical.A systematic analysis of metabolomics utility in mitochondrial neurodevelopmental disorders has not been performed, but the broader range of metabolites surveyed may provide superior sensitivity over traditional metabolic testing.A distinctive metabolomic signature for NAXE and NAXD deficiency including abnormalities in nicotinamide branched-chain amino acids, and lipids was reported [86,87].Resolution of these abnormalities in NAXE and NAXD deficiency with niacin supplementation suggests metabolomics may also be a useful therapeutic biomarker in these conditions and perhaps other mitochondrial neurodevelopmental disorders.

Multiplexed assays of variant effect and saturation genome editing
Traditional approaches to study genetic variation's impact on mitochondrial function utilize patient samples or genetic manipulation of cell lines including transgene expression or gene editing.While these approaches offer important insights into mitochondrial gene function, they are not readily translated to high throughput screens capable of resolving the consequences of all possible genetic variation within a single gene, Fig. 5. Saturation Genome Editing and Multiplexed Assays of Variant Effect.Saturation Genome Editing (SGE) involves the introduction of all possible single nucleotide variants (SNVs) and indels into a gene or locus of interest (top half) into cell lines using CRISPR-Cas9.The experimental process is described in the bottom half of the figure.The cell line HAP1 is often used for genome editing since it contains a haploid genome.An SNV or indel library and CRISPR-Cas9 is introduced into cells to induce homology directed repair.Each cell acquires a single mutation.If the gene is essential for cellular survival, genetic variants resulting in a null allele will cause cell death.The cell population is subsequent sequenced by ES or GS.The dropout of variants within the population indicates that the variant is pathogenic.The degree of variant depletion is calculated and used to generate a functional score.Multiplexed assays of variant effect (MAVE) thus consist of the functional score of each variant assayed in an SGE study.locus, or gene family.The development and optimization of highly efficient genome editing techniques like CRISPR-Cas9 has finally made this a reality through multiplexed assays of variant effect (MAVE) and saturation genome editing (SGE) (Fig. 5) [88][89][90].SGE utilizes CRISPR/Cas9-mediated homology directed repair to introduce genetic variants into a gene or locus at scale.To date, SGE has primarily been performed using the haploid cell line HAP1.The introduction of gene-disrupting genetic variants into essential genes leads to their depletion in culture over time, whereas benign genetic variants persist.Each variant can thus be scored based on its persistence in culture.The approach is most sensitive for null alleles (e.g., loss-of-function variants), although hypomorphic alleles may cause slower but appreciable variant depletion [90].The functional scores for each variant form an atlas which helps resolve variants of uncertain significance (VUS) into benign versus pathogenic.SGE also provides a map of essential domains, amino acids, and nucleotides within genes.Although most SGE studies have focused on cellular viability, other approaches using high-throughput imaging readouts of mitochondrial content or morphology can be conceived [91][92][93].While MAVE studies of mitochondrial genes have not been published, several are ongoing (https://www.mavedb.org/).As many nuclear mitochondrial genes are essential, they should be good targets for SGE [94].
Unfortunately, the implementation of SGE for mitochondrial genome is technically challenging.Challenges in mitochondrial genome editing include 1) multiple copies of the mitochondrial genome per cell, 2) lack of homology directed repair in mitochondria, 3) the mitochondrial membrane as a barrier, and 4) RNA import into mitochondria [95].Recent reports of successful base editing of the mitochondrial genome, an approach which does not require homology directed repair, provide an important step toward mitochondrial genome SGE [96].

Complete functional annotation of the human genome requires global data aggregation
While the completion of the Human Genome Project in April 2003 was a major milestone in human genetics and medicine, it was only the beginning of vast, ongoing global effort to catalog and functionally annotate all possible human genetic variations.Due to rapid population growth over the last century, it is estimated nearly all possible single nucleotide variants exist within the human population [97].Thus, humanity has engaged in a massive saturation mutagenesis experiment [97].While models of human genetic variation in non-human organisms and cell lines have been helpful, the translation of human genetics into these systems is imperfect at best.The ultimate proof of a genetic variant's phenotypic impact is the recurrent association of a variant with a phenotype in humans (Fig. 6).
Allelic seriesthe collection of variant alleles at a gene or genomic locus which generate a range of phenotypic traitsare an essential component of functional genome annotation [98].Examples of mitochondrial genes with evolving and complex allelic series include HPDL, NAXD, and NAXE.HPDL encodes the mitochondrial protein 4-hydroxyphenylpyruvate dioxygenase-like [99][100][101].Although it shares significant sequence similarly with the tyrosine catabolic enzyme 4-hydroxyphenylpyruvate dioxygenase, it is not involved in tyrosine catabolism but instead functions within the CoQ10 biosynthesis pathway [102].Biallelic variants in HPDL cause a mitochondrial neurodevelopmental disease spectrum ranging from neonatal mitochondrial encephalopathy and Leigh syndrome to adolescent-onset, isolated hereditary spastic paraplegia (HSP) [99][100][101].While genotype-phenotype correlations were not initially apparent, subsequent expansion of the allelic series demonstrates HPDL loss-of-function variants associate with more severe phenotypes in keeping with a recessive disease model [101].Furthermore, a likely Latino founder allele, HPDL: c.859T>C p.Y287H, shows a strong genotype-phenotype correlation with simple HSP in the homozygous state (D.Calame, personal observation).NAXD and NAXE encode NAD(P)HX dehydratase and NAD(P)HX epimerase, respectively [103,104].The two enzymes form a two-step pathway for the detoxification of hydrated, toxic forms of NAD(P)H.As hydrated NAD(P)H accumulates during cellular stress, deficiency of either enzyme results in metabolic crises triggered by fever: progressive, early-onset encephalopathy with brain edema and/or leukoencephalopathy (PEBEL) [103][104][105].While initial publications described severe metabolic decompensation and death in infants or toddlers, the evolving NAXE and NAXD allelic series demonstrates a wider phenotypic spectrum with some individuals surviving crises, lacking cardinal clinical features like skin rashes, or even having adolescent-onset [86,87].While genotype-phenotype correlations remain under study, there is preliminary evidence that pathogenic variants early in NAXD may associate with milder phenotypes [86].A potential explanation for this genotype-phenotype correlation is NAXD encodes distinct cytosolic and mitochondrial isoforms; early mutations should only impact the cytosolic isoform and therefore preserve mitochondrial NAXD activity [86].
Functional human genome annotation is also uncovering genetic mimics of mitochondrial neurodevelopmental disorders.An illustrative example is SNAPC4 [106].SNAPC4 encodes a component of the small nuclear RNA (snRNA)-activating protein complex.The snRNA-activating complex directs snRNA transcription and is consequently required for spliceosome assembly and proper RNA splicing.Biallelic pathogenic variation in SNAPC4 decreases snRNA expression and leads to global splicing dysregulation.While most individuals with SNAPC4-related disorders have non-specific neurodevelopmental disorders characterized by developmental delay and movement disorders, some have severe infantile-onset neurologic dysfunction with bilateral basal ganglia T2 hyperintensity and a clinical diagnosis of Leigh syndrome.Other causes of bilateral infantile striatal necrosis without evidence of mitochondrial dysfunction include NUP62, VAC14, and ADAR1 [107][108][109].Identification of these mimics is critical.Phenotypic and genetic heterogeneity are likely a barrier to success in mitochondrial disease clinical trials; the exclusion of non-mitochondrial genetic mimics is needed to maximize the likelihood of trial success.Furthermore, treatments optimizing mitochondrial function or targeting mitochondrial dysfunction are unlikely to benefit patients with non-mitochondrial disorders.
Realizing the full promise of the human genomics at scale remains a massive undertaking.The functional annotation of common genetic variation through genome-wide association studies (GWAS) is nearly complete.As common variant studies have concentrated in Europe and the United States due to global inequities, remaining common variant studies must concentrate on the Global South and other understudied populations [110].In contrast, the study of rare genetic variation has only just begun.An explosion of neurodevelopmental disease-gene associations has resulted from family-based studies, next generation sequencing, and online gene matchmaking [4,[111][112][113][114][115][116][117][118][119][120][121][122][123][124][125][126][127][128], but genome-wide functional annotation requires greater international integration of human genomic and phenotypic data.National programs like UK Biobank and All of Us represent an important step in this direction [129,130], but a wealth of genomic and phenotypic data remains locked away in data silos: research databases, clinical diagnostic laboratories, and electronic health records.Ultimately, the complete functional annotation of the human genome requires greater global data sharing with societal commitment for the benefit of all.

Challenges inherent to multi-omic approaches
While considerable progress has been made towards the identification of the genes and genetic variants responsible for mitochondrial neurodevelopmental disorders, a definitive genetic diagnosis is obtained through ES or GS in less than half of patients [12,131].There are several potential explanations for inconclusive ES/GS [41].First, many disease-gene associations clearly remain to be discovered.The pace of gene discovery is rapid, yet as of December 8th, 2023, only 4868 of the approximately 20,000 computationally-annotated protein coding genes have phenotypic associations in the Online Mendelian Inheritance in Man (OMIM) database (https://omim.org/).Second, perceived phenotypic mismatch and phenotypic expansions can result in a failure to obtain a definitive diagnosis.Phenotypic data guides the analysis of genomic data as the genetic diagnosis is expected to match the clinical synopsis of the associated disease trait.In some cases, clinical phenotyping data is extracted into structured ontology terms (e.g., Human Phenotype Ontology) which then feeds into variant filtering algorithms [132][133][134].Therefore, incorrect or limited clinical phenotypic data may result in the filtering out of a pathogenic variant and therefore a "negative" result.Misinterpretation of phenotypic data can also play a role.For example, hypotonia is often misconstrued as a sign of neuromuscular (i.e., muscle or peripheral nerve) disease but is just as likely to have a central nervous system origin in mitochondrial neurodevelopmental disorders.Phenotypic expansion refers to the extension of the phenotypic spectrum associated with a particular gene or locus [135].It is now recognized most genes associate with a phenotypic spectrum ranging from mild to severe disease due to gene dosage effects in haploinsufficient or recessive disease models (e.g., hypomorphic missense variant allele versus null alleles) or allele-specific disease mechanisms (e.g., loss-of-function variants and haploinsufficiency versus missense variants and gain of function).However, the phenotypic spectrum can be so broad as to limit recognition of allelic disorders, e.g., biallelic HPDL variants associating with Leigh syndrome on the severe end of the spectrum and simple hereditary spastic paraplegia on the mild end.Third, technical limitations may prevent a definitive diagnosis.For example, ES may struggle to detect small copy number variants (e.g., single or multi-exon deletions and duplications) and will generally not detect pathogenic noncoding variants.Forth, there may be inadequate evidence to interpret the pathogenicity of a variant.This can reflect challenges in variant interpretation due to limitations of in silico predictive programs but also may be due to absence of family genomic data to determine de novo status or to test whether variants segregate according to Mendelian expectations.This is especially true for noncoding variation as in silico tools capable of evaluating their functional impact remain in a nascent stage.Furthermore, as population genomic databases which do not contain individual-level phenotypic data increase in size by ingesting data from large-scale programs like UK Biobank, All of US, 100, 000 Genomes Project, and Regeneron Genomics Center which do not screen out individuals with severe neurodevelopmental disorders, there may be new challenges in the interpretation of ultra-rare genetic variants in a Mendelian context [129,130,136,137].
There are many limitations to multi-omics including cost, a requirement for highly specialized knowledge, and logistical challenges.Collaboration and team science are a requirement for large-scale multiomic studies.Fortunately, costs continue to drop with technical advances.The cost of GS has dramatically declined from the time of the Human Genome Project (~$3 billion US dollars) to the present era of $100 genomes [27,138].Another important limitation is a lack of platforms for the integrated analysis of multi-omic data.The open source seqr software platform allows for integrative analysis of genomic and transcriptomic data; efforts to incorporate proteomic data are ongoing [139].Further advances are needed to fully realize the potential of functional genomics in mitochondrial neurodevelopmental disorders.

Small molecule therapeutics in mitochondrial neurodevelopmental disorders
We live in an exciting era.RNA/DNA therapies like antisense oligonucleotides (ASOs), small interfering RNAs (siRNAs), gene replacement therapies, and genome editing are entering the clinical realm with proven therapeutic benefit and life extension for rare, previously fatal Mendelian diseases [140,141].Yet with ~20,000 protein-coding genes and the considerable genetic heterogeneity of mitochondrial neurodevelopmental disorders, RNA/DNA therapeutics will remain inaccessible to most patients.However, it is important to recognize other therapeutic strategies like small molecules can be highly efficacious (Table 1).Indeed, small molecule therapeutics should be considered for all patients with mitochondrial neurodevelopmental disorders, even those without established small molecule treatments, by conceptualizing the gene defect within the context of the underlying biological pathway.Here we provide several examples of small molecule therapeutics for mitochondrial neurodevelopmental disorders.
Coenzyme Q10 deficiencies and HPDL-related disorders Coenzyme Q10 (CoQ10) is an integral cofactor in mitochondrial oxidative phosphorylation.Several genes are involved in CoQ10 biosynthesis (Table 1).Primary CoQ10 deficiency typically manifests with brain, muscle, and kidney involvement like encephalomyopathy, cerebellar ataxia, severe multisystem infantile disease, or steroidresistant nephrotic syndrome [142].CoQ10 supplementation results in clinical stabilization and even symptom reversal in some but not all individuals with primary CoQ10 deficiency [142].Intriguingly, hydroxyphenylpyruvate dioxygenase-like (HPDL) was recently found to produce the CoQ10 intermediate 4-hydroxymandelate (4-HMA) [102].As biallelic pathogenic variation in HPDL causes mitochondrial neurodevelopmental disorders ranging from Leigh syndrome to simple hereditary spastic paraplegia, supplementation of CoQ10 or 4-HMA might have therapeutic benefit in these conditions [99][100][101].A natural history study of HPDL-related disorders (ClinicalTrials.govID: NCT05848271) is ongoing with the end goal of a clinical trial of 4-HMA supplementation.
Thiamine, biotin, and SLC19A3/TPK1 deficiency SLC19A3 is a member of the solute carrier (SLC) family and encodes a thiamine transporter [143].Individuals with biallelic pathogenic SLC19A3 variants develop infantile or early childhood encephalopathy and Leigh syndrome [143].TPK1 encodes thiamine pyrophosphokinase, an enzyme which converts thiamine into the active cofactor thiamine pyrophosphate [144].Thiamine pyrophosphate deficiency similarly causes episodic encephalopathy and Leigh-like syndromes [144].While some patients with SLC19A3 or TPK1 deficiency respond to biotin supplementation, most evidence supports thiamine supplementation as the definitive treatment [144,145].Early intervention and long-term supplementation is critical.Neurologic injury following encephalopathic crisis may be irreversible, and recurrence after discontinuation of thiamine therapy has been reported [146].Substrate supplementation should be considered for other disorders of SLC transporters; increased substrate concentrations may boost substrate uptake through hypomorphic transporters or transport through alternative, low-affinity transporters.
Niacin and NAXE/NAXD deficiency NAXE and NAXD provide a two-step mechanism for the detoxification of hydrated NAD(P)H (NAD(P)HX) [103,104].Children with NAXE or NAXD deficiency develop encephalopathy and brain edema triggered by febrile illness or other metabolic stressors.Lactate in CSF or plasma may be elevated or normal.Pellagra-like skin lesions have observed in NAXE-related crises.The morbidity and mortality of crises is high.As NAXE/NAXD deficiency is expected to cause NAD(P)HX accumulation and NADþ pool depletion, trials of prophylactic niacin supplementation to replete NADþ pools has been performed in several individuals with prolonged survival without crisis recurrence and improvement in neurologic symptoms [86,87].Metabolomic testing during NAXE/NAXD crises supports the NADþ pool depletion hypothesis, and repeat testing post-niacin supplementation showed metabolic normalization [86,87].NADþ deficiency has also been identified in adults with mitochondrial myopathies, and niacin repletion increased blood and muscle NADþ levels and improved muscle strength and mitochondrial biogenesis [147].Systematic study of NADþ levels in individuals with mitochondrial neurodevelopmental disorders is needed to determine whether niacin supplementation might be beneficial more broadly in this population.

Jak inhibition and type 1 interferonopathies
ATAD3A encodes the ATPase family AAA-domain containing protein 3A, a nuclear-encoded mitochondrial membrane protein involved in mtDNA organization and mitochondrial dynamics [148].Monoallelic and biallelic pathogenic variants in ATAD3A cause mitochondrial neurodevelopmental disorders: autosomal dominant and recessive Harel-Yoon syndrome and autosomal recessive pontocerebellar hypoplasia [148][149][150][151]. Recent expansion of the ATAD3A allelic series demonstrates some individuals with Harel-Yoon syndrome due to dominant negative ATAD3A variants exhibit autoinflammation/autoimmunity and elevated type 1 interferon activity [152].Loss of ATAD3A in Harel-Yoon syndrome results in autoinflammation due to cytosolic mtDNA leakage.Cytosolic mtDNA is recognized as 'foreign' by DNA sensors cGAS/STING which activate the type 1 interferon pathway.A similar phenomenon has been identified in individuals with PNPT1-related disorders [153].Furthermore, pathogenic variation in SUPV3L1, a gene encoding an RNA helicase involved in the mitochondrial degradosome, results in double-stranded RNA (dsRNA) accumulation, type 1 interferon activation, mitochondrial dysfunction, and white matter disease [154,155].Like cytosolic mtDNA, dsRNA is recognized by cellular sensors including MDA5, RIG-1, TLR7/8.TLR3, PKR, and OAS1-3 which also trigger the type 1 interferon pathway.
Excess type 1 interferon is the hallmark of Aicardi-Gouiterres syndrome (AGS), a genetic neurodevelopmental autoinflammatory disorder resulting from dysregulated type 1 interferon expression and signaling [156].A direct overlap between AGS and severe mitochondrial neurodevelopmental disorders was suggested in a case report of a patient with AGS and a homozygous PNPT1 missense variant [157].As excess type 1 interferon drives neurodegeneration in AGS, it is possible but yet unproven that type 1 interferon may also contribute to the pathogenesis of mitochondrial neurodevelopmental disorders.This could potentially have therapeutic implications as Janus kinase inhibitors suppress type 1 interferon signaling and have demonstrated therapeutic benefit in AGS.Jak inhibition may also be beneficial for autoinflammation/autoimmunity in mitochondrial neurodevelopmental disorders.

Conclusions and perspectives
Progress in mitochondrial disease mechanisms and therapeutics has been hampered by genetic heterogeneity, an ever-widening clinical spectrum, and non-mitochondrial disease genetic mimics.The diagnostic rate in suspected mitochondrial neurodevelopmental disorders remains low [158].However, functional genomic approachesshort and long-read genomic sequencing, transcriptomics, proteomics, metabolomics, saturation genome editing, and global genotype-phenotype data aggregationrepresent a critical step towards the comprehensive functional annotation of all genes and genomic regions involved in mitochondrial disease.These advances suggest the bold prediction of the United States National Human Genome Research Institute (NHGRI) that "the biological function(s) of every human gene will be known" and "the clinical relevance of all encountered genomics variants will be readily predictable, rendering the diagnostic designation 'variant of uncertain significance (VUS)' obsolete" by 2030 may be achievable [159].Furthermore, the insights into mitochondrial disease pathogenesis resulting from these studies will provide a platform for development of better disease biomarkers, small molecule medicines, and RNA/DNA therapies.
While clinical genomics has brought precision diagnosis to many individuals with mitochondrial neurodevelopmental disorders, the implementation of precision medicine has been slow.Precision medicine refers to the development of individualized medical plan based on genotype, phenotype, and environment.A major barrier to precision medicine in mitochondrial neurodevelopmental disorders is their rarity.With many conditions identified only in a single family or a few families, it is not yet possible to identify genotype-phenotype correlations which may inform prognostication.Their rarity also provides a limit to the natural history data available to inform prognostication.Ascertainment bias also plays a role as comprehensive genetic investigations like ES only became available in the clinic in 2011 and are primarily performed on children in most healthcare systems; therefore, there is a paucity of clinical records available for adults to inform our understanding of the long-term trajectory and sequelae of most mitochondrial neurodevelopmental disorders.
One way to address precision medicine in mitochondrial neurodevelopmental disorders is through the implementation of artificial intelligence (AI) and machine learning (ML).The integration of large datasets like electronic health records, genomic databases, and highthroughput multi-omic investigations provides an ideal opportunity for the application of AI/ML to mitochondrial neurodevelopmental disorders.AI/ML approaches have already been applied to metabolomic data within UK Biobank to identify correlations between metabolomic states and common diseases [160].AI programs have also been shown to accurately predict risk of cardiac disease and cancer [161].The rich datasets generated for individuals with mitochondrial neurodevelopmental disordersneuroimaging, clinical phenotyping, genotyping, and potentially multi-omic technologiesmay provide an opportunity to use AI to improve precision medicine.Finally, AI may provide opportunities to quickly identify biological pathways and FDA-approved drugs for drug repurposing in mitochondrial neurodevelopmental disorders [162].A future in which precision medicine is a reality for all patients with mitochondrial neurodevelopmental disorders is on the horizon.

Fig. 1 .
Fig. 1.Physiological functions of mitochondria and examples of mitochondrial dysfunction in neurodevelopmental disorders.Physiologic functions of mitochondria are shown in the left column.Examples of mitochondrial dysfunction responsible for mitochondrial neurodevelopmental disorders are shown in the right column.Figure created using Biorender.com.
Fig.3.Transcriptomic approaches to mitochondrial neurodevelopmental disorders.RNA sequencing (RNA-seq) can identify gene expression outliers, splicing outliers, and allelic imbalance/ dropout within cells, tissues, or other biospecimens.A) Normal gene expression within cells or tissues is established by performing RNA-seq on control samples (blue bars).Each bar represents gene expression levels within a sample.Gene expression within the sample (red bar) is found to be an outlier.The lower expression within the sample indicates that a genetic variant reducing gene expression may be present within the gene or gene regulatory elements.B) RNA splicing results in the removal of intronic sequences and the splicing of exons (colored rectangles).Abnormal splicing patterns (splicing outliers) may result from pathogenic splicing variants; these can be detected in comparison to control samples.The pathogenic splicing variant in the figure results in exon skipping in half of transcripts (yellow box).C) The relative expression of two alleles at a particular locus can be quantified by examining how many reads contain a single nucleotide variant or polymorphism versus the reference sequence.Pathogenic variation within a gene or gene regulatory element that reduce one allele's expression will result in allelic imbalance or dropout as shown in the figure.Figure created using Biorender.com.

Fig. 4 .
Fig.4.Proteomic approaches to mitochondrial neurodevelopmental disorders.Summary of a bottom-up proteomic study.Proteomics can be performed on biospecimens, tissue samples, or cell lines (sample selection).Samples are then enzymatically or chemically digested to generate short peptides.These peptides then undergo MS, LC-MS, or HPLC-MS.Raw data then undergoes peptide-spectrum matching to identify the sequence of each peptide.Individual peptides are then quantified, and peptide fragment sequences are compared with the sequences of known proteins; the quantified proteins are then biologically interpreted.Figure created using Biorender.com.

Fig. 6 .
Fig. 6.Functional annotation of the human genome through global data aggregation.The integration and aggregation of individual-level genomic and phenotypic data derived from health records (top half) on a global scale is necessary to provide full functional annotation of the human genome.Such efforts are necessary to fully realize the promise of precision medicine through accurate prediction of mitochondrial dysfunction and mitochondrial neurodevelopmental disorders from an individual's genotype (bottom half).Figure created using Biorender.com.

Table 1
Small molecule therapeutics for mitochondrial neurodevelopmental disorders.
OMIM, Online Mendelian Inheritance in Man (www.omim.org);?Theoretical small molecule therapeutics with limited supporting clinical data.