Microsatellite and Minisatellite Analysis of Leptosphaeria maculans in Australia Reveals Regional Genetic Differentiation

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

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

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

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

Genetic and physical mapping of loci for resistance to blackleg disease in canola (Brassica napus L.)

Sustainable canola production is essential to meet growing human demands for vegetable oil, biodiesel, and meal for stock feed markets. Blackleg, caused by the fungal pathogen, Leptosphaeria maculans is a devastating disease that can lead to significant yield loss in many canola production regions worldwide. Breakdown of race-specific resistance to L. maculans in commercial cultivars poses a constant threat to the canola industry. To identify new alleles, especially for quantitative resistance (QR), we analyzed 177 doubled haploid (DH) lines derived from an RP04/Ag-Outback cross. DH lines were evaluated for QR under field conditions in three experiments conducted at Wagga Wagga (2013, 2014) and Lake Green (2015), and under shade house conditions using the 'ascospore shower' test. DH lines were also characterized for qualitative R gene-mediated resistance via cotyledon tests with two differential single spore isolates, IBCN17 and IBCN76, under glasshouse conditions. Based on 18,851 DArTseq markers, a linkage map representing 2,019 unique marker bins was constructed and then utilized for QTL detection. Marker regression analysis identified 22 significant marker associations for resistance, allowing identification of two race-specific resistance R genes, Rlm3 and Rlm4, and 21 marker associations for QR loci. At least three SNP associations for QR were repeatedly detected on chromosomes A03, A07 and C04 across phenotyping environments. Physical mapping of markers linked with these consistent QR loci on the B. napus genome assembly revealed their localization in close proximity of the candidate genes of B. napus BnaA03g26760D (A03), BnaA07g20240D (A07) and BnaC04g02040D (C04). Annotation of these candidate genes revealed their association with protein kinase and jumonji proteins implicated in defense resistance. Both Rlm3 and Rlm4 genes identified in this DH population did not show any association with resistance loci detected under either field and/or shade house conditions (ascospore shower) suggesting that both genes are ineffective in conferring resistance to L. maculans in Australian field conditions. Taken together, our study identified sequence-based molecular markers for dissecting R and QR loci to L. maculans in a canola DH population from the RP04/Ag-Outback cross.

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.

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.

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.

Genome-wide Association Study Identifies New Loci for Resistance to Leptosphaeria maculans in Canola

Key message "We identified both quantitative and quantitative resistance loci to Leptosphaeria maculans, a fungal pathogen, causing blackleg disease in canola. Several genome-wide significant associations were detected at known and new loci for blackleg resistance. We further validated statistically significant associations in four genetic mapping populations, demonstrating that GWAS marker loci are indeed associated with resistance to L. maculans. One of the novel loci identified for the first time, Rlm12, conveys adult plant resistance in canola." Blackleg, caused by Leptosphaeria maculans, is a significant disease which affects the sustainable production of canola (Brassica napus). This study reports a genome-wide association study based on 18,804 polymorphic SNPs to identify loci associated with qualitative and quantitative resistance to L. maculans. Genomic regions delimited with 694 significant SNP markers, that are associated with resistance evaluated using 12 single spore isolates and pathotypes from four canola stubble were identified. Several significant associations were detected at known disease resistance loci including in the vicinity of recently cloned Rlm2/LepR3 genes, and at new loci on chromosomes A01/C01, A02/C02, A03/C03, A05/C05, A06, A08, and A09. In addition, we validated statistically significant associations on A01, A07, and A10 in four genetic mapping populations, demonstrating that GWAS marker loci are indeed associated with resistance to L. maculans. One of the novel loci identified for the first time, Rlm12, conveys adult plant resistance and mapped within 13.2 kb from Arabidopsis R gene of TIR-NBS class. We showed that resistance loci are located in the vicinity of R genes of Arabidopsis thaliana and Brassica napus on the sequenced genome of B. napus cv. Darmor-bzh. Significantly associated SNP markers provide a valuable tool to enrich germplasm for favorable alleles in order to improve the level of resistance to L. maculans in canola.

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.

Genetic discontinuities and disequilibria in recently established populations of the plant pathogenic fungus Mycosphaerella fijiensis

Dispersal processes of fungal plant pathogens can be inferred from analysis of spatial genetic structures resulting from recent range expansion. The relative importance of long-distance dispersal (LDD) events vs. gradual dispersal in shaping population structures depends on the geographical scale considered. The fungus Mycosphaerella fijiensis, pathogenic on banana, is an example of a recent worldwide epidemic. Founder effects in this species were detected at both global and continental scale, suggesting stochastic spread of the disease through LDD events. In this study, we analysed the structure of M. fijiensis populations in two recently (∼1979-1980) colonized areas in Costa Rica and Cameroon. Isolates collected in 10-15 sites distributed along a ∼250- to 300- km-long transect in each country were analysed using 19 microsatellite markers. We detected low-to-moderate genetic differentiation among populations in both countries and isolation by distance in Cameroon. Combined with historical data, these observations suggest continuous range expansion at the scale of banana-production area through gradual dispersal of spores. However, both countries displayed specific additional signatures of colonization: a sharp discontinuity in gene frequencies was observed along the Cameroon transect, while the Costa Rican populations seemed not yet to have reached genetic equilibrium. These differences in the genetic characteristics of M. fijiensis populations in two recently colonized areas are discussed in the light of historical data on disease spread and ecological data on landscape features.