IL-27 maintains cytotoxic Ly6C+ γδ T cells that arise from immature precursors.

In mice, γδ-T lymphocytes that express the co-stimulatory molecule, CD27, are committed to the IFNγ-producing lineage during thymic development. In the periphery, these cells play a critical role in host defense and anti-tumor immunity. Unlike αß-T cells that rely on MHC-presented peptides to drive their terminal differentiation, it is unclear whether MHC-unrestricted γδ-T cells undergo further functional maturation after exiting the thymus. Here, we provide evidence of phenotypic and functional diversity within peripheral IFNγ-producing γδ T cells. We found that CD27+ Ly6C- cells convert into CD27+Ly6C+ cells, and these CD27+Ly6C+ cells control cancer progression in mice, while the CD27+Ly6C- cells cannot. The gene signatures of these two subsets were highly analogous to human immature and mature γδ-T cells, indicative of conservation across species. We show that IL-27 supports the cytotoxic phenotype and function of mouse CD27+Ly6C+ cells and human Vδ2+ cells, while IL-27 is dispensable for mouse CD27+Ly6C- cell and human Vδ1+ cell functions. These data reveal increased complexity within IFNγ-producing γδ-T cells, comprising immature and terminally differentiated subsets, that offer new insights into unconventional T-cell biology.

Pathway level subtyping identifies a slow-cycling biological phenotype associated with poor clinical outcomes in colorectal cancer.

Molecular stratification using gene-level transcriptional data has identified subtypes with distinctive genotypic and phenotypic traits, as exemplified by the consensus molecular subtypes (CMS) in colorectal cancer (CRC). Here, rather than gene-level data, we make use of gene ontology and biological activation state information for initial molecular class discovery. In doing so, we defined three pathway-derived subtypes (PDS) in CRC: PDS1 tumors, which are canonical/LGR5+ stem-rich, highly proliferative and display good prognosis; PDS2 tumors, which are regenerative/ANXA1+ stem-rich, with elevated stromal and immune tumor microenvironmental lineages; and PDS3 tumors, which represent a previously overlooked slow-cycling subset of tumors within CMS2 with reduced stem populations and increased differentiated lineages, particularly enterocytes and enteroendocrine cells, yet display the worst prognosis in locally advanced disease. These PDS3 phenotypic traits are evident across numerous bulk and single-cell datasets, and demark a series of subtle biological states that are currently under-represented in pre-clinical models and are not identified using existing subtyping classifiers.

PTEN deficiency exposes a requirement for an ARF GTPase module for integrin-dependent invasion in ovarian cancer.

Dysregulation of the PI3K/AKT pathway is a common occurrence in high-grade serous ovarian carcinoma (HGSOC), with the loss of the tumour suppressor PTEN in HGSOC being associated with poor prognosis. The cellular mechanisms of how PTEN loss contributes to HGSOC are largely unknown. We here utilise time-lapse imaging of HGSOC spheroids coupled to a machine learning approach to classify the phenotype of PTEN loss. PTEN deficiency induces PI(3,4,5)P3 -rich and -dependent membrane protrusions into the extracellular matrix (ECM), resulting in a collective invasion phenotype. We identify the small GTPase ARF6 as a crucial vulnerability of HGSOC cells upon PTEN loss. Through a functional proteomic CRISPR screen of ARF6 interactors, we identify the ARF GTPase-activating protein (GAP) AGAP1 and the ECM receptor ß1-integrin (ITGB1) as key ARF6 interactors in HGSOC regulating PTEN loss-associated invasion. ARF6 functions to promote invasion by controlling the recycling of internalised, active ß1-integrin to maintain invasive activity into the ECM. The expression of the CYTH2-ARF6-AGAP1 complex in HGSOC patients is inversely associated with outcome, allowing the identification of patient groups with improved versus poor outcome. ARF6 may represent a therapeutic vulnerability in PTEN-depleted HGSOC.

