Comprehensive molecular characterization of mitochondrial genomes in human cancers

HMGA2 drives the IGFBP1/AKT pathway to counteract the increase in P27KIP1 protein levels in mtDNA/RNA-less cancer cells

Recent comprehensive analyses of mtDNA and orthogonal RNA-sequencing data revealed that in numerous human cancers, mtDNA copy numbers and mtRNA amounts are significantly reduced, followed by low respiratory gene expression. Under such conditions (called mt-Low), cells encounter severe cell proliferation defects; therefore, they must acquire countermeasures against this fatal disadvantage during malignant transformation. This study elucidated a countermeasure against the mt-Low condition-induced antiproliferative effects in hepatocellular carcinoma (HCC) cells. The mechanism relied on the architectural transcriptional regulator HMGA2, which was preferably expressed in HCC cells of the mt-Low type in vitro and in vivo. Detailed in vitro analyses suggest that HMGA2 regulates insulin-like growth factor binding protein 1 (IGFBP1) expression, leading to AKT activation, which then phosphorylates the cyclin-dependent kinase inhibitor (CKI), P27KIP1, and facilitates its ubiquitin-mediated degradation. Accordingly, intervention in the HMGA2 function by RNAi resulted in an increase in P27KIP1 levels and an induction of senescence-like cell proliferation inhibition in mt-Low-type HCC cells. Conclusively, the HMGA2/IGFBP1/AKT axis has emerged as a countermeasure against P27KIP1 CKI upregulation under mt-Low conditions, thereby circumventing cell proliferation inhibition and supporting the tumorigenic state. Notably, similar to in vitro cell lines, HMGA2 was likely to regulate IGFBP1 expression in HCC in vivo, thereby contributing to poor patient prognosis. Considering the significant number of cases under mt-Low or the threat of CKI upregulation cancer-wide, the axis is noteworthy as a vulnerability of cancer cells or target for tumor-agnostic therapy inducing irreversible cell proliferation inhibition via CKI upregulation in a large population with cancer.

MicroRNA and Metabolic Profiling of a Primary Ovarian Neuroendocrine Carcinoma Pulmonary-Type Reveals a High Degree of Similarity with Small Cell Lung Cancer

Small cell neuroendocrine carcinoma is most frequently found in the lung (SCLC), but it has been also reported, albeit with a very low incidence, in the ovary. Here, we analyze a case of primary small cell carcinoma of the ovary of pulmonary type (SCCOPT), a rare and aggressive tumor with poor prognosis, whose biology and molecular features have not yet been thoroughly investigated. The patient affected by SCCOPT had a residual tumor following chemotherapy which displayed pronounced similarity with neuroendocrine tumors and lung cancer in terms of its microRNA expression profile and mTOR-downstream activation. By analyzing the metabolic markers of the neoplastic lesion, we established a likely glycolytic signature. In conclusion, this in-depth characterization of SCCOPT could be useful for future diagnoses, possibly aided by microRNA profiling, allowing clinicians to adopt the most appropriate therapeutic strategy.

Implications of mitochondrial DNA variants in pediatric B-cell acute lymphoblastic leukemia

Background

Research on the role of variations in the mitochondrial genome in pathogenesis of acute lymphoblastic leukemia (ALL) has been unfolding at a rapid rate. Our laboratory has previously described higher number of copies of the mitochondrial genomes per cell in pediatric ALL patients as compared to the healthy controls. In the current study, we evaluated the pattern of mitochondrial genome variations in 20 de-novo pediatric B-ALL cases and seven controls. Quantitative real-time Polymerase Chain Reaction was used for estimation of mitochondrial genomes’ copy number in bone marrow samples of each ALL patient and peripheral blood samples of controls. The complete mitochondrial genomes of all samples were sequenced using the Illumina platform.

Results

Sequencing data analysis using multiple mitochondrial genome databases revealed 325 variants in all 27 samples, out of which 221 variants were previously known while 104 were unassigned, new variants. The 325 variants consisted of 7 loss-of-function variants, 131 synonymous variants, 75 missense variants, and 112 non-coding variants. New, missense variants (n = 21) were identified in genes encoding the electron transport chain complexes with most of them encoding ND4, ND5 of complex I. Missense and loss-of-function variants were found to be deleterious by many predictor databases of pathogenicity. MuTect2 identified true somatic variants present only in tumors between patient-sibling pairs and showed overlap with missense and loss-of-function variants. Online MtDNA-server showed heteroplasmic and homoplasmic variants in mitochondrial genome.

Conclusions

The data suggest that some of these variations might have a deleterious impact on the expression of mitochondrial encoded genes with a possible functional relevance in leukemia.

Development and Validation of a Novel Mitochondrion and Ferroptosis-Related Long Non-Coding RNA Prognostic Signature in Hepatocellular Carcinoma

Mitochondrion and ferroptosis are related to tumorigenesis and tumor progression of hepatocellular carcinoma (HCC). Therefore, this study focused on exploring the participation of lncRNAs in mitochondrial dysfunction and ferroptosis using public datasets from The Cancer Genome Atlas (TCGA) database. We identified the mitochondrion- and ferroptosis-related lncRNAs by Pearson's analysis and lasso-Cox regression. Moreover, real-time quantitative reverse transcription PCR (RT-qPCR) was utilized to further confirm the abnormal expression of these lncRNAs. Based on eight lncRNAs, the MF-related lncRNA prognostic signature (LPS) with outstanding stratification ability and prognostic prediction capability was constructed. In addition, functional enrichment analysis and immune cell infiltration analysis were performed to explore the possible functions of lncRNAs and their impact on the tumor microenvironment. The pathways related to G2M checkpoint and MYC were activated, and the infiltration ratio of regulatory T cells and M0 and M2 macrophages was higher in the high-risk group. In conclusion, these lncRNAs may affect mitochondria functions, ferroptosis, and immune cell infiltration in HCC through specific pathways, which may provide valuable insight into the progression and therapies of HCC.

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.

