Tumor-Extrinsic Axl Expression Shapes an Inflammatory Microenvironment Independent of Tumor Cell Promoting Axl Signaling in Hepatocellular Carcinoma

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.

Mapping functional to morphological variation reveals the basis of regional extracellular matrix subversion and nerve invasion in pancreatic cancer

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.

Increased enhancer-promoter interactions during developmental enhancer activation in mammals

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.

Long-term Multimodal Recording Reveals Epigenetic Adaptation Routes in Dormant Breast Cancer Cells

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.

Dynamic enhancer landscapes in human craniofacial development

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.

Logical design of synthetic cis-regulatory DNA for genetic tracing of cell identities and state changes

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.

A spatio-temporally constrained gene regulatory network directed by PBX1/2 acquires limb patterning specificity via HAND2

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.

MYC determines lineage commitment in KRAS-driven primary liver cancer development

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.

Cell Type- and Tissue-specific Enhancers in Craniofacial Development

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 development1-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.

Topologically associating domain boundaries are required for normal genome function

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.

The co-evolution of the genome and epigenome in colorectal cancer

Colorectal malignancies are a leading cause of cancer-related death1 and have undergone extensive genomic study2,3. However, DNA mutations alone do not fully explain malignant transformation4-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.

Cis-regulatory architecture and functions of a gene desert controlling pleiotropic expression and cardiac pacemaker development

 

Gene deserts are extensive genomic regions spanning more than 500 kilobases and enriched near developmental genes. Despite their proposed critical roles in development, the enhancer architecture, and biological functions of most gene deserts in mammalian genomes remain unknown. The Shox2 transcription factor is a key regulator of cardiac pacemaker differentiation and at the genomic level is flanked by a centromeric gene desert harboring a multitude of predicted cis-regulatory modules (CRMs). Using CRISPR-Cas9 we demonstrate a requirement of this gene desert for pleiotropic expression and embryonic survival through full control of Shox2 in the sinoatrial node (SAN). To decode the underlying chromatin architecture and CRM logic, we performed region-specific chromatin conformation capture, epigenomic analysis and transgenic in vivo reporter assays. While the Shox2 regulatory landscape was found partitioned into largely tissue-invariant chromatin loops, we identified a surprising cardiac-specific contact domain within the gene desert. Subsequent transgenic analysis revealed clustering of non-cardiac enhancers with activities in limb, craniofacial or neuronal populations in this domain, pointing to a potential tissue-specific 3D-mechanism for enhancer attenuation. In addition, SAN-chromatin assessment and H3K27ac profiling from human fetal cardiac compartments enabled the identification of a unique cardiac enhancer located outside of the contact domain and essential for robust Shox2 levels during cardiac pacemaker development. In summary, our results serve as a blueprint to investigate gene desert function in context to cardiac defects and identify the Shox2 gene desert as a robust cis-regulatory hub indispensable for pleiotropic patterning, cardiac pacemaker differentiation and embryonic survival.

Funding Acknowledgement

Type of funding sources: Public grant(s) – National budget only. Main funding source(s): Swiss National Science Foundation (SNSF),National Institutes of Health (NIH)

Genome-wide fetalization of enhancer architecture in heart disease

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.

Interferons reshape the 3D conformation and accessibility of macrophage chromatin

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.

Transcriptional network orchestrating regional patterning of cortical progenitors

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.

SMAD4 target genes are part of a transcriptional network that integrates the response to BMP and SHH signaling during early limb bud patterning

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.

Summary: The transcriptional targets of SMAD4 in early limb buds are identified and the largely opposing impact of BMP and SHH signaling on early digit patterning and outgrowth is revealed.

Long non-coding RNA TINCR suppresses metastatic melanoma dissemination by preventing ATF4 translation

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.

Author Correction: An atlas of dynamic chromatin landscapes in mouse fetal development

A Correction to this paper has been published: https://doi.org/10.1038/s41586-020-03089-4.

Phenotypic Mapping of Pathologic Cross-Talk between Glioblastoma and Innate Immune Cells by Synthetic Genetic Tracing

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.

Evolution of Advanced Chronic Lymphoid Leukemia Unveiled by Single-Cell Transcriptomics: A Case Report

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.

Abstract PO-104: Phenotypic mapping of pathological crosstalk between glioblastoma and innate immune cells by synthetic genetic tracing

