Skip to main content

Molecular mechanism of modulating miR482b level in tomato with botrytis cinerea infection

Abstract

Background

Plant miRNAs are involved in the response to biotic and abiotic stresses by altering their expression levels, and they play an important role in the regulation of plant resistance to stress. However, the molecular mechanism that regulates the expression levels of miRNAs in plants with biotic and abiotic stress still needs to be explored. Previously, we found that the expression of the miR482 family was changed in tomato infected by Botrytis cinerea. In this study, we investigated and uncovered the mechanism underlying the response of miR482 to B. cinerea infection in tomato.

Results

First, RT-qPCR was employed to detect the expression patterns of miR482b in tomato infected by B. cinerea, and results showed that miR482b primary transcripts (pri-miR482b) were up-regulated in B. cinerea-infected leaves, but the mature miR482b was down-regulated. Subsequently, we used rapid amplification cDNA end method to amplify the full-length of pri-miR482b. Result showed that the pri-miR482b had two isoforms, with the longer one (consisting 300 bp) having an extra fragment of 53 bp in the 3’-end compared with the shorter one. In vitro Dicer assay indicated that the longer isoform pri-miR482b-x1 had higher efficiency in the post-transcriptional splicing of miRNA than the shorter isoform pri-miR482b-x2. In addition, the transcription level of mature miR482b was much higher in transgenic Arabidopsis overexpressing pri-miR482b-x1 than that in OE pri-miR482b-x2 Arabidopsis. These results confirmed that this extra 53 bp in pri-miR482b-x1 might play a key role in the miR482b biogenesis of post-transcription processing.

Conclusions

Extra 53 bp in pri-miR482b-x1 enhanced miR482b biogenesis, which elevated the transcription level of miR482b. This study clarified the response of miR482 to B. cinerea infection in tomato, thereby helping us further understand the molecular mechanisms that regulate the expression levels of other miRNAs.

Peer Review reports

Background

MicroRNAs (miRNAs) are a class of sRNAs with length of 20 ~ 24 nt. In plants, the miRNA gene is transcribed by RNA polymerase II [1]. The primary transcripts (pri-miRNAs) fold into an incomplete stem-loop structure, which is cleaved by RNase III type enzyme Dicer-like 1 (DCL1) to form an incomplete hairpin precursor miRNA, termed pre-miRNAs. The pre-miRNAs are then cleaved by DCL1 or DCL4 [2, 3] to generate the miRNA duplex, namely miRNA: miRNA* double-stranded dimer. Methylation of the 3’-end of the duplex is by methyltransferase HEN1. The plant exportin-5 homologous protein (HASTY, HST) transports the miRNA duplex into the cytoplasm, and they are incorporated into the RISC (RNA-induced silencing complex) together with an Argonaute (AGO) protein, where one strand is selected to become the mature miRNA [4], binding to the target mRNA complements to cleave or inhibit translation to achieve negative regulation of the target gene. In addition to the above classical mode of miRNA synthesis, it can be synthesized by the action of the DCL enzyme or the like [5].

MiRNAs act as negative regulators of gene expression in eukaryotes and participate to regulate the growth and development of plants [6,7,8,9,10] and disease resistance [11,12,13,14]. To date, many miRNAs exhibit different expression profiles in plants responding to B. cinerea infection. Zhao et al. found that miR5254, miR165a-3p, miR3897-3p and miR6450a are involved in the defense response of tree peony during the invasion of B. cinerea [15]. Soto-Suarez et al. found that in addition to its role in controlling development, miR396 contributes to the dynamic defense response to necrotrophic (B. cinerea) and hemibiotrophic (P. cucumerina) fungal pathogens in Arabidopsis [16]. In strawberry, miR5290a negatively regulates its target gene PIRL to increase resistance to B. cinerea [17]. The 31 miRNAs in lily are differentially expressed in leaves infected by Botrytis ellipsoidea and respond to the stress of B. ellipsoidea [18]. In our previous study, we found that the expression levels of multiple miRNAs in tomato are correlated with B. cinerea infection via miRNA microarray and high-throughput sequencing technology [19, 20]. Moreover, miR319 and miR394 act as positive and negative regulators in resistance to B. cinerea infection, respectively [21, 22].

Although the expression levels and the biological functions of many miRNAs have been well studied in plants in response to environmental changes, the mechanisms that regulate the expression of the environmentally responsive miRNAs are still poorly understood. Bielewicz et al. found that introns are crucial for the expression levels of two miRNA genes: MIR163 and MIR161. Removal of their introns leads to a drop-off in the level of both miRNAs [23]. In addition, Schwab et al. showed that the introns located in the 3’ -end of the stem-loop structure can promote mature miRNA accumulation [24]. Inhibition of splicing at the 93 − 25 3’ SS (between the sequences of pre-miRNA 93 and 25) results in elevated miRNA levels [25]. In the previous study, we found that the expression level of miR482 was changed in tomato leaves with B. cinerea infection [26]. In the present study, we found that miR482b was significantly down-regulated in tomato leaves infected by B. cinerea, but its primary transcript was up-regulated. This result suggested that the different expression patterns between miR482b and its primary might be involved in the processing efficiency of post-transcriptional levels rather than that in transcription levels. Therefore, in this study, we aimed to reveal the mechanism that regulated the expression of miR482b in tomato’s response to B. cinerea infection.