Insights regarding mitochondrial DNA copy number alterations in human cancer (Review)

Mitochondria are the critical organelles involved in various cellular functions. Mitochondrial biogenesis is activated by multiple cellular mechanisms which require a synchronous regulation between mitochondrial DNA (mtDNA) and nuclear DNA (nDNA). The mitochondrial DNA copy number (mtDNA-CN) is a proxy indicator for mitochondrial activity, and its alteration reflects mitochondrial biogenesis and function. Despite the precise mechanisms that modulate the amount and composition of mtDNA, which have not been fully elucidated, mtDNA-CN is known to influence numerous cellular pathways that are associated with cancer and as well as multiple other diseases. In addition, the utility of current technology in measuring mtDNA-CN contributes to its extensive assessment of diverse traits and tumorigenesis. The present review provides an overview of mtDNA-CN variations across human cancers and an extensive summary of the existing knowledge on the regulation and machinery of mtDNA-CN. The current information on the advanced methods used for mtDNA-CN assessment is also presented.

Mitochondrial Genetic and Epigenetic Regulations in Cancer: Therapeutic Potential

Mitochondria are dynamic organelles managing crucial processes of cellular metabolism and bioenergetics. Enabling rapid cellular adaptation to altered endogenous and exogenous environments, mitochondria play an important role in many pathophysiological states, including cancer. Being under the control of mitochondrial and nuclear DNA (mtDNA and nDNA), mitochondria adjust their activity and biogenesis to cell demands. In cancer, numerous mutations in mtDNA have been detected, which do not inactivate mitochondrial functions but rather alter energy metabolism to support cancer cell growth. Increasing evidence suggests that mtDNA mutations, mtDNA epigenetics and miRNA regulations dynamically modify signalling pathways in an altered microenvironment, resulting in cancer initiation and progression and aberrant therapy response. In this review, we discuss mitochondria as organelles importantly involved in tumorigenesis and anti-cancer therapy response. Tumour treatment unresponsiveness still represents a serious drawback in current drug therapies. Therefore, studying aspects related to genetic and epigenetic control of mitochondria can open a new field for understanding cancer therapy response. The urgency of finding new therapeutic regimens with better treatment outcomes underlines the targeting of mitochondria as a suitable candidate with new therapeutic potential. Understanding the role of mitochondria and their regulation in cancer development, progression and treatment is essential for the development of new safe and effective mitochondria-based therapeutic regimens.

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.

Biomarkers of mitochondrial origin: a futuristic cancer diagnostic

Cancer is a highly fatal disease without effective early-stage diagnosis and proper treatment. Along with the oncoproteins and oncometabolites, several organelles from cancerous cells are also emerging as potential biomarkers. Mitochondria isolated from cancer cells are one such biomarker candidates. Cancerous mitochondria exhibit different profiles compared with normal ones in morphology, genomic, transcriptomic, proteomic and metabolic landscape. Here, the possibilities of exploring such characteristics as potential biomarkers through single-cell omics and Artificial Intelligence (AI) are discussed. Furthermore, the prospects of exploiting the biomarker-based diagnosis and its futuristic utilization through circulatory tumor cell technology are analyzed. A successful alliance of circulatory tumor cell isolation protocols and a single-cell omics platform can emerge as a next-generation diagnosis and personalized treatment procedure.

Long read mitochondrial genome sequencing using Cas9-guided adaptor ligation

The mitochondrial genome (mtDNA) is an important source of disease-causing genetic variability, but existing sequencing methods limit understanding, precluding phased measurement of mutations and clear detection of large sporadic deletions. We adapted a method for amplification-free sequence enrichment using Cas9 cleavage to obtain full length nanopore reads of mtDNA. We then utilized the long reads to phase mutations in a patient with an mtDNA-linked syndrome and demonstrated that this method can map age-induced mtDNA deletions. We believe this method will offer deeper insight into our understanding of mtDNA variation.

Identification of Novel Cytochrome C1 (CYC1) Gene Expression in Oral Squamous Cell Carcinoma- An Evaluative Study

Introduction

Cytochrome C1 (CYC1) is an important subunit of mitochondrial complex III and plays a vital role in oxidative phosphorylation (OXPHOS) and reactive oxygen species generation. Overexpression of the CYC1 gene has been implicated in cancer development and its prognosis previously, but unexplored in head-and-neck squamous cell carcinomas (HNSCC), especially oral squamous cell carcinoma (OSCC).

Materials and Methods

CYC1 m-RNA expression and gene alterations were assessed using the Cancer Genome Atlas dataset in HNSCC and validated in OSCC tissues using real-time polymerase chain reaction (RT-PCR). The protein-protein interaction (PPI) network and functional enrichment pathways were also analysed.

Results

A thorough analysis of the TCGA (The Cancer Genome Atlas) database revealed that CYC1 was overexpressed in the HNSCC cases and the increased expression correlated with several parameters which involve the prediction of advanced diseases such as histopathological grade, tumour-node-metastasis staging, and nodal metastases (P < 0.05). The expression of CYC1 was validated using RT-PCR showing significant upregulation (P < 0.05) in OSCC tissue samples compared to the normal tissue counterparts. PPI network and functional analysis show the prominent role of CYC1 in OXPHOS, especially in electron transport chain III complex regulation.

Discussion

The study revealed that CYC1 is highly expressed in HNSCC, and is validated in the OSCC patient tissue samples compared to the normal counterparts and associated with advanced disease stages and grade of the tumour. CYC1 could be a novel promising therapeutic and prognostic marker in HNSCC, especially in OSCC.

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.

Metabolic reprogramming: a bridge between aging and tumorigenesis

