Genome-wide Association Study Identifies New Loci for Resistance to Sclerotinia Stem Rot in Brassica napus

Integrated Assays of Genome-Wide Association Study, Multi-Omics Co-Localization, and Machine Learning Associated Calcium Signaling Genes with Oilseed Rape Resistance to Sclerotinia sclerotiorum

Sclerotinia sclerotiorum (Ss) is one of the most devastating fungal pathogens, causing huge yield loss in multiple economically important crops including oilseed rape. Plant resistance to Ss pertains to quantitative disease resistance (QDR) controlled by multiple minor genes. Genome-wide identification of genes involved in QDR to Ss is yet to be conducted. In this study, we integrated several assays including genome-wide association study (GWAS), multi-omics co-localization, and machine learning prediction to identify, on a genome-wide scale, genes involved in the oilseed rape QDR to Ss. Employing GWAS and multi-omics co-localization, we identified seven resistance-associated loci (RALs) associated with oilseed rape resistance to Ss. Furthermore, we developed a machine learning algorithm and named it Integrative Multi-Omics Analysis and Machine Learning for Target Gene Prediction (iMAP), which integrates multi-omics data to rapidly predict disease resistance-related genes within a broad chromosomal region. Through iMAP based on the identified RALs, we revealed multiple calcium signaling genes related to the QDR to Ss. Population-level analysis of selective sweeps and haplotypes of variants confirmed the positive selection of the predicted calcium signaling genes during evolution. Overall, this study has developed an algorithm that integrates multi-omics data and machine learning methods, providing a powerful tool for predicting target genes associated with specific traits. Furthermore, it makes a basis for further understanding the role and mechanisms of calcium signaling genes in the QDR to Ss.

Natural variation in BnaA07.MKK9 confers resistance to Sclerotinia stem rot in oilseed rape

Sclerotinia stem rot (SSR), caused by the necrotrophic fungus Sclerotinia sclerotiorum, is one of the most devastating diseases for several major oil-producing crops. Despite its impact, the genetic basis of SSR resistance in plants remains poorly understood. Here, through a genome-wide association study, we identify a key gene, BnaA07. MKK9, that encodes a mitogen-activated protein kinase kinase that confers SSR resistance in oilseed rape. Our functional analyses reveal that BnaA07.MKK9 interacts with BnaC03.MPK3 and BnaC03.MPK6 and phosphorylates them at the TEY activation motif, triggering a signaling cascade that initiates biosynthesis of ethylene, camalexin, and indole glucosinolates, and promotes accumulation of H2O2 and the hypersensitive response, ultimately conferring resistance. Furthermore, variations in the coding sequence of BnaA07.MKK9 alter its kinase activity and improve SSR resistance by ~30% in cultivars carrying the advantageous haplotype. These findings enhance our understanding of SSR resistance and may help engineer novel diversity for future breeding of oilseed rape.

LESION MIMIC MUTANT 1 confers basal resistance to Sclerotinia sclerotiorum in rapeseed via a salicylic acid-dependent pathway

Rapeseed (Brassica napus) is a major edible oilseed crop consumed worldwide. However, its yield is seriously affected by infection from the broad-spectrum non-obligate pathogen Sclerotinia sclerotiorum due to a lack of highly resistant germplasm. Here, we identified a Sclerotinia-resistant and light-dependent lesion mimic mutant from an ethyl methanesulfonate-mutagenized population of the rapeseed inbred Zhongshuang 11 (ZS11) named lesion mimic mutant 1 (lmm1). The phenotype of lmm1 is controlled by a single recessive gene, named LESION MIMIC MUTANT 1 (LMM1), which mapped onto chromosome C04 by bulked segregant analysis within a 2.71-Mb interval. Histochemical analysis indicated that H2O2 strongly accumulated and cell death occurred around the lesion mimic spots. Among 877 differentially expressed genes (DEGs) between ZS11 and lmm1 leaves, 188 DEGs were enriched in the defense response, including 95 DEGs involved in systemic acquired resistance, which is consistent with the higher salicylic acid levels in lmm1. Combining bulked segregant analysis and transcriptome analysis, we identified a significantly up-regulated gene, BnaC4.PR2, which encodes β-1,3-glucanase, as the candidate gene for LMM1. Overexpression of BnaC4.PR2 may induce a reactive oxygen species burst to trigger partial cell death and systemic acquired resistance. Our study provides a new genetic resource for S. sclerotiorum resistance as well as new insights into disease resistance breeding in B. napus.

Genetic and molecular analysis of stem rot (Sclerotinia sclerotiorum) resistance in Brassica napus (canola type)

Identifying the molecular and genetic basis of resistance to Sclerotinia stem rot (Sclerotinia sclerotiorum) is critical for developing long-term and cost-effective management of this disease in rapeseed/canola (Brassica napus). Current cultural or chemical management options provide, at best, only partial and/or sporadic control. Towards this, a B. napus breeding population (Mystic x Rainbow), including the parents, F1, F2, BC1P1 and BC1P2, was utilized in a field study to determine the inheritance pattern of Sclerotinia stem rot resistance (based on stem lesion length, SLL). Broad sense heritability was 0.58 for SLL and 0.44 for days to flowering (DTF). There was a significant negative correlation between SLL and stem diameter (SD) (r = -0.39) and between SLL and DTF (r = -0.28), suggesting co-selection of SD and DTF traits, along with SLL, should assist in improving overall resistance. Non-additive genetic variance was evident for SLL, DTF, and SD. In a genome wide association study (GWAS), a significant quantitative trait locus (QTL) was identified for SLL. Several putative candidate marker trait associations (MTA) were located within this QTL region. Overall, this study has provided valuable new understanding of inheritance of resistance to S. sclerotiorum, and has identified QTL, MTAs and transgressive segregants with high-level resistances. Together, these will foster more rapid selection for multiple traits associated with Sclerotinia stem rot resistance, by enabling breeders to make critical choices towards selecting/developing cultivars with enhanced resistance to this devastating pathogen.

Genetic determinants of lettuce resistance to drop caused by Sclerotinia minor identified through genome-wide association mapping frequently co-locate with loci regulating anthocyanin content

KEY MESSAGE

GWAS identified 19 QTLs for resistance to Sclerotinia minor, 11 of them co-locating with red leaf color. Lower disease incidence was observed in red and dark red accessions. Lettuce (Lactuca sativa L.), one of the most economically important vegetables grown primarily in moderate climates around the world, is susceptible to many diseases including lettuce drop caused by the soilborne fungus Sclerotinia minor. Complete resistance to S. minor has not been identified in cultivated lettuce or its wild relatives. We conducted five experiments over 4 years with the diversity panel of almost 500 lettuce accessions to evaluate their response to the pathogen in an artificially infested field. The lowest disease incidence (DI) was observed in cultivars Eruption, Infantry, and Annapolis (median DI of 12.1-17.5%), while the highest DI was recorded for cultivars Reine des Glaces, Wayahead, and line FL. 43007 (median DI of 81.0-95.2%). Overall, significantly lower DI was observed in red and dark red accessions compared to those with a lower anthocyanin content. Genome-wide association mapping identified 19 QTLs for resistance to S. minor, 21 for the presence of red leaf color or its variations caused by the anthocyanin content, and one for the green color intensity. Eleven of the QTLs for disease resistance were located within 10 Mb of the loci associated with red color or anthocyanin content identified in this diversity panel. The frequent, non-random co-location of QTLs, together with the lower DI observed in red and dark red accessions suggests that lettuce interaction with S. minor may be partly influenced by anthocyanins. We have identified RLL2 and ANS, the genes of the anthocyanin biosynthesis pathway that co-locate with resistance QTLs, as candidates for functional studies to ascertain the involvement of anthocyanins in lettuce resistance against S. minor. Resistance QTLs closely linked with QTLs for anthocyanin content could be used to develop lettuce with a relatively high partial resistance and red color, while those not associated with anthocyanins could be used to develop partially resistant cultivars of green color.

Subgenome-dominant expression and alternative splicing in response to Sclerotinia infection in polyploid Brassica napus and progenitors

Polyploidy has played an extensive role in the evolution of flowering plants. Allopolyploids, with subgenomes containing duplicated gene pairs called homeologs, can show rapid transcriptome changes including novel alternative splicing (AS) patterns. The extent to which abiotic stress modulates AS of homeologs is a nascent topic in polyploidy research. We subjected both resynthesized and natural lines of polyploid Brassica napus, along with the progenitors Brassica rapa and Brassica oleracea, to infection with the fungal pathogen Sclerotinia sclerotiorum. RNA-sequencing analyses revealed widespread divergence between polyploid subgenomes in both gene expression and AS patterns. Resynthesized B. napus displayed significantly more A and C subgenome biased homeologs under pathogen infection than during uninfected growth. Differential AS (DAS) in response to infection was highest in natural B. napus (12 709 DAS events) and lower in resynthesized B. napus (8863 DAS events). Natural B. napus had more upregulated events and fewer downregulated events. There was a global expression bias towards the B. oleracea-derived (C) subgenome in both resynthesized and natural B. napus, enhanced by widespread non-parental downregulation of the B. rapa-derived (A) homeolog. In the resynthesized B. napus, this resulted in a disproportionate C subgenome contribution to the pathogen defense response, characterized by biases in both transcript expression levels and the proportion of induced genes. Our results elucidate the complex ways in which Sclerotinia infection affects expression and AS of homeologous genes in resynthesized and natural B. napus.

