Nuclear Fractionation Reveals Thousands of Chromatin-Tethered Noncoding RNAs Adjacent to Active Genes

Role of long noncoding RNAs in diabetes-associated peripheral arterial disease

Diabetes mellitus (DM) is a metabolic disease that heightens the risks of many vascular complications, including peripheral arterial disease (PAD). Various types of cells, including but not limited to endothelial cells (ECs), vascular smooth muscle cells (VSMCs), and macrophages (MΦs), play crucial roles in the pathogenesis of DM-PAD. Long non-coding RNAs (lncRNAs) are epigenetic regulators that play important roles in cellular function, and their dysregulation in DM can contribute to PAD. This review focuses on the developing field of lncRNAs and their emerging roles in linking DM and PAD. We review the studies investigating the role of lncRNAs in crucial cellular processes contributing to DM-PAD, including those in ECs, VSMCs, and MΦ. By examining the intricate molecular landscape governed by lncRNAs in these relevant cell types, we hope to shed light on the roles of lncRNAs in EC dysfunction, inflammatory responses, and vascular remodeling contributing to DM-PAD. Additionally, we provide an overview of the research approach and methodologies, from identifying disease-relevant lncRNAs to characterizing their molecular and cellular functions in the context of DM-PAD. We also discuss the potential of leveraging lncRNAs in the diagnosis and therapeutics for DM-PAD. Collectively, this review provides a summary of lncRNA-regulated cell functions contributing to DM-PAD and highlights the translational potential of leveraging lncRNA biology to tackle this increasingly prevalent and complex disease.

Transcription factors ERα and Sox2 have differing multiphasic DNA- and RNA-binding mechanisms

Many transcription factors (TFs) have been shown to bind RNA, leading to open questions regarding the mechanism(s) of this RNA binding and its role in regulating TF activities. Here, we use biophysical assays to interrogate the k on, k off, and K d for DNA and RNA binding of two model human TFs, ERα and Sox2. Unexpectedly, we found that both proteins exhibit multiphasic nucleic acid-binding kinetics. We propose that Sox2 RNA and DNA multiphasic binding kinetics can be explained by a conventional model for sequential Sox2 monomer association and dissociation. In contrast, ERα nucleic acid binding exhibited biphasic dissociation paired with novel triphasic association behavior, in which two apparent binding transitions are separated by a 10-20 min "lag" phase depending on protein concentration. We considered several conventional models for the observed kinetic behavior, none of which adequately explained all the ERα nucleic acid-binding data. Instead, simulations with a model incorporating sequential ERα monomer association, ERα nucleic acid complex isomerization, and product "feedback" on isomerization rate recapitulated the general kinetic trends for both ERα DNA and RNA binding. Collectively, our findings reveal that Sox2 and ERα bind RNA and DNA with previously unappreciated multiphasic binding kinetics, and that their reaction mechanisms differ with ERα binding nucleic acids via a novel reaction mechanism.

Linker histone H1 regulates homeostasis of heterochromatin-associated cRNAs

Chromatin-associated RNAs (cRNAs) are a poorly characterized fraction of cellular RNAs that co-purify with chromatin. Their full complexity and the mechanisms regulating their packaging and chromatin association remain poorly understood. Here, we address these questions in Drosophila. We find that cRNAs constitute a heterogeneous group of RNA species that is abundant in heterochromatic transcripts. We show that heterochromatic cRNAs interact with the heterogeneous nuclear ribonucleoproteins (hnRNP) hrp36/hrp48 and that depletion of linker histone dH1 impairs this interaction. dH1 depletion induces the accumulation of RNA::DNA hybrids (R-loops) in heterochromatin and, as a consequence, increases retention of heterochromatic cRNAs. These effects correlate with increased RNA polymerase II (RNAPII) occupancy at heterochromatin. Notably, impairing cRNA assembly by depletion of hrp36/hrp48 mimics heterochromatic R-loop accumulation induced by dH1 depletion. We also show that dH1 depletion alters nucleosome organization, increasing accessibility of heterochromatin. Altogether, these perturbations facilitate annealing of cRNAs to the DNA template, enhancing R-loop formation and cRNA retention at heterochromatin.

Porcine cis-acting lnc-CAST positively regulates CXCL8 expression through histone H3K27ac

The chemokine CXCL8, also known as the neutrophil chemotactic factor, plays a crucial role in mediating inflammatory responses and managing cellular immune reactions during viral infections. Porcine reproductive and respiratory syndrome virus (PRRSV) primarily infects pulmonary alveolar macrophages (PAMs), leading to acute pulmonary infections. In this study, we explored a novel long non-coding RNA (lncRNA), termed lnc-CAST, situated within the Cxcl8 gene locus. This lncRNA was found to be highly expressed in porcine macrophages. We observed that both lnc-CAST and CXCL8 were significantly upregulated in PAMs following PRRSV infection, and after treatments with lipopolysaccharide (LPS) or lipoteichoic acid (LTA). Furthermore, we noticed a concurrent upregulation of lnc-CAST and CXCL8 expression in lungs of PRRSV-infected pigs. We then determined that lnc-CAST positively influenced CXCL8 expression in PAMs. Overexpression of lnc-CAST led to an increase in CXCL8 production, which in turn enhanced the migration of epithelial cells and the recruitment of neutrophils. Conversely, inhibiting lnc-CAST expression resulted in reduced CXCL8 production in PAMs, leading to decreased migration levels of epithelial cells and neutrophils. From a mechanistic perspective, we found that lnc-CAST, localized in the nucleus, facilitated the enrichment of histone H3K27ac in CXCL8 promoter region, thereby stimulating CXCL8 transcription in a cis-regulatory manner. In conclusion, our study underscores the pivotal critical role of lnc-CAST in regulating CXCL8 production, offering valuable insights into chemokine regulation and lung damage during PRRSV infection.

Transcription factors ERα and Sox2 have differing multiphasic DNA and RNA binding mechanisms

Many transcription factors (TFs) have been shown to bind RNA, leading to open questions regarding the mechanism(s) of this RNA binding and its role in regulating TF activities. Here we use biophysical assays to interrogate the k o n , k o f f , and K d for DNA and RNA binding of two model human transcription factors, ERα and Sox2. Unexpectedly, we found that both proteins exhibited multiphasic nucleic acid binding kinetics. We propose that Sox2 RNA and DNA multiphasic binding kinetics could be explained by a conventional model for sequential Sox2 monomer association and dissociation. In contrast, ERα nucleic acid binding exhibited biphasic dissociation paired with novel triphasic association behavior, where two apparent binding transitions are separated by a 10-20 min "lag" phase depending on protein concentration. We considered several conventional models for the observed kinetic behavior, none of which adequately explained all the ERα nucleic acid binding data. Instead, simulations with a model incorporating sequential ERα monomer association, ERα nucleic acid complex isomerization, and product "feedback" on isomerization rate recapitulated the general kinetic trends for both ERα DNA and RNA binding. Collectively, our findings reveal that Sox2 and ERα bind RNA and DNA with previously unappreciated multiphasic binding kinetics, and that their reaction mechanisms differ with ERα binding nucleic acids via a novel reaction mechanism.

Genomic profiling of six human somatic histone H1 variants denotes that H1X accumulates at recently incorporated transposable elements

