Global Hierarchical Gene Diversity Analysis Suggests the Fertile Crescent Is Not the Center of Origin of the Barley Scald Pathogen Rhynchosporium secalis

Assessment and modeling using machine learning of resistance to scald (Rhynchosporium commune) in two specific barley genetic resources subsets

Barley production worldwide is limited by several abiotic and biotic stresses and breeding of highly productive and adapted varieties is key to overcome these challenges. Leaf scald, caused by Rhynchosporium commune is a major disease of barley that requires the identification of novel sources of resistance. In this study two subsets of genebank accessions were used: one extracted from the Reference set developed within the Generation Challenge Program (GCP) with 191 accessions, and the other with 101 accessions selected using the filtering approach of the Focused Identification of Germplasm Strategy (FIGS). These subsets were evaluated for resistance to scald at the seedling stage under controlled conditions using two Moroccan isolates, and at the adult plant stage in Ethiopia and Morocco. The results showed that both GCP and FIGS subsets were able to identify sources of resistance to leaf scald at both plant growth stages. In addition, the test of independence and goodness of fit showed that FIGS filtering approach was able to capture higher percentages of resistant accessions compared to GCP subset at the seedling stage against two Moroccan scald isolates, and at the adult plant stage against four field populations of Morocco and Ethiopia, with the exception of Holetta nursery 2017. Furthermore, four machine learning models were tuned on training sets to predict scald reactions on the test sets based on diverse metrics (accuracy, specificity, and Kappa). All models efficiently identified resistant accessions with specificities higher than 0.88 but showed different performances between isolates at the seedling and to field populations at the adult plant stage. The findings of our study will help in fine-tuning FIGS approach using machine learning for the selection of best-bet subsets for resistance to scald disease from the large number of genebank accessions.

Population structure of Erysiphe necator on domesticated and wild vines in the Middle East raises questions on the origin of the grapevine powdery mildew pathogen

Plant pathogens usually originate and diversify in geographical regions where hosts and pathogens co-evolve. Erysiphe necator, the causal agent of grape powdery mildew, is a destructive pathogen of grapevines worldwide. Although Eastern US is considered the centre of origin and diversity of E. necator, previous reports on resistant native wild and domesticated Asian grapevines suggest Asia as another possible origin of the pathogen. By using multi-locus sequencing, microsatellites and a novel application of amplicon sequencing (AmpSeq), we show that the population of E. necator in Israel is composed of three genetic groups: Groups A and B that are common worldwide, and a new group IL, which is genetically differentiated from any known group in Europe and Eastern US. Group IL showed distinguished ecological characteristics: it was dominant on wild and traditional vines (95%); its abundance increased along the season; and was more aggressive than A and B isolates on both wild and domesticated vines. The low genetic diversity within group IL suggests that it has invaded Israel from another origin. Therefore, we suggest that the Israeli E. necator population was founded by at least two invasions, of which one could be from a non-East American source, possibly from Asian origin.

Recent insights into barley and Rhynchosporium commune interactions

Rhynchosporium commune is the causal pathogen of scald in barley (Hordeum vulgare), a foliar disease that can reduce yield by up to 40% in susceptible cultivars. R. commune is found worldwide in all temperate growing regions and is regarded as one of the most economically important barley pathogens. It is a polycyclic pathogen with the ability to rapidly evolve new virulent strains in response to resistance genes deployed in commercial cultivars. Hence, introgression and pyramiding of different loci for resistance (qualitative or quantitative) through marker-assisted selection is an effective way to improve scald resistance in barley. This review summarizes all 148 resistance quantitative trait loci reported at the date of submission of this review and projects them onto the barley physical map, where it is clear many loci co-locate on chromosomes 3H and 7H. We have summarized the major named resistance loci and reiterated the renaming of Rrs15 (CI8288) to Rrs17. This review provides a comprehensive resource for future discovery and breeding efforts of qualitative and quantitative scald resistance loci.

Fine mapping QSc.VR4, an effective and stable scald resistance locus in barley (Hordeum vulgare L.), to a 0.38-Mb region enriched with LRR-RLK and GLP genes

An effective and stable quantitative resistance locus, QSc.VR4, was fine mapped, characterized and physically anchored to the short arm of 4H, conferring adult plant resistance to the fungus Rhynchosporium commune in barley. Scald caused by Rhynchosporium commune is one of the most destructive barley diseases worldwide. Accumulation of adult plant resistance (APR) governed by multiple resistance alleles is predicted to be effective and long-lasting against a broad spectrum of pathotypes. However, the molecular mechanisms that control APR remain poorly understood. Here, quantitative trait loci (QTL) analysis of APR and fine mapping were performed on five barley populations derived from a common parent Vlamingh, which expresses APR to scald. Two QTLs, designated QSc.VR4 and QSc.BR7, were detected from a cross between Vlamingh and Buloke. Our data confirmed that QSc.VR4 is an effective and stable APR locus, residing on the short arm of chromosome 4H, and QSc.BR7 derived from Buloke may be an allele of reported Rrs2. High-resolution fine mapping revealed that QSc.VR4 is located in a 0.38 Mb genomic region between InDel markers 4H2282169 and 4H2665106. The gene annotation analysis and sequence comparison suggested that a gene cluster containing two adjacent multigene families encoding leucine-rich repeat receptor kinase-like proteins (LRR-RLKs) and germin-like proteins (GLPs), respectively, is likely contributing to scald resistance. Adult plant resistance (APR) governed by QSc.VR4 may confer partial levels of resistance to the fungus Rhynchosporium commune and, furthermore, be an important resource for gene pyramiding that may contribute broad-based and more durable resistance.

