LongTarget: a tool to predict lncRNA DNA-binding motifs and binding sites via Hoogsteen base-pairing analysis

Human-specific protein-coding and lncRNA genes cast sex-biased genes in the brain and their relationships with brain diseases

BACKGROUND

Gene expression shows sex bias in the brain as it does in other organs. Since female and male humans exhibit noticeable differences in emotions, logical thinking, movement, spatial orientation, and even the incidence of neurological disorders, sex biases in the brain are especially interesting, but how they are determined, whether they are conserved or lineage specific, and what the consequences of the biases are, remain poorly explored and understood.

METHODS

Based on RNA-seq datasets from 16  and 14 brain regions in humans and macaques across developmental periods and from patients with brain diseases, we used linear mixed models (LMMs) to differentiate variations in gene expression caused by factors of interest and confounding factors and identify four types of sex-biased genes. Effect size and confidence in each effect were measured upon the local false sign rate (LFSR). We utilized the biomaRt R package to acquire orthologous genes in humans and macaques from the BioMart Ensembl website. Transcriptional regulation of sex-biased genes by sex hormones and lncRNAs were analyzed using the CellOracle, GENIE3, and Longtarget programs. Sex-biased genes' functions were revealed by gene set enrichment analysis using multiple methods.

RESULTS

Lineage-specific sex-biased genes greatly determine the distinct sex biases in human and macaque brains. In humans, those encoding proteins contribute directly to immune-related functions, and those encoding lncRNAs intensively regulate the expression of other sex-biased genes, especially genes with immune-related functions. The identified sex-specific differentially expressed genes (ssDEGs) upon gene expression in disease and normal samples also indicate that protein-coding ssDEGs are conserved in humans and macaques but that lncRNA ssDEGs are not conserved. The results answer the above questions, reveal an intrinsic relationship between sex biases in the brain and sex-biased susceptibility to brain diseases, and will help researchers investigate human- and sex-specific ncRNA targets for brain diseases.

CONCLUSIONS

Human-specific genes greatly cast sex-biased genes in the brain and their relationships with brain diseases, with protein-coding genes contributing to immune response related functions and lncRNA genes critically regulating sex-biased genes. The high proportions of lineage-specific lncRNAs in mammalian genomes indicate that sex biases may have evolved rapidly in not only the brain but also other organs.

Alternative splicing of CARM1 regulated by LincGET-guided paraspeckles biases the first cell fate in mammalian early embryos

The heterogeneity of CARM1 controls first cell fate bias during early mouse development. However, how this heterogeneity is established is unknown. Here, we show that Carm1 mRNA is of a variety of specific exon-skipping splicing (ESS) isoforms in mouse two-cell to four-cell embryos that contribute to CARM1 heterogeneity. Disruption of paraspeckles promotes the ESS of Carm1 precursor mRNAs (pre-mRNAs). LincGET, but not Neat1, is required for paraspeckle assembly and inhibits the ESS of Carm1 pre-mRNAs in mouse two-cell to four-cell embryos. We further find that LincGET recruits paraspeckles to the Carm1 gene locus through HNRNPU. Interestingly, PCBP1 binds the Carm1 pre-mRNAs and promotes its ESS in the absence of LincGET. Finally, we find that the ESS seen in mouse two-cell to four-cell embryos decreases CARM1 protein levels and leads to trophectoderm fate bias. Our findings demonstrate that alternative splicing of CARM1 has an important role in first cell fate determination.

Targeting and engineering long non-coding RNAs for cancer therapy

RNA therapeutics (RNATx) aim to treat diseases, including cancer, by targeting or employing RNA molecules for therapeutic purposes. Amongst the most promising targets are long non-coding RNAs (lncRNAs), which regulate oncogenic molecular networks in a cell type-restricted manner. lncRNAs are distinct from protein-coding genes in important ways that increase their therapeutic potential yet also present hurdles to conventional clinical development. Advances in genome editing, oligonucleotide chemistry, multi-omics and RNA engineering are paving the way for efficient and cost-effective lncRNA-focused drug discovery pipelines. In this Review, we present the emerging field of lncRNA therapeutics for oncology, with emphasis on the unique strengths and challenges of lncRNAs within the broader RNATx framework. We outline the necessary steps for lncRNA therapeutics to deliver effective, durable, tolerable and personalized treatments for cancer.

RNA-DNA triplexes: molecular mechanisms and functional relevance

Interactions of RNA with DNA are principles of gene expression control that have recently gained considerable attention. Among RNA-DNA interactions are R-loops and RNA-DNA hybrid G-quadruplexes, as well as RNA-DNA triplexes. It is proposed that RNA-DNA triplexes guide RNA-associated regulatory proteins to specific genomic locations, influencing transcription and epigenetic decision making. Although triplex formation initially was considered solely an in vitro event, recent progress in computational, biochemical, and biophysical methods support in vivo functionality with relevance for gene expression control. Here, we review the central methodology and biology of triplexes, outline paradigms required for triplex function, and provide examples of physiologically important triplex-forming long non-coding RNAs.

PVT1 interacts with polycomb repressive complex 2 to suppress genomic regions with pro-apoptotic and tumour suppressor functions in multiple myeloma

Multiple myeloma is a heterogeneous hematological disease that originates from the bone marrow and is characterized by the monoclonal expansion of malignant plasma cells. Despite novel therapies, multiple myeloma remains clinically challenging. A common feature among patients with poor prognosis is the increased activity of the epigenetic silencer EZH2, which is the catalytic subunit of the PRC2. Interestingly, the recruitment of PRC2 lacks sequence specificity and, to date, the molecular mechanisms that define which genomic locations are destined for PRC2-mediated silencing remain unknown. The presence of a long non-coding RNA (lncRNA)-binding pocket on EZH2 suggests that lncRNA could potentially mediate PRC2 recruitment to specific genomic regions. Here, we coupled RNA immunoprecipitation sequencing, RNA-sequencing and chromatin immunoprecipitation-sequencing analysis of human multiple myeloma primary cells and cell lines to identify potential lncRNA partners to EZH2. We found that the lncRNA plasmacytoma variant translocation 1 (PVT1) directly interacts with EZH2 and is overexpressed in patients with a poor prognosis. Moreover, genes predicted to be targets of PVT1 exhibited H3K27me3 enrichment and were associated with pro-apoptotic and tumor suppressor functions. In fact, PVT1 inhibition independently promotes the expression of the PRC2 target genes ZBTB7C, RNF144A and CCDC136. Altogether, our work suggests that PVT1 is an interacting partner in PRC2-mediated silencing of tumor suppressor and pro-apoptotic genes in multiple myeloma, making it a highly interesting potential therapeutic target.

LincR-PPP2R5C Promotes Th2 Cell Differentiation Through PPP2R5C/PP2A by Forming an RNA-DNA Triplex in Allergic Asthma

PURPOSE

The roles and mechanisms of long noncoding RNAs (lncRNAs) in T helper 2 (Th2) differentiation from allergic asthma are poorly understood. We aimed to explore a novel lncRNA, LincR-protein phosphatase 2 regulatory subunit B' gamma (PPP2R5C), in Th2 differentiation in a mouse model of asthma.

METHODS

LincR-PPP2R5C from RNA-seq data of CD4+ T cells of asthma-like mice were validated and confirmed by quantitative reverse transcription polymerase chain reaction, northern blotting, nuclear and cytoplasmic separation, and fluorescence in situ hybridization (FISH). Lentiviruses encoding LincR-PPP2R5C or shRNA were used to overexpress or silence LincR-PPP2R5C in CD4+ T cells. The interactions between LincR-PPP2R5C and PPP2R5C were explored with western blotting, chromatin isolation by RNA purification assay, and fluorescence resonance energy transfer. An ovalbumin-induced acute asthma model in knockout (KO) mice (LincR-PPP2R5C KO, CD4 conditional LincR-PPP2R5C KO) was established to explore the roles of LincR-PPP2R5C in Th2 differentiation.

RESULTS

LncR-PPP2R5C was significantly higher in CD4+ T cells from asthmatic mice ex vivo and Th2 cells in vitro. The lentivirus encoding LincR-PPP2R5C suppressed Th1 differentiation; in contrast, the short hairpin RNA (shRNA) lentivirus decreased LincR-PPP2R5C and Th2 differentiation. Mechanistically, LincR-PPP2R5C deficiency suppressed the phosphatase activity of the protein phosphatase 2A (PP2A) holocomplex, resulting in a decline in Th2 differentiation. The formation of an RNA-DNA triplex between LincR-PPP2R5C and the PPP2R5C promoter enhanced PPP2R5C expression and activated PP2A. LincR-PPP2R5C KO and CD4 conditional KO decreased Th2 differentiation, airway hyperresponsiveness and inflammatory responses.