Glioblastoma is a lethal brain tumor with a high degree of heterogeneity and resistance to therapy. The lack of genetic tracing approaches to selectively identify tumor states and fate transitions limited our understanding of tumor homeostasis. By translating glioblastoma subtype signatures into synthetic genetic tracing cassettes, we investigated tumor heterogeneity at cellular and molecular level, in vitro and in vivo. Synthetic genetic tracing demonstrated that proneural glioblastoma is a hardwired identity, whereas the mesenchymal glioblastoma is an adaptive and metastable state downstream pro-inflammatory and differentiation cues, similarly wired in breast and lung tumors. Importantly, we discovered that innate immune cells divert glioblastoma cells into a proneural-to-mesenchymal transition causal to therapeutic resistance. Systematic phenotypic mapping using synthetic locus control regions (sLCRs) is a simple, automated, and scalable strategy for genetic tracing, generally applicable to study developmental and disease homeostasis. In glioblastoma, it uncovered causality between various (micro)environmental, genetic and pharmacological perturbations and mesenchymal commitment.

Citation Format: Matthias Jürgen Schmitt, Carlos Company, Yuliia Dramaretska, Iros Barozzi, Andreas Göhrig, Sonia Kertalli, Melanie Großmann, Heike Naumann, Jikke Wierikx, Danielle Hulsman, Rainer Glass, Massimo Squatrito, Michela Serresi, Gaetano Gargiulo. Phenotypic mapping of pathological crosstalk between glioblastoma and innate immune cells by synthetic genetic tracing [abstract]. In: Proceedings of the AACR Virtual Special Conference on Tumor Heterogeneity: From Single Cells to Clinical Impact; 2020 Sep 17-18. Philadelphia (PA): AACR; Cancer Res 2020;80(21 Suppl):Abstract nr PO-104.

Exploiting evolutionary steering to induce collateral drug sensitivity in cancer

Drug resistance mediated by clonal evolution is arguably the biggest problem in cancer therapy today. However, evolving resistance to one drug may come at a cost of decreased fecundity or increased sensitivity to another drug. These evolutionary trade-offs can be exploited using 'evolutionary steering' to control the tumour population and delay resistance. However, recapitulating cancer evolutionary dynamics experimentally remains challenging. Here, we present an approach for evolutionary steering based on a combination of single-cell barcoding, large populations of 108-109 cells grown without re-plating, longitudinal non-destructive monitoring of cancer clones, and mathematical modelling of tumour evolution. We demonstrate evolutionary steering in a lung cancer model, showing that it shifts the clonal composition of the tumour in our favour, leading to collateral sensitivity and proliferative costs. Genomic profiling revealed some of the mechanisms that drive evolved sensitivity. This approach allows modelling evolutionary steering strategies that can potentially control treatment resistance.

Targeting the scaffolding role of LSD1 (KDM1A) poises acute myeloid leukemia cells for retinoic acid-induced differentiation

The histone demethylase LSD1 is deregulated in several tumors, including leukemias, providing the rationale for the clinical use of LSD1 inhibitors. In acute promyelocytic leukemia (APL), pharmacological doses of retinoic acid (RA) induce differentiation of APL cells, triggering degradation of the PML-RAR oncogene. APL cells are resistant to LSD1 inhibition or knockout, but targeting LSD1 sensitizes them to physiological doses of RA without altering of PML-RAR levels, and extends survival of leukemic mice upon RA treatment. The combination of RA with LSD1 inhibition (or knockout) is also effective in other non-APL, acute myeloid leukemia (AML) cells. Nonenzymatic activities of LSD1 are essential to block differentiation, while RA with targeting of LSD1 releases a differentiation gene expression program, not strictly dependent on changes in histone H3K4 methylation. Integration of proteomic/epigenomic/mutational studies showed that LSD1 inhibitors alter the recruitment of LSD1-containing complexes to chromatin, inhibiting the interaction between LSD1 and the transcription factor GFI1.

Abstract ES4-2: Single cell transcriptomics reveals multi-step adaptations to endocrine therapy

Resistant tumours are thought to arise from the action of Darwinian selection on genetically heterogenous cancer cell populations. However, simple clonal selection is inadequate to describe the late relapses often characterising luminal breast cancers treated with endocrine therapy (ET), suggesting a more complex interplay between genetic and non-genetic factors. During this talk, I will briefly introduce a series of molecular tools that can be used to study transcriptional changes at the single cell level in relevant models. I will then show how used these tools to dissect the contributions of clonal genetic diversity and transcriptional plasticity during the early and late phases of ET at single-cell resolution. Using single-cell RNA-sequencing and imaging we could disentangle the transcriptional variability of plastic cells and define a rare subpopulation of pre-adapted (PA) cells which undergoes further transcriptomic reprogramming and copy number changes to acquire full resistance. We find evidence for sub-clonal expression of a PA signature in primary tumours and for dominant expression in clustered circulating tumour cells. I will then update the audience on how it is possible to integrate these technologies directly into clinical studies and what developments are happening in the field of single cell technologies. I will finally propose a multi-step model for ET resistance development and advocate the use of stage-specific biomarkers.

