Comparison of Genetic and Virulence Diversity of Melampsora larici-populina Populations on Wild and Cultivated Poplar and Influence of the Alternate Host

Approximate Bayesian Computation applied to time series of population genetic data disentangles rapid genetic changes and demographic variations in a pathogen population

Adaptation can occur at remarkably short timescales in natural populations, leading to drastic changes in phenotypes and genotype frequencies over a few generations only. The inference of demographic parameters can allow understanding how evolutionary forces interact and shape the genetic trajectories of populations during rapid adaptation. Here we propose a new Approximate Bayesian Computation (ABC) framework that couples a forward and individual-based model with temporal genetic data to disentangle genetic changes and demographic variations in a case of rapid adaptation. We test the accuracy of our inferential framework and evaluate the benefit of considering a dense versus sparse sampling. Theoretical investigations demonstrate high accuracy in both model and parameter estimations, even if a strong thinning is applied to time series data. Then, we apply our ABC inferential framework to empirical data describing the population genetic changes of the poplar rust pathogen following a major event of resistance overcoming. We successfully estimate key demographic and genetic parameters, including the proportion of resistant hosts deployed in the landscape and the level of standing genetic variation from which selection occurred. Inferred values are in accordance with our empirical knowledge of this biological system. This new inferential framework, which contrasts with coalescent-based ABC analyses, is promising for a better understanding of evolutionary trajectories of populations subjected to rapid adaptation.

The Melampsora americana Population on Salix purpurea in the Great Lakes Region Is Highly Diverse with a Contributory Influence of Clonality

Shrub willows (Salix spp.) are emerging as a viable lignocellulosic, second-generation bioenergy crop with many growth characteristics favorable for marginal lands in New York State and surrounding areas. Willow rust, caused by members of the genus Melampsora, is the most limiting disease of shrub willow in this region and remains extremely understudied. In this study, genetic diversity, genetic structure, and pathogen clonality were examined in Melampsora americana over two growing seasons via genotyping-by-sequencing to identify single-nucleotide polymorphism markers. In conjunction with this project, a reference genome of rust isolate R15-033-03 was generated to aid in variant discovery. Sampling between years allowed regional and site-specific investigation into population dynamics, in the context of both wild and cultivated hosts within high-density plantings. This work revealed that this pathogen is largely panmictic over the sampled areas, with few sites showing moderate genetic differentiation. These data support the hypothesis of sexual recombination between growing seasons because no genotype persisted across the two years of sampling. Additionally, clonality was determined as a driver of pathogen populations within cultivated fields and single shrubs; however, there is also evidence of high genetic diversity of rust isolates in all settings. This work provides a framework for M. americana population structure in the Great Lakes region, providing crucial information that can aid in future resistance breeding efforts.

Genomic signatures of a major adaptive event in the pathogenic fungus Melampsora larici-populina

The recent availability of genome-wide sequencing techniques has allowed systematic screening for molecular signatures of adaptation, including in nonmodel organisms. Host–pathogen interactions constitute good models due to the strong selective pressures that they entail. We focused on an adaptive event which affected the poplar rust fungus Melampsora larici-populina when it overcame a resistance gene borne by its host, cultivated poplar. Based on 76 virulent and avirulent isolates framing narrowly the estimated date of the adaptive event, we examined the molecular signatures of selection. Using an array of genome scan methods based on different features of nucleotide diversity, we detected a single locus exhibiting a consistent pattern suggestive of a selective sweep in virulent individuals (excess of differentiation between virulent and avirulent samples, linkage disequilibrium, genotype–phenotype statistical association, and long-range haplotypes). Our study pinpoints a single gene and further a single amino acid replacement which may have allowed the adaptive event. Although our samples are nearly contemporary to the selective sweep, it does not seem to have affected genome diversity further than the immediate vicinity of the causal locus, which can be explained by a soft selective sweep (where selection acts on standing variation) and by the impact of recombination in mitigating the impact of selection. Therefore, it seems that properties of the life cycle of M. larici-populina, which entails both high genetic diversity and outbreeding, has facilitated its adaptation.

A point mutation and large deletion at the candidate avirulence locus AvrMlp7 in the poplar rust fungus correlate with poplar RMlp7 resistance breakdown

