Skip to main content

High throughput sequencing of two celery varieties small RNAs identifies microRNAs involved in temperature stress response



MicroRNAs (miRNAs) are small, non-coding RNAs of 20 to 24 nucleotides that regulate gene expression and responses to biotic and abiotic stress. Till now, no reports have previously been published concerning miRNAs in celery.


Two small RNAs libraries were constructed from two celery varieties, ‘Jinnan Shiqin’ and ‘Ventura’, and characterized by deep sequencing. A total of 431 (418 known and 13 novel) and 346 (341 known and five novel) miRNAs were identified in celery varieties ‘Jinnan Shiqin’ and ‘Ventura’, respectively. Potential miRNA-target genes were predicted and annotated by screening diverse protein databases, including Gene Ontology, Cluster of Orthologous Groups and Kyoto Encyclopedia of Genes and Genomes. Significant differential expression between the two varieties was seen for 221 miRNAs. qRT-PCR was used to analyze the abundance of six miRNAs under cold and heat stress conditions. The results showed that miRNAs may have important functions in controlling temperature stress in celery.


A large number of miRNAs were identified in celery, and their target genes, functional annotations, and gene expression patterns have been explored.

These findings provide the first information on celery miRNAs and enhance understanding of celery miRNA regulatory mechanisms under extreme temperature stress.


MicroRNAs (miRNAs) are small, non-coding RNAs found in animals and plants, and mainly function in regulating gene expression at the post-transcriptional level [1]. miRNAs are highly conserved in eukaryotes, and are an important component in the evolution of genetic regulation [2, 3]. Mature miRNAs are generated from primary miRNAs (pri-miRNAs) via two steps. The pri-miRNAs are cut by RNA polymerase II into precursor miRNAs (pre-miRNAs) that have a stem-loop structure and are 70 to 100 nucleotides (nt) in length. The pre-miRNAs are then cleaved into mature miRNAs of 19 to 23 nt by Dicer-Like1 in the cytoplasm [4].

The first two characterized miRNAs, lin-4 and let-7, were identified as regulators of the juvenile-to-adult phase transition in Caenorhabditis elegans[5, 6]. The latest miRBase database (miRBase19) contains 21,264 entries representing hairpin pre-miRNAs, expressing 25,141 mature miRNA products in 193 species. Numerous studies have revealed that miRNAs are involved in diverse biological and metabolic processes, such as regulation of cell growth and development, abiotic stress response, pathogen defense and gene translational repression [710]. Levels of miR156, miR167, and miR164 increase during virus intrusion [11, 12]. Furthermore, miR172 and miR159 affect flowering time regulation in plants [13, 14], whereas miR160 has important roles in the regulation of plant development and hormonal signals [15]. In cold conditions, miR393 over-expression is induced [16].

Celery (Apium graveolens L.) is an annual or perennial herb belonging to the Apiaceae family. Although celery originated from the Mediterranean and the Middle East, it is now cultivated worldwide. Celery is rich in carotene, vitamins, carbohydrates, and volatile aroma compounds, and has excellent medicinal properties to regulate the digestive system and blood lipids [1719]. However, celery is one of the most common allergenic foods in many European countries [20]. Numerous allergenic compounds contained in celery tissue can induce allergic symptoms [21]. miRNAs have significant biological functions in metabolic and immune responses [8, 22], and exploring the miRNAs present in celery will provide a foundation for the study of their functions. Till now, there have been no previously published reports concerning miRNA s in Apiaceae species.

Extreme temperature is one of the major limiting factors on celery growth and yield. Low temperatures can lead to early bolting, and high temperatures can induce various diseases. To investigate the roles of celery miRNAs under temperature stress, two varieties, ‘Jinnan Shiqin’ and ‘Ventura’, were selected for small RNA and transcriptome sequencing. ‘Jinnan Shiqin’ was bred by the Institute of Tianjin City, China, while ‘Ventura’ originated in the United States and was introduced to China. These two varieties have similar phenotypes, although they originated from different geographical areas. Here, the significant differential expression miRNAs were detected between two varieties, and the results provide useful information for miRNA that can response to temperature stress.


Sequence analysis of small RNAs

To identify the low-abundance candidate miRNAs in the ‘Jinnan Shiqin’ and ‘Ventura’ celery varieties, Solexa sequencing was performed on the small RNAs libraries. A total of 6,951,840 and 6,057,911 reads were generated from ‘Jinnan Shiqin’ and ‘Ventura’ libraries, respectively. After filtering the reads for those of low quality, 3’ insert null, poly(A), length < 16 or length > 30, the majority of small RNAs were 21 to 24 nt in length (Figure 1). The length distributions of small RNAs were similar for both varieties. The 24 nt small RNAs were the most abundant, representing 35.75% and 24.50% of small RNAs in ‘Jinnan Shiqin’ and ‘Ventura’, respectively. This result is consistent with reports for other plants, such as Arabidopsis thaliana, trifoliate oranges, peanuts and olives [2326]. Some studies have reported that 24 nt siRNAs are involved in heterochromatin modification, particularly in genomes with high repetitive sequence content [27, 28].