Citation Format: L Magnani, I Barozzi, S Hong. Single cell transcriptomics reveals multi-step adaptations to endocrine therapy [abstract]. In: Proceedings of the 2019 San Antonio Breast Cancer Symposium; 2019 Dec 10-14; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2020;80(4 Suppl):Abstract nr ES4-2.

Abstract P4-07-06: Breast cancer subtype classification using NanoString and RNAseq technologies

Background: Immunohistochemistry (IHC) is conventionally used in clinical practice to define Breast Cancer (BC) subtypes, while PAM50 is useful in specific cases. The former is partly subjective and semi-quantitative, the latter is based on gene expression profile. We compared gene expression data obtained by NanoString BC 360™ panel with those obtained by RNAseq on 12 Formalin-Fixed Paraffin-Embedded (FFPE) BC samples.

Methods: RNA was isolated from FFPE tumors and the quality was checked before performing NanoString BC 360™ assay and preparing RNAseq library (NEBNext Ultra II RNA Library Prep Kit, Illumina) following the manufacturers’ instructions. Libraries were sequenced on NextSeq500 (Illumina). Reads were aligned using Kallisto and raw read counts were normalized as Transcripts Per Million (TPM). PAM50 and Tumor Inflammation Signature (TIS) were determined on NanoString data. Two-sided t-test on TPM and on normalized NanoString counts was used to compare luminal and Triple Negative (TN) BCs (p <0.05 statistically significant), considering the 734 genes shared by the 2 assays. Deseq2 (1.22.1) was used to perform a differential gene expression analysis on raw RNAseq counts. Enrichment tests were performed on KEGG pathways (FDR < 0.05). Python and R softwares were used for statistical analyses.

Results: Of the 12 BC samples, 3 were luminal A-like, 3 luminal B (HER2+)-like and 6 TN by IHC (Table 1). Pearson’s correlation coefficient between NanoString counts and RNAseq TPM was 0.72.

Out of the 734 genes shared by both methods, those differentially expressed between luminal and TN BCs, defined based on IHC, were 135 according to RNAseq TPM, 155 according to NanoString counts and 88 according to both methods. These 88 genes belonged to “cell cycle”, “pathways in cancer”, “HIF-1 signaling”, “progesterone-mediated oocyte maturation”, “prostate cancer” and “basal cell carcinoma” KEGG pathways.

When considering the whole transcriptome analyzed by RNAseq, the differential gene expression analysis between luminal and TN/basal BCs was performed according to both IHC and PAM50 classifications. Considering the former, there were 1143 differentially expressed genes out of a total of 20936 (5.5%) and the pathways involved were “fatty acid elongation”, “cell cycle”, “oocyte meiosis”, “p53 signaling” and “alcoholism”. The analysis on the latter classification showed 1608 differentially expressed genes out of 21048 (7.6%), and the enriched pathways were “cell cycle”, “p53 signaling”, “DNA replication”, “alcoholism” and “systemic lupus erythematosus”.

Regarding the TIS, no significant differences were found between TN/basal and luminal BCs, comparing both IHC and the PAM50 subtypes. An expected heterogeneity was found among TNBCs, that showed 3 distinct signatures: basal-like immune activated, basal-like immune-suppressed and luminal androgen receptor.

Conclusions: This preliminary analysis shows a good concordance between NanoString and RNAseq data. The choice of the method must take into account important factors such as the costs, the number of genes analyzed, the requirement of bioinformatic support, the working time. Further analyses on pathway enrichment will be presented.

Table 1. Breast cancer samples characteristics.Sample ID Age Grade T N IHC subtype PAM50 Subtype TIS scoreTNBC signature 0172222aLumALumB7,0NA026031c2aTNBasal7,2BLIA037221cxTNLumA7,4LAR045231c0LumALumA7,5NA053531c1aTNBasal5,5BLIS0675320TNBasal6,2BLIS077931cxLumB HER2+LumB6,7NA0878220LumALumB7,7NA097431a2aTNBasal8,5BLIA106031c1LumB HER2+LumA8,2NA116621b0LumB HER2+LumA6,1NA1241321aTNBasal7,5BLIATIS: Tumor Inflammation Signature; BLIA: Basal-like Immune Activated; BLIS: Basal-like Immune-Suppressed; LAR: Luminal Androgen Receptor.

Citation Format: Andrea Rocca, Sara Ravaioli, Eugenio Fonzi, Iros Barozzi, Ylenia Perone, Luca Magnani, Francesca Pirini, Giovanni Martinelli, Sara Bravaccini. Breast cancer subtype classification using NanoString and RNAseq technologies [abstract]. In: Proceedings of the 2019 San Antonio Breast Cancer Symposium; 2019 Dec 10-14; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2020;80(4 Suppl):Abstract nr P4-07-06.

