Architectural and Functional Commonalities between Enhancers and Promoters

ASH2L Mediates Epidermal Differentiation and Hair Follicle Morphogenesis through H3K4me3 Modification

The processes of epidermal development in mammals are regulated by complex molecular mechanisms, such as histone modifications. Histone H3 lysine K4 methylation mediated by COMPASS (complex of proteins associated with Set1) methyltransferase is associated with gene activation, but its effect on epidermal lineage development remains unclear. Therefore, we constructed a mouse model of specific ASH2L (COMPASS methyltransferase core subunit) deletion in epidermal progenitor cells and investigated its effect on the development of mouse epidermal lineage. Furthermore, downstream target genes regulated by H3K4me3 were screened using RNA sequencing combined with Cleavage Under Targets and Tagmentation sequencing. Deletion of ASH2L in epidermal progenitor cells caused thinning of the suprabasal layer of the epidermis and delayed hair follicle morphogenesis in newborn mice. These phenotypes may be related to the reduced proliferative capacity of epidermal and hair follicle progenitor cells. ASH2L depletion may also lead to depletion of the epidermal stem cell pools in late mouse development. Finally, genes related to hair follicle development (Shh, Edar, and Fzd6), Notch signaling pathway (Notch2, Notch3, Hes5, and Nrarp), and ΔNp63 were identified as downstream target genes regulated by H3K4me3. Collectively, ASH2L-dependent H3K4me3 modification served as an upstream epigenetic regulator in epidermal differentiation and hair follicle morphogenesis in mice.

Transcriptional cascades during fasting amplify gluconeogenesis and instigate a secondary wave of ketogenic gene transcription

BACKGROUND AND AIMS

During fasting, bodily homeostasis is maintained due to hepatic production of glucose (gluconeogenesis) and ketone bodies (ketogenesis). The main hormones governing hepatic fuel production are glucagon and glucocorticoids that initiate transcriptional programs aimed at supporting gluconeogenesis and ketogenesis.

METHODS

Using primary mouse hepatocytes as an ex vivo model, we employed transcriptomic analysis (RNA-seq), genome-wide profiling of enhancer dynamics (ChIP-seq), perturbation experiments (inhibitors, shRNA), hepatic glucose production measurements and computational analyses.

RESULTS

We found that in addition to the known metabolic genes transcriptionally induced by glucagon and glucocorticoids, these hormones induce a set of genes encoding transcription factors (TFs) thereby initiating transcriptional cascades. Upon activation by glucocorticoids, the glucocorticoid receptor (GR) induced the genes encoding two TFs: CCAAT/enhancer-binding protein beta (C/EBPβ) and peroxisome proliferator-activated receptor alpha (PPARα). We found that the GR-C/EBPβ cascade mainly serves as a secondary amplifier of primary hormone-induced gene programs. C/EBPβ augmented gluconeogenic gene expression and hepatic glucose production. Conversely, the GR-PPARα cascade initiated a secondary transcriptional wave of genes supporting ketogenesis. The cascade led to synergistic induction of ketogenic genes which is dependent on protein synthesis. Genome-wide analysis of enhancer dynamics revealed numerous enhancers activated by the GR-PPARα cascade. These enhancers were proximal to ketogenic genes, enriched for the PPARα response element and showed increased PPARα binding.

CONCLUSION

This study reveals abundant transcriptional cascades occurring during fasting. These cascades serve two separated purposes: the amplification of the gluconeogenic transcriptional program and the induction of a gene program aimed at enhancing ketogenesis.

An RNA-centric view of transcription and genome organization

Foundational models of transcriptional regulation involve the assembly of protein complexes at DNA elements associated with specific genes. These assemblies, which can include transcription factors, cofactors, RNA polymerase, and various chromatin regulators, form dynamic spatial compartments that contribute to both gene regulation and local genome architecture. This DNA-protein-centric view has been modified with recent evidence that RNA molecules have important roles to play in gene regulation and genome structure. Here, we discuss evidence that gene regulation by RNA occurs at multiple levels that include assembly of transcriptional complexes and genome compartments, feedback regulation of active genes, silencing of genes, and control of protein kinases. We thus provide an RNA-centric view of transcriptional regulation that must reside alongside the more traditional DNA-protein-centric perspectives on gene regulation and genome architecture.

Dynamic chromatin architecture identifies new autoimmune-associated enhancers for IL2 and novel genes regulating CD4+ T cell activation

Genome-wide association studies (GWAS) have identified hundreds of genetic signals associated with autoimmune disease. The majority of these signals are located in non-coding regions and likely impact cis-regulatory elements (cRE). Because cRE function is dynamic across cell types and states, profiling the epigenetic status of cRE across physiological processes is necessary to characterize the molecular mechanisms by which autoimmune variants contribute to disease risk. We localized risk variants from 15 autoimmune GWAS to cRE active during TCR-CD28 co-stimulation of naïve human CD4+ T cells. To characterize how dynamic changes in gene expression correlate with cRE activity, we measured transcript levels, chromatin accessibility, and promoter-cRE contacts across three phases of naive CD4+ T cell activation using RNA-seq, ATAC-seq, and HiC. We identified ~1200 protein-coding genes physically connected to accessible disease-associated variants at 423 GWAS signals, at least one-third of which are dynamically regulated by activation. From these maps, we functionally validated a novel stretch of evolutionarily conserved intergenic enhancers whose activity is required for activation-induced IL2 gene expression in human and mouse, and is influenced by autoimmune-associated genetic variation. The set of genes implicated by this approach are enriched for genes controlling CD4+ T cell function and genes involved in human inborn errors of immunity, and we pharmacologically validated eight implicated genes as novel regulators of T cell activation. These studies directly show how autoimmune variants and the genes they regulate influence processes involved in CD4+ T cell proliferation and activation.

The role of circular RNA during the urological cancer metastasis: exploring regulatory mechanisms and potential therapeutic targets

Metastasis is a major contributor to treatment failure and death in urological cancers, representing an important biomedical challenge at present. Metastases form as a result of cancer cells leaving the primary site, entering the vasculature and lymphatic vessels, and colonizing clones elsewhere in the body. However, the specific regulatory mechanisms of action underlying the metastatic process of urological cancers remain incompletely elucidated. With the deepening of research, circular RNAs (circRNAs) have been found to not only play a significant role in tumor progression and prognosis but also show aberrant expression in various tumor metastases, consequently impacting tumor metastasis through multiple pathways. Therefore, circRNAs are emerging as potential tumor markers and treatment targets. This review summarizes the research progress on elucidating how circRNAs regulate the urological cancer invasion-metastasis cascade response and related processes, as well as their role in immune microenvironment remodeling and circRNA vaccines. This body of work highlights circRNA regulation as an emerging therapeutic target for urological cancers, which should motivate further specific research in this regard.

Enhancers associated with unstable RNAs are rare in plants

Unstable transcripts have emerged as markers of active enhancers in vertebrates and shown to be involved in many cellular processes and medical disorders. However, their prevalence and role in plants is largely unexplored. Here, we comprehensively captured all actively initiating (nascent) transcripts across diverse crops and other plants using capped small (cs)RNA sequencing. We discovered that unstable transcripts are rare in plants, unlike in vertebrates, and when present, often originate from promoters. In addition, many 'distal' elements in plants initiate tissue-specific stable transcripts and are likely bona fide promoters of as-yet-unannotated genes or non-coding RNAs, cautioning against using reference genome annotations to infer putative enhancer sites. To investigate enhancer function, we integrated data from self-transcribing active regulatory region (STARR) sequencing. We found that annotated promoters and other regions that initiate stable transcripts, but not those marked by unstable or bidirectional unstable transcripts, showed stronger enhancer activity in this assay. Our findings underscore the blurred line between promoters and enhancers and suggest that cis-regulatory elements can encompass diverse structures and mechanisms in eukaryotes, including humans.