Characterisation of barley landraces from Syria and Jordan for resistance to rhynchosporium and identification of diagnostic markers for Rrs1Rh4

Diagnostic markers for Rrs1_Rh4 have been identified by testing for associations between SNPs within the Rrs1 interval in 150 barley genotypes and their resistance to Rhynchosporium commune isolates recognised by lines containing Rrs1. Rhynchosporium or barley scald, caused by the destructive fungal pathogen Rhynchosporium commune, is one of the most economically important diseases of barley in the world. Barley landraces from Syria and Jordan demonstrated high resistance to rhynchosporium in the field. Genotyping of a wide range of barley cultivars and landraces, including known sources of different Rrs1 genes/alleles, across the Rrs1 interval, followed by association analysis of this genotypic data with resistance phenotypes to R. commune isolates recognised by Rrs1, allowed the identification of diagnostic markers for Rrs1_Rh4. These markers are specific to Rrs1_Rh4 and do not detect other Rrs1 genes/alleles. The Rrs1_Rh4 diagnostic markers represent a resource that can be exploited by breeders for the sustainable deployment of varietal resistance in new cultivars. Thirteen out of the 55 most resistant Syrian and Jordanian landraces were shown to contain markers specific to Rrs1_Rh4. One of these lines came from Jordan, with the remaining 12 lines from different locations in Syria. One of the Syrian landraces containing Rrs1_Rh4 was also shown to have Rrs2. The remaining landraces that performed well against rhynchosporium in the field are likely to contain other resistance genes and represent an important novel resource yet to be exploited by European breeders.

Characterization of the Barley Net Blotch Pathosystem at the Center of Origin of Host and Pathogen

Net blotch (NB) is a major disease of barley caused by the fungus Pyrenophora teres f. teres (Ptt), and P. teres f. maculata (Ptm). Ptt and Ptm infect the cultivated crop (Hordeum vulgare) and its wild relatives (H. vulgare ssp. spontaneum and H. murinum ssp. glaucum). The main goal of this research was to study the NB-causing pathogen in the crop center of origin. To address this, we have constructed a Ptt (n = 15) and Ptm (n = 12) collection isolated from three barley species across Israel. Isolates were characterized genetically and phenotypically. Aggressiveness of the isolates was determined based on necrotrophic growth rate on detached leaves of barley. In addition, isolates were genetically characterized by the mating type, followed by phylogenetic analysis, clustering them into seven groups. The analysis showed no significant differentiation of isolates based on either geographic origin, host of origin or form (Ptt vs. Ptm). Nevertheless, there was a significant difference in aggressiveness among the isolates regardless of host species, geographic location or sampling site. Moreover, it was apparent that the isolates derived from wild hosts were more variable in their necrotrophic growth rate, compared to isolates sampled from cultivated hosts, thereby suggesting that NB plays a major role in epidemiology at the center of barley origin where most of the diversity lies. Ptm has significantly higher necrotrophic and saprotrophic growth rates than Ptt, and for both a significant negative correlation was found between light intensity exposure and growth rates.

The emergence of the multi-species NIP1 effector in Rhynchosporium was accompanied by high rates of gene duplications and losses

Plant pathogens secrete effector proteins to manipulate the host and facilitate infection. Cognate hosts trigger strong defence responses upon detection of these effectors. Consequently, pathogens and hosts undergo rapid coevolutionary arms races driven by adaptive evolution of effectors and receptors. Because of their high rate of turnover, most effectors are thought to be species-specific and the evolutionary trajectories are poorly understood. Here, we investigate the necrosis-inducing protein 1 (NIP1) effector in the multihost pathogen genus Rhynchosporium. We retraced the evolutionary history of the NIP1 locus using whole-genome assemblies of 146 strains covering four closely related species. NIP1 orthologues were present in all species but the locus consistently segregated presence-absence polymorphisms suggesting long-term balancing selection. We also identified previously unknown paralogues of NIP1 that were shared among multiple species and showed substantial copy-number variation within R. commune. The NIP1A paralogue was under significant positive selection suggesting that NIP1A is the dominant effector variant coevolving with host immune receptors. Consistent with this prediction, we found that copy number variation at NIP1A had a stronger effect on virulence than NIP1B. Our analyses unravelled the origins and diversification mechanisms of a pathogen effector family shedding light on how pathogens gain adaptive genetic variation.

Resistance to Rhynchosporium commune in a collection of European spring barley germplasm

Association analyses of resistance to Rhynchosporium commune in a collection of European spring barley germplasm detected 17 significant resistance quantitative trait loci. The most significant association was confirmed as Rrs1. Rhynchosporium commune is a fungal pathogen of barley which causes a highly destructive and economically important disease known as rhynchosporium. Genome-wide association mapping was used to investigate the genetic control of host resistance to R. commune in a collection of predominantly European spring barley accessions. Multi-year disease nursery field trials revealed 8 significant resistance quantitative trait loci (QTL), whilst a separate association mapping analysis using historical data from UK national and recommended list trials identified 9 significant associations. The most significant association identified in both current and historical data sources, collocated with the known position of the major resistance gene Rrs1. Seedling assays with R. commune single-spore isolates expressing the corresponding avirulence protein NIP1 confirmed that this locus is Rrs1. These results highlight the significant and continuing contribution of Rrs1 to host resistance in current elite spring barley germplasm. Varietal height was shown to be negatively correlated with disease severity, and a resistance QTL was identified that co-localised with the semi-dwarfing gene sdw1, previously shown to contribute to disease escape. The remaining QTL represent novel resistances that are present within European spring barley accessions. Associated markers to Rrs1 and other resistance loci, identified in this study, represent a set of tools that can be exploited by breeders for the sustainable deployment of varietal resistance in new cultivars.

