LnCompare: gene set feature analysis for human long non-coding RNAs

Exploring the landscape of tools and resources for the analysis of long non-coding RNAs

In recent years, research on long non-coding RNAs (lncRNAs) has gained considerable attention due to the increasing number of newly identified transcripts. Several characteristics make their functional evaluation challenging, which called for the urgent need to combine molecular biology with other disciplines, including bioinformatics. Indeed, the recent development of computational pipelines and resources has greatly facilitated both the discovery and the mechanisms of action of lncRNAs. In this review, we present a curated collection of the most recent computational resources, which have been categorized into distinct groups: databases and annotation, identification and classification, interaction prediction, and structure prediction. As the repertoire of lncRNAs and their analysis tools continues to expand over the years, standardizing the computational pipelines and improving the existing annotation of lncRNAs will be crucial to facilitate functional genomics studies.

lncHUB2: aggregated and inferred knowledge about human and mouse lncRNAs

Long non-coding ribonucleic acids (lncRNAs) account for the largest group of non-coding RNAs. However, knowledge about their function and regulation is limited. lncHUB2 is a web server database that provides known and inferred knowledge about the function of 18 705 human and 11 274 mouse lncRNAs. lncHUB2 produces reports that contain the secondary structure fold of the lncRNA, related publications, the most correlated coding genes, the most correlated lncRNAs, a network that visualizes the most correlated genes, predicted mouse phenotypes, predicted membership in biological processes and pathways, predicted upstream transcription factor regulators, and predicted disease associations. In addition, the reports include subcellular localization information; expression across tissues, cell types, and cell lines, and predicted small molecules and CRISPR knockout (CRISPR-KO) genes prioritized based on their likelihood to up- or downregulate the expression of the lncRNA. Overall, lncHUB2 is a database with rich information about human and mouse lncRNAs and as such it can facilitate hypothesis generation for many future studies. The lncHUB2 database is available at https://maayanlab.cloud/lncHUB2. Database URL: https://maayanlab.cloud/lncHUB2.

A comprehensive survey of the approaches for pathway analysis using multi-omics data integration

Pathway analysis has been widely used to detect pathways and functions associated with complex disease phenotypes. The proliferation of this approach is due to better interpretability of its results and its higher statistical power compared with the gene-level statistics. A plethora of pathway analysis methods that utilize multi-omics setup, rather than just transcriptomics or proteomics, have recently been developed to discover novel pathways and biomarkers. Since multi-omics gives multiple views into the same problem, different approaches are employed in aggregating these views into a comprehensive biological context. As a result, a variety of novel hypotheses regarding disease ideation and treatment targets can be formulated. In this article, we review 32 such pathway analysis methods developed for multi-omics and multi-cohort data. We discuss their availability and implementation, assumptions, supported omics types and databases, pathway analysis techniques and integration strategies. A comprehensive assessment of each method's practicality, and a thorough discussion of the strengths and drawbacks of each technique will be provided. The main objective of this survey is to provide a thorough examination of existing methods to assist potential users and researchers in selecting suitable tools for their data and analysis purposes, while highlighting outstanding challenges in the field that remain to be addressed for future development.

Parsing the role of NSP1 in SARS-CoV-2 infection

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) leads to shutoff of protein synthesis, and nsp1, a central shutoff factor in coronaviruses, inhibits cellular mRNA translation. However, the diverse molecular mechanisms employed by nsp1 as well as its functional importance are unresolved. By overexpressing various nsp1 mutants and generating a SARS-CoV-2 mutant, we show that nsp1, through inhibition of translation and induction of mRNA degradation, targets translated cellular mRNA and is the main driver of host shutoff during infection. The propagation of nsp1 mutant virus is inhibited exclusively in cells with intact interferon (IFN) pathway as well as in vivo, in hamsters, and this attenuation is associated with stronger induction of type I IFN response. Therefore, although nsp1's shutoff activity is broad, it plays an essential role, specifically in counteracting the IFN response. Overall, our results reveal the multifaceted approach nsp1 uses to shut off cellular protein synthesis and uncover nsp1's explicit role in blocking the IFN response.

