Transcription factor antagonism regulates heterogeneity in embryonic stem cell states.

Gene expression heterogeneity underlies cell states and contributes to developmental robustness. While heterogeneity can arise from stochastic transcriptional processes, the extent to which it is regulated is unclear. Here, we characterize the regulatory program underlying heterogeneity in murine embryonic stem cell (mESC) states. We identify differentially active and transcribed enhancers (DATEs) across states. DATEs regulate differentially expressed genes and are distinguished by co-binding of transcription factors Klf4 and Zfp281. In contrast to other factors that interact in a positive feedback network stabilizing mESC cell-type identity, Klf4 and Zfp281 drive opposing transcriptional and chromatin programs. Abrogation of factor binding to DATEs dampens variation in gene expression, and factor loss alters kinetics of switching between states. These results show antagonism between factors at enhancers results in gene expression heterogeneity and formation of cell states, with implications for the generation of diverse cell types during development.

Systematic functional characterization of antisense eRNA of protocadherin α composite enhancer.

Enhancers generate bidirectional noncoding enhancer RNAs (eRNAs) that may regulate gene expression. At present, the eRNA function remains enigmatic. Here, we report a 5' capped antisense eRNA PEARL (Pcdh eRNA associated with R-loop formation) that is transcribed from the protocadherin (Pcdh) α HS5-1 enhancer region. Through loss- and gain-of-function experiments with CRISPR/Cas9 DNA fragment editing, CRISPRi, and CRISPRa, as well as locked nucleic acid strategies, in conjunction with ChIRP, MeDIP, DRIP, QHR-4C, and HiChIP experiments, we found that PEARL regulates Pcdhα gene expression by forming local RNA-DNA duplexes (R-loops) in situ within the HS5-1 enhancer region to promote long-distance chromatin interactions between distal enhancers and target promoters. In particular, increased levels of eRNA PEARL via perturbing transcription elongation factor SPT6 lead to strengthened local three-dimensional chromatin organization within the Pcdh superTAD. These findings have important implications regarding molecular mechanisms by which the HS5-1 enhancer regulates stochastic Pcdhα promoter choice in single cells in the brain.

The Integrator Complex Attenuates Promoter-Proximal Transcription at Protein-Coding Genes.

The transition of RNA polymerase II (Pol II) from initiation to productive elongation is a central, regulated step in metazoan gene expression. At many genes, Pol II pauses stably in early elongation, remaining engaged with the 25- to 60-nt-long nascent RNA for many minutes while awaiting signals for release into the gene body. However, 15%-20% of genes display highly unstable promoter Pol II, suggesting that paused polymerase might dissociate from template DNA at these promoters and release a short, non-productive mRNA. Here, we report that paused Pol II can be actively destabilized by the Integrator complex. Specifically, we present evidence that Integrator utilizes its RNA endonuclease activity to cleave nascent RNA and drive termination of paused Pol II. These findings uncover a previously unappreciated mechanism of metazoan gene repression, akin to bacterial transcription attenuation, wherein promoter-proximal Pol II is prevented from entering productive elongation through factor-regulated termination.

R-Loops Promote Antisense Transcription across the Mammalian Genome.

Widespread antisense long noncoding RNA (lncRNA) overlap with many protein-coding genes in mammals and emanate from gene promoter, enhancer, and termination regions. However, their origin and biological purpose remain unclear. We show that these antisense lncRNA can be generated by R-loops that form when nascent transcript invades the DNA duplex behind elongating RNA polymerase II (Pol II). Biochemically, R-loops act as intrinsic Pol II promoters to induce de novo RNA synthesis. Furthermore, their removal across the human genome by RNase H1 overexpression causes the selective reduction of antisense transcription. Consequently, we predict that R-loops act to facilitate the synthesis of many gene proximal antisense lncRNA. Not only are R-loops widely associated with DNA damage and repair, but we now show that they have the capacity to promote de novo transcript synthesis that may have aided the evolution of gene regulation.

