Search

Signaling between cancer and nonmalignant (stromal) cells in the tumor microenvironment (TME) is a key to tumor progression. Here, we deconvoluted bulk tumor transcriptomes to infer cross-talk between ligands and receptors on cancer and stromal cells in the TME of 20 solid tumor types. This approach recovered known transcriptional hallmarks of cancer and stromal cells and was concordant with single-cell, in situ hybridization and IHC data. Inferred autocrine cancer cell interactions varied between tissues but often converged on Ephrin, BMP, and FGFR-signaling pathways. Analysis of immune checkpoints nominated interactions with high levels of cancer-to-immune cross-talk across distinct tumor types. Strikingly, PD-L1 was found to be highly expressed in stromal rather than cancer cells. Overall, our study presents a new resource for hypothesis generation and exploration of cross-talk in the TME.Significance:This study provides deconvoluted bulk tumor transcriptomes across multiple cancer types to infer cross-talk in the tumor microenvironment.

Supplementary results and methods

Chromatin is organised in the form of domains known to be insulated from neighbouring regions in the genome. The colocalisation of genes sharing similar functions in the topologically associated domains (TAD) and the orchestration of their activity in cancer-cells and their drug-response have not been comprehensively explored. Here we analyzed patterns in the activity of TADs in cancer-cells along with chromatin interaction profiles to understand their modes of drug-response. Using our approach on the transcriptome of 819 cancer cell lines revealed that for multiple drugs, IC50 values were substantially more correlated with the activity of a few TADs than any of the genes. Such a result signifies the potential utility of TAD activity in highlighting common mechanisms of drug response in spite of heterogeneity in the association of genes with drug-resistance. Further we performed analysis of 9014 cancer sample transcriptomes from 20 cancer types to find the relationship between TAD-activity and survival to highlight context of cause and actionable mechanism for better therapy. Overall our analysis highlighted higher importance and influence of combinatorics of genes in many TADs in comparison to independent effect of their own genes, on the patient-survival and drug-response of cancer-cells.

Solid tumors are generally characterized by an acidic tumor microenvironment (TME) that favors cancer progression, therapy resistance and immune evasion. By single-cell RNA-sequencing analysis in individuals with pancreatic ductal adenocarcinoma (PDAC), we reveal solute carrier family 4 member 4 (SLC4A4) as the most abundant bicarbonate transporter, predominantly expressed by epithelial ductal cells. Functionally, SLC4A4 inhibition in PDAC cancer cells mitigates the acidosis of the TME due to bicarbonate accumulation in the extracellular space and a decrease in lactate production by cancer cells as the result of reduced glycolysis. In PDAC-bearing mice, genetic or pharmacological SLC4A4 targeting improves T cell-mediated immune response and breaches macrophage-mediated immunosuppression, thus inhibiting tumor growth and metastases. In addition, Slc4a4 targeting in combination with immune checkpoint blockade is able to overcome immunotherapy resistance and prolong survival. Overall, our data propose SLC4A4 as a therapeutic target to unleash an antitumor immune response in PDAC. ispartof: NATURE CANCER vol:3 issue:12 ispartof: location:England status: Published online

Cancer treatment is gradually evolving from the classical use of nonspecific cytotoxic drugs targeting generic mechanisms of cell growth and proliferation. Instead, new "patient-specific treatment paradigms" that are based on an individual patient's tumor-specific molecular features are emerging, and these include "druggable" genomic alterations such as oncogenic driver mutations, downstream activities of cancer-signaling pathways, and the expression of specific genes involved in tumorigenesis and cancer progression. This evolving landscape of making evidence-based treatment decisions forms the foundation of precision oncology, which aims to deliver "the right drug, to the right patient and at the right time". The long-term vision for this approach is to maximize the treatment efficacy while minimizing exposure to ineffective therapy and reducing co-morbidity-related side effects. Successful clinical translation and implementation of this vision have the potential to revolutionize treatment paradigms from predominantly reactive, to more evidence-based, proactive and predictive care. In this article, we review the past and current approaches in precision oncology, and describe their remarkable power and limitations. We also speculate on the evolution of newly emerging methodologies of the future that can be used to address some of the key challenges associated with the existing paradigms. This article is categorized under: CancerGenetics/Genomics/Epigenetics CancerMolecular and Cellular Physiology CancerComputational Models.

