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.     

Evolution of linked avirulence effectors in Leptosphaeria maculans is affected by genomic environment and exposure to resistance genes in host plants

Brassica napus (canola) cultivars and isolates of the blackleg fungus, Leptosphaeria maculans interact in a 'gene for gene' manner whereby plant resistance (R) genes are complementary to pathogen avirulence (Avr) genes. Avirulence genes encode proteins that belong to a class of pathogen molecules known as effectors, which includes small secreted proteins that play a role in disease. In Australia in 2003 canola cultivars with the Rlm1 resistance gene suffered a breakdown of disease resistance, resulting in severe yield losses. This was associated with a large increase in the frequency of virulence alleles of the complementary avirulence gene, AvrLm1, in fungal populations. Surprisingly, the frequency of virulence alleles of AvrLm6 (complementary to Rlm6) also increased dramatically, even though the cultivars did not contain Rlm6. In the L. maculans genome, AvrLm1 and AvrLm6 are linked along with five other genes in a region interspersed with transposable elements that have been degenerated by Repeat-Induced Point (RIP) mutations. Analyses of 295 Australian isolates showed deletions, RIP mutations and/or non-RIP derived amino acid substitutions in the predicted proteins encoded by these seven genes. The degree of RIP mutations within single copy sequences in this region was proportional to their proximity to the degenerated transposable elements. The RIP alleles were monophyletic and were present only in isolates collected after resistance conferred by Rlm1 broke down, whereas deletion alleles belonged to several polyphyletic lineages and were present before and after the resistance breakdown. Thus, genomic environment and exposure to resistance genes in B. napus has affected the evolution of these linked avirulence genes in L. maculans.     

Crystal structure of the effector AvrLm4–7 of Leptosphaeria maculans reveals insights into its translocation into plant cells and recognition by resistance proteins

The avirulence gene AvrLm4–7 of Leptosphaeria maculans, the causal agent of stem canker in Brassica napus (oilseed rape), confers a dual specificity of recognition by two resistance genes (Rlm4 and Rlm7) and is strongly involved in fungal fitness. In order to elucidate the biological function of AvrLm4–7 and understand the specificity of recognition by Rlm4 and Rlm7, the AvrLm4–7 protein was produced in Pichia pastoris and its crystal structure was determined. It revealed the presence of four disulfide bridges, but no close structural analogs could be identified. A short stretch of amino acids in the C terminus of the protein, (R/N)(Y/F)(R/S)E(F/W), was well‐conserved among AvrLm4–7 homologs. Loss of recognition of AvrLm4–7 by Rlm4 is caused by the mutation of a single glycine to an arginine residue located in a loop of the protein. Loss of recognition by Rlm7 is governed by more complex mutational patterns, including gene loss or drastic modifications of the protein structure. Three point mutations altered residues in the well‐conserved C–terminal motif or close to the glycine involved in Rlm4‐mediated recognition, resulting in the loss of Rlm7‐mediated recognition. Transient expression in Nicotiana benthamiana (tobacco) and particle bombardment experiments on leaves from oilseed rape suggested that AvrLm4–7 interacts with its cognate R proteins inside the plant cell, and can be translocated into plant cells in the absence of the pathogen. Translocation of AvrLm4–7 into oilseed rape leaves is likely to require the (R/N)(Y/F)(R/S)E(F/W) motif as well as an RAWG motif located in a nearby loop that together form a positively charged region.     

One gene‐one name: the AvrLmJ1 avirulence gene of Leptosphaeria maculans is AvrLm5

Leptosphaeria maculans, the causal agent of blackleg disease, interacts with Brassica napus (oilseed rape, canola) and other Brassica hosts in a gene‐for‐gene manner. The avirulence gene AvrLmJ1 has been cloned previously and shown to interact with an unidentified Brassica juncea resistance gene. In this study, we show that the AvrLmJ1 gene maps to the same position as the AvrLm5 locus. Furthermore, isolates complemented with the AvrLmJ1 locus confer avirulence towards B. juncea genotypes harbouring Rlm5. These findings demonstrate that AvrLmJ1 is AvrLm5 and highlight the need for shared resources to characterize accurately avirulence and/or resistance genes.     

Rapid identification of the Leptosphaeria maculans avirulence gene AvrLm2 using an intraspecific comparative genomics approach