Genome-Wide Detection of Genes Under Positive Selection in Worldwide Populations of the Barley Scald Pathogen

Coevolution between hosts and pathogens generates strong selection pressures to maintain resistance and infectivity, respectively. Genomes of plant pathogens often encode major effect loci for the ability to successfully infect specific host genotypes. Hence, spatial heterogeneity in host genotypes coupled with abiotic factors could lead to locally adapted pathogen populations. However, the genetic basis of local adaptation is poorly understood. Rhynchosporium commune, the pathogen causing barley scald disease, interacts at least partially in a gene-for-gene manner with its host. We analyzed global field populations of 125 R. commune isolates to identify candidate genes for local adaptation. Whole genome sequencing data showed that the pathogen is subdivided into three genetic clusters associated with distinct geographic and climatic regions. Using haplotype-based selection scans applied independently to each genetic cluster, we found strong evidence for selective sweeps throughout the genome. Comparisons of loci under selection among clusters revealed little overlap, suggesting that ecological differences associated with each cluster led to variable selection regimes. The strongest signals of selection were found predominantly in the two clusters composed of isolates from Central Europe and Ethiopia. The strongest selective sweep regions encoded protein functions related to biotic and abiotic stress responses. Selective sweep regions were enriched in genes encoding functions in cellular localization, protein transport activity, and DNA damage responses. In contrast to the prevailing view that a small number of gene-for-gene interactions govern plant pathogen evolution, our analyses suggest that the evolutionary trajectory is largely determined by spatially heterogeneous biotic and abiotic selection pressures.

Genetic structure of Cercospora beticola populations on Beta vulgaris in New York and Hawaii

Cercospora leaf spot (CLS), caused by Cercospora beticola, is a major disease of Beta vulgaris worldwide. No sexual stage is known for C. beticola but in its asexual form it overwinters on infected plant debris as pseudostromata, and travels short distances by rain splash-dispersed conidiospores. Cercospora beticola infects a broad range of host species and may be seedborne. The relative contribution of these inoculum sources to CLS epidemics on table beet is not well understood. Pathogen isolates collected from table beet, Swiss chard and common lambsquarters in mixed-cropping farms and monoculture fields in New York and Hawaii, USA, were genotyped (n = 600) using 12 microsatellite markers. All isolates from CLS symptoms on lambsquarters were identified as C. chenopodii. Sympatric populations of C. beticola derived from Swiss chard and table beet were not genetically differentiated. Results suggested that local (within field) inoculum sources may be responsible for the initiation of CLS epidemics in mixed-cropping farms, whereas external sources of inoculum may be contributing to CLS epidemics in the monoculture fields in New York. New multiplex PCR assays were developed for mating-type determination for C. beticola. Implications of these findings for disease management are discussed.

Population Genetic Structure of Venturia effusa, Cause of Pecan Scab, in the Southeastern United States

Venturia effusa is the most important pathogen of pecan in the southeastern United States. Little information exists on the population biology and genetic diversity of the pathogen. A hierarchical sampling of 784 isolates from 63 trees in 11 pecan orchards in the southeastern United States were screened against a set of 30 previously characterized microsatellite markers. Populations were collected from Georgia (n = 2), Florida (n = 1), Alabama (n = 2), Mississippi (n = 1), Louisiana (n = 1), Illinois (n = 1), Oklahoma (n = 1), Texas (n = 1), and Kansas (n = 1). Clonality was low in all orchard populations (≤10.1% of isolates), and there were consistently high levels of genotypic diversity (Shannon-Weiner's index = 3.49 to 4.59) and gene diversity (Nei's measure = 0.513 to 0.713). Analysis of molecular variance showed that, although 81% of genetic diversity occurred at the scale of the individual tree, 16% occurred between orchards and only 3% between trees within orchards. All populations could be differentiated from each other (P = 0.01), and various cluster analyses indicated that some populations were more closely related compared with other pairs of populations. This is indicative of some limited population differentiation in V. effusa in the southeastern United States. Bayesian and nearest-neighbor methods suggested eight clusters, with orchards from Georgia and Florida being grouped together. A minimum spanning tree of all 784 isolates also indicated some isolate identification with source population. Linkage disequilibrium was detected in all but one population (Kansas), although 8 of the 11 populations had <20% of loci at disequilibrium. A Mantel test demonstrated a relationship between physical and genetic distance between populations (Z = 11.9, r = 0.559, P = 0.001). None of the populations were at mutation-drift equilibrium. All but 3 of the 11 populations had a deficiency of gene diversity compared with that expected at mutation-drift equilibrium (indicating population expansion); the remaining populations had an excess of gene diversity compared with that expected at mutation-drift equilibrium (indicating a recent bottleneck). These observations are consistent with the known history of pecan and pecan scab, which is that V. effusa became an issue on cultivated pecan in the last approximately 120 years (recent population expansion). Recently reported mating type genes and the sexual stage of this fungus may help explain the observed population characteristics, which bear a strong resemblance to those of other well-characterized sexually reproducing ascomycete pathogens.

