The TERT promoter is polycomb-repressed in neuroblastoma cells with long telomeres.

Acquiring a telomere maintenance mechanism is a hallmark of high-risk neuroblastoma and commonly occurs by expressing telomerase (TERT). Telomerase-negative neuroblastoma has long telomeres and utilize the telomerase-independent alternative lengthening of telomeres (ALT) mechanism. Conversely, no discernable telomere maintenance mechanism is detected in a fraction of neuroblastoma with long telomeres. Here, we show, unlike most cancers, DNA of the TERT promoter is broadly hypomethylated in neuroblastoma. In telomerase-positive neuroblastoma cells, the hypomethylated DNA promoter is approximately 1.5-kb. The TERT locus shows active chromatin marks with low enrichment for the repressive mark, H3K27me3. MYCN, a commonly amplified oncogene in neuroblstoma, binds to the promoter and induces TERT expression. Strikingly, in neuroblastoma with long telomeres, the hypomethylated region spans the entire TERT locus, including multiple nearby genes with enrichment for the repressive H3K27me3 chromatin mark. Furthermore, subtelomeric regions showed enrichment of repressive chromatin marks in neuroblastomas with long telomeres relative to those with short telomeres. These repressive marks were even more evident at the genic loci, suggesting a telomere position effect. Inhibiting H3K27 methylation by three different EZH2 inhibitors induced the expression of TERT in cell lines with long telomeres and H3K27me3 marks in the promoter region. EZH2 inhibition facilitated MYCN binding to the TERT promoter in neuroblastoma cells with long telomeres. Taken together, these data suggest that epigenetic regulation of TERT expression differs in neuroblastoma depending on the telomere maintenance status, and H3K27 methylation is important in repressing TERT expression in neuroblastoma with long telomeres.

A single-vector intersectional AAV strategy for interrogating cellular diversity and brain function.

As discovery of cellular diversity in the brain accelerates, so does the need for tools that target cells based on multiple features. Here we developed Conditional Viral Expression by Ribozyme Guided Degradation (ConVERGD), an adeno-associated virus-based, single-construct, intersectional targeting strategy that combines a self-cleaving ribozyme with traditional FLEx switches to deliver molecular cargo to specific neuronal subtypes. ConVERGD offers benefits over existing intersectional expression platforms, such as expanded intersectional targeting with up to five recombinase-based features, accommodation of larger and more complex payloads and a vector that is easy to modify for rapid toolkit expansion. In the present report we employed ConVERGD to characterize an unexplored subpopulation of norepinephrine (NE)-producing neurons within the rodent locus coeruleus that co-express the endogenous opioid gene prodynorphin (Pdyn). These studies showcase ConVERGD as a versatile tool for targeting diverse cell types and reveal Pdyn-expressing NE+ locus coeruleus neurons as a small neuronal subpopulation capable of driving anxiogenic behavioral responses in rodents.

Functional characterization of cooperating MGA mutations in RUNX1::RUNX1T1 acute myeloid leukemia.

MGA (Max-gene associated) is a dual-specificity transcription factor that negatively regulates MYC-target genes to inhibit proliferation and promote differentiation. Loss-of-function mutations in MGA have been commonly identified in several hematological neoplasms, including acute myeloid leukemia (AML) with RUNX1::RUNX1T1, however, very little is known about the impact of these MGA alterations on normal hematopoiesis or disease progression. We show that representative MGA mutations identified in patient samples abolish protein-protein interactions and transcriptional activity. Using a series of human and mouse model systems, including a newly developed conditional knock-out mouse strain, we demonstrate that loss of MGA results in upregulation of MYC and E2F targets, cell cycle genes, mTOR signaling, and oxidative phosphorylation in normal hematopoietic cells, leading to enhanced proliferation. The loss of MGA induces an open chromatin state at promoters of genes involved in cell cycle and proliferation. RUNX1::RUNX1T1 expression in Mga-deficient murine hematopoietic cells leads to a more aggressive AML with a significantly shortened latency. These data show that MGA regulates multiple pro-proliferative pathways in hematopoietic cells and cooperates with the RUNX1::RUNX1T1 fusion oncoprotein to enhance leukemogenesis.

Histone lysine demethylase 4 family proteins maintain the transcriptional program and adrenergic cellular state of MYCN-amplified neuroblastoma.