Spatial regulation of the glycocalyx component podocalyxin is a switch for prometastatic function.

The glycocalyx component and sialomucin podocalyxin (PODXL) is required for normal tissue development by promoting apical membranes to form between cells, triggering lumen formation. Elevated PODXL expression is also associated with metastasis and poor clinical outcome in multiple tumor types. How PODXL presents this duality in effect remains unknown. We identify an unexpected function of PODXL as a decoy receptor for galectin-3 (GAL3), whereby the PODXL-GAL3 interaction releases GAL3 repression of integrin-based invasion. Differential cortical targeting of PODXL, regulated by ubiquitination, is the molecular mechanism controlling alternate fates. Both PODXL high and low surface levels occur in parallel subpopulations within cancer cells. Orthotopic intraprostatic xenograft of PODXL-manipulated cells or those with different surface levels of PODXL define that this axis controls metastasis in vivo. Clinically, interplay between PODXL-GAL3 stratifies prostate cancer patients with poor outcome. Our studies define the molecular mechanisms and context in which PODXL promotes invasion and metastasis.

PD-1 and TIM-3 differentially regulate subsets of mouse IL-17A-producing γδ T cells.

IL-17A-producing γδ T cells in mice consist primarily of Vγ6+ tissue-resident cells and Vγ4+ circulating cells. How these γδ T cell subsets are regulated during homeostasis and cancer remains poorly understood. Using single-cell RNA sequencing and flow cytommetry, we show that lung Vγ4+ and Vγ6+ cells from tumor-free and tumor-bearing mice express contrasting cell surface molecules as well as distinct co-inhibitory molecules, which function to suppress their expansion. Vγ6+ cells express constitutively high levels of PD-1, whereas Vγ4+ cells upregulate TIM-3 in response to tumor-derived IL-1ß and IL-23. Inhibition of either PD-1 or TIM-3 in mammary tumor-bearing mice increased Vγ6+ and Vγ4+ cell numbers, respectively. We found that genetic deletion of γδ T cells elicits responsiveness to anti-PD-1 and anti-TIM-3 immunotherapy in a mammary tumor model that is refractory to T cell checkpoint inhibitors, indicating that IL-17A-producing γδ T cells instigate resistance to immunotherapy. Together, these data demonstrate how lung IL-17A-producing γδ T cell subsets are differentially controlled by PD-1 and TIM-3 in steady-state and cancer.

Traject3d allows label-free identification of distinct co-occurring phenotypes within 3D culture by live imaging.

Single cell profiling by genetic, proteomic and imaging methods has expanded the ability to identify programmes regulating distinct cell states. The 3-dimensional (3D) culture of cells or tissue fragments provides a system to study how such states contribute to multicellular morphogenesis. Whether cells plated into 3D cultures give rise to a singular phenotype or whether multiple biologically distinct phenotypes arise in parallel is largely unknown due to a lack of tools to detect such heterogeneity. Here we develop Traject3d (Trajectory identification in 3D), a method for identifying heterogeneous states in 3D culture and how these give rise to distinct phenotypes over time, from label-free multi-day time-lapse imaging. We use this to characterise the temporal landscape of morphological states of cancer cell lines, varying in metastatic potential and drug resistance, and use this information to identify drug combinations that inhibit such heterogeneity. Traject3d is therefore an important companion to other single-cell technologies by facilitating real-time identification via live imaging of how distinct states can lead to alternate phenotypes that occur in parallel in 3D culture.

Cancer-associated fibroblasts require proline synthesis by PYCR1 for the deposition of pro-tumorigenic extracellular matrix.