Parsing the role of NSP1 in SARS-CoV-2 infection

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the cause of the ongoing coronavirus disease 19 (COVID-19) pandemic. Despite its urgency, we still do not fully understand the molecular basis of SARS-CoV-2 pathogenesis and its ability to antagonize innate immune responses. SARS-CoV-2 leads to shutoff of cellular protein synthesis and over-expression of nsp1, a central shutoff factor in coronaviruses, inhibits cellular gene translation. However, the diverse molecular mechanisms nsp1 employs as well as its functional importance in infection are still unresolved. By overexpressing various nsp1 mutants and generating a SARS-CoV-2 mutant in which nsp1 does not bind ribosomes, we untangle the effects of nsp1. We uncover that nsp1, through inhibition of translation and induction of mRNA degradation, is the main driver of host shutoff during SARS-CoV-2 infection. Furthermore, we find the propagation of nsp1 mutant virus is inhibited specifically in cells with intact interferon (IFN) response as well as in-vivo , in infected hamsters, and this attenuation is associated with stronger induction of type I IFN response. This illustrates that nsp1 shutoff activity has an essential role mainly in counteracting the IFN response. Overall, our results reveal the multifaceted approach nsp1 uses to shut off cellular protein synthesis and uncover the central role it plays in SARS-CoV-2 pathogenesis, explicitly through blockage of the IFN response.

Transcriptomic analysis reveals myometrial topologically associated domains linked to onset of human term labor

Changes in cell phenotype are thought to occur through the expression of groups of co-regulated genes within topologically associated domains (TADs). In this paper we allocate genes expressed within the myometrium of the human uterus during the onset of term labor into TADs. Transformation of the myometrial cells of the uterus into a contractile phenotype during term human labor is the result of a complex interaction of different epigenomic and genomic layers. Recent work suggests that the transcription factor RELA lies at the top of this regulatory network. Using deep RNA sequencing (RNAseq) analysis of myometrial samples (n = 16) obtained at term from women undergoing Caesarean section prior to or after the onset of labor we have identified evidence for how other gene expression regulatory elements interact with transcription factors in the labor phenotype transition. Gene set enrichment analysis of our RNAseq data identified three modules of enriched genes (M1, M2 and M3), which in gene ontology studies are linked to matrix degradation, smooth muscle and immune gene signatures, respectively. These genes were predominantly located within chromosomal TADs suggesting co-regulation of expression. Our transcriptomic analysis also identified significant differences in the expression of long non-coding RNAs (lncRNA), microRNAs (miRNA) and transcription factors that were predicted to target genes within the TADs. Additionally, network analysis revealed 15 new lncRNA (MCM3AP-AS1, TUG1, MIR29B2CHG, HCG18, LINC00963, KCNQ1OT1, NEAT1, HELLPAR, SNHG16, NUTM2B-AS1, MALAT1, PSMA3-AS1, GABPB1-AS1, NORAD, NKILA) and four miRNA (mir-145, mir-223, mir-let-7a, mir-132) as top gene hubs with three transcription factors (NFKB1, RELA, ESR1) as master regulators. Together, these factors are likely to be involved in co-regulatory networks driving a myometrial transformation to generate an estrogen sensitive phenotype. We conclude that lncRNA and miRNA targeting the estrogen receptor 1 and nuclear factor kappa B pathways play a key role in the initiation of human labor. For the first time we perform an integrative analysis to present a multi-level genomic signature made of mRNA, ncRNA and transcription factors in the myometrium for spontaneous term labor.

Identification of lncRNAs based on different patterns of immune infiltration in gastric cancer

BACKGROUND

Gastric cancer is one of the most common malignant tumors in the world, which brings great challenges to people's life and health. The purpose of this study was to investigate immune related-lncRNAs and identify new biomarkers for the prognosis of gastric cancer (GC).

METHODS

