ceRNA network development and tumor-infiltrating immune cell analysis in hepatocellular carcinoma

Cuproptosis-Related Genes MTF1 and LIPT1 as Novel Prognostic Biomarker in Acute Myeloid Leukemia

Acute myeloid leukemia (AML) is a life-threatening hematologic malignant disease with high morbidity and mortality in both adults and children. Cuproptosis, a novel mode of cell death, plays an important role in tumor development, but the functional mechanisms of cuproptosis-related genes (CRGs) in AML are unclear. The differential expression of CRGs between tumors such as AML and normal tissues in UCSC XENA, TCGA and GTEx was verified using R (version: 3.6.3). Lasso regression, Cox regression and Nomogram were used to screen for prognostic biomarkers of AML and to construct corresponding prognostic models. Kaplan-Meier analysis, ROC analysis, clinical correlation analysis, immune infiltration analysis and enrichment analysis were used to further investigate the correlation and functional mechanisms of CRGs with AML. The ceRNA regulatory network was used to identify the mRNA-miRNA-lncRNA regulatory axis. Cuproptosis-related genes LIPT1, MTF1, GLS and CDKN2A were highly expressed in AML, while FDX1, LIAS, DLD, DLAT, PDHA1, SLC31A1 and ATP7B were lowly expressed in AML. Lasso regression, Cox regression, Nomogram and calibration curve finally identified MTF1 and LIPT1 as two novel prognostic biomarkers of AML and constructed the corresponding prognostic models. In addition, all 12 CRGs had predictive power for AML, with MTF1, LIAS, SLC31A1 and CDKN2A showing more reliable results. Further analysis showed that ATP7B was closely associated with mutation types such as FLT3, NPM1, RAS and IDH1 R140 in AML, while the expression of MTF1, LIAS and ATP7B in AML was closely associated with immune infiltration. Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Set Enrichment Analysis (GSEA) revealed that biological functions such as metal ion transmembrane transporter activity, haptoglobin binding and oxygen carrier activity, pathways such as interferon alpha response, coagulation, UV response DN, apoptosis, hypoxia and heme metabolism all play a role in the development of AML. The ceRNA regulatory network revealed that 6 lncRNAs such as MALAT1, interfere with MTF1 expression through 6 miRNAs such as hsa-miR-32-5p, which in turn affect the development and progression of AML. In addition, APTO-253 has the potential to become an AML-targeted drug. The cuproptosis-related genes MTF1 and LIPT1 can be used as prognostic biomarkers in AML. A total of six lncRNAs, including MALAT1, are involved in the expression and regulation of MTF1 in AML through six miRNAs such as hsa-miR-32-5p.

Screening immune-related blood biomarkers for DKD-related HCC using machine learning

Background

Diabetes mellitus is a significant health problem worldwide, often leading to diabetic kidney disease (DKD), which may also influence the occurrence of hepatocellular carcinoma (HCC). However, the relationship and diagnostic biomarkers between DKD and HCC are unclear.

Methods

Using public database data, we screened DKD secretory RNAs and HCC essential genes by limma and WGCNA. Potential mechanisms, drugs, and biomarkers for DKD-associated HCC were identified using PPI, functional enrichment, cMAP, and machine learning algorithms, and a diagnostic nomogram was constructed. Then, ROC, calibration, and decision curves were used to evaluate the diagnostic performance of the nomograms. In addition, immune cell infiltration in HCC was explored using CIBERSORT. Finally, the detectability of critical genes in blood was verified by qPCR.

Results

104 DEGs associated with HCC using WGCNA were identified. 101 DEGs from DKD were predicated on secreting into the bloodstream with Exorbase datasets. PPI analysis identified three critical modules considered causative genes for DKD-associated HCC, primarily involved in inflammation and immune regulation. Using lasso and RM, four hub genes associated with DKD-associated HCC were identified, and a diagnostic nomogram confirmed by DCA curves was established. The results of immune cell infiltration showed immune dysregulation in HCC, which was associated with the expression of four essential genes. PLVAP was validated by qPCR as a possible blood-based diagnostic marker for DKD-related HCC.

Conclusion

We revealed the inflammatory immune pathways of DKD-related HCC and developed a diagnostic nomogram for HCC based on PLVAP, C7, COL15A1, and MS4A6A. We confirmed with qPCR that PLVAP can be used as a blood marker to assess the risk of HCC in DKD patients.

Developing four cuproptosis-related lncRNAs signature to predict prognosis and immune activity in ovarian cancer

BACKGROUND

There has been a recent discovery of a new type of cell death produced by copper-iron ions, called Cuproptosis (copper death). The purpose of this study was to identify LncRNA signatures associated with Cuproptosis in ovarian cancer that could be used as prognostic indicators.

