Deep sequencing of the mitochondrial genome reveals common heteroplasmic sites in NADH dehydrogenase genes

Association analysis of mitochondrial DNA heteroplasmic variants: Methods and application

We rigorously assessed a comprehensive association testing framework for heteroplasmy, employing both simulated and real-world data. This framework employed a variant allele fraction (VAF) threshold and harnessed multiple gene-based tests for robust identification and association testing of heteroplasmy. Our simulation studies demonstrated that gene-based tests maintained an appropriate type I error rate at α = 0.001. Notably, when 5 % or more heteroplasmic variants within a target region were linked to an outcome, burden-extension tests (including the adaptive burden test, variable threshold burden test, and z-score weighting burden test) outperformed the sequence kernel association test (SKAT) and the original burden test. Applying this framework, we conducted association analyses on whole-blood derived heteroplasmy in 17,507 individuals of African and European ancestries (31 % of African Ancestry, mean age of 62, with 58 % women) with whole genome sequencing data. We performed both cohort- and ancestry-specific association analyses, followed by meta-analysis on both pooled samples and within each ancestry group. Our results suggest that mtDNA-encoded genes/regions are likely to exhibit varying rates in somatic aging, with the notably strong associations observed between heteroplasmy in the RNR1 and RNR2 genes (p < 0.001) and advance aging by the Original Burden test. In contrast, SKAT identified significant associations (p < 0.001) between diabetes and the aggregated effects of heteroplasmy in several protein-coding genes. Further research is warranted to validate these findings. In summary, our proposed statistical framework represents a valuable tool for facilitating association testing of heteroplasmy with disease traits in large human populations.

MitoH3: Mitochondrial Haplogroup and Homoplasmic/Heteroplasmic Variant Calling Pipeline for Alzheimer's Disease Sequencing Project

Background

Mitochondrial DNA (mtDNA) is a double-stranded circular DNA and has multiple copies in each cell. Excess heteroplasmy, the coexistence of distinct variants in copies of mtDNA within a cell, may lead to mitochondrial impairments. Accurate determination of heteroplasmy in whole-genome sequencing (WGS) data has posed a significant challenge because mitochondria carrying heteroplasmic variants cannot be distinguished during library preparation. Moreover, sequencing errors, contamination, and nuclear mtDNA segments can reduce the accuracy of heteroplasmic variant calling.

Objective

To efficiently and accurately call mtDNA homoplasmic and heteroplasmic variants from the large-scale WGS data generated from the Alzheimer's Disease Sequencing Project (ADSP), and test their association with Alzheimer's disease (AD).

Methods

In this study, we present MitoH3-a comprehensive computational pipeline for calling mtDNA homoplasmic and heteroplasmic variants and inferring haplogroups in the ADSP WGS data. We first applied MitoH3 to 45 technical replicates from 6 subjects to define a threshold for detecting heteroplasmic variants. Then using the threshold of 5% ≤variant allele fraction≤95%, we further applied MitoH3 to call heteroplasmic variants from a total of 16,113 DNA samples with 6,742 samples from cognitively normal controls and 6,183 from AD cases.

Results

This pipeline is available through the Singularity container engine. For 4,311 heteroplasmic variants identified from 16,113 samples, no significant variant count difference was observed between AD cases and controls.

Conclusions

Our streamlined pipeline, MitoH3, enables computationally efficient and accurate analysis of a large number of samples.

A naturally occurring variant of SHLP2 is a protective factor in Parkinson's disease

Mitochondrial DNA single nucleotide polymorphisms (mtSNPs) have been associated with a reduced risk of developing Parkinson's disease (PD), yet the underlying mechanisms remain elusive. In this study, we investigate the functional role of a PD-associated mtSNP that impacts the mitochondrial-derived peptide (MDP) Small Humanin-like Peptide 2 (SHLP2). We identify m.2158 T > C, a mtSNP associated with reduced PD risk, within the small open reading frame encoding SHLP2. This mtSNP results in an alternative form of SHLP2 (lysine 4 replaced with arginine; K4R). Using targeted mass spectrometry, we detect specific tryptic fragments of SHLP2 in neuronal cells and demonstrate its binding to mitochondrial complex 1. Notably, we observe that the K4R variant, associated with reduced PD risk, exhibits increased stability compared to WT SHLP2. Additionally, both WT and K4R SHLP2 show enhanced protection against mitochondrial dysfunction in in vitro experiments and confer protection against a PD-inducing toxin, a mitochondrial complex 1 inhibitor, in a mouse model. This study sheds light on the functional consequences of the m.2158 T > C mtSNP on SHLP2 and provides insights into the potential mechanisms by which this mtSNP may reduce the risk of PD.