The deployment of plant varieties carrying resistance genes (R) exerts strong selection pressure on pathogen populations. Rapidly evolving avirulence genes (Avr) allow pathogens to escape R-mediated plant immunity through a variety of mechanisms, leading to virulence. The poplar rust fungus Melampsora larici-populina is a damaging pathogen of poplars in Europe. It underwent a major adaptive event in 1994, with the breakdown of the poplar RMlp7 resistance gene. Population genomics studies identified a locus in the genome of M. larici-populina that probably corresponds to the candidate avirulence gene AvrMlp7. Here, to further characterize this effector, we used a population genetics approach on a comprehensive set of 281 individuals recovered throughout a 28-year period encompassing the resistance breakdown event. Using two dedicated molecular tools, genotyping at the candidate locus highlighted two different alterations of a predominant allele found mainly before the resistance breakdown: a nonsynonymous mutation and a complete deletion of this locus. This results in six diploid genotypes: three genotypes related to the avirulent phenotype and three related to the virulent phenotype. The temporal survey of the candidate locus revealed that both alterations were found in association during the resistance breakdown event. They pre-existed before the breakdown in a heterozygous state with the predominant allele cited above. Altogether, these results suggest that the association of both alterations at the candidate locus AvrMlp7 drove the poplar rust adaptation to RMlp7-mediated immunity. This study demonstrates for the first time a case of adaptation from standing genetic variation in rust fungi during a qualitative resistance breakdown.

Genetic analysis reveals long-standing population differentiation and high diversity in the rust pathogen Melampsora lini

A priority for research on infectious disease is to understand how epidemiological and evolutionary processes interact to influence pathogen population dynamics and disease outcomes. However, little is understood about how population adaptation changes across time, how sexual vs. asexual reproduction contribute to the spread of pathogens in wild populations and how diversity measured with neutral and selectively important markers correlates across years. Here, we report results from a long-term study of epidemiological and genetic dynamics within several natural populations of the Linum marginale-Melampsora lini plant-pathogen interaction. Using pathogen isolates collected from three populations of wild flax (L. marginale) spanning 16 annual epidemics, we probe links between pathogen population dynamics, phenotypic variation for infectivity and genomic polymorphism. Pathogen genotyping was performed using 1567 genome-wide SNP loci and sequence data from two infectivity loci (AvrP123, AvrP4). Pathogen isolates were phenotyped for infectivity using a differential set. Patterns of epidemic development were assessed by conducting surveys of infection prevalence in one population (Kiandra) annually. Bayesian clustering analyses revealed host population and ecotype as key predictors of pathogen genetic structure. Despite strong fluctuations in pathogen population size and severe annual bottlenecks, analysis of molecular variance revealed that pathogen population differentiation was relatively stable over time. Annually, varying levels of clonal spread (0–44.8%) contributed to epidemics. However, within populations, temporal genetic composition was dynamic with rapid turnover of pathogen genotypes, despite the dominance of only four infectivity phenotypes across the entire study period. Furthermore, in the presence of strong fluctuations in population size and migration, spatial selection may maintain pathogen populations that, despite being phenotypically stable, are genetically highly dynamic.

Melampsora lini is a rust fungus that infects native flax, Linum marginale in south-eastern Australia where its epidemiology and evolution have been intensively studied since 1987. Over that time, substantial diversity in the pathotypic structure of M. lini has been demonstrated but an understanding of how genetic diversity in pathogen populations is maintained through space and time is lacking. Here we integrated phenotypic, genotypic and epidemiological datasets spanning 16 annual epidemics across three host populations to examine long-term pathogen genetic dynamics. The results show that host ecotype is the dominant selective force in the face of strong bottlenecks and annual patterns of genetic turnover. Results from previous studies indicate that in this geographic region, M. lini lacks the capacity to reproduce sexually–we thus expected to find limited genetic diversity and evidence for strong clonality influencing genetic dynamics within growing seasons. However, the breadth of genomic coverage provided by the SNP markers revealed high levels of genotypic variation within M. lini populations. This discovery contrasts with observed phenotypic dynamics as the epidemics of this pathogen were largely dominated by four pathotypes across the study period. Based on a detailed assessment and comparison of pathotypic and genotypic patterns, our study increases the understanding of how genetic diversity is generated and maintained through space and time within wild pathogen populations. The implications for the management of resistance to pathogens in agricultural or conservation contexts are significant: the appearance of clonality may be hiding high levels of pathogen diversity and recombination. Understanding how this diversity is generated could provide new and unique ways to mitigate or suppress the emergence of infectious strains, allowing to efficiently combat harmful diseases.

