The murine Microenvironment Cell Population counter method to estimate abundance of tissue-infiltrating immune and stromal cell populations in murine samples using gene expression

Overexpression of miR-96 leads to retinal degeneration in mice

Double knockout of miR-183 and miR-96 results in retinal degeneration in mice; however, single knockout of miR-96 leads to developmental delay but not substantial retinal degeneration. To further explore the role of miR-96, we overexpressed this miRNA in mouse retinas. Interestingly, we found that overexpression of miR-96 at a safe dose results in retinal degeneration in the mouse retina. The retinal photoreceptors dramatically degenerated in the miR-96-overexpressing group, as shown by OCT, ERG and cryosectioning at one month after subretinal injection. Degenerative features such as TUNEL signals and reactive gliosis were observed in the miR-96-overexpressing retina. RNA-seq data revealed that immune responses and microglial activation occurred in the degenerating retina. Further qRT‒PCR and immunostaining experiments verified the microglial activation. Moreover, the number of microglia in the miR-96-overexpressing retinas was significantly increased. Our findings demonstrate that appropriate miR-96 expression is required for mouse retinal homeostasis.

Enhancing antitumor immunity and achieving tumor eradication with IL11RA mRNA immunotherapy

Current methods for delivering genes to target tumors face significant challenges, including off-target effects and immune responses against delivery vectors. In this study, we developed a novel approach using messenger RNA (mRNA) to encode IL11RA for local immunotherapy, aiming to harness the immune system to combat tumors. Our research uncovered a compelling correlation between IL11RA expression and CD8 + T cell levels across multiple tumor types, with elevated IL11RA expression correlating with improved overall survival. Examination of the Pan-Cancer Atlas dataset showed a significant reduction in IL11RA expression in various cancer types compared to normal tissue, raising questions about its potential role in tumorigenesis. To achieve efficient in vivo expression of IL11RA, we synthesized two mRNA sequences mimicking the wild-type protein. These mRNA sequences were formulated and capped to ensure effective delivery, resulting in robust expression within tumor sites. Our investigation into IL11RA mRNA therapy demonstrated its effectiveness in controlling tumor growth when administered both intratumorally and intravenously in mouse models. Additionally, IL11RA mRNA treatment significantly stimulated the expansion of CD8 + T cells within tumors, draining lymph nodes, and the spleen. Transcriptome analysis revealed distinct transcriptional patterns associated with T cell functions. Using multiple deconvolution algorithms, we found substantial infiltration of CD8 + T cells following IL11RA mRNA treatment, highlighting its immunomodulatory effects within the tumor microenvironment. In conclusion, IL11RA mRNA therapy presents a promising strategy for tumor regression with potential immunomodulatory effects and clinical implications for improved survival outcomes.

Obesogenic High-Fat Diet and MYC Cooperate to Promote Lactate Accumulation and Tumor Microenvironment Remodeling in Prostate Cancer

Cancer cells exhibit metabolic plasticity to meet oncogene-driven dependencies while coping with nutrient availability. A better understanding of how systemic metabolism impacts the accumulation of metabolites that reprogram the tumor microenvironment (TME) and drive cancer could facilitate development of precision nutrition approaches. Using the Hi-MYC prostate cancer mouse model, we demonstrated that an obesogenic high-fat diet (HFD) rich in saturated fats accelerates the development of c-MYC-driven invasive prostate cancer through metabolic rewiring. Although c-MYC modulated key metabolic pathways, interaction with an obesogenic HFD was necessary to induce glycolysis and lactate accumulation in tumors. These metabolic changes were associated with augmented infiltration of CD206+ and PD-L1+ tumor-associated macrophages (TAM) and FOXP3+ regulatory T cells, as well as with the activation of transcriptional programs linked to disease progression and therapy resistance. Lactate itself also stimulated neoangiogenesis and prostate cancer cell migration, which were significantly reduced following treatment with the lactate dehydrogenase inhibitor FX11. In patients with prostate cancer, high saturated fat intake and increased body mass index were associated with tumor glycolytic features that promote the infiltration of M2-like TAMs. Finally, upregulation of lactate dehydrogenase, indicative of a lactagenic phenotype, was associated with a shorter time to biochemical recurrence in independent clinical cohorts. This work identifies cooperation between genetic drivers and systemic metabolism to hijack the TME and promote prostate cancer progression through oncometabolite accumulation. This sets the stage for the assessment of lactate as a prognostic biomarker and supports strategies of dietary intervention and direct lactagenesis blockade in treating advanced prostate cancer.

SIGNIFICANCE

Lactate accumulation driven by high-fat diet and MYC reprograms the tumor microenvironment and promotes prostate cancer progression, supporting the potential of lactate as a biomarker and therapeutic target in prostate cancer. See related commentary by Frigo, p. 1742.

A computational pipeline for identifying gene targets and signalling pathways in cancer cells to improve lymphocyte infiltration and immune checkpoint therapy efficacy

BACKGROUND

Tumour-infiltrating lymphocytes (TILs) are crucial for effective immune checkpoint blockade (ICB) therapy in solid tumours. However, ∼70% of these tumours exhibit poor lymphocyte infiltration, rendering ICB therapies less effective.

METHODS

We developed a bioinformatics pipeline integrating multiple previously unconsidered factors or datasets, including tumour cell immune-related pathways, copy number variation (CNV), and single tumour cell sequencing data, as well as tumour mRNA-seq data and patient survival data, to identify targets that can potentially improve T cell infiltration and enhance ICB efficacy. Furthermore, we conducted wet-lab experiments and successfully validated one of the top-identified genes.

FINDINGS

We applied this pipeline in solid tumours of the Cancer Genome Atlas (TCGA) and identified a set of genes in 18 cancer types that might potentially improve lymphocyte infiltration and ICB efficacy, providing a valuable drug target resource to be further explored. Importantly, we experimentally validated SUN1, which had not been linked to T cell infiltration and ICB therapy previously, but was one of the top-identified gene targets among 3 cancer types based on the pipeline, in a mouse colon cancer syngeneic model. We showed that Sun1 KO could significantly enhance antigen presentation, increase T-cell infiltration, and improve anti-PD1 treatment efficacy. Moreover, with a single-cell multiome analysis, we identified subgene regulatory networks (sub-GRNs) showing Stat proteins play important roles in enhancing the immune-related pathways in Sun1-KO cancer cells.

INTERPRETATION

This study not only established a computational pipeline for discovering new gene targets and signalling pathways in cancer cells that block T-cell infiltration, but also provided a gene target pool for further exploration in improving lymphocyte infiltration and ICB efficacy in solid tumours.

FUNDING

A full list of funding bodies that contributed to this study can be found in the Acknowledgements section.

Dynamic changes of immunocyte subpopulations in thermogenic activation of adipose tissues

Objectives

The effects of cold exposure on whole-body metabolism in humans have gained increasing attention. Brown or beige adipose tissues are crucial in cold-induced thermogenesis to dissipate energy and thus have the potential to combat metabolic disorders. Despite the immune regulation of thermogenic adipose tissues, the overall changes in vital immune cells during distinct cold periods remain elusive. This study aimed to discuss the overall changes in immune cells under different cold exposure periods and to screen several potential immune cell subpopulations on thermogenic regulation.

Methods

Cibersort and mMCP-counter algorithms were employed to analyze immune infiltration in two (brown and beige) thermogenic adipose tissues under distinct cold periods. Changes in some crucial immune cell populations were validated by reanalyzing the single-cell sequencing dataset (GSE207706). Flow cytometry, immunofluorescence, and quantitative real-time PCR assays were performed to detect the proportion or expression changes in mouse immune cells of thermogenic adipose tissues under cold challenge.

Results

The proportion of monocytes, naïve, and memory T cells increased, while the proportion of NK cells decreased under cold exposure in brown adipose tissues.

Conclusion

Our study revealed dynamic changes in immune cell profiles in thermogenic adipose tissues and identified several novel immune cell subpopulations, which may contribute to thermogenic activation of adipose tissues under cold exposure.

Identification of potential immunotherapeutic targets and prognostic biomarkers in Graves' disease using weighted gene co-expression network analysis

Graves' disease (GD) is an autoimmune disorder characterized by hyperthyroidism resulting from autoantibody-induced stimulation of the thyroid gland. Despite recent advancements in understanding GD's pathogenesis, the molecular processes driving disease progression and treatment response remain poorly understood. In this study, we aimed to identify crucial immunogenic factors associated with GD prognosis and immunotherapeutic response. To achieve this, we implemented a comprehensive screening strategy that combined computational immunogenicity-potential scoring with multi-parametric cluster analysis to assess the immunomodulatory genes in GD-related subtypes involving stromal and immune cells. Utilizing weighted gene co-expression network analysis (WGCNA), we identified co-expressed gene modules linked to cellular senescence and immune infiltration in CD4+ and CD8+ GD samples. Additionally, gene set enrichment analysis enabled the identification of hallmark pathways distinguishing high- and low-immune subtypes. Our WGCNA analysis revealed 21 gene co-expression modules comprising 1,541 genes associated with immune infiltration components in various stages of GD, including T cells, M1 and M2 macrophages, NK cells, and Tregs. These genes primarily participated in T cell proliferation through purinergic signaling pathways, particularly neuroactive ligand-receptor interactions, and DNA binding transcription factor activity. Three genes, namely PRSS1, HCRTR1, and P2RY4, exhibited robustness in GD patients across multiple stages and were involved in immune cell infiltration during the late stage of GD (p < 0.05). Importantly, HCRTR1 and P2RY4 emerged as potential prognostic signatures for predicting overall survival in high-immunocore GD patients (p < 0.05). Overall, our study provides novel insights into the molecular mechanisms driving GD progression and highlights potential key immunogens for further investigation. These findings underscore the significance of immune infiltration-related cellular senescence in GD therapy and present promising targets for the development of new immunotherapeutic strategies.

An inducible Cd79b mutation confers ibrutinib sensitivity in mouse models of Myd88-driven diffuse large B-cell lymphoma

ABSTRACT