Five avirulence genes from Leptosphaeria maculans, the causal agent of blackleg of canola (Brassica napus), have been identified previously through map‐based cloning. In this study, a comparative genomic approach was used to clone the previously mapped AvrLm2. Given the lack of a presence–absence gene polymorphism coincident with the AvrLm2 phenotype, 36 L. maculans isolates were resequenced and analysed for single‐nucleotide polymorphisms (SNPs) in predicted small secreted protein‐encoding genes present within the map interval. Three SNPs coincident with the AvrLm2 phenotype were identified within LmCys1, previously identified as a putative effector‐coding gene. Complementation of a virulent isolate with LmCys1, as the candidate AvrLm2 allele, restored the avirulent phenotype on Rlm2‐containing B. napus lines. AvrLm2 encodes a small cysteine‐rich protein with low similarity to other proteins in the public databases. Unlike other avirulence genes, AvrLm2 resides in a small GC island within an AT‐rich isochore of the genome, and was never found to be deleted completely in virulent isolates.     

A new avirulence gene of Leptosphaeria maculans,AvrLm14, identifies a resistance source in American broccoli (Brassica oleracea) genotypes

In many cultivated crops, sources of resistance to diseases are sparse and rely on introgression from wild relatives. Agricultural crops often are allopolyploids resulting from interspecific crosses between related species, which are sources of diversity for resistance genes. This is the case for Brassica napus (oilseed rape, canola), an interspecific hybrid between Brassica rapa (turnip) and Brassica oleracea (cabbage). B. napus has a narrow genetic basis and few effective resistance genes against stem canker (blackleg) disease, caused by the fungus Leptosphaeria maculans, are currently available. B. rapa diversity has proven to be a valuable source of resistance (Rlm, LepR) genes, while B. oleracea genotypes were mostly considered susceptible. Here we identified a new resistance source in B. oleracea genotypes from America, potentially effective against French L. maculans isolates under both controlled and field conditions. Genetic analysis of fungal avirulence and subsequent cloning and validation identified a new avirulence gene termed AvrLm14 and suggested a typical gene-for-gene interaction between AvrLm14 and the postulated Rlm14 gene. AvrLm14 shares all the usual characteristics of L. maculans avirulence genes: it is hosted in a genomic region enriched in transposable elements and heterochromatin marks H3K9me3, its expression is repressed during vegetative growth but shows a strong overexpression 5–9 days following cotyledon infection, and it encodes a small secreted protein enriched in cysteine residues with few matches in databases. Similar to the previously cloned AvrLm10-A, AvrLm14 contributes to reduce lesion size on susceptible cotyledons, pointing to a complex interplay between effectors promoting or reducing lesion development.     

Use of high resolution melting analysis to genotype the avirulence AvrLm4‐7 gene of Leptosphaeria maculans,a fungal pathogen of Brassica napus

The oilseed rape (Brassica napus) stem canker disease, due to the fungal pathogen Leptosphaeria maculans, is mainly controlled by host genetic resistance. Since 2004, the specific resistance gene Rlm7 is widely used in France. Specific resistance is effective when fungal populations are mainly composed of avirulent isolates. The development of molecular tools for the identification of virulent isolates towards Rlm7 was needed to undertake large‐scale surveys and to monitor the emergence of virulent populations in fields. Previous studies have described a large diversity of molecular events leading to virulence towards Rlm7, rendering conventional polymerase chain reaction (PCR) methods inapplicable to identify virulent isolates. Interestingly, a very limited nucleotide polymorphism was observed for avirulent, AvrLm7, alleles. Such characteristics were exploited here to develop a diagnostic method based on high resolution melting (HRM) analysis of the AvrLm4‐7 gene. High resolution melting analysis of a collection of 206 reference isolates revealed only four different profiles within 100 avirulent isolates and 87% of virulent isolates showed either no amplification or HRM curves distinct from those of avirulent isolates. The reliability of the method was confirmed using a second set of 119 unknown isolates, for which biological phenotyping and HRM genotyping were in agreement for 93% of the isolates. HRM combined with the PCR amplification of a larger fragment encompassing AvrLm4‐7 led to a correct diagnostic for 97.5% of the isolates.     

Leptosphaeria maculans AvrLm9: a new player in the game of hide and seek with AvrLm4‐7