Elevated production of collagen-rich extracellular matrix is a hallmark of cancer-associated fibroblasts (CAFs) and a central driver of cancer aggressiveness. Here we find that proline, a highly abundant amino acid in collagen proteins, is newly synthesized from glutamine in CAFs to make tumour collagen in breast cancer xenografts. PYCR1 is a key enzyme for proline synthesis and highly expressed in the stroma of breast cancer patients and in CAFs. Reducing PYCR1 levels in CAFs is sufficient to reduce tumour collagen production, tumour growth and metastatic spread in vivo and cancer cell proliferation in vitro. Both collagen and glutamine-derived proline synthesis in CAFs are epigenetically upregulated by increased pyruvate dehydrogenase-derived acetyl-CoA levels. PYCR1 is a cancer cell vulnerability and potential target for therapy; therefore, our work provides evidence that targeting PYCR1 may have the additional benefit of halting the production of a pro-tumorigenic extracellular matrix. Our work unveils new roles for CAF metabolism to support pro-tumorigenic collagen production.

Serum Proteomics and Plasma Fibulin-3 in Differentiation of Mesothelioma From Asbestos-Exposed Controls and Patients With Other Pleural Diseases.

INTRODUCTION: Malignant pleural mesothelioma (MPM) is difficult to diagnose. An accurate blood biomarker could prompt specialist referral or be deployed in future screening. In earlier retrospective studies, SOMAscan proteomics (Somalogic, Boulder, CO) and fibulin-3 seemed highly accurate, but SOMAscan has not been validated prospectively and subsequent fibulin-3 data have been contradictory. METHODS: A multicenter prospective observational study was performed in 22 centers, generating a large intention-to-diagnose cohort. Blood sampling, processing, and diagnostic assessment were standardized, including a 1-year follow-up. Plasma fibulin-3 was measured using two enzyme-linked immunosorbent assays (CloudClone [used in previous studies] and BosterBio, Pleasanton, CA). Serum proteomics was measured using the SOMAscan assay. Diagnostic performance (sensitivity at 95% specificity, area under the curve [AUC]) was benchmarked against serum mesothelin (Mesomark, Fujirebio Diagnostics, Malvern, PA). Biomarkers were correlated against primary tumor volume, inflammatory markers, and asbestos exposure. RESULTS: A total of 638 patients with suspected pleural malignancy (SPM) and 110 asbestos-exposed controls (AECs) were recruited. SOMAscan reliably differentiated MPM from AECs (75% sensitivity, 88.2% specificity, validation cohort AUC 0.855) but was not useful in patients with differentiating non-MPM SPM. Fibulin-3 (by BosterBio after failed CloudClone validation) revealed 7.4% and 11.9% sensitivity at 95% specificity in MPM versus non-MPM SPM and AECs, respectively (associated AUCs 0.611 [0.557-0.664], p = 0.0015) and 0.516 [0.443-0.589], p = 0.671), both inferior to mesothelin. SOMAscan proteins correlated with inflammatory markers but not with asbestos exposure. Neither biomarker correlated with tumor volume. CONCLUSIONS: SOMAscan may prove useful as a future screening test for MPM in asbestos-exposed persons. Neither fibulin-3 nor SOMAscan should be used for diagnosis or pathway stratification.

A model of k-mer surprisal to quantify local sequence information content surrounding splice regions.

Molecular sequences carry information. Analysis of sequence conservation between homologous loci is a proven approach with which to explore the information content of molecular sequences. This is often done using multiple sequence alignments to support comparisons between homologous loci. These methods therefore rely on sufficient underlying sequence similarity with which to construct a representative alignment. Here we describe a method using a formal metric of information, surprisal, to analyse biological sub-sequences without alignment constraints. We applied our model to the genomes of five different species to reveal similar patterns across a panel of eukaryotes. As the surprisal of a sub-sequence is inversely proportional to its occurrence within the genome, the optimal size of the sub-sequences was selected for each species under consideration. With the model optimized, we found a strong correlation between surprisal and CG dinucleotide usage. The utility of our model was tested by examining the sequences of genes known to undergo splicing. We demonstrate that our model can identify biological features of interest such as known donor and acceptor sites. Analysis across all annotated coding exon junctions in Homo sapiens reveals the information content of coding exons to be greater than the surrounding intron regions, a consequence of increased suppression of the CG dinucleotide in intronic space. Sequences within coding regions proximal to exon junctions exhibited novel patterns within DNA and coding mRNA that are not a function of the encoded amino acid sequence. Our findings are consistent with the presence of secondary information encoding features such as DNA and RNA binding sites, multiplexed through the coding sequence and independent of the information required to define the corresponding amino-acid sequence. We conclude that surprisal provides a complementary methodology with which to locate regions of interest in the genome, particularly in situations that lack an appropriate multiple sequence alignment.

