Search

Abstract Background High-throughput sequencing is generating massive amounts of data at a pace that largely exceeds the throughput of data analysis routines. Here we introduce Fish the ChIPs (FC), a computational pipeline aimed at a broad public of users and designed to perform complete ChIP-Seq data analysis of an unlimited number of samples, thus increasing throughput, reproducibility and saving time. Results Starting from short read sequences, FC performs the following steps: 1) quality controls, 2) alignment to a reference genome, 3) peak calling, 4) genomic annotation, 5) generation of raw signal tracks for visualization on the UCSC and IGV genome browsers. FC exploits some of the fastest and most effective tools today available. Installation on a Mac platform requires very basic computational skills while configuration and usage are supported by a user-friendly graphic user interface. Alternatively, FC can be compiled from the source code on any Unix machine and then run with the possibility of customizing each single parameter through a simple configuration text file that can be generated using a dedicated user-friendly web-form. Considering the execution time, FC can be run on a desktop machine, even though the use of a computer cluster is recommended for analyses of large batches of data. FC is perfectly suited to work with data coming from Illumina Solexa Genome Analyzers or ABI SOLiD and its usage can potentially be extended to any sequencing platform. Conclusions Compared to existing tools, FC has two main advantages that make it suitable for a broad range of users. First of all, it can be installed and run by wet biologists on a Mac machine. Besides it can handle an unlimited number of samples, being convenient for large analyses. In this context, computational biologists can increase reproducibility of their ChIP-Seq data analyses while saving time for downstream analyses. Reviewers This article was reviewed by Gavin Huttley, George Shpakovski and Sarah Teichmann.

Immunity suffers a function deficit during aging, and the incidence of cancer is increased in the elderly. However, most cancer models employ young mice, which are poorly representative of adult cancer patients. We have previously reported that Triple-Therapy (TT), involving antigen-presenting-cell activation by vinorelbine and generation of TCF1 + -stem-cell-like T cells (scTs) by cyclophosphamide significantly improved anti-PD-1 efficacy in anti-PD1-resistant models like Triple-Negative Breast Cancer (TNBC) and Non-Hodgkin's Lymphoma (NHL), due to T-cell-mediated tumor killing. Here, we describe the effect of TT on TNBC growth and on tumor-microenvironment (TME) of young (6-8w, representative of human puberty) versus adult (12 m, representative of 40y-humans) mice. TT-efficacy was similar in young and adults, as CD8 + scTs were only marginally reduced in adults. However, single-cell analyses revealed major differences in the TME: adults had fewer CD4 + scTs, B-naïve and NK-cells, and more memory-B-cells. Cancer-associated-fibroblasts (CAF) with an Extracellular Matrix (ECM) deposition-signature (Matrix-CAFs) were more common in young mice, while pro-inflammatory stromal populations and myofibroblasts were more represented in adults. Matrix-CAFs in adult mice displayed decreased ECM-remodeling abilities, reduced collagen deposition, and a different pattern of interactions with the other cells of the TME. Taken together, our results suggest that age-dependent differences in the TME should be considered when designing preclinical studies., Competing Interests: Competing interests: The authors have declared that no conflict of interest exists. Ethics: Experiments involving animals were approved by the Italian Ministry of Health and have been done in accordance with the applicable Italian laws (D.L.vo 26/14 and following amendments), the Institutional Animal Care and Use Committee, and the institutional guidelines at the European Institute of Oncology, following ARRIVE guidelines (projects 267/2020 and 219/2024)., (© 2025. The Author(s).)

Approximately a quarter of the human genome consists of gene deserts, large regions devoid of genes often located adjacent to developmental genes and thought to contribute to their regulation. However, defining the regulatory functions embedded within these deserts is challenging due to their large size. Here, we explore the cis-regulatory architecture of a gene desert flanking the Shox2 gene, which encodes a transcription factor indispensable for proximal limb, craniofacial, and cardiac pacemaker development. We identify the gene desert as a regulatory hub containing more than 15 distinct enhancers recapitulating anatomical subdomains of Shox2 expression. Ablation of the gene desert leads to embryonic lethality due to Shox2 depletion in the cardiac sinus venosus, caused in part by the loss of a specific distal enhancer. The gene desert is also required for stylopod morphogenesis, mediated via distributed proximal limb enhancers. In summary, our study establishes a multi-layered role of the Shox2 gene desert in orchestrating pleiotropic developmental expression through modular arrangement and coordinated dynamics of tissue-specific enhancers., (© 2024. The Author(s).)

