Oncogenic KRAS Regulates Tumor Cell Signaling via Stromal Reciprocation.

Oncogenic mutations regulate signaling within both tumor cells and adjacent stromal cells. Here, we show that oncogenic KRAS (KRAS(G12D)) also regulates tumor cell signaling via stromal cells. By combining cell-specific proteome labeling with multivariate phosphoproteomics, we analyzed heterocellular KRAS(G12D) signaling in pancreatic ductal adenocarcinoma (PDA) cells. Tumor cell KRAS(G12D) engages heterotypic fibroblasts, which subsequently instigate reciprocal signaling in the tumor cells. Reciprocal signaling employs additional kinases and doubles the number of regulated signaling nodes from cell-autonomous KRAS(G12D). Consequently, reciprocal KRAS(G12D) produces a tumor cell phosphoproteome and total proteome that is distinct from cell-autonomous KRAS(G12D) alone. Reciprocal signaling regulates tumor cell proliferation and apoptosis and increases mitochondrial capacity via an IGF1R/AXL-AKT axis. These results demonstrate that oncogene signaling should be viewed as a heterocellular process and that our existing cell-autonomous perspective underrepresents the extent of oncogene signaling in cancer. VIDEO ABSTRACT.

Cancer-associated protein kinase C mutations reveal kinase's role as tumor suppressor.

Protein kinase C (PKC) isozymes have remained elusive cancer targets despite the unambiguous tumor promoting function of their potent ligands, phorbol esters, and the prevalence of their mutations. We analyzed 8% of PKC mutations identified in human cancers and found that, surprisingly, most were loss of function and none were activating. Loss-of-function mutations occurred in all PKC subgroups and impeded second-messenger binding, phosphorylation, or catalysis. Correction of a loss-of-function PKCß mutation by CRISPR-mediated genome editing in a patient-derived colon cancer cell line suppressed anchorage-independent growth and reduced tumor growth in a xenograft model. Hemizygous deletion promoted anchorage-independent growth, revealing that PKCß is haploinsufficient for tumor suppression. Several mutations were dominant negative, suppressing global PKC signaling output, and bioinformatic analysis suggested that PKC mutations cooperate with co-occurring mutations in cancer drivers. These data establish that PKC isozymes generally function as tumor suppressors, indicating that therapies should focus on restoring, not inhibiting, PKC activity.

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.

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.

In silico error correction improves cfDNA mutation calling.

MOTIVATION: Circulating-free DNA (cfDNA) profiling by sequencing is an important minimally invasive protocol for monitoring the mutation profile of solid tumours in cancer patients. Since the concentration of available cfDNA is limited, sample library generation relies on multiple rounds of PCR amplification, during which the accumulation of errors results in reduced sensitivity and lower accuracy. RESULTS: We present PCR Error Correction (PEC), an algorithm to identify and correct errors in short read sequencing data. It exploits the redundancy that arises from multiple rounds of PCR amplification. PEC is particularly well suited to applications such as single-cell sequencing and circulating tumour DNA (ctDNA) analysis, in which many cycles of PCR are used to generate sufficient DNA for sequencing from small amounts of starting material. When applied to ctDNA analysis, PEC significantly improves mutation calling accuracy, achieving similar levels of performance to more complex strategies that require additional protocol steps and access to calibration DNA datasets. AVAILABILITY AND IMPLEMENTATION: PEC is available under the GPL-v3 Open Source licence, and is freely available from: https://github.com/CRUKMI-ComputationalBiology/PCR_Error_Correction.git. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Mouse RUNX1C regulates premegakaryocytic/erythroid output and maintains survival of megakaryocyte progenitors.

