Transcriptome Network Analysis Identifies CXCL13-CXCR5 Signaling Modules in the Prostate Tumor Immune Microenvironment

The role of proinflammatory cytokines and CXC chemokines (CXCL1-CXCL16) in the progression of prostate cancer: insights on their therapeutic management

Reproductive cancers are malignancies that develop in the reproductive organs. One of the leading cancers affecting the male reproductive system on a global scale is prostate cancer (PCa). The negative consequences of PCa metastases endure and are severe, significantly affecting mortality and life quality for those who are affected. The association between inflammation and PCa has captured interest for a while. Inflammatory cells, cytokines, CXC chemokines, signaling pathways, and other elements make up the tumor microenvironment (TME), which is characterized by inflammation. Inflammatory cytokines and CXC chemokines are especially crucial for PCa development and prognosis. Cytokines (interleukins) and CXC chemokines such as IL-1, IL-6, IL-7, IL-17, TGF-β, TNF-α, CXCL1-CXCL6, and CXCL8-CXCL16 are thought to be responsible for the pleiotropic effects of PCa, which include inflammation, progression, angiogenesis, leukocyte infiltration in advanced PCa, and therapeutic resistance. The inflammatory cytokine and CXC chemokines systems are also promising candidates for PCa suppression and immunotherapy. Therefore, the purpose of this work is to provide insight on how the spectra of inflammatory cytokines and CXC chemokines evolve as PCa develops and spreads. We also discussed recent developments in our awareness of the diverse molecular signaling pathways of these circulating cytokines and CXC chemokines, as well as their associated receptors, which may one day serve as PCa-targeted therapies. Moreover, the current status and potential of theranostic PCa therapies based on cytokines, CXC chemokines, and CXC receptors (CXCRs) are examined.

B cells and the coordination of immune checkpoint inhibitor response in patients with solid tumors

Immunotherapy profoundly changed the landscape of cancer therapy by providing long-lasting responses in subsets of patients and is now the standard of care in several solid tumor types. However, immunotherapy activity beyond conventional immune checkpoint inhibition is plateauing, and biomarkers are overall lacking to guide treatment selection. Most studies have focused on T cell engagement and response, but there is a growing evidence that B cells may be key players in the establishment of an organized immune response, notably through tertiary lymphoid structures. Mechanisms of B cell response include antibody-dependent cellular cytotoxicity and phagocytosis, promotion of CD4+ and CD8+ T cell activation, maintenance of antitumor immune memory. In several solid tumor types, higher levels of B cells, specific B cell subpopulations, or the presence of tertiary lymphoid structures have been associated with improved outcomes on immune checkpoint inhibitors. The fate of B cell subpopulations may be widely influenced by the cytokine milieu, with versatile roles for B-specific cytokines B cell activating factor and B cell attracting chemokine-1/CXCL13, and a master regulatory role for IL-10. Roles of B cell-specific immune checkpoints such as TIM-1 are emerging and could represent potential therapeutic targets. Overall, the expanding field of B cells in solid tumors of holds promise for the improvement of current immunotherapy strategies and patient selection.

TP63-TRIM29 axis regulates enhancer methylation and chromosomal instability in prostate cancer

BACKGROUND

Prostate adenocarcinoma (PRAD) is the second leading cause of cancer-related deaths in men. High variability in DNA methylation and a high rate of large genomic rearrangements are often observed in PRAD.

RESULTS

To investigate the reasons for such high variance, we integrated DNA methylation, RNA-seq, and copy number alterations datasets from The Cancer Genome Atlas (TCGA), focusing on PRAD, and employed weighted gene co-expression network analysis (WGCNA). Our results show that only single cluster of co-expressed genes is associated with genomic and epigenomic instability. Within this cluster, TP63 and TRIM29 are key transcription regulators and are downregulated in PRAD. We discovered that TP63 regulates the level of enhancer methylation in prostate basal epithelial cells. TRIM29 forms a complex with TP63 and together regulates the expression of genes specific to the prostate basal epithelium. In addition, TRIM29 binds DNA repair proteins and prevents the formation of the TMPRSS2:ERG gene fusion typically observed in PRAD.

CONCLUSION

