Identification of Mitochondrial-Related Prognostic Biomarkers Associated With Primary Bile Acid Biosynthesis and Tumor Microenvironment of Hepatocellular Carcinoma

Mitochondrial-related genes markers that predict survival in patients with head and neck squamous cell carcinoma affect immunomodulation through hypoxia, glycolysis, and angiogenesis pathways

Mitochondria play a crucial role in the occurrence and development of tumors. We used mitochondria-related genes for consistent clustering to identify three stable molecular subtypes of head and neck squamous cell carcinoma (HNSCC) with different prognoses, mutations, and immune characteristics. Significant differences were observed in clinical characteristics, immune microenvironment, immune cell infiltration, and immune cell scores. TP53 was the most significantly mutated; cell cycle-related pathways and tumorigenesis-related pathways were activated in different subtypes. Risk modeling was conducted using a multifactor stepwise regression method, and nine genes were identified as mitochondria-related genes affecting prognosis (DKK1, EFNB2, ITGA5, AREG, EPHX3, CHGB, P4HA1, CCND1, and JCHAIN). Risk score calculations revealed significant differences in prognosis, immune cell scores, immune cell infiltration, and responses to conventional chemotherapy drugs. Glycolysis, angiogenesis, hypoxia, and tumor-related pathways were positively correlated with the RiskScore. Clinical samples were subjected to qPCR to validate the results. In this work, we constructed a prognostic model based on the mitochondrial correlation score, which well reflects the risk and positive factors for the prognosis of patients with HNSCC. This model can be used to guide individualized adjuvant and immunotherapy in patients with HNSCC.

ATAD3A: A Key Regulator of Mitochondria-Associated Diseases

Mitochondrial membrane protein ATAD3A is a member of the AAA-domain-containing ATPases superfamily. It is important for the maintenance of mitochondrial DNA, structure, and function. In recent years, an increasing number of ATAD3A mutations have been identified in patients with neurological symptoms. Many of these mutations disrupt mitochondrial structure, function, and dynamics and are lethal to patients at a young age. Here, we summarize the current understanding of the relationship between ATAD3A and mitochondria, including the interaction of ATAD3A with mitochondrial DNA and mitochondrial/ER proteins, the regulation of ATAD3A in cholesterol mitochondrial trafficking, and the effect of known ATAD3A mutations on mitochondrial function. In the current review, we revealed that the oligomerization and interaction of ATAD3A with other mitochondrial/ER proteins are vital for its various functions. Despite affecting different domains of the protein, nearly all documented mutations observed in ATAD3A exhibit either loss-of-function or dominant-negative effects, potentially leading to disruption in the dimerization of ATAD3A; autophagy; mitophagy; alteration in mitochondrial number, size, and cristae morphology; and diminished activity of mitochondrial respiratory chain complexes I, IV, and V. These findings imply that ATAD3A plays a critical role in mitochondrial dynamics, which can be readily perturbed by ATAD3A mutation variants.

FAM family gene prediction model reveals heterogeneity, stemness and immune microenvironment of UCEC

Background: Endometrial cancer (UCEC) is a highly heterogeneous gynecologic malignancy that exhibits variable prognostic outcomes and responses to immunotherapy. The Familial sequence similarity (FAM) gene family is known to contribute to the pathogenesis of various malignancies, but the extent of their involvement in UCEC has not been systematically studied. This investigation aimed to develop a robust risk profile based on FAM family genes (FFGs) to predict the prognosis and suitability for immunotherapy in UCEC patients. Methods: Using the TCGA-UCEC cohort from The Cancer Genome Atlas (TCGA) database, we obtained expression profiles of FFGs from 552 UCEC and 35 normal samples, and analyzed the expression patterns and prognostic relevance of 363 FAM family genes. The UCEC samples were randomly divided into training and test sets (1:1), and univariate Cox regression analysis and Lasso Cox regression analysis were conducted to identify the differentially expressed genes (FAM13C, FAM110B, and FAM72A) that were significantly associated with prognosis. A prognostic risk scoring system was constructed based on these three gene characteristics using multivariate Cox proportional risk regression. The clinical potential and immune status of FFGs were analyzed using CiberSort, SSGSEA, and tumor immune dysfunction and rejection (TIDE) algorithms. qRT-PCR and IHC for detecting the expression levels of 3-FFGs. Results: Three FFGs, namely, FAM13C, FAM110B, and FAM72A, were identified as strongly associated with the prognosis of UCEC and effective predictors of UCEC prognosis. Multivariate analysis demonstrated that the developed model was an independent predictor of UCEC, and that patients in the low-risk group had better overall survival than those in the high-risk group. The nomogram constructed from clinical characteristics and risk scores exhibited good prognostic power. Patients in the low-risk group exhibited a higher tumor mutational load (TMB) and were more likely to benefit from immunotherapy. Conclusion: This study successfully developed and validated novel biomarkers based on FFGs for predicting the prognosis and immune status of UCEC patients. The identified FFGs can accurately assess the prognosis of UCEC patients and facilitate the identification of specific subgroups of patients who may benefit from personalized treatment with immunotherapy and chemotherapy.