High-resolution label-free 3D mapping of extracellular pH of single living cells

Dynamic mapping of extracellular pH (pHe) at the single-cell level is critical for understanding the role of H+ in cellular and subcellular processes, with particular importance in cancer. While several pHe sensing techniques have been developed, accessing this information at the single-cell level requires improvement in sensitivity, spatial and temporal resolution. We report on a zwitterionic label-free pH nanoprobe that addresses these long-standing challenges. The probe has a sensitivity > 0.01 units, 2 ms response time, and 50 nm spatial resolution. The platform was integrated into a double-barrel nanoprobe combining pH sensing with feedback-controlled distance dependance via Scanning Ion Conductance Microscopy. This allows for the simultaneous 3D topographical imaging and pHe monitoring of living cancer cells. These classes of nanoprobes were used for real-time high spatiotemporal resolution pHe mapping at the subcellular level and revealed tumour heterogeneity of the peri-cellular environments of melanoma and breast cancer cells.

Single-cell transcriptomics reveals multi-step adaptations to endocrine therapy

Resistant tumours are thought to arise from the action of Darwinian selection on genetically heterogenous cancer cell populations. However, simple clonal selection is inadequate to describe the late relapses often characterising luminal breast cancers treated with endocrine therapy (ET), suggesting a more complex interplay between genetic and non-genetic factors. Here, we dissect the contributions of clonal genetic diversity and transcriptional plasticity during the early and late phases of ET at single-cell resolution. Using single-cell RNA-sequencing and imaging we disentangle the transcriptional variability of plastic cells and define a rare subpopulation of pre-adapted (PA) cells which undergoes further transcriptomic reprogramming and copy number changes to acquire full resistance. We find evidence for sub-clonal expression of a PA signature in primary tumours and for dominant expression in clustered circulating tumour cells. We propose a multi-step model for ET resistance development and advocate the use of stage-specific biomarkers.

Dissection of acute stimulus-inducible nucleosome remodeling in mammalian cells

Accessibility of the genomic regulatory information is largely controlled by the nucleosome-organizing activity of transcription factors (TFs). While stimulus-induced TFs bind to genomic regions that are maintained accessible by lineage-determining TFs, they also increase accessibility of thousands of cis-regulatory elements. Nucleosome remodeling events underlying such changes and their interplay with basal positioning are unknown. Here, we devised a novel quantitative framework discriminating different types of nucleosome remodeling events in micrococcal nuclease ChIP-seq (chromatin immunoprecipitation [ChIP] combined with high-throughput sequencing) data sets and used it to analyze nucleosome dynamics at stimulus-regulated cis-regulatory elements. At enhancers, remodeling preferentially affected poorly positioned nucleosomes while sparing well-positioned nucleosomes flanking the enhancer core, indicating that inducible TFs do not suffice to overrule basal nucleosomal organization maintained by lineage-determining TFs. Remodeling events appeared to be combinatorially driven by multiple TFs, with distinct TFs showing, however, different remodeling efficiencies. Overall, these data provide a systematic view of the impact of stimulation on nucleosome organization and genome accessibility in mammalian cells.

Author Correction: Differential epigenetic reprogramming in response to specific endocrine therapies promotes cholesterol biosynthesis and cellular invasion

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

Control of inducible gene expression links cohesin to hematopoietic progenitor self-renewal and differentiation

Cohesin is important for 3D genome organization. Nevertheless, even the complete removal of cohesin has surprisingly little impact on steady-state gene transcription and enhancer activity. Here we show that cohesin is required for the core transcriptional response of primary macrophages to microbial signals, and for inducible enhancer activity that underpins inflammatory gene expression. Consistent with a role for inflammatory signals in promoting myeloid differentiation of hematopoietic stem and progenitor cells (HPSCs), cohesin mutations in HSPCs led to reduced inflammatory gene expression and increased resistance to differentiation-inducing inflammatory stimuli. These findings uncover an unexpected dependence of inducible gene expression on cohesin, link cohesin with myeloid differentiation, and may help explain the prevalence of cohesin mutations in human acute myeloid leukemia.

Enhancer mapping uncovers phenotypic heterogeneity and evolution in patients with luminal breast cancer

