Estimating frequencies of virulent isolates in field populations of a plant pathogenic fungus, Leptosphaeria maculans, using high-throughput pyrosequencing

Twenty Years of Leptosphaeria maculans Population Survey in France Suggests Pyramiding Rlm3 and Rlm7 in Rapeseed Is a Risky Resistance Management Strategy

Strategies for plant resistance gene deployment aim to preserve their durability to highly adaptable fungal pathogens. While the pyramiding of resistance genes is often proposed as an effective way to increase their durability, molecular mechanisms by which the pathogen can overcome the resistance also are important aspects to take into account. Here, we report a counterexample where pyramiding of two resistance genes of Brassica napus, Rlm3 and Rlm7, matching the Leptosphaeria maculans avirulence genes AvrLm3 and AvrLm4-7, respectively, favored the selection of double-virulent isolates. We previously demonstrated that the presence of a functional AvrLm4-7 gene in an isolate masks the Rlm3-AvrLm3 recognition. Rlm7 was massively deployed in France since 2004. L. maculans populations were surveyed on a large scale (>7,600 isolates) over a period of 20 years, and resistance gene deployment at the regional scale was determined. Mutations in isolates overcoming both resistance genes were analyzed. All data indicated that the simultaneous success of Rlm7, the deployment of varieties pyramiding Rlm3 and Rlm7, along with the decrease in areas cultivated with Rlm3 only, contributed to the success of virulent isolates toward Rlm7, and more recently to both Rlm3 and Rlm7. Experimental field assays proved that resistance gene alternation was a better strategy compared with pyramiding in this context. Our study also illustrated an unusually high sequence diversification of AvrLm3 and AvrLm4-7 under such a selection pressure, and identified a few regions of the AvrLm4-7 protein involved in both its recognition by Rlm7 and in its AvrLm3-Rlm3 masking ability. [Formula: see text] Copyright © 2022 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

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.

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

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

Comparison of non-subjective relative fungal biomass measurements to quantify the Leptosphaeria maculans-Brassica napus interaction

BACKGROUND

Blackleg disease, caused by the fungal pathogen Leptosphaeria maculans, is a serious threat to canola (Brassica napus) production worldwide. Quantitative resistance to this disease is a highly desirable trait but is difficult to precisely phenotype. Visual scores can be subjective and are prone to assessor bias. Methods to assess variation in quantitative resistance more accurately were developed based on quantifying in planta fungal biomass, including the Wheat Germ Agglutinin Chitin Assay (WAC), qPCR and ddPCR assays.

RESULTS

Disease assays were conducted by inoculating a range of canola cultivars with L. maculans isolates in glasshouse experiments and assessing fungal biomass in cotyledons, petioles and stem tissue harvested at different timepoints post-inoculation. PCR and WAC assay results were well correlated, repeatable across experiments and host tissues, and able to differentiate fungal biomass in different host-isolate treatments. In addition, the ddPCR assay was shown to differentiate between L. maculans isolates.

CONCLUSIONS

The ddPCR assay is more sensitive in detecting pathogens and more adaptable to high-throughput methods by using robotic systems than the WAC assay. Overall, these methods proved accurate and non-subjective, providing alternatives to visual assessments to quantify the L. maculans-B. napus interaction in all plant tissues throughout the progression of the disease in seedlings and mature plants and have potential for fine-scale blackleg resistance phenotyping in canola.

Technological Advances in Phytopathogen Detection and Metagenome Profiling Techniques

The use of advanced molecular methods in plant pathology and applied microbiology has necessitated for more accurate, rapid detection and identification of plant pathogens. This is particularly significant given accelerated emergence of virulence that leads to increased prevalence of plant pathogens. Thus, the capacity to contain plant pathogens and ultimately disease progression is key to ensuring crop biosecurity and overall food security. Of recent, research on pathogens utilizes a holistic approach focusing on elucidating growth dynamics within the entire biome rather than studying individual or closely related isolates in unison. This has advanced knowledge and information of microbial ecosystem within natural environments in the twenty first century. Applied technological platforms used for rapid detection and profiling microbial biomes in this regard include digital PCR, pyrosequencing, Illumina, DNA microarray and barcoding, Ion torrent, and nanopore. These technologies have been applied in various fields including human health and medicine, marine and animal biology, crop production and water quality research, to mention but a few. Although much has been done and achieved through the development of several technologies, more accuracy is required to circumvent the shortfalls still experienced. This includes integrating existing methods with new applications such as viability PCRs and microbial viability testing. Hence, this review provides critical analysis of some widely used latest technologies in rapid detection and identification of plant pathogens, and profiling plant associated microbiomes that reveal growth dynamics and population diversity. The advantages and limitations of the technologies are also discussed.

Analysis of Repeat Induced Point (RIP) Mutations in Leptosphaeria maculans Indicates Variability in the RIP Process Between Fungal Species