Our study demonstrates that TRIM29 and TP63 are important regulators in maintaining the identity of the basal epithelium under physiological conditions. Furthermore, we uncover the role of TRIM29 in PRAD development.

Transcriptome Analysis Identifies Tumor Immune Microenvironment Signaling Networks Supporting Metastatic Castration-Resistant Prostate Cancer

Prostate cancer (PCa) is the second most common cause of cancer death in American men. Metastatic castration-resistant prostate cancer (mCRPC) is the most lethal form of PCa and preferentially metastasizes to the bones through incompletely understood molecular mechanisms. Herein, we processed RNA sequencing data from patients with mCRPC (n = 60) and identified 14 gene clusters (modules) highly correlated with mCRPC bone metastasis. We used a novel combination of weighted gene co-expression network analysis (WGCNA) and upstream regulator and gene ontology analyses of clinically annotated transcriptomes to identify the genes. The cyan module (M14) had the strongest positive correlation (0.81, p = 4 × 10-15) with mCRPC bone metastasis. It was associated with two significant biological pathways through KEGG enrichment analysis (parathyroid hormone synthesis, secretion, and action and protein digestion and absorption). In particular, we identified 10 hub genes (ALPL, PHEX, RUNX2, ENPP1, PHOSPHO1, PTH1R, COL11A1, COL24A1, COL22A1, and COL13A1) using cytoHubba of Cytoscape. We also found high gene expression for collagen formation, degradation, absorption, cell-signaling peptides, and bone regulation processes through Gene Ontology (GO) enrichment analysis.

MetaTiME integrates single-cell gene expression to characterize the meta-components of the tumor immune microenvironment

Recent advances in single-cell RNA sequencing have shown heterogeneous cell types and gene expression states in the non-cancerous cells in tumors. The integration of multiple scRNA-seq datasets across tumors can indicate common cell types and states in the tumor microenvironment (TME). We develop a data driven framework, MetaTiME, to overcome the limitations in resolution and consistency that result from manual labelling using known gene markers. Using millions of TME single cells, MetaTiME learns meta-components that encode independent components of gene expression observed across cancer types. The meta-components are biologically interpretable as cell types, cell states, and signaling activities. By projecting onto the MetaTiME space, we provide a tool to annotate cell states and signature continuums for TME scRNA-seq data. Leveraging epigenetics data, MetaTiME reveals critical transcriptional regulators for the cell states. Overall, MetaTiME learns data-driven meta-components that depict cellular states and gene regulators for tumor immunity and cancer immunotherapy.

Exosome-derived miRNAs regulate macrophage-colorectal cancer cell cross-talk during aggressive tumor development

Colorectal cancer is one of the leading causes of death worldwide. Its incidence and mortality have significantly increased during the past few years. Colorectal cancer cells cross-talk with other cells through exosomes in their tumor microenvironment. The miRNAs containing exosomes are responsible for tumor growth, invasion, and metastasis. Multiple studies have shown that exosomal miRNAs are key players in the crosstalk between cancerous, immune, and stromal cells during colorectal cancer development. They help in the establishment of the tumorigenic microenvironment by reprogramming macrophages towards a pro-tumorigenic phenotype. In this review, we discussed various exosomal miRNAs derived both from colorectal cancer cells and macrophages that promote or inhibit cancer aggression. We also discussed various miRNA-based therapeutic approaches to inhibit cancer progression.

Metabolic Networks: Weighted Gene Correlation Network Analysis

Weighted Gene Correlation Network Analysis (WGCNA) is used to build weighted gene networks representing direct interconnections among genes. This method is useful to identify gene modules associated with biological functions, revealing core functional differences between samples. Here, we describe a step-by-step methodology to build a WGCNA network associated with a phenotype of interest. The results can be visualized using Cytoscape and other available software or used as a basis for further functional enrichment analyses.

Data mining and network analysis reveals C-X-C chemokine receptor type 5 is involved in the pathophysiology of age-related macular degeneration

