Non-coding Transcripts from Enhancers: New Insights into Enhancer Activity and Gene Expression Regulation

The long non-coding RNAs (lncRNA) in the pathogenesis of gastric cancer cells: molecular mechanisms and involvement miRNAs

A complex sequence of occurrences, including host genetic vulnerability, Helicobacter pylori infection, and other environmental variables, culminate in gastric cancer (GC). The development of several genetic and epigenetic changes in oncogenes and tumor suppressor genes causes dysregulation of several signaling pathways, which upsets the cell cycle and the equilibrium between cell division and apoptosis, leading to GC. Developments in computational biology and RNA-seq technology enable quick detection and characterization of long non-coding RNAs (lncRNAs). Recent studies have shown that long non-coding RNAs (lncRNAs) have multiple roles in the development of gastric cancer. These lncRNAs interact with molecules of protein, RNA, DNA, and/or combinations. This review article explores several gastric cancer-associated lncRNAs, such as ADAMTS9-AS2, UCA1, XBP-1, and LINC00152. These various lncRNAs could change GC cell apoptosis, migration, and invasion features in the tumor microenvironment. This review provides an overview of the most recent research on lncRNAs and GC cell apoptosis, migration, invasion, and drug resistance, focusing on studies conducted in cancer cells and healthy cells during differentiation.

Identification of key eRNAs for intervertebral disc degeneration by integrated multinomial bioinformatics analysis

BACKGROUND

Intervertebral disc degeneration (IVDD) is a common degenerative condition leading to abnormal stress distribution under load, causing intervertebral stenosis, facet joint degeneration, and foraminal stenosis. Very little is known about the molecular mechanism of eRNAs in IVDD.

METHODS

Gene expression profiles of 38 annulus disc samples composed of 27 less degenerated discs (LDs) and 11 more degenerated discs (MDs) were retrieved from the GEO database. Then, differentially expressed enhancer RNAs (DEeRNAs), differentially expressed target genes (DETGs), and differentially expressed transcription factors (DETFs), hallmark of cancer signalling pathways according to GSVA; the types and quantity of immune cells according to CIBERSORT; and immune gene sets according to ssGSEA were analysed to construct an IVDD-related eRNA network. Then, multidimensional validation was performed to explore the interactions among DEeRNAs, DETFs and DEGs in space.

RESULTS

A total of 53 components, 14 DETGs, 15 DEeRNAs, 3 DETFs, 5 immune cells, 9 hallmarks, and 7 immune gene sets, were selected to construct the regulatory network. After validation by online multidimensional databases, 21 interactive DEeRNA-DEG-DETF axes related to IVDD exacerbation were identified, among which the C1S-CTNNB1-CHD4 axis was the most significant.

CONCLUSION

Based upon the results of our study, we theorize that the C1S-CTNNB1-CHD4 axis plays a vital role in IVDD exacerbation. Specifically, C1S recruits CTNNB1 and upregulates the expression of CHD4 in IVDD, and subsequently, CHD4 suppresses glycolysis and activates oxidative phosphorylation, thus generating insoluble collagen fibre deposits and leading to the progression of IVDD. Overall, these DEeRNAs could comprise promising therapeutic targets for IVDD due to their high tissue specificity.

Identification of enhancer RNA AC005515.1 as a novel biomarker for prognosis in esophageal cancer and predictors of immunotherapy response

Background

The enhancer RNA (eRNA) signature shows excellent promise in the prognostic role of many malignancies, but its value has not been fully explored in esophageal cancer (ESCA).

Methods

We comprehensively analyzed 33 oncogene expression matrices and clinical data from The Cancer Genome Atlas (TCGA) and identified ESCA prognostic-related key eRNAs by Kaplan-Meier and co-expression analysis. We also investigated the prognostic role of the key eRNA using a series of bioinformatics approaches, including immune infiltration, immune function, immune subtypes, and the tumor microenvironment. Finally, the tumor immune dysfunction and exclusion (TIDE) score was used to predict the immune response to immune checkpoint blockade (ICB) therapy.

Results

We identified eRNA AC005515.1, AC012368.1, AP003469.2, Clorf61, and WDFY3-AS2 were associated with the prognosis of ESCA. AC005515.1 is a critical prognostic-related eRNA in ESCA and was significantly co-expressed with immune checkpoint genes (CTLA4, CD274, etc.). In the pan-cancer analysis, AC005515.1 was also associated with the prognosis of seven cancers, including kidney renal papillary cell carcinoma and low-grade brain glioma. It was also found to be co-expressed with immune checkpoint genes in these tumors. Moreover, high expression of AC005515.1 was associated with CD8+ T cells and M1 macrophages infiltration, and the AC005515.1 high-expression group had a higher TIDE score in ESCA.

Conclusions

Overall, eRNA AC005515.1 is associated with the local immune environment of ESCA and may become a new biomarker of ESCA prognosis and immunotherapy response.

Enhancer Function in the 3D Genome

In this review, we consider various aspects of enhancer functioning in the context of the 3D genome. Particular attention is paid to the mechanisms of enhancer-promoter communication and the significance of the spatial juxtaposition of enhancers and promoters in 3D nuclear space. A model of an activator chromatin compartment is substantiated, which provides the possibility of transferring activating factors from an enhancer to a promoter without establishing direct contact between these elements. The mechanisms of selective activation of individual promoters or promoter classes by enhancers are also discussed.

Single-cell RNA sequencing identifies distinct transcriptomic signatures between PMA/ionomycin- and αCD3/αCD28-activated primary human T cells

Immunologists have activated T cells in vitro using various stimulation methods, including phorbol myristate acetate (PMA)/ionomycin and αCD3/αCD28 agonistic antibodies. PMA stimulates protein kinase C, activating nuclear factor-κB, and ionomycin increases intracellular calcium levels, resulting in activation of nuclear factor of activated T cell. In contrast, αCD3/αCD28 agonistic antibodies activate T cells through ZAP-70, which phosphorylates linker for activation of T cell and SH2-domain-containing leukocyte protein of 76 kD. However, despite the use of these two different in vitro T cell activation methods for decades, the differential effects of chemical-based and antibody-based activation of primary human T cells have not yet been comprehensively described. Using single-cell RNA sequencing (scRNA-seq) technologies to analyze gene expression unbiasedly at the single-cell level, we compared the transcriptomic profiles of the non-physiological and physiological activation methods on human peripheral blood mononuclear cell-derived T cells from four independent donors. Remarkable transcriptomic differences in the expression of cytokines and their respective receptors were identified. We also identified activated CD4 T cell subsets (CD55+) enriched specifically by PMA/ionomycin activation. We believe this activated human T cell transcriptome atlas derived from two different activation methods will enhance our understanding, highlight the optimal use of these two in vitro T cell activation assays, and be applied as a reference standard when analyzing activated specific disease-originated T cells through scRNA-seq.

Integrative profiling of gene expression and chromatin accessibility elucidates specific transcriptional networks in porcine neutrophils

Neutrophils are vital components of the immune system for limiting the invasion and proliferation of pathogens in the body. Surprisingly, the functional annotation of porcine neutrophils is still limited. The transcriptomic and epigenetic assessment of porcine neutrophils from healthy pigs was performed by bulk RNA sequencing and transposase accessible chromatin sequencing (ATAC-seq). First, we sequenced and compared the transcriptome of porcine neutrophils with eight other immune cell transcriptomes to identify a neutrophil-enriched gene list within a detected neutrophil co-expression module. Second, we used ATAC-seq analysis to report for the first time the genome-wide chromatin accessible regions of porcine neutrophils. A combined analysis using both transcriptomic and chromatin accessibility data further defined the neutrophil co-expression network controlled by transcription factors likely important for neutrophil lineage commitment and function. We identified chromatin accessible regions around promoters of neutrophil-specific genes that were predicted to be bound by neutrophil-specific transcription factors. Additionally, published DNA methylation data from porcine immune cells including neutrophils were used to link low DNA methylation patterns to accessible chromatin regions and genes with highly enriched expression in porcine neutrophils. In summary, our data provides the first integrative analysis of the accessible chromatin regions and transcriptional status of porcine neutrophils, contributing to the Functional Annotation of Animal Genomes (FAANG) project, and demonstrates the utility of chromatin accessible regions to identify and enrich our understanding of transcriptional networks in a cell type such as neutrophils.

