Multi-year linkage and association mapping confirm the high number of genomic regions involved in oilseed rape quantitative resistance to blackleg

Using genotype imputation to integrate Canola populations for genome-wide association and genomic prediction of blackleg resistance

BACKGROUND

Integrating germplasm populations genotyped by different genotyping platforms via genotype imputation is a way to utilize accumulated genetic resources. In this study, we used 278 canola samples genotyped via whole-genome sequencing (WGS) at 10× coverage to evaluate the imputation accuracy of three imputation approaches. The optimal imputation methods were used to impute and integrate two Canola genotype datasets: a diverse canola collection genotyped by genotyping-by-sequencing via transcriptome (GBS-t) and a double haploid (DH) line collection genotyped with low-coverage WGS (skim-WGS). The genomic predictive ability (GP) and detection power of marker‒trait association (GWAS) of the combined population for blackleg resistance were evaluated.

RESULTS

The empirical imputation accuracy (r2) measured as the squared correlation between observed and imputed genotypes was moderate for Minimac3 when imputing from the GBS-t density to the WGS. The accuracy dramatically improved from 0.64 to 0.82 by removing SNPs with poor Minimac3-reported Rsq (Rsq < 0.2) quality statistics. The r2 for GLIMPSE was higher than that for Beagle when imputing from different low-coverage to full-coverage WGS. We imputed and integrated the diverse canola collection and the DH lines, and the combined population showed similar or slightly greater predictive ability (PA) for blackleg resistance traits than did each of the single populations with ~ 921 K SNPs. Higher marker-trait association (MTA) detection powers were indicated with the combined population; however, similar numbers of MTAs were discovered when each single population was combined in a meta-GWAS.

CONCLUSION

It is feasible to impute and integrate germplasms from different sequencing platforms for downstream analyses. However, genetic heterogeneity across populations could add complexity to the analysis. Increasing the sample size by combining datasets showed slightly greater predictive ability and greater detection power in GWASs in the present study.

Quantifying Genetic Parameters for Blackleg Resistance in Rapeseed: A Comparative Study

Selection is a fundamental part of the plant breeding process, enabling the identification and development of varieties with desirable traits. Thanks to advances in genetics and biotechnology, the selection process has become more precise and efficient, resulting in faster breeding progress and better adaptation of crops to environmental challenges. Genetic parameters related to gene additivity and epistasis play a key role and can influence decisions on the suitability of breeding material. In this study, 188 rapeseed doubled haploid lines were assessed in field conditions for resistance to Leptosphaeria spp. Through next-generation sequencing, a total of 133,764 molecular markers (96,121 SilicoDArT and 37,643 SNP) were obtained. The similarity of the DH lines at the phenotypic and genetic levels was calculated. The results indicate that the similarity at the phenotypic level was markedly different from the similarity at the genetic level. Genetic parameters related to additive gene action effects and epistasis (double and triple) were calculated using two methods: based on phenotypic observations only and using molecular marker observations. All evaluated genetic parameters (additive, additive-additive and additive-additive-additive) were statistically significant for both estimation methods. The parameters associated with the interaction (double and triple) had opposite signs depending on the estimation method.

DArTseq-Based, High-Throughput Identification of Novel Molecular Markers for the Detection of Blackleg (Leptosphaeria Spp.) Resistance in Rapeseed

Blackleg disease, caused by Leptosphaeria spp. fungi, is one of the most important diseases of Brassica napus, responsible for severe yield losses worldwide. Blackleg resistance is controlled by major R genes and minor quantitative trait loci (QTL). Due to the high adaptation ability of the pathogen, R-mediated resistance can be easily broken, while the resistance mediated via QTL is believed to be more durable. Thus, the identification of novel molecular markers linked to blackleg resistance for B. napus breeding programs is essential. In this study, 183 doubled haploid (DH) rapeseed lines were assessed in field conditions for resistance to Leptosphaeria spp. Subsequently, DArTseq-based Genome-Wide Association Study (GWAS) was performed to identify molecular markers linked to blackleg resistance. A total of 133,764 markers (96,121 SilicoDArT and 37,643 SNP) were obtained. Finally, nine SilicoDArT and six SNP molecular markers were associated with plant resistance to Leptosphaeria spp. at the highest significance level, p < 0.001. Importantly, eleven of these fifteen markers were found within ten genes located on chromosomes A06, A07, A08, C02, C03, C06 and C08. Given the immune-related functions of the orthologues of these genes in Arabidopsis thaliana, the identified markers hold great promise for application in rapeseed breeding programs.

