A genome-wide relay of signalling-responsive enhancers drives hematopoietic specification

HDI-STARR-seq: Condition-specific enhancer discovery in mouse liver in vivo

STARR-seq and other massively-parallel reporter assays are widely used to discover functional enhancers in transfected cell models, which can be confounded by plasmid vector-induced type-I interferon immune responses and lack the multicellular environment and endogenous chromatin state of complex mammalian tissues. Here, we describe HDI-STARR-seq, which combines STARR-seq plasmid library delivery to the liver, by hydrodynamic tail vein injection (HDI), with reporter RNA transcriptional initiation driven by a minimal Albumin promoter, which we show is essential for mouse liver STARR-seq enhancer activity assayed 7 days after HDI. Importantly, little or no vector-induced innate type-I interferon responses were observed. Comparisons of HDI-STARR-seq activity between male and female mouse livers and in livers from males treated with an activating ligand of the transcription factor CAR ( Nr1i3 ) identified many condition-dependent enhancers linked to condition-specific gene expression. Further, thousands of active liver enhancers were identified using a high complexity STARR-seq library comprised of ∼50,000 genomic regions released by DNase-I digestion of mouse liver nuclei. When compared to stringently inactive library sequences, the active enhancer sequences identified were highly enriched for liver open chromatin regions with activating histone marks (H3K27ac, H3K4me1, H3K4me3), were significantly closer to gene transcriptional start sites, and were significantly depleted of repressive (H3K27me3, H3K9me3) and transcribed region histone marks (H3K36me3). HDI-STARR-seq offers substantial improvements over current methodologies for large scale, functional profiling of enhancers, including condition-dependent enhancers, in liver tissue in vivo, and can be adapted to characterize enhancer activities in a variety of species and tissues by selecting suitable tissue- and species-specific promoter sequences.

Multiplex profiling of developmental cis-regulatory elements with quantitative single-cell expression reporters

The inability to scalably and precisely measure the activity of developmental cis-regulatory elements (CREs) in multicellular systems is a bottleneck in genomics. Here we develop a dual RNA cassette that decouples the detection and quantification tasks inherent to multiplex single-cell reporter assays. The resulting measurement of reporter expression is accurate over multiple orders of magnitude, with a precision approaching the limit set by Poisson counting noise. Together with RNA barcode stabilization via circularization, these scalable single-cell quantitative expression reporters provide high-contrast readouts, analogous to classic in situ assays but entirely from sequencing. Screening >200 regions of accessible chromatin in a multicellular in vitro model of early mammalian development, we identify 13 (8 previously uncharacterized) autonomous and cell-type-specific developmental CREs. We further demonstrate that chimeric CRE pairs generate cognate two-cell-type activity profiles and assess gain- and loss-of-function multicellular expression phenotypes from CRE variants with perturbed transcription factor binding sites. Single-cell quantitative expression reporters can be applied in developmental and multicellular systems to quantitatively characterize native, perturbed and synthetic CREs at scale, with high sensitivity and at single-cell resolution.

A novel cis-regulatory element regulates αD and αA-globin gene expression in chicken erythroid cells

Background

Cis-regulatory elements (CREs) play crucial roles in regulating gene expression during erythroid cell differentiation. Genome-wide erythroid-specific CREs have not been characterized in chicken erythroid cells, which is an organism model used to study epigenetic regulation during erythropoiesis.

Methods

Analysis of public genome-wide accessibility (ATAC-seq) maps, along with transcription factor (TF) motif analysis, CTCF, and RNA Pol II occupancy, as well as transcriptome analysis in fibroblasts and erythroid HD3 cells, were used to characterize erythroid-specific CREs. An α-globin CRE was identified, and its regulatory activity was validated in vitro and in vivo by luciferase activity and genome-editing assays in HD3 cells, respectively. Additionally, circular chromosome conformation capture (UMI-4C) assays were used to distinguish its role in structuring the α-globin domain in erythroid chicken cells.

Results

Erythroid-specific CREs displayed occupancy by erythroid TF binding motifs, CTCF, and RNA Pol II, as well as an association with genes involved in hematopoiesis and cell differentiation. An α-globin CRE, referred to as CRE-2, was identified as exhibiting enhancer activity over αD and αA genes in vitro and in vivo. Induction of terminal erythroid differentiation showed that α-globin CRE-2 is required for the induction of αD and αA. Analysis of TF binding motifs at α-globin CRE-2 shows apparent regulation mediated by GATA-1, YY1, and CTCF binding.

Conclusion

Our findings demonstrate that cell-specific CREs constitute a key mechanism that contributes to the fine-tuning gene regulation of erythroid cell differentiation and provide insights into the annotation and characterization of CREs in chicken cells.

Chromatin priming elements direct tissue-specific gene activity before hematopoietic specification

