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

Mitochondrial DNA content and methylation in sperm of patients with asthenozoospermia

PURPOSE

The aim of the current study was to investigate the mtDNA methylation levels and mtDNA copy numbers in the sperm of patients with asthenozoospermia and compare them to those observed in controls with normozoospermia.

METHODS

Pyrosequencing analysis of the methylation levels of the mitochondrial D-loop and MT-CO1/chr1:631,907-632083/chrX:26,471,887-126,472,063 (hereinafter referred to as "MT-CO1-AVG") region and quantitative PCR analysis of the mtDNA copy number were performed on sperm from 30 patients with asthenozoospermia and 30 controls with normozoospermia.

RESULTS

Compared with those of controls with normozoospermia, the methylation levels of D-loop and MT-CO1-AVG regions and mtDNA copy number were significantly higher in patients with asthenozoospermia. The methylation level of the D-loop region in patients with asthenozoospermia and controls with normozoospermia and that of MT-CO1-AVG region in patients with asthenozoospermia showed a decreasing tendency with increasing total sperm motility. A significant inverse correlation between the mtDNA copy number and total sperm motility was observed in patients with asthenozoospermia but not in controls with normozoospermia. In patients with asthenozoospermia, but not in controls with normozoospermia, we observed a significant inverse correlation between D-loop methylation levels and mtDNA copy number, while no significant correlation was observed between MT-CO1-AVG methylation levels and mtDNA copy number.

CONCLUSION

These results reveal the occurrence of mtDNA methylation in human sperm and altered D-loop and MT-CO1-AVG methylation levels in patients with asthenozoospermia. Additional research is needed to determine the function of these features in the etiology and course of asthenozoospermia.

Immune infiltration and prognosis in gastric cancer: role of NAD+ metabolism-related markers

Background

This study endeavored to develop a nicotinamide adenine dinucleotide (NAD+) metabolism-related biomarkers in gastric cancer (GC), which could provide a theoretical foundation for prognosis and therapy of GC patients.

Methods

In this study, differentially expressed genes (DEGs1) between GC and paraneoplastic tissues were overlapped with NAD+ metabolism-related genes (NMRGs) to identify differentially expressed NMRGs (DE-NMRGs). Then, GC patients were divided into high and low score groups by gene set variation analysis (GSVA) algorithm for differential expression analysis to obtain DEGs2, which was overlapped with DEGs1 for identification of intersection genes. These genes were further analyzed using univariate Cox and least absolute shrinkage and selection operator (LASSO) regression analyses to obtain prognostic genes for constructing a risk model. Enrichment and immune infiltration analyses further investigated investigate the different risk groups, and qRT-PCR validated the prognostic genes.

Results

Initially, we identified DE-NMRGs involved in NAD biosynthesis, with seven (DNAJB13, CST2, THPO, CIDEA, ONECUT1, UPK1B and SNCG) showing prognostic significance in GC. Subsequent, a prognostic model was constructed in which the risk score, derived from the expression profiles of these genes, along with gender, emerged as robust independent predictors of patient outcomes in GC. Enrichment analysis linked high-risk patients to synaptic membrane pathways and low-risk to the CMG complex pathway. Tumor immune infiltration analysis revealed correlations between risk scores and immune cell abundance, suggesting a relationship between NAD+ metabolism and immune response in GC. The prognostic significance of our identified genes was validated by qRT-PCR, which confirmed their upregulated expression in GC tissue samples.

Conclusion

In this study, seven NAD+ metabolism-related markers were established, which is of great significance for the development of prognostic molecular biomarkers and clinical prognosis prediction for gastric cancer patients.

Translational insights into the hormetic potential of carbon dioxide: from physiological mechanisms to innovative adjunct therapeutic potential for cancer

Background

Carbon dioxide (CO2), traditionally viewed as a mere byproduct of cellular respiration, plays a multifaceted role in human physiology beyond simple elimination through respiration. CO2 may regulate the tumor microenvironment by significantly affecting the release of oxygen (O2) to tissues through the Bohr effect and by modulating blood pH and vasodilation. Previous studies suggest hypercapnia (elevated CO2 levels) might trigger optimized cellular mechanisms with potential therapeutic benefits. The role of CO2 in cellular stress conditions within tumor environments and its impact on O2 utilization offers a new investigative area in oncology.

Objectives

