Transcribed enhancers lead waves of coordinated transcription in transitioning mammalian cells

Long Non-Coding RNAs and RNA-Binding Proteins in Pancreatic Cancer Development and Progression

Pancreatic ductal adenocarcinoma (PDAC) is the most common type of pancreatic cancer and is responsible for about 467,000 cancer deaths annually. An oftentimes asymptomatic early phase of this disease results in a delayed diagnosis, and patients often present with advanced disease. Current treatment options have limited survival benefits, and only a minor patient population carries actionable genomic alterations. Hence, innovative personalized treatment strategies that consider molecular, cellular and functional analyses are urgently needed for pancreatic cancer patients. However, the majority of the genetic alterations found in PDAC are currently undruggable, or patients' response is not as expected. Therefore, non-genomic biomarkers and alternative molecular targets should be considered in order to advance the clinical management of PDAC patients. In line with this, recent gene expression and single-cell transcriptome analyses have identified molecular subtypes and transcriptional cell states that affect disease progression and drug efficiency. In this review, we will introduce long non-coding RNAs (lncRNAs) as well as RNA-binding proteins (RBPs) that are able to modulate the transcriptome of a cell through diverse mechanisms, thereby contributing to disease progression. We will provide a brief overview about the general functions of lncRNAs and RBPs, respectively. Subsequently, we will highlight selected lncRNAs and RBPs that have been shown to play a role in PDAC development, progression and drug response. Finally, we will present strategies aiming to interfere with the expression and function of lncRNAs and RBPs.

Drug-induced cis-regulatory elements in human hepatocytes affect molecular phenotypes associated with adverse reactions

Genomic variation drives phenotypic diversity, including individual differences in drug response. While coding polymorphisms linked to drug efficacy and adverse reactions are well characterized, the contribution of noncoding regulatory elements remains underexplored. Using CAGE (Cap Analysis of Gene Expression), profiling transcription initiations of mRNAs and enhancer RNAs, we identify candidate cis-regulatory elements (CREs) and assessed their activities simultaneously in HepG2 cells expressing the drug-responsive transcription factor pregnane X receptor (PXR). Comparison with GWAS data reveals strong enrichment of the drug-induced CREs near variants associated with bilirubin and vitamin D levels. Among those bound by PXR in primary hepatocytes, we identify enhancers of UGT1A1, TSKU, and CYP24A1 and functional alleles that alter regulatory activities. We also find that TSKU influences expression of vitamin D-metabolizing enzymes. This study expands the landscape of PXR-mediated regulatory elements and uncovers noncoding variants impacting drug response, providing insights into the genomic basis of adverse drug reactions.

Recent advances in designing synthetic plant regulatory modules

Introducing novel functions in plants through synthetic multigene circuits requires strict transcriptional regulation. Currently, the use of natural regulatory modules in synthetic circuits is hindered by our limited knowledge of complex plant regulatory mechanisms, the paucity of characterized promoters, and the possibility of crosstalk with endogenous circuits. Synthetic regulatory modules can overcome these limitations. This article introduces an integrative de novo approach for designing plant synthetic promoters by utilizing the available online tools and databases. The recent achievements in designing and validating synthetic plant promoters, enhancers, transcription factors, and the challenges of establishing synthetic circuits in plants are also discussed.

Validation of selective catalytic BmCBP inhibitors that regulate the Bm30K-24 protein expression in silkworm, Bombyx mori

The cAMP response element binding protein (CREB)-binding protein (CBP) is a histone acetyltransferase that plays an indispensable role in regulating the acetylation of histone and non-histone proteins. Recently, it has been discovered that chemical inhibitors A485 and C646 can bind to Bombyx mori's CBP (BmCBP) and inhibit its acetyltransferase activity. Notably, the binding ability of A485 with BmCBP showed a very low Kd value of 48 nM by surface plasmon resonance (SPR) test. Further identification showed that both A485 and C646 can decrease the acetylation level of known substrate H3K27 and only 1 μM of A485 can almost completely inhibit the acetylation of H3K27, suggesting that A485 is an effective inhibitor of BmCBP's acetyltransferase activity. Moreover, it was confirmed that A485 could downregulate the expression of acetylated Bm30K-24 protein at a post-translational level through acetylation modification by BmCBP. Additionally, it was found that A485 can downregulate the stability of Bm30K-24 and improve its ubiquitination level, suggesting that the acetylation modification by BmCBP could compete with ubiquitination modification at the same lysine site on Bm30K-24, thereby affecting its protein stability. Here, we predict that A485 may be a potent CBP acetyltransferase inhibitor which could be utilized to inhibit acetyltransferase activity in insects, including silkworms.

A transcription coupling model for how enhancers communicate with their target genes

How enhancers communicate with their target genes to influence transcription is an unresolved question of fundamental importance. Current models of the mechanism of enhancer-target gene or enhancer-promoter (E-P) communication are transcription-factor-centric and underappreciate major findings, including that enhancers are themselves transcribed by RNA polymerase II, which correlates with enhancer activity. In this Perspective, we posit that enhancer transcription and its products, enhancer RNAs, are elementary components of enhancer-gene communication. Specifically, we discuss the possibility that transcription at enhancers and at their cognate genes are linked and that this coupling is at the basis of how enhancers communicate with their targets. This model of transcriptional coupling between enhancers and their target genes is supported by growing experimental evidence and represents a synthesis of recent key discoveries.

True length of diverse capped RNA sequencing (TLDR-seq): 5'-3'-end sequencing of capped RNAs regardless of 3'-end status

Analysis of transcript function is greatly aided by knowledge of the full-length RNA sequence. New long-read sequencing enabled by Oxford Nanopore and PacBio devices have the potential to provide full-length transcript information; however, standard methods still lack the ability to capture true RNA 5' ends and select for polyadenylated (pA+) transcripts only. Here, we present a method that, by utilizing cap trapping and 3'-end adapter ligation, sequences transcripts between their exact 5' and 3' ends regardless of polyadenylation status and without the need for ribosomal RNA depletion, with the ability to characterize polyadenylation length of RNAs, if any. The method shows high reproducibility, can faithfully detect 5' ends, 3' ends and splice junctions, and produces gene-expression estimates that are highly correlated to those of short-read sequencing techniques. We also demonstrate that the method can detect and sequence full-length nonadenylated (pA-) RNAs, including long noncoding RNAs, promoter upstream transcripts, and enhancer RNAs, and present cases where pA+ and pA- RNAs show preferences for different but closely located transcription start sites. Our method is therefore useful for the characterization of diverse capped RNA species and analysis of relationships between transcription initiation, termination, and RNA processing.

Transcriptional enhancers in human neuronal differentiation provide clues to neuronal disorders

Genome-wide association studies (GWASs) have identified thousands of variants associated with complex phenotypes, including neuropsychiatric disorders. To better understand their pathogenesis, it is necessary to identify the functional roles of these variants, which are largely located in non-coding DNA regions. Here, we employ a human mesencephalic neuronal cell differentiation model, LUHMES, with sensitive and high-resolution methods to discover enhancers (NET-CAGE), perform DNA conformation analysis (Capture Hi-C) to link enhancers to their target genes, and finally validate selected interactions. We expand the number of known enhancers active in differentiating human LUHMES neurons to 47,350, and find overlap with GWAS variants for Parkinson's disease and schizophrenia. Our findings reveal a fine-tuned regulation of human neuronal differentiation, even between adjacent developmental stages; provide a valuable resource for further studies on neuronal development, regulation, and disorders; and emphasize the importance of exploring the vast regulatory potential of non-coding DNA and enhancers.