Brain microenvironment-driven resistance to immune and targeted therapies in acral melanoma.

BACKGROUND: Combination treatments targeting the MEK-ERK pathway and checkpoint inhibitors have improved overall survival in melanoma. Resistance to treatment especially in the brain remains challenging, and rare disease subtypes such as acral melanoma are not typically included in trials. Here we present analyses from longitudinal sampling of a patient with metastatic acral melanoma that became resistant to successive immune and targeted therapies. METHODS: We performed whole-exome sequencing and RNA sequencing on an acral melanoma that progressed on successive immune (nivolumab) and targeted (dabrafenib) therapy in the brain to identify resistance mechanisms. In addition, we performed growth inhibition assays, reverse phase protein arrays and immunoblotting on patient-derived cell lines using dabrafenib in the presence or absence of cerebrospinal fluid (CSF) in vitro. Patient-derived xenografts were also developed to analyse response to dabrafenib. RESULTS: Immune escape following checkpoint blockade was not due to loss of tumour cell recognition by the immune system or low neoantigen burden, but was associated with distinct changes in the microenvironment. Similarly, resistance to targeted therapy was not associated with acquired mutations but upregulation of the AKT/phospho-inositide 3-kinase pathway in the presence of CSF. CONCLUSION: Heterogeneous tumour interactions within the brain microenvironment enable progression on immune and targeted therapies and should be targeted in salvage treatments.

Publisher Correction: Pulmonary venous circulating tumor cell dissemination before tumor resection and disease relapse.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Profiling of Circulating Free DNA Using Targeted and Genome-wide Sequencing in Patients with SCLC.

INTRODUCTION: SCLC accounts for approximately 250,000 deaths worldwide each year. Acquisition of adequate tumor biopsy samples is challenging, and liquid biopsies present an alternative option for patient stratification and response monitoring. METHODS: We applied whole genome next-generation sequencing to circulating free DNA (cfDNA) from 39 patients with limited-stage (LS) SCLC and 30 patients with extensive-stage SCLC to establish genome-wide copy number aberrations and also performed targeted mutation analysis of 110 SCLC associated genes. Quantitative metrics were calculated for copy number aberrations, including percent genome amplified (PGA [the percentage of genomic regions amplified]), Z-score (a measure of standard deviation), and Moran's I (a measure of spatial autocorrelation). In addition CellSearch, an epitope-dependent enrichment platform, was used to enumerate circulating tumor cells (CTCs) from a parallel blood sample. RESULTS: Genome-wide and targeted cfDNA sequencing data identified tumor-related changes in 94% of patients with LS SCLC and 100% of patients with extensive-stage SCLC. Parallel analysis of CTCs based on at least 1 CTC/7.5 mL of blood increased tumor detection frequencies to 95% for LS SCLC. Both CTC counts and cfDNA readouts correlated with disease stage and overall survival. CONCLUSIONS: We demonstrate that a simple cfDNA genome-wide copy number approach provides an effective means of monitoring patients through treatment and show that targeted cfDNA sequencing identifies potential therapeutic targets in more than 50% of patients. We are now incorporating this approach into additional studies and trials of targeted therapies.

Publisher Correction: In silico prediction of housekeeping long intergenic non-coding RNAs reveals HKlincR1 as an essential player in lung cancer cell survival.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Pulmonary venous circulating tumor cell dissemination before tumor resection and disease relapse.

