Transcriptional enhancers: from properties to genome-wide predictions

High-resolution dissection of human cell type-specific enhancers in cis and trans activities

The spatiotemporal-specific gene expression is regulated by cell type-specific regulatory elements. Here we selected the H3K4me1-associated DNA sequences as candidate enhancers in two different human cell lines and performed ChIP-STARR-seq to quantify the cell-type-specific enhancer activities with high-resolution. We investigated how the activity landscape of enhancers would change when transferred from native cells (cis activity) to another cell lines (trans activity). We obtained enhancers cis activity maps and trans activity maps in two different cell lines. The cis and trans activity maps enabled us to identify cell type-specific active enhancers, with enrichment of motifs of differentially expressed TFs. Comparisons between the cis and trans activity maps revealed general consistent regulatory property with different levels of activity in two cell lines, suggesting sequence intrinsic regulatory properties remain similar in different types of cells. This study provides a new perspective on sequence intrinsic enhancer activities in different types of cells.

Twenty-first century mouse genetics is again at an inflection point

The laboratory mouse has been the premier model organism for biomedical research owing to the availability of multiple well-characterized inbred strains, its mammalian physiology and its homozygous genome, and because experiments can be performed under conditions that control environmental variables. Moreover, its genome can be genetically modified to assess the impact of allelic variation on phenotype. Mouse models have been used to discover or test many therapies that are commonly used today. Mouse genetic discoveries are often made using genome-wide association study methods that compare allelic differences in panels of inbred mouse strains with their phenotypic responses. Here we examine changes in the methods used to analyze mouse genetic models of biomedical traits during the twenty-first century. To do this, we first examine where mouse genetics was before the first inflection point, which was just before the revolution in genome sequencing that occurred 20 years ago, and then describe the factors that have accelerated the pace of mouse genetic discovery. We focus on mouse genetic studies that have generated findings that either were translated to humans or could impact clinical medicine or drug development. We next explore how advances in computational capabilities and in DNA sequencing methodology during the past 20 years could enhance the ability of mouse genetics to produce solutions for twenty-first century public-health problems.

A novel geminivirus-derived 3' flanking sequence of terminator mediates the gene expression enhancement

Exploring the new elements to re-design the expression cassette is crucial in synthetic biology. Viruses are one of the most important sources for exploring gene expression elements. In this study, we found that the DNA sequence of the SBG51 deltasatellite from the Sweet potato leaf curl virus (SPLCV) greatly enhanced the gene expression when flanked downstream of the terminator. The SBG51 sequence increased transient GFP gene expression in Nicotiana benthamiana leaves by up to ~6 times and ~10 times compared to the gene expression controlled by the UBQ10 promoter and 35S promoter alone, respectively. The increased GFP gene expression level contributed to the continuous accumulation of GFP protein and GFP fluorescence until 8 days post-inoculation (dpi). The SBG51 sequence also enhanced the gene expression in the transgenic Arabidopsis plants and maintained the spatio-temporal pattern of the FLOWERING LOCUS T (FT) and TOO MANY MOUTHS (TMM) promoters. We identified a 123 bp of AT-rich sequence containing seven "ATAAA" or "TTAAA" elements from the SBG51 DNA, which had the gene expression enhancement effect. Furthermore, the artificial synthetic sequences containing tandem repeated "ATAAA" or "TTAAA" elements were sufficient to increase the gene expression but did not alter the polyadenylation of mRNA, similar to the function of matrix attachment regions (MAR). Additionally, the compact artificial synthetic sequence also had an effect on yeast when the expression cassette was integrated into the genome. We conclude that the geminivirus deltasatellite-derived sequence and the "ATAAA"/"TTAAA" elements are powerful tools for enhancing gene expression.

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

BACKGROUND

STARR-seq and other massively-parallel reporter assays are widely used to discover functional enhancers in transfected cell models, which can be confounded by plasmid vector-induced type-I interferon immune responses and lack the multicellular environment and endogenous chromatin state of complex mammalian tissues.

RESULTS

We describe HDI-STARR-seq, which combines STARR-seq plasmid library delivery to the liver, by hydrodynamic tail vein injection (HDI), with reporter RNA transcriptional initiation driven by a minimal Albumin promoter, which we show is essential for mouse liver STARR-seq enhancer activity assayed 7 days after HDI. Importantly, little or no vector-induced innate type-I interferon responses were observed. Comparisons of HDI-STARR-seq activity between male and female mouse livers and in livers from males treated with an activating ligand of the transcription factor (TF) CAR (Nr1i3) identified many condition-dependent enhancers linked to condition-specific gene expression. Further, thousands of active liver enhancers were identified using a high complexity STARR-seq library comprised of ~ 50,000 genomic regions released by DNase-I digestion of mouse liver nuclei. When compared to stringently inactive library sequences, the active enhancer sequences identified were highly enriched for liver open chromatin regions with activating histone marks (H3K27ac, H3K4me1, H3K4me3), were significantly closer to gene transcriptional start sites, and were significantly depleted of repressive (H3K27me3, H3K9me3) and transcribed region histone marks (H3K36me3).

CONCLUSION

HDI-STARR-seq offers substantial improvements over current methodologies for large scale, functional profiling of enhancers, including condition-dependent enhancers, in liver tissue in vivo, and can be adapted to characterize enhancer activities in a variety of species and tissues by selecting suitable tissue- and species-specific promoter sequences.

LDB1 establishes multi-enhancer networks to regulate gene expression

How specific enhancer-promoter pairing is established remains mostly unclear. Besides the CTCF/cohesin machinery, few nuclear factors have been studied for a direct role in physically connecting regulatory elements. Using a murine erythroid cell model, we show via acute degradation experiments that LDB1 directly and broadly promotes connectivity among regulatory elements. Most LDB1-mediated contacts, even those spanning hundreds of kb, can form in the absence of CTCF, cohesin, or YY1 as determined using multiple degron systems. Moreover, an engineered LDB1-driven chromatin loop is cohesin independent. Cohesin-driven loop extrusion does not stall at LDB1-occupied sites but aids the formation of a subset of LDB1-anchored loops. Leveraging the dynamic reorganization of nuclear architecture during the transition from mitosis to G1 phase, we observe that loop formation and de novo LDB1 occupancy correlate and can occur independently of structural loops. Tri-C and Region Capture Micro-C reveal that LDB1 organizes multi-enhancer networks to activate transcription. These findings establish LDB1 as a driver of spatial connectivity.

Identification, characterization, and design of plant genome sequences using deep learning

Due to its excellent performance in processing large amounts of data and capturing complex non-linear relationships, deep learning has been widely applied in many fields of plant biology. Here we first review the application of deep learning in analyzing genome sequences to predict gene expression, chromatin interactions, and epigenetic features (open chromatin, transcription factor binding sites, and methylation sites) in plants. Then, current motif mining and functional component design and synthesis based on generative adversarial networks, large models, and attention mechanisms are elaborated in detail. The progress of protein structure and function prediction, genomic prediction, and large model applications based on deep learning is also discussed. Finally, this work provides prospects for the future development of deep learning in plants with regard to multiple omics data, algorithm optimization, large language models, sequence design, and intelligent breeding.

A genome-wide screen identifies silencers with distinct chromatin properties and mechanisms of repression