Comprehensive investigation of the expression profiles of common long noncoding RNAs during microglial activation

Microglia, similar to peripheral macrophages, are the primary immune cells of the central nervous system (CNS). Microglia exist in the resting state in the healthy CNS, but can be activated and polarized into either M1 or M2 subtypes for immune defense and the maintenance of CNS homeostasis by multiple stimuli. Several long noncoding RNAs (lncRNAs) mediate human inflammatory diseases and neuropathologies by regulating their target genes. However, the function of common lncRNAs that contribute to microglial activation remains unclear. Thus, we used bioinformatic approaches to identify common lncRNAs involved in microglial activation in vitro. Our study identified several lncRNAs as common regulators of microglial activation. We identified 283 common mRNAs and 53 common lncRNAs during mouse M1 microglial activation processes, whereas 26 common mRNAs and five common lncRNAs were identified during mouse M2 microglial activation processes. A total of 648 common mRNAs and 274 common lncRNAs were identified during the activation of human M1 microglia. In addition, we identified 1,920 common co-expressed pairs in mouse M1 activation processes and 25 common co-expressed pairs in mouse M2 activation processes. Our study provides a comprehensive understanding of common lncRNA expression profiles in microglial activation processes in vitro. The list of common lncRNAs identified in this study provides novel evidence and clues regarding the molecular mechanisms underlying microglial activation.

Interplay between mammalian heat shock factors 1 and 2 in physiology and pathology

The heat-shock factors (HSFs) belong to an evolutionary conserved family of transcription factors that were discovered already over 30 years ago. The HSFs have been shown to a have a broad repertoire of target genes, and they also have crucial functions during normal development. Importantly, HSFs have been linked to several disease states, such as neurodegenerative disorders and cancer, highlighting their importance in physiology and pathology. However, it is still unclear how HSFs are regulated and how they choose their specific target genes under different conditions. Posttranslational modifications and interplay among the HSF family members have been shown to be key regulatory mechanisms for these transcription factors. In this review, we focus on the mammalian HSF1 and HSF2, including their interplay, and provide an updated overview of the advances in understanding how HSFs are regulated and how they function in multiple processes of development, aging, and disease. We also discuss HSFs as therapeutic targets, especially the recently reported HSF1 inhibitors.

Oxidative Stress and Its Modulation by Ladostigil Alter the Expression of Abundant Long Non-Coding RNAs in SH-SY5Y Cells

Neurodegenerative disorders, brain injury, and the decline in cognitive function with aging are accompanied by a reduced capacity of cells in the brain to cope with oxidative stress and inflammation. In this study, we focused on the response to oxidative stress in SH-SY5Y, a human neuroblastoma cell line. We monitored the viability of the cells in the presence of oxidative stress. Such stress was induced by hydrogen peroxide or by Sin1 (3-morpholinosydnonimine) that generates reactive oxygen and nitrogen species (ROS and RNS). Both stressors caused significant cell death. Our results from the RNA-seq experiments show that SH-SY5Y cells treated with Sin1 for 24 h resulted in 94 differently expressed long non-coding RNAs (lncRNAs), including many abundant ones. Among the abundant lncRNAs that were upregulated by exposing the cells to Sin1 were those implicated in redox homeostasis, energy metabolism, and neurodegenerative diseases (e.g., MALAT1, MIAT, GABPB1-AS1, NEAT1, MIAT, GABPB1-AS1, and HAND2-AS1). Another group of abundant lncRNAs that were significantly altered under oxidative stress included cancer-related SNHG family members. We tested the impact of ladostigil, a bifunctional reagent with antioxidant and anti-inflammatory properties, on the lncRNA expression levels. Ladostigil was previously shown to enhance learning and memory in the brains of elderly rats. In SH-SY5Y cells, several lncRNAs involved in transcription regulation and the chromatin structure were significantly induced by ladostigil. We anticipate that these poorly studied lncRNAs may act as enhancers (eRNA), regulating transcription and splicing, and in competition for miRNA binding (ceRNA). We found that the induction of abundant lncRNAs, such as MALAT1, NEAT-1, MIAT, and SHNG12, by the Sin1 oxidative stress paradigm specifies only the undifferentiated cell state. We conclude that a global alteration in the lncRNA profiles upon stress in SH-SY5Y may shift cell homeostasis and is an attractive in vitro system to characterize drugs that impact the redox state of the cells and their viability.

ISR8/IRF1-AS1 Is Relevant for IFNα and NF-κB Responses

The study of the interferon (IFN) α-induced cell transcriptome has shown altered expression of several long non-coding RNAs (lncRNAs). ISR8/IRF1-AS1 (IFN stimulated RNA 8), located close to IFN regulatory factor 1 (IRF1) coding gene, transcribes a lncRNA induced at early times after IFNα treatment or IRF1 or NF-κB activation. Depletion or overexpression of ISR8 RNA does not lead to detected deregulation of the IFN response. Surprisingly, disruption of ISR8 locus with CRISPR-Cas9 genome editing results in cells that fail to induce several key ISGs and pro-inflammatory cytokines after a trigger with IFNα or overexpression of IRF1 or the NF-κB subunit RELA. This suggests that the ISR8 locus may play a relevant role in IFNα and NF-κB pathways. Interestingly, IFNα, IRFs and NF-κB-responding luciferase reporters are normally induced in ISR8-disrupted cells when expressed from a plasmid but not when integrated into the genome. Therefore, IFNα and NF-κB pathways are functional to induce the expression of exogenous episomic transcripts but fail to activate transcription from genomic promoters. Transcription from these promoters is not restored with silencing inhibitors, by decreasing the levels of several negative regulators or by overexpression of inducers. Transcriptome analyses indicate that ISR8-disrupted cells have a drastic increase in the levels of negative regulators such as XIST and Zinc finger proteins. Our results agree with ISR8 loci being an enhancer region that is fundamental for proper antiviral and proinflammatory responses. These results are relevant because several SNPs located in the ISR8 region are associated with chronic inflammatory and autoimmune diseases including Crohn’s disease, inflammatory bowel disease, ulcerative colitis or asthma.

Epigenetic regulation of epithelial to mesenchymal transition: a trophoblast perspective

Epigenetic changes alter expression of genes at both pre- and post-transcriptional levels without changing their DNA sequence. Accumulating evidence suggests that such changes can modify cellular behaviour and characteristics required during development and in response to various extracellular stimuli. Trophoblast cells develop from the outermost trophectoderm layer of the blastocyst and undergo many phenotypic changes as the placenta develops. One such phenotypic change is differentiation of the epithelial natured cytotrophoblasts into the mesenchymal natured extravillous trophoblasts. The extravillous trophoblasts are primarily responsible for invading into the maternal decidua and thus establishing connection with the maternal spiral arteries. Any dysregulation of this process can have adverse effects on the pregnancy outcome. Hence, tight regulation of this epithelial-mesenchymal transition is critical for successful pregnancy. This review summarizes the recent research on the epigenetic regulation of the epithelial-mesenchymal transition occurring in the trophoblast cells during placental development. The functional significance of chemical modifications of DNA and histone, which regulate transcription, as well as non-coding RNAs, which control gene expression post-transcriptionally, is discussed in relation to trophoblast biology.

Progress in Expression Pattern and Molecular Regulation Mechanism of LncRNA in Bovine Mastitis

Simple Summary

Bovine mastitis is an inflammatory disease of the mammary glands that causes serious harm to cow health and huge economic losses. Susceptibility or resistance to mastitis in individual cows is mainly determined by genetic factors, including coding genes and non-coding genes. Long non-coding RNAs (lncRNAs) are non-coding RNA molecules with a length of more than 200 nucleotides (nt) that have recently been discovered. They can regulate a variety of diseases of humans and animals, especially the immune response and inflammatory disease process. This paper reviews the role of long non-coding RNA (lncRNA) in inflammatory diseases, emphasizes on the latest research progress of lncRNA expression and the molecular regulatory mechanism in bovine mastitis, and looks forward to the research and application prospect of lncRNA in bovine mastitis, intending to provide a reference for scientific researchers to systematically understand this research field.

