New Avirulence Genes in the Phytopathogenic Fungus Leptosphaeria maculans

DArTseq-Based, High-Throughput Identification of Novel Molecular Markers for the Detection of Blackleg (Leptosphaeria Spp.) Resistance in Rapeseed

Blackleg disease, caused by Leptosphaeria spp. fungi, is one of the most important diseases of Brassica napus, responsible for severe yield losses worldwide. Blackleg resistance is controlled by major R genes and minor quantitative trait loci (QTL). Due to the high adaptation ability of the pathogen, R-mediated resistance can be easily broken, while the resistance mediated via QTL is believed to be more durable. Thus, the identification of novel molecular markers linked to blackleg resistance for B. napus breeding programs is essential. In this study, 183 doubled haploid (DH) rapeseed lines were assessed in field conditions for resistance to Leptosphaeria spp. Subsequently, DArTseq-based Genome-Wide Association Study (GWAS) was performed to identify molecular markers linked to blackleg resistance. A total of 133,764 markers (96,121 SilicoDArT and 37,643 SNP) were obtained. Finally, nine SilicoDArT and six SNP molecular markers were associated with plant resistance to Leptosphaeria spp. at the highest significance level, p < 0.001. Importantly, eleven of these fifteen markers were found within ten genes located on chromosomes A06, A07, A08, C02, C03, C06 and C08. Given the immune-related functions of the orthologues of these genes in Arabidopsis thaliana, the identified markers hold great promise for application in rapeseed breeding programs.

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.

Regulation of effector gene expression as concerted waves in Leptosphaeria maculans: a two-player game

Effector genes, encoding molecules involved in disease establishment, are concertedly expressed throughout the lifecycle of plant-pathogenic fungi. However, little is known about how effector gene expression is regulated. Since many effector genes are located in repeat-rich regions, the role of chromatin remodeling in their regulation was recently investigated, notably establishing that the repressive histone modification H3K9me3, deposited by KMT1, was involved in several fungal species including Leptosphaeria maculans. Nevertheless, previous data suggest that a second regulatory layer, probably involving a specific transcription factor (TF), might be required. In L. maculans, a Dothideomycete causing stem canker of oilseed rape, we identified the ortholog of Pf2, a TF belonging to the Zn2Cys6 fungal-specific family, and described as essential for pathogenicity and effector gene expression. We investigated its role together with KMT1, by inactivating and over-expressing LmPf2 in a wild-type strain and a ∆kmt1 mutant. Functional analyses of the corresponding transformants highlighted an essential role of LmPf2 in the establishment of pathogenesis and we found a major effect of LmPf2 on the induction of effector gene expression once KMT1 repression is lifted. Our results show, for the first time, a dual control of effector gene expression.

Transcriptomics of temperature-sensitive R gene-mediated resistance identifies a WAKL10 protein interaction network

Understanding temperature-sensitivity of R gene-mediated resistance against apoplastic pathogens is important for sustainable food production in the face of global warming. Here, we show that resistance of Brassica napus cotyledons against Leptosphaeria maculans was temperature-sensitive in introgression line Topas-Rlm7 but temperature-resilient in Topas-Rlm4. A set of 1,646 host genes was differentially expressed in Topas-Rlm4 and Topas-Rlm7 in response to temperature. Amongst these were three WAKL10 genes, including BnaA07g20220D, representing the temperature-sensitive Rlm7-1 allele and Rlm4. Network analysis identified a WAKL10 protein interaction cluster specifically for Topas-Rlm7 at 25 °C. Diffusion analysis of the Topas-Rlm4 network identified WRKY22 as a putative regulatory target of the ESCRT-III complex-associated protein VPS60.1, which belongs to the WAKL10 protein interaction community. Combined enrichment analysis of gene ontology terms considering gene expression and network data linked vesicle-mediated transport to defence. Thus, dysregulation of effector-triggered defence in Topas-Rlm7 disrupts vesicle-associated resistance against the apoplastic pathogen L. maculans.

Molecular Interactions Between Leptosphaeria maculans and Brassica Species

Canola is an important oilseed crop, providing food, feed, and fuel around the world. However, blackleg disease, caused by the ascomycete Leptosphaeria maculans, causes significant yield losses annually. With the recent advances in genomic technologies, the understanding of the Brassica napus-L. maculans interaction has rapidly increased, with numerous Avr and R genes cloned, setting this system up as a model organism for studying plant-pathogen associations. Although the B. napus-L. maculans interaction follows Flor's gene-for-gene hypothesis for qualitative resistance, it also puts some unique spins on the interaction. This review discusses the current status of the host-pathogen interaction and highlights some of the future gaps that need addressing moving forward.

In silico integration of disease resistance QTL, genes and markers with the Brassica juncea physical map

Brassica juncea (AABB), Indian mustard, is a source of disease resistance genes for a wide range of pathogens. The availability of reference genome sequences for B. juncea has made it possible to characterise the genomic structure and distribution of these disease resistance genes. Potentially functional disease resistance genes can be identified by co-localization with genetically mapped disease resistance quantitative trait loci (QTL). Here we identify and characterise disease resistance gene analogs (RGAs), including nucleotide-binding site-leucine-rich repeat (NLR), receptor-like kinase (RLK) and receptor-like protein (RLP) classes, and investigate their association with disease resistance QTL intervals. The molecular genetic marker sequences for four white rust (Albugo candida) disease resistance QTL, six blackleg (Leptosphaeria maculans) disease resistance QTL and BjCHI1, a gene cloned from B. juncea for hypocotyl rot disease, were extracted from previously published studies and used to compare with candidate RGAs. Our results highlight the complications for the identification of functional resistance genes, including the duplicated appearance of genetic markers for several resistance loci, including Ac2(t), AcB1-A4.1, AcB1-A5.1, Rlm6 and PhR2 in both the A and B genomes, due to the presence of homoeologous regions. Furthermore, the white rust loci, Ac2(t) and AcB1-A4.1, mapped to the same position on chromosome A04 and may be different alleles of the same gene. Despite these challenges, a total of nine candidate genomic regions hosting 14 RLPs, 28 NLRs and 115 RLKs were identified. This study facilitates the mapping and cloning of functional resistance genes for applications in crop improvement programs.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-022-01309-5.

A new family of structurally conserved fungal effectors displays epistatic interactions with plant resistance proteins