Only a handful of somatic alterations have been linked to endocrine therapy resistance in hormone-dependent breast cancer, potentially explaining ∼40% of relapses. If other mechanisms underlie the evolution of hormone-dependent breast cancer under adjuvant therapy is currently unknown. In this work, we employ functional genomics to dissect the contribution of cis-regulatory elements (CRE) to cancer evolution by focusing on 12 megabases of noncoding DNA, including clonal enhancers, gene promoters, and boundaries of topologically associating domains. Parallel epigenetic perturbation (CRISPRi) in vitro reveals context-dependent roles for many of these CREs, with a specific impact on dormancy entrance and endocrine therapy resistance. Profiling of CRE somatic alterations in a unique, longitudinal cohort of patients treated with endocrine therapies identifies a limited set of noncoding changes potentially involved in therapy resistance. Overall, our data uncover how endocrine therapies trigger the emergence of transient features which could ultimately be exploited to hinder the adaptive process. Significance: This study shows that cells adapting to endocrine therapies undergo changes in the usage or regulatory regions. Dormant cells are less vulnerable to regulatory perturbation but gain transient dependencies which can be exploited to decrease the formation of dormant persisters., (©2024 The Authors; Published by the American Association for Cancer Research.)

Patients with estrogen receptor-positive breast cancer receive adjuvant endocrine therapies (ET) that delay relapse by targeting clinically undetectable micrometastatic deposits. Yet, up to 50% of patients relapse even decades after surgery through unknown mechanisms likely involving dormancy. To investigate genetic and transcriptional changes underlying tumor awakening, we analyzed late relapse patients and longitudinally profiled a rare cohort treated with long-term neoadjuvant ETs until progression. Next, we developed an in vitro evolutionary study to record the adaptive strategies of individual lineages in unperturbed parallel experiments. Our data demonstrate that ETs induce nongenetic cell state transitions into dormancy in a stochastic subset of cells via epigenetic reprogramming. Single lineages with divergent phenotypes awaken unpredictably in the absence of recurrent genetic alterations. Targeting the dormant epigenome shows promising activity against adapting cancer cells. Overall, this study uncovers the contribution of epigenetic adaptation to the evolution of resistance to ETs., Significance: This study advances the understanding of therapy-induced dormancy with potential clinical implications for breast cancer. Estrogen receptor-positive breast cancer cells adapt to endocrine treatment by entering a dormant state characterized by strong heterochromatinization with no recurrent genetic changes. Targeting the epigenetic rewiring impairs the adaptation of cancer cells to ETs. See related commentary by Llinas-Bertran et al., p. 704. This article is featured in Selected Articles from This Issue, p. 695., (©2024 The Authors; Published by the American Association for Cancer Research.)

The activation of the receptor tyrosine kinase Axl by Gas6 is a major driver of tumorigenesis. Despite recent insights, tumor cell-intrinsic and -extrinsic Axl functions are poorly understood in hepatocellular carcinoma (HCC). Thus, we analyzed the cell-specific aspects of Axl in liver cancer cells and in the tumor microenvironment. We show that tumor-intrinsic Axl expression decreased the survival of mice and elevated the number of pulmonary metastases in a model of resection-based tumor recurrence. Axl expression increased the invasion of hepatospheres by the activation of Akt signaling and a partial epithelial-to-mesenchymal transition (EMT). However, the liver tumor burden of Axl +/+ mice induced by diethylnitrosamine plus carbon tetrachloride was reduced compared to systemic Axl -/- mice. Tumors of Axl +/+ mice were highly infiltrated with cytotoxic cells, suggesting a key immune-modulatory role of Axl. Interestingly, hepatocyte-specific Axl deficiency did not alter T cell infiltration, indicating that these changes are independent of tumor cell-intrinsic Axl. In this context, we observed an upregulation of multiple chemokines in Axl +/+ compared to Axl -/- tumors, correlating with HCC patient data. In line with this, Axl is associated with a cytotoxic immune signature in HCC patients. Together these data show that tumor-intrinsic Axl expression fosters progression, while tumor-extrinsic Axl expression shapes an inflammatory microenvironment.

