Homoeologous duplicated regions are involved in quantitative resistance of Brassica napus to stem canker

Genetic control of root morphology in response to nitrogen across rapeseed diversity

Rapeseed (Brassica napus L.) is an oil-containing crop of great economic value but with considerable nitrogen requirement. Breeding root systems that efficiently absorb nitrogen from the soil could be a driver to ensure genetic gains for more sustainable rapeseed production. The aim of this study is to identify genomic regions that regulate root morphology in response to nitrate availability. The natural variability offered by 300 inbred lines was screened at two experimental locations. Seedlings grew hydroponically with low or elevated nitrate levels. Fifteen traits related to biomass production and root morphology were measured. On average across the panel, a low nitrate level increased the root-to-shoot biomass ratio and the lateral root length. A large phenotypic variation was observed, along with important heritability values and genotypic effects, but low genotype-by-nitrogen interactions. Genome-wide association study and bulk segregant analysis were used to identify loci regulating phenotypic traits. The first approach nominated 319 SNPs that were combined into 80 QTLs. Three QTLs identified on the A07 and C07 chromosomes were stable across nitrate levels and/or experimental locations. The second approach involved genotyping two groups of individuals from an experimental F2 population created by crossing two accessions with contrasting lateral root lengths. These individuals were found in the tails of the phenotypic distribution. Co-localized QTLs found in both mapping approaches covered a chromosomal region on the A06 chromosome. The QTL regions contained some genes putatively involved in root organogenesis and represent selection targets for redesigning the root morphology of rapeseed.

Key Advances in the New Era of Genomics-Assisted Disease Resistance Improvement of Brassica Species

Disease resistance improvement remains a major focus in breeding programs as diseases continue to devastate Brassica production systems due to intensive cultivation and climate change. Genomics has paved the way to understand the complex genomes of Brassicas, which has been pivotal in the dissection of the genetic underpinnings of agronomic traits driving the development of superior cultivars. The new era of genomics-assisted disease resistance breeding has been marked by the development of high-quality genome references, accelerating the identification of disease resistance genes controlling both qualitative (major) gene and quantitative resistance. This facilitates the development of molecular markers for marker assisted selection and enables genome editing approaches for targeted gene manipulation to enhance the genetic value of disease resistance traits. This review summarizes the key advances in the development of genomic resources for Brassica species, focusing on improved genome references, based on long-read sequencing technologies and pangenome assemblies. This is further supported by the advances in pathogen genomics, which have resulted in the discovery of pathogenicity factors, complementing the mining of disease resistance genes in the host. Recognizing the co-evolutionary arms race between the host and pathogen, it is critical to identify novel resistance genes using crop wild relatives and synthetic cultivars or through genetic manipulation via genome-editing to sustain the development of superior cultivars. Integrating these key advances with new breeding techniques and improved phenotyping using advanced data analysis platforms will make disease resistance improvement in Brassica species more efficient and responsive to current and future demands.

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.

Functional divergence of Brassica napus BnaABI1 paralogs in the structurally conserved PP2CA gene subfamily of Brassicaceae

Group A PP2C (PP2CA) genes form a gene subfamily whose members play an important role in regulating many biological processes by dephosphorylation of target proteins. In this study we examined the effects of evolutionary changes responsible for functional divergence of BnaABI1 paralogs in Brassica napus against the background of the conserved PP2CA gene subfamily in Brassicaceae. We performed comprehensive phylogenetic analyses of 192 PP2CA genes in 15 species in combination with protein structure homology modeling. Fundamentally, the number of PP2CA genes remained relatively constant in these taxa, except in the Brassica genus and Camelina sativa. The expansion of this gene subfamily in these species has resulted from whole genome duplication. We demonstrated a high degree of structural conservation of the PP2CA genes, with a few minor variations between the different PP2CA groups. Furthermore, the pattern of conserved sequence motifs in the PP2CA proteins and their secondary and 3D structures revealed strong conservation of the key ion-binding sites. Syntenic analysis of triplicated regions including ABI1 paralogs revealed significant structural rearrangements of the Brassica genomes. The functional and syntenic data clearly show that triplication of BnaABI1 in B. napus has had an impact on its functions, as well as the positions of adjacent genes in the corresponding chromosomal regions. The expression profiling of BnaABI1 genes showed functional divergence, i.e. subfunctionalization, potentially leading to neofunctionalization. These differences in expression are likely due to changes in the promoters of the BnaABI1 paralogs. Our results highlight the complexity of PP2CA gene subfamily evolution in Brassicaceae.

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.