We downloaded data from The Cancer Genome Atlas (TCGA) and used R software to determine the ESTIMATEScore, ImmuneScore, and StromalScore of each tumor sample. We performed prognostic analysis and identified the differentially expressed lnRNAs, which were then used to construct a prognostic model. Among the 44 hub genes in the competitive endogenous RNA (ceRNA) network, 3 differentially expressed genes were verified by qPCR.

RESULTS

Based on the degree of immune infiltration, cluster A had a higher ESTIMATEScore, ImmuneScore, and StromalScore and higher expression levels of PD-L1 (CD274) and CTLA4 than cluster B. Univariate Cox analysis was conducted for these differential lncRNAs, and 57 lncRNAs were found to have prognostic value (P<0.05). gene cluster A had a worse prognosis than gene cluster B (P=0.021). Then, a prognostic model was constructed. The low-risk group had a significantly higher survival rate. Finally, the qPCR results showed that the expression levels of BMPER, PRUNE2, and RBPMS2 were low in GC cell lines.

CONCLUSIONS

We identified a risk score of 19 lncRNAs as a prognostic marker of GC. There was a relationship between these 19 prognostic-related lncRNAs and the subtypes of infiltrating immune cells. An approach for predicting the prognosis of GC was therefore provided in this study.

Development and Validation of Genome Instability-Associated lncRNAs to Predict Prognosis and Immunotherapy of Patients With Hepatocellular Carcinoma

The pathophysiology of hepatocellular carcinoma (HCC) is prevalently related to genomic instability. However, research on the association of extensive genome instability lncRNA (GILnc) with the prognosis and immunotherapy of HCC remains scarce. We placed the top 25% of somatic mutations into the genetically unstable group and placed the bottom 25% of somatic mutations into the genetically stable group, and then to identify different expression of GILnc between the two groups. Then, LASSO was used to identify the most powerful prognostic GILnc, and a risk score for each patient was calculated according to the formula. Based on a computational frame, 245 different GILncs in HCC were identified. An eight GILnc model was successfully established to predict overall survival in HCC patients based on LASSO, then we divided HCC patients into high-risk and low-risk groups, and a significantly shorter overall survival in the high-risk group was observed compared to those in the low-risk group, and this was validated in GSE76427 and Tongji cohorts. GSEA revealed that the high-risk group was more likely to be enriched in cancer-specific pathways. Besides, the GILnc signature has greater prognostic significance than TP53 mutation status alone, and it is capable of identifying intermediate subtype groups existing with partial TP53 functionality in TP53 wild-type patients. Importantly, the high-risk group was associated with the therapeutic efficacy of PD-L1 blockade, suggesting that the development of potential drugs targeting these GILnc could aid the clinical benefits of immunotherapy. Finally, the GILnc signature model is better than the prediction performance of two recently published lncRNA signatures. In summary, we applied bioinformatics approaches to suggest that an eight GILnc model could serve as prognostic biomarkers to provide a novel direction to explore the pathogenesis of HCC.

A novel 8-genome instability-associated lncRNAs signature predicting prognosis and drug sensitivity in gastric cancer

BACKGROUND

Genome instability lncRNA (GILnc) is prevalently related with gastric cancer (GC) pathophysiology. However, the study on the relationship GILnc and prognosis and drug sensitivity of GC remains scarce.

METHOD

We extracted expression data of 375 GC patients from TCGA cohort and 205 GC patients from GSE26942 cohort. Then, lncRNA was separated from expression data, and systematically characterized the 8 marker lncRNAs using the LASSO method. Next, we constructed a GILnc model (GILnc score) to quantify the GILnc index of each GC patient. Finally, we analyzed the relationship between GILnc score and clinical traits including survival outcomes, TP53, and drug sensitivity of GC.

RESULTS