Large-scale population survey of Leptosphaeria maculans in France highlights both on-going breakdowns and potentially effective resistance genes in oilseed rape

BACKGROUND

Leptosphaeria maculans, the cause of stem canker of oilseed rape, develops gene-for-gene interactions with its host and shows a high evolutionary potential to 'break down' novel resistance genes (R, Rlm) deployed in cultivars over large areas. For optimal management of R genes, updated knowledge of the population structure of the pathogen is needed. In France, large-scale surveys have been done at 10-year intervals since 2000. Here we report the characterization of a large L. maculans population collected in France in 2019-2020.

RESULTS

A total of 844 isolates were collected from 11 sites in ten French departments and were phenotyped for their virulence against nine Brassica napus R genes. All isolates were virulent toward Rlm2 and Rlm9. Very few isolates were avirulent on Rlm1 (1.8%) and Rlm4 (0.6%). Avirulent isolates toward Rlm7 ('AvrLm7') varied from 67% to 11.3%, depending on the site sampled, illustrating the ongoing breakdown of Rlm7. The decrease of AvrLm7 isolates (29.2% at the national level) compared to the 2010 survey (96.5%) was accompanied by an increase of avirulent isolates on Rlm3 (0% in 2010; 54% in 2019-2020). However, virulent isolates on both Rlm3 and Rlm7, previously rarely detected, were found in all sites with a frequency of 17.3%. Finally, most or all isolates were avirulent on Rlm11 (96.1%), LepR2 (RlmS, 99.8%), and Rlm6 (100%), suggesting these three genes still effectively control the disease.

CONCLUSION

These data will help guide strategies for breeding and deploying resistant oilseed rape varieties against L. maculans in France. © 2023 The Authors. Pest Management Science published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

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.

Brassica napus Haploid and Double Haploid Production and Its Latest Applications

Rapeseed is one of the most important oil crops in the world. Increasing demand for oil and limited agronomic capabilities of present-day rapeseed result in the need for rapid development of new, superior cultivars. Double haploid (DH) technology is a fast and convenient approach in plant breeding as well as genetic research. Brassica napus is considered a model species for DH production based on microspore embryogenesis; however, the molecular mechanisms underlying microspore reprogramming are still vague. It is known that morphological changes are accompanied by gene and protein expression patterns, alongside carbohydrate and lipid metabolism. Novel, more efficient methods for DH rapeseed production have been reported. This review covers new findings and advances in Brassica napus DH production as well as the latest reports related to agronomically important traits in molecular studies employing the double haploid rapeseed lines.

Molecular Interactions Between Leptosphaeria maculans and Brassica Species

Canola is an important oilseed crop, providing food, feed, and fuel around the world. However, blackleg disease, caused by the ascomycete Leptosphaeria maculans, causes significant yield losses annually. With the recent advances in genomic technologies, the understanding of the Brassica napus-L. maculans interaction has rapidly increased, with numerous Avr and R genes cloned, setting this system up as a model organism for studying plant-pathogen associations. Although the B. napus-L. maculans interaction follows Flor's gene-for-gene hypothesis for qualitative resistance, it also puts some unique spins on the interaction. This review discusses the current status of the host-pathogen interaction and highlights some of the future gaps that need addressing moving forward.

"Late" effectors from Leptosphaeria maculans as tools for identifying novel sources of resistance in Brassica napus

The Dothideomycete Leptosphaeria maculans, causing stem canker (blackleg) of Brassica napus, secretes different cocktails of effectors at specific infection stages. Some effectors ("Late" effectors) are specifically produced during the long asymptomatic phase of stem colonization. By manipulating their expression so that they are overexpressed during cotyledon infection (OEC transformants of the fungus), we previously postulated that resistance genes operating in the stem may be involved in gene-for-gene relationship and thus contribute to quantitative disease resistance (QDR). Here, we selected 10 relevant new effector genes, and we generated OEC transformants to screen a collection of 130 B. napus genotypes, representative of the available diversity in the species. Five B. napus accessions showed a typical hypersensitive response when challenged with effectors LmSTEE98 or LmSTEE6826 at the cotyledon stage, and all belong to the semi-winter type of the diversity panel. In addition, five winter-type genotypes displayed an intermediate response to another late effector, LmSTEE7919. These new interactions now have to be genetically validated to check that they also correspond to gene-for-gene interactions. In all cases, they potentially provide novel resources, easy to breed for, and accounting for part of the quantitative resistance in a species for which we are currently facing limited resistance sources.

