Data mining combines bioinformatics discover immunoinfiltration-related gene SERPINE1 as a biomarker for diagnosis and prognosis of stomach adenocarcinoma

Pan-cancer analysis of SERPINE1 with a concentration on immune therapeutic and prognostic in gastric cancer

The serine protease inhibitor clade E member 1 (SERPINE1) is a key modulator of the plasminogen/plasminase system and has been demonstrated to promote tumor progression and metastasis in various tumours. However, although much literature has explored the cancer-promoting mechanism of SERPINE1, the pan-cancer analyses of its predictive value and immune response remain unexplored. The differential expression, and survival analysis of SERPINE1 expression in multiple cancers were analysed using The Cancer Genome Atlas and Genotype-Tissue Expression database. Kaplan-Meier (K-M) plotter and survival data analysis were used to analyze the prognostic value of SERPINE1 expression, including overall survival (OS), disease-specific survival, disease-free interval and progression-free interval and investigated the relationship of SERPINE1 expression with microsatellite instability. We further analysed the correlation between the expression of SERPINE1 and immune infiltration. The Kyoto Encyclopaedia of Genes and Genomes pathway was used for enrichment analysis, and the Gene Set Enrichment Analysis (GSEA) database was used to perform pathway analysis. Finally, in vitro experiments demonstrated that knockdown or overexpression of SERPINE1 could alter the proliferation and migration of gastric cancer (GC) cells. The results indicated that SERPINE1 expression levels different significantly between cancer and normal tissues, meanwhile, it was highly expressed in various cancers. By analysing online data, it has been observed that the gene SERPINE1 exhibits heightened expression levels across a variety of human cancers, significantly impacting patient survival rates. Notably, the presence of SERPINE1 was strongly associated with decrease OS and disease-free survival in individuals diagnosed with GC. Furthermore, an observed link indicates that higher levels of SERPINE expression are associated with increased infiltration of immune cells in GC. Finally, in vitro experiments showed that knockdown or overexpression of SERPINE1 inhibited the growth, and migration, of GC cells. SERPINE1expression potentially represents a novel prognostic biomarker due to its significant association with immune cell infiltration in GC. This study shows that SERPINE1 is an oncogene that participates in regulating the immune infiltration and affecting the prognosis of patients in multiple cancers, especially in GC. These findings underscore the importance of further investigating the role of SERPINE1 in cancer progression and offer a promising direction for the development of new therapeutic strategies.

Identification and functional analysis of lactic acid metabolism-related differentially expressed genes in hepatocellular carcinoma

Background: Hepatocellular carcinoma (HCC) is a malignant tumor with high morbidity and mortality rate that seriously threatens human health. We aimed to investigate the expression, prognostic value, and immune cell infiltration of lactic acid metabolism-related genes (LAMRGs) in HCC using bioinformatics. Methods: The HCC database (The Cancer Genome Atlas-Liver Hepatocellular Carcinoma) was downloaded from the Cancer Genome Atlas (TCGA). Differentially expressed genes (DEGs) between normal and tumor groups were identified. The LAMRGs were obtained from literature and GeneCards and MSigDB databases. Lactic acid metabolism-related differentially expressed genes (LAMRDEGs) in HCC were screened from the DEGs and LAMRGs. Functional enrichment analyses of the screened LAMRDEGs were further conducted using Gene Ontology (GO) analysis, Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, and Gene Set Enrichment Analysis (GSEA). The genes were used in multivariate Cox regression and least absolute shrinkage and selection operator (LASSO) analyses to construct a prognostic model. Then, a protein-protein interaction network was constructed using STRING and CTD databases. Furthermore, the CIBERSORTx online database was used to assess the relationship between immune cell infiltration and hub genes. Results: Twenty-eight lactic acid metabolism-related differentially expressed genes (LAMRDEGs) were identified. The GO and KEGG analyses showed that the LAMRDEGs were related to the prognosis of HCC. The GSEA indicated that the LAMRDEGs were significantly enriched in tumor related pathways. In the multivariate Cox regression analysis, 14 key genes (E2F1, SERPINE1, GYS2, SPP1, PCK1, CCNB1, CYP2C9, IGFBP3, KDM8, RCAN1, ALPL, FBP1, NQO1, and LCAT) were found to be independent prognostic factors of HCC. Finally, the LASSO and Cox regression analyses showed that six key genes (SERPINE1, SPP1, CCNB1, CYP2C9, NQO1, and LCAT) were associated with HCC prognosis. Moreover, the correlation analyses revealed that the expression of the six key genes were associated with immune infiltrates of HCC. Conclusion: The LAMRDEGs can predict the prognosis and may be associated with immune cells infiltration in patients with HCC. These genes might be the promising biomarkers for the prognosis and treatment of HCC.

