Comprehensive molecular characterization of mitochondrial genomes in human cancers

Single-mitochondrion sequencing uncovers distinct mutational patterns and heteroplasmy landscape in mouse astrocytes and neurons

BACKGROUND

Mitochondrial (mt) heteroplasmy can cause adverse biological consequences when deleterious mtDNA mutations accumulate disrupting "normal" mt-driven processes and cellular functions. To investigate the heteroplasmy of such mtDNA changes, we developed a moderate throughput mt isolation procedure to quantify the mt single-nucleotide variant (SNV) landscape in individual mouse neurons and astrocytes. In this study, we amplified mt-genomes from 1645 single mitochondria isolated from mouse single astrocytes and neurons to (1) determine the distribution and proportion of mt-SNVs as well as mutation pattern in specific target regions across the mt-genome, (2) assess differences in mtDNA SNVs between neurons and astrocytes, and (3) study co-segregation of variants in the mouse mtDNA.

RESULTS

(1) The data show that specific sites of the mt-genome are permissive to SNV presentation while others appear to be under stringent purifying selection. Nested hierarchical analysis at the levels of mitochondrion, cell, and mouse reveals distinct patterns of inter- and intra-cellular variation for mt-SNVs at different sites. (2) Further, differences in the SNV incidence were observed between mouse neurons and astrocytes for two mt-SNV 9027:G > A and 9419:C > T showing variation in the mutational propensity between these cell types. Purifying selection was observed in neurons as shown by the Ka/Ks statistic, suggesting that neurons are under stronger evolutionary constraint as compared to astrocytes. (3) Intriguingly, these data show strong linkage between the SNV sites at nucleotide positions 9027 and 9461.

CONCLUSIONS

This study suggests that segregation as well as clonal expansion of mt-SNVs is specific to individual genomic loci, which is important foundational data in understanding of heteroplasmy and disease thresholds for mutation of pathogenic variants.

Non-Coding RNAs of Mitochondrial Origin: Roles in Cell Division and Implications in Cancer

Non-coding RNAs (ncRNAs) are a heterogeneous group, in terms of structure and sequence length, consisting of RNA molecules that do not code for proteins. These ncRNAs have a central role in the regulation of gene expression and are virtually involved in every process analyzed, ensuring cellular homeostasis. Although, over the years, much research has focused on the characterization of non-coding transcripts of nuclear origin, improved bioinformatic tools and next-generation sequencing (NGS) platforms have allowed the identification of hundreds of ncRNAs transcribed from the mitochondrial genome (mt-ncRNA), including long non-coding RNA (lncRNA), circular RNA (circRNA), and microRNA (miR). Mt-ncRNAs have been described in diverse cellular processes such as mitochondrial proteome homeostasis and retrograde signaling; however, the function of the majority of mt-ncRNAs remains unknown. This review focuses on a subgroup of human mt-ncRNAs whose dysfunction is associated with both failures in cell cycle regulation, leading to defects in cell growth, cell proliferation, and apoptosis, and the development of tumor hallmarks, such as cell migration and metastasis formation, thus contributing to carcinogenesis and tumor development. Here we provide an overview of the mt-ncRNAs/cancer relationship that could help the future development of new biomedical applications in the field of oncology.

NeMu: a comprehensive pipeline for accurate reconstruction of neutral mutation spectra from evolutionary data

The recognized importance of mutational spectra in molecular evolution is yet to be fully exploited beyond human cancer studies and model organisms. The wealth of intraspecific polymorphism data in the GenBank repository, covering a broad spectrum of genes and species, presents an untapped opportunity for detailed mutational spectrum analysis. Existing methods fall short by ignoring intermediate substitutions on the inner branches of phylogenetic trees and lacking the capability for cross-species mutational comparisons. To address these challenges, we present the NeMu pipeline, available at https://nemu-pipeline.com, a tool grounded in phylogenetic principles designed to provide comprehensive and scalable analysis of mutational spectra. Utilizing extensive sequence data from numerous available genome projects, NeMu rapidly and accurately reconstructs the neutral mutational spectrum. This tool, facilitating the reconstruction of gene- and species-specific mutational spectra, contributes to a deeper understanding of evolutionary mechanisms across the broad spectrum of known species.

Decoding mitochondria's role in immunity and cancer therapy

The functions of mitochondria, including energy production and biomolecule synthesis, have been known for a long time. Given the rising incidence of cancer, the role of mitochondria in cancer has become increasingly popular. Activated by components released by mitochondria, various pathways interact with each other to induce immune responses to protect organisms from attack. However, mitochondria play dual roles in the progression of cancer. Abnormalities in proteins, which are the elementary structures of mitochondria, are closely linked with oncogenesis. Both the aberrant accumulation of intermediates and mutations in enzymes result in the generation and progression of cancer. Therefore, targeting mitochondria to treat cancer may be a new strategy. Several drugs aimed at inhibiting mutated enzymes and accumulated intermediates have been tested clinically. Here, we discuss the current understanding of mitochondria in cancer and the interactions between mitochondrial functions, immune responses, and oncogenesis. Furthermore, we discuss mitochondria as hopeful targets for cancer therapy, providing insights into the progression of future therapeutic strategies.

Mitochondrial genetics through the lens of single-cell multi-omics

Mitochondria carry their own genetic information encoding for a subset of protein-coding genes and translational machinery essential for cellular respiration and metabolism. Despite its small size, the mitochondrial genome, its natural genetic variation and molecular phenotypes have been challenging to study using bulk sequencing approaches, due to its variation in cellular copy number, non-Mendelian modes of inheritance and propensity for mutations. Here we highlight emerging strategies designed to capture mitochondrial genetic variation across individual cells for lineage tracing and studying mitochondrial genetics in primary human cells and clinical specimens. We review recent advances surrounding single-cell mitochondrial genome sequencing and its integration with functional genomic readouts, including leveraging somatic mitochondrial DNA mutations as clonal markers that can resolve cellular population dynamics in complex human tissues. Finally, we discuss how single-cell whole mitochondrial genome sequencing approaches can be utilized to investigate mitochondrial genetics and its contribution to cellular heterogeneity and disease.

Our current understanding of the biological impact of endometrial cancer mtDNA genome mutations and their potential use as a biomarker

Endometrial cancer (EC) is a devastating and common disease affecting women's health. The NCI Surveillance, Epidemiology, and End Results Program predicted that there would be >66,000 new cases in the United States and >13,000 deaths from EC in 2023, and EC is the sixth most common cancer among women worldwide. Regulation of mitochondrial metabolism plays a role in tumorigenesis. In proliferating cancer cells, mitochondria provide the necessary building blocks for biosynthesis of amino acids, lipids, nucleotides, and glucose. One mechanism causing altered mitochondrial activity is mitochondrial DNA (mtDNA) mutation. The polyploid human mtDNA genome is a circular double-stranded molecule essential to vertebrate life that harbors genes critical for oxidative phosphorylation plus mitochondrial-derived peptide genes. Cancer cells display aerobic glycolysis, known as the Warburg effect, which arises from the needs of fast-dividing cells and is characterized by increased glucose uptake and conversion of glucose to lactate. Solid tumors often contain at least one mtDNA substitution. Furthermore, it is common for cancer cells to harbor mixtures of wild-type and mutant mtDNA genotypes, known as heteroplasmy. Considering the increase in cancer cell energy demand, the presence of functionally relevant carcinogenesis-inducing or environment-adapting mtDNA mutations in cancer seems plausible. We review 279 EC tumor-specific mtDNA single nucleotide variants from 111 individuals from different studies. Many transition mutations indicative of error-prone DNA polymerase γ replication and C to U deamination events were present. We examine the spectrum of mutations and their heteroplasmy and discuss the potential biological impact of recurrent, non-synonymous, insertion, and deletion mutations. Lastly, we explore current EC treatments, exploiting cancer cell mitochondria for therapy and the prospect of using mtDNA variants as an EC biomarker.

Single-Mitochondrion Sequencing Uncovers Distinct Mutational Patterns and Heteroplasmy Landscape in Mouse Astrocytes and Neurons

Background

