RNA-Seq of Early-Infected Poplar Leaves by the Rust Pathogen Melampsora larici-populina Uncovers PtSultr3;5, a Fungal-Induced Host Sulfate Transporter

Abstract

Biotroph pathogens establish intimate interactions with their hosts that are conditioned by the successful secretion of effectors in infected tissues and subsequent manipulation of host physiology. The identification of early-expressed pathogen effectors and early-modulated host functions is currently a major goal to understand the molecular basis of biotrophy. Here, we report the 454-pyrosequencing transcriptome analysis of early stages of poplar leaf colonization by the rust fungus Melampsora larici-populina. Among the 841,301 reads considered for analysis, 616,879 and 649 were successfully mapped to Populus trichocarpa and M. larici-populina genome sequences, respectively. From a methodological aspect, these results indicate that this single approach is not appropriate to saturate poplar transcriptome and to follow transcript accumulation of the pathogen. We identified 19 pathogen transcripts encoding early-expressed small-secreted proteins representing candidate effectors of interest for forthcoming studies. Poplar RNA-Seq data were validated by oligoarrays and quantitatively analysed, which revealed a highly stable transcriptome with a single transcript encoding a sulfate transporter (herein named PtSultr3;5, POPTR_0006s16150) showing a dramatic increase upon colonization by either virulent or avirulent M. larici-populina strains. Perspectives connecting host sulfate transport and biotrophic lifestyle are discussed.

About the Authors:

Benjamin Petre

Affiliation: Unité Mixte de Recherche 1136 ‘Interactions Arbres/Microorganismes’, INRA (Institut National de la Recherche Agronomique)/Université de Lorraine, Centre INRA de Nancy, Champenoux, France

Emmanuelle Morin

Contributed equally to this work with: Emmanuelle Morin, Emilie Tisserant

Affiliation: Unité Mixte de Recherche 1136 ‘Interactions Arbres/Microorganismes’, INRA (Institut National de la Recherche Agronomique)/Université de Lorraine, Centre INRA de Nancy, Champenoux, France

Emilie Tisserant

Contributed equally to this work with: Emmanuelle Morin, Emilie Tisserant

Affiliation: Unité Mixte de Recherche 1136 ‘Interactions Arbres/Microorganismes’, INRA (Institut National de la Recherche Agronomique)/Université de Lorraine, Centre INRA de Nancy, Champenoux, France

Stéphane Hacquard

Affiliation: Unité Mixte de Recherche 1136 ‘Interactions Arbres/Microorganismes’, INRA (Institut National de la Recherche Agronomique)/Université de Lorraine, Centre INRA de Nancy, Champenoux, France

Corinne Da Silva

Affiliation: CEA-Genoscope, Centre National de Séquençage, Evry, France

Julie Poulain

Affiliation: CEA-Genoscope, Centre National de Séquençage, Evry, France

Christine Delaruelle

Affiliation: Unité Mixte de Recherche 1136 ‘Interactions Arbres/Microorganismes’, INRA (Institut National de la Recherche Agronomique)/Université de Lorraine, Centre INRA de Nancy, Champenoux, France

Francis Martin

Affiliation: Unité Mixte de Recherche 1136 ‘Interactions Arbres/Microorganismes’, INRA (Institut National de la Recherche Agronomique)/Université de Lorraine, Centre INRA de Nancy, Champenoux, France

Nicolas Rouhier

Affiliation: Unité Mixte de Recherche 1136 ‘Interactions Arbres/Microorganismes’, INRA (Institut National de la Recherche Agronomique)/Université de Lorraine, Centre INRA de Nancy, Champenoux, France

Annegret Kohler

Affiliation: Unité Mixte de Recherche 1136 ‘Interactions Arbres/Microorganismes’, INRA (Institut National de la Recherche Agronomique)/Université de Lorraine, Centre INRA de Nancy, Champenoux, France

Sébastien Duplessis

* E-mail: [email protected]

Affiliation: Unité Mixte de Recherche 1136 ‘Interactions Arbres/Microorganismes’, INRA (Institut National de la Recherche Agronomique)/Université de Lorraine, Centre INRA de Nancy, Champenoux, France

Introduction

Parasitic microbes continuously challenge plants but constitutive and inducible host immunity is generally sufficient to confer broad-spectrum resistance to non-adapted pathogens [1]. A successful infection requires that microbes secrete effectors into infected tissues to repress and manipulate host defense and physiology, resulting in the promotion of biotrophic growth and the proper setting of compatible interaction [2]. In turn, plants possess several hundreds of resistance (R) proteins that, upon direct or indirect recognition of pathogen effectors, trigger strong defense responses and lead to incompatibility [3]. The issue of an interaction as compatible or incompatible is determined within hours following the formation of the first structures of infection and subsequent effector secretion, even after the very first physical contacts [4]. At these early stages, the ability of pathogen effectors to rapidly and efficiently manipulate host functions and simultaneously escape R protein recognition is though to be the key of compatibility [5]. Identification of both effectors and pathogen-modulated host-functions is currently a major goal to unravel molecular basis of plant-microbe interactions [6]. To date, most fungal effectors identified are lineage-specific small-secreted proteins (SSP) of unknown function and hundreds of such SSP coding-genes (i.e. candidate effectors) have been identified in fungal genomes [7]–[11]. Host functions targeted by bacterial effectors are mostly restricted to the plant immune system signaling network, vesicle trafficking, protein degradation systems and transcription of specific transporters [12]. Only a handful of eukaryote effector targets are currently known but the larger effector candidate repertoires compared with those of prokaryotes suggests specialized host manipulations likely linked to nutrient acquisition strategies, infection structures (e.g. haustoria) formation and host reprogramming [13]–[15].

Genome sequences of Populus trichocarpa and Melampsora larici-populina are available, making the poplar-poplar rust interaction an emerging model for post-genomic studies of tree immunity and fungal biotrophy [16]–[17]. Several genome-wide analyses uncovered the huge expansion and diversification of defense- or virulence-related gene families in poplar and the poplar rust fungus, respectively (for a recent review, see [18]). Several transcriptomic studies identified defense response genes in poplar upon infection by Melampsora spp. between 2 and 10 days after inoculation during compatible or incompatible interactions (for a review, see [19]). M. larici-populina transcriptome analyses were also conducted during a compatible time-course infection as well as in specific infection structures captured by laser microdissection at late sporulating stages [11], [20]–[21]. A major output of these studies is the timely coordinated and structure-specific expression of hundreds of genes encoding SSPs all along the infection process, suggesting an early and sustained host manipulation [18]. Contrasting with the flow of data gained on both the host and the pathogen transcriptomes at late stages of infection, our knowledge of the first two days of infection remains limited, whereas they encompass a crucial step in setting the outcome of the interaction.

Next-generation sequencing technologies have been an important breakthrough for sensitive, quantitative, annotation-independent and high-throughput transcriptome analyses referred to as RNA-Seq [22]–[24]. Recent works reported the use of RNA-Seq (Illumina or 454-pyrosequencing) to get insights into the transcriptomes of plant-microbe associations [25]–[27]. Using original or phylogenetically related reference genomes, the authors were able to discriminate microbial and host transcripts. In the poplar-poplar rust pathosystem, the availability of both well-annotated genome sequences would greatly facilitate this kind of approach.

