Pancreatic ductal adenocarcinoma (PDAC) tumors can originate either from acinar or ductal cells in the adult pancreas. We re-analyze multiple pancreas and PDAC single-cell RNA-seq datasets and find a subset of nonmalignant acinar cells, which we refer to as acinar edge (AE) cells, whose transcriptomes highly diverge from a typical acinar cell in each dataset. Genes upregulated among AE cells are enriched for transcriptomic signatures of pancreatic progenitors, acinar dedifferentiation, and several oncogenic programs. AE-upregulated genes are upregulated in human PDAC tumors, and consistently, their promoters are hypomethylated. High expression of these genes is associated with poor patient survival. The fraction of AE-like cells increases with age in healthy pancreatic tissue, which is not explained by clonal mutations, thus pointing to a nongenetic source of variation. The fraction of AE-like cells is also significantly higher in human pancreatitis samples. Finally, we find edge-like states in lung, liver, prostate, and colon tissues, suggesting that subpopulations of healthy cells across tissues can exist in pre-neoplastic states. Significance: These findings show “edge” epithelial cell states with oncogenic transcriptional activity in human organs without oncogenic mutations. In the pancreas, the fraction of acinar cells increases with age.

Recent studies exploring the impact of methylation in tumor evolution suggest that while the methylation status of many of the CpG sites are preserved across distinct lineages, others are altered as the cancer progresses. Since changes in methylation status of a CpG site may be retained in mitosis, they could be used to infer the progression history of a tumor via single-cell lineage tree reconstruction. In this work, we introduce the first principled distance-based computational method, Sgootr, for inferring a tumor's single methylation lineage tree and jointly identifying lineage-informative CpG sites which harbor changes in methylation status that are retained along the lineage. We apply Sgootr on single-cell bisulfite-treated whole genome sequencing data of multi-regionally-sampled tumor cells from 9 metastatic colorectal cancer patients, as well as multi-regionally-sampled single-cell reduced-representation bisulfite sequencing data from a glioblastoma patient. We demonstrate that the tumor lineages constructed reveal a simple model underlying tumor progression and metastatic seeding. A comparison of Sgootr against alternative approaches shows that Sgootr can construct lineage trees with fewer migration events and more in concordance with the sequential-progression model of tumor evolution, with a running time a fraction of that used in prior studies. Lineage-informative CpG sites identified by Sgootr are in inter-CpG island (CGI) regions, as opposed to intra-CGIs, which have been the main regions of interest in genomic methylation-related analyses.

Melanoma plasticity, driven by phenotype state switching, underlies clinically relevant traits such as metastasis and therapy resistance. As melanoma progression is thought to recapitulate aspects of neural crest cell (NCC) development, understanding embryonic melanocyte specification and lineage fate decisions of closely related NCCs may illuminate the pathways co-opted during disease evolution. Here, we use a mouse model to isolate and sequence Dopachrome tautomerase (Dct) expressing NCCs, the precursors of melanocytes, at two key developmental stages. We classify these lineages and devise a Developmental Gene Module (DGM) scoring system to interrogate lineage state switching in melanoma samples. In bulk transcriptomes, activation of DGMs representing embryonic Schwann Cell Precursors (SCPs)—multipotent stem cells—in patient tumors predicts poor response to immune checkpoint inhibitors (ICI). Co-activation of SCP and Mesenchymal-like (Mes.) modules further correlates with resistance to MAPK inhibitors. Notably, single-cell analyses reveal that melanoma cells can simultaneously express multiple DGMs, forming “hybrid” states. Cells in a hybrid Neural/SCP state are enriched in early metastasis and ICI-resistant tumors and are insensitive to inflammatory stimuli. We demonstrate that targeting Hdac2 , a histone deacetylase associated with this Neural/SCP hybrid state, promotes a mesenchymal-like state switch, remodels the tumor microenvironment, and sensitizes melanoma cells to TNFα and tumors to ICI therapy. Our methodology thus reveals dynamic patterns of lineage state switching correlated with melanoma tumor evolution to drive insight into new therapeutic targets. Teaser: Newly identified melanoblast cell states are reawakened in metastasizing and therapy-resistant melanomas.

