Splice-Break: exploiting an RNA-seq splice junction algorithm to discover mitochondrial DNA deletion breakpoints and analyses of psychiatric disorders

High number of mitochondrial DNA alterations in postmortem brain tissue of patients with schizophrenia compared to healthy controls

Previous studies have shown mitochondrial dysfunction in schizophrenia (SZ) patients, which may be caused by mitochondrial DNA (mtDNA) alterations. However, there are few studies in SZ that have analyzed mtDNA in brain samples by next-generation sequencing (NGS). To address this gap, we used mtDNA-targeted NGS and qPCR to characterize mtDNA alterations in brain samples from patients with SZ (n = 40) and healthy controls (HC) (n = 40). 35 % of SZ patients showed mtDNA alterations, a significantly higher prevalence compared to 10 % of HC. Specifically, SZ patients had a significantly higher frequency of deletions (35 vs. 5 in HC), with a mean number of deletions of 3.8 in SZ vs. 1.0 in HC. Likely pathogenic missense variants were also significantly more frequent in patients with SZ than in HC (10 vs. three HC), encompassing 14 variants in patients and three in HC. The pathogenic tRNA variant m.3243A>G was identified in one SZ patient with a high heteroplasmy level of 32.2 %. While no significant differences in mtDNA copy number (mtDNA-CN) were observed between SZ and HC, antipsychotic users had significantly higher mtDNA-CN than non-users. These findings suggest a potential role for mtDNA alterations in the pathophysiology of SZ that require further validation and functional studies.

Common mitochondrial deletions in RNA-Seq: evaluation of bulk, single-cell, and spatial transcriptomic datasets

Common mitochondrial DNA (mtDNA) deletions are large structural variants in the mitochondrial genome that accumulate in metabolically active tissues with age and have been investigated in various diseases. We applied the Splice-Break2 pipeline (designed for high-throughput quantification of mtDNA deletions) to human RNA-Seq datasets and describe the methodological considerations for evaluating common deletions in bulk, single-cell, and spatial transcriptomics datasets. A robust evaluation of 1570 samples from 14 RNA-Seq studies showed: (i) the abundance of some common deletions detected in PCR-amplified mtDNA correlates with levels observed in RNA-Seq data; (ii) RNA-Seq library preparation method has a strong effect on deletion detection; (iii) deletions had a significant, positive correlation with age in brain and muscle; (iv) deletions were enriched in cortical grey matter, specifically in layers 3 and 5; and (v) brain regions with dopaminergic neurons (i.e., substantia nigra, ventral tegmental area, and caudate nucleus) had remarkable enrichment of common mtDNA deletions.

RNA G-quadruplex folding is a multi-pathway process driven by conformational entropy

The kinetics of folding is crucial for the function of many regulatory RNAs including RNA G-quadruplexes (rG4s). Here, we characterize the folding pathways of a G-quadruplex from the telomeric repeat-containing RNA by combining all-atom molecular dynamics and coarse-grained simulations with circular dichroism experiments. The quadruplex fold is stabilized by cations and thus, the ion atmosphere forming a double layer surrounding the highly charged quadruplex guides the folding process. To capture the ionic double layer in implicit solvent coarse-grained simulations correctly, we develop a matching procedure based on all-atom simulations in explicit water. The procedure yields quantitative agreement between simulations and experiments as judged by the populations of folded and unfolded states at different salt concentrations and temperatures. Subsequently, we show that coarse-grained simulations with a resolution of three interaction sites per nucleotide are well suited to resolve the folding pathways and their intermediate states. The results reveal that the folding progresses from unpaired chain via hairpin, triplex and double-hairpin constellations to the final folded structure. The two- and three-strand intermediates are stabilized by transient Hoogsteen interactions. Each pathway passes through two on-pathway intermediates. We hypothesize that conformational entropy is a hallmark of rG4 folding. Conformational entropy leads to the observed branched multi-pathway folding process for TERRA25. We corroborate this hypothesis by presenting the free energy landscapes and folding pathways of four rG4 systems with varying loop length.

Enhanced mitochondrial genome analysis: bioinformatic and long-read sequencing advances and their diagnostic implications

INTRODUCTION

Primary mitochondrial diseases (PMDs) comprise a large and heterogeneous group of genetic diseases that result from pathogenic variants in either nuclear DNA (nDNA) or mitochondrial DNA (mtDNA). Widespread adoption of next-generation sequencing (NGS) has improved the efficiency and accuracy of mtDNA diagnoses; however, several challenges remain.

