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.

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.

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.

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.

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.

Mutation pattern analysis reveals polygenic mini-drivers associated with relapse after surgery in lung adenocarcinoma.

The genomic lesions found in malignant tumours exhibit a striking degree of heterogeneity. Many tumours lack a known driver mutation, and their genetic basis is unclear. By mapping the somatic mutations identified in primary lung adenocarcinomas onto an independent coexpression network derived from normal tissue, we identify a critical gene network enriched for metastasis-associated genes. While individual genes within this module were rarely mutated, a significant accumulation of mutations within this geneset was predictive of relapse in lung cancer patients that have undergone surgery. Since it is the density of mutations within this module that is informative, rather than the status of any individual gene, these data are in keeping with a 'mini-driver' model of tumorigenesis in which multiple mutations, each with a weak effect, combine to form a polygenic driver with sufficient power to significantly alter cell behaviour and ultimately patient outcome. These polygenic mini-drivers therefore provide a means by which heterogeneous mutation patterns can generate the consistent hallmark changes in phenotype observed across tumours.

Signaling pathway screening platforms are an efficient approach to identify therapeutic targets in cancers that lack known driver mutations: a case report for a cancer of unknown primary origin.

Precision medicine aims to tailor cancer therapies to target specific tumor-promoting aberrations. For tumors that lack actionable drivers, which occurs frequently in the clinic, extensive molecular characterization and pre-clinical drug efficacy studies will be required. A cell line maintained at low passage and a patient- derived xenograft model (PDX) were generated using a fresh biopsy from a patient with a poorly-differentiated neuroendocrine tumor of unknown primary origin. Next-generation sequencing, high throughput signaling network analysis, and drug efficacy trials were then conducted to identify actionable targets for therapeutic intervention. No actionable mutations were identified after whole exome sequencing of the patient's DNA. However, whole genome sequencing revealed amplification of the 3q and 5p chromosomal arms, that include the PIK3CA and RICTOR genes, respectively. We then conducted pathway analysis, which revealed activation of the AKT pathway. Based on this analysis, efficacy of PIK3CA and AKT inhibitors were evaluated in the tumor biopsy-derived cell culture and PDX, and response to the AKT inhibitor AZD5363 was observed both in vitro and in vivo indicating the patient would benefit from targeted therapies directed against the serine/threonine kinase AKT. In conclusion, our study demonstrates that high throughput signaling pathway analysis will significantly aid in identifying actionable alterations in rare tumors and guide patient stratification into early-phase clinical trials.

Truncation- and motif-based pan-cancer analysis reveals tumor-suppressing kinases.

A major challenge in cancer genomics is identifying "driver" mutations from the many neutral "passenger" mutations within a given tumor. To identify driver mutations that would otherwise be lost within mutational noise, we filtered genomic data by motifs that are critical for kinase activity. In the first step of our screen, we used data from the Cancer Cell Line Encyclopedia and The Cancer Genome Atlas to identify kinases with truncation mutations occurring within or before the kinase domain. The top 30 tumor-suppressing kinases were aligned, and hotspots for loss-of-function (LOF) mutations were identified on the basis of amino acid conservation and mutational frequency. The functional consequences of new LOF mutations were biochemically validated, and the top 15 hotspot LOF residues were used in a pan-cancer analysis to define the tumor-suppressing kinome. A ranked list revealed MAP2K7, an essential mediator of the c-Jun N-terminal kinase (JNK) pathway, as a candidate tumor suppressor in gastric cancer, despite its mutational frequency falling within the mutational noise for this cancer type. The majority of mutations in MAP2K7 abolished its catalytic activity, and reactivation of the JNK pathway in gastric cancer cells harboring LOF mutations in MAP2K7 or the downstream kinase JNK suppressed clonogenicity and growth in soft agar, demonstrating the functional relevance of inactivating the JNK pathway in gastric cancer. Together, our data highlight a broadly applicable strategy to identify functional cancer driver mutations and define the JNK pathway as tumor-suppressive in gastric cancer.

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.

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.

Recurrent MLK4 Loss-of-Function Mutations Suppress JNK Signaling to Promote Colon Tumorigenesis.