QTL mapping and transcriptome analysis identify novel QTLs and candidate genes in Brassica villosa for quantitative resistance against Sclerotinia sclerotiorum

Novel QTLs and candidate genes for Sclerotinia-resistance were identified in B. villosa, a wild Brassica species, which represents a new genetic source for improving oilseed rape resistance to SSR. Sclerotinia stem rot (SSR), caused by Sclerotinia sclerotiorum, is one of the most destructive diseases in oilseed rape growing regions. To date, there is no effective genetic resistance against S. sclerotiorum in the B. napus germplasm and knowledge of the molecular plant-fungal interaction is also limited. To identify new resistance resources, we screened a set of wild Brassica species and identified B. villosa (BRA1896) with a high level of Sclerotinia-resistance. Two segregating F₂ populations for Sclerotinia-resistance, generated by interspecific crosses between the resistant B. villosa (BRA1896) and the wild susceptible B. oleracea (BRA1909) were assessed for Sclerotinia-resistance. Genetic mapping using a 15-k Illumina Infinium SNP-array resulted in a high-density genetic map containing 1,118 SNP markers and spanning a total genetic length of 792.2 cM. QTL analysis revealed seven QTLs explaining 3.8% to 16.5% of phenotypic variance. Intriguingly, RNAseq-based transcriptome analysis identified genes and pathways specific to B. villosa, of which a cluster of five genes encoding putative receptor-like kinases (RLKs) and two pathogenesis-related (PR) proteins are co-localized within a QTL on chromosome C07. Furthermore, transcriptomic analysis revealed enhanced ethylene (ET)-activated signaling in the resistant B. villosa, which is associated with a stronger plant immune response, depressed cell death, and enhanced phytoalexin biosynthesis compared to the susceptible B. oleracea. Our data demonstrates that B. villosa represents a novel and unique genetic source for improving oilseed rape resistance against SSR.

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

Quantitative disease resistance (QDR) controls the association of the light leaf spot pathogen with Brassica napus; four QDR loci that were in linkage disequilibrium and eight gene expression markers were identified. 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 195 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.

Identification of candidate genes for LepR1 resistance against Leptosphaeria maculans in Brassica napus

Utilising resistance (R) genes, such as LepR1, against Leptosphaeria maculans, the causal agent of blackleg in canola (Brassica napus), could help manage the disease in the field and increase crop yield. Here we present a genome wide association study (GWAS) in B. napus to identify LepR1 candidate genes. Disease phenotyping of 104 B. napus genotypes revealed 30 resistant and 74 susceptible lines. Whole genome re-sequencing of these cultivars yielded over 3 million high quality single nucleotide polymorphisms (SNPs). GWAS in mixed linear model (MLM) revealed a total of 2,166 significant SNPs associated with LepR1 resistance. Of these SNPs, 2108 (97%) were found on chromosome A02 of B. napus cv. Darmor bzh v9 with a delineated LepR1_mlm1 QTL at 15.11-26.08 Mb. In LepR1_mlm1, there are 30 resistance gene analogs (RGAs) (13 nucleotide-binding site-leucine rich repeats (NLRs), 12 receptor-like kinases (RLKs), and 5 transmembrane-coiled-coil (TM-CCs)). Sequence analysis of alleles in resistant and susceptible lines was undertaken to identify candidate genes. This research provides insights into blackleg resistance in B. napus and assists identification of the functional LepR1 blackleg resistance gene.

Comparative transcriptome profiling reveals the importance of GmSWEET15 in soybean susceptibility to Sclerotinia sclerotiorum

Soybean sclerotinia stem rot (SSR) is a disease caused by Sclerotinia sclerotiorum that causes incalculable losses in soybean yield each year. Considering the lack of effective resistance resources and the elusive resistance mechanisms, we are urged to develop resistance genes and explore their molecular mechanisms. Here, we found that loss of GmSWEET15 enhanced the resistance to S. sclerotiorum, and we explored the molecular mechanisms by which gmsweet15 mutant exhibit enhanced resistance to S. sclerotiorum by comparing transcriptome. At the early stage of inoculation, the wild type (WT) showed moderate defense response, whereas gmsweet15 mutant exhibited more extensive and intense transcription reprogramming. The gmsweet15 mutant enriched more biological processes, including the secretory pathway and tetrapyrrole metabolism, and it showed stronger changes in defense response, protein ubiquitination, MAPK signaling pathway-plant, plant-pathogen interaction, phenylpropanoid biosynthesis, and photosynthesis. The more intense and abundant transcriptional reprogramming of gmsweet15 mutant may explain how it effectively delayed colonization by S. sclerotiorum. In addition, we identified common and specific differentially expressed genes between WT and gmsweet15 mutant after inoculation with S. sclerotiorum, and gene sets and genes related to gmsweet15_24 h were identified through Gene Set Enrichment Analysis. Moreover, we constructed the protein-protein interaction network and gene co-expression networks and identified several groups of regulatory networks of gmsweet15 mutant in response to S. sclerotiorum, which will be helpful for the discovery of candidate functional genes. Taken together, our results elucidate molecular mechanisms of delayed colonization by S. sclerotiorum after loss of GmSWEET15 in soybean, and we propose novel resources for improving resistance to SSR.

Quantitative trait locus mapping and improved resistance to sclerotinia stem rot in a backbone parent of rapeseed (Brassica napus L.)

There are three main challenges to improving sclerotinia stem rot (SSR) resistance in rapeseed (Brassica napus L.). First, breeding materials such as the backbone parents have not been extensively investigated, making the findings of previous studies difficult to directly implement. Second, SSR resistance and flowering time (FT) loci are typically linked; thus, use of these loci requires sacrifice of the rapeseed growth period. Third, the SSR resistance loci in susceptible materials are often neglected, thereby reducing the richness of resistant resources. This study was conducted to investigate the stem resistance, disease index, and FT of a doubled haploid population consisting of 151 lines constructed from the backbone parent 19514A and conventional rapeseed cultivar ZY50 within multiple environments. Quantitative trait locus (QTL) mapping revealed 13 stem resistance QTLs, 9 disease index QTLs, and 20 FT QTLs. QTL meta-analysis showed that uqA04, uqC03.1, and uqC03.2 were repeatable SSR resistance QTLs derived from different parents but not affected by the FT. Based on these three QTLs, we proposed a strategy for improving the SSR resistance of 19514A and ZY50. This study improves the understanding of the resistance to rapeseed SSR and genetic basis of FT and demonstrates that SSR resistance QTLs can be mined from parents with a minimal resistance level difference, thereby supporting the application of backbone parents in related research and resistance improvement.

The Sclerotinia sclerotiorum-inducible promoter pBnGH17D7 in Brassica napus: isolation, characterization, and application in host-induced gene silencing

Sclerotinia stem rot (SSR), caused by Sclerotinia sclerotiorum, is among the most devastating diseases in Brassica napus worldwide. Conventional breeding for SSR resistance in Brassica species is challenging due to the limited availability of resistant germplasm. Therefore, genetic engineering is an attractive approach for developing SSR-resistant Brassica crops. Compared with the constitutive promoter, an S. sclerotiorum-inducible promoter would avoid ectopic expression of defense genes that may cause plant growth deficits. In this study, we generated a S. sclerotiorum-inducible promoter. pBnGH17D7, from the promoter of B. napus glycosyl hydrolase 17 gene (pBnGH17). Specifically, 5'-deletion and promoter activity analyses in transgenic Arabidopsis thaliana plants defined a 189 bp region of pBnGH17 which was indispensable for S. sclerotiorum-induced response. Compared with pBnGH17, pBnGH17D7 showed a similar response upon S. sclerotiorum infection, but lower activity in plant tissues in the absence of S. sclerotiorum infection. Moreover, we revealed that the transcription factor BnTGA7 directly binds to the TGACG motif in pBnGH17D7 to activate BnGH17. Ultimately, pBnGH17D7 was exploited for engineering Sclerotinia-resistant B. napus via host-induced gene silencing. It induces high expression of siRNAs against the S. sclerotiorum pathogenic factor gene specifically during infection, leading to increased resistance.

