Skip to main content
  • Research article
  • Open access
  • Published:

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

Abstract

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.

Background

Rapeseed (also known as canola and oilseed rape, Brassica napus L., 2n = 4× = 38; genomes AACC) is the second largest oilseed crop after soybean and provides a valuable rotational crop for farmers in many parts of the world. During the past 15 years, global oilseed rape production has doubled to >60 million tonnes (http://faostat.fao.org/; data sourced April 2012), due to the high demand for healthy vegetable oil, feedstock and a renewable source for the biodiesel industry. This amphidiploid Brassica species most likely originated as a result of repeated natural hybridization and genome doubling between the monogenomic diploid species Brassica rapa (2n = 2× = 20, genome AA) and Brassica oleracea (2n = 2× = 18, genome CC) along the Mediterranean coastline in Southern Europe, and was probably selected as an oilseed crop only 300–400 years ago [1–3].

Selection of useful variation in crop plants has been a major thrust of early farmers since the dawn of agriculture. In modern Brassica breeding programs, classical genetic analysis and molecular genetic approaches have been used to improve our understanding of the inheritance of various qualitative and quantitative loci and causative genes, to estimate their number, position and genetic effects, and to identify DNA-based markers associated with traits of agronomic importance [4].

Several linkage maps have been constructed from recombination data in B. napus mapping populations, which depict the distances between loci, as well as their order on a chromosome. These linkage maps have been based on a range of marker systems such as restriction fragment length polymorphisms (RFLPs), randomly amplified polymorphic DNAs (RAPD), amplified fragment length polymorphisms (AFLPs), simple sequence repeats (SSRs), sequence- tagged sites (STSs), sequence-related amplified polymorphisms (SRAPs) and single nucleotide polymorphisms (SNPs) [5–14]. These maps have been used to uncover qualitative and quantitative trait loci (QTLs) underlying traits of agricultural importance, for linkage disequilibrium assessment, candidate gene identification, marker-assisted selection, and genetic evolutionary studies [14–19]. The number of marker loci in these genetic linkage maps varied from 219 to 13,551, where the marker density was dependent upon the level of polymorphism between the parental lines of mapping populations and type of marker system used for detecting polymorphisms.

In order to increase marker density for predictive molecular breeding and map-based cloning of genes, dense consensus genetic maps have been developed for various crops including B. napus[14, 20, 21]. Such consensus maps are developed by mapping several common markers from each linkage group, in different mapping populations, and anchoring the linkage maps with consensus markers. However, most B. napus genetic maps have been based either upon low-throughput marker assays or on different marker systems, rendering them difficult to utilise to compare the map locations of markers and therefore their linkages with traits of interest across populations.

Diversity Array Technology (DArT) markers do not require prior DNA sequence information and can simultaneously genotype a large number of SNPs and insertion/deletion polymorphisms across the genome in a single assay, providing a low-cost option for genotyping [22]. The DArT marker system was first used in rice, and has since been extensively utilised for genetic applications such as genetic diversity analysis, construction of genetic linkage maps, and QTL/linkage analysis in different crops [23, 24]. Recently, a 3072-clone B. napus DArT array was used to generate a linkage map using a B. napus doubled haploid (DH) population derived from Lynx-/Monty-028DH [13]. This information was used to identify QTLs for flowering time, resistance to blackleg disease caused by Leptosphaeria maculans, and resistance to pod shatter [25–27].

Comparative mapping studies have shown extensive co-linearity between the genomes of Brassica and Arabidopsis, which diverged from a common ancestor approximately 20 MYA [28–30]. The genome of B. rapa has recently been sequenced [31], whilst the genomes of B. oleracea and B. napus are currently being annotated and analysed and are expected to become publicly available in 2013. These advances in structural and functional genomics have provided us with an opportunity to align the genetic DArT markers with the sequenced genomes of B. rapa and B. oleracea. This allows the identification of markers for marker-assisted selection in rapeseed improvement programs.

In the present study, we constructed a consensus map of rapeseed based on DArT marker datasets produced using six DH B. napus populations. These populations were derived from Ag-Castle/Topas (AT), BLN2762/Surpass400 (BS), Lynx-037DH/Monty-028DH (LM), Maxol*1/Westar-10 (MW), Skipton/Ag-Spectrum (SAS), and Tapidor/Ningyou7 (TN). The parents of the mapping populations originated from Australia, Canada, China and Europe, and are currently being used internationally in rapeseed germplasm enhancement programs. We physically mapped the B. napus DArT sequences using genome sequences from scaffolds of B. rapa and B. oleracea. DArT marker loci associated with agronomic traits, including various components of phenology, plant architecture, seed and oil quality, nutrient uptake and mobilisation, and resistance to blackleg disease and abiotic stresses were also identified.

Results

Map construction from individual mapping populations

Individual component maps were constructed from six populations that were genotyped with DArT, along with a selection of SSR markers (Table 1, Additional file 1). The number of DArT marker loci on each map ranged from 217 (AT population) to 437 (LM population). The TN map consisted of the most markers (971); 403 DArT and 568 non-DArT (381 SSRs, 69 STSs including candidate genes, 68 SNPs, 34 RFLPs, 14 centromeric sequence-related and two AFLP markers). The genetic linkage map length for the six populations varied from 759.9 to 2,288 cM, with a mean length of 1446.5 cM.

Table 1 The distribution and density of markers in the linkage maps of individual doubled haploid populations

The order of the DArT marker loci was generally similar among the different maps (Additional file 2). However, there was some evidence of chromosomal segmental rearrangements, such as inversions or translocations, on some genomic regions on chromosomes A7, A9, A10, C1, C2, and C6. For example, XbrPb-658333 and XbrPb-659113 were mapped within 0.5 cM on chromosome A10 in the TN population; however, these markers were not linked in the SAS population and were located on different linkage groups (A10-II and A10-III on the SAS map). In the LM population, a high frequency of homoeologous recombination was observed between chromosomes A7 and C6 [13]. The map lengths and marker densities were positively correlated among the populations (r = 0.73).

Segregation distortion was much more pronounced on chromosome A7 in the populations derived from AT (original data not shown), MW and SAS [32, 33], where either targeted selection for desirable alleles (for example for different blackleg resistance genes Rlm1, Rlm3 and Rlm4) might have occurred within the breeding programs or could have been a low level of homoeologous interaction with chromosome C6.

The DArT markers were generally well distributed throughout the genome, although some chromosomes exhibited either no or low polymorphism for the DArT markers (Additional file 2). For example, chromosome A3 had low DArT marker density (4 DArT/35 total markers) in the TN map. A number of DArT markers co-segregated and therefore mapped on as the same locus. DArT marker redundancy (markers that map to the identical map position as their respective neighbouring markers) ranged from 9.2% in the TN to 47.5% in the AT mapping population (Table 1). In order to estimate redundancy of the DArT markers at the sequence level, we aligned the DArT sequences using ClustalW program and found that the majority of co-segregating DArT clones (0.01 cM = ~6 kb) had high levels (77% to 99.4%) of sequence identities (Additional file 3). However, a marker cluster on chromosome C3 had lower sequence identities (53 to 80%) and hence may not represent the same locus. Co-segregation of markers may have occurred because recombination is difficult to observe between loci that are very close to each other, particularly in the small populations used here (Table 2).

Table 2 Attributes of doubled haploid genetic mapping populations used for consensus map construction

Consensus map construction

Individual component maps were aligned to the reference linkage map of the TN-DH population to construct a consensus map. A total of 1359 marker loci [791 DArTs, 381 SSRs, 69 STSs, 68 SNPs, 34 RFLPs, 14 centromeric region specific (CS) and two AFLPs] were integrated in the consensus map, covering 1987.2 cM, with an average interval distance of 1.46 cM, across the 19 chromosomes of B. napus (Table 3, Additional file 4). The A and C genomes had genetic lengths of 1062.5 cM and 924.6 cM, with the mean distances between neighbouring markers of 1.08 cM and 2.48 cM, respectively. Approximately 2.6 times more markers were mapped on the A genome (986) than the C genome (373), suggesting that the A genome was more polymorphic than the C genome in B. napus.

Table 3 The distribution and density of markers across 19 rapeseed chromosomes of the consensus map

The length of individual chromosomes varied from 59.5 cM (chromosome C7) to 153.9 cM (chromosome C3), whereas the number of marker loci per chromosome ranged from 23 (chromosome C7) to 156 (chromosome A3). Chromosomes A5 and A9 had the highest marker density (1 marker/0.89 cM), and chromosome C5 had the lowest marker density (1 marker/3.91 cM). Both map length and marker density across chromosomes were moderately correlated (r2 = 0.56, Additional file 5).

In general, the consensus map had good coverage of DArT and non-DArT markers, however there were 16 genomic sites (accounting for 6% of the genetic map) where adjacent markers were >10 cM apart: on chromosomes A1, A2, A7, A10, C2, C4, C5, C6, C8 and C9 (Table 3, Additional file 4). Of the 791 genome-wide DArT markers, 22% showed redundancy, whereas only 0.9% of non-DArT markers showed redundancy. Chromosomes A6 and C1 exhibited the maximum redundancy of 3.7% and 40.9% for non-DArT and DArT markers, respectively (Table 3). The overall marker density on the consensus map (1 marker/1.46 cM) was higher compared to that observed in individual populations (1 marker/1.89 - 4.20 cM) (Tables 1 and 3). A total of 10 major clusters of DArT markers (≥5 loci/0.01 cM) were identified on chromosomes A1, A2, A3, A6, A9, and C1 (Additional file 4).

Comparison of individual maps and the consensus map

The congruency of marker order and map positions was consistent with an earlier map of TN [9]. The positions of the DArT markers on the consensus map were generally consistent with their locations on the six individual linkage maps from the biparental crosses (Additional file 2). The order of markers was moderately correlated in most pair-wise comparisons of populations, and highly correlated between the AT and SAS populations (Additional file 6). However, there was a poor correlation between the BS and LM populations (r < 0.4). Due to the diverse backgrounds of the DH populations used in this study, the number of shared DArT markers was extremely variable and ranged from 0 to 27 per chromosome. The largest number of shared markers per chromosome (27) was between the LM and BS maps of chromosome A3, covering a genetic distance of ~123 cM. Spearman’s rank correlation of the marker orders also varied between different component maps (Additional file 6). Therefore, it was difficult to reliably compare the marker orders across individual maps.

Colinearity of the consensus map and the sequence assemblies of A/C ancestral species

The genomic locations of DArT markers and their homology (sequence identities and bit scores) are provided in Additional files 7 and 8. Approximately 80% (2,457/3,072) of DArT markers were physically mapped to the sequenced genomes of B. rapa and B. oleracea. This enabled evaluation of the correspondence of the DArT markers on the genetic and physical maps. The location of 615 DArT clones on the genome scaffolds of B. rapa and B. oleracea (256 on the A genome and 359 on the C genome) could not be determined, as they did not return any hit. These sequences may represent unassembled regions of the current genome assemblies or alternatively may be divergent B. napus sequences (Additional file 9). The alignment of the integrated genetic and physical maps was generally in agreement (Additional file 10); however, in a number of instances marker order as determined by sequence homology suggested chromosomal rearrangements on the linkage groups A3, A9, C2 and C4 and or due to inaccuracies in mapping. For two linkage groups, C7 and C8, the alignment of the DArT markers to the genome scaffolds suggested that significant portions of these two chromosomes were not polymorphic in the six crosses.

Distribution of marker duplications across A and C genomes

Eighty-three markers (18 DArT and 65 non-DArT) showed inter- and intra-genomic duplicated loci on the homoeologous chromosomes of the A and C genomes (Additional files 4 and 11). The frequency of duplicated (multi-locus) markers varied from 2.27% (18/791) to 11.4% (65/568) for DArT and non-DArT markers, respectively, in the consensus map. Chromosome A9 had the most duplicated marker loci (13), followed by A2 (10), whereas C1 and C5 did not contain any duplicated marker loci. In addition, up to five alleles for the marker Xpw123 could be mapped on chromosomes A6, A9 and C8. Several DArT markers also showed multiple sequence identities with the sequences of B. rapa and B. oleracea localised on different chromosomes (Additional files 7 and 8).

Association between DArT markers and agricultural traits

We surveyed all genetically mapped markers, including DArT, that showed association with traits of agricultural importance in the mapping populations. This included traits from the published literature, along with fatty acid quality QTLs, published for the first time here. In total 139 marker intervals were identified, associated with resistance to Sclerotinia and L. maculans (causing blackleg disease), yield related traits (seed number, seed weight, seed yield, oil yield, branch number, and pod number), oil quality traits (glucosinolate content, fatty acid content, and vitamin E), biomass yield, phenology traits (plant height, flowering time, vernalisation requirement, maturity time), male sterility, boron efficiency, and uptake of nutrients, including boron, calcium, iron, magnesium, copper, zinc, and phosphorus (Additional files 4 and 11) [9, 27, 32–47]. We did not attempt to locate all (246) QTLs that were reported to be involved in the glucosinolate metabolic pathway in the TN population [35]. Details of four newly identified QTLs governing 16:0, 18:1, 18:2 and 18:3 fatty acids along with five correlated QTLs for 20:1, 20:2, 22:0, 22:1 and 24:0 in the LM population are presented in Table 4. These highly significant QTLs accounted for 13.9% to 75.6% of the explained variation. The marker interval XbrPb-X657955 - dFAD2a explained the highest (75.6%) variation, whilst the marker interval brPb-662948 - brPb-660893 (13.9%) explained the lowest amount of variation in linolenic acid and palmitic acid content, respectively. The fatty acid desaturase-2 (FAD2) gene based marker was mapped within two of the QTLs on chromosome A5. A majority of the qualitative and quantitative trait related loci identified in this study, resided within 10 cM of DArT marker intervals. Some previously published QTLs could not be precisely identified due to the lack of common shared markers between the consensus map and the original mapping population and hence they were defined as approximate (Additional file 12).

Table 4 Newly identified QTL for fatty acid seed content in the LM population

Discussion

In this study, we developed the first consensus genetic map of the rapeseed genome based on DArT markers. The development of this consensus map will provide a platform to compare chromosomal locations of markers across rapeseed populations and facilitate identification of simple and complex inherited trait-marker associations, comprehensive assessment of genetic diversity, and whole genome selection in rapeseed breeding programs. The consensus map consisted of 1,359 markers spanning all 19 B. napus chromosomes, with an average marker density of one marker per 1.46 cM. This corresponds to approximately 880 Mbp of the B. napus genome, 94% of the mapped AC genome was contained within intervals <10 cM. These observations suggest that this consensus map will be suitable for various applications including detection of quantitative traits in rapeseed improvement programs, as in most QTL analyses a 10-cM interval between marker loci is commonly used for regression analysis. This consensus map based on the sequenced DArT markers will also allow the positional cloning of the causative genes controlling phenotypic variation, estimation of linkage disequilibrium at the individual chromosome/genome level, and identification of genomic rearrangements such as translocations that occurred during the period between the genome triplication of Brassica and the divergence of B. rapa and B. oleracea and their distant relative model plant, Arabidopsis.

Different marker systems were used in each population, for example, the LM population was mapped with DArT markers, it had the highest density of SSR markers. These SSR markers (developed by the Agiculture and Agri-Food Canada Consortium) were not used in any other population undertaken in this study except for the TN population. The TN population was also mapped with a variety of markers based on SNPs and CS, which were not used in any other population. Furthermore, DArT and SSR markers amplify intra- and inter-homoalleles as a result of chromosomal duplications between homoeologous sequences among subgenomes and between paralogous duplicated sequences and other rearrangements such as in the TN (Additional files 4 and 11) and SAS populations [35, 48], which makes it difficult to compare individual component maps. Of the DArT clones that were useful for mapping across the six segregating populations in this study, relatively few loci were common between populations. This could be due to utilisation of a range of DH populations derived from genetically diverse parental lines used in order to map the DArT clones across the B. napus genome. For example, the cultivar Surpass 400 is derived from several backcrosses into cultivated B. napus from B. rapa ssp. sylvestris. The introgression of chromosomal segments from B. rapa into Ningyou7 – one of the parental lines of the TNDH population - has also been recently reported [49]. We could not ascertain the map positions of some DArT loci or their order on certain chromosomes (for example in the MW population, LG I to VI, Additional file 1) due to the absence of shared markers among mapping populations. However, marker order was generally consistent within some chromosomes where more than four shared markers were present for assessment of statistical significance [21] and there was a global conservation of co-linearity between the corresponding linkage groups of integrated maps and the A (B. rapa) and C (B. oleracea) genome scaffolds. On some chromosomes no polymorphism was found between the parental lines of the DH populations, for example on chromosome C8 in the populations MW and BS. These findings suggest that the DArT sequences in the genomic representations were conserved.

The results presented here demonstrate that the consensus map was longer (1,987.2 cM) than the six individual maps with the exception of the LM population. The LM map contained a higher proportion of distorted and distinct loci [13], than most of the populations. The recombination frequency can also vary with genetic background [48]. In addition, scoring errors tend to be higher with SSR markers which show extensive stuttering and multiple bands (alleles), as compared to DArT markers. These biological and non-biological factors could have skewed the map estimates.

In general, the DArT markers were well-distributed across the genome, however, certain chromosomal regions showed extensive clustering such as on chromosomes A9 and C2 (Additional file 4). This may indicate the presence of gene-rich regions and uneven distribution of recombination events along chromosomes, or it may suggest that the DArT markers are preferentially surveying DNA polymorphisms that are unevenly distributed along chromosomes. Similar results have been reported in other crops including in B. napus[14, 49, 50]. Uneven distribution of DArT markers between the A and C genomes may indicate genome evolution events of B. napus from its diploid progenitor species or uneven representation of DArT clones from source species used in the array development. The DArT markers derived from genomic representations are known to introduce some degree of redundancy [22]. However, in this study, genomic representations were prepared using the CNG methylation sensitive restriction enzyme PstI, which generates low and single copy DNA fractions in plants.

In comparison to the non-DArT markers (11.4%), a low proportion (2.27%) of multi-locus DArT markers was found. Hybridisation-based markers (such as DArT) select against multi-locus markers, because hybridisation intensities contributed by different loci are difficult to resolve in the DArT allele calling process and such markers are normally scored as monomorphic. A similar frequency (~1.8%) of multi-locus DArT markers has been reported in rye, sorghum, and barley [22, 50, 51]. Existence of multi-locus markers as a result of DNA (DArT) hybridisation from orthologous and non-homologous regions of rapeseed, and multiple hits between DArT sequences and the sequenced scaffolds of B. rapa and B. oleracea, provide strong evidence for the presence of intra- and inter-chromosomal duplicated loci in B. napus. Previous studies have reported that most of the approximately 1.2 Gb genome of B. napus[52] comprises genome sequences from the two progenitor species, which exhibit significant co-linearity. However, homoeologous recombination plays a major role in chromosome rearrangements, such as duplications and reciprocal translocations [53–57]. In certain cases, the physical position of DArT markers on the B. rapa and B. oleracea scaffolds did not correspond with their genetic positions on the linkage maps of segregating populations as reported previously [58]. This could have been caused by large translocations within the mapping populations, mapping inaccuracies and/or errors in assembling genome scaffolds. However, there was a good correspondence between genetic ordering and the current genome sequence assemblies. Approximately 20% of the DArT marker sequences did not have a significant match to B. rapa and B. oleracea. This may be due to the incomplete genome coverage in the current genomic scaffolds or some of them may represent novel loci that may have evolved in B. napus as a result of inter-specific recombination between B. rapa and B. oleracea. However, a more likely explanation is biological (structural) variation among ancestral genotypes involved in the development of B. napus cultivars used for mapping, as compared to B. rapa and B. oleracea genotypes used for sequencing the A and C genomes [58]. It is also possible that the reference accessions used to build current A and C genome scaffold assemblies lack some genomic regions that are present in the wider gene-pools.

Discrepancies such as inversions and translocations were identified when comparing the consensus map and the individual component maps, as have been reported previously in B. napus[59] and in comparisons of the A genome in B. napus and B. rapa[58]. Apart from genome rearrangements, this could be due to the smaller size of the populations used here [60], and variation in the number of recombination events in different regions of the plant genomes [61].

A recent study has shown good co-linearity of DArT and non-DArT markers between the genetic linkage map of the SAS population and the assembled pseudomolecules of B. rapa[27]. As the sequenced genome scaffolds of other A and C genome Brassica species become available in the public domain, DArT sequences (available on http://www.diversityarrays.com/) can be used for comparative genomic analysis, identification of QTLs for predictive breeding (as shown in this study) and for the identification of candidate genes linked to traits of agricultural importance.

In this study, we related DArT markers with the majority of QTLs that have been identified in six breeding populations, however this can be extended to other traits and populations of B. napus, as well as to related A and C genome Brassica species. For example, several DArT markers were related to the positions of QTLs associated with flowering and quality components in the TN and SAS populations [9, 15, 27, 37, 40]. Two QTLs for flowering time were mapped in the genomic regions harbouring homoeologues of FLOWERING TIME LOCUS C of Arabidopsis, BnFLC.A3 on chromosomes A3 and Bn.FLC.A10 on chromosome A10 [15, 27, 40]. In the LM population, we identified four QTLs associated with 16:0, 18:1, 18:2 and 18:3 fatty acids on chromosomes A1, A5 and C4 (Table 4). Among them, two major QTLs associated with 18:1 (oleic acid) and 18:2 (linoleic acid) were mapped on A5 within a region containing FAD2 gene (Table 4). Localisation of these QTLs suggested that FAD2 controls variation for both oleic acid and linoleic acid content in this population. The FAD2 of the endoplasmic reticulum encodes ω-6 desaturases which is responsible for conversion of oleic acid to linoleic acid by inserting a double bond at the ω-6 position [62].

The parents of mapping populations used in this study such as Ag-Castle, Maxol, Monty, Skipton, Ag-Spectrum, and Surpass 400 are known to harbour resistance both at the seedling (race-specific, qualitative) and adult plant stages (race-specific and/or race-non-specific) to L. maculans[26, 32, 63]. DArT markers were mapped in the vicinity of a genomic region associated with Rlm1, Rlm3, and Rlm4 genes for resistance to L. maculans on the chromosome A7 [32, 33]. It was interesting to note that some of the loci associated with correlated traits such as seed yield, pod number, seed weight, flowering time, seed number and plant height were localised on the same genomic regions (A1, A2, A3, A7, C3 and C6) on the consensus map, which may have been due to pleiotropic effects of a QTL on different traits and/or tight linkage between multiple traits and QTLs [64].

Conclusions

We constructed a high-density consensus map of B. napus utilising 1,359 DArT and non-DArT based markers. This consensus map was useful in locating loci associated with various traits of agronomic importance. The sequences of the DArT clones are publicly available (http://www.diversityarrays.com/) and provide a valuable resource in predicted genomic breeding, map-based gene cloning, and comparative analysis studies of A and C genomes Brassica species. The development of the integrated consensus map described in conjunction with the physical locations of B. napus DArT markers on the genome scaffolds, and the identification of molecular markers flanking genomic regions associated with agronomic traits, will empower rapeseed breeding programs to identify candidate/causative genes controlling genetic variation for such traits. This will enhance selection efficiency, especially of quantitative traits governed by a large number of genes and influenced with G × E interaction via marker-assisted selection.

Methods

Mapping populations

Six DH populations derived from 12 parents were used for consensus map construction. The structure of these populations and the traits that were studied are summarised in Table 2. We used the TN-DH population from parents Tapidor (Winter-type, European cultivar) and Ningyou7 (semi-winter, Chinese cultivar) as a reference, because this population has been extensively used by the international rapeseed community for the genetic/comparative mapping of a range of traits of agronomic importance such as plant architecture, flowering time, seed and oil quality attributes, boron efficiency, and glucosinolate content [9, 15, 34–38]. The TN-DH population has been used to map 614 non-DArT markers based upon SSRs, SNP, STSs, single strand conformation polymorphisms (SSCPs), RFLPs, cleaved amplified polymorphic sequences (CAPSs), and AFLPs, as described previously [9, 35, 39] and with SNPs identified from the rapeseed transcriptome [49]. This population has resulted in the discovery and mapping of candidate genes for sucrose transporter, α-tocopherol, fatty acid elongase, indehiscent (IND) gene for pod shatter, and FLOWERING LOCUS C (FLC), FLOWERING LOCUS T, APETALA2 for flowering time [15, 34, 40–42, 65]. The TN-DH published map was further saturated with DArT markers in this study.

We used the published map data of Lynx-037DH/Monty-028DH (LM) population [13]. The LM map was constructed utilising DArTs, SSRs, intron polymorphism (IP), and candidate gene-based markers for fatty acid desaturase and FLC genes. Other populations derived from AT, BS, MW, and SAS were genotyped predominantly with DArT markers and their individual component maps were developed in this study.

Trait phenotyping

Seed fatty acid quality (specifically 14:0, 16:0, 16:1, 18:0, 18:1, 18:2, 18:3, 20:0, 20:1, 20:2, 22:0, 22:1, 24:0 and 24:1 fatty acids) was measured in the LM DH population using gas chromatography (GC; Perkin-Elmer Auto System Gas Chromatograph; Waltham, MA, USA) using 1.5 g of seed per line as described previously [43]. DH lines were grown as spaced plants in two partially replicated spaced plant nurseries at Manjimup Horticultural Research Centre (Manjimup, WA, Australia) in the summer or 2004/05 and at Shenton Park Field Station (Perth, WA, Australia) in the winter of 2005, with 12 DH lines included in both nurseries. Self-pollination was enforced by enclosing racemes in perforated plastic bread bags.

DNA isolation

DNA was isolated from the fresh leaf samples of 2-3-week-old seedlings from parental lines and their progenies using a standard phenol/chloroform extraction method [66].

DArT marker analyses

All six populations were genotyped with the B. napus version 1.0 DArT microarray containing 3,072 markers that were selected after genotyping several hundreds of B. napus accessions at Diversity Arrays Technology Pty Ltd (DArT P/L, Yarralumla, ACT, Australia) as described previously [13]. The PstI/BstNI genomic representations of individual samples (parental lines and their segregating derivatives) were generated, labelled with fluorescent dyes (Cy3 and Cy5) by random priming [67], and hybridised with the B. napus DArT array. Images of microarrays were acquired using a scanner (Tecan LS300; Grodig, Salzburg, Austria) and further analysed with the DArTsoft software version 7.4.7 (DArT P/L). The same software was used to score polymorphisms among parental lines and their DH lines in a binary format (for the presence of marker in the representation as ‘1’ and for the absence as ‘0’, as described previously [13, 22]. For quality control, two measures were used: the first was based on the Q value (a quality parameter measuring bimodality of signal distribution between ‘0’ and ‘1’ clusters) and the second was based on call rate (P, the percentage of DNA samples with defined ‘0’ or ‘1’ allele calls). Only high quality DArT clones with Q > 77%, a call rate >97% and 100% allele-calling consistency (reproducibility) across the technical replicates were selected as markers for genetic mapping of different DH populations. However, some lower quality (those with less well supported map positions) were integrated as attached markers into the most appropriate positions in framework component maps, as described previously in the LM population [13].

PCR-based marker analysis

Both publicly and privately available SSR primer-pairs with prefixes ‘BRMS’, ‘BRAS’, ‘CB’, ‘MR’, and ‘MD’ were obtained from the literature [5, 68, 69]. Two assays; capillary electrophoresis and agarose gel electrophoresis were performed for allele sizing of SSR amplicons generated with fluorescently labelled and non-fluorescent primers, respectively. The 5’ ends of the forward primers were tailed with 19-bp long M13 sequences and labelled with a fluorescent dye (D2, D3 or D4; Beckman Coulter Inc., Fullerton, USA) as previously described [64, 70]. SSR primers were synthesised by Sigma-Aldrich (Castle Hill, NSW, Australia) and were diluted to 1 to 10 pmol/μl depending upon the assay used. PCRs and amplifications were carried out as described previously [64]. Amplified DNA fragments (>600 bp) were separated on a CEQ8000 DNA sequencer (Beckman Coulter Inc., Brea, CA, USA) and their sizes measured using fragment analysis software. Fragments were scored by allele sizes in bp and then converted into binary format as described above. We adopted the standard nomenclature for the DArT markers as described previously [13].

Linkage map construction

Individual genetic (component) maps of the six populations were constructed separately using DArT linkage group and marker ordering software [71]. Prior to construction of the map, the redundant markers were ‘binned’, and within each bin the marker with the highest quality was used for map construction.

The DArT marker ordering system followed a three-step process similar to the one described previously [72]. In the first step, the ‘first’ markers were grouped. A complete graph of all the markers was constructed where the weight assigned to the edges corresponded to the recombination frequency. Edges corresponding to a recombination frequency above a threshold (typically around 0.2) were removed and the resultant connected components were collected into linkage groups. In the second step, the markers in each group were ordered. The optimum marker order was posed as the travelling salesman path within the group [73]. The optimal solution was based on the Concorde solver. Finally, the marker positions were converted to centiMorgans (cM) by applying the Kosambi function [74] to the recombination frequencies.

The reference map of TN was compared with the map constructed using software JoinMap 4.0 [75] with a threshold LOD score of 3.0, however the minimum LOD score of 1.4 was used to incorporate some markers. Non-DArT based markers (SSRs, SNPs, STSs, SSCPs, RFLPs, CAPSs, and AFLPs) with known chromosomal locations on the TN genetic linkage map were used to assign linkage groups to chromosomes A1-A10 and C1-C9 representing AA and CC genomes, respectively. Subsequently, this reference map was used to assign the chromosomal positions for DArT markers in other described mapping populations.

Allele segregation ratios for each marker locus were determined using χ2 tests to determine whether they conformed to expected Mendelian ratios for one locus (1:1 allelic ratio) or two locus (1:3 or 3:1 allelic ratio) models [13]. Markers that showed significant segregation distortion (P < 0.01) were discarded and were not used for linkage mapping.

Consensus map construction

The segregation marker data and their order determined for each population were integrated into the consensus map with the DArT consensus map software [71]. The software required the following inputs: (i) seed (reference) map with one linkage group and marker positions per chromosome; and (ii) set of linkage groups from individual populations with marker positions and chromosome assignment for each group. The construction procedure used the following pseudo-code algorithm: (i) initialised consensus map with the seed map; and (ii) for each chromosome, we found a subset of linkage groups for this chromosome and repeated until subset was empty for each group. The process required at least three markers in common with the consensus map and correlated the positions of common markers with the consensus map. The group with the highest commonality [correlation × log (number of common markers)] was identified, and if the correlation was larger than 0.5, all markers were joined to the consensus map by linear interpolation. The group was removed from the subset and the process repeated. The individual and consensus maps were visualised graphically using the software MapChart [76].

Sequence analysis of DArT clones

All 3,072 DArT markers that represented the current B. napus DArT array (version 1) were sequenced using ABI Prism Big Dye Terminator Cycle Sequencing Ready Reaction Kit (Perkin Elmer, MA, USA) and analysed on an ABI DNA Sequencer. The clones were sequenced using universal M13 forward primer. The sequences were analysed with the dedicated pipeline developed by Diversity Arrays Technology Pty Ltd (unpublished).

In silico mapping of DArT markers on the B. rapa and B. oleraceagenomes

Sequences of DArT clones were used to search the assembled B. rapa[31] and B. oleracea (Isobel Parkin and Andrew Sharpe, planned for publication) genome scaffolds using BLASTN [77] and the predicted marker positions were compared with their locations on the consensus genetic map. The vector and adapter sequences were trimmed from DNA sequences before determining the percentage sequence identities. We recorded the percentage identity using an E-value 10-5 threshold of all scoring pairs rather than recording only the best hit, in order to capture genome-wide sequence similarities.

Detection of DArT sequence redundancy

In order to estimate redundancy at the sequence level, we aligned the DArT sequences with a multiple alignment algorithm CLUSTALW [78] implemented in Geneious software [79], with the default option for analysis of gaps. Although it is difficult to establish a threshold to declare a DArT marker redundant, we used a pragmatic approach and assumed that markers showing <80% sequence homology were less likely to be identical, because they should not hybridise to the same targets under stringent hybridisation conditions.

Identifying trait-marker associations and chromosomal (marker) rearrangements

The consensus map developed in this study was used to (i) localise QTLs and major genes that were identified in the six DH populations; and (ii) identify chromosomal rearrangements such as duplications and inversions on the basis of markers. The location of genetic markers on the consensus map was delimited by identifying published common markers flanking the qualitative and quantitative loci in the individual genetic mapping population. QTL analysis of seed fatty acid content in the LM population was performed using MultiQTL v2.5 (MultiQTL Ltd, Haifa, Israel) using default parameters. The goodness of fit of one-gene models was compared to those of two-gene models.

References

  1. Song K, Osborne TC: Polyphyletic origins of Brassica napus: new evidence based on organelle and nuclear RFLP analyses. Genome. 1992, 35: 992-1001. 10.1139/g92-152.

    Article  Google Scholar 

  2. Gómez-Campo C, Prakash S: Origin and domestication. Biology of Brassica coenospecies. Edited by: Gómez-Campo C. 1999, Netherlands: Elsevier, 33-58.

    Chapter  Google Scholar 

  3. U N: Genomic analysis in Brassica with special reference to the experimental formation of B. napus and peculiar mode of fertilisation. Jpn J Bot. 1935, 7: 389-452.

    Google Scholar 

  4. Snowdon RJ, Friedt W: Molecular markers in Brassica oilseed breeding: current status and future possibilities. Plant Breeding. 2004, 123 (1): 1-8. 10.1111/j.1439-0523.2003.00968.x.

    Article  CAS  Google Scholar 

  5. Piquemal J, Cinquin E, Couton F, Rondeau C, Seignoret E, Doucet I, Perret D, Villeger MJ, Vincourt P, Blanchard P: Construction of an oilseed rape (Brassica napus L.) genetic map with SSR markers. Theor Appl Genet. 2005, 111 (8): 1514-1523. 10.1007/s00122-005-0080-6.

    Article  CAS  PubMed  Google Scholar 

  6. Ferreira ME, Williams PH, Osborn TC: RFLP mapping of Brassica napus using doubled haploid lines. Theor Appl Genet. 1994, 89 (5): 615-621.

    Article  CAS  PubMed  Google Scholar 

  7. Parkin IAP, Sharpe AG, Keith DJ, Lydiate DJ: Identification of the A and C genomes of the amphidiploid Brassica napus (oilseed rape). Genome. 1995, 38: 1122-1131. 10.1139/g95-149.

    Article  CAS  PubMed  Google Scholar 

  8. Uzunova M, Ecke W, WeiBleder K, Robbelen G: Mapping the genome of rapeseed (Brassica napus L). I. Construction of an RFLP linkage map and localisation of QTLs for seed glucosinolate content. Theor Appl Genet. 1995, 90: 194-204.

    Article  CAS  PubMed  Google Scholar 

  9. Qiu D, Morgan C, Shi J, Long Y, Liu J, Li R, Zhuang X, Wang Y, Tan X, Dietrich E: A comparative linkage map of oilseed rape and its use for QTL analysis of seed oil and erucic acid content. Theor Appl Genet. 2006, 114: 67-80. 10.1007/s00122-006-0411-2.

    Article  CAS  PubMed  Google Scholar 

  10. Cheng X, Xu J, Xia S, Gu J, Yang Y, Fu J, Qian X, Zhang S, Wu J, Liu K: Development and genetic mapping of microsatellite markers from genome survey sequences in Brassica napus. Theor Appl Genet. 2009, 118 (6): 1121-1131. 10.1007/s00122-009-0967-8.

    Article  CAS  PubMed  Google Scholar 

  11. Li G, Quiros CF: Sequence-related amplified polymorphism (SRAP), a new marker system based on a simple PCR reaction: its application to mapping and gene tagging in Brassica. Theor Appl Genet. 2001, 103 (2): 455-461. 10.1007/s001220100570.

    Article  CAS  Google Scholar 

  12. Sun Z, Wang Z, Tu J, Zhang J, Yu F, McVetty P, Li G: An ultradense genetic recombination map for Brassica napus consisting of 13551 SRAP markers. Theor Appl Genet. 2007, 114 (8): 1305-1317. 10.1007/s00122-006-0483-z.

    Article  CAS  PubMed  Google Scholar 

  13. Raman H, Raman R, Nelson MN, Aslam MN, Rajasekaran R, Wratten N, Cowling WA, Kilian A, Sharpe AG, Schondelmaier J: Diversity Array Technology markers: Genetic diversity analyses and linkage map construction in rapeseed (Brassica napus L.). DNA Res. 2012, 19 (1): 51-65. 10.1093/dnares/dsr041.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  14. Delourme R, Falentin C, Fomeju B, Boillot M, Lassalle G, Andre I, Duarte J, Gauthier V, Lucante N, Marty A, Pauchon M: High-density SNP-based genetic map development and linkage disequilibrium assessment in Brassica napus L. BMC Genomics. 2013, 14: 120-10.1186/1471-2164-14-120.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  15. Long Y, Shi J, Qiu D, Li R, Zhang C, Wang J, Hou J, Zhao J, Shi L, Park B-S: Flowering time quantitative trait loci analysis of oilseed brassica in multiple environments and genomewide alignment with Arabidopsis. Genetics. 2007, 177 (4): 2433-2444.

    PubMed Central  CAS  PubMed  Google Scholar 

  16. Ecke W, Uzunova M, Weißleder K: Mapping the genome of rapeseed (Brassica napus L.). II. Localization of genes controlling erucic acid synthesis and seed oil content. Theor Appl Genet. 1995, 91 (6): 972-977.

    CAS  PubMed  Google Scholar 

  17. Hasan M, Seyis F, Badani A, Pons-Kühnemann J, Friedt W, Lühs W, Snowdon R: Analysis of genetic diversity in the Brassica napus L. gene pool using SSR markers. Genet Resour Crop Evol. 2006, 53: 793-802. 10.1007/s10722-004-5541-2.

    Article  CAS  Google Scholar 

  18. Hu X, Sullivan-Gilbert M, Gupta M, Thompson S: Mapping of the loci controlling oleic and linolenic acid contents and development of fad2 and fad3 allele-specific markers in canola (Brassica napus L.). Theor Appl Genet. 2006, 113 (3): 497-507. 10.1007/s00122-006-0315-1.

    Article  CAS  PubMed  Google Scholar 

  19. Pilet ML, Delourme R, Foisset N, Renard M: Identification of loci contributing to quantitative field resistance to blackleg disease, causal agent Leptosphaeria maculans (Desm.) Ces. et de Not., in winter rapeseed (Brassica napus L.). Theor Appl Genet. 1998, 96: 23-30. 10.1007/s001220050704.

    Article  Google Scholar 

  20. Lombard V, Delourme R: A consensus linkage map for rapeseed (Brassica napus L.): construction and integration of three individual maps from DH populations. Theor Appl Genet. 2001, 103 (4): 491-507. 10.1007/s001220100560.

    Article  CAS  Google Scholar 

  21. Wang J, Lydiate D, Parkin I, Falentin C, Delourme R, Carion P, King G: Integration of linkage maps for the amphidiploid Brassica napus and comparative mapping with Arabidopsis and Brassica rapa. BMC Genomics. 2011, 12: 101-10.1186/1471-2164-12-101.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  22. Wenzl P, Li H, Carling J, Zhou M, Raman H, Paul E, Hearnden P, Maier C, Xia L, Caig V: A high-density consensus map of barley linking DArT markers to SSR. RFLP and STS loci and agricultural traits. BMC Genomics. 2006, 7: 206-10.1186/1471-2164-7-206.

    Article  PubMed Central  PubMed  Google Scholar 

  23. Jaccoud D, Peng K, Feinstein D, Kilian A: Diversity arrays: a solid state technology for sequence information independent genotyping. Nucleic Acids Res. 2001, 29: e25-10.1093/nar/29.4.e25.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  24. Gupta PK, Rustgi S, Mir RR: Array-based high-throughput DNA markers for crop improvement. Heredity. 2008, 101 (1): 5-18. 10.1038/hdy.2008.35.

    Article  CAS  PubMed  Google Scholar 

  25. Raman R, Raman H, Kadkol GP, Coombes N, Taylor B, Luckett D: Genome-wide association analyses of loci for shatter resistance in Brassicas. Proceedings of the 17th Australian Research Assembly on Brassicas: 2011. 2011, Wagga Wagga, NSW, Australia, 36-41. (http://www.australianoilseeds.com/__data/assets/pdf_file/0010/8299/S2-P3-Raman.pdf)

    Google Scholar 

  26. Raman H, Raman R, Taylor B, Lindbeck K, Coombes N, Eckermann P, Batley J, Edwards D, Price A, Rehman A: Blackleg resistance in rapeseed: Phenotypic screen, molecular markers, and genome-wide linkage and association mapping. Proceedings of the 17th Australian Research Assembly on Brassicas. 2011, Wagga Wagga, NSW, Australia, 61-64. (http://www.australianoilseeds.com/__data/assets/pdf_file/0013/8311/S3-P1-Raman.pdf)

    Google Scholar 

  27. Raman H, Raman R, Eckermann P, Coombes N, Manoli S, Zou X, Edwards D, Meng J, Prangnell R, Stiller J: Genetic and physical mapping of flowering time loci in oilseed rape (Brassica napus L.). Theor Appl Genet. 2013, 126 (1): 119-132. 10.1007/s00122-012-1966-8.

    Article  CAS  PubMed  Google Scholar 

  28. Yang Y-W, Lai K-N, Tai P-Y, Li W-H: Rates of nucleotide substitution in angiosperm mitochondrial DNA sequences and dates of divergence between Brassica and other angiosperm lineages. J Mol Evol. 1999, 48 (5): 597-604. 10.1007/PL00006502.

    Article  CAS  PubMed  Google Scholar 

  29. Schranz ME, Lysak MA, Mitchell-Olds T: The ABC's of comparative genomics in the Brassicaceae: building blocks of crucifer genomes. Trends Plant Sci. 2006, 11 (11): 535-542. 10.1016/j.tplants.2006.09.002.

    Article  CAS  PubMed  Google Scholar 

  30. Cheung F, Trick M, Drou N, Lim YP, Park J-Y, Kwon S-J, Kim J-A, Scott RJ, Pires C, Paterson AH: Comparative analysis between homoeologous genome segments of Brassica napus and its progenitor species reveals extensive sequence-level divergence. Plant Cell. 2009, 21: 1912-1928. 10.1105/tpc.108.060376.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  31. Wang X, Wang H, Wang J, Sun R, Wu J, Liu S, Bai Y, Mun J-H, Bancroft I, Cheng F: The genome of the mesopolyploid crop species Brassica rapa. Nat Gen. 2011, 43: 1035-1039. 10.1038/ng.919.

    Article  CAS  Google Scholar 

  32. Raman R, Taylor B, Marcroft S, Stiller J, Eckermann P, Coombes N, Rehman A, Lindbeck K, Luckett D, Wratten N: Molecular mapping of qualitative and quantitative loci for resistance to Leptosphaeria maculans; causing blackleg disease in canola (Brassica napus L.). Theor Appl Genet. 2012, 125 (2): 405-418. 10.1007/s00122-012-1842-6.

    Article  CAS  PubMed  Google Scholar 

  33. Raman R, Taylor B, Lindbeck K, Coombes N, Barbulescu D, Salisbury P, Raman H: Molecular mapping and validation of Rlm1 gene for resistance to Leptosphaeria maculans in canola (Brassica napus L.). Crop & Pasture Science. 2012, 63: 1007-1017. 10.1071/CP12255.

    Article  CAS  Google Scholar 

  34. Shi JQ, Li RY, Qiu D, Jiang CC, Long Y, Morgan C, Bancroft I, Zhao JY, Meng JL: Unraveling the complex trait of crop yield with quantitative trait loci mapping in Brassica napus. Genetics. 2009, 182 (3): 851-861. 10.1534/genetics.109.101642.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  35. Feng J, Long Y, Shi L, Shi J, Barker G, Meng J: Characterisation of metabolite quantitative trait loci and metabolic networks that control glucosinolate concentration in the seeds and leaves of Brassica napus. New Phytol. 2012, 193: 96-108. 10.1111/j.1469-8137.2011.03890.x.

    Article  CAS  PubMed  Google Scholar 

  36. Zhao H, Shi L, Duan X, Xu F, Wang Y, Meng J: Mapping and validation of chromsome regions conferring a new boron-efficient locus in Brassica napus. Mol breeding. 2008, 22: 495-506. 10.1007/s11032-008-9193-3.

    Article  CAS  Google Scholar 

  37. Wang J, Long Y, Wu BD, Liu J, Jiang CC, Shi L, Zhao JW, King GJ, Meng JL: The evolution of Brassica napus FLOWERING LOCUS T paralogues in the context of inverted chromosomal duplication blocks. BMC Evol Biol. 2009, 9: 271-10.1186/1471-2148-9-271.

    Article  PubMed Central  PubMed  Google Scholar 

  38. Liu J, Yang J, Li R, Shi L, Zhang C, Long Y, Xu F, Meng J: Analysis of genetic factors that control shoot mineral concentrations in rapeseed (Brassica napus) in different boron environments. Plant Soil. 2009, 320 (1–2): 255-266.

    Article  CAS  Google Scholar 

  39. Long Y, Xia W, Li R, Wang J, Shao M, Feng J, King GJ, Meng J: Epigenetic QTL mapping in Brassica napus. Genetics. 2011, 189: 1093-1102. 10.1534/genetics.111.131615.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  40. Zou X, Suppanz I, Raman H, Hou J, Wang J, Long Y, Jung C, Meng J: Comparative analysis of FLC homologues in Brassicaceae provides insight into their role in the evolution of oilseed rape. PLoS One. 2012, 7 (9): e45751-10.1371/journal.pone.0045751.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  41. Endrigkeit J, Wang X, Cai D, Zhang C, Long Y, Meng J, Jung C: Genetic mapping, cloning, and functional characterization of the BnaX.VTE4 gene encoding a gamma-tocopherol methyltransferase from oilseed rape. Theor Appl Genet. 2009, 119: 567-575. 10.1007/s00122-009-1066-6.

    Article  CAS  PubMed  Google Scholar 

  42. Li F, Ma C, Wang X, Gao C, Zhang J, Wang Y, Cong N, Li X, Wen J, Yi B: Characterization of sucrose transporter alleles and their association with seed yield-related traits in Brassica napus L. BMC Plant Biol. 2011, 11 (1): 168-10.1186/1471-2229-11-168.

    Article  PubMed Central  PubMed  Google Scholar 

  43. Aslam MN: Inheritance of molecular markers for high oleic and low linolenic acids in canola (Brassica napus) and impact of genotype and environment on these traits. PhD Thesis. School of Plant Biology. 2007, The University of Western Australia

    Google Scholar 

  44. Zhao J, Meng J: Detection of loci controlling seed glucosinolate content and their association with Sclerotinia resistance in Brassica napus. Plant Breed. 2003, 122: 19-23. 10.1046/j.1439-0523.2003.00784.x.

    Article  CAS  Google Scholar 

  45. Wang N, Qian W, Suppanz I, Wei L, Mao B, Long Y, Meng J, Müller AE, Jung C: Flowering time variation in oilseed rape (Brassica napus L.) is associated with allelic variation in the FRIGIDA homologue BnaA.FRI.a. J Exp Bot. 2011, 62 (15): 5641-5658. 10.1093/jxb/err249.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  46. Fritsche S, Wang X, Li J, Stich B, Kopisch-Obuch FJ, Endrigkeit J, Leckband G, Dreyer F, Friedt W, Meng J: A candidate gene-based association study of tocopherol content and composition in rapeseed (Brassica napus). Front. Plant Science. 2012, 3: 129-

    Article  CAS  Google Scholar 

  47. Li YY, Shen J, Wang T, Chen Q, Zhang X, Fu T, Meng J, Tu J, Ma C: QTL analysis of yield-related traits and their association with functional markers in Brassica napus L. Aust J Agric Res. 2007, 58 (8): 759-766. 10.1071/AR06350.

    Article  CAS  Google Scholar 

  48. Nelson MN, Nixon J, Lydiate DJ: Genome-wide analysis of the frequency and distribution of crossovers at male and female meiosis in Sinapis alba L. (white mustard). Theor Appl Genet. 2005, 111: 31-43. 10.1007/s00122-005-1961-4.

    Article  CAS  PubMed  Google Scholar 

  49. Ramchiary N, Nguyen VD, Li X, Hong CP, Dhandapani V, Choi SR, Yu G, Piao ZY, Lim YP: Genic microsatellite markers in Brassica rapa: Development, characterization, mapping, and their utility in other cultivated and wild Brassica relatives. DNA Res. 2011, 18 (5): 305-320. 10.1093/dnares/dsr017.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  50. Milczarski P, Bolibok-Bragoszewska H, Myskow B, Stojalowski S, Heller-Uszynska K, Goralska M, Bragoszewski P, Uszynski G, Kilian A, Rakoczy-Trojanowska M: A high density consensus map of rye (Secale cereale L.) based on DArT markers. PLoS One. 2011, 6 (12): e28495-10.1371/journal.pone.0028495.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  51. Mace E, Xia L, Jordan D, Halloran K, Parh D, Huttner E, Wenzl P, Kilian A: DArT markers: diversity analyses and mapping in Sorghum bicolor. BMC Genomics. 2008, 9: 26-10.1186/1471-2164-9-26.

    Article  PubMed Central  PubMed  Google Scholar 

  52. Arumuganathan K, Earle ED: Nuclear DNA content of some important plant species. Plant Mol Biol Reporter. 1991, 9: 208-218. 10.1007/BF02672069.

    Article  CAS  Google Scholar 

  53. Bohuon EJR, Keith DJ, Parkin IAP, Sharpe A, Lydiate DJ: Alignment of the conserved C genome of Brassica oleracea and Brassica napus. Theor Appl Genet. 1996, 93: 833-839. 10.1007/BF00224083.

    Article  CAS  PubMed  Google Scholar 

  54. Nicolas SD, Mignon GL, Eber F, Coriton O, Monod H, Clouet V, Huteau V, Lostanlen A, Delourme R, Chalhoub B: Homeologous recombination plays a major role in chromosome rearrangements that occur during meiosis of Brassica napus haploids. Genetics. 2007, 175 (2): 487-503. 10.1534/genetics.106.062968.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  55. Nicolas SD, Monod H, Eber F, Che`vre A-M, Jenczewski E: Non-random distribution of extensive chromosome rearrangements in Brassica napus depends on genome organization. Plant J. 2012, 70 (4): 691-703. 10.1111/j.1365-313X.2012.04914.x.

    Article  CAS  PubMed  Google Scholar 

  56. Parkin IA, Sharpe AG, Lydiate DJ: Patterns of genome duplication within the Brassica napus genome. Genome. 2003, 46 (2): 291-303. 10.1139/g03-006.

    Article  CAS  PubMed  Google Scholar 

  57. Sharpe A, Parkin IA, Keith D, Lydiate DJ: Frequent reciprocal translocations in the amphidiploid genome of oilseed rape (Brassica napus). Genome. 1995, 38: 1112-1121. 10.1139/g95-148.

    Article  CAS  PubMed  Google Scholar 

  58. Jiang C, Ramchiary N, Ma Y, Jin M, Feng J, Li R, Wang H, Long Y, Choi S, Zhang C: Structural and functional comparative mapping between the Brassica A genomes in allotetraploid Brassica napus and diploid Brassica rapa. Theor Appl Genet. 2011, 123: 927-941. 10.1007/s00122-011-1637-1.

    Article  CAS  PubMed  Google Scholar 

  59. Mayerhofer R, Wilde K, Mayerhofer M, Lydiate D, Bansal VK, Good AG, Parkin IAP: Complexities of chromosome landing in a highly duplicated genome: Toward map-based cloning of a gene controlling blackleg resistance in Brassica napus. Genetics. 2005, 171 (4): 1977-1988. 10.1534/genetics.105.049098.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  60. Hauge BM, Hanley SM, Cartinhour S, Cherry JM, Goodman HM, Koornneef M, Stam P, Chang C, Kempin S, Medrano L: An integrated genetic/RFLP map of the Arabidopsis thaliana genome. Plant J. 1993, 3 (5): 745-754. 10.1111/j.1365-313X.1993.00745.x.

    Article  CAS  Google Scholar 

  61. Akhunov ED, Goodyear AW, Geng S, Qi L, Echalier B, Gill BS, Miftahudin , Gustafson JP, Lazo G, Chao S: The organization and rate of evolution of wheat genomes are correlated with recombination rates along chromosome arms. Genome Res. 2003, 13: 753-763. 10.1101/gr.808603.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  62. Zhang J, Liu H, Sun J, Li B, Zhu Q, Chen S, Hongxia Zhang H: Arabidopsis fatty scid desaturase FAD2 is required for salt tolerance during seed germination and early seedling growth. PLoS One. 2012, 7 (1): e30355-10.1371/journal.pone.0030355.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  63. Rouxel T, Willner E, Coudard L, Balesdent M-H: Screening and identification of resistance to Leptosphaeria maculans (stem canker) in Brassica napus accessions. Euphytica. 2003, 133 (2): 219-231. 10.1023/A:1025597622490.

    Article  CAS  Google Scholar 

  64. Raman R, Allen H, Diffey S, Raman H, Martin P, McKelvie K: Localisation of quantitative trait loci for quality attributes in a doubled haploid population of wheat (Triticum aestivum L). Genome. 2009, 52 (8): 701-715. 10.1139/G09-045.

    Article  CAS  PubMed  Google Scholar 

  65. Girin T, Stephenson P, Goldsack CMP, Kempin SA, Perez A, Pires N, Sparrow PA, Wood TA, Yanofsky MF, Østergaard L: Brassicaceae INDEHISCENT genes specify valve margin cell fate and repress replum formation. Plant J. 2010, 63 (2): 329-338. 10.1111/j.1365-313X.2010.04244.x.

    Article  CAS  PubMed  Google Scholar 

  66. Davis GL, Dibner DM, Battey FJ: Preparation of DNA from Eukaryotic cells: General method. Basic Methods in Molecular Biology. Edited by: Davis LG. 1989, New York: Elsevier, 44-46.

    Google Scholar 

  67. Wenzl P, Carling J, Kudrna D, Jaccoud D, Huttner E, Kleinhofs A, Kilian A: Diversity arrays technology (DArT) for whole-genome profiling of barley. Proc Nat Acad Sci, USA. 2004, 101: 9915-9920. 10.1073/pnas.0401076101.

    Article  CAS  Google Scholar 

  68. Suwabe K, Morgan C, Bancroft I: Integration of Brassica A genome genetic linkage map between Brassica napus and B. rapa. Genome. 2008, 41: 169-176.

    Google Scholar 

  69. Suwabe K, Tsukazaki H, Iketani H, Hatakeyama K, Kondo M, Fujimura M, Nunome T, Fukuoka H, Hirai M, Matsumoto S: Simple sequence repeat-based comparative genomics between Brassica rapa and Arabidopsis thaliana: The genetic origin of clubroot resistance. Genetics. 2006, 173 (1): 309-319. 10.1534/genetics.104.038968.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  70. Schuelke M: An economic method for the fluorescent labeling of PCR fragments. Nat Biotechnol. 2000, 18: 233-234. 10.1038/72708.

    Article  CAS  PubMed  Google Scholar 

  71. Hudson CJ, Freeman JS, Kullan AR, Petroli CD, Sansaloni CP, Kilian A, Detering F, Grattapaglia D, Potts BM, Myburg AA: A reference linkage map for Eucalyptus. BMC Genomics. 2012, 13: 240-10.1186/1471-2164-13-240.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  72. Cheema J, Dicks J: Computational approaches and software tools for genetic linkage map estimation in plants. Brief Bioinform. 2009, 10: 595-608. 10.1093/bib/bbp045.

    Article  CAS  PubMed  Google Scholar 

  73. Wu Y, Bhat PR, Close TJ, Lonardi S: Efficient and accurate construction of genetic linkage maps from the minimum spanning tree of a graph. PLoS Genet. 2008, 4 (10): e1000212-10.1371/journal.pgen.1000212.

    Article  PubMed Central  PubMed  Google Scholar 

  74. Kosambi DD: The estimation of map distances from recombination values. Ann Eugen. 1944, 12: 172-175.

    Article  Google Scholar 

  75. Van Ooijen JW: JoinMap 4.0: Software for the calculation of genetic linkage maps in experimental populations. 2006, Wageningen, Netherlands: Kyazma B.V

    Google Scholar 

  76. Voorrips RE: MapChart: Software for the graphical presentation of linkage maps and QTLs. J Hered. 2002, 93 (1): 77-78. 10.1093/jhered/93.1.77.

    Article  CAS  PubMed  Google Scholar 

  77. Altschul SF, Madden TL, Schäffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25 (17): 3389-3402. 10.1093/nar/25.17.3389.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  78. Thompson JD, Higgins DG, Gibson TJ: CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994, 22 (22): 4673-4680. 10.1093/nar/22.22.4673.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  79. Drummond AJ, Ashton B, Cheung M, Heled J, Kearse M, Moir R, Stones-Havas S, Thierer T, Wilson A: Geneious Pro v 5.4.6. 2005–2011, Available at http://www.geneious.com

    Google Scholar 

Download references

Acknowledgments

Authors acknowledge the International Rapeseed Consortium, NSW Agricultural Genomics Centre, Australia and the Grains Research and Development Corporation, NSWDPI, and UQ for supporting the Australian Canola Molecular Marker Program (DAN 117). M.N.A. was supported by Australian Research Council project LP0210571, with industry partners Council of Grain Grower Organizations Ltd (Australia) and Norddeutsche Pflanzenzucht Hans-Georg Lembke KG (Germany). The authors are thankful to the Victorian Department of Primary Industries for providing the doubled haploid population from Maxol/Westar, and Dr Neil Coombes (NSWDPI) for statistical analysis, Dr Sahana Manoli and Mr. Pradeep Ruparao (UQ) for in silico mapping of DArT clones.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Harsh Raman.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

HR and AK coordinated and designed the study. HR, RR, MNN, and FD, wrote the manuscript. RR, YL and MNA isolated DNA samples and generated the PCR marker data of mapping populations. AK, FD and HR analyzed the DArT data and constructed genetic and consensus maps. HR, JZ, YL and JM compiled published trait-marker associations. DE, IP, JB and AS performed in silico physical mapping with A and C genome scaffolds. HR investigated DArT sequence analysis. MNA, WAC and MNN analysed fatty acid profiles and conducted QTL mapping in the LM population. NL and DL provided AT mapping population. All the authors have commented, read and approved the final manuscript.

Electronic supplementary material

12864_2012_4970_MOESM1_ESM.xls

Additional file 1: Molecular markers and their map positions on the genetic linkage map of the doubled haploid populations and consensus map.(XLS 846 KB)

12864_2012_4970_MOESM2_ESM.rtf

Additional file 2: Comparative linkage map of different chromosomes constructed from the individual six doubled haploid populations derived from Ag-Castle/Topas (AT), BLN2762/Surpass400 (BS), Maxol*1/Westar-10 (MW), Lynx-037DH/Monty-028DH (LM), Skipton/Ag-Spectrum (SAS) and Tapidor/Ningyou7 (TN), respectively (for enlarged view × 1.5).(RTF 14 MB)

12864_2012_4970_MOESM3_ESM.rtf

Additional file 3: Percent sequence identities of DArT clones that map in a cluster on different chromosomes of the consensus map.(RTF 172 KB)

12864_2012_4970_MOESM4_ESM.rtf

Additional file 4: Consensus map of B. napus constructed from the six doubled haploid populations derived from Ag-Castle/Topas, BLN2762/Surpass400, Lynx-037DH/Monty-028DH, Maxol*1/Westar-10, Skipton/Ag-Spectrum, and Tapidor/Ningyou7, respectively (for enlarged view × 5).(RTF 14 MB)

12864_2012_4970_MOESM5_ESM.rtf

Additional file 5: Comparisons between marker density and map length (in cM) in the consensus map developed from six component individual maps.(RTF 109 KB)

Additional file 6: Shared markers across different mapping populations.(RTF 427 KB)

Additional file 7: Sequence identities of DArT clones with the genome scaffold of Brassica rapa.(XLS 3 MB)

Additional file 8: Sequence identities of DArT clones with the genome scaffold of Brassica oleracea.(XLS 7 MB)

12864_2012_4970_MOESM9_ESM.xls

Additional file 9: List of DArT clones that did not show any 'hit' with the genome scaffolds of Brassica rapa and Brassica oleracea.(XLS 44 KB)

12864_2012_4970_MOESM10_ESM.rtf

Additional file 10: Relationship between genetic map distance and physical map distance for all the 19 Brassica A and C genome chromosomes. Genetic distances are given in cM, derived from the B. napus consensus map. Map order of linkage groups was orientated according to genome scaffolds. DArT sequences which were genetically mapped on to linkage groups were only aligned against the corresponding scaffolds. DArT sequences those showed multiple sequence alignments within the linkage group (chromosome) were also included. (RTF 2 MB)

12864_2012_4970_MOESM11_ESM.xls

Additional file 11: List of inter-and intra homoalleles of marker loci mapped on the consensus map developed from six mapping populations.(XLS 78 KB)

12864_2012_4970_MOESM12_ESM.xls

Additional file 12: Chromosomal position of molecular markers associated with various traits segregating in the mapping populations used for construction of consensus map.(XLS 189 KB)

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Raman, H., Raman, R., Kilian, A. et al. A consensus map of rapeseed (Brassica napusL.) based on diversity array technology markers: applications in genetic dissection of qualitative and quantitative traits. BMC Genomics 14, 277 (2013). https://0-doi-org.brum.beds.ac.uk/10.1186/1471-2164-14-277

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/1471-2164-14-277

Keywords