Histone H1, a vital component in chromatin structure, binds to linker DNA and regulates nuclear processes. We have investigated the distribution of histone H1 variants in a breast cancer cell line using ChIP-Seq. Two major groups of variants are identified: H1.2, H1.3, H1.5 and H1.0 are abundant in low GC regions (B compartment), while H1.4 and H1X preferentially localize in high GC regions (A compartment). Examining their abundance within transposable elements (TEs) reveals that H1X and H1.4 are enriched in recently-incorporated TEs (SVA and SINE-Alu), while H1.0/H1.2/H1.3/H1.5 are more abundant in older elements. Notably, H1X is particularly enriched in SVA families, while H1.4 shows the highest abundance in young AluY elements. Although low GC variants are generally enriched in LINE, LTR and DNA repeats, H1X and H1.4 are also abundant in a subset of recent LINE-L1 and LTR repeats. H1X enrichment at SVA and Alu is consistent across multiple cell lines. Further, H1X depletion leads to TE derepression, suggesting its role in maintaining TE repression. Overall, this study provides novel insights into the differential distribution of histone H1 variants among repetitive elements, highlighting the potential involvement of H1X in repressing TEs recently incorporated within the human genome.

Noncoding snoRNA host genes are a distinct subclass of long noncoding RNAs

Mammalian genomes are pervasively transcribed into different noncoding (nc)RNA classes, each one with its own hallmarks and exceptions. Some of them are nested into each other, such as host genes for small nucleolar RNAs (snoRNAs), which were long believed to simply act as molecular containers strictly facilitating snoRNA biogenesis. However, recent findings show that noncoding snoRNA host genes (ncSNHGs) display features different from those of 'regular' long ncRNAs (lncRNAs) and, more importantly, they can exert independent and unrelated functions to those of the encoded snoRNAs. Here, we review and summarize past and recent evidence that ncSNHGs form a defined subclass among the plethora of lncRNAs, and discuss future research that can further elucidate their biological relevance.

Co-expression analysis suggests lncRNA-mRNA interactions enhance antiviral immune response during acute Chikungunya fever in whole blood of pediatric patients

Chikungunya virus is an arbovirus that causes the neglected tropical disease chikungunya fever, common in tropical areas worldwide. There is evidence that arboviruses alter host transcriptome and modulate immune response; this modulation may involve transcriptional and post-transcriptional control mechanisms mediated by long non-coding RNAs (lncRNAs). Herein, we employed bioinformatic analysis to evaluate co-expression of lncRNAs and their putative target mRNAs in whole blood during natural Chikungunya infection in adolescent boys. Sequencing data from GSE99992 was uploaded to the Galaxy web server, where data was aligned with HISAT2, gene counts were estimated with HTSeq-count, and differential expression was run with DESeq2. After gene classification with Biomart, Pearson's correlation was applied to identify potential interactions between lncRNAs and mRNAs, which were later classified into cis and trans according to genomic location (FEELnc) and binding potential (LncTar), respectively. We identified 1,975 mRNAs and 793 lncRNAs that were differentially expressed between the acute and convalescent stages of infection in the blood. Of the co-expressed lncRNAs and mRNAs, 357 potentially interact in trans and 9 in cis; their target mRNAs enriched pathways related to immune response and viral infections. Out of 52 enriched KEGG pathways, the RIG-I like receptor signaling is enriched by the highest number of target mRNAs. This pathway starts with the recognition of viral pathogens, leading to innate immune response mediated by the production of IFN-I and inflammatory cytokines. Our findings indicate that alterations in lncRNA expression in adolescent boys, induced by acute Chikungunya infection, potentially modulate mRNAs that contribute to antiviral immune responses.

Hypoxic regulation of hypoxia inducible factor 1 alpha via antisense transcription

Impaired oxygen homeostasis is a frequently encountered pathophysiological factor in multiple complex diseases, including cardiovascular disease and cancer. While the canonical hypoxia response pathway is well characterized, less is known about the role of noncoding RNAs in this process. Here, we investigated the nascent and steady-state noncoding transcriptional responses in endothelial cells and their potential roles in regulating the hypoxic response. Notably, we identify a novel antisense long noncoding RNA that convergently overlaps the majority of the hypoxia inducible factor 1 alpha (HIF1A) locus, which is expressed across several cell types and elevated in atherosclerotic lesions. The antisense (HIF1A-AS) is produced as a stable, unspliced, and polyadenylated nuclear retained transcript. HIF1A-AS is highly induced in hypoxia by both HIF1A and HIF2A and exhibits anticorrelation with the coding HIF1A transcript and protein expression. We further characterized this functional relationship by CRISPR-mediated bimodal perturbation of the HIF1A-AS promoter. We provide evidence that HIF1A-AS represses the expression of HIF1a in cis by repressing transcriptional elongation and deposition of H3K4me3, and that this mechanism is dependent on the act of antisense transcription itself. Overall, our results indicate a critical regulatory role of antisense mediated transcription in regulation of HIF1A expression and cellular response to hypoxia.

The chromatin - triple helix connection

Mammalian genomes are extensively transcribed, producing a large number of coding and non-coding transcripts. A large fraction of the nuclear RNAs is physically associated with chromatin, functioning in gene activation and silencing, shaping higher-order genome organisation, such as involvement in long-range enhancer-promoter interactions, transcription hubs, heterochromatin, nuclear bodies and phase transitions. Different mechanisms allow the tethering of these chromatin-associated RNAs (caRNA) to chromosomes, including RNA binding proteins, the RNA polymerases and R-loops. In this review, we focus on the sequence-specific targeting of RNA to DNA by forming triple helical structures and describe its interplay with chromatin. It turns out that nucleosome positioning at triple helix target sites and the nucleosome itself are essential factors in determining the formation and stability of triple helices. The histone H3-tail plays a critical role in triple helix stabilisation, and the role of its epigenetic modifications in this process is discussed.

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.

Merging short and stranded long reads improves transcript assembly

Long-read RNA sequencing has arisen as a counterpart to short-read sequencing, with the potential to capture full-length isoforms, albeit at the cost of lower depth. Yet this potential is not fully realized due to inherent limitations of current long-read assembly methods and underdeveloped approaches to integrate short-read data. Here, we critically compare the existing methods and develop a new integrative approach to characterize a particularly challenging pool of low-abundance long noncoding RNA (lncRNA) transcripts from short- and long-read sequencing in two distinct cell lines. Our analysis reveals severe limitations in each of the sequencing platforms. For short-read assemblies, coverage declines at transcript termini resulting in ambiguous ends, and uneven low coverage results in segmentation of a single transcript into multiple transcripts. Conversely, long-read sequencing libraries lack depth and strand-of-origin information in cDNA-based methods, culminating in erroneous assembly and quantitation of transcripts. We also discover a cDNA synthesis artifact in long-read datasets that markedly impacts the identity and quantitation of assembled transcripts. Towards remediating these problems, we develop a computational pipeline to "strand" long-read cDNA libraries that rectifies inaccurate mapping and assembly of long-read transcripts. Leveraging the strengths of each platform and our computational stranding, we also present and benchmark a hybrid assembly approach that drastically increases the sensitivity and accuracy of full-length transcript assembly on the correct strand and improves detection of biological features of the transcriptome. When applied to a challenging set of under-annotated and cell-type variable lncRNA, our method resolves the segmentation problem of short-read sequencing and the depth problem of long-read sequencing, resulting in the assembly of coherent transcripts with precise 5' and 3' ends. Our workflow can be applied to existing datasets for superior demarcation of transcript ends and refined isoform structure, which can enable better differential gene expression analyses and molecular manipulations of transcripts.