Approximately 50% of patients with early-stage non-small-cell lung cancer (NSCLC) who undergo surgery with curative intent will relapse within 5 years1,2. Detection of circulating tumor cells (CTCs) at the time of surgery may represent a tool to identify patients at higher risk of recurrence for whom more frequent monitoring is advised. Here we asked whether CellSearch-detected pulmonary venous CTCs (PV-CTCs) at surgical resection of early-stage NSCLC represent subclones responsible for subsequent disease relapse. PV-CTCs were detected in 48% of 100 patients enrolled into the TRACERx study3, were associated with lung-cancer-specific relapse and remained an independent predictor of relapse in multivariate analysis adjusted for tumor stage. In a case study, genomic profiling of single PV-CTCs collected at surgery revealed higher mutation overlap with metastasis detected 10 months later (91%) than with the primary tumor (79%), suggesting that early-disseminating PV-CTCs were responsible for disease relapse. Together, PV-CTC enumeration and genomic profiling highlight the potential of PV-CTCs as early predictors of NSCLC recurrence after surgery. However, the limited sensitivity of PV-CTCs in predicting relapse suggests that further studies using a larger, independent cohort are warranted to confirm and better define the potential clinical utility of PV-CTCs in early-stage NSCLC.

In silico prediction of housekeeping long intergenic non-coding RNAs reveals HKlincR1 as an essential player in lung cancer cell survival.

Prioritising long intergenic noncoding RNAs (lincRNAs) for functional characterisation is a significant challenge. Here we applied computational approaches to discover lincRNAs expected to play a critical housekeeping (HK) role within the cell. Using the Illumina Human BodyMap RNA sequencing dataset as a starting point, we first identified lincRNAs ubiquitously expressed across a panel of human tissues. This list was then further refined by reference to conservation score, secondary structure and promoter DNA methylation status. Finally, we used tumour expression and copy number data to identify lincRNAs rarely downregulated or deleted in multiple tumour types. The resulting list of candidate essential lincRNAs was then subjected to co-expression analyses using independent data from ENCODE and The Cancer Genome Atlas (TCGA). This identified a substantial subset with a predicted role in DNA replication and cell cycle regulation. One of these, HKlincR1, was selected for further characterisation. Depletion of HKlincR1 affected cell growth in multiple lung cancer cell lines, and led to disruption of genes involved in cell growth and viability. In addition, HKlincR1 expression was correlated with overall survival in lung adenocarcinoma patients. Our in silico studies therefore reveal a set of housekeeping noncoding RNAs of interest both in terms of their role in normal homeostasis, and their relevance in tumour growth and maintenance.

Genomic Evaluation of Multiparametric Magnetic Resonance Imaging-visible and -nonvisible Lesions in Clinically Localised Prostate Cancer.

