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.     

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.     

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.     

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.     

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.     

Independent breakdown events of the Brassica napus Rlm7 resistance gene including via the off-target impact of a dual-specificity avirulence interaction

Protection of many crops is achieved through the use of genetic resistance. Leptosphaeria maculans, the causal agent of blackleg disease of Brassica napus, has emerged as a model for understanding gene-for-gene interactions that occur between plants and pathogens. Whilst many of the characterized avirulence effector genes interact with a single resistance gene in the host, the AvrLm4-7 avirulence gene is recognized by two resistance genes, Rlm4 and Rlm7. Here, we report the “breakdown” of the Rlm7 resistance gene in Australia, under two different field conditions. The first, and more typical, breakdown probably resulted from widescale use of Rlm7-containing cultivars whereby selection has led to an increase of individuals in the L. maculans population that have undergone repeat-induced point (RIP) mutations at the AvrLm4-7 locus. This has rendered the AvrLm4-7 gene ineffective and therefore these isolates have become virulent towards both Rlm4 and Rlm7. The second, more atypical, situation was the widescale use of Rlm4 cultivars. Whilst a single-nucleotide polymorphism is the more common mechanism of virulence towards Rlm4, in this field situation, RIP mutations have been selected leading to the breakdown of resistance for both Rlm4 and Rlm7. This is an example of a resistance gene being rendered ineffective without having grown cultivars with the corresponding resistance gene due to the dual specificity of the avirulence gene. These findings highlight the value of pathogen surveillance in the context of expanded knowledge about potential complexities for Avr–R interactions for the deployment of appropriate resistance gene strategies.     

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.     

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.     

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.     

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.     

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.     

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.     

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.     

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.     

Identifying genetic diversity of avirulence genes in Leptosphaeria maculans using whole genome sequencing

Next generation sequencing technology allows rapid re-sequencing of individuals, as well as the discovery of single nucleotide polymorphisms (SNPs), for genomic diversity and evolutionary analyses. By sequencing two isolates of the fungal plant pathogen Leptosphaeria maculans, the causal agent of blackleg disease in Brassica crops, we have generated a resource of over 76 million sequence reads aligned to the reference genome. We identified over 21,000 SNPs with an overall SNP frequency of one SNP every 2,065 bp. Sequence validation of a selection of these SNPs in additional isolates collected throughout Australia indicates a high degree of polymorphism in the Australian population. In preliminary phylogenetic analysis, isolates from Western Australia clustered together and those collected from Brassica juncea stubble were identical. These SNPs provide a novel marker resource to study the genetic diversity of this pathogen. We demonstrate that re-sequencing provides a method of validating previously characterised SNPs and analysing differences in important genes, such as the disease related avirulence genes of L. maculans. Understanding the genetic characteristics of this devastating pathogen is vital in developing long-term solutions to managing blackleg disease in Brassica crops.     

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.