Whole-Genome Sequence of Synthesized Allopolyploids in Cucumis Reveals Insights into the Genome Evolution of Allopolyploidization

The importance of allopolyploidy in plant evolution has been widely recognized. The genetic changes triggered by allopolyploidy, however, are not yet fully understood due to inconsistent phenomena reported across diverse species. The construction of synthetic polyploids offers a controlled approach to systematically reveal genomic changes that occur during the process of polyploidy. This study reports the first fully sequenced synthetic allopolyploid constructed from a cross between Cucumis sativus and C. hystrix, with high-quality assembly. The two subgenomes are confidently partitioned and the C. sativus-originated subgenome predominates over the C. hystrix-originated subgenome, retaining more sequences and showing higher homeologous gene expression. Most of the genomic changes emerge immediately after interspecific hybridization. Analysis of a series of genome sequences from several generations (S0, S4-S13) of C. ×hytivus confirms that genomic changes occurred in the very first generations, subsequently slowing down as the process of diploidization is initiated. The duplicated genome of the allopolyploid with double genes from both parents broadens the genetic base of C. ×hytivus, resulting in enhanced phenotypic plasticity. This study provides novel insights into plant polyploid genome evolution and demonstrates a promising strategy for the development of a wide array of novel plant species and varieties through artificial polyploidization.

Identification of resistance loci against new pathotypes of Plasmodiophora brassicae in Brassica napus based on genome-wide association mapping

Genetic resistance is a successful strategy for management of clubroot (Plasmodiophora brassicae) of brassica crops, but resistance can break down quickly. Identification of novel sources of resistance is especially important when new pathotypes arise. In the current study, the reaction of 177 accessions of Brassica napus to four new, virulent pathotypes of P. brassicae was assessed. Each accession was genotyped using genotyping by sequencing to identify and map novel sources of clubroot resistance using mixed linear model (MLM) analysis. The majority of accessions were highly susceptible (70–100 DSI), but a few accessions exhibited strong resistance (0–20 DSI) to pathotypes 5X (21 accessions), 3A (8), 2B (7), and 3D (15), based on the Canadian Clubroot Differential system. In total, 301,753 SNPs were mapped to 19 chromosomes. Population structure analysis indicated that the 177 accessions belong to seven major populations. SNPs were associated with resistance to each pathotype using MLM. In total, 13 important SNP loci were identified, with 9 SNPs mapped to the A-genome and 4 to the C-genome. The SNPs were associated with resistance to pathotypes 5X (2 SNPs), 3A (4), 2B (5) and 3D (6). A Blast search of 1.6 Mb upstream and downstream from each SNP identified 13 disease-resistance genes or domains. The distance between a SNP locus and the nearest resistance gene ranged from 0.04 to 0.74 Mb. The resistant lines and SNP markers identified in this study can be used to breed for resistance to the most prevalent new pathotypes of P. brassicae in Canada.

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.

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.

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.

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.

Detection of First Marker Trait Associations for Resistance Against Sclerotinia sclerotiorum in Brassica juncea-Erucastrum cardaminoides Introgression Lines

A set of 96 Brassica juncea-Erucastrum cardaminoides introgression lines (ILs) were developed with genomic regions associated with Sclerotinia stem rot (Sclerotinia sclerotiorum) resistance from a wild Brassicaceous species E. cardaminoides. ILs were assessed for their resistance responses to stem inoculation with S. sclerotiorum, over three crop seasons (season I, 2011/2012; II, 2014/2015; III, 2016-2017). Initially, ILs were genotyped with transferable SSR markers and subsequently through genotyping by sequencing. SSR based association mapping identified six marker loci associated to resistance in both A and B genomes. Subsequent genome-wide association analysis (GWAS) of 84 ILs recognized a large number of SNPs associated to resistance, in chromosomes A03, A06, and B03. Chromosomes A03 and A06 harbored the maximum number of resistance related SNPs. Annotation of linked genomic regions highlighted an array of resistance mechanisms in terms of signal transduction pathways, hypersensitive responses and production of anti-fungal proteins and metabolites. Of major importance was the clustering of SNPs, encoding multiple resistance genes on small regions spanning approximately 885 kb region on chromosome A03 and 74 kb on B03. Five SNPs on chromosome A03 (6,390,210-381) were associated with LRR-RLK (receptor like kinases) genes that encode LRR-protein kinase family proteins. Genetic factors associated with pathogen-associated molecular patterns (PAMPs) and effector-triggered immunity (ETI) were predicted on chromosome A03, exhibiting 11 SNPs (6,274,763-994). These belonged to three R-Genes encoding TIR-NBS-LRR proteins. Marker trait associations (MTAs) identified will facilitate marker assisted introgression of these critical resistances, into new cultivars of B. juncea initially and, subsequently, into other crop Brassica species.

