Relationship between a 7-mRNA signature of the pancreatic adenocarcinoma microenvironment and patient prognosis (a STROBE-compliant article)

Comprehensive modular analyses of scar subtypes illuminate underlying molecular mechanisms and potential therapeutic targets

Pathological scarring resulting from traumas and wounds, such as hypertrophic scars and keloids, pose significant aesthetic, functional and psychological challenges. This study provides a comprehensive transcriptomic analysis of these conditions, aiming to illuminate underlying molecular mechanisms and potential therapeutic targets. We employed a co-expression and module analysis tool to identify significant gene clusters associated with distinct pathophysiological processes and mechanisms, notably lipid metabolism, sebum production, cellular energy metabolism and skin barrier function. This examination yielded critical insights into several skin conditions including folliculitis, skin fibrosis, fibrosarcoma and congenital ichthyosis. Particular attention was paid to Module Cluster (MCluster) 3, encompassing genes like BLK, TRPV1 and GABRD, all displaying high expression and potential implications in immune modulation. Preliminary immunohistochemistry validation supported these findings, showing elevated expression of these genes in non-fibrotic samples rich in immune activity. The complex interplay of different cell types in scar formation, such as fibroblasts, myofibroblasts, keratinocytes and mast cells, was also explored, revealing promising therapeutic strategies. This study underscores the promise of targeted gene therapy for pathological scars, paving the way for more personalised therapeutic approaches. The results necessitate further research to fully ascertain the roles of these identified genes and pathways in skin disease pathogenesis and potential therapeutics. Nonetheless, our work forms a strong foundation for a new era of personalised medicine for patients suffering from pathological scarring.

Immunological role and clinical prognostic significance of P2RY6 in lung adenocarcinoma: a multi-omics studies and single-cell sequencing analysis

BACKGROUND

There is increasing evidence that recombinant human P2Y purinoceptor 6 (P2RY6) may be involved in inflammatory responses. However, the role of P2RY6 in lung adenocarcinoma (LUAD) remains unknown.

METHODS

We used transcriptomic, genomic, single-cell transcriptomic, and methylation sequencing data from The Cancer Genome Atlas database to analyze the aberrant status and prognostic value of P2RY6 in a variety of tumors. The LUAD single-cell sequencing dataset was used to explore the effect of P2RY6 on the tumor microenvironment. Cell-type identification by estimating relative subsets of RNA transcripts (CIBERSORT) was used to quantify immune cells in the tumor microenvironment. We also analyzed the correlation of P2RY6 with immune checkpoints and immune regulation-related genes. The correlation of between tumor mutation burden (TMB), microsatellite instability (MSI), and P2RY6 expression was also analyzed simultaneously. Tissue microarray and immunohistochemistry were employed to assess the expression of P2RY6 in internal tumor samples.

RESULTS

Our findings indicate that P2RY6 exhibits significantly higher expression levels in various cancer tissues, particularly in LUAD. High expression of P2RY6 was closely associated with a poor prognosis for patients, and it plays a role in regulating immune-related pathways, such as cytokine-cytokine receptor interaction. Notably, P2RY6 expression is closely linked to the abundance of CD8 + T cells. Furthermore, we have developed a P2RY6-related inflammation prediction model that demonstrates promising results in predicting the prognosis of LUAD patients, with an AUC (area under the curve) value of 0.83. This performance is significantly better than the traditional TNM staging system. Through single-cell transcriptome sequencing analysis, we observed that high P2RY6 expression is associated with increased intercellular communication. Additionally, pathway enrichment analysis revealed that P2RY6 influences antigen presentation and processing pathways within the LUAD microenvironment.

CONCLUSIONS

This study suggests that P2RY6 would be a new target for immunotherapy in LUAD.

GBP5 Expression Predicted Prognosis of Immune Checkpoint Inhibitors in Small Cell Lung Cancer and Correlated with Tumor Immune Microenvironment