Tissue-specific gene regulation during development involves the interplay between transcription factors and epigenetic regulators binding to enhancer and promoter elements. The pattern of active enhancers defines the cellular differentiation state. However, developmental gene activation involves a previous step called chromatin priming which is not fully understood. We recently developed a genome-wide functional assay that allowed us to functionally identify enhancer elements integrated in chromatin regulating five stages spanning the in vitro differentiation of embryonic stem cells to blood. We also measured global chromatin accessibility, histone modifications, and transcription factor binding. The integration of these data identified and characterised cis-regulatory elements which become activated before the onset of gene expression, some of which are primed in a signalling-dependent fashion. Deletion of such a priming element leads to a delay in the up-regulation of its associated gene in development. Our work uncovers the details of a complex network of regulatory interactions with the dynamics of early chromatin opening being at the heart of dynamic tissue-specific gene expression control.

Lhx3/4 initiates a cardiopharyngeal-specific transcriptional program in response to widespread FGF signaling

Individual signaling pathways, such as fibroblast growth factors (FGFs), can regulate a plethora of inductive events. According to current paradigms, signal-dependent transcription factors (TFs), such as FGF/MapK-activated Ets family factors, partner with lineage-determining factors to achieve regulatory specificity. However, many aspects of this model have not been rigorously investigated. One key question relates to whether lineage-determining factors dictate lineage-specific responses to inductive signals or facilitate these responses in collaboration with other inputs. We utilize the chordate model Ciona robusta to investigate mechanisms generating lineage-specific induction. Previous studies in C. robusta have shown that cardiopharyngeal progenitor cells are specified through the combined activity of FGF-activated Ets1/2.b and an inferred ATTA-binding transcriptional cofactor. Here, we show that the homeobox TF Lhx3/4 serves as the lineage-determining TF that dictates cardiopharyngeal-specific transcription in response to pleiotropic FGF signaling. Targeted knockdown of Lhx3/4 leads to loss of cardiopharyngeal gene expression. Strikingly, ectopic expression of Lhx3/4 in a neuroectodermal lineage subject to FGF-dependent specification leads to ectopic cardiopharyngeal gene expression in this lineage. Furthermore, ectopic Lhx3/4 expression disrupts neural plate morphogenesis, generating aberrant cell behaviors associated with execution of incompatible morphogenetic programs. Based on these findings, we propose that combinatorial regulation by signal-dependent and lineage-determinant factors represents a generalizable, previously uncategorized regulatory subcircuit we term "cofactor-dependent induction." Integration of this subcircuit into theoretical models will facilitate accurate predictions regarding the impact of gene regulatory network rewiring on evolutionary diversification and disease ontogeny.

From Genotype to Phenotype: How Enhancers Control Gene Expression and Cell Identity in Hematopoiesis

Blood comprises a wide array of specialized cells, all of which share the same genetic information and ultimately derive from the same precursor, the hematopoietic stem cell (HSC). This diversity of phenotypes is underpinned by unique transcriptional programs gradually acquired in the process known as hematopoiesis. Spatiotemporal regulation of gene expression depends on many factors, but critical among them are enhancers-sequences of DNA that bind transcription factors and increase transcription of genes under their control. Thus, hematopoiesis involves the activation of specific enhancer repertoires in HSCs and their progeny, driving the expression of sets of genes that collectively determine morphology and function. Disruption of this tightly regulated process can have catastrophic consequences: in hematopoietic malignancies, dysregulation of transcriptional control by enhancers leads to misexpression of oncogenes that ultimately drive transformation. This review attempts to provide a basic understanding of enhancers and their role in transcriptional regulation, with a focus on normal and malignant hematopoiesis. We present examples of enhancers controlling master regulators of hematopoiesis and discuss the main mechanisms leading to enhancer dysregulation in leukemia and lymphoma.

Genome-wide transcription factor-binding maps reveal cell-specific changes in the regulatory architecture of human HSPCs

Hematopoietic stem and progenitor cells (HSPCs) rely on a complex interplay among transcription factors (TFs) to regulate differentiation into mature blood cells. A heptad of TFs (FLI1, ERG, GATA2, RUNX1, TAL1, LYL1, LMO2) bind regulatory elements in bulk CD34+ HSPCs. However, whether specific heptad-TF combinations have distinct roles in regulating hematopoietic differentiation remains unknown. We mapped genome-wide chromatin contacts (HiC, H3K27ac, HiChIP), chromatin modifications (H3K4me3, H3K27ac, H3K27me3) and 10 TF binding profiles (heptad, PU.1, CTCF, STAG2) in HSPC subsets (stem/multipotent progenitors plus common myeloid, granulocyte macrophage, and megakaryocyte erythrocyte progenitors) and found TF occupancy and enhancer-promoter interactions varied significantly across cell types and were associated with cell-type-specific gene expression. Distinct regulatory elements were enriched with specific heptad-TF combinations, including stem-cell-specific elements with ERG, and myeloid- and erythroid-specific elements with combinations of FLI1, RUNX1, GATA2, TAL1, LYL1, and LMO2. Furthermore, heptad-occupied regions in HSPCs were subsequently bound by lineage-defining TFs, including PU.1 and GATA1, suggesting that heptad factors may prime regulatory elements for use in mature cell types. We also found that enhancers with cell-type-specific heptad occupancy shared a common grammar with respect to TF binding motifs, suggesting that combinatorial binding of TF complexes was at least partially regulated by features encoded in DNA sequence motifs. Taken together, this study comprehensively characterizes the gene regulatory landscape in rare subpopulations of human HSPCs. The accompanying data sets should serve as a valuable resource for understanding adult hematopoiesis and a framework for analyzing aberrant regulatory networks in leukemic cells.