Differential gene transcription enables development and homeostasis in all animals and is regulated by two major classes of distal cis-regulatory DNA elements (CREs): enhancers and silencers. Although enhancers have been thoroughly characterized, the properties and mechanisms of silencers remain largely unknown. By an unbiased genome-wide functional screen in Drosophila melanogaster S2 cells, we discover a class of silencers that bind one of three transcription factors (TFs) and are generally not included in chromatin-defined CRE catalogs as they mostly lack detectable DNA accessibility. The silencer-binding TF CG11247, which we term Saft, safeguards cell fate decisions in vivo and functions via a highly conserved domain we term zinc-finger-associated C-terminal (ZAC) and the corepressor G9a, independently of G9a's H3K9-methyltransferase activity. Overall, our identification of silencers with unexpected properties and mechanisms has important implications for the understanding and future study of repressive CREs, as well as the functional annotation of animal genomes.

Benchmarking tools for transcription factor prioritization

Spatiotemporal regulation of gene expression is controlled by transcription factor (TF) binding to regulatory elements, resulting in a plethora of cell types and cell states from the same genetic information. Due to the importance of regulatory elements, various sequencing methods have been developed to localise them in genomes, for example using ChIP-seq profiling of the histone mark H3K27ac that marks active regulatory regions. Moreover, multiple tools have been developed to predict TF binding to these regulatory elements based on DNA sequence. As altered gene expression is a hallmark of disease phenotypes, identifying TFs driving such gene expression programs is critical for the identification of novel drug targets. In this study, we curated 84 chromatin profiling experiments (H3K27ac ChIP-seq) where TFs were perturbed through e.g., genetic knockout or overexpression. We ran nine published tools to prioritize TFs using these real-world datasets and evaluated the performance of the methods in identifying the perturbed TFs. This allowed the nomination of three frontrunner tools, namely RcisTarget, MEIRLOP and monaLisa. Our analyses revealed opportunities and commonalities of tools that will help to guide further improvements and developments in the field.

Synthetic Promoters in Gene Therapy: Design Approaches, Features and Applications

Gene therapy is a promising approach to the treatment of various inherited diseases, but its development is complicated by a number of limitations of the natural promoters used. The currently used strong ubiquitous natural promoters do not allow for the specificity of expression, while natural tissue-specific promoters have lowactivity. These limitations of natural promoters can be addressed by creating new synthetic promoters that achieve high levels of tissue-specific target gene expression. This review discusses recent advances in the development of synthetic promoters that provide a more precise regulation of gene expression. Approaches to the design of synthetic promoters are reviewed, including manual design and bioinformatic methods using machine learning. Examples of successful applications of synthetic promoters in the therapy of hereditary diseases and cancer are presented, as well as prospects for their clinical use.

Integrative pan-cancer analysis reveals a common architecture of dysregulated transcriptional networks characterized by loss of enhancer methylation

Aberrant DNA methylation contributes to gene expression deregulation in cancer. However, these alterations' precise regulatory role and clinical implications are still not fully understood. In this study, we performed expression-methylation Quantitative Trait Loci (emQTL) analysis to identify deregulated cancer-driving transcriptional networks linked to CpG demethylation pan-cancer. By analyzing 33 cancer types from The Cancer Genome Atlas, we identified and confirmed significant correlations between CpG methylation and gene expression (emQTL) in cis and trans, both across and within cancer types. Bipartite network analysis of the emQTL revealed groups of CpGs and genes related to important biological processes involved in carcinogenesis including proliferation, metabolism and hormone-signaling. These bipartite communities were characterized by loss of enhancer methylation in specific transcription factor binding regions (TFBRs) and the CpGs were topologically linked to upregulated genes through chromatin loops. Penalized Cox regression analysis showed a significant prognostic impact of the pan-cancer emQTL in many cancer types. Taken together, our integrative pan-cancer analysis reveals a common architecture where hallmark cancer-driving functions are affected by the loss of enhancer methylation and may be epigenetically regulated.

The Evolution of Complex Multicellularity in Land Plants

The evolution of complex multicellularity in land plants represents a pivotal event in the history of life on Earth, characterized by significant increases in biological complexity. This transition, classified as a Major Evolutionary Transition (MET), is best understood through the framework of Evolutionary Transitions in Individuality (ETIs), which focuses on formerly independent entities forming higher-level units that lose their reproductive autonomy. While much of the ETI literature has concentrated on the early stages of multicellularity, such as the formation and maintenance stages, this paper seeks to address the less explored transformation stage. To do so, we apply an approach that we call Transitions in Structural Complexity (TSCs), which focuses on the emergence of new units of organization via the three key evolutionary processes of modularization, subfunctionalization, and integration to the evolution of land plants. To lay the groundwork, we first explore the relationships between sex, individuality, and units of selection to highlight a sexual life cycle-based perspective on ETIs by examining the early stages of the transition to multicellularity (formation) in the sexual life cycle of the unicellular common ancestor of land plants, emphasizing the differences between the transition to multicellularity in eumetazoans and land plants. We then directly apply the TSC approach in this group, identifying key evolutionary events such as the distinct evolutionary innovations like shoot, root, vascular systems, and specialized reproductive structures, arguing that bringing these under the broader rubric of TSCs affords a degree of explanatory unification. By examining these evolutionary processes, this paper provides a new perspective on the evolution of multicellularity in land plants, highlighting both parallels and distinctions with the animal kingdom.

Large-scale analysis of the integration of enhancer-enhancer signals by promoters

Genes are often regulated by multiple enhancers. It is poorly understood how the individual enhancer activities are combined to control promoter activity. Anecdotal evidence has shown that enhancers can combine sub-additively, additively, synergistically, or redundantly. However, it is not clear which of these modes are more frequent in mammalian genomes. Here, we systematically tested how pairs of enhancers activate promoters using a three-way combinatorial reporter assay in mouse embryonic stem cells. By assaying about 69,000 enhancer-enhancer-promoter combinations we found that enhancer pairs generally combine near-additively. This behaviour was conserved across seven developmental promoters tested. Surprisingly, these promoters scale the enhancer signals in a non-linear manner that depends on promoter strength. A housekeeping promoter showed an overall different response to enhancer pairs, and a smaller dynamic range. Thus, our data indicate that enhancers mostly act additively, but promoters transform their collective effect non-linearly.

Decoding resilience: ecology, regulation, and evolution of biosynthetic gene clusters

Secondary metabolism is crucial for plant survival and can generate chemistry with nutritional, therapeutic, and industrial value. Biosynthetic genes of selected secondary metabolites cluster within localised chromosomal regions. The arrangement of these biosynthetic gene clusters (BGCs) challenges the long-held model of random gene order in eukaryotes, raising questions about their regulation, ecological significance, and evolution. In this review, we address these questions by exploring the contribution of BGCs to ecologically relevant plant-biotic interactions, while also evaluating the molecular-(epi)genetic mechanisms controlling their coordinated stress- and tissue-specific expression. Based on evidence that BGCs have distinct chromatin signatures and are enriched with transposable elements (TEs), we integrate emerging hypotheses into an updated evolutionary model emphasising how stress-induced epigenetic processes have shaped BGC formation.

Predicting genome-wide tissue-specific enhancers via combinatorial transcription factor genomic occupancy analysis

