Mitochondrial DNA copy number variation across human cancers

Physiological Cell Culture Media Tune Mitochondrial Bioenergetics and Drug Sensitivity in Cancer Cell Models

Simple Summary

Cell biologists trust in standard media for analyzing cellular functions and for the specification of target-oriented drug efficacies in cell culture settings. Here, we present a general applicable workflow for the constant monitoring of bioenergetic states of cells grown in 2D cell models to accompany tailored drug discovery efforts. Using in-depth high-resolution respirometry analyses (HRR) of mitochondrial function, we unveiled that the human-plasma-like media (HPLM) altered cellular energetic states. In a systematic HRR setup for drug profiling experiments, we revealed an unexpected side effect of an FDA-approved cancer drug on mitochondrial function, exclusively in HPLM. Thus, we believe that both the recordings of bioenergetic states and the use of more physiological media would improve and reshape cell-based drug discovery ventures.

Abstract

Two-dimensional cell cultures are established models in research for studying and perturbing cell-type specific functions. However, many limitations apply to the cell growth in a monolayer using standard cell culture media. Although they have been used for decades, their formulations do not mimic the composition of the human cell environment. In this study, we analyzed the impact of a newly formulated human plasma-like media (HPLM) on cell proliferation, mitochondrial bioenergetics, and alterations of drug efficacies using three distinct cancer cell lines. Using high-resolution respirometry, we observed that cells grown in HPLM displayed significantly altered mitochondrial bioenergetic profiles, particularly related to mitochondrial density and mild uncoupling of respiration. Furthermore, in contrast to standard media, the growth of cells in HPLM unveiled mitochondrial dysfunction upon exposure to the FDA-approved kinase inhibitor sunitinib. This seemingly context-dependent side effect of this drug highlights that the selection of the cell culture medium influences the assessment of cancer drug sensitivities. Thus, we suggest to prioritize media with a more physiological composition for analyzing bioenergetic profiles and to take it into account for assigning drug efficacies in the cell culture model of choice.

Mitochondria dysfunction in circulating tumor cells

Circulating tumor cells (CTCs) represent a subset of heterogeneous cells, which, once released from a tumor site, have the potential to give rise to metastasis in secondary sites. Recent research focused on the attempt to detect and characterize these rare cells in the circulation, and advancements in defining their molecular profile have been reported in diverse tumor species, with potential implications for clinical applications. Of note, metabolic alterations, involving mitochondria, have been implicated in the metastatic process, as key determinants in the transition of tumor cells to a mesenchymal or stemness-like phenotype, in drug resistance, and in induction of apoptosis. This review aimed to briefly analyse the most recent knowledge relative to mitochondria dysfunction in CTCs, and to envision implications of altered mitochondria in CTCs for a potential utility in clinics.

Mitochondrial DNA mutations in aging and cancer

Advancing age is a major risk factor for malignant transformation and the development of cancer. As such, over 50% of neoplasms occur in individuals over the age of 70. The pathologies of both ageing and cancer have been characterized by respective groups of molecular hallmarks, and while some features are divergent between the two pathologies, several are shared. Perturbed mitochondrial function is one such common hallmark, and this observation therefore suggests that mitochondrial alterations may be of significance in age‐related cancer development. There is now considerable evidence documenting the accumulation of somatic mitochondrial DNA (mtDNA) mutations in ageing human postmitotic and replicative tissues. Similarly, mutations of the mitochondrial genome have been reported in human cancers for decades. The plethora of functions in which mitochondria partake, such as oxidative phosphorylation, redox balance, apoptosis and numerous biosynthetic pathways, manifests a variety of ways in which alterations in mtDNA may contribute to tumour growth. However, the specific mechanisms by which mtDNA mutations contribute to tumour progression remain elusive and often contradictory. This review aims to consolidate current knowledge and describe future direction within the field.

Hypersensitivity to ferroptosis in chromophobe RCC is mediated by a glutathione metabolic dependency and cystine import via solute carrier family 7 member 11

Chromophobe (Ch) renal cell carcinoma (RCC) arises from the intercalated cell in the distal nephron. There are no proven treatments for metastatic ChRCC. A distinguishing characteristic of ChRCC is strikingly high levels of reduced (GSH) and oxidized (GSSG) glutathione. Here, we demonstrate that ChRCC-derived cells exhibit higher sensitivity to ferroptotic inducers compared with clear-cell RCC. ChRCC-derived cells are critically dependent on cystine via the cystine/glutamate antiporter xCT to maintain high levels of glutathione, making them sensitive to inhibitors of cystine uptake and cyst(e)inase. Gamma-glutamyl transferase 1 (GGT1), a key enzyme in glutathione homeostasis, is markedly suppressed in ChRCC relative to normal kidney. Importantly, GGT1 overexpression inhibits the proliferation of ChRCC cells in vitro and in vivo, suppresses cystine uptake, and decreases levels of GSH and GSSG. Collectively, these data identify ferroptosis as a metabolic vulnerability in ChRCC, providing a potential avenue for targeted therapy for these distinctive tumors.

Decrease of MtDNA copy number affects mitochondrial function and involves in the pathological consequences of ischaemic stroke

The mtDNA copy number can affect the function of mitochondria and play an important role in the development of diseases. However, there are few studies on the mechanism of mtDNA copy number variation and its effects in IS. The specific mechanism of mtDNA copy number variation is still unclear. In this study, mtDNA copy number of 101 IS patients and 101 normal controls were detected by qRT‐PCR, the effect of D‐loop variation on mtDNA copy number of IS patients was explored. Then, a TFAM gene KD‐OE PC12 cell model was constructed to explore the effect of mtDNA copy number variation on mitochondrial function. The results showed that the mtDNA copy number level of the IS group was significantly lower than that of the normal control group (p < 0.05). The relative expression of TFAM gene mRNA in the cells of the OGD/R treatment group was significantly lower than that of the control group (p < 0.05). In addition, after TFAM gene knockdown and over‐expression plasmids were transfected into HEK 293T cells, mtDNA copy number and ATP production level of Sh‐TFAM transfection group was significantly decreased (p < 0.05), while mtDNA copy number and ATP production level of OE‐TFAM transfected group were significantly higher than that of blank control group and OE‐ctrl negative control group (p < 0.01). Our study demonstrated that mitochondrial D‐loop mutation and TFAM gene dysfunction can cause the decrease of mtDNA copy number, thus affecting the mitochondrial metabolism and function of nerve cells, participating in the pathological damage mechanism of IS.

The impact of rare germline variants on human somatic mutation processes

Somatic mutations are an inevitable component of ageing and the most important cause of cancer. The rates and types of somatic mutation vary across individuals, but relatively few inherited influences on mutation processes are known. We perform a gene-based rare variant association study with diverse mutational processes, using human cancer genomes from over 11,000 individuals of European ancestry. By combining burden and variance tests, we identify 207 associations involving 15 somatic mutational phenotypes and 42 genes that replicated in an independent data set at a false discovery rate of 1%. We associate rare inherited deleterious variants in genes such as MSH3, EXO1, SETD2, and MTOR with two phenotypically different forms of DNA mismatch repair deficiency, and variants in genes such as EXO1, PAXIP1, RIF1, and WRN with deficiency in homologous recombination repair. In addition, we identify associations with other mutational processes, such as APEX1 with APOBEC-signature mutagenesis. Many of the genes interact with each other and with known mutator genes within cellular sub-networks. Considered collectively, damaging variants in the identified genes are prevalent in the population. We suggest that rare germline variation in diverse genes commonly impacts mutational processes in somatic cells.

Potential of Mitochondrial Ribosomal Genes as Cancer Biomarkers Demonstrated by Bioinformatics Results

Next-generation sequencing and bioinformatics analyses have clearly revealed the roles of mitochondrial ribosomal genes in cancer development. Mitochondrial ribosomes are composed of three RNA components encoded by mitochondrial DNA and 82 specific protein components encoded by nuclear DNA. They synthesize mitochondrial inner membrane oxidative phosphorylation (OXPHOS)-related proteins and participate in various biological activities via the regulation of energy metabolism and apoptosis. Mitochondrial ribosomal genes are strongly associated with clinical features such as prognosis and foci metastasis in patients with cancer. Accordingly, mitochondrial ribosomes have become an important focus of cancer research. We review recent advances in bioinformatics research that have explored the link between mitochondrial ribosomes and cancer, with a focus on the potential of mitochondrial ribosomal genes as biomarkers in cancer.

ESR1 Regulates the Obesity- and Metabolism-Differential Gene MMAA to Inhibit the Occurrence and Development of Hepatocellular Carcinoma