DNA Transposons Favor De Novo Transcript Emergence Through Enrichment of Transcription Factor Binding Motifs

De novo genes emerge from noncoding regions of genomes via succession of mutations. Among others, such mutations activate transcription and create a new open reading frame (ORF). Although the mechanisms underlying ORF emergence are well documented, relatively little is known about the mechanisms enabling new transcription events. Yet, in many species a continuum between absent and very prominent transcription has been reported for essentially all regions of the genome. In this study, we searched for de novo transcripts by using newly assembled genomes and transcriptomes of seven inbred lines of Drosophila melanogaster, originating from six European and one African population. This setup allowed us to detect sample specific de novo transcripts, and compare them to their homologous nontranscribed regions in other samples, as well as genic and intergenic control sequences. We studied the association with transposable elements (TEs) and the enrichment of transcription factor motifs upstream of de novo emerged transcripts and compared them with regulatory elements. We found that de novo transcripts overlap with TEs more often than expected by chance. The emergence of new transcripts correlates with regions of high guanine-cytosine content and TE expression. Moreover, upstream regions of de novo transcripts are highly enriched with regulatory motifs. Such motifs are more enriched in new transcripts overlapping with TEs, particularly DNA TEs, and are more conserved upstream de novo transcripts than upstream their 'nontranscribed homologs'. Overall, our study demonstrates that TE insertion is important for transcript emergence, partly by introducing new regulatory motifs from DNA TE families.

An oncogenic enhancer promotes melanoma progression via regulating ETV4 expression

BACKGROUND

Enhancers are important gene regulatory elements that promote the expression of critical genes in development and disease. Aberrant enhancer can modulate cancer risk and activate oncogenes that lead to the occurrence of various cancers. However, the underlying mechanism of most enhancers in cancer remains unclear. Here, we aim to explore the function and mechanism of a crucial enhancer in melanoma.

METHODS

Multi-omics data were applied to identify an enhancer (enh17) involved in melanoma progression. To evaluate the function of enh17, CRISPR/Cas9 technology were applied to knockout enh17 in melanoma cell line A375. RNA-seq, ChIP-seq and Hi-C data analysis integrated with luciferase reporter assay were performed to identify the potential target gene of enh17. Functional experiments were conducted to further validate the function of the target gene ETV4. Multi-omics data integrated with CUT&Tag sequencing were performed to validate the binding profile of the inferred transcription factor STAT3.

RESULTS

An enhancer, named enh17 here, was found to be aberrantly activated and involved in melanoma progression. CRISPR/Cas9-mediated deletion of enh17 inhibited cell proliferation, migration, and tumor growth of melanoma both in vitro and in vivo. Mechanistically, we identified ETV4 as a target gene regulated by enh17, and functional experiments further support ETV4 as a target gene that is involved in cancer-associated phenotypes. In addition, STAT3 acts as a transcription factor binding with enh17 to regulate the transcription of ETV4.

CONCLUSIONS

Our findings revealed that enh17 plays an oncogenic role and promotes tumor progression in melanoma, and its transcriptional regulatory mechanisms were fully elucidated, which may open a promising window for melanoma prevention and treatment.

Inverted Alu repeats: friends or foes in the human transcriptome

Alu elements are highly abundant primate-specific short interspersed nuclear elements that account for ~10% of the human genome. Due to their preferential location in gene-rich regions, especially in introns and 3' UTRs, Alu elements can exert regulatory effects on the expression of both host and neighboring genes. When two Alu elements with inverse orientations are positioned in close proximity, their transcription results in the generation of distinct double-stranded RNAs (dsRNAs), known as inverted Alu repeats (IRAlus). IRAlus are key immunogenic self-dsRNAs and post-transcriptional cis-regulatory elements that play a role in circular RNA biogenesis, as well as RNA transport and stability. Recently, IRAlus dsRNAs have emerged as regulators of transcription and activators of Z-DNA-binding proteins. The formation and activity of IRAlus can be modulated through RNA editing and interactions with RNA-binding proteins, and misregulation of IRAlus has been implicated in several immune-associated disorders. In this review, we summarize the emerging functions of IRAlus dsRNAs, the regulatory mechanisms governing IRAlus activity, and their relevance in the pathogenesis of human diseases.

Dissecting the cis-regulatory syntax of transcription initiation with deep learning

Despite extensive characterization of mammalian Pol II transcription, the DNA sequence determinants of transcription initiation at a third of human promoters and most enhancers remain poorly understood. Hence, we trained and interpreted a neural network called ProCapNet that accurately models base-resolution initiation profiles from PRO-cap experiments using local DNA sequence. ProCapNet learns sequence motifs with distinct effects on initiation rates and TSS positioning and uncovers context-specific cryptic initiator elements intertwined within other TF motifs. ProCapNet annotates predictive motifs in nearly all actively transcribed regulatory elements across multiple cell-lines, revealing a shared cis-regulatory logic across promoters and enhancers mediated by a highly epistatic sequence syntax of cooperative and competitive motif interactions. ProCapNet models of RAMPAGE profiles measuring steady-state RNA abundance at TSSs distill initiation signals on par with models trained directly on PRO-cap profiles. ProCapNet learns a largely cell-type-agnostic cis-regulatory code of initiation complementing sequence drivers of cell-type-specific chromatin state critical for accurate prediction of cell-type-specific transcription initiation.

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.

Evolution of the enhancer-rich regulatory region of the gene for the cell-type specific transcription factor POU1F1

Precise spatio-temporal expression of genes in organogenesis is regulated by the coordinated interplay of DNA elements such as promoter and enhancers present in the regulatory region of a given locus. POU1F1 transcription factor plays a crucial role in the development of somatotrophs, lactotrophs and thyrotrophs in the anterior pituitary gland, and in maintaining high expression of growth hormone, prolactin and TSH. In mouse, expression of POU1F1 is controlled by a region fenced by two CTCF sites, containing 5 upstream enhancer elements, designated E-A (5' to 3'). Elements C, B and A correspond to elements shown previously to play a role in pituitary development and hormonal expression; functional roles for elements E and D have not been reported. We performed comparative sequence analysis of this regulatory region and discovered that three elements, B, C and E, are present in all vertebrate groups except Agnatha. One very long (>2 kb) element (A) is unique to mammals suggesting a specific change in regulation of the gene in this group. Using DNA accessibility assay (ATAC-seq) we showed that conserved elements in anterior pituitary of four non-mammals are open, suggesting functionality as regulatory elements. We showed that, in many non-mammalian vertebrates, an additional upstream exon closely follows element E, leading to alternatively spliced transcripts. Here, element E functions as an alternative promoter, but in mammals this feature is lost, suggesting conversion of alternative promoter to enhancer. Our work shows that regulation of POU1F1 changed markedly during the course of vertebrate evolution, use of a low number of enhancer elements combined with alternative promoters in non-mammalian vertebrates being replaced by use of a unique combination of regulatory units in mammals. Most importantly, our work suggests that evolutionary conversion of alternate promoter to enhancer could be one of the evolutionary mechanisms of enhancer birth.