Mitochondrial (mt) heteroplasmy can cause adverse biological consequences when deleterious mtDNA mutations accumulate disrupting 'normal' mt-driven processes and cellular functions. To investigate the heteroplasmy of such mtDNA changes we developed a moderate throughput mt isolation procedure to quantify the mt single-nucleotide variant (SNV) landscape in individual mouse neurons and astrocytes In this study we amplified mt-genomes from 1,645 single mitochondria (mts) isolated from mouse single astrocytes and neurons to 1. determine the distribution and proportion of mt-SNVs as well as mutation pattern in specific target regions across the mt-genome, 2. assess differences in mtDNA SNVs between neurons and astrocytes, and 3. Study cosegregation of variants in the mouse mtDNA.

Results

1. The data show that specific sites of the mt-genome are permissive to SNV presentation while others appear to be under stringent purifying selection. Nested hierarchical analysis at the levels of mitochondrion, cell, and mouse reveals distinct patterns of inter- and intra-cellular variation for mt-SNVs at different sites. 2. Further, differences in the SNV incidence were observed between mouse neurons and astrocytes for two mt-SNV 9027:G>A and 9419:C>T showing variation in the mutational propensity between these cell types. Purifying selection was observed in neurons as shown by the Ka/Ks statistic, suggesting that neurons are under stronger evolutionary constraint as compared to astrocytes. 3. Intriguingly, these data show strong linkage between the SNV sites at nucleotide positions 9027 and 9461.

Conclusion

This study suggests that segregation as well as clonal expansion of mt-SNVs is specific to individual genomic loci, which is important foundational data in understanding of heteroplasmy and disease thresholds for mutation of pathogenic variants.

RESPIRATION DEFECTS LIMIT SERINE SYNTHESIS REQUIRED FOR LUNG CANCER GROWTH AND SURVIVAL

Mitochondrial function is important for both energetic and anabolic metabolism. Pathogenic mitochondrial DNA (mtDNA) mutations directly impact these functions, resulting in the detrimental consequences seen in human mitochondrial diseases. The role of pathogenic mtDNA mutations in human cancers is less clear; while pathogenic mtDNA mutations are observed in some cancer types, they are almost absent in others. We report here that the proofreading mutant DNA polymerase gamma ( PolG D256A ) induced a high mtDNA mutation burden in non-small-cell lung cancer (NSCLC), and promoted the accumulation of defective mitochondria, which is responsible for decreased tumor cell proliferation and viability and increased cancer survival. In NSCLC cells, pathogenic mtDNA mutations increased glycolysis and caused dependence on glucose. The glucose dependency sustained mitochondrial energetics but at the cost of a decreased NAD+/NADH ratio that inhibited de novo serine synthesis. Insufficient serine synthesis, in turn, impaired the downstream synthesis of GSH and nucleotides, leading to impaired tumor growth that increased cancer survival. Unlike tumors with intact mitochondrial function, NSCLC with pathogenic mtDNA mutations were sensitive to dietary serine and glycine deprivation. Thus, mitochondrial function in NSCLC is required specifically to sustain sufficient serine synthesis for nucleotide production and redox homeostasis to support tumor growth, explaining why these cancers preserve functional mtDNA.

In brief

High mtDNA mutation burden in non-small-cell lung cancer (NSCLC) leads to the accumulation of respiration-defective mitochondria and dependency on glucose and glycolytic metabolism. Defective respiratory metabolism causes a massive accumulation of cytosolic nicotinamide adenine dinucleotide + hydrogen (NADH), which impedes serine synthesis and, thereby, glutathione (GSH) and nucleotide synthesis, leading to impaired tumor growth and increased survival.

Highlights

Proofreading mutations in Polymerase gamma led to a high burden of mitochondrial DNA mutations, promoting the accumulation of mitochondria with respiratory defects in NSCLC.Defective respiration led to reduced proliferation and viability of NSCLC cells increasing survival to cancer.Defective respiration caused glucose dependency to fuel elevated glycolysis.Altered glucose metabolism is associated with high NADH that limits serine synthesis, leading to impaired GSH and nucleotide production.Mitochondrial respiration defects sensitize NSCLC to dietary serine/glycine starvation, further increasing survival.

Abstract Figure

DNA mismatch and damage patterns revealed by single-molecule sequencing

Mutations accumulate in the genome of every cell of the body throughout life, causing cancer and other diseases1,2. Most mutations begin as nucleotide mismatches or damage in one of the two strands of the DNA before becoming double-strand mutations if unrepaired or misrepaired3,4. However, current DNA-sequencing technologies cannot accurately resolve these initial single-strand events. Here we develop a single-molecule, long-read sequencing method (Hairpin Duplex Enhanced Fidelity sequencing (HiDEF-seq)) that achieves single-molecule fidelity for base substitutions when present in either one or both DNA strands. HiDEF-seq also detects cytosine deamination-a common type of DNA damage-with single-molecule fidelity. We profiled 134 samples from diverse tissues, including from individuals with cancer predisposition syndromes, and derive from them single-strand mismatch and damage signatures. We find correspondences between these single-strand signatures and known double-strand mutational signatures, which resolves the identity of the initiating lesions. Tumours deficient in both mismatch repair and replicative polymerase proofreading show distinct single-strand mismatch patterns compared to samples that are deficient in only polymerase proofreading. We also define a single-strand damage signature for APOBEC3A. In the mitochondrial genome, our findings support a mutagenic mechanism occurring primarily during replication. As double-strand DNA mutations are only the end point of the mutation process, our approach to detect the initiating single-strand events at single-molecule resolution will enable studies of how mutations arise in a variety of contexts, especially in cancer and ageing.

Rational combination platform trial design for children and young adults with diffuse midline glioma: A report from PNOC

Background Diffuse midline glioma (DMG) is a devastating pediatric brain tumor unresponsive to hundreds of clinical trials. Approximately 80% of DMGs harbor H3K27M oncohistones, which reprogram the epigenome to increase the metabolic profile of the tumor cells. Methods We have previously shown preclinical efficacy of targeting both oxidative phosphorylation and glycolysis through treatment with ONC201, which activates the mitochondrial protease ClpP, and paxalisib, which inhibits PI3K/mTOR, respectively. Results ONC201 and paxalisib combination treatment aimed at inducing metabolic distress led to the design of the first DMG-specific platform trial PNOC022 (NCT05009992). Conclusions Here, we expand on the PNOC022 rationale and discuss various considerations, including liquid biome, microbiome, and genomic biomarkers, quality-of-life endpoints, and novel imaging modalities, such that we offer direction on future clinical trials in DMG.

Single-cell mtDNA dynamics in tumors is driven by coregulation of nuclear and mitochondrial genomes

The extent of cell-to-cell variation in tumor mitochondrial DNA (mtDNA) copy number and genotype, and the phenotypic and evolutionary consequences of such variation, are poorly characterized. Here we use amplification-free single-cell whole-genome sequencing (Direct Library Prep (DLP+)) to simultaneously assay mtDNA copy number and nuclear DNA (nuDNA) in 72,275 single cells derived from immortalized cell lines, patient-derived xenografts and primary human tumors. Cells typically contained thousands of mtDNA copies, but variation in mtDNA copy number was extensive and strongly associated with cell size. Pervasive whole-genome doubling events in nuDNA associated with stoichiometrically balanced adaptations in mtDNA copy number, implying that mtDNA-to-nuDNA ratio, rather than mtDNA copy number itself, mediated downstream phenotypes. Finally, multimodal analysis of DLP+ and single-cell RNA sequencing identified both somatic loss-of-function and germline noncoding variants in mtDNA linked to heteroplasmy-dependent changes in mtDNA copy number and mitochondrial transcription, revealing phenotypic adaptations to disrupted nuclear/mitochondrial balance.

Approaches and challenges in identifying, quantifying, and manipulating dynamic mitochondrial genome variations

Mitochondria, the cellular powerhouses, possess their own unique genetic system, including replication, transcription, and translation. Studying these processes is crucial for comprehending mitochondrial disorders, energy production, and their related diseases. Over the past decades, various approaches have been applied in detecting and quantifying mitochondrial genome variations with also the purpose of manipulation of mitochondria or mitochondrial genome for therapeutics. Understanding the scope and limitations of above strategies is not only fundamental to the understanding of basic biology but also critical for exploring disease-related novel target(s), as well to develop innovative therapies. Here, this review provides an overview of different tools and techniques for accurate mitochondrial genome variations identification, quantification, and discuss novel strategies for the manipulation of mitochondria to develop innovative therapeutic interventions, through combining the insights gained from the study of mitochondrial genetics with ongoing single cell omics combined with advanced single molecular tools.

