Comprehensive molecular characterization of mitochondrial genomes in human cancers

Mitochondrial-Encoded Complex I Impairment Induces a Targetable Dependency on Aerobic Fermentation in Hürthle Cell Carcinoma of the Thyroid

A metabolic hallmark of cancer identified by Warburg is the increased consumption of glucose and secretion of lactate, even in the presence of oxygen. Although many tumors exhibit increased glycolytic activity, most forms of cancer rely on mitochondrial respiration for tumor growth. We report here that Hürthle cell carcinoma of the thyroid (HTC) models harboring mitochondrial DNA-encoded defects in complex I of the mitochondrial electron transport chain exhibit impaired respiration and alterations in glucose metabolism. CRISPR-Cas9 pooled screening identified glycolytic enzymes as selectively essential in complex I-mutant HTC cells. We demonstrate in cultured cells and a patient-derived xenograft model that small-molecule inhibitors of lactate dehydrogenase selectively induce an ATP crisis and cell death in HTC. This work demonstrates that complex I loss exposes fermentation as a therapeutic target in HTC and has implications for other tumors bearing mutations that irreversibly damage mitochondrial respiration.

SIGNIFICANCE

HTC is enriched in somatic mtDNA mutations predicted to affect complex I of the electron transport chain (ETC). We demonstrate that these mutations impair respiration and induce a therapeutically tractable reliance on aerobic fermentation for cell survival. This work provides a rationale for targeting fermentation in cancers harboring irreversible genetically encoded ETC defects. See related article by Gopal et al., p. 1904. This article is highlighted in the In This Issue feature, p. 1749.

Effectors Enabling Adaptation to Mitochondrial Complex I Loss in Hürthle Cell Carcinoma

Oncocytic (Hürthle cell) carcinoma of the thyroid (HCC) is genetically characterized by complex I mitochondrial DNA mutations and widespread chromosomal losses. Here, we utilize RNA sequencing and metabolomics to identify candidate molecular effectors activated by these genetic drivers. We find glutathione biosynthesis, amino acid metabolism, mitochondrial unfolded protein response, and lipid peroxide scavenging to be increased in HCC. A CRISPR-Cas9 knockout screen in a new HCC model reveals which pathways are key for fitness, and highlights loss of GPX4, a defense against lipid peroxides and ferroptosis, as a strong liability. Rescuing complex I redox activity with the yeast NADH dehydrogenase (NDI1) in HCC cells diminishes ferroptosis sensitivity, while inhibiting complex I in normal thyroid cells augments ferroptosis induction. Our work demonstrates unmitigated lipid peroxide stress to be an HCC vulnerability that is mechanistically coupled to the genetic loss of mitochondrial complex I activity.

SIGNIFICANCE

HCC harbors abundant mitochondria, mitochondrial DNA mutations, and chromosomal losses. Using a CRISPR-Cas9 screen inspired by transcriptomic and metabolomic profiling, we identify molecular effectors essential for cell fitness. We uncover lipid peroxide stress as a vulnerability coupled to mitochondrial complex I loss in HCC. See related article by Frank et al., p. 1884. This article is highlighted in the In This Issue feature, p. 1749.

Role of mitochondrial alterations in human cancer progression and cancer immunity

Dysregulating cellular metabolism is one of the emerging cancer hallmarks. Mitochondria are essential organelles responsible for numerous physiologic processes, such as energy production, cellular metabolism, apoptosis, and calcium and redox homeostasis. Although the "Warburg effect," in which cancer cells prefer aerobic glycolysis even under normal oxygen circumstances, was proposed a century ago, how mitochondrial dysfunction contributes to cancer progression is still unclear. This review discusses recent progress in the alterations of mitochondrial DNA (mtDNA) and mitochondrial dynamics in cancer malignant progression. Moreover, we integrate the possible regulatory mechanism of mitochondrial dysfunction-mediated mitochondrial retrograde signaling pathways, including mitochondrion-derived molecules (reactive oxygen species, calcium, oncometabolites, and mtDNA) and mitochondrial stress response pathways (mitochondrial unfolded protein response and integrated stress response) in cancer progression and provide the possible therapeutic targets. Furthermore, we discuss recent findings on the role of mitochondria in the immune regulatory function of immune cells and reveal the impact of the tumor microenvironment and metabolism remodeling on cancer immunity. Targeting the mitochondria and metabolism might improve cancer immunotherapy. These findings suggest that targeting mitochondrial retrograde signaling in cancer malignancy and modulating metabolism and mitochondria in cancer immunity might be promising treatment strategies for cancer patients and provide precise and personalized medicine against cancer.

Mitochondrial DNA mutations in Medulloblastoma

To date, several studies on genomic events underlying medulloblastoma (MB) biology have expanded our understanding of this tumour entity and led to its division into four groups-WNT, SHH, group 3 (G3) and group 4 (G4). However, there is little information about the relevance of pathogenic mitochondrial DNA (mtDNA) mutations and their consequences across these. In this report, we describe the case of a female patient with MB and a mitochondriopathy, followed by a study of mtDNA variants in MB groups. After being diagnosed with G4 MB, the index patient was treated in line with the HIT 2000 protocol with no indications of relapse after five years. Long-term side effects of treatment were complemented by additional neurological symptoms and elevated lactate levels ten years later, resulting in suspected mitochondrial disease. This was confirmed by identifying a mutation in the MT-TS1 gene which appeared homoplasmic in patient tissue and heteroplasmic in the patient's mother. Motivated by this case, we explored mtDNA mutations across 444 patients from ICGC and HIT cohorts. While there was no statistically significant enrichment of mutations in one MB group, both cohorts encompassed a small group of patients harbouring potentially deleterious mtDNA variants. The case presented here highlights the possible similarities between sequelae caused by MB treatment and neurological symptoms of mitochondrial dysfunction, which may apply to patients across all MB groups. In the context of the current advances in characterising and interpreting mtDNA aberrations, recognising affected patients could enhance our future knowledge regarding the mutations' impact on carcinogenesis and cancer treatment.

The human placenta exhibits a unique transcriptomic void

The human placenta exhibits a unique genomic architecture with an unexpectedly high mutation burden and many uniquely expressed genes. The aim of this study is to identify transcripts that are uniquely absent or depleted in the placenta. Here, we show that 40 of 46 of the other organs have no selectively depleted transcripts and that, of the remaining six, the liver has the largest number, with 26. In contrast, the term placenta has 762 depleted transcripts. Gene Ontology analysis of this depleted set highlighted multiple pathways reflecting known unique elements of placental physiology. For example, transcripts associated with neuronal function are in the depleted set-as expected given the lack of placental innervation. However, this demonstrated overrepresentation of genes involved in mitochondrial function (p = 5.8 × 10-10), including PGC-1α, the master regulator of mitochondrial biogenesis, and genes involved in polyamine metabolism (p = 2.1 × 10-4).

Mutations in Structural Genes of the Mitochondrial Complex IV May Influence Breast Cancer

Although it has gained more attention in recent years, the relationship between breast cancer (BC) and mitochondrial oxidative phosphorylation (OXPHOS) is still not well understood. Importantly, Complex IV or Cytochrome C Oxidase (COX) of OXPHOS is one of the key players in mitochondrial balance. An in silico investigation of mutations in structural genes of Complex IV was conducted in BC, comprising 2107 samples. Our findings show four variants (rs267606614, rs753969142, rs199476128 and rs267606884) with significant pathogenic potential. Moreover, we highlight nine genes (MT-CO1, MT-CO2, MT-CO3, CO4I2, COX5A, COX5B, COX6A2, COX6C and COX7B2) with a potential impact on BC.

Next-generation sequencing reveals mitogenome diversity in plasma extracellular vesicles from colorectal cancer patients

BACKGROUND

Recent reports have demonstrated that the entire mitochondrial genome can be secreted in extracellular vesicles (EVs), but the biological attributes of this cell-free mitochondrial DNA (mtDNA) remain insufficiently understood. We used next-generation sequencing to compare plasma EV-derived mtDNA to that of whole blood (WB), peripheral blood mononuclear cells (PBMCs), and formalin-fixed paraffin-embedded (FFPE) tumor tissue from eight rectal cancer patients and WB and fresh-frozen (FF) tumor tissue from eight colon cancer patients.

METHODS

Total DNA was isolated before the mtDNA was enriched by PCR with either two primer sets generating two long products or multiple primer sets (for the FFPE tumors), prior to the sequencing. mtDNA diversity was assessed as the total variant number, level of heteroplasmy (mutant mtDNA copies mixed with wild-type copies), variant distribution within the protein-coding genes, and the predicted functional effect of the variants in the different sample types. Differences between groups were compared by paired Student's t-test or ANOVA with Dunnett's multiple comparison tests when comparing matched samples from patients. Mann-Whitney U test was used when comparing differences between the cancer types and patient groups. Pearson correlation analysis was performed.