In Vivo Screening Unveils Pervasive RNA-Binding Protein Dependencies in Leukemic Stem Cells and Identifies ELAVL1 as a Therapeutic Target

Acute myeloid leukemia (AML) is fueled by leukemic stem cells (LSC) whose determinants are challenging to discern from hematopoietic stem cells (HSC) or uncover by approaches focused on general cell properties. We have identified a set of RNA-binding proteins (RBP) selectively enriched in human AML LSCs. Using an in vivo two-step CRISPR-Cas9 screen to assay stem cell functionality, we found 32 RBPs essential for LSCs in MLL-AF9;NrasG12D AML. Loss-of-function approaches targeting key hit RBP ELAVL1 compromised LSC-driven in vivo leukemic reconstitution, and selectively depleted primitive malignant versus healthy cells. Integrative multiomics revealed differentiation, splicing, and mitochondrial metabolism as key features defining the leukemic ELAVL1-mRNA interactome with mitochondrial import protein, TOMM34, being a direct ELAVL1-stabilized target whose repression impairs AML propagation. Altogether, using a stem cell-adapted in vivo CRISPR screen, this work demonstrates pervasive reliance on RBPs as regulators of LSCs and highlights their potential as therapeutic targets in AML.

SIGNIFICANCE

LSC-targeted therapies remain a significant unmet need in AML. We developed a stem-cell-adapted in vivo CRISPR screen to identify key LSC drivers. We uncover widespread RNA-binding protein dependencies in LSCs, including ELAVL1, which we identify as a novel therapeutic vulnerability through its regulation of mitochondrial metabolism. This article is highlighted in the In This Issue feature, p. 171.

The pan-cancer analysis of the oncogenic role of FAM72A as a BRCA prognostic biomarker and immunotherapeutic target

In this study, we first comprehensively investigated the expression profile, mutation status, and survival analysis of FAM72A as well as the correlation between FAM72A and DNA damage repair, methylation, and cell stemness analysis using bioinformatics techniques. In addition, we also analyzed the relationship between FAM72A and immune cell infiltration and pathway enrichment. The role of FAM72A in breast cancer (BC) was so conspicuous that we analyzed the prognostic significance and clinicopathological parameter's relevance of FAM72A in BC. We also validated biological functions by applying in vitro experiments. FAM72A was highly expressed in 26 types of a total of 31 cancers, while it expressed low levels in only five cancers. FAM72A expression was relative to clinical stages in nine cancers and has a significant difference in disease-free survival among 31 kinds of cancers. In addition, FAM72A has negatively correlated with cancer-associated fibroblast and endothelial cells in BC but positively correlated with follicular helper T cells. Univariate and multivariate cox regression analyses identified T, N, M, age, and FAM72A expression as independent influences on BC prognosis, so we created a nomogram to predict patient survival benefits. In in vitro experiments, we verified that downregulation of FAM72A not only inhibited cell proliferation, colony formation, cell migration, cell invasion, and G2/M cell cycle transition but also promoted apoptosis of breast invasive carcinoma cells. Our study discovered FAM72A as a clinically meaningful biomarker for prognostic predicting and a guiding target for immune treatment in BC.

Fam72a functions as a cell-cycle-controlled gene during proliferation and antagonizes apoptosis through reprogramming PP2A substrates