MLK4 is a member of the mixed-lineage family of kinases that regulate the JNK, p38, and ERK kinase signaling pathways. MLK4 mutations have been identified in various human cancers, including frequently in colorectal cancer, where their function and pathobiological importance have been uncertain. In this study, we assessed the functional consequences of MLK4 mutations in colon tumorigenesis. Biochemical data indicated that a majority of MLK4 mutations are loss-of-function (LOF) mutations that can exert dominant-negative effects. In seeking to understand the abrogated activity of these mutants, we elucidated a new MLK4 catalytic domain structure. To determine whether MLK4 is required to maintain tumorigenic phenotypes, we reconstituted its signaling axis in colon cancer cells harboring MLK4-inactivating mutations. We found that restoring MLK4 activity reduced cell viability, proliferation, and colony formation in vitro and delayed tumor growth in vivo. Mechanistic investigations established that restoring the function of MLK4 selectively induced the JNK pathway and its downstream targets, cJUN, ATF3, and the cyclin-dependent kinase inhibitors CDKN1A and CDKN2B. Our work indicates that MLK4 is a novel tumor-suppressing kinase harboring frequent LOF mutations that lead to diminished signaling in the JNK pathway and enhanced proliferation in colon cancer.

Hypoxia-driven splicing into noncoding isoforms regulates the DNA damage response.

Tumour hypoxia is associated with poor patient outcome and resistance to therapy. It is accompanied by widespread changes in gene expression mediated largely through the transcription factors HIF1/2/3α. Hypoxia impacts on multiple pathways throughout the cell and has widespread effects on phenotype. Here we use sample-specific annotation approaches to determine the changes in transcript architecture that arise as result of alternative splicing in hypoxic cells. Using in vivo data generated from a time course in reduced oxygenation we identified genome-wide switching between coding and noncoding isoforms, including a significant number of components of the DNA damage response pathway. Notably, HDAC6, a master regulator of the cytotoxic response, and TP53BP1, which sits at the nexus of the double-strand break repair pathway, both underwent a marked transition towards an intron-retention pattern with a concomitant decline in protein levels. These transitions from coding to noncoding isoforms were recapitulated in a large and independent cohort of 499 colorectal samples taken from The Cancer Genome Atlas (TCGA). The set of altered genes was enriched for multiple components of the Fanconi Anaemia, nucleotide excision and double-strand break repair pathways, and together correlating with tumour status at last contact. Altogether, these data demonstrate a new role for hypoxia-driven alternative splicing in regulating DNA damage response, and highlight the importance of considering alternative splicing as a critical factor in our understanding of human disease.

Using large-scale genomics data to identify driver mutations in lung cancer: methods and challenges.

Lung cancer is the commonest cause of cancer death in the world and carries a poor prognosis for most patients. While precision targeting of mutated proteins has given some successes for never- and light-smoking patients, there are no proven targeted therapies for the majority of smokers with the disease. Despite sequencing hundreds of lung cancers, known driver mutations are lacking for a majority of tumors. Distinguishing driver mutations from inconsequential passenger mutations in a given lung tumor is extremely challenging due to the high mutational burden of smoking-related cancers. Here we discuss the methods employed to identify driver mutations from these large datasets. We examine different approaches based on bioinformatics, in silico structural modeling and biological dependency screens and discuss the limitations of these approaches.

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.

Discrepancies in cancer genomic sequencing highlight opportunities for driver mutation discovery.

Cancer genome sequencing is being used at an increasing rate to identify actionable driver mutations that can inform therapeutic intervention strategies. A comparison of two of the most prominent cancer genome sequencing databases from different institutes (Cancer Cell Line Encyclopedia and Catalogue of Somatic Mutations in Cancer) revealed marked discrepancies in the detection of missense mutations in identical cell lines (57.38% conformity). The main reason for this discrepancy is inadequate sequencing of GC-rich areas of the exome. We have therefore mapped over 400 regions of consistent inadequate sequencing (cold-spots) in known cancer-causing genes and kinases, in 368 of which neither institute finds mutations. We demonstrate, using a newly identified PAK4 mutation as proof of principle, that specific targeting and sequencing of these GC-rich cold-spot regions can lead to the identification of novel driver mutations in known tumor suppressors and oncogenes. We highlight that cross-referencing between genomic databases is required to comprehensively assess genomic alterations in commonly used cell lines and that there are still significant opportunities to identify novel drivers of tumorigenesis in poorly sequenced areas of the exome. Finally, we assess other reasons for the observed discrepancy, such as variations in dbSNP filtering and the acquisition/loss of mutations, to give explanations as to why there is a discrepancy in pharmacogenomic studies, given recent concerns with poor reproducibility of data.

Epithelial and stromal microRNA signatures of columnar cell hyperplasia linking Let-7c to precancerous and cancerous breast cancer cell proliferation.