AREAS COVERED

In this review, we briefly summarize the current state of the art in molecular diagnostics for mtDNA and consider the implications of improved whole genome sequencing (WGS), bioinformatic techniques, and the adoption of long-read sequencing, for PMD diagnostics.

EXPERT OPINION

We anticipate that the application of PCR-free WGS from blood DNA will increase in diagnostic laboratories, while for adults with myopathic presentations, WGS from muscle DNA may become more widespread. Improved bioinformatic strategies will enhance WGS data interrogation, with more accurate delineation of mtDNA and NUMTs (nuclear mitochondrial DNA segments) in WGS data, superior coverage uniformity, indirect measurement of mtDNA copy number, and more accurate interpretation of heteroplasmic large-scale rearrangements (LSRs). Separately, the adoption of diagnostic long-read sequencing could offer greater resolution of complex LSRs and the opportunity to phase heteroplasmic variants.

Single-cell multi-omics of mitochondrial DNA disorders reveals dynamics of purifying selection across human immune cells

Pathogenic mutations in mitochondrial DNA (mtDNA) compromise cellular metabolism, contributing to cellular heterogeneity and disease. Diverse mutations are associated with diverse clinical phenotypes, suggesting distinct organ- and cell-type-specific metabolic vulnerabilities. Here we establish a multi-omics approach to quantify deletions in mtDNA alongside cell state features in single cells derived from six patients across the phenotypic spectrum of single large-scale mtDNA deletions (SLSMDs). By profiling 206,663 cells, we reveal the dynamics of pathogenic mtDNA deletion heteroplasmy consistent with purifying selection and distinct metabolic vulnerabilities across T-cell states in vivo and validate these observations in vitro. By extending analyses to hematopoietic and erythroid progenitors, we reveal mtDNA dynamics and cell-type-specific gene regulatory adaptations, demonstrating the context-dependence of perturbing mitochondrial genomic integrity. Collectively, we report pathogenic mtDNA heteroplasmy dynamics of individual blood and immune cells across lineages, demonstrating the power of single-cell multi-omics for revealing fundamental properties of mitochondrial genetics.

Mitochondria-wide association study observed significant interactions of mitochondrial respiratory and the inflammatory in the development of anxiety and depression

The aim of this study was to investigate the possible interaction of mitochondrial dysfunction and inflammatory cytokines in the risk of anxiety and depression. We utilized the UK Biobank for the sample of this study. A mitochondria-wide association(MiWAS) and interaction analysis was performed to investigate the interaction effects of mitochondrial DNA (mtDNA)×C-reactive protein (CRP) on the risks of self-reported anxiety (N = 72,476), general anxiety disorder (GAD-7) scores (N = 80,853), self-reported depression (N = 80,778), Patient Health Questionnaire (PHQ-9) scores (N = 80,520) in total samples, females and males, respectively, adjusting for sex, age, Townsend deprivation index (TDI), education score, alcohol intake, smoking and 10 principal components. In all, 25 mtSNPs and 10 mtSNPs showed significant level of association with self-reported anxiety and GAD-7 score respectively. A total of seven significant mtDNA × CRP interactions were found for anxiety, such as m.3915G>A(MT-ND1) for self-reported anxiety in total subjects (P = 6.59 × 10^-3), m.4561T>C(MT-ND2) (P = 3.04 × 10^-3) for GAD-7 score in total subjects. For depression, MiWAS identified 17 significant mtSNPs for self-reported depression and 14 significant mtSNPs for PHQ-9 scores. 17 significant mtDNA associations (2 for self-reported depression and 15 for PHQ-9 score) was identified, such as m.14869G>A(MT-CYB; P = 2.22 × 10^-3) associated with self-reported depression and m.4561T>C (MT-ND2; P value = 3.02 × 10^-8) associated with PHQ-9 score in all subjects. In addition, 5 common mtDNA shared with anxiety and depression were found in MiWAS, and 4 common mtDNA variants were detected to interact with CRP for anxiety and depression, such as m.9899T>C(MT-CO3). Our study suggests the important interaction effects of mitochondrial function and CRP on the risks of anxiety and depression.

Assessing the joint effects of mitochondrial function and human behavior on the risks of anxiety and depression

BACKGROUND

