Stem-cell states converge in multistage cutaneous squamous cell carcinoma development.

Stem cells play a critical role in cancer development by contributing to cell heterogeneity, lineage plasticity, and drug resistance. We created gene expression networks from hundreds of mouse tissue samples (both normal and tumor) and integrated these with lineage tracing and single-cell RNA-seq, to identify convergence of cell states in premalignant tumor cells expressing markers of lineage plasticity and drug resistance. Two of these cell states representing multilineage plasticity or proliferation were inversely correlated, suggesting a mutually exclusive relationship. Treatment of carcinomas in vivo with chemotherapy repressed the proliferative state and activated multilineage plasticity whereas inhibition of differentiation repressed plasticity and potentiated responses to cell cycle inhibitors. Manipulation of this cell state transition point may provide a source of potential combinatorial targets for cancer therapy.

The PTPN2/PTPN1 inhibitor ABBV-CLS-484 unleashes potent anti-tumour immunity.

Immune checkpoint blockade is effective for some patients with cancer, but most are refractory to current immunotherapies and new approaches are needed to overcome resistance1,2. The protein tyrosine phosphatases PTPN2 and PTPN1 are central regulators of inflammation, and their genetic deletion in either tumour cells or immune cells promotes anti-tumour immunity3-6. However, phosphatases are challenging drug targets; in particular, the active site has been considered undruggable. Here we present the discovery and characterization of ABBV-CLS-484 (AC484), a first-in-class, orally bioavailable, potent PTPN2 and PTPN1 active-site inhibitor. AC484 treatment in vitro amplifies the response to interferon and promotes the activation and function of several immune cell subsets. In mouse models of cancer resistant to PD-1 blockade, AC484 monotherapy generates potent anti-tumour immunity. We show that AC484 inflames the tumour microenvironment and promotes natural killer cell and CD8+ T cell function by enhancing JAK-STAT signalling and reducing T cell dysfunction. Inhibitors of PTPN2 and PTPN1 offer a promising new strategy for cancer immunotherapy and are currently being evaluated in patients with advanced solid tumours (ClinicalTrials.gov identifier NCT04777994 ). More broadly, our study shows that small-molecule inhibitors of key intracellular immune regulators can achieve efficacy comparable to or exceeding that of antibody-based immune checkpoint blockade in preclinical models. Finally, to our knowledge, AC484 represents the first active-site phosphatase inhibitor to enter clinical evaluation for cancer immunotherapy and may pave the way for additional therapeutics that target this important class of enzymes.

Gene networks reveal stem-cell state convergence during preneoplasia and progression to malignancy in multistage skin carcinogenesis.

Adult mammalian stem cells play critical roles in normal tissue homeostasis, as well as in tumor development, by contributing to cell heterogeneity, plasticity, and development of drug resistance. The relationship between different types of normal and cancer stem cells is highly controversial and poorly understood. Here, we carried out gene expression network analysis of normal and tumor samples from genetically heterogeneous mice to create network metagenes for visualization of stem-cell networks, rather than individual stem-cell markers, at the single-cell level during multistage carcinogenesis. We combined this approach with lineage tracing and single-cell RNASeq of stem cells and their progeny, identifying a previously unrecognized hierarchy in which Lgr6+ stem cells from tumors generate progeny that express a range of other stem-cell markers including Sox2, Pitx1, Foxa1, Klf5, and Cd44. Our data identify a convergence of multiple stem-cell and tumor-suppressor pathways in benign tumor cells expressing markers of lineage plasticity and oxidative stress. This same single-cell population expresses network metagenes corresponding to markers of cancer drug resistance in human tumors of the skin, lung and prostate. Treatment of mouse squamous carcinomas in vivo with the chemotherapeutic cis-platin resulted in elevated expression of the genes that mark this cell population. Our data have allowed us to create a simplified model of multistage carcinogenesis that identifies distinct stem-cell states at different stages of tumor progression, thereby identifying networks involved in lineage plasticity, drug resistance, and immune surveillance, providing a rich source of potential targets for cancer therapy.