Comparative genomics to explore phylogenetic relationship, cryptic sexual potential and host specificity of Rhynchosporium species on grasses

BACKGROUND

The Rhynchosporium species complex consists of hemibiotrophic fungal pathogens specialized to different sweet grass species including the cereal crops barley and rye. A sexual stage has not been described, but several lines of evidence suggest the occurrence of sexual reproduction. Therefore, a comparative genomics approach was carried out to disclose the evolutionary relationship of the species and to identify genes demonstrating the potential for a sexual cycle. Furthermore, due to the evolutionary very young age of the five species currently known, this genus appears to be well-suited to address the question at the molecular level of how pathogenic fungi adapt to their hosts.

RESULTS

The genomes of the different Rhynchosporium species were sequenced, assembled and annotated using ab initio gene predictors trained on several fungal genomes as well as on Rhynchosporium expressed sequence tags. Structures of the rDNA regions and genome-wide single nucleotide polymorphisms provided a hypothesis for intra-genus evolution. Homology screening detected core meiotic genes along with most genes crucial for sexual recombination in ascomycete fungi. In addition, a large number of cell wall-degrading enzymes that is characteristic for hemibiotrophic and necrotrophic fungi infecting monocotyledonous hosts were found. Furthermore, the Rhynchosporium genomes carry a repertoire of genes coding for polyketide synthases and non-ribosomal peptide synthetases. Several of these genes are missing from the genome of the closest sequenced relative, the poplar pathogen Marssonina brunnea, and are possibly involved in adaptation to the grass hosts. Most importantly, six species-specific genes coding for protein effectors were identified in R. commune. Their deletion yielded mutants that grew more vigorously in planta than the wild type.

CONCLUSION

Both cryptic sexuality and secondary metabolites may have contributed to host adaptation. Most importantly, however, the growth-retarding activity of the species-specific effectors suggests that host adaptation of R. commune aims at extending the biotrophic stage at the expense of the necrotrophic stage of pathogenesis. Like other apoplastic fungi Rhynchosporium colonizes the intercellular matrix of host leaves relatively slowly without causing symptoms, reminiscent of the development of endophytic fungi. Rhynchosporium may therefore become an object for studying the mutualism-parasitism transition.

QTL mapping of temperature sensitivity reveals candidate genes for thermal adaptation and growth morphology in the plant pathogenic fungus Zymoseptoria tritici

Different thermal environments impose strong, differential selection on populations, leading to local adaptation, but the genetic basis of thermal adaptation is poorly understood. We used quantitative trait locus (QTL) mapping in the fungal wheat pathogen Zymoseptoria tritici to study the genetic architecture of thermal adaptation and identify candidate genes. Four wild-type strains originating from the same thermal environment were crossed to generate two mapping populations with 263 (cross 1) and 261 (cross 2) progeny. Restriction site-associated DNA sequencing was used to genotype 9745 (cross 1) and 7333 (cross 2) single-nucleotide polymorphism markers segregating within the mapping population. Temperature sensitivity was assessed using digital image analysis of colonies growing at two different temperatures. We identified four QTLs for temperature sensitivity, with unique QTLs found in each cross. One QTL had a logarithm of odds score >11 and contained only six candidate genes, including PBS2, encoding a mitogen-activated protein kinase kinase associated with low temperature tolerance in Saccharomyces cerevisiae. This and other QTLs showed evidence for pleiotropy among growth rate, melanization and growth morphology, suggesting that many traits can be correlated with thermal adaptation in fungi. Higher temperatures were highly correlated with a shift to filamentous growth among the progeny in both crosses. We show that thermal adaptation has a complex genetic architecture, with natural populations of Z. tritici harboring significant genetic variation for this trait. We conclude that Z. tritici populations have the potential to adapt rapidly to climate change and expand into new climatic zones.

Are PECTIN ESTERASE INHIBITOR Genes Involved in Mediating Resistance to Rhynchosporium commune in Barley?

A family of putative PECTIN ESTERASE INHIBITOR (PEI) genes, which were detected in the genomic region co-segregating with the resistance gene Rrs2 against scald caused by Rhynchosporium commune in barley, were characterized and tested for their possible involvement in mediating resistance to the pathogen by complementation and overexpression analysis. The sequences of the respective genes were derived from two BAC contigs originating from the susceptible cultivar ‘Morex’. For the genes HvPEI2, HvPEI3, HvPEI4 and HvPEI6, specific haplotypes for 18 resistant and 23 susceptible cultivars were detected after PCR-amplification and haplotype-specific CAPS-markers were developed. None of the tested candidate genes HvPEI2, HvPEI3 and HvPEI4 alone conferred a high resistance level in transgenic over-expression plants, though an improvement of the resistance level was observed especially with OE-lines for gene HvPEI4. These results do not confirm but also do not exclude an involvement of the PEI gene family in the response to the pathogen. A candidate for the resistance gene Rrs2 could not be identified yet. It is possible that Rrs2 is a PEI gene or another type of gene which has not been detected in the susceptible cultivar ‘Morex’ or the full resistance reaction requires the presence of several PEI genes.

The influence of genetic drift and selection on quantitative traits in a plant pathogenic fungus

