- Research article
- Open Access
MicroRNA-mediated susceptible poplar gene expression regulation associated with the infection of virulent Melampsora larici-populina
BMC Genomics volume 17, Article number: 59 (2016)
Rust caused by Melampsora larici-populina is one of the most damaging diseases of poplars. Rust is considered to be a model pathogen for genetic studies because both pathogen and host genomes are available. The poplar ‘Robusta’, whose general rust resistance is defeated by virulent rust E4, provides suitable host material for studies of the gene regulation involved in rust resistance/susceptibility. In this study, we investigated the microRNA-mediated susceptible poplar gene expression regulation associated with the infection of virulent rust. We were particularly interested in delineating the host-pathogen interactions with a specific focus on microRNAs (miRNAs).
To study the susceptibility of poplar to M. larici-populina, small RNA (sRNA) libraries, a degradome cDNA library and digital gene expression libraries were constructed for rust-inoculated and rust-free susceptible poplar ‘Robusta’ leaves through high-throughput sequencing. Altogether, 12,722 regulating interactions were identified. The results delineated the framework of post-transcriptional regulation of gene expression in the susceptible poplar, which was infected by the virulent rust. The results indicated that pathogen-associated molecular patterns (PAMPs) and PAMP-triggered immunity were induced by the infection of virulent rust E4 and that miRNAs still functioned at this stage. After this stage, miRNA-regulated R genes, such as TIR-NBS-LRR and CC-NBS-LRR, were not fully functional. Additionally, the rust-responsive miRNAs did not regulate the signaling component genes related to the salicylic acid pathway or the hypersensitive response.
We found that the defense-related post-transcriptional regulation of the susceptible poplar ‘Robusta’ functions normally only at the stage of PAMPs and PAMP-triggered immunity (PTI). More importantly, the miRNA-mediated post-transcriptional regulation of defense signal pathway genes were inactivated by the infection of virulent rust at the stage of effector-triggered susceptibility and during the following stages of salicylic acid and hypersensitive responses. This inactivation was the major characteristic of ‘Robusta’ susceptibility.
Rust caused by Melampsora larici-populina is one of the most severe diseases of poplars. Rust is indigenous to Eurasia and later spread to Australia  and North America . In the Far East, M. larici-populina infection occurs in a wide range of poplar species and hybrids and often causes severe damages [3, 4]. It has a very complex life cycle, producing 5 spore stages during its life cycle and alternating on larch (Larix spp.).
Poplars (Populus spp.) occur naturally in most parts of the northern hemisphere from subarctic to subtropical regions  and are widely planted beyond their natural range in both the southern and northern hemispheres . Poplars and their hybrids vary in rust resistance. Although hybrids of P. deltoides and P. nigra or P. trichocarpa were selected for their immunity to rust , breakdown of the resistance to M. larici-populina in some of these cultivars was detected upon the appearance of races E1 and E2 [8, 9]. Subsequently, many cultivars that had been highly resistant to rust were infected by a new race, E4 . Resistance and susceptibility are at opposite ends of the same spectrum, and research on disease resistance cannot be conducted without reference to susceptibility . Although great advances have been made in researching the genetic mechanisms of plant disease resistance, the mechanisms underlying plant disease susceptibility remain unclear.
Populus deltoides possesses qualitative resistance to M. larici-populina . Although hybrid vigor is generally high in these cultivars, most of the hybrids exhibit high rust susceptibility once the major resistance genes inherited from P. deltoides are defeated [12, 13]. In a previous study, hybrid poplar P. nigra × P. deltoides ‘Robusta’ was found to have non-race-specific resistance to M. larici-populina isolates collected from England , but the resistance was defeated by E4. Therefore, ‘Robusta’ provides suitable host material for studies of gene regulation involved in rust resistance/susceptibility.
MicroRNAs (miRNAs) play important roles in plant post-transcriptional gene regulation by repressing the translation of target mRNAs or targeting them for cleavage . In plants, miRNAs regulate diverse processes, such as development , abiotic stress tolerance , and biotic stress . Several studies have indicated that miRNAs play critical roles in disease resistance or susceptibility responses . For example, several plant miRNAs, such as miR160, miR167, miR398 and miR1885, are responsive to biotic stress during both resistance and susceptibility interactions [17–19]. Additionally, in susceptible poplar, all fungi-responsive miRNAs have been found to be up-regulated in response to the infection of a canker pathogen .
The genome sequences of both poplar (P. trichocarpa) and its rust pathogen M. larici-populina have now been published . Thus, the poplar–poplar rust pathogen system is the only system in which the genome sequences of both tree host and rust pathogen are available. Although poplar genomic analyses have revealed the genetic potential of susceptibility and transcriptomics has allowed for the deciphering of gene expression in space and time, little is known about the post-transcriptional regulation of resistance gene expression. Recent evidence has suggested that post-transcriptional silencing involving small RNAs of cellular rather than pathogenic origin might have broad implications in potentiating basal defense and race-specific resistance to microbes in plants . The present study was conducted to decipher the regulatory functions of miRNA in a susceptible poplar and to unravel the complexity of post-transcriptional regulation of resistance genes in the susceptible poplar ‘Robusta’.
Small RNA sequencing profile
To study the post-transcriptional regulation associated with poplar susceptibility to rust, miRNA accumulation was investigated by small RNA sequencing. Two small RNA (sRNA) libraries constructed from a combination of seven selected time points of the ‘Robusta’ leaves of -rust (rust-free) and + rust (E4 rust-inoculated) were sequenced on an Illumina GAIIx system. To acquire more details of the sRNA background in one sequencing, seven time points covering the major developmental transitions of the urediniospores in poplar leaves were selected, consistently with previous studies [23, 24]. A total of 8,840,983 (Additional file 1: Figure S1 A) and 9,736,262 (Additional file 1: Figure S1 B) raw sequences were generated from the -rust and + rust libraries, respectively. Blast homology searches were performed against M. larici-populina (version 1.0, JGI website http://genome.jgi-psf.org/Mellp1/), and the mapped reads were removed from the + rust library.
The length distribution of the mappable reads showed that 15-nt reads (18 %) were the most abundant sequences in the + rust sRNA library and that 21-nt reads (17.5 %) were the second most abundant sequences (Additional file 1: Figure S2). In contrast, the -rust sRNA library showed that 21-nt reads were the most abundant sequences (Additional file 1: Figure S2), constituting 19.3 % of all reads. This latter result is in agreement with results from sRNA studies of P. balsamifera  and Vitis vinifera , in which the 21-nt sRNAs showed the highest abundance. However, this result differs from the results of studies of other plants, in which 24-nt sRNAs were most abundant [27–28].
By removing the redundant sequences from mappable reads, 3,820 and 4,070 unique sequences were identified from the -rust and + rust libraries, respectively. After removing the consensus sequences, a total of 5,138 unique miRNAs were identified in these two sRNA libraries. These unique miRNAs were classified into two groups: known miRNAs and predicted miRNAs (Additional file 1: Figure S1 C and D, Additional file 1: Table S1). Known miRNAs were further divided into two subgroups: Gp1a and Gp1b (Additional file 1: Text S1). Predicted miRNAs were further divided into six subgroups: Gp2a, Gp2b, Gp3a, Gp3b, Gp4a and Gp4b (Additional file 1: Text S1). By removing the no-hairpin subgroups and the subgroups that did not map to the poplar genome, Gp1a, Gp1b, Gp2a and Gp4a were identified as candidate miRNA subgroups. Among these candidate miRNA subgroups, 1,474 and 1,475 miRNAs were obtained from the -rust and + rust sRNA libraries, respectively. The potential novel miRNAs identified by deep sequencing were deposited in The European Nucleotide Archive (WEBIN ID number: Hx2000050201).
According to the eleven features of the miRNA hairpins (LC Sciences, Houston, TX, USA, Additional file 1: Text S2), the secondary structure of genes corresponding to the potential novel miRNAs (Gp4a) were identified (Table 1), and 4 predicted secondary structures of the potential novel miRNAs are shown as examples (Additional file 1: Figure S3).
Identification of rust-responsive miRNA in P. nigra × P. deltoides
MiRNAs regulate gene expression in plants by promoting the degradation of or repressing the translation of target mRNAs, both transcriptionally and post-transcriptionally [18, 28]. To associate rust fungi M. larici-populina infection with an alteration of poplar miRNA accumulation, the differential expression of miRNAs in the two libraries was compared using the log2- (+rust/-rust) fold changes of the normalized miRNA copy numbers (Fig. 1a). Log2- (+rust/-rust) fold changes >1 were designated as indicating ‘up-regulated miRNA’ , and changes < −1 were designated as indicating ‘down-regulated miRNA’ (Fig. 1b).
Analysis of the differential expression of 5,138 miRNAs showed that 2,126 (41.38 %) were expressed in both the -rust and + rust libraries (co-expression, Fig. 1c), whereas 1,335 (25.98 %) and 1,677 (32.64 %) were preferentially expressed in -rust and + rust libraries, respectively (Fig. 1c). Analysis of co-expressed miRNAs resulted in the identification of 541 miRNAs (out of 2,126; 25.45 %) that showed significant differential expression between the -rust and + rust libraries. Of these, 376 (73.15 %) were up-regulated, and 165 (32.10 %) were down-regulated after the inoculation of rust.
In the four candidate miRNA subgroups, 1,081 miRNA members were responsive to rust infection (log2-fold changes >1 or < −1). Restricting the analysis to only those miRNAs with more than 5 raw reads in both of the two sRNA libraries revealed changes of only 61 miRNAs as shown in Fig. 2. In the subgroup Gp1a, 7 miRNAs were up-regulated and 7 miRNAs were down-regulated in response to rust infection (Fig. 2a). With the exception of 3 up-regulated miRNAs, no differentially expressed miRNAs were found in the Gp1b subgroup (Fig. 2b). In the subgroup Gp2a, 5 miRNAs were up-regulated and 8 miRNAs were down-regulated in response to rust infection (Fig. 2c). In total, 12 miRNAs were up-regulated and 7 were down-regulated in response to rust infection in the Gp4a subgroup and were potentially novel miRNAs (Fig. 2d). The most robustly changed miRNA was PC-5p-310376_3 (Gp4a, Fig. 2d), with a peak log2-fold change of 12.641.
Identification of miRNA targets
To identify genes targeted by miRNAs, a degradome cDNA library from + rust poplar leaves was constructed. In total, 8,409,874 raw, short sequencing reads were obtained from this library (Additional file 1: Table S2), and 5,077,704 unique reads mapped to the Populus mRNA database (http://www.ncbi.nlm.nih.gov/nuccore/). A list of all of the confidently detected mRNA targets along with the corresponding alignments for the miRNA–mRNA interactions were generated by the pipeline. All five categories (Additional file 1: Text S3)  of the sliced-target transcripts were found. According to the 1,475 miRNAs identified in the + rust sRNA library, 947 (64.20 %) miRNAs potentially targeted to 8,463 poplar genes, and a total of 12,722 regulating interactions were found.
These targets included 31 plant defense genes that were regulated by 26 miRNA (Table 2), such as the disease resistance protein TIR-NBS-LRR, regulated by miR482a (Fig. 3f); the HXXXD-type acyl-transferase family protein gene, regulated by miR6468 (Fig. 3j); and the wound-responsive family protein, regulated by miR394-p3. Additionally, these targets included 27 transcription factors that were regulated by 27 miRNAs (Additional file 1: Table S3), such as the AP2/ERF superfamily transcription factor gene, regulated by miR172g-3p (Fig. 3a); the GRAS family transcription factor gene, regulated by miR171i-p3 (Fig. 3 b); the MYB transcription factor gene, regulated by miR828a (Fig. 3d); and the WRKY70 gene, regulated by miR172b-5p (Fig. 3e).
Digital gene expression
To characterize the target gene expression pattern, two digital gene expression (DGE) libraries of -rust and + rust were constructed and sequenced. There were 16,409,618 and 18,186,321 reads sequenced from the -rust and + rust libraries, respectively. More than 99 % of the raw reads passed the filter cut-off, resulting in 18,117,197 and 16,338,785 clean reads. In total, 73,013 unique genes were obtained, and 65,822 of those unique genes were mappable. Of the unique reads, 1,331,698 (25.07 %) and 1,419,523 (24.85 %) sense reads were mapped to one gene for the -rust and + rust libraries, respectively. The remaining 73.38 % of the unique reads for the –rust library and the remaining 73.60 % for the + rust library matched to multiple genes. The remaining reads were unique antisense or both sense and antisense reads that were mappable to genes.
To compare the differential expression patterns between the -rust and + rust libraries, we normalized (by the means of TPM) the read distribution for the gene expression level in each library to construct an effective library size and determined the significance of differentially expressed transcripts with x 2 (P ≤ 0.01). According to the annotation of these genes, the rust infection-regulated gene expression was analyzed (Fig. 4).
According to the results of the DGE analysis, the targets of the potential novel miRNAs with altered expression were analyzed (Tables 1 and Additional file 1: S4). Then, the log2-fold changes of plant defense-related genes and their regulating miRNAs were analyzed (Table 2).
On the basis of differential expression analysis  and control of the false discovery rate (FDR) , 282 differentially expressed target genes were identified. To compare these results with the real-time quantitative polymerase chain reaction (RT-qPCR) results, the log2- (+rust/-rust) fold changes >1 were designated as indicating ‘up-regulated genes’, and the changes < −1 were designated as indicating ‘down-regulated genes’.
Gene Ontology (GO) analysis was performed for all of the differentially expressed genes (Additional file 1: Table S5). Altogether, 282 differentially expressed target genes were linked to 56 GO terms. Once the GO terms of interest were identified for the unigenes, the differentially expressed genes and their associated miRNAs were analyzed (Additional file 1: Text S5). From the GO terms, we found that many processes, such as response to stress (GO:0006950), heat shock protein binding (GO:0031072), unfolded protein binding (GO:0051082), oxidation-reduction process(GO:0055114), were prominently down-regulated and that the processes of response to oxidative stress (GO:0006979), metabolic process (GO:0008152), biosynthetic process (GO:0009058), and cell wall macromolecule catabolic process (GO:0016998) were up-regulated.
The impact of rust infection on post-transcriptional regulation in the susceptible poplar
From the analysis of the miRNA, degradome cDNA and DGE libraries, 12,722 regulatory interactions (947 miRNAs to 8,463 target genes) were identified (Fig. 5a). Among the 947 miRNAs involved, 250 down-regulated (log2-fold changes <1) and 301 up-regulated (log2-fold changes >1) rust-responsive miRNAs were associated with 4,099 and 4,781 regulatory interactions (Fig. 6), respectively. Among the 4,099 regulatory interactions targeted by 250 down-regulated miRNAs, 59 (1.439 %) target genes were up-regulated (log2-fold changes >1), 61 (1.488 %) target genes were down-regulated (log2-fold changes < −0.5), 3,967 (96.77 %) target genes were unchanged (log2-fold changes > −1 and <1), and the other target genes (12 regulating interactions) were undetected by DGE analysis. Among the 4,781 regulatory interactions targeted by 301 up-regulated miRNAs, 44 (0.92 %) target genes were down-regulated, 73 (1.53 %) target genes were up-regulated, 4,648 (97.22 %) target genes were unchanged, and the other target genes (16 regulating interactions) were undetected by DGE analysis. The resistance-related miRNAs and their target genes were found in poplar ‘Robusta’, but not all of those miRNAs and their target genes were in compliance with this negative regulation (Additional file 1: Figure S4, Additional file 1: Table S4, transcription model M2b).
Furthermore, we analyzed the post-transcriptional regulation of 82 miRNAs to 171 resistance genes, such as the genes encoding intracellular nucleotide-binding-site (NBS) -leucine-rich repeat (LRR) proteins (Fig. 5b, A1 and A2), signaling cascades (Fig. 5b, B1 to B5), transcription factors (Fig. 5b, C1 to C7), reactive oxygen (Fig. 5b, D1 to D7) and pathogenesis-related (Fig. 5b, E1 to E10) (Additional file 1: Table S4). Out of these 82 miRNAs, 28 (34 %) regulated only one target gene, whereas the other 54 (66 %) regulated more than one target gene. In contrast, 132 (77 %) resistance genes were regulated by one miRNA, whereas the remaining genes (23 %) were regulated by more than one miRNA.
Altogether, 226 resistance regulating interactions were identified, and these regulating interactions were further divided into five regulation modes. The first mode included 73 regulating interactions (32 %). The expression of target genes was unchanged (Additional file 1: Table S4, regulation mode M1b) or up-regulated (Additional file 1: Table S4, regulation mode M1a) in the first mode, although the accumulation of miRNA was increased. Only 9 regulating interactions (4 %) were identified in the second mode. In this mode, the accumulation of miRNA was unchanged, but the expression of target genes was up-regulated (Additional file 1: Table S4, regulation mode M2a) or strongly suppressed (Additional file 1: Table S4, regulation mode M2c). The maximum amount of regulating interactions, a total of 89 interactions (39 %), was identified in the third mode. In this mode, the accumulation of miRNA was decreased, but the expression of target genes was down-regulated (Additional file 1: Table S4, regulation mode M3c) or unchanged (Additional file 1: Table S4, regulation mode M3b). The fourth mode included 31 regulating interactions (14 %, Additional file 1: Table S4, regulation mode M2b), and both the miRNA and the target gene were unresponsive to the rust infection. The last mode, which is consistent with the classical post-transcriptional regulation theory, included 24 regulating interactions (11 %). In the last mode, the accumulation of miRNA and the expression of target genes were in compliance with negative correlation (Additional file 1: Table S4, regulation modes M1c and M3a).
Validation and temporal pattern of miRNA expression and target genes
To avoid missing significant differential changes in the mixed samples sequencing, RT-qPCR was performed for 10 differentially expressed miRNAs (log2 > 1 or < −1) and target genes at 2 h post inoculation (hpi) of rust, 6 hpi, 12 hpi, 1 day post inoculation (dpi), 2 dpi, 4 dpi and 7 dpi. These miRNAs were selected based on principles Additional file 1: Text S5 and are listed in Additional file 1: Table S6. The expression patterns of the RT-qPCR results showed details that had not been found in the sRNA libraries and DGE sequencing.
According to previous studies, before 1 dpi, infection hyphae extend into the mesophyll and differentiate into the first haustorial structures [32, 33]. Then, biotrophic growth continues, and the fungal biomass greatly increases between 2 and 4 dpi [34, 35]. During this time period, the expression of miRNA volatility continually changes over time, and we found 5 expression patterns related to the development of rust occurring on poplar leaves (Additional file 1: Figure S5; α, β, γ, δ and ε patterns, Additional file 1: Text S6). The expression of targets (Additional file 1: Figure S6) did not accompany changes in miRNA for all 7 time periods (Fig. 7), but several negative interactions were found for the α-pattern miRNA (miR172b-5p targeting the WRKY70 gene, miR171i-p3 targeting the GRAS family transcription factor gene, miR393_R + 2 targeting the F-box/RNI-like superfamily gene, miR828a targeting the MYB transcription factor gene and miR172b-5p targeting the WRKY70 gene) and the γ-pattern miRNA (miR171e-p5 targeting the epsin N-terminal homology domain-containing protein gene) at 1 dpi; in these cases, the targets were up-regulated, whereas the miRNAs were down-regulated (Fig. 7, boxed with red dashed lines).
Framework of the susceptible poplar defense signal pathway
The miRNA and target global maps indicated that most of the related KEGG (Kyoto Encyclopedia of Genes and Genomes) pathways (Additional file 1: Figure S7) were down-regulated or up-regulated after rust infection, whereas few target genes were responsive to rust infection in the susceptible poplar. The other pathways, such as plant hormone signal transduction, cell cycle and plant-pathogen interaction, showed the same pattern; i.e., the miRNAs reacted strongly to the infection, but the target genes reacted mildly. The pathway of Plant-pathogen Interaction indicated that the miRNAs related to fungal pathogen-associated molecular patterns (PAMPs, including reactive oxygen species), Hypersensitive response (HR) and fungal effectors were responsive to rust infection. However, only 4 of 31 genes were up-regulated in the same pathway.
According to the regulatory interactions, an miRNA-mediated susceptible poplar defense signal pathway related to rust infection was revealed (Additional file 1: Table S7). The receptor-like kinase genes RLK 25, RLK 26, RLK 29, RLK 34 and RLK 42 were up-regulated, whereas the regulating miRNAs were down-regulated by rust infection. RLK 14 was up-regulated, but the regulating miRNAs remained unchanged. RLK 3 was up-regulated, as were the regulating miRNAs up-regulated.
The NBS-LRR family of resistance (R) proteins induced effector-triggered susceptibility (ETS). The results of degradome cDNA library sequencing indicated that 25 TIR-NBS-LRR family genes were regulated by 9 miRNAs. Although 7 miRNAs were down-regulated significantly, none of the corresponding 18 target genes were up-regulated. Two miRNAs were up-regulated, whereas the target TIR-NBS-LRR family genes were unchanged. Only 1 TIR-NBS-LRR gene family was up-regulated after the infection of rust fungi, whereas the regulating miRNA was unchanged. Additionally, 4 unchanging CC-NBS-LRR family genes were regulated by 2 down-regulated miRNAs. The suppression of the TIR-NBS-LRR-regulating miRNA indicated that recognition occurred between the rust effectors and the poplar R proteins but that the signal promoting the expression of R proteins was blocked. This process explains why ‘Robusta’ cannot prevent rust infection at the ETS stage.
Many downstream defense signal pathways were also found. First, pathogenesis-related (PR) genes, which were triggered by effector signals and transcription factors, were investigated (Additional file 1: Table S7). MiR5517-p5 was up-regulated, and its target PR gene was suppressed. Two pairs of PR genes and their regulating miRNA were unchanged. MiR3447-p3 was strongly up-regulated, yet its target PR gene was unchanged.
Second, the pathway related to salicylic acid (SA), which is the main signal molecule for systemic acquired resistance (SAR) and HR, did not respond to rust infection, although many of the regulating miRNAs were found to be down-regulated or up-regulated after infection. The enhanced disease susceptibility 1 gene (EDS1, regulated by the promotion of miR3434-p3) and the phytoalexin-deficient 4 gene (PAD4, regulated by the suppression of miR4248b-p3), which act to induce genes upstream of SA accumulation, did not respond to rust infection. Ankyrin-repeat protein NPR1, which acts downstream of SA and promotes the expression of the SAR-associated genes PR-1, BGL2, and PR-5 in Arabidopsis , also did not respond to rust infection. The expression of the ‘Robusta’ enhanced disease susceptibility 5 (EDS5) gene was up-regulated slightly (it is also regulated by the suppression of miR4248b-p3), whereas SA induction–deficient 2 (SID2) was unchanged (regulated by the up-regulated PC-5p-1150271_1) after rust infection, although EDS5 and SID2 function upstream of SA accumulation in Arabidopsis . The programmed cell death-related gene also did not respond to rust infection. These results suggest that the signal transduction pathways are blocked, which typically leads to an HR-related accumulation of SA and programmed cell death to halt pathogen invasion. Third, with the exception of the suppression of mitogen-activated protein kinase 4 (MPK4, regulated by the down-re down-regulated PC-3p-272434_3) down-regulated and the promotion of 12-oxpophytodienoic acid reductase 3 (OPR3, regulated by the down-regulated PC-5p-1312398_1) up-regulated, no genes related to jasmonic acid (JA) or ethylene (ET) responded to rust infection, including coronatine-insensitive protein 1 (COI1, undetected), jasmonic acid-amido synthetase 1 (JAR1, regulated by the down-regulated miR393b-p5_1) and Ethylene-insensitive protein 2 (EIN2, undetected). Normally, JA- and ET-dependent resistance is up-regulated by lesion formation in A. thaliana. MPK4, COI1 and JAR1 are required to transduce the JA signal, whereas EIN2 is required to transduce the ET signal .
The miRNA-mediated post-transcriptional regulation of plant resistance to biotic stress has been described in several plant pathosystems [22, 36–41]. In contrast, the molecular mechanisms of plant susceptibility to biotic stress and the post-transcriptional regulation of such susceptibility are poorly studied. Identifying the genetic basis of host susceptibility/resistance is a prerequisite for understanding the interaction between pathogen and plant host. Based on miRNA, degradome cDNA and DGE library analyses, we investigated the regulatory functions of miRNA in a susceptible poplar under rust attack. We found that the genes involved in the miRNA-mediated post-transcriptional regulation of the defense signal pathway were inactivated after infection by E4 in ‘Robusta’, whereas this regulation was functional during E1 infection. This inactivation was the major characteristic of ‘Robusta’ susceptibility.
The small RNA sequencing revealed a significant peak in abundance of 15-nt reads, with 835 unique 15-nt miRNAs were found in the + rust sRNA library. Moreover, 120 miRNAs targeted 113 unigenes, which included many disease-related genes, such as the AP2/B3 transcription factor protein family, the disease resistance CC-NBS-LRR class protein family, the disease resistance TIR-NBS-LRR class protein family, cellulose synthase 6 genes, LRR and NB-ARC domain-containing disease-resistance protein genes and stress-inducible protein genes. However, unexpectedly, none of these disease-related target genes responded to rust infection. These results indicate that the 15-nt miRNAs responded to rust infection and showed a significant increase but that they failed to regulate their target genes in the susceptible poplar.
Basal defense and R gene-mediated resistance are the two branches of the defense response carried out by plants under pathogen stress . Basal defense marks the first line of defense and consists of a set of defined receptors referred to as pattern recognition receptors (PRRs), which recognize conserved, slowly evolving MAMPs or PAMPs . The receptor-like kinases, which are part of the initial stages of PAMPs and PTI, were up-regulated in this compatible (susceptible poplar vs. virulent rust) pathosystem, and the miRNAs related to these stages were responsive to rust infection. These results indicate that the poplar defense commences with pathogen-associated molecular patterns (PAMPs) and PAMP-triggered immunity (PTI) and that the miRNAs are still functional at this stage.
The second branch of the defense response acts largely inside the cell and uses the polymorphic NB-LRR protein products encoded mostly by R genes . Although the basal defense of the poplar was almost completely destroyed in this compatible pathosystem, R gene-mediated resistance was not fully functional. Moreover, the post-transcriptional regulation of the R gene did not occur. Only 1 TIR-NBS-LRR gene family, regulated by miR472b, was up-regulated after rust fungi infection, although poplar genomes contain 64 TIR-NBS-LRR and 119 CC-NBS-LRR immune receptors encoded by R genes that recognize the pathogen infection and trigger resistance responses .
Twelve miRNAs that guide the cleavage of 4 CC-NBS-LRR and 25 TIR-NBS-LRR transcripts were identified (Additional file 1: Table S7). Two TIR-NBS-LRR and 2 CC-NBS-LRR genes were regulated by miRNAs that belong to the mir482 superfamily; this result is consistent with those from previous studies [46, 47]. Both miR482b and its target CC-NBS-LRR were unresponsive to rust infection. In contrast, although miR482 was up-regulated (log2-fold change =1.77), 6 TIR-NBS-LRR target genes were unchanged (log2-fold change values ranging from 0.04 to 0.49). As one of the most important plant immune system branches , NBS-LRR-mediated disease resistance is effective against obligate biotrophic pathogens and hemi-biotrophic pathogens . Our results showed that the susceptible poplar NBS-LRR-related genes were not up-regulated by the rust biotrophic pathogen. Interestingly, another study showed that NBS-LRRs are not substantially induced by an incompatible rust strain in P. trichocarpa × P. deltoides . These findings suggest that the lack of substantial induction of NBS-LRRs by infection is not a characteristic specific to susceptible poplar but is common to both susceptible and resistant poplars.
All of the signaling component genes, including the genes related to ET, SA, JA and HR, were unresponsive to rust infection, which suggests that the regulation of miRNA-based HR signaling was blocked. Normally, plants produce phytohormones to cope with infection by diverse pathogens . SA is involved in resistance to biotrophic pathogens, whereas JA and ET mediate resistance primarily to necrotrophic pathogens . These results are consistent with our expectations and the gene-for-gene hypothesis, because the heterozygous resistant (Rr) poplar and the homozygous virulent (avr) rust were used in this study. Most of the studied plant pathosystems are related to the R-Avr gene-mediated disease response during pathogen invasion . However, in our study, the plant pathosystem of susceptible poplar vs. virulent rust (an R-avr gene-mediated disease response) provides a new perspective for understanding plant resistance/susceptibility. However, our findings provide only initial insight into this topic, and a great effort will be required for a full understanding of the post-transcriptional regulation related to plant resistance/susceptibility.
The results of the sRNA, degradome cDNA and DGE library analyses indicated that the defense-related post-transcriptional regulation of the susceptible poplar ‘Robusta’ only function normally during the initial stages of PAMPs and PTI. More importantly, the miRNA-mediated post-transcriptional regulation of defense signal pathway genes were inactivated after infection by virulent rust at the stage of effector-triggered susceptibility and at the subsequent stages of salicylic acid and HR response. This inactivation is the major characteristic of ‘Robusta’ susceptibility.
The filial generation of virulent E4 and avirulent E1 was used in this study. Rust-infected poplar leaves were collected from P. trichocarpa cv. trichobelat Markington (northern England) and Alice Holt, Surrey (southern England). Rust isolates were derived from single uredinial pustules according previously reported methods . The rust spores were stored at −20 °C.
Plant materials and inoculation procedure
The 1-year-old dormant cuttings of the hybrid poplar (P. nigra × P. deltoides) ‘Robusta’ were used as a source of plant tissue. ‘Robusta’ was chosen because a previous study has found that this cultivar showed non-race-specific resistance to M. larici-populina isolates collected from England .
‘Robusta’ plants were grown in pots containing a sand-peat (50:50, v/v) mixture and were watered daily with deionized water under 16-h/8-h photoperiod in greenhouse conditions and 60-70 % relative humidity for 10 weeks. At the time of inoculation, young trees were 70–100 cm tall and bore 10 to 14 fully expanded leaves.
The inoculation procedure was performed as described by Pei et al.  with some modifications. In brief, fully expanded leaves from leaf plastochrony index (LPI) 5 to 9 were detached from ‘Robusta’ plants and spray-inoculated on their abaxial surface with a rust spore suspension in deionized water containing 0.004 % Tween 20 (1 drop in 100 ml) adjusted to 100,000 spores ml−1, or with deionized water containing 0.004 % Tween 20 as a control. After inoculation, the inoculated leaves were incubated in a growth chamber at 16 °C with 16 h day−1 illumination (80 uE m−2 s−1) for various durations. Each treatment was replicated three times. The samples harvested at different time points (2 hpi, 6 hpi, 12 hpi, 1 dpi, 2 dpi, 4 dpi, and 7 dpi) in the different treatments were then immediately snap-frozen in liquid nitrogen and maintained at −80 °C for further nucleic acid isolation.
Total RNA preparation
Total RNA was extracted from frozen leaves of ‘Robusta’ with the CTAB method. Two sets of total RNA were prepared, with one derived from the original RNA pool prepared from E4 rust-inoculated leaves (+rust) at the 2 hpi, 6 hpi, 12 hpi, 1 dpi, 2 dpi, 4 dpi, and 7 dpi time points and the other from the RNA pool derived from rust-free leaves (−rust, the control) at the same time points.
sRNA, degradome cDNA and DGE library preparation and sequencing
Two sRNA libraries were generated from the -rust and + rust RNA pools using the Illumina TruSeq Small RNA Preparation Kit according to Illumina’s TruSeq Small RNA Sample Preparation Guide 1.
The degradome cDNA library of + rust was prepared following procedures previously described (Additional file 1: Text S7) [48,51-53]. Two DGE libraries representing -rust and + rust were prepared according to the Illumina/Solexa standard protocol (directional mRNASeq sample preparation part #15018460 Rev. A, October 2010. Illumina, San Diego, CA, USA) (Additional file 1: Text S8).
The purified sRNA, degradome cDNA and DGE libraries were used for cluster generation on Illumina’s Cluster Station and then sequenced on the Illumina GAIIx platform by LC Sciences (Houston, TX, USA). To analyze the differential expression, three biological replicates were performed with the sRNA, degradome cDNA and DGE libraries. Each sample of degradome cDNA was sequenced for 80 cycles on one lane of the Illumina GAIIx platform (Illumina, San Diego, CA, USA) with the 2 × 100 bp module. DGE libraries were sequenced using the 1 × 36 bp module.
Raw sRNA library sequencing reads were obtained using Illumina’s Sequencing Control Studio software version 2.8 (SCS v2.8) following real-time sequencing image analysis and base-calling by Illumina’s Real-Time Analysis version 1.8.70 (RTA v1.8.70).
Raw sequencing reads of the degradome cDNA library were obtained using Illumina’s Pipeline version 1.5 software following sequencing image analysis by the Pipeline Firecrest Module and base-calling by the Pipeline Bustard Module.
A proprietary pipeline script, ACGT101-miR version 4.2 (LC Sciences, Houston, TX, USA), was used for the sRNA library sequencing data analysis to obtain mappable sequences from raw sequencing data and to map the miRNA-mappable unique sequences to pre-miRNA, mature miRNA or genome.
A global normalization was used to correct sRNA copy numbers among the -rust and + rust libraries (Additional file 1: Text S10). The differential gene expression of the DGE libraries was analyzed according to the method of Simon Anders, .
RT-qPCR of mature miRNAs
RT-qPCR was performed with the One Step Prime-Script miRNA cDNA Synthesis Kit and the SYBR Premix ExTag II reagent kit (TaKaRa, Dalian, China) using the Stratagene Mx3000P qPCR system (Agilent, USA) to validate the expression levels of miRNAs after the rust inoculation at 2 hpi, 6 hpi, 12 hpi, 1 dpi, 2 dpi, 4 dpi and 7 dpi. All of the primers used in this study are listed in Additional file 1: Table S8. The 5.8S ribosomal RNA was used as the internal control . This experiment was performed on three biological replicates. Quantification of the RT-qPCR results was conducted as previously described . To compare the RT-qPCR results with the results of the high-throughput sequencing, the RT-qPCR results were normalized as log2- (+rust miRNA copy/-rust miRNA copy) fold changes.
Availability of supporting data
All the supporting data are included as additional files as “Additional file 1”.
Coronatine-insensitive protein 1
Digital gene expression
Days post inoculation
Enhanced disease susceptibility 1
Enhanced disease susceptibility 5
Ethylene-insensitive protein 2
Hours post inoculation
Jasmonic acid-amido synthetase 1
Kyoto Encyclopedia of Genes and Genomes
Leaf plastochrony index
Pathogen-associated molecular patterns
Pattern recognition receptors
The PAMP-triggered immunity
Quantitative real time polymerase chain reaction
Systemic acquired resistance
Van Kraayenoord CWS, Laundon GF, Spiers AG. Poplar rusts invade New Zealand. Plant Dis Reporter. 1974;58:423–7.
Newcombe G, Chastagner GA. First report of the Eurasian poplar leaf rust fungus, Melampsora larici-populina, in North America. Plant Dis. 1993;77:532–5.
Pei MH, Shang YZ. Study on the leaf rust of cathay poplar caused by Melampsora larici-tremulae Kleb [in Chinese]. J Northeastern Forestry Institute. 1984;12:40–9.
Pei MH, Ruiz C, Bayon C, Hunter T, Lonsdale D. Pathogenic variation in poplar rust Melampsora larici-populina from England. Eur J Plant Pathol. 2005;111:147–55.
Bean WJ. Trees and shrubs hardy in the British Isles (Vol. III, 8th edn). London: John Murray Ltd.; 1976.
Boyle J, Winjum J, Kavanagh K, Jensen E. Planted forests: contributions to the quest for sustainable societies. The Netherlands: Kluwer Academic Publishers; 1999. p. 90.
Steenackers M. Breeding poplars for rust resistance -recent advances. Communications of the Faculty of Agricultural Sciences of the State University of Ghent, Belgium. 1988; 53:417–22
Steenackers V. Nouvelle race physiologique de Melampsora larici-populina en Belgique (communication provisoire). In: FAO/CIP, 22e Réunion de Groupe de Travail des Maladies. Italy: Casale Monferrato; 1982. p. 6–10.
Pinon J, van Dam BC, Genetet I, De Kam M. Two pathogenic races of Melampsora larici-populina in north-western Europe. Eur J Forest Pathol. 1987;17:47–53.
Steenackers M, Steenackers V, Delporte T. A new race of M. larici-populina in Belgium. In: International Poplar Commission, Turkey: Izmit. 1994.p.3-7.
Nancy AE. Plant disease susceptibility genes? Plant Cell. 2002;14:1983–6.
Dowkiw A, Jorge V, Villar M, Voisin E,Guérin V, Faivre-Rampant P, et al. Breeding poplars with durable resistance to Melampsora larici-populina leaf rust: a multidisciplinary approach to understand and delay pathogen adaptation. In: Sniezko RA, Yanchuk AD, Kliejunas JT, et al., editors. Proceedings of the fourth international workshop on the genetics of host-parasite interactions in forestry: Disease and insect resistance in forest trees. Albany: Department of Agriculture; 2012. p. 3–40.
Baohong Z, Xiaoping P, Cobb GP, Anderson TA. Plant microRNA: A small regulatory molecule with big impact. Dev Biol. 2006;289:3–6.
Jones-Rhoades MWW, Bartel DPP, Bartel B. MicroRNAs and their regulatory roles in plants. Ann Rev Plant Biol. 2006;57:19–53.
Jones-Rhoades MWW, David PB. Computational identification of plant microRNAs and their targets, including a stress-induced miRNA. Mol Cell. 2004;14:787–99.
Padmanabhan C, Zhang X, Jin H. Host small RNAs are big contributors to plant innate immunity. Curr Opin Plant Biol. 2009;12:465–72.
Fahlgren N, Howell MD, Kasschau KD, Chapman EJ, Sullivan CM, Cumbie JS, et al. High-throughput sequencing of Arabidopsis microRNAs: evidence for frequent birth and death of MIRNA genes. PLoS One. 2007;2:e219. doi:10.1371/journal.pone.0000219.
Jagadeeswaran G, Saini A, Sunkar R. Biotic and abiotic stress down-regulate miR398 expression in Arabidopsis. Planta. 2009;229:1009–14.
He XF, Fang YY, Feng L, Guo HS. Characterization of conserved and novel microRNAs and their targets, including a TuMV-induced TIR-NBS-LRR class R gene-derived novel miRNA in Brassica. FEBS Lett. 2008;582:2445–52.
Zhao JP, Jiang XL, Zhang BY, Su XH. Involvement of microRNA-mediated gene expression regulation in the pathological development of stem canker disease in Populus trichocarpa. PLoS One. 2012;9:e44968. doi:10.1371/journal.pone.0044968.
Petre B, Morin E, Tisserant E, Hacquard S, Da Silva C, Poulain J, et al. RNA-seq of early-infected poplar leaves by the rust pathogen Melampsora larici-populina uncovers PtSultr3;5, a fungal-induced host sulfate transporter. PLoS One. 2012;8:e44408. doi:10.1371/journal.pone.0044408.
Olivier V. Post-transcriptional RNA, silencing in plant-microbe interactions: a touch of robustness and versatility. Curr Opin Plant Biol. 2008;11:464–70.
Hacquard S, Delaruelle C, Legué V, Tisserant E, Kohler A, Frey P, et al. Laser capture microdissection of uredinia formed by Melampsora larici-populina revealed a transcriptional switch between biotrophy and sporulation. Mol Plant-Microbe Interact. 2010;23:1275–86.
Duplessis S, Hacquard S, Delaruelle C, et al. Melampsora larici-populinatranscript profilingduring germination and timecourse infection of poplar leaves reveals dynamic expression patternsassociated with virulence and biotrophy. Mol Plant-Microbe Interact. 2011;7:808–18.
Abdelali B, Phillip KW, Scott DL, Claude WD, John EC. Conservation and divergence of microRNAs in Populus. BMC Genomics. 2007;8:481.
Vitantonio P, Gyorgy S, Simon M, Laura M, Vincent M, Tamas D, et al. Identification of grapevine microRNAs and their targets using high-throughput sequencing and degradome analysis. Plant J. 2010;62:960–76.
Rajagopalan R, Vaucheret H, Trejo J, Bartel DP. A diverse and evolutionarily fluid set of microRNAs in Arabidopsis thaliana. Genes Dev. 2006;20:3407–25.
Moxon S, Jing R, Szittya G, Schwach F, Rusholme PRL, Moulton V, et al. Deep sequencing of tomato short RNAs identifies microRNAs targeting genes involved in fruit ripening. Genome Res. 2008;18:1602–9.
Jinghua Y, Xunyan L, Baochen X, Na Z, Xiaodong Y, Mingfang Z. Identification of miRNAs and their targets using high-throughput sequencing and degradome analysis in cytoplasmic male-sterile and its maintainer fertile lines of brassica juncea. BMC Genomics. 2013;4:9.
Simon A, Wolfgang H. Differential expression analysis for sequence count data. Genome Biol. 2010;11:R106.
By Yoav B, Daniel Y. The control of the false discovery rate in multiple testing under dependency. Annals Statistics. 2001;29:1165–88.
Stéphane H, Benjamin P, Pascal F, Arnaud H, Nicolas R, Sébastien D. The poplar-poplar rust interaction: insights from genomics and transcriptomics. Journal of Pathogens. 2011;2011:1–11
Laurans F, Pilate G. Histological aspects of a hyper sensitive response in poplar to Melampsora larici-populina. Phytopathology. 1999;89:233–8.
Rinaldi C, Kohler A, Frey P, Duchaussoy F, Ningre N, Couloux A, et al. Transcript profiling of poplar leaves upon infection with compatible and incompatible strains of the foliar rust Melampsora larici-populina. Plant Physiol. 2007;144:347–66.
Hacquard S, Veneault-Fourrey C, Delaruelle C, Frey P, Martin F, Duplessis S. Validation of Melampsora larici-populina reference genes for in planta RT-quantitative PCR expression profiling during time-course infection of poplar leaves. Physiol Mol Plant Pathol. 2011;75:106–12.
Glazebrook J. Genes controlling expression of defense responses in Arabidopsis - 2001 status. Curr Opin Plant Biol. 2001;4:301–8.
Nawrath C, Metraux JP. Salicylic acid induction-deficient mutants of Arabidopsis express PR-2 and PR-5 and accumulate high levels of camalexin. Plant Cell. 1999;11:1393–404.
Navarro L, Dunoyer P, Jay F, Arnold B, Dharmasiri N, Estelle M, et al. A plant miRNA contributes to antibacterial resistance by repressing auxin signaling. Science. 2006;312:436–9.
Monk CE, Hutvagner G, Arthur JSC. Regulation of miRNA transcription in macrophages in response to Candida albicans. PLoS One. 2010;10:e13669.
Liu PT, Wheelwright M, Teles R, Komisopoulou E, Edfeldt K, Ferguson B, et al. MicroRNA-21 targets the vitamin D-dependent antimicrobial pathway in leprosy. Nat Med. 2012;18:267–73.
Ma F, Xu S, Liu X, Zhang Q, Xu X, Liu M et al. The microRNA miR-29 controls innate and adaptive immune responses to intracellular bacterial infection by targeting interferon-gamma. Nat Immunol. 2011;12:861–9.
Om PG, Pradeep S, Raj KG, Indu S. Current status on role of miRNAs during plant fungus interaction. Physiol Mol Plant Pathol. 2014;85:1–7.
Jonathan DGJ, Jeffery LD. The plant immune system. Nature. 2006;444:323–9.
Dangl JL, Jones JDG. Plant pathogens and integrated defence responses to infection. Nature. 2001;411:826–33.
Annegret K, Cécile R, Sébastien D, Marie B, Danny G, Frédéric D, et al. Genome-wide identification of NBS resistance genes in Populus trichocarpa. Plant MolBiol. 2008;66:619–36.
Padubidri VS, Ho-Ming C, Kanu P, Donna MB, Bruno AS, David CB. A microRNA superfamily regulates nucleotide binding site-leucine-rich repeats and other mRNAs. Plant Cell. 2012;3:859–74.
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;12:e84390. doi:10.1371/journal.pone.0084390.
Glazebrook J. Contrasting mechanisms of defense against biotrophic and necrotrophic pathogens. Annu Rev Phytopathol. 2005;43:205–27.
Robert-Seilaniantz A, Navarro L, Bari R, Jones JDG. Pathological hormone imbalances. Curr Opin Plant Biol. 2007;10:372e9.
Grant MR, Jones JDG. Hormone (dis) hormonymoulds plant health and disease. Science. 2009;324:750e2.
Zhaorong M, Ceyda C, Michael JA. Arabidopsis lyrata small RNAs: transient MIRNA and small interfering RNA loci within the Arabidopsis genus. Plant Cell. 2010;22:1090–103.
Addo-Quaye C, Eshoo TW, Bartel DP, Axtell MJ. Endogenous siRNA and miRNA targets identified by sequencing of the Arabidopsis degradome. Curr Biol. 2008;18:758–62.
Addo-Quaye C, Miller W, Axtell MJ. CleaveLand: a pipeline for using degradome data to find cleaved small RNA targets. Bioinformatics. 2009;1:130–1.
Lu S, Sun YH. Chiang VL Stress-responsive microRNAs in Populus. Plant J. 2008;1:131–51.
Peng S, Dan L, Zhoujia Z, Weilun Y, Xinli X. Identification of drought-responsive and novel Populus trichocarpa microRNAs by high-throughput sequencing and their targets using degradome analysis. BMC Genomics. 2013;14:233.
We are grateful to Dr. Minghao Pei (Rothamsted Research, Harpenden, Hertfordshire, UK) for his helpful suggestions and critical reading of the manuscript.
This work was supported by the National Natural Science Foundation of China (31500524, 31570637), the Fundamental Research Funds for the Central Universities (2572015EY02, DL13BA01 and 2572015CA09), the National High Technology Research and Development Program of China (863 Program) and the National Natural Science Foundation of China (31300571).
The authors declare that they have no competing interests.
Conceived and designed the experiments: DL, FW and ZW. Performed the experiments: DL and FW. Analyzed the data: FW and DL. Contributed reagents/materials/analytical tools: FW, DL, LZ, CW, QC, CN, RZ, YL and BW. Wrote the paper: FW, DL and QC. All authors have read and approved the manuscript.
Danlei Li and Feng Wang contributed equally to this work.
Unique miRNA classification. Text S2. Eleven features of miRNAs hairpin. Text S3. Five categories of the sliced-target transcripts. Text S4. Gene Ontology (GO) analysis. Text S5. Principles of the miRNA selected for RT-qPCR. Text S6. The 5 expression patterns of miRNA volatility change related to the development of rust. Text S7. Methods of + rust degradome cDNA libraries construction. Text S8. Detailed methods of RNAseq libraries construction. Text S9. The functions of CleaveLand 3.0 used in this study. Text S10. Procedure of a modified global normalization. Figure S1. Pie plots of data filtering and database mapping of the sRNA libraries from P. nigra × deltoids infected (+rust) and uninfected (−rust) with rust fungi M. larici-populina. Figure S2. Length distribution of mappable reads in sRNA libraries from P. nigra × deltoids infected (+rust) and uninfected (−rust) with rust fungi M. larici-populina. Figure S3. Predicted secondary structures of potential novel miRNAs from P. nigra × deltoids. Sequences indicated in red correspond to predicted miRNA. Figure S4. Targets of rust-responsive miRNAs among rust-regulated expressed genes. Figure S5. miRNA expression level validated by RT-qPCR at different time points post inoculation. Figure S6. Target gene expression level validated by RT-qPCR at different time points post inoculation. Figure S7. KEGG pathway and atlas. Table S1. A summary of standard data analysis results of the sRNA libraries from P. nigra × deltoids ‘Robusta’ infected (+rust) and uninfected (−rust) with rust fungi M. larici-populina. Table S2. Data summary of degradome library. Table S3. Transcription factors and their regulating miRNA. Table S4. post-transcriptional analyses of miRNA and resistance genes response to the infection of rust in ‘Robusta’. Table S5. The Gene Ontology of the rust-responsive gene and their regulating miRNA. Table S6. miRNA and target genes expression level and symptoms at different time points post inoculation. Table S7. Poplar defense pathway miRNA and target genes. Table S8. Primers used in the RT-qPCR. (DOC 3467 kb)
About this article
Cite this article
Li, D., Wang, F., Wang, C. et al. MicroRNA-mediated susceptible poplar gene expression regulation associated with the infection of virulent Melampsora larici-populina . BMC Genomics 17, 59 (2016). https://doi.org/10.1186/s12864-015-2286-6
- Melampsora larici-populina