iBRIDGE: A Data Integration Method to Identify Inflamed Tumors from Single-cell RNA-Seq Data and Differentiate Cell Type-Specific Markers of Immune-Cell Infiltration.

The development of immune checkpoint-based immunotherapies has been a major advancement in the treatment of cancer, with a subset of patients exhibiting durable clinical responses. A predictive biomarker for immunotherapy response is the preexisting T-cell infiltration in the tumor immune microenvironment (TIME). Bulk transcriptomics-based approaches can quantify the degree of T-cell infiltration using deconvolution methods and identify additional markers of inflamed/cold cancers at the bulk level. However, bulk techniques are unable to identify biomarkers of individual cell types. Although single-cell RNA sequencing (scRNA-seq) assays are now being used to profile the TIME, to our knowledge there is no method of identifying patients with a T-cell inflamed TIME from scRNA-seq data. Here, we describe a method, iBRIDGE, which integrates reference bulk RNA-seq data with the malignant subset of scRNA-seq datasets to identify patients with a T-cell inflamed TIME. Using two datasets with matched bulk data, we show iBRIDGE results correlated highly with bulk assessments (0.85 and 0.9 correlation coefficients). Using iBRIDGE, we identified markers of inflamed phenotypes in malignant cells, myeloid cells, and fibroblasts, establishing type I and type II interferon pathways as dominant signals, especially in malignant and myeloid cells, and finding the TGFß-driven mesenchymal phenotype not only in fibroblasts but also in malignant cells. Besides relative classification, per-patient average iBRIDGE scores and independent RNAScope quantifications were used for threshold-based absolute classification. Moreover, iBRIDGE can be applied to in vitro grown cancer cell lines and can identify the cell lines that are adapted from inflamed/cold patient tumors.

Chemoradiation therapy alters the PD-L1 score in locoregional recurrent squamous cell carcinomas of the head and neck.

PD-L1 testing guides therapeutic decision-making for head and neck squamous cell carcinoma (HNSCC). We sought to understand whether chemoradiation therapy (CRT) influences the PD-L1 combined positive score (CPS) and other biomarkers of response to immunotherapy. PD-L1 expression was assessed using immunohistochemistry, and bulk RNA sequencing was performed on 146 HNSCC patients (65 primary sites, 50 paired local recurrences, and 31 paired regional recurrences). PD-L1 was scored using the CPS of ≥1, ≥20, and ≥50. Overall, 98 %, 54 %, and 17 % of HNSCCs had a CPS ≥1, ≥20, and ≥50, respectively. When using a cut-off of ≥1, CRT did not significantly change CPS at the locoregional recurrent site. However, there were significant changes when using CPS ≥20 or ≥50. The CPS changed for 32 % of patients when using a CPS ≥20 (p < 0.001). When using a CPS ≥50, there was a 20-23 % (p = 0.0058-0.00067) discordance rate at the site of locoregional recurrence. Oral cavity cancers had a significantly higher discordant rate than other primary sites for CPS ≥50, 44 % (8/18, p = 0.0058) and 58 % (7/12, p = 0.00067) discordance at the site of local and regional recurrence, respectively. When evaluating the 18 gene IFN-É£ signature predictive of response to anti-PD-1 blockade, there was a statistically significant increase in the IFN-É£ signature in recurrent larynx cancer (p = 0.02). Our study demonstrates that when using a higher cut-off of CPS ≥20 and ≥50, a repeat biopsy may be warranted after CRT for local and regional recurrent HNSCCs.

An anti-PD-1-GITR-L bispecific agonist induces GITR clustering-mediated T cell activation for cancer immunotherapy.