Global mapping of RNA-chromatin contacts reveals a proximity-dominated connectivity model for ncRNA-gene interactions

Non-coding RNAs (ncRNAs) are transcribed throughout the genome and provide regulatory inputs to gene expression through their interaction with chromatin. Yet, the genomic targets and functions of most ncRNAs are unknown. Here we use chromatin-associated RNA sequencing (ChAR-seq) to map the global network of ncRNA interactions with chromatin in human embryonic stem cells and the dynamic changes in interactions during differentiation into definitive endoderm. We uncover general principles governing the organization of the RNA-chromatin interactome, demonstrating that nearly all ncRNAs exclusively interact with genes in close three-dimensional proximity to their locus and provide a model predicting the interactome. We uncover RNAs that interact with many loci across the genome and unveil thousands of unannotated RNAs that dynamically interact with chromatin. By relating the dynamics of the interactome to changes in gene expression, we demonstrate that activation or repression of individual genes is unlikely to be controlled by a single ncRNA.

Transcriptome and methylome sequencing reveals altered long non-coding RNA genes expression and their aberrant DNA methylation in equine sarcoids

Recent publications confirmed that long non-coding RNAs (lncRNAs) perform an essential function in gene-specific transcription regulation. Nevertheless, despite its important role, lncRNA has not yet been described in equine sarcoids, the skin neoplasia of horses. Therefore, the aim of this study is to deepen the knowledge about lncRNA expression in the pathogenesis of equine sarcoids and provide new insight into the regulatory function of lncRNA in the bovine papillomavirus-dependent neoplasia of horse dermal tissues. RNA sequencing (RNA-seq) data from 12 equine sarcoid samples and the corresponding controls were reanalyzed in this study. A total of 3396 differentially expressed (DE) lncRNAs and 128 DElncRNA-DE genes (DEGs) pairs were identified. Differentially expressed lncRNAs predicted target genes were enriched in pathways associated with inter alia the extracellular matrix disassembly and cancer pathways. Furthermore, methylation data from the same samples were integrated into the analysis, and 12 DElncRNAs were described as potentially disturbed by aberrant methylation. In conclusion, this study presents novel data about lncRNA's role in the pathogenesis of equine sarcoids.

A model for cis-regulation of transcriptional condensates and gene expression by proximal lncRNAs

Long noncoding RNAs (lncRNAs) perform several important functions in cells including cis-regulation of transcription. Barring a few specific cases, the mechanisms underlying transcriptional regulation by lncRNAs remain poorly understood. Transcriptional proteins can form condensates via phase separation at protein-binding loci (BL) on the genome (e.g., enhancers and promoters). lncRNA-coding genes are present at loci in close genomic proximity of these BL and these RNAs can interact with transcriptional proteins via attractive heterotypic interactions mediated by their net charge. Motivated by these observations, we propose that lncRNAs can dynamically regulate transcription in cis via charge-based heterotypic interactions with transcriptional proteins in condensates. To study the consequences of this mechanism, we developed and studied a dynamical phase-field model. We find that proximal lncRNAs can promote condensate formation at the BL. Vicinally localized lncRNA can migrate to the BL to attract more protein because of favorable interaction free energies. However, increasing the distance beyond a threshold leads to a sharp decrease in protein recruitment to the BL. This finding could potentially explain why genomic distances between lncRNA-coding genes and protein-coding genes are conserved across metazoans. Finally, our model predicts that lncRNA transcription can fine-tune transcription from neighboring condensate-controlled genes, repressing transcription from highly expressed genes and enhancing transcription of genes expressed at a low level. This nonequilibrium effect can reconcile conflicting reports that lncRNAs can enhance or repress transcription from proximal genes.

Palmitic Acid-Induced Long Noncoding RNA PARAIL Regulates Inflammation via Interaction With RNA-Binding Protein ELAVL1 in Monocytes and Macrophages

BACKGROUND

Obesity and diabetes are associated with elevated free fatty acids like palmitic acid (PA), which promote chronic inflammation and impaired inflammation resolution associated with cardiometabolic disorders. Long noncoding RNAs (lncRNAs) are implicated in inflammatory processes; however, their roles in PA-regulated inflammation and resolution are unclear.

METHODS

We performed RNA-sequencing analysis to identify PA-regulated coding genes and novel lncRNAs in CD14+ monocytes from healthy volunteers. We investigated the regulation and function of an uncharacterized PA-induced lncRNA PARAIL (PA-regulated anti-inflammatory lncRNA). We examined its role in inflammation resolution by employing knockdown and overexpression strategies in human and mouse macrophages. We also used RNA pulldown coupled with mass spectrometry to identify PARAIL interacting nuclear proteins and their mechanistic involvement in PARAIL functions in human macrophages.

RESULTS

Treatment of human CD14+ monocytes with PA-induced several lncRNAs and genes associated with inflammatory phenotype. PA strongly induced lncRNA PARAIL expressed near RIPK2. PARAIL was also induced by cytokines and infectious agents in human monocytes/macrophages and was regulated by NF-κB (nuclear factor-kappa B). Time course studies showed PARAIL was induced during inflammation resolution phase in PA-treated macrophages. PARAIL knockdown with antisense oligonucleotides upregulated key inflammatory genes and vice versa with PARAIL overexpression. We found that PARAIL interacts with ELAVL1 (ELAV-like RNA-binding protein 1) protein via adenylate/uridylate-rich elements (AU-rich elements; AREs). ELAVL1 knockdown inhibited the anti-inflammatory functions of PARAIL. Moreover, PARAIL knockdown increased cytosolic localization of ELAVL1 and increased the stability of ARE-containing inflammatory genes. Mouse orthologous Parail was downregulated in macrophages from mice with diabetes and atherosclerosis. Parail overexpression attenuated proinflammatory genes in mouse macrophages.

CONCLUSIONS

Upregulation of PARAIL under acute inflammatory conditions contributes to proresolution mechanisms via PARAIL-ELAVL1 interactions. Conversely, PARAIL downregulation in cardiometabolic diseases enhances ELAVL1 function and impairs inflammation resolution to further augment inflammation. Thus, inflammation-resolving lncRNAs like PARAIL represent novel targets to combat inflammatory cardiometabolic diseases.

Profiling lariat intermediates reveals genetic determinants of early and late co-transcriptional splicing

Long introns with short exons in vertebrate genes are thought to require spliceosome assembly across exons (exon definition), rather than introns, thereby requiring transcription of an exon to splice an upstream intron. Here, we developed CoLa-seq (co-transcriptional lariat sequencing) to investigate the timing and determinants of co-transcriptional splicing genome wide. Unexpectedly, 90% of all introns, including long introns, can splice before transcription of a downstream exon, indicating that exon definition is not obligatory for most human introns. Still, splicing timing varies dramatically across introns, and various genetic elements determine this variation. Strong U2AF2 binding to the polypyrimidine tract predicts early splicing, explaining exon definition-independent splicing. Together, our findings question the essentiality of exon definition and reveal features beyond intron and exon length that are determinative for splicing timing.