Our previous studies found that the C-X-C motif chemokine receptor 5 (CXCR5) loss leads to retinal pigment epithelium (RPE) dysfunction and AMD pathogenesis. The current study aimed to characterize the G protein-coupled receptor (GPCR) structure of CXCR5 and analyze its interactions with AMD-related risk genes. The sequence alignments, homology model of CXCR5 and structural assessment analysis were performed. Data and text mining were then performed to identify AMD-related risk genes and their interaction with CXCR5 using statistical and mathematical algorithms. Sequence alignment and phylogenetic tree analysis revealed that human CXCR5 was highly similar (85.4839%) to the rabbit. The least similarity (33.871%) was found to be in zebrafish compared to the other species. The CXCR5 model structural assessment and secondary structure analysis exhibited an excellent model. Network analysis revealed that IL10, TNF, ICAM1, CXCL1, CXCL8, APP, TLR4, SELL, C3, IL17A and CCR2 were the most connected genes CXCR5. These findings suggest that CXCR5 signaling may regulate the biological function of RPE and modulate AMD pathophysiology via GPCR signaling and interacting with identified AMD risk genes. In summary, the data presented here provide novel and crucial insights into the molecular mechanisms of CXCR5 involvement in AMD.Communicated by Ramaswamy H. Sarma.

Chemokines and cytokines: Axis and allies in prostate cancer pathogenesis

Chemokines are recognized as the major contributor to various tumorigenesis, tumor heterogeneity, and failures of current cancer therapies. The tumor microenvironment (TME) is enriched with chemokines and cytokines and plays a pivotal role in cancer progression. Chronic inflammation is also considered an instructive process of cancer progression, where chemokines are spatiotemporally secreted by malignant cells and leukocyte subtypes that initiate cell trafficking into the TME. In various cancers, prostate cancer (PCa) is reported as one of the leading cancers in the worldwide male population. The chemokines-mediated signaling pathways are intensively involved in PCa progression and metastasis. Emerging evidence suggests that chemokines and cytokines are responsible for the pleiotropic actions in cancer, including the growth, angiogenesis, endothelial mesenchymal transition, leukocyte infiltration, and hormone escape for advanced PCa and therapy resistance. Chemokine's system and immune cells represent a promising target to suppress tumorigenic environments and serve as potential therapy/immunotherapy for the PCa. In this review, an attempt has been made to shed light on the alteration of chemokine and cytokine profiles during PCa progression and metastasis. We also discussed the recent findings of the diverse molecular signaling of these circulating chemokines and their corresponding receptors that could become future targets for therapeutic management of PCa.

TIGIT is the central player in T-cell suppression associated with CAR T-cell relapse in mantle cell lymphoma

BACKGROUND

Chimeric antigen receptor (CAR) T-cell therapy using brexucabtagene autoleucel (BA) induces remission in many patients with mantle cell lymphoma (MCL), and BA is the only CAR T-cell therapy approved by the FDA for MCL. However, development of relapses to BA is recognized with poor patient outcomes. Multiple CAR T-cell therapies have been approved for other lymphomas and the resistance mechanisms have been investigated. However, the mechanisms underlying BA relapse in MCL have not been investigated and whether any previously reported resistance mechanisms apply to BA-relapsed patients with MCL is unknown.

METHODS

To interrogate BA resistance mechanisms in MCL, we performed single-cell RNA sequencing on 39 longitudinally collected samples from 15 BA-treated patients, and multiplex cytokine profiling on 80 serial samples from 20 patients.

RESULTS

We demonstrate that after BA relapse, the proportion of T cells, especially cytotoxic T cells (CTLs), decreased among non-tumor cells, while the proportion of myeloid cells correspondingly increased. TIGIT, LAG3, and CD96 were the predominant checkpoint molecules expressed on exhausted T cells and CTLs; only TIGIT was significantly increased after relapse. CTLs expanded during remission, and then contracted during relapse with upregulated TIGIT expression. Tumor cells also acquired TIGIT expression after relapse, leading to the enhanced interaction of tumor cell TIGIT with monocyte CD155/PVR. In myeloid cells, post-relapse HLA-II expression was reduced relative to pretreatment and during remission. Myeloid-derived suppressor cells (MDSCs) were enriched after relapse with elevated expression of activation markers, including CLU (clusterin) and VCAN (versican). Extracellular chemokines (CCL4, CXCL9, CXCL13), soluble checkpoint inhibitors (sPD-L1, sTIM3, s4-1BB), and soluble receptors (sIL-2R, sTNFRII) were decreased during remission but elevated after relapse.

CONCLUSIONS