Results

1. Expression patterns of miR482b and its primary in B. cinerea-infected tomato leaves

To understand the expression patterns of miR482b, the abundance of the mature miRNA and its primary (pri-miR482b) was measured by quantitative reverse transcription PCR (RT-qPCR) in B. cinerea-infected tomato at different time points. The results showed that the abundance of pri-miR482b was significantly up-regulated at 24 and 72 hpi (Fig. 1 A). Interestingly, the expression level of miR482b was significantly down-regulated at all three time points (Fig. 1B). Correspondingly, four nucleotide binding site-leucine-rich repeat genes (NBS-LRR), namely Solyc02g036270.2, Solyc04g009070.1, Solyc12g016220.2, and Solyc05g008070.2, which were experimentally confirmed as the target genes of miR482b [27], were up-regulated in B. cinerea-infected leaves at all three time points compared with 0 hpi (Fig. 1 C-F), showing negative regulation by miR482b. These results proposed that miR482b was induced at the transcriptional level but inhibited the splicing process after transcription in tomato leaves infected by B. cinerea.

Fig. 1
figure 1

Expression patterns of miR482b and its target genes in mock- and B. cinerea-infected leaves at 0, 12, 24 and 72 hpi. A-B) The expressions of pri-miR482b (A) and miR482b (B). SlyU6 were used as the internal control. Results are expressed as means ± SD of three biological replicates. Asterisks indicate a significant difference (P < 0.05) compared with the corresponding 0 hpi leaves. C-F) The expressions of four NBS-LRR genes (Solyc02g036270.2, Solyc04g009070.1, Solyc12g016220.2, and Solyc05g008070.2) targeted by miR482b. SlyUbq3 were used as the internal control. Results are expressed as means ± SD of three biological replicates. Asterisks indicate a significant difference (P < 0.05) compared with the corresponding 0 hpi leaves

2. Cloning and chromosomal location of pri-miR482b

To detect whether the down-regulation of miR482b is due to post-transcriptional inhibition, the full-length cDNA of the miR482b primary was first cloned by the RACE method. The results showed that a 256 bp fragment was amplified by 5’-RACE (Fig. 2 A). Interestingly, 3’-RACE results showed two different amplicons in 3’-RACE (Fig. 2 A). DNA sequencing revealed that the longer one was 254 bp in length, and the shorter one was 201 bp, with a 53 bp deletion at the 3’-end of the longer one. The complete cDNA of the miR482b gene was compiled by overlapping the sequences of the cloned cDNA and the 5’-RACE and 3’-RACE PCR products. Two isoforms of pri-miR482b transcripts without poly(A) consisted of 300 and 247 bp, respectively. The longer isoform was named pri-miR482b-x1, and the shorter one was named pri-miR482b-x2 (Fig. 2B C). Compared with pri-miR482b-x2, an extra small hairpin structure was presented in the 3’-end of pri-miR482b-x1 (Fig. 2D). On the basis of the full-length cDNA sequence of pri-miR482b, miR482b was located in the intergenic region, and the transcription start site was at 37,497,067 of the minus strand of the tomato chromosome 6 which had a full length of 49,794,276 bp. Both sequences have been deposited into the GenBank database (accession numbers: MW590251 and MW590252).

Fig. 2
figure 2

Cloning and sequence analysis of full-length pri-miR482b. A) RACE Amplification of pri-miR482b. Lane M: DL 2000 DNA marker; Lane 1: 5’ end of pri-miR482b; Lane2: 3’ end of pri-miR482b. B) Schematic of two isoforms of pri-miR482b. C) Alignment of two isoforms of pri-miR482b. D) Secondary structures of two pri-miR482b isoforms

3. Expression patterns of pri-miR482b-x1 and pri-miR482b-x2 in B. cinerea-infected tomato leaves

To understand the expression patterns of both isoforms of pri-miR482b, the abundance of pri-miR482b-x1 and pri-miR482b-x2 was further quantified in B. cinerea-infected tomato at different time points. The transcript level of pri-miR482b-x1 was significantly down-regulated in B. cinerea-inoculated leaves at 12, 24 and 72 hpi compared with that at 0 hpi, whereas the transcript level of pri-miR482b-x2 was significantly up-regulated at 24 and 72 hpi (Fig. 3). These results indicated that miR482b was down-regulated in B. cinerea-inoculated leaves due to the inhibited expression of pri-miR482b-x1. Thus, the extra stem-loop structure in pri-miR482b-x1 might play a key role in the miR482b biogenesis of post-transcription processing.

Fig. 3
figure 3

Expression patterns of two pri-miR482b isoforms in mock- and B. cinerea-infected leaves at 0, 12, 24 and 72 hpi. SlyUbq3 was used as the internal control. Results are expressed as means ± SD of three biological replicates. Various letters indicate a significant difference among samples (P < 0.05)

