- Research article
- Open Access
Identification of cucumber circular RNAs responsive to salt stress
BMC Plant Biology volume 19, Article number: 164 (2019)
Circular RNAs (circRNAs) are 3′-5′ head-to-tail covalently closed non-coding RNA that have been proved to play essential roles in many cellular and developmental processes. However, no information relate to cucumber circRNAs is available currently, especially under salt stress condition.
In this study, we sequenced circRNAs in cucumber and a total of 2787 were identified, with 1934 in root and 44 in leaf being differentially regulated under salt stress. Characteristics analysis of these circRNAs revealed following features: most of them are exon circRNAs (79.51%) and they prefer to arise from middle exon(s) of parent genes (2035/2516); moreover, most of circularization events (88.3%) use non-canonical-GT/AG splicing signals; last but not least, pairing-driven circularization is not the major way to generate cucumber circRNAs since very few circRNAs (18) contain sufficient flanking complementary sequences. Annotation and enrichment analysis of both parental genes and target mRNAs were launched to uncover the functions of differentially expressed circRNAs induced by salt stress. The results showed that circRNAs may be paly roles in salt stress response by mediating transcription, signal transcription, cell cycle, metabolism adaptation, and ion homeostasis related pathways. Moreover, circRNAs may function to regulate proline metabolisms through regulating associated biosynthesis and degradation genes.
The present study identified large number of cucumber circRNAs and function annotation revealed their possible biological roles in response to salt stress. Our findings will lay a solid foundation for further structure and function studies of cucumber circRNAs.
Circular RNAs (circRNAs) have recently emerged as a new form of non-coding RNA in a wide variety of organisms . Different from the traditional linear RNAs, they are special endogenous RNA molecules that are predominantly formed by 3′-5′ head-to-tail ligation in a backsplicing reaction. Although observed for decades in eukaryotic cells, circRNAs had been perceived as splicing errors in the past. Recently, knowledge of circRNAs has changed substantially. Thanks to the technological breakthroughs in high-throughput sequencing and high-efficiency bioinformatics analyzing strategies, a vast number of circRNAs have been identified in mice, humans, Archaea, and Caenorhabditis elegans, etc. . Characterization analyses indicate that circRNAs can arise from exons (exonic circRNA), introns (intronic circRNA), and intergenic region and their expression patterns are often cell, tissue and developmental stage-specific . Furthermore, evidence showed that circRNAs are much more stable than traditional linear RNAs and therefore might be involved in biological processes through a distinct way .
In recent years, the biological importance of circRNAs in transcriptional and post-transcriptional regulation have been increasingly illuminated. One of the most studied functions of circRNAs in animals is that they can function as microRNA (miRNA) sponge. For example, Zheng, et al.  reported that human circHIPK3 has 18 potential binding sites and can sponge to 9 miRNAs. Wang, et al.  attested that a heart-related circRNA (HRCR) acts as a ceRNA to sponge miR-223 and inhibit its activity. Another remarkable role of circRNAs is to regulate parental gene expression. For example, in human, some circRNAs have been reported to promote the expression of their parental genes in cis via specific RNA-RNA interactions . Besides, some circRNAs are even able to regulate cell cycle progression or ageing by directly interacting with protein [7, 8]. Despite the relevant functional roles of several circRNAs have been revealed, the majority of them are still mysterious to us, especially their functions in biological processes and associated regulating pathways.
In contrast to animals, to date, plant circRNAs have not been comprehensively identified. Wang, et al.  firstly carried out circRNA research in plant and found clear evidence for RNA circles in Arabidopsis. Meanwhile, Ye, et al.  identified 12,037 and 6012 circRNAs in rice (Oryza sativa) and Arabidopsis through genome-wide level inspecting, which confirmed the widespread of circRNAs in both monocot and dicot plants. Later, circRNAs have been identified in tomato [1, 11], potato , kiwifruit , wheat , soybean , tea , barely , and maize . It is worth noting that although circRNA-miRNA-mRNA network has been constructed in many of these studies in plant, unlike in animals, miRNA sponge roles of circRNAs in plants has not been proven yet. Until now, functional studies on circRNAs are limited, especially in plants. Some studies suggested that circRNAs can regulate gene expression at the transcriptional and/or posttranscriptional levels by interacting with miRNAs [6, 14]. A transgenic study revealed that overexpression of a circRNA construct could reduce the expression level of its parental gene in transgenic rice, suggesting that circRNAs can act as a negative regulator for their parental genes. Similarly, in tomato, Tan, et al.  overexpressed a circRNA derived from a key gene in carotenoid biosynthesis (Phytoene Synthase 1), resulting in a decreasement of Phytoene Synthase 1 mRNA abundance, as well as the lycopene and β-carotene accumulation. More recently, Cheng, et al.  reported that over-expression of a lariat-derived circRNA, generated from the first intron of AT5G37720, alters expression of approximately 800 genes and influence development of Arabidopsis. Moreover, an increasing number of studies in plants has also proposed that the circRNAs play important roles in regulating the response to environmental stresses. In rice, 27 exonic circRNAs were found to be differentially expressed under phosphate-sufficient and -starvation conditions, suggesting that circRNAs play roles in stress-specific biological processes in plants . Similarly, in wheat, kiwifruit, and potato, circRNAs have been found to be specific responded to plant dehydration and bacterial stresses [12,13,14]. In Arabidopsis, Pan, et al.  revealed that heat stress induced more circRNAs than in normal condition, and these circRNAs may participate in plant response to heat stress through the circRNA-mediated ceRNA networks. These observations confirmed that circRNAs involved in a wide range of biological processes in plants, including growth and development, abiotic and biotic stress response.
Several characteristics of animal circRNAs, such as conservation, tissue-specific expression profiles, and longer flanking introns, are shared by plant circRNAs. Some distinct features of plant circRNAs have also been reported. For instance, plant seems to possess much less repetitive elements and reverse complementary sequences in their flanking introns . Recently, Ye, et al.  found that large number of rice circRNAs (92.7%) are flanked by diverse non-GT/AG splicing signals while most human exonic circRNAs are flanked by canonical GT/AG splicing signals. The first achievement towards understanding the molecular mechanisms of plant circRNAs has been obtained in Arabidopsis. Conn, et al.  demonstrated that Arabidopsis circRNA generated from the exon 6 of the SEPALLATA3 (SEP3) genes regulates the expression of its parental gene (SEP3) by binding strongly to its cognate DNA locus (forming an RNA:DNA hybrid or R-loop), revealing the ability of specific circRNAs to skew splicing preference and favour the cognate alternative splicing (AS) mRNA variant. Furthermore, in maize, Chen et al.  hypothesized that transposons may be involved in the formation of circRNAs and they may further modulate phenotypic variation. Taken together, these studies provide novel insight into the biological function and generation mechanism of circRNAs in plant. Genome-wide discovery and characterization of circRNAs should be performed in more plants, which could provide a basis for further studies of the evolution, biogenesis, and function of circRNAs.
Salt stress is a major abiotic stress that limits plant growth and development. Currently, more than 20% of the world’s agricultural irrigated land is adversely affected by levels of salt that could markedly decrease the harvested agriculture yield . Cucumber is an important controlled-environment horticultural crop and is fairly sensitive to salinity [24, 25]. Previously, non-coding RNAs, such as microRNA, have been shown to play critical roles in regulating gene expression in response to abiotic stress including salt stress in plant [26,27,28]. However, to our knowledge, there is still no evidence showing the role of circRNAs, a new type of non-coding RNAs, in the salt stress resistance of cucumber. Besides, although genome-wide identification of circRNAs has been performed in multiple plant species, the characteristics and functions of circRNAs are not clear in cucumber. Our study firstly focused on identification and characterization of circRNAs in cucumber, and uncovered a potential functional role of circRNAs in the regulation of salt stress response.
Plant materials and treatments
Cucumber (Cucumis sativus L. ‘JinYou 1’, Xintiandi Co., Yangling, Shannxi, China) seeds were rinsed thoroughly in distilled water and germinated on moist gauze in an incubator at 28 °C for 2 days. The germinated seeds were sown in quartz sands in a greenhouse with an average day/night (12 h/12 h) temperature of 28 °C/18 °C. The seedlings were transferred to 15-L plastic containers filled with 1/4 strength of Hoagland nutrient solution at two-leaf stage, and each container had 10 plants. Three days later, the strength of Hoagland solution was increased to 1/2. Seven days after transplanting, salt treatment was started by adding sodium chloride (NaCl 75 mM) to the nutrient solution. The pH of nutrient solution was adjusted to 6.0 using 0.2 M H2SO4 or 1 M KOH.
Growth parameters and root morphology traits
Dry weight was determined by collecting shoots and roots after 2, 4, 6, 8, 10, 12, 14, 16, 18 days’ of treatment, washing with distilled water, dried in an oven at 75 °C for 72 h, and then weighted. Roots were collected and determined by scanning roots with the root scanner platform (Model MRS-9600TFU2L, Microtek, Shanghai, China). The scanned images were analyzed with a root image analysis software (MRS-9600TFU2L, Zhongjing Technology Co. Ltd., Hangzhou, China). Each treatment included three replications and was successively performed in 2 independent experiments.
After 3 days of treatment, the second fully expanded leaves and roots were collected from 12 independent plants of the CK and NaCl-treated seedlings with every 6 of them being mixed as one biological replicate for each treatment and then frozen in liquid nitrogen immediately and stored at − 80 °C. All the treatments had 2 independent experiment replicates. The total RNA extraction and quality confirmation were performed according to the methods reported by Yin, et al. . Then, following the manufacturer’s procedure, 2 μg of total RNA from each sample was treated with Ribo-Zero rRNA Removal Kit (for plant, Epicentre, Madison, WI, USA) to obtain ribosomal RNA (rRNA)-depleted RNA, and was further treated with RNase R (Epicentre, Madison, WI, USA) to eliminate linear RNAs. The remaining RNAs were used to generate circRNA-seq libraries according to Yin, et al. . Then the libraries were sequenced on Illumina HiSeq™ X-ten platform at Biomarker Technologies Co., Ltd. (Beijing, China), and 150 bp pair-end reads were generated. The raw sequencing data were deposited in the US National Center for Biotechnology Information (NCBI) Sequence Read Archive (https://www.ncbi.nlm.nih.gov/sra) with a Bioproject ID: PRJNA449130.
Identification of circRNAs
These low quality reads, including that (1) unknown (N) bases are more than 5%, (2) contain adaptor sequences, and (3) contain low quality bases (Q ≤ 20) more than 50%, were removed. The remaining high-quality reads were aligned against the previously reported C. sativus reference genome v2 using BWA (v0.7.15, mem -T 19). The output SAM files were then inspected by CIRI2 (v2.05), the considerable reliability and sensitive tool, to identify circRNAs . Unanchored scaffolds in cucumber genome were merged as ChrUn (detailed information can be found in Additional file 1: Table S1). CircRNA expression levels were normalized to be SRPBM (SRPBM = number of circular reads/number of mapped reads (units in billion)/read length).
Annotation and GO enrichment of circRNAs
Databases including Nr (NCBI non-redundant protein sequences), GO (Gene Ontology), KEGG (Kyoto Encyclopedia of Genes and Genomes), KOG/COG (Clusters of Orthologous Groups of proteins), Swiss-Prot (A manually annotated and reviewed protein sequence database, https://www.uniprot.org/), and Pfam (Protein family) were used to perform gene function annotation. The identification of differentially expressed circRNAs were performed using R package EBSeq. The cut-off criterion were set as |log2(Fold change)| > 1 and false discovery rate (FDR) adjusted p-value ≤0.05. GO annotations of the parent genes of differentially expressed circRNAs were collected and used to perform GO enrichment by topGO with default parameters. Veen and Weblogo diagrams were draw using online tools (http://bioinformatics.psb.ugent.be/webtools/Venn/, http://weblogo.berkeley.edu/logo.cgi). Circos and Violin diagrams were draw using correspondence R packages.
Prediction for circRNA-miRNA-mRNA relationships
To construct the circRNA-miRNA-mRNA network, we collected the mature sequences of cucumber miRNAs from the previously reported publications [30,31,32,33,34,35] (Additional file 2: Table S2). Then the miRNAs identical alignment regions in circRNAs were predicted by TargetFinder using default parameters . Meanwhile, cucumber coding genes (mRNAs) targeted by collected miRNAs were predicted by searching the selected library [Cucumis sativus (cucumber), cds, Cucumber genome sequencing project, version 2] with online tool (http://plantgrn.noble.org/psRNATarget/) . The circRNA-miRNA-mRNA regulating network was generated by Cytoscape (v3.6.0) . The eggNOG class annotations of these miRNAs targeted mRNAs were collected and used to interpret the possible regulating roles of circRNAs.
Validation of circRNAs
To further confirm the circRNAs identified from cucumber, total RNA was extracted from root and leaf tissues using Plant RNA Kit (Omega, London, UK), and then was subjected to reverse transcription reaction with the RevertAid First Strand cDNA Synthesis Kit (Thermo Scientific, USA) to synthesize the first cDNA strand according to the manufacturer’s protocol. The divergent PCR primers were designed using the “out-facing” strategy to guarantee the amplifications were from circle template. Primers were synthesized by Sangon Biotech Co., Ltd. (Shanghai, China). Sanger sequencing (BGI Co., Ltd., Wuhan, Hubei, China) were performed to validate the presence of the back-spliced junction sites. Quantitative real-time PCR (qRT-PCR) was performed using SYBR Green Master Mix (Vazyme, Nanjing, China) on a CFX 96 Real-Time PCR system (Bio-Rad) to validate the expression levels of circRNAs and correspondence parental genes according to reported methods . The relative quantity was calculated with 2–ΔΔCt method . Each sample contains three replications, and each replication contains two technical replications. Primers for these genes were designed using the Primer Premier 5 software and showed in Additional file 3: Table S3.
Proline measurements, enzymatic activity and qRT-PCR assay
Proline content was measured spectrophotometrically using the ninhydrin assay . The activities of proline dehydrogenase (PDH) and D1-pyrroline-5-carboxylate synthetase (P5CS) were determined according to Singh, et al. . Protein in the dialyzed extracts was measured colorimetrically at 595 nm based on the formation of protein-dye complex (the binding of Coomassie Brilliant blue G-250 to protein) as described by Bradford  using bovine serum albumin (BSA) as the standard. RNA extraction and quantitative qPCR analysis of the expression level of PDH (Csa5G643300) P5CS (Csa3G185140) were performed as described above with gene-specific primers (Additional file 3: Table S3).
Plant growth parameters under control and salt stress conditions
As can been seen in Fig. 1, salt stress significantly decreased plant dry weight, root superficial area, and root volume after 8 days of treatment. Root length was decreased by salt stress after 4 days of treatment compared with control (Fig. 1). Since the growth parameters have shown a downward trend in NaCl-treated seedlings compared with CK after 4 days of treatment, the 3 days’ treatment seedlings were used to identify circRNAs in cucumber because we were particularly interested in the early responses to salinity.
Identification of circRNAs in cucumber
To identify cucumber circRNAs, we extracted total RNA from leaf (L) and root (R) tissues under both control and salt stress conditions, constructed ribosomal RNA- and liner RNA-depleted RNA libraries, and deeply sequenced the libraries using Illumina Hiseq X-ten platform. In total, about 450 million raw reads were obtained from eight samples (CK-L1, CK-L2, Na-L1, Na-L2, CK-R1, CK-R2, Na-R1, and Na-R2). After removing low quality reads, clean reads were aligned to reference genome using BWA-MEM, about 334 million reads were uniquely mapped (Table 1). Then CIRI2, the circRNA identification tool, was used to inspect the produced SAM files, which resulted in the detection of 154,210 junction reads. Based on the back-spliced reads, 2787 circRNAs were finally detected, with 827 in leaf and 2420 in root (Fig. 2a, Additional file 4: Table S4). Among these 2787 circRNAs, 367 and 1960 circRNAs were specifically expressed in leaf and root, respectively, indicating that circRNAs are tissues specific (Fig. 2b). Besides, by divergent polymerase chain reaction (PCR), 35 out of 40 (87.5%) selected circRNAs were confirmed to be processed from back-splicing, demonstrating the reliability of our high-throughput sequencing based circRNA identification (Fig. 2d, e. Additional file 3: Table S3). The overall abundances of circRNAs in cucumber were positively correlated with their parental gene abundances both in leaf and root (r = 0.713, p < 0.005, Fig. 2c).
Characterization of cucumber circRNAs
Genomic origin analysis showed that the largest majority (2216, 79.51%) of the 2787 identified circRNAs were exonic circRNAs, 434 (15.57%) were intergenic circRNAs, and the rest 137 (4.92%) were intronic circRNAs (Fig. 3a). Genomic density analysis suggested that these cucumber circRNAs were distributed across the whole genome. Nevertheless, similar to coding genes, circRNAs were more commonly found distribution at both ends of chromosomes (Fig. 3b). As to chromosome distribution, the most circRNAs were found in chromosome 3, followed by chromosome 6 and 1 (Fig. 3c). Considering the length, most circRNAs (n = 2391; only known splice lengths without introns) are shorter than 1600 nucleotide (nt), the median length is ~ 570 nt, and large majority of long circRNAs are intergenic type (Fig. 3d. Additional file 4: Table S4). Alternative circularization analysis showed that 1131 alternative back-splicing circularization events were detected from 434 unique chromosome loci, of which, 319 from exon region, 66 from intergenic region, 6 from intron region, and 43 from mixture regions (including exon, intron and/or intergenic types). Statistically speaking, 295 of 434 loci produced two distinct circRNA isoforms, 75 produced three isoforms, 35 produced four isoforms, 12 produced five isoforms, and 2 regions produced nine different circRNAs isoforms (Fig. 3e, Additional file 4: Table S4). Then, we investigated the expression abundance of the circRNA isoforms within one gene locus and found that about 74.71% parental genes (845 in 1131, circRNAs without isoforms were excluded in the analysis) produced a significantly higher expressed circRNA (at least twofold higher than other isoforms) (Fig. 3f). These results indicate that there is often a predominantly expressed circRNA isoform from one gene locus.
To explore conserved features of circRNAs in different plant species, the circRNAs originated from the cognate orthologous genes of Arabidopsis thaliana, Oryza sativa, Glycine max, Zea mays, Solanum lycopersicum, Solanum tuberosum, Triticum aestivum and cucumber were compared. circRNAs belonging to these species were downloaded from the PlantcircBase . The orthologous gene sets between cucumber and these species were generated by BioMart from EnsemblPlants database. Among the 8362, 11,147, 3876, 2209, 984, 260, 31, and 1642 parent genes that produced exonic circRNAs in Arabidopsis, rice, soybean, maize, tomato, potato, barley, and cucumber, 1493 gene pairs between cucumber (1210 unigenes) and Arabidopsis (1462 unigenes), 1459 gene pairs between cucumber (1075 unigenes) and rice (1422 unigenes), 897 gene pairs between cucumber (656 unigenes) and soybean (877 unigenes), 526 gene pairs between cucumber (451 unigenes) and maize (502 unigenes), 241 gene pairs between cucumber (234 unigenes) and tomato (234 unigenes), 24 gene pairs between cucumber (24 unigenes) and potato (24 unigenes), and 6 gene pairs between cucumber (6 unigenes) and barley (5 unigenes) were orthologs (Additional file 5: Table S5), accounting for 17.48, 12.76, 22.63, 22.73, 22.78, 9.23, 16.13, and 95.83% of Arabidopsis, rice, soybean, maize, tomato, potato, barley, and cucumber genes producing exonic circRNAs, respectively. Furthermore, orthologous genes of 1175 cucumber parental genes (1175/1642, 72%) were found in more than 2 species, while 286 of that (286/1642, 17%) were found in 1 species and no orthologous genes of 181 cucumber parental genes were detected from other species (Additional file 6: Figure S1). Meanwhile, according to reciprocal BLAST search analysis, 257 cucumber circRNAs had significant alignments with 1540 circRNAs obtained from 12 species in PlantcircBase (e value <1e-5). The nucleotide identity levels between 257 cucumber circRNAs and significant alignments were > 70%. These results further suggested the conservation feature of cucumber circRNAs.
Orientation-opposite complementary sequences are not preferentially flanked by circularized exons
To further reveal the possible biosynthesis mechanisms of exon circRNAs, we then investigated the genomic features of more confidence exon circRNAs with high expression levels (circRNAs supported by at least five junction reads) (Fig. 4a). As can been seen in Fig. 4b, among the highly expressed exonic circular RNAs identified in this study, the vast majority were processed from those exons distributed in the middle of parent genes (2035), and only a few circRNAs contain the first (144) or the last exons (337) (Fig. 4b), implying that the formation of circular RNAs is commonly coupled with RNA splicing process. Moreover, most circular RNAs have multiple exons (generally 2–5 exons), while 101 out of 2035 highly expressed circular RNAs contain only one annotated exon (Fig. 4c). Interestingly, these exons from circular RNAs comprised by only one circularized exon were much longer than exons from circular RNAs containing multiple circularized exons, revealing that a certain length of exon(s) was probably preferred by exon circRNAs to maximize the circularization (Fig. 4d). Several studies proposed that, compared to animals, plant (e.g. Arabidopsis thaliana, rice and soybean) exonic circRNAs are less likely to be generated from exon(s) flanked by introns containing repetitive and/or reverse complementary sequences [10, 15, 43]. To explore whether or not cucumber is the same, we analyzed the flanking sequences of exon circRNAs and found that only 18 circRNAs contain at least 18-bp flanking complementary sequences at both end of 1000-bp flanking sequences (Additional file 7: Table S6). Meanwhile, miniature invertedrepeat transposable elements (MITEs) searching found that only 6 circRNAs have two flanking sequences that overlapped with MITEs. These results confirmed that the circularization of exon circRNAs in cucumber maybe not mediated by flanking orientation-opposite complementary sequences. In both plants and animals, U2-dependent spliceosome is involved in the splicing of vast majority of introns with 5’-GT and 3’-AG terminal dinucleotides . To confirm whether the 5’-GT and 3’-AG are the circRNA splicing signal in cucumber, the splice sites of 2787 cucumber circRNAs were further analyzed and we found that the splicing sites were dominated by C and T nucleotide in 5’ termini and A and G nucleotide in 3’ termini (Fig. 4e). And only a minority portion (325 or 11.7%) were flanked by the common GT/AG (including CT/AC equivalent) splicing signals (Fig. 4f). Moreover, a diverse set of non-GT/AG splicing signals were also found and the top five non-GT/AG terminal dinucleotides pairs were GA/AG (115), CT/AT (113), AT/AG (104), CT/TC (95), and CT/AG (75) (Fig. 4f). Accordingly, our analysis demonstrates that, except the canonical GT/AG, cucumber circRNAs also use diverse non-GT/AG splicing signals for circularization.
Differentially expressed circRNAs induced by salt stress
To overview the effect of salt stress on cucumber circRNAs, we pair-wised compared the whole expression levels of all detected circular RNAs. Wilcoxon rank-sum test indicated that circRNAs are significantly up-regulated by salt stress in both tissues (leaf, P < 6.07e-06; root, P < 2.2e-16. Figure 5a). In order to reveal cucumber circRNAs that have biological functions in responding to salt stress, we compared the expression profiles of differentially expressed circRNAs between control and salt stress treatments. Among the 2787 circRNAs, 1966 circRNAs were significant differential expressed under salt stress, including 1934 in root (991 up- and 943 down-regulated, Additional file 8: Table S7) and 44 in leaf (19 up- and 25 down-regulated, Additional file 9: Table S8) (Fig. 5b).
Functional categorization of stress-regulated parent genes
Functional classification by GO annotation
Since circRNA function could be related to the known function of the host linear transcripts [1, 3], we predicted and annotated the parent genes of regulated circRNAs. GO enrichment analysis showed that, in root, GO terms related to: (1) Biological Process (BP), including positive regulation of Ras GTPase activity (GO:0032320), xylan metabolic process (GO:0045491), regulation of multi-organism process (GO:0043900), vernalization response (GO:0010048), cell wall macromolecule biosynthetic process (GO:0044038), cell division (GO:0051301), and microtubule nucleation (GO:0007020); (2) Cellular Component (CC), including mitochondrial outer membrane (GO:0005741), and endoplasmic reticulum (GO:0005783); (3) Molecular Function (MF), including Ras GTPase activator activity (GO:0005099), zinc ion binding (GO:0008270), identical protein binding (GO:0042802), transition metal ion transmembrane transporter activity (GO:0046915), endonuclease activity, active with either ribo- or deoxyribonucleic acids and producing 5′-phosphomonoesters (GO:0016893), endoribonuclease activity (GO:0004521), protein serine/threonine kinase activity (GO:0004674), and oxidoreductase activity, acting on paired donors, with incorporation or reduction of molecular oxygen, 2-oxoglutarate as one donor, and incorporation of one atom each of oxygen into both donors (GO:0016706), were specifically enriched (Fig. 5c, Additional file 8: Table S7). Similar results were also reported in previous studies, which suggested that plant salt tolerance is a complex regulating network involving in signal transduction, and substance and energy metabolism [44, 45].
Genes related to abiotic stress responses and plant development
In leaf, the expression of seven annotated circRNAs, including Chr1:3244941|3,245,447, Chr1:20336460|20,339,620 Chr1:21205717|21,206,424, Chr1:25970551|25,973,020, Chr3:10198472|10,200,676, Chr4:23091126|23,091,996, and Chr6:9785683|9,788,312, were significantly down-regulated (Table 2). Simultaneously, the parent genes of these seven circRNAs have been proved to be vital for regulating abiotic stress responses and plant development. For example, protein phosphatase 2C (PP2C), members of major phosphatase class, have been reported to regulate abiotic stress triggered signaling pathways. In Arabidopsis, overexpression of OsPP108, a group A type PP2C from rice, increased the plant’s insensitive to ABA and tolerance to high salt and mannitol stresses during seed germination, root growth, and overall seedling growth . Plant growth is mainly the results of cell division and elongation, while salt stress inhibits cell division and cell wall metabolism [44, 46]. A decrease in cell wall extensibility has been suggested to be responsible for leaf growth reduction under salt stress [47, 48]. In this research, parent gene of circRNA Chr2:226834|229,565, annotated as type II inositol 1,4,5-trisphosphate 5-phosphatase FRA3-like, has been reported to plays an essential role in the secondary wall synthesis in fiber cells and xylem vessels . Another gene, annotated as molybdenum cofactor sulfurase (parent gene of circRNA Chr6:9785683|9,788,312), were found to be an important factor in regulating the ABA levels in plant tissues in response to drought, salt, or ABA treatment [50,51,52].
Similar to leaf, in the roots, several parent genes of differentially expressed circRNAs that belong to protein serine/threonine kinase activity (GO:0004674) including SNF1-related kinases (SnRKs) and PP2C, and cell process like 7 metabolic process (GO:0045491) and cell division (GO:0051301) were obtained (Table 3). As signal transductors, serine/threonine kinases involve in multiple stress tolerance in plants, including salt, drought, and freezing stresses [53, 54]. Moreover, a number of genes related to cell biogenesis including xyloglucan glycosyltransferase and sucrose synthase were identified (Table 3, Additional file 8: Table S7). Thus, a hypothesis was raised that these circRNAs were valuable parts in the signaling transduction pathways and cell wall biosynthesis and modification process through regulating the expression level of their parent genes.
Genes related to the ion transport
The ability of a plant to maintain a high cytosolic K+/Na+ ratio is critical to plant salt tolerance. Na+/H+ (NHX) antiporters and K+ transporter have important role in cellular Na+, K+ homeostasis [19, 55]. In this study, the expression level of four circRNAs, Chr3:25629954|25,631,297, Chr2:16186442|16,191,680, Chr3:7346435|7,347,924 and Chr5:2818033|2,823,737, were significantly affected by salinity (Table 3). Their parent genes were predicted to encode a vacuolar cation/H+ exchanger 2-like (CHX2-like), K+ efflux antiporter 6-like (KEA6-like), K+ efflux antiporter 2-like protein (KEA2-like), and sodium/hydrogen exchanger 7-like protein, (NHX7-like), respectively. NHX7, also known as SOS1, is a member of vacuolar NHX, and has been proved to improve salt tolerance through reducing Na+ accumulation in the shoot of Arabidopsis thaliana . The functions of NHX and KEA transporters are involved in K+ and H+ movements across endomembranes, and it is becoming more and more accepted that the expression level of these transporters are regulated by NaCl and osmotic stress [57, 58]. However, the detailed function of these antiporters in K+ homeostasis is poorly understood [57, 59]. These results in our study might have implications for the role of circRNAs related to Na+/H+ (NHX) antiporters and K+ transporters in cucumber and could be helpful to further discovered the mechanisms of the maintenance of optimal K+/Na+ ratios under saline conditions.
circRNA-miRNA-mRNA regulating network
Functional classification by eggNOG annotations
Besides cis-regulating parent genes, circRNAs showed an ability to act as competing endogenous RNAs . To explore the potential functional importance of circRNAs by binding to miRNA and further affect the post-transcriptional regulation of genes, we construct the circRNA-miRNA-mRNA network in cucumber. Among the 1966 differentially expressed circRNAs, 30 circRNAs were found to contain miRNA identical regions, implying that these circRNAs may have miRNA binding sites and could function as miRNA sponges, which need to be validated experimentally. Target prediction showed 45 miRNAs to be the potentially binding targets of these 30 circRNAs (Additional file 10: Table S9). Moreover, 269 mRNAs were predicted to be the potential targets of 26 miRNAs (Additional file 11: Table S11, Additional file 12: Table S12). The functions of 269 cucumber genes which transcript the potential target mRNAs were further analyzed by eggNOG class annotations analysis. Just like the GO terms analysis to circRNA parent genes, iron transporting and cell cycle controlling were obtained (Fig. 6).
Genes related to transcription factors
Besides, a lot of genes involved in transcription regulations and signal transduction were obtained, which has been reported to play important roles in response to environmental stresses . To address complicated environmental stresses, many transcription factor (TFs) genes are activated in plant cells through signal perception and subsequent signal transduction . Accordingly, in this study, many of the predicted miRNA targets are various TFs, including ethylene-responsive TFs, bZIP, MYB, SBP and ZF-HD (Additional file 10: Table S9). In plants, Ethylene is an important inducer of defence-related genes in both biotic and abiotic stresses [63, 64]. In this study, four RAP2-like ethylene-responsive transcription factors (Csa3G736760, Csa4G292470, Csa5G175970, and Csa6G296960) were predicted to be the targets of scaffold01066:88458_88477. It is reported that RAP2-like ethylene-responsive transcription factor participates in ABA, salt and osmotic stress responses in Arabidopsis . Moreover, ERF054-like ethylene-responsive transcription factor (Csa6G128610) and CRF4-like ethylene-responsive transcription factor (Csa5G139630) were predicted to be the targets of csa-miRn6-3p (Table 4. Additional file 10: Table S9). ERF054-like ethylene-responsive transcription factor is found to involved in salt stress response in radish (Raphanus sativus L.) and cold stress response in Rosa multiflora [66, 67]. Further studies are required to explore their exact roles in cucumber salt stress adaptation.
MYB transcription factors play impotent roles in diverse processes including developmental control, cell cycle control, secondary metabolism and responses to biotic and abiotic stresses in plant . In this study, several mRNA annotated as MYB-like were identified (Table 4. Additional file 10: Table S9). To be specific, targets of miR159b is functional annotated as MYB29-like transcription factor (Csa4G022940). In Arabidopsis, MYB29 has been revealed to function in the complex interplay of ethylene, jasmonic acid, and salicylic acid, and ROS signaling by regulating the expression of various ethylene responsive factors and WRKY transcription factors . In summary, combining with these former studies, our results reveal that those miRNAs and their corresponding circRNAs playing roles in salt stress responding by remolding the transcript of TFs (Fig. 6).
This circRNA-miRNA-mRNA network has been constructed in many studies related to plant circRNAs. The circRNA-miRNA interaction network could provide valuable candidates for experimentally texting the significant regulatory function by non-coding RNA regulatory network involving circRNAs in plants. However, this was associated with two problems. First, since only a few experimentally proven miRNA targets could be deciphered, it is hard to predict the functions of circRNAs more accurately. Second, unlike in animal, in plant, there is still no experimental evidence proves circRNAs can act as miRNA sponge and Chu, et al.  suggested that miRNA sponges might not be the major functional mechanism of plant circRNAs. Thus, we further mined conserved miRNAs, which have been well known to extensively regulate analogous targets and play critical roles in plant development and abiotic stresses . Some well-known miRNA families like miR156, miR159, miR172, miR398 miR399 and miR396 were predicted to be targeted by certain cucumber circRNAs. And most of these miRNAs were proved to be play roles in diverse interconnected biological programs and plays critical roles in responses to biotic and abiotic stresses. (Additional file 13: Table S10). In further studies, these conserved miRNAs can be important candidates to understand the circRNA-miRNA interactions as well as specific functions of these circRNAs in cucumber plant.
circRNAs mediate proline (pro) metabolisms by regulating corresponding biosynthesis and degradation genes
In this study, the cluster for ‘amino acid transport and metabolism’ was one of the largest groups among the annotated COG categories. In this category, circRNA Chr3:13041746|13,044,121 (corresponding parent gene Csa6G008780) and Chr6:915006|915,916 (corresponding parent gene Csa3G185140) are both significantly up-regulated under salt stress. Csa6G008780 encodes a rate-limiting enzyme involved in the biosynthesis of proline in higher plants, D1-pyrroline-5-carboxylate synthetase (P5CS). Csa3G185140 (homolog to Arabidopsis Drought and Freezing Responsive gene 1) has been reported to positively regulate proline accumulation through interacts with proline degradation enzymes PDH and P5CDH and compromises their activities in Arabidopsis . Further, we test whether circRNAs participate in proline synthesis and degradation regulation in physiological level. To this end, the expression level of Csa6G008780, Csa3G185140, and the activity of P5CS, PDH were determined. As a result, in the synthesis pathway, circRNA Chr6:915006|915,916 is significantly down-regulated whereas the expression level of P5CS (Csa6G008780) is significantly increased under salt stress. Accordingly, proline content and P5CS activity are significantly increased. Proline content is determined by proline synthesis and degradation. In the degradation pathway, the expression level of Chr6:915006|915,916 and Csa3G185140 is significantly up-regulated while PDH (Csa5G643300) expression level and PDH activity were significantly downregulated by salt stress, which is in consistent with the incensement of proline content under salt stress. These results suggested that Chr6:915006|915,916 may be response to salt stress and positively regulate Csa3G185140, which may also be responsible for the decrease of PDH (Csa5G643300) expression level and eventually PDH activity. Our results suggested that Chr3:13041746|13,044,121 and Chr6:915006|915,916 may act as a negative and positive regulator of their parent genes, respectively, which need to be further confirmed in the molecular level.
Circular RNAs (circRNAs) are a class of special RNAs generated by a non-linear backsplicing event between a downstream splice donor and an upstream splice acceptor. Recently, with the development of circRNA-seq, a huge number of circRNAs have been identified in both animals and plants. The discovering of these widely expressed and highly conserved circRNAs has increased the potential impact of non-coding RNAs on cellular function . Compared with animals, the biogenesis, regulation, and function of circRNAs are largely unknown in plants . Ye, et al.  characterized circRNAs in Oryza sativa and Arabidopsis thaliana and reported that the biogenesis mechanisms of circRNAs from plants could be different from those in animals. Further studies showed that plant circRNAs are conserved, expressed at low concentrations, and in a tissue-specific manner. Moreover, recent studies in rice, wheat, kiwifruit, and potato demonstrated that the circRNAs play an important role in regulating the response to environmental stress such as phosphate imbalance, dehydration stress, and pathogen invasion [10, 12,13,14]. Cucumber is an important horticultural crop and is fairly sensitive to salinity. To our knowledge, there is still no evidence showing the role of circRNAs in salt-stressed cucumber.
In the present study, we reported the first genome-wide identification and characterization of circRNAs in cucumber and the potential regulation roles of circRNAs under salt stress condition. circRNAs can derive from both root and leaf. In total, 2787 circRNAs were detected, including 827 from leaf and 2420 from root (Table 1. Additional file 4: Table S4). In consistent with previous studies, we found that circRNAs are tissues specific, with 367 and 1960 specifically expressed in leaf, and root, respectively [15, 74]. The circRNAs identified in cucumber was less than that in Oryza sativa (26,160), Arabidopsis thaliana (31,079), soybean (5372), and kiwifruit (3582), but more than that in wheat (88), maize (496), barely (47), tea (342), and tomato (854), which may due to the difference in species, methods (e.g. experimental strategy or bioinformatic method), tissues (e.g. leaf, root, fruit, stem, shoot), and circRNA prediction tools (e.g. CIRI2, fnd_circ, CIRCexplorer), that are used in different researches . Besides, by divergent PCR, 35 of these circRNAs were validated to be processed from back-splicing (Fig. 2d, e and Additional file 2: Table S2), confirming the reliability of our circRNA-seq results (Additional file 4: Table S4). This large amount of circRNAs detected in this study indicated that circRNAs may contain one of the largest RNA families in cucumber transcription. It also increased the potential impact of non-coding RNAs on cell function and the complexity of regulatory processes . Examining the roles of individual circRNAs will be the subject of future investigations.
CircRNAs can be produced from exons (exonic circRNA), introns (intronic circRNA), and intergenic regions of the genome . In rice, Arabidopsis, tea and tomato, most identified circRNAs were exonic circular RNAs [10, 16, 43]. Similarly, in this study, the majority of the identified circRNAs (2216, 79.51%) were exonic circRNAs (Fig. 3a). However, in wheat and kiwifruit, 51 and 60.2% of the circRNAs were intergenic circRNAs, respectively. In soybean, most of the total circRNAs (in leaf, root, and stem) were intronic circRNAs, while there are more intronic in the root and more exonic in the stem . Low percentage of exonic circRNAs in wheat might due to the huge genome size with comparable fewer available genes and higher percentage of introns circRNAs in soybean might be attributed to the large and duplicated genome and multiple copies of genes [14, 15]. Integrating, these results indicated that the molecular basis of circular RNA biogenesis in plants is quite complicated, which is species- and tissue-specific [13,14,15]. Since biogenesis of circRNAs have close relationship with their structure, we then explored the characteristics of cucumber circular RNAs.
The different isoforms of circRNAs are results of alternative back-spliced circularization. In human, studies have shown that exon circularization depends on flanking intronic complementary sequences. Further evidence suggested that alternative formation of inverted repeated Alu pairs and the competition between them can lead to alternative circularization. Our results revealed that most alternative circularization events were located in the exon regions, but, interestingly, most alternative circularization with more than 3 isoforms were located in the intergenic regions (Fig. 3e). Besides, chromosome 3, the longest one in cucumber genome, produced more circRNAs than others. Further analysis revealed that circRNAs numbers are significantly correlated with chromosome length (Additional file 14: Figure S2, R2 = 0.9298 ***). Interestingly, these exons detected in circRNAs with only one circularized exon are much longer than these exons in circRNAs with multiple circularized exons (Fig. 4d), revealing that circRNAs processing may prefer a certain length to maximize exon(s) circularization. Similarly, most circular RNAs in H9 human embryonic stem cells contain two to three exons and circular RNAs with only one circularized exon were much longer than those from circular RNAs with multiple circularized exons . But whether it is a common character need to be further confirmed in both animal and plant.
Compared to animals, plant exonic circRNAs are less likely to be generated from exon(s) flanked by introns containing repetitive and/or reverse complementary sequences. In rice, Lu, et al.  showed that only a very small part of rice circular RNAs contain two flanking sequences overlapping with miniature invertedrepeat transposable elements (MITEs) through searching 500-bp flanking sequences of the backsplice sites. Furthermore, Ye, et al.  found that, unlike human exonic circRNAs, most of which are flanked by canonical GT/AG splicing signals, the vast majority of rice circRNAs contain very diverse non-GT/AG splicing signals. To determine whether exon circularization depends on flanking intronic complementary sequences in cucumber, we detected putative complementary sequences by testing 1000-bp flanking sequences of the backsplice sites from circRNAs and only 18 circRNAs were detected to contain flanking complementary sequence longer than 17-bp. Furthermore, we compared the cucumber circRNAs flanking sequences with plant miniature invertedrepeat transposable elements (MITEs) and only 6 circRNAs were found to have two flanking sequences that overlapped with MITEs. To further validate this observation, we expanded the flanking sequences to 2-kb, and, as expected, still a small number of circRNAs (1.5%) have flanking complementary sequences longer than 17-bp, suggesting the biosynthesis of cucumber circRNAs possibly not dependent on flanking complementary sequences (Additional file 3: Table S3). Similar result was also observed in rice. Therefore, in cucumber, other mechanisms may be involved in circRNA generation in addition to pairing-driven circularization.
The responses of plant to salt stress are very complex processes and involve many salt-inducible genes and signaling transduction pathways. In this study, we explored the roles of circRNAs in cucumber in responses to salt stress and found 1966 circRNAs significantly differentially expressed under salt stress compared with CK (Fig. 5. Additional file 8: Table S7, Additional file 9: Table S8). Among these 1966 circRNAs, a large number of them (1934) belong to the root and only 44 belong to the leaf, which may be attributable to: (1) as suggested by previous studies, levels of circular RNAs vary across tissues, species and treatments. For example, circRNAs exhibited specific expression patterns in different tissues of soybean, with 776, 3171 and 2165 being identified from leaves, roots and stems, respectively (Zhao et al. 2017). circRNAs are also treatment-specific. In Arabidopsis, heat stress induced more circRNAs than in normal condition . In barely, circRNAs behaved differently in leaf and seed tissues when comparing the treatments with foliar application of micronutrients . (2) Root is the first organ that respond to salt stress, which senses salt signals to initial reestablishment of cellular ionic homeostasis and osmotic equilibrium under stress conditions . Some stress-responsive genes function primarily in the root. The biological functions of circRNA could be related to the known function of the host linear transcripts [10, 43]. GO enrichment analysis on parent genes of these differentially expressed circRNAs showed that GO terms related to salt stress response including amino acid transport and metabolism, energy production and conversion, and ion transporting were specifically enriched (Fig. 5c. Additional file 8: Table S7) . Thus, a hypothesis was raised that these circRNAs were valuable parts in cucumber in response to salt stress through regulating the expression level of their parent genes.
Amino acids have various prominent functions in plants including stress response . Proline (Pro) biosynthesis and degradation are known to undergo alteration in response to diverse environmental stresses. Moreover, proline has been considered to be one of the most commonly occurring osmolytes and scavenger of reactive oxygen species (ROS) under salt stress . In plants, proline is synthesized from glutamate by the action of two enzymes, P5CS and ∆1-pyrroline-5-carboxylate reductase (P5CR). P5CS converts glutamate into pyrroline-5-carboxylate (P5C) and has been proposed as a rate-limiting enzyme involved in the biosynthesis of proline in higher plants. Conversely, proline is oxidized to glutamate via PDH and ∆1-pyrroline-5-carboxylate dehydrogenase (P5CDH) . In this study, circRNA Chr3:13041746|13,044,121 and Chr6:915006|915,916 may act as a negative and positive regulator of their parent genes Csa5G643300 and Csa3G185140 respectively, which were reported to participate in proline biosynthesis and degradation (Fig. 7). The production of circRNAs has been reported to have a negative or positive effect on expression of their parental gene. For example, in rice, Lu, et al.  suggested that circRNA and its linear form might act as a negative regulator of its parental gene. Li, et al.  proposed that intron-containing circRNAs (EIciRNAs) interact with U1 snRNP and promote transcription of their parental genes. To further explore the relationship between circRNAs and their parent genes, the expression level of 6 random selected circRNAs and their parent mRNAs in different time points were examined. Results suggested that in most cases, the expression patterns of parent mRNAs showed an opposite trend with circRNAs. But there are exceptions, the expression levels of CsCir12 were coincident with the increased expression of their parental genes after 3, 6 and 9 days of treatment (Additional file 15: Figure S3). In barley, weak and negative correlations between cellular levels of circular RNAs and the linear-transcript levels of their parental genes were reported . However, expression profiles of circRNAs in some plant including tea and Arabidopsis showed a positive correlation with their parental genes (Tong et al. 2018) [16, 19]. Comprehensive identification and quantification of circRNAs in more plant species are needed to reveal the relationship between the abundances of circRNAs and their corresponding parental genes. The complicate relationship between circRNAs and their parent genes need to be further explored by molecular biology approaches.
We first demonstrated the existence of abundant circRNAs in cucumber and their possible regulating roles in response to salt stress. Characterization analysis revel two key findings in our work: (1) circRNAs prefer to arise from middle exon(s) of parent genes and specific length of exons are favorited by circularization; (2) Unlike animals, non-GT/AG splicing signals are common in circRNAs and pairing-driven circularization is not the major way to generate circRNAs. Annotation and enrichment analysis on differentially expressed circRNAs reveal that these circRNAs are associated with numerous cucumber coding genes generally functional related to salt stress response, including transcription, signal transcription, cell wall biogenesis, metabolism adaptation (e.g. carbohydrate transport and metabolism, and amino acid transport and metabolism) and ion homeostasis, suggesting circRNAs’ wide involvement of participating in cucumber responding to salt stress (Fig. 7). Our study expands our understanding about the characters of plant circRNAs and the regulating diversity of cucumber transcriptome, which will facilitate the future regulatory mechanism analyses of cucumber circRNA in responding to salt stress.
Cation/H+ exchanger 2-like
False discovery rate
K+ efflux antiporter 2-like protein
K+ efflux antiporter 6-like
Kyoto Encyclopedia of Genes and Genomes
Clusters of Orthologous Groups of proteins
Miniature invertedrepeat transposable elements
National Center for Biotechnology Information
Sodium/hydrogen exchanger 7-like protein
NCBI non-redundant protein sequences
Protein phosphatase 2C
Quality real time PCR
Number of circular reads/number of mapped reads (units in billion)/read length
A manually annotated and reviewed protein sequence database
Yin JL, Liu MY, Ma DF, Wu JW, Li SL, Zhu YX, Han B. Identification of circular RNAs and their targets during tomato fruit ripening. Postharvest Biol Tec. 2018;136:90–8.
Salzman J. Circular RNA expression: its potential regulation and function. Trends Genet. 2016;32(5):309–16.
Li L, Guo J, Chen Y, Chang C, Xu C. Comprehensive CircRNA expression profile and selection of key CircRNAs during priming phase of rat liver regeneration. BMC Genomics. 2017;18(1):80.
Zheng Q, Bao C, Guo W, Li S, Chen J, Chen B, Luo Y, Lyu D, Li Y, Shi G, et al. Circular RNA profiling reveals an abundant circHIPK3 that regulates cell growth by sponging multiple miRNAs. Nat Commun. 2016;7:11215.
Wang K, Bo L, Fang L, Wang JX, Liu CY, Bing Z, Zhou LY, Teng S, Man W, Tao Y. A circular RNA protects the heart from pathological hypertrophy and heart failure by targeting miR-223. Eur Heart J. 2016;37(33):2602–11.
Li Z, Huang C, Bao C, Chen L, Lin M, Wang X, Zhong G, Yu B, Hu W, Dai L, et al. Exon-intron circular RNAs regulate transcription in the nucleus. Nat Struct Mol Biol. 2015;22:256.
Du WW, Yang W, Liu E, Yang Z, Dhaliwal P, Yang BB. Foxo3 circular RNA retards cell cycle progression via forming ternary complexes with p21 and CDK2. Nucleic Acids Res. 2016;44(6):2846–58.
Memczak S, Jens M, Elefsinioti A, Torti F, Krueger J, Rybak A, Maier L, Mackowiak SD, Gregersen LH, Munschauer M, et al. Circular RNAs are a large class of animal RNAs with regulatory potency. Nature. 2013;495(7441):333–8.
Wang PL, Bao Y, Yee MC, Barrett SP, Hogan GJ, Olsen MN, Dinneny JR, Brown PO, Salzman J. Circular RNA is expressed across the eukaryotic tree of life. PLoS One. 2014;9(3):e90859.
Ye CY, Chen L, Liu C, Zhu QH, Fan LJ. Widespread noncoding circular RNAs in plants. New Phytol. 2015;208(1):88–95.
Tan J, Zhou Z, Niu Y, Sun X, Deng Z. Identification and functional characterization of tomato circrnas derived from genes involved in fruit pigment accumulation. Sci Rep. 2017;7:8594.
Zhou R, Zhu YX, Zhao J, Fang ZW, Wang SP, Yin JL, Chu ZH, Ma DF. Transcriptome-wide identification and characterization of potato circular RNAs in response to Pectobacterium carotovorum subspecies brasiliense infection. Int J Mol Sci. 2018;19(1):71.
Wang Z, Liu Y, Li D, Li L, Zhang Q, Wang S, Huang H. Identification of circular RNAs in kiwifruit and their species-specific response to bacterial canker pathogen invasion. Front Plant Sci. 2017;8:413.
Wang Y, Yang M, Wei S, Qin F, Zhao H, Suo B. Identification of circular RNAs and their targets in leaves of Triticum aestivum L. under dehydration stress. Front Plant Sci. 2016;7:2024.
Zhao W, Cheng Y, Zhang C, You Q, Shen X, Guo W, Jiao Y. Genome-wide identification and characterization of circular RNAs by high throughput sequencing in soybean. Sci Rep. 2017;7(1):5636.
Wei T, Jie Y, Yan H, Li F, Zhou Q, Wei C, Bennetzen JL. Circular RNA architecture and differentiation during leaf bud to young leaf development in tea (Camellia sinensis). Planta. 2018;248(10):1–13.
Darbani B, Noeparvar S, Borg S. Identification of circular RNAs from the parental genes involved in multiple aspects of cellular metabolism in barley. Front Plant Sci. 2016;7:776.
Chen L, Zhang P, Fan Y, Lu Q, Li Q, Yan J, Muehlbauer GJ, Schnable PS, Dai M, Li L. Circular RNAs mediated by transposons are associated with transcriptomic and phenotypic variation in maize. New Phytol. 2018;217(3):3.
Cheng J, Zhang Y, Li Z, Wang T, Zhang X, Zheng B. A lariat-derived circular RNA is required for plant development in Arabidopsis. Sci China Life Sci. 2018;61(2):204–13.
Pan T, Sun X, Liu Y, Li H, Deng G, Lin H, Wang S. Heat stress alters genome-wide profiles of circular RNAs in Arabidopsis. Plant Mol Biol. 2018;96(3):217–29.
Ye CY, Zhang X, Chu Q, Liu C, Yu Y, Jiang W, Zhu Q-H, Fan LJ, Guo L. Full-length sequence assembly reveals circular RNAs with diverse non-GT/AG splicing signals in rice. RNA Biol. 2017;14(8):1055–63.
Conn VM, Hugouvieux V, Nayak A, Conos SA, Capovilla G, Cildir G, Jourdain A, Tergaonkar V, Schmid M, Zubieta C, et al. A circRNA from SEPALLATA3 regulates splicing of its cognate mRNA through R-loop formation. NAT PLANTS. 2017;3(5):17053.
Yang Y, Guo Y. Elucidating the molecular mechanisms mediating plant salt-stress responses. New Phytol. 2018;217(2):523–39.
Yin JL, Jia J, Lian Z, Hu Y, Guo J, Huo H, Zhu YX, Gong H. Silicon enhances the salt tolerance of cucumber through increasing polyamine accumulation and decreasing oxidative damage. Ecotox Environ Safe. 2019;169:8–17.
Zhu YX, Xu XB, Hu YH, Han WH, Yin JL, Li HL, Gong HJ. Silicon improves salt tolerance by increasing root water uptake in Cucumis sativus L. Plant Cell Rep. 2015;34(9):1629–46.
Cao X, Wu Z, Jiang F, Zhou R, Yang Z. Identification of chilling stress-responsive tomato microRNAs and their target genes by high-throughput sequencing and degradome analysis. BMC Genomics. 2014;15(1):1130.
Xia JH, Liu P, Liu F, Lin G, Sun F, Tu R, Yue GH. Analysis of stress-responsive transcriptome in the intestine of Asian seabass (Lates calcarifer) using RNA-Seq. DNA Res. 2013;20(5):449–60.
Li Y, Li C, Bai L, He C, Yu X. MicroRNA and target gene responses to salt stress in grafted cucumber seedlings. Acta Physiol Plant. 2016;38(2):1–12.
Gao Y, Zhang J, Zhao F. Circular RNA identification based on multiple seed matching. Brief Bioinform. 2017;1:8.
Burkhardt A, Day B. Transcriptome and small RNAome dynamics during a resistant and susceptible interaction between cucumber and downy mildew. Plant Genome. 2016;9(1):1.
Jin W, Wu F. Identification and characterization of cucumber microRNAs in response to Pseudoperonospora cubensis infection. Gene. 2015;569(2):225–32.
Li C, Li Y, Bai L, Zhang T, He C, Yan Y, Yu X. Grafting-responsive miRNAs in cucumber and pumpkin seedlings identified by high-throughput sequencing at whole genome level. Physiol Plantarum. 2014;151(4):406–22.
Liu HW, Luo LX, Liang CQ, Jiang N, Liu PF, Li JQ. High-throughput sequencing identifies novel and conserved cucumber (Cucumis sativus L.) microRNAs in response to cucumber green mottle mosaic virus infection. PLoS One. 2015;10(6):e0129002.
Mao W, Li Z, Xia X, Li Y, Yu J. A combined approach of high-throughput sequencing and degradome analysis reveals tissue specific expression of microRNAs and their targets in cucumber. PLoS One. 2012;7(3):e33040.
Martinez G, Forment J, Llave C, Pallas V, Gomez G. High-throughput sequencing, characterization and detection of new and conserved cucumber miRNAs. PLoS One. 2011;6(5):e19523.
Bo X, Wang S. TargetFinder: a software for antisense oligonucleotide target site selection based on MAST and secondary structures of target mRNA. Bioinformatics. 2005;21(8):1401–2.
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13(11):2498–504.
Yin JL, Fang ZW, Sun C, Zhang P, Zhang X, Lu C, Wang SP, Ma DF, Zhu YX. Rapid identification of a stripe rust resistant gene in a space-induced wheat mutant using specific locus amplified fragment (SLAF) sequencing. Sci Rep. 2018;8(1):3086.
Bates LS, Waldren RP, Teare ID. Rapid determination of free proline for water-stress studies. Plant Soil. 1973;39(1):205–7.
Singh A, Jha SK, Bagri J, Pandey GK. ABA inducible rice protein phosphatase 2C confers ABA insensitivity and abiotic stress tolerance in Arabidopsis. PLoS One. 2015;10(4):e0125168.
Bradford MM. A rapid and sensitive method for the quantitation of microgram quantities of protein utilizing the principle of protein-dye binding. Anal Biochem. 1976;72(s1–2):248–54.
Chu Q, Zhang X, Zhu X, Liu C, Mao L, Ye C, Zhu Q, Fan L. PlantcircBase: a database for plant circular RNAs. Mol Plant. 2017;10(8):1126–8.
Lu T, Cui L, Zhou Y, Zhu C, Fan D, Gong H, Zhao Q, Zhou C, Zhao Y, Lu D. Transcriptome-wide investigation of circular RNAs in rice. RNA. 2015;21(12):2076–87.
Kim SW, Jeong HJ, Jung KH. Integrating omics analysis of salt stress-responsive genes in rice. Genes Genom. 2015;37(8):645–55.
Sun W, Xu X, Zhu H, Liu A, Liu L, Li J, Hua X. Comparative transcriptomic profiling of a salt-tolerant wild tomato species and a salt-sensitive tomato cultivar. Plant Cell Physiol. 2010;51(6):997–1006.
Burssens S, Himanen K, Cotte BVD, Beeckman T, Montagu MV, Inzé D, Verbruggen N. Expression of cell cycle regulatory genes and morphological alterations in response to salt stress in Arabidopsis thaliana. Planta. 2000;211(5):632–40.
Liu L, Si L, Meng X, Luo L. Comparative transcriptomic analysis reveals novel genes and regulatory mechanisms of Tetragenococcus halophilus in response to salt stress. J Ind Microbiol Biot. 2015;42(4):601–16.
Uddin MN, Hanstein S, Leubner R, Schubert S. Leaf cell-wall components as influenced in the first phase of salt stress in three maize (Zea mays L.) hybrids differing in salt resistance. J Agron Crop Sci. 2013;199(6):405–15.
Zhong R, Burk DH, Morrison WH, Ye Z-H. FRAGILE FIBER3, an Arabidopsis gene encoding a type II inositol polyphosphate 5-phosphatase, is required for secondary wall synthesis and actin organization in FIBER cells. Plant Cell. 2004;16(12):3242–59.
Huang PM, Chen JY, Wang SJ. Tissue-specific regulation of rice molybdenum cofactor sulfurase gene in response to salt stress and ABA. Acta Physiol Plant. 2009;31(3):545–51.
Lu Y, Li Y, Zhang J, Xiao Y, Yue Y, Duan L, Zhang M, Li Z. Overexpression of Arabidopsis molybdenum cofactor sulfurase gene confers drought tolerance in maize (Zea mays L.). PLoS One. 2013;8(1):e52126.
Xiong L, Ishitani M, Lee H, Zhu JK. The Arabidopsis LOS5/ABA3 locus encodes a molybdenum cofactor sulfurase and modulates cold stress- and osmotic stress-responsive gene expression. Plant Cell. 2001;13(9):2063–83.
Mao X, Zhang H, Tian S, Chang X, Jing R. TaSnRK2.4, an SNF1-type serine/threonine protein kinase of wheat (Triticum aestivum L.), confers enhanced multistress tolerance in Arabidopsis. J Exp Bot. 2010;61(3):683–96.
Sun XL, Yu QY, Tang LL, Ji W, Bai X, Cai H, Liu XF, Ding XD, Zhu YM. GsSRK, a G-type lectin S-receptor-like serine/threonine protein kinase, is a positive regulator of plant tolerance to salt stress. J Plant Physiol. 2013;170(5):505–15.
Zhu JK. Salt and drought stress signal transduction in plants. Annu Rev Plant Biol. 2002;53(53):247–73.
Shi HZ, Lee BH, Wu SJ, Zhu JK. Overexpression of a plasma membrane Na+/H+ antiporter gene improves salt tolerance in Arabidopsis thaliana. Nat Biotechnol. 2002;21(1):81–5.
Cuin TA, Shabala S. Compatible solutes mitigate damaging effects of salt stress by reducing the impact of stress-induced reactive oxygen species. Plant Signal Behav. 2008;3(3):207–8.
Mottaleb SA, Rodríguez-Navarro A, Haro R. Knockouts of physcomitrella patens CHX1 and CHX2 transporters reveal high complexity of potassium homeostasis. Plant Cell Physiol. 2013;54(9):1455–68.
Han L, Li JL, Wang L, Shi WM, Su YH. Identification and localized expression of putative K+/H+ antiporter genes in Arabidopsis. Acta Physiol Plant. 2015;37(5):101.
Hansen TB, Jensen TI, Clausen BH, Bramsen JB, Finsen B, Damgaard CK, Kjems J. Natural RNA circles function as efficient microRNA sponges. Nature. 2013;495(7441):384–8.
Du X, Wang G, Ji J, Shi L, Guan C, Jin C. Comparative transcriptome analysis of transcription factors in different maize varieties under salt stress conditions. Plant Growth Regu. 2017;81(1):183–95.
Gao F, Zhou J, Deng RY, Zhao HX, Li CL, Chen H, Suzuki T, Park SU, Wu Q. Overexpression of a tartary buckwheat R2R3-MYB transcription factor gene, FtMYB9, enhances tolerance to drought and salt stresses in transgenic Arabidopsis. J Plant Physiol. 2017;214:81.
Pan Y, Seymour GB, Lu C, Hu Z, Chen X, Chen G. An ethylene response factor (ERF5) promoting adaptation to drought and salt tolerance in tomato. Plant Cell Rep. 2012;31(2):349–60.
Wang KLC, Li H, Ecker JR. Ethylene biosynthesis and signaling networks. Plant Cell. 2002;14(suppl 1):S131–51.
Zhu Q, Zhang J, Gao X, Tong J, Xiao L, Li W, Zhang H. The Arabidopsis AP2/ERF transcription factor RAP2.6 participates in ABA, salt and osmotic stress responses. Gene. 2010;457(1):1–12.
Sun X, Xu L, Wang Y, Luo X, Zhu X, Kinuthia KB, Nie S, Feng H, Li C, Liu L. Transcriptome-based gene expression profiling identifies differentially expressed genes critical for salt stress response in radish (Raphanus sativus L.). Plant Cell Rep. 2016;35(2):329–46.
Zhang X, Zhang J, Zhang W, Yang T, Xiong Y, Che D. Transcriptome sequencing and de novo analysis of Rosa multiflora under cold stress. Acta Physiol Plant. 2016;38(7):164.
Dubos C, Stracke R, Grotewold E, Weisshaar B, Martin C, Lepiniec L. MYB transcription factors in Arabidopsis. Trends Plant Sci. 2010;15(10):573–81.
Zhang X, Ivanova A, Vandepoele K, Radomiljac JD, Van de Velde J, Berkowitz O, Willems P, Xu Y, Ng S, Van Aken O, et al. The transcription factor MYB29 is a regulator of ALTERNATIVE OXIDASE 1. Plant Physiol. 2017;173(3):1824–43.
Chu Q, Bai P, Zhu X, Zhang X, Mao L, Zhu Q, Fan L, Ye C. Characteristics of plant circular RNAs. Brief Bioinform. 2018. https://doi.org/10.1093/bib/bby111.
Qin Z, Li C, Mao L, Wu L. Novel insights from non-conserved microRNAs in plants. Front Plant Sci. 2014;5:586.
Ren Y, Miao M, Meng Y, Cao J, Fan T, Yue J, Xiao F, Liu Y, Cao S. DFR1-mediated inhibition of proline degradation pathway regulates drought and freezing tolerance in Arabidopsis. Cell Rep. 2018;23(13):3960–74.
Errichelli L, Dini Modigliani S, Laneve P, Colantoni A, Legnini I, Capauto D, Rosa A, De Santis R, Scarfo R, Peruzzi G, et al. FUS affects circular RNA expression in murine embryonic stem cell-derived motor neurons. Nat Commun. 2017;8:14741.
Shen Y, Guo X, Wang W. Identification and characterization of circular RNAs in zebrafish. FEBS Lett. 2016;591(1):213.
Yin JL, Ma DF, Liu LC, Xia YC, Zhu YX. Biology features of circular RNAs and their research progress in plants. Acta Bot Boreali-Oc. 2017;37(12):2510–8.
Zhang XO, Wang HB, Zhang Y, Lu X, Chen LL, Yang L. Complementary sequence-mediated exon circularization. Cell. 2014;159(1):134–47.
Hildebrandt TM, Nunes Nesi A, Araújo WL, Braun HP. Amino acid catabolism in plants. Mol Plant. 2015;8(11):1563–79.
We thank Prof. Yongli Qiao for beneficial comments on the initial project design and data analysis.
This work was supported by the National Natural Science Foundation of China (No. 31701911), the Open Project Program of Engineering Research Center of Ecology and Agricultural Use of Wetland, Ministry of Education (No. KF201707). The funders had no role in the experiment design, data analysis, decision to publish, or preparation of the manuscript.
Availability of data and materials
We have deposited our data in the National Center for Biotechnology Information (NCBI) Sequence Read Archive with a Bioproject ID: PRJNA449130.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Physical position of unanchored scaffolds in merged ChrUn. (XLSX 17 kb)
Table S2. The mature sequences of reported cucumber miRNAs. (XLSX 87 kb)
Table S3. Primers used for circRNAs sequences validation. (XLSX 13 kb)
Table S4. Detailed information of identified circRNAs. (XLSX 726 kb)
Table S5. circRNAs derived from orthologous genes of cucumber and seven other species. (XLSX 235 kb)
Figure S1. The number of gene pairs in the parent genes of cucumber circRNAs and other species for orthologs. (DOCX 30 kb)
Table S6. circRNAs with complementary or MITE flanking sequences. (XLSX 18 kb)
Table S7. Significantly differentially expressed circRNAs in root. (XLSX 474 kb)
Table S8. Significantly differentially expressed circRNAs in leaf. (XLSX 19 kb)
Table S9. miRNAs sponged by circRNAs. (XLSX 17 kb)
Table S11. miRNAs sponged by significantly differentially expressed circRNAs and their mRNA targets in root. (XLSX 90 kb)
Table S12. miRNAs sponged by significantly differentially expressed circRNAs and their mRNA targets in leaf. (XLSX 39 kb)
Table S10. Evolutionary conserved miRNAs sponged by circRNAs. (XLSX 13 kb)
Figure S2. Chromosome length versus circRNA number. (DOCX 81 kb)
Figure S3. Expression analysis of circRNAs and their parent genes. (DOCX 308 kb)
About this article
Cite this article
Zhu, YX., Jia, JH., Yang, L. et al. Identification of cucumber circular RNAs responsive to salt stress. BMC Plant Biol 19, 164 (2019). https://doi.org/10.1186/s12870-019-1712-3
- Cucumis sativus
- Cucumber circRNAs
- Salt stress
- Parent genes
- GO enrichment