Identification prognosis-associated immune genes in colon adenocarcinoma

An integrated framework for prognosis prediction and drug response modeling in colorectal liver metastasis drug discovery

BACKGROUND

Colorectal cancer (CRC) is the third most prevalent cancer globally, and liver metastasis (CRLM) is the primary cause of death. Hence, it is essential to discover novel prognostic biomarkers and therapeutic drugs for CRLM.

METHODS

This study developed two liver metastasis-associated prognostic signatures based on differentially expressed genes (DEGs) in CRLM. Additionally, we employed an interpretable deep learning model utilizing drug sensitivity databases to identify potential therapeutic drugs for high-risk CRLM patients. Subsequently, in vitro and in vivo experiments were performed to verify the efficacy of these compounds.

RESULTS

These two prognostic models exhibited superior performance compared to previously reported ones. Obatoclax, a BCL-2 inhibitor, showed significant differential responses between high and low risk groups classified by prognostic models, and demonstrated remarkable effectiveness in both Transwell assay and CT26 colorectal liver metastasis mouse model.

CONCLUSIONS

This study highlights the significance of developing specialized prognostication approaches and investigating effective therapeutic drugs for patients with CRLM. The application of a deep learning drug response model provides a new drug discovery strategy for translational medicine in precision oncology.

Identification of disulfidptosis-related subtypes, characterization of tumor microenvironment infiltration, and development of a prognosis model in breast cancer

Introduction

Breast cancer (BC) is now the most common type of cancer in women. Disulfidptosis is a new regulation of cell death (RCD). RCD dysregulation is causally linked to cancer. However, the comprehensive relationship between disulfidptosis and BC remains unknown. This study aimed to explore the predictive value of disulfidptosis-related genes (DRGs) in BC and their relationship with the TME.

Methods

This study obtained 11 disulfidptosis genes (DGs) from previous research by Gan et al. RNA sequencing data of BC were downloaded from the Cancer Genome Atlas (TCGA) and Gene Expression Omnibus database (GEO) databases. First, we examined the effect of DG gene mutations and copy number changes on the overall survival of breast cancer samples. We then used the expression profile data of 11 DGs and survival data for consensus clustering, and BC patients were divided into two clusters. Survival analysis, gene set variation analysis (GSVA) and ss GSEA were used to compare the differences between them. Subsequently, DRGs were identified between the clusters used to perform Cox regression and least absolute shrinkage and selection operator regression (LASSO) analyses to construct a prognosis model. Finally, the immune cell infiltration pattern, immunotherapy response, and drug sensitivity of the two subtypes were analyzed. CCK-8 and a colony assay obtained by knocking down genes and gene sequencing were used to validate the model.

Result

Two DG clusters were identified based on the expression of 11DGs. Then, 225 DRGs were identified between them. RS, composed of six genes, showed a significant relationship with survival, immune cell infiltration, clinical characteristics, immune checkpoints, immunotherapy response, and drug sensitivity. Low-RS shows a better prognosis and higher immunotherapy response than high-RS. A nomogram with perfect stability constructed using signature and clinical characteristics can predict the survival of each patient. CCK-8 and colony assay obtained by knocking down genes have demonstrated that the knockdown of high-risk genes in the RS model significantly inhibited cell proliferation.

Discussion

This study elucidates the potential relationship between disulfidptosis-related genes and breast cancer and provides new guidance for treating breast cancer.

A novel gene signature related to oxidative stress predicts the prognosis in clear cell renal cell carcinoma

