Molecular mapping of qualitative and quantitative loci for resistance to Leptosphaeria maculans causing blackleg disease in canola (Brassica napus L.)

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.

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

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

The Global Assessment of Oilseed Brassica Crop Species Yield, Yield Stability and the Underlying Genetics

The global demand for oilseeds is increasing along with the human population. The family of Brassicaceae crops are no exception, typically harvested as a valuable source of oil, rich in beneficial molecules important for human health. The global capacity for improving Brassica yield has steadily risen over the last 50 years, with the major crop Brassica napus (rapeseed, canola) production increasing to ~72 Gt in 2020. In contrast, the production of Brassica mustard crops has fluctuated, rarely improving in farming efficiency. The drastic increase in global yield of B. napus is largely due to the demand for a stable source of cooking oil. Furthermore, with the adoption of highly efficient farming techniques, yield enhancement programs, breeding programs, the integration of high-throughput phenotyping technology and establishing the underlying genetics, B. napus yields have increased by >450 fold since 1978. Yield stability has been improved with new management strategies targeting diseases and pests, as well as by understanding the complex interaction of environment, phenotype and genotype. This review assesses the global yield and yield stability of agriculturally important oilseed Brassica species and discusses how contemporary farming and genetic techniques have driven improvements.

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.

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.

Characterization of non-specific lipid transfer protein (nsLtp) gene families in the Brassica napus pangenome reveals abundance variation

BACKGROUND

Brassica napus is an important agricultural species, improving stress resistance was one of the main breeding goals at present. Non-specific lipid transfer proteins (nsLTPs) are small, basic proteins which are involved in some biotic or abiotic stress responses. B. napus is susceptible to a variety of fungal diseases, so identify the BnLTPs and their expression in disease responses is very important. The common reference genome of B. napus does not contain all B. napus genes because of gene presence/absence variations between individuals. Therefore, it was necessary to search for candidate BnLTP genes in the B. napus pangenome.

RESULTS

In the present study, the BnLTP genes were identified throughout the pangenome, and different BnLTP genes were presented among varieties. Totally, 246 BnLTP genes were identified and could be divided into five types (1, 2, C, D, and G). The classification, phylogenetic reconstruction, chromosome distribution, functional annotation, and gene expression were analyzed. We also identified potential cis-elements that respond to biotic and abiotic stresses in the 2 kb upstream regions of all BnLTP genes. RNA sequencing analysis showed that the BnLTP genes were involved in the response to Sclerotinia sclerotiorum infection. We identified 32 BnLTPs linked to blackleg resistance quantitative trait locus (QTL).

CONCLUSION

The identification and analysis of LTP genes in the B. napus pangenome could help to elucidate the function of BnLTP family members and provide new information for future molecular breeding in B. napus.

Comparison of non-subjective relative fungal biomass measurements to quantify the Leptosphaeria maculans-Brassica napus interaction

BACKGROUND

Blackleg disease, caused by the fungal pathogen Leptosphaeria maculans, is a serious threat to canola (Brassica napus) production worldwide. Quantitative resistance to this disease is a highly desirable trait but is difficult to precisely phenotype. Visual scores can be subjective and are prone to assessor bias. Methods to assess variation in quantitative resistance more accurately were developed based on quantifying in planta fungal biomass, including the Wheat Germ Agglutinin Chitin Assay (WAC), qPCR and ddPCR assays.

RESULTS

Disease assays were conducted by inoculating a range of canola cultivars with L. maculans isolates in glasshouse experiments and assessing fungal biomass in cotyledons, petioles and stem tissue harvested at different timepoints post-inoculation. PCR and WAC assay results were well correlated, repeatable across experiments and host tissues, and able to differentiate fungal biomass in different host-isolate treatments. In addition, the ddPCR assay was shown to differentiate between L. maculans isolates.

CONCLUSIONS

The ddPCR assay is more sensitive in detecting pathogens and more adaptable to high-throughput methods by using robotic systems than the WAC assay. Overall, these methods proved accurate and non-subjective, providing alternatives to visual assessments to quantify the L. maculans-B. napus interaction in all plant tissues throughout the progression of the disease in seedlings and mature plants and have potential for fine-scale blackleg resistance phenotyping in canola.

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 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.

Genetic Analysis of a Horizontal Resistance Locus BLMR2 in Brassica napus

