Skip to main content
Download PDF
  • Research article
  • Open access
  • Published:

Identification and characterization of CircRNAs involved in the regulation of wheat root length

Biological Research volume 52, Article number: 19 (2019) Cite this article

Abstract

Background

Recent studies indicate that circular RNAs (circRNAs) may play important roles in the regulation of plant growth and development. Plant roots are the main organs of nutrient and water uptake. However, whether circRNAs involved in the regulation of plant root growth remains to be elucidated.

Methods

LH9, XN979 and YN29 are three Chinese wheat varieties with contrasting root lengths. Here, the root circRNA expression profiles of LH9, XN979 and YN29 were examined by using high-throughput sequencing technology.

Results

Thirty-three and twenty-two differentially expressed circRNAs (DECs) were identified in the YN29-LH9 comparison and YN29-XN979 comparison, respectively. Among them, ten DECs coexisted in both comparisons. As the roots of both LH9 and XN979 were significantly larger and deeper than YN29, the ten DECs coexisting in the two comparisons were highly likely to be involved in the regulation of wheat root length. Moreover, three of the ten DECs have potential miRNA binding sites. Real-time PCR analysis showed that the expression levels of the potential binding miRNAs exhibited significant differences between the long root plants and the short root plants.

Conclusions

The expression levels of some circRNAs exhibited significant differences in wheat varieties with contrasting root phenotypes. Ten DECs involved in the regulation of wheat root length were successfully identified in which three of them have potential miRNAs binding sites. The expression levels of putative circRNA-binding miRNAs were correlated with their corresponding circRNAs. Our results provide new clues for studying the potential roles of circRNAs in the regulation of wheat root length.

Background

Plant roots are the main organs of nutrient and water uptake from soil. Genetic improvement of root traits is of vital importance for improving crop nutrient and water use efficiencies [1,2,3,4,5]. Understanding the molecular mechanisms controlling plant root growth and development is beneficial to molecular breeding aimed at improving root traits and resource utilization efficiency of crops. Previous studies have shown that plant root growth and development are controlled by phytohormones and sustained by the root apical meristem (RAM) [6,7,8,9,10,11,12]. Besides, redox regulation and the balance of ROS also have important roles in maintaining RAM activity [13,14,15]. Therefore, the growth and development of plant roots are regulated by complex hormone signals and other related pathways.

CircRNAs are endogenous non-coding RNAs produced by back-splicing of pre-mRNA [16]. The 5′ and 3′ ends of the circRNA are linked together to form a covalent closed loop structure [17]. Currently, circRNAs have been widely studied in animals [18,19,20,21]. However, the role of circRNA in plants has not attracted enough attention [22]. Until recent years, several literatures reported circRNAs identified in plants such as Arabidopsis [23,24,25,26,27], soybean [28, 29], rice [25, 30], maize [31,32,33], tomato [34, 35], barley [36], tea [37], cotton [38, 39] and wheat [40, 41]. For example, Ye et al. performed genome-wide identification of circRNAs in rice and Arabidopsis using available public RNA-Seq data and identified 12,037 and 6012 circRNAs, respectively. Moreover, the parent genes of over 700 exonic circRNAs were orthologues between rice and Arabidopsis, suggesting the conservation of circRNAs in plants [25]. Chen et al. performed circRNA-Seq on maize seedling leaves and uncovered 2804 circRNAs. They found that sequences related to LINE1-like elements (LLEs) and their reverse complementary pairs (LLERCPs) are significantly enriched in the flanking regions of circRNAs [31]. Furthermore, genes with LLERCP-mediated circRNAs are enriched among loci that are associated with phenotypic variation. Therefore, circRNAs are likely to be involved in the modulation of phenotypic variation by LLERCPs [31]. Moreover, studies showed that circRNAs may also play roles in response to Verticillium wilt in cotton [39], maize iranian mosaic virus infection in maize [32], TYLCV infection in tomato [35], drought stress and nitrogen deficiency in wheat [38, 41]. These findings indicate that circRNAs are present in different plant species and may play important roles in the regulation of growth and development, and stress response.