Association analysis of mitochondrial DNA heteroplasmic variants: methods and application

We rigorously assessed a comprehensive association testing framework for heteroplasmy, employing both simulated and real-world data. This framework employed a variant allele fraction (VAF) threshold and harnessed multiple gene-based tests for robust identification and association testing of heteroplasmy. Our simulation studies demonstrated that gene-based tests maintained an appropriate type I error rate at α=0.001. Notably, when 5% or more heteroplasmic variants within a target region were linked to an outcome, burden-extension tests (including the adaptive burden test, variable threshold burden test, and z-score weighting burden test) outperformed the sequence kernel association test (SKAT) and the original burden test. Applying this framework, we conducted association analyses on whole-blood derived heteroplasmy in 17,507 individuals of African and European ancestries (31% of African Ancestry, mean age of 62, with 58% women) with whole genome sequencing data. We performed both cohort- and ancestry-specific association analyses, followed by meta-analysis on both pooled samples and within each ancestry group. Our results suggest that mtDNA-encoded genes/regions are likely to exhibit varying rates in somatic aging, with the notably strong associations observed between heteroplasmy in the RNR1 and RNR2 genes (p<0.001) and advance aging by the Original Burden test. In contrast, SKAT identified significant associations (p<0.001) between diabetes and the aggregated effects of heteroplasmy in several protein-coding genes. Further research is warranted to validate these findings. In summary, our proposed statistical framework represents a valuable tool for facilitating association testing of heteroplasmy with disease traits in large human populations.

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.

A Pilot Mitochondrial Genome-Wide Association on Migraine Among Saudi Arabians

Background

Mitochondrial DNA (mtDNA) mutations have been reported in multiple neurological diseases and helped to explain the pathophysiology of these diseases. Similarly, variations in mtDNA might exist in migraine and can explain the effect of low ATP production in the neurons on the initiation of migraine attack. Therefore, in the current study we aim to explore the association of mtDNA mutations on migraine in the Saudi population.

Subjects and Methods

Over 1950 young Saudi female students were screened for migraine, among that a total of 103 satisfied the ICHD-3 criteria. However,  20 migraine cases confirmed in the neurology clinic and gave consent to participate in the study. Another 20 age-matched healthy controls were also recruited. Mitochondrial sequence variations were filtered from exome sequencing using NCBI GenBank Reference Sequence: NC_012920.1 and analysed using MITOMAP. Genes with significant single nucleotide polymorphisms (SNPs) were investigated by the gene functional classification tool DAVID and functional enrichment analysis of protein-protein interaction networks through STRING 11.5 for the most significant associated genes.

Results

Genome wide analysis of the mitochondrial sequence variations between the patients with migraine and control revealed the association of 30 SNPs (p < 0.05) in the mitochondrial genome. The highest significance (p = 0.001033) was observed in a coding SNP (rs1603225278) in the CYTB gene and rs386829281 in the region of origin of replication. Twenty-four significant SNPs were in the coding region of nine (ND5, ND4, COX2, COX1, ND3, CYTB, COX3, ND2 and ND1) genes.

Conclusion

This is the first study to demonstrate the association of mtDNA variations with migraine in the Saudi population. The current findings will help to highlight the significance of mtDNA mutations to migraine pathophysiology and will serve as a reference data for larger national and international studies.

Mitochondrial DNA Is a Vital Driving Force in Ischemia-Reperfusion Injury in Cardiovascular Diseases

According to the latest Global Burden of Disease Study, cardiovascular disease (CVD) is the leading cause of death, and ischemic heart disease and stroke are the cause of death in approximately half of CVD patients. In CVD, mitochondrial dysfunction following ischemia-reperfusion (I/R) injury results in heart failure. The proper functioning of oxidative phosphorylation (OXPHOS) and the mitochondrial life cycle in cardiac mitochondria are closely related to mitochondrial DNA (mtDNA). Following myocardial I/R injury, mitochondria activate multiple repair and clearance mechanisms to repair damaged mtDNA. When these repair mechanisms are insufficient to restore the structure and function of mtDNA, irreversible mtDNA damage occurs, leading to mtDNA mutations. Since mtDNA mutations aggravate OXPHOS dysfunction and affect mitophagy, mtDNA mutation accumulation leads to leakage of mtDNA and proteins outside the mitochondria, inducing an innate immune response, aggravating cardiovascular injury, and leading to the need for external interventions to stop or slow the disease course. On the other hand, mtDNA released into the circulation after cardiac injury can serve as a biomarker for CVD diagnosis and prognosis. This article reviews the pathogenic basis and related research findings of mtDNA oxidative damage and mtDNA leak-triggered innate immune response associated with I/R injury in CVD and summarizes therapeutic options that target mtDNA.