Neuroblastoma with MYCN amplification (MNA) is a high-risk disease that has a poor survival rate. Neuroblastoma displays cellular heterogeneity, including more differentiated (adrenergic) and more primitive (mesenchymal) cellular states. Here, we demonstrate that MYCN oncoprotein promotes a cellular state switch in mesenchymal cells to an adrenergic state, accompanied by induction of histone lysine demethylase 4 family members (KDM4A-C) that act in concert to control the expression of MYCN and adrenergic core regulatory circulatory (CRC) transcription factors. Pharmacologic inhibition of KDM4 blocks expression of MYCN and the adrenergic CRC transcriptome with genome-wide induction of transcriptionally repressive H3K9me3, resulting in potent anticancer activity against neuroblastomas with MNA by inducing neuroblastic differentiation and apoptosis. Furthermore, a short-term KDM4 inhibition in combination with conventional, cytotoxic chemotherapy results in complete tumor responses of xenografts with MNA. Thus, KDM4 blockade may serve as a transformative strategy to target the adrenergic CRC dependencies in MNA neuroblastomas.

THUMPD2 catalyzes the N2-methylation of U6 snRNA of the spliceosome catalytic center and regulates pre-mRNA splicing and retinal degeneration.

The mechanisms by which the relatively conserved spliceosome manages the enormously large number of splicing events that occur in humans (∼200 000 versus ∼300 in yeast) are poorly understood. Here, we show deposition of one RNA modification-N2-methylguanosine (m2G) on the G72 of U6 snRNA (the catalytic center of the spliceosome) promotes efficient pre-mRNA splicing activity in human cells. This modification was identified to be conserved among vertebrates. Further, THUMPD2 was demonstrated as the methyltransferase responsible for U6 m2G72 by explicitly recognizing the U6-specific sequences and structural elements. The knock-out of THUMPD2 eliminated U6 m2G72 and impaired the pre-mRNA splicing activity, resulting in thousands of changed alternative splicing events of endogenous pre-mRNAs in human cells. Notably, the aberrantly spliced pre-mRNA population elicited the nonsense-mediated mRNA decay pathway. We further show that THUMPD2 was associated with age-related macular degeneration and retinal function. Our study thus demonstrates how an RNA epigenetic modification of the major spliceosome regulates global pre-mRNA splicing and impacts physiology and disease.

RNA-binding protein RBM5 plays an essential role in acute myeloid leukemia by activating the oncogenic protein HOXA9.

BACKGROUND: The oncogenic protein HOXA9 plays a critical role in leukemia transformation and maintenance, and its aberrant expression is a hallmark of most aggressive acute leukemia. Although inhibiting the upstream regulators of HOXA9 has been proven as a significant therapeutic intervention, the comprehensive regulation network controlling HOXA9 expression in leukemia has not been systematically investigated. RESULTS: Here, we perform genome-wide CRISPR/Cas9 screening in the HOXA9-driven reporter acute leukemia cells. We identify a poorly characterized RNA-binding protein, RBM5, as the top candidate gene required to maintain leukemia cell fitness. RBM5 is highly overexpressed in acute myeloid leukemia (AML) patients compared to healthy individuals. RBM5 loss triggered by CRISPR knockout and shRNA knockdown significantly impairs leukemia maintenance in vitro and in vivo. Through domain CRISPR screening, we reveal that RBM5 functions through a noncanonical transcriptional regulation circuitry rather than RNA splicing, such an effect depending on DNA-binding domains. By integrative analysis and functional assays, we identify HOXA9 as the downstream target of RBM5. Ectopic expression of HOXA9 rescues impaired leukemia cell proliferation upon RBM5 loss. Importantly, acute protein degradation of RBM5 through auxin-inducible degron system immediately reduces HOXA9 transcription. CONCLUSIONS: We identify RBM5 as a new upstream regulator of HOXA9 and reveal its essential role in controlling the survival of AML. These functional and molecular mechanisms further support RBM5 as a promising therapeutic target for myeloid leukemia treatment.

Acute myeloid leukemias with UBTF tandem duplications are sensitive to menin inhibitors.

