A proteogenomic portrait of lung squamous cell carcinoma.

Lung squamous cell carcinoma (LSCC) remains a leading cause of cancer death with few therapeutic options. We characterized the proteogenomic landscape of LSCC, providing a deeper exposition of LSCC biology with potential therapeutic implications. We identify NSD3 as an alternative driver in FGFR1-amplified tumors and low-p63 tumors overexpressing the therapeutic target survivin. SOX2 is considered undruggable, but our analyses provide rationale for exploring chromatin modifiers such as LSD1 and EZH2 to target SOX2-overexpressing tumors. Our data support complex regulation of metabolic pathways by crosstalk between post-translational modifications including ubiquitylation. Numerous immune-related proteogenomic observations suggest directions for further investigation. Proteogenomic dissection of CDKN2A mutations argue for more nuanced assessment of RB1 protein expression and phosphorylation before declaring CDK4/6 inhibition unsuccessful. Finally, triangulation between LSCC, LUAD, and HNSCC identified both unique and common therapeutic vulnerabilities. These observations and proteogenomics data resources may guide research into the biology and treatment of LSCC.

An early cell shape transition drives evolutionary expansion of the human forebrain.

The human brain has undergone rapid expansion since humans diverged from other great apes, but the mechanism of this human-specific enlargement is still unknown. Here, we use cerebral organoids derived from human, gorilla, and chimpanzee cells to study developmental mechanisms driving evolutionary brain expansion. We find that neuroepithelial differentiation is a protracted process in apes, involving a previously unrecognized transition state characterized by a change in cell shape. Furthermore, we show that human organoids are larger due to a delay in this transition, associated with differences in interkinetic nuclear migration and cell cycle length. Comparative RNA sequencing (RNA-seq) reveals differences in expression dynamics of cell morphogenesis factors, including ZEB2, a known epithelial-mesenchymal transition regulator. We show that ZEB2 promotes neuroepithelial transition, and its manipulation and downstream signaling leads to acquisition of nonhuman ape architecture in the human context and vice versa, establishing an important role for neuroepithelial cell shape in human brain expansion.

Regulation of Collective Metastasis by Nanolumenal Signaling.

Collective metastasis is defined as the cohesive migration and metastasis of multicellular tumor cell clusters. Disrupting various cell adhesion genes markedly reduces cluster formation and colonization efficiency, yet the downstream signals transmitted by clustering remain largely unknown. Here, we use mouse and human breast cancer models to identify a collective signal generated by tumor cell clusters supporting metastatic colonization. We show that tumor cell clusters produce the growth factor epigen and concentrate it within nanolumina-intercellular compartments sealed by cell-cell junctions and lined with microvilli-like protrusions. Epigen knockdown profoundly reduces metastatic outgrowth and switches clusters from a proliferative to a collective migratory state. Tumor cell clusters from basal-like 2, but not mesenchymal-like, triple-negative breast cancer cell lines have increased epigen expression, sealed nanolumina, and impaired outgrowth upon nanolumenal junction disruption. We propose that nanolumenal signaling could offer a therapeutic target for aggressive metastatic breast cancers.

Proteogenomic Characterization of Endometrial Carcinoma.

We undertook a comprehensive proteogenomic characterization of 95 prospectively collected endometrial carcinomas, comprising 83 endometrioid and 12 serous tumors. This analysis revealed possible new consequences of perturbations to the p53 and Wnt/ß-catenin pathways, identified a potential role for circRNAs in the epithelial-mesenchymal transition, and provided new information about proteomic markers of clinical and genomic tumor subgroups, including relationships to known druggable pathways. An extensive genome-wide acetylation survey yielded insights into regulatory mechanisms linking Wnt signaling and histone acetylation. We also characterized aspects of the tumor immune landscape, including immunogenic alterations, neoantigens, common cancer/testis antigens, and the immune microenvironment, all of which can inform immunotherapy decisions. Collectively, our multi-omic analyses provide a valuable resource for researchers and clinicians, identify new molecular associations of potential mechanistic significance in the development of endometrial cancers, and suggest novel approaches for identifying potential therapeutic targets.