Genetic dissection of heterotic loci associated with plant weight by Graded pool-seq in heading Chinese cabbage (Brassica rapa)

Four heterotic QTL and a heterozygous segment for plant weight were identified by Graded Pool-Seq, QTL-seq and traditional genetic linkage analysis in heading Chinese cabbage. Heading Chinese cabbage (Brassica rapa L. spp. pekinensis) is a cross-pollinated leafy vegetable with significant heterosis. The use of heterosis is important for breeding high-yield Chinese cabbage hybrids. However, the formation and mechanism of heterosis have not been studied. We dissected the molecular mechanism of heterosis of yield-related traits in Chinese cabbage. An F₁ hybrid with high-parent heterosis of yield-related traits was selected and self-pollinated to generate segregating F₂ populations. QTL-seq, Graded Pool-seq (GPS), and traditional genetic linkage analysis were used to identify four heterotic quantitative trait loci (QTL) for plant weight: qPW1.1, qPW5.1, qPW7.1, and qPW8.1. Traditional genetic linkage analysis over two years showed that qPW8.1, located in marker A08_S45 (18,172,719) and A08_S85 (18,196,752), was mapped to a 23.5 kb genomic region. QTL qPW8.1 explained 8.6% and 23.6% of the phenotypic variation in plant weight and the total numbers of head leaves, respectively, and contained a heterozygous segment that might control the heterosis of plant weight. The qPW1.1 made an 11.7% phenotypic contribution to plant weight. The qPW7.1 was sensitive to environmental influence and explained 10.7% of the phenotypic variance. QTL qPW5.1 had a significant signal and was located in a genetic region near the centromere showing high heterozygosity. The "pseudo-overdominance" and "synergistic allelic" effects from parent line "XJD4" appear to play an important role in heterosis for plant weight in Chinese cabbage. These results provide a basis for an improved understanding of the molecular mechanism of yield-related traits and their heterosis.

Epigenomic and structural events preclude recombination in Brassica napus

Meiotic recombination is a major evolutionary process generating genetic diversity at each generation in sexual organisms. However, this process is highly regulated, with the majority of crossovers lying in the distal chromosomal regions that harbor low DNA methylation levels. Even in these regions, some islands without recombination remain, for which we investigated the underlying causes. Genetic maps were established in two Brassica napus hybrids to detect the presence of such large nonrecombinant islands. The role played by DNA methylation and structural variations in this local absence of recombination was determined by performing bisulfite sequencing and whole genome comparisons. Inferred structural variations were validated using either optical mapping or oligo fluorescence in situ hybridization. Hypermethylated or inverted regions between Brassica genomes were associated with the absence of recombination. Pairwise comparisons of nine B. napus genome assemblies revealed that such inversions occur frequently and may contain key agronomic genes such as resistance to biotic stresses. We conclude that such islands without recombination can have different origins, such as DNA methylation or structural variations in B. napus. It is thus essential to take into account these features in breeding programs as they may hamper the efficient combination of favorable alleles in elite varieties.

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.

Metabolic Capacity Differentiates Plenodomus lingam from P. biglobosus Subclade 'brassicae', the Causal Agents of Phoma Leaf Spotting and Stem Canker of Oilseed Rape (Brassica napus) in Agricultural Ecosystems

In contrast to the long-lasting taxonomic classification of Plenodomus lingam and P. biglobosus as one species, formerly termed Leptosphaeria maculans, both species form separate monophyletic groups, comprising sub-classes, differing considerably with epidemiology towards Brassicaceae plants. Considering the great differences between P. lingam and P. biglobosus, we hypothesized their metabolic capacities vary to a great extent. The experiment was done using the FF microplates (Biolog Inc., Hayward, CA, USA) containing 95 carbon sources and tetrazolium dye. The fungi P. lingam and P. biglobosus subclade 'brassicae' (3 isolates per group) were cultured on PDA medium for 6 weeks at 20 °C and then fungal spores were used as inoculum of microplates. The test was carried out in triplicate. We have demonstrated that substrate richness, calculated as the number of utilized substrates (measured at λ490 nm), and the number of substrates allowing effective growth of the isolates (λ750 nm), showed significant differences among tested species. The most efficient isolate of P. lingam utilized 36 carbon sources, whereas P. biglobosus utilized 60 substrates. Among them, 25-29 carbon sources for P. lingam and 34-48 substrates for P. biglobosus were efficiently used, allowing their growth. Cluster analysis based on Senath criteria divided P. biglobosus into two groups and P. lingam isolates formed one group (33% similarity). We deduce the similarities between the tested species help them coexist on the same host plant and the differences greatly contribute to their different lifestyles, with P. biglobosus being less specialized and P. lingam coevolving more strictly with the host plant.

