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

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.


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][8][9][10][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][14][15].
Genome sequences of Populus trichocarpa and Melampsora laricipopulina 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][23][24]. Recent works reported the use of RNA-Seq (Illumina or 454-pyrosequencing) to get insights into the transcriptomes of plant-microbe associations [25][26][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 wellannotated 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.

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. laricipopulina has been established and tissues were harvested at 18, 21 and 24 hours post-inoculation (hpi) for the incompatible interaction (I 18 -I 21 -I 24 ) and 18, 24 and 48 hpi for the compatible interaction (C 18 -C 24 -C 48 ) (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][31][32][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 earlyregulated poplar genes not previously reported. A later time-point in the compatible interaction (C 48 ) 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 'contigblast 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. laricipopulina transcripts.

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 48 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 18 [37][38][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.
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 5265 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 flyhead 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   (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 nonnegligeable technical biase between these approaches discussed in Text S1 (Fig. S6).

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.
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 48 (515 reads) with an almost 20-fold up-regulation compared with the other five conditions that remain stable (26615 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 consid- ering the specificity and amplitude of the induction of PtSultr3;5, we hypothesized that this specific up-regulation could result from a pathogen manipulation.
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 hosttranscriptome 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][56][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. presymbiotic 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 plantmicrobe interactions.

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 mg of starting total RNA and 19 cycles of PCR amplification in 25 mL. 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 mg of total RNA and following manufacturer instructions.

454-pyrosequencing and Oligoarrays
Sequences were generated from 5 mg 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 18 , I 21 , I 24 conditions mixed 1/3-1/3-1/3 together with specific tags; and another one for the C 18 , C 24 , C 48 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) 2Ct ]/ [(reference gene primer efficiency) 2Ct ], 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 reannotation 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/). 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 doublehead arrow marks additional predicted amino acids in the P. trichocarpa 'Nisqually-1' genome discussed in Test S3. (TIF) Text S1 Uncorrelated oligoarray/RNA-Seq transcript expression.