Aging is the most robust risk factor for cancer development, with more than 60% of cancers occurring in those aged 60 and above. However, how aging and tumorigenesis are intertwined is poorly understood and a matter of significant debate. Metabolic changes are hallmarks of both aging and tumorigenesis. The deleterious consequences of aging include dysfunctional cellular processes, the build‐up of metabolic byproducts and waste molecules in circulation and within tissues, and stiffer connective tissues that impede blood flow and oxygenation. Collectively, these age‐driven changes lead to metabolic reprogramming in different cell types of a given tissue that significantly affects their cellular functions. Here, we put forward the idea that metabolic changes that happen during aging help create a favorable environment for tumorigenesis. We review parallels in metabolic changes that happen during aging and how these changes function both as adaptive mechanisms that enable the development of malignant phenotypes in a cell‐autonomous manner and as mechanisms that suppress immune surveillance, collectively creating the perfect environment for cancers to thrive. Hence, antiaging therapeutic strategies that target the metabolic reprogramming that occurs as we age might provide new opportunities to prevent cancer initiation and/or improve responses to standard‐of‐care anticancer therapies.

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].

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.

Computational analysis of cancer genome sequencing data

Distilling biologically meaningful information from cancer genome sequencing data requires comprehensive identification of somatic alterations using rigorous computational methods. As the amount and complexity of sequencing data have increased, so has the number of tools for analysing them. Here, we describe the main steps involved in the bioinformatic analysis of cancer genomes, review key algorithmic developments and highlight popular tools and emerging technologies. These tools include those that identify point mutations, copy number alterations, structural variations and mutational signatures in cancer genomes. We also discuss issues in experimental design, the strengths and limitations of sequencing modalities and methodological challenges for the future.

Mitochondrial-Related Transcriptome Feature Correlates with Prognosis, Vascular Invasion, Tumor Microenvironment, and Treatment Response in Hepatocellular Carcinoma

Background

Hepatocellular carcinoma (HCC) is the most common subtype of primary liver cancer, which was highly correlated with metabolic dysfunction. Nevertheless, the association between nuclear mitochondrial-related transcriptome and HCC remained unclear.

Materials and Methods

A total of 147 nuclear mitochondrial-related genes (NMRGs) were downloaded from the MITOMAP: A Human Mitochondrial Genome Database. The training dataset was downloaded from The Cancer Genome Atlas (TCGA), while validation datasets were retrieved from the International Cancer Genome Consortium (ICGC) and Gene Expression Omnibus (GEO). The univariate and multivariate, and least absolute shrinkage and selection operator (LASSO) Cox regression analyses were applied to construct a NMRG signature, and the value of area under receiver operating characteristic curve (AUC) was utilized to assess the signature and nomogram. Then, data from the Genomics of Drug Sensitivity in Cancer (GDSC) were used for the evaluation of chemotherapy response in HCC.

Results

Functional enrichment of differentially expressed genes (DEGs) between HCC and paired normal tissue samples demonstrated that mitochondrial dysfunction was significantly associated with HCC development. Survival analysis showed a total of 35 NMRGs were significantly correlated with overall survival (OS) of HCC, and the LASSO Cox regression analysis further identified a 25-NMRG signature and corresponding prognosis score based on their transcriptional profiling. HCC patients were divided into high- and low-risk groups according to the median prognosis score, and high-risk patients had significantly worse OS (median OS: 27.50 vs. 83.18 months, P < 0.0001). The AUC values for OS at 1, 3, and 5 years were 0.79, 0.77, and 0.77, respectively. The prognostic capacity of NMRG signature was verified in the GSE14520 dataset and ICGC-HCC cohort. Besides, the NMRG signature outperformed each NMRG and clinical features in prognosis prediction and could also differentiate whether patients presented with vascular invasions (VIs) or not. Subsequently, a prognostic nomogram (C-index: 0.753, 95% CI: 0.703~0.804) by the integration of age, tumor metastasis, and NMRG prognosis score was constructed with the AUC values for OS at 1, 3, and 5 years were 0.82, 0.81, and 0.82, respectively. Notably, significant enrichment of regulatory and follicular helper T cells in high-risk group indicated the potential treatment of immune checkpoint inhibitors for these patients. Interestingly, the NMRG signature could also identify the potential responders of sorafenib or transcatheter arterial chemoembolization (TACE) treatment. Additionally, HCC patients in high-risk group appeared to be more sensitive to cisplatin, vorinostat, and methotrexate, reversely, patients in low-risk group had significantly higher sensitivity to paclitaxel and bleomycin instead.

Conclusions

In summary, the development of NMRG signature provided a more comprehensive understanding of mitochondrial dysfunction in HCC, helped predict prognosis and tumor microenvironment, and provided potential targeted therapies for HCC patients with different NMRG prognosis scores.

Comprehensive profiling of 1015 patients' exomes reveals genomic-clinical associations in colorectal cancer

The genetic basis of colorectal cancer (CRC) and its clinical associations remain poorly understood due to limited samples or targeted genes in current studies. Here, we perform ultradeep whole-exome sequencing on 1015 patients with CRC as part of the ChangKang Project. We identify 46 high-confident significantly mutated genes, 8 of which mutate in 14.9% of patients: LYST, DAPK1, CR2, KIF16B, NPIPB15, SYTL2, ZNF91, and KIAA0586. With an unsupervised clustering algorithm, we propose a subtyping strategy that classisfies CRC patients into four genomic subtypes with distinct clinical characteristics, including hypermutated, chromosome instability with high risk, chromosome instability with low risk, and genome stability. Analysis of immunogenicity uncover the association of immunogenicity reduction with genomic subtypes and poor prognosis in CRC. Moreover, we find that mitochondrial DNA copy number is an independent factor for predicting the survival outcome of CRCs. Overall, our results provide CRC-related molecular features for clinical practice and a valuable resource for translational research.

Identification of Somatic Mitochondrial DNA Mutations, Heteroplasmy, and Increased Levels of Catenanes in Tumor Specimens Obtained from Three Endometrial Cancer Patients