Hallmarks of Brain Plasticity

Cerebral plasticity is the ability of the brain to change and adapt in response to experience or learning. Its hallmarks are developmental flexibility, complex interactions between genetic and environmental influences, and structural-functional changes comprising neurogenesis, axonal sprouting, and synaptic remodeling. Studies on brain plasticity have important practical implications. The molecular characteristics of changes in brain plasticity may reveal disease course and the rehabilitative potential of the patient. Neurological disorders are linked with numerous cerebral non-coding RNAs (ncRNAs), in particular, microRNAs; the discovery of their essential role in gene regulation was recently recognized and awarded a Nobel Prize in Physiology or Medicine in 2024. Herein, we review the association of brain plasticity and its homeostasis with ncRNAs, which make them putative targets for RNA-based diagnostics and therapeutics. New insight into the concept of brain plasticity may provide additional perspectives on functional recovery following brain damage. Knowledge of this phenomenon will enable physicians to exploit the potential of cerebral plasticity and regulate eloquent networks with timely interventions. Future studies may reveal pathophysiological mechanisms of brain plasticity at macro- and microscopic levels to advance rehabilitation strategies and improve quality of life in patients with neurological diseases.

Application Strategies of Super-Enhancer RNA in Cardiovascular Diseases

Cardiovascular diseases (CVDs) are a leading cause of death worldwide, and new therapeutic strategies are urgently needed. In recent years, enhancer RNAs (eRNAs) have gradually attracted attention because they offer new directions for the treatment of CVDs. Super-enhancer RNAs (seRNAs) are a subset of non-coding RNAs that are transcribed from regions of the genome known as super enhancers, which are large clusters of enhancers with a high density of transcription factors and cofactors. These regions play a pivotal role in regulating genes involved in cell identity and disease progression. This article reviews the characteristics of seRNAs, their expression patterns, and regulatory mechanisms in the cardiovascular system. We also explore their role in the occurrence and development of CVDs, as well as their potential as diagnostic biomarkers and therapeutic targets. Currently, therapies targeting seRNAs are a research hotspot. The development of specific inhibitors or activators is expected to facilitate precise interventions for CVDs. In addition, the use of gene editing techniques to modify relevant eRNA introduces new possibilities for disease treatment. This review aims to provide a comprehensive overview of seRNAs in CVDs and discusses their potential as a novel class of therapeutic targets.

Update of the FANTOM web resource: enhancement for studying noncoding genomes