Genome-wide association study, haplotype analysis, and genomic prediction reveal the genetic basis of yield-related traits in soybean (Glycine max L.)

Identifying the genetic components underlying yield-related traits in soybean is crucial for improving its production and productivity. Here, 211 soybean genotypes were evaluated across six environments for four yield-related traits, including seed yield per plant (SYP), number of pods per plant number of seeds per plant and 100-seed weight (HSW). Genome-wide association study (GWAS) and genomic prediction (GP) analyses were performed using 12,617 single nucleotide polymorphism markers from NJAU 355K SoySNP Array. A total of 57 SNPs were significantly associated with four traits across six environments and a combined environment using five Genome-wide association study models. Out of these, six significant SNPs were consistently identified in more than three environments using multiple GWAS models. The genomic regions (±670 kb) flanking these six consistent SNPs were considered stable QTL regions. Gene annotation and in silico expression analysis revealed 15 putative genes underlying the stable QTLs that might regulate soybean yield. Haplotype analysis using six significant SNPs revealed various allelic combinations regulating diverse phenotypes for the studied traits. Furthermore, the GP analysis revealed that accurate breeding values for the studied soybean traits is attainable at an earlier generation. Our study paved the way for increasing soybean yield performance within a short breeding cycle.

The evolutionary and molecular features of the broad-host-range plant pathogen Sclerotinia sclerotiorum

Sclerotinia sclerotiorum is a pathogenic fungus that infects hundreds of plant species, including many of the world's most important crops. Key features of S. sclerotiorum include its extraordinary host range, preference for dicotyledonous plants, relatively slow evolution, and production of protein effectors that are active in multiple host species. Plant resistance to this pathogen is highly complex, typically involving numerous polymorphisms with infinitesimally small effects, which makes resistance breeding a major challenge. Due to its economic significance, S. sclerotiorum has been subjected to a large amount of molecular and evolutionary research. In this updated pathogen profile, we review the evolutionary and molecular features of S. sclerotiorum and discuss avenues for future research into this important species.

Unraveling the tripartite interaction of volatile compounds of Streptomyces rochei with grain mold pathogens infecting sorghum

Sorghum is a major grain crop used in traditional meals and health drinks, and as an efficient fuel. However, its productivity, value, germination, and usability are affected by grain mold, which is a severe problem in sorghum production systems, which reduces the yield of harvested grains for consumer use. The organic approach to the management of the disease is essential and will increase consumer demand. Bioactive molecules like mVOC (volatile organic compound) identification are used to unravel the molecules responsible for antifungal activity. The Streptomyces rochei strain (ASH) has been reported to be a potential antagonist to many pathogens, with high levels of VOCs. The present study aimed to study the inhibitory effect of S. rochei on sorghum grain mold pathogens using a dual culture technique and via the production of microbial volatile organic compounds (mVOCs). mVOCs inhibited the mycelial growth of Fusarium moniliforme by 63.75 and Curvularia lunata by 68.52%. mVOCs suppressed mycelial growth and inhibited the production of spores by altering the structure of mycelia in tripartite plate assay. About 45 mVOCs were profiled when Streptomyces rochei interacted with these two pathogens. In the present study, several compounds were upregulated or downregulated by S. rochei, including 2-methyl-1-butanol, methanoazulene, and cedrene. S. rochei emitted novel terpenoid compounds with peak areas, such as myrcene (1.14%), cymene (6.41%), and ç-terpinene (7.32%) upon interaction with F. moniliforme and C. lunata. The peak area of some of the compounds, including furan 2-methyl (0.70%), benzene (1.84%), 1-butanol, 2-methyl-(8.25%), and myrcene (1.12)%, was increased during tripartite interaction with F. moniliforme and C. lunata, which resulted in furan 2-methyl (6.60%), benzene (4.43%), butanol, 2-methyl (18.67%), and myrcene (1.14%). These metabolites were implicated in the sesquiterpenoid and alkane biosynthetic pathways and the oxalic acid degradation pathway. The present study shows how S. rochei exhibits hyperparasitism, competition, and antibiosis via mVOCs. In addition to their antimicrobial functions, these metabolites could also enhance plant growth.

Genome-wide identification and functional exploration of the legume lectin genes in Brassica napus and their roles in Sclerotinia disease resistance

As one of the largest classes of lectins, legume lectins have a variety of desirable features such as antibacterial and insecticidal activities as well as anti-abiotic stress ability. The Sclerotinia disease (SD) caused by the soil-borne fungus Sclerotinia sclerotiorum is a devastating disease affecting most oil crops such as Brassica napus. Here, we identified 130 legume lectin (LegLu) genes in B. napus, which could be phylogenetically classified into seven clusters. The BnLegLu gene family has been significantly expanded since the whole-genome duplication (WGD) or segmental duplication. Gene structure and conserved motif analysis suggested that the BnLegLu genes were well conserved in each cluster. Moreover, relative to those genes only containing the legume lectin domain in cluster VI-VII, the genes in cluster I-V harbored a transmembrane domain and a kinase domain linked to the legume lectin domain in the C terminus. The expression of most BnLegLu genes was relatively low in various tissues. Thirty-five BnLegLu genes were responsive to abiotic stress, and 40 BnLegLu genes were strongly induced by S. sclerotiorum, with a most significant up-regulation of 715-fold, indicating their functional roles in SD resistance. Four BnLegLu genes were located in the candidate regions of genome-wide association analysis (GWAS) results which resulted from a worldwide rapeseed population consisting of 324 accessions associated with SD. Among them, the positive role of BnLegLus-16 in SD resistance was validated by transient expression in tobacco leaves. This study provides important information on BnLegLu genes, particularly about their roles in SD resistance, which may help targeted functional research and genetic improvement in the breeding of B. napus.

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

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

Supplementary Information

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

Genetic mapping and genomic prediction of sclerotinia stem rot resistance to rapeseed/canola (Brassica napus L.) at seedling stage

GWAS detected ninety-eight significant SNPs associated with Sclerotinia sclerotiorum resistance. Six statistical models resulted in medium to high predictive ability, depending on trait, indicating potential of genomic prediction for disease resistance breeding. The lack of complete host resistance and a complex resistance inheritance nature between rapeseed/canola and Sclerotinia sclerotiorum often limits the development of functional molecular markers that enable breeding for sclerotinia stem rot (SSR) resistance. However, genomics-assisted selection has the potential to accelerate the breeding for SSR resistance. Therefore, genome-wide association (GWA) mapping and genomic prediction (GP) were performed using a diverse panel of 337 rapeseed/canola genotypes. Three-week-old seedlings were screened using the petiole inoculation technique (PIT). Days to wilt (DW) up to 2 weeks and lesion phenotypes (LP) at 3, 4, and 7 days post-inoculation (dpi) were recorded. A strong correlation (r = - 0.90) between DW and LP_4dpi implied that a single time point scoring at four days could be used as a proxy trait. GWA analyses using single-locus (SL) and multi-locus (ML) models identified a total of 41, and 208 significantly associated SNPs, respectively. Out of these, ninety-eight SNPs were identified by a combination of the SL model and any of the ML models, at least two ML models, or two traits. These SNPs explained 1.25-12.22% of the phenotypic variance and considered as significant, could be associated with SSR resistance. Eighty-three candidate genes with a function in disease resistance were associated with the significant SNPs. Six GP models resulted in moderate to high (0.42-0.67) predictive ability depending on SSR resistance traits. The resistant genotypes and significant SNPs will serve as valuable resources for future SSR resistance breeding. Our results also highlight the potential of genomic selection to improve rapeseed/canola breeding for SSR resistance.

Mining of Cloned Disease Resistance Gene Homologs (CDRHs) in Brassica Species and Arabidopsis thaliana

Simple Summary

Developing cultivars with resistance genes (R genes) is an effective strategy to support high yield and quality in Brassica crops. The availability of clone R gene and genomic sequences in Brassica species and Arabidopsis thaliana provide the opportunity to compare genomic regions and survey R genes across genomic databases. In this paper, we aim to identify genes related to cloned genes through sequence identity, providing a repertoire of species-wide related R genes in Brassica crops. The comprehensive list of candidate R genes can be used as a reference for functional analysis.

Abstract

