{'fraction_non_alphanumeric': 0.04021523647691872, 'fraction_numerical': 0.020390824129141887, 'mean_word_length': 5.41671714112659, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 0, 'https://': 4, 'lorem ipsum': 0, 'www.': 4, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 6, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
RiceLncPedia: a comprehensive database of rice long non‐coding RNAs

Long noncoding RNAs (lncRNAs) are referred as RNA molecules with length of at least 200 nucleotides (nt) and usually have low protein-coding potential (Chekanova 2015). In plants, emerging evidences indicate that lncRNAs function as key modulators in development and stress response at the epigenetic, transcriptional and post-transcriptional levels (Chekanova 2015; Lucero et al. 2021). Multiple comprehensive databases have been constructed for human or animals, such as LncBook (Ma et al. 2019).

Long non-coding RNAs (lncRNAs) are referred as RNA molecules with length of at least 200 nucleotides (nt) and usually have low protein-coding potential (Chekanova, 2015). In plants, emerging evidence indicate that lncRNAs function as key modulators in development and stress response at the epigenetic, transcriptional and post-transcriptional levels (Chekanova, 2015;Lucero et al., 2021). Multiple comprehensive databases have been constructed for human or animals, such as LncBook (Ma et al., 2019). Comparatively, the small amount of lncRNAs, the small sample scale or less comprehensive were the main limitations for current plant lncRNA databases, especially for rice, one of the most widely staple food and model crops (Table 1). For example, GREENC provides lncRNAs in 45 plant species, but without their expression profile and genomic features (Paytuvi-Gallart et al., 2019). A rice genome re-annotation database IC4R 2.0 harbours 3215 lncRNA loci, 6259 transcripts but without the relevant multi-omic features of them (Sang et al., 2020). PLncDB V2.0 was updated very recently, containing information for 11565 rice lncRNAs identified from 98 RNA-Seq libraries (Jin et al., 2021). Here, we developed a database, RiceLncPedia (http://3dgenome.hzau.edu.cn/RiceLnc Pedia), to systematically characterize rice lncRNAs with expression profile and multi-omic features to facilitate the understanding and research of rice lncRNAs, including as follows: (i) lncRNA expression profiles in various tissues, development stages and stress treatments; (ii) lncRNA associations with genome variations; (iii) the linkage of lncRNAs with phenotypes; (iv) the overlap information of lncRNAs and transposon elements; and (v) the lncRNAs predicted as miRNA targets or miRNA precursors.
We identified high-confidence rice lncRNAs based on 2312 publicly available RNA-seq libraries following the unified pipeline ( Figure 1a). Briefly, low-quality reads and adapter sequences were trimmed and the clean RNA-seq reads were mapped to rice reference genome Os-Nipponbare-Reference-IRGSP-1.0.41. Transcripts were assembled and merged to acquire comprehensive nonredundant transcripts for subsequent analysis: (i) filter out known protein-coding transcripts, rRNA and tRNA; (ii) transcripts with lengths less than 200nt and with FPKM scores smaller than 0.5 in all samples were discarded successively; (iii) protein-coding potential was predicted using Coding Potential Calculator (CPC2), Plant Long Non-Coding RNA Prediction by Random fOrests (PlncPRO) and PfamScan software. As a consequence, RiceLncPedia accommodates 6925 rice lncRNAs in 5812 gene loci ( Figure 1b). In addition, 40 experimental validated rice lncRNAs were also integrated into RiceLncPedia (Table 1). In the transcript section, each lncRNA transcript is assigned to a unique accession number and shown with molecular features, including location, length, GC content (%), exon number and category. LncRNAs are classified into intergenic lncRNA (lincRNA), intronic lncRNA, sense lncRNA, antisense lncRNA and long non-coding isoforms of known genes according their positions relative to coding genes ( Figure 1b).
For each given lncRNA transcript, a specific page was linked to incorporate sequence, coding score, genome browser view, expression profile, variation, overlapped transposons, small RNA associations,QTLandGWASinformationrelevanttothelncRNA( Figure 1c). Because the specific expression in a specific tissue or under a specific condition indicates the function association (Yanai et al., 2005), the expressionprofilesofanygivenlncRNAcanbevisualizedinbarcharts, for a few represented projects, covering diverse tissues such as leaf, shoot, root, seed, glume and panicle callus, samples from phosphate starvation, salt, cadmium, drought, cold, osmotic and flood stresses as well as samples grown under JA and ABA treatments. An interactive graphic was presented for further visualization of each lncRNA expression in genome browser ( Figure 1c). To explore the relationship of samples, we clustered the represented 339 libraries mentioned above according to all 6925 lncRNAs expression values ( Figure 1d). The resultant clusters were well matched between the indica and japonica groups, basically indicating the reliability of lncRNA expression profiles in RiceLncPedia.
The multi-omics page provides different molecular features for all lncRNAs. The lncRNA expression profiles were provided across all 2312 collected RNA-seq libraries, which are available for download in RiceLncPedia. We calculated the maximum, average and median (FPKM) as well as expression breadth, coefficient of variance (CV), tissue-specificity index and stress-responsive index (s-Value) for each lncRNA transcript, which were harboured in expression section. The specific expressed lncRNAs in given tissues or growth conditions can be screened by selecting a dataset such as 'ABA treatment' and defining a specified range of CV, s-value or expression breadth (Figure 1c, Help section).
Genome variation section contained 50441 SNPs in 4883 lncRNA transcripts with an average of about 10 SNPs per lncRNA transcript (Figure 1e,f), by comparing the position of lncRNAs with SNPs based on 3000 genome projects (http://snpseek.irri.
Precursor of osa-MIR156j (c) (e) Note: a indicates the number include both microRNAs and lncRNAs. b represents the unknown number of rice lncRNAs since the database is not accessible following the link in the reference. c means the experimentally validated lncRNAs were integrated from EVLncRNAs database. All the data item features in RiceLncPedia are shown in bold font. The rice lncRNA database 1493 org/_download.zul). This information will facilitate the research of lncRNA variation association with their structures, expressions, interactions and functions. In plants, some lncRNA-SNPs were implicated to play potential roles in regulating agricultural traits through GWAS or QTL analysis. We, therefore, predicted the lncRNA-SNP-phenotype association if any rice agricultural GWAS tag SNP co-located with a specific lncRNA. Similarly, a lncRNA resided in any rice QTL was also thought of being associated with the relevant trait. The QTL section shows 6684 rice lncRNAs co-located with 513 QTLs, such as 1000 grain weight, drought tolerance and so on, belonging to 25 tissues, development stages or stress tolerance. The GWAS section presents 384 GWAS SNPs residing in 66 lncRNAs transcripts, which refers to 11 agricultural traits (Figure 1e,f). A specified trait-related lncRNAs can be retrieved by selecting the trait in left menu.
A number of lncRNAs were reported to be originated from transposons in plants, and it was demonstrated that TE-associated lncRNAs show tissue-specific transcription and play vital roles in plant abiotic stress responses (Wang et al., 2017). We overall identified 82 transposons overlapped with 448 lncRNA transcripts, involving 474 transposon and lncRNA transcript relations (Figure 1e,f) by comparing the positions of lncRNAs with transposons (Genomic coordinates of Japonica transposon elements, https://www.genome.arizona.edu/cgi-bin/rite/index.cgi). All lncRNAs overlapping with TE were contained in RiceLncPedia as TE-lncRNAs associations in transposon section.
To facilitate the function prediction of rice lncRNAs, we predicted lncRNA targets of microRNAs with psRNATarget software (Dai et al., 2018) and screened the precursors of miRNAs by comparing lncRNAs sequences with rice miRNA precursor (pre-miRNA) hairpin sequences (http://www.mirbase.org/). Blast 2.7 was used with the threshold e-value ≤ 10 À5 , coverage per cent bigger than 90% and -max_hsps as 1. The secondary structures of lncRNAs and relevant pre-miRNAs were built with the localized RNAfold program with default parameters (Gruber et al., 2008). In the end, the Small RNA targets section contains 6060 lncRNAs targets of 713 Osa-miRNAs, building up 64153 lncRNA and Osa-miRNA interactions (Figure 1e). Pre-miRNA section harbours 312 lncRNAs with high homology with 48 pre-miRNAs, involving 554 relations of lncRNAs and miRNAs (Figure 1e,f). The homology information of lncRNAs with rice pre-miRNAs, the optimal secondary structures in both dot-bracket notation and graphical style with the minimum free energy were available to be downloaded in RiceLncPedia.
RiceLncPedia database was constructed using Django as backend Web framework and PostgreSQL (https://www.postgresql.org/ ) as the database engine. JQuery and AJAX (Asynchronous JavaScript and XML) were used to develop Web interfaces. As for the front-end framework, we employed Bootstrap (https://getbootstrap.com) to supply a series of templates to design Web pages with consistent interface components. We adopted the icon in Font Awesome in the RiceLncPedia website (http://www.fontawesome.com.cn/). Data visualization was powered by Pyecharts (https://github.com/pyecha rts/pyecharts) to add interactive diagrams to our website. All the data and methods can be downloaded in download page.
In summary, RiceLncPedia houses a comprehensive collection of rice lncRNAs from the widest samples and with systematic annotation through integrating multi-omics data, covering molecular features, expression profiles, sequences variations, lncRNA-miRNA association, lncRNAs-transposon association and agricultural traits association. All the methods and data are available in help or download page. Future development of RiceLncPedia will refer to regular updates of newly discovered rice lncRNAs, integration of differentially expressed lncRNAs in more diverse tissues and environments, epigenetic features of lncRNAs and the association of lncRNAs with protein-coding genes, experimentally validated lncRNAs and more lncRNA-phenotype associations. We are looking forward to any reasonable suggestions from worldwide scientists, with the aim to provide a continually updated and rich knowledge reservoir of rice lncRNAs and serve as a valuable resource for rice research communities.