Leptosphaeria maculans causes blackleg disease in Brassica napus. The blackleg disease is mainly controlled by resistance genes in B. napus. Previous studies have shown that the blackleg resistant BLMR2 locus that conferred horizontal resistance under field conditions, is located on chromosome A10 of B. napus. The purpose of this study is to fine map this locus and hence identify a candidate gene underlying horizontal resistance. The spectrum of resistance to L. maculans isolates of the resistance locus BLMR2 was analyzed using near isogenic lines, resistant, and susceptible cultivars. The results showed that this locus was horizontally resistant to all isolates tested. Sequence characterized amplified regions (SCAR), simple sequence repeats (SSR), and single nucleotide polymorphism (SNP) markers were developed in the chromosome region of BLMR2 and a fine genetic map was constructed. Two molecular markers narrowed BLMR2 in a 53.37 kb region where six genes were annotated. Among the six annotated genes, BnaA10g11280D/BnaA10g11290D encoding a cytochrome P450 protein were predicted as the candidate of BLMR2. Based on the profiling of pathogen induced transcriptome, three expressed genes in the six annotated genes were identified while only cytochrome P450 showed upregulation. The candidate corresponds to the gene involved in the indole glucosinolate biosynthesis pathway and plant basal defense in Arabidopsis thaliana. The molecular markers identified in this study will allow the quick incorporation of the BLMR2 allele in rapeseed cultivars to enhance blackleg 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.

Comparison and evolutionary analysis of Brassica nucleotide binding site leucine rich repeat (NLR) genes and importance for disease resistance breeding

The Brassica genus contains many agriculturally significant oilseed and vegetable crops, however the crop yield is threatened by a range of fungal and bacterial pathogens. Nucleotide Binding Site Leucine Rich Repeat (NLR) genes play important roles in plant innate immunity. The evolution of NLR genes is influenced by genomic processes and pathogen selection. At the whole genome level, whole genome duplications (WGDs) generate abundant gene copies, most of which are lost during genome fractionation. At sub-genomic levels, some retained copies undergo duplication forming clusters which facilitate rapid evolution through recombination. The number, distribution and genetic variations of the NLR genes vary among Brassica species and within populations suggesting differential selection pressure exerted by pathogen populations throughout the evolutionary history. A study of the evolution of disease resistance genes in agriculturally important plants such as Brassicas helps gain insights into their function and inform the identification of resistance genes for breeding of resistant lines.

Genome-wide transcriptome reveals mechanisms underlying Rlm1-mediated blackleg resistance on canola

Genetic resistance to blackleg (Leptosphaeria maculans, Lm) of canola (Brassica napus, Bn) has been extensively studied, but the mechanisms underlying the host-pathogen interaction are still not well understood. Here, a comparative transcriptome analysis was performed on a resistant doubled haploid Bn line carrying the resistance gene Rlm1 following inoculation with a virulent (avrLm1) or avirulent (AvrLm1) Lm isolate on cotyledons. A total of 6999 and 3015 differentially expressed genes (DEGs) were identified, respectively, in inoculated local tissues with compatible (susceptible) and incompatible (resistant) interactions. Functional enrichment analysis found several biological processes, including protein targeting to membrane, ribosome and negative regulation of programmed cell death, were over-represented exclusively among up-regulated DEGs in the resistant reaction, whereas significant enrichment of salicylic acid (SA) and jasmonic acid (JA) pathways observed for down-regulated DEGs occurred only in the susceptible reaction. A heat-map analysis showed that both biosynthesis and signaling of SA and JA were induced more significantly in the resistant reaction, implying that a threshold level of SA and JA signaling is required for the activation of Rlm1-mediated resistance. Co-expression network analysis revealed close correlation of a gene module with the resistance, involving DEGs regulating pathogen-associated molecular pattern recognition, JA signaling and transcriptional reprogramming. Substantially fewer DEGs were identified in mock-inoculated (control) cotyledons, relative to those in inoculated local tissues, including those involved in SA pathways potentially contributing to systemic acquired resistance (SAR). Pre-inoculation of cotyledon with either an avirulent or virulent Lm isolate, however, failed to induce SAR on remote tissues of same plant despite elevated SA and PR1 protein. This study provides insights into the molecular mechanism of Rlm1-mediated resistance to blackleg.

Current Status of the Disease-Resistant Gene(s)/QTLs, and Strategies for Improvement in Brassica juncea

Brassica juncea is a major oilseed crop in tropical and subtropical countries, especially in south-east Asia like India, China, Bangladesh, and Pakistan. The widespread cultivation of genetically similar varieties tends to attract fungal pathogens which cause heavy yield losses in the absence of resistant sources. The conventional disease management techniques are often expensive, have limited efficacy, and cause additional harm to the environment. A substantial approach is to identify and use of resistance sources within the Brassica hosts and other non-hosts to ensure sustainable oilseed crop production. In the present review, we discuss six major fungal pathogens of B. juncea: Sclerotinia stem rot (Sclerotinia sclerotiorum), Alternaria blight (Alternaria brassicae), White rust (Albugo candida), Downy mildew (Hyaloperonospora parasitica), Powdery mildew (Erysiphe cruciferarum), and Blackleg (Leptoshaeria maculans). From discussing studies on pathogen prevalence in B. juncea, the review then focuses on highlighting the resistance sources and quantitative trait loci/gene identified so far from Brassicaceae and non-filial sources against these fungal pathogens. The problems in the identification of resistance sources for B. juncea concerning genome complexity in host subpopulation and pathotypes were addressed. Emphasis has been laid on more elaborate and coordinated research to identify and deploy R genes, robust techniques, and research materials. Examples of fully characterized genes conferring resistance have been discussed that can be transformed into B. juncea using advanced genomics tools. Lastly, effective strategies for B. juncea improvement through introgression of novel R genes, development of pre-breeding resistant lines, characterization of pathotypes, and defense-related secondary metabolites have been provided suggesting the plan for the development of resistant B. juncea.

