The Rust Fungus Melampsora larici-populina Expresses a Conserved Genetic Program and Distinct Sets of Secreted Protein Genes During Infection of Its Two Host Plants, Larch and Poplar.

Mechanisms required for broad-spectrum or specific host colonization of plant parasites are poorly understood. As a perfect illustration, heteroecious rust fungi require two alternate host plants to complete their life cycles. Melampsora larici-populina infects two taxonomically unrelated plants, larch, on which sexual reproduction is achieved, and poplar, on which clonal multiplication occurs, leading to severe epidemics in plantations. We applied deep RNA sequencing to three key developmental stages of M. larici-populina infection on larch: basidia, pycnia, and aecia, and we performed comparative transcriptomics of infection on poplar and larch hosts, using available expression data. Secreted protein was the only significantly overrepresented category among differentially expressed M. larici-populina genes between the basidial, the pycnial, and the aecial stages, highlighting their probable involvement in the infection process. Comparison of fungal transcriptomes in larch and poplar revealed a majority of rust genes were commonly expressed on the two hosts and a fraction exhibited host-specific expression. More particularly, gene families encoding small secreted proteins presented striking expression profiles that highlight probable candidate effectors specialized on each host. Our results bring valuable new information about the biological cycle of rust fungi and identify genes that may contribute to host specificity.

Patterns of genomic variation in the poplar rust fungus Melampsora larici-populina identify pathogenesis-related factors.

Melampsora larici-populina is a fungal pathogen responsible for foliar rust disease on poplar trees, which causes damage to forest plantations worldwide, particularly in Northern Europe. The reference genome of the isolate 98AG31 was previously sequenced using a whole genome shotgun strategy, revealing a large genome of 101 megabases containing 16,399 predicted genes, which included secreted protein genes representing poplar rust candidate effectors. In the present study, the genomes of 15 isolates collected over the past 20 years throughout the French territory, representing distinct virulence profiles, were characterized by massively parallel sequencing to assess genetic variation in the poplar rust fungus. Comparison to the reference genome revealed striking structural variations. Analysis of coverage and sequencing depth identified large missing regions between isolates related to the mating type loci. More than 611,824 single-nucleotide polymorphism (SNP) positions were uncovered overall, indicating a remarkable level of polymorphism. Based on the accumulation of non-synonymous substitutions in coding sequences and the relative frequencies of synonymous and non-synonymous polymorphisms (i.e., PN/PS ), we identify candidate genes that may be involved in fungal pathogenesis. Correlation between non-synonymous SNPs in genes encoding secreted proteins (SPs) and pathotypes of the studied isolates revealed candidate genes potentially related to virulences 1, 6, and 8 of the poplar rust fungus.

Transcriptome analysis of poplar rust telia reveals overwintering adaptation and tightly coordinated karyogamy and meiosis processes.

Most rust fungi have a complex life cycle involving up to five different spore-producing stages. The telial stage that produces melanized overwintering teliospores is one of these and plays a fundamental role for generating genetic diversity as karyogamy and meiosis occur at that stage. Despite the importance of telia for the rust life cycle, almost nothing is known about the fungal genetic programs that are activated in this overwintering structure. In the present study, the transcriptome of telia produced by the poplar rust fungus Melampsora larici-populina has been investigated using whole genome exon oligoarrays and RT-qPCR. Comparative expression profiling at the telial and uredinial stages identifies genes specifically expressed or up-regulated in telia including osmotins/thaumatin-like proteins (TLPs) and aquaporins that may reflect specific adaptation to overwintering as well numerous lytic enzymes acting on plant cell wall, reflecting extensive cell wall remodeling at that stage. The temporal dynamics of karyogamy was followed using combined RT-qPCR and DAPI-staining approaches. This reveals that fusion of nuclei and induction of karyogamy-related genes occur simultaneously between the 25 and 39 days post inoculation time frame. Transcript profiling of conserved meiosis genes indicates a preferential induction right after karyogamy and corroborates that meiosis begins prior to overwintering and is interrupted in Meiosis I (prophase I, diplonema stage) until teliospore germination in early spring.

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.

A comprehensive analysis of genes encoding small secreted proteins identifies candidate effectors in Melampsora larici-populina (poplar leaf rust).