Recognition of a pathogen avirulence (AVR) effector protein by a cognate plant resistance (R) protein triggers a set of immune responses that render the plant resistant. Pathogens can escape this so-called Effector-Triggered Immunity (ETI) by different mechanisms including the deletion or loss-of-function mutation of the AVR gene, the incorporation of point mutations that allow recognition to be evaded while maintaining virulence function, and the acquisition of new effectors that suppress AVR recognition. The Dothideomycete Leptosphaeria maculans, causal agent of oilseed rape stem canker, is one of the few fungal pathogens where suppression of ETI by an AVR effector has been demonstrated. Indeed, AvrLm4-7 suppresses Rlm3- and Rlm9-mediated resistance triggered by AvrLm3 and AvrLm5-9, respectively. The presence of AvrLm4-7 does not impede AvrLm3 and AvrLm5-9 expression, and the three AVR proteins do not appear to physically interact. To decipher the epistatic interaction between these L. maculans AVR effectors, we determined the crystal structure of AvrLm5-9 and obtained a 3D model of AvrLm3, based on the crystal structure of Ecp11-1, a homologous AVR effector candidate from Fulvia fulva. Despite a lack of sequence similarity, AvrLm5-9 and AvrLm3 are structural analogues of AvrLm4-7 (structure previously characterized). Structure-informed sequence database searches identified a larger number of putative structural analogues among L. maculans effector candidates, including the AVR effector AvrLmS-Lep2, all produced during the early stages of oilseed rape infection, as well as among effector candidates from other phytopathogenic fungi. These structural analogues are named LARS (for Leptosphaeria AviRulence and Suppressing) effectors. Remarkably, transformants of L. maculans expressing one of these structural analogues, Ecp11-1, triggered oilseed rape immunity in several genotypes carrying Rlm3. Furthermore, this resistance could be suppressed by AvrLm4-7. These results suggest that Ecp11-1 shares a common activity with AvrLm3 within the host plant which is detected by Rlm3, or that the Ecp11-1 structure is sufficiently close to that of AvrLm3 to be recognized by Rlm3.

A new family of structurally conserved fungal effectors displays epistatic interactions with plant resistance proteins

Recognition of a pathogen avirulence (AVR) effector protein by a cognate plant resistance (R) protein triggers a set of immune responses that render the plant resistant. Pathogens can escape this so-called Effector-Triggered Immunity (ETI) by different mechanisms including the deletion or loss-of-function mutation of the AVR gene, the incorporation of point mutations that allow recognition to be evaded while maintaining virulence function, and the acquisition of new effectors that suppress AVR recognition. The Dothideomycete Leptosphaeria maculans, causal agent of oilseed rape stem canker, is one of the few fungal pathogens where suppression of ETI by an AVR effector has been demonstrated. Indeed, AvrLm4-7 suppresses Rlm3- and Rlm9-mediated resistance triggered by AvrLm3 and AvrLm5-9, respectively. The presence of AvrLm4-7 does not impede AvrLm3 and AvrLm5-9 expression, and the three AVR proteins do not appear to physically interact. To decipher the epistatic interaction between these L. maculans AVR effectors, we determined the crystal structure of AvrLm5-9 and obtained a 3D model of AvrLm3, based on the crystal structure of Ecp11-1, a homologous AVR effector candidate from Fulvia fulva. Despite a lack of sequence similarity, AvrLm5-9 and AvrLm3 are structural analogues of AvrLm4-7 (structure previously characterized). Structure-informed sequence database searches identified a larger number of putative structural analogues among L. maculans effector candidates, including the AVR effector AvrLmS-Lep2, all produced during the early stages of oilseed rape infection, as well as among effector candidates from other phytopathogenic fungi. These structural analogues are named LARS (for Leptosphaeria AviRulence and Suppressing) effectors. Remarkably, transformants of L. maculans expressing one of these structural analogues, Ecp11-1, triggered oilseed rape immunity in several genotypes carrying Rlm3. Furthermore, this resistance could be suppressed by AvrLm4-7. These results suggest that Ecp11-1 shares a common activity with AvrLm3 within the host plant which is detected by Rlm3, or that the Ecp11-1 structure is sufficiently close to that of AvrLm3 to be recognized by Rlm3.

Two independent approaches converge to the cloning of a new Leptosphaeria maculans avirulence effector gene, AvrLmS-Lep2

Brassica napus (oilseed rape, canola) seedling resistance to Leptosphaeria maculans, the causal agent of blackleg (stem canker) disease, follows a gene-for-gene relationship. The avirulence genes AvrLmS and AvrLep2 were described to be perceived by the resistance genes RlmS and LepR2, respectively, present in B. napus 'Surpass 400'. Here we report cloning of AvrLmS and AvrLep2 using two independent methods. AvrLmS was cloned using combined in vitro crossing between avirulent and virulent isolates with sequencing of DNA bulks from avirulent or virulent progeny (bulked segregant sequencing). AvrLep2 was cloned using a biparental cross of avirulent and virulent L. maculans isolates and a classical map-based cloning approach. Taking these two approaches independently, we found that AvrLmS and AvrLep2 are the same gene. Complementation of virulent isolates with this gene confirmed its role in inducing resistance on Surpass 400, Topas-LepR2, and an RlmS-line. The gene, renamed AvrLmS-Lep2, encodes a small cysteine-rich protein of unknown function with an N-terminal secretory signal peptide, which is a common feature of the majority of effectors from extracellular fungal plant pathogens. The AvrLmS-Lep2/LepR2 interaction phenotype was found to vary from a typical hypersensitive response through intermediate resistance sometimes towards susceptibility, depending on the inoculation conditions. AvrLmS-Lep2 was nevertheless sufficient to significantly slow the systemic growth of the pathogen and reduce the stem lesion size on plant genotypes with LepR2, indicating the potential efficiency of this resistance to control the disease in the field.

Efficacy of Blackleg Major Resistance Genes in B. napus in Germany

Leptosphaeria maculans is one of the major pathogens of oilseed rape (B. napus). It causes blackleg disease, which accounts for significant yield losses worldwide. Using cultivars that harbor major resistance (R) genes is one of the most effective control methods. However, the efficacy of major R genes is related to the frequency of the corresponding avirulence (Avr) genes in a L. maculans population. In this paper, we report the Avr profiles of L. maculans populations and the ratio of its mating types in Northern and Central regions of Germany. Eleven Avr genes in five-hundred and seventy-four isolates were characterized either by applying cotyledon tests on a B. napus differential set or by amplifying avirulence gene-specific PCR markers. Fifty-two races were determined, among which the most dominant race was Avrlm6, -7, -11, AvrlepR1, -R2. Results showed that the resistance gene Rlm2 is 100% ineffective, some other major R genes such as Rlm1, Rlm3, Rlm4 and LepR3 are partially effective (with corresponding Avr frequencies ≤ 42%), while LepR1, LepR2, Rlm6, Rlm11 and Rlm7 can still provide relatively effective resistance in the German fields investigated (with corresponding Avr frequencies of 63-100%). Sexual reproduction is a factor that enhances the potential of L. maculans to evolve under selection pressure. Mating types of the L. maculans populations did not deviate from the ratio of 1:1 in the examined regions, indicating that sexual reproduction and ascospores play central roles in the L. maculans lifecycle. Overall, this study provides an important dataset for the establishment of a strategic plan to preserve the efficacies of major R genes in Germany by applying cultivar rotations of oilseed rape.

