Syntenic quantitative trait loci and genomic divergence for Sclerotinia resistance and flowering time in Brassica 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.

Functional and evolutionary study of MLO gene family in the regulation of Sclerotinia stem rot resistance in Brassica napus L

BACKGROUND

Oilseed rape (Brassica napus L.) is known as one of the most important oilseed crops cultivated around the world. However, its production continuously faces a huge challenge of Sclerotinia stem rot (SSR), a destructive disease caused by the fungus Sclerotinia sclerotiorum, resulting in huge yield loss annually. The SSR resistance in B. napus is quantitative and controlled by a set of minor genes. Identification of these genes and pyramiding them into a variety are a major strategy for SSR resistance breeding in B. napus.

RESULTS

Here, we performed a genome-wide association study (GWAS) using a natural population of B. napus consisting of 222 accessions to identify BnaA08g25340D (BnMLO2_2) as a candidate gene that regulates the SSR resistance. BnMLO2_2 was a member of seven homolog genes of Arabidopsis Mildew Locus O 2 (MLO2) and the significantly SNPs were mainly distributed in the promoter of BnMLO2_2, suggesting a role of BnMLO2_2 expression level in the regulation of SSR resistance. We expressed BnMLO2_2 in Arabidopsis and the transgenic plants displayed an enhanced SSR resistance. Transcriptome profiling of different tissues of B. napus revealed that BnMLO2_2 had the most expression level in leaf and silique tissues among all the 7 BnMLO2 members and also expressed higher in the SSR resistant accession than in the susceptible accession. In Arabidopsis, mlo2 plants displayed reduced resistance to SSR, whereas overexpression of MLO2 conferred plants an enhanced SSR resistance. Moreover, a higher expression level of MLO2 showed a stronger SSR resistance in the transgenic plants. The regulation of MLO2 in SSR resistance may be associated with the cell death. Collinearity and phylogenetic analysis revealed a large expansion of MLO family in Brassica crops.

CONCLUSION

Our study revealed an important role of BnMLO2 in the regulation of SSR resistance and provided a new gene candidate for future improvement of SSR resistance in B. napus and also new insights into understanding of MLO family evolution in Brassica crops.

Integrating GWAS, linkage mapping and gene expression analyses reveal the genetic control of first branch height in Brassica napus L

Rapeseed (Brassica napus L.) is a crucial oil crop cultivated worldwide. First branch height, an essential component of rapeseed plant architecture, has an important effect on yield and mechanized harvesting; however, the underlying genetic mechanism remains unclear. In this study, based on the 60K single nucleotide polymorphism array and a recombinant inbred lines population derived from M083 and 888-5, a total of 19 QTLs were detected in five environments, distributed on linkage groups A02, A09, A10, C06, and C07, which explained phenotypic variation ranging from 4.87 to 29.87%. Furthermore, 26 significant SNPs were discovered on Chr.A02 by genome-wide association study in a diversity panel of 324 re-sequencing accessions. The major QTL of the first branch height trait was co-located on Chr.A02 by integrating linkage mapping and association mapping. Eleven candidate genes were screened via allelic variation analysis, inter-subgenomic synteny analysis, and differential expression of genes in parental shoot apical meristem tissues. Among these genes, BnaA02g13010D, which encodes a TCP transcription factor, was confirmed as the target gene according to gene function annotation, haplotype analysis, and full-length gene sequencing, which revealed that TATA insertion/deletion in the promoter region was closely linked to significantly phenotypic differences BnaA02.TCP1 ^M083 overexpression resulted in decreased branch height and increased branch number in Arabidopsis. These results provide a genetic basis for first branch height and the ideal architecture of B. napus.

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.

Multienvironment QTL analysis delineates a major locus associated with homoeologous exchanges for water-use efficiency and seed yield in canola

Canola varieties exhibit variation in drought avoidance and drought escape traits, reflecting adaptation to water-deficit environments. Our understanding of underlying genes and their interaction across environments in improving crop productivity is limited. A doubled haploid population was analysed to identify quantitative trait loci (QTL) associated with water-use efficiency (WUE) related traits. High WUE in the vegetative phase was associated with low seed yield. Based on the resequenced parental genome data, we developed sequence-capture-based markers and validated their linkage with carbon isotope discrimination (Δ¹³ C) in an F₂ population. RNA sequencing was performed to determine the expression of candidate genes underlying Δ¹³ C QTL. QTL contributing to main and QTL × environment interaction effects for Δ¹³ C and yield were identified. One multiple-trait QTL for Δ¹³ C, days to flower, plant height, and seed yield was identified on chromosome A09. Interestingly, this QTL region overlapped with a homoeologous exchange (HE) event, suggesting its association with the multiple traits. Transcriptome analysis revealed 121 significantly differentially expressed genes underlying Δ¹³ C QTL on A09 and C09, including in HE regions. Sorting out the negative relationship between vegetative WUE and seed yield is a priority. Genetic and genomic resources and knowledge so developed could improve canola WUE and yield.