Our data demonstrate that multiple tumor-intrinsic and -extrinsic factors are associated with T-cell suppression and BA relapse. Among these, TIGIT appears to be the central player given its elevated expression after BA relapse in not only CTLs but also MCL cells. The acquisition of TIGIT expression on tumor cells is MCL-specific and has not been reported in other CAR T-treated diseases. Together, our data suggest that co-targeting TIGIT may prevent CAR T relapses and thus promote long-term progression-free survival in MCL patients.

Serum Levels of CXCL13 Are an Independent Predictor of Survival Following Resection of Biliary Tract Cancer

Simple Summary

Biliary tract cancer (BTC) is a primary liver malignancy with poor outcome. The identification of the ideal surgical candidates is often challenging and stratification algorithms comprising the parameters of individual tumor biology are missing. Here, we investigated a potential role of circulating CXCL1, CXCL10 and CXCL13 in patients with resectable BTC as novel biomarkers and could show that elevated levels of CXCL13 both before and after tumor resection identified a subgroup of patients with significantly impaired outcomes following tumor resection. Thus, the present study supports a fundamental role of the CXC chemokine family in BTC and identifies circulating levels of CXCL13 as a previously unrecognized marker for predicting outcomes following the resection of BTC.

Abstract

Background: The prognosis of biliary tract cancer (BTC) has remained very poor. Although tumor resection represents a potentially curative therapy for selected patients, tumor recurrence is common, and 5-year survival rates have remained below 50%. As stratification algorithms comprising the parameters of individual tumor biology are missing, the identification of ideal patients for extensive tumor surgery is often challenging. The CXC chemokine family exerts decisive functions in cell–cell interactions and has only recently been associated with cancer, but little is known about their function in BTC. Here, we aim to evaluate a potential role of circulating CXCL1, CXCL10 and CXCL13 in patients with resectable BTC. Methods: Serum levels of CXCL1, CXCL10 and CXCL13 were measured by multiplex immunoassay in a cohort of 119 BTC patients undergoing tumor resection and 50 control samples. Results: Circulating levels of CXCL1, CXCL10 and CXCL13 were all significantly elevated in BTC patients compared to healthy controls and increased the diagnostic power of established tumor markers such as CA19-9 when used in combination. Importantly, elevated levels of CXCL13 both before and after tumor resection identified a subgroup of patients with significantly impaired outcomes following tumor resection. As such, BTC patients with initial CXCL13 levels above the ideal prognostic cut-off value (25.01 pg/mL) had a median overall survival (OS) of 290 days compared to 969 days for patients with low initial CXCL13 levels. The prognostic value of circulating CXCL13 was further confirmed by uni- and multivariate Cox regression analyses. Finally, the individual kinetics of CXCL13 before and after tumor resection were also indicative of patient outcomes. Conclusion: Our data support a fundamental role of the CXC chemokine family in BTC and identified circulating levels of CXCL13 as a previously unrecognized marker for predicting outcomes following the resection of BTC.

Chemokine receptors differentially expressed by race category and molecular subtype in the breast cancer TCGA cohort

Racial disparities in mortality due to metastasis remain significant among breast cancer patients. Chemokine receptors contribute to breast tumors and metastatic outcome. We explored for significant differences in chemokine receptor expression in breast tumors from Black, Asian, and White patients in The Cancer Genome Atlas. We show that despite sharing the same molecular subtype, expression of the chemokine receptors ACKR1, CCR3, CCR6, CCRL1, CCRL2, CXCR1, CXCR2, CXCR4, CXCR6, and CXC3CR1 was significantly different depending on racial group. For patients with triple negative breast cancer, CCR3 was higher in Black versus White and CCRL2 was higher in Asian versus White. In luminal A tumors, ACKR1 was lower in Asian versus White, CCR3 was higher in Black versus White, and CCR6 and CXC3CR1 were lower in Black versus White. In luminal B tumors, CCRL2 was lower in Black versus White, CXCR1 and CXC3CR1 were lower in Asian versus White, and CXCR2 was lower in Black and Asian versus White. In HER2 enriched tumors, CCR3 was higher in Black versus White and CXCR4 lower in Asian versus White. CCR3, CCR6, and CXCR6 associated with worse patient survival. These findings can inform improved treatment strategies to decrease racial disparities in breast cancer burden.

