|Home | About | Journals | Submit | Contact Us | Français|
The yellow dwarf disease associated with phytoplasmas is one of the most devastating diseases of mulberry and the pathogenesis involved in the disease is poorly understood. To analyze the molecular mechanisms mediating gene expression in mulberry-phytoplasma interaction, the comprehensive sRNA changes of mulberry leaf in response to phytoplasma-infection were examined. A total of 164 conserved miRNAs and 23 novel miRNAs were identified, and 62 conserved miRNAs and 13 novel miRNAs were found to be involved in the response to phytoplasma-infection. Meanwhile, target genes of the responsive miRNAs were identified by sequencing of the degradome library. In addition, the endogenous siRNAs were sequenced, and their expression profiles were characterized. Interestingly, we found that phytoplasma infection induced the accumulation of mul-miR393-5p which was resulted from the increased transcription of MulMIR393A, and mul-miR393-5p most likely initiate the biogenesis of siRNAs from TIR1 transcript. Based on the results, we can conclude that phytoplasma-responsive sRNAs modulate multiple hormone pathways and play crucial roles in the regulation of development and metabolism. These responsive sRNAs may work cooperatively in the response to phytoplasma-infection and be responsible for some symptoms in the infected plants.
Mulberry tree is a perennial woody plant that has long been planted, and is of considerable economic importance to sericulture1. It is susceptible to a number of diseases2, among which yellow dwarf disease associated with aster yellows phytoplasma is one of the most devastating diseases3. Phytoplasmas are obligate endocellular parasites lacking cell walls, and have devastating effects on more than several hundred plant species worldwide4. Infected plants show dramatic changes in morphology including yellowing, phyllody, stunting, proliferation, and witches' broom5. These abnormal developments in plants can cause devastating losses in agriculture and forestry6. Because phytoplasmas were unculturable in pure cultures in laboratory broth, these organisms remain one of the most poorly characterized plant pathogens, and the mechanisms underlying their pathogenicities are not yet understood7.
When faced with pathogen invasion, plants activate sophisticated response mechanisms to reprogram gene expression8. MicroRNAs (miRNAs) and endogenous small interfering RNAs (siRNAs) are two main categories of small regulatory RNAs which globally regulate plant immunity by inhibiting target gene expression at transcriptional or post-transcriptional level9. Increasing evidence indicated that miRNAs serve as one important mechanism for mediating gene expression in plant-microorganism interactions10,11,12,13,14. As for the phytoplasma-responsive miRNAs, only one study has explored in Mexican lime (Citrus aurantifolia L.) infected by ‘Candidatus Phytoplasma aurantifolia’15. Though some of the pathogen-responsive miRNA families are deeply conserved among various plant species, their function may be specie specific16. Moreover, every plant has its own specie-specific miRNAs which play important roles in the regulatory networks associated with the stress resistance, and individual miRNAs of the same family may be expressed differentially and have different functions in response to the same stress13,16. Therefore, the phytoplasma responsive miRNAs have not been fully explored, and how their gene expression reprogrammings achieved are largely unknown. As far as we know, there was no miRNA information reported in mulberry so far. Characterization of the phytoplasma responsive miRNAs in mulberry may provide a novel platform to better understand the plant-phytoplasma interactions.
SiRNAs are distinguished with miRNAs by the structure of their precursors and by their targets, and they arise from a long double-stranded RNAs and typically direct the cleavage of transcripts to which they are completely complementary17. Though thousands of endogenous siRNAs have been sequenced18,19,20, their biological roles are largely unknown except for the functions of some chromatin-associated siRNAs in DNA methylation and transcriptional gene silencing and the roles of transacting siRNAs (ta-siRNA) in plant development and hormone signaling21. Except for the nat-siRNAATGB2, which was specifically induced by the effector avrRpt2 of Pseudomonas syringae pv. tomato DC3000 (Pst DC3000)21, the AtlsiRNA1, which was strongly induced by Pst (avrRpt2), and some siRNAs were found differentially expressed in response to Heterodera schachtii infection22, limited siRNAs involved in plant biotic stress has been reported. Therefore, identification of entire set of siRNAs and their expressions in response to biotic stress is needed to unravel the complex siRNA-mediated regulatory networks controlling gene expression in plant-microorganism interactions.
In the present study, based on the transcriptome information attained, we employed the next generation high throughput sequencing technology to characterize the miRNAs and siRNAs in mulberry. By comparing the normalized expression levels of sRNAs from the healthy and infected leaf libraries, the phytoplasma-responsive miRNAs and siRNAs were identified. In addition, the potential target genes of the responsive miRNAs were identified by sequencing of the degradome library and their functions involved in the response of mulberry to phytoplasma infection were discussed.
To examine the phytoplasma-responsive sRNAs, two sRNA libraries were constructed from phytoplasma-infected leaves and healthy leaves, respectively, and were then subjected to Solexa deep sequencing. The resulting raw sequence reads were cleaned and yielded 28 208 781 and 26 814 061 transcriptome-matching reads from the healthy and phytoplasma-infected leaf libraries, respectively. More than 60% of these mapped small RNAs were 20–24 nt in length with 24 and 21 nt as the major size groups (Fig. 1), consistent with the size of Dicer-like cleavage products. The matched sRNA sequences were clustered into several RNA categories such as known miRNAs, repeats, rRNAs, tRNAs, snRNAs, snoRNAs and unannotatated sRNAs. The known miRNAs account for 41.25% of all sequence reads in the healthy leaf library and 41.85% in the phytoplasma-infected leaf library. However, the proportion of sRNA unique sequences derived from known miRNAs represented only a very small fraction (0.64% and 0.79%, respectively) of the total unique sequences (Table 1). The highest proportion of unique sequences was unannotatated sRNA sequence, which probably include some novel miRNA candidates and lots of siRNAs as well as other type RNAs.
When the sRNA sequences obtained from our libraries were aligned to plant miRNAs in the miRBase database, 164 miRNA members showed perfect matches to the known miRNAs. These miRNAs identified have been shown conserved in many plant species. These conserved miRNAs varied significantly in expression abundances in both libraries (Table 2), and a total of 62 conserved miRNAs identified were found to be significantly responsive to phytoplasma infection, among which, 37 miRNAs were significantly decreased and 25 miRNAs were increased in the infected leaf (P < 0.05, fold 2.0). These differentially expressed miRNAs included some highly expressed miRNAs such as mul-miR156a-5p and mul-miR157a-5p, also some low abundance miRNAs such as mul-miR394a-5p and mul-miR4343a. Moreover, many differentially expressed miRNA-3p sequences, such as mul-miR1023b-3p, mul-miR1157-3p and mul-miR157d-3p, were identified, and these miRNAs (3p) may be the authentic miRNAs since the corresponding miRNA-5p sequences were not detected (Table 2). Real-time PCR analysis for 24 individual miRNAs covering different expression patterns were performed, and the results demonstrated that there was a very strong correlation between PCR data and read frequencies, which indicated that the sequencing profiles are quantitative and reliable (Fig. 2A).
All the non-annotated clean tag sequences were mapped to our mulberry transcriptome database, and only 23 sequences perfectly matched transcriptome sequences able to fold into hairpin structures. These sequences comprise 20 putative new miRNA families (Table 3). The lengths of novel miRNAs identified vary from 20 to 23 nt with a peak at 21 nt, and most begin with a 5′ uridine which is a characteristic feature of miRNAs. These newly identified mulberry miRNA precursors have negative folding free energies (from −21.0 to −105.0 kcal mol−1 with an average of about −46.9 kcal mol−1) according to Mfold, which are similar to the free energy values of rice and Arabidopsis23 and much lower than the folding free energies of tRNA (−27.5 kcal mol−1) and rRNA (−33 kcal mol−1)24. Of these miRNAs, 3 candidates contained both miRNA-5p and miRNA-3p sequences. Moreover, other 5 miRNA-3p sequences were also detected (Table 3). The detection of miRNA-3p is a strong clue for the formation of precursor hairpin structures and represents further evidence supporting them as novel miRNAs. Compared to the conserved miRNAs, most of the novel miRNAs exhibited relatively low abundance as indicated by their frequencies (Table 2, ,4).4). Among the 23 novel miRNAs, only 7 had at least 20 transcripts per million (TPM) in both libraries with the highest abundance 506 TPM for mul-miRn12-5p. Therefore, the expression levels of the majority of the novel miRNAs identified were low.
In order to identify the novel miRNAs involved in response to phytoplasma-infection, normalized expression profiles of the novel miRNAs in healthy and infected leaves were compared (Table 4). Among 13 differential novel miRNAs, 7 miRNAs were significantly decreased and 6 miRNAs were up-regulated more than 2.0 fold changes (P < 0.05) in the infected leaf. Real-time PCR analysis to determine the expression patterns of the 13 differential miRNAs (Fig. 2B) reveals a very strong correlation between qRT-PCR data and read frequencies, demonstrating that the profiles of these novel miRNAs detected by Illumina sequencing are reliable.
To understand the functions of phytoplasma-responsive miRNAs, target genes of the phytoplasma-responsive miRNAs including conserved and novel miRNAs were predicted. In total, 42 differentially expressed miRNAs were predicted to match 122 target genes (Table 5). The predicted miRNA targets were further experimentally verified by mRNA degradome sequencing. In total, we obtained 25 850 468 raw reads from the degradome libraries. After removing the reads without the adaptor, 25 721 355 clean reads were obtained, among which the percent of 20–21 nt reads was 99.39%. The distinct sequences were aligned to the mulberry transcriptome database, and 6913 912 sequences could be matched to the transcriptome without mismatch. The matched sequences were further analyzed by the CleaveLand pipeline, and 69 target genes of 28 phytoplasma-responsive miRNAs were qualified (Table 5), of which 8 genes were cleaved by two unconserved miRNAs and 61 genes were cleaved by 26 conserved miRNAs. Among these identified targets, 27 belonged to category 0, 5 were in category I, 13 were in category II, 15 were in category III and 9 were in category VI. These results indicated that most of these targets were efficiently cleaved by miRNA. The representative ‘target plots’ (tplots) of identified targets included all five categories were showed in Figure 3.
In most cases, the identified miRNAs were predicted to cleave two or more different targets. For example, three members of SPL transcription factor family genes were predicted to be cleaved by mul-miR157a-5p, and mul-miR160a-5p was predicted to slice three genes belonging to the auxin response factor. Inconstantly, some miRNAs were only predicted to cleave one target such as mul-miR1030a, mul-miR164c-3p, mul-miR6463, and mul-miR6475. Unfortunately, we could not detect the cleavage signature for most of miRNAs in this degradome library. The low identification frequency might be because that their target mRNAs have not been annotated in our mulberry transcriptome data.
The expressions of the targets of phytoplasma-responsive miRNAs were investigated by real-time PCR in response to phytoplasma-infection (Fig. 4). The results indicated that most of the target genes showed a negative correlation with the corresponding miRNA expression, and this is consistent with miRNA function in guiding the cleavage of target mRNAs. But it should be noted that a similar expression pattern between the target gene and its corresponding miRNAs was also observed. For instance, mul-miR169a-3p was positively correlated with its target gene encoding phosphatidylinositol 4-kinase. Of course, the putative target genes were bioinformatically predicted, and are subjected to experimental verification.
To better understand miRNA functions, we subjected the identified target genes to Gene Ontology (GO) analysis. The result of GO analysis demonstrated that the target genes of the phytoplasma-responsive miRNAs were classified into 10 categories according to their ontologies in Arabidopsis based on KEGG functional annotations (Fig. 5). The first category of the predicted target genes encoding an array of transcription factors involved in regulation of gene expression. In the second category, target genes associated with signaling pathway, indicating that there were many signaling pathways involved in phytoplasma infection. Target genes involved in metabolic process belong to the third category, suggesting that phytoplasma infection may change diverse metabolic processes in infected plant. The other categories included stress response, development, DNA and RNA processing, and so on. In most cases, targets of the same miRNA were with similar functions (Table 5), however, some miRNAs, such as mul-miR394a-5p, target many genes with different functions. So these miRNAs might be play potential roles in the expression control of genes related to diverse biologic processes, and the regulation networks of miRNAs involved in the response of mulberry to phytoplasma-infection were intricate.
According to the standard that siRNAs were a pair of perfectly complementary sRNAs with 2 nt overhangs at the 3′-end, 315 910 and 271 367 unique siRNA sequences were identified in the healthy and infected leaf small RNA libraries, respectively. The siRNAs identified varied in length from 20 to 24 nt, and the 24 nt siRNAs, which were considered as long siRNAs, were the major size groups (Fig. 6). In order to determine the siRNAs responsive to phytoplasma infection, we looked for differential siRNAs by comparing the normalized expression levels of the siRNAs between the healthy and infected leaf libraries. A total of 14 598 siRNAs were identified to be significantly responsive to phytoplasma infection. Among these differential siRNAs, 10 745 siRNAs were significantly decreased and 3853 siRNAs were up-regulated in the infected leaf (P < 0.05, fold 2.0) (Fig. 6; Supplementary Text).
It was reported that Arabidopsis miR393 plays a role in plant PAMP-triggered immunity in Arabidopsis10. Though phytoplasma lacks cell wall and flagellum, and has no PAMPs like flg224,6, our data showed that mul-miR393 was also induced by phytoplasma infection (Table 2). Therefore, mul-miR393 was induced by phytoplasma through some pathway other than flg22. MiR393 is potentially encoded by two genes, AtMIR393A and AtMIR393B, in Arabidopsis18. To study mul-miR393 function in the response to phytoplasma infection, we cloned the two genes encoding mul-miR393, MulMIR393A and MulMIR393B, in mulberry. Sequence analysis showed that limited homology exists between the MulMIR393A and MulMIR393B, but the mature mul-miRNA393-5p sequences produced are identical (Fig. 7A). Meanwhile, the promoters of MulMIR393A and MulMIR393B were cloned, and cis-acting regulatory elements analysis revealed that the two promoters were different, such as the promoter of MulMIR393A has a cis-acting element involved in salicylic acid responsiveness, but this was not found in the promoter of MulMIR393B. GUS activity analysis indicated that MulMIR393A::GUS showed induced GUS activity, but in contrast MulMIR393B::GUS showed no induced GUS activity after application of salicylic acid (SA) (Fig. 7B), confirming the cis-acting regulatory elements analysis results above. Therefore, MulMIR393A and MulMIR393B have different regulators. Previous studies showed that miR393 is a developmentally regulated miRNA which arises predominantly from AtMIR393B in aerial organs of Arabidopsis25. But when the plants were treated with flg22, the expressed level of miR393 was up-regulated mainly from activated transcription of AtMIR393A10. Northern blot and real-time PCR analysis results showed that the transcript level of MulMIR393A was increased in the infected leaves, but the transcript level of MulMIR393B was unaltered (Fig. 7C, D). Therefore, the mul-miR393 accumulation resulted from the increased transcripts of MulMIR393A in responding to phytoplasma infection.
So far, miRNAs and their regulatory functions have been extensively characterized in model plant species, but to our knowledge, it was the first report about the sRNAs of mulberry. In present study, solexa high-throughput sequencing of small RNAs revealed a diverse and complex small RNA population. Among the 62 conserved phytoplasma-responsive miRNAs identified in mulberry, only 6 miRNAs, including miR157a, miR172a-3p, miR4414, miR473, miR477a, and miR479, were also differentially expressed in Candidatus Phytoplasma aurantifolia infected Mexican lime (Citrus aurantifolia L.)15. Therefore, the effect of phytoplasma infection on the expression of miRNAs depends on the pathogen and host types. Furthermore, it was found that the miRNAs of the same family were differentially regulated in the response to phytoplasma-infection. For instance, the expression level of mul-miR160a-5p was up-regulated by phytoplasma-infection, but the expression level of mul-miR160b-3p showed no obvious change in the infected leaves (Table 2). So the gene expression programmings of the miRNAs, which share a high degree of sequence similarity and belong to the same family, were controlled by different regulator16.
Phytoplasma-associated disorders are believed to be the auxinic diseases of plants26. It was reported that many auxin response factors were down-regulated by the effector (tengu) of Onion yellows phytoplasma, and suggested that the symptoms caused by phytoplasma may be involved in an auxin-related pathway5. In plants, a group of miRNAs regulate and fine-tune multiple hormone signaling pathways coordinately27. Our results showed some miRNAs, which target the genes associated with auxin signaling and metabolism, were differentially expressed between the healthy and infected leaves. MiR393 was the first miRNA identified which targets auxin receptor mRNAs and subsequently represses auxin signaling in Arabidopsis10, and we found that mul-miR393 was up-regulated in the infected leaves (Table 2). In addition to miR393, target analyses showed that mul-miRn10-3p was also involved in auxin signaling pathway by regulating auxin response factor (Table 5). Moreover, miR160, miR164 and miR166, which were associated with auxin signaling pathway by regulating different transcription factors13, which were found differentially expressed in healthy and infected leaves (Table 2). Furthermore, the differential mul-miR169a-3p and mul-miR6204 were predicted to target the genes coding for O-fucosyltransferase family protein and SAUR-like auxin-responsive protein family, respectively, which were associated with auxin metabolism (Table 5). Therefore, these differentially expressed miRNAs may play a role in altering auxin signaling and metabolism in infected plants.
Beside the auxin signaling pathway, miRNAs were also involved in regulating other hormone signaling pathways. For instance, miR164, miR166 and miR398 were reported involved in abscisic acid signaling pathway28,29,30. MiR166 and miR172 were found to be down-regulated by gibberellin31 and cytokinin32, respectively. The expression levels of these miRNAs were differential between healthy and infected leaves (Table 2). In addition, mul-miR2111a-5p was predicted to target the gene coding ACC oxidase 1 associated with ethylene metabolism, and mul-miR2595 and mul-miRn20-5p were predicted to target the pentatricopeptide repeat-containing protein gene involved in the response to SA stimulus (Table 5). Therefore, it is more likely that the phytoplasma infection-induced alteration in hormonal signaling triggers reprogramming of the plant growth and developmental pattern, leading to some symptoms in infected plants. This is in accord with early reports that phytoplasma infection disturbs hormonal balance in the host plants4,33,34,35,36.
Symptoms induced in infected plants suggest that phytoplasma infection may modulate developmental processes within the plant host37. This was confirmed by the reports that the expression of effector protein (phyllogen) of phytoplasmas in Arabidopsis can alter floral development, resulting in the production of leaf-like flowers38,39. In Arabidopsis, many Squamosa Promoter-Binding Protein-Like (SPL) genes are post-transcriptionally regulated by miR156. Meanwhile, some SPL genes, such as AtSPL9, in turn positively regulate the expression of miR172. This forms the miR156-AtSPL9-miR172 regulatory pathway40. Our results showed that mul-miR156 was up-regulated and its target gene SPL9 was down-regulated, meanwhile the mul-miR172 was down-regulated (Table 2; Fig. 2A). Therefore, the miR156-SPL9-miR172 regulatory pathway may be also conserved in mulberry and responsive to phytoplasma infection. SPL family of transcription factors play important roles in flower and fruit development, plant architecture and phase transition40. AP2 domain-containing transcription factor, which is the target of mul-miR172, takes part in regulating flowering time and floral organ identity41. Therefore, the expression changes of miR156 and miR172 might disturb the plant development process and be responsible for the mulberry yellow dwarf symptoms such as development of green leaf-like structures instead of flowers and sterility of flowers.
In addition, our results showed some other miRNAs, which target the genes associated with plant growth and development, were also differentially expressed. For example, the expression level of miR166h-3p, which targets the class III HD-Zip protein 6 (Table 2, ,5)5) involved in the regulation of shoot meristem initiation, vascular development42, was down-regulated by phytoplasma infection, and this probably resulted in an abnormal plant development such as meristem defects and was partially responsible for the symptoms such as dwarf and witches' broom. Meanwhile, miR160, miR162, miR164, miR165, miR168, and miR319 were reported to modulate plant development through different pathway41,43,44,45,46,47,48. These up or down-regulation of these miRNAs (Table 2) may disorder the expressions of many genes involved in diverse development processes and change the architecture of infected plants causing some symptoms in the infected plants.
Our data showed that some differentially expressed miRNAs target the genes associated with CHO metabolism, protein metabolism, lipid metabolism, and so on. For example, the up-regulated mul-miR160a-5p may down-regulate its target, beta-amylase gene (Table 2, ,5),5), and responsible for the accumulation of starch in the infected leaves. Similarly, the up-regulated mul-miR529b may down-regulate its target, chlorophyll synthase gene (Table 2, ,5),5), resulting in the blockage of the biosynthesis of chlorophyll and the decoloration symptoms. Thus, the differential expressions of these miRNAs might disturb the normal metabolism process and partly responsible for some of the mulberry yellow dwarf symptoms. Based on the information obtained, the complex regulatory networks of miRNAs involved in the response of mulberry to phytoplasma-infection is outlined in Figure 8.
Unlike miRNAs, the lengths of siRNAs are more variable, ranging from 21 to 30 nt. The biogenesis of the siRNAs depends on specific RNA-dependent RNA polymerases and Dicer-like (DCL) proteins and RNA polymerase IV49. The siRNAs produced by different DCLs may have different lengths and functions50. In addition to the relatively limited number of miRNAs, thousands of differentially expressed siRNAs were involved in the response of mulberry to phytoplasma-infection. The lengths of these siRNAs are more variable (Fig. 6), and the targeted genes of these differentially expressed siRNAs were associated with a wide range of functions in hormone networks, development and metabolism. It was reported that the 22-nt forms of miRNAs are important for triggering secondary siRNA production in plants51,52. Previous research showed that miR393 is necessary to initiate the production of the secondary siRNAs from the transport inhibitor response 1/auxin signaling F-box 2 (TIR1/AFB2) clade of auxin receptors (TAARs) transcripts in Arabidopsis25. Sequencing results showed that mul-miR393 is 22-nt long, so we wondered whether it can trigger the production of the secondary siRNAs. Our small RNA library sequencing data showed there were some siRNAs matching to the TIR1 mRNA downstream of the mul-miR393-5p binding sites (Figure 9). Thus, mul-miR393 might also trigger the production of the secondary siRNAs. However, we can't predict which one is necessary to trigger the production of these siRNAs, mul-miR393a or mul-miR393b. Future studies will be carried out to explore the functions of these phytoplasma-responsive siRNAs, and identify the pathways connected with the response against phytoplasma-infection. In conclusion, the information provided here would be particularly useful for understanding the function of sRNAs in mulberry and be aid to reveal the mechanisms underlying of phytoplasma pathogenicity and the symptoms they cause in host plants.
One-year-old cutting seedlings derived from the same mother tree, Husang 32 (Morus multicaulis Perr.), were used as rootstocks which were grafted with the scions collected from healthy and phytoplasma-infected mulberry trees (Husang 32), respectively. All the grafts were incubated in a growth chamber at 26°C, humidity 90% and under 12 h of light. Six weeks after infection, the plant showing typical symptoms as yellowing of the leaves, stunting, and witches'-broom were confirmed by PCR assay with an amplified fragment of the 16S rRNA gene (GenBank Accession No. EF532410) and TEM detection as described before1. The plants that were phytoplasma-positive by PCR assay and TEM detection were used as treatments while the plants grafted with the scions collected from healthy trees were used as controls.
Leaves were collected from the infected plants and the controls, and RNAs were isolated from the leaves of healthy and infected plants separately using TRIzol kit (Invitrogen) following the manufacturer's instructions. Total RNAs were separated on 15% denaturing PAGE and small RNAs between 14 and 30 nt were enriched and ligated with 5′ and 3′ adaptors. The ligation products were purified and amplified after reverse transcription. The amplification products were purified and used for sequencing by Solexa machine.
The raw sequences tags from HiSeq sequencing were cleaned by getting rid of adaptors, low-quality tags, as well as contaminants to generate clean reads. Length distribution of clean tags is then summarized, and the trimmed sequences longer than 18 nt were used for further analyses. The remaining clean reads were annotatated into different categories through standard bioinformatics analysis. Subsequently, the sequences were compared with a miRBase database v19.0 (http://www.mirbase.org/). The sequences with identical or related sequences from other plants were identified as conserved miRNAs, and the other sequences were analyzed for novel miRNA candidates.
To discover novel miRNAs, all non-annotated clean tag sequences were mapped to our mulberry transcriptome database. The characteristic hairpin structure of miRNA precursor was used to predict novel miRNA using the Mireap algorithm (http://sourceforge.net/projects/mireap/) by exploring the secondary structure, the dicer cleavage site, the minimum free energy, and other criteria provided by before53.
To identify siRNAs, tags from clean reads were aligned with each other. A pair of perfectly complementary sRNAs with 2 nt overhangs at the 3′-end were considered as siRNAs.
The expression of miRNAs and siRNAs in healthy and infected leaf samples were normalized to obtain the expression of transcripts per million, and only the sRNA with an expression value of ≥1, it was considered in future differential expression analysis. Based on the normalized expression analysis, the fold-change and P-value were calculated according to the following equation. Normalized expression = (Actual sRNA sequencing reads count/Total clean reads count) ×1,000,000; Fold change = Log2 (x/y); P-value formula:
The N1 and N2 represent total count of clean reads of a given sRNA in the sRNA library of phytoplasma infected and healthy leaves, respectively. The x and y represent normalized expression level of a given sRNA in the sRNA library of phytoplasma infected and healthy leaves, respectively.
The sequences from the mulberry transcriptome data were chosen to predict the miRNA targets, and the plant small RNA target analysis server, psRNATarget (http://plantgrn.noble.org/psRNATarget/), was used to predict the target genes following the rules suggested by Allen et al.54. All the predicted miRNA targets were experimentally verified by mRNA degradome sequencing. A degradome sequencing library was constructed using mRNAs isolated from the healthy and infected leaves according to the protocols described previously55,56 and then were submitted for parallel analysis of RNA end (PARE) sequencing on a Solexa Analyzer. All the adaptor sequences and low quality sequencing reads were removed from the raw sequencing reads to obtain the clean reads with the length of 20 and 21 nt, which were used to identify potentially cleaved targets with the software package, CleaveLand 3.0, as previously described55,56. The clean reads were mapped to the mulberry transcriptome unigene sequences obtained, and only the perfect matching reads were kept and extend to 35–36 nt by adding 15 nt of upstream of the sequence. All resulting reads were reverse-complemented and aligned to the miRNA identified and scored by previously described method54. All targets identified were classified into five categories based on the abundance of the resulting mRNA tag relative to the overall profile of degradome reads that matched the target57. When a target mRNA was predicted or identified, a BlastN search against a reference Arabidopsis thaliana database downloaded from TAIR (http://www.arabidopsis.org/; version TAIR10) was used to provide the gene ontologies, and then GO analysis was performed for the matched Arabidopsis accession entries based on their TAIR GO categories.
RNA was extracted using the TRIzol® reagent following the manufacturer's recommendations (Invitrogen) and digested with DNase I. Real-time PCR was performed using the PrimeScript™ miRNA qPCR Starter Kit Ver.2.0 for miRNAs and the SYBR Premix Ex TaqTM kit for mRNA according to the manufacturer's protocol on the Rotor-Gene 3000A system. The U6 and actin genes were amplified as reference genes for miRNA and mRNA normalization, respectively. The primers (Supplementary table 1, 2) used to amplify the genes were designed based on our available mulberry transcriptome data. Comparative cycle threshold (Ct) method58 was used to evaluate the relative gene expression level. All samples were assayed in triplicate.
Twenty micrograms of total RNA of each sample isolated as above were separated on 1.2% formaldehyde denatured agarose gel, and then were blotted onto nylon Hybond N membrane (Hybond N, Amersham). The blots were hybridized with digoxigenin-labeled RNA probes prepared using the PCR DIG Probe Synthesis Kit (Roche, Germany). Prehybridization, hybridization, membrane washing, and detection were carried out as described previously59.
The promoter sequences of MulMIR393A and MulMIR393B genes were cloned using Tail-PCR method and analyzed using the program online (http://bioinformatics.psb.ugent.be/webtools/plantcare/html/) to search for cis-acting motifs. Then they were cloned into the vector pBI121, which harbors the promoter 35S and the GUS gene, to replace 35S, respectively. The pBI121 vector containing MulMIR393A::GUS or MulMIR393B::GUS was introduced into Agrobacteriun tumefaciens strain GV3101. Tobacco (Nicotiana benthamiana) plants were used for transformation by infection of leaf sections with A. tumefaciens strain harboring the constructs described above. Selected transgenic shoots were treated with SA by dipping the leaves into a solution containing 1 mM SA, and the leaves treated with water were used as control. Histochemical staining for GUS activity was performed 24 h later using the method described by Jefferson60 with X-gluc as a substrate.
This work was funded by the national natural science foundation of China (No. 30972366, 31070573, 31100478) and science foundation for the excellent youth scholars of Shandong province of China (No. BS2009NY024, BS2010NY015).
The authors declare no competing financial interests.
Author Contributions Y.P. and X.L. were responsible for design, sample preparations, and contributed to write the manuscript. Y.Q. and F.Y. were responsible for high-throughput deep sequencing analysis. C.Z. and Y.Y. were responsible for PCR analysis. H.L. and H. were responsible for mul-miR393 analysis. All authors edited the manuscript and reviewed the manuscript.