Identification of Virulences of the Rust Fungus Melampsora larici-populina Occurring in Chile

Melampsora larici-populina Kleb. is a Eurasian species that causes leaf rust of many species of Populus from sections Tacamahaca and Aigeiros. This rust has been inadvertently introduced to many parts around the world affecting poplar species that grow naturally, which may be susceptible to this fungus. In Chile, early European settlers introduced poplars; rust, attributed to M. larici-populina, has been reported in these trees since 1918. However, a modern confirmation has been lacking, and pathogenic variation of the Chilean population of M. larici-populina has not been investigated. Using a morphological analysis of urediniospores and sequencing of internal transcribed spacer (ITS) regions, we confirmed the presence of M. larici-populina in Chile. ITS regions exhibited 100% homology with M. larici-populina. Scanning electron microscopy showed that spores were 26 to 47 µm in length, 13 to 16 µm in width, and echinulate except for apices, which are smooth, characteristics described for this species of rust. The variability of M. larici-populina is characterized by the presence of pathotypes, which allows the fungus to infect despite the resistance of certain poplar hybrids. We concluded that the identified spores belong to M. larici-populina, with virulences 1, 2, 3, 4, 5, and 6. These results describe variation in virulence of M. larici-populina, which suggests the presence of a sexual stage in Chile.

Inferring pathogen dynamics from temporal count data: the emergence of Xylella fastidiosa in France is probably not recent

Unravelling the ecological structure of emerging plant pathogens persisting in multi-host systems is challenging. In such systems, observations are often heterogeneous with respect to time, space and host species, and may lead to biases of perception. The biased perception of pathogen ecology may be exacerbated by hidden fractions of the whole host population, which may act as infection reservoirs. We designed a mechanistic-statistical approach to help understand the ecology of emerging pathogens by filtering out some biases of perception. This approach, based on SIR (Susceptible-Infected-Removed) models and a Bayesian framework, disentangles epidemiological and observational processes underlying temporal counting data. We applied our approach to French surveillance data on Xylella fastidiosa, a multi-host pathogenic bacterium recently discovered in Corsica, France. A model selection led to two diverging scenarios: one scenario without a hidden compartment and an introduction around 2001, and the other with a hidden compartment and an introduction around 1985. Thus, Xylella fastidiosa was probably introduced into Corsica much earlier than its discovery, and its control could be arduous under the hidden compartment scenario. From a methodological perspective, our approach provides insights into the dynamics of emerging plant pathogens and, in particular, the potential existence of infection reservoirs.

Rust Infection of Black Poplar Trees Reduces Photosynthesis but Does Not Affect Isoprene Biosynthesis or Emission

Poplar (Populus spp.) trees are widely distributed and play an important role in ecological communities and in forestry. Moreover, by releasing high amounts of isoprene, these trees impact global atmospheric chemistry. One of the most devastating diseases for poplar is leaf rust, caused by fungi of the genus Melampsora. Despite the wide distribution of these biotrophic pathogens, very little is known about their effects on isoprene biosynthesis and emission. We therefore infected black poplar (P. nigra) trees with the rust fungus M. larici-populina and monitored isoprene emission and other physiological parameters over the course of infection to determine the underlying mechanisms. We found an immediate and persistent decrease in photosynthesis during infection, presumably caused by decreased stomatal conductance mediated by increased ABA levels. At the same time, isoprene emission remained stable during the time course of infection, consistent with the stability of its biosynthesis. There was no detectable change in the levels of intermediates or gene transcripts of the methylerythritol 4-phosphate (MEP) pathway in infected compared to control leaves. Rust infection thus does not affect isoprene emission, but may still influence the atmosphere via decreased fixation of CO₂.

Testing Differences Between Pathogen Compositions with Small Samples and Sparse Data

The structure of pathogen populations is an important driver of epidemics affecting crops and natural plant communities. Comparing the composition of two pathogen populations consisting of assemblages of genotypes or phenotypes is a crucial, recurrent question encountered in many studies in plant disease epidemiology. Determining whether there is a significant difference between two sets of proportions is also a generic question for numerous biological fields. When samples are small and data are sparse, it is not straightforward to provide an accurate answer to this simple question because routine statistical tests may not be exactly calibrated. To tackle this issue, we built a computationally intensive testing procedure, the generalized Monte Carlo plug-in test with calibration test, which is implemented in an R package available at https://doi.org/10.5281/zenodo.635791 . A simulation study was carried out to assess the performance of the proposed methodology and to make a comparison with standard statistical tests. This study allows us to give advice on how to apply the proposed method, depending on the sample sizes. The proposed methodology was then applied to real datasets and the results of the analyses were discussed from an epidemiological perspective. The applications to real data sets deal with three topics in plant pathology: the reproduction of Magnaporthe oryzae, the spatial structure of Pseudomonas syringae, and the temporal recurrence of Puccinia triticina.

