A thousand-genome panel retraces the global spread and adaptation of a major fungal crop pathogen

Human activity impacts the evolutionary trajectories of many species worldwide. Global trade of agricultural goods contributes to the dispersal of pathogens reshaping their genetic makeup and providing opportunities for virulence gains. Understanding how pathogens surmount control strategies and cope with new climates is crucial to predicting the future impact of crop pathogens. Here, we address this by assembling a global thousand-genome panel of Zymoseptoria tritici, a major fungal pathogen of wheat reported in all production areas worldwide. We identify the global invasion routes and ongoing genetic exchange of the pathogen among wheat-growing regions. We find that the global expansion was accompanied by increased activity of transposable elements and weakened genomic defenses. Finally, we find significant standing variation for adaptation to new climates encountered during the global spread. Our work shows how large population genomic panels enable deep insights into the evolutionary trajectory of a major crop pathogen.

Complete Genome Sequences of Septoria linicola: A Resource for Studying a Damaging Flax Pathogen

Fungal genus Septoria causes diseases in a wide range of plants. Here, we report the first genome sequences of two strains of Septoria linicola, the causal agent of the pasmo disease of flax (Linum usitatissimum). The genome of the first strain, SE15195, was fully assembled in 16 chromosomes, while 35 unitigs were obtained for a second strain, SE14017. Structural annotations predicted 13,096 and 13,085 protein-encoding genes and transposable elements content of 19.0 and 18.1% of the genome for SE15195 and SE14017, respectively. The four smaller chromosomes 13 to 16 show genomics features of potential accessory chromosomes. The assembly of these two genomes is a new resource for studying S. linicola and improving management of pasmo. [Formula: see text] Copyright © 2022 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Low Amplitude Boom-and-Bust Cycles Define the Septoria Nodorum Blotch Interaction

INTRODUCTION

Septoria nodorum blotch (SNB) is a complex fungal disease of wheat caused by the Dothideomycete fungal pathogen Parastagonospora nodorum. The fungus infects through the use of necrotrophic effectors (NEs) that cause necrosis on hosts carrying matching dominant susceptibility genes. The Western Australia (WA) wheatbelt is a SNB "hot spot" and experiences significant under favorable conditions. Consequently, SNB has been a major target for breeders in WA for many years.

MATERIALS AND METHODS

In this study, we assembled a panel of 155 WA P. nodorum isolates collected over a 44-year period and compared them to 23 isolates from France and the USA using 28 SSR loci.

RESULTS

The WA P. nodorum population was clustered into five groups with contrasting properties. 80% of the studied isolates were assigned to two core groups found throughout the collection location and time. The other three non-core groups that encompassed transient and emergent populations were found in restricted locations and time. Changes in group genotypes occurred during periods that coincided with the mass adoption of a single or a small group of widely planted wheat cultivars. When introduced, these cultivars had high scores for SNB resistance. However, the field resistance of these new cultivars often declined over subsequent seasons prompting their replacement with new, more resistant varieties. Pathogenicity assays showed that newly emerged isolates non-core are more pathogenic than old isolates. It is likely that the non-core groups were repeatedly selected for increased virulence on the contemporary popular cultivars.

DISCUSSION

The low level of genetic diversity within the non-core groups, difference in virulence, low abundance, and restriction to limited locations suggest that these populations more vulnerable to a population crash when the cultivar was replaced by one that was genetically different and more resistant. We characterize the observed pattern as a low-amplitude boom-and-bust cycle in contrast with the classical high amplitude boom-and-bust cycles seen for biotrophic pathogens where the contrast between resistance and susceptibility is typically much greater. Implications of the results are discussed relating to breeding strategies for more sustainable SNB resistance and more generally for pathogens with NEs.

Genome structure and reproductive behaviour influence the evolutionary potential of a fungal phytopathogen