Genome-Wide Association Analysis Combined With Quantitative Trait Loci Mapping and Dynamic Transcriptome Unveil the Genetic Control of Seed Oil Content in Brassica napus L

Rapeseed, an allotetraploid oil crop, provides vegetable oil for human consumption. The growing demand for oilseeds has necessitated the development of rapeseed varieties with improved quality. Therefore, a clear understanding of the genetic basis underlying the seed oil content (SOC) is required. In this study, a natural population comprising 204 diverse accessions and recombinant inbred lines (RILs) derived from Brassica napus and Sinapis alba via distant hybridization were collected for genome-wide association analysis (GWAS) and quantitative trait loci (QTL) mapping of the SOC trait, respectively. The variable coefficient of the RIL and natural populations ranged from 7.43 to 10.43% and 8.40 to 10.91%. Then, a high-density linkage map was constructed based on whole genome re-sequencing (WGS); the map harbored 2,799 bin markers and covered a total distance of 1,835.21 cM, with an average marker interval of 0.66 cM. The QTLs for SOC on chromosome A07 were stably detected in both single and multiple environments. Finally, a novel locus qA07.SOC was identified as the major QTL for SOC based on the GWAS and RIL populations. In addition, the RNA-seq results showed that photosynthesis, lipid biosynthesis proteins, fatty acid metabolism, and unsaturated fatty acid biosynthesis were significantly different between the developed seeds of the two parents of the RIL population. By comparing the variation information and expression levels of the syntenic genes within qA07.SOC and its syntenic genomic regions, as well as through haplotype analysis via GWAS, BnaA07.STR18, BnaA07.NRT1, and BnaA07g12880D were predicted as candidate genes in the qA07.SOC interval. These stable QTLs containing candidate genes and haplotypes can potentially provide a reliable basis for marker-assisted selection in B. napus breeding for SOC.

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.

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.

BnaA03.MKK5-BnaA06.MPK3/BnaC03.MPK3 Module Positively Contributes to Sclerotinia sclerotiorum Resistance in Brassica napus

Brassica napus (oilseed rape) is one of the most important oil crops worldwide, but its growth is seriously threatened by Sclerotinia sclerotiorum. The mechanism of oilseed rape response to this pathogen has rarely been studied. Here, it was identified that BnaA03.MKK5 whose expression was induced by S. sclerotiorum infection was involved in plant immunity. BnaA03.MKK5 overexpression lines exhibited decreased disease symptoms compared to wild-type plants, accompanied by the increased expression of camalexin-biosynthesis-related genes, including BnPAD3 and BnCYP71A13. In addition, two copies of BnMPK3 (BnA06.MPK3 and BnC03.MPK3) were induced by Sclerotinia incubation, and BnaA03.MKK5 interacted with BnaA06.MPK3/BnaC03.MPK3 in yeast. These interactions were confirmed using in vivo co-immunoprecipitation assays. In vitro phosphorylation assays showed that BnaA06.MPK3 and BnaC03.MPK3 were the direct phosphorylation substrates of BnaA03.MKK5. The transgenic oilseed rape plants including BnaA06.MPK3 and BnaC03.MPK3 overexpression lines and BnMPK3 gene editing lines mediated by CRISPR/Cas9 were generated; the results of the genetic transformation of BnaA06.MPK3/BnaC03.MPK3 indicate that BnMPK3 also has a positive role in Sclerotinia resistance. This study provides information about the potential mechanism of B. napus defense against S. Sclerotiorum mediated by a detailed BnaA03.MKK5-BnaA06.MPK3/BnaC03.MPK3 module.

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.

Different Shades of Kale-Approaches to Analyze Kale Variety Interrelations

Brassica oleracea is a vegetable crop with an amazing morphological diversity. Among the various crops derived from B. oleracea, kale has been in the spotlight globally due to its various health-benefitting compounds and many different varieties. Knowledge of the existing genetic diversity is essential for the improved breeding of kale. Here, we analyze the interrelationships, population structures, and genetic diversity of 72 kale and cabbage varieties by extending our previous diversity analysis and evaluating the use of summed potential lengths of shared haplotypes (SPLoSH) as a new method for such analyses. To this end, we made use of the high-density Brassica 60K SNP array, analyzed SNPs included in an available Brassica genetic map, and used these resources to generate and evaluate the information from SPLoSH data. With our results we could consistently differentiate four groups of kale across all analyses: the curly kale varieties, Italian, American, and Russian varieties, as well as wild and cultivated types. The best results were achieved by using SPLoSH information, thus validating the use of this information in improving analyses of interrelations in kale. In conclusion, our definition of kale includes the curly varieties as the kales in a strict sense, regardless of their origin. These results contribute to a better understanding of the huge diversity of kale and its interrelations.

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.