Background

The discovery and development of immune checkpoint inhibitors (ICIs) has significantly enhanced the arsenal of immunotherapy treatments available for cancer patients. The identification of biomarkers that are indicative of an individual's sensitivity to treatment with ICIs is useful for screening SCLC patients prior to commencement of any ICIs based immunotherapy. However, the relationship between GBP5 and the prognosis of SCLC immunotherapy is still unclear and requires further study.

Methods

We downloaded two SCLC datasets, namely the George-SCLC and Jiang-SCLC cohorts. We used the TIDE algorithm to predict the efficacy of immunotherapy for SCLC patients. The QuanTIseq, MCPcounter, and EPIC algorithms are used to calculate the proportions of immune cells in SCLC patients. Additionally, we retrospectively collected 35 SCLC samples from the first affiliated hospital of the Hengyang Medical school.

Results

Patients in each cohort were devided into two groups with high (GBP5-High) and low (GBP5-Low) expression of GBP5. In both cohorts, the GBP5-High population had a higher proportion of patients that responded well to immunotherapy (responders) (p < 0.05). In addition, both GBP5-High subgroups had significantly increased cytotoxicity, chemokines, antigen presenting, and TNF family related genes. We also determined that GBP5 was related to high-level infiltration of B cells, CD4+T cells, CD8+T cells and NK cells.

Conclusion

In this study, we found that GBP5 has the potential to be used as a biomarker of ICIs efficacy for SCLC patients. GBP5 is related to the quantity of inflammatory molecules, a high level of immune infiltration, and a highly activated immune response pathway.

Vertebrate Cell Differentiation, Evolution, and Diseases: The Vertebrate-Specific Developmental Potential Guardians VENTX/NANOG and POU5/OCT4 Enter the Stage

During vertebrate development, embryonic cells pass through a continuum of transitory pluripotent states that precede multi-lineage commitment and morphogenesis. Such states are referred to as “refractory/naïve” and “competent/formative” pluripotency. The molecular mechanisms maintaining refractory pluripotency or driving the transition to competent pluripotency, as well as the cues regulating multi-lineage commitment, are evolutionarily conserved. Vertebrate-specific “Developmental Potential Guardians” (vsDPGs; i.e., VENTX/NANOG, POU5/OCT4), together with MEK1 (MAP2K1), coordinate the pluripotency continuum, competence for multi-lineage commitment and morphogenesis in vivo. During neurulation, vsDPGs empower ectodermal cells of the neuro-epithelial border (NEB) with multipotency and ectomesenchyme potential through an “endogenous reprogramming” process, giving rise to the neural crest cells (NCCs). Furthermore, vsDPGs are expressed in undifferentiated-bipotent neuro-mesodermal progenitor cells (NMPs), which participate in posterior axis elongation and growth. Finally, vsDPGs are involved in carcinogenesis, whereby they confer selective advantage to cancer stem cells (CSCs) and therapeutic resistance. Intriguingly, the heterogenous distribution of vsDPGs in these cell types impact on cellular potential and features. Here, we summarize the findings about the role of vsDPGs during vertebrate development and their selective advantage in evolution. Our aim to present a holistic view regarding vsDPGs as facilitators of both cell plasticity/adaptability and morphological innovation/variation. Moreover, vsDPGs may also be at the heart of carcinogenesis by allowing malignant cells to escape from physiological constraints and surveillance mechanisms.

RAB39B as a Chemosensitivity-Related Biomarker for Diffuse Large B-Cell Lymphoma

Background: Diffuse large B-cell lymphoma (DLBCL) is the most common aggressive lymphoma with an increased tendency to relapse or refractoriness. RAB39B, a member of the Ras-oncogene superfamily, is associated with a variety of tumors. Nevertheless, the role of RAB39B in DLBCL is still unknown. This study aimed to identify the role of RAB39B in DLBCL using integrated bioinformatics analysis.