Obesity is often regarded as a factor that promotes tumorigenesis, but the role of obesity in promoting hepatocellular carcinoma (HCC) is still controversial. We compared the trend change of 14 obesity-related genes in the formation and development of HCC in normal, adjacent, and HCC tissues. Mendelian randomization (MR) analysis was used to verify the relationship between obesity and HCC occurrence. Metabolism of cobalamin-associated A (MMAA) was discovered as an obesity- and metabolism-differential gene, and its function in HCC was tested in vitro and in vivo. Finally, we explored how obese female patients with an originally high expression of female estrogen receptor (ESR1) directly upregulated MMAA to interfere with the progression of HCC. Fourteen obesity-related genes were downregulated in adjacent and tumoral tissues compared with normal liver tissues, which indicated that obesity may be inversely related to the occurrence of HCC and was consistent with the results of MR analysis. We also discovered that MMAA is a metabolic gene closely related to the occurrence and development of HCC by mining the TCGA database, and it functioned an anti-tumor-promoting role in HCC by damaging the mitochondrial function and preserving the redox balance. We further verified that obese females with a high expression of ESR1 can regulate MMAA to protect HCC from progression. This study elucidates that obesity might be a protective factor for female HCC patients, as they originally highly expressed ESR1, which could upregulate MMAA to suppress tumor growth and participate in metabolic reprogramming.

Mutational profiling of mtDNA control region reveals tumor-specific evolutionary selection involved in mitochondrial dysfunction

Background

Mitochondrial DNA (mtDNA) mutations alter mitochondrial function in oxidative metabolism and play an important role in tumorigenesis. A series of studies have demonstrated that the mtDNA control region (mtCTR), which is essential for mtDNA replication and transcription, represents a mutational hotspot in human tumors. However, a comprehensive pan-cancer evolutionary pattern analysis of mtCTR mutations is urgently needed.

Methods

We generated a comprehensive combined dataset containing 10026 mtDNA somatic mutations from 4664 patients, covering 20 tumor types based on public and private next-generation sequencing data.

Findings

Our results demonstrated a significantly higher and much more variable mutation rate in mtCTR than in the coding region across different tumor types. Moreover, our data showed a remarkable distributional bias of tumor somatic mutations between the hypervariable segment (HVS) and non-HVS, with a significantly higher mutation density and average mutation sites in HVS. Importantly, the tumor-specific mutational pattern between mtCTR HVS and non-HVS was identified, which was classified into three evolutionary selection types (relaxed, moderate, and strict constraint types). Analysis of substitution patterns revealed that the prevalence of C_H > T_H in non-HVS greatly contributed to the mutational selection pattern of mtCTR across different tumor types. Furthermore, we found that the mutational pattern of mtCTR in the four tumor types was clearly associated with mitochondrial biogenesis, mitochondrial oxidative metabolism, and the overall survival of patients.

Interpretation

Our results suggest that somatic mutations in mtCTR may be shaped by tumor-specific selective pressure and are involved in tumorigenesis.

Fundings

National Natural Science Foundation of China [grants 82020108023, 81830070, 81872302], and Autonomous Project of State Key Laboratory of Cancer Biology, China [grants CBSKL2019ZZ06, CBSKL2019ZZ27].

Mitochondrial health quality control: measurements and interpretation in the framework of predictive, preventive, and personalized medicine

Mitochondria are the “gatekeeper” in a wide range of cellular functions, signaling events, cell homeostasis, proliferation, and apoptosis. Consequently, mitochondrial injury is linked to systemic effects compromising multi-organ functionality. Although mitochondrial stress is common for many pathomechanisms, individual outcomes differ significantly comprising a spectrum of associated pathologies and their severity grade. Consequently, a highly ambitious task in the paradigm shift from reactive to predictive, preventive, and personalized medicine (PPPM/3PM) is to distinguish between individual disease predisposition and progression under circumstances, resulting in compromised mitochondrial health followed by mitigating measures tailored to the individualized patient profile. For the successful implementation of PPPM concepts, robust parameters are essential to quantify mitochondrial health sustainability. The current article analyses added value of Mitochondrial Health Index (MHI) and Bioenergetic Health Index (BHI) as potential systems to quantify mitochondrial health relevant for the disease development and its severity grade. Based on the pathomechanisms related to the compromised mitochondrial health and in the context of primary, secondary, and tertiary care, a broad spectrum of conditions can significantly benefit from robust quantification systems using MHI/BHI as a prototype to be further improved. Following health conditions can benefit from that: planned pregnancies (improved outcomes for mother and offspring health), suboptimal health conditions with reversible health damage, suboptimal life-style patterns and metabolic syndrome(s) predisposition, multi-factorial stress conditions, genotoxic environment, ischemic stroke of unclear aetiology, phenotypic predisposition to aggressive cancer subtypes, pathologies associated with premature aging and neuro/degeneration, acute infectious diseases such as COVID-19 pandemics, among others.

Mutant KRAS drives metabolic reprogramming and autophagic flux in premalignant pancreatic cells

Mutational activation of the KRAS gene occurs in almost all pancreatic ductal adenocarcinoma (PDAC) and is the earliest molecular event in their carcinogenesis. Evidence has accumulated of the metabolic reprogramming in PDAC, such as amino acid homeostasis and autophagic flux. However, the biological effects of KRAS mutation on metabolic reprogramming at the earlier stages of PDAC carcinogenesis are unclear. Here we report dynamic metabolic reprogramming in immortalized human non-cancerous pancreatic ductal epithelial cells, in which a KRAS mutation was induced by gene-editing, which may mimic early pancreatic carcinogenesis. Similar to the cases of PDAC, KRAS gene mutation increased the dependency on glucose and glutamine for maintaining the intracellular redox balance. In addition, the intracellular levels of amino acids were significantly decreased because of active protein synthesis, and the cells required greater autophagic flux to maintain their viability. The lysosomal inhibitor chloroquine significantly inhibited cell proliferation. Therefore, metabolic reprogramming is an early event in carcinogenesis initiated by KRAS gene mutation, suggesting a rationale for the development of nutritional interventions that suppress or delay the development of PDAC.

Cancer/Testis Antigen 55 is required for cancer cell proliferation and mitochondrial DNA maintenance

Cancer/Testis Antigens (CTAs) represent a group of proteins whose expression under physiological conditions is restricted to testis but activated in many human cancers. Also, it was observed that co-expression of multiple CTAs worsens the patient prognosis. Five CTAs were reported acting in mitochondria and we recently reported 147 transcripts encoded by 67 CTAs encoding for proteins potentially targeted to mitochondria. Among them, we identified the two isoforms encoded by CT55 for whom the function is poorly understood. First, we found that patients with tumors expressing wild-type CT55 are associated with poor survival. Moreover, CT55 silencing decreases dramatically cell proliferation. Second, to investigate the role of CT55 on mitochondria, we first show that CT55 is localized to both mitochondria and endoplasmic reticulum (ER) due to the presence of an ambiguous N-terminal targeting signal. Then, we show that CT55 silencing decreases mtDNA copy number and delays mtDNA recovery after an acute depletion. Moreover, demethylation of CT55 promotor increases its expression, which in turn increases mtDNA copy number. Finally, we measured the mtDNA copy number in NCI-60 cell lines and screened for genes whose expression is strongly correlated to mtDNA amount. We identified CT55 as the second highest correlated hit. Also, we show that compared to siRNA scrambled control (siCtrl) treatment, CT55 specific siRNA (siCT55) treatment down-regulates aerobic respiration, indicating that CT55 sustains mitochondrial respiration. Altogether, these data show for first time that CT55 acts on mtDNA copy number, modulates mitochondrial activity to sustain cancer cell proliferation.

Single cell transcriptomics reveals unique metabolic profiles of ependymoma subgroups

OBJECTIVE

Ependymomas are biologically diverse tumors with five major genomic subgroups. However, intratumor heterogeneity continues to be poorly understood. The present study characterized the metabolic landscapes of ependymoma subgroups at the single-cell level.

METHODS

Expression profiles from 11,200 ependymoma single cells derived from the five major subgroups and 7,200 ependymoma-derived non-neoplastic cells were computationally analyzed using a robust workflow to elucidate relative differences in metabolic pathway activities.

RESULTS

Dimensionality reduction using metabolic expression profiles exhibited clustering corresponding to each tumor subgroup, but non-neoplastic cells exhibited no discernable differences between subgroups. From the 80 metabolic pathways examined, over 75 pathways had significantly different activity scores between ependymoma subgroups. Further analysis of metabolic heterogeneity suggests that mitochondrial oxidative phosphorylation accounts for considerable metabolic variation within tumor subgroups and non-neoplastic cells of the same cell type. Drug metabolism pathways, specifically those involving cytochromes P450, were also found to be major contributors to heterogeneity.

CONCLUSIONS

Ependymoma subgroups display distinct metabolic differences as evaluated through gene expression profiles with certain pathways contributing greatly to intra-subgroup variation. These results may account for variation in tumor metabolism, treatment response, and potential targeting approaches that disrupt metabolic signalling.

The Potential Role of Mitochondrial Acetaldehyde Dehydrogenase 2 in Urological Cancers From the Perspective of Ferroptosis and Cellular Senescence