Regulatory mechanisms of dopamine metabolism in a marine Meyerozyma guilliermondii GXDK6 under NaCl stress as revealed by integrative multi-omics analysis

Dopamine can be used to treat depression, myocardial infarction, and other diseases. However, few reports are available on the de novo microbial synthesis of dopamine from low-cost substrate. In this study, integrated omics technology was used to explore the dopamine metabolism of a novel marine multi-stress-tolerant aromatic yeast Meyerozyma guilliermondii GXDK6. GXDK6 was found to have the ability to biosynthesize dopamine when using glucose as the substrate. 14 key genes for the biosynthesis of dopamine were identified by whole genome-wide analysis. Transcriptomic and proteomic data showed that the expression levels of gene AAT2 encoding aspartate aminotransferase (regulating dopamine anabolism) were upregulated, while gene AO-I encoding copper amine oxidase (involved in dopamine catabolism) were downregulated under 10 % NaCl stress compared with non-NaCl stress, thereby contributing to biosynthesis of dopamine. Further, the amount of dopamine under 10 % NaCl stress was 2.51-fold higher than that of zero NaCl, which was consistent with the multi-omics results. Real-time fluorescence quantitative PCR (RT-qPCR) and high-performance liquid chromatography (HPLC) results confirmed the metabolic model of dopamine. Furthermore, by overexpressing AAT2, AST enzyme activity was increased by 24.89 %, the expression of genes related to dopamine metabolism was enhanced, and dopamine production was increased by 56.36 % in recombinant GXDK6AAT2. In conclusion, Meyerozyma guilliermondii GXDK6 could utilize low-cost carbon source to synthesize dopamine, and NaCl stress promoted the biosynthesis of dopamine.

Bioinformatics analysis identifies coagulation factor II receptor as a potential biomarker in stomach adenocarcinoma

Coagulation factor 2 thrombin receptor (F2R), a member of the G protein-coupled receptor family, plays an important role in regulating blood clotting through protein hydrolytic cleavage mediated receptor activation. However, the underlying biological mechanisms by which F2R affects the development of gastric adenocarcinoma are not fully understood. This study aimed to systematically analyze the role of F2R in gastric adenocarcinoma. Stomach adenocarcinoma (STAD)-related gene microarray data and corresponding clinicopathological information were downloaded from the Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases. Differential expression genes (DEGs) associated with F2R were analyzed using Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), gene set enrichment analysis (GSEA), and protein-protein interaction (PPI) networks. F2R mRNA expression data were utilized to estimate stromal cell and immune cell scores in gastric cancer tissue samples, including stromal score, immune score, and ESTIMATE score, derived from single-sample enrichment studies. Analysis of TCGA and GEO databases revealed significantly higher F2R expression in STAD tissues compared to normal tissues. Patients with high F2R expression had shorter survival times than those with low F2R expression. F2R expression was significantly correlated with tumor (T) stage, node (N) stage, histological grade and pathological stage. Enrichment analysis of F2R-related genes showed that GO terms were mainly related to circulation-mediated human immune response, immunoglobulin, cell recognition and phagocytosis. KEGG analysis indicated associations to extracellular matrix (ECM) receptor interactions, neuroactive ligand-receptor interactions, the phosphoinositide-3-kinase-protein kinase B/Akt (PI3K-AKT) signaling pathway, the Wnt signaling pathway and the transforming growth factor-beta (TGF-β) signaling pathway. GSEA revealed connections to DNA replication, the Janus kinase/signal transducers and activators of transcription (JAK-STAT) signaling pathway, the mitogen-activated protein kinase (MAPK) signaling pathway and oxidative phosphorylation. Drug sensitivity analysis demonstrated positive correlations between F2R and several drugs, including BEZ235, CGP-60474, Dasatinib, HG-6-64-1, Aazopanib, Rapamycin, Sunitinib and TGX221, while negative correlation with CP724714, FH535, GSK1904529A, JNK-9L, LY317615, pyrimidine, rTRAIL and Vinorelbine. Knocking down F2R in GC cell lines resulted in slowed proliferation, migration, and invasion. All statistical analyses were performed using R software (version 4.2.1) and GraphPad Prism 9.0. p < 0.05 was considered statistically significant. In conclusion, this study underscores the significance of F2R as a potential biomarker in gastric adenocarcinoma, shedding light on its molecular mechanisms in tumorigenesis. F2R holds promise for aiding in the diagnosis, prognosis, and targeted therapy of STAD.