Enhancers are non-coding cis-regulatory elements crucial for transcriptional regulation. Mutations in enhancers can disrupt gene regulation, leading to disease phenotypes. Identifying enhancers and their tissue-specific activity is challenging due to their lack of stereotyped sequences. This study presents a sequence-based computational model that uses combinatorial transcription factor (TF) genomic occupancy to predict tissue-specific enhancers. Trained on diverse datasets, including ENCODE and Vista enhancer browser data, the model predicted 25 000 forebrain-specific cis-regulatory modules (CRMs) in the human genome. Validation using biochemical features, disease-associated SNPs, and in vivo zebrafish analysis confirmed its effectiveness. This model aids in predicting enhancers lacking well-characterized chromatin features, complementing experimental approaches in tissue-specific enhancer discovery.

Genes with differential expression across ancestries are enriched in ancestry-specific disease effects likely due to gene-by-environment interactions

Multi-ancestry genome-wide association studies (GWASs) have highlighted the existence of variants with ancestry-specific effect sizes. Understanding where and why these ancestry-specific effects occur is fundamental to understanding the genetic basis of human diseases and complex traits. Here, we characterized genes differentially expressed across ancestries (ancDE genes) at the cell-type level by leveraging single-cell RNA-sequencing data in peripheral blood mononuclear cells for 21 individuals with East Asian (EAS) ancestry and 23 individuals with European (EUR) ancestry (172,385 cells); then, we tested whether variants surrounding those genes were enriched in disease variants with ancestry-specific effect sizes by leveraging ancestry-matched GWASs of 31 diseases and complex traits (average n ∼ 90,000 and ∼ 267,000 in EAS and EUR, respectively). We observed that ancDE genes tended to be cell-type specific and enriched in genes interacting with the environment and in variants with ancestry-specific disease effect sizes, which suggests cell-type-specific, gene-by-environment interactions shared between regulatory and disease architectures. Finally, we illustrated how different environments might have led to ancestry-specific myeloid cell leukemia 1 (MCL1) expression in B cells and ancestry-specific allele effect sizes in lymphocyte count GWASs for variants surrounding MCL1. Our results imply that large single-cell and GWAS datasets from diverse ancestries are required to improve our understanding of human diseases.

Integration of multi-omics summary data reveals the role of N6-methyladenosine in neuropsychiatric disorders

N6-methyladenosine (m6A) methylation regulates gene expression/protein by influencing numerous aspects of mRNA metabolism and contributes to neuropsychiatric diseases. Here, we integrated multi-omics data and genome-wide association study summary data of schizophrenia (SCZ), bipolar disorder (BP), attention deficit hyperactivity disorder (ADHD), autism spectrum disorder (ASD), major depressive disorder (MDD), Alzheimer's disease (AD), and Parkinson's disease (PD) to reveal the role of m6A in neuropsychiatric disorders by using transcriptome-wide association study (TWAS) tool and Summary-data-based Mendelian randomization (SMR). Our investigation identified 86 m6A sites associated with seven neuropsychiatric diseases and then revealed 7881 associations between m6A sites and gene expressions. Based on these results, we discovered 916 significant m6A-gene associations involving 82 disease-related m6A sites and 606 genes. Further integrating the 58 disease-related genes from TWAS and SMR analysis, we obtained 61, 8, 7, 3, and 2 associations linking m6A-disease, m6A-gene, and gene-disease for SCZ, BP, AD, MDD, and PD separately. Functional analysis showed the m6A mapped genes were enriched in "response to stimulus" pathway. In addition, we also analyzed the effect of gene expression on m6A and the post-transcription effect of m6A on protein. Our study provided new insights into the genetic component of m6A in neuropsychiatric disorders and unveiled potential pathogenic mechanisms where m6A exerts influences on disease through gene expression/protein regulation.

UVB-induced TRPS1 regulates MITF transcription activity to promote skin pigmentation

Hyperpigmented dermatoses are characterized by increased skin pigmentation caused by genetic, environmental factors and inflammation, which lasts a long time and is difficult to treat. Ultraviolet (UV), especially ultraviolet B (UVB), is the primary external factor inducing skin pigmentation. However, the specific regulatory mechanisms are not fully understood. Through analysis of GEO datasets from four UV-exposed skin cell/tissue samples, we found that TRPS1 is the only gene differentially expressed in multiple datasets (GSE22083, GSE67098 and GSE70280) and highly positively correlated with the expression of key melanogenesis genes. Consistently, we observed that TRPS1 is highly expressed in sun-exposed skin tissues compared to non-exposed skin. Additionally, the expression of TRPS1 was also significantly upregulated after UVB irradiation in isolated skin tissues and melanocytes, while knockdown of TRPS1 expression inhibited the UVB-induced melanogenesis. Further research revealed that overexpression of TRPS1 increased melanin content and tyrosinase activity in MNT1 cells, as well as upregulated the expression levels of key melanogenesis genes (MITF, TYR, TYRP1, DCT). In contrast, inhibition of TRPS1 expression showed the opposite effect. Moreover, we found that TRPS1 can bind to the promoter region of MITF, inhibiting the expression of MITF can antagonize the melanogenesis induced by TRPS1. In conclusion, UVB-induced TRPS1 promotes melanogenesis by activating the transcriptional activity of MITF.

Evaluating the representational power of pre-trained DNA language models for regulatory genomics

The emergence of genomic language models (gLMs) offers an unsupervised approach to learning a wide diversity of cis-regulatory patterns in the non-coding genome without requiring labels of functional activity generated by wet-lab experiments. Previous evaluations have shown that pre-trained gLMs can be leveraged to improve predictive performance across a broad range of regulatory genomics tasks, albeit using relatively simple benchmark datasets and baseline models. Since the gLMs in these studies were tested upon fine-tuning their weights for each downstream task, determining whether gLM representations embody a foundational understanding of cis-regulatory biology remains an open question. Here we evaluate the representational power of pre-trained gLMs to predict and interpret cell-type-specific functional genomics data that span DNA and RNA regulation. Our findings suggest that probing the representations of pre-trained gLMs do not offer substantial advantages over conventional machine learning approaches that use one-hot encoded sequences. This work highlights a major gap with current gLMs, raising potential issues in conventional pre-training strategies for the non-coding genome.

A comparative analysis of planarian genomes reveals regulatory conservation in the face of rapid structural divergence

The planarian Schmidtea mediterranea is being studied as a model species for regeneration, but the assembly of planarian genomes remains challenging. Here, we report a high-quality haplotype-phased, chromosome-scale genome assembly of the sexual S2 strain of S. mediterranea and high-quality chromosome-scale assemblies of its three close relatives, S. polychroa, S. nova, and S. lugubris. Using hybrid gene annotations and optimized ATAC-seq and ChIP-seq protocols for regulatory element annotation, we provide valuable genome resources for the planarian research community and a first comparative perspective on planarian genome evolution. Our analyses reveal substantial divergence in protein-coding sequences and regulatory regions but considerable conservation within promoter and enhancer annotations. We also find frequent retrotransposon-associated chromosomal inversions and interchromosomal translocations within the genus Schmidtea and, remarkably, independent and nearly complete losses of ancestral metazoan synteny in Schmidtea and two other flatworm groups. Overall, our results suggest that platyhelminth genomes can evolve without syntenic constraints.

Oncogenic Enhancers in Leukemia