The cell cycle is key to life. After decades of research, it is unclear whether any parts of this process have yet to be identified. Fam72a is a poorly characterized gene and is evolutionarily conserved across multicellular organisms. Here, we have found that Fam72a is a cell-cycle-regulated gene that is transcriptionally and post-transcriptionally regulated by FoxM1 and APC/C, respectively. Functionally, Fam72a directly binds to tubulin and both the Aα and B56 subunits of PP2A-B56 to modulate tubulin and Mcl1 phosphorylation, which in turn affects the progression of the cell cycle and signaling of apoptosis. Moreover, Fam72a is involved in early responses to chemotherapy, and it efficiently antagonizes various anticancer compounds such as CDK and Bcl2 inhibitors. Thus, Fam72a switches the tumor-suppressive PP2A to be oncogenic by reprogramming its substrates. These findings identify a regulatory axis of PP2A and a protein member in the cell cycle and tumorigenesis regulatory network in human cells.

Novel Gene Signatures Promote Epithelial-Mesenchymal Transition (EMT) in Glucose Deprivation-Based Microenvironment to Predict Recurrence-Free Survival in Hepatocellular Carcinoma

Background

Current research studies have suggested that glucose deprivation (GD)-based tumor microenvironment (TME) can promote epithelial-mesenchymal transition (EMT) of tumor cells, leading to tumor invasion and metastasis. However, no one has yet studied detailedly the synthetic studies that include GD features in TME with EMT status. In our research, we comprehensively developed and validated a robust signature regarding GD and EMT status to provide prognostic value for patients with liver cancer.

Methods

GD and EMT status were estimated with transcriptomic profiles based on WGCNA and t-SNE algorithms. Two cohorts of training (TCGA_LIHC) and validation (GSE76427) datasets were analyzed with the Cox regression and logistic regression analyses. We identified a 2-mRNA signature to establish a GD-EMT-based gene risk model for the prediction of HCC relapse.

Results

Patients with significant GD-EMT status were divided into two subgroups: GD^low/EMT^low and GD^high/EMT^high, with the latter having significantly worse recurrence-free survival (P < 0.01). We employed the least absolute shrinkage and selection operator (LASSO) technique as a method for HNF4A and SLC2A4 filtering and constructing a risk score for risk stratification. In the multivariate analysis, this risk score predicted recurrence-free survival (RFS) in both the discovery and validation cohorts and remained valid in patients stratified by TNM stage and age at diagnosis. The nomogram that combines risk score and TNM stage as well as age produces improved performance and net benefits in the analysis of calibration and decision curves in training and validation groups.

Conclusions

The GD-EMT-based signature predictive model may provide a prognosis classifier for HCC patients with a high risk of postoperative recurrence to decrease the relapse rate.

Mitochondria dysfunction in airway epithelial cells is associated with type 2-low asthma

Background: Type 2 (T2)-low asthma can be severe and corticosteroid-resistant. Airway epithelial cells play a pivotal role in the development of asthma, and mitochondria dysfunction is involved in the pathogenesis of asthma. However, the role of epithelial mitochondria dysfunction in T2-low asthma remains unknown. Methods: Differentially expressed genes (DEGs) were identified using gene expression omnibus (GEO) dataset GSE4302, which is originated from airway epithelial brushings from T2-high (n = 22) and T2-low asthma patients (n = 20). Gene set enrichment analysis (GSEA) was implemented to analyze the potential biological pathway involved between T2-low and T2-high asthma. T2-low asthma related genes were identified using weighted gene co-expression network analysis (WGCNA). The mitochondria-related genes (Mito-RGs) were referred to the Molecular Signatures Database (MSigDB). T2-low asthma related mitochondria (T2-low-Mito) DEGs were obtained by intersecting the DEGs, T2-low asthma related genes, and Mito-RGs. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) was performed to further explore the potential function of the T2-low-Mito DEGs. In addition, the hub genes were further identified by protein-protein interaction (PPI), and the expressions of hub genes were verified in another GEO dataset GSE67472 and bronchial brushings from patients recruited at Tongji Hospital. Results: Six hundred and ninety-two DEGs, including 107 downregulated genes and 585 upregulated genes were identified in airway epithelial brushings from T2-high and T2-low asthma patients included in GSE4302 dataset. GSEA showed that mitochondrial ATP synthesis coupled electron transport is involved in T2-low asthma. Nine hundred and four T2-low asthma related genes were identified using WGCNA. Twenty-two T2-low-Mito DEGs were obtained by intersecting the DEGs, T2-low asthma and Mito-RGs. The GO enrichment analysis of the T2-low-Mito DEGs showed significant enrichment of mitochondrial respiratory chain complex assembly, and respiratory electron transport chain. PPI network was constructed using 22 T2-low-Mito DEGs, and five hub genes, ATP5G1, UQCR10, NDUFA3, TIMM10, and NDUFAB1, were identified. Moreover, the expression of these hub genes was validated in another GEO dataset, and our cohort of asthma patients. Conclusion: This study suggests that mitochondria dysfunction contributes to T2-low asthma.