A Muscle-Specific Enhancer RNA Mediates Cohesin Recruitment and Regulates Transcription In trans.

The enhancer regions of the myogenic master regulator MyoD give rise to at least two enhancer RNAs. Core enhancer eRNA (CEeRNA) regulates transcription of the adjacent MyoD gene, whereas DRReRNA affects expression of Myogenin in trans. We found that DRReRNA is recruited at the Myogenin locus, where it colocalizes with Myogenin nascent transcripts. DRReRNA associates with the cohesin complex, and this association correlates with its transactivating properties. Despite being expressed in undifferentiated cells, cohesin is not loaded on Myogenin until the cells start expressing DRReRNA, which is then required for cohesin chromatin recruitment and maintenance. Functionally, depletion of either cohesin or DRReRNA reduces chromatin accessibility, prevents Myogenin activation, and hinders muscle cell differentiation. Thus, DRReRNA ensures spatially appropriate cohesin loading in trans to regulate gene expression.

JMJD6 Licenses ERα-Dependent Enhancer and Coding Gene Activation by Modulating the Recruitment of the CARM1/MED12 Co-activator Complex.

Whereas the actions of enhancers in gene transcriptional regulation are well established, roles of JmjC-domain-containing proteins in mediating enhancer activation remain poorly understood. Here, we report that recruitment of the JmjC-domain-containing protein 6 (JMJD6) to estrogen receptor alpha (ERα)-bound active enhancers is required for RNA polymerase II recruitment and enhancer RNA production on enhancers, resulting in transcriptional pause release of cognate estrogen target genes. JMJD6 is found to interact with MED12 in the mediator complex to regulate its recruitment. Unexpectedly, JMJD6 is necessary for MED12 to interact with CARM1, which methylates MED12 at multiple arginine sites and regulates its chromatin binding. Consistent with its role in transcriptional activation, JMJD6 is required for estrogen/ERα-induced breast cancer cell growth and tumorigenesis. Our data have uncovered a critical regulator of estrogen/ERα-induced enhancer coding gene activation and breast cancer cell potency, providing a potential therapeutic target of ER-positive breast cancers.

Spirits in the Material World: Enhancer RNAs in Transcriptional Regulation.

Responses to developmental and environmental cues depend on precise spatiotemporal control of gene transcription. Enhancers, which comprise DNA elements bound by regulatory proteins, can activate target genes in response to these external signals. Recent studies have shown that enhancers are transcribed to produce enhancer RNAs (eRNAs). Do eRNAs play a functional role in activating gene expression or are they non-functional byproducts of nearby transcription machinery? The unstable nature of eRNAs and over-reliance on knockdown approaches have made elucidating the possible functions of eRNAs challenging. We focus here on studies using cloned eRNAs to study their function as transcripts, revealing roles for eRNAs in enhancer-promoter looping, recruiting transcriptional machinery, and facilitating RNA polymerase pause-release to regulate gene expression.

Intragenic Enhancers Attenuate Host Gene Expression.

Eukaryotic gene transcription is regulated at many steps, including RNA polymerase II (Pol II) recruitment, transcription initiation, promoter-proximal Pol II pause release, and transcription termination; however, mechanisms regulating transcription during productive elongation remain poorly understood. Enhancers, which activate gene transcription, themselves undergo Pol II-mediated transcription, but our understanding of enhancer transcription and enhancer RNAs (eRNAs) remains incomplete. Here we show that transcription at intragenic enhancers interferes with and attenuates host gene transcription during productive elongation. While the extent of attenuation correlates positively with nascent eRNA expression, the act of intragenic enhancer transcription alone, but not eRNAs, explains the attenuation. Through CRISPR/Cas9-mediated deletions, we demonstrate a physiological role for intragenic enhancer-mediated transcription attenuation in cell fate determination. We propose that intragenic enhancers not only enhance transcription of one or more genes from a distance but also fine-tune transcription of their host gene through transcription interference, facilitating differential utilization of the same regulatory element for disparate functions.