This study aims to explore CO2's role in the tumor environment, particularly how its physiological properties and adaptive responses can influence therapeutic strategies.

Methods

By applying a structured translational approach using the Work Breakdown Structure method, the study divided the analysis into six interconnected work packages to comprehensively analyze the interactions between carbon dioxide and the tumor microenvironment. Methods included systematic literature reviews, data analyses, data integration for identifying critical success factors and exploring extracellular environment modulation. The research used SMART criteria for assessing innovation and the applicability of results.

Results

The research revealed that the human body's adaptability to hypercapnic conditions could potentially inform innovative strategies for manipulating the tumor microenvironment. This could enhance O2 utilization efficiency and manage adaptive responses to cellular stress. The study proposed that carbon dioxide's hormetic potential could induce beneficial responses in the tumor microenvironment, prompting clinical protocols for experimental validation. The research underscored the importance of pH regulation, emphasizing CO2 and carbonic acid's role in modulating metabolic and signaling pathways related to cancer.

Conclusion

The study underscores CO2 as vital to our physiology and suggests potential therapeutic uses within the tumor microenvironment. pH modulation and cellular oxygenation optimization via CO2 manipulation could offer innovative strategies to enhance existing cancer therapies. These findings encourage further exploration of CO2's therapeutic potential. Future research should focus on experimental validation and exploration of clinical applications, emphasizing the need for interdisciplinary and collaborative approaches to tackle current challenges in cancer treatment.

A mitochondria-related genes associated neuroblastoma signature - based on bulk and single-cell transcriptome sequencing data analysis, and experimental validation

Background

Neuroblastoma (NB), characterized by its marked heterogeneity, is the most common extracranial solid tumor in children. The status and functionality of mitochondria are crucial in regulating NB cell behavior. While the significance of mitochondria-related genes (MRGs) in NB is still missing in key knowledge.

Materials and methods

This study leverages consensus clustering and machine learning algorithms to construct and validate an MRGs-related signature in NB. Single-cell data analysis and experimental validation were employed to characterize the pivotal role of FEN1 within NB cells.

Results

MRGs facilitated the classification of NB patients into 2 distinct clusters with considerable differences. The constructed MRGs-related signature and its quantitative indicators, mtScore and mtRisk, effectively characterize the MRGs-related patient clusters. Notably, the MRGs-related signature outperformed MYCN in predicting NB patient prognosis and was adept at representing the tumor microenvironment (TME), tumor cell stemness, and sensitivity to the chemotherapeutic agents Cisplatin, Topotecan, and Irinotecan. FEN1, identified as the most contributory gene within the MRGs-related signature, was found to play a crucial role in the communication between NB cells and the TME, and in the developmental trajectory of NB cells. Experimental validations confirmed FEN1's significant influence on NB cell proliferation, apoptosis, cell cycle, and invasiveness.

Conclusion

The MRGs-related signature developed in this study offers a novel predictive tool for assessing NB patient prognosis, immune infiltration, stemness, and chemotherapeutic sensitivity. Our findings unveil the critical function of FEN1 in NB, suggesting its potential as a therapeutic target.

Identification of VEGFs-related gene signature for predicting microangiogenesis and hepatocellular carcinoma prognosis

Microangiogenesis is an important prognostic factor in various cancers, including hepatocellular carcinoma (HCC). The Vascular Endothelial Growth Factor (VEGF) has been shown to contribute to tumor angiogenesis. Recently, several studies have investigated the regulation of VEGF production by a single gene, with few researchers exploring all genes that affect VEGF production. In this study, we comprehensively analyzed all genes affecting VEGF production in HCC and developed a risk model and gene-based risk score based on VEGF production. Moreover, the model's predictive capacity on prognosis of HCCs was verified using training and validation datasets. The developed model showed good prediction of the overall survival rate. Patients with a higher risk score experienced poor outcomes compared to those with a lower risk score. Furthermore, we identified the immunological causes of the poor prognosis of patients with high-risk scores comparing with those with low-risk scores.

Construction and Analysis of a Mitochondrial Metabolism-Related Prognostic Model for Breast Cancer to Evaluate Survival and Immunotherapy