Based on a computational frame, 205 GILncs in GC has been identified. Then, a 8 GILncs was successfully established to predict overall survival in GC patients based on LASSO analysis, divided GC samples into high GILnc score and low GILnc score groups with significantly different outcome and was validated in multiple independent patient cohorts. Furthermore, GILnc model is better than the prediction performance of two recently published lncRNA signatures, and the high GILnc score group was more sensitive to mitomycin. Besides, the GILnc score has greater prognostic significance than TP53 mutation status alone and is capable of identifying intermediate subtype group existing with partial TP53 functionality in TP53 wild-type patients. Finally, GILnc signature as verified in GSE26942.

CONCLUSION

We applied bioinformatics approaches to suggest that a 8 GILnc signature could serve as prognostic biomarkers, and provide a novel direction to explore the pathogenesis of GC.

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

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

Cancer LncRNA Census 2 (CLC2): an enhanced resource reveals clinical features of cancer lncRNAs

Long non-coding RNAs (lncRNAs) play key roles in cancer and are at the vanguard of precision therapeutic development. These efforts depend on large and high-confidence collections of cancer lncRNAs. Here, we present the Cancer LncRNA Census 2 (CLC2). With 492 cancer lncRNAs, CLC2 is 4-fold greater in size than its predecessor, without compromising on strict criteria of confident functional/genetic roles and inclusion in the GENCODE annotation scheme. This increase was enabled by leveraging high-throughput transposon insertional mutagenesis screening data, yielding 92 novel cancer lncRNAs. CLC2 makes a valuable addition to existing collections: it is amongst the largest, contains numerous unique genes (not found in other databases) and carries functional labels (oncogene/tumour suppressor). Analysis of this dataset reveals that cancer lncRNAs are impacted by germline variants, somatic mutations and changes in expression consistent with inferred disease functions. Furthermore, we show how clinical/genomic features can be used to vet prospective gene sets from high-throughput sources. The combination of size and quality makes CLC2 a foundation for precision medicine, demonstrating cancer lncRNAs’ evolutionary and clinical significance.

LncSEA: a platform for long non-coding RNA related sets and enrichment analysis

Long non-coding RNAs (lncRNAs) have been proven to play important roles in transcriptional processes and various biological functions. Establishing a comprehensive collection of human lncRNA sets is urgent work at present. Using reference lncRNA sets, enrichment analyses will be useful for analyzing lncRNA lists of interest submitted by users. Therefore, we developed a human lncRNA sets database, called LncSEA, which aimed to document a large number of available resources for human lncRNA sets and provide annotation and enrichment analyses for lncRNAs. LncSEA supports >40 000 lncRNA reference sets across 18 categories and 66 sub-categories, and covers over 50 000 lncRNAs. We not only collected lncRNA sets based on downstream regulatory data sources, but also identified a large number of lncRNA sets regulated by upstream transcription factors (TFs) and DNA regulatory elements by integrating TF ChIP-seq, DNase-seq, ATAC-seq and H3K27ac ChIP-seq data. Importantly, LncSEA provides annotation and enrichment analyses of lncRNA sets associated with upstream regulators and downstream targets. In summary, LncSEA is a powerful platform that provides a variety of types of lncRNA sets for users, and supports lncRNA annotations and enrichment analyses. The LncSEA database is freely accessible at http://bio.liclab.net/LncSEA/index.php.

Gliadin, through the Activation of Innate Immunity, Triggers lncRNA NEAT1 Expression in Celiac Disease Duodenal Mucosa

Celiac disease (CD) is an autoimmune enteropathy arising in genetically predisposed subjects exposed to gluten, which activates both innate and adaptive immunity. Although the pathogenesis is common to all patients, the clinical spectrum is quite variable, and differences could be explained by gene expression variations. Among the factors able to affect gene expression, there are lncRNAs. We evaluated the expression profile of 87 lncRNAs in CD vs. healthy control (HC) intestinal biopsies by RT-qPCR array. Nuclear enriched abundant transcript 1 (NEAT1) and taurine upregulated gene 1 (TUG1) were detected as downregulated in CD patients at diagnosis, but their expression increased in biopsies of patients on a gluten-free diet (GFD) exposed to gluten. The increase in NEAT1 expression after gluten exposure was mediated by IL-15 and STAT3 activation and binding to the NEAT1 promoter, as demonstrated by gel shift assay. NEAT1 is localized in the nucleus and can regulate gene expression by sequestering transcription factors, and it has been implicated in immune regulation and control of cell proliferation. The demonstration of its regulation by gluten thus also supports the role of lncRNAs in CD and prompts further research on these RNAs as gene expression regulators.