Recent Findings Unravel Genes and Genetic Factors Underlying Leptosphaeria maculans Resistance in Brassica napus and Its Relatives

Among the Brassica oilseeds, canola (Brassica napus) is the most economically significant globally. However, its production can be limited by blackleg disease, caused by the fungal pathogen Lepstosphaeria maculans. The deployment of resistance genes has been implemented as one of the key strategies to manage the disease. Genetic resistance against blackleg comes in two forms: qualitative resistance, controlled by a single, major resistance gene (R gene), and quantitative resistance (QR), controlled by numerous, small effect loci. R-gene-mediated blackleg resistance has been extensively studied, wherein several genomic regions harbouring R genes against L. maculans have been identified and three of these genes were cloned. These studies advance our understanding of the mechanism of R gene and pathogen avirulence (Avr) gene interaction. Notably, these studies revealed a more complex interaction than originally thought. Advances in genomics help unravel these complexities, providing insights into the genes and genetic factors towards improving blackleg resistance. Here, we aim to discuss the existing R-gene-mediated resistance, make a summary of candidate R genes against the disease, and emphasise the role of players involved in the pathogenicity and resistance. The comprehensive result will allow breeders to improve resistance to L. maculans, thereby increasing yield.

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.

Characterization of disease resistance genes in the Brassica napus pangenome reveals significant structural variation

Methods based on single nucleotide polymorphism (SNP), copy number variation (CNV) and presence/absence variation (PAV) discovery provide a valuable resource to study gene structure and evolution. However, as a result of these structural variations, a single reference genome is unable to cover the entire gene content of a species. Therefore, pangenomics analysis is needed to ensure that the genomic diversity within a species is fully represented. Brassica napus is one of the most important oilseed crops in the world and exhibits variability in its resistance genes across different cultivars. Here, we characterized resistance gene distribution across 50 B. napus lines. We identified a total of 1749 resistance gene analogs (RGAs), of which 996 are core and 753 are variable, 368 of which are not present in the reference genome (cv. Darmor-bzh). In addition, a total of 15 318 SNPs were predicted within 1030 of the RGAs. The results showed that core R-genes harbour more SNPs than variable genes. More nucleotide binding site-leucine-rich repeat (NBS-LRR) genes were located in clusters than as singletons, with variable genes more likely to be found in clusters. We identified 106 RGA candidates linked to blackleg resistance quantitative trait locus (QTL). This study provides a better understanding of resistance genes to target for genomics-based improvement and improved disease resistance.

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.

Inheritance Pattern and Molecular Markers for Resistance to Blackleg Disease in Cabbage

The inheritance and causal loci for resistance to blackleg, a devastating disease of Brassicaceous crops, are yet to be known in cabbage (Brassica oleracea L.). Here, we report the pattern of inheritance and linked molecular marker for this trait. A segregating BC₁ population consisting of 253 plants was raised from resistant and susceptible parents, L29 (♀) and L16 (♂), respectively. Cotyledon resistance bioassay of BC₁ population, measured based on a scale of 0-9 at 12 days after inoculation with Leptosphaeria maculans isolate 03-02 s, revealed the segregation of resistance and ratio, indicative of dominant monogenic control of the trait. Investigation of potential polymorphism in the previously identified differentially expressed genes within the collinear region of 'B. napus blackleg resistant loci Rlm1' in B. oleracea identified two insertion/deletion (InDel) mutations in the intron and numerous single nucleotide polymorphisms (SNPs) throughout the LRR-RLK gene Bol040029, of which six SNPs in the first exon caused the loss of two LRR domains in the susceptible line. An InDel marker, BLR-C-InDel based on the InDel mutations, and a high resolution melting (HRM) marker, BLR-C-2808 based on the SNP C²⁸⁰⁸T in the second exon were developed, which predicated the resistance status of the BC₁ population with 80.24%, and of 24 commercial inbred lines with 100% detection accuracy. This is the first report of inheritance and molecular markers linked with blackleg resistance in cabbage. This study will enhance our understanding of the trait, and will be helpful in marker assisted breeding aiming at developing resistant cabbage varieties.

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.

GWAS hints at pleiotropic roles for FLOWERING LOCUS T in flowering time and yield-related traits in canola

Background

Transition to flowering at the right time is critical for local adaptation and to maximize grain yield in crops. Canola is an important oilseed crop with extensive variation in flowering time among varieties. However, our understanding of underlying genes and their role in canola productivity is limited.

Results