Single-cell analysis reveals context-dependent, cell-level selection of mtDNA

Heteroplasmy occurs when wild-type and mutant mitochondrial DNA (mtDNA) molecules co-exist in single cells1. Heteroplasmy levels change dynamically in development, disease and ageing2,3, but it is unclear whether these shifts are caused by selection or drift, and whether they occur at the level of cells or intracellularly. Here we investigate heteroplasmy dynamics in dividing cells by combining precise mtDNA base editing (DdCBE)4 with a new method, SCI-LITE (single-cell combinatorial indexing leveraged to interrogate targeted expression), which tracks single-cell heteroplasmy with ultra-high throughput. We engineered cells to have synonymous or nonsynonymous complex I mtDNA mutations and found that cell populations in standard culture conditions purge nonsynonymous mtDNA variants, whereas synonymous variants are maintained. This suggests that selection dominates over simple drift in shaping population heteroplasmy. We simultaneously tracked single-cell mtDNA heteroplasmy and ancestry, and found that, although the population heteroplasmy shifts, the heteroplasmy of individual cell lineages remains stable, arguing that selection acts at the level of cell fitness in dividing cells. Using these insights, we show that we can force cells to accumulate high levels of truncating complex I mtDNA heteroplasmy by placing them in environments where loss of biochemical complex I activity has been reported to benefit cell fitness. We conclude that in dividing cells, a given nonsynonymous mtDNA heteroplasmy can be harmful, neutral or even beneficial to cell fitness, but that the 'sign' of the effect is wholly dependent on the environment.

A PCR-independent approach for mtDNA enrichment and next-generation sequencing: comprehensive evaluation and clinical application

BACKGROUND

Sequencing the mitochondrial genome has been increasingly important for the investigation of primary mitochondrial diseases (PMD) and mitochondrial genetics. To overcome the limitations originating from PCR-based mtDNA enrichment, we set out to develop and evaluate a PCR-independent approach in this study, named Pime-Seq (PCR-independent mtDNA enrichment and next generation Sequencing).

RESULTS

By using the optimized mtDNA enrichment procedure, the mtDNA reads ratio reached 88.0 ± 7.9% in the sequencing library when applied on human PBMC samples. We found the variants called by Pime-Seq were highly consistent among technical repeats. To evaluate the accuracy and reliability of this method, we compared Pime-Seq with lrPCR based NGS by performing both methods simultaneously on 45 samples, yielding 1677 concordant variants, as well as 146 discordant variants with low-level heteroplasmic fraction, in which Pime-Seq showed higher reliability. Furthermore, we applied Pime-Seq on 4 samples of PMD patients retrospectively, and successfully detected all the pathogenic mtDNA variants. In addition, we performed a prospective study on 192 apparently healthy pregnant women during prenatal screening, in which Pime-Seq identified pathogenic mtDNA variants in 4 samples, providing extra information for better health monitoring in these cases.

CONCLUSIONS

Pime-Seq can obtain highly enriched mtDNA in a PCR-independent manner for high quality and reliable mtDNA deep-sequencing, which provides us an effective and promising tool for detecting mtDNA variants for both clinical and research purposes.

IFNα-induced BST2+ tumor-associated macrophages facilitate immunosuppression and tumor growth in pancreatic cancer by ERK-CXCL7 signaling

Pancreatic ductal adenocarcinoma (PDAC) features an immunosuppressive tumor microenvironment (TME) that resists immunotherapy. Tumor-associated macrophages, abundant in the TME, modulate T cell responses. Bone marrow stromal antigen 2-positive (BST2+) macrophages increase in KrasG12D/+; Trp53R172H/+; Pdx1-Cre mouse models during PDAC progression. However, their role in PDAC remains elusive. Our findings reveal a negative correlation between BST2+ macrophage levels and PDAC patient prognosis. Moreover, an increased ratio of exhausted CD8+ T cells is observed in tumors with up-regulated BST2+ macrophages. Mechanistically, BST2+ macrophages secrete CXCL7 through the ERK pathway and bind with CXCR2 to activate the AKT/mTOR pathway, promoting CD8+ T cell exhaustion. The combined blockade of CXCL7 and programmed death-ligand 1 successfully decelerates tumor growth. Additionally, cGAS-STING pathway activation in macrophages induces interferon (IFN)α synthesis leading to BST2 overexpression in the PDAC TME. This study provides insights into IFNα-induced BST2+ macrophages driving an immune-suppressive TME through ERK-CXCL7 signaling to regulate CD8+ T cell exhaustion in PDAC.

Mitochondria: A source of potential biomarkers for non-communicable diseases

Mitochondria, as an endosymbiont of eukaryotic cells, controls multiple cellular activities, including respiration, reactive oxygen species production, fatty acid synthesis, and death. Though the majority of functional mitochondrial proteins are translated through a nucleus-controlled process, very few of them (∼10%) are translated within mitochondria through their own machinery. Germline and somatic mutations in mitochondrial and nuclear DNA significantly impact mitochondrial homeostasis and function. Such modifications disturbing mitochondrial biogenesis, metabolism, or mitophagy eventually resulted in cellular pathophysiology. In this chapter, we discussed the impact of mitochondria and its dysfunction on several non-communicable diseases like cancer, diabetes, neurodegenerative, and cardiovascular problems. Mitochondrial dysfunction and its outcome could be screened by currently available omics-based techniques, flow cytometry, and high-resolution imaging. Such characterization could be evaluated as potential biomarkers to assess the disease burden and prognosis.

Mitochondrial-Stem Cell Connection: Providing Additional Explanations for Understanding Cancer

The cancer paradigm is generally based on the somatic mutation model, asserting that cancer is a disease of genetic origin. The mitochondrial-stem cell connection (MSCC) proposes that tumorigenesis may result from an alteration of the mitochondria, specifically a chronic oxidative phosphorylation (OxPhos) insufficiency in stem cells, which forms cancer stem cells (CSCs) and leads to malignancy. Reviewed evidence suggests that the MSCC could provide a comprehensive understanding of all the different stages of cancer. The metabolism of cancer cells is altered (OxPhos insufficiency) and must be compensated by using the glycolysis and the glutaminolysis pathways, which are essential to their growth. The altered mitochondria regulate the tumor microenvironment, which is also necessary for cancer evolution. Therefore, the MSCC could help improve our understanding of tumorigenesis, metastases, the efficiency of standard treatments, and relapses.

Mitochondrial DNA copy number and cancer risks: A comprehensive Mendelian randomization analysis

Mitochondrial DNA plays a critical role in the pathophysiology of cancer. However, the associations between mitochondrial DNA copy number (mtDNA-CN) and cancer risk are controversial. Mendelian randomization (MR) analyses were performed using three independent instrumental variables (IVs) to explore potential associations between mtDNA-CN and 20 types of cancer. The three sets of IVs were primarily obtained from participants in the UK Biobank and the Cohorts for Heart and Aging Research in Genomic Epidemiology consortium using different methods. The outcome data of cancers were investigated using summary statistics from the FinnGen cohort. The potential causal associations were evaluated using the MR-Egger regression, weighted median, inverse-variance weighted (IVW), and weighted mode methods. The robustness of IVW estimates was validated using leave-one-out sensitivity analysis. Additionally, a meta-analysis was conducted to pool results from three sets of IVs. The results revealed that genetically predicted mtDNA-CN was not associated with cancer risk (odds ratio = 1.02; 95% confidence interval: 0.95-1.10). Subgroup analyses indicated no causal association between mtDNA-CN and breast, lung, prostate, skin, colorectal, gastric, liver, cervical uteri, esophageal, thyroid, bladder, pancreas, kidney, corpus uteri, ovary, brain, larynx, and anus cancers. It was observed that mtDNA-CN was associated with lip, oral cavity, and testis cancers. However, these results should be interpreted with caution because a small number of patients with lip and oral cavity or testis cancers were included. The comprehensive MR analysis demonstrated that mtDNA-CN is not a suitable biomarker for tumor risk assessment.

Leveraging new methods for comprehensive characterization of mitochondrial DNA in esophageal squamous cell carcinoma

BACKGROUND