Modern agriculture favours the selection and spread of novel plant diseases. Furthermore, crop genetic resistance against pathogens is often rendered ineffective within a few years of its commercial deployment. Leptosphaeria maculans, the cause of phoma stem canker of oilseed rape, develops gene-for-gene interactions with its host plant, and has a high evolutionary potential to render ineffective novel sources of resistance in crops. Here, we established a four-year field experiment to monitor the evolution of populations confronted with the newly released Rlm7 resistance and to investigate the nature of the mutations responsible for virulence against Rlm7. A total of 2551 fungal isolates were collected from experimental crops of a Rlm7 cultivar or a cultivar without Rlm7. All isolates were phenotyped for virulence and a subset was genotyped with neutral genetic markers. Virulent isolates were investigated for molecular events at the AvrLm4-7 locus. Whilst virulent isolates were not found in neighbouring crops, their frequency had reached 36% in the experimental field after four years. An extreme diversity of independent molecular events leading to virulence was identified in populations, with large-scale Repeat Induced Point mutations or complete deletion of AvrLm4-7 being the most frequent. Our data suggest that increased mutability of fungal genes involved in the interactions with plants is directly related to their genomic environment and reproductive system. Thus, rapid allelic diversification of avirulence genes can be generated in L. maculans populations in a single field provided that large population sizes and sexual reproduction are favoured by agricultural practices.

Genes under positive selection in a model plant pathogenic fungus, Botrytis

The rapid evolution of particular genes is essential for the adaptation of pathogens to new hosts and new environments. Powerful methods have been developed for detecting targets of selection in the genome. Here we used divergence data to compare genes among four closely related fungal pathogens adapted to different hosts to elucidate the functions putatively involved in adaptive processes. For this goal, ESTs were sequenced in the specialist fungal pathogens Botrytis tulipae and Botrytis ficariarum, and compared with genome sequences of Botrytis cinerea and Sclerotinia sclerotiorum, responsible for diseases on over 200 plant species. A maximum likelihood-based analysis of 642 predicted orthologs detected 21 genes showing footprints of positive selection. These results were validated by resequencing nine of these genes in additional Botrytis species, showing they have also been rapidly evolving in other related species. Twenty of the 21 genes had not previously been identified as pathogenicity factors in B. cinerea, but some had functions related to plant-fungus interactions. The putative functions were involved in respiratory and energy metabolism, protein and RNA metabolism, signal transduction or virulence, similarly to what was detected in previous studies using the same approach in other pathogens. Mutants of B. cinerea were generated for four of these genes as a first attempt to elucidate their functions.

Genomic analysis of the necrotrophic fungal pathogens Sclerotinia sclerotiorum and Botrytis cinerea

Sclerotinia sclerotiorum and Botrytis cinerea are closely related necrotrophic plant pathogenic fungi notable for their wide host ranges and environmental persistence. These attributes have made these species models for understanding the complexity of necrotrophic, broad host-range pathogenicity. Despite their similarities, the two species differ in mating behaviour and the ability to produce asexual spores. We have sequenced the genomes of one strain of S. sclerotiorum and two strains of B. cinerea. The comparative analysis of these genomes relative to one another and to other sequenced fungal genomes is provided here. Their 38-39 Mb genomes include 11,860-14,270 predicted genes, which share 83% amino acid identity on average between the two species. We have mapped the S. sclerotiorum assembly to 16 chromosomes and found large-scale co-linearity with the B. cinerea genomes. Seven percent of the S. sclerotiorum genome comprises transposable elements compared to <1% of B. cinerea. The arsenal of genes associated with necrotrophic processes is similar between the species, including genes involved in plant cell wall degradation and oxalic acid production. Analysis of secondary metabolism gene clusters revealed an expansion in number and diversity of B. cinerea-specific secondary metabolites relative to S. sclerotiorum. The potential diversity in secondary metabolism might be involved in adaptation to specific ecological niches. Comparative genome analysis revealed the basis of differing sexual mating compatibility systems between S. sclerotiorum and B. cinerea. The organization of the mating-type loci differs, and their structures provide evidence for the evolution of heterothallism from homothallism. These data shed light on the evolutionary and mechanistic bases of the genetically complex traits of necrotrophic pathogenicity and sexual mating. This resource should facilitate the functional studies designed to better understand what makes these fungi such successful and persistent pathogens of agronomic crops.

Effector diversification within compartments of the Leptosphaeria maculans genome affected by Repeat-Induced Point mutations

Fungi are of primary ecological, biotechnological and economic importance. Many fundamental biological processes that are shared by animals and fungi are studied in fungi due to their experimental tractability. Many fungi are pathogens or mutualists and are model systems to analyse effector genes and their mechanisms of diversification. In this study, we report the genome sequence of the phytopathogenic ascomycete Leptosphaeria maculans and characterize its repertoire of protein effectors. The L. maculans genome has an unusual bipartite structure with alternating distinct guanine and cytosine-equilibrated and adenine and thymine (AT)-rich blocks of homogenous nucleotide composition. The AT-rich blocks comprise one-third of the genome and contain effector genes and families of transposable elements, both of which are affected by repeat-induced point mutation, a fungal-specific genome defence mechanism. This genomic environment for effectors promotes rapid sequence diversification and underpins the evolutionary potential of the fungus to adapt rapidly to novel host-derived constraints.