Integrative approaches based on genomic techniques in the functional studies on enhancers

With the development of sequencing technology and the dramatic drop in sequencing cost, the functions of noncoding genes are being characterized in a wide variety of fields (e.g. biomedicine). Enhancers are noncoding DNA elements with vital transcription regulation functions. Tens of thousands of enhancers have been identified in the human genome; however, the location, function, target genes and regulatory mechanisms of most enhancers have not been elucidated thus far. As high-throughput sequencing techniques have leapt forwards, omics approaches have been extensively employed in enhancer research. Multidimensional genomic data integration enables the full exploration of the data and provides novel perspectives for screening, identification and characterization of the function and regulatory mechanisms of unknown enhancers. However, multidimensional genomic data are still difficult to integrate genome wide due to complex varieties, massive amounts, high rarity, etc. To facilitate the appropriate methods for studying enhancers with high efficacy, we delineate the principles, data processing modes and progress of various omics approaches to study enhancers and summarize the applications of traditional machine learning and deep learning in multi-omics integration in the enhancer field. In addition, the challenges encountered during the integration of multiple omics data are addressed. Overall, this review provides a comprehensive foundation for enhancer analysis.

Heterozygous Kmt2d loss diminishes enhancers to render medulloblastoma cells vulnerable to combinatory inhibition of lysine demethylation and oxidative phosphorylation

The histone H3 lysine 4 (H3K4) methyltransferase KMT2D (also called MLL4) is one of the most frequently mutated epigenetic modifiers in medulloblastoma (MB) and other types of cancer. Notably, heterozygous loss of KMT2D is prevalent in MB and other cancer types. However, what role heterozygous KMT2D loss plays in tumorigenesis has not been well characterized. Here, we show that heterozygous Kmt2d loss highly promotes MB driven by heterozygous loss of the MB suppressor gene Ptch in mice. Heterozygous Kmt2d loss upregulated tumor-promoting programs, including oxidative phosphorylation and G-protein-coupled receptor signaling, in Ptch-mutant-driven MB genesis. Mechanistically, both downregulation of the transcription-repressive tumor suppressor gene NCOR2 by heterozygous Kmt2d loss and upregulation of the oncogene MycN by heterozygous Ptch loss increased the expression of tumor-promoting genes. Moreover, heterozygous Kmt2d loss extensively diminished enhancer signals (e.g., H3K27ac) and H3K4me3 signature, including those for tumor suppressor genes (e.g., Ncor2). Combinatory pharmacological inhibition of oxidative phosphorylation and the H3K4 demethylase LSD1 drastically reduced tumorigenicity of MB cells bearing heterozygous Kmt2d loss. These findings reveal the mechanistic basis underlying the MB-promoting effect of heterozygous KMT2D loss, provide a rationale for a therapeutic strategy for treatment of KMT2D-deficient MB, and have mechanistic implications for the molecular pathogenesis of other types of cancer bearing heterozygous KMT2D loss.

Prognosis of colon cancer patients based on enhancer RNAs-related genes

PURPOSE

Colon cancer (CC) is a cancer of the large intestine with high prevalence and poor prognosis. enhancer RNAs. Therefore, valuable tools or biomarkers for predicting patient status, directing clinical practice, and reducing overtreatment are needed. Enhancer RNAs (eRNAs), a class of noncoding RNAs transcribed from enhancers, have been shown to function as regulators of oncogene or tumor suppressor gene expression. The aim of our study was to explore the potential roles of eRNAs and their target enhancer-related genes (ERGs) in the prognosis of CC.

METHODS

Selected CC cases (stage I-III) from The Cancer Genome Atlas database were used as a training set, and cases from the Gene Expression Omnibus were used as the validation set. ERGs associated with prognosis were screened through three steps: potential, candidate, and prognosis ERGs. Multivariate Cox proportional hazards analysis was used to identify independent prognostic factors, and a nomogram was created. Calibration curves were drawn by comparing predicted and observed survival probability. For validation, the calibration curves and ROC analysis were also applied to two external validation sets. The biological significance and clinical application of the genes obtained were investigated.

RESULTS

Based on the multiple tiers of strict screening, 11 prognostic ERGs were obtained, which were combined to obtain a prognosis signature. A compound nomogram integrating age, TNM classification, and the prognostic signature was constructed. The model was reliable in distinguishing the risk of patients with stage I-III CC, with AUCs of 0.78 and 0.70 at 5 and 7 years, respectively. There was good reproducibility in calibration curves. The prognostic model also yielded good prediction capability in the validation sets.

CONCLUSION

In this study, the usefulness and specificity of the ERGs in prognosis were described, which should be considered a key feature in the clinical guidance of CC patients with early stage. We concluded that the major implications of the eRNAs and ERGs should be valued, which would be an emerging hallmark in the prognosis of cancer.

Epromoters are new players in the regulatory landscape with potential pleiotropic roles

Precise spatiotemporal control of gene expression during normal development and cell differentiation is achieved by the combined action of proximal (promoters) and distal (enhancers) cis-regulatory elements. Recent studies have reported that a subset of promoters, termed Epromoters, works also as enhancers to regulate distal genes. This new paradigm opened novel questions regarding the complexity of our genome and raises the possibility that genetic variation within Epromoters has pleiotropic effects on various physiological and pathological traits by differentially impacting multiple proximal and distal genes. Here, we discuss the different observations pointing to an important role of Epromoters in the regulatory landscape and summarize the evidence supporting a pleiotropic impact of these elements in disease. We further hypothesize that Epromoter might represent a major contributor to phenotypic variation and disease.

Enhancers associated with unstable RNAs are rare in plants

Unstable transcripts have emerged as markers of active enhancers in vertebrates and shown to be involved in many cellular processes and medical disorders. However, their prevalence and role in plants is largely unexplored. Here, we comprehensively captured all actively initiating ("nascent") transcripts across diverse crops and other plants using capped small (cs)RNA-seq. We discovered that unstable transcripts are rare, unlike in vertebrates, and often originate from promoters. Additionally, many "distal" elements in plants initiate tissue-specific stable transcripts and are likely bone fide promoters of yet-unannotated genes or non-coding RNAs, cautioning against using genome annotations to infer "enhancers" or transcript stability. To investigate enhancer function, we integrated STARR-seq data. We found that annotated promoters, and other regions that initiate stable transcripts rather than unstable transcripts, function as stronger enhancers in plants. Our findings underscore the blurred line between promoters and enhancers and suggest that cis-regulatory elements encompass diverse structures and mechanisms in eukaryotes.

Retinal organoid and gene editing for basic and translational research

The rapid evolution of two technologies has greatly transformed the basic, translational, and clinical research in the mammalian retina. One is the retinal organoid (RO) technology. Various induction methods have been created or adapted to generate species-specific, disease-specific, and experimental-targeted retinal organoids (ROs). The process of generating ROs can highly mimic the in vivo retinal development, and consequently, the ROs resemble the retina in many aspects including the molecular and cellular profiles. The other technology is the gene editing, represented by the classical CRISPR-Cas9 editing and its derivatives such as prime editing, homology independent targeted integration (HITI), base editing and others. The combination of ROs and gene editing has opened up countless possibilities in the study of retinal development, pathogenesis, and therapeutics. We review recent advances in the ROs, gene editing methodologies, delivery vectors, and related topics that are particularly relevant to retinal studies.

Using Synthetic DNA Libraries to Investigate Chromatin and Gene Regulation