4. Detection of the splicing efficiency for the two pri-miR482b isoforms

To detect whether the two isoforms of pri-miR482b have different efficiencies in the biogenesis of miR482b, both biotin-labeled pri-miR482b-x1 and -x2 were transcribed in vitro via T7 RNA polymerase with NTPs and biotin-labeled UTP and then incubated with 2 µL of miRNA splicing proteins (5 mg/mL). The results showed that pri-miR482b-x1 was spliced to produce pre-miR482b-x1 and mature miR482b, whereas pri-miR482b-x2 only produced a pre-miR482b-like fragment; no mature miR482b was produced (Fig. 4 A; Additional file 1 (Figure S1B-C)).

Fig. 4
figure 4

Dicer splicing activities of two pri-miR482b isoforms. (A) Splicing efficiencies of two pri-miR482b isoforms in vitro. Land 1: pri-miR482b-x1; Lane 2: pri-miR482b-x2; Lane 3: pri-miR482b-x1 with 2 mg/mL protein extraction; Lane 4: pri-miR482b-x2 with 2 mg/mL protein extraction. (B) Expression levels of pri-miR482b-x1 in pri-miR482b-x1 OE plants and the expression level of pri-miR482b-x2 in pri-miR482b-x2 OE plants via semi-quantification RT-PCR. Lanes 1–3: three OE lines of pri-miR482b-x1; Lanes 4–6: three OE lines of pri-miR482b-x2. C-D) Expression levels of pre-miR482b (C) and miR482b (D) in pri-miR482b-x1 and pri-miR482b-x2 OE plants via quantification RT-PCR. AtUBQ10 was used as the internal control. Results are expressed as means ± SD of three biological replicates. Asterisks indicates a significant difference (P < 0.05) compared with the corresponding OE_pri-miR482b-x1-1

To investigate the splicing efficiency of both isoforms in vivo, Arabidopsis, which does not have an endogenous miR482 gene was selected for overexpressing pri-miR482b-x1 and pri-miR482b-x2 by transgenic technology. RT-PCR results showed that pri-miR482b-x1 and -x2 were expressed in transgenic plants overexpressing pri-miR482b-x1 and -x2, respectively. The expression levels of pri-miR482b were lower in pri-miR482b-x1 overexpression (OE) plants than in pri-miR482b-x2 OE plant (Fig. 4B). Moreover, RT-qPCR results revealed that the expression level of pre-miR482b was lower in pri-miR482b-x1 OE plants than in pri-miR482b-x2 OE plants (Fig. 4 C), but mature miR482b had higher abundance in pri-miR482b-x1 OE plants than in pri-miR482b-x2 OE plants (Fig. 4D). These results confirmed that the 53 nt extra sequence in the 3’-end of pri-miR482b-x1 played a critical role in miR482b biogenesis after transcription.

5. Function of miR482b in plants infected by B. cinerea

To understand the role of miR482b in plants infected by B. cinerea, WT and transgenic plants were inoculated with B. cinerea for 48 h, and the physical appearance of the plants was assessed. The transgenic plants had larger necrotic spots than the WT plants (Fig. 5 A-C; Additional file 1 (Figure S2)). Moreover, the necrotic spots of pri-miR482b-x1 OE plants were ~ 2 fold larger than those of pri-miR482b-x2 OE plants. In addition, two marker genes, plant defensing 1.2 (PDF1.2) and pathogen related protein 4 (PR4) [28], of the ERF branch of the jasmonic acid signaling pathway were quantitative detected in WT and transgenic plants. The results showed that the expression of both genes was downregulated in transgenic pri-miR482b OE plants compared with wide type plants. Moreover, the inhibition of both genes is more significant in pri-miR482b-x1 OE plants than that in pri-miR482b-x2 OE plants (Fig. 6).

Fig. 5
figure 5

Resistance analysis of transgenic Arabidopsis overexpressing pri-miR482b-x1 and pri-miR482b-x2 against B. cinerea. (A) Disease symptoms on the B. cinerea-infected leaves of WT and transgenic Arabidopsis overexpressing pri-miR482b-x1and pri-miR482b-x2 for 48 h. (B) Trypan blue staining confirmed the pathological spots. (C) Statistical analysis of the pathogenic spot diameters. Asterisks indicate a significant difference (P < 0.01) compared with the corresponding WT

Fig. 6
figure 6

Expression levels of two genes (PDF1.2 and PR4) of JA-signal pathway in WT and transgenic Arabidopsis. SlyUbq3 was used as the internal control. Results are expressed as means ± SD of three biological replicates. Various letters indicate a significant difference among samples (P < 0.05)

Discussion