As one of the most prevalent malignancies among women, breast cancer (BC) is tightly linked to metabolic dysfunction. However, the correlation between mitochondrial metabolism-related genes (MMRGs) and BC remains unclear. The training and validation datasets for BC were obtained from The Cancer Genome Atlas and Gene Expression Omnibus databases, respectively. MMRG-related data were obtained from the Molecular Signatures Database. A risk score prognostic model incorporating MMRGs was established based on univariate, LASSO, and multivariate Cox regression analyses. Independent factors affecting BC prognosis were identified through regression analysis and presented in a nomogram. Single-sample gene set enrichment analysis was employed to assess the immune levels of high-risk (HR) and low-risk (LR) groups. The sensitivity of BC patients in the two groups to common anti-tumor drugs was evaluated by utilizing the Genomics of Drug Sensitivity in Cancer database. 12 MMRGs significantly associated with survival were selected from 1234 MMRGs. A 12-gene risk score prognostic model was built. In the multivariate regression analysis incorporating classical clinical factors, the MMRG-related risk score remained an independent prognostic factor. As revealed by tumor immune microenvironment analysis, the LR group with higher survival rates had elevated immune levels. The drug sensitivity results unmasked that the LR group demonstrated higher sensitivity to Irinotecan, Nilotinib, and Oxaliplatin, while the HR group demonstrated higher sensitivity to Lapatinib. The development of MMRG characteristics provides a comprehensive understanding of mitochondrial metabolism in BC, aiding in the prediction of prognosis and tumor microenvironment, and offering promising therapeutic choices for BC patients with different MMRG risk scores.

Identification of cuproptosis and ferroptosis-related subgroups and development of a signature for predicting prognosis and tumor microenvironment landscape in hepatocellular carcinoma

Background

Ferroptosis and cuproptosis play a crucial role in the progression and dissemination of hepatocellular carcinoma (HCC). The primary objective of this study was to develop a unique scoring system for predicting the prognosis and immunological landscape of HCC based on ferroptosis-related genes (FRGs) and cuproptosis-related genes (CRGs).

Methods

As the training cohort, we assembled a novel HCC cohort by merging gene expression data and clinical data from The Cancer Genome Atlas (TCGA) database, and Gene Expression Omnibus (GEO) database. The validation cohort consisted of 230 HCC cases taken from the International Cancer Genome Consortium (ICGC) database. Multiple genomic characteristics, such as tumor mutation burden (TMB), and copy number variations were analyzed concurrently. On the basis of the expression of CRGs and FRGs, patients were classified into cuproptosis and ferroptosis subtypes. Then, we constructed a risk model using least absolute shrinkage and selection operator (LASSO) analysis and Cox regression analysis based on ferroptosis and cuproptosis-related differentially expressed genes (DEGs). Patients were separated into two groups according to median risk score. We compared the immunophenotype, tumor microenvironment (TME), cancer stem cell index, and treatment sensitivity of two groups.

Results

Three subtypes of ferroptosis and two subtypes of cuproptosis were identified among the patients. A greater likelihood of survival (P<0.05) was expected for patients in FRGcluster B and CRGcluster B. After that, a confirmed risk signature for ferroptosis and cuproptosis was developed and tested. Patients in the low-risk group had significantly higher survival rates than those in the high-risk group, according to our study (P<0.001). There was also a strong correlation between the signature and other variables including immunophenoscore, TMB, cancer stem cell index, immunological checkpoint genes, and sensitivity to chemotherapeutics.

Conclusions

Through this comprehensive research, we identified a unique risk signature associated with HCC patients' treatment status and prognosis. Our findings highlight FRGs' and CRGs' significance in clinical practice and imply ferroptosis and cuproptosis may be therapeutic targets for HCC patients.

Diagnostic role and immune correlates of programmed cell death-related genes in hepatocellular carcinoma

Programmed cell death (PCD) is thought to have multiple roles in tumors. Here, the roles of PCD-related genes were comprehensively analyzed to evaluate their values in hepatocellular carcinoma (HCC) diagnosis and prognosis. Gene expression and single-cell data of HCC patients, and PCD-related genes were collected from public databases. The diagnostic and prognostic roles of differentially expressed PCD-related genes in HCC were explored by univariate and multivariate Cox regression analyses. Single-cell data were further analyzed for the immune cells and expression of feature genes. Finally, we evaluated the expression of genes by quantitative real-time polymerase chain reaction and Western blot, and the proportion of immune cells was detected by flow cytometry in HCC samples. We obtained 52 differentially expressed PCD-related genes in HCC, based on which the consensus clustering analysis cluster 2 was found to have a worse prognosis than cluster 1. Then 10 feature genes were identified using LASSO analysis, and programmed cell death index (PCDI) was calculated to divided HCC patients into high-PCDI and low-PCDI groups. Worse prognosis was observed in high-PCDI group. Cox regression analysis showed that PCDI is an independent prognostic risk factor for HCC patients. Additionally, SERPINE1 and G6PD of feature genes significantly affect patient survival. Macrophages and Tregs were significantly positively correlated with PCDI. G6PD mainly expressed in macrophages, SERPINE1 mainly expressed in fibroblast. The experimental results confirmed the high expression of SERPINE1 and G6PD in HCC compared with the control, and the infiltration level of macrophages and Treg in HCC was also obviously elevated. PCDI may be a new predictor for the diagnosis of patients with HCC. The association of SERPINE1 and G6PD with the immune environment will provide new clues for HCC therapy.

