The C-type lectin DCIR contributes to the immune response and pathogenesis of colorectal cancer

Development and progression of malignancies are accompanied and influenced by alterations in the surrounding immune microenvironment. Understanding the cellular and molecular interactions between immune cells and cancer cells has not only provided important fundamental insights into the disease, but has also led to the development of new immunotherapies. The C-type lectin Dendritic Cell ImmunoReceptor (DCIR) is primarily expressed by myeloid cells and is an important regulator of immune homeostasis, as demonstrated in various autoimmune, infectious and inflammatory contexts. Yet, the impact of DCIR on cancer development remains largely unknown. Analysis of available transcriptomic data of colorectal cancer (CRC) patients revealed that high DCIR gene expression is associated with improved patients' survival, immunologically "hot" tumors and high immunologic constant of rejection, thus arguing for a protective and immunoregulatory role of DCIR in CRC. In line with these correlative data, we found that deficiency of DCIR1, the murine homologue of human DCIR, leads to the development of significantly larger tumors in an orthotopic murine model of CRC. This phenotype is accompanied by an altered phenotype of tumor-associated macrophages (TAMs) and a reduction in the percentage of activated effector CD4+ and CD8+ T cells in CRC tumors of DCIR1-deficient mice. Overall, our results show that DCIR promotes antitumor immunity in CRC, making it an attractive target for the future development of immunotherapies to fight the second deadliest cancer in the world.

Location matters: spatial dynamics of tumor-infiltrating T cell subsets is prognostic in colon cancer

Background

Colon cancer is a heterogeneous disease and consists of various molecular subtypes. Despite advances in high-throughput expression profiling, limitations remain in predicting clinical outcome and assigning specific treatment to individual cases. Tumor-immune interactions play a critical role, with tumors that activate the immune system having better outcome for the patient. The localization of T cells within tumor epithelium, to enable direct contact, is essential for antitumor function, but bulk DNA/RNA sequencing data lacks spatial distribution information. In this study, we provide spatial T cell tumor distribution and connect these data with previously determined genomic data in the AC-ICAM colon cancer patient cohort.

Methods

Colon cancer patients (n=90) with transcriptome data available were selected. We used a custom multiplex immunofluorescence assay on colon tumor tissue sections for quantifying T cell subsets spatial distribution in the tumor microenvironment, in terms of cell number, location, mutual distance, and distance to tumor cells. Statistical analyses included the previously determined Immunologic Constant of Rejection (ICR) transcriptome correlation and patient survival, revealing potential prognostic value in T cell spatial distribution.

Results

T cell phenotypes were characterized and CD3+CD8-FoxP3- T cells were found to be the predominant tumor-infiltrating subtype while CD3+FoxP3+ T cells and CD3+CD8+ T cells showed similar densities. Spatial distribution analysis elucidated that proliferative T cells, characterized by Ki67 expression, and Granzyme B-expressing T cells were predominantly located within the tumor epithelium. We demonstrated an increase in immune cell density and a decrease in the distance of CD3+CD8+ T cells to the nearest tumor cell, in the immune active, ICR High, immune subtypes. Higher densities of stromal CD3+FoxP3+ T cells showed enhanced survival outcomes, and patients exhibited superior clinical benefits when greater spatial distances were observed between CD3+CD8-FoxP3- or CD3+CD8+ T cells and CD3+FoxP3+ T cells.

Conclusion

Our study's in-depth analysis of the spatial distribution and densities of major T cell subtypes within the tumor microenvironment has provided valuable information that paves the way for further research into the intricate relationships between immune cells and colon cancer development.

Organizing training workshops on gene literature retrieval, profiling, and visualization for early career researchers

Early-career researchers must acquire the skills necessary to effectively search and extract information from biomedical literature. This ability is for instance crucial for evaluating the novelty of experimental results, and assessing potential publishing opportunities. Given the rapidly growing volume of publications in the field of biomedical research, new systematic approaches need to be devised and adopted for the retrieval and curation of literature relevant to a specific theme. In this context, we present a hands-on training curriculum aimed at retrieval, profiling, and visualization of literature associated with a given topic. The curriculum was implemented in a workshop in January 2021. Here we provide supporting material and step-by-step implementation guidelines with the ISG15 gene literature serving as an illustrative use case. Workshop participants can learn several skills, including: 1) building and troubleshoot PubMed queries in order to retrieve the literature associated with a gene of interest; 2) identifying key concepts relevant to given themes (such as cell types, diseases, and biological processes); 3) measuring the prevalence of these concepts in the gene literature; 4) extracting key information from relevant articles, and 5) developing a background section or summary on the basis of this information. Finally, trainees can learn to consolidate the structured information captured through this process for presentation via an interactive web application.

An integrated tumor, immune and microbiome atlas of colon cancer

The lack of multi-omics cancer datasets with extensive follow-up information hinders the identification of accurate biomarkers of clinical outcome. In this cohort study, we performed comprehensive genomic analyses on fresh-frozen samples from 348 patients affected by primary colon cancer, encompassing RNA, whole-exome, deep T cell receptor and 16S bacterial rRNA gene sequencing on tumor and matched healthy colon tissue, complemented with tumor whole-genome sequencing for further microbiome characterization. A type 1 helper T cell, cytotoxic, gene expression signature, called Immunologic Constant of Rejection, captured the presence of clonally expanded, tumor-enriched T cell clones and outperformed conventional prognostic molecular biomarkers, such as the consensus molecular subtype and the microsatellite instability classifications. Quantification of genetic immunoediting, defined as a lower number of neoantigens than expected, further refined its prognostic value. We identified a microbiome signature, driven by Ruminococcus bromii, associated with a favorable outcome. By combining microbiome signature and Immunologic Constant of Rejection, we developed and validated a composite score (mICRoScore), which identifies a group of patients with excellent survival probability. The publicly available multi-omics dataset provides a resource for better understanding colon cancer biology that could facilitate the discovery of personalized therapeutic approaches.