Pan-cancer analysis of SERPINE family genes as biomarkers of cancer prognosis and response to therapy

Background: Serine protease inhibitor E (SERPINE) family genes participate in the tumor growth, cancer cell survival and metastasis. However, the SERPINE family members role in the prognosis and their clinical therapeutic potentials in various human cancer types have not been elaborately explored. Methods: We preliminarily analyzed expression levels and prognostic values of SERPINE family genes, and investigated the correlation between SERPINEs expression and tumor microenvironment (TME), Stemness score, clinical characteristic, immune infiltration, tumor mutational burden (TMB), immune subtype, and drug sensitivity in pan-cancer, which based on updated public databases and integrated some bioinformatics analysis methods. In addition, we conducted the enrichment analysis of SERPINEs from DAVID and KOBAS databases. Results: SERPINE1, SERPINE2, and SERPINE3 expression were upregulated in nine cancers, twelve cancers, and six cancers, respectively. The expression of SERPINE family genes was associated with the prognosis in several cancers from The Cancer Genome Atlas (TCGA). Furthermore, SERPINE family genes expression also had a significant relation to stromal and immune scores, and RNA stemness score and DNA stemness score in pan-cancer. SERPINE1 and SERPINE2 expression significantly increased in tumor advanced stage in colon adenocarcinoma (COAD). Results showed that SERPINE1 and SERPINE2 expression were negatively related with B cells and Monocytes, respectively. SERPINE2 expression had a significantly positive relation with B cells and Macrophages. In terms of TMB, SERPINE1, SERPINE2, and SERPINE3 were found to associated with TMB in seven cancers, fourteen cancers, and four cancers, respectively. Moreover, all SERPINE gene family members were significantly correlated with immune subtypes. SERPINE1 expression had a significantly positive or negative correlation with drug sensitivity. Conclusion: The study indicated the great potential of SERPINE family genes as biomarkers for prognosis and provided valuable strategies for further investigation of SERPINE family genes as potential targets in cancer.

Characterization of the Prognosis and Tumor Microenvironment of Cellular Senescence-related Genes through scRNA-seq and Bulk RNA-seq Analysis in GC

BACKGROUND

Cellular senescence (CS) is thought to be the primary cause of cancer development and progression. This study aimed to investigate the prognostic role and molecular subtypes of CS-associated genes in gastric cancer (GC).

MATERIALS AND METHODS

The CellAge database was utilized to acquire CS-related genes. Expression data and clinical information of GC patients were obtained from The Cancer Genome Atlas (TCGA) database. Patients were then grouped into distinct subtypes using the "Consesus- ClusterPlus" R package based on CS-related genes. An in-depth analysis was conducted to assess the gene expression, molecular function, prognosis, gene mutation, immune infiltration, and drug resistance of each subtype. In addition, a CS-associated risk model was developed based on Cox regression analysis. The nomogram, constructed on the basis of the risk score and clinical factors, was formulated to improve the clinical application of GC patients. Finally, several candidate drugs were screened based on the Cancer Therapeutics Response Portal (CTRP) and PRISM Repurposing dataset.