Intratumoral heterogeneity (ITH) can promote cancer progression and treatment failure, but the complexity of the regulatory programs and contextual factors involved complicates its study. To understand the specific contribution of ITH to immune checkpoint blockade (ICB) response, we generated single cell-derived clonal sublines from an ICB-sensitive and genetically and phenotypically heterogeneous mouse melanoma model, M4. Genomic and single cell transcriptomic analyses uncovered the diversity of the sublines and evidenced their plasticity. Moreover, a wide range of tumor growth kinetics were observed in vivo , in part associated with mutational profiles and dependent on T cell-response. Further inquiry into melanoma differentiation states and tumor microenvironment (TME) subtypes of untreated tumors from the clonal sublines demonstrated correlations between highly inflamed and differentiated phenotypes with the response to anti-CTLA-4 treatment. Our results demonstrate that M4 sublines generate intratumoral heterogeneity at both levels of intrinsic differentiation status and extrinsic TME profiles, thereby impacting tumor evolution during therapeutic treatment. These clonal sublines proved to be a valuable resource to study the complex determinants of response to ICB, and specifically the role of melanoma plasticity in immune evasion mechanisms.

Background Exosomal microRNAs (exomiRs), transported via exosomes, play a pivotal role in intercellular communication. In cancer, exomiRs influence tumor progression by regulating key cellular processes such as proliferation, angiogenesis, and metastasis. Their role in mediating communication between cancer cells and the tumor microenvironment highlights their significance as potential diagnostic and therapeutic targets. Methodology In this study, we aimed to characterize the role of exomiRs in influencing the pre-metastatic niche (PMN). Across 7 tumor types, including 4 cell lines and three tumors, we extracted high confidence exomiRs (Log FC >= 2 in exosomes relative to control) and their targets (experimentally identified and targeted by at least 2 exomiRs). Subsequently, we identified enriched pathways and selected the top 100 high-confidence exomiR targets based on the frequency of their appearance in the enriched pathways. These top 100 targets were consistently used throughout the analysis. Results Cancer cell line and tumor derived ExomiRs have significantly higher GC content relative to genomic background. Pathway enriched among the top exomiR targets included general cancer-associated processes such as “wound healing” and “regulation of epithelial cell proliferation”, as well as cancer-specific processes, such as “regulation of angiogenesis in kidney” (KIRC), “ossification” in lung (LUAD), and “positive regulation of cytokine production” in pancreatic cancer (PAAD). Similarly, ‘Pathways in cancer’ and ‘MicroRNAs in cancer’ ranked among the top 10 enriched KEGG pathways in all cancer types. ExomiR targets were not only enriched for cancer-specific tumor suppressor genes (TSG) but are also downregulated in pre-metastatic niche formed in lungs compared to normal lung. Motif analysis shows high similarity among motifs identified from exomiRs across cancer types. Our analysis recapitulates exomiRs associated with M2 macrophage differentiation and chemoresistance such as miR-21 and miR-222-3p, regulating signaling pathways such as PTEN/PI3/Akt, NF-κB, etc. Cox regression indicated that exomiR targets are significantly associated with overall survival of patients in TCGA. Lastly, a Support Vector Machine (SVM) model using exomiR target gene expression classified responders and non-responders to neoadjuvant chemotherapy with an AUROC of 0.96 (in LUAD), higher than other previously reported gene signatures. Conclusion Our study characterizes the pivotal role of exomiRs in shaping the PMN in diverse cancers, underscoring their diagnostic and therapeutic potential.

Aim: Growing evidence points to non-genetic mechanisms underlying long-term resistance to cancer therapies. These mechanisms involve pre-existing or therapy-induced transcriptional cell states that confer resistance. However, the relationship between early transcriptional responses to treatment and the eventual emergence of resistant states remains poorly understood. Furthermore, it is unclear whether such early resistance-associated transcriptional responses are evolutionarily conserved. In this study, we examine the similarity between early transcriptional responses and long-term resistant states, assess their clinical relevance, and explore their evolutionary conservation across species. Methods: We integrated datasets on early drug responses and long-term resistance from multiple cancer cell lines, bacteria, and yeast to identify early transcriptional changes predictive of long-term resistance and assess their evolutionary conservation. Using genome-wide CRISPR-Cas9 knockout screens, we evaluated the impact of genes associated with resistant transcriptional states on drug sensitivity. Clinical datasets were analyzed to explore the prognostic value of the identified resistance-associated gene signatures. Results: We found that transcriptional states observed in drug-naive cells and shortly after treatment overlapped with those seen in fully resistant populations. Some of these shared features appear to be evolutionarily conserved. Knockout of genes marking resistant states sensitized ovarian cancer cells to Prexasertib. Moreover, early resistance gene signatures effectively distinguished therapy responders from non-responders in multiple clinical cancer trials and differentiated premalignant breast lesions that progressed to malignancy from those that remained benign. Conclusion: Early cellular transcriptional responses to therapy exhibit key similarities to fully resistant states across different drugs, cancer types, and species. Gene signatures defining these early resistance states have prognostic value in clinical settings.