BACKGROUND: The prostate cancer (PCa) diagnostic pathway is undergoing a radical change with the introduction of multiparametric magnetic resonance imaging (mpMRI), genomic testing, and different prostate biopsy techniques. It has been proposed that these tests should be used in a sequential manner to optimise risk stratification. OBJECTIVE: To characterise the genomic, epigenomic, and transcriptomic features of mpMRI-visible and -nonvisible PCa in clinically localised disease. DESIGN, SETTING, AND PARTICIPANTS: Multicore analysis of fresh prostate tissue sampled immediately after radical prostatectomy was performed for intermediate- to high-risk PCa. INTERVENTION: Low-pass whole-genome, exome, methylation, and transcriptome profiling of patient tissue cores taken from microscopically benign and cancerous areas in the same prostate. Circulating free and germline DNA was assessed from the blood of five patients. OUTCOME MEASUREMENT AND STATISTICAL ANALYSIS: Correlations between preoperative mpMRI and genomic characteristics of tumour and benign prostate samples were assessed. Gene profiles for individual tumour cores were correlated with existing genomic classifiers currently used for prognostication. RESULTS AND LIMITATIONS: A total of 43 prostate cores (22 tumour and 21 benign) were profiled from six whole prostate glands. Of the 22 tumour cores, 16 were tumours visible and six were tumours nonvisible on mpMRI. Intratumour genomic, epigenomic, and transcriptomic heterogeneity was found within mpMRI-visible lesions. This could potentially lead to misclassification of patients using signatures based on copy number or RNA expression. Moreover, three of the six cores obtained from mpMRI-nonvisible tumours harboured one or more genetic alterations commonly observed in metastatic castration-resistant PCa. No circulating free DNA alterations were found. Limitations include the small cohort size and lack of follow-up. CONCLUSIONS: Our study supports the continued use of systematic prostate sampling in addition to mpMRI, as avoidance of systematic biopsies in patients with negative mpMRI may mean that clinically significant tumours harbouring genetic alterations commonly seen in metastatic PCa are missed. Furthermore, there is inconsistency in individual genomics when genomic classifiers are applied. PATIENT SUMMARY: Our study shows that tumour heterogeneity within prostate tumours visible on multiparametric magnetic resonance imaging (mpMRI) can lead to misclassification of patients if only one core is used for genomic analysis. In addition, some cancers that were missed by mpMRI had genomic aberrations that are commonly seen in advanced metastatic prostate cancer. Avoiding biopsies in mpMRI-negative cases may mean that such potentially lethal cancers are missed.

Utility of ctDNA to support patient selection for early phase clinical trials: the TARGET study.

Next-generation sequencing (NGS) of circulating tumor DNA (ctDNA) supports blood-based genomic profiling but is not yet routinely implemented in the setting of a phase I trials clinic. TARGET is a molecular profiling program with the primary aim to match patients with a broad range of advanced cancers to early phase clinical trials on the basis of analysis of both somatic mutations and copy number alterations (CNA) across a 641 cancer-associated-gene panel in a single ctDNA assay. For the first 100 TARGET patients, ctDNA data showed good concordance with matched tumor and results were turned round within a clinically acceptable timeframe for Molecular Tumor Board (MTB) review. When a 2.5% variant allele frequency (VAF) threshold was applied, actionable mutations were identified in 41 of 100 patients, and 11 of these patients received a matched therapy. These data support the application of ctDNA in this early phase trial setting where broad genomic profiling of contemporaneous tumor material enhances patient stratification to novel therapies and provides a practical template for bringing routinely applied blood-based analyses to the clinic.

Ex vivo culture of cells derived from circulating tumour cell xenograft to support small cell lung cancer research and experimental therapeutics.

BACKGROUND AND PURPOSE: Small cell lung cancer (SCLC) is an aggressive disease with median survival of <2 years. Tumour biopsies for research are scarce, especially from extensive-stage patients, with repeat sampling at disease progression rarely performed. We overcame this limitation for relevant preclinical models by developing SCLC circulating tumour cell derived explants (CDX), which mimic the donor tumour pathology and chemotherapy response. To facilitate compound screening and identification of clinically relevant biomarkers, we developed short-term ex vivo cultures of CDX tumour cells. EXPERIMENTAL APPROACH: CDX tumours were disaggregated, and the human tumour cells derived were cultured for a maximum of 5 weeks. Phenotypic, transcriptomic and pharmacological characterization of these cells was performed. KEY RESULTS: CDX cultures maintained a neuroendocrine phenotype, and most changes in the expression of protein-coding genes observed in cultures, for up to 4 weeks, were reversible when the cells were re-implanted in vivo. Moreover, the CDX cultures exhibited a similar sensitivity to chemotherapy compared to the corresponding CDX tumour in vivo and were able to predict in vivo responses to therapeutic candidates. CONCLUSIONS AND IMPLICATIONS: Short-term cultures of CDX provide a tractable platform to screen new treatments, identify predictive and pharmacodynamic biomarkers and investigate mechanisms of resistance to better understand the progression of this recalcitrant tumour.