Diffuse large B-cell lymphoma (DLBCL) is the most common aggressive lymphoma and constitutes a highly heterogenous disease. Recent comprehensive genomic profiling revealed the identity of numerous molecularly defined DLBCL subtypes, including a cluster which is characterized by recurrent aberrations in MYD88, CD79B, and BCL2, as well as various lesions promoting a block in plasma cell differentiation, including PRDM1, TBL1XR1, and SPIB. Here, we generated a series of autochthonous mouse models to mimic this DLBCL cluster and specifically focused on the impact of Cd79b mutations in this setting. We show that canonical Cd79b immunoreceptor tyrosine-based activation motif (ITAM) mutations do not accelerate Myd88- and BCL2-driven lymphomagenesis. Cd79b-mutant murine DLBCL were enriched for IgM surface expression, reminiscent of their human counterparts. Moreover, Cd79b-mutant lymphomas displayed a robust formation of cytoplasmic signaling complexes involving MYD88, CD79B, MALT1, and BTK. These complexes were disrupted upon pharmacological BTK inhibition. The BTK inhibitor-mediated disruption of these signaling complexes translated into a selective ibrutinib sensitivity of lymphomas harboring combined Cd79b and Myd88 mutations. Altogether, this in-depth cross-species comparison provides a framework for the development of molecularly targeted therapeutic intervention strategies in DLBCL.

Making mouse transcriptomics deconvolution accessible with immunedeconv

Summary

Transcriptome deconvolution has emerged as a reliable technique to estimate cell-type abundances from bulk RNA sequencing data. Unlike their human equivalents, methods to quantify the cellular composition of complex tissues from murine transcriptomics are sparse and sometimes not easy to use. We extended the immunedeconv R package to facilitate the deconvolution of mouse transcriptomics, enabling the quantification of murine immune-cell types using 13 different methods. Through immunedeconv, we further offer the possibility of tweaking cell signatures used by deconvolution methods, providing custom annotations tailored for specific cell types and tissues. These developments strongly facilitate the study of the immune-cell composition of mouse models and further open new avenues in the investigation of the cellular composition of other tissues and organisms.

Availability and implementation

The R package and the documentation are available at https://github.com/omnideconv/immunedeconv.

Advancing the Boundaries of Immunotherapy in Lung Adenocarcinoma with Idiopathic Pulmonary Fibrosis by a Biomimetic Proteinoid Enabling Selective Endocytosis

The coexistence of lung adenocarcinoma (LUAD) with idiopathic pulmonary fibrosis (IPF), which has been extensively documented as a prominent risk factor for checkpoint inhibitor-related pneumonitis (CIP) in patients undergoing immunotherapy, has long been considered a restricted domain for the use of immune checkpoint inhibitors (ICIs). To overcome it, an approach was employed herein to specifically target PD-L1 within the cellular interior, surpassing the conventional focus solely on the cytomembrane, thereby facilitating the development of ICIs capable of distinguishing between LUAD cells and noncancerous cells based on their distinctive endocytic propensities. By exploiting the aurophilicity-driven self-assembly of a PD-L1 binding peptide (PDBP) and subsequently encapsulating it within erythrocyte membranes (EM), the resulting biomimetic ICIs protein EMS-PDBP exhibited extraordinary selectivity in internalizing LUAD cells, effectively targeting PD-L1 within cancer cells while hindering its membrane translocation. The EMS-PDBP treatment not only reactivated the antitumor immune response in the LUAD orthotopic allograft mouse model but also demonstrated a favorable safety profile by effectively eliminating any immune-related adverse events (irAEs). Most significantly, EMS-PDBP successfully and safely restored the antitumor immune response in a mouse model of LUAD with coexistent IPF, thus shattering the confines of ICIs immunotherapy. The reported EMS-PDBP collectively offers a potential strategy for immune reactivation to overcome the limitations of immunotherapy in LUAD coexisting with IPF.

SKP2 Knockout in Rb1/p53-Deficient Mouse Models of Osteosarcoma Induces Immune Infiltration and Drives a Transcriptional Program with a Favorable Prognosis

Osteosarcoma is an aggressive bone malignancy with a poor prognosis. One putative proto-oncogene in osteosarcoma is SKP2, encoding a substrate recognition factor of the SCF E3 ubiquitin ligase. We previously demonstrated that Skp2 knockout in murine osteosarcoma improved survival and delayed tumorigenesis. Here, we performed RNA sequencing (RNA-seq) on tumors from a transgenic osteosarcoma mouse model with conditional Trp53 and Rb1 knockouts in the osteoblast lineage ("DKO": Osx1-Cre;Rb1lox/lox;p53lox/lox) and a triple-knockout model with additional Skp2 germline knockout ("TKO": Osx1-Cre;Rb1lox/lox;p53lox/lox;Skp2-/-), followed by qPCR and immunohistochemistry validation. To investigate the clinical implications of our results, we analyzed a human osteosarcoma patient cohort ("NCI-TARGET OS") with RNA-seq and clinical data. We found large differences in gene expression after SKP2 knockout. Surprisingly, we observed increased expression of genes related to immune microenvironment infiltration in TKO tumors, especially the signature genes for macrophages and to a lesser extent, T cells, B cells, and vascular cells. We also uncovered a set of relevant transcription factors that may mediate these changes. In osteosarcoma patient cohorts, high expression of genes upregulated in TKO was correlated with favorable overall survival, which was largely explained by the macrophage gene signatures. This relationship was further supported by our finding that SKP2 expression was negatively correlated with macrophage infiltration in the NCI-TARGET osteosarcoma and the TCGA Sarcoma cohorts. Overall, our findings indicate that SKP2 may mediate immune exclusion from the osteosarcoma tumor microenvironment, suggesting that SKP2 modulation in osteosarcoma may induce antitumor immune activation.

Transcriptomic profiling of Polycystic Kidney Disease identifies paracrine factors in the early cyst microenvironment

Initial cysts that are formed upon Pkd1 loss in mice impose persistent stress on surrounding tissue and trigger a cystic snowball effect, in which local aberrant PKD-related signaling increases the likelihood of new cyst formation, ultimately leading to accelerated disease progression. Although many pathways have been associated with PKD progression, the knowledge of early changes near initial cysts is limited. To perform an unbiased analysis of transcriptomic alterations in the cyst microenvironment, microdomains were collected from kidney sections of iKsp-Pkd1del mice with scattered Pkd1-deletion using Laser Capture Microdissection. These microdomains were defined as F4/80-low cystic, representing early alterations in the cyst microenvironment, F4/80-high cystic, with more advanced alterations, or non-cystic. RNA sequencing and differential gene expression analysis revealed 953 and 8088 dysregulated genes in the F4/80-low and F4/80-high cyst microenvironment, respectively, when compared to non-cystic microdomains. In the early cyst microenvironment, several injury-repair, growth, and tissue remodeling-related pathways were activated, accompanied by mild metabolic changes. In the more advanced F4/80-high microdomains, these pathways were potentiated and the metabolism was highly dysregulated. Upstream regulator analysis revealed a series of paracrine factors with increased activity in the early cyst microenvironment, including TNFSF12 and OSM. In line with the upstream regulator analysis, TWEAK and Oncostatin-M promoted cell proliferation and inflammatory gene expression in renal epithelial cells and fibroblasts in vitro. Collectively, our data provide an overview of molecular alterations that specifically occur in the cyst microenvironment and identify paracrine factors that may mediate early and advanced alterations in the cyst microenvironment.

Analysis of multiple programmed cell death-related prognostic genes and functional validations of necroptosis-associated genes in oesophageal squamous cell carcinoma

BACKGROUND

Oesophageal squamous cell carcinoma (ESCC) is a lethal malignancy. Immune checkpoint inhibitors (ICIs) showed great clinical benefits for patients with ESCC. We aimed to construct a model predicting prognosis and response to ICIs by integrating diverse programmed cell death (PCD) forms.

METHODS

Genes related to 14 PCDs were collected to generate multi-gene signatures, including apoptosis, necroptosis, pyroptosis, ferroptosis, and cuproptosis. Bulk and single-cell RNA transcriptome datasets were used to develop and validate the model. We assessed the functions of two necroptosis-related genes in ESCC cells by Western blot, co-immunoprecipitation (Co-IP), LDH release assay, CCK-8, and migration assay, followed by immunohistochemistry (IHC) staining on samples of patients with ESCC (n = 67).

FINDINGS

We built and validated a 16-gene prognostic combined cell death index (CCDI) by combining immunogenic cell death (ICD) and necroptosis signatures. The CCDI could also predict response to ICIs in cancer, as shown by Tumour Immune Dysfunction and Exclusion (TIDE) analysis, confirmed in four independent ICI clinical trials. Trajectory analysis revealed that HOOK1 and CUL4A might affect ESCC cell fate. We found that HOOK1 induced necroptosis and inhibited the proliferation and migration of ESCC cells, while CUL4A exhibited the opposite effects. Co-IP assay confirmed that HOOK1 and CUL4A promoted and reduced necrosome formation in ESCC cells. Data from patients with ESCC further supported that HOOK1 and CUL4A might be a tumour suppressor and oncogene, respectively.

INTERPRETATION

We constructed a CCDI model with potential in predicting prognosis and response to ICIs in cancer. HOOK1 and CUL4A in the CCDI model are crucial prognostic biomarkers in ESCC.

FUNDING

The Natural Science Foundation of China [82172786], The National Cancer Center Climbing Fund of China [NCC201908B06], The Natural Science Foundation of Heilongjiang Province [LH2021H077].

Antigen presentation deficiency, mesenchymal differentiation, and resistance to immunotherapy in the murine syngeneic CT2A tumor model

Background

The GL261 and CT2A syngeneic tumor lines are frequently used as immunocompetent orthotopic mouse models of human glioblastoma (huGBM) but demonstrate distinct differences in their responses to immunotherapy.

Methods

To decipher the cell-intrinsic mechanisms that drive immunotherapy resistance in CT2A-luc and to define the aspects of human cancer biology that these lines can best model, we systematically compared their characteristics using whole exome and transcriptome sequencing, and protein analysis through immunohistochemistry, Western blot, flow cytometry, immunopeptidomics, and phosphopeptidomics.