Clear cell renal cell carcinoma (ccRCC) is considered to be related to the worse prognosis, which might in part be attributed to the early recurrence and metastasis, compared with other type of kidney cancer. Oxidative stress refers to an imbalance between production of oxidants and antioxidant defense. Accumulative studies have indicated that oxidative stress genes contribute to the tumor invasion, metastasis and drug sensitivity. However, the biological functions of oxidative stress genes in ccRCC remain largely unknown. In this study, we identified 1,399 oxidative stress genes from GeneCards with a relevance score ≥7. Data for analysis were accessed from The Cancer Genome Atlas (TCGA) and the International Cancer Genome Consortium (ICGC) database, and were utilized as training set and validation set respectively. Univariate Cox analysis, least absolute shrinkage and selection operator (LASSO) Cox regression and multivariate Cox were employed to construct a prognostic signature in ccRCC. Finally, a prognostic signature including four different oxidative stress genes was constructed from 1,399 genes, and its predictive performance was verified through Kaplan-Meier survival analysis and the receiver operating characteristic (ROC) curve. Interestingly, we found that there was significant correlation between the expression of oxidative stress genes and the immune infiltration and the sensitivity of tumor cells to chemotherapeutics. Moreover, the highest hazard ratio gene urocortin (UCN) was chosen for further study; some necessary vitro experiments proved that the UCN could promote the ability of ccRCC proliferation and migration and contribute to the degree of oxidative stress. In conclusion, it was promising to predict the prognosis of ccRCC through the four oxidative stress genes signature. UCN played oncogenic roles in ccRCC by influencing proliferation and oxidative stress pathway, which was expected to be the novel therapeutic target for ccRCC.

Regularized survival learning and cross-database analysis enabled identification of colorectal cancer prognosis-related immune genes

Colon adenocarcinoma is the most common type of colorectal cancer. The prognosis of advanced colorectal cancer patients who received treatment is still very poor. Therefore, identifying new biomarkers for prognosis prediction has important significance for improving treatment strategies. However, the power of biomarker analyses was limited by the used sample size of individual database. In this study, we combined Genotype-Tissue Expression (GTEx) and The Cancer Genome Atlas (TCGA) databases to expand the number of healthy tissue samples. We screened differentially expressed genes between the GTEx healthy samples and TCGA tumor samples. Subsequently, we applied least absolute shrinkage and selection operator (LASSO) regression and multivariate Cox analysis to identify nine prognosis-related immune genes: ANGPTL4, IDO1, NOX1, CXCL3, LTB4R, IL1RL2, CD72, NOS2, and NUDT6. We computed the risk scores of samples based on the expression levels of these genes and divided patients into high- and low-risk groups according to this risk score. Survival analysis results showed a significant difference in survival rate between the two risk groups. The high-risk group had a significantly lower overall survival rate and poorer prognosis. We found the receiver operating characteristic based on the risk score was showed to accurately predict patients' prognosis. These prognosis-related immune genes may be potential biomarkers for colorectal cancer diagnosis and treatment. Our open-source code is freely available from GitHub at https://github.com/gutmicrobes/Prognosis-model.git.

Establishment and validation of an aging-related risk signature associated with prognosis and tumor immune microenvironment in breast cancer

BACKGROUND

Breast cancer (BC) is a highly malignant and heterogeneous tumor which is currently the cancer with the highest incidence and seriously endangers the survival and prognosis of patients. Aging, as a research hotspot in recent years, is widely considered to be involved in the occurrence and development of a variety of tumors. However, the relationship between aging-related genes (ARGs) and BC has not yet been fully elucidated.

MATERIALS AND METHODS

The expression profiles and clinicopathological data were acquired in the Cancer Genome Atlas (TCGA) and the gene expression omnibus (GEO) database. Firstly, the differentially expressed ARGs in BC and normal breast tissues were investigated. Based on these differential genes, a risk model was constructed composed of 11 ARGs via univariate and multivariate Cox analysis. Subsequently, survival analysis, independent prognostic analysis, time-dependent receiver operating characteristic (ROC) analysis and nomogram were performed to assess its ability to sensitively and specifically predict the survival and prognosis of patients, which was also verified in the validation set. In addition, functional enrichment analysis and immune infiltration analysis were applied to reveal the relationship between the risk scores and tumor immune microenvironment, immune status and immunotherapy. Finally, multiple datasets and real-time polymerase chain reaction (RT-PCR) were utilized to verify the expression level of the key genes.

RESULTS