Prognostic and Immunological Implications of FAM72A in Pan-Cancer and Functional Validations

The family with sequence similarity 72 Member A (FAM72A) is overexpressed in several types of cancer. However, its contributions to tumorigenesis remain largely unknown. Based on The Cancer Genome Atlas (TCGA) database, FAM72A was upregulated across 33 types of cancer. Accordingly, high levels of FAM72A predicted inferior outcomes in half of the cancer types using survival analysis (the Kaplan-Meier curve and univariate Cox regression model). Receiver operating characteristic (ROC) analysis demonstrated that FAM72A showed high accuracy in distinguishing cancerous tissues from normal ones. FAM72A was correlated with immune and stromal scores and immune cell infiltrations in various tumors. Moreover, FAM72A was also associated with tumor mutation burden (TMB), microsatellite instability (MSI), and immune checkpoint genes. Immunophenoscore (IPS) further validated that the FAM72A^low tumor showed high immunogenicity and tended to respond to anti-PD1/PDL1/PDL2, anti-CTLA4 treatment, and combined immunotherapies. We also investigated the functional role of FAM72A in lung adenocarcinoma (LUAD). In vitro studies demonstrated that the ectopic expression of FAM72A accelerated the proliferation and migration of NSCLC cells, whereas silencing FAM72A showed the opposite effects on them. In short, FAM72A had prognostic potential and correlated with tumor immunogenicity in various tumors. Functional analysis indicated that FAM72A is an oncogene in LUAD.

Integrated systemic analysis of FAM72A to identify its clinical relevance, biological function, and relationship to drug sensitivity in hepatocellular carcinoma

Background

The family with sequence similarity 72 member A (FAM72A) protein has been identified as an effector of multiple pathological processes in many cancers. The value of FAM72A in HCC remains largely unknown.

Methods

Data from TCGA-LIHC, ICGC-LIRI-JP, IMvigor210, cBioPortal, GeneMANIA, and TIMER were processed and visualized to explore the association between FAM72A and the prognosis, stemness phenotype, mutational burden, immune cell infiltration, and drug sensitivity in HCC patients. Potential pathways were also revealed. Furthermore, we experimentally verified the results in vivo and in vitro using immunohistochemistry, western blotting, and CCK-8 assays.

Results

First, FAM72A mRNA expression was significantly upregulated in HCC. High FAM72A expression was independently associated with a poor prognosis. Experimental validation confirmed that FAM72A was remarkably overexpressed in HCC patients and mice. Moreover, FAM72A knockdown suppressed HCC cell proliferation. In addition, the frequency of TP53 mutations was significantly higher in the high FAM72A expression group. Subsequently, the enrichment analysis revealed that FAM72A was closely related to immune processes and mTOR pathways. Silencing FAM72A increased the expression levels of mTOR in HCC cell lines. The FAM72A-mTOR pathway was strongly associated with a poor prognosis for patients with HCC. Patients with high FAM72A expression levels might be more resistant to sorafenib. Furthermore, the expression of FAM72A and mTOR was significantly associated with the abundance of some tumor-infiltrating immune cells, especially CD4+ T cells. Finally, patients with high levels of FAM72A and mTOR were more sensitive to immunotherapy.

Conclusions

FAM72A, a member of the FAM72 family, might be a prognostic and immunotherapeutic target for HCC patients.

Fuzheng Huayu Recipe and its active compounds inhibited HBeAg production by promoting TOMM34 gene expression in HBV-infected hepatocytes

Background and aim: Fuzheng Huayu Recipe (FZHY) is a Chinese patent medicine (approval No. Z20020074) included in the national medical insurance catalogue, which is mainly used for anti-hepatic fibrosis treatment of hepatitis B virus (HBV) induced liver fibrosis and liver cirrhosis. In clinical practice, we discovered that FZHY might also have a direct anti-HBV effect on inhibiting HBeAg production, but the mechanism underlying was unclear. This study aimed to clarify the molecular mechanism of the inhibition effect of FZHY on HBeAg production.