Finding invisible quantitative trait loci with missing data

Evolutionary processes during plant polyploidization and speciation have led to extensive presence-absence variation (PAV) in crop genomes, and there is increasing evidence that PAV associates with important traits. Today, high-resolution genetic analysis in major crops frequently implements simple, cost-effective, high-throughput genotyping from single nucleotide polymorphism (SNP) hybridization arrays; however, these are normally not designed to distinguish PAV from failed SNP calls caused by hybridization artefacts. Here, we describe a strategy to recover valuable information from single nucleotide absence polymorphisms (SNaPs) by population-based quality filtering of SNP hybridization data to distinguish patterns associated with genuine deletions from those caused by technical failures. We reveal that including SNaPs in genetic analyses elucidate segregation of small to large-scale structural variants in nested association mapping populations of oilseed rape (Brassica napus), a recent polyploid crop with widespread structural variation. Including SNaP markers in genomewide association studies identified numerous quantitative trait loci, invisible using SNP markers alone, for resistance to two major fungal diseases of oilseed rape, Sclerotinia stem rot and blackleg disease. Our results indicate that PAV has a strong influence on quantitative disease resistance in B. napus and that SNaP analysis using cost-effective SNP array data can provide extensive added value from 'missing data'. This strategy might also be applicable for improving the precision of genetic mapping in many important crop species.

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.

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

A repertoire of the genomic regions involved in quantitative resistance to Leptosphaeria maculans in winter oilseed rape was established from combined linkage-based QTL and genome-wide association (GWA) mapping. Linkage-based mapping of quantitative trait loci (QTL) and genome-wide association studies are complementary approaches for deciphering the genomic architecture of complex agronomical traits. In oilseed rape, quantitative resistance to blackleg disease, caused by L. maculans, is highly polygenic and is greatly influenced by the environment. In this study, we took advantage of multi-year data available on three segregating populations derived from the resistant cv Darmor and multi-year data available on oilseed rape panels to obtain a wide overview of the genomic regions involved in quantitative resistance to this pathogen in oilseed rape. Sixteen QTL regions were common to at least two biparental populations, of which nine were the same as previously detected regions in a multi-parental design derived from different resistant parents. Eight regions were significantly associated with quantitative resistance, of which five on A06, A08, A09, C01 and C04 were located within QTL support intervals. Homoeologous Brassica napus genes were found in eight homoeologous QTL regions, which corresponded to 657 pairs of homoeologous genes. Potential candidate genes underlying this quantitative resistance were identified. Genomic predictions and breeding are also discussed, taking into account the highly polygenic nature of this resistance.

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.

Mapping of homoeologous chromosome exchanges influencing quantitative trait variation in Brassica napus

Genomic rearrangements arising during polyploidization are an important source of genetic and phenotypic variation in the recent allopolyploid crop Brassica napus. Exchanges among homoeologous chromosomes, due to interhomoeologue pairing, and deletions without compensating homoeologous duplications are observed in both natural B. napus and synthetic B. napus. Rearrangements of large or small chromosome segments induce gene copy number variation (CNV) and can potentially cause phenotypic changes. Unfortunately, complex genome restructuring is difficult to deal with in linkage mapping studies. Here, we demonstrate how high-density genetic mapping with codominant, physically anchored SNP markers can detect segmental homoeologous exchanges (HE) as well as deletions and accurately link these to QTL. We validated rearrangements detected in genetic mapping data by whole-genome resequencing of parental lines along with cytogenetic analysis using fluorescence in situ hybridization with bacterial artificial chromosome probes (BAC-FISH) coupled with PCR using primers specific to the rearranged region. Using a well-known QTL region influencing seed quality traits as an example, we confirmed that HE underlies the trait variation in a DH population involving a synthetic B. napus trait donor, and succeeded in narrowing the QTL to a small defined interval that enables delineation of key candidate genes.

The Impact of Open Pollination on the Structural Evolutionary Dynamics, Meiotic Behavior, and Fertility of Resynthesized Allotetraploid Brassica napus L