RESULTS

In both cancer types, EV mtDNA presented twice as many variants and had significantly more low-level heteroplasmy than WB mtDNA. The EV mtDNA variants were clustered in the coding regions, and the proportion of EV mtDNA variants that were missense mutations (i.e., estimated to moderately affect the mitochondrial protein function) was significantly higher than in WB and tumor tissues. Nonsense mutations (i.e., estimated to highly affect the mitochondrial protein function) were only observed in the tumor tissues and EVs.

CONCLUSION

Taken together, plasma EV mtDNA in CRC patients exhibits a high degree of diversity.

TRIAL REGISTRATION

ClinicalTrials.gov: NCT01816607 . Registered 22 March 2013.

TERT expression is associated with metastasis from thin primaries, exhausted CD4+ T cells in melanoma and with DNA repair across cancer entities

Telomerase reverse transcriptase (TERT) promoter mutations occur frequently in cancer, have been associated with increased TERT expression and cell proliferation, and could potentially influence therapeutic regimens for melanoma. As the role of TERT expression in malignant melanoma and the non-canonical functions of TERT remain understudied, we aimed to extend the current knowledge on the impact of TERT promoter mutations and expression alterations in tumor progression by analyzing several highly annotated melanoma cohorts. Using multivariate models, we found no consistent association for TERT promoter mutations or TERT expression with the survival rate in melanoma cohorts under immune checkpoint inhibition. However, the presence of CD4+ T cells increased with TERT expression and correlated with the expression of exhaustion markers. While the frequency of promoter mutations did not change with Breslow thickness, TERT expression was increased in metastases arising from thinner primaries. As single-cell RNA-sequencing (RNA-seq) showed that TERT expression was associated with genes involved in cell migration and dynamics of the extracellular matrix, this suggests a role of TERT during invasion and metastasis. Co-regulated genes found in several bulk tumors and single-cell RNA-seq cohorts also indicated non-canonical functions of TERT related to mitochondrial DNA stability and nuclear DNA repair. This pattern was also evident in glioblastoma and across other entities. Hence, our study adds to the role of TERT expression in cancer metastasis and potentially also immune resistance.

Mitochondrial DNA is a target of HBV integration

Hepatitis B virus (HBV) may integrate into the genome of infected cells and contribute to hepatocarcinogenesis. However, the role of HBV integration in hepatocellular carcinoma (HCC) development remains unclear. In this study, we apply a high-throughput HBV integration sequencing approach that allows sensitive identification of HBV integration sites and enumeration of integration clones. We identify 3339 HBV integration sites in paired tumour and non-tumour tissue samples from 7 patients with HCC. We detect 2107 clonally expanded integrations (1817 in tumour and 290 in non-tumour tissues), and a significant enrichment of clonal HBV integrations in mitochondrial DNA (mtDNA) preferentially occurring in the oxidative phosphorylation genes (OXPHOS) and D-loop region. We also find that HBV RNA sequences are imported into the mitochondria of hepatoma cells with the involvement of polynucleotide phosphorylase (PNPASE), and that HBV RNA might have a role in the process of HBV integration into mtDNA. Our results suggest a potential mechanism by which HBV integration may contribute to HCC development.

The contribution of mitochondrial DNA alterations to aging, cancer, and neurodegeneration

Mitochondrial DNA (mtDNA) is as a double-stranded molecule existing in hundreds to thousands copies in cells depending on cell metabolism and exposure to endogenous and/or environmental stressors. The coordination of mtDNA replication and transcription regulates the pace of mitochondrial biogenesis to guarantee the minimum number of organelles per cell. mtDNA inheritance follows a maternal lineage, although bi-parental inheritance has been reported in some species and in the case of mitochondrial diseases in humans. mtDNA mutations (e.g., point mutations, deletions, copy number variations) have been identified in the setting of several human diseases. For instance, sporadic and inherited rare disorders involving the nervous system as well higher risk of developing cancer and neurodegenerative conditions, including Parkinson's and Alzheimer's disease, have been associated with polymorphic mtDNA variants. An accrual of mtDNA mutations has also been identified in several tissues and organs, including heart and muscle, of old experimental animals and humans, which may contribute to the development of aging phenotypes. The role played by mtDNA homeostasis and mtDNA quality control pathways in human health is actively investigated for the possibility of developing targeted therapeutics for a wide range of conditions.

Metastatic pattern of ovarian cancer delineated by tracing the evolution of mitochondrial DNA mutations

Ovarian cancer (OC) is the most lethal gynecologic tumor and is characterized by a high rate of metastasis. Challenges in accurately delineating the metastatic pattern have greatly restricted the improvement of treatment in OC patients. An increasing number of studies have leveraged mitochondrial DNA (mtDNA) mutations as efficient lineage-tracing markers of tumor clonality. We applied multiregional sampling and high-depth mtDNA sequencing to determine the metastatic patterns in advanced-stage OC patients. Somatic mtDNA mutations were profiled from a total of 195 primary and 200 metastatic tumor tissue samples from 35 OC patients. Our results revealed remarkable sample-level and patient-level heterogeneity. In addition, distinct mtDNA mutational patterns were observed between primary and metastatic OC tissues. Further analysis identified the different mutational spectra between shared and private mutations among primary and metastatic OC tissues. Analysis of the clonality index calculated based on mtDNA mutations supported a monoclonal tumor origin in 14 of 16 patients with bilateral ovarian cancers. Notably, mtDNA-based spatial phylogenetic analysis revealed distinct patterns of OC metastasis, in which a linear metastatic pattern exhibited a low degree of mtDNA mutation heterogeneity and a short evolutionary distance, whereas a parallel metastatic pattern showed the opposite trend. Moreover, a mtDNA-based tumor evolutionary score (MTEs) related to different metastatic patterns was defined. Our data showed that patients with different MTESs responded differently to combined debulking surgery and chemotherapy. Finally, we observed that tumor-derived mtDNA mutations were more likely to be detected in ascitic fluid than in plasma samples. Our study presents an explicit view of the OC metastatic pattern, which sheds light on efficient treatment for OC patients.

Comparative analyses define differences between BHD-associated renal tumour and sporadic chromophobe renal cell carcinoma

BACKGROUND

Birt-Hogg-Dubé (BHD) syndrome, caused by germline alteration of folliculin (FLCN) gene, develops hybrid oncocytic/chromophobe tumour (HOCT) and chromophobe renal cell carcinoma (ChRCC), whereas sporadic ChRCC does not harbor FLCN alteration. To date, molecular characteristics of these similar histological types of tumours have been incompletely elucidated.

METHODS

To elucidate renal tumourigenesis of BHD-associated renal tumours and sporadic renal tumours, we conducted whole genome sequencing (WGS) and RNA-sequencing (RNA-seq) of sixteen BHD-associated renal tumours from nine unrelated BHD patients, twenty-one sporadic ChRCCs and seven sporadic oncocytomas. We then compared somatic mutation profiles with FLCN variants and RNA expression profiles between BHD-associated renal tumours and sporadic renal tumours.

FINDINGS

RNA-seq analysis revealed that BHD-associated renal tumours and sporadic renal tumours have totally different expression profiles. Sporadic ChRCCs were clustered into two distinct clusters characterized by L1CAM and FOXI1 expressions, molecular markers for renal tubule subclasses. Increased mitochondrial DNA (mtDNA) copy number with fewer variants was observed in BHD-associated renal tumours compared to sporadic ChRCCs. Cell-of-origin analysis using WGS data demonstrated that BHD-associated renal tumours and sporadic ChRCCs may arise from different cells of origin and second hit FLCN alterations may occur in early third decade of life in BHD patients.

INTERPRETATION

These data further our understanding of renal tumourigenesis of these two different types of renal tumours with similar histology.

FUNDING

This study was supported by JSPS KAKENHI Grants, RIKEN internal grant, and the Intramural Research Program of the National Institutes of Health (NIH), National Cancer Institute (NCI), Center for Cancer Research.

Targeting Mitochondrial DNA Transcription by POLRMT Inhibition or Depletion as a Potential Strategy for Cancer Treatment

Transcription of the mitochondrial genome is essential for the maintenance of oxidative phosphorylation (OXPHOS) and other functions directly related to this unique genome. Considerable evidence suggests that mitochondrial transcription is dysregulated in cancer and cancer metastasis and contributes significantly to cancer cell metabolism. Recently, inhibitors of the mitochondrial DNA-dependent RNA polymerase (POLRMT) were identified as potentially attractive new anti-cancer compounds. These molecules (IMT1, IMT1B) inactivate cancer cell metabolism through reduced transcription of mitochondrially-encoded OXPHOS subunits such as ND1-5 (Complex I) and COI-IV (Complex IV). Studies from our lab have discovered small molecule regulators of the mitochondrial matrix caseinolytic protease (ClpP) as probable inhibitors of mitochondrial transcription. These compounds activate ClpP proteolysis and lead to the rapid depletion of POLRMT and other matrix proteins, resulting in inhibition of mitochondrial transcription and growth arrest. Herein we present a comparison of POLRMT inhibition and ClpP activation, both conceptually and experimentally, and evaluate the results of these treatments on mitochondrial transcription, inhibition of OXPHOS, and ultimately cancer cell growth. We discuss the potential for targeting mitochondrial transcription as a cancer cell vulnerability.

