Chang-Qi Zhu

Cancer development requires a permissive microenvironment that is shaped by interactions between tumor cells, stroma, and the surrounding matrix. As collagen receptors, the leukocyte-associated immunoglobulin-like receptor (LAIR) family allows the immune system to interact with the extracellular matrix. However, little is known about their role in regulating tumor immunity and cancer progression. Methods: Genetic analysis of resected human lung adenocarcinoma was correlated to clinical-pathological characteristics, gene ontologies, and single cell RNA sequencing (scRNASeq). LAIR2 production was determined in subsets of immune cells isolated from blood leukocytes and lung adenocarcinoma tumor. Functional assays were used to determine the role of LAIR2 in tumorigenesis. Results: LAIR2 expression was adversely prognostic in lung adenocarcinoma. LAIR2 was preferentially produced by activated CD4⁺ T cells and enhanced in vitro tumor invasion into collagen. scRNASeq analysis of tumor infiltrating T cells revealed that LAIR2 expression co-localized with FOXP3 expressing cells and shared a transcriptional signature with tumor-associated regulatory T (T_reg) cells. A CD4⁺ LAIR2⁺ T_reg gene signature was prognostically significant in the TCGA dataset (n = 439; hazard ratio (HR) = 1.37; 95% confidence interval (CI), 1.05–1.77, p = 0.018) and validated in NCI Director’s Challenge lung adenocarcinoma dataset (n = 488; HR = 1.54; 95% CI, 1.14–2.09, p = 0.0045). Conclusions: Our data support a role for LAIR2 in lung adenocarcinoma tumorigenesis and identify a CD4⁺ LAIR2⁺ T_reg gene signature in lung adenocarcinoma prognosis. LAIR2 provides a novel target for development of immunotherapies.

Pancreatic ductal adenocarcinoma (PDAC) remains one of the most lethal human cancers, with 5-year patient survival rates of KRAS are the predominant oncogenic drivers of PDAC but are accompanied by additional lower frequency genetic alterations. Our group previously identified the guanine nucleotide exchange factor ARHGEF10 in a genomic screen for genes with copy number alterations that may synergize with oncogenic KRAS to promote PDAC carcinogenesis. In the present study we show that ARHGEF10 possesses putative tumor suppressor function in PDAC. ARHGEF10 expression is reduced in over 70% of PDAC cell lines, and copy number loss is documented in more than 30% of PDAC patient-derived xenografts. Loss of ARHGEF10 expression enhanced subcutaneous tumor growth in mouse models, while its exogenous expression greatly impaired tumorigenesis. Loss of ARHGEF10 expression also increased in vitro proliferation, invasion, and motility of PDAC cell lines, and enhanced their metastatic spread in orthotopic mouse models. Treatment of ARHGEF10-depleted cells with the inhibitor dasatinib reduced levels of phospho Src kinase and attenuated motility and invasion in vitro. Together, our data indicate that ARHGEF10 may function as a tumor suppressor in PDAC.

The burden of somatic mutations and neoantigens has been associated with improved survival in cancer treated with immunotherapies, especially non-small cell lung cancer (NSCLC). However, there is uncertainty about their effect on outcome in early-stage untreated cases. We posited that the burden of mutations in a specific set of genes may also contribute to the prognosis of early NSCLC patients. From a small cohort of 36 NSCLC cases, we were able to identify somatic mutations and copy number alterations in 865 genes that contributed to patient overall survival. Simply, the number of altered genes (NAG) among these 865 genes was associated with longer disease-free survival (hazard ratio (HR) = 0.153, p = 1.48 × 10⁻⁴). The gene expression signature distinguishing patients with high/low NAG was also prognostic in three independent datasets. Patients with a high NAG could be further stratified based on the presence of immunogenic mutations, revealing a further subgroup of stage I NSCLC with even better prognosis (85% with >5 years survival), and associated with cytotoxic T-cell expression. Importantly, 95% of the highly-altered genes lacked direct relation to cancer, but were implicated in pathways regulating cell proliferation, motility and immune response.