Dissection of Quantitative Blackleg Resistance Reveals Novel Variants of Resistance Gene Rlm9 in Elite Brassica napus

Blackleg is one of the major fungal diseases in oilseed rape/canola worldwide. Most commercial cultivars carry R gene-mediated qualitative resistances that confer a high level of race-specific protection against Leptosphaeria maculans, the causal fungus of blackleg disease. However, monogenic resistances of this kind can potentially be rapidly overcome by mutations in the pathogen's avirulence genes. To counteract pathogen adaptation in this evolutionary arms race, there is a tremendous demand for quantitative background resistance to enhance durability and efficacy of blackleg resistance in oilseed rape. In this study, we characterized genomic regions contributing to quantitative L. maculans resistance by genome-wide association studies in a multiparental mapping population derived from six parental elite varieties exhibiting quantitative resistance, which were all crossed to one common susceptible parental elite variety. Resistance was screened using a fungal isolate with no corresponding avirulence (AvrLm) to major R genes present in the parents of the mapping population. Genome-wide association studies revealed eight significantly associated quantitative trait loci (QTL) on chromosomes A07 and A09, with small effects explaining 3-6% of the phenotypic variance. Unexpectedly, the qualitative blackleg resistance gene Rlm9 was found to be located within a resistance-associated haploblock on chromosome A07. Furthermore, long-range sequence data spanning this haploblock revealed high levels of single-nucleotide and structural variants within the Rlm9 coding sequence among the parents of the mapping population. The results suggest that novel variants of Rlm9 could play a previously unknown role in expression of quantitative disease resistance in oilseed rape.

Status and advances in mining for blackleg (Leptosphaeria maculans) quantitative resistance (QR) in oilseed rape (Brassica napus)

KEY MESSAGE

Quantitative resistance (QR) loci discovered through genetic and genomic analyses are abundant in the Brassica napus genome, providing an opportunity for their utilization in enhancing blackleg resistance. Quantitative resistance (QR) has long been utilized to manage blackleg in Brassica napus (canola, oilseed rape), even before major resistance genes (R-genes) were extensively explored in breeding programmes. In contrast to R-gene-mediated qualitative resistance, QR reduces blackleg symptoms rather than completely eliminating the disease. As a polygenic trait, QR is controlled by numerous genes with modest effects, which exerts less pressure on the pathogen to evolve; hence, its effectiveness is more durable compared to R-gene-mediated resistance. Furthermore, combining QR with major R-genes has been shown to enhance resistance against diseases in important crops, including oilseed rape. For these reasons, there has been a renewed interest among breeders in utilizing QR in crop improvement. However, the mechanisms governing QR are largely unknown, limiting its deployment. Advances in genomics are facilitating the dissection of the genetic and molecular underpinnings of QR, resulting in the discovery of several loci and genes that can be potentially deployed to enhance blackleg resistance. Here, we summarize the efforts undertaken to identify blackleg QR loci in oilseed rape using linkage and association analysis. We update the knowledge on the possible mechanisms governing QR and the advances in searching for the underlying genes. Lastly, we lay out strategies to accelerate the genetic improvement of blackleg QR in oilseed rape using improved phenotyping approaches and genomic prediction tools.

A Modified Meiotic Recombination in Brassica napus Largely Improves Its Breeding Efficiency

Simple Summary

The selection of varieties more resilient to disease and climate change requires generating new genetic diversity for breeding. The main mechanism for reshuffling genetic information is through the recombination of chromosomes during meiosis. We showed in oilseed rape (Brassica napus, AACC, 2n = 4x = 38), which is a natural hybrid formed from a cross between turnip (B. rapa, AA, 2n = 2x = 20) and cabbage (B. oleracea, CC, 2n = 2x = 18), that there is significantly more crossovers occurring along the entire A chromosomes in allotriploid AAC (crossbetween B. napus and B. rapa) than in diploid AA or allotetraploid AACC hybrids. We demonstrated that these allotriploid AAC hybrids are highly efficient to introduce new variability within oilseed rape varieties, notably by enabling the introduction of small genomic regions carrying genes controlling agronomically interesting traits.

Abstract