Results

The transcriptional profiles of GL261-luc2 and CT2A-luc tumors resembled those of some huGBMs, despite neither line sharing the essential genetic or histologic features of huGBM. Both models exhibited striking hypermutation, with clonal hotspot mutations in RAS genes (Kras p.G12C in GL261-luc2 and Nras p.Q61L in CT2A-luc). CT2A-luc distinctly displayed mesenchymal differentiation, upregulated angiogenesis, and multiple defects in antigen presentation machinery (e.g. Tap1 p.Y488C and Psmb8 p.A275P mutations) and interferon response pathways (e.g. copy number losses of loci including IFN genes and reduced phosphorylation of JAK/STAT pathway members). The defect in MHC class I expression could be overcome in CT2A-luc by interferon-γ treatment, which may underlie the modest efficacy of some immunotherapy combinations. Additionally, CT2A-luc demonstrated substantial baseline secretion of the CCL-2, CCL-5, and CCL-22 chemokines, which play important roles as myeloid chemoattractants.

Conclusion

Although the clinical contexts that can be modeled by GL261 and CT2A for huGBM are limited, CT2A may be an informative model of immunotherapy resistance due to its deficits in antigen presentation machinery and interferon response pathways.

Molecular and clinical characterization of ICOS expression in breast cancer through large-scale transcriptome data

ICOS (Inducible T Cell Costimulator), one of the co-stimulatory B7 superfamily members, was characterized as a co-stimulatory receptor for T-cell enhancement. However, the role of ICOS in breast cancer remains largely unknown. The present study systematically investigated the expression pattern and its relation to clinical characteristics and immunotherapy by integrating multiple clinical cohorts and large-scale gene expression data. This study included 2994 breast tumor samples with transcriptome data and matched clinical data. To make our findings more reliable, we set the TCGA cohort as the discovery set and the METABRIC cohort as the validation set. The expression of ICOS in breast cancer is strongly associated with major clinical and molecular characteristics. There is an association between higher ICOS expression and malignant subtypes and grades of tumors. In addition, gene ontology analysis based on genes significantly correlated with ICOS expression indicated that the expression of ICOS is mainly associated with immune responses and inflammation. We also observed strong correlations between ICOS and other promising immune-checkpoint molecules, including PD1, PDL1, CTLA4, and IDO1. Furthermore, we found that ICOS expression is associated with the response to anti-PDL1 immunotherapy and may serve as a biomarker for immunotherapy prediction. Our results indicated higher ICOS expression is significantly associated with favorable survival in triple-negative breast cancer (TNBC) patients, but not for all subtypes of breast cancer patients. In summary, ICOS correlates with higher malignant breast cancers, and it contributes to the regulation of the immune microenvironment of breast tumors, making it a potential biomarker and immunotherapy target.

Targeting the lipid kinase PIKfyve upregulates surface expression of MHC class I to augment cancer immunotherapy

Despite the remarkable clinical success of immunotherapies in a subset of cancer patients, many fail to respond to treatment and exhibit resistance. Here, we found that genetic or pharmacologic inhibition of the lipid kinase PIKfyve, a regulator of autophagic flux and lysosomal biogenesis, upregulated surface expression of major histocompatibility complex class I (MHC-I) in cancer cells via impairing autophagic flux, resulting in enhanced cancer cell killing mediated by CD8+ T cells. Genetic depletion or pharmacologic inhibition of PIKfyve elevated tumor-specific MHC-I surface expression, increased intratumoral functional CD8+ T cells, and slowed tumor progression in multiple syngeneic mouse models. Importantly, enhanced antitumor responses by Pikfyve-depletion were CD8+ T cell- and MHC-I-dependent, as CD8+ T cell depletion or B2m knockout rescued tumor growth. Furthermore, PIKfyve inhibition improved response to immune checkpoint blockade (ICB), adoptive cell therapy, and a therapeutic vaccine. High expression of PIKFYVE was also predictive of poor response to ICB and prognostic of poor survival in ICB-treated cohorts. Collectively, our findings show that targeting PIKfyve enhances immunotherapies by elevating surface expression of MHC-I in cancer cells, and PIKfyve inhibitors have potential as agents to increase immunotherapy response in cancer patients.

ICBcomb: a comprehensive expression database for immune checkpoint blockade combination therapy