Endometrial carcinoma (EC) is the most common type of gynecologic malignant epithelial tumor, with the death rate from this disease doubling over the past 20 years. Mitochondria provide cancer cells with necessary anabolic building blocks such as amino acids, lipids, and nucleotides, and EC samples have been shown to increase mitochondrial biogenesis. In cancer, mitochondrial DNA (mtDNA) heteroplasmy studies suggest that heteroplasmic variants encode predicted pathogenic proteins. We investigated the mtDNA genotypes within peri-normal and tumor specimens obtained from three individuals diagnosed with EC. DNA extracts from peri-normal and tumor tissues were used for mtDNA-specific next-generation sequencing and analyses of mtDNA content and topoisomers. The three tumors harbor heteroplasmic somatic mutations, and at least one mutation in each carcinoma is predicted to deleteriously alter a mtDNA-encoded protein. Somatic heteroplasmy linked to two mtDNA tRNA genes was found in separate tumors, and two heteroplasmic non-coding variants were identified in a single EC tumor. While two tumors had altered mtDNA content, all three displayed increased mtDNA catenanes. Our findings support that EC cells require wild-type mtDNA, but heteroplasmic mutations may alter mitochondrial metabolism to help promote cancer cell growth and proliferation.

The Entanglement between Mitochondrial DNA and Tumor Metastasis

Simple Summary

Mitochondrial dysfunction is one of the main features of cancer cells. As genetic material in mitochondria, mitochondrial DNA (mtDNA) variations and dysregulation of mitochondria-encoded genes have been shown to correlate with survival outcomes in cancer patients. Cancer metastasis is often a major cause of treatment failure, which is a multi-step cascade process. With the development of gene sequencing and in vivo modeling technology, the role of mtDNA in cancer metastasis has been continuously explored. Our review systematically provides a summary of the multiple roles of mtDNA in cancer metastasis and presents the broad prospects for mtDNA in cancer prediction and therapy.

Abstract

Mitochondrial DNA, the genetic material in mitochondria, encodes essential oxidative phosphorylation proteins and plays an important role in mitochondrial respiration and energy transfer. With the development of genome sequencing and the emergence of novel in vivo modeling techniques, the role of mtDNA in cancer biology is gaining more attention. Abnormalities of mtDNA result in not only mitochondrial dysfunction of the the cancer cells and malignant behaviors, but regulation of the tumor microenvironment, which becomes more aggressive. Here, we review the recent progress in the regulation of cancer metastasis using mtDNA and the underlying mechanisms, which may identify opportunities for finding novel cancer prediction and therapeutic targets.

CRISPR-free base editors with enhanced activity and expanded targeting scope in mitochondrial and nuclear DNA

The all-protein cytosine base editor DdCBE uses TALE proteins and a double-stranded DNA-specific cytidine deaminase (DddA) to mediate targeted C•G-to-T•A editing. To improve editing efficiency and overcome the strict TC sequence-context constraint of DddA, we used phage-assisted non-continuous and continuous evolution to evolve DddA variants with improved activity and expanded targeting scope. Compared to canonical DdCBEs, base editors with evolved DddA6 improved mitochondrial DNA (mtDNA) editing efficiencies at TC by 3.3-fold on average. DdCBEs containing evolved DddA11 offered a broadened HC (H = A, C or T) sequence compatibility for both mitochondrial and nuclear base editing, increasing average editing efficiencies at AC and CC targets from less than 10% for canonical DdCBE to 15-30% and up to 50% in cell populations sorted to express both halves of DdCBE. We used these evolved DdCBEs to efficiently install disease-associated mtDNA mutations in human cells at non-TC target sites. DddA6 and DddA11 substantially increase the effectiveness and applicability of all-protein base editing.

Bladder cancer cells shift rapidly and spontaneously to cisplatin-resistant oxidative phosphorylation that is trackable in real time

Genetic mutations have long been recognized as drivers of cancer drug resistance, but recent work has defined additional non-genetic mechanisms of plasticity, wherein cancer cells assume a drug resistant phenotype marked by altered epigenetic and transcriptional states. Currently, little is known about the real-time, dynamic nature of this phenotypic shift. Using a bladder cancer model of nongenetic plasticity, we discovered that rapid transition to drug resistance entails upregulation of mitochondrial gene expression and a corresponding metabolic shift towards the tricarboxylic acid cycle and oxidative phosphorylation. Based on this distinction, we were able to track cancer cell metabolic profiles in real time using fluorescence lifetime microscopy (FLIM). We observed single cells transitioning spontaneously to an oxidative phosphorylation state over hours to days, a trend that intensified with exposure to cisplatin chemotherapy. Conversely, pharmacological inhibition of oxidative phosphorylation significantly reversed the FLIM metabolic signature and reduced cisplatin resistance. These rapid, spontaneous metabolic shifts offer a new means of tracking nongenetic cancer plasticity and forestalling the emergence of drug resistance.

Heteroplasmic mitochondrial DNA variants in cardiovascular diseases

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

A computational workflow for predicting cancer neo-antigens

Neo-antigens presented on cell surface play a pivotal role in the success of immunotherapies. Peptides derived from mutant proteins are thought to be the primary source of neo-antigens presented on the surface of cancer cells. Mutation data from cancer genome sequencing is often used to predict cancer neo-antigens. However, this strategy is associated with significant false positives as many coding mutations may not be expressed at the protein level. Hence, we describe a computational workflow to integrate genomic and proteomic data to predictpotential neo-antigens.

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.

Metabolic tricks of cancer cells

One of the characteristics of cancer cells important for tumorigenesis is their metabolic plasticity. Indeed, in various stress conditions, cancer cells can reshape their metabolic pathways to support the increased energy request due to continuous growth and rapid proliferation. Moreover, selective pressures in the tumor microenvironment, such as hypoxia, acidosis, and competition for resources, force cancer cells to adapt by complete reorganization of their metabolism. In this review, we highlight the characteristics of cancer metabolism and discuss its clinical significance, since overcoming metabolic plasticity of cancer cells is a key objective of modern cancer therapeutics and a better understanding of metabolic reprogramming may lead to the identification of possible targets for cancer therapy.

