A Comprehensive Characterization of the Function of LincRNAs in Transcriptional Regulation Through Long-Range Chromatin Interactions

Decoding the complexity of coding and non-coding RNAs across maize anther development at the isoform level

Anther is a key male reproductive organ that is essential for the plant life cycle, from the sporophyte to the gametophyte generation. To explore isoform-level transcriptional landscape of developing anthers in maize (Zea mays L.), we analyzed Iso-Seq data from anthers collected at 10 developmental stages, together with strand-specific RNA-seq, CAGE-seq, and PAS-seq data. Of the 152,026 high-confidence full-length isoforms identified, 68.8% have not been described; these include 22,365 isoforms that originate from previously unannotated loci and 82,167 novel isoforms that originate from annotated protein-coding genes. Using our newly developed strategy to detect dynamic expression patterns of isoforms, we identified 13,899 differentially variable regions (DVRs); surprisingly, 1,275 genes contain more than two DVRs, revealing highly efficient utilization of limited genic regions. We identified 7,876 long non-coding RNAs (lncRNAs) from 4,098 loci, most of which were preferentially expressed during cell differentiation and meiosis. We also detected 371 long-range interactions involving intergenic lncRNAs (lincRNAs); interestingly, 243 were lincRNA-gene ones, and the interacting genes were highly expressed in anthers, suggesting that many potential lncRNA regulators of key genes are required for anther development. This study provides valuable resources and fundamental information for studying the essential transcripts of key genes during anther development.

Simultaneous profiling of chromatin-associated RNA at targeted DNA loci and RNA-RNA Interactions through TaDRIM-seq

Eukaryotic genomes are extensively transcribed into various types of RNAs, many of which are physically associated with chromatin in cis at their transcription sites or in trans to other genomic loci. Emerging roles have been uncovered for these chromatin-associated RNAs (caRNAs) in gene regulation and genome organization, yet they remain challenging to interrogate. Here, we present TaDRIM-seq, a technique employing Protein G (PG)-Tn5-targeted DNA elements and in situ proximity ligation to concurrently probe caRNAs across diverse genomic regions as well as global RNA-RNA interactions within intact nuclei. Notably, this approach diminishes required cell inputs, minimizes hands-on time compared to established methodologies, and is compatible in both mammalian cells and plants. Using this technique, we identify extensive caRNAs at DNA anchor regions associated with chromatin loops and reveal diurnal variation in RNA-DNA and RNA-RNA connectivity networks within rice.

The regulatory activities of MALAT1 in the development of bone and cartilage diseases

Long non-coding RNAs (lncRNAs) have been comprehensively implicated in various cellular functions by mediating transcriptional or post-transcriptional activities. MALAT1 is involved in the differentiation, proliferation, and apoptosis of multiple cell lines, including BMSCs, osteoblasts, osteoclasts, and chondrocytes. Interestingly, MALAT1 may interact with RNAs or proteins, regulating cellular processes. Recently, MALAT1 has been reported to be associated with the development of bone and cartilage diseases by orchestrating the signaling network. The involvement of MALAT1 in the pathological development of bone and cartilage diseases makes it available to be a potential biomarker for clinical diagnosis or prognosis. Although the potential mechanisms of MALAT1 in mediating the cellular processes of bone and cartilage diseases are still needed for further elucidation, MALAT1 shows great promise for drug development.

Targeting the Expression of Long Noncoding RNAs in Murine Satellite Cells from Single Myofibers