Despite the recent explosion in genome-wide studies in chromatin and gene regulation, we are still far from extracting a set of genetic rules that can predict the function of the regulatory genome. One major reason for this deficiency is that gene regulation is a multi-layered process that involves an enormous variable space, which cannot be fully explored using native genomes. This problem can be partially solved by introducing synthetic DNA libraries into cells, a method that can test the regulatory roles of thousands to millions of sequences with limited variables. Here, we review recent applications of this method to study transcription factor (TF) binding, nucleosome positioning, and transcriptional activity. We discuss the design principles, experimental procedures, and major findings from these studies and compare the pros and cons of different approaches.

High-sensitive nascent transcript sequencing reveals BRD4-specific control of widespread enhancer and target gene transcription

Gene transcription by RNA polymerase II (Pol II) is under control of promoters and distal regulatory elements known as enhancers. Enhancers are themselves transcribed by Pol II correlating with their activity. How enhancer transcription is regulated and coordinated with transcription at target genes has remained unclear. Here, we developed a high-sensitive native elongating transcript sequencing approach, called HiS-NET-seq, to provide an extended high-resolution view on transcription, especially at lowly transcribed regions such as enhancers. HiS-NET-seq uncovers new transcribed enhancers in human cells. A multi-omics analysis shows that genome-wide enhancer transcription depends on the BET family protein BRD4. Specifically, BRD4 co-localizes to enhancer and promoter-proximal gene regions, and is required for elongation activation at enhancers and their genes. BRD4 keeps a set of enhancers and genes in proximity through long-range contacts. From these studies BRD4 emerges as a general regulator of enhancer transcription that may link transcription at enhancers and genes.

3D organization of regulatory elements for transcriptional regulation in Arabidopsis

BACKGROUND

Although spatial organization of compartments and topologically associating domains at large scale is relatively well studied, the spatial organization of regulatory elements at fine scale is poorly understood in plants.

RESULTS

Here we perform high-resolution chromatin interaction analysis using paired-end tag sequencing approach. We map chromatin interactions tethered with RNA polymerase II and associated with heterochromatic, transcriptionally active, and Polycomb-repressive histone modifications in Arabidopsis. Analysis of the regulatory repertoire shows that distal active cis-regulatory elements are linked to their target genes through long-range chromatin interactions with increased expression of the target genes, while poised cis-regulatory elements are linked to their target genes through long-range chromatin interactions with depressed expression of the target genes. Furthermore, we demonstrate that transcription factor MYC2 is critical for chromatin spatial organization, and propose that MYC2 occupancy and MYC2-mediated chromatin interactions coordinately facilitate transcription within the framework of 3D chromatin architecture. Analysis of functionally related gene-defined chromatin connectivity networks reveals that genes implicated in flowering-time control are functionally compartmentalized into separate subdomains via their spatial activity in the leaf or shoot apical meristem, linking active mark- or Polycomb-repressive mark-associated chromatin conformation to coordinated gene expression.

CONCLUSION

The results reveal that the regulation of gene transcription in Arabidopsis is not only by linear juxtaposition, but also by long-range chromatin interactions. Our study uncovers the fine scale genome organization of Arabidopsis and the potential roles of such organization in orchestrating transcription and development.

Stochastic motion and transcriptional dynamics of pairs of distal DNA loci on a compacted chromosome

Chromosomes in the eukaryotic nucleus are highly compacted. However, for many functional processes, including transcription initiation, the pairwise motion of distal chromosomal elements such as enhancers and promoters is essential and necessitates dynamic fluidity. Here, we used a live-imaging assay to simultaneously measure the positions of pairs of enhancers and promoters and their transcriptional output while systematically varying the genomic separation between these two DNA loci. Our analysis reveals the coexistence of a compact globular organization and fast subdiffusive dynamics. These combined features cause an anomalous scaling of polymer relaxation times with genomic separation leading to long-ranged correlations. Thus, encounter times of DNA loci are much less dependent on genomic distance than predicted by existing polymer models, with potential consequences for eukaryotic gene expression.

Association between MGMT Enhancer Methylation and MGMT Promoter Methylation, MGMT Protein Expression, and Overall Survival in Glioblastoma

The repair protein O6-methylguanine-DNA methyltransferase (MGMT) is regulated epigenetically, mainly by the methylation of the MGMT promoter. MGMT promoter methylation status has emerged as a prognostic and predictive biomarker for patients with newly diagnosed glioblastoma (GBM). However, a strong negative correlation between MGMT promoter methylation and MGMT protein expression cannot be applied as a rule for all GBM patients. In order to investigate if the DNA methylation status of MGMT enhancers is associated with MGMT promoter methylation, MGMT expression, and the overall survival (OS) of GBM patients, we established assays based on high-resolution melting analysis and pyrosequencing for one intragenic and three intergenic MGMT enhancers. For CpGs in an enhancer located 560 kb upstream of the MGMT promoter, we found a significant negative correlation between the methylation status and MGMT protein levels of GBM samples expressing MGMT. The methylation status of CpGs in the intragenic enhancer (hs696) was strongly negatively correlated with MGMT promoter methylation and was significantly higher in MGMT-expressing GBM samples than in MGMT-non-expressing GBM samples. Moreover, low methylation of CpGs 01-03 and CpGs 09-13 was associated with the longer OS of the GBM patients. Our findings indicate an association between MGMT enhancer methylation and MGMT promoter methylation, MGMT protein expression, and/or OS.

ASH2L Controls Ureteric Bud Morphogenesis through the Regulation of RET/GFRA1 Signaling Activity in a Mouse Model

SIGNIFICANCE STATEMENT

Causes of congenital anomalies of the kidney and urinary tract (CAKUT) remain unclear. The authors investigated whether and how inactivation of Ash2l -which encodes a subunit of the COMPASS methyltransferase responsible for genome-wide histone H3 lysine K4 (H3K4) methylation-might contribute to CAKUT. In a mouse model, inactivation of Ash2l in the ureteric bud (UB) lineage led to CAKUT-like phenotypes. Removal of ASH2L led to deficient H3K4 trimethylation, which slowed cell proliferation at the UB tip, delaying budding and impairing branching morphogenesis. The absence of ASH2L also downregulated the expression of Ret , Gfra1 , and Wnt11 genes involved in RET/GFRA1 signaling. These findings identify ASH2L-mediated H3K4 methylation as an upstream epigenetic regulator of signaling crucial for UB morphogenesis and indicate that deficiency or dysregulation of these processes may lead to CAKUT.

BACKGROUND

Ureteric bud (UB) induction and branching morphogenesis are fundamental to the establishment of the renal architecture and are key determinants of nephron number. Defective UB morphogenesis could give rise to a spectrum of malformations associated with congenital anomalies of the kidney and urinary tract (CAKUT). Signaling involving glial cell line-derived neurotrophic factor and its receptor rearranged during transfection (RET) and coreceptor GFRA1 seems to be particularly important in UB development. Recent epigenome profiling studies have uncovered dynamic changes of histone H3 lysine K4 (H3K4) methylation during metanephros development, and dysregulated H3K4 methylation has been associated with a syndromic human CAKUT.

METHODS

To investigate whether and how inactivation of Ash2l , which encodes a subunit of the COMPASS methyltransferase responsible for genome-wide H3K4 methylation, might contribute to CAKUT, we inactivated Ash2l specifically from the UB lineage in C57BL/6 mice and examined the effects on genome-wide H3K4 methylation and metanephros development. Genes and epigenome changes potentially involved in these effects were screened using RNA-seq combined with Cleavage Under Targets and Tagmentation sequencing.