Costimulatory receptors such as glucocorticoid-induced tumor necrosis factor receptor-related protein (GITR) play key roles in regulating the effector functions of T cells. In human clinical trials, however, GITR agonist antibodies have shown limited therapeutic effect, which may be due to suboptimal receptor clustering-mediated signaling. To overcome this potential limitation, a rational protein engineering approach is needed to optimize GITR agonist-based immunotherapies. Here we show a bispecific molecule consisting of an anti-PD-1 antibody fused with a multimeric GITR ligand (GITR-L) that induces PD-1-dependent and FcγR-independent GITR clustering, resulting in enhanced activation, proliferation and memory differentiation of primed antigen-specific GITR+PD-1+ T cells. The anti-PD-1-GITR-L bispecific is a PD-1-directed GITR-L construct that demonstrated dose-dependent, immunologically driven tumor growth inhibition in syngeneic, genetically engineered and xenograft humanized mouse tumor models, with a dose-dependent correlation between target saturation and Ki67 and TIGIT upregulation on memory T cells. Anti-PD-1-GITR-L thus represents a bispecific approach to directing GITR agonism for cancer immunotherapy.

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.

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.

Mutational signatures in tumours induced by high and low energy radiation in Trp53 deficient mice.

Ionising radiation (IR) is a recognised carcinogen responsible for cancer development in patients previously treated using radiotherapy, and in individuals exposed as a result of accidents at nuclear energy plants. However, the mutational signatures induced by distinct types and doses of radiation are unknown. Here, we analyse the genetic architecture of mammary tumours, lymphomas and sarcomas induced by high (56Fe-ions) or low (gamma) energy radiation in mice carrying Trp53 loss of function alleles. In mammary tumours, high-energy radiation is associated with induction of focal structural variants, leading to genomic instability and Met amplification. Gamma-radiation is linked to large-scale structural variants and a point mutation signature associated with oxidative stress. The genomic architecture of carcinomas, sarcomas and lymphomas arising in the same animals are significantly different. Our study illustrates the complex interactions between radiation quality, germline Trp53 deficiency and tissue/cell of origin in shaping the genomic landscape of IR-induced tumours.

Immune oncology, immune responsiveness and the theory of everything.

Anti-cancer immunotherapy is encountering its own checkpoint. Responses are dramatic and long lasting but occur in a subset of tumors and are largely dependent upon the pre-existing immune contexture of individual cancers. Available data suggest that three landscapes best define the cancer microenvironment: immune-active, immune-deserted and immune-excluded. This trichotomy is observable across most solid tumors (although the frequency of each landscape varies depending on tumor tissue of origin) and is associated with cancer prognosis and response to checkpoint inhibitor therapy (CIT). Various gene signatures (e.g. Immunological Constant of Rejection - ICR and Tumor Inflammation Signature - TIS) that delineate these landscapes have been described by different groups. In an effort to explain the mechanisms of cancer immune responsiveness or resistance to CIT, several models have been proposed that are loosely associated with the three landscapes. Here, we propose a strategy to integrate compelling data from various paradigms into a "Theory of Everything". Founded upon this unified theory, we also propose the creation of a task force led by the Society for Immunotherapy of Cancer (SITC) aimed at systematically addressing salient questions relevant to cancer immune responsiveness and immune evasion. This multidisciplinary effort will encompass aspects of genetics, tumor cell biology, and immunology that are pertinent to the understanding of this multifaceted problem.

Lgr6 is a stem cell marker in mouse skin squamous cell carcinoma.

