A whole genome scan of SNP data suggests a lack of abundant hard selective sweeps in the genome of the broad host range plant pathogenic fungus Sclerotinia sclerotiorum

Genetic and Pathogenic Variability among Isolates of Sporisorium reilianum Causing Sorghum Head Smut

Sporisorium reilianum, the causal agent of sorghum (Sorghum bicolor (L.) Moench) head smut, is present in most sorghum-producing regions. This seed replacement fungal disease can reduce yield by up to 80% in severely infected fields. Management of this disease can be challenging due to the appearance of different pathotypes within the pathogenic population. In this research, the genetic variability and pathogenicity of isolates collected from five Texas Counties was conducted. Due to the lack of available space, 21 out of 32 sequenced isolates were selected and evaluated for virulence patterns on the six sorghum differentials, Tx7078, BTx635, SC170-6-17 (TAM2571), SA281 (Early Hegari), Tx414, and BTx643. The results reveal the occurrence of a new pathotype, 1A, and four previously documented US pathotypes when the 21 isolates were evaluated for virulence patterns on the differentials. The most prevalent was pathotype 5, which was recovered from Brazos, Hidalgo, Nueces, and Willacy Counties, Texas. This pathotype was followed by 1A and 6 in frequency of recovery. Pathotype 4 was identified only from isolates collected from Hidalgo County, while pathotype 1 was from Burleson County, Texas. It appeared that the previous US head smut pathotypes (2 and 3) are no longer common, and the new pathotypes, 1A, 5, and 6, are now predominant. The phylogenetic tree constructed from the single-nucleotide polymorphism (SNP) data through the neighbor-joining method showed high genetic diversity among the tested isolates. Some of the diverse clades among the tested isolates were independent of their sampled locations. Notably, HS37, HS49, and HS65 formed a clade and were classified as 1A in the virulence study, while HS 61 and HS 66, which were collected from Nueces County, were grouped and identified as pathotype 5.

The broad host range pathogen Sclerotinia sclerotiorum produces multiple effector proteins that induce host cell death intracellularly

Sclerotinia sclerotiorum is a broad host range necrotrophic fungal pathogen, which causes disease on many economically important crop species. S. sclerotiorum has been shown to secrete small effector proteins to kill host cells and acquire nutrients. We set out to discover novel necrosis-inducing effectors and characterize their activity using transient expression in Nicotiana benthamiana leaves. Five intracellular necrosis-inducing effectors were identified with differing host subcellular localization patterns, which were named intracellular necrosis-inducing effector 1-5 (SsINE1-5). We show for the first time a broad host range pathogen effector, SsINE1, that uses an RxLR-like motif to enter host cells. Furthermore, we provide preliminary evidence that SsINE5 induces necrosis via an NLR protein. All five of the identified effectors are highly conserved in globally sourced S. sclerotiorum isolates. Taken together, these results advance our understanding of the virulence mechanisms employed by S. sclerotiorum and reveal potential avenues for enhancing genetic resistance to this damaging fungal pathogen.

The use of evolutionary analyses to predict functionally relevant traits in filamentous plant pathogens

Identifying traits involved in plant-pathogen interactions is one of the major objectives in molecular plant pathology. Evolutionary analyses may assist in the identification of genes encoding traits that are involved in virulence and local adaptation, including adaptation to agricultural intervention strategies. In the past decades, the number of available genome sequences of fungal plant pathogens has rapidly increased, providing a rich source for the discovery of functionally important genes as well as inference of species histories. Positive selection in the form of diversifying or directional selection leaves particular signatures in genome alignments and can be identified with statistical genetics methods. This review summarises the concepts and approaches used in evolutionary genomics and lists major discoveries related to plant-pathogen adaptative evolution. We underline the significant contribution of evolutionary genomics in discovering virulence-related traits and the study of plant-pathogen ecology and adaptive evolution.

Population Genomics Provide Insights into the Global Genetic Structure of Colletotrichum graminicola, the Causal Agent of Maize Anthracnose