The success of immune checkpoint blockade (ICB) promotes the immunotherapy to be a new pillar in cancer treatment. However, the low response rate of the ICB therapy limits its application. To increase the response rate and enhance efficacy, the ICB combination therapy has emerged and its clinical trials are increasing. Nevertheless, the gene expression profile and its pattern of ICB combination were not comprehensively studied, which limits the understanding of the ICB combination therapy and the identification of new drugs. Here, we constructed ICBcomb (http://bioinfo.life.hust.edu.cn/ICBcomb/), a comprehensive database, by analyzing the human and mouse expression data of the ICB combination therapy and comparing them between groups treated with ICB, other drugs or their combinations. ICBcomb contains 1399 samples across 29 cancer types involving 52 drugs. It provides a user-friendly web interface for demonstrating the results of the available comparisons in the ICB combination therapy datasets with five functional modules: [1, 2] the 'Dataset/Disease' modules for browsing the expression, enrichment and comparison results in each dataset or disease; [3] the 'Gene' module for inputting a gene symbol and displaying its expression and comparison results across datasets/diseases; [4] the 'Gene Set' module for GSVA/GSEA enrichment analysis on the built-in gene sets and the user-input gene sets in different comparisons; [5] the 'Immune Cell' module for immune cell infiltration comparison between different groups by immune cell abundance analysis. The ICBcomb database provides the first resource for gene expression profile and comparison in ICB combination therapy, which may provide clues for discovering the mechanism of effective combination strategies and new combinatory drugs.

Ly6C-high monocytes alleviate brain injury in experimental subarachnoid hemorrhage in mice

BACKGROUND

Subarachnoid hemorrhage (SAH) is an uncommon type of potentially fatal stroke. The pathophysiological mechanisms of brain injury remain unclear, which hinders the development of drugs for SAH. We aimed to investigate the pathophysiological mechanisms of SAH and to elucidate the cellular and molecular biological response to SAH-induced injury.

METHODS

A cross-species (human and mouse) multiomics approach combining high-throughput data and bioinformatic analysis was used to explore the key pathophysiological processes and cells involved in SAH-induced brain injury. Patient data were collected from the hospital (n = 712). SAH was established in adult male mice via endovascular perforation, and flow cytometry, a bone marrow chimera model, qPCR, and microglial depletion experiments were conducted to explore the origin and chemotaxis mechanism of the immune cells. To investigate cell effects on SAH prognosis, murine neurological function was evaluated based on a modified Garcia score, pole test, and rotarod test.

RESULTS

The bioinformatics analysis confirmed that inflammatory and immune responses were the key pathophysiological processes after SAH. Significant increases in the monocyte levels were observed in both the mouse brains and the peripheral blood of patients after SAH. Ly6C-high monocytes originated in the bone marrow, and the skull bone marrow contribute a higher proportion of these monocytes than neutrophils. The mRNA level of Ccl2 was significantly upregulated after SAH and was greater in CD11b-positive than CD11b-negative cells. Microglial depletion, microglial inhibition, and CCL2 blockade reduced the numbers of Ly6C-high monocytes after SAH. With CCR2 antagonization, the neurological function of the mice exhibited a slow recovery. Three days post-SAH, the monocyte-derived dendritic cell (moDC) population had a higher proportion of TNF-α-positive cells and a lower proportion of IL-10-positive cells than the macrophage population. The ratio of moDCs to macrophages was higher on day 3 than on day 5 post-SAH.

CONCLUSIONS

Inflammatory and immune responses are significantly involved in SAH-induced brain injury. Ly6C-high monocytes derived from the bone marrow, including the skull bone marrow, infiltrated into mouse brains via CCL2 secreted from microglia. Moreover, Ly6C-high monocytes alleviated neurological dysfunction after SAH.

Rat Deconvolution as Knowledge Miner for Immune Cell Trafficking from Toxicogenomics Databases

Toxicogenomics databases are useful for understanding biological responses in individuals because they include a diverse spectrum of biological responses. Although these databases contain no information regarding immune cells in the liver, which are important in the progression of liver injury, deconvolution that estimates cell-type proportions from bulk transcriptome could extend immune information. However, deconvolution has been mainly applied to humans and mice and less often to rats, which are the main target of toxicogenomics databases. Here, we developed a deconvolution method for rats to retrieve information regarding immune cells from toxicogenomics databases. The rat-specific deconvolution showed high correlations for several types of immune cells between spleen and blood, and between liver treated with toxicants compared with those based on human and mouse data. Additionally, we found 4 clusters of compounds in Open TG-GATEs database based on estimated immune cell trafficking, which are different from those based on transcriptome data itself. The contributions of this work are three-fold. First, we obtained the gene expression profiles of 6 rat immune cells necessary for deconvolution. Second, we clarified the importance of species differences on deconvolution. Third, we retrieved immune cell trafficking from toxicogenomics databases. Accumulated and comparable immune cell profiles of massive data of immune cell trafficking in rats could deepen our understanding of enable us to clarify the relationship between the order and the contribution rate of immune cells, chemokines and cytokines, and pathologies. Ultimately, these findings will lead to the evaluation of organ responses in Adverse Outcome Pathway.

SETD2 Loss and ATR Inhibition Synergize to Promote cGAS Signaling and Immunotherapy Response in Renal Cell Carcinoma

PURPOSE

Immune checkpoint blockade (ICB) demonstrates durable clinical benefits in a minority of patients with renal cell carcinoma (RCC). We aimed to identify the molecular features that determine the response and develop approaches to enhance it.

EXPERIMENTAL DESIGN

We investigated the effects of SET domain-containing protein 2 (SETD2) loss on the DNA damage response pathway, the cytosolic DNA-sensing pathway, the tumor immune microenvironment, and the response to ataxia telangiectasia and rad3-related (ATR) and checkpoint inhibition in RCC.

RESULTS

ATR inhibition activated the cyclic GMP-AMP synthase (cGAS)-interferon regulatory factor 3 (IRF3)-dependent cytosolic DNA-sensing pathway, resulting in the concurrent expression of inflammatory cytokines and immune checkpoints. Among the common RCC genotypes, SETD2 loss is associated with preferential ATR activation and sensitizes cells to ATR inhibition. SETD2 knockdown promoted the cytosolic DNA-sensing pathway in response to ATR inhibition. Treatment with the ATR inhibitor VE822 concurrently upregulated immune cell infiltration and immune checkpoint expression in Setd2 knockdown Renca tumors, providing a rationale for ATR inhibition plus ICB combination therapy. Setd2-deficient Renca tumors demonstrated greater vulnerability to ICB monotherapy or combination therapy with VE822 than Setd2-proficient tumors. Moreover, SETD2 mutations were associated with a higher response rate and prolonged overall survival in patients with ICB-treated RCC but not in patients with non-ICB-treated RCC.

CONCLUSIONS

SETD2 loss and ATR inhibition synergize to promote cGAS signaling and enhance immune cell infiltration, providing a mechanistic rationale for the combination of ATR and checkpoint inhibition in patients with RCC with SETD2 mutations.

Photodynamic augmentation of oncolytic virus therapy for central nervous system malignancies

Oncolytic viruses (OVs) have emerged as a clinical therapeutic modality potentially effective for cancers that evade conventional therapies, including central nervous system malignancies. Rationally designed combinatorial strategies can augment the efficacy of OVs by boosting tumor-selective cytotoxicity and modulating the tumor microenvironment (TME). Photodynamic therapy (PDT) of cancer not only mediates direct neoplastic cell death but also primes the TME to sensitize the tumor to secondary therapies, allowing for the combination of two potentially synergistic therapies with broader targets. Here, we created G47Δ-KR, clinical oncolytic herpes simplex virus G47Δ that expresses photosensitizer protein KillerRed (KR). Optical properties and cytotoxic effects of G47Δ-KR infection followed by amber LED illumination (peak wavelength: 585-595 nm) were examined in human glioblastoma (GBM) and malignant meningioma (MM) models in vitro. G47Δ-KR infection of tumor cells mediated KR expression that was activated by LED and produced reactive oxygen species, leading to cell death that was more robust than G47Δ-KR without light. In vivo, we tested photodynamic-oncolytic virus (PD-OV) therapy employing intratumoral injection of G47Δ-KR followed by laser light tumor irradiation (wavelength: 585 nm) in GBM and MM xenografts. PD-OV therapy was feasible in these models and resulted in potent anti-tumor effects that were superior to G47Δ-KR alone (without laser light) or laser light alone. RNA sequencing analysis of post-treatment tumor samples revealed PD-OV therapy-induced increases in TME infiltration of variable immune cell types. This study thus demonstrated the proof-of-concept that G47Δ-KR enables PD-OV therapy for neuro-oncological malignancies and warrants further research to advance potential clinical translation.

Exploring the role of Hmox1 in ferroptosis and immune infiltration during renal ischemia-reperfusion injury

Ischemia-reperfusion injury (IRI) is inevitable in kidney transplantations and, as a complex pathophysiological process, it can be greatly impacted by ferroptosis and immune inflammation. Our study aimed to identify the biomarkers of renal IRI (RIRI) and elucidate their relationship with immune infiltration. In this study, the GSE148420 database was used as a training set to analyze differential genes and overlap them with ferroptosis-related genes to identify hub genes using a protein-protein interaction (PPI) network, the least absolute shrinkage and selection operator (LASSO), and random forest algorithm (RFA). We verified the hub gene and ferroptosis-related phenotypes in a verification set and animal experiments involving unilateral IRI with contralateral nephrectomy in rats. Gene set enrichment analysis (GSEA) of single genes was conducted according to the hub gene to predict related endogenous RNAs (ceRNAs) and drugs to establish a network. Finally, we used the Cibersort to analyze immunological infiltration and conducted Spearman's correlation analysis. We identified 5456 differential genes and obtained 26 ferroptosis-related differentially expressed genes. Through PPI, LASSO, and RFA, Hmox1 was identified as the only hub gene and its expression levels were verified using verification sets. In animal experiments, Hmox1 was verified as a key biomarker. GSEA of single genes revealed the seven most related pathways, and the ceRNAs network included 138 mRNAs and miRNAs. We predicted 11 related drugs and their three-dimensional structural maps. Thus, Hmox1 was identified as a key biomarker and regulator of ferroptosis in RIRI and its regulation of ferroptosis was closely related to immune infiltration.

Sexually dimorphic effect of H-ferritin genetic manipulation on survival and tumor microenvironment in a mouse model of glioblastoma

PURPOSE

Iron plays a crucial role in various biological mechanisms and has been found to promote tumor growth. Recent research has shown that the H-ferritin (FTH1) protein, traditionally recognized as an essential iron storage protein, can transport iron to GBM cancer stem cells, reducing their invasion activity. Moreover, the binding of extracellular FTH1 to human GBM tissues, and brain iron delivery in general, has been found to have a sex bias. These observations raise questions, addressed in this study, about whether H-ferritin levels extrinsic to the tumor can affect tumor cell pathways and if this impact is sex-specific.

METHODS

To interrogate the role of systemic H-ferritin in GBM we introduce a mouse model in which H-ferritin levels are genetically manipulated. Mice that were genetically manipulated to be heterozygous for H-ferritin (Fth1+/-) gene expression were orthotopically implanted with a mouse GBM cell line (GL261). Littermate Fth1 +/+ mice were used as controls. The animals were evaluated for survival and the tumors were subjected to RNA sequencing protocols. We analyzed the resulting data utilizing the murine Microenvironment Cell Population (mMCP) method for in silico immune deconvolution. mMCP analysis estimates the abundance of tissue infiltrating immune and stromal populations based on cell-specific gene expression signatures.

RESULTS

There was a clear sex bias in survival. Female Fth1+/- mice had significantly poorer survival than control females (Fth1+/+). The Fth1 genetic status did not affect survival in males. The mMCP analysis revealed a significant reduction in T cells and CD8 + T cell infiltration in the tumors of females with Fth1+/- background as compared to the Fth1+/+. Mast and fibroblast cell infiltration was increased in females and males with Fth1+/- background, respectively, compared to Fth1+/+ mice.

CONCLUSION

Genetic manipulation of Fth1 which leads to reduced systemic levels of FTH1 protein had a sexually dimorphic impact on survival. Fth1 heterozygosity significantly worsened survival in females but did not affect survival in male GBMs. Furthermore, the genetic manipulation of Fth1 significantly affected tumor infiltration of T-cells, CD8 + T cells, fibroblasts, and mast cells in a sexually dimorphic manner. These results demonstrate a role for FTH1 and presumably iron status in establishing the tumor cellular landscape that ultimately impacts survival and further reveals a sex bias that may inform the population studies showing a sex effect on the prevalence of brain tumors.

Effectively utilizing publicly available databases for cancer target evaluation

The majority of compounds designed against cancer drug targets do not progress to become approved drugs, mainly due to lack of efficacy and/or unmanageable toxicity. Robust target evaluation is therefore required before progressing through the drug discovery process to reduce the high attrition rate. There are a wealth of publicly available databases that can be mined to generate data as part of a target evaluation. It can, however, be challenging to learn what databases are available, how and when they should be used, and to understand the associated limitations. Here, we have compiled and present key, freely accessible and easy-to-use databases that house informative datasets from in vitro, in vivo and clinical studies. We also highlight comprehensive target review databases that aim to bring together information from multiple sources into one-stop portals. In the post-genomics era, a key objective is to exploit the extensive cell, animal and patient characterization datasets in order to deliver precision medicine on a patient-specific basis. Effective utilization of the highlighted databases will go some way towards supporting the cancer research community achieve these aims.

The histone chaperone function of Daxx is dispensable for embryonic development

Daxx functions as a histone chaperone for the histone H3 variant, H3.3, and is essential for embryonic development. Daxx interacts with Atrx to form a protein complex that deposits H3.3 into heterochromatic regions of the genome, including centromeres, telomeres, and repeat loci. To advance our understanding of histone chaperone activity in vivo, we developed two Daxx mutant alleles in the mouse germline, which abolish the interactions between Daxx and Atrx (DaxxY130A), and Daxx and H3.3 (DaxxS226A). We found that the interaction between Daxx and Atrx is dispensable for viability; mice are born at the expected Mendelian ratio and are fertile. The loss of Daxx-Atrx interaction, however, does cause dysregulated expression of endogenous retroviruses. In contrast, the interaction between Daxx and H3.3, while not required for embryonic development, is essential for postnatal viability. Transcriptome analysis of embryonic tissues demonstrates that this interaction is important for silencing endogenous retroviruses and for maintaining proper immune cell composition. Overall, these results clearly demonstrate that Daxx has both Atrx-dependent and independent functions in vivo, advancing our understanding of this epigenetic regulatory complex.

HPVTIMER: A shiny web application for tumor immune estimation in human papillomavirus-associated cancers

The tumor immune microenvironment (TIME) is closely associated with tumor formation, particularly linked to the human papillomavirus (HPV), and regulates tumor initiation, proliferation, infiltration, and metastasis. With the rise of immunotherapy, an increasing amount of sample data used for TIME exploration is available in databases. However, no currently available web tool enables a comprehensive exploration of the TIME of HPV-associated cancers by leveraging these data. We have developed a web tool called HPV-associated Tumor Immune MicroEnvironment ExploreR (HPVTIMER), which provides a comprehensive analysis platform that integrates over 10,000 genes and 2290 tumor samples from 65 transcriptome data sets across 8 cancer types sourced from the Gene Expression Omnibus (GEO) database. The tool features four built-in analysis modules, namely, the differential expression analysis module, correlation analysis module, immune infiltration analysis module, and pathway analysis module. These modules enable users to perform systematic and vertical analyses. We used several analytical modules in HPVTIMER to briefly explore the role of CDKN2A in head and neck squamous cell carcinomas. We expect that HPVTIMER will help users explore the immune microenvironment of HPV-associated cancers and uncover potential immune regulatory mechanisms and immunotherapeutic targets. HPVTIMER is available at http://www.hpvtimer.com/.

CXCR2 expression during melanoma tumorigenesis controls transcriptional programs that facilitate tumor growth

BACKGROUND

Though the CXCR2 chemokine receptor is known to play a key role in cancer growth and response to therapy, a direct link between expression of CXCR2 in tumor progenitor cells during induction of tumorigenesis has not been established.

METHODS

To characterize the role of CXCR2 during melanoma tumorigenesis, we generated tamoxifen-inducible tyrosinase-promoter driven BrafV600E/Pten-/-/Cxcr2-/- and NRasQ61R/INK4a-/-/Cxcr2-/- melanoma models. In addition, the effects of a CXCR1/CXCR2 antagonist, SX-682, on melanoma tumorigenesis were evaluated in BrafV600E/Pten-/- and NRasQ61R/INK4a-/- mice and in melanoma cell lines. Potential mechanisms by which Cxcr2 affects melanoma tumorigenesis in these murine models were explored using RNAseq, mMCP-counter, ChIPseq, and qRT-PCR; flow cytometry, and reverse phosphoprotein analysis (RPPA).

RESULTS

Genetic loss of Cxcr2 or pharmacological inhibition of CXCR1/CXCR2 during melanoma tumor induction resulted in key changes in gene expression that reduced tumor incidence/growth and increased anti-tumor immunity. Interestingly, after Cxcr2 ablation, Tfcp2l1, a key tumor suppressive transcription factor, was the only gene significantly induced with a log2 fold-change greater than 2 in these three different melanoma models.

CONCLUSIONS

Here, we provide novel mechanistic insight revealing how loss of Cxcr2 expression/activity in melanoma tumor progenitor cells results in reduced tumor burden and creation of an anti-tumor immune microenvironment. This mechanism entails an increase in expression of the tumor suppressive transcription factor, Tfcp2l1, along with alteration in the expression of genes involved in growth regulation, tumor suppression, stemness, differentiation, and immune modulation. These gene expression changes are coincident with reduction in the activation of key growth regulatory pathways, including AKT and mTOR.

SKP2 knockout in Rb1/p53 deficient mouse models of osteosarcoma induces immune infiltration and drives a transcriptional program with a favorable prognosis

Purpose

Osteosarcoma (OS) is an aggressive bone malignancy with a poor prognosis. One putative proto-oncogene in OS is SKP2, encoding a substrate recognition factor of the SCF E3 ubiquitin ligase. We previously demonstrated that SKP2 knockout in murine OS improved survival and delayed tumorigenesis. Here we aim to define the SKP2 drives transcriptional program and its clinical implication in OS.

Experimental Design

We performed RNA-sequencing (RNA-seq) on tumors from a transgenic OS mouse model with conditional Trp53 and Rb1 knockouts in the osteoblast lineage ("DKO": Osx1-Cre;Rb1lox/lox;p53lox/lox) and a triple-knockout model with additional Skp2 germline knockout ("TKO": Osx1-Cre;Rb1lox/lox;p53lox/lox;SKP2-/-). We validated our RNA-seq findings using qPCR and immunohistochemistry. To investigate the clinical implications of our results, we analyzed a human OS patient cohort ("NCI-TARGET OS") with RNA-seq and clinical data.

Results

We found large differences in gene expression after SKP2 knockout. Strikingly, we observed increased expression of genes related to immune microenvironment infiltration in TKO tumors. We observed significant increases in signature genes for macrophages and to a lesser extent, T cells, B cells and vascular cells. We also uncovered a set of relevant transcription factors that may mediate the changes. In OS patient cohorts, high expression of genes upregulated in TKO was correlated with favorable overall survival, which was largely explained by the macrophage gene signatures. This relationship was further supported by our finding that SKP2 expression was negatively correlated with macrophage infiltration in the NCI-TARGET OS and the TCGA Sarcoma cohort.

Conclusion

Our findings indicate that SKP2 may mediate immune exclusion from the OS tumor microenvironment, suggesting that SKP2 modulation in OS may induce anti-tumor immune activation.

Potential clinical use of azacitidine and MEK inhibitor combination therapy in PTPN11-mutated juvenile myelomonocytic leukemia

Juvenile myelomonocytic leukemia (JMML) is a rare myeloproliferative neoplasm of childhood. The molecular hallmark of JMML is hyperactivation of the Ras/MAPK pathway with the most common cause being mutations in the gene PTPN11, encoding the protein tyrosine phosphatase SHP2. Current strategies for treating JMML include using the hypomethylating agent, 5-azacitidine (5-Aza) or MEK inhibitors trametinib and PD0325901 (PD-901), but none of these are curative as monotherapy. Utilizing an Shp2E76K/+ murine model of JMML, we show that the combination of 5-Aza and PD-901 modulates several hematologic abnormalities often seen in JMML patients, in part by reducing the burden of leukemic hematopoietic stem and progenitor cells (HSC/Ps). The reduced JMML features in drug-treated mice were associated with a decrease in p-MEK and p-ERK levels in Shp2E76K/+ mice treated with the combination of 5-Aza and PD-901. RNA-sequencing analysis revealed a reduction in several RAS and MAPK signaling-related genes. Additionally, a decrease in the expression of genes associated with inflammation and myeloid leukemia was also observed in Shp2E76K/+ mice treated with the combination of the two drugs. Finally, we report two patients with JMML and PTPN11 mutations treated with 5-Aza, trametinib, and chemotherapy who experienced a clinical response because of the combination treatment.

CXCR2 expression during melanoma tumorigenesis controls transcriptional programs that facilitate tumor growth

Background

Though the CXCR2 chemokine receptor is known to play a key role in cancer growth and response to therapy, a direct link between expression of CXCR2 in tumor progenitor cells during induction of tumorigenesis has not been established.

Methods

To characterize the role of CXCR2 during melanoma tumorigenesis, we generated tamoxifen-inducible tyrosinase-promoter driven Braf ^V600E /Pten ^-/- /Cxcr2 ^-/- and NRas ^Q61R /INK4a ^-/- /Cxcr2 ^-/- melanoma models. In addition, the effects of a CXCR1/CXCR2 antagonist, SX-682, on melanoma tumorigenesis were evaluated in Braf ^V600E /Pten ^-/- and NRas ^Q61R /INK4a ^-/- mice and in melanoma cell lines. Potential mechanisms by which Cxcr2 affects melanoma tumorigenesis in these murine models were explored using RNAseq, mMCP-counter, ChIPseq, and qRT-PCR; flow cytometry, and reverse phosphoprotein analysis (RPPA).

Results

Genetic loss of Cxcr2 or pharmacological inhibition of CXCR1/CXCR2 during melanoma tumor induction resulted in key changes in gene expression that reduced tumor incidence/growth and increased anti-tumor immunity. Interestingly, after Cxcr2 ablation, Tfcp2l1 , a key tumor suppressive transcription factor, was the only gene significantly induced with a log ₂ fold-change greater than 2 in these three different melanoma models.

Conclusions

Here, we provide novel mechanistic insight revealing how loss of Cxcr2 expression/activity in melanoma tumor progenitor cells results in reduced tumor burden and creation of an anti-tumor immune microenvironment. This mechanism entails an increase in expression of the tumor suppressive transcription factor, Tfcp2l1, along with alteration in the expression of genes involved in growth regulation, tumor suppression, stemness, differentiation, and immune modulation. These gene expression changes are coincident with reduction in the activation of key growth regulatory pathways, including AKT and mTOR.

COL11A1-driven positive feedback loop modulates fibroblast transformation and activates pancreatic cancer progression

Pancreatic ductal adenocarcinoma (PDAC) is one of the most common leading causes of cancer death. The cancer-associated fibroblasts (CAFs) in the tumor microenvironment (TME) aggravate the malignant behavior of PDAC. However, it is still unknown how PDAC induces normal fibroblasts (NFs) to CAFs. In present research, we found that PDAC-derived collagen type XI alpha 1 (COL11A1) promoted the conversion of NFs to CAF-like cells. It included morphological and corresponding molecular marker changes. Activation of the nuclear factor-κB (NF-κB) pathway was involved in this process. Corresponding, CAFs cells could secrete interleukin 6 (IL-6), which promoted the invasion and the epithelial-mesenchymal transition of PDAC cells. Furthermore, IL-6 promoted the expression of transcription factor Activating Transcription Factor 4 by activating the Mitogen-Activated Protein Kinase/extracellular-signal-regulated kinase pathway. The latter directly promotes the expression of COL11A1. This way, a feedback loop of mutual influence was constructed between PDAC and CAFs. Our research proposed a novel concept for PDAC-educated NFs. PDAC-COL11A1-fibroblast-IL-6-PDAC axis might contribute to the cascade between PDAC and TME.

MmCMS: mouse models' consensus molecular subtypes of colorectal cancer

BACKGROUND

Colorectal cancer (CRC) primary tumours are molecularly classified into four consensus molecular subtypes (CMS1-4). Genetically engineered mouse models aim to faithfully mimic the complexity of human cancers and, when appropriately aligned, represent ideal pre-clinical systems to test new drug treatments. Despite its importance, dual-species classification has been limited by the lack of a reliable approach. Here we utilise, develop and test a set of options for human-to-mouse CMS classifications of CRC tissue.

METHODS

Using transcriptional data from established collections of CRC tumours, including human (TCGA cohort; n = 577) and mouse (n = 57 across n = 8 genotypes) tumours with combinations of random forest and nearest template prediction algorithms, alongside gene ontology collections, we comprehensively assess the performance of a suite of new dual-species classifiers.

RESULTS

We developed three approaches: MmCMS-A; a gene-level classifier, MmCMS-B; an ontology-level approach and MmCMS-C; a combined pathway system encompassing multiple biological and histological signalling cascades. Although all options could identify tumours associated with stromal-rich CMS4-like biology, MmCMS-A was unable to accurately classify the biology underpinning epithelial-like subtypes (CMS2/3) in mouse tumours.

CONCLUSIONS

When applying human-based transcriptional classifiers to mouse tumour data, a pathway-level classifier, rather than an individual gene-level system, is optimal. Our R package enables researchers to select suitable mouse models of human CRC subtype for their experimental testing.

General transcription factor TAF4 antagonizes epigenetic silencing by Polycomb to maintain intestine stem cell functions

Taf4 (TATA-box binding protein-associated factor 4) is a subunit of the general transcription factor TFIID, a component of the RNA polymerase II pre-initiation complex that interacts with tissue-specific transcription factors to regulate gene expression. Properly regulated gene expression is particularly important in the intestinal epithelium that is constantly renewed from stem cells. Tissue-specific inactivation of Taf4 in murine intestinal epithelium during embryogenesis compromised gut morphogenesis and the emergence of adult-type stem cells. In adults, Taf4 loss impacted the stem cell compartment and associated Paneth cells in the stem cell niche, epithelial turnover and differentiation of mature cells, thus exacerbating the response to inflammatory challenge. Taf4 inactivation ex vivo in enteroids prevented budding formation and maintenance and caused broad chromatin remodeling and a strong reduction in the numbers of stem and progenitor cells with a concomitant increase in an undifferentiated cell population that displayed high activity of the Ezh2 and Suz12 components of Polycomb Repressive Complex 2 (PRC2). Treatment of Taf4-mutant enteroids with a specific Ezh2 inhibitor restored buddings, cell proliferation and the stem/progenitor compartment. Taf4 loss also led to increased PRC2 activity in cells of adult crypts associated with modification of the immune/inflammatory microenvironment that potentiated Apc-driven tumorigenesis. Our results reveal a novel function of Taf4 in antagonizing PRC2-mediated repression of the stem cell gene expression program to assure normal development, homeostasis, and immune-microenvironment of the intestinal epithelium.

Deep dive into the immune response against murine mesothelioma permits design of novel anti-mesothelioma therapeutics

Given the need to improve the efficacy of standard-of-care immunotherapy (anti-CTLA-4 + anti-PD-1) in human malignant pleural mesothelioma (hMPM), we thoroughly characterized the immunobiology of the AB12 murine mesothelioma (MM) model, aiming to increase its accuracy in predicting the response of hMPM to immunotherapy and in designing novel anti-hMPM treatments. Specifically, we used immunologic, transcriptomic and survival analyses, to synchronize the MM tumor growth phases and immune evolution with the histo-molecular and immunological characteristics of hMPM while also determining the anti-MM efficacy of standard-of-care anti-hMPM immunotherapy as a benchmark that novel therapeutics should meet. We report that early-, intermediate- and advanced- AB12 tumors are characterized by a bell-shaped anti-tumor response that peaks in intermediate tumors and decays in advanced tumors. We further show that intermediate- and advanced- tumors match with immune active ("hot") and immune inactive ("cold") hMPM respectively, and that they respond to immunotherapy in a manner that corresponds well with its performance in real-life settings. Finally, we show that in advanced tumors, addition of cisplatin to anti CTLA-4 + anti PD-1 can extend mice survival and invigorate the decaying anti-tumor response. Therefore, we highlight this triple combination as a worthy candidate to improve clinical outcomes in hMPM.

PCSK9 facilitates melanoma pathogenesis via a network regulating tumor immunity

BACKGROUND

PCSK9 regulates cholesterol homeostasis and promotes tumorigenesis. However, the relevance of these two actions and the mechanisms underlying PCSK9's oncogenic roles in melanoma and other cancers remain unclear.

METHODS

PCSK9's association with melanoma was analysed using the TCGA dataset. Empty vector (EV), PCSK9, gain-of-function (D374Y), and loss-of-function (Q152H) PCSK9 mutant were stably-expressed in murine melanoma B16 cells and studied for impact on B16 cell-derived oncogenesis in vitro and in vivo using syngeneic C57BL/6 and Pcsk9^-/- mice. Intratumoral accumulation of cholesterol was determined. RNA-seq was performed on individual tumor types. Differentially-expressed genes (DEGs) were derived from the comparisons of B16 PCSK9, B16 D374Y, or B16 Q152H tumors to B16 EV allografts and analysed for pathway alterations.

RESULTS

PCSK9 expression and its network negatively correlated with the survival probability of patients with melanoma. PCSK9 promoted B16 cell proliferation, migration, and growth in soft agar in vitro, formation of tumors in C57BL/6 mice in vivo, and accumulation of intratumoral cholesterol in a manner reflecting its regulation of the low-density lipoprotein receptor (LDLR): Q152H, EV, PCSK9, and D374Y. Tumor-associated T cells, CD8 + T cells, and NK cells were significantly increased in D374Y tumors along with upregulations of multiple immune checkpoints, IFNγ, and 143 genes associated with T cell dysfunction. Overlap of 36 genes between the D374Y DEGs and the PCSK9 DEGs predicted poor prognosis of melanoma and resistance to immune checkpoint blockade (ICB) therapy. CYTH4, DENND1C, AOAH, TBC1D10C, EPSTI1, GIMAP7, and FASL (FAS ligand) were novel predictors of ICB therapy and displayed high level of correlations with multiple immune checkpoints in melanoma and across 30 human cancers. We observed FAS ligand being among the most robust biomarkers of ICB treatment and constructed two novel and effective multigene panels predicting response to ICB therapy. The profiles of allografts produced by B16 EV, PCSK9, D374Y, and Q152H remained comparable in C57BL/6 and Pcsk9^-/- mice.

CONCLUSIONS

Tumor-derived PCSK9 plays a critical role in melanoma pathogenesis. PCSK9's oncogenic actions are associated with intratumoral cholesterol accumulation. PCSK9 systemically affects the immune system, contributing to melanoma immune evasion. Novel biomarkers derived from the PCSK9-network effectively predicted ICB therapy responses.

Transcriptomic reveals the ferroptosis features of host response in a mouse model of Zika virus infection

Zika virus (ZIKV) is a neurotropic flavivirus. The outbreak of ZIKV in 2016 created a global health emergency. However, the underlying pathogenic mechanisms remain elusive. We investigated the host response features of in vivo replication in a mouse model of ZIKV infection, by performing a series of transcriptomic and bioinformatic analyses of ZIKV and mock-infected brain tissue. Tissue damage, inflammatory cells infiltration and high viral replication were observed in the brain tissue of ZIKV infected mice. RNA-Seq of the brain indicated the activation of ferroptosis pathways. Enrichment analysis of ferroptosis regulators revealed their involvement in pathways such as mineral absorption, fatty acid biosynthesis, fatty acid degradation, PPAR signaling pathway, peroxidase, and adipokinesine signalling pathway. We then identified 12 interacted hub ferroptosis regulators (CYBB, HMOX1, CP, SAT1, TF, SLC39A14, FTL, LPCAT3, FTH1, SLC3A2, TP53, and SLC40A1) that were related to the differential expression of CD8⁺ T cells, microglia and monocytes. CYBB, HMOX1, SALT, and SLAC40A1 were selected as potential biomarkers of ZIKV infection. Finally, we validated our results using RT-qPCR and outside available datasets. For the first time, we proposed a possible mechanism of ferroptosis in brain tissue infected by ZIKV in mice and identified the four key ferroptosis regulators.

Hippo Signaling Pathway Regulates Cancer Cell-Intrinsic MHC-II Expression

MHC-II is known to be mainly expressed on the surface of antigen-presenting cells. Evidence suggests MHC-II is also expressed by cancer cells and may be associated with better immunotherapy responses. However, the role and regulation of MHC-II in cancer cells remain unclear. In this study, we leveraged data mining and experimental validation to elucidate the regulation of MHC-II in cancer cells and its role in modulating the response to immunotherapy. We collated an extensive collection of omics data to examine cancer cell-intrinsic MHC-II expression and its association with immunotherapy outcomes. We then tested the functional relevance of cancer cell-intrinsic MHC-II expression using a syngeneic transplantation model. Finally, we performed data mining to identify pathways potentially involved in the regulation of MHC-II expression, and experimentally validated candidate regulators. Analyses of preimmunotherapy clinical samples in the CheckMate 064 trial revealed that cancer cell-intrinsic MHC-II protein was positively correlated with more favorable immunotherapy outcomes. Comprehensive meta-analyses of multiomics data from an exhaustive collection of data revealed that MHC-II is heterogeneously expressed in various solid tumors, and its expression is particularly high in melanoma. Using a syngeneic transplantation model, we further established that melanoma cells with high MHC-II responded better to anti-PD-1 treatment. Data mining followed by experimental validation revealed the Hippo signaling pathway as a potential regulator of melanoma MHC-II expression. In summary, we identified the Hippo signaling pathway as a novel regulator of cancer cell-intrinsic MHC-II expression. These findings suggest modulation of MHC-II in melanoma could potentially improve immunotherapy response.

Positive regulators of T cell proliferation as biomarkers for predicting prognosis and characterizing the immune landscape in lung adenocarcinoma

Lung adenocarcinoma (LUAD) is the one of the most prevalent and fatal form of malignant tumors worldwide. Recently, immunotherapy is widely used in the treatment of patients with LUAD and has proved to be clinically effective in improve the prognosis of patients. But there still has been a tremendous thrust to further improve the efficacy of immunotherapy in individual patients with LUAD. The suppression of T cells and their effector functions in the tumor microenvironment (TME) of LUAD is one of the primary reasons for the low efficacy of immunotherapy in some patients with LUAD. Therefore, identifying positive regulators of T cell proliferation (TPRs) may offer novel avenues for LUAD immunotherapy. In this study, we comprehensively evaluated the infiltration patterns of TPRs in 1,066 patients with LUAD using unsupervised consensus clustering and identified correlations with genomic and clinicopathological characteristics. Three infiltrating TPR clusters were defined, and a TPR-related risk signature composed of nine TPRs was constructed using least absolute shrinkage and selection operator-Cox regression algorithms to classify the individual TPR infiltration patterns. Cluster 1 exhibited high levels of T cell infiltration and activation of immune-related signaling pathways, whereas cluster 2 was characterized by robust T cell immune infiltration and enrichment of pathways associated with carcinogenic gene sets and tumor immunity. Cluster 3 was characterized as an immune-desert phenotype. Moreover, the TPR signature was confirmed as an independent prognostic biomarker for drug sensitivity in patients with LUAD. In conclusion, the TPR signature may serve as a novel tool for effectively characterizing immune characteristics and evaluating the prognosis of patients with LUAD.

Sexually dimorphic activation of innate antitumor immunity prevents adrenocortical carcinoma development

Unlike most cancers, adrenocortical carcinomas (ACCs) are more frequent in women than in men, but the underlying mechanisms of this sexual dimorphism remain elusive. Here, we show that inactivation of Znrf3 in the mouse adrenal cortex, recapitulating the most frequent alteration in ACC patients, is associated with sexually dimorphic tumor progression. Although female knockouts develop metastatic carcinomas at 18 months, adrenal hyperplasia regresses in male knockouts. This male-specific phenotype is associated with androgen-dependent induction of senescence, recruitment, and differentiation of highly phagocytic macrophages that clear out senescent cells. In contrast, in females, macrophage recruitment is delayed and dampened, which allows for aggressive tumor progression. Consistently, analysis of TCGA-ACC data shows that phagocytic macrophages are more prominent in men and are associated with better prognosis. Together, these data show that phagocytic macrophages are key players in the sexual dimorphism of ACC that could be previously unidentified allies in the fight against this devastating cancer.

Comprehensive analysis for cellular senescence-related immunogenic characteristics and immunotherapy prediction of acute myeloid leukemia

In malignancies, cellular senescence is critical for carcinogenesis, development, and immunological regulation. Patients with acute myeloid leukemia (AML) have not investigated a reliable cellular senescence-associated profile and its significance in outcomes and therapeutic response. Cellular senescence-related genes were acquired from the CellAge database, while AML data were obtained from the GEO and TCGA databases. The TCGA-AML group served as a training set to construct a prognostic risk score signature, while the GSE71014 set was used as a testing set to validate the accuracy of the signature. Through exploring the expression profiles of cellular senescence-related genes (SRGs) in AML patients, we used Lasso and Cox regression analysis to establish the SRG-based signature (SRGS), which was validated as an independent prognostic predictor for AML patients via clinical correlation. Survival analysis showed that AML patients in the low-risk score group had a longer survival time. Tumor immune infiltration and functional enrichment analysis demonstrated that AML patients with low-risk scores had higher immune infiltration and active immune-related pathways. Meanwhile, drug sensitivity analysis and the TIDE algorithm showed that the low-risk score group was more susceptible to chemotherapy and immunotherapy. Cell line analysis in vitro further confirmed that the SRGs in the proposed signature played roles in the susceptibility to cytarabine and YM155. Our results indicated that SRGS, which regulates the immunological microenvironment, is a reliable predictor of the clinical outcome and immunotherapeutic response in AML.

SimBu: bias-aware simulation of bulk RNA-seq data with variable cell-type composition

MOTIVATION

As complex tissues are typically composed of various cell types, deconvolution tools have been developed to computationally infer their cellular composition from bulk RNA sequencing (RNA-seq) data. To comprehensively assess deconvolution performance, gold-standard datasets are indispensable. Gold-standard, experimental techniques like flow cytometry or immunohistochemistry are resource-intensive and cannot be systematically applied to the numerous cell types and tissues profiled with high-throughput transcriptomics. The simulation of 'pseudo-bulk' data, generated by aggregating single-cell RNA-seq expression profiles in pre-defined proportions, offers a scalable and cost-effective alternative. This makes it feasible to create in silico gold standards that allow fine-grained control of cell-type fractions not conceivable in an experimental setup. However, at present, no simulation software for generating pseudo-bulk RNA-seq data exists.

RESULTS

We developed SimBu, an R package capable of simulating pseudo-bulk samples based on various simulation scenarios, designed to test specific features of deconvolution methods. A unique feature of SimBu is the modeling of cell-type-specific mRNA bias using experimentally derived or data-driven scaling factors. Here, we show that SimBu can generate realistic pseudo-bulk data, recapitulating the biological and statistical features of real RNA-seq data. Finally, we illustrate the impact of mRNA bias on the evaluation of deconvolution tools and provide recommendations for the selection of suitable methods for estimating mRNA content. SimBu is a user-friendly and flexible tool for simulating realistic pseudo-bulk RNA-seq datasets serving as in silico gold-standard for assessing cell-type deconvolution methods.

AVAILABILITY AND IMPLEMENTATION

SimBu is freely available at https://github.com/omnideconv/SimBu as an R package under the GPL-3 license.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Identification of key monocytes/macrophages related gene set of the early-stage abdominal aortic aneurysm by integrated bioinformatics analysis and experimental validation

Objective

Abdominal aortic aneurysm (AAA) is a lethal peripheral vascular disease. Inflammatory immune cell infiltration is a central part of the pathogenesis of AAA. It’s critical to investigate the molecular mechanisms underlying immune infiltration in early-stage AAA and look for a viable AAA marker.

Methods

In this study, we download several mRNA expression datasets and scRNA-seq datasets of the early-stage AAA models from the NCBI-GEO database. mMCP-counter and CIBERSORT were used to assess immune infiltration in early-stage experimental AAA. The scRNA-seq datasets were then utilized to analyze AAA-related gene modules of monocytes/macrophages infiltrated into the early-stage AAA by Weighted Correlation Network analysis (WGCNA). After that, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) functional enrichment analysis for the module genes was performed by ClusterProfiler. The STRING database was used to create the protein-protein interaction (PPI) network. The Differentially Expressed Genes (DEGs) of the monocytes/macrophages were explored by Limma-Voom and the key gene set were identified. Then We further examined the expression of key genes in the human AAA dataset and built a logistic diagnostic model for distinguishing AAA patients and healthy people. Finally, real-time quantitative polymerase chain reaction (RT-qPCR) and Enzyme Linked Immunosorbent Assay (ELISA) were performed to validate the gene expression and serum protein level between the AAA and healthy donor samples in our cohort.