CONCLUSIONS

LincR-PPP2R5C regulated PPP2R5C expression and PP2A activity by forming an RNA-DNA triplex with the PPP2R5C promoter, leading to Th2 polarization in a mouse model of acute asthma. Our data presented the first definitive evidence of lncRNAs in the regulation of Th2 cells in asthma.

A PPARγ/long noncoding RNA axis regulates adipose thermoneutral remodeling in mice

Interplay between energy-storing white adipose cells and thermogenic beige adipocytes contributes to obesity and insulin resistance. Irrespective of specialized niche, adipocytes require the activity of the nuclear receptor PPARγ for proper function. Exposure to cold or adrenergic signaling enriches thermogenic cells though multiple pathways that act synergistically with PPARγ; however, the molecular mechanisms by which PPARγ licenses white adipose tissue to preferentially adopt a thermogenic or white adipose fate in response to dietary cues or thermoneutral conditions are not fully elucidated. Here, we show that a PPARγ/long noncoding RNA (lncRNA) axis integrates canonical and noncanonical thermogenesis to restrain white adipose tissue heat dissipation during thermoneutrality and diet-induced obesity. Pharmacologic inhibition or genetic deletion of the lncRNA Lexis enhances uncoupling protein 1-dependent (UCP1-dependent) and -independent thermogenesis. Adipose-specific deletion of Lexis counteracted diet-induced obesity, improved insulin sensitivity, and enhanced energy expenditure. Single-nuclei transcriptomics revealed that Lexis regulates a distinct population of thermogenic adipocytes. We systematically map Lexis motif preferences and show that it regulates the thermogenic program through the activity of the metabolic GWAS gene and WNT modulator TCF7L2. Collectively, our studies uncover a new mode of crosstalk between PPARγ and WNT that preserves white adipose tissue plasticity.

Pan-glioma analyses reveal species- and tumor-specific regulation of neuron-glioma synapse genes by lncRNAs

Gliomas are highly heterogeneous and aggressive. Malignant cells in gliomas can contact normal neurons through a synapse-like structure (called neuron-to-glioma synapse, NGS) to promote their proliferation, but it is unclear whether NGS gene expression and regulation show species- and tumor-specificity. This question is important in that many anti-cancer drugs are developed upon mouse models. To address this question, we conducted a pan-glioma analysis using nine scRNA-seq datasets from humans and mice. We also experimentally validated the key element of our methods and verified a key result using TCGA datasets of the same glioma types. Our analyses revealed that NGS gene expression and regulation by lncRNAs are highly species- and tumor-specific. Importantly, simian-specific lncRNAs are more involved in NGS gene regulation than lncRNAs conserved in mammals, and transgenic mouse gliomas have little in common with PDX mouse models and human gliomas in terms of NGS gene regulation. The analyses suggest that simian-specific lncRNAs are a new and rich class of potential targets for tumor-specific glioma treatment, and provide pertinent data for further experimentally and clinically exmining the targets.

TFOFinder: Python program for identifying purine-only double-stranded stretches in the predicted secondary structure(s) of RNA targets

Nucleic acid probes are valuable tools in biology and chemistry and are indispensable for PCR amplification of DNA, RNA quantification and visualization, and downregulation of gene expression. Recently, triplex-forming oligonucleotides (TFO) have received increased attention due to their improved selectivity and sensitivity in recognizing purine-rich double-stranded RNA regions at physiological pH by incorporating backbone and base modifications. For example, triplex-forming peptide nucleic acid (PNA) oligomers have been used for imaging a structured RNA in cells and inhibiting influenza A replication. Although a handful of programs are available to identify triplex target sites (TTS) in DNA, none are available that find such regions in structured RNAs. Here, we describe TFOFinder, a Python program that facilitates the identification of intramolecular purine-only RNA duplexes that are amenable to forming parallel triple helices (pyrimidine/purine/pyrimidine) and the design of the corresponding TFO(s). We performed genome- and transcriptome-wide analyses of TTS in Drosophila melanogaster and found that only 0.3% (123) of total unique transcripts (35,642) show the potential of forming 12-purine long triplex forming sites that contain at least one guanine. Using minimization algorithms, we predicted the secondary structure(s) of these transcripts, and using TFOFinder, we found that 97 (79%) of the identified 123 transcripts are predicted to fold to form at least one TTS for parallel triple helix formation. The number of transcripts with potential purine TTS increases when the strict search conditions are relaxed by decreasing the length of the probe or by allowing up to two pyrimidine inversions or 1-nucleotide bulge in the target site. These results are encouraging for the use of modified triplex forming probes for live imaging of endogenous structured RNA targets, such as pre-miRNAs, and inhibition of target-specific translation and viral replication.

Long noncoding RNA AI662270 promotes kidney fibrosis through enhancing METTL3-mediated m6 A modification of CTGF mRNA

The sustained release of profibrotic cytokines, mainly transforming growth factor-β (TGF-β), leads to the occurrence of kidney fibrosis and chronic kidney disease (CKD). Connective tissue growth factor (CTGF) appears to be an alternative target to TGF-β for antifibrotic therapy in CKD. In this study, we found that long noncoding RNA AI662270 was significantly increased in various renal fibrosis models. In vivo, ectopic expression of AI662270 alone was sufficient to activate interstitial fibroblasts and drive kidney fibrosis, whereas inhibition of AI662270 blocked the activation of interstitial fibroblasts and ameliorated kidney fibrosis in various murine models. Mechanistic studies revealed that overexpression of AI662270 significantly increased CTGF product, which was required for the role of AI662270 in driving kidney fibrosis. Furthermore, AI662270 binds to the CTGF promoter and directly interacts with METTL3, the methyltransferase of RNA N6 -methyladenosine (m6 A) modification. Functionally, AI662270-mediated recruitment of METTL3 increased the m6 A methylation of CTGF mRNA and consequently enhanced CTGF mRNA stability. In conclusion, our results support that AI662270 promotes CTGF expression at the posttranscriptional stage by recruiting METTL3 to the CTGF promoter and depositing m6 A modifications on the nascent mRNA, thereby, uncovering a novel regulatory mechanism of CTGF in the pathogenesis of kidney fibrosis.

Small extracellular vesicles-transported lncRNA TDRKH-AS1 derived from AOPPs-treated trophoblasts initiates endothelial cells pyroptosis through PDIA4/DDIT4 axis in preeclampsia

BACKGROUND

Substantial studies have demonstrated that oxidative stress placenta and endothelial injury are considered to inextricably critical events in the pathogenesis of preeclampsia (PE). Systemic inflammatory response and endothelial dysfunction are induced by the circulating factors released from oxidative stress placentae. As a novel biomarker of oxidative stress, advanced oxidation protein products (AOPPs) levels are strongly correlated with PE characteristics. Nevertheless, the molecular mechanism underlying the effect of factors is still largely unknown.

METHODS

With the exponential knowledge on the importance of placenta-derived extracellular vesicles (pEVs), we carried out lncRNA transcriptome profiling on small EVs (sEVs) secreted from AOPPs-treated trophoblast cells and identified upregulated lncRNA TDRKH-AS1 as a potentially causative factor for PE. We isolated and characterized sEVs from plasma and trophoblast cells by transmission electron microscopy (TEM), nanoparticle tracking analysis (NTA) and western blotting. The expression and correlation of lncRNA TDRKH-AS1 were evaluated using qRT-PCR in plasmatic sEVs and placentae from patients. Pregnant mice injected with TDRKH-AS1-riched trophoblast sEVs was performed to detect the TDRKH-AS1 function in vivo. To investigate the potential effect of sEVs-derived TDRKH-AS1 on endothelial function in vitro, transcriptome sequencing, scanning electron Microscopy (SEM), immunofluorescence, ELISA and western blotting were conducted in HUVECs. RNA pulldown, mass spectrometry, RNA immunoprecipitation (RIP), chromatin isolation by RNA purification (ChIRP) and coimmunoprecipitation (Co-IP) were used to reveal the latent mechanism of TDRKH-AS1 on endothelial injury.