RUNX1 is crucial for the regulation of megakaryocyte specification, maturation, and thrombopoiesis. Runx1 possesses 2 promoters: the distal P1 and proximal P2 promoters. The major protein isoforms generated by P1 and P2 are RUNX1C and RUNX1B, respectively, which differ solely in their N-terminal amino acid sequences. RUNX1C is the most abundantly expressed isoform in adult hematopoiesis, present in all RUNX1-expressing populations, including the cKit+ hematopoietic stem and progenitor cells. RUNX1B expression is more restricted, being highly expressed in the megakaryocyte lineage but downregulated during erythropoiesis. We generated a Runx1 P1 knock-in of RUNX1B, termed P1-MRIPV This mouse line lacks RUNX1C expression but has normal total RUNX1 levels, solely comprising RUNX1B. Using this mouse line, we establish a specific requirement for the P1-RUNX1C isoform in megakaryopoiesis, which cannot be entirely compensated for by RUNX1B overexpression. P1 knock-in megakaryocyte progenitors have reduced proliferative capacity and undergo increased cell death, resulting in thrombocytopenia. P1 knock-in premegakaryocyte/erythroid progenitors demonstrate an erythroid-specification bias, evident from increased erythroid colony-forming ability and decreased megakaryocyte output. At a transcriptional level, multiple erythroid-specific genes are upregulated and megakaryocyte-specific transcripts are downregulated. In addition, proapoptotic pathways are activated in P1 knock-in premegakaryocyte/erythroid progenitors, presumably accounting for the increased cell death in the megakaryocyte progenitor compartment. Unlike in the conditional adult Runx1 null models, megakaryocytic maturation is not affected in the P1 knock-in mice, suggesting that RUNX1B can regulate endomitosis and thrombopoiesis. Therefore, despite the high degree of structural similarity, RUNX1B and RUNX1C isoforms have distinct and specific roles in adult megakaryopoiesis.

Correction: RUNX1B Expression Is Highly Heterogeneous and Distinguishes Megakaryocytic and Erythroid Lineage Fate in Adult Mouse Hematopoiesis.

[This corrects the article DOI: 10.1371/journal.pgen.1005814.].

RUNX1B Expression Is Highly Heterogeneous and Distinguishes Megakaryocytic and Erythroid Lineage Fate in Adult Mouse Hematopoiesis.

The Core Binding Factor (CBF) protein RUNX1 is a master regulator of definitive hematopoiesis, crucial for hematopoietic stem cell (HSC) emergence during ontogeny. RUNX1 also plays vital roles in adult mice, in regulating the correct specification of numerous blood lineages. Akin to the other mammalian Runx genes, Runx1 has two promoters P1 (distal) and P2 (proximal) which generate distinct protein isoforms. The activities and specific relevance of these two promoters in adult hematopoiesis remain to be fully elucidated. Utilizing a dual reporter mouse model we demonstrate that the distal P1 promoter is broadly active in adult hematopoietic stem and progenitor cell (HSPC) populations. By contrast the activity of the proximal P2 promoter is more restricted and its upregulation, in both the immature Lineage- Sca1high cKithigh (LSK) and bipotential Pre-Megakaryocytic/Erythroid Progenitor (PreMegE) populations, coincides with a loss of erythroid (Ery) specification. Accordingly the PreMegE population can be prospectively separated into "pro-erythroid" and "pro-megakaryocyte" populations based on Runx1 P2 activity. Comparative gene expression analyses between Runx1 P2+ and P2- populations indicated that levels of CD34 expression could substitute for P2 activity to distinguish these two cell populations in wild type (WT) bone marrow (BM). Prospective isolation of these two populations will enable the further investigation of molecular mechanisms involved in megakaryocytic/erythroid (Mk/Ery) cell fate decisions. Having characterized the extensive activity of P1, we utilized a P1-GFP homozygous mouse model to analyze the impact of the complete absence of Runx1 P1 expression in adult mice and observed strong defects in the T cell lineage. Finally, we investigated how the leukemic fusion protein AML1-ETO9a might influence Runx1 promoter usage. Short-term AML1-ETO9a induction in BM resulted in preferential P2 upregulation, suggesting its expression may be important to establish a pre-leukemic environment.

twoddpcr: an R/Bioconductor package and Shiny app for Droplet Digital PCR analysis.

SUMMARY: Droplet Digital PCR (ddPCR) is a sensitive platform used to quantify specific nucleic acid molecules amplified by polymerase chain reactions. Its sensitivity makes it particularly useful for the detection of rare mutant molecules, such as those present in a sample of circulating free tumour DNA obtained from cancer patients. ddPCR works by partitioning a sample into individual droplets for which the majority contain only zero or one target molecule. Each droplet then becomes a reaction chamber for PCR, which through the use of fluorochrome labelled probes allows the target molecules to be detected by measuring the fluorescence intensity of each droplet. The technology supports two channels, allowing, for example, mutant and wild type molecules to be detected simultaneously in the same sample. As yet, no open source software is available for the automatic gating of two channel ddPCR experiments in the case where the droplets can be grouped into four clusters. Here, we present an open source R package 'twoddpcr', which uses Poisson statistics to estimate the number of molecules in such two channel ddPCR data. Using the Shiny framework, an accompanying graphical user interface (GUI) is also included for the package, allowing users to adjust parameters and see the results in real-time. AVAILABILITY AND IMPLEMENTATION: twoddpcr is available from Bioconductor (3.5) at https://bioconductor.org/packages/twoddpcr/ . A Shiny-based GUI suitable for non-R users is available as a standalone application from within the package and also as a web application at http://shiny.cruk.manchester.ac.uk/twoddpcr/ . CONTACT: ged.brady@cruk.manchester.ac.uk or crispin.miller@cruk.manchester.ac.uk. PACKAGE MAINTAINER: anthony.chiu@cruk.manchester.ac.uk.