Mitochondrial Proteins as Source of Cancer Neoantigens

In the past decade, anti-tumour immune responses have been successfully exploited to improve the outcome of patients with different cancers. Significant progress has been made in taking advantage of different types of T cell functions for therapeutic purposes. Despite these achievements, only a subset of patients respond favorably to immunotherapy. Therefore, there is a need of novel approaches to improve the effector functions of immune cells and to recognize the major targets of anti-tumour immunity. A major hallmark of cancer is metabolic rewiring associated with switch of mitochondrial functions. These changes are a consequence of high energy demand and increased macromolecular synthesis in cancer cells. Such adaptations in tumour cells might generate novel targets of tumour therapy, including the generation of neoantigens. Here, we review the most recent advances in research on the immune response to mitochondrial proteins in different cellular conditions.

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.

Mitochondria in tumour progression: a network of mtDNA variants in different types of cancer

Background

Mitochondrial participation in tumorigenesis and metastasis has been studied for many years, but several aspects of this mechanism remain unclear, such as the association of mitochondrial DNA (mtDNA) with different cancers. Here, based on two independent datasets, we modelled an mtDNA mutation-cancer network by systematic integrative analysis including 37 cancer types to identify the mitochondrial variants found in common among them.

Results

Our network showed mtDNA associations between gastric cancer and other cancer types, particularly kidney, liver, and prostate cancers, which is suggestive of a potential role of such variants in the metastatic processes among these cancer types. A graph-based interactive web tool was made available at www2.lghm.ufpa.br/mtdna. We also highlighted that most shared variants were in the MT-ND4, MT-ND5 and D-loop, and that some of these variants were nonsynonymous, indicating a special importance of these variants and regions regarding cancer progression, involving genomic and epigenomic alterations.

Conclusions

This study reinforces the importance of studying mtDNA in cancer and offers new perspectives on the potential involvement of different mitochondrial variants in cancer development and metastasis.

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.

Circulating Mitochondrial DNA and Inter-Organelle Contact Sites in Aging and Associated Conditions

Mitochondria are primarily involved in cell bioenergetics, regulation of redox homeostasis, and cell death/survival signaling. An immunostimulatory property of mitochondria has also been recognized which is deployed through the extracellular release of entire or portioned organelle and/or mitochondrial DNA (mtDNA) unloading. Dynamic homo- and heterotypic interactions involving mitochondria have been described. Each type of connection has functional implications that eventually optimize mitochondrial activity according to the bioenergetic demands of a specific cell/tissue. Inter-organelle communications may also serve as molecular platforms for the extracellular release of mitochondrial components and subsequent ignition of systemic inflammation. Age-related chronic inflammation (inflamm-aging) has been associated with mitochondrial dysfunction and increased extracellular release of mitochondrial components—in particular, cell-free mtDNA. The close relationship between mitochondrial dysfunction and cellular senescence further supports the central role of mitochondria in the aging process and its related conditions. Here, we provide an overview of (1) the mitochondrial genetic system and the potential routes for generating and releasing mtDNA intermediates; (2) the pro-inflammatory pathways elicited by circulating mtDNA; (3) the participation of inter-organelle contacts to mtDNA homeostasis; and (4) the link of these processes with senescence and age-associated conditions.

mitoDataclean: A machine learning approach for the accurate identification of cross-contamination-derived tumor mitochondrial DNA mutations

Next-generation sequencing (NGS) of mitochondrial DNA (mtDNA) has widespread applications in aging and cancer studies. However, cross-contamination of mtDNA constitutes a major concern. Previous methods for the detection of mtDNA contamination mainly focus on haplogroup-level phylogeny, but neglect haplotype-level differences, leading to limited sensitivity and accuracy. In this study, we present mitoDataclean, a random-forest-based machine learning package for accurate identification of cross-contamination, evaluation of contamination levels and detection of contamination-derived variants in mtDNA NGS data. Comprehensive optimization of mitoDataclean revealed that training simulation with mixtures of small haplogroup distance and low polymorphic difference was critical for optimal modeling. Compared with existing methods, mitoDataclean exhibited significantly improved sensitivity and accuracy for the detection of sample contamination in simulated data. In addition, mitoDataclean achieved area under the curve values of 0.91 and 0.97 for discerning genuine and contamination-derived mtDNA variants in a simulated Western dataset and private sequencing contamination data, respectively, suggesting that this tool may be applicable for different populations and samples with different sources of contamination. Finally, mitoDataclean was further evaluated in several private and public datasets and showed a robust ability for contamination detection. Altogether, our study demonstrates that mitoDataclean may be used for accurate detection of contaminated samples and contamination-derived variants in mtDNA NGS data. This article is protected by copyright. All rights reserved.

Dynamic Cancer Cell Heterogeneity: Diagnostic and Therapeutic Implications

Though heterogeneity of cancers is recognized and has been much discussed in recent years, the concept often remains overlooked in different routine examinations. Indeed, in clinical or biological articles, reviews, and textbooks, cancers and cancer cells are generally presented as evolving distinct entities rather than as an independent heterogeneous cooperative cell population with its self-oriented biology. There are, therefore, conceptual gaps which can mislead the interpretations/diagnostic and therapeutic approaches. In this short review, we wish to summarize and discuss various aspects of this dynamic evolving heterogeneity and its biological, pathological, clinical, diagnostic, and therapeutic implications, using thyroid carcinoma as an illustrative example.

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.

Hypermethylation of mitochondrial DNA facilitates bone metastasis of renal cell carcinoma

