EVLncRNAs: a manually curated database for long non-coding RNAs validated by low-throughput experiments.

Long non-coding RNAs (lncRNAs) play important functional roles in various biological processes. Early databases were utilized to deposit all lncRNA candidates produced by high-throughput experimental and/or computational techniques to facilitate classification, assessment and validation. As more lncRNAs are validated by low-throughput experiments, several databases were established for experimentally validated lncRNAs. However, these databases are small in scale (with a few hundreds of lncRNAs only) and specific in their focuses (plants, diseases or interactions). Thus, it is highly desirable to have a comprehensive dataset for experimentally validated lncRNAs as a central repository for all of their structures, functions and phenotypes. Here, we established EVLncRNAs by curating lncRNAs validated by low-throughput experiments (up to 1 May 2016) and integrating specific databases (lncRNAdb, LncRANDisease, Lnc2Cancer and PLNIncRBase) with additional functional and disease-specific information not covered previously. The current version of EVLncRNAs contains 1543 lncRNAs from 77 species that is 2.9 times larger than the current largest database for experimentally validated lncRNAs. Seventy-four percent lncRNA entries are partially or completely new, comparing to all existing experimentally validated databases. The established database allows users to browse, search and download as well as to submit experimentally validated lncRNAs. The database is available at http://biophy.dzu.edu.cn/EVLncRNAs.

INTRODUCTION

The ENCODE project found that three-quarters of the human genome is capable of being transcribed (1) but <3% of which code for proteins (2). The vast majority of the genome transcribe into non-coding RNAs (lncRNAs), many of which are shown to have important biological functions (3). Small non-coding RNAs (<200 nt), such as microRNAs, snoRNAs, siRNAs and piRNAs, have been extensively studied (4,5). The discoveries of two lncRNAs, H19 and Xist (6,7), inspired the search for additional lncRNAs and determination of their biological functions (8). LncRNAs have been shown to be involved in regulating gene expression and other cellular processes and implicated in disease etiology (8–10) including cancers (9,11–13). As more and more lncRNAs are being discovered, a number of databases (14–16) have been established for different purposes.

Early databases were utilized to deposit all lncRNA candidates from high-throughput techniques to facilitate classification, assessment and validation. Examples include LNCipedia for human lncRNA (17), PNRD (Plant non-coding RNA database) (18) and PlantNATsDB for plants (19), and NONCODE for 16 different species (20). Some databases such as NPInter 3.0 (21) contain both manually curated as well as high-throughput experiments. Recognizing potential false positives in high-throughput experiments, five recent databases were built to contain lncRNA validated by low-throughput experiments only. They are lncRNAdb, a manually curated database with 287 functional lncRNAs from eukaryotic organisms (22), LncRNADisease with 1102 entries of experimentally validated lncRNA–disease associations (23), Lnc2Cancer with manually curated 1057 associations between 531 lncRNAs and 86 human cancers (24), lncRInter with 922 lncRNA interaction pairs between 276 lncRNAs and 597 partners in 15 organisms (25), and PLNlncRbase with manually collected 420 plant lncRNAs in 43 plant species after excluding those from high-throughput experiments (26). All of these datasets were well designed for their respective specific purposes [eukaryotic only without legible disease associations (lncRNAdb), disease-focused (Lnc2Cancer and LncRNADisease), interaction-focused (lncRInter) and plant-focused (PLNlncRbase)]. Thus, it is highly desirable to have a central repository for all lncRNAs validated by low-throughput experiments along with their sequence, structure, function and related disease information.

Here, we manually curated all experimentally validated lncRNAs (EVlncRNAs) published prior to 1 May 2016. We compared these lncRNAs with those collected in several specific databases (lncRNAdb, Lnc2Cancer, LncRNADisease and PLNlncRbase) and curated additional functional and disease-specific information not covered previously. The final dataset contains 1543 lncRNAs from 77 species, 2324 associations between 886 lncRNAs and 338 diseases, 793 entries of functional descriptions of 664 lncRNAs, and 1163 interactions between 445 lncRNAs and their interacting partners. More than 70% entries are partially or completely new, comparing to existing specific databases of experimentally validated lncRNAs.

DATABASE CONSTRUCTION

For this database, we only kept lncRNAs confirmed by low-throughput experiments such as quantitative reverse transcription-polymerase chain reaction (qRT-PCR), knock down, western blot, northern blot and luciferase reporter assays. Western blotting validates the interaction between lncRNAs and proteins often after knocking down (27,28). qRT-PCR and northern blot confirm the expression level of lncRNAs (29,30). Luciferase reporter assays examine the promoter activity of lncRNA or their interacting partners and confirm the interaction between lncRNAs and their partners (31,32). We excluded lncRNAs derived from RNA-seq only because these RNAs were transcribed but not necessarily functional. Moreover, RNA-seq data suffers from potential problems of partial transcripts, sequencing biases and low sequencing depth for lowly expressed transcripts, in particular (33). Furthermore, it is still challenging to separate peptide-coding from non-coding RNAs. Given the uncertainty, we have excluded lncRNAs if RNA-seq data is the only evidence. For similar reasons, lncRNAs with only microarray evidence are also excluded.