Nuclear PKM2 binds pre-mRNA at folded G-quadruplexes and reveals their gene regulatory role.

Nuclear localization of the metabolic enzyme PKM2 is widely observed in various cancer types. We identify nuclear PKM2 as a non-canonical RNA-binding protein (RBP) that specifically interacts with folded RNA G-quadruplex (rG4) structures in precursor mRNAs (pre-mRNAs). PKM2 occupancy at rG4s prevents the binding of repressive RBPs, such as HNRNPF, and promotes the expression of rG4-containing pre-mRNAs (the "rG4ome"). We observe an upregulation of the rG4ome during epithelial-to-mesenchymal transition and a negative correlation of rG4 abundance with patient survival in different cancer types. By preventing the nuclear accumulation of PKM2, we could repress the rG4ome in triple-negative breast cancer cells and reduce migration and invasion of cancer cells in vitro and in xenograft mouse models. Our data suggest that the balance of folded and unfolded rG4s controlled by RBPs impacts gene expression during tumor progression.

Cellular adaptation to cancer therapy along a resistance continuum.

Advancements in precision oncology over the past decades have led to new therapeutic interventions, but the efficacy of such treatments is generally limited by an adaptive process that fosters drug resistance1. In addition to genetic mutations2, recent research has identified a role for non-genetic plasticity in transient drug tolerance3 and the acquisition of stable resistance4,5. However, the dynamics of cell-state transitions that occur in the adaptation to cancer therapies remain unknown and require a systems-level longitudinal framework. Here we demonstrate that resistance develops through trajectories of cell-state transitions accompanied by a progressive increase in cell fitness, which we denote as the 'resistance continuum'. This cellular adaptation involves a stepwise assembly of gene expression programmes and epigenetically reinforced cell states underpinned by phenotypic plasticity, adaptation to stress and metabolic reprogramming. Our results support the notion that epithelial-to-mesenchymal transition or stemness programmes-often considered a proxy for phenotypic plasticity-enable adaptation, rather than a full resistance mechanism. Through systematic genetic perturbations, we identify the acquisition of metabolic dependencies, exposing vulnerabilities that can potentially be exploited therapeutically. The concept of the resistance continuum highlights the dynamic nature of cellular adaptation and calls for complementary therapies directed at the mechanisms underlying adaptive cell-state transitions.

In vivo single-cell CRISPR uncovers distinct TNF programmes in tumour evolution.

The tumour evolution model posits that malignant transformation is preceded by randomly distributed driver mutations in cancer genes, which cause clonal expansions in phenotypically normal tissues. Although clonal expansions can remodel entire tissues1-3, the mechanisms that result in only a small number of clones transforming into malignant tumours remain unknown. Here we develop an in vivo single-cell CRISPR strategy to systematically investigate tissue-wide clonal dynamics of the 150 most frequently mutated squamous cell carcinoma genes. We couple ultrasound-guided in utero lentiviral microinjections, single-cell RNA sequencing and guide capture to longitudinally monitor clonal expansions and document their underlying gene programmes at single-cell transcriptomic resolution. We uncover a tumour necrosis factor (TNF) signalling module, which is dependent on TNF receptor 1 and involving macrophages, that acts as a generalizable driver of clonal expansions in epithelial tissues. Conversely, during tumorigenesis, the TNF signalling module is downregulated. Instead, we identify a subpopulation of invasive cancer cells that switch to an autocrine TNF gene programme associated with epithelial-mesenchymal transition. Finally, we provide in vivo evidence that the autocrine TNF gene programme is sufficient to mediate invasive properties and show that the TNF signature correlates with shorter overall survival of patients with squamous cell carcinoma. Collectively, our study demonstrates the power of applying in vivo single-cell CRISPR screening to mammalian tissues, unveils distinct TNF programmes in tumour evolution and highlights the importance of understanding the relationship between clonal expansions in epithelia and tumorigenesis.

The epithelial-mesenchymal plasticity landscape: principles of design and mechanisms of regulation.

