Identification of transcriptional programs using dense vector representations defined by mutual information with GeneVector.

Deciphering individual cell phenotypes from cell-specific transcriptional processes requires high dimensional single cell RNA sequencing. However, current dimensionality reduction methods aggregate sparse gene information across cells, without directly measuring the relationships that exist between genes. By performing dimensionality reduction with respect to gene co-expression, low-dimensional features can model these gene-specific relationships and leverage shared signal to overcome sparsity. We describe GeneVector, a scalable framework for dimensionality reduction implemented as a vector space model using mutual information between gene expression. Unlike other methods, including principal component analysis and variational autoencoders, GeneVector uses latent space arithmetic in a lower dimensional gene embedding to identify transcriptional programs and classify cell types. In this work, we show in four single cell RNA-seq datasets that GeneVector was able to capture phenotype-specific pathways, perform batch effect correction, interactively annotate cell types, and identify pathway variation with treatment over time.

Ovarian cancer mutational processes drive site-specific immune evasion.

High-grade serous ovarian cancer (HGSOC) is an archetypal cancer of genomic instability1-4 patterned by distinct mutational processes5,6, tumour heterogeneity7-9 and intraperitoneal spread7,8,10. Immunotherapies have had limited efficacy in HGSOC11-13, highlighting an unmet need to assess how mutational processes and the anatomical sites of tumour foci determine the immunological states of the tumour microenvironment. Here we carried out an integrative analysis of whole-genome sequencing, single-cell RNA sequencing, digital histopathology and multiplexed immunofluorescence of 160 tumour sites from 42 treatment-naive patients with HGSOC. Homologous recombination-deficient HRD-Dup (BRCA1 mutant-like) and HRD-Del (BRCA2 mutant-like) tumours harboured inflammatory signalling and ongoing immunoediting, reflected in loss of HLA diversity and tumour infiltration with highly differentiated dysfunctional CD8+ T cells. By contrast, foldback-inversion-bearing tumours exhibited elevated immunosuppressive TGFß signalling and immune exclusion, with predominantly naive/stem-like and memory T cells. Phenotypic state associations were specific to anatomical sites, highlighting compositional, topological and functional differences between adnexal tumours and distal peritoneal foci. Our findings implicate anatomical sites and mutational processes as determinants of evolutionary phenotypic divergence and immune resistance mechanisms in HGSOC. Our study provides a multi-omic cellular phenotype data substrate from which to develop and interpret future personalized immunotherapeutic approaches and early detection research.

Single-cell genomic variation induced by mutational processes in cancer.

How cell-to-cell copy number alterations that underpin genomic instability1 in human cancers drive genomic and phenotypic variation, and consequently the evolution of cancer2, remains understudied. Here, by applying scaled single-cell whole-genome sequencing3 to wild-type, TP53-deficient and TP53-deficient;BRCA1-deficient or TP53-deficient;BRCA2-deficient mammary epithelial cells (13,818 genomes), and to primary triple-negative breast cancer (TNBC) and high-grade serous ovarian cancer (HGSC) cells (22,057 genomes), we identify three distinct 'foreground' mutational patterns that are defined by cell-to-cell structural variation. Cell- and clone-specific high-level amplifications, parallel haplotype-specific copy number alterations and copy number segment length variation (serrate structural variations) had measurable phenotypic and evolutionary consequences. In TNBC and HGSC, clone-specific high-level amplifications in known oncogenes were highly prevalent in tumours bearing fold-back inversions, relative to tumours with homologous recombination deficiency, and were associated with increased clone-to-clone phenotypic variation. Parallel haplotype-specific alterations were also commonly observed, leading to phylogenetic evolutionary diversity and clone-specific mono-allelic expression. Serrate variants were increased in tumours with fold-back inversions and were highly correlated with increased genomic diversity of cellular populations. Together, our findings show that cell-to-cell structural variation contributes to the origins of phenotypic and evolutionary diversity in TNBC and HGSC, and provide insight into the genomic and mutational states of individual cancer cells.

Rejuvenating dysfunctional T cells in ovarian cancer: CD28 is the license to kill.

High-grade serous ovarian cancers (HGSOCs) exhibit limited response to immune checkpoint blockade. In a new study in Cancer Cell, Duraiswamy et al. highlight intratumoral CD28 co-stimulation by myeloid-antigen-presenting cells as a key mechanism required for activation of programmed cell death receptor 1 (PD-1)+ tumor-infiltrating T lymphocytes during PD-1 blockade in HGSOC.

