Sources of genomic diversity in the self-fertile plant pathogen, Sclerotinia sclerotiorum, and consequences for resistance breeding.

The ascomycete, Sclerotinia sclerotiorum, has a broad host range and causes yield loss in dicotyledonous crops world wide. Genomic diversity was determined in a population of 127 isolates obtained from individual canola (Brassica napus) fields in western Canada. Genotyping with 39 simple sequence repeat (SSR) markers revealed each isolate was a unique haplotype. Analysis of molecular variance showed 97% was due to isolate and 3% due to geographical location. Testing of mycelium compatibility among 133 isolates identified clones of mutually compatible isolates with 86-95% similar SSR haplotype, whereas incompatible isolates were highly diverse. In the Province of Manitoba, 61% of isolates were compatible forming clones and stings of pairwise compatible isolates not described before. In contrast, only 35% of isolates were compatible in Alberta without forming clones and strings, while 39% were compatible in Saskatchewan with a single clone, but no strings. These difference can be explained by wetter growing seasons and more susceptible crop species in Manitoba favouring frequent mycelium interaction and more life cycles over time, which might also explain similar differences observed in other geographical areas and host crops. Analysis of linkage disequilibrium rejected random recombination, consistent with a self-fertile fungus, restricted outcrossing due to mycelium incompatibility, and only a single annual opportunity for genomic recombination during meiosis in the ascospore stage between non-sister chromatids in the rare event nuclei from different isolates come together. More probable sources of genomic diversity is slippage during DNA replication and point mutation affecting single nucleotides that accumulate and likely increase mycelium incompatibility in a population over time. A phylogenetic tree based on SSR haplotype grouped isolates into 17 sub-populations. Aggressiveness was tested by inoculating one isolate from each sub-population onto B. napus lines with quantitative resistance. Analysis of variance was significant for isolate, line, and isolate by line interaction. These isolates represent the genomic and pathogenic diversity in western Canada, and are suitable for resistance screening in canola breeding programs.

Modeling first order additive × additive epistasis improves accuracy of genomic prediction for sclerotinia stem rot resistance in canola.

The fungus Sclerotinia sclerotiorum infects hundreds of plant species including many crops. Resistance to this pathogen in canola (Brassica napus L. subsp. napus) is controlled by numerous quantitative trait loci (QTL). For such polygenic traits, genomic prediction may be useful for breeding as it can capture many QTL at once while also considering nonadditive genetic effects. Here, we test application of common regression models to genomic prediction of S. sclerotiorum resistance in canola in a diverse panel of 218 plants genotyped at 24,634 loci. Disease resistance was scored by infection with an aggressive isolate and monitoring over 3 wk. We found that including first-order additive × additive epistasis in linear mixed models (LMMs) improved accuracy of breeding value estimation between 3 and 40%, depending on method of assessment, and correlation between phenotypes and predicted total genetic values by 14%. Bayesian models performed similarly to or worse than genomic relationship matrix-based models for estimating breeding values or overall phenotypes from genetic values. Bayesian ridge regression, which is most similar to the genomic relationship matrix-based approach in the amount of shrinkage it applies to marker effects, was the most accurate of this family of models. This confirms several studies indicating the highly polygenic nature of sclerotinia stem rot resistance. Overall, our results highlight the use of simple epistasis terms for prediction of breeding values and total genetic values for a complex disease resistance phenotype in canola.

A whole genome scan of SNP data suggests a lack of abundant hard selective sweeps in the genome of the broad host range plant pathogenic fungus Sclerotinia sclerotiorum.