Methods: RAB39B expression data were examined using TIMER, UCSC, and GEO databases. The LinkedOmics database was used to study the genes and signaling pathways related to RAB39B expression. A Protein–protein interaction network was performed in STRING. TIMER was used to analyze the correlation between RAB39B and infiltrating immune cells. The correlation between RAB39B and m6A-related genes in DLBCL was analyzed using TCGA data. The RAB39B ceRNA network was constructed based on starBase and miRNet2.0 databases. Drug sensitivity information was obtained from the GSCA database.

Results: RAB39B was highly expressed in multiple tumors including DLBCL. The protein–protein interaction network showed enrichment of autophagy and RAS family proteins. Functional enrichment analysis of RAB39B co-expression genes revealed that RAB39B was closely related to DNA replication, protein synthesis, cytokine–cytokine receptor interaction, JAK-STAT signaling pathway, NF-kappa B signaling pathway, and autophagy. Immune infiltrate analysis showed that the amount of RAB39B was negatively correlated with iDC, Tem, and CD8 T-cell infiltration. CD4⁺ T cell and DC were negatively correlated with CNV of RAB39B. DLBCL cohort analysis found that RAB39B expression was related to 14 m6A modifier genes, including YTHDC1, YTHDC2, YTHDF1, YTHDF2, YTHDF3, RBMX, ZC3H13, METTL14, METTL3, RBM15, RBM15B, VIRMA, FTO, and ALKBH5. We constructed 14 possible ceRNA networks of RAB39B in DLBCL. The RAB39B expression was associated with decreased sensitivity of chemotherapy drugs such as dexamethasone, doxorubicin, etoposide, vincristine, and cytarabine and poor overall survival in DLBCL. In vitro experiments showed that RAB39B was associated with proliferation, apoptosis, and drug sensitivity of DLBCL cells.

Conclusion: RAB39B is abnormally elevated and related to drug resistance and poor OS in DLBCL, which may be due to its involvement in immune infiltration, m6A modification, and regulation by multiple non-coding RNAs. RAB39B may be used as an effective biomarker for the diagnosis and treatment of DLBCL.

A novel immune-related ceRNA network that predicts prognosis and immunotherapy response in lung adenocarcinoma

Background

The tumor microenvironment plays an important role in the progression and malignancy of lung adenocarcinoma and affects the immunotherapy response. There is increasing evidence that long non-coding RNAs (lncRNAs) as competing endogenous RNAs (ceRNAs) have significant functions in the development and treatment response of various kinds of cancer. This study aimed to explore the association between immune-related lncRNA-microRNA (miRNA)-messenger RNA (mRNA)-ceRNA networks, and the prognosis of and immunotherapy response in lung adenocarcinoma.

Methods

RNA-sequencing (RNA-seq) and miRNA-seq data from The Cancer Genome Atlas (TCGA) were used to evaluate the infiltration of immune cells in lung adenocarcinoma samples by undertaking a single-sample gene set enrichment analysis (ssGSEA) to divide the cells into high and low immune cell infiltration groups. The differentially expressed mRNA (DEmRNA) was further analyzed by a weighted gene co-expression network analysis (WGCNA), search tool for recurring instances of neighboring genes (STRING), and Cytoscape to select hub genes. The ceRNA network was constructed using Cytoscape. Additionally, survival analyses were conducted to screen out prognostic candidate genes.

Results

Seven thousand five hundred and thirty-eight mRNAs, 540 lncRNAs, and 138 miRNAs were found to be differentially expressed between the high and low immune cell infiltration groups. The two DEmRNA modules most significantly associated with immune cell infiltration were further analyzed, and four clusters, including 179 DEmRNAs, were selected based on Molecular Complex Detection (MCODE) scores. The selected DEmRNAs in the four clusters were mainly enriched in pathways involved in regulating the immune response. Ultimately, a ceRNA network of SNHG6-hsa-miR-30e-5p-CYSLTR1 was identified as being associated with the prognosis of and immunotherapy response in lung adenocarcinoma.

