Mitochondrial respiratory gene expression is suppressed in many cancers

A computational study of gene expression patterns in head and neck squamous cell carcinoma using TCGA data

Aim: Head and Neck squamous cell carcinoma (HNSCC) is the second most prevalent cancer in Pakistan. Methods: Gene expression data from TCGA and GETx for normal genes to analyze Differentially Expressed Genes (DEGs). Data was further investigated using the Enrichr tool to perform Gene Ontology (GO). Results: Our analysis identified most significantly differentially expressed genes and explored their established cellular functions as well as their potential involvement in tumor development. We found that the highly expressed Keratin family and S100A9 genes. The under-expressed genes KRT4 and KRT13 provide instructions for the production of keratin proteins. Conclusion: Our study suggests that factors such as poor oral hygiene and smokeless tobacco can result in oral stress and cellular damage and cause cancer.

Mitochondrial complex I promotes kidney cancer metastasis

Most kidney cancers are metabolically dysfunctional1-4, but how this dysfunction affects cancer progression in humans is unknown. We infused 13C-labelled nutrients in over 80 patients with kidney cancer during surgical tumour resection. Labelling from [U-13C]glucose varies across subtypes, indicating that the kidney environment alone cannot account for all tumour metabolic reprogramming. Compared with the adjacent kidney, clear cell renal cell carcinomas (ccRCCs) display suppressed labelling of tricarboxylic acid (TCA) cycle intermediates in vivo and in ex vivo organotypic cultures, indicating that suppressed labelling is tissue intrinsic. [1,2-13C]acetate and [U-13C]glutamine infusions in patients, coupled with measurements of respiration in isolated human kidney and tumour mitochondria, reveal lower electron transport chain activity in ccRCCs that contributes to decreased oxidative and enhanced reductive TCA cycle labelling. However, ccRCC metastases unexpectedly have enhanced TCA cycle labelling compared with that of primary ccRCCs, indicating a divergent metabolic program during metastasis in patients. In mice, stimulating respiration or NADH recycling in kidney cancer cells is sufficient to promote metastasis, whereas inhibiting electron transport chain complex I decreases metastasis. These findings in humans and mice indicate that metabolic properties and liabilities evolve during kidney cancer progression, and that mitochondrial function is limiting for metastasis but not growth at the original site.

Lack of Mitochondrial DNA Provides Metabolic Advantage in Yeast Osmoadaptation

Alterations in mitochondrial function have been linked to a variety of cellular and organismal stress responses including apoptosis, aging, neurodegeneration and tumorigenesis. However, adaptation to mitochondrial dysfunction can occur through the activation of survival pathways, whose mechanisms are still poorly understood. The yeast Saccharomyces cerevisiae is an invaluable model organism for studying how mitochondrial dysfunction can affect stress response and adaptation processes. In this study, we analyzed and compared in the absence and in the presence of osmostress wild-type cells with two models of cells lacking mitochondrial DNA: ethidium bromide-treated cells (ρ0) and cells lacking the mitochondrial pyrimidine nucleotide transporter RIM2 (ΔRIM2). Our results revealed that the lack of mitochondrial DNA provides an advantage in the kinetics of stress response. Additionally, wild-type cells exhibited higher osmosensitivity in the presence of respiratory metabolism. Mitochondrial mutants showed increased glycerol levels, required in the short-term response of yeast osmoadaptation, and prolonged oxidative stress. The involvement of the mitochondrial retrograde signaling in osmoadaptation has been previously demonstrated. The expression of CIT2, encoding the peroxisomal isoform of citrate synthase and whose up-regulation is prototypical of RTG pathway activation, appeared to be increased in the mutants. Interestingly, selected TCA cycle genes, CIT1 and ACO1, whose expression depends on RTG signaling upon stress, showed a different regulation in ρ0 and ΔRIM2 cells. These data suggest that osmoadaptation can occur through different mechanisms in the presence of mitochondrial defects and will allow us to gain insight into the relationships among metabolism, mitochondria-mediated stress response, and cell adaptation.

FSHD muscle shows perturbation in fibroadipogenic progenitor cells, mitochondrial function and alternative splicing independently of inflammation

Facioscapulohumeral muscular dystrophy (FSHD) is a prevalent, incurable myopathy. FSHD is highly heterogeneous, with patients following a variety of clinical trajectories, complicating clinical trials. Skeletal muscle in FSHD undergoes fibrosis and fatty replacement that can be accelerated by inflammation, adding to heterogeneity. Well controlled molecular studies are thus essential to both categorize FSHD patients into distinct subtypes and understand pathomechanisms. Here, we further analyzed RNA-sequencing data from 24 FSHD patients, each of whom donated a biopsy from both a non-inflamed (TIRM-) and inflamed (TIRM+) muscle, and 15 FSHD patients who donated peripheral blood mononucleated cells (PBMCs), alongside non-affected control individuals. Differential gene expression analysis identified suppression of mitochondrial biogenesis and up-regulation of fibroadipogenic progenitor (FAP) gene expression in FSHD muscle, which was particularly marked on inflamed samples. PBMCs demonstrated suppression of antigen presentation in FSHD. Gene expression deconvolution revealed FAP expansion as a consistent feature of FSHD muscle, via meta-analysis of 7 independent transcriptomic datasets. Clustering of muscle biopsies separated patients in an unbiased manner into clinically mild and severe subtypes, independently of known disease modifiers (age, sex, D4Z4 repeat length). Lastly, the first genome-wide analysis of alternative splicing in FSHD muscle revealed perturbation of autophagy, BMP2 and HMGB1 signalling. Overall, our findings reveal molecular subtypes of FSHD with clinical relevance and identify novel pathomechanisms for this highly heterogeneous condition.

Comparison of Transcriptomic Analysis of the Conjunctiva in Glaucoma-Treated Eyes with Dry Eyes and Healthy Controls

Ocular surface disease (OSD) associated with topical glaucoma drugs is a common issue impacting treatment adherence. We aimed to identify conjunctival transcriptomic changes in glaucoma and dry eye patients, comparing them to healthy controls. Bulbar conjunctival specimens were collected via impression cytology from 33 patients treated for glaucoma, 9 patients with dry eye, and 14 healthy controls. RNA extraction and bulk RNA sequencing were performed, followed by bioinformatics analysis to detect gene dysregulation. Ingenuity pathways analysis (IPA) identified pathways and biological processes associated with these transcriptomic changes. Sequencing analysis revealed 200 modified genes in glaucoma patients compared to healthy individuals, 233 differentially expressed genes in dry eye patients versus controls, and 650 genes in treated versus dry eye samples. In glaucoma patients, 79% of altered pathways were related to host defense, while dry eye patients showed a 39% involvement of host response, 15% in cellular proliferation and integrity, and 16% of mitochondrial dysfunction. These findings were validated through qRT-PCR. Glaucoma patients showed an intensified conjunctival immune response as a potential cause of OSD, whereas in dry eye patients, in addition to the immune response, other mechanisms such as mitochondrial dysfunction or reduced cellular proliferation were observed.

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.