A chromosome-scale genome assembly of the tomato pathogen Cladosporium fulvum reveals a compartmentalized genome architecture and the presence of a dispensable chromosome

Cladosporium fulvum is a fungal pathogen that causes leaf mould of tomato. The reference genome of this pathogen was released in 2012 but its high repetitive DNA content prevented a contiguous assembly and further prohibited the analysis of its genome architecture. In this study, we combined third generation sequencing technology with the Hi-C chromatin conformation capture technique, to produce a high-quality and near complete genome assembly and gene annotation of a Race 5 isolate of C. fulvum. The resulting genome assembly contained 67.17 Mb organized into 14 chromosomes (Chr1-to-Chr14), all of which were assembled telomere-to-telomere. The smallest of the chromosomes, Chr14, is only 460 kb in size and contains 25 genes that all encode hypothetical proteins. Notably, PCR assays revealed that Chr14 was absent in 19 out of 24 isolates of a world-wide collection of C. fulvum, indicating that Chr14 is dispensable. Thus, C. fulvum is currently the second species of Capnodiales shown to harbour dispensable chromosomes. The genome of C. fulvum Race 5 is 49.7 % repetitive and contains 14 690 predicted genes with an estimated completeness of 98.9%, currently one of the highest among the Capnodiales. Genome structure analysis revealed a compartmentalized architecture composed of gene-dense and repeat-poor regions interspersed with gene-sparse and repeat-rich regions. Nearly 39.2 % of the C. fulvum Race 5 genome is affected by Repeat-Induced Point (RIP) mutations and evidence of RIP leakage toward non-repetitive regions was observed in all chromosomes, indicating the RIP plays an important role in the evolution of this pathogen. Finally, 345 genes encoding candidate effectors were identified in C. fulvum Race 5, with a significant enrichment of their location in gene-sparse regions, in accordance with the 'two-speed genome' model of evolution. Overall, the new reference genome of C. fulvum presents several notable features and is a valuable resource for studies in plant pathogens.

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.

Dissection of Quantitative Blackleg Resistance Reveals Novel Variants of Resistance Gene Rlm9 in Elite Brassica napus

Blackleg is one of the major fungal diseases in oilseed rape/canola worldwide. Most commercial cultivars carry R gene-mediated qualitative resistances that confer a high level of race-specific protection against Leptosphaeria maculans, the causal fungus of blackleg disease. However, monogenic resistances of this kind can potentially be rapidly overcome by mutations in the pathogen's avirulence genes. To counteract pathogen adaptation in this evolutionary arms race, there is a tremendous demand for quantitative background resistance to enhance durability and efficacy of blackleg resistance in oilseed rape. In this study, we characterized genomic regions contributing to quantitative L. maculans resistance by genome-wide association studies in a multiparental mapping population derived from six parental elite varieties exhibiting quantitative resistance, which were all crossed to one common susceptible parental elite variety. Resistance was screened using a fungal isolate with no corresponding avirulence (AvrLm) to major R genes present in the parents of the mapping population. Genome-wide association studies revealed eight significantly associated quantitative trait loci (QTL) on chromosomes A07 and A09, with small effects explaining 3-6% of the phenotypic variance. Unexpectedly, the qualitative blackleg resistance gene Rlm9 was found to be located within a resistance-associated haploblock on chromosome A07. Furthermore, long-range sequence data spanning this haploblock revealed high levels of single-nucleotide and structural variants within the Rlm9 coding sequence among the parents of the mapping population. The results suggest that novel variants of Rlm9 could play a previously unknown role in expression of quantitative disease resistance in oilseed rape.

Remote homology clustering identifies lowly conserved families of effector proteins in plant-pathogenic fungi

Plant diseases caused by fungal pathogens are typically initiated by molecular interactions between 'effector' molecules released by a pathogen and receptor molecules on or within the plant host cell. In many cases these effector-receptor interactions directly determine host resistance or susceptibility. The search for fungal effector proteins is a developing area in fungal-plant pathology, with more than 165 distinct confirmed fungal effector proteins in the public domain. For a small number of these, novel effectors can be rapidly discovered across multiple fungal species through the identification of known effector homologues. However, many have no detectable homology by standard sequence-based search methods. This study employs a novel comparison method (RemEff) that is capable of identifying protein families with greater sensitivity than traditional homology-inference methods, leveraging a growing pool of confirmed fungal effector data to enable the prediction of novel fungal effector candidates by protein family association. Resources relating to the RemEff method and data used in this study are available from https://figshare.com/projects/Effector_protein_remote_homology/87965.

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.

Validating the Strategic Deployment of Blackleg Resistance Gene Groups in Commercial Canola Fields on the Canadian Prairies

Blackleg, caused by the fungal pathogen Leptosphaeria maculans, is a serious threat to canola (Brassica napus L.) production in western Canada. Crop scouting and extended crop rotation, along with the use of effective genetic resistance, have been key management practices available to mitigate the impact of the disease. In recent years, new pathogen races have reduced the effectiveness of some of the resistant cultivars deployed. Strategic deployment and rotation of major resistance (R) genes in cultivars have been used in France and Australia to help increase the longevity of blackleg resistance. Canada also introduced a grouping system in 2017 to identify blackleg R genes in canola cultivars. The main objective of this study was to examine and validate the concept of R gene deployment through monitoring the avirulence (Avr) profile of L. maculans population and disease levels in commercial canola fields within the Canadian prairies. Blackleg disease incidence and severity was collected from 146 cultivars from 53 sites across Manitoba, Saskatchewan, and Alberta in 2018 and 2019, and the results varied significantly between gene groups, which is likely influenced by the pathogen population. Isolates collected from spring and fall stubble residues were examined for the presence of Avr alleles AvrLm1, AvrLm2, AvrLm3, AvrLm4, AvrLm5, AvrLm6, AvrLm7, AvrLm9, AvrLm10, AvrLm11, AvrLepR1, AvrLepR2, AvrLep3, and AvrLmS using a set of differential host genotypes carrying known resistance genes or PCR-based markers. The Simpson's evenness index was very low, due to two dominant L. maculans races (AvrLm2-4-5-6-7-10-11 and AvrLm2-5-6-7-10-11) representing 49% of the population, but diversity of the population was high from the 35 L. maculans races isolated in Manitoba. AvrLm6 and AvrLm11 were found in all 254 L. maculans isolates collected in Manitoba. Knowledge of the blackleg disease levels in relation to the R genes deployed, along with the L. maculans Avr profile, helps to measure the effectiveness of genetic resistance.