METHODS

RNA sequencing (RNA-seq) profiles with clinicopathological data from TCGA database were used to select prognostic CRLs and then constructed prognostic risk model using multivariate regression analysis and LASSO algorithms. An independent dataset from GEO database was used to validate the prognostic performance. Combined with clinical factors, we further constructed a prognostic nomogram. In addition, tumor immune microenvironment, somatic mutation and drug sensitivity were analyzed using ssGSEA, GSVA, ESTIMATE and CIBERSORT algorithms.

RESULT

A total of 129 CRLs were selected whose expression levels were significantly related to expression levels of 10 cuproptosis-related genes. The univariate Cox regression analysis showed that 12 CRLs were associated with overall survival (OS). Using LASSO algorithms and multivariate regression analysis, we constructed a four-CRLs prognostic signature in the training dataset. Patients in the training dataset could be classified into high- or low-risk subgroups with significantly different OS (log-rank p < 0.001). The prognostic performance was confirmed in TCGA-OC cohort (log-rank p < 0.001) and an independent GEO cohort (log-rank p = 0.023). Multivariate cox regression analysis proved the four-CRLs signature was an independent prognostic factor for OC. Additionally, different risk subtypes showed significantly different levels of immune cells, signal pathways, and drug response.

CONCLUSION

We established a prognostic signature based on cuproptosis-related lncRNAs for OC patients, which will be of great value in predicting the prognosis patients and may provide a new perspective for research and individualized treatment.

Delineation of a SMARCA4-specific competing endogenous RNA network and its function in hepatocellular carcinoma

BACKGROUND

Hepatocellular carcinoma (HCC) is a common malignancy worldwide, and the mortality rate continues to rise each year. SMARCA4 expression has been associated with poor prognosis in various types of cancer; however, the specific mechanism of action of SMARCA4 in HCC needs to be fully elucidated.

AIM

To explore the specific mechanism of action of SMARCA4 in HCC.

METHODS

Herein, the expression level of SMARCA4 as well as its association with HCC prognosis were evaluated using transcriptome profiling and clinical data of 18 different types of cancer collected from The Cancer Genome Atlas database. Furthermore, SMARCA4-high and -low groups were identified. Thereafter, gene ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analyses were performed to identify the function of SMARCA4, followed by construction of a SMARCA4-specific competing endogenous RNA (ceRNA) network using starBase database. The role of SMARCA4 in immunotherapy and its association with immune cells were assessed using correlation analysis.

RESULTS

It was observed that SMARCA4 was overexpressed and negatively correlated with prognosis in HCC. Further, SMARCA4 expression was positively associated with tumor mutational burden, microsatellite stability, and immunotherapy efficacy. The SNHG3/THUMP3-AS1-miR-139-5p-SMARCA4 ceRNA network was established and could be assumed to serve as a stimulatory mechanism in HCC.

CONCLUSION

The findings of this study demonstrated that SMARCA4 plays a significant role in progression and immune infiltration in HCC. Moreover, a ceRNA network was detected, which was found to be correlated with poor prognosis in HCC. The findings of this study could contribute towards the identification of predictive markers for immunotherapy and a novel mechanism of action for HCC treatment.

Identification and Validation of Cuproptosis-Related Prognostic Signature and Associated Regulatory Axis in Uterine Corpus Endometrial Carcinoma

Background: Uterine corpus endometrial carcinoma (UCEC) is a common gynecological malignancy globally with high recurrence and mortality rates. Cuproptosis is a new type of programmed cell death involved in tumor cell proliferation and growth, angiogenesis, and metastasis.

Methods: The difference in cuproptosis-related genes (CRGs) between UCEC tissues and normal tissues deposited in The Cancer Genome Atlas database was calculated using the “limma” R package. LASSO Cox regression analysis was conducted to construct a prognostic cuproptosis–related signature. Kaplan–Meier analysis was conducted to compare the survival of UCEC patients. A ceRNA network was constructed to identify the lncRNA–miRNA–mRNA regulatory axis. Quantitative reverse transcription-polymerase chain reaction (qRT-PCR) was performed to verify CRG expression in UCEC.