Various diseases severely affect Brassica crops, leading to significant global yield losses and a reduction in crop quality. In this study, we used the complete protein sequences of 49 cloned resistance genes (R genes) that confer resistance to fungal and bacterial diseases known to impact species in the Brassicaceae family. Homology searches were carried out across Brassica napus, B. rapa, B. oleracea, B. nigra, B. juncea, B. carinata and Arabidopsis thaliana genomes. In total, 660 cloned disease R gene homologs (CDRHs) were identified across the seven species, including 431 resistance gene analogs (RGAs) (248 nucleotide binding site-leucine rich repeats (NLRs), 150 receptor-like protein kinases (RLKs) and 33 receptor-like proteins (RLPs)) and 229 non-RGAs. Based on the position and distribution of specific homologs in each of the species, we observed a total of 87 CDRH clusters composed of 36 NLR, 16 RLK and 3 RLP homogeneous clusters and 32 heterogeneous clusters. The CDRHs detected consistently across the seven species are candidates that can be investigated for broad-spectrum resistance, potentially providing resistance to multiple pathogens. The R genes identified in this study provide a novel resource for the future functional analysis and gene cloning of Brassicaceae R genes towards crop improvement.

Construction of a Quantitative Genomic Map, Identification and Expression Analysis of Candidate Genes for Agronomic and Disease-Related Traits in Brassica napus

Rapeseed is the second most important oil crop in the world. Improving seed yield and seed oil content are the two main highlights of the research. Unfortunately, rapeseed development is frequently affected by different diseases. Extensive research has been made through many years to develop elite cultivars with high oil, high yield, and/or disease resistance. Quantitative trait locus (QTL) analysis has been one of the most important strategies in the genetic deciphering of agronomic characteristics. To comprehend the distribution of these QTLs and to uncover the key regions that could simultaneously control multiple traits, 4,555 QTLs that have been identified during the last 25 years were aligned in one unique map, and a quantitative genomic map which involved 128 traits from 79 populations developed in 12 countries was constructed. The present study revealed 517 regions of overlapping QTLs which harbored 2,744 candidate genes and might affect multiple traits, simultaneously. They could be selected to customize super-rapeseed cultivars. The gene ontology and the interaction network of those candidates revealed genes that highly interacted with the other genes and might have a strong influence on them. The expression and structure of these candidate genes were compared in eight rapeseed accessions and revealed genes of similar structures which were expressed differently. The present study enriches our knowledge of rapeseed genome characteristics and diversity, and it also provided indications for rapeseed molecular breeding improvement in the future.

QTL Identification for Stem Fiber, Strength and Rot Resistance in a DH Population from an Alien Introgression of Brassica napus

Stem fiber, stem strength and stem-rot resistance are important agronomic traits in Brassica napus. To understand the molecular mechanism that controls the stem-related traits, we investigated the stem lignin (ADL), cellulose (Cel), hemicellulose (Hem) content, S/G monolignol ratio (SG), stem breaking force (BF), breaking strength (F) and Sclerotinia sclerotiorum resistance (SSR). Each trait was significantly positively or negatively correlated with more than three of the other six traits. QTL mapping for ADL, Cel, Hem, SG, BF, F and SSR were performed using a doubled haploid population derived from an intertribal B. napus introgression line 'Y689' crossed with B. napus cv. 'Westar'. A total of 67 additive QTL were identified and integrated into 55 consensus QTL by meta-analysis. Among the 55 consensus QTL, 23 (41.8%) QTL were co-located and were integrated into 11 unique QTL. The QTL by environment (Q × E) interactions were analyzed and 22 combined QTL were identified. In addition, candidate genes within the QTL intervals were proposed based on the known function of Arabidopsis orthologs. These results provided valuable information for improving lodging resistance, S. sclerotiorum resistance and mechanized harvesting of B. napus.

Two decades of association mapping: Insights on disease resistance in major crops

Climate change across the globe has an impact on the occurrence, prevalence, and severity of plant diseases. About 30% of yield losses in major crops are due to plant diseases; emerging diseases are likely to worsen the sustainable production in the coming years. Plant diseases have led to increased hunger and mass migration of human populations in the past, thus a serious threat to global food security. Equipping the modern varieties/hybrids with enhanced genetic resistance is the most economic, sustainable and environmentally friendly solution. Plant geneticists have done tremendous work in identifying stable resistance in primary genepools and many times other than primary genepools to breed resistant varieties in different major crops. Over the last two decades, the availability of crop and pathogen genomes due to advances in next generation sequencing technologies improved our understanding of trait genetics using different approaches. Genome-wide association studies have been effectively used to identify candidate genes and map loci associated with different diseases in crop plants. In this review, we highlight successful examples for the discovery of resistance genes to many important diseases. In addition, major developments in association studies, statistical models and bioinformatic tools that improve the power, resolution and the efficiency of identifying marker-trait associations. Overall this review provides comprehensive insights into the two decades of advances in GWAS studies and discusses the challenges and opportunities this research area provides for breeding resistant varieties.

Quantitative Inheritance of Sclerotinia Stem Rot Resistance in Brassica napus and Relationship to Cotyledon and Leaf Resistances

Sclerotinia sclerotiorum is a necrotrophic fungus causing devastating stem rot and associated yield losses of canola/rapeseed (Brassica napus) worldwide, including in Australia. Developing host resistance against Sclerotinia stem rot is critical if this disease in canola/rapeseed is to be successfully managed, as cultural or chemical control options provide only partial or sporadic control. Three B. napus breeding populations, C2, C5 and C6, including the parents, F₁, F₂, BC₁P1, and BC₂P2, were used in a field study with an objective of exploring the inheritance pattern of disease resistance (based on stem lesion length [SLL]) and the genetic relationships of disease with stem diameter (SD) or days to first flowering (DTF), and to compare these new adult plant stem resistances against S. sclerotiorum with those of seedling (cotyledon and leaf) resistances in earlier studies. Heritability (broad sense) of SLL was 0.57 and 0.73 for population C2 at 3 and 5 weeks postinoculation and 0.21 for population C5 at 5 weeks postinoculation. Additive genetic variance was evident within all 3 populations for DTF but not for SD. Narrow-sense heritability for DTF was 0.48 (C2), 0.42 (C5), and 0.32 (C6). SD, DTF, and SLL were all inherited independently, with no significant genetic covariance between traits in bivariate analysis. Genetic variance for SLL in populations C2 and C5 was entirely nonadditive, and there was significant nonadditive genetic covariance of SLL at 3 and 5 weeks postinoculation. Generation means analysis in population C2 supported the conclusion that complex epistatic interactions controlled SLL. Several C2 and C5 progeny showed high adult plant stem resistance, which may be critical in developing enhanced stem resistance in canola/rapeseed. Although population C6 showed no genetic variation for SLL resistance in this study, it showed significant nonadditive genetic variance at the cotyledon and leaf stages in earlier studies. We conclude that host resistance varies across different plant growth stages, and breeding must be targeted for resistance at each growth stage. In populations C2, C5, and C6, resistance to S. sclerotiorum in stem, leaf, and cotyledon was always controlled by nonadditive effects such as complex epistasis or dominance. Overall, our findings in relation to the quantitative inheritance of Sclerotinia stem rot resistance, together with the new high-level resistances identified, will enable breeders to select/develop genotypes with enhanced resistances to S. sclerotiorum.

Utilization of a Sugarcane100K Single Nucleotide Polymorphisms Microarray-Derived High-Density Genetic Map in Quantitative Trait Loci Mapping and Function Role Prediction of Genes Related to Chlorophyll Content in Sugarcane

Chlorophyll is the most important pigment for plant photosynthesis that plays an important role in crop growth and production. In this study, the chlorophyll content trait was explored to improve sugarcane yield. Two hundred and eighty-five F₁ progenies from the cross YT93-159 × ROC22 with significantly different chlorophyll contents were included as test materials. The chlorophyll content of the +1 leaves during elongation phase was measured using a SPAD-502 meter through a three-crop cycle (plant cane, first ratoon, and second ratoon). Linkage analysis was conducted on a high-density genetic map constructed based on the sugarcane 100K SNP chip. In addition, Fv/Fm, plant height, stalk diameter, brix data were collected on plant cane during the elongation and maturation phases. The results showed that the +1 leaf SPAD values, which can be used as an important reference to evaluate the growth potential of sugarcane, were significantly and positively correlated with the Fv/Fm during elongation phase, as well as with plant height, stalk diameter, and brix during maturity phase (P < 0.01). The broad sense heritability (H ²) of the chlorophyll content trait was 0.66 for plant cane crop, 0.67 for first ratoon crop, and 0.73 for second ratoon crop, respectively, indicating that this trait was mainly controlled by genetic factors. Thirty-one quantitative trait loci (QTL) were detected by QTL mapping. Among them, a major QTL, qCC-R1, could account for 12.95% of phenotypic variation explained (PVE), and the other 30 minor QTLs explained 2.37-7.99% PVE. Twenty candidate genes related to chlorophyll content were identified in the QTLs plus a 200-Kb extension region within either sides, of which four were homologous genes involved in the chlorophyll synthesis process and the remaining 16 played a certain role in chlorophyll catabolic pathway, chloroplast organization, or photosynthesis. These results provide a theoretical reference for analyzing the genetic mechanism of chlorophyll synthesis and subsequent improvement of photosynthetic characteristics in sugarcane.