The escalatory Red Queen: Population extinction and replacement following arms race dynamics in poplar rust

Host-parasite systems provide convincing examples of Red Queen co-evolutionary dynamics. Yet, a key process underscored in Van Valen's theory - that arms race dynamics can result in extinction - has never been documented. One reason for this may be that most sampling designs lack the breadth needed to illuminate the rapid pace of adaptation by pathogen populations. In this study, we used a 25-year temporal sampling to decipher the demographic history of a plant pathogen: the poplar rust fungus, Melampsora larici-populina. A major adaptive event occurred in 1994 with the breakdown of R7 resistance carried by several poplar cultivars widely planted in Western Europe since 1982. The corresponding virulence rapidly spread in M. larici-populina populations and nearly reached fixation in northern France, even on susceptible hosts. Using both temporal records of virulence profiles and temporal population genetic data, our analyses revealed that (i) R7 resistance breakdown resulted in the emergence of a unique and homogeneous genetic group, the so-called cultivated population, which predominated in northern France for about 20 years, (ii) selection for Vir7 individuals brought with it multiple other virulence types via hitchhiking, resulting in an overall increase in the population-wide number of virulence types and (iii) - above all - the emergence of the cultivated population superseded the initial population which predominated at the same place before R7 resistance breakdown. Our temporal analysis illustrates how antagonistic co-evolution can lead to population extinction and replacement, hence providing direct evidence for the escalation process which is at the core of Red Queen dynamics.

Weeds, as ancillary hosts, pose disproportionate risk for virulent pathogen transfer to crops

Background

The outcome of the arms race between hosts and pathogens depends heavily on the interactions between their genetic diversity, population size and transmission ability. Theory predicts that genetically diverse hosts will select for higher virulence and more diverse pathogens than hosts with low genetic diversity. Cultivated hosts typically have lower genetic diversity and thus small effective population sizes, but can potentially harbour large pathogen population sizes. On the other hand, hosts, such as weeds, which are genetically more diverse and thus have larger effective population sizes, usually harbour smaller pathogen population sizes. Large pathogen population sizes may lead to more opportunities for mutation and hence more diverse pathogens. Here we test the predictions that pathogen neutral genetic diversity will increase with large pathogen population sizes and host diversity, whereas diversity under selection will increase with host diversity. We assessed and compared the diversity of a fungal pathogen, Rhynchosporium commune, on weedy barley grass (which have a large effective population size) and cultivated barley (low genetic diversity) using microsatellites, effector locus nip1 diversity and pathogen aggressiveness in order to assess the importance of weeds in the evolution of the neutral and selected diversity of pathogens.

Results

The findings indicated that the large barley acreage and low host diversity maintains higher pathogen neutral genetic diversity and lower linkage disequilibrium, while the weed maintains more pathotypes and higher virulence diversity at nip1. Strong evidence for more pathogen migration from barley grass to barley suggests transmission of virulence from barley grass to barley is common.

Conclusions

Pathogen census population size is a better predictor for neutral genetic diversity than host diversity. Despite maintaining a smaller pathogen census population size, barley grass acts as an important ancillary host to R. commune, harbouring highly virulent pathogen types capable of transmission to barley. Management of disease on crops must therefore include management of weedy ancillary hosts, which may harbour disproportionate supplies of virulent pathogen strains.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0680-6) contains supplementary material, which is available to authorized users.

Detection and phylogenetic relationships of Puccinia emaculata and Uromyces graminicola (Pucciniales) on switchgrass in New York State using rDNA sequence information

