Predicting gene expression state and prioritizing putative enhancers using 5hmC signal.

BACKGROUND: Like its parent base 5-methylcytosine (5mC), 5-hydroxymethylcytosine (5hmC) is a direct epigenetic modification of cytosines in the context of CpG dinucleotides. 5hmC is the most abundant oxidized form of 5mC, generated through the action of TET dioxygenases at gene bodies of actively-transcribed genes and at active or lineage-specific enhancers. Although such enrichments are reported for 5hmC, to date, predictive models of gene expression state or putative regulatory regions for genes using 5hmC have not been developed. RESULTS: Here, by using only 5hmC enrichment in genic regions and their vicinity, we develop neural network models that predict gene expression state across 49 cell types. We show that our deep neural network models distinguish high vs low expression state utilizing only 5hmC levels and these predictive models generalize to unseen cell types. Further, in order to leverage 5hmC signal in distal enhancers for expression prediction, we employ an Activity-by-Contact model and also develop a graph convolutional neural network model with both utilizing Hi-C data and 5hmC enrichment to prioritize enhancer-promoter links. These approaches identify known and novel putative enhancers for key genes in multiple immune cell subsets. CONCLUSIONS: Our work highlights the importance of 5hmC in gene regulation through proximal and distal mechanisms and provides a framework to link it to genome function. With the recent advances in 6-letter DNA sequencing by short and long-read techniques, profiling of 5mC and 5hmC may be done routinely in the near future, hence, providing a broad range of applications for the methods developed here.

Development of a New Model System to Study Long-Distance Interactions Supported by Architectural Proteins.

Chromatin architecture is critical for the temporal and tissue-specific activation of genes that determine eukaryotic development. The functional interaction between enhancers and promoters is controlled by insulators and tethering elements that support specific long-distance interactions. However, the mechanisms of the formation and maintenance of long-range interactions between genome regulatory elements remain poorly understood, primarily due to the lack of convenient model systems. Drosophila became the first model organism in which architectural proteins that determine the activity of insulators were described. In Drosophila, one of the best-studied DNA-binding architectural proteins, Su(Hw), forms a complex with Mod(mdg4)-67.2 and CP190 proteins. Using a combination of CRISPR/Cas9 genome editing and attP-dependent integration technologies, we created a model system in which the promoters and enhancers of two reporter genes are separated by 28 kb. In this case, enhancers effectively stimulate reporter gene promoters in cis and trans only in the presence of artificial Su(Hw) binding sites (SBS), in both constructs. The expression of the mutant Su(Hw) protein, which cannot interact with CP190, and the mutation inactivating Mod(mdg4)-67.2, lead to the complete loss or significant weakening of enhancer-promoter interactions, respectively. The results indicate that the new model system effectively identifies the role of individual subunits of architectural protein complexes in forming and maintaining specific long-distance interactions in the D. melanogaster model.

Genomic context sensitizes regulatory elements to genetic disruption.

Genomic context critically modulates regulatory function but is difficult to manipulate systematically. The murine insulin-like growth factor 2 (Igf2)/H19 locus is a paradigmatic model of enhancer selectivity, whereby CTCF occupancy at an imprinting control region directs downstream enhancers to activate either H19 or Igf2. We used synthetic regulatory genomics to repeatedly replace the native locus with 157-kb payloads, and we systematically dissected its architecture. Enhancer deletion and ectopic delivery revealed previously uncharacterized long-range regulatory dependencies at the native locus. Exchanging the H19 enhancer cluster with the Sox2 locus control region (LCR) showed that the H19 enhancers relied on their native surroundings while the Sox2 LCR functioned autonomously. Analysis of regulatory DNA actuation across cell types revealed that these enhancer clusters typify broader classes of context sensitivity genome wide. These results show that unexpected dependencies influence even well-studied loci, and our approach permits large-scale manipulation of complete loci to investigate the relationship between regulatory architecture and function.

LINC00887 Acts as an Enhancer RNA to Promote Medullary Thyroid Carcinoma Progression by Binding with FOXQ1.