Abstract

Bovine mastitis is an inflammatory disease caused by pathogenic microbial infection, trauma, or other factors. Its morbidity is high, and it is difficult to cure, causing great harm to the health of cows and the safety of dairy products. Susceptibility or resistance to mastitis in individual cows is mainly determined by genetic factors, including coding genes and non-coding genes. Long non-coding RNAs (lncRNAs) are a class of endogenous non-coding RNA molecules with a length of more than 200 nucleotides (nt) that have recently been discovered. They can regulate the immune response of humans and animals on three levels (transcription, epigenetic modification, and post-transcription), and are widely involved in the pathological process of inflammatory diseases. Over the past few years, extensive findings revealed basic roles of lncRNAs in inflammation, especially bovine mastitis. This paper reviews the expression pattern and mechanism of long non-coding RNA (lncRNA) in inflammatory diseases, emphasizes on the latest research progress of the lncRNA expression pattern and molecular regulatory mechanism in bovine mastitis, analyzes the molecular regulatory network of differentially expressed lncRNAs, and looks forward to the research and application prospect of lncRNA in bovine mastitis, laying a foundation for molecular breeding and the biological therapy of bovine mastitis.

An Immune-Related Long Noncoding RNA Signature as a Prognostic Biomarker for Human Endometrial Cancer

Background

Endometrial cancer is among the most common malignant tumors threatening the health of women. Recently, immunity and long noncoding RNA (lncRNA) have been widely examined in oncology and shown to play important roles in oncology. Here, we searched for immune-related lncRNAs as prognostic biomarkers to predict the outcome of patients with endometrial cancer.

Methods

RNA sequencing data for 575 endometrial cancer samples and immune-related genes were downloaded from The Cancer Genome Atlas (TCGA) database and gene set enrichment analysis (GSEA) gene sets, respectively. Immune-related lncRNAs showing a coexpression relationship with immune-related genes were obtained, and Cox regression analysis was performed to construct the prognostic model. Survival, independent prognostic, and clinical correlation analyses were performed to evaluate the prognostic model. Immune infiltration of endometrial cancer samples was also evaluated. Functional annotation of 12 immune-related lncRNAs was performed using GSEA software. Prognostic nomogram and survival analysis for independent prognostic risk factors were performed to evaluate the prognostic model and calculate the survival time based on the prognostic model.

Results

Twelve immune-related lncRNAs (ELN-AS1, AC103563.7, PCAT19, AF131215.5, LINC01871, AC084117.1, NRAV, SCARNA9, AL049539.1, POC1B-AS1, AC108134.4, and AC019080.5) were obtained, and a prognostic model was constructed. The survival rate in the high-risk group was significantly lower than that in the low-risk group. Patient age, pathological grade, the International Federation of Gynecology and Obstetrics (FIGO) stage, and risk status were the risk factors. The 12 immune-related lncRNAs correlated with patient age, pathological grade, and FIGO stage. Principal component analysis and functional annotation showed that the high-risk and low-risk groups separated better, and the immune status of the high-risk and low-risk groups differed. Nomogram and receiver operating characteristic (ROC) curves effectively predicted the prognosis of endometrial cancer. Additionally, age, pathological grade, FIGO stage, and risk status were all related to patient survival.

Conclusion

We identified 12 immune-related lncRNAs affecting the prognosis of endometrial cancer, which may be useful as therapeutic targets and molecular biomarkers.

EMX2OS plays a prognosis-associated enhancer RNA role in gastric cancer

Enhancer RNAs (eRNAs), a subclass of lncRNAs, are derived from enhancer regions. The function of eRNAs has been reported by many previous studies. However, the role of eRNAs in gastric cancer, especially the prognosis-associated eRNAs, has not been studied yet.In this study, we have used a novel approach to screened key eRNAs in gastric cancer. Kaplan-Meier correlation analysis and Co-expression analysis were used to find the most significant survival-associated eRNAs. Enrichment analysis is applied to explore the key functions and pathways of screened eRNAs. The correlation and survival analysis are used to evaluate targeted genes in the pan-cancer analysisA total of 63 prognostic-associated eRNAs in gastric cancer were identified, the top 6 eRNAs were LINC01714, ZNF192P1, AC079760.2, LINC01645, EMX2OS, and AC114489.2. The correlation analysis demonstrated the top 10 screened eRNAs and their targeted genes. The results demonstrated that EMX2OS was ranked as the top eRNA according to the results of the Kaplan-Meier analysis. The correlation analysis demonstrated that eRNA EMX2OS is correlated with age, grade, stage, and cancer status. The pan-cancer analysis demonstrated that EMX2OS was associated with poor survival outcomes in adrenocortical carcinoma, cervical squamous cell carcinoma and endocervical adenocarcinoma, kidney renal clear cell carcinoma, stomach adenocarcinoma, and uveal melanoma.In this study, survival-related eRNAs were screened and the correlation between survival-related eRNAs and their targeted genes was demonstrated. EMX2OS plays a prognosis-associated eRNA role in gastric cancer, which might be a novel therapeutic target in clinical practice.

Animal-eRNAdb: a comprehensive animal enhancer RNA database