An 11-gene signature (including FABP7, IGHD, SPIB, CTSW, IGKC, SEZ6, S100B, CXCL1, IGLV6-57, CPLX2 and CCL19) was established to predict the survival of BC patients, which was validated by the GEO cohort. Based on the risk model, the BC patients were divided into high- and low-risk groups, and the high-risk patients showed worse survival. Stepwise ROC analysis and Cox analyses demonstrated the good performance and independence of the model. Moreover, a nomogram combined with the risk score and clinical parameters was built for prognostic prediction. Functional enrichment analysis revealed the robust relationship between the risk model with immune-related functions and pathways. Subsequent immune microenvironment analysis, immunotherapy, etc., indicated that the immune status of patients in the high-risk group decreased, and the anti-tumor immune function was impaired, which was significantly different with those in the low-risk group. Eventually, the expression level of FABP7, IGHD, SPIB, CTSW, IGKC, SEZ6, S100B, CXCL1, IGLV6-57 and CCL19 was identified as down-regulated in tumor cell line, while CPLX2 up-regulated, which was mostly similar with the results in TCGA and Human Protein Atlas (HPA) via RT-PCR.

CONCLUSIONS

In summary, our study constructed a risk model composed of ARGs, which could be used as a solid model for predicting the survival and prognosis of BC patients. Moreover, this model also played an important role in tumor immunity, providing a new direction for patient immune status assessment and immunotherapy selection.

Development and validation of a chromatin regulator prognostic signature in colon adenocarcinoma

Aberrant expression of chromatin regulators (CRs) could lead to the development of various diseases including cancer. However, the biological function and prognosis role of CRs in colon adenocarcinoma (COAD) remains unclear. We performed the clustering analyses for expression profiling of COAD downloaded from The Cancer Genome Atlas. We developed a chromatin regulator prognostic model, which was validated in an independent cohort data. Time-intendent receiver operating characteristics curve was used to evaluate predict ability of model. Univariate and multivariate cox regression were used to assess independence of risk score. Nomogram was established to assess individual risk. Gene ontology, and Kyoto Encyclopedia of genes and genomes, gene set variation analysis and gene set enrichment analysis were performed to explore the function of CRs. Immune infiltration and drug sensitivity were also performed to assess effect of CRs on treatment in COAD. COAD can be separated into two subtypes with different clinical characteristics and prognosis. The C2 had elevated immune infiltration levels and low tumor purity. Using 12 chromatin regulators, we developed and validated a prognostic model that can predict the overall survival of COAD patients. We built a risk score that can be an independent prognosis predictor of COAD. The nomogram score system achieved the best predict ability and were also confirmed by decision curve analysis. There were significantly different function and pathway enrichment, immune infiltration levels, and tumor mutation burden between high-risk and low-risk group. The external validation data also indicated that high-risk group had higher stable disease/progressive disease response rate and poorer prognosis than low-risk group. Besides, the signature genes included in the model could cause chemotherapy sensitivity to some small molecular compounds. Our integrative analyses for chromatin regulators could provide new insights for the risk management and individualized treatment in COAD.

The role of a cuproptosis-related prognostic signature in colon cancer tumor microenvironment and immune responses

Background: Colon adenocarcinoma (COAD) is a common malignant tumor of the digestive tract with poor clinical outcomes. Cuproptosis is a novel cell death mechanism and linked to mitochondrial respiration. However, the role of cuproptosis in colon cancer tumor microenvironment (TME) and immune responses remains unknown.

Methods: We conducted difference analysis to identify the differential expressed cuproptosis-related genes (CRGs). According to the CRGs, the TCGA-COAD samples were categorized using consensus clustering. The LASSO regression analysis was utilized to develop the cuproptosis-related signature. We then verified the model reliability by Kaplan–Meier, PCA, and ROC analysis. The GES39582 cohort served as the validation set. GO and KEGG functional analyses were conducted to investigate the underlying mechanism. We compared the infiltration levels of immune cells, the expression levels of immune checkpoints, and microsatellite instability (MSI) status between the high- and low-risk groups. Additionally, the relationships between the risk signature and immune cells and cancer stem cell (CSC) were analyzed.

