Genetic mapping of the Leptosphaeria maculans avirulence gene corresponding to the LepR1 resistance gene of Brassica napus

Breeding and management of major resistance genes to stem canker/blackleg in Brassica crops

Blackleg (also known as Phoma or stem canker) is a major, worldwide disease of Brassica crop species, notably B. napus (rapeseed, canola), caused by the ascomycete fungus Leptosphaeria maculans. The outbreak and severity of this disease depend on environmental conditions and management practices, as well as a complex interaction between the pathogen and its hosts. Genetic resistance is a major method to control the disease (and the only control method in some parts of the world, such as continental Europe), but efficient use of genetic resistance is faced with many difficulties: (i) the scarcity of germplasm/genetic resources available, (ii) the different history of use of resistance genes in different parts of the world and the different populations of the fungus the resistance genes are exposed to, (iii) the complexity of the interactions between the plant and the pathogen that expand beyond typical gene-for-gene interactions, (iv) the incredible evolutionary potential of the pathogen and the importance of knowing the molecular processes set up by the fungus to "breakdown' resistances, so that we may design high-throughput diagnostic tools for population surveys, and (v) the different strategies and options to build up the best resistances and to manage them so that they are durable. In this paper, we aim to provide a comprehensive overview of these different points, stressing the differences between the different continents and the current prospects to generate new and durable resistances to blackleg disease.

A gene-for-gene interaction involving a 'late' effector contributes to quantitative resistance to the stem canker disease in Brassica napus

The control of stem canker disease of Brassica napus (rapeseed), caused by the fungus Leptosphaeria maculans is based largely on plant genetic resistance: single-gene specific resistance (Rlm genes) or quantitative, polygenic, adult-stage resistance. Our working hypothesis was that quantitative resistance partly obeys the gene-for-gene model, with resistance genes 'recognizing' fungal effectors expressed during late systemic colonization. Five LmSTEE (stem-expressed effector) genes were selected and placed under the control of the AvrLm4-7 promoter, an effector gene highly expressed at the cotyledon stage of infection, for miniaturized cotyledon inoculation test screening of a gene pool of 204 rapeseed genotypes. We identified a rapeseed genotype, 'Yudal', expressing hypersensitive response to LmSTEE98. The LmSTEE98-RlmSTEE98 interaction was further validated by inactivation of the LmSTEE98 gene with a CRISPR-Cas9 approach. Isolates with mutated versions of LmSTEE98 induced more severe stem symptoms than the wild-type isolate in 'Yudal'. This single-gene resistance was mapped in a 0.6 cM interval of the 'Darmor_bzh' × 'Yudal' genetic map. One typical gene-for-gene interaction contributes partly to quantitative resistance when L. maculans colonizes the stems of rapeseed. With numerous other effectors specific to stem colonization, our study provides a new route for resistance gene discovery, elucidation of quantitative resistance mechanisms and selection for durable resistance.

Candidate Rlm6 resistance genes against Leptosphaeria. maculans identified through a genome-wide association study in Brassica juncea (L.) Czern

One hundred and sixty-seven B. juncea varieties were genotyped on the 90K Brassica assay (42,914 SNPs), which led to the identification of sixteen candidate genes for Rlm6. Brassica species are at high risk of severe crop loss due to pathogens, especially Leptosphaeria maculans (the causal agent of blackleg). Brassica juncea (L.) Czern is an important germplasm resource for canola improvement, due to its good agronomic traits, such as heat and drought tolerance and high blackleg resistance. The present study is the first using genome-wide association studies to identify candidate genes for blackleg resistance in B. juncea based on genome-wide SNPs obtained from the Illumina Infinium 90 K Brassica SNP array. The verification of Rlm6 in B. juncea was performed through a cotyledon infection test. Genotyping 42,914 single nucleotide polymorphisms (SNPs) in a panel of 167 B. juncea lines revealed a total of seven SNPs significantly associated with Rlm6 on chromosomes A07 and B04 in B. juncea. Furthermore, 16 candidate Rlm6 genes were found in these regions, defined as nucleotide binding site leucine-rich-repeat (NLR), leucine-rich repeat RLK (LRR-RLK) and LRR-RLP genes. This study will give insights into the blackleg resistance in B. juncea and facilitate identification of functional blackleg resistance genes which can be used in Brassica breeding.

The Rlm13 Gene, a New Player of Brassica napus-Leptosphaeria maculans Interaction Maps on Chromosome C03 in Canola