An Extended DNA Binding Domain of the Estrogen Receptor Alpha Directly Interacts with RNAs in Vitro

Estrogen receptor alpha (ERα) is a ligand-responsive transcription factor critical for sex determination and development. Recent reports challenge the canonical view of ERα function by suggesting an activity beyond binding dsDNA at estrogen-responsive promotor elements: association with RNAs in vivo. Whether these interactions are direct or indirect remains unknown, which limits the ability to understand the extent, specificity, and biological role of ERα-RNA binding. Here we demonstrate that an extended DNA-binding domain of ERα directly binds a wide range of RNAs in vitro with structural specificity. ERα binds RNAs that adopt a range of hairpin-derived structures independent of sequence, while interacting poorly with single- and double-stranded RNA. RNA affinities are only 4-fold weaker than consensus dsDNA and significantly tighter than nonconsensus dsDNA sequences. Moreover, RNA binding is competitive with DNA binding. Together, these data show that ERα utilizes an extended DNA-binding domain to achieve a high-affinity/low-specificity mode for interacting with RNA.

HSV-1 and influenza infection induce linear and circular splicing of the long NEAT1 isoform

The herpes simplex virus 1 (HSV-1) virion host shut-off (vhs) protein cleaves both cellular and viral mRNAs by a translation-initiation-dependent mechanism, which should spare circular RNAs (circRNAs). Here, we show that vhs-mediated degradation of linear mRNAs leads to an enrichment of circRNAs relative to linear mRNAs during HSV-1 infection. This was also observed in influenza A virus (IAV) infection, likely due to degradation of linear host mRNAs mediated by the IAV PA-X protein and cap-snatching RNA-dependent RNA polymerase. For most circRNAs, enrichment was not due to increased circRNA synthesis but due to a general loss of linear RNAs. In contrast, biogenesis of a circRNA originating from the long isoform (NEAT1_2) of the nuclear paraspeckle assembly transcript 1 (NEAT1) was induced both in HSV-1 infection-in a vhs-independent manner-and in IAV infection. This was associated with induction of novel linear splicing of NEAT1_2 both within and downstream of the circRNA. NEAT1_2 forms a scaffold for paraspeckles, nuclear bodies located in the interchromatin space, must likely remain unspliced for paraspeckle assembly and is up-regulated in HSV-1 and IAV infection. We show that NEAT1_2 splicing and up-regulation can be induced by ectopic co-expression of the HSV-1 immediate-early proteins ICP22 and ICP27, potentially linking increased expression and splicing of NEAT1_2. To identify other conditions with NEAT1_2 splicing, we performed a large-scale screen of published RNA-seq data. This uncovered both induction of NEAT1_2 splicing and poly(A) read-through similar to HSV-1 and IAV infection in cancer cells upon inhibition or knockdown of CDK7 or the MED1 subunit of the Mediator complex phosphorylated by CDK7. In summary, our study reveals induction of novel circular and linear NEAT1_2 splicing isoforms as a common characteristic of HSV-1 and IAV infection and highlights a potential role of CDK7 in HSV-1 or IAV infection.

LncRNA RUS shapes the gene expression program towards neurogenesis

The chromatin-associated lncRNA RUS binds in the vicinity to neural differentiation-associated genes and regulates them in a context-dependent manner to enable proper neuron development.

The evolution of brain complexity correlates with an increased expression of long, noncoding (lnc) RNAs in neural tissues. Although prominent examples illustrate the potential of lncRNAs to scaffold and target epigenetic regulators to chromatin loci, only few cases have been described to function during brain development. We present a first functional characterization of the lncRNA LINC01322, which we term RUS for “RNA upstream of Slitrk3.” The RUS gene is well conserved in mammals by sequence and synteny next to the neurodevelopmental gene Slitrk3. RUS is exclusively expressed in neural cells and its expression increases during neuronal differentiation of mouse embryonic cortical neural stem cells. Depletion of RUS locks neuronal precursors in an intermediate state towards neuronal differentiation resulting in arrested cell cycle and increased apoptosis. RUS associates with chromatin in the vicinity of genes involved in neurogenesis, most of which change their expression upon RUS depletion. The identification of a range of epigenetic regulators as specific RUS interactors suggests that the lncRNA may mediate gene activation and repression in a highly context-dependent manner.

Histone H1 regulates non-coding RNA turnover on chromatin in a m6A-dependent manner

Linker histones are highly abundant chromatin-associated proteins with well-established structural roles in chromatin and as general transcriptional repressors. In addition, it has been long proposed that histone H1 exerts context-specific effects on gene expression. Here, we identify a function of histone H1 in chromatin structure and transcription using a range of genomic approaches. In the absence of histone H1, there is an increase in the transcription of non-coding RNAs, together with reduced levels of m6A modification leading to their accumulation on chromatin and causing replication-transcription conflicts. This strongly suggests that histone H1 prevents non-coding RNA transcription and regulates non-coding transcript turnover on chromatin. Accordingly, altering the m6A RNA methylation pathway rescues the replicative phenotype of H1 loss. This work unveils unexpected regulatory roles of histone H1 on non-coding RNA turnover and m6A deposition, highlighting the intimate relationship between chromatin conformation, RNA metabolism, and DNA replication to maintain genome performance.

Non-coding RNAs associated with Prader-Willi syndrome regulate transcription of neurodevelopmental genes in human induced pluripotent stem cells

Mutations and aberrant gene expression during cellular differentiation lead to neurodevelopmental disorders, such as Prader-Willi syndrome (PWS), which results from the deletion of an imprinted locus on paternally inherited chromosome 15. We analyzed chromatin-associated RNA in human induced pluripotent cells (iPSCs) upon depletion of hybrid small nucleolar long non-coding RNAs (sno-lncRNAs) and 5' snoRNA capped and polyadenylated long non-coding RNAs (SPA-lncRNAs) transcribed from the locus deleted in PWS. We found that rapid ablation of these lncRNAs affects transcription of specific gene classes. Downregulated genes contribute to neurodevelopment and neuronal maintenance, while upregulated genes are predominantly involved in the negative regulation of cellular metabolism and apoptotic processes. Our data reveal the importance of SPA-lncRNAs and sno-lncRNAs in controlling gene expression in iPSCs and provide a platform for synthetic experimental approaches in PWS studies. We conclude that ncRNAs transcribed from the PWS locus are critical regulators of a transcriptional signature, which is important for neuronal differentiation and development.

Chromatin Hubs: A biological and computational outlook

This review discusses our current understanding of chromatin biology and bioinformatics under the unifying concept of “chromatin hubs.” The first part reviews the biology of chromatin hubs, including chromatin–chromatin interaction hubs, chromatin hubs at the nuclear periphery, hubs around macromolecules such as RNA polymerase or lncRNAs, and hubs around nuclear bodies such as the nucleolus or nuclear speckles. The second part reviews existing computational methods, including enhancer–promoter interaction prediction, network analysis, chromatin domain callers, transcription factory predictors, and multi-way interaction analysis. We introduce an integrated model that makes sense of the existing evidence. Understanding chromatin hubs may allow us (i) to explain long-unsolved biological questions such as interaction specificity and redundancy of mechanisms, (ii) to develop more realistic kinetic and functional predictions, and (iii) to explain the etiology of genomic disease.

Emerging Functions of lncRNA Loci beyond the Transcript Itself