RUNX1 positively regulates a cell adhesion and migration program in murine hemogenic endothelium prior to blood emergence.

During ontogeny, the transcription factor RUNX1 governs the emergence of definitive hematopoietic cells from specialized endothelial cells called hemogenic endothelium (HE). The ultimate consequence of this endothelial-to-hematopoietic transition is the concomitant activation of the hematopoietic program and downregulation of the endothelial program. However, due to the rare and transient nature of the HE, little is known about the initial role of RUNX1 within this population. We, therefore, developed and implemented a highly sensitive DNA adenine methyltransferase identification-based methodology, including a novel data analysis pipeline, to map early RUNX1 transcriptional targets in HE cells. This novel transcription factor binding site identification protocol should be widely applicable to other low abundance cell types and factors. Integration of the RUNX1 binding profile with gene expression data revealed an unexpected early role for RUNX1 as a positive regulator of cell adhesion- and migration-associated genes within the HE. This suggests that RUNX1 orchestrates HE cell positioning and integration prior to the release of hematopoietic cells. Overall, our genome-wide analysis of the RUNX1 binding and transcriptional profile in the HE provides a novel comprehensive resource of target genes that will facilitate the precise dissection of the role of RUNX1 in early blood development.

Targeted genetic dependency screen facilitates identification of actionable mutations in FGFR4, MAP3K9, and PAK5 in lung cancer.

Approximately 70% of patients with non-small-cell lung cancer present with late-stage disease and have limited treatment options, so there is a pressing need to develop efficacious targeted therapies for these patients. This remains a major challenge as the underlying genetic causes of ~50% of non-small-cell lung cancers remain unknown. Here we demonstrate that a targeted genetic dependency screen is an efficient approach to identify somatic cancer alterations that are functionally important. By using this approach, we have identified three kinases with gain-of-function mutations in lung cancer, namely FGFR4, MAP3K9, and PAK5. Mutations in these kinases are activating toward the ERK pathway, and targeted depletion of the mutated kinases inhibits proliferation, suppresses constitutive activation of downstream signaling pathways, and results in specific killing of the lung cancer cells. Genomic profiling of patients with lung cancer is ushering in an era of personalized medicine; however, lack of actionable mutations presents a significant hurdle. Our study indicates that targeted genetic dependency screens will be an effective strategy to elucidate somatic variants that are essential for lung cancer cell viability.

Evaluation and validation of a robust single cell RNA-amplification protocol through transcriptional profiling of enriched lung cancer initiating cells.

BACKGROUND: Although profiling of RNA in single cells has broadened our understanding of development, cancer biology and mechanisms of disease dissemination, it requires the development of reliable and flexible methods. Here we demonstrate that the EpiStem RNA-Amp™ methodology reproducibly generates microgram amounts of cDNA suitable for RNA-Seq, RT-qPCR arrays and Microarray analysis. RESULTS: Initial experiments compared amplified cDNA generated by three commercial RNA-Amplification protocols (Miltenyi µMACS™ SuperAmp™, NuGEN Ovation® One-Direct System and EpiStem RNA-Amp™) applied to single cell equivalent levels of RNA (25-50 pg) using Affymetrix arrays. The EpiStem RNA-Amp™ kit exhibited the highest sensitivity and was therefore chosen for further testing. A comparison of Affymetrix array data from RNA-Amp™ cDNA generated from single MCF7 and MCF10A cells to reference controls of unamplified cDNA revealed a high degree of concordance. To assess the flexibility of the amplification system single cell RNA-Amp™ cDNA was also analysed using RNA-Seq and high-density qPCR, and showed strong cross-platform correlations. To exemplify the approach we used the system to analyse RNA profiles of small populations of rare cancer initiating cells (CICs) derived from a NSCLC patient-derived xenograft. RNA-Seq analysis was able to identify transcriptional differences in distinct subsets of CIC, with one group potentially enriched for metastasis formation. Pathway analysis revealed that the distinct transcriptional signatures demonstrated in the CIC subpopulations were significantly correlated with published stem-cell and epithelial-mesenchymal transition signatures. CONCLUSIONS: The combined results confirm the sensitivity and flexibility of the RNA-Amp™ method and demonstrate the suitability of the approach for identifying clinically relevant signatures in rare, biologically important cell populations.

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.