Genome-wide association mapping and genomic prediction for adult stage sclerotinia stem rot resistance in Brassica napus (L) under field environments

Sclerotinia stem rot (SSR) is a fungal disease of rapeseed/canola that causes significant seed yield losses and reduces its oil content and quality. In the present study, the reaction of 187 diverse canola genotypes to SSR was characterized at full flowering stage using the agar plug to stem inoculation method in four environments. Genome-wide association study (GWAS) using three different algorithms identified 133 significant SNPs corresponding with 123 loci for disease traits like stem lesion length (LL), lesion width (LW), and plant mortality at 14 (PM_14D) and 21 (PM_21D) days. The explained phenotypic variation of these SNPs ranged from 3.6 to 12.1%. Nineteen significant SNPs were detected in two or more environments, disease traits with at least two GWAS algorithms. The strong correlations observed between LL and other three disease traits evaluated, suggest they could be used as proxies for SSR resistance phenotyping. Sixty-nine candidate genes associated with disease resistance mechanisms were identified. Genomic prediction (GP) analysis with all the four traits employing genome-wide markers resulted in 0.41-0.64 predictive ability depending on the model specifications. The highest predictive ability for PM_21D with three models was about 0.64. From our study, the identified resistant genotypes and stable significant SNP markers will serve as a valuable resource for future SSR resistance breeding. Our study also suggests that genomic selection holds promise for accelerating canola breeding progress by enabling breeders to select SSR resistance genotypes at the early stage by reducing the need to phenotype large numbers of genotypes.

The Bnapus50K array: a quick and versatile genotyping tool for Brassica napus genomic breeding and research

Rapeseed is a globally cultivated commercial crop, primarily grown for its oil. High-density single nucleotide polymorphism (SNP) arrays are widely used as a standard genotyping tool for rapeseed research, including for gene mapping, genome-wide association studies, germplasm resource analysis, and cluster analysis. Although considerable rapeseed genome sequencing data have been released, DNA arrays are still an attractive choice for providing additional genetic data in an era of high-throughput whole-genome sequencing. Here, we integrated re-sequencing DNA array data (32,216, 304 SNPs) from 505 inbred rapeseed lines, allowing us to develop a sensitive and efficient genotyping DNA array, Bnapus50K, with a more consistent genetic and physical distribution of probes. A total of 42,090 high-quality probes were filtered and synthesized, with an average distance between adjacent SNPs of 8 kb. To improve the practical application potential of this array in rapeseed breeding, we also added 1,618 functional probes related to important agronomic traits such as oil content, disease resistance, male sterility, and flowering time. The additional probes also included those specifically for detecting genetically modified material. These probes show a good detection efficiency and are therefore useful for gene mapping, along with crop variety improvement and identification. The novel Bnapus50K DNA array developed in this study could prove to be a quick and versatile genotyping tool for B. napus genomic breeding and research.

Brassica napus Mediator Subunit16 Induces BnMED25- and BnWRKY33-Activated Defense Signaling to Confer Sclerotinia sclerotiorum Resistance

The plant mediator is a highly conserved protein complex that interacts with transcription factors (TFs) and RNA polymerase II (RNAP II) to relay regulatory information during transcription. Plant immune response is one of the biological processes that is orchestrated by this regulatory mechanism. Brassica napus, an important oil crop, is severely attacked by a devastating disease Sclerotinia stem rot. Here, we explored broad-spectrum disease resistant roles of B. napus mediator subunit 16 (BnMED16) and its host defense mechanism against fugal pathogen Sclerotinia sclerotiorum. We found that BnMED16 expression was significantly increased by S. sclerotiorum infection, and its homologous overexpression resulted in rapid and comprehensive defense responses from the beginning to the end. This affected signal transduction with multiple channels including pathogen recognition, intracellular Ca²⁺ concentration, reactive oxygen species (ROS) accumulation and clearance, and activation of mitogen-activated protein kinase (MAPK) signaling cascades initially. Subsequently, pathogen-/defense-related genes and hormone-responsive pathways were highly activated, which resulted in enhanced cell wall and secretion of defense proteases. Furthermore, the biochemical analysis showed that BnMED16 interacts with BnMED25 and BnWRKY33. Additionally, BnMED25 also interacts with TFs BnMYC2, BnCOI1, and BnEIN3 of the JA/ET signal transduction pathway. Taken together, we proposed a hypothetical model that BnMED16 confers S. sclerotiorum resistance by enhancing BnMED25-mediated JA/ET defense pathways and BnWRKY33-activated defense signaling in B. napus. The BnMED16 overexpressing lines with enhanced broad-spectrum disease resistance could be useful for breeding Sclerotinia-resistant oilseed rape varieties, as well as serving as basis for further strategy development in resistance breeding.

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.

Sclerotinia Stem Rot Resistance in Rapeseed: Recent Progress and Future Prospects

Sclerotinia stem rot (SSR) of rapeseed (Brassica napus), caused by the soil-borne fungus Sclerotinia sclerotiorum, is one of the main diseases seriously affecting the yield and oil quality of infected rapeseed crops. The complexity of the inheritance of resistance and of the interaction mechanisms between rapeseed and S. sclerotiorum limits resistance gene identification and molecular breeding. In this review, the latest progress of research into resistance to SSR in B. napus is summarized from the following three directions: the pathogenesis mechanisms of S. sclerotiorum, the resistance mechanisms of B. napus toward S. sclerotiorum, and rapeseed breeding for resistance to SSR. This review aims to provide a theoretical basis and useful reference for analyzing the mechanism of the interaction between B. napus and S. sclerotiorum, searching for gene loci associated with the resistance response, and for achieving disease-resistance genetic manipulation and molecular design breeding in rapeseed.

Understanding Host-Pathogen Interactions in Brassica napus in the Omics Era

Brassica napus (canola/oilseed rape/rapeseed) is an economically important crop, mostly found in temperate and sub-tropical regions, that is cultivated widely for its edible oil. Major diseases of Brassica crops such as Blackleg, Clubroot, Sclerotinia Stem Rot, Downy Mildew, Alternaria Leaf Spot and White Rust have caused significant yield and economic losses in rapeseed-producing countries worldwide, exacerbated by global climate change, and, if not remedied effectively, will threaten global food security. To gain further insights into the host-pathogen interactions in relation to Brassica diseases, it is critical that we review current knowledge in this area and discuss how omics technologies can offer promising results and help to push boundaries in our understanding of the resistance mechanisms. Omics technologies, such as genomics, proteomics, transcriptomics and metabolomics approaches, allow us to understand the host and pathogen, as well as the interaction between the two species at a deeper level. With these integrated data in multi-omics and systems biology, we are able to breed high-quality disease-resistant Brassica crops in a more holistic, targeted and accurate way.

Exploring sunflower responses to Sclerotinia head rot at early stages of infection using RNA-seq analysis

Sclerotinia head rot (SHR), caused by the necrotrophic fungus Sclerotinia sclerotiorum, is one of the most devastating sunflower crop diseases. Despite its worldwide occurrence, the genetic determinants of plant resistance are still largely unknown. Here, we investigated the Sclerotinia-sunflower pathosystem by analysing temporal changes in gene expression in one susceptible and two tolerant inbred lines (IL) inoculated with the pathogen under field conditions. Differential expression analysis showed little overlapping among ILs, suggesting genotype-specific control of cell defense responses possibly related to differences in disease resistance strategies. Functional enrichment assessments yielded a similar pattern. However, all three ILs altered the expression of genes involved in the cellular redox state and cell wall remodeling, in agreement with current knowledge about the initiation of plant immune responses. Remarkably, the over-representation of long non-coding RNAs (lncRNA) was another common feature among ILs. Our findings highlight the diversity of transcriptional responses to SHR within sunflower breeding lines and provide evidence of lncRNAs playing a significant role at early stages of defense.

Differential Alternative Splicing Genes and Isoform Regulation Networks of Rapeseed (Brassica napus L.) Infected with Sclerotinia sclerotiorum

Alternative splicing (AS) is a post-transcriptional level of gene expression regulation that increases transcriptome and proteome diversity. How the AS landscape of rapeseed (Brassica napus L.) changes in response to the fungal pathogen Sclerotinia sclerotiorum is unknown. Here, we analyzed 18 RNA-seq libraries of mock-inoculated and S. sclerotiorum-inoculated susceptible and tolerant B. napus plants. We found that infection increased AS, with intron retention being the main AS event. To determine the key genes functioning in the AS response, we performed a differential AS (DAS) analysis. We identified 79 DAS genes, including those encoding splicing factors, defense response proteins, crucial transcription factors and enzymes. We generated coexpression networks based on the splicing isoforms, rather than the genes, to explore the genes’ diverse functions. Using this weighted gene coexpression network analysis alongside a gene ontology enrichment analysis, we identified 11 modules putatively involved in the pathogen defense response. Within these regulatory modules, six DAS genes (ascorbate peroxidase 1, ser/arg-rich protein 34a, unknown function 1138, nitrilase 2, v-atpase f, and amino acid transporter 1) were considered to encode key isoforms involved in the defense response. This study provides insight into the post-transcriptional response of B. napus to S. sclerotiorum infection.