Intratumor morphological heterogeneity of pancreatic ductal adenocarcinoma (PDAC) predicts clinical outcomes but is only partially understood at the molecular level. To elucidate the gene expression programs underpinning intratumor morphological variation in PDAC, we investigated and deconvoluted at single cell level the molecular profiles of histologically distinct clusters of PDAC cells. We identified three major morphological and functional variants that co-exist in varying proportions in all PDACs, display limited genetic diversity, and are associated with a distinct organization of the extracellular matrix: a glandular variant with classical ductal features; a transitional variant displaying abortive ductal structures and mixed endodermal and myofibroblast-like gene expression; and a poorly differentiated variant lacking ductal features and basement membrane, and showing neuronal lineage priming. Ex vivo and in vitro evidence supports the occurrence of dynamic transitions among these variants in part influenced by extracellular matrix composition and stiffness and associated with local, specifically neural, invasion., Competing Interests: Declaration of interests L.M. reports research funding from Bristol-Myers Squibb and Roche, and consulting fees from Tethis, which are all unrelated to the content of this study., (Copyright © 2024 Elsevier Inc. All rights reserved.)

Remote enhancers are thought to interact with their target promoters via physical proximity, yet the importance of this proximity for enhancer function remains unclear. Here we investigate the three-dimensional (3D) conformation of enhancers during mammalian development by generating high-resolution tissue-resolved contact maps for nearly a thousand enhancers with characterized in vivo activities in ten murine embryonic tissues. Sixty-one percent of developmental enhancers bypass their neighboring genes, which are often marked by promoter CpG methylation. The majority of enhancers display tissue-specific 3D conformations, and both enhancer-promoter and enhancer-enhancer interactions are moderately but consistently increased upon enhancer activation in vivo. Less than 14% of enhancer-promoter interactions form stably across tissues; however, these invariant interactions form in the absence of the enhancer and are likely mediated by adjacent CTCF binding. Our results highlight the general importance of enhancer-promoter physical proximity for developmental gene activation in mammals., (© 2024. The Author(s), under exclusive licence to Springer Nature America, Inc.)

The genetic basis of human facial variation and craniofacial birth defects remains poorly understood. Distant-acting transcriptional enhancers control the fine-tuned spatiotemporal expression of genes during critical stages of craniofacial development. However, a lack of accurate maps of the genomic locations and cell type-resolved activities of craniofacial enhancers prevents their systematic exploration in human genetics studies. Here, we combine histone modification, chromatin accessibility, and gene expression profiling of human craniofacial development with single-cell analyses of the developing mouse face to define the regulatory landscape of facial development at tissue- and single cell-resolution. We provide temporal activity profiles for 14,000 human developmental craniofacial enhancers. We find that 56% of human craniofacial enhancers share chromatin accessibility in the mouse and we provide cell population- and embryonic stage-resolved predictions of their in vivo activity. Taken together, our data provide an expansive resource for genetic and developmental studies of human craniofacial development., (© 2024. The Author(s).)