Kidney cancers including clear cell carcinoma (RCC) are identified with very vulnerable mitochondria DNA (mtDNA) and frequent epigenetic aberrations. Bone metastasis from RCC is prevalent and destructive. Bone marrow contains a quite hypoxic microenvironment that usually insitigate 50% of hypermethylation events in conferring a selective advantage for tumor growth. We hypothesized that hypermethylation of mtDNA in RCC cells would significantly contribute to bone metastatic tumor progression. Methylation-specific polymerase chain reaction assay (MSP) was adopted to measure the methylation status of D-loop region of mtDNA in 15 pairs of bone metastatic and primary RCC as well as tumor adjescent normal kidney tissues. mtDNA copy number was examined by the real-time quantitative polymerase chain reaction (qPCR). Western blotting analysis was used to measure the accumulation of several DNA methyltransferases (DNMTs) in the mitochondria and nucleus fractions of bone metastatic RCC cells. mRNA expression of mitochondria encoded genes was examined by RT-PCR. Reactive oxygen species (ROS), mitochondrial membrane potential and ATP content were measured using in vitro cells treated with de-methylation drug 5-Azacytidine (5-Aza). Non-invasive bioluminescent imaging was performed to monitor tumor occurrence in skeleton in mice. Our results showed that the D-loop region in bone metastatic tumor cells was markedly hypermethylated than those in primary RCC tumor cells, that is associated with a decreased mtDNA copy number and accumulation of DNMT1 in the mitochondria. The bone-tropism tumor colonization and progression of RCC cells was significantly suppressed by demethylating the D-loop region of mtDNA and reducing the intracellular level of ROS and ATP by 5-Aza treatment. In conclusion, our study provided a direct association between hypermethylation of mtDNA in RCC with bone metastastic tumor growth.

The frequency of the known mitochondrial variants associated with drug-induced toxicity in a Korean population

BACKGROUND

Few studies have annotated the whole mitochondrial DNA (mtDNA) genome associated with drug responses in Asian populations. This study aimed to characterize mtDNA genetic profiles, especially the distribution and frequency of well-known genetic biomarkers associated with diseases and drug-induced toxicity in a Korean population.

METHOD

Whole mitochondrial genome was sequenced for 118 Korean subjects by using a next-generation sequencing approach. The bioinformatic pipeline was constructed for variant calling, haplogroup classification and annotation of mitochondrial mutation.

RESULTS

A total of 681 variants was identified among all subjects. The MT-TRNP gene and displacement loop showed the highest numbers of variants (113 and 74 variants, respectively). The m.16189T > C allele, which is known to reduce the mtDNA copy number in human cells was detected in 25.4% of subjects. The variants (m.2706A > G, m.3010A > G, and m.1095T > C), which are associated with drug-induced toxicity, were observed with the frequency of 99.15%, 30.51%, and 0.08%, respectively. The m.2150T > A, a genotype associated with highly disruptive effects on mitochondrial ribosomes, was identified in five subjects. The D and M groups were the most dominant groups with the frequency of 34.74% and 16.1%, respectively.

CONCLUSIONS

Our finding was consistent with Korean Genome Project and well reflected the unique profile of mitochondrial haplogroup distribution. It was the first study to annotate the whole mitochondrial genome with drug-induced toxicity to predict the ADRs event in clinical implementation for Korean subjects. This approach could be extended for further study for validation of the potential ethnic-specific mitochondrial genetic biomarkers in the Korean population.

Mitochondrial metabolism as a potential therapeutic target in myeloid leukaemia

While the understanding of the genomic aberrations that underpin chronic and acute myeloid leukaemia (CML and AML) has allowed the development of therapies for these diseases, limitations remain. These become apparent when looking at the frequency of treatment resistance leading to disease relapse in leukaemia patients. Key questions regarding the fundamental biology of the leukaemic cells, such as their metabolic dependencies, are still unresolved. Even though a majority of leukaemic cells are killed during initial treatment, persistent leukaemic stem cells (LSCs) and therapy-resistant cells are still not eradicated with current treatments, due to various mechanisms that may contribute to therapy resistance, including cellular metabolic adaptations. In fact, recent studies have shown that LSCs and treatment-resistant cells are dependent on mitochondrial metabolism, hence rendering them sensitive to inhibition of mitochondrial oxidative phosphorylation (OXPHOS). As a result, rewired energy metabolism in leukaemic cells is now considered an attractive therapeutic target and the significance of this process is increasingly being recognised in various haematological malignancies. Therefore, identifying and targeting aberrant metabolism in drug-resistant leukaemic cells is an imperative and a relevant strategy for the development of new therapeutic options in leukaemia. In this review, we present a detailed overview of the most recent studies that present experimental evidence on how leukaemic cells can metabolically rewire, more specifically the importance of OXPHOS in LSCs and treatment-resistant cells, and the current drugs available to target this process. We highlight that uncovering specific energy metabolism dependencies will guide the identification of new and more targeted therapeutic strategies for myeloid leukaemia.

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).

Racial disparities in the genetic landscape of lung cancer

Lung cancer has the highest cancer-related mortality worldwide and in the United States. Although reduced tobacco consumption and advancement in therapies have led to a modest decline in lung cancer death rates over the past two decades; the overall survival rate is still disappointing. Moreover, race-associated disparities are also observed, especially in the clinical outcomes. Socioeconomic factors are considered major contributors in cancer health disparities, however, the differences in the genetic landscape of lung cancer among different racial groups have also been reported. In this review, we shed light on the genetic heterogeneity of lung cancer and race-associated differences in genetic alterations to build a framework for future studies to understand the biological basis of lung cancer disparities.

Identifying Somatic Mitochondrial DNA Mutations