In the present study, we used 454-pyrosequencing to analyse the transcriptome of compatible and incompatible poplar-poplar rust interactions within the first two days of infection. Methodological aspects such as discrimination of reads origin (poplar vs. poplar rust) and poplar RNA-Seq data validation by oligoarrays were considered. This work identified five early-expressed candidate effectors from M. larici-populina as well as an early-induced sulfate transporter from poplar (PtSultr3;5). These data complete our knowledge of the poplar-poplar rust interaction and will help prioritizing functional studies in this pathosystem.

Results and Discussion

454-pyrosequencing Transcriptome of Poplar Leaves Infected by M. larici-populina

A time-course infection of poplar leaf disks inoculated with either the virulent strain 98AG31 (compatible interaction) or the avirulent strain 93ID6 (incompatible interaction) of M. larici-populina has been established and tissues were harvested at 18, 21 and 24 hours post-inoculation (hpi) for the incompatible interaction (I₁₈-I₂₁-I₂₄) and 18, 24 and 48 hpi for the compatible interaction (C₁₈-C₂₄-C₄₈) (see Fig. S1 for a summary of all experiments and analyses presented in this study). M. larici-populina haustoria are formed in both compatible and incompatible interactions as soon as 17 hpi [28], and the first poplar defense responses are observed only in incompatible interactions at 48 hpi [29]. Studies conducted on the flax rust and the bean rust model systems showed that haustoria are key structures for rust fungi to secrete effectors into their host [30]–[33]. Thus, our experiment was designed to identify early-expressed putative rust effectors (i.e. after haustoria formation, between 18 and 24 hpi) as well as early-regulated poplar genes not previously reported. A later time-point in the compatible interaction (C₄₈) was included to define the dynamic of SSP expression and to detect modulation in the host transcriptome. For each time-point, cDNA were sequenced by 454-pyrosequencing (GS-FLX titanium), generating a total of 841,301 high-quality nucleotide reads (∼304 Mb) with an average length of 356 bp (Table 1, Fig. S2). For the six conditions, we independently assembled reads into contigs. As presented in Table 1, a total of 713,505 reads were assembled into 90,398 contigs. The average contig size is 659 pb whereas the average number of read per contigs is around 8 (Fig. S2). Homology searches against both P. trichocarpa and M. larici-populina genome sequences attributed more than 97% of the assembled reads (694,138) to poplar (693,489) or the poplar rust (649) genes (Table 1). This approach will hereafter be referred to as ‘contig-blast method’. We used a second method for read processing, which consisted in the mapping of all the reads onto the reference genomes and current available annotations after reassembly using the Program to Assemble Spliced Alignments (PASA; [34]) (herein called ‘read-mapping method’). This method identified 616,879 and 359 reads corresponding to P. trichocarpa and M. larici-populina genes, respectively, which is comparable to the numbers gained by the contig-blast method (Table 1). Moreover, expression levels of poplar transcripts (i.e. number of reads per transcripts) estimated by both approaches correlated well, indicating that these methods provide similar outputs (Fig. S3). The data processed through the read-mapping method were more stringent and thus retained for the quantitative analysis of poplar transcript levels, whereas the contig-blast method was retained for the analysis of M. larici-populina transcripts.

[Figure omitted. See PDF.]

Table 1. Summary of 454-pyrosequencing transcriptome data from poplar leaves infected by M. larici-populina.

https://doi.org/10.1371/journal.pone.0044408.t001

M. larici-populina Early-expressed Transcripts and Candidate Effectors

A total of 649 reads were aligned to 361 M. larici-populina genes by the contig-blast method, whereas 359 reads were mapped to 280 M. larici-populina genes by the read-mapping method (Table 1). In total, 97% of the genes identified by the read-mapping method were also retrieved by the contig-blast method. The lower number of genes identified by the read-mapping method are due to intrinsic differences in assembly parameters. The C₄₈ condition itself accounted for almost half of the reads whatever method used, which is in accordance with the fungal growth previously reported at this stage of the infection [29], [35]. The low number of reads does not allow proper quantitative comparison of fungal transcripts abundance between time-points. Thus, we focused our attention on the 40 genes cumulating a total of three reads or more with any of the two methods (Table 2). Among those, a total of 21 genes encode secreted proteins, of which 19 correspond to small secreted proteins of less than 300 amino acids, previously reported as putative candidate effectors of the poplar rust fungus [11]. Despite a very low level of fungal transcripts detected in infected leaf samples, the high number of transcripts encoding putative effectors is striking and in accordance with the levels reported in purified haustoria of other fungal biotroph [32], [36] or in M. larici-populina biotrophy-related structures microdissected from infected poplar tissues [20]. Noteworthy, Mlp-124067 has been reported among the most highly expressed fungal transcript at 24 hpi during the compatible poplar-poplar rust interaction [21]. The most expressed transcript (Mlp-94736, 92 reads) encodes a SSP of 158 amino acids with no functional annotation and only a few homologs in databases (Table 2). Strinkingly, expression of the transcript is mostly explained by the I₁₈ time-point that cumulates 86 reads. Mlp-94736 belongs to a small cluster of genes within 40 kb on the scaffold 72 of which 4 genes (Mlp-94735, Mlp-94736, Mlp-94740 and Mlp-94741) were expressed in the dataset. The second most expressed transcript (Mlp-95625, 36 reads) encodes a SSP of 89 amino acids that possesses a short inositol phosphatase domain (conserved domain cd01640). Recent studies have reported that entry of oomycete and fungal effectors into plant cells is mediated by binding to phosphatidylinositol phosphates (PIP) present on the extracellular face of the host plasmalemma [37]–[39]. The early expression of a secreted inositol phosphatase by the poplar rust fungus is puzzling and it is tempting to speculate a role related to the PIP-binding of effectors. Several genes showing homology to ribosomal proteins and elongation factors were also identified (Table 2), which is in accordance with the profile of haustoria expressed transcripts of the obligate biotrophic fungus Golovinomyces orontii [36]. Noticeable differences in reads numbers could be observed between the two read assignment methods, particularly for the genes cumulating the largest numbers of reads with the contig-blast method (Table 2). We manually validated the alignment of these contigs onto the M. larici-populina genome and their absence with the second method is due to differences in alignment parameters.

[Figure omitted. See PDF.]

Table 2. List of selected M. larici-populina (Mlp) genes expressed at early stages of poplar leaf infection and cumulating three reads for all time-points assessed in the study with any of the read assignment method.

https://doi.org/10.1371/journal.pone.0044408.t002