RESULTS

The expression level of TDRKH-AS1 was significantly increased in plasmatic sEVs and placentae from patients, and elevated TDRKH-AS1 in plasmatic sEVs was positively correlated with clinical severity of the patients. Moreover, pregnant mice injected with TDRKH-AS1-riched trophoblast sEVs exhibited a hallmark feature of PE with increased blood pressure and systemic inflammatory responses. Pyroptosis, an inflammatory form of programmed cell death, is involved in the development of PE. Indeed, our in vitro study indicated that sEVs-derived TDRKH-AS1 secreted from AOPPs-induced trophoblast elevated DDIT4 expression levels to trigger inflammatory response of pyroptosis in endothelial cells through interacting with PDIA4.

CONCLUSIONS

Herein, results in the present study supported that TDRKH-AS1 in sEVs isolated from oxidative stress trophoblast may be implicated in the pathogenesis of PE via inducing pyroptosis and aggravating endothelial dysfunction.

Interactome of Long Non-Coding RNAs: Transcriptomic Expression Patterns and Shaping Cancer Cell Phenotypes

RNA biology has gained extensive recognition in the last two decades due to the identification of novel transcriptomic elements and molecular functions. Cancer arises, in part, due to the accumulation of mutations that greatly contribute to genomic instability. However, the identification of differential gene expression patterns of wild-type loci has exceeded the boundaries of mutational study and has significantly contributed to the identification of molecular mechanisms that drive carcinogenic transformation. Non-coding RNA molecules have provided a novel avenue of exploration, providing additional routes for evaluating genomic and epigenomic regulation. Of particular focus, long non-coding RNA molecule expression has been demonstrated to govern and direct cellular activity, thus evidencing a correlation between aberrant long non-coding RNA expression and the pathological transformation of cells. lncRNA classification, structure, function, and therapeutic utilization have expanded cancer studies and molecular targeting, and understanding the lncRNA interactome aids in defining the unique transcriptomic signatures of cancer cell phenotypes.

TRIPBASE: a database for identifying the human genomic DNA and lncRNA triplexes

Long-non-coding RNAs (lncRNAs) are defined as RNA sequences which are >200 nt with no coding capacity. These lncRNAs participate in various biological mechanisms, and are widely abundant in a diversity of species. There is well-documented evidence that lncRNAs can interact with genomic DNAs by forming triple helices (triplexes). Previously, several computational methods have been designed based on the Hoogsteen base-pair rule to find theoretical RNA-DNA:DNA triplexes. While powerful, these methods suffer from a high false-positive rate between the predicted triplexes and the biological experiments. To address this issue, we first collected the experimental data of genomic RNA-DNA triplexes from antisense oligonucleotide (ASO)-mediated capture assays and used Triplexator, the most widely used tool for lncRNA-DNA interaction, to reveal the intrinsic information on true triplex binding potential. Based on the analysis, we proposed six computational attributes as filters to improve the in-silico triplex prediction by removing most false positives. Further, we have built a new database, TRIPBASE, as the first comprehensive collection of genome-wide triplex predictions of human lncRNAs. In TRIPBASE, the user interface allows scientists to apply customized filtering criteria to access the potential triplexes of human lncRNAs in the cis-regulatory regions of the human genome. TRIPBASE can be accessed at https://tripbase.iis.sinica.edu.tw/.

3plex enables deep computational investigation of triplex forming lncRNAs

Long non-coding RNAs (lncRNAs) regulate gene expression through different molecular mechanisms, including DNA binding via the formation of RNA:DNA:DNA triple helices (TPXs). Despite the increasing amount of experimental evidence, TPXs investigation remains challenging. Here we present 3plex, a software able to predict TPX interactions in silico. Given an RNA sequence and a set of DNA sequences, 3plex integrates 1) Hoogsteen pairing rules that describe the biochemical interactions between RNA and DNA nucleotides, 2) RNA secondary structure prediction and 3) determination of the TPX thermal stability derived from a collection of TPX experimental evidences. We systematically collected and uniformly re-analysed published experimental lncRNA binding sites on human and mouse genomes. We used these data to evaluate 3plex performance and showed that its specific features allow a reliable identification of TPX interactions. We compared 3plex with the other available software and obtained comparable or even better accuracy at a fraction of the computation time. Interestingly, by inspecting collected data with 3plex we found that TPXs tend to be shorter and more degenerated than previously expected and that the majority of analysed lncRNAs can directly bind to the genome by TPX formation. Those results suggest that an important fraction of lncRNAs can exert its biological function through this mechanism. The software is available at https://github.com/molinerisLab/3plex.

A hybrid approach to assess the structural impact of long noncoding RNA mutations uncovers key NEAT1 interactions in colorectal cancer

Long noncoding RNAs (lncRNAs) are emerging players in cancer and they entail potential as prognostic biomarkers or therapeutic targets. Earlier studies have identified somatic mutations in lncRNAs that are associated with tumor relapse after therapy, but the underlying mechanisms behind these associations remain unknown. Given the relevance of secondary structure for the function of some lncRNAs, some of these mutations may have a functional impact through structural disturbance. Here, we examined the potential structural and functional impact of a novel A > G point mutation in NEAT1 that has been recurrently observed in tumors of colorectal cancer patients experiencing relapse after treatment. Here, we used the nextPARS structural probing approach to provide first empirical evidence that this mutation alters NEAT1 structure. We further evaluated the potential effects of this structural alteration using computational tools and found that this mutation likely alters the binding propensities of several NEAT1-interacting miRNAs. Differential expression analysis on these miRNA networks shows upregulation of Vimentin, consistent with previous findings. We propose a hybrid pipeline that can be used to explore the potential functional effects of lncRNA somatic mutations.

Aid or Antagonize: Nuclear Long Noncoding RNAs Regulate Host Responses and Outcomes of Viral Infections

Long noncoding RNAs (lncRNAs) are transcripts measuring >200 bp in length and devoid of protein-coding potential. LncRNAs exceed the number of protein-coding mRNAs and regulate cellular, developmental, and immune pathways through diverse molecular mechanisms. In recent years, lncRNAs have emerged as epigenetic regulators with prominent roles in health and disease. Many lncRNAs, either host or virus-encoded, have been implicated in critical cellular defense processes, such as cytokine and antiviral gene expression, the regulation of cell signaling pathways, and the activation of transcription factors. In addition, cellular and viral lncRNAs regulate virus gene expression. Viral infections and associated immune responses alter the expression of host lncRNAs regulating immune responses, host metabolism, and viral replication. The influence of lncRNAs on the pathogenesis and outcomes of viral infections is being widely explored because virus-induced lncRNAs can serve as diagnostic and therapeutic targets. Future studies should focus on thoroughly characterizing lncRNA expressions in virus-infected primary cells, investigating their role in disease prognosis, and developing biologically relevant animal or organoid models to determine their suitability for specific therapeutic targeting. Many cellular and viral lncRNAs localize in the nucleus and epigenetically modulate viral transcription, latency, and host responses to infection. In this review, we provide an overview of the role of nuclear lncRNAs in the pathogenesis and outcomes of viral infections, such as the Influenza A virus, Sendai Virus, Respiratory Syncytial Virus, Hepatitis C virus, Human Immunodeficiency Virus, and Herpes Simplex Virus. We also address significant advances and barriers in characterizing lncRNA function and explore the potential of lncRNAs as therapeutic targets.

Dynamic chromatin accessibility tuning by the long noncoding RNA ELDR accelerates chondrocyte senescence and osteoarthritis

Epigenetic reprogramming plays a critical role in chondrocyte senescence during osteoarthritis (OA) pathology, but the underlying molecular mechanisms remain to be elucidated. Here, using large-scale individual datasets and genetically engineered (Col2a1-CreER^T2;Eldr^flox/flox and Col2a1-CreER^T2;ROSA26-LSL-Eldr^+/+ knockin) mouse models, we show that a novel transcript of long noncoding RNA ELDR is essential for the development of chondrocyte senescence. ELDR is highly expressed in chondrocytes and cartilage tissues of OA. Mechanistically, exon 4 of ELDR physically mediates a complex consisting of hnRNPL and KAT6A to regulate histone modifications of the promoter region of IHH, thereby activating hedgehog signaling and promoting chondrocyte senescence. Therapeutically, GapmeR-mediated silencing of ELDR in the OA model substantially attenuates chondrocyte senescence and cartilage degradation. Clinically, ELDR knockdown in cartilage explants from OA-affected individuals decreased the expression of senescence markers and catabolic mediators. Taken together, these findings uncover an lncRNA-dependent epigenetic driver in chondrocyte senescence, highlighting that ELDR could be a promising therapeutic avenue for OA.