Analytic pipelines to assess the relationship between immune response and germline genetics in human tumors

Germline genetic variants modulate human immune response. We present analytical pipelines for assessing the contribution of hosts' genetic background to the immune landscape of solid tumors using harmonized data from more than 9,000 patients in The Cancer Genome Atlas (TCGA). These include protocols for heritability, genome-wide association studies (GWAS), colocalization, and rare variant analyses. These workflows are developed around the structure of TCGA but can be adapted to explore other repositories or in the context of cancer immunotherapy. For complete details on the use and execution of this protocol, please refer to Sayaman et al. (2021).

The immune landscape of solid pediatric tumors

Background

Large immunogenomic analyses have demonstrated the prognostic role of the functional orientation of the tumor microenvironment in adult solid tumors, this variable has been poorly explored in the pediatric counterpart.

Methods

We performed a systematic analysis of public RNAseq data (TARGET) for five pediatric tumor types (408 patients): Wilms tumor (WLM), neuroblastoma (NBL), osteosarcoma (OS), clear cell sarcoma of the kidney (CCSK) and rhabdoid tumor of the kidney (RT). We assessed the performance of the Immunologic Constant of Rejection (ICR), which captures an active Th1/cytotoxic response. We also performed gene set enrichment analysis (ssGSEA) and clustered more than 100 well characterized immune traits to define immune subtypes and compared their outcome.

Results

A higher ICR score was associated with better survival in OS and high risk NBL without MYCN amplification but with poorer survival in WLM. Clustering of immune traits revealed the same five principal modules previously described in adult tumors (TCGA). These modules divided pediatric patients into six immune subtypes (S1-S6) with distinct survival outcomes. The S2 cluster showed the best overall survival, characterized by low enrichment of the wound healing signature, high Th1, and low Th2 infiltration, while the reverse was observed in S4. Upregulation of the WNT/Beta-catenin pathway was associated with unfavorable outcomes and decreased T-cell infiltration in OS.

Conclusions

We demonstrated that extracranial pediatric tumors could be classified according to their immune disposition, unveiling similarities with adults’ tumors. Immunological parameters might be explored to refine diagnostic and prognostic biomarkers and to identify potential immune-responsive tumors.

Transcriptomic and immunophenotypic profiling reveals molecular and immunological hallmarks of colorectal cancer tumourigenesis

OBJECTIVE

Biological insights into the stepwise development and progression of colorectal cancer (CRC) are imperative to develop tailored approaches for early detection and optimal clinical management of this disease. Here, we aimed to dissect the transcriptional and immunologic alterations that accompany malignant transformation in CRC and to identify clinically relevant biomarkers through spatial profiling of pT1 CRC samples.

DESIGN

We employed digital spatial profiling (GeoMx) on eight pT1 CRCs to study gene expression in the epithelial and stromal segments across regions of distinct histology, including normal mucosa, low-grade and high-grade dysplasia and cancer. Consecutive histology sections were profiled by imaging mass cytometry to reveal immune contextures. Finally, publicly available single-cell RNA-sequencing data was analysed to determine the cellular origin of relevant transcripts.

RESULTS

Comparison of gene expression between regions within pT1 CRC samples identified differentially expressed genes in the epithelium (n=1394 genes) and the stromal segments (n=1145 genes) across distinct histologies. Pathway analysis identified an early onset of inflammatory responses during malignant transformation, typified by upregulation of gene signatures such as innate immune sensing. We detected increased infiltration of myeloid cells and a shift in macrophage populations from pro-inflammatory HLA-DR⁺CD204^- macrophages to HLA-DR^-CD204⁺ immune-suppressive subsets from normal tissue through dysplasia to cancer, accompanied by the upregulation of the CD47/SIRPα 'don't eat me signal'.

CONCLUSION

Spatial profiling revealed the molecular and immunological landscape of CRC tumourigenesis at early disease stage. We identified biomarkers with strong association with disease progression as well as targetable immune processes that are exploitable in a clinical setting.

Correlation of Immunological and Histopathological Features with Gene Expression-Based Classifiers in Colon Cancer Patients

The purpose of this study was to evaluate the association between four distinct histopathological features: (1) tumor infiltrating lymphocytes, (2) mucinous differentiation, (3) tumor-stroma ratio, plus (4) tumor budding and two gene expression-based classifiers—(1) consensus molecular subtypes (CMS) plus (2) colorectal cancer intrinsic subtypes (CRIS). All four histopathological features were retrospectively scored on hematoxylin and eosin sections of the most invasive part of the primary tumor in 218 stage II and III colon cancer patients from two independent cohorts (AMC-AJCC-90 and AC-ICAM). RNA-based CMS and CRIS assignments were independently obtained for all patients. Contingency tables were constructed and a χ2 test was used to test for statistical significance. Odds ratios with 95% confidence intervals were calculated. The presence of tumor infiltrating lymphocytes and a mucinous phenotype (>50% mucinous surface area) were strongly correlated with CMS1 (p < 0.001 and p = 0.008) and CRIS-A (p = 0.006 and p < 0.001). The presence of mucus (≥ 10%) was associated with CMS3: mucus was present in 64.1% of all CMS3 tumors (p < 0.001). Although a clear association between tumor-stroma ratio and CMS4 was established in this study (p = 0.006), still 32 out of 61 (52.5%) CMS4 tumors were scored as stroma-low, indicating that CMS4 tumors cannot be identified solely based on stromal content. Higher budding counts were seen in CMS4 and CRIS-B tumors (p = 0.045 and p = 0.046). No other associations of the measured parameters were seen for any of the other CRIS subtypes. Our analysis revealed clear associations between histopathologic features and CMS or CRIS subtypes. However, identification of distinct molecular subtypes solely based on histopathology proved to be infeasible. Combining both molecular and morphologic features could potentially improve patient stratification.