Wheat (Triticum aestivum L.) is one of the most important food crops in the world. In wheat, circRNAs have been shown to be involved in response to drought stress. Wang et al. isolated 88 circRNAs and found that 62 circRNAs were differentially expressed under drought stress conditions compared with control [38]. In a recent study, Ren et al. identified six circRNAs involved in the common response to nitrogen deficiency stress and 23 circRNAs involved in the regulation of low nitrogen-promoted root growth in wheat [41]. These studies show that circRNAs may play important roles in responding to abiotic stresses. However, it is unclear whether circRNAs are involved in the regulation of wheat root growth under normal growth conditions. To explore this question, the root circRNA expression profiles of three wheat varieties (including two long root varieties and one short root variety) were obtained using high-throughput sequencing technology. Differentially expressed circRNAs (DECs) were identified and further validated using real-time PCR technology. Target miRNAs of the DECs were predicted and the expression levels of these miRNAs in the roots of LH9, XN979 and YN29 were also examined. This is the first report on the identification of differentially expressed circRNAs in wheat varieties with contrasting root length. The results are helpful for further investigate of the potential roles of circRNAs in regulating wheat root length.

Methods

Plant materials

XN979, LH9 and YN29 are three Chinese wheat varieties. The roots of LH9 and XN979 are significantly longer than those of YN29. Here, XN979, LH9 and YN29 were selected as materials to identify circRNAs involved in the regulation of wheat root length.

Plant growth conditions and evaluation of root phenotype

Seed sterilization, germination and the growth conditions of wheat plants were conducted as previously described by Ren et al. [42]. Plants were randomly placed and grown in a greenhouse with six replications each. The maximum root length (MRL) and total root length (TRL) of XN979, LH9 and YN29 were analyzed by using WinRHIZO software (Regent Instruments, Canada) after 15 days of transferring. The developmental stages of wheat plants were Zadoks growth scale 13 [43]. The roots of XN979, LH9 and YN29 were fast-frozen using liquid nitrogen for RNA extraction.

Libraries construction and sequencing

Root total RNA was extracted using Trizol reagent following the manufacturer’s procedure. The concentration and purity of total RNA were measured by NanoDrop ND-1000 spectrophotometer (NanoDrop Technologies, Wilmington, DE, USA). Approximately 10 µg total RNA was used to deplete the ribosomal RNA (rRNA) following the instructions of the Epicentre Ribo-Zero Gold Kit (Illumina, San Diego, USA). The remaining RNAs were used as templates for the construction of cDNA libraries following the instructions of the RNA-Seq sample preparation kit (Illumina, San Diego, USA) [41]. Three biological replicates were analyzed. The samples were named as XN979-1, -2, -3; LH9-1, -2, -3; YN29-1, -2, -3, respectively. Paired-end sequencing was performed on an Illumina Hiseq 2500 platform (Hangzhou Shangyi biotechnology company, Hangzhou, China).

Identification of circRNAs

The obtained clean reads were tempted to align with wheat reference genome (Triticum aestivum TGACv1.0) using the bowtie2 (bowtie2-2.2.2) alignment method. The reads of the linear RNA can be mapped appropriately to wheat reference genome, while the reads at the loop forming junctions of circRNAs cannot be directly aligned to the wheat reference genome [41]. Then the find_circ software was employed to detect head-to-tail splicing (back-spliced) of the remaining unmapped RNA-seq reads (default setting). The detected back-spliced reads were further filtered to predict circRNAs following the recommended setting rules (GU/AG appears on both sides of the splice site, clear breakpoint can be detected; ≤ 2 mismatches; The length of the circRNA junctions ≤ 100 kb) [41, 44].

Differential expression analysis of circRNAs

Differential expression analysis of circRNAs between different varieties was performed using the DEseq R package [41]. Only the circRNAs with p values ≤ 0.05 and |log2 (foldchange)| ≥ 1 were defined as DECs [45, 46]. psRNATarget software was used to predict circRNA-miRNA interactions of the DECs [47].

Real-time PCR analysis

To confirm and quantify the predicted circRNAs, divergent primers of five randomly selected DECs were designed based on the flanking sequences of the head-to-tail splicing site of each circRNA (Additional file 1: Table S1). The cDNA samples were used as templates and mixed with primers and SYBR Green PCR Real Master Mix (Tiangen, China) for real-time PCR analyses. In the quantitative analysis of selected miRNAs, miRcute Plus miRNA First-Strand cDNA Synthesis Kit (Tiangen, China) was used for cDNA synthesis and miRcute Plus miRNA qPCR Detection Kit (Tiangen) was used for quantitative analysis according to the kit instructions. The reverse primers were provided in the miRcute Plus miRNA qPCR Detection Kit and the forward ones were designed according to the instructions of the kit (Additional file 1: Table S1). Real-time PCRs were performed on a Thermal Cycler CFX96 Real-Time System (Bio-Rad, USA). The program settings are as follows: 95 °C for 5 min, then 40 cycles of 95 °C for 15 s, 60 °C for 15 s, 72 °C for 15 s. TaActin was used as an internal reference gene to normalize the expression levels of the investigated DECs and miRNAs. SPSS 21.0 software was used to analyze the statistical significance of the data.