hnRNPA0 promotes MYB expression by interacting with enhancer lncRNA MY34UE-AS in human leukemia cells.

MYB is a key regulator of hematopoiesis and erythropoiesis, and dysregulation of MYB is closely involved in the development of leukemia, however the mechanism of MYB regulation remains still unclear so far. Our previous study identified a long noncoding RNA (lncRNA) derived from the -34 kb enhancer of the MYB locus, which can promote MYB expression, the proliferation and migration of human leukemia cells, and is therefore termed MY34UE-AS. Then the interacting partner proteins of MY34UE-AS were identified and studied in the present study. hnRNPA0 was identified as a binding partner of MY34UE-AS through RNA pulldown assay, which was further validated through RNA immunoprecipitation (RIP). hnRNPA0 interacted with MY34UE-AS mainly through its RRM2 domain. hnRNPA0 overexpression upregulated MYB and increased the proliferation and migration of K562 cells, whereas hnRNPA0 knockdown showed opposite effects. Rescue experiments showed MY34UE-AS was required for above mentioned functions of hnRNPA0. These results reveal that hnRNPA0 is involved in leukemia through upregulating MYB expression by interacting with MY34UE-AS, suggesting that the hnRNPA0/MY34UE-AS axis could serve as a potential target for leukemia treatment.

A mutator-derived prognostic eRNA signature provides insight into the pathogenesis of breast cancer.

Abundant evidence suggests that enhancer RNA (eRNA) is closely related to tumorigenesis, and the role of eRNA transcription in promoting genomic instability in cancers is gradually unveiled. However, research on the evaluation of the prognostic value and molecular mechanisms of genomic instability associated eRNAs in breast cancer is long overdue. Here, we integratively analyzed eRNA expression and somatic mutation profiles in breast cancer genome. We identified genomic instability associated eRNAs and developed a prognostic signature based on these eRNAs with the area under the curve (AUC) around 0.8 at 9-year survival. We further found the prognostic value of this signature is independent of common clinical factors and is better than TP53 status. Higher expression of genomic instability associated genes in the high-risk group was observed, suggesting that this eRNA signature may serve as an indicator of genomic instability in breast cancer. We found prognostic eRNA co-expressed genes are mainly enriched in Gene set 'Breast Cancer 8P12-P11 Amplicon', Gene set 'Metabolism of lipids' and GO process 'Ubiquitin protein ligase binding'. Furthermore, 11 eRNA-signature regulated genes are identified by assessing promoter-enhancer interaction. Among these genes, F11R, BHLHE40, and NECTIN4 are previously reported oncogenes and EGOT is a tumor suppressor gene, indicating the direct roles of eRNAs in tumorigenesis.

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.

CRED9: A differentially expressed elncRNA regulates expression of transcription factor CEBPA.

Enhancer RNAs (eRNA) are non-coding transcripts produced from active enhancers and have potential gene regulatory function. CCAAT enhancer-binding protein alpha (CEBPA) is a transcription factor generally involved in metabolism, cell cycle inhibition, hematopoiesis, adipogenesis, hepatogenesis, and is associated with tumorigenesis. In this study, we demonstrate that an enhancer-associated long non-coding RNA (elncRNA), transcribed from an enhancer located 9kb downstream from the transcriptional start site (TSS) of CEBPA, positively regulates the expression of CEBPA. As a result, we named this elncRNA 'CEBPA regulatory elncRNA downstream 9kb' or 'CRED9'. CRED9 expression level positively correlates with CEBPA mRNA expression across multiple cell lines as detected by RT droplet digital PCR. Knockdown of CRED9 resulted in a reduction of CEBPA mRNA expression in Hep3B cells. Additionally, CRED9 knockdown in Hep3B and HepG2 cells resulted in lower CEBPA protein expression. We also found that knockdown of CRED9 in Hep3B cells caused a 57.8% reduction in H3K27ac levels at the +9kb CEBPA enhancer. H3K27ac has previously been described as a marker of active enhancers. Taken together, the evidence presented here supports a previously proposed model whereby, in some contexts, eRNA transcripts are necessary to amplify and maintain H3K27ac levels at a given enhancer. Ultimately, this study adds to the growing body of evidence that elncRNA transcripts have important roles in enhancer function and gene regulation.