The degree of intrinsic and interpatient phenotypic heterogeneity and its role in tumor evolution is poorly understood. Phenotypic drifts can be transmitted via inheritable transcriptional programs. Cell-type specific transcription is maintained through the activation of epigenetically defined regulatory regions including promoters and enhancers. Here we have annotated the epigenome of 47 primary and metastatic estrogen-receptor (ERα)-positive breast cancer clinical specimens and inferred phenotypic heterogeneity from the regulatory landscape, identifying key regulatory elements commonly shared across patients. Shared regions contain a unique set of regulatory information including the motif for transcription factor YY1. We identify YY1 as a critical determinant of ERα transcriptional activity promoting tumor growth in most luminal patients. YY1 also contributes to the expression of genes mediating resistance to endocrine treatment. Finally, we used H3K27ac levels at active enhancer elements as a surrogate of intra-tumor phenotypic heterogeneity to track the expansion and contraction of phenotypic subpopulations throughout breast cancer progression. By tracking the clonality of SLC9A3R1-positive cells, a bona fide YY1-ERα-regulated gene, we show that endocrine therapies select for phenotypic clones under-represented at diagnosis. Collectively, our data show that epigenetic mechanisms significantly contribute to phenotypic heterogeneity and evolution in systemically treated breast cancer patients.

Thrombopoietin signaling to chromatin elicits rapid and pervasive epigenome remodeling within poised chromatin architectures

Thrombopoietin (TPO) is a critical cytokine regulating hematopoietic stem cell maintenance and differentiation into the megakaryocytic lineage. However, the transcriptional and chromatin dynamics elicited by TPO signaling are poorly understood. Here, we study the immediate early transcriptional and cis-regulatory responses to TPO in hematopoietic stem/progenitor cells (HSPCs) and use this paradigm of cytokine signaling to chromatin to dissect the relation between cis- regulatory activity and chromatin architecture. We show that TPO profoundly alters the transcriptome of HSPCs, with key hematopoietic regulators being transcriptionally repressed within 30 minutes of TPO. By examining cis-regulatory dynamics and chromatin architectures, we demonstrate that these changes are accompanied by rapid and extensive epigenome remodeling of cis-regulatory landscapes that is spatially coordinated within topologically associating domains (TADs). Moreover, TPO-responsive enhancers are spatially clustered and engage in preferential homotypic intra- and inter-TAD interactions that are largely refractory to TPO signaling. By further examining the link between cis-regulatory dynamics and chromatin looping, we show that rapid modulation of cis-regulatory activity is largely independent of chromatin looping dynamics. Finally, we show that, although activated and repressed cis-regulatory elements share remarkably similar DNA sequence compositions, transcription factor binding patterns accurately predict rapid cis-regulatory responses to TPO.

Ultraconserved Enhancers Are Required for Normal Development

Non-coding "ultraconserved" regions containing hundreds of consecutive bases of perfect sequence conservation across mammalian genomes can function as distant-acting enhancers. However, initial deletion studies in mice revealed that loss of such extraordinarily constrained sequences had no immediate impact on viability. Here, we show that ultraconserved enhancers are required for normal development. Focusing on some of the longest ultraconserved sites genome wide, located near the essential neuronal transcription factor Arx, we used genome editing to create an expanded series of knockout mice lacking individual or combinations of ultraconserved enhancers. Mice with single or pairwise deletions of ultraconserved enhancers were viable and fertile but in nearly all cases showed neurological or growth abnormalities, including substantial alterations of neuron populations and structural brain defects. Our results demonstrate the functional importance of ultraconserved enhancers and indicate that remarkably strong sequence conservation likely results from fitness deficits that appear subtle in a laboratory setting.

Distal Limb Patterning Requires Modulation of cis-Regulatory Activities by HOX13

The combinatorial expression of Hox genes along the body axes is a major determinant of cell fate and plays a pivotal role in generating the animal body plan. Loss of HOXA13 and HOXD13 transcription factors (HOX13) leads to digit agenesis in mice, but how HOX13 proteins regulate transcriptional outcomes and confer identity to the distal-most limb cells has remained elusive. Here, we report on the genome-wide profiling of HOXA13 and HOXD13 in vivo binding and changes of the transcriptome and chromatin state in the transition from the early to the late-distal limb developmental program, as well as in Hoxa13^−/−; Hoxd13^−/−limbs. Our results show that proper termination of the early limb transcriptional program and activation of the late-distal limb program are coordinated by the dual action of HOX13 on cis-regulatory modules.

Limb-Enhancer Genie: An accessible resource of accurate enhancer predictions in the developing limb