Methods: The decrease degree of serum HBeAg titer in FZHY + entecavir (ETV) group patients were analyzed through clinical data. C57BL/6N-Tg (1.28HBV)/Vst HBV transgenic mice were used for in vivo experiments. HepG2. 2.15 cells (wild-type HBV replication cells) were used for in vitro experiments.

Results: The clinical study results showed that the decrease degree of serum HBeAg titer in FZHY+ETV group was significantly higher than that in ETV group after 48 weeks treatment. In vivo experiments results showed that FZHY could significantly reduce the serum HBeAg titer in HBV transgenic mice, and promote HBeAg seroconversion. In vitro experiments results showed that FZHY could reduce HBeAg titer dependently, but it did not significantly inhibit the expression of HBsAg and HBV-DNA. Further cell experiments in vitro discovered that TOMM34 might be the key target for FZHY to inhibit HBeAg production. The subsequent pharmacological screening experiment of 20 active compounds in FZHY showed that quercetin, baicalin and cordycepin could promote the expression of TOMM34 gene and reduce the production of HBeAg.

Conclusion: In conclusion, FZHY and its active compounds quercetin, baicalin and cordycepin could inhibit HBeAg production by promoting the expression of TOMM34 gene in HBV-infected hepatocytes.

Identification and analysis of mitochondria-related key genes of heart failure

Mitochondria-induced cell death is a vital mechanism of heart failure (HF). Thus, identification of mitochondria-related genes (Mito-RGs) based on transcriptome sequencing data of HF might provide novel diagnostic markers and therapeutic targets for HF. First, bioinformatics analysis was conducted on the GSE57338, GSE76701, GSE136547, and GSE77399 datasets in the Gene Expression Omnibus. Next, we analyzed HF-Mito differentially expressed genes (DEGs) using the protein-protein interaction (PPI) network for obtaining critical genes and exploring their functions. Subsequently, immune cell scores of the HF and normal groups were compared. The potential alteration mechanisms of the key genes were investigated by constructing a competing endogenous RNA network. Finally, we predicted potential therapeutic agents and validated the expression levels of the key genes. Twenty-three HF-Mito DEGs were acquired in the GSE57338 dataset, and the PPI network obtained four key genes, including IFIT3, XAF1, RSAD2, and MX1. According to gene set enrichment analysis, the key genes showed high enrichment in myogenesis and hypoxia. Immune cell analysis demonstrated that aDCs, B cells, and 20 other immune cell types varied between the HF and normal groups. Moreover, we observed that H19 might affect the expression of IFIT3, AXF1, and RSAD2. PCGEM1 might regulate RSAD2 expression. A total of 515 potential therapeutic drugs targeting the key genes, such as tretinoin, silicon dioxide, and bisphenol A, were acquired. Finally, IFIT3, RSAD2, and MX1 expression increased in HF samples compared with normal samples in the GSE76701 dataset, conforming to the GSE57338 dataset analysis. This work screened four key genes, namely, IFIT3, XAF1, RSAD2, and MX1, which can be further explored in subsequent studies for their specific molecular mechanisms in HF.

RNA methylation regulators contribute to poor prognosis of hepatocellular carcinoma associated with the suppression of bile acid metabolism: a multi-omics analysis

RNA methylation has been known to promote the initiation and progression of many types of cancer, including hepatocellular carcinoma (HCC). To fully understand the importance of this post-transcriptional modification in HCC, a thorough investigation that combines different patterns of RNA methylation is urgently needed. In this study, we investigated the regulators of the three most common types of RNA methylation: m6A, N1-methyladenosine (m1A) and 5-methylcytosine (m5C). Based on the genomic and proteomic data, we constructed a classifier consisting of seven RNA methylation regulators. This classifier performed well and robustly predicted the prognosis of HCC patients. By analysis using this classifier, we found that the primary bile acid biosynthesis pathway was mostly downregulated in high-risk HCC patients. Furthermore, we found that the gene expression patterns regulated by several bile acids were similar to those regulated by some well-defined anti-tumor compounds, indicating that bile acid metabolism plays a crucial role in the progression of HCC, and the related metabolites can be used as the potential agents for HCC treatments. Moreover, our study revealed a crosstalk between RNA methylation and bile acid regulators, demonstrating a novel mechanism of the downregulation of bile acid metabolism in HCC and providing new insights into how RNA methylation regulators affect the oncogenesis of HCC.

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