Occurrence of a new subclade of Leptosphaeria biglobosa in Western Australia

Stem canker of crucifers is caused by an ascomycete species complex comprising of two main species, Leptosphaeria maculans and L. biglobosa. These are composed of at least seven distinct subclades based on biochemical data or on sequences of internal transcribed spacer (ITS), the mating type MAT1-2 or fragments of actin or beta-tubulin genes. In the course of a wide-scale characterization of the race structure of L. maculans from Western Australia, a few isolates from two locations failed to amplify specific sequences of L. maculans, i.e., the mating-type or minisatellite alleles. Based on both pathogenicity tests and ITS size, these isolates were classified as belonging to the L. biglobosa species. Parsimony and distance analyses performed on ITS, actin and beta-tubulin sequences revealed that these isolates formed a new L. biglobosa subclade, more related to the Canadian L. biglobosa 'canadensis' subclade than to the L. biglobosa 'australensis' isolates previously described in Australia (Victoria). They are termed here as L. biglobosa 'occiaustralensis'. These isolates were mainly recovered from resistant oilseed rape cultivars that included the Brassica rapa sp. sylvestris-derived resistance source, but not from the susceptible cv. Westar. The pathogenicity of L. biglobosa 'occiaustralensis' to cotyledons of most oilseed rape genotypes was higher than that of L. biglobosa 'canadensis' or L. biglobosa 'australensis' isolates.

Heterochromatin-like regions as ecological niches for avirulence genes in the Leptosphaeria maculans genome: map-based cloning of AvrLm6

Map-based cloning of avirulence genes of the AvrLml-2-6 cluster was recently undertaken in Leptosphaeria maculans and led to the identification of AvrLm1. The ensuing chromosome walk toward AvrLm6 resulted in the delineation of a 562-kb bacterial artificial chromosome (BAC) clone contig in an avirulent isolate. Following sequencing of the contig and sequence comparison with a virulent isolate, four AvrLm6 candidate genes were identified. Complementation of the virulent isolate with the four candidates was performed and one gene was found to fully restore the avirulent phenotype on Rlm6 oilseed rape genotypes. AvrLm6 was found to be located in the same genome context as AvrLml, because it is a solo gene surrounded by 85 and 48 kb of degenerated repeats on its 5' and 3' sides, respectively. AvrLm6 is an orphan gene encoding a small, potentially secreted, cysteine-rich protein. Comparison of AvrLm1 and AvrLm6 expressions by quantitative reverse-transcription polymerase chain reaction revealed that both genes are highly overexpressed during primary leaf infection. Using RNA interference, decreasing expression of AvrLm6 was shown to result in virulence toward Rlm6 genotypes whenever the expression was reduced by more than 60% compared with the wild-type isolate.

Lost in the middle of nowhere: theAvrLm1avirulence gene of the DothideomyceteLeptosphaeria maculans

Leptosphaeria maculans, a Dothideomycete causing stem canker on oilseed rape (Brassica napus), develops gene-for-gene interactions with its host plants. To date, nine resistance genes (Rlm1-9) have been identified in Brassica spp. The corresponding nine avirulence genes (AvrLm1-9) in L. maculans have been mapped at four independent loci, thereby revealing two clusters of three and four linked avirulence genes. Here, we report the completion of map-based cloning of AvrLm1. AvrLm1 was genetically delineated within a 7.3 centimorgan interval corresponding to a 439 kb BAC contig. AvrLm1 is a single copy gene isolated within a 269 kb non-coding, heterochromatin-like region. The region comprised a number of degenerated, nested copies of four long-terminal repeat (LTR) retrotransposons, including Pholy and three novel Gypsy-like retrotransposons. AvrLm1 restored the avirulent phenotype on Rlm1 cultivars following functional complementation of virulent isolates. AvrLm1 homologues were not detected in other Leptosphaeria species or in known fungal genomes including the closely related species Stagonospora nodorum. The predicted AvrLm1 protein is composed of 205 amino acids, of which only one is a cysteine residue. It contains a peptide signal suggesting extracellular localization. Unlike most other fungal avirulence genes, AvrLm1 is constitutively expressed, with a probable increased level of expression upon plant infection, suggesting the absence of tight regulation of AvrLm1 expression.