Meiotic recombination is the main tool used by breeders to generate biodiversity, allowing genetic reshuffling at each generation. It enables the accumulation of favorable alleles while purging deleterious mutations. However, this mechanism is highly regulated with the formation of one to rarely more than three crossovers, which are not randomly distributed. In this study, we showed that it is possible to modify these controls in oilseed rape (Brassica napus, AACC, 2n = 4x = 38) and that it is linked to AAC allotriploidy and not to polyploidy per se. To that purpose, we compared the frequency and the distribution of crossovers along A chromosomes from hybrids carrying exactly the same A nucleotide sequence, but presenting three different ploidy levels: AA, AAC and AACC. Genetic maps established with 202 SNPs anchored on reference genomes revealed that the crossover rate is 3.6-fold higher in the AAC allotriploid hybrids compared to AA and AACC hybrids. Using a higher SNP density, we demonstrated that smaller and numerous introgressions of B. rapa were present in AAC hybrids compared to AACC allotetraploid hybrids, with 7.6 Mb vs. 16.9 Mb on average and 21 B. rapa regions per plant vs. nine regions, respectively. Therefore, this boost of recombination is highly efficient to reduce the size of QTL carried in cold regions of the oilseed rape genome, as exemplified here for a QTL conferring blackleg resistance.

A gene-for-gene interaction involving a 'late' effector contributes to quantitative resistance to the stem canker disease in Brassica napus

The control of stem canker disease of Brassica napus (rapeseed), caused by the fungus Leptosphaeria maculans is based largely on plant genetic resistance: single-gene specific resistance (Rlm genes) or quantitative, polygenic, adult-stage resistance. Our working hypothesis was that quantitative resistance partly obeys the gene-for-gene model, with resistance genes 'recognizing' fungal effectors expressed during late systemic colonization. Five LmSTEE (stem-expressed effector) genes were selected and placed under the control of the AvrLm4-7 promoter, an effector gene highly expressed at the cotyledon stage of infection, for miniaturized cotyledon inoculation test screening of a gene pool of 204 rapeseed genotypes. We identified a rapeseed genotype, 'Yudal', expressing hypersensitive response to LmSTEE98. The LmSTEE98-RlmSTEE98 interaction was further validated by inactivation of the LmSTEE98 gene with a CRISPR-Cas9 approach. Isolates with mutated versions of LmSTEE98 induced more severe stem symptoms than the wild-type isolate in 'Yudal'. This single-gene resistance was mapped in a 0.6 cM interval of the 'Darmor_bzh' × 'Yudal' genetic map. One typical gene-for-gene interaction contributes partly to quantitative resistance when L. maculans colonizes the stems of rapeseed. With numerous other effectors specific to stem colonization, our study provides a new route for resistance gene discovery, elucidation of quantitative resistance mechanisms and selection for durable resistance.

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

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

Enhancing backcross programs through increased recombination

Background

Introgression of a quantitative trait locus (QTL) by successive backcrosses is used to improve elite lines (recurrent parent) by introducing alleles from exotic material (donor parent). In the absence of selection, the proportion of the donor genome decreases by half at each generation. However, since selection is for the donor allele at the QTL, elimination of the donor genome around that QTL will be much slower than in the rest of the genome (i.e. linkage drag). Using markers to monitor the genome around the QTL and in the genetic background can accelerate the return to the recurrent parent genome. Successful introgression of a locus depends partly on the occurrence of crossovers at favorable positions. However, the number of crossovers per generation is limited and their distribution along the genome is heterogeneous. Recently, techniques have been developed to modify these two recombination parameters.

Results

In this paper, we assess, by simulations in the context of Brassicaceae, the effect of increased recombination on the efficiency of introgression programs by studying the decrease in linkage drag and the recovery of the recurrent genome. The simulated selection schemes begin by two generations of foreground selection and continue with one or more generations of background selection. Our results show that, when the QTL is in a region that initially lacked crossovers, an increase in recombination rate can decrease linkage drag by nearly ten-fold after the foreground selection and improves the return to the recurrent parent. However, if the QTL is in a region that is already rich in crossovers, an increase in recombination rate is detrimental.

Conclusions

Depending on the recombination rate in the region targeted for introgression, increasing it can be beneficial or detrimental. Thus, the simulations analysed in this paper help us understand how an increase in recombination rate can be beneficial. They also highlight the best methods that can be used to increase recombination rate, depending on the situation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-021-00619-0.

Long-read assembly of the Brassica napus reference genome Darmor-bzh

BACKGROUND