Given the important role of miRNAs in the regulation of plant growth and development and stress response, the abundance of miRNAs is tightly regulated at multiple levels, including transcriptional and post transcriptional steps [29]. In transcriptional steps, evidence suggested that general RNA-binding proteins (RBPs) binding to the terminal loop or stem of an miRNA stem-loop structure can positively or negatively affect microprocessor mediated pri-miRNA and/or pre-miRNA processing [30]. For example, the RBP protein LIN28 can specifically recognize and bind to the terminal loop of let-7 and then affect the processing of the let-7 precursor by blocking the activity of Drosha and Dicer [31,32,33]. HnRNP A1 binds to the terminal loop of pri-mir-18a and induces relaxation at the stem–loop structure near the DROSHA cleavage site, resulting in increased efficiency of miRNA processing [29]. In addition, the sequence variation, such as single nucleotide polymorphisms, may play a role in the biogenesis of miRNA. A rare genetic variation in the terminal loop of pri-miR-30c-1 (G27 to A), which directly affects the processing of pri-mir-30c-1 by inducing a secondary RNA structure rearrangement and facilitates binding of the trans-acting factor SRSF3 [34], results in increased levels of mature miR-30c [34, 35]. Finally, the introns are crucial for the expression levels of two miRNA genes, namely, MIR163 and MIR161. Removal of their introns leads to a drop-off in the level of both miRNAs [23]. Schwab et al. also showed that the introns located at the 3’-end of the stem-loop can promote mature miRNA accumulation [24]. In this study, we found that the processing efficiency of pri-miR482b was inhibited when a 53 nt fragment was absent in the 3’-end of pri-miR482b.

The miR482 family comprises plant-specific small RNAs that have been found in 23 plants including Populus trichocarpa [36], Pinus taeda [37, 38], Glycinemax [39], Malus domestica [40], Phaseolus vulgaris [40], Medicago truncatula [41], Solanum lycopersicum [42], and so on. Evidence showed that miR482 is involved in response to biotic stress and acts as a negative regulator in plant-pathogen interaction through inhibited R genes. In tomato and potato, members of the miR482 family are expected to target approximately 20 % of resistance genes [43]. Lu et al. [36, 37] showed that Ptc-miR482 can cleave anti-pathogenic protein genes involved in plant resistance to biotic and abiotic stresses. Therefore, overexpression of miR482 can decrease the resistance of host plants against pathogens. Yang et al. [44] found that the overexpression of potato miR482e enhances the sensitivity of plants to Verticillium dahliae by targeting a series of NBS-LRR genes. Jiang et al. [27] found that miR482b negatively regulates the infection of Phytophthora infestans by targeting the NBS-LRR resistance genes. Feng et al. [45] found that ghr-miR482a plays a role in resistance to G. barbadense by targeting the R gene. By contrast, to increase resistance, host plants would like to reduce the expression of miR482 and increase the expression of the NBS-LRR genes in the pathogen-host interaction. Zhu et al. [46] found that the expression levels of the NBS-LRR genes were induced in V. dahlia-infected cottons by inhibiting the expression of ghr-miR482c, ghr-miR482d.2 and ghr-miR482b/miR482b.2. Ouyang et al. [47] found that tomato plants up-regulated an R gene by inhibiting the expression of sly-miR482f to enhance resistance to Fusarium oxysporum. In this study, we also found that the expression levels of miR482b and its targets were down- and up-regulated in tomato with B. cinerea infection, respectively, resulting in the up-regulation of its target genes (NBS-LRRs). Moreover, we confirmed that increasing the transcript level of miR482b could decrease resistance to B. cinerea infection in plants. This result indicated that miR482b played a negative role in plant resistance to B. cinerea and suggested that the plants can resist the infection of B. cinerea by changing the isoform levels of miR482b.

The jasmonic acid-mediated defense pathway has an important role in plant resistance against necrotrophic pathogens [48]. There are two distinct branches in Arabidopsis JA pathway including ethylene response factor (ERF) branch and the MYC branch [28]. The ERF branch of the JA pathway is typically activated upon infection by necrotrophic pathogens [49]. Induction of the ERF branch results in the activation of a large set of JA/ET-responsive genes, including the marker gene PDF1.2 [50, 51]. In the other branch of the JA pathway, MYC branch, is typically activated upon wounding or feeding by herbivorous insects. Therefore, we tried to understand the effect of miR482b over-expression on the JA pathway through investigating the expression of two marker genes in ERF branch of JA pathway. Results showed that the PDF1.2, a marker genes of JA pathway representative ERF branch, was significantly decreased in pri-miR482b OE plants compared to WT. Moreover, PDF1.2 has lower expression level in pri-miR482b-x1 OE plants than that in pri-miR482b-x2 OE plants (Fig. 6). On the contrary, PDF1.2 would be induced by B. cinerea to activate JA pathway for increasing the resistance in Arabidopsis [28]. In addition, a similar expression pattern was observed for PR4, which is another JA-responsive gene of the ERF branch and its expression levels could be induced by necrotrophic pathogen [52], but was also decreased in pri-miR482b OE plants compared to WT. Our results suggested that the overexpression of pri-miR482b enhanced the susceptibility of the transgenic plants to B. cinerea infection through inhibiting the ERF branch of JA signaling pathway.

Conclusions