Figure 1

Length distributions of small RNAs in two celery varieties.

The small RNAs were further classified into different categories by performing BLAST searches against Rfam. Non-coding RNAs included miRNA, rRNA, tRNA, snoRNA, snRNA, and other unannotated RNAs. More categories of small RNAs were present in ‘Jinnan Shiqin’ than in ‘Ventura’ (Table 1).

Table 1 Distribution of small RNAs among different categories in ‘Jinnan Shiqin’ and ‘Ventura’

Identification of conserved miRNAs

To identify conserved miRNAs in the two celery varieties, small RNA sequences were compared with currently known mature plant miRNAs in the miRBase. Total of 418 unique miRNA sequences belonging to 28 miRNA families were identified in ‘Jinnan Shiqin’, and total of 341 unique miRNA sequences belonging to 34 miRNA families were identified in ‘Ventura’ (Figure 2). Most of the identified miRNAs showed high sequence similarities to miRNAs in other plants. The distributions of miRNA families were similar between the two libraries. Several families, such as miR156, miR159, miR166, miR167-1 and miR396, were relatively abundant, whereas some families were not. The length of miRNAs and nucleotide preference distributions are shown in Figure 3. The majority of miRNAs tended to start with 5’-U and not with 5’-G, which is consistent with typical miRNA sequence patterns.

Figure 2

Numbers of identical miRNA members in each family in ‘Jinnan Shiqin’ and ‘Ventura’.

Figure 3

First nucleotide bias of miRNAs in ‘Jinnan Shiqin’ and ‘Ventura’.

Identification of novel miRNAs

To identify novel miRNA sequences, all unannotated small RNAs were searched against transcriptome sequence data (SRX326597). After searching for potential pre-miRNAs and predicting their hairpin-like structures, 13 and five unique sequences were identified as novel miRNAs in ‘Jinnan Shiqin’ and ‘Ventura’, respectively (Table 2 and Table 3). Only one novel miRNA, celery-miR-27, was identified in both celery varieties. The novel miRNA sequences were 20 to 24 nt in length, and 24 nt reads were the most abundant. The range in length of pre-miRNAs was 77 to 102 nt in ‘Jinnan Shiqin’, with an average of 94 nt. The average minimum free energy values in ‘Jinnan Shiqin’ and ‘Ventura’ were −23.18 and −31.12 kcal/mol, respectively, which are much higher than values for other plant pre-miRNAs (−59.5 and −71.0 kcal/mol in Arabidopsis and rice, respectively) [29].

Table 2 Novel miRNA candidates in ‘Jinnan Shiqin’
Table 3 Novel miRNA candidates in ‘Ventura’

Prediction and annotation of miRNA target genes

To better understand the functions of the identified miRNAs, putative target genes were predicted using the psRNA Target program. A total of 503 and 408 potential target genes were identified in ‘Jinnan Shiqin’ and ‘Ventura’, respectively. To evaluate the potential functions of these miRNA target genes, GO analysis was used [30]. The miRNA target genes were categorized according to biological process, cellular component, and molecular function (Figure 4). Fifty-six genes were categorized as cellular components. Based on molecular function, 111 genes were classified into nine categories, among which the most over-represented GO terms were binding and catalytic activities. Eight different biological processes were found, among which the three most frequent terms were metabolic process, cellular process, and biological regulation. To further evaluate the completeness of the transcriptome and the effectiveness of the annotation process, the annotated sequences were searched for genes that could be assigned to Clusters of Orthologous Groups (COG) classifications [31]. Overall, 57 sequences were assigned to COG classifications (Figure 5), and COG-annotated putative proteins were functionally classified into at least 25 molecular families. The cluster for general function prediction was the largest group, followed in descending order by translation, secondary metabolite biosynthesis, and signal transduction. The biological interpretation of these target genes was further examined using KEGG pathway analysis. Twenty-six different pathways were found, and the most frequently represented pathways included those for transcription factors (ten enzymes), plant hormone signal transduction (four enzymes), chromosomes (four enzymes), and protein kinase inhibitors (3 enzymes).

Figure 4

GO categories and distribution of miRNA targets in celery.

Figure 5

COG function classification of consensus sequences in celery.

Differential expression of miRNAs between celery varieties