Genetic variability and distribution of mating type alleles in field populations of Leptosphaeria maculans from France

Leptosphaeria maculans is the most ubiquitous fungal pathogen of Brassica crops and causes the devastating stem canker disease of oilseed rape worldwide. We used minisatellite markers to determine the genetic structure of L. maculans in four field populations from France. Isolates were collected at three different spatial scales (leaf, 2-m2 field plot, and field) enabling the evaluation of spatial distribution of the mating type alleles and of genetic variability within and among field populations. Within each field population, no gametic disequilibrium between the minisatellite loci was detected and the mating type alleles were present at equal frequencies. Both sexual and asexual reproduction occur in the field, but the genetic structure of these populations is consistent with annual cycles of randomly mating sexual reproduction. All L. maculans field populations had a high level of gene diversity (H = 0.68 to 0.75) and genotypic diversity. Within each field population, the number of genotypes often was very close to the number of isolates. Analysis of molecular variance indicated that >99.5% of the total genetic variability was distributed at a small spatial scale, i.e., within 2-m2 field plots. Population differentiation among the four field populations was low (GST < 0.02), suggesting a high degree of gene exchange between these populations. The high gene flow evidenced here in French populations of L. maculans suggests a rapid countrywide diffusion of novel virulence alleles whenever novel resistance sources are used.

Analysis of Leptosphaeria maculans Race Structure in a Worldwide Collection of Isolates

ABSTRACT Leptosphaeria maculans, the causal agent of stem canker of oilseed rape, develops gene-for-gene interactions with its hosts. To date, eight L. maculans avirulence (Avr) genes, AvrLm1 to AvrLm8, have been genetically characterized. An additional Avr gene, AvrLm9, that interacts with the resistance gene Rlm9, was genetically characterized here following in vitro crosses of the pathogen. A worldwide collection of 63 isolates, including the International Blackleg of Crucifers Network collection, was genotyped at these nine Avr loci. In a first step, isolates were classified into pathogenicity groups (PGs) using two published differential sets. This analysis revealed geographical disparities as regards the proportion of each PG. Genotyping of isolates at all Avr loci confirmed the disparities between continents, in terms of Avr allele frequencies, particularly for AvrLm2, AvrLm3, AvrLm7, AvrLm8, and AvrLm9, or in terms of race structure, diversity, and complexity. Twenty-six distinct races were identified in the collection. A larger number of races (n = 18) was found in Australia than in Europe (n = 8). Mean number of virulence alleles per isolate was also higher in Australia (5.11 virulence alleles) than in Europe (4.33) and Canada (3.46). Due to the diversity of populations of L. maculans evidenced here at the race level, a new, open terminology is proposed for L. maculans race designation, indicating all Avr loci for which the isolate is avirulent.

Truncated and RIP-degenerated copies of the LTR retrotransposon are clustered in a pericentromeric region of the genome

The LMR1 5.2 kb interspersed repeat of Leptosphaeria maculans was described by Taylor and Borgmann [Mol. Plant Microbe Interact. 7 (1994) 181] as an uncharacterized repeated element sharing homologies with both LINEs and SINEs. Here, we used the LMR1 sequence as a template to identify the full-length element within a 184-kb genomic sequence corresponding to the pericentromeric region of the 2.80 Mb chromosome of isolate v23.1.3. This region comprises (i) one 6980-bp full-sized Pholy element bordered by two 275- to 280-bp long terminal repeats (LTRs), (ii) five Pholy-related sequences, usually truncated at their 3' ends, and (iii) five solo-LTRs. Structural features strongly suggested that Pholy corresponds to an ancient copia-like retrotransposon, sharing strong homologies with the Elsa retrotransposon of Stagonospora nodorum. Pholy was also suggested to be specific to pericentromeric regions. Comparative analysis of the structure of the Pholy-like sequences occurring in the 184-kb contig and in other parts of the genome showed that this family of repeats is highly degenerated following extensive repeat induced point mutation (RIP).

Identification and characterization of polymorphic minisatellites in the phytopathogenic ascomycete Leptosphaeria maculans