RESULTS

According to the cluster result, patients were categorized into two molecular subtypes (C1 and C2). The two subtypes revealed distinct expression levels, overall survival (OS) and clinical presentations, mutation profiles, tumor microenvironment (TME), and drug resistance. A risk model was developed by selecting eight genes from the differential expression genes (DEGs) between two molecular subtypes. Patients with GC were categorized into two risk groups, with the high-risk group exhibiting a poor prognosis, a higher TME level, and increased expression of immune checkpoints. Function enrichment results suggested that genes were enriched in DNA repaired pathway in the low-risk group. Moreover, the Tumor Immune Dysfunction and Exclusion (TIDE) analysis indicated that immunotherapy is likely to be more beneficial for patients in the low-risk group. Drug analysis results revealed that several drugs, including ML210, ML162, dasatinib, idronoxil, and temsirolimus, may contribute to the treatment of GC patients in the high-risk group. Moreover, the risk model genes presented a distinct expression in single-cell levels in the GSE150290 dataset.

CONCLUSION

The two molecular subtypes, with their own individual OS rate, expression patterns, and immune infiltration, lay the foundation for further exploration into the GC molecular mechanism. The eight gene signatures could effectively predict the GC prognosis and can serve as reliable markers for GC patients.

9S1R nullomer peptide induces mitochondrial pathology, metabolic suppression, and enhanced immune cell infiltration, in triple-negative breast cancer mouse model

Nullomers are the shortest strings of absent amino acid (aa) sequences in a species or group of species. Primes are those nullomers that have not been detected in the genome of any species. 9S1R is a 5-aa peptide prime sequence attached to 5-arginine aa, used to treat triple negative breast cancer (TNBC) in an in vivo mouse model. This unique peptide, administered with a trehalose carrier (9S1R-NulloPT), offers enhanced solubility and exhibits distinct anti-cancer effects against TNBC. In our study, we investigated the effect of 9S1R-NulloPT on tumor growth, metabolism, metastatic burden, tumor immune-microenvironment (TME), and transcriptome of aggressive mouse TNBC tumors. Notably, treated mice had smaller tumors in the initial phase of the treatment, as compared to untreated control, and diminished in vivo and ex vivo bioluminescence at later-stages - indicative of metabolically quiescent, dying tumors. The treatment also caused changes in TME with increased infiltration of immune cells and altered tumor transcriptome, with 365 upregulated genes and 710 downregulated genes. Consistent with in vitro data, downregulated genes were enriched in cellular metabolic processes (179), specifically mitochondrial TCA cycle/oxidative phosphorylation (44), and translation machinery/ribosome biogenesis (45). The upregulated genes were associated with the developmental (13), ECM organization (12) and focal adhesion pathways (7). In conclusion, our study demonstrates that 9S1R-NulloPT effectively reduced tumor growth during its initial phase, altering the TME and tumor transcriptome. The treatment induced mitochondrial pathology which led to a metabolic deceleration in tumors, aligning with in vitro observations.

Integrated bioinformatics analysis of noncoding RNAs with tumor immune microenvironment in gastric cancer

In recent years, molecular and genetic research hotspots of gastric cancer have been investigated, including microRNAs, long noncoding RNAs (lncRNAs) and messenger RNA (mRNAs). The study on the role of lncRNAs may help to develop personalized treatment and identify potential prognostic biomarkers in gastric cancer. The RNA-seq and miRNA-seq data of gastric cancer were downloaded from the TCGA database. Differential analysis of RNA expression between gastric cancer samples and normal samples was performed using the edgeR package. The ceRNA regulatory network was visualized using Cytoscape. KEGG pathway analysis of mRNAs in the ceRNA network was performed using the clusterProfiler package. CIBERSORT was used to distinguish 22 immune cell types and the prognosis-related genes and immune cells were determined using Kaplan-Meier and Cox proportional hazard analyses. To estimate these nomograms, we used receiver operating characteristic and calibration curve studies. The ceRNA regulation network of gastric cancer was built in this study, and the genes in the network were analyzed for prognosis. A total of 980 lncRNAs were differentially expressed, of which 774 were upregulated and 206 were downregulated. A survival study identified 15 genes associated with gastric cancer prognosis, including VCAN-AS1, SERPINE1, AL139002.1, LINC00326, AC018781.1, C15orf54, hsa-miR-145. Monocytes and Neutrophils were associated with the survival rate of gastric cancer. Our research uncovers new ceRNA network for the detection, treatment, and monitoring of gastric cancer.