ABSTRACT: UBTF tandem duplications (UBTF-TDs) have recently emerged as a recurrent alteration in pediatric and adult acute myeloid leukemia (AML). UBTF-TD leukemias are characterized by a poor response to conventional chemotherapy and a transcriptional signature that mirrors NUP98-rearranged and NPM1-mutant AMLs, including HOX-gene dysregulation. However, the mechanism by which UBTF-TD drives leukemogenesis remains unknown. In this study, we investigated the genomic occupancy of UBTF-TD in transformed cord blood CD34+ cells and patient-derived xenograft models. We found that UBTF-TD protein maintained genomic occupancy at ribosomal DNA loci while also occupying genomic targets commonly dysregulated in UBTF-TD myeloid malignancies, such as the HOXA/HOXB gene clusters and MEIS1. These data suggest that UBTF-TD is a gain-of-function alteration that results in mislocalization to genomic loci dysregulated in UBTF-TD leukemias. UBTF-TD also co-occupies key genomic loci with KMT2A and menin, which are known to be key partners involved in HOX-dysregulated leukemias. Using a protein degradation system, we showed that stemness, proliferation, and transcriptional signatures are dependent on sustained UBTF-TD localization to chromatin. Finally, we demonstrate that primary cells from UBTF-TD leukemias are sensitive to the menin inhibitor SNDX-5613, resulting in markedly reduced in vitro and in vivo tumor growth, myeloid differentiation, and abrogation of the UBTF-TD leukemic expression signature. These findings provide a viable therapeutic strategy for patients with this high-risk AML subtype.

ABCC4 impacts megakaryopoiesis and protects megakaryocytes against 6-mercaptopurine induced cytotoxicity.

The role of ABCC4, an ATP-binding cassette transporter, in the process of platelet formation, megakaryopoiesis, is unknown. Here, we show that ABCC4 is highly expressed in megakaryocytes (MKs). Mining of public genomic data (ATAC-seq and genome wide chromatin interactions, Hi-C) revealed that key megakaryopoiesis transcription factors (TFs) interacted with ABCC4 regulatory elements and likely accounted for high ABCC4 expression in MKs. Importantly these genomic interactions for ABCC4 ranked higher than for genes with known roles in megakaryopoiesis suggesting a role for ABCC4 in megakaryopoiesis. We then demonstrate that ABCC4 is required for optimal platelet formation as in vitro differentiation of fetal liver derived MKs from Abcc4-/- mice exhibited impaired proplatelet formation and polyploidization, features required for optimal megakaryopoiesis. Likewise, a human megakaryoblastic cell line, MEG-01 showed that acute ABCC4 inhibition markedly suppressed key processes in megakaryopoiesis and that these effects were related to reduced cAMP export and enhanced dissociation of a negative regulator of megakaryopoiesis, protein kinase A (PKA) from ABCC4. PKA activity concomitantly increased after ABCC4 inhibition which was coupled with significantly reduced GATA-1 expression, a TF needed for optimal megakaryopoiesis. Further, ABCC4 protected MKs from 6-mercaptopurine (6-MP) as Abcc4-/- mice show a profound reduction in MKs after 6-MP treatment. In total, our studies show that ABCC4 not only protects the MKs but is also required for maximal platelet production from MKs, suggesting modulation of ABCC4 function might be a potential therapeutic strategy to regulate platelet production.

PAX3-FOXO1 dictates myogenic reprogramming and rhabdomyosarcoma identity in endothelial progenitors.

Fusion-positive rhabdomyosarcoma (FP-RMS) driven by the expression of the PAX3-FOXO1 (P3F) fusion oncoprotein is an aggressive subtype of pediatric rhabdomyosarcoma. FP-RMS histologically resembles developing muscle yet occurs throughout the body in areas devoid of skeletal muscle highlighting that FP-RMS is not derived from an exclusively myogenic cell of origin. Here we demonstrate that P3F reprograms mouse and human endothelial progenitors to FP-RMS. We show that P3F expression in aP2-Cre expressing cells reprograms endothelial progenitors to functional myogenic stem cells capable of regenerating injured muscle fibers. Further, we describe a FP-RMS mouse model driven by P3F expression and Cdkn2a loss in endothelial cells. Additionally, we show that P3F expression in TP53-null human iPSCs blocks endothelial-directed differentiation and guides cells to become myogenic cells that form FP-RMS tumors in immunocompromised mice. Together these findings demonstrate that FP-RMS can originate from aberrant development of non-myogenic cells driven by P3F.

Systematic characterization of the HOXA9 downstream targets in MLL-r leukemia by noncoding CRISPR screens.