Development and investigation of metabolism-associated risk assessment models for patients with viral hepatitis

Dysregulation of metabolism plays an important role in the onset and progression of multiple pathogenic diseases, including viral hepatitis. However, a model to predict viral hepatitis risk by metabolic pathways is still lacking. Thus, we developed two risk assessment models for viral hepatitis based on metabolic pathways identified through univariate and least absolute shrinkage and selection operator (LASSO) Cox regression analysis. The first model is designed to assess the progression of the disease by evaluating changes in the Child–Pugh class, hepatic decompensation, and the development of hepatocellular carcinoma. The second model is focused on determining the prognosis of the illness, taking into account the patient’s cancer status. Our models were further validated by Kaplan–Meier plots of survival curves. In addition, we investigated the contribution of immune cells in metabolic processes and identified three distinct subsets of immune cells—CD8+ T cells, macrophages, and NK cells—that have significantly affected metabolic pathways. Specifically, our findings suggest that resting or inactive macrophages and NK cells contribute to maintaining metabolic homeostasis, particularly with regard to lipid and α-amino acid metabolism, thereby potentially reducing the risk of viral hepatitis progression. Moreover, maintaining metabolic homeostasis ensures a balance between killer-proliferative and exhausted CD8+ T cells, which helps in mitigating CD8+ T cell-mediated liver damage while preserving energy reserves. In conclusion, our study offers a useful tool for early disease detection in viral hepatitis patients through metabolic pathway analysis and sheds light on the immunological understanding of the disease through the examination of immune cell metabolic disorders.

Ferroptosis and cuproptosis prognostic signature for prediction of prognosis, immunotherapy and drug sensitivity in hepatocellular carcinoma: development and validation based on TCGA and ICGC databases

Background

Hepatocellular carcinoma (HCC) is a common malignancy. Ferroptosis and cuproptosis promote HCC spread and proliferation. While fewer studies have combined ferroptosis and cuproptosis to construct prognostic signature of HCC. This work attempts to establish a novel scoring system for predicting HCC prognosis, immunotherapy, and medication sensitivity based on ferroptosis-related genes (FRGs) and cuproptosis-related genes (CRGs).

Methods

FerrDb and previous literature were used to identify FRGs. CRGs came from original research. The Cancer Genome Atlas (TCGA) and International Cancer Genome Consortium (ICGC) databases included the HCC transcriptional profile and clinical information [survival time, survival status, age, gender, Tumor Node Metastasis (TNM) stage, etc.]. Correlation, Cox, and least absolute shrinkage and selection operator (LASSO) regression analyses were used to narrow down prognostic genes and develop an HCC risk model. Using "caret", R separated TCGA-HCC samples into a training risk set and an internal test risk set. As external validation, we used ICGC samples. We employed Kaplan-Meier analysis and receiver operating characteristic (ROC) curve to evaluate the model's clinical efficacy. CIBERSORT and TIMER measured immunocytic infiltration in high- and low-risk populations.

Results