The species of rust fungi (Pucciniales) inciting disease on switchgrass (Panicum virgatum) grown in bioenergy feedstock systems across the north-central and eastern United States remain unclear. In the present study, the species number and phylogenetic relationships of rust species affecting switchgrass were examined in 2011-2013 at two sites in New York State as well as selected sites in Alabama, Iowa, Nebraska, Pennsylvania, South Dakota, and West Virginia using ribosomal RNA gene data (partial internal transcribed spacer [ITS] 1, complete 5.8 subunit [S] and ITS2, and partial 28S). Uredinial group and teliospore morphology were also utilized to delimit taxa in collection years 2012 and 2013. Maximum likelihood, maximum parsimony, and Bayesian analyses demonstrated two monophyletic clades. Clade I consisted of Puccinia emaculata and included the majority of single-sorus samples across sites, whereas, Clade II included multiple samples from Iowa, Nebraska, and South Dakota. Single-telial samples for Clade I possessed only two-celled teliospores while Clade II samples possessed only one-celled teliospores, and hence, were readily diagnosed morphologically to P. emaculata and Uromyces graminicola, respectively. No U. graminicola sequences exist in GenBank to compare with our Clade II samples; however, based on teliospore morphology, the identity of Clade II taxa is U. graminicola.

Genome-wide patterns of segregation and linkage disequilibrium: the construction of a linkage genetic map of the poplar rust fungus Melampsora larici-populina

The poplar rust fungus Melampsora larici-populina causes significant yield reduction and severe economic losses in commercial poplar plantations. After several decades of breeding for qualitative resistance and subsequent breakdown of the released resistance genes, breeders now focus on quantitative resistance, perceived to be more durable. But quantitative resistance also can be challenged by an increase of aggressiveness in the pathogen. Thus, it is of primary importance to better understand the genetic architecture of aggressiveness traits. To this aim, our goal is to build a genetic linkage map for M. larici-populina in order to map quantitative trait loci related to aggressiveness. First, a large progeny of M. larici-populina was generated through selfing of the reference strain 98AG31 (which genome sequence is available) on larch plants, the alternate host of the poplar rust fungus. The progeny's meiotic origin was validated through a segregation analysis of 115 offspring with 14 polymorphic microsatellite markers, of which 12 segregated in the expected 1:2:1 Mendelian ratio. A microsatellite-based linkage disequilibrium analysis allowed us to identify one potential linkage group comprising two scaffolds. The whole genome of a subset of 47 offspring was resequenced using the Illumina HiSeq 2000 technology at a mean sequencing depth of 6X. The reads were mapped onto the reference genome of the parental strain and 144,566 SNPs were identified across the genome. Analysis of distribution and polymorphism of the SNPs along the genome led to the identification of 2580 recombination blocks. A second linkage disequilibrium analysis, using the recombination blocks as markers, allowed us to group 81 scaffolds into 23 potential linkage groups. These preliminary results showed that a high-density linkage map could be constructed by using high-quality SNPs based on low-coverage resequencing of a larger number of M. larici-populina offspring.

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.

Association analysis identifies Melampsora ×columbiana poplar leaf rust resistance SNPs

Populus species are currently being domesticated through intensive time- and resource-dependent programs for utilization in phytoremediation, wood and paper products, and conversion to biofuels. Poplar leaf rust disease can greatly reduce wood volume. Genetic resistance is effective in reducing economic losses but major resistance loci have been race-specific and can be readily defeated by the pathogen. Developing durable disease resistance requires the identification of non-race-specific loci. In the presented study, area under the disease progress curve was calculated from natural infection of Melampsora ×columbiana in three consecutive years. Association analysis was performed using 412 P. trichocarpa clones genotyped with 29,355 SNPs covering 3,543 genes. We found 40 SNPs within 26 unique genes significantly associated (permutated P<0.05) with poplar rust severity. Moreover, two SNPs were repeated in all three years suggesting non-race-specificity and three additional SNPs were differentially expressed in other poplar rust interactions. These five SNPs were found in genes that have orthologs in Arabidopsis with functionality in pathogen induced transcriptome reprogramming, Ca²⁺/calmodulin and salicylic acid signaling, and tolerance to reactive oxygen species. The additive effect of non-R gene functional variants may constitute high levels of durable poplar leaf rust resistance. Therefore, these findings are of significance for speeding the genetic improvement of this long-lived, economically important organism.

A single, plastic population of Mycosphaerella pinodes causes ascochyta blight on winter and spring peas (Pisum sativum) in France