Mitochondrial remodeling in colorectal cancer initiation, progression, metastasis, and therapy: A review

Colorectal cancer (CRC) is a major health concern with multifactorial pathophysiology representing intense therapeutic challenges. It is well known that deregulation of spatiotemporally-controlled signaling pathways and their metabolic reprogramming effects play a pivotal role in the development and progression of CRC. As such, the mitochondrial role in CRC initiation gained a lot of attention recently, as it is considered the powerhouse that regulates the bioenergetics in CRC. In addition, the crosstalk between microRNAs (miRNAs) and mitochondrial dysfunction has become a newfangled passion for deciphering CRC molecular mechanisms. This review sheds light on the relationship between different signaling pathways involved in metabolic reprogramming and their therapeutic targets, alterations in mitochondrial DNA content, mitochondrial biogenesis, and mitophagy, and the role of polymorphisms in mitochondrial genes as well as miRNAs regulating mitochondrial proteins in CRC initiation, progression, metastasis, and resistance to various therapies.

Quantitative haplotype-resolved analysis of mitochondrial DNA heteroplasmy in Human single oocytes, blastoids, and pluripotent stem cells

Maternal mitochondria are the sole source of mtDNA for every cell of the offspring. Heteroplasmic mtDNA mutations inherited from the oocyte are a common cause of metabolic diseases and associated with late-onset diseases. However, the origin and dynamics of mtDNA heteroplasmy remain unclear. We used our individual Mitochondrial Genome sequencing (iMiGseq) technology to study mtDNA heterogeneity, quantitate single nucleotide variants (SNVs) and large structural variants (SVs), track heteroplasmy dynamics, and analyze genetic linkage between variants at the individual mtDNA molecule level in single oocytes and human blastoids. Our study presented the first single-mtDNA analysis of the comprehensive heteroplasmy landscape in single human oocytes. Unappreciated levels of rare heteroplasmic variants well below the detection limit of conventional methods were identified in healthy human oocytes, of which many are reported to be deleterious and associated with mitochondrial disease and cancer. Quantitative genetic linkage analysis revealed dramatic shifts of variant frequency and clonal expansions of large SVs during oogenesis in single-donor oocytes. iMiGseq of a single human blastoid suggested stable heteroplasmy levels during early lineage differentiation of naïve pluripotent stem cells. Therefore, our data provided new insights of mtDNA genetics and laid a foundation for understanding mtDNA heteroplasmy at early stages of life.

Evaluation of Mitochondrial Function in Blood Samples Shows Distinct Patterns in Subjects with Thyroid Carcinoma from Those with Hyperplasia

Metabolic adaptations are a hallmark of cancer and may be exploited to develop novel diagnostic and therapeutic tools. Only about 50% of the patients who undergo thyroidectomy due to suspicion of thyroid cancer actually have the disease, highlighting the diagnostic limitations of current tools. We explored the possibility of using non-invasive blood tests to accurately diagnose thyroid cancer. We analyzed blood and thyroid tissue samples from two independent cohorts of patients undergoing thyroidectomy at the Hospital Universitario 12 de Octubre (Madrid, Spain). As expected, histological comparisons of thyroid cancer and hyperplasia revealed higher proliferation and apoptotic rates and enhanced vascular alterations in the former. Notably, they also revealed increased levels of membrane-bound phosphorylated AKT, suggestive of enhanced glycolysis, and alterations in mitochondrial sub-cellular distribution. Both characteristics are common metabolic adaptations in primary tumors. These data together with reduced mtDNA copy number and elevated levels of the mitochondrial antioxidant PRX3 in cancer tissue samples suggest the presence of mitochondrial oxidative stress. In plasma, cancer patients showed higher levels of cfDNA and mtDNA. Of note, mtDNA plasma levels inversely correlated with those in the tissue, suggesting that higher death rates were linked to lower mtDNA copy number. In PBMCs, cancer patients showed higher levels of PGC-1α, a positive regulator of mitochondrial function, but this increase was not associated with a corresponding induction of its target genes, suggesting a reduced activity in cancer patients. We also observed a significant difference in the PRDX3/PFKFB3 correlation at the gene expression level, between carcinoma and hyperplasia patients, also indicative of increased systemic metabolic stress in cancer patients. The correlation of mtDNA levels in tissue and PBMCs further stressed the interconnection between systemic and tumor metabolism. Evaluation of the mitochondrial gene ND1 in plasma, PBMCs and tissue samples, suggested that it could be a good biomarker for systemic oxidative metabolism, with ND1/mtDNA ratio positively correlating in PBMCs and tissue samples. In contrast, ND4 evaluation would be informative of tumor development, with ND4/mtDNA ratio specifically altered in the tumor context. Taken together, our data suggest that metabolic dysregulation in thyroid cancer can be monitored accurately in blood samples and might be exploited for the accurate discrimination of cancer from hyperplasia.

Tumour mitochondrial DNA mutations drive aerobic glycolysis to enhance checkpoint blockade

The mitochondrial genome encodes essential machinery for respiration and metabolic homeostasis but is paradoxically among the most common targets of somatic mutation in the cancer genome, with truncating mutations in respiratory complex I genes being most over-represented1. While mitochondrial DNA (mtDNA) mutations have been associated with both improved and worsened prognoses in several tumour lineages1-3, whether these mutations are drivers or exert any functional effect on tumour biology remains controversial. Here we discovered that complex I-encoding mtDNA mutations are sufficient to remodel the tumour immune landscape and therapeutic resistance to immune checkpoint blockade. Using mtDNA base editing technology4 we engineered recurrent truncating mutations in the mtDNA-encoded complex I gene, Mt-Nd5, into murine models of melanoma. Mechanistically, these mutations promoted utilisation of pyruvate as a terminal electron acceptor and increased glycolytic flux without major effects on oxygen consumption, driven by an over-reduced NAD pool and NADH shuttling between GAPDH and MDH1, mediating a Warburg-like metabolic shift. In turn, without modifying tumour growth, this altered cancer cell-intrinsic metabolism reshaped the tumour microenvironment in both mice and humans, promoting an anti-tumour immune response characterised by loss of resident neutrophils. This subsequently sensitised tumours bearing high mtDNA mutant heteroplasmy to immune checkpoint blockade, with phenocopy of key metabolic changes being sufficient to mediate this effect. Strikingly, patient lesions bearing >50% mtDNA mutation heteroplasmy also demonstrated a >2.5-fold improved response rate to checkpoint inhibitor blockade. Taken together these data nominate mtDNA mutations as functional regulators of cancer metabolism and tumour biology, with potential for therapeutic exploitation and treatment stratification.

MYC-driven increases in mitochondrial DNA copy number occur early and persist throughout prostatic cancer progression

Increased mitochondrial function may render some cancers vulnerable to mitochondrial inhibitors. Since mitochondrial function is regulated partly by mitochondrial DNA copy number (mtDNAcn), accurate measurements of mtDNAcn could help reveal which cancers are driven by increased mitochondrial function and may be candidates for mitochondrial inhibition. However, prior studies have employed bulk macrodissections that fail to account for cell type-specific or tumor cell heterogeneity in mtDNAcn. These studies have often produced unclear results, particularly in prostate cancer. Herein, we developed a multiplex in situ method to spatially quantify cell type specific mtDNAcn. We show that mtDNAcn is increased in luminal cells of high-grade prostatic intraepithelial neoplasia (HGPIN), is increased in prostatic adenocarcinomas (PCa), and is further elevated in metastatic castration-resistant prostate cancer. Increased PCa mtDNAcn was validated by two orthogonal methods and is accompanied by increases in mtRNAs and enzymatic activity. Mechanistically, MYC inhibition in prostate cancer cells decreases mtDNA replication and expression of several mtDNA replication genes, and MYC activation in the mouse prostate leads to increased mtDNA levels in the neoplastic prostate cells. Our in situ approach also revealed elevated mtDNAcn in precancerous lesions of the pancreas and colon/rectum, demonstrating generalization across cancer types using clinical tissue samples.

Mitochondrial Damage Induced by T-2 Mycotoxin on Human Skin-Fibroblast Hs68 Cell Line