Although the study of leukemogenesis has traditionally focused on protein-coding genes, the role of enhancer dysregulation is becoming increasingly recognized. The advent of high-throughput sequencing, together with a better understanding of enhancer biology, has revealed how various genetic and epigenetic lesions produce oncogenic enhancers that drive transformation. These aberrations include translocations that lead to enhancer hijacking, point mutations that modulate enhancer activity, and copy number alterations that modify enhancer dosage. In this review, we describe these mechanisms in the context of leukemia and discuss potential therapeutic avenues to target these regulatory elements. Significance: Large-scale sequencing projects have uncovered recurrent gene mutations in leukemia, but the picture remains incomplete: some patients harbor no such aberrations, whereas others carry only a few that are insufficient to bring about transformation on their own. One of the missing pieces is enhancer dysfunction, which only recently has emerged as a critical driver of leukemogenesis. Knowledge of the various mechanisms of enhancer dysregulation is thus key for a complete understanding of leukemia and its causes, as well as the development of targeted therapies in the era of precision medicine.

W2V-repeated index: Prediction of enhancers and their strength based on repeated fragments

Enhancers are crucial in gene expression regulation, dictating the specificity and timing of transcriptional activity, which highlights the importance of their identification for unravelling the intricacies of genetic regulation. Therefore, it is critical to identify enhancers and their strengths. Repeated sequences in the genome are repeats of the same or symmetrical fragments. There has been a great deal of evidence that repetitive sequences contain enormous amounts of genetic information. Thus, We introduce the W2V-Repeated Index, designed to identify enhancer sequence fragments and evaluates their strength through the analysis of repeated K-mer sequences in enhancer regions. Utilizing the word2vector algorithm for numerical conversion and Manta Ray Foraging Optimization for feature selection, this method effectively captures the frequency and distribution of K-mer sequences. By concentrating on repeated K-mer sequences, it minimizes computational complexity and facilitates the analysis of larger K values. Experiments indicate that our method performs better than all other advanced methods on almost all indicators.

Relaxin Modulates the Genomic Actions and Biological Effects of Estrogen in the Myometrium

Estradiol (E2) and relaxin (Rln) are steroid and polypeptide hormones, respectively, with important roles in the female reproductive tract, including myometrium. Some actions of Rln, which are mediated by its membrane receptor RXFP1, require or are augmented by E2 signaling through its cognate nuclear steroid receptor, estrogen receptor alpha (ERα). In contrast, other actions of Rln act in opposition to the effects of E2. Here we explored the molecular and genomic mechanisms that underlie the functional interplay between E2 and Rln in the myometrium. We used both ovariectomized female mice and immortalized human myometrial cells expressing wild-type or mutant ERα (hTERT-HM-ERα cells). Our results indicate that Rln modulates the genomic actions and biological effects of estrogen in the myometrium and myometrial cells by reducing phosphorylation of ERα on serine 118 (S118), as well as by reducing the E2-dependent binding of ERα across the genome. These effects were associated with changes in the hormone-regulated transcriptome, including a decrease in the E2-dependent expression of some genes and enhanced expression of others. The inhibitory effects of Rln cotreatment on the E2-dependent phosphorylation of ERα required the nuclear dual-specificity phosphatases DUSP1 and DUSP5. Moreover, the inhibitory effects of Rln were reflected in a concomitant inhibition of the E2-dependent contraction of myometrial cells. Collectively, our results identify a pathway that integrates Rln/RXFP1 and E2/ERα signaling, resulting in a convergence of membrane and nuclear signaling pathways to control genomic and biological outcomes.

LDB1 establishes multi-enhancer networks to regulate gene expression

How specific enhancer-promoter pairing is established is still mostly unclear. Besides the CTCF/cohesin machinery, only a few nuclear factors have been studied for a direct role in physically connecting regulatory elements. Here, we show via acute degradation experiments that LDB1 directly and broadly promotes enhancer-promoter loops. Most LDB1-mediated contacts, even those spanning hundreds of kb, can form in the absence of CTCF, cohesin, or YY1 as determined via the use of multiple degron systems. Moreover, an engineered LDB1-driven chromatin loop is cohesin independent. Cohesin-driven loop extrusion does not stall at LDB1 occupied sites but may aid the formation of a subset of LDB1 anchored loops. Leveraging the dynamic reorganization of nuclear architecture during the transition from mitosis to G1-phase, we establish a relationship between LDB1-dependent interactions in the context of TAD organization and gene activation. Lastly, Tri-C and Region Capture Micro-C reveal that LDB1 organizes multi-enhancer networks to activate transcription. This establishes LDB1 as a direct driver of regulatory network inter-connectivity.

Artificially inserted strong promoter containing multiple G-quadruplexes induces long-range chromatin modification

Although the role of G-quadruplex (G4) DNA structures has been suggested in chromosomal looping this was not tested directly. Here, to test causal function, an array of G4s, or control sequence that does not form G4s, were inserted within chromatin in cells. In vivo G4 formation of the inserted G4 sequence array, and not the control sequence, was confirmed using G4-selective antibody. Compared to the control insert, we observed a remarkable increase in the number of 3D chromatin looping interactions from the inserted G4 array. This was evident within the immediate topologically associated domain (TAD) and throughout the genome. Locally, recruitment of enhancer histone marks and the transcriptional coactivator p300/Acetylated-p300 increased in the G4-array, but not in the control insertion. Resulting promoter-enhancer interactions and gene activation were clear up to 5 Mb away from the insertion site. Together, these show the causal role of G4s in enhancer function and long-range chromatin interactions. Mechanisms of 3D topology are primarily based on DNA-bound architectural proteins that induce/stabilize long-range interactions. Involvement of the underlying intrinsic DNA sequence/structure in 3D looping shown here therefore throws new light on how long-range chromosomal interactions might be induced or maintained.

Quantitative analysis of cis-regulatory elements in transcription with KAS-ATAC-seq

Cis-regulatory elements (CREs) are pivotal in orchestrating gene expression throughout diverse biological systems. Accurate identification and in-depth characterization of functional CREs are crucial for decoding gene regulation networks during cellular processes. In this study, we develop Kethoxal-Assisted Single-stranded DNA Assay for Transposase-Accessible Chromatin with Sequencing (KAS-ATAC-seq) to quantitatively analyze the transcriptional activity of CREs. A main advantage of KAS-ATAC-seq lies in its precise measurement of ssDNA levels within both proximal and distal ATAC-seq peaks, enabling the identification of transcriptional regulatory sequences. This feature is particularly adept at defining Single-Stranded Transcribing Enhancers (SSTEs). SSTEs are highly enriched with nascent RNAs and specific transcription factors (TFs) binding sites that define cellular identity. Moreover, KAS-ATAC-seq provides a detailed characterization and functional implications of various SSTE subtypes. Our analysis of CREs during mouse neural differentiation demonstrates that KAS-ATAC-seq can effectively identify immediate-early activated CREs in response to retinoic acid (RA) treatment. Our findings indicate that KAS-ATAC-seq provides more precise annotation of functional CREs in transcription. Future applications of KAS-ATAC-seq would help elucidate the intricate dynamics of gene regulation in diverse biological processes.

BertSNR: an interpretable deep learning framework for single-nucleotide resolution identification of transcription factor binding sites based on DNA language model

MOTIVATION

Transcription factors are pivotal in the regulation of gene expression, and accurate identification of transcription factor binding sites (TFBSs) at high resolution is crucial for understanding the mechanisms underlying gene regulation. The task of identifying TFBSs from DNA sequences is a significant challenge in the field of computational biology today. To address this challenge, a variety of computational approaches have been developed. However, these methods face limitations in their ability to achieve high-resolution identification and often lack interpretability.

