Integrating genetic and gene co-expression analysis identifies gene networks involved in alcohol and stress responses

Although the link between stress and alcohol is well recognized, the underlying mechanisms of how they interplay at the molecular level remain unclear. The purpose of this study is to identify molecular networks underlying the effects of alcohol and stress responses, as well as their interaction on anxiety behaviors in the hippocampus of mice using a systems genetics approach. Here, we applied a gene co-expression network approach to transcriptomes of 41 BXD mouse strains under four conditions: stress, alcohol, stress-induced alcohol and control. The co-expression analysis identified 14 modules and characterized four expression patterns across the four conditions. The four expression patterns include up-regulation in no restraint stress and given an ethanol injection (NOE) but restoration in restraint stress followed by an ethanol injection (RSE; pattern 1), down-regulation in NOE but rescue in RSE (pattern 2), up-regulation in both restraint stress followed by a saline injection (RSS) and NOE, and further amplification in RSE (pattern 3), and up-regulation in RSS but reduction in both NOE and RSE (pattern 4). We further identified four functional subnetworks by superimposing protein-protein interactions (PPIs) to the 14 co-expression modules, including g-aminobutyric acid receptor (GABA) signaling, glutamate signaling, neuropeptide signaling, cAMP-dependent signaling. We further performed module specificity analysis to identify modules that are specific to stress, alcohol, or stress-induced alcohol responses. Finally, we conducted causality analysis to link genetic variation to these identified modules, and anxiety behaviors after stress and alcohol treatments. This study underscores the importance of integrative analysis and offers new insights into the molecular networks underlying stress and alcohol responses

SPP1 is associated with adverse prognosis and predicts immunotherapy efficacy in penile cancer

Abstract Background The effect of SPP1 in squamous cell carcinoma of the penis (PSCC) remained unknown. We attempted to clarify the function of the SPP1 gene in PSCC. Method Eight paired penile cancer specimens (including penile cancer tissue, paracancerous tissue, and positive lymph node tissue) subjected to whole transcriptome sequencing were analysed to identify differentially expressed genes. We used immunohistochemistry to detect the expression of SPP1 protein and immune cell related proteins in penile cancer tissue. Then, we performed weighted gene coexpression network analysis (WGCNA) to identify the genes related to SPP1 in penile cancer tissue and positive lymph node tissue. Based on the GSE57955 dataset, the CIBERSORT and ssGSEA algorithms were carried out to investigate the immune environment of PSCC. GSVA analysis was conducted to identify the signaling pathways related to SPP1 subgroups. Enzyme-linked immunosorbent assay (ELISA) method was adopted to detect SPP1 level in the serum of 60 patients with penile cancer. Results Differential analysis indicated that SPP1 was the most differentially upregulated gene in both penile cancer tissues and positive lymph node tissues. Survival analysis suggested that the prognosis of the low-SPP1 group was significantly poorer than that of the high-SPP1 group. Subsequently, immune-related bioinformatics showed that SPP1 was significantly associated with B cells, CD8 + T cells, CD4 + T cells, macrophages, helper T cells, neutrophils and dendritic cells. The immunohistochemical results showed that the high-SPP1 group was characterized by relatively high expression of CD16 and relatively low expression of CD4. GSVA analysis indicated that high-SPP1 group was significantly associated with immune-related pathways such as PD-L1 expression and the PD-1 checkpoint pathway in cancer and the TNF signaling pathway. ELISA demonstrated that the serum level of SPP1 in patients with positive lymph node metastasis of penile cancer was significantly higher than that in patients with negative lymph node metastasis of penile cancer. Conclusion Our study shows that the SPP1 gene might be an effective biomarker for predicting the prognosis and the efficacy of immunotherapy in PSCC patients

Image1_Immune-related gene risk score predicting the effect of immunotherapy and prognosis in bladder cancer patients.TIF

Background: Immune checkpoint inhibitor therapy has changed the treatment model of metastatic bladder cancer. However, only approximately 20% of patients benefit from this therapy, and robust biomarkers to predict the effect of immunotherapy are still lacking. In this study, we aimed to investigate whether immune-related genes could be indicators for the prognosis of bladder cancer patients and the effect of immunotherapy.Methods: Based on bladder cancer dataset from the Cancer Genome Atlas (TCGA) and GSE48075, 22 immune microenvironment-related cells were identified by CIBERSORT. After performing a series of bioinformatic and machine learning approaches, we identified distinct tumor microenvironment clusters and three bladder cancer specific immune-related genes (EGFR, OAS1 and MST1R). Then, we constructed immune-related gene risk score (IRGRS) by using the Cox regression method and validated it with the IMvigor210 dataset.Results: IRGRS-high patients had a worse overall survival than IRGRS-low patients, which was consistent with the result in the IMvigor210 dataset. Comprehensive analysis shows that patients with high IRGRS scores are mainly enriched in basal/squamous type (Ba/Sq), and tumor metabolism-related pathways are more Active, with higher TP53 and RB1 gene mutation rates, lower CD4+/CD8+ T cell infiltration, higher M0 macrophage infiltration, and lower immunotherapy efficacy. In contrast, Patients with low IRGRS scores are mainly enriched in the luminal papillary type (LumP), which is associated with the activation of IL-17 and TNF signaling pathways, higher mutation rates of FGFR3 and CDKN1A genes, higher CD4+/CD8+ T cell infiltration content, and The level of M0 macrophage infiltration was relatively low, and the immunotherapy was more probably effective.Conclusion: Our study constructed an IRGRS for bladder cancer and clarified the immune and molecular characteristics of IRGRS-defined subgroups of bladder cancer to investigate the association between IRGRS and its potential implications for prognosis and immunotherapy.</p