Overproduction of reactive oxygen species (ROS) and superlative lipid peroxidation promote tumorigenesis, and mitochondrial aldehyde dehydrogenase 2 (ALDH2) is associated with the detoxification of ROS-mediated lipid peroxidation-generated reactive aldehydes such as 4-hydroxy-2-nonenal (4-HNE), malondialdehyde, and acrolein due to tobacco smoking. ALDH2 has been demonstrated to be highly associated with the prognosis and chemoradiotherapy sensitivity of many types of cancer, including leukemia, lung cancer, head and neck cancer, esophageal cancer, hepatocellular cancer, pancreatic cancer, and ovarian cancer. In this study, we explored the possible relationship between ALDH2 and urological cancers from the aspects of ferroptosis, epigenetic alterations, proteostasis, mitochondrial dysfunction, and cellular senescence.

Exogenous Glutathione Protects IPEC-J2 Cells against Oxidative Stress through a Mitochondrial Mechanism

The accumulation of reactive oxygen species (ROS) triggers oxidative stress in cells by oxidizing and modifying various cellular components, preventing them from performing their inherent functions, ultimately leading to apoptosis and autophagy. Glutathione (GSH) is a ubiquitous intracellular peptide with multiple functions. In this study, a hydrogen peroxide (H2O2)-induced oxidative damage model in IPEC-J2 cells was used to investigate the cellular protection mechanism of exogenous GSH against oxidative stress. The results showed that GSH supplement improved the cell viability reduced by H2O2-induced oxidative damage model in IPEC-J2 cells in a dose-dependent manner. Moreover, supplement with GSH also attenuated the H2O2-induced MMP loss, and effectively decreased the H2O2-induced mitochondrial dysfunction by increasing the content of mtDNA and upregulating the expression TFAM. Exogenous GSH treatment significantly decreased the ROS and MDA levels, improved SOD activity in H2O2-treated cells and reduced H2O2-induced early apoptosis in IPEC-J2 cells. This study showed that exogenous GSH can protect IPEC-J2 cells against apoptosis induced by oxidative stress through mitochondrial mechanisms.

PGC1α/β Expression Predicts Therapeutic Response to Oxidative Phosphorylation Inhibition in Ovarian Cancer

Ovarian cancer is the deadliest gynecologic cancer, and novel therapeutic options are crucial to improve overall survival. Here we provide evidence that impairment of oxidative phosphorylation (OXPHOS) can help control ovarian cancer progression, and this benefit correlates with expression of the two mitochondrial master regulators PGC1α and PGC1β. In orthotopic patient-derived ovarian cancer xenografts (OC-PDX), concomitant high expression of PGC1α and PGC1β (PGC1α/β) fostered a unique transcriptional signature, leading to increased mitochondrial abundance, enhanced tricarboxylic acid cycling, and elevated cellular respiration that ultimately conferred vulnerability to OXPHOS inhibition. Treatment with the respiratory chain complex I inhibitor IACS-010759 caused mitochondrial swelling and ATP depletion that consequently delayed malignant progression and prolonged the lifespan of high PGC1α/β-expressing OC-PDX-bearing mice. Conversely, low PGC1α/β OC-PDXs were not affected by IACS-010759, thus pinpointing a selective antitumor effect of OXPHOS inhibition. The clinical relevance of these findings was substantiated by analysis of ovarian cancer patient datasets, which showed that 25% of all cases displayed high PGC1α/β expression along with an activated mitochondrial gene program. This study endorses the use of OXPHOS inhibitors to manage ovarian cancer and identifies the high expression of both PGC1α and β as biomarkers to refine the selection of patients likely to benefit most from this therapy.

SIGNIFICANCE

OXPHOS inhibition in ovarian cancer can exploit the metabolic vulnerabilities conferred by high PGC1α/β expression and offers an effective approach to manage patients on the basis of PGC1α/β expression.

Bioenergetic and Metabolic Adaptation in Tumor Progression and Metastasis

The ability of cancer cells to adjust their metabolism in response to environmental changes is a well-recognized hallmark of cancer. Diverse cancer and non-cancer cells within tumors compete for metabolic resources. Metabolic demands change frequently during tumor initiation, progression and metastasis, challenging our quest to better understand tumor biology and develop novel therapeutics. Vascularization, physical constraints, immune responses and genetic instability promote tumor evolution resulting in immune evasion, opportunities to breach basement membrane barriers and spread through the circulation and lymphatics. In addition, the unfolded protein response linked to the ubiquitin proteasome system is a key player in addressing stoichiometric imbalances between nuclear and mitochondrially-encoded protein subunits of respiratory complexes, and nuclear-encoded mitochondrial ribosomal protein subunits. While progressive genetic changes, some of which affect metabolic adaptability, contribute to tumorigenesis and metastasis through clonal expansion, epigenetic changes are also important and more dynamic in nature. Understanding the role of stromal and immune cells in the tumor microenvironment in remodeling cancer cell energy metabolism has become an increasingly important area of research. In this perspective, we discuss the adaptations made by cancer cells to balance mitochondrial and glycolytic energy metabolism. We discuss how hypoxia and nutrient limitations affect reductive and oxidative stress through changes in mitochondrial electron transport activity. We propose that integrated responses to cellular stress in cancer cells are central to metabolic flexibility in general and bioenergetic adaptability in particular and are paramount in tumor progression and metastasis.

Chemotherapy Resistance: Role of Mitochondrial and Autophagic Components

Simple Summary

Chemotherapy resistance is a common occurrence during cancer treatment that cancer researchers are attempting to understand and overcome. Mitochondria are a crucial intracellular signaling core that are becoming important determinants of numerous aspects of cancer genesis and progression, such as metabolic reprogramming, metastatic capability, and chemotherapeutic resistance. Mitophagy, or selective autophagy of mitochondria, can influence both the efficacy of tumor chemotherapy and the degree of drug resistance. Regardless of the fact that mitochondria are well-known for coordinating ATP synthesis from cellular respiration in cellular bioenergetics, little is known its mitophagy regulation in chemoresistance. Recent advancements in mitochondrial research, mitophagy regulatory mechanisms, and their implications for our understanding of chemotherapy resistance are discussed in this review.

Abstract

Cancer chemotherapy resistance is one of the most critical obstacles in cancer therapy. One of the well-known mechanisms of chemotherapy resistance is the change in the mitochondrial death pathways which occur when cells are under stressful situations, such as chemotherapy. Mitophagy, or mitochondrial selective autophagy, is critical for cell quality control because it can efficiently break down, remove, and recycle defective or damaged mitochondria. As cancer cells use mitophagy to rapidly sweep away damaged mitochondria in order to mediate their own drug resistance, it influences the efficacy of tumor chemotherapy as well as the degree of drug resistance. Yet despite the importance of mitochondria and mitophagy in chemotherapy resistance, little is known about the precise mechanisms involved. As a consequence, identifying potential therapeutic targets by analyzing the signal pathways that govern mitophagy has become a vital research goal. In this paper, we review recent advances in mitochondrial research, mitophagy control mechanisms, and their implications for our understanding of chemotherapy resistance.

Ultra-deep whole genome bisulfite sequencing reveals a single methylation hotspot in human brain mitochondrial DNA

While DNA methylation is established as a major regulator of gene expression in the nucleus, the existence of mitochondrial DNA (mtDNA) methylation remains controversial. Here, we characterized the mtDNA methylation landscape in the prefrontal cortex of neurological healthy individuals (n=26) and patients with Parkinson’s disease (n=27), using a combination of whole-genome bisulphite sequencing (WGBS) and bisulphite-independent methods. Accurate mtDNA mapping from WGBS data required alignment to an mtDNA reference only, to avoid misalignment to nuclear mitochondrial pseudogenes. Once correctly aligned, WGBS data provided ultra-deep mtDNA coverage (16,723 ± 7,711) and revealed overall very low levels of cytosine methylation. The highest methylation levels (5.49 ± 0.97%) were found on CpG position m.545, located in the heavy-strand promoter 1 region. The m.545 methylation was validated using a combination of methylation-sensitive DNA digestion and quantitative PCR analysis. We detected no association between mtDNA methylation profile and Parkinson’s disease. Interestingly, m.545 methylation correlated with the levels of mtDNA transcripts, suggesting a putative role in regulating mtDNA gene expression. In addition, we propose a robust framework for methylation analysis of mtDNA from WGBS data, which is less prone to false-positive findings due to misalignment of nuclear mitochondrial pseudogene sequences.

Role of Mitochondrial Stress Response in Cancer Progression

Mitochondria are subcellular organelles that are a hub for key biological processes, such as bioenergetic, biosynthetic, and signaling functions. Mitochondria are implicated in all oncogenic processes, from malignant transformation to metastasis and resistance to chemotherapeutics. The harsh tumor environment constantly exposes cancer cells to cytotoxic stressors, such as nutrient starvation, low oxygen, and oxidative stress. Excessive or prolonged exposure to these stressors can cause irreversible mitochondrial damage, leading to cell death. To survive hostile microenvironments that perturb mitochondrial function, cancer cells activate a stress response to maintain mitochondrial protein and genome integrity. This adaptive mechanism, which is closely linked to mitochondrial function, enables rapid adjustment and survival in harsh environmental conditions encountered during tumor dissemination, thereby promoting cancer progression. In this review, we describe how the mitochondria stress response contributes to the acquisition of typical malignant traits and highlight the potential of targeting the mitochondrial stress response as an anti-cancer therapeutic strategy.