Accumulating evidence indicates that HOXA9 dysregulation is necessary and sufficient for leukemic transformation and maintenance. However, it remains largely unknown how HOXA9, as a homeobox transcriptional factor, binds to noncoding regulatory sequences and controls the downstream genes. Here, we conduct dropout CRISPR screens against 229 HOXA9-bound peaks identified by ChIP-seq. Integrative data analysis identifies reproducible noncoding hits, including those located in the distal enhancer of FLT3 and intron of CDK6. The Cas9-editing and dCas9-KRAB silencing of the HOXA9-bound sites significantly reduce corresponding gene transcription and impair cell proliferation in vitro, and in vivo by transplantation into NSG female mice. In addition, RNA-seq, Q-PCR analysis, chromatin accessibility change, and chromatin conformation evaluation uncover the noncoding regulation mechanism of HOXA9 and its functional downstream genes. In summary, our work improves our understanding of how HOXA9-associated transcription programs reconstruct the regulatory network specifying MLL-r dependency.

Functional Characterization of Cooperating MGA Mutations in RUNX1::RUNX1T1 Acute Myeloid Leukemia.

MGA (Max-gene associated) is a dual-specificity transcription factor that negatively regulates MYC-target genes to inhibit proliferation and promote differentiation. Loss-of-function mutations in MGA have been commonly identified in several hematological neoplasms, including acute myeloid leukemia (AML) with RUNX1::RUNX1T1, however, very little is known about the impact of these MGA alterations on normal hematopoiesis or disease progression. We show that representative MGA mutations identified in patient samples abolish protein-protein interactions and transcriptional activity. Using a series of human and mouse model systems, including a newly developed conditional knock-out mouse strain, we demonstrate that loss of MGA results in upregulation of MYC and E2F targets, cell cycle genes, mTOR signaling, and oxidative phosphorylation in normal hematopoietic cells, leading to enhanced proliferation. The loss of MGA induces an open chromatin state at promotors of genes involved in cell cycle and proliferation. RUNX1::RUNX1T1 expression in Mga-deficient murine hematopoietic cells leads to a more aggressive AML with a significantly shortened latency. These data show that MGA regulates multiple pro-proliferative pathways in hematopoietic cells and cooperates with the RUNX1::RUNX1 T1 fusion oncoprotein to enhance leukemogenesis.

Latent Epigenetic Programs in Müller Glia Contribute to Stress, Injury, and Disease Response in the Retina.

Previous studies have demonstrated the dynamic changes in chromatin structure during retinal development that correlate with changes in gene expression. However, a major limitation of those prior studies was the lack of cellular resolution. Here, we integrate single-cell (sc) RNA-seq and scATAC-seq with bulk retinal data sets to identify cell type-specific changes in the chromatin structure during development. Although most genes' promoter activity is strongly correlated with chromatin accessibility, we discovered several hundred genes that were transcriptionally silent but had accessible chromatin at their promoters. Most of those silent/accessible gene promoters were in the Müller glial cells. The Müller cells are radial glia of the retina and perform a variety of essential functions to maintain retinal homeostasis and respond to stress, injury, or disease. The silent/accessible genes in Müller glia are enriched in pathways related to inflammation, angiogenesis, and other types of cell-cell signaling and were rapidly activated when we tested 15 different physiologically relevant conditions to mimic retinal stress, injury, or disease in human and murine retinae. We refer to these as "pliancy genes" because they allow the Müller glia to rapidly change their gene expression and cellular state in response to different types of retinal insults. The Müller glial cell pliancy program is established during development, and we demonstrate that pliancy genes are necessary and sufficient for regulating inflammation in the murine retina in vivo. In zebrafish, Müller glia can de-differentiate and form retinal progenitor cells that replace lost neurons. The pro-inflammatory pliancy gene cascade is not activated in zebrafish Müller glia following injury, and we propose a model in which species-specific pliancy programs underly the differential response to retinal damage in species that can regenerate retinal neurons (zebrafish) versus those that cannot (humans and mice).

Identification of Evolutionarily Conserved VSX2 Enhancers in Retinal Development.

Super-enhancers (SEs) are expansive regions of genomic DNA that regulate the expression of genes involved in cell identity and cell fate. Recently, we found that distinct modules within a murine SE regulate gene expression of master regulatory transcription factor Vsx2 in a developmental stage- and cell-type specific manner. Vsx2 is expressed in retinal progenitor cells as well as differentiated bipolar neurons and Müller glia. Mutations in VSX2 in humans and mice lead to microphthalmia due to a defect in retinal progenitor cell proliferation. Deletion of a single module within the Vsx2 SE leads to microphthalmia. Deletion of a separate module within the SE leads to a complete loss of bipolar neurons, yet the remainder of the retina develops normally. Furthermore, the Vsx2 SE is evolutionarily conserved in vertebrates, suggesting that these modules are important for retinal development across species. In the present study, we examine the ability of these modules to drive retinal development between species. By inserting the human build of one Vsx2 SE module into a mouse with microphthalmia, eye size was rescued. To understand the implications of these SE modules in a model of human development, we generated human retinal organoids. Deleting one module results in small organoids, recapitulating the small-eyed phenotype of mice with microphthalmia, while deletion of the other module leads to a complete loss of ON cone bipolar neurons. This prototypical SE serves as a model for uncoupling developmental stage- and cell-type specific effects of neurogenic transcription factors with complex expression patterns. Moreover, by elucidating the gene regulatory mechanisms, we can begin to examine how dysregulation of these mechanisms contributes to phenotypic diversity and disease.