Psychiatric disorders have great health hazards and the exact pathogeny remains elusive now. We aim to explore the potential interaction effects of mitochondrial function and human behavior on the risks of anxiety and depression.

METHODS

The genome-wide association study (GWAS) data of mitochondrial function (N = 383,476-982,072) were obtained from published studies. Individual level genotype and phenotype data of anxiety, depression and behavioral factors (including drinking, smoking and physical activity) were all from the UK Biobank (N = 84,805-85,164). We first calculated the polygenic risk scores (PRS) of mitochondrial function as the instrumental variables, and then constructed linear regression analyses to systematically explore the potential interaction effects of mitochondrial function and human behavior on anxiety and depression.

RESULTS

In total samples, we observed mitochondrial heteroplasmy (MtHz) vs. Drinking (P_GAD-7 = 6.49 × 10^-3; P_PHQ-9 = 1.89 × 10^-3) was positively associated with both anxiety and depression. In males, MtHz vs. Drinking (P_Male = 3.46 × 10^-5) was positively correlated with depression. In females, blood mitochondrial DNA copy number (mtDNA-CN) vs. Drinking (P_Female = 8.63 × 10^-3) was negatively related to anxiety. Furthermore, we identified additional 6 suggestive interaction effects (P < 0.05) for anxiety and depression.

LIMITATIONS

Considering all subjects were from UK Biobank, it should be careful to extrapolate our findings to other populations with different genetic background.

CONCLUSIONS

Our results suggest the significant impacts of mitochondrial function and human behavior interactions on the development of anxiety and depression, providing new clues for clarifying the pathogenesis of anxiety and depression.

Large Common Mitochondrial DNA Deletions Are Associated with a Mitochondrial SNP T14798C Near the 3' Breakpoints

Introduction

Large somatic deletions of mitochondrial DNA (mtDNA) accumulate with aging in metabolically active tissues such as the brain. We have cataloged the breakpoints and frequencies of large mtDNA deletions in the human brain.

Methods

We quantified 112 high-frequency mtDNA somatic deletions across four human brain regions with the Splice-Break2 pipeline. In addition, we utilized PLINK/Seq to test the association of mitochondrial genotypes with the abundance of these high-frequency mtDNA deletions. A conservative p value threshold of 5E-08 was used to find the significant loci.

Results

One mtDNA SNP (T14798C) was significantly associated with mtDNA deletions in two brain regions, the dorsolateral prefrontal cortex (DLPFC) and the superior temporal gyrus. Since the DLPFC showed the most robust association between T14798C and two deletion breakpoints (7816-14807 and 5462-14807), this association was tested in the DLPFC of a replication sample and validated the first results. Incorporating the C allele at 14,798 bp increased the perfect/imperfect length of the repeat at the 3' breakpoint of the two associated deletions.

Conclusion

This is the first study to identify the association of mtDNA SNP with large mtDNA deletions in the human brain. The T14798C allele located in the MT-CYB gene is a common polymorphism that occurs in several mitochondrial haplogroups. We hypothesize that the T14798C association with two deletions occurs by extending the repeat length around the 3' deletion breakpoints. This simple mechanism suggests that mtDNA SNPs can affect the mitochondrial genome structure, especially in brain where high levels of reactive oxygen species lead to deletion accumulation with aging.

Mitochondria DNA copy number, mitochondria DNA total somatic deletions, Complex I activity, synapse number, and synaptic mitochondria number are altered in schizophrenia and bipolar disorder

Mitochondrial dysfunction is a neurobiological phenomenon implicated in the pathophysiology of schizophrenia and bipolar disorder that can synergistically affect synaptic neurotransmission. We hypothesized that schizophrenia and bipolar disorder share molecular alterations at the mitochondrial and synaptic levels. Mitochondria DNA (mtDNA) copy number (CN), mtDNA common deletion (CD), mtDNA total deletion, complex I activity, synapse number, and synaptic mitochondria number were studied in the postmortem human dorsolateral prefrontal cortex (DLPFC), superior temporal gyrus (STG), primary visual cortex (V1), and nucleus accumbens (NAc) of controls (CON), and subjects with schizophrenia (SZ), and bipolar disorder (BD). The results showed (i) the mtDNA CN is significantly higher in DLPFC of both SZ and BD, decreased in the STG of BD, and unaltered in V1 and NAc of both SZ and BD; (ii) the mtDNA CD is significantly higher in DLPFC of BD while unaltered in STG, V1, and NAc of both SZ and BD; (iii) The total deletion burden is significantly higher in DLPFC in both SZ and BD while unaltered in STG, V1, and NAc of SZ and BD; (iv) Complex I activity is significantly lower in DLPFC of both SZ and BD, which is driven by the presence of medications, with no alteration in STG, V1, and NAc. In addition, complex I protein concentration, by ELISA, was decreased across three cortical regions of SZ and BD subjects; (v) The number of synapses is decreased in DLPFC of both SZ and BD, while the synaptic mitochondria number was significantly lower in female SZ and female BD compared to female controls. Overall, these findings will pave the way to understand better the pathophysiology of schizophrenia and bipolar disorder for therapeutic interventions.