Characterization of tumour microenvironment reprogramming reveals invasion in epithelial ovarian carcinoma

BACKGROUND

Patients with epithelial ovarian carcinoma (EOC) are usually diagnosed at an advanced stage with tumour cell invasion. However, identifying the underlying molecular mechanisms and biomarkers of EOC proliferation and invasion remains challenging.

RESULTS

Herein, we explored the relationship between tumour microenvironment (TME) reprogramming and tissue invasion based on single-cell RNA sequencing (scRNA-seq) datasets. Interestingly, hypoxia, oxidative phosphorylation (OXPHOS) and glycolysis, which have biologically active trajectories during epithelial mesenchymal transition (EMT), were positively correlated. Moreover, energy metabolism and anti-apoptotic activity were found to be critical contributors to intratumor heterogeneity. In addition, HMGA1, EGR1 and RUNX1 were found to be critical drivers of the EMT process in EOC. Experimental validation revealed that suppressing EGR1 expression inhibited tumour cell invasion, significantly upregulated the expression of E-cadherin and decreased the expression of N-cadherin. In cell components analysis, cancer-associated fibroblasts (CAFs) were found to significantly contribute to immune infiltration and tumour invasion, and the accumulation of CAFs was associated with poorer patient survival.

CONCLUSION

We revealed the molecular mechanism and biomarkers of tumour invasion and TME reprogramming in EOC, which provides effective targets for the suppression of tumour invasion.

Integrated transcriptomic and metabolomic analysis reveals the metabolic programming of GM-CSF- and M-CSF- differentiated mouse macrophages

Macrophages play a critical role in the inflammatory response and tumor development. Macrophages are primarily divided into pro-inflammatory M1-like and anti-inflammatory M2-like macrophages based on their activation status and functions. In vitro macrophage models could be derived from mouse bone marrow cells stimulated with two types of differentiation factors: GM-CSF (GM-BMDMs) and M-CSF (M-BMDMs), to represent M1- and M2-like macrophages, respectively. Since macrophage differentiation requires coordinated metabolic reprogramming and transcriptional rewiring in order to fulfill their distinct roles, we combined both transcriptome and metabolome analysis, coupled with experimental validation, to gain insight into the metabolic status of GM- and M-BMDMs. The data revealed higher levels of the tricarboxylic acid cycle (TCA cycle), oxidative phosphorylation (OXPHOS), fatty acid oxidation (FAO), and urea and ornithine production from arginine in GM-BMDMs, and a preference for glycolysis, fatty acid storage, bile acid metabolism, and citrulline and nitric oxide (NO) production from arginine in M-BMDMs. Correlation analysis with the proteomic data showed high consistency in the mRNA and protein levels of metabolic genes. Similar results were also obtained when compared to RNA-seq data of human monocyte derived macrophages from the GEO database. Furthermore, canonical macrophage functions such as inflammatory response and phagocytosis were tightly associated with the representative metabolic pathways. In the current study, we identified the core metabolites, metabolic genes, and functional terms of the two distinct mouse macrophage populations. We also distinguished the metabolic influences of the differentiation factors GM-CSF and M-CSF, and wish to provide valuable information for in vitro macrophage studies.

Inter-organellar nucleic acid communication by a mitochondrial tRNA regulates nuclear metabolic transcription

Efficient communication between mitochondria and the nucleus underlies homoeostatic metabolic control, though the involved mitochondrial factors and their mechanisms are poorly defined. Here, we report the surprising detection of multiple mitochondrial-derived transfer RNAs (mito-tRNAs) within the nuclei of human cells. Focused studies of nuclear-transported mito-tRNA-asparagine (mtAsn) revealed that its cognate charging enzyme (NARS2) is also present in the nucleus. MtAsn promoted interaction of NARS2 with histone deacetylase 2 (HDAC2), and repressed HDAC2 association with specific chromatin loci. Perturbation of this axis using antisense oligonucleotides promoted nucleotide biogenesis and enhanced breast cancer growth, and RNA and nascent transcript sequencing demonstrated specific alterations in the transcription of nuclear genes. These findings uncover nucleic-acid mediated communication between two organelles and the existence of a machinery for nuclear gene regulation by a mito-tRNA that restricts tumor growth through metabolic control.

Highlights

Multiple mitochondrial-derived tRNAs are detected in human cell nucleiMtAsn promotes binding between NARS2 and HDAC2Metabolic alterations driven by mtAsn impact cell proliferationMtAsn inhibition releases HDAC2 to bind and transcriptionally regulate multiple nuclear genes.

REFINED-CNN framework for survival prediction with high-dimensional features

Robust and accurate survival prediction of clinical trials using high-throughput genomics data is a fundamental challenge in pharmacogenomics. Current machine learning tools often provide limited predictive performance and model interpretation in these settings. In the present study, we extend the application of REFINED-CNN from regression tasks to making survival predictions, by mapping high-dimensional RNA sequencing data into REFINED images which are conducive to CNN modeling. We show that the REFINED-CNN survival model can be easily adapted to new tasks of a similar nature (e.g., predicting on new cancer types) using transfer learning with a low number of patients. Furthermore, the model can also be interpreted both locally and globally through risk score back propagation that quantifies each feature (e.g., gene) importance in survival prediction task for the patient or cancer type of interest.

Impact of supraphysiologic MDM2 expression on chromatin networks and therapeutic responses in sarcoma

Amplification of MDM2 on supernumerary chromosomes is a common mechanism of P53 inactivation across tumors. Here, we investigated the impact of MDM2 overexpression on chromatin, gene expression, and cellular phenotypes in liposarcoma. Three independent regulatory circuits predominate in aggressive, dedifferentiated tumors. RUNX and AP-1 family transcription factors bind mesenchymal gene enhancers. P53 and MDM2 co-occupy enhancers and promoters associated with P53 signaling. When highly expressed, MDM2 also binds thousands of P53-independent growth and stress response genes, whose promoters engage in multi-way topological interactions. Overexpressed MDM2 concentrates within nuclear foci that co-localize with PML and YY1 and could also contribute to P53-independent phenotypes associated with supraphysiologic MDM2. Importantly, we observe striking cell-to-cell variability in MDM2 copy number and expression in tumors and models. Whereas liposarcoma cells are generally sensitive to MDM2 inhibitors and their combination with pro-apoptotic drugs, MDM2-high cells tolerate them and may underlie the poor clinical efficacy of these agents.

Chronic airway epithelial hypoxia exacerbates injury in muco-obstructive lung disease through mucus hyperconcentration

