A five-mRNA signature associated with post-translational modifications can better predict recurrence and survival in cervical cancer

Diminazene aceturate-induced cytotoxicity is associated with the deregulation of cell cycle signaling and downregulation of oncogenes Furin, c-MYC, and FOXM1 in human cervical carcinoma Hela cells

Diminazene aceturate (DIZE) is an FDA-listed small molecule known for the treatment of African sleeping sickness. In vivo studies showed that DIZE may be beneficial for a range of human ailments. However, there is very limited information on the effects of DIZE on human cancer cells. The current study aimed to investigate the cytotoxic responses of DIZE, using the human carcinoma Hela cell line. WST-1 cell proliferation assay showed that DIZE inhibited the viability of Hela cells in a dose-dependent manner and the observed response was associated with the downregulation of Ki67 and PCNA cell proliferation markers. DIZE-treated cells stained with acridine orange-ethidium and JC-10 dye revealed cell death and loss of mitochondrial membrane potential (Ψm), compared with DMSO (vehicle) control, respectively. Cellular immunofluorescence staining of DIZE-treated cells showed upregulation of caspase 3 activities. DIZE-treated cells showed downregulation of mRNA for G1/S genes CCNA2 and CDC25A, S-phase genes MCM3 and PLK4, and G2/S phase transition/mitosis genes Aurka and PLK1. These effects were associated with decreased mRNA expression of Furin, c-Myc, and FOXM1 oncogenes. These results suggested that DIZE may be considered for its effects on other cancer types. To the best of our knowledge, this is the first study to evaluate the effect of DIZE on human cervical cancer cells.

Revealing the immune cell subtype reconstitution profile in patients from the CLARITY study using deconvolution algorithms after cladribine tablets treatment

Immune Cell Deconvolution methods utilizing gene expression profiling to quantify immune cells in tissues and blood are an appealing alternative to flow cytometry. Our objective was to investigate the applicability of deconvolution approaches in clinical trial settings to better investigate the mode of action of drugs for autoimmune diseases. Popular deconvolution methods CIBERSORT and xCell were validated using gene expression from the publicly available GSE93777 dataset that has comprehensive matching flow cytometry. As shown in the online tool, ~ 50% of signatures show strong correlation (r > 0.5) with the remainder showing moderate correlation, or in a few cases, no correlation. Deconvolution methods were then applied to gene expression data from the phase III CLARITY study (NCT00213135) to evaluate the immune cell profile of relapsing multiple sclerosis patients treated with cladribine tablets. At 96 weeks after treatment, deconvolution scores showed the following changes vs placebo: naïve, mature, memory CD4⁺ and CD8⁺ T cells, non-class switched, and class switched memory B cells and plasmablasts were significantly reduced, naïve B cells and M2 macrophages were more abundant. Results confirm previously described changes in immune cell composition following cladribine tablets treatment and reveal immune homeostasis of pro- vs anti-inflammatory immune cell subtypes, potentially supporting long-term efficacy.

Immune landscape and risk prediction based on pyroptosis-related molecular subtypes in triple-negative breast cancer

The survival outcome of triple-negative breast cancer (TNBC) remains poor, with difficulties still existing in prognosis assessment and patient stratification. Pyroptosis, a newly discovered form of programmed cell death, is involved in cancer pathogenesis and progression. The role of pyroptosis in the tumor microenvironment (TME) of TNBC has not been fully elucidated. In this study, we disclosed global alterations in 58 pyroptosis-related genes at somatic mutation and transcriptional levels in TNBC samples collected from The Cancer Genome Atlas and Gene Expression Omnibus databases. Based on the expression patterns of genes related to pyroptosis, we identified two molecular subtypes that harbored different TME characteristics and survival outcomes. Then, based on differentially expressed genes between two subtypes, we established a 12-gene score with robust efficacy in predicting short- and long-term overall survival of TNBC. Patients at low risk exhibited a significantly better prognosis, more antitumor immune cell infiltration, and higher expression of immune checkpoints including PD-1, PD-L1, CTLA-4, and LAG3. The comprehensive analysis of the immune landscape in TNBC indicated that alterations in pyroptosis-related genes were closely related to the formation of the immune microenvironment and the intensity of the anticancer response. The 12-gene score provided new information on the risk stratification and immunotherapy strategy for highly heterogeneous patients with TNBC.

A Novel Prognostic Risk Model for Cervical Cancer Based on Immune Checkpoint HLA-G-Driven Differentially Expressed Genes