Role of chemokine systems in cancer and inflammatory diseases

Chemokines are a large family of small secreted proteins that have fundamental roles in organ development, normal physiology, and immune responses upon binding to their corresponding receptors. The primary functions of chemokines are to coordinate and recruit immune cells to and from tissues and to participate in regulating interactions between immune cells. In addition to the generally recognized antimicrobial immunity, the chemokine/chemokine receptor axis also exerts a tumorigenic function in many different cancer models and is involved in the formation of immunosuppressive and protective tumor microenvironment (TME), making them potential prognostic markers for various hematologic and solid tumors. In fact, apart from its vital role in tumors, almost all inflammatory diseases involve chemokines and their receptors in one way or another. Modulating the expression of chemokines and/or their corresponding receptors on tumor cells or immune cells provides the basis for the exploitation of new drugs for clinical evaluation in the treatment of related diseases. Here, we summarize recent advances of chemokine systems in protumor and antitumor immune responses and discuss the prevailing understanding of how the chemokine system operates in inflammatory diseases. In this review, we also emphatically highlight the complexity of the chemokine system and explore its potential to guide the treatment of cancer and inflammatory diseases.

JW. Multivariant Transcriptome Analysis Identifies Modules and Hub Genes Associated with Poor Outcomes in Newly Diagnosed Multiple Myeloma Patients

The molecular mechanisms underlying chemoresistance in some newly diagnosed multiple myeloma (MM) patients receiving standard therapies (lenalidomide, bortezomib, and dexamethasone) are poorly understood. Identifying clinically relevant gene networks associated with death due to MM may uncover novel mechanisms, drug targets, and prognostic biomarkers to improve the treatment of the disease. This study used data from the MMRF CoMMpass RNA-seq dataset (N = 270) for weighted gene co-expression network analysis (WGCNA), which identified 21 modules of co-expressed genes. Genes differentially expressed in patients with poor outcomes were assessed using two independent sample t-tests (dead and alive MM patients). The clinical performance of biomarker candidates was evaluated using overall survival via a log-rank Kaplan-Meier and ROC test. Four distinct modules (M10, M13, M15, and M20) were significantly correlated with MM vital status and differentially expressed between the dead (poor outcomes) and the alive MM patients within two years. The biological functions of modules positively correlated with death (M10, M13, and M20) were G-protein coupled receptor protein, cell-cell adhesion, cell cycle regulation genes, and cellular membrane fusion genes. In contrast, a negatively correlated module to MM mortality (M15) was the regulation of B-cell activation and lymphocyte differentiation. MM biomarkers CTAG2, MAGEA6, CCND2, NEK2, and E2F2 were co-expressed in positively correlated modules to MM vital status, which was associated with MM's lower overall survival.

CXCL13 in Cancer and Other Diseases: Biological Functions, Clinical Significance, and Therapeutic Opportunities

The development of cancer is a multistep and complex process involving interactions between tumor cells and the tumor microenvironment (TME). C-X-C chemokine ligand 13 (CXCL13) and its receptor, CXCR5, make crucial contributions to this process by triggering intracellular signaling cascades in malignant cells and modulating the sophisticated TME in an autocrine or paracrine fashion. The CXCL13/CXCR5 axis has a dominant role in B cell recruitment and tertiary lymphoid structure formation, which activate immune responses against some tumors. In most cancer types, the CXCL13/CXCR5 axis mediates pro-neoplastic immune reactions by recruiting suppressive immune cells into tumor tissues. Tobacco smoke and haze (smohaze) and the carcinogen benzo(a)pyrene induce the secretion of CXCL13 by lung epithelial cells, which contributes to environmental lung carcinogenesis. Interestingly, the knockout of CXCL13 inhibits benzo(a)pyrene-induced lung cancer and azoxymethane/dextran sodium sulfate-induced colorectal cancer in mice. Thus, a better understanding of the context-dependent functions of the CXCL13/CXCR5 axis in tumor tissue and the TME is required to design an efficient immune-based therapy. In this review, we summarize the molecular events and TME alterations caused by CXCL13/CXCR5 and briefly discuss the potentials of agents targeting this axis in different malignant tumors.

Identification of Gene Co-Expression Networks Associated with Consensus Molecular Subtype-1 of Colorectal Cancer