Single-cell RNA sequencing integrated with bulk RNA sequencing analysis reveals diagnostic and prognostic signatures and immunoinfiltration in gastric cancer

BACKGROUND

Early diagnosis and prognostic predication of gastric cancer (GC) pose significant challenges in current clinical practice of GC treatments. Therefore, our aim was to explore relevant gene signatures that can predict the prognosis of GC patients.

METHODS

Here, we established a single-cell transcriptional atlas of GC, focusing on the expression of T-cell-related genes for cell-cell communication analysis, trajectory analysis, and transcription factor regulatory network analysis. Additionally, we conducted validation and prediction of immune-related prognostic gene signatures in GC patients using TCGA and GEO data. Based on these prognostic gene signatures, we predicted the immune infiltration status of GC patients by grouping the patient samples into high or low-risk groups.

RESULTS

Based on 10 tumor samples and corresponding normal samples from GC patients, we selected 18,416 cells for subsequent analysis using single-cell sequencing. From these, we identified 3,284 T-cells and obtained 641 differentially expressed genes related to T-cells from 5 different T-cell subtypes. By integrating bulk RNA sequencing data, we identified prognostic signatures associated with T-cells. Stratifying patients based on these prognostic signatures into high-risk or low-risk groups allowed us to effectively predict their survival rates and the immunoinfiltration status of the tumor microenvironment.

CONCLUSION

This study explored prognostic gene signatures associated with T-cells in GC patients, providing insights into predicting patients' survival rates and immunoinfiltration levels.

Nullomer peptide increases immune cell infiltration and reduces tumor metabolism in triple negative breast cancer mouse model

Background

Nullomers are the shortest strings of absent amino acid (aa) sequences in a species or group of species. Primes are those nullomers that have not been detected in the genome of any species. 9S1R is a 5-aa peptide derived from a prime sequence that is tagged with 5 arginine aa, used to treat triple negative breast cancer (TNBC) in an in vivo TNBC mouse model. 9S1R is administered in trehalose (9S1R-NulloPT), which enhances solubility and exhibits some independent effects against tumor growth and is thus an important component in the drug preparation.

Method

We examined the effect of 9S1R-NulloPT on tumor growth, metabolism, metastatic burden, necrosis, tumor immune microenvironment, and the transcriptome of aggressive mouse TNBC tumors.

Results

The peptide-treated mice had smaller tumors in the initial phase of the treatment, as compared to the untreated control, and reduced in vivo bioluminescence at later stages, which is indicative of metabolically inactive tumors. A decrease in ex vivo bioluminescence was also observed in the excised tumors of treated mice, but not in the secondary metastasis in the lungs. The treatment also caused changes in tumor immune microenvironment with increased infiltration of immune cells and margin inflammation. The treatment upregulated 365 genes and downregulated 710 genes in tumors compared to the untreated group. Consistent with in vitro findings in breast cancer cell lines, downregulated genes in the treated TNBC tumors include Cellular Metabolic Process Related genes (179), specifically mitochondrial genes associated with TCA cycle/oxidative phosphorylation (44), and translation machinery/ribosome biogenesis genes (45). Among upregulated genes, the Developmental Pathway (13), ECM Organization (12) and Focal Adhesion Related Pathways (7) were noteworthy. We also present data from a pilot study using a bilateral BC mouse model, which supports our findings.

Conclusion

In conclusion, although 9S1R-NulloPT was moderate at reducing the tumor volume, it altered the tumor immune microenvironment as well as the tumor transcriptome, rendering tumors metabolically less active by downregulating the mitochondrial function and ribosome biogenesis. This corroborates previously published in vitro findings.