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.     

Seed transmission of turnip yellow mosaic virus in winter turnip and winter oilseed rapes

This is the first record of seed transmission of turnip yellow mosaic virus (TYMV) in oilseed and turnip rapes. The seed transmission of TYMV in a naturally infected winter turnip rape (Brassica napus var. silvestris) cultivar Perko PVH was investigated. By ELISA 1.6%, 3.2% and 8.3% seed transmission of the virus was found in seed of plants from three localities. The proportion of infected seeds produced by artificially infected plants of winter oilseed rape (Brassica napus ssp. oleifera) and winter turnip rape cultivars was determined. The virus transmission rate, expressed as the proportion of virus-infected plants which germinated from the seed was for the oilseed rape cvs Jet Neuf 0.1%, Solida 0.4%, Silesia 0.8%, Darmor 1.2%, SL-507 0.2%, SL-509 0.0% and for the winter turnip rape cv. Perko 1.5%. ELISA cannot be used in direct tests on bulk seed lots to estimate proportion of infected seed, but must be used on germinated seedlings.     

Direct nitrous oxide emissions from oilseed rape cropping – a meta‐analysis

Oilseed rape is one of the leading feedstocks for biofuel production in Europe. The climate change mitigation effect of rape methyl ester (RME) is particularly challenged by the greenhouse gas (GHG) emissions during crop production, mainly as nitrous oxide (N₂O) from soils. Oilseed rape requires high nitrogen fertilization and crop residues are rich in nitrogen, both potentially causing enhanced N₂O emissions. However, GHG emissions of oilseed rape production are often estimated using emission factors that account for crop‐type specifics only with respect to crop residues. This meta‐analysis therefore aimed to assess annual N₂O emissions from winter oilseed rape, to compare them to those of cereals and to explore the underlying reasons for differences. For the identification of the most important factors, linear mixed effects models were fitted with 43 N₂O emission data points deriving from 12 different field sites. N₂O emissions increased exponentially with N‐fertilization rates, but interyear and site‐specific variability were high and climate variables or soil parameters did not improve the prediction model. Annual N₂O emissions from winter oilseed rape were 22% higher than those from winter cereals fertilized at the same rate. At a common fertilization rate of 200 kg N ha^?1 yr^?1, the mean fraction of fertilizer N that was lost as N₂O‐N was 1.27% for oilseed rape compared to 1.04% for cereals. The risk of high yield‐scaled N₂O emissions increased after a critical N surplus of about 80 kg N ha^?1 yr^?1. The difference in N₂O emissions between oilseed rape and cereal cultivation was especially high after harvest due to the high N contents in oilseed rape's crop residues. However, annual N₂O emissions of winter oilseed rape were still lower than predicted by the Stehfest and Bouwman model. Hence, the assignment of oilseed rape to the crop‐type classes of cereals or other crops should be reconsidered.     

Sequence‐level comparative analysis of the Brassica napus genome around two stearoyl‐ACP desaturase loci

We conducted a sequence‐level comparative analyses, at the scale of complete bacterial artificial chromosome (BAC) clones, between the genome of the most economically important Brassica species, Brassica napus (oilseed rape), and those of Brassica rapa, the genome of which is currently being sequenced, and Arabidopsis thaliana. We constructed a new B. napus BAC library and identified and sequenced clones that contain homoeologous regions of the genome including stearoyl‐ACP desaturase‐encoding genes. We sequenced the orthologous region of the genome of B. rapa and conducted comparative analyses between the Brassica sequences and those of the orthologous region of the genome of A. thaliana. The proportion of genes conserved (～56%) is lower than has been reported previously between A. thaliana and Brassica (～66%). The gene models for sets of conserved genes were used to determine the extent of nucleotide conservation of coding regions. This was found to be 84.2 ± 3.9% and 85.8 ± 3.7% between the B. napus A and C genomes, respectively, and that of A. thaliana, which is consistent with previous results for other Brassica species, and 97.5 ± 3.1% between the B. napus A genome and B. rapa, and 93.1 ± 4.9% between the B. napus C genome and B. rapa. The divergence of the B. napus genes from the A genome and the B. rapa genes was greater than anticipated and indicates that the A genome ancestor of the B. napus cultivar studied was relatively distantly related to the cultivar of B. rapa selected for genome sequencing.     