We report our analyses of a diverse GWAS panel (300–368 accessions) of canola and identify SNPs that are significantly associated with variation in flowering time and response to photoperiod across multiple locations. We show that several of these associations map in the vicinity of FLOWERING LOCUS T (FT) paralogs and its known transcriptional regulators. Complementary QTL and eQTL mapping studies, conducted in an Australian doubled haploid population, also detected consistent genomic regions close to the FT paralogs associated with flowering time and yield-related traits. FT sequences vary between accessions. Expression levels of FT in plants grown in field (or under controlled environment cabinets) correlated with flowering time. We show that markers linked to the FT paralogs display association with variation in multiple traits including flowering time, plant emergence, shoot biomass and grain yield.

Conclusions

Our findings suggest that FT paralogs not only control flowering time but also modulate yield-related productivity traits in canola.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5964-y) contains supplementary material, which is available to authorized users.

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.

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.

The high-quality genome of Brassica napus cultivar 'ZS11' reveals the introgression history in semi-winter morphotype

Allotetraploid oilseed rape (Brassica napus L.) is an agriculturally important crop. Cultivation and breeding of B. napus by humans has resulted in numerous genetically diverse morphotypes with optimized agronomic traits and ecophysiological adaptation. To further understand the genetic basis of diversification and adaptation, we report a draft genome of an Asian semi-winter oilseed rape cultivar 'ZS11' and its comprehensive genomic comparison with the genomes of the winter-type cultivar 'Darmor-bzh' as well as two progenitors. The integrated BAC-to-BAC and whole-genome shotgun sequencing strategies were effective in the assembly of repetitive regions (especially young long terminal repeats) and resulted in a high-quality genome assembly of B. napus 'ZS11'. Within a short evolutionary period (~6700 years ago), semi-winter-type 'ZS11' and the winter-type 'Darmor-bzh' maintained highly genomic collinearity. Even so, certain genetic differences were also detected in two morphotypes. Relative to 'Darmor-bzh', both two subgenomes of 'ZS11' are closely related to its progenitors, and the 'ZS11' genome harbored several specific segmental homoeologous exchanges (HEs). Furthermore, the semi-winter-type 'ZS11' underwent potential genomic introgressions with B. rapa (A_r ). Some of these genetic differences were associated with key agronomic traits. A key gene of A03.FLC3 regulating vernalization-responsive flowering time in 'ZS11' was first experienced HE, and then underwent genomic introgression event with A_r , which potentially has led to genetic differences in controlling vernalization in the semi-winter types. Our observations improved our understanding of the genetic diversity of different B. napus morphotypes and the cultivation history of semi-winter oilseed rape in Asia.

Leptosphaeria maculans Alters Glucosinolate Profiles in Blackleg Disease-Resistant and -Susceptible Cabbage Lines

Blackleg, a fungal disease caused by Leptosphaeria maculans, is one of the most devastating diseases of Brassica crops worldwide. Despite notable progress elucidating the roles of glucosinolates in pathogen defense, the complex interaction between B. oleracea (cabbage) and L. maculans infection that leads to the selective induction of genes involved in glucosinolate production and subsequent modulation of glucosinolate profiles remains to be fully understood. The current study was designed to identify glucosinolate-biosynthesis genes induced by L. maculans and any associated alterations in glucosinolate profiles to explore their roles in blackleg resistance in 3-month-old cabbage plants. The defense responses of four cabbage lines, two resistant and two susceptible, were investigated using two L. maculans isolates, 03-02 s and 00-100 s. A simultaneous increase in the aliphatic glucosinolates glucoiberverin (GIV) and glucoerucin (GER) and the indolic glucosinolates glucobrassicin (GBS) and neoglucobrassicin (NGBS) was associated with complete resistance. An increase in either aliphatic (GIV) or indolic (GBS and MGBS) glucosinolates was associated with moderate resistance. Indolic glucobrassicin (GBS) and neoglucobrassicin (NGBS) were increased in both resistant and susceptible interactions. Pearson correlation showed positive association between GER content with GSL-OH (Bol033373) expression. Expressions of MYB34 (Bol007760), ST5a (Bol026200), and CYP81F2 (Bol026044) were positively correlated with the contents of both GBS and MGBS. Our results confirm that L. maculans infection induces glucosinolate-biosynthesis genes in cabbage, with concomitant changes in individual glucosinolate contents. In resistant lines, both aliphatic and indolic glucosinolates are associated with resistance, with aliphatic GIV and GER and indolic MGBS glucosinolates particularly important. The association between the genes, the corresponding glucosinolates, and plant resistance broaden our molecular understanding of glucosinolate mediated defense against L. maculans in cabbage.

Single R Gene Introgression Lines for Accurate Dissection of the Brassica - Leptosphaeria Pathosystem

Seven blackleg resistance (R) genes (Rlm1, Rlm2, Rlm3, Rlm4, LepR1, LepR2 & LepR3) were each introgressed into a common susceptible B. napus doubled-haploid (DH) line through reciprocal back-crossing, producing single-R gene introgression lines (ILs) for use in the pathological and molecular study of Brassica-Leptosphaeria interactions. The genomic positions of the R genes were defined through molecular mapping and analysis with transgenic L. maculans isolates was used to confirm the identity of the introgressed genes where possible. Using L. maculans isolates of contrasting avirulence gene (Avr) profiles, we preformed extensive differential pathology for phenotypic comparison of the ILs to other B. napus varieties, demonstrating the ILs can provide for the accurate assessment of Avr-R gene interactions by avoiding non-Avr dependant alterations to resistance responses which can occur in some commonly used B. napus varieties. Whole-genome SNP-based assessment allowed us to define the donor parent introgressions in each IL and provide a strong basis for comparative molecular dissection of the pathosystem.

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.