Drug resistance is a major obstacle to the treatment of most human tumors. In this study, we find that dual-specificity phosphatase 16 (DUSP16) regulates resistance to chemotherapy in nasopharyngeal carcinoma, colorectal cancer, gastric and breast cancer. Cancer cells expressing higher DUSP16 are intrinsically more resistant to chemotherapy-induced cell death than cells with lower DUSP16 expression. Overexpression of DUSP16 in cancer cells leads to increased resistance to cell death upon chemotherapy treatment. In contrast, knockdown of DUSP16 in cancer cells increases their sensitivity to treatment. Mechanistically, DUSP16 inhibits JNK and p38 activation, thereby reducing BAX accumulation in mitochondria to reduce apoptosis. Analysis of patient survival in head & neck cancer and breast cancer patient cohorts supports DUSP16 as a marker for sensitivity to chemotherapy and therapeutic outcome. This study therefore identifies DUSP16 as a prognostic marker for the efficacy of chemotherapy, and as a therapeutic target for overcoming chemoresistance in cancer., Chemoresistance is one of the main challenges for cancer therapy success. Here, the authors show that dual-specificity phosphatase 16 (DUSP16) expression is associated with chemoresistance in several types of cancer through impairing mitochondria-associated apoptosis.

The potential to migrate is one of the most fundamental functions for various epithelial, mesenchymal, and immune cells. Image analysis of motile cell populations, both primary and cultured, typically reveals an intercellular variability in migration speeds. However, cell migration chromatography, the sorting of large populations of cells based on their migratory characteristics, cannot be easily performed. The lack of such methods has hindered our understanding of the direct correlation between the capacity to migrate and other cellular properties. Here, we report two novel, easily implementable and readily scalable methods to sort millions of live migratory cancer and immune cells based on their spontaneous migration in two-dimensional and three-dimensional microenvironments, respectively. Correlative downstream transcriptomic, molecular and functional tests reveal marked differences between the fast and slow subpopulations in patient-derived cancer cells. We further employ our method to reveal that sorting the most migratory cytotoxic T lymphocytes yields a pool of cells with enhanced cytotoxicity against cancer cells. This phenotypic assay opens new avenues for the precise characterization of the mechanisms underlying hither to unexplained heterogeneities in migratory phenotypes within a cell population, and for the targeted enrichment of the most potent migratory leukocytes in immunotherapies.

Transcriptional reactivation of hTERT is the limiting step in tumorigenesis. While mutations in hTERT promoter present in 19% of cancers are recognized as key drivers of hTERT reactivation, mechanisms by which wildtype hTERT (WT-hTERT) promoter is reactivated, in majority of human cancers, remain unknown. Using primary colorectal cancers (CRC) we identified Tert INTeracting region 2 (T-INT2), the critical chromatin region essential for reactivating WT-hTERT promoter in CRCs. Elevated β-catenin and JunD level in CRC facilitates chromatin interaction between hTERT promoter and T-INT2 that is necessary to turn on hTERTexpression. Pharmacological screens uncovered salinomycin, which inhibits JunD mediated hTERT-T-INT2 interaction that is required for the formation of a stable transcription complex on the hTERT promoter. Our results showed for the first time how known CRC alterations, such as APC, lead to WT-hTERT promoter reactivation during stepwise-tumorigenesis and provide a new perspective for developing cancer-specific drugs.Healthy and cancer cells harbor the same DNA sequence, but reactivation of the Human Telomerase Reverse Transcriptase (hTERT) gene is observed only in cancer cells. How does that happen was not known for over three decades of research? This study identifies a specific DNA structure that forms only in cancer cells and brings the necessary molecular machinery into the correct position to activate the hTERT gene. The detailed mechanism of hTERT activation provided in this study will be instrumental in designing cancer cell-specific hTERT inhibitors, especially since all the other ways of inhibiting telomerase failed in the clinic.