Transcriptomic profile investigations highlight a putative role for NUDT16 in sepsis

Sepsis is an aberrant systemic inflammatory response mediated by the acute activation of the innate immune system. Neutrophils are important contributors to the innate immune response that controls the infection, but harbour the risk of collateral tissue damage such as thrombosis and organ dysfunction. A better understanding of the modulations of cellular processes in neutrophils and other blood cells during sepsis is needed and can be initiated via transcriptomic profile investigations. To that point, the growing repertoire of publicly accessible transcriptomic datasets serves as a valuable resource for discovering and/or assessing the robustness of biomarkers. We employed systematic literature mining, reductionist approach to gene expression profile and empirical in vitro work to highlight the role of a Nudix hydrolase family member, NUDT16, in sepsis. The relevance and implication of the expression of NUDT16 under septic conditions and the putative functional roles of this enzyme are discussed.

673 (Re-) Solving the biology of colorectal cancer onset and progression to improve treatment and prevention

Background

Colorectal cancer (CRC) development is accompanied by the gradual accumulation of genetic alterations in epithelial cells of the colon and rectum.1 2 The paradigm of the adenoma-carcinoma sequence was originally centered around cancer cells; however, it is now clear that the tumor microenvironment plays a substantial role in cancer progression and patient outcome.3 In recent years, technologies have evolved rapidly, allowing the multiplexed quantification of gene expression while preserving spatial context.4 Furthermore, some spatial transcriptomic technologies also allow the parallel interrogation of different cell populations in the tumor microenvironment. Here, we performed digital spatial profiling on early-stage CRC samples to elucidate the biological processes that are at the basis of malignant transformation and to identify novel therapeutic targets and (immune) biomarkers.

Methods

Endoscopically resected early-stage CRC samples were obtained at Leiden University Medical Center. In total, 144 areas of illumination were interrogated with GeoMx digital spatial profiling using the Cancer Transcriptome Atlas (>1,800 genes). In each of eight samples, nine regions of interest with different levels of cancer progression were selected, including normal epithelium, transition areas, low-, and high-grade dysplasia, and invasive carcinoma (figure 1A). We segmented each region based on cytokeratin and vimentin protein expression (figure 1B). Immunohistochemical detection was performed on these samples and 26 additional samples to validate targets associated with disease progression.

Results

Digital spatial profiling allowed us to dissect transcriptional alterations in epithelial and stromal fractions between different regions from healthy tissue, different degrees of dysplasia, and cancer. Gene expression data revealed a clear separation of profiled areas by histologic category. Interestingly, gene expression features in the stromal compartment provided a better data-driven separation of histologic categories than the epithelial fraction (figure 1C). Substantial changes in immune-related pathways were identified, including differential expression of specific immunomodulators. We validated the expression of several candidate biomarkers/targets that demonstrated consistent alterations from normal tissue to cancer by immunohistochemistry. Several proteins were identified that could clearly discriminate benign from malignant tissue.

Conclusions

We here demonstrated the unique biological insights that are provided by spatial examination of early-stage CRC by digital spatial profiling. We identified specific genes that were altered during CRC tumorigenesis, in epithelial and stromal/immune fractions. Furthermore, our results indicate an essential role for innate immunity in colorectal cancer onset and progression. The genes identified by this approach could potentially serve as novel biomarkers and targets for early interception or prevention of CRC development.

Acknowledgements

This work was supported by the European Research Council (ERC) Starting grant awarded to Dr. Noel F. de Miranda and the Stichting Management Apothekers en de Gezondheidszorg (STIMAG) Research grant awarded to Jessica Roelands.

Trial Registration

N/A

References

Fearon ER, Vogelstein B. A genetic model for colorectal tumorigenesis. Cell 1990;61(5):759–767. doi: 10.1016/0092-8674(90)90186-I.

Nowell PC. The clonal evolution of tumor cell populations. Science 1976;194(4260):23–28. doi: 10.1126/science.959840.

Hanahan D, Weinberg RA. Hallmarks of cancer: the next generation. Cell 2011;144(5):646–674. doi: 10.1016/j.cell.2011.02.013.

Merritt CR, et al. Multiplex digital spatial profiling of proteins and RNA in fixed tissue. Nat Biotechnol 2020;38(5):586–599. doi: 10.1038/s41587-020-0472-9.

Ethics Approval

This study was approved by the METC Leiden-Den Haag-Delft (protocol B20.039). Patient samples were anonymised and handled according to the medical ethical guidelines described in the Code of Conduct for Proper Secondary Use of Human Tissue of the Dutch Federation of Biomedical Scientific Societies.

Abstract 673 Figure 1