Results

Monocytes/macrophages were identified as the major immune cells infiltrating the early-stage experimental AAA. After pseudocell construction of monocytes/macrophages from scRNA-seq datasets and WGCNA analysis, four gene modules from two datasets were identified positively related to AAA, mainly enriched in Myeloid Leukocyte Migration, Collagen-Containing Extracellular matrix, and PI3K-Akt signaling pathway by functional enrichment analysis. Thbs1, Clec4e, and Il1b were identified as key genes among the hub genes in the modules, and the high expression of Clec4e, Il1b, and Thbs1 was confirmed in the other datasets. Then, in human AAA transcriptome datasets, the high expression of CLEC4E, IL1B was confirmed and a logistic regression model based on the two gene expressions was built, with an AUC of 0.9 in the train set and 0.79 in the validated set. Additionally, in our cohort, we confirmed the increased serum protein levels of IL-1β and CLEC4E in AAA patients as well as the increased expression of these two genes in AAA aorta samples.

Conclusion

This study identified monocytes/macrophages as the main immune cells infiltrated into the early-stage AAA and constructed a logistic regression model based on monocytes/macrophages related gene set. This study could aid in the early diagnostic of AAA.

Loss of MHC-I antigen presentation correlated with immune checkpoint blockade tolerance in MAPK inhibitor-resistant melanoma