SNIP1 and PRC2 coordinate cell fates of neural progenitors during brain development.

Stem cell survival versus death is a developmentally programmed process essential for morphogenesis, sizing, and quality control of genome integrity and cell fates. Cell death is pervasive during development, but its programming is little known. Here, we report that Smad nuclear interacting protein 1 (SNIP1) promotes neural progenitor cell survival and neurogenesis and is, therefore, integral to brain development. The SNIP1-depleted brain exhibits dysplasia with robust induction of caspase 9-dependent apoptosis. Mechanistically, SNIP1 regulates target genes that promote cell survival and neurogenesis, and its activities are influenced by TGFß and NFκB signaling pathways. Further, SNIP1 facilitates the genomic occupancy of Polycomb complex PRC2 and instructs H3K27me3 turnover at target genes. Depletion of PRC2 is sufficient to reduce apoptosis and brain dysplasia and to partially restore genetic programs in the SNIP1-depleted brain in vivo. These findings suggest a loci-specific regulation of PRC2 and H3K27 marks to toggle cell survival and death in the developing brain.

PP2Ac Deficiency Enhances Tumor Immunogenicity by Activating STING-Type I Interferon Signaling in Glioblastoma.

Glioblastoma (GBM) is an immunologically "cold" tumor that does not respond to current immunotherapy. Here, we demonstrate a fundamental role for the α-isoform of the catalytic subunit of protein phosphatase-2A (PP2Ac) in regulating glioma immunogenicity. Genetic ablation of PP2Ac in glioma cells enhanced double-stranded DNA (dsDNA) production and cGAS-type I IFN signaling, MHC-I expression, and tumor mutational burden. In coculture experiments, PP2Ac deficiency in glioma cells promoted dendritic cell (DC) cross-presentation and clonal expansion of CD8+ T cells. In vivo, PP2Ac depletion sensitized tumors to immune-checkpoint blockade and radiotherapy treatment. Single-cell analysis demonstrated that PP2Ac deficiency increased CD8+ T-cell, natural killer cell, and DC accumulation and reduced immunosuppressive tumor-associated macrophages. Furthermore, loss of PP2Ac increased IFN signaling in myeloid and tumor cells and reduced expression of a tumor gene signature associated with worse patient survival in The Cancer Genome Atlas. Collectively, this study establishes a novel role for PP2Ac in inhibiting dsDNA-cGAS-STING signaling to suppress antitumor immunity in glioma. SIGNIFICANCE: PP2Ac deficiency promotes cGAS-STING signaling in glioma to induce a tumor-suppressive immune microenvironment, highlighting PP2Ac as a potential therapeutic target to enhance tumor immunogenicity and improve response to immunotherapy.

Phosphorylation of 53BP1 by ATM enforce neurodevelopmental programs in cortical organoids.

53BP1 is a well-established DNA damage repair factor recently shown to regulate gene expression and critically influence tumor suppression and neural development. For gene regulation, how 53BP1 is regulated remains unclear. Here, we showed that 53BP1-serine 25 phosphorylation by ATM is required for neural progenitor cell proliferation and neuronal differentiation in cortical organoids. 53BP1-serine 25 phosphorylation dynamics controls 53BP1 target genes for neuronal differentiation and function, cellular response to stress, and apoptosis. Beyond 53BP1, ATM is required for phosphorylation of factors in neuronal differentiation, cytoskeleton, p53 regulation, and ATM, BNDF, and WNT signaling pathways for cortical organoid differentiation. Overall, our data suggest that 53BP1 and ATM control key genetic programs required for human cortical development.

Interrogating bromodomain inhibitor resistance in KMT2A-rearranged leukemia through combinatorial CRISPR screens.

