Immune status is prognostic for poor survival in colorectal cancer patients and is associated with tumour hypoxia

Background Immunohistochemical quantification of the immune response is prognostic for colorectal cancer (CRC). Here, we evaluate the suitability of alternative immune classifiers on prognosis and assess whether they relate to biological features amenable to targeted therapy. Methods Overall survival by immune (CD3, CD4, CD8, CD20 and FOXP3) and immune-checkpoint (ICOS, IDO-1 and PD-L1) biomarkers in independent CRC cohorts was evaluated. Matched mutational and transcriptomic data were interrogated to identify associated biology. Results Determination of immune-cold tumours by combined low-density cell counts of CD3, CD4 and CD8 immunohistochemistry constituted the best prognosticator across stage II–IV CRC, particularly in patients with stage IV disease (HR 1.98 [95% CI: 1.47–2.67]). These immune-cold CRCs were associated with tumour hypoxia, confirmed using CAIX immunohistochemistry (P = 0.0009), which may mediate disease progression through common biology (KRAS mutations, CRIS-B subtype and SPP1 mRNA overexpression). Conclusions Given the significantly poorer survival of immune-cold CRC patients, these data illustrate that assessment of CD4-expressing cells complements low CD3 and CD8 immunohistochemical quantification in the tumour bulk, potentially facilitating immunophenotyping of patient biopsies to predict prognosis. In addition, we found immune-cold CRCs to associate with a difficult-to-treat, poor prognosis hypoxia signature, indicating that these patients may benefit from hypoxia-targeting clinical trials

β-defensin 1 expression in HCV infected liver/liver cancer: an important role in protecting HCV progression and liver cancer development

Abstract β-defensin family plays a role in host defense against viral infection, however its role in HCV infection is still unknown. In this study, we demonstrated that β-defensin 1 was significantly reduced in HCV-infected liver specimens. Treatment with interferon and ribavirin upregulated β-defensin-1, but not other β-defensin tested, with the extent and duration of upregulation associated with treatment response. We investigated β-defensin family expression in liver cancer in publicly available datasets and found that among all the β-defensins tested, only β-defensin 1 was significantly downregulated, suggesting β-defensin 1 plays a crucial role in liver cancer development. Further analysis identified E-cadherin as the top positive correlated gene, while hepatocyte growth factor-regulated tyrosine kinase substrate as the top negative correlated gene. Expression of two proteoglycans were also positively correlated with that of β-defensin 1. We have also identified small molecules as potential therapeutic agents to reverse β-defensin 1-associated gene signature. Furthermore, the downregulation of β-defensin 1 and E-cadherin, and upregulation of hepatocyte growth factor-regulated tyrosine kinase substrate, were further confirmed in liver cancer and adjacent normal tissue collected from in-house Chinese liver cancer patients. Together, our results suggest β-defensin 1 plays an important role in protecting HCV progression and liver cancer development

Rational drug repurposing using sscMap analysis in a HOX-TALE model of leukemia

Drug discovery and development are often hampered by lack of target identification and clinical tractability. Repurposing of approved drugs to life-threatening diseases such as leukemia is emerging as a promising alternative approach. Connectivity mapping systems link approved drugs with disease-related gene signatures. Relevant preclinical models provide essential tools for system validation and proof-of-concept studies. Herein we describe procedures aimed at generating disease-based gene signatures and applying them to established cross-referencing databases of potential candidate drugs. As a proof of principle, we present the identification of Entinostat as a candidate drug for the treatment of HOX TALE-related leukemia

Transcriptional Profiling of a Patient-Matched Cohort of Glioblastoma (IDH-Wildtype) for Therapeutic Target and Repurposing Drug Identification.

Glioblastoma (GBM) is the most prevalent and aggressive adult brain tumor. Despite multi-modal therapies, GBM recurs, and patients have poor survival (~14 months). Resistance to therapy may originate from a subpopulation of tumor cells identified as glioma-stem cells (GSC), and new treatments are urgently needed to target these. The biology underpinning GBM recurrence was investigated using whole transcriptome profiling of patient-matched initial and recurrent GBM (recGBM). Differential expression analysis identified 147 significant probes. In total, 24 genes were validated using expression data from four public cohorts and the literature. Functional analyses revealed that transcriptional changes to recGBM were dominated by angiogenesis and immune-related processes. The role of MHC class II proteins in antigen presentation and the differentiation, proliferation, and infiltration of immune cells was enriched. These results suggest recGBM would benefit from immunotherapies. The altered gene signature was further analyzed in a connectivity mapping analysis with QUADrATiC software to identify FDA-approved repurposing drugs. Top-ranking target compounds that may be effective against GSC and GBM recurrence were rosiglitazone, nizatidine, pantoprazole, and tolmetin. Our translational bioinformatics pipeline provides an approach to identify target compounds for repurposing that may add clinical benefit in addition to standard therapies against resistant cancers such as GBM