All EVLncRNAs entries were curated as follows. We first searched the PubMed with keywords matching ‘lncRNA’, ‘long non coding RNA’, ‘long ncRNA’ or ‘long non-code RNA’ along with their plural forms. We then separated those articles covered in previous databases (lncRNAdb, LncRNADisease, Lnc2Cancer and PLNlncRbase) from those completely new ones, which led to 1267 new publications. From these literatures, we manually extracted information on newly validated lncRNAs, their related diseases and interactions, and found 155 lncRNAs not listed in any existing databases for experimentally validated lncRNAs.

Next, the articles referenced as well as not referenced by previous databases were investigated for locating additional information required by our database. We manually examined all entries in lncRNAdb, LncRNADisease, Lnc2Cancer and PLNlncRbase. We excluded 767 lncRNAs collected by PLNlncRbase that was obtained from high-throughput microarray or RNA-seq. For lncRNAs collected by above databases but lacking function annotation, disease association or other information, we manually curated information from literature. In short, we added disease associated to entries in lncRNAdb, function annotations to entries in LncRNADisease, other interaction and disease information to entries in Lnc2Cancer, and interaction information to entries in PLNlncRbase. Only 399 lncRNAs indexed in previous databases were incorporated into our database without finding new information, compared to the fact that 989 entries incorporated were modified with new information. This modification along with 155 novel lncRNAs leads to 1144 (74%, 1144/1543) lncRNAs with new information not covered by existing databases of experimentally validated lncRNAs.

DATABASE CONTENT

Our database EVLncRNAs is intended to be a comprehensive database that includes all species and covers functional and disease-specific roles for all lncRNAs validated by one or more low-throughput experiments. In EVLncRNAs, lncRNAs were named following the HUGO Gene nomenclature Committee (HGNC) (34) with alias included. Sequence and position information of lncRNAs were annotated according to NCBI with links pointing to accessions in the NCBI (35) and Ensembl (36) provided if known. All data in EVLncRNAs were organized using MySQL. The web interface was designed by PHP. The web services were built on Apache Tomcat. The EVLncRNAs database is available at http://biophy.dzu.edu.cn/EVLncRNAs/.

The database contains a total of 1543 lncRNAs from 77 species along with their annotated functions, interaction partners and relevant diseases, if known. Table 1 provides statistics of experimentally validated datasets for comparison. EVlncRNAs has nearly three times more lncRNAs than the largest existing dataset and covers 6 new species, 89 new diseases and 405 new interactions. Each entry contains general information, such as lncRNA official name, alias, species, chromosome, start site, end site, chain, exon number, assembly version, lncRNA class (such as lincRNA, antisense, exonic, intronic etc.) and accession numbers to NCBI and Ensembl. Moreover, we have included 14 peptide-coding lncRNAs.

Table 1.

Comparison of EVLncRNAs with other experimentally validated lncRNA databases

Database	#LncRNAs	#Species	#Diseases	#Disease associations	#Functions (Excl. Interactions)	#Interactions
lncRNAdb	287	71	N/A	N/A	287	307
LncRNADisease	321	14	221	1102	N/A	475
Lnc2Cancer	531	1	86	1057	N/A	N/A
PLNlncRbase	420*	43	N/A	N/A	420*	N/A
lncRInter	276	15	N/A	N/A	N/A	922
EVLncRNAs	1543	77	338	2324	793	1163

*After excluding 767 lncRNA from high-throughput microarray and RNA-seq.

EVLncRNAs contains 2324 associations between 886 lncRNAs and 338 diseases along with information on experimental methods, experimental samples, expression pattern, dysfunction types, detailed description and original articles (PubMed ID and web link). Among the entries, 89 diseases and 346 lncRNAs–disease associations were manually curated from recent literatures that never appeared in other databases. The 346 lncRNAs–disease associations comprise of 76 associations between previously known lncRNAs and diseases, and 270 associations on new lncRNAs or diseases. Experimentally validated functions of lncRNAs were classified according to expression, mutation, interaction, locus and epigenetics (23).

For those lncRNAs that are not associated with diseases, we collected their functional information such as gene regulation, cell differentiation or embryonic development (interaction information is collected separately, see below). There are 166 novel entries in 793 entries from recent literature.