Recent Development in Biomedical Applications of Oligonucleotides with Triplex-Forming Ability

A DNA structure, known as triple-stranded DNA, is made up of three oligonucleotide chains that wind around one another to form a triple helix (TFO). Hoogsteen base pairing describes how triple-stranded DNA may be built at certain conditions by the attachment of the third strand to an RNA, PNA, or DNA, which might all be employed as oligonucleotide chains. In each of these situations, the oligonucleotides can be employed as an anchor, in conjunction with a specific bioactive chemical, or as a messenger that enables switching between transcription and replication through the triplex-forming zone. These data are also considered since various illnesses have been linked to the expansion of triplex-prone sequences. In light of metabolic acidosis and associated symptoms, some consideration is given to the impact of several low-molecular-weight compounds, including pH on triplex production in vivo. The review is focused on the development of biomedical oligonucleotides with triplexes.

Human Macrophage Long Intergenic Noncoding RNA, SIMALR, Suppresses Inflammatory Macrophage Apoptosis via NTN1 (Netrin-1)

BACKGROUND

Long noncoding RNAs (lncRNAs) have emerged as novel regulators of macrophage biology and inflammatory cardiovascular diseases. However, studies focused on lncRNAs in human macrophage subtypes, particularly human lncRNAs that are not conserved in rodents, are limited.

METHODS

Through RNA-sequencing of human monocyte-derived macrophages, we identified suppressor of inflammatory macrophage apoptosis lncRNA (SIMALR). Lipopolysaccharide/IFNγ (interferon γ) stimulated human macrophages were treated with SIMALR antisense oligonucleotides and subjected to RNA-sequencing to investigate the function of SIMALR. Western blots, luciferase assay, and RNA immunoprecipitation were performed to validate function and potential mechanism of SIMALR. RNAscope was performed to identify SIMALR expression in human carotid atherosclerotic plaques.

RESULTS

RNA-sequencing of human monocyte-derived macrophages identified SIMALR, a human macrophage-specific long intergenic noncoding RNA that is highly induced in lipopolysaccharide/IFNγ-stimulated macrophages. SIMALR knockdown in lipopolysaccharide/IFNγ stimulated THP1 human macrophages induced apoptosis of inflammatory macrophages, as shown by increased protein expression of cleaved PARP (poly[ADP-ribose] polymerase), caspase 9, caspase 3, and Annexin V+. RNA-sequencing of control versus SIMALR knockdown in lipopolysaccharide/IFNγ-stimulated macrophages showed Netrin-1 (NTN1) to be significantly decreased upon SIMALR knockdown. We confirmed that NTN1 knockdown in lipopolysaccharide/IFNγ-stimulated macrophages induced apoptosis. The SIMALR knockdown-induced apoptotic phenotype was rescued by adding recombinant NTN1. NTN1 promoter-luciferase reporter activity was increased in HEK293T (human embryonic kidney 293) cells treated with lentiviral overexpression of SIMALR. NTN1 promoter activity is known to require HIF1α (hypoxia-inducible factor 1 subunit alpha), and our studies suggest that SIMALR may interact with HIF1α to regulate NTN1 transcription, thereby regulating macrophages apoptosis. SIMALR was found to be expressed in macrophages in human carotid atherosclerotic plaques of symptomatic patients.

CONCLUSIONS

SIMALR is a nonconserved, human macrophage lncRNA expressed in atherosclerosis that suppresses macrophage apoptosis. SIMALR partners with HIF1α (hypoxia-inducible factor 1 subunit alpha) to regulate NTN1, which is a known macrophage survival factor. This work illustrates the importance of interrogating the functions of human lncRNAs and exploring their translational and therapeutic potential in human atherosclerosis.

Computational Methods to Study DNA:DNA:RNA Triplex Formation by lncRNAs

Long non-coding RNAs (lncRNAs) impact cell function via numerous mechanisms. In the nucleus, interactions between lncRNAs and DNA and the consequent formation of non-canonical nucleic acid structures seems to be particularly relevant. Along with interactions between single-stranded RNA (ssRNA) and single-stranded DNA (ssDNA), such as R-loops, ssRNA can also interact with double-stranded DNA (dsDNA) to form DNA:DNA:RNA triplexes. A major challenge in the study of DNA:DNA:RNA triplexes is the identification of the precise RNA component interacting with specific regions of the dsDNA. As this is a crucial step towards understanding lncRNA function, there exist several computational methods designed to predict these sequences. This review summarises the recent progress in the prediction of triplex formation and highlights important DNA:DNA:RNA triplexes. In particular, different prediction tools (Triplexator, LongTarget, TRIPLEXES, Triplex Domain Finder, TriplexFFP, TriplexAligner and Fasim-LongTarget) will be discussed and their use exemplified by selected lncRNAs, whose DNA:DNA:RNA triplex forming potential was validated experimentally. Collectively, these tools revealed that DNA:DNA:RNA triplexes are likely to be numerous and make important contributions to gene expression regulation.

Identification of Long Noncoding RNAs That Exert Transcriptional Regulation by Forming RNA-DNA Triplexes in Prostate Cancer

Long noncoding RNAs (lncRNAs) are involved in transcriptional regulation, and their deregulation is associated with the development of various human cancers, including prostate cancer (PCa). However, their underlying mechanisms remain unclear. In this study, lncRNAs that interact with DNA and regulate mRNA transcription in PCa were screened and identified to promote PCa development. First, 4195 protein-coding genes (PCGs, mRNAs) were obtained from the The Cancer Genome Atlas (TCGA) database, in which 1148 lncRNAs were differentially expressed in PCa. Then, 44,270 pairs of co-expression relationships were calculated between 612 lncRNAs and 2742 mRNAs, of which 42,596 (96%) were positively correlated. Among the 612 lncRNAs, 392 had the potential to interact with the promoter region to form DNA:DNA:RNA triplexes, from which lncRNA AD000684.2(AC002128.1) was selected for further validation. AC002128.1 was highly expressed in PCa. Furthermore, AD000684.2 positively regulated the expression of the correlated genes. In addition, AD000684.2 formed RNA-DNA triplexes with the promoter region of the regulated genes. Functional assays also demonstrated that lncRNA AD000684.2 promotes cell proliferation and motility, as well as inhibits apoptosis, in PCa cell lines. The results suggest that AD000684.2 could positively regulate the transcription of target genes via triplex structures and serve as a candidate prognostic biomarker and target for new therapies in human PCa.

An overview of long noncoding RNAs: Biology, functions, therapeutics, analysis methods, and bioinformatics tools

Long noncoding RNAs (lncRNAs) are a diverse class of RNAs whose functions are widespread in all branches of life and have been the focus of attention in the last decade. While a huge number of lncRNAs have been identified, there is still much work to be done and plenty to be learned. In the current review, we begin with the biogenesis and function of lncRNAs as they are involved in the different cellular processes from regulating the architecture of chromosomes to controlling translation and post-translation modifications. Questions on how overexpression, mutations, or deficiency of lncRNAs can affect the cellular status and result in the pathogenesis of various human diseases are responded to. Besides, we allocate an overview of several studies, concerning the application of lncRNAs either as diagnostic and prognostic biomarkers or novel therapeutics. We also introduce the currently available techniques to explore details of lncRNAs such as their function, cellular localization, and structure. In the last section, as exponentially growing data in this area need to be gathered and organized in comprehensive databases, we have a particular focus on presenting general and specialized databases. Taken together, with this review, we aim to provide the latest information on different aspects of lncRNAs to highlight their importance in physiopathologic states and take a step towards helping future studies.

Implementation of a combined bioinformatics and experimental approach to address lncRNA mechanism of action: The example of NRIR