RESULTS

UB-specific inactivation of Ash2l caused CAKUT-like phenotypes mainly involving renal dysplasia at birth, which were associated with deficient H3K4 trimethylation. Ash2l inactivation slowed proliferation of cells at the UB tip, delaying budding and impairing UB branching morphogenesis. These effects were associated with downregulation of Ret , Gfra1 , and Wnt11 , which participate in RET/GFRA1 signaling.

CONCLUSIONS

These experiments identify ASH2L-dependent H3K4 methylation in the UB lineage as an upstream epigenetic regulator of RET/GFRA1 signaling in UB morphogenesis, which, if deficient, may lead to CAKUT.

An in Silico Approach to Identifying TF Binding Sites: Analysis of the Regulatory Regions of BUSCO Genes from Fungal Species in the Ceratocystidaceae Family

Transcriptional regulation controls gene expression through regulatory promoter regions that contain conserved sequence motifs. These motifs, also known as regulatory elements, are critically important to expression, which is driving research efforts to identify and characterize them. Yeasts have been the focus of such studies in fungi, including in several in silico approaches. This study aimed to determine whether in silico approaches could be used to identify motifs in the Ceratocystidaceae family, and if present, to evaluate whether these correspond to known transcription factors. This study targeted the 1000 base-pair region upstream of the start codon of 20 single-copy genes from the BUSCO dataset for motif discovery. Using the MEME and Tomtom analysis tools, conserved motifs at the family level were identified. The results show that such in silico approaches could identify known regulatory motifs in the Ceratocystidaceae and other unrelated species. This study provides support to ongoing efforts to use in silico analyses for motif discovery.

Enhancer profiling identifies epigenetic markers of endocrine resistance and reveals therapeutic options for metastatic castration-resistant prostate cancer patients

Androgen Receptor (AR) signaling inhibitors, including enzalutamide, are treatment options for patients with metastatic castration-resistant prostate cancer (mCRPC), but resistance inevitably develops. Using metastatic samples from a prospective phase II clinical trial, we epigenetically profiled enhancer/promoter activities with H3K27ac chromatin immunoprecipitation followed by sequencing, before and after AR-targeted therapy. We identified a distinct subset of H3K27ac-differentially marked regions that associated with treatment responsiveness. These data were successfully validated in mCRPC patient-derived xenograft models (PDX). In silico analyses revealed HDAC3 as a critical factor that can drive resistance to hormonal interventions, which we validated in vitro . Using cell lines and mCRPC PDX tumors in vitro , we identified drug-drug synergy between enzalutamide and the pan-HDAC inhibitor vorinostat, providing therapeutic proof-of-concept. These findings demonstrate rationale for new therapeutic strategies using a combination of AR and HDAC inhibitors to improve patient outcome in advanced stages of mCRPC.

Synthetic analysis of chromatin tracing and live-cell imaging indicates pervasive spatial coupling between genes

The role of the spatial organization of chromosomes in directing transcription remains an outstanding question in gene regulation. Here, we analyze two recent single-cell imaging methodologies applied across hundreds of genes to systematically analyze the contribution of chromosome conformation to transcriptional regulation. Those methodologies are (1) single-cell chromatin tracing with super-resolution imaging in fixed cells; and (2) high-throughput labeling and imaging of nascent RNA in living cells. Specifically, we determine the contribution of physical distance to the coordination of transcriptional bursts. We find that individual genes adopt a constrained conformation and reposition toward the centroid of the surrounding chromatin upon activation. Leveraging the variability in distance inherent in single-cell imaging, we show that physical distance - but not genomic distance - between genes on individual chromosomes is the major factor driving co-bursting. By combining this analysis with live-cell imaging, we arrive at a corrected transcriptional correlation of [Formula: see text] for genes separated by < 400 nm. We propose that this surprisingly large correlation represents a physical property of human chromosomes and establishes a benchmark for future experimental studies.

Stochastic motion and transcriptional dynamics of pairs of distal DNA loci on a compacted chromosome

Chromosomes in the eukaryotic nucleus are highly compacted. However, for many functional processes, including transcription initiation, the 3D pair-wise motion of distal chromosomal elements, such as enhancers and promoters, is essential and necessitates dynamic fluidity. Therefore, the interplay of chromosome organization and dynamics is crucial for gene regulation. Here, we use a live imaging assay to simultaneously measure the positions of pairs of enhancers and promoters and their transcriptional output in the developing fly embryo while systematically varying the genomic separation between these two DNA loci. Our analysis reveals a combination of a compact globular organization and fast subdiffusive dynamics. These combined features cause an anomalous scaling of polymer relaxation times with genomic separation and lead to long-ranged correlations compared to existing polymer models. This scaling implies that encounter times of DNA loci are much less dependent on genomic separation than predicted by existing polymer models, with potentially significant consequences for eukaryotic gene expression.

Complex Analysis of Single-Cell RNA Sequencing Data

Single-cell RNA sequencing (scRNA-seq) is a revolutionary tool for studying the physiology of normal and pathologically altered tissues. This approach provides information about molecular features (gene expression, mutations, chromatin accessibility, etc.) of cells, opens up the possibility to analyze the trajectories/phylogeny of cell differentiation and cell-cell interactions, and helps in discovery of new cell types and previously unexplored processes. From a clinical point of view, scRNA-seq facilitates deeper and more detailed analysis of molecular mechanisms of diseases and serves as a basis for the development of new preventive, diagnostic, and therapeutic strategies. The review describes different approaches to the analysis of scRNA-seq data, discusses the advantages and disadvantages of bioinformatics tools, provides recommendations and examples of their successful use, and suggests potential directions for improvement. We also emphasize the need for creating new protocols, including multiomics ones, for the preparation of DNA/RNA libraries of single cells with the purpose of more complete understanding of individual cells.

Plant Synthetic Promoters: Advancement and Prospective

Native/endogenous promoters have several fundamental limitations in terms of their size, Cis-elements distribution/patterning, and mode of induction, which is ultimately reflected in their insufficient transcriptional activity. Several customized synthetic promoters were designed and tested in plants during the past decade to circumvent such constraints. Such synthetic promoters have a built-in capacity to drive the expression of the foreign genes at their maximum amplitude in plant orthologous systems. The basic structure and function of the promoter has been discussed in this review, with emphasis on the role of the Cis-element in regulating gene expression. In addition to this, the necessity of synthetic promoters in the arena of plant biology has been highlighted. This review also provides explicit information on the two major approaches for developing plant-based synthetic promoters: the conventional approach (by utilizing the basic knowledge of promoter structure and Cis-trans interaction) and the advancement in gene editing technology. The success of plant genetic manipulation relies on the promoter efficiency and the expression level of the transgene. Therefore, advancements in the field of synthetic promoters has enormous potential in genetic engineering-mediated crop improvement.

TGFβ1-Induced EMT in the MCF10A Mammary Epithelial Cell Line Model Is Executed Independently of SNAIL1 and ZEB1 but Relies on JUNB-Coordinated Transcriptional Regulation