Non-coding RNAs in immunoregulation and autoimmunity: Technological advances and critical limitations.

Immune cell function is critically dependent on precise control over transcriptional output from the genome. In this respect, integration of environmental signals that regulate gene expression, specifically by transcription factors, enhancer DNA elements, genome topography and non-coding RNAs (ncRNAs), are key components. The first three have been extensively investigated. Even though non-coding RNAs represent the vast majority of cellular RNA species, this class of RNA remains historically understudied. This is partly because of a lag in technological and bioinformatic innovations specifically capable of identifying and accurately measuring their expression. Nevertheless, recent progress in this domain has enabled a profusion of publications identifying novel sub-types of ncRNAs and studies directly addressing the function of ncRNAs in human health and disease. Many ncRNAs, including circular and enhancer RNAs, have now been demonstrated to play key functions in the regulation of immune cells and to show associations with immune-mediated diseases. Some ncRNAs may function as biomarkers of disease, aiding in diagnostics and in estimating response to treatment, while others may play a direct role in the pathogenesis of disease. Importantly, some are relatively stable and are amenable to therapeutic targeting, for example through gene therapy. Here, we provide an overview of ncRNAs and review technological advances that enable their study and hold substantial promise for the future. We provide context-specific examples by examining the associations of ncRNAs with four prototypical human autoimmune diseases, specifically rheumatoid arthritis, psoriasis, inflammatory bowel disease and multiple sclerosis. We anticipate that the utility and mechanistic roles of these ncRNAs in autoimmunity will be further elucidated in the near future.

Comprehensive Pan-Cancer Mutation Density Patterns in Enhancer RNA.

Significant advances have been achieved in understanding the critical role of enhancer RNAs (eRNAs) in the complex field of gene regulation. However, notable uncertainty remains concerning the biology of eRNAs, highlighting the need for continued research to uncover their exact functions in cellular processes and diseases. We present a comprehensive study to scrutinize mutation density patterns, mutation strand bias, and mutation burden in eRNAs across multiple cancer types. Our findings reveal that eRNAs exhibit mutation strand bias akin to that observed in protein-coding RNAs. We also identified a novel pattern, in which mutation density is notably diminished around the central region of the eRNA, but conspicuously elevated towards both the beginning and end. This pattern can be potentially explained by a mechanism involving heightened transcriptional activity and the activation of transcription-coupled repair. The central regions of the eRNAs appear to be more conserved, hinting at a potential mechanism preserving their structural and functional integrity, while the extremities may be more susceptible to mutations due to increased exposure. The evolutionary trajectory of this mutational pattern suggests a nuanced adaptation in eRNAs, where stability at their core coexists with flexibility at their extremities, potentially facilitating their diverse interactions with other genetic entities.

Genomic determinants implicated in the glucocorticoid-mediated induction of KLF9 in pulmonary epithelial cells.