Genetic drift and selection are ubiquitous evolutionary forces acting to shape genetic variation in populations. While their relative importance has been well studied in plants and animals, less is known about their relative importance in fungal pathogens. Because agro-ecosystems are more homogeneous environments than natural ecosystems, stabilizing selection may play a stronger role than genetic drift or diversifying selection in shaping genetic variation among populations of fungal pathogens in agro-ecosystems. We tested this hypothesis by conducting a Q_ST/F_ST analysis using agricultural populations of the barley pathogen Rhynchosporium commune. Population divergence for eight quantitative traits (Q_ST) was compared with divergence at eight neutral microsatellite loci (F_ST) for 126 pathogen strains originating from nine globally distributed field populations to infer the effects of genetic drift and types of selection acting on each trait. Our analyses indicated that five of the eight traits had Q_ST values significantly lower than F_ST, consistent with stabilizing selection, whereas one trait, growth under heat stress (22°C), showed evidence of diversifying selection and local adaptation (Q_ST>F_ST). Estimates of heritability were high for all traits (means ranging between 0.55–0.84), and average heritability across traits was negatively correlated with microsatellite gene diversity. Some trait pairs were genetically correlated and there was significant evidence for a trade-off between spore size and spore number, and between melanization and growth under benign temperature. Our findings indicate that many ecologically and agriculturally important traits are under stabilizing selection in R. commune and that high within-population genetic variation is maintained for these traits.

Rhynchosporium commune: a persistent threat to barley cultivation

Rhynchosporium commune is a haploid fungus causing scald or leaf blotch on barley, other Hordeum spp. and Bromus diandrus.

TAXONOMY

Rhynchosporium commune is an anamorphic Ascomycete closely related to the teleomorph Helotiales genera Oculimacula and Pyrenopeziza.

DISEASE SYMPTOMS

Rhynchosporium commune causes scald-like lesions on leaves, leaf sheaths and ears. Early symptoms are generally pale grey oval lesions. With time, the lesions acquire a dark brown margin with the centre of the lesion remaining pale green or pale brown. Lesions often merge to form large areas around which leaf yellowing is common. Infection frequently occurs in the leaf axil, which can lead to chlorosis and eventual death of the leaf.

LIFE CYCLE

Rhynchosporium commune is seed borne, but the importance of this phase of the disease is not fully understood. Debris from previous crops and volunteers, infected from the stubble from previous crops, are considered to be the most important sources of the disease. Autumn-sown crops can become infected very soon after sowing. Secondary spread of disease occurs mainly through splash dispersal of conidia from infected leaves. Rainfall at the stem extension growth stage is the major environmental factor in epidemic development. DETECTION AND QUANTIFICATION: Rhynchosporium commune produces unique beak-shaped, one-septate spores both on leaves and in culture. The development of a specific polymerase chain reaction (PCR) and, more recently, quantitative PCR (qPCR) has allowed the identification of asymptomatic infection in seeds and during the growing season.

DISEASE CONTROL

The main measure for the control of R. commune is the use of fungicides with different modes of action, in combination with the use of resistant cultivars. However, this is constantly under review because of the ability of the pathogen to adapt to host plant resistance and to develop fungicide resistance.

Pathogenicity Spectra and Screening for Resistance in Barley Against Tunisian Pyrenophora teres f. teres

This work aimed to determine patterns of pathogenicity in Pyrenophora teres f. teres and to identify potentially effective resistance sources that could be used as breeding material to control net blotch in Tunisia. Extensive pathogenic variability was detected in 85 isolates of P. teres causing net blotch of barley in Tunisia. Based on unweighted pair-group method with arithmetic averaging clustering and mean disease rating scores, three distinct virulence groups were identified. The isolates were classified into 23 pathotypes. Pathogenic variability within the groups was higher than that between the groups, a finding that can guide a rational choice of isolates for screening lines as part of a breeding program. Conversely, studying the relationship between geographic and pathotypic structure allowed us to detect a significant isolation by distance pattern, suggesting a regular and gradual dispersal of the pathogen over this spatial scale. Using specific resistance properties of individual barley genotypes as virulence markers, all the differential barley genotypes were shown to be distinct, and no single source of resistance was totally effective against all isolates.

Necrosis-inducing proteins of Rhynchosporium commune, effectors in quantitative disease resistance

The barley pathogen Rhynchosporium commune secretes necrosis-inducing proteins NIP1, NIP2, and NIP3. Expression analysis revealed that NIP1 transcripts appear to be present in fungal spores already, whereas NIP2 and NIP3 are synthesized after inoculation of host plants. To assess the contribution of the three effector proteins to disease development, deletion mutants were generated. The development of these fungal mutants on four barley cultivars was quantified in comparison with that of the parent wild-type strain and with two fungal strains failing to secrete an "active" NIP1 avirulence protein, using quantitative polymerase chain reaction as well as microscopic imaging after fungal green fluorescent protein tagging. The impact of the three deletions varied quantitatively depending on the host genotype, suggesting that the activities of the fungal effectors add up to produce stronger growth patterns and symptom development. Alternatively, recognition events of differing intensities may be converted into defense gene expression in a quantitative manner.

Presence and functionality of mating type genes in the supposedly asexual filamentous fungus Aspergillus oryzae