Leptosphaeria maculans Isolates Reveal Their Allele Frequency in Western Canada

Blackleg, caused by Leptosphaeria maculans, is a major disease of canola in Canada, Australia, and Europe. For effective deployment of resistant varieties and disease management, it is crucial to understand the population structure of L. maculans. In this study, we analyzed L. maculans isolates from commercial fields in western Canada from 2014 to 2016 for the presence and frequency of avirulence (Avr) genes. A total of 1,584 isolates were examined for the presence of Avr genes AvrLm1, AvrLm2, AvrLm3, AvrLm4, AvrLm6, AvrLm7, AvrLm9, AvrLepR1, AvrLepR2, and AvrLmS via a set of differential host genotypes carrying known resistance genes and a PCR assay. Several Avr genes showed a higher frequency in the pathogen population, such as AvrLm6 and AvrLm7, which were present in >90% of isolates, whereas AvrLm3, AvrLm9, and AvrLepR2 showed frequencies of <10%. A total of 189 races (different combinations of Avr genes) were detected, with Avr-2-4-6-7-S, Avr-1-4-6-7, and Avr-2-4-6-7 as the three predominant races. When the effect of crop rotation was assessed, only a 3-year rotation showed a significantly higher frequency of AvrLm2 relative to shorter rotations. This study provides the information for producers to select effective canola varieties for blackleg management and for breeders to deploy new R genes in disease resistance breeding in western Canada.

Large-scale transcriptomics to dissect 2 years of the life of a fungal phytopathogen interacting with its host plant

BACKGROUND

The fungus Leptosphaeria maculans has an exceptionally long and complex relationship with its host plant, Brassica napus, during which it switches between different lifestyles, including asymptomatic, biotrophic, necrotrophic, and saprotrophic stages. The fungus is also exemplary of "two-speed" genome organisms in the genome of which gene-rich and repeat-rich regions alternate. Except for a few stages of plant infection under controlled conditions, nothing is known about the genes mobilized by the fungus throughout its life cycle, which may last several years in the field.

RESULTS

We performed RNA-seq on samples corresponding to all stages of the interaction of L. maculans with its host plant, either alive or dead (stem residues after harvest) in controlled conditions or in field experiments under natural inoculum pressure, over periods of time ranging from a few days to months or years. A total of 102 biological samples corresponding to 37 sets of conditions were analyzed. We show here that about 9% of the genes of this fungus are highly expressed during its interactions with its host plant. These genes are distributed into eight well-defined expression clusters, corresponding to specific infection lifestyles or to tissue-specific genes. All expression clusters are enriched in effector genes, and one cluster is specific to the saprophytic lifestyle on plant residues. One cluster, including genes known to be involved in the first phase of asymptomatic fungal growth in leaves, is re-used at each asymptomatic growth stage, regardless of the type of organ infected. The expression of the genes of this cluster is repeatedly turned on and off during infection. Whatever their expression profile, the genes of these clusters are enriched in heterochromatin regions associated with H3K9me3 or H3K27me3 repressive marks. These findings provide support for the hypothesis that part of the fungal genes involved in niche adaptation is located in heterochromatic regions of the genome, conferring an extreme plasticity of expression.

CONCLUSION

This work opens up new avenues for plant disease control, by identifying stage-specific effectors that could be used as targets for the identification of novel durable disease resistance genes, or for the in-depth analysis of chromatin remodeling during plant infection, which could be manipulated to interfere with the global expression of effector genes at crucial stages of plant infection.

Genome-Wide Identification and Analysis of the Valine-Glutamine Motif-Containing Gene Family in Brassica napus and Functional Characterization of BnMKS1 in Response to Leptosphaeria maculans

Proteins containing valine-glutamine (VQ) motifs play important roles in plant growth and development as well as in defense responses to both abiotic and biotic stresses. Blackleg disease, which is caused by Leptosphaeria maculans, is the most important disease in canola (Brassica napus) worldwide; however, the identification of Brassica napus VQs and their functions in response to blackleg disease have not yet been reported. In this study, we conducted a genome-wide identification and characterization of the VQ gene family in Brassica napus, including chromosome location, phylogenetic relations, gene structure, motif domain, synteny analysis, and cis-elements categorization of their promoter regions. To understand Brassica napus VQ gene function in response to blackleg disease, we overexpressed BnVQ7 (BnaA01g36880D, also known as the mitogen-activated protein kinase 4 substrate 1 [MKS1] gene) in a blackleg-susceptible canola variety, Westar. Overexpression of BnMKS1 in canola did not improve its resistance to blackleg disease at the seedling stage; however, transgenic canola plants overexpressing BnMKS1 displayed an enhanced resistance to L. maculans infection at the adult plant stage. Expression levels of downstream and defense marker genes in cotyledons increased significantly at the necrotrophic stage of L. maculans infection in the overexpression line of BnMKS1, suggesting that the salicylic acid- and jasmonic acid-mediated signaling pathways were both involved in the defense responses. Together, these results suggest that BnMKS1 might play an important role in defense against L. maculans.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Recent Findings Unravel Genes and Genetic Factors Underlying Leptosphaeria maculans Resistance in Brassica napus and Its Relatives

Among the Brassica oilseeds, canola (Brassica napus) is the most economically significant globally. However, its production can be limited by blackleg disease, caused by the fungal pathogen Lepstosphaeria maculans. The deployment of resistance genes has been implemented as one of the key strategies to manage the disease. Genetic resistance against blackleg comes in two forms: qualitative resistance, controlled by a single, major resistance gene (R gene), and quantitative resistance (QR), controlled by numerous, small effect loci. R-gene-mediated blackleg resistance has been extensively studied, wherein several genomic regions harbouring R genes against L. maculans have been identified and three of these genes were cloned. These studies advance our understanding of the mechanism of R gene and pathogen avirulence (Avr) gene interaction. Notably, these studies revealed a more complex interaction than originally thought. Advances in genomics help unravel these complexities, providing insights into the genes and genetic factors towards improving blackleg resistance. Here, we aim to discuss the existing R-gene-mediated resistance, make a summary of candidate R genes against the disease, and emphasise the role of players involved in the pathogenicity and resistance. The comprehensive result will allow breeders to improve resistance to L. maculans, thereby increasing yield.