Understanding the genetic diversity and mechanisms underlying genetic variation in pathogen populations is crucial to the development of effective control strategies. We investigated the genetic diversity and reproductive biology of Colletotrichum graminicola isolates which infect maize by sequencing the genomes of 108 isolates collected from 14 countries using restriction site-associated DNA sequencing (RAD-seq) and whole-genome sequencing (WGS). Clustering analyses based on single-nucleotide polymorphisms revealed three genetic groups delimited by continental origin, compatible with short-dispersal of the pathogen and geographic subdivision. Intra- and intercontinental migration was observed between Europe and South America, likely associated with the movement of contaminated germplasm. Low clonality, evidence of genetic recombination, and high phenotypic diversity were detected. We show evidence that, although it is rare (possibly due to losses of sexual reproduction- and meiosis-associated genes) C. graminicola can undergo sexual recombination. Our results support the hypotheses that intra- and intercontinental pathogen migration and genetic recombination have great impacts on the C. graminicola population structure. IMPORTANCE Plant pathogens cause significant reductions in yield and crop quality and cause enormous economic losses worldwide. Reducing these losses provides an obvious strategy to increase food production without further degrading natural ecosystems; however, this requires knowledge of the biology and evolution of the pathogens in agroecosystems. We employed a population genomics approach to investigate the genetic diversity and reproductive biology of the maize anthracnose pathogen (Colletotrichum graminicola) in 14 countries. We found that the populations are correlated with their geographical origin and that migration between countries is ongoing, possibly caused by the movement of infected plant material. This result has direct implications for disease management because migration can cause the movement of more virulent and/or fungicide-resistant genotypes. We conclude that genetic recombination is frequent (in contrast to the traditional view of C. graminicola being mainly asexual), which strongly impacts control measures and breeding programs aimed at controlling this disease.

Phylogenomics and gene selection in Aspergillus welwitschiae: Possible implications in the pathogenicity in Agave sisalana

Aspergillus welwitschiae causes bole rot disease in sisal (Agave sisalana and related species) which affects the production of natural fibers in Brazil, the main worldwide producer of sisal fibers. This fungus is a saprotroph with a broad host range. Previous research established A. welwitschiae as the only causative agent of bole rot in the field, but little is known about the evolution of this species and its strains. In this work, we performed a comparative genomics analysis of 40 Aspergillus strains. We show the conflicting molecular identity of this species, with one sisal-infecting strain sharing its last common ancestor with Aspergillus niger, having diverged only 833 thousand years ago. Furthermore, our analysis of positive selection reveals sites under selection in genes coding for siderophore transporters, Sodium‑calcium exchangers, and Phosphatidylethanolamine-binding proteins (PEBPs). Herein, we discuss the possible impacts of these gene functions on the pathogenicity in sisal.

Transcriptional response to host chemical cues underpins the expansion of host range in a fungal plant pathogen lineage

The host range of parasites is an important factor in assessing the dynamics of disease epidemics. The evolution of pathogens to accommodate new hosts may lead to host range expansion, a process the molecular bases of which are largely enigmatic. The fungus Sclerotinia sclerotiorum has been reported to parasitize more than 400 plant species from diverse eudicot families while its close relative, S. trifoliorum, is restricted to plants from the Fabaceae family. We analyzed S. sclerotiorum global transcriptome reprogramming on hosts from six botanical families and reveal a flexible, host-specific transcriptional program. We generated a chromosome-level genome assembly for S. trifoliorum and found near-complete gene space conservation in two representative strains of broad and narrow host range Sclerotinia species. However, S. trifoliorum showed increased sensitivity to the Brassicaceae defense compound camalexin. Comparative analyses revealed a lack of transcriptional response to camalexin in the S. trifoliorum strain and suggest that regulatory variation in detoxification and effector genes at the population level may associate with the genetic accommodation of Brassicaceae in the Sclerotinia host range. Our work proposes transcriptional plasticity and the co-existence of signatures for generalist and polyspecialist adaptive strategies in the genome of a plant pathogen.

Population Genomics Reveals Molecular Determinants of Specialization to Tomato in the Polyphagous Fungal Pathogen Botrytis cinerea in France

Many fungal plant pathogens encompass multiple populations specialized on different plant species. Understanding the factors underlying pathogen adaptation to their hosts is a major challenge of evolutionary microbiology, and it should help to prevent the emergence of new specialized pathogens on novel hosts. Previous studies have shown that French populations of the gray mold pathogen Botrytis cinerea parasitizing tomato and grapevine are differentiated from each other, and have higher aggressiveness on their host of origin than on other hosts, indicating some degree of host specialization in this polyphagous pathogen. Here, we aimed at identifying the genomic features underlying the specialization of B. cinerea populations to tomato and grapevine. Based on whole genome sequences of 32 isolates, we confirmed the subdivision of B. cinerea pathogens into two genetic clusters on grapevine and another, single cluster on tomato. Levels of genetic variation in the different clusters were similar, suggesting that the tomato-specific cluster has not recently emerged following a bottleneck. Using genome scans for selective sweeps and divergent selection, tests of positive selection based on polymorphism and divergence at synonymous and nonsynonymous sites, and analyses of presence and absence variation, we identified several candidate genes that represent possible determinants of host specialization in the tomato-associated population. This work deepens our understanding of the genomic changes underlying the specialization of fungal pathogen populations.