Hepatitis B surface antigen (HBsAg) loss or functional cure (FC), is considered the desirable therapeutic outcome for chronic hepatitis B (CHB) patients. However, the immuno-pathological biomarkers and underlying mechanisms remain unclear. Here we present a comprehensive single cell-transcriptomic atlas together with immune-phenotyping of disease-associated cell states (DACS) isolated from intra-hepatic tissue and matched PBMCs of either CHB or FC patients. We find that the intra-hepatic environment displays specific cell identities and molecular signatures that are distinct from PBMCs. FC is associated with emergence of an altered adaptive immune response marked by CD4 cytotoxic T lymphocytes (CD4-CTLs), and an activated innate response represented by liver-resident natural killer (LR-NK) cells. Overall, these findings provide novel insights into immuno-pathological cell states associated with FC that could serve as prognostic biomarkers.

T-cell acute lymphoblastic leukemia (T-ALL) is an aggressive hematological malignancy, and novel therapeutics are much needed. Profiling patient leukemia’ drug sensitivities ex vivo, we discovered that 44.4% of childhood and 16.7% of adult T-ALL cases exquisitely respond to dasatinib. Applying network-based systems pharmacology analyses to examine signal circuitry, we identified preTCR-LCK activation as the driver of dasatinib sensitivity, and T-ALL-specific LCK dependency was confirmed in genome-wide CRISPR-Cas9 screens. Dasatinib-sensitive T-ALLs exhibited high BCL-XL and low BCL2 activity and venetoclax resistance. Discordant sensitivity of T-ALL to dasatinib and venetoclax is strongly correlated with T-cell differentiation, particularly with the dynamic shift in LCK vs. BCL2 activation. Finally, single-cell analysis identified leukemia heterogeneity in LCK and BCL2 signaling and T-cell maturation stage, consistent with dasatinib response. In conclusion, our results indicate that developmental arrest in T-ALL drives differential activation of preTCR-LCK and BCL2 signaling in this leukemia, providing unique opportunities for targeted therapy.

Significance Mutations in hTERT promoter are seen in over 19% of human cancers, irrespective of the cancer type. Understanding the molecular players that regulate Mut-hTERT promoters may help the design of effective targeting strategies to inhibit telomerase reactivation specifically in cancer cells. Our work uses genome-wide functional screens to identify 30 specific regulators of Mut-hTERT promoters. These candidates identified from the screening serve as an excellent resource to understand how telomerase is reactivated and as targets for making inhibitors to telomerase, a key driver of cancer., Cancer-specific hTERT promoter mutations reported in 19% of cancers result in enhanced telomerase activity. Understanding the distinctions between transcriptional regulation of wild-type (WT) and mutant (Mut) hTERT promoters may open up avenues for development of inhibitors which specially block hTERT expression in cancer cells. To comprehensively identify physiological regulators of WT- or Mut-hTERT promoters, we generated several isogenic reporter cells driven by endogenous hTERT loci. Genome-wide CRISPR-Cas9 and small interfering RNA screens using these isogenic reporter lines identified specific regulators of Mut-hTERT promoters. We validate and characterize one of these hits, namely, MED12, a kinase subunit of mediator complex. We demonstrate that MED12 specifically drives expression of hTERT from the Mut-hTERT promoter by mediating long-range chromatin interaction between the proximal Mut-hTERT promoter and T-INT1 distal regulatory region 260 kb upstream. Several hits identified in our screens could serve as potential therapeutic targets, inhibition of which may specifically block Mut-hTERT promoter driven telomerase reactivation in cancers.

Long non-coding RNAs (lncRNAs) regulate numerous biological processes involved in both development and carcinogenesis. Hippo-YAP/TAZ signaling, a critical pathway responsible for organ size control, is often dysregulated in a variety of cancers. However, the nature and function of YAP/TAZ-regulated lncRNAs during tumorigenesis remain largely unexplored. By profiling YAP/TAZ-regulated lncRNAs, we identified SFTA1P as a novel transcriptional target and a positive feedback regulator of YAP/TAZ signaling. Using non-small cell lung cancer (NSCLC) cell lines, we show that SFTA1P is transcriptionally activated by YAP/TAZ in a TEAD-dependent manner. Functionally, knockdown of SFTA1P in NSCLC cell lines inhibited proliferation, induced programmed cell death, and compromised their tumorigenic potential. Mechanistically, SFTA1P knockdown decreased TAZ protein abundance and consequently, the expression of YAP/TAZ transcriptional targets. We provide evidence that this phenomenon could potentially be mediated via its interaction with TAZ mRNA to regulate TAZ translation. Our results reveal SFTA1P as a positive feedback regulator of Hippo-YAP/TAZ signaling, which may serve as the molecular basis for lncRNA-based therapies against YAP/TAZ-driven cancers.