Heteroplasmic mitochondrial DNA variants in cardiovascular diseases

Mitochondria are implicated in the pathogenesis of cardiovascular diseases (CVDs) but the reasons for this are not well understood. Maternally-inherited population variants of mitochondrial DNA (mtDNA) which affect all mtDNA molecules (homoplasmic) are associated with cardiometabolic traits and the risk of developing cardiovascular disease. However, it is not known whether mtDNA mutations only affecting a proportion of mtDNA molecules (heteroplasmic) also play a role. To address this question, we performed a high-depth (~1000-fold) mtDNA sequencing of blood DNA in 1,399 individuals with hypertension (HTN), 1,946 with ischemic heart disease (IHD), 2,146 with ischemic stroke (IS), and 723 healthy controls. We show that the per individual burden of heteroplasmic single nucleotide variants (mtSNVs) increases with age. The age-effect was stronger for low-level heteroplasmies (heteroplasmic fraction, HF, 5-10%), likely reflecting acquired somatic events based on trinucleotide mutational signatures. After correcting for age and other confounders, intermediate heteroplasmies (HF 10-95%) were more common in hypertension, particularly involving non-synonymous variants altering the amino acid sequence of essential respiratory chain proteins. These findings raise the possibility that heteroplasmic mtSNVs play a role in the pathophysiology of hypertension.

Comparison of mitochondrial DNA sequences from whole blood and lymphoblastoid cell lines

Lymphoblastoid cell lines (LCLs) provide an unlimited source of genomic DNA for genetic studies. Here, we compared mtDNA sequence variants, heteroplasmic or homplasmic, between LCL (sequenced by mitoRCA-seq method) and whole blood samples (sequenced through whole genome sequencing approach) of the same 130 participants in the Framingham Heart Study. We applied harmonization of sequence coverages and consistent quality control to mtDNA sequences. We identified 866 variation sites in the 130 LCL samples and 666 sites in the 130 blood samples. More than 94% of the identified homoplasmies were present in both LCL and blood samples while more than 70% of heteroplasmic sites were uniquely present either in LCL or in blood samples. The LCL and whole blood samples carried a similar number of homoplasmic variants (p = 0.45) per sample while the LCL carried a greater number of heteroplasmic variants than whole blood per sample (p < 2.2e-16). Furthermore, the LCL samples tended to accumulate low level heteroplasmies (heteroplasmy level in 3-25%) than their paired blood samples (p = 0.001). These results suggest that cautions should be taken in the interpretation and comparison of findings when different tissues/cell types or different sequencing technologies are applied to obtain mtDNA sequences.

Genetic landscape of human mitochondrial genome using whole genome sequencing