The pathogenic fungus Sclerotinia sclerotiorum infects over 600 species of plant. It is present in numerous environments throughout the world and causes significant damage to many agricultural crops. Fragmentation and lack of gene flow between populations may lead to population sub-structure. Within discrete recombining populations, positive selection may lead to a 'selective sweep'. This is characterised by an increase in frequency of a favourable allele leading to reduction in genotypic diversity in a localised genomic region due to the phenomenon of genetic hitchhiking. We aimed to assess whether isolates of S. sclerotiorum from around the world formed genotypic clusters associated with geographical origin and to determine whether signatures of population-specific positive selection could be detected. To do this, we sequenced the genomes of 25 isolates of S. sclerotiorum collected from four different continents-Australia, Africa (north and south), Europe and North America (Canada and the northen United States) and conducted SNP based analyses of population structure and selective sweeps. Among the 25 isolates, there was evidence for two major population clusters. One of these consisted of 11 isolates from Canada, the USA and France (population 1), and the other consisted of nine isolates from Australia and one from Morocco (population 2). The rest of the isolates were genotypic outliers. We found that there was evidence of outcrossing in these two populations based on linkage disequilibrium decay. However, only a single candidate selective sweep was observed, and it was present in population 2. This sweep was close to a Major Facilitator Superfamily transporter gene, and we speculate that this gene may have a role in nutrient uptake from the host. The low abundance of selective sweeps in the S. sclerotiorum genome contrasts the numerous examples in the genomes of other fungal pathogens. This may be a result of its slow rate of evolution and low effective recombination rate due to self-fertilisation and vegetative reproduction.

Microsatellite markers used for genome-wide association mapping of partial resistance to Sclerotinia sclerotiorum in a world collection of Brassica napus.

The fungal pathogen Sclerotinia sclerotiorum causes stem rot of oilseed rape (Brassica napus) worldwide. In preparation for genome-wide association mapping (GWAM) of sclerotinia resistance in B. napus, 152 accessions from diverse geographical regions were screened with a single Canadian isolate, #321. Plants were inoculated by attaching mycelium plugs to the main stem at full flower. Lesion lengths measured 7, 14 and 21 days after inoculation were used to calculate the area under the disease progress curve (AUDPC). Depth of penetration was noted and used to calculate percent soft and collapsed lesions (% s + c). The two disease traits were highly correlated (r = 0.93). Partially resistant accessions (AUDPC <7 and % s + c <2) were identified primarily from South Korea and Japan with a few from Pakistan, China and Europe. Genotyping of accessions with 84 simple sequence repeat markers provided 690 polymorphic loci for GWAM. The general linear model in TASSEL best fitted the data when adjusted for population structure (STRUCTURE), GLM + Q. After correction for positive false discovery rate, 34 loci were significantly associated with both disease traits of which 21 alleles contributed to resistance, while the remaining enhanced susceptibility. The phenotypic variation explained by the loci ranged from 6 to 25 %. Five loci mapped to published quantitative trait loci conferring sclerotinia resistance in Chinese lines.

IGS Minisatellites Useful for Race Differentiation in Colletotrichum lentis and a Likely Site of Small RNA Synthesis Affecting Pathogenicity.

Colletotrichum lentis is a fungal pathogen of lentil in Canada but rarely reported elsewhere. Two races, Ct0 and Ct1, have been identified using differential lines. Our objective was to develop a PCR-probe differentiating these races. Sequences of the translation elongation factor 1α (tef1α), RNA polymerase II subunit B2 (rpb2), ATP citrate lyase subunit A (acla), and internal transcribed spacer (ITS) regions were monomorphic, while the intergenic spacer (IGS) region showed length polymorphisms at two minisatellites of 23 and 39 nucleotides (nt). A PCR-probe (39F/R) amplifying the 39 nt minisatellite was developed which subsequently revealed 1-5 minisatellites with 1-12 repeats in C. lentis. The probe differentiated race Ct1 isolates having 7, 9 or 7+9 repeats from race Ct0 having primarily 2 or 4 repeats, occasionally 5, 6, or 8, but never 7 or 9 repeats. These isolates were collected between 1991 and 1999. In a 2012 survey isolates with 2 and 4 repeats increased from 34% to 67%, while isolated with 7 or 9 repeats decreased from 40 to 4%, likely because Ct1 resistant lentil varieties had been grown. The 39 nt repeat was identified in C. gloeosporioides, C. trifolii, Ascochyta lentis, Sclerotinia sclerotiorum and Botrytis cinerea. Thus, the 39F/R PCR probe is not species specific, but can differentiate isolates based on repeat number. The 23 nt minisatellite in C. lentis exists as three length variants with ten sequence variations differentiating race Ct0 having 14 or 19 repeats from race Ct1 having 17 repeats, except for one isolate. RNA-translation of 23 nt repeats forms hairpins and has the appropriate length to suggest that IGS could be a site of small RNA synthesis, a hypothesis that warrants further investigation. Small RNA from fungal plant pathogens able to silence genes either in the host or pathogen thereby aiding infection have been reported.