Enhancer RNAs (eRNAs) are a class of non-coding RNAs transcribed from enhancers. As the markers of active enhancers, eRNAs play important roles in gene regulation and are associated with various complex traits and characteristics. With increasing attention to eRNAs, numerous eRNAs have been identified in different human tissues. However, the expression landscape, regulatory network and potential functions of eRNAs in animals have not been fully elucidated. Here, we systematically characterized 185 177 eRNAs from 5085 samples across 10 species by mapping the RNA sequencing data to the regions of known enhancers. To explore their potential functions based on evolutionary conservation, we investigated the sequence similarity of eRNAs among multiple species. In addition, we identified the possible associations between eRNAs and transcription factors (TFs) or nearby genes to decipher their possible regulators and target genes, as well as characterized trait-related eRNAs to explore their potential functions in biological processes. Based on these findings, we further developed Animal-eRNAdb (http://gong_lab.hzau.edu.cn/Animal-eRNAdb/), a user-friendly database for data searching, browsing and downloading. With the comprehensive characterization of eRNAs in various tissues of different species, Animal-eRNAdb may greatly facilitate the exploration of functions and mechanisms of eRNAs.

Identification and Validation of an Immune-Related eRNA Prognostic Signature for Hepatocellular Carcinoma

Background

Enhancer RNAs (eRNAs) are intergenic long non-coding RNAs (lncRNAs) that participate in the progression of malignancies by targeting tumor-related genes and immune checkpoints. However, the potential role of eRNAs in hepatocellular carcinoma (HCC) is unclear. In this study, we aimed to construct an immune-related eRNA prognostic model that could be used to prospectively assess the prognosis of patients with HCC.

Methods

Gene expression profiles of patients with HCC were downloaded from The Cancer Genome Atlas (TCGA). The eRNAs co-expressed from immune genes were identified as immune-related eRNAs. Cox regression analyses were applied in a training cohort to construct an immune-related eRNA signature (IReRS), that was subsequently used to analyze a testing cohort and combination of the two cohorts. Kaplan-Meier and receiver operating characteristic (ROC) curves were used to validate the predictive effect in the three cohorts. Gene Set Enrishment Analysis (GSEA) computation was used to identify an IReRS-related signaling pathway. A web-based cell type identification by estimating relative subsets of RNA transcripts (CIBERSORT) computation was used to evaluate the relationship between the IReRS and infiltrating immune cells.

Results

A total of sixty-four immune-related eRNAs (IReRNAs) was identified in HCC, and 14 IReRNAs were associated with overall survival (OS). Five IReRNAs were used for constructing an immune-related eRNA signature (IReRS), which was shown to correlate with poor survival and to be an independent prognostic biomarker for HCC. The GSEA results showed that the IReRS was correlated to cancer-related and immune-related pathways. Moreover, we found that IReRS was correlated to infiltrating immune cells, including CD8⁺ T cells and M0 macrophages. Finally, differential expressions of the five risk IReRNAs in tumor tissues vs. adjacent normal tissues and their prognostic values were verified, in which the AL445524.1 may function as an oncogene that affects prognosis partly by regulating CD4-CLTA4 related genes.

Conclusion

Our results suggest that the IReRS could serve as a biomarker for predicting prognosis in patients with HCC. Additionally, it may be correlated to the tumor immune microenvironment and could also be used as a biomarker in immunotherapy for HCC.

Unveiling the long non-coding RNA profile of porcine reproductive and respiratory syndrome virus-infected porcine alveolar macrophages

Background

Long noncoding RNA (lncRNA) is highly associated with inflammatory response and virus-induced interferon production. By far the majority of studies have focused on the immune-related lncRNAs of mice and humans, but the function of lncRNAs in porcine immune cells are poorly understood. Porcine reproductive and respiratory syndrome virus (PRRSV) impairs local immune responses in the lungs of nursery and growing pigs, whereas the virus triggers the inflammatory responses. Porcine alveolar macrophage (PAM) is the primary target cell of PRRSV, thus PRRSV is used as an in vitro model of inflammation. Here, we profiled lncRNA and mRNA repertories from PRRSV-infected PAMs to explore the underlying mechanism of porcine lncRNAs in regulating host immune responses.

Results

In this study, a total of 350 annotated lncRNAs and 1792 novel lncRNAs in PAMs were identified through RNA-seq analysis. Among them 86 differentially expressed (DE) lncRNAs and 406 DE protein-coding mRNAs were identified upon PRRSV incubation. GO category and KEGG pathway enrichment analyses revealed that these DE lncRNAs and mRNAs were mainly involved in inflammation- and pathogen infection-induced pathways. The results of dynamic correlated expression networks between lncRNAs and their predicted target genes uncovered that numerous lncRNAs, such as XLOC-022175, XLOC-019295, and XLOC-017089, were correlated with innate immune genes. Further analysis validated that these three lncRNAs were positively correlated with their predicted target genes including CXCL2, IFI6, and CD163. This study suggests that porcine lncRNAs affect immune responses against PRRSV infection through regulating their target genes in PAMs.

Conclusion

This study provides both transcriptomic and epigenetic status of porcine macrophages. In response to PRRSV infection, comprehensive DE lncRNAs and mRNAs were profiled from PAMs. Co-expression analysis demonstrated that lncRNAs are emerging as the important modulators of immune gene activities through their critical influence upon PRRSV infection in porcine macrophages.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07482-9.

Long noncoding RNA LINC00518 induces radioresistance by regulating glycolysis through an miR-33a-3p/HIF-1α negative feedback loop in melanoma

The long noncoding RNA, LINC00518, is highly expressed in various types of cancers and is involved in cancer progression. Although LINC00518 promotes the metastasis of cutaneous malignant melanoma (CMM), the mechanism underlaying its effects on CMM radiosensitivity remains unclear. In this study, LINC00518 expression was significantly upregulated in CMM samples, and LINC00518 levels were associated with poor prognosis of patients with CMM. Knockdown of LINC00518 in CMM cells significantly inhibited cell invasion, migration, proliferation, and clonogenicity. LINC00518-mediated invasion, migration, proliferation, and clonogenicity were negatively regulated by the microRNA, miR-33a-3p, in vitro, which increased sensitivity to radiotherapy via inhibition of the hypoxia-inducible factor 1α (HIF-1α)/lactate dehydrogenase A glycolysis axis. Additionally, HIF-1α recognized the miR-33a-3p promoter region and recruited histone deacetylase 2, which decreased the expression of miR-33a-3p and formed an LINC00518/miR-33a-3p/HIF-1α negative feedback loop. Furthermore, signaling with initially activated glycolysis and radioresistance in CMM cells was impaired by Santacruzamate A, a histone deacetylase inhibitor, and 2-deoxy-D-glucose, a glycolytic inhibitor. Lastly, knockdown of LINC00518 expression sensitized CMM cancer cells to radiotherapy in an in vivo subcutaneously implanted tumor model. In conclusion, LINC00518 was confirmed to be an oncogene in CMM, which induces radioresistance by regulating glycolysis through an miR-33a-3p/HIF-1α negative feedback loop. Our study, may provide a potential strategy to improve the treatment outcome of radiotherapy in CMM.

The Biological Roles of lncRNAs and Future Prospects in Clinical Application

Chemo and radiation therapies are the most commonly used therapies for cancer, but they can induce DNA damage, resulting in the apoptosis of host cells. DNA double-stranded breaks (DSBs) are the most lethal form of DNA damage in cells, which are constantly caused by a wide variety of genotoxic agents, both environmentally and endogenously. To maintain genomic integrity, eukaryotic organisms have developed a complex mechanism for the repair of DNA damage. Researches reported that many cellular long noncoding RNAs (lncRNAs) were involved in the response of DNA damage. The roles of lncRNAs in DNA damage response can be regulated by the dynamic modification of N6-adenosine methylation (m6A). The cellular accumulation of DNA damage can result in various diseases, including cancers. Additionally, lncRNAs also play roles in controlling the gene expression and regulation of autophagy, which are indirectly involved with individual development. The dysregulation of these functions can facilitate human tumorigenesis. In this review, we summarized the origin and overview function of lncRNAs and highlighted the roles of lncRNAs involved in the repair of DNA damage.

An integrated analysis of enhancer RNAs in glioma and a validation of their prognostic values

Glioma, a highly aggressive neuroepithelial malignant brain tumor, is associated with high disability and recurrence rates. Enhancer RNA (eRNA) plays a significant role in tumor proliferation and metastasis; however, their functions in gliomas need further evaluation. We used the computational pipeline, PreSTIGE, to predict tissue-specific enhancer-derived RNAs and the underlying regulatory genes. Using data retrieved from the TCGA and CGGA databases, a LASSO regression analysis and multiCox proportional hazards regression analyses were performed to determine the hub eRNAs associated with glioma prognosis. Quantitative reverse transcription PCR was performed on the glioma samples to evaluate the expression characteristics of the identified hub eRNAs. To construct a risk signature, we selected three eRNAs, including CRNDE, MRPS31P5, and LINC00844, for their significant prognostic values. The predictive value of the risk signature was validated using the CGGA and Rembrandt cohorts. Apart from the risk signature, the nomogram performed well at predicting OS in glioma patients. An eRNA-target gene regulatory network was established, which we evaluated using a target gene enrichment analysis. Pathway and gene ontology (GO) analyses demonstrated that the risk signature is associated with mRNA processing and spliceosome in glioma. Furthermore, we found that hub eRNAs potentially regulate the expressions of numerous splicing factors, such as MOV10 and SEC31B, and are correlated with prognosis-associated alteration splicing (AS). In conclusion, we established a risk signature that comprises three eRNAs, which can accurately be utilized as targets to predict prognosis in glioma patients.

What do we know about the role of lncRNAs in multiple sclerosis?

Multiple sclerosis is a chronic, inflammatory and degenerative disease of the central nervous system of unknown aetiology although well-defined evidence supports an autoimmune pathogenesis. So far, the exact mechanisms leading to autoimmune diseases are still only partially understood. We know that genetic, epigenetic, molecular, and cellular factors resulting in pathogenic inflammatory responses are certainly involved. Long non-coding RNAs (lncRNAs) are non-protein coding transcripts longer than 200 nucleotides that play an important role in both innate and acquired immunity, so there is great interest in lncRNAs involved in autoimmune diseases. The research on multiple sclerosis has been enriched with many studies on the molecular role of lncRNAs in the pathogenesis of the disease and their potential application as diagnostic and prognostic biomarkers. In particular, many multiple sclerosis fields of research are based on the identification of lncRNAs as possible biomarkers able to predict the onset of the disease, its activity degree, its progression phase and the response to disease-modifying drugs. Last but not least, studies on lncRNAs can provide a new molecular target for new therapies, missing, so far, a cure for multiple sclerosis. While our knowledge on the role of lncRNA in multiple sclerosis has recently improved, further studies are required to better understand the specific role of lncRNAs in this neurological disease. In this review, we present the most recent studies on molecular characterization of lncRNAs in multiple sclerosis disorder discussing their clinical relevance as biomarkers for diagnosis and treatments.

Characterization of complete lncRNAs transcriptome reveals the functional and clinical impact of lncRNAs in multiple myeloma

Multiple myeloma (MM) is an incurable disease, whose clinical heterogeneity makes its management challenging, highlighting the need for biological features to guide improved therapies. Deregulation of specific long non-coding RNAs (lncRNAs) has been shown in MM, nevertheless, the complete lncRNA transcriptome has not yet been elucidated. In this work, we identified 40,511 novel lncRNAs in MM samples. lncRNAs accounted for 82% of the MM transcriptome and were more heterogeneously expressed than coding genes. A total of 10,351 overexpressed and 9,535 downregulated lncRNAs were identified in MM patients when compared with normal bone-marrow plasma cells. Transcriptional dynamics study of lncRNAs in the context of normal B-cell maturation revealed 989 lncRNAs with exclusive expression in MM, among which 89 showed de novo epigenomic activation. Knockdown studies on one of these lncRNAs, SMILO (specific myeloma intergenic long non-coding RNA), resulted in reduced proliferation and induction of apoptosis of MM cells, and activation of the interferon pathway. We also showed that the expression of lncRNAs, together with clinical and genetic risk alterations, stratified MM patients into several progression-free survival and overall survival groups. In summary, our global analysis of the lncRNAs transcriptome reveals the presence of specific lncRNAs associated with the biological and clinical behavior of the disease.

Silencing of lncRNA MIAT alleviates LPS-induced pneumonia via regulating miR-147a/NKAP/NF-κB axis

Purpose: Pneumonia is a respiratory disease with an increasing incidence in recent years. More and more studies have revealed that lncRNAs can regulate the transcriptional expression of target genes at different stage. Herein, we aimed to explore the effect of lncRNA MIAT in LPS-induced pneumonia, and further illuminate the possible underlying mechanisms.

Method and results: Mice were intraperitoneally injected with LPS, and the lung inflammation was evaluated. Microarray showed lncRNA MIAT was up-regulated in LPS-induced pulmonary inflammation. And qRT-PCR and FISH assay indicated that MIAT was increased in mice with LPS injection. Functional analysis showed sh-MIAT inhibited LPS-induced inflammation response, inhibited apoptosis level and protected lung function. As well, si-MIAT removed the injury of LPS on mouse lung epithelial TC-1 cells, and inhibited the activation of NF-κB signaling. Furthermore, MIAT acted as a sponge of miR-147a, and miR-147a directly targeted NKAP. Functionally, AMO-147a or NKAP remitted the beneficial effects of si-MIAT on LPS-induced inflammation response of TC-1 cells.

Conclusion: Deletion of MIAT protected against LPS-induced lung inflammation via regulating miR-147a/NKAP, which might provide new insight for pneumonia treatment.

Downregulation of enhancer RNA EMX2OS is associated with poor prognosis in kidney renal clear cell carcinoma

Enhancer RNAs are a subclass of long non-coding RNAs transcribed from enhancer regions that play an important role in the transcriptional regulation of genes. However, their role in kidney renal clear cell carcinoma (KIRC) is largely unknown. Herein, we identified the key enhancer RNAs in KIRC via an integrated data analysis method. Gene expression profiles and clinical data of KIRC and 32 other cancer types were acquired using the University of California Santa Cruz Xena platform. Reported enhancer RNAs and genes regulated by them were selected as putative enhancer RNA-target pairs. Kaplan-Meier survival and correlation analyses were performed to identify the key enhancer RNAs. Finally, EMX2OS was identified as the enhancer RNA most associated with survival, with EMX2 as its target. EMX2OS downregulation was significantly associated with higher histological grade, advanced stage, and poorer prognosis. The results were validated in pan-cancer data from The Cancer Genome Atlas and RT-qPCR analysis of 12 pairs of KIRC and normal real-world samples. Functional enrichment analysis indicated that several metabolism-associated signaling pathways were enriched. This study demonstrated that EMX2OS is a key metabolism-associated enhancer RNA in KIRC with a favorable impact on survival and may be a novel therapeutic target in KIRC.

Cis-acting lnc-eRNA SEELA directly binds histone H4 to promote histone recognition and leukemia progression

Background

Long noncoding enhancer RNAs (lnc-eRNAs) are a subset of stable eRNAs identified from annotated lncRNAs. They might act as enhancer activity-related therapeutic targets in cancer. However, the underlying mechanism of epigenetic activation and their function in cancer initiation and progression remain largely unknown.

Results

We identify a set of lncRNAs as lnc-eRNAs according to the epigenetic signatures of enhancers. We show that these lnc-eRNAs are broadly activated in MLL-rearranged leukemia (MLL leukemia), an aggressive leukemia caused by a chromosomal translocation, through a mechanism by which the HOXA cluster initiates enhancer activity, and the epigenetic reader BRD4 cooperates with the coregulator MLL fusion oncoprotein to induce transcriptional activation. To demonstrate the functional roles of lnc-eRNAs, two newly identified lnc-eRNAs transcribed from the SEELA eRNA cluster (SEELA), SEELA1 and SEELA2, are chosen for further studies. The results show that SEELA mediated cis-activated transcription of the nearby oncogene Serine incorporate 2 (SERINC2) by directly binding to the K31 amino acid (aa) of histone H4. Chromatin-bound SEELA strengthens the interaction between chromatin and histone modifiers to promote histone recognition and oncogene transcription. Further studies show that the SEELA-SERINC2 axis regulated aspects of cancer metabolism, such as sphingolipid synthesis, to affect leukemia progression.

Conclusions

This study shows that lnc-eRNAs are epigenetically activated by cancer-initiating oncoproteins and uncovers a cis-activating mechanism of oncogene transcription control based on lnc-eRNA-mediated epigenetic regulation of enhancer activity, providing insights into the critical roles of lnc-eRNAs in cancer initiation and progression.

Supplementary information

Supplementary information accompanies this paper at 10.1186/s13059-020-02186-x.

The emerging role of super enhancer-derived noncoding RNAs in human cancer

Super enhancers (SEs) are large clusters of adjacent enhancers that drive the expression of genes which regulate cellular identity; SE regions can be enriched with a high density of transcription factors, co-factors, and enhancer-associated epigenetic modifications. Through enhanced activation of their target genes, SEs play an important role in various diseases and conditions, including cancer. Recent studies have shown that SEs not only activate the transcriptional expression of coding genes to directly regulate biological functions, but also drive the transcriptional expression of non-coding RNAs (ncRNAs) to indirectly regulate biological functions. SE-derived ncRNAs play critical roles in tumorigenesis, including malignant proliferation, metastasis, drug resistance, and inflammatory response. Moreover, the abnormal expression of SE-derived ncRNAs is closely related to the clinical and pathological characterization of tumors. In this review, we summarize the functions and roles of SE-derived ncRNAs in tumorigenesis and discuss their prospective applications in tumor therapy. A deeper understanding of the potential mechanism underlying the action of SE-derived ncRNAs in tumorigenesis may provide new strategies for the early diagnosis of tumors and targeted therapy.

eRNAs and Superenhancer lncRNAs Are Functional in Human Prostate Cancer

Prostate cancer (PCa) is one of the most commonly diagnosed cancers in males worldwide. lncRNAs (long noncoding RNAs) play a significant role in the occurrence and development of PCa. eRNAs (enhancer RNAs) and SE-lncRNAs (superenhancer lncRNAs) are important elements of lncRNAs, but the role of eRNAs and SE-lncRNAs in PCa remains largely unclear. In this work, we identified 681 eRNAs and 292 SE-lncRNAs that were expressed differentially in PCa using a microarray. We also found that eRNAs transcribed from active open chromatin had significantly higher expression than those from active closed chromatin, and SE-lncRNAs had a little higher expression than eRNAs. Next, we constructed a transcriptional regulation network that eRNA-related enhancer and the target genes shared the same TF-binding motifs. Further, we investigated whether CTCF played a role in mediating the transcriptional regulation network. eRNAs, especially those that regulate androgen response genes, may be candidates for prognostic biomarkers and therapy targets. Our work provides a new perspective for developing medical treatments and therapies for prostate cancer.

Biogenesis mechanisms of circular RNA can be categorized through feature extraction of a machine learning model

Motivation

In recent years, multiple circular RNAs (circRNA) biogenesis mechanisms have been discovered. Although each reported mechanism has been experimentally verified in different circRNAs, no single biogenesis mechanism has been proposed that can universally explain the biogenesis of all tens of thousands of discovered circRNAs. Under the hypothesis that human circRNAs can be categorized according to different biogenesis mechanisms, we designed a contextual regression model trained to predict the formation of circular RNA from a random genomic locus on human genome, with potential biogenesis factors of circular RNA as the features of the training data.

Results

After achieving high prediction accuracy, we found through the feature extraction technique that the examined human circRNAs can be categorized into seven subgroups, according to the presence of the following sequence features: RNA editing sites, simple repeat sequences, self-chains, RNA binding protein binding sites and CpG islands within the flanking regions of the circular RNA back-spliced junction sites. These results support all of the previously reported biogenesis mechanisms of circRNA and solidify the idea that multiple biogenesis mechanisms co-exist for different subset of human circRNAs. Furthermore, we uncover a potential new links between circRNA biogenesis and flanking CpG island. We have also identified RNA binding proteins putatively correlated with circRNA biogenesis.

Availability and implementation

Scripts and tutorial are available at http://wanglab.ucsd.edu/star/circRNA. This program is under GNU General Public License v3.0.

Supplementary information

Supplementary data are available at Bioinformatics online.

Molecular Basis of the Function of Transcriptional Enhancers

Transcriptional enhancers are major genomic elements that control gene activity in eukaryotes. Recent studies provided deeper insight into the temporal and spatial organization of transcription in the nucleus, the role of non-coding RNAs in the process, and the epigenetic control of gene expression. Thus, multiple molecular details of enhancer functioning were revealed. Here, we describe the recent data and models of molecular organization of enhancer-driven transcription.

Long non-coding RNAs and TGF-β signaling in cancer

Cancer is driven by genetic mutations in oncogenes and tumor suppressor genes and by cellular events that develop a misregulated molecular microenvironment in the growing tumor tissue. The tumor microenvironment is guided by the excessive action of specific cytokines including transforming growth factor‐β (TGF‐β), which normally controls embryonic development and the homeostasis of young or adult tissues. As a consequence of the genetic alterations generating a given tumor, TGF‐β can preserve its homeostatic function and attempt to limit neoplastic expansion, whereas, once the tumor has progressed to an aggressive stage, TGF‐β can synergize with various oncogenic stimuli to facilitate tumor invasiveness and metastasis. TGF‐β signaling mechanisms via Smad proteins, various ubiquitin ligases, and protein kinases are relatively well understood. Such mechanisms regulate the expression of genes encoding proteins or non‐coding RNAs. Among non‐coding RNAs, much has been understood regarding the regulation and function of microRNAs, whereas the role of long non‐coding RNAs is still emerging. This article emphasizes TGF‐β signaling mechanisms leading to the regulation of non‐coding genes, the function of such non‐coding RNAs as regulators of TGF‐β signaling, and the contribution of these mechanisms in specific hallmarks of cancer.

Structure encoding in DNA

It is proposed that transposons and related long non-coding RNA define the fine structure of body parts. Although morphogens have long been known to direct the formation of many gross structures in early embryonic development, they do not have the necessary precision to define a structure down to the individual cellular level. Using the distinction between procedural and declarative knowledge in information processing as an analogy, it is hypothesized that DNA encodes fine structure in a manner that is different from the genetic code for proteins. The hypothesis states that repeated or near-repeated sequences that are in transposons and non-coding RNA define body part structures. As the cells in a body part go through the epigenetic process of differentiation, the action of methylation serves to inactivate all but the relevant structure definitions and some associated cell type genes. The transposons left active will then physically modify the DNA sequence in the heterochromatin to establish the local context in the three-dimensional body part structure. This brings the encoded definition of the cell type to the histone. The histone code for that cell type starts the regulatory cascade that turns on the genes associated with that particular type of cell, transforming it from a multipotent cell to a fully differentiated cell. This mechanism creates structures in the musculoskeletal system, the organs of the body, the major parts of the brain, and other systems.

lncRNA HOXB-AS3 exacerbates proliferation, migration, and invasion of lung cancer via activating the PI3K-AKT pathway

Lung cancer remains the leading cause of cancer-related death all over the world. In spite of the great advances made in surgery and chemotherapy, the prognosis of lung cancer patients is poor. A substantial fraction of long noncoding RNAs (lncRNAs) can regulate various cancers. A recent study has reported that lncRNA HOXB-AS3 plays a critical role in cancers. However, its biological function remains unclear in lung cancer progression. In the current research, we found HOXB-AS3 was obviously elevated in NSCLC tissues and cells. Functional assays showed that inhibition of HOXB-AS3 was able to repress A549 and H1975 cell proliferation, cell colony formation ability and meanwhile, triggered cell apoptosis. Furthermore, the lung cancer cell cycle was mostly blocked in the G1 phase whereas the cell ratio in the S phase was reduced. Also, A549 and H1975 cell migration and invasion capacity were significantly repressed by the loss of HOXB-AS3. The PI3K/AKT pathway has been implicated in the carcinogenesis of multiple cancers. Here, we displayed that inhibition of HOXB-AS3 suppressed lung cancer cell progression via inactivating the PI3K/AKT pathway. Subsequently, in vivo experiments were utilized in our study and it was demonstrated that HOXB-AS3 contributed to lung cancer tumor growth via modulating the PI3K/AKT pathway. Overall, we implied that HOXB-AS3 might provide a new perspective for lung cancer treatment via targeting PI3K/AKT.

A Single Cell but Many Different Transcripts: A Journey into the World of Long Non-Coding RNAs

In late 2012 it was evidenced that most of the human genome is transcribed but only a small percentage of the transcripts are translated. This observation supported the importance of non-coding RNAs and it was confirmed in several organisms. The most abundant non-translated transcripts are long non-coding RNAs (lncRNAs). In contrast to protein-coding RNAs, they show a more cell-specific expression. To understand the function of lncRNAs, it is fundamental to investigate in which cells they are preferentially expressed and to detect their subcellular localization. Recent improvements of techniques that localize single RNA molecules in tissues like single-cell RNA sequencing and fluorescence amplification methods have given a considerable boost in the knowledge of the lncRNA functions. In recent years, single-cell transcription variability was associated with non-coding RNA expression, revealing this class of RNAs as important transcripts in the cell lineage specification. The purpose of this review is to collect updated information about lncRNA classification and new findings on their function derived from single-cell analysis. We also retained useful for all researchers to describe the methods available for single-cell analysis and the databases collecting single-cell and lncRNA data. Tables are included to schematize, describe, and compare exposed concepts.

New insights into transcriptional reprogramming during cellular stress

Cellular stress triggers reprogramming of transcription, which is required for the maintenance of homeostasis under adverse growth conditions. Stress-induced changes in transcription include induction of cyto-protective genes and repression of genes related to the regulation of the cell cycle, transcription and metabolism. Induction of transcription is mediated through the activation of stress-responsive transcription factors that facilitate the release of stalled RNA polymerase II and so allow for transcriptional elongation. Repression of transcription, in turn, involves components that retain RNA polymerase II in a paused state on gene promoters. Moreover, transcription during stress is regulated by a massive activation of enhancers and complex changes in chromatin organization. In this Review, we highlight the latest research regarding the molecular mechanisms of transcriptional reprogramming upon stress in the context of specific proteotoxic stress responses, including the heat-shock response, unfolded protein response, oxidative stress response and hypoxia response.

GATA2 and PU.1 Collaborate To Activate the Expression of the Mouse Ms4a2 Gene, Encoding FcεRIβ, through Distinct Mechanisms

GATA factors GATA1 and GATA2 and ETS factor PU.1 are known to function antagonistically during hematopoietic development. In mouse mast cells, however, these factors are coexpressed and activate the expression of the Ms4a2 gene encoding the β chain of the high-affinity IgE receptor (FcεRI). The present study showed that these factors cooperatively regulate Ms4a2 gene expression through distinct mechanisms. Although GATA2 and PU.1 contributed almost equally to Ms4a2 gene expression, gene ablation experiments revealed that simultaneous knockdown of both factors showed neither a synergistic nor an additive effect. A chromatin immunoprecipitation analysis showed that they shared DNA binding to the +10.4-kbp region downstream of the Ms4a2 gene with chromatin looping factor LDB1, whereas the proximal -60-bp region was exclusively bound by GATA2 in a mast cell-specific manner. Ablation of PU.1 significantly reduced the level of GATA2 binding to both the +10.4-kbp and -60-bp regions. Surprisingly, the deletion of the +10.4-kbp region by genome editing completely abolished the Ms4a2 gene expression as well as the cell surface expression of FcεRI. These results suggest that PU.1 and LDB1 play central roles in the formation of active chromatin structure whereas GATA2 directly activates the Ms4a2 promoter.

A novel enhancer RNA, Hmrhl, positively regulates its host gene, phkb, in chronic myelogenous leukemia

Noncoding RNAs are increasingly being accredited with key roles in gene regulation during development and disease. Here we report the discovery and characterization of a novel long noncoding RNA, Hmrhl, which shares synteny and partial sequence similarity with the mouse lncRNA, Mrhl. The human homolog, Hmrhl, transcribed from intron 14 of phkb gene, is 5.5 kb in size, expressed in all tissues examined and is associated with chromatin. Analysis of Hmrhl locus using ENCODE database revealed that it exhibits hallmarks of enhancers like the open chromatin configuration, binding of transcription factors, enhancer specific histone signature etc. in the K562 Chronic Myelogenous Leukemia (CML) cells. We compared the expression of Hmrhl in the normal lymphoblast cell line, GM12878, with that of K562 cells and lymphoma samples and show that it is highly upregulated in leukemia as well as several cases of lymphoma. Further, we validated the enhancer properties of Hmrhl locus in K562 cells with the help of ChIP-qPCR and Luciferase assay. Moreover, siRNA mediated down-regulation of Hmrhl in K562 cells leads to a concomitant down regulation of its parent gene, phkb, showing that Hmrhl functions as an enhancer RNA and positively regulates its host gene, phkb, in chronic myelogenous leukemia. This study is significant in view of the fact that a better understanding of mechanism of gene regulation under normal conditions and its perturbation in cancer could in turn help in its therapeutic intervention through molecular medicine/RNA based drug discovery.

Control of LDL Uptake in Human Cells by Targeting the LDLR Regulatory Long Non-coding RNA BM450697

Hypercholesterolemia is a condition that is characterized by very high levels of cholesterol in the blood and is a major correlating factor with heart disease. Indeed, high levels of the low-density lipoprotein (LDL) have been causally linked to the development of atherosclerotic cardiovascular disease (ASCVD). A method to specifically reduce cholesterol in the blood in a long-term, stable manner could prove therapeutically relevant. Cholesterol is removed from the blood by the LDL receptor (LDLR) in the liver. Others and we have discovered that a long non-coding RNA (lncRNA; BM450697) functions as an endogenous epigenetic regulator of LDLR and that the repression of this lncRNA by the action of small interfering RNAs (siRNAs) results in the activation of LDLR. We found here, through the interrogation of two siRNAs that can target this lncRNA, both in a transcriptional and post-transcriptional manner, that BM450697 functions as a local scaffold for modulating LDLR transcription. Moreover, we found that conjugation of α-N-acetylgalactosamine (GalNAc) with two lncRNA-directed siRNAs allows for direct liver cell targeting of this lncRNA and functional enhanced uptake of cholesterol. Collectively, these data suggest that targeting the BM450697 lncRNA regulator of LDLR may result in a more specific, long-term, targeted approach to regulating cholesterol in the blood.

The AP-1 transcriptional complex: Local switch or remote command?

The ubiquitous family of AP-1 dimeric transcription complexes is involved in virtually all cellular and physiological functions. It is paramount for cells to reprogram gene expression in response to cues of many sorts and is involved in many tumorigenic processes. How AP-1 controls gene transcription has largely remained elusive till recently. The advent of the "omics" technologies permitting genome-wide studies of transcription factors has however changed and improved our view of AP-1 mechanistical actions. If these studies confirm that AP-1 can sometimes act as a local transcriptional switch operating in the vicinity of transcription start sites (TSS), they strikingly indicate that AP-1 principally operates as a remote command binding to distal enhancers, placing chromatin architecture dynamics at the heart of its transcriptional actions. They also unveil novel constraints operating on AP-1, as well as novel mechanisms used to regulate gene expression via transcription-pioneering-, chromatin-remodeling- and chromatin accessibility maintenance effects.

Long Non-Coding RNAs and the Innate Immune Response

Innate immunity provides the initial defence against infection and it is now clear that long non-coding RNAs (lncRNAs) are important regulators of this response. Following activation of the innate response, we commonly see rapid induction of these lncRNAs and this is often mediated via the pro-inflammatory transcription factor, nuclear factor-κB (NF-κB). Knockdown studies have shown that lncRNAs tend to act in trans to regulate the expression of multiple inflammatory mediators and other responses. Mechanistically, many lncRNAs have demonstrated acting through heterogeneous nuclear ribonucleoproteins, complexes that are implicated chromatin re-modelling, transcription process and translation. In addition, these lncRNAs have also been shown to interact with multiple other proteins involved in the regulation of chromatin re-modelling, as well as those proteins involved in intracellular immune signalling, which include NF-κB. In this review, we will describe the evidence that supports this emerging role of lncRNA in the innate immune response.

Regulatory Non-coding RNAs Network in Non-alcoholic Fatty Liver Disease

Non-alcoholic fatty liver disease (NAFLD) spectrum comprises simple steatosis and non-alcoholic steatohepatitis (NASH) that can lead to fibrosis and cirrhosis. The patients usually have no history of excessive alcohol consumption and other etiologies that can cause fatty liver. Understanding of the pathophysiology of NAFLD has revealed that non-coding RNAs (ncRNAs) play significant roles in modulating the disease susceptibility, pathogenesis and progression. Currently, the ncRNAs are grouped according to their sizes and their regulatory or housekeeping functions. Each of these ncRNAs has a wide range of involvement in the regulation of the genes and biological pathways. Here, we briefly review the current literature the regulatory ncRNAs in NAFLD pathogenesis and progression, mainly the microRNAs, long non-coding RNAs and circular RNAs. We also discuss the co-regulatory functions and interactions between these ncRNAs in modulating the disease pathogenesis. Elucidation of ncRNAs in NAFLD may facilitate the identification of early diagnostic biomarkers and development of therapeutic strategies for NAFLD.

AP001056.1, A Prognosis-Related Enhancer RNA in Squamous Cell Carcinoma of the Head and Neck

A growing number of long non-coding RNAs (lncRNAs) have been linked to squamous cell carcinoma of the head and neck (SCCHN). A subclass of lncRNAs, termed enhancer RNAs (eRNAs), are derived from enhancer regions and could contribute to enhancer function. In this study, we developed an integrated data analysis approach to identify key eRNAs in SCCHN. Tissue-specific enhancer-derived RNAs and their regulated genes previously predicted using the computational pipeline PreSTIGE, were considered as putative eRNA-target pairs. The interactive web servers, TANRIC (the Atlas of Noncoding RNAs in Cancer) and cBioPortal, were used to explore the RNA levels and clinical data from the Cancer Genome Atlas (TCGA) project. Requiring that key eRNAs should show significant associations with overall survival (Kaplan⁻Meier log-rank test, p < 0.05) and the predicted target (correlation coefficient r > 0.4, p < 0.001), we identified five key eRNA candidates. The most significant survival-associated eRNA was AP001056.1 with ICOSLG encoding an immune checkpoint protein as its regulated target. Another 1640 genes also showed significant correlation with AP001056.1 (r > 0.4, p < 0.001), with the "immune system process" being the most significantly enriched biological process (adjusted p < 0.001). Our results suggest that AP001056.1 is a key immune-related eRNA in SCCHN with a positive impact on clinical outcome.

Characterization and identification of long non-coding RNAs based on feature relationship

Motivation

The significance of long non-coding RNAs (lncRNAs) in many biological processes and diseases has gained intense interests over the past several years. However, computational identification of lncRNAs in a wide range of species remains challenging; it requires prior knowledge of well-established sequences and annotations or species-specific training data, but the reality is that only a limited number of species have high-quality sequences and annotations.

Results

Here we first characterize lncRNAs in contrast to protein-coding RNAs based on feature relationship and find that the feature relationship between open reading frame length and guanine-cytosine (GC) content presents universally substantial divergence in lncRNAs and protein-coding RNAs, as observed in a broad variety of species. Based on the feature relationship, accordingly, we further present LGC, a novel algorithm for identifying lncRNAs that is able to accurately distinguish lncRNAs from protein-coding RNAs in a cross-species manner without any prior knowledge. As validated on large-scale empirical datasets, comparative results show that LGC outperforms existing algorithms by achieving higher accuracy, well-balanced sensitivity and specificity, and is robustly effective (>90% accuracy) in discriminating lncRNAs from protein-coding RNAs across diverse species that range from plants to mammals. To our knowledge, this study, for the first time, differentially characterizes lncRNAs and protein-coding RNAs based on feature relationship, which is further applied in computational identification of lncRNAs. Taken together, our study represents a significant advance in characterization and identification of lncRNAs and LGC thus bears broad potential utility for computational analysis of lncRNAs in a wide range of species.

Availability and implementation

LGC web server is publicly available at http://bigd.big.ac.cn/lgc/calculator. The scripts and data can be downloaded at http://bigd.big.ac.cn/biocode/tools/BT000004.

Supplementary information

Supplementary data are available at Bioinformatics online.

Enhancer-Promoter Interaction of SELF PRUNING 5G Shapes Photoperiod Adaptation

Tomato (Solanum lycopersicum) is a major vegetable fruit grown and consumed worldwide. Modern cultivated tomatoes are derived from their wild relative, Solanum pimpinellifolium, a short-day plant that originated from the Andean region of South America. The molecular underpinnings of the regional adaptation and expansion of domesticated tomato remain largely unclear. In this study, we examined flowering time in wild and cultivated tomatoes under both long-day and short-day conditions. Using quantitative trait locus mapping in a recombinant inbred line population, we identified SELF PRUNING 5G (SP5G) as a major locus influencing daylength adaptation in tomato. Genetic diversity analysis revealed that the genomic region harboring SP5G shows signatures of a domestication sweep. We found that a 52-bp sequence within the 3' untranslated region of SP5G is essential for the enhanced expression of this gene, leading to delayed flowering time in tomatoes through a promoter-enhancer interaction that occurs only under long-day conditions. We further demonstrate that the absence of the 52-bp sequence attenuates the promoter-enhancer interaction and reduces SP5G expression in cultivated tomatoes, making their flowering time insensitive to daylength. Our findings demonstrate that cis-regulatory variation at the enhancer region of the SP5G 3' untranslated region confers reduced photoperiodic response in cultivated tomatoes, uncovering a regulatory mechanism that could potentially be used to manipulate flowering time in tomato through novel biotechnological approaches.

The enhancer RNA lnc-SLC4A1-1 epigenetically regulates unexplained recurrent pregnancy loss (URPL) by activating CXCL8 and NF-kB pathway

Background

Enhancer RNAs (eRNAs) are a group of lncRNAs transcribed from enhancers, whose regulatory effects on gene expression are an emerging area of interest. However, the role of eRNAs in regulating trophoblast cells and unexplained recurrent pregnancy loss (URPL) remains elusive.

Methods

We profiled eRNAs in villi from URPL patients and matched controls by RNA-seq. Functions of URPL-related eRNAs were further investigated in vitro.

Results

We identified lnc-SLC4A1-1, which was transcribed from an active enhancer marked with H3K27ac and H3K4me1 and so-called eRNA, highly expressed in URPL patients. Gain-of-function experiments indicated that lnc-SLC4A1-1 facilitated trophoblast cell migration and apoptosis. Mechanistically, as an eRNA, lnc-SLC4A1-1 was retained in the nuclei and recruited transcription factor NF-κB to bind to CXCL8, resulting in increased H3K27ac in the CXCL8 promoter and subsequent elevation of CXCL8 expression. Activation of CXCL8 exacerbated inflammatory reactions in trophoblast cells by inducing TNF-α and IL-1β, which could be blocked by an antagonist of lnc-SLC4A1-1.

Interpretation

These findings indicate that an eRNA, lnc-SLC4A1-1, alters trophoblast function via activation of immune responses and by regulating the NF-κB/CXCL8 axis. Our study provides new insights in understanding lncRNA/eRNA function in pathological pregnancy, potentially informing on therapeutic strategies for URPL.

Fund

National Natural Science Foundation of China, Natural Science Foundation of Jiangsu Province, National Key Research and Development Program, the Priority Academic Program for the Development of Jiangsu Higher Education Institutions.

Cross-species inference of long non-coding RNAs greatly expands the ruminant transcriptome

Background

mRNA-like long non-coding RNAs (lncRNAs) are a significant component of mammalian transcriptomes, although most are expressed only at low levels, with high tissue-specificity and/or at specific developmental stages. Thus, in many cases lncRNA detection by RNA-sequencing (RNA-seq) is compromised by stochastic sampling. To account for this and create a catalogue of ruminant lncRNAs, we compared de novo assembled lncRNAs derived from large RNA-seq datasets in transcriptional atlas projects for sheep and goats with previous lncRNAs assembled in cattle and human. We then combined the novel lncRNAs with the sheep transcriptional atlas to identify co-regulated sets of protein-coding and non-coding loci.

Results

Few lncRNAs could be reproducibly assembled from a single dataset, even with deep sequencing of the same tissues from multiple animals. Furthermore, there was little sequence overlap between lncRNAs that were assembled from pooled RNA-seq data. We combined positional conservation (synteny) with cross-species mapping of candidate lncRNAs to identify a consensus set of ruminant lncRNAs and then used the RNA-seq data to demonstrate detectable and reproducible expression in each species. In sheep, 20 to 30% of lncRNAs were located close to protein-coding genes with which they are strongly co-expressed, which is consistent with the evolutionary origin of some ncRNAs in enhancer sequences. Nevertheless, most of the lncRNAs are not co-expressed with neighbouring protein-coding genes.

Conclusions

Alongside substantially expanding the ruminant lncRNA repertoire, the outcomes of our analysis demonstrate that stochastic sampling can be partly overcome by combining RNA-seq datasets from related species. This has practical implications for the future discovery of lncRNAs in other species.

Electronic supplementary material

The online version of this article (10.1186/s12711-018-0391-0) contains supplementary material, which is available to authorized users.

A family of long intergenic non-coding RNA genes in human chromosomal region 22q11.2 carry a DNA translocation breakpoint/AT-rich sequence

FAM230C, a long intergenic non-coding RNA (lincRNA) gene in human chromosome 13 (chr13) is a member of lincRNA genes termed family with sequence similarity 230. An analysis using bioinformatics search tools and alignment programs was undertaken to determine properties of FAM230C and its related genes. Results reveal that the DNA translocation element, the Translocation Breakpoint Type A (TBTA) sequence, which consists of satellite DNA, Alu elements, and AT-rich sequences is embedded in the FAM230C gene. Eight lincRNA genes related to FAM230C also carry the TBTA sequences. These genes were formed from a large segment of the 3’ half of the FAM230C sequence duplicated in chr22, and are specifically in regions of low copy repeats (LCR22)s, in or close to the 22q.11.2 region. 22q11.2 is a chromosomal segment that undergoes a high rate of DNA translocation and is prone to genetic deletions. FAM230C-related genes present in other chromosomes do not carry the TBTA motif and were formed from the 5’ half region of the FAM230C sequence. These findings identify a high specificity in lincRNA gene formation by gene sequence duplication in different chromosomes.