Advancements in Genomic and Behavioral Neuroscience Analysis for the Study of Normal and Pathological Brain Function

Psychiatric and neurological disorders are influenced by an undetermined number of genes and molecular pathways that may differ among afflicted individuals. Functionally testing and characterizing biological systems is essential to discovering the interrelationship among candidate genes and understanding the neurobiology of behavior. Recent advancements in genetic, genomic, and behavioral approaches are revolutionizing modern neuroscience. Although these tools are often used separately for independent experiments, combining these areas of research will provide a viable avenue for multidimensional studies on the brain. Herein we will briefly review some of the available tools that have been developed for characterizing novel cellular and animal models of human disease. A major challenge will be openly sharing resources and datasets to effectively integrate seemingly disparate types of information and how these systems impact human disorders. However, as these emerging technologies continue to be developed and adopted by the scientific community, they will bring about unprecedented opportunities in our understanding of molecular neuroscience and behavior.

Comprehensive summary of mitochondrial DNA alterations in the postmortem human brain: A systematic review

BACKGROUND

Mitochondrial DNA (mtDNA) encodes 37 genes necessary for synthesizing 13 essential subunits of the oxidative phosphorylation system. mtDNA alterations are known to cause mitochondrial disease (MitD), a clinically heterogeneous group of disorders that often present with neuropsychiatric symptoms. Understanding the nature and frequency of mtDNA alterations in health and disease could be a cornerstone in disentangling the relationship between biochemical findings and clinical symptoms of brain disorders. This systematic review aimed to summarize the mtDNA alterations in human brain tissue reported to date that have implications for further research on the pathophysiological significance of mtDNA alterations in brain functioning.

METHODS

We searched the PubMed and Embase databases using distinct terms related to postmortem human brain and mtDNA up to June 10, 2021. Reports were eligible if they were empirical studies analysing mtDNA in postmortem human brains.

FINDINGS

A total of 158 of 637 studies fulfilled the inclusion criteria and were clustered into the following groups: MitD (48 entries), neurological diseases (NeuD, 55 entries), psychiatric diseases (PsyD, 15 entries), a miscellaneous group with controls and other clinical diseases (5 entries), ageing (20 entries), and technical issues (5 entries). Ten entries were ascribed to more than one group. Pathogenic single nucleotide variants (pSNVs), both homo- or heteroplasmic variants, have been widely reported in MitD, with heteroplasmy levels varying among brain regions; however, pSNVs are rarer in NeuD, PsyD and ageing. A lower mtDNA copy number (CN) in disease was described in most, but not all, of the identified studies. mtDNA deletions were identified in individuals in the four clinical categories and ageing. Notably, brain samples showed significantly more mtDNA deletions and at higher heteroplasmy percentages than blood samples, and several of the deletions present in the brain were not detected in the blood. Finally, mtDNA heteroplasmy, mtDNA CN and the deletion levels varied depending on the brain region studied.

INTERPRETATION

mtDNA alterations are well known to affect human tissues, including the brain. In general, we found that studies of MitD, NeuD, PsyD, and ageing were highly variable in terms of the type of disease or ageing process investigated, number of screened individuals, studied brain regions and technology used. In NeuD and PsyD, no particular type of mtDNA alteration could be unequivocally assigned to any specific disease or diagnostic group. However, the presence of mtDNA deletions and mtDNA CN variation imply a role for mtDNA in NeuD and PsyD. Heteroplasmy levels and threshold effects, affected brain regions, and mitotic segregation patterns of mtDNA alterations may be involved in the complex inheritance of NeuD and PsyD and in the ageing process. Therefore, more information is needed regarding the type of mtDNA alteration, the affected brain regions, the heteroplasmy levels, and their relationship with clinical phenotypes and the ageing process.