Results: The expression of FDX1, LIAS, DLAT, and CDKN2A were upregulated, whereas the expression of LIPT1, DLD, PDHB, MTF1, and GLS were downregulated in UCEC versus normal tissues. The genetic mutation landscape of CRGs in UCEC was also summarized. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analyses revealed that these CRGs were enriched in the tricarboxylic acid (TCA) cycle, glycolysis, and HIF-1 signaling pathway. LASSO Cox regression analysis was performed and identified a cuproptosis-related prognostic signature including these three prognostic biomarkers (CDKN2A, GLS, and LIPT1). UCEC patients with high risk scores had a poor prognosis with an area under the curve of 0.782 and 0.764 on 3- and 5-year receiver operating characteristic curves. Further analysis demonstrated a significant correlation between CDKN2A and pTNM stage, tumor grade, immune cell infiltration, drug sensitivity, tumor mutational burden (TMB) score, and microsatellite instable (MSI) score. The data validation of qRT-PCR further demonstrated the upregulation of CDKN2A and the downregulation of LIPT1 and GLS in UCEC versus normal tissues. The ceRNA network also identified lncRNA XIST/miR-125a-5p/CDKN2A regulatory axis for UCEC.

Conclusion: The current study identified a cuproptosis-related prognostic signature including these three prognostic biomarkers (CDKN2A, GLS, and LIPT1) for UCEC. The ceRNA network also identified that lncRNA XIST/miR-125a-5p/CDKN2A regulatory axis may be involved in the progression of UCEC. Further in vivo and in vitro studies should be conducted to verify these results.

Mechanism of a Herbal Formula Associated with Prognosis and Immune Infiltration in LIHC: Transcriptomics Analysis and Molecular Dynamics Simulations

Background

The aim of this study is to explore the interactions between effective monomers of herbal formulas and their therapeutic targets using systems biology approaches which may be a promising approach to unraveling their underlying mechanisms. Shentao Ruangan decoction (STRGD), which has been experimentally, clinically demonstrated to be effective in treating liver hepatocellular carcinoma (LIHC), was selected.

Methods

Bioactive ingredients and drug targets of STRGD were retrieved from the traditional Chinese medicine systems pharmacology database and analysis platform and BATMAN-TCM databases. LIHC-related differentially expressed genes (DEGs) and key modules were identified by a weighted gene coexpression network analysis using The Cancer Genome Atlas data. The Kaplan–Meier analysis was used to investigate the relationship between STRGD tumor targets and patients survival. The CIBERSORT deconvolution algorithm was used to analyze the correlation between STRGD tumor targets and infiltrating immune cells. Enrichment analysis was used to analyze biological functions. Interactions between STRGD compounds and LIHC-immune-related genes were investigated using molecular docking and MDS.

Results

We identified 24 STRGD tumor targets, which were found to be correlated with survival and the level of immune cell infiltration in LIHC patients. Immune infiltration, gene set enrichment, and Kyoto Encyclopedia of Genes and Genomes analyses highlighted the roles of T and B cell subsets, which were both related to activator protein 1 (AP1), in STRGD action. Docking studies and HPLC indicated that tanshinone IIA is the main compound of STRGD in LIHC treatment, and MDS showed that the potential LIHC-immune-related targets 1FOS and 1JUN firmly bind to tanshinone IIA.

Conclusions

The mechanisms of STRGD in improving the immune and survival status of LIHC patients include interactions between STRGD compounds and LIHC-immune-related targets. The findings of this study can guide research studies on the potential usefulness of tanshinone IIA in the development of drugs targeting 1JUN and 1FOS for the treatment of LIHC.

Single-Cell Sequencing Identifies the Heterogeneity of CD8+ T Cells and Novel Biomarker Genes in Hepatocellular Carcinoma

CD8+ T cells are required for the establishment of antitumor immunity, and their substantial infiltration is associated with a good prognosis. However, CD8+ T cell subsets in the tumor microenvironment may play distinct roles in tumor progression, prognosis, and immunotherapy. In this study, we used the scRNA-seq data of hepatocellular carcinoma (HCC) to reveal the heterogeneity of different CD8+ T cell subsets. The scRNA-seq data set GSE149614 was obtained from the GEO database, and the transcriptome and sample phenotypic data of TCGA-LIHC were obtained from the TCGA database. CD8+ T cell subtypes and metabolic gene sets were obtained from published reports. The data processing and analysis of CD8+ T cell groups was performed by R language. The PPI network was constructed to obtain the hub genes, and the KM survival curve of the hub genes was further plotted to determine the hub genes with differences in survival. CD8+ T cells in HCC were divided into 7 subsets, and the cytotoxic CD8 T cells 4 subset showed considerable differences between the TP53-mutant and nonmutant groups, as well as between different degrees of cirrhosis, HCC grades, stages, ages, and body weights. Cytotoxic CD8 T cells 4 differential genes were analyzed by TCGA-LIHC data and single-cell sequencing data set. 10 hub genes were found: FGA, ApoA1, ApoH, AHSG, FGB, HP, TTR, TF, HPX, and APOC3. Different subsets of CD8+ T cells were found to contribute to heterogeneous prognosis and pathway activity in HCC. Alterations in the cytotoxic and immune checkpoint gene expression during CD8+ T cell differentiation were also identified. We found that cytotoxic CD8 T cells 4 is closely associated with survival and prognosis of HCC and identified four differential genes that can be used as biological markers for survival, prognosis, and clinically relevant characteristics of HCC. Results of this study could help finding targets for immunotherapy of HCC and aid in the accelerated development of immunotherapy for HCC.