Epithelial-mesenchymal plasticity (EMP) enables cells to interconvert between several states across the epithelial-mesenchymal landscape, thereby acquiring hybrid epithelial/mesenchymal phenotypic features. This plasticity is crucial for embryonic development and wound healing, but also underlies the acquisition of several malignant traits during cancer progression. Recent research using systems biology and single-cell profiling methods has provided novel insights into the main forces that shape EMP, which include the microenvironment, lineage specification and cell identity, and the genome. Additionally, key roles have emerged for hysteresis (cell memory) and cellular noise, which can drive stochastic transitions between cell states. Here, we review these forces and the distinct but interwoven layers of regulatory control that stabilize EMP states or facilitate epithelial-mesenchymal transitions (EMTs) and discuss the therapeutic potential of manipulating the EMP landscape.

Inhibition of YTHDF2 triggers proteotoxic cell death in MYC-driven breast cancer.

RNA-binding proteins (RBPs) are critical regulators of post-transcriptional gene expression, and aberrant RBP-RNA interactions can promote cancer progression. Here, we interrogate the function of RBPs in cancer using pooled CRISPR-Cas9 screening and identify 57 RBP candidates with distinct roles in supporting MYC-driven oncogenic pathways. We find that disrupting YTHDF2-dependent mRNA degradation triggers apoptosis in triple-negative breast cancer (TNBC) cells and tumors. eCLIP and m6A sequencing reveal that YTHDF2 interacts with mRNAs encoding proteins in the MAPK pathway that, when stabilized, induce epithelial-to-mesenchymal transition and increase global translation rates. scRibo-STAMP profiling of translating mRNAs reveals unique alterations in the translatome of single cells within YTHDF2-depleted solid tumors, which selectively contribute to endoplasmic reticulum stress-induced apoptosis in TNBC cells. Thus, our work highlights the therapeutic potential of RBPs by uncovering a critical role for YTHDF2 in counteracting the global increase of mRNA synthesis in MYC-driven breast cancers.

GREM1 is required to maintain cellular heterogeneity in pancreatic cancer.

Pancreatic ductal adenocarcinoma (PDAC) shows pronounced epithelial and mesenchymal cancer cell populations1-4. Cellular heterogeneity in PDAC is an important feature in disease subtype specification3-5, but how distinct PDAC subpopulations interact, and the molecular mechanisms that underlie PDAC cell fate decisions, are incompletely understood. Here we identify the BMP inhibitor GREM16,7 as a key regulator of cellular heterogeneity in pancreatic cancer in human and mouse. Grem1 inactivation in established PDAC in mice resulted in a direct conversion of epithelial into mesenchymal PDAC cells within days, suggesting that persistent GREM1 activity is required to maintain the epithelial PDAC subpopulations. By contrast, Grem1 overexpression caused an almost complete 'epithelialization' of highly mesenchymal PDAC, indicating that high GREM1 activity is sufficient to revert the mesenchymal fate of PDAC cells. Mechanistically, Grem1 was highly expressed in mesenchymal PDAC cells and inhibited the expression of the epithelial-mesenchymal transition transcription factors Snai1 (also known as Snail) and Snai2 (also known as Slug) in the epithelial cell compartment, therefore restricting epithelial-mesenchymal plasticity. Thus, constant suppression of BMP activity is essential to maintain epithelial PDAC cells, indicating that the maintenance of the cellular heterogeneity of pancreatic cancer requires continuous paracrine signalling elicited by a single soluble factor.

ST3GAL5-catalyzed gangliosides inhibit TGF-ß-induced epithelial-mesenchymal transition via TßRI degradation.