TXNRD1 [hazard ratio (HR) =1.477, P<0.001], FTL (HR =1.373, P=0.001), GPX4 (HR =1.650, P=0.004), PRDX1 (HR =1.576, P=0.002), VDAC2 (HR =1.728, P=0.008), OTUB1 (HR =1.826, P=0.002), NRAS (HR =1.596, P=0.005), SLC38A1 (HR =1.290, P=0.002), and SLC1A5 (HR =1.306, P<0.001) were distinguished to build predictive model. In both the model cohort (P<0.001) and the validation cohort (P<0.05), low-risk patients had superior overall survival (OS). The areas under the curve (AUCs) of the ROC curves in the training cohort (1-, 3-, and 5-year AUCs: 0.751, 0.727, and 0.743), internal validation cohort (1-, 3-, and 5-year AUCs: 0.826, 0.624, and 0.589), and ICGC cohort (1-, 3-, and 5-year AUCs: 0.699, 0.702, and 0.568) were calculated. Infiltration of immune cells and immunological checkpoints were also connected with our signature. Treatments with BI.2536, Epothilone.B, Gemcitabine, Mitomycin.C, Obatoclax. Mesylate, and Sunitinib may profit high-risk patients.

Conclusions

We analyzed FRGs and CRGs profiles in HCC and established a unique risk model for treatment and prognosis. Our data highlight FRGs and CRGs in clinical practice and suggest ferroptosis and cuproptosis may be therapeutic targets for HCC patients. To validate the model's clinical efficacy, more HCC cases and prospective clinical assessments are needed.

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

Background

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

Methods

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

Results

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

Conclusion

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

Supplementary Information

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

Overweight/obesity-related transcriptomic signature as a correlate of clinical outcome, immune microenvironment, and treatment response in hepatocellular carcinoma

BACKGROUNDS

The pandemic of overweight and obesity (quantified by body mass index (BMI) ≥ 25) has rapidly raised the patient number of non-alcoholic fatty hepatocellular carcinoma (HCC), and several clinical trials have shown that BMI is associated with the prognosis of HCC. However, whether overweight/obesity is an independent prognostic factor is arguable, and the role of overweight/obesity-related metabolisms in the progression of HCC is scarcely known.

MATERIALS AND METHODS

In the present study, clinical information, mRNA expression profile, and genomic data were downloaded from The Cancer Genome Atlas (TCGA) as a training cohort (TCGA-HCC) for the identification of overweight/obesity-related transcriptome. Machine learning and the Cox regression analysis were conducted for the construction of the overweight/obesity-associated gene (OAG) signature. The Kaplan-Meier curve, receiver operating characteristic (ROC) curve, and the Cox regression analysis were performed to assess the prognostic value of the OAG signature, which was further validated in two independent retrospective cohorts from the International Cancer Genome Consortium (ICGC) and Gene Expression Omnibus (GEO). Subsequently, functional enrichment, genomic profiling, and tumor microenvironment (TME) evaluation were utilized to characterize biological activities associated with the OAG signature. GSE109211 and GSE104580 were retrieved to evaluate the underlying response of sorafenib and transcatheter arterial chemoembolization (TACE) treatment, respectively. The Genomics of Drug Sensitivity in Cancer (GDSC) database was employed for the evaluation of chemotherapeutic response.

RESULTS

Overweight/obesity-associated transcriptome was mainly involved in metabolic processes and noticeably and markedly correlated with prognosis and TME of HCC. Afterward, a novel established OAG signature (including 17 genes, namely, GAGE2D, PDE6A, GABRR1, DCAF8L1, DPYSL4, SLC6A3, MMP3, RIBC2, KCNH2, HTRA3, PDX1, ATHL1, PRTG, SHC4, C21orf29, SMIM32, and C1orf133) divided patients into high and low OAG score groups with distinct prognosis (median overall survival (OS): 24.87 vs. 83.51 months, p < 0.0001), and the values of area under ROC curve (AUC) in predicting 1-, 2-, 3-, and 4-year OS were 0.81, 0.80, 0.83, and 0.85, respectively. Moreover, the OAG score was independent of clinical features and also exhibited a good ability for prognosis prediction in the ICGC-LIHC-JP cohort and GSE54236 dataset. Expectedly, the OAG score was also highly correlated with metabolic processes, especially oxidative-related signaling pathways. Furthermore, abundant enrichment of chemokines, receptors, MHC molecules, and other immunomodulators as well as PD-L1/PD-1 expression among patients with high OAG scores indicated that they might have better responses to immunotherapy. However, probably exclusion of T cells from infiltrating tumors resulting in lower infiltration of effective T cells would restrict immunotherapeutic effects. In addition, the OAG score was significantly associated with the response of sorafenib and TACE treatment.

CONCLUSIONS

Overall, this study comprehensively disclosed the relationship between BMI-guided transcriptome and HCC. Moreover, the OAG signature had the potential clinical applications in the future to promote clinical management and precision medicine of HCC.