Epigenomic mapping of enhancer-associated chromatin modifications facilitates the genome-wide discovery of tissue-specific enhancers in vivo. However, reliance on single chromatin marks leads to high rates of false-positive predictions. More sophisticated, integrative methods have been described, but commonly suffer from limited accessibility to the resulting predictions and reduced biological interpretability. Here we present the Limb-Enhancer Genie (LEG), a collection of highly accurate, genome-wide predictions of enhancers in the developing limb, available through a user-friendly online interface. We predict limb enhancers using a combination of >50 published limb-specific datasets and clusters of evolutionarily conserved transcription factor binding sites, taking advantage of the patterns observed at previously in vivo validated elements. By combining different statistical models, our approach outperforms current state-of-the-art methods and provides interpretable measures of feature importance. Our results indicate that including a previously unappreciated score that quantifies tissue-specific nuclease accessibility significantly improves prediction performance. We demonstrate the utility of our approach through in vivo validation of newly predicted elements. Moreover, we describe general features that can guide the type of datasets to include when predicting tissue-specific enhancers genome-wide, while providing an accessible resource to the general biological community and facilitating the functional interpretation of genetic studies of limb malformations.

High constitutive activity of a broad panel of housekeeping and tissue-specific cis-regulatory elements depends on a subset of ETS proteins

Here, Curina et al. found that a subgroup of ETS family transcription factors (the ELF proteins) shows a strong bias in its genomic distribution by binding in very close proximity (<60 nt) to the transcription start sites of housekeeping genes. They show that a limited number of highly active transcription factors can equip cis-regulatory elements with disparate functional roles and cell type specificity with the ability to efficiently promote transcription.

Enhancers and promoters that control the transcriptional output of terminally differentiated cells include cell type-specific and broadly active housekeeping elements. Whether the high constitutive activity of these two groups of cis-regulatory elements relies on entirely distinct or instead also on shared regulators is unknown. By dissecting the cis-regulatory repertoire of macrophages, we found that the ELF subfamily of ETS proteins selectively bound within 60 base pairs (bp) from the transcription start sites of highly active housekeeping genes. ELFs also bound constitutively active, but not poised, macrophage-specific enhancers and promoters. The role of ELFs in promoting high-level constitutive transcription was suggested by multiple evidence: ELF sites enabled robust transcriptional activation by endogenous and minimal synthetic promoters, ELF recruitment was stabilized by the transcriptional machinery, and ELF proteins mediated recruitment of transcriptional and chromatin regulators to core promoters. These data suggest that the co-optation of a limited number of highly active transcription factors represents a broadly adopted strategy to equip both cell type-specific and housekeeping cis-regulatory elements with the ability to efficiently promote transcription.

Acquired CYP19A1 amplification is an early specific mechanism of aromatase inhibitor resistance in ERα metastatic breast cancer

Tumor evolution is shaped by many variables, potentially involving external selective pressures induced by therapies. After surgery, patients with estrogen receptor (ERα)-positive breast cancer are treated with adjuvant endocrine therapy, including selective estrogen receptor modulators (SERMs) and/or aromatase inhibitors (AIs). However, more than 20% of patients relapse within 10 years and eventually progress to incurable metastatic disease. Here we demonstrate that the choice of therapy has a fundamental influence on the genetic landscape of relapsed diseases. We found that 21.5% of AI-treated, relapsed patients had acquired CYP19A1 (encoding aromatase) amplification (CYP19A1amp). Relapsed patients also developed numerous mutations targeting key breast cancer-associated genes, including ESR1 and CYP19A1. Notably, CYP19A1amp cells also emerged in vitro, but only in AI-resistant models. CYP19A1 amplification caused increased aromatase activity and estrogen-independent ERα binding to target genes, resulting in CYP19A1amp cells showing decreased sensitivity to AI treatment. These data suggest that AI treatment itself selects for acquired CYP19A1amp and promotes local autocrine estrogen signaling in AI-resistant metastatic patients.

Reciprocal insulation analysis of Hi-C data shows that TADs represent a functionally but not structurally privileged scale in the hierarchical folding of chromosomes

Understanding how regulatory sequences interact in the context of chromosomal architecture is a central challenge in biology. Chromosome conformation capture revealed that mammalian chromosomes possess a rich hierarchy of structural layers, from multi-megabase compartments to sub-megabase topologically associating domains (TADs) and sub-TAD contact domains. TADs appear to act as regulatory microenvironments by constraining and segregating regulatory interactions across discrete chromosomal regions. However, it is unclear whether other (or all) folding layers share similar properties, or rather TADs constitute a privileged folding scale with maximal impact on the organization of regulatory interactions. Here, we present a novel algorithm named CaTCH that identifies hierarchical trees of chromosomal domains in Hi-C maps, stratified through their reciprocal physical insulation, which is a single and biologically relevant parameter. By applying CaTCH to published Hi-C data sets, we show that previously reported folding layers appear at different insulation levels. We demonstrate that although no structurally privileged folding level exists, TADs emerge as a functionally privileged scale defined by maximal boundary enrichment in CTCF and maximal cell-type conservation. By measuring transcriptional output in embryonic stem cells and neural precursor cells, we show that the likelihood that genes in a domain are coregulated during differentiation is also maximized at the scale of TADs. Finally, we observe that regulatory sequences occur at genomic locations corresponding to optimized mutual interactions at the same scale. Our analysis suggests that the architectural functionality of TADs arises from the interplay between their ability to partition interactions and the specific genomic position of regulatory sequences.