In this study, we demonstrate the benefit of applying combined strategies to analyze lncRNA action based on bioinformatics and experimental information. This strategy was developed to identify the molecular function of negative regulator of interferon response (NRIR), a type I interferon-stimulated gene (ISG), that we have previously demonstrated to be involved in the upregulation of a subset of ISGs in LPS-stimulated human monocytes. In this study, we provide experimental evidence that NRIR is localized in cellular nuclei, enriched on the chromatin fraction, and upregulates ISGs acting at the transcriptional level. In silico analysis of secondary structures identified distinct NRIR structural domains, comprising putative DNA- and protein-binding regions. In parallel, the presence of a putative DNA-binding domain in NRIR and the five putative NRIR-binding sites in the promoter of NRIR-target genes support the function of NRIR as a transcriptional regulator of its target genes. By use of integrated experimental/bioinformatics approaches, comprising database and literature mining together with in silico analysis of putative NRIR-binding proteins, we identified a list of eight transcription factors (TFs) shared by the majority of NRIR-target genes and simultaneously able to bind TF binding sites enriched in the NRIR-target gene promoters. Among these TFs, the predicted NRIR:STAT interactions were experimentally validated by RIP assay.

A universal model of RNA.DNA:DNA triplex formation accurately predicts genome-wide RNA-DNA interactions

RNA.DNA:DNA triple helix (triplex) formation is a form of RNA-DNA interaction which regulates gene expression but is difficult to study experimentally in vivo. This makes accurate computational prediction of such interactions highly important in the field of RNA research. Current predictive methods use canonical Hoogsteen base pairing rules, which whilst biophysically valid, may not reflect the plastic nature of cell biology. Here, we present the first optimization approach to learn a probabilistic model describing RNA-DNA interactions directly from motifs derived from triplex sequencing data. We find that there are several stable interaction codes, including Hoogsteen base pairing and novel RNA-DNA base pairings, which agree with in vitro measurements. We implemented these findings in TriplexAligner, a program that uses the determined interaction codes to predict triplex binding. TriplexAligner predicts RNA-DNA interactions identified in all-to-all sequencing data more accurately than all previously published tools in human and mouse and also predicts previously studied triplex interactions with known regulatory functions. We further validated a novel triplex interaction using biophysical experiments. Our work is an important step towards better understanding of triplex formation and allows genome-wide analyses of RNA-DNA interactions.

Fasim-LongTarget enables fast and accurate genome-wide lncRNA/DNA binding prediction

Many long noncoding RNAs (lncRNAs) can bind to DNA sequences proximal and distal to abundant genes, thereby regulating gene expression by recruiting epigenomic modification enzymes to binding sites. Because a lncRNA's target genes scattering in a genome have correlated functions, epigenetic analyses should often be genome-wide on both genome and transcriptome levels. Multiple tools have been developed for predicting lncRNA/DNA binding, but fast and accurate genome-wide prediction remains a challenge. Here we report Fasim-LongTarget (a revised version of LongTarget), compare its performance with TDF and LongTarget using the experimental data of the lncRNA MEG3, NEAT1, and MALAT1, and describe a case of genome-wide prediction. Fasim-LongTarget is as accurate as LongTarget and more accurate than TDF and is 200 times faster than LongTarget, making accurate genome-wide prediction feasible. The code is available on the Github website (https://github.com/LongTarget/Fasim-LongTarget), and the online service is available on the LongTarget website (https://lncRNA.smu.edu.cn).

Pluripotency exit is guided by the Peln1-mediated disruption of intrachromosomal architecture

The molecular circuitry that causes stem cells to exit from pluripotency remains largely uncharacterized. Using chromatin RNA in situ reverse transcription sequencing, we identified Peln1 as a novel chromatin RNA component in the promoter complex of Oct4, a stem cell master transcription factor gene. Peln1 was negatively associated with pluripotent status during somatic reprogramming. Peln1 overexpression caused E14 cells to exit from pluripotency, while Peln1 downregulation induced robust reprogramming. Mechanistically, we discovered that Peln1 interacted with the Oct4 promoter and recruited the DNA methyltransferase DNMT3A. By de novo altering the epigenotype in the Oct4 promoter, Peln1 dismantled the intrachromosomal loop that is required for the maintenance of pluripotency. Using RNA reverse transcription-associated trap sequencing, we showed that Peln1 targets multiple pathway genes that are associated with stem cell self-renewal. These findings demonstrate that Peln1 can act as a new epigenetic player and use a trans mechanism to induce an exit from the pluripotent state in stem cells.

A feedforward circuit between KLF5 and lncRNA KPRT4 contributes to basal-like breast cancer

Basal-like breast cancer (BLBC) is the most aggressive subtype of breast cancer with a poor prognosis. Long noncoding RNAs (lncRNAs) play critical roles in human cancers. Krüppel-like Factor 5 (KLF5) is a key oncogenic transcription factor in BLBC. However, the underlying mechanism of mutual regulation between KLF5 and lncRNA remains largely unknown. Here, we demonstrate that lncRNA KPRT4 promotes BLBC cell proliferation in vitro and in vivo. Mechanistically, KLF5 directly binds to the promoter of KPRT4 to promote KPRT4 transcription. Reciprocally, KPRT4 recruits the YB-1 transcription factor to the KLF5 promoter by interacting with YB-1 at its 5' domain and forming an RNA-DNA-DNA triplex structure at its 3' domain, resulting in enhanced transcription of KLF5 and ultimately establishing a feedforward circuit to promote cell proliferation. Moreover, the antisense oligonucleotide (ASO)-based therapy targeting KPRT4 substantially attenuated tumor growth in vivo. Clinically, the expression levels of YB-1, KLF5 and KPRT4 are positively correlated in clinical breast specimens. Together, our data suggest that KPRT4 is a major molecule for BLBC progression and that the feedforward circuit between KLF5 and KPRT4 may represent a potential therapeutic target in BLBC.

E2F1-induced lncRNA, EMSLR regulates lncRNA LncPRESS1

E2F1 induces hundreds of protein-coding genes influencing diverse signaling pathways but much less is known about its non-coding RNA targets. For identifying E2F1-dependent oncogenic long non-coding RNAs (lncRNAs), we carried out genome-wide transcriptome analysis and discovered an lncRNA, EMSLR, which is induced both in lung adenocarcinoma (LUAD) and lung squamous cell carcinoma (LUSC). EMSLR depletion blocks the cells in G1 phase and inhibits the clonogenic ability indicating that it is essential for the tumor-related phenotypes. We discovered that EMSLR represses the promoter activity of another lncRNA, LncPRESS1, which is located 6.9 kb upstream of EMSLR and they display an inverse expression pattern in lung cancer cell lines. Depletion of C-MYC results in downregulation of EMSLR and simultaneous upregulation of EMSLR target LncPRESS1, exemplifying how C-MYC and E2F1 signal transduction pathways control the network of lncRNA genes to modulate cell proliferation and differentiation.

The acetyltransferase p300 is recruited in trans to multiple enhancer sites by lncSmad7

The histone acetyltransferase p300 (also known as KAT3B) is a general transcriptional coactivator that introduces the H3K27ac mark on enhancers triggering their activation and gene transcription. Genome-wide screenings demonstrated that a large fraction of long non-coding RNAs (lncRNAs) plays a role in cellular processes and organ development although the underlying molecular mechanisms remain largely unclear (1,2). We found 122 lncRNAs that interacts directly with p300. In depth analysis of one of these, lncSmad7, is required to maintain ESC self-renewal and it interacts to the C-terminal domain of p300. lncSmad7 also contains predicted RNA-DNA Hoogsteen forming base pairing. Combined Chromatin Isolation by RNA precipitation followed by sequencing (ChIRP-seq) together with CRISPR/Cas9 mutagenesis of the target sites demonstrate that lncSmad7 binds and recruits p300 to enhancers in trans, to trigger enhancer acetylation and transcriptional activation of its target genes. Thus, these results unveil a new mechanism by which p300 is recruited to the genome.

Illuminating lncRNA Function Through Target Prediction

Most of the transcribed human genome codes for noncoding RNAs (ncRNAs), and long noncoding RNAs (lncRNAs) make for the lion's share of the human ncRNA space. Despite growing interest in lncRNAs, because there are so many of them, and because of their tissue specialization and, often, lower abundance, their catalog remains incomplete and there are multiple ongoing efforts to improve it. Consequently, the number of human lncRNA genes may be lower than 10,000 or higher than 200,000. A key open challenge for lncRNA research, now that so many lncRNA species have been identified, is the characterization of lncRNA function and the interpretation of the roles of genetic and epigenetic alterations at their loci. After all, the most important human genes to catalog and study are those that contribute to important cellular functions-that affect development or cell differentiation and whose dysregulation may play a role in the genesis and progression of human diseases. Multiple efforts have used screens based on RNA-mediated interference (RNAi), antisense oligonucleotide (ASO), and CRISPR screens to identify the consequences of lncRNA dysregulation and predict lncRNA function in select contexts, but these approaches have unresolved scalability and accuracy challenges. Instead-as was the case for better-studied ncRNAs in the past-researchers often focus on characterizing lncRNA interactions and investigating their effects on genes and pathways with known functions. Here, we focus most of our review on computational methods to identify lncRNA interactions and to predict the effects of their alterations and dysregulation on human disease pathways.