The differences between the miRNA profiles for the two celery varieties are possibly related to differences in their responses to the surrounding environment. IDEG6 software was used to analyze the miRNAs of both varieties. A false discovery rate (FDR) < 0.001 and an absolute threshold value of the log2 ratio fold-change > 1 were used to determine the statistical significance of miRNA abundance. A total of 221 differentially expressed miRNAs were selected from the two libraries (Figure 6 and Additional file 1: Table S1). In ‘Jinnan Shiqin’, 55 genes were up-regulated, whereas 166 genes were down-regulated. Some families showed large differences in abundance between the two varieties. For example, mtr-miR166e and osa-miR396a were more abundant in ‘Jinnan Shiqin’, whereas sbi-miR166b and ghr-miR396b were more abundant in ‘Ventura’.

Figure 6

Expression distribution of miRNAs between ‘Jinnan Shiqin’ and ‘Ventura’.

qRT-PCR analysis of miRNAs following cold or heat stress

In this study, six miRNAs (miR160, miR164, miR168, miR394, miR395 and miR408), whose expression was previously reported to be regulated by cold or heat stress, were analyzed after temperature treatment. The changed profiles of these miRNAs are shown in Figure 7. Almost all of these miRNAs were up-regulated by cold or heat stress in both celery varieties. Aside from miR168, the abundance of the miRNAs was higher under cold stress in ‘Jinnan Shiqin’ than in ‘Ventura’, with the biggest change in miR164. After heat treatment, miR395 and miR408 were highly abundant in ‘Jinnan Shiqin’, whereas miR164, miR168, and miR394 were more abundant in ‘Ventura’. These results indicate that these six miRNAs are involved in temperature stress responses during a relatively short period of time.

Figure 7

qRT-PCR analysis of miRNAs in ‘Jinnan Shiqin’ and ‘Ventura’. CK: control check. Error bars were calculated based on three replicates.


High throughput sequencing technology has been extensively applied in small RNA research. Thousands of miRNAs and their functions have been identified in higher plants. To date, only a single report detailing transcriptome data in celery has been published, and thus, comprehensive miRNA information is not available for this plant [32]. Moreover, there are no reports concerning miRNAs in other species of Apiaceae. In the present study, miRNAs were identified and characterized from two celery varieties, namely, ‘Jinnan Shiqin’ and ‘Ventura’, which come from different geographical sources but have similar morphology. Till now, this study is the first to identify and investigate small RNAs in celery, and the results provide new information for further research into the functions, biological pathways and evolution of target genes related to temperature stresses in celery.

Over six million reads of 16 to 30 nt were obtained from each library. Based on the sequence conservation of mature miRNAs, 418 and 341 conserved miRNAs were identified in ‘Jinnan Shiqin’ and ‘Ventura’ libraries, respectively. Most conserved miRNAs showed high sequence similarities to other plants, and were distributed in 37 conserved miRNA families. Axtell et al. [33] concluded that miRNA families were highly conserved in various plants, and that these miRNAs performed analogous regulatory functions. Some highly conserved miRNA families, including miR156, miR159, miR166 and miR396, showed relatively high numbers of reads in the two celery varieties. Several studies have reported that these miRNAs have crucial roles in biotic and abiotic stress responses, plant development and fertility, and cell proliferation [14, 3436]. Other families, such as miR395 and miR812, had very low abundance in both celery varieties. These findings suggest that over long evolutionary timescales, the miRNAs evolved at different rates and have different roles in plant development. Putative targets gene were predicted to provide more information concerning the identified miRNAs, in which most target genes were involved in molecular functions and transcription pathways.

After searching for potential pre-miRNAs and predicting their hairpin-like structures, 13 and five novel miRNAs were identified in ‘Jinnan Shiqin’ and ‘Ventura’, respectively. Surprisingly, among all the newly identified miRNAs, only one (celery-miR-27) was present in both varieties. Therefore, the functions of these novel celery miRNAs should be analyzed in future studies.

Significantly different abundance was seen for 221 miRNAs between ‘Jinnan Shiqin’ and ‘Ventura’. Different miRNA genes accounted for a large proportion of all identified miRNAs. The two varieties have similar phenotypes, although they originate from different geographical areas: ‘Jinnan Shiqin’ from Asia and ‘Ventura’ from North America. This finding suggests that miRNA evolution in the same plant species is more conservative, and that variations in miRNA machinery in the two celery varieties may be critical for the differences in miRNA expression.