Results

Root phenotypes of different wheat varieties

The root phenotypes of LH9, XN979 and YN29 are shown in Fig. 1. There are significant differences between the long root varieties (LH9, XN979) and the short root variety (YN29). The MRL of LH9 and XN979 were 13.7 cm and 9.2 cm longer than that of YN29 (Fig. 1a), and the TRL of LH9 and XN979 were 302.4 cm and 211.9 cm longer than that of YN29 respectively (Fig. 1b). Overall, the roots of LH9 and XN979 were significantly larger and deeper than those of YN29 (Fig. 1c).

Fig. 1
figure 1

The root length and morphology of LH9, XN979 and YN29. The maximum root length (a), total root length (b) and root morphology (c) of LH9, XN979 and YN29. Bar = 5 cm

Full size image

CircRNAs identification

To determine circRNAs involved in the regulation of root length, nine cDNA libraries from wheat roots of LH9, XN979 and YN29 were constructed and sequenced. Over 100 million raw reads were generated in each library (among 100.9 to 104.7 million raw reads in each library). About 76% to 85% reads were successfully mapped to wheat genome (TGAC v1.0). Find_circ software was employed to detect head-to-tail splicing (back-spliced) of the remaining 15–24% unmapped RNA-seq reads. The number of circRNAs identified in each sample ranged from 285 to 478. Among them, more than 70% of circRNAs are exonic circRNAs (Fig. 2a). The proportion of intronic circRNAs is the lowest, no more than 6.5% in each sample. The proportion of intergenic circRNAs is between the two, accounting for 18.7–24.9% (Fig. 2a).

Fig. 2
figure 2

Statistical analysis of identified circRNAs and differentially expressed circRNAs (DECs) in the roots of LH9, XN979 and YN29. a The number of exonic circRNAs, intergenic circRNAs, and intronic circRNAs in each sequenced sample. b Venn diagram analysis of DECs in the YN29-LH9 and YN29-XN979 comparisons. c The number of up-regulated and down-regulated DECs in each comparison

Full size image

Identification of differentially expressed circRNAs

We compared the circRNAs expression profiles in the long root varieties (XN979 and LH9) and the short root variety (YN29), and obtained DECs between them. In the YN29-LH9 comparison, 33 circRNAs were identified as DECs (p value ≤ 0.05 along with |log2 (foldchange)| ≥ 1) (Fig. 2b; Additional file 2: Table S2). Among them, 18 circRNAs were up-regulated and 15 were down-regulated in YN29 compared with LH9 (Fig. 2c). We totally identified 22 DECs in the one to one comparison between YN29 and XN979 (Fig. 2b; Additional file 3: Table S3). Among them, 15 circRNAs showed up-regulation and seven circRNAs exhibited down-regulation in YN29 compared with XN979 (Fig. 2c). Among the circRNAs identified above, ten circRNAs existed in both the YN29-LH9 comparison and the YN29-XN979 comparison (Fig. 2b; Table 1). Twenty-three DECs were specifically found in the YN29-LH9 comparison and 12 DECs were specifically found in the YN29-XN979 comparison (Fig. 2b). Since the roots of both LH9 and XN979 are significantly larger and deeper than those of YN29 (Fig. 1c), the ten DECs coexisting in the two comparisons are highly likely to be involved in the regulation of wheat root length.

Table 1 CircRNAs involved in the regulation of wheat root length
Full size table

Real-time PCR verification of DECs

Five randomly selected DECs were selected for data verification of RNA-seq by using real-time PCR technology. Results showed that the expression levels of the five DECs detected using real-time PCR technology matched well with the results of RNA-seq (Fig. 3), indicating that the results of RNA-seq are reliable.

Fig. 3
figure 3

Relative expression analysis of differentially expressed circRNAs (DECs). a The results of relative expression analysis of circRNAs by using real-time PCR technology (*p < 0.05). b The results of relative expression analysis of circRNAs by using RNA-sequence technology

Full size image

The regulation of miRNAs levels by circRNAs