Ligand-activated glucocorticoid receptor (GR) elicits variable glucocorticoid-modulated transcriptomes in different cell types. However, some genes, including Krüppel-like factor 9 (KLF9), a putative transcriptional repressor, demonstrate conserved responses. We show that glucocorticoids induce KLF9 expression in the human airways in vivo and in differentiated human bronchial epithelial (HBE) cells grown at air-liquid interface (ALI). In A549 and BEAS-2B pulmonary epithelial cells, glucocorticoids induce KLF9 expression with similar kinetics to primary HBE cells in submersion culture. A549 and BEAS-2B ChIP-seq data reveal four common glucocorticoid-induced GR binding sites (GBSs). Two GBSs mapped to the 5'-proximal region relative to KLF9 transcription start site (TSS) and two occurred at distal sites. These were all confirmed in primary HBE cells. Global run-on (GRO) sequencing indicated robust enhancer RNA (eRNA) production from three of these GBSs in BEAS-2B cells. This was confirmed in A549 cells, plus submersion, and ALI culture of HBE cells. Cloning each GBS into luciferase reporters revealed glucocorticoid-induced activity requiring a glucocorticoid response element (GRE) within each distal GBS. While the proximal GBSs drove modest reporter induction by glucocorticoids, this region exhibited basal eRNA production, RNA polymerase II enrichment, and looping to the TSS, plausibly underlying constitutive KLF9 expression. Post glucocorticoid treatment, interactions between distal and proximal GBSs and the TSS correlated with KLF9 induction. CBP/P300 silencing reduced proximal GBS activity, but negligibly affected KLF9 expression. Overall, a model for glucocorticoid-mediated regulation of KLF9 involving multiple GBSs is depicted. This work unequivocally demonstrates that mechanistic insights gained from cell lines can translate to physiologically relevant systems.

Identification of cell cycle-associated and -unassociated regulators for expression of a hepatocellular carcinoma oncogene cyclin-dependent kinase inhibitor 3.

Human cyclin-dependent kinase inhibitor 3 (CDKN3) is a known oncogene in hepatocellular carcinoma (HCC) and its expression is promoted during tumor development. CDKN3 serves as a cell cycle regulator and its dysregulation is considered to be a causal factor for tumor progression. However, the molecular basis of the regulation of CDKN3 expression remains largely elusive. Using in silico approach, we identified CDKN3SE, a super enhancer (SE), and enhancer RNA (eRNA) candidates transcribed from this SE. Among the eRNA candidates, the expression of CDKN3eRNA was detected in the human HCC model cell line HepG2, and was found to facilitate the expression of CDKN3 without affecting the cell proliferation rate. In silico screening revealed two DNA-binding transcription factors, upstream stimulatory factor (USF) 1 and 2, involved in the regulation of CDKN3eRNA expression on CDKN3SE. A knock-down of USF1/USF2 expression in the HepG2 cells did not affect CDKN3eRNA expression, while the expression of CDKN3 was down-regulated. In a USF2 dominant negative HepG2 cell line generated by genome editing, a drastically altered cell shape and lowered cell proliferation rate were found; however, the expression of CDKN3eRNA appeared unaffected. Thus, the present study illustrated two regulators for CDKN3 expression: USF2, as a cell cycle-associated protein regulator, and CDKN3eRNA, as a cell cycle-unassociated RNA regulator.

A comprehensive prognostic and immune analysis of enhancer RNA identifies IGFBP7-AS1 as a novel prognostic biomarker in Uterine Corpus Endometrial Carcinoma.

BACKGROUND: Long non-coding RNAs (lncRNA) have been implicated in a hand of studies that supported an involvement and co-operation in Uterine Corpus Endometrial Carcinoma (UCEC). Enhancer RNAs (eRNA), a functional subtype of lncRNA, have a key role throughout the genome to guide protein production, thus potentially associated with diseases. METHODS: In this study, we mainly applied the Cancer Genome Atlas (TCGA) dataset to systematically discover crucial eRNAs involving UCEC. For the key eRNAs in UCEC, we employed RT-qPCR to compare eRNA expression levels in tumor tissues and paired normal adjacent tissues from UCEC patients for validation. Furthermore, the relationships between the key eRNAs and immune activities were measured from several aspects, including the analysis for tumor microenvironment, immune infiltration cells, immune check point genes, tumor mutation burden, and microsatellite instability, as well as m6A related genes. Finally, the key eRNAs were verified by a comprehensive pan-cancer analysis. RESULTS: IGFBP7 Antisense RNA 1 (IGFBP7-AS1) was identified as the key eRNA for its expression patterns of low levels in tumor tissues and favorable prognostic value in UCEC correlated with its target gene IGFBP7. In RT-qPCR analysis, IGFBP7-AS1 and IGFBP7 had down-regulated expression in tumor tissues, which was consistent with previous analysis. Moreover, IGFBP7-AS1 was found closely related with immune response in relevant immune analyses. Besides, IGFBP7-AS1 and its target gene IGFBP7 correlated with a multi-omics pan-cancer analysis. CONCLUSIONS: Finally, we suggested that IGFBP7-AS1 played a key role in impacting on clinical outcomes of UCEC patients for its possible influence on immune activity.