RESULTS

We propose BertSNR, an interpretable deep learning framework for identifying TFBSs at single-nucleotide resolution. BertSNR integrates sequence-level and token-level information by multi-task learning based on pre-trained DNA language models. Benchmarking comparisons show that our BertSNR outperforms the existing state-of-the-art methods in TFBS predictions. Importantly, we enhanced the interpretability of the model through attentional weight visualization and motif analysis, and discovered the subtle relationship between attention weight and motif. Moreover, BertSNR effectively identifies TFBSs in promoter regions, facilitating the study of intricate gene regulation.

AVAILABILITY AND IMPLEMENTATION

The BertSNR source code can be found at https://github.com/lhy0322/BertSNR.

TRIM24 Cooperates with Ras Mutation to Drive Glioma Progression through snoRNA Recruitment of PHAX and DNA-PKcs

The factors driving glioma progression remain poorly understood. Here, the epigenetic regulator TRIM24 is identified as a driver of glioma progression, where TRIM24 overexpression promotes HRasV12 anaplastic astrocytoma (AA) progression into epithelioid GBM (Ep-GBM)-like tumors. Co-transfection of TRIM24 with HRasV12 also induces Ep-GBM-like transformation of human neural stem cells (hNSCs) with tumor protein p53 gene (TP53) knockdown. Furthermore, TRIM24 is highly expressed in clinical Ep-GBM specimens. Using single-cell RNA-sequencing (scRNA-Seq), the authors show that TRIM24 overexpression impacts both intratumoral heterogeneity and the tumor microenvironment. Mechanically, HRasV12 activates phosphorylated adaptor for RNA export (PHAX) and upregulates U3 small nucleolar RNAs (U3 snoRNAs) to recruit Ku-dependent DNA-dependent protein kinase catalytic subunit (DNA-PKcs). Overexpressed TRIM24 is also recruited by PHAX to U3 snoRNAs, thereby facilitating DNA-PKcs phosphorylation of TRIM24 at S767/768 residues. Phosphorylated TRIM24 induces epigenome and transcription factor network reprogramming and promotes Ep-GBM-like transformation. Targeting DNA-PKcs with the small molecule inhibitor NU7441 synergizes with temozolomide to reduce Ep-GBM tumorigenicity and prolong animal survival. These findings provide new insights into the epigenetic regulation of Ep-GBM-like transformation and suggest a potential therapeutic strategy for patients with Ep-GBM.

Single-Cell Chromatin Accessibility Analysis Reveals Subgroup-Specific TF-NTR Regulatory Circuits in Medulloblastoma

Medulloblastoma (MB) stands as one of the prevalent malignant brain tumors among pediatric patients. Despite its prevalence, the intricate interplay between the regulatory program driving malignancy in MB cells and their interactions with the microenvironment remains insufficiently understood. Leveraging the capabilities of single-cell Assay for Transposase-Accessible Chromatin sequencing (scATAC-seq), the chromatin accessibility landscape is unveiled across 59,015 distinct MB cells. This expansive dataset encompasses cells belonging to discrete molecular subgroups, namely SHH, WNT, Group3, and Group4. Within these chromatin accessibility profiles, specific regulatory elements tied to individual subgroups are uncovered, shedding light on the distinct activities of transcription factors (TFs) that likely orchestrate the tumorigenesis process. Moreover, it is found that certain neurotransmitter receptors (NTRs) are subgroup-specific and can predict MB subgroup classification when combined with their associated transcription factors. Notably, targeting essential NTRs within tumors influences both the in vitro sphere-forming capability and the in vivo tumorigenic capacity of MB cells. These findings collectively provide fresh insights into comprehending the regulatory networks and cellular dynamics within MBs. Furthermore, the significance of the TF-NTR regulatory circuits is underscored as prospective biomarkers and viable therapeutic targets.

Epigenomic identification of vernalization cis-regulatory elements in winter wheat

BACKGROUND

Winter wheat undergoes vernalization, a process activated by prolonged exposure to low temperatures. During this phase, flowering signals are generated and transported to the apical meristems, stimulating the transition to the inflorescence meristem while inhibiting tiller bud elongation. Although some vernalization genes have been identified, the key cis-regulatory elements and precise mechanisms governing this process in wheat remain largely unknown.

RESULTS

In this study, we construct extensive epigenomic and transcriptomic profiling across multiple tissues-leaf, axillary bud, and shoot apex-during the vernalization of winter wheat. Epigenetic modifications play a crucial role in eliciting tissue-specific responses and sub-genome-divergent expressions during vernalization. Notably, we observe that H3K27me3 primarily regulates vernalization-induced genes and has limited influence on vernalization-repressed genes. The integration of these datasets enables the identification of 10,600 putative vernalization-related regulatory elements including distal accessible chromatin regions (ACRs) situated 30Kb upstream of VRN3, contributing to the construction of a comprehensive regulatory network. Furthermore, we discover that TaSPL7/15, integral components of the aging-related flowering pathway, interact with the VRN1 promoter and VRN3 distal regulatory elements. These interactions finely regulate their expressions, consequently impacting the vernalization process and flowering.

CONCLUSIONS

Our study offers critical insights into wheat vernalization's epigenomic dynamics and identifies the putative regulatory elements crucial for developing wheat germplasm with varied vernalization characteristics. It also establishes a vernalization-related transcriptional network, and uncovers that TaSPL7/15 from the aging pathway participates in vernalization by directly binding to the VRN1 promoter and VRN3 distal regulatory elements.

AI-Assisted Rational Design and Activity Prediction of Biological Elements for Optimizing Transcription-Factor-Based Biosensors

The rational design, activity prediction, and adaptive application of biological elements (bio-elements) are crucial research fields in synthetic biology. Currently, a major challenge in the field is efficiently designing desired bio-elements and accurately predicting their activity using vast datasets. The advancement of artificial intelligence (AI) technology has enabled machine learning and deep learning algorithms to excel in uncovering patterns in bio-element data and predicting their performance. This review explores the application of AI algorithms in the rational design of bio-elements, activity prediction, and the regulation of transcription-factor-based biosensor response performance using AI-designed elements. We discuss the advantages, adaptability, and biological challenges addressed by the AI algorithms in various applications, highlighting their powerful potential in analyzing biological data. Furthermore, we propose innovative solutions to the challenges faced by AI algorithms in the field and suggest future research directions. By consolidating current research and demonstrating the practical applications and future potential of AI in synthetic biology, this review provides valuable insights for advancing both academic research and practical applications in biotechnology.

Neph1 is required for neurite branching and is negatively regulated by the PRRXL1 homeodomain factor in the developing spinal cord dorsal horn