Conclusions

The present study extends understandings of immune-related lncRNA-miRNA-mRNA-ceRNA networks and identifies novel targets and a regulatory pathway for anti-tumor immunotherapy.

Identification of Prognostic Genes in the Tumor Microenvironment of Hepatocellular Carcinoma

Background: Hepatocellular carcinoma (HCC) is one of the most common malignant tumors in the world. The efficacy of immunotherapy usually depends on the interaction of immunomodulation in the tumor microenvironment (TME). This study aimed to explore the potential stromal-immune score-based prognostic genes related to immunotherapy in HCC through bioinformatics analysis.

Methods: ESTIMATE algorithm was applied to calculate the immune/stromal/Estimate scores and tumor purity of HCC using the Cancer Genome Atlas (TCGA) transcriptome data. Functional enrichment analysis of differentially expressed genes (DEGs) was analyzed by the Database for Annotation, Visualization, and Integrated Discovery database (DAVID). Univariate and multivariate Cox regression analysis and least absolute shrinkage and selection operator (LASSO) regression analysis were performed for prognostic gene screening. The expression and prognostic value of these genes were further verified by KM-plotter database and the Human Protein Atlas (HPA) database. The correlation of the selected genes and the immune cell infiltration were analyzed by single sample gene set enrichment analysis (ssGSEA) algorithm and Tumor Immune Estimation Resource (TIMER).

Results: Data analysis revealed that higher immune/stromal/Estimate scores were significantly associated with better survival benefits in HCC within 7 years, while the tumor purity showed a reverse trend. DEGs based on both immune and stromal scores primarily affected the cytokine–cytokine receptor interaction signaling pathway. Among the DEGs, three genes (CASKIN1, EMR3, and GBP5) were found most significantly associated with survival. Moreover, the expression levels of CASKIN1, EMR3, and GBP5 genes were significantly correlated with immune/stromal/Estimate scores or tumor purity and multiple immune cell infiltration. Among them, GBP5 genes were highly related to immune infiltration.

Conclusion: This study identified three key genes which were related to the TME and had prognostic significance in HCC, which may be promising markers for predicting immunotherapy outcomes.

P2Y receptors for extracellular nucleotides: Contributions to cancer progression and therapeutic implications

Purinergic receptors for extracellular nucleotides and nucleosides contribute to a vast array of cellular and tissue functions, including cell proliferation, intracellular and transmembrane ion flux, immunomodulation and thrombosis. In mammals, the purinergic receptor system is composed of G protein-coupled P1 receptors A₁, A_2A, A_2B and A₃ for extracellular adenosine, P2X1-7 receptors that are ATP-gated ion channels and G protein-coupled P2Y_{1,2,4,6,11,12,13 and 14} receptors for extracellular ATP, ADP, UTP, UDP and/or UDP-glucose. Recent studies have implicated specific P2Y receptor subtypes in numerous oncogenic processes, including cancer tumorigenesis, metastasis and chemotherapeutic drug resistance, where G protein-mediated signaling cascades modulate intracellular ion concentrations and activate downstream protein kinases, Src family kinases as well as numerous mitogen-activated protein kinases. We are honored to contribute to this special issue dedicated to the founder of the field of purinergic signaling, Dr. Geoffrey Burnstock, by reviewing the diverse roles of P2Y receptors in the initiation, progression and metastasis of specific cancers with an emphasis on pharmacological and genetic strategies employed to delineate cell-specific and P2Y receptor subtype-specific responses that have been investigated using in vitro and in vivo cancer models. We further highlight bioinformatic and empirical evidence on P2Y receptor expression in human clinical specimens and cover clinical perspectives where P2Y receptor-targeting interventions may have therapeutic relevance to cancer treatment.