Introgression and pyramiding of genetic loci from wild Brassica oleracea into B. napus for improving Sclerotinia resistance of rapeseed

Resistant rapeseed lines pyramided with multiple resistant QTLs derived from Brassica oleracea were developed via a hexaploidy strategy. Rapeseed (Brassica napus L.) suffers heavily from Sclerotinia stem rot, but the breeding of Sclerotinia-resistant rapeseed cultivar has been unsuccessful. During the study, interspecific hexaploids were generated between rapeseed variety 'Zhongshuang 9' and a wild B. oleracea which was highly resistant to S. sclerotiorum, followed by backcrossing with Zhongshuang 9 and successive selfing. By molecular marker-assisted selection, three major resistant QTLs were transferred and pyramided from B. oleracea into two BC1F8 lines which exhibited ~ 35% higher resistance level than Zhongshuang 9 and produced good seed yield and seed quality. It is the first report on successful development of Sclerotinia-resistant rapeseed lines by introducing multiple resistant loci from wild B. oleracea. This study revealed the effectiveness of pyramiding multiple QTLs in improving Sclerotinia resistance in rapeseed and provided a novel breeding strategy on utilization of B. oleracea in rapeseed improvement.

An update on the arsenal: mining resistance genes for disease management of Brassica crops in the genomic era

Brassica species include many economically important crops that provide nutrition and health-promoting substances to humans worldwide. However, as with all crops, their production is constantly threatened by emerging viral, bacterial, and fungal diseases, whose incidence has increased in recent years. Traditional methods of control are often costly, present limited effectiveness, and cause environmental damage; instead, the ideal approach is to mine and utilize the resistance genes of the Brassica crop hosts themselves. Fortunately, the development of genomics, molecular genetics, and biological techniques enables us to rapidly discover and apply resistance (R) genes. Herein, the R genes identified in Brassica crops are summarized, including their mapping and cloning, possible molecular mechanisms, and application in resistance breeding. Future perspectives concerning how to accurately discover additional R gene resources and efficiently utilize these genes in the genomic era are also discussed.

BnaMPK6 is a determinant of quantitative disease resistance against Sclerotinia sclerotiorum in oilseed rape

Sclerotinia sclerotiorum causes a devastating disease in oilseed rape (Brassica napus), resulting in major economic losses. Resistance response of B. napus against S. sclerotiorum exhibits a typical quantitative disease resistance (QDR) characteristic, but the molecular determinants of this QDR are largely unknown. In this study, we isolated a B. napus mitogen-activated protein kinase gene, BnaMPK6, and found that BnaMPK6 expression is highly responsive to infection by S. sclerotiorum and treatment with salicylic acid (SA) or jasmonic acid (JA). Moreover, overexpression (OE) of BnaMPK6 significantly enhances resistance to S. sclerotiorum, whereas RNAi in BnaMPK6 significantly reduces this resistance. These results showed that BnaMPK6 plays an important role in defense to S. sclerotiorum. Furthermore, expression of defense genes associated with SA-, JA- and ethylene (ET)-mediated signaling was investigated in BnaMPK6-RNAi, WT and BnaMPK6-OE plants after S. sclerotiorum infection, and consequently, it was indicated that the activation of ET signaling by BnaMPK6 may play a role in the defense. Further, four BnaMPK6-encoding homologous loci were mapped in the B. napus genome. Using the allele analysis and expression analysis on the four loci, we demonstrated that the locus BnaA03.MPK6 makes an important contribution to QDR against S. sclerotiorum. Our data indicated that BnaMPK6 is a previously unknown determinant of QDR against S. sclerotiorum in B. napus.

Identification of QTLs Containing Resistance Genes for Sclerotinia Stem Rot in Brassica napus Using Comparative Transcriptomic Studies

Sclerotinia stem rot is a major disease in Brassica napus that causes yield losses of 10-20% and reaching 80% in severely infected fields. SSR not only causes yield reduction but also causes low oil quality by reducing fatty acid content. There is a need to identify resistant genetic sources with functional significance for the breeding of SSR-resistant cultivars. In this study, we identified 17 QTLs involved in SSR resistance in three different seasons using SNP markers and disease lesion development after artificial inoculation. There were no common QTLs in all 3 years, but there were three QTLs that appeared in two seasons covering all seasons with a shared QTL. The QTLs identified in the 2 years were SRA9a, SRC2a and SRC3a with phenotypic effect variances of 14.75 and 11.57% for SRA9a, 7.49 and 10.38% for SRC3a and 7.73 and 6.81% for SRC2a in their 2 years, respectively. The flowering time was also found to have a negative correlation with disease resistance, i.e., early-maturing lines were more susceptible to disease. The stem width has shown a notably weak effect on disease development, causing researchers to ignore its effect. Given that flowering time is an important factor in disease resistance, we used comparative RNA-sequencing analysis of resistant and susceptible lines with consistent performance in 3 years with almost the same flowering time to identify the resistance genes directly involved in resistance within the QTL regions. Overall, there were more genes differentially expressed in resistant lines 19,970 than in susceptible lines 3936 compared to their mock-inoculated lines, demonstrating their tendency to cope with disease. We identified 36 putative candidate genes from the resistant lines that were upregulated in resistant lines compared to resistant mock and susceptible lines that might be involved in resistance to SSR.

Molecular and genetic analysis of defensive responses of Brassica juncea - B. fruticulosa introgression lines to Sclerotinia infection

Sclerotinia stem rot caused by Sclerotinia sclerotiorum is a major disease of crop brassicas, with inadequate variation for resistance in primary gene pools. We utilized a wild Brassicaceae species with excellent resistance against stem rot to develop a set of B. juncea - B. fruticulosa introgression lines (ILs). These were assessed for resistance using a highly reproducible stem inoculation technique against a virulent pathogen isolate. Over 40% of ILs showed higher levels of resistance. IL-43, IL-175, IL-215, IL-223 and IL-277 were most resistant ILs over three crop seasons. Sequence reads (21x) from the three most diverse ILs were then used to create B. juncea pseudomolecules, by replacing SNPs of reference B. juncea with those of re-sequenced ILs. Genotyping by sequencing (GBS) was also carried out for 88 ILs. Resultant sequence tags were then mapped on to the B. juncea pseudomolecules, and SNP genotypes prepared for each IL. Genome wide association studies helped to map resistance responses to stem rot. A total of 13 significant loci were identified on seven B. juncea chromosomes (A01, A03, A04, A05, A08, A09 and B05). Annotation of the genomic region around identified SNPs allowed identification of 20 candidate genes belonging to major disease resistance protein families, including TIR-NBS-LRR class, Chitinase, Malectin/receptor-like protein kinase, defensin-like (DEFL), desulfoglucosinolate sulfotransferase protein and lipoxygenase. A majority of the significant SNPs could be validated using whole genome sequences (21x) from five advanced generation lines being bred for Sclerotinia resistance as compared to three susceptible B. juncea germplasm lines. Our findings not only provide critical new understanding of the defensive pathway of B. fruticulosa resistance, but will also enable development of marker candidates for assisted transfer of introgressed resistant loci in to agronomically superior cultivars of crop Brassica.

Comparative Reaction of Camelina sativa to Sclerotinia sclerotiorum and Leptosphaeria maculans

Sclerotinia sclerotiorum and Leptosphaeria maculans are two of the most important pathogens of many cruciferous crops. The reaction of 30 genotypes of Camelina sativa (false flax) was determined against both pathogens. C. sativa genotypes were inoculated at seedling and adult stages with two pathotypes of S. sclerotiorum, highly virulent MBRS-1 and less virulent WW-1. There were significant differences (P < 0.001) among genotypes, between pathotypes, and a significant interaction between genotypes and pathotypes in relation to percent cotyledon disease index (% CDI) and stem lesion length. Genotypes 370 (% CDI 20.5, stem lesion length 1.8 cm) and 253 (% CDI 24.8, stem lesion length 1.4 cm) not only consistently exhibited cotyledon and stem resistance, in contrast to susceptible genotype 2305 (% CDI 37.7, stem lesion length 7.2 cm), but their resistance was independent to S. sclerotiorum pathotype. A F₅-recombinant inbred line population was developed from genotypes 370 × 2305 and responses characterized. Low broad-sense heritability indicated a complex pattern of inheritance of resistance to S. sclerotiorum. Six isolates of L. maculans, covering combinations of five different avirulent loci (i.e., five different races), were tested on C. sativa cotyledons across two experiments. There was a high level of resistance, with % CDI < 17, and including development of a hypersensitive reaction. This is the first report of variable reaction of C. sativa to different races of L. maculans and the first demonstrating comparative reactions of C. sativa to S. sclerotiorum and L. maculans. This study not only provides new understanding of these comparative resistances in C. sativa, but highlights their potential as new sources of resistance, both for crucifer disease-resistance breeding in general and to enable broader adoption of C. sativa as a more sustainable oilseed crop in its own right.