Mitochondria are cellular organelles that play an essential role in eukaryotes, producing the energy needed for a cell to survive. Beyond the ~3.2 Gb of nuclear genomic DNA, each human cell has hundreds of mitochondria which carry one or a few copies of the 16.5 kb circular mitochondria DNA (mtDNA). Despite its small size, the circular genome encodes 37 genes, including 13 proteins that generate respiratory chain complexes together with other proteins of nuclear origin. Similar to nuclear genome, mtDNA in cancer cells frequently harbor somatically acquired alterations. Whole-genome or whole-exome sequencing of the tumor and its matched normal tissues (frequently blood or adjacent non-tumor tissues) enables sensitive and efficient detection of somatic mtDNA mutations. Because each cancer cell commonly carries hundreds to thousands of mtDNA copies, detection of mtDNA mutations is dependent on the heteroplasmic level of each mutation. Here, we describe strategies to accurately identify somatic mtDNA mutations in cancer genome studies.

Efficient exosome subpopulation isolation and proteomic profiling using a Sub-ExoProfile chip towards cancer diagnosis and treatment

Deconstruction of the heterogeneity of surface marker-dependent exosome subpopulations by the Sub-ExoProfile chip.

Genetic landscape of human mitochondrial genome using whole genome sequencing

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

Metabolic therapy and bioenergetic analysis: The missing piece of the puzzle

BACKGROUND

Aberrant metabolism is recognized as a hallmark of cancer, a pillar necessary for cellular proliferation. Regarding bioenergetics (ATP generation), most cancers display a preference not only toward aerobic glycolysis ("Warburg effect") and glutaminolysis (mitochondrial substrate level-phosphorylation) but also toward other metabolites such as lactate, pyruvate, and fat-derived sources. These secondary metabolites can assist in proliferation but cannot fully cover ATP demands.

SCOPE OF REVIEW

The concept of a static metabolic profile is challenged by instances of heterogeneity and flexibility to meet fuel/anaplerotic demands. Although metabolic therapies are a promising tool to improve therapeutic outcomes, either via pharmacological targets or press-pulse interventions, metabolic plasticity is rarely considered. Lack of bioenergetic analysis in vitro and patient-derived models is hindering translational potential. Here, we review the bioenergetics of cancer and propose a simple analysis of major metabolic pathways, encompassing both affordable and advanced techniques. A comprehensive compendium of Seahorse XF bioenergetic measurements is presented for the first time.

MAJOR CONCLUSIONS

Standardization of principal readouts might help researchers to collect a complete metabolic picture of cancer using the most appropriate methods depending on the sample of interest.

Rationale and design of the CRAFT (Continuous ReAssessment with Flexible ExTension in Rare Malignancies) multicenter phase II trial

BACKGROUND

Approvals of cancer therapeutics are primarily disease entity specific. Current molecular diagnostic approaches frequently identify actionable alterations in rare cancers or rare subtypes of common cancers for which the corresponding treatments are not approved and unavailable within clinical trials due to entity-related eligibility criteria. Access may be negotiated with health insurances. However, approval rates vary, and critical information required for a scientific evaluation of treatment-associated risks and benefits is not systematically collected. Thus clinical trials with optimized patient selection and comprehensive molecular characterization are essential for translating experimental treatments into standard care.

PATIENTS AND METHODS

Continuous ReAssessment with Flexible ExTension in Rare Malignancies (CRAFT) is an open-label phase II trial for adults with pretreated, locally advanced, or metastatic solid tumors. Based on the evaluation by a molecular tumor board, patients are assigned to combinations of six molecularly targeted agents and a programmed death-ligand 1 (PD-L1) antagonist within seven study arms focusing on (i) BRAF V600 mutations; (ii) ERBB2 amplification and/or overexpression, activating ERBB2 mutations; (iii) ALK rearrangements, activating ALK mutations; (iv and v) activating PIK3CA and AKT mutations, other aberrations predicting increased PI3K-AKT pathway activity; (vi) aberrations predicting increased RAF-MEK-ERK pathway activity; (vii) high tumor mutational burden and other alterations predicting sensitivity to PD-L1 inhibition. The primary endpoint is the disease control rate (DCR) at week 16; secondary and exploratory endpoints include the progression-free survival ratio, overall survival, and patient-reported outcomes. Using Simon's optimal two-stage design, 14 patients are accrued for each study arm. If three or fewer patients achieve disease control, the study arm is stopped. Otherwise, 11 additional patients are accrued. If the DCR exceeds 7 of 25 patients, the null hypothesis is rejected for the respective study arm.

CONCLUSIONS

CRAFT was activated in October 2021 and will recruit at 10 centers in Germany.

TRIAL REGISTRATION NUMBERS

EudraCT: 2019-003192-18; ClinicalTrials.gov: NCT04551521.

Longitudinal Single-Cell Dynamics of Chromatin Accessibility and Mitochondrial Mutations in Chronic Lymphocytic Leukemia Mirror Disease History

While cancers evolve during disease progression and in response to therapy, temporal dynamics remain difficult to study in humans due to the lack of consistent barcodes marking individual clones in vivo. We employ mitochondrial single-cell assay for transposase-accessible chromatin with sequencing to profile 163,279 cells from 9 patients with chronic lymphocytic leukemia (CLL) collected across disease course and utilize mitochondrial DNA (mtDNA) mutations as natural genetic markers of cancer clones. We observe stable propagation of mtDNA mutations over years in the absence of strong selective pressure, indicating clonal persistence, but dramatic changes following tight bottlenecks, including disease transformation and relapse posttherapy, paralleled by acquisition of copy-number variants and changes in chromatin accessibility and gene expression. Furthermore, we link CLL subclones to distinct chromatin states, providing insight into nongenetic sources of relapse. mtDNA mutations thus mirror disease history and provide naturally occurring genetic barcodes to enable patient-specific study of cancer subclonal dynamics.

SIGNIFICANCE

Single-cell multi-omic profiling of CLL reveals the utility of somatic mtDNA mutations as in vivo barcodes, which mark subclones that can evolve over time along with changes in accessible chromatin and gene expression profiles to capture dynamics of disease evolution. See related commentary by Hilton and Scott, p. 2965. This article is highlighted in the In This Issue feature, p. 2945.

Impaired TFAM expression promotes mitochondrial damage to drive fibroblast activation and fibrosis in systemic sclerosis