It has been reported that circRNAs can act as ceRNAs (competitive endogenous RNA) of miRNAs, inhibiting their binding to mRNA molecules, thereby regulating gene expression [21, 44]. To detect whether the ten DECs coexisted in both comparisons can perform the function, psRNATarget software was employed to predict potential miRNA binding sites of the circRNAs. Results showed that three out of the ten DECs were predicted to have two to nine corresponding miRNAs binding sites (Table 1). CircRNA23 have putative miRNA binding sites for tae-miR1134 and tae-miR5085. CircRNA3645 have putative miRNA binding sites for tae-miR1133, tae-miR1137a and tae-miR9773. According to the psRNATarget-based prediction, circRNA2473 can bind to nine different miRNAs. To investigate whether these circRNAs regulate the expression of the corresponding miRNAs, eight potential circRNA-binding miRNAs were selected to check their expression levels in the long root plants (LH9 and XN979) and the short root plants (YN29). Results showed that six of them exhibited significant differences between the long root plants (LH9 and XN979) and the short root plants (YN29) (Fig. 4).

Fig. 4
figure 4

Relative expression analysis of potential binding miRNAs of the DECs in the roots of LH9, XN979 and YN29. Significant difference at p < 0.05 is indicated by different letters above the columns

Full size image

Discussion

Deeper rooting is beneficial for plant production and survival under water and nitrogen deficient conditions [48,49,50]. Therefore, genetic improvement of root traits is essential for improving crop nutrient and water use efficiencies. Recent studies have shown that circRNAs are widespread in plants and play important roles in regulating plant growth and development, and stress tolerance [23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41]. However, whether circRNAs, as a new type of non-coding RNA molecules, participate in the regulation of root growth is not clear. Here, an integrated comparative root transcriptomic study of long root plants (LH9 and XN979) and short root plants (YN29) was conducted to explore this question. The number of circRNAs identified in each sample ranged from 285 to 478. More than 70% circRNAs belong to exonic circRNAs category. The intronic circRNAs only account for 3.9–6.5% of the total number of circRNAs in different samples (Fig. 2a). The proportion of intronic circRNAs is similar to the results of previous reports in Arabidopsis (3.8% intronic circRNAs), tomato (3.6%) and wheat (less than 6.5%) [16, 19, 41]. In the YN29-LH9 comparison and YN29-XN979 comparison, 33 and 22 DECs were identified, respectively (Fig. 2b). In total, ten circRNAs existed in both the YN29-LH9 comparison and the YN29-XN979 comparison. As has been mentioned above, both the roots of LH9 and XN979 were significantly larger and deeper than those of YN29 (Fig. 1c), therefore, the ten coexisted DECs in both comparisons should play more critical roles in the regulation of root length.

Previous studies in animal and human have confirmed that circRNAs can act as miRNA sponges to capture miRNAs from their target genes via ceRNA (competing endogenous RNAs) networks [51, 52]. Similar reports have recently been made in plants [26, 34, 38, 41]. To reveal whether the ten DECs can target miRNAs and participate in transcriptional regulation of genes, psRNATarget software was used to identify potential miRNA binding sites. As expected, three of the ten circRNAs had two to nine putative miRNA-binding sites (Table 1), which was consist with another study in wheat [38]. Moreover, the results of real-time PCR analysis showed that the expression levels of six out of the eight investigated potential binding miRNAs exhibited significant differences between the long root plants (LH9 and XN979) and the short root plants (YN29) (Fig. 4). Therefore, the expression levels of these miRNAs were correlated with their corresponding circRNAs. These results indicate that the three circRNAs may bind to the target miRNAs and modulate the transcript levels of these miRNAs. Interestingly, it has been reported that some putative target miRNAs are involved in the regulation of plant growth and development [53,54,55,56,57]. For example, miR164 has been reported controlling root development in Arabidopsis, maize and potato [53,54,55]. In Arabidopsis, one of miR164 target genes is NAC1, which transduces auxin signals for lateral root emergence [58, 59]. In mir164 mutants, plants express less miR164 and more NAC1 mRNA, and produced more lateral roots. Moreover, the mutant phenotypes can be complemented by expression miR164a and miR164b genomic sequences [53]. Besides, the target genes of the potential circRNA binding miRNAs such as miR1122, miR1125, miR1134 and miR1133 are also involved in plant growth, development, metabolism and stress response [56, 57]. Therefore, these circRNAs may regulate wheat root length by modulating target miRNAs levels. We also noticed that different circRNAs may have the same miRNA binding site, and one circRNA may have several corresponding miRNA binding sites (Table 1). Since the mechanism of how circRNAs regulate their target miRNA are not clear in plants so far and these miRNA binding sites are predicted by psRNATarget software, whether circRNA can bind to each of them still requires further experimental validation. Furthermore, how do the circRNAs that do not have miRNAs binding sites work? What are the mechanisms by which different circRNAs coordinately regulate the length of wheat roots? There is still a lot of work to be done to clarify the regulatory roles of these circRNAs.