Descriptive data are rapidly expanding in biomedical research. Instead, functional validation methods with sufficient complexity remain underdeveloped. Transcriptional reporters allow experimental characterization and manipulation of developmental and disease cell states, but their design lacks flexibility. Here, we report logical design of synthetic cis-regulatory DNA (LSD), a computational framework leveraging phenotypic biomarkers and trans-regulatory networks as input to design reporters marking the activity of selected cellular states and pathways. LSD uses bulk or single-cell biomarkers and a reference genome or custom cis-regulatory DNA datasets with user-defined boundary regions. By benchmarking validated reporters, we integrate LSD with a computational ranking of phenotypic specificity of putative cis-regulatory DNA. Experimentally, LSD-designed reporters targeting a wide range of cell states are functional without minimal promoters. Applied to broadly expressed genes from human and mouse tissues, LSD generates functional housekeeper-like sLCRs compatible with size constraints of AAV vectors for gene therapy applications. A mesenchymal glioblastoma reporter designed by LSD outperforms previously validated ones and canonical cell surface markers. In genome-scale CRISPRa screens, LSD facilitates the discovery of known and novel bona fide cell-state drivers. Thus, LSD captures core principles of cis-regulation and is broadly applicable to studying complex cell states and mechanisms of transcriptional regulation., (© 2024. The Author(s).)

A lingering question in developmental biology has centered on how transcription factors with widespread distribution in vertebrate embryos can perform tissue-specific functions. Here, using the murine hindlimb as a model, we investigate the elusive mechanisms whereby PBX TALE homeoproteins, viewed primarily as HOX cofactors, attain context-specific developmental roles despite ubiquitous presence in the embryo. We first demonstrate that mesenchymal-specific loss of PBX1/2 or the transcriptional regulator HAND2 generates similar limb phenotypes. By combining tissue-specific and temporally controlled mutagenesis with multi-omics approaches, we reconstruct a gene regulatory network (GRN) at organismal-level resolution that is collaboratively directed by PBX1/2 and HAND2 interactions in subsets of posterior hindlimb mesenchymal cells. Genome-wide profiling of PBX1 binding across multiple embryonic tissues further reveals that HAND2 interacts with subsets of PBX-bound regions to regulate limb-specific GRNs. Our research elucidates fundamental principles by which promiscuous transcription factors cooperate with cofactors that display domain-restricted localization to instruct tissue-specific developmental programs., (© 2023. The Author(s).)

Background & Aims: Primary liver cancer (PLC) comprises hepatocellular carcinoma (HCC) and intrahepatic cholangiocarcinoma (iCCA), two frequent and lethal tumour types that differ regarding their tumour biology and responses to cancer therapies. Liver cells harbour a high degree of cellular plasticity and can give rise to either HCC or iCCA. However, little is known about the cell-intrinsic mechanisms directing an oncogenically transformed liver cell to either HCC or iCCA. The scope of this study was to identify cell-intrinsic factors determining lineage commitment in PLC., Methods: Cross-species transcriptomic and epigenetic profiling was applied to murine HCCs and iCCAs and to two human PLC cohorts. Integrative data analysis comprised epigenetic Landscape In Silico deletion Analysis (LISA) of transcriptomic data and Hypergeometric Optimization of Motif EnRichment (HOMER) analysis of chromatin accessibility data. Identified candidate genes were subjected to functional genetic testing in non-germline genetically engineered PLC mouse models (shRNAmir knockdown or overexpression of full-length cDNAs)., Results: Integrative bioinformatic analyses of transcriptomic and epigenetic data pinpointed the Forkhead-family transcription factors FOXA1 and FOXA2 as MYC-dependent determination factors of the HCC lineage. Conversely, the ETS family transcription factor ETS1 was identified as a determinant of the iCCA lineage, which was found to be suppressed by MYC during HCC development. Strikingly, shRNA-mediated suppression of FOXA1 and FOXA2 with concomitant ETS1 expression fully switched HCC to iCCA development in PLC mouse models., Conclusions: The herein reported data establish MYC as a key determinant of lineage commitment in PLC and provide a molecular explanation why common liver-damaging risk factors such as alcoholic or non-alcoholic steatohepatitis can lead to either HCC or iCCA., Impact and Implications: Liver cancer is a major health problem and comprises hepatocellular carcinoma (HCC) and intrahepatic cholangiocarcinoma (iCCA), two frequent and lethal tumour types that differ regarding their morphology, tumour biology, and responses to cancer therapies. We identified the transcription factor and oncogenic master regulator MYC as a switch between HCC and iCCA development. When MYC levels are high at the time point when a hepatocyte becomes a tumour cell, an HCC is growing out. Conversely, if MYC levels are low at this time point, the result is the outgrowth of an iCCA. Our study provides a molecular explanation why common liver-damaging risk factors such as alcoholic or non-alcoholic steatohepatitis can lead to either HCC or iCCA. Furthermore, our data harbour potential for the development of better PLC therapies., (Copyright © 2023 The Authors. Published by Elsevier B.V. All rights reserved.)