Simple Summary

Colorectal cancer (CRC) is a frequently lethal disease with heterogenous outcomes. Alterations in the Wnt signaling pathways have been shown to promote activation of signaling pathways such as MAPK and PI3K-Akt. Consensus molecular subtyping (CMS) provides a cohesive structure to classify the heterogeneity of CRC using gene expression analysis. CMS is categorized into four subtypes: CMS1, immune; CMS2, canonical; CMS3, metabolic; and CMS4, mesenchymal. Here, we identify co-expressed gene networks associated with CMS1. Our findings distinguish co-expressed gene networks that play a pivotal role in key features specific for CMS1, such as immune infiltration and activation. The co-expressed gene networks for CMS1 were significantly and positively correlated with the TNF, WNT, and ERK1 and ERK2 signaling pathways. This study highlights the relevance of CMS1 gene networks relating to oncogenic signaling cascades, cell activation, and positive regulation of immune responses, promoting CRC progressiveness.

Abstract

Colorectal cancer (CRC) is driven in part by dysregulated Wnt, Ras-Raf-MAPK, TGF-β, and PI3K-Akt signaling. The progression of CRC is also promoted by molecular alterations and heterogeneous—yet interconnected—gene mutations, chromosomal instability, transcriptomic subtypes, and immune signatures. Genomic alterations of CRC progression lead to changes in RNA expression, which support CRC metastasis. An RNA-based classification system used for CRC, known as consensus molecular subtyping (CMS), has four classes. CMS1 has the lowest survival after relapse of the four CRC CMS phenotypes. Here, we identify gene signatures and associated coding mRNAs that are co-expressed during CMS1 CRC progression. Using RNA-seq data from CRC primary tumor samples, acquired from The Cancer Genome Atlas (TCGA), we identified co-expression gene networks significantly correlated with CMS1 CRC progression. CXCL13, CXCR5, IL10, PIK3R5, PIK3AP1, CCL19, and other co-expressed genes were identified to be positively correlated with CMS1. The co-expressed eigengene networks for CMS1 were significantly and positively correlated with the TNF, WNT, and ERK1 and ERK2 signaling pathways, which together promote cell proliferation and survival. This network was also aligned with biological characteristics of CMS1 CRC, being positively correlated to right-sided tumors, microsatellite instability, chemokine-mediated signaling pathways, and immune responses. CMS1 also differentially expressed genes involved in PI3K-Akt signaling. Our findings reveal CRC gene networks related to oncogenic signaling cascades, cell activation, and positive regulation of immune responses distinguishing CMS1 from other CRC subtypes.

microRNA-425 loss mediates amyloid plaque microenvironment heterogeneity and promotes neurodegenerative pathologies

Different cellular and molecular changes underlie the pathogenesis of Alzheimer's disease (AD). Among these, neuron‐specific dysregulation is a necessary event for accumulation of classic pathologies including amyloid plaques. Here, we show that AD‐associated pathophysiology including neuronal cell death, inflammatory signaling, and endolysosomal dysfunction is spatially colocalized to amyloid plaques in regions with abnormal microRNA‐425 (miR‐425) levels and this change leads to focal brain microenvironment heterogeneity, that is, an amyloid plaque‐associated microenvironment (APAM). APAM consists of multiple specific neurodegenerative signature pathologies associated with senile plaques that contribute to the heterogeneity and complexity of AD. Remarkably, miR‐425, a neuronal‐specific regulator decreased in AD brain, maintains a normal spatial transcriptome within brain neurons. We tested the hypothesis that miR‐425 loss correlates with enhanced levels of mRNA targets downstream, supporting APAM and AD progression. A miR‐425‐deficient mouse model has enhanced APP amyloidogenic processing, neuroinflammation, neuron loss, and cognitive impairment. In the APP/PS1 mouse model, intervening with miR‐425 supplementation ameliorated APAM changes and memory deficits. This study reveals a novel mechanism of dysregulation of spatial transcriptomic changes in AD brain, identifying a probable neuronal‐specific microRNA regulator capable of staving off amyloid pathogenesis. Moreover, our findings provide new insights for developing AD treatment strategies with miRNA oligonucleotide(s).

Characteristics of B lymphocyte infiltration in HPV+ head and neck squamous cell carcinoma