Background: Adjuvant chemotherapy (ACT) provides modest benefit in resected non-small cell lung cancer (NSCLC) patients. Genome-wide studies have identified gene copy number aberrations (CNA), but their prognostic implication is unknown. Methods: DNA from 1,013 FFPE tumor samples from three pivotal multicenter randomized trials (ACT vs. control) in the LACE-Bio consortium (median follow-up: 5.2 years) was successfully extracted, profiled using a molecular inversion probe SNP assay, normalized relative to a pool of normal tissues and segmented. Minimally recurrent regions were identified. P values were adjusted to control the false discovery rate (Q values). Results: A total of 976 samples successfully profiled, 414 (42%) adenocarcinoma (ADC), 430 (44%) squamous cell carcinoma (SCC) and 132 (14%) other NSCLC; 710 (73%) males. We identified 431 recurrent regions, with on average 51 gains and 43 losses; 253 regions (59%) were ≤3 Mb. Most frequent gains (up to 48%) were on chr1, 3q, 5p, 6p, 8q, 22q; most frequent losses (up to 40%) on chr3p, 8p, 9p. CNA frequency of 195 regions was significantly different (Q≤0.05) between ADC and SCC. Fourteen regions (7p11–12, 9p21, 18q12, and 19p11–13) were associated with disease-free survival (DFS) (univariate P≤0.005, QCDKN2A/B [hazard ratio (HR) for 2-fold lower CN: 1.5 (95% CI: 1.2–1.9), PSTK11 [HR =2.4 (1.3–4.3), P=0.005, Q=0.15]. Chromosomal instability was associated with poorer DFS (HR =1.5, P=0.015), OS (HR =1.2, P=0.189) and lung-cancer specific survival (HR =1.7, P=0.003). Conclusions: These large-scale genome-wide analyses of gene CNA provide new candidate prognostic markers for stage I–III NSCLC.

Cancer cell repopulation during chemotherapy or radiotherapy is a major factor limiting the efficacy of treatment. Cancer stem cells (CSC) may play critical roles during this process. We aim to demonstrate the role of mesothelioma stem cells (MSC) in treatment failure and eventually to design specific target therapies against MSC to improve the efficacy of treatment in malignant mesothelioma. Murine mesothelioma AB12 and RN5 cells were used to compare tumorigenicity in mice. The expression of CSC-associated genes was evaluated by quantitative real-time PCR in both cell lines treated with chemo-radiation. Stemness properties of MSC-enriched RN5-EOS-Puro2 cells were characterized with flow cytometry and immunostaining. A MSC-specific gene profile was screened by microarray assay and confirmed thereafter. Gene Ontology analysis of the selected genes was performed by GOMiner. Tumor growth delay of murine mesothelioma AB12 cells was achieved after each cycle of cisplatin treatment, however, tumors grew back rapidly due to cancer cell repopulation between courses of chemotherapy. Strikingly, a 10-times lower number of irradiated cells in both cell lines led to a similar tumor incidence and growth rate as with untreated cells. The expression of CSC-associated genes such as CD24, CD133, CD90 and uPAR was dramatically up-regulated, while others did not change significantly after chemoradiation. Highly enriched MSC after selection with puromycin displayed an increasing GFP-positive population and showed typical properties of stemness. Comparatively, the proportion of MSC significantly increased after RN5-EOS parental cells were treated with either chemotherapy, γ-ray radiation, or a combination of the two, while MSC showed more resistance to the above treatments. A group of identified genes are most likely MSC-specific, and major pathways related to regulation of cell growth or apoptosis are involved. Upregulation of the gene transcripts Tnfsf18, Serpinb9b, Ly6a, and Nppb were confirmed. Putative MSC possess the property of stemness showing more resistance to chemoradiation, suggesting that MSC may play critical roles in cancer cell repopulation. Further identification of selected genes may be used to design novel target therapies against MSC, so as to eliminate cancer cell repopulation in mesothelioma.

// <![CDATA[ $('.header-date').hide();$('#titleAuthors').hide();$('#abstractHeader').hide(); // ]]> Tomas Tokar1, Chiara Pastrello1, Varune R. Ramnarine1, 2, Chang-Qi Zhu1, Kenneth J. Craddock1, Larrisa A. Pikor3, Emily A. Vucic3, Simon Vary1, 4, 5, Frances A. Shepherd1, Ming-Sound Tsao1, 6, 7, Wan L. Lam3 and Igor Jurisica1, 6, 8, 9 1Princess Margaret Cancer Centre, University Health Network, Toronto, Canada 2The Vancouver Prostate Centre, Vancouver General Hospital, Vancouver, Canada 3Department of Integrative Oncology, British Columbia Cancer Research Centre, Vancouver, Canada 4Mathematical Institute, University of Oxford, Oxford, United Kingdom 5Faculty of Mathematics, Physics and Informatics, Comenius University, Bratislava, Slovakia 6Department of Medical Biophysics, University of Toronto, Toronto, Canada 7Department of Laboratory Medicine and Pathobiology, University of Toronto, Toronto, Canada 8Department of Computer Science, University of Toronto, Toronto, Canada 9Institute of Neuroimmunology, Slovak Academy of Sciences, Bratislava, Slovakia Correspondence to: Igor Jurisica, email: juris@ai.utoronto.ca Keywords: lung adenocarcinoma; copy number aberrations; microRNA; gene regulatory network; prognostic signature Received: August 25, 2017 Accepted: January 02, 2018 Published: January 08, 2018 ABSTRACT In many cancers, significantly down- or upregulated genes are found within chromosomal regions with DNA copy number alteration opposite to the expression changes. Generally, this paradox has been overlooked as noise, but can potentially be a consequence of interference of epigenetic regulatory mechanisms, including microRNA-mediated control of mRNA levels. To explore potential associations between microRNAs and paradoxes in non-small-cell lung cancer (NSCLC) we curated and analyzed lung adenocarcinoma (LUAD) data, comprising gene expressions, copy number aberrations (CNAs) and microRNA expressions. We integrated data from 1,062 tumor samples and 241 normal lung samples, including newly-generated array comparative genomic hybridization (aCGH) data from 63 LUAD samples. We identified 85 “paradoxical” genes whose differential expression consistently contrasted with aberrations of their copy numbers. Paradoxical status of 70 out of 85 genes was validated on sample-wise basis using The Cancer Genome Atlas (TCGA) LUAD data. Of these, 41 genes are prognostic and form a clinically relevant signature, which we validated on three independent datasets. By meta-analysis of results from 9 LUAD microRNA expression studies we identified 24 consistently-deregulated microRNAs. Using TCGA-LUAD data we showed that deregulation of 19 of these microRNAs explains differential expression of the paradoxical genes. Our results show that deregulation of paradoxical genes is crucial in LUAD and their expression pattern is maintained epigenetically, defying gene copy number status.

Cancer development and biology is influenced by the host immune system. Emerging data indicate that the context of immune cell infiltrates may contribute to cancer prognosis. However, the types of infiltrating immune cells that are critical for cancer development remain controversial. In attempts to gain insights into the immune networks that regulate and/or predict tumor progression, gene expression analysis was conducted on microarray datasets of resected tumor samples from 128 early-stage non–small cell lung cancer (NSCLC) adenocarcinoma patients. By limiting analysis to immune-related genes, we identified a 9-gene signature using MAximizing R Square Algorithm that selected for the greatest separation between favorable and adverse prognostic patient subgroups. The prognostic value of this 9-gene signature was validated in 10 additional independently published microarray datasets of lung adenocarcinoma [n = 1,097; overall survival hazard ratio (HR), 2.05; 95% confidence interval, 1.64–2.56; P < 0.0001] and was found to be an independent prognostic indicator relative to tumor stage (overall survival HR, 2.09, 95% confidence interval, 1.65–2.66; P < 0.0001). Network analysis revealed that genes associated with Fcε complex (FCER1, MS4A2) formed the largest and most significant pathway of the signature. Using immunohistochemistry, we validated that MS4A2, the β subunit of the IgE receptor expressed on mast cells, is a favorable prognostic indicator and show that MS4A2 gene expression is an independent prognostic marker for early-stage lung cancer patient survival. Cancer Immunol Res; 5(9); 821–9. ©2017 AACR.

Non-small-cell lung carcinoma (NSCLC) accounts for 85% of malignant lung tumors and is the leading cause of cancer deaths. Our group previously identified Tripartite Motif 14 (TRIM14) as a component of a prognostic multigene expression signature for NSCLC. Little is known about the function of TRIM14 protein in normal or disease states. We investigated the functional and prognostic role of TRIM14 in NSCLC using in vitro and in vivo perturbation model systems. Firstly, a pooled RNAi screen identified TRIM14 to effect cell proliferation/survival in NSCLC cells. Secondly, silencing of TRIM14 expression significantly enhanced tumor growth in NSCLC xenograft mouse models, while exogenous TRIM14 expression attenuated tumorigenesis. In addition, differences in apoptotic activity between TRIM14-deficient and control tumors suggests that TRIM14 tumor suppressor activity may depend on cell death signaling pathways. TRIM14-deficient cell lines showed both resistance to hypoxia-induced cell death and attenuation of interferon response via STAT1 signaling. Consistent with these phenotypes, multivariate analyses on published mRNA expression datasets of over 600 primary NSCLCs demonstrated that low TRIM14 mRNA levels are significantly associated with poorer prognosis in early stage NSCLC patients. Our functional data therefore establish a novel tumor suppressive role for TRIM14 in NSCLC progression.

Availability of lung cancer models that closely mimic human tumors remains a significant gap in cancer research, as tumor cell lines and mouse models may not recapitulate the spectrum of lung cancer heterogeneity seen in patients. We aimed to establish a patient‐derived tumor xenograft (PDX) resource from surgically resected non‐small cell lung cancer (NSCLC). Fresh tumor tissue from surgical resection was implanted and grown in the subcutaneous pocket of non‐obese severe combined immune deficient (NOD SCID) gamma mice. Subsequent passages were in NOD SCID mice. A subset of matched patient and PDX tumors and non‐neoplastic lung tissues were profiled by whole exome sequencing, single nucleotide polymorphism (SNP) and methylation arrays, and phosphotyrosine (pY)‐proteome by mass spectrometry. The data were compared to published NSCLC datasets of NSCLC primary and cell lines. 127 stable PDXs were established from 441 lung carcinomas representing all major histological subtypes: 52 adenocarcinomas, 62 squamous cell carcinomas, one adeno‐squamous carcinoma, five sarcomatoid carcinomas, five large cell neuroendocrine carcinomas, and two small cell lung cancers. Somatic mutations, gene copy number and expression profiles, and pY‐proteome landscape of 36 PDXs showed greater similarity with patient tumors than with established cell lines. Novel somatic mutations on cancer associated genes were identified but only in PDXs, likely due to selective clonal growth in the PDXs that allows detection of these low allelic frequency mutations. The results provide the strongest evidence yet that PDXs established from lung cancers closely mimic the characteristics of patient primary tumors.