Thousands of long noncoding RNAs (lncRNAs) are actively transcribed in mammalian genomes. This class of RNAs has important regulatory functions in a broad range of cellular processes and diseases. Numerous lncRNAs have been demonstrated to mediate gene regulation through RNA-based mechanisms. Simultaneously, non-functional lncRNA transcripts derived from the activity of lncRNA loci have been identified, which underpin the notion that a considerable fraction of lncRNA loci exert regulatory functions through mechanisms associated with the production or the activity of lncRNA loci beyond the synthesized transcripts. We particularly distinguish two main RNA-independent components associated with regulatory effects; the act of transcription and the activity of DNA regulatory elements. We describe the experimental approaches to distinguish and understand the functional mechanisms derived from lncRNA loci. These scenarios reveal emerging mechanisms important to understanding the lncRNA implications in genome biology.

Long non-coding RNAs and exosomal lncRNAs: Potential functions in lung cancer progression, drug resistance and tumor microenvironment remodeling

Among the different kinds of tumors threatening human life, lung cancer is one that is commonly observed in both males and females. The aggressive behavior of lung cancer and interactions occurring in tumor microenvironment enhances the malignancy of this tumor. The lung tumor cells have demonstrated capacity in developing chemo- and radio-resistance. LncRNAs are a category of non-coding RNAs that do not encode proteins, but their aberrant expression is responsible for tumor development, especially lung cancer. In the present review, we focus on both lncRNAs and exosomal lncRNAs in lung cancer, and their ability in regulating proliferation and metastasis. Cell cycle progression and molecular mechanisms related to lung cancer metastasis such as EMT and MMPs are regulated by lncRNAs. LncRNAs interact with miRNAs, STAT, Wnt, EZH2, PTEN and PI3K/Akt signaling pathways to affect progression of lung cancer cells. LncRNAs demonstrate both tumor-suppressor and tumor-promoting functions in lung cancer. They can be considered as biomarkers in lung cancer and especially exosomal lncRNAs present in body fluids are potential tools for minimally invasive diagnosis. Furthermore, we discuss regulation of lncRNAs by anti-cancer drugs and genetic tools as well as the role of these factors in therapy response of lung cancer cells.

Ribosomal RNA regulates chromosome clustering during mitosis

Noncoding RNAs are known to associate with mitotic chromosomes, but the identities and functions of chromosome-associated RNAs in mitosis remain elusive. Here, we show that rRNA species associate with condensed chromosomes during mitosis. In particular, pre-rRNAs such as 45S, 32S, and 30S are highly enriched on mitotic chromosomes. Immediately following nucleolus disassembly in mitotic prophase, rRNAs are released and associate with and coat each condensed chromosome at prometaphase. Using unbiased mass spectrometry analysis, we further demonstrate that chromosome-bound rRNAs are associated with Ki-67. Moreover, the FHA domain and the repeat region of Ki-67 recognize and anchor rRNAs to chromosomes. Finally, suppression of chromosome-bound rRNAs by RNA polymerase I inhibition or by using rRNA-binding-deficient Ki-67 mutants impair mitotic chromosome dispersion during prometaphase. Our study thus reveals an important role of rRNAs in preventing chromosome clustering during mitosis.

Long intergenic non-protein coding RNA 847 promotes laryngeal squamous cell carcinoma progression through the microRNA-181a-5p/zinc finger E-box binding homeobox 2 axis

The present study is targeted at investigating the effects of long intergenic non-protein coding RNA 847 (LINC00847) on the malignant biological behaviors of laryngeal squamous cell carcinoma (LSCC) cells, and the mechanisms. Quantitative real-time PCR and Western blotting were conducted for detecting the expressions of LINC00847, microRNA-181a-5p (miR-181a-5p) and zinc finger E-box binding homeobox 2 (ZEB2) in LSCC cell lines and tissue samples. BrdU, cell counting kit-8, scratch wound healing, Transwell and flow cytometry assays were utilized for detecting cell proliferation, migration, invasion, and cell cycle progression. Dual-luciferase reporter gene, RNA binding protein immunoprecipitation (RIP), and RNA pull-down assays were utilized to investigate the interaction among LINC00847, miR-181a-5p, and ZEB2. The subcellular location of LINC00847 was determined by RNA fluorescence in situ hybridization (RNA-FISH) assay. Tumor growth was evaluated using a xenograft model of nude mice. It was revealed that LINC00847 expression was increased in LSCC tissues, and its high expression was associated with lymph node metastasis and poor differentiation. LINC00847 was mainly located in the cytoplasm of LSCC cells, and LINC00847 overexpression promoted LSCC cell proliferation, migration, invasion, and accelerated the cell cycle progression while knocking down LINC00847 had the opposite effects in vitro and inhibited the tumor growth in vivo. LINC00847 directly targeted miR-181a-5p and negatively modulated miR-181a-5p expression. ZEB2 was a target gene of miR-181a-5p, and was positively and indirectly modulated by LINC00847. Our data suggest that LINC00847 promotes LSCC progression by regulating the miR-181a-5p/ZEB2 axis.

Deciphering the RNA universe in sperm in its role as a vertical information carrier

The inheritance of neurophysiologic and neuropsychologic complex diseases can only partly be explained by the Mendelian concept of genetic inheritance. Previous research showed that both psychological disorders like post-traumatic stress disorder and metabolic diseases are more prevalent in the progeny of affected parents. This could suggest an epigenetic mode of transmission. Human studies give first insight into the scope of intergenerational influence of stressors but are limited in exploring the underlying mechanisms. Animal models have elucidated the mechanistic underpinnings of epigenetic transmission. In this review, we summarize progress on the mechanisms of paternal intergenerational transmission by means of sperm RNA in mouse models. We discuss relevant details for the modelling of RNA-mediated transmission, point towards currently unanswered questions and propose experimental considerations for tackling these questions.

Non-coding RNA LEVER sequestration of PRC2 can mediate long range gene regulation

Polycomb Repressive Complex 2 (PRC2) is an epigenetic regulator required for gene silencing during development. Although PRC2 is a well-established RNA-binding complex, the biological function of PRC2-RNA interaction has been controversial. Here, we study the gene-regulatory role of the inhibitory PRC2-RNA interactions. We report a nuclear long non-coding RNA, LEVER, which mapped 236 kb upstream of the β-globin cluster as confirmed by Nanopore sequencing. LEVER RNA interacts with PRC2 in its nascent form, and this prevents the accumulation of the H3K27 repressive histone marks within LEVER locus. Interestingly, the accessible LEVER chromatin, in turn, suppresses the chromatin interactions between the ε-globin locus and β-globin locus control region (LCR), resulting in a repressive effect on ε-globin gene expression. Our findings validate that the nascent RNA-PRC2 interaction inhibits local PRC2 function in situ. More importantly, we demonstrate that such a local process can in turn regulate the expression of neighboring genes.

Identification of a long non-coding RNA LEVER, that inhibits the Polycomb Repressive Complex 2 (PRC2) and controls nearby embryonic form of beta-globin gene, provides additional evidence for PRC2-RNA functional interaction.

Chromatin-associated orphan snoRNA regulates DNA damage-mediated differentiation via a non-canonical complex