Transcriptional alterations in early-stage colorectal cancer. Digital spatial profiling defines transcriptional alterations in early-stage colorectal cancer. (A) Schematic representation of an early-stage CRC sample containing regions with different levels of cancer progression, including normal epithelium, transition areas, low-, and high-grade dysplasia, and invasive carcinoma. (B) Segmentation based immunofluorescent labelling with antibodies directed against PanCK and Vimentin in one of the early-stage CRC samples. Artificial overlay of implemented segmentation is indicated for each ROI, visualizing Vimentin+ (pink) and PanCK+ (orange) segments. Inset: higher magnification of an individual ROI. (C) Dimension reduction of expression of all quantified genes by t-Distributed Stochastic Neighbor Embedding (tSNE). tSNE plots are annotated by segment (left), and histological region (right).

Combined Assessment of the Tumor-Stroma Ratio and Tumor Immune Cell Infiltrate for Immune Checkpoint Inhibitor Therapy Response Prediction in Colon Cancer

The best current biomarker strategies for predicting response to immune checkpoint inhibitor (ICI) therapy fail to account for interpatient variability in response rates. The histologic tumor-stroma ratio (TSR) quantifies intratumoral stromal content and was recently found to be predictive of response to neoadjuvant therapy in multiple cancer types. In the current work, we predicted the likelihood of ICI therapy responsivity of 335 therapy-naive colon adenocarcinoma tumors from The Cancer Genome Atlas, using bioinformatics approaches. The TSR was scored on diagnostic tissue slides, and tumor-infiltrating immune cells (TIICs) were inferred from transcriptomic data. Tumors with high stromal content demonstrated increased T regulatory cell infiltration (p = 0.014) but failed to predict ICI therapy response. Consequently, we devised a hybrid tumor microenvironment classification of four stromal categories, based on histological stromal content and transcriptomic-deconvoluted immune cell infiltration, which was associated with previously established transcriptomic and genomic biomarkers for ICI therapy response. By integrating these biomarkers, stroma-low/immune-high tumors were predicted to be most responsive to ICI therapy. The framework described here provides evidence for expansion of current histological TIIC quantification to include the TSR as a novel, easy-to-use biomarker for the prediction of ICI therapy response.

Myeloid Cells Are Enriched in Tonsillar Crypts, Providing Insight into the Host Tropism of Human Papillomavirus

Viruses are the second leading cause of cancer worldwide, and human papillomavirus (HPV)-associated head and neck cancers are increasing in incidence in the United States. HPV preferentially infects the crypts of the tonsils rather than the surface epithelium. The present study sought to characterize the unique microenvironment within the crypts to better understand the host tropism of HPV to a lymphoid-rich organ. Laser-capture microdissection of distinct anatomic areas (crypts, surface epithelium, and germinal centers) of the tonsil, coupled with transcriptional analysis and multiparameter immunofluorescence staining, was performed and demonstrated that the tonsillar crypts are enriched with myeloid populations that co-express multiple canonical and noncanonical immune checkpoints, including PD-L1, CTLA-4, HAVCR2 (TIM-3), ADORA2A, IDO1, BTLA, LGALS3, CDH1, CEACAM1, PVR, and C10orf54 (VISTA). The resident monocytes may foster a permissive microenvironment that facilitates HPV infection and persistence. Furthermore, the myeloid populations within HPV-associated tonsil cancers co-express the same immune checkpoints, providing insight into potential novel immunotherapeutic targets for HPV-associated head and neck cancers.

Integrated transcriptional-phenotypic analysis captures systemic immunomodulation following antiangiogenic therapy in renal cell carcinoma patients

BACKGROUND

The combination of immune checkpoint blockade (ICB) with standard therapies is becoming a common approach for overcoming resistance to cancer immunotherapy in most human malignancies including metastatic renal cell carcinoma (mRCC). In this regard, insights into the immunomodulatory properties of antiangiogenic agents may help designing multidrug schedules based on specific immune synergisms.

METHODS

We used orthogonal transcriptomic and phenotyping platforms combined with functional analytic pipelines to elucidate the immunomodulatory effect of the antiangiogenic agent pazopanib in mRCC patients. Nine patients were studied longitudinally over a period of 6 months. We also analyzed transcriptional data from The Cancer Genome Atlas (TCGA) RCC cohort (N = 571) to assess the prognostic implications of our findings. The effect of pazopanib was assessed in vitro on NK cells and T cells. Additionally, myeloid-derived suppressor (MDSC)-like cells were generated from CD14+ monocytes transfected with mimics of miRNAs associated with MDSC function in the presence or absence of pazopanib.

RESULTS

Pazopanib administration caused a rapid and dramatic reshaping in terms of frequency and transcriptional activity of multiple blood immune cell subsets, with a downsizing of MDSC and regulatory T cells in favor of a strong enhancement in PD-1 expressing cytotoxic T and Natural Killer effectors. These changes were paired with an increase of the expression of transcripts reflecting activation of immune-effector functions. This immunomodulation was marked but transient, peaking at the third month of treatment. Moreover, the intratumoral expression level of a MDSC signature (MDSC INT) was strongly associated with poor prognosis in RCC patients. In vitro experiments indicate that the observed immunomodulation might be due to an inhibitory effect on MDSC-mediated suppression, rather than a direct effect on NK and T cells.

CONCLUSIONS

The marked but transient nature of this immunomodulation, peaking at the third month of treatment, provides the rationale for the use of antiangiogenics as a preconditioning strategy to improve the efficacy of ICB.

