Novel gene loci associated with susceptibility or cryptic quantitative resistance to Pyrenopeziza brassicae in Brassica napus

doi:10.21203/rs.3.rs-1327286/v1

PDF

Research Article

Novel gene loci associated with susceptibility or cryptic quantitative resistance to Pyrenopeziza brassicae in Brassica napus

https://doi.org/10.21203/rs.3.rs-1327286/v1

This work is licensed under a CC BY 4.0 License

You are reading this latest preprint version

Quantitative disease resistance (QDR) can provide durable control of pathogens in crops in contrast to resistance (R) gene-mediated resistance which can break down due to pathogen evolution. QDR is therefore a desirable trait in crop improvement but little is known about the causative genes and so it is difficult to incorporate into breeding programmes. Light leaf spot, caused by Pyrenopeziza brassicae, is an important disease of oilseed rape (canola, Brassica napus). To identify new QDR gene loci, we used a high-throughput screening pathosystem with P. brassicae on 196 lines of B. napus combined with an association transcriptomics platform. We show that all resistance against P. brassicae was associated with QDR and not R gene-mediated. We used genome-wide association analysis with an improved B. napus population structure to reveal four gene loci significantly (P = 0.0001) associated with QDR in regions showing linkage disequilibrium. On chromosome A09, enhanced resistance was associated with heterozygosity for a cytochrome P450 gene co-localising with a previously-described locus for seed glucosinolate content. In addition, eight significant gene expression markers with a false discovery rate of 0.001 were associated with QDR against P. brassicae. For seven of these, expression was positively correlated with resistance whereas for one, a HXXXD-type acyl-transferase, negative correlation indicated a potential susceptibility gene. The study identifies novel QDR loci for susceptibility and resistance, including novel cryptic QDR genes associated with heterozygosity, that will inform future crop improvement.

Apoplast

Cylindrosporium concentricum

host-pathogen interaction

quantitative trait locus

rapeseed

subcuticular niche

Quantitative resistance (QR) controls the association of the light leaf spot pathogen with Brassica napus; four QR loci that were in linkage disequilibrium and eight gene expression markers were identified.

Quantitative disease resistance (QDR) is predominant in natural plant populations and provides robust and durable protection from pathogens in ecosystems (Delplace et al. 2020). Although QDR could potentially provide more durable resistance in crop plants, it is not routinely selected for in breeding since the underlying gene loci are poorly characterised and contribute partial but additive resistance which is difficult to track in breeding programmes (Nelson et al. 2018). Understanding the molecular basis of QDR will provide novel opportunities for introducing durable disease resistance into crops and reducing the use of pesticides. Our investigation uses a pathosystem which, combined with an association genetics platform, establishes a novel approach to identify candidate QDR genes for crop improvement in Brassica napus.

Plants have evolved layers of immunity for defence against various pathogens with different modes of infection and life styles. Although the layers of plant immunity are not distinct and may overlap to some extent (Ngou et al. 2021; Yuan et al. 2021), the first is considered to be the rapid detection of microbial patterns usually described as pattern-triggered immunity (PTI). Adapted pathogens secrete effectors that suppress PTI as they colonise the host plant and cause disease (Jones and Dangl 2006). Effectors can be recognised by resistance proteins, typically nucleotide-binding, leucine-rich repeat (NLR) immune receptors, which provide the second layer of immunity to protect the plant. When operating against cell-penetrating pathogens, this effector-triggered immunity (ETI) is the principal mechanism in qualitative, or R gene-mediated, resistance and provides a powerful defence response (Jones & Dangl 2006). However, operating against extra-cellular pathogens, R gene-mediated resistance, termed effector-triggered defence (ETD) acts more slowly and is not an immune but a resistance response (Stotz et al. 2014). Pathogen populations can mutate or lose effectors so that they are no longer recognised, leading to resistance breakdown in crops. This contrasts with QDR, which is predicted to provide durable protection against pathogens. Although less well-studied, mechanisms of QDR can include enhanced cell wall thickening, modified signalling processes and enhanced secondary metabolism (Cowger and Brown 2019). Interestingly, studies with the apoplastic (extracellular) fungal pathogen Leptosphaeria maculans have shown that ETD may also contribute to QDR in B. napus, reinforcing the concept of overlap between layers of defence in plants (Jiquel et al. 2021).

After soybean, oilseed rape (OSR; Brassica napus) is economically the second most important vegetable oil crop in the world. Amongst the biotic threats that challenge OSR production, light leaf spot (LLS), caused by the apoplastic ascomycete Pyrenopeziza brassicae (anamorph Cylindrosporium concentricum), ranks in the top 10 most damaging diseases on the crop in Europe (Zheng et al. 2020). This pathogen has a widespread geographic distribution, occurring in the United Kingdom (UK) and, increasingly, continental Europe, the Pacific Northwest of the United States, Asia (Japan, Philippines) and New Zealand (Carmody et al. 2019; Karandeni Dewage et al. 2018). P. brassicae epidemics are initiated by wind-dispersed ascospores. Acervuli then produce asexual conidia in infected plant parts; these conidia are rain splash-dispersed to establish the polycyclic stage of the LLS disease epidemic.

LLS has become the most damaging disease of OSR in the UK. The disease accounts for up to £160 million yield loss annually in England, despite expenditure of £20M on fungicides. The severity of the disease is greater in Scotland than in England (Karandeni Dewage et al. 2018). However, LLS has progressively become a greater problem in parts of the UK other than Scotland over more than a decade and may now be considered a national emergency. Simultaneously, P. brassicae has become a problem on brassica vegetables in the UK and other places, causing market losses due to surface blemishes on crops like Brussels sprouts. Besides B. napus and B. oleracea, P. brassicae has been observed on mustard rape (B. juncea) and B. rapa (Carmody et al. 2019).

One of the limitations of studying QDR in crop species is that resistance of individual QDR loci is only partial and can be masked by R-gene mediated resistance. P. brassicae on B. napus is an ideal pathosystem to study QDR because there are no known R genes effective against the pathogen. Moreover, one study has indicated the presence of QDR for LLS. Six environmentally stable quantitative trait loci (QTL) for resistance against P. brassicae were mapped in a doubled haploid (DH) population of B. napus derived from a cross between moderately resistant Darmor-bzh and susceptible Yudal cultivars (Pilet et al. 1998). Four new QTL for resistance against P. brassicae were identified on linkage groups C01, C06 and C09 (Karandeni Dewage et al. 2022). Two major loci for resistance against P. brassicae were identified (Bradburne et al. 1999) and one was mapped to the bottom of chromosome A1 of B. napus (Boys et al. 2012). This locus could be a QDR locus since it was associated with substantial decrease in LLS hyphal growth rather than death of the pathogen typical of R gene-mediated resistance (Boys et al. 2012; Stotz et al. 2014). However, apart from these examples, little is known about QDR for resistance against P. brassicae in B. napus.