Single-cell RNA sequencing reveals distinct tumor microenvironmental patterns in lung adenocarcinoma.

Recent developments in immuno-oncology demonstrate that not only cancer cells, but also the tumor microenvironment can guide precision medicine. A comprehensive and in-depth characterization of the tumor microenvironment is challenging since its cell populations are diverse and can be important even if scarce. To identify clinically relevant microenvironmental and cancer features, we applied single-cell RNA sequencing to ten human lung adenocarcinomas and ten normal control tissues. Our analyses revealed heterogeneous carcinoma cell transcriptomes reflecting histological grade and oncogenic pathway activities, and two distinct microenvironmental patterns. The immune-activated CP²E microenvironment was composed of cancer-associated myofibroblasts, proinflammatory monocyte-derived macrophages, plasmacytoid dendritic cells and exhausted CD8+ T cells, and was prognostically unfavorable. In contrast, the inert N³MC microenvironment was characterized by normal-like myofibroblasts, non-inflammatory monocyte-derived macrophages, NK cells, myeloid dendritic cells and conventional T cells, and was associated with a favorable prognosis. Microenvironmental marker genes and signatures identified in single-cell profiles had progonostic value in bulk tumor profiles. In summary, single-cell RNA profiling of lung adenocarcinoma provides additional prognostic information based on the microenvironment, and may help to predict therapy response and to reveal possible target cell populations for future therapeutic approaches.

Mitogen-activated protein kinase activity drives cell trajectories in colorectal cancer.

In colorectal cancer, oncogenic mutations transform a hierarchically organized and homeostatic epithelium into invasive cancer tissue lacking visible organization. We sought to define transcriptional states of colorectal cancer cells and signals controlling their development by performing single-cell transcriptome analysis of tumors and matched non-cancerous tissues of twelve colorectal cancer patients. We defined patient-overarching colorectal cancer cell clusters characterized by differential activities of oncogenic signaling pathways such as mitogen-activated protein kinase and oncogenic traits such as replication stress. RNA metabolic labeling and assessment of RNA velocity in patient-derived organoids revealed developmental trajectories of colorectal cancer cells organized along a mitogen-activated protein kinase activity gradient. This was in contrast to normal colon organoid cells developing along graded Wnt activity. Experimental targeting of EGFR-BRAF-MEK in cancer organoids affected signaling and gene expression contingent on predictive KRAS/BRAF mutations and induced cell plasticity overriding default developmental trajectories. Our results highlight directional cancer cell development as a driver of non-genetic cancer cell heterogeneity and re-routing of trajectories as a response to targeted therapy.

Mutation-specific effects of NRAS oncogenes in colorectal cancer cells.

In colorectal cancer (CRC), the prevalence of NRAS mutations (5-9%) is inferior to that of KRAS mutations (40-50%). NRAS mutations feature lately during tumour progression and drive resistance to anti-EGFR therapy in KRAS wild-type tumours. To elucidate specific functions of NRAS mutations in CRC, we expressed doxycycline-inducible G12D and Q61K mutations in the CRC cell line Caco-2. A focused phospho-proteome analysis based on the Bio-Plex platform, which interrogated the activity of MAPK, PI3K, mTOR, STAT, p38, JNK and ATF2, did not reveal significant differences between Caco-2 cells expressing NRASG12D, NRASQ61K and KRASG12V. However, phenotypic read-outs were different. The NRAS Q61K mutation promoted anchorage-independent proliferation and tumorigenicity, similar to features driven by canonical KRAS mutations. In contrast, expression of NRASG12D resulted in reduced proliferation and apoptosis. At the transcriptome level, we saw upregulation of cytokines and chemokines. IL1A, IL11, CXCL8 (IL-8) and CCL20 exhibited enhanced secretion into the culture medium. In addition, RNA sequencing results indicated activation of the IL1-, JAK/STAT-, NFκB- and TNFα signalling pathways. These results form the basis for an NRASG12D-driven inflammatory phenotype in CRC.

SFPQ Depletion Is Synthetically Lethal with BRAFV600E in Colorectal Cancer Cells.