Progressive Loss of Function in a Limb Enhancer during Snake Evolution

The evolution of body shape is thought to be tightly coupled to changes in regulatory sequences, but specific molecular events associated with major morphological transitions in vertebrates have remained elusive. We identified snake-specific sequence changes within an otherwise highly conserved long-range limb enhancer of Sonic hedgehog (Shh). Transgenic mouse reporter assays revealed that the in vivo activity pattern of the enhancer is conserved across a wide range of vertebrates, including fish, but not in snakes. Genomic substitution of the mouse enhancer with its human or fish ortholog results in normal limb development. In contrast, replacement with snake orthologs caused severe limb reduction. Synthetic restoration of a single transcription factor binding site lost in the snake lineage reinstated full in vivo function to the snake enhancer. Our results demonstrate changes in a regulatory sequence associated with a major body plan transition and highlight the role of enhancers in morphological evolution. PAPERCLIP.

Genome-wide compendium and functional assessment of in vivo heart enhancers

Whole-genome sequencing is identifying growing numbers of non-coding variants in human disease studies, but the lack of accurate functional annotations prevents their interpretation. We describe the genome-wide landscape of distant-acting enhancers active in the developing and adult human heart, an organ whose impairment is a predominant cause of mortality and morbidity. Using integrative analysis of >35 epigenomic data sets from mouse and human pre- and postnatal hearts we created a comprehensive reference of >80,000 putative human heart enhancers. To illustrate the importance of enhancers in the regulation of genes involved in heart disease, we deleted the mouse orthologs of two human enhancers near cardiac myosin genes. In both cases, we observe in vivo expression changes and cardiac phenotypes consistent with human heart disease. Our study provides a comprehensive catalogue of human heart enhancers for use in clinical whole-genome sequencing studies and highlights the importance of enhancers for cardiac function.

Identification of non-coding variants has outstripped our ability to annotate and interpret them. Dickel et al. present a compendium of over 80,000 putative human heart enhancers and demonstrate that two conserved enhancers are required for proper cardiac function in mice.

Fishing for Function in the Human Gene Pool

Identification and characterization of causal non-coding variants in human genomes is challenging and requires substantial experimental resources. A new study by Tehranchi et al. describes a cost-effective approach for accurate mapping of molecular quantitative trait loci (QTLs) from pooled samples, a powerful way to link disease-associated changes to molecular functions.

TET2 Regulates Mast Cell Differentiation and Proliferation through Catalytic and Non-catalytic Activities

Dioxygenases of the TET family impact genome functions by converting 5-methylcytosine (5mC) in DNA to 5-hydroxymethylcytosine (5hmC). Here, we identified TET2 as a crucial regulator of mast cell differentiation and proliferation. In the absence of TET2, mast cells showed disrupted gene expression and altered genome-wide 5hmC deposition, especially at enhancers and in the proximity of downregulated genes. Impaired differentiation of Tet2-ablated cells could be relieved or further exacerbated by modulating the activity of other TET family members, and mechanistically it could be linked to the dysregulated expression of C/EBP family transcription factors. Conversely, the marked increase in proliferation induced by the loss of TET2 could be rescued exclusively by re-expression of wild-type or catalytically inactive TET2. Our data indicate that, in the absence of TET2, mast cell differentiation is under the control of compensatory mechanisms mediated by other TET family members, while proliferation is strictly dependent on TET2 expression.

Differential epigenetic reprogramming in response to specific endocrine therapies promotes cholesterol biosynthesis and cellular invasion

Endocrine therapies target the activation of the oestrogen receptor alpha (ERα) via distinct mechanisms, but it is not clear whether breast cancer cells can adapt to treatment using drug-specific mechanisms. Here we demonstrate that resistance emerges via drug-specific epigenetic reprogramming. Resistant cells display a spectrum of phenotypical changes with invasive phenotypes evolving in lines resistant to the aromatase inhibitor (AI). Orthogonal genomics analysis of reprogrammed regulatory regions identifies individual drug-induced epigenetic states involving large topologically associating domains (TADs) and the activation of super-enhancers. AI-resistant cells activate endogenous cholesterol biosynthesis (CB) through stable epigenetic activation in vitro and in vivo. Mechanistically, CB sparks the constitutive activation of oestrogen receptors alpha (ERα) in AI-resistant cells, partly via the biosynthesis of 27-hydroxycholesterol. By targeting CB using statins, ERα binding is reduced and cell invasion is prevented. Epigenomic-led stratification can predict resistance to AI in a subset of ERα-positive patients.