Synergistic anti-tumor efficacy of mutant isocitrate dehydrogenase 1 inhibitor SYC-435 with standard therapy in patient-derived xenograft mouse models of glioma

•

A novel pair of orthotopic PDX models of glioma bearing IDH1-R132H/R132C mutations.

•

New mutant IDH1i (SY-435) with standard therapy led to strong therapeutic efficacy.

•

H3K4/K9 methylation/mtDNA-encoded molecules mediate anti-tumor activity of SYC-435.

•

Discovered MYO1F, CTC1 and BCL9 as novel genes that mediated SYC-435 resistance.

Clinical outcomes in patients with WHO grade II/III astrocytoma, oligodendroglioma or secondary glioblastoma remain poor. Isocitrate dehydrogenase 1 (IDH1) is mutated in > 70% of these tumors, making it an attractive therapeutic target. To determine the efficacy of our newly developed mutant IDH1 inhibitor, SYC-435 (1-hydroxypyridin-2-one), we treated orthotopic glioma xenograft model (IC-BT142AOA) carrying R132H mutation and our newly established orthotopic patient-derived xenograft (PDX) model of recurrent anaplastic oligoastrocytoma (IC-V0914AOA) bearing R132C mutation. In addition to suppressing IDH1 mutant cell proliferation in vitro, SYC-435 (15 mg/kg, daily x 28 days) synergistically prolonged animal survival times with standard therapies (Temozolomide + fractionated radiation) mediated by reduction of H3K4/H3K9 methylation and expression of mitochondrial DNA (mtDNA)-encoded molecules. Furthermore, RNA-seq of the remnant tumors identified genes (MYO1F, CTC1 and BCL9) and pathways (base excision repair, TCA cycle II, sirtuin signaling, protein kinase A, eukaryotic initiation factor 2 and α-adrenergic signaling) as mediators of therapy resistance. Our data demonstrated the efficacy SYC-435 in targeting IDH1 mutant gliomas when combined with standard therapy and identified a novel set of genes that should be prioritized for future studies to overcome SYC-435 resistance.

Single-cell mtDNA heteroplasmy in colorectal cancer

Human mitochondria can be genetically distinct within the same individual, a phenomenon known as heteroplasmy. In cancer, this phenomenon seems exacerbated, and most mitochondrial mutations seem to be heteroplasmic. How this genetic variation is arranged within and among normal and tumor cells is not well understood. To address this question, here we sequenced single-cell mitochondrial genomes from multiple normal and tumoral locations in four colorectal cancer patients. Our results suggest that single cells, both normal and tumoral, can carry various mitochondrial haplotypes. Remarkably, this intra-cell heteroplasmy can arise before tumor development and be maintained afterward in specific tumoral cell subpopulations. At least in the colorectal patients studied here, the somatic mutations in the single-cells do not seem to have a prominent role in tumorigenesis.

Mutations in TP53 and PIK3CA genes in hepatocellular carcinoma patients are associated with chronic Schistosomiasis

The purpose of this study was to evaluate the genetic variation of the PIK3CA gene and the histopathological changes in liver tissue of patients with chronic Schistosomiasis to predict hepatocellular carcinoma. In this retrospective, the study samples were taken from 20 patients, divided into chronic schistosomiasis infected group of people (S) and chronic schistosomiasis uninfected group of people (C). The liver tissue biopsy samples for histological examinations were obtained only from chronic Schistosomiasis patients (n = 9). The blood samples were obtained from groups S and C for the mutational analysis of the PIK3CA and TP53 genes. The results suggest that the patients diagnosed with chronic Schistosomiasis were 9 (55%), and healthy patients without Schistosomiasis were 11 (45%). Histological results found that proliferation of fibrosis was observed in the hepatocytes of schistosomiasis patients. A total of 8 mutations (5 male, 3 female) were detected in PIK3CA and TP53 genes. Including 1634 A > G substitution mutations in PIK3CA, which was the only mutation found in males and females among the 8 mutations, accounting 22.22%. PIK3CA gene mutations were found more predominant in male groups as compared to other TP53 gene mutations. In conclusion, this study found that patients with chronic Schistosomiasis are at risk of PIK3CA gene mutations, eventually leading to hepatocytes fibrosis and liver cancer.

Arsenic exposure and biomarkers for oxidative stress and telomere length in indigenous populations in Bolivia

BACKGROUND

Women living in the Bolivian Andes are environmentally exposed to arsenic, yet there is scarce information about arsenic-related effects in this region. Several biomarkers for telomere length and oxidative stress (mitochondrial DNA copy number, mtDNAcn; 8-Oxo-2'-deoxyguanosine, 8-oxo-dG; and 4-hydroxy nonenal mercapturic acid, 4-HNE-MA) have been previously linked to arsenic, and some of which are prospective biomarkers for cancer risk.

OBJECTIVE AND HYPOTHESIS

To evaluate associations between arsenic exposure and telomere length, mtDNAcn, 8-oxo-dG, and 4-HNE-MA in Bolivians. Arsenic exposure was hypothesized to be positively associated with all four toxicity biomarkers, particularly in individuals with a less efficient arsenic metabolism.

METHODS

The study encompassed 193 indigenous women. Arsenic exposure was assessed in urine as the sum of inorganic arsenic metabolite concentrations (U-As) measured by HPLC-HG-ICP-MS, and in whole blood as total arsenic (B-As) measured by ICP-MS. Efficiency of arsenic metabolism was evaluated by a polymorphism (rs3740393) in the main arsenic methylating gene AS3MT measured by TaqMan allelic discrimination, and by the relative fractions of urinary inorganic arsenic metabolites. Telomere length and mtDNAcn were determined in peripheral blood leukocytes by quantitative PCR, and urinary 8-oxo-dG and 4-HNE-MA by LC-MS/MS.

RESULTS

U-As and B-As were associated with longer telomeres and higher mtDNAcn, particularly in women with a less efficient arsenic metabolism. Urinary 8-oxo-dG and 4-HNE-MA were positively associated with U-As, but only 4-HNE-MA was associated with B-As. Arsenic metabolism efficiency did not have a clear effect on the concentrations of either of these biomarkers.

CONCLUSION

Bolivian women showed indications of arsenic toxicity, measured by four different biomarkers. Telomere length, mtDNAcn, and 4-HNE-MA were positively associated with both U-As and B-As. The association of arsenic exposure with telomere length and mtDNAcn was only present in Bolivian women with a less efficient metabolism. These findings call for additional efforts to evaluate and reduce arsenic exposure in Bolivia.

Discovery and Lead Optimization of Benzene-1,4-disulfonamides as Oxidative Phosphorylation Inhibitors

Inhibition of oxidative phosphorylation (OXPHOS) is a promising therapeutic strategy for select cancers that are dependent on aerobic metabolism. Here, we report the discovery, optimization, and structure-activity relationship (SAR) study of a series of novel OXPHOS inhibitors. The hit compound, benzene-1,4-disulfonamide 1, was discovered in a phenotypic screen selective for cytotoxicity in a galactose-containing medium. Our multi-parameter optimization campaign led to the discovery of 65 (DX3-235), showing nanomolar inhibition of complex I function and adenosine triphosphate (ATP) production in a galactose-containing medium resulting in significant cytotoxicity. Importantly, 64 (DX3-234), a close analogue of 65, is well tolerated in mice and shows significant single agent efficacy in a Pan02 syngeneic pancreatic cancer model, suggesting that highly potent and selective OXPHOS inhibitors can be useful for the treatment of pancreatic cancer.

Genetic Alterations in Mitochondrial DNA Are Complementary to Nuclear DNA Mutations in Pheochromocytomas

Simple Summary

Mitochondrial DNA (mtDNA) alterations have been reported to play important roles in cancer development and metastasis. However, there is scarce information about pheochromocytomas and paragangliomas (PCCs/PGLs) formation. To determine the potential roles of mtDNA alterations in PCCs/PGLs, we analyzed a panel of 26 nuclear susceptibility genes and the entire mtDNA sequence of 77 human tumors, using NGS. We also performed an analysis of copy-number alterations, large mtDNA deletion, and gene/protein expression. Our results revealed that 53.2% of the tumors harbor a mutation in the susceptibility genes and 16.9% harbor complementary mitochondrial mutations. Large deletions and depletion of mtDNA were found in 26% and 87% of tumors, respectively, accompanied by a reduced expression of the mitochondrial biogenesis markers (PCG1α, NRF1, and TFAM). Furthermore, P62 and LC3a gene expression suggested increased mitophagy, which is linked to mitochondrial dysfunction. These finding suggest a complementarity and a potential contributing role in PCCs/PGLs tumorigenesis.

Abstract

Background: Somatic mutations, copy-number variations, and genome instability of mitochondrial DNA (mtDNA) have been reported in different types of cancers and are suggested to play important roles in cancer development and metastasis. However, there is scarce information about pheochromocytomas and paragangliomas (PCCs/PGLs) formation. Material: To determine the potential roles of mtDNA alterations in sporadic PCCs/PGLs, we analyzed a panel of 26 nuclear susceptibility genes and the entire mtDNA sequence of seventy-seven human tumors, using next-generation sequencing, and compared the results with normal adrenal medulla tissues. We also performed an analysis of copy-number alterations, large mtDNA deletion, and gene and protein expression. Results: Our results revealed that 53.2% of the tumors harbor a mutation in at least one of the targeted susceptibility genes, and 16.9% harbor complementary mitochondrial mutations. More than 50% of the mitochondrial mutations were novel and predicted pathogenic, affecting mitochondrial oxidative phosphorylation. Large deletions were found in 26% of tumors, and depletion of mtDNA occurred in more than 87% of PCCs/PGLs. The reduction of the mitochondrial number was accompanied by a reduced expression of the regulators that promote mitochondrial biogenesis (PCG1α, NRF1, and TFAM). Further, P62 and LC3a gene expression suggested increased mitophagy, which is linked to mitochondrial dysfunction. Conclusion: The pathogenic role of these finding remains to be shown, but we suggest a complementarity and a potential contributing role in PCCs/PGLs tumorigenesis.

Metabolic resistance to the inhibition of mitochondrial transcription revealed by CRISPR-Cas9 screen

Cancer cells depend on mitochondria to sustain their increased metabolic need and mitochondria therefore constitute possible targets for cancer treatment. We recently developed small-molecule inhibitors of mitochondrial transcription (IMTs) that selectively impair mitochondrial gene expression. IMTs have potent antitumor properties in vitro and in vivo, without affecting normal tissues. Because therapy-induced resistance is a major constraint to successful cancer therapy, we investigated mechanisms conferring resistance to IMTs. We employed a CRISPR-Cas9 (clustered regularly interspaced short palindromic repeats)-(CRISP-associated protein 9) whole-genome screen to determine pathways conferring resistance to acute IMT1 treatment. Loss of genes belonging to von Hippel-Lindau (VHL) and mammalian target of rapamycin complex 1 (mTORC1) pathways caused resistance to acute IMT1 treatment and the relevance of these pathways was confirmed by chemical modulation. We also generated cells resistant to chronic IMT treatment to understand responses to persistent mitochondrial gene expression impairment. We report that IMT1-acquired resistance occurs through a compensatory increase of mitochondrial DNA (mtDNA) expression and cellular metabolites. We found that mitochondrial transcription factor A (TFAM) downregulation and inhibition of mitochondrial translation impaired survival of resistant cells. The identified susceptibility and resistance mechanisms to IMTs may be relevant for different types of mitochondria-targeted therapies.

Genome-wide analysis of mitochondrial DNA copy number reveals loci implicated in nucleotide metabolism, platelet activation, and megakaryocyte proliferation

Mitochondrial DNA copy number (mtDNA-CN) measured from blood specimens is a minimally invasive marker of mitochondrial function that exhibits both inter-individual and intercellular variation. To identify genes involved in regulating mitochondrial function, we performed a genome-wide association study (GWAS) in 465,809 White individuals from the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) consortium and the UK Biobank (UKB). We identified 133 SNPs with statistically significant, independent effects associated with mtDNA-CN across 100 loci. A combination of fine-mapping, variant annotation, and co-localization analyses was used to prioritize genes within each of the 133 independent sites. Putative causal genes were enriched for known mitochondrial DNA depletion syndromes (p = 3.09 × 10-15) and the gene ontology (GO) terms for mtDNA metabolism (p = 1.43 × 10-8) and mtDNA replication (p = 1.2 × 10-7). A clustering approach leveraged pleiotropy between mtDNA-CN associated SNPs and 41 mtDNA-CN associated phenotypes to identify functional domains, revealing three distinct groups, including platelet activation, megakaryocyte proliferation, and mtDNA metabolism. Finally, using mitochondrial SNPs, we establish causal relationships between mitochondrial function and a variety of blood cell-related traits, kidney function, liver function and overall (p = 0.044) and non-cancer mortality (p = 6.56 × 10-4).