Unlike solid organs, human airway epithelia derive their oxygen from inspired air rather than the vasculature. Many pulmonary diseases are associated with intraluminal airway obstruction caused by aspirated foreign bodies, virus infection, tumors, or mucus plugs intrinsic to airway disease, including cystic fibrosis (CF). Consistent with requirements for luminal O2, airway epithelia surrounding mucus plugs in chronic obstructive pulmonary disease (COPD) lungs are hypoxic. Despite these observations, the effects of chronic hypoxia (CH) on airway epithelial host defense functions relevant to pulmonary disease have not been investigated. Molecular characterization of resected human lungs from individuals with a spectrum of muco-obstructive lung diseases (MOLDs) or COVID-19 identified molecular features of chronic hypoxia, including increased EGLN3 expression, in epithelia lining mucus-obstructed airways. In vitro experiments using cultured chronically hypoxic airway epithelia revealed conversion to a glycolytic metabolic state with maintenance of cellular architecture. Chronically hypoxic airway epithelia unexpectedly exhibited increased MUC5B mucin production and increased transepithelial Na+ and fluid absorption mediated by HIF1α/HIF2α-dependent up-regulation of β and γENaC (epithelial Na+ channel) subunit expression. The combination of increased Na+ absorption and MUC5B production generated hyperconcentrated mucus predicted to perpetuate obstruction. Single-cell and bulk RNA sequencing analyses of chronically hypoxic cultured airway epithelia revealed transcriptional changes involved in airway wall remodeling, destruction, and angiogenesis. These results were confirmed by RNA-in situ hybridization studies of lungs from individuals with MOLD. Our data suggest that chronic airway epithelial hypoxia may be central to the pathogenesis of persistent mucus accumulation in MOLDs and associated airway wall damage.

The Interplay between Dysregulated Metabolism and Epigenetics in Cancer

Cellular metabolism (or energetics) and epigenetics are tightly coupled cellular processes. It is arguable that of all the described cancer hallmarks, dysregulated cellular energetics and epigenetics are the most tightly coregulated. Cellular metabolic states regulate and drive epigenetic changes while also being capable of influencing, if not driving, epigenetic reprogramming. Conversely, epigenetic changes can drive altered and compensatory metabolic states. Cancer cells meticulously modify and control each of these two linked cellular processes in order to maintain their tumorigenic potential and capacity. This review aims to explore the interplay between these two processes and discuss how each affects the other, driving and enhancing tumorigenic states in certain contexts.

Mitochondria-derived small RNAs as diagnostic biomarkers in lung cancer patients through a novel ratio-based expression analysis methodology

Small non-coding RNAs are potential diagnostic biomarkers for lung cancer. Mitochondria-derived small RNA (mtRNA) is a novel regulatory small non-coding RNA that only recently has been identified and cataloged. Currently, there are no reports of studies of mtRNA in human lung cancer. Currently, normalization methods are unstable, and they often fail to identify differentially expressed small non-coding RNAs (sncRNAs). In order to identify reliable biomarkers for lung cancer screening, we used a ratio-based method using mtRNAs newly discovered in human peripheral blood mononuclear cells. In the discovery cohort (AUC = 0.981) and independent validation cohort (AUC = 0.916) the prediction model of eight mtRNA ratios distinguished lung cancer patients from controls. The prediction model will provide reliable biomarkers that will allow blood-based screening to become more feasible and will help make lung cancer diagnosis more accurate in clinical practice.

Inhibition of mitochondrial translocase SLC25A5 and histone deacetylation is an effective combination therapy in neuroblastoma

The mitochondrion is a gatekeeper of apoptotic processes, and mediates drug resistance to several chemotherapy agents used to treat cancer. Neuroblastoma is a common solid cancer in young children with poor clinical outcomes following conventional chemotherapy. We sought druggable mitochondrial protein targets in neuroblastoma cells. Among mitochondria-associated gene targets, we found that high expression of the mitochondrial adenine nucleotide translocase 2 (SLC25A5/ANT2), was a strong predictor of poor neuroblastoma patient prognosis and contributed to a more malignant phenotype in pre-clinical models. Inhibiting this transporter with PENAO reduced cell viability in a panel of neuroblastoma cell lines in a TP53-status-dependant manner. We identified the histone deacetylase inhibitor, suberanilohydroxamic acid (SAHA), as the most effective drug in clinical use against mutant TP53 neuroblastoma cells. SAHA and PENAO synergistically reduced cell viability, and induced apoptosis, in neuroblastoma cells independent of TP53-status. The SAHA and PENAO drug combination significantly delayed tumour progression in pre-clinical neuroblastoma mouse models, suggesting that these clinically advanced inhibitors may be effective in treating the disease.

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

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

Mitochondrial redox adaptations enable alternative aspartate synthesis in SDH-deficient cells

The oxidative tricarboxylic acid (TCA) cycle is a central mitochondrial pathway integrating catabolic conversions of NAD +to NADH and anabolic production of aspartate, a key amino acid for cell proliferation. Several TCA cycle components are implicated in tumorigenesis, including loss-of-function mutations in subunits of succinate dehydrogenase (SDH), also known as complex II of the electron transport chain (ETC), but mechanistic understanding of how proliferating cells tolerate the metabolic defects of SDH loss is still lacking. Here, we identify that SDH supports human cell proliferation through aspartate synthesis but, unlike other ETC impairments, the effects of SDH inhibition are not ameliorated by electron acceptor supplementation. Interestingly, we find aspartate production and cell proliferation are restored to SDH-impaired cells by concomitant inhibition of ETC complex I (CI). We determine that the benefits of CI inhibition in this context depend on decreasing mitochondrial NAD+/NADH, which drives SDH-independent aspartate production through pyruvate carboxylation and reductive carboxylation of glutamine. We also find that genetic loss or restoration of SDH selects for cells with concordant CI activity, establishing distinct modalities of mitochondrial metabolism for maintaining aspartate synthesis. These data therefore identify a metabolically beneficial mechanism for CI loss in proliferating cells and reveal how compartmentalized redox changes can impact cellular fitness.

Mitochondrial metabolism in primary and metastatic human kidney cancers

Most kidney cancers display evidence of metabolic dysfunction1-4 but how this relates to cancer progression in humans is unknown. We used a multidisciplinary approach to infuse 13C-labeled nutrients during surgical tumour resection in over 70 patients with kidney cancer. Labeling from [U-13C]glucose varies across cancer subtypes, indicating that the kidney environment alone cannot account for all metabolic reprogramming in these tumours. Compared to the adjacent kidney, clear cell renal cell carcinomas (ccRCC) display suppressed labelling of tricarboxylic acid (TCA) cycle intermediates in vivo and in organotypic slices cultured ex vivo, indicating that suppressed labeling is tissue intrinsic. Infusions of [1,2-13C]acetate and [U-13C]glutamine in patients, coupled with respiratory flux of mitochondria isolated from kidney and tumour tissue, reveal primary defects in mitochondrial function in human ccRCC. However, ccRCC metastases unexpectedly have enhanced labeling of TCA cycle intermediates compared to primary ccRCCs, indicating a divergent metabolic program during ccRCC metastasis in patients. In mice, stimulating respiration in ccRCC cells is sufficient to promote metastatic colonization. Altogether, these findings indicate that metabolic properties evolve during human kidney cancer progression, and suggest that mitochondrial respiration may be limiting for ccRCC metastasis but not for ccRCC growth at the site of origin.