Blackleg disease of Brassica napus caused by Leptosphaeria maculans (Lm) is largely controlled by the deployment of race‐specific resistance (R) genes. However, selection pressure exerted by R genes causes Lm to adapt and give rise to new virulent strains through mutation and deletion of effector genes. Therefore, a knowledge of effector gene function is necessary for the effective management of the disease. Here, we report the cloning of Lm effector AvrLm9 which is recognized by the resistance gene Rlm9 in B. napus cultivar Goéland. AvrLm9 was mapped to scaffold 7 of the Lm genome, co‐segregating with the previously reported AvrLm5 (previously known as AvrLmJ1). Comparison of AvrLm5 alleles amongst the 37 re‐sequenced Lm isolates and transgenic complementation identified a single point mutation correlating with the AvrLm9 phenotype. Therefore, we renamed this gene as AvrLm5‐9 to reflect the dual specificity of this locus. Avrlm5‐9 transgenic isolates were avirulent when inoculated on the B. napus cultivar Goéland. The expression of AvrLm5‐9 during infection was monitored by RNA sequencing. The recognition of AvrLm5‐9 by Rlm9 is masked in the presence of AvrLm4‐7, another Lm effector. AvrLm5‐9 and AvrLm4‐7 do not interact, and AvrLm5‐9 is expressed in the presence of AvrLm4‐7. AvrLm5‐9 is the second Lm effector for which host recognition is masked by AvrLm4‐7. An understanding of this complex interaction will provide new opportunities for the engineering of broad‐spectrum recognition.     

Leptosphaeria maculans avirulence gene AvrLm4-7 confers a dual recognition specificity by the Rlm4 and Rlm7 resistance genes of oilseed rape, and circumvents Rlm4-mediated recognition through a single amino acid change

Leptosphaeria maculans is the ascomycete responsible for one of the most damaging diseases of oilseed rape ( Brassica napus ), stem canker of crucifers. Both avirulence ( AvrLm ) genes in the fungus and resistance ( Rlm ) genes in the plant are genetically clustered. Using a map-based cloning strategy, we delineated a 238 kb region containing the AvrLm7 locus. Structural features of the region were reminiscent of those previously found on another chromosome for genomic regions encompassing AvrLm1 and AvrLm6 , i.e. GC-equilibrated, gene-rich isochores alternating with AT-rich, recombination-deficient, gene-poor isochores. These latter corresponded to mosaics of degenerated and truncated transposable elements. AvrLm7 is the only gene located within a 60 kb AT-rich isochore. It induced resistance responses in plants harbouring either Rlm7 or Rlm4 , and was thus renamed AvrLm4-7 . It encodes a 143-amino-acid cysteine-rich protein, predicted to be secreted, and strongly induced during early stages of plant infection. Sequencing and restriction analyses of AvrLm4–AvrLm7 or avrLm4–AvrLm7 alleles in L. maculans field isolates, and targeted point mutagenesis strongly suggested that one single base mutation, leading to the change of a glycine to an arginine residue, was responsible for the loss of AvrLm4 specificity whereas AvrLm7 recognition was unaltered.     

Molecular screening for avirulence alleles AvrLm1 and AvrLm6 in airborne inoculum of Leptosphaeria maculans and winter oilseed rape (Brassica napus) plants from Poland and the UK

A combination of staining, light microscopy and SYBR green- and dual-labelled fluorescent probe-based qPCR chemistries with species- and gene-specific primers was employed to evaluate fluctuations in the aerial biomass of Leptosphaeria maculans spores captured by volumetric spore trappings in Poznan, Poland (2006, 2008) and Harpenden, UK (2002, 2006). Arising from these surveys, DNA samples extracted from Burkard spore-trap tapes were screened for fluctuation patterns in the frequencies of AvrLm1 and AvrLm6, the most prominent of the 15 genes that code for avirulence effectors in this Dothideomycete cause of the destructive phoma stem canker disease of oilseed rape worldwide. In Poznan, very low frequencies of AvrLm1 allele were found in the autumn of both 2006 and 2008, reflecting significantly increased cultivation of rape seed with Rlm1-based resistance. In contrast, at least six folds-higher frequencies of AvrLm6, which were also confirmed by end-point PCR bioassays on phoma-infected leaves from the same region of Poland, were obtained during both years. In the UK, however, relatively higher AvrLm1 allele titres were found in L. maculans spores captured in air samples from the autumn of 2002 on the experimental fields of Rothamsted Research, Harpenden, that were historically sown to genetically heterogeneous B. napus cultivars. In the 2006 screen these levels had plummeted, to a 1:4 ratio, in favour of frequencies of the AvrLm6 allele. Patterns of fluctuations in erg11 (CYP51) fragments coding for sterol 14α-demethylase suggest October as the month with the most viable wind-dispersed L. maculans propagules of each season of the screens.     