Small nucleolar RNAs (snoRNAs) are commonly acknowledged as a class of homogeneous non-coding RNAs that guide ribosomal RNA modifications. However, snoRNAs referred to as orphans have largely unknown functions. Here, we systematically profile chromatin-associated snoRNAs (casnoRNAs) in mammalian cells and identify a subgroup of orphan casnoRNAs responding to DNA damage stress, among which SNORA73 shows the most marked reduction in chromatin enrichment. Downregulated SNORA73 maintains cancer genome stability and differentiation block in hematopoietic malignancy. Mechanistically, casnoRNA the 5' end non-canonical structure of SNORA73 is critical for its function and binding to poly (ADP-ribose) polymerase 1 (PARP1). SNORA73 inhibits PARP1 auto-PARylation to affect cancer genome stability by forming a small nucleolar ribonucleoprotein (snoRNP) with PARP1 and canonical H/ACA proteins DKC1/NHP2. Our findings reveal the role of an orphan snoRNA serving as casnoRNA and highlights a link between non-canonical structure of snoRNA and their functional diversity.

Subcellular Transcriptomics and Proteomics: A Comparative Methods Review

The internal environment of cells is molecularly crowded, which requires spatial organization via subcellular compartmentalization. These compartments harbor specific conditions for molecules to perform their biological functions, such as coordination of the cell cycle, cell survival, and growth. This compartmentalization is also not static, with molecules trafficking between these subcellular neighborhoods to carry out their functions. For example, some biomolecules are multifunctional, requiring an environment with differing conditions or interacting partners, and others traffic to export such molecules. Aberrant localization of proteins or RNA species has been linked to many pathological conditions, such as neurological, cancer, and pulmonary diseases. Differential expression studies in transcriptomics and proteomics are relatively common, but the majority have overlooked the importance of subcellular information. In addition, subcellular transcriptomics and proteomics data do not always colocate because of the biochemical processes that occur during and after translation, highlighting the complementary nature of these fields. In this review, we discuss and directly compare the current methods in spatial proteomics and transcriptomics, which include sequencing- and imaging-based strategies, to give the reader an overview of the current tools available. We also discuss current limitations of these strategies as well as future developments in the field of spatial -omics.

Genome-wide analysis and functional annotation of chromatin-enriched noncoding RNAs in rice during somatic cell regeneration

BACKGROUND

Plants have the remarkable ability to generate callus, a pluripotent cell mass that acquires competence for subsequent tissue regeneration. Global chromatin remodeling is required for this cell fate transition, but how the process is regulated is not fully understood. Chromatin-enriched noncoding RNAs (cheRNAs) are thought to play important roles in maintaining chromatin state. However, whether cheRNAs participate in somatic cell regeneration in plants has not yet been clarified.

RESULTS

To uncover the characteristics and functions of cheRNAs during somatic cell reprogramming in plants, we systematically investigate cheRNAs during callus induction, proliferation and regeneration in rice. We identify 2284 cheRNAs, most of which are novel long non-coding RNAs or small nucleolar RNAs. These cheRNAs, which are highly conserved across plant species, shuttle between chromatin and the nucleoplasm during somatic cell regeneration. They positively regulate the expression of neighboring genes via specific RNA motifs, which may interact with DNA motifs around cheRNA loci. Large-scale mutant analysis shows that cheRNAs are associated with plant size and seed morphology. Further detailed functional investigation of two che-lncRNAs demonstrates that their loss of function impairs cell dedifferentiation and plant regeneration, highlighting the functions of cheRNAs in regulating the expression of neighboring genes via specific motifs. These findings support cis- regulatory roles of cheRNAs in influencing a variety of rice traits.

CONCLUSIONS

cheRNAs are a distinct subclass of regulatory non-coding RNAs that are required for somatic cell regeneration and regulate rice traits. Targeting cheRNAs has great potential for crop trait improvement and breeding in future.

The METTL5-TRMT112 N6-methyladenosine methyltransferase complex regulates mRNA translation via 18S rRNA methylation

Ribosomal RNAs (rRNAs) have long been known to carry chemical modifications, including 2'O-methylation, pseudouridylation, N⁶-methyladenosine (m⁶A), and N^6,6-dimethyladenosine. While the functions of many of these modifications are unclear, some are highly conserved and occur in regions of the ribosome critical for mRNA decoding. Both 28S rRNA and 18S rRNA carry single m⁶A sites, and while the methyltransferase ZCCHC4 has been identified as the enzyme responsible for the 28S rRNA m⁶A modification, the methyltransferase responsible for the 18S rRNA m⁶A modification has remained unclear. Here, we show that the METTL5-TRMT112 methyltransferase complex installs the m⁶A modification at position 1832 of human 18S rRNA. Our work supports findings that TRMT112 is required for METTL5 stability and reveals that human METTL5 mutations associated with microcephaly and intellectual disability disrupt this interaction. We show that loss of METTL5 in human cancer cell lines and in mice regulates gene expression at the translational level; additionally, Mettl5 knockout mice display reduced body size and evidence of metabolic defects. While recent work has focused heavily on m⁶A modifications in mRNA and their roles in mRNA processing and translation, we demonstrate here that deorphanizing putative methyltransferase enzymes can reveal previously unappreciated regulatory roles for m⁶A in noncoding RNAs.

The genetic basis of urate control and gout: Insights into molecular pathogenesis from follow-up study of genome-wide association study loci

This review focuses on the post-genome-wide association study (GWAS) era in gout, i.e., the translation of GWAS genetic association signals into biologically informative knowledge. Analytical and experimental follow-up of individual loci, based on the identification of causal genetic variants, reveals molecular pathogenic pathways. We summarize in detail the largest GWAS in urate to date, then we review follow-up studies and molecular insights from ABCG2, HNF4A, PDZK1, MAF, GCKR, ALDH2, ALDH16A1, SLC22A12, SLC2A9, ABCC4, and SLC22A13, including the role of insulin signaling. One common factor in these pathways is the importance of transcriptional control, including the HNF4α transcription factor. The new molecular knowledge reveals new targets for intervention to manage urate levels and prevent gout.

Chromatin lncRNA Platr10 controls stem cell pluripotency by coordinating an intrachromosomal regulatory network

BACKGROUND

A specific 3-dimensional intrachromosomal architecture of core stem cell factor genes is required to reprogram a somatic cell into pluripotency. As little is known about the epigenetic readers that orchestrate this architectural remodeling, we used a novel chromatin RNA in situ reverse transcription sequencing (CRIST-seq) approach to profile long noncoding RNAs (lncRNAs) in the Oct4 promoter.

RESULTS

We identify Platr10 as an Oct4 - Sox2 binding lncRNA that is activated in somatic cell reprogramming. Platr10 is essential for the maintenance of pluripotency, and lack of this lncRNA causes stem cells to exit from pluripotency. In fibroblasts, ectopically expressed Platr10 functions in trans to activate core stem cell factor genes and enhance pluripotent reprogramming. Using RNA reverse transcription-associated trap sequencing (RAT-seq), we show that Platr10 interacts with multiple pluripotency-associated genes, including Oct4, Sox2, Klf4, and c-Myc, which have been extensively used to reprogram somatic cells. Mechanistically, we demonstrate that Platr10 helps orchestrate intrachromosomal promoter-enhancer looping and recruits TET1, the enzyme that actively induces DNA demethylation for the initiation of pluripotency. We further show that Platr10 contains an Oct4 binding element that interacts with the Oct4 promoter and a TET1-binding element that recruits TET1. Mutation of either of these two elements abolishes Platr10 activity.