Human leukocyte antigen G (HLA-G) is a potential checkpoint molecule that plays a key role in cervical carcinogenesis. The purpose of this study was to construct and validate a prognostic risk model to predict the overall survival (OS) of cervical cancer patients, providing a reference for individualized clinical treatment that may lead to better clinical outcomes. HLA-G-driven differentially expressed genes (DEGs) were obtained from two cervical carcinoma cell lines, namely, SiHa and HeLa, with stable overexpression of HLA-G by RNA sequencing (RNA-seq). The biological functions of these HLA-G-driven DEGs were analysed by GO enrichment and KEGG pathway using the “clusterProfiler” package. The protein-protein interactions (PPIs) were assessed using the STRING database. The prognostic relevance of each DEG was evaluated by univariate Cox regression using the TCGA-CESC dataset. After the TCGA-CESC cohort was randomly divided into training set and testing set, and a prognostic risk model was constructed by LASSO and stepwise multivariate Cox regression analysis in training set and validated in testing set or in different types of cervical cancer set. The predictive ability of the prognostic risk model or nomogram was evaluated by a series of bioinformatics methods. A total of 1108 candidate HLA-G-driven DEGs, including 391 upregulated and 717 downregulated genes, were obtained and were enriched mostly in the ErbB pathway, steroid biosynthesis, and MAPK pathway. Then, an HLA-G-driven DEG signature consisting of the eight most important prognostic genes CD46, LGALS9, PGM1, SPRY4, CACNB3, PLIN2, MSMO1, and DAGLB was identified as a key predictor of cervical cancer. Multivariate Cox regression analysis showed that this signature is an independent risk factor for the overall survival of CESC patients. Kaplan-Meier survival analysis showed that the 5-year overall survival rate is 23.0% and 84.6% for the high-risk and low-risk patients, respectively (P<0.001). The receiver operating characteristic (ROC) curve of this prognostic model with an area under the curve (AUC) was 0.896 for 5 years, which was better than that of other clinical traits. This prognostic risk model was also successfully validated in different subtypes of cervical cancer, including the keratinizing squamous cell carcinoma, non-keratinizing squamous cell carcinoma, squamous cell neoplasms, non-squamous cell neoplasms set. Single-sample gene set enrichment (ssGSEA) algorithm and Tumor Immune Dysfunction and Exclusion (TIDE) analysis confirmed that this signature influence tumour microenvironment and immune checkpoint blockade. A nomogram that integrated risk score, age, clinical stage, histological grade, and pathological type was then built to predict the overall survival of CESC patients and evaluated by calibration curves, AUC, concordance index (C-index) and decision curve analysis (DCA). To summarize, we developed and validated a novel prognostic risk model for cervical cancer based on HLA-G-driven DEGs, and the prognostic signature showed great ability in predicting the overall survival of patients with cervical cancer.

Identifying a cervical cancer survival signature based on mRNA expression and genome-wide copy number variations

Cervical cancer mortality is the second highest in gynecological cancers. This study developed a new model based on copy number variation data and mRNA data for overall survival prediction of cervical cancer. Differentially expressed genes from The Cancer Genome Atlas dataset detected by univariate Cox regression analysis were further simplified to six by least absolute shrinkage and selection operator (Lasso) and stepwise Akaike information criterion (stepAIC). The study developed a six-gene signature, which was further verified in independent dataset. Association between immune infiltration and risk score was investigated by immune score. The relation between the signature and functional pathways was examined by gene set enrichment analysis. Ninety-nine differentially expressed genes were detected, and C11orf80, FOXP3, GSN, HCCS, PGAM5, and RIBC2 were identified as key genes to construct a six-gene signature. The prognostic signature showed a significant correlation with overall survival (hazard ratio, HR = 3.45, 95% confidence interval (CI) = 2.08-5.72, p < 0.00001). Immune score showed a negative correlation with the risk score calculated by the signature (p < 0.05). Four immune-related pathways were closely associated with risk score (p < 0.0001). The six-gene prognostic signature was an effective tool to predict overall survival of cervical cancer. In conclusion, the newly identified six genes may be considered as new drug targets for cervical cancer treatment.

A novel lncRNA-mRNA-miRNA signature predicts recurrence and disease-free survival in cervical cancer

