Abstract

Circular RNAs (circRNAs) are endogenous non-coding RNAs that form a covalently closed continuous loop, are widely distributed and play important roles in a series of developmental processes. In plants, an increasing number of studies have found that circRNAs can regulate plant metabolism and are involved in plant responses to biotic or abiotic stress. Acting as miRNA decoys is a critical way for circRNAs to perform their functions. Therefore, we developed GreenCircRNA—a database for plant circRNAs acting as miRNA decoys that is dedicated to providing a plant-based platform for detailed exploration of plant circRNAs and their potential decoy functions. This database includes over 210 000 circRNAs from 69 species of plants; the main data sources of circRNAs in this database are NCBI, EMBL-EBI and Phytozome. To investigate the function of circRNAs as competitive endogenous RNAs, the possibility of circRNAs from 38 plants to act as miRNA decoys was predicted. Moreover, we provide basic information for the circRNAs in the database, including their locations, host genes and relative expression levels, as well as full-length sequences, host gene GO (Gene Ontology) numbers and circRNA visualization. GreenCircRNA is the first database for the prediction of circRNAs that act as miRNA decoys and contains the largest number of plant species.

Database URL: http://greencirc.cn

Introduction

Circular RNAs (circRNAs) have been a hot topic in non-coding RNA research in recent years. First discovered in the 1970s (1, 2), circRNAs are characterized by covalently closed-loop structures with neither a 5′ cap nor a 3′ polyadenylated tail and are generated by back splicing (3). Because of their special cyclic structure, circRNAs are insusceptible to degradation by RNA exonuclease or RNase R, which suggests that these RNAs have important functions in vivo (4, 5). With the development of high-throughput RNA sequencing, thousands of circRNAs have been identified in humans, mammals and fungi (6–8). Moreover, many circRNAs have been identified in plants, including Arabidopsis thaliana (9), maize (10), tomato (11), wheat (12) and rice (13). Recent studies found that plant circRNAs play regulatory roles in stress response (14). For instance, Vv-circATS1 responds to cold stress in grape (15), and overexpression of circR5g05160 can improve the resistance of rice to Magnaporthe oryzae (16). MicroRNAs (miRNAs) are a class of non-coding RNAs that are approximately 20–24 nucleotides in length (17) and interact with mRNAs to regulate gene expression (18). Studies have shown that circRNAs are able to sequester and inactivate miRNAs by acting as miRNA decoys (also known as miRNA sponges) (19–21), for which function circRNAs carry a short stretch of sequences that share homology with miRNA binding sites in endogenous mRNA targets. For example, circRNA ciRS-7 has 73 conserved miR-7-binding sites and strongly suppresses miR-7 activity by acting as a miR-7 decoy, affecting expression levels of miR-7 targets (22–24). Currently, many studies are also focusing on the function of circRNA as a miRNA decoy in plants. For instance, 102 circRNAs were found to act as decoys for 24 corresponding miRNAs in tomato (11), 6 circRNAs that respond to dehydration were found to act as decoys for 26 corresponding miRNAs in wheat (12) and 346 circRNAs were found to act as decoys in Zea mays (25). Thus, inferring the function of plant circRNAs based on circRNAs as miRNA decoys is an effective approach.

GreenCircRNA framework.
Figure 1

GreenCircRNA framework.