The obligate biotrophic rust fungus Melampsora larici-populina is the most devastating and widespread pathogen of poplars. Studies over recent years have identified various small secreted proteins (SSP) from plant biotrophic filamentous pathogens and have highlighted their role as effectors in host-pathogen interactions. The recent analysis of the M. larici-populina genome sequence has revealed the presence of 1,184 SSP-encoding genes in this rust fungus. In the present study, the expression and evolutionary dynamics of these SSP were investigated to pinpoint the arsenal of putative effectors that could be involved in the interaction between the rust fungus and poplar. Similarity with effectors previously described in Melampsora spp., richness in cysteines, and organization in large families were extensively detailed and discussed. Positive selection analyses conducted over clusters of paralogous genes revealed fast-evolving candidate effectors. Transcript profiling of selected M. laricipopulina SSP showed a timely coordinated expression during leaf infection, and the accumulation of four candidate effectors in distinct rust infection structures was demonstrated by immunolocalization. This integrated and multifaceted approach helps to prioritize candidate effector genes for functional studies.

Melampsora larici-populina transcript profiling during germination and timecourse infection of poplar leaves reveals dynamic expression patterns associated with virulence and biotrophy.

Melampsora larici-populina is responsible for poplar leaf rust disease and causes severe epidemics in poplar plantations in Europe. The poplar rust genome has been recently sequenced and, in order to find the genetic determinants associated with its biotrophic lifestyle, we generated a whole-genome custom oligoarray and analyzed transcript profiles of M. larici-populina during the infection timecourse in poplar leaves. Different stages were investigated during the asexual development of the rust fungus, including resting and germinating urediniospores and seven in planta stages in the telial host. In total, 76% of the transcripts were detected during leaf infection as well as in urediniospores, whereas 20% were only detected in planta, including several transporters and many small secreted proteins (SSP). We focused our analysis on gene categories known to be related to plant colonization and biotrophic growth in rust pathogens, such as SSP, carbohydrate active enzymes (CAZymes), transporters, lipases, and proteases. Distinct sets of SSP transcripts were expressed all along the infection process, suggesting highly dynamic expression of candidate rust effectors. In contrast, transcripts encoding transporters and proteases were mostly expressed after 48 h postinoculation, when numerous haustoria are already formed in the leaf mesophyll until uredinia formation, supporting their role in nutrient acquisition during biotrophic growth. Finally, CAZymes and lipase transcripts were predominantly expressed at late stages of infection, highlighting their importance during sporulation.

Laser capture microdissection of uredinia formed by Melampsora larici-populina revealed a transcriptional switch between biotrophy and sporulation.

The foliar rust caused by the basidiomycete Melampsora larici-populina is the main disease affecting poplar plantations in Europe. The biotrophic status of rust fungi is a major limitation to study gene expression of cell or tissue types during host infection. At the uredinial stage, infected poplar leaves contain distinct rust tissues such as haustoria, infection hyphae, and uredinia with sporogenous hyphae and newly formed asexual urediniospores. Laser capture microdissection (LCM) was used to isolate three areas corresponding to uredinia and subjacent zones in the host mesophyll for expression analysis with M. larici-populina whole-genome exon oligoarrays. Optimization of tissue preparation prior to LCM allowed isolation of RNA of good integrity for genome-wide expression profiling. Our results indicate that the poplar rust uredinial stage is marked by distinct genetic programs related to biotrophy in the host palisade mesophyll and to sporulation in the uredinium. A strong induction of transcripts encoding small secreted proteins, likely containing rust effectors, is observed in the mesophyll, suggesting a late maintenance of suppression of host defense in the tissue containing haustoria and infection hyphae. On the other hand, cell cycle and cell defense rescue transcripts are strongly accumulated in the sporulation area. This combined LCM-transcriptomic approach brings new insights on the molecular mechanisms underlying urediniospore formation in rust fungi.

Morphological and molecular analyses in Scleroderma species associated with some Caesalpinioid legumes, Dipterocarpaceae and Phyllanthaceae trees in southern Burkina Faso.

A combination of morphotypes, polymerase chain reaction/restriction fragment length polymorphism (PCR/RFLP) analyses and internal transcribed spacer (ITS) sequencing was used to investigate Scleroderma species that were collected from woodlands in Burkina Faso. We harvested 52 specimens from 20 sites during rainy seasons between 1997 and 2000. According to their morphological features, these specimens were initially characterised, and we then identified six species of Scleroderma. Two of the species were clearly determined as Scleroderma dictyosporum Pat. and S. verrucosum Pers. The four remaining species were characteristically described as (1) displaying big spores with spines up to 2 microm (Scleroderma sp1), (2) producing spores without ornamentation (Scleroderma sp2), (3) spores with very small spines (Scleroderma sp3) and (4) with yellow sporocarps and sub-spherical spores (Scleroderma sp4). The specimens were then analysed using PCR/RFLP of the intergenic regions of rDNA, ITS and IGS1 and ITS sequencing. The restriction fragments obtained with two endonucleases, HinfI and MboI on ITS and IGS1 regions, showed that some isolates of S. dictyosporum had the same patterns as isolates and basidiocarps of Scleroderma sp4 (IR265, IR408, SP4-2903). Isolates of Scleroderma sp3 (IR252) had common restriction fragments as isolates of S. verrucosum (IR500, IR600). Intraspecific differences were observed in the two previously determined species, as well as in Scleroderma sp2. The ITS sequencing and phylogenetic analyses showed that the ribotypes identified by PCR/RFLP within these species might be phylogenetic species. Combining these molecular results allowed regrouping the six morphological species in three sets of cryptic species: a first set with two species including S. dictyosporum Pat., a second set with four species, including both S. verrucosum Pers. and Scleroderma sp1 and a third set with two species, including Scleroderma sp2. These investigations and the combined morphological and molecular analyses used to sort out species paved the way for identifying larger populations of Scleroderma species in Burkina Faso and other tropical zones.