Network-based identification of key master regulators associated with an immune-silent cancer phenotype

A cancer immune phenotype characterized by an active T-helper 1 (Th1)/cytotoxic response is associated with responsiveness to immunotherapy and favorable prognosis across different tumors. However, in some cancers, such an intratumoral immune activation does not confer protection from progression or relapse. Defining mechanisms associated with immune evasion is imperative to refine stratification algorithms, to guide treatment decisions and to identify candidates for immune-targeted therapy. Molecular alterations governing mechanisms for immune exclusion are still largely unknown. The availability of large genomic datasets offers an opportunity to ascertain key determinants of differential intratumoral immune response. We follow a network-based protocol to identify transcription regulators (TRs) associated with poor immunologic antitumor activity. We use a consensus of four different pipelines consisting of two state-of-the-art gene regulatory network inference techniques, regularized gradient boosting machines and ARACNE to determine TR regulons, and three separate enrichment techniques, including fast gene set enrichment analysis, gene set variation analysis and virtual inference of protein activity by enriched regulon analysis to identify the most important TRs affecting immunologic antitumor activity. These TRs, referred to as master regulators (MRs), are unique to immune-silent and immune-active tumors, respectively. We validated the MRs coherently associated with the immune-silent phenotype across cancers in The Cancer Genome Atlas and a series of additional datasets in the Prediction of Clinical Outcomes from Genomic Profiles repository. A downstream analysis of MRs specific to the immune-silent phenotype resulted in the identification of several enriched candidate pathways, including NOTCH1, TGF-$\beta $, Interleukin-1 and TNF-$\alpha $ signaling pathways. TGFB1I1 emerged as one of the main negative immune modulators preventing the favorable effects of a Th1/cytotoxic response.

Organizing gene literature retrieval, profiling, and visualization training workshops for early career researchers

Developing the skills needed to effectively search and extract information from biomedical literature is essential for early-career researchers. It is, for instance, on this basis that the novelty of experimental results, and therefore publishing opportunities, can be evaluated. Given the unprecedented volume of publications in the field of biomedical research, new systematic approaches need to be devised and adopted for the retrieval and curation of literature relevant to a specific theme. Here we describe a hands-on training curriculum aimed at retrieval, profiling, and visualization of literature associated with a given topic. This curriculum was implemented in a workshop in January 2021. We provide supporting material and step-by-step implementation guidelines with the ISG15 gene literature serving as an illustrative use case. Through participation in such a workshop, trainees can learn: 1) to build and troubleshoot PubMed queries in order to retrieve the literature associated with a gene of interest; 2) to identify key concepts relevant to given themes (such as cell types, diseases, and biological processes); 3) to measure the prevalence of these concepts in the gene literature; 4) to extract key information from relevant articles, and 5) to develop a background section or summary on the basis of this information. Finally, trainees can learn to consolidate the structured information captured through this process for presentation via an interactive web application.

BloodGen3Module: Blood transcriptional module repertoire analysis and visualization using R

Motivation

We previously described the construction and characterization of fixed reusable blood transcriptional module repertoires. More recently we released a third iteration (‘BloodGen3’ module repertoire) that comprises 382 functionally annotated modules and encompasses 14 168 transcripts. Custom bioinformatic tools are needed to support downstream analysis, visualization and interpretation relying on such fixed module repertoires.

Results

We have developed and describe here an R package, BloodGen3Module. The functions of our package permit group comparison analyses to be performed at the module-level, and to display the results as annotated fingerprint grid plots. A parallel workflow for computing module repertoire changes for individual samples rather than groups of samples is also available; these results are displayed as fingerprint heatmaps. An illustrative case is used to demonstrate the steps involved in generating blood transcriptome repertoire fingerprints of septic patients. Taken together, this resource could facilitate the analysis and interpretation of changes in blood transcript abundance observed across a wide range of pathological and physiological states.

Availability and implementation

The BloodGen3Module package and documentation are freely available from Github: https://github.com/Drinchai/BloodGen3Module.

Supplementary information

Supplementary data are available at Bioinformatics online.

Ancestry-associated transcriptomic profiles of breast cancer in patients of African, Arab, and European ancestry

Breast cancer largely dominates the global cancer burden statistics; however, there are striking disparities in mortality rates across countries. While socioeconomic factors contribute to population-based differences in mortality, they do not fully explain disparity among women of African ancestry (AA) and Arab ancestry (ArA) compared to women of European ancestry (EA). In this study, we sought to identify molecular differences that could provide insight into the biology of ancestry-associated disparities in clinical outcomes. We applied a unique approach that combines the use of curated survival data from The Cancer Genome Atlas (TCGA) Pan-Cancer clinical data resource, improved single-nucleotide polymorphism-based inferred ancestry assignment, and a novel breast cancer subtype classification to interrogate the TCGA and a local Arab breast cancer dataset. We observed an enrichment of BasalMyo tumors in AA patients (38 vs 16.5% in EA, p = 1.30E - 10), associated with a significant worse overall (hazard ratio (HR) = 2.39, p = 0.02) and disease-specific survival (HR = 2.57, p = 0.03). Gene set enrichment analysis of BasalMyo AA and EA samples revealed differences in the abundance of T-regulatory and T-helper type 2 cells, and enrichment of cancer-related pathways with prognostic implications (AA: PI3K-Akt-mTOR and ErbB signaling; EA: EGF, estrogen-dependent and DNA repair signaling). Strikingly, AMPK signaling was associated with opposing prognostic connotation (AA: 10-year HR = 2.79, EA: 10-year HR = 0.34). Analysis of ArA patients suggests enrichment of BasalMyo tumors with a trend for differential enrichment of T-regulatory cells and AMPK signaling. Together, our findings suggest that the disparity in the clinical outcome of AA breast cancer patients is likely related to differences in cancer-related and microenvironmental features.