Integrated clinical characteristics and omics analysis identifies a ferroptosis and iron-metabolism-related lncRNA signature for predicting prognosis and therapeutic responses in ovarian cancer

Background

Ferroptosis and iron-metabolism are regulated by Long non-coding RNAs (lncRNAs) in ovarian cancer (OC). Therefore, a comprehensive analysis of ferroptosis and iron-metabolism related lncRNAs (FIRLs) in OC is crucial for proposing therapeutic strategies and survival prediction.

Methods

In multi-omics data from OC patients, FIRLs were identified by calculating Pearson correlation coefficients with ferroptosis and iron-metabolism related genes (FIRGs). Cox-Lasso regression analysis was performed on the FIRLs to screen further the lncRNAs participating in FIRLs signature. In addition, all patients were divided into two robust risk subtypes using the FIRLs signature. Receiver operator characteristic (ROC) curve, Kaplan–Meier analysis, decision curve analysis (DCA), Cox regression analysis and calibration curve were used to confirm the clinical benefits of FIRLs signature. Meanwhile, two nomograms were constructed to facilitate clinical application. Moreover, the potential biological functions of the signature were investigated by genes function annotation. Finally, immune microenvironment, chemotherapeutic sensitivity, and the response of PARP inhibitors were compared in different risk groups using diversiform bioinformatics algorithms.

Results

The raw data were randomized into a training set (n = 264) and a testing set (n = 110). According to Pearson coefficients between FIRGs and lncRNAs, 1075 FIRLs were screened for univariate Cox regression analysis, and then LASSO regression analysis was used to construct 8-FIRLs signature. It is worth mentioning that a variety of analytical methods indicated excellent predictive performance for overall survival (OS) of FIRLs signature (p < 0.05). The multivariate Cox regression analysis showed that FIRLs signature was an independent prognostic factor for OS (p < 0.05). Moreover, significant differences in the abundance of immune cells, immune-related pathways, and drug response were excavated in different risk subtypes (p < 0.05).

Conclusion

The FIRLs signature can independently predict overall survival and therapeutic effect in OC patients.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13048-022-00944-y.

Bioinformatics Analysis Using ATAC-seq and RNA-seq for the Identification of 15 Gene Signatures Associated With the Prediction of Prognosis in Hepatocellular Carcinoma

Background

Gene expression (RNA-seq) and overall survival (OS) in TCGA were combined using chromosome accessibility (ATAC-seq) to search for key molecules affecting liver cancer prognosis.

Methods

We used the assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) to analyse chromatin accessibility in the promoter regions of whole genes in liver hepatocellular carcinoma (LIHC) and then screened differentially expressed genes (DEGs) at the mRNA level by transcriptome sequencing technology (RNA-seq). We obtained genes significantly associated with overall survival (OS) by a one-way Cox analysis. The three were screened by taking intersection and further using a Kaplan–Meier (KM) for validation. A prognostic model was constructed using the obtained genes by LASSO regression analysis.The expression of these genes in hepatocellular carcinomas was then analysed. The protein expression of these genes was verified using the Human Protein Atlas(HPA) online datasets and immunohistochemistry.

Results

ATAC-seq, RNA-seq and survival analysis, combined with a LASSO prediction model, identified signatures of 15 genes (PRDX6, GCLM, HTATIP2, SEMA3F, UCK2, NOL10, KIF18A, RAP2A, BOD1, GDI2, ZIC2, GTF3C6 SLC1A5, ERI3 and SAC3D1), all of which were highly expressed in hepatocellular carcinoma. The LASSO prognostic model showed that this risk score had high predictive accuracy for the survival prognosis at 1, 3 and 5 years. A KM curve analysis showed that high expression of all 15 gene signatures was significantly associated with a poor prognosis in LIHC patients. HPA analysis of protein expression showed that PRDX6, GCLM, HTATIP2, NOL10, KIF18A, RAP2A and GDI2 were highly expressed in the hepatocellular carcinoma tissues compared with normal control tissues.

Conclusions

PRDX6, GCLM, HTATIP2, SEMA3F, UCK2, NOL10, KIF18A, RAP2A, BOD1, GDI2, ZIC2, GTF3C6, SLC1A5, ERI3 and SAC3D1 may affect the prognosis of LIHC.