Epigenome-wide association study of mitochondrial genome copy number

We conducted cohort- and race-specific epigenome-wide association analyses of mitochondrial deoxyribonucleic acid (mtDNA) copy number (mtDNA CN) measured in whole blood from participants of African and European origins in five cohorts (n = 6182, mean age = 57-67 years, 65% women). In the meta-analysis of all the participants, we discovered 21 mtDNA CN-associated DNA methylation sites (CpG) (P < 1 × 10-7), with a 0.7-3.0 standard deviation increase (3 CpGs) or decrease (18 CpGs) in mtDNA CN corresponding to a 1% increase in DNA methylation. Several significant CpGs have been reported to be associated with at least two risk factors (e.g. chronological age or smoking) for cardiovascular disease (CVD). Five genes [PR/SET domain 16, nuclear receptor subfamily 1 group H member 3 (NR1H3), DNA repair protein, DNA polymerase kappa and decaprenyl-diphosphate synthase subunit 2], which harbor nine significant CpGs, are known to be involved in mitochondrial biosynthesis and functions. For example, NR1H3 encodes a transcription factor that is differentially expressed during an adipose tissue transition. The methylation level of cg09548275 in NR1H3 was negatively associated with mtDNA CN (effect size = -1.71, P = 4 × 10-8) and was positively associated with the NR1H3 expression level (effect size = 0.43, P = 0.0003), which indicates that the methylation level in NR1H3 may underlie the relationship between mtDNA CN, the NR1H3 transcription factor and energy expenditure. In summary, the study results suggest that mtDNA CN variation in whole blood is associated with DNA methylation levels in genes that are involved in a wide range of mitochondrial activities. These findings will help reveal molecular mechanisms between mtDNA CN and CVD.

Monitoring and Modulating mtDNA G-Quadruplex Dynamics Reveal Its Close Relationship to Cell Glycolysis

The mitochondrial DNA G-quadruplex (mtDNA G4) is a potential regulatory element for the regulation of mitochondrial functions; however, its relevance and specific roles in diseases remain largely unknown. Here, we engineered a set of chemical probes, including MitoISCH, an mtDNA G4-specific fluorescent probe, together with MitoPDS, a mitochondria-targeted G4-stabilizing agent, to thoroughly investigate mtDNA G4s. Using MitoISCH to monitor previously intractable dynamics of mtDNA G4s, we surprisingly found that their formation was prevalent only in endothelial and cancer cells that rely on glycolysis for energy production. Consistent with this, promotion of mtDNA G4 folding by MitoPDS in turn caused glycolysis-related gene activation and glycolysis enhancement. Remarkably, this close relationship among mtDNA G4s, glycolysis, and cancer cells further allowed MitoISCH to accumulate in tumors and label them in vivo. Our work reveals an unprecedented link between mtDNA G4s and cell glycolysis, suggesting that mtDNA G4s may be a novel cancer biomarker and therapeutic target deserving further exploration.

A mitochondrial response to oxidative stress mediated by unscheduled RNA-DNA hybrids (R-loops)

How oxidative stress promotes aging-related human diseases like cancer and neurodegeneration remains unclear. Here, we discuss the origins and implications of an oxidative-stress response recently reported to destabilize the mitochondrial (mt) genome via unscheduled RNA/DNA hybrid (R-loop) accumulation, by impairing the recruitment of RNAseH1 to the regulatory regions of mtDNA.

Malignancy and NF-κB signalling strengthen coordination between expression of mitochondrial and nuclear-encoded oxidative phosphorylation genes

BACKGROUND

Mitochondria are ancient endosymbiotic organelles crucial to eukaryotic growth and metabolism. The mammalian mitochondrial genome encodes for 13 mitochondrial proteins, and the remaining mitochondrial proteins are encoded by the nuclear genome. Little is known about how coordination between the expression of the two sets of genes is achieved.

RESULTS

Correlation analysis of RNA-seq expression data from large publicly available datasets is a common method to leverage genetic diversity to infer gene co-expression modules. Here we use this method to investigate nuclear-mitochondrial gene expression coordination. We identify a pitfall in correlation analysis that results from the large variation in the proportion of transcripts from the mitochondrial genome in RNA-seq data. Commonly used normalisation techniques based on total read counts, such as FPKM or TPM, produce artefactual negative correlations between mitochondrial- and nuclear-encoded transcripts. This also results in artefactual correlations between pairs of nuclear-encoded genes, with important consequences for inferring co-expression modules beyond mitochondria. We show that these effects can be overcome by normalizing using the median-ratio normalisation (MRN) or trimmed mean of M values (TMM) methods. Using these normalisations, we find only weak and inconsistent correlations between mitochondrial and nuclear-encoded mitochondrial genes in the majority of healthy human tissues from the GTEx database.