The G-protein-coupled receptors LGR4, LGR5 and LGR6 are Wnt signaling mediators, but their functions in squamous cell carcinoma (SCC) are unclear. Using lineage tracing in Lgr5-EGFP-CreERT2/Rosa26-Tomato and Lgr6-EGFP-CreERT2/Rosa26-Tomato reporter mice, we demonstrate that Lgr6, but not Lgr5, acts as an epithelial stem cell marker in SCCs in vivo. We identify, by single-molecule in situ hybridization and cell sorting, rare cells positive for Lgr6 expression in immortalized keratinocytes and show that their frequency increases in advanced SCCs. Lgr6 expression is enriched in cells with stem cell characteristics, and Lgr6 downregulation in vivo causes increased epidermal proliferation with expanded lineage tracing from epidermal stem cells positive for Lgr6 expression. Surprisingly, mice with germline knockout of Lgr6 are predisposed to SCC development, through a mechanism that includes compensatory upregulation of Lgr5. These data provide a model for human patients with germline loss-of-function mutations in Wnt pathway genes, including RSPO1 or LGR4, who show increased susceptibility to squamous tumor development.

Panx3 links body mass index and tumorigenesis in a genetically heterogeneous mouse model of carcinogen-induced cancer.

BACKGROUND: Body mass index (BMI) has been implicated as a primary factor influencing cancer development. However, understanding the relationship between these two complex traits has been confounded by both environmental and genetic heterogeneity. METHODS: In order to gain insight into the genetic factors linking BMI and cancer, we performed chemical carcinogenesis on a genetically heterogeneous cohort of interspecific backcross mice ((Mus Spretus × FVB/N) F1 × FVB/N). Using this cohort, we performed quantitative trait loci (QTL) analysis to identify regions linked to BMI. We then performed an integrated analysis incorporating gene expression, sequence comparison between strains, and gene expression network analysis to identify candidate genes influencing both tumor development and BMI. RESULTS: Analysis of QTL linked to tumorigenesis and BMI identified several loci associated with both phenotypes. Exploring these loci in greater detail revealed a novel relationship between the Pannexin 3 gene (Panx3) and both BMI and tumorigenesis. Panx3 is positively associated with BMI and is strongly tied to a lipid metabolism gene expression network. Pre-treatment Panx3 gene expression levels in normal skin are associated with tumor susceptibility and inhibition of Panx function strongly influences inflammation. CONCLUSIONS: These studies have identified several genetic loci that influence both BMI and carcinogenesis and implicate Panx3 as a candidate gene that links these phenotypes through its effects on inflammation and lipid metabolism.

Gene Expression Architecture of Mouse Dorsal and Tail Skin Reveals Functional Differences in Inflammation and Cancer.

Inherited germline polymorphisms can cause gene expression levels in normal tissues to differ substantially between individuals. We present an analysis of the genetic architecture of normal adult skin from 470 genetically unique mice, demonstrating the effect of germline variants, skin tissue location, and perturbation by exogenous inflammation or tumorigenesis on gene signaling pathways. Gene networks related to specific cell types and signaling pathways, including sonic hedgehog (Shh), Wnt, Lgr family stem cell markers, and keratins, differed at these tissue sites, suggesting mechanisms for the differential susceptibility of dorsal and tail skin to development of skin diseases and tumorigenesis. The Pten tumor suppressor gene network is rewired in premalignant tumors compared to normal tissue, but this response to perturbation is lost during malignant progression. We present a software package for expression quantitative trait loci (eQTL) network analysis and demonstrate how network analysis of whole tissues provides insights into interactions between cell compartments and signaling molecules.

Evolution of metastasis revealed by mutational landscapes of chemically induced skin cancers.

Human tumors show a high level of genetic heterogeneity, but the processes that influence the timing and route of metastatic dissemination of the subclones are unknown. Here we have used whole-exome sequencing of 103 matched benign, malignant and metastatic skin tumors from genetically heterogeneous mice to demonstrate that most metastases disseminate synchronously from the primary tumor, supporting parallel rather than linear evolution as the predominant model of metastasis. Shared mutations between primary carcinomas and their matched metastases have the distinct A-to-T signature of the initiating carcinogen dimethylbenzanthracene, but non-shared mutations are primarily G-to-T, a signature associated with oxidative stress. The existence of carcinomas that either did or did not metastasize in the same host animal suggests that there are tumor-intrinsic factors that influence metastatic seeding. We also demonstrate the importance of germline polymorphisms in determining allele-specific mutations, and we identify somatic genetic alterations that are specifically related to initiation of carcinogenesis by Hras or Kras mutations. Mouse tumors that mimic the genetic heterogeneity of human cancers can aid our understanding of the clonal evolution of metastasis and provide a realistic model for the testing of novel therapies.