Background

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

Materials and Methods

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

Results

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

Conclusions

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

Diagnostic Value, Prognostic Value, and Immune Infiltration of LOX Family Members in Liver Cancer: Bioinformatic Analysis

Background

Liver cancer (LC) is well known for its prevalence as well as its poor prognosis. The aberrant expression of lysyl oxidase (LOX) family is associated with liver cancer, but their function and prognostic value in LC remain largely unclear. This study aimed to explore the function and prognostic value of LOX family in LC through bioinformatics analysis and meta-analysis.

Results

The expression levels of all LOX family members were significantly increased in LC. Area under the receiver operating characteristic curve (AUC) of LOXL2 was 0.946 with positive predictive value (PPV) of 0.994. LOX and LOXL3 were correlated with worse prognosis. Meta-analysis also validated effect of LOX on prognosis. Nomogram of these two genes and other predictors was also plotted. There was insufficient data from original studies to conduct meta-analysis on LOXL3. The functions of LOX family members in LC were mostly involved in extracellular and functions and structures. The expressions of LOX family members strongly correlated with various immune infiltrating cells and immunomodulators in LC.

Conclusions

For LC patients, LOXL2 may be a potential diagnostic biomarker, while LOX and LOXL3 have potential prognostic and therapeutic values. Positive correlation between LOX family and infiltration of various immune cells and immunomodulators suggests the need for exploration of their roles in the tumor microenvironment and for potential immunotherapeutic to target LOX family proteins.

Hepatotoxicity of cantharidin is associated with the altered bile acid metabolism

Cantharidin (CTD) is an effective antitumor agent. However, it exhibits significant hepatotoxicity, the mechanism of which remains unclear. In this study, biochemical and histopathological analyses complemented with ultra-high-performance liquid chromatography-tandem mass spectrometry (UHPLC-MS/MS)-based targeted metabolomic analysis of bile acids (BAs) were employed to investigate CTD-induced hepatotoxicity in rats. Sixteen male and female Sprague-Dawley rats were randomly divided into two groups: control and CTD (1.0 mg/kg) groups. Serum and liver samples were collected after 28 days of intervention. Biochemical, histopathological, and BA metabolomic analyses were performed for all samples. Further, the key biomarkers of CTD-induced hepatotoxicity were identified via multivariate and metabolic pathway analyses. In addition, metabolite-gene-enzyme network and Kyoto Encyclopedia of Genes and Genomes pathway analyses were used to identify the signaling pathways related to CTD-induced hepatotoxicity. The results revealed significantly increased levels of biochemical indices (alanine aminotransferase, aspartate aminotransferase, and total bile acid). Histopathological analysis revealed that the hepatocytes were damaged. Further, 20 endogenous BAs were quantitated via UHPLC-MS/MS, and multivariate and metabolic pathway analyses of BAs revealed that hyocholic acid, cholic acid, and chenodeoxycholic acid were the key biomarkers of CTD-induced hepatotoxicity. Meanwhile, primary and secondary BA biosynthesis and taurine and hypotaurine metabolism were found to be associated with the mechanism by which CTD induced hepatotoxicity in rats. This study provides useful insights for research on the mechanism of CTD-induced hepatotoxicity.

Identification of a Chemotherapeutic Lead Molecule for the Potential Disruption of the FAM72A-UNG2 Interaction to Interfere with Genome Stability, Centromere Formation, and Genome Editing

Family with sequence similarity 72 A (FAM72A) is a pivotal mitosis-promoting factor that is highly expressed in various types of cancer. FAM72A interacts with the uracil-DNA glycosylase UNG2 to prevent mutagenesis by eliminating uracil from DNA molecules through cleaving the N-glycosylic bond and initiating the base excision repair pathway, thus maintaining genome integrity. In the present study, we determined a specific FAM72A-UNG2 heterodimer protein interaction using molecular docking and dynamics. In addition, through in silico screening, we identified withaferin B as a molecule that can specifically prevent the FAM72A-UNG2 interaction by blocking its cell signaling pathways. Our results provide an excellent basis for possible therapeutic approaches in the clinical treatment of cancer.

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

INTRODUCTION

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

MATERIALS AND METHODS

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

RESULTS

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

CONCLUSIONS

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