The cell-adhesion molecule NEPH1 is required for maintaining the structural integrity and function of the glomerulus in the kidneys. In the nervous system of Drosophila and C. elegans, it is involved in synaptogenesis and axon branching, which are essential for establishing functional circuits. In the mammalian nervous system, the expression regulation and function of Neph1 has barely been explored. In this study, we provide a spatiotemporal characterization of Neph1 expression in mouse dorsal root ganglia (DRGs) and spinal cord. After the neurogenic phase, Neph1 is broadly expressed in the DRGs and in their putative targets at the dorsal horn of the spinal cord, comprising both GABAergic and glutamatergic neurons. Interestingly, we found that PRRXL1, a homeodomain transcription factor that is required for proper establishment of the DRG-spinal cord circuit, prevents a premature expression of Neph1 in the superficial laminae of the dorsal spinal cord at E14.5, but has no regulatory effect on the DRGs or on either structure at E16.5. By chromatin immunoprecipitation analysis of the dorsal spinal cord, we identified four PRRXL1-bound regions within the Neph1 introns, suggesting that PRRXL1 directly regulates Neph1 transcription. We also showed that Neph1 is required for branching, especially at distal neurites. Together, our work showed that Prrxl1 prevents the early expression of Neph1 in the superficial dorsal horn, suggesting that Neph1 might function as a downstream effector gene for proper assembly of the DRG-spinal nociceptive circuit.

Ptip safeguards the epigenetic control of skeletal stem cell quiescence and potency in skeletogenesis

Stem cells remain in a quiescent state for long-term maintenance and preservation of potency; this process requires fine-tuning regulatory mechanisms. In this study, we identified the epigenetic landscape along the developmental trajectory of skeletal stem cells (SSCs) in skeletogenesis governed by a key regulator, Ptip (also known as Paxip1, Pax interaction with transcription-activation domain protein-1). Our results showed that Ptip is required for maintaining the quiescence and potency of SSCs, and loss of Ptip in type II collagen (Col2)+ progenitors causes abnormal activation and differentiation of SSCs, impaired growth plate morphogenesis, and long bone dysplasia. We also found that Ptip suppressed the glycolysis of SSCs through downregulation of phosphoglycerate kinase 1 (Pgk1) by repressing histone H3 lysine 27 acetylation (H3K27ac) at the promoter region. Notably, inhibition of glycolysis improved the function of SSCs despite Ptip deficiency. To the best of our knowledge, this is the first study to establish an epigenetic framework based on Ptip, which safeguards skeletal stem cell quiescence and potency through metabolic control. This framework is expected to improve SSC-based treatments of bone developmental disorders.

Characterization of activating cis-regulatory elements from the histone genes of Chlamydomonas reinhardtii

Regulation of gene expression in eukaryotes is controlled by cis-regulatory modules (CRMs). A major class of CRMs are enhancers which are composed of activating cis-regulatory elements (CREs) responsible for upregulating transcription. To date, most enhancers and activating CREs have been studied in angiosperms; in contrast, our knowledge about these key regulators of gene expression in green algae is limited. In this study, we aimed at characterizing putative activating CREs/CRMs from the histone genes of the unicellular model alga Chlamydomonas reinhardtii. To test the activity of four candidates, reporter constructs consisting of a tetramerized CRE, an established promoter, and a gene for the mCerulean3 fluorescent protein were incorporated into the nuclear genome of C. reinhardtii, and their activity was quantified by flow cytometry. Two tested candidates, Eupstr and Ehist cons, significantly upregulated gene expression and were characterized in detail. Eupstr, which originates from highly expressed genes of C. reinhardtii, is an orientation-independent CRE capable of activating both the RBCS2 and β2-tubulin promoters. Ehist cons, which is a CRM from histone genes of angiosperms, upregulates the β2-tubulin promoter in C. reinhardtii over a distance of at least 1.5 kb. The octamer motif present in Ehist cons was identified in C. reinhardtii and the related green algae Chlamydomonas incerta, Chlamydomonas schloesseri, and Edaphochlamys debaryana, demonstrating its high evolutionary conservation. The results of this investigation expand our knowledge about the regulation of gene expression in green algae. Furthermore, the characterized activating CREs/CRMs can be applied as valuable genetic tools.

Enhancer selectivity in space and time: from enhancer-promoter interactions to promoter activation

The primary regulators of metazoan gene expression are enhancers, originally functionally defined as DNA sequences that can activate transcription at promoters in an orientation-independent and distance-independent manner. Despite being crucial for gene regulation in animals, what mechanisms underlie enhancer selectivity for promoters, and more fundamentally, how enhancers interact with promoters and activate transcription, remain poorly understood. In this Review, we first discuss current models of enhancer-promoter interactions in space and time and how enhancers affect transcription activation. Next, we discuss different mechanisms that mediate enhancer selectivity, including repression, biochemical compatibility and regulation of 3D genome structure. Through 3D polymer simulations, we illustrate how the ability of 3D genome folding mechanisms to mediate enhancer selectivity strongly varies for different enhancer-promoter interaction mechanisms. Finally, we discuss how recent technical advances may provide new insights into mechanisms of enhancer-promoter interactions and how technical biases in methods such as Hi-C and Micro-C and imaging techniques may affect their interpretation.

The evolutionary consequences of interactions between the epigenome, the genome and the environment

The epigenome is the suite of interacting chemical marks and molecules that helps to shape patterns of development, phenotypic plasticity and gene regulation, in part due to its responsiveness to environmental stimuli. There is increasing interest in understanding the functional and evolutionary importance of this sensitivity under ecologically realistic conditions. Observations that epigenetic variation abounds in natural populations have prompted speculation that it may facilitate evolutionary responses to rapid environmental perturbations, such as those occurring under climate change. A frequent point of contention is whether epigenetic variants reflect genetic variation or are independent of it. The genome and epigenome often appear tightly linked and interdependent. While many epigenetic changes are genetically determined, the converse is also true, with DNA sequence changes influenced by the presence of epigenetic marks. Understanding how the epigenome, genome and environment interact with one another is therefore an essential step in explaining the broader evolutionary consequences of epigenomic variation. Drawing on results from experimental and comparative studies carried out in diverse plant and animal species, we synthesize our current understanding of how these factors interact to shape phenotypic variation in natural populations, with a focus on identifying similarities and differences between taxonomic groups. We describe the main components of the epigenome and how they vary within and between taxa. We review how variation in the epigenome interacts with genetic features and environmental determinants, with a focus on the role of transposable elements (TEs) in integrating the epigenome, genome and environment. And we look at recent studies investigating the functional and evolutionary consequences of these interactions. Although epigenetic differentiation in nature is likely often a result of drift or selection on stochastic epimutations, there is growing evidence that a significant fraction of it can be stably inherited and could therefore contribute to evolution independently of genetic change.

Single-cell nascent RNA sequencing unveils coordinated global transcription

Transcription is the primary regulatory step in gene expression. Divergent transcription initiation from promoters and enhancers produces stable RNAs from genes and unstable RNAs from enhancers1,2. Nascent RNA capture and sequencing assays simultaneously measure gene and enhancer activity in cell populations3. However, fundamental questions about the temporal regulation of transcription and enhancer-gene coordination remain unanswered, primarily because of the absence of a single-cell perspective on active transcription. In this study, we present scGRO-seq-a new single-cell nascent RNA sequencing assay that uses click chemistry-and unveil coordinated transcription throughout the genome. We demonstrate the episodic nature of transcription and the co-transcription of functionally related genes. scGRO-seq can estimate burst size and frequency by directly quantifying transcribing RNA polymerases in individual cells and can leverage replication-dependent non-polyadenylated histone gene transcription to elucidate cell cycle dynamics. The single-nucleotide spatial and temporal resolution of scGRO-seq enables the identification of networks of enhancers and genes. Our results suggest that the bursting of transcription at super-enhancers precedes bursting from associated genes. By imparting insights into the dynamic nature of global transcription and the origin and propagation of transcription signals, we demonstrate the ability of scGRO-seq to investigate the mechanisms of transcription regulation and the role of enhancers in gene expression.