P. brassicae and Rhynchosporium commune are closely related discomycetes (Goodwin 2002; Penselin et al. 2016) that occupy a subcuticular apoplastic niche in their respective brassica and barley hosts (Stotz et al. 2014). It is therefore relevant to compare resistance mechanisms that operate against these related pathogens. The major resistance locus Rrs1, containing wall-associated kinases, controls resistance against Rhynchosporium commune in barley (Looseley et al. 2020). Rrs1 is the only R gene against R. commune that has a corresponding Avr gene, NIP1, that encodes a necrosis-inducing protein (Rohe et al. 1995). Additionally, multiple quantitative resistance loci were found to be involved in resistance of barley against R. commune (Buttner et al. 2020).

The development of association genetics and transcriptomics has enabled the identification of pathogen resistance loci that was not possible in bi-parental mapping populations (Bartoli and Roux 2017). The approaches take advantage of recombination events that have accumulated in natural populations to identify genetic polymorphisms associated with phenotypes of interest. The method requires genomic or transcriptomic sequences from a diversity collection and phenotypic data for the trait of interest. The approach has been used to identify new resistance loci operating against Sclerotinia sclerotiorum stem rot and potential R genes for resistance in B. napus against Plasmodiophora brassicae (clubroot) (Hejna et al. 2019; Wu et al. 2016). However, no such studies have been done to identify QDR against Pyrenopeziza brassicae.