LncRNAs have been recently implicated in the epigenetic control of muscle differentiation and their functional characterization has traditionally relied upon in vitro models of myogenic differentiation. However, the use of experimental paradigms to specifically target lncRNAs expression in muscle stem cells (MuSCs), also known as satellite cells, represents an important requisite to interrogate their function in more physiological contexts. Since isolation and culture of single myofibers preserves satellite cells within their physiological niche underneath the surrounding basal lamina, this procedure represents the optimal approach to follow satellite cell dynamics ex-vivo, such as activation from quiescence, expansion of committed progenitors, differentiation, and self-renewal. Here, we detail an optimized protocol to isolate viable single myofibers from the extensor digitorum longus (EDL) skeletal muscle of adult mice and to manipulate the expression of lncRNAs by antisense LNA GapmeRs-mediated knock-down (KD). Furthermore, we describe a method of EdU incorporation that, coupled to lncRNA KD and subsequent immunofluorescence analysis of proliferating, differentiating, and satellite cell-specific markers, permits the inference of lncRNAs function on muscle stem cells dynamics. Graphic abstract: Graphical representation of the single myofiber isolation method. Experimental workflow showing the main steps of the protocol procedure: EDL muscle harvesting from the mouse hindlimb; EDL digestion into single myofibers; transfection with antisense oligos and culture for 96h; immunofluorescence protocol and image outcome.

Non-Coding Variants in Cancer: Mechanistic Insights and Clinical Potential for Personalized Medicine

The cancer genome is characterized by extensive variability, in the form of Single Nucleotide Polymorphisms (SNPs) or structural variations such as Copy Number Alterations (CNAs) across wider genomic areas. At the molecular level, most SNPs and/or CNAs reside in non-coding sequences, ultimately affecting the regulation of oncogenes and/or tumor-suppressors in a cancer-specific manner. Notably, inherited non-coding variants can predispose for cancer decades prior to disease onset. Furthermore, accumulation of additional non-coding driver mutations during progression of the disease, gives rise to genomic instability, acting as the driving force of neoplastic development and malignant evolution. Therefore, detection and characterization of such mutations can improve risk assessment for healthy carriers and expand the diagnostic and therapeutic toolbox for the patient. This review focuses on functional variants that reside in transcribed or not transcribed non-coding regions of the cancer genome and presents a collection of appropriate state-of-the-art methodologies to study them.

Identification and analysis of short open reading frames (sORFs) in the initially annotated noncoding RNA LINC00493 from human cells

Whole transcriptome analyses have revealed that mammalian genomes are massively transcribed, resulting in the production of huge numbers of transcripts with unknown functions (TUFs). Previous research has categorized most TUFs as noncoding RNAs (ncRNAs) because most previously studied TUFs do not encode open reading frames (ORFs) with biologically significant lengths [>100 amino acids (AAs)]. Recent studies, however, have reported that several transcripts harbouring small ORFs that encode peptides shorter than 100 AAs are translated and play important biological functions. Here, we examined the translational capacity of transcripts annotated as ncRNAs in human cells, and identified several hundreds of ribosome-associated transcripts previously annotated as ncRNAs. Ribosome footprinting and polysome profiling analyses revealed that 61 of them are potentially translatable. Among them, 45 were nonnonsense-mediated mRNA decay targets, suggesting that they are productive mRNAs. We confirmed the translation of one ncRNA, LINC00493, by luciferase reporter assaying and western blotting of a FLAG-tagged LINC00493 peptide. While proteomic analysis revealed that the LINC00493 peptide interacts with many mitochondrial proteins, immunofluorescence assays showed that its peptide is mitochondrially localized. Our findings indicate that some transcripts annotated as ncRNAs encode peptides and that unannotated peptides may perform important roles in cells.

Wnt-regulated lncRNA discovery enhanced by in vivo identification and CRISPRi functional validation

BACKGROUND

Wnt signaling is an evolutionarily conserved developmental pathway that is frequently hyperactivated in cancer. While multiple protein-coding genes regulated by Wnt signaling are known, the functional lncRNAs regulated by Wnt signaling have not been systematically characterized.

METHODS

We comprehensively mapped Wnt-regulated lncRNAs from an orthotopic Wnt-addicted pancreatic cancer model and examined the response of lncRNAs to Wnt inhibition between in vivo and in vitro cancer models. We further annotated and characterized these Wnt-regulated lncRNAs using existing genomic classifications (using data from FANTOM5) in the context of Wnt signaling and inferred their role in cancer pathogenesis (using GWAS and expression data from the TCGA). To functionally validate Wnt-regulated lncRNAs, we performed CRISPRi screens to assess their role in cancer cell proliferation both in vivo and in vitro.