Mitochondria play essential roles in tumorigenesis; however, little is known about the contribution of mitochondrial DNA (mtDNA) to esophageal squamous cell carcinoma (ESCC). Whole-genome sequencing (WGS) is by far the most efficient technology to fully characterize the molecular features of mtDNA; however, due to the high redundancy and heterogeneity of mtDNA in regular WGS data, methods for mtDNA analysis are far from satisfactory.

METHODS

Here, we developed a likelihood-based method dMTLV to identify low-heteroplasmic mtDNA variants. In addition, we described fNUMT, which can simultaneously detect non-reference nuclear sequences of mitochondrial origin (non-ref NUMTs) and their derived artifacts. Using these new methods, we explored the contribution of mtDNA to ESCC utilizing the multi-omics data of 663 paired tumor-normal samples.

RESULTS

dMTLV outperformed the existing methods in sensitivity without sacrificing specificity. The verification using Nanopore long-read sequencing data showed that fNUMT has superior specificity and more accurate breakpoint identification than the current methods. Leveraging the new method, we identified a significant association between the ESCC overall survival and the ratio of mtDNA copy number of paired tumor-normal samples, which could be potentially explained by the differential expression of genes enriched in pathways related to metabolism, DNA damage repair, and cell cycle checkpoint. Additionally, we observed that the expression of CBWD1 was downregulated by the non-ref NUMTs inserted into its intron region, which might provide precursor conditions for the tumor cells to adapt to a hypoxic environment. Moreover, we identified a strong positive relationship between the number of mtDNA truncating mutations and the contribution of signatures linked to tumorigenesis and treatment response.

CONCLUSIONS

Our new frameworks promote the characterization of mtDNA features, which enables the elucidation of the landscapes and roles of mtDNA in ESCC essential for extending the current understanding of ESCC etiology. dMTLV and fNUMT are freely available from https://github.com/sunnyzxh/dMTLV and https://github.com/sunnyzxh/fNUMT , respectively.

Advancing nanotechnology for neoantigen-based cancer theranostics

Neoantigens play a pivotal role in the field of tumour therapy, encompassing the stimulation of anti-tumour immune response and the enhancement of tumour targeting capability. Nonetheless, numerous factors directly influence the effectiveness of neoantigens in bolstering anti-tumour immune responses, including neoantigen quantity and specificity, uptake rates by antigen-presenting cells (APCs), residence duration within the tumour microenvironment (TME), and their ability to facilitate the maturation of APCs for immune response activation. Nanotechnology assumes a significant role in several aspects, including facilitating neoantigen release, promoting neoantigen delivery to antigen-presenting cells, augmenting neoantigen uptake by dendritic cells, shielding neoantigens from protease degradation, and optimizing interactions between neoantigens and the immune system. Consequently, the development of nanotechnology synergistically enhances the efficacy of neoantigens in cancer theranostics. In this review, we provide an overview of neoantigen sources, the mechanisms of neoantigen-induced immune responses, and the evolution of precision neoantigen-based nanomedicine. This encompasses various therapeutic modalities, such as neoantigen-based immunotherapy, phototherapy, radiotherapy, chemotherapy, chemodynamic therapy, and other strategies tailored to augment precision in cancer therapeutics. We also discuss the current challenges and prospects in the application of neoantigen-based precision nanomedicine, aiming to expedite its clinical translation.

Mitochondrial DNA mutations drive aerobic glycolysis to enhance checkpoint blockade response in melanoma

The mitochondrial genome (mtDNA) encodes essential machinery for oxidative phosphorylation and metabolic homeostasis. Tumor mtDNA is among the most somatically mutated regions of the cancer genome, but whether these mutations impact tumor biology is debated. We engineered truncating mutations of the mtDNA-encoded complex I gene, Mt-Nd5, into several murine models of melanoma. These mutations promoted a Warburg-like metabolic shift that reshaped tumor microenvironments in both mice and humans, consistently eliciting an anti-tumor immune response characterized by loss of resident neutrophils. Tumors bearing mtDNA mutations were sensitized to checkpoint blockade in a neutrophil-dependent manner, with induction of redox imbalance being sufficient to induce this effect in mtDNA wild-type tumors. Patient lesions bearing >50% mtDNA mutation heteroplasmy demonstrated a response rate to checkpoint blockade that was improved by ~2.5-fold over mtDNA wild-type cancer. These data nominate mtDNA mutations as functional regulators of cancer metabolism and tumor biology, with potential for therapeutic exploitation and treatment stratification.

Increases in computationally predicted deleterious variants of unknown significance and sperm mtDNA copy numbers may be associated with semen quality

BACKGROUND

Mitochondria are essential for sperm motility because they provide the energy required for the movement. Changes in sperm mtDNA, such as point mutations, large-scale deletions, or copy number variations, may interfere with ATP production and reduce sperm motility. However, it is not clear if changes in mtDNA are linked to semen quality.

OBJECTIVES

To explore the association between sperm mitochondrial DNA (mtDNA) changes and semen quality.

MATERIALS AND METHODS

Sixty-five oligo and/or astheno and/or terato patients (O/A/T) patients and 41 controls were recruited from couples undergoing assisted reproduction. Semen and blood samples were collected from the same individual on the day of oocyte retrieval to extract, isolate and purify mtDNA for next-generation sequencing. mtDNA copy numbers were assessed in 64 patient and 39 control sperm DNA samples using quantitative real-time PCR. The 4977 bp deletion was assessed in 20 patient and 20 control sperm DNA samples using polymerase chain reaction.

RESULTS

The mtDNA of patients was more likely to carry pathogenic variants or variants of unknown significance (VUSs) (P = 0.091) with higher heteroplasmy levels (P < 0.05) than that of controls. Interestingly, 33.85% of O/A/T patients (22 out of 65) lacked unique variants in their spermatozoa. but presented an exceptionally high mtDNA copy number (P < 0.0001). Moreover, we observed a decrease in the heteroplasmy level of common mtDNA variants shared by somatic and gamete cells (P < 0.0001) and the emergence of a very large number of de novo mtDNA variants with low-level heteroplasmy in spermatozoa.

DISCUSSION AND CONCLUSION

The increases in the number of computationally predicted deleterious VUS and mtDNA copies in spermatozoa may be associated with semen quality. Exposure to environmental mutation pressure that causes novel mtDNA variants with low-level heteroplasmy may occur during spermatogenesis. Furthermore, when a certain harmful threshold is reached, male germ cells may degrade mtDNA with mutations and replicate the correct mtDNA sequence to maintain the mitochondrial function in spermatozoa.

Forward Genetic Screens Identify Mechanisms of Resistance to Small-Molecule Lactate Dehydrogenase Inhibitors

Altered metabolism is a hallmark of cancer; however, it has been difficult to specifically target metabolism in cancer for therapeutic benefit. Cancers with genetically defined defects in metabolic enzymes constitute a subset of cancers where targeting metabolism is potentially accessible. Hürthle cell carcinoma of the thyroid (HTC) tumors frequently harbor deleterious mitochondrial DNA (mtDNA) mutations in subunits of complex I of the mitochondrial electron transport chain (ETC). Previous work has shown that HTC models with deleterious mtDNA mutations exhibit mitochondrial ETC defects that expose lactate dehydrogenase (LDH) as a therapeutic vulnerability. Here, we performed forward genetic screens to identify mechanisms of resistance to small-molecule LDH inhibitors. We identified two distinct mechanisms of resistance: upregulation of an LDH isoform and a compound-specific resistance mutation. Using these tools, we demonstrate that the anticancer activity of LDH inhibitors in cell line and xenograft models of complex I mutant HTC is through on-target LDH inhibition.

The known unknowns of mitochondrial carcinogenesis: de novo NUMTs and intercellular mitochondrial transfer

The translocation of mitochondrial DNA (mtDNA) sequences into the nuclear genome, resulted in the occurrence of nuclear sequences of mitochondrial origin (NUMTs) which can be detected in nearly all sequenced eukaryotes. However, de novo mtDNA insertions can contribute to the development of pathological conditions including cancer. Recent data indicate that de novo mtDNA translocation into chromosomes can occur due to genotoxic influence of DNA double-strand break-inducing environmental mutagens. This confirms the hypothesis of the involvement of genome instability in the occurrence of mtDNA fragments in chromosomes. Mounting evidence indicates that mitochondria can be transferred from normal cells to cancer cells and recover cellular respiration. These exchanged mitochondria can facilitate cancer progression and metastasis. This review article provides a comprehensive overview of the potential carcinogenicity of mtDNA insertions, and the relevance of mtDNA escape in cancer progression, metastasis, and treatment resistance in humans. Potential molecular targets involved in mtDNA escape and exchange of mitochondria that can be of possible clinical benefits are presented and discussed. Understanding these processes could lead to improved diagnostic approaches, novel therapeutic strategies, and a deeper understanding of the intricate relationship between mitochondria, nuclear DNA, and cancer biology.