The aim of this study was to characterise genomic regions associated with QDR against P. brassicae in OSR under glasshouse conditions. To achieve this, we performed association genetic analysis with a B. napus diversity set developed through the OREGIN initiative (https://www.herts.ac.uk/oregin) combined with genotype and expression data (Havlickova et al. 2018) and phenotypic measurement of P. brassicae infection that we developed for the study. This study provides new insights into QDR mechanisms, and supports breeding efforts to generate durable disease-resistant crops.

Glasshouse growth conditions

Glasshouses at the Bayfordbury campus, University of Hertfordshire, and Rothamsted Research were utilised for experiments scored for partial resistance against P. brassicae with baseline temperatures set to 16°C during the day and 14°C during the night (see diurnal cycles below). Actual temperatures recorded were outside this baseline range, however, due to fluctuations in temperature and light during the experimental period (Supporting Information Table S1). At Bayfordbury, supplemental lighting was used for 12 h per day using sodium high pressure lamps (Sylvania SHP-TS 400 W GroLux), which automatically switched on once natural daylight decreased to < 115 µmol m⁻² s⁻¹. Supplemental lighting was used for 12 h and 14.5 h at Rothamsted Research for experiments 8 and experiments 9 and 10, respectively (Supporting Information Table S1).

Isolation of P. brassicae populations

Two populations of P. brassicae were used. The first population was amplified having obtained infected leaves from KWS Saat SE & Co. KGaA (Einbeck, Germany) from field experiments on the island in Fehmarn, Germany (54.4701, 11.1329) in 2016. The second population was obtained from infected leaves of KWS-grown B. napus genotypes, including reference cultivars Cuillin and Express, at a Rothamsted Research field site (51.813125, -0.382005) in 2019. Populations obtained in 2016 and 2019 were used for seven and three glasshouse experiments at Bayfordbury and Rothamsted Research, respectively (Supporting Information Table S2).

Leaves of B. napus with the greatest amount of P. brassicae sporulation already present were sampled from the field experiments and placed into polyethylene bags for transportation back to the University of Hertfordshire. The spores were dislodged from the leaves by pipetting sterile distilled water over the leaf and the run-off was collected in a glass beaker. This spore suspension was filtered using Miracloth (Merck Millipore, Watford, U.K.) and the spore concentration determined using a haemocytometer slide under a stereo-microscope (GX Microscopes, XTC-3A1). The concentration was adjusted to 10⁵ spores ml⁻¹ and dispensed into 50 ml aliquots for use as inoculum in the experiments or stored at -20^oC.

Quantification of P. brassicae sporulation on leaves of B. napus accessions

Each of 10 experiments tested 24 accessions, each with five inoculated replicates with the exception of experiments 1 and 3, which each tested only 23 accessions (Supporting Information Table S2). In total there were 195 accessions in these experiments. The replicates were arranged in a randomised alpha block design, so that each block included one replicate plant from each of the 24 accessions. All of the lines screened could not be grown at the same time; therefore seven experiments were done at Bayfordbury from 2016-2018 with the cultivars Tapidor, Imola (resistant), Bristol (susceptible) and Temple included in each experiment to act as references for normalising symptom severity between experiments. In 2019, another three experiments were done at Rothamsted Research with the susceptible one of the replicated reference cultivars being Cabriolet instead of Bristol.

Seeds were stratified for 2 days at 4^oC prior to germination on damp filter paper in Petri dishes in the dark. The successfully germinated seeds were then planted into 8 x 5 cell seed trays with a 50:50 mixture of John Innes No.3 compost and multipurpose compost (Miracle Gro). The seedlings were left to grow in the glasshouse for approximately 2 weeks until the first and second true leaves had emerged. Established plants were transplanted into 9 cm diameter individual pots using the same potting mix.

After a further 2 weeks, a total of 120 plants at growth stage 1,5 (Sylvester-Bradley 1985) were sprayed with 100 ml of P. brassicae spore suspension, at a concentration of 10⁵ spores ml⁻¹ water containing 0.05% Silwet (Spiess-Urania Chemicals GmbH, Hamburg, Germany) as a wetting agent. Subsequently, plants were incubated for 48 h under a clear polyethylene sheet over a plastic frame, 2 x 1 x 0.5 m in dimension. The plants were grown for another 3 weeks before the whole plant was sampled by cutting at the base of the stem just above the soil surface, wrapped in a paper towel, and placed into a labelled polyethylene freezer bag to preserve the whole plant. The inside of the bag and the paper towel were sprayed with sterile distilled water to create humidity and then sampled plants were incubated at 4°C in a cold room. After 5 and 10 days of incubation, the sampled plants were assessed visually for the presence of P. brassicae sporulation using a modified LLS 1-6 severity scale (Karandeni Dewage et al. 2021). Specifically, score = 1 (no sporulation), score = 2 (< 10% leaf area with sporulation), score = 3 (10-25% leaf area with sporulation), score = 4 (25-50% leaf area with sporulation), score = 5 (50-75% leaf area with sporulation) and score = 6 (75-100% leaf area with sporulation). Each of three to six leaves per plant was given a score and then these scores were averaged to give a disease severity rating.

Statistical analysis of phenotypic traits

The quantitative resistance work consisted of 10 glasshouse experiments. A total of 195 accessions were scored for P. brassicae sporulation. Reference genotypes replicated in each of the ten experiments were used to determine normality of the data with the Shapiro-Wilk test and lack of effect of experiment for the entire study. Phenotypic data were analysed in R using a package for non-linear mixed-effects models (nlme). The function for the model had the following structure: model=lme(Score~Genotype, data=LLS, random =~1|Experiment/Genotype). Adjusted means were generated and ordered according to the scores of the genotypes. A Wilcoxon signed-rank test was used for pair-wise comparisons of the adjusted means of genotypes. A significance threshold of P < 0.05 was used to assess pair-wise differences in disease score. Data were visualised in R using ordered means dotchart and hist functions.

Association transcriptomics

The phenotyping datasets of the diversity set for LLS resistance scores were analysed using an association transcriptomics pipeline based on programs used to map traits in B. napus previously with minor modifications (Harper et al. 2012; Wells et al. 2013). Genotype and expression level datasets used were published (Havlickova et al. 2018) and available from York Knowledgebase (http://yorknowledgebase.info). This dataset was reduced to include only the accessions within this study.

Gene expression marker (GEM) associations were determined by linear regression using Reads Per Kilobase of transcript, per Million mapped reads (RPKM) to predict a quantitative outcome of the trait value. All markers with an average expression less than 0.5 RPKM were removed before analysis.

An updated population structure was calculated for the accessions used within this study using a Bayesian clustering approach. A Markov Chain Monte Carlo (MCMC) algorithm was implemented in the population-genetic software STRUCTURE V2.3.1. One of the requirements of STRUCTURE is unlinked markers; therefore, the single nucleotide polymorphism (SNP) file was adjusted before the analysis using the following criteria: SNPs were required to be biallelic, with a minor allele frequency (MAF) > 0.05 and a minimum distance of 500 kb between markers. Markers within 100 kb of the centromeres, based on published findings (Cheng et al. 2013; Mason et al. 2016), were excluded. STRUCTURE was run using the admixture model with uncorrelated allele frequencies with a burn-in period of 100,000 iterations and MCMC analyses of 100,000 permutations. The accessions were not assigned to a given population. Ten iterations were run for each value of K, the number of subpopulations estimated to make up the total population. STRUCTURE HARVESTER (Earl and vonHoldt 2012) was used to determine the optimal K value, by generating a series of ΔK values, which represent the mean likelihood of K divided by the standard deviation of K, for the population. To further investigate population clusters TASSEL v5 was used to construct a phylogenetic tree, using the Neighbour Joining method and all SNPs with MAF > 0.05. A cluster matching and permutation program (CLUMPP) was used from STRUCTURE HARVESTER for K = 6 to generate the Q matrix input (Jakobsson and Rosenberg 2007).

Genome-wide association (GWA) mapping was performed using TASSEL v5 using SNP markers with an allele frequency > 0.05. Analysis was conducted using generalised linear models (GLM) and mixed linear models (MLM) to determine the optimal model. Kinship data were calculated using TASSEL’s ‘centered IBS’ method. Optimum compression level and P3D variance component estimation were used as MLM options. The false discovery rate (FDR) was determined using the Shiny implementation of the q-value R package (Storey et al. 2020). GEM and GWA results were visualised using R (https://github.com/BRAVO-research-project/pyrenopeziza-resistance). The most significantly associated SNP markers were selected for further analysis, including distribution within the population and allelic effect.

The level of linkage disequilibrium (LD) varies between and across chromosomes depending on the position and level of selection. To determine the specific level of LD at each locus, the mean pair-wise r² for all markers on a chromosome to each of the 11 significantly associated markers was calculated within TASSEL v5 using the site by all analysis option. Markers were considered in LD when r² > 0.15.

LLS scoring system

Percentage of leaf area with P. brassicae sporulation observed on infected leaves was the most consistent measurement of disease severity across all accessions (Karandeni Dewage et al. 2021) and was therefore used as the basis for the scoring system. Disease severity was scored on a scale of 1 to 6, with a score of 1 for no sporulation and 6 for the most sporulation (Fig. 1A). The image of part of an infected cv. Tapidor leaf (Fig. 1A, score 2) shows why P. brassicae is referred as to C. concentricum in its imperfect stage; this pathogen produces concentric rings of acervuli on its hosts. Patchy sporulation was observed on leaf laminas of ‘couve-nabiça’ and cv. Capitol with scores of 3 to 4, respectively (Fig. 1A). The entire leaf laminas of cv. Musette and cv. Daichousen were covered with acervuli and scored 5 and 6, respectively. Whilst the resistance response black necrotic flecking was sometimes observed during the screening process, it was not present on all of the accessions tested, and so this phenotype was not considered suitable for scoring QDR.

Variation in sporulation of P. brassicae on leaves of diverse B. napus accessions

A diversity set of 195 accessions was tested in glasshouse experiments for the amount of pathogen sporulation after spray inoculation of B. napus seedlings with local P. brassicae populations. The distribution of disease scores showed wide variation in P. brassicae sporulation among diverse B. napus accessions with an approximately normal distribution of this trait (Fig. 1B).

The glasshouse screen consisted of 10 independent experiments, each with 23 or 24 B. napus accessions with four reference cultivars per experiment (Supporting Information Table S2). The first seven experiments were done at the Bayfordbury campus of the University of Hertfordshire and the last three experiments at Rothamsted Research. Instead of cv. Bristol, cv. Cabriolet was used as a susceptible reference cultivar for the experiments at Rothamsted Research. All other reference cultivars were identical between both sites. Irrespective of the location, there were significant differences between reference cultivars in disease score. No significant effects of experiment and experiment-by-cultivar interaction on disease score were observed (Supporting Information Notes S1). Combination of all 10 experiments and comparison of the shared reference cultivars Imola, Tapidor and Temple resulted in significant effects of both cultivar and experiment, but no significant experiment-by-cultivar interactions on disease score (Supporting Information Table S2). The disease scores of all three reference cultivars were less at Rothamsted Research than at Bayfordbury, with cv. Temple differing the most between these two environments (Supporting Information Figure S1). Restricting the analysis to cv. Imola and cv. Tapidor eliminated the significant effect of experiment. Differences in disease scores between Rothamsted Research and Bayfordbury could have resulted from different environmental conditions, pathogen inoculum and assessors at the two sites.

Non-linear mixed model analysis established that B. napus cv. Cabriolet and cv. Imola scored as the most susceptible and most resistant cultivars amongst the five reference cultivars tested, respectively (Supporting Information Table S2). Intermediate scores were observed for the other three cultivars Temple, Tapidor and Bristol. While cv. Cabriolet scored as the most susceptible cultivar amongst all 195 accessions tested, 60 accessions scored less sporulation than cv. Imola; half of these 60 accessions scored significantly less than cv. Imola including all 11 accessions that supported the least P. brassicae sporulation. These data clearly show that quantitative resistance present in 30 diverse accessions resulted in less P. brassicae sporulation than in cv. Imola, which has a major QDR locus against this pathogen.

GWA mapping of quantitative resistance against P. brassicae

SNP data for the RIPR genotype dataset from the resources page of York Knowledgebase (http://yorknowledgebase.info) were used for this study. Following analysis with STRUCTURE, calculation of ΔK divided the population into two clusters (Figure 2); cluster one mainly comprising of winter OSR and fodder types and cluster two comprising of other crop types (Supporting Information Table S3). It was shown that K = 2 is a common outcome when using the ΔK method (Janes et al. 2017). ΔK frequently identifies K = 2 as the top level of hierarchical structure, and further analysis is required to determine whether more subpopulations are present. Our analysis identified a further maximum in ΔK at K = 6. This divided the population into groups comprising the different crop types; cluster one - winter and fodder; cluster two – swede; cluster three – spring OSR; clusters four and five – Siberian kale types, and six semi-winter (Chinese) OSR (Supporting Information Table S3). Subsequent phylogenetic analysis showed a delineation of these crop types (Figure 2), with the different crop types forming clear subgroups within the tree. A number of accessions did not cluster with their given crop types. However, bar diagram outputs from STRUCTURE (Supporting Information Figure S3) showed a level of admixture within these accessions. Given the evidence for population substructure beyond K = 2, a structure of K = 6 was taken forward for use in association mapping.

GWA mapping using TASSEL identified a generalised linear modelling approach as an optimal fit for the phenotypic data (Supporting Information Figure S4). Eleven significant marker associations with LLS infection score were observed at P < 0.0001 (Table 1); however, none of these were significant at the FDR 0.05. This is not unexpected as resistance against P. brassicae is a highly quantitative trait, with no known R gene loci.

The allelic effects of identified GWA maxima were determined (Table 2). The distribution of alleles contributing to resistance was not determined by crop type or phylogenetic relationship. Due to sequence similarities, cross alignment of transcriptome reads occurs between homeologous loci in the A and C genomes. This means that allelic calls can be the same in each genome or carry alternate alleles in the A and C genomes, referred to as a hemi-SNP, resulting in an ambiguity call during SNP calling. Three loci, LLSC01, LLSA09a and LLSC02, showed the strongest resistance when present as hemi-SNPs, suggesting a resistance benefit linked to carrying an alternative allele at the homeologous locus. Loci LLSA01 and LLSC04 showed that the majority of OSR in the panel carried the resistant allele; therefore these may have already been selected for during breeding. Lines carrying alternate alleles at the homeologous position were also present, suggesting some breeding lines may not be optimised for these potential resistance loci. Loci LLSA02, LLSA07, LLSA09b and LLSC08b carried both A and C resistance alleles in a small number of lines, with most lines carrying alternate alleles at the homeologous loci. For locus LLSC08b, only one line carried susceptible alleles in both sub-genomes, resulting in an elevated disease score. For loci LLSC05 and LLSC08a, resistant alleles at the two homeologous loci were not present within winter-OSR (WOSR) or in the case of LLSC08a, within the panel tested. Hemi-SNPs associated with enhanced resistance indicate currently un-described cryptic loci that could be further optimised for resistance in breeding unused loci for resistance breeding against P. brassicae.

Although the 11 GWA markers detail the genetic variation most significantly associated with resistance against P. brassicae within this analysis, transcriptome sequencing does not provide all potential genetic variants present across the panel. Genetic polymorphisms close to the causal variation will be associated on the basis of genetic linkage or LD. In the case of LD, the GWA markers may not be indicative of the causal gene. Instead, causal genes will be more or less closely linked to the GWA marker. Each GWA marker was tested for LD against all other markers on a chromosome. Four of the 11 significantly associated markers, loci LLSA02, LLSA09, LLSC02, LLSC08a, exhibited LD with r² > 0.15, thus defining the region where the causal gene is likely to be situated (Table 3, Figure 3).

Eight gene expression markers are associated with resistance or susceptibility to P. brassicae infection

The transcriptomes of the 195 accessions were used to analyse the disease scores as a function of the expression of all gene models using a general linear model with a false discovery rate of 0.001. A list of eight genes with P-values < 1.6 x 10⁻⁷ was selected for further analysis (Table 4). Linear regressions of the disease score versus the expression of each gene generated different slopes (Supporting Information Figure S2); the expression of seven genes was negatively correlated with the disease score. Thus, expression was greatest when the least sporulation was observed. These seven GEMs may therefore contribute to partial resistance against P. brassicae. Another gene, encoding an HXXXD-type acyl-transferase (Table 4), had the opposite expression pattern with an expression positively correlated with the disease score. This is therefore a candidate gene for susceptibility to P. brassicae.

QDR is the defence mechanism against P. brassicae in B. napus

Our pathosystem combined with association genetics enabled identification of candidate QDR gene loci for resistance or susceptibility to P. brassicae. GWA mapping clearly demonstrated the existence of multiple genes contributing to QDR, confirming that combinations of genes with relatively minor effects are the predominant mechanism of resistance against P. brassicae. This contrasts to resistance against L. maculans, which is controlled by both R genes (Larkan et al. 2013; Larkan et al. 2020) and quantitative resistance loci (Huang et al. 2019). We could associate enhanced resistance with hemi-SNPs at some loci, indicating the presence of previously-unknown cryptic QDR. Importantly, accessions were identified that were more resistant to P. brassicae than the well-characterised cv. Imola, which contains a single major QDR locus for resistance against this pathogen (Boys et al. 2012). Our investigation reveals novel QDR loci that are the primary mechanisms of resistance against P. brassicae. Our findings could also provide insight into mechanisms of QDR against other closely related apoplastic fungal pathosystems, including R. commune and Venturia inaequalis for durable disease control (Stotz et al. 2014).

Population structure analysis reflects history of OSR cultivation

The origin of cultivated OSR has been considered to be Europe (Lu et al. 2019). Spring cultivars that are commonly grown in North America and Australia were developed in the late 18th century and semi-winter cultivars were introduced into China in the 20th century. Moreover, distinct morphotypes and subspecies of swedes (B. napus subsp. rapifera) and kale (B. napus subsp. pabularia) were developed as root and leaf crops (An et al. 2019).

The population structure analysis using more than a thousand SNPs suggested genetic subdivisions into winter, spring, semi-winter and fodder OSR as well as swedes and Siberian kale types are consistent with other published data based on B. napus genome sequences (An et al. 2019; Lu et al. 2019). The phylogenetic comparison reported here may be biased by the B. napus accessions that were used, but it suggests WOSR to be ancestral to kales, with spring OSR, swedes and semi-spring OSR being more derived (Figure 2). The established population structure was used with associative transcriptomics to identify GWA markers and GEMs linked to quantitative resistance against P. brassicae.

GWA mapping indicates four QDR loci against P. brassicae

GWA mapping identified four loci in LD, showing multiple markers within a region associated with the infection score, and therefore more likely linked to QDR against P. brassicae than single associated SNP markers. The loci on chromosomes A02, A09, C02 and C08 were not located in regions previously identified for resistance against P. brassicae (Boys et al. 2012; Bradburne et al. 1999; Karandeni Dewage et al. 2022; Pilet et al. 1998). The LD observed on chromosome A09 coincides with a homeologous QTL for seed glucosinolate content (Qian et al. 2014). The spread of the LD on chromosomes A02, A09 and C02 was too large to identify candidate gene loci that might be responsible for QDR against P. brassicae. In contrast, the LD spread on chromosome C08 was narrow. The GWA marker Bo8g108400, corresponding to BnaC08g41550D, encodes an Cys-rich protein of unknown function. The hemi-SNP results in Gly204Cys substitution, which could be functionally significant (Perry et al. 2009). Of note, this codon is otherwise conserved in B. oleracea, B. rapa and for the homeolog BnaA09g47370D. Collectively, this could point to an important novel gene involved in host-pathogen interactions.

An alternative explanation could be that this gene was located right next to an ortholog of the Arabidopsis thaliana flavin-containing monooxygenase gene FMO_GS−OX5 (At1g12140) that converts methylthioalkyl to methylsulfinylalkyl (MS) glucosinolates. MS glucosinolates are precursors to MS isothiocyanates that are toxic to phytopathogenic fungi (Stotz et al. 2011b). Importantly, FMO_GS−OX5 has a preference for long-chain aliphatic glucosinolates (Li et al. 2008), which release isothiocyanates that are most toxic to phytopathogenic fungi like S. sclerotiorum (Stotz et al. 2011b). Notably, a mutation in the FMO_GS−OX5 gene also altered cytokinin and jasmonate levels (Garrido et al. 2020), which is of significance considering that P. brassicae is a cytokinin-producing pathogen (Ashby 1997). Although the impact of jasmonates on infections by P. brassicae has not yet been studied, jasmonates are known to influence many pathosystems (Stotz et al. 2011a; Zheng et al. 2012). The FMO_GS−OX5 gene is also part of an FMO cluster and the paralog At1g12200 that corresponds to BnaC08g41500D is induced after inoculation with S. sclerotiorum or B. cinerea (Stotz et al. 2011b). The role of this chrC08 locus in resistance against P. brassicae is therefore worth investigating.

The hemi-SNP for GWA marker Bo2g163990, encoding a nuclear localised DEAD box ATP-dependent RNA helicase, would result in a Gly565Asp substitution. However, in B. rapa and B. napus, a Ser and Asn are found in the same position, respectively. The amino acid substitutions would occur in the GUCT domain of these types of helicases. As these helicases are involved in development and abiotic stress responses (Liu et al. 2016; Perroud et al. 2021), it is not immediately obvious whether this helicase or another gene within LD are responsible for QDR at this locus.

The disclosed GWA markers will assist with the development of molecular markers for marker-assisted selection, plant breeding and crop improvement.

Putative function of GEMs in QDR

The expression of eight GEMs was correlated with resistance against or susceptibility to P. brassicae. Amongst them were three genes that were previously reported to be involved in host pathogen interactions, including a gene encoding a cinnamate-4-hydroxylase (C4H). C4H genes are induced after infection of B. napus with S. sclerotiorum or L. maculans (Becker et al. 2017; Wu et al. 2016) and infection of the liverwort Marchantia polymorpha with the oomycete Phytophthora palmivora (Carella et al. 2019). The product 4-coumarate feeds into phenylpropanoid, flavonoid and lignin biosynthesis. Mutations in the C4H gene of A. thaliana have pleiotropic developmental defects (El Houari et al. 2021). Multiple corresponding C4H genes occur in B. napus, which may have developed specialised functions, including involvement in resistance against P. brassicae.

Another gene, PLC4, is involved in resistance against pathogens. Expression of the tomato gene SlPLC4 is tightly regulated in response to the apoplastic fungal pathogen Cladosporium fulvum with a 10-fold increase in expression 7 days after infection (Vossen et al. 2010). SlPLC4 contributes to R gene-mediated resistance against Cladosporium fulvum (Vossen et al. 2010). Unlike SlPLC6, SlPLC4 is not involved in resistance against Verticillium dahliae or Pseudomonas syringae. PLC4 of A. thaliana interacts with PEN3 (Campe et al. 2016), which is involved in resistance to penetration and export of protective metabolites against microbial invasion (Lu et al. 2015; Stein et al. 2006).

Moreover, KIN10 recently emerged as a gene involved in resistance against the clubroot pathogen in A. thaliana (Chen et al. 2021). KIN10 is a central regulator in response to energy deprivation, facilitating survival under stress conditions (Baena-Gonzalez et al. 2007) including pathogen challenge (Chen et al. 2021). These findings are consistent with the correlated of KIN10 expression and resistance against P. brassicae that we observed (Supporting Information Figure S2). All of these above examples suggest that the GEMs we identified are of significance and potential breeding targets. Importantly, it also implicates other, less characterised genes, in resistance against pathogens (Table 1).

All resistance against P. brassicae found in 196 accessions of B. napus results from QDR with no evidence of R gene-mediated resistance. For genome-wide association studies, we developed an improved population structure with six phylogenetic groups based on B. napus crop types. We identified four gene loci significantly (P = 0.0001) associated with QDR and in LD. On chromosome A09, enhanced resistance was associated with a cryptic heterozygous locus for a cytochrome P450 gene co-localising with a previously-described QTL for seed glucosinolate content. The expression of seven gene expression markers were positively correlated with resistance whereas one, a HXXXD-type acyl-transferase, was negatively correlated and so is a potential susceptibility gene. The results provide new insight into QDR against P. brassicae in B. napus and can be used for marker-assisted breeding in crop improvement.

Funding

We thank KWS for their shipment of infected leaf material from their field sites. Funding by the UK Department for Environment, Food and Rural Affairs (Defra, CH0104) enabled associative transcriptomics. ERA-CAPS/BBSRC grant BB/N005007/1 (MAQBAT) provided funding for HUS, CJR, RW, BDLF, GKM, HJS and HF. HUS, CJR and RW are also supported by BBSRC grant BB/V01725X/1. CJR and RW received additional funding from BBSRC grants BBS/E/J/000PR9797 and BB/P003095/1 (BRAVO), respectively.

Competing Interests

The authors have no relevant financial or non-financial interests to disclose.

Author contributions

HF, AMA, RW, GKM, BDLF, CJR, HJS and HUS contributed to the conception and design of the study. HF, AMA, RW, GKM and HJS contributed to plant and pathogen material preparation. HF, AMA, RW, HW and HUS were involved in data collection and analysis. The first draft of the manuscript was written by HUS and all authors commented on previous versions of the manuscript. All authors read and approved the final manuscript.

Data Availability

The datasets generated during and/or analysed during the current study are available from the corresponding author on reasonable request.

An H, Qi X, Gaynor ML, Hao Y, Gebken SC, Mabry ME, McAlvay AC, Teakle GR, Conant GC, Barker MS, Fu T, Yi B, Pires JC (2019) Transcriptome and organellar sequencing highlights the complex origin and diversification of allotetraploid Brassica napus. Nat Commun 10:2878
Ashby AM (1997) A molecular view through the looking glass: The Pyrenopeziza brassicae Brassica interaction. Advances in Botanical Research 24:31-70
Baena-Gonzalez E, Rolland F, Thevelein JM, Sheen J (2007) A central integrator of transcription networks in plant stress and energy signalling. Nature 448:938-942
Bartoli C, Roux F (2017) Genome-wide association studies In plant pathosystems: Toward an ecological genomics approach. Front Plant Sci 8:763
Becker MG, Zhang X, Walker PL, Wan JC, Millar JL, Khan D, Granger MJ, Cavers JD, Chan AC, Fernando DWG, Belmonte MF (2017) Transcriptome analysis of the Brassica napus-Leptosphaeria maculans pathosystem identifies receptor, signaling and structural genes underlying plant resistance. Plant J 90:573-586
Boys EF, Roques SE, West JS, Werner CP, Kinga GJ, Dyer PS, Fitt BDL (2012) Effects of R gene-mediated resistance in Brassica napus (oilseed rape) on asexual and sexual sporulation of Pyrenopeziza brassicae (light leaf spot). Plant Pathol 61:543–554
Bradburne R, Majer D, Magrath R, Werner CP, Lewis B, Mithen R (1999) Winter oilseed rape with high levels of resistance to Pyrenopeziza brassicae derived from wild Brassica species. Plant Pathol 48:550-558
Buttner B, Draba V, Pillen K, Schweizer G, Maurer A (2020) Identification of QTLs conferring resistance to scald (Rhynchosporium commune) in the barley nested association mapping population HEB-25. BMC Genomics 21:837
Campe R, Langenbach C, Leissing F, Popescu GV, Popescu SC, Goellner K, Beckers GJ, Conrath U (2016) ABC transporter PEN3/PDR8/ABCG36 interacts with calmodulin that, like PEN3, is required for Arabidopsis nonhost resistance. New Phytol 209:294-306
Carella P, Gogleva A, Hoey DJ, Bridgen AJ, Stolze SC, Nakagami H, Schornack S (2019) Conserved biochemical defenses underpin host responses to oomycete infection in an early-divergent land plant lineage. Curr Biol 29:2282-2294 e2285
Carmody SM, King KM, Ocamb CM, Fraaije BA, West JS, du Toit LJ (2019) A phylogenetically distinct lineage of Pyrenopeziza brassicae associated with chlorotic leaf spot of Brassicaceae in North America. Plant Pathol 69:518-537
Chen W, Li Y, Yan R, Ren L, Liu F, Zeng L, Sun S, Yang H, Chen K, Xu L, Liu L, Fang X, Liu S (2021) SnRK1.1-mediated resistance of Arabidopsis thaliana to clubroot disease is inhibited by the novel Plasmodiophora brassicae effector PBZF1. Mol Plant Pathol 22:1057-1069
Cheng F, Mandakova T, Wu J, Xie Q, Lysak MA, Wang X (2013) Deciphering the diploid ancestral genome of the Mesohexaploid Brassica rapa. Plant Cell 25:1541-1554
Cowger C, Brown JKM (2019) Durability of quantitative resistance in crops: greater than we know? Annu Rev Phytopathol 57:253-277
Delplace F, Huard-Chauveau C, Dubiella U, Khafif M, Alvarez E, Langin G, Roux F, Peyraud R, Roby D (2020) Robustness of plant quantitative disease resistance is provided by a decentralized immune network. Proc Natl Acad Sci U S A 117:18099-18109
Earl DA, vonHoldt BM (2012) STRUCTURE HARVESTER: a website and program for visualizing STRUCTURE output and implementing the Evanno method. Conserv Genet Res 4:359–361
El Houari I, Van Beirs C, Arents HE, Han H, Chanoca A, Opdenacker D, Pollier J, Storme V, Steenackers W, Quareshy M, Napier R, Beeckman T, Friml J, De Rybel B, Boerjan W, Vanholme B (2021) Seedling developmental defects upon blocking CINNAMATE-4-HYDROXYLASE are caused by perturbations in auxin transport. New Phytol 230:2275-2291
Garrido AN, Supijono E, Boshara P, Douglas SJ, Stronghill PE, Li B, Nambara E, Kliebenstein DJ, Riggs CD (2020) flasher, a novel mutation in a glucosinolate modifying enzyme, conditions changes in plant architecture and hormone homeostasis. Plant J 103:1989-2006
Goodwin SB (2002) The barley scald pathogen Rhynchosporium secalis is closely related to the discomycetes Tapesia and Pyrenopeziza. Mycol Res 106:645-654
Harper AL, Trick M, Higgins J, Fraser F, Clissold L, Wells R, Hattori C, Werner P, Bancroft I (2012) Associative transcriptomics of traits in the polyploid crop species Brassica napus. Nat Biotechnol 30:798-802
Havlickova L, He Z, Wang L, Langer S, Harper AL, Kaur H, Broadley MR, Gegas V, Bancroft I (2018) Validation of an updated Associative Transcriptomics platform for the polyploid crop species Brassica napus by dissection of the genetic architecture of erucic acid and tocopherol isoform variation in seeds. Plant J 93:181-192
Hejna O, Havlickova L, He Z, Bancroft I, Curn V (2019) Analysing the genetic architecture of clubroot resistance variation in Brassica napus by associative transcriptomics. Mol Breed 39:112
Huang YJ, Paillard S, Kumar V, King GJ, Fitt BDL, Delourme R (2019) Oilseed rape (Brassica napus) resistance to growth of Leptosphaeria maculans in leaves of young plants contributes to quantitative resistance in stems of adult plants. PLoS One 14:e0222540
Jakobsson M, Rosenberg NA (2007) CLUMPP: a cluster matching and permutation program for dealing with label switching and multimodality in analysis of population structure. Bioinformatics 23:1801-1806
Janes JK, Miller JM, Dupuis JR, Malenfant RM, Gorrell JC, Cullingham CI, Andrew RL (2017) The K = 2 conundrum. Mol Ecol 26:3594-3602
Jiquel A, Gervais J, Geistodt-Kiener A, Delourme R, Gay EJ, Ollivier B, Fudal I, Faure S, Balesdent MH, Rouxel T (2021) A gene-for-gene interaction involving a 'late' effector contributes to quantitative resistance to the stem canker disease in Brassica napus. New Phytol 231:1510-1524
Jones JDG, Dangl JL (2006) The plant immune system. Nature 444:323-329
Karandeni Dewage CS, Klöppel CA, Stotz HU, Fitt BDL (2018) Host–pathogen interactions in relation to management of light leaf spot disease (caused by Pyrenopeziza brassicae) on Brassica species. Crop Pasture Sci 69:9-19
Karandeni Dewage CS, Qi A, Stotz HU, Huang Y-J, Fitt BDL (2021) Interactions in the Brassica napus–Pyrenopeziza brassicae pathosystem and sources of resistance to P. brassicae (light leaf spot). Plant Pathol 70:2104-2114
Karandeni Dewage CS, Qi A, Stotz HU, Huang Y-J, Fitt BDL (2022) QTL mapping for resistance against Pyrenopeziza brassicae derived from a Brassica napus secondary gene pool. Front Plant Sci accepted
Larkan NJ, Lydiate DJ, Parkin IA, Nelson MN, Epp DJ, Cowling WA, Rimmer SR, Borhan MH (2013) The Brassica napus blackleg resistance gene LepR3 encodes a receptor-like protein triggered by the Leptosphaeria maculans effector AVRLM1. New Phytol 197:595-605
Larkan NJ, Ma L, Haddadi P, Buchwaldt M, Parkin IAP, Djavaheri M, Borhan MH (2020) The Brassica napus wall-associated kinase-like (WAKL) gene Rlm9 provides race-specific blackleg resistance. Plant J 104:892-900
Li J, Hansen BG, Ober JA, Kliebenstein DJ, Halkier BA (2008) Subclade of flavin-monooxygenases involved in aliphatic glucosinolate biosynthesis. Plant Physiol 148:1721-1733
Liu Y, Tabata D, Imai R (2016) A cold-Inducible DEAD-Box RNA helicase from Arabidopsis thaliana regulates plant growth and development under low temperature. PLoS One 11:e0154040
Looseley ME, Griffe LL, Buttner B, Wright KM, Bayer MM, Coulter M, Thauvin JN, Middlefell-Williams J, Maluk M, Okpo A, Kettles N, Werner P, Byrne E, Avrova A (2020) Characterisation of barley landraces from Syria and Jordan for resistance to rhynchosporium and identification of diagnostic markers for Rrs1Rh4. Theor Appl Genet 133:1243-1264
Lu K, Wei L, Li X, Wang Y, Wu J, Liu M, Zhang C, Chen Z, Xiao Z, Jian H, Cheng F, Zhang K, Du H, Cheng X, Qu C, Qian W, Liu L, Wang R, Zou Q, Ying J, Xu X, Mei J, Liang Y, Chai YR, Tang Z, Wan H, Ni Y, He Y, Lin N, Fan Y, Sun W, Li NN, Zhou G, Zheng H, Wang X, Paterson AH, Li J (2019) Whole-genome resequencing reveals Brassica napus origin and genetic loci involved in its improvement. Nat Commun 10:1154
Lu X, Dittgen J, Pislewska-Bednarek M, Molina A, Schneider B, Svatos A, Doubsky J, Schneeberger K, Weigel D, Bednarek P, Schulze-Lefert P (2015) Mutant allele-specific uncoupling of PENETRATION3 functions reveals engagement of the ATP-binding cassette transporter in distinct tryptophan metabolic pathways. Plant Physiol 168:814-827
Mason AS, Rousseau-Gueutin M, Morice J, Bayer PE, Besharat N, Cousin A, Pradhan A, Parkin IA, Chevre AM, Batley J, Nelson MN (2016) Centromere locations in Brassica A and C genomes revealed through half-tetrad analysis. Genetics 202:513-523
Nelson R, Wiesner-Hanks T, Wisser R, Balint-Kurti P (2018) Navigating complexity to breed disease-resistant crops. Nat Rev Genet 19:21-33
Ngou BPM, Ahn HK, Ding P, Jones JDG (2021) Mutual potentiation of plant immunity by cell-surface and intracellular receptors. Nature 592:110-115
Penselin D, Munsterkotter M, Kirsten S, Felder M, Taudien S, Platzer M, Ashelford K, Paskiewicz KH, Harrison RJ, Hughes DJ, Wolf T, Shelest E, Graap J, Hoffmann J, Wenzel C, Woltje N, King KM, Fitt BD, Guldener U, Avrova A, Knogge W (2016) Comparative genomics to explore phylogenetic relationship, cryptic sexual potential and host specificity of Rhynchosporium species on grasses. BMC Genomics 17:953
Perroud PF, Demko V, Ako AE, Khanal R, Bokor B, Pavlovic A, Jasik J, Johansen W (2021) The nuclear GUCT domain-containing DEAD-box RNA helicases govern gametophytic and sporophytic development in Physcomitrium patens. Plant Mol Biol 107:307-325
Perry J, Brachmann A, Welham T, Binder A, Charpentier M, Groth M, Haage K, Markmann K, Wang TL, Parniske M (2009) TILLING in Lotus japonicus identified large allelic series for symbiosis genes and revealed a bias in functionally defective ethyl methanesulfonate alleles toward glycine replacements. Plant Physiol 151:1281-1291
Pilet ML, Delourme R, Foisset N, Renard M (1998) Identification of QTL involved in field resistance to light leaf spot (Pyrenopeziza brassicae) and blackleg resistance (Leptosphaeria maculans) in winter rapeseed (Brassica napus L.). Theor Appl Genet 97:398-406
Qian L, Qian W, Snowdon RJ (2014) Sub-genomic selection patterns as a signature of breeding in the allopolyploid Brassica napus genome. BMC Genomics 15:1170
Rohe M, Gierlich A, Hermann H, Hahn M, Schmidt B, Rosahl S, Knogge W (1995) The race-specific elicitor, NIP1, from the barley pathogen, Rhynchosporium secalis, determines avirulence on host plants of the Rrs1 resistance genotype. EMBO J 14:4168-4177
Stein M, Dittgen J, Sanchez-Rodriguez C, Hou BH, Molina A, Schulze-Lefert P, Lipka V, Somerville S (2006) Arabidopsis PEN3/PDR8, an ATP binding cassette transporter, contributes to nonhost resistance to inappropriate pathogens that enter by direct penetration. Plant Cell 18:731-746
Storey JD, Bass AJ, Dabney A, Robinson D (2020) Q-value estimation for false discovery rate control. 2.22.0 edn, Bioconductor open source software for bioinformatics, p R package
Stotz HU, Jikumaru Y, Shimada Y, Sasaki E, Stingl N, Mueller MJ, Kamiya Y (2011a) Jasmonate-dependent and COI1-independent defense responses against Sclerotinia sclerotiorum in Arabidopsis thaliana: auxin is part of COI1-independent defense signaling. Plant Cell Physiol 52:1941-1956
Stotz HU, Mitrousia GK, de Wit PJ, Fitt BDL (2014) Effector-triggered defence against apoplastic fungal pathogens. Trends Plant Sci 19:491-500
Stotz HU, Sawada Y, Shimada Y, Hirai MY, Sasaki E, Krischke M, Brown PD, Saito K, Kamiya Y (2011b) Role of camalexin, indole glucosinolates, and side chain modification of glucosinolate-derived isothiocyanates in defense of Arabidopsis against Sclerotinia sclerotiorum. The Plant journal : for cell and molecular biology 67:81-93
Sylvester-Bradley R (1985) Revision of a code for stages of development in oilseed rape (Brassica napus L.). Asp Appl Biol 10:395–400
Vossen JH, Abd-El-Haliem A, Fradin EF, van den Berg GC, Ekengren SK, Meijer HJ, Seifi A, Bai Y, ten Have A, Munnik T, Thomma BP, Joosten MH (2010) Identification of tomato phosphatidylinositol-specific phospholipase-C (PI-PLC) family members and the role of PLC4 and PLC6 in HR and disease resistance. Plant J 62:224-239
Wells R, Trick M, Fraser F, Soumpourou E, Clissold L, Morgan C, Pauquet J, Bancroft I (2013) Sequencing-based variant detection in the polyploid crop oilseed rape. BMC Plant Biol 13:111
Wu J, Zhao Q, Liu S, Shahid M, Lan L, Cai G, Zhang C, Fan C, Wang Y, Zhou Y (2016) Genome-wide association study identifies new loci for resistance to sclerotinia stem rot in Brassica napus. Front Plant Sci 7:1418
Yuan M, Jiang Z, Bi G, Nomura K, Liu M, Wang Y, Cai B, Zhou JM, He SY, Xin XF (2021) Pattern-recognition receptors are required for NLR-mediated plant immunity. Nature 592:105-109
Zheng X, Koopmann B, Ulber B, von Tiedemann A (2020) A global survey on diseases and pests in oilseed rape—current challenges and innovative strategies of control. Front Agron 2:590908
Zheng XY, Spivey NW, Zeng W, Liu PP, Fu ZQ, Klessig DF, He SY, Dong X (2012) Coronatine promotes Pseudomonas syringae virulence in plants by activating a signaling cascade that inhibits salicylic acid accumulation. Cell Host Microbe 11:587-596

Table 1 to 4 is available in the Supplemental Files section.

PDF

Reviews received at journal
15 Feb, 2022
Reviewers invited by journal
10 Feb, 2022
Editor assigned by journal
05 Feb, 2022
First submitted to journal
04 Feb, 2022

You are reading this latest preprint version

Novel gene loci associated with susceptibility or cryptic quantitative resistance to Pyrenopeziza brassicae in Brassica napus

Status:

Version 1

Abstract

Figures

Key Messages

Introduction

Methods

Glasshouse growth conditions

Statistical analysis of phenotypic traits

Association transcriptomics

Results

LLS scoring system

Discussion

Population structure analysis reflects history of OSR cultivation

Putative function of GEMs in QDR

Conclusion

Declarations

References

Tables

Supplementary Files

Status:

Version 1