Combining TSS-MPRA and sensitive TSS profile dissimilarity scoring to study the sequence determinants of transcription initiation

Cis-regulatory elements (CREs) can be classified by the shapes of their transcription start site (TSS) profiles, which are indicative of distinct regulatory mechanisms. Massively parallel reporter assays (MPRAs) are increasingly being used to study CRE regulatory mechanisms, yet the degree to which MPRAs replicate individual endogenous TSS profiles has not been determined. Here, we present a new low-input MPRA protocol (TSS-MPRA) that enables measuring TSS profiles of episomal reporters as well as after lentiviral reporter chromatinization. To sensitively compare MPRA and endogenous TSS profiles, we developed a novel dissimilarity scoring algorithm (WIP score) that outperforms the frequently used earth mover's distance on experimental data. Using TSS-MPRA and WIP scoring on 500 unique reporter inserts, we found that short (153 bp) MPRA promoter inserts replicate the endogenous TSS patterns of ∼60% of promoters. Lentiviral reporter chromatinization did not improve fidelity of TSS-MPRA initiation patterns, and increasing insert size frequently led to activation of extraneous TSS in the MPRA that are not active in vivo. We discuss the implications of our findings, which highlight important caveats when using MPRAs to study transcription mechanisms. Finally, we illustrate how TSS-MPRA and WIP scoring can provide novel insights into the impact of transcription factor motif mutations and genetic variants on TSS patterns and transcription levels.

Identification and characterization of enhancer elements controlling cell type-specific and signalling dependent chromatin programming during hematopoietic development

The development of multi-cellular organisms from a single fertilized egg requires to differentially execute the information encoded in our DNA. This complex process is regulated by the interplay of transcription factors with a chromatin environment, both of which provide the epigenetic information maintaining cell-type specific gene expression patterns. Moreover, transcription factors and their target genes form vast interacting gene regulatory networks which can be exquisitely stable. However, all developmental processes originate from pluripotent precursor cell types. The production of terminally differentiated cells from such cells, therefore, requires successive changes of cell fates, meaning that genes relevant for the next stage of differentiation must be switched on and genes not relevant anymore must be switched off. The stimulus for the change of cell fate originates from extrinsic signals which set a cascade of intracellular processes in motion that eventually terminate at the genome leading to changes in gene expression and the development of alternate gene regulatory networks. How developmental trajectories are encoded in the genome and how the interplay between intrinsic and extrinsic processes regulates development is one of the major questions in developmental biology. The development of the hematopoietic system has long served as model to understand how changes in gene regulatory networks drive the differentiation of the various blood cell types. In this review, we highlight the main signals and transcription factors and how they are integrated at the level of chromatin programming and gene expression control. We also highlight recent studies identifying the cis-regulatory elements such as enhancers at the global level and explain how their developmental activity is regulated by the cooperation of cell-type specific and ubiquitous transcription factors with extrinsic signals.

Unraveling the therapeutic mechanisms of dichloroacetic acid in lung cancer through integrated multi-omics approaches: metabolomics and transcriptomics

Objective: The aim of this study was to investigate the molecular mechanisms underlying the therapeutic effects of dichloroacetic acid (DCA) in lung cancer by integrating multi-omics approaches, as the current understanding of DCA's role in cancer treatment remains insufficiently elucidated. Methods: We conducted a comprehensive analysis of publicly available RNA-seq and metabolomic datasets and established a subcutaneous xenograft model of lung cancer in BALB/c nude mice (n = 5 per group) treated with DCA (50 mg/kg, administered via intraperitoneal injection). Metabolomic profiling, gene expression analysis, and metabolite-gene interaction pathway analysis were employed to identify key pathways and molecular players involved in the response to DCA treatment. In vivo evaluation of DCA treatment on tumor growth and MIF gene expression was performed in the xenograft model. Results: Metabolomic profiling and gene expression analysis revealed significant alterations in metabolic pathways, including the Warburg effect and citric acid cycle, and identified the MIF gene as a potential therapeutic target in lung cancer. Our analysis indicated that DCA treatment led to a decrease in MIF gene expression and an increase in citric acid levels in the treatment group. Furthermore, we observed a potential interaction between citric acid and the MIF gene, suggesting a novel mechanism underlying the therapeutic effects of DCA in lung cancer. Conclusion: This study underscores the importance of integrated omics approaches in deciphering the complex molecular mechanisms of DCA treatment in lung cancer. The identification of key metabolic pathways and the novel finding of citric acid elevation, together with its interaction with the MIF gene, provide promising directions for the development of targeted therapeutic strategies and improving clinical outcomes for lung cancer patients.