BACKGROUND: Medullary thyroid carcinoma (MTC) is a rare but aggressive endocrine malignancy that originates from the parafollicular C cells of the thyroid gland. Enhancer RNAs (eRNAs) are non-coding RNAs transcribed from enhancer regions, which are critical regulators of tumorigenesis. However, the roles and regulatory mechanisms of eRNAs in MTC remain poorly understood. This study aims to identify key eRNAs regulating the malignant phenotype of MTC and to uncover transcription factors involved in the regulation of key eRNAs. METHODS: GSE32662 and GSE114068 were used for the identification of differentially expressed genes, eRNAs, enhancers and enhancer-regulated genes in MTC. Metascape and the transcription factor affinity prediction method were used for gene function enrichment and transcription factor prediction, respectively. qRT-PCR was used to detect gene transcription levels. ChIP-qPCR was used to assess the binding of histone H3 lysine 27 acetylation (H3K27ac)-enriched regions to anti- H3K27ac. RIP-qPCR was used to detect the binding between FOXQ1 and LINC00887. CCK8 and Transwell were performed to measure the proliferation and invasion of MTC cells, respectively. Intracellular reactive oxygen species (ROS) levels were quantified using a ROS assay kit. RESULTS: Four eRNAs (H1FX-AS1, LINC00887, MCM3AP-AS1 and A1BG-AS1) were screened, among which LINC00887 was the key eRNA promoting the proliferation and invasion of MTC cells. A total of 135 genes controlled by LINC00887-regulated enhancers were identified; among them, BCL2, PRDX1, SFTPD, TPO, GSS, RAD52, ZNF580, and ZFP36L1 were significantly enriched in the "ROS metabolic process" term. As a transcription factor regulating genes enriched in the "ROS metabolic process" term, FOXQ1 could recruit LINC00887. Overexpression of FOXQ1 restored LINC00887 knockdown-induced downregulation of GSS and ZFP36L1 transcription in MTC cells. Additionally, FOXQ1 overexpression counteracted the inhibitory effects of LINC00887 knockdown on the proliferation and invasion of MTC cells and the promotion of intracellular ROS accumulation induced by LINC00887 knockdown. CONCLUSION: LINC00887 was identified as a key eRNA promoting the malignant phenotype of MTC cells. The involvement of FOXQ1 was essential for LINC00887 to play a pro-tumorigenic role in MTC. Our findings suggest that the FOXQ1/LINC00887 axis is a potential therapeutic target for MTC.

Shaping gene expression and its evolution by chromatin architecture and enhancer activity.

Transcriptional regulation plays a pivotal role in orchestrating the intricate genetic programs governing embryonic development. The expression of developmental genes relies on the combined activity of several cis-regulatory elements (CREs), such as enhancers and silencers, which can be located at long linear distances from the genes that they regulate and that interact with them through establishment of chromatin loops. Mutations affecting their activity or interaction with their target genes can lead to developmental disorders and are thought to have importantly contributed to the evolution of the animal body plan. The income of next-generation-sequencing approaches has allowed identifying over a million of sequences with putative regulatory potential in the human genome. Characterizing their function and establishing gene-CREs maps is essential to decode the logic governing developmental gene expression and is one of the major challenges of the post-genomic era. Chromatin 3D organization plays an essential role in determining how CREs specifically contact their target genes while avoiding deleterious off-target interactions. Our understanding of these aspects has greatly advanced with the income of chromatin conformation capture techniques and fluorescence microscopy approaches to visualize the organization of DNA elements in the nucleus. Here we will summarize relevant aspects of how the interplay between CRE activity and chromatin 3D organization regulates developmental gene expression and how it relates to pathological conditions and the evolution of animal body plan.

Identification and genetic validation of leukemia inhibitory factor super-enhancers in acute respiratory distress syndrome and lung cancer.

Super-enhancers play prominent roles in driving robust pathological gene expression, but they are hidden in human genome at noncoding regions, making them difficult to explore. Leukemia inhibitory factor (LIF) is a multifunctional cytokine crucially involved in acute respiratory distress syndrome (ARDS) and lung cancer progression. However, the mechanisms governing LIF regulation in disease contexts remain largely unexplored. In this study, we observed elevated levels of LIF in the bronchoalveolar lavage fluid (BALF) of patients with sepsis-related ARDS compared to those with nonsepsis-related ARDS. Furthermore, both basal and LPS-induced LIF expression were under the control of super-enhancers. Through analysis of H3K27Ac ChIP-seq data, we pinpointed three potential super-enhancers (LIF-SE1, LIF-SE2, and LIF-SE3) located proximal to the LIF gene in cells. Notably, genetic deletion of any of these three super-enhancers using CRISPR-Cas9 technology led to a significant reduction in LIF expression. Moreover, in cells lacking these super-enhancers, both cell growth and invasion capabilities were substantially impaired. Our findings highlight the critical role of three specific super-enhancers in regulating LIF expression and offer new insights into the transcriptional regulation of LIF in ARDS and lung cancer.