Complex adaptive architecture underlies adaptation to quantitative host resistance in a fungal plant pathogen

Plant pathogens often adapt to plant genetic resistance so characterization of the architecture underlying such an adaptation is required to understand the adaptive potential of pathogen populations. Erosion of banana quantitative resistance to a major leaf disease caused by polygenic adaptation of the causal agent, the fungus Pseudocercospora fijiensis, was recently identified in the northern Caribbean region. Genome scan and quantitative genetics approaches were combined to investigate the adaptive architecture underlying this adaptation. Thirty-two genomic regions showing host selection footprints were identified by pool sequencing of isolates collected from seven plantation pairs of two cultivars with different levels of quantitative resistance. Individual sequencing and phenotyping of isolates from one pair revealed significant and variable levels of correlation between haplotypes in 17 of these regions with a quantitative trait of pathogenicity (the diseased leaf area). The multilocus pattern of haplotypes detected in the 17 regions was found to be highly variable across all the population pairs studied. These results suggest complex adaptive architecture underlying plant pathogen adaptation to quantitative resistance with a polygenic basis, redundancy, and a low level of parallel evolution between pathogen populations. Candidate genes involved in quantitative pathogenicity and host adaptation of P. fijiensis were identified in genomic regions by combining annotation analysis with available biological data.

Modeling first order additive × additive epistasis improves accuracy of genomic prediction for sclerotinia stem rot resistance in canola

The fungus Sclerotinia sclerotiorum infects hundreds of plant species including many crops. Resistance to this pathogen in canola (Brassica napus L. subsp. napus) is controlled by numerous quantitative trait loci (QTL). For such polygenic traits, genomic prediction may be useful for breeding as it can capture many QTL at once while also considering nonadditive genetic effects. Here, we test application of common regression models to genomic prediction of S. sclerotiorum resistance in canola in a diverse panel of 218 plants genotyped at 24,634 loci. Disease resistance was scored by infection with an aggressive isolate and monitoring over 3 wk. We found that including first-order additive × additive epistasis in linear mixed models (LMMs) improved accuracy of breeding value estimation between 3 and 40%, depending on method of assessment, and correlation between phenotypes and predicted total genetic values by 14%. Bayesian models performed similarly to or worse than genomic relationship matrix-based models for estimating breeding values or overall phenotypes from genetic values. Bayesian ridge regression, which is most similar to the genomic relationship matrix-based approach in the amount of shrinkage it applies to marker effects, was the most accurate of this family of models. This confirms several studies indicating the highly polygenic nature of sclerotinia stem rot resistance. Overall, our results highlight the use of simple epistasis terms for prediction of breeding values and total genetic values for a complex disease resistance phenotype in canola.

Impact of Preconditioning Temperature and Duration Period on Carpogenic Germination of Diverse Sclerotinia sclerotiorum Populations in Southwestern Australia

The soilborne pathogen Sclerotinia sclerotiorum (Lib.) de Bary is the causal agent of Sclerotinia stem rot, a severe disease of broad-leaf crops including canola/rapeseed (Brassica napus) that can result in significant yield losses. Sclerotia, the hard melanized resting structure of the pathogen, requires preconditioning before carpogenic germination can occur. We investigated the effect of preconditioning temperature (4, 20, 35, 50°C, and field conditions) and duration (0, 30, 60, 120, 179, 240, and 301 days) on germination of S. sclerotiorum sclerotia collected from five canola fields in the southwestern Australia grain belt. The ecological diversity of each population was characterized using mycelial compatibility group (MCG) typing. No response was observed for isolates conditioned at 4°C at any time period, indicating that chilling is not a preconditioning requirement for these isolates. Sclerotia required preconditioning for a minimum of 60 days before any significant increase in germination occurred, with no further increases in germination recorded in response to longer conditioning after 60 days. The highest germination was observed in sclerotia conditioned at 50°C. The MCG results indicated significant diversity within and between populations, suggesting local adaptation to different environments as well as ensuring the ability to respond to seasonal variation between years.

Large-Scale Genome Scanning within Exonic Regions Revealed the Contributions of Selective Sweep Prone Genes to Host Divergence and Adaptation in Magnaporthe oryzae Species Complex