Cervical cancer (CC) patients have a poor prognosis due to the high recurrence rate. However, there are still no effective molecular signatures to predict the recurrence and survival rates for CC patients. Here, we aimed to identify a novel signature based on three types of RNAs [messenger RNA (mRNAs), microRNA (miRNAs), and long non-coding RNAs (lncRNAs)]. A total of 763 differentially expressed mRNAs (DEMs), 46 lncRNAs (DELs), and 22 miRNAs (DEMis) were identified between recurrent and non-recurrent CC patients using the datasets collected from the Gene Expression Omnibus (GSE44001; training) and The Cancer Genome Atlas (RNA- and miRNA-sequencing; testing) databases. A competing endogenous RNA network was constructed based on 23 DELs, 15 DEMis, and 426 DEMs, in which 15 DELs, 13 DEMis, and 390 DEMs were significantly associated with disease-free survival (DFS). A prognostic signature, containing two DELs (CD27-AS1, LINC00683), three DEMis (hsa-miR-146b, hsa-miR-1238, hsa-miR-4648), and seven DEMs (ARMC7, ATRX, FBLN5, GHR, MYLIP, OXCT1, RAB39A), was developed after LASSO analysis. The built risk score could effectively separate the recurrence rate and DFS of patients in the high- and low-risk groups. The accuracy of this risk score model for DFS prediction was better than that of the FIGO (International Federation of Gynecology and Obstetrics) staging (the area under receiver operating characteristic curve: training, 0.954 vs 0.501; testing, 0.882 vs 0.656; and C-index: training, 0.855 vs 0.539; testing, 0.711 vs 0.508). In conclusion, the high predictive accuracy of our signature for DFS indicated its potential clinical application value for CC patients.

Development of a Predictive Immune-Related Gene Signature Associated With Hepatocellular Carcinoma Patient Prognosis

Background:

Hepatocellular carcinoma (HCC) remains the third leader cancer-associated cause of death globally, but the etiological basis for this complex disease remains poorly clarified. The present study was thus conceptualized to define a prognostic immune-related gene (IRG) signature capable of predicting immunotherapy responsiveness and overall survival (OS) in patients with HCC.

Methods:

Five differentially expressed IRG associated with HCC were established the immune-related risk model through univariate Cox regression and least absolute shrinkage and selection operator (LASSO) regression analyses. Patients were separated at random into training and testing cohorts, after which the association between the identified IRG signature and OS was evaluated using the “survival” R package. In addition, maftools was leveraged to assess mutational data, with tumor mutation burden (TMB) scores being calculated as follows: (total mutations/total bases) × 10⁶. Immune-related risk term abundance was quantified via “ssGSEA” algorithm using the “gsva” R package.

Results:

HCC patients were successfully stratified into low-risk and high-risk groups based upon a signature composed of 5 differentially expressed IRGs, with overall survival being significantly different between these 2 groups in training cohort, testing cohort and overall patient cohort (P = 1.745e-06, P = 1.888e-02, P = 4.281e-07). No association was observed between TMB and this IRG risk score in the overall patient cohort (P = 0.461). Notably, 19 out of 29 immune-related risk terms differed substantially in the overall patient dataset. These risk terms mainly included checkpoints, human leukocyte antigens, natural killer cells, dendritic cells, and major histocompatibility complex class I.

Conclusion:

In summary, an immune-related prognostic gene signature was successfully developed and used to predict survival outcomes and immune system status in patients with HCC. This signature has the potential to help guide immunotherapeutic treatment planning for patients affected by this deadly cancer.

Development and Validation of a Prognostic Classifier Based on Lipid Metabolism-Related Genes in Gastric Cancer

Background: Dysregulation of lipid metabolism plays important roles in the tumorigenesis and progression of gastric cancer (GC). The present study aimed to establish a prognostic model based on the lipid metabolism–related genes in GC patients.

Materials and Methods: Two GC datasets from the Gene Expression Atlas, GSE62254 (n = 300) and GSE26942 (n = 217), were used as training and validation cohorts to establish a risk predictive scoring model. The efficacy of this model was assessed by ROC analysis. The association of the risk predictive scores with patient characteristics and immune cell subtypes was evaluated. A nomogram was constructed based on the risk predictive score model and other prognostic factors.

Results: A risk predictive score model was established based on the expression of 19 lipid metabolism–related genes (LPL, IPMK, PLCB3, CDIPT, PIK3CA, DPM2, PIGZ, GPD2, GPX3, LTC4S, CYP1A2, GALC, SGMS1, SMPD2, SMPD3, FUT6, ST3GAL1, B4GALNT1, and ACADS). The time-dependent ROC analysis revealed that the risk predictive score model was stable and robust. Patients with high risk scores had significantly unfavorable overall survival compared with those with low risk scores in both the training and validation cohorts. A higher risk score was associated with more aggressive features, including a higher tumor grade, a more advanced TNM stage, and diffuse type of Lauren classification of GC. Moreover, distinct immune cell subtypes and signaling pathways were found between the high–risk and low–risk score groups. A nomogram containing patients’ age, tumor stage, adjuvant chemotherapy, and the risk predictive score could accurately predict the survival probability of patients at 1, 3, and 5 years.