The molecular mechanism of aging and the role in neurodegenerative diseases

Aging is a complex and inevitable biological process affected by a combination of external environmental and genetic factors. Humans are currently living longer than ever before, accompanied with aging-related alterations such as diminished autophagy, decreased immunological function, mitochondrial malfunction, stem cell failure, accumulation of somatic and mitochondrial DNA mutations, loss of telomere, and altered nutrient metabolism. Aging leads to a decline in body functions and age-related diseases, for example, Alzheimer's disease, which adversely affects human health and longevity. The quality of life of the elderly is greatly diminished by the increase in their life expectancy rather than healthy life expectancy. With the rise in the age of the global population, aging and related diseases have become the focus of attention worldwide. In this review, we discuss several major mechanisms of aging, including DNA damage and repair, free radical oxidation, telomeres and telomerase, mitochondrial damage, inflammation, and their role in neurodegenerative diseases to provide a reference for the prevention of aging and its related diseases.

Effects of hexavalent chromium on mitochondria and their implications in carcinogenesis

Hexavalent chromium (Cr(VI)) is a well-known occupational and environmental human carcinogen. The cellular effect of Cr(VI) is complex and often nonspecific due to its ability to modulate multiple cellular targets. The toxicity of Cr(VI) is strongly linked to the generation of reactive oxygen species (ROS) during its reduction process. ROS can cause oxidation of cellular macromolecules, such as proteins, lipids, and DNA, thereby altering their functions. A major genotoxic effect of Cr(VI) that contributes to carcinogenesis is the formation of DNA adducts, which can lead to DNA damage. Modulations of cellular signaling pathways and epigenetics may also contribute to the carcinogenic effects of Cr(VI). Cr(VI) has a major impact on many aspects of mitochondrial biology, including oxidative phosphorylation, mitophagy, and mitochondrial biogenesis. These effects have the potential to alter the trajectory of Cr(VI)-induced carcinogenic process. This perspective article summarizes current understandings of the effect of Cr(VI) on mitochondria and discusses the future directions of research in this area, particularly with regard to carcinogenesis.

Yeast EndoG prevents genome instability by degrading cytoplasmic DNA

In metazoans release of mitochondrial DNA or retrotransposon cDNA to cytoplasm can cause sterile inflammation and disease 1. Cytoplasmic nucleases degrade these DNA species to limit inflammation 2,3. It remains unknown whether degradation these DNA also prevents nuclear genome instability. To address this question, we decided to identify the nuclease regulating transfer of these cytoplasmic DNA species to the nucleus. We used an amplicon sequencing-based method in yeast enabling analysis of millions of DSB repair products. Nuclear mtDNA (NUMTs) and retrotransposon cDNA insertions increase dramatically in nondividing stationary phase cells. Yeast EndoG (Nuc1) nuclease limits insertions of cDNA and transfer of very long mtDNA (>10 kb) that forms unstable circles or rarely insert in the genome, but it promotes formation of short NUMTs (~45-200 bp). Nuc1 also regulates transfer of cytoplasmic DNA to nucleus in aging or during meiosis. We propose that Nuc1 preserves genome stability by degrading retrotransposon cDNA and long mtDNA, while short NUMTs can originate from incompletely degraded mtDNA. This work suggests that nucleases eliminating cytoplasmic DNA play a role in preserving genome stability.

Reanalysis of mtDNA mutations of human primordial germ cells (PGCs) reveals NUMT contamination and suggests that selection in PGCs may be positive

The resilience of the mitochondrial genome (mtDNA) to a high mutational pressure depends, in part, on negative purifying selection in the germline. A paradigm in the field has been that such selection, at least in part, takes place in primordial germ cells (PGCs). Specifically, Floros et al. (Nature Cell Biology 20: 144-51) reported an increase in the synonymity of mtDNA mutations (a sign of purifying selection) between early-stage and late-stage PGCs. We re-analyzed Floros' et al. data and determined that their mutational dataset was significantly contaminated with single nucleotide variants (SNVs) derived from a nuclear sequence of mtDNA origin (NUMT) located on chromosome 5. Contamination was caused by co-amplification of the NUMT sequence by cross-specific PCR primers. Importantly, when we removed NUMT-derived SNVs, the evidence of purifying selection was abolished. In addition to bulk PGCs, Floros et al. reported the analysis of single-cell late-stage PGCs, which were amplified with different sets of PCR primers that cannot amplify the NUMT sequence. Accordingly, there were no NUMT-derived SNVs among single PGC mutations. Interestingly, single PGC mutations show adecreaseof synonymity with increased intracellular mutant fraction. More specifically, nonsynonymous mutations show faster intracellular genetic drift towards higher mutant fraction than synonymous ones. This pattern is incompatible with predominantly negative selection. This suggests that germline selection of mtDNA mutations is a complex phenomenon and that the part of this process that takes place in PGCs may be predominantly positive. However counterintuitive, positive germline selection of detrimental mtDNA mutations has been reported previously andpotentially may be evolutionarily advantageous.

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

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

Mitochondrial response of glioma cells to temozolomide

Metabolic adaptations are central for carcinogenesis and response to therapy, but little is known about the contribution of mitochondrial dynamics to the response of glioma cells to the standard treatment with temozolomide (TMZ). Glioma cells responded to TMZ with mitochondrial mass increased and the production of round structures of dysfunctional mitochondria. At single-cell level, asymmetric mitosis contributed to the heterogeneity of mitochondrial levels. It affected the fitness of cells in control and treated condition, indicating that the mitochondrial levels are relevant for glioma cell fitness in the presence of TMZ.

MRPL48 is a novel prognostic and predictive biomarker of hepatocellular carcinoma

BACKGROUND

Hepatocellular carcinoma (HCC) is one of the most prevalent forms of cancer and poses a threat to the health and survival of humans. Mitochondrial ribosomal protein L48 (MRPL48) belongs to the mitochondrial ribosomal protein family, which participates in energy production. Studies have shown that MRPL48 can predict osteosarcoma incidence and prognosis, as well as promotes colorectal cancer progression. However, the role of MRPL48 in HCC remains unknown.

METHODS

TCGA, GEO, HCCDB, CPTAC, SMART, UALCAN, Kaplan-Meier plotter, cBioPortal, and MethSurv were performed for bioinformatics purposes. Quantitative RT-PCR, immunoblotting, and functional studies were conducted to validate the methodology in vitro.

RESULTS

MRPL48 was greatly overexpressed in HCC tissues, compared with healthy tissue, which was subsequently demonstrated in vitro as well. The survival and regression analyses showed that MRPL48 expression is of significant clinical prognostic value in HCC. The ROC curve and nomogram analysis indicated that MRPL48 is a powerful predictor of HCC. MRPL48 methylation was adversely associated with the expression of MRPL48, and patients with a low level of methylation had poorer overall survival than those with a high level of methylation. GSEA showed that the expression of the MRPL48 was correlated with Resolution of Sister Chromatid Cohesion, Mitotic Prometaphase, Retinoblastoma Gene in Cancer, RHO Gtpases Activate Formins, Mitotic Metaphase and Anaphase, and Cell Cycle Checkpoints. An analysis of immune cell infiltration showed a significant association between MRPL48 and immune cell infiltration subsets, which impacted the survival of HCC patients. Additionally, MRPL48 knockdown reduced HCC cell proliferation, migration, and invasion in vitro.

CONCLUSIONS

We demonstrated that MRPL48 expression may be associated with HCC development and prognosis. These findings may open up new research directions and opportunities for the development of HCC treatments.

Yeast EndoG prevents genome instability by degrading cytoplasmic DNA