Enhancer selectivity across cell types delineates three functionally distinct enhancer-promoter regulation patterns.

BACKGROUND: Multiple enhancers co-regulating the same gene is prevalent and plays a crucial role during development and disease. However, how multiple enhancers coordinate the same gene expression across various cell types remains largely unexplored at genome scale. RESULTS: We develop a computational approach that enables the quantitative assessment of enhancer specificity and selectivity across diverse cell types, leveraging enhancer-promoter (E-P) interactions data. We observe two well-known gene regulation patterns controlled by enhancer clusters, which regulate the same gene either in a limited number of cell types (Specific pattern, Spe) or in the majority of cell types (Conserved pattern, Con), both of which are enriched for super-enhancers (SEs). We identify a previously overlooked pattern (Variable pattern, Var) that multiple enhancers link to the same gene, but rarely coexist in the same cell type. These three patterns control the genes associating with distinct biological function and exhibit unique epigenetic features. Specifically, we discover a subset of Var patterns contains Shared enhancers with stable enhancer-promoter interactions in the majority of cell types, which might contribute to maintaining gene expression by recruiting abundant CTCF. CONCLUSIONS: Together, our findings reveal three distinct E-P regulation patterns across different cell types, providing insights into deciphering the complexity of gene transcriptional regulation.

Copper-binding proteins and exonic splicing enhancers and silencers.

Eukaryotic DNA codes not only for proteins but contains a wealth of information required for accurate splicing of messenger RNA precursors and inclusion of constitutively or alternatively spliced exons in mature transcripts. This "auxiliary" splicing code has been characterized as exonic splicing enhancers and silencers (ESE and ESS). The exact interplay between protein and splicing codes is, however, poorly understood. Here, we show that exons encoding copper-coordinating amino acids in human cuproproteins lack ESEs and/or have an excess of ESSs, yet RNA sequencing and expressed sequence tags data show that they are more efficiently included in mature transcripts by the splicing machinery than average exons. Their largely constitutive inclusion in messenger RNA is facilitated by stronger splice sites, including polypyrimidine tracts, consistent with an important role of the surrounding intron architecture in ensuring high expression of metal-binding residues during evolution. ESE/ESS profiles of codons and entire exons that code for copper-coordinating residues were very similar to those encoding residues that coordinate zinc but markedly different from those that coordinate calcium. Together, these results reveal how the traditional and auxiliary splicing motifs responded to constraints of metal coordination in proteins.

A Bidirectional Non-Coding RNA Promoter Mediates Long-Range Gene Expression Regulation.

Recent evidence suggests that human gene promoters display gene expression regulatory mechanisms beyond the typical single gene local transcription modulation. In mammalian genomes, genes with an associated bidirectional promoter are abundant; bidirectional promoter architecture serves as a regulatory hub for a gene pair expression. However, it has been suggested that its contribution to transcriptional regulation might exceed local transcription initiation modulation. Despite their abundance, the functional consequences of bidirectional promoter architecture remain largely unexplored. This work studies the long-range gene expression regulatory role of a long non-coding RNA gene promoter using chromosome conformation capture methods. We found that this particular bidirectional promoter contributes to distal gene expression regulation in a target-specific manner by establishing promoter-promoter interactions. In particular, we validated that the promoter-promoter interactions of this regulatory element with the promoter of distal gene BBX contribute to modulating the transcription rate of this gene; removing the bidirectional promoter from its genomic context leads to a rearrangement of BBX promoter-enhancer interactions and to increased gene expression. Moreover, long-range regulatory functionality is not directly dependent on its associated non-coding gene pair expression levels.

Multiple cis-regulatory elements control prox1a expression in distinct lymphatic vascular beds.