T-2 toxin is produced by different Fusarium species and belongs to the group of type A trichothecene mycotoxins. T-2 toxin contaminates various grains, such as wheat, barley, maize, or rice, thus posing a risk to human and animal health. The toxin has toxicological effects on human and animal digestive, immune, nervous and reproductive systems. In addition, the most significant toxic effect can be observed on the skin. This in vitro study focused on T-2 toxicity on human skin fibroblast Hs68 cell line mitochondria. In the first step of this study, T-2 toxin's effect on the cell mitochondrial membrane potential (MMP) was determined. The cells were exposed to T-2 toxin, which resulted in dose- and time-dependent changes and a decrease in MMP. The obtained results revealed that the changes of intracellular reactive oxygen species (ROS) in the Hs68 cells were not affected by T-2 toxin. A further mitochondrial genome analysis showed that T-2 toxin in a dose- and time-dependent manner decreased the number of mitochondrial DNA (mtDNA) copies in cells. In addition, T-2 toxin genotoxicity causing mtDNA damage was evaluated. It was found that incubation of Hs68 cells in the presence of T-2 toxin, in a dose- and time-dependent manner, increased the level of mtDNA damage in both tested mtDNA regions: NADH dehydrogenase subunit 1 (ND1) and NADH dehydrogenase subunit 5 (ND5). In conclusion, the results of the in vitro study revealed that T-2 toxin shows adverse effects on Hs68 cell mitochondria. T-2 toxin induces mitochondrial dysfunction and mtDNA damage, which may cause the disruption of adenosine triphosphate (ATP) synthesis and, in consequence, cell death.

Single-strand mismatch and damage patterns revealed by single-molecule DNA sequencing

Mutations accumulate in the genome of every cell of the body throughout life, causing cancer and other genetic diseases1-4. Almost all of these mosaic mutations begin as nucleotide mismatches or damage in only one of the two strands of the DNA prior to becoming double-strand mutations if unrepaired or misrepaired5. However, current DNA sequencing technologies cannot resolve these initial single-strand events. Here, we developed a single-molecule, long-read sequencing method that achieves single-molecule fidelity for single-base substitutions when present in either one or both strands of the DNA. It also detects single-strand cytosine deamination events, a common type of DNA damage. We profiled 110 samples from diverse tissues, including from individuals with cancer-predisposition syndromes, and define the first single-strand mismatch and damage signatures. We find correspondences between these single-strand signatures and known double-strand mutational signatures, which resolves the identity of the initiating lesions. Tumors deficient in both mismatch repair and replicative polymerase proofreading show distinct single-strand mismatch patterns compared to samples deficient in only polymerase proofreading. In the mitochondrial genome, our findings support a mutagenic mechanism occurring primarily during replication. Since the double-strand DNA mutations interrogated by prior studies are only the endpoint of the mutation process, our approach to detect the initiating single-strand events at single-molecule resolution will enable new studies of how mutations arise in a variety of contexts, especially in cancer and aging.

Mitochondrial Metabolism in Pancreatic Ductal Adenocarcinoma: From Mechanism-Based Perspectives to Therapy

Pancreatic ductal adenocarcinoma (PDAC), the fourteenth most common malignancy, is a major contributor to cancer-related death with the utmost case fatality rate among all malignancies. Functional mitochondria, regardless of their complex ecosystem relative to normal cells, are essential in PDAC progression. Tumor cells' potential to produce ATP as energy, despite retaining the redox potential optimum, and allocating materials for biosynthetic activities that are crucial for cell growth, survival, and proliferation, are assisted by mitochondria. The polyclonal tumor cells with different metabolic profiles may add to carcinogenesis through inter-metabolic coupling. Cancer cells frequently possess alterations in the mitochondrial genome, although they do not hinder metabolism; alternatively, they change bioenergetics. This can further impart retrograde signaling, educate cell signaling, epigenetic modifications, chromatin structures, and transcription machinery, and ultimately satisfy cancer cellular and nuclear demands. To maximize the tumor microenvironment (TME), tumor cells remodel nearby stromal cells and extracellular matrix. These changes initiate polyclonality, which is crucial for growth, stress response, and metastasis. Here, we evaluate all the intrinsic and extrinsic pathways drawn by mitochondria in carcinogenesis, emphasizing the perspectives of mitochondrial metabolism in PDAC progression and treatment.

EFHD1, a novel mitochondrial regulator of tumor metastasis in clear cell renal cell carcinoma

The biological function of many mitochondrial proteins in mechanistic detail has not been well investigated in clear cell renal cell carcinoma (ccRCC). A seven-mitochondrial-gene signature was generated by Lasso regression analysis to improve the prediction of prognosis of patients with ccRCC, using The Cancer Genome Atlas and Clinical Proteomic Tumor Analysis Consortium cohort. Among those seven genes, EFHD1 is less studied and its role in the progression of ccRCC remains unknown. The decreased expression of EFHD1 was validated in clinical samples and was correlated with unfavorable outcome. Overexpression of EFHD1 in ccRCC cells resulted in the reduction of mitochondrial Ca²⁺ , and the inhibition of cell migration and invasion in vitro and tumor metastasis in vivo. Mechanistically, EFHD1 physically bound to the core mitochondrial calcium transporter (mitochondrial calcium uniporter, MCU) through its N-terminal domain. The interaction between EFHD1 and MCU suppressed the uptake of Ca²⁺ into mitochondria, and deactivated the Hippo/YAP signaling pathway. Further data revealed that the ectopic expression of EFHD1 upregulated STARD13 to enhance the phosphorylation of YAP protein at Ser-127. The knockdown of STARD13 or the overexpression of MCU partly abrogated the EFHD1-mediated induction of phosphorylation of YAP at Ser-127 and suppression of cell migration. Taken together, the newly identified EFHD1-MCU-STARD13 axis participates in the modulation of the Hippo/YAP pathway and serves as a novel regulator in the progression of ccRCC.

Variants in the Control Region of Mitochondrial Genome Associated with type 2 Diabetes in a Cohort of Mexican Mestizos

BACKGROUND

According to the International Diabetes Federation, Mexico is seventh place in the prevalence of type 2 diabetes (T2D) worldwide. Mitochondrial DNA variant association studies in multifactorial diseases like T2D are scarce in Mexican populations.

AIM OF THE STUDY

The objective of this study was to analyze the association between 18 variants in the mtDNA control region and T2D and related metabolic traits in a Mexican mestizo population from Mexico City.

METHODS

This study included 1001 participants divided into 477 cases with T2D and 524 healthy controls aged between 42 and 62 years and 18 mtDNA variants with frequencies >15%.

RESULTS

Association analyses matched by age and sex showed differences in the distribution between cases and controls for variants m.315_316insC (p = 1.18 × 10^-6), m.489T>C (p = 0.009), m.16362T>C (p = 0.001), and m.16519T>C (p = 0.004). The associations between T2D and variants m.315_316ins (OR = 6.13, CI = 3.42-10.97, p = 1.97 × 10^-6), m.489T>C (OR = 1.45, CI = 1.00-2.11, p = 0.006), m.16362T>C (OR = 2.17, CI = 1.57-3.00, p = 0.001), and m.16519T>C (OR = 1.69, CI = 1.23-2.33, p = 0.006) were significant after performing logistic regression models adjusted for age, sex, and diastolic blood pressure. Metabolic traits in the control group through linear regressions, adjusted for age, sex and BMI, and corrected for multiple comparisons showed nominal association between glucose and variants m.263A>G (p <0.050), m.16183A>C (p <0.010), m.16189T>C (p <0.020), and m.16223C>T (p <0.024); triglycerides, and cholesterol and variant m.309_310insC (p <0.010 and p <0.050 respectively); urea, and creatinine, and variant m.315_316insC (p <0.007, and p <0.004 respectively); diastolic blood pressure and variants m.235A>G (p <0.016), m.263A>G (p <0.013), m.315_316insC (p <0.043), and m.16111C>T (p <0.022).

CONCLUSION

These results demonstrate a strong association between variant m.315_316insC and T2D and a nominal association with T2D traits.

Spectrum of germline and somatic mitochondrial DNA variants in Tuberous Sclerosis Complex

Tuberous Sclerosis Complex (TSC) is caused by loss of function variants in either TSC1 or TSC2 and is characterized by broad phenotypic heterogeneity. Currently, there is limited knowledge regarding the role of the mitochondrial genome (mtDNA) in TSC pathogenesis. In this study, we aimed to determine the prevalence and spectrum of germline and somatic mtDNA variants in TSC and identify potential disease modifiers. Analysis of mtDNA amplicon massively parallel sequencing (aMPS) data, off-target mtDNA from whole-exome sequencing (WES), and/or qPCR, revealed mtDNA alterations in 270 diverse tissues (139 TSC-associated tumors and 131 normal tissue samples) from 199 patients and six healthy individuals. Correlation of clinical features to mtDNA variants and haplogroup analysis was done in 102 buccal swabs (age: 20-71 years). No correlation was found between clinical features and either mtDNA variants or haplogroups. No pathogenic variants were identified in the buccal swab samples. Using in silico analysis, we identified three predicted pathogenic variants in tumor samples: MT-ND4 (m.11742G>A, p. Cys328Tyr, VAF: 43%, kidney angiomyolipoma), MT-CYB (m.14775T>C, p. Leu10Pro, VAF: 43%, LAM abdominal tumor) and MT-CYB (m.15555C>T, p. Pro270Leu, VAF: 7%, renal cell carcinoma). Large deletions of the mitochondrial genome were not detected. Analysis of tumors from 23 patients with corresponding normal tissue did not reveal any recurrent tumor-associated somatic variants. The mtDNA/gDNA ratio between tumors and corresponding normal tissue was also unchanged. Overall, our findings demonstrate that the mitochondrial genome is highly stable across tissues and within TSC-associated tumors.