The genetics of resistance to lettuce drop (Sclerotinia spp.) in lettuce in a recombinant inbred line population from Reine des Glaces × Eruption

KEY MESSAGE

Two QTLs for resistance to lettuce drop, qLDR1.1 and qLDR5.1, were identified. Associated SNPs will be useful in breeding for lettuce drop and provide the foundation for future molecular analysis. Lettuce drop, caused by Sclerotinia minor and S. sclerotiorum, is an economically important disease of lettuce. The association of resistance to lettuce drop with the commercially undesirable trait of fast bolting has hindered the integration of host resistance in control of this disease. Eruption is a slow-bolting cultivar that exhibits a high level of resistance to lettuce drop. Eruption also is completely resistant to Verticillium wilt caused by race 1 of Verticillium dahliae. A recombinant inbred line population from the cross Reine des Glaces × Eruption was genotyped by sequencing and evaluated for lettuce drop and bolting in separate fields infested with either S. minor or V. dahliae. Two quantitative trait loci (QTLs) for lettuce drop resistance were consistently detected in at least two experiments, and two other QTLs were identified in another experiment; the alleles for resistance at all four QTLs originated from Eruption. A QTL for lettuce drop resistance on linkage group (LG) 5, qLDR5.1, was consistently detected in all experiments and explained 11 to 25% of phenotypic variation. On LG1, qLDR1.1 was detected in two experiments explaining 9 to 12% of the phenotypic variation. Three out of four resistance QTLs are distinct from QTLs for bolting; qLDR5.1 is pleiotropic or closely linked with a QTL for early bolting; however, the rate of bolting shows only a small effect on the variance in resistance observed at this locus. The SNP markers linked with these QTLs will be useful in breeding for resistance through marker-assisted selection.

CropSNPdb: a database of SNP array data for Brassica crops and hexaploid bread wheat