Considering the instability of the environment, plants are frequently challenged by various biotic and abiotic stresses, including temperature stress. Numerous studies have confirmed that small RNAs have important roles in biotic and abiotic stress responses in plants. Some miRNAs are temperature sensitive [37]. For instance, in wheat, miR393 and Ta-miR2002 were up-regulated after 0.5 h of heat treatment [38], while ten cold-regulated miRNAs were detected in Arabidopsis by microarray [39]. In castor bean, 41 miRNAs were down-regulated and four others were up-regulated under cold stress compared with normal growth conditions [40]. In the present study, qRT-PCR method was used to investigate six miRNAs under cold and heat stress conditions. All miRNAs in both celery varieties showed sensitivity to temperature change, and abundance was remarkably up-regulated. This finding is consistent with results in Arabidopsis, poplar, and wheat [38, 39, 41]. The stronger response of miRNAs to cold compared to heat stress may be related to celery being a cool-season biennial plant. The abundance of most miRNAs was higher in ‘Jinnan Shiqin’ than in ‘Ventura’, suggesting that ‘Jinnan Shiqin’ is more resistant to temperature stress. This result is consistent with the phenotypes of the two plants in relation to temperature stress. Thus, these miRNAs may have important roles in temperature stress defense. However, understanding the molecular mechanisms underlying these roles requires further study.


This is the first report to investigate small RNAs in celery; a large number of small RNAs were characterized as known and novel miRNAs. Putative target genes were predicted and then annotated by GO, COG, and KEGG databases to explore gene functions. The differential expression of miRNAs under temperature stress conditions between the two varieties suggests that the miRNA machinery in the two celery varieties may be different. Taken together, this study provides useful information for understanding the functions and regulatory mechanisms of miRNAs.


Plant materials

Celery plants (A. graveolens L. cvs. ‘Jinnan Shiqin’ and ‘Ventura’) were grown in pots containing soil:vermiculite mixture (3:1) in a controlled-environment growth chamber programmed for 16 h/8 h at 25°C/16°C for day/night and 3000 lux of light intensity. Two month-old plants were transferred to growth chambers set at 4°C or 38°C under the same light intensity and day length as above. Leaves and stems were collected 1 h after transfer to the hot or cold chambers, immediately frozen in liquid nitrogen, and then stored at −80°C until use.

RNA isolation, small RNA library development and Solexa sequencing

Total RNA was extracted using an RNA kit (RNAsimple Total RNA Kit, Tiangen, Beijing, China) according to the manufacturer’s instructions. RNA quantity and quality were examined using gel electrophoresis and a Nanodrop ND1000 spectrophotometer. Small RNAs of 18 to 30 nt were separated from total RNAs by polyacrylamide gel electrophoresis. Small RNA molecules were then ligated to Solexa adaptors at both 5’- and 3’-ends, and then converted to cDNA by RT-PCR. The purified DNA samples were sequenced using an Illumina Cluster Station and Illumina Genome Analyzer. The data were submitted to the National Center for Biotechnology Information under accession number, SRR1125031.

Prediction of conserved and novel miRNAs