The potential for sexual reproduction in Aspergillus oryzae was assessed by investigating the presence and functionality of MAT genes. Previous genome studies had identified a MAT1-1 gene in the reference strain RIB40. We now report the existence of a complementary MAT1-2 gene and the sequencing of an idiomorphic region from A. oryzae strain AO6. This allowed the development of a PCR diagnostic assay, which detected isolates of the MAT1-1 and MAT1-2 genotypes among 180 strains assayed, including industrial tane-koji isolates. Strains used for sake and miso production showed a near-1:1 ratio of the MAT1-1 and MAT1-2 mating types, whereas strains used for soy sauce production showed a significant bias toward the MAT1-2 mating type. MAT1-1 and MAT1-2 isogenic strains were then created by genetic manipulation of the resident idiomorph, and gene expression was compared by DNA microarray and quantitative real-time PCR (qRT-PCR) methodologies under conditions in which MAT genes were expressed. Thirty-three genes were found to be upregulated more than 10-fold in either the MAT1-1 host strain or the MAT1-2 gene replacement strain relative to each other, showing that both the MAT1-1 and MAT1-2 genes functionally regulate gene expression in A. oryzae in a mating type-dependent manner, the first such report for a supposedly asexual fungus. MAT1-1 expression specifically upregulated an α-pheromone precursor gene, but the functions of most of the genes affected were unknown. The results are consistent with a heterothallic breeding system in A. oryzae, and prospects for the discovery of a sexual cycle are discussed.

Population-based resequencing analysis of wild and cultivated barley revealed weak domestication signal of selection and bottleneck in the Rrs2 scald resistance gene region

Many plant disease resistance (R) genes have been cloned, but the potential of utilizing these plant R-gene genomic resources for genetic inferences of plant domestication history remains unexplored. A population-based resequencing analysis of the genomic region near the Rrs2 scald resistance gene was made in 51 accessions of wild and cultivated barley from 41 countries. Fifteen primer pairs were designed to sample the genomic region with a total length of 10 406 bp. More nucleotide diversity was found in wild (π = 0.01846) than cultivated (π = 0.01507) barley samples. Three distinct groups of 29 haplotypes were detected for all 51 samples, and they were well mixed with wild and cultivated barley samples from different countries and regions. The neutrality tests by Tajima's D were not significant, but a significant (P < 0.05) case by Fu and Li's D* and F* was found in the barley cultivar samples. The estimate of selection intensity by K(a)/K(s) was 0.691 in wild barley and 0.580 in cultivated barley. The estimate of the minimum recombination events was 16 in wild barley and 19 in cultivated barley. A coalescence simulation revealed a bottleneck intensity of 1.5 to 2 since barley domestication. Together, the domestication signal in the genomic region was weak both in human selection and domestication bottleneck.

Contrasting introduction scenarios among continents in the worldwide invasion of the banana fungal pathogen Mycosphaerella fijiensis

Reconstructing and characterizing introduction routes is a key step towards understanding the ecological and evolutionary factors underlying successful invasions and disease emergence. Here, we aimed to decipher scenarios of introduction and stochastic demographic events associated with the global spread of an emerging disease of bananas caused by the destructive fungal pathogen Mycosphaerella fijiensis. We analysed the worldwide population structure of this fungus using 21 microsatellites and 8 sequence-based markers on 735 individuals from 37 countries. Our analyses designated South-East Asia as the source of the global invasion and supported the location of the centre of origin of M. fijiensis within this area. We confirmed the occurrence of bottlenecks upon introduction into other continents followed by widespread founder events within continents. Furthermore, this study suggested contrasting introduction scenarios of the pathogen between the African and American continents. While potential signatures of admixture resulting from multiple introductions were detected in America, all the African samples examined seem to descend from a single successful founder event. In combination with historical information, our study reveals an original and unprecedented global scenario of invasion for this recently emerging disease caused by a wind-dispersed pathogen.

Sexual recombinants make a significant contribution to epidemics caused by the wheat pathogen Phaeosphaeria nodorum

We conducted a 2-year mark-release-recapture field experiment to quantify the relative contributions of immigration and sexual and asexual reproduction to epidemics of Stagonospora nodorum blotch caused by Phaeosphaeria nodorum. The epidemic was initiated using nine genetically distinct P. nodorum isolates. Infected plants were sampled four times across two growing seasons. In total, 1,286 isolates were recovered and assayed with 10 microsatellite markers and 1 minisatellite marker. The proportion of isolates having multilocus haplotypes (MLHTs) identical to the inoculated isolates decreased steadily from 86% in the first collection to 25% in the fourth collection. The novel isolates that had different MLHTs compared with the marked inoculants originated through immigration and sexual recombination. By the end of the experiment, nearly three-quarters of the novel isolates originated from sexual recombination. Our results indicate that recombinant offspring and airborne immigrant ascospores can make significant contributions to epidemics of Stagonospora nodorum blotch during a growing season.

Two new species of Rhynchosporium

Rhynchosporium consists of two species, R. secalis and R. orthosporum. Both are pathogens of grasses with R. secalis infecting a variety of Poaceae hosts and R. orthosporum infecting Dactylis glomerata. Phylogenetic analyses of multilocus DNA sequence data on R. secalis isolates originating from cultivated barley, rye, triticale and other grasses, including Agropyron spp., Bromus diandrus and Hordeum spp., resolved the monophyletic groups into three species according to their respective hosts. Host specificity according to phylogenetic lineages was confirmed with pathogenicity studies. Because R. secalis was described first on rye this name is retained for Rhynchosporium isolates infecting rye and triticale. Rhynchosporium isolates infecting cultivated barley and other Hordeum spp. and Bromus diandrus belong to a distinct species, R. commune. Similarly isolates infecting Agropyron spp. represent a distinct species of Rhynchosporium, namely R. agropyri. A PCR-RFLP assay was developed as a rapid tool for species identification of R. secalis and R. commune.