With the huge progress in animal circRNAs research, some animal circRNA databases have been established, such as CircNet, Circ2Traits and CircRNADb (26–28). In contrast, with the advance of plant circRNAs, there are only two databases of plant circRNAs that have been built, namely, PlantCircNet and PlantcircBase (29, 30), but there is no database for identifying circRNAs as miRNA decoys in plants. In this study, we established Green CircRNA—a database for plant circRNAs that act as miRNA decoys. This database contains 69 plants in total, with 38 plants having relevant miRNA information. Based on the information available, we investigated circRNAs as miRNA decoys and identified mRNAs as miRNA targets in these 38 plants. These findings will facilitate further analysis of the function of circRNAs as competitive endogenous RNAs. GreenCircRNA is a comprehensive plant circRNA database containing 213 494 circRNAs from 69 plant species using 4116 transcriptome sequencing data sets and provides relevant information, full-length sequences and regulatory networks for those circRNAs. We believe that GreenCircRNA is a comprehensive and valuable resource and an important platform for further research on plant circRNAs. The data set used by our database can be downloaded freely. Henceforth, we will continue to supplement the data in GreenCircRNA (http://greencirc.cn).

Aims of the database

CircRNAs play essential roles in regulating plant development and metabolism. The mechanism of formation, function and conservation of circRNAs in plants is the focus of recent circRNA research (6, 8, 31). To gain a deeper understanding of the mechanism and functions of circRNAs as miRNA decoys in plants, we constructed a comprehensive database for plant circRNAs that act as miRNA decoys. This database was created with three main goals: (i) use of publicly available high-throughput transcriptome sequencing data to identify circRNAs in various plants and archive related information for these circRNAs, (ii) analysis of the potential function of circRNAs by identifying circRNAs that act as miRNA decoys and (iii) provision of a user-friendly website with useful web-based tools for the investigation of plant circRNAs.

Web interface.
Figure 2

Web interface.

Materials and Methods Data collection and implementation

The GreenCircRNA database website is built on a `LAMP' (Linux, Apache, MySQL and PHP) open-source architecture. The principal data sources of GreenCircRNA are NCBI (https://www.ncbi.nlm.nih.gov/) (32), EMBL-EBI (https://www.ebi.ac.uk/), SRA database, Phytozome (https://phytozome.jgi.doe.gov/pz/portal.html), miRBase (33) (www.mirbase.org/) and Plant MicroRNA Database (PMRD: http://bioinformatics.cau.edu.cn/PMRD/) (34) (Figure 1). We downloaded plant genome data and annotation information from Phytozome12; transcriptome sequence data, miRNA sequences and mRNA sequences were collected from NCBI, EMBL-EBI, miRBase andPMRD.

Identification of circRNAs and related information

After obtaining all types of data, we identified circRNAs from various plants using CIRI (v2.0.6) and CIRCexplorer2 (v2.3.3) (35–38). In addition, we provide related information regarding circRNAs in the database, including the following:

  • (1) Location of the circRNA in the genome;

  • (2) Junction read number of circRNAs and junction read IDs;

  • (3) Relative expression of the circRNA calculated by [2*junction_reads/(2*junction_reads + non-junction_reads)];

  • (4) Classification of circRNAs (exon, intron and intergenic region);

  • (5) Host gene of circRNAs and GO (Gene Ontology) number of the host gene (excluding intergenic region circRNAs);

  • (6) SRA accession number that generated the circRNAs;

  • (7) Sequence of the circRNA extracted from the start position to the end position (named circ-genome-seq).

Table 1

Statistics of GreenCircRNA database

SpeciesNumber of circRNAsNumber of exon circRNAsNumber of intergenic region circRNAsNumber of intron circRNAsNumber of circRNAs full-length
Amaranthus hypochondriacus15365478721170
Ananas comosus72322042875178
Arabidopsis lyrata2886623103212310
Arabidopsis thaliana10 707525455448992469
Asparagus officinalis1303379767157962
Botryococcus braunii3729023646123
Brachypodium distachyon55522929432227
Brachypodium hybridum288941904104
Brachypodium stacei752249424
Brachypodium sylvaticum25906175261447541
Brassica oleracea capitata21134156781020399
Brassica rapa FPsc1241054605
Capsella grandiflora25171885597351357
Capsella rubella4121257554610001605
Carica papaya3984120080119831089
Chenopodium quinoa37554103183162611
Chlamydomonas reinhardtii28251142258142593
Chromochloris zofingiensis2538915597
Cicer arietinum134579651138253
Citrus clementina1016162440414145
Citrus sinensis18292083551266176
Cucumis sativus5432399211692711472
Daucus carota105349751838327
Dunaliella salina42191679
Eucalyptus grandis3253287290165284
Eutrema salsugineum2481031450119
Fragaria vesca63521341102539861080
Glycine max7938295292540611850
Gossypium hirsutum39882598129199846
Gossypium raimondii652189269194195
Helianthus annuus11 858662847534775141
Hordeum vulgare1162504341317163
Kalanchoe fedtschenkoi11061599074035
Lactuca sativa5540399614291152936
Linum usitatissimum202246249510650
Malus domestica31416535651923219
Manihot esculenta454325591906782190
Marchantia polymorpha7074982045479
Medicago truncatula56682244167017541580
Micromonas pusilla CCMP1545611931022936815180
Micromonas sp.RCC2991518758648
Mimulus guttatus129674753514664
Musa acuminata25041017140087995
Olea europaea860507862477
Oryza sativa17 4796684140593903931
Oryza sativa Kitaake105050947863340
Panicum hallii157856196156254
Panicum virgatum32964412743112408
Phaseolus vulgaris8866624125031225856
Physcomitrella patens43026711053283
Populus deltoides WV9466112252712127
Populus trichocarpa2090229183427248
Porphyra umbilicalis6634472106378
Prunus persica29251976916331420
Ricinus communis28618591264738976
Salix purpurea2310589163784301
Setaria italica2099138866546909
Setaria viridis1514103344041833
Solanum lycopersicum2426122510901111485
Solanum tuberosum6644113634472061880
Sorghum bicolor116554058639436
Spirodela polyrhiza831566238
Theobroma cacao44433169696578664
Trifolium pratense1792759919114390
Triticum aestivum318222496792542330
Vigna unguiculata576746361097343778
Vitis vinifera7412506795113942974
Zea mays12 0353111111878062630
Zostera marina171106632111
SpeciesNumber of circRNAsNumber of exon circRNAsNumber of intergenic region circRNAsNumber of intron circRNAsNumber of circRNAs full-length
Amaranthus hypochondriacus15365478721170
Ananas comosus72322042875178
Arabidopsis lyrata2886623103212310
Arabidopsis thaliana10 707525455448992469
Asparagus officinalis1303379767157962
Botryococcus braunii3729023646123
Brachypodium distachyon55522929432227
Brachypodium hybridum288941904104
Brachypodium stacei752249424
Brachypodium sylvaticum25906175261447541
Brassica oleracea capitata21134156781020399
Brassica rapa FPsc1241054605
Capsella grandiflora25171885597351357
Capsella rubella4121257554610001605
Carica papaya3984120080119831089
Chenopodium quinoa37554103183162611
Chlamydomonas reinhardtii28251142258142593
Chromochloris zofingiensis2538915597
Cicer arietinum134579651138253
Citrus clementina1016162440414145
Citrus sinensis18292083551266176
Cucumis sativus5432399211692711472
Daucus carota105349751838327
Dunaliella salina42191679
Eucalyptus grandis3253287290165284
Eutrema salsugineum2481031450119
Fragaria vesca63521341102539861080
Glycine max7938295292540611850
Gossypium hirsutum39882598129199846
Gossypium raimondii652189269194195
Helianthus annuus11 858662847534775141
Hordeum vulgare1162504341317163
Kalanchoe fedtschenkoi11061599074035
Lactuca sativa5540399614291152936
Linum usitatissimum202246249510650
Malus domestica31416535651923219
Manihot esculenta454325591906782190
Marchantia polymorpha7074982045479
Medicago truncatula56682244167017541580
Micromonas pusilla CCMP1545611931022936815180
Micromonas sp.RCC2991518758648
Mimulus guttatus129674753514664
Musa acuminata25041017140087995
Olea europaea860507862477
Oryza sativa17 4796684140593903931
Oryza sativa Kitaake105050947863340
Panicum hallii157856196156254
Panicum virgatum32964412743112408
Phaseolus vulgaris8866624125031225856
Physcomitrella patens43026711053283
Populus deltoides WV9466112252712127
Populus trichocarpa2090229183427248
Porphyra umbilicalis6634472106378
Prunus persica29251976916331420
Ricinus communis28618591264738976
Salix purpurea2310589163784301
Setaria italica2099138866546909
Setaria viridis1514103344041833
Solanum lycopersicum2426122510901111485
Solanum tuberosum6644113634472061880
Sorghum bicolor116554058639436
Spirodela polyrhiza831566238
Theobroma cacao44433169696578664
Trifolium pratense1792759919114390
Triticum aestivum318222496792542330
Vigna unguiculata576746361097343778
Vitis vinifera7412506795113942974
Zea mays12 0353111111878062630
Zostera marina171106632111
Table 1

Statistics of GreenCircRNA database

SpeciesNumber of circRNAsNumber of exon circRNAsNumber of intergenic region circRNAsNumber of intron circRNAsNumber of circRNAs full-length
Amaranthus hypochondriacus15365478721170
Ananas comosus72322042875178
Arabidopsis lyrata2886623103212310
Arabidopsis thaliana10 707525455448992469
Asparagus officinalis1303379767157962
Botryococcus braunii3729023646123
Brachypodium distachyon55522929432227
Brachypodium hybridum288941904104
Brachypodium stacei752249424
Brachypodium sylvaticum25906175261447541
Brassica oleracea capitata21134156781020399
Brassica rapa FPsc1241054605
Capsella grandiflora25171885597351357
Capsella rubella4121257554610001605
Carica papaya3984120080119831089
Chenopodium quinoa37554103183162611
Chlamydomonas reinhardtii28251142258142593
Chromochloris zofingiensis2538915597
Cicer arietinum134579651138253
Citrus clementina1016162440414145
Citrus sinensis18292083551266176
Cucumis sativus5432399211692711472
Daucus carota105349751838327
Dunaliella salina42191679
Eucalyptus grandis3253287290165284
Eutrema salsugineum2481031450119
Fragaria vesca63521341102539861080
Glycine max7938295292540611850
Gossypium hirsutum39882598129199846
Gossypium raimondii652189269194195
Helianthus annuus11 858662847534775141
Hordeum vulgare1162504341317163
Kalanchoe fedtschenkoi11061599074035
Lactuca sativa5540399614291152936
Linum usitatissimum202246249510650
Malus domestica31416535651923219
Manihot esculenta454325591906782190
Marchantia polymorpha7074982045479
Medicago truncatula56682244167017541580
Micromonas pusilla CCMP1545611931022936815180
Micromonas sp.RCC2991518758648
Mimulus guttatus129674753514664
Musa acuminata25041017140087995
Olea europaea860507862477
Oryza sativa17 4796684140593903931
Oryza sativa Kitaake105050947863340
Panicum hallii157856196156254
Panicum virgatum32964412743112408
Phaseolus vulgaris8866624125031225856
Physcomitrella patens43026711053283
Populus deltoides WV9466112252712127
Populus trichocarpa2090229183427248
Porphyra umbilicalis6634472106378
Prunus persica29251976916331420
Ricinus communis28618591264738976
Salix purpurea2310589163784301
Setaria italica2099138866546909
Setaria viridis1514103344041833
Solanum lycopersicum2426122510901111485
Solanum tuberosum6644113634472061880
Sorghum bicolor116554058639436
Spirodela polyrhiza831566238
Theobroma cacao44433169696578664
Trifolium pratense1792759919114390
Triticum aestivum318222496792542330
Vigna unguiculata576746361097343778
Vitis vinifera7412506795113942974
Zea mays12 0353111111878062630
Zostera marina171106632111
SpeciesNumber of circRNAsNumber of exon circRNAsNumber of intergenic region circRNAsNumber of intron circRNAsNumber of circRNAs full-length
Amaranthus hypochondriacus15365478721170
Ananas comosus72322042875178
Arabidopsis lyrata2886623103212310
Arabidopsis thaliana10 707525455448992469
Asparagus officinalis1303379767157962
Botryococcus braunii3729023646123
Brachypodium distachyon55522929432227
Brachypodium hybridum288941904104
Brachypodium stacei752249424
Brachypodium sylvaticum25906175261447541
Brassica oleracea capitata21134156781020399
Brassica rapa FPsc1241054605
Capsella grandiflora25171885597351357
Capsella rubella4121257554610001605
Carica papaya3984120080119831089
Chenopodium quinoa37554103183162611
Chlamydomonas reinhardtii28251142258142593
Chromochloris zofingiensis2538915597
Cicer arietinum134579651138253
Citrus clementina1016162440414145
Citrus sinensis18292083551266176
Cucumis sativus5432399211692711472
Daucus carota105349751838327
Dunaliella salina42191679
Eucalyptus grandis3253287290165284
Eutrema salsugineum2481031450119
Fragaria vesca63521341102539861080
Glycine max7938295292540611850
Gossypium hirsutum39882598129199846
Gossypium raimondii652189269194195
Helianthus annuus11 858662847534775141
Hordeum vulgare1162504341317163
Kalanchoe fedtschenkoi11061599074035
Lactuca sativa5540399614291152936
Linum usitatissimum202246249510650
Malus domestica31416535651923219
Manihot esculenta454325591906782190
Marchantia polymorpha7074982045479
Medicago truncatula56682244167017541580
Micromonas pusilla CCMP1545611931022936815180
Micromonas sp.RCC2991518758648
Mimulus guttatus129674753514664
Musa acuminata25041017140087995
Olea europaea860507862477
Oryza sativa17 4796684140593903931
Oryza sativa Kitaake105050947863340
Panicum hallii157856196156254
Panicum virgatum32964412743112408
Phaseolus vulgaris8866624125031225856
Physcomitrella patens43026711053283
Populus deltoides WV9466112252712127
Populus trichocarpa2090229183427248
Porphyra umbilicalis6634472106378
Prunus persica29251976916331420
Ricinus communis28618591264738976
Salix purpurea2310589163784301
Setaria italica2099138866546909
Setaria viridis1514103344041833
Solanum lycopersicum2426122510901111485
Solanum tuberosum6644113634472061880
Sorghum bicolor116554058639436
Spirodela polyrhiza831566238
Theobroma cacao44433169696578664
Trifolium pratense1792759919114390
Triticum aestivum318222496792542330
Vigna unguiculata576746361097343778
Vitis vinifera7412506795113942974
Zea mays12 0353111111878062630
Zostera marina171106632111

CircRNAs as miRNA decoys

Functional research on circRNAs is the main challenge in the study of plant circRNAs. To infer the potential functions of circRNAs, we established a circRNA–miRNA–mRNA network in which circRNAs act as miRNA decoys and mRNAs act as miRNA targets. First, we downloaded miRNA sequences from miRBase and PMRD and mRNA sequences from Phytozome12 and extracted circRNA sequences by using an in-house Perl script. Next, to establish the relationship between miRNAs and other RNAs, we used RNAplex with default parameters to predict RNA–RNA hybridization sites (39, 40). We then identified miRNA decoys (circRNAs) and targets (mRNA) following a method proposed in our previous report (41–43, 25). The criteria used to define a miRNA decoy were as follows: no more than six mismatched or inserted bases present between the 9th and 20th nucleotides of the miRNA 5′ end, perfect matching of the second to eighth bases of the miRNA 5′ end sequence and no more than four mismatches or indels in other regions. The criteria used to define a miRNA target were as follows: at most, one mismatch or indel was allowed between the 9th and 12th positions of the 5′ end of miRNA sequences, the total number of bulges or mismatches in the other regions was not allowed to exceed 4 nt and no continuous mismatches were allowed. Finally, a picture of the circRNA–miRNA–mRNA network for a species was generated using Cytoscape (v3.7.2) (44). The influence of circRNAs on other mRNAs via miRNAs can be assessed according to the circRNA–miRNA–mRNA network, and the potential functions of the circRNAs can be inferred.

Full-length sequences of circRNAs

The full-length sequences of circRNAs are important for subsequent analysis of the internal structural features and functions of circRNAs, which can help in evaluating the translation potential of circRNAs. However, most circRNA databases do not provide the full-length sequences of circRNAs. For our database, we used downstream programs of CIRI, named CIRI-full, to assemble the full-length sequences of circRNAs (45). After removing redundant results with the in-house Perl script, we obtained final circRNA full-length sequences.

CircRNA visualization

CircRNAs form covalently closed-loop structures with neither 5′-3′ polarities nor polyadenylated tails. As providing the position of circRNAs in the genome does not lead to an intuitive understanding of its structure for users, a schematic of each circRNA is available in the database. The schematic consists of two parts: one is a line of the host gene of the circRNA, with the exons and introns in the gene labeled; the other is a circle of the circular structure of circRNA, which is marked with corresponding colors that are the same as the position on the gene. Thus, the position and structure of the circRNA can be observed visually, which makes the information easy to understand.

(A) Distribution of three types of circRNAs in maize. (B) Visualization of the circRNA Zmays_10:10199025|10 199 243.
Figure 3

(A) Distribution of three types of circRNAs in maize. (B) Visualization of the circRNA Zmays_10:10199025|10 199 243.

(A) Relative expression of circRNAs in Zea mays. (B) Length distribution of full-length circRNAs in Zea mays.
Figure 4

(A) Relative expression of circRNAs in Zea mays. (B) Length distribution of full-length circRNAs in Zea mays.

(A) CircRNA–miRNA–mRNA sub-network, taking circRNA Zmays_3:225352281|225 387 878 as an example. Blue nodes, miRNA. Pink node, circRNA as a miRNA decoy. Green nodes, mRNAs as miRNA targets. Gray edges, correlations. (B) Sequence list for Figure 5A, taking circRNA Zmays_3: 225352281|225 387 878 as an example. The first column is a circRNA that acts as a miRNA decoy, the second column shows miRNAs and the third column contains mRNAs that act as miRNA targets. This circRNA may act as a decoy for three miRNAs: zma-miR171k-3p/zma-miR171h-3p, zma-miR395a-5p and zma-miR395e-5p/zma-miR395h-5p/zma-miR395j-5p/zma-miR395p-5p.
Figure 5

(A) CircRNA–miRNA–mRNA sub-network, taking circRNA Zmays_3:225352281|225 387 878 as an example. Blue nodes, miRNA. Pink node, circRNA as a miRNA decoy. Green nodes, mRNAs as miRNA targets. Gray edges, correlations. (B) Sequence list for Figure 5A, taking circRNA Zmays_3: 225352281|225 387 878 as an example. The first column is a circRNA that acts as a miRNA decoy, the second column shows miRNAs and the third column contains mRNAs that act as miRNA targets. This circRNA may act as a decoy for three miRNAs: zma-miR171k-3p/zma-miR171h-3p, zma-miR395a-5p and zma-miR395e-5p/zma-miR395h-5p/zma-miR395j-5p/zma-miR395p-5p.

Results

Usage and access

GreenCircRNA mainly includes the following modules: Home, Species, Search, Download and FAQ (Figure 2). Users can browse, search and download circRNA information through the web interface of GreenCircRNA.

Species

Basic information of many circRNAs belonging to 69 plant species is included in this module. There is an individual interface of each plant containing detailed circRNA associated information, including a list of all circRNAs in this plant, a histogram showing the relative expression levels of all circRNAs, a length distribution histogram of the circRNA full-length sequences and a circRNA–miRNA–mRNA network illustrating the potential relationships of these RNAs. Furthermore, there is a separate page for every circRNA that displays detailed information for the circRNA, including the location of the circRNA in the genome, the relative expression level, the full-length sequence, the host gene and annotation information. Users can obtain all information for a given circRNA on this page, and this page also provides a circRNA network that shows circRNAs as miRNA decoys in a tabular and graphical manner.

Search

The module enables users to search for circRNAs by host gene, miRNA ID, circRNA ID and SRA ID. In addition, a subset search is available in this module, and users can search circRNAs by a series of criteria, such as plant species name, chromosome and relative expression level.

Download

Related information for circRNAs, including basic information, circ-genome-seq, full-length sequences and network, can all be downloaded for free in CSV format or fasta format from the `download’ module.

Data summary

The GreenCircRNA database covers 69 plant species, which is the highest number of plant species in all plant circRNA databases to date. We downloaded a total of 4116 transcriptome sequencing data sets from SRA and EMBL-EBI for circRNA identification and eventually obtained 213 494 circRNAs. These circRNAs were classified into three categories: exon circRNAs, intron circRNAs and intergenic region circRNAs. Among all the identified circRNAs, 95 010 (44.50%) belong to the exon category, 65 175 (30.53%) belong to the intergenic region category and 53 309 (24.97%) belong to the intron category. We then extracted the sequences of all circRNAs from genome sequences (circ-genome-seq). Furthermore, we assembled 68 237 full-length sequences of circRNAs using CIRI-full and predicted the conditions that circRNAs act as miRNA decoys for 38 plant species (Table 1).

A case study

Taking maize as an example, we downloaded 181 transcriptome sequencing data sets from the SRA and EMBL-EBI databases for circRNA identification and obtained 12 035 circRNAs in total, including 3111 (25.85%) belonging to exons, 7806 (64.86%) belonging to introns and 1118 (9.29%) belonging to intergenic regions (Figure 3A). By contrast, the results of software CIRI predicted more exon circRNAs, and the result of software CIRCexplorer predicted more intron circRNAs. In addition, 2630 full-length sequences of circRNAs in maize were assembled by CIRI-full; the relative expression levels and length distributions of the circRNA full-length sequences are displayed in histograms (Figure 4A and B). Moreover, an individual page shows detailed information for each circRNA, including the location in the genome, the relative expression level and visualization. For instance, the exon circRNA Zmays_10:10199025|10 199 243 is located in the middle of the third exon of the gene `GRMZM2G046284.v6a’ (Figure 3B). We analyzed circRNAs that may act as miRNA decoys in maize and showed their relationships by the circRNA–miRNA–mRNA network and a table that includes these RNAs involved in the network on the single circRNA page. For example, circRNA Zmays_3:225352281|225 387 878 may act as decoys for three miRNAs: zma-miR171k-3p/zma-miR171h-3p, zma-miR395a-5p and zma-miR395e-5p/zma-miR395h-5p/zma-miR395j-5p/zma-miR395p-5p (Figure 5).

Discussion and future prospects

Increasing evidence has proven that circRNAs play important roles in various biological processes, but few studies have examined circRNAs in plants. Furthermore, specific circRNA data for analyzing the sequence and structure of circRNAs are not available for most plant species, and the functions and mechanisms of most circRNAs are unclear. Although two databases of plant circRNAs have been built, these databases cover relatively fewer plant species. In this study, we developed GreenCircRNA, a comprehensive database of plant circRNAs that includes circRNAs in 69 plants, and analyzed the potential decoy function of these circRNAs in 38 plants. Users can freely search and download information related to circRNAs. We hope that as a platform, GreenCircRNA will help researchers to study the basic properties and characteristics of plant circRNAs and will be useful for further research on the internal structure, translation function and mechanism of circRNAs. This database can still be improved, and we will continuously identify and collect more circRNAs and update and improve GreenCircRNA to provide accurate information regarding plant circRNAs and circRNAs as miRNA decoys.

Funding

National Science Foundation of China (31770333, 31370329 and 11631012); Program for New Century Excellent Talents in University (NCET-12-0896); Fundamental Research Funds for the Central Universities (GK201403004).

Conflict of interest. The authors declare no conflicts of interest.

References

1.

Sanger
,
H.L.
,
Klotz
,
G.
,
Riesner
,
D.
et al.  (
1976
)
Viroids are single-stranded covalently closed circular RNA molecules existing as highly base-paired rod-like structures
.
Proc. Natl. Acad. Sci. U. S. A.
,
73
,
3852
3856
.

2.

Houseley
,
J.M.
,
Garcia-Casado
,
Z.
,
Pascual
,
M.
et al.  (
2006
)
Noncanonical RNAs from transcripts of the Drosophila muscleblind gene
.
J. Hered.
,
97
,
253
260
.

3.

Chen
,
L.L.
and
Yang
,
L.
(
2015
)
Regulation of circRNA biogenesis
.
RNA Biol.
,
12
,
381
388
.

4.

Suzuki
,
H.
and
Tsukahara
,
T.
(
2014
)
A view of pre-mRNA splicing from RNase R resistant RNAs
.
Int. J. Mol. Sci.
,
15
,
9331
9342
.

5.

Xiao
,
M.S.
,
Ai
,
Y.
and
Wilusz
,
J.E.
(
2020
)
Biogenesis and functions of circular RNAs come into focus
.
Trends Cell Biol.
,
30
,
226
240
.

6.

Andranik
,
I.S.
Andranik
,
I.
Sebastian
,
M.
Emanuel
,
W.
et al.  (
2015
)
Analysis of intron sequences reveals hallmarks of circular RNA biogenesis in animals
.
Cell Rep.
,
10
,
170
177
.

7.

Salzman
,
J.
,
Chen
,
R.E.
,
Olsen
,
M.N.
et al.  (
2013
)
Cell-type specific features of circular RNA expression
.
PLoS Genet.
,
9
,
e1003777
.

8.

Jeck
,
W.R.
,
Sorrentino
,
J.A.
,
Wang
,
K.
et al.  (
2013
)
Circular RNAs are abundant, conserved, and associated with ALU repeats
.
RNA
,
19
,
141
157
.

9.

Wang
,
P.L.
,
Bao
,
Y.
,
Yee
,
M.C.
et al.  (
2014
)
Circular RNA is expressed across the eukaryotic tree of life
.
PLoS One
,
9
,
e90859
.

10.

Chen
,
L.
,
Zhang
,
P.
,
Fan
,
Y.
et al.  (
2018
)
Circular RNAs mediated by transposons are associated with transcriptomic and phenotypic variation in maize
.
New Phytol.
,
217
,
1292
1306
.

11.

Zuo
,
J.
,
Wang
,
Q.
,
Zhu
,
B.
et al.  (
2016
)
Deciphering the roles of circRNAs on chilling injury in tomato
.
Biochem. Biophys. Res. Commun.
,
479
,
132
138
.

12.

Wang
,
Y.
,
Yang
,
M.
,
Wei
,
S.
et al.  (
2016
)
Identification of circular RNAs and their targets in leaves of Triticum aestivum L. under dehydration stress
.
Front. Plant Sci.
,
7
,
2024
doi: .

13.

Lu
,
T.
,
Cui
,
L.
,
Zhou
,
Y.
et al.  (
2015
)
Transcriptome-wide investigation of circular RNAs in rice
.
RNA
,
21
,
2076
2087
.

14.

Chu
,
Q.J.
,
Shen
,
E.
,
Ye
,
C.Y.
et al.  (
2018
)
Emerging roles of plant circular RNAs
.
J. Plant Cell Dev.
,
1
,
1
14

15.

Gao
,
Z.
,
Li
,
J.
,
Luo
,
M.
et al.  (
2019
)
Characterization and cloning of grape circular RNAs identified the cold resistance-related Vv-circATS1
Plant Physiol.
,
180
,
966
985
.

16.

Fan
,
J.
,
Quan
,
W.
,
Li
,
G.B.
et al.  (
2020
)
circRNAs are involved in the rice-Magnaporthe oryzae interaction
Plant Physiol.
,
182
,
272
286
.

17.

Pasquinelli
,
A.E.
(
2002
)
MicroRNAs: deviants no longer
.
Trends Genet.
,
18
,
171
173
.

18.

Fabian
,
M.R.
and
Sonenberg
,
N.
(
2012
)
The mechanics of miRNA-mediated gene silencing: a look under the hood of miRISC
.
Nat. Struct. Mol. Biol.
,
19
,
586
593
.

19.

Memczak
,
S.
,
Jens
,
M.
,
Elefsinioti
,
A.
et al.  (
2013
)
Circular RNAs are a large class of animal RNAs with regulatory potency
.
Nature
,
495
,
333
338
.

20.

Franco-Zorrilla
,
J.M.
,
Valli
,
A.
,
Todesco
,
M.
et al.  (
2007
)
Target mimicry provides a new mechanism for regulation of microRNA activity
.
Nat. Genet.
,
39
,
1033
1037
.

21.

Hansen
,
T.B.
,
Jensen
,
T.I.
,
Clausen
,
B.H.
et al.  (
2013
)
Natural RNA circles function as efficient microRNA sponges
.
Nature
,
495
,
384
388
.

22.

Peng
,
L.
,
Yuan
,
X.Q.
and
Li
,
G.C.
(
2015
)
The emerging landscape of circular RNA ciRS-7 in cancer (review)
.
Oncol. Rep.
,
33
,
2669
2674
.

23.

Weng
,
W.
,
Wei
,
Q.
,
Toden
,
S.
et al.  (
2017
)
Circular RNA ciRS-7-a promising prognostic biomarker and a potential therapeutic target in colorectal cancer
.
Clin. Cancer Res.
,
23
,
2541
2016
.

24.

Yu
,
L.
,
Gong
,
X.
,
Sun
,
L.
et al.  (
2016
)
The circular RNA Cdr1as act as an oncogene in hepatocellular carcinoma through targeting miR-7 expression
.
PLoS One
,
11
,
e0158347
.

25.

Tang
,
B.
,
Hao
,
Z.
,
Zhu
,
Y.
et al.  (
2018
)
Genome-wide identification and functional analysis of circRNAs in Zea mays.
PLoS One
,
13
,
e0202375
.

26.

Liu
,
Y.C.
,
Li
,
J.R.
,
Sun
,
C.H.
et al.  (
2016
)
CircNet: a database of circular RNAs derived from transcriptome sequencing data
.
Nucleic Acids Res.
,
44
,
D209
D215
.

27.

Ghosal
,
S.
,
Das
,
S.
,
Sen
,
R.
et al.  (
2013
)
Circ2Traits: a comprehensive database for circular RNA potentially associated with disease and traits
.
Front. Genet.
,
4
,
283
.

28.

Chen
,
X.
,
Han
,
P.
,
Zhou
,
T.
et al.  (
2016
)
circRNADb: a comprehensive database for human circular RNAs with protein-coding annotations
.
Sci. Rep.
,
6
,
34985
.

29.

Chu
,
Q.
,
Zhang
,
X.
,
Zhu
,
X.
et al.  (
2017
)
PlantcircBase: a database for plant circular RNAs
.
Mol. Plant
,
10
,
1126
1128
.

30.

Zhang
,
P.
,
Meng
,
X.
,
Chen
,
H.
et al.  (
2017
)
PlantCircNet: a database for plant circRNA–miRNA–mRNA regulatory networks
.
Database (Oxford)
,
1
. doi: .

31.

Li
,
Q.F.
,
Zhang
,
Y.C.
,
Chen
,
Y.Q.
et al.  (
2017
)
Circular RNAs roll into the regulatory network of plants
.
Biochem. Biophys. Res. Commun.
,
488
,
382
386
.

32.

Pruitt
,
K.D.
,
Tatusova
,
T.
,
Brown
,
G.R.
et al.  (
2012
)
NCBI Reference Sequences (RefSeq): current status, new features and genome annotation policy
.
Nucleic Acids Res.
,
40
,
D130
D135
.

33.

Griffiths-Jones
,
S.
(
2006
)
miRBase: the microRNA sequence database
.
Methods Mol. Biol.
,
342
,
129
138
.

34.

Zhang
,
Z.
,
Yu
,
J.
,
Li
,
D.
et al.  (
2010
)
PMRD: plant microRNA database
.
Nucleic Acids Res.
,
38
,
D806
D813
.

35.

Gao
,
Y.
,
Wang
,
J.
and
Zhao
,
F.
(
2015
)
CIRI: an efficient and unbiased algorithm for de novo circular RNA identification
.
Genome Biol.
,
16
,
4
. doi: .

36.

Gao
,
Y.
and
Zhao
,
F.
(
2018
)
Computational strategies for exploring circular RNAs
.
Trends Genet.
,
34
,
389
400
.

37.

Xiangxiang
,
Z.
,
Wei
,
L.
,
Maozu
,
G.
et al.  (
2017
)
A comprehensive overview and evaluation of circular RNA detection tools
.
PLoS Comput. Biol.
,
13
,
e1005420
.

38.

Zhang
,
X.O.
,
Wang
,
H.B.
,
Zhang
,
Y.
et al.  (
2014
)
Complementary sequence-mediated exon circularization
.
Cell
,
159
,
134
147
.

39.

Tafer
,
H.
and
Hofacker
,
I.L.
(
2008
)
RNAplex: a fast tool for RNA-RNA interaction search
.
Bioinformatics
,
24
,
2657
2663
.

40.

Lorenz
,
R.
,
Bernhart
,
S.H.
,
Christian
,
H.S.
et al.  (
2011
)
ViennaRNA package 2.0
.
Algorithms Mol. Biol.
,
6
,
26
. doi: .

41.

Rajewsky
,
N.
Jurga
,
S.
and
Barciszewski
,
J.
(
2017
) Genome-wide function analysis of lincRNAs as miRNA targets or decoys in plant.
Plant Epigenetics
.
RNA Technologies (Chapter 8)
,
149
162
. doi: .

42.

Hua-Jun
,
W.
,
Zhi-Min
,
W.
,
Meng
,
W.
et al.  (
2013
)
Widespread long noncoding RNAs as endogenous target mimics for microRNAs in plants
.
Plant Physiol.
,
161
,
1875
1884
.

43.

Fan
,
C.
,
Hao
,
Z.
,
Yan
,
J.
and
Li
,
G.
(
2015
)
Genome-wide identification and functional analysis of lincRNAs acting as miRNA targets or decoys in maize
.
BMC Genomics
,
16
,
793
.

44.

Shannon
,
P.
,
Markiel
,
A.
,
Ozier
,
O.
et al.  (
2003
)
Cytoscape: a software environment for integrated models of biomolecular interaction networks
.
Genome Res.
,
13
,
2498
2504
.

45.

Zheng
,
Y.
,
Ji
,
P.
,
Chen
,
S.
et al.  (
2019
)
Reconstruction of full-length circular RNAs enables isoform-level quantification
.
Genome Med.
,
11
,
2
. doi: .

Author notes

These authors contributed equally to this work.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted reuse, distribution, and reproduction in any medium, provided the original work is properly cited.