Increasing evidences suggest that mitochondrial dysfunction is implicated in diseases and aging, and whole-genome sequencing (WGS) is the most unbiased method in analyzing the mitochondrial genome (mtDNA). However, the genetic landscape of mtDNA in the Chinese population has not been fully examined. Here, we described the genetic landscape of mtDNA using WGS data from Chinese individuals (n = 3241). We identified 3892 mtDNA variants, of which 3349 (86%) were rare variants. Interestingly, we observed a trend toward extreme heterogeneity of mtDNA variants. Our study observed a distinct purifying selection on mtDNA, which inhibits the accumulation of harmful heteroplasmies at the individual level: (1) mitochondrial dN/dS ratios were much less than 1; (2) the dN/dS ratio of heteroplasmies was higher than homoplasmies; (3) heteroplasmies had more indels and predicted deleterious variants than homoplasmies. Furthermore, we found that haplogroup M (20.27%) and D (20.15%) had the highest frequencies in the Chinese population, followed by B (18.51%) and F (16.45%). The number of variants per individual differed across haplogroup groups, with a higher number of homoplasmies for the M lineage. Meanwhile, mtDNA copy number was negatively correlated with age but positively correlated with the female sex. Finally, we developed an mtDNA variation database of Chinese populations called MTCards (http://genemed.tech/mtcards/) to facilitate the query of mtDNA variants in this study. In summary, these findings contribute to different aspects of understanding mtDNA, providing a better understanding of the genetic basis of mitochondrial-related diseases.

Presence and transmission of mitochondrial heteroplasmic mutations in human populations of European and African ancestry

We investigated the concordance of mitochondrial DNA heteroplasmic mutations (heteroplasmies) in 6745 maternal pairs of European (EA, n = 4718 pairs) and African (AA, n = 2027 pairs) Americans in whole blood. Mother-offspring pairs displayed the highest concordance rate, followed by sibling-sibling and more distantly-related maternal pairs. The allele fractions of concordant heteroplasmies exhibited high correlation (R2 = 0.8) between paired individuals. Discordant heteroplasmies were more likely to be in coding regions, be nonsynonymous or nonsynonymous-deleterious (p < 0.001). The number of deleterious heteroplasmies was significantly correlated with advancing age (20-44, 45-64, and ≥65 years, p-trend = 0.01). One standard deviation increase in heteroplasmic burden (i.e., the number of heteroplasmies carried by an individual) was associated with 0.17 to 0.26 (p < 1e - 23) standard deviation decrease in mtDNA copy number, independent of age. White blood cell count and differential count jointly explained 0.5% to 1.3% (p ≤ 0.001) variance in heteroplasmic burden. A genome-wide association and meta-analysis identified a region at 11p11.12 (top signal rs779031139, p = 2.0e - 18, minor allele frequency = 0.38) associated with the heteroplasmic burden. However, the 11p11.12 region is adjacent to a nuclear mitochondrial DNA (NUMT) corresponding to a 542 bp area of the D-loop. This region was no longer significant after excluding heteroplasmies within the 542 bp from the heteroplasmic burden. The discovery that blood mtDNA heteroplasmies were both inherited and somatic origins and that an increase in heteroplasmic burden was strongly associated with a decrease in average number of mtDNA copy number in blood are important findings to be considered in association studies of mtDNA with disease traits.

Comparison of whole genome sequencing and targeted sequencing for mitochondrial DNA

Mitochondrial dysfunction has emerged to be associated with a broad spectrum of diseases, and there is an increasing demand for accurate detection of mitochondrial DNA (mtDNA) variants. Whole genome sequencing (WGS) has been the dominant sequencing approach to identify genetic variants in recent decades, but most studies focus on variants on the nuclear genome. Whole genome sequencing is also costly and time consuming. Sequencing specifically targeted for mtDNA is commonly used in the diagnostic settings and has lower costs. However, there is a lack of pairwise comparisons between these two sequencing approaches for calling mtDNA variants on a population basis. In this study, we compared WGS and mtDNA-targeted sequencing (targeted-seq) in analyzing mitochondrial DNA from 1499 participants recruited into the Severe Asthma Research Program (SARP). Our study reveals that targeted-sequencing and WGS have comparable capacity to determine genotypes and to call haplogroups and homoplasmies on mtDNA. However, there exists a large variability in calling heteroplasmies, especially for low-frequency heteroplasmies, which indicates that investigators should be cautious about heteroplasmies acquired from different sequencing methods. Further research is highly desired to improve variant detection methods for mitochondrial DNA.

Nuclear genome-wide associations with mitochondrial heteroplasmy

The role of the nuclear genome in maintaining the stability of the mitochondrial genome (mtDNA) is incompletely known. mtDNA sequence variants can exist in a state of heteroplasmy, which denotes the coexistence of organellar genomes with different sequences. Heteroplasmic variants that impair mitochondrial capacity cause disease, and the state of heteroplasmy itself is deleterious. However, mitochondrial heteroplasmy may provide an intermediate state in the emergence of novel mitochondrial haplogroups. We used genome-wide genotyping data from 982,072 European ancestry individuals to evaluate variation in mitochondrial heteroplasmy and to identify the regions of the nuclear genome that affect it. Age, sex, and mitochondrial haplogroup were associated with the extent of heteroplasmy. GWAS identified 20 loci for heteroplasmy that exceeded genome-wide significance. This included a region overlapping mitochondrial transcription factor A (TFAM), which has multiple roles in mtDNA packaging, replication, and transcription. These results show that mitochondrial heteroplasmy has a heritable nuclear component.

Metabolism in Pulmonary Hypertension

Pulmonary arterial hypertension (PAH) is characterized by impaired regulation of pulmonary hemodynamics and vascular growth. Alterations of metabolism and bioenergetics are increasingly recognized as universal hallmarks of PAH, as metabolic abnormalities are identified in lungs and hearts of patients, animal models of the disease, and cells derived from lungs of patients. Mitochondria are the primary organelle critically mediating the complex and integrative metabolic pathways in bioenergetics, biosynthetic pathways, and cell signaling. Here, we review the alterations in metabolic pathways that are linked to the pathologic vascular phenotype of PAH, including abnormalities in glycolysis and glucose oxidation, fatty acid oxidation, glutaminolysis, arginine metabolism, one-carbon metabolism, the reducing and oxidizing cell environment, and the tricarboxylic acid cycle, as well as the effects of PAH-associated nuclear and mitochondrial mutations on metabolism. Understanding of the metabolic mechanisms underlying PAH provides important knowledge for the design of new therapeutics for treatment of patients.

A novel m.11406 T > A mutation in mitochondrial ND4 gene causes MELAS syndrome

Pathogenic point mutations of mitochondrial DNA (mtDNA) are associated with a large number of heterogeneous diseases involving multiple systems with which patients may present with a wide range of clinical phenotypes. In this study, we describe a novel heteroplasmic missense mutation, m.11406 T > A, of the ND4 gene encoding the subunit 4 of mitochondrial complex I in a 32-year-old woman with recurrent epileptic seizure, headache and bilateral hearing loss. Skeletal muscle histochemistry demonstrated that approximately 20% of fibers were cytochrome C oxidase (COX) deficient with increased activity of succinate dehydrogenase (SDH). Further investigations in muscle specimens showed significantly reduced level of ND4 protein. It is interesting that the subunits of complex I (ND1 and NDFUB8) and complex IV(CO1) were also remarkably decreased. These findings indicate that ND1, NDFUB8 and CO1 are more susceptible than other subunits to mutations in the mitochondrial ND4 gene.

Genetic and phenotypic landscape of the mitochondrial genome in the Japanese population

The genetic landscape of mitochondrial DNA (mtDNA) has been elusive. By analyzing mtDNA using the whole genome sequence (WGS) of Japanese individuals (n = 1928), we identified 2023 mtDNA variants and high-resolution haplogroups. Frequency spectra of the haplogroups were population-specific and were heterogeneous among geographic regions within Japan. Application of machine learning methods could finely classify the subjects corresponding to the high-digit mtDNA sub-haplogroups. mtDNA had distinct genetic structures from that of nuclear DNA (nDNA), characterized by no distance-dependent linkage disequilibrium decay, sparse tagging of common variants, and the existence of common haplotypes spanning the entire mtDNA. We did not detect any evidence of mtDNA–nDNA (or mtDNA copy number–nDNA) genotype associations. Together with WGS-based mtDNA variant imputation, we conducted a phenome-wide association study of 147,437 Japanese individuals with 99 clinical phenotypes. We observed pleiotropy of mtDNA genetic risk on the five late-onset human complex traits including creatine kinase (P = 1.7 × 10⁻¹²).

Kenichi Yamamoto et al. report a genetic analysis of mitochondrial DNA (mtDNA) and a phenome-wide association study in Japanese individuals from the BioBank Japan Project. They describe the genetic landscape of the mitochondria and identify pleiotropic mtDNA variants associated with 5 late-onset complex traits.

microRNA-222 Attenuates Mitochondrial Dysfunction During Transmissible Gastroenteritis Virus Infection

Transmissible gastroenteritis virus (TGEV) is a member of Coronaviridae family. Our previous research showed that TGEV infection could induce mitochondrial dysfunction and upregulate miR-222 level. Therefore, we presumed that miR-222 might be implicated in regulating mitochondrial dysfunction induced by TGEV infection. To verify the hypothesis, the effect of miR-222 on mitochondrial dysfunction was tested and we showed that miR-222 attenuated TGEV-induced mitochondrial dysfunction. To investigate the underlying molecular mechanism of miR-222 in TGEV-induced mitochondrial dysfunction, a quantitative proteomic analysis of PK-15 cells that were transfected with miR-222 mimics and infected with TGEV was performed. In total, 4151 proteins were quantified and 100 differentially expressed proteins were obtained (57 upregulated, 43 downregulated), among which thrombospondin-1 (THBS1) and cluster of differentiation 47 (CD47) were downregulated. THBS1 was identified as the target of miR-222. Knockdown of THBS1 and CD47 decreased mitochondrial Ca2+ level and increased mitochondrial membrane potential (MMP) level. Reversely, overexpression of THBS1 and CD47 elevated mitochondrial Ca2+ level and reduced mitochondrial membrane potential (MMP) level. Together, our data establish a significant role of miR-222 in regulating mitochondrial dysfunction in response to TGEV infection.