Origins of the amphiploid species <Emphasis Type="Italic">Brassica napus</Emphasis>L. investigated by chloroplast and nuclear molecular markers

Background  

The amphiploid species Brassica napus (oilseed rape, Canola) is a globally important oil crop yielding food, biofuels and industrial compounds such as lubricants and surfactants. Identification of the likely ancestors of each of the two genomes (designated A and C) found in B. napus would facilitate incorporation of novel alleles from the wider Brassica genepool in oilseed rape crop genetic improvement programmes. Knowledge of the closest extant relatives of the genotypes involved in the initial formation of B. napus would also allow further investigation of the genetic factors required for the formation of a stable amphiploid and permit the more efficient creation of fully fertile re-synthesised B. napus. We have used a combination of chloroplast and nuclear genetic markers to investigate the closest extant relatives of the original maternal progenitors of B. napus. This was based on a comprehensive sampling of the relevant genepools, including 83 accessions of A genome B. rapa L. (both wild and cultivated types), 94 accessions of B. napus and 181 accessions of C genome wild and cultivated B. oleracea L. and related species.     

Genetic diversity of oilseedBrassica napus germ plasm based on restriction fragment length polymorphisms

Oilseed rape (Brassica napus) is an important oilseed crop worldwide. Cultivars have been developed for many growing regions, however little is known about genetic diversity inB. napus germ plasm. The purpose of the research presented here was to study the genetic diversity and relationships ofB. napus accessions using restriction fragment length polymorphisms (RFLPs). Eighty threeB. napus accessions were screened using 43 genomic DNA clones which revealed 161 polymorphic fragments. Each accession was uniquely identified by the markers with the exception of the near-isogenic cvs Triton and Tower. The RFLP data were analyzed by cluster analysis of similarity coefficients and by principal component analysis. Overall, there were three major groups of cultivars. The first group included only spring accessions, the second mostly winter accessions and the third, rutabagas and oilseed rape accessions from China and Japan. These results indicate that withinB. napus, winter and spring cultivars represent genetically distinct groups. The grouping of accessions by cluster analysis was generally consistent with known pedigrees. This consistency included the grouping of lines derived both by backcrossing or self-pollination with their parents.     

Analysis of leaf appearance,leaf death and phoma leaf spot,caused by Leptosphaeria maculans,on oilseed rape (Brassica napus) cultivars

Development of phoma leaf spot (caused by Leptosphaeria maculans) on winter oilseed rape (canola, Brassica napus) was assessed in two experiments at Rothamsted in successive years (2003–04 and 2004–05 growing seasons). Both experiments compared oilseed rape cultivars Eurol, Darmor, Canberra and Lipton, which differ in their resistance to L. maculans. Data were analysed to describe disease development in terms of increasing numbers of leaves affected over thermal time from sowing. The cultivars showed similar patterns of leaf spot development in the 2003–04 experiment when inoculum concentration was relatively low (up to 133 ascospores m⁻³ air), Darmor developing 5.3 diseased leaves per plant by 5 May 2004, Canberra 6.6, Eurol 6.8 and Lipton 7.5. Inoculum concentration was up to sevenfold greater in 2004–05, with Eurol and Darmor developing 2.4 diseased leaves per plant by 16 February 2005, whereas Lipton and Canberra developed 2.8 and 3.0 diseased leaves, respectively. Based on three defined periods of crop development, a piece-wise linear statistical model was applied to the progress of the leaf spot disease (cumulative diseased leaves) in relation to appearance (‘birth’) and death of leaves for individual plants of each cultivar. Estimates of the thermal time from sowing until appearance of the first leaf or death of the first leaf, the rate of increase in number of diseased leaves and the area under the disease progress line (AUDPL) for the first time period were made. In 2004–05, Canberra (1025 leaves ×°C days) and Lipton (879) had greater AUDPL values than Eurol (427) and Darmor (598). For Darmor and Lipton, the severity of leaf spotting could be related to the severity of stem canker at harvest. Eurol had less leaf spotting but severe stem canker, whereas Canberra had more leaf spotting but less severe canker.     