Development of a specific marker for detection of a functional AvrLm9 allele and validating the interaction between AvrLm7 and AvrLm9 in Leptosphaeria maculans

Blackleg, which is caused by the fungus Leptosphaeria maculans (L. maculans), is a major disease of canola in western Canada and worldwide. Long-term use of one source of resistance could cause the breakdown of its effectiveness. Therefore, appropriate use of R genes is very important, and knowledge about the distribution of avirulence genes is a prerequisite for effectively deploying resistance. Of the 14 avirulence genes identified in L. maculans, AvrLm5 and AvrLm9 were recognized as the two alleles of the same gene based on two single nucleotide polymorphisms, C⁸⁵T and G¹⁶⁴A/C. In this study, a specific marker was developed to identify AvrLm5 and AvrLm9 based on two single nucleotide polymorphisms, C⁸⁵T and G¹⁶⁴A/C, which are responsible for the function of AvrLm9. The specific marker can be used to discriminate the AvrLm9 from avrLm9 accurately in L. maculans isolates, which is consistent with inoculation tests in isolates without AvrLm4-7. This specific marker was used to screen 1229 isolates collected from fields in the years 2014 through 2016 in Manitoba. From 68 to 84% of the isolates were found to contain the AvrLm9 allele; while 4-7% of them were avirulent on the variety Goéland with Rlm9 loci. Furthermore, no isolates having both AvrLm9 and AvrLm7 were detected using a cotyledon test, while 67% to 84% of isolates contained both avirulence genes via PCR detection, implying suppression of AvrLm9 by AvrLm7. In addition, avirulence gene profiles of the other 10 avirulence alleles were examined with the 1229 isolates using cotyledon tests or PCR amplifications. Taken together, this research enables the fast identification of AvrLm5/9, provides the Avr genes' landscape of western Canada and elaborates the relationship between AvrLm9 and AvrLm7 using isolates from grower fields.

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.

Generation and Characterization of a Virulent Leptosphaeria maculans Isolate Carrying a Mutated AvrLm7 Gene Using the CRISPR/Cas9 System

Blackleg, caused by the fungal pathogen Leptosphaeria maculans, is the most important disease affecting canola (Brassica napus) crops worldwide. We employed the clustered regularly interspaced short palindromic repeat (CRISPR)/CRISPR-associated (Cas) system to generate the mutant isolate umavr7 from a point mutation of the AvrLm7 coding region in a L. maculans isolate (UMAvr7). Reverse transcription PCR and transcriptome data confirmed that the AvrLm7 gene was knocked out in the mutant isolate. Pathogenicity tests indicated that umavr7 can cause large lesions on a set of Brassica differential genotypes that express different resistance (R) genes. Comparative pathogenicity tests between UMAvr7 (wild type) and umavr7 on the corresponding B. napus genotype 01-23-2-1 (with Rlm7) showed that umavr7 is a mutant isolate, producing large gray/green lesions on cotyledons. The pathogenicity of the mutant isolate was shifted from avirulent to virulent on the B. napus Rlm7 genotype. Therefore, this mutant is virulence on the identified resistant genes to blackleg disease in B. napus genotypes. Superoxide accumulated differently in cotyledons in response to infection with UMAvr7 and umavr7, especially in resistant B. napus genotype 01-23-2-1. Resistance/susceptibility was further evaluated on 123 B. napus genotypes with the mutant isolate, umavr7. Only 6 of the 123 genotypes showed resistance to umavr7. The identification of these six resistant B. napus genotypes will lead to further studies on the development of blackleg disease resistance through breeding and the identification of novel R genes.

LmCBP1, a secreted chitin-binding protein, is required for the pathogenicity of Leptosphaeria maculans on Brassica napus

Leptosphaeria maculans is the causal agent of blackleg disease on Brassica napus. Determining the underlying functions of genes required for pathogenesis is essential for understanding the infection process. A chitin-binding protein (LmCBP1) was discovered as a pathogenicity factor for the infection of B. napus by L. maculans through gene knockout using the CRISPR-Cas9 system. Chitin-binding activity was demonstrated through a chitin-protein binding assay. A secreted signal peptide was detected using a yeast secreted-signal peptide trap assay. An increased expression level during the infection stage was also observed, suggesting that LmCBP1 is a secreted protein. The knockout mutants showed decreased infection on B. napus, with reduced pathogenicity on ten cultivars with/without diverse R genes. The mutants were more sensitive to H₂O₂ compared to wild type L. maculans isolate JN3. This study provides evidence of the virulence of a novel chitin-binding protein LmCBP1 on B. napus through mutants created via the CRISPR-Cas9 system.

Inheritance Pattern and Molecular Markers for Resistance to Blackleg Disease in Cabbage

The inheritance and causal loci for resistance to blackleg, a devastating disease of Brassicaceous crops, are yet to be known in cabbage (Brassica oleracea L.). Here, we report the pattern of inheritance and linked molecular marker for this trait. A segregating BC₁ population consisting of 253 plants was raised from resistant and susceptible parents, L29 (♀) and L16 (♂), respectively. Cotyledon resistance bioassay of BC₁ population, measured based on a scale of 0-9 at 12 days after inoculation with Leptosphaeria maculans isolate 03-02 s, revealed the segregation of resistance and ratio, indicative of dominant monogenic control of the trait. Investigation of potential polymorphism in the previously identified differentially expressed genes within the collinear region of 'B. napus blackleg resistant loci Rlm1' in B. oleracea identified two insertion/deletion (InDel) mutations in the intron and numerous single nucleotide polymorphisms (SNPs) throughout the LRR-RLK gene Bol040029, of which six SNPs in the first exon caused the loss of two LRR domains in the susceptible line. An InDel marker, BLR-C-InDel based on the InDel mutations, and a high resolution melting (HRM) marker, BLR-C-2808 based on the SNP C²⁸⁰⁸T in the second exon were developed, which predicated the resistance status of the BC₁ population with 80.24%, and of 24 commercial inbred lines with 100% detection accuracy. This is the first report of inheritance and molecular markers linked with blackleg resistance in cabbage. This study will enhance our understanding of the trait, and will be helpful in marker assisted breeding aiming at developing resistant cabbage varieties.

Effect of Elevated CO2 Concentration on the Disease Severity of Compatible and Incompatible Interactions of Brassica napus-Leptosphaeria maculans Pathosystem