A Novel Effector Protein SsERP1 Inhibits Plant Ethylene Signaling to Promote Sclerotinia sclerotiorum Infection

Sclerotinia sclerotiorum is one of the most devastating pathogens in Brassica napus and causes huge economic loss worldwide. Though around one hundred putative effectors have been predicted in Sclerotinia sclerotiorum genome, their functions are largely unknown. In this study, we cloned and characterized a novel effector, SsERP1 (ethylene pathway repressor protein 1), in Sclerotinia sclerotiorum. SsERP1 is a secretory protein highly expressed at the early stages of Sclerotinia sclerotiorum infection. Ectopic overexpression of SsERP1 in plant leaves promoted Sclerotinia sclerotiorum infection, and the knockout mutants of SsERP1 showed reduced pathogenicity but retained normal mycelial growth and sclerotium formation, suggesting that SsERP1 specifically contributes to the pathogenesis of Sclerotinia sclerotiorum. Transcriptome analysis indicated that SsERP1 promotes Sclerotinia sclerotiorum infection by inhibiting plant ethylene signaling pathway. Moreover, we showed that knocking down SsERP1 by in vitro synthesized double-strand RNAs was able to effectively inhibit Sclerotinia sclerotiorum infection, which verifies the function of SsERP1 in Sclerotinia sclerotiorum pathogenesis and further suggests a potential strategy for Sclerotinia disease control.

Mapping of a major QTL controlling plant height using a high-density genetic map and QTL-seq methods based on whole-genome resequencing in Brassica napus

Plant height is a crucial element related to plant architecture that influences the seed yield of oilseed rape (Brassica napus L.). In this study, we isolated a natural B. napus mutant, namely a semi-dwarf mutant (sdw-e), which exhibits a 30% reduction in plant height compared with Zhongshuang 11-HP (ZS11-HP). Quantitative trait locus sequencing (QTL-seq) was conducted using two extreme DNA bulks in F₂ populations in Wuchang-2017 derived from ZS11-HP × sdw-e to identify QTLs associated with plant height. The result suggested that two QTL intervals were located on chromosome A10. The F₂ population consisting of 200 individuals in Yangluo-2018 derived from ZS11-HP × sdw-e was used to construct a high-density linkage map using whole-genome resequencing. The high-density linkage map harbored 4323 bin markers and covered a total distance of 2026.52 cM with an average marker interval of 0.47 cM. The major QTL for plant height named qPHA10 was identified on linkage group A10 by interval mapping and composite interval mapping methods. The major QTL qPHA10 was highly consistent with the QTL-seq results. And then, we integrated the variation sites and expression levels of genes in the major QTL interval to predict the candidate genes. Thus, the identified QTL and candidate genes could be used in marker-assisted selection for B. napus breeding in the future.

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.

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.

BnaMPK3 Is a Key Regulator of Defense Responses to the Devastating Plant Pathogen Sclerotinia sclerotiorum in Oilseed Rape

The disease caused by Sclerotinia sclerotiorum has traditionally been difficult to control, resulting in tremendous economic losses in oilseed rape (Brassica napus). Identification of important genes in the defense responses is critical for molecular breeding, an important strategy for controlling the disease. Here, we report that a B. napus mitogen-activated protein kinase gene, BnaMPK3, plays an important role in the defense against S. sclerotiorum in oilseed rape. BnaMPK3 is highly expressed in the stems, flowers and leaves, and its product is localized in the nucleus. Furthermore, BnaMPK3 is highly responsive to infection by S. sclerotiorum and treatment with jasmonic acid (JA) or the biosynthesis precursor of ethylene (ET), but not to treatment with salicylic acid (SA) or abscisic acid. Moreover, overexpression (OE) of BnaMPK3 in B. napus and Nicotiana benthamiana results in significantly enhanced resistance to S. sclerotiorum, whereas resistance is diminished in RNAi transgenic plants. After S. sclerotiorum infection, defense responses associated with ET, JA, and SA signaling are intensified in the BnaMPK3-OE plants but weakened in the BnaMPK3-RNAi plants when compared to those in the wild type plants; by contrast the level of both H₂O₂ accumulation and cell death exhibits a reverse pattern. The candidate gene association analyses show that the BnaMPK3-encoding BnaA06g18440D locus is a cause of variation in the resistance to S. sclerotiorum in natural B. napus population. These results suggest that BnaMPK3 is a key regulator of multiple defense responses to S. sclerotiorum, which may guide the resistance improvement of oilseed rape and related economic crops.