The consensus molecular subtype (CMS) classification of colorectal cancer is based on bulk transcriptomics. The underlying epithelial cell diversity remains unclear. We analyzed 373,058 single-cell transcriptomes from 63 patients, focusing on 49,155 epithelial cells. We identified a pervasive genetic and transcriptomic dichotomy of malignant cells, based on distinct gene expression, DNA copy number and gene regulatory network. We recapitulated these subtypes in bulk transcriptomes from 3,614 patients. The two intrinsic subtypes, iCMS2 and iCMS3, refine CMS. iCMS3 comprises microsatellite unstable (MSI-H) cancers and one-third of microsatellite-stable (MSS) tumors. iCMS3 MSS cancers are transcriptomically more similar to MSI-H cancers than to other MSS cancers. CMS4 cancers had either iCMS2 or iCMS3 epithelium; the latter had the worst prognosis. We defined the intrinsic epithelial axis of colorectal cancer and propose a refined 'IMF' classification with five subtypes, combining intrinsic epithelial subtype (I), microsatellite instability status (M) and fibrosis (F). ispartof: NATURE GENETICS vol:54 issue:7 pages:963-+ ispartof: location:United States status: published

The human genome contains more than 200,000 gene isoforms. However, different isoforms can be highly similar, and with an average length of 1.5kb remain difficult to study with short read sequencing. To systematically evaluate the ability to study the transcriptome at a resolution of individual isoforms we profiled 5 human cell lines with short read cDNA sequencing and Nanopore long read direct RNA, amplification-free direct cDNA, PCR-cDNA sequencing. The long read protocols showed a high level of consistency, with amplification-free RNA and cDNA sequencing being most similar. While short and long reads generated comparable gene expression estimates, they differed substantially for individual isoforms. We find that increased read length improves read-to-transcript assignment, identifies interactions between alternative promoters and splicing, enables the discovery of novel transcripts from repetitive regions, facilitates the quantification of full-length fusion isoforms and enables the simultaneous profiling of m6A RNA modifications when RNA is sequenced directly. Our study demonstrates the advantage of long read RNA sequencing and provides a comprehensive resource that will enable the development and benchmarking of computational methods for profiling complex transcriptional events at isoform-level resolution.

The recent coronavirus disease 2019 (COVID-19) pandemic is caused by severe acute respiratory syndrome coronavirus 2. COVID-19 was first reported in China (December 2019) and is now prevalent across the globe. Entry of severe acute respiratory syndrome coronavirus 2 into mammalian cells requires the binding of viral Spike (S) proteins to the angiotensin-converting enzyme 2 receptor. Once entered, the S protein is primed by a specialized serine protease, transmembrane serine protease 2 in the host cell. Importantly, besides the respiratory symptoms that are consistent with other common respiratory virus infections when patients become viremic, a significant number of COVID-19 patients also develop liver comorbidities. We explored whether a specific target cell-type in the mammalian liver could be implicated in disease pathophysiology other than the general deleterious response to cytokine storms. Here, we used single-cell RNA-seq to survey the human liver and identified potentially implicated liver cell-type for viral ingress. We analyzed ~300,000 single cells across five different (i.e., human fetal, healthy, cirrhotic, tumor, and adjacent normal) liver tissue types. This study reports on the co-expression of angiotensin-converting enzyme 2 and transmembrane serine protease 2 in a TROP2+ liver progenitor population. Importantly, we detected enrichment of this cell population in the cirrhotic liver when compared with tumor tissue. These results indicated that in COVID-19-associated liver dysfunction and cell death, a viral infection of TROP2+ progenitors in the liver might significantly impair liver regeneration in patients with liver cirrhosis.