The combination of long reads and long-range information to produce genome assemblies is now accepted as a common standard. This strategy not only allows access to the gene catalogue of a given species but also reveals the architecture and organization of chromosomes, including complex regions such as telomeres and centromeres. The Brassica genus is not exempt, and many assemblies based on long reads are now available. The reference genome for Brassica napus, Darmor-bzh, which was published in 2014, was produced using short reads and its contiguity was extremely low compared with current assemblies of the Brassica genus.

FINDINGS

Herein, we report the new long-read assembly of Darmor-bzh genome (Brassica napus) generated by combining long-read sequencing data and optical and genetic maps. Using the PromethION device and 6 flowcells, we generated ∼16 million long reads representing 93× coverage and, more importantly, 6× with reads longer than 100 kb. This ultralong-read dataset allows us to generate one of the most contiguous and complete assemblies of a Brassica genome to date (contig N50 > 10 Mb). In addition, we exploited all the advantages of the nanopore technology to detect modified bases and sequence transcriptomic data using direct RNA to annotate the genome and focus on resistance genes.

CONCLUSION

Using these cutting-edge technologies, and in particular by relying on all the advantages of the nanopore technology, we provide the most contiguous Brassica napus assembly, a resource that will be valuable to the Brassica community for crop improvement and will facilitate the rapid selection of agronomically important traits.

Meta-analysis of GWAS in canola blackleg (Leptosphaeria maculans) disease traits demonstrates increased power from imputed whole-genome sequence

Blackleg disease causes yield losses in canola (Brassica napus L.). To identify resistance genes and genomic regions, genome-wide association studies (GWAS) of 585 diverse winter and spring canola accessions were performed using imputed whole-genome sequence (WGS) and transcriptome genotype-by-sequencing (GBSt). Blackleg disease phenotypes were collected across three years in six trials. GWAS were performed in several ways and their respective power was judged by the number of significant single nucleotide polymorphisms (SNP), the false discovery rate (FDR), and the percentage of SNP that validated in additional field trials in two subsequent years. WGS GWAS with 1,234,708 million SNP detected a larger number of significant SNP, achieved a lower FDR and a higher validation rate than GBSt with 64,072 SNP. A meta-analysis combining survival and average internal infection resulted in lower FDR but also lower validation rates. The meta-analysis GWAS identified 79 genomic regions (674 SNP) conferring potential resistance to L. maculans. While several GWAS signals localised in regions of known Rlm genes, fifty-three new potential resistance regions were detected. Seventeen regions had underlying genes with putative functions related to disease defence or stress response in Arabidopsis thaliana. This study provides insight into the genetic architecture and potential molecular mechanisms underlying canola L. maculans resistance.

Genome-Wide Association Mapping Identifies Novel Loci for Quantitative Resistance to Blackleg Disease in Canola

Blackleg disease, caused by the fungal pathogen Leptosphaeria maculans, continues to be a major concern for sustainable production of canola (Brassica napus L.) in many parts of the world. The deployment of effective quantitative resistance (QR) is recognized as a durable strategy in providing natural defense to pathogens. Herein, we uncover loci for resistance to blackleg in a genetically diverse panel of canola accessions by exploiting historic recombination events which occurred during domestication and selective breeding by genome-wide association analysis (GWAS). We found extensive variation in resistance to blackleg at the adult plant stage, including for upper canopy infection. Using the linkage disequilibrium and genetic relationship estimates from 12,414 high quality SNPs, GWAS identified 59 statistically significant and "suggestive" SNPs on 17 chromosomes of B. napus genome that underlie variation in resistance to blackleg, evaluated under field and shade-house conditions. Each of the SNP association accounted for up to 25.1% of additive genetic variance in resistance among diverse panel of accessions. To understand the homology of QR genomic regions with Arabidopsis thaliana genome, we searched the synteny between QR regions with 22 ancestral blocks of Brassicaceae. Comparative analyses revealed that 25 SNP associations for QR were localized in nine ancestral blocks, as a result of genomic rearrangements. We further showed that phenological traits such as flowering time, plant height, and maturity confound the genetic variation in resistance. Altogether, these findings provided new insights on the complex genetic control of the blackleg resistance and further expanded our understanding of its genetic architecture.

Identification of resistance loci in Chinese and Canadian canola/rapeseed varieties against Leptosphaeria maculans based on genome-wide association studies

Background