Expression and regulation of Sclerotinia sclerotiorum necrosis and ethylene-inducing peptides (NEPs).

Successful host colonization by necrotrophic plant pathogens requires the induction of plant cell death to provide the nutrients needed for infection establishment and progression. We have cloned two genes encoding necrosis and ethylene-inducing peptides from Sclerotinia sclerotiorum, which we named SsNep1 and SsNep2. The peptides encoded by these genes induce necrosis when expressed transiently in tobacco leaves. SsNep1 is expressed at a very low level relative to SsNep2 during infection. The expression of SsNep2 was induced by contact with solid surfaces and occurred in both the necrotic zone and at the leading margin of the infection. SsNep2 expression was dependent on calcium and cyclic adenosine monophosphate signalling, as compounds affecting these pathways reduced or abolished SsNep2 expression coincident with a partial or total loss of virulence.

Patterns of differential gene expression in Brassica napus cultivars infected with Sclerotinia sclerotiorum.

SUMMARY The fungal pathogen Sclerotinia sclerotiorum infects a broad range of dicotyledonous plant species and causes stem rot in Brassica napus. To elucidate the mechanisms underlying the defence response, the patterns of gene expression in the partially resistant B. napus cultivar ZhongYou 821 (ZY821) and the susceptible cultivar Westar were studied using a B. napus oligonucleotide microarray. Although maximum differential gene expression was observed at 48 h post-inoculation (hpi) in both cultivars, increased transcript levels were detected in cv. ZY821 at the earlier stages of infection (6-12 hpi) for many genes, including those encoding defence-associated proteins, such as chitinases, glucanases, osmotins and lectins, as well as genes encoding transcription factors belonging to the zinc finger, WRKY, APETALA2 (AP2) and MYB classes. In both cultivars, genes encoding enzymes involved in jasmonic acid, ethylene and auxin synthesis were induced, as were those for gibberellin degradation. In addition, changes in the expression of genes encoding enzymes involved in carbohydrate and energy metabolism appeared to be directed towards shuttling carbon reserves to the tricarboxylic acid cycle and generating reactive oxygen species. Transcripts from genes encoding enzymes involved in glucosinolate and phenylpropanoid biosynthesis were highly elevated in both cultivars, suggesting that secondary metabolites are also components of the response to S. sclerotiorum in B. napus.

Differential expression of duplicated peroxidase genes in the allotetraploid Brassica napus.

Gene redundancy due to polyploidization provides a selective advantage for plant adaptation. We examined the expression patterns of two peroxidase genes (BnPOX1 and BnPOX2) in the natural allotetraploid Brassica napus and the model diploid progenitors Brassica rapa (Br) and Brassica oleracea (Bo) in response to the fungal pathogen Sclerotinia sclerotiorum. We demonstrated that the Bo homeolog of BnPOX1 was up-regulated after infection, while both BnPOX2 homeologs were down-regulated. A bias toward reciprocal expression of the homeologs of BnPOX1 in different organs in the natural allotetraploid of B. napus was also observed. These results suggest that subfunctionalization of the duplicated BnPOX genes after B. napus polyploidization as well as subneofunctionalization of the homeologs in response to this specific biotic stress has occurred. Retention of expression patterns in the diploid progenitors and the natural allotetraploid in some organs indicates that the function of peroxidase genes has been conserved during evolution.