Columnar cell hyperplasia (CCH) is the earliest histologically identifiable breast lesion linked to cancer progression and is characterized by increased proliferation, decreased apoptosis and elevated oestrogen receptor α (ERα) expression. The mechanisms underlying the initiation of these lesions have not been clarified but might involve early and fundamental changes in cancer progression. MiRNAs are key regulators of several biological processes, acting by influencing the post-transcriptional regulation of numerous targets, thus making miRNAs potential candidates in cancer initiation. Here we have defined novel epithelial as well as stromal miRNA signatures from columnar cell hyperplasia lesions compared to normal terminal duct lobular units by using microdissection and miRNA microarrays. Let-7c were among the identified downregulated epithelial miRNAs and its functions were delineated in unique CCH derived cells and breast cancer cell line MCF-7 suggesting anti-proliferative traits potentially due to effects on Myb and ERα. MiR-132 was upregulated in the stroma surrounding CCH compared to stoma surrounding normal terminal duct lobular units (TDLUs), and overexpression of miR-132 in immortalized fibroblasts and in fibroblasts co-cultured with epithelial CCH cells caused substantial expression changes of genes involved in metabolism, DNA damage and cell motility. The miRNA signatures identified in CCH indicate early changes in the epithelial and stromal compartment of CCH and could represent early key alterations in breast cancer progression that potentially could be targeted in novel prevention or treatment schedules.

Tumorigenicity and genetic profiling of circulating tumor cells in small-cell lung cancer.

Small-cell lung cancer (SCLC), an aggressive neuroendocrine tumor with early dissemination and dismal prognosis, accounts for 15-20% of lung cancer cases and ∼200,000 deaths each year. Most cases are inoperable, and biopsies to investigate SCLC biology are rarely obtainable. Circulating tumor cells (CTCs), which are prevalent in SCLC, present a readily accessible 'liquid biopsy'. Here we show that CTCs from patients with either chemosensitive or chemorefractory SCLC are tumorigenic in immune-compromised mice, and the resultant CTC-derived explants (CDXs) mirror the donor patient's response to platinum and etoposide chemotherapy. Genomic analysis of isolated CTCs revealed considerable similarity to the corresponding CDX. Most marked differences were observed between CDXs from patients with different clinical outcomes. These data demonstrate that CTC molecular analysis via serial blood sampling could facilitate delivery of personalized medicine for SCLC. CDXs are readily passaged, and these unique mouse models provide tractable systems for therapy testing and understanding drug resistance mechanisms.

A global non-coding RNA system modulates fission yeast protein levels in response to stress.

Non-coding RNAs (ncRNAs) are frequent and prevalent across the taxa. Although individual non-coding loci have been assigned a function, most are uncharacterized. Their global biological significance is unproven and remains controversial. Here we investigate the role played by ncRNAs in the stress response of Schizosaccharomyces pombe. We integrate global proteomics and RNA sequencing data to identify a systematic programme in which elevated antisense RNA arising both from ncRNAs and from 3'-overlapping convergent gene pairs is directly associated with substantial reductions in protein levels throughout the genome. We describe an extensive array of ncRNAs with trans associations that have the potential to influence multiple pathways. Deletion of one such locus reduces levels of atf1, a transcription factor downstream of the stress-activated mitogen-activated protein kinase (MAPK) pathway, and alters sensitivity to oxidative stress. These non-coding transcripts therefore regulate specific stress responses, adding unanticipated information-processing capacity to the MAPK signalling system.

Investigation of radiosensitivity gene signatures in cancer cell lines.

Intrinsic radiosensitivity is an important factor underlying radiotherapy response, but there is no method for its routine assessment in human tumours. Gene signatures are currently being derived and some were previously generated by expression profiling the NCI-60 cell line panel. It was hypothesised that focusing on more homogeneous tumour types would be a better approach. Two cell line cohorts were used derived from cervix [n = 16] and head and neck [n = 11] cancers. Radiosensitivity was measured as surviving fraction following irradiation with 2 Gy (SF2) by clonogenic assay. Differential gene expression between radiosensitive and radioresistant cell lines (SF2 median) was investigated using Affymetrix GeneChip Exon 1.0ST (cervix) or U133A Plus2 (head and neck) arrays. There were differences within cell line cohorts relating to tissue of origin reflected by expression of the stratified epithelial marker p63. Of 138 genes identified as being associated with SF2, only 2 (1.4%) were congruent between the cervix and head and neck carcinoma cell lines (MGST1 and TFPI), and these did not partition the published NCI-60 cell lines based on SF2. There was variable success in applying three published radiosensitivity signatures to our cohorts. One gene signature, originally trained on the NCI-60 cell lines, did partially separate sensitive and resistant cell lines in all three cell line datasets. The findings do not confirm our hypothesis but suggest that a common transcriptional signature can reflect the radiosensitivity of tumours of heterogeneous origins.