RAC2, AEP, and ICAM1 expression are associated with CNS disease in a mouse model of pre-B childhood acute lymphoblastic leukemia.

We developed a murine model of CNS disease to obtain a better understanding of the pathogenesis of CNS involvement in pre-B-cell acute lymphoblastic leukemia (ALL). Semiquantitative proteomic discovery-based approaches identified unique expression of asparaginyl endopeptidase (AEP), intercellular adhesion molecule 1 (ICAM1), and ras-related C3 botulinum toxin substrate 2 (RAC2), among others, in an invasive pre-B-cell line that produced CNS leukemia in NOD-SCID mice. Targeting RAC2 significantly inhibited in vitro invasion and delayed disease onset in mice. Induced expression of RAC2 in cell lines with low/absent expression of AEP and ICAM1 did not result in an invasive phenotype or murine CNS disease. Flow cytometric analysis identified an enriched population of blast cells expressing ICAM1/lymphocyte function associated antigen-1 (LFA-1)/CD70 in the CD10(+)/CD19(+) fraction of bone marrow aspirates obtained from relapsed compared with normal controls and those with primary disease. CD10(+)/CD19(+) fractions obtained from relapsed patients also express RAC2 and give rise to CNS disease in mice. Our data suggest that combinations of processes are involved in the pathogenesis of CNS disease in pre-B-cell ALL, support a model in which CNS disease occurs as a result of external invasion, and suggest that targeting the processes of adhesion and invasion unique to pre-B cells may prevent recurrences within the CNS.

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.

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.

Programmed fluctuations in sense/antisense transcript ratios drive sexual differentiation in S. pombe.

Strand-specific RNA sequencing of S. pombe revealed a highly structured programme of ncRNA expression at over 600 loci. Waves of antisense transcription accompanied sexual differentiation. A substantial proportion of ncRNA arose from mechanisms previously considered to be largely artefactual, including improper 3' termination and bidirectional transcription. Constitutive induction of the entire spk1+, spo4+, dis1+ and spo6+ antisense transcripts from an integrated, ectopic, locus disrupted their respective meiotic functions. This ability of antisense transcripts to disrupt gene function when expressed in trans suggests that cis production at native loci during sexual differentiation may also control gene function. Consistently, insertion of a marker gene adjacent to the dis1+ antisense start site mimicked ectopic antisense expression in reducing the levels of this microtubule regulator and abolishing the microtubule-dependent 'horsetail' stage of meiosis. Antisense production had no impact at any of these loci when the RNA interference (RNAi) machinery was removed. Thus, far from being simply 'genome chatter', this extensive ncRNA landscape constitutes a fundamental component in the controls that drive the complex programme of sexual differentiation in S. pombe.

Gastrointestinal differentiation marker Cytokeratin 20 is regulated by homeobox gene CDX1.

CDX1 is a transcription factor that plays a key role in intestinal development and differentiation. However, the downstream targets of CDX1 are less well defined than those of its close homologue, CDX2. We report here the identification of downstream targets of CDX1 using microarray gene-expression analysis and other approaches. Keratin 20 (KRT20), a member of the intermediate filament and a well-known marker of intestinal differentiation, was initially identified as one of the genes likely to be directly regulated by CDX1. CDX1 and KRT20 mRNA expression were significantly correlated in a panel of 38 colorectal cancer cell lines. Deletion and mutation analysis of the KRT20 promoter showed that the minimum regulatory region for the control of KRT20 expression by CDX1 is within 246 bp upstream of the KRT20 transcription start site. ChIP analysis confirmed that CDX1 binds to the predicted CDX elements in this region of the KRT20 promoter in vivo. In addition, immunohistochemistry showed expression of CDX1 parallels that of KRT20 in the normal crypt, which further supports their close relationship. In summary, our observations strongly imply that KRT20 is directly regulated by CDX1, and therefore suggest a role for CDX1 in maintaining differentiation in intestinal epithelial cells. Because a key feature of the development of a cancer is an unbalanced program of proliferation and differentiation, dysregulation of CDX1 may be an advantage for the development of a colorectal carcinoma. This could, therefore, explain the relatively frequent down regulation of CDX1 in colorectal carcinomas by hypermethylation.

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.