Germline genetic contribution to the immune landscape of cancer

Understanding the contribution of the host's genetic background to cancer immunity may lead to improved stratification for immunotherapy and to the identification of novel therapeutic targets. We investigated the effect of common and rare germline variants on 139 well-defined immune traits in ∼9000 cancer patients enrolled in TCGA. High heritability was observed for estimates of NK cell and T cell subset infiltration and for interferon signaling. Common variants of IFIH1, TMEM173 (STING1), and TMEM108 were associated with differential interferon signaling and variants mapping to RBL1 correlated with T cell subset abundance. Pathogenic or likely pathogenic variants in BRCA1 and in genes involved in telomere stabilization and Wnt-β-catenin also acted as immune modulators. Our findings provide evidence for the impact of germline genetics on the composition and functional orientation of the tumor immune microenvironment. The curated datasets, variants, and genes identified provide a resource toward further understanding of tumor-immune interactions.

A balance score between immune stimulatory and suppressive microenvironments identifies mediators of tumour immunity and predicts pan-cancer survival

BACKGROUND

The balance between immune-stimulatory and immune-suppressive mechanisms in the tumour microenvironment is associated with tumour rejection and can predict the efficacy of immune checkpoint-inhibition therapies.

METHODS

We consider the observed differences between the transcriptional programmes associated with cancer types where the levels of immune infiltration predict a favourable prognosis versus those in which the immune infiltration predicts an unfavourable prognosis and defined a score named Mediators of Immune Response Against Cancer in soLid microEnvironments (MIRACLE). MIRACLE deconvolves T cell infiltration, from inhibitory mechanisms, such as TGFβ, EMT and PI3Kγ signatures.

RESULTS

Our score outperforms current state-of-the-art immune signatures as a predictive marker of survival in TCGA (n = 9305, HR: 0.043, p value: 6.7 × 10^-36). In a validation cohort (n = 7623), MIRACLE predicts better survival compared to other immune metrics (HR: 0.1985, p value: 2.73 × 10^-38). MIRACLE also predicts response to checkpoint-inhibitor therapies (n = 333). The tumour-intrinsic factors inversely associated with the reported score such as EGFR, PRKAR1A and MAP3K1 are frequently associated with immune-suppressive phenotypes.

CONCLUSIONS

The association of cancer outcome with the level of infiltrating immune cells is mediated by the balance of activatory and suppressive factors. MIRACLE accounts for this balance and predicts favourable cancer outcomes.

Expression of NK cell receptor ligands in primary colorectal cancer tissue in relation to the phenotype of circulating NK- and NKT cells, and clinical outcome

INTRODUCTION

Natural killer (NK) cells and natural killer T (NKT) cells are implicated in the development and progression of colorectal cancer (CRC). Tumor cells express NK cell receptor ligands that modulate their function. This study aimed to investigate the expression of such ligands in CRC in relation to the phenotype of circulating NK- and NKT cells, and clinical outcome.

METHODS

Primary tumor tissues were analyzed for protein expression of NK cell ligands using immunohistochemistry with automated image analysis in a cohort of 78 CRC patients. For 24 of the 78 patients, RNA expression of NK cell ligands was analyzed in primary tumor tissue using RNA sequencing. Receptor expression on circulating NK- and NKT cells was previously measured by us in 71 of the 78 patients using flow cytometry.

RESULTS

High Proliferating Cell Nuclear Antigen (PCNA) protein expression in the primary tumor associated with shorter disease-free survival (DFS) of CRC patients (P = 0.026). A trend was observed towards shorter DFS in CRC patients with above-median galectin-3 protein expression in the primary tumor (P = 0.055). High protein expression of galectin-3, CD1d, and human leukocyte antigen (HLA) class I, and high RNA expression of UL16-binding protein (ULBP)-1, -2, and -5, and HLA-E in the tumor tissue correlated with low expression of the corresponding receptors on circulating NK- or NKT cells (P < 0.05).

CONCLUSIONS

Galectin-3 and PCNA expression in the primary tumor may be prognostic biomarkers in CRC patients. Furthermore, our results suggest that NK cell receptor ligands expressed by tumor cells may modulate the phenotype of circulating NK- and NKT cells, and facilitate immune escape of metastasizing cells.

Annexin A3 in sepsis: novel perspectives from an exploration of public transcriptome data

According to publicly available transcriptome datasets, the abundance of Annexin A3 (ANXA3) is robustly increased during the course of sepsis; however, no studies have examined the biological significance or clinical relevance of ANXA3 in this pathology. Here we explored this interpretation gap and identified possible directions for future research. Based on reference transcriptome datasets, we found that ANXA3 expression is restricted to neutrophils, is upregulated in vitro after exposure to plasma obtained from septic patients, and is associated with adverse clinical outcomes. Secondly, an increase in ANXA3 transcript abundance was also observed in vivo, in the blood of septic patients in multiple independent studies. ANXA3 is known to mediate calcium-dependent granules-phagosome fusion in support of microbicidal activity in neutrophils. More recent work has also shown that ANXA3 enhances proliferation and survival of tumour cells via a Caspase-3-dependent mechanism. And this same molecule is also known to play a critical role in regulation of apoptotic events in neutrophils. Thus, we posit that during sepsis ANXA3 might either play a beneficial role, by facilitating microbial clearance and resolution of the infection; or a detrimental role, by prolonging neutrophil survival, which is known to contribute to sepsis-mediated organ damage.