Development of iFOX‐hunting as a functional genomic tool and demonstration of its use to identify early senescence‐related genes in the polyploid Brassica napus

Functional genomic studies of many polyploid crops, including rapeseed (Brassica napus), are constrained by limited tool sets. Here we report development of a gain‐of‐function platform, termed ‘iFOX (inducible Full‐length cDNA OvereXpressor gene)‐Hunting’, for inducible expression of B. napus seed cDNAs in Arabidopsis. A Gateway‐compatible plant gene expression vector containing a methoxyfenozide‐inducible constitutive promoter for transgene expression was developed. This vector was used for cloning of random cDNAs from developing B. napus seeds and subsequent Agrobacterium‐mediated transformation of Arabidopsis. The inducible promoter of this vector enabled identification of genes upon induction that are otherwise lethal when constitutively overexpressed and to control developmental timing of transgene expression. Evaluation of a subset of the resulting ~6000 Arabidopsis transformants revealed a high percentage of lines with full‐length B. napus transgene insertions. Upon induction, numerous iFOX lines with visible phenotypes were identified, including one that displayed early leaf senescence. Phenotypic analysis of this line (rsl‐1327) after methoxyfenozide induction indicated high degree of leaf chlorosis. The integrated B. napuscDNA was identified as a homolog of an Arabidopsis acyl‐CoA binding protein (ACBP) gene designated BnACBP1‐like. The early senescence phenotype conferred by BnACBP1‐like was confirmed by constitutive expression of this gene in Arabidopsis and B. napus. Use of the inducible promoter in the iFOX line coupled with RNA‐Seq analyses allowed mechanistic clues and a working model for the phenotype associated with BnACBP1‐like expression. Our results demonstrate the utility of iFOX‐Hunting as a tool for gene discovery and functional characterization of Brassica napus genome.     

Reconstituting the genome of a young allopolyploid crop,Brassica napus,with its related species

Brassica napus (AⁿAⁿCⁿCⁿ) is an important worldwide oilseed crop, but it is a young allotetraploid with a short evolutionary history and limited genetic diversity. To significantly broaden its genetic diversity and create a novel heterotic population for sustainable rapeseed breeding, this study reconstituted the genome of B. napus by replacing it with the subgenomes from 122 accessions of Brassica rapa (A^rA^r) and 74 accessions of Brassica carinata (B^cB^cC^cC^c) and developing a novel gene pool of B. napus through five rounds of extensive recurrent selection. When compared with traditional B. napus using SSR markers and high‐throughput SNP/Indel markers through genotyping by sequencing, the newly developed gene pool and its homozygous progenies exhibited a large genetic distance, rich allelic diversity, new alleles and exotic allelic introgression across all 19 AC chromosomes. In addition to the abundant genomic variation detected in the AC genome, we also detected considerable introgression from the eight chromosomes of the B genome. Extensive trait variation and some genetic improvements were present from the early recurrent selection to later generations. This novel gene pool produced equally rich phenotypic variation and should be valuable for rapeseed genetic improvement. By reconstituting the genome of B. napus by introducing subgenomic variation within and between the related species using intense selection and recombination, the whole genome could be substantially reorganized. These results serve as an example of the manipulation of the genome of a young allopolyploid and provide insights into its rapid genome evolution affected by interspecific and intraspecific crosses.     

Patterns of molecular variation in a species-wide germplasm set of <Emphasis Type="Italic">Brassica napus</Emphasis>