Epithelial-mesenchymal transition (EMT) is pivotal in the initiation and development of cancer cell metastasis. We observed that the abundance of glycosphingolipids (GSLs), especially ganglioside subtypes, decreased significantly during TGF-ß-induced EMT in NMuMG mouse mammary epithelial cells and A549 human lung adenocarcinoma cells. Transcriptional profiling showed that TGF-ß/SMAD response genes and EMT signatures were strongly enriched in NMuMG cells, along with depletion of UDP-glucose ceramide glucosyltransferase (UGCG), the enzyme that catalyzes the initial step in GSL biosynthesis. Consistent with this finding, genetic or pharmacological inhibition of UGCG promoted TGF-ß signaling and TGF-ß-induced EMT. UGCG inhibition promoted A549 cell migration, extravasation in the zebrafish xenograft model, and metastasis in mice. Mechanistically, GSLs inhibited TGF-ß signaling by promoting lipid raft localization of the TGF-ß type I receptor (TßRI) and by increasing TßRI ubiquitination and degradation. Importantly, we identified ST3GAL5-synthesized a-series gangliosides as the main GSL subtype involved in inhibition of TGF-ß signaling and TGF-ß-induced EMT in A549 cells. Notably, ST3GAL5 is weakly expressed in lung cancer tissues compared to adjacent nonmalignant tissues, and its expression correlates with good prognosis.

LncRNA LITATS1 suppresses TGF-ß-induced EMT and cancer cell plasticity by potentiating TßRI degradation.

Epithelial cells acquire mesenchymal phenotypes through epithelial-mesenchymal transition (EMT) during cancer progression. However, how epithelial cells retain their epithelial traits and prevent malignant transformation is not well understood. Here, we report that the long noncoding RNA LITATS1 (LINC01137, ZC3H12A-DT) is an epithelial gatekeeper in normal epithelial cells and inhibits EMT in breast and non-small cell lung cancer cells. Transcriptome analysis identified LITATS1 as a TGF-ß target gene. LITATS1 expression is reduced in lung adenocarcinoma tissues compared with adjacent normal tissues and correlates with a favorable prognosis in breast and non-small cell lung cancer patients. LITATS1 depletion promotes TGF-ß-induced EMT, migration, and extravasation in cancer cells. Unbiased pathway analysis demonstrated that LITATS1 knockdown potently and selectively potentiates TGF-ß/SMAD signaling. Mechanistically, LITATS1 enhances the polyubiquitination and proteasomal degradation of TGF-ß type I receptor (TßRI). LITATS1 interacts with TßRI and the E3 ligase SMURF2, promoting the cytoplasmic retention of SMURF2. Our findings highlight a protective function of LITATS1 in epithelial integrity maintenance through the attenuation of TGF-ß/SMAD signaling and EMT.

ARID1A-BAF coordinates ZIC2 genomic occupancy for epithelial-to-mesenchymal transition in cranial neural crest specification.

The BAF chromatin remodeler regulates lineage commitment including cranial neural crest cell (CNCC) specification. Variants in BAF subunits cause Coffin-Siris syndrome (CSS), a congenital disorder characterized by coarse craniofacial features and intellectual disability. Approximately 50% of individuals with CSS harbor variants in one of the mutually exclusive BAF subunits, ARID1A/ARID1B. While Arid1a deletion in mouse neural crest causes severe craniofacial phenotypes, little is known about the role of ARID1A in CNCC specification. Using CSS-patient-derived ARID1A+/- induced pluripotent stem cells to model CNCC specification, we discovered that ARID1A-haploinsufficiency impairs epithelial-to-mesenchymal transition (EMT), a process necessary for CNCC delamination and migration from the neural tube. Furthermore, wild-type ARID1A-BAF regulates enhancers associated with EMT genes. ARID1A-BAF binding at these enhancers is impaired in heterozygotes while binding at promoters is unaffected. At the sequence level, these EMT enhancers contain binding motifs for ZIC2, and ZIC2 binding at these sites is ARID1A-dependent. When excluded from EMT enhancers, ZIC2 relocates to neuronal enhancers, triggering aberrant neuronal gene activation. In mice, deletion of Zic2 impairs NCC delamination, while ZIC2 overexpression in chick embryos at post-migratory neural crest stages elicits ectopic delamination from the neural tube. These findings reveal an essential ARID1A-ZIC2 axis essential for EMT and CNCC delamination.

Machine learning identifies activation of RUNX/AP-1 as drivers of mesenchymal and fibrotic regulatory programs in gastric cancer.