HLA-G: A New Immune Checkpoint in Cancer?

Human leukocyte antigen G (HLA-G), known as a central protein in providing immune tolerance to the fetus in pregnant women, is also studied for a possible role in tumor development. Many studies have claimed HLA-G as a new immune checkpoint in cancer. Therefore, HLA-G and its receptors might be targets for immune checkpoint blockade in cancer immunotherapy. In order to substantiate that HLA-G is indeed an immune checkpoint in cancer, two important questions need to be answered: (1) To what extent is HLA-G expressed in the tumor by cancer cells? and (2) What is the function of HLA-G in cancer immune evasion? In this review, we discuss these questions. We agree that HLA-G is a potentially new immune checkpoint in cancer, but additional evidence is required to show the extent of intra-tumor and inter-tumor expression. These studies should focus on tumor expression patterns of the seven different HLA-G isoforms and of the receptors for HLA-G. Furthermore, specific roles for the different HLA-G isoforms should be established.

Oncogenic states dictate the prognostic and predictive connotations of intratumoral immune response

BACKGROUND

An immune active cancer phenotype typified by a T helper 1 (Th-1) immune response has been associated with increased responsiveness to immunotherapy and favorable prognosis in some but not all cancer types. The reason of this differential prognostic connotation remains unknown.

METHODS

To explore the contextual prognostic value of cancer immune phenotypes, we applied a multimodal pan-cancer analysis among 31 different histologies (9282 patients), encompassing immune and oncogenic transcriptomic analysis, mutational and neoantigen load and copy number variations.

RESULTS

We demonstrated that the favorable prognostic connotation conferred by the presence of a Th-1 immune response was abolished in tumors displaying specific tumor-cell intrinsic attributes such as high transforming growth factor-beta (TGF-β) signaling and low proliferation capacity. This observation was independent of mutation rate. We validated this observation in the context of immune checkpoint inhibition. WNT-β catenin, barrier molecules, Notch, hedgehog, mismatch repair, telomerase activity and AMPK signaling were the pathways most coherently associated with an immune silent phenotype together with mutations of driver genes including IDH1/2, FOXA2, HDAC3, PSIP1, MAP3K1, KRAS, NRAS, EGFR, FGFR3, WNT5A and IRF7.

CONCLUSIONS

This is the first systematic study demonstrating that the prognostic and predictive role of a bona fide favorable intratumoral immune response is dependent on the disposition of specific oncogenic pathways. This information could be used to refine stratification algorithms and prioritize hierarchically relevant targets for combination therapies.

Genetic Variation in CCL5 Signaling Genes and Triple Negative Breast Cancer: Susceptibility and Prognosis Implications

Triple-negative breast cancer (TNBC) accounts for ~15–20% of breast cancer (BC) and has a higher rate of early relapse and mortality compared to other subtypes. The Chemokine (C-C motif) ligand 5 (CCL5) and its signaling pathway have been linked to TNBC. We aimed to investigate the susceptibility and prognostic implications of genetic variation in CCL5 signaling genes in TNBC in the present study. We characterized variants in CCL5 and that of six other CCL5 signaling genes (CCND1, ZMIZ1, CASP8, NOTCH2, MAP3K21, and HS6ST3) among 1,082 unrelated Tunisian subjects (544 BC patients, including 196 TNBC, and 538 healthy controls), assessed the association of the variants with BC-specific overall survival (OVS) and progression-free survival (PFS), and correlated CCL5 mRNA and serum levels with CCL5 genotypes. We found a highly significant association between the CCND1 rs614367-TT genotype (OR = 5.14; P = 0.004) and TNBC risk, and identified a significant association between the rs614367-T allele and decreased PFS in TNBC. A decreased risk of lymph node metastasis was associated with the MAP3K21 rs1294255-C allele, particularly in rs1294255-GC (OR = 0.47; P = 0.001). CCL5 variants (rs2107538 and rs2280789) were linked to CCL5 serum and mRNA levels. In the TCGA TNBC/Basal-like cohort the MAP3K21 rs1294255-G allele was associated with a decreased OVS. High expression of CCL5 in breast tumors was significantly associated with an increased OVS in all BC patients, but particularly in TNBC/Basal-like patients. In conclusion, genetic variation in CCL5 signaling genes may predict not only TNBC risk but also disease aggressiveness.

Long-Chain Acyl-CoA Synthetase 1 Role in Sepsis and Immunity: Perspectives From a Parallel Review of Public Transcriptome Datasets and of the Literature

A potential role for the long-chain acyl-CoA synthetase family member 1 (ACSL1) in the immunobiology of sepsis was explored during a hands-on training workshop. Participants first assessed the robustness of the potential gap in biomedical knowledge identified via an initial screen of public transcriptome data and of the literature associated with ACSL1. Increase in ACSL1 transcript abundance during sepsis was confirmed in several independent datasets. Querying the ACSL1 literature also confirmed the absence of reports associating ACSL1 with sepsis. Inferences drawn from both the literature (via indirect associations) and public transcriptome data (via correlation) point to the likely participation of ACSL1 and ACSL4, another family member, in inflammasome activation in neutrophils during sepsis. Furthermore, available clinical data indicate that levels of ACSL1 and ACSL4 induction was significantly higher in fatal cases of sepsis. This denotes potential translational relevance and is consistent with involvement in pathways driving potentially deleterious systemic inflammation. Finally, while ACSL1 expression was induced in blood in vitro by a wide range of pathogen-derived factors as well as TNF, induction of ACSL4 appeared restricted to flagellated bacteria and pathogen-derived TLR5 agonists and IFNG. Taken together, this joint review of public literature and omics data records points to two members of the acyl-CoA synthetase family potentially playing a role in inflammasome activation in neutrophils. Translational relevance of these observations in the context of sepsis and other inflammatory conditions remain to be investigated.