Global warming by increased atmospheric CO2 concentration has been widely accepted. Yet, there has not been any consistent conclusion on the doubled CO2 concentration that in the future will affect plant disease incidence and severity. Blackleg disease, mainly caused by Leptosphaeria maculans, is a major disease on canola production globally. Brassica napus and L. maculans have a gene-for-gene interaction, which causes an incompatible reaction between canola plants carrying resistance genes and L. maculans isolates carrying corresponding avirulence genes. In this study, B. napus varieties and lines inoculated with different Leptosphaeria isolates were subjected to simulated growth conditions, namely, growth chambers with normal environments and with controlled CO2 concentrations of 400, 600, and 800 ppm. The results indicated that the elevated CO2 concentrations have no noticeable effect on the inferred phenotypes of the canola-blackleg interactions. However, the disease severity decreased in most of the B. napus-L. maculans interactions at extremely high CO2 concentration (800 ppm). The varied pathogenicity changes of the B. napus-L. maculans pathosystem under elevated CO2 concentrations at 400 or 600 ppm may be due to the genetic background or physiological differences in plants and pathogenicity differences in L. maculans isolates having different Avr gene profiles. The mechanisms by which elevated CO2 concentrations affect the B. napus-L. maculans pathosystem will help us understand how climate change will impact crops and diseases.

Characterization of the Race Structure of Leptosphaeria maculans Causing Blackleg of Winter Canola in Oklahoma and Kansas

Blackleg, caused by the fungus Leptosphaeria maculans, is a widespread disease of winter canola (Brassica napus) in Oklahoma and Kansas. Deployment of genetic resistance is the primary strategy for managing blackleg. Resistance genes (Rlm) in canola interact with avirulence genes in the fungus (AvrLm) in a gene-for-gene manner. Little is known about the diversity and frequency of avirulence genes and the race structure in the region. Isolates of Leptosphaeria spp. were collected from diseased leaves in nine counties in Oklahoma and one county in Kansas from 2009 to 2013. Based on pathogenicity and PCR amplification of mating type and species-specific internal transcribed spacer loci, most isolates (n = 90) were L. maculans. The presence of avirulence genes was evaluated using phenotypic interactions on cotyledons of differential cultivars with Rlm1, Rlm2, Rlm3, and Rlm4 and amplification of AvrLm1, AvrLm4-7, and AvrLm6 by PCR. The avirulence alleles AvrLm6 and AvrLm7 were present in the entire L. maculans population. AvrLm1 was found in 34% of the population, AvrLm2 in 4%, and AvrLm4 in only 1%. A total of five races, defined as combinations of avirulence alleles, were identified that included AvrLm1-2-6-7, AvrLm2-6-7, AvrLm4-6-7, AvrLm1-6-7, and AvrLm6-7. Races virulent on the most Rlm genes, AvrLm1-6-7 at 32% and AvrLm6-7 at 62%, were predominant. Defining the avirulence allele frequency and race structure of L. maculans should be useful for the identification and development of resistant cultivars and hybrids for blackleg management in the region. The results suggest that Rlm6 and Rlm7 would be effective, although their deployment should be integrated with quantitative resistance and cultural practices, such as crop rotation, that limit selection pressure on Rlm genes.

A two genes - for - one gene interaction between Leptosphaeria maculans and Brassica napus

Interactions between Leptosphaeria maculans, causal agent of stem canker of oilseed rape, and its Brassica hosts are models of choice to explore the multiplicity of 'gene-for-gene' complementarities and how they diversified to increased complexity in the course of plant-pathogen co-evolution. Here, we support this postulate by investigating the AvrLm10 avirulence that induces a resistance response when recognized by the Brassica nigra resistance gene Rlm10. Using genome-assisted map-based cloning, we identified and cloned two AvrLm10 candidates as two genes in opposite transcriptional orientation located in a subtelomeric repeat-rich region of the genome. The AvrLm10 genes encode small secreted proteins and show expression profiles in planta similar to those of all L. maculans avirulence genes identified so far. Complementation and silencing assays indicated that both genes are necessary to trigger Rlm10 resistance. Three assays for protein-protein interactions showed that the two AvrLm10 proteins interact physically in vitro and in planta. Some avirulence genes are recognized by two distinct resistance genes and some avirulence genes hide the recognition specificities of another. Our L. maculans model illustrates an additional case where two genes located in opposite transcriptional orientation are necessary to induce resistance. Interestingly, orthologues exist for both L. maculans genes in other phytopathogenic species, with a similar genome organization, which may point to an important conserved effector function linked to heterodimerization of the two proteins.

A Novel Effector Gene SCRE2 Contributes to Full Virulence of Ustilaginoidea virens to Rice

Ustilaginoidea virens, the causal agent of rice false smut (RFS), has become one of the most devastating rice pathogens worldwide. As a group of essential virulence factors, the effectors in the filamentous fungus might play central roles in the interaction between plants and pathogens. However, little is known about the roles of individual effectors in U. virens virulence. In this study, we identified and characterized a small secreted cysteine-rich effector, SCRE2, in U. virens. SCRE2 was first confirmed as an effector through yeast secretion, protein localization and translocation assays, as well as its expression pattern during U. virens infection. Transient expression of SCRE2 in Nicotiana benthamiana suppressed necrosis-like defense symptoms triggered by the mammalian BAX and oomycete elicitin INF1 proteins. The ability of SCRE2 to inhibit immunity-associated responses in N. benthamiana, including elicitor-triggered cell death and oxidative burst, is further defined to a small peptide region SCRE2^68-85 through expressing a series of truncated proteins. Convincingly, ectopic expression of SCRE2 in the transgenic rice cells significantly inhibited pathogen-associated molecular pattern-triggered immunity including flg22- and chitin-induced defense gene expression and oxidative burst. Furthermore, the scre2 knockout mutant generated by the CRISPR/Cas9 system greatly attenuated in U. virens virulence to rice. Collectively, this study indicates that the effector SCRE2 is able to inhibit plant immunity and is required for full virulence of U. virens.

Characterization of Callose Deposition and Analysis of the Callose Synthase Gene Family of Brassica napus in Response to Leptosphaeria maculans

Callose plays a critical role in different biological processes including development as well as in the response to multiple biotic and abiotic stresses. In this study, we characterized the callose deposition in cotyledons of different Brassica napus varieties post-inoculated with different Leptosphaeria maculans isolates. Further, members of the callose synthase gene were identified from the whole genome of B. napus using the 12 Arabidopsis thaniana callose synthase protein sequences, and were then classified into three groups based on their phylogenetic relationships. Chromosomal location and duplication patterns indicated uneven distribution and segmental duplication patterns of BnCalS genes in the B. napus genome. Subsequently, gene structures, conserved domains analysis, and protein properties were analyzed for BnCalS genes. In addition, 12 B. napus orthologs of the AtCalS were selected for investigating the tissue expression pattern, indicating diverse expression patterns for these BnCalS genes. Responses of the selected 12 orthologs and all the BnCalS genes were characterized in the different types (AvrLm1-Rlm1, AvrLm4-Rlm4, AvrLepR1-LepR1) of B. napus–L. maculans interactions and B. napus-Leptosphaeria biglobosa interactions, implying their potential roles in response to Leptosphaeria infection.