Expression of M. larici-populina transcripts during poplar leaf colonization (compatible interaction) and urediniospore germination was previously assessed by the mean of Genome oligoarrays [21]. Among the set of 40 fungal genes, 30 were represented on the oligoarray, of which 27 are expressed in the situations tested and particularly, 24 are expressed at 24 and 48 hpi, which provides a nice support for the transcripts detected by the 454-pyrosequencing approach. Overall, the low number of fungal transcripts detected in our study reflects the limited quantity of fungal structures in early-infected leaves [29], [35]. In comparison, Fernandez and collaborators [27] identified almost one fungal for two host transcripts in the coffee-coffee rust fungus interaction, however this study targeted a late stage of infection when the fungus has completely invaded host tissues. Although providing valuable insights into early-expressed fungal transcripts, our analysis remains scale-limited and argues for the use of other NGS technologies such as Illumina for a deeper sampling of the pathogen transcriptome.

454-pyrosequencing does not Saturate the Poplar Transcriptome

A total of 616,879 reads were mapped onto 21,311 poplar transcripts (i.e. predicted protein coding genes in poplar genome), considered as expressed when at least one read was mapped (Table 1). Among these 21,311 transcripts only 6,755 were expressed in all conditions. The average and median of total reads number per expressed transcript are 28.95 and 8 respectively, ranging from 1 to 6,513 (Fig. S4A). Saturation curves indicate that 52±5 new expressed transcripts are detected for each 1000 supplemental reads after 100,000 (Fig. S4B). This demonstrates that the poplar transcriptome is not saturated and consequently unsuitable for whole-genome quantitative expression analysis. As specified recently by Malone and Oliver [40], an important concern in RNA-Seq experiments is the depth of sequencing required to efficiently sample the transcriptome. By using Illumina sequencing, Graveley and collaborators [41] reported that even 50 millions of mapped reads were not enough to fully saturate the fly-head transcriptome. Despite huge progresses in sequencing, it seems that full saturation of transcriptomes still cannot be considered as an easy purpose.

In order to isolate a subset of transcript with a sufficient coverage to enable quantitative analysis of expression, we dramatically reduced our set from 21,311 to 1,909 highly expressed (HE) transcripts, with a minimal arbitrary cut-off of 60 total reads per HE transcript (Table 1, Fig. S4).

Poplar RNA-Seq Expression Levels are Validated by Oligoarrays and RT-qPCR

In order to validate poplar RNA-Seq data, we performed whole-genome poplar oligoarrays and compared expression values gained by both methods. The DNA-chip approach detected the expression of 34,964 transcripts (76% of the coding genome) of which 31,644 were expressed in all conditions (Fig. 1a). Compared with the 21,311 expressed transcripts revealed by RNA-Seq, oligoarrays clearly demonstrate that the RNA-Seq approach was not deep enough to cover the entire poplar transcriptome. Oligoarrays confirmed the expression of 20,105 of the 21,311 transcripts detected as expressed by RNA-Seq (Fig. 1a). Among the 1,206 unconfirmed transcripts, 440 were detected below the oligoarrays background threshold whereas 766 were not spotted on the chip. Otherwise, 1,794 of the 1,909 HE transcripts were detected as highly-expressed on oligoarrays and are hereafter referred to as validated highly expressed (VHE) transcripts. Among the 115 non-validated HE transcripts, only 5 were detected below the threshold whereas 110 were not included in the chip design. The expression levels of the 1,794 VHE transcripts assessed by oligoarray and RNA-Seq showed a good correlation (r² = 0.37) that increases with expression level (Fig. 1b). To further validate the RNA-Seq expression profiles, we tested 12 selected VHE transcripts by RT-qPCR. The expression profiles were in agreement with those obtained with the two other methods (Fig. S5). Otherwise, it is noteworthy to mention that 20% of the most highly expressed transcripts on oligoarrays were not supported by RNA-Seq, revealing a unidirectional and non-negligeable technical biase between these approaches discussed in Text S1 (Fig. S6).

[Figure omitted. See PDF.]

Figure 1. Ln-transformed expression levels of poplar transcripts.

(A) Venn-diagram of expressed poplar transcripts measured by whole-genome oligoarray and RNA-Seq. Numbers in brackets correspond to the percentage of the coding genome; numbers in parentheses correspond to the highly expressed (HE) transcripts and numbers in parentheses and in bold correspond to HE transcripts validated by oligoarrays (i.e. VHE transcripts). (B) RNA-Seq and oligoarrays average expression levels correlation. The black rectangle indicates the area containing the 1,794 VHE transcripts discussed in the text.

https://doi.org/10.1371/journal.pone.0044408.g001

Transcriptome of Poplar Leaves: Overall Stability at Early Stages of Infection

In order to identify early-regulated transcripts in the poplar transcriptome, expression levels were compared pair-by-pair between the 6 different conditions tested (Fig. 2). Among the VHE transcripts, 1,793 showed overall stable RNA-Seq expression levels, with average relative standard deviation percentage (%RSD, see methods) values of 35, ranging from 2 to 101 (Fig. S7a). The only highly variable transcript with a %RSD of 186 is a predicted sulfate transporter homologous to Arabidopsis thaliana Sultr3;5 (AtSultr3;5) [42] and hybrid poplar Populus tremula x Populus alba Sultr3;5 (PtaSultr3;5) [43] and thus was herein named P. trichocarpa Sultr3;5 (PtSultr3;5). The overall stability of VHE transcript expression levels was confirmed by oligoarrays, which present %RSD values from 1 to 62 (Fig. S7b). Moreover, expression levels of the 34,964 oligoarray-expressed transcripts were equally stable, with only 18 transcripts with %RSD above 100, which all correspond to low and barely expressed genes (Fig. 2, Table S1). Functional classification of the 1,794 VHE transcripts as well as detailed analysis of the stress responsive phenylpropanoid pathway further confirmed the absence of defense responses (Text S2, Fig. S8 and S9). Hence, except PtSultr3;5, the leaf transcriptome remains extremely stable for the time-points investigated.

[Figure omitted. See PDF.]

Figure 2. RNA-Seq and oligoarray two-by-two expression levels correlations along time-course infection of poplar leaves by M. larici-populina.

Ln-transformed expression levels of the 1,794 validated highly expressed (VHE) transcripts are plotted for RNA-Seq on the up-right, and expression levels of the 34,964 transcripts detected on oligoarrays are plotted on the bottom-left. Arrows indicate values corresponding to the sulfate transporter PtSultr3;5 (Poptr_0006s16150) transcript discussed in the text.

https://doi.org/10.1371/journal.pone.0044408.g002

Many host transcriptomes have been reported to be highly dynamic during the first events of plant-microbe interactions [44]–[45], as highlighed recently in soybean leaves as early as 12 hours after inoculation by the rust fungus Phakospora pachyrhizi [46]. Nevertheless, in the poplar-poplar rust system, typical defense-response genes encoding pathogenesis-related (PR) proteins such as PR1 or PR5, glutathion S-transferase 18 (GST18) or the rust-induced secreted protein (RISP) are not regulated and present a basal expression at early stages of leaf infection (Table S1, [29], [47]). Indeed, the first defense responses are reported at 48 hpi in the incompatible interaction, concomitant with an arrest of M. larici-populina avirulent strain progression in leaves [28]–[29]. Considering the global stability of the poplar trancriptome and the fact that sulfate transporters are not reported as typical defense-related genes in the literature, we conclude that no defense responses could be detected in early rust-infected poplar leaves.