In metazoans release of mitochondrial DNA or retrotransposon cDNA to cytoplasm can cause sterile inflammation and disease. Cytoplasmic nucleases degrade these DNA species to limit inflammation. It remains unknown whether degradation these DNA also prevents nuclear genome instability. To address this question, we decided to identify the nuclease regulating transfer of these cytoplasmic DNA species to the nucleus. We used an amplicon sequencing-based method in yeast enabling analysis of millions of DSB repair products. Nu clear mt DNA (NUMTs) and retrotransposon cDNA insertions increase dramatically in nondividing stationary phase cells. Yeast EndoG (Nuc1) nuclease limits insertions of cDNA and transfer of very long mtDNA (>10 kb) that forms unstable circles or rarely insert in the genome, but it promotes formation of short NUMTs (∼45-200 bp). Nuc1 also regulates transfer of cytoplasmic DNA to nucleus in aging or during meiosis. We propose that Nuc1 preserves genome stability by degrading retrotransposon cDNA and long mtDNA, while short NUMTs can originate from incompletely degraded mtDNA. This work suggests that nucleases eliminating cytoplasmic DNA play a role in preserving genome stability.

Mitochondrial DNA Variants at Low-Level Heteroplasmy and Decreased Copy Numbers in Chronic Kidney Disease (CKD) Tissues with Kidney Cancer

The human mitochondrial genome (mtDNA) is a circular DNA molecule with a length of 16.6 kb, which contains a total of 37 genes. Somatic mtDNA mutations accumulate with age and environmental exposure, and some types of mtDNA variants may play a role in carcinogenesis. Recent studies observed mtDNA variants not only in kidney tumors but also in adjacent kidney tissues, and mtDNA dysfunction results in kidney injury, including chronic kidney disease (CKD). To investigate whether a relationship exists between heteroplasmic mtDNA variants and kidney function, we performed ultra-deep sequencing (30,000×) based on long-range PCR of DNA from 77 non-tumor kidney tissues of kidney cancer patients with CKD (stages G1 to G5). In total, this analysis detected 697 single-nucleotide variants (SNVs) and 504 indels as heteroplasmic (0.5% ≤ variant allele frequency (VAF) < 95%), and the total number of detected SNVs/indels did not differ between CKD stages. However, the number of deleterious low-level heteroplasmic variants (pathogenic missense, nonsense, frameshift and tRNA) significantly increased with CKD progression (p < 0.01). In addition, mtDNA copy numbers (mtDNA-CNs) decreased with CKD progression (p < 0.001). This study demonstrates that mtDNA damage, which affects mitochondrial genes, may be involved in reductions in mitochondrial mass and associated with CKD progression and kidney dysfunction.

Resolving therapy resistance mechanisms in multiple myeloma by multiomics subclone analysis

Intratumor heterogeneity as a clinical challenge becomes most evident after several treatment lines, when multidrug-resistant subclones accumulate. To address this challenge, the characterization of resistance mechanisms at the subclonal level is key to identify common vulnerabilities. In this study, we integrate whole-genome sequencing, single-cell (sc) transcriptomics (scRNA sequencing), and chromatin accessibility (scATAC sequencing) together with mitochondrial DNA mutations to define subclonal architecture and evolution for longitudinal samples from 15 patients with relapsed or refractory multiple myeloma. We assess transcriptomic and epigenomic changes to resolve the multifactorial nature of therapy resistance and relate it to the parallel occurrence of different mechanisms: (1) preexisting epigenetic profiles of subclones associated with survival advantages, (2) converging phenotypic adaptation of genetically distinct subclones, and (3) subclone-specific interactions of myeloma and bone marrow microenvironment cells. Our study showcases how an integrative multiomics analysis can be applied to track and characterize distinct multidrug-resistant subclones over time for the identification of molecular targets against them.

Centuries of genome instability and evolution in soft-shell clam, Mya arenaria, bivalve transmissible neoplasia

Transmissible cancers are infectious parasitic clones that metastasize to new hosts, living past the death of the founder animal in which the cancer initiated. We investigated the evolutionary history of a cancer lineage that has spread though the soft-shell clam (Mya arenaria) population by assembling a chromosome-scale soft-shell clam reference genome and characterizing somatic mutations in transmissible cancer. We observe high mutation density, widespread copy-number gain, structural rearrangement, loss of heterozygosity, variable telomere lengths, mitochondrial genome expansion and transposable element activity, all indicative of an unstable cancer genome. We also discover a previously unreported mutational signature associated with overexpression of an error-prone polymerase and use this to estimate the lineage to be >200 years old. Our study reveals the ability for an invertebrate cancer lineage to survive for centuries while its genome continues to structurally mutate, likely contributing to the evolution of this lineage as a parasitic cancer.

Somatic evolution of marine transmissible leukemias in the common cockle, Cerastoderma edule

Transmissible cancers are malignant cell lineages that spread clonally between individuals. Several such cancers, termed bivalve transmissible neoplasia (BTN), induce leukemia-like disease in marine bivalves. This is the case of BTN lineages affecting the common cockle, Cerastoderma edule, which inhabits the Atlantic coasts of Europe and northwest Africa. To investigate the evolution of cockle BTN, we collected 6,854 cockles, diagnosed 390 BTN tumors, generated a reference genome and assessed genomic variation across 61 tumors. Our analyses confirmed the existence of two BTN lineages with hemocytic origins. Mitochondrial variation revealed mitochondrial capture and host co-infection events. Mutational analyses identified lineage-specific signatures, one of which likely reflects DNA alkylation. Cytogenetic and copy number analyses uncovered pervasive genomic instability, with whole-genome duplication, oncogene amplification and alkylation-repair suppression as likely drivers. Satellite DNA distributions suggested ancient clonal origins. Our study illuminates long-term cancer evolution under the sea and reveals tolerance of extreme instability in neoplastic genomes.

13C tracer analysis reveals the landscape of metabolic checkpoints in human CD8+ T cell differentiation and exhaustion

Introduction

Naïve T cells remain in an actively maintained state of quiescence until activation by antigenic signals, upon which they start to proliferate and generate effector cells to initiate a functional immune response. Metabolic reprogramming is essential to meet the biosynthetic demands of the differentiation process, and failure to do so can promote the development of hypofunctional exhausted T cells.

Methods

Here we used 13C metabolomics and transcriptomics to study the metabolism of CD8+ T cells in their complete course of differentiation from naïve over stem-like memory to effector cells and in exhaustion-inducing conditions.

Results

The quiescence of naïve T cells was evident in a profound suppression of glucose oxidation and a decreased expression of ENO1, downstream of which no glycolytic flux was detectable. Moreover, TCA cycle activity was low in naïve T cells and associated with a downregulation of SDH subunits. Upon stimulation and exit from quiescence, the initiation of cell growth and proliferation was accompanied by differential expression of metabolic enzymes and metabolic reprogramming towards aerobic glycolysis with high rates of nutrient uptake, respiration and lactate production. High flux in anabolic pathways imposed a strain on NADH homeostasis, which coincided with engagement of the proline cycle for mitochondrial redox shuttling. With acquisition of effector functions, cells increasingly relied on glycolysis as opposed to oxidative phosphorylation, which was, however, not linked to changes in mitochondrial abundance. In exhaustion, decreased effector function concurred with a reduction in mitochondrial metabolism, glycolysis and amino acid import, and an upregulation of quiescence-associated genes, TXNIP and KLF2, and the T cell suppressive metabolites succinate and itaconate.

Discussion

Overall, these results identify multiple metabolic features that regulate quiescence, proliferation and effector function, but also exhaustion of CD8+ T cells during differentiation. Thus, targeting these metabolic checkpoints may be a promising therapeutic strategy for both prevention of exhaustion and promotion of stemness of anti-tumor T cells.

VRK3 depletion induces cell cycle arrest and metabolic reprogramming of pontine diffuse midline glioma - H3K27 altered cells

We previously identified VRK3 as a specific vulnerability in DMG-H3K27M cells in a synthetic lethality screen targeting the whole kinome. The aim of the present study was to elucidate the mechanisms by which VRK3 depletion impact DMG-H3K27M cell fitness. Gene expression studies after VRK3 knockdown emphasized the inhibition of genes involved in G1/S transition of the cell cycle resulting in growth arrest in G1. Additionally, a massive modulation of genes involved in chromosome segregation was observed, concomitantly with a reduction in the level of phosphorylation of serine 10 and serine 28 of histone H3 supporting the regulation of chromatin condensation during cell division. This last effect could be partly due to a concomitant decrease of the chromatin kinase VRK1 in DMG following VRK3 knockdown. Furthermore, a metabolic switch specific to VRK3 function was observed towards increased oxidative phosphorylation without change in mitochondria content, that we hypothesized would represent a cell rescue mechanism. This study further explored the vulnerability of DMG-H3K27M cells to VRK3 depletion suggesting potential therapeutic combinations, e.g. with the mitochondrial ClpP protease activator ONC201.