The FANTOM web resource (https://fantom.gsc.riken.jp/) has been a unique resource for studying mammalian genomes, which is built on the research activities conducted in the international collaborative project FANTOM (Functional ANnoTation Of the Mammalian genome). In recent updates, we expanded annotations for long non-coding RNAs (lncRNAs) and transcribed cis-regulatory elements (CREs). The former was derived from the large-scale lncRNA perturbations in induced pluripotent stem cells (iPSCs) and integrative analysis of Hi-C data conducted in the sixth iteration of the project (FANTOM6). The resulting annotations of lncRNAs, according to the impact on cellular and molecular phenotypes and the potential RNA-chromatin interactions, are accessible via the interactive ZENBU-Reports framework. The latter involves a new platform, fanta.bio (https://fanta.bio/), which collects transcribed CREs identified via use of an extended dataset of CAGE profiles. The CREs, with their annotations including genetic and epigenetic information, are accessible via a dedicated interface as well as the UCSC Genome Browser Database. These updates offer enhanced opportunities to investigate the functions of non-coding regions within mammalian genomes.

Super-enhancer-driven ameboidal-type cell migration-related MMP14 expression in tongue squamous cell carcinoma switched by BATF and ATF3

BACKGROUND

Tongue squamous cell carcinoma (TSCC) exhibits an aggressive biological behavior of lymph node and distant metastasis, which contributes to poorer prognosis and results in tongue function loss or death. In addition to known regulators and pathways of cell migration in TSCC, it is important to uncover pivotal switches governing tumor metastasis.

METHODS

Cancer cell migration-associated transcriptional and epigenetic characteristics were profiled in TSCC, and the specific super-enhancers (SEs) were identified. Molecular function and mechanism studies were used to investigate the pivotal switches in TSCC metastasis.

RESULTS

Ameboidal-type cell migration-related genes accompanied by transcriptional and epigenetic activity were enriched in TSCC. Meanwhile, the higher-ranked SE-related genes showed significant differences between 43 paired tumor and normal samples from the TCGA TSCC cohort. In addition, key motifs were detected in SE regions, and transcription factor-related expression levels were significantly associated with TSCC survival status. Notably, BATF and ATF3 regulated the expression of ameboidal-type cell migration-related MMP14 by switching the interaction with the SE region.

CONCLUSION

SEs and related key motifs transcriptional regulate tumor metastasis-associated MMP14 and might be potential therapeutic targets for TSCC.

Long noncoding RNA VENTHEART is required for ventricular cardiomyocyte specification and function

RATIONALE

Cardiac-expressed long noncoding RNAs (lncRNAs) are important for cardiomyocyte (CM) differentiation and function. Several lncRNAs have been identified and characterized for early CM lineage commitment, however those in later CM lineage specification and maturation remain less well studied. Moreover, unique atrial / ventricular lncRNA expression has never been studied in detail.

OBJECTIVES

Here, we characterized a novel ventricular myocyte-restricted lncRNA, not expressed in atrial myocytes, and conserved only in primates.

METHODS AND RESULTS

First, we performed single cell RNA-seq on human pluripotent stem cell derived cardiomyocytes (hPSC-CM) at the late stages of 2, 6 and 12 weeks of differentiation. Weighted correlation network analysis identified core gene modules, including a set of lncRNAs highly abundant and predominantly expressed in the human heart. A lncRNA (we call VENTHEART, VHRT) co-expressed with cardiac maturation and ventricular-specific genes MYL2 and MYH7, and was expressed in fetal and adult human ventricles, but not atria. CRISPR-mediated deletion of the VHRT gene led to impaired CM sarcomere formation and significant disruption of the ventricular CM gene program. Indeed, a similar disruption was not observed in VHRT KO hPSC-derived atrial CM, suggesting that VHRT exhibits only ventricular myocyte subtype-specific effects. Optical recordings validated that loss of VHRT significantly prolonged action potential duration at 90 % repolarization (APD90) for ventricular-like, but not atrial-like, CMs.

CONCLUSION

This reports the first lncRNA that is exclusively required for proper ventricular, and not atrial, CM specification and function.

Transcriptional machinery as an architect of genome structure

Chromatin organization, facilitated by compartmentalization and loop extrusion, is crucial for proper gene expression and cell viability. Transcription has long been considered important for shaping genome architecture due to its pervasive activity across the genome and impact on the local chromatin environment. Although earlier studies suggested a minimal contribution of transcription to shaping global genome structure, recent insights from high-resolution chromatin contact mapping, polymer simulations, and acute perturbations have revealed its critical role in dynamic chromatin organization at the level of active genes and enhancer-promoter interactions. In this review, we discuss these latest advances, highlighting the direct interplay between transcriptional machinery and loop extrusion. Finally, we explore how transcription of genes and non-coding regulatory elements may contribute to the specificity of gene regulation, focusing on enhancers as sites of targeted cohesin loading.

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. 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 and a highly epistatic sequence syntax of cooperative and competitive motif interactions. ProCapNet models of steady-state RAMPAGE profiles 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.

SCAR-6 elncRNA locus epigenetically regulates PROZ and modulates coagulation and vascular function

In this study, we characterize a novel lncRNA-producing gene locus that we name Syntenic Cardiovascular Conserved Region-Associated lncRNA-6 (scar-6) and functionally validate its role in coagulation and cardiovascular function. A 12-bp deletion of the scar-6 locus in zebrafish (scar-6gib007Δ12/Δ12) results in cranial hemorrhage and vascular permeability. Overexpression, knockdown and rescue with the scar-6 lncRNA modulates hemostasis in zebrafish. Molecular investigation reveals that the scar-6 lncRNA acts as an enhancer lncRNA (elncRNA), and controls the expression of prozb, an inhibitor of factor Xa, through an enhancer element in the scar-6 locus. The scar-6 locus suppresses loop formation between prozb and scar-6 sequences, which might be facilitated by the methylation of CpG islands via the prdm14-PRC2 complex whose binding to the locus might be stabilized by the scar-6 elncRNA transcript. Binding of prdm14 to the scar-6 locus is impaired in scar-6gib007Δ12/Δ12 zebrafish. Finally, activation of the PAR2 receptor in scar-6gib007Δ12/Δ12 zebrafish triggers NF-κB-mediated endothelial cell activation, leading to vascular dysfunction and hemorrhage. We present evidence that the scar-6 locus plays a role in regulating the expression of the coagulation cascade gene prozb and maintains vascular homeostasis.

Transcriptional coactivator MED15 is required for beta cell maturation

Mediator, a co-regulator complex required for RNA Polymerase II activity, interacts with tissue-specific transcription factors to regulate development and maintain homeostasis. We observe reduced Mediator subunit MED15 expression in endocrine hormone-producing pancreatic islets isolated from people living with type 2 diabetes and sought to understand how MED15 and Mediator control gene expression programs important for the function of insulin-producing β-cells. Here we show that Med15 is expressed during mouse β-cell development and maturation. Knockout of Med15 in mouse β-cells causes defects in β-cell maturation without affecting β-cell mass or insulin expression. ChIP-seq and co-immunoprecipitation analyses found that Med15 binds β-cell transcription factors Nkx6-1 and NeuroD1 to regulate key β-cell maturation genes. In support of a conserved role during human development, human embryonic stem cell-derived β-like cells, genetically engineered to express high levels of MED15, express increased levels of maturation markers. We provide evidence of a conserved role for Mediator in β-cell maturation and demonstrate an additional layer of control that tunes β-cell transcription factor function.

RNA-driven phase transitions in biomolecular condensates

RNAs and RNA-binding proteins can undergo spontaneous or active condensation into phase-separated liquid-like droplets. These condensates are cellular hubs for various physiological processes, and their dysregulation leads to diseases. Although RNAs are core components of many cellular condensates, the underlying molecular determinants for the formation, regulation, and function of ribonucleoprotein condensates have largely been studied from a protein-centric perspective. Here, we highlight recent developments in ribonucleoprotein condensate biology with a particular emphasis on RNA-driven phase transitions. We also present emerging future directions that might shed light on the role of RNA condensates in spatiotemporal regulation of cellular processes and inspire bioengineering of RNA-based therapeutics.

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.

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.

Acetalax (Oxyphenisatin Acetate, NSC 59687) and Bisacodyl Cause Oncosis in Triple-Negative Breast Cancer Cell Lines by Poisoning the Ion Exchange Membrane Protein TRPM4

Triple-negative breast cancer (TNBC) is clinically aggressive and relatively unresponsive to current therapies. Therefore, the development of new anticancer agents is needed to satisfy clinical needs. Oxyphenisatin acetate (Acetalax), which had been used as a laxative, has recently been reported to have anticancer activity in murine models. In this study, we demonstrate that Acetalax and its diphenolic laxative structural analogue bisacodyl (Dulcolax) exhibit potent antiproliferative activity in TNBC cell lines and cause oncosis, a nonapoptotic cell death characterized by cellular and nuclear swelling and cell membrane blebbing, leading to mitochondrial dysfunction, ATP depletion, and enhanced immune and inflammatory responses. Mechanistically, we provide evidence that transient receptor potential melastatin member 4 (TRPM4) is poisoned by Acetalax and bisacodyl in MDA-MB468, BT549, and HS578T TNBC cells. MDA-MB231 and MDA-MB436 TNBC cells without endogenous TRPM4 expression as well as TRPM4-knockout TNBC cells were found to be Acetalax- and bisacodyl-resistant. Conversely, ectopic expression of TRPM4 sensitized MDA-MB231 and MDA-MB436 cells to Acetalax. TRPM4 was also lost in cells with acquired Acetalax resistance. Moreover, TRPM4 is rapidly degraded by the ubiquitin-proteasome system upon acute exposure to Acetalax and bisacodyl. Together, these results demonstrate that TRPM4 is a previously unknown target of Acetalax and bisacodyl and that TRPM4 expression in cancer cells is a predictor of Acetalax and bisacodyl efficacy and could be used for the clinical development of these drugs as anticancer agents.

SIGNIFICANCE

Acetalax and bisacodyl kill cancer cells by causing oncosis following poisoning of the plasma membrane sodium transporter TRPM4 and represent a new therapeutic approach for TNBC.

Pan-cancer transcriptome analysis reveals widespread regulation through alternative tandem transcription initiation

Abnormal transcription initiation from alternative first exon has been reported to promote tumorigenesis. However, the prevalence and impact of gene expression regulation mediated by alternative tandem transcription initiation were mostly unknown in cancer. Here, we developed a robust computational method to analyze alternative tandem transcription start site (TSS) usage from standard RNA sequencing data. Applying this method to pan-cancer RNA sequencing datasets, we observed widespread dysregulation of tandem TSS usage in tumors, many of which were independent of changes in overall expression level or alternative first exon usage. We showed that the dynamics of tandem TSS usage was associated with epigenomic modulation. We found that significant 5' untranslated region shortening of gene TIMM13 contributed to increased protein production, and up-regulation of TIMM13 by CRISPR-mediated transcriptional activation promoted proliferation and migration of lung cancer cells. Our findings suggest that dysregulated tandem TSS usage represents an addtional layer of cancer-associated transcriptome alterations.

An atlas of transcribed enhancers across helper T cell diversity for decoding human diseases

Transcribed enhancer maps can reveal nuclear interactions underpinning each cell type and connect specific cell types to diseases. Using a 5' single-cell RNA sequencing approach, we defined transcription start sites of enhancer RNAs and other classes of coding and noncoding RNAs in human CD4+ T cells, revealing cellular heterogeneity and differentiation trajectories. Integration of these datasets with single-cell chromatin profiles showed that active enhancers with bidirectional RNA transcription are highly cell type-specific and that disease heritability is strongly enriched in these enhancers. The resulting cell type-resolved multimodal atlas of bidirectionally transcribed enhancers, which we linked with promoters using fine-scale chromatin contact maps, enabled us to systematically interpret genetic variants associated with a range of immune-mediated diseases.

Uncovering the dynamics and consequences of RNA isoform changes during neuronal differentiation

Static gene expression programs have been extensively characterized in stem cells and mature human cells. However, the dynamics of RNA isoform changes upon cell-state-transitions during cell differentiation, the determinants and functional consequences have largely remained unclear. Here, we established an improved model for human neurogenesis in vitro that is amenable for systems-wide analyses of gene expression. Our multi-omics analysis reveals that the pronounced alterations in cell morphology correlate strongly with widespread changes in RNA isoform expression. Our approach identifies thousands of new RNA isoforms that are expressed at distinct differentiation stages. RNA isoforms mainly arise from exon skipping and the alternative usage of transcription start and polyadenylation sites during human neurogenesis. The transcript isoform changes can remodel the identity and functions of protein isoforms. Finally, our study identifies a set of RNA binding proteins as a potential determinant of differentiation stage-specific global isoform changes. This work supports the view of regulated isoform changes that underlie state-transitions during neurogenesis.

linc-ADAIN, a human adipose lincRNA, regulates adipogenesis by modulating KLF5 and IL-8 mRNA stability

Adipose tissue remodeling and dysfunction, characterized by elevated inflammation and insulin resistance, play a central role in obesity-related development of type 2 diabetes (T2D) and cardiovascular diseases. Long intergenic non-coding RNAs (lincRNAs) are important regulators of cellular functions. Here, we describe the functions of linc-ADAIN (adipose anti-inflammatory), an adipose lincRNA that is downregulated in white adipose tissue of obese humans. We demonstrate that linc-ADAIN knockdown (KD) increases KLF5 and interleukin-8 (IL-8) mRNA stability and translation by interacting with IGF2BP2. Upregulation of KLF5 and IL-8, via linc-ADAIN KD, leads to an enhanced adipogenic program and adipose tissue inflammation, mirroring the obese state, in vitro and in vivo. KD of linc-ADAIN in human adipose stromal cell (ASC) hTERT adipocytes implanted into mice increases adipocyte size and macrophage infiltration compared to implanted control adipocytes, mimicking hallmark features of obesity-induced adipose tissue remodeling. linc-ADAIN is an anti-inflammatory lincRNA that limits adipose tissue expansion and lipid storage.

Enhancer switching in cell lineage priming is linked to eRNA, Brg1's AT-hook, and SWI/SNF recruitment

RNA transcribed from enhancers, i.e., eRNA, has been suggested to directly activate transcription by recruiting transcription factors and co-activators. Although there have been specific examples of eRNA functioning in this way, it is not clear how general this may be. We find that the AT-hook of SWI/SNF preferentially binds RNA and, as part of the esBAF complex, associates with eRNA transcribed from intronic and intergenic regions. Our data suggest that SWI/SNF is globally recruited in cis by eRNA to cell-type-specific enhancers, representative of two distinct stages that mimic early mammalian development, and not at enhancers that are shared between the two stages. In this manner, SWI/SNF facilitates recruitment and/or activation of MLL3/4, p300/CBP, and Mediator to stage-specific enhancers and super-enhancers that regulate the transcription of metabolic and cell lineage priming-related genes. These findings highlight a connection between ATP-dependent chromatin remodeling and eRNA in cell identity and typical- and super-enhancer activation.

Helper T cell immunity in humans with inherited CD4 deficiency

CD4+ T cells are vital for host defense and immune regulation. However, the fundamental role of CD4 itself remains enigmatic. We report seven patients aged 5-61 years from five families of four ancestries with autosomal recessive CD4 deficiency and a range of infections, including recalcitrant warts and Whipple's disease. All patients are homozygous for rare deleterious CD4 variants impacting expression of the canonical CD4 isoform. A shorter expressed isoform that interacts with LCK, but not HLA class II, is affected by only one variant. All patients lack CD4+ T cells and have increased numbers of TCRαβ+CD4-CD8- T cells, which phenotypically and transcriptionally resemble conventional Th cells. Finally, patient CD4-CD8- αβ T cells exhibit intact responses to HLA class II-restricted antigens and promote B cell differentiation in vitro. Thus, compensatory development of Th cells enables patients with inherited CD4 deficiency to acquire effective cellular and humoral immunity against an unexpectedly large range of pathogens. Nevertheless, CD4 is indispensable for protective immunity against at least human papillomaviruses and Trophyrema whipplei.

Reproducible single cell annotation of programs underlying T-cell subsets, activation states, and functions

T-cells recognize antigens and induce specialized gene expression programs (GEPs) enabling functions including proliferation, cytotoxicity, and cytokine production. Traditionally, different classes of helper T-cells express mutually exclusive responses - for example, Th1, Th2, and Th17 programs. However, new single-cell RNA sequencing (scRNA-Seq) experiments have revealed a continuum of T-cell states without discrete clusters corresponding to these subsets, implying the need for new analytical frameworks. Here, we advance the characterization of T-cells with T-CellAnnoTator (TCAT), a pipeline that simultaneously quantifies pre-defined GEPs capturing activation states and cellular subsets. From 1,700,000 T-cells from 700 individuals across 38 tissues and five diverse disease contexts, we discover 46 reproducible GEPs reflecting the known core functions of T-cells including proliferation, cytotoxicity, exhaustion, and T helper effector states. We experimentally characterize several novel activation programs and apply TCAT to describe T-cell activation and exhaustion in Covid-19 and cancer, providing insight into T-cell function in these diseases.

Genome organization and botanical diversity

The rich diversity of angiosperms, both the planet's dominant flora and the cornerstone of agriculture, is integrally intertwined with a distinctive evolutionary history. Here, we explore the interplay between angiosperm genome organization and botanical diversity, empowered by genomic approaches ranging from genetic linkage mapping to analysis of gene regulation. Commonality in the genetic hardware of plants has enabled robust comparative genomics that has provided a broad picture of angiosperm evolution and implicated both general processes and specific elements in contributing to botanical diversity. We argue that the hardware of plant genomes-both in content and in dynamics-has been shaped by selection for rather substantial differences in gene regulation between plants and animals such as maize and human, organisms of comparable genome size and gene number. Their distinctive genome content and dynamics may reflect in part the indeterminate development of plants that puts strikingly different demands on gene regulation than in animals. Repeated polyploidization of plant genomes and multiplication of individual genes together with extensive rearrangement and differential retention provide rich raw material for selection of morphological and/or physiological variations conferring fitness in specific niches, whether natural or artificial. These findings exemplify the burgeoning information available to employ in increasing knowledge of plant biology and in modifying selected plants to better meet human needs.

Msx genes delineate a novel molecular map of the developing cerebellar neuroepithelium

In the early cerebellar primordium, there are two progenitor zones, the ventricular zone (VZ) residing atop the IVth ventricle and the rhombic lip (RL) at the lateral edges of the developing cerebellum. These zones give rise to the several cell types that form the GABAergic and glutamatergic populations of the adult cerebellum, respectively. Recently, an understanding of the molecular compartmentation of these zones has emerged. To add to this knowledge base, we report on the Msx genes, a family of three transcription factors, that are expressed downstream of Bone Morphogenetic Protein (BMP) signaling in these zones. Using fluorescent RNA in situ hybridization, we have characterized the Msx (Msh Homeobox) genes and demonstrated that their spatiotemporal pattern segregates specific regions within the progenitor zones. Msx1 and Msx2 are compartmentalized within the rhombic lip (RL), while Msx3 is localized within the ventricular zone (VZ). The relationship of the Msx genes with an early marker of the glutamatergic lineage, Atoh1, was examined in Atoh1-null mice and it was found that the expression of Msx genes persisted. Importantly, the spatial expression of Msx1 and Msx3 altered in response to the elimination of Atoh1. These results point to the Msx genes as novel early markers of cerebellar progenitor zones and more importantly to an updated view of the molecular parcellation of the RL with respect to the canonical marker of the RL, Atoh1.

OVCH1 Antisense RNA 1 is differentially expressed between non-frail and frail old adults

While some old adults stay healthy and non-frail up to late in life, others experience multimorbidity and frailty often accompanied by a pro-inflammatory state. The underlying molecular mechanisms for those differences are still obscure. Here, we used gene expression analysis to understand the molecular underpinning between non-frail and frail individuals in old age. Twenty-four adults (50% non-frail and 50% frail) from InCHIANTI study were included. Total RNA extracted from whole blood was analyzed by Cap Analysis of Gene Expression (CAGE). CAGE identified transcription start site (TSS) and active enhancer regions. We identified a set of differentially expressed (DE) TSS and enhancer between non-frail and frail and male and female participants. Several DE TSSs were annotated as lncRNA (XIST and TTTY14) and antisense RNAs (ZFX-AS1 and OVCH1 Antisense RNA 1). The promoter region chr6:366,786,54-366,787,97;+ was DE and overlapping the longevity CDKN1A gene. GWAS-LD enrichment analysis identifies overlapping LD-blocks with the DE regions with reported traits in GWAS catalog (isovolumetric relaxation time and urinary tract infection frequency). Furthermore, we used weighted gene co-expression network analysis (WGCNA) to identify changes of gene expression associated with clinical traits and identify key gene modules. We performed functional enrichment analysis of the gene modules with significant trait/module correlation. One gene module is showing a very distinct pattern in hub genes. Glycogen Phosphorylase L (PYGL) was the top ranked hub gene between non-frail and frail. We predicted transcription factor binding sites (TFBS) and motif activity. TF involved in age-related pathways (e.g., FOXO3 and MYC) shows different expression patterns between non-frail and frail participants. Expanding the study of OVCH1 Antisense RNA 1 and PYGL may help understand the mechanisms leading to loss of homeostasis that ultimately causes frailty.

Assembling the RNA therapeutics toolbox

From the approval of COVID-19 mRNA vaccines to the 2023 Nobel Prize awarded for nucleoside base modifications, RNA therapeutics have entered the spotlight and are transforming drug development. While the term "RNA therapeutics" has been used in various contexts, this review focuses on treatments that utilize RNA as a component or target RNA for therapeutic effects. We summarize the latest advances in RNA-targeting tools and RNA-based technologies, including but not limited to mRNA, antisense oligos, siRNAs, small molecules and RNA editors. We focus on the mechanisms of current FDA-approved therapeutics but also provide a discussion on the upcoming workforces. The clinical utility of RNA-based therapeutics is enabled not only by the advances in RNA technologies but in conjunction with the significant improvements in chemical modifications and delivery platforms, which are also briefly discussed in the review. We summarize the latest RNA therapeutics based on their mechanisms and therapeutic effects, which include expressing proteins for vaccination and protein replacement therapies, degrading deleterious RNA, modulating transcription and translation efficiency, targeting noncoding RNAs, binding and modulating protein activity and editing RNA sequences and modifications. This review emphasizes the concept of an RNA therapeutic toolbox, pinpointing the readers to all the tools available for their desired research and clinical goals. As the field advances, the catalog of RNA therapeutic tools continues to grow, further allowing researchers to combine appropriate RNA technologies with suitable chemical modifications and delivery platforms to develop therapeutics tailored to their specific clinical challenges.

Are non-protein coding RNAs junk or treasure?: An attempt to explain and reconcile opposing viewpoints of whether the human genome is mostly transcribed into non-functional or functional RNAs

The human genome project's lasting legacies are the emerging insights into human physiology and disease, and the ascendance of biology as the dominant science of the 21st century. Sequencing revealed that >90% of the human genome is not coding for proteins, as originally thought, but rather is overwhelmingly transcribed into non-protein coding, or non-coding, RNAs (ncRNAs). This discovery initially led to the hypothesis that most genomic DNA is "junk", a term still championed by some geneticists and evolutionary biologists. In contrast, molecular biologists and biochemists studying the vast number of transcripts produced from most of this genome "junk" often surmise that these ncRNAs have biological significance. What gives? This essay contrasts the two opposing, extant viewpoints, aiming to explain their bases, which arise from distinct reference frames of the underlying scientific disciplines. Finally, it aims to reconcile these divergent mindsets in hopes of stimulating synergy between scientific fields.

An Emerging Role for Enhancer RNAs in Brain Disorders

Noncoding DNA undergoes widespread context-dependent transcription to produce noncoding RNAs. In recent decades, tremendous advances in genomics and transcriptomics have revealed important regulatory roles for noncoding DNA elements and the RNAs that they produce. Enhancers are one such element that are well-established drivers of gene expression changes in response to a variety of factors such as external stimuli, cellular responses, developmental cues, and disease states. They are known to act at long distances, interact with multiple target gene loci simultaneously, synergize with other enhancers, and associate with dynamic chromatin architectures to form a complex regulatory network. Recent advances in enhancer biology have revealed that upon activation, enhancers transcribe long noncoding RNAs, known as enhancer RNAs (eRNAs), that have been shown to play important roles in enhancer-mediated gene regulation and chromatin-modifying activities. In the brain, enhancer dysregulation and eRNA transcription has been reported in numerous disorders from acute injuries to chronic neurodegeneration. Because this is an emerging area, a comprehensive understanding of eRNA function has not yet been achieved in brain disorders; however, the findings to date have illuminated a role for eRNAs in activity-driven gene expression and phenotypic outcomes. In this review, we highlight the breadth of the current literature on eRNA biology in brain health and disease and discuss the challenges as well as focus areas and strategies for future in-depth research on eRNAs in brain health and disease.

Genetic regulatory effects in response to a high-cholesterol, high-fat diet in baboons

Steady-state expression quantitative trait loci (eQTLs) explain only a fraction of disease-associated loci identified through genome-wide association studies (GWASs), while eQTLs involved in gene-by-environment (GxE) interactions have rarely been characterized in humans due to experimental challenges. Using a baboon model, we found hundreds of eQTLs that emerge in adipose, liver, and muscle after prolonged exposure to high dietary fat and cholesterol. Diet-responsive eQTLs exhibit genomic localization and genic features that are distinct from steady-state eQTLs. Furthermore, the human orthologs associated with diet-responsive eQTLs are enriched for GWAS genes associated with human metabolic traits, suggesting that context-responsive eQTLs with more complex regulatory effects are likely to explain GWAS hits that do not seem to overlap with standard eQTLs. Our results highlight the complexity of genetic regulatory effects and the potential of eQTLs with disease-relevant GxE interactions in enhancing the understanding of GWAS signals for human complex disease using non-human primate models.

KLF7 regulates super-enhancer-driven IGF2BP2 overexpression to promote the progression of head and neck squamous cell carcinoma

BACKGROUND

Head and neck squamous carcinoma (HNSCC) is known for its high aggressiveness and susceptibility to cervical lymph node metastasis, which greatly contributes to its poor prognosis. During tumorigenesis, many types of cancer cells acquire oncogenic super-enhancers (SEs) that drive the overexpression of oncogenes, thereby maintaining malignant progression. This study aimed to identify and validate the role of oncogenic SE-associated genes in the malignant progression of HNSCC.

METHODS

We identified HNSCC cell-specific SE-associated genes through H3K27Ac ChIP-seq and overlapped them with HNSCC-associated genes obtained from The Cancer Genome Atlas (TCGA) dataset and Gene Expression Omnibus (GEO) datasets using weighted gene coexpression network analysis (WGCNA) to identify hub genes. The expression of IGF2BP2 and KLF7 in HNSCC was detected using clinical samples. To determine the biological role of IGF2BP2, we performed CCK-8, colony formation assay, Transwell migration assay, invasion assay, and orthotopic xenograft model experiments. Furthermore, we utilized a CRISPR/Cas9 gene-editing system, small-molecule inhibitors, ChIP-qPCR, and dual-luciferase reporter assays to investigate the molecular mechanisms of IGF2BP2 and its upstream transcription factors.

RESULTS

Our study identified IGF2BP2 as a hub SE-associated gene that exhibited aberrant expression in HNSCC tissues. Increased expression of IGF2BP2 was observed to be linked with malignant progression and unfavorable prognosis in HNSCC patients. Both in vitro and in vivo experiments confirmed that IGF2BP2 promotes the tumorigenicity and metastasis of HNSCC by promoting cell proliferation, migration, and invasion. Mechanistically, the IGF2BP2-SE region displayed enrichment for H3K27Ac, BRD4, and MED1, which led to the inhibition of IGF2BP2 transcription and expression through deactivation of the SE-associated transcriptional program. Additionally, KLF7 was found to induce the transcription of IGF2BP2 and directly bind to its promoter and SE regions. Moreover, the abundance of KLF7 exhibited a positive correlation with the abundance of IGF2BP2 in HNSCC. Patients with high expression of both KLF7 and IGF2BP2 showed poorer prognosis. Lastly, we demonstrated that the small molecule inhibitor JQ1, targeting BRD4, attenuated the proliferation and metastatic abilities of HNSCC cells.

CONCLUSIONS

Our study reveals the critical role of IGF2BP2 overexpression mediated by SE and KLF7 in promoting HNSCC progression. Targeting SE-associated transcriptional programs may represent a potential therapeutic strategy in managing HNSCC.

The relationship between extreme inter-individual variation in macrophage gene expression and genetic susceptibility to inflammatory bowel disease

The differentiation of resident intestinal macrophages from blood monocytes depends upon signals from the macrophage colony-stimulating factor receptor (CSF1R). Analysis of genome-wide association studies (GWAS) indicates that dysregulation of macrophage differentiation and response to microorganisms contributes to susceptibility to chronic inflammatory bowel disease (IBD). Here, we analyzed transcriptomic variation in monocyte-derived macrophages (MDM) from affected and unaffected sib pairs/trios from 22 IBD families and 6 healthy controls. Transcriptional network analysis of the data revealed no overall or inter-sib distinction between affected and unaffected individuals in basal gene expression or the temporal response to lipopolysaccharide (LPS). However, the basal or LPS-inducible expression of individual genes varied independently by as much as 100-fold between subjects. Extreme independent variation in the expression of pairs of HLA-associated transcripts (HLA-B/C, HLA-A/F and HLA-DRB1/DRB5) in macrophages was associated with HLA genotype. Correlation analysis indicated the downstream impacts of variation in the immediate early response to LPS. For example, variation in early expression of IL1B was significantly associated with local SNV genotype and with subsequent peak expression of target genes including IL23A, CXCL1, CXCL3, CXCL8 and NLRP3. Similarly, variation in early IFNB1 expression was correlated with subsequent expression of IFN target genes. Our results support the view that gene-specific dysregulation in macrophage adaptation to the intestinal milieu is associated with genetic susceptibility to IBD.

Pig-eRNAdb: a comprehensive enhancer and eRNA dataset of pigs

Enhancers and the enhancer RNAs (eRNAs) have been strongly implicated in regulations of transcriptions. Based the multi-omics data (ATAC-seq, ChIP-seq and RNA-seq) from public databases, Pig-eRNAdb is a dataset that comprehensively integrates enhancers and eRNAs for pigs using the machine learning strategy, which incorporates 82,399 enhancers and 37,803 eRNAs from 607 samples across 15 tissues of pigs. This user-friendly dataset covers a comprehensive depth of enhancers and eRNAs annotation for pigs. The coordinates of enhancers and the expression patterns of eRNAs are downloadable. Besides, thousands of regulators on eRNAs, the target genes of eRNAs, the tissue-specific eRNAs, and the housekeeping eRNAs are also accessible as well as the sequence similarity of eRNAs with humans. Moreover, the tissue-specific eRNA-trait associations encompass 652 traits are also provided. It will crucially facilitate investigations on enhancers and eRNAs with Pig-eRNAdb as a reference dataset in pigs.

Molecular profiling of high-level athlete skeletal muscle after acute endurance or resistance exercise - A systems biology approach

OBJECTIVE

Long-term high-level exercise training leads to improvements in physical performance and multi-tissue adaptation following changes in molecular pathways. While skeletal muscle baseline differences between exercise-trained and untrained individuals have been previously investigated, it remains unclear how training history influences human multi-omics responses to acute exercise.

METHODS

We recruited and extensively characterized 24 individuals categorized as endurance athletes with >15 years of training history, strength athletes or control subjects. Timeseries skeletal muscle biopsies were taken from M. vastus lateralis at three time-points after endurance or resistance exercise was performed and multi-omics molecular analysis performed.

RESULTS

Our analyses revealed distinct activation differences of molecular processes such as fatty- and amino acid metabolism and transcription factors such as HIF1A and the MYF-family. We show that endurance athletes have an increased abundance of carnitine-derivates while strength athletes increase specific phospholipid metabolites compared to control subjects. Additionally, for the first time, we show the metabolite sorbitol to be substantially increased with acute exercise. On transcriptional level, we show that acute resistance exercise stimulates more gene expression than acute endurance exercise. This follows a specific pattern, with endurance athletes uniquely down-regulating pathways related to mitochondria, translation and ribosomes. Finally, both forms of exercise training specialize in diverging transcriptional directions, differentiating themselves from the transcriptome of the untrained control group.

CONCLUSIONS

We identify a "transcriptional specialization effect" by transcriptional narrowing and intensification, and molecular specialization effects on metabolomic level Additionally, we performed multi-omics network and cluster analysis, providing a novel resource of skeletal muscle transcriptomic and metabolomic profiling in highly trained and untrained individuals.

Stress-induced nuclear speckle reorganization is linked to activation of immediate early gene splicing

Current models posit that nuclear speckles (NSs) serve as reservoirs of splicing factors and facilitate posttranscriptional mRNA processing. Here, we discovered that ribotoxic stress induces a profound reorganization of NSs with enhanced recruitment of factors required for splice-site recognition, including the RNA-binding protein TIAR, U1 snRNP proteins and U2-associated factor 65, as well as serine 2 phosphorylated RNA polymerase II. NS reorganization relies on the stress-activated p38 mitogen-activated protein kinase (MAPK) pathway and coincides with splicing activation of both pre-existing and newly synthesized pre-mRNAs. In particular, ribotoxic stress causes targeted excision of retained introns from pre-mRNAs of immediate early genes (IEGs), whose transcription is induced during the stress response. Importantly, enhanced splicing of the IEGs ZFP36 and FOS is accompanied by relocalization of the corresponding nuclear mRNA foci to NSs. Our study reveals NSs as a dynamic compartment that is remodeled under stress conditions, whereby NSs appear to become sites of IEG transcription and efficient cotranscriptional splicing.

The chromatin signatures of enhancers and their dynamic regulation

Enhancers are cis-regulatory elements that can stimulate gene expression from distance, and drive precise spatiotemporal gene expression profiles during development. Functional enhancers display specific features including an open chromatin conformation, Histone H3 lysine 27 acetylation, Histone H3 lysine 4 mono-methylation enrichment, and enhancer RNAs production. These features are modified upon developmental cues which impacts their activity. In this review, we describe the current state of knowledge about enhancer functions and the diverse chromatin signatures found on enhancers. We also discuss the dynamic changes of enhancer chromatin signatures, and their impact on lineage specific gene expression profiles, during development or cellular differentiation.

Genetic regulatory effects in response to a high cholesterol, high fat diet in baboons

Steady-state expression quantitative trait loci (eQTLs) explain only a fraction of disease-associated loci identified through genome-wide association studies (GWAS), while eQTLs involved in gene-by-environment (GxE) interactions have rarely been characterized in humans due to experimental challenges. Using a baboon model, we found hundreds of eQTLs that emerge in adipose, liver, and muscle after prolonged exposure to high dietary fat and cholesterol. Diet-responsive eQTLs exhibit genomic localization and genic features that are distinct from steady-state eQTLs. Furthermore, the human orthologs associated with diet-responsive eQTLs are enriched for GWAS genes associated with human metabolic traits, suggesting that context-responsive eQTLs with more complex regulatory effects are likely to explain GWAS hits that do not seem to overlap with standard eQTLs. Our results highlight the complexity of genetic regulatory effects and the potential of eQTLs with disease-relevant GxE interactions in enhancing the understanding of GWAS signals for human complex disease using nonhuman primate models.

Transposable element-initiated enhancer-like elements generate the subgenome-biased spike specificity of polyploid wheat

Transposable elements (TEs) comprise ~85% of the common wheat genome, which are highly diverse among subgenomes, possibly contribute to polyploid plasticity, but the causality is only assumed. Here, by integrating data from gene expression cap analysis and epigenome profiling via hidden Markov model in common wheat, we detect a large proportion of enhancer-like elements (ELEs) derived from TEs producing nascent noncoding transcripts, namely ELE-RNAs, which are well indicative of the regulatory activity of ELEs. Quantifying ELE-RNA transcriptome across typical developmental stages reveals that TE-initiated ELE-RNAs are mainly from RLG_famc7.3 specifically expanded in subgenome A. Acquisition of spike-specific transcription factor binding likely confers spike-specific expression of RLG_famc7.3-initiated ELE-RNAs. Knockdown of RLG_famc7.3-initiated ELE-RNAs resulted in global downregulation of spike-specific genes and abnormal spike development. These findings link TE expansion to regulatory specificity and polyploid developmental plasticity, highlighting the functional impact of TE-driven regulatory innovation on polyploid evolution.

Enhancer target prediction: state-of-the-art approaches and future prospects

Enhancers are genomic regions that regulate gene transcription and are located far away from the transcription start sites of their target genes. Enhancers are highly enriched in disease-associated variants and thus deciphering the interactions between enhancers and genes is crucial to understanding the molecular basis of genetic predispositions to diseases. Experimental validations of enhancer targets can be laborious. Computational methods have thus emerged as a valuable alternative for studying enhancer-gene interactions. A variety of computational methods have been developed to predict enhancer targets by incorporating genomic features (e.g. conservation, distance, and sequence), epigenomic features (e.g. histone marks and chromatin contacts) and activity measurements (e.g. covariations of enhancer activity and gene expression). With the recent advances in genome perturbation and chromatin conformation capture technologies, data on experimentally validated enhancer targets are becoming available for supervised training of these methods and evaluation of their performance. In this review, we categorize enhancer target prediction methods based on their rationales and approaches. Then we discuss their merits and limitations and highlight the future directions for enhancer targets prediction.

ZEB2 and MEIS1 independently contribute to hematopoiesis via early hematopoietic enhancer activation

Cell differentiation is achieved by acquiring a cell type-specific transcriptional program and epigenetic landscape. While the cell type-specific patterning of enhancers has been shown to precede cell fate decisions, it remains unclear how regulators of these enhancers are induced to initiate cell specification and how they appropriately restrict cells that differentiate. Here, using embryonic stem cell-derived hematopoietic cell differentiation cultures, we show the activation of some hematopoietic enhancers during arterialization of hemogenic endothelium, a prerequisite for hematopoiesis. We further reveal that ZEB2, a factor involved in the transcriptional regulation of arterial endothelial cells, and a hematopoietic regulator MEIS1 are independently required for activating these enhancers. Concomitantly, ZEB2 or MEIS1 deficiency impaired hematopoietic cell development. These results suggest that multiple regulators expressed from an earlier developmental stage non-redundantly contribute to the establishment of hematopoietic enhancer landscape, thereby restricting cell differentiation despite the unrestricted expression of these regulators to hematopoietic cells.

Genome-wide measurement of RNA dissociation from chromatin classifies transcripts by their dynamics and reveals rapid dissociation of enhancer lncRNAs

Long non-coding RNAs (lncRNAs) are involved in gene expression regulation in cis. Although enriched in the cell chromatin fraction, to what degree this defines their regulatory potential remains unclear. Furthermore, the factors underlying lncRNA chromatin tethering, as well as the molecular basis of efficient lncRNA chromatin dissociation and its impact on enhancer activity and target gene expression, remain to be resolved. Here, we developed chrTT-seq, which combines the pulse-chase metabolic labeling of nascent RNA with chromatin fractionation and transient transcriptome sequencing to follow nascent RNA transcripts from their transcription on chromatin to release and allows the quantification of dissociation dynamics. By incorporating genomic, transcriptomic, and epigenetic metrics, as well as RNA-binding protein propensities, in machine learning models, we identify features that define transcript groups of different chromatin dissociation dynamics. Notably, lncRNAs transcribed from enhancers display reduced chromatin retention, suggesting that, in addition to splicing, their chromatin dissociation may shape enhancer activity.

Understanding blood development and leukemia using sequencing-based technologies and human cell systems

Our current understanding of human hematopoiesis has undergone significant transformation throughout the years, challenging conventional views. The evolution of high-throughput technologies has enabled the accumulation of diverse data types, offering new avenues for investigating key regulatory processes in blood cell production and disease. In this review, we will explore the opportunities presented by these advancements for unraveling the molecular mechanisms underlying normal and abnormal hematopoiesis. Specifically, we will focus on the importance of enhancer-associated regulatory networks and highlight the crucial role of enhancer-derived transcription regulation. Additionally, we will discuss the unprecedented power of single-cell methods and the progression in using in vitro human blood differentiation system, in particular induced pluripotent stem cell models, in dissecting hematopoietic processes. Furthermore, we will explore the potential of ever more nuanced patient profiling to allow precision medicine approaches. Ultimately, we advocate for a multiparameter, regulatory network-based approach for providing a more holistic understanding of normal hematopoiesis and blood disorders.

Mapping the epigenomic landscape of human monocytes following innate immune activation reveals context-specific mechanisms driving endotoxin tolerance

BACKGROUND

Monocytes are key mediators of innate immunity to infection, undergoing profound and dynamic changes in epigenetic state and immune function which are broadly protective but may be dysregulated in disease. Here, we aimed to advance understanding of epigenetic regulation following innate immune activation, acutely and in endotoxin tolerant states.

METHODS

We exposed human primary monocytes from healthy donors (n = 6) to interferon-γ or differing combinations of endotoxin (lipopolysaccharide), including acute response (2 h) and two models of endotoxin tolerance: repeated stimulations (6 + 6 h) and prolonged exposure to endotoxin (24 h). Another subset of monocytes was left untreated (naïve). We identified context-specific regulatory elements based on epigenetic signatures for chromatin accessibility (ATAC-seq) and regulatory non-coding RNAs from total RNA sequencing.

RESULTS

We present an atlas of differential gene expression for endotoxin and interferon response, identifying widespread context specific changes. Across assayed states, only 24-29% of genes showing differential exon usage are also differential at the gene level. Overall, 19.9% (6,884 of 34,616) of repeatedly observed ATAC peaks were differential in at least one condition, the majority upregulated on stimulation and located in distal regions (64.1% vs 45.9% of non-differential peaks) within which sequences were less conserved than non-differential peaks. We identified enhancer-derived RNA signatures specific to different monocyte states that correlated with chromatin accessibility changes. The endotoxin tolerance models showed distinct chromatin accessibility and transcriptomic signatures, with integrated analysis identifying genes and pathways involved in the inflammatory response, detoxification, metabolism and wound healing. We leveraged eQTL mapping for the same monocyte activation states to link potential enhancers with specific genes, identifying 1,946 unique differential ATAC peaks with 1,340 expression associated genes. We further use this to inform understanding of reported GWAS, for example involving FCHO1 and coronary artery disease.

CONCLUSION

This study reports context-specific regulatory elements based on transcriptomic profiling and epigenetic signatures for enhancer-derived RNAs and chromatin accessibility in immune tolerant monocyte states, and demonstrates the informativeness of linking such elements and eQTL to inform future mechanistic studies aimed at defining therapeutic targets of immunosuppression and diseases.

A body map of super-enhancers and their function in pig

Introduction

Super-enhancers (SEs) are clusters of enhancers that act synergistically to drive the high-level expression of genes involved in cell identity and function. Although SEs have been extensively investigated in humans and mice, they have not been well characterized in pigs.

Methods

Here, we identified 42,380 SEs in 14 pig tissues using chromatin immunoprecipitation sequencing, and statistics of its overall situation, studied the composition and characteristics of SE, and explored the influence of SEs characteristics on gene expression.

Results

We observed that approximately 40% of normal enhancers (NEs) form SEs. Compared to NEs, we found that SEs were more likely to be enriched with an activated enhancer and show activated functions. Interestingly, SEs showed X chromosome depletion and short interspersed nuclear element enrichment, implying that SEs play an important role in sex traits and repeat evolution. Additionally, SE-associated genes exhibited higher expression levels and stronger conservation than NE-associated genes. However, genes with the largest SEs had higher expression levels than those with the smallest SEs, indicating that SE size may influence gene expression. Moreover, we observed a negative correlation between SE gene distance and gene expression, indicating that the proximity of SEs can affect gene activity. Gene ontology enrichment and motif analysis revealed that SEs have strong tissue-specific activity. For example, the CORO2B gene with a brain-specific SE shows strong brain-specific expression, and the phenylalanine hydroxylase gene with liver-specific SEs shows strong liver-specific expression.

Discussion

In this study, we illustrated a body map of SEs and explored their functions in pigs, providing information on the composition and tissue-specific patterns of SEs. This study can serve as a valuable resource of gene regulatory and comparative analyses to the scientific community and provides a theoretical reference for genetic control mechanisms of important traits in pigs.

Liver endothelial cells in NAFLD and transition to NASH and HCC

Non-alcoholic fatty liver disease (NAFLD) is considered as the hepatic manifestation of metabolic syndrome, which is characterised by obesity, insulin resistance, hypercholesterolemia and hypertension. NAFLD is the most frequent liver disease worldwide and more than 10% of NAFLD patients progress to the inflammatory and fibrotic stage of non-alcoholic steatohepatitis (NASH), which can lead to end-stage liver disease including hepatocellular carcinoma (HCC), the most frequent primary malignant liver tumor. Liver sinusoidal endothelial cells (LSEC) are strategically positioned at the interface between blood and hepatic parenchyma. LSECs are highly specialized cells, characterised by the presence of transcellular pores, called fenestrae, and exhibit anti-inflammatory and anti-fibrotic characteristics under physiological conditions. However, during NAFLD development they undergo capillarisation and acquire a phenotype similar to vascular endothelial cells, actively promoting all pathophysiological aspects of NAFLD, including steatosis, inflammation, and fibrosis. LSEC dysfunction is critical for the progression to NASH and HCC while restoring LSEC homeostasis appears to be a promising approach to prevent NAFLD progression and its complications and even reverse tissue damage. In this review we present current information on the role of LSEC throughout the progressive phases of NAFLD, summarising in vitro and in vivo experimental evidence and data from human studies.

Transcriptional signature of islet neogenesis-associated protein peptide-treated rat pancreatic islets reveals induction of novel long non-coding RNAs

Background

Diabetes mellitus is characterized by chronic hyperglycemia with loss of β-cell function and mass. An attractive therapeutic approach to treat patients with diabetes in a non-invasive way is to harness the innate regenerative potential of the pancreas. The Islet Neogenesis-Associated Protein pentadecapeptide (INGAP-PP) has been shown to induce β-cell regeneration and improve their function in rodents. To investigate its possible mechanism of action, we report here the global transcriptional effects induced by the short-term INGAP-PP in vitro treatment of adult rat pancreatic islets.

Methods and findings

Rat pancreatic islets were cultured in vitro in the presence of INGAP-PP for 4 days, and RNA-seq was generated from triplicate treated and control islet samples. We performed a de novo rat gene annotation based on the alignment of RNA-seq reads. The list of INGAP-PP-regulated genes was integrated with epigenomic data. Using the new gene annotation generated in this work, we quantified RNA-seq data profiled in INS-1 cells treated with IL1β, IL1β+Calcipotriol (a vitamin D agonist) or vehicle, and single-cell RNA-seq data profiled in rat pancreatic islets. We found 1,669 differentially expressed genes by INGAP-PP treatment, including dozens of previously unannotated rat transcripts. Genes differentially expressed by the INGAP-PP treatment included a subset of upregulated transcripts that are associated with vitamin D receptor activation. Supported by epigenomic and single-cell RNA-seq data, we identified 9 previously unannotated long noncoding RNAs (lncRNAs) upregulated by INGAP-PP, some of which are also differentially regulated by IL1β and vitamin D in β-cells. These include Ri-lnc1, which is enriched in mature β-cells.

Conclusions

Our results reveal the transcriptional program that could explain the enhancement of INGAP-PP-mediated physiological effects on β-cell mass and function. We identified novel lncRNAs that are induced by INGAP-PP in rat islets, some of which are selectively expressed in pancreatic β-cells and downregulated by IL1β treatment of INS-1 cells. Our results suggest a relevant function for Ri-lnc1 in β-cells. These findings are expected to provide the basis for a deeper understanding of islet translational results from rodents to humans, with the ultimate goal of designing new therapies for people with diabetes.