Epithelial-mesenchymal transition (EMT) fosters cancer cell invasion and metastasis, the main cause of cancer-related mortality. Growing evidence that SNAIL and ZEB transcription factors, typically portrayed as master regulators of EMT, may be dispensable for this process, led us to re-investigate its mechanistic underpinnings. For this, we used an unbiased computational approach that integrated time-resolved analyses of chromatin structure and differential gene expression, to predict transcriptional regulators of TGFβ1-inducible EMT in the MCF10A mammary epithelial cell line model. Bioinformatic analyses indicated comparatively minor contributions of SNAIL proteins and ZEB1 to TGFβ1-induced EMT, whereas the AP-1 subunit JUNB was anticipated to have a much larger impact. CRISPR/Cas9-mediated loss-of-function studies confirmed that TGFβ1-induced EMT proceeded independently of SNAIL proteins and ZEB1. In contrast, JUNB was necessary and sufficient for EMT in MCF10A cells, but not in A549 lung cancer cells, indicating cell-type-specificity of JUNB EMT-regulatory capacity. Nonetheless, the JUNB-dependence of EMT-associated transcriptional reprogramming in MCF10A cells allowed to define a gene expression signature which was regulated by TGFβ1 in diverse cellular backgrounds, showed positively correlated expression with TGFβ signaling in multiple cancer transcriptomes, and was predictive of patient survival in several cancer types. Altogether, our findings provide novel mechanistic insights into the context-dependent control of TGFβ1-driven EMT and thereby may lead to improved diagnostic and therapeutic options.

Efficient Selection of Enhancers and Promoters from MIA PaCa-2 Pancreatic Cancer Cells by ChIP-lentiMPRA

A library of active genome regulatory elements (putative promoters and enhancers) from MIA PaCa-2 pancreatic adenocarcinoma cells was constructed using a specially designed lentiviral vector and a massive parallel reporter assay (ChIP-lentiMPRA). Chromatin immunoprecipitation of the cell genomic DNA by H3K27ac antibodies was used for primary enrichment of the library for regulatory elements. Totally, 11,264 unique genome regions, many of which are capable of enhancing the expression of the CopGFP reporter gene from the minimal CMV promoter, were identified. The regions tend to be located near promoters. Based on the proximity assay, we found an enrichment of highly expressed genes among those associated with three or more mapped distal regions (2 kb distant from the 5'-ends of genes). It was shown significant enrichment of genes related to carcinogenesis or Mia PaCa-2 cell identity genes in this group. In contrast, genes associated with 1-2 distal regions or only with proximal regions (within 2 kbp of the 5'-ends of genes) are more often related to housekeeping functions. Thus, ChIP-lentiMPRA is a useful strategy for creating libraries of regulatory elements for the study of tumor-specific gene transcription.

A Myb enhancer-guided analysis of basophil and mast cell differentiation

The transcription factor MYB is a crucial regulator of hematopoietic stem and progenitor cells. However, the nature of lineage-specific enhancer usage of the Myb gene is largely unknown. We identify the Myb -68 enhancer, a regulatory element which marks basophils and mast cells. Using the Myb -68 enhancer activity, we show a population of granulocyte-macrophage progenitors with higher potential to differentiate into basophils and mast cells. Single cell RNA-seq demonstrates the differentiation trajectory is continuous from progenitors to mature basophils in vivo, characterizes bone marrow cells with a gene signature of mast cells, and identifies LILRB4 as a surface marker of basophil maturation. Together, our study leads to a better understanding of how MYB expression is regulated in a lineage-associated manner, and also shows how a combination of lineage-related reporter mice and single-cell transcriptomics can overcome the rarity of target cells and enhance our understanding of gene expression programs that control cell differentiation in vivo.

Cell fate decisions, transcription factors and signaling during early retinal development

The development of the vertebrate eyes is a complex process starting from anterior-posterior and dorso-ventral patterning of the anterior neural tube, resulting in the formation of the eye field. Symmetrical separation of the eye field at the anterior neural plate is followed by two symmetrical evaginations to generate a pair of optic vesicles. Next, reciprocal invagination of the optic vesicles with surface ectoderm-derived lens placodes generates double-layered optic cups. The inner and outer layers of the optic cups develop into the neural retina and retinal pigment epithelium (RPE), respectively. In vitro produced retinal tissues, called retinal organoids, are formed from human pluripotent stem cells, mimicking major steps of retinal differentiation in vivo. This review article summarizes recent progress in our understanding of early eye development, focusing on the formation the eye field, optic vesicles, and early optic cups. Recent single-cell transcriptomic studies are integrated with classical in vivo genetic and functional studies to uncover a range of cellular mechanisms underlying early eye development. The functions of signal transduction pathways and lineage-specific DNA-binding transcription factors are dissected to explain cell-specific regulatory mechanisms underlying cell fate determination during early eye development. The functions of homeodomain (HD) transcription factors Otx2, Pax6, Lhx2, Six3 and Six6, which are required for early eye development, are discussed in detail. Comprehensive understanding of the mechanisms of early eye development provides insight into the molecular and cellular basis of developmental ocular anomalies, such as optic cup coloboma. Lastly, modeling human development and inherited retinal diseases using stem cell-derived retinal organoids generates opportunities to discover novel therapies for retinal diseases.

An Expanded Interplay Network between NF-κB p65 (RelA) and E2F1 Transcription Factors: Roles in Physiology and Pathology

Transcription Factors (TFs) are the main regulators of gene expression, controlling among others cell homeostasis, identity, and fate. TFs may either act synergistically or antagonistically on nearby regulatory elements and their interplay may activate or repress gene expression. The family of NF-κB TFs is among the most important TFs in the regulation of inflammation, immunity, and stress-like responses, while they also control cell growth and survival, and are involved in inflammatory diseases and cancer. The family of E2F TFs are major regulators of cell cycle progression in most cell types. Several studies have suggested the interplay between these two TFs in the regulation of numerous genes controlling several biological processes. In the present study, we compared the genomic binding landscape of NF-κB RelA/p65 subunit and E2F1 TFs, based on high throughput ChIP-seq and RNA-seq data in different cell types. We confirmed that RelA/p65 has a binding profile with a high preference for distal enhancers bearing active chromatin marks which is distinct to that of E2F1, which mostly generates promoter-specific binding. Moreover, the RelA/p65 subunit and E2F1 cistromes have limited overlap and tend to bind chromatin that is in an active state even prior to immunogenic stimulation. Finally, we found that a fraction of the E2F1 cistrome is recruited by NF-κΒ near pro-inflammatory genes following LPS stimulation in immune cell types.

Characterization of sequence determinants of enhancer function using natural genetic variation

Sequence variation in enhancers that control cell-type-specific gene transcription contributes significantly to phenotypic variation within human populations. However, it remains difficult to predict precisely the effect of any given sequence variant on enhancer function due to the complexity of DNA sequence motifs that determine transcription factor (TF) binding to enhancers in their native genomic context. Using F₁-hybrid cells derived from crosses between distantly related inbred strains of mice, we identified thousands of enhancers with allele-specific TF binding and/or activity. We find that genetic variants located within the central region of enhancers are most likely to alter TF binding and enhancer activity. We observe that the AP-1 family of TFs (Fos/Jun) are frequently required for binding of TEAD TFs and for enhancer function. However, many sequence variants outside of core motifs for AP-1 and TEAD also impact enhancer function, including sequences flanking core TF motifs and AP-1 half sites. Taken together, these data represent one of the most comprehensive assessments of allele-specific TF binding and enhancer function to date and reveal how sequence changes at enhancers alter their function across evolutionary timescales.

A multi-omic dissection of super-enhancer driven oncogenic gene expression programs in ovarian cancer