During embryonic development, lymphatic endothelial cell (LEC) precursors are distinguished from blood endothelial cells by the expression of Prospero-related homeobox 1 (Prox1), which is essential for lymphatic vasculature formation in mouse and zebrafish. Prox1 expression initiation precedes LEC sprouting and migration, serving as the marker of specified LECs. Despite its crucial role in lymphatic development, Prox1 upstream regulation in LECs remains to be uncovered. SOX18 and COUP-TFII are thought to regulate Prox1 in mice by binding its promoter region. However, the specific regulation of Prox1 expression in LECs remains to be studied in detail. Here, we used evolutionary conservation and chromatin accessibility to identify enhancers located in the proximity of zebrafish prox1a active in developing LECs. We confirmed the functional role of the identified sequences through CRISPR/Cas9 mutagenesis of a lymphatic valve enhancer. The deletion of this region results in impaired valve morphology and function. Overall, our results reveal an intricate control of prox1a expression through a collection of enhancers. Ray-finned fish-specific distal enhancers drive pan-lymphatic expression, whereas vertebrate-conserved proximal enhancers refine expression in functionally distinct subsets of lymphatic endothelium.

Massively parallel characterization of regulatory elements in the developing human cortex.

Nucleotide changes in gene regulatory elements are important determinants of neuronal development and diseases. Using massively parallel reporter assays in primary human cells from mid-gestation cortex and cerebral organoids, we interrogated the cis-regulatory activity of 102,767 open chromatin regions, including thousands of sequences with cell type-specific accessibility and variants associated with brain gene regulation. In primary cells, we identified 46,802 active enhancer sequences and 164 variants that alter enhancer activity. Activity was comparable in organoids and primary cells, suggesting that organoids provide an adequate model for the developing cortex. Using deep learning we decoded the sequence basis and upstream regulators of enhancer activity. This work establishes a comprehensive catalog of functional gene regulatory elements and variants in human neuronal development.

Identification of drug responsive enhancers by predicting chromatin accessibility change from perturbed gene expression profiles.

Individual may response to drug treatment differently due to their genetic variants located in enhancers. These variants can alter transcription factor's (TF) binding strength, affect enhancer's chromatin activity or interaction, and eventually change expression level of downstream gene. Here, we propose a computational framework, PERD, to Predict the Enhancers Responsive to Drug. A machine learning model was trained to predict the genome-wide chromatin accessibility from transcriptome data using the paired expression and chromatin accessibility data collected from ENCODE and ROADMAP. Then the model was applied to the perturbed gene expression data from Connectivity Map (CMAP) and Cancer Drug-induced gene expression Signature DataBase (CDS-DB) and identify drug responsive enhancers with significantly altered chromatin accessibility. Furthermore, the drug responsive enhancers were related to the pharmacogenomics genome-wide association studies (PGx GWAS). Stepping on the traditional drug-associated gene signatures, PERD holds the promise to enhance the causality of drug perturbation by providing candidate regulatory element of those drug associated genes.

Pervasive structural heterogeneity rewires glioblastoma chromosomes to sustain patient-specific transcriptional programs.

Glioblastoma multiforme (GBM) encompasses brain malignancies marked by phenotypic and transcriptional heterogeneity thought to render these tumors aggressive, resistant to therapy, and inevitably recurrent. However, little is known about how the spatial organization of GBM genomes underlies this heterogeneity and its effects. Here, we compile a cohort of 28 patient-derived glioblastoma stem cell-like lines (GSCs) known to reflect the properties of their tumor-of-origin; six of these were primary-relapse tumor pairs from the same patient. We generate and analyze 5 kbp-resolution chromosome conformation capture (Hi-C) data from all GSCs to systematically map thousands of standalone and complex structural variants (SVs) and the multitude of neoloops arising as a result. By combining Hi-C, histone modification, and gene expression data with chromatin folding simulations, we explain how the pervasive, uneven, and idiosyncratic occurrence of neoloops sustains tumor-specific transcriptional programs via the formation of new enhancer-promoter contacts. We also show how even moderately recurrent neoloops can relate to patient-specific vulnerabilities. Together, our data provide a resource for dissecting GBM biology and heterogeneity, as well as for informing therapeutic approaches.

Regulation of Myc transcription by an enhancer cluster dedicated to pluripotency and early embryonic expression.