Regulation of 3D genome organization during T cell activation

Within the three-dimensional (3D) nuclear space, the genome organizes into a series of orderly structures that impose important influences on gene regulation. T lymphocytes, crucial players in adaptive immune responses, undergo intricate transcriptional remodeling upon activation, leading to differentiation into specific effector and memory T cell subsets. Recent evidence suggests that T cell activation is accompanied by dynamic changes in genome architecture at multiple levels, providing a unique biological context to explore the functional relevance and molecular mechanisms of 3D genome organization. Here, we summarize recent advances that link the reorganization of genome architecture to the remodeling of transcriptional programs and conversion of cell fates during T cell activation and differentiation. We further discuss how various chromatin architecture regulators, including CCCTC-binding factor and several transcription factors, collectively modulate the genome architecture during this process.

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

Background

STARR-seq and other massively-parallel reporter assays are widely used to discover functional enhancers in transfected cell models, which can be confounded by plasmid vector-induced type-I interferon immune responses and lack the multicellular environment and endogenous chromatin state of complex mammalian tissues.

Results

Here, we describe HDI-STARR-seq, which combines STARR-seq plasmid library delivery to the liver, by hydrodynamic tail vein injection (HDI), with reporter RNA transcriptional initiation driven by a minimal Albumin promoter, which we show is essential for mouse liver STARR-seq enhancer activity assayed 7 days after HDI. Importantly, little or no vector-induced innate type-I interferon responses were observed. Comparisons of HDI-STARR-seq activity between male and female mouse livers and in livers from males treated with an activating ligand of the transcription factor CAR (Nr1i3) identified many condition-dependent enhancers linked to condition-specific gene expression. Further, thousands of active liver enhancers were identified using a high complexity STARR-seq library comprised of ~ 50,000 genomic regions released by DNase-I digestion of mouse liver nuclei. When compared to stringently inactive library sequences, the active enhancer sequences identified were highly enriched for liver open chromatin regions with activating histone marks (H3K27ac, H3K4me1, H3K4me3), were significantly closer to gene transcriptional start sites, and were significantly depleted of repressive (H3K27me3, H3K9me3) and transcribed region histone marks (H3K36me3).

Conclusions

HDI-STARR-seq offers substantial improvements over current methodologies for large scale, functional profiling of enhancers, including condition-dependent enhancers, in liver tissue in vivo, and can be adapted to characterize enhancer activities in a variety of species and tissues by selecting suitable tissue- and species-specific promoter sequences.

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

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

Genome-wide enhancer RNA profiling adds molecular links between genetic variation and human cancers

BACKGROUND

Dysregulation of enhancer transcription occurs in multiple cancers. Enhancer RNAs (eRNAs) are transcribed products from enhancers that play critical roles in transcriptional control. Characterizing the genetic basis of eRNA expression may elucidate the molecular mechanisms underlying cancers.

METHODS

Initially, a comprehensive analysis of eRNA quantitative trait loci (eRNAQTLs) was performed in The Cancer Genome Atlas (TCGA), and functional features were characterized using multi-omics data. To establish the first eRNAQTL profiles for colorectal cancer (CRC) in China, epigenomic data were used to define active enhancers, which were subsequently integrated with transcription and genotyping data from 154 paired CRC samples. Finally, large-scale case-control studies (34,585 cases and 69,544 controls) were conducted along with multipronged experiments to investigate the potential mechanisms by which candidate eRNAQTLs affect CRC risk.

RESULTS

