Identification and characterization of Prunus persica miRNAs in response to UVB radiation in greenhouse through high-throughput sequencing
BMC Genomics volume 18, Article number: 938 (2017)
MicroRNAs (miRNAs) are small non-coding RNAs that regulate gene expression of target mRNAs involved in plant growth, development, and abiotic stress. As one of the most important model plants, peach (Prunus persica) has high agricultural significance and nutritional values. It is well adapted to be cultivated in greenhouse in which some auxiliary conditions like temperature, humidity, and UVB etc. are needed to ensure the fruit quality. However, little is known about the genomic information of P. persica under UVB supplement. Transcriptome and expression profiling data for this species are therefore important resources to better understand the biological mechanism of seed development, formation and plant adaptation to environmental change. Using a high-throughput miRNA sequencing, followed by qRT-PCR tests and physiological properties determination, we identified the responsive-miRNAs under low-dose UVB treatment and described the expression pattern and putative function of related miRNAs and target genes in chlorophyll and carbohydrate metabolism.
A total of 164 known peach miRNAs belonging to 59 miRNA families and 109 putative novel miRNAs were identified. Some of these miRNAs were highly conserved in at least four other plant species. In total, 1794 and 1983 target genes for known and novel miRNAs were predicted, respectively. The differential expression profiles of miRNAs between the control and UVB-supplement group showed that UVB-responsive miRNAs were mainly involved in carbohydrate metabolism and signal transduction. UVB supplement stimulated peach to synthesize more chlorophyll and sugars, which was verified by qRT-PCR tests of related target genes and metabolites’ content measurement.
The high-throughput sequencing data provided the most comprehensive miRNAs resource available for peach study. Our results identified a series of differentially expressed miRNAs/target genes that were predicted to be low-dose UVB-responsive. The correlation between transcriptional profiles and metabolites contents in UVB supplement groups gave novel clues for the regulatory mechanism of miRNAs in Prunus. Low-dose UVB supplement could increase the chlorophyll and sugar (sorbitol) contents via miRNA-target genes and therefore improve the fruit quality in protected cultivation of peaches.
As an important environment signal, sunlight provides energy for the growth and development of plants , but its ultraviolet (UV) radiation part causes abiotic stress potentially influence the biological processes of plants. Since the late 1980s when awareness of stratospheric ozone layer depletion triggered concerns about the potentially harmful effects of increased UVB radiation, many studies have shown that UVB causes non-specific damage to DNA, proteins and lipids [2,3,4]. On the other hand, there is overwhelming evidence that other than substantially impeding plant growth, low-dose UVB is an environmental regulator affecting gene expression, cellular and metabolic activities, and growth and development [5,6,7,8]. Whether UVB radiation is a stressor or a regulator is determined by the fluence rate and exposure time . Nevertheless, the regulatory mechanism of plants responding to the UVB-lack environment, for example in the greenhouse where the UVB radiation level is 30%–70% lower than outdoors, were rarely reported .
Most of the photomorphogenic responses to low-dose UVB are mediated by the photoreceptor UV RESISTANCE LOCUS8 (UVR8). Subsequent structural and functional characterization revealed that the UVR8 has a unique regulatory mechanism in photoreception [10,11,12]. After UVB treatment, UVR8 interacts with the E3 ubiquitin ligase (transducin/WD40 repeat-like superfamily protein) CONSTITUTIVELY PHOTOMORPHOGENIC1 (COP1), following the ubiquitination of the basic leucine-zipper (bZIP family) transcription factor ELONGATED HYPOCOTYL5 (HY5) which is primarily in the initiation of photo-morphogenesis [13,14,15,16].
MiRNAs, small endogenous non-coding RNAs approximately 21–24 nucleotides (nt) in length, play an important role in regulating gene expression at the post-transcriptional level [17,18,19]. A large number of miRNAs have been recently identified in plants via high-throughput sequencing, and numerous miRNAs have been entered into the miRBase 21. MiRNAs are involved in regulating growth, development, root initiation and development, hormone balance, floral morphogenesis and reproductive performance [20, 21]. Stress-regulated miRNAs in plants confer resistance to the extreme conditions, including UVB, drought, salt, cold and heat. In addition, the expression of miRNAs can alter the behavior of plants in response to both abiotic and biotic stresses [21, 22]. Previous reports have shown that miRNA induction was involved in regulating auxin signaling via miR160, miR167 and miR393 thus, becoming an important strategy for photomorphism in plants .
China is the largest producer of both outside-grown and inside-grown peaches and nectarines in the world. There are nearly 16,000 ha of protected peach and nectarines cultivation, 2.3% of the total area . As a new agricultural form, protected production has been rapidly developed. Peach (Prunus persica), which has been cultivated for more than 4000 years, is one of the most important fruits in the world . Peach has a small genome and it reaches reproductive maturity in a relatively short time. In 2010, the Genome Sequencing Project of the peach double haploid cultivar ‘Lovell’ was completed, which generated, 230 Mb genome sequence and 202 assembly scaffolds . Therefore, peach is considered to be a useful forest model species for genetic and ecological research. Following these findings, several reports on the identification of miRNAs in different peach tissues have been published [27,28,29]. Meanwhile, peach has some unique biological features not commonly found in other model organisms, such as a 3–5 year juvenile period before blossom, the formation of fleshy fruit with a hardened endocarp and chilling-requirement dormancy mechanism [30,31,32,33,34]. Nectarine, because of its low-chill and nutritional value, is selected as one of the most important fruits in the protected cultivation industry which targeted early and high markets. In previous research, we found that the supplement of UVB radiation can improve the fruit quality and the ability to compete for C-assimilate [35, 36]. Considering the distinct environment especially the light condition in greenhouse, more deep studies related to the molecular and metabolic mechanism under UVB irradiance are needed.
In this study, we generated over 2 billion bases of high-quality RNA sequence with Illumina platform. In a single run, we identified 31,763,592 raw sequences including thousands of seed target and metabolism genes. Our results identified a series of differentially expressed miRNAs/target genes that were predicted to be low-dose UVB-responsive. The correlation between transcriptional profiles and metabolites’ contents helped elucidate the regulatory mechanisms of peach under the UVB supplement. Our findings of correlation among miRNA, target genes and metabolites provided clues for breeding with high-quality fruit and other properties which are suitable for greenhouse cultivation.
Analysis of miRNA sequences
Using Illumina sequencing, a total of 74,119,581 and 76,123,247 raw reads were obtained from control and UVB-treatment groups, respectively (Table 1). After discarding 3′ adapter deletions, insertion deletions, 5′ adapter contaminants, poly-A sequences and sequences less than 18 nt from the high-quality reads, 52,375,087 and 56,244,851 clean reads were used for further analysis. The proportions of clean reads were 70.79% and 73.85% of the total reads obtained from the two libraries, respectively.
The small RNA (sRNA) reads were typically 18 to 30 nt in length (Fig. 1). Among these sequences, 21 nt sRNAs were the most abundant in the two libraries, accounting for 16.56% and 13.22% of the total reads, followed by 24 nt sRNAs, which accounted for 11.79% and 10.26% of the total reads, respectively. Furthermore, we observed that the number of less than 24 nt length sequences in the control libraries was more abundant than that in the treatment libraries (74.08% and 64.14%, respectively). Additionally, a large proportion of unique sequences (>85% in both libraries) were unclassified sRNAs, suggesting a broad existence of miRNAs in peach.
Known miRNAs in peach
Known miRNAs in peach affected by UVB were identified through homologous alignment analysis using the plant miRNA in miRBase 21. A total of 164 known miRNAs belonging to 59 families were obtained from the deep sequencing. The dominant miRNA families are shown in Additional file 1: Table S1, and most of these miRNAs were largely conserved in various plant species. The expression levels of a few miRNA families, such as miR166, miR1511 and miR398, were evidentially high in both libraries. Some conserved miRNAs were reported only in few species, such as miR3627 in five (Vitis vinifera, Populus trichocarpa, Malus domestica, Solanum tuberosum, and P. persica), miR1511 in three (M. domestica, S. tuberosum, and P. persica), while miR8133 only in P. persica.
MiRNAs have a broad range of expression levels, varying from several to millions of reads. Most of the conserved miRNAs were identified from the two libraries, and certain miRNAs were abundant in some samples but scarce or even lacking in other samples. For example, the expression of miR159 generated 323,238 and 128,812 reads in the CK and UVB libraries respectively. Moreover, the number of reads for different members of the same family varied widely. For example, the expression of miR7122a in both libraries generated 21,291 reads while only 417 reads of miR7122b in both libraries.
Novel miRNAs in peach
In the present study, 109 novel miRNAs from peach deep sequencing were identified (Additional file 2: Table S2). The most abundant miRNA was Pp04_27840-3p, with 33,630 and 42,292 reads in the UVB libraries and CK libraries, respectively. Many novel miRNAs from our database were conserved with miRNAs from other species to a certain degree (<50%), such as Pp02_15663-5p and Pp03–22,312-3p corresponding to miR172d and miR2950-5p in grape (Vitis vinifera), respectively.
Differential expression of miRNA
The comparison of the miRNA expression levels in CK and UVB groups showed that 164 known and 109 novel miRNAs were identified in the two libraries. The analysis of the differential expression of miRNAs in the UVB treatment and control libraries showed that 9 known and 6 novel differentially expressed miRNAs from the two libraries might play important roles in the UVB response(Fig. 2). In brief, 8 miRNAs were up-regulated, including Pp03–22,312-3p, Pp03–22,312-5p, Pp05–19,842-3p, Pp06–35,148-3p, Pp06–35,148-5p, miR397, miRNA171d-3p and miRNA3627-5p, and 7 miRNAs were down-regulated, including miRNA395d, miRNA395e, miRNA7122b-5p, miRNA399a, miRNA399b, miRNA8133-3p, and Pp05–28,899-3p in UVB treatment.
Target prediction and functional analysis
A large number of targets were predicted for most differentially expressed miRNAs. As for the known miRNAs, there were 1928 pathways that accounted for the largest percentage of the total targets. The results of Gene ontology classification and top 30 preferential KEGG pathway analysis revealed that target genes of these miRNAs were involved in various biological and biochemical processes in plant growth and development (Figs. 3 and 4) (Additional file 3: Table S3, Additional file 4: Table S4), such as porphyrin and chlorophyll metabolism (pper00860), pentose and glucuronate interconversions (pper00040), citrate cycle (TCA cycle, pper04712 and pper00020), circadian rhythm - plant and metabolic pathways (pper01100).
Expression validation of UVB-responsive miRNA
Totally, we verified the expression of 31 miRNAs (Table 2) via qRT-PCR in which 15 miRNAs were retrieved in the database (Fig. 5). We selected miR5059 as reference miRNA  and 15 other miRNAs that related to UVB irradiance. The primers were listed in Additional file 5: Table S5. The qRT-PCR results showed the consistency with sequencing data except that miR397 levels under UVB treatment did not have significant difference with CK. The results showed some difference with the database. Many miRNAs were up-regulated other than down-regulation from high-throughput sequencing, such as miR398a-5p (6.67-fold), miR398a-3p (2.09-fold), miR6263 (5.21-fold), miR6260 (2.31-fold), and miR319a (1.75-fold). On the contrary, miR1511 (2.31-fold), miR171c (2.31-fold), and miR3627-3p (1.61-fold) were down-regulated under UVB treatment (Table 2; Fig. 5). Some conserved UVB-responsive miRNAs were not remarkably expressed such as miR156a, miR160a, miR166a, miR393a, miR402a which were considered to be involved in plant acclimation and adaptation of biotic or abiotic stress.
Chlorophyll metabolism under UVB treatment
KEGG pathway analysis of miRNA target genes (Fig. 4) showed that genes involved in chlorophyll metabolism might be significantly regulated by UVB. We determined chlorophyll content and crucial miRNA/genes expression levels to give clues for UVB effects on chlorophyll metabolism. Spectrophotometry results showed that the content of chlorophyll was gradually increased during the whole development period (Fig. 6a). Its content reached up to 5.24 mg/g Fw in the mature stage of UVB-group, higher than 4.17 mg/g Fw of the control group. The expression levels of miR171c in UVB-group were lower than the untreated samples, which was 2.3-fold difference in the mature stage (Fig. 6b). Fig. 6c showed the different expression profiles of genes involved in chlorophyll metabolism. Scarecrow-like protein (SCL), the target gene of miR171c, was highly up-regulated during the developmental period. On the contrary, pchlide oxidoreductase C (PORC), the downstream gene of SCL was remarkably down-regulated. The others did not show significant difference between CK and UVB groups. Specifically, SCL was up-regulated by 2.2-fold after UVB treatment in fruit mature phase. PORC was down-regulated by 2.8-fold (Fig. 6c). The gene information and sequences of primers were listed in Additional file 6: Table S6.
Carbohydrate metabolism under UVB treatment
Carbohydrate metabolism was another pathway targeted by miRNAs predicted by KEGG analysis. The contents of carbohydrates in leaves were significantly affected by UVB treatment (Fig. 7). Sorbitol, the main form of sugar, was up-regulated by 1.24-fold with UVB stimulation in the mature period (Fig. 7a), while the content of sucrose was dramatically decreased and did not show significant difference with/without UVB treatment (Fig. 7b). Fructose and glucose contents were increased with the developmental process, and the UVB-groups were higher than the control samples (Fig. 7c and d).
We screened 19 genes related with sugar metabolism from P. persica genome, and applied qRT-PCR assay to investigate the correlation between sugar synthesis and miRNA under UVB treatment (Fig. 8). The crucial genes involved in sucrose pathway did not show significant difference between UVB and control groups (Fig. 8a). However, genes in sorbitol metabolism were expressed distinctly after UVB stimulation (Fig. 8b). The expression levels of Sorbitol −3-orbitol −6-phosphate dehydrogenase (S6PDH), NADP dependent sorbitol dehydrogenase (NADP-SDH) and Sorbitol transporter (SOT) were increased by 2.9-fold, 6.1-fold and 2.0-fold respectively under UVB treatment compared to the control groups. However, NAD dependent sorbitol dehydrogenase (NAD-SDH) was rarely expressed under UVB treatment. Moreover, the expression levels of Pyrophosphate--fructose 6-phosphate 1-phosphotransferase subunit alpha (PFP-α) and Phosphofructokinase (PFK) in hexose pathway were1.8-fold and 1.6-fold higher in UVB group than those in control group (Fig. 8c). Hexose carrier protein (HEX6) showed a 3.2-fold decrease after UVB stimulation (Fig. 8c). These results indicated that UVB had a major influence on the interconversion and transportation of monosaccharides. The gene information and sequences of primers were listed in Additional file 6: Table S6.
There are many problems in the protected cultivation of peach trees, such as the basic theory and growth/development pattern, the environmental intelligent control, the standardization system with high quality and efficiency etc. . It is mainly ascribed to the weak light irradiance and quality, and the quality of fruits in the greenhouse is worse than that in the outdoor cultivation . Our investigation on the regulatory mechanism of complementary UVB will provide significant effects for elucidating the growth and development pattern of P. persica, thus increasing the fruit quality.
High-dose UVB may produce phenols, ROS and cause damage to DNA, proteins, and membranes [40,41,42]. However, low-dose UVB may trigger early adaptation to environment and regulate plants seeding, development and growth such as the induction of alterations in antioxidant status, the regulation of phenylpropanoids, cinnamates, or flavonoids pathways, chlorophyll and pyridoxine biosynthesis pathways . In our previous study, Chen et al. selected three UVB radiation levels and 1.44 Kj m−2 d−1 showed the most effective function on the improvement of total soluble solid, anthocyanin and the repression of total acid in peach cultivated in greenhouse . Therefore the UVB radiation level (1.44 Kj m−2 d−1) was applied in this study as the most suitable dose to regulate the growth and fruit quality of P. persica in the greenhouse environment.
MiRNAs play important roles in plant response and adaptation to environmental change. Thus, understanding the miRNA-mediated regulatory network of UVB supplement will lay the foundation for unraveling the complex molecular genetic mechanism of positive effects on fruit’s agronomic traits improvement. A growing evidences suggested that miRNA-guided gene regulation could play a vital role in plant response to UVB radiation [42, 45, 46]. In the present study, a total of 4.02 M and 3.83 M unique sRNA sequences were obtained from the control and UVB-treatment libraries, respectively, suggesting adequate sequencing depth for further analysis. The majority of total sRNA reads ranged from 18 nt to 30 nt in length (Fig. 1), which was consistent with the typical size for Dicer-derived products . The most abundant length is 21 nt followed by 24 nt which was consistent with previous studies in peach [48,49,50,51].
In our study, the expression levels of 2 known miRNAs were highly up-regulated (miR171d-3p, miR3627-5p) and 6 known miRNAs showed significant down-regulation (miR395d, miR395e, miR399a, miR399b, miR7122b-5p, miR8133-3p). These miRNAs were predicted to be involved in distinct metabolic pathways.
MiR171c was predicted to target SCL6, SCL22 and SCL27, a family of transcription factors which were involved in the morphogenesis, proliferation of meristematic cells, polar organization and chlorophyll synthesis [52,53,54]. Further study found that its target gene tomato (Solamum lycopersicum) gras transcription factor gene (SlGRAS24) impacts multiple agronomical traits, such as plant height, flowering time, leaf architecture, lateral branch number, root length, fruit set and development, via regulating gibberellin and auxin homeostasis [55, 56]. It is indicated that the UVB-responsive miRNA171 plays an important role in regulating the growth and development of plant. In the present study, we investigated the potential relation between miR171c and chlorophyll synthesis. The content of chlorophyll was gradually increased during the whole development period under UVB treatment, which was in accordance with the studies by Chen . Although previous reports showed that miRNA171 was up-regulated under light , our results found that miRNA171c was less expressed after UVB stimulation. Further analysis on genes involved in miR171-SCL model showed that SCL and PORC were up-regulated and down-regulated respectively, which was in accordance with this model. PORC had positive effects as an upstream gene of chlorophyll synthetic pathway. However, our results showed an opposite regulation pattern between PORC and chlorophyll content. Interestingly, the expression levels of Chlorophyll a-b binding protein (CBR) and Red chlorophyll catabolite reductase (PCCR) (genes involved in chlorophyll degradation) had no change after UVB treatment, while the CHLG was up-regulated by 1.6-fold. According to these results, we presumed that the increase of chlorophyll content under UVB might be not only related with miR171-SCL model but also regulated by non miR171-SCL pathway like the regulation of CHLG. Further verification is needed to illustrate the regulatory mechanism of chlorophyll metabolism under UVB.
A conserved miRNA family (miR3627) was reported only in five species (Solanum tuberosum, Malus domestica, Populus trichocarpa, Citrus trifolia, and Prunus persica), and was identified as a chilling responsive miRNA in P. persica . In our study, miRNA3627-5p was up-regulated under UVB treatment, and we found an interesting physiological phenomenon that the germination rate of treatment group (89.2%) was 12% higher than the CK group (77.5%), which occurred 8 months after the termination of UVB treatment. It implied that UVB radiation may have a long-term and sustainable influence on plants via modification called “UVB memory” , and this phenomenon involved a series of miRNAs such as miR3627-5p. Bioinformatics analysis of miR3627-5p in our libraries showed that it had 415 target genes which referred to many pathways such as metabolic, protein and amino acid metabolism, RNA transport and circadian rhythm. However, the molecular mechanism of miR3627-5p regulation in chilling still needs more experimental studies.
The target genes of miR3627-5p: ppa025787mg, ppa007623mg, ppa016917mg, ppa021860mg were involved in pentose and glucuronate interconversions, starch and sucrose metabolism. Besides, miR3627-5p target gene ppa007934mg is involved in carbon metabolism and TCA cycle. The expression levels of these five target genes were all down-regulated (Additional file 7: Figure S1), which were consistent with the up-regulation of miR3627-5p. Also the predicted target genes of novel miRNA Pp03_22,312-5p (upregulated) were involved in sugar metabolism. This indicated that miRNA may participate in the formation of fruit sweetness, an essential characteristic of fruit quality which was in accordance with previous reports .
Different sugar components in leaves were analyzed by High Performance Liquid Chromatography (HPLC) and we found that the carbohydrate contents changed variously under UVB treatment (Fig. 7). Accordingly, we speculated some UVB-responsive miRNAs which targeted sugar-metabolizing pathway, such as miR3627-5p. Further qRT-PCR showed relative expression of the key genes in different sugar metabolic pathways (Fig. 8). Sorbitol is the main form of sugar in Rosaccae leaves and is transferred to fruits for the storage of carbohydrates in the mature stage. The content of sorbitol in our study was gradually increased with the developmental process, which was assistant with previous reports . qPCR results showed that sorbitol-related genes were distinctly expressed with UVB treatment. The relative expression levels of S6PDH was remarkably up-regulated in UVB-group, indicating the increase of S6PDH activities and therefore stimulating the synthesis of sorbitol. NAD-SDH1 and NAD-SDH2 were down-regulated while the levels of NADP-SDH were increased after UVB treatment. This suggested that the plants restricted NAD-SDHs expression to reduce sorbitol degradation while increased NADP-SDH transcripts for the synthesis of fructose from sorbitol under UVB treatment. SOT was up-regulated, responsible for more active transportation of sorbitol, which accorded with the increase of sorbitol contents.
The content of sucrose did not change obviously during the continuous growth period, which could be explained by the slight change of expression levels of genes involved in sucrose metabolism, such as Sucrose synthase (SUS), Sucrose transporter (SUT), Tonoplast monosaccharide transporter (TMT) etc. In hexose pathway, the crucial genes PFP-a and PFK were up-regulated while Fructose-1,6-bisphosphatase (FBPase) and HEX6 were down-regulated, implying the UVB treatment had an influence on hexose metabolism.
In our study, the UVB treatment led to different expression patterns of genes related with sugar metabolism. Previous research have shown that the regulation was mainly mediated by miRNA. Our results were essentially in agreement with these reports, and more investigations on functional verification need to be performed.
In this study, we constructed 2 miRNA libraries for low-dose UVB radiation groups and control groups of P. persica in greenhouse. A total of 164 known and 109 novel miRNAs were identified. In brief, 8 miRNAs were highly up-regulated and 7 miRNAs were significantly down-regulated in the UVB treatment groups, which were mainly predicted to be targeted in signal transduction, carbohydrate metabolism and stress response etc. Combined with qRT-PCR tests and the measurement of the related metabolites, our results showed that low-dose UVB radiation could regulate the expression patterns of some miRNAs e.g. miR3627-5p and Pp_22,312, and cause the expression levels of genes in carbohydrate (sorbitol, fructose, and glucose etc.) and chlorophyll pathway, and therefore indirectly affect sugar contents and fruit quality (Fig. 9). Our study provided a comprehensive database of miRNA for P. persica and a theoretical basis for further investigations of the function of miRNA in regulating the biological features of peach in greenhouse.
Plant material and tissue collection
The experiment was carried out on 7-year-old peach trees (Prunus persica var. nectarine Zhongyou No.5) cultivated in experimental station of Shandong Agriculture university, Taian, China. The trees were treated with 1.44 Kj m−2 d−1 UVB radiation during the whole growing period. UVB was provided by the dedicated UV lamp of 40 W, 297 nm (Nanjing Kazhi), hanging at the position of 1.5 m above the plants. UVB-type single-channel UV irradiator (Beijing Normal University Photoelectric Instrument Factory) equipped with 297 nm probe was used to determine the UVB radiation dose of 1.44 Kj m−2 d−1. The function leaves within the range of 80–120 cm below the lamp were selected. The on/off time of UV light was controlled through the electronic automatic control device, from 7 days after blossom to the fruit mature stage. UV light was kept on from 9 h30 to 10 h30 every day and stopped at cloudy, rainy and snowy days. Function leaves were sampled on the 7th day after blossom and every week after, until the mature period. These samples were washed with DEPC-treated H2O, immediately frozen in liquid nitrogen and stored at −80 °C until use. The function leaves of fruit mature stage were selected for high-quality deep-sequencing.
Small RNA isolation and Illumina sequencing
Total RNA in peach was isolated using the mirVan miRNA Isolation kit (Ambion; Thermo Fisher Scientific, USA) and purified using the miRNeasy Mini kit (Qiagen, Germany), following the manufacturer’s instructions. RNA was quantified using a spectrophotometer (NanoDrop, Thermo Fisher Scientific, USA). Purified RNA was frozen in liquid nitrogen and then stored at −80 °C until required.
We constructed an RNA library using the NEBNext Ultra RNA Library Prep kit for Illumina (New England Biolabs, USA). There were four steps: Firstly, total RNA (approx. 1 μg) was spliced into shorter fragments (200–500 bp) in the NEBNext First Strand Synthesis Reaction Buffer and the fragments were used to produce the double-stranded cDNA; secondly, the cDNA was end-repaired and ligated with Illumina-specific adaptors; thirdly, we used 200 bp inserts from the library and selected suitable fragments for PCR amplification; last, we performed PCR using Phusion High-Fidelity DNA polymerase (New England Biolabs, USA) and purified the products with a QIAquick Nucleotide Removal kit (Qiagen, Germany). Then we sequenced the new RNA library on an Illumina HiSeq 2500 system using 2 × 150 base pairs paired-end sequencing.
Analysis of sequencing data
Clean reads were obtained from raw reads after removing low-quality and adapter reads. SOAP software was used for the mapping of clean reads to the peach genome. The non-coding RNAs, including rRNAs, scRNAs, snoRNAs, snRNAs, and tRNAs deposited in the NCBI GenBank database and Rfam (11.0) database, were removed. We also excluded the small RNAs corresponding to the exons and introns of mRNA and repeat sequences. The remaining sRNA sequences were aligned to the miRBase 21 database, with a maximum of two mismatches, to identify known miRNAs in P. persica. The obtained sequences were used to predict hairpin structures using the perpl program. The remaining unannotated sRNAs were used to predict novel miRNAs using Mireap software .
Target prediction of miRNAs and functional analysis
Target prediction of miRNAs followed rules referring to Allen et al. : a. ≤two adjacent mismatches in the miRNA/target duplex; b. ≤ four mismatches between the sRNA and target gene; c. ≤ 2.5 mismatches at positions 1–12 of the 5′-miRNA/target duplex; d. no mismatches at 10–11 of the miRNA/target duplex; e. no adjacent mismatches at 2–12 of the miRNA/target duplex; f. the minimum free energy (MFE) of the miRNA/target duplex should be 75% of the MFE bound to the perfect complement. Target genes were searched using peach genome information. To better understand the roles of miRNAs in peach under low-dose UVB treatment, the potential target functions were annotated using the Gene Ontology and KEGG pathway database.
Differential expression analysis of miRNAs
The miRNA reads were used to analyze differential expression and determine significant differences between the control and treatment libraries. The frequency of miRNAs was normalized to one million to reduce potential errors before calculating the fold-change, P-value and ratio. Normalized expression = (actual miRNA counts/total counts of clean reads) × 1,000,000. Fold-change = log2 (miRNA normalized read counts in the treatment library/ in the control library).
A fold change larger than 1 or less than −1 and a P-value less than 0.01 suggested the highly significant difference in the miRNA expression between two libraries. When the fold-change was greater than 1 or less than −1 and the P-value was between 0.01 and 0.05, the expression of the miRNA was significantly different between the two libraries. The ratio of miRNA normalized read counts in treatment library/in control library was used to determine changes in the expression of an miRNA in the treatment samples compared with the control samples. When the ratio was more than 2, the miRNA was indicated as up-regulated, and when the ratio was less than 1/2, the miRNA was down-regulated .
Relative expression of miRNA and target genes
qRT-PCR was performed to verify the expression levels of identified miRNAs and targets genes using the IQ5 Quantitative Real-time PCR Detection System (Bio-Rad, California, USA) with the SYBR® PrimeScript™ miRNA RT-PCR Kit (TaKaRa, Dalian, China). The reactions were performed in a total volume of 25 μL containing 2.0 μL of diluted cDNA (100 ng/μL), 1 μL of each primer (10 μM), and 12.5 μL of SYBR Green premix Ex Taq II with the following reaction conditions: 95 °C for 30 s, followed by 40 cycles of 95 °C for 5 s and 60 °C for 20 s, then dissociation curve with 95 °C for 60 s, 55 °C for 30 s and 95 °C for 30 s. The reference genes for qRT-PCR of miRNA and target genes were miR5059 and beta-actin, respectively [65, 66]. Each sample was processed in triplicate. All validated primer sequences of miRNAs/target genes are listed in Table S5 and Table S6.
Metabolites analysis by HPLC
Four sugar components: sucrose, sorbitol, glucose and fructose of leaves, were analyzed according to the method of Karkacier et al. . In a mortar pre-cooled in the −20 °C refrigerator, 100 mg of fresh leaves were ground with liquid nitrogen and extracted with 1 mL NANO pure water into a 2 ml tube. The tubes were then vigorously shaken for 15 s, sonicated for 15 min and centrifuged at 12,000 rpm for 15 min. The supernatant was sterilized by filtering through a 0.45 μm membrane filter and stored at −20 °C prior to sugar components’ measurement using HPLC. The HPLC system was programmed to inject 50 μL crude extracts automatically. Online detection was performed using a Waters 410 differential refractrometer detector and the data were analyzed by Oirigin75 software. The whole program used a MetaCarb 87 °C equipped with a guard column as the analytical column, and the deionized water as the mobile phase with a 0.5 mL min−1 flow rate. Glucose, sorbitol, fructose and sucrose purchased from company were used as standards .
Function Leaves were taken every 7 days after flowering to measure chlorophyll content via spectrophotometry. The maximum UV absorption wavelength of chlorophyll a and chlorophyll b is 645 nm and 663 nm respectively. The total chlorophyll content was analyzed with the following formula :
CT = Ca + Cb = 20.29 A645 + 8.05 A663.
Chlorophyll a-b binding protein
Constitutively photomorphogenic 1
Hexose carrier protein
High Performance Liquid Chromatography
Elongated hypocotyl 5
NAD-sorbitol dehydrogenase 2
Neutral invertase 8
Pyrophosphate--fructose 6-phosphate 1-phosphotransferase subunit alpha
Pyrophosphate--fructose 6-phosphate 1-phosphotransferase subunit beta
Pchlide oxidoreductase C
Red chlorophyll catabolite reductase
Sorbitol −6- phosphate dehydrogenase
A tomato (solamum lycopersicum) GRAS transcription factor gene
Sucrose phosphate synthase
Sucrose synthase 4
Sucrose transporter 2
Sucrose transporter 4
Tonoplast monosaccharide transporter
Ultraviolet radiation B
UV resistance locus 8
Kami C, Lorrain S, Hornitschek P, Fankhauser C. Light-regulated plant growth and development. Curr Top Dev Biol. 2010;91:29–66.
Lidon FJC, Teixeira M, Ramalho JC. Decay of the chloroplast pool of ascorbate switches on the oxidative burst in UV-B-irradiated rice. J Agron Crop Sci. 2012;198:130–44.
Pitzschke A, Forzani C, Hirt H. Reactive oxygen species signaling in plants. Antioxid Redox Signal. 2006;8:1757–64.
Gill SS, Tuteja N. Reactive oxygen species and antioxidant machinery in abiotic stress tolerance in crop plants. Plant Physiol Biochem. 2010;48:909–30.
Jenkins GI. Signal transduction in responses to UVB radiation. Annu Rev Plant Biol. 2009;60:407–31.
Brosche M, Strid A. Molecular events following perception of ultraviolet-B radiation by plants: UVB induced signal transduction pathways and changes in gene expression. Physiol Plant. 2003;117:1–10.
Ballare CL, Caldwell MM, Flint SD, Robinson SA, Bornman JF. Effects of solar ultraviolet radiation on terrestrial ecosystems. patterns, mechanisms, and interactions with climate change. Photochem Photobiol Sci. 2011;10:226–41.
Li FR, Peng SL, Chen BM, Hou YP. A meta-analysis of the responses of woody and herbaceous plants to elevated ultraviolet-B radiation. Acta Oecol. 2010;36:1–9.
Deckmyn G, Gaeyenberghs E, Ceulemans R. Reduced UV-B in greenhouses decreases white clover response to enhance CO2. Environ Exp Bot. 2001;46(2):109–17.
Rizzini L, Favory JJ, Cloix C, Faggionato D, O’Hara A, Kaiserli E, Baumeister R, Schäfer E, Nagy F, Jenkins GI, et al. Perception of UV-B by the Arabidopsis UVR8 protein. Science. 2011;332(6025):103–6.
Christie JM, Arvai AS, Baxter KJ, Heilmann M, Pratt AJ. Plant UVR8 photoreceptor senses UVB by tryptophan-mediated disruption of cross-dimer salt bridges. Science. 2012;335:1492–6.
Wu D, Hu Q, Yan Z, Chen W, Yan C. Structural basis of ultraviolet-B perception by UVR8. Nature. 2012;484:214–9.
Heijde M, Ulm R. UVB photoreceptor-mediated signaling in plants. Trends Plant Sci. 2012;17:230–7.
Yi C, Deng XW. COP1-from plant photomorphogenesis to mammalian tumorigenesis. Trends Cell Biol. 2005;15:618–25.
Saijo Y, Sullivan JA, Wang H, Yang J, Shen Y. The COP1-SPA1 interaction defines a critical step in phytochrome A-mediated regulation of HY5 activity. Genes Dev. 2003;17:2642–7.
Osterlund MT, Hardtke CS, Wei N, Deng XW. Targeted destabilization of HY5 during light-regulated development of Arabidopsis. Nature. 2000;405:462–6.
Jones-Rhoades MW, Bartel DP, Bartel B. MicroRNAs and their regulatory roles in plants. Annu Rev Plant Biol. 2006;57:19–53.
Ramachandran V, Chen X. Small RNA metabolism in Arabidopsis. Trends Plant Sci. 2008;13:368–74.
Voinnet O. Origin, biogenesis, and activity of plant microRNAs. Cell. 2009;136:669–87.
Chen X. MicroRNA biogenesis through Dicer-like 1 protein functions. Proc Natl Acad Sci U S A. 2004;101:12753–8.
Sunkar R, Li Y, Jagadeeswaran G. Functions of microRNAs in plant stress responses. Trends Plant Sci. 2012;17:196–203.
Sunkar R, Chinnusamy V, Zhu J, Zhu JK. Small RNAs as big players in plant abotic stress responses and nutrient deprivation. Trends Plant Sci. 2007;12:301–9.
Pérez-Quintero AL, Quintero A, Urrego O, Vanegas P, López C. Bioinformatic identification of cassava miRNAs differentially expressed in response to infection by Xanthomonas axonopodis pv. manihotis. BMC Plant Biol. 2012;12:1–11.
Wang Z, Niu L. Peach industry status and recommendations. Fruit Growers’ Friend (in Chinese). 2012;11:37–8.
Shulaev V, Korban SS, Sosinski B, Abbott AG, Aldwinckle HS, Folta KM, Lezzoni A, Main D, Arus P, Dandekar AM, et al. Multiple models for Rosaceae genomic. Plant Physiol. 2008;147(3):985–1003.
Verde I, Abbott AG, Scalabrin S, Jung S, Shu S, Marroni F, Zhebentyayeva T, Dettori MT, Grimwood J, Cattonaro F, et al. The high-quality draft genome of peach (Prunus persica) identifies unique patterns of genetic diversity, domestication and genome evolution. Nat Genet. 2013;45(5):487–94.
Zhu H, Xia R, Zhao BY, An YQ, Dardick DC, Callahan MA, Liu ZR. Unique expression, processing regulation, and regulatory network of peach (prunus persica) miRNAs. BMC Plant Biol. 2012;12(1):149.
Gao Z, Luo X, Shi T, Cai B, Zhang Z, Cheng ZM. Identification and validation of potential conserved microRNAs and their targets in peach (Prunus persica). Mol Cells. 2012;34:239–49.
Luo X, Gao Z, Shi T, Cheng Z, Zhang Z, Ni ZJ. Identification of miRNAs and their target genes in peach (Prunus persica L.) using high-throughput sequencing and degradome analysis. PLoS One. 2013;8(11):e79090.
Reig G, Alegre S, Gatius F, Iglesias I. Adaptability of peach cultivars [Prunus persica (L.) Batsch] to the climatic conditions of the Ebro Valley, with special focus on fruit quality. Sci Hortic-Amsterdam. 2015;190:149–60.
Byrne DH. Peach breeding trends: a worldwide perspective. Acta Hort. 2002;592:49–59.
Horvath DP, Anderson JV, Chao WS, Foley ME. Knowing when to grow: signals regulating bud dormancy. Trends Plant Sci. 2003;8(11):534–40.
Cirilli M, Bassi D, Ciacciulli A. Sugars in peach fruit: a breeding perspective. Hortic Res. 2016;2:15067.
Crisosto C, Costa G. The peach: botany, production, and uses. In: Layne DR, Bassi D, editors. . Cambridge: CAB International; 2008. p. 536–49.
Yu N, Li D, Tan Q, Zhang H, Gao D. Effect of UVB radiation on assimilate translocation and distribution in fruiting shoot of protected peach. Chin J Appl Environment Biol. 2013;19(1):157–63.
Yu N, Tan Q, Tan Y, Zhang H, Gao D. Effects of UVB radiation on 15N urea absorption, utilization and distribution in fruiting shoot of peach under protected culture. Plant Nutrition Fertilizer Sci. 2012;18(2):491–8.
Luo X, Shi T, Sun H. Selection of suitable inner reference genes for normalization of microRNAs expression response to abiotic stresses by RT-qPCR in leaves, flowers and young stems of peach. Sci Hortic-Amsterdam. 2014;165(3):281–7.
Gao D. The current conditions and developing tendency of protected cultivation of fruit trees in China. Deciduous Fruits (in Chinese). 2016;48(1):1–4.
Li ZY, Gao DS, Qian S, Zhang JH, Li ZJ, Wang C. Effects of different light environments on the fruit quality of peach in greenhouse. J Anhui Agri Sci. 2009;37(21):9933V9934–63.
Brown BA, Jenkins GI. UV-B signaling pathways with different fluence-rate response profiles are distinguished in mature Arabidopsis leaf tissue by requirement for UVR8, HY5, and HYH. Plant Physiol. 2008;146:576–88.
Parul P, Samiksha S, Rachana S. Changing scenario in plant UV-B research: UV-B from a generic stressor to a specific regulator. J Photoch Photobio B. 2015;153:334–43.
Ulm R, Nagy F. Signalling and gene regulation in response to ultraviolet light. Curr Opin Plant Biol. 2005;8:477–82.
Hideg E, Jansen MA, Strid A. UV-B exposure, ROS, and stress: inseparable companions or loosely linked associates? Trends Plant Sci. 2013;18:107–15.
Chen XD. The effects of ultraviolet-B radiation intensity and different plasic film on development characteristics of peach flower and fruit in protected culture. Shandong Agriculture University, Horticulture and Engineering college: Master Thesis; 2009.
Hectors K, Prinsen E, De CW, Jansen MAK, Guisez Y. Arabidopsis thaliana plants acclimated to low dose rates of ultraviolet B radiation show specific changes in morphology and gene expression in the absence of stress symptoms. New Phytol. 2007;175:255–70.
Wang B, Sun Y, Song N. Identification of UVB-induced microRNAs in wheat. Genet Mol Res. 2013;12(4):4213–21.
Henderson IR, Zhang X, Lu C, Johnson L, Meyers BC, Green PJ. Dissecting Arabidopsis thaliana DICER function in small RNA processing, gene silencing and DNA methylation patterning. Nat Genet. 2006;38:721–5.
Chen L, Ren Y, Zhang Y, Xu J, Sun F, Zhang Z, Wang Y. Genome-wide identification and expression analysis of heat-responsive and novel microRNAs in Populus tomentosa. Gene. 2012;504:160–5.
Chen L, Zhang Y, Ren Y, Xu J, Zhang Z, Wang Y. Genome-wide identification of cold-responsive and new microRNAs in Populus tomentosa by high-throughput sequencing. Biochem Biophys Res Commun. 2012;417:892–6.
Li B, Duan H, Li J, Deng X, Yin W. Global identification of miRNAs and targets in Populus euphratica under salt stress. Plant Mol Biol. 2013;81:525–39.
Ren Y, Chen L, Zhang Y, Kang X, Zhang Z, Wang Y. Identification of novel and conserved Populus tomentosa microRNA as components of a response to water stress. Funct Integr Genomics. 2012;12:327–39.
Wang L, Mai YX, Zhang YC, Luo Q, Yang HQ. MicroRNA171c-targeted SCL6-II, SCL6-III, and SCL6-IV genes regulate shoot branching in Arabidopsis. Mol Plant. 2010;3:794–806.
Curaba J, Talbot M, Li Z, Helliwell C. Over-expression of microRNA171 affects phase transitions and floral meristem determinancy in barley. BMC Plant Biol. 2013;13:6.
Fan T, Li X, Yang W, Xia K, Ouyang J, Zhang M. Rice osa-mir171c mediates phase change from vegetative to reproductive development and shoot apical meristem maintenance by repressin four OsHAM transcription factors. PLoS One. 2015;10(5):e0125833.
Huang W, Peng S, Xian Z, Lin D, Hu G, Yang L, Ren M, Li Z. Overexpression of a tomato mir171 target gene SlGRAS24 impacts multiple agronomical traits via regulating gibberellin and auxin homeostasis. Plant Biotechnol J. 2017;15(4):472–88.
Huang W, Xian Z, Kang X, Tang N, Li Z. Genome-wide identification, phylogeny and expression analysis of GRAS gene family in tomato. BMC Plant Biol. 2015;15:209.
Ma Z, Hu X, Cai W, Huang W, Xhou X, Luo Q, Yang H, Wang J, Huang J. Arabidopsis miR171-targeted scarecrow-like proteins bind to GT cis-elements and mediate gibberellin-regulated chlorophyll biosynthesis under light conditions. PLoS Genet. 2014;10(8):E1004519.
Barakat A, Sriram A, Park J, Zhebentyayeva T, Main D, Abbott A. Genome wide identification of chilling responsive microRNAs in Prunus persica. BMC Genomics. 2012;13:481.
Müller-Xing R, Xing Q, Goodrich J. Footprints of the sun: memory of UV and light stress in plants. Front Plant Sci. 2014;5:474.
Kroger M, Meister K, Kava R. Low-calorie sweeteners and other sugar substitutes: a review of the safety issues. Compr Rev Food Sci F. 2006;5(2):35–47.
Guo X, Li S, Liu G, Fu Z, Li S. Seasonal changes in carbohydrate content and reated enzyme activity in fruit and leaves of “Yanfengyihao” peach variety. J Fruit Science (in Chinese). 2004;21(3):196–200.
Ding D, Li W, Han M, Wang Y, Fu Z, Wang B, Tang J. Identification and characterization of maize microRNAs involved in developing ears. Plant Biol. 2013;16(1):9–15.
Allen E, Xie Z, Gustafson AM, Sung GH. Evolution of microRNA genes by inverted duplication of target gene sequences in Arabidopsis thaliana. Nat Genet. 2004;36(12):1282.
Robinson MD, McCarthy DJ, Smyth GK. EdgeR: a bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26(1):139–40.
Luo X, Shi T, Sun H, Song J, Ni Z, Gao Z. Selection of suitable inner reference genes for normalisation of microRNA expression response to abiotic stresses by RT-qPCR in leaves, flowers and young stems of peach. Sci Hortic-Amsterdam. 2014;165(3):281–7.
Wang D, Gao Z, Du P, Xiao W, Tan Q, Chen X, Li L, Gao D. Expression of ABA metabolism-related genes suggests similarities and differences between seed dormancy and bud dormancy of peach (Prunus persica). Front Plant Sci. 2016;6:1248.
Karkacier M, Erbas M, Uslu MK, Aksu M. Comparison of different extraction and detection methods for sugars using amino-bonded phase HPLC. J Chromatogr Sci. 2003;41(6):331–43.
Sornkanok V, Zheng H, Peng Q, Jiang Q, Wang H, Fang T, Liao L, Wang L, He H, Han Y. Assessment of sugar components and genes involved in the regulation of sucrose accumulation in peach fruit. J Agric Food Chem. 2016;64:6723–9.
Shabala SN, Shabala SI, Martynenko AI, Babourina O, Newman IA. Salinity effect on bioelectric activity, growth, Na+ accumulation and chlorophyll fluorescence of maize leaves: a comparative survey and prospects for screening. Aust J Plant Physiol. 1998;25:609–16.
We thank WX in the lab for providing the trees and BBW and ZJZ for DNA extraction. We thank GW from the CapitalBio Technology for sequencing performed and bioinformatics support. We are grateful to all co-authors who participated in the studies mentioned in the text that were published by our groups. All authors declared no conflict of interest.
This work was supported by the National Natural Science Foundation of China [Grant number 31601706], Natural Science Foundation of Shandong Province [Grant number ZR2016CM09] and Science and Technology Innovation Team of Shandong Agriculture University-Facility Horticulture Advantages Team (SYL2017YSTD07). The funding institutions had no direct role in study design, sample collection, analysis, and interpretation of date, nor in manuscript writing. Annual reports were submitted to the funding institutions tracking the progress of the projects.
Availability of data and materials
The dataset supporting the conclusions of this article is available in the NCBI database, with the BioProject: ID PRJNA391378 and Locus tag prefix CEW84.
Ethics approval and consent to participate
The peach trees and fruits were cultivated and collected in the No. 15 greenhouse of Scientific and Technological Innovation Park in Shandong Agriculture University. The authors declared the compliance with institutional, national and international guidelines. The experiment was environmentally friendly and had no negative effects on the local ecosystem.
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.
The summary of known miRNAs prediction and expression in control and UVB supplement libraries. (PDF 250 kb)
The summary of novel miRNAs prediction and expression in control and UVB supplement libraries. (PDF 162 kb)
Details of targets genes and their annotation, GO classification, and KEGG pathway for the known miRNAs. (PDF 2499 kb)
Details of targets genes and their annotation, GO classification, and KEGG pathway for the known miRNAs. (PDF 2823 kb)
Primers for qRT-PCR verification of miRNAs. (PDF 17 kb)
The names, Genbank IDs, involved pathways of verified target genes and their primers for qRT-PCR tests. (PDF 79 kb)
qRT-PCR analysis of five target genes predicted for miR3627-5p. Beta-actin was the internal control. Each experiment was performed with three biological replicates. (PDF 66 kb)
About this article
Cite this article
Li, S., Shao, Z., Fu, X. et al. Identification and characterization of Prunus persica miRNAs in response to UVB radiation in greenhouse through high-throughput sequencing. BMC Genomics 18, 938 (2017). https://doi.org/10.1186/s12864-017-4347-5