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

BACKGROUND: 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. RESULTS: 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. CONCLUSION: 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.