Canola exhibits an extensive genetic variation for resistance to blackleg disease, caused by the fungal pathogen Leptosphaeria maculans. Despite the identification of several Avr effectors and R (race-specific) genes, specific interactions between Avr-R genes are not yet fully understood in the Brassica napus-L. maculans pathosystem. In this study, we investigated the genetic basis of resistance in an F_{2 : 3} population derived from Australian canola varieties CB-Telfer (Rlm4)/ATR-Cobbler (Rlm4) using a single-spore isolate of L. maculans, PHW1223. A genetic linkage map of the CB-Telfer/ATR-Cobbler population was constructed using 7,932 genotyping-by-sequencing-based DArTseq markers and subsequently utilized for linkage and haplotype analyses. Genetic linkage between DArTseq markers and resistance to PHW1223 isolate was also validated using the B. napus 60K Illumina Infinium array. Our results revealed that a major locus for resistance, designated as Rlm13, maps on chromosome C03. To date, no R gene for resistance to blackleg has been reported on the C subgenome in B. napus. Twenty-four candidate R genes were predicted to reside within the quantitative trait locus (QTL) region. We further resequenced both the parental lines of the mapping population (CB-Telfer and ATR-Cobbler, > 80 × coverage) and identified several structural sequence variants in the form of single-nucleotide polymorphisms (SNPs), insertions/deletions (InDels), and presence/absence variations (PAVs) near Rlm13. Comparative mapping revealed that Rlm13 is located within the homoeologous A03/C03 region in ancestral karyotype block "R" of Brassicaceae. Our results provide a "target" for further understanding the Avr-Rlm13 gene interaction as well as a valuable tool for increasing resistance to blackleg in canola germplasm.

Meta-analysis of GWAS in canola blackleg (Leptosphaeria maculans) disease traits demonstrates increased power from imputed whole-genome sequence

Blackleg disease causes yield losses in canola (Brassica napus L.). To identify resistance genes and genomic regions, genome-wide association studies (GWAS) of 585 diverse winter and spring canola accessions were performed using imputed whole-genome sequence (WGS) and transcriptome genotype-by-sequencing (GBSt). Blackleg disease phenotypes were collected across three years in six trials. GWAS were performed in several ways and their respective power was judged by the number of significant single nucleotide polymorphisms (SNP), the false discovery rate (FDR), and the percentage of SNP that validated in additional field trials in two subsequent years. WGS GWAS with 1,234,708 million SNP detected a larger number of significant SNP, achieved a lower FDR and a higher validation rate than GBSt with 64,072 SNP. A meta-analysis combining survival and average internal infection resulted in lower FDR but also lower validation rates. The meta-analysis GWAS identified 79 genomic regions (674 SNP) conferring potential resistance to L. maculans. While several GWAS signals localised in regions of known Rlm genes, fifty-three new potential resistance regions were detected. Seventeen regions had underlying genes with putative functions related to disease defence or stress response in Arabidopsis thaliana. This study provides insight into the genetic architecture and potential molecular mechanisms underlying canola L. maculans resistance.

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 effector AvrLm4-7 affects salicylic acid (SA) and ethylene (ET) signalling and hydrogen peroxide (H2 O2 ) accumulation in Brassica napus

To achieve host colonization, successful pathogens need to overcome plant basal defences. For this, (hemi)biotrophic pathogens secrete effectors that interfere with a range of physiological processes of the host plant. AvrLm4-7 is one of the cloned effectors from the hemibiotrophic fungus Leptosphaeria maculans 'brassicaceae' infecting mainly oilseed rape (Brassica napus). Although its mode of action is still unknown, AvrLm4-7 is strongly involved in L. maculans virulence. Here, we investigated the effect of AvrLm4-7 on plant defence responses in a susceptible cultivar of B. napus. Using two isogenic L. maculans isolates differing in the presence of a functional AvrLm4-7 allele [absence ('a4a7') and presence ('A4A7') of the allele], the plant hormone concentrations, defence-related gene transcription and reactive oxygen species (ROS) accumulation were analysed in infected B. napus cotyledons. Various components of the plant immune system were affected. Infection with the 'A4A7' isolate caused suppression of salicylic acid- and ethylene-dependent signalling, the pathways regulating an effective defence against L. maculans infection. Furthermore, ROS accumulation was decreased in cotyledons infected with the 'A4A7' isolate. Treatment with an antioxidant agent, ascorbic acid, increased the aggressiveness of the 'a4a7' L. maculans isolate, but not that of the 'A4A7' isolate. Together, our results suggest that the increased aggressiveness of the 'A4A7' L. maculans isolate could be caused by defects in ROS-dependent defence and/or linked to suppressed SA and ET signalling. This is the first study to provide insights into the manipulation of B. napus defence responses by an effector of L. maculans.

Elucidating the Role of Effectors in Plant-Fungal Interactions: Progress and Challenges

Pathogenic fungi have diverse growth lifestyles that support fungal colonization on plants. Successful colonization and infection for all lifestyles depends upon the ability to modify living host plants to sequester the necessary nutrients required for growth and reproduction. Secretion of virulence determinants referred to as “effectors” is assumed to be the key governing factor that determines host infection and colonization. Effector proteins are capable of suppressing plant defense responses and alter plant physiology to accommodate fungal invaders. This review focuses on effector molecules of biotrophic and hemibiotrophic plant pathogenic fungi, and the mechanism required for the release and uptake of effector molecules by the fungi and plant cells, respectively. We also place emphasis on the discovery of effectors, difficulties associated with predicting the effector repertoire, and fungal genomic features that have helped promote effector diversity leading to fungal evolution. We discuss the role of specific effectors found in biotrophic and hemibiotrophic fungi and examine how CRISPR/Cas9 technology may provide a new avenue for accelerating our ability in the discovery of fungal effector function.

