Introduction

Human mitochondrial DNA (mtDNA) consists of a circular genome of 16,569 base pairs (bp) that contains 13 protein-encoding genes traceable to primordial maternal haplogroups (Anderson et al. 1981; Taanman 1999). The subunits ND1-3, ND4, ND4-6 of Complex I, cytochrome b, subunits CO1-3 of Complex IV, and subunits A6 and A8 of Complex V (ATP synthase) expressed from mtDNA represent essential components of the electron transport chain. Nuclear diploid gene sequences encode the vast majority of cognate protein subunits, accessory proteins, and assembly-related protein chaperones. Complex regulatory factors driving mtDNA gene expression in response to physiological or pathological changes in cellular energy demands are reciprocally dynamic, and they functionally extend to a multiplicity of nuclear target genes (Quiros et al. 2016). Within this mechanistic context, heteroplasmy of mtDNA specifically refers to dynamically determined co-expression of wild type (WT) inherited polymorphisms and collective time-dependent somatic mutations within individual mtDNA genomes (see Stefano et al. 2017).

The total number of mitochondria and mtDNA copy numbers vary by at least three orders of magnitude within the individual cell and tissue types (Allen 2015; He et al. 2010; Jacobi et al. 2001; Stefano et al. 2017). The patterns of expression and distribution of heteroplasmic mtDNA may represent a regulatory mechanism functionally linked to the maintenance of ATP levels in response to physiological demand (Allen 2015; He et al. 2010; Jacobi et al 2001; Stefano et al. 2017). The review presents interrelated topics that integrate recently elucidated mechanistic aspects of mtDNA heteroplasmy with special reference to mitochondrial associated immune processes across multiple cell types that include myeloid cells, monocytes, macrophages, and antigen-presenting dendritic cells. Our overall goal is to critically discuss and integrate convergent lines of empirical evidence that potentially link expression of mtDNA heteroplasmy by functionally interactive CNS cell types to the extent and severity of acute and chronic post-COVID-19 neurological disorders.

Functional Relevance of mtDNA Heteroplasmy to Human Health and Disease

Previously we have presented a functionally-oriented, multifaceted review of mitochondrial heteroplasmy from cellular, physiological, and pathophysiological perspectives (Stefano et al. 2017). Many studies have focused on identifying the potential etiological roles of heritable mutations in protein-encoding genes and mtDNA genes with sporadic somatic mutations (Davila and Zamorano 2013; He et al. 2010; Stefano et al. 2017). These mutations have been identified in congenital cardiomyopathies, neurological diseases, and human malignancy (Davila and Zamorano 2013; He et al. 2010; Stefano et al. 2017). The maintenance of adequate cell energy within dynamically changing physiological environments may be compromised following the loss of function mutations in protein-encoding mtDNA genes (Aryaman et al. 2019; Hoitzing et al. 2019). Accordingly, we have surmised that cell-specific and tissue-specific heteroplasmic mtDNA expression in healthy individuals is functionally associated with bidirectional intracellular communication involving mitochondrial signaling pathways and nuclear DNA gene expression (Stefano et al. 2017).

Stochastic modeling algorithms have been developed to provide theoretical frameworks that functionally link variability in cellular distributions of heteroplasmic mtDNA within mitochondrial networks (Aryaman et al. 2019; Hoitzing et al. 2019). These analyses suggest that complex regulatory factors underlying fusion and fission ratios are functionally associated with state-dependent variance in the mean distribution of heteroplasmic mtDNA copy numbers (Aryaman et al. 2019; Hoitzing et al. 2019). Accordingly, empirically determined tissue-specific variances in the rates of formation of heteroplasmic mtDNA distributions may represent biomarkers of optimized mitochondrial dynamics (Aryaman et al. 2019; Hoitzing et al. 2019). Within a clinical context, prospective studies on the homeostatic regulatory factors that underlie the observed intrinsic differences in cell-specific and organ-specific distribution of heteroplasmic mtDNA may result in novel treatment approaches for mitochondrial disease or enhanced quality of life for aging populations.