Oncoproteins such as the BRAFV600E kinase endow cancer cells with malignant properties, but they also create unique vulnerabilities. Targeting of BRAFV600E-driven cytoplasmic signaling networks has proved ineffective, as patients regularly relapse with reactivation of the targeted pathways. We identify the nuclear protein SFPQ to be synthetically lethal with BRAFV600E in a loss-of-function shRNA screen. SFPQ depletion decreases proliferation and specifically induces S-phase arrest and apoptosis in BRAFV600E-driven colorectal and melanoma cells. Mechanistically, SFPQ loss in BRAF-mutant cancer cells triggers the Chk1-dependent replication checkpoint, results in decreased numbers and reduced activities of replication factories, and increases collision between replication and transcription. We find that BRAFV600E-mutant cancer cells and organoids are sensitive to combinations of Chk1 inhibitors and chemically induced replication stress, pointing toward future therapeutic approaches exploiting nuclear vulnerabilities induced by BRAFV600E.

Cell type-dependent differential activation of ERK by oncogenic KRAS in colon cancer and intestinal epithelium.

Oncogenic mutations in KRAS or BRAF are frequent in colorectal cancer and activate the ERK kinase. Here, we find graded ERK phosphorylation correlating with cell differentiation in patient-derived colorectal cancer organoids with and without KRAS mutations. Using reporters, single cell transcriptomics and mass cytometry, we observe cell type-specific phosphorylation of ERK in response to transgenic KRASG12V in mouse intestinal organoids, while transgenic BRAFV600E activates ERK in all cells. Quantitative network modelling from perturbation data reveals that activation of ERK is shaped by cell type-specific MEK to ERK feed forward and negative feedback signalling. We identify dual-specificity phosphatases as candidate modulators of ERK in the intestine. Furthermore, we find that oncogenic KRAS, together with ß-Catenin, favours expansion of crypt cells with high ERK activity. Our experiments highlight key differences between oncogenic BRAF and KRAS in colorectal cancer and find unexpected heterogeneity in a signalling pathway with fundamental relevance for cancer therapy.

Perturbation-response genes reveal signaling footprints in cancer gene expression.

Aberrant cell signaling can cause cancer and other diseases and is a focal point of drug research. A common approach is to infer signaling activity of pathways from gene expression. However, mapping gene expression to pathway components disregards the effect of post-translational modifications, and downstream signatures represent very specific experimental conditions. Here we present PROGENy, a method that overcomes both limitations by leveraging a large compendium of publicly available perturbation experiments to yield a common core of Pathway RespOnsive GENes. Unlike pathway mapping methods, PROGENy can (i) recover the effect of known driver mutations, (ii) provide or improve strong markers for drug indications, and (iii) distinguish between oncogenic and tumor suppressor pathways for patient survival. Collectively, these results show that PROGENy accurately infers pathway activity from gene expression in a wide range of conditions.

A compendium of ERK targets.

The RAF-MEK-ERK cascade is one of the most studied signaling pathways as it controls many vital cellular programs. There has been an immense amount of effort to determine ERK target proteins involved in regulating these programs. Classical biochemical and genetic approaches have elicited hundreds of direct ERK substrates, and with the advent of phospho-proteomic technologies, numerous studies have expanded the number of ERK target proteins. Here, we compile a comprehensive ERK target phospho-site archive, in which we gathered information from various research studies, and we provide this archive as an online database to form a searchable compendium of ERK targets.

An immediate-late gene expression module decodes ERK signal duration.

The RAF-MEK-ERK signalling pathway controls fundamental, often opposing cellular processes such as proliferation and apoptosis. Signal duration has been identified to play a decisive role in these cell fate decisions. However, it remains unclear how the different early and late responding gene expression modules can discriminate short and long signals. We obtained both protein phosphorylation and gene expression time course data from HEK293 cells carrying an inducible construct of the proto-oncogene RAF By mathematical modelling, we identified a new gene expression module of immediate-late genes (ILGs) distinct in gene expression dynamics and function. We find that mRNA longevity enables these ILGs to respond late and thus translate ERK signal duration into response amplitude. Despite their late response, their GC-rich promoter structure suggested and metabolic labelling with 4SU confirmed that transcription of ILGs is induced immediately. A comparative analysis shows that the principle of duration decoding is conserved in PC12 cells and MCF7 cells, two paradigm cell systems for ERK signal duration. Altogether, our findings suggest that ILGs function as a gene expression module to decode ERK signal duration.