Conclusion: A novel 19-gene risk predictive score model was developed based on the lipid metabolism–related genes, which could be a potential prognostic indicator and therapeutic target of GC.

Pan-Cancer Analysis of FURIN as a Potential Prognostic and Immunological Biomarker

Background

Furin is a calcium-dependent protease that processes various precursor proteins through diverse secretory pathways. The deregulation of FURIN correlated with the prognosis of patients in numerous diseases. However, the role of FURIN in human pan-cancer is still largely unknown.

Methods

Multiple bioinformatic methods were employed to comprehensively analyze the correlation of FURIN expression with prognosis, mismatch repair (MMR), microsatellite instability (MSI), tumor mutational burden (TMB), DNA methylation, tumor immune infiltration, and common immune checkpoint inhibitors (ICIs) from the public database, and aim to evaluate the potential prognostic value of FURIN across cancers.

Results

FURIN was aberrantly expressed and was strongly correlated with MMR, MSI, TMB, and DNA methylation in multiple types of cancer. Moreover, survival analysis across cancers revealed that FURIN expression was correlated with overall survival (OS) in four cancers, disease-specific survival (DSS) in five cancers, progression-free interval (PFI) in seven cancers, and disease-free interval (DFI) in two cancers. Also, FURIN expression was related to immune cell infiltration in 6 cancers and ImmuneScore/StromalScore in 10 cancers, respectively. In addition, FURIN expression also showed strong association between expression levels and immune checkpoint markers in three cancers.

Conclusion

FURIN can serve as a significant prognostic biomarker and correlate with tumor immunity in human pan-cancer.

Construction of Gene Modules and Analysis of Prognostic Biomarkers for Cervical Cancer by Weighted Gene Co-Expression Network Analysis

Background

Despite advances in the understanding of neoplasm, patients with cervical cancer still have a poor prognosis. Identifying prognostic markers of cervical cancer may enable early detection of recurrence and more effective treatment.

Methods

Gene expression profiling data were acquired from the Gene Expression Omnibus database. After data normalization, genes with large variation were screened out. Next, we built co-expression modules by using weighted gene co-expression network analysis to investigate the relationship between the modules and clinical traits related to cervical cancer progression. Functional enrichment analysis was also applied on these co-expressed genes. We integrated the genes into a human protein-protein interaction (PPI) network to expand seed genes and build a co-expression network. For further analysis of the dataset, the Cancer Genome Atlas (TCGA) database was used to identify seed genes and their correlation to cervical cancer prognosis. Verification was further conducted by qPCR and the Human Protein Atlas (HPA) database to measure the expression of hub genes.

Results

Using WGCNA, we identified 25 co-expression modules from 10,016 genes in 128 human cervical cancer samples. After functional enrichment analysis, the magenta, brown, and darkred modules were selected as the three most correlated modules for cancer progression. Additionally, seed genes in the three modules were combined with a PPI network to identify 31 tumor-specific genes. Hierarchical clustering and Gepia results indicated that the expression quantity of hub genes NDC80, TIPIN, MCM3, MCM6, POLA1, and PRC1 may determine the prognosis of cervical cancer. Finally, TIPIN and POLA1 were further filtered by a LASSO model. In addition, their expression was identified by immunohistochemistry in HPA database as well as a biological experiment.

Conclusion

Our research provides a co-expression network of gene modules and identifies TIPIN and POLA1 as stable potential prognostic biomarkers for cervical cancer.

Identification of an immune-based mRNA-lncRNA signature for overall survival in cervical squamous cell carcinoma

Aim: To better predict the survival of cervical squamous cell carcinoma (CESC) patients, we aimed to construct a signature according to different immune infiltration. Methods: We downloaded the RNA sequences of CESC patients from the Cancer Genome Atlas database. By using single-sample gene set enrichment analysis, we separated the samples into high- and low-immunity groups. Then we separated the samples into training and testing datasets and performed the following analyses: univariate, least absolute shrinkage and selection operator analysis, multivariate Cox regression analyses and weighted gene coexpression network analysis using R software. Gene ontology and Kyoto Encyclopedia of Genes and Genomes studies were performed using the Database for Annotation, Visualization and Integrated Discovery website. Results & conclusion: We finally identified a signature with three mRNAs and two lncRNAs: ADGRG5, HSH2D, ZMAT4, RBAKDN and LINC00200. In short, our study constructed an mRNA-lncRNA signature related to immune infiltration to better predict the survival of CESC patients.