Neoantigens: promising targets for cancer therapy

Recent advances in neoantigen research have accelerated the development and regulatory approval of tumor immunotherapies, including cancer vaccines, adoptive cell therapy and antibody-based therapies, especially for solid tumors. Neoantigens are newly formed antigens generated by tumor cells as a result of various tumor-specific alterations, such as genomic mutation, dysregulated RNA splicing, disordered post-translational modification, and integrated viral open reading frames. Neoantigens are recognized as non-self and trigger an immune response that is not subject to central and peripheral tolerance. The quick identification and prediction of tumor-specific neoantigens have been made possible by the advanced development of next-generation sequencing and bioinformatic technologies. Compared to tumor-associated antigens, the highly immunogenic and tumor-specific neoantigens provide emerging targets for personalized cancer immunotherapies, and serve as prospective predictors for tumor survival prognosis and immune checkpoint blockade responses. The development of cancer therapies will be aided by understanding the mechanism underlying neoantigen-induced anti-tumor immune response and by streamlining the process of neoantigen-based immunotherapies. This review provides an overview on the identification and characterization of neoantigens and outlines the clinical applications of prospective immunotherapeutic strategies based on neoantigens. We also explore their current status, inherent challenges, and clinical translation potential.

Mitochondrial Dynamics during Development

Mitochondria are dynamic membrane-bound organelles in eukaryotic cells. These are important for the generation of chemical energy needed to power various cellular functions and also support metabolic, energetic, and epigenetic regulation in various cells. These organelles are also important for communication with the nucleus and other cellular structures, to maintain developmental sequences and somatic homeostasis, and for cellular adaptation to stress. Increasing information shows mitochondrial defects as an important cause of inherited disorders in different organ systems. In this article, we provide an extensive review of ontogeny, ultrastructural morphology, biogenesis, functional dynamics, important clinical manifestations of mitochondrial dysfunction, and possibilities for clinical intervention. We present information from our own clinical and laboratory research in conjunction with information collected from an extensive search in the databases PubMed, EMBASE, and Scopus.

Nuclear and mitochondrial DNA alterations in pheochromocytomas and paragangliomas, and their potential treatment

Mitochondrial DNA (mtDNA) alterations have been reported in different types of cancers and are suggested to play important roles in cancer development and metastasis. However, there is little information about its involvement in pheochromocytomas and paragangliomas (PCCs/PGLs) formation. PCCs and PGLs are rare endocrine tumors of the chromaffin cells in the adrenal medulla and extra-adrenal paraganglia that can synthesize and secrete catecholamines. Over the last 3 decades, the genetic background of about 60% of PCCs/PGLs involving nuclear DNA alterations has been determined. Recently, a study showed that mitochondrial alterations can be found in around 17% of the remaining PCCs/PGLs. In this review, we summarize recent knowledge regarding both nuclear and mitochondrial alterations and their involvement in PCCs/PGLs. We also provide brief insights into the genetics and the molecular pathways associated with PCCs/PGLs and potential therapeutical targets.

Roles of mitochondrial genetics in cancer metastasis

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

Mitochondrial transfer/transplantation: an emerging therapeutic approach for multiple diseases

Mitochondria play a pivotal role in energy generation and cellular physiological processes. These organelles are highly dynamic, constantly changing their morphology, cellular location, and distribution in response to cellular stress. In recent years, the phenomenon of mitochondrial transfer has attracted significant attention and interest from biologists and medical investigators. Intercellular mitochondrial transfer occurs in different ways, including tunnelling nanotubes (TNTs), extracellular vesicles (EVs), and gap junction channels (GJCs). According to research on intercellular mitochondrial transfer in physiological and pathological environments, mitochondrial transfer hold great potential for maintaining body homeostasis and regulating pathological processes. Multiple research groups have developed artificial mitochondrial transfer/transplantation (AMT/T) methods that transfer healthy mitochondria into damaged cells and recover cellular function. This paper reviews intercellular spontaneous mitochondrial transfer modes, mechanisms, and the latest methods of AMT/T. Furthermore, potential application value and mechanism of AMT/T in disease treatment are also discussed.

Respiratory Complex I dysfunction in cancer: from a maze of cellular adaptive responses to potential therapeutic strategies

Mitochondria act as key organelles in cellular bioenergetics and biosynthetic processes producing signals that regulate different molecular networks for proliferation and cell death. This ability is also preserved in pathologic contexts such as tumorigenesis, during which bioenergetic changes and metabolic reprogramming confer flexibility favoring cancer cell survival in a hostile microenvironment. Although different studies epitomize mitochondrial dysfunction as a protumorigenic hit, genetic ablation or pharmacological inhibition of respiratory complex I causing a severe impairment is associated with a low-proliferative phenotype. In this scenario, it must be considered that despite the initial delay in growth, cancer cells may become able to resume proliferation exploiting molecular mechanisms to overcome growth arrest. Here, we highlight the current knowledge on molecular responses activated by complex I-defective cancer cells to bypass physiological control systems and to re-adapt their fitness during microenvironment changes. Such adaptive mechanisms could reveal possible novel molecular players in synthetic lethality with complex I impairment, thus providing new synergistic strategies for mitochondrial-based anticancer therapy.

Establishment of a platform for measuring mitochondrial oxygen consumption rate for cardiac mitochondrial toxicity

The heart has an abundance of mitochondria since cardiac muscles require copious amounts of energy for providing continuous blood through the circulatory system, thereby implying that myocardial function is largely reliant on mitochondrial energy. Thus, cardiomyocytes are susceptible to mitochondrial dysfunction and are likely targets of mitochondrial toxic drugs. Various methods have been developed to evaluate mitochondrial toxicity by evaluating toxicological mechanisms, but an optimized and standardized assay for cardiomyocytes remains unmet. We have therefore attempted to standardize the evaluation system for determining cardiac mitochondrial toxicity, using AC16 human and H9C2 rat cardiomyocytes. Three clinically administered drugs (acetaminophen, amiodarone, and valproic acid) and two anticancer drugs (doxorubicin and tamoxifen) which are reported to have mitochondrial effects, were applied in this study. The oxygen consumption rate (OCR), which directly reflects mitochondrial function, and changes in mRNA levels of mitochondrial respiratory complex I to complex V, were analyzed. Our results reveal that exposure to all five drugs results in a concentration-dependent decrease in the basal and maximal levels of OCR in AC16 cells and H9C2 cells. In particular, marked reduction in the OCR was observed after treatment with doxorubicin. The reduction in OCR after exposure to mitochondrial toxic drugs was found to be associated with reduced mRNA expression in the mitochondrial respiratory complexes, suggesting that the cardiac mitochondrial toxicity of drugs is majorly due to dysfunction of mitochondrial respiration. Based on the results of this study, we established and standardized a protocol to measure OCR in cardiomyocytes. We expect that this standardized evaluation system for mitochondrial toxicity can be applied as basic data for establishing a screening platform to evaluate cardiac mitochondrial toxicity of drugs, during the developmental stage of new drug discovery.

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.

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.

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.