Gene duplication contributes to evolutionary potential, yet many duplications in a genome arise from the activity of "selfish" genetic elements such as transposable elements. Fungi have a number of mechanisms by which they limit the expansion of transposons, including Repeat Induced Point mutation (RIP). RIP has been best characterized in the Sordariomycete Neurospora crassa, wherein duplicated DNA regions are recognized after cell fusion, but before nuclear fusion during the sexual cycle, and then mutated. While "signatures" of RIP appear in the genome sequences of many fungi, the species most distant from N. crassa in which the process has been experimentally demonstrated to occur is the Dothideomycete Leptosphaeria maculans In the current study, we show that similar to N. crassa, nonlinked duplications can trigger RIP; however, the frequency of the generated RIP mutations is extremely low in L maculans (< 0.1%) and requires a large duplication to initiate RIP, and that multiple premeiotic mitoses are involved in the RIP process. However, a single sexual cycle leads to the generation of progeny with unique haplotypes, despite progeny pairs being generated from mitosis. We hypothesize that these different haplotypes may be the result of the deamination process occurring post karyogamy, leading to unique mutations within each of the progeny pairs. These findings indicate that the RIP process, while common to many fungi, differs between fungi and that this impacts on the fate of duplicated DNA.

Genome scans on experimentally evolved populations reveal candidate regions for adaptation to plant resistance in the potato cyst nematode Globodera pallida

Improving resistance durability involves to be able to predict the adaptation speed of pathogen populations. Identifying the genetic bases of pathogen adaptation to plant resistances is a useful step to better understand and anticipate this phenomenon. Globodera pallida is a major pest of potato crop for which a resistance QTL, GpaV_vrn , has been identified in Solanum vernei. However, its durability is threatened as G. pallida populations are able to adapt to the resistance in few generations. The aim of this study was to investigate the genomic regions involved in the resistance breakdown by coupling experimental evolution and high-density genome scan. We performed a whole-genome resequencing of pools of individuals (Pool-Seq) belonging to G. pallida lineages derived from two independent populations having experimentally evolved on susceptible and resistant potato cultivars. About 1.6 million SNPs were used to perform the genome scan using a recent model testing for adaptive differentiation and association to population-specific covariables. We identified 275 outliers and 31 of them, which also showed a significant reduction in diversity in adapted lineages, were investigated for their genic environment. Some candidate genomic regions contained genes putatively encoding effectors and were enriched in SPRYSECs, known in cyst nematodes to be involved in pathogenicity and in (a)virulence. Validated candidate SNPs will provide a useful molecular tool to follow frequencies of virulence alleles in natural G. pallida populations and define efficient strategies of use of potato resistances maximizing their durability.

A pyrosequencing-based method to quantify genetic substitutions associated with resistance to succinate dehydrogenase inhibitor fungicides in Botrytis spp. populations

BACKGROUND

The genetic underlying resistance mechanisms in the population of the phytopathogenic fungus Botrytis cinerea are well documented. Specifically, several genetic substitutions associated with succinate dehydrogenase inhibitor (SDHI)-based fungicide resistance have been identified in the succinate dehydrogenase gene. The objective of the present work was to develop a molecular tool for accurate quantification of these genetic substitutions within Botrytis populations. A test using the PyroMark Q24 instrument was designed to detect and quantify five genetic substitutions associated with SDHI resistance.

RESULTS

The technique is based on sequencing by synthesis, and it generated quantitative and accurate data with a limit of quantification of a minimum of 500 spores. There was a linear relationship between the known and estimated percentages of spores with the targeted genetic substitutions and wild-type strains at ratios of 0-100%, with a 20% increment.

CONCLUSION

With the pyrosequencing assay developed in this study, a large number of Botrytis spp. individuals can be characterised in a timely fashion with greater accuracy than by commonly used methods. Hence, pyrosequencing-based methods will be useful for improving our understanding of fungicide resistance, detecting the arrival of new genetic substitutions, monitoring shifts in fungal populations and assessing the effectiveness of antiresistance strategies, and for routine monitoring of fungicide resistance.

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.

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.

Epigenetic control of effectors in plant pathogens

Plant pathogens display impressive versatility in adapting to host immune systems. Pathogen effector proteins facilitate disease but can become avirulence (Avr) factors when the host acquires discrete recognition capabilities that trigger immunity. The mechanisms that lead to changes to pathogen Avr factors that enable escape from host immunity are diverse, and include epigenetic switches that allow for reuse or recycling of effectors. This perspective outlines possibilities of how epigenetic control of Avr effector gene expression may have arisen and persisted in filamentous plant pathogens, and how it presents special problems for diagnosis and detection of specific pathogen strains or pathotypes.

Eukaryotic diversity in premise drinking water using 18S rDNA sequencing: implications for health risks

The goal of this study was to characterize microbial eukaryotes over a 12-month period to provide insight into the occurrence of potential bacterial predators and hosts in premise plumbing. Nearly 6,300 partial 18S rRNA gene sequences from 24 hot (36.9-39.0 °C) and cold (6.8-29.1 °C) drinking water samples were analyzed and classified into major eukaryotic groups. Each major group, consisting of free-living amoebae (FLA)/protozoa, algae, copepods, dinoflagellates, fungi, nematodes, and unique uncultured eukaryotic sequences, showed limited diversity dominated by a few distinct populations, which may be characteristic of oligotrophic environments. Changes in the relative abundance of predators such as nematodes, copepods, and FLA appear to be related to temperature and seasonal changes in water quality. Sequences nearly identical to FLA such as Hartmannella vermiformis, Echinamoeba thermarmum, Pseudoparamoeba pagei, Protacanthamoeba bohemica, Platyamoeba sp., and Vannella sp. were obtained. In addition to FLA, various copepods, rotifers, and nematodes have been reported to internalize viral and bacterial pathogens within drinking water systems thus potentially serving as transport hosts; implications of which are discussed further. Increasing the knowledge of eukaryotic occurrence and their relationship with potential pathogens should aid in assessing microbial risk associated with various eukaryotic organisms in drinking water.

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

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

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

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