Conclusions

Our study revealed that the expression levels of some circRNAs in roots exhibited significant differences between the long root plants (LH9 and XN979) and the short root plants (YN29). Ten DECs involved in the regulation of wheat root length were successfully identified, three of which have potential miRNAs binding sites. The expression levels of putative circRNA-binding miRNAs were correlated with their corresponding circRNAs. These results provide new clues for investigating the functions of circRNAs in the regulation of wheat root length.

Abbreviations

BR:

brassinosteroid

ceRNA:

competing endogenous RNA

circRNAs:

circular RNAs

DECs:

differentially expressed circRNAs

MRL:

maximum root length

rRNA:

ribosomal RNAs

TRL:

total root length

References

  1. de Dorlodot S, Forster B, Pagès L, Price A, Tuberosa R, Draye X. Root system architecture: opportunities and constraints for genetic improvement of crops. Trends Plant Sci. 2007;12:474–81.

    Article  Google Scholar 

  2. Lynch J. Roots of the second green revolution. Aust J Bot. 2007;55:493–512.

    Article  Google Scholar 

  3. Den Herder G, Van Isterdael G, Beeckman T, De Smet I. The roots of a new green revolution. Trends Plant Sci. 2010;15:600–7.

    Article  Google Scholar 

  4. White PJ, George TS, Gregory PJ, Bengough AG, Hallett PD, McKenzie BM. Matching roots to their environment. Ann Bot. 2013;112:207–22.

    Article  CAS  Google Scholar 

  5. Yu P, Li XX, White PJ, Li CJ. A large and deep root system underlies high nitrogen-use efficiency in maize production. PLoS ONE. 2015;10:e0126293.

    Article  Google Scholar 

  6. Pacifici E, Polverari L, Sabatini S. Plant hormone crosstalk: the pivot of root growth. J Exp Bot. 2015;66:1113–21.

    Article  CAS  Google Scholar 

  7. Dello Ioio R, Nakamura K, Moubayidin L, Perilli S, Taniguchi M, Morita MT, Aoyama T, Costantino P, Sabatini S. A genetic frame work for the control of cell division and differentiation in the root meristem. Science. 2008;322:1380–4.

    Article  Google Scholar 

  8. Vragovic K, Sela A, Friedlander-Shani L, Fridman Y, Hacham Y, Holland N, Bartom E, Mockler TC, Savaldi-Goldstein S. Translatome analyses capture of opposing tissue-specific brassinosteroid signals orchestrating root meristem differentiation. Proc Natl Acad Sci USA. 2015;112:923–8.

    Article  CAS  Google Scholar 

  9. Li J, Nagpal P, Vitart V, McMorris TC, Chory J. A role for brassinosteroids in light-dependent development of Arabidopsis. Science. 1996;272:398–401.

    Article  CAS  Google Scholar 

  10. Müssig C, Shin GH, Altmann T. Brassinosteroids promote root growth in Arabidopsis. Plant Physiol. 2003;133:1261–71.

    Article  Google Scholar 

  11. Roddick JG, Rijnenberg AL, Ikekawa N. Developmental effects of 24-epibrassinolide in excised roots of tomato grown in vitro. Physiol Plant. 1993;87:453–8.

    Article  CAS  Google Scholar 

  12. Clouse SD, Langford M, McMorris TC. A brassinosteroid-insensitive mutant in Arabidopsis thaliana exhibits multiple defects in growth and development. Plant Physiol. 1996;111:671–8.

    Article  CAS  Google Scholar 

  13. He X, Fang JJ, Li JJ, Qu BY, Ren YZ, Ma WY, Zhao X, Li B, Wang DW, Li ZS, et al. A genotypic difference in primary root length is associated with the inhibitory role of transforming growth factor-beta receptor-interacting protein-1 on root meristem size in wheat. Plant J. 2014;77:931–43.

    Article  CAS  Google Scholar 

  14. De Tullio MC, Jiang K, Feldman LJ. Redox regulation of root apical meristem organization: connecting root development to its environment. Plant Physiol Biochem. 2010;48:328–36.

    Article  Google Scholar 

  15. Xia XJ, Chen Z, Yu JQ. ROS mediate brassinosteroids-induced plant stress responses. Plant Signal Behav. 2010;5:532–4.

    Article  CAS  Google Scholar 

  16. Lasda E, Parker R. Circular RNAs: diversity of form and function. RNA. 2014;20:1829–42.

    Article  CAS  Google Scholar 

  17. Chen LL, Yang L. Regulation of circRNA biogenesis. RNA Biol. 2015;12:381–8.

    Article  Google Scholar 

  18. Salzman J, Gawad C, Wang PL, Lacayo N, Brown PO. Circular RNAs are the predominant transcript isoform from hundreds of human genes in diverse cell types. PLoS ONE. 2012;7:e30733.

    Article  CAS  Google Scholar 

  19. Jeck WR, Sharpless NE. Detecting and characterizing circular RNAs. Nat Biotechnol. 2014;32:453–61.

    Article  CAS  Google Scholar 

  20. Westholm JO, Miura P, Olson S, Shenker S, Joseph B, Sanfilippo P, Celniker SE, Graveley BR, Lai EC. Genome-wide analysis of Drosophila circular RNAs reveals their structural and sequence properties and age-dependent neural accumulation. Cell Rep. 2014;9:1966–80.

    Article  CAS  Google Scholar 

  21. Fan XY, Zhang XN, Wu XL, Guo HS, Hu YQ, Tang FC, Huang YY. Single-cell RNA-seq transcriptome analysis of linear and circular RNAs in mouse preimplantation embryos. Genome Biol. 2015;16:148.

    Article  Google Scholar 

  22. Sablok G, Zhao H, Sun X. Plant Circular RNAs (circRNAs): transcriptional regulation beyond miRNAs in plants. Mol. Plant. 2016;9:192–4.

    Article  CAS  Google Scholar 

  23. Pan T, Sun XQ, Liu YX, Li H, Deng GB, Lin HH, Wang SH. Heat stress alters genome-wide profiles of circular RNAs in Arabidopsis. Plant Mol Biol. 2018;96:217–29.

    Article  CAS  Google Scholar 

  24. Liu TF, Zhang L, Chen G, Shi TL. Identifying and characterizing the circular RNAs during the lifespan of Arabidopsis leaves. Front Plant Sci. 2017;8:1278.

    Article  Google Scholar 

  25. Ye CY, Chen L, Liu C, Zhu QH, Fan L. Widespread noncoding circular RNAs in plants. New Phytol. 2015;208:88–95.

    Article  CAS  Google Scholar 

  26. Chen G, Cui JW, Wang L, Zhu YF, Lu ZG, Jin B. Genome-wide identification of circular RNAs in Arabidopsis thaliana. Front Plant Sci. 2017;8:1678.

    Article  Google Scholar 

  27. Meng X, Zhang P, Chen Q, Wang J, Chen M. Identification and characterization of ncRNA-associated ceRNA networks in Arabidopsis leaf development. BMC Genomics. 2018;19:607.

    Article  Google Scholar 

  28. Zhao W, Cheng YH, Zhang C, You QB, Shen XJ, Guo W, Jiao YQ. Genome-wide identification and characterization of circular RNAs by high throughput sequencing in soybean. Sci. Rep. 2017;7:5636.

    Article  Google Scholar 

  29. Chen L, Ding X, Zhang H, He T, Li Y, Wang T, Li X, Jin L, Song Q, Yang S, et al. Comparative analysis of circular RNAs between soybean cytoplasmic male-sterile line NJCMS1A and its maintainer NJCMS1B by high-throughput sequencing. BMC Genomics. 2018;19:663.

    Article  Google Scholar 

  30. Lu TT, Cui LL, Zhou Y, Zhu CR, Fan DL, Gong H, Zhao Q, Zhou CC, Zhao Z, Lu DF, et al. Transcriptome-wide investigation of circular RNAs in rice. RNA. 2015;21:2076–87.

    Article  CAS  Google Scholar 

  31. 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. 2017;217:1292–306.

    Article  Google Scholar 

  32. Ghorbani A, Izadpanah K, Peters JR, Dietzgen RG, Mitter N. Detection and profiling of circular RNAs in uninfected and maize Iranian mosaic virus-infected maize. Plant Sci. 2018;274:402–9.

    Article  CAS  Google Scholar 

  33. Tang B, Hao Z, Zhu Y, Zhang H, Li G. Genome-wide identification and functional analysis of circRNAs in Zea mays. PLoS ONE. 2018;13:e0202375.

    Article  CAS  Google Scholar 

  34. Zuo J, Wang Q, Zhu B, Luo Y, Gao L. Deciphering the roles of circRNAs on chilling injury in tomato. Biochem Biophys Res Commun. 2016;479:132–8.

    Article  CAS  Google Scholar 

  35. Wang J, Yang Y, Jin L, Ling X, Liu T, Chen T, Ji Y, Yu W, Zhang B. Re-analysis of long non-coding RNAs and prediction of circRNAs reveal their novel roles in susceptible tomato following TYLCV infection. BMC Plant Biol. 2018;18:104.

    Article  Google Scholar 

  36. 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.

    Article  Google Scholar 

  37. Tong W, Yu J, Hou Y, 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:1417–29.

    Article  CAS  Google Scholar 

  38. Zhao T, Wang L, Li S, Xu M, Guan X, Zhou B. Characterization of conserved circular RNA in polyploid Gossypium species and their ancestors. FEBS Lett. 2017;591:3660–9.

    Article  CAS  Google Scholar 

  39. Xiang L, Cai C, Cheng J, Wang L, Wu C, Shi Y, Luo J, He L, Deng Y, Zhang X, et al. Identification of circularRNAs and their targets in Gossypium under Verticillium wilt stress based on RNA-seq. PeerJ. 2018;6:e4500.

    Article  Google Scholar 

  40. Wang YX, Yang M, Wei SM, Qin FJ, Zhao HJ, Suo B. Identification of circular RNAs and their targets in leaves of Triticum aestivum L. under dehydration stress. Front. Plant Sci. 2017;7:2024.

    PubMed  PubMed Central  Google Scholar 

  41. Ren Y, Yue H, Li L, Xu Y, Wang Z, Xin Z, Lin T. Identification and characterization of circRNAs involved in the regulation of low nitrogen-promotedroot growth in hexaploid wheat. Biol Res. 2018;51:43.

    Article  Google Scholar 

  42. Ren YZ, He X, Liu DC, Li JJ, Zhao XQ, Li B, Tong YP, Zhang AM, Li ZS. Major quantitative trait loci for seminal root morphology of wheat seedlings. Mol Breeding. 2012;30:139–48.

    Article  Google Scholar 

  43. Zadoks JC, Chang TT, Konzak CF. Decimal code for the growth stages of cereals. Weed Res. 1974;4:415–21.

    Article  Google Scholar 

  44. 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:333–8.

    Article  CAS  Google Scholar 

  45. Zhou L, Chen JH, Li ZZ, Li XX, Hu XD, Huang Y, Zhao XK, Liang CZ, Wang Y, Sun L, et al. Integrated profiling of microRNAs and mRNAs: microRNAs located on Xq27.3 associate with clear cell renal cell carcinoma. PloS ONE. 2010;5:e15224.

    Article  CAS  Google Scholar 

  46. Anders S, Huber W. Differential expression analysis for sequence count data. Genome Biol. 2010;11:R106.

    Article  CAS  Google Scholar 

  47. Dai X, Zhao PX. psRNATarget: a plant small RNA target analysis server. Nucleic Acids Res. 2011;39:W155.

    Article  CAS  Google Scholar 

  48. Uga Y, Sugimoto K, Ogawa S, Rane J, Ishitani M, Hara N, Kitomi Y, Inukai Y, Ono K, Kanno N, et al. Control of root system architecture by DEEPER ROOTING 1 increases rice yield under drought conditions. Nat Genet. 2013;45:1097–102.

    Article  CAS  Google Scholar 

  49. Comas LH, Becker SR, Cruz VMV, Byrne PF, Dierig DA. Root traits contributing to plant productivity under drought. Front Plant Sci. 2013;4:442.

    Article  Google Scholar 

  50. Koevoets IT, Venema JH, Elzenga JT, Testerink C. Roots withstanding their environment: exploiting root system architecture responses to abiotic stress to improve crop tolerance. Front Plant Sci. 2016;31:1335.

    Google Scholar 

  51. 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:384–8.

    Article  CAS  Google Scholar 

  52. Zhang S, Zhu D, Li H, Li H, Feng C, Zhang W. Characterization of circRNA-associated-ceRNA networks in a senescence-accelerated mouse prone 8 brain. Mol Ther. 2017;25:2053–61.

    Article  CAS  Google Scholar 

  53. Guo HS, Xie Q, Fei JF, Chua NH. MicroRNA directs mRNA cleavage of the transcription factor NAC1 to downregulate auxin signals for Arabidopsis lateral root development. Plant Cell. 2005;17:1376–86.

    Article  CAS  Google Scholar 

  54. Li J, Guo GH, Guo WW, Guo GG, Tong D, Ni ZF, Sun QX, Yao YY. miRNA164-directed cleavage of ZmNAC1 confers lateral root development in maize (Zea mays L.). BMC Plant Biol. 2012;12:220.

    Article  CAS  Google Scholar 

  55. Zhang L, Yao L, Zhang N, Yang J, Zhu X, Tang X, Calderón-Urrea A, Si H. Lateral root development in potato is mediated by stu-mi164 regulation of NAC transcription factor. Front Plant Sci. 2018;9:383.

    Article  Google Scholar 

  56. Pandey B, Gupta OP, Pandey DM, Sharma I, Sharma P. Identification of new stress-induced microRNA and their targets in wheat using computational approach. Plant Signal Behav. 2013;8:e23932.

    Article  Google Scholar 

  57. Zhou L, Liu Y, Liu Z, Kong D, Duan M, Luo L. Genome-wide identification and analysis of drought-responsive microRNAs in Oryza sativa. J Exp Bot. 2010;61:4157–68.

    Article  CAS  Google Scholar 

  58. Rhoades MW, Reinhart BJ, Lim LP, Burge CB, Bartel B, Bartel DP. Prediction of plant microRNA targets. Cell. 2002;110:513–20.

    Article  CAS  Google Scholar 

  59. Xie Q, Guo HS, Dallman G, Fang S, Weissman AM, Chua NH. SINAT5 promotes ubiquitin-related degradation of NAC1 to attenuate auxin signals. Nature. 2002;419:167–70.

    Article  CAS  Google Scholar 