CONCLUSIONS

We show that a subset of healthy tissues with high expression of NF-κB show significant coordination, suggesting a role for NF-κB in ensuring balanced expression between mitochondrial and nuclear genes. Contrastingly, most cancer types show robust coordination of nuclear and mitochondrial OXPHOS gene expression, identifying this as a feature of gene regulation in cancer.

Hypoxia Regulates Endogenous Double-Stranded RNA Production via Reduced Mitochondrial DNA Transcription

Hypoxia is a common phenomenon in solid tumours strongly linked to the hallmarks of cancer. Hypoxia promotes local immunosuppression and downregulates type I interferon (IFN) expression and signalling, which contribute to the success of many cancer therapies. Double-stranded RNA (dsRNA), transiently generated during mitochondrial transcription, endogenously activates the type I IFN pathway. We report the effects of hypoxia on the generation of mitochondrial dsRNA (mtdsRNA) in breast cancer. We found a significant decrease in dsRNA production in different cell lines under hypoxia. This effect was HIF1α/2α-independent. mtdsRNA was responsible for induction of type I IFN and significantly decreased after hypoxia. Mitochondrially encoded gene expression was downregulated and mtdsRNA bound by the dsRNA-specific J2 antibody was decreased during hypoxia. These findings reveal a new mechanism of hypoxia-induced immunosuppression that could be targeted by hypoxia-activated therapies.

Micropeptide ASAP encoded by LINC00467 promotes colorectal cancer progression by directly modulating ATP synthase activity

Emerging evidence has shown that open reading frames inside long noncoding RNAs (lncRNAs) could encode micropeptides. However, their roles in cellular energy metabolism and tumor progression remain largely unknown. Here, we identified a 94 amino acid-length micropeptide encoded by lncRNA LINC00467 in colorectal cancer. We also characterized its conservation across higher mammals, localization to mitochondria, and the concerted local functions. This peptide enhanced the ATP synthase construction by interacting with the subunits α and γ (ATP5A and ATP5C), increased ATP synthase activity and mitochondrial oxygen consumption rate, and thereby promoted colorectal cancer cell proliferation. Hence, this micropeptide was termed ATP synthase-associated peptide (ASAP). Furthermore, loss of ASAP suppressed patient-derived xenograft growth with attenuated ATP synthase activity and mitochondrial ATP production. Clinically, high expression of ASAP and LINC00467 predicted poor prognosis of colorectal cancer patients. Taken together, our findings revealed a colorectal cancer-associated micropeptide as a vital player in mitochondrial metabolism and provided a therapeutic target for colorectal cancer.

A Novel Mitochondrial-Related Nuclear Gene Signature Predicts Overall Survival of Lung Adenocarcinoma Patients

Background: Lung cancer is the leading cause of cancer-related death worldwide, of which lung adenocarcinoma (LUAD) is one of the main histological subtypes. Mitochondria are vital for maintaining the physiological function, and their dysfunction has been found to be correlated with tumorigenesis and disease progression. Although, some mitochondrial-related genes have been found to correlate with the clinical outcomes of multiple tumors solely. The integrated relationship between nuclear mitochondrial genes (NMGs) and the prognosis of LUAD remains unclear.

Methods: The list of NMGs, gene expression data, and related clinical information of LUAD were downloaded from public databases. Bioinformatics methods were used and obtained 18 prognostic related NMGs to construct a risk signature.

Results: There were 18 NMGs (NDUFS2, ATP8A2, SCO1, COX14, COA6, RRM2B, TFAM, DARS2, GARS, YARS2, EFG1, GFM1, MRPL3, MRPL44, ISCU, CABC1, HSPD1, and ETHE1) identified by LASSO regression analysis. The mRNA expression of these 18 genes was positively correlated with their relative linear copy number alteration (CNA). Meanwhile, the established risk signature could effectively distinguish high- and low-risk patients, and its predictive capacity was validated in three independent gene expression omnibus (GEO) cohorts. Notably, a significantly lower prevalence of actionable EGFR alterations was presented in patients with high-risk NMGs signature but accompanied with a more inflame immune tumor microenvironment. Additionally, multicomponent Cox regression analysis showed that the model was stable when risk score, tumor stage, and lymph node stage were considered, and the 1-, 3-, and 5-year AUC were 0.74, 0.75, and 0.70, respectively.

Conclusion: Together, this study established a signature based on NMGs that is a prognostic biomarker for LUAD patients and has the potential to be widely applied in future clinical settings.

Baicalein resensitizes tamoxifen-resistant breast cancer cells by reducing aerobic glycolysis and reversing mitochondrial dysfunction via inhibition of hypoxia-inducible factor-1α

Drug resistance is a major hurdle for the effectiveness of tamoxifen (TAM) to provide clinical benefit. Therefore, it is essential to identify a sensitizer that could be used to improve TAM efficacy in treating TAM-resistant breast cancer. Here, we investigated the ability of baicalein to reverse TAM resistance. We found that baicalein increased the efficacy of TAM in inhibiting proliferation and inducing apoptosis of TAM-resistant cells. It also enhanced the TAM-induced growth reduction of resistant cells from NOD/SCID mouse mammary fat pads, without causing obvious systemic toxicity. Analyses using the CellMiner tool and the Kaplan-Meier plotter database showed that HIF-1α expression was inversely correlated with TAM therapeutic response in NCI-60 cancer cells and breast cancer patients. HIF-1α expression was increased in TAM-resistant cells due to an increase in mRNA levels and reduced ubiquitin-mediated degradation. Baicalein reduced HIF-1α expression by promoting its interaction with PHD2 and pVHL, thus facilitating ubiquitin ligase-mediated proteasomal degradation and thereby suppressing the nuclear translocation, binding to the hypoxia-response element, and transcriptional activity of HIF-1α. As a result, baicalein downregulated aerobic glycolysis by restricting glucose uptake, lactate production, ATP generation, lactate/pyruvate ratio and expression of HIF-1α-targeted glycolytic genes, thereby enhancing the antiproliferative efficacy of TAM. Furthermore, baicalein interfered with HIF-1α inhibition of mitochondrial biosynthesis, which increased mitochondrial DNA content and mitochondrial numbers, restored the generation of reactive oxygen species in mitochondria, and thus enhanced the TAM-induced mitochondrial apoptotic pathway. The HIF-1α stabilizer dimethyloxallyl glycine prevented the baicalein-induced downregulation of glycolysis and mitochondrial biosynthesis and reduced the effects of baicalein on reversing TAM resistance. Our results indicate that baicalein is a promising candidate to help overcome TAM resistance by sensitizing resistant cells to TAM-induced growth inhibition and apoptosis. The mechanism underlying the effects of baicalein consists of inhibition of HIF-1α-mediated aerobic glycolysis and mitochondrial dysfunction.

Beyond base excision repair: an evolving picture of mitochondrial DNA repair

Mitochondria are highly specialised organelles required for key cellular processes including ATP production through cellular respiration and controlling cell death via apoptosis. Unlike other organelles, mitochondria contain their own DNA genome which encodes both protein and RNA required for cellular respiration. Each cell may contain hundreds to thousands of copies of the mitochondrial genome, which is essential for normal cellular function - deviation of mitochondrial DNA (mtDNA) copy number is associated with cellular ageing and disease. Furthermore, mtDNA lesions can arise from both endogenous or exogenous sources and must either be tolerated or corrected to preserve mitochondrial function. Importantly, replication of damaged mtDNA can lead to stalling and introduction of mutations or genetic loss, mitochondria have adapted mechanisms to repair damaged DNA. These mechanisms rely on nuclear-encoded DNA repair proteins that are translocated into the mitochondria. Despite the presence of many known nuclear DNA repair proteins being found in the mitochondrial proteome, it remains to be established which DNA repair mechanisms are functional in mammalian mitochondria. Here, we summarise the existing and emerging research, alongside examining proteomic evidence, demonstrating that mtDNA damage can be repaired using Base Excision Repair (BER), Homologous Recombination (HR) and Microhomology-mediated End Joining (MMEJ). Critically, these repair mechanisms do not operate in isolation and evidence for interplay between pathways and repair associated with replication is discussed. Importantly, characterising non-canonical functions of key proteins and understanding the bespoke pathways used to tolerate, repair or bypass DNA damage will be fundamental in fully understanding the causes of mitochondrial genome mutations and mitochondrial dysfunction.

Association between telomere length and mitochondrial copy number and cancer risk in humans: A meta-analysis on more than 300,000 individuals

In the last decades the association of leukocyte telomere length (LTL) and mitochondrial copy number (mtDNAcn) with cancer risk has been the focus of many reports, however the relation is not yet completely understood. A meta-analysis of 112 studies including 64,184 cancer cases and 278,641 controls that analysed LTL and mtDNAcn in relation to cancer risk has been conducted to further our understanding of the topic. Stratified analyses for tumor type were also performed. Overall, no association was observed for all cancer combined neither for LTL nor mtDNAcn. Significant associations were detected for these biomarkers and specific cancer type; however, a large degree of heterogeneity was present, even within the same tumor type. Alternatives approaches based on polymorphic variants, such as polygenic risk scores and mendelian randomization, could be adopted to unravel the causal correlation of telomere length and mitochondrial copy number with cancer risk.