Human papillomavirus (HPV) is an important etiological factor of head and neck squamous cell carcinoma (HNSCC). HPV⁺ HNSCC patients usually have a better prognosis, which probably results from the higher infiltration of B lymphocytes. This study was purposed to detect the infiltration of B lymphocyte subsets and the correlation between B lymphocyte subsets and the prognosis in HPV‐related HNSCC. In this study, 124 HPV⁺ and 513 HPV⁻ HNSCC samples were obtained from the Gene Expression Omnibus (GEO) database and The Cancer Genome Atlas (TCGA) database for transcriptomic analysis. Infiltration of B lymphocytes subsets was detected with 7 HPV⁺ HNSCC and 13 HPV⁻ HNSCC tissues through immunohistochemistry and immunofluorescence. One HPV⁻ HNSCC sample was detected with single‐cell sequencing for chemokine analysis. In the results, the infiltration of plasma cells (CD19⁺CD38⁺) and memory B cells (MS4A1⁺CD27⁺) was higher in HPV⁺ HNSCC samples. High infiltration of plasma cells and memory B cells was related to a better prognosis. High density of B lymphocytes was positively correlated with high CXCL13 production mainly from CD4⁺ T lymphocytes in HNSCC. These results indicated that a high density of plasma cells and memory B cells could predict excellent prognosis. CD4⁺ T lymphocytes might affect B lymphocytes and their subsets through the CXCL13/CXCR5 axis in HNSCC.

Bacillus amyloliquefaciens TL Downregulates the Ileal Expression of Genes Involved in Immune Responses in Broiler Chickens to Improve Growth Performance

Bacillus amyloliquefaciens TL promotes broiler chicken performance by improving nutrient absorption and utilization and reducing intestinal inflammation. In this study, RNA-sequencing (RNA-seq)-based transcriptomes of ileal tissues collected from probiotic-fed and control broiler chickens were analyzed to elucidate the effects of the probiotic B. amyloliquefaciens TL, as a feed additive, on the gut immune function. In total, 475 genes were significantly differentially expressed between the ileum of probiotic-fed and control birds. The expression of genes encoding pyruvate kinase, prothymosin-α, and heat stress proteins was high in the ileum of probiotic-fed birds (FPKM > 500), but not in the control group. The gene ontology functional enrichment and pathway enrichment analyses revealed that the uniquely expressed genes in the control group were mostly involved in immune responses, whereas those in the probiotic group were involved in fibroblast growth factor receptor signaling pathways and positive regulation of cell proliferation. Bacillus amyloliquefaciens TL downregulated the expression of certain proinflammatory factors and affected the cytokine–cytokine receptor interaction pathway. Furthermore, B. amyloliquefaciens TL in broiler diets altered the expression of genes involved in immune functions in the ileum. Thus, it might contribute to improved broiler growth by regulating the immune system and reducing intestinal damage in broilers.

Hub gene identification and prognostic model construction for isocitrate dehydrogenase mutation in glioma

•

We identified ten hub genes which were driving IDH status in GBM and LGG.

•

We constructed a prognostic model for IDH-mutant patients.

•

Our findings have important clinical implications for accurate treatment in glioma.

Our study attempted to identify hub genes related to isocitrate dehydrogenase (IDH) mutation in glioma and develop a prognostic model for IDH-mutant glioma patients. In a first step, ten hub genes significantly associated with the IDH status were identified by weighted gene coexpression analysis (WGCNA). The functional enrichment analysis demonstrated that the most enriched terms of these hub genes were cadherin binding and glutathione metabolism. Three of these hub genes were significantly linked with the survival of glioma patients. 328 samples of IDH-mutant glioma were separated into two datasets: a training set (N = 228) and a test set (N = 100). Based on the training set, we identified two IDH-mutant subtypes with significantly different pathological features by using consensus clustering. A 31 gene-signature was identified by the least absolute shrinkage and selection operator (LASSO) algorithm and used for establishing a differential prognostic model for IDH-mutant patients. In addition, the test set was employed for validating the prognostic model, and the model was proven to be of high value in classifying prognostic information of samples. The functional annotation revealed that the genes related to the model were mainly enriched in nuclear division, DNA replication, and cell cycle. Collectively, this study provided novel insights into the molecular mechanism of IDH mutation in glioma, and constructed a prognostic model which can be effective for predicting prognosis of glioma patients with IDH-mutation, which might promote the development of IDH target agents in glioma therapies and contribute to accurate prognostication and management in IDH-mutant glioma patients.