The human genome contains regulatory elements, such as enhancers, that are often rewired by cancer cells for the activation of genes that promote tumorigenesis and resistance to therapy. This is especially true for cancers that have little or no known driver mutations within protein coding genes, such as ovarian cancer. Herein, we utilize an integrated set of genomic and epigenomic datasets to identify clinically relevant super-enhancers that are preferentially amplified in ovarian cancer patients. We systematically probe the top 86 super-enhancers, using CRISPR-interference and CRISPR-deletion assays coupled to RNA-sequencing, to nominate two salient super-enhancers that drive proliferation and migration of cancer cells. Utilizing Hi-C, we construct chromatin interaction maps that enable the annotation of direct target genes for these super-enhancers and confirm their activity specifically within the cancer cell compartment of human tumors using single-cell genomics data. Together, our multi-omic approach examines a number of fundamental questions about how regulatory information encoded into super-enhancers drives gene expression networks that underlie the biology of ovarian cancer.

Dynamic enhancer transcription associates with reprogramming of immune genes during pattern triggered immunity in Arabidopsis

BACKGROUND

Enhancers are cis-regulatory elements present in eukaryote genomes, which constitute indispensable determinants of gene regulation by governing the spatiotemporal and quantitative expression dynamics of target genes, and are involved in multiple life processes, for instance during development and disease states. The importance of enhancer activity has additionally been highlighted for immune responses in animals and plants; however, the dynamics of enhancer activities and molecular functions in plant innate immunity are largely unknown. Here, we investigated the involvement of distal enhancers in early innate immunity in Arabidopsis thaliana.

RESULTS

A group of putative distal enhancers producing low-abundance transcripts either unidirectionally or bidirectionally are identified. We show that enhancer transcripts are dynamically modulated in plant immunity triggered by microbe-associated molecular patterns and are strongly correlated with open chromatin, low levels of methylated DNA, and increases in RNA polymerase II targeting and acetylated histone marks. Dynamic enhancer transcription is correlated with target early immune gene expression patterns. Cis motifs that are bound by immune-related transcription factors, such as WRKYs and SARD1, are highly enriched within upregulated enhancers. Moreover, a subset of core pattern-induced enhancers are upregulated by multiple patterns from diverse pathogens. The expression dynamics of putative immunity-related enhancers and the importance of WRKY binding motifs for enhancer function were also validated.

CONCLUSIONS

Our study demonstrates the general occurrence of enhancer transcription in plants and provides novel information on the distal regulatory landscape during early plant innate immunity, providing new insights into immune gene regulation and ultimately improving the mechanistic understanding of the plant immune system.

Analysis of the landscape of human enhancer sequences in biological databases

The process of gene regulation extends as a network in which both genetic sequences and proteins are involved. The levels of regulation and the mechanisms involved are multiple. Transcription is the main control mechanism for most genes, being the downstream steps responsible for refining the transcription patterns. In turn, gene transcription is mainly controlled by regulatory events that occur at promoters and enhancers. Several studies are focused on analyzing the contribution of enhancers in the development of diseases and their possible use as therapeutic targets. The study of regulatory elements has advanced rapidly in recent years with the development and use of next generation sequencing techniques. All this information has generated a large volume of information that has been transferred to a growing number of public repositories that store this information. In this article, we analyze the content of those public repositories that contain information about human enhancers with the aim of detecting whether the knowledge generated by scientific research is contained in those databases in a way that could be computationally exploited. The analysis will be based on three main aspects identified in the literature: types of enhancers, type of evidence about the enhancers, and methods for detecting enhancer-promoter interactions. Our results show that no single database facilitates the optimal exploitation of enhancer data, most types of enhancers are not represented in the databases and there is need for a standardized model for enhancers. We have identified major gaps and challenges for the computational exploitation of enhancer data.

Enhancer RNAs stimulate Pol II pause release by harnessing multivalent interactions to NELF

Enhancer RNAs (eRNAs) are long non-coding RNAs that originate from enhancers. Although eRNA transcription is a canonical feature of activated enhancers, the molecular features required for eRNA function and the mechanism of how eRNAs impinge on target gene transcription have not been established. Thus, using eRNA-dependent RNA polymerase II (Pol II) pause release as a model, we here investigate the requirement of sequence, structure and length of eRNAs for their ability to stimulate Pol II pause release by detaching NELF from paused Pol II. We find eRNAs not to exert their function through common structural or sequence motifs. Instead, eRNAs that exhibit a length >200 nucleotides and that contain unpaired guanosines make multiple, allosteric contacts with NELF subunits -A and -E to trigger efficient NELF release. By revealing the molecular determinants of eRNA function, our study establishes eRNAs as an important player in Pol II pause release, and it provides new insight into the regulation of metazoan transcription.

Skeletal and gene-regulatory functions of nuclear sex steroid hormone receptors

The wide variety of sex hormone actions underlie bone growth and health, and their actions mediate gene regulation by the cognate nuclear receptors. Nuclear androgen and estrogen receptors (AR, and ERα/ERβ) are hormone-dependent and DNA binding- transcription regulatory factors, and gene regulation by sex hormones often accompany with chromatin remodeling under aid of a number of co-regulators. As sex hormone biosynthesis is under highly regulated systemic and local regulations, the skeletal actions of sex hormones could be inferred from only the phenotypic abnormalities in skeleton in mouse genetic models deficient of nuclear receptors selectively in specific types of bone cells as well as at specific cell differentiation stages. Anabolic androgen actions and anti-bone resorptive estrogen actions are discussed here from the phenotypic abnormalities in such model mice. Though rapid gene regulation by sex hormones may not require chromatin reorganization, dynamic chromatin reconfiguration looks to facilitate profound and long-term hormonal actions. In this review, we focus the recent findings in gene regulation at a chromatin level, particularly of the function of enhancer RNAs transcribed from strong enhancers, and in the role of liquid-liquid phase separation state in transcription initiation through chromatin reconfiguration.

Enhancer transcription detected in the nascent transcriptomic landscape of bread wheat

The precise spatiotemporal gene expression is orchestrated by enhancers that lack general sequence features and thus are difficult to be computationally identified. By nascent RNA sequencing combined with epigenome profiling, we detect active transcription of enhancers from the complex bread wheat genome. We find that genes associated with transcriptional enhancers are expressed at significantly higher levels, and enhancer RNA is more precise and robust in predicting enhancer activity compared to chromatin features. We demonstrate that sub-genome-biased enhancer transcription could drive sub-genome-biased gene expression. This study highlights enhancer transcription as a hallmark in regulating gene expression in wheat.

Decoding Human Genome Regulatory Features That Influence HIV-1 Proviral Expression and Fate Through an Integrated Genomics Approach

Fundamental principles of HIV-1 integration into the human genome have been revealed in the past 2 decades. However, the impact of the integration site on proviral transcription and expression remains poorly understood. Solving this problem requires the analysis of multiple genomic datasets for thousands of proviral integration sites. Here, we generated and combined large-scale datasets, including epigenetics, transcriptome, and 3-dimensional genome architecture to interrogate the chromatin states, transcription activity, and nuclear sub-compartments around HIV-1 integrations in Jurkat CD4⁺ T cells to decipher human genome regulatory features shaping the transcription of proviral classes based on their position and orientation in the genome. Through a Hidden Markov Model and ranked informative values prior to a machine learning logistic regression model, we defined nuclear sub-compartments and chromatin states contributing to genomic architecture, transcriptional activity, and nucleosome density of regions neighboring the integration site, as additive features influencing HIV-1 expression. Our integrated genomics approach also allows for a robust experimental design, in which HIV-1 can be genetically introduced into precise genomic locations with known regulatory features to assess the relationship of integration positions to viral transcription and fate.