Download references

Authors’ contributions

YX wrote the paper. YR performed most of the experiments; YX took part in partial work of this research; YR, DC and TL conceived and designed the experiments; DC gave many helpful suggestions on the writing. All authors read and approved the final manuscript.

Acknowledgements

We thank Huifang Yue for her help in root phenotypic analysis and data processing.

Competing interests

The authors declare that they have no competing interests.

Availability of data and materials

The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.

Consent for publication

Not applicable.

Ethics approval and consent to participate

Not applicable.

Funding

This work was supported by the National Natural Science Foundation of China (31572197), the National Natural Science Foundation of China (31201678), the Natural Science Foundation of Henan province (162300410133) and the State Key Laboratory Program (PCCE-KF-2017-04).

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Author information

Author notes

  1. Yanhua Xu and Yongzhe Ren contributed equally to this work

Authors and Affiliations

  1. College of Agronomy, Henan Agricultural University, Zhengzhou, 450002, China

    Yanhua Xu, Yongzhe Ren, Tongbao Lin & Dangqun Cui

  2. State Key Laboratory of Wheat and Maize Crop Science, Henan Agricultural University, Zhengzhou, 450002, China

    Yanhua Xu, Yongzhe Ren, Tongbao Lin & Dangqun Cui

  3. Collaborative Innovation Center of Henan Grain Crops, Henan Agricultural University, Zhengzhou, 450002, China

    Yanhua Xu, Yongzhe Ren, Tongbao Lin & Dangqun Cui

  4. College of life science, Shangqiu Normal University, Shangqiu, 476000, China

    Yanhua Xu

Authors
  1. Yanhua Xu
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Yongzhe Ren
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Tongbao Lin
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Dangqun Cui
    View author publications

    You can also search for this author in PubMed Google Scholar

Corresponding authors

Correspondence to Yongzhe Ren or Dangqun Cui.

Additional files

Additional file 1: Table S1.

Primers used for real-time PCR.

Additional file 2: Table S2.

List of differentially expressed circRNAs (DECs) in the YN29-LH9 comparison.

Additional file 3: Table S3.

List of differentially expressed circRNAs (DECs) in the YN29-XN979 comparison.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Xu, Y., Ren, Y., Lin, T. et al. Identification and characterization of CircRNAs involved in the regulation of wheat root length. Biol Res 52, 19 (2019). https://doi.org/10.1186/s40659-019-0228-5

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s40659-019-0228-5

Keywords

Biological Research

ISSN: 0717-6287

Contact us