Rapeseed (Brassica napus L.) is the leading European oilseed crop serving as source for edible oil and renewable energy. The objectives of our study were to (i) examine the population structure of a large and diverse set of B. napus inbred lines, (ii) investigate patterns of genetic diversity within and among different germplasm types, (iii) compare the two genomes of B. napus with regard to genetic diversity, and (iv) assess the extent of linkage disequilibrium (LD) between simple sequence repeat (SSR) markers. Our study was based on 509 B. napus inbred lines genotyped with 89 genome-specific SSR primer combinations. Both a principal coordinate analysis and software STRUCTURE revealed that winter types, spring types, and swedes were assigned to three major clusters. The genetic diversity of winter oilseed rape was lower than the diversity found in other germplasm types. Within winter oilseed rape types, a decay of genetic diversity with more recent release dates and reduced levels of erucic acid and glucosinolates was observed. The percentage of linked SSR loci pairs in significant (r ² > Q _{95 unlinked loci pairs}) LD was 6.29% for the entire germplasm set. Furthermore, LD decayed rapidly with distance, which will allow a relatively high mapping resolution in genome-wide association studies using our germplasm set, but, on the other hand, will require a high number of markers.     

Comparative mitochondrial genome analysis reveals the evolutionary rearrangement mechanism in Brassica

The genus Brassica has many species that are important for oil, vegetable and other food products. Three mitochondrial genome types (mitotype) originated from its common ancestor. In this paper, a B. nigra mitochondrial main circle genome with 232,407 bp was generated through de novo assembly. Synteny analysis showed that the mitochondrial genomes of B. rapa and B. oleracea had a better syntenic relationship than B. nigra. Principal components analysis and development of a phylogenetic tree indicated maternal ancestors of three allotetraploid species in Us triangle of Brassica. Diversified mitotypes were found in allotetraploid B. napus, in which napus‐type B. napus was derived from B. oleracea, while polima‐type B. napus was inherited from B. rapa. In addition, the mitochondrial genome of napus‐type B. napus was closer to botrytis‐type than capitata‐type B. oleracea. The sub‐stoichiometric shifting of several mitochondrial genes suggested that mitochondrial genome rearrangement underwent evolutionary selection during domestication and/or plant breeding. Our findings clarify the role of diploid species in the maternal origin of allotetraploid species in Brassica and suggest the possibility of breeding selection of the mitochondrial genome.     

Molecular and phenotypic characterization of near isogenic lines at QTL for quantitative resistance to Leptosphaeria maculans in oilseed rape (Brassica napus L.)

The most common and effective way to control phoma stem canker (blackleg) caused by Leptosphaeria maculans in oilseed rape (Brassica napus) is by breeding resistant cultivars. Specific resistance genes have been identified in B. napus and related species but in some B. napus cultivars resistance is polygenic [mediated by quantitative trait loci (QTL)], postulated to be race non-specific and durable. The genetic basis of quantitative resistance in the French winter oilseed rape ‘Darmor’, which was derived from ‘Jet Neuf’, was previously examined in two genetic backgrounds. Stable QTL involved in blackleg resistance across year and genetic backgrounds were identified. In this study, near isogenic lines (NILs) were produced in the susceptible background ‘Yudal’ for four of these QTL using marker-assisted selection. Various strategies were used to develop new molecular markers, which were mapped in these QTL regions. These were used to characterize the length and homozygosity of the ‘Darmor-bzh’ introgressed segment in the NILs. Individuals from each NIL were evaluated in blackleg disease field trials and assessed for their level of stem canker in comparison to the recurrent line ‘Yudal’. The effect of QTL LmA2 was clearly validated and to a lesser extent, QTL LmA9 also showed an effect on the disease level. This work provides valuable material that can be used to study the mode of action of genetic factors involved in L. maculans quantitative resistance.     

Effects of non‐host plant odour on Meligethes aeneus during immigration to oilseed rape