Multi-environment QTL studies suggest a role for cysteine-rich protein kinase genes in quantitative resistance to blackleg disease in Brassica napus

BACKGROUND

Resistance to the blackleg disease of Brassica napus (canola/oilseed rape), caused by the hemibiotrophic fungal pathogen Leptosphaeria maculans, is determined by both race-specific resistance (R) genes and quantitative resistance loci (QTL), or adult-plant resistance (APR). While the introgression of R genes into breeding material is relatively simple, QTL are often detected sporadically, making them harder to capture in breeding programs. For the effective deployment of APR in crop varieties, resistance QTL need to have a reliable influence on phenotype in multiple environments and be well defined genetically to enable marker-assisted selection (MAS).

RESULTS

Doubled-haploid populations produced from the susceptible B. napus variety Topas and APR varieties AG-Castle and AV-Sapphire were analysed for resistance to blackleg in two locations over 3 and 4 years, respectively. Three stable QTL were detected in each population, with two loci appearing to be common to both APR varieties. Physical delineation of three QTL regions was sufficient to identify candidate defense-related genes, including a cluster of cysteine-rich receptor-like kinases contained within a 49 gene QTL interval on chromosome A01. Individual L. maculans isolates were used to define the physical intervals for the race-specific R genes Rlm3 and Rlm4 and to identify QTL common to both field studies and the cotyledon resistance response.

CONCLUSION

Through multi-environment QTL analysis we have identified and delineated four significant and stable QTL suitable for MAS of quantitative blackleg resistance in B. napus, and identified candidate genes which potentially play a role in quantitative defense responses to L. maculans.

Computational Prediction of Effector Proteins in Fungi: Opportunities and Challenges

Effector proteins are mostly secretory proteins that stimulate plant infection by manipulating the host response. Identifying fungal effector proteins and understanding their function is of great importance in efforts to curb losses to plant diseases. Recent advances in high-throughput sequencing technologies have facilitated the availability of several fungal genomes and 1000s of transcriptomes. As a result, the growing amount of genomic information has provided great opportunities to identify putative effector proteins in different fungal species. There is little consensus over the annotation and functionality of effector proteins, and mostly small secretory proteins are considered as effector proteins, a concept that tends to overestimate the number of proteins involved in a plant-pathogen interaction. With the characterization of Avr genes, criteria for computational prediction of effector proteins are becoming more efficient. There are 100s of tools available for the identification of conserved motifs, signature sequences and structural features in the proteins. Many pipelines and online servers, which combine several tools, are made available to perform genome-wide identification of effector proteins. In this review, available tools and pipelines, their strength and limitations for effective identification of fungal effector proteins are discussed. We also present an exhaustive list of classically secreted proteins along with their key conserved motifs found in 12 common plant pathogens (11 fungi and one oomycete) through an analytical pipeline.

Comparative Transcriptomic Analysis of Virulence Factors in Leptosphaeria maculans during Compatible and Incompatible Interactions with Canola

Leptosphaeria maculans is a hemibiotrophic fungus that causes blackleg of canola (Brassica napus), one of the most devastating diseases of this crop. In the present study, transcriptome profiling of L. maculans was performed in an effort to understand and define the pathogenicity genes that govern both the biotrophic and the necrotrophic phase of the fungus, as well as those that separate a compatible from an incompatible interaction. For this purpose, comparative RNA-seq analyses were performed on L. maculans isolate D5 at four different time points following inoculation on susceptible cultivar Topas-DH16516 or resistant introgression line Topas-Rlm2. Analysis of 1.6 billion Illumina reads readily identified differentially expressed genes that were over represented by candidate secretory effector proteins, CAZymes, and other pathogenicity genes. Comparisons between the compatible and incompatible interactions led to the identification of 28 effector proteins whose chronology and level of expression suggested a role in the establishment and maintenance of biotrophy with the plant. These included all known Avr genes of isolate D5 along with eight newly characterized effectors. In addition, another 15 effector proteins were found to be exclusively expressed during the necrotrophic phase of the fungus, which supports the concept that L. maculans has a separate and distinct arsenal contributing to each phase. As for CAZymes, they were often highly expressed at 3 dpi but with no difference in expression between the compatible and incompatible interactions, indicating that other factors were necessary to determine the outcome of the interaction. However, their significantly higher expression at 11 dpi in the compatible interaction confirmed that they contributed to the necrotrophic phase of the fungus. A notable exception was LysM genes whose high expression was singularly observed on the susceptible host at 7 dpi. In the case of TFs, their higher expression at 7 and 11 dpi on susceptible Topas support an important role in regulating the genes involved in the different pathogenic phases of L. maculans. In conclusion, comparison of the transcriptome of L. maculans during compatible and incompatible interactions has led to the identification of key pathogenicity genes that regulate not only the fate of the interaction but also lifestyle transitions of the fungus.

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.

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.