Bromo- and extra-terminal domain inhibitors (BETi) have exhibited therapeutic activities in many cancers. However, the mechanisms controlling BETi response and resistance are not well understood. We conducted genome-wide loss-of-function CRISPR screens using BETi-treated KMT2A-rearranged (KMT2A-r) cell lines. We revealed that Speckle-type POZ protein (SPOP) gene (Speckle Type BTB/POZ Protein) deficiency caused significant BETi resistance, which was further validated in cell lines and xenograft models. Proteomics analysis and a kinase-vulnerability CRISPR screen indicated that cells treated with BETi are sensitive to GSK3 perturbation. Pharmaceutical inhibition of GSK3 reversed the BETi-resistance phenotype. Based on this observation, a combination therapy regimen inhibiting both BET and GSK3 was developed to impede KMT2A-r leukemia progression in patient-derived xenografts in vivo. Our results revealed molecular mechanisms underlying BETi resistance and a promising combination treatment regimen of ABBV-744 and CHIR-98014 by utilizing unique ex vivo and in vivo KMT2A-r PDX models.

How Estrogen, Testosterone, and Sex Differences Influence Serum Immunoglobulin Isotype Patterns in Mice and Humans.

Females often exhibit superior immune responses compared to males toward vaccines and pathogens such as influenza viruses and SARS-CoV-2. To help explain these differences, we first studied serum immunoglobulin isotype patterns in C57BL/6 male and female mice. We focused on IgG2b, an isotype that lends to virus control and that has been previously shown to be elevated in murine females compared to males. Improvements in IgG2b serum levels, and/or IgG2b ratios with other non-IgM isotypes, were observed when: (i) wildtype (WT) female mice were compared to estrogen receptor knockout mice (IgG2b, IgG2b/IgG3, IgG2b/IgG1, and IgG2b/IgA were all higher in WT mice), (ii) unmanipulated female mice were compared to ovariectomized mice (IgG2b/IgA was higher in unmanipulated animals), (iii) female mice were supplemented with estrogen in the context of an inflammatory insult (IgG2b and IgG2b/IgG3 were improved by estrogen supplementation), and (iv) male mice were supplemented with testosterone, a hormone that can convert to estrogen in vivo (IgG2b, IgG2b/IgG3, IgG2b/IgG1, and IgG2b/IgA were all improved by supplementation). We next examined data from three sets of previously described male and female human blood samples. In each case, there were higher IgG2 levels, and/or ratios of IgG2 with non-IgM isotypes, in human females compared to males. The effects of sex and sex hormones in the mouse and human studies were subtle, but frequent, suggesting that sex hormones represent only a fraction of the factors that influence isotype patterns. Examination of the gene loci suggested that upregulation of murine IgG2b or human IgG2 could be mediated by estrogen receptor binding to estrogen response elements and cytosine-adenine (CA) repeats upstream of respective Cγ genes. Given that murine IgG2b and human IgG2 lend to virus control, the isotype biases in females may be sufficient to improve outcomes following vaccination or infection. Future attention to sex hormone levels, and consequent immunoglobulin isotype patterns, in clinical trials are encouraged to support the optimization of vaccine and drug products for male and female hosts.

Spatiotemporal regulation of cholangiocarcinoma growth and dissemination by peritumoral myofibroblasts in a Vcam1-dependent manner.

Intrahepatic cholangiocarcinoma (iCCA) is characterized by its highly desmoplastic stroma. Myofibroblasts (MFs) are present both within the tumor mass (intratumoral MFs, iMFs) and at the tumor border (peritumoral MFs, pMFs). Using a spheroid-based coculture system, we show that the initial iCCA-pMF contact is growth suppressive to the tumor cells. However, prolonged iCCA-pMF interaction elicits significant tumor cell invasion and dissemination. We find that vascular cell adhesion molecule-1 (Vcam1) level is elevated in tumor cells in contact with pMFs but low in disseminated tumor cells both in vitro and in vivo. A gene regulatory network analysis of mouse and patient iCCA tumors and Vcam1 knockout (Vcam1KO) demonstrate a heavy involvement of Vcam1 in epithelial-to-mesenchymal transition. While Vcam1KO has only a limited impact on tumor cell growth in their monoculture, Vcam1KO spheroids exhibit instant dissemination and a severe growth defect when cocultured with pMFs. When transplanted into the liver, Vcam1KO iCCA cells show a similar increase in dissemination but a significant defect in establishing primary and metastatic tumors. Incomplete blocking of Vcam1 in vivo reduces the size but increase the number of metastatic lesions. Overall, our study shows a spatiotemporal regulation of iCCA growth and dissemination by pMFs in a Vcam1-dependent manner.