Identification of a Nuclear Mitochondrial-Related Multi-Genes Signature to Predict the Prognosis of Bladder Cancer

INTRODUCTION

Bladder cancer (BC) is one of the most prevalent urinary cancers, and its management is still a problem causing recurrence and progression, elevating mortality.

MATERIALS AND METHODS

We aimed at the nuclear mitochondria-related genes (MTRGs), collected from the MITOMAP: A Human Mitochondrial Genome Database. Meanwhile, the expression profiles and clinical information of BC were downloaded from the Cancer Genome Atlas (TCGA) as a training group. The univariate, multivariate, and the least absolute shrinkage and selection operator (LASSO) Cox regression analyses were used to construct a nuclear mitochondrial-related multi-genes signature and the prognostic nomogram.

RESULTS

A total of 17 nuclear MTRGs were identified to be correlated with the overall survival (OS) of BC patients, and a nuclear MTRGs signature based on 16 genes expression was further determined by the LASSO Cox regression analysis. Based on a nuclear MTRGs scoring system, BC patients from the TCGA cohort were divided into high- and low- nuclear MTRGs score groups. Patients with a high nuclear MTRGs score exhibited a significantly poorer outcome (median OS: 92.90 vs 20.20 months, p<0.0001). The nuclear MTRGs signature was further verified in three independent datasets, namely, GSE13507, GSE31684, and GSE32548, from the Gene Expression Omnibus (GEO). The BC patients with a high nuclear MTRGs score had significantly worse survival (median OS in GSE13507: 31.52 vs 98.00 months, p<0.05; GSE31684: 32.85 months vs unreached, p<0.05; GSE32548: unreached vs unreached, p<0.05). Furthermore, muscle-invasive bladder cancer (MIBC) patients had a significantly higher nuclear MTRGs score (p<0.05) than non-muscle-invasive bladder cancer (NMIBC) patients. The integrated signature outperformed each involved MTRG. In addition, a nuclear MTRGs-based nomogram was constructed as a novel prediction prognosis model, whose AUC values for OS at 1, 3, 5 years were 0.76, 0.75, and 0.75, respectively, showing the prognostic nomogram had good and stable predicting ability. Enrichment analyses of the hallmark gene set and KEGG pathway revealed that the E2F targets, G2M checkpoint pathways, and cell cycle had influences on the survival of BC patients. Furthermore, the analysis of tumor microenvironment indicated more CD8+ T cells and higher immune score in patients with high nuclear MTRGs score, which might confer sensitivity to immune checkpoint inhibitors.

CONCLUSIONS

Not only could the signature and prognostic nomogram predict the prognosis of BC, but it also had potential therapeutic guidance.

Circulating Cell-Free mtDNA Content as a Non-invasive Prognostic Biomarker in HCC Patients Receiving TACE and Traditional Chinese Medicine

Hepatocellular carcinoma (HCC) accounts for 70-85% of liver cancer, and about 85% of HCC are hepatitis B virus-related (HBV-HCC) in China. Transarterial chemoembolization (TACE) combined with traditional Chinese medicine (TCM) has been reported as an effective treatment. Potential biomarkers to stratify patients who may benefit from this treatment are needed. In this study, we aimed to evaluate whether circulating cell-free mitochondrial DNA (ccf-mtDNA) content was associated with the outcome of HCC patients, especially of those who received the combination treatment of TACE and TCM. Univariate and multivariate Cox analyses were conducted to evaluate the association between ccf-mtDNA content and the overall survival of HBV-HCC patients. Kaplan-Meier analysis was used to compare the survival differences between patients with low and high ccf-mtDNA content. In a hospital-based cohort with 141 HBV-HCC patients, there was no statistically significant association between the ccf-mtDNA content and the overall survival of HBV-HCC patients in the univariate analysis, but a borderline significant association was found in the multivariate analyses. In a subcohort of 50 HBV-HCC patients who received TACE and TCM treatment, high ccfDNA content conferred an increased death risk with a hazard ratio of 4.01 (95% confidence interval: 1.25-12.84, p = 0.019) in the multivariate analysis. Kaplan-Meier survival analysis also showed that patients with high ccf-mtDNA content had unfavorable survival (log rank p = 0.097). Our findings suggest that ccf-mtDNA content is a potential non-invasive prognostic biomarker in HCC patients receiving TACE and TCM treatment.

Mitochondrial Metabolism Targeted Nanoplatform for Efficient Triple-Negative Breast Cancer Combination Therapy

Tumor reprogram pathway of mitochondrial metabolism is an emerging approach for malignant tumor treatment, such as triple-negative breast cancer. In this study, a tumor/mitochondria cascaded targeting, adenosine-triphosphate (ATP) responsive nanocarrier of zeolitic imidazolate framework-90 (ZIF-90) for breast cancer combination therapy is reported. Atovaquone (AVO) and hemin are loaded into ZIF-90, then a peptide iRGD with tumor-targeting ability is modified on the ZIF-90 nanoplatform. Hemin can specifically degrade BTB and CNC homology1 (BACH1), resulting in the changes of mitochondrial metabolism, and AVO acts as the inhibitor of the electron transport chain (ETC). The degradation of BACH1 using hemin can effectively improve the anti-tumor efficiency of mitochondrial metabolism inhibitor AVO, by increasing dependency on mitochondrial respiration. This nanoplatform displays both tumor-targeting and mitochondria-targeting capacity with high level of ATP responsive drug release behavior. The specific characteristic of mitochondria-targeting ability of this nanoplatform can increase the accumulation of AVO in the mitochondria, and in turn, can effectively improve the inhibition of the ETC. Both in vitro and in vivo results reveal that this composite nanocarrier has excellent tumor inhibition ability with limited side effects. Accordingly, this study provides an attractive strategy in the mitochondrial metabolism for cancer targeted therapy.

EVT-701 is a novel selective and safe mitochondrial complex 1 inhibitor with potent anti-tumor activity in models of solid cancers

Targeting the first protein complex of the mitochondrial electron transport chain (MC1) in cancer has become an attractive therapeutic approach in the recent years, given the metabolic vulnerabilities of cancer cells. The anticancer effect exerted by the pleiotropic drug metformin and the associated reduction in hypoxia-inducible factor 1α (HIF-1α) levels putatively mediated by MC1 inhibition led to the development of HIF-1α inhibitors, such as BAY87-2243, with a more specific MC1 targeting. However, the development of BAY87-2243 was stopped early in phase 1 due to dose-independent emesis and thus there is still no clinical proof of concept for the approach. Given the importance of mitochondrial metabolism during cancer progression, there is still a strong therapeutic need to develop specific and safe MC1 inhibitors. We recently reported the synthesis of compounds with a novel chemotype and potent action on HIF-1α degradation and MC1 inhibition. We describe here the selectivity, safety profile and anti-cancer activity in solid tumors of lead compound EVT-701. In addition, using murine models of lung cancer and of Non-Hodgkin's B cell lymphoma we demonstrated that EVT-701 reduced tumor growth and lymph node invasion when used as a single agent therapy. LKB1 deficiency in lung cancer was identified as a potential indicator of accrued sensitivity to EVT-701, allowing stratification and selection of patients in clinical trials. Altogether these results support further evaluation of EVT-701 alone or in combination in preclinical models and eventually in patients.

Emerging methods for and novel insights gained by absolute quantification of mitochondrial DNA copy number and its clinical applications

The past thirty years have seen a surge in interest in pathophysiological roles of mitochondria, and the accurate quantification of mitochondrial DNA copy number (mCN) in cells and tissue samples is a fundamental aspect of assessing changes in mitochondrial health and biogenesis. Quantification of mCN between studies is surprisingly variable due to a combination of physiological variability and diverse protocols being used to measure this endpoint. The advent of novel methods to quantify nucleic acids like digital polymerase chain reaction (dPCR) and high throughput sequencing offer the ability to measure absolute values of mCN. We conducted an in-depth survey of articles published between 1969 -- 2020 to create an overview of mCN values, to assess consensus values of tissue-specific mCN, and to evaluate consistency between methods of assessing mCN. We identify best practices for methods used to assess mCN, and we address the impact of using specific loci on the mitochondrial genome to determine mCN. Current data suggest that clinical measurement of mCN can provide diagnostic and prognostic value in a range of diseases and health conditions, with emphasis on cancer and cardiovascular disease, and the advent of means to measure absolute mCN should improve future clinical applications of mCN measurements.

Molecular differences in mitochondrial DNA genomes of dogs with malignant mammary tumours