FUNDING

Hospital Universitari Institut Pere Mata; Institut d'Investigació Sanitària Pere Virgili; Instituto de Salud Carlos III, Ministerio de Ciencia e Innovación (PI18/00514).

Rapid structure-function insights via hairpin-centric analysis of big RNA structure probing datasets

The functions of RNA are often tied to its structure, hence analyzing structure is of significant interest when studying cellular processes. Recently, large-scale structure probing (SP) studies have enabled assessment of global structure-function relationships via standard data summarizations or local folding. Here, we approach structure quantification from a hairpin-centric perspective where putative hairpins are identified in SP datasets and used as a means to capture local structural effects. This has the advantage of rapid processing of big (e.g. transcriptome-wide) data as RNA folding is circumvented, yet it captures more information than simple data summarizations. We reformulate a statistical learning algorithm we previously developed to significantly improve precision of hairpin detection, then introduce a novel nucleotide-wise measure, termed the hairpin-derived structure level (HDSL), which captures local structuredness by accounting for the presence of likely hairpin elements. Applying HDSL to data from recent studies recapitulates, strengthens and expands on their findings which were obtained by more comprehensive folding algorithms, yet our analyses are orders of magnitude faster. These results demonstrate that hairpin detection is a promising avenue for global and rapid structure-function analysis, furthering our understanding of RNA biology and the principal features which drive biological insights from SP data.

Triplex and other DNA motifs show motif-specific associations with mitochondrial DNA deletions and species lifespan

The "theory of resistant biomolecules" posits that long-lived species show resistance to molecular damage at the level of their biomolecules. Here, we test this hypothesis in the context of mitochondrial DNA (mtDNA) as it implies that predicted mutagenic DNA motifs should be inversely correlated with species maximum lifespan (MLS). First, we confirmed that guanine-quadruplex and direct repeat (DR) motifs are mutagenic, as they associate with mtDNA deletions in the human major arc of mtDNA, while also adding mirror repeat (MR) and intramolecular triplex motifs to a growing list of potentially mutagenic features. What is more, triplex motifs showed disease-specific associations with deletions and an apparent interaction with guanine-quadruplex motifs. Surprisingly, even though DR, MR and guanine-quadruplex motifs were associated with mtDNA deletions, their correlation with MLS was explained by the biased base composition of mtDNA. Only triplex motifs negatively correlated with MLS even after adjusting for body mass, phylogeny, mtDNA base composition and effective number of codons. Taken together, our work highlights the importance of base composition for the comparative biogerontology of mtDNA and suggests that future research on mitochondrial triplex motifs is warranted.

Accurate mapping of mitochondrial DNA deletions and duplications using deep sequencing

Deletions and duplications in mitochondrial DNA (mtDNA) cause mitochondrial disease and accumulate in conditions such as cancer and age-related disorders, but validated high-throughput methodology that can readily detect and discriminate between these two types of events is lacking. Here we establish a computational method, MitoSAlt, for accurate identification, quantification and visualization of mtDNA deletions and duplications from genomic sequencing data. Our method was tested on simulated sequencing reads and human patient samples with single deletions and duplications to verify its accuracy. Application to mouse models of mtDNA maintenance disease demonstrated the ability to detect deletions and duplications even at low levels of heteroplasmy.

Unbiased PCR-free spatio-temporal mapping of the mtDNA mutation spectrum reveals brain region-specific responses to replication instability

Background

The accumulation of mtDNA mutations in different tissues from various mouse models has been widely studied especially in the context of mtDNA mutation-driven ageing but has been confounded by the inherent limitations of the most widely used approaches. By implementing a method to sequence mtDNA without PCR amplification prior to library preparation, we map the full unbiased mtDNA mutation spectrum across six distinct brain regions from mice.

Results

We demonstrate that ageing-induced levels of mtDNA mutations (single nucleotide variants and deletions) reach stable levels at 50 weeks of age but can be further elevated specifically in the cortex, nucleus accumbens (NAc), and paraventricular thalamic nucleus (PVT) by expression of a proof-reading-deficient mitochondrial DNA polymerase, Polg^D181A. The increase in single nucleotide variants increases the fraction of shared SNVs as well as their frequency, while characteristics of deletions remain largely unaffected. In addition, Polg^D181A also induces an ageing-dependent accumulation of non-coding control-region multimers in NAc and PVT, a feature that appears almost non-existent in wild-type mice.