Allopolyploidy, which results from the merger and duplication of two divergent genomes, has played a major role in the evolution and diversification of flowering plants. The genomic changes that occur in resynthesized or natural neopolyploids have been extensively studied, but little is known about the effects of the reproductive mode in the initial generations that may precede its successful establishment. To truly reflect the early generations of a nascent polyploid, two resynthesized allotetraploid Brassica napus populations were obtained for the first time by open pollination. In these populations, we detected a much lower level of aneuploidy (third generation) compared with those previously published populations obtained by controlled successive selfing. We specifically studied 33 resynthesized B. napus individuals from our two open pollinated populations, and showed that meiosis was affected in both populations. Their genomes were deeply shuffled after allopolyploidization: up to 8.5 and 3.5% of the C and A subgenomes were deleted in only two generations. The identified deletions occurred mainly at the distal part of the chromosome, and to a significantly greater extent on the C rather than the A subgenome. Using Fluorescent In Situ Hybridization (BAC-FISH), we demonstrated that four of these deletions corresponded to fixed translocations (via homeologous exchanges). We were able to evaluate the size of the structural variations and their impact on the whole genome size, gene content, and allelic diversity. In addition, the evolution of fertility was assessed, to better understand the difficulty encountered by novel polyploid individuals before the putative formation of a novel stable species.

Genome-Wide Association Studies In Plant Pathosystems: Toward an Ecological Genomics Approach

The emergence and re-emergence of plant pathogenic microorganisms are processes that imply perturbations in both host and pathogen ecological niches. Global change is largely assumed to drive the emergence of new etiological agents by altering the equilibrium of the ecological habitats which in turn places hosts more in contact with pathogen reservoirs. In this context, the number of epidemics is expected to increase dramatically in the next coming decades both in wild and crop plants. Under these considerations, the identification of the genetic variants underlying natural variation of resistance is a pre-requisite to estimate the adaptive potential of wild plant populations and to develop new breeding resistant cultivars. On the other hand, the prediction of pathogen's genetic determinants underlying disease emergence can help to identify plant resistance alleles. In the genomic era, whole genome sequencing combined with the development of statistical methods led to the emergence of Genome Wide Association (GWA) mapping, a powerful tool for detecting genomic regions associated with natural variation of disease resistance in both wild and cultivated plants. However, GWA mapping has been less employed for the detection of genetic variants associated with pathogenicity in microbes. Here, we reviewed GWA studies performed either in plants or in pathogenic microorganisms (bacteria, fungi and oomycetes). In addition, we highlighted the benefits and caveats of the emerging joint GWA mapping approach that allows for the simultaneous identification of genes interacting between genomes of both partners. Finally, based on co-evolutionary processes in wild populations, we highlighted a phenotyping-free joint GWA mapping approach as a promising tool for describing the molecular landscape underlying plant - microbe interactions.

Genome-wide Association Study Identifies New Loci for Resistance to Leptosphaeria maculans in Canola

Key message "We identified both quantitative and quantitative resistance loci to Leptosphaeria maculans, a fungal pathogen, causing blackleg disease in canola. Several genome-wide significant associations were detected at known and new loci for blackleg resistance. We further validated statistically significant associations in four genetic mapping populations, demonstrating that GWAS marker loci are indeed associated with resistance to L. maculans. One of the novel loci identified for the first time, Rlm12, conveys adult plant resistance in canola." Blackleg, caused by Leptosphaeria maculans, is a significant disease which affects the sustainable production of canola (Brassica napus). This study reports a genome-wide association study based on 18,804 polymorphic SNPs to identify loci associated with qualitative and quantitative resistance to L. maculans. Genomic regions delimited with 694 significant SNP markers, that are associated with resistance evaluated using 12 single spore isolates and pathotypes from four canola stubble were identified. Several significant associations were detected at known disease resistance loci including in the vicinity of recently cloned Rlm2/LepR3 genes, and at new loci on chromosomes A01/C01, A02/C02, A03/C03, A05/C05, A06, A08, and A09. In addition, we validated statistically significant associations on A01, A07, and A10 in four genetic mapping populations, demonstrating that GWAS marker loci are indeed associated with resistance to L. maculans. One of the novel loci identified for the first time, Rlm12, conveys adult plant resistance and mapped within 13.2 kb from Arabidopsis R gene of TIR-NBS class. We showed that resistance loci are located in the vicinity of R genes of Arabidopsis thaliana and Brassica napus on the sequenced genome of B. napus cv. Darmor-bzh. Significantly associated SNP markers provide a valuable tool to enrich germplasm for favorable alleles in order to improve the level of resistance to L. maculans in canola.