The aim of this study was to determine molecular defects in mitochondrial DNA (mtDNA) with the use of large-scale genome analysis in malignant canine mammary gland tumours and indicate whether these changes were linked with the carcinogenesis process. With the use of the NGS technology, we sequenced 27 samples of mtDNA isolated from blood and tumours obtained from 13 dogs with mammary gland tumours. The total number of mutations and polymorphisms in the analysed mitochondrial genomes was 557. We identified 383 single nucleotide polymorphisms (SNP), 32 indels (or length polymorphisms), 4 mutations, 137 heteroplasmic positions and 1 indel mutation. The highest variability (132 changes) was observed in the variable number of tandem repeats (VNTR) region. The heteroplasmy rate in VNTR varied among individuals and even between two tumours in one organism. Our previous study resulted in determination of a probable CpG island in this region, thus it is not excluded that these changes might alter mtDNA methylation. Only the ATP8 gene was not affected by any polymorphisms or mutations, whereas the COX1 gene had the highest number of polymorphisms from all protein-coding genes. One change m.13594G>A was detected in a region spanning two genes: ND5 and ND6, from which a deleterious effect was observed for the ND5 protein. Molecular changes were frequently observed in the TΨC loop, which is thought to interact with ribosomal RNA.

Mitochondrial DNA and MitomiR Variations in Pancreatic Cancer: Potential Diagnostic and Prognostic Biomarkers

Pancreatic cancer is an aggressive disease with poor prognosis. Only about 15-20% of patients diagnosed with pancreatic cancer can undergo surgical resection, while the remaining 80% are diagnosed with locally advanced or metastatic pancreatic ductal adenocarcinoma (PDAC). In these cases, chemotherapy and radiotherapy only confer marginal survival benefit. Recent progress has been made in understanding the pathobiology of pancreatic cancer, with a particular effort in discovering new diagnostic and prognostic biomarkers, novel therapeutic targets, and biomarkers that can predict response to chemo- and/or radiotherapy. Mitochondria have become a focus in pancreatic cancer research due to their roles as powerhouses of the cell, important subcellular biosynthetic factories, and crucial determinants of cell survival and response to chemotherapy. Changes in the mitochondrial genome (mtDNA) have been implicated in chemoresistance and metastatic progression in some cancer types. There is also growing evidence that changes in microRNAs that regulate the expression of mtDNA-encoded mitochondrial proteins (mitomiRs) or nuclear-encoded mitochondrial proteins (mitochondria-related miRs) could serve as diagnostic and prognostic cancer biomarkers. This review discusses the current knowledge on the clinical significance of changes of mtDNA, mitomiRs, and mitochondria-related miRs in pancreatic cancer and their potential role as predictors of cancer risk, as diagnostic and prognostic biomarkers, and as molecular targets for personalized cancer therapy.

Genomic and evolutionary classification of lung cancer in never smokers

Lung cancer in never smokers (LCINS) is a common cause of cancer mortality but its genomic landscape is poorly characterized. Here high-coverage whole-genome sequencing of 232 LCINS showed 3 subtypes defined by copy number aberrations. The dominant subtype (piano), which is rare in lung cancer in smokers, features somatic UBA1 mutations, germline AR variants and stem cell-like properties, including low mutational burden, high intratumor heterogeneity, long telomeres, frequent KRAS mutations and slow growth, as suggested by the occurrence of cancer drivers' progenitor cells many years before tumor diagnosis. The other subtypes are characterized by specific amplifications and EGFR mutations (mezzo-forte) and whole-genome doubling (forte). No strong tobacco smoking signatures were detected, even in cases with exposure to secondhand tobacco smoke. Genes within the receptor tyrosine kinase-Ras pathway had distinct impacts on survival; five genomic alterations independently doubled mortality. These findings create avenues for personalized treatment in LCINS.

Mitochondrial Dysfunction in Diseases, Longevity, and Treatment Resistance: Tuning Mitochondria Function as a Therapeutic Strategy

Mitochondria are very important intracellular organelles because they have various functions. They produce ATP, are involved in cell signaling and cell death, and are a major source of reactive oxygen species (ROS). Mitochondria have their own DNA (mtDNA) and mutation of mtDNA or change the mtDNA copy numbers leads to disease, cancer chemo/radioresistance and aging including longevity. In this review, we discuss the mtDNA mutation, mitochondrial disease, longevity, and importance of mitochondrial dysfunction in cancer first. In the later part, we particularly focus on the role in cancer resistance and the mitochondrial condition such as mtDNA copy number, mitochondrial membrane potential, ROS levels, and ATP production. We suggest a therapeutic strategy employing mitochondrial transplantation (mtTP) for treatment-resistant cancer.

The Role of Mitochondrial DNA Variation in Drug Response: A Systematic Review

Background: The triad of drug efficacy, toxicity and resistance underpins the risk-benefit balance of all therapeutics. The application of pharmacogenomics has the potential to improve the risk-benefit balance of a given therapeutic via the stratification of patient populations based on DNA variants. A growth in the understanding of the particulars of the mitochondrial genome, alongside the availability of techniques for its interrogation has resulted in a growing body of literature examining the impact of mitochondrial DNA (mtDNA) variation upon drug response. Objective: To critically evaluate and summarize the available literature, across a defined period, in a systematic fashion in order to map out the current landscape of the subject area and identify how the field may continue to advance. Methods: A systematic review of the literature published between January 2009 and December 2020 was conducted using the PubMed database with the following key inclusion criteria: reference to specific mtDNA polymorphisms or haplogroups, a core objective to examine associations between mtDNA variants and drug response, and research performed using human subjects or human in vitro models. Results: Review of the literature identified 24 articles reporting an investigation of the association between mtDNA variant(s) and drug efficacy, toxicity or resistance that met the key inclusion criteria. This included 10 articles examining mtDNA variations associated with antiretroviral therapy response, 4 articles examining mtDNA variants associated with anticancer agent response and 4 articles examining mtDNA variants associated with antimicrobial agent response. The remaining articles covered a wide breadth of medications and were therefore grouped together and referred to as "other." Conclusions: Investigation of the impact of mtDNA variation upon drug response has been sporadic to-date. Collective assessment of the associations identified in the articles was inconclusive due to heterogeneous methods and outcomes, limited racial/ethnic groups, lack of replication and inadequate statistical power. There remains a high degree of idiosyncrasy in drug response and this area has the potential to explain variation in drug response in a clinical setting, therefore further research is likely to be of clinical benefit.

Association of Mitochondrial DNA Copy Number and Telomere Length with Prevalent and Incident Cancer and Cancer Mortality in Women: A Prospective Swedish Population-Based Study

Changes in mitochondrial DNA copy number (mtDNA-CN) and telomere length have, separately, been proposed as risk factors for various cancer types. However, those results are conflicting. Here, mtDNA-CN and relative telomere length were measured in 3225 middle-aged women included in a large population-based prospective cohort. The baseline mtDNA-CN in patients with prevalent breast cancer was significantly higher (12.39 copies/µL) than cancer-free individuals. During an average of 15.2 years of follow-up, 520 patients were diagnosed with cancer. Lower mtDNA-CN was associated with decreased risk of genital organ cancer (hazard ratio (HR), 0.84), and shorter telomere length was associated with increased risk of urinary system cancer (HR, 1.79). Furthermore, mtDNA-CN was inversely associated with all-cause (HR, 1.20) and cancer-specific mortality (HR, 1.21) when considering all cancer types. Surprisingly, shorter telomere length was associated with decreased risk of cancer-specific mortality when considering all cancer types (HR, 0.85). Finally, lower mtDNA-CN and shorter telomere length were associated with increased risk of both all-cause and cancer-specific mortality in genital organ cancer patients. In this study population, we found that mtDNA-CN and telomere length were significantly associated with prevalent and incident cancer and cancer mortality. However, these associations were cancer type specific and need further investigation.

Mitochondrial Heteroplasmy Shifting as a Potential Biomarker of Cancer Progression

Cancer is a serious health problem with a high mortality rate worldwide. Given the relevance of mitochondria in numerous physiological and pathological mechanisms, such as adenosine triphosphate (ATP) synthesis, apoptosis, metabolism, cancer progression and drug resistance, mitochondrial genome (mtDNA) analysis has become of great interest in the study of human diseases, including cancer. To date, a high number of variants and mutations have been identified in different types of tumors, which coexist with normal alleles, a phenomenon named heteroplasmy. This mechanism is considered an intermediate state between the fixation or elimination of the acquired mutations. It is suggested that mutations, which confer adaptive advantages to tumor growth and invasion, are enriched in malignant cells. Notably, many recent studies have reported a heteroplasmy-shifting phenomenon as a potential shaper in tumor progression and treatment response, and we suggest that each cancer type also has a unique mitochondrial heteroplasmy-shifting profile. So far, a plethora of data evidencing correlations among heteroplasmy and cancer-related phenotypes are available, but still, not authentic demonstrations, and whether the heteroplasmy or the variation in mtDNA copy number (mtCNV) in cancer are cause or consequence remained unknown. Further studies are needed to support these findings and decipher their clinical implications and impact in the field of drug discovery aimed at treating human cancer.