MYC plays various roles in pluripotent stem cells, including the promotion of somatic cell reprogramming to pluripotency, the regulation of cell competition and the control of embryonic diapause. However, how Myc expression is regulated in this context remains unknown. The Myc gene lies within a ~ 3-megabase gene desert with multiple cis-regulatory elements. Here we use genomic rearrangements, transgenesis and targeted mutation to analyse Myc regulation in early mouse embryos and pluripotent stem cells. We identify a topologically-associated region that homes enhancers dedicated to Myc transcriptional regulation in stem cells of the pre-implantation and early post-implantation embryo. Within this region, we identify elements exclusively dedicated to Myc regulation in pluripotent cells, with distinct enhancers that sequentially activate during naive and formative pluripotency. Deletion of pluripotency-specific enhancers dampens embryonic stem cell competitive ability. These results identify a topologically defined enhancer cluster dedicated to early embryonic expression and uncover a modular mechanism for the regulation of Myc expression in different states of pluripotency.

CUX1 regulates human hematopoietic stem cell chromatin accessibility via the BAF complex.

CUX1 is a homeodomain-containing transcription factor that is essential for the development and differentiation of multiple tissues. CUX1 is recurrently mutated or deleted in cancer, particularly in myeloid malignancies. However, the mechanism by which CUX1 regulates gene expression and differentiation remains poorly understood, creating a barrier to understanding the tumor-suppressive functions of CUX1. Here, we demonstrate that CUX1 directs the BAF chromatin remodeling complex to DNA to increase chromatin accessibility in hematopoietic cells. CUX1 preferentially regulates lineage-specific enhancers, and CUX1 target genes are predictive of cell fate in vivo. These data indicate that CUX1 regulates hematopoietic lineage commitment and homeostasis via pioneer factor activity, and CUX1 deficiency disrupts these processes in stem and progenitor cells, facilitating transformation.

STR mutations on chromosome 15q cause thyrotropin resistance by activating a primate-specific enhancer of MIR7-2/MIR1179.

Thyrotropin (TSH) is the master regulator of thyroid gland growth and function. Resistance to TSH (RTSH) describes conditions with reduced sensitivity to TSH. Dominantly inherited RTSH has been linked to a locus on chromosome 15q, but its genetic basis has remained elusive. Here we show that non-coding mutations in a (TTTG)4 short tandem repeat (STR) underlie dominantly inherited RTSH in all 82 affected participants from 12 unrelated families. The STR is contained in a primate-specific Alu retrotransposon with thyroid-specific cis-regulatory chromatin features. Fiber-seq and RNA-seq studies revealed that the mutant STR activates a thyroid-specific enhancer cluster, leading to haplotype-specific upregulation of the bicistronic MIR7-2/MIR1179 locus 35 kb downstream and overexpression of its microRNA products in the participants' thyrocytes. An imbalance in signaling pathways targeted by these micro-RNAs provides a working model for this cause of RTSH. This finding broadens our current knowledge of genetic defects altering pituitary-thyroid feedback regulation.

A CTCF-binding site in the Mdm1-Il22-Ifng locus shapes cytokine expression profiles and plays a critical role in early Th1 cell fate specification.

Cytokine expression during T cell differentiation is a highly regulated process that involves long-range promoter-enhancer and CTCF-CTCF contacts at cytokine loci. Here, we investigated the impact of dynamic chromatin loop formation within the topologically associating domain (TAD) in regulating the expression of interferon gamma (IFN-γ) and interleukin-22 (IL-22); these cytokine loci are closely located in the genome and are associated with complex enhancer landscapes, which are selectively active in type 1 and type 3 lymphocytes. In situ Hi-C analyses revealed inducible TADs that insulated Ifng and Il22 enhancers during Th1 cell differentiation. Targeted deletion of a 17 bp boundary motif of these TADs imbalanced Th1- and Th17-associated immunity, both in vitro and in vivo, upon Toxoplasma gondii infection. In contrast, this boundary element was dispensable for cytokine regulation in natural killer cells. Our findings suggest that precise cytokine regulation relies on lineage- and developmental stage-specific interactions of 3D chromatin architectures and enhancer landscapes.

A new suite of reporter vectors and a novel landing site survey system to study cis-regulatory elements in diverse insect species.