After removing the impure sequences (the low quality reads, adaptor reads, and reads with length < 16 or length > 30), unique reads were queried against ribosomal RNAs (rRNAs) and transfer RNAs (tRNAs) from GenBank ( rRNA, tRNA, small nucleolar RNA (snoRNA), and small nuclear RNA (snRNA) sequences were obtained from Rfam ( Unique reads were also used for a nucleotide–nucleotide Basic Local Alignment Search Tool (BLASTn) search against the miRNA database (miRBase 16.0) to identify conserved miRNAs. To identify novel miRNAs, the Mireap program was used to obtain all candidate precursors with hairpin-like structures that were perfectly mapped by sequencing tags ( The secondary structures of putative pre-miRNAs were checked using Mfold [42]. The criteria chosen for stem-loop hairpins were described by Meyers et al. [43].

Prediction of potential miRNA target genes

Target genes of miRNAs were predicted using the psRNA Target program ( The rules used for target prediction were based on those suggested by Allen et al. [44] and Schwab et al. [45]. BLASTn hits with < 4 mismatches were chosen as candidate targets, and then nucleotide 6-frame translation-protein (blastx) was used to obtain their putative functions.

Functional annotation of the potential miRNA target genes

To investigate the putative functions of potential target genes, the target sequences were annotated using diverse protein databases, including Gene Ontology (GO), Cluster of Orthologous Groups (COG) and Kyoto Encyclopedia of Genes and Genomes (KEGG) [30, 31, 46]. The GO categorization results are expressed as three independent hierarchies for biological processes, cellular components, and molecular functions.

Differential expression of miRNAs between celery varieties

To select differentially expressed miRNAs between the two libraries, the frequency of miRNAs was normalized to one million of the total number of miRNA reads in each sample. The selection method used was according to Audic and Claverie [47]. IDEG6 software [48] was used to analyze the differential expression.

qRT-PCR analysis of miRNA abundance under temperature stress conditions

Total RNA was extracted using an RNAsimple Total RNA Kit (Tiangen) according to the manufacturer’s instructions. Small RNAs were reverse-transcribed into cDNA using the One Step PrimeScript® miRNA cDNA Synthesis Kit (TaKaRa, Dalian, China). Quantitative real-time PCR (qRT-PCR) was performed using the MyiQ single-color real-time PCR detection system (Bio-Rad, Hercules, CA, USA). The reactions were carried out in a total volume of 20 μL containing 2.0 μL of diluted cDNA, 0.8 μL of each primer, and 10 μL of SYBR GreenI Mix with the following cycling profile: 95°C for 30 s; followed by 40 cycles at 95°C for 5 s, 60°C for 20 s; Melting curve analysis was performed (61 cycles at 65°C for 10 s) to verify specific amplification. Each sample was processed in triplicate, and 5.8S rRNA was used as an internal control. The qRT-PCR primers are listed in Additional file 2: Table S2.


  1. 1.

    Obernosterer G, Leuschner PJ, Alenius M, Martinez J: Post-transcriptional regulation of microRNA expression. RNA (New York, NY). 2006, 12 (7): 1161-1167. 10.1261/rna.2322506.

    CAS  Article  Google Scholar 

  2. 2.

    Kren BT, Wong PY, Sarver A, Zhang X, Zeng Y, Steer CJ: MicroRNAs identified in highly purified liver-derived mitochondria may play a role in apoptosis. RNA Biol. 2009, 6 (1): 65-72. 10.4161/rna.6.1.7534.

    CAS  PubMed Central  PubMed  Article  Google Scholar 

  3. 3.

    Lee CT, Risom T, Strauss WM: Evolutionary conservation of microRNA regulatory circuits: an examination of microRNA gene complexity and conserved microRNA-target interactions through metazoan phylogeny. DNA Cell Biol. 2007, 26 (4): 209-218. 10.1089/dna.2006.0545.

    CAS  PubMed  Article  Google Scholar 

  4. 4.

    Voinnet O: Origin, biogenesis, and activity of plant MicroRNAs. Cell. 2009, 136 (4): 669-687. 10.1016/j.cell.2009.01.046.

    CAS  PubMed  Article  Google Scholar 

  5. 5.

    Lee RC, Feinbaum RL, Ambros V: The C. elegans heterochronic gene lin-4 encodes small RNAs with antisense complementarity to lin-14. Cell. 1993, 75 (5): 843-854. 10.1016/0092-8674(93)90529-Y.

    CAS  PubMed  Article  Google Scholar 

  6. 6.

    Reinhart BJ, Slack FJ, Basson M, Pasquinelli AE, Bettinger JC, Rougvie AE, Horvitz HR, Ruvkun G: The 21-nucleotide let-7 RNA regulates developmental timing in Caenorhabditis elegans. Nature. 2000, 403 (6772): 901-906. 10.1038/35002607.

    CAS  PubMed  Article  Google Scholar 

  7. 7.

    Brennecke J, Hipfner DR, Stark A, Russell RB, Cohen SM: bantam encodes a developmentally regulated microRNA that controls cell proliferation and regulates the proapoptotic gene hid in Drosophila. Cell. 2003, 113 (1): 25-36. 10.1016/S0092-8674(03)00231-9.

    CAS  PubMed  Article  Google Scholar 

  8. 8.

    De Lima JC, Loss-Morais G, Margis R: MicroRNAs play critical roles during plant development and in response to abiotic stresses. Genet Mol Biol. 2012, 35 (4 (suppl)): 1069-1077.

    CAS  PubMed Central  PubMed  Article  Google Scholar 

  9. 9.

    Chen L, Ren Y, Zhang Y, Xu J, Zhang Z, Wang Y: Genome-wide profiling of novel and conserved Populus microRNAs involved in pathogen stress response by deep sequencing. Planta. 2012, 235 (5): 873-883. 10.1007/s00425-011-1548-z.

    CAS  PubMed  Article  Google Scholar 

  10. 10.

    Huntzinger E, Izaurralde E: Gene silencing by microRNAs: contributions of translational repression and mRNA decay. Nat Rev Genet. 2011, 12 (2): 99-110. 10.1038/nrg2936.

    CAS  PubMed  Article  Google Scholar 

  11. 11.

    Chen J, Li WX, Xie D, Peng JR, Ding SW: Viral virulence protein suppresses RNA silencing-mediated defense but upregulates the role of microrna in host gene expression. Plant Cell. 2004, 16 (5): 1302-1313. 10.1105/tpc.018986.

    CAS  PubMed Central  PubMed  Article  Google Scholar 

  12. 12.

    Bazzini AA, Hopp HE, Beachy RN, Asurmendi S: Infection and coaccumulation of tobacco mosaic virus proteins alter microRNA levels, correlating with symptom and plant development. Proc Natl Acad Sci U S A. 2007, 104 (29): 12157-12162. 10.1073/pnas.0705114104.

    CAS  PubMed Central  PubMed  Article  Google Scholar 

  13. 13.

    Aukerman MJ, Sakai H: Regulation of flowering time and floral organ identity by a MicroRNA and its APETALA2-like target genes. Plant Cell. 2003, 15 (11): 2730-2741. 10.1105/tpc.016238.

    CAS  PubMed Central  PubMed  Article  Google Scholar 

  14. 14.

    Li X, Bian H, Song D, Ma S, Han N, Wang J, Zhu M: Flowering time control in ornamental gloxinia (Sinningia speciosa) by manipulation of miR159 expression. Annals of botany. 2013, 111 (5): 791-799. 10.1093/aob/mct034.

    CAS  PubMed Central  PubMed  Article  Google Scholar 

  15. 15.

    Mallory AC, Bartel DP, Bartel B: MicroRNA-directed regulation of Arabidopsis AUXIN RESPONSE FACTOR17 is essential for proper development and modulates expression of early auxin response genes. The Plant cell. 2005, 17 (5): 1360-1375. 10.1105/tpc.105.031716.

    CAS  PubMed Central  PubMed  Article  Google Scholar 

  16. 16.

    Sunkar R, Zhu JK: Novel and stress-regulated microRNAs and other small RNAs from Arabidopsis. Plant Cell. 2004, 16 (8): 2001-2019. 10.1105/tpc.104.022830.

    CAS  PubMed Central  PubMed  Article  Google Scholar 

  17. 17.

    Raffo A, Sinesio F, Moneta E, Nardo N, Peparaio M, Paoletti F: Internal quality of fresh and cold stored celery petioles described by sensory profile, chemical and instrumental measurements. Eur Food Res Technol. 2005, 222 (5–6): 590-599.

    Google Scholar 

  18. 18.

    Mongeau R, Siddiqui IR, Emery J, Brassard R: Effect of dietary fiber concentrated from celery, parsnip, and rutabaga on intestinal function, serum cholesterol, and blood glucose response in rats. J Agric Food Chem. 1990, 38 (1): 195-200. 10.1021/jf00091a043.

    CAS  Article  Google Scholar 

  19. 19.

    Song X, Wang Y: Research progress in the medicinal function of celery. J Anhui Agric Sci. 2008, 36 (15): 6360-6361.

    Google Scholar 

  20. 20.

    Eriksson NE, Moller C, Werner S, Magnusson J, Bengtsson U, Zolubas M: Self-reported food hypersensitivity in Sweden, Denmark, Estonia, Lithuania, and Russia. J Investig Allergol Clin Immunol. 2004, 14 (1): 70-79.

    CAS  PubMed  Google Scholar 

  21. 21.

    Hoffmann-Sommergruber K, Demoly P, Crameri R, Breiteneder H, Ebner C, Laimer Da Camara Machado M, Blaser K, Ismail C, Scheiner O, Bousquet J, Menz G: IgE reactivity to Api g 1, a major celery allergen, in a Central European population is based on primary sensitization by Bet v 1. J Allergy Clin Immunol. 1999, 104 (2 Pt 1): 478-484.

    CAS  PubMed  Article  Google Scholar 

  22. 22.

    Bi Y, Liu G, Yang R: MicroRNAs: novel regulators during the immune response. J Cell Physiol. 2009, 218 (3): 467-472. 10.1002/jcp.21639.

    CAS  PubMed  Article  Google Scholar 

  23. 23.

    Rajagopalan R, Vaucheret H, Trejo J, Bartel DP: A diverse and evolutionarily fluid set of microRNAs in Arabidopsis thaliana. Genes Dev. 2006, 20 (24): 3407-3425. 10.1101/gad.1476406.

    CAS  PubMed Central  PubMed  Article  Google Scholar 

  24. 24.

    Song C, Wang C, Zhang C, Korir NK, Yu H, Ma Z, Fang J: Deep sequencing discovery of novel and conserved microRNAs in trifoliate orange (Citrus trifoliata). BMC Genomics. 2010, 11: 431-10.1186/1471-2164-11-431.

    PubMed Central  PubMed  Article  Google Scholar 

  25. 25.

    Chi X, Yang Q, Chen X, Wang J, Pan L, Chen M, Yang Z, He Y, Liang X, Yu S: Identification and characterization of microRNAs from peanut (Arachis hypogaea L.) by high-throughput sequencing. PloS One. 2011, 6 (11): e27530-10.1371/journal.pone.0027530.

    CAS  PubMed Central  PubMed  Article  Google Scholar 

  26. 26.

    Donaire L, Pedrola L, Rosa Rde L, Llave C: High-throughput sequencing of RNA silencing-associated small RNAs in olive (Olea europaea L.). PloS One. 2011, 6 (11): e27916-10.1371/journal.pone.0027916.

    CAS  PubMed Central  PubMed  Article  Google Scholar 

  27. 27.

    Pontes O, Costa-Nunes P, Vithayathil P, Pikaard CS: RNA polymerase V functions in Arabidopsis interphase heterochromatin organization independently of the 24-nt siRNA-directed DNA methylation pathway. Mol Plant. 2009, 2 (4): 700-710. 10.1093/mp/ssp006.

    CAS  PubMed Central  PubMed  Article  Google Scholar 

  28. 28.

    Jia Y, Lisch DR, Ohtsu K, Scanlon MJ, Nettleton D, Schnable PS: Loss of RNA-dependent RNA polymerase 2 (RDR2) function causes widespread and unexpected changes in the expression of transposons, genes, and 24-nt small RNAs. PLoS Genet. 2009, 5 (11): e1000737-10.1371/journal.pgen.1000737.

    PubMed Central  PubMed  Article  Google Scholar 

  29. 29.

    Bonnet E, Wuyts J, Rouze P, Van de Peer Y: Evidence that microRNA precursors, unlike other non-coding RNAs, have lower folding free energies than random sequences. Bioinformatics. 2004, 20 (17): 2911-2917. 10.1093/bioinformatics/bth374.

    CAS  PubMed  Article  Google Scholar 

  30. 30.

    Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, Harris MA, Hill DP, Issel-Tarver L, Kasarskis A, Lewis S, Matese JC, Richardson JE, Ringwald M, Rubin GM, Sherlock G: Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. 2000, 25 (1): 25-29. 10.1038/75556.

    CAS  PubMed Central  PubMed  Article  Google Scholar 

  31. 31.

    Tatusov RL, Fedorova ND, Jackson JD, Jacobs AR, Kiryutin B, Koonin EV, Krylov DM, Mazumder R, Mekhedov SL, Nikolskaya AN, Rao BS, Smirnov S, Sverdlov AV, Vasudevan S, Wolf YI, Yin JJ, Natale DA: The COG database: an updated version includes eukaryotes. BMC Bioinforma. 2003, 4: 41-10.1186/1471-2105-4-41.

    Article  Google Scholar 

  32. 32.

    Fu N, Wang Q, Shen HL: De novo assembly, gene annotation and marker development using illumina paired-End transcriptome sequences in celery (apium graveolens L.). PloS One. 2013, 8 (2): e57686-10.1371/journal.pone.0057686.

    CAS  PubMed Central  PubMed  Article  Google Scholar 

  33. 33.

    Axtell MJ, Snyder JA, Bartel DP: Common functions for diverse small RNAs of land plants. Plant Cell. 2007, 19 (6): 1750-1769. 10.1105/tpc.107.051706.

    CAS  PubMed Central  PubMed  Article  Google Scholar 

  34. 34.

    Khraiwesh B, Zhu JK, Zhu J: Role of miRNAs and siRNAs in biotic and abiotic stress responses of plants. Biochim Biophys Acta. 2012, 1819 (2): 137-148. 10.1016/j.bbagrm.2011.05.001.

    CAS  PubMed Central  PubMed  Article  Google Scholar 

  35. 35.

    Li Y, Li C, Ding G, Jin Y: Evolution of MIR159/319 microRNA genes and their post-transcriptional regulatory link to siRNA pathways. BMC Evol Biol. 2011, 11 (1): 122-10.1186/1471-2148-11-122.

    CAS  PubMed Central  PubMed  Article  Google Scholar 

  36. 36.

    Rodriguez RE, Mecchia MA, Debernardi JM, Schommer C, Weigel D, Palatnik JF: Control of cell proliferation in Arabidopsis thaliana by microRNA miR396. Development. 2010, 137 (1): 103-112. 10.1242/dev.043067.

    CAS  PubMed Central  PubMed  Article  Google Scholar 

  37. 37.

    Zhou X, Wang G, Sutoh K, Zhu JK, Zhang W: Identification of cold-inducible microRNAs in plants by transcriptome analysis. Biochim Biophys Acta. 2008, 1779 (11): 780-788. 10.1016/j.bbagrm.2008.04.005.

    CAS  PubMed  Article  Google Scholar 

  38. 38.

    Xin M, Wang Y, Yao Y, Xie C, Peng H, Ni Z, Sun Q: Diverse set of microRNAs are responsive to powdery mildew infection and heat stress in wheat (Triticum aestivum L.). BMC Plant Biol. 2010, 10: 123-10.1186/1471-2229-10-123.

    PubMed Central  PubMed  Article  Google Scholar 

  39. 39.

    Liu HH, Tian X, Li YJ, Wu CA, Zheng CC: Microarray-based analysis of stress-regulated microRNAs in Arabidopsis thaliana. RNA (New York, NY). 2008, 14 (5): 836-843. 10.1261/rna.895308.

    CAS  Article  Google Scholar 

  40. 40.

    Zeng C, Wang W, Zheng Y, Chen X, Bo W, Song S, Zhang W, Peng M: Conservation and divergence of microRNAs and their functions in Euphorbiaceous plants. Nucleic Acids Res. 2010, 38 (3): 981-995. 10.1093/nar/gkp1035.

    CAS  PubMed Central  PubMed  Article  Google Scholar 

  41. 41.

    Lu S, Sun YH, Chiang VL: Stress-responsive microRNAs in Populus. Plant J. 2008, 55 (1): 131-151. 10.1111/j.1365-313X.2008.03497.x.

    CAS  PubMed  Article  Google Scholar 

  42. 42.

    Zuker M: Mfold web server for nucleic acid folding and hybridization prediction. Nucleic Acids Res. 2003, 31 (13): 3406-3415. 10.1093/nar/gkg595.

    CAS  PubMed Central  PubMed  Article  Google Scholar 

  43. 43.

    Meyers BC, Axtell MJ, Bartel B, Bartel DP, Baulcombe D, Bowman JL, Cao X, Carrington JC, Chen X, Green PJ, Griffiths-Jones S, Jacobsen SE, Mallory AC, Martienssen RA, Poethig RS, Qi Y, Vaucheret H, Voinnet O, Watanabe Y, Weigel D, Zhu JK: Criteria for annotation of plant MicroRNAs. Plant Cell. 2008, 20 (12): 3186-3190. 10.1105/tpc.108.064311.

    CAS  PubMed Central  PubMed  Article  Google Scholar 

  44. 44.

    Allen E, Xie Z, Gustafson AM, Carrington JC: microRNA-directed phasing during trans-acting siRNA biogenesis in plants. Cell. 2005, 121 (2): 207-221. 10.1016/j.cell.2005.04.004.

    CAS  PubMed  Article  Google Scholar 

  45. 45.

    Schwab R, Palatnik JF, Riester M, Schommer C, Schmid M, Weigel D: Specific effects of MicroRNAs on the plant transcriptome. Dev Cell. 2005, 8 (4): 517-527. 10.1016/j.devcel.2005.01.018.

    CAS  PubMed  Article  Google Scholar 

  46. 46.

    Kanehisa M, Goto S: KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000, 28 (1): 27-30. 10.1093/nar/28.1.27.

    CAS  PubMed Central  PubMed  Article  Google Scholar 

  47. 47.

    Audic S, Claverie JM: The significance of digital gene expression profiles. Genome Res. 1997, 7 (10): 986-995.

    CAS  PubMed  Google Scholar 

  48. 48.

    Romualdi C, Bortoluzzi S, d’Alessi F, Danieli GA: IDEG6: a web tool for detection of differentially expressed genes in multiple tag sampling experiments. Physiol Genomics. 2003, 12 (2): 159-162.

    CAS  PubMed  Article  Google Scholar 

Download references


The research was supported by National Natural Science Foundation of China (31272175); New Century Excellent Talents in University (NCET-11-0670); Graduate Educated Innovation Project of Jiangsu (CXZZ13_0297); Jiangsu Natural Science Foundation (BK20130027); Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD) and Jiangsu Shuangchuang Project.

Author information



Corresponding author

Correspondence to Ai-Sheng Xiong.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

ASX conceived and designed the experiments. MYL, FW, ZSX, QJ, GFT and JM performed the experiments. MYL, FW, JQ and ASX analyzed the data. ASX contributed reagents/materials/analysis tools. MYL wrote the paper. MYL and ASX revised the paper. All authors read and approved the final version of manuscript.

Electronic supplementary material

Authors’ original submitted files for images

Rights and permissions

Reprints and Permissions

About this article

Cite this article

Li, M., Wang, F., Xu, Z. et al. High throughput sequencing of two celery varieties small RNAs identifies microRNAs involved in temperature stress response. BMC Genomics 15, 242 (2014).

Download citation


  • microRNAs
  • Temperature stress
  • Deep sequencing
  • qRT-PCR
  • Celery