The genetic basis of craniofacial birth defects and general variation in human facial shape remains poorly understood. Distant-acting transcriptional enhancers are a major category of non-coding genome function and have been shown to control the fine-tuned spatiotemporal expression of genes during critical stages of craniofacial development 1-3 . However, a lack of accurate maps of the genomic location and cell type-specific in vivo activities of all craniofacial enhancers prevents their systematic exploration in human genetics studies. Here, we combined histone modification and chromatin accessibility profiling from different stages of human craniofacial development with single-cell analyses of the developing mouse face to create a comprehensive catalogue of the regulatory landscape of facial development at tissue- and single cell-resolution. In total, we identified approximately 14,000 enhancers across seven developmental stages from weeks 4 through 8 of human embryonic face development. We used transgenic mouse reporter assays to determine the in vivo activity patterns of human face enhancers predicted from these data. Across 16 in vivo validated human enhancers, we observed a rich diversity of craniofacial subregions in which these enhancers are active in vivo . To annotate the cell type specificities of human-mouse conserved enhancers, we performed single-cell RNA-seq and single-nucleus ATAC-seq of mouse craniofacial tissues from embryonic days e11.5 to e15.5. By integrating these data across species, we find that the majority (56%) of human craniofacial enhancers are functionally conserved in mice, providing cell type- and embryonic stage-resolved predictions of their in vivo activity profiles. Using retrospective analysis of known craniofacial enhancers in combination with single cell-resolved transgenic reporter assays, we demonstrate the utility of these data for predicting the in vivo cell type specificity of enhancers. Taken together, our data provide an expansive resource for genetic and developmental studies of human craniofacial development., Competing Interests: Declaration of Interests Bing Ren is a co-founder of Arima Genomics, Inc, and Epigenome Technologies, Inc.

Topologically associating domain (TAD) boundaries partition the genome into distinct regulatory territories. Anecdotal evidence suggests that their disruption may interfere with normal gene expression and cause disease phenotypes 1-3 , but the overall extent to which this occurs remains unknown. Here we demonstrate that targeted deletions of TAD boundaries cause a range of disruptions to normal in vivo genome function and organismal development. We used CRISPR genome editing in mice to individually delete eight TAD boundaries (11-80 kb in size) from the genome. All deletions examined resulted in detectable molecular or organismal phenotypes, which included altered chromatin interactions or gene expression, reduced viability, and anatomical phenotypes. We observed changes in local 3D chromatin architecture in 7 of 8 (88%) cases, including the merging of TADs and altered contact frequencies within TADs adjacent to the deleted boundary. For 5 of 8 (63%) loci examined, boundary deletions were associated with increased embryonic lethality or other developmental phenotypes. For example, a TAD boundary deletion near Smad3/Smad6 caused complete embryonic lethality, while a deletion near Tbx5/Lhx5 resulted in a severe lung malformation. Our findings demonstrate the importance of TAD boundary sequences for in vivo genome function and reinforce the critical need to carefully consider the potential pathogenicity of noncoding deletions affecting TAD boundaries in clinical genetics screening., (© 2023. The Author(s).)