Anti-diabetic rosiglitazone remodels the adipocyte transcriptome by redistributing transcription to PPARγ-driven enhancers.

Rosiglitazone (rosi) is a powerful insulin sensitizer, but serious toxicities have curtailed its widespread clinical use. Rosi functions as a high-affinity ligand for peroxisome proliferator-activated receptor γ (PPARγ), the adipocyte-predominant nuclear receptor (NR). The classic model, involving binding of ligand to the NR on DNA, explains positive regulation of gene expression, but ligand-dependent repression is not well understood. We addressed this issue by studying the direct effects of rosi on gene transcription using global run-on sequencing (GRO-seq). Rosi-induced changes in gene body transcription were pronounced after 10 min and correlated with steady-state mRNA levels as well as with transcription at nearby enhancers (enhancer RNAs [eRNAs]). Up-regulated eRNAs occurred almost exclusively at PPARγ-binding sites, to which rosi treatment recruited coactivators, including MED1, p300, and CBP. In contrast, transcriptional repression by rosi involved a loss of coactivators from eRNA sites devoid of PPARγ and enriched for other transcription factors, including AP-1 factors and C/EBPs. Thus, rosi activates and represses transcription by fundamentally different mechanisms that could inform the future development of anti-diabetic drugs.

Enhancer RNA (eRNA) in Human Diseases.

Enhancer RNAs (eRNAs), a class of non-coding RNAs (ncRNAs) transcribed from enhancer regions, serve as a type of critical regulatory element in gene expression. There is increasing evidence demonstrating that the aberrant expression of eRNAs can be broadly detected in various human diseases. Some studies also revealed the potential clinical utility of eRNAs in these diseases. In this review, we summarized the recent studies regarding the pathological mechanisms of eRNAs as well as their potential utility across human diseases, including cancers, neurodegenerative disorders, cardiovascular diseases and metabolic diseases. It could help us to understand how eRNAs are engaged in the processes of diseases and to obtain better insight of eRNAs in diagnosis, prognosis or therapy. The studies we reviewed here indicate the enormous therapeutic potency of eRNAs across human diseases.

Additional evidence of HS5-1 enhancer eRNA PEARL for protocadherin alpha gene regulation.

The regulation of target genes by distal enhancers usually determines the fate and function of cells. Active enhancers in specific regions of chromatin may transcribe bidirectionally to produce long non-coding enhancer RNA (eRNA) to regulate gene expression. We recently found that an antisense enhancer eRNA PEARL (Pcdh eRNA associated with R-loop formation) regulates gene expression of members of the Pcdhα cluster via R-loop formation. To further explore the biological function of eRNA, we performed additional genetic and molecular experiments such as CRISPR (clustered regularly interspaced short palindromic repeats) DNA-fragment editing, RT-PCR, and qPCR. First, we performed expression analyses of the HS5-1 eRNA PEARL and found that it was expressed in a tissue-specific manner. In addition, upon CRISPR DNA-fragment deletion or inversion of the CTCF sites in the HS5-1 enhancer region, the expression of eRNA PEARL was reduced to 2%-10% and the expression of Pcdhα gene cluster was also reduced to 13%-68% of the original levels. Finally, deletion of the bidirectional transcription start site (TSS) of HS5-1 eRNA or inversion of TSS of the eRNA PEARL resulted in approximately 60% or 40% decrease of levels of Pcdhα gene expression. In summary, these data suggested a functional role of the HS5-1 eRNA in gene regulation of the Pcdhα cluster, providing a new direction for future researches on the regulatory mechanisms of clustered Pcdh gene expression in the brain.