- Research article
- Open Access
Identification and profiling of miRNAs during herbivory reveals jasmonate-dependent and -independent patterns of accumulation in Nicotiana attenuata
BMC Plant Biology volume 12, Article number: 209 (2012)
Plant microRNAs (miRNAs) play key roles in the transcriptional responses to environmental stresses. However, the role of miRNAs in responses to insect herbivory has not been thoroughly explored. To identify herbivory-responsive miRNAs, we identified conserved miRNAs in the ecological model plant Nicotiana attenuata whose interactions with herbivores have been well-characterized in both laboratory and field studies.
We identified 59 miRNAs from 36 families, and two endogenous trans-acting small interfering RNAs (tasiRNA) targeted by miRNAs. We characterized the response of the precursor and mature miRNAs to simulated attack from the specialist herbivore Manduca sexta by quantitative PCR analysis and used ir-aoc RNAi transformants, deficient in jasmonate biosynthesis, to identify jasmonate-dependent and -independent miRNA regulation. Expression analysis revealed that groups of miRNAs and tasiRNAs were specifically regulated by either mechanical wounding or wounding plus oral secretions from M. sexta larvae, and these small RNAs were accumulated in jasmonate-dependent or -independent manners. Moreover, cDNA microarray analysis indicated that the expression patterns of the corresponding target genes were correlated with the accumulation of miRNAs and tasiRNAs.
We show that a group of miRNAs and tasiRNAs orchestrates the expression of target genes involved in N. attenuata’s responses to herbivore attack.
A group of non-coding small RNAs (smRNAs) plays an important role in transcript regulation by binding to their target sequences, resulting in transcriptional degradation, transcriptional or translational inhibition of the targets [1–5]. The smRNAs are classified into two major classes: microRNAs (miRNAs), and small interfering RNAs (siRNAs). Primary transcripts of miRNAs are processed into precursors of miRNAs that form secondary stem-and-loop structures, which are processed by the ribonuclease DICER-like 1 (DCL1) into miRNA/miRNA* duplexes, which are subsequently incorporated into the RNA induced silencing complex (RISC) [5–7]. The siRNAs are further classified into trans acting siRNAs (tasiRNAs), chromatin-associated cis acting siRNAs, and natural antisense siRNAs, based on their biogenesis [8–10]. Biogenesis of tasiRNAs is regulated by miRNAs, which direct cleavage of primary tasiRNA (TAS) transcripts encoding tasiRNAs, resulting in second-strand RNA synthesis by RNA-dependent RNA polymerases (RdRs). The double-stranded RNAs are diced by DCL4 to generate tasiRNAs in Arabidopsis thaliana[9–11].
Plant miRNAs and siRNAs are involved in several developmental processes : embryogenesis , organ polarity , leaf formation , root development [11, 16], phytohormone signaling [17, 18], and flowering time . Plant defense signaling is also regulated by miRNAs in response to different abiotic stresses [20, 21] including heat, cold, drought [22, 23], and UV-B radiation . For example, A. thaliana miR399 (Ath-miR399), induced during phosphate starvation, targets the ubiquitin-conjugating E2 enzyme involved in phosphate uptake from the soil . Under drought stress, Ath-miR159 regulates MYB33 and MYB101 transcription factors, which activate abscisic acid responses during seed germination . Ath-miR398 regulates Cu/Zn-superoxide dismutase genes, which detoxify superoxide radicals . A recent study reported that several miRNAs are induced upon mechanical wounding in tobacco leaves and roots .
Plant miRNAs are also involved in biotic interactions. Ath-miR393 is induced by flagellin-derived PAMP peptide 22, and targets the F-box protein and transport inhibitor response 1, which plays a key role in antibacterial responses . Ath-miR160, Ath-miR167, and Ath-miR825 are induced in response to infection by Pst DC3000 hrcC , and A. thaliana and Nicotiana tabacum plants infected by TYMVp69 virus accumulate high levels of miR156, miR160, and miR164 [30, 31]. Plant miRNAs are also involved in beneficial interactions with bacteria: miR482, miR1512, and miR1515 play a role during rhizobial infection in Glycine max nodulation with Bradyrhizobium japonicum. However, little is known about the role of plant miRNAs in the response to insect herbivores.
The wild tobacco Nicotiana attenuata and its herbivore community have become an ecological model system for the study of plant-herbivore interactions. During attack by insect herbivores, N. attenuata rapidly induces jasmonate-mediated defense responses, which reconfigure primary and secondary metabolism [33, 34]. Jasmonates comprise jasmonic acid (JA), its derivatives and conjugates; the jasmonates and in particular, the active hormone jasmonoyl-isoleucine (JA-Ile) regulate most defenses against chewing herbivores . Fatty acid amino acid conjugates (FACs) in oral secretion (OS) from larvae of the specialist herbivore, Manduca sexta, trigger jasmonate-mediated direct and indirect defenses in N. attenuata, such as nicotine accumulation, proteinase inhibitor production, diterpene glycoside biosynthesis, and emission of green leaf volatiles [36–38]. Transgenic plants impaired in jasmonate biosynthesis or signaling show increased susceptibility to herbivory in both glasshouse and field studies [35, 39–41].
OS-elicitation dramatically changes the smRNA population in N. attenuata, and two major components of the smRNA pathway, RdRs and DCL proteins, function in biotic and abiotic stress responses [42–45]. Silencing of N. attenuata RdR1, DCL3, and DCL4 results in impaired defense responses against M. sexta herbivory [42, 44, 45]. Silencing either NaRdR1 or NaDCL4 impairs jasmonic acid (JA) accumulation, and co-silencing NaDCL3 and NaDCL4 reduces JA levels, indicating that RdR1/DCL4-mediated smRNAs are critical regulators of responses to insect herbivory.
To deepen our understanding of the roles that smRNAs play in plant-insect interactions, we identified primary miRNA (MIR) transcripts and TAS transcripts encoding tasiRNAs in a transcriptome database of N. attenuata, and computationally analyzed secondary stem-and-loop structures of MIR transcripts. To understand the role of jasmonates in regulating miRNAs, we examined miRNA accumulation in jasmonate-deficient allene oxide cyclase (AOC) RNAi lines. The AOC protein provides a precursor for JA biosynthesis . Expression analysis of miRNAs and tasiRNAs with their putative target genes provides evidence for a key role of plant smRNAs in the response to herbivory.
Results and discussion
Identification of conserved miRNAs and their precursors in N. attenuata
To identify conserved miRNAs in N. attenuata, we used a 454-transcriptome database of N. attenuata to conduct a BLAST search against conserved plant miRNAs in the miRBase (http://www.mirbase.org) (Figures 1A and 1B). This search identified 59 potential miRNAs distributed in 36 families (Table 1). We used the BLASTX algorithm against the NCBI protein database to check that the putative primary transcripts of miRNAs were non-coding. Web-based mFOLD software (http://mfold.rna.albany.edu/) was used to predict secondary stem-and-loop structures. Of the identified miRNA-precursors, 52 had stem-and-loop structures (Figure 1C and Additional file 1), which were created with minimum free energies (MFE) ranging from ΔG = −97.5 kcal mol-1 to −33.3 kcal mol-1 (Table 1) with an average MFE of −62.1 kcal mol-1. This average MFE is comparable to that found in A. thaliana (−59.5 kcal mol-1), higher than in the red alga Porphyra yezoensis (−41.7 kcal mol-1) and lower than in the monocots rice (−71.0 kcal mol-1) and wheat (−72.4 kcal mol-1) [22, 48]. Only seven predicted miRNA-precursors transcripts did not form stem-and-loop structures or were not stable (Table 1). We identified several N. attenuata (Nat) miRNA families (Nat-miR403, Nat-miR478, Nat-miR482, Nat-miR1128, Nat-miR1133, Nat-miR1446, Nat-miR1863, Nat-miR2911, and Nat-miR5281) which were not reported in N. tabaccum[28, 49]. Among these, Nat-miR478, Nat-miR482, Nat-miR1128, Nat-miR1133, Nat-miR1446, Nat-miR1863, and Nat-miR5281are absent in A. thaliana but are close homologues to those in other plant species (Table 1).
Next, we designed probes to detect N. attenuata miRNAs on RNA blots ( Additional file 2). We performed northern blot hybridization using 40 μg of total RNA extracted from rosette leaves to detect selected miRNAs. Accumulation of miRNAs varied ( Additional file 3). Accumulation of Nat-miR159, Nat-miR171, Nat-miR172, and Nat-miR319 was high compared to Nat-miR157, Nat-miR393, Nat-miR396, and Nat-miR828 in leaves from rosette-stage plants. For further analyses of precursor and mature miRNA abundance, we used real-time quantitative PCR (qPCR) with specific primer sets ( Additional file 4 and 5).
Identification of conserved tasiRNAs in N. attenuata
Four families of endogenous tasiRNAs (TAS1, TAS2, TAS3, and TAS4) identified in A. thaliana are regulated by miRNAs [9, 50]. We found three TAS3 transcripts and one TAS4 transcript in N. attenuata (Figure 2), and constructed a phylogenetic tree with their homologs from different plant species to examine the evolutionary relationships of TAS3 expressed in dicotyledonous and monocotyledonous plant species [50, 51]. Not surprisingly, NaTAS3 members were grouped amongst members of the dicotyledonous plant species (Figure 2A).
Because TAS3 and TAS4 transcripts contain the binding sites of miR390 and miR828, respectively [9, 11, 52], we blasted Nat-miR390 and Nat-miR828 against the TAS transcripts (Figure 2B), and compared the N. attenuata tasiRNAs with those in A. thaliana (Figure 2C). Sequence analysis indicated that binding sites of miRNA and tasiRNA sequences are highly conserved in N. attenuata and A. thaliana.
Wound- and OS-inducible miRNAs
Next, we examined the abundance of miRNAs and their target genes in leaves of rosette-stage WT plants changed after wounding and the application of diluted M. sexta OS (W+OS) or water to wounds (W+W, as a control for W+OS) versus no treatment (control) (Figures 3 and 4). W+OS treatment faithfully mimics the majority of responses elicited by M. sexta feeding [36, 38], but has the distinct advantage of allowing the time of elicitation to be precisely controlled (as opposed to the sporadic nature of M. sexta larval feeding behavior) and hence greatly increases the reproducibility of transcriptional analyses.
The abundance of mature Nat-miR159, Nat-miR160, Nat-miR167, Nat-miR396, Nat-miR403, Nat-miR408, and Nat-miR828 increased after both W+W and W+OS treatments, abundances of which did not differ (Figure 3A). We categorized these miRNAs wound-inducible miRNAs. Twelve conserved miRNA families (Nat-miR156, Nat-miR164, Nat-miR166, Nat-miR168, Nat-miR171, Nat-miR172, Nat-miR319, Nat-miR390, Nat-miR393, Nat-miR394, Nat-miR398, and Nat-miR1446) were significantly increased in W+OS treated leaves compared to control and W+W treated leaves (Figure 3B) and these were classified as OS-inducible miRNAs.
It is well-established that miRNAs and tasiRNAs bind to their target mRNAs by perfect or imperfect complementarity . Such complementarity permits the identification of miRNA targets in databases. We blasted miRNA sequences against an in-house transriptomic database using BLASTN algorithms with default parameters allowing 1–4 mistmatches. We identified 111 potential targets, inluding targets of tasiRNAs ( Additional file 6).
Ath-miR164 is known to negatively regulate ORE1/NAC2 transcription factors which are involved in age-dependent senescence  and mutations in ORE1/NAC2 delays chlorophyll loss in old leaves of A. thaliana. We found a NAC transcription factor containing a Nat-miR164 binding site in N. attenuata (Figure 4). The NAC transcripts were significantly reduced in W+W and W+OS elicited leaves (Figure 4), which correlated with elevated Nat-miR164 levels after OS-elicitation (Figure 3B). This result may explain how M. sexta attack delays senescence in N. attenuata.
In addition, ARF6 and ARF8 were putative targets of wound-inducible Nat-miR167, and their transcripts decreased concurrently with the accumulation of Nat-miR167 (Figure 4). Potential targets of OS-inducible Nat-miR171, GRAS domain transcription factors, were down-regulated in OS-elicited plants (Figure 4). Transcripts of several OS-inducible Nat-miR156 targets annotated as SQUAMOSA promoter binding proteins (SBPs) were significantly down-regulated after W+OS treatment. Abundance of Nat-miR172 increased significantly in W+OS treated leaves (Figure 3B) and consistently, transcripts of its putative targets, AP2-like proteins, were significantly down-regulated in W+W and W+OS treated leaves (Figure 4).
The most significant change in miRNA transcripts after W+OS treatment was for Nat-miR390, which showed a 75-fold increase compared to control leaves and a 30-fold increase compared to W+W treated leaves (Figure 3B). Ath-miR390 cleaves TAS transcripts, resulting in production of tasiRNA3 mediated by the RdR6/DCL4 complex [11, 52]. TasiRNA3 controls the transcription of auxin response factors (ARFs) 2, 3, and 4, which regulate leaf morphology and lateral root growth in A. thaliana[11, 16, 52]. Overexpression of TAS3 leads to an increased number of lateral roots, and the knock-out mutant shows impaired lateral root growth . Abundance of mature Nat-tasiRNA3 in W+OS treated leaves was increased (Figure 5A) and Nat-tasiRNA3 was regulated in a JA-independent manner (Figure 5B). The targets of tasiRNA3, transcripts homologous to A. thaliana ARF2, 3 and 4 were significantly reduced after W+OS elicitation (Figure 4), suggesting that Nat-miR390 could affect the architecture of roots and thereby regulate the production of nicotine, which is synthesized in the roots, or promote tolerance of herbivory by increasing mineral uptake or sugar storage in roots. Our previous study shows that silencing NaDCL4, which in turn reduces the accumulation of tasiRNA3, dramatically impairs root growth and nicotine accumulation in N. attenuata. We validated the microarray data by qPCR analysis ( Additional file 7), and the result showed similar expression patterns ( Additional file 8).
The key role of jasmonates in mediating responses to herbivory is well established and it was not surprising to find that the levels of jasmonates dramatically increased in W+OS treated leaves compared to W+W treatment ( Additional file 9) . Silencing of NaAOC reduces JA accumulation in W+OS treated leaves by 90-100% ( Additional file 9) and is known to silence the production of JA-mediated defenses . OS-elicitation of this genotype allowed us to understand which miRNAs are regulated by jasmonates.
Abundance of 10 miRNAs (Nat-miR159, Nat-miR160, Nat-miR164, Nat-miR166, Nat-miR168, Nat-miR171, Nat-miR172, Nat-miR393, Nat-miR403, and Nat-miR408) increased after W+W and W+OS treatments, but did not differ between WT and ir-aoc plants (Figure 6A). We considered these miRNAs as JA-independent miRNAs. The abundance of Nat-miR156, Nat-miR167, Nat-miR390, Nat-miR396, Nat-miR398, and Nat-miR1446 were significantly higher in both control and W+OS elicited leaves of ir-aoc than of WT plants (Figure 6B), indicating that jasmonates or JA-signaling negatively influences levels of these miRNAs, regardless of treatment. Nat-miR390 abundance was higher in ir-aoc only after W+OS treatment (Figure 6B). While abundance of Nat-miR319 and Nat-miR394 were induced by W+OS treatment in WT (Figure 3B), they were not induced by W+OS treatment in ir-aoc plants (Figure 6C). We considered these as OS-inducible JA-dependent miRNAs. Wound-inducible Nat-miR828 was induced in both W+W and W+OS treatments in WT (Figure 3A), but their levels were not dramatically induced in elicited JA-deficient ir-aoc plant compared to WT (Figure 6C).
One main target of miR319 in A. thaliana is TCP (TEOSINTEBRANCHED/CYCLOIDEA/PCF) transcription factors, which positively regulates jasmonate biosynthesis . The TCP4 protein can bind a promoter of LIPOXYGENASE2 (LOX2) in A. thaliana. We do not yet know the orthologs of A. thaliana TCPs in N. attenuata, but the timing of Nat-miR319 induction was similar with that of W+OS elicited JA accumulation and the induction of NaLOX3, which is the functional homolog of A. thaliana LOX2[40, 55]. W+OS treatment amplifies the wound-induced JA accumulation ( Additional file 9) [36, 38] and does the same for Nat-miR319 (Figure 3B). This suggests that Nat-miR319 could play a role in the fine-tuning regulation of jasmonate biosynthesis  during herbivory.
Nat-miR828 and Nat-tasiRNA4 were increased by W+W and W+OS treatments in WT (Figure 3A), but not in ir-aoc plants (Figures 5 and 6C). Levels of Nat-miR828 dramatically decreased in ir-aoc plants compared to WT and its primary transcript as well (Figure 6C and Additional file 10). One target of miR828 is a TAS4 transcript, which encodes tasiRNA4 processed by RdR6/DCL4 proteins [9, 56]. TasiRNA4 targets several MYB transcription factors which regulate phenylpropanoid biosynthesis [9, 56]. Several phenylpropanoid defense metabolites increase in N. attenuata during herbivory and silencing of DCL4 alters the accumulation of dicaffeoyl spermidine and caffeoylputrescine in W+OS-elicited leaves . A jasmonate-inducible MYB transcription factor, NaMYB8, is involved in plant defense against herbivory and phenylpropanoid biosynthesis [57, 58]. Although we were not able to find a MYB transcription factor containing a Nat-tasiRNA4 binding site in our current cDNA library, sequence conservation of miR828 and tasiRNA4 in A. thaliana and N. attenuata will guide the identification of targets which regulate secondary metabolite production.
Temporal effect of single elicitation on miRNAs and their targets accumulation
To examine the kinetics of miRNAs in response to insect herbivory, we measured the accumulation of miRNAs and their targets at 1 h and 5 h after treatments. Most OS-induced miRNAs increased in abundance within 1 h after OS-elicitation and continue to increase until 5 h, except miR166 and miR390 ( Additional file 11). The miRNA390 increased at 1 h but decreased again by 5 h post-elicitation ( Additional file 11). OS-inducible miRNAs were also induced by wounding but the increase after wounding was less (miR164, miR166, miR172, miR319, and miR398) or slower (miR168, miR390, and miR393) than that after OS-elicitation. Interestingly, most of the miRNA and tasiRNA target genes were reduced at 1 h after OS-elicitation but return to control levels (SBPa, ARF2, ARF3, TAS3a, TAS3b, HD-ZIPIII, TIR1) or even increased to levels greater than prior to elicitation (MYB, ARF6, ARF8), suggesting a complex relationship between the regulator mature miRNAs and their target genes.
Comparison between primary and mature miRNA transcript accumulation
We checked the accumulation of MIR transcripts using previously reported microarray data  ( Additional file 12) and confirmed these responses by qPCR ( Additional file 10). Accumulation of MIR156, MIR159, MIR164, MIR166, MIR167, MIR168, MIR172, MIR393, MIR396, MIR398, MIR403, and MIR1446 was not correlated to the abundance of their mature miRNAs. Changes in transcript levels of MIR160, MIR171, MIR319, MIR390, MIR408, and MIR828 reflected the transcript accumulation of their mature miRNAs. Generally, a weak relationship between MIR and mature miRNA accumulation has been reported . Our data also suggests that the processing of MIR transcripts is more important than the transcription of MIR genes in herbivory-induced miRNA regulation.
In order to build a database of plant miRNAs functionally involved in plant-insect interactions, we investigated the accumulations of miRNAs and their targets in N. attenuata after OS-elicitation. We classified W+OS-induced miRNAs and tasiRNAs into four groups: JA-dependent or –independent wound-inducible and OS-inducible miRNAs or tasiRNAs (Figure 7). Herbivore-attacked N. attenuata plants induce defense metabolites and tune their physiology to tolerate insect attack. This study shows that W+OS elicitation, a rigorous means of mimicking herbivore attack, rapidly changed the expression of miRNAs involved in flowering time, root morphology, senescence, hormone regulation, and metabolite synthesis. In future work, we will experimentally characterize the function of these W+OS regulated JA-dependent and JA-independent miRNAs and their targets.
Plant material and growth conditions
WT and JA-deficient ir-aoc RNAi lines of the 30th inbred generation of the N. attenuata (Torrey ex. Watson) (originally collected in Southwestern Utah, USA) were used for the experiments. Seeds were germinated on MS basal medium supplemented with vitamin B5 (GB5, Duchefa, http://www.duchefa.com). For treatment, the leaf at the +1 node of the rosette (youngest fully-expanded leaf) was wounded with a fabric pattern wheel, and either 20 uL MilliQ water (W+W) or 20 μL M. sexta OS (diluted 5X with sterile water) (W+OS) were applied. Leaves at the same position from untreated plants served as controls. Plants were grown in a glasshouse under 16/8 h (long day) (26-28°C) supplemental light from Master Suns-U PIA Agro 400-W sodium lights.
Identification of miRNAs and precursors, and prediction of miRNA targets
To identify conserved miRNA and tasiRNA, we blasted all conserved plant miRNAs present in miRBase (http://www.mirbase.org) against N. attenuata’s in house 454-transriptome database using default search parameters allowing 1 or 2 mismatches (Figure 1A). Hits of 20 to 24 nt sequences against non-coding transcripts with up to four nucleotide mismatches were selected as candidates to check for predicted miRNA secondary structure. The public web-based mFOLD server (http://mfold.rna.albany.edu/) was used to predict secondary stem-and-loop structures using default parameters, folding temperature (37°C) and ionic conditions (1M NaCl) with minimum free energy (MFE) formations (Table 1). Predicted miRNA-precursors were additionally depicted using RNAshapes (http://bibiserv.techfak.uni-bielefeld.de/rapidshapes)  to create stem-and-loop structures ( Additional file 1). For TAS3 identification we used the following TAS3 orthologs: Saccharum officinarum (EU327139), Sorghum bicolor (EU327137), Zea mays (EU327127, EU293143), Oryza barthii (GQ420228), Triticum aestivum (EU327134), A. thaliana TAS3a (AT3G17185), TAS3b (AT5G49615), TAS3c (At5g57735), Nicotiana tabacum (FJ804751, FJ804743), Nicotiana benthamiana (FJ804742), Datura stramonium (FJ804744). A neighbor-joining tree was built for TAS3 transcripts using MEGA4 with group evaluation and 1000 bootstrap replicates .
RNA extraction and RNA blot hybridization
RNA blot hybridization was performed as described by Molnar et al. (2007) . Extracted total RNA was treated with DNase I (Fermentas; http://www.fermentas.com) according to the manufacturer’s protocol. Enzymes were removed by phenol-chloroform extraction and total RNA was re-isolated by ethanol precipitation. A denaturing 15% polyacrylamid gel containing 7M urea was prepared using the BIORAD Mini-Protean 3 Cell system (http://www.bio-rad.com). Forty μg of total RNA were denatured with an equal volume of 2x gel-loading dye (Fermentas, http:/http://www.fermentas.de) at 65°C for 5 minutes. Denatured RNA samples were loaded into the prerun 15% denaturing polyacrylamid gel at 50 V. A ZR small-RNATM Ladder (Zymo Research, http://www.zymoresearch.com) was also loaded. The gel was run at 150 V until the bromphenol blue in the loading dye reached the bottom of the gel. Next, the portion of the gel containing small RNA (between the bromphenol blue and xylene cyanol) was cut out and blotted onto a nylon membrane (Nytran Supercharge; http://www.whatman.com) by capillary wet transfer overnight. Transferred RNA was crosslinked to the membrane with UV light at 1200 μJ (UV Stratalinker 2400, Stratagene; http://www.stratagene.com). The membrane was pre-hybridized in 5–10 mL of ULTRAhyb-Oligo Hybridization Buffer (http://www.ambion.com) at 40°C for at least 30 min. DNA oligo probes ( Additional file 2) were labeled using 32P with U4 polynucleotide kinase (Fermentas; http://www.fermentas.com). The reaction mix was incubated at 37°C for 1 h. Unincorporated nucleotides were separated using a Microspin G-25 column according to the manufacturer’s instructions (GE Healthcare, http://www.gelifesciences.com). The labeled probe was mixed with the hybridization solution and hybridized to the membrane at 40°C overnight. The membrane was washed with a washing solution (2xSSC, 0.1% SDS) at 40°C for 10 minutes, wrapped with plastic saran wrap, and exposed to phospho-imaging plates.
Quantitative real-time PCR (qPCR) and cDNA microarray
One μg of total RNA were reverse-transcribed using SuperScript® II reverse transcriptase following the manufacturer’s protocol (Invitrogen; http://www.invitrogen.com). Twenty μg of cDNA was used to perform quantitative real-time PCR with SYBR Green using gene-specific primers ( Additional file 4 and 10) designed for MIR and target genes. Elongation factor (NaEF) was used as reference house-keeping gene for analysis.
For reverse transcription of miRNA and tasiRNAs into cDNA, we used the miScript Reverse Transcription Kit (Qiagen; http://www.qiagen.com). 10 μg of total RNA were used for qPCR ( Additional file 5) with the miScript SYBR Green PCR Kit (Qiagen; http://www.qiagen.com) to quantify miRNAs and tasiRNAs. All qPCR data were analyzed using the 2ΔΔCt calculation method .
We used a cDNA microarray NCBI GEO database (Platform GPL13527, accession number GSE30287) . For data analyses, raw data was normalized to the 75th percentile and log2-transformed. Comparisons with greater than a 2-fold change were tested by Fisher’s PLSD test following an ANOVA.
JA and JA-Ile were co-extracted from leaf tissue as previously described . One hundred to 150 mg of lamina tissue from control and treated plants were used for phytohormone extraction with 1 ml of ethyl acetate spiked with 200 ng of D2-JA, and 40 ng of 13C6-JA-Ile as internal standards. Fifteen μL of the supernatant were analyzed on a Varian 1200 L Triple-Quadrapol-MS with a ProntoSIL column (C18; 5 μm, 50 × 2 mm).
Data were calculated with the StatView Software using the one-way analysis of variance ANOVA (means were compared by the lowest standard deviation (LSD)) algorithm (SAS Institute Inc., Cary, NC, USA).
Wu L, Zhou H, Zhang Q, Zhang J, Ni F, Liu C, Qi Y: DNA methylation mediated by a microRNA pathway. Mol Cell. 2010, 38 (3): 465-475. 10.1016/j.molcel.2010.03.008.
Ronemus M, Martienssen R: RNA interference: methylation mystery. Nature. 2005, 433 (7025): 472-473. 10.1038/433472a.
Chellappan P, Xia J, Zhou X, Gao S, Zhang X, Coutino G, Vazquez F, Zhang W, Jin H: siRNAs from miRNA sites mediate DNA methylation of target genes. Nucleic Acids Res. 2010, 38 (20): 6883-6894. 10.1093/nar/gkq590.
Khraiwesh B, Arif MA, Seumel GI, Ossowski S, Weigel D, Reski R, Frank W: Transcriptional control of gene expression by microRNAs. Cell. 2010, 140 (1): 111-122. 10.1016/j.cell.2009.12.023.
Bartel DP: MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004, 116 (2): 281-297. 10.1016/S0092-8674(04)00045-5.
Henderson IR, Zhang X, Lu C, Johnson L, Meyers BC, Green PJ, Jacobsen SE: Dissecting Arabidopsis thaliana DICER function in small RNA processing, gene silencing and DNA methylation patterning. Nat Genet. 2006, 38 (6): 721-725. 10.1038/ng1804.
Kurihara Y, Watanabe Y: Arabidopsis micro-RNA biogenesis through Dicer-like 1 protein functions. Proc Natl Acad Sci USA. 2004, 101 (34): 12753-12758. 10.1073/pnas.0403115101.
Ruiz-Ferrer V, Voinnet O: Roles of plant small RNAs in biotic stress responses. Annu Rev Plant Biol. 2009, 60: 485-510. 10.1146/annurev.arplant.043008.092111.
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.
Xie Z, Allen E, Wilken A, Carrington JC: DICER-LIKE 4 functions in trans-acting small interfering RNA biogenesis and vegetative phase change in Arabidopsis thaliana. Proc Natl Acad Sci USA. 2005, 102 (36): 12984-12989. 10.1073/pnas.0506426102.
Marin E, Jouannet V, Herz A, Lokerse AS, Weijers D, Vaucheret H, Nussaume L, Crespi MD, Maizel A: miR390, Arabidopsis TAS3 tasiRNAs, and their AUXIN RESPONSE FACTOR targets define an autoregulatory network quantitatively regulating lateral root growth. Plant Cell. 2010, 22 (4): 1104-1117. 10.1105/tpc.109.072553.
Chen XM: Small RNAs and their roles in plant development. Annu Rev Cell Dev Bi. 2009, 25: 21-44. 10.1146/annurev.cellbio.042308.113417.
Willmann MR, Mehalick AJ, Packer RL, Jenik PD: MicroRNAs regulate the timing of embryo maturation in Arabidopsis. Plant Physiol. 2011, 155 (4): 1871-1884. 10.1104/pp.110.171355.
Kidner CA, Martienssen RA: Spatially restricted microRNA directs leaf polarity through ARGONAUTE1. Nature. 2004, 428 (6978): 81-84. 10.1038/nature02366.
Palatnik JF, Allen E, Wu X, Schommer C, Schwab R, Carrington JC, Weigel D: Control of leaf morphogenesis by microRNAs. Nature. 2003, 425 (6955): 257-263. 10.1038/nature01958.
Yoon EK, Yang JH, Lim J, Kim SH, Kim SK, Lee WS: Auxin regulation of the microRNA390-dependent transacting small interfering RNA pathway in Arabidopsis lateral root development. Nucleic Acids Res. 2010, 38 (4): 1382-1391. 10.1093/nar/gkp1128.
Liu Q, Chen YQ: Insights into the mechanism of plant development: interactions of miRNAs pathway with phytohormone response. Biochem Biophys Res Commun. 2009, 384 (1): 1-5. 10.1016/j.bbrc.2009.04.028.
Schommer C, Palatnik JF, Aggarwal P, Chetelat A, Cubas P, Farmer EE, Nath U, Weigel D: Control of jasmonate biosynthesis and senescence by miR319 targets. PLoS Biol. 2008, 6 (9): e230-10.1371/journal.pbio.0060230.
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.
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.
Kulcheski FR, de Oliveira LF, Molina LG, Almerao MP, Rodrigues FA, Marcolino J, Barbosa JF, Stolf-Moreira R, Nepomuceno AL, Marcelino-Guimaraes FC, et al: Identification of novel soybean microRNAs involved in abiotic and biotic stresses. BMC Genomics. 2011, 12: 307-10.1186/1471-2164-12-307.
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.
Zhou L, Liu Y, Liu Z, Kong D, Duan M, Luo L: Genome-wide identification and analysis of drought-responsive microRNAs in Oryza sativa. J Exp Bot. 2010, 61 (15): 4157-4168. 10.1093/jxb/erq237.
Zhou X, Wang G, Zhang W: UV-B responsive microRNA genes in Arabidopsis thaliana. Mol Syst Biol. 2007, 3: 103-
Chiou TJ, Aung K, Lin SI, Wu CC, Chiang SF, Su CL: Regulation of phosphate homeostasis by MicroRNA in Arabidopsis. Plant Cell. 2006, 18 (2): 412-421. 10.1105/tpc.105.038943.
Reyes JL, Chua NH: ABA induction of miR159 controls transcript levels of two MYB factors during Arabidopsis seed germination. Plant J. 2007, 49 (4): 592-606. 10.1111/j.1365-313X.2006.02980.x.
Jones-Rhoades MW, Bartel DP: Computational identification of plant microRNAs and their targets, including a stress-induced miRNA. Mol Cell. 2004, 14 (6): 787-799. 10.1016/j.molcel.2004.05.027.
Tang S, Wang Y, Li Z, Gui Y, Xiao B, Xie J, Zhu QH, Fan L: Identification of wounding and topping responsive small RNAs in tobacco (Nicotiana tabacum). BMC Plant Biol. 2012, 12 (1): 28-10.1186/1471-2229-12-28.
Zhang W, Gao S, Zhou X, Chellappan P, Chen Z, Zhang X, Fromuth N, Coutino G, Coffey M, Jin H: Bacteria-responsive microRNAs regulate plant innate immunity by modulating plant hormone networks. Plant Mol Biol. 2011, 75 (1–2): 93-105.
Navarro L, Jay F, Nomura K, He SY, Voinnet O: Suppression of the microRNA pathway by bacterial effector proteins. Science. 2008, 321 (5891): 964-967. 10.1126/science.1159505.
Kasschau KD, Xie Z, Allen E, Llave C, Chapman EJ, Krizan KA, Carrington JC: P1/HC-Pro, a viral suppressor of RNA silencing, interferes with Arabidopsis development and miRNA function. Dev Cell. 2003, 4 (2): 205-217. 10.1016/S1534-5807(03)00025-X.
Li H, Deng Y, Wu T, Subramanian S, Yu O: Misexpression of miR482, miR1512, and miR1515 increases soybean nodulation. Plant Physiol. 2010, 153 (4): 1759-1770. 10.1104/pp.110.156950.
Schwachtje J, Baldwin IT: Why does herbivore attack reconfigure primary metabolism?. Plant Physiol. 2008, 146 (3): 845-851. 10.1104/pp.107.112490.
Kessler A, Baldwin IT: Herbivore-induced plant vaccination. Part I. The orchestration of plant defenses in nature and their fitness consequences in the wild tobacco Nicotiana attenuata. Plant J. 2004, 38 (4): 639-649. 10.1111/j.1365-313X.2004.02076.x.
Kang JH, Wang L, Giri A, Baldwin IT: Silencing threonine deaminase and JAR4 in Nicotiana attenuata impairs jasmonic acid-isoleucine-mediated defenses against Manduca sexta. Plant Cell. 2006, 18 (11): 3303-3320. 10.1105/tpc.106.041103.
Halitschke R, Schittko U, Pohnert G, Boland W, Baldwin IT: Molecular interactions between the specialist herbivore Manduca sexta (Lepidoptera, Sphingidae) and its natural host Nicotiana attenuata. III. Fatty acid-amino acid conjugates in herbivore oral secretions are necessary and sufficient for herbivore-specific plant responses. Plant Physiol. 2001, 125 (2): 711-717. 10.1104/pp.125.2.711.
Bonaventure G, VanDoorn A, Baldwin IT: Herbivore-associated elicitors: FAC signaling and metabolism. Trends Plant Sci. 2011, 16 (6): 294-299. 10.1016/j.tplants.2011.01.006.
McCloud ES, Baldwin IT: Herbivory and caterpillar regurgitants amplify the wound-induced increases in jasmonic acid but not nicotine in Nicotiana sylvestris. Planta. 1997, 203 (4): 430-435. 10.1007/s004250050210.
Halitschke R, Ziegler J, Keinanen M, Baldwin IT: Silencing of hydroperoxide lyase and allene oxide synthase reveals substrate and defense signaling crosstalk in Nicotiana attenuata. Plant J. 2004, 40 (1): 35-46. 10.1111/j.1365-313X.2004.02185.x.
Kallenbach M, Alagna F, Baldwin IT, Bonaventure G: Nicotiana attenuata SIPK, WIPK, NPR1, and fatty acid-amino acid conjugates participate in the induction of jasmonic acid biosynthesis by affecting early enzymatic steps in the pathway. Plant Physiol. 2010, 152 (1): 96-106. 10.1104/pp.109.149013.
Paschold A, Halitschke R, Baldwin IT: Co(i)-ordinating defenses: NaCOI1 mediates herbivore- induced resistance in Nicotiana attenuata and reveals the role of herbivore movement in avoiding defenses. Plant J. 2007, 51 (1): 79-91. 10.1111/j.1365-313X.2007.03119.x.
Pandey SP, Shahi P, Gase K, Baldwin IT: Herbivory-induced changes in the small-RNA transcriptome and phytohormone signaling in Nicotiana attenuata. Proc Natl Acad Sci USA. 2008, 105 (12): 4559-4564. 10.1073/pnas.0711363105.
Liu Y, Gao Q, Wu B, Ai T, Guo X: NgRDR1, an RNA-dependent RNA polymerase isolated from Nicotiana glutinosa, was involved in biotic and abiotic stresses. Plant Physiol Biochem. 2009, 47 (5): 359-368. 10.1016/j.plaphy.2008.12.017.
Pandey SP, Baldwin IT: RNA-directed RNA polymerase 1 (RdR1) mediates the resistance of Nicotiana attenuata to herbivore attack in nature. Plant J. 2007, 50 (1): 40-53. 10.1111/j.1365-313X.2007.03030.x.
Bozorov TA, Prakash Pandey S, Dinh ST, Kim SG, Heinrich M, Gase K, Baldwin IT: DICER-like proteins and their role in plant-herbivore interactions in Nicotiana attenuata. J Integr Plant Biol. 2012, 54 (3): 189-206. 10.1111/j.1744-7909.2012.01104.x.
Kim SG, Yon F, Gaquerel E, Gulati J, Baldwin IT: Tissue specific diurnal rhythms of metabolites and their regulation during herbivore attack in a native tobacco, Nicotiana attenuata. PLoS One. 2011, 6 (10): e26214-10.1371/journal.pone.0026214.
Kallenbach M, Bonaventure G, Gilardoni PA, Wissgott A, Baldwin IT: Empoasca leafhoppers attack wild tobacco plants in a jasmonate-dependent manner and identify jasmonate mutants in natural populations. Proc Natl Acad Sci USA. 2012, 109 (24): 1548-1557. 10.1073/pnas.1200363109.
Liang C, Zhang X, Zou J, Xu D, Su F, Ye N: Identification of miRNA from Porphyra yezoensis by high-throughput sequencing and bioinformatics analysis. PLoS One. 2010, 5 (5): e10698-10.1371/journal.pone.0010698.
Frazier TP, Xie F, Freistaedter A, Burklew CE, Zhang B: Identification and characterization of microRNAs and their target genes in tobacco (Nicotiana tabacum). Planta. 2010, 232 (6): 1289-1308. 10.1007/s00425-010-1255-1.
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.
Shen D, Wang S, Chen H, Zhu QH, Helliwell C, Fan LJ: Molecular phylogeny of miR390-guided trans-acting siRNA genes (TAS3) in the grass family. Plant Syst Evol. 2009, 283 (1–2): 125-132.
Adenot X, Elmayan T, Lauressergues D, Boutet S, Bouche N, Gasciolli V, Vaucheret H: DRB4-dependent TAS3 trans-acting siRNAs control leaf morphology through AGO7. Curr Biol. 2006, 16 (9): 927-932. 10.1016/j.cub.2006.03.035.
Kim JH, Woo HR, Kim J, Lim PO, Lee IC, Choi SH, Hwang D, Nam HG: Trifurcate feed-forward regulation of age-dependent cell death involving miR164 in Arabidopsis. Science. 2009, 323 (5917): 1053-1057. 10.1126/science.1166386.
Schwachtje J, Minchin PE, Jahnke S, van Dongen JT, Schittko U, Baldwin IT: SNF1-related kinases allow plants to tolerate herbivory by allocating carbon to roots. Proc Natl Acad Sci USA. 2006, 103 (34): 12935-12940. 10.1073/pnas.0602316103.
Halitschke R, Baldwin IT: Antisense LOX expression increases herbivore performance by decreasing defense responses and inhibiting growth-related transcriptional reorganization in Nicotiana attenuata. Plant J. 2003, 36 (6): 794-807. 10.1046/j.1365-313X.2003.01921.x.
Luo QJ, Mittal A, Jia F, Rock CD: An autoregulatory feedback loop involving PAP1 and TAS4 in response to sugars in Arabidopsis. Plant Mol Biol. 2012, 80 (1): 117-129. 10.1007/s11103-011-9778-9.
Kaur H, Heinzel N, Schottner M, Baldwin IT, Galis I: R2R3-NaMYB8 regulates the accumulation of phenylpropanoid-polyamine conjugates, which are essential for local and systemic defense against insect herbivores in Nicotiana attenuata. Plant Physiol. 2010, 152 (3): 1731-1747. 10.1104/pp.109.151738.
Onkokesung N, Gaquerel E, Kotkar H, Kaur H, Baldwin IT, Galis I: MYB8 controls inducible phenolamide levels by activating three novel hydroxycinnamoyl-coenzyme A:polyamine transferases in Nicotiana attenuata. Plant Physiol. 2012, 158 (1): 389-407. 10.1104/pp.111.187229.
Schwab R, Voinnet O: miRNA processing turned upside down. EMBO J. 2009, 28 (23): 3633-3634. 10.1038/emboj.2009.334.
Steffen P, Voss B, Rehmsmeier M, Reeder J, Giegerich R: RNAshapes: an integrated RNA analysis package based on abstract shapes. Bioinformatics. 2006, 22 (4): 500-503. 10.1093/bioinformatics/btk010.
Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Mol Biol Evol. 2007, 24 (8): 1596-1599. 10.1093/molbev/msm092.
Molnar A, Schwach F, Studholme DJ, Thuenemann EC, Baulcombe DC: miRNAs control gene expression in the single-cell alga Chlamydomonas reinhardtii. Nature. 2007, 447 (7148): 1126-1129. 10.1038/nature05903.
Bubner B, Gase K, Baldwin IT: Two-fold differences are the detection limit for determining transgene copy numbers in plants by real-time PCR. BMC Biotechnol. 2004, 4: 14-10.1186/1472-6750-4-14.
Wang L, Halitschke R, Kang JH, Berg A, Harnisch F, Baldwin IT: Independently silencing two JAR family members impairs levels of trypsin proteinase inhibitors but not nicotine. Planta. 2007, 226 (1): 159-167. 10.1007/s00425-007-0477-3.
We thank Dr. Tamara Krügel and the gardening staff for growing and maintaining plants, Dr. Meredith C. Schuman for language editing, and the Max Planck Society for funding.
The authors declare that they have no competing interests.
TAB performed the experimental work and analyzed the data. TAB, ITB and SGK participated in the design of the study. ITB and SGK conceived of the study and edited the manuscript. TAB drafted the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 8: Transcript abundance of miRNA target genes. Abundance of miRNA targets after W+W and W+OS treatments. For each sample, one leaf on the rosette of a 32-day-old WT plant was left untreated (control) or treated with wounding plus water (W+W) or wounding plus OS (W+OS) and harvested 1 h or 5 h post treatment. Shown are mean (± SE) levels of three replicates per line. Letters indicate significant differences (P < 0.05) in Fisher’s PLSD tests following an ANOVA. (TIFF 755 KB)
Additional file 9: JA and JA-Ile levels were impaired in ir-aoc plants after elicitation. Accumulation of Jasmonic acid (JA) and jasmonoyl isoleucine (JA-Ile) increased in W+W and W+OS treated leaves and these accumulations were dramatically altered in ir-aoc plants. For each sample, one leaf on the rosette of a 32-day-old plant was treated and harvested 1 h or 5 h post treatment. Untreated plants served as control. Asterisks indicate significant differences (***, P < 0.001) in Fisher’ s PLSD tests following an ANOVA. (TIFF 900 KB)
Additional file 10: qPCR data showing accumulation of primary miRNA transcripts in W+W and W+OS treated leaves in wild type and ir- aoc. For each sample, one leaf on the rosette of a 32-day-old plants was treated with wounding plus water (W+W) or wounding plus OS (W+OS) and harvested 1 h or 5 h post treatment. Untreated plants served as control (CTRL). Shown are mean (± SE) levels of three replicates. Letters indicate significant differences (P < 0.05) in Fisher’s PLSD test following an ANOVA. (TIFF 1 MB)
Additional file 11: Time course expression of miRNAs and their targets in W+W and W+OS treated leaves. (A) qPCR data showing the abundance of mature miRNAs. (B) Microarray data showing the accumulation of miRNAs targets. For each sample, one leaf on the rosette of a 32-day-old WT plant was left untreated (control) or treated with wounding plus water (W+W) or wounding plus OS (W+OS) and harvested 1 h or 5 h post treatment. Shown are mean (± SE) levels of three replicates per line. Lowercase letters (W+OS) and italic letters (W+W) indicate significant differences (P < 0.05) in Fisher’s PLSD tests following an ANOVA. (TIFF 6 MB)
Additional file 12: Microarray data showing accumulation of primary miRNAs in W+W and W+OS treated leaves of N. attenuata . Mean (± SE) levels of normalized transcript abundance of miRNAs and tasiRNAs target genes after W+W and W+OS treatments. Untreated plants served as control (CTRL). Signal intensities in microarray data was normalized using the 75th percentile value and log2-transformation. Letters indicate significant differences (P < 0.05) in Fisher’s PLSD test following an ANOVA. (TIFF 1 MB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.