The database also contains 1163 entries of 445 lncRNAs interacting with other molecules, where we included the name of interaction target, level of interaction, type of interaction, description of interaction and original publication (PubMed ID and linked). Here, ‘level of interaction’ represents levels of interaction partner including DNA, RNA and protein, and ‘types of interactions’ were classified into binding, regulation and co-expression; ‘description of interaction’ means descriptions of lncRNA interacting process from the literature. Among these entries, 500 entries were curated from recent literatures. In other words, 43% (500/1163) of interactions entries are new.

USER INTERFACE

The interface of EVLncRNAs allows users to browse, search and download. As shown Figure 1, there are four pull-down menus: Species, Disease, Interaction and Peptide-related in the ‘Browse’ menu (located in labeled position 1). Each menu contains sub-menus. Users can browse relevant entries in EVLncRNAs by clicking any pull-down menus or their sub-menus. Alternatively, users can click navigation bar (located in position 2) on the left to view the corresponding page or entry. Using ‘HOXA11-AS’ as an example (position 3), the page of this entry displays the lncRNA information (position 4), its related disease (position 5), function (position 6) and interaction information (position 7). If an lncRNA involves the function of interacting with other molecules, its interaction function will be described in more details in a following table. In the ‘Search’ page, EVLncRNAs enables users to search by any keywords, such as lncRNA name, alias, disease name, experimental methods, associated components and level of interaction. EVLncRNAs offers ‘fuzzy’ searching capabilities and returns all matching records. Figure 2 shows the search result of ‘HIF1A-AS1’ with specific link to the lncRNA found (labeled in position 1) and original articles (with PMID, labeled in 2) along with its disease association (position 3), function (position 4) and interaction information (position 5) as labeled. EVLncRNAs supports users to download the database data in the ‘Download’ page. In addition, EVLncRNAs provides links to the tools for lncRNA prediction (37,38) in the ‘Tools’ page. EVLncRNAs also permits users to submit novel experimentally validated lncRNAs and related diseases or associated components in the ‘Submit’ page. The submitted record will be included in the next release after manually checked by our curators. A detailed user tutorial is available in the ‘Help’ page.

Figure 1.

Browse page of EVLncRNAs.

Figure 2.

Search result page of EVLncRNAs.

FUTURE EXTENSIONS

It is expected that the number of experimentally validated lncRNAs will continue to increase. EVLncRNAs will be manually curated and updated in a regular basis. We will also gradually add new tools for analysis and prediction of lncRNA, when they are available.

DISCUSSION AND CONCLUSION

LncRNAs have been found to play an essential role in gene regulation and diseases (3,11,12). However, the majority of databases for lncRNAs relied on computational prediction and/or high-throughput experiments with uncertain accuracy. A few experimentally validated lncRNA databases (22–26) have limited information available. The purpose of this work is to build a comprehensive dataset for all lncRNAs that were validated by low-throughput experiments.

The current version of EVLncRNAs (up to 1 May 2016) contains the largest collection of experimentally validated lncRNAs across 77 species. This is close to three times larger than any other experimentally validated datasets. There are 869 lncRNAs for human, 220 for animal, 428 for plant and 26 for microbe (Supplementary Figure S1). The top 20 most studied species on lncRNA were shown in Figure 3. The human and mouse used as a human disease model are top two species with the most identified lncRNAs. The next species is Arabidopsis thaliana, which is a model organism for plant research. About 34 lncRNAs have orthologs in human, mouse and/or rat, which mean the same functional lncRNAs in different species, such as XIST of human and Xist of mouse, and 23 lncRNAs have orthologs in plant species. This should provide important information for evolution studies.

Figure 3.

Top 20 species for the number of experimentally validated lncRNAs.

Disease association is an important component of EVLncRNAs. Diseases were classified into cancers (147 types) and others. Top 15 cancers and top 15 other diseases associated with the highest number of lncRNAs are shown in Supplementary Figures S2 and 3, respectively. That the number of cancer-related lncRNAs is greater than those related with other diseases is likely due to the fact that more researchers are working on cancers because even in LncRNADisease (specifically designed to include lncRNAs relevant to human disease) collections, the majority (627/1070 = 59%) are related to neoplastic conditions. That is, the dominance of cancer-related lncRNAs is not caused by the contribution from the entries from Lnc2Cancer collecting lncRNAs relevant to cancers. Currently, the disease associated with the most number of lncRNAs (>100) is hepatocelluar carcinoma. Many lncRNAs are associated with more than one disease. Those lncRNAs associated with more than 10 diseases are shown in Supplementary Figure S4. H19, for example, is associated with 40 cancers and 24 other diseases.

In addition, EVLncRNAs contains 1163 entries of interaction information between 445 lncRNAs and their associated partners. As shown in Supplementary Figure S5, the highest number of lncRNAs is involved in regulation. This is followed by binding and co-expression.

Click here for additional data file.