Tumor-derived exosomal miR-934 induces macrophage M2 polarization to promote liver metastasis of colorectal cancer

Background

Mounting evidence has demonstrated the vital importance of tumor-associated macrophages (TAMs) and exosomes in the formation of the premetastatic niche. However, the molecular mechanisms by which tumor-derived exosomal miRNAs interact with TAMs underlying premetastatic niche formation and colorectal cancer liver metastasis (CRLM) remain largely unknown.

Methods

Transmission electron microscopy and differential ultracentrifugation were used to verify the existence of exosomes. In vivo and in vitro assays were used to identify roles of exosomal miR-934. RNA pull-down assay, dual-luciferase reporter assay, etc. were applied to clarify the mechanism of exosomal miR-934 regulated the crosstalk between CRC cells and M2 macrophages.

Results

In the present study, we first demonstrated the aberrant overexpression of miR-934 in colorectal cancer (CRC), especially in CRLM, and its correlation with the poor prognosis of CRC patients. Then, we verified that CRC cell-derived exosomal miR-934 induced M2 macrophage polarization by downregulating PTEN expression and activating the PI3K/AKT signaling pathway. Moreover, we revealed that hnRNPA2B1 mediated miR-934 packaging into exosomes of CRC cells and then transferred exosomal miR-934 into macrophages. Interestingly, polarized M2 macrophages could induce premetastatic niche formation and promote CRLM by secreting CXCL13, which activated a CXCL13/CXCR5/NFκB/p65/miR-934 positive feedback loop in CRC cells.

Conclusions

These findings indicate that tumor-derived exosomal miR-934 can promote CRLM by regulating the crosstalk between CRC cells and TAMs. These findings reveal a tumor and TAM interaction in the metastatic microenvironment mediated by tumor-derived exosomes that affects CRLM. The present study also provides a theoretical basis for secondary liver cancer.

Database Mining of Genes of Prognostic Value for the Prostate Adenocarcinoma Microenvironment Using the Cancer Gene Atlas

Background

Prostate adenocarcinoma (PRAD) is a common malignant tumor in elderly men. Our research uses The Cancer Gene Atlas (TCGA) database to find potential related genes for predicting the prognosis of patients with PRAD.

Methods

We downloaded gene expression profiles and clinical sample information from TCGA for 490 patients with PRAD (patient age: 41-78 years). We calculated stromal and immune scores using the ESTIMATE algorithm to predict the level of stromal and immune cell infiltration. We categorized patients with PRAD in TCGA into high and low score arrays according to their median immune/stromal scores and identified differentially expressed genes (DEGs) that were significantly correlated with the prognosis of PRAD. Then, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were performed. The association between DEGs and overall survival was investigated by weighted Kaplan–Meier survival analysis and multivariate analysis. Furthermore, the protein-protein interaction network (PPI) of DEGs was constructed using the STRING tool. Finally, the hub genes were identified by analyzing the degree of association of PPI networks.

Results

We found that 8 individual DEGs, C6, S100A12, MLC1, PAX5, C7, FAM162B, CAMK1G, and TCEAL5, were significantly predictive of favorable overall survival and one DEG, EPYC, was associated with poor overall survival. GO and KEGG pathway analyses revealed that the DEGs were associated with immune responses. Moreover, 30 hub genes were obtained using the PPI network of DEGs: ITGAM, CD4, CD3E, IL-10, LCP2, ITGB2, ZAP-70, C3, CCL19, CXCL13, CXCL9, BTK, CCL21, CD247, CD28, CD3D, FCER1G, PTPRC, TYROBP, CCR5, ITK, CCL13, CCR1, CCR2, CD79B, CYBB, IL2RG, JAK3, PLCG2, and CD19. These prominent nodes had the most associations with other genes, indicating that they might play crucial roles in the prognosis of PRAD.

Conclusions

We extracted a list of genes associated with the prostate adenocarcinoma microenvironment, which might contribute to the prediction and interpretation of PRAD prognosis.