The fungal pathogen Leptosphaeria maculans (Lm). causes blackleg disease on canola/rapeseed in many parts of the world. It is important to use resistant cultivars to manage the disease and minimize yield losses. In this study, twenty-two Lm isolates were used to identify resistance genes in a collection of 243 canola/rapeseed (Brassica napus L.) accessions from Canada and China. These Lm isolates carry different compliments of avirulence genes, and the investigation was based on a genome-wide association study (GWAS) and genotype-by-sequencing (GBS).

Results

Using the CROP-SNP pipeline, a total of 81,471 variants, including 78,632 SNPs and 2839 InDels, were identified. The GWAS was performed using TASSEL 5.0 with GLM + Q model. Thirty-two and 13 SNPs were identified from the Canadian and Chinese accessions, respectively, tightly associated with blackleg resistance with P values < 1 × 10^− 4. These SNP loci were distributed on chromosomes A03, A05, A08, A09, C01, C04, C05, and C07, with the majority of them on A08 followed by A09 and A03. The significant SNPs identified on A08 were all located in a 2010-kb region and associated with resistance to 12 of the 22 Lm isolates. Furthermore, 25 resistance gene analogues (RGAs) were identified in these regions, including two nucleotide binding site (NBS) domain proteins, fourteen RLKs, three RLPs and six TM-CCs. These RGAs can be the potential candidate genes for blackleg resistance.

Conclusion

This study provides insights into potentially new genomic regions for discovery of additional blackleg resistance genes. The identified regions associated with blackleg resistance in the germplasm collection may also contribute directly to the development of canola varieties with novel resistance genes against blackleg of canola.

Exploring Mechanisms of Quantitative Resistance to Leptosphaeria maculans (Blackleg) in the Cotyledons of Canola (Brassica napus) Based on Transcriptomic and Microscopic Analyses

Using resistant cultivars is a common approach to managing blackleg of canola/rapeseed caused by Leptosphaeria maculans (Lm). Quantitative resistance (QR), as opposed to major-gene resistance, is of interest because it is generally more durable, due to its multi-genetic basis. However, the mechanisms and genes underlying QR are mostly unknown. In this study, potential QR modes of action in “74-44 BL” was explored. This Canadian canola cultivar showed moderate but consistent race-nonspecific resistance at the cotyledon and adult-plant stages. A susceptible cultivar, “Westar”, was used as a control. After inoculation, the lesions developed more slowly on the cotyledons of 74-44 BL than those of Westar. We used RNA sequencing (-RNA-seq) to identify genes and their functions, putatively related to this resistance, and found that genes involved in programmed cell death (PCD), reactive oxygen species (ROS), signal transduction or intracellular endomembrane transport were most differentially expressed. ROS production was assessed in relation to Lm hyphal growth and lesion size; it occurred beyond the tissue colonized by Lm in 74-44 BL and appeared to trigger rapid cell death, limiting cotyledon colonization by Lm. In contrast, Lm grew more rapidly in Westar, often catching up with the ring of ROS and surpassing lesion boundaries. It appears that QR in 74-44 BL cotyledons is associated with limited colonization by Lm possibly mediated via ROS. The RNA-seq data also showed a link between ROS, signal transduction, and endomembrane vesicle trafficking, as well as PCD in the resistance. These results provide a starting point for a better understanding of the mechanisms behind QR against Lm in canola.

Quantitative phenotyping of shell suture strength in walnut (Juglans regia L.) enhances precision for detection of QTL and genome-wide association mapping

Walnut shell suture strength directly impacts the ability to maintain shell integrity during harvest and processing, susceptibility to insect damage and other contamination, and the proportion of kernel halves recovered during cracking. Suture strength is therefore an important breeding objective. Here, two methods of phenotyping this trait were investigated: 1) traditional, qualitative and rather subjective scoring on an interval scale by human observers, and; 2) quantitative and continuous measurements captured by a texturometer. The aim of this work was to increase the accuracy of suture strength phenotyping and to then apply two mapping approaches, quantitative trait loci (QTL) mapping and genome wide association (GWAS) models, in order to dissect the genetic basis of the walnut suture trait. Using data collected on trees within the UC Davis Walnut Improvement Program (n = 464), the genetic correlation between the texturometer method and qualitatively scored method was high (0.826). Narrow sense heritability calculated using quantitative measurements was 0.82. A major QTL for suture strength was detected on LG05, explaining 34% of the phenotypic variation; additionally, two minor QTLs were identified on LG01 and LG11. All three QTLs were confirmed with GWAS on corresponding chromosomes. The findings reported in this study are relevant for application towards a molecular breeding program in walnut.

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

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