Co-localisation of the blackleg resistance genes Rlm2 and LepR3 on Brassica napus chromosome A10

BACKGROUND

The protection of canola (Brassica napus) crops against blackleg disease, caused by the fungal pathogen Leptosphaeria maculans, is largely mediated by race-specific resistance genes (R-genes). While many R-genes effective against blackleg disease have been identified in Brassica species, information of the precise genomic locations of the genes is limited.

RESULTS

In this study, the Rlm2 gene for resistance to blackleg, located on chromosome A10 of the B. napus cultivar 'Glacier', was targeted for fine mapping. Molecular markers tightly linked to the gene were developed for use in mapping the resistance locus and defining the physical interval in B. napus. Rlm2 was localised to a 5.8 cM interval corresponding to approximately 873 kb of the B. napus chromosome A10.

CONCLUSION

The recently-cloned B. napus R-gene, LepR3, occupies the same region of A10 as Rlm2 and analysis of the putative B. napus and B. rapa genes in the homologous region identified several additional candidate defense-related genes that may control Rlm2 function.

Gene loss in the fungal canola pathogen Leptosphaeria maculans

Recent comparisons of the increasing number of genome sequences have revealed that variation in gene content is considerably more prevalent than previously thought. This variation is likely to have a pronounced effect on phenotypic diversity and represents a crucial target for the assessment of genomic diversity. Leptosphaeria maculans, a causative agent of phoma stem canker, is the most devastating fungal pathogen of Brassica napus (oilseed rape/canola). A number of L. maculans genes are known to be present in some isolates but lost in the others. We analyse gene content variation within three L. maculans isolates using a hybrid mapping and genome assembly approach and identify genes which are present in one of the isolates but missing in the others. In total, 57 genes are shown to be missing in at least one isolate. The genes encode proteins involved in a range of processes including oxidative processes, DNA maintenance, cell signalling and sexual reproduction. The results demonstrate the effectiveness of the method and provide new insight into genomic diversity in L. maculans.

SNP markers-based map construction and genome-wide linkage analysis in Brassica napus

An Illumina Infinium array comprising 5306 single nucleotide polymorphism (SNP) markers was used to genotype 175 individuals of a doubled haploid population derived from a cross between Skipton and Ag-Spectrum, two Australian cultivars of rapeseed (Brassica napus L.). A genetic linkage map based on 613 SNP and 228 non-SNP (DArT, SSR, SRAP and candidate gene markers) covering 2514.8 cM was constructed and further utilized to identify loci associated with flowering time and resistance to blackleg, a disease caused by the fungus Leptosphaeria maculans. Comparison between genetic map positions of SNP markers and the sequenced Brassica rapa (A) and Brassica oleracea (C) genome scaffolds showed several genomic rearrangements in the B. napus genome. A major locus controlling resistance to L. maculans was identified at both seedling and adult plant stages on chromosome A07. QTL analyses revealed that up to 40.2% of genetic variation for flowering time was accounted for by loci having quantitative effects. Comparative mapping showed Arabidopsis and Brassica flowering genes such as Phytochrome A/D, Flowering Locus C and agamous-Like MADS box gene AGL1 map within marker intervals associated with flowering time in a DH population from Skipton/Ag-Spectrum. Genomic regions associated with flowering time and resistance to L. maculans had several SNP markers mapped within 10 cM. Our results suggest that SNP markers will be suitable for various applications such as trait introgression, comparative mapping and high-resolution mapping of loci in B. napus.

Genome-wide delineation of natural variation for pod shatter resistance in Brassica napus

Resistance to pod shattering (shatter resistance) is a target trait for global rapeseed (canola, Brassica napus L.), improvement programs to minimise grain loss in the mature standing crop, and during windrowing and mechanical harvest. We describe the genetic basis of natural variation for shatter resistance in B. napus and show that several quantitative trait loci (QTL) control this trait. To identify loci underlying shatter resistance, we used a novel genotyping-by-sequencing approach DArT-Seq. QTL analysis detected a total of 12 significant QTL on chromosomes A03, A07, A09, C03, C04, C06, and C08; which jointly account for approximately 57% of the genotypic variation in shatter resistance. Through Genome-Wide Association Studies, we show that a large number of loci, including those that are involved in shattering in Arabidopsis, account for variation in shatter resistance in diverse B. napus germplasm. Our results indicate that genetic diversity for shatter resistance genes in B. napus is limited; many of the genes that might control this trait were not included during the natural creation of this species, or were not retained during the domestication and selection process. We speculate that valuable diversity for this trait was lost during the natural creation of B. napus. To improve shatter resistance, breeders will need to target the introduction of useful alleles especially from genotypes of other related species of Brassica, such as those that we have identified.