Magnaporthe oryzae, one of the most notorious plant pathogens in the agronomic ecosystem, causes a destructive rice blast disease around the world. The blast fungus infects wide arrays of cultivated and non-cultivated plants within the Poaceae. Studies have shown that host speciation exerts selection pressure that drives the evolution and divergence of the M. oryzae population. Population genetic relationship deducted by genome-wide single nucleotide polymorphisms showed that M. oryzae differentiation is highly consistent with the host speciation process. In particular, the rice-infecting population of M. oryzae is distinct from populations from other hosts. However, how genome regions prone to host-mediated selection pressures associated with speciation in M. oryzae, especially at a large-scale population level, has not been extensively characterized. Here, we detected strong evidence of sweep selection throughout the genomes of rice and non-rice pathotypes of M. oryzae population using integrated haplotype score (iHS), cross population extended haplotype homozygosity (XPEHH), and cross population composite likelihood ratio (XPCLR) tests. Functional annotation analyses of the genes associated with host-mediated selection pressure showed that 14 pathogenicity-related genes are under positive selection pressure. Additionally, we showed that 17 candidate effector proteins are under positive and divergent selection among the blast fungus population through sweep selection analysis. Specifically, we find that a divergent selective gene, MGG_13871, is experiencing host-directed mutation in two amino acid residues in rice and non-rice infecting populations. These results provide a crucial insight into the impact of selective sweeping on the differentiation of M. oryzae populations and the dynamic influences of genomic regions in promoting host adaptation and speciation among M. oryzae species.

Carpogenic Germinability of Diverse Sclerotinia sclerotiorum Populations Within the Southwestern Australian Grain Belt

Sclerotinia stem rot, caused by the necrotrophic plant pathogen Sclerotinia sclerotiorum (Lib.) de Bary, is a major disease of canola and pulses in Australia. Current disease management relies greatly on cultural and chemical means of control. Timing of fungicide applications remains a challenge, because efficacy is dependent on accurate prediction of ascospore release and presence on the plant. The aims of this study were to determine the optimal temperature for carpogenic germination of S. sclerotiorum populations sampled from canola and lupin fields in southwestern Australia and characterize diversity using mycelial compatibility groupings (MCGs). Sclerotia were collected from four diseased canola and one diseased lupin field from across southwestern Australia. Forty sclerotia from each population were incubated at four alternating temperatures of 30/15, 20/15, 20/4, and 15/4°C (12-h/12-h light/dark cycle) and assessed every 2 to 3 days for a 180-day period. MCG groupings for populations were characterized using 12 reference isolates. Results indicated the time to initial carpogenic germination decreased as diurnal temperature fluctuations decreased, with a fluctuation of 5°C (20/15°C) having the most rapid initial germination followed by 11°C (15/4°C) followed by 16°C (20/4°C). Optimal germination temperature for all five populations was 20/15°C; however, population responses to other diurnal temperature regimes varied considerably. No germination was observed at 30/15°C. MCG results indicate extensive diversity within and between populations, with at least 40% of sclerotia within each population unable to be characterized. We suggest that this diversity has enabled S. sclerotiorum populations to adapt to varying environmental conditions within southwestern Australia.

Bioinformatic Detection of Positive Selection Pressure in Plant Pathogens: The Neutral Theory of Molecular Sequence Evolution in Action

The genomes of plant pathogenic fungi and oomycetes are often exposed to strong positive selection pressure. During speciation, shifts in host range and preference can lead to major adaptive changes. Furthermore, evolution of total host resistance to most isolates can force rapid evolutionary changes in host-specific pathogens. Crop pathogens are subjected to particularly intense selective pressures from monocultures and fungicides. Detection of the footprints of positive selection in plant pathogen genomes is a worthwhile endeavor as it aids understanding of the fundamental biology of these important organisms. There are two main classes of test for detection of positively selected alleles. Tests based on the ratio of non-synonymous to synonymous substitutions per site detect the footprints of multiple fixation events between divergent lineages. Thus, they are well-suited to the study of ancient adaptation events spanning speciations. On the other hand, tests that scan genomes for local fluctuations in allelic diversity within populations are suitable for detection of recent positive selection in populations. In this review, I briefly describe some of the more widely used tests of positive selection and the theory underlying them. I then discuss various examples of their application to plant pathogen genomes, emphasizing the types of genes that are associated with signatures of positive selection. I conclude with a discussion of the practicality of such tests for identification of pathogen genes of interest and the important features of pathogen ecology that must be taken into account for accurate interpretation.