A game of hide and seek between avirulence genes AvrLm4-7 and AvrLm3 in Leptosphaeria maculans

Extending the durability of plant resistance genes towards fungal pathogens is a major challenge. We identified and investigated the relationship between two avirulence genes of Leptosphaeria maculans, AvrLm3 and AvrLm4-7. When an isolate possesses both genes, the Rlm3-mediated resistance of oilseed rape (Brassica napus) is not expressed due to the presence of AvrLm4-7 but virulent isolates toward Rlm7 recover the AvrLm3 phenotype. Combining genetic and genomic approaches (genetic mapping, RNA-seq, BAC (bacterial artificial chromosome) clone sequencing and de novo assembly) we cloned AvrLm3, a telomeric avirulence gene of L. maculans. AvrLm3 is located in a gap of the L. maculans reference genome assembly, is surrounded by repeated elements, encodes for a small secreted cysteine-rich protein and is highly expressed at early infection stages. Complementation and silencing assays validated the masking effect of AvrLm4-7 on AvrLm3 recognition by Rlm3 and we showed that the presence of AvrLm4-7 does not impede AvrLm3 expression in planta. Y2H assays suggest the absence of physical interaction between the two avirulence proteins. This unusual interaction is the basis for field experiments aiming to evaluate strategies that increase Rlm7 durability.

Comparative Transcriptomic Analysis of Virulence Factors in Leptosphaeria maculans during Compatible and Incompatible Interactions with Canola

Leptosphaeria maculans is a hemibiotrophic fungus that causes blackleg of canola (Brassica napus), one of the most devastating diseases of this crop. In the present study, transcriptome profiling of L. maculans was performed in an effort to understand and define the pathogenicity genes that govern both the biotrophic and the necrotrophic phase of the fungus, as well as those that separate a compatible from an incompatible interaction. For this purpose, comparative RNA-seq analyses were performed on L. maculans isolate D5 at four different time points following inoculation on susceptible cultivar Topas-DH16516 or resistant introgression line Topas-Rlm2. Analysis of 1.6 billion Illumina reads readily identified differentially expressed genes that were over represented by candidate secretory effector proteins, CAZymes, and other pathogenicity genes. Comparisons between the compatible and incompatible interactions led to the identification of 28 effector proteins whose chronology and level of expression suggested a role in the establishment and maintenance of biotrophy with the plant. These included all known Avr genes of isolate D5 along with eight newly characterized effectors. In addition, another 15 effector proteins were found to be exclusively expressed during the necrotrophic phase of the fungus, which supports the concept that L. maculans has a separate and distinct arsenal contributing to each phase. As for CAZymes, they were often highly expressed at 3 dpi but with no difference in expression between the compatible and incompatible interactions, indicating that other factors were necessary to determine the outcome of the interaction. However, their significantly higher expression at 11 dpi in the compatible interaction confirmed that they contributed to the necrotrophic phase of the fungus. A notable exception was LysM genes whose high expression was singularly observed on the susceptible host at 7 dpi. In the case of TFs, their higher expression at 7 and 11 dpi on susceptible Topas support an important role in regulating the genes involved in the different pathogenic phases of L. maculans. In conclusion, comparison of the transcriptome of L. maculans during compatible and incompatible interactions has led to the identification of key pathogenicity genes that regulate not only the fate of the interaction but also lifestyle transitions of the fungus.

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.

The Brassica napus blackleg resistance gene LepR3 encodes a receptor-like protein triggered by the Leptosphaeria maculans effector AVRLM1

LepR3, found in the Brassica napus cv 'Surpass 400', provides race-specific resistance to the fungal pathogen Leptosphaeria maculans, which was overcome after great devastation in Australia in 2004. We investigated the LepR3 locus to identify the genetic basis of this resistance interaction. We employed a map-based cloning strategy, exploiting collinearity with the Arabidopsis thaliana and Brassica rapa genomes to enrich the map and locate a candidate gene. We also investigated the interaction of LepR3 with the L. maculans avirulence gene AvrLm1 using transgenics. LepR3 was found to encode a receptor-like protein (RLP). We also demonstrated that avirulence towards LepR3 is conferred by AvrLm1, which is responsible for both the Rlm1 and LepR3-dependent resistance responses in B. napus. LepR3 is the first functional B. napus disease resistance gene to be cloned. AvrLm1's interaction with two independent resistance loci, Rlm1 and LepR3, highlights the need to consider redundant phenotypes in 'gene-for-gene' interactions and offers an explanation as to why LepR3 was overcome so rapidly in parts of Australia.