Assessing quantitative resistance against Leptosphaeria maculans (phoma stem canker) in Brassica napus (oilseed rape) in young plants

Quantitative resistance against Leptosphaeria maculans in Brassica napus is difficult to assess in young plants due to the long period of symptomless growth of the pathogen from the appearance of leaf lesions to the appearance of canker symptoms on the stem. By using doubled haploid (DH) lines A30 (susceptible) and C119 (with quantitative resistance), quantitative resistance against L. maculans was assessed in young plants in controlled environments at two stages: stage 1, growth of the pathogen along leaf veins/petioles towards the stem by leaf lamina inoculation; stage 2, growth in stem tissues to produce stem canker symptoms by leaf petiole inoculation. Two types of inoculum (ascospores; conidia) and three assessment methods (extent of visible necrosis; symptomless pathogen growth visualised using the GFP reporter gene; amount of pathogen DNA quantified by PCR) were used. In stage 1 assessments, significant differences were observed between lines A30 and C119 in area of leaf lesions, distance grown along veins/petioles assessed by visible necrosis or by viewing GFP and amount of L. maculans DNA in leaf petioles. In stage 2 assessments, significant differences were observed between lines A30 and C119 in severity of stem canker and amount of L. maculans DNA in stem tissues. GFP-labelled L. maculans spread more quickly from the stem cortex to the stem pith in A30 than in C119. Stem canker symptoms were produced more rapidly by using ascospore inoculum than by using conidial inoculum. These results suggest that quantitative resistance against L. maculans in B. napus can be assessed in young plants in controlled conditions. Development of methods to phenotype quantitative resistance against plant pathogens in young plants in controlled environments will help identification of stable quantitative resistance for control of crop diseases.

Identifying genetic diversity of avirulence genes in Leptosphaeria maculans using whole genome sequencing

Next generation sequencing technology allows rapid re-sequencing of individuals, as well as the discovery of single nucleotide polymorphisms (SNPs), for genomic diversity and evolutionary analyses. By sequencing two isolates of the fungal plant pathogen Leptosphaeria maculans, the causal agent of blackleg disease in Brassica crops, we have generated a resource of over 76 million sequence reads aligned to the reference genome. We identified over 21,000 SNPs with an overall SNP frequency of one SNP every 2,065 bp. Sequence validation of a selection of these SNPs in additional isolates collected throughout Australia indicates a high degree of polymorphism in the Australian population. In preliminary phylogenetic analysis, isolates from Western Australia clustered together and those collected from Brassica juncea stubble were identical. These SNPs provide a novel marker resource to study the genetic diversity of this pathogen. We demonstrate that re-sequencing provides a method of validating previously characterised SNPs and analysing differences in important genes, such as the disease related avirulence genes of L. maculans. Understanding the genetic characteristics of this devastating pathogen is vital in developing long-term solutions to managing blackleg disease in Brassica crops.

A consensus map of rapeseed (Brassica napus L.) based on diversity array technology markers: applications in genetic dissection of qualitative and quantitative traits

Background

Dense consensus genetic maps based on high-throughput genotyping platforms are valuable for making genetic gains in Brassica napus through quantitative trait locus identification, efficient predictive molecular breeding, and map-based gene cloning. This report describes the construction of the first B. napus consensus map consisting of a 1,359 anchored array based genotyping platform; Diversity Arrays Technology (DArT), and non-DArT markers from six populations originating from Australia, Canada, China and Europe. We aligned the B. napus DArT sequences with genomic scaffolds from Brassica rapa and Brassica oleracea, and identified DArT loci that showed linkage with qualitative and quantitative loci associated with agronomic traits.

Results

The integrated consensus map covered a total of 1,987.2 cM and represented all 19 chromosomes of the A and C genomes, with an average map density of one marker per 1.46 cM, corresponding to approximately 0.88 Mbp of the haploid genome. Through in silico physical mapping 2,457 out of 3,072 (80%) DArT clones were assigned to the genomic scaffolds of B. rapa (A genome) and B. oleracea (C genome). These were used to orientate the genetic consensus map with the chromosomal sequences. The DArT markers showed linkage with previously identified non-DArT markers associated with qualitative and quantitative trait loci for plant architecture, phenological components, seed and oil quality attributes, boron efficiency, sucrose transport, male sterility, and race-specific resistance to blackleg disease.

Conclusions

The DArT markers provide increased marker density across the B. napus genome. Most of the DArT markers represented on the current array were sequenced and aligned with the B. rapa and B. oleracea genomes, providing insight into the Brassica A and C genomes. This information can be utilised for comparative genomics and genomic evolution studies. In summary, this consensus map can be used to (i) integrate new generation markers such as SNP arrays and next generation sequencing data; (ii) anchor physical maps to facilitate assembly of B. napus genome sequences; and (iii) identify candidate genes underlying natural genetic variation for traits of interest.

The Brassica napus blackleg resistance gene LepR3 encodes a receptor-like protein triggered by the Leptosphaeria maculans effector AVRLM1