A Poplar Sulfate Transporter is Specifically Induced by the Rust Fungus

The sulfate transporter PtSultr3;5 displayed the highest %RSD among all VHE transcripts (Fig. S7). Interestingly, the variability observed is explained by the condition C₄₈ (515 reads) with an almost 20-fold up-regulation compared with the other five conditions that remain stable (26±15 reads) (Table S1). The different methods used in this study to quantify gene expression provided highly similar profiles for this gene, validating its effective regulation (Fig. 3a). Moreover, complementary infection time-points included in the RT-qPCR assay, as well as comparison to mock-inoculated control leaves, provided a more accurate expression pattern (Fig. 3b). A strong induction of PtSultr3;5 was observed at 48 hpi in both compatible and incompatible poplar-polar rust interactions compared with all other conditions that showed an overall stable expression. Thus, PtSultr3;5 presents a fungal-specific strain-independent induction in poplar leaves. Since we did not detect neither defense responses nor large-scale transcriptome regulations, and considering the specificity and amplitude of the induction of PtSultr3;5, we hypothesized that this specific up-regulation could result from a pathogen manipulation.

[Figure omitted. See PDF.]

Figure 3. Expression profile of PtSultr3;5.

(A) Expression pattern of PtSultr3;5 (Poptr_0006s16150) assessed by oligoarray, RT-qPCR and RNA-Seq. RNA-Seq and RT-qPCR values are Ln-transformed. For RNA-Seq, both results from the contig-blast (cb) and the read-mapping (rm) methods are presented. Values are in arbitrary units. (B) Detailed RT-qPCR profile of PtSultr3;5 between 0 and 48 hours post-inoculation (hpi). Expression values are normalized to the ubiquitin reference gene expression (see methods).

https://doi.org/10.1371/journal.pone.0044408.g003

PtSultr3;5 is the ortholog of AtSultr3;5 (At5g19600) that belongs to the family of low-affinity sulfate transporter and is localized at plasma membranes [48]. Interestingly, PtSultr3;5 is also the ortholog of the Lotus japonicus symbiotic sulfate transporter 1 (SST1) with 57% of identity (Fig. S10). In L. japonicus intracellular rhizobia symbiotic interactions, SST1 is located onto the symbiosome membrane, the interface for host-symbiont nutrient exchanges [49]. Moreover, Krusell and colleagues [50] demonstrated that sst1 expression is nodule-specific and that a mutant lacking sst1 was defective for symbiosis. The authors concluded that SST1 is essential for symbiosis, likely by translocating sulfate from the host cell to the symbiont. Interestingly, Felten and collaborators [51] reported that PtSultr3;5 is one of the 15 most induced transcripts in poplar roots during early interaction with the symbiotic fungus Laccaria bicolor. In these experiments, the host and its symbiont where not in contact, suggesting that indirect and diffusible signals were directly or indirectly responsible of host-transcriptome reprogramming.

No obvious link has been established so far between host sulfur metabolism and mechanisms of pathogenicity in biotrophs. However, in a recent groundbreaking study, Mukhtar and collaborators [52] showed that multiple effectors from divergent biotroph pathogens target response to low sulfur proteins 1 and 2 (LSU1 and LSU2). As mentioned by the authors, the identification of sulfur-related proteins as conserved targets of effectors is puzzling since one of the hallmarks of biotrophy is the loss of key genes necessary to assimilate sulfur. Indeed, this has been reported for several biotroph pathogens belonging to distinct phyla such as the oomycete Hyaloperonospora arabidopsidis, the ascomycete Blumeria gaminis and the basidiomycete Puccinia graminis f. sp. tritici [9], [17], [53]–[54]. In the case of M. larici-populina, genes required to perform the primary sulfate assimilation were nevertheless identified, although a transketolase domain of the beta-subunit of the sulfide reductase is missing [17]. Interestingly, 2 out of 3 M. larici-populina sulfate transporter transcripts are expressed during poplar leaf infection and in urediniospores, however sulfate metabolism transcripts are barely detected [21].

It is tempting to speculate that M. larici-populina could perform a targeted-manipulation of poplar physiology by inducing PtSultr3;5 to derivate sulfate via a host-encoded mechanism. Such a manipulation of host transporter expression by pathogens has been described in interactions between Xanthomonas spp. and rice or pepper, where host genes encoding sugar transporters are specifically and necessarily induced by transcription activator-like (TAL) bacterial effectors to settle compatibility [55]–[57]. Interestingly, fungal pathogens also induce sugar exporters during plant infection [57] and although effectors have been brought up, the molecular mechanism underpining this phenomenon remains to be determined [58]. To our knowledge, no TAL effectors have been reported outside of the bacterial genus Xanthomonas and Ralstonia. Interestingly, a recent report by Plett and colleagues [14] has shown that a small-secreted effector protein from the symbiotic fungus L. bicolor targets the host nucleus to settle compatibility through poplar roots transcriptome modulation. Moreover, in plant-oomycete interactions, the massive delivery of effector proteins into the host nucleus has been observed [15], [59]. Altogether, these studies raise interesting perspectives on the ability of biotrophic eukaryote plant symbionts to perform targeted modulations of host genetic programs.

Conclusions

In this study, we performed a transcriptome analysis of early stages of poplar leaf infection by the rust fungus M. larici-populina both during compatible and incompatible interactions. We covered critical steps of the fungus biotrophic growth in the host, when the first haustorial structures are formed ahead of typical host defense responses. Complementary large-scale expression approaches were used including 454 pyrosequencing and oligoarrays to assess the accuracy of the expression data. At these early stages, only a few fungal transcripts were detected in planta indicating that 454-pyrosequencing could not support investigation of fungal function. Nevertheless, early-expressed candidate effectors corresponding to probable rust pathogenicity determinants were detected and will be targeted in future investigations. We report the absence of induction of typical defense response genes in infected poplar leaves before 48 hpi, which is contrasted by a highly induced poplar sulfate transporter. This last observation, in the light of our knowledge of biotrophy, is puzzling and establishes a new correlative link between biotrophic growth and host sulfate. PtSultr3;5 is also early induced in poplar roots upon indirect interaction with a symbiotic fungus (i.e. pre-symbiotic contact) and could represent a conserved mechanism targeted by symbiotic and pathogenic biotrophs to settle compatibility with poplar. Further functional investigations are now needed to understand the role of PtSultr3;5 and more generally to underpin how plant sulfate transport could play a role in plant-microbe interactions.

Materials and Methods

Biological Material, Infection Procedures and Time-course Infection

Melampsora larici-populina isolates 98AG31 (pathotype 3-4-7) and 93ID6 (pathotype 3-4), respectively virulent and avirulent on the hybrid P. trichocarpa x P. deltoides poplar cultivar ‘Beaupré’, were used in this study. Urediniospores multiplication and poplar leaf inoculation procedures were carried out as previously described [29] except that inoculations were performed on 5 cm² leaf disks as reported in other studies [60]–[62]. For time-course infection analyses, leaf disks were harvested at 0, 12, 15, 18, 21, 24 and 48 hpi in compatible (C) and incompatible (I) interactions as well as in mock-inoculated (water agar) condition (T). The 6 following conditions were used for 454-pyrosequencing and oligoarrays: I₁₈, I₂₁, I₂₄, C₁₈, C₂₄, C₄₈ (Fig. S1). All of these conditions were used for RT-qPCR analyses. No biological duplication was performed.