Plant diseases are caused by pathogen populations continuously subjected to evolutionary forces (genetic flow, selection, and recombination). Ascochyta blight, caused by Mycosphaerella pinodes, is one of the most damaging necrotrophic pathogens of field peas worldwide. In France, both winter and spring peas are cultivated. Although these crops overlap by about 4 months (March to June), primary Ascochyta blight infections are not synchronous on the two crops. This suggests that the disease could be due to two different M. pinodes populations, specialized on either winter or spring pea. To test this hypothesis, 144 pathogen isolates were collected in the field during the winter and spring growing seasons in Rennes (western France), and all the isolates were genotyped using amplified fragment length polymorphism (AFLP) markers. Furthermore, the pathogenicities of 33 isolates randomly chosen within the collection were tested on four pea genotypes (2 winter and 2 spring types) grown under three climatic regimes, simulating winter, late winter, and spring conditions. M. pinodes isolates from winter and spring peas were genetically polymorphic but not differentiated according to the type of cultivars. Isolates from winter pea were more pathogenic than isolates from spring pea on hosts raised under winter conditions, while isolates from spring pea were more pathogenic than those from winter pea on plants raised under spring conditions. These results show that disease developed on winter and spring peas was initiated by a single population of M. pinodes whose pathogenicity is a plastic trait modulated by the physiological status of the host plant.

Variation in infectivity and aggressiveness in space and time in wild host-pathogen systems: causes and consequences

Variation in host resistance and in the ability of pathogens to infect and grow (i.e. pathogenicity) is important as it provides the raw material for antagonistic (co)evolution and therefore underlies risks of disease spread, disease evolution and host shifts. Moreover, the distribution of this variation in space and time may inform us about the mode of coevolutionary selection (arms race vs. fluctuating selection dynamics) and the relative roles of G × G interactions, gene flow, selection and genetic drift in shaping coevolutionary processes. Although variation in host resistance has recently been reviewed, little is known about overall patterns in the frequency and scale of variation in pathogenicity, particularly in natural systems. Using 48 studies from 30 distinct host-pathogen systems, this review demonstrates that variation in pathogenicity is ubiquitous across multiple spatial and temporal scales. Quantitative analysis of a subset of extensively studied plant-pathogen systems shows that the magnitude of within-population variation in pathogenicity is large relative to among-population variation and that the distribution of pathogenicity partly mirrors the distribution of host resistance. At least part of the variation in pathogenicity found at a given spatial scale is adaptive, as evidenced by studies that have examined local adaptation at scales ranging from single hosts through metapopulations to entire continents and - to a lesser extent - by comparisons of pathogenicity with neutral genetic variation. Together, these results support coevolutionary selection through fluctuating selection dynamics. We end by outlining several promising directions for future research.

Variation and selection of quantitative traits in plant pathogens

The first section presents the quantitative traits of pathogenicity that are most commonly measured by plant pathologists, how the expression of those traits is influenced by environmental factors, and why the traits must be taken into account for understanding pathogen evolution in agricultural systems. Particular attention is given to the shared genetic control of these traits by the host and the pathogen. Next, the review discusses how quantitative traits account for epidemic development and how they can be related to pathogen fitness. The main constraints that influence the evolution of quantitative traits in pathogen populations are detailed. Finally, possible directions for research on the management of pathogen virulence (as defined by evolutionists) and host quantitative resistance are presented. The review evaluates how the theoretical corpus developed by epidemiologists and evolutionists may apply to plant pathogens in the context of agriculture. The review also analyzes theoretical papers and compares the modeling hypotheses to the biological characteristics of plant pathogens.

Obligate biotrophy features unraveled by the genomic analysis of rust fungi

Rust fungi are some of the most devastating pathogens of crop plants. They are obligate biotrophs, which extract nutrients only from living plant tissues and cannot grow apart from their hosts. Their lifestyle has slowed the dissection of molecular mechanisms underlying host invasion and avoidance or suppression of plant innate immunity. We sequenced the 101-Mb genome of Melampsora larici-populina, the causal agent of poplar leaf rust, and the 89-Mb genome of Puccinia graminis f. sp. tritici, the causal agent of wheat and barley stem rust. We then compared the 16,399 predicted proteins of M. larici-populina with the 17,773 predicted proteins of P. graminis f. sp tritici. Genomic features related to their obligate biotrophic lifestyle include expanded lineage-specific gene families, a large repertoire of effector-like small secreted proteins, impaired nitrogen and sulfur assimilation pathways, and expanded families of amino acid and oligopeptide membrane transporters. The dramatic up-regulation of transcripts coding for small secreted proteins, secreted hydrolytic enzymes, and transporters in planta suggests that they play a role in host infection and nutrient acquisition. Some of these genomic hallmarks are mirrored in the genomes of other microbial eukaryotes that have independently evolved to infect plants, indicating convergent adaptation to a biotrophic existence inside plant cells.

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.