A detailed in silico analysis of secondary metabolite biosynthesis clusters in the genome of the broad host range plant pathogenic fungus Sclerotinia sclerotiorum

BACKGROUND

The broad host range pathogen Sclerotinia sclerotiorum infects over 400 plant species and causes substantial yield losses in crops worldwide. Secondary metabolites are known to play important roles in the virulence of plant pathogens, but little is known about the secondary metabolite repertoire of S. sclerotiorum. In this study, we predicted secondary metabolite biosynthetic gene clusters in the genome of S. sclerotiorum and analysed their expression during infection of Brassica napus using an existing transcriptome data set. We also investigated their sequence diversity among a panel of 25 previously published S. sclerotiorum isolate genomes.

RESULTS

We identified 80 putative secondary metabolite clusters. Over half of the clusters contained at least three transcriptionally coregulated genes. Comparative genomics revealed clusters homologous to clusters in the closely related plant pathogen Botrytis cinerea for production of carotenoids, hydroxamate siderophores, DHN melanin and botcinic acid. We also identified putative phytotoxin clusters that can potentially produce the polyketide sclerin and an epipolythiodioxopiperazine. Secondary metabolite clusters were enriched in subtelomeric genomic regions, and those containing paralogues showed a particularly strong association with repeats. The positional bias we identified was borne out by intraspecific comparisons that revealed putative secondary metabolite genes suffered more presence / absence polymorphisms and exhibited a significantly higher sequence diversity than other genes.

CONCLUSIONS

These data suggest that S. sclerotiorum produces numerous secondary metabolites during plant infection and that their gene clusters undergo enhanced rates of mutation, duplication and recombination in subtelomeric regions. The microevolutionary regimes leading to S. sclerotiorum secondary metabolite diversity have yet to be elucidated. Several potential phytotoxins documented in this study provide the basis for future functional analyses.

The Notorious Soilborne Pathogenic Fungus Sclerotinia sclerotiorum: An Update on Genes Studied with Mutant Analysis

Ascomycete Sclerotinia sclerotiorum (Lib.) de Bary is one of the most damaging soilborne fungal pathogens affecting hundreds of plant hosts, including many economically important crops. Its genomic sequence has been available for less than a decade, and it was recently updated with higher completion and better gene annotation. Here, we review key molecular findings on the unique biology and pathogenesis process of S. sclerotiorum, focusing on genes that have been studied in depth using mutant analysis. Analyses of these genes have revealed critical players in the basic biological processes of this unique pathogen, including mycelial growth, appressorium establishment, sclerotial formation, apothecial and ascospore development, and virulence. Additionally, the synthesis has uncovered gaps in the current knowledge regarding this fungus. We hope that this review will serve to build a better current understanding of the biology of this under-studied notorious soilborne pathogenic fungus.

Small RNAs from the plant pathogenic fungus Sclerotinia sclerotiorum highlight host candidate genes associated with quantitative disease resistance

Fungal plant pathogens secrete effector proteins and metabolites to cause disease. Additionally, some species transfer small RNAs (sRNAs) into plant cells to silence host mRNAs through complementary base pairing and suppress plant immunity. The fungus Sclerotinia sclerotiorum infects over 600 plant species, but little is known about the molecular processes that govern interactions with its many hosts. In particular, evidence for the production of sRNAs by S. sclerotiorum during infection is lacking. We sequenced sRNAs produced by S. sclerotiorum in vitro and during infection of two host species, Arabidopsis thaliana and Phaseolus vulgaris. We found that S. sclerotiorum produces at least 374 distinct highly abundant sRNAs during infection, mostly originating from repeat-rich plastic genomic regions. We predicted the targets of these sRNAs in A. thaliana and found that these genes were significantly more down-regulated during infection than the rest of the genome. Predicted targets of S. sclerotiorum sRNAs in A. thaliana were enriched for functional domains associated with plant immunity and were more strongly associated with quantitative disease resistance in a genome-wide association study (GWAS) than the rest of the genome. Mutants in A. thaliana predicted sRNA target genes SERK2 and SNAK2 were more susceptible to S. sclerotiorum than wild-type, suggesting that S. sclerotiorum sRNAs may contribute to the silencing of immune components in plants. The prediction of fungal sRNA targets in plant genomes can be combined with other global approaches, such as GWAS, to assist in the identification of plant genes involved in quantitative disease resistance.