Results: Finally, we identified 9 differentially expressed CRGs in COAD. According to the expression of CRGs, the TCGA-COAD samples were separated into two clusters. The 11-gene signature was established by LASSO, and it had excellent predictive power for COAD prognosis. Besides, we used the GSE39582 cohort to validate the prognostic value of the model. GO and KEGG results demonstrated that the survival differences between two risk groups was mainly linked to the extracellular matrix (ECM). Further immune characterization analysis showed the significant differences in the immune cell infiltration and immune responses between two risk groups.

Conclusion: Overall, the novel cuproptosis-related signature was able to accurately predict COAD prognosis and played important roles in COAD tumor microenvironment and immune responses.

Novel six-gene prognostic signature based on colon adenocarcinoma immune-related genes

Background

Colon adenocarcinoma (COAD) is one of the most common gastrointestinal tumors worldwide, and immunotherapy is one of the most promising treatments for it. Identifying immune genes involved in the development and maintenance of cancer is key to the use of tumor immunotherapy. This study aimed to determine the prognostic value of immune genes in patients with COAD and to establish an immune-related gene signature. Differentially expressed genes, immune-related genes (DEIGs), and transcription factors (DETFs) were screened using the following databases: Cistrome, The Cancer Genome Atlas (TCGA), the Immunology Database and Analysis Portal, and InnateDB. We constructed a network showing the regulation of DEIGs by DETFs. Using weighted gene co-expression network analysis, we prepared 5 co-expressed gene modules; 6 hub genes (CD1A, CD1B, FGF9, GRP, SERPINE1, and F2RL2) obtained using univariate and multivariate regression analysis were used to construct a risk model. Patients from TCGA database were divided into high- and low-risk groups based on whether their risk score was greater or less than the mean; the public dataset GSE40967, which contains gene expression profiles of 566 colon cancer patients, was used for validation.

Results

Survival analysis, somatic gene mutations, and tumor-infiltrating immune cells differed significantly between the high- and low-risk groups.

Conclusions

This immune-related gene signature could play an important role in guiding treatment, making prognoses, and potentially developing future clinical applications.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-022-04909-2.

Identification of Prognostic Fatty Acid Metabolism lncRNAs and Potential Molecular Targeting Drugs in Uveal Melanoma

Background

The aim of this study was to identify prognostic fatty acid metabolism lncRNAs and potential molecular targeting drugs in uveal melanoma through integrated bioinformatics analysis.

Methods

In the present study, we obtained the expression matrix of 309 FAM-mRNAs and identified 225 FAM-lncRNAs by coexpression network analysis. We then performed univariate Cox analysis, LASSO regression analysis, and cross-validation and finally obtained an optimized UVM prognosis prediction model composed of four PFAM-lncRNAs (AC104129.1, SOS1-IT1, IDI2-AS1, and DLGAP1-AS2).

Results

The survival curves showed that the survival time of UVM patients in the high-risk group was significantly lower than that in the low-risk group in the train cohort, test cohort, and all patients in the prognostic prediction model (P < 0.05). We further performed risk prognostic assessment, and the results showed that the risk scores of the high-risk group in the train cohort, test cohort, and all patients were significantly higher than those of the low-risk group (P < 0.05), patient survival decreased and the number of deaths increased with increasing risk scores, and AC104129.1, SOS1-IT1, and DLGAP1-AS2 were high-risk PFAM-lncRNAs, while IDI2-AS1 were low-risk PFAM-lncRNAs. Afterwards, we further verified the accuracy and the prognostic value of our model in predicting prognosis by PCA analysis and ROC curves.

Conclusion

We identified 24 potential molecularly targeted drugs with significant sensitivity differences between high- and low-risk UVM patients, of which 13 may be potential targeted drugs for high-risk patients. Our findings have important implications for early prediction and early clinical intervention in high-risk UVM patients.

Development of a novel hypoxia-immune–related LncRNA risk signature for predicting the prognosis and immunotherapy response of colorectal cancer

Background

Colorectal cancer (CRC) is one of the most common digestive system tumors worldwide. Hypoxia and immunity are closely related in CRC; however, the role of hypoxia-immune–related lncRNAs in CRC prognosis is unknown.

Methods