SNP-based association study of kernel architecture in a worldwide collection of durum wheat germplasm

Durum wheat, genetic resource with favorable alleles is considered as natural gene pool for wheat breeding. Kernel size and weight are important factors affecting grain yield in crops. Here, association analysis was performed to dissect the genetic constitution of kernel-related traits in 150 lines collected from 46 countries and regions using a set of EST-derived and genome-wide SNP markers with five consecutive years of data. Total 109 significant associations for eight kernel-related traits were detected under a mix linear model, generating 54 unique SNP markers distributed on 13 of 14 chromosomes. Of which, 19 marker-trait associations were identified in two or more environments, including one stable and pleiotropic SNP BE500291_5_A_37 on chromosome 5A correlated with six kernel traits. Although most of our SNP loci were overlapped with the previously known kernel weight QTLs, several novel loci for kernel traits in durum were reported. Correlation analysis implied that the moderate climatic variables during growth and development of durum are needed for the large grain size and high grain weight. Combined with our previous studies, we found that chromosome 5A might play an important role in durum growth and development.

Oilseed rape (Brassica napus) resistance to growth of Leptosphaeria maculans in leaves of young plants contributes to quantitative resistance in stems of adult plants

Key message: One QTL for resistance against Leptosphaeria maculans growth in leaves of young plants in controlled environments overlapped with one QTL detected in adult plants in field experiments. The fungal pathogen Leptosphaeria maculans initially infects leaves of oilseed rape (Brassica napus) in autumn in Europe and then grows systemically from leaf lesions along the leaf petiole to the stem, where it causes damaging phoma stem canker (blackleg) in summer before harvest. Due to the difficulties of investigating resistance to L. maculans growth in leaves and petioles under field conditions, identification of quantitative resistance typically relies on end of season stem canker assessment on adult plants. To investigate whether quantitative resistance can be detected in young plants, we first selected nine representative DH (doubled haploid) lines from an oilseed rape DY ('Darmor-bzh' × 'Yudal') mapping population segregating for quantitative resistance against L. maculans for controlled environment experiment (CE). We observed a significant correlation between distance grown by L. maculans along the leaf petiole towards the stem (r = 0.91) in CE experiments and the severity of phoma stem canker in field experiments. To further investigate quantitative trait loci (QTL) related to resistance against growth of L. maculans in leaves of young plants in CE experiments, we selected 190 DH lines and compared the QTL detected in CE experiments with QTL related to stem canker severity in stems of adult plants in field experiments. Five QTL for resistance to L. maculans growth along the leaf petiole were detected; collectively they explained 35% of the variance. Two of these were also detected in leaf lesion area assessments and each explained 10-12% of the variance. One QTL on A02 co-localized with a QTL detected in stems of adult plants in field experiments. This suggests that resistance to the growth of L. maculans from leaves along the petioles towards the stems contributes to the quantitative resistance assessed in stems of adult plants in field experiments at the end of the growing season.

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.

Stable Quantitative Resistance Loci to Blackleg Disease in Canola (Brassica napus L.) Over Continents

The hemibiotrophic fungus, Leptosphaeria maculans is the most devastating pathogen, causing blackleg disease in canola (Brassica napus L). To study the genomic regions involved in quantitative resistance (QR), 259-276 DH lines from Darmor-bzh/Yudal (DYDH) population were assessed for resistance to blackleg under shade house and field conditions across 3 years. In different experiments, the broad sense heritability varied from 43 to 95%. A total of 27 significant quantitative trait loci (QTL) for QR were detected on 12 chromosomes and explained between 2.14 and 10.13% of the genotypic variance. Of the significant QTL, at least seven were repeatedly detected across different experiments on chromosomes A02, A07, A09, A10, C01, and C09. Resistance alleles were mainly contributed by 'Darmor-bzh' but 'Yudal' also contributed few of them. Our results suggest that plant maturity and plant height may have a pleiotropic effect on QR in our conditions. We confirmed that Rlm9 which is present in 'Darmor-bzh' is not effective to confer resistance in our Australian field conditions. Comparative mapping showed that several R genes coding for nucleotide-binding leucine-rich repeat (LRR) receptors map in close proximity (within 200 Kb) of the significant trait-marker associations on the reference 'Darmor-bzh' genome assembly. More importantly, eight significant QTL regions were detected across diverse growing environments: Australia, France, and United Kingdom. These stable QTL identified herein can be utilized for enhancing QR in elite canola germplasm via marker- assisted or genomic selection strategies.