Identification of a Six-Gene Signature for Predicting the Overall Survival of Cervical Cancer Patients

Background

Although the incidence of cervical cancer has decreased in recent decades with the development of human papillomavirus vaccines and cancer screening, cervical cancer remains one of the leading causes of cancer-related death worldwide. Identifying potential biomarkers for cervical cancer treatment and prognosis prediction is necessary.

Methods

Samples with mRNA sequencing, copy number variant, single nucleotide polymorphism and clinical follow-up data were downloaded from The Cancer Genome Atlas database and randomly divided into a training dataset (N=146) and a test dataset (N=147). We selected and identified a prognostic gene set and mutated gene set and then integrated the two gene sets with the random survival forest algorithm and constructed a prognostic signature. External validation and immunohistochemical staining were also performed.

Results

We obtained 1416 differentially expressed prognosis-related genes, 624 genes with copy number amplification, 1038 genes with copy number deletion, and 163 significantly mutated genes. A total of 75 candidate genes were obtained after overlapping the differentially expressed genes and the genes with genomic variations. Subsequently, we obtained six characteristic genes through the random survival forest algorithm. The results showed that high expression of SLC19A3, FURIN, SLC22A3, and DPAGT1 and low expression of CCL17 and DES were associated with a poor prognosis in cervical cancer patients. We constructed a six-gene signature that can separate cervical cancer patients according to their different overall survival rates, and it showed robust performance for predicting survival (training set: p ˂ 0.001, AUC = 0.82; testing set: p ˂ 0.01, AUC = 0.59).

Conclusion

Our study identified a novel six-gene signature and nomogram for predicting the overall survival of cervical cancer patients, which may be beneficial for clinical decision-making for individualized treatment.

Tumor infiltrating immune cells (TIICs) as a biomarker for prognosis benefits in patients with osteosarcoma

BACKGROUND

Osteosarcoma is a rare malignant bone tumor in adolescents and children. Poor prognosis has always been a difficult problem for patients with osteosarcoma. Recent studies have shown that tumor infiltrating immune cells (TIICs) are associated with the clinical outcome of osteosarcoma patients. The aim of our research was to construct a risk score model based on TIICs to predict the prognosis of patients with osteosarcoma.

METHODS

CIBERSORTX algorithm was used to calculate the proportion of 22 TIIC types in osteosarcoma samples. Kaplan-Meier curves were drawn to investigate the prognostic value of 22 TIIC types. Forward stepwise approach was used to screen a minimal set of immune cell types. Multivariate Cox PHR analysis was performed to construct an immune risk score model.

RESULTS

Osteosarcoma samples with CIBERSORTX output p value less than 0.05 were selected for research. Kaplan-Meier curves indicated that naive B cells (p = 0.047) and Monocytes (p = 0.03) in osteosarcoma are associated with poor prognosis. An immune risk score model was constructed base on eight immune cell types, and the ROC curve showed that the immune risk score model is reliable in predicting the prognosis of patients with osteosarcoma (AUC = 0.724). Besides, a nomogram model base on eight immune cell types was constructed to predict the survival rate of patients with osteosarcoma.

CONCLUSIONS

TIICs are closely related to the prognosis of osteosarcoma. The immune risk score model based on TIICs is reliable in predicting the prognosis of osteosarcoma.

A five-mRNA signature associated with post-translational modifications can better predict recurrence and survival in cervical cancer

High mortality of patients with cervical cancer (CC) stresses the imperative of prognostic biomarkers for CC patients. Additionally, the vital status of post‐translational modifications (PTMs) in the progression of cancers has been reported by numerous researches. Therefore, the purpose of this research was to dig a prognostic signature correlated with PTMs for CC. We built a five‐mRNA (GALNTL6, ARSE, DPAGT1, GANAB and FURIN) prognostic signature associated with PTMs to predict both disease‐free survival (DFS) (hazard ratio [HR] = 3.967, 95% CI = 1.985‐7.927; P < .001) and overall survival (HR = 2.092, 95% CI = 1.138‐3.847; P = .018) for CC using data from The Cancer Genome Atlas database. Then, the robustness of the signature was validated using GSE44001 and the Human Protein Atlas (HPA) database. CIBERSORT algorithm analysis displayed that activated CD4 memory T cell was also an independent indicator for DFS (HR = 0.426, 95% CI = 0.186‐0.978; P = .044) which could add additional prognostic value to the signature. Collectively, the PTM‐related signature and activated CD4 memory T cell can provide new avenues for the prognostic predication of CC. These findings give further insights into effective treatment strategies for CC, providing opportunities for further experimental and clinical validations.