CONCLUSION

These data suggest that Platr10 functions as a novel chromatin RNA molecule to control pluripotency in trans by modulating chromatin architecture and regulating DNA methylation in the core stem cell factor network.

Nascent RNA scaffolds contribute to chromosome territory architecture and counter chromatin compaction

Nuclear chromosomes transcribe far more RNA than required to encode protein. Here we investigate whether non-coding RNA broadly contributes to cytological-scale chromosome territory architecture. We develop a procedure that depletes soluble proteins, chromatin, and most nuclear RNA from the nucleus but does not delocalize XIST, a known architectural RNA, from an insoluble chromosome "scaffold." RNA-seq analysis reveals that most RNA in the nuclear scaffold is repeat-rich, non-coding, and derived predominantly from introns of nascent transcripts. Insoluble, repeat-rich (C₀T-1) RNA co-distributes with known scaffold proteins including scaffold attachment factor A (SAF-A), and distribution of these components inversely correlates with chromatin compaction in normal and experimentally manipulated nuclei. We further show that RNA is required for SAF-A to interact with chromatin and for enrichment of structurally embedded "scaffold attachment regions" prevalent in euchromatin. Collectively, the results indicate that long nascent transcripts contribute a dynamic structural role that promotes the open architecture of active chromosome territories.

Nascent RNA antagonizes the interaction of a set of regulatory proteins with chromatin

A number of regulatory factors are recruited to chromatin by specialized RNAs. Whether RNA has a more general role in regulating the interaction of proteins with chromatin has not been determined. We used proteomics methods to measure the global impact of nascent RNA on chromatin in embryonic stem cells. Surprisingly, we found that nascent RNA primarily antagonized the interaction of chromatin modifiers and transcriptional regulators with chromatin. Transcriptional inhibition and RNA degradation induced recruitment of a set of transcriptional regulators, chromatin modifiers, nucleosome remodelers, and regulators of higher-order structure. RNA directly bound to factors, including BAF, NuRD, EHMT1, and INO80 and inhibited their interaction with nucleosomes. The transcriptional elongation factor P-TEFb directly bound pre-mRNA, and its recruitment to chromatin upon Pol II inhibition was regulated by the 7SK ribonucleoprotein complex. We postulate that by antagonizing the interaction of regulatory proteins with chromatin, nascent RNA links transcriptional output with chromatin composition.

Brain Long Noncoding RNAs: Multitask Regulators of Neuronal Differentiation and Function

The extraordinary cellular diversity and the complex connections established within different cells types render the nervous system of vertebrates one of the most sophisticated tissues found in living organisms. Such complexity is ensured by numerous regulatory mechanisms that provide tight spatiotemporal control, robustness and reliability. While the unusual abundance of long noncoding RNAs (lncRNAs) in nervous tissues was traditionally puzzling, it is becoming clear that these molecules have genuine regulatory functions in the brain and they are essential for neuronal physiology. The canonical view of RNA as predominantly a 'coding molecule' has been largely surpassed, together with the conception that lncRNAs only represent 'waste material' produced by cells as a side effect of pervasive transcription. Here we review a growing body of evidence showing that lncRNAs play key roles in several regulatory mechanisms of neurons and other brain cells. In particular, neuronal lncRNAs are crucial for orchestrating neurogenesis, for tuning neuronal differentiation and for the exact calibration of neuronal excitability. Moreover, their diversity and the association to neurodegenerative diseases render them particularly interesting as putative biomarkers for brain disease. Overall, we foresee that in the future a more systematic scrutiny of lncRNA functions will be instrumental for an exhaustive understanding of neuronal pathophysiology.

A myocardin-adjacent lncRNA balances SRF-dependent gene transcription in the heart

In this study, Anderson et al. identify a cardiac lncRNA transcribed adjacent to myocardin, named CARDINAL, that antagonizes SRF-dependent mitogenic gene transcription in the heart. They also show that CARDINAL forms a nuclear complex with SRF and inhibits TCF-mediated transactivation of the promitogenic gene c-fos, suggesting CARDINAL functions as an essential RNA cofactor for SRF in the heart.

Myocardin, a potent coactivator of serum response factor (SRF), competes with ternary complex factor (TCF) proteins for SRF binding to balance opposing mitogenic and myogenic gene programs in cardiac and smooth muscle. Here we identify a cardiac lncRNA transcribed adjacent to myocardin, named CARDINAL, which antagonizes SRF-dependent mitogenic gene transcription in the heart. CARDINAL-deficient mice show ectopic TCF/SRF-dependent mitogenic gene expression and decreased cardiac contractility in response to age and ischemic stress. CARDINAL forms a nuclear complex with SRF and inhibits TCF-mediated transactivation of the promitogenic gene c-fos, suggesting CARDINAL functions as an RNA cofactor for SRF in the heart.

New insights into the functional role of retrotransposon dynamics in mammalian somatic cells

Retrotransposons are genetic elements present across all eukaryotic genomes. While their role in evolution is considered as a potentially beneficial natural source of genetic variation, their activity is classically considered detrimental due to their potentially harmful effects on genome stability. However, studies are increasingly shedding light on the regulatory function and beneficial role of somatic retroelement reactivation in non-pathological contexts. Here, we review recent findings unveiling the regulatory potential of retrotransposons, including their role in noncoding RNA transcription, as modulators of mammalian transcriptional and epigenome landscapes. We also discuss technical challenges in deciphering the multifaceted activity of retrotransposable elements, highlighting an unforeseen central role of this neglected portion of the genome both in early development and in adult life.

Intergenic RNA mainly derives from nascent transcripts of known genes

BACKGROUND

Eukaryotic genomes undergo pervasive transcription, leading to the production of many types of stable and unstable RNAs. Transcription is not restricted to regions with annotated gene features but includes almost any genomic context. Currently, the source and function of most RNAs originating from intergenic regions in the human genome remain unclear.

RESULTS

We hypothesize that many intergenic RNAs can be ascribed to the presence of as-yet unannotated genes or the "fuzzy" transcription of known genes that extends beyond the annotated boundaries. To elucidate the contributions of these two sources, we assemble a dataset of more than 2.5 billion publicly available RNA-seq reads across 5 human cell lines and multiple cellular compartments to annotate transcriptional units in the human genome. About 80% of transcripts from unannotated intergenic regions can be attributed to the fuzzy transcription of existing genes; the remaining transcripts originate mainly from putative long non-coding RNA loci that are rarely spliced. We validate the transcriptional activity of these intergenic RNAs using independent measurements, including transcriptional start sites, chromatin signatures, and genomic occupancies of RNA polymerase II in various phosphorylation states. We also analyze the nuclear localization and sensitivities of intergenic transcripts to nucleases to illustrate that they tend to be rapidly degraded either on-chromatin by XRN2 or off-chromatin by the exosome.

CONCLUSIONS

We provide a curated atlas of intergenic RNAs that distinguishes between alternative processing of well-annotated genes from independent transcriptional units based on the combined analysis of chromatin signatures, nuclear RNA localization, and degradation pathways.

Sequence deeper without sequencing more: Bayesian resolution of ambiguously mapped reads