A total of 300,112 eRNAQTLs were identified across 30 different cancer types, which exert their influence on eRNA transcription by modulating chromatin status, binding affinity to transcription factors and RNA-binding proteins. These eRNAQTLs were found to be significantly enriched in cancer risk loci, explaining a substantial proportion of cancer heritability. Additionally, tumor-specific eRNAQTLs exhibited high responsiveness to the development of cancer. Moreover, the target genes of these eRNAs were associated with dysregulated signaling pathways and immune cell infiltration in cancer, highlighting their potential as therapeutic targets. Furthermore, multiple ethnic population studies have confirmed that an eRNAQTL rs3094296-T variant decreases the risk of CRC in populations from China (OR = 0.91, 95%CI 0.88-0.95, P = 2.92 × 10-7) and Europe (OR = 0.92, 95%CI 0.88-0.95, P = 4.61 × 10-6). Mechanistically, rs3094296 had an allele-specific effect on the transcription of the eRNA ENSR00000155786, which functioned as a transcriptional activator promoting the expression of its target gene SENP7. These two genes synergistically suppressed tumor cell proliferation. Our curated list of variants, genes, and drugs has been made available in CancereRNAQTL ( http://canernaqtl.whu.edu.cn/#/ ) to serve as an informative resource for advancing this field.

CONCLUSION

Our findings underscore the significance of eRNAQTLs in transcriptional regulation and disease heritability, pinpointing the potential of eRNA-based therapeutic strategies in cancers.

Mechanistic analysis of enhancer sequences in the estrogen receptor transcriptional program

Estrogen Receptor α (ERα) is a major lineage determining transcription factor (TF) in mammary gland development. Dysregulation of ERα-mediated transcriptional program results in cancer. Transcriptomic and epigenomic profiling of breast cancer cell lines has revealed large numbers of enhancers involved in this regulatory program, but how these enhancers encode function in their sequence remains poorly understood. A subset of ERα-bound enhancers are transcribed into short bidirectional RNA (enhancer RNA or eRNA), and this property is believed to be a reliable marker of active enhancers. We therefore analyze thousands of ERα-bound enhancers and build quantitative, mechanism-aware models to discriminate eRNAs from non-transcribing enhancers based on their sequence. Our thermodynamics-based models provide insights into the roles of specific TFs in ERα-mediated transcriptional program, many of which are supported by the literature. We use in silico perturbations to predict TF-enhancer regulatory relationships and integrate these findings with experimentally determined enhancer-promoter interactions to construct a gene regulatory network. We also demonstrate that the model can prioritize breast cancer-related sequence variants while providing mechanistic explanations for their function. Finally, we experimentally validate the model-proposed mechanisms underlying three such variants.

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.

Mapping putative enhancers in mouse oocytes and early embryos reveals TCF3/12 as key folliculogenesis regulators

Dynamic epigenomic reprogramming occurs during mammalian oocyte maturation and early development. However, the underlying transcription circuitry remains poorly characterized. By mapping cis-regulatory elements using H3K27ac, we identified putative enhancers in mouse oocytes and early embryos distinct from those in adult tissues, enabling global transitions of regulatory landscapes around fertilization and implantation. Gene deserts harbour prevalent putative enhancers in fully grown oocytes linked to oocyte-specific genes and repeat activation. Embryo-specific enhancers are primed before zygotic genome activation and are restricted by oocyte-inherited H3K27me3. Putative enhancers in oocytes often manifest H3K4me3, bidirectional transcription, Pol II binding and can drive transcription in STARR-seq and a reporter assay. Finally, motif analysis of these elements identified crucial regulators of oogenesis, TCF3 and TCF12, the deficiency of which impairs activation of key oocyte genes and folliculogenesis. These data reveal distinctive regulatory landscapes and their interacting transcription factors that underpin the development of mammalian oocytes and early embryos.

Molecular insights into regulatory RNAs in the cellular machinery

It is apparent that various functional units within the cellular machinery are derived from RNAs. The evolution of sequencing techniques has resulted in significant insights into approaches for transcriptome studies. Organisms utilize RNA to govern cellular systems, and a heterogeneous class of RNAs is involved in regulatory functions. In particular, regulatory RNAs are increasingly recognized to participate in intricately functioning machinery across almost all levels of biological systems. These systems include those mediating chromatin arrangement, transcription, suborganelle stabilization, and posttranscriptional modifications. Any class of RNA exhibiting regulatory activity can be termed a class of regulatory RNA and is typically represented by noncoding RNAs, which constitute a substantial portion of the genome. These RNAs function based on the principle of structural changes through cis and/or trans regulation to facilitate mutual RNA‒RNA, RNA‒DNA, and RNA‒protein interactions. It has not been clearly elucidated whether regulatory RNAs identified through deep sequencing actually function in the anticipated mechanisms. This review addresses the dominant properties of regulatory RNAs at various layers of the cellular machinery and covers regulatory activities, structural dynamics, modifications, associated molecules, and further challenges related to therapeutics and deep learning.

The roles of nuclear orphan receptor NR2F6 in anti-viral innate immunity

Proper transcription regulation by key transcription factors, such as IRF3, is critical for anti-viral defense. Dynamics of enhancer activity play important roles in many biological processes, and epigenomic analysis is used to determine the involved enhancers and transcription factors. To determine new transcription factors in anti-DNA-virus response, we have performed H3K27ac ChIP-Seq and identified three transcription factors, NR2F6, MEF2D and MAFF, in promoting HSV-1 replication. NR2F6 promotes HSV-1 replication and gene expression in vitro and in vivo, but not dependent on cGAS/STING pathway. NR2F6 binds to the promoter of MAP3K5 and activates AP-1/c-Jun pathway, which is critical for DNA virus replication. On the other hand, NR2F6 is transcriptionally repressed by c-Jun and forms a negative feedback loop. Meanwhile, cGAS/STING innate immunity signaling represses NR2F6 through STAT3. Taken together, we have identified new transcription factors and revealed the underlying mechanisms involved in the network between DNA viruses and host cells.

Construction of cytomegalovirus promoter-driven gene expression system in Laodelphax striatellus

The small brown planthopper (SBPH, Laodelphax striatellus) is a significant rice pest, responsible for transmitting rice stripe virus (RSV) in a persistent and propagative manner. RSV is one of the most detrimental rice viruses, causing rice stripe disease, which results in considerable loss of rice grain yield. While RNA interference and gene knockout techniques have enabled gene downregulation in SBPH, no system currently exists for the overexpression of endogenous or exogenous genes. Consequently, the development of a protein expression system for SBPH is imperative to serve as a technical foundation for pest control and gene function investigations. This study aimed to construct an expression vector using the promoter of the constitutive-expressed tubulin gene of SBPH, and promoter of human cytomegalovirus (CMV). Fluorescence experiments demonstrated that both tubulin and CMV promoter could drive green fluorescent protein (GFP) expression in SBPH, and could also facilitate the expression of a nucleocapsid protein (NP) -GFP fusion protein containing viral NP with comparable efficiency. Through expression vector optimization, we have identified that the 3 tandem CMV promoters display a significantly higher promoter activity compared with both the 2 tandem CMV promoters and the single CMV promoter. In addition, the incorporation of Star polycation nanoparticles significantly enhanced the expression efficiency in SBPH. These results provide a promising technical platform for investigating gene functions in SBPH.

Conserved Cis-Acting Range Extender Element Mediates Extreme Long-Range Enhancer Activity in Mammals

While most mammalian enhancers regulate their cognate promoters over moderate distances of tens of kilobases (kb), some enhancers act over distances in the megabase range. The sequence features enabling such extreme-distance enhancer-promoter interactions remain elusive. Here, we used in vivo enhancer replacement experiments in mice to show that short- and medium-range enhancers cannot initiate gene expression at extreme-distance range. We uncover a novel conserved cis-acting element, Range EXtender (REX), that confers extreme-distance regulatory activity and is located next to a long-range enhancer of Sall1. The REX element itself has no endogenous enhancer activity. However, addition of the REX to other short- and mid-range enhancers substantially increases their genomic interaction range. In the most extreme example observed, addition of the REX increased the range of an enhancer by an order of magnitude, from its native 71kb to 840kb. The REX element contains highly conserved [C/T]AATTA homeodomain motifs. These motifs are enriched around long-range limb enhancers genome-wide, including the ZRS, a benchmark long-range limb enhancer of Shh. Mutating the [C/T]AATTA motifs within the ZRS does not affect its limb-specific enhancer activity at short range, but selectively abolishes its long-range activity, resulting in severe limb reduction in knock-in mice. In summary, we identify a sequence signature globally associated with long-range enhancer-promoter interactions and describe a prototypical REX element that is necessary and sufficient to confer extreme-distance gene activation by remote enhancers.

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.

H3K4me1 facilitates promoter-enhancer interactions and gene activation during embryonic stem cell differentiation

Histone H3 lysine 4 mono-methylation (H3K4me1) marks poised or active enhancers. KMT2C (MLL3) and KMT2D (MLL4) catalyze H3K4me1, but their histone methyltransferase activities are largely dispensable for transcription during early embryogenesis in mammals. To better understand the role of H3K4me1 in enhancer function, we analyze dynamic enhancer-promoter (E-P) interactions and gene expression during neural differentiation of the mouse embryonic stem cells. We found that KMT2C/D catalytic activities were only required for H3K4me1 and E-P contacts at a subset of candidate enhancers, induced upon neural differentiation. By contrast, a majority of enhancers retained H3K4me1 in KMT2C/D catalytic mutant cells. Surprisingly, H3K4me1 signals at these KMT2C/D-independent sites were reduced after acute depletion of KMT2B, resulting in aggravated transcriptional defects. Our observations therefore implicate KMT2B in the catalysis of H3K4me1 at enhancers and provide additional support for an active role of H3K4me1 in enhancer-promoter interactions and transcription in mammalian cells.

Paramutation at the maize pl1 locus is associated with RdDM activity at distal tandem repeats

Exceptions to Mendelian inheritance often highlight novel chromosomal behaviors. The maize Pl1-Rhoades allele conferring plant pigmentation can display inheritance patterns deviating from Mendelian expectations in a behavior known as paramutation. However, the chromosome features mediating such exceptions remain unknown. Here we show that small RNA production reflecting RNA polymerase IV function within a distal downstream set of five tandem repeats is coincident with meiotically-heritable repression of the Pl1-Rhoades transcription unit. A related pl1 haplotype with three, but not one with two, repeat units also displays the trans-homolog silencing typifying paramutations. 4C interactions, CHD3a-dependent small RNA profiles, nuclease sensitivity, and polyadenylated RNA levels highlight a repeat subregion having regulatory potential. Our comparative and mutant analyses show that transcriptional repression of Pl1-Rhoades correlates with 24-nucleotide RNA production and cytosine methylation at this subregion indicating the action of a specific DNA-dependent RNA polymerase complex. These findings support a working model in which pl1 paramutation depends on trans-chromosomal RNA-directed DNA methylation operating at a discrete cis-linked and copy-number-dependent transcriptional regulatory element.