Targeting Mitochondria to Control Ageing and Senescence

Ageing is accompanied by a progressive impairment of cellular function and a systemic deterioration of tissues and organs, resulting in increased vulnerability to multiple diseases. Here, we review the interplay between two hallmarks of ageing, namely, mitochondrial dysfunction and cellular senescence. The targeting of specific mitochondrial features in senescent cells has the potential of delaying or even reverting the ageing process. A deeper and more comprehensive understanding of mitochondrial biology in senescent cells is necessary to effectively face this challenge. Here, we discuss the main alterations in mitochondrial functions and structure in both ageing and cellular senescence, highlighting the differences and similarities between the two processes. Moreover, we describe the treatments available to target these pathways and speculate on possible future directions of anti-ageing and anti-senescence therapies targeting mitochondria.

Genomic features of renal cell carcinoma developed during end-stage renal disease and dialysis

Patients with end-stage renal disease (ESRD) or receiving dialysis have a much higher risk for renal cell carcinoma (RCC), but carcinogenic mechanisms and genomic features remain little explored and undefined. This study's goal was to identify the genomic features of ESRD RCC and characterize them for associations with tumor histology and dialysis exposure. In this study, we obtained 33 RCCs, with various histological subtypes, that developed in ESRD patients receiving dialysis and performed whole-genome sequencing and transcriptome analyses. Driver events, copy-number alteration (CNA) analysis and mutational signature profiling were performed using an analysis pipeline that integrated data from germline and somatic SNVs, Indels and structural variants as well as CNAs, while transcriptome data were analyzed for differentially expressed genes and through gene set enrichment analysis. ESRD related clear cell RCCs' driver genes and mutations mirrored those in sporadic ccRCCs. Longer dialysis periods significantly correlated with a rare mutational signature SBS23, whose etiology is unknown, and increased mitochondrial copy number. All acquired cystic disease (ACD)-RCCs, which developed specifically in ESRD patients, showed chromosome 16q amplification. Gene expression analysis suggests similarity between certain ACD-RCCs and papillary RCCs and in TCGA papillary RCCs with chromosome 16 gain identified enrichment for genes related to DNA repair, as well as pathways related to reactive oxygen species, oxidative phosphorylation and targets of Myc. This analysis suggests that ESRD or dialysis could induce types of cellular stress that impact some specific types of genomic damage leading to oncogenesis.

Multiple metals exposure and blood mitochondrial DNA copy number: A cross-sectional study from the Dongfeng-Tongji cohort

OBJECTIVE

Mitochondria are essential organelles that execute fundamental biological processes, while mitochondrial DNA is vulnerable to environmental insults. The aim of this study was to investigate the individual and mixture effect of plasma metals on blood mitochondria DNA copy number (mtDNAcn).

METHODS

This study involved 1399 randomly selected subcohort participants from the Dongfeng-Tongji cohort. The blood mtDNAcn and plasma levels of 23 metals were determined by using quantitative real-time polymerase chain reaction (qPCR) and inductively coupled plasma mass spectrometer (ICP-MS), respectively. The multiple linear regression was used to explore the association between each metal and mtDNAcn, and the LASSO penalized regression was performed to select the most significant metals. We also used the quantile g-computation analysis to assess the mixture effect of multiple metals.

RESULTS

Based on multiple linear regression models, each 1% increase in plasma concentration of copper (Cu), rubidium (Rb), and titanium (Ti) was associated with a separate 0.16% [β(95% CI) = 0.158 (0.066, 0.249), P = 0.001], 0.20% [β(95% CI) = 0.196 (0.073, 0.318), P = 0.002], and 0.25% [β(95% CI) = 0.245 (0.081, 0.409), P = 0.003] increase in blood mtDNAcn. The LASSO regression also confirmed Cu, Rb, and Ti as significant predictors for mtDNAcn. There was a significant mixture effect of multiple metals on increasing mtDNAcn among the elder participants (aged ≥65), with an approximately 11% increase in mtDNAcn for each quartile increase in all metal concentrations [β(95% CI) = 0.146 (0.048, 0.243), P = 0.004].

CONCLUSIONS

Our results show that plasma Cu, Rb and Ti were associated with increased blood mtDNA, and we further revealed a significant mixture effect of all metals on mtDNAcn among elder population. These findings may provide a novel perspective on the effect of metals on mitochondrial dysfunction.

MicroRNA-Mediated Mitochondrial Dysfunction Is Involved in the Anti-triple-Negative Breast Cancer Cell Activity of Phytosesquiterpene Lactones

Aims: Emerging evidence suggests that modulating redox homeostasis through targeting mitochondrial functions may be a useful strategy for suppressing triple-negative breast cancer (TNBC) activities. However, whether there are specific microRNAs (miRNAs) involved in regulating oxidative stress-associated mitochondrial functions that can act as therapeutic targets to suppress TNBC activities remains unclear. Here, we aimed to identify the role of redox-associated miRNAs in TNBC and investigated their potential as therapeutic targets. Results: We identified oxidative stress-responsive differentially expressed miRNAs (DEMs) regulated by phytosesquiterpene lactone deoxyelephantopin (DET) and its novel derivative DETD-35, which are known to inhibit TNBC growth and metastasis in vitro and in vivo, using comparative miRNA microarray analysis and reactive oxygen species (ROS) scavenging approaches. Mitochondrial dysfunction was identified as a major biological function regulated by a few specific DEMs. In particular, miR-4284 was identified to play a role in DET- and DETD-35-mediated ROS production, mitochondrial basal proton leak, and antiproliferation activity in TNBC cells. Moreover, DET- and DETD-35-induced mitochondrial DNA damage was observed in TNBC cells and xenograft tumors. miR-4284 was also identified to play a role in oxidative DNA damage in TNBC tumors. Innovation: We identified a novel role for miR-4284 in regulating mitochondrial basal proton leak in TNBC cells, and highlighted its significance in TNBC tumor oxidative DNA damage, and its direct correlation with TNBC patient survival. Conclusion: We used DET and DETD-35 as proof of concept to demonstrate that activities of anticancer agents can involve regulation of multiple miRNAs playing different roles in cancer progression. Antioxid. Redox Signal. 38, 198-214.

Integrated bioinformatic analysis of mitochondrial metabolism-related genes in acute myeloid leukemia

Background

Acute myeloid leukemia (AML) is a common hematologic malignancy characterized by poor prognoses and high recurrence rates. Mitochondrial metabolism has been increasingly recognized to be crucial in tumor progression and treatment resistance. The purpose of this study was to examined the role of mitochondrial metabolism in the immune regulation and prognosis of AML.

Methods

In this study, mutation status of 31 mitochondrial metabolism-related genes (MMRGs) in AML were analyzed. Based on the expression of 31 MMRGs, mitochondrial metabolism scores (MMs) were calculated by single sample gene set enrichment analysis. Differential analysis and weighted co-expression network analysis were performed to identify module MMRGs. Next, univariate Cox regression and the least absolute and selection operator regression were used to select prognosis-associated MMRGs. A prognosis model was then constructed using multivariate Cox regression to calculate risk score. We validated the expression of key MMRGs in clinical specimens using immunohistochemistry (IHC). Then differential analysis was performed to identify differentially expressed genes (DEGs) between high- and low-risk groups. Functional enrichment, interaction networks, drug sensitivity, immune microenvironment, and immunotherapy analyses were also performed to explore the characteristic of DEGs.

Results

Given the association of MMs with prognosis of AML patients, a prognosis model was constructed based on 5 MMRGs, which could accurately distinguish high-risk patients from low-risk patients in both training and validation datasets. IHC results showed that MMRGs were highly expressed in AML samples compared to normal samples. Additionally, the 38 DEGs were mainly related to mitochondrial metabolism, immune signaling, and multiple drug resistance pathways. In addition, high-risk patients with more immune-cell infiltration had higher Tumor Immune Dysfunction and Exclusion scores, indicating poor immunotherapy response. mRNA-drug interactions and drug sensitivity analyses were performed to explore potential druggable hub genes. Furthermore, we combined risk score with age and gender to construct a prognosis model, which could predict the prognosis of AML patients.