Next-generation sequencing (NGS) has transformed molecular biology and contributed to many seminal insights into genomic regulation and function. Apart from whole-genome sequencing, an NGS workflow involves alignment of the sequencing reads to the genome of study, after which the resulting alignments can be used for downstream analyses. However, alignment is complicated by the repetitive sequences; many reads align to more than one genomic locus, with 15-30% of the genome not being uniquely mappable by short-read NGS. This problem is typically addressed by discarding reads that do not uniquely map to the genome, but this practice can lead to systematic distortion of the data. Previous studies that developed methods for handling ambiguously mapped reads were often of limited applicability or were computationally intensive, hindering their broader usage. In this work, we present SmartMap: an algorithm that augments industry-standard aligners to enable usage of ambiguously mapped reads by assigning weights to each alignment with Bayesian analysis of the read distribution and alignment quality. SmartMap is computationally efficient, utilizing far fewer weighting iterations than previously thought necessary to process alignments and, as such, analyzing more than a billion alignments of NGS reads in approximately one hour on a desktop PC. By applying SmartMap to peak-type NGS data, including MNase-seq, ChIP-seq, and ATAC-seq in three organisms, we can increase read depth by up to 53% and increase the mapped proportion of the genome by up to 18% compared to analyses utilizing only uniquely mapped reads. We further show that SmartMap enables the analysis of more than 140,000 repetitive elements that could not be analyzed by traditional ChIP-seq workflows, and we utilize this method to gain insight into the epigenetic regulation of different classes of repetitive elements. These data emphasize both the dangers of discarding ambiguously mapped reads and their power for driving biological discovery.

Insight Into Chromatin-Enriched RNA: A Key Chromatin Regulator in Tumors

Chromatin-enriched RNAs (cheRNAs) constitute a special class of long noncoding RNAs (lncRNAs) that are enriched around chromatin and function to activate neighboring or distal gene transcription. Recent studies have shown that cheRNAs affect chromatin structure and gene expression by recruiting chromatin modifiers or acting as bridges between distal enhancers and promoters. The abnormal transcription of cheRNAs plays an important role in the occurrence of many diseases, particularly tumors. The critical effect of cancer stem cells (CSCs) on the formation and development of tumors is well known, but the function of cheRNAs in tumorigenesis, especially in CSC proliferation and stemness maintenance, is not yet fully understood. This review focuses on the mechanisms of cheRNAs in epigenetic regulation and chromatin conformation and discusses the way cheRNAs function in CSCs to deepen the understanding of tumorigenesis and provide novel insight to advance tumor-targeting therapy.

RNAs as Regulators of Cellular Matchmaking

RNA molecules are increasingly being identified as facilitating or impeding the interaction of proteins and nucleic acids, serving as so-called scaffolds or decoys. Long non-coding RNAs have been commonly implicated in such roles, particularly in the regulation of nuclear processes including chromosome topology, regulation of chromatin state and gene transcription, and assembly of nuclear biomolecular condensates such as paraspeckles. Recently, an increased awareness of cytoplasmic RNA scaffolds and decoys has begun to emerge, including the identification of non-coding regions of mRNAs that can also function in a scaffold-like manner to regulate interactions of nascently translated proteins. Collectively, cytoplasmic RNA scaffolds and decoys are now implicated in processes such as mRNA translation, decay, protein localization, protein degradation and assembly of cytoplasmic biomolecular condensates such as P-bodies. Here, we review examples of RNA scaffolds and decoys in both the nucleus and cytoplasm, illustrating common themes, the suitability of RNA to such roles, and future challenges in identifying and better understanding RNA scaffolding and decoy functions.

Modulation of Brain Pathology by Enhancer RNAs in Cerebral Ischemia

Recent studies have reported widespread stimulus-dependent transcription of mammalian enhancers into noncoding enhancer RNAs (eRNAs), some of which have central roles in the enhancer-mediated induction of target genes and modulation of phenotypic outcomes during development and disease. In cerebral ischemia, the expression and functions of eRNAs are virtually unknown. Here, we applied genome-wide H3K27ac ChIP-seq and genome-wide RNA-seq to identify enhancer elements and stroke-induced eRNAs, respectively, in the mouse cerebral cortex during transient focal ischemia. Following a 1-h middle cerebral artery occlusion (MCAO) and 6 h of reperfusion, we identified 77 eRNAs that were significantly upregulated in stroke as compared to sham, of which 55 were exclusively expressed in stroke. The knockdown of two stroke-induced eRNAs in the mouse brain resulted in significantly larger infarct volumes as compared to controls, suggesting that these eRNAs are involved in the post-stroke neuroprotective response. A preliminary comparison of eRNA expression in the male versus female cortices revealed sex-dependent patterns that may underlie the physiological differences in response to stroke between the two sexes. Together, this study is the first to illuminate the eRNA landscape in the post-stroke cortex and demonstrate the significance of an eRNA in modulating post-stroke cortical brain damage.

Global analysis of expression, maturation and subcellular localization of mouse liver transcriptome identifies novel sex-biased and TCPOBOP-responsive long non-coding RNAs

BACKGROUND

While nuclear transcription and RNA processing and localization are well established for protein coding genes (PCGs), these processes are poorly understood for long non-coding (lnc)RNAs. Here, we characterize global patterns of transcript expression, maturation and localization for mouse liver RNA, including more than 15,000 lncRNAs. PolyA-selected liver RNA was isolated and sequenced from four subcellular fractions (chromatin, nucleoplasm, total nucleus, and cytoplasm), and from the chromatin-bound fraction without polyA selection.

RESULTS

Transcript processing, determined from normalized intronic to exonic sequence read density ratios, progressively increased for PCG transcripts in going from the chromatin-bound fraction to the nucleoplasm and then on to the cytoplasm. Transcript maturation was similar for lncRNAs in the chromatin fraction, but was significantly lower in the nucleoplasm and cytoplasm. LncRNA transcripts were 11-fold more likely to be significantly enriched in the nucleus than cytoplasm, and 100-fold more likely to be significantly chromatin-bound than nucleoplasmic. Sequencing chromatin-bound RNA greatly increased the sensitivity for detecting lowly expressed lncRNAs and enabled us to discover and localize hundreds of novel regulated liver lncRNAs, including lncRNAs showing sex-biased expression or responsiveness to TCPOBOP a xenobiotic agonist ligand of constitutive androstane receptor (Nr1i3).

CONCLUSIONS

Integration of our findings with prior studies and lncRNA annotations identified candidate regulatory lncRNAs for a variety of hepatic functions based on gene co-localization within topologically associating domains or transcription divergent or antisense to PCGs associated with pathways linked to hepatic physiology and disease.

Biomolecular Condensates and Gene Activation in Development and Disease

Activating the right gene at the right time and place is essential for development. Emerging evidence suggests that this process is regulated by the mesoscale compartmentalization of the gene-control machinery, RNA polymerase II and its cofactors, within biomolecular condensates. Coupling gene activity to the reversible and dynamic process of condensate formation is proposed to enable the robust and precise changes in gene-regulatory programs during signaling and development. The macromolecular features that enable condensates and the regulatory pathways that control them are dysregulated in disease, highlighting their importance for normal physiology. In this review, we will discuss the role of condensates in gene activation; the multivalent features of protein, RNA, and DNA that enable reversible condensate formation; and how these processes are utilized in normal and disease biology. Understanding the regulation of condensates promises to provide novel insights into how organization of the gene-control machinery regulates development and disease.