The use of semiochemicals for the manipulation of the pollen beetle, Meliethes aeneus (Fabricius) (Coleoptera: Nitidulidae), is being investigated for potential incorporation into a push‐pull strategy for this pest, which damages oilseed rape, Brassica napus L. (Brassicaceae), throughout Europe. Previous laboratory behavioural studies using volatiles from non‐host plants showed that M. aeneus is repelled by the odour of lavender, Lavendula angustifolia Mill. (Lamiaceae), essential oil. This article reports on semi‐field and field trials to investigate this behaviour under more realistic conditions. Semi‐field experiments were conducted to assess the relative importance of olfaction at different points in host location behaviour by M. aeneus. The results showed that oilseed rape plants treated with lavender odour were less colonised by M. aeneus in comparison with an untreated control, but that the treatment effect was much reduced if the lavender odour was applied after colonisation. The field experiment demonstrated that lavender odour caused a significant reduction in the number of adult M. aeneus infesting the oilseed rape plants in the treatment plots compared to the control plots. Overall, these findings are very encouraging for the future development of a push‐pull pest control system.     

Genome of ‘Charleston Gray’, the principal American watermelon cultivar,and genetic characterization of 1,365 accessions in the U.S. National Plant Germplasm System watermelon collection

Years of selection for desirable fruit quality traits in dessert watermelon (Citrullus lanatus) has resulted in a narrow genetic base in modern cultivars. Development of novel genomic and genetic resources offers great potential to expand genetic diversity and improve important traits in watermelon. Here, we report a high‐quality genome sequence of watermelon cultivar ‘Charleston Gray’, a principal American dessert watermelon, to complement the existing reference genome from ‘97103’, an East Asian cultivar. Comparative analyses between genomes of ‘Charleston Gray’ and ‘97103’ revealed genomic variants that may underlie phenotypic differences between the two cultivars. We then genotyped 1365 watermelon plant introduction (PI) lines maintained at the U.S. National Plant Germplasm System using genotyping‐by‐sequencing (GBS). These PI lines were collected throughout the world and belong to three Citrullus species, C. lanatus, C. mucosospermus and C. amarus. Approximately 25 000 high‐quality single nucleotide polymorphisms (SNPs) were derived from the GBS data using the ‘Charleston Gray’ genome as the reference. Population genomic analyses using these SNPs discovered a close relationship between C. lanatus and C. mucosospermus and identified four major groups in these two species correlated to their geographic locations. Citrullus amarus was found to have a distinct genetic makeup compared to C. lanatus and C. mucosospermus. The SNPs also enabled identification of genomic regions associated with important fruit quality and disease resistance traits through genome‐wide association studies. The high‐quality ‘Charleston Gray’ genome and the genotyping data of this large collection of watermelon accessions provide valuable resources for facilitating watermelon research, breeding and improvement.     

Genetic changes in a novel breeding population of Brassica napus synthesized from hundreds of crosses between B. rapa and B. carinata

Introgression of genomic variation between and within related crop species is a significant evolutionary approach for population differentiation, genome reorganization and trait improvement. Using the Illumina Infinium Brassica 60K SNP array, we investigated genomic changes in a panel of advanced generation new‐type Brassica napus breeding lines developed from hundreds of interspecific crosses between 122 Brassica rapa and 74 Brassica carinata accessions, and compared them with representative accessions of their three parental species. The new‐type B. napus population presented rich genetic diversity and abundant novel genomic alterations, consisting of introgressions from B. rapa and B. carinata, novel allelic combinations, reconstructed linkage disequilibrium patterns and haplotype blocks, and frequent deletions and duplications (nonrandomly distributed), particularly in the C subgenome. After a much shorter, but very intensive, selection history compared to traditional B. napus, a total of 15 genomic regions with strong selective sweeps and 112 genomic regions with putative signals of selective sweeps were identified. Some of these regions were associated with important agronomic traits that were selected for during the breeding process, while others were potentially associated with restoration of genome stability and fertility after interspecific hybridization. Our results demonstrate how a novel method for population‐based crop genetic improvement can lead to rapid adaptation, restoration of genome stability and positive responses to artificial selection.     