Immune checkpoint blockade and MAPK-targeted combined therapy is a promising regimen for advanced melanoma patients. However, the clinical benefit from this combo regimen remains limited, especially in patients who acquired resistance to MAPK-targeted therapy. Here, we systematically characterized the immune landscape during MAPK-targeted therapy in patients and mouse melanoma models. We observed that both the abundance of tumor-infiltrated T cells and the expression of immune-related genes were upregulated in the drug-responsive period, but downregulated in the resistance period, implying that acquired drug resistance dampens the antitumor immune response. Further transcriptomic dissection indicated that loss of MHC-I antigen presentation on tumor cells plays a critical role in the reduction of T cell infiltration during drug resistance. Survival analysis demonstrates that loss of antigen presentation and reduction of T-cell infiltration during acquired drug resistance are associated with poorer clinical response and prognosis of anti-PD-1 therapy in melanoma patients. In addition, we identified that alterations in the MAPK inhibitor resistance-related oncogenic signaling pathway closely correlated with deficiency of MHC-I antigen presentation, including activation of the PI3K-mTOR, MAPK, and Wnt pathways. In conclusion, our research illuminates that decreased infiltration of T cells is associated with acquired drug resistance during MAPK-targeted therapy, which may underlie the cross-resistance to immune checkpoint blockade.