Gastric cancer (GC) is the fifth most common cancer worldwide and is a heterogeneous disease. Among GC subtypes, the mesenchymal phenotype (Mes-like) is more invasive than the epithelial phenotype (Epi-like). Although gene expression of the epithelial-to-mesenchymal transition (EMT) has been studied, the regulatory landscape shaping this process is not fully understood. Here we use ATAC-seq and RNA-seq data from a compendium of GC cell lines and primary tumors to detect drivers of regulatory state changes and their transcriptional responses. Using the ATAC-seq data, we developed a machine learning approach to determine the transcription factors (TFs) regulating the subtypes of GC. We identified TFs driving the mesenchymal (RUNX2, ZEB1, SNAI2, AP-1 dimer) and the epithelial (GATA4, GATA6, KLF5, HNF4A, FOXA2, GRHL2) states in GC. We identified DNA copy number alterations associated with dysregulation of these TFs, specifically deletion of GATA4 and amplification of MAPK9 Comparisons with bulk and single-cell RNA-seq data sets identified activation toward fibroblast-like epigenomic and expression signatures in Mes-like GC. The activation of this mesenchymal fibrotic program is associated with differentially accessible DNA cis-regulatory elements flanking upregulated mesenchymal genes. These findings establish a map of TF activity in GC and highlight the role of copy number driven alterations in shaping epigenomic regulatory programs as potential drivers of GC heterogeneity and progression.

Cadherin adhesion complexes direct cell aggregation in the epithelial transition of Wnt-induced nephron progenitor cells.

In the developing mammalian kidney, nephron formation is initiated by a subset of nephron progenitor cells (NPCs). Wnt input activates a ß-catenin (Ctnnb1)-driven, transcriptional nephrogenic program and the mesenchymal to epithelial transition (MET) of NPCs. Using an in vitro mouse NPC culture model, we observed that activation of the Wnt pathway results in the aggregation of induced NPCs, which is an initiating step in the MET program. Genetic removal showed aggregation was dependent on ß-catenin. Modulating extracellular Ca2+ levels showed cell-cell contacts were Ca2+ dependent, suggesting a role for cadherin (Cdh)-directed cell adhesion. Molecular analysis identified Cdh2, Cdh4 and Cdh11 in NPCs, and the ß-catenin directed upregulation of Cdh3 and Cdh4 accompanying the MET of induced NPCs. Mutational analysis of ß-catenin supported a role for a Lef/Tcf-ß-catenin-mediated transcriptional response in the cell aggregation process. Genetic removal of all four cadherins, and independent removal of α-catenin or of ß-catenin-α-catenin interactions, abolished aggregation, but not the inductive response to Wnt pathway activation. These findings, and data in an accompanying article highlight the role of ß-catenin in linking transcriptional programs to the morphogenesis of NPCs in mammalian nephrogenesis.

Spatial transition tensor of single cells.

Spatial transcriptomics and messenger RNA splicing encode extensive spatiotemporal information for cell states and transitions. The current lineage-inference methods either lack spatial dynamics for state transition or cannot capture different dynamics associated with multiple cell states and transition paths. Here we present spatial transition tensor (STT), a method that uses messenger RNA splicing and spatial transcriptomes through a multiscale dynamical model to characterize multistability in space. By learning a four-dimensional transition tensor and spatial-constrained random walk, STT reconstructs cell-state-specific dynamics and spatial state transitions via both short-time local tensor streamlines between cells and long-time transition paths among attractors. Benchmarking and applications of STT on several transcriptome datasets via multiple technologies on epithelial-mesenchymal transitions, blood development, spatially resolved mouse brain and chicken heart development, indicate STT's capability in recovering cell-state-specific dynamics and their associated genes not seen using existing methods. Overall, STT provides a consistent multiscale description of single-cell transcriptome data across multiple spatiotemporal scales.

SAMD1 suppresses epithelial-mesenchymal transition pathways in pancreatic ductal adenocarcinoma.