Polymorphism at the ribosomal DNA ITS and its relation to postglacial re-colonization routes of the Perigord truffle Tuber melanosporum.

• Glaciations and postglacial migrations are major factors responsible for the present patterns of genetic variation we see in natural populations in Europe. For ectomycorrhizal fungi, escape from refugia can only follow range expansion by their specific hosts. • To infer phylogeographic relationships within Tuber melanosporum, sequences of internal transcribed spacers (ITS) and the 5.8S coding region of the ribosomal DNA repeat were obtained for 188 individuals sampled over the entire distribution of this species in France, and in north-western Italy and north-eastern Spain. • Ten distinct ITS haplotypes were distinguished, mapped and treated using F- and NST -statistics and nested clade (NCA) analyses. They showed a significant genetic differentiation between regional populations. NCA revealed a geographical association of ITS haplotypes, an old fragmentation into two major groups of populations, which likely colonized regions on different sides of the French Central Massif. • This re-colonization pattern is reminiscent of the one observed for host trees of the Perigord truffle, such as oaks and hazelnut trees. This suggests that host postglacial expansion was one of the major factors that shaped the mycobiont population structure.

The poplar root transcriptome: analysis of 7000 expressed sequence tags.

To date, most poplar expressed sequence tags (ESTs) are from above-ground tissues such as wood, leaf and buds. Here, we present a large-scale production of ESTs from roots of the hybrid cottonwood, Populus trichocarpaxdeltoides. cDNA libraries were generated from the root system of 2-month-old rooted cuttings, and roots of 2.5-month-old cuttings water-stressed for 19 days. Partial sequences obtained from 7013 clones were assembled into 1347 clusters and 3527 singletons. This set of ESTs represents 4874 unique transcripts expressed in roots. Putative functions could be assigned to 3021 (62%) of the transcripts. A significant portion of the ESTs encode proteins of common metabolic pathways; energy and metabolism represented 5% and 8% of total transcripts, respectively. Of specific interest to root functions are the 6% of ESTs involved in signalling pathways and hormone metabolism, and 4% encoding transporters and channels. The current poplar root ESTs and the aspen root ESTs present in public databases represent 6700 unique transcripts. The Unigene set was selected from the ESTs and used to generate nylon microarrays. Changes in aquaporins and transporter transcripts were then studied during adventitious root development.

Analysis of expressed sequence tags from the ectomycorrhizal basidiomycetes Laccaria bicolor and Pisolithus microcarpus.

• In an effort to discover genes that are expressed in the ectomycorrhizal basidiomycetes Laccaria bicolor and Pisolithus microcarpus, and in P. microcarpus/Eucalyptus globulus ectomycorrhizas, we have sequenced 1519 and 1681 expressed sequence tags (ESTs) from L. bicolor and P. microcarpus cDNA libraries. • Contig analysis resulted in 905 and 806 tentative consensus sequences (unique transcripts) in L. bicolor and P. microcarpus, respectively. For 36% of the ESTs, significant similarities to sequences in databases were detected. The most abundant transcripts showed no similarity to previously identified genes. Sequence redundancy analysis between different developmental stages indicated that several genes were differentially expressed in free-living mycelium and symbiotic tissues of P. microcarpus. • Based on sequence similarity, 11% of L. bicolor unique transcripts were also detected in P. microcarpus. Similarly, L. bicolor and P. microcarpus shared only a low proportion of common transcripts with other basidiomycetous fungi, such as Pleurotus ostreatus and Agaricus bisporus. Such a low proportion of shared transcripts between basidiomycetes suggests, on the one hand, that the variability of expressed transcripts in different fungi and fungal tissues is considerably high. On the other hand, it might reflect the low number of GenBank entries of basidiomycetous origin and stresses the necessity of an additional sequencing effort. • The present ESTs provide a valuable resource for future research on the development and functioning of ectomycorrhizas.