Theoretical modeling studies are complemented and empirically validated by in vivo analyses of engineered heteroplasmic mouse strains with nonpathological mtDNA polymorphisms (Lechuga-Vieco et al. 2020). In heteroplasmic mice, integrative analysis showed non-random segregation of mtDNA copy numbers in all tissues studied (Lechuga-Vieco et al. 2020). The authors concluded that mtDNA heteroplasmy reflected an intracellular sorting process based on organelle selection that was critically dependent on the functional coupling of the rate of mitochondrial oxidative phosphorylation to cellular metabolic requirements (Lechuga-Vieco et al. 2020).

Functional Relevance of mtDNA Heteroplasmy to Mitochondrial Transfer Mechanisms

In light of the above, the combination of stochastic modeling and complementary studies of heteroplasmic mouse strains with nonpathological mtDNA polymorphisms has our increased understanding of the potential depth of biological roles linked to mtDNA heteroplasmy (Aryaman et al. 2019; Hoitzing et al. 2019; Lechuga-Vieco et al. 2020). Within the last decade, key studies have supported the physiological importance of state-dependent transfer of functional mitochondria from healthy donor cells to metabolically compromised cells (Stefano and Kream 2016). Operationally, advantageous mitochondrial transfer indicates a superimposition of functionally competent heteroplasmic mtDNA copy numbers within the cell to restore mitochondrial dynamics within multiple organ systems (Jackson et al. 2016; Kitani et al. 2014; McCoy-Simandle et al. 2016; Tatsuta et al. 2014; Wang and Gerdes 2015). In cell-specific or organ-specific systems, intercellular transfer of functional mitochondria involves synergistic mechanisms involving tunneling nanotubes and/or cellular-derived exosomes as delivery systems (Jackson et al. 2016; Kitani et al. 2014; McCoy-Simandle et al. 2016; Tatsuta et al. 2014; Wang and Gerdes 2015). We contend that the summated effects of mass transfer of physiologically competent mitochondria in response to physiological demands may involve restorative transcriptional and translational processes mediated by functional heteroplasmic mtDNA copy numbers within cytoplasmic microdomains.

From translational perspectives, continuation of integrated biochemical and molecular studies designed to critically probe complex mechanisms of intercellular transport of functional mitochondria are required to critically validate the safety and efficacy of using functionally viable mitochondria as a novel therapeutic approach to restore cellular and organ function in disease and healthy aging (Stefano and Kream 2016; McCully et al. 2016). Recent advances in the use of selective CRISPR gene editing and functional reversal of heritable mtDNA mutations in mitochondrial diseases have supported the possibility of using these molecular bioengineering technologies (Mok et al. 2020). Accordingly, future therapeutic usage of mitochondrial replacement therapy employing highly purified preparations of autologous mitochondria previously engineered to provide optimal oxidative phosphorylation and ATP production represents a potentially novel approach for treating heritable mitochondrial diseases and well-studied metabolic disorders such as Type II diabetes. Similar considerations may apply to the potential future use of mitochondrial replacement therapy with selectively engineered mtDNA copy numbers for targeted therapy or personalized medicine in genetic, infectious, and degenerative diseases (Pacheu-Grau et al. 2013). In these therapeutic applications, tissue-specific and organ-specific distribution of heteroplasmic mtDNA copy numbers may provide quality control parameters for the successful application of mitochondrial replacement therapy strategies.

Finally, intact cell-free mitochondria with normal energy production have been identified in healthy human blood samples (Al Amir Dache et al. 2020). In this study, the authors estimated that between 200,000 and 3,700,000 functional mitochondria per ml of extracted plasma were present (Al Amir Dache et al. 2020). They proposed that circulating cell-free functional mitochondria were representative of a novel class of signaling organelles involved in complex regulatory activities and intercellular communication (Al Amir Dache et al. 2020). However, it remains to be confirmed whether blood-borne functional mitochondria serve as sentinels for maintaining homeostatic metabolic activities and a reserve for essential cellular functions.

Modulation of Interactive Immune Responses by mtDNA Heteroplasmy