RESULTS

We identified 3633 lncRNAs, of which 1503 were regulated by Wnt signaling in an orthotopic Wnt-addicted pancreatic cancer model. These lncRNAs were much more sensitive to changes in Wnt signaling in xenografts than in cultured cells. Our analysis suggested that Wnt signaling inhibition could influence the co-expression relationship of Wnt-regulated lncRNAs and their eQTL-linked protein-coding genes. Wnt-regulated lncRNAs were also implicated in specific gene networks involved in distinct biological processes that contribute to the pathogenesis of cancers. Consistent with previous genome-wide lncRNA CRISPRi screens, around 1% (13/1503) of the Wnt-regulated lncRNAs were found to modify cancer cell growth in vitro. This included CCAT1 and LINC00263, previously reported to regulate cancer growth. Using an in vivo CRISPRi screen, we doubled the discovery rate, identifying twice as many Wnt-regulated lncRNAs (25/1503) that had a functional effect on cancer cell growth.

CONCLUSIONS

Our study demonstrates the value of studying lncRNA functions in vivo, provides a valuable resource of lncRNAs regulated by Wnt signaling, and establishes a framework for systematic discovery of functional lncRNAs.

Small RNAs, Big Diseases

The past two decades have seen extensive research done to pinpoint the role of microRNAs (miRNAs) that have led to discovering thousands of miRNAs in humans. It is not, therefore, surprising to see many of them implicated in a number of common as well as rare human diseases. In this review article, we summarize the progress in our understanding of miRNA-related research in conjunction with different types of cancers and neurodegenerative diseases, as well as their potential in generating more reliable diagnostic and therapeutic approaches.

MALAT1 Long Non-Coding RNA: Functional Implications

The mammalian genome is pervasively transcribed and the functional significance of many long non-coding RNA (lncRNA) transcripts are gradually being elucidated. Metastasis Associated Lung Adenocarcinoma Transcript 1 (MALAT1) is one of the most well-studied lncRNAs. MALAT1 is a highly conserved nuclear retained lncRNA that is abundantly expressed in cells and tissues and has been shown to play a role in regulating genes at both the transcriptional and post-transcriptional levels in a context-dependent manner. However, Malat1 has been shown to be dispensable for normal development and viability in mice. Interestingly, accumulating evidence suggests that MALAT1 plays an important role in numerous diseases including cancer. Here, we discuss the current state-of-knowledge in regard to MALAT1 with respect to its function, role in diseases, and the potential therapeutic opportunities for targeting MALAT1 using antisense oligonucleotides and small molecules.

Architectural RNAs for Membraneless Nuclear Body Formation

Long noncoding RNAs (lncRNAs) are fundamental regulators of various cellular processes. A subset of lncRNAs, termed architectural RNAs (arcRNAs), function in the formation and maintenance of phase-separated membraneless organelles in multiple eukaryotic species. These membraneless organelles represent an important type of compartmentalization in the crowded cellular environment and have several distinct features. The NEAT1_2 lncRNA is a well-characterized arcRNA that functions as an essential scaffold of paraspeckle nuclear bodies. Here, we describe the biogenesis of paraspeckles on arcRNAs through phase separation, focusing on the specific functions of multiple NEAT1_2 RNA domains and their partner RNA-binding proteins. Finally, we present an updated model of paraspeckle formation and discuss future perspectives of research into arcRNA-instructed architectures of phase-separated nuclear bodies.

Systemic Actions of Breast Cancer Facilitate Functional Limitations