This study was the first to reveal the molecular mechanism underlying the suppression of miR482b expression in tomato infected by B. cinerea as follows: a 53 nt fragment was spliced in the 3’-end of the normal pri-miR482b, namely pri-miR482b-x1, in B. cinerea-infected tomato. The truncated pri-miR482b isoform, namely, pri-miR482b-x2, has lower processing efficiency in miR482b biogenesis, leading to the decrease in transcript level of miR482b in B. cinerea-infected tomato.

Methods

1. Plants and B. cinerea inoculation

Seeds of tomato cv. MicroTom, purchased from Nanjing Fengshuo Yuanyi Co., Ltd in China. Arabidopsis thaliana (Col-0) and Botrytis cinerea were provided by plant pathology laboratory of Zhejiang Sci-tech University (Hangzhou, Zhejiang). Tomato and B. cinerea were cultured according to Meng et al. [53]. Seeds of tomato were grown with a 12:12 h photoperiod at ~ 22 °C. Six-week-old plants were inoculated with B. cinerea solution containing 5 × 106 conidiospores/mL. The B. cinerea- and mock-inoculated leaves were harvested at 0, 12, 24 and 72 h post inoculation (hpi). The samples were frozen in liquid nitrogen and stored at -70 °C for transcript level analyses.

2. RNA extraction, reverse transcription and real-time PCR (RT-qPCR)

Total RNAs were extracted and quantified of total RNA were performed according to Meng et al. [53]. Total RNAs were extracted using TRIzol reagent, treated by RNase-free DNase, and quantified by using a NanoDrop ND-1000 spectrophotometer. For poly(A) RNAs, equal quantities of total RNA (1 µg) were reverse-transcribed at 42 °C using SuperScript III Reverse Transcriptase (Invitrogen) and 2.5 µM Oligo(dT18). A similar reaction without reverse transcriptase was also performed as a control to confirm the absence of genomic DNA in subsequent steps. For miRNAs, reverse transcription was performed using the One Step PrimeScript miRNA cDNA Synthesis Kit (TaKaRa, Dalian, China) according to the manufacturer’s protocol.

SYBR Green PCR was performed according to Meng et al. [53]. In brief, 2 µL of cDNA template was added to 12.5 µL of 2× SYBR Green PCR master mix (Takara), 1 µM specific primers and ddH2O to a final volume of 25 µL. The reactions were amplified for 10 s at 95 °C, followed by 40 cycles of 95 °C for 10 s and 60 °C for 30 s. All reactions were performed in triplicate, and controls (no template and no RT) were included for each gene. The threshold cycle (Ct) values were automatically determined by the ABI Prism 7300 Sequence Detection System (PE Applied Biosystems, USA). The fold-changes were calculated using the 2−ΔΔCT method, where ΔΔCt = (Ct,target - Ct,inner)Infection - (Ct,target − Ct,inner)Mock [54].

3. Rapid amplification of cDNA ends (RACE)

The full-length of pri-miR482b was obtained by rapid amplification of cDNA ends (RACE) with a SMART RACE cDNA Amplification Kit (Invitrogen) according to the manufacturer’s protocol. Total RNA (1 µg) obtained from tomato leaves was used for cDNA synthesis. The gene-specific primers pri-miR482b-R and pri-miR482b-F were used in the RACE of the 5’ and 3’ ends, respectively (Additional file 2). The 5’ and 3’ cDNA fragments obtained from RACE were cloned into a pMD19-T vector (Takara) and sequenced. The RNA secondary structure was predicted by the RNAfold program [55].

4. In vitro transcription and Dicer splicing assay of pri-miR482b

The transcription and splicing assay for the pri-miR482b in vitro were performed according to Qi et al. [56]. The DNA template of pri-miR482b was amplified by using T7 promoter anchored primers (Additional file 2). Resulting DNAs were used for in vitro transcription under the presence of Biotin-UTP according to the manufacturer’s protocol (Roche). For Dicer activity assay, RNAs were incubated with 10 µg of protein in 20 µL of reaction buffer containing 100 mM NaCl, 1 mM ATP, 0.2 mM GTP, 1.2 mM MgCl2, 25 mM creatine phosphate, 30 µg/mL creatine kinase, and 4 U RNase Inhibitor at room temperature for 10 min. RNAs were extracted, precipitated, and dissolved in water.

The RNAs were electrophoresed in 8 M Urea 10 % PAGE gel and then transferred to a Hybond N + membrane (Amersham). The membrane was cross-linked twice with 254 nm UV light at 120 mJ/cm2. Biotin signal was detected by a Chemiluminescent Biotin-labeled Nucleic Acid Detection Kit (Beyotime, China) according to the manufacturer’s protocol. Images of the membranes were captured with a chemiluminescence detection system (Chemi-Doc; Bio-Rad Laboratories, Ltd., Hemel Hempstead, UK) with a CCD camera after 10 s of exposure.

5. Gene constructs and the generation of transgenic Arabidopsis plants

Two isoforms of pri-miR482b were amplified from the cDNA of tomato leaves using specific primers (Additional file 2) and then cloned into the pBIN438 expression vector downstream of the CaMV 35 S promoter region through the Pst I and Xba I restriction sites. The construct was introduced into A. tumefaciens GV3101 and transformed into Arabidopsis “Columbia” (Col-0) according to Wu et al. [53].