Advances in sequencing technology have led to a rapid rise in the genomic data available for plants, driving new insights into the evolution, domestication and improvement of crops. Single nucleotide polymorphisms (SNPs) are a major component of crop genomic diversity, and are invaluable as genetic markers in research and breeding programs. High-throughput SNP arrays, or 'SNP chips', can generate reproducible sets of informative SNP markers and have been broadly adopted. Although there are many public repositories for sequencing data, which are routinely uploaded, there are no formal repositories for crop SNP array data. To make SNP array data more easily accessible, we have developed CropSNPdb (http://snpdb.appliedbioinformatics.com.au), a database for SNP array data produced by the Illumina Infinium™ hexaploid bread wheat (Triticum aestivum) 90K and Brassica 60K arrays. We currently host SNPs from datasets covering 526 Brassica lines and 309 bread wheat lines, and provide search, download and upload utilities for users. CropSNPdb provides a useful repository for these data, which can be applied for a range of genomics and molecular crop-breeding activities.

NBS-Encoding Genes in Brassica napus Evolved Rapidly After Allopolyploidization and Co-localize With Known Disease Resistance Loci

Genes containing nucleotide-binding sites (NBS) play an important role in pathogen resistance in plants. However, the evolutionary fate of NBS-encoding genes after formation of allotetraploid Brassica napus (A_nA_nC_nC_n, 2n = 38) is still unknown. We performed a genome-wide comparison of putatively functional NBS-encoding genes in B. napus and its progenitor species Brassica rapa (A_rA_r, 2n = 20) and Brassica oleracea (C_oC_o, 2n = 18), identifying 464, 202, and 146 putatively functional NBS-encoding genes respectively, with genes unevenly distributed in several clusters. The A_n-subgenome of B. napus possessed similar numbers of NBS-encoding genes (191 genes) to the A_r genome of B. rapa (202 genes) and similar clustering patterns. However, the C_n genome of B. napus had many more genes (273) than the B. oleracea C_o genome (146), with different clustering trends. Only 97 NBS-encoding genes (66.4%) in B. oleracea were homologous with NBS-encoding genes in B. napus, while 176 NBS-encoding genes (87.1%) were homologous between B. rapa and B. napus. These results suggest a greater diversification of NBS-encoding genes in the C genome may have occurred after formation of B. napus. Although most NBS-encoding genes in B. napus appeared to derive from the progenitors, the birth and death of several NBS-encoding genes was also putatively mediated by non-homologous recombination. The Ka/Ks values of most homologous pairs between B. napus and the progenitor species were less than 1, suggesting purifying selection during B. napus evolution. The majority of NBS-encoding genes (60% in all species) showed higher expression levels in root tissue (out of root, leaf, stem, seed and flower tissue types). Comparative analysis of NBS-encoding genes with mapped resistance QTL against three major diseases of B. napus (blackleg, clubroot and Sclerotinia stem rot) found 204 NBS-encoding genes in B. napus located within 71 resistance QTL intervals. The majority of NBS-encoding genes were co-located with resistance QTLs against a single disease, while 47 genes were co-located with QTLs against two diseases and 3 genes were co-located with QTLs against all three. Our results revealed significant variation as well as interesting evolutionary trajectories of NBS-encoding genes in the different Brassica subgenomes, while co-localization of NBS-encoding genes and resistance QTL may facilitate resistance breeding in oilseed rape.

Recent Advances in Mechanisms of Plant Defense to Sclerotinia sclerotiorum

Sclerotinia sclerotiorum (Lib.) de Bary is an unusual pathogen which has the broad host range, diverse infection modes, and potential double feeding lifestyles of both biotroph and necrotroph. It is capable of infecting over 400 plant species found worldwide and more than 60 names have agriculturally been used to refer to diseases caused by this pathogen. Plant defense to S. sclerotiorum is a complex biological process and exhibits a typical quantitative disease resistance (QDR) response. Recent studies using Arabidopsis thaliana and crop plants have obtained new advances in mechanisms used by plants to cope with S. sclerotiorum infection. In this review, we focused on our current understanding on plant defense mechanisms against this pathogen, and set up a model for the defense process including three stages: recognition of this pathogen, signal transduction and defense response. We also have a particular interest in defense signaling mediated by diverse signaling molecules. We highlight the current challenges and unanswered questions in both the defense process and defense signaling. Essentially, we discussed candidate resistance genes newly mapped by using high-throughput experiments in important crops, and classified these potential gene targets into different stages of the defense process, which will broaden our understanding of the genetic architecture underlying quantitative resistance to S. sclerotiorum. We proposed that more powerful mapping population(s) will be required for accurate and reliable QDR gene identification.

Genome-wide association study identifies favorable SNP alleles and candidate genes for waterlogging tolerance in chrysanthemums

Chrysanthemums are sensitive to waterlogging stress, and the development of screening methods for tolerant germplasms or genes and the breeding of tolerant new varieties are of great importance in chrysanthemum breeding. To understand the genetic basis of waterlogging tolerance (WT) in chrysanthemums, we performed a genome-wide association study (GWAS) using 92,811 single nucleotide polymorphisms (SNPs) in a panel of 88 chrysanthemum accessions, including 64 spray cut and 24 disbud chrysanthemums. The results showed that the average MFVW (membership function value of waterlogging) of the disbud type (0.65) was significantly higher than that of the spray type (0.55) at P < 0.05, and the MFVW of the Asian accessions (0.65) was significantly higher than that of the European accessions (0.48) at P < 0.01. The GWAS performed using the general linear model (GLM) and mixed linear model (MLM) identified 137 and 14 SNP loci related to WT, respectively, and 11 associations were commonly predicted. By calculating the phenotypic effect values for 11 common SNP loci, six highly favorable SNP alleles that explained 12.85-21.85% of the phenotypic variations were identified. Furthermore, the dosage-pyramiding effects of the favorable alleles and the significant linear correlations between the numbers of highly favorable alleles and phenotypic values were identified (r ² = 0.45; P < 0.01). A major SNP locus (Marker6619-75) was converted into a derived cleaved amplified polymorphic sequence (dCAPS) marker that cosegregated with WT with an average efficiency of 78.9%. Finally, four putative candidate genes in the WT were identified via quantitative real-time PCR (qRT-PCR). The results presented in this study provide insights for further research on WT mechanisms and the application of molecular marker-assisted selection (MAS) in chrysanthemum WT breeding programs.

Genome-wide association study identifies favorable SNP alleles and candidate genes for waterlogging tolerance in chrysanthemums

Chrysanthemums are sensitive to waterlogging stress, and the development of screening methods for tolerant germplasms or genes and the breeding of tolerant new varieties are of great importance in chrysanthemum breeding. To understand the genetic basis of waterlogging tolerance (WT) in chrysanthemums, we performed a genome-wide association study (GWAS) using 92,811 single nucleotide polymorphisms (SNPs) in a panel of 88 chrysanthemum accessions, including 64 spray cut and 24 disbud chrysanthemums. The results showed that the average MFVW (membership function value of waterlogging) of the disbud type (0.65) was significantly higher than that of the spray type (0.55) at P < 0.05, and the MFVW of the Asian accessions (0.65) was significantly higher than that of the European accessions (0.48) at P < 0.01. The GWAS performed using the general linear model (GLM) and mixed linear model (MLM) identified 137 and 14 SNP loci related to WT, respectively, and 11 associations were commonly predicted. By calculating the phenotypic effect values for 11 common SNP loci, six highly favorable SNP alleles that explained 12.85-21.85% of the phenotypic variations were identified. Furthermore, the dosage-pyramiding effects of the favorable alleles and the significant linear correlations between the numbers of highly favorable alleles and phenotypic values were identified (r ² = 0.45; P < 0.01). A major SNP locus (Marker6619-75) was converted into a derived cleaved amplified polymorphic sequence (dCAPS) marker that cosegregated with WT with an average efficiency of 78.9%. Finally, four putative candidate genes in the WT were identified via quantitative real-time PCR (qRT-PCR). The results presented in this study provide insights for further research on WT mechanisms and the application of molecular marker-assisted selection (MAS) in chrysanthemum WT breeding programs.

CRISPR/Cas9-Mediated Multiplex Genome Editing of the BnWRKY11 and BnWRKY70 Genes in Brassica napus L

Targeted genome editing is a desirable means of basic science and crop improvement. The clustered, regularly interspaced, palindromic repeat (CRISPR)/Cas9 (CRISPR-associated 9) system is currently the simplest and most commonly used system in targeted genomic editing in plants. Single and multiplex genome editing in plants can be achieved under this system. In Arabidopsis, AtWRKY11 and AtWRKY70 genes were involved in JA- and SA-induced resistance to pathogens, in rapeseed (Brassica napus L.), BnWRKY11 and BnWRKY70 genes were found to be differently expressed after inoculated with the pathogenic fungus, Sclerotinia sclerotiorum (Lib.) de Bary. In this study, two Cas9/sgRNA constructs targeting two copies of BnWRKY11 and four copies of BnWRKY70 were designed to generate BnWRKY11 and BnWRKY70 mutants respectively. As a result, twenty-two BnWRKY11 and eight BnWRKY70 independent transformants (T₀) were obtained, with the mutation ratios of 54.5% (12/22) and 50% (4/8) in BnWRKY11 and BnWRKY70 transformants respectively. Eight and two plants with two copies of mutated BnWRKY11 and BnWRKY70 were obtained respectively. In T₁ generation of each plant examined, new mutations on target genes were detected with high efficiency. The vast majority of BnWRKY70 mutants showed editing in three copies of BnWRKY70 in examined T₁ plants. BnWRKY70 mutants exhibited enhanced resistance to Sclerotinia, while BnWRKY11 mutants showed no significant difference in Sclerotinia resistance when compared to non-transgenic plants. In addition, plants that overexpressed BnWRKY70 showed increased sensitivity when compared to non-transgenic plants. Altogether, our results demonstrated that BnWRKY70 may function as a regulating factor to negatively control the Sclerotinia resistance and CRISPR/Cas9 system could be used to generate germplasm in B. napus with high resistance against Sclerotinia.

Genomic evidence for genes encoding leucine-rich repeat receptors linked to resistance against the eukaryotic extra- and intracellular Brassica napus pathogens Leptosphaeria maculans and Plasmodiophora brassicae

Genes coding for nucleotide-binding leucine-rich repeat (LRR) receptors (NLRs) control resistance against intracellular (cell-penetrating) pathogens. However, evidence for a role of genes coding for proteins with LRR domains in resistance against extracellular (apoplastic) fungal pathogens is limited. Here, the distribution of genes coding for proteins with eLRR domains but lacking kinase domains was determined for the Brassica napus genome. Predictions of signal peptide and transmembrane regions divided these genes into 184 coding for receptor-like proteins (RLPs) and 121 coding for secreted proteins (SPs). Together with previously annotated NLRs, a total of 720 LRR genes were found. Leptosphaeria maculans-induced expression during a compatible interaction with cultivar Topas differed between RLP, SP and NLR gene families; NLR genes were induced relatively late, during the necrotrophic phase of pathogen colonization. Seven RLP, one SP and two NLR genes were found in Rlm1 and Rlm3/Rlm4/Rlm7/Rlm9 loci for resistance against L. maculans on chromosome A07 of B. napus. One NLR gene at the Rlm9 locus was positively selected, as was the RLP gene on chromosome A10 with LepR3 and Rlm2 alleles conferring resistance against L. maculans races with corresponding effectors AvrLm1 and AvrLm2, respectively. Known loci for resistance against L. maculans (extracellular hemi-biotrophic fungus), Sclerotinia sclerotiorum (necrotrophic fungus) and Plasmodiophora brassicae (intracellular, obligate biotrophic protist) were examined for presence of RLPs, SPs and NLRs in these regions. Whereas loci for resistance against P. brassicae were enriched for NLRs, no such signature was observed for the other pathogens. These findings demonstrate involvement of (i) NLR genes in resistance against the intracellular pathogen P. brassicae and a putative NLR gene in Rlm9-mediated resistance against the extracellular pathogen L. maculans.

Exploring the genetics of lesion and nodal resistance in pea (Pisum sativum L.) to Sclerotinia sclerotiorum using genome-wide association studies and RNA-Seq

The disease white mold caused by the fungus Sclerotinia sclerotiorum is a significant threat to pea production, and improved resistance to this disease is needed. Nodal resistance in plants is a phenomenon where a fungal infection is prevented from passing through a node, and the infection is limited to an internode region. Nodal resistance has been observed in some pathosystems such as the pea (Pisum sativum L.)-S. sclerotiorum pathosystem. In addition to nodal resistance, different pea lines display different levels of stem lesion size restriction, referred to as lesion resistance. It is unclear whether the genetics of lesion resistance and nodal resistance are identical or different. This study applied genome-wide association studies (GWAS) and RNA-Seq to understand the genetic makeup of these two types of resistance. The time series RNA-Seq experiment consisted of two pea lines (the susceptible 'Lifter' and the partially resistant PI 240515), two treatments (mock inoculated samples and S. sclerotiorum-inoculated samples), and three time points (12, 24, and 48 hr post inoculation). Integrated results from GWAS and RNA-Seq analyses identified different redox-related transcripts for lesion and nodal resistances. A transcript encoding a glutathione S-transferase was the only shared resistance variant for both phenotypes. There were more leucine rich-repeat containing transcripts found for lesion resistance, while different candidate resistance transcripts such as a VQ motif-containing protein and a myo-inositol oxygenase were found for nodal resistance. This study demonstrated the robustness of combining GWAS and RNA-Seq for identifying white mold resistance in pea, and results suggest different genetics underlying lesion and nodal resistance.

Development and Applications of a High Throughput Genotyping Tool for Polyploid Crops: Single Nucleotide Polymorphism (SNP) Array

Polypoid species play significant roles in agriculture and food production. Many crop species are polyploid, such as potato, wheat, strawberry, and sugarcane. Genotyping has been a daunting task for genetic studies of polyploid crops, which lags far behind the diploid crop species. Single nucleotide polymorphism (SNP) array is considered to be one of, high-throughput, relatively cost-efficient and automated genotyping approaches. However, there are significant challenges for SNP identification in complex, polyploid genomes, which has seriously slowed SNP discovery and array development in polyploid species. Ploidy is a significant factor impacting SNP qualities and validation rates of SNP markers in SNP arrays, which has been proven to be a very important tool for genetic studies and molecular breeding. In this review, we (1) discussed the pros and cons of SNP array in general for high throughput genotyping, (2) presented the challenges of and solutions to SNP calling in polyploid species, (3) summarized the SNP selection criteria and considerations of SNP array design for polyploid species, (4) illustrated SNP array applications in several different polyploid crop species, then (5) discussed challenges, available software, and their accuracy comparisons for genotype calling based on SNP array data in polyploids, and finally (6) provided a series of SNP array design and genotype calling recommendations. This review presents a complete overview of SNP array development and applications in polypoid crops, which will benefit the research in molecular breeding and genetics of crops with complex genomes.

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.