Intrahepatic CD69+Vδ1 T cells re-circulate in the blood of patients with metastatic colorectal cancer and limit tumor progression

Background

More than 50% of all patients with colorectal cancer (CRC) develop liver metastases (CLM), a clinical condition characterized by poor prognosis and lack of reliable prognostic markers. Vδ1 cells are a subset of tissue-resident gamma delta (γδ) T lymphocytes endowed with a broad array of antitumor functions and showing a natural high tropism for the liver. However, little is known about their impact in the clinical outcomes of CLM.

Methods

We isolated human γδ T cells from peripheral blood (PB) and peritumoral (PT) tissue of 93 patients undergone surgical procedures to remove CLM. The phenotype of freshly purified γδ T cells was assessed by multiparametric flow cytometry, the transcriptional profiles by single cell RNA-sequencing, the functional annotations by Gene Ontology enrichment analyses and the clonotype by γδ T cell receptor (TCR)-sequencing.

Results

The microenvironment of CLM is characterized by a heterogeneous immune infiltrate comprising different subsets of γδ tumor-infiltrating lymphocytes (TILs) able to egress the liver and re-circulate in PB. Vδ1 T cells represent the largest population of γδ TILs within the PT compartment of CLM that is greatly enriched in Vδ1 T effector (T_EF) cells expressing constitutive high levels of CD69. These Vδ1 CD69⁺ TILs express a distinct phenotype and transcriptional signature, show high antitumor potential and correlate with better patient clinical outcomes in terms of lower numbers of liver metastatic lesions and longer overall survival (OS). Moreover, intrahepatic CD69⁺ Vδ1 TILs can egress CLM tissue to re-circulate in PB, where they retain a phenotype, transcriptional signature and TCR clonal repertoires resembling their liver origin. Importantly, even the increased frequencies of the CD69⁺ terminally differentiated (T_EMRA) Vδ1 cells in PB of patients with CLM significantly correlate with longer OS. The positive prognostic score of high frequencies of CD69⁺ T_EMRA Vδ1 cells in PB is independent from the neoadjuvant chemotherapy and immunotherapy regimens administered to patients with CLM prior surgery.

Conclusions

The enrichment of tissue-resident CD69⁺ Vδ1 T_EMRA cells re-circulating at high frequencies in PB of patients with CLM limits tumor progression and represents a new important clinical tool to either predict the natural history of CLM or develop alternative therapeutic protocols of cellular therapies.

Krüppel-like factor 15 integrated autophagy and gluconeogenesis to maintain glucose homeostasis under 20-hydroxyecdysone regulation

The regulation of glycometabolism homeostasis is vital to maintain health and development of animal and humans; however, the molecular mechanisms by which organisms regulate the glucose metabolism homeostasis from a feeding state switching to a non-feeding state are not fully understood. Using the holometabolous lepidopteran insect Helicoverpa armigera, cotton bollworm, as a model, we revealed that the steroid hormone 20-hydroxyecdysone (20E) upregulated the expression of transcription factor Krüppel-like factor (identified as Klf15) to promote macroautophagy/autophagy, apoptosis and gluconeogenesis during metamorphosis. 20E via its nuclear receptor EcR upregulated Klf15 transcription in the fat body during metamorphosis. Knockdown of Klf15 using RNA interference delayed pupation and repressed autophagy and apoptosis of larval fat body during metamorphosis. KLF15 promoted autophagic flux and transiting to apoptosis. KLF15 bound to the KLF binding site (KLF bs) in the promoter of Atg8 (autophagy-related gene 8/LC3) to upregulate Atg8 expression. Knockdown Atg8 reduced free fatty acids (FFAs), glycerol, free amino acids (FAAs) and glucose levels. However, knockdown of Klf15 accumulated FFAs, glycerol, and FAAs. Glycolysis was switched to gluconeogenesis, trehalose and glycogen synthesis were changed to degradation during metamorphosis, which were accompanied by the variation of the related genes expression. KLF15 upregulated phosphoenolpyruvate carboxykinase (Pepck) expression by binding to KLF bs in the Pepck promoter for gluconeogenesis, which utilised FFAs, glycerol, and FAAs directly or indirectly to increase glucose in the hemolymph. Taken together, 20E via KLF15 integrated autophagy and gluconeogenesis by promoting autophagy-related and gluconeogenesis-related genes expression.

Glucose is the direct substrate for energy production in animal and humans. Autophagy and gluconeogenesis are known to help organisms maintaining energy substrates; however, the mechanism of integration of autophagy and gluconeogenesis is unclear. Holometabolous insects stop feeding during metamorphosis under steroid hormone 20-hydroxyecdysone (20E) regulation, providing a good model for the study. Using lepidopteran insect Helicoverpa armigera, cotton bollworm, as a model, we revealed that Krüppel-like factor 15 (KLF15) integrated autophagy and gluconeogenesis to maintain glucose homeostasis under 20E regulation. 20E increased Klf15 expression, and KLF15 in turn promoted autophagy-related and gluconeogenesis-related genes expression during metamorphosis. Autophagy and apoptosis of the fat body provided substrates for gluconeogenesis. This work clarified the important functions and mechanisms of KLF15 in autophagy and glycometabolism reprogramming for glucose homeostasis after feeding stop during insect metamorphosis.

Single cell transcriptomics reveals unique metabolic profiles of ependymoma subgroups

OBJECTIVE

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

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.

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.

RoDiCE: robust differential protein co-expression analysis for cancer complexome

MOTIVATION

The full spectrum of abnormalities in cancer-associated protein complexes remains largely unknown. Comparing the co-expression structure of each protein complex between tumor and healthy cells may provide insights regarding cancer-specific protein dysfunction. However, the technical limitations of mass spectrometry-based proteomics, including contamination with biological protein variants, causes noise that leads to non-negligible over- (or under-) estimating co-expression.

RESULTS