Colorectal malignancies are a leading cause of cancer-related death 1  and have undergone extensive genomic study 2,3 . However, DNA mutations alone do not fully explain malignant transformation 4-7 . Here we investigate the co-evolution of the genome and epigenome of colorectal tumours at single-clone resolution using spatial multi-omic profiling of individual glands. We collected 1,370 samples from 30 primary cancers and 8 concomitant adenomas and generated 1,207 chromatin accessibility profiles, 527 whole genomes and 297 whole transcriptomes. We found positive selection for DNA mutations in chromatin modifier genes and recurrent somatic chromatin accessibility alterations, including in regulatory regions of cancer driver genes that were otherwise devoid of genetic mutations. Genome-wide alterations in accessibility for transcription factor binding involved CTCF, downregulation of interferon and increased accessibility for SOX and HOX transcription factor families, suggesting the involvement of developmental genes during tumourigenesis. Somatic chromatin accessibility alterations were heritable and distinguished adenomas from cancers. Mutational signature analysis showed that the epigenome in turn influences the accumulation of DNA mutations. This study provides a map of genetic and epigenetic tumour heterogeneity, with fundamental implications for understanding colorectal cancer biology., (© 2022. The Author(s).)

Heart disease is associated with re-expression of key transcription factors normally active only during prenatal development of the heart. However, the impact of this reactivation on the regulatory landscape in heart disease is unclear. Here, we use RNA-seq and ChIP-seq targeting a histone modification associated with active transcriptional enhancers to generate genome-wide enhancer maps from left ventricle tissue from up to 26 healthy controls, 18 individuals with idiopathic dilated cardiomyopathy (DCM), and five fetal hearts. Healthy individuals have a highly reproducible epigenomic landscape, consisting of more than 33,000 predicted heart enhancers. In contrast, we observe reproducible disease-associated changes in activity at 6,850 predicted heart enhancers. Combined analysis of adult and fetal samples reveals that the heart disease epigenome and transcriptome both acquire fetal-like characteristics, with 3,400 individual enhancers sharing fetal regulatory properties. We also provide a comprehensive data resource (http://heart.lbl.gov) for the mechanistic exploration of DCM etiology., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2022 The Authors. Published by Elsevier Inc. All rights reserved.)

Engagement of macrophages in innate immune responses is directed by type I and type II interferons (IFN-I and IFN-γ, respectively). IFN triggers drastic changes in cellular transcriptomes, executed by JAK-STAT signal transduction and the transcriptional control of interferon-stimulated genes (ISG) by STAT transcription factors. Here, we study the immediate-early nuclear response to IFN-I and IFN-γ in murine macrophages. We show that the mechanism of gene control by both cytokines includes a rapid increase of DNA accessibility and rearrangement of the 3D chromatin contacts particularly between open chromatin of ISG loci. IFN-stimulated gene factor 3 (ISGF3), the major transcriptional regulator of ISG, controlled homeostatic and, most notably, induced-state DNA accessibility at a subset of ISG. Increases in DNA accessibility correlated with the appearance of activating histone marks at surrounding nucleosomes. Collectively our data emphasize changes in the three-dimensional nuclear space and epigenome as an important facet of transcriptional control by the IFN-induced JAK-STAT pathway., Competing Interests: The authors declare no competing interests., (© 2022 The Author(s).)

We uncovered a transcription factor (TF) network that regulates cortical regional patterning in radial glial stem cells. Screening the expression of hundreds of TFs in the developing mouse cortex identified 38 TFs that are expressed in gradients in the ventricular zone (VZ). We tested whether their cortical expression was altered in mutant mice with known patterning defects ( Emx2, Nr2f1 , and Pax6 ), which enabled us to define a cortical regionalization TF network (CRTFN). To identify genomic programming underlying this network, we performed TF ChIP-seq and chromatin-looping conformation to identify enhancer-gene interactions. To map enhancers involved in regional patterning of cortical progenitors, we performed assays for epigenomic marks and DNA accessibility in VZ cells purified from wild-type and patterning mutant mice. This integrated approach has identified a CRTFN and VZ enhancers involved in cortical regional patterning in the mouse., Competing Interests: Competing interest statement: J.L.R. is a cofounder, stockholder, and currently on the scientific board of Neurona, a company studying the potential therapeutic use of interneuron transplantation., (Copyright © 2021 the Author(s). Published by PNAS.)

SMAD4 regulates gene expression in response to BMP and TGFβ signal transduction, and is required for diverse morphogenetic processes, but its target genes have remained largely elusive. Here, we identify the SMAD4 target genes in mouse limb buds using an epitope-tagged Smad4 allele for ChIP-seq analysis in combination with transcription profiling. This analysis shows that SMAD4 predominantly mediates BMP signal transduction during early limb bud development. Unexpectedly, the expression of cholesterol biosynthesis enzymes is precociously downregulated and intracellular cholesterol levels are reduced in Smad4-deficient limb bud mesenchymal progenitors. Most importantly, our analysis reveals a predominant function of SMAD4 in upregulating target genes in the anterior limb bud mesenchyme. Analysis of differentially expressed genes shared between Smad4- and Shh-deficient limb buds corroborates this function of SMAD4 and also reveals the repressive effect of SMAD4 on posterior genes that are upregulated in response to SHH signaling. This analysis uncovers opposing trans-regulatory inputs from SHH- and SMAD4-mediated BMP signal transduction on anterior and posterior gene expression during the digit patterning and outgrowth in early limb buds., Competing Interests: Competing interests The authors declare no competing or financial interests., (© 2021. Published by The Company of Biologists Ltd.)

Interactions between the splicing machinery and RNA polymerase II increase protein-coding gene transcription. Similarly, exons and splicing signals of enhancer-generated long noncoding RNAs (elncRNAs) augment enhancer activity. However, elncRNAs are inefficiently spliced, suggesting that, compared with protein-coding genes, they contain qualitatively different exons with a limited ability to drive splicing. We show here that the inefficiently spliced first exons of elncRNAs as well as promoter-antisense long noncoding RNAs (pa-lncRNAs) in human and mouse cells trigger a transcription termination checkpoint that requires WDR82, an RNA polymerase II-binding protein, and its RNA-binding partner of previously unknown function, ZC3H4. We propose that the first exons of elncRNAs and pa-lncRNAs are an intrinsic component of a regulatory mechanism that, on the one hand, maximizes the activity of these cis-regulatory elements by recruiting the splicing machinery and, on the other, contains elements that suppress pervasive extragenic transcription.

Transition from proliferative-to-invasive phenotypes promotes metastasis and therapy resistance in melanoma. Reversion of the invasive phenotype, however, is challenged by the poor understanding of mechanisms underlying its maintenance. Here, we report that the lncRNA TINCR is down-regulated in metastatic melanoma and its silencing increases the expression levels of invasive markers, in vitro migration, in vivo tumor growth, and resistance to BRAF and MEK inhibitors. The critical mediator is ATF4, a central player of the integrated stress response (ISR), which is activated in TINCR-depleted cells in the absence of starvation and eIF2α phosphorylation. TINCR depletion increases global protein synthesis and induces translational reprogramming, leading to increased translation of mRNAs encoding ATF4 and other ISR proteins. Strikingly, re-expression of TINCR in metastatic melanoma suppresses the invasive phenotype, reduces numbers of tumor-initiating cells and metastasis formation, and increases drug sensitivity. Mechanistically, TINCR interacts with mRNAs associated with the invasive phenotype, including ATF4, preventing their binding to ribosomes. Thus, TINCR is a suppressor of the melanoma invasive phenotype, which functions in nutrient-rich conditions by repressing translation of selected ISR RNAs., (© 2021 The Authors. Published under the terms of the CC BY NC ND 4.0 license.)

Glioblastoma is a lethal brain tumor that exhibits heterogeneity and resistance to therapy. Our understanding of tumor homeostasis is limited by a lack of genetic tools to selectively identify tumor states and fate transitions. Here, we use glioblastoma subtype signatures to construct synthetic genetic tracing cassettes and investigate tumor heterogeneity at cellular and molecular levels, in vitro and in vivo . Through synthetic locus control regions, we demonstrate that proneural glioblastoma is a hardwired identity, whereas mesenchymal glioblastoma is an adaptive and metastable cell state driven by proinflammatory and differentiation cues and DNA damage, but not hypoxia. Importantly, we discovered that innate immune cells divert glioblastoma cells to a proneural-to-mesenchymal transition that confers therapeutic resistance. Our synthetic genetic tracing methodology is simple, scalable, and widely applicable to study homeostasis in development and diseases. In glioblastoma, the method causally links distinct (micro)environmental, genetic, and pharmacologic perturbations and mesenchymal commitment. SIGNIFICANCE: Glioblastoma is heterogeneous and incurable. Here, we designed synthetic reporters to reflect the transcriptional output of tumor cell states and signaling pathways' activity. This method is generally applicable to study homeostasis in normal tissues and diseases. In glioblastoma, synthetic genetic tracing causally connects cellular and molecular heterogeneity to therapeutic responses. This article is highlighted in the In This Issue feature, p. 521 ., (©2020 American Association for Cancer Research.)

Genetic and transcriptional heterogeneity of Chronic lymphocytic leukaemia (CLL) limits prevention of disease progression. Longitudinal single-cell transcriptomics represents the state-of-the-art method to profile the disease heterogeneity at diagnosis and to inform about disease evolution. Here, we apply single-cell RNA-seq to a CLL case, sampled at diagnosis and relapse, that was treated with FCR (Fludarabine, Cyclophosphamide, Rituximab) and underwent a dramatic decrease in CD19 expression during disease progression. Computational analyses revealed a major switch in clones' dominance during treatment. The clone that expanded at relapse showed 17p and 3p chromosomal deletions, and up-regulation of pathways related to motility, cytokine signaling and antigen presentation. Single-cell RNA-seq uniquely revealed that this clone was already present at low frequency at diagnosis, and it displays feature of plasma cell differentiation, consistent with a more aggressive phenotype. This study shows the benefit of single-cell profiling of CLL heterogeneity at diagnosis, to identify clones that might otherwise not be recognized and to determine the best treatment options., (Copyright © 2020 Ostasov, Robertson, Piazza, Datta, Apperley, Houdova, Lysak, Holubova, Tesarova, Caputo and Barozzi.)

Enhancers are important non-coding elements, but they have traditionally been hard to characterize experimentally. The development of massively parallel assays allows the characterization of large numbers of enhancers for the first time. Here, we developed a framework using Drosophila STARR-seq to create shape-matching filters based on meta-profiles of epigenetic features. We integrated these features with supervised machine-learning algorithms to predict enhancers. We further demonstrated that our model could be transferred to predict enhancers in mammals. We comprehensively validated the predictions using a combination of in vivo and in vitro approaches, involving transgenic assays in mice and transduction-based reporter assays in human cell lines (153 enhancers in total). The results confirmed that our model can accurately predict enhancers in different species without re-parameterization. Finally, we examined the transcription factor binding patterns at predicted enhancers versus promoters. We demonstrated that these patterns enable the construction of a secondary model that effectively distinguishes enhancers and promoters.

The Encyclopedia of DNA Elements (ENCODE) project has established a genomic resource for mammalian development, profiling a diverse panel of mouse tissues at 8 developmental stages from 10.5 days after conception until birth, including transcriptomes, methylomes and chromatin states. Here we systematically examined the state and accessibility of chromatin in the developing mouse fetus. In total we performed 1,128 chromatin immunoprecipitation with sequencing (ChIP-seq) assays for histone modifications and 132 assay for transposase-accessible chromatin using sequencing (ATAC-seq) assays for chromatin accessibility across 72 distinct tissue-stages. We used integrative analysis to develop a unified set of chromatin state annotations, infer the identities of dynamic enhancers and key transcriptional regulators, and characterize the relationship between chromatin state and accessibility during developmental gene regulation. We also leveraged these data to link enhancers to putative target genes and demonstrate tissue-specific enrichments of sequence variants associated with disease in humans. The mouse ENCODE data sets provide a compendium of resources for biomedical researchers and achieve, to our knowledge, the most comprehensive view of chromatin dynamics during mammalian fetal development to date.