NY-ESO-1 Based Immunotherapy of Cancer: Current Perspectives

NY-ESO-1 or New York esophageal squamous cell carcinoma 1 is a well-known cancer-testis antigen (CTAs) with re-expression in numerous cancer types. Its ability to elicit spontaneous humoral and cellular immune responses, together with its restricted expression pattern, have rendered it a good candidate target for cancer immunotherapy. In this review, we provide background information on NY-ESO-1 expression and function in normal and cancerous tissues. Furthermore, NY-ESO-1-specific immune responses have been observed in various cancer types; however, their utility as biomarkers are not well determined. Finally, we describe the immune-based therapeutic options targeting NY-ESO-1 that are currently in clinical trial. We will highlight the recent advancements made in NY-ESO-1 cancer vaccines, adoptive T cell therapy, and combinatorial treatment with checkpoint inhibitors and will discuss the current trends for future NY-ESO-1 based immunotherapy. Cancer treatment has been revolutionized over the last few decades with immunotherapy emerging at the forefront. Immune-based interventions have shown promising results, providing a new treatment avenue for durable clinical responses in various cancer types. The majority of successful immunotherapy studies have been reported in liquid cancers, whereas these approaches have met many challenges in solid cancers. Effective immunotherapy in solid cancers is hampered by the complex, dynamic tumor microenvironment that modulates the extent and phenotype of the antitumor immune response. Furthermore, many solid tumor-associated antigens are not private but can be found in normal somatic tissues, resulting in minor to detrimental off-target toxicities. Therefore, there is an ongoing effort to identify tumor-specific antigens to target using various immune-based modalities. CTAs are considered good candidate targets for immunotherapy as they are characterized by a restricted expression in normal somatic tissues concomitant with a re-expression in solid epithelial cancers. Moreover, several CTAs have been found to induce a spontaneous immune response, NY-ESO-1 being the most immunogenic among the family members. Hence, this review will focus on NY-ESO-1 and discuss the past and current NY-ESO-1 targeted immunotherapeutic strategies.

The MAPK hypothesis: immune-regulatory effects of MAPK-pathway genetic dysregulations and implications for breast cancer immunotherapy

With the advent of checkpoint inhibition, immunotherapy has revolutionized the clinical management of several cancers, but has demonstrated limited efficacy in mammary carcinoma. Transcriptomic profiling of cancer samples defined distinct immunophenotypic categories characterized by different prognostic and predictive connotations. In breast cancer, genomic alterations leading to the dysregulation of mitogen-activated protein kinase (MAPK) pathways have been linked to an immune-silent phenotype associated with poor outcome and treatment resistance. These aberrations include mutations of MAP3K1 and MAP2K4, amplification of KRAS, BRAF, and RAF1, and truncations of NF1. Anticancer therapies targeting MAPK signaling by BRAF and MEK inhibitors have demonstrated clear immunologic effects. These off-target properties could be exploited to convert the immune-silent tumor phenotype into an immune-active one. Preclinical evidence supports that MAPK-pathway inhibition can dramatically increase the efficacy of immunotherapy. In this review, we provide a detailed overview of the immunomodulatory impact of MAPK-pathway blockade through BRAF and MEK inhibitions. While BRAF inhibition might be relevant in melanoma only, MEK inhibition is potentially applicable to a wide range of tumors. Context-dependent similarities and differences of MAPK modulation will be dissected, in light of the complexity of the MAPK pathways. Therapeutic strategies combining the favorable effects of MAPK-oriented interventions on the tumor microenvironment while maintaining T-cell function will be presented. Finally, we will discuss recent studies highlighting the rationale for the implementation of MAPK-interference approaches in combination with checkpoint inhibitors and immune agonists in breast cancer.

A collection of annotated and harmonized human breast cancer transcriptome datasets, including immunologic classification

The increased application of high-throughput approaches in translational research has expanded the number of publicly available data repositories. Gathering additional valuable information contained in the datasets represents a crucial opportunity in the biomedical field. To facilitate and stimulate utilization of these datasets, we have recently developed an interactive data browsing and visualization web application, the Gene Expression Browser (GXB). In this note, we describe a curated compendium of 13 public datasets on human breast cancer, representing a total of 2142 transcriptome profiles. We classified the samples according to different immune based classification systems and integrated this information into the datasets. Annotated and harmonized datasets were uploaded to GXB. Study samples were categorized in different groups based on their immunologic tumor response profiles, intrinsic molecular subtypes and multiple clinical parameters. Ranked gene lists were generated based on relevant group comparisons. In this data note, we demonstrate the utility of GXB to evaluate the expression of a gene of interest, find differential gene expression between groups and investigate potential associations between variables with a specific focus on immunologic classification in breast cancer. This interactive resource is publicly available online at: http://breastcancer.gxbsidra.org/dm3/geneBrowser/list.