RAB20 Promotes Proliferation via G2/M Phase through the Chk1/cdc25c/cdc2-cyclinB1 Pathway in Penile Squamous Cell Carcinoma

RAB20, a member of the RAS GTPase oncogene family, is overexpressed in several cancers with poor outcomes, promoting tumorigenesis and inducing genomic instability. Here, we performed comprehensive genomic sequencing on eight penile squamous cell carcinoma (PSCC) and normal tissue pairs and found that RAB20 was upregulated in tumors, especially in metastatic lymph nodes. RAB20 overexpression in tumors was further verified by qPCR, Western blotting, and immunohistochemistry of our newly established PSCC cell lines and paired tissues. The clinical significance of RAB20 was validated in 259 PSCC patients, the largest cohort to date, and high RAB20 expression positively correlated with the T, N, M status, extranodal extension, and clinical stage (all p p = 0.011, HR = 2.090; 95% Cl: 1.183–4.692), and PSCC patients with high RAB20 expression experienced shorter 5-year cancer-specific survival times (p < 0.001). Furthermore, tumorigenesis assays demonstrated that RAB20 knockdown inhibited cell proliferation, migration, and colony formation in vitro and tumor growth in vivo. RAB20 depletion also induced PSCC cell cycle arrest at G2/M by increasing Chk1 expression and promoting cdc25c phosphorylation to reduce cdc2-cyclinB1 complex formation. Our study revealed an oncogenic role for RAB20 in promoting PSCC cell proliferation at the G2/M phase via the Chk1/cdc25c/cdc2-cyclinB1 pathway. Thus, RAB20 could be a promising prognostic biomarker of advanced PSCC with poor patient survival outcomes and could be a potential therapeutic target

Table4_Immune-related gene risk score predicting the effect of immunotherapy and prognosis in bladder cancer patients.XLS

Background: Immune checkpoint inhibitor therapy has changed the treatment model of metastatic bladder cancer. However, only approximately 20% of patients benefit from this therapy, and robust biomarkers to predict the effect of immunotherapy are still lacking. In this study, we aimed to investigate whether immune-related genes could be indicators for the prognosis of bladder cancer patients and the effect of immunotherapy.Methods: Based on bladder cancer dataset from the Cancer Genome Atlas (TCGA) and GSE48075, 22 immune microenvironment-related cells were identified by CIBERSORT. After performing a series of bioinformatic and machine learning approaches, we identified distinct tumor microenvironment clusters and three bladder cancer specific immune-related genes (EGFR, OAS1 and MST1R). Then, we constructed immune-related gene risk score (IRGRS) by using the Cox regression method and validated it with the IMvigor210 dataset.Results: IRGRS-high patients had a worse overall survival than IRGRS-low patients, which was consistent with the result in the IMvigor210 dataset. Comprehensive analysis shows that patients with high IRGRS scores are mainly enriched in basal/squamous type (Ba/Sq), and tumor metabolism-related pathways are more Active, with higher TP53 and RB1 gene mutation rates, lower CD4+/CD8+ T cell infiltration, higher M0 macrophage infiltration, and lower immunotherapy efficacy. In contrast, Patients with low IRGRS scores are mainly enriched in the luminal papillary type (LumP), which is associated with the activation of IL-17 and TNF signaling pathways, higher mutation rates of FGFR3 and CDKN1A genes, higher CD4+/CD8+ T cell infiltration content, and The level of M0 macrophage infiltration was relatively low, and the immunotherapy was more probably effective.Conclusion: Our study constructed an IRGRS for bladder cancer and clarified the immune and molecular characteristics of IRGRS-defined subgroups of bladder cancer to investigate the association between IRGRS and its potential implications for prognosis and immunotherapy.</p

Table1_Immune-related gene risk score predicting the effect of immunotherapy and prognosis in bladder cancer patients.XLS