Conclusion

Our study provided a prognostic predictor for AML patients and revealed that mitochondrial metabolism is associated with immune regulation and drug resistant in AML, providing vital clues for immunotherapies.

Mutational signature of mtDNA confers mechanistic insight into oxidative metabolism remodeling in colorectal cancer

Rationale: Mitochondrial dysfunction caused by mitochondrial DNA (mtDNA) mutations and subsequent metabolic defects are closely involved in tumorigenesis and progression in a cancer-type specific manner. To date, the mutational pattern of mtDNA somatic mutations in colorectal cancer (CRC) tissues and its clinical implication are still not completely clear. Methods: In the present study, we generated a large mtDNA somatic mutation dataset from three CRC cohorts (432, 1,015, and 845 patients, respectively) and then most comprehensively characterized the CRC-specific evolutionary pattern and its clinical implication. Results: Our results showed that the mtDNA control region (mtCTR) with a high mutation density exhibited a distinct mutation spectrum characterizing a high enrichment of L-strand C > T mutations, which was contrary to the H-strand C > T mutational bias observed in the mtDNA coding region (mtCDR) (P < 0.001). Further analysis clearly confirmed the relaxed evolutionary selection of mtCTR mutations, which was mainly characterized by the similar distribution of hypervariable region (HVS) and non-HVS mutation density. Moreover, significant negative selection was identified in mutations of mtDNA complex V (ATP6/ATP8) and tRNA loop regions. Although our data showed that oxidative metabolism was commonly increased in CRC cells, mtDNA somatic mutations in CRC tissues were not closely associated with mitochondrial biogenesis, oxidative metabolism, and clinical progression, suggesting a cancer-type specific relationship between mtDNA mutations and mitochondrial metabolic functions in CRC cells. Conclusion: Our study identified the CRC-specific evolutionary mode of mtDNA mutations, which is possibly matched to specific mitochondrial metabolic remodeling and confers new mechanic insight into CRC tumorigenesis.

Acetyl-CoA regulates lipid metabolism and histone acetylation modification in cancer

Acetyl-CoA, as an important molecule, not only participates in multiple intracellular metabolic reactions, but also affects the post-translational modification of proteins, playing a key role in the metabolic activity and epigenetic inheritance of cells. Cancer cells require extensive lipid metabolism to fuel for their growth, while also require histone acetylation modifications to increase the expression of cancer-promoting genes. As a raw material for de novo lipid synthesis and histone acetylation, acetyl-CoA has a major impact on lipid metabolism and histone acetylation in cancer. More importantly, in cancer, acetyl-CoA connects lipid metabolism with histone acetylation, forming a more complex regulatory mechanism that influences cancer growth, proliferation, metastasis.

Mitochondrial DNA Mutations and Ageing

Mitochondria are subcellular organelles present in most eukaryotic cells which play a significant role in numerous aspects of cell biology. These include carbohydrate and fatty acid metabolism to generate cellular energy through oxidative phosphorylation, apoptosis, cell signalling, haem biosynthesis and reactive oxygen species production. Mitochondrial dysfunction is a feature of many human ageing tissues, and since the discovery that mitochondrial DNA mutations were a major underlying cause of changes in oxidative phosphorylation capacity, it has been proposed that they have a role in human ageing. However, there is still much debate on whether mitochondrial DNA mutations play a causal role in ageing or are simply a consequence of the ageing process. This chapter describes the structure of mammalian mitochondria, and the unique features of mitochondrial genetics, and reviews the current evidence surrounding the role of mitochondrial DNA mutations in the ageing process. It then focusses on more recent discoveries regarding the role of mitochondrial dysfunction in stem cell ageing and age-related inflammation.

Maximizing Small Biopsy Patient Samples: Unified RNA-Seq Platform Assessment of over 120,000 Patient Biopsies

Despite its wide-ranging benefits, whole-transcriptome or RNA exome profiling is challenging to implement in a clinical diagnostic setting. The Unified Assay is a comprehensive workflow wherein exome-enriched RNA-sequencing (RNA-Seq) assays are performed on clinical samples and analyzed by a series of advanced machine learning-based classifiers. Gene expression signatures and rare and/or novel genomic events, including fusions, mitochondrial variants, and loss of heterozygosity were assessed using RNA-Seq data generated from 120,313 clinical samples across three clinical indications (thyroid cancer, lung cancer, and interstitial lung disease). Since its implementation, the data derived from the Unified Assay have allowed significantly more patients to avoid unnecessary diagnostic surgery and have played an important role in guiding follow-up decisions regarding treatment. Collectively, data from the Unified Assay show the utility of RNA-Seq and RNA expression signatures in the clinical laboratory, and their importance to the future of precision medicine.

Roles of Mitochondrial DNA Damage in Kidney Diseases: A New Biomarker

The kidney is a mitochondria-rich organ, and kidney diseases are recognized as mitochondria-related pathologies. Intact mitochondrial DNA (mtDNA) maintains normal mitochondrial function. Mitochondrial dysfunction caused by mtDNA damage, including impaired mtDNA replication, mtDNA mutation, mtDNA leakage, and mtDNA methylation, is involved in the progression of kidney diseases. Herein, we review the roles of mtDNA damage in different setting of kidney diseases, including acute kidney injury (AKI) and chronic kidney disease (CKD). In a variety of kidney diseases, mtDNA damage is closely associated with loss of kidney function. The level of mtDNA in peripheral serum and urine also reflects the status of kidney injury. Alleviating mtDNA damage can promote the recovery of mitochondrial function by exogenous drug treatment and thus reduce kidney injury. In short, we conclude that mtDNA damage may serve as a novel biomarker for assessing kidney injury in different causes of renal dysfunction, which provides a new theoretical basis for mtDNA-targeted intervention as a therapeutic option for kidney diseases.

Roles of mitochondrial genetics in cancer metastasis

The contributions of mitochondria to cancer have been recognized for decades. However, the focus on the metabolic role of mitochondria and the diminutive size of the mitochondrial genome compared to the nuclear genome have hindered discovery of the roles of mitochondrial genetics in cancer. This review summarizes recent data demonstrating the contributions of mitochondrial DNA (mtDNA) copy-number variants (CNVs), somatic mutations, and germline polymorphisms to cancer initiation, progression, and metastasis. The goal is to summarize accumulating data to establish a framework for exploring the contributions of mtDNA to neoplasia and metastasis.

Mitochondrial DNA Copy Number: Linking Diabetes and Cancer

Recent Advances: Various studies have suggested that mitochondrial DNA copy number (mtDNA-CN), a surrogate biomarker of mitochondrial dysfunction, is an easily quantifiable biomarker for chronic diseases, including diabetes and cancer. However, current knowledge is limited, and the results are controversial. This has been attributed mainly to methodology and study design. Critical Issues: The incidence of diabetes and cancer has increased significantly in recent years. Moreover, type 2 diabetes (T2D) has been shown to be a risk factor for cancer. mtDNA-CN has been associated with both T2D and cancer. However, it is not known whether mtDNA-CN plays any role in the association between T2D and cancer. Significance: In this review, we have discussed mtDNA-CN in diabetes and cancer, and reviewed the literature and methodology used in published studies so far. Based on the literature review, we have speculated how mtDNA-CN may act as a link between diabetes and cancer. Furthermore, we have provided some recommendations for reliable translation of mtDNA-CN as a biomarker. Future Directions: Further research is required to elucidate the role of mtDNA-CN in the association between T2D and cancer. If established, early lifestyle interventions, such as physical activity and diet control that improve mitochondrial function, may help preventing cancer in patients with T2D. Antioxid. Redox Signal. 37, 1168-1190.

Mitochondrial DNA is a major source of driver mutations in cancer

Mitochondrial DNA (mtDNA) mutations are among the most common genetic events in all tumors and directly impact metabolic homeostasis. Despite the central role mitochondria play in energy metabolism and cellular physiology, the role of mutations in the mitochondrial genomes of tumors has been contentious. Until recently, genomic and functional studies of mtDNA variants were impeded by a lack of adequate tumor mtDNA sequencing data and available methods for mitochondrial genome engineering. These barriers and a conceptual fog surrounding the functional impact of mtDNA mutations in tumors have begun to lift, revealing a path to understanding the role of this essential metabolic genome in cancer initiation and progression. Here we discuss the history, recent developments, and challenges that remain for mitochondrial oncogenetics as the impact of a major new class of cancer-associated mutations is unveiled.

Tobacco smoke exposure and mitochondrial DNA copy number on neurobehavioural performance: A community study