Genetic mapping of the <Emphasis Type="Italic">Leptosphaeria maculans</Emphasis> avirulence gene corresponding to the <Emphasis Type="Italic">LepR1</Emphasis> resistance gene of <Emphasis Type="Italic">Brassica napus</Emphasis>

AvrLepR1 of the fungal pathogen Leptosphaeria maculans is the avirulence gene that corresponds to Brassica LepR1, a plant gene controlling dominant, race-specific resistance to this pathogen. An in vitro cross between the virulent L. maculans isolate, 87-41, and the avirulent isolate, 99-56, was performed in order to map the AvrLepR1 gene. The disease reactions of the 94 of the resulting F₁ progenies were tested on the canola line ddm-12-6s-1, which carries LepR1. There were 44 avirulent progenies and 50 virulent progenies suggesting a 1:1 segregation ratio and that the avirulence of 99-56 on ddm-12-6s-1 is controlled by a single gene. Tetrad analysis also indicated a 1:1 segregation ratio. The AvrLepR1 gene was positioned on a genetic map of L. maculans relative to 259 sequence-related amplified polymorphism (SRAP) markers, two cloned avirulence genes (AvrLm1 and AvrLm4-7) and the mating type locus (MAT1). The genetic map consisted of 36 linkage groups, ranging in size from 13.1 to 163.7 cM, and spanned a total of 2,076.4 cM. The AvrLepR1 locus was mapped to linkage group 4, in the 13.1 cM interval flanked by the SRAP markers SBG49-110 and FT161-223. The AvrLm4-7 locus was also positioned on linkage group 4, close to but distinct from the AvrLepR1 locus, in the 5.4 cM interval flanked by FT161-223 and P1314-300. This work will make possible the further characterization and map-based cloning of AvrLepR1. A combination of genetic mapping and pathogenicity tests demonstrated that AvrLepR1 is different from each of the L. maculans avirulence genes that have been characterized previously.     

Lost in the middle of nowhere: the AvrLm1 avirulence gene of the Dothideomycete Leptosphaeria 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.     

Detection,introgression and localization of genes conferring specific resistance to <Emphasis Type="Italic">Leptosphaeria maculans</Emphasis> from <Emphasis Type="Italic">Brassica rapa</Emphasis> into <Emphasis Type="Italic">B. napus</Emphasis>

Blackleg (stem canker) caused by the fungus Leptosphaeria maculans is one of the most damaging diseases of oilseed rape (Brassica napus). Crop relatives represent a valuable source of “new” resistance genes that could be used to diversify cultivar resistance. B. rapa, one of the progenitors of B. napus, is a potential source of new resistance genes. However, most of the accessions are heterozygous so it is impossible to directly detect the plant genes conferring specific resistance due to the complex patterns of avirulence genes in L. maculans isolates. We developed a strategy to simultaneously characterize and introgress resistance genes from B. rapa, by homologous recombination, into B. napus. One B. rapa plant resistant to one L. maculans isolate was used to produce B. rapa backcross progeny and a resynthesized B. napus plant from which a population of doubled haploid lines was derived after crossing with natural B. napus. We then used molecular analyses and resistance tests on these populations to identify and map the resistance genes and to characterize their introgression from B. rapa into B. napus. Three specific genes conferring resistance to L. maculans (Rlm1, Rlm2 and Rlm7) were identified in B. rapa. Comparisons of genetic maps showed that two of these genes were located on the R7 linkage group, in a region homologous to the region on linkage group N7 in B. napus, where these genes have been reported previously. The results of our study offer new perspectives for gene introgression and cloning in Brassicas.     

An avirulence gene,AvrLmJ1, from the blackleg fungus,Leptosphaeria maculans,confers avirulence to Brassica juncea cultivars

The fungus Leptosphaeria maculans causes blackleg of Brassica species. Here, we report the mapping and subsequent cloning of an avirulence gene from L. maculans. This gene, termed AvrLmJ1, confers avirulence towards all three Brassica juncea cultivars tested. Analysis of RNA‐seq data showed that AvrLmJ1 is housed in a region of the L. maculans genome which contains only one gene that is highly expressed in planta. The closest genes are 57 and 33 kb away and, like other avirulence genes of L. maculans, AvrLmJ1 is located within an AT‐rich, gene‐poor region of the genome. The encoded protein is 141 amino acids, has a predicted signal peptide and is cysteine rich. Two virulent isolates contain a premature stop codon in AvrLmJ1. Complementation of an isolate that forms cotyledonary lesions on B. juncea with the wild‐type allele of AvrLmJ1 confers avirulence towards all three B. juncea cultivars tested, suggesting that the gene may confer species‐specific avirulence activity.     