We propose a robust algorithm for identifying protein complex aberrations in cancer based on differential protein co-expression testing. Our method based on a copula is sufficient for improving identification accuracy with noisy data compared to conventional linear correlation-based approaches. As an application, we use large-scale proteomic data from renal cancer to show that important protein complexes, regulatory signaling pathways and drug targets can be identified. The proposed approach surpasses traditional linear correlations to provide insights into higher-order differential co-expression structures.

AVAILABILITY AND IMPLEMENTATION

https://github.com/ymatts/RoDiCE.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

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

Simple Summary

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

Abstract

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

OUP accepted manuscript

Background

The glioblastoma-amplified sequence (GBAS) is a newly identified gene that is amplified in approximately 40% of glioblastomas. This article probes into the expression, prognostic significance, and possible pathways of GBAS in ovarian cancer (OC).

Method

Immunohistochemical methods were used to evaluate the expression level of GBAS in OC and its relationship with clinicopathological characteristics and prognosis. Glioblastoma-amplified sequence shRNA was designed to transfect into OC cell lines to silence GBAS expression, then detect the proliferation, apoptosis, and migration ability of the cell. Furthermore, an in vitro tumor formation experiment in mice was constructed to prove the effect of GBAS expression on the growth of OC in vivo. To further study the regulation mechanism of GBAS, we performed co-immunoprecipitation (Co-IP) and shotgun LC-MS mass spectrometry identification.

Results

Immunohistochemistry indicated that GBAS was markedly overexpressed in OC compared with normal ovarian tissue and was associated with lymph node metastasis. Inhibition of GBAS expression can significantly reduce OC cell proliferation, colony formation, promote cell apoptosis, and reduce the ability of cell migration and invasion. In vivo tumor formation experiments showed that the size and weight of tumors in mice after GBAS expression knockdown was significantly smaller. Glioblastoma-amplified sequence may be combined with elongation factor 1 alpha 1 (eEF1A1) to achieve its regulation in OC. Bioinformatics analysis data indicate that GBAS may be a key regulator of mitochondria-associated pathways, therefore controlling cancer progression. MicroRNA-27b, MicroRNA-23a, and MicroRNA-590 may directly targeting GBAS affects the biological behavior of OC cells.

Conclusion

The glioblastoma-amplified sequence may regulate the proliferation and metastasis of OC cells by combining with eEF1A1.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

The genetic source tracking of human urinary exosomes

The genetic origins of nanoscale extracellular vesicles in our body fluids remains unclear. Here, we perform a tracking analysis of urinary exosomes via RNA sequencing, revealing that urine exosomes mostly express tissue-specific genes for the bladder and have close cell-genetic relationships to the endothelial cell, basal cell, monocyte, and dendritic cell. Tracking the differentially expressed genes of cancers and corresponding enrichment analysis show urine exosomes are intensively involved in immune activities, indicating that they may be harnessed as reliable biomarkers of noninvasive liquid biopsy in cancer genomic diagnostics and precision medicine.

New insights into genetic characteristics between multiple myeloma and COVID-19: An integrative bioinformatics analysis of gene expression omnibus microarray and the cancer genome atlas data

Background

Multiple myeloma (MM) is a hematological malignancy. Coronavirus disease 2019 (COVID‐19) infection correlates with MM features. This study aimed to identify MM prognostic biomarkers with potential association with COVID‐19.

Methods

Differentially expressed genes (DEGs) in five MM data sets (GSE47552, GSE16558, GSE13591, GSE6477, and GSE39754) with the same expression trends were screened out. Functional enrichment analysis and the protein‐protein interaction network were performed for all DEGs. Prognosis‐associated DEGs were screened using the stepwise Cox regression analysis in the cancer genome atlas (TCGA) MMRF‐CoMMpass cohort and the GSE24080 data set. Prognosis‐associated DEGs associated with COVID‐19 infection in the GSE164805 data set were also identified.

Results

A total of 98 DEGs with the same expression trends in five data sets were identified, and 83 DEGs were included in the protein‐protein interaction network. Cox regression analysis identified 16 DEGs were associated with MM prognosis in the TCGA cohort, and only the cytochrome c oxidase subunit 6C (COX6C) gene (HR = 1.717, 95% CI 1.231–2.428, p = .002) and the nucleotide‐binding oligomerization domain containing 2 (NOD2) gene (HR = 0.882, 95% CI 0.798–0.975, p = .014) were independent factors related to MM prognosis in the GSE24080 data set. Both of them were downregulated in patients with mild COVID‐19 infection compared with controls but were upregulated in patients with severe COVID‐19 compared with patients with mild illness.

Conclusions

The NOD2 and COX6C genes might be used as prognostic biomarkers in MM. The two genes might be associated with the development of COVID‐19 infection.

NKX3.1 Localization to Mitochondria Suppresses Prostate Cancer Initiation

Mitochondria provide the first line of defense against the tumor-promoting effects of oxidative stress. Here we show that the prostate-specific homeoprotein NKX3.1 suppresses prostate cancer initiation by protecting mitochondria from oxidative stress. Integrating analyses of genetically engineered mouse models, human prostate cancer cells, and human prostate cancer organotypic cultures, we find that, in response to oxidative stress, NKX3.1 is imported to mitochondria via the chaperone protein HSPA9, where it regulates transcription of mitochondrial-encoded electron transport chain (ETC) genes, thereby restoring oxidative phosphorylation and preventing cancer initiation. Germline polymorphisms of NKX3.1 associated with increased cancer risk fail to protect from oxidative stress or suppress tumorigenicity. Low expression levels of NKX3.1 combined with low expression of mitochondrial ETC genes are associated with adverse clinical outcome, whereas high levels of mitochondrial NKX3.1 protein are associated with favorable outcome. This work reveals an extranuclear role for NKX3.1 in suppression of prostate cancer by protecting mitochondrial function. SIGNIFICANCE: Our findings uncover a nonnuclear function for NKX3.1 that is a key mechanism for suppression of prostate cancer. Analyses of the expression levels and subcellular localization of NKX3.1 in patients at risk of cancer progression may improve risk assessment in a precision prevention paradigm, particularly for men undergoing active surveillance.See related commentary by Finch and Baena, p. 2132.This article is highlighted in the In This Issue feature, p. 2113.

miQC: An adaptive probabilistic framework for quality control of single-cell RNA-sequencing data