Delineation of the DNA Structural Features of Eukaryotic Core Promoter Classes

The eukaryotic transcription is orchestrated from a chunk of the DNA region stated as the core promoter. Multifarious and punctilious core promoter signals, viz., TATA-box, Inr, BREs, and Pause Button, are associated with a subset of genes and regulate their spatiotemporal expression. However, the core promoter architecture linked with these signals has not been investigated exhaustively for several species. In this study, we attempted to envisage the adaptive binding landscape of the transcription initiation machinery as a function of DNA structure. To this end, we deployed a set of k-mer based DNA structural estimates and regular expression models derived from experiments, molecular dynamic simulations, and theoretical frameworks, and high-throughout promoter data sets retrieved from the eukaryotic promoter database. We categorized protein-coding gene core promoters based on characteristic motifs at precise locations and analyzed the B-DNA structural properties and non-B-DNA structural motifs for 15 different eukaryotic genomes. We observed that Inr, BREd, and no-motif classes display common patterns of DNA sequence and structural environment. TATA-containing, BREu, and Pause Button classes show a deviant behavior with the TATA class displaying varied axial and twisting flexibility while BREu and Pause Button leaned toward G-quadruplex motif enrichment. Intriguingly, DNA meltability and shape signals are conserved irrespective of the presence or absence of distinct core promoter motifs in the majority of species. Altogether, here we delineated the conserved DNA structural signals associated with several promoter classes that may contribute to the chromatin configuration, orchestration of transcription machinery, and DNA duplex melting during the transcription process.

Interplay between nuclear factor-κB, p38 MAPK and glucocorticoid receptor signaling synergistically induces functional TLR2 in lung epithelial cells

While glucocorticoids act via the glucocorticoid receptor (GR; NR3C1) to reduce the expression of many inflammatory genes, repression is not an invariable outcome. Here, we explore synergy occurring between synthetic glucocorticoids (dexamethasone and budesonide) and proinflammatory cytokines (IL1B and TNF) on the expression of the toll-like receptor 2 (TLR2). This effect is observed in epithelial cell lines and both undifferentiated and differentiated primary human bronchial epithelial cells (pHBECs). In A549 cells, IL1B-plus-glucocorticoid–induced TLR2 expression required nuclear factor (NF)-κB and GR. Likewise, in A549 cells, BEAS-2B cells, and pHBECs, chromatin immunoprecipitation identified GR- and NF-κB/p65-binding regions ∼32 kb (R1) and ∼7.3 kb (R2) upstream of the TLR2 gene. Treatment of BEAS-2B cells with TNF or/and dexamethasone followed by global run-on sequencing confirmed transcriptional activity at these regions. Furthermore, cloning R1 or R2 into luciferase reporters revealed transcriptional activation by budesonide or IL1B, respectively, while R1+R2 juxtaposition enabled synergistic activation by IL1B and budesonide. In addition, small-molecule inhibitors and siRNA knockdown showed p38α MAPK to negatively regulate both IL1B-induced TLR2 expression and R1+R2 reporter activity. Finally, agonism of IL1B-plus-dexamethasone–induced TLR2 in A549 cells and pHBECs stimulated NF-κB- and interferon regulatory factor-dependent reporter activity and chemokine release. We conclude that glucocorticoid-plus-cytokine-driven synergy at TLR2 involves GR and NF-κB acting via specific enhancer regions, which combined with the inhibition of p38α MAPK promotes TLR2 expression. Subsequent inflammatory effects that occur following TLR2 agonism may be pertinent in severe neutrophilic asthma or chronic obstructive pulmonary disease, where glucocorticoid-based therapies are less efficacious.

Single-cell-resolved dynamics of chromatin architecture delineate cell and regulatory states in zebrafish embryos

DNA accessibility of cis-regulatory elements (CREs) dictates transcriptional activity and drives cell differentiation during development. While many genes regulating embryonic development have been identified, the underlying CRE dynamics controlling their expression remain largely uncharacterized. To address this, we produced a multimodal resource and genomic regulatory map for the zebrafish community, which integrates single-cell combinatorial indexing assay for transposase-accessible chromatin with high-throughput sequencing (sci-ATAC-seq) with bulk histone PTMs and Hi-C data to achieve a genome-wide classification of the regulatory architecture determining transcriptional activity in the 24-h post-fertilization (hpf) embryo. We characterized the genome-wide chromatin architecture at bulk and single-cell resolution, applying sci-ATAC-seq on whole 24-hpf stage zebrafish embryos, generating accessibility profiles for ∼23,000 single nuclei. We developed a genome segmentation method, ScregSeg (single-cell regulatory landscape segmentation), for defining regulatory programs, and candidate CREs, specific to one or more cell types. We integrated the ScregSeg output with bulk measurements for histone post-translational modifications and 3D genome organization and identified new regulatory principles between chromatin modalities prevalent during zebrafish development. Sci-ATAC-seq profiling of npas4l/cloche mutant embryos identified novel cellular roles for this hematovascular transcriptional master regulator and suggests an intricate mechanism regulating its expression. Our work defines regulatory architecture and principles in the zebrafish embryo and establishes a resource of cell-type-specific genome-wide regulatory annotations and candidate CREs, providing a valuable open resource for genomics, developmental, molecular, and computational biology.

Exploring the role of eRNA in regulating gene expression

eRNAs as the products of enhancers can regulate gene expression via various possible ways, but which regulation way is more reasonable is debatable in biology, and in particular, how eRNAs impact gene expression remains unclear. Here we introduce a mechanistic model of gene expression to address these issues. This model considers three possible regulation ways of eRNA: Type-I by which eRNA regulates transcriptional activity by facilitating the formation of enhancer-promoter (E-P) loop, Type-II by which eRNA directly promotes the mRNA production rate, and mixed regulation (i.e., the combination of Type-I and Type-II). We show that with the increase of the E-P loop length, mRNA distribution can transition from unimodality to bimodality or vice versa in all the three regulation cases. However, in contrast to the other two regulations, Type-II regulation can lead to the highest mean mRNA level and the lowest mRNA noise, independent of the E-P loop length. These results would not only reveal the essential mechanism of how eRNA regulates gene expression, but also imply a new mechanism for phenotypic switching, namely the E-P loop can induce phenotypic switching.

Constructing gene regulatory networks using epigenetic data

The biological processes that drive cellular function can be represented by a complex network of interactions between regulators (transcription factors) and their targets (genes). A cell's epigenetic state plays an important role in mediating these interactions, primarily by influencing chromatin accessibility. However, how to effectively use epigenetic data when constructing a gene regulatory network remains an open question. Almost all existing network reconstruction approaches focus on estimating transcription factor to gene connections using transcriptomic data. In contrast, computational approaches for analyzing epigenetic data generally focus on improving transcription factor binding site predictions rather than deducing regulatory network relationships. We bridged this gap by developing SPIDER, a network reconstruction approach that incorporates epigenetic data into a message-passing framework to estimate gene regulatory networks. We validated SPIDER's predictions using ChIP-seq data from ENCODE and found that SPIDER networks are both highly accurate and include cell-line-specific regulatory interactions. Notably, SPIDER can recover ChIP-seq verified transcription factor binding events in the regulatory regions of genes that do not have a corresponding sequence motif. The networks estimated by SPIDER have the potential to identify novel hypotheses that will allow us to better characterize cell-type and phenotype specific regulatory mechanisms.