Pan-Cancer Analysis of the Oncogenic and Immunological Role of Solute Carrier Family 6 Member 8 (SLC6A8)

There is mounting evidence on the implication of SLC6A8 in the initiation and progression of human cancers. However, a comprehensive understanding of the role of SLC6A8 in pan-cancer remains elusive yet. Bioinformatics analysis was performed to investigate the expression and mutation profiles of SLC6A8 in cancers, and the association of SLC6A8 expression with cancer patients’ survival and immune cell infiltration. In general, SLC6A8 is significantly upregulated across multiple cancers. SLC6A8 expression was inconsistently prognostic in different types of cancer, albeit associated with favorable survival in the vast majority of cancers. The receiver operating characteristic curves showed that SLC6A8 was relatively accurate in identifying possible cancer patients. The genetic alteration of SLC6A8, including mutation, amplification and deletion, was frequently present across various types of cancer. Mechanistically, SLC6A8 might be involved in tumorigenesis through “carbon metabolism” and “HIF-1 signaling pathway.” Besides, SLC6A8 expression had significant correlation with immune checkpoints genes and tumor-infiltrating immune cell abundances. The present study offers a novel insight into the roles of SLC6A8 in the oncogenesis and development of multiple common cancers.

GM-CSF drives myelopoiesis, recruitment and polarisation of tumour-associated macrophages in cholangiocarcinoma and systemic blockade facilitates antitumour immunity

OBJECTIVE