Effectiveness of Rlm7 resistance against Leptosphaeria maculans (phoma stem canker) in UK winter oilseed rape cultivars

The Rlm7 gene in Brassica napus is an important source of resistance for control of phoma stem canker on oilseed rape caused by the fungus Leptosphaeria maculans. This study shows the first report of L. maculans isolates virulent against Rlm7 in the UK. Leptosphaeria maculans isolates virulent against Rlm7 represented 3% of the pathogen population when cultivars with the Rlm7 gene represented 5% of the UK oilseed rape area in 2012/13. However, the Rlm7 gene has been widely used since then, representing >15% of the UK oilseed rape area in 2015/16. Winter oilseed rape field experiments included cultivars with the Rlm7 gene, with the Rlm4 gene or without Rlm genes and took place at five sites in the UK over four cropping seasons. An increase in phoma leaf spotting severity on Rlm7 cultivars in successive seasons was observed. Major resistance genes played a role in preventing severe phoma leaf spotting at the beginning of the cropping season and, in addition, quantitative resistance (QR) in the cultivars examined made an important contribution to control of phoma stem canker development at the end of the cropping season. Deployment of the Rlm7 resistance gene against L. maculans in cultivars with QR in combination with sustainable disease management practices will prolong the use of this gene for effective control of phoma stem canker epidemics.

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.

Current Status and Challenges in Identifying Disease Resistance Genes in Brassica napus

Brassica napus is an economically important crop across different continents including temperate and subtropical regions in Europe, Canada, South Asia, China and Australia. Its widespread cultivation also brings setbacks as it plays host to fungal, oomycete and chytrid pathogens that can lead to serious yield loss. For sustainable crop production, identification of resistance (R) genes in B. napus has become of critical importance. In this review, we discuss four key pathogens affecting Brassica crops: Clubroot (Plasmodiophora brassicae), Blackleg (Leptosphaeria maculans and L. biglobosa), Sclerotinia Stem Rot (Sclerotinia sclerotiorum), and Downy Mildew (Hyaloperonospora parasitica). We first review current studies covering prevalence of these pathogens on Brassica crops and highlight the R genes and QTL that have been identified from Brassica species against these pathogens. Insights into the relationships between the pathogen and its Brassica host, the unique host resistance mechanisms and how these affect resistance outcomes is also presented. We discuss challenges in identification and deployment of R genes in B. napus in relation to highly specific genetic interactions between host subpopulations and pathogen pathotypes and emphasize the need for common or shared techniques and research materials or tighter collaboration between researchers to reconcile the inconsistencies in the research outcomes. Using current genomics tools, we provide examples of how characterization and cloning of R genes in B. napus can be carried out more effectively. Lastly, we put forward strategies to breed resistant cultivars through introgressions supported by genomic approaches and suggest prospects that can be implemented in the future for a better, pathogen-resistant B. napus.

Different waves of effector genes with contrasted genomic location are expressed by Leptosphaeria maculans during cotyledon and stem colonization of oilseed rape

Leptosphaeria maculans, the causal agent of stem canker disease, colonizes oilseed rape (Brassica napus) in two stages: a short and early colonization stage corresponding to cotyledon or leaf colonization, and a late colonization stage during which the fungus colonizes systemically and symptomlessly the plant during several months before stem canker appears. To date, the determinants of the late colonization stage are poorly understood; L. maculans may either successfully escape plant defences, leading to stem canker development, or the plant may develop an 'adult-stage' resistance reducing canker incidence. To obtain an insight into these determinants, we performed an RNA-sequencing (RNA-seq) pilot project comparing fungal gene expression in infected cotyledons and in symptomless or necrotic stems. Despite the low fraction of fungal material in infected stems, sufficient fungal transcripts were detected and a large number of fungal genes were expressed, thus validating the feasibility of the approach. Our analysis showed that all avirulence genes previously identified are under-expressed during stem colonization compared with cotyledon colonization. A validation RNA-seq experiment was then performed to investigate the expression of candidate effector genes during systemic colonization. Three hundred and seven 'late' effector candidates, under-expressed in the early colonization stage and over-expressed in the infected stems, were identified. Finally, our analysis revealed a link between the regulation of expression of effectors and their genomic location: the 'late' effector candidates, putatively involved in systemic colonization, are located in gene-rich genomic regions, whereas the 'early' effector genes, over-expressed in the early colonization stage, are located in gene-poor regions of the genome.

Unusual evolutionary mechanisms to escape effector-triggered immunity in the fungal phytopathogen Leptosphaeria maculans

Leptosphaeria maculans is the fungus responsible for the stem canker disease of oilseed rape (Brassica napus). AvrLm3 and AvrLm4-7, two avirulence effector genes of L. maculans, are involved in an unusual relationship: AvrLm4-7 suppresses the Rlm3-mediated resistance. Here, we assessed AvrLm3 polymorphism in a collection of 235 L. maculans isolates. No field isolates exhibited deletion or inactivating mutations in AvrLm3, as observed for other L. maculans avirulence genes. Eleven isoforms of the AvrLm3 protein were found. In isolates virulent towards both Rlm3 and Rlm7 (a3a7), the loss of the Rlm3-mediated resistance response was due to two distinct mechanisms. First, when AvrLm4-7 was inactivated (deletion or inactivating mutations), amino acid substitutions in AvrLm3 generated virulent isoforms of the protein. Second, when only point mutations were observed in AvrLm4-7, a3a7 isolates still contained an avirulent allele of AvrLm3. Directed mutagenesis confirmed that some point mutations in AvrLm4-7 were sufficient for the fungus to escape Rlm7-mediated resistance while maintaining the suppression of the AvrLm3 phenotype. Signatures of positive selection were also identified in AvrLm3. The complex evolutionary mechanisms enabling L. maculans to escape Rlm3-mediated resistance while preserving AvrLm3 integrity, along with observed reduced aggressiveness of isolates silenced for AvrLm3, serves to emphasize the importance of this effector in pathogenicity towards B. napus. While the common response to resistance gene pressure is local selection of isolates depleted in the cognate avirulence gene, this example contributes to complexify the gene-for-gene concept of plant-pathogen evolution with a 'camouflaged' model allowing retention of nondispensable avirulence effectors.