Genome‐wide selection footprints and deleterious variations in young Asian allotetraploid rapeseed

Brassica napus (AACC, 2n = 38) is an important oilseed crop grown worldwide. However, little is known about the population evolution of this species, the genomic difference between its major genetic groups, such as European and Asian rapeseed, and the impacts of historical large‐scale introgression events on this young tetraploid. In this study, we reported the de novo assembly of the genome sequences of an Asian rapeseed (B. napus), Ningyou 7, and its four progenitors and compared these genomes with other available genomic data from diverse European and Asian cultivars. Our results showed that Asian rapeseed originally derived from European rapeseed but subsequently significantly diverged, with rapid genome differentiation after hybridization and intensive local selective breeding. The first historical introgression of B. rapa dramatically broadened the allelic pool but decreased the deleterious variations of Asian rapeseed. The second historical introgression of the double‐low traits of European rapeseed (canola) has reshaped Asian rapeseed into two groups (double‐low and double‐high), accompanied by an increase in genetic load in the double‐low group. This study demonstrates distinctive genomic footprints and deleterious SNP (single nucleotide polymorphism) variants for local adaptation by recent intra‐ and interspecies introgression events and provides novel insights for understanding the rapid genome evolution of a young allopolyploid crop.     

Oilseed rape: learning about ancient and recent polyploid evolution from a recent crop species

Oilseed rape (Brassica napus) is one of our youngest crop species, arising several times under cultivation in the last few thousand years and completely unknown in the wild. Oilseed rape originated from hybridisation events between progenitor diploid species B. rapa and B. oleracea, both important vegetable species. The diploid progenitors are also ancient polyploids, with remnants of two previous polyploidisation events evident in the triplicated genome structure. This history of polyploid evolution and human agricultural selection makes B. napus an excellent model with which to investigate processes of genomic evolution and selection in polyploid crops. The ease of de novo interspecific hybridisation, responsiveness to tissue culture, and the close relationship of oilseed rape to the model plant Arabidopsis thaliana, coupled with the recent availability of reference genome sequences and suites of molecular cytogenetic and high‐throughput genotyping tools, allow detailed dissection of genetic, genomic and phenotypic interactions in this crop. In this review we discuss the past and present uses of B. napus as a model for polyploid speciation and evolution in crop species, along with current and developing analysis tools and resources. We further outline unanswered questions that may now be tractable to investigation.     

Evaluation of Linkage Disequilibrium Pattern and Association Study on Seed Oil Content in Brassica napus Using ddRAD Sequencing

High-density genetic markers are the prerequisite for understanding linkage disequilibrium (LD) and genome-wide association studies (GWASs) of complex traits in crops. To evaluate the LD pattern in oilseed rape, we sequenced a previous association panel containing 189 B. napus inbred lines using double-digested restriction-site associated DNA (ddRAD) and genotyped 19,327 RAD tags. A total of 15,921 RAD tags were assigned to a published genetic linkage map and the majority (71.1%) of these tags was uniquely mapped to the draft reference genome “Darmor-bzh.” The distance of LD decay was 1,214 kb across the genome at the background level (r² = 0.26), with the distances of LD decay being 405 kb and 2,111 kb in the A and C subgenomes, respectively. A total of 361 haplotype blocks with length > 100 kb were identified in the entire genome. The association panel could be classified into two groups, P1 and P2, which are essentially consistent with the geographical origins of varieties. A large number of group-specific haplotypes were identified, reflecting that varieties in the P1 and P2 groups experienced distinct selection in breeding programs to adapt their different growth habitats. GWAS repeatedly detected two loci significantly associated with oil content of seeds based on the developed SNPs, suggesting that the high-density SNPs were useful for understanding the genetic determinants of complex traits in GWAS.