Pan-tissue mitochondrial phenotyping reveals lower OXPHOS expression and function across cancer types

Targeting mitochondrial oxidative phosphorylation (OXPHOS) to treat cancer has been hampered due to serious side-effects potentially arising from the inability to discriminate between non-cancerous and cancerous mitochondria. Herein, comprehensive mitochondrial phenotyping was leveraged to define both the composition and function of OXPHOS across various murine cancers and compared to both matched normal tissues and other organs. When compared to both matched normal tissues, as well as high OXPHOS reliant organs like heart, intrinsic expression of the OXPHOS complexes, as well as OXPHOS flux were discovered to be consistently lower across distinct cancer types. Assuming intrinsic OXPHOS expression/function predicts OXPHOS reliance in vivo, these data suggest that pharmacologic blockade of mitochondrial OXPHOS likely compromises bioenergetic homeostasis in healthy oxidative organs prior to impacting tumor mitochondrial flux in a clinically meaningful way. Although these data caution against the use of indiscriminate mitochondrial inhibitors for cancer treatment, considerable heterogeneity was observed across cancer types with respect to both mitochondrial proteome composition and substrate-specific flux, highlighting the possibility for targeting discrete mitochondrial proteins or pathways unique to a given cancer type.

Mitochondrial heteroplasmic shifts reveal a positive selection of breast cancer

BACKGROUND

Breast cancer is, despite screening, not always detected early enough and is together with other tumor types known to shed genetic information in circulation. Unlike single-copy nuclear DNA, mitochondrial DNA (mtDNA) copies range from 100s to 10,000s per cell, thus providing a potentially alternative to identify potential missing cancer information in circulation at an early stage.

METHODS

To characterize mitochondrial mutation landscapes in breast cancer, whole mtDNA sequencing and bioinformatics analyses were performed on 86 breast cancer biopsies and 50 available matched baseline cancer-free whole blood samples from the same individuals, selected from a cohort of middle-aged women in Sweden. To determine whether the mutations can be detected in blood plasma prior to cancer diagnosis, we further designed a nested case-control study (n = 663) and validated the shortlisted mutations using droplet digital PCR.

RESULTS

We detected different mutation landscapes between biopsies and matched whole blood samples. Compared to whole blood samples, mtDNA from biopsies had higher heteroplasmic mutations in the D-loop region (P = 0.02), RNR2 (P = 0.005), COX1 (P = 0.037) and CYTB (P = 0.006). Furthermore, the germline mtDNA mutations had higher heteroplasmy level than the lost (P = 0.002) and de novo mutations (P = 0.04). The nonsynonymous to synonymous substitution ratio (dN/dS) was higher for the heteroplasmic mutations (P = 7.25 × 10-12) than that for the homoplasmic mutations, but the de novo (P = 0.06) and lost mutations (P = 0.03) had lower dN/dS than the germline mutations. Interestingly, we found that the critical regions for mitochondrial transcription: MT-HSP1 (odds ratio [OR]: 21.41), MT-TFH (OR: 7.70) and MT-TAS2 (OR: 3.62), had significantly higher heteroplasmic mutations than the rest of the D-loop sub-regions. Finally, we found that the presence of mt.16093T > C mutation increases 67% risk of developing breast cancer.

CONCLUSIONS

Our findings show that mitochondrial genetic landscape changes during cancer pathogenesis and positive selection of mtDNA heteroplasmic mutations in breast cancer. Most importantly, the mitochondrial mutations identified in biopsies can be traced back in matched plasma samples and could potentially be used as early breast cancer diagnostic biomarkers.

Forward genetic screens identify mechanisms of resistance to small molecule lactate dehydrogenase inhibitors

Altered metabolism is a hallmark of cancer; however, it has been difficult to specifically target metabolism in cancer for therapeutic benefit. Cancers with genetically defined defects in metabolic enzymes constitute a subset of cancers where targeting metabolism is potentially accessible. Hürthle cell carcinoma of the thyroid (HTC) tumors frequently harbor deleterious mitochondrial DNA (mtDNA) mutations in subunits of complex I of the mitochondrial electron transport chain (ETC). Previous work has shown that HTC models with deleterious mtDNA mutations exhibit mitochondrial ETC defects that expose lactate dehydrogenase (LDH) as a therapeutic vulnerability. Here, we performed forward genetic screens to identify mechanisms of resistance to small molecule LDH inhibitors. We identified two distinct mechanisms of resistance: upregulation of an LDH isoform and a compound-specific resistance mutation. Using these tools, we demonstrate that the anti-cancer activity of LDH inhibitors in cell line and xenograft models of complex I-mutant HTC is through on-target LDH inhibition.

Deleterious heteroplasmic mitochondrial mutations are associated with an increased risk of overall and cancer-specific mortality

Mitochondria carry their own circular genome and disruption of the mitochondrial genome is associated with various aging-related diseases. Unlike the nuclear genome, mitochondrial DNA (mtDNA) can be present at 1000 s to 10,000 s copies in somatic cells and variants may exist in a state of heteroplasmy, where only a fraction of the DNA molecules harbors a particular variant. We quantify mtDNA heteroplasmy in 194,871 participants in the UK Biobank and find that heteroplasmy is associated with a 1.5-fold increased risk of all-cause mortality. Additionally, we functionally characterize mtDNA single nucleotide variants (SNVs) using a constraint-based score, mitochondrial local constraint score sum (MSS) and find it associated with all-cause mortality, and with the prevalence and incidence of cancer and cancer-related mortality, particularly leukemia. These results indicate that mitochondria may have a functional role in certain cancers, and mitochondrial heteroplasmic SNVs may serve as a prognostic marker for cancer, especially for leukemia.

Mitochondrial metabolism-related signature depicts immunophenotype and predicts therapeutic response in testicular germ cell tumors

In recent years, there has been growing evidence linking mitochondrial dysfunction to the development and progression of cancer. However, the role of mitochondrial metabolism-related genes (MMRGs) in testicular germ cell tumor (TGCT) remains unclear. We downloaded clinical pathology, transcriptome, and somatic mutation data for TGCT from public databases and conducted univariate Cox regression analysis to investigate prognostic correlations. We also used consensus clustering to identify molecular subtypes, comparing differential expression genes, biological processes, Kyoto Encyclopedia of Genes and Genomes pathways, mutations, prognosis, immune infiltration, drug sensitivity, and immune therapeutic response between these subtypes. We constructed multi-gene risk features and nomograms for TGCT prognosis. Fifteen MMRGs were significantly correlated with progression-free survival in TGCT patients. Based on these genes, we identified 2 molecular subtypes which showed significant differences in somatic mutations, prognosis, and immune cell infiltration. These subtypes could also indicate drug sensitivity and immune therapeutic response; the subtype with poor prognosis showed a higher potential benefit from some drugs and immunotherapy. Abnormalities in immune-related biological processes and extracellular matrix as well as Kyoto Encyclopedia of Genes and Genomes pathways such as PI3K-AKT signaling pathway, pat5hways in cancer, primary immunodeficiency, and neutrophil extracellular trap formation were associated with significant differences in phenotypes among subtypes. Finally, we constructed an 8-gene TGCT risk feature based on differential expression genes between subtypes which performed well in TGCT patient prognostic evaluation. Our study elucidated the prognostic correlation between MMRGs and TGCT and established MMRG-derived molecular subtypes and risk features for personalized treatment of TGCT which have potential clinical application value.

A systematic comparison of human mitochondrial genome assembly tools

BACKGROUND