The influence of tobacco smoke has been a controversial and very questionable subject within the field of neurological behaviours. To examine the dose-response relationships between tobacco smoke and neurological performance, we investigated whether mitochondrial DNA copy number (mtDNAcn) mediates these relationships. We used restricted cubic spline models to estimate the dose-response relationships. A mediation model was also used to detect the mediating effect. Increased cotinine was negatively associated with auditory memory scores and a 0.51 decrease in mtDNAcn. MtDNAcn acts as a mediator between cotinine and auditory memory. Tobacco smoke levels were inversely associated with mtDNAcn and neurobehavioural changes, and there was a mediation effect between cotinine levels and auditory memory by mtDNAcn.

Detection of mitochondrial DNA mutations in circulating mitochondria-originated extracellular vesicles for potential diagnostic applications in pancreatic adenocarcinoma

There is a complete lack of highly sensitive and specific biomarkers for early pancreatic ductal adenocarcinoma (PDAC) diagnosis, limiting multi-modal therapeutic options. Mitochondrial DNA (mtDNA) is an excellent resource for biomarker discovery because of its high copy number and increased mutational frequency in cancer cells. We examined if mtDNA mutations can be detected in circulating extracellular vesicles (EVs) of PDAC patients and used for discerning between cancer and non-cancer subjects. A greater yield of circulating EVs (~ 1.4 fold; p = 0.002) was obtained in PDAC patients (n = 20) than non-cancer (NC) individuals (n = 10). PDAC-EVs contained a higher quantity of total DNA (~ 5.5 folds; p = 0.0001) than NC-EVs and had greater enrichment of mtDNA (~ 14.02-fold; p = 0.0001). PDAC-EVs also had higher levels of cardiolipin (a mitochondrial inner-membrane phospholipid), suggestive of their mitochondrial origin. All mtDNA mutations in PDAC-EVs were unique and frequency was remarkably higher. Most mtDNA mutations (41.5%) in PDAC-EVs were in the respiratory complex-I (RCI) (ND1-ND6), followed by the RCIII gene (CYTB; 11.2%). Among the non-coding genes, D-Loop and RNR2 exhibited the most mutations (15.2% each). Altogether, our study establishes, for the first time, that mtDNA mutations can be detected in circulating EVs and potentially serve as a tool for reliable PDAC diagnosis.

Nuclear-embedded mitochondrial DNA sequences in 66,083 human genomes

DNA transfer from cytoplasmic organelles to the cell nucleus is a legacy of the endosymbiotic event-the majority of nuclear-mitochondrial segments (NUMTs) are thought to be ancient, preceding human speciation^1-3. Here we analyse whole-genome sequences from 66,083 people-including 12,509 people with cancer-and demonstrate the ongoing transfer of mitochondrial DNA into the nucleus, contributing to a complex NUMT landscape. More than 99% of individuals had at least one of 1,637 different NUMTs, with 1 in 8 individuals having an ultra-rare NUMT that is present in less than 0.1% of the population. More than 90% of the extant NUMTs that we evaluated inserted into the nuclear genome after humans diverged from apes. Once embedded, the sequences were no longer under the evolutionary constraint seen within the mitochondrion, and NUMT-specific mutations had a different mutational signature to mitochondrial DNA. De novo NUMTs were observed in the germline once in every 10⁴ births and once in every 10³ cancers. NUMTs preferentially involved non-coding mitochondrial DNA, linking transcription and replication to their origin, with nuclear insertion involving multiple mechanisms including double-strand break repair associated with PR domain zinc-finger protein 9 (PRDM9) binding. The frequency of tumour-specific NUMTs differed between cancers, including a probably causal insertion in a myxoid liposarcoma. We found evidence of selection against NUMTs on the basis of size and genomic location, shaping a highly heterogenous and dynamic human NUMT landscape.

Precise and simultaneous quantification of mitochondrial DNA heteroplasmy and copy number by digital PCR

Mitochondrial DNA (mtDNA) is present in multiple copies and phenotypic consequences of mtDNA mutations depend on the mutant load surpassing a specific threshold. Additionally, changes in mtDNA copy number can impact mitochondrial ATP production, resulting in disease. Therefore, the precise determination of mtDNA heteroplasmy and copy number is crucial to the study of mitochondrial diseases. However, current methods can be imprecise, and quantifying small changes in either heteroplasmy or copy number is challenging. We developed a new approach to measure mtDNA heteroplasmy using a single digital PCR (dPCR) probe. This method is based on the observation that fluorescent-labeled probes in dPCR exhibit different intensities depending on the presence of a single nucleotide change in the sequence bound by the probe. This finding allowed us to precisely and simultaneously determine mtDNA copy number and heteroplasmy levels using duplex dPCR. We tested this approach in two different models (human and mouse), which proved faster and more internally controlled when compared to other published methods routinely used in the mitochondrial genetics field. We believe this approach could be broadly applicable to the detection and quantification of other mixed genetic variations.

Modulating mitochondrial DNA mutations: factors shaping heteroplasmy in the germ line and somatic cells

Until recently it was thought that most humans only harbor one type of mitochondrial DNA (mtDNA), however, deep sequencing and single-cell analysis has shown the converse - that mixed populations of mtDNA (heteroplasmy) are the norm. This is important because heteroplasmy levels can change dramatically during transmission in the female germ line, leading to high levels causing severe mitochondrial diseases. There is also emerging evidence that low level mtDNA mutations contribute to common late onset diseases such as neurodegenerative disorders and cardiometabolic diseases because the inherited mutation levels can change within developing organs and non-dividing cells over time. Initial predictions suggested that the segregation of mtDNA heteroplasmy was largely stochastic, with an equal tendency for levels to increase or decrease. However, transgenic animal work and single-cell analysis have shown this not to be the case during germ-line transmission and in somatic tissues during life. Mutation levels in specific mtDNA regions can increase or decrease in different contexts and the underlying molecular mechanisms are starting to be unraveled. In this review we provide a synthesis of recent literature on the mechanisms of selection for and against mtDNA variants. We identify the most pertinent gaps in our understanding and suggest ways these could be addressed using state of the art techniques.

Lung mitochondrial DNA copy number, inflammatory biomarkers, gene transcription and gene methylation in vapers and smokers

BACKGROUND

Mitochondrial DNA copy number (mtCN) maintains cellular function and homeostasis, and is linked to nuclear DNA methylation and gene expression. Increased mtCN in the blood is associated with smoking and respiratory disease, but has received little attention for target organ effects for smoking or electronic cigarette (EC) use.

METHODS

Bronchoscopy biospecimens from healthy EC users, smokers (SM), and never-smokers (NS) were assessed for associations of mtCN with mtDNA point mutations, immune responses, nuclear DNA methylation and gene expression using linear regression. Ingenuity pathway analysis was used for enriched pathways. GEO and TCGA respiratory disease datasets were used to explore the involvement of mtCN-associated signatures.

FINDINGS

mtCN was higher in SM than NS, but EC was not statistically different from either. Overall there was a negative association of mtCN with a point mutation in the D-loop but no difference within groups. Positive associations of mtCN with IL-2 and IL-4 were found in EC only. mtCN was significantly associated with 71,487 CpGs and 321 transcripts. 263 CpGs were correlated with nearby transcripts for genes enriched in the immune system. EC-specific mtCN-associated-CpGs and genes were differentially expressed in respiratory diseases compared to controls, including genes involved in cellular movement, inflammation, metabolism, and airway hyperresponsiveness.

INTERPRETATION

Smoking may elicit a lung toxic effect through mtCN. While the impact of EC is less clear, EC-specific associations of mtCN with nuclear biomarkers suggest exposure may not be harmless. Further research is needed to understand the role of smoking and EC-related mtCN on lung disease risks.

FUNDING

The National Cancer Institute, the National Heart, Lung, and Blood Institute, the Food and Drug Administration Center for Tobacco Products, the National Center For Advancing Translational Sciences, and Pelotonia Intramural Research Funds.

A novel signature model based on mitochondrial-related genes for predicting survival of colon adenocarcinoma

Background

Colon cancer is the foremost reason of cancer-related mortality worldwide. Colon adenocarcinoma constitutes 90% of colon cancer, and most patients with colon adenocarcinoma (COAD) are identified until advanced stage. With the emergence of an increasing number of novel pathogenic mechanisms and treatments, the role of mitochondria in the development of cancer, has been studied and reported with increasing frequency.

Methods

We systematically analyzed the effect of mitochondria-related genes in COAD utilizing RNA sequencing dataset from The Cancer Genome Atlas database and 1613 mitochondrial function-related genes from MitoMiner database. Our approach consisted of differentially expressed gene, gene set enrichment analysis, gene ontology terminology, Kyoto Encyclopedia of Genes and Genomes, independent prognostic analysis, univariate and multivariate analysis, Kaplan–Meier survival analysis, immune microenvironment correlation analysis, and Cox regression analysis.

Results