Comparative analyses between traditional model organisms, such as the fruit fly Drosophila melanogaster, and more recent model organisms, such as the red flour beetle Tribolium castaneum, have provided a wealth of insight into conserved and diverged aspects of gene regulation. While the study of trans-regulatory components is relatively straightforward, the study of cis-regulatory elements (CREs, or enhancers) remains challenging outside of Drosophila. A central component of this challenge has been finding a core promoter suitable for enhancer-reporter assays in diverse insect species. Previously, we demonstrated that a Drosophila Synthetic Core Promoter (DSCP) functions in a cross-species manner in Drosophila and Tribolium. Given the over 300 million years of divergence between the Diptera and Coleoptera, we reasoned that DSCP-based reporter constructs will be useful when studying cis-regulation in a variety of insect models across the holometabola and possibly beyond. To this end, we sought to create a suite of new DSCP-based reporter vectors, leveraging dual compatibility with piggyBac and PhiC31-integration, the 3xP3 universal eye marker, GATEWAY cloning, different colors of reporters and markers, as well as Gal4-UAS binary expression. While all constructs functioned properly with a Tc-nub enhancer in Drosophila, complications arose with tissue-specific Gal4-UAS binary expression in Tribolium. Nevertheless, the functionality of these constructs across multiple holometabolous orders suggests a high potential compatibility with a variety of other insects. In addition, we present the piggyLANDR (piggyBac-LoxP AttP Neutralizable Destination Reporter) platform for the establishment of proper PhiC31 landing sites free from position effects. As a proof-of-principle, we demonstrated the workflow for piggyLANDR in Drosophila. The potential utility of these tools ranges from molecular biology research to pest and disease-vector management, and will help advance the study of gene regulation beyond traditional insect models.

TAD border deletion at the Kit locus causes tissue-specific ectopic activation of a neighboring gene.

Topologically associated domains (TADs) restrict promoter-enhancer interactions, thereby maintaining the spatiotemporal pattern of gene activity. However, rearrangements of the TADs boundaries do not always lead to significant changes in the activity pattern. Here, we investigated the consequences of the TAD boundaries deletion on the expression of developmentally important genes encoding tyrosine kinase receptors: Kit, Kdr, Pdgfra. We used genome editing in mice to delete the TADs boundaries at the Kit locus and characterized chromatin folding and gene expression in pure cultures of fibroblasts, mast cells, and melanocytes. We found that although Kit is highly active in both mast cells and melanocytes, deletion of the TAD boundary between the Kit and Kdr genes results in ectopic activation only in melanocytes. Thus, the epigenetic landscape, namely the mutual arrangement of enhancers and actively transcribing genes, is important for predicting the consequences of the TAD boundaries removal. We also found that mice without a TAD border between the Kit and Kdr genes have a phenotypic manifestation of the mutation - a lighter coloration. Thus, the data obtained shed light on the principles of interaction between the 3D chromatin organization and epigenetic marks in the regulation of gene activity.

Estrogen receptor 1 chromatin profiling in human breast tumors reveals high inter-patient heterogeneity with enrichment of risk SNPs and enhancer activity at most-conserved regions.

Estrogen Receptor 1 (ESR1; also known as ERα, encoded by ESR1 gene) is the main driver and prime drug target in luminal breast cancer. ESR1 chromatin binding is extensively studied in cell lines and a limited number of human tumors, using consensi of peaks shared among samples. However, little is known about inter-tumor heterogeneity of ESR1 chromatin action, along with its biological implications. Here, we use a large set of ESR1 ChIP-seq data from 70 ESR1+ breast cancers to explore inter-patient heterogeneity in ESR1 DNA binding to reveal a striking inter-tumor heterogeneity of ESR1 action. Of note, commonly shared ESR1 sites show the highest estrogen-driven enhancer activity and are most engaged in long-range chromatin interactions. In addition, the most commonly shared ESR1-occupied enhancers are enriched for breast cancer risk SNP loci. We experimentally confirm SNVs to impact chromatin binding potential for ESR1 and its pioneer factor FOXA1. Finally, in the TCGA breast cancer cohort, we can confirm these variations to associate with differences in expression for the target gene. Cumulatively, we reveal a natural hierarchy of ESR1-chromatin interactions in breast cancers within a highly heterogeneous inter-tumor ESR1 landscape, with the most common shared regions being most active and affected by germline functional risk SNPs for breast cancer development.