Brassica napus possesses an expanded set of polygalacturonase inhibitor protein genes that are differentially regulated in response to Sclerotinia sclerotiorum infection, wounding and defense hormone treatment.

Most plants encode a limited set of polygalacturonase inhibitor (PGIP) genes that may be involved in aspects of plant development, but more importantly in the inactivation of polygalacturonases (PG) secreted by pathogens. Previously, we characterized two Brassica napus PGIP genes, BnPgip1 and BnPgip2, which were differentially expressed in response to pathogen infection and wounding. Here we report that the B. napus genome encodes a set of at least 16 PGIP genes that are similar to BnPgip1 or BnPgip2. This is the largest Pgip gene family reported to date. Comparison of the BnPGIPs revealed several sites within the xxLxLxx region of leucine rich repeats that form beta-sheets along the interacting face of the PGIP that are hypervariable and represent good candidates for generating PGIP diversity. Characterization of the regulatory regions and RT-PCR studies with gene-specific primers revealed that individual genes were differentially responsive to pathogen infection, mechanical wounding and signaling molecules. Many of the BnPgip genes responded to infection by the necrotic pathogen, Sclerotinia sclerotiorum; however, these genes were also induced either by jasmonic acid, wounding and salicylic acid or some combination thereof. The large number of PGIPs and the differential manner in which they are regulated likely ensures that B. napus can respond to attack from a broad spectrum of pathogens and pests.

Interaction of Sclerotinia sclerotiorum with a resistant Brassica napus cultivar: expressed sequence tag analysis identifies genes associated with fungal pathogenesis.

Sclerotinia sclerotiorum is a ubiquitous necrotrophic fungal pathogen capable of infecting a wide range of plants. To identify genes involved in fungal development and pathogenesis we generated 2232 expressed sequence tags (ESTs) from two cDNA libraries constructed using either mycelia grown in pectin medium or tissues from infected Brassica napus stems. A total of 774 individual fungal genes were identified of which 39 were represented only among the infected plant EST collection. Annotation of 534 unigenes was possible following the categories applied to Saccharomyces cerevisiae and the Universal Gene Ontology scheme. cDNAs were identified that encoded potential pathogenicity factors including four endopolygalacturonases, two exopolygalacturonases, and several metabolite transporters. The potential role of these genes, as well as those encoding signal transduction factors, in the infection process is discussed.

Interaction of Sclerotinia sclerotiorum with Brassica napus: cloning and characterization of endo- and exo-polygalacturonases expressed during saprophytic and parasitic modes.

Five major and several minor PG isoenzymes were identified in a Sclerotinia sclerotiorum isolate from Brassica napus by isoelectric focusing and pectin gel overlays. Using a combination of degenerate PCR and expressed sequence tags (ESTs) four endo-polygalacturonase (PG) genes, designated as sspg1d, sspg3, sspg5, and sspg6, and two exo-PG genes, ssxpg1 and ssxpg2, were identified. SSPG1d is a member of the PG gene family previously described by Fraissinet-Tachet et al. [Curr. Genet. 29 (1995) 96]. The mature SSPG1d is a neutral PG, whereas fully processed SSPG3, SSPG5, and SSPG6 are acidic enzymes. Under saprophytic growth conditions, sspg1d, sspg3, sspg5, and ssxpg1 expression was induced by pectin and galacturonic acid and subject to catabolite repression by glucose. Conditions could not be identified under which sspg6 or ssxpg2 were expressed well. Transfer of mycelia from liquid media to solid substrates induced expression of sspg1d suggesting that it may also be regulated by thigmotrophic interactions. Under pathogenic conditions, sspg1d was highly expressed during infection. sspg3 was also expressed during infection, albeit at lower levels than sspg1d, whereas sspg5, sspg6, and ssxpg1 were expressed only weakly.