Evolutionary genetic algorithm identifies <i>IL2RB</i> as a potential predictive biomarker for immune-checkpoint therapy in colorectal cancer.

Identifying robust predictive biomarkers to stratify colorectal cancer (CRC) patients based on their response to immune-checkpoint therapy is an area of unmet clinical need. Our evolutionary algorithm Atlas Correlation Explorer (ACE) represents a novel approach for mining The Cancer Genome Atlas (TCGA) data for clinically relevant associations. We deployed ACE to identify candidate predictive biomarkers of response to immune-checkpoint therapy in CRC. We interrogated the colon adenocarcinoma (COAD) gene expression data across nine immune-checkpoints (PDL1, PDCD1, CTLA4, LAG3, TIM3, TIGIT, ICOS, IDO1 and BTLA). IL2RB was identified as the most common gene associated with immune-checkpoint genes in CRC. Using human/murine single-cell RNA-seq data, we demonstrated that IL2RB was expressed predominantly in a subset of T-cells associated with increased immune-checkpoint expression (P n = 115) revealed sensitive, specific staining of a subset of lymphocytes and a strong association with FOXP3+ lymphocytes (P IL2RB mRNA positively correlated with three previously-published gene signatures of response to immune-checkpoint therapy (P IL2RB to be extensively linked to immune-checkpoints in CRC; its expression should be investigated for clinical utility as a potential predictive biomarker for CRC patients receiving immune-checkpoint blockade

Automated tumour recognition and digital pathology scoring unravels new role for PD-L1 in predicting good outcome in ER-/HER2+ breast cancer

The role of PD-L1 as a prognostic and predictive biomarker is an area of great interest. However, there is a lack of consensus on how to deliver PD-L1 as a clinical biomarker. At the heart of this conundrum is the subjective scoring of PD-L1 IHC in most studies to date. Current standard scoring systems involve separation of epithelial and inflammatory cells and find clinical significance in different percentages of expression, e.g., above or below 1%. Clearly, an objective, reproducible and accurate approach to PD-L1 scoring would bring a degree of necessary consistency to this landscape. Using a systematic comparison of technologies and the application of QuPath, a digital pathology platform, we show that high PD-L1 expression is associated with improved clinical outcome in Triple Negative breast cancer in the context of standard of care (SoC) chemotherapy, consistent with previous findings. In addition, we demonstrate for the first time that high PD-L1 expression is also associated with better outcome in ER- disease as a whole including HER2+ breast cancer. We demonstrate the influence of antibody choice on quantification and clinical impact with the Ventana antibody (SP142) providing the most robust assay in our hands. Through sampling different regions of the tumour, we show that tumour rich regions display the greatest range of PD-L1 expression and this has the most clinical significance compared to stroma and lymphoid rich areas. Furthermore, we observe that both inflammatory and epithelial PD-L1 expression are associated with improved survival in the context of chemotherapy. Moreover, as seen with PD-L1 inhibitor studies, a low threshold of PD-L1 expression stratifies patient outcome. This emphasises the importance of using digital pathology and precise biomarker quantitation to achieve accurate and reproducible scores that can discriminate low PD-L1 expression

Connectivity Mapping for Candidate Therapeutics Identification Using Next Generation Sequencing RNA-Seq Data

The advent of next generation sequencing technologies (NGS) has expanded the area of genomic research, offering high coverage and increased sensitivity over older microarray platforms. Although the current cost of next generation sequencing is still exceeding that of microarray approaches, the rapid advances in NGS will likely make it the platform of choice for future research in differential gene expression. Connectivity mapping is a procedure for examining the connections among diseases, genes and drugs by differential gene expression initially based on microarray technology, with which a large collection of compound-induced reference gene expression profiles have been accumulated. In this work, we aim to test the feasibility of incorporating NGS RNA-Seq data into the current connectivity mapping framework by utilizing the microarray based reference profiles and the construction of a differentially expressed gene signature from a NGS dataset. This would allow for the establishment of connections between the NGS gene signature and those microarray reference profiles, alleviating the associated incurring cost of re-creating drug profiles with NGS technology. We examined the connectivity mapping approach on a publicly available NGS dataset with androgen stimulation of LNCaP cells in order to extract candidate compounds that could inhibit the proliferative phenotype of LNCaP cells and to elucidate their potential in a laboratory setting. In addition, we also analyzed an independent microarray dataset of similar experimental settings. We found a high level of concordance between the top compounds identified using the gene signatures from the two datasets. The nicotine derivative cotinine was returned as the top candidate among the overlapping compounds with potential to suppress this proliferative phenotype. Subsequent lab experiments validated this connectivity mapping hit, showing that cotinine inhibits cell proliferation in an androgen dependent manner. Thus the results in this study suggest a promising prospect of integrating NGS data with connectivity mapping

Combination of a proteomics approach and reengineering of meso scale network models for prediction of mode-of-action for tyrosine kinase inhibitors

In drug discovery, the characterisation of the precise modes of action (MoA) and of unwanted off-target effects of novel molecularly targeted compounds is of highest relevance. Recent approaches for identification of MoA have employed various techniques for modeling of well defined signaling pathways including structural information, changes in phenotypic behavior of cells and gene expression patterns after drug treatment. However, efficient approaches focusing on proteome wide data for the identification of MoA including interference with mutations are underrepresented. As mutations are key drivers of drug resistance in molecularly targeted tumor therapies, efficient analysis and modeling of downstream effects of mutations on drug MoA is a key to efficient development of improved targeted anti-cancer drugs. Here we present a combination of a global proteome analysis, reengineering of network models and integration of apoptosis data used to infer the mode-of-action of various tyrosine kinase inhibitors (TKIs) in chronic myeloid leukemia (CML) cell lines expressing wild type as well as TKI resistance conferring mutants of BCR-ABL. The inferred network models provide a tool to predict the main MoA of drugs as well as to grouping of drugs with known similar kinase inhibitory activity patterns in comparison to drugs with an additional MoA. We believe that our direct network reconstruction approach, demonstrated on proteomics data, can provide a complementary method to the established network reconstruction approaches for the preclinical modeling of the MoA of various types of targeted drugs in cancer treatment. Hence it may contribute to the more precise prediction of clinically relevant on- and off-target effects of TKIs

TAZ Expression as a Prognostic Indicator in Colorectal Cancer

Agency of Science, Technology, and Research (A*STAR), SingaporeThe Hippo pathway restricts the activity of transcriptional coactivators TAZ (WWTR1) and YAP. TAZ and YAP are reported to be overexpressed in various cancers, however, their prognostic significance in colorectal cancers remains unstudied. The expression levels of TAZ and YAP, and their downstream transcriptional targets, AXL and CTGF, were extracted from two independent colon cancer patient datasets available in the Gene Expression Omnibus database, totaling 522 patients. We found that mRNA expressions of both TAZ and YAP were positively correlated with those of AXL and CTGF (p<0.05). High level mRNA expression of TAZ, AXL or CTGF significantly correlated with shorter survival. Importantly, patients co-overexpressing all 3 genes had a significantly shorter survival time, and combinatorial expression of these 3 genes was an independent predictor for survival. The downstream target genes for TAZ-AXL-CTGF overexpression were identified by Java application MyStats. Interestingly, genes that are associated with colon cancer progression (ANTXR1, EFEMP2, SULF1, TAGLN, VCAN, ZEB1 and ZEB2) were upregulated in patients co-overexpressing TAZ-AXL-CTGF. This TAZ-AXL-CTGF gene expression signature (GES) was then applied to Connectivity Map to identify small molecules that could potentially be utilized to reverse this GES. Of the top 20 small molecules identified by connectivity map, amiloride (a potassium sparing diuretic,) and tretinoin (all-trans retinoic acid) have shown therapeutic promise in inhibition of colon cancer cell growth. Using MyStats, we found that low level expression of either ANO1 or SQLE were associated with a better prognosis in patients who co-overexpressed TAZ-AXL-CTGF, and that ANO1 was an independent predictor of survival together with TAZ-AXL-CTGF. Finally, we confirmed that TAZ regulates Axl, and plays an important role in clonogenicity and non-adherent growth in vitro and tumor formation in vivo. These data suggest that TAZ could be a therapeutic target for the treatment of colon cancer