Genetic basis of nitrogen use efficiency and yield stability across environments in winter rapeseed

Background

Nitrogen use efficiency is an important breeding trait that can be modified to improve the sustainability of many crop species used in agriculture. Rapeseed is a major oil crop with low nitrogen use efficiency, making its production highly dependent on nitrogen input. This complex trait is suspected to be sensitive to genotype × environment interactions, especially genotype × nitrogen interactions. Therefore, phenotyping diverse rapeseed populations under a dense network of trials is a powerful approach to study nitrogen use efficiency in this crop. The present study aimed to determine the quantitative trait loci (QTL) associated with yield in winter oilseed rape and to assess the stability of these regions under contrasting nitrogen conditions for the purpose of increasing nitrogen use efficiency.

Results

Genome-wide association studies and linkage analyses were performed on two diversity sets and two doubled-haploid populations. These populations were densely genotyped, and yield-related traits were scored in a multi-environment design including seven French locations, six growing seasons (2009 to 2014) and two nitrogen nutrition levels (optimal versus limited). Very few genotype × nitrogen interactions were detected, and a large proportion of the QTL were stable across nitrogen nutrition conditions. In contrast, strong genotype × trial interactions in which most of the QTL were specific to a single trial were found. To obtain further insight into the QTL × environment interactions, genetic analyses of ecovalence were performed to identify the genomic regions contributing to the genotype × nitrogen and genotype × trial interactions. Fifty-one critical genomic regions contributing to the additive genetic control of yield-associated traits were identified, and the structural organization of these regions in the genome was investigated.

Conclusions

Our results demonstrated that the effect of the trial was greater than the effect of nitrogen nutrition levels on seed yield-related traits under our experimental conditions. Nevertheless, critical genomic regions associated with yield that were stable across environments were identified in rapeseed.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-016-0432-z) contains supplementary material, which is available to authorized users.

Genome-wide association analysis and differential expression analysis of resistance to Sclerotinia stem rot in Brassica napus

Brassica napus is one of the most important oil crops in the world, and stem rot caused by the fungus Sclerotinia sclerotiorum results in major losses in yield and quality. To elucidate resistance genes and pathogenesis-related genes, genome-wide association analysis of 347 accessions was performed using the Illumina 60K Brassica SNP (single nucleotide polymorphism) array. In addition, the detached stem inoculation assay was used to select five highly resistant (R) and susceptible (S) B. napus lines, 48 h postinoculation with S. sclerotiorum for transcriptome sequencing. We identified 17 significant associations for stem resistance on chromosomes A8 and C6, five of which were on A8 and 12 on C6. The SNPs identified on A8 were located in a 409-kb haplotype block, and those on C6 were consistent with previous QTL mapping efforts. Transcriptome analysis suggested that S. sclerotiorum infection activates the immune system, sulphur metabolism, especially glutathione (GSH) and glucosinolates in both R and S genotypes. Genes found to be specific to the R genotype related to the jasmonic acid pathway, lignin biosynthesis, defence response, signal transduction and encoding transcription factors. Twenty-four genes were identified in both the SNP-trait association and transcriptome sequencing analyses, including a tau class glutathione S-transferase (GSTU) gene cluster. This study provides useful insight into the molecular mechanisms underlying the plant's response to S. sclerotiorum.

Identification of environmentally stable QTL for resistance against Leptosphaeria maculans in oilseed rape (Brassica napus)

Six stable QTL for resistance against L. maculans (phoma stem canker) have been identified by QTL × environment interaction analysis using data from five winter oilseed rape field experiments. Phoma stem canker, caused by Leptosphaeria maculans, is a disease of worldwide importance on oilseed rape (Brassica napus). Quantitative trait loci (QTL)-mediated resistance against L. maculans in B. napus is considered to be race non-specific and potentially durable. Identification and evaluation of QTL for resistance to L. maculans is important for breeding oilseed rape cultivars with durable resistance. An oilseed rape mapping population was used to detect QTL for resistance against L. maculans in five winter oilseed rape field experiments under different environments. A total of 17 QTL involved in 'field' quantitative resistance against L. maculans were detected and collectively explained 51% of the phenotypic variation. The number of QTL detected in each experiment ranged from two to nine and individual QTL explained 2-25% of the phenotypic variation. QTL × environment interaction analysis suggested that six of these QTL were less sensitive to environmental factors, so they were considered to be stable QTL. Markers linked to these stable QTL will be valuable for selection to breed for effective resistance against L. maculans in different environments, which will contribute to sustainable management of the disease.