Leptosphaeria maculans causes phoma stem canker, the most serious disease of oilseed rape world-wide. Sexual recombination is important in the pathogen life cycle and increases the risk of plant resistance genes being overcome rapidly. Thus, there is a need to develop easy-to-use molecular markers suitable for large-scale population genetic studies. The minisatellite MinLm1, showing six alleles in natural populations, has previously been used as a marker to survey populations. Here, we report the characterization of five new minisatellites (MinLm2-MinLm6), of which four were identified by a systematic search for tandemly repeated polymorphic regions in BAC-end sequencing data from L. maculans. Of 782 BAC-end sequences analysed, 43 possessed putative minisatellite-type repeats and four of these (MinLm3-MinLm6) displayed both consistent PCR amplification and size polymorphism in a collection of L. maculans isolates of diverse origins. Cloning and sequencing of each allele confirmed that polymorphism was due to variation in the repeat number of a core motif ranging from 11 bp (MinLm3) to 51 bp (MinLm4). The number of alleles found for each minisatellite ranged from three (MinLm4) to nine (MinLm2), with eight, five and six for MinLm3, MinLm5 and MinLm6, respectively. MinLm2-MinLm6 are all single locus markers specific to L. maculans and share some common features, such as conservation of core motifs and incomplete direct repeats in the flanking regions. To our knowledge, L. maculans is the first fungal species for which six polymorphic single locus minisatellite markers have been reported.

Analysis of molecular markers genetically linked to the Leptosphaeria maculans avirulence gene AvrLm1 in field populations indicates a highly conserved event leading to virulence on Rlm1 genotypes

Map-based cloning of the avirulence gene AvrLm1 of Leptosphaeria maculans was initiated utilizing a genetic map of the fungus and a BAC library constructed from an AvrLm1 isolate. Seven polymorphic DNA markers closely linked to AvrLm1 were identified. Of these, two were shown to border the locus on its 5' end and were present, with size polymorphism, in both the virulent and the avirulent isolates. In contrast, three markers, J19-1.1, J53-1.3 (in coupling phase with avirulence), and Vir1 (in repulsion phase with avirulence), cosegregated with AvrLm1 in 312 progeny from five in vitro crosses. J19-1.1 and J53-1.3 were never amplified in the virulent parents or progeny, whereas Vir1 was never amplified in the avirulent parents or progeny. J19-1.1 and J53-1.3 were shown to be separated by 40 kb within a 184-kb BAC contig. In addition, the 1.6-cM genetic distance between J53-1.3 and the nearest recombinant marker corresponded to a 121-kb physical distance. When analyzing a European Union-wide collection of 192 isolates, J53-1.3, J19-1.1, and Vir1 were found to be closely associated with the AvrLm1 locus. The results of polymerase chain reaction amplification with primers for the three markers were in accordance with the interaction phenotype for 92.2% (J53-1.3), 90.6% (J19-1.1), and 88.0% (Vir1) of the isolates. In addition, genome organization of the AvrLm1 region was highly conserved in field isolates, because 89.1% of the avirulent isolates and 79.0% of the virulent isolates showed the same association of markers as that of the parents of in vitro crosses. The large-scale analysis of field isolates with markers originating from the genetic map therefore confirms (i) the physical proximity between the markers and the target locus and (ii) that AvrLm1 is located in (or close to) a recombination-deficient genome region. As a consequence, map-based markers provided us with high-quality markers for an overview of the occurrence of race "AvrLm1" at the field scale. These data were used to propose hypotheses on evolution towards virulence in field isolates.

Molecular characterisation and polymorphism of MinLm1, a minisatellite from the phytopathogenic ascomycete Leptosphaeria maculans

A sequence-characterised amplified region marker was identified in the phytopathogenic fungus Leptosphaeria maculans, which generated a single-banding pattern corresponding to six alleles showing size polymorphism between L. maculans field isolates. The size polymorphism was due to 2-7 tandem repeats of the 23-bp motif 5' TCTTACTTACATACACACCTCCC 3'. The repeated sequence, termed MinLm1, shares many features specific to minisatellites, e.g. a very strong G/C strand asymmetry, the presence of 6-bp direct repeats at both ends of the sequence and its occurrence in a region rich in microsatellites such as (CT)n, (ATG)n, (GTG)n and (CAT)n. MinLm1 shows a very high degree of conservation of the bases from one repeat to another and from one isolate to another (percent match range: 99.6-100%), whatever their geographical or temporal relatedness. MinLm1 is a single-locus minisatellite located on chromosomes sized 2.79 Mb and 2.48 Mb, of L. maculans isolates a.2 and H5, respectively. In agricultural populations of L. maculans, two alleles of MinLm1 were prevalent, corresponding to 2x and 5x repeats of the core motif. Differences in allele frequencies were observed in some cropping conditions, suggesting that MinLm1 is an informative marker for epidemiological studies of the pathogen.