Data used in the current study were sourced from the Gene Expression Omnibus and The Cancer Genome Atlas (TCGA) databases. CRC patients were divided into low- and high-hypoxia groups using the single-sample gene set enrichment analysis (ssGSEA) algorithm and into low- and high-immune groups using the Estimation of STromal and Immune cells in MAlignant Tumours using Expression data (ESTIMATE) algorithm. Differentially expressed lncRNAs (DElncRNAs) between low- and high-hypoxia groups, low- and high-immune groups, and tumor and control samples were identified using the limma package. Hypoxia-immune–related lncRNAs were obtained by intersecting these DElncRNAs. A hypoxia-immune–related lncRNA risk signature was developed using univariate Cox regression and least absolute shrinkage and selection operator (LASSO) analyses. The tumor microenvironments in the low- and high-risk groups were evaluated using ssGSEA, ESTIMATE, and the expression of immune checkpoints. The therapeutic response in the two groups was assessed using TIDE, IPS, and IC50. A ceRNA network based on signature lncRNAs was constructed. Finally, we used RT-qPCR to verify the expression of hypoxia-immune–related lncRNA signatures in normal and cancer tissues.

Results

Using differential expression analysis, and univariate Cox and LASSO regression analyses, ZNF667-AS1, LINC01354, LINC00996, DANCR, CECR7, and LINC01116 were selected to construct a hypoxia-immune–related lncRNA signature. The performance of the risk signature in predicting CRC prognosis was validated in internal and external datasets, as evidenced by receiver operating characteristic curves. In addition, we observed significant differences in the tumor microenvironment and immunotherapy response between low- and high-risk groups and constructed a CECR7–miRNA–mRNA regulatory network in CRC. Furthermore, RT-qPCR results confirmed that the expression patterns of the six lncRNA signatures were consistent with those in TCGA-CRC cohort.

Conclusion

Our study identified six hypoxia-immune–related lncRNAs for predicting CRC survival and sensitivity to immunotherapy. These findings may enrich our understanding of CRC and help improve CRC treatment. However, large-scale long-term follow-up studies are required for verification.

Therapeutic targets and functions of curcumol against COVID-19 and colon adenocarcinoma

Since 2019, the coronavirus disease (COVID-19) has caused 6,319,395 deaths worldwide. Although the COVID-19 vaccine is currently available, the latest variant of the virus, Omicron, spreads more easily than earlier strains, and its mortality rate is still high in patients with chronic diseases, especially cancer patients. So, identifying a novel compound for COVID-19 treatment could help reduce the lethal rate of the viral infection in patients with cancer. This study applied network pharmacology and systematic bioinformatics analysis to determine the possible use of curcumol for treating colon adenocarcinoma (COAD) in patients infected with COVID-19. Our results showed that COVID-19 and COAD in patients shared a cluster of genes commonly deregulated by curcumol. The clinical pathological analyses demonstrated that the expression of gamma-aminobutyric acid receptor subunit delta (GABRD) was associated with the patients' hazard ratio. More importantly, the high expression of GABRD was associated with poor survival rates and the late stages of COAD in patients. The network pharmacology result identified seven-core targets, including solute carrier family 6 member 3, gamma-aminobutyric acid receptor subunit pi, butyrylcholinesterase, cytochrome P450 3A4, 17-beta-hydroxysteroid dehydrogenase type 2, progesterone receptor, and GABRD of curcumol for treating patients with COVID-19 and COAD. The bioinformatic analysis further highlighted their importance in the biological processes and molecular functions in gland development, inflammation, retinol, and steroid metabolism. The findings of this study suggest that curcumol could be an alternative compound for treating patients with COVID-19 and COAD.

Identification and validation of prognosis-associated DNA repair gene signatures in colorectal cancer