Single-cell RNA-sequencing (scRNA-seq) has made it possible to profile gene expression in tissues at high resolution. An important preprocessing step prior to performing downstream analyses is to identify and remove cells with poor or degraded sample quality using quality control (QC) metrics. Two widely used QC metrics to identify a 'low-quality' cell are (i) if the cell includes a high proportion of reads that map to mitochondrial DNA (mtDNA) encoded genes and (ii) if a small number of genes are detected. Current best practices use these QC metrics independently with either arbitrary, uniform thresholds (e.g. 5%) or biological context-dependent (e.g. species) thresholds, and fail to jointly model these metrics in a data-driven manner. Current practices are often overly stringent and especially untenable on certain types of tissues, such as archived tumor tissues, or tissues associated with mitochondrial function, such as kidney tissue [1]. We propose a data-driven QC metric (miQC) that jointly models both the proportion of reads mapping to mtDNA genes and the number of detected genes with mixture models in a probabilistic framework to predict the low-quality cells in a given dataset. We demonstrate how our QC metric easily adapts to different types of single-cell datasets to remove low-quality cells while preserving high-quality cells that can be used for downstream analyses. Our software package is available at https://bioconductor.org/packages/miQC.

A pig BodyMap transcriptome reveals diverse tissue physiologies and evolutionary dynamics of transcription

A comprehensive transcriptomic survey of pigs can provide a mechanistic understanding of tissue specialization processes underlying economically valuable traits and accelerate their use as a biomedical model. Here we characterize four transcript types (lncRNAs, TUCPs, miRNAs, and circRNAs) and protein-coding genes in 31 adult pig tissues and two cell lines. We uncover the transcriptomic variability among 47 skeletal muscles, and six adipose depots linked to their different origins, metabolism, cell composition, physical activity, and mitochondrial pathways. We perform comparative analysis of the transcriptomes of seven tissues from pigs and nine other vertebrates to reveal that evolutionary divergence in transcription potentially contributes to lineage-specific biology. Long-range promoter–enhancer interaction analysis in subcutaneous adipose tissues across species suggests evolutionarily stable transcription patterns likely attributable to redundant enhancers buffering gene expression patterns against perturbations, thereby conferring robustness during speciation. This study can facilitate adoption of the pig as a biomedical model for human biology and disease and uncovers the molecular bases of valuable traits.

A comprehensive transcriptomic survey of the pig could enable mechanistic understanding of tissue specialization and accelerate its use as a biomedical model. Here the authors characterize four distinct transcript types in 31 adult pig tissues to dissect their distinct structural and transcriptional features and uncover transcriptomic variability related to tissue physiology.

Single-cell analyses of renal cell cancers reveal insights into tumor microenvironment, cell of origin, and therapy response

Renal cell carcinomas (RCCs) are heterogeneous malignancies thought to arise from kidney tubular epithelial cells, and clear cell RCC is the most common entity. This study demonstrates that cell atlases generated from benign kidney and two common RCCs using single-cell RNA sequencing can predict putative cells of origin for more than 10 RCC subtypes. A focused analysis of distinct cell-type compartments reveals the potential role of tumor epithelia in promoting immune infiltration and other molecular attributes of the tumor microenvironment. Finally, an observed association between the lack of immunotherapy response and endothelial cell fraction has important clinical implications. The current study, therefore, significantly contributes toward understanding disease ontogenies and the molecular dynamics of tumor epithelia and the microenvironment.

Diverse subtypes of renal cell carcinomas (RCCs) display a wide spectrum of histomorphologies, proteogenomic alterations, immune cell infiltration patterns, and clinical behavior. Delineating the cells of origin for different RCC subtypes will provide mechanistic insights into their diverse pathobiology. Here, we employed single-cell RNA sequencing (scRNA-seq) to develop benign and malignant renal cell atlases. Using a random forest model trained on this cell atlas, we predicted the putative cell of origin for more than 10 RCC subtypes. scRNA-seq also revealed several attributes of the tumor microenvironment in the most common subtype of kidney cancer, clear cell RCC (ccRCC). We elucidated an active role for tumor epithelia in promoting immune cell infiltration, potentially explaining why ccRCC responds to immune checkpoint inhibitors, despite having a low neoantigen burden. In addition, we characterized an association between high endothelial cell types and lack of response to immunotherapy in ccRCC. Taken together, these single-cell analyses of benign kidney and RCC provide insight into the putative cell of origin for RCC subtypes and highlight the important role of the tumor microenvironment in influencing ccRCC biology and response to therapy.

The Clinicopathologic and Molecular Landscape of Clear Cell Papillary Renal Cell Carcinoma: Implications in Diagnosis and Management

BACKGROUND

Clear cell papillary renal cell carcinoma (CCPRCC) is a recently described tumor entity. Several questions remain about its epidemiology, molecular features, and clinical behavior.

OBJECTIVE

To comprehensively evaluate clinicopathologic and molecular features of CCPRCC, and compare it with more common kidney cancer subtypes.

DESIGN, SETTING, AND PARTICIPANTS

We identified 89 CCPRCC patients and compared their clinicopathologic features with 1120 localized clear cell renal cell carcinoma (ccRCC) and 129 type 1 papillary renal cell carcinoma (pRCC) patients.

OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS

Nonparametric statistical testing was used to compare relevant features between tumor types. Overall, cancer-specific survival (CSS) and metastasis-free survival estimates were calculated from initial diagnosis using the Kaplan-Meier method. Patients with ipsilateral multifocal disease were explored further. A subset of CCPRCC tumors underwent genomic analysis and were compared with other RCC subtypes.

RESULTS AND LIMITATIONS

A higher proportion of female (45% vs 32%) and African-American (19% vs 3%) patients were observed in the CCPRCC cohort than in the ccRCC and pRCC cohorts. CCPRCC tumors also had increased odds of presenting with additional ipsilateral masses (odds ratio [OR]: 4.41 [confidence interval {CI}: 2.34, 8.15], p < 0.001) and bilateral disease (OR: 4.80 [CI: 2.40, 9.59], p < 0.001) compared with ccRCC tumors. On molecular analysis, CCPRCC tumors showed fewer somatic aberrations and a greater degree of mitochondrial DNA depletion. In multifocal CCPRCC tumors, histologic concordance among the different renal cell carcinoma masses was estimated at 44% (7/16), and none of the individuals presenting exclusively with CCPRCC tumors developed metastatic disease after 5 yr. In contrast, multifocal tumors with CCPRCC and other nonconcordant histologies were more likely to experience adverse outcomes (CSS, log rank p = 0.034).

CONCLUSIONS

CCPRCC is characterized by distinct molecular and epidemiologic features that could be used to refine current diagnostic approaches. Although their clinical course is generally indolent, multifocal CCPRCC tumors represent a unique diagnostic challenge. In this context, single-mass biopsies could miss concomitant aggressive disease, with a potential negative impact on patient outcomes. Furthermore, high discordance rates in multifocal CCPRCC tumors have important clinical implications in management.

PATIENT SUMMARY

We explored the molecular and clinical features of clear cell papillary renal cell carcinoma (CCPRCC) relative to other kidney cancer subtypes. While CCPRCC generally conveys a good prognosis, additional caution should be taken when it is diagnosed using biopsy if multiple kidney masses are present.

Mitochondrial miR-181a-5p promotes glucose metabolism reprogramming in liver cancer by regulating the electron transport chain