Conclusions

Our data provide a novel view of the spatio-temporal accumulation of mtDNA mutations using very limited tissue input. The differential response of brain regions to a state of replication instability provides insight into a possible heterogenic mitochondrial landscape across the brain that may be involved in the ageing phenotype and mitochondria-associated disorders.

Ultrasensitive deletion detection links mitochondrial DNA replication, disease, and aging

BACKGROUND

Acquired human mitochondrial genome (mtDNA) deletions are symptoms and drivers of focal mitochondrial respiratory deficiency, a pathological hallmark of aging and late-onset mitochondrial disease.

RESULTS

To decipher connections between these processes, we create LostArc, an ultrasensitive method for quantifying deletions in circular mtDNA molecules. LostArc reveals 35 million deletions (~ 470,000 unique spans) in skeletal muscle from 22 individuals with and 19 individuals without pathogenic variants in POLG. This nuclear gene encodes the catalytic subunit of replicative mitochondrial DNA polymerase γ. Ablation, the deleted mtDNA fraction, suffices to explain skeletal muscle phenotypes of aging and POLG-derived disease. Unsupervised bioinformatic analyses reveal distinct age- and disease-correlated deletion patterns.

CONCLUSIONS

These patterns implicate replication by DNA polymerase γ as the deletion driver and suggest little purifying selection against mtDNA deletions by mitophagy in postmitotic muscle fibers. Observed deletion patterns are best modeled as mtDNA deletions initiated by replication fork stalling during strand displacement mtDNA synthesis.

Psychiatric drugs impact mitochondrial function in brain and other tissues

Mitochondria have been linked to the etiology of schizophrenia (SZ). However, studies of mitochondria in SZ might be confounded by the effects of pharmacological treatment with antipsychotic drugs (APDs) and other common medications. This review summarizes findings on relevant mitochondria mechanisms underlying SZ, and the potential impact of psychoactive drugs including primarily APDs, but also antidepressants and anxiolytics. The summarized data suggest that APDs impair mitochondria function by decreasing Complex I activity and ATP production and dissipation of the mitochondria membrane potential. At the same time, in the brains of patients with SZ, antipsychotic drug treatment normalizes gene expression modules enriched in mitochondrial genes that are decreased in SZ. This indicates that APDs may have both positive and negative effects on mitochondria. The available evidence suggests three conclusions i) alterations in mitochondria functions in SZ exist prior to APD treatment, ii) mitochondria alterations in SZ can be reversed by APD treatment, and iii) APDs directly cause impairment of mitochondria function. Overall, the mechanisms of action of psychiatric drugs on mitochondria are both direct and indirect; we conclude the effects of APDs on mitochondria may contribute to both their therapeutic and metabolic side effects. These studies support the hypothesis that neuronal mitochondria are an etiological factor in SZ. Moreover, APDs and other drugs must be considered in the evaluation of this pathophysiological role of mitochondria in SZ. Considering these effects, pharmacological actions on mitochondria may be a worthwhile target for further APD development.

Novel Complex Interactions between Mitochondrial and Nuclear DNA in Schizophrenia and Bipolar Disorder

Mitochondrial dysfunction has been associated with schizophrenia (SZ) and bipolar disorder (BD). This review examines recent publications and novel associations between mitochondrial genes and SZ and BD. Associations of nuclear-encoded mitochondrial variants with SZ were found using gene- and pathway-based approaches. Two control region mitochondrial DNA (mtDNA) SNPs, T16519C and T195C, both showed an association with SZ and BD. A review of 4 studies of A15218G located in the cytochrome B oxidase gene (CYTB, SZ = 11,311, control = 35,735) shows a moderate association with SZ (p = 2.15E-03). Another mtDNA allele A12308G was nominally associated with psychosis in BD type I subjects and SZ. The first published study testing the epistatic interaction between nuclear-encoded and mitochondria-encoded genes demonstrated evidence for potential interactions between mtDNA and the nuclear genome for BD. A similar analysis for the risk of SZ revealed significant joint effects (34 nuclear-mitochondria SNP pairs with joint effect p ≤ 5E-07) and significant enrichment of projection neurons. The mitochondria-encoded gene CYTB was found in both the epistatic interactions for SZ and BD and the single SNP association of SZ. Future efforts considering population stratification and polygenic risk scores will test the role of mitochondrial variants in psychiatric disorders.