Photosynthetic and respiratory changes in leaves of poplar elicited by rust infection

Poplars are challenged by a wide range of pathogens during their lifespan, and have an innate immunity system that activates defence responses to restrict pathogen growth. Large-scale expression studies of poplar-rust interactions have shown concerted transcriptional changes during defence responses, as in other plant pathosystems. Detailed analysis of expression profiles of metabolic pathways in these studies indicates that photosynthesis and respiration are also important components of the poplar response to rust infection. This is consistent with our current understanding of plant pathogen interactions as defence responses impose substantive demands for resources and energy that are met by reorganization of primary metabolism. This review applies the results of poplar transcriptome analyses to current research describing how plants divert energy from plant primary metabolism for resistance mechanisms.

Molecular and histochemical characterisation of two distinct poplar Melampsora leaf rust pathosystems

In this study, we compared interactions of two Melampsora foliar rust species with poplar, which resulted in either limited or abundant pathogen proliferation. In the pathosystem exhibiting limited pathogen growth, a defence response was observed after invasion of poplar leaf tissues by the biotroph, with late and clear production of reactive oxygen species (ROS) and other products. Characterisation of the histological, biochemical and transcriptional events occurring in both pathosystems showed striking similarity with components of plant defence reactions observed during qualitative resistance. Key components associated with development of an active defence response, such as up-regulation of pathogenesis-related (PR) genes, were observed during infection. Moreover, the time course and strength of gene induction appear to be critical determinants for the outcome of the tree-pathogen interaction. This work provides basic biochemical characterisation and expression data for the study of so-called partial resistance in the poplar-rust pathosystem, which is also applicable to other plant-pathogen interactions resulting in quantitative disease resistance.

Population structure of Potebniamyces pyri in the U.S. Pacific Northwest and evidence of outcrossing inferred with sequence-characterized amplified region markers

Potebniamyces pyri is the cause of Phacidiopycnis rot of d'Anjou pear, which is grown primarily in Washington and Oregon. To estimate the population structure of P. pyri, 146 single-spore isolates were sampled from five major pear-production areas and scored for variation at eight sequence-characterized amplified region (SCAR) loci. Significant genetic differentiation was detected among five subpopulations and a total of 54 multilocus genotypes were identified, with significant genotypic diversity in each subpopulation. No genotype was shared among more than three subpopulations. To estimate the relationship between phenotype and multilocus SCAR genotype, four to five representative isolates of each dominant SCAR genotype in each subpopulation were assayed for growth rate on oatmeal agar and virulence on d'Anjou pear fruit. Significant differences in daily growth rates and virulence were detected among genotypes; however, genotype was not predictive of virulence. To assess the mating system of the pathogen, 10 ascospores were sampled from each of 20 apothecia from a commercial orchard and scored for five SCAR markers. Segregation of alleles at one or more SCAR loci was detected among 18 of 20 ascospore progeny sets, indicating that P. pyri is likely a heterothallic fungus with a predominantly outcrossing mating system.

Transcriptome profiling in hybrid poplar following interactions with Melampsora rust fungi

In natural conditions, plants are subjected to a combination of biotic stresses and often have to cope with simultaneous pathogen infections. In this report, we aim to understand the global transcriptional response of hybrid poplar NM6 (Populus nigra x P. maximowiczii) to infection by two biotrophic Melampsora fungi, Melampsora larici-populina and M. medusae f. sp. deltoidae. These pathogens triggered different responses after inoculation of poplar leaves. Transcript profiling using the GeneChip Poplar Genome Array revealed a total of 416 differentially expressed transcripts whose expression level was > or = twofold relative to controls. Interestingly, approximately half of the differentially expressed genes in infected leaves showed altered expression following interaction with either of the Melampsora spp. We also infected poplar leaves simultaneously with both Melampsora spp. to investigate potential interaction between the responses to the individual pathogens during a mixed infection. For this mixed inoculation, the number of differentially expressed transcripts increased to 648 and our analysis showed that infection with both fungi also induced a common set of genes. The genes induced after Melampsora spp. infection were mainly related to primary and secondary metabolic processes, cell-wall reinforcement and lignification, defense and stress-related mechanisms, and signal perception and transduction.