Consequently, 8 genes were identified to construct 8 mitochondrial-related gene model by applying Cox regression analysis, CDC25C, KCNJ11, NOL3, P4HA1, QSOX2, Trap1, DNAJC28, and ATCAY. Meanwhile, we assessed the connection between this model and clinical parameters or immune microenvironment. Risk score was an independent predictor for COAD patients’ survival with an AUC of 0.687, 0.752 and 0.762 at 1-, 3- and 5-year in nomogram, respectively. The group with the highest risk score had the lowest survival rate and the worst clinical stages. Additionally, its predictive capacity was validated in GSE39582 cohort.

Conclusion

In summary, we established a prognostic pattern of mitochondrial-related genes, which can predict overall survival in COAD, which may enable a more optimized approach for the clinical treatment and scientific study of COAD. This gene signature model has the potential to improve prognosis and treatment for COAD patients in the future, and to be widely implemented in clinical settings. The utilization of this mitochondrial-related gene signature model may be benefit in the treatments and medical decision-making of COAD.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12911-022-02020-3.

A mitochondria-specific mutational signature of aging: increased rate of A &gt; G substitutions on the heavy strand

The mutational spectrum of the mitochondrial DNA (mtDNA) does not resemble any of the known mutational signatures of the nuclear genome and variation in mtDNA mutational spectra between different organisms is still incomprehensible. Since mitochondria are responsible for aerobic respiration, it is expected that mtDNA mutational spectrum is affected by oxidative damage. Assuming that oxidative damage increases with age, we analyse mtDNA mutagenesis of different species in regards to their generation length. Analysing, (i) dozens of thousands of somatic mtDNA mutations in samples of different ages (ii) 70053 polymorphic synonymous mtDNA substitutions reconstructed in 424 mammalian species with different generation lengths and (iii) synonymous nucleotide content of 650 complete mitochondrial genomes of mammalian species we observed that the frequency of AH > GH substitutions (H: heavy strand notation) is twice bigger in species with high versus low generation length making their mtDNA more AH poor and GH rich. Considering that AH > GH substitutions are also sensitive to the time spent single-stranded (TSSS) during asynchronous mtDNA replication we demonstrated that AH > GH substitution rate is a function of both species-specific generation length and position-specific TSSS. We propose that AH > GH is a mitochondria-specific signature of oxidative damage associated with both aging and TSSS.

A method for multiplexed full-length single-molecule sequencing of the human mitochondrial genome

Methods to reconstruct the mitochondrial DNA (mtDNA) sequence using short-read sequencing come with an inherent bias due to amplification and mapping. They can fail to determine the phase of variants, to capture multiple deletions and to cover the mitochondrial genome evenly. Here we describe a method to target, multiplex and sequence at high coverage full-length human mitochondrial genomes as native single-molecules, utilizing the RNA-guided DNA endonuclease Cas9. Combining Cas9 induced breaks, that define the mtDNA beginning and end of the sequencing reads, as barcodes, we achieve high demultiplexing specificity and delineation of the full-length of the mtDNA, regardless of the structural variant pattern. The long-read sequencing data is analysed with a pipeline where our custom-developed software, baldur, efficiently detects single nucleotide heteroplasmy to below 1%, physically determines phase and can accurately disentangle complex deletions. Our workflow is a tool for studying mtDNA variation and will accelerate mitochondrial research.

svaRetro and svaNUMT: modular packages for annotating retrotransposed transcripts and nuclear integration of mitochondrial DNA in genome sequencing data

Nuclear integration of mitochondrial genomes and retrocopied transcript insertion are biologically important but often-overlooked aspects of structural variant (SV) annotation. While tools for their detection exist, these typically rely on reanalysis of primary data using specialised detectors rather than leveraging calls from general purpose structural variant callers. Such reanalysis potentially leads to additional computational expense and does not take advantage of advances in general purpose structural variant calling. Here, we present svaRetro and svaNUMT; R packages that provide functions for annotating novel genomic events, such as nonreference retrocopied transcripts and nuclear integration of mitochondrial DNA. The packages were developed to work within the Bioconductor framework. We evaluate the performance of these packages to detect events using simulations and public benchmarking datasets, and annotate processed transcripts in a public structural variant database. svaRetro and svaNUMT provide modular, SV-caller agnostic tools for downstream annotation of structural variant calls.

Aging-associated accumulation of mitochondrial DNA mutations in tumor origin

The majority of cancer patients are among aged population, suggesting an urgent need to advance our knowledge on complicated relationship between aging and cancer. It has been hypothesized that metabolic changes during aging could act as a driver for tumorigenesis. Given the fact that mitochondrial DNA (mtDNA) mutations are common in both tumors and aged tissues, it is interesting to contemplate possible role of age-related mtDNA mutations in tumorigenesis. MtDNA encodes genes essential for mitochondrial metabolism, and mtDNA mutates at a much higher rate than nuclear genome. Random drifting of somatic mtDNA mutations, as a result of cell division or mitochondrial turnover during aging, may lead to more and more cells harboring high-frequency pathogenic mtDNA mutations, albeit at different loci, in single-cells. Such mutations can induce metabolic reprogramming, nuclear genome instability and immune response, which might increase the likelihood of tumorigenesis. In this review, we summarize current understanding of how mtDNA mutations accumulate with aging and how these mutations could mechanistically contribute to tumor origin. We also discuss potential prevention strategies for mtDNA mutation-induced tumorigenesis, and future works needed in this direction.

Single nucleotide polymorphism patterns associated with a cancer resistant phenotype

BACKGROUND AND AIMS

In this study ten mouse strains representing ~90% of genetic diversity in laboratory mice (B6C3F1/J, C57BL/6J, C3H/HeJ, A/J, NOD.B1oSnH2/J, NZO/HILtJ, 129S1/SvImJ, WSB/EiJ, PWK/PhJ, CAST/EiJ) were examined to identify the mouse strain with the lowest incidence of cancer. The unique single polymorphisms (SNPs) associated with this low cancer incidence are reported.

METHODS

Evaluations of cancer incidence in the 10 mouse strains were based on gross and microscopic diagnosis of tumors. Single nucleotide polymorphisms (SNPs) in the coding regions of the genome were derived from the respective mouse strains located in the Sanger mouse sequencing database and the B6C3F1/N genome from the National Toxicology Program (NTP).

RESULTS

The WSB strain had an overall lower incidence of both benign and malignant tumors compared to the other mouse strains. At 2 years, the incidence of total malignant tumors (Poly-3 incidence rate) ranged from 2% (WSB) to 92% (C3H) in males, and 14% (WSB) to 93% (NZO) in females, and the total incidence of benign and malignant tumor incidence ranged from 13% (WSB) to 99% (C3H) in males and 25% (WSB) to 96% (NOD) in females. Single nucleotide polymorphism (SNP) patterns were examined in the following strains: B6C3F1/N, C57BL/6J, C3H/HeJ, 129S1/SvImJ, A/J, NZO/HILtJ, CAST/EiJ, PWK/PhJ, and WSB/EiJ. We identified 7519 SNPs (involving 5751 Ensembl transcripts of 3453 Ensembl Genes) that resulted in a unique amino acid change in the coding region of the WSB strain.

CONCLUSIONS

The inherited genetic patterns in the WSB cancer-resistant mouse strain occurred in genes involved in multiple cell functions including mitochondria, metabolic, immune, and membrane-related cell functions. The unique SNP patterns in a cancer resistant mouse strain provides insights for understanding and developing strategies for cancer prevention.

Targeting Mitochondrial Function with Chemoptogenetics

Mitochondria are ATP-generating organelles in eukaryotic cells that produce reactive oxygen species (ROS) during oxidative phosphorylation (OXPHOS). Mitochondrial DNA (mtDNA) is packaged within nucleoids and, due to its close proximity to ROS production, endures oxidative base damage. This damage can be repaired by base excision repair (BER) within the mitochondria, or it can be degraded via exonucleases or mitophagy. Persistent mtDNA damage may drive the production of dysfunctional OXPHOS components that generate increased ROS, or OXPHOS components may be directly damaged by ROS, which then can cause more mtDNA damage and create a vicious cycle of ROS production and mitochondrial dysfunction. If mtDNA damage is left unrepaired, mtDNA mutations including deletions can result. The accumulation of mtDNA mutations has been associated with conditions ranging from the aging process to cancer and neurodegenerative conditions, but the sequence of events leading to mtDNA mutations and deletions is yet unknown. Researchers have utilized many systems and agents for generating ROS in mitochondria to observe the downstream effects on mtDNA, ROS, and mitochondrial function; yet, there are various drawbacks to these methodologies that limit their precision. Here, we describe a novel chemoptogenetic approach to target oxidative damage to mitochondria and mtDNA with a high spatial and temporal resolution so that the downstream effects of ROS-induced damage can be measured with a high precision in order to better understand the mechanism of mitochondrial dysfunction in aging, cancer, and neurodegenerative diseases.