Immune evasion is key to cancer initiation and later at metastasis, but its dynamics at intermediate stages, where potential therapeutic interventions could be applied, is undefined. Here we show, using multi-dimensional analyses of resected tumours, their adjacent non-tumour tissues and peripheral blood, that extensive immune remodelling takes place in patients with stage I to III hepatocellular carcinoma (HCC). We demonstrate the depletion of anti-tumoural immune subsets and accumulation of immunosuppressive or exhausted subsets along with reduced tumour infiltration of CD8 T cells peaking at stage II tumours. Corresponding transcriptomic modification occur in the genes related to antigen presentation, immune responses, and chemotaxis. The progressive immune evasion is validated in a murine model of HCC. Our results show evidence of ongoing tumour-immune co-evolution during HCC progression and offer insights into potential interventions to reverse, prevent or limit the progression of the disease.

Simple Summary Lung cancers lead cancer-related mortalities, with Non-Small Cell Lung Cancer (NSCLC) representing a substantial proportion of these cases. Perturbation of Hippo-YAP/TAZ signaling in NSCLC could be mainly attributed to post-transcriptional regulators since genetic alterations to the signaling pathway are known to be rare. In this study, we identified miR-582-5p as a novel, post-transcriptional regulator of Hippo-YAP/TAZ signaling. Our work revealed an inhibitory function of miR-582-5p on YAP/TAZ signaling in NSCLC cells, whereby the tumorigenic potential is diminished upon the overexpression of miR-582-5p. We also uncovered the regulation of miR-582-5p expression by YAP/TAZ, suggesting a potential feedback loop of YAP/TAZ signaling mediated by miR-582-5p. Mechanistically, we discovered NCKAP1 and PIP5K1C, regulators of actin polymerization, as novel and direct targets of miR-582-5p. Restoring miR-582-5p expression in NSCLC cells resulted in deficient F-actin, which potentially mediates the miR-582-5p-driven tumor suppression in a YAP/TAZ-dependent manner. Our findings underscore the anti-tumor function of miR-582-5p in NSCLC, positing its therapeutic potential in YAP/TAZ-driven lung cancers. Abstract The Hippo-YAP/TAZ signaling pathway is an evolutionarily conserved signaling pathway involved in a broad spectrum of biological processes, including tumorigenesis. Whilst aberrant Hippo-YAP/TAZ signaling is frequently reported in various cancers, the genetic alterations of this pathway are relatively rare, suggesting regulation at the post-transcriptional level. MicroRNAs play key role in tumorigenesis by regulating gene expression post-transcriptionally. Amongst the cancer-relevant microRNAs, miR-582-5p suppresses cell growth and tumorigenesis by inhibiting the expression of oncogenes, including AKT3, MAP3K2 and NOTCH1. Given the oncogenic role of YAP/TAZ in solid tumors, we scrutinized the possible interplay between miR-582-5p and Hippo-YAP/TAZ signaling. Correlation analysis in NSCLC cells revealed a positive relationship between the expression of mature miR-582-5p and the proportion of phosphorylated YAP/TAZ. Intriguingly, YAP/TAZ knockdown reduced the expression of mature miR-582-5p but increased that of primary miR-582. Overexpression of miR-582-5p resulted in increased phosphorylation of YAP/TAZ with a concomitant reduction in cell proliferation and enhanced apoptosis. Mechanistically, we find that miR-582-5p targets actin regulators NCKAP1 and PIP5K1C, which may be responsible for the observed alteration in F-actin, known to modulate YAP/TAZ. We postulate that regulation of the actin cytoskeleton by miR-582-5p may attenuate YAP/TAZ activity. Altogether, this study reveals a novel mechanism of YAP/TAZ regulation by miR-582-5p in a cytoskeleton-dependent manner and suggests a negative feedback loop, highlighting the therapeutic potential of restoring miR-582-5p expression in treating NSCLC.