OBJECTIVES

The transcription factor TFAM is controlling the transcription of core proteins required for mitochondrial homeostasis. The aim of the current study was to investigate changes in TFAM expression in systemic sclerosis (SSc), to analyze mitochondrial function and to evaluate the consequences for fibroblast activation.

METHODS

The expression of TFAM was analyzed by immunofluorescence and Western blot. The effects of TFAM knockout were investigated in cultured fibroblasts and in bleomycin-induced skin and lung fibrosis and in TβRI^act -induced skin fibrosis.

RESULTS

The expression of TFAM was downregulated in fibroblasts in SSc skin and in cultured SSc fibroblasts. The downregulation of TFAM was associated with decreased mitochondrial number and accumulation of damaged mitochondria with release of mtDNA, accumulation of deletions in mtDNA, metabolic alterations with impaired OXPHOS and release of the mitokine GDF15. Chronic, but not acute, exposure of normal fibroblasts to TGFβ mimicked the finding in SSc fibroblasts with downregulation of TFAM and accumulation of mitochondrial damage. Downregulation of TFAM promotes fibroblast activation with upregulation of fibrosis-relevant GO-terms in RNASeq, partially in a ROS-dependent manner. Mice with fibroblast-specific knockout of TFAM are prone to fibrotic tissue remodeling with fibrotic responses even to NaCl instillation and enhanced sensitivity to bleomycin injection and TβRIact-overexpression. TFAM knockout fosters SMAD3 signaling to promote fibroblast activation.

CONCLUSIONS

Alterations in the key mitochondrial transcription factor TFAM in response to prolonged activation of TGFβ and associated mitochondrial damage induce transcriptional programs that promote fibroblast-to-myofibroblast transition and drive tissue fibrosis.

MitoScape: A big-data, machine-learning platform for obtaining mitochondrial DNA from next-generation sequencing data

The growing number of next-generation sequencing (NGS) data presents a unique opportunity to study the combined impact of mitochondrial and nuclear-encoded genetic variation in complex disease. Mitochondrial DNA variants and in particular, heteroplasmic variants, are critical for determining human disease severity. While there are approaches for obtaining mitochondrial DNA variants from NGS data, these software do not account for the unique characteristics of mitochondrial genetics and can be inaccurate even for homoplasmic variants. We introduce MitoScape, a novel, big-data, software for extracting mitochondrial DNA sequences from NGS. MitoScape adopts a novel departure from other algorithms by using machine learning to model the unique characteristics of mitochondrial genetics. We also employ a novel approach of using rho-zero (mitochondrial DNA-depleted) data to model nuclear-encoded mitochondrial sequences. We showed that MitoScape produces accurate heteroplasmy estimates using gold-standard mitochondrial DNA data. We provide a comprehensive comparison of the most common tools for obtaining mtDNA variants from NGS and showed that MitoScape had superior performance to compared tools in every statistically category we compared, including false positives and false negatives. By applying MitoScape to common disease examples, we illustrate how MitoScape facilitates important heteroplasmy-disease association discoveries by expanding upon a reported association between hypertrophic cardiomyopathy and mitochondrial haplogroup T in men (adjusted p-value = 0.003). The improved accuracy of mitochondrial DNA variants produced by MitoScape will be instrumental in diagnosing disease in the context of personalized medicine and clinical diagnostics.

Recent studies have highlighted the importance of mitochondrial DNA variation in both primary mitochondrial disease and complex, human pathology including COVID-19, and space-flight stress. The vast amount of existing, next-generation sequencing (NGS) data can be leveraged to interrogate both nuclear and mitochondrial DNA (mtDNA) sequence simultaneously, allowing for analysis of the interplay between mitochondrial and nuclear encoded genes in mitochondrial function. Identifying mtDNA sequence accurately is complicated by the presence of nuclear encoded mitochondrial sequences (NUMTs), which are homologous to mtDNA. Current software for analyzing mtDNA from NGS do not accurately model the unique characteristics of mitochondrial genetics. We introduce MitoScape, a novel, big-data, software which models mitochondrial genetics through machine learning to accurately identify mtDNA sequence from NGS data. MitoScape takes advantage of rho-zero cell data to model the characteristics of NUMTs. We show that MitoScape produces more accurate heteroplasmy estimates compared to published software. We provide an example of applying MitoScape in replicating an association between hypertrophic cardiomyopathy and mitochondrial haplogroup T in men. MitoScape is an important contribution to mitochondrial genomics allowing for accurate mtDNA variants, and the ability to tailor mtDNA analysis in different population and disease contexts, which is not available in other software.

A replication-linked mutational gradient drives somatic mutation accumulation and influences germline polymorphisms and genome composition in mitochondrial DNA

Mutations in mitochondrial DNA (mtDNA) cause maternally inherited diseases, while somatic mutations are linked to common diseases of aging. Although mtDNA mutations impact health, the processes that give rise to them are under considerable debate. To investigate the mechanism by which de novo mutations arise, we analyzed the distribution of naturally occurring somatic mutations across the mouse and human mtDNA obtained by Duplex Sequencing. We observe distinct mutational gradients in G→A and T→C transitions delimited by the light-strand origin and the mitochondrial Control Region (mCR). The gradient increases unequally across the mtDNA with age and is lost in the absence of DNA polymerase γ proofreading activity. In addition, high-resolution analysis of the mCR shows that important regulatory elements exhibit considerable variability in mutation frequency, consistent with them being mutational ‘hot-spots’ or ‘cold-spots’. Collectively, these patterns support genome replication via a deamination prone asymmetric strand-displacement mechanism as the fundamental driver of mutagenesis in mammalian DNA. Moreover, the distribution of mtDNA single nucleotide polymorphisms in humans and the distribution of bases in the mtDNA across vertebrate species mirror this gradient, indicating that replication-linked mutations are likely the primary source of inherited polymorphisms that, over evolutionary timescales, influences genome composition during speciation.