Pathotype and microsatellite analyses reveal new sources of resistance to barley scald in Tunisia

We examined the variation and relationships between pathogenicity and a microsatellite-based haplotype in 79 Tunisian Rhynchosporium secalis isolates that were collected from the most commonly cultivated barley populations in Tunisia, Rihane cv. and local landraces, with the goal of finding genes that might be used to monitor resistance to scald. Isolates could be classified into three distinct virulence groups based on artificial inoculation of 19 differential cultivars with known scald resistance genes. The resistance gene BRR2 carried by the Astrix differential cultivar appeared to be the most effective in Tunisia. Pathotypes sampled from the Rihane host were more virulent than those sampled from local barley landraces. Because some differential cultivars that carried the same resistance genes showed different reaction patterns to 48 of the isolates, we postulated that other unknown resistance gene(s) specific to Tunisian isolates may be prevalent and could be used in Tunisian barley breeding programs. Microsatellite fingerprinting allowed the detection of 11 alleles linked to the virulence and pathogenic identification of 52% of the tested isolates. Thus, microsatellite analysis may provide a rapid tool for pathogen detection, without an inoculation step that requires long incubation periods before ultimate disease assessment.

Fine mapping, physical mapping and development of diagnostic markers for the Rrs2 scald resistance gene in barley

The Rrs2 gene confers resistance to the fungal pathogen Rhynchosporium secalis which causes leaf scald, a major barley disease. The Rrs2 gene was fine mapped to an interval of 0.08 cM between markers 693M6_6 and P1D23R on the distal end of barley chromosome 7HS using an Atlas (resistant) x Steffi (susceptible) mapping population of 9,179 F(2)-plants. The establishment of a physical map of the Rrs2 locus led to the discovery that Rrs2 is located in an area of suppressed recombination within this mapping population. The analysis of 58 barley genotypes revealed a large linkage block at the Rrs2 locus extending over several hundred kb which is present only in Rrs2 carrying cultivars. Due to the lack of recombination in the mapping population and the presence of a Rrs2-specific linkage block, we assume a local chromosomal rearrangement (alien introgression or inversion) in Rrs2 carrying varieties. The variety analysis led to the discovery of eight SNPs which were diagnostic for the Rrs2 phenotype. Based on these SNPs diagnostic molecular markers (CAPS and pyrosequencing markers) were developed which are highly useful for marker-assisted selection in resistance gene pyramiding programmes for Rhynchosporium secalis resistance in barley.

Phylogeographical analyses reveal global migration patterns of the barley scald pathogen Rhynchosporium secalis

A phylogeographical analysis of the scald pathogen Rhynchosporium secalis was conducted using nuclear DNA sequences from two neutral restriction fragment length polymorphism loci and the mating-type idiomorphs. Approximately 500 isolates sampled from more than 60 field populations from five continents were analysed to infer migration patterns and the demographic history of the fungus. Migration rates among continents were generally low, consistent with earlier reports of significant population subdivision among continents. Northern Europe was mainly a source population for global migration. We hypothesize that the pathogen only recently moved out of its centre of origin, resulting in founder populations that are reproductively isolated due to the contemporary absence of long-distance gene flow.

Ustilago maydis populations tracked maize through domestication and cultivation in the Americas

The domestication of crops and the development of agricultural societies not only brought about major changes in human interactions with the environment but also in plants' interactions with the diseases that challenge them. We evaluated the impact of the domestication of maize from teosinte and the widespread cultivation of maize on the historical demography of Ustilago maydis, a fungal pathogen of maize. To determine the evolutionary response of the pathogen's populations, we obtained multilocus genotypes for 1088 U. maydis diploid individuals from two teosinte subspecies in Mexico and from maize in Mexico and throughout the Americas. Results identified five major U. maydis populations: two in Mexico; two in South America; and one in the United States. The two populations in Mexico diverged from the other populations at times comparable to those for the domestication of maize at 6000-10000 years before present. Maize domestication and agriculture enforced sweeping changes in U. maydis populations such that the standing variation in extant pathogen populations reflects evolution only since the time of the crop's domestication.

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.

Rapid speciation following recent host shifts in the plant pathogenic fungus Rhynchosporium

Agriculture played a significant role in increasing the number of pathogen species and in expanding their geographic range during the last 10,000 years. We tested the hypothesis that a fungal pathogen of cereals and grasses emerged at the time of domestication of cereals in the Fertile Crescent and subsequently speciated after adaptation to its hosts. Rhynchosporium secalis, originally described from rye, causes an important disease on barley called scald, although it also infects other species of Hordeum and Agropyron. Phylogenetic analyses based on four DNA sequence loci identified three host-associated lineages that were confirmed by cross-pathogenicity tests. Bayesian analyses of divergence time suggested that the three lineages emerged between approximately 1200 to 3600 years before present (B.P.) with a 95% highest posterior density ranging from 100 to 12,000 years B.P. depending on the implemented clock models. The coalescent inference of demographic history revealed a very recent population expansion for all three pathogens. We propose that Rhynchosporium on barley, rye, and Agropyron host species represent three cryptic pathogen species that underwent independent evolution and ecological divergence by host-specialization. We postulate that the recent emergence of these pathogens followed host shifts. The subsequent population expansions followed the expansion of the cultivated host populations and accompanying expansion of the weedy Agropyron spp. found in fields of cultivated cereals. Hence, agriculture played a major role in the emergence of the scald diseases, the adaptation of the pathogens to new hosts and their worldwide dissemination.