Pancreatic ductal adenocarcinoma (PDAC) poses a significant threat due to its tendency to evade early detection, frequent metastasis, and the subsequent challenges in devising effective treatments. Processes that govern epithelial-mesenchymal transition (EMT) in PDAC hold promise for advancing novel therapeutic strategies. SAMD1 (SAM domain-containing protein 1) is a CpG island-binding protein that plays a pivotal role in the repression of its target genes. Here, we revealed that SAMD1 acts as a repressor of genes associated with EMT. Upon deletion of SAMD1 in PDAC cells, we observed significantly increased migration rates. SAMD1 exerts its effects by binding to specific genomic targets, including CDH2, encoding N-cadherin, which emerged as a driver of enhanced migration upon SAMD1 knockout. Furthermore, we discovered the FBXO11-containing E3 ubiquitin ligase complex as an interactor and negative regulator of SAMD1, which inhibits SAMD1 chromatin-binding genome-wide. High FBXO11 expression in PDAC is associated with poor prognosis and increased expression of EMT-related genes, underlining an antagonistic relationship between SAMD1 and FBXO11. In summary, our findings provide insights into the regulation of EMT-related genes in PDAC, shedding light on the intricate role of SAMD1 and its interplay with FBXO11 in this cancer type.

Fat1 deletion promotes hybrid EMT state, tumour stemness and metastasis.

FAT1, which encodes a protocadherin, is one of the most frequently mutated genes in human cancers1-5. However, the role and the molecular mechanisms by which FAT1 mutations control tumour initiation and progression are poorly understood. Here, using mouse models of skin squamous cell carcinoma and lung tumours, we found that deletion of Fat1 accelerates tumour initiation and malignant progression and promotes a hybrid epithelial-to-mesenchymal transition (EMT) phenotype. We also found this hybrid EMT state in FAT1-mutated human squamous cell carcinomas. Skin squamous cell carcinomas in which Fat1 was deleted presented increased tumour stemness and spontaneous metastasis. We performed transcriptional and chromatin profiling combined with proteomic analyses and mechanistic studies, which revealed that loss of function of FAT1 activates a CAMK2-CD44-SRC axis that promotes YAP1 nuclear translocation and ZEB1 expression that stimulates the mesenchymal state. This loss of function also inactivates EZH2, promoting SOX2 expression, which sustains the epithelial state. Our comprehensive analysis identified drug resistance and vulnerabilities in FAT1-deficient tumours, which have important implications for cancer therapy. Our studies reveal that, in mouse and human squamous cell carcinoma, loss of function of FAT1 promotes tumour initiation, progression, invasiveness, stemness and metastasis through the induction of a hybrid EMT state.

Reconstruction of single-cell lineage trajectories and identification of diversity in fates during the epithelial-to-mesenchymal transition.

Exploring the complexity of the epithelial-to-mesenchymal transition (EMT) unveils a diversity of potential cell fates; however, the exact timing and mechanisms by which early cell states diverge into distinct EMT trajectories remain unclear. Studying these EMT trajectories through single-cell RNA sequencing is challenging due to the necessity of sacrificing cells for each measurement. In this study, we employed optimal-transport analysis to reconstruct the past trajectories of different cell fates during TGF-beta-induced EMT in the MCF10A cell line. Our analysis revealed three distinct trajectories leading to low EMT, partial EMT, and high EMT states. Cells along the partial EMT trajectory showed substantial variations in the EMT signature and exhibited pronounced stemness. Throughout this EMT trajectory, we observed a consistent downregulation of the EED and EZH2 genes. This finding was validated by recent inhibitor screens of EMT regulators and CRISPR screen studies. Moreover, we applied our analysis of early-phase differential gene expression to gene sets associated with stemness and proliferation, pinpointing ITGB4, LAMA3, and LAMB3 as genes differentially expressed in the initial stages of the partial versus high EMT trajectories. We also found that CENPF, CKS1B, and MKI67 showed significant upregulation in the high EMT trajectory. While the first group of genes aligns with findings from previous studies, our work uniquely pinpoints the precise timing of these upregulations. Finally, the identification of the latter group of genes sheds light on potential cell cycle targets for modulating EMT trajectories.