Genetic structure of the poplar rust fungus Melampsora larici-populina: evidence for isolation by distance in Europe and recent founder effects overseas

Dispersal has a great impact on the genetic structure of populations, but remains difficult to estimate by direct measures. In particular, gradual and stochastic dispersal are often difficult to assess and to distinguish, although they have different evolutionary consequences. Plant pathogens, especially rust fungi, are suspected to display both dispersal modes, though on different spatial scales. In this study, we inferred dispersal capacities of the poplar rust fungus Melampsora larici-populina by examining the genetic diversity and structure of 13 populations from eight European and two overseas countries in the Northern hemisphere. M. larici-populina was sampled from both cultivated hybrid poplars and on the wild host, Populus nigra. The populations were analyzed with 11 microsatellite and 8 virulence markers. Although isolates displayed different virulence profiles according to the host plant, neutral markers revealed little population differentiation with respect to the type of host. This suggests an absence of reproductive isolation between populations sampled from cultivated and wild poplars. Conversely, studying the relationship between geographic and genetic structure allowed us to distinguish between isolation by distance (IBD) patterns and long distance dispersal (LDD) events. The European populations exhibited a significant IBD pattern, suggesting a regular and gradual dispersal of the pathogen over this spatial scale. Nonetheless, the genetic differentiation between these populations was low, suggesting an important gene flow on a continental scale. The two overseas populations from Iceland and Canada were shown to result from rare LDD events, and exhibited signatures of strong founder effects. Furthermore, the high genetic differentiation between both populations suggested that these two recent introductions were independent. This study illustrated how the proper use of population genetics methods can enable contrasted dispersal modes to be revealed.

Landscape epidemiology of plant diseases

Many agricultural landscapes are characterized by a high degree of heterogeneity and fragmentation. Landscape ecology focuses on the influence of habitat heterogeneity in space and time on ecological processes. Landscape epidemiology aims at applying concepts and approaches originating from landscape ecology to the study of pathogen dynamics at the landscape scale. However, despite the strong influence that the landscape properties may have on the spread of plant diseases, landscape epidemiology has still received little attention from plant pathologists. Some recent methodological and technological progress provides new and powerful tools to describe and analyse the spatial patterns of host-pathogen interactions. Here, we review some important topics in plant pathology that may benefit from a landscape perspective. These include the influence of: landscape composition on the global inoculum pressure; landscape heterogeneity on pathogen dynamics; landscape structure on pathogen dispersal; and landscape properties on the emergence of pathogens and on their evolution.

Transcript profiling of poplar leaves upon infection with compatible and incompatible strains of the foliar rust Melampsora larici-populina

To understand key processes governing defense mechanisms in poplar (Populus spp.) upon infection with the rust fungus Melampsora larici-populina, we used combined histological and molecular techniques to describe the infection of Populus trichocarpa x Populus deltoides 'Beaupré' leaves by compatible and incompatible fungal strains. Striking differences in host-tissue infection were observed after 48-h postinoculation (hpi) between compatible and incompatible interactions. No reactive oxygen species production could be detected at infection sites, while a strong accumulation of monolignols occurred in the incompatible interaction after 48 hpi, indicating a late plant response once the fungus already penetrated host cells to form haustorial infection structures. P. trichocarpa whole-genome expression oligoarrays and sequencing of cDNAs were used to determine changes in gene expression in both interactions at 48 hpi. Temporal expression profiling of infection-regulated transcripts was further compared by cDNA arrays and reverse transcription-quantitative polymerase chain reaction. Among 1,730 significantly differentially expressed transcripts in the incompatible interaction, 150 showed an increase in concentration > or =3-fold, whereas 62 were decreased by > or =3-fold. Regulated transcripts corresponded to known genes targeted by R genes in plant pathosystems, such as inositol-3-P synthase, glutathione S-transferases, and pathogenesis-related proteins. However, the transcript showing the highest rust-induced up-regulation encodes a putative secreted protein with no known function. In contrast, only a few transcripts showed an altered expression in the compatible interaction, suggesting a delay in defense response between incompatible and compatible interactions in poplar. This comprehensive analysis of early molecular responses of poplar to M. larici-populina infection identified key genes that likely contain the fungus proliferation in planta. Sequences from the SSH library described in this article can be retrieved in GenBank under accession numbers CT 027996 to CT 029994 and CT 033829.