LepR3, found in the Brassica napus cv 'Surpass 400', provides race-specific resistance to the fungal pathogen Leptosphaeria maculans, which was overcome after great devastation in Australia in 2004. We investigated the LepR3 locus to identify the genetic basis of this resistance interaction. We employed a map-based cloning strategy, exploiting collinearity with the Arabidopsis thaliana and Brassica rapa genomes to enrich the map and locate a candidate gene. We also investigated the interaction of LepR3 with the L. maculans avirulence gene AvrLm1 using transgenics. LepR3 was found to encode a receptor-like protein (RLP). We also demonstrated that avirulence towards LepR3 is conferred by AvrLm1, which is responsible for both the Rlm1 and LepR3-dependent resistance responses in B. napus. LepR3 is the first functional B. napus disease resistance gene to be cloned. AvrLm1's interaction with two independent resistance loci, Rlm1 and LepR3, highlights the need to consider redundant phenotypes in 'gene-for-gene' interactions and offers an explanation as to why LepR3 was overcome so rapidly in parts of Australia.

Genetic and physical mapping of flowering time loci in canola (Brassica napus L.)

We identified quantitative trait loci (QTL) underlying variation for flowering time in a doubled haploid (DH) population of vernalisation-responsive canola (Brassica napus L.) cultivars Skipton and Ag-Spectrum and aligned them with physical map positions of predicted flowering genes from the Brassica rapa genome. Significant genetic variation in flowering time and response to vernalisation were observed among the DH lines from Skipton/Ag-Spectrum. A molecular linkage map was generated comprising 674 simple sequence repeat, sequence-related amplified polymorphism, sequence characterised amplified region, Diversity Array Technology, and candidate gene based markers loci. QTL analysis indicated that flowering time is a complex trait and is controlled by at least 20 loci, localised on ten different chromosomes. These loci each accounted for between 2.4 and 28.6% of the total genotypic variation for first flowering and response to vernalisation. However, identification of consistent QTL was found to be dependant upon growing environments. We compared the locations of QTL with the physical positions of predicted flowering time genes located on the sequenced genome of B. rapa. Some QTL associated with flowering time on A02, A03, A07, and C06 may represent homologues of known flowering time genes in Arabidopsis; VERNALISATION INSENSITIVE 3, APETALA1, CAULIFLOWER, FLOWERING LOCUS C, FLOWERING LOCUS T, CURLY LEAF, SHORT VEGETATIVE PHASE, GA3 OXIDASE, and LEAFY. Identification of the chromosomal location and effect of the genes influencing flowering time may hasten the development of canola varieties having an optimal time for flowering in target environments such as for low rainfall areas, via marker-assisted selection.

Comparative analysis of FLC homologues in Brassicaceae provides insight into their role in the evolution of oilseed rape

We identified nine FLOWERING LOCUS C homologues (BnFLC) in Brassica napus and found that the coding sequences of all BnFLCs were relatively conserved but the intronic and promoter regions were more divergent. The BnFLC homologues were mapped to six of 19 chromosomes. All of the BnFLC homologues were located in the collinear region of FLC in the Arabidopsis genome except BnFLC.A3b and BnFLC.C3b, which were mapped to noncollinear regions of chromosome A3 and C3, respectively. Four of the homologues were associated significantly with quantitative trait loci for flowering time in two mapping populations. The BnFLC homologues showed distinct expression patterns in vegetative and reproductive organs, and at different developmental stages. BnFLC.A3b was differentially expressed between the winter-type and semi-winter-type cultivars. Microsynteny analysis indicated that BnFLC.A3b might have been translocated to the present segment in a cluster with other flowering-time regulators, such as a homologue of FRIGIDA in Arabidopsis. This cluster of flowering-time genes might have conferred a selective advantage to Brassica species in terms of increased adaptability to diverse environments during their evolution and domestication process.

Identification and characterization of candidate Rlm4 blackleg resistance genes in Brassica napus using next-generation sequencing

A thorough understanding of the relationships between plants and pathogens is essential if we are to continue to meet the agricultural needs of the world's growing population. The identification of genes underlying important quantitative trait loci is extremely challenging in complex genomes such as Brassica napus (canola, oilseed rape or rapeseed). However, recent advances in next-generation sequencing (NGS) enable much quicker identification of candidate genes for traits of interest. Here, we demonstrate this with the identification of candidate disease resistance genes from B. napus for its most devastating fungal pathogen, Leptosphaeria maculans (blackleg fungus). These two species are locked in an evolutionary arms race whereby a gene-for-gene interaction confers either resistance or susceptibility in the plant depending on the genotype of the plant and pathogen. Preliminary analysis of the complete genome sequence of Brassica rapa, the diploid progenitor of B. napus, identified numerous candidate genes with disease resistance characteristics, several of which were clustered around a region syntenic with a major locus (Rlm4) for blackleg resistance on A7 of B. napus. Molecular analyses of the candidate genes using B. napus NGS data are presented, and the difficulties associated with identifying functional gene copies within the highly duplicated Brassica genome are discussed.