Genome-Wide Analysis of Seed Acid Detergent Lignin (ADL) and Hull Content in Rapeseed (Brassica napus L.)

A stable yellow-seeded variety is the breeding goal for obtaining the ideal rapeseed (Brassica napus L.) plant, and the amount of acid detergent lignin (ADL) in the seeds and the hull content (HC) are often used as yellow-seeded rapeseed screening indices. In this study, a genome-wide association analysis of 520 accessions was performed using the Q + K model with a total of 31,839 single-nucleotide polymorphism (SNP) sites. As a result, three significant associations on the B. napus chromosomes A05, A09, and C05 were detected for seed ADL content. The peak SNPs were within 9.27, 14.22, and 20.86 kb of the key genes BnaA.PAL4, BnaA.CAD2/BnaA.CAD3, and BnaC.CCR1, respectively. Further analyses were performed on the major locus of A05, which was also detected in the seed HC examination. A comparison of our genome-wide association study (GWAS) results and previous linkage mappings revealed a common chromosomal region on A09, which indicates that GWAS can be used as a powerful complementary strategy for dissecting complex traits in B. napus. Genomic selection (GS) utilizing the significant SNP markers based on the GWAS results exhibited increased predictive ability, indicating that the predictive ability of a given model can be substantially improved by using GWAS and GS.

Comparative genomic analysis of duplicated homoeologous regions involved in the resistance of Brassica napus to stem canker

All crop species are current or ancient polyploids. Following whole genome duplication, structural and functional modifications result in differential gene content or regulation in the duplicated regions, which can play a fundamental role in the diversification of genes underlying complex traits. We have investigated this issue in Brassica napus, a species with a highly duplicated genome, with the aim of studying the structural and functional organization of duplicated regions involved in quantitative resistance to stem canker, a disease caused by the fungal pathogen Leptosphaeria maculans. Genome-wide association analysis on two oilseed rape panels confirmed that duplicated regions of ancestral blocks E, J, R, U, and W were involved in resistance to stem canker. The structural analysis of the duplicated genomic regions showed a higher gene density on the A genome than on the C genome and a better collinearity between homoeologous regions than paralogous regions, as overall in the whole B. napus genome. The three ancestral sub-genomes were involved in the resistance to stem canker and the fractionation profile of the duplicated regions corresponded to what was expected from results on the B. napus progenitors. About 60% of the genes identified in these duplicated regions were single-copy genes while less than 5% were retained in all the duplicated copies of a given ancestral block. Genes retained in several copies were mainly involved in response to stress, signaling, or transcription regulation. Genes with resistance-associated markers were mainly retained in more than two copies. These results suggested that some genes underlying quantitative resistance to stem canker might be duplicated genes. Genes with a hydrolase activity that were retained in one copy or R-like genes might also account for resistance in some regions. Further analyses need to be conducted to indicate to what extent duplicated genes contribute to the expression of the resistance phenotype.

Karyotype and gene order evolution from reconstructed extinct ancestors highlight contrasts in genome plasticity of modern rosid crops

We used nine complete genome sequences, from grape, poplar, Arabidopsis, soybean, lotus, apple, strawberry, cacao, and papaya, to investigate the paleohistory of rosid crops. We characterized an ancestral rosid karyotype, structured into 7/21 protochomosomes, with a minimal set of 6,250 ordered protogenes and a minimum physical coding gene space of 50 megabases. We also proposed ancestral karyotypes for the Caricaceae, Brassicaceae, Malvaceae, Fabaceae, Rosaceae, Salicaceae, and Vitaceae families with 9, 8, 10, 6, 12, 9, 12, and 19 protochromosomes, respectively. On the basis of these ancestral karyotypes and present-day species comparisons, we proposed a two-step evolutionary scenario based on allohexaploidization involving the newly characterized A, B, and C diploid progenitors leading to dominant (stable) and sensitive (plastic) genomic compartments in any modern rosid crops. Finally, a new user-friendly online tool, “DicotSyntenyViewer” (available from http://urgi.versailles.inra.fr/synteny-dicot), has been made available for accurate translational genomics in rosids.