Impact of biotic and abiotic factors on the expression of fungal effector-encoding genes in axenic growth conditions

In phytopathogenic fungi, the expression of hundreds of small secreted protein (SSP)-encoding genes is induced upon primary infection of plants while no or a low level of expression is observed during vegetative growth. In some species such as Leptosphaeria maculans, this coordinated in-planta upregulation of SSP-encoding genes expression relies on an epigenetic control but the signals triggering gene expression in-planta are unknown. In the present study, biotic and abiotic factors that may relieve suppression of SSP-encoding gene expression during axenic growth of L. maculans were investigated. Some abiotic factors (temperature, pH) could have a limited effect on SSP gene expression. In contrast, two types of cellular stresses induced by antibiotics (cycloheximide, phleomycin) activated strongly the transcription of SSP genes. A transcriptomic analysis to cycloheximide exposure revealed that biological processes such as ribosome biosynthesis and rRNA processing were induced whereas important metabolic pathways such as glycogen and nitrogen metabolism, glycolysis and tricarboxylic acid cycle activity were down-regulated. A quantitatively different expression of SSP-encoding genes compared to plant infection was also detected. Interestingly, the same physico-chemical parameters as those identified here for L. maculans effectors were identified to regulate positively or negatively the expression of bacterial effectors. This suggests that apoplastic phytopathogens may react to similar physiological parameters for regulation of their effector genes.

Single R Gene Introgression Lines for Accurate Dissection of the Brassica - Leptosphaeria Pathosystem

Seven blackleg resistance (R) genes (Rlm1, Rlm2, Rlm3, Rlm4, LepR1, LepR2 & LepR3) were each introgressed into a common susceptible B. napus doubled-haploid (DH) line through reciprocal back-crossing, producing single-R gene introgression lines (ILs) for use in the pathological and molecular study of Brassica-Leptosphaeria interactions. The genomic positions of the R genes were defined through molecular mapping and analysis with transgenic L. maculans isolates was used to confirm the identity of the introgressed genes where possible. Using L. maculans isolates of contrasting avirulence gene (Avr) profiles, we preformed extensive differential pathology for phenotypic comparison of the ILs to other B. napus varieties, demonstrating the ILs can provide for the accurate assessment of Avr-R gene interactions by avoiding non-Avr dependant alterations to resistance responses which can occur in some commonly used B. napus varieties. Whole-genome SNP-based assessment allowed us to define the donor parent introgressions in each IL and provide a strong basis for comparative molecular dissection of the pathosystem.

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.

A Large Family of AvrLm6-like Genes in the Apple and Pear Scab Pathogens, Venturia inaequalis and Venturia pirina

Venturia inaequalis and V. pirina are Dothideomycete fungi that cause apple scab and pear scab disease, respectively. Whole genome sequencing of V. inaequalis and V. pirina isolates has revealed predicted proteins with sequence similarity to AvrLm6, a Leptosphaeria maculans effector that triggers a resistance response in Brassica napus and B. juncea carrying the resistance gene, Rlm6. AvrLm6-like genes are present as large families (>15 members) in all sequenced strains of V. inaequalis and V. pirina, while in L. maculans, only AvrLm6 and a single paralog have been identified. The Venturia AvrLm6-like genes are located in gene-poor regions of the genomes, and mostly in close proximity to transposable elements, which may explain the expansion of these gene families. An AvrLm6-like gene from V. inaequalis with the highest sequence identity to AvrLm6 was unable to trigger a resistance response in Rlm6-carrying B. juncea. RNA-seq and qRT-PCR gene expression analyses, of in planta- and in vitro-grown V. inaequalis, has revealed that many of the AvrLm6-like genes are expressed during infection. An AvrLm6 homolog from V. inaequalis that is up-regulated during infection was shown (using an eYFP-fusion protein construct) to be localized to the sub-cuticular stroma during biotrophic infection of apple hypocotyls.

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 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.

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.

Co-localisation of the blackleg resistance genes Rlm2 and LepR3 on Brassica napus chromosome A10

BACKGROUND

The protection of canola (Brassica napus) crops against blackleg disease, caused by the fungal pathogen Leptosphaeria maculans, is largely mediated by race-specific resistance genes (R-genes). While many R-genes effective against blackleg disease have been identified in Brassica species, information of the precise genomic locations of the genes is limited.

RESULTS

In this study, the Rlm2 gene for resistance to blackleg, located on chromosome A10 of the B. napus cultivar 'Glacier', was targeted for fine mapping. Molecular markers tightly linked to the gene were developed for use in mapping the resistance locus and defining the physical interval in B. napus. Rlm2 was localised to a 5.8 cM interval corresponding to approximately 873 kb of the B. napus chromosome A10.

CONCLUSION

The recently-cloned B. napus R-gene, LepR3, occupies the same region of A10 as Rlm2 and analysis of the putative B. napus and B. rapa genes in the homologous region identified several additional candidate defense-related genes that may control Rlm2 function.

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.

Molecular cytogenetic identification of B genome chromosomes linked to blackleg disease resistance in Brassica napus × B. carinata interspecific hybrids

Provide evidence that the Brassica B genome chromosome B3 carries blackleg resistance gene, and also the B genome chromosomes were inherited several generations along with B. napus chromosomes. Blackleg disease caused by fungus Leptosphaeria maculans causes significant yield losses in Brassica napus. Brassica carinata possesses excellent resistance to this disease. To introgress blackleg resistance, crosses between B. napus cv. Westar and B. carinata were done. The interspecific-hybrids were backcrossed twice to Westar and self-pollinated three times to produce BC2S3 families. Doubled haploid lines (DH1) were produced from one blackleg resistant family. SSR markers were used to study the association between B genome chromosome(s) and blackleg resistance. The entire B3 chromosome of B. carinata was associated with blackleg resistance in DH1. A second DH population (DH2) was produced from F1s of resistant DH1 lines crossed to blackleg susceptible B. napus cv. Polo where resistance was found to be associated with SSR markers from the middle to bottom of the B3 and top of the B8 chromosomes. The results demonstrated that the B3 chromosome carried gene(s) for blackleg resistance. Genomic in situ hybridization (GISH) and GISH-like analysis of the DH2 lines revealed that susceptible lines, in addition to B. napus chromosomes, possessed one pair of B genome chromosomes (2n = 40), while resistant lines had either one (2n = 40) or two pairs (2n = 42) of B chromosomes. The molecular and GISH data suggested that the B chromosome in the susceptible lines was B7, while it was difficult to confirm the identity of the B chromosomes in the resistant lines. Also, B chromosomes were found to be inherited over several generations along with B. napus chromosomes.

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.