The Brassica napus receptor‐like protein RLM2 is encoded by a second allele of the LepR3/Rlm2 blackleg resistance locus

Leucine‐rich repeat receptor‐like proteins (LRR‐RLPs) are highly adaptable parts of the signalling apparatus for extracellular detection of plant pathogens. Resistance to blackleg disease of Brassica spp. caused by Leptosphaeria maculans is largely governed by host race‐specific R‐genes, including the LRR‐RLP gene LepR3. The blackleg resistance gene Rlm2 was previously mapped to the same genetic interval as LepR3. In this study, the LepR3 locus of the Rlm2 Brassica napus line ‘Glacier DH24287’ was cloned, and B. napus transformants were analysed for recovery of the Rlm2 phenotype. Multiple B. napus, B. rapa and B. juncea lines were assessed for sequence variation at the locus. Rlm2 was found to be an allelic variant of the LepR3 LRR‐RLP locus, conveying race‐specific resistance to L. maculans isolates harbouring AvrLm2. Several defence‐related LRR‐RLPs have previously been shown to associate with the RLK SOBIR1 to facilitate defence signalling. Bimolecular fluorescence complementation (BiFC) and co‐immunoprecipitation of RLM2‐SOBIR1 studies revealed that RLM2 interacts with SOBIR1 of Arabidopsis thaliana when co‐expressed in Nicotiana benthamiana. The interaction of RLM2 with AtSOBIR1 is suggestive of a conserved defence signalling pathway between B. napus and its close relative A. thaliana.     

Genetic diversity and population structure of Leptosphaeria maculans isolates in Western Canada

Leptosphaeria maculans is a serious concern for canola production worldwide. For effective disease management, knowledge of the pathogen's genetic variability and population structure is a prerequisite. In this study, whole-genome sequencing was performed for 162 of 1590 L. maculans isolates collected in the years 2007–2008 and 2012–2014 in Western Canada. DNA variants in genome-wide and specific regions including avirulence (Avr) genes were characterized. A total of 31,870 high-quality polymorphic DNA variants were used to study L. maculans genetic diversity and population structure. Cluster analysis showed that 150 isolates were clustered into 2 main groups and 4 subgroups by DNA variants located in either Avr or small secreted protein-encoding genes and into 2 main groups and 6 subgroups by genome-wide variants. The analysis of nucleotide diversity and differentiation also confirmed genetic variation within a population and among populations. Principal component analysis with genome-wide variants showed that the isolates collected in 2012–2014 were more genetically diverse than those collected in 2007–2008. Population structure analysis discovered three distinct sub-populations. Although isolates from Saskatchewan and Alberta were of similar genetic composition, Manitoba isolates were highly diverse. Genome-wide association study detected DNA variants in genes AvrLm4-7, Lema_T86300, and Lema_T86310 associated with the years of collection.     

Temperature and leaf wetness duration affect phenotypic expression of Rlm6-mediated resistance to Leptosphaeria maculans in Brassica napus

Near-isogenic Brassica napus lines carrying/lacking resistance gene Rlm6 were used to investigate the effects of temperature and leaf wetness duration on phenotypic expression of Rlm6-mediated resistance. Leaves were inoculated with ascospores or conidia of Leptosphaeria maculans carrying the effector gene AvrLm6. Incubation period to the onset of lesion development, number of lesions and lesion diameter were assessed. Symptomless growth of L. maculans from leaf lesions to stems was investigated using a green fluorescent protein (GFP) expressing isolate carrying AvrLm6. L. maculans produced large grey lesions on Darmor (lacking Rlm6) at 5-25 degrees C and DarmorMX (carrying Rlm6) at 25 degrees C, but small dark spots and 'green islands' on DarmorMX at 5-20 degrees C. With increasing temperature/wetness duration, numbers of lesions/spots generally increased. GFP-expressing L. maculans grew from leaf lesions down leaf petioles to stems on DarmorMX at 25 degrees C but not at 15 degrees C. We conclude that temperature and leaf wetness duration affect the phenotypic expression of Rlm6-mediated resistance in leaves and subsequent L. maculans spread down petioles to produce stem cankers.