Background: Immune checkpoint inhibitor therapy has changed the treatment model of metastatic bladder cancer. However, only approximately 20% of patients benefit from this therapy, and robust biomarkers to predict the effect of immunotherapy are still lacking. In this study, we aimed to investigate whether immune-related genes could be indicators for the prognosis of bladder cancer patients and the effect of immunotherapy.Methods: Based on bladder cancer dataset from the Cancer Genome Atlas (TCGA) and GSE48075, 22 immune microenvironment-related cells were identified by CIBERSORT. After performing a series of bioinformatic and machine learning approaches, we identified distinct tumor microenvironment clusters and three bladder cancer specific immune-related genes (EGFR, OAS1 and MST1R). Then, we constructed immune-related gene risk score (IRGRS) by using the Cox regression method and validated it with the IMvigor210 dataset.Results: IRGRS-high patients had a worse overall survival than IRGRS-low patients, which was consistent with the result in the IMvigor210 dataset. Comprehensive analysis shows that patients with high IRGRS scores are mainly enriched in basal/squamous type (Ba/Sq), and tumor metabolism-related pathways are more Active, with higher TP53 and RB1 gene mutation rates, lower CD4+/CD8+ T cell infiltration, higher M0 macrophage infiltration, and lower immunotherapy efficacy. In contrast, Patients with low IRGRS scores are mainly enriched in the luminal papillary type (LumP), which is associated with the activation of IL-17 and TNF signaling pathways, higher mutation rates of FGFR3 and CDKN1A genes, higher CD4+/CD8+ T cell infiltration content, and The level of M0 macrophage infiltration was relatively low, and the immunotherapy was more probably effective.Conclusion: Our study constructed an IRGRS for bladder cancer and clarified the immune and molecular characteristics of IRGRS-defined subgroups of bladder cancer to investigate the association between IRGRS and its potential implications for prognosis and immunotherapy.</p

Image2_Immune-related gene risk score predicting the effect of immunotherapy and prognosis in bladder cancer patients.TIF

Background: Immune checkpoint inhibitor therapy has changed the treatment model of metastatic bladder cancer. However, only approximately 20% of patients benefit from this therapy, and robust biomarkers to predict the effect of immunotherapy are still lacking. In this study, we aimed to investigate whether immune-related genes could be indicators for the prognosis of bladder cancer patients and the effect of immunotherapy.Methods: Based on bladder cancer dataset from the Cancer Genome Atlas (TCGA) and GSE48075, 22 immune microenvironment-related cells were identified by CIBERSORT. After performing a series of bioinformatic and machine learning approaches, we identified distinct tumor microenvironment clusters and three bladder cancer specific immune-related genes (EGFR, OAS1 and MST1R). Then, we constructed immune-related gene risk score (IRGRS) by using the Cox regression method and validated it with the IMvigor210 dataset.Results: IRGRS-high patients had a worse overall survival than IRGRS-low patients, which was consistent with the result in the IMvigor210 dataset. Comprehensive analysis shows that patients with high IRGRS scores are mainly enriched in basal/squamous type (Ba/Sq), and tumor metabolism-related pathways are more Active, with higher TP53 and RB1 gene mutation rates, lower CD4+/CD8+ T cell infiltration, higher M0 macrophage infiltration, and lower immunotherapy efficacy. In contrast, Patients with low IRGRS scores are mainly enriched in the luminal papillary type (LumP), which is associated with the activation of IL-17 and TNF signaling pathways, higher mutation rates of FGFR3 and CDKN1A genes, higher CD4+/CD8+ T cell infiltration content, and The level of M0 macrophage infiltration was relatively low, and the immunotherapy was more probably effective.Conclusion: Our study constructed an IRGRS for bladder cancer and clarified the immune and molecular characteristics of IRGRS-defined subgroups of bladder cancer to investigate the association between IRGRS and its potential implications for prognosis and immunotherapy.</p

Volatiles and Transcriptome Profiling Revealed the Formation of &lsquo;Taro-like&rsquo; Aroma in the Leaf of Pumpkin (Cucurbita moschata)

&lsquo;Taro-like&rsquo; odor is an important economic trait of pumpkin species. The metabolic and molecular bases of this aromatic trait remain largely unexplored. Therefore, in this study, gas chromatography-mass spectrometry, GC-Olfactometry, and RNA-seq technology were used to illuminate the differential volatile compounds, the key volatile compounds, and differentially expressed genes (DEGs) in leaves from two pumpkin samples. Eight volatile compounds, including (E)-2-nonenal, 3-octanol, 2-ethyl-1-hexanol, 1-nonanol, &alpha;-terpineol, 2,3-pentanedione, caryophyllene, and 2-acetyl-1-pyrroline, were only detected in the sample with &lsquo;taro-like&rsquo; aroma. Moreover, the variable importance in projection scores of all the above eight volatile compounds were &gt;1.0 using PLS-DA analysis. The compounds 2-acetyl-1-pyrroline, 3-octanol, 1-nonanol, and (E)-3,7-dimethyl-2,6-octadienal were identified as the key contributors using GC-Olfactometry analysis. It was determined that 2-acetyl-1-pyrroline might play a significant role in &lsquo;taro-like&rsquo; aroma. Furthermore, most of the differential volatile compounds were derived from fatty acids, and the DEGs were also involved in the pathways related to degradation, metabolism, and biosynthesis of fatty acids. Moreover, five genes involved in the accumulation of 2-acetyl-1-pyrroline showed differential expression, and their expression trends were consistent with 2-acetyl-1-pyrroline. This study offers the basis for further studies on the mechanism of &lsquo;taro-like&rsquo; aroma in pumpkins