Bioinformatics Approaches for Functional Prediction of Long Noncoding RNAs

There is accumulating evidence that long noncoding RNAs (lncRNAs) play crucial roles in biological processes and diseases. In recent years, computational models have been widely used to predict potential lncRNA-disease relations. In this chapter, we systematically describe various computational algorithms and prediction tools that have been developed to elucidate the roles of lncRNAs in diseases, coding potential/functional characterization, or ascertaining their involvement in critical biological processes as well as provide a comprehensive summary of these applications.

Annotation and curation of the causality information in LncRNADisease

Disease causative non-coding RNAs (ncRNAs) are of great importance in understanding a disease, for they directly contribute to the development or progress of a disease. Identifying the causative ncRNAs can provide vital implications for biomedical researches. In this work, we updated the long non-coding RNA disease database (LncRNADisease) with long non-coding RNA (lncRNA) causality information with manual annotations of the causal associations between lncRNAs/circular RNAs (circRNAs) and diseases by reviewing related publications. Of the total 11 568 experimental associations, 2297 out of 10 564 lncRNA-disease associations and 198 out of 1004 circRNA-disease associations were identified to be causal, whereas 635 lncRNAs and 126 circRNAs were identified to be causative for the development or progress of at least one disease. The updated information and functions of the database can offer great help to future researches involving lncRNA/circRNA-disease relationship. The latest LncRNADisease database is available at http://www.rnanut.net/lncrnadisease.

Target Enrichment Enables the Discovery of lncRNAs with Somatic Mutations or Altered Expression in Paraffin-Embedded Colorectal Cancer Samples

Simple Summary

Alterations in long noncoding RNAs and their mutations have been increasingly recognized in tumorogenesis and cancer progression awakening especial interest as potential novel cancer biomarkers and therapeutic targets. The use of adjuvant chemotherapy in stage II colorectal cancer patients is challenging, and new biomarkers are required to identify patients with high probability of relapse. We focused on translational potential of non-coding RNAs in colorectal cancer. In this study, we aim to validate a new tool which couples target enrichment and RNAseq for transcriptomics studies of lncRNAs in formalin-fixed paraffin embedded (FFPE) tissue samples. Our results show that this new approach efficiently detects lncRNAs and differences in their expression between healthy and tumor FFPE tissues, as well as somatic mutations in expressed lncRNAs, identifying novel lncRNAs as potential candidates for colorectal cancer. This new approach could represent a promising avenue that would reduce costs and enable more efficient translational research.

Abstract

Long non-coding RNAs (lncRNAs) play important roles in cancer and are potential new biomarkers or targets for therapy. However, given the low and tissue-specific expression of lncRNAs, linking these molecules to particular cancer types and processes through transcriptional profiling is challenging. Formalin-fixed, paraffin-embedded (FFPE) tissues are abundant resources for research but are prone to nucleic acid degradation, thereby complicating the study of lncRNAs. Here, we designed and validated a probe-based enrichment strategy to efficiently profile lncRNA expression in FFPE samples, and we applied it for the detection of lncRNAs associated with colorectal cancer (CRC). Our approach efficiently enriched targeted lncRNAs from FFPE samples, while preserving their relative abundance, and enabled the detection of tumor-specific mutations. We identified 379 lncRNAs differentially expressed between CRC tumors and matched healthy tissues and found tumor-specific lncRNA variants. Our results show that numerous lncRNAs are differentially expressed and/or accumulate variants in CRC tumors, thereby suggesting a role in CRC progression. More generally, our approach unlocks the study of lncRNAs in FFPE samples, thus enabling the retrospective use of abundant, well documented material available in hospital biobanks.