A major goal of our research has been to investigate the functional roles of key mitochondrial regulatory molecules in the pathophysiology of acute and chronic inflammation, with special reference to altered expression of proinflammatory and anti-inflammatory cytokines and their protein receptor complexes (Esch and Stefano 2002; Esch et al. 2020; Stefano et al. 2015). A physiological state of proinflammation is a key autoregulatory mechanism that maintains immune surveillance of infective microorganisms, including viruses (Esch and Stefano 2002; Esch et al. 2020). In support of this contention, myeloblast-derived immune cells, monocytes, macrophages, and dendritic antigen-presenting cells undergo distinct programmed metabolic changes in mitochondrial bioenergetics in innate and adaptive immunological processes (Angajala et al. 2018; Zuo and Wan 2019). Polarized M1 macrophages utilize uncoupled mitochondrial oxidative phosphorylation functionally coupled to a disrupted cell cycle to maintain a proinflammatory phenotype (Angajala et al. 2018; Zuo and Wan 2019). M2 macrophages retain concerted mitochondrial systems to produce anti-inflammatory responses (Angajala et al. 2018; Zuo and Wan 2019). Immunologically-driven mitochondrial processes are associated with morphological changes and altered ratios of fusion and fission linked to mitophagy rates and junctional coupling to the endoplasmic reticulum (Angajala et al. 2018; Zuo and Wan 2019). The potential regulatory roles of heteroplasmic mtDNA copy numbers to successfully drive these integrative cellular processes have been discussed above ((Aryaman et al. 2019; Hoitzing et al. 2019; Lechuga-Vieco et al. 2020).

Recent studies have observed multiple proinflammatory effects from extracellular transport and uptake of mitochondrial degradation products from sites of cell and tissue damage or high intracellular mitophagy rates (Pollara et al. 2018; Zhu et al. 2018; Nagakawa et al. 2020). Diverse mitochondrial degradation products have been qualitatively termed mitochondrial damage-associated molecular patterns (mtDAMPs) and include N-formylated peptide fragments derived from mitochondrial proteins, intact full-length mtDNA copy numbers, mtDNA fragments, and small organic molecules such as ATP and cardiolipin that are functionally linked to ongoing innate immune responses to pathogens (Al Amir Dache et al. 2020; Esch et al. 2020; Chandel 2014; Tiku et al. 2020; West et al. 2011). In healthy individuals, ongoing cell turnover may result in a basal level of blood-borne mtDAMPs and circulating mitochondria to maintain a constitutive proinflammatory state that is significantly enhanced after immunological challenges by infectious pathogens. Interestingly, intact mitochondria released from apoptotic or necrotic dying cells have been observed to directly mediate proinflammatory processes with downstream effects on adaptive immune responses that may promote long term pathophysiological effects under chronic conditions (Pollara et al. 2018; Zhu et al. 2018; Nagakawa et al. 2020).

In light of the above, a recent study has provided empirical evidence in support of novel signaling activities of blood-borne mitochondria to mediate innate immune cell-to-cell communication (Song et al. 2020). Intact and cell-free mitochondria from human plasma were shown to present membrane-associated, immunologically active proteins, including CD270 and programmed cell death-ligand 1 (PD-L1) (Song et al. 2020). Furthermore, intact mitochondria from human plasma and immunologically active surface proteins were associated with upregulation of activated CD4 + T-cells and CD8 + T-cells and reduced concentrations of proinflammatory cytokines (Song et al. 2020). Additional studies are warranted to critically evaluate the potential regulatory roles of differences in the distribution of heteroplasmic mtDNA copy numbers within blood-borne metabolically competent mitochondria in the partial mediation of innate immune responses.

Modulation of the Immune Response Involving the Gut-Brain Axis by Heteroplasmic mtDNA

Communication between intrinsic luminal and immune cell populations and commensal strains of gut microbiota maintain homeostatic innate and adaptive immune processes within the gastrointestinal (GI) tract (Stefano et al. 2018). Recent studies have identified a pivotal role of host cell mitochondria and immune signaling to maintain optimal diversity and functional integrity of gut bacteria commensal strains (Hirose et al. 2017; Saint-Georges-Chaumet and Edeas 2016; Yardeni et al. 2019). Mitochondrial oxidative phosphorylation is functionally coupled to a disrupted mitochondrial tricarboxylic acid cycle (TCA) cycle and is associated with a proinflammatory phenotype in enteric immune cells, mediated by reactive oxygen species (ROS) (Saint-Georges-Chaumet and Edeas 2016; Yardeni et al. 2019). Enteric microglia found in Auerbach's plexus and Meisner’s plexus of the enteric nervous system are key players in maintaining innate immunity in the GI tract (Verheijden et al. 2015). The sustained release of proinflammatory mediators within the gut converts enteric microglia into active immune cells with an M1 macrophage-like phenotype (Sonetti et al. 1994; Stefano et al. 1994).

As previously discussed, activation of innate immune responses to infection by pathogens is dependent on functionally competent mitochondria. Within the enteric nervous system and GI tract, a complementary array of normative mitochondrial functions is adversely affected during infection by bacterial pathogens (Lobet et al. 2015). Repeated exposure to outer membrane vesicles from Gram-negative bacteria may induce mitochondrial dysfunction and apoptosis functionally linked to the activation and sustained release of interleukin-1β to enhance proinflammatory processes in the GI tract (Deo et al. 2020). Recently reported studies from a mouse model evaluated the downstream effects of selective single nucleotide mutations in mitochondrial respiratory complex protein genes expressed by mtDNA open reading frames (ORFs) on changes in mitochondrial function and immune processes and differences in populations of gut microbiota (Hirose et al. 2017, 2019; Schilf et al. 2021). In sum, these data gleaned from several complementary studies on heteroplasmic mouse strains strongly suggest an important regulatory role of mtDNA heteroplasmy within mitochondria of luminal and immune gut cells to maintain mitochondrial bioenergy homeostasis and innate immune signaling processes.

Bidirectional gut-brain signaling pathways are reciprocally regulated by metabolic and immune processes of the host GI cell mitochondria (Diaz Heijtz et al. 2011; Zhou and Foster 2015). Conversely, dysbiosis of major strains of commensal gut microbiota has been functionally associated with degenerative CNS disorders that include Parkinson’s disease (Bedarf et al. 2017; Devos et al. 2013; Forsyth et al. 2011; Malkki 2017; Mulak and Bonaz 2015; Rietdijk et al. 2017; Scheperjans et al. 2015), Alzheimer's disease, Huntington's disease, and amyotrophic lateral sclerosis (ALS) (Main and Minter 2017; Tremlett et al. 2017). Several comorbid pathophysiological events have been shown to precede classic motor symptoms in Parkinson’s disease and include dysbiosis of the gut microbiota, loss of intrinsic enteric nervous system and lymphoid coupling, and progressive colonic inflammation (Bodea et al. 2014; Braak et al. 2003; Dobbs et al. 1999; Forsyth et al. 2011; Kosikowska et al. 2016; Lebouvier et al. 2009; Rietdijk et al. 2017; Villaran et al. 2010). It is therefore possible that the pathophysiology of Parkinson’s disease involves the slow progression of degenerative changes in the gut transmitted by peripheral nerves and humoral regulatory loops to the CNS with degeneration of the nigrostriatal CNS pathway (Braak and Del Tredici 2017; Braak et al. 2003; Rietdijk et al. 2017). More recently, the possible functional roles of heritable polymorphisms in mtDNA open reading frames (ORFs) have been investigated as a source of mitochondrial dysfunction in Parkinson’s disease-associated neurodegenerative processes (Müller-Nedebock et al. 2021).

Modulation of Mitochondrial, Metabolic, and Immune Responses in COVID-19 Neurological Disorders by Heteroplasmic mtDNA

A review of complementary in vivo, in vitro and neuropathological studies has highlighted convergent and potentially synergistic strategies used by SARS-CoV-2 for neuro-invasion, infection, and replication within CNS neurons and glia with subsequent cellular and tissue damage linked to multifaceted neurological dysfunction (Crunfli et al. 2020; Matschke et al. 2020; Song et al. 2021). Song and coworkers utilized a three pronged self-validating approach to evaluate the mechanistic capacity of SARS-CoV-2 to infect the CNS via in vitro usage of human brain organoids, in vivo modeling employing transgenic mice overexpressing human ACE2, and analysis of human autopsy specimens from patients who died from COVID-19-associated pathologies. Overall, the authors provided evidence for the neuroinvasive capacity of SARS-CoV-2 in an apparent ACE2-dependent mechanism for cellular entry and replication. Furthermore, they empirically demonstrated that CNS structures are sites for high replicative potential for SARS-CoV-2 and that SARS-CoV-2 infection subsequently engenders significant neuronal death in the in vitro human brain organoids model. Examination of dying cells in the organoid cultures indicated that neuronal death did not directly colocalize with viral infection, thereby suggesting that infected cells are capable of engendering regions of hypoxia within microdomains that potentially lead to tissue damage due to metabolic insufficiencies. Importantly, these observations are line with current critical thinking that links COVID-19 progression with successful takeover or hijacking of essential mitochondrial processes, including oxidative phosphorylation and ATP production (Burtscher et al. 2020; Stefano et al. 2020, 2021; Wang et al. 2020). Accordingly, successful SARS-CoV-2 infection is proposed to be functionally associated with a functional override of antiviral innate immune responses associated with multiple mitochondrial processes (Burtscher et al. 2020; Shenoy, 2020).

As discussed above, enteric microglia found in Auerbach's plexus and Meisner’s plexus of the enteric nervous system are key players in maintaining innate immunity in the GI tract (Verheijden et al. 2015) and sustained release of proinflammatory mediators within the gut converts enteric microglia into active immune cells with an M1 macrophage-like phenotype (Sonetti et al. 1994; Stefano et al. 1994). In similar fashion, neuroinflammatory conditions within the CNS induce brain microglia to undergo mitochondrial-dependent metabolic reprogramming from oxidative glycolysis towards aerobic glycolysis (Bernier et al. 2020; Lynch 2020). Aberrant neuroimmune responses engendered by polarized microglia that are functionally linked to widespread metabolic dysregulation appear to represent common components of neurodegenerative diseases, COVID-19 neurological disorders and coordinate damage to CNS neurons and astrocytes (Crunfli et al. 2020; Matschke et al. 2020; Song et al. 2021). Interestingly, extracellular release of functionally competent mitochondria from astrocytes coupled to uptake by damaged neurons utilizing an in vivo rodent stroke model has been observed to promote enhanced cell survival signaling and neurological outcome (Hayakawa et al. 2016). As discussed above, additional studies are warranted to critically evaluate the potential regulatory roles of differences in the distribution of heteroplasmic mtDNA copy numbers within metabolically competent mitochondria in damaged CNS neurons and glia subsequent to COVID-19 infection.

A high resolution predictive index of the extent of SARS-CoV-2 neuro-invasion may be obtained from a computational modeling study that localized full-length SARS-CoV-2 genomic RNA and relatively small sub-genomic RNA (sgRNA) transcripts across eight subcellular regions, with enrichment, observed in the mitochondrial matrix and nucleolus (Wu et al. 2020). Comparative analysis showed that SARS-CoV-2 RNA exhibited significantly stronger mitochondrial matrix and nucleolar residency signals than SAR-CoV-1 and related coronaviruses (Wu et al. 2020). Notably, it was hypothesized that the presence of SARS-CoV-2 in the mitochondrial matrix was a strong indicator of local translation of sgRNA transcripts and the expression of inhibitory proteins, including the open reading frame ORF9b (Wu et al. 2020). ORF9b was shown to directly interact with mitochondrial TOMM70 to downregulate host cell immunity and facilitate viral replication (Wu et al. 2020). These authors hypothesized that the high degree of similarity between endoplasmic reticulum-derived double-membrane vesicles (DMVs) required for SARS-CoV-2 replication and the double membrane of the mitochondrial matrix might be indicative of a common mechanism of intracellular viral RNA transfer and functional sequestration (Wu et al. 2020). Subsequent reviews have hypothesized a link between SARS-CoV-2 infection and replication that requires successful takeover of mitochondrial processes (Singh et al. 2020; Wang et al. 2020). The findings from several studies that have been reviewed support multiple sites of action of viral proteins expressed by SARS-CoV-2 sgRNA transcripts to suppress essential mitochondrial functions resulting in down-regulation of innate and adaptive immune processes (Singh et al. 2020).

Accordingly, a unified hypothesis of the mechanism of action for functional entrapment, or takeover, of host mitochondria, involves targeted, potentially synergistic, effects of small viral proteins expressed from sgORFs or sgRNA transcripts at critical regulatory sites within the mitochondrial matrix. The scope of potential inhibitory effects mediated by small viral proteins within the mitochondrial matrix may extend to protein complexes composed of subunits expressed from both normal DNA and mtDNA. Whereas the findings of the computational analyses by Wu and coworkers await secondary validation by complementary empirical studies, these observations have implications for advancing the understanding of the systemic harmful effects of SARS-CoV-2 infection in patients with COVID-19. Infection with SARS-CoV-2 is both organ-specific and systemic. The presence of intact mitochondria that contain SARS-CoV-2 RNA and viral proteins may act as infectious units or as complex inflammatory vehicles. The subsequent regulatory roles of intercellular and blood-borne transport and uptake of functional mitochondria and mtDNA may maintain state-dependent levels of proinflammatory effects.

In COVID -19, a paradoxical viral infection of circulating platelets lacking the required angiotensin-converting enzyme-2 (ACE-2) receptor has been reported in severely ill patients presenting with thrombotic complications associated with SARS-CoV-2 infection (Manne et al. 2020). ACE2-independent uptake of fully functional SARS-CoV-2 genomic mRNA in circulating platelets as an alternate mechanism of viral infection appears to be consistent with the predictive existence of circulating blood-borne mitochondria carrying resident full-length genomic RNA sequestered within the mitochondrial matrix (Wu et al. 2020). Patients with COVID-19 who have thrombotic complications of SARS-CoV-2 infection have been reported to transiently express antiphospholipid autoantibodies, including mitochondrial-derived anticardiolipin immunoglobulins (Dotan et al. 2021). This functional observation supports that intact circulating mitochondria with resident biologically active species of SARS-CoV-2 RNA or small viral proteins act as self-contained complex inflammatory vehicles. The multiple effects of heteroplasmic mtDNA copy numbers to regulate normative mitochondrial dynamics may predict the alteration of the residency patterns of SARS-CoV-2 RNA strains and the expressed viral proteins within the mitochondrial matrix. The functional consequences of interactive heritable and acquired mtDNA polymorphisms on the progression of COVID -19 may then be traced to intrinsic differences in the maternal haplotype associated with bidirectional communication of essential mitochondrial processes with mtDNA gene expression.

Conclusions

This review has presented a critically oriented discussion of potential functional implications of mtDNA heteroplasmy in human health and disease with special reference to its putative role in dysregulation of metabolic and immunological processes associated with rundown of mitochondrial function in COVID-19 neurological disorders. We contend that mtDNA heteroplasmy serves as a dynamic informational reservoir and key regulatory mediator that facilitates bidirectional communication pathways involved in temporal linkage of nuclear and mitochondrial gene expression found in viruses, prokaryotic, and eukaryotic cells. Chronic dysregulation of mitochondrial function leads to the initiation and persistence of diverse pathophysiological states. It may result from the temporal loss of ongoing restorative metabolic and immunological processes inherently linked to the normative cell-specific patterns of mtDNA heteroplasmy.

We have assembled a working mechanistic model proposed for the progression of COVID-19, which may involve extracellular and blood-borne transport and uptake of intact mitochondria with resident biologically active species of SARS-CoV-2 RNA and viral proteins that function as infectious units or complex inflammatory vehicles. The resulting multiple effects of mtDNA heteroplasmy to regulate normative mitochondrial dynamics may alter the patterns of SARS-CoV-2 RNA species and expressed viral proteins within the mitochondrial matrix. This hypothesis supports the need for empirical studies to test the relationships between the mtDNA heteroplasmy and the immune response to SARS-CoV-2 infection, including in the CNS and in the GI tract.

Complex mitochondrial functions from primordial prokaryotic endosymbiont processes may be dependent on mtDNA heteroplasmy as an informational reservoir. Further understanding of mtDNA copy number heteroplasmic distribution may represent a dynamic algorithm for targeting bidirectional communication processes between nuclear and mitochondrial-genomes in SARS-CoV-2 infection. Increased understanding of mtDNA heteroplasmy in the ongoing COVID-19 pandemic involves successful innate and adaptive immunological defense mechanisms dependent on the conformational matching of shared genetic vocabularies. Heteroplasmy of mtDNA provides a molecular substrate for creating novel enhanced immunological selection relevant to SARS-CoV-2 and future emerging pathogenic viruses.