Intrahepatic cholangiocarcinoma (iCCA) is rising in incidence, and at present, there are limited effective systemic therapies. iCCA tumours are infiltrated by stromal cells, with high prevalence of suppressive myeloid populations including tumour-associated macrophages (TAMs) and myeloid-derived suppressor cells (MDSCs). Here, we show that tumour-derived granulocyte-macrophage colony-stimulating factor (GM-CSF) and the host bone marrow is central for monopoiesis and potentiation of TAMs, and abrogation of this signalling axis facilitates antitumour immunity in a novel model of iCCA.

METHODS

Blood and tumours were analysed from iCCA patients and controls. Treatment and correlative studies were performed in mice with autochthonous and established orthotopic iCCA tumours treated with anti-GM-CSF monoclonal antibody.

RESULTS

Systemic elevation in circulating myeloid cells correlates with poor prognosis in patients with iCCA, and patients who undergo resection have a worse overall survival if tumours are more infiltrated with CD68+ TAMs. Mice with spontaneous iCCA demonstrate significant elevation of monocytic myeloid cells in the tumour microenvironment and immune compartments, and tumours overexpress GM-CSF. Blockade of GM-CSF with a monoclonal antibody decreased tumour growth and spread. Mice bearing orthotopic tumours treated with anti-GM-CSF demonstrate repolarisation of immunosuppressive TAMs and MDSCs, facilitating T cell response and tumour regression. GM-CSF blockade dampened inflammatory gene networks in tumours and TAMs. Human tumours with decreased GM-CSF expression exhibit improved overall survival after resection.

CONCLUSIONS

iCCA uses the GM-CSF-bone marrow axis to establish an immunosuppressive tumour microenvironment. Blockade of the GM-CSF axis promotes antitumour T cell immunity.

Cancer Stem Cell-Associated Immune Microenvironment in Recurrent Glioblastomas

Glioblastoma multiforme (GBM) is the most incurable tumor (due to the difficulty in complete surgical resection and the resistance to conventional chemo/radiotherapies) that displays a high relapse frequency. Cancer stem cells (CSCs) have been considered as a promising target responsible for therapy resistance and cancer recurrence. CSCs are known to organize a self-advantageous microenvironment (niche) for their maintenance and expansion. Therefore, understanding how the microenvironment is reconstructed by the remaining CSCs after conventional treatments and how it eventually causes recurrence should be essential to inhibit cancer recurrence. However, the number of studies focusing on recurrence is limited, particularly those related to tumor immune microenvironment, while numerous data have been obtained from primary resected samples. Here, we summarize recent investigations on the immune microenvironment from the viewpoint of recurrent GBM (rGBM). Based on the recurrence-associated immune cell composition reported so far, we will discuss how CSCs manipulate host immunity and create the special microenvironment for themselves to regrow. An integrated understanding of the interactions between CSCs and host immune cells at the recurrent phase will lead us to develop innovative therapies and diagnoses to achieve GBM eradication.

Extracellular vimentin mimics VEGF and is a target for anti-angiogenic immunotherapy

Anti-angiogenic cancer therapies possess immune-stimulatory properties by counteracting pro-angiogenic molecular mechanisms. We report that tumor endothelial cells ubiquitously overexpress and secrete the intermediate filament protein vimentin through type III unconventional secretion mechanisms. Extracellular vimentin is pro-angiogenic and functionally mimics VEGF action, while concomitantly acting as inhibitor of leukocyte-endothelial interactions. Antibody targeting of extracellular vimentin shows inhibition of angiogenesis in vitro and in vivo. Effective and safe inhibition of angiogenesis and tumor growth in several preclinical and clinical studies is demonstrated using a vaccination strategy against extracellular vimentin. Targeting vimentin induces a pro-inflammatory condition in the tumor, exemplified by induction of the endothelial adhesion molecule ICAM1, suppression of PD-L1, and altered immune cell profiles. Our findings show that extracellular vimentin contributes to immune suppression and functions as a vascular immune checkpoint molecule. Targeting of extracellular vimentin presents therefore an anti-angiogenic immunotherapy strategy against cancer.

The pro-tumorigenic effects of vimentin have been attributed to intracellular functions in tumour cells so far. Here, the authors show that tumour endothelial cells can secrete vimentin as a pro-angiogenic factor and that targeting of vimentin can be used as an immunotherapeutic strategy.

Increased Expression of Homeobox 5 Predicts Poor Prognosis: A Potential Prognostic Biomarker for Glioma

Background

The homeobox gene 5 (HOXB5) encodes a transcription factor that regulates the embryonic development of the central nervous system. Notably, its expression pattern and prognostic role in glioma remain unelucidated.

Methods

This study identified the relationship between HOXB5 and glioma by investigating HOXB5 expression data from The Cancer Genome Atlas and the Genotype Tissue Expression databases and validating the obtained data using the Chinese Glioma Genome Atlas database. Western blots were used to identify HOXB5 expression levels in glioma cells and clinical samples. Kaplan-Meier and multivariate Cox regression analyses were performed to assess the prognostic value of HOXB5. The key functions and signaling pathways related to HOXB5 were analyzed using GO, KEGG, and GSEA. Immune infiltration was calculated using the microenvironment cell populations-counter, estimate the proportion of immune and cancer, and ESTIMATE algorithms.

Results

The expression of HOXB5 was upregulated in glioma and generally increased with malignancy. HOXB5 was an independent prognostic factor for glioma patients. A nomogram was further built that integrated HOXB5, and it showed stratifying prediction accuracy and efficiency. HOXB5 was associated with the regulation of cell growth, endothelial cell growth, and the IL-6/JAK-STAT3 pathway, and was determined to possibly promote stomatal specimen enrichment and angiogenesis.

Conclusion

HOXB5 protein is overexpressed in glioma and might serve as a good predictive factor of this disease.

Molecular Subtyping Resource: a user-friendly tool for rapid biological discovery from transcriptional data

Generation of transcriptional data has dramatically increased in the past decade, driving the development of analytical algorithms that enable interrogation of the biology underpinning the profiled samples. However, these resources require users to have expertise in data wrangling and analytics, reducing opportunities for biological discovery by 'wet-lab' users with a limited programming skillset. Although commercial solutions exist, costs for software access can be prohibitive for academic research groups. To address these challenges, we have developed an open source and user-friendly data analysis platform for on-the-fly bioinformatic interrogation of transcriptional data derived from human or mouse tissue, called Molecular Subtyping Resource (MouSR). This internet-accessible analytical tool, https://mousr.qub.ac.uk/, enables users to easily interrogate their data using an intuitive 'point-and-click' interface, which includes a suite of molecular characterisation options including quality control, differential gene expression, gene set enrichment and microenvironmental cell population analyses from RNA sequencing. The MouSR online tool provides a unique freely available option for users to perform rapid transcriptomic analyses and comprehensive interrogation of the signalling underpinning transcriptional datasets, which alleviates a major bottleneck for biological discovery. This article has an associated First Person interview with the first author of the paper.

Identification of Diagnostic Biomarkers Associated with Stromal and Immune Cell Infiltration in Fatty Infiltration After Rotator Cuff Tear by Integrating Bioinformatic Analysis and Machine-Learning

Purpose

The present study aimed to explore potential diagnostic biomarkers for fatty infiltration (FI) of the rotator cuff muscles after rotator cuff tear (RCT) and investigate the influence of stromal and immune cell infiltration on this pathology.

Methods

The GSE130447 and GSE103266 datasets were downloaded from the Gene Expression Omnibus database. Differentially expressed genes (DEGs) were identified, and gene set enrichment analyses were performed by R software. Two machine learning algorithms, random forest and multiple support vector machine recursive feature elimination (mSVM-RFE), were used to screen candidate biomarkers. The diagnostic value of the screened biomarkers was further validated by the area under the ROC curve (AUC) in the GSE103266 dataset. Murine microenvironment cell population counter (mMCP-counter) method was employed to estimate stromal and immune cell infiltration of FI. The correlation between biomarkers and infiltrated immune and stromal cell subsets was further analyzed.

Results

A total of 2123 DEGs were identified. The identified DEGs were predominantly linked to immune system process, extracellular matrix organization and PPAR signalling pathway. FABP5 (AUC = 0.958) and MGP (AUC = 1) were screened as diagnostic biomarkers of FI. Stromal and immune cell infiltration analysis showed that monocytes, mast cells, vessels, endothelial cells and fibroblasts may be related to the process of FI. FABP5 and MGP were positively correlated with vessels whereas negatively correlated with monocytes and mast cells.

Conclusion

FABP5 and MGP can serve as diagnostic biomarkers of FI after RCT, and stromal and immune cell infiltration may play a crucial role in this pathology.

Network Pharmacology and Inflammatory Microenvironment Strategy Approach to Finding the Potential Target of Siraitia grosvenorii (Luo Han Guo) for Glioblastoma

Background: Glioblastoma (GBM) is the most common and aggressive primary intracranial tumor of the central nervous system, and the prognosis of GBM remains a challenge using the standard methods of treatment—TMZ, radiation, and surgical resection. Traditional Chinese medicine (TCM) is a helpful complementary and alternative medicine. However, there are relatively few studies on TCM for GBM.

Purpose: We aimed to find the connection between TCM and anti-GBM.

Study design: Network pharmacology and inflammatory microenvironment strategy were used to predict Siraitia grosvenorii (Luo Han Guo) target for treating glioblastoma.

Methods: We mainly used network pharmacology and bioinformatics.

Results: CCL5 was significantly highly expressed in GBM with poor prognostics. Uni-cox and randomForest were used to determine that CCL5 was especially a biomarker in GBM. CCL5 was also the target for SG and TMZ. The active ingredient of Luo Han Guo — squalene and CCL5 —showed high binding efficiency. CCL5, a chemotactic ligand, was enriched and positively correlated in eosinophils. CCL5 was also the target of Luo Han Guo, and its effective active integrate compound –— squalene — might act on CCL5.

Conclusion: SG might be a new complementary therapy of the same medicine and food, working on the target CCL5 and playing an anti-GBM effect. CCL5 might affect the immune microenvironment of GBM.