Methods and plant material described above complied with relevant institutional, national and international guidelines and legislation.

Availability of data and materials

All data generated during this study are included in this published article and its supplementary information files, except for the two isoform sequences of pri-miR482b which had been deposited into GenBank database (accession numbers: MW590251 and MW590252).

Abbreviations

miRNAs:

microRNAs

DCL1:

Dicer-like 1

NBS-LRR:

Nucleotide binding site-leucine-rich repeats

RISC:

RNA-induced silencing complex

RT-qPCR:

quantitative reverse transcription PCR

WT:

wild type

OE:

overexpression

RBPs:

RNA-binding proteins

CaMV:

Cauliflower mosaic virus

RACE:

Rapid amplification of cDNA ends

PDF1.2:

plant defensing 1.2

PR4:

pathogen related protein 4

References

  1. Sun X, Lin L, Sui N. Regulation mechanism of microRNA in plant response to abiotic stress and breeding. Mol Biol Rep. 2019;46(1):1447–57.

    Article  CAS  PubMed  Google Scholar 

  2. Katiyar-Agarwal S, Gao S, Vivian-Smith A, Jin H. A novel class of bacteria-induced small RNAs in Arabidopsis. Genes Dev. 2007;21(23):3123–34.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Jin HL. Endogenous small RNAs and antibacterial immunity in plants. FEBS Lett. 2008;582(18):2679–84.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Budak H, Akpinar BA. Plant miRNAs: biogenesis, organization and origins. Funct Integr Genomics. 2015;15(5):523–31.

    Article  CAS  PubMed  Google Scholar 

  5. Kim YK, Kim B, Kim VN. Re-evaluation of the roles of DROSHA, Exportin 5, and DICER in microRNA biogenesis. Proc Natl Acad Sci U S A. 2016;113(13):E1881–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Duan P, Ni S, Wang J, Zhang B, Xu R, Wang Y, et al. Regulation of OsGRF4 by OsmiR396 controls grain size and yield in rice. Nat Plants. 2015;2:15203.

    Article  PubMed  CAS  Google Scholar 

  7. Gao F, Wang K, Liu Y, Chen Y, Chen P, Shi Z, et al. Blocking miR396 increases rice yield by shaping inflorescence architecture. Nat Plants. 2015;2:15196.

    Article  PubMed  CAS  Google Scholar 

  8. Zhang YC, Yu Y, Wang CY, Li ZY, Liu Q, Xu J, et al. Overexpression of microRNA OsmiR397 improves rice yield by increasing grain size and promoting panicle branching. Nat Biotechnol. 2013;31(9):848–52.

    Article  CAS  PubMed  Google Scholar 

  9. Wu G, Park MY, Conway SR, Wang JW, Weigel D, Poethig RS. The sequential action of miR156 and miR172 regulates developmental timing in Arabidopsis. Cell. 2009;138(4):750–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Li N, Yang TX, Guo ZY, Wang QS, Chai M, Wu MB, et al. Maize microRNA166 Inactivation Confers Plant Development and Abiotic Stress Resistance. Int J Mol Sci. 2020;21(24).

  11. Hong YH, Meng J, He XL, Zhang YY, Luan YS. Overexpression of MiR482c in Tomato Induces Enhanced Susceptibility to Late Blight. Cells-Basel. 2019;8(8).

  12. Jiang N, Meng J, Cui J, Sun GX, Luan YS. Function identification of miR482b, a negative regulator during tomato resistance to Phytophthora infestans. Hortic Res-England. 2018;5.

  13. Yang L, Mu XY, Liu C, Cai JH, Shi K, Zhu WJ, et al. Overexpression of potato miR482e enhanced plant sensitivity to Verticillium dahliae infection. J Integr Plant Biol. 2015;57(12):1078–88.

    Article  CAS  PubMed  Google Scholar 

  14. Zhang YJ, Wang W, Chen J, Liu JB, Xia MX, Shen FF. Identification of miRNAs and Their Targets in Cotton Inoculated with Verticillium dahliae by High- Throughput Sequencing and Degradome Analysis. Int J Mol Sci. 2015;16(7):14749–68.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Zhao DQ, Gong SJ, Hao ZJ, Tao J. Identification of miRNAs Responsive to Botrytis cinerea in Herbaceous Peony (Paeonia lactiflora Pall.) by High-Throughput Sequencing. Genes. 2015;6(3):918–34.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Soto-Suarez M, Baldrich P, Weigel D, Rubio-Somoza I, San Segundo B. The Arabidopsis miR396 mediates pathogen-associated molecular pattern-triggered immune responses against fungal pathogens. Sci Rep. 2017;7.

  17. Liang YX, Guan YH, Wang SX, Li YJ, Zhang ZH, Li H. Identification and characterization of known and novel microRNAs in strawberry fruits induced by Botrytis cinerea. Sci Rep. 2018;8.

  18. Gao X, Cui Q, Cao QZ, Liu Q, He HB, Zhang DM, et al. Transcriptome-Wide Analysis of Botrytis elliptica Responsive microRNAs and Their Targets in Lilium Regale Wilson by High-Throughput Sequencing and Degradome Analysis. Front Plant Sci. 2017;8.

  19. Jin WB, Wu FL. Characterization of miRNAs associated with Botrytis cinerea infection of tomato leaves. BMC Plant Biol. 2015;15.

  20. Jin WB, Wu FL, Xiao L, Liang GW, Zhen YX, Guo ZK, et al. Microarray-based Analysis of Tomato miRNA Regulated by Botrytis cinerea. J Plant Growth Regul. 2012;31(1):38–46.

    Article  CAS  Google Scholar 

  21. Tian X, Song LP, Wang Y, Jin WB, Tong FD, Wu FL. miR394 Acts as a Negative Regulator of Arabidopsis Resistance to B. cinerea Infection by Targeting LCR. Front Plant Sci. 2018;9.

  22. Wu FL, Qi JY, Meng X, Jin WB. miR319c acts as a positive regulator of tomato against Botrytis cinerea infection by targeting TCP29. Plant Sci. 2020;300.

  23. Bielewicz D, Kalak M, Kalyna M, Windels D, Barta A, Vazquez F, et al. Introns of plant pri-miRNAs enhance miRNA biogenesis. EMBO Rep. 2013;14(7):622–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Schwab R, Speth C, Laubinger S, Voinnet O. Enhanced microRNA accumulation through stemloop-adjacent introns. EMBO Rep. 2013;14(7):615–21.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Agranat-Tamir L, Shomron N, Sperling J, Sperling R. Interplay between pre-mRNA splicing and microRNA biogenesis within the supraspliceosome. Nucleic Acids Res. 2014;42(7):4640–51.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Jin W, Wu F. Characterization of miRNAs associated with Botrytis cinerea infection of tomato leaves. BMC Plant Biol. 2015;15:1.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  27. Jiang N, Meng J, Cui J, Sun G, Luan Y. Function identification of miR482b, a negative regulator during tomato resistance to Phytophthora infestans. Hortic Res. 2018;5:9.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  28. Castillo N, Pastor V, Chávez N, Arró M, Boronat A, Flors V, et al. Inactivation of UDP-Glucose Sterol Glucosyltransferases Enhances Arabidopsis Resistance to Botrytis cinerea. other. 2019;10.

  29. Michlewski G, Caceres JF. Post-transcriptional control of miRNA biogenesis. Rna. 2019;25(1):1–16.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Ratnadiwakara M, Mohenska M, Anko ML. Splicing factors as regulators of miRNA biogenesis - links to human disease. Semin Cell Dev Biol. 2018;79:113–22.

    Article  CAS  PubMed  Google Scholar 

  31. Newman MA, Thomson JM, Hammond SM. Lin-28 interaction with the Let-7 precursor loop mediates regulated microRNA processing. Rna. 2008;14(8):1539–49.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Piskounova E, Viswanathan SR, Janas M, LaPierre RJ, Daley GQ, Sliz P, et al. Determinants of microRNA processing inhibition by the developmentally regulated RNA-binding protein Lin28. J Biol Chem. 2008;283(31):21310–4.

    Article  CAS  PubMed  Google Scholar 

  33. Viswanathan SR, Daley GQ, Gregory RI. Selective blockade of MicroRNA processing by Lin28. Science. 2008;320(5872):97–100.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Fernandez N, Cordiner RA, Young RS, Hug N, Macias S, Caceres JF. Genetic variation and RNA structure regulate microRNA biogenesis. Nat Commun. 2017;8.

  35. Shen J, Ambrosone CB, Zhao H. Novel genetic variants in microRNA genes and familial breast cancer. Int J Cancer. 2009;124(5):1178–82.

    Article  CAS  PubMed  Google Scholar 

  36. Lu S, Sun YH, Shi R, Clark C, Li L, Chiang VL. Novel and mechanical stress-responsive MicroRNAs in Populus trichocarpa that are absent from Arabidopsis. Plant Cell. 2005;17(8):2186–203.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Lu S, Sun YH, Amerson H, Chiang VL. MicroRNAs in loblolly pine (Pinus taeda L.) and their association with fusiform rust gall development. Plant J. 2007;51(6):1077–98.

    Article  CAS  PubMed  Google Scholar 

  38. Lu S, Sun YH, Chiang VL. Stress-responsive microRNAs in Populus. Plant J. 2008;55(1):131–51.

    Article  CAS  PubMed  Google Scholar 

  39. Subramanian S, Fu Y, Sunkar R, Barbazuk WB, Zhu JK, Yu O. Novel and nodulation-regulated microRNAs in soybean roots. BMC Genomics. 2008;9:160.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  40. Arenas-Huertero C, Pérez B, Rabanal F, Blanco-Melo D, De la Rosa C, Estrada-Navarrete G, et al. Conserved and novel miRNAs in the legume Phaseolus vulgaris in response to stress. Plant Mol Biol. 2009;70(4):385–401.

    Article  CAS  PubMed  Google Scholar 

  41. Devers EA, Branscheid A, May P, Krajinski F. Stars and symbiosis: microRNA- and microRNA*-mediated transcript cleavage involved in arbuscular mycorrhizal symbiosis. Plant Physiol. 2011;156(4):1990–2010.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Mohorianu I, Schwach F, Jing R, Lopez-Gomollon S, Moxon S, Szittya G, et al. Profiling of short RNAs during fleshy fruit development reveals stage-specific sRNAome expression patterns. Plant J. 2011;67(2):232–46.

    Article  CAS  PubMed  Google Scholar 

  43. de Vries S, Kloesges T, Rose LE. Evolutionarily Dynamic, but Robust, Targeting of Resistance Genes by the miR482/2118 Gene Family in the Solanaceae. Genome Biol Evol. 2015;7(12):3307–21.

    Article  PubMed  PubMed Central  Google Scholar 

  44. Yang L, Mu X, Liu C, Cai J, Shi K, Zhu W, et al. Overexpression of potato miR482e enhanced plant sensitivity to Verticillium dahliae infection. J Integr Plant Biol. 2015;57(12):1078–88.

    Article  CAS  PubMed  Google Scholar 

  45. Feng J, Liu S, Wang M, Lang Q, Jin C. Identification of microRNAs and their targets in tomato infected with Cucumber mosaic virus based on deep sequencing. Planta. 2014;240(6):1335–52.

    Article  CAS  PubMed  Google Scholar 

  46. Zhu QH, Fan L, Liu Y, Xu H, Llewellyn D. Wilson I: miR482 regulation of NBS-LRR defense genes during fungal pathogen infection in cotton. PLoS One. 2013;8(12):e84390.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  47. Ouyang S, Park G, Atamian HS, Han CS, Stajich JE, Kaloshian I, et al. MicroRNAs suppress NB domain genes in tomato that confer resistance to Fusarium oxysporum. PLoS Pathog. 2014;10(10):e1004464.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  48. Rowe HC, Walley JW, Jason C, Chan KF, Katayoon D, Kliebenstein DJ, et al. Deficiencies in Jasmonate-Mediated Plant Defense Reveal Quantitative Variation in Botrytis cinerea Pathogenesis. PLoS Pathog. 2010;6(4):e1000861.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  49. Proietti S, Caarls L, Coolen S, Van Pelt JA, Van Wees SCM, Pieterse CMJ. Genome-wide association study reveals novel players in defense hormone crosstalk in Arabidopsis. Plant Cell Environ. 2018;41(10):2342–56.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Lorenzo O, Piqueras R, Sánchez-Serrano J, Solano R. ETHYLENE RESPONSE FACTOR1 Integrates Signals from Ethylene and Jasmonate Pathways in Plant Defense. Plant Cell. 2003;15(1):165–78.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Penninckx IAMA, Thomma BPHJ, Buchala, Metraux JP, Broekaert WF. Concomitant activation of jasmonate and ethylene response pathways is required for induction of a plant defensin gene in Arabidopsis. Plant Cell. 1998.

  52. Bertini L, Leonardi L, Caporale C, Tucci M, Cascone N, Berardino ID, et al. Pathogen-responsive wheat PR4 genes are induced by activators of systemic acquired resistance and wounding. Plant Sci. 2003;164(6):1067–78.

    Article  CAS  Google Scholar 

  53. Meng X, Jin W, Wu F. Novel tomato miRNA miR1001 initiates cross-species regulation to suppress the conidiospore germination and infection virulence of Botrytis cinerea in vitro. Gene. 2020;759:145002.

    Article  CAS  PubMed  Google Scholar 

  54. Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001;25(4):402–8.

    Article  CAS  PubMed  Google Scholar 

  55. Lorenz R, Bernhart SH, Siederdissen CHZ, Tafer H, Flamm C, Stadler PF, et al. ViennaRNA Package 2.0. Algorithm Mol Biol. 2011;6.

  56. Qi Y, Denli AM, Hannon GJ. Biochemical specialization within Arabidopsis RNA silencing pathways. Mol Cell. 2005;19(3):421–8.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

We would like to thank Xing Tian, Jingyi Qi and Yi Wang for their excellent technical assistance.

Funding

This study was supported by the Natural Science Foundation of Zhejiang Province (grant no. LY21C140004 and LY18C150010), and the 521 Talent Foundation of Zhejiang Sci-Tech University.

Author information

Authors and Affiliations

Authors

Contributions

FW and WJ designed this research. FW, JX, TG and DH performed the experiments and analyzed the data. FW and WJ wrote and revised this paper. All Authors have read and approved the final manuscript.

Corresponding author

Correspondence to Weibo Jin.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interest

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Wu, F., Xu, J., Gao, T. et al. Molecular mechanism of modulating miR482b level in tomato with botrytis cinerea infection. BMC Plant Biol 21, 496 (2021). https://doi.org/10.1186/s12870-021-03203-2

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12870-021-03203-2

Keywords