Background Overexpression of epidermal growth factor receptor (EGFR), and downstream pathway activation appears to be a common oncogenic driver in the majority of head and neck squamous cell cancers (HNSCCs); yet targeting EGFR for the treatment of HNSCC has met with limited success. Apart from the anti-EGFR antibody cetuximab, no small molecule EGFR/tyrosine kinase inhibitors (TKIs) have progressed to routine clinical use. The aim of this study was to determine factors contributing to the lack of response to TKIs and identify alternative therapeutic vulnerabilities. Methods Genomic and transcriptomic sequencing, high-throughput compound screens, overexpression and siRNA knockdown, western blot, in vivo xenograft studies. Findings We derived three pairs of isogenic gefitinib (TKI)-sensitive and resistant patient-derived HNSCC cell lines. Genomic sequencing of gefitinib-resistant cell lines identified a lack of activating and resistance-associated EGFR mutations. Instead, transcriptomic sequencing showed upregulated EMT gene signature in the gefitinib-resistant cells with a corresponding increase in their migratory phenotype. Additionally, the resistant cell displayed reduced growth rate. Surprisingly, while gefitinib-resistant cells were independent of EGFR for survival, they nonetheless displayed activation of downstream ERK and AKT signalling. High-throughput screening (HTS) of druggable, small molecule libraries revealed that the gefitinib-resistant cells were particularly sensitive to inhibitors of genes involved in cell cycle and mitosis, such as Aurora kinase inhibitors (AKIs), cyclin-dependent kinase (CDK) inhibitors, and microtubule inhibitors. Notably our results showed that in the EGFR inhibited state, Aurora kinases are essential for cell survival. Interpretation Our study demonstrates that in the absence of activating EGFR mutations, HNSCCs may gain resistance to gefitinib through decreased cell proliferation, which makes them exceptionally vulnerable to cell-cycle inhibitors. Funding Agency for Science, Technology, and Research (A*STAR), National Medical Research Council (NMRC), and the National Institutes of Health (NIH)/National Cancer Institute (NCI).

SummaryWhile understanding heterogeneity in molecular signatures across patients underpins precision oncology, there is increasing appreciation for taking intra-tumor heterogeneity into account. Single-cell RNA-seq (scRNA-seq) technologies have facilitated investigations into the role of intra-tumor transcriptomic heterogeneity (ITTH) in tumor biology and evolution, but their application to in silico models of drug response has not been explored. Based on large-scale analysis of cancer omics datasets, we highlight the utility of ITTH for predicting clinical outcomes. We then show that heterogeneous gene expression signatures obtained from scRNA-seq data can be accurately analyzed (80%) in a recommender system framework (CaDRReS-Sc) for in silico drug response prediction. Patient-derived cell lines capturing transcriptomic heterogeneity from primary and metastatic tumors were used as in vitro proxies for validating monotherapy predictions (Pearson r>0.6), as well as optimal drug combinations to target different subclonal populations (>10% improvement). Applying CaDRReS-Sc to the increasing number of publicly available tumor scRNA-seq datasets can serve as an in silico screen for further in vitro and in vivo drug repurposing studies.Graphical abstractHighlightsLarge-scale analysis to establish the impact of transcriptomic heterogeneity within tumors on clinical outcomesCalibrated recommender system for drug response prediction based on single-cell RNA-seq data (CaDRReS-Sc)Prediction of drug response in patient-derived cell lines with transcriptomic heterogeneityIn silico identification of drug combinations that work based on clonal vulnerabilities

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

Transfer RNA (tRNA)-derived fragments (tRFs) are an emerging class of conserved small non-coding RNAs that play important roles in post-transcriptional gene regulation. High-throughput sequencing of multiple biological samples have identified heterogeneous species of tRFs with distinct functionalities. These small RNAs have garnered a lot of scientific attention due to their ubiquitous expression and versatility in regulating various biological processes. In this review, we highlight our current understanding of tRF biogenesis and their regulatory functions. We summarize the diverse modes of biogenesis through which tRFs are generated and discuss the mechanism through which different tRF species regulate gene expression and the biological implications. Finally, we conceptualize research areas that require focus to strengthen our understanding of the biogenesis and function of tRFs.

Signaling between cancer and nonmalignant (stromal) cells in the tumor microenvironment (TME) is a key to tumor progression. Here, we deconvoluted bulk tumor transcriptomes to infer cross-talk between ligands and receptors on cancer and stromal cells in the TME of 20 solid tumor types. This approach recovered known transcriptional hallmarks of cancer and stromal cells and was concordant with single-cell, in situ hybridization and IHC data. Inferred autocrine cancer cell interactions varied between tissues but often converged on Ephrin, BMP, and FGFR-signaling pathways. Analysis of immune checkpoints nominated interactions with high levels of cancer-to-immune cross-talk across distinct tumor types. Strikingly, PD-L1 was found to be highly expressed in stromal rather than cancer cells. Overall, our study presents a new resource for hypothesis generation and exploration of cross-talk in the TME. Significance: This study provides deconvoluted bulk tumor transcriptomes across multiple cancer types to infer cross-talk in the tumor microenvironment.