Spatiotemporal expression of long noncoding RNA Moshe modulates heart cell lineage commitment

The roles of long non-coding RNA (LncRNA) have been highlighted in various development processes including congenital heart defects (CHD). Here, we characterized the molecular function of LncRNA, Moshe (1010001N08ik-203), one of the Gata6 antisense transcripts located upstream of Gata6, which is involved in both heart development and the most common type of congenital heart defect, atrial septal defect (ASD). During mouse embryonic development, Moshe was first detected during the cardiac mesoderm stage (E8.5 to E9.5) where Gata6 is expressed and continues to increase at the atrioventricular septum (E12.5), which is involved in ASD. Functionally, the knock-down of Moshe during cardiogenesis caused significant repression of Nkx2.5 in cardiac progenitor stages and resulted in the increase in major SHF lineage genes, such as cardiac transcriptional factors (Isl1, Hand2, Tbx2), endothelial-specific genes (Cd31, Flk1, Tie1, vWF), a smooth muscle actin (a-Sma) and sinoatrial node-specific genes (Shox2, Tbx18). Chromatin Isolation by RNA Purification showed Moshe activates Nkx2.5 gene expression via direct binding to its promoter region. Of note, Moshe was conserved across species, including human, pig and mouse. Altogether, this study suggests that Moshe is a heart-enriched lncRNA that controls a sophisticated network of cardiogenesis by repressing genes in SHF via Nkx2.5 during cardiac development and may play an important role in ASD.

Long Noncoding RNA MIAT Controls Advanced Atherosclerotic Lesion Formation and Plaque Destabilization

Supplemental Digital Content is available in the text.

Long noncoding RNAs (lncRNAs) are important regulators of biological processes involved in vascular tissue homeostasis and disease development. The present study assessed the functional contribution of the lncRNA myocardial infarction-associated transcript (MIAT) to atherosclerosis and carotid artery disease.

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.

Tumor-derived exosomal BCYRN1 activates WNT5A/VEGF-C/VEGFR3 feedforward loop to drive lymphatic metastasis of bladder cancer

BACKGROUND

Patients with lymph node (LN) metastatic bladder cancer (BCa) present with extremely poor prognosis. BCa-derived exosomes function as crucial bioactive cargo carriers to mediate the signal transduction in tumor microenvironment triggering tumor metastasis. However, the mechanisms underlying exosome-mediated LN metastasis in BCa are unclear.

METHODS

We conducted the high-throughput sequencing to explore the expression profile of long noncoding RNA (lncRNA) in urinary exosomes (urinary-EXO) from patients with BCa and further evaluated the clinical relevance of exosomal lncRNA BCYRN1 in a larger 210-case cohort. The functional role of exosomal BCYRN1 was evaluated through the migration and tube formation assays in vitro and the footpad-popliteal LN metastasis model in vivo. RNA pull-down assays, luciferase assays, and actinomycin assays were conducted to detect the regulatory mechanism of exosomal BCYRN1.

RESULTS

LncRNA BCYRN1 was substantially upregulated in urinary-EXO from patients with BCa, and associated with the LN metastasis of BCa. We demonstrated that exosomal BCYRN1 markedly promoted tube formation and migration of human lymphatic endothelial cells (HLECs) in vitro and lymphangiogenesis and LN metastasis of BCa in vivo. Mechanistically, BCYRN1 epigenetically upregulated WNT5A expression by inducing hnRNPA1-associated H3K4 trimethylation in WNT5A promoter, which activated Wnt/β-catenin signaling to facilitate the secretion of VEGF-C in BCa. Moreover, exosomal BCYRN1 was transmitted to HLECs to stabilize the VEGFR3 mRNA and thus formed an hnRNPA1/WNT5A/VEGFR3 feedforward regulatory loop, ultimately promoting the lymphatic metastasis of BCa. Importantly, blocking VEGFR3 with specific inhibitor, SAR131675 significantly impaired exosomal BCYRN1-induced the LN metastasis in vivo. Clinically, exosomal BCYRN1 was positively associated with the shorter survival of BCa patients and identified as a poor prognostic factor of patients.

CONCLUSION

Our results uncover a novel mechanism by which exosomal BCYRN1 synergistically enhances VEGF-C/VEGFR3 signaling-induced lymphatic metastasis of BCa, indicating that BCYRN1 may serve as an encouraging therapeutic target for patients with BCa.

LncRNA LINC00525 suppresses p21 expression via mRNA decay and triplex-mediated changes in chromatin structure in lung adenocarcinoma

BACKGROUND

Emerging evidence suggests that long noncoding RNAs (lncRNAs) play crucial roles in various cancers. In the present study, we aim to investigate the function and molecular mechanism of an up-regulated and survival-associated lncRNA, LINC00525, in lung adenocarcinoma (LUAD).

METHODS

The expression level of LINC00525 in tissues was determined by quantitative reverse transcription polymerase chain reaction (RT-qPCR) and in situ hybridization (ISH). The functional role of LINC00525 in LUAD was investigated using gain-and loss-of-function approaches, both in vivo and in vitro. RNA pull-down, RNA immunoprecipitation (RIP), chromatin immunoprecipitation (ChIP), triplex-capture assay, dual-luciferase assay, gene expression microarray, and bioinformatics analysis were used to investigate the potential underlying mechanisms involved.

RESULTS

LINC00525 is highly expressed in LUAD cells and tissues. Survival analysis indicated that upregulation of LINC00525 was associated with poor prognosis in patients with LUAD patients. Knockdown of LINC00525 inhibited cell proliferation and cell cycle progression in vitro. In xenograft models, LINC00525 knockdown suppressed tumor growth and tumorigenesis of tumor-bearing mice. Mechanistically, LINC00525 epigenetically suppressed p21 transcription by guiding Enhancer Of Zeste 2 Polycomb Repressive Complex 2 Subunit (EZH2) to the p21 promoter through an formation of RNA-DNA triplex with the p21 promoter, leading to increased trimethylation of lysine 27 on histone 3 (H3K27me3) of the p21 promoter. In addition, LINC00525 repressed p21 expression post-transcriptionally by enhancing p21 mRNA decay. LINC00525 promoted p21 mRNA decay by competitively binding to RNA Binding Motif Single Stranded Interacting Protein 2 (RBMS2).

CONCLUSION

Our findings demonstrate that LINC00525 promotes the progression of LUAD by reducing the transcription and stability of p21 mRNA in concert with EZH2 and RBMS2, thus suggesting that LINC00525 may be a potential therapeutic target for clinical intervention in LUAD.

Long non-coding RNA histone deacetylase 4 antisense RNA 1 (HDAC4-AS1) inhibits HDAC4 expression in human ARPE-19 cells with hypoxic stress

Age-related macular degeneration (AMD) is resulted from choroidal neovascularization (CNV)-mediated cicatrization and vision loss. The sustained retinal hypoxia in retinal pigment epithelium (RPE) cells was reported to contribute to CNV. However, the underlying genetic regulatory network of hypoxia response in RPE is not fully understood. In this study, human ARPE-19 RPE cells were cultured under the anoxia for 24 h and later re-oxygenated in normoxia. Then the transcriptome was investigated via high throughput sequencing. We observed that long non-coding RNA (lncRNA) histone deacetylase 4 antisense RNA 1 (HDAC4-AS1) was increased in hypoxic condition compared to normal control and decreased after re-oxygenation addition, while the change of HDAC4 expression was reduced in hypoxic condition compared to normal control and up-regulated after re-oxygenation addition in ARPE-19 cells. Furthermore, HDAC4-AS1 knockdown could suppress the transcription activity of HDAC4 only in hypoxia condition, and fluorescence in situ hybridization and pull down assay indicated that transcripts of HDAC4-AS1 could substantially bind to the promoter of HDAC4 and facilitate the recruitment of HIF-1α. Finally, we also determined the specific regions of HDAC4-AS1 that contribute to the interaction with HIF-1α and the promoter of HDAC4. Taken together, these outcomes declared that HDAC4-AS1 could inhibit HDAC4 expression through regulating HIF-1α in human ARPE-19 cells with hypoxic stress.

Aging-associated lncRNAs are evolutionarily conserved and participate in NFκB signaling

The transcriptome undergoes global changes during aging, including both protein-coding and noncoding RNAs. Using comparative genomics, we identify aging-associated long noncoding RNAs (lncRNAs) that are under evolutionary constraint and are more conserved than lncRNAs that do not change with age. Aging-associated lncRNAs are enriched for functional elements, including binding sites for RNA-binding proteins and transcription factors, in particular nuclear factor kappa B (NFκB). Using CRISPR screening, we discovered that 13 of the aging-associated lncRNAs were regulators of the NFκB pathway, and we named this family 'NFκB modulating aging-related lncRNAs (NFKBMARLs)'. Further characterization of NFκBMARL-1 reveals it can be traced to 29 Ma before humans and is induced by NFκB during aging, inflammation and senescence. Reciprocally, NFκBMARL-1 directly regulates transcription of the NFκB inhibitor NFKBIZ in cis within the same topologically associated domain by binding to the NFKBIZ enhancer and recruiting RELA to the NFKBIZ promoter. These findings reveal many aging-associated lncRNAs are evolutionarily conserved components of the NFκB pathway.

SUMOylation promotes extracellular vesicle-mediated transmission of lncRNA ELNAT1 and lymph node metastasis in bladder cancer

Small ubiquitin-like modifier (SUMO) binding (termed SUMOylation) emerged as the inducer for the sorting of bioactive molecules into extracellular vesicles (EVs), triggering lymphangiogenesis and further driving tumor lymph node (LN) metastasis, but the precise mechanisms remain largely unclear. Here, we show that bladder cancer (BCa) cell-secreted EVs mediated intercellular communication with human lymphatic endothelial cells (HLECs) through transmission of the long noncoding RNA ELNAT1 and promoted lymphangiogenesis and LN metastasis in a SUMOylation-dependent manner in both cultured BCa cell lines and mouse models. Mechanistically, ELNAT1 induced UBC9 overexpression to catalyze the SUMOylation of hnRNPA1 at the lysine 113 residue, which mediated recognition of ELNAT1 by the endosomal sorting complex required for transport (ESCRT) and facilitated its packaging into EVs. EV-mediated ELNAT1 was specifically transmitted into HLECs and epigenetically activated SOX18 transcription to induce lymphangiogenesis. Importantly, blocking the SUMOylation of tumor cells by downregulating UBC9 expression markedly reduced lymphatic metastasis in EV-mediated, ELNAT1-treated BCa in vivo. Clinically, EV-mediated ELNAT1 was correlated with LN metastasis and a poor prognosis for patients with BCa. These findings highlight a molecular mechanism whereby the EV-mediated ELNAT1/UBC9/SOX18 regulatory axis promotes lymphangiogenesis and LN metastasis in BCa in a SUMOylation-dependent manner and implicate ELNAT1 as an attractive therapeutic target for LN metastatic BCa.

RNA:DNA triple helices: from peculiar structures to pervasive chromatin regulators

The genomes of complex eukaryotes largely contain non-protein-coding DNA, which is pervasively transcribed into a plethora of non-coding RNAs (ncRNAs). The functional importance of many of these ncRNAs has been investigated in the last two decades, revealing their crucial and multifaceted roles in chromatin regulation. A common mode of action of ncRNAs is the recruitment of chromatin modifiers to specific regions in the genome. Whereas many ncRNA-protein interactions have been characterised in detail, binding of ncRNAs to their DNA target sites is much less understood. Recently developed RNA-centric methods have mapped the genome-wide distribution of ncRNAs, however, how ncRNAs achieve locus-specificity remains mainly unresolved. In terms of direct RNA-DNA interactions, two kinds of triple-stranded structures can be formed: R-loops consisting of an RNA:DNA hybrid and a looped out DNA strand, and RNA:DNA triple helices (triplexes), in which the RNA binds to the major groove of the DNA double helix by sequence-specific Hoogsteen base pairing. In this essay, we will review the current knowledge about RNA:DNA triplexes, summarising triplex formation rules, detection methods, and ncRNAs reported to engage in triplexes. While the functional characterisation of RNA:DNA triplexes is still anecdotal, recent advances in high-throughput and computational analyses indicate their widespread distribution in the genome. Thus, we are witnessing a paradigm shift in the appreciation of RNA:DNA triplexes, away from exotic structures towards a prominent mode of ncRNA-chromatin interactions.

DeepLGP: a novel deep learning method for prioritizing lncRNA target genes

MOTIVATION

Although long non-coding RNAs (lncRNAs) have limited capacity for encoding proteins, they have been verified as biomarkers in the occurrence and development of complex diseases. Recent wet-lab experiments have shown that lncRNAs function by regulating the expression of protein-coding genes (PCGs), which could also be the mechanism responsible for causing diseases. Currently, lncRNA-related biological data are increasing rapidly. Whereas, no computational methods have been designed for predicting the novel target genes of lncRNA.

RESULTS

In this study, we present a graph convolutional network (GCN) based method, named DeepLGP, for prioritizing target PCGs of lncRNA. First, gene and lncRNA features were selected, these included their location in the genome, expression in 13 tissues and miRNA-mediated lncRNA-gene pairs. Next, GCN was applied to convolve a gene interaction network for encoding the features of genes and lncRNAs. Then, these features were used by the convolutional neural network for prioritizing target genes of lncRNAs. In 10-cross validations on two independent datasets, DeepLGP obtained high area under curves (0.90-0.98) and area under precision-recall curves (0.91-0.98). We found that lncRNA pairs with high similarity had more overlapped target genes. Further experiments showed that genes targeted by the same lncRNA sets had a strong likelihood of causing the same diseases, which could help in identifying disease-causing PCGs.

AVAILABILITY AND IMPLEMENTATION

https://github.com/zty2009/LncRNA-target-gene.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

linc-AAM Facilitates Gene Expression Contributing to Macrophage Activation and Adaptive Immune Responses

Although various long noncoding RNAs (lncRNAs) are specifically expressed in activated macrophages, their in vivo functions and mechanisms of action are largely unexplored. Here, we identify a long intergenic noncoding RNA associated with activated macrophage (linc-AAM) and elucidate its function and mechanisms. linc-AAM is highly expressed in activated macrophages. In vitro function analysis reveals that linc-AAM facilitates macrophage activation and promotes the expression of immune response genes (IRGs). In mechanisms, linc-AAM interacts with heterogeneous nuclear ribonucleoprotein L (hnRNPL) via two CACACA motifs, resulting in its dissociation from histone H3 to activate chromatin and facilitate transcription of IRGs. Of note, linc-AAM knockout (KO) mice manifest impaired antigen-specific cellular and humoral immune responses to ovalbumin (OVA) in vivo. Altogether, the results uncover a mechanism of lncRNA in modulating hnRNPL function and confirm that linc-AAM acts as a transcription enhancer to activate macrophages and promote adaptive immunity.

Deep Learning in LncRNAome: Contribution, Challenges, and Perspectives

Long non-coding RNAs (lncRNA), the pervasively transcribed part of the mammalian genome, have played a significant role in changing our protein-centric view of genomes. The abundance of lncRNAs and their diverse roles across cell types have opened numerous avenues for the research community regarding lncRNAome. To discover and understand lncRNAome, many sophisticated computational techniques have been leveraged. Recently, deep learning (DL)-based modeling techniques have been successfully used in genomics due to their capacity to handle large amounts of data and produce relatively better results than traditional machine learning (ML) models. DL-based modeling techniques have now become a choice for many modeling tasks in the field of lncRNAome as well. In this review article, we summarized the contribution of DL-based methods in nine different lncRNAome research areas. We also outlined DL-based techniques leveraged in lncRNAome, highlighting the challenges computational scientists face while developing DL-based models for lncRNAome. To the best of our knowledge, this is the first review article that summarizes the role of DL-based techniques in multiple areas of lncRNAome.

Deep learning based DNA:RNA triplex forming potential prediction

BACKGROUND

Long non-coding RNAs (lncRNAs) can exert functions via forming triplex with DNA. The current methods in predicting the triplex formation mainly rely on mathematic statistic according to the base paring rules. However, these methods have two main limitations: (1) they identify a large number of triplex-forming lncRNAs, but the limited number of experimentally verified triplex-forming lncRNA indicates that maybe not all of them can form triplex in practice, and (2) their predictions only consider the theoretical relationship while lacking the features from the experimentally verified data.

RESULTS

In this work, we develop an integrated program named TriplexFPP (Triplex Forming Potential Prediction), which is the first machine learning model in DNA:RNA triplex prediction. TriplexFPP predicts the most likely triplex-forming lncRNAs and DNA sites based on the experimentally verified data, where the high-level features are learned by the convolutional neural networks. In the fivefold cross validation, the average values of Area Under the ROC curves and PRC curves for removed redundancy triplex-forming lncRNA dataset with threshold 0.8 are 0.9649 and 0.9996, and these two values for triplex DNA sites prediction are 0.8705 and 0.9671, respectively. Besides, we also briefly summarize the cis and trans targeting of triplexes lncRNAs.

CONCLUSIONS

The TriplexFPP is able to predict the most likely triplex-forming lncRNAs from all the lncRNAs with computationally defined triplex forming capacities and the potential of a DNA site to become a triplex. It may provide insights to the exploration of lncRNA functions.

Towards a comprehensive pipeline to identify and functionally annotate long noncoding RNA (lncRNA)

Long noncoding RNAs (lncRNAs) are implicated in various genetic diseases and cancer, attributed to their critical role in gene regulation. They are a divergent group of RNAs and are easily differentiated from other types with unique characteristics, functions, and mechanisms of action. In this review, we provide a list of some of the prominent data repositories containing lncRNAs, their interactome, and predicted and validated disease associations. Next, we discuss various wet-lab experiments formulated to obtain the data for these repositories. We also provide a critical review of in silico methods available for the identification purpose and suggest techniques to further improve their performance. The bulk of the methods currently focus on distinguishing lncRNA transcripts from the coding ones. Functional annotation of these transcripts still remains a grey area and more efforts are needed in that space. Finally, we provide details of current progress, discuss impediments, and illustrate a roadmap for developing a generalized computational pipeline for comprehensive annotation of lncRNAs, which is essential to accelerate research in this area.

Silencing of Long Non-coding RNA GAS5 Suppresses Neuron Cell Apoptosis and Nerve Injury in Ischemic Stroke Through Inhibiting DNMT3B-Dependent MAP4K4 Methylation

Ischemic stroke is associated with various physiological and pathological processes including neuronal apoptosis. Growth-arrest-specific transcript 5 (GAS5), a long non-coding RNA (lncRNA), has been recently reported to affect ischemic stroke-induced neuron apoptosis, while its mechanisms remain largely undefined. Through in silico analysis, GAS5 was predicted to interact with the promoter of MAP4K4. The aim of the present study was therefore to investigate the possible role of GAS5 in the progression of ischemic stroke via regulation of mitogen-activated protein kinase kinase kinase kinase 4 (MAP4K4) methylation. The expression of MAP4K4 was found to be lowly expressed in the clinical samples collected from 55 patients. MAP4K4 was suggested to be methylated in an in vitro model of oxygen-glucose deprivation (OGD)-treated mouse primary cortical neurons, while its overexpression could inhibit OGD-induced neuronal apoptosis. A series of dual-luciferase reporter, RIP, RNA pull-down, ChIP MSP, and BSP assays confirmed that GAS5 significantly induced MAP4K4 methylation and downregulated MAP4K4 expression through the recruitment of DNA methyltransferase 3B (DNMT3B). An in vivo ischemic stroke model was developed using middle cerebral artery occlusion (MCAO). Upregulation of GAS5 promoted OGD-induced neuronal apoptosis in the in vitro model and increased cerebral infarction size and neurological score in the in vivo model by reducing MAP4K4 expression. Collectively, the present study highlights that silencing GAS5 may inhibit neuronal apoptosis and improve neurological function in ischemic stroke by suppressing DNMT3B-mediated MAP4K4 methylation, which contributes to better understanding of the pathologies of ischemic stroke and development of novel therapeutic options for this disease.

Long noncoding RNA functionality in imprinted domain regulation

Genomic imprinting is a parent-of-origin dependent phenomenon that restricts transcription to predominantly one parental allele. Since the discovery of the first long noncoding RNA (lncRNA), which notably was an imprinted lncRNA, a body of knowledge has demonstrated pivotal roles for imprinted lncRNAs in regulating parental-specific expression of neighboring imprinted genes. In this Review, we will discuss the multiple functionalities attributed to lncRNAs and how they regulate imprinted gene expression. We also raise unresolved questions about imprinted lncRNA function, which may lead to new avenues of investigation. This Review is dedicated to the memory of Denise Barlow, a giant in the field of genomic imprinting and functional lncRNAs. With her passion for understanding the inner workings of science, her indominable spirit and her consummate curiosity, Denise blazed a path of scientific investigation that made many seminal contributions to genomic imprinting and the wider field of epigenetic regulation, in addition to inspiring future generations of scientists.

LncRNA-SLC16A1-AS1 induces metabolic reprogramming during Bladder Cancer progression as target and co-activator of E2F1

Long non-coding RNAs (lncRNAs) have emerged as integral components of E2F1-regulated gene regulatory networks (GRNs), but their implication in advanced or treatment-refractory malignancy is unknown. Methods: We combined high-throughput transcriptomic approaches with bioinformatics and structure modeling to search for lncRNAs that participate in E2F1-activated prometastatic GRNs and their phenotypic targets in the highly-relevant case of E2F1-driven aggressive bladder cancer (BC). RNA immunoprecipitation was performed to verify RNA-protein interactions. Functional analyses including qRT-PCR, immunoblotting, luciferase assays and measurement of extracellular fluxes were conducted to validate expression and target gene regulation. Results: We identified E2F1-responsive lncRNA-SLC16A1-AS1 and its associated neighboring protein-coding gene, SLC16A1/MCT1, which both promote cancer invasiveness. Mechanistically, upon E2F1-mediated co-transactivation of the gene pair, SLC16A1-AS1 associates with E2F1 in a structure-dependent manner and forms an RNA-protein complex that enhances SLC16A1/MCT1 expression through binding to a composite SLC16A1-AS1:E2F1-responsive promoter element. Moreover, SLC16A1-AS1 increases aerobic glycolysis and mitochondrial respiration and fuels ATP production by fatty acid β-oxidation. These metabolic changes are accompanied by alterations in the expression of the SLC16A1-AS1:E2F1-responsive gene PPARA, a key mediator of fatty acid β-oxidation. Conclusions: Our results unveil a new gene regulatory program by which E2F1-induced lncRNA-SLC16A1-AS1 forms a complex with its transcription factor that promotes cancer metabolic reprogramming towards the acquisition of a hybrid oxidative phosphorylation/glycolysis cell phenotype favoring BC invasiveness.

Exosomal long noncoding RNA LNMAT2 promotes lymphatic metastasis in bladder cancer

Patients with bladder cancer (BCa) with clinical lymph node (LN) metastasis have an extremely poor prognosis. VEGF-C has been demonstrated to play vital roles in LN metastasis in BCa. However, approximately 20% of BCa with LN metastasis exhibits low VEGF-C expression, suggesting a VEGF-C-independent mechanism for LN metastasis of BCa. Herein, we demonstrate that BCa cell-secreted exosome-mediated lymphangiogenesis promoted LN metastasis in BCa in a VEGF-C-independent manner. We identified an exosomal long noncoding RNA (lncRNA), termed lymph node metastasis-associated transcript 2 (LNMAT2), that stimulated human lymphatic endothelial cell (HLEC) tube formation and migration in vitro and enhanced tumor lymphangiogenesis and LN metastasis in vivo. Mechanistically, LNMAT2 was loaded to BCa cell-secreted exosomes by directly interacting with heterogeneous nuclear ribonucleoprotein A2B1 (hnRNPA2B1). Subsequently, exosomal LNMAT2 was internalized by HLECs and epigenetically upregulated prospero homeobox 1 (PROX1) expression by recruitment of hnRNPA2B1 and increasing the H3K4 trimethylation level in the PROX1 promoter, ultimately resulting in lymphangiogenesis and lymphatic metastasis. Therefore, our findings highlight a VEGF-C-independent mechanism of exosomal lncRNA-mediated LN metastasis and identify LNMAT2 as a therapeutic target for LN metastasis in BCa.