Redox-Mediated Suberoylanilide Hydroxamic Acid Sensitivity in Breast Cancer

AIMS

Vorinostat (suberoylanilide hydroxamic acid; SAHA) is a histone deacetylase inhibitor (HDACi) approved in the clinics for the treatment of T-cell lymphoma and with the potential to be effective also in breast cancer. We investigated the responsiveness to SAHA in human breast primary tumors and cancer cell lines.

RESULTS

We observed a differential response to drug treatment in both human breast primary tumors and cancer cell lines. Gene expression analysis of the breast cancer cell lines revealed that genes involved in cell adhesion and redox pathways, especially glutathione metabolism, were differentially expressed in the cell lines resistant to SAHA compared with the sensitive ones, indicating their possible association with drug resistance mechanisms. Notably, such an association was also observed in breast primary tumors. Indeed, addition of buthionine sulfoximine (BSO), a compound capable of depleting cellular glutathione, significantly enhanced the cytotoxicity of SAHA in both breast cancer cell lines and primary breast tumors.

INNOVATION

We identify and validate transcriptional differences in genes involved in redox pathways, which include potential predictive markers of sensitivity to SAHA.

CONCLUSION

In breast cancer, it could be relevant to evaluate the expression of antioxidant genes that may favor tumor resistance as a factor to consider for potential clinical application and treatment with epigenetic drugs (HDACis).

A dual cis-regulatory code links IRF8 to constitutive and inducible gene expression in macrophages

The transcription factor IRF8 controls both developmental and inflammatory stimulus-inducible genes in macrophages. Mancino et al. found that IRF8 is recruited to distinct sets of DNA consensus sequences before and after lipopolysaccharide (LPS) stimulation. While constitutively expressed IRF8-dependent genes contained only sites mediating basal IRF8/PU.1 recruitment, inducible IRF8-dependent genes contained variable combinations of constitutive and inducible sites.

The transcription factor (TF) interferon regulatory factor 8 (IRF8) controls both developmental and inflammatory stimulus-inducible genes in macrophages, but the mechanisms underlying these two different functions are largely unknown. One possibility is that these different roles are linked to the ability of IRF8 to bind alternative DNA sequences. We found that IRF8 is recruited to distinct sets of DNA consensus sequences before and after lipopolysaccharide (LPS) stimulation. In resting cells, IRF8 was mainly bound to composite sites together with the master regulator of myeloid development PU.1. Basal IRF8–PU.1 binding maintained the expression of a broad panel of genes essential for macrophage functions (such as microbial recognition and response to purines) and contributed to basal expression of many LPS-inducible genes. After LPS stimulation, increased expression of IRF8, other IRFs, and AP-1 family TFs enabled IRF8 binding to thousands of additional regions containing low-affinity multimerized IRF sites and composite IRF–AP-1 sites, which were not premarked by PU.1 and did not contribute to the basal IRF8 cistrome. While constitutively expressed IRF8-dependent genes contained only sites mediating basal IRF8/PU.1 recruitment, inducible IRF8-dependent genes contained variable combinations of constitutive and inducible sites. Overall, these data show at the genome scale how the same TF can be linked to constitutive and inducible gene regulation via distinct combinations of alternative DNA-binding sites.

Comparative evaluation of DNase-seq footprint identification strategies

DNase I is an enzyme preferentially cleaving DNA in highly accessible regions. Recently, Next-Generation Sequencing has been applied to DNase I assays (DNase-seq) to obtain genome-wide maps of these accessible chromatin regions. With high-depth sequencing, DNase I cleavage sites can be identified with base-pair resolution, revealing the presence of protected regions ("footprints"), corresponding to bound molecules on the DNA. Integrating footprint positions close to transcription start sites with motif analysis can reveal the presence of regulatory interactions between specific transcription factors (TFs) and genes. However, this inference heavily relies on the accuracy of the footprint call and on the sequencing depth of the DNase-seq experiment. Using ENCODE data, we comprehensively evaluate the performances of two recent footprint callers (Wellington and DNaseR) and one metric (the Footprint Occupancy Score, or FOS), and assess the consequences of different footprint calls on the reconstruction of TF-TF regulatory networks. We rate Wellington as the method of choice among those tested: not only its predictions are the best in terms of accuracy, but also the properties of the inferred networks are robust against sequencing depth.