On the origin and spread of the Scab disease of apple: out of central Asia

BACKGROUND

Venturia inaequalis is an ascomycete fungus responsible for apple scab, a disease that has invaded almost all apple growing regions worldwide, with the corresponding adverse effects on apple production. Monitoring and predicting the effectiveness of intervention strategies require knowledge of the origin, introduction pathways, and population biology of pathogen populations. Analysis of the variation of genetic markers using the inferential framework of population genetics offers the potential to retrieve this information.

METHODOLOGY/PRINCIPAL FINDINGS

Here, we present a population genetic analysis of microsatellite variation in 1,273 strains of V. inaequalis representing 28 orchard samples from seven regions in five continents. Analysis of molecular variance revealed that most of the variation (88%) was distributed within localities, which is consistent with extensive historical migrations of the fungus among and within regions. Despite this shallow population structure, clustering analyses partitioned the data set into separate groups corresponding roughly to geography, indicating that each region hosts a distinct population of the fungus. Comparison of the levels of variability among populations, along with coalescent analyses of migration models and estimates of genetic distances, was consistent with a scenario in which the fungus emerged in Central Asia, where apple was domesticated, before its introduction into Europe and, more recently, into other continents with the expansion of apple growing. Across the novel range, levels of variability pointed to multiple introductions and all populations displayed signatures of significant post-introduction increases in population size. Most populations exhibited high genotypic diversity and random association of alleles across loci, indicating recombination both in native and introduced areas.

CONCLUSIONS/SIGNIFICANCE

Venturia inaequalis is a model of invasive phytopathogenic fungus that has now reached the ultimate stage of the invasion process with a broad geographic distribution and well-established populations displaying high genetic variability, regular sexual reproduction, and demographic expansion.

The origins of plant pathogens in agro-ecosystems

Plant pathogens can emerge in agricultural ecosystems through several mechanisms, including host-tracking, host jumps, hybridization and horizontal gene transfer. High-throughput DNA sequencing coupled with new analytical approaches make it possible to differentiate among these mechanisms and to infer the time and place where pathogens first emerged. We present several examples to illustrate the different mechanisms and timescales associated with the origins of important plant pathogens. In some cases pathogens were domesticated along with their hosts during the invention of agriculture approximately 10,000 years ago. In other cases pathogens appear to have emerged very recently and almost instantaneously following horizontal gene transfer or hybridization. The predominant unifying feature in these examples is the environmental and genetic uniformity of the agricultural ecosystem in which the pathogens emerged. We conclude that agro-ecosystems will continue to select for new pathogens unless they are re-engineered to make them less conducive to pathogen emergence.

The origin and colonization history of the barley scald pathogen Rhynchosporium secalis

The origins of pathogens and their past and present migration patterns are often unknown. We used phylogenetic haplotype clustering in conjunction with model-based coalescent approaches to reconstruct the genetic history of the barley leaf pathogen Rhynchosporium secalis using the avirulence gene NIP1 and its flanking regions. Our results falsify the hypothesis that R. secalis emerged in association with its host during the domestication of barley 10,000 to 15,000 years ago in the Fertile Crescent and was introduced into Europe through the migration of Neolithic farmers. Estimates of time since most recent common ancestor (2500-5000 BP) placed the emergence of R. secalis clearly after the domestication of barley. We propose that modern populations of R. secalis originated in northern Europe following a host switch, most probably from a wild grass onto cultivated barley shortly after barley was introduced into northern Europe. R. secalis subsequently spread southwards into already established European barley-growing areas.

Genetic diversity of Rhynchosporium secalis in Tunisia as revealed by pathotype, AFLP, and microsatellite analyses

Genetic variability among 122 Rhynchosporium secalis isolates collected from barley in three regions of Tunisia was investigated using host differentials, amplified fragment length polymorphism (AFLP), and microsatellite markers. The isolates were collected from a widely grown scald-susceptible barley cultivar Rihane and a range of local landrace cultivars in geographically distinct regions with different agroclimatic conditions. Pathotypic diversity (the proportion of unique pathotypes) was high in R. secalis populations from the high (100% diversity), moderate (95%), and low (100%) rainfall areas of Tunisia, and from both Rihane (which is the sole variety grown in the high rainfall region) and local landraces (which predominate in the low rainfall area). This may reflect a general adaptability for aggressiveness and suggests that the widely grown cultivar Rihane has exerted little or no selection pressure on the pathogen population since its release in 1983. Genotypic diversity (GD), defined as the probability that two individuals taken at random had different genotypes, was high for populations from Rihane, local landraces, and different agro-ecological zones (GD = 0.96-0.99). There was low genetic differentiation among pathogen populations from different host populations (G(ST) < or = 0.08, theta < or = 0.12) and agro-ecological zones (G(ST) < or = 0.05, theta < or = 0.04), which may be partly explained by gene flow due to the movement of infected stubble around the country. There was no correlation (r = 0.06, P = 0.39) between virulence phenotype and AFLP haplotype. A phenetic tree revealed groups with low bootstrap values that did not reflect the grouping of isolates based on host, pathotype, or agro-ecological region. The implications of these findings for R. secalis evolutionary potential and scald-resistance breeding in Tunisia are discussed.