Whereas phenotypic assays such as Boyden chambers and wound healing assays can easily be employed to characterize the migratory potential of cells at the population level, few methods exist that can sort subpopulations of cells based on their migratory behaviour from an initial heterogeneous pool. In this paper, we present an approach to sort migratory cancer and immune cells based on their spontaneous migration in 2D and 3D microenvironments. Using this method, which is easy to implement and readily scalable, millions of live cells can be sorted based on their migratory characteristics and then subjected to downstream genomic, molecular and functional tests. We reveal that enrichment of the most migratory cytotoxic T lymphocytes yields a pool of cells with enhanced cytotoxicity against cancer cells. This new functional sorting method opens new avenues for the precise characterization of the mechanisms underlying hitherto unexplained heterogeneities in migratory phenotypes within a cell population, and for the targeted enrichment of the most potent migratory leukocytes in immunotherapies.

SummaryThe recent pandemic of coronavirus disease 2019 (COVID-19) is caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). COVID-19 was first reported in China (December 2019) and now prevalent in ∼170 countries across the globe. Entry of SARS-CoV-2 into mammalian cells require the binding of viral Spike (S) proteins to the ACE2 (angiotensin converting enzyme 2) receptor. Once entered the S protein is primed by a specialised serine protease, TMPRSS2 (Transmembrane Serine Protease 2) in the host cell. Importantly, beside respiratory symptoms, consistent with other common respiratory virus infection when patients become viraemic, a significant number of COVID-19 patients also develop liver comorbidities. We explored if specific target cell-type in the mammalian liver, could be implicated in disease pathophysiology other than the general deleterious response to cytokine storms. Here we employed single-cell RNA-seq (scRNA-seq) to survey the human liver and identified potentially implicated liver cell-type for viral ingress. We report the co-expression of ACE2 and TMPRSS2 in a TROP2+ liver progenitor population. Importantly, we fail to detect the expression of ACE2 in hepatocyte or any other liver (immune and stromal) cell types. These results indicated that in COVID-19 associated liver dysfunction and cell death, viral infection of TROP2+ progenitors in liver may significantly impaired liver regeneration and could lead to pathology.Highlights- EPCAM+ Liver progenitors co-express ACE2 and TMPRSS2- ACE2 and TMPRSS2 expression is highest in TROP2high progenitors- ACE2 and TMPRSS2 cells express cholangiocyte biased fate markers- ACE2 and TMPRSS2 positive cells are absent in human fetal liver

Chemo-resistance is one of the major causes of cancer-related deaths. Here we used single-cell transcriptomics to investigate divergent modes of chemo-resistance in tumor cells. We observed that higher degree of phenotypic intra-tumor heterogeneity (ITH) favors selection of pre-existing drug-resistant cells, whereas phenotypically homogeneous cells engage covert epigenetic mechanisms to trans-differentiate under drug-selection. This adaptation was driven by selection-induced gain of H3K27ac marks on bivalently poised resistance-associated chromatin, and therefore not expressed in the treatment-naïve setting. Mechanistic interrogation of this phenomenon revealed that drug-induced adaptation was acquired upon the loss of stem factor SOX2, and a concomitant gain of SOX9. Strikingly we observed an enrichment of SOX9 at drug-induced H3K27ac sites, suggesting that tumor evolution could be driven by stem cell-switch-mediated epigenetic plasticity. Importantly, JQ1 mediated inhibition of BRD4 could reverse drug-induced adaptation. These results provide mechanistic insights into the modes of therapy-induced cellular plasticity and underscore the use of epigenetic inhibitors in targeting tumor evolution.