The mutational landscapes of genetic and chemical models of Kras-driven lung cancer.

Next-generation sequencing of human tumours has refined our understanding of the mutational processes operative in cancer initiation and progression, yet major questions remain regarding the factors that induce driver mutations and the processes that shape mutation selection during tumorigenesis. Here we performed whole-exome sequencing on adenomas from three mouse models of non-small-cell lung cancer, which were induced either by exposure to carcinogens (methyl-nitrosourea (MNU) and urethane) or by genetic activation of Kras (Kras(LA2)). Although the MNU-induced tumours carried exactly the same initiating mutation in Kras as seen in the Kras(LA2) model (G12D), MNU tumours had an average of 192 non-synonymous, somatic single-nucleotide variants, compared with only six in tumours from the Kras(LA2) model. By contrast, the Kras(LA2) tumours exhibited a significantly higher level of aneuploidy and copy number alterations compared with the carcinogen-induced tumours, suggesting that carcinogen-induced and genetically engineered models lead to tumour development through different routes. The wild-type allele of Kras has been shown to act as a tumour suppressor in mouse models of non-small-cell lung cancer. We demonstrate that urethane-induced tumours from wild-type mice carry mostly (94%) Kras Q61R mutations, whereas those from Kras heterozygous animals carry mostly (92%) Kras Q61L mutations, indicating a major role for germline Kras status in mutation selection during initiation. The exome-wide mutation spectra in carcinogen-induced tumours overwhelmingly display signatures of the initiating carcinogen, while adenocarcinomas acquire additional C > T mutations at CpG sites. These data provide a basis for understanding results from human tumour genome sequencing, which has identified two broad categories of tumours based on the relative frequency of single-nucleotide variations and copy number alterations, and underline the importance of carcinogen models for understanding the complex mutation spectra seen in human cancers.

The Drosophila surface glia transcriptome: evolutionary conserved blood-brain barrier processes.

Central nervous system (CNS) function is dependent on the stringent regulation of metabolites, drugs, cells, and pathogens exposed to the CNS space. Cellular blood-brain barrier (BBB) structures are highly specific checkpoints governing entry and exit of all small molecules to and from the brain interstitial space, but the precise mechanisms that regulate the BBB are not well understood. In addition, the BBB has long been a challenging obstacle to the pharmacologic treatment of CNS diseases; thus model systems that can parse the functions of the BBB are highly desirable. In this study, we sought to define the transcriptome of the adult Drosophila melanogaster BBB by isolating the BBB surface glia with fluorescence activated cell sorting (FACS) and profiling their gene expression with microarrays. By comparing the transcriptome of these surface glia to that of all brain glia, brain neurons, and whole brains, we present a catalog of transcripts that are selectively enriched at the Drosophila BBB. We found that the fly surface glia show high expression of many ATP-binding cassette (ABC) and solute carrier (SLC) transporters, cell adhesion molecules, metabolic enzymes, signaling molecules, and components of xenobiotic metabolism pathways. Using gene sequence-based alignments, we compare the Drosophila and Murine BBB transcriptomes and discover many shared chemoprotective and small molecule control pathways, thus affirming the relevance of invertebrate models for studying evolutionary conserved BBB properties. The Drosophila BBB transcriptome is valuable to vertebrate and insect biologists alike as a resource for studying proteins underlying diffusion barrier development and maintenance, glial biology, and regulation of drug transport at tissue barriers.