RNA Extraction and cDNA Preparation

Total RNA extraction and quantity/quality control were performed as previously described [29]. For 454-pyrosequencing and oligoarray experiments, tagged double-strands (ds) cDNA were synthesized with the SMARTer PCR cDNA Synthesis Kit (Clontech, Saint-Germain-en-Laye, France), following manufacturer instructions, with 0.5 µg of starting total RNA and 19 cycles of PCR amplification in 25 µL. PCR products were purified with the Quiaquick PCR Purification Kit (Qiagen, Couteboeuf, France), quantified with a nanodrop 1000 spectrophotometer (Labtech, Palaiseau, France) and their quality was assessed by electrophoretic RNA profiling with an Experion analyzer using the Experion RNA Standard-Sens analysis kit (Bio-Rad, Marnes-la-Coquette, France). For RT-qPCR experiments, cDNA were synthesized using the iScript cDNA Synthesis Kit (Bio-Rad,) with 1 µg of total RNA and following manufacturer instructions.

454-pyrosequencing and Oligoarrays

Sequences were generated from 5 µg of full-length cDNA using the 454-GS-FLEX Titanium pyrosequencing technology on a Roche 454 Titanium sequencer (single read strategy) at the CEA-Genoscope (Centre National de Séquençage, Evry, France; http://www.genoscope.cns.fr/spip) following standard procedures recommended by Roche. Two half-plates were used: one for the I₁₈, I₂₁, I₂₄ conditions mixed 1/3-1/3-1/3 together with specific tags; and another one for the C₁₈, C₂₄, C₄₈ conditions mixed 1/3-1/3-1/3 together with specific tags. Raw 454-sequences are available upon request. The interest and use of 454-sequences for complementary annotation of poplar and poplar rust genomes are discussed in Text S3. Whole-genome poplar oligoarrays were generated and hybridized at the NimbleGen facilities (NimbleGen systems, Reykjavik, Iceland) following standard protocols [29], [51]. Microarray probe intensities were quantile normalized across chips, using the Arraystar software (DNAStar, Inc., Madison, Wi, USA). The background level was fixed at 150 (average 3-fold random probe intensities across oligoarray hybridizations). The complete oligoarray expression datasets are available at the NCBI GEO website as series number GSE34802.

RT-qPCR

Specific primers were designed for targeted ‘Beaupré’ transcripts and RT-qPCR were carried out as previously detailed on technical duplicates [20]. Sequences gained by RNA-Seq were used for primer design in order to avoid polymorphism bias during alleles amplifications of the ‘Beaupré’ hybrid cultivar. Transcript expression was normalized to the poplar ubiquitin reference transcript [29] and calculated using the following equation: Relative Expression  =  [(target gene primer efficiency)^−Ct]/[(reference gene primer efficiency)^−Ct], where Ct is the number of PCR amplification cycles necessary for signal detection [63]. Primer sequences are presented in Table S2.

Bioinformatic Procedures for RNA-Seq Data Analysis

Raw reads were processed to remove SmartER adaptor sequences and reads shorter than 70 pb were discarded. For the contig-blast method, reads were assembled into contigs with the Mimicking Intelligent Read Assembly (MIRA) program [64] and blast homology searches were performed against both M. larici-populina (version 1.0, JGI website http://genome.jgi-psf.org/Mellp1/) and P. trichocarpa (version 2.2, Phytozome website http://www.phytozome.net/poplar) genome sequences. Contigs were assigned to the reference genomes based on blast scores and e-values, and ambiguous assignments were eventually resolved by searching for plant or fungal homologs in international databases as previously described [27]. All assignments were validated by genome assignments with GenomeThreader [65]. For the read-mapping method, reads were aligned directly on P. trichocarpa and M. larici-populina genome sequences using Program to Assemble Spliced Alignments (PASA, [34], http://pasa.sourceforge.net/) with default parameters and GMAP as the program to align transcripts and annotated coding sequences available at the P. trichocarpa and M. larici-populina genomes websites. PASA included de novo assembly on the genome considering available annotations and proposed predictions for gene correction, gene creation as well as identification of putative splicing variants (see Text S3). Numbers of reads attributed to annotated genes at the six time-points were retrieved using in-house scripts and reads’ numbers were normalized between the six conditions prior further quantitative comparisons (Table 1). Expression analyses and graphics were done with Excel for Mac (Microsoft) or R version 2.9.2 (http://cran.r-project.org/). For a given transcript, %RSD of the six expression values was calculated according to the following equation: %RSD  =  [(standard deviation/average) x 100]. Sequences were analyzed in Text Wrangler version 2.3 (Bare Bones Software, Inc.) and Tablet version 1.10.05.21 [66]. Alignments were performed with the Multalin online software (http://multalin.toulouse.inra.fr/multalin/). The Functional Catalogue (FunCat) version 2.1 was used for classification of genes with slight modifications within categories [67]. Read and contig sequences as well as PASA outputs (i.e. gene model re-annotation as well as identification of alternative splicing and of allelic polymorphism) are available upon request. Saturation curves were performed with Analytic Rarefaction software version 2.0 (http://www.huntmountainsoftware.com/).

Supporting Information

[Figure omitted. See PDF.]

Figure S1.

Experimental design and bioinformatic procedure summary. (A) Experimental design of the time-course infection. All conditions and time-points were used for standard cDNA synthesis and subsequent RT-qPCR, whereas the 6 conditions highlighted in black (I₁₈, I₂₁, I₂₄, C₁₈, C₂₄, C₄₈) were used for hybridization on poplar oligoarrays and for tagged cDNA synthesis and subsequent 454-pyrosequencing as depicted in b. (B) Summary of the bioinformatic procedure and analyses performed in the study.

https://doi.org/10.1371/journal.pone.0044408.s001

(TIF)

Figure S2.

Read and contig length distribution.

https://doi.org/10.1371/journal.pone.0044408.s002

(TIF)

Figure S3.

RNA-Seq data analysis methods comparison. RNA-Seq mean expression levels are Ln-transformed.

https://doi.org/10.1371/journal.pone.0044408.s003

(TIF)

Figure S4.

RNA-Seq coverage of the poplar transcriptome. (A) Distribution of RNA-Seq reads among the 21,311 expressed poplar transcripts. The black rectangle highlights the subset of 1,909 highly expressed (HE) transcripts considered for quantitative analyses. (B) Average saturation curve of the poplar transcriptome by RNA-Seq.

https://doi.org/10.1371/journal.pone.0044408.s004

(TIF)

Figure S5.

Correlation between RNA-Seq, oligoarrays and RT-qPCR normalized expression levels for 12 selected validated highly expressed (VHE) transcripts (listed in Table S2).

https://doi.org/10.1371/journal.pone.0044408.s005

(TIF)

Figure S6.

Dot-plot of the 285 most-expressed transcripts on oligoarrays (average expression values above 20,000) and corresponding RNA-Seq expression values.

https://doi.org/10.1371/journal.pone.0044408.s006

(TIF)

Figure S7.

Percentage-relative standard deviation (%RSD) distributions. (A) %RSD for RNA-Seq expression levels of the 1,794 validated highly expressed (VHE) transcripts. (B) %RSD of oligoarray expression levels of the 1,794 VHE transcripts. Arrows in (A) and (B) indicate values of the fungal-induced sulfate transporter (PtSultr3;5, Poptr_0006s16150).

https://doi.org/10.1371/journal.pone.0044408.s007

(TIF)

Figure S8.

Functional Catalogue (FunCat) classification of the 1,794 validated highly expressed (VHE) transcripts. (A) Distribution of FunCat categories. Values represent percentage of the 1,794 VHE transcripts. (B) Average expression levels in FunCat categories. (C) Average percentage-relative standard deviation (%RSD) in FunCat categories.

https://doi.org/10.1371/journal.pone.0044408.s008

(TIF)

Figure S9.

Transcriptional levels of the phenylpropanoid pathway transcripts in poplar leaves infected by M. larici-populina .

https://doi.org/10.1371/journal.pone.0044408.s009

(TIF)

Figure S10.

Alignment of selected plant sulfate transporters homologous to PtSultr3;5. Amino acid sequences of PtSultr3;5 (Populus trichocarpa, POPTR_0006s16150.1), Beaupré-9 (contig C48-lrc9, P. trichocarpa x Populus deltoides ‘Beaupré’), Beaupré-124 (contig C48-lrc-124, P. trichocarpa x P. deltoides ‘Beaupré’), AtSultr3;5 (Arabidopsis thaliana, At5g19600), NCBI ID 117557144 (PtaSultr3;5, Populus tremula x Populus alba), NCBI ID 255549068 (Ricinus communis), NCBI ID 225445290 (Vitis vinifera), NCBI ID 45720463 (Brassica oleracea), NCBI ID 297812143 (Arabidopsis lyrata) and SST1 (Lotus japonicus). Asterisks indicate polymorphic residues between ‘Beaupré’ alleles and the double-head arrow marks additional predicted amino acids in the P. trichocarpa ‘Nisqually-1’ genome discussed in Test S3.

https://doi.org/10.1371/journal.pone.0044408.s010

(TIF)

Table S1.

Summary table of expression, annotation and sequence data used in the study for poplar genes.

https://doi.org/10.1371/journal.pone.0044408.s011

(ZIP)

Table S2.

List of primers used in this study.

https://doi.org/10.1371/journal.pone.0044408.s012

(XLS)

Text S1.

Uncorrelated oligoarray/RNA-Seq transcript expression.

https://doi.org/10.1371/journal.pone.0044408.s013

(DOC)

Text S2.

Functional classification of the 1,794 poplar validated highly expressed (VHE) genes.

https://doi.org/10.1371/journal.pone.0044408.s014

(DOC)

Text S3.

Interest and use of 454-sequences for complementary annotation of poplar and poplar rust genomes.

https://doi.org/10.1371/journal.pone.0044408.s015

(DOC)

Acknowledgments

The authors thank Patrice Vion, Jean-louis Churin and Olivier Cael (INRA Nancy) for technical support and advices on poplar culture in greenhouses; Pascal Frey and Bénédicte Fabre (INRA Nancy) for the access to the M. larici-populina strains collection and discussions on the poplar-poplar rust pathosystem; Aurore Coince (INRA Nancy) for advices on R graphical tools and saturation curves design; Diana Fernandez for fruitful discussions in the course of the collaboration on plant-rust fungi transcriptome analyses (collaborative CEA-Génoscope/INRA/IRD sequencing project).

Author Contributions

Conceived and designed the experiments: BP SD. Performed the experiments: BP SH CD CDS JP. Analyzed the data: BP EM ET SD. Contributed reagents/materials/analysis tools: FM AK. Wrote the paper: BP SH FM NR SD.

Citation: Petre B, Morin E, Tisserant E, Hacquard S, Da Silva C, Poulain J, et al. (2012) 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 7(8): e44408. https://doi.org/10.1371/journal.pone.0044408

References

1. Schulze-Lefert P, Panstruga R (2011) A molecular evolutionary concept connecting nonhost resistance, pathogen host range, and pathogen speciation. Trends Plant Sci 16: 117–125.

2. Hok S, Attard A, Keller H (2010) Getting the Most from the Host: How Pathogens Force Plants to Cooperate in Disease. Mol Plant Microbe Interact 23: 1253–1259.

3. Dodds PN, Rathjen JP (2010) Plant immunity: towards an integrated view of plant-pathogen interactions. Nat Rev Genet 11: 539–548.

4. Nirmala J, Draeder T, Chen X, Steffenson B, Kleinhofs A (2010) Stem rust spores elicit rapid RPG1 phosphorylation. Mol Plant Microbe Interact 23: 1635–1642.

5. Schneider DJ, Collmer A (2010) Studying Plant–Pathogen Interactions in the Genomics Era: Beyond Molecular Koch’s Postulates to Systems Biology. Annu Rev Phytopathol 48: 457–479.

6. Panstruga R, Dodds PN (2009) Terrific Protein Traffic: The Mystery of Effector Protein Delivery by Filamentous Plant Pathogens. Science 324: 748–750.

7. Stergiopoulos I, de Wit PJGM (2009) Fungal effector proteins. Annu Rev Phytopathol 47: 233–263.

8. Jonge (de) R, Bolton MD, Thomma BPHJ (2011) How filamentous pathogen co-opt plants: the ins and outs of fungal effectors. Curr Opin Plant Biol 14: 1–7.

9. Schmidt SM, Panstruga R (2011) Pathogenomics of fungal plant parasites: what have we learnt about pathogenesis? Curr Opin Plant Biol 14: 1–8.

10. Duplessis S, Joly DL, Doods PN (2012) Rust effectors. In: Martin F, Kamoun S, editors. Effectors in Plant-Microbe Interactions (1st edition). Chichester: John Wiley & Sons, Ltd. 155–193.

11. Hacquard S, Joly DL, Lin Y-C, Tisserant E, Feau N, et al.. (2012) A comprehensive analysis of genes encoding small secreted proteins identifies candidate effectors in Melampsora larici-populina (poplar leaf rust). Mol Plant Microbe Interact doi:10.1094/MPMI-09–11–0238.

12. Block A, Alfano JR (2011) Plant targets for Pseudomonas syringae type III effectors: virulence targets or guarded decoys? Curr Opin Microbiol 14: 39–46.

13. Koeck M, Hardham AR, Dodds PN (2011) The role of effectors in biotrophic and hemibiotrophic fungi in infection. Cell Microbiol 13: 1847–1859.

14. Plett JM, Kemppainen M, Kale SD, Kohler A, Legué V, et al. (2011) A Secreted Effector Protein of Laccaria bicolor Is Required for Symbiosis Development. Curr Biol 21: 1197–1203.

15. Caillaud M-C, Piquerez SJM, Fabro G, Steinbrenner J, Ishaque N, et al. (2012) Subcellular localization of the Hpa RxLR effector repertoire identifies a tonoplast-associated protein HaRxL17 that confers enhanced plant susceptibility. Plant J 69: 252–265.

16. Tuskan GA, Difazio S, Jansson S, Bohlmann J, Grigoriev I, et al. (2006) The genome of black cottonwood, Populus trichocarpa (Torr. & Gray). Science 313: 1596–1604.

17. Duplessis S, Cuomo CA, Lin Y-C, Aerts A, Tisserand E, et al.. (2011a) Obligate biotrophy features unraveled by the genomic analysis of rust fungi. Proc. Natl. Acad. Sci. U.S.A. doi/10.1073/pnas.1019315108.

18. Hacquard S, Petre B, Frey P, Hecker A, Rouhier N, et al.. (2011a) The Poplar-Poplar Rust Interaction: Insights from Genomics and Transcriptomics. Journal of Pathogens doi:10.4061/2011/716041.

19. Duplessis S, Major I, Martin F, Séguin A (2009) Poplar and Pathogen Interactions: Insights from Populus Genome-Wide Analyses of Resistance and Defense Gene Families and Gene Expression Profiling. Crit Rev Plant Sci 28: 309–334.

20. Hacquard S, Delaruelle C, Legué V, Tisserant E, Kohler A, et al. (2010) Laser Capture Microdissection of Uredinia Formed by Melampsora larici-populina Revealed a Transcriptional Switch Between Biotrophy and Sporulation. Mol Plant Microbe Interact 23: 1275–1286.

21. Duplessis S, Hacquard S, Delaruelle C, Tisserant E, Frey P, et al. (2011b) Melampsora larici-populina Transcripts Profiling During Germination and Time-course Infection of Poplar Leaves Reveals Dynamic Expression Patterns Associated with Virulence and Biotrophy. Mol Plant Microbe Interact 24: 808–818.

22. Rothberg JM, Leamon JH (2008) The development and impact of 454 sequencing. Nat Biotechnol 26: 1117–1124.

23. Simon SA, Zhai J, Nandety RS, McCormick KP, Zeng J, et al. (2009) Short-Read Sequencing Technologies for Transcriptional Analyses. Annu Rev Plant Biol 60: 305–333.

24. Metzker ML (2010) Sequencing technologies – the next generation. Nat Rev Genet 11: 31–46.

25. Kemen E, Gardiner A, Schultz-Larsen T, Kemen AC, Balmuth AL, et al. (2011) Gene Gain and Loss during Evolution of Obligate Parasitism in the White Rust Pathogen of Arabidopsis thaliana. Plos Biol 9(7): e1001094.

26. Tisserant E, Kohler A, Dozolme-Seddas P, Balestrini R, Benabdellah K, et al.. (2011) The transcriptome of the arbuscular mycorrhizal fungus Glomus intraradices (DAOM 197198) reveals functional tradeoffs in an obligat symbiont. New Phytol doi: 10.1111/j.1469–8137.2011.03948.x.

27. Fernandez D, Tisserand E, Talhinhas P, Azinheira H, Vieira A, et al. (2012) 454-pyrosequencing of Coffea arabica leaves infected by the rust fungus Hemeleia vastatrix reveals in planta-expressed pathogen-secreted proteins and plant functions in a late compatible plant-rust interaction. Mol Plant Pathol 5: 515–595.

28. Laurans F, Pilate G (1999) Histological Aspects of a Hypersensitive Response in Poplar to Melampsora larici-populina. Phytopathology 89: 223–238.

29. Rinaldi C, Kohler A, Frey P, Duchaussoy F, Ningre N, et al. (2007) Transcript Profiling of Poplar Leaves upon Infection with Compatible and Incompatible Strains of the Foliar Rust Melampsora larici-populina. Plant Physiol 144: 347–366.

30. Hahn M, Mendgen K (1997) Characterization of In Planta-Induced Rust Genes Isolated from a Hautorium-Specific cDNA Library. Mol Plant Microbe Interact 4: 427–437.

31. Kemen E, Kemen AC, Rafiqi M, Hempel U, Mendgen K, et al. (2005) Identification of a Protein from Rust Fungi Transferred from Haustoria into Infected Plant Cells. Mol Plant Microbe Interact 11: 1030–1039.

32. Catanzariti A-M, Dodds PN, Lawrence GJ, Ayliffe MA, Ellis JG (2006) Haustorially Expressed Secreted Proteins from Flax Rust Are Highly Enriched for Avirulence Elicitors. Plant Cell 18: 243–256.

33. Rafiqi M, Gan PHP, Ravensdale M, Laxrence GJ, Ellis JG, et al. (2010) Internalization of Flax Rust Avirulence Proteins into Flax and Tobacco Cell Can Occur in the Absence of the Pathogen. Plant Cell 22: 2017–2032.

34. Haas BJ, Delcher AL, Mount SM, Wortman JR, Smith RK, et al. (2003) Improving the Arabidopsis genome annotation using maximal transcript alignment assemblies. Nucl Acid Res 31: 5654–5666.

35. Hacquard S, Veneault-Fourrey C, Delaruelle C, Frey P, Martin F, et al. (2011b) 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 75: 106–112.

36. Weßling R, Schmidt SM, Micali CO, Knaust F, Reinhardt R, et al.. (2012) Transcriptome analysis of enriched Golovinomyces orontii haustoria by deep 454 pyrosequencing. Fung Genet Biol, doi:10.1016/j.fgb.2012.04.001.

37. Kale SD, Gu B, Capelluto DGS, Dou D, Feldman E, et al. (2010) External Lipid PI3P Mediates Entry of Eukaryotic Pathogen Effectors into Plant and Animal Host Cells. Cell 142: 284–295.

38. Kale SD, Tyler BM (2011) Entry of oomycete and fungal effectors into plant and animal host cells. Cell Microbiol 13: 1839–1848.

39. Yaeno T, Li H, Chaparro-Garcia A, Schornack S, Koshiba S, et al.. (2011) Phosphatidylinositol monophosphate-binding interface in the oomycete RXLR effector AVR3a is required for its stability in host cells to modulate plant immunity. Proc. Natl. Acad. Sci. U.S.A. doi/10.1073/pnas.1106002108.

40. Malone JH, Oliver B (2011) Microarrays, deep sequencing and the true measure of the transcriptome. BMC Biol 9: 34.

41. Graveley BR, Brooks AN, Carlson JW, Duff MO, Landolin JM, et al. (2011) The developmental transcriptome of Drosophila melanogaster. Nature 471: 473–479.

42. Kumar S, Hasan Asif M, Chakrabarty D, Deo Tripathi R, Kumar Trivedi P (2011) Differential expression and alternative splicing in rice sulphate transporter family members regulate sulphur status during plant growth, development and stress conditions. Funct Integr Genomic doi: 10.1007/s10142-010-0207-y.

43. Dürr J, Bücking H, Mult S, Wildhagen H, Palme K, et al. (2010) Seasonal and cell type specific expression of sulfate transporters in the phloem of Populus reveals tree specific characteristics for SO42− storage and mobilization. Plant Mol Biol 72: 499–517.

44. Wise RP, Moscou MJ, Bogdanove AJ, Whitham SA (2007) Transcript Profiling in Host-Pathogen Interactions. Annu Rev Phytopathol 45: 329–369.

45. Zipfel C (2009) Early molecular events in PAMP-triggered immunity. Curr Opin Plant Biol 12: 414–420.

46. Schneider KT, van de Mortel M, Bancroft TJ, Braun E, Nettleton D, et al. (2011) Biphasic gene expression changes elicited by Phakopsora pachyrhizi in soybean correlate with fungal penetration and haustoria formation. Plant Physiol 157: 355–371.

47. Petre B, Major I, Rouhier N, Duplessis S (2011) Genome-wide analysis of eukaryote thaumatin-like proteins (TLPs) with an emphasis on poplar. BMC Plant Biol 11: 33.

48. Kataoka T, Hayashi N, Yamaya T, Takahashi H (2004) Root-to-Shoot Transport of Sulfate in Arabidopsis. Evidence for the Role of SULTR3;5 as a Component of Low-Affinity Sulfate Transport System in the Root Vasculature. Plant Physiol 136: 4198–4204.

49. Wienkoop S, Saalbach G (2003) Proteome analysis: Novel proteins identified at the peribacteroid membrane from Lotus japonicus root nodules. Plant Physiol 131: 1080–1090.

50. Krusell L, Krause K, Ott T, Desbrosses G, Kramer U, et al. (2005) The Sulfate Transporter SST1 Is Crucial for Symbiotic Nitrogen Fixation in Lotus japonicus Root Nodules. Plant Cell 17: 1625–1636.

51. Felten J, Kohler A, Morin E, Bhalerao RP, Palme K, et al. (2009) The Ectomycorrhizal Fungus Laccaria bicolor Stimulates Lateral Roots Formation in Poplar and Arabidopsis through Auxin Transport and Signaling. Plant Physiol 151: 1991–2005.

52. Mukhtar MS, Carvunis A-R, Dreze M, Epple P, Steinbrenner J, et al. (2011) Independantly Evolved Virulence Effectors Converge onto Hubs in a Plant Immune System Network. Science 333: 596–601.

53. Baxter L, Tripathy S, Ishaque N, Boot N, Cabral A, et al. (2010) Signatures of Adaptation to Obligate Biotrophy in the Hyaloperonospora arabidopsidis Genome. Science 330: 1549–1551.

54. Spanu PD, Abbott JC, Amselem J, Burgis TA, Soanes DM, et al. (2010) Genome Expansion and Gene Loss in Powdery Mildew Fungi Reveal Tradeoffs in Extreme Parasitism. Science 330: 1543–1546.

55. Antony G, Zhou J, Huang S, Li T, Liu B, et al. (2010) Rice xa13 Recessive Resistance to Bacterial Blight Is Defeated by Induction of the Disease Susceptibility Gene Os-11N3. Plant Cell 22: 3864–3876.

56. Boch J, Bonas U (2010) Xanthomonas AvrBs3 Family-Type III Effectors: Discovery and Function. Annu Rev Phytopathol 48: 419–436.

57. Chen L-Q, Hou B-H, Lalonde S, Takanaga H, Hartung ML, et al. (2010) Sugar transporters for intercellular exchange and nutrition of pathogens. Nature 468: 527–532.

58. Talbot NJ (2010) Raiding the sweet shop. Nature 510: 468–469.

59. Schornack S, van Damme M, Bozkurt TO, Cano LM, Smoker M, et al. (2010) Ancient class of translocated oomycete effectors targets the host nucleus. Proc. Natl. Acad. Sci. U.S.A. 107: 17421–17426.

60. Dowkiw A, Husson C, Frey P, Pinon J, Bastien C (2003) Partial resistance to Melampsora larici-populina leaf rust in hybrid poplars: genetic variability in inoculated excised leaf disk bioassay and relationship with complete resistance. Phytopathology 93: 421–427.

61. Boyle B, Hamelin RC, Séguin A (2005) In vivo monitoring of obligate biotrophic pathogen growth by kinetic PCR. Appl Environ Microbiol 71: 1546–1552.

62. Boyle B, Levée V, Hamel L-P, Nicole M-C, Séguin A (2010) Molecular and histochemical characterisation of two distinct poplar Melampsora leaf rust pathosystems. Plant Biol 12: 364–376.

63. Pfaffl MW (2001) A new mathematical model for relative quantification in real-time RT-PCR. Nucl Acid Res 29: 2003–2007.

64. Chevreux B, Pfisterer T, Drescher B, Driesel AJ, Muller WEG, et al. (2004) Using the miraEST Assembler for Reliable and Automated mRNA Transcript Assembly and SNP Detection in Sequenced ESTs. Genome Res 14: 1147–1159.

65. Gremme G, Brendel V, Sparks ME, Kurtz S (2005) Engineering a software tool for gene structure prediction in higher organisms. Inform Software Technol 47: 965–978.

66. Milne I, Bayer M, Cardle L, Shaw P, Stephen G, et al. (2010) Tablet – next generation sequence assembly visualization. Bioinformatics 26: 401–402.

67. Ruepp A, Zollner A, Maier D, Albermann K, Hani J, et al. (2004) The FunCat, a functional annotation scheme for systematic classification of proteins from whole genomes. Nucl Acid Res 32: 5539–5545.

68. Marioni JC, Mason CE, Mane SM, Stephens M, Gilad Y (2008) RNA-seq: An assessment of technical reproductibility and comparison with gene expression arrays. Genome Res 18: 1509–1517.

69. Tsai C-J, Harding SA, Tschaplinski TJ, Lindroth RL, Yuan Y (2006) Genome-wide analysis of the structural genes regulating defense phenylpropanoid metabolism in Populus. New Phytol 172: 47–62.

70. Hamberger B, Ellis M, Friedmann M, de Azevedo Souza C, Barbazuk B, et al. (2007) Genome-wide analyses of phenylpropanoid-related genes in Populus trichocarpa, Arabidopsis thaliana, and Oryza sativa: the Populus lignin toolbox and conservation and diversification of angiosperm gene families. Can J Bot 85: 1182–1201.

71. Miranda M, Ralph SG, Mellway R, White R, Heath MC, et al. (2007) The Transcriptional Response of Hybrid Poplar (Populus trichocarpa x P. deltoides) to Infection by Melampsora medusae Leaf Rust Involves Induction of Flavonoid Pathway Genes Leading to the Accumulation of Proanthocyanidins. Mol Plant Microbe Interact 20: 816–831.

72. Azaiez A, Boyle B, Levée V, Séguin A (2009) Transcriptome profiling in hybrid poplar following interactions with Melampsora rust fungi. Mol Plant Microbe Interact 22: 190–200.