Mitochondria are the cell organelles that produce most of the chemical energy required to power the cell's biochemical reactions. Despite being a part of a eukaryotic host cell, the mitochondria contain a separate genome whose origin is linked with the endosymbiosis of a prokaryotic cell by the host cell and encode independent genomic information throughout their genomes. Mitochondrial genomes accommodate essential genes and are regularly utilized in biotechnology and phylogenetics. Various assemblers capable of generating complete mitochondrial genomes are being continuously developed. These tools often use whole-genome sequencing data as an input containing reads from the mitochondrial genome. Till now, no published work has explored the systematic comparison of all the available tools for assembling human mitochondrial genomes using short-read sequencing data. This evaluation is required to identify the best tool that can be well-optimized for small-scale projects or even national-level research.

RESULTS

In this study, we have tested the mitochondrial genome assemblers for both simulated datasets and whole genome sequencing (WGS) datasets of humans. For the highest computational setting of 16 computational threads with the simulated dataset having 1000X read depth, MitoFlex took the least execution time of 69 s, and IOGA took the longest execution time of 1278 s. NOVOPlasty utilized the least computational memory of approximately 0.098 GB for the same setting, whereas IOGA utilized the highest computational memory of 11.858 GB. In the case of WGS datasets for humans, GetOrganelle and MitoFlex performed the best in capturing the SNPs information with a mean F1-score of 0.919 at the sequencing depth of 10X. MToolBox and NOVOPlasty performed consistently across all sequencing depths with a mean F1 score of 0.897 and 0.890, respectively.

CONCLUSIONS

Based on the overall performance metrics and consistency in assembly quality for all sequencing data, MToolBox performed the best. However, NOVOPlasty was the second fastest tool in execution time despite being single-threaded, and it utilized the least computational resources among all the assemblers when tested on simulated datasets. Therefore, NOVOPlasty may be more practical when there is a significant sample size and a lack of computational resources. Besides, as long-read sequencing gains popularity, mitochondrial genome assemblers must be developed to use long-read sequencing data.

Exosomal transfer leads to chemoresistance through oxidative phosphorylation-mediated stemness phenotype in colorectal cancer

Background: Recently years have seen the increasing evidence identifying that OXPHOS is involved in different processes of tumor progression and metastasis and has been proposed to be a potential therapeutical target for cancer treatment. However, the exploration in oxidative phosphorylation-mediated chemoresistance is still scarce. In our study, we identify exosomal transfer leads to chemoresistance by reprogramming metabolic phenotype in recipient cells. Methods: RNA sequencing analysis was used to screen altered targets mediating exosome transfer-induced chemoresistance. Seahorse assay allowed us to measure mitochondrial respiration. Stemness was measured by spheroids formation assay. Serum exosomes were isolated for circ_0001610 quantification. Results: The induced oxidative phosphorylation leads to more stem-like properties, which is dependent on the transfer of exosomal circ_0001610. Exosome transfer results in the removal of miR-30e-5p-mediated suppression of PGC-1a, a master of mitochondrial biogenesis and function. Consequently, increased PGC-1a reshapes cellular metabolism towards oxidative phosphorylation, leading to chemoresistance. Inhibition of OXPHOS or exosomal si-circ_0001610 increases the sensitivity of chemotherapy by decreasing cell stemness in vitro and in vivo. Conclusion: Our data suggests that exosomal circ_0001610-induced OXPHOS plays an important role in chemoresistance and supports a therapeutical potential of circ_0001610 inhibitors in the treatment of oxaliplatin-resistant colorectal cancer by manipulating cell stemness.

Pathogenic mitochondrial DNA mutations inhibit melanoma metastasis

Mitochondrial DNA (mtDNA) mutations are frequently observed in cancer, but their contribution to tumor progression is controversial. To evaluate the impact of mtDNA variants on tumor growth and metastasis, we created human melanoma cytoplasmic hybrid (cybrid) cell lines transplanted with wildtype mtDNA or pathogenic mtDNA encoding variants that partially or completely inhibit oxidative phosphorylation. Homoplasmic pathogenic mtDNA cybrids reliably established tumors despite dysfunctional oxidative phosphorylation. However, pathogenic mtDNA variants disrupted spontaneous metastasis of subcutaneous tumors and decreased the abundance of circulating melanoma cells in the blood. Pathogenic mtDNA did not induce anoikis or inhibit organ colonization of melanoma cells following intravenous injections. Instead, migration and invasion were reduced, indicating that limited circulation entry functions as a metastatic bottleneck amidst mtDNA dysfunction. Furthermore, analysis of selective pressure exerted on the mitochondrial genomes of heteroplasmic cybrid lines revealed a suppression of pathogenic mtDNA allelic frequency during melanoma growth. Collectively, these findings demonstrate that functional mtDNA is favored during melanoma growth and enables metastatic entry into the blood.

High expression of novel biomarker TBRG4 promotes the progression and invasion of oral squamous cell carcinoma

BACKGROUND

Transforming growth factor β regulator 4 (TBRG4) is a potential prognostic indicator in various cancers, especially squamous cell carcinomas, and is associated with disease amelioration and poor outcomes. The study aimed to assess the expression pattern of TBRG4 in patients with operable oral squamous cell carcinoma (OSCC) to understand its role in tumour progression using indicators of disease outcome like tumour stage, grade, nodal metastasis, and pattern of invasion.

METHODS

TBRG4 expression was assessed by analyzing 51 cancer and adjacent non-cancerous tissues of OSCC patients using quantitative real-time PCR, and Western blot. TBRG4 expression was also analysed in The Cancer Genome Atlas Head-Neck Squamous Cell Carcinoma (TCGA-HNSC) dataset using the UALCAN tool (http://ualcan.path.uab.edu/). The relationship between TBRG4 expression and the patient's prognosis was analysed using Kaplan-Meier plotter.

RESULTS

Both mRNA and protein levels of TBRG4 were significantly increased in OSCC tissues. The TBGR4 expression was significantly associated with advanced stages (III and IV) and the worst pattern of invasion (WPOI-4 and 5). High TBRG4 expression was also significantly associated with reduced overall survival (p = 0.011). In addition, the analysis of TBRG4 gene expression and clinical data from TCGA, identified that TBRG4 was highly expressed in HPV negative OSCC patients and positively correlated with worst overall survival.

CONCLUSION

The present study suggests that the high expression of TBRG4 might serve as a novel prognostic biomarker for HPV-negative OSCC, which can be validated by future additional investigations in larger cohorts along with functional studies.

Mitochondrial heterogeneity in diseases

As key organelles involved in cellular metabolism, mitochondria frequently undergo adaptive changes in morphology, components and functions in response to various environmental stresses and cellular demands. Previous studies of mitochondria research have gradually evolved, from focusing on morphological change analysis to systematic multiomics, thereby revealing the mitochondrial variation between cells or within the mitochondrial population within a single cell. The phenomenon of mitochondrial variation features is defined as mitochondrial heterogeneity. Moreover, mitochondrial heterogeneity has been reported to influence a variety of physiological processes, including tissue homeostasis, tissue repair, immunoregulation, and tumor progression. Here, we comprehensively review the mitochondrial heterogeneity in different tissues under pathological states, involving variant features of mitochondrial DNA, RNA, protein and lipid components. Then, the mechanisms that contribute to mitochondrial heterogeneity are also summarized, such as the mutation of the mitochondrial genome and the import of mitochondrial proteins that result in the heterogeneity of mitochondrial DNA and protein components. Additionally, multiple perspectives are investigated to better comprehend the mysteries of mitochondrial heterogeneity between cells. Finally, we summarize the prospective mitochondrial heterogeneity-targeting therapies in terms of alleviating mitochondrial oxidative damage, reducing mitochondrial carbon stress and enhancing mitochondrial biogenesis to relieve various pathological conditions. The possibility of recent technological advances in targeted mitochondrial gene editing is also discussed.

Chromophobe renal cell carcinoma

Chromophobe renal cell carcinoma (ChRCC) is the second most common variant histology (non-clear cell) RCC. ChRCC is distinct from clear cell RCC (ccRCC) in terms of genetics, genomics, metabolism, cell of origin, and response to targeted and immune therapies. The pathogenesis of ChRCC remains unclear, but current data suggest two potential mechanisms: mTORC1 hyperactivation through PTEN pathway mutations and mitochondrial dysfunction leading to oxidative stress. There are no specific approved treatments for ChRCC, although some responses to tyrosine kinase and mTOR inhibitors have been observed. Response to immunotherapy is generally limited. Targetable pathways involving innate lymphoid cells/IL-15 and cysteine homeostasis/ferroptosis have recently been identified.