Colorectal cancer (CRC) is the third most common malignant tumor. DNA damage plays a crucial role in tumorigenesis, and abnormal DNA repair pathways affect the occurrence and progression of CRC. In the current study, we aimed to construct a DNA repair-related gene (DRG) signature to predict the overall survival (OS) of patients with CRC patients. The differentially expressed DRGs (DE-DRGs) were analyzed using The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases. The prognostic gene signature was identified by univariate Cox regression and least absolute shrinkage and selection operator (LASSO)-penalized Cox proportional hazards regression analysis. The predictive ability of the model was evaluated using the Kaplan–Meier curves and time-dependent receiver operating characteristic (ROC) curves. The gene set enrichment analysis (GSEA) was performed to explore the underlying biological processes and signaling pathways. ESTIMATE and CIBERSORT were implemented to estimate the tumor immune score and immune cell infiltration status between the different risk group. The half-maximal inhibitory concentration (IC50) was evaluated to representing the drug response of this signature. Nine DE-DRGs (ESCO2, AXIN2, PLK1, CDC25C, IGF1, TREX2, ALKBH2, ESR1 and MC1R) signatures was constructed to classify patients into high- and low-risk groups. The risk score was an independent prognostic indicator of OS (hazard ratio > 1, P < 0.001). The genetic alteration analysis indicated that the nine DE-DRGs in the signature were changed in 63 required samples (100%), and the major alteration was missense mutation. Function enrichment analysis revealed that the immune response and mtotic sister chromatid segregation were the main biological processes. The high-risk group had higher immune score than the low-risk group. What’s more, low-risk patients were more sensitive to selumetinib and dasatinib. The nine DE-DRGs signature was significantly associated with OS and provided a new insight for the diagnosis and treatment of CRC.

Centromere Protein A (CENPA) Regulates Metabolic Reprogramming in the Colon Cancer Cells by Transcriptionally Activating Karyopherin Subunit Alpha 2 (KPNA2)

The karyopherin α₂ subunit gene (KPNA2) has been reported as an oncogene and is involved in metabolic reprogramming in cancer. This study aimed to explore the function of KPNα₂ in the growth and glycolysis in colon cancer (CC) cells. Genes from the Oncomine database that were differentially expressed in multiple CC types were screened. Bioinformatics analysis suggested that KPNA2 was highly expressed in CC: High expression of KPNA2 was detected in the CC cell lines. Down-regulation of KPNA2 reduced viability and DNA-replication ability, and it increased apoptosis of HCT116 and LoVo cells. It also reduced glucose consumption, extracellular acidification rate, and the ATP production in cells. Centromere protein A (CENPA) was confirmed as an upstream transcription activator of KPNA2. There was significant H3K27ac modification in the promoter region of KPNA2. CENPA mainly recruited histone acetyltransferase general control of amino acid synthesis (GCN)-5 to the promoter region of KPNA2 to induce transcription activation. Overexpression of either CENPA or GCN-5 blocked the role of short hairpin KPNα₂ and restored growth and glycolysis in CC cells. To conclude, the findings from this study suggest that CENPA recruits GCN-5 to the promoter region of KPNA2 to induce KPNα₂ activation, which strengthens growth and glycolysis in, and augments the development of, CC.

Immune cell infiltration-associated signature in colon cancer and its prognostic implications

Tumor immune cell infiltration (ICI) has been reported in various studies to be correlated with tumor diagnosis, clinical treatment sensitivity and prognosis. It is an important direction to study the characteristics of immune cell infiltration and develop new prognostic markers to improve the treatment of colon cancer. In this paper, we systematically analyzed the ICI characteristics and obtained three ICI clusters. Then, the ICI scores were constructed and its prognostic implications were discussed. From the results, the ICI score patterns were linked to a great survival difference (p<0.001). A high ICI score was characterized by a higher fraction of plasma cells, CD8⁺ T cells, memory resting CD4⁺ T cells, monocytes, eosinophils and dendritic cells, which had better prognosis. Macrophages and neutrophils were increased in low ICI score patients with decreased overall survival. Immune checkpoint molecules (PDCD1, CD274, LAG3, IDO1, CTLA-4, TIGHT and HAVCR2) were found to be significantly overexpressed in the low ICI score subgroup. In addition, we also studied the correlation between the tumor mutation burden (TMB) and ICI score. This study indicated the ICI score could serve as a potential prognostic biomarker for colon cancer patients’ immunotherapy.

Immune-Related Gene Expression Analysis Revealed Three lncRNAs as Prognostic Factors for Colon Cancer