Liver cancer and other malignant tumor cells rely on the glycolytic pathway to obtain energy (i.e. the Warburg effect); however, the underlying mechanism is unclear. Mitochondria are sites of oxidative phosphorylation and adenosine triphosphate (ATP) production. The 13 constituent respiratory chain proteins encoded by the mitochondrial genome (namely, mtDNA) play essential roles. We found that in human hepatocellular carcinoma (HCC) tissues, 11 out of the 13 mtDNA-encoded genes exhibited decreased mRNA levels and 5 genes displayed decreased protein levels, including the cytochrome B (mt-CYB) and cytochrome C oxidase II (mt-CO2) genes. Mitochondrial gene sequencing revealed abnormalities in the levels of a large number of mitochondrial miRNAs (mitomiRs). MicroRNA-181a-5p (mir-181a-5p), which potentially targets genes encoding mt-CYB and mt-CO2 protein, was screened out from 549 downregulated mitomiRs via bioinformatic analysis. After overexpression of mitomiR-181a-5p, mt-CYB and mt-CO2 levels were reduced in HCC cells, and the mitochondrial membrane potential (MMP) maintained by the electron transport chain (ETC) was decreased. Furthermore, the expression of hexokinase 2 (HK2) and glucose transporter type 1 (GLUT1) was upregulated, accompanied by elevated glucose, lactic acid release, and activity of lactate dehydrogenase (LDH). In vivo experiments confirmed that constitutive mitomiR-181a-5p expression caused reprogramming of glucose metabolism and promoted tumor growth and early lung metastasis in liver cancer. In summary, the present study reveals the important role of mitomiRs in glucose metabolism reprogramming in liver cancer, which is of considerable value in exploring new therapeutic targets for HCC.

Mitochondrial DNA alterations underlie an irreversible shift to aerobic glycolysis in fumarate hydratase-deficient renal cancer

Understanding the mechanisms of the Warburg shift to aerobic glycolysis is critical to defining the metabolic basis of cancer. Hereditary leiomyomatosis and renal cell carcinoma (HLRCC) is an aggressive cancer characterized by biallelic inactivation of the gene encoding the Krebs cycle enzyme fumarate hydratase, an early shift to aerobic glycolysis, and rapid metastasis. We observed impairment of the mitochondrial respiratory chain in tumors from patients with HLRCC. Biochemical and transcriptomic analyses revealed that respiratory chain dysfunction in the tumors was due to loss of expression of mitochondrial DNA (mtDNA)-encoded subunits of respiratory chain complexes, caused by a marked decrease in mtDNA content and increased mtDNA mutations. We demonstrated that accumulation of fumarate in HLRCC tumors inactivated the core factors responsible for replication and proofreading of mtDNA, leading to loss of respiratory chain components, thereby promoting the shift to aerobic glycolysis and disease progression in this prototypic model of glucose-dependent human cancer.

Reprogramming of aerobic glycolysis in non-transformed mouse liver with pyruvate dehydrogenase complex deficiency

The Pyruvate Dehydrogenase Complex (PDC), a key enzyme in glucose metabolism, catalyzes an irreversible oxidative decarboxylation reaction of pyruvate to acetyl‐CoA, linking the cytosolic glycolytic pathway to mitochondrial tricarboxylic acid cycle and oxidative phosphorylation. Earlier we reported a down‐regulation of several key hepatic lipogenic enzymes and their upstream regulators in liver‐specific PDC‐deficient mouse (L‐PDCKO model by deleting the Pdha1 gene). In this study we investigated gene expression profiles of key glycolytic enzymes and other proteins that respond to various metabolic stresses in liver from L‐PDCKO mice. Transcripts of several, such as hexokinase 2, phosphoglycerate kinase 1, pyruvate kinase muscle‐type 2, and lactate dehydrogenase B as well as those for the nonglycolysis‐related proteins, CD‐36, C/EBP homologous protein, and peroxisome proliferator‐activated receptor γ, were up‐regulated in L‐PDCKO liver whereas hypoxia‐induced factor‐1α, pyruvate dehydrogenase kinase 1 and Sirtuin 1 transcripts were down‐regulated. The protein levels of pyruvate kinase muscle‐type 2 and lactate dehydrogenase B were increased whereas that of lactate dehydrogenase A was decreased in PDC‐deficient mouse liver. Analysis of endoplasmic reticulum and oxidative stress indicators suggests that the L‐PDCKO liver showed evidence of the former but not the latter. These findings indicate that (i) liver‐specific PDC deficiency is sufficient to induce “aerobic glycolysis characteristic” in mouse liver, and (ii) the mechanism(s) responsible for these changes appears distinct from that which induces the Warburg effect in some cancer cells.

Targeting Metabolic Plasticity and Flexibility Dynamics for Cancer Therapy

Cancer cells continuously rewire their metabolism to fulfill their need for rapid growth and survival while subject to changes in environmental cues. Thus, a vital component of a cancer cell lies in its metabolic adaptability. The constant demand for metabolic alterations requires flexibility, that is, the ability to utilize different metabolic substrates; as well as plasticity, that is, the ability to process metabolic substrates in different ways. In this review, we discuss how dynamic changes in cancer metabolism affect tumor progression and the consequential implications for cancer therapy. SIGNIFICANCE: Recognizing cancer dynamic metabolic adaptability as an entity can lead to targeted therapy that is expected to decrease drug resistance.

Mitophagy-associated genes PINK1 and PARK2 are independent prognostic markers of survival in papillary renal cell carcinoma and associated with aggressive tumor behavior

The aim of this study was to investigate the mitophagy-related genes PINK1 and PARK2 in papillary renal cell carcinoma and their association with prognosis. In silico data of PINK1 and PARK2 were analyzed in TCGA cohorts of papillary renal cell carcinoma comprising 290 tumors and 33 corresponding non-neoplastic renal tissues. Protein expression data from a cohort of 95 papillary renal cell carcinoma patients were analyzed and associated with clinical-pathological parameters including survival. PINK1 and PARK2 were significantly downregulated in papillary renal cell carcinoma at transcript and protein levels. Reduced transcript levels of PINK1 and PARK2 were negatively associated with overall survival (p < 0.05). At the protein level, PARK2 and PINK1 expression were positively correlated (correlation coefficient 0.286, p = 0.04) and reduced PINK1 protein expression was prognostic for shorter survival. Lower PINK1 protein levels were found in tumors with metastases at presentation and in tumors of higher pT-stages. The multivariate analysis revealed mRNA expression of PINK1 and PARK2 as well as PINK1 protein expression as independent prognostic factors for shorter overall survival. The downregulation of PINK1 is a strong predictor of poor survival in papillary renal cell carcinoma. Immunohistochemical PINK1 expression in resected pRCC should be considered as an additional prognostic marker for routine practice.