Breast cancer is a disease of a specific organ, but its effects are felt throughout the body. The systemic effects of breast cancer can lead to functional limitations in patients who suffer from muscle weakness, fatigue, pain, fibromyalgia, or many other dysfunctions, which hasten cancer-associated death. Mechanistic studies have identified quite a few molecular defects in skeletal muscles that are associated with functional limitations in breast cancer. These include circulating cytokines such as TNF-α, IL-1, IL-6, and TGF-β altering the levels or function of myogenic molecules including PAX7, MyoD, and microRNAs through transcriptional regulators such as NF-κB, STAT3, and SMADs. Molecular defects in breast cancer may also include reduced muscle mitochondrial content and increased extracellular matrix deposition leading to energy imbalance and skeletal muscle fibrosis. This review highlights recent evidence that breast cancer-associated molecular defects mechanistically contribute to functional limitations and further provides insights into therapeutic interventions in managing functional limitations, which in turn may help to improve quality of life in breast cancer patients.

Identifying lncRNA-mediated regulatory modules via ChIA-PET network analysis

Background

Although several studies have provided insights into the role of long non-coding RNAs (lncRNAs), the majority of them have unknown function. Recent evidence has shown the importance of both lncRNAs and chromatin interactions in transcriptional regulation. Although network-based methods, mainly exploiting gene-lncRNA co-expression, have been applied to characterize lncRNA of unknown function by means of ’guilt-by-association’, no strategy exists so far which identifies mRNA-lncRNA functional modules based on the 3D chromatin interaction graph.

Results

To better understand the function of chromatin interactions in the context of lncRNA-mediated gene regulation, we have developed a multi-step graph analysis approach to examine the RNA polymerase II ChIA-PET chromatin interaction network in the K562 human cell line. We have annotated the network with gene and lncRNA coordinates, and chromatin states from the ENCODE project. We used centrality measures, as well as an adaptation of our previously developed Markov State Models (MSM) clustering method, to gain a better understanding of lncRNAs in transcriptional regulation. The novelty of our approach resides in the detection of fuzzy regulatory modules based on network properties and their optimization based on co-expression analysis between genes and gene-lncRNA pairs. This results in our method returning more bona fide regulatory modules than other state-of-the art approaches for clustering on graphs.

Conclusions

Interestingly, we find that lncRNA network hubs tend to be significantly enriched in evolutionary conserved lncRNAs and enhancer-like functions. We validated regulatory functions for well known lncRNAs, such as MALAT1 and the enhancer-like lncRNA FALEC. In addition, by investigating the modular structure of bigger components we mine putative regulatory functions for uncharacterized lncRNAs.

Electronic supplementary material

The online version of this article (10.1186/s12859-019-2900-8) contains supplementary material, which is available to authorized users.

Discovery, Identification, and Functional Characterization of Plant Long Intergenic Noncoding RNAs After Virus Infection

Long intergenic noncoding RNAs (lincRNAs), which possess diverse features such as remodeling chromatin and genome architecture, RNA stabilization, and genome architecture, are important regulatory factors in plant genomes. They serve to fine-tune the expression of neighboring genes. Here, we describe a procedure of discovery, identification, and functional characterization of plant lincRNAs after virus infection. From high-throughput RNA-Seq transcriptome analysis, the noncoding RNA transcripts with significant fold changes (upregulation or downregulation) will be discovered and identified. The lincRNA of interest will be further confirmed and validated using rapid amplification of cDNA ends (RACE). In addition, functional characterization of the lincRNA will be followed up through overexpression and knockdown strategies.

Into the Wild: GWAS Exploration of Non-coding RNAs