Colorectal cancer (CRC) is one of the most common cancers. Almost 80% of CRC cases are colon adenocarcinomas (COADs). Several studies have indicated the role of immunotherapy in the treatment of various cancers. Our study aimed to identify immune-related long non-coding RNAs (lncRNAs) and to use them to construct a risk assessment model for evaluating COAD prognosis. Using differential expression, correlation, and Cox regression analyses, we identified three immune-related differentially expressed lncRNAs (IR-DELs) and used them to construct a risk assessment model. The area under the curve (AUC) for each receiver operating characteristic (ROC) curve at 3-, 5-, and 10-years were greater than 0.6. In addition, the risk assessment model was correlated with several immune cells and factors. The three IR-DELs (AC124067.4, LINC02604, and MIR4435-2HG) identified in this study can be used to predict outcomes for patients with COAD and might help in identifying those who can benefit from anti-tumor immunotherapy.

An integrative microenvironment approach for laryngeal carcinoma: the role of immune/methylation/autophagy signatures on disease clinical prognosis and single-cell genotypes

The effects of methylation/autophagy-related genes (MARGs) and immune infiltration in the tumor microenvironment on the prognosis of laryngeal cancer were comprehensively explored in this study. Survival analysis screened out 126 MARGs and 10 immune cells potentially associated with the prognosis of laryngeal carcinoma. Cox and lasso regression analyses were then used to select 8 MARGs (CAPN10, DAPK2, MBTPS2, ST13, CFLAR, FADD, PEX14 and TSC2) and 2 immune cells (Eosinophil and Mast cell) to obtain the prognostic risk scoring system (pRS). The pRS was used to establish a risk prediction model for the prognosis of laryngeal cancer. The predictive ability of the prediction model was evaluated by GEO datasets and our clinical samples. Further analysis revealed that pRS is highly associated with single nucleotide polymorphism (SNP), copy number variation (CNV), immune checkpoint blockade (ICB) therapy and tumor microenvironment. Moreover, the screened pRS-related ceRNA network and circ_0002951/miR-548k/HAS2 pathway provide potential therapeutic targets and biomarkers of laryngocarcinoma. Based on the clustering results of pRS-related genes, single cells were then genotyped and revealed by integrated scRNA-seq in laryngeal cancer samples. Fibroblasts were found enriched in high risk cell clusters at the scRNA-seq level. Fibroblast-related ligand-receptor interactions were then exposed and a neural network-based deep learning model based on these pRS-related hub gene signatures was also established with a high accuracy in cell type prediction. In conclusion, the combination of single-cell and transcriptome laryngeal carcinoma landscape analyses can investigate the link between the tumor microenvironmental and prognostic characteristics.

Identification of the hub genes in gastric cancer through weighted gene co-expression network analysis

Background

Gastric cancer is one of the most lethal tumors and is characterized by poor prognosis and lack of effective diagnostic or therapeutic biomarkers. The aim of this study was to find hub genes serving as biomarkers in gastric cancer diagnosis and therapy.

Methods

GSE66229 from Gene Expression Omnibus (GEO) was used as training set. Genes bearing the top 25% standard deviations among all the samples in training set were performed to systematic weighted gene co-expression network analysis (WGCNA) to find candidate genes. Then, hub genes were further screened by using the “least absolute shrinkage and selection operator” (LASSO) logistic regression. Finally, hub genes were validated in the GSE54129 dataset from GEO by supervised learning method artificial neural network (ANN) algorithm.

Results

Twelve modules with strong preservation were identified by using WGCNA methods in training set. Of which, five modules significantly related to gastric cancer were selected as clinically significant modules, and 713 candidate genes were identified from these five modules. Then, ADIPOQ, ARHGAP39, ATAD3A, C1orf95, CWH43, GRIK3, INHBA, RDH12, SCNN1G, SIGLEC11 and LYVE1 were screened as the hub genes. These hub genes successfully differentiated the tumor samples from the healthy tissues in an independent testing set through artificial neural network algorithm with the area under the receiver operating characteristic curve at 0.946.

Conclusions

These hub genes bearing diagnostic and therapeutic values, and our results may provide a novel prospect for the diagnosis and treatment of gastric cancer in the future.