Summary Here we report a molecular docking-based approach to identify small molecules that can target the β-catenin (β-cat)-TCF4 protein-protein interaction (PPI), a key effector complex for nuclear Wnt signaling activity. Specifically, we developed and optimized a computational model of β-cat using publicly available β-cat protein crystal structures, and existing β-cat-TCF4 interaction inhibitors as the training set. Using our computational model to an in silico screen predicted 27 compounds as good binders to β-cat, of which 3 were identified to be effective against a Wnt-responsive luciferase reporter. In vitro functional validation experiments revealed GB1874 as an inhibitor of the Wnt pathway that targets the β-cat-TCF4 PPI. GB1874 also affected the proliferation and stemness of Wnt-addicted colorectal cancer (CRC) cells in vitro. Encouragingly, GB1874 inhibited the growth of CRC tumor xenografts in vivo, thus demonstrating its potential for further development into therapeutics against Wnt-associated cancer indications., Graphical abstract, Highlights • A computational docking model for β-catenin was developed and optimized. • The model was trained using protein crystal structures and known inhibitor compounds. • Compound GB1874 was identified by in silico screening and functionally validated. • GB1874 inhibited colon cancer growth in vitro and in vivo by blocking Wnt activity., Pharmacology; Biochemistry; Structural biology; Cancer; Omics

Genomics-driven cancer therapeutics has gained prominence in personalized cancer treatment. However, its utility in indications lacking biomarker-driven treatment strategies remains limited. Here we present a “phenotype-driven precision-oncology” approach, based on the notion that biological response to perturbations, chemical or genetic, in ex vivo patient-individualized models can serve as predictive biomarkers for therapeutic response in the clinic. We generated a library of “screenable” patient-derived primary cultures (PDCs) for head and neck squamous cell carcinomas that reproducibly predicted treatment response in matched patient-derived-xenograft models. Importantly, PDCs could guide clinical practice and predict tumour progression in two n = 1 co-clinical trials. Comprehensive “-omics” interrogation of PDCs derived from one of these models revealed YAP1 as a putative biomarker for treatment response and survival in ~24% of oral squamous cell carcinoma. We envision that scaling of the proposed PDC approach could uncover biomarkers for therapeutic stratification and guide real-time therapeutic decisions in the future., Treatment response in patient-derived models may serve as a biomarker for response in the clinic. Here, the authors use paired patient-derived mouse xenografts and patient-derived primary culture models from head and neck squamous cell carcinomas, including metastasis, as models for high-throughput screening of anti-cancer drugs.

Signaling between cancer and nonmalignant (stromal) cells in the tumor microenvironment (TME) is key to tumorigenesis yet challenging to decipher from tumor transcriptomes. Here, we report an unbiased, data-driven approach to deconvolute bulk tumor transcriptomes and predict crosstalk between ligands and receptors on cancer and stromal cells in the TME of 20 solid tumor types. Our approach recovers known transcriptional hallmarks of cancer and stromal cells and is concordant with single-cell and immunohistochemistry data, underlining its robustness. Pan-cancer analysis reveals previously unrecognized features of cancer-stromal crosstalk. We find that autocrine cancer cell cross-talk varied between tissues but often converged on known cancer signaling pathways. In contrast, many stromal cross-talk interactions were highly conserved across tumor types. Interestingly, the immune checkpoint ligand PD-L1 was overexpressed in stromal rather than cancer cells across all tumor types. Moreover, we predicted and experimentally validated aberrant ligand and receptor expression in cancer cells of basal and luminal breast cancer, respectively. Collectively, our findings validate a data-driven method for tumor transcriptome deconvolution and establishes a new resource for hypothesis generation and downstream functional interrogation of the TME in tumorigenesis and disease progression.

Transfer RNA (tRNA)‐derived small RNAs (tsRNAs) have recently emerged as important regulators of protein translation and shown to have diverse biological functions. However, the underlying cellular and molecular mechanisms of tsRNA function in the context of dynamic cell‐state transitions remain unclear. Expression analysis of tsRNAs in distinct heterologous cell and tissue models of stem vs. differentiated states revealed a differentiation‐dependent enrichment of 5′‐tsRNAs. We report the identification of a set of 5′‐tsRNAs that is upregulated in differentiating mouse embryonic stem cells (mESCs). Notably, interactome studies with differentially enriched 5′‐tsRNAs revealed a switch in their association with “effector” RNPs and “target” mRNAs in different cell states. We demonstrate that specific 5′‐tsRNAs can preferentially interact with the RNA‐binding protein, Igf2bp1, in the RA‐induced differentiated state. This association influences the transcript stability and thereby translation of the pluripotency‐promoting factor, c‐Myc, thus providing a mechanistic basis for how 5′‐tsRNAs can modulate stem cell states in mESCs. Together our study highlights the role of 5′‐tsRNAs in defining distinct cell states.