Genome-wide association studies (GWAS) have proven a fundamental tool to identify common variants associated to complex traits, thus contributing to unveil the genetic components of human disease. Besides, the advent of GWAS contributed to expose unexpected findings that urged to redefine the framework of population genetics. First, loci identified by GWAS had small effect sizes and could only explain a fraction of the predicted heritability of the traits under study. Second, the majority of GWAS hits mapped within non-coding regions (such as intergenic or intronic regions) where new functional RNA species (such as lncRNAs or circRNAs) have started to emerge. Bigger cohorts, meta-analysis and technical improvements in genotyping allowed identification of an increased number of genetic variants associated to coronary artery disease (CAD) and cardiometabolic traits. The challenge remains to infer causal mechanisms by which these variants influence cardiovascular disease development. A tendency to assign potential causal variants preferentially to coding genes close to lead variants contributed to disregard the role of non-coding elements. In recent years, in parallel to an increased knowledge of the non-coding genome, new studies started to characterize disease-associated variants located within non-coding RNA regions. The upcoming of databases integrating single-nucleotide polymorphisms (SNPs) and non-coding RNAs together with novel technologies will hopefully facilitate the discovery of causal non-coding variants associated to disease. This review attempts to summarize the current knowledge of genetic variation within non-coding regions with a focus on long non-coding RNAs that have widespread impact in cardiometabolic diseases.

Genome-wide epigenetic landscape of pig lincRNAs and their evolution during porcine domestication

AIM

We aimed to identify previously unreported long intergenic noncoding RNAs (lincRNAs) in the porcine liver, an important metabolic tissue, and further illustrate the epigenomic landscapes and the evolution of lincRNAs.

MATERIALS & METHODS

We used porcine omics data and comprehensively analyzed and identified lincRNAs and their methylation, expression and evolutionary patterns during pig domestication.

RESULTS

LincRNAs exhibit highly methylated promoter and downstream regions, as well as lower expression levels and higher tissue specificity than protein-coding genes. We identified a batch of lincRNAs with selection signals that are associated with pig domestication, which are more highly expressed in the liver than in other tissues (19:10/8/6/3/2/1/1). Interestingly, the lincRNA linc-sscg1779 and its target gene C6, which is crucial in liver metabolism, are differentially expressed during pig domestication.

CONCLUSION

Although they may originate from noisy transcripts, lincRNAs may be subjected to artificial selection. This phenomenon implies the functional importance of lincRNAs in pig domestication.

LncRNAs in TGF-β-Driven Tissue Fibrosis

Transforming growth factor-β (TGF-β) is a crucial mediator in tissue fibrosis that promotes accumulation of extracellular matrix (ECM), myofibroblasts to epithelial–mesenchymal transition (EMT), endothelial-mesenchymal transition (EndoMT), and apoptosis via canonical and noncanonical signaling pathways. In the past decades, a number of microRNAs have been reported to participate in TGF-β-mediated tissue scarring; however, the roles of long noncoding RNAs (lncRNAs) in fibrogenesis remain largely unknown. Recently, emerging evidence has shown that lncRNAs are involved in the development of different diseases, including cancer, autoimmune diseases, cardiovascular diseases, and fibrotic diseases. In this review, we summarize the current updates of lncRNAs in TGF-β1-driven tissue fibrosis and discuss their therapeutic potential for the treatment of chronic fibrotic diseases.

CisPi: a transcriptomic score for disclosing cis-acting disease-associated lincRNAs

Motivation

Long intergenic noncoding RNAs (lincRNAs) have risen to prominence in cancer biology as new biomarkers of disease. Those lincRNAs transcribed from active cis-regulatory elements (enhancers) have provided mechanistic insight into cis-acting regulation; however, in the absence of an enhancer hallmark, computational prediction of cis-acting transcription of lincRNAs remains challenging. Here, we introduce a novel transcriptomic method: a cis-regulatory lincRNA-gene associating metric, termed 'CisPi'. CisPi quantifies the mutual information between lincRNAs and local gene expression regarding their response to perturbation, such as disease risk-dependence. To predict risk-dependent lincRNAs in neuroblastoma, an aggressive pediatric cancer, we advance this scoring scheme to measure lincRNAs that represent the minority of reads in RNA-Seq libraries by a novel side-by-side analytical pipeline.

Results

Altered expression of lincRNAs that stratifies tumor risk is an informative readout of oncogenic enhancer activity. Our CisPi metric therefore provides a powerful computational model to identify enhancer-templated RNAs (eRNAs), eRNA-like lincRNAs, or active enhancers that regulate the expression of local genes. First, risk-dependent lincRNAs revealed active enhancers, over-represented neuroblastoma susceptibility loci, and uncovered novel clinical biomarkers. Second, the prioritized lincRNAs were significantly prognostic. Third, the predicted target genes further inherited the prognostic significance of these lincRNAs. In sum, RNA-Seq alone is sufficient to identify disease-associated lincRNAs using our methodologies, allowing broader applications to contexts in which enhancer hallmarks are not available or show limited sensitivity.

Availability and implementation

The source code is available on request. The prioritized lincRNAs and their target genes are in the Supplementary Material.

Supplementary information

Supplementary data are available at Bioinformatics online.

Genome-Wide Association and Replication Study of Hepatotoxicity Induced by Antiretrovirals Alone or with Concomitant Anti-Tuberculosis Drugs

Drug-induced hepatotoxicity (DIH) is a common adverse event that is associated with both antiretroviral (ARV) and anti-tuberculosis drugs (ATD). Moreover, the genetic variations predisposing ARV- and ARV-ATD-induced liver toxicity in African populations are not well investigated, despite the two diseases being the major global health problems in sub-Saharan Africa. We performed a genome-wide association study (GWAS) and replication study to identify the genetic variants linked to the risk of developing DIH due to ARV drugs alone, and ARV-ATD co-treatment in Ethiopian HIV-positive patients. Treatment-naïve newly diagnosed HIV patients (n = 719) with or without tuberculosis (TB) co-infection were enrolled prospectively and received efavirenz-based ARV therapy with or without rifampicin-based short course ATD, respectively. Whole-genome genotyping was performed by using the Illumina Omni Express Exome Bead Chip genotyping array with 951,117 single nucleotide polymorphisms (SNPs) on a total of 41 cases of DIH, and 452 people without DIH (treatment tolerants). The replication study was carried out for 100 SNPs with the lowest p-values (top SNPs) by using an independent cohort consisting of 18 DIH cases and 208 treatment tolerants. We identified a missense SNP rs199650082 (2756G→A, R919Q, p = 1.4 × 10^-6, odds ratio [OR] = 18.2, 95% confidence interval [CI] = 7.1-46.9) in an endoplasmic reticulum to the nucleus signaling-1 (ERN1) gene on chromosome 17 to be associated with DIH in the ARV-only cohort. In the ARV-ATD co-treatment groups, rs4842407, a long intergenic noncoding RNAs (lincRNAs) transcript variant on chromosome 12, was associated with DIH (p = 5.3 × 10^-7, OR = 5.4, 95% CI = 2.8-10.3). These genetic variants that are putatively associated with DIH due to ARV drugs alone and ARV-ATD co-treatment establish a foundation for future personalized medicine in people with HIV and TB and call for larger studies in independent populations.

Non-coding RNAs in skeletal muscle regeneration

Following injury, skeletal muscles can regenerate from muscle specific stem cells, called satellite cells. Quiescent in uninjured muscles, satellite cells become activated, proliferate and differentiate into myotubes. Muscle regeneration occurs following distinct main overlapping phases, including inflammation, regeneration and maturation of the regenerated myofibers. Each step of muscle regeneration is orchestrated through complex signaling networks and gene regulatory networks, leading to the expression of specific set of genes in each concerned cell type. Apart from the well-established transcriptional mechanisms involving the myogenic regulatory factors of the MyoD family, increasing data indicate that each step of muscle regeneration is controlled by a wide range of non-coding RNAs. In this review, we discuss the role of two classes of non-coding RNAs (microRNAs and long non-coding RNAs) in the inflammatory, regeneration and maturation steps of muscle regeneration.

Biological and chemical guided isolation of 3,4-secograyanane diterpenoids from the roots of Pieris formosa

Seventeen new 3,4-secograyanoids (1–17), together with seven known compounds (18–24), were isolated from the roots of Pieris formosa.