- Research article
- Open Access
High-density SNP-based genetic map development and linkage disequilibrium assessment in Brassica napusL
BMC Genomics volume 14, Article number: 120 (2013)
High density genetic maps built with SNP markers that are polymorphic in various genetic backgrounds are very useful for studying the genetics of agronomical traits as well as genome organization and evolution. Simultaneous dense SNP genotyping of segregating populations and variety collections was applied to oilseed rape (Brassica napus L.) to obtain a high density genetic map for this species and to study the linkage disequilibrium pattern.
We developed an integrated genetic map for oilseed rape by high throughput SNP genotyping of four segregating doubled haploid populations. A very high level of collinearity was observed between the four individual maps and a large number of markers (>59%) was common to more than two maps. The precise integrated map comprises 5764 SNP and 1603 PCR markers. With a total genetic length of 2250 cM, the integrated map contains a density of 3.27 markers (2.56 SNP) per cM. Genotyping of these mapped SNP markers in oilseed rape collections allowed polymorphism level and linkage disequilibrium (LD) to be studied across the different collections (winter vs spring, different seed quality types) and along the linkage groups. Overall, polymorphism level was higher and LD decayed faster in spring than in “00” winter oilseed rape types but this was shown to vary greatly along the linkage groups.
Our study provides a valuable resource for further genetic studies using linkage or association mapping, for marker assisted breeding and for Brassica napus sequence assembly and genome organization analyses.
Genetic linkage maps are highly valuable tools for comparative genome analyses and the identification of genomic regions carrying major genes and quantitative trait loci (QTL) controlling agronomical traits. They are a prerequisite for further map-based cloning or marker-assisted breeding programs. In recent years, the establishment of genetic maps have benefited from the development of new types of molecular markers which take advantage of automated sequencing and genotyping technologies. While the first marker-based genetic maps were built with restriction fragment length polymorphisms (RFLPs), random amplified polymorphic DNAs (RAPDs) and amplified fragment length polymorphisms (AFLPs), dense genetic maps now include simple sequence repeats (SSRs) and more recently single nucleotide polymorphisms (SNPs). Dense genetic maps based on sequence-derived markers allow finer comparative genome analyses to be performed based on comparisons with sequenced related genomes and to accelerate the process of map-based cloning of major genes and QTL. They are also very useful tools to assist sequence assembly in whole de novo genome sequencing projects [1–3]. Moreover, by integrating genetic map data with genotyping data generated from collections of accessions/varieties linkage disequilibrium (LD) pattern along the genome of a given species can be investigated, which is a prerequisite for precise genome wide association studies (GWAS). GWAS performed with a large number of SNPs have been reported in a number of crop species such as maize [4, 5], Arabidopsis , barley [7–11], and rice [12–15] The success of GWAS to locate genes responsible for complex traits depends on the extent of LD, the number, the distribution and the diversity of markers and the underlying structure in the studied collections. Since the diversity of markers and the extent of LD may vary depending on the history of the collections [7, 15], they should be investigated prior to GWAS design.
Oilseed rape (Brassica napus) is a prominent oilseed crop in most world continents including America, Europe, Australia and Asia and is cultivated for food (oil) and feed (meal) as well as for non-food uses such as biofuels or lubricants. It is the second world oil crop after soybean (http://faostat3.fao.org/home/index.html; May 2011) with a world production of more than 60 million tonnes per year. B. napus is an amphidiploid species (AC genome, n = 19) that arose from hybridization between B. rapa (A genome, n = 10) and B. oleracea (C genome, n = 9) diploid species  within the past 10,000 years . B. napus includes spring and winter oilseed rape, rutabaga or swede, and some fodder crops. It likely originated from a few interspecific hybridization events  and has only a short domestication history of about 400–500 years [17, 19]. For these reasons, the genetic diversity within B. napus germplasm is rather low compared to that of its two progenitor species B. rapa and B. oleracea. Moreover, two bottlenecks have occurred during breeding of modern oilseed rape varieties through the selection for low erucic acid content in the oil and low glucosinolate content in the seeds, which reduced the genetic diversity in modern varieties .
Over the last 20 years, many genetic B. napus maps have been built, which have been progressively integrating various types of markers [21–27]. These maps have been used for genetic studies of various agronomical traits including development traits [28, 29], seed quality [30–34], yield components [35–38] and disease resistance [39–44] as well as for genetic study of chromosome pairing . The establishment of genetic maps of diploid and amphidiploid Brassica species, and their comparison and alignment to Arabidopsis genome sequence provided insights into Brassica genome organization and evolution after the different rounds of polyploidization and diploidization occurring in these species history. Extensive collinearity observed between A. thaliana and B. napus led to the description of a genomic block system determined by Parkin et al. , who demonstrated that the structure of the Brassica A and C genomes could be described with approximately 21 conserved blocks. A framework built of 24 genomic blocks (A-X) within the ancestral karyotype was then proposed that represents an extension of the above mentioned study . This conserved block structure was then further investigated in related species such as B. juncea or B. oleracea and led to a block arrangement comparison in the A, B and C genomes . It was also recently confirmed in B. napus using dense genetic maps with SSR , and SNP  markers.
The availability of high numbers of markers now makes it possible to investigate more precisely genome wide diversity and the extent of LD in oilseed rape (OSR). To date, the published studies relied on either a low number of lines or low number of markers. Ecke et al.  used 845 AFLP markers to examine the extent of LD in a set of 85 winter oilseed rape lines. Bus et al.  investigated patterns of genetic diversity and the extent of LD in 509 inbred lines corresponding to different germplasms of oilseed rape with 89 SSR markers. Xiao et al.  assessed the genetic diversity and the extent of LD in a panel of 192 inbred lines from all over the world but with a great proportion originating from China using 451 SSR markers. Harper et al.  carried out associative transcriptomics on 53 B. napus lines using >60 K SNPs and confirmed the low overall level of LD in B. napus.
In this context, our objectives were to (i) obtain a dense integrated SNP genetic map of B. napus built from four segregating populations (ii) investigate the polymorphism of these SNPs within and among different germplasm types of a large B. napus collection and (iii) assess the extent and pattern of LD between densely mapped SNP markers.
A total of 7322 SNPs was selected for Infinium genotyping according to the criteria described in methods section. Of these, 5986 were retained to build the four individual maps as they exhibited clear segregation patterns after Genome Studio analysis. The numbers of SNP markers mapped were 2664, 2763, 3385 and 2301 for the DYDH (‘Darmor-bzh’ × ‘Yudal’ doubled haploid), TNDH (‘Tapidor’ × ‘Ningyou7’ doubled haploid), AADH (‘Aviso’ × ‘Aburamasari’ doubled haploid) and AMDH (‘Aviso’ × ‘Montego’ doubled haploid) populations, respectively (Table 1). In addition, 833 and 831 PCR markers (SSR and sequence-derived markers) were mapped in the DYDH and TNDH populations. Individual maps covered 1947 cM for the TNDH and AMDH populations, 2049 cM for the DYDH and 3495 cM for the AADH ones and the numbers of markers per cM were 1.85, 1.18, 1.71 and 0.97, respectively (Table 1 and Additional file 1: Table S1). The percentage of mapped SNPs that showed segregation distortion was estimated at 46.2%, 28.3%, 16.1% and 16.3% in the DYDH, TNDH, AADH and AMDH populations, respectively. Most of the linkage groups in the DYDH and TNDH maps showed segregation distortion except A1, A4, A8, A10, C1, C4, C8 on DYDH and A1, C1, C2, C4, C5, C9 on the TNDH maps (Additional file 2: Table S2). On the AMDH map, only A2, A3, A9, C1 and C9 showed segregation distortion. Many regions were distorted in more than one map such as regions on A2, A3, A4, A9, A10, C3, C6, C7, C8 and C9.
The integrated map comprised 7367 markers (5764 SNP and 1603 PCR markers) and covered 2250 cM, which corresponds to a density of 3.27 markers (2.56 SNPs) per cM (Table 2). Twice as many SNP markers were assigned a position on the A genome compared to the C genome: 3942 (68.4%) and 1822 (31.6%) were mapped on the A and C genomes, respectively, which correspond to 4.6 and 2.17 SNPs per cM. This difference between genomes A and C was less pronounced for the AMDH population. Of these 5764 SNP markers, 2350 (41%), 2150 (37%), 1112 (19%) and 152 (3%) were mapped in only one or were common to two, three or four populations, respectively. The AMDH population had the smallest number of markers (28%) in common with the other populations whereas 35 to 40% of the markers were common between the AADH, DYDH and TNDH populations. These percentages were similar if we consider the SNPs that were mapped to the A and C genomes, respectively.
The recombination rate on the AADH map was higher than on the others (Additional file 3: Figures S1 and Additional file 4: Figure S2) as expected due to the mating scheme used to produce the population. This increase was relatively homogeneous across all the linkage groups, except for the bottom of A2 and the top of C2. An increase in recombination rate was also observed at the top of A3 on the TNDH map. Overall, very good collinearity was observed between all the maps with the exception of three inversions that were observed at the bottom of A2 on the TNDH map, at the top of C8 on the AADH map and at the bottom of C8 on the AMDH map. The AMDH population was the least polymorphic with the lowest number of markers and an overall lower marker density: many regions were not polymorphic at all or had a very scarce number of markers. Non-polymorphic regions were observed on almost all the linkage groups and in particular very low coverage was obtained for linkage groups A8, A10, C5 and C8.
We anchored the genetic maps onto the Arabidopsis genome using homology searches with the SNP context sequences and the previous anchorage method described in Wang et al.  for SSR and other PCR markers. Out of the 7367 mapped markers, 5725 gave hits with Arabidopsis genes. From these hits, 119 collinearity blocks were identified and represented in relation to the 24 blocks defined by Schranz et al.  (Figure 1). Since the collinearity between the different individual maps was very good, the conservation of the block organization between these individual maps was also very good with some additional blocks on some maps due to the different number of markers mapped in these regions (Additional file 3: Figure S1). However, few polymorphic markers were identified in some blocks especially on the AMDH map, which, as mentioned above, showed the lowest polymorphism rate. Indeed, on this map some blocks were totally or partially monomorphic (e.g. F on A1; E on A2 and C2; W, J and I on A3; V on A6; E on A7; B and A on A8; A on A9; W on A10; A and N on C8; R on C9). The A block was also partially monomorphic or missing on A8 of the DYDH map as well as the R block on A2 and C2 of the TNDH map.
In addition to these 7367 markers included in the integrated map, 222 markers were assigned to different linkage groups across the individual maps (Additional file 5: Table S3). Further investigation of these cases showed that these markers were included in the same synteny blocks and mapped to duplicated regions present on different linkage groups. Thus they appear to be homoeologous or duplicated loci.
Validation of a 1536 SNP public set in a GoldenGate assay
A set of 1536 SNPs was selected for even distribution over the integrated linkage map and a good level of polymorphism in the OSR collection. The characteristics of these 1536 SNPs including their anchorage onto Arabidopsis and B. rapa sequences are described in Additional file 6: Table S4. Of these 1536 SNPs, 1104 (72%) were successfully genotyped on the Illumina BeadXpress platform. Twenty-nine could not be localized on the integrated map as they were assigned to different linkage groups on the individual maps. The distribution of the remaining 1507 SNPs, that were located on the integrated linkage map (Table 3), showed that the mean distances between two SNP markers were 1.32 and 1.79 cM for the original set and the GoldenGate validated set, respectively, which on average corresponded to a 35% increase. Consequently, the number of gaps over 10 cM increased from 22 to 32, out of 1488 and 1069 intervals, respectively in these two sets. Out of the 1536 SNPs, 1412 and 1527 showed a significant hit with Arabidopsis genes or B. rapa sequences, respectively. A total of 1461 SNPs had a significant hit on the B. rapa pseudo-chromosomes: 883 and 578 were located on the A and C linkage groups, respectively. Ninety-five percent of the SNPs located on A linkage groups showed a hit on the corresponding pseudo-chromosomes. Eighty-six percent of the SNPs located on the C linkage groups showed a hit on the expected pseudo-chromosomes according to the known collinearity between the A and C genomes.
Polymorphism in the B. napuscollection
A total of 5685 SNPs were validated and scored on the B. napus collection. On average, only 1.1% of the data (0–10.3) was missing and 1.4% (0–21.9) was heterozygous (or of a mix of homozygous and heterozygous) per variety. Only 13 varieties showed more than 10% of heterozygosity. A total of 4881 SNPs with a minor allele frequency (MAF) greater than 5% on the whole collection was retained for further polymorphism study and analysis of molecular variance (AMOVA). An AMOVA was carried out to assess genetic differentiation between fodder, spring oilseed rape and winter oilseed rape varieties, and between the three seed quality subgroups (“++”, “0+” and “00”) within the spring and winter oilseed rape types. The level of within subgroup variation was 72.5 and 68.7% whether we considered fodder varieties or not, with FST indices of 0.275 and 0.312 and differentiation indices between types of 0.181 and 0.223, respectively. The FST indices were 0.089 and 0.170 within winter and spring oilseed rape types, respectively.
Of the 4881 SNPs above, 4363 were localised on the integrated map and were further used for principal component analyses (PCA). Three PCAs were performed using either i) the whole set of mapped SNPs; ii) the 2854 SNPs that were mapped on the A genome or iii) the 1509 SNPs that were mapped on the C genome. The first two axes accounted for 19.2%, 30.6% and 20.1% of the variation, respectively in these three PCAs (Figure 2). Out of the 1536 SNPs selected subset, 1507 SNPs were located on the integrated map. In the PCA performed with these 1507 SNPs or with the 908 and 599 SNPs mapped on the A and C genomes, respectively, the first two axes accounted for 18.5%, 18.3% and 20.0% of the variation. In each case, the first axis differentiated the spring oilseed rape from the winter oilseed rape accessions and the fodder rape varieties were between the two groups. The second axis mainly discriminated European and Canadian (“0+” and “00”) spring oilseed rape from Asian spring “++” oilseed rape. No clear differentiation between “++”, “+0” and “00” winter oilseed rape subgroups was observed.
The 4363 SNPs that were localised on the integrated map were further used for investigating polymorphism variation across the linkage groups. Table 4 and Additional file 7: Figure S3 show the mean polymorphism information content (PIC) values for each linkage groups in the whole collection, the fodder, spring and winter oilseed rape types, as well as in the “++”, “0+” and “00” winter oilseed rape subgroups. The polymorphism level was significantly higher in spring oilseed rape than in fodder and winter oilseed rape and was significantly lower in winter oilseed rape than in fodder rape. The exception was the C1 linkage group where the polymorphism was much lower in the spring oilseed rape types. On average, the polymorphism level was slightly lower for the A than for the C linkage groups except for the spring oilseed rape types and there was a great variability between linkage groups. A1, A2, A9 and A10 were the least polymorphic linkage groups for the A genome, especially in winter oilseed rape. For the C genome, C2, C8, C9 and C2, C9 were the least polymorphic in winter oilseed rape and fodder types, respectively. For winter oilseed rape, the polymorphism level was quite similar between the “++”, “0+” and “00” subgroups for the A linkage groups whereas there was a greater variability for the C linkage groups with some linkage groups such as C2, C4, C5, C9 showing lower PIC values in the “0+” and “00” subgroups and other linkage groups such as C1, C3, C7 showing higher PIC values in the “00” subgroup. The mean PIC values were also estimated for the whole collection with the mapped SNPs from the 1536 SNP subset. These were significantly higher than the ones estimated with the full SNP set (Table 4), due to the criteria used for their selection.
Additional file 7: Figure S3 shows that the PIC values varied differently along the different linkage groups. When we compared fodder rape, spring oilseed rape (SOSR) and winter oilseed rape (WOSR), there were large regions on most A linkage groups except A8 where the level of polymorphism was higher in the spring oilseed rape types. The polymorphism level in winter oilseed rape was either lower or similar to the one observed in fodder rape, depending on the A regions. On the C linkage groups, the level of polymorphism was lower in some regions in winter oilseed rape e.g. on C2, C4, C6 and C8 linkage groups but elsewhere the variation in PIC values between the different types was more erratic. The C1 and to a lesser extent the C5 linkage groups showed a contrasted situation with a lower polymorphism in spring oilseed rape along most of their length. When we compared the three winter oilseed rape subgroups (“++”, “0+” and “00”), variation in PIC values along the different linkage groups was very contrasted since it decreased or increased in the “0+” and/or “00” subgroups depending on regions. Nevertheless, regions showing a decrease in polymorphism were more numerous.
The 4329 SNPs that were localised on the integrated map and had a MAF above 5% were then used for mean LD estimation and LD pattern study along the different linkage groups depending on the oilseed rape types. The whole collection or spring, winter and “00” winter oilseed rape were considered so that a sufficient number of varieties was included to estimate the LD.
The mean pairwise r2 was estimated at 0.037, 0.057, 0.017 and 0.021 in the whole, the spring, the winter and the “00” winter oilseed rape collections, respectively. This corresponded to 2.9%, 6.6%, 0.4% and 0.7% of the SNP pairs showing a r2 value higher than 0.2 (Table 5). In the whole collection and the spring types, a high percentage of the pairs with a significant LD were between SNP markers located on different chromosomes (85-90%) whereas in the winter types, most were intra-chromosomic pairs. The r2 value was greater than 0.5 for 0.12% to 0.25% of the pairs and greater than 0.8 for less than 0.1% of the pairs. Very few pairs with r2 > 0.8 were observed between SNP markers located on different chromosomes.
LD decay was estimated globally and for each linkage groups from the four collections. The non-linear regression of the LD measure (r2) relative to genetic map distance and the genetic distance at which the estimated r2 fell below 0.2, as well as the effective size, were estimated (Table 6). The trend lines of these non-linear regressions are shown in Additional file 8: Figure S4. The genetic distance at which the estimated r2 fell below 0.2 was 0.6-0.7 cM for the whole collection, the spring and winter types whether it was estimated on the whole genome or the A and C genomes. However on the C genome for the spring types, this distance was estimated at 1.2 cM. For the “00” winter types, it was estimated at 1.2 cM on the whole genome as well as on the genomes A or C. This value varied depending on the linkage group and the collection, ranging from 0.2 to 3.4 cM. On all the linkage groups, the extent of LD was overall higher for the “00” winter types than for the winter types that included the “++”, “0+” and “00” varieties. Different linkage groups showed a higher extent of LD depending on whether it was estimated from the spring or the winter types e.g. A9, C2 and C9 for the spring types and A2, A6, A8 and C6 for the winter types. The effective size varied accordingly. The very different LD patterns between spring and winter types was also evident from the LD plots obtained for each linkage group. This was true for the linkage groups cited above such as for C9 (Figure 3) but some other differences could be seen on most of the groups (Additional file 9: Figure S5).
In this study, we could built a high density SNP integrated B. napus map and depict polymorphism level and LD decay over the linkage groups across different B. napus collections by integrating genotyping data of a large set of SNPs in both segregating populations and diverse collections.
As reviewed by Kaur et al. , SNP discovery is challenging in allopolyploid species such as B. napus. SNPs may arise both between allelic (homologous) sequences within subgenomes and between homoeologous sequences among subgenomes but also from polymorphisms between paralogous duplicated sequences. SNP discovery has been based on B. napus ESTs sequence analysis [55–57] or on second generation high throughput sequencing [50, 58–60]. SNP genotyping using Illumina GoldenGate assays was shown to be possible in B. napus species after careful selection of the SNP . Here, we report the first high throughput genotyping study in oilseed rape using Illumina Infinium and GoldenGate technology. From the 1536 SNPs tested on the two platforms, 1104 were validated on both, suggesting that it should be possible to use one platform or the other depending on the required number of SNPs.
In this study, we generated an integrated map with 7367 markers including 5764 SNPs, which corresponds to 3.3 markers every cM. This large number of mapped markers was obtained by integrating four mapping populations. Marker density actually increased from one marker (one SNP) every 0.54 (0.7-0.8) cM on the individual maps to one marker (SNP) every 0.3 (0.39) cM on the integrated map. This density was comparable to that obtained in Wang et al.  (one marker every 0.34 cM) but the number of markers in common between populations was much higher in our study (60% of the markers were common to at least two maps in our study compared to only 20% in Wang et al. ). The size of the population was also larger in our study, which led to a more accurate ordering of the markers on the integrated map. Bancroft et al.  built a B. napus map with 23027 SNPs by transcriptome sequencing of 37 TNDH lines. These SNP were distributed in 527 recombination bins (one bin corresponded to SNPs having the same scoring genotyping data). Our integrated maps exhibited 3177 bin loci (one bin corresponded to SNPs mapped within 0.1 cM), thus the mapping resolution increased considerably with the increase in population size. The collinearity blocks we identified were compared to those reported in Panjabi et al. , Parkin , Wang et al.  and Bancroft et al. . Compared to Bancroft et al. , which is at present the most complete study, five and 10 small blocks were missing on the A and C genomes, respectively. Nevertheless, three new blocks were identified in our study: S on A8, V-W on C2 and B on C7. In Bancroft’s study, two markers corresponding to these blocks were identified but not declared as a block. Moreover, the V-W and B blocks were identified in the homoeologous regions on A2 and A7 linkage groups, which supports their occurrence in these genomic regions.
A higher number of markers was mapped on the A genome than on the C genome, as previously reported by Bancroft et al. . As in their study, here, this difference between the A and C genomes was more pronounced on the three crosses involving Asiatic parental lines. The hypothesis is that Asiatic cultivars are partly derived from crosses involving its progenitor species B. rapa[61, 62], which increased the genetic diversity of the A genome. This introgression of B. rapa genetic information was recently shown for ‘Nignyou7’ (one of the parent of the TNDH map) by Bancroft et al. . Indeed, there was less difference between the number of markers mapped on the A and C genomes in the AMDH population which is derived from a cross between two French winter OSR varieties. No such difference was observed in Wang et al.  but the markers mapped in their study were mainly contributed by two crosses involving resynthetised B. napus, which is enriched in polymorphisms in both genomes. We probably succeeded in capturing this high level of polymorphism on the A genome because the original sequencing to identify the SNP was performed on material that included Asiatic varieties.
We observed a high level of segregation distortion especially in the DYDH and TNDY populations. Such segregation distortions were reported in many B. napus maps (e.g. ). Segregation distortion and clustering of the skewed loci are common features of microspore-derived DH populations in various species ( for a review), including oilseed rape and may be related to differential responsiveness to microspore culture between the two parental lines, which leads to skewed loci in regions involved in the microspore culture responsiveness.
Very good collinearity was observed between all the individual maps, which made it easy to integrate the four individual maps accurately. Only three inversions were identified on the A2 and C8 linkage groups. These inversions could be due to mapping inaccuracies and need to be confirmed. Due to careful selection at the beginning of the study to target homologous SNPs, very few SNPs were not assigned a position on the same linkage groups among the different maps. Those SNPs that did map to different linkage groups were located in duplicated regions within or between the A and C genomes where there is a high level sequence similarity. The map derived from the AADH population was 75% bigger than those derived from the other crosses, which was expected from the way the DH population was obtained. DH lines were produced after intermating F2 plants, which increased the number of recombination events. This type of highly recombinant population is of interest for obtaining better mapping resolution [64, 65]. The map derived from the AMDH population was the least dense with some regions missing, due to the lack of polymorphism between the two parental lines. In many cases, the monomorphic regions corresponded to quite complete B. napus/Arabidopsis collinearity blocks. The lower marker density can be related to the lower level of polymorphism revealed within winter oilseed rape or “00” winter oilseed rape compared with spring oilseed rape.
A moderate level of differentiation was observed between the different B. napus types (WOSR, SOSR, fodder rape) as revealed by the estimated differentiation indices (0.18 between the three groups and 0.22 between SOSR and WOSR). A similar result was reported by Bus et al. . They used 89 SSR primer combinations to assess the diversity in a set of 509 oilseed rape lines which included WOSR, SOSR, fodder and swede rape lines from diverse origins. Xiao et al.  found a lower level of differentiation between their groups defined from a collection of 192 oilseed rape lines genotyped with 451 SSR markers. However, each of their groups, as defined by STRUCTURE analysis, was constituted of lines from China, Europe, Canada and Australia. Within the SOSR and WOSR groups, a very low level of differentiation was observed between the three subgroups which corresponded to the three quality types (“++”, “0+” and “00”) with 91% and 83% of the variation present within the subgroups in WOSR and SOSR, respectively. Examination of the collection structure with PCA showed distinct clustering of WOSR, SOSR and fodder rape lines while the first two axes did not account for a large part of the variation (19.5%). This differentiation was previously reported by Diers and Osborn , Hasan et al.  and Bus et al.  and can be related to the relatively distinct breeding history between these pools and their adaptation to different environments or uses. However, our data allowed two groups to be differentiated within SOSR which mainly contained European and Canadian or Asian spring oilseed rape, thus corresponding to a differentiation due to geographic origin. Some SOSR lines were located at an intermediate position between these two groups such as ‘Grouse’ and ‘Marnoo’ (Australian cvs), ‘Chine Wuhan’ and ‘Yeong Dang’ (Asiatic cvs) or ‘Industry’ (a European cv with high erucic acid and low glucosinolate content). No such differentiation was observed within the WOSR lines which all originated from Europe, as previously reported by Ecke et al.  and Bus et al. . Fodder rape lines were located at an intermediate position between SOSR and WOSR although closer to WOSR than to SOSR. The exception was ‘Liho’ which is a spring fodder rape and grouped with the SOSR. This result is consistent with the fact that fodder types and oilseed types were derived from the same ancestral spring and winter pools and were bred for different uses, as reported by Bus et al. . The same structuration was obtained whether we consider SNPs from the A or the C genomes. The percentage of variation accounted for by the first two axes was higher with the SNPs from the A than from the C genome when we used the whole SNP set where the markers from the A genome were overrepresented. When we considered the 1507 chosen SNP markers, the percentage of variation accounted for by the first two axes was similar with the SNPs from both genomes, indicating that an overrepresentation of markers in some regions can biased the results. A selection of evenly spaced markers over all the linkage groups is thus recommended for genetic structuration. The optimal number has to be assessed by testing different sample sizes.
To assess the genetic diversity in the collection and within the different germplasm groups, we examined the mean PIC and its evolution along each linkage group. The PIC values ranged between 0.1 and 0.35 depending on the position on the linkage groups and depending on the collections used. Similar observations were made in barley and maize with the same type of markers [4, 10, 68]. Our set of SNPs was only derived from exonic sequences, which could have lowered the level of revealed diversity compared to intronic SNPs . On average, PIC values were lower in WOSR than in SOSR. This difference was more important for SNPs derived from the A genome rather than the C genome. This might be because the SOSR lines have more diverse geographic origins than the WOSR lines and have undergone differential selection to adapt to the different continents. In WOSR, the mean PIC values were not very different between the three seed quality types. The mean PIC value was only slightly higher in “++” WOSR type but this result should be taken with caution due to the different sizes of the three subgroups. Bus et al.  reported a lower genetic diversity for “00” seed quality WOSR varieties and the difference was also quite low. The variation in PIC values along the different linkage groups between the three WOSR seed quality types was very contrasted, with more numerous regions showing a PIC decrease for the “0+” and/or “00” types. These variations could be related to potential selection signatures within and between these types. In barley, specific chromosomal regions exhibited contrasting levels of diversity in different germplasm subgroups. A region of reduced diversity in winter barley in the central part of chromosome 5H was attributed to the small number of founding genotypes that contributed to the winter seasonal growth habit locus Hrn-1. Similarly, an abrupt decrease in diversity on the short arm of chromosome 3H observed in all groups coincided with the locus of non-shattering of ears after ripening . We therefore investigated whether the regions with reduced diversity in either the “0+” and /or the “00” subgroups could be related to the position of genes controlling erucic acid and glucosinolate content. Erucic acid genes are located on A8 and C3. Numerous genes are involved in the glucosinolate pathway [69, 70] but the major QTL controlling total glucosinolate content are located on A2, A9, C2, C7 and C9 [71, 72], Delourme, unpubl. data. Their positions are indicated in Additional file 2: Table 2. However, there was no decrease in diversity clearly surrounding erucic acid genes in the “0+” and the “00” subgroups or total glucosinolate QTL in the “00” subgroup. The only exception could be at the top of C9. It can be hypothesized that many recombination events have occurred around the selected genes in these regions during different rounds of intercrossing between the varieties since the original crosses were made with the genitors of low erucic and low glucosinolate content. Significant decreases in diversity were observed in other regions, which could be related to breeding for other agronomical traits but a more precise investigation of QTL located in these regions should be made before drawing any conclusions.
The mean pairwise r2 values are close to previous estimates of 0.027  or 0.0247  and confirms the low overall level of linkage disequilibrium in B. napus. LD was observed to decay below a critical level (r2 value 0.2) within a map distance between 0.6 and 1.2 cM among the subgroups. This value is in accordance with previous studies performed either on smaller oilseed rape collections or with a smaller number of markers [20, 52]. This level is lower than that detected in a collection of 85 WOSR lines genotyped with 845 AFLP markers, where the LD decayed within 2–3 cM at r2 < 0.2 . This could be due to the fact that this latter collection comprised only “00” seed quality types. In our study, the extent of LD was also higher in the “00” WOSR collection than in the whole WOSR or the SOSR collections. From the size of the genetic maps and of the genomes in Brassica species, on average 1 cM can be roughly estimated to correspond to ~500 kb [51, 73]. This means that on average for the whole genome, the extent of LD is between 300–1000 kb depending on the collections but great variations were observed across the linkage groups (from ~100 kb on A10 to ~1700 kb on C9 in SOSR). In many species, LD decay varied across the germplasms used. For example, LD extends less than 1 kb for maize landraces  and roughly 2 kb for diverse inbred maize lines  but can be as high as 100 kb for commercial elite inbred lines . LD decay can also vary considerably from locus to locus 1–4 kb  up to 800 kb . Similar differences were observed in rice, 50–500 kb, [12, 15, 79], Arabidopsis, 10–250 kb [80, 81] or barley, 90–210 kb  or 4–8 cM . Generally, the extent of LD is related to the mating system of the species, the breeding history of the species (e.g. the occurrence of bottlenecks) and the genetic diversity of the different germplasms . LD decay is more rapid in outcrossing species and in pools with higher genetic diversity. Oilseed rape is bred as a selfing species and has undergone two bottlenecks for seed quality improvement (to eliminate erucic acid and decrease glucosinolate content), which led to a LD decay similar to other selfing species and to a higher LD extent in the “00” WOSR collection.
The LD patterns varied greatly among the linkage groups and these variations were different between the SOSR or WOSR types. Similar results were reported in maize [4, 5], barley  or rice . Different patterns of LD along the chromosomes in various pools can be related to variation in recombination rate and in the history of recombination for specific chromosome regions within these pools. The centromeres were localised approximately on the integrated map by mapping centromeric markers developed by Pouilly et al.  to the DYDH map. These were then included on the integrated map and on the LD heatmaps (Additional file 2: Table S2 and Additional file 9: Figure S5). On many linkage groups, LD seems to be extended across the centromeric regions as reported in barley  but other regions which do not correspond to centromeres also showed extended LD such as on A2, A8, C8 or C9. Differences in allele frequencies  could have also influenced the distribution of LD along the chromosomes. In barley, larger LD extent in some chromosomic regions was caused by markers with low allele frequencies . In rice, LD decay rates in indica and japonica subspecies were only weakly correlated across the genome in relation to a relatively long history of partial reproductive isolation of these self-fertilized subspecies . Since SOSR and WOSR have a relatively distinct breeding history, a similar hypothesis can also be proposed in our case.
With high throughput SNP genotyping on four segregating DH populations, we developed an integrated genetic map for oilseed rape that comprises 5764 SNP and 1603 PCR markers. This significantly improves the marker density or mapping accuracy compared to previously published genetic maps. The genotyping of these mapped SNP markers in collections allowed polymorphism level and linkage disequilibrium to be studied in oilseed rape. Both were shown to vary across the different collections (winter vs spring, seed quality types) and across the linkage groups. Taking into account the length of the genetic map (~2500 cM) and the mean LD extent (0.7 – 1.2 cM for r2 > 0.2), a relatively low number of evenly spaced SNPs (few thousands) would be necessary to perform genome wide association studies in oilseed rape. However, this number should be adjusted to obtain a sufficient SNP density throughout the genome and to take into account the variation in LD along the linkage groups. A set of 1536 public SNPs was set up, of which 72% were validated on a GoldenGate platform. They provide evenly spaced SNPs showing a good level of polymorphism in oilseed rape. Information regarding the other SNPs can be requested from J Pauquet (Jerome.Pauquet@biogemma.com). Our study provides a valuable resource for further genetic studies through linkage or association mapping, marker assisted breeding and Brassica sequence assembly and comparative mapping.
Four doubled haploid (DH) B. napus populations were used. Two have already been described in previous studies: DYDH [23, 31] and TNDH . Sets of 280 and 94 DH lines were used for these two populations, respectively. The third population, referred to as the AMDH population, was derived from the cross between two French winter oilseed rape varieties ‘Aviso’ and ‘Montego’ and consisted in 87 DH lines produced from the F1 between these two parents. The fourth population, referred to as the AADH population, was derived from the cross between a French winter oilseed rape ‘Aviso’ and a Japanese oilseed rape ‘Aburamasari’ and consisted in 96 DH derived from 192 intermated F2 plants: each F2 plant was used once, as male or female, in a cross with another F2 plant so that 96 hybrids were generated and one DH was derived per hybrid.
A B. napus collection of 313 inbred lines from different geographical origins was used for diversity analyses and linkage disequilibrium assessment. It consisted of 65 spring oilseed rape (SOSR) lines, 223 winter oilseed rape (WOSR) lines and 25 fodder rape lines from Europe, Australia, Asia and Canada. The SOSR and WOSR groups were divided in three subgroups depending on their seed quality types: “++” for high erucic acid and glucosinolate content, “0+” for low erucic acid and high glucosinolate content and “00” for double low types. A description of this collection is presented in Additional file 10: Table S5.
SNP origin and selection
Two sets of SNPs were used. The first set was obtained in previous internal research programs in Biogemma using sequence capture technology. Publicly available Brassica ESTs contigs (http://brassica.nbi.ac.uk/array_info.html) corresponding to a wide range of gene function were used to capture the corresponding genomic DNA in OSR genotypes including ‘Aviso’, ‘Montego’ and ‘Aburamasari’ as parent of mapping populations and Asiatic lines [88, 89]. A custom 2.1 M probes sequence capture Nimblegen (Roche NimbleGen, Inc., Madison, USA) microarray designed from those contigs was used with a protocol adapted from Albert et al. . Briefly, 454 sequencing libraries were synthesized, hybridized on microarray and subsequently, specifically hybridized library fragments were eluted and sequenced on a 454 GS-FLX sequencer. Reads were mapped against the targeted contigs and then assembled within each cluster using MIRA software  and SNP detection was performed using stringent criteria based on base quality, absence of heterozygosity within genotypes and 2X minimal allele coverage. The second set corresponds to SNPs identified between ‘Tapidor’ and ‘Ningyou7’ . All the SNPs were submitted to the Illumina Assay Design Tool (ADT) (Illumina, San Diego, CA), and only SNPs with designability scores > 0.4 for both Infinium and GoldenGate chemistries were included in further analyses.
A total of 7322 SNPs was selected for Infinium genotyping (4703 from the first set and 2619 from the second set). They were well distributed in silico over the Arabidopsis genome in order to obtain an even distribution of markers in the B. napus genome. The 7322 SNPs targeted 4190 EST contigs . To facilitate their use in a later GoldenGate genotyping assay, we also took care to select those SNPs that were at least 60 bp from another polymorphism. For SNPs derived from ‘Tapidor’ and ‘Ningyou7’, we also only considered SNPs that were at least 60 bp from an intron.
DNA was isolated from young leaves and DNA extracted using the DNeasy 96 Plant Kit (Qiagen, Courtaboeuf, France). DNA was quantified with the Quant-iT™ PicoGreen® Assay (Invitrogen, Carlsbad, USA), using the Appliskan multiplate reader (Thermo Scientific, Courtaboeuf, France). Concentrations were adjusted to a minimum of 50 ng/μL and were submitted to a provider, where the Infinium® assay was performed following the manufacturer's protocol (Illumina Inc., San Diego, USA). The automatic allele calling for each locus was accomplished using the Genome Studio software (Illumina Inc., San Diego, USA). The clusters were manually edited when necessary. Technical replicates and signal intensities were controlled and only the most reliable calls were retained.
Genetic maps construction
Individual genetic maps
For DYDH and TNDH, PCR markers that were previously genotyped [27, 28, 36] were added to the genotyping SNP matrix. Segregation of each marker was tested by Chi-square test for goodness of fit (1:1; P = 0.01) for the DYDH, AADH, AMDH and TNDH populations. At first, Mapmaker Exp/3.0b  was used to build a framework map for each individual genetic map. A minimum LOD score of 4.0 with a maximum genetic distance of 30 cM was first used to associate loci into initial linkage groups. A full multipoint linkage analysis was performed to determine the most probable locus order of highly informative markers (order with a LOD of 3.0 and with the highest log-likelihood ratio) for each linkage group. The remaining markers from each linkage group were manually integrated at their most likely position using the ‘try’ command. Double-crossover events were examined and the original scores rechecked for potential scoring errors. The order of the final set of framework loci within the linkage groups was re-verified using the ‘ripple’ command with a sliding window of five loci and a LOD score threshold of 3.0. Once the framework maps were built, the remaining markers were mapped using the ActionMap software . Each locus was mapped independently to the framework map, so low-quality mapping of some loci did not alter subsequent mapping of other loci and results referred to a stable reference map. All genetic distances were expressed in centimorgans using the Kosambi mapping function . Linkage groups were named according to the international Brassica nomenclature with A1- A10 and C1- C9 corresponding to the linkage groups of the A and C genomes, respectively. They were oriented as in Parkin et al. .
Integrated genetic map
Once the individual genetic maps were obtained, an integrated map was constructed through a projection process using the BioMercator V3.2.2 software . Using the iterative projection method, and in order to limit error propagation, the projection process started with the map that presents the more stable framework and well ordered and spread loci i.e. the map derived from the DYDH population.
Selection of a 1536 SNP set and validation in a GoldenGate assay
A set of 1536 SNPs were further chosen from the previous Infinium® assay to design four custom VeraCode OPA sets for the Illumina BeadXpress Reader. They were evenly distributed on the integrated map and showed a low number of missing data (<5%) and high minor allele frequencies (MAF > 10%) on the OSR panel. For each OPA run, a plate of 96 samples with 5 μL of genomic DNA normalized to 50 ng/μL was genotyped using the “GoldenGate Genotyping Assay for VeraCode Manual Protocol” (Illumina Inc., San Diego, USA). The 96 samples were a subset of the collection previously used for Infinium genotyping. Automatic allele calling for each locus was accomplished using the Genome Studio software (Illumina Inc., San Diego, USA). The clusters were manually edited when necessary.
Homology search with arabidopsis and B. Rapa
The sequences associated with each set of genetic markers were used as queries in homology searches against the Arabidopsis thaliana pseudo-chromosomes (TAIR10 release, ftp://ftp.arabidopsis.org/home/tair/Sequences/blast_datasets/TAIR10_blastsets/, version date: 16/04/2012), and against B. rapa Chiifu-401 pseudo-chromosomes in EnsemblPlants (IVF-CAASv1.14, ftp://ftp.ensemblgenomes.org/pub/plants/release-14/fasta/Brassica_rapa/dna/, version date: 27/05/2012). For homology searches on Arabidopsis thaliana, following parameters were used: TBlastX with match = 1, mismatch = −3, gap open penalty = −1, gap extension penalty = −1, word size = 3, and low complexity sequences filtered. A fairly low expect value (E-value) was used as the exclusion cutoff (1E-06). At least five consecutive homologous loci were required to define a collinearity block. Collinearity blocks were colour- coded according to the convention of Schranz et al. . For homology searches on B. rapa pseudo-chromosomes, following parameters were used: BlastN with match = 1, mismatch = −3, gap open penalty = 1, gap extension penalty = 2, word size = 7, and low complexity sequences filtered. A fairly low expect value (E-value) was used as the exclusion cutoff (1E-06).
The minor allele frequencies (MAF), percentage of heterozygosity and polymorphism information content (PIC) were estimated for each SNP marker within the different collections using PowerMarker v3.25 software . Mean PIC values were compared with Wilcoxon test (α = 5%). Population differentiation was studied using analysis of molecular variance (AMOVA) performed with Arlequin v3.1 software  and principal component analyses (PCA) performed with Darwin v5.0.158 . LD was estimated as the correlation coefficient r2 between all pairs of SNPs (with MAF > 5%) within and between linkage groups using PLINK v1.07 program . The overall decay of LD in relation to genetic distance was evaluated with R software (R development Core team, 2011) using the non linear regression of r2 according to Gaut and Long  with E[r2] = 1/(1 + 4Nec) where c is the recombination rate in Morgans and Ne the effective population size. LD heatmaps were built for each linkage group with the R package LDheatmap implemented in R software .
Paterson AH, Bowers JE, Bruggmann R, Dubchak I, Grimwood J, Gundlach H, Haberer G, Hellsten U, Mitros T, Poliakov A, Schmutz J, Spannagl M, Tang H, Wang X, Wicker T, Bharti AK, Chapman J, Feltus FA, Gowik U, Grigoriev IV, Lyons E, Maher CA, Martis M, Narechania A, Otillar RP, Penning BW, Salamov AA, Wang Y, Zhang L, Carpita NC, Freeling M, Gingle AR, Hash CT, Keller B, Klein P, Kresovich S, McCann MC, Ming R, Peterson DG, Rahman M, Ware D, Westhoff P, Mayer KFX, Messing J, Rokhsar DS: The Sorghum bicolor genome and the diversification of grasses. Nature. 2009, 457: 551-556.
Wei F, Zhang J, Zhou S, He R, Schaeffer M, Collura K, Kudrna D, Faga BP, Wissotski M, Golser W, Rock SM, Graves TA, Fulton RS, Coe E, Schnable PS, Schwartz DC, Ware D, Clifton SW, Wilson RK, Wing RA: The physical and genetic framework of the maize B73 genome. PLoS Genet. 2009, 5: e1000715-
Schmutz J, Cannon SB, Schlueter J, Ma J, Mitros T, Nelson W, Hyten DL, Song Q, Thelen JJ, Cheng J, Xu D, Hellsten U, May GD, Yu Y, Sakurai T, Umezawa T, Bhattacharyya MK, Sandhu D, Valliyodan B, Lindquist E, Peto M, Grant D, Shu S, Goodstein D, Barry K, Futrell-Griggs M, Abernathy B, Du J, Tian Z, Zhu L, Gill N, Joshi T, Libault M, Sethuraman A, Zhang XC, Shinozaki K, Nguyen HT, Wing RA, Cregan P, Specht J, Grimwood J, Rokhsar D, Stacey G, Shoemaker RC, Jackson SA: Genome sequence of the paleopolyploid soybean. Nature. 2010, 463: 178-183.
Lu Y, Shah T, Hao Z, Taba S, Zhang S, Gao S, Liu J, Cao M, Wang J, Prakash AB, Rong T, Xu Y: Comparative SNP and Haplotype analysis reveals a higher genetic diversity and rapider LD decay in tropical than temperate germplasm in maize. PLoS One. 2011, 6: e24861-
Van Inghelandt D, Reif JC, Dhillon BS, Flament P, Melchinger AE: Extent and genome-wide distribution of linkage disequilibrium in commercial maize germplasm. Theor Appl Genet. 2011, 123: 11-20.
Atwell S, Huang YS, Vilhja’lmsson BJ, Willems G, Horton M, Li Y, Meng D, Platt A, Tarone AM, Hu TT, Jiang R, Muliyati NW, Zhang X, Amer MA, Baxter I, Brachi B, Chory J, Dean C, Debieu M, de Meaux J, Ecker JR, Faure N, Kniskern JM, Jones JDG, Michael T, Nemri A, Roux F, Salt DE, Tang C, Todesco M: Genome-wide association study of 107 phenotypes in Arabidopsis thaliana inbred lines. Nature. 2010, 465: 627-631.
Hamblin MT, Close TJ, Bhat PR, Chao S, Kling JG, Abraham KJ, Blake T, Brooks WS, Cooper B, Griffey CA, Hayes PM, Hole DJ, Horsley RD, Obert DE, Smith KP, Ullrich SE, Muehlbauer GJ, Jannink J-L: Population structure and linkage disequilibrium in U.S. Barley germplasm: implications for association mapping. Crop Sci. 2010, 50: 556-566.
Lorenz AJ, Hamblin MT, Jannink J-L: Performance of single nucleotide polymorphisms versus haplotypes for genome-wide association analysis in barley. PLoS One. 2010, 5: e14079-
Bradbury P, Parker T, Hamblin MT, Jannink J-L: Assessment of power and false discovery rate in genome-wide association studies using the BarleyCAP germplasm. Crop Sci. 2011, 51: 52-59.
Comadran J, Russell JR, Booth A, Pswarayi A, Ceccarelli S, Grando S, Stanca AM, Pecchioni N, Akar T, Al-Yassin A, Benbelkacem A, Ouabbou H, Bort J, van Eeuwijk FA, Thomas WTB, Romagosa I: Mixed model association scans of multi-environmental trial data reveal major loci controlling yield and yield related traits in Hordeum vulgare in Mediterranean environments. Theor Appl Genet. 2011, 122: 1363-1373.
Ehrenreich IM, Hanzawa Y, Chou L, Roe JL, Kover PX, Purugganan MD: Candidate gene association mapping of Arabidopsis flowering time. Genetics. 2011, 183: 325-335.
McNally KL, Childs KL, Bohnert R, Davidson RM, Zhao K, Ulat VJ, Zeller G, Clark RM, Hoen DR, Bureau TE, Stokowski R, Ballinger DG, Frazer KA, Cox DR, Padhukasahasram B, Bustamant CD, Weigel D, Mackill DJ, Bruskiewich RM, Rätsch G, Buell CR, Leung H, Leach JE: Genomewide SNP variation reveals relationships among landraces and modern varieties of rice. Proc Natl Acad Sci USA. 2009, 106: 12273-12278.
Huang X, Wei X, Sang T, Zhao Q, Feng Q, Zhao Y, Li C, Zhu C, Lu T, Zhang Z, Li M, Fan D, Guo Y, Wang A, Wang L, Deng L, Li W, Lu Y, Weng Q, Liu K, Huang T, Zhou T, Jing Y, Li W, Lin Z, Buckler ES, Qian Q, Zhang Q-F, Li J, Han B: Genome-wide association studies of 14 agronomic traits in rice landraces. Nat Genet. 2010, 42: 961-969.
Famoso AN, Zhao K, Clark RT, Tung C-W, Wright MH, Bustamante C, Kochian LV, McCouch SR: Genetic architecture of aluminum tolerance in rice (Oryza sativa) determined through genome-wide association analysis and QTL mapping. PLoS Genet. 2011, 7: e1002221-
Zhao K, Tung C-W, Eizenga GC, Wright MH, Ali ML, Price AH, Norton GJ, Islam MR, Reynolds A, Mezey J, McClung AM, Bustamante CD, McCouch SR: Genome-wide association mapping reveals a rich genetic architecture of complex traits in Oryza sativa. Nature Comm. 2011, 2: 467-
U N: Genomic analysis in Brassica with special reference to the experimental formation of B. napus and peculiar mode of fertilization. Japan J Bot. 1935, 7: 389-452.
Gomez-Campo C, Prakash S: Origin and Domestication. Biology of Brassica coenospecies. Edited by: Gomez-Campo C. 1999, Amsterdam, The Netherlands: Elsevier Science B.V, 33-58.
Allender CJ, King GJ: Origin of the amphiploid species Brassica napus L. Investigated by chloroplast and nuclear molecular markers. BMC Plant Biol. 2010, 10: 54-
Toxopeus H: The domestication of Brassica crops. 1979, Wageningen: Proc Eucarpia Conference on the breeding of Cruciferous crops, 47-56.
Bus A, Körber N, Snowdon RJ, Stich B: Patterns of molecular variation in a species-wide germplasm set of Brassica napus. Theor Appl Genet. 2011, 123: 1413-1423.
Parkin IA, Sharpe AG, Keith DJ, Lydiate DJ: Identification of the A and C genomes of amphidiploid Brassica napus (oilseed rape). Genome. 1995, 38: 1122-1131.
Parkin IAP, Gulden SM, Sharpe AG, Lukens L, Trick M, Osborn TC, Lydiate DJ: Segmental structure of the Brassica napus genome based on comparative analysis with Arabidopsis thaliana. Genetics. 2005, 171: 765-781.
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: 491-507.
Lowe AJ, Moule C, Trick M, Edwards KJ: Efficient large-scale development of microsatellites for marker and mapping applications in Brassica crop species. Theor Appl Genet. 2004, 108: 1103-1112.
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: 1514-1523.
Suwabe K, Morgan C, Bancroft I: Integration of Brassica A genome genetic linkage map between Brassica napus and B. rapa. Genome. 2008, 51: 169-176.
Wang J, Lydiate DJ, Parkin IAP, Falentin C, Delourme R, Carion PWC, King GJ: Integration of linkage maps for the amphidiploid Brassica napusand comparative mapping with Arabidopsis and Brassica rapa. BMC Genomics. 2011, 12: 101-
Long Y, Shi J, Qiu D, Li R, Zhang C, Wang J, Hou J, Zhao J, Shi L, Park B-S, Choi SR, Lim YP, Meng J: Flowering time quantitative trait loci analysis of oilseed Brassica in multiple environments and genomewide alignment with Arabidopsis. Genetics. 2007, 177: 2433-2444.
Mei DS, Wang HZ, Hu Q, Li YD, Xu YS, Li YC: QTL analysis on plant height and flowering time in Brassica napus. Plant Breeding. 2009, 128: 458-465.
Uzunova M, Ecke W, Weissleder K, Robbelen G: Mapping the genome of rapeseed (Brassica napusL).1. Construction of an RFLP linkage map and localization of QTLs for seed glucosinolate content. Theor Appl Genet. 1995, 90: 194-204.
Delourme R, Falentin C, Huteau V, Clouet V, Horvais R, Gandon B, Specel S, Hanneton L, Dheu JE, Deschamps M, Margale E, Vincourt P, Renard M: Genetic control of oil content in oilseed rape (Brassica napus L.). Theor Appl Genet. 2006, 113: 1331-1345.
Zhao J, Becker HC, Zhang D, Zhang Y, Ecke W: Conditional QTL mapping of oil content in rapeseed with respect to protein content and traits related to plant development and grain yield. Theor Appl Genet. 2006, 113: 33-38.
Qiu D, Morgan C, Shi J, Long Y, Liu J, Li R, Zhuang X, Wang Y, Tan X, Dietrich E, Weihmann T, Everett C, Vanstraelen S, Beckett P, Fraser F, Trick M, Barnes S, Wilmer J, Schmidt R, Li J, Li D, Meng J, Bancroft I: 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.
Chen G, Geng J, Rahman M, Liu X, Tu J, Fu T, Li G, McVetty PBE, Tahir M: Identification of QTL for oil content, seed yield, and flowering time in oilseed rape (Brassica napus). Euphytica. 2010, 175: 161-174.
Radoev M, Becker HC, Ecke W: Genetic analysis of heterosis for yield and yield components in rapeseed (Brassica napus L.) by quantitative trait locus mapping. Genetics. 2008, 179: 1547-1558.
Shi J, Li R, Qiu D, Jiang C, Long Y, Morgan C, Bancroft I, Zhao J, Meng J: Unraveling the complex trait of crop yield with quantitative trait loci mapping in Brassica napus. Genetics. 2009, 182: 851-861.
Basunanda P, Radoev M, Ecke W, Friedt W, Becker HC, Snowdon RJ: Comparative mapping of quantitative trait loci involved in heterosis for seedling and yield traits in oilseed rape (Brassica napus L.). Theor Appl Genet. 2010, 120: 271-281.
Wei C, Yongshan Z, Jinbo Y, Chaozhi M, Jinxing T, Fu T: Quantitative trait loci mapping for two seed yield component traits in an oilseed rape (Brassica napus) cross. Plant Breeding. 2011, 130: 640-646.
Pilet ML, Duplan G, Archipiano M, Barret P, Baron C, Horvais R, Tanguy X, Lucas MO, Renard M, Delourme R: Stability of QTL for field resistance to blackleg across two genetic backgrounds in oilseed rape. Crop Sci. 2001, 41: 197-205.
Manzanares-Dauleux MJ, Delourme R, Baron F, Thomas G: Mapping of one major gene and of QTLs involved in resistance to clubroot in Brassica napus. Theor Appl Genet. 2000, 101: 885-891.
Zhao J, Udall JA, Quijada PA, Grau CR, Meng J, Osborn TC: Quantitative trait loci for resistance to Sclerotinia sclerotiorum and its association with a homeologous non-reciprocal transposition in Brassica napus L. Theor Appl Genet. 2006, 112: 509-516.
Werner S, Diederichsen E, Frauen M, Schondelmaier J, Jung C: Genetic mapping of clubroot resistance genes in oilseed rape. Theor Appl Genet. 2008, 116: 363-372.
Rygulla W, Showdon RJ, Friedt W, Happstadius I, Cheung WY, Chen D: Identification of quantitative trait loci for resistance against Verticillium longisporum in oilseed rape (Brassica napus). Phytopathology. 2008, 98: 215-221.
Kaur S, Cogan NOI, Ye G, Baillie RC, Hand ML, Ling AE, Mcgearey AK, Kaur J, Hopkins CJ, Todorovic M, Mountford H, Edwards D, Batley J, Burton W, Salisbury P, Gororo N, Marcroft S, Kearney G, Smith KF, Forster JW, Spangenberg GC: Genetic map construction and QTL mapping of resistance to blackleg (Leptosphaeria maculans) disease in Australian canola (Brassica napus L.) cultivars. Theor Appl Genet. 2009, 120: 71-83.
Liu Z, Adamczyk K, Manzanares-Dauleux MJ, Eber F, Lucas M-O, Delourme R, Chèvre A-M, Jenczewski E: Mapping PrBn and other quantitative trait loci responsible for the control of homeologous chromosome pairing in oilseed rape (Brassica napus L.) haploids. Genetics. 2006, 174: 1583-1596.
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: 535-542.
Panjabi P, Jagannath A, Bisht NC, Padmaja KL, Sharma S, Gupta V, Pradhan AK, Pental D: Comparative mapping of Brassica juncea and Arabidopsis thaliana using intron polymorphism (IP) markers: homoeologous relationships, diversification and evolution of the A, B and C Brassica genomes. BMC Genomics. 2008, 9: 113-
Kaczmarek M, Koczyk G, Ziolkowski PA, Babula-Skowronska D, Sadowski J: Comparative analysis of the Brassica oleracea genetic map and the Arabidopsis thaliana genome. Genome. 2009, 52: 620-633.
Parkin IAP: Chasing ghosts: Comparative mapping in the Brassicaceae. Genetics and Genomics of the Brassicaceae, Plant Genetics and Genomics: Crop and Models 9. Edited by: Schmidt R, Bancroft I. 2011, New-York: Springer Science+Business Media, 153-170.
Bancroft I, Morgan C, Fraser F, Higgins J, Wells R, Clissold L, Baker D, Long Y, Meng J, Wang X, Liu S, Trick M: Dissecting the genome of the polyploid crop oilseed rape by transcriptome sequencing. Nat Biotechnol. 2011, 29: 762-768.
Ecke W, Clemens R, Honsdorf N, Becker HC: Extent and structure of linkage disequilibrium in canola quality winter rapeseed (Brassica napus L.). Theor Appl Genetics. 2010, 120: 921-931.
Xiao Y, Cai D, Yang W, Ye W, Younas M, Wu J, Liu K: Genetic structure and linkage disequilibrium pattern of a rapeseed (Brassica napus L.) association mapping panel revealed by microsatellites. Theor Appl Genet Online First. 2012, 125: 437-447.
Harper AL, Trick M, Higgins J, Fraser F, Clissold L, Wells R, Hattori C, Werner P, Bancroft I: Associative transcriptomics of traits in the polyloid crop species Brassica napus. Nat Biotechnol. 2012, 30: 798-802.
Kaur S, Francki MG, Forster JW: Identification, characterization and interpretation of single-nucleotide sequence variation in allopolyploid crop species. Plant Biotechnol J. 2012, 10: 125-138.
Durstewitz G, Polley A, Plieske J, Luerssen H, Graner EM, Wieseke R, Ganal MW: SNP discovery by amplicon sequencing and multiplex SNP genotyping in the allopolyploid species Brassica napus. Genome. 2010, 53: 948-956.
Wiggins M, Tang S, Bai Y, Lu F, Powers C, Pita F, Ubayasena L, Ehlert Z, Kubik T, Gingera G, Stoll C, Ripley V, Greene T, Thompson S, Kumpatla S: High-throughput single nucleotide polymorphism (SNP) discovery and marker validation in Brassica napus. 2010, Saskatoon, Canada: Proceedings of the 17th Crucifer Genetics Workshop (Brassica 2010), 140-Abstracts,
Hu Z, Huang S, Sun M, Wang H, Hua W: Development and application of single nucleotide polymorphism markers in the polyploid Brassica napus by 454 sequencing of expressed sequence tags. Plant Breeding. 2012, 131: 293-299.
Sidebottom C, Clarke W, Li R, Higgins E, Links M, Parkin I, Sharpe A: SNP discovery in Brassica napus using 454 titanium sequencing. 2010, San Diego, California, USA: Proceedings of Plant and Animal Genome XVIII, P378-
Tang S, Wiggins M, Nipper R, Gribbin J, Johnson E, Lillegard N, Greene T, Thompson S, Kumpatla S: SNP discovery using restriction site-associated DNA (RAD) LongRead sequencing in Brassica napus, a polyploid species. 2010, Saskatoon, Canada: Proceedings of the 17th Crucifer Genetics Workshop (Brassica 2010), 109-Abstracts
Bus A, Hecht J, Huettel B, Reinhardt R, Stich B: High-throughput polymorphism detection and genotyping in Brassica napus using next-generation RAD sequencing. BMC Genomics. 2012, 13: 281-
Shiga T: Rape breeding by interspecific crossing between Brassica napus and Brassica campestris in Japan. Japan Agric Res Quaterly. 1970, 5: 5-10.
Qian W, Meng J, Li M, Frauen M, Sass O, Noack J, Jung C: Introgression of genomic components from Chinese B. rapa contributes to widening the genetic diversity in rapeseed (B. napus L.), with emphasis on the evolution of Chinese rapeseed. Theor Appl Genet. 2006, 113: 49-54.
Foisset N, Delourme R: Segregation distortion in androgenic plants. In vitro haploid production in higher plants. Edited by: Mohan Jain S, Sopory S, Veilleux R. 1996, Dordrecht, The Netherlands: Kluwer Academic Publishers, 189-201.
Liu S-C, Kowalski SP, Lan T-H: Genome-wide high-resolution mapping by recurrent intermating using Arabidopsis thaliana as a model. Genetics. 1996, 142: 247-258.
Lee M, Sharopova N, Beavis WD, Grant D, Katt M, Blair D, Hallauer A: Expanding the genetic map of maize with the intermated B73 × Mo17 (IBM) population. Plant Mol Biol. 2002, 48: 453-461.
Diers BW, McVetty PBE, Osborn TC: Relationship between heterosis and genetic distance based on restriction fragment length polymorphism markers in oilseed rape (Brassica napus L). Crop Sci. 1996, 36: 79-83.
Hasan M, Seyis F, Badani AG, Pons-Kühnemann J, Friedt W, Lühs W, Snowdon RJ: Analysis of genetic diversity in the Brassica napus L. Gene pool using SSR markers. Genet Res Crop Evol. 2006, 53: 793-802.
Rostoks N, Ramsay L, MacKenzie K, Cardle L, Bhat PR, Roose ML, Svensson JT, Stein N, Varshney RK, Marshall DF, Graner A, Close TJ, Waugh R: Recent history of artificial outcrossing facilitates whole-genome association mapping in elite inbred crop varieties. Proc Natl Acad Sci U S A. 2006, 103: 18656-18661.
Halkier BA, Gershenzon J: Biology and biochemistry of glucosinolates. Annu Rev Plant Biol, Book Series: Annu Rev Plant Biol. 2006, 57: 303-333.
Feng J, Long Y, Shi L, Shi J, Barker G, Meng J: Characterization of metabolite quantitative trait loci and metabolic networks that control glucosinolate concentration in the seeds and leaves of Brassica napus. New Phytol. 2011, 193: 96-108.
Howell PM, Sharpe AG, Lydiate DJ: Homoeologous loci control the accumulation of seed glucosinolates in oilseed rape (Brassica napus). Genome. 2003, 46: 454-460.
Quijada PA, Udall JA, Lambert B, Osborn TC: Quantitative trait analysis of seed yield and other complex traits in hybrid spring rapeseed (Brassica napus L.): 1. Identification of genomic regions from winter germplasm. Theor Appl Genet. 2006, 113: 549-561.
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: 309-319.
Tenaillon MI, Sawkins MC, Long AD, Gaut RL, Doebley JF, Gaut BS: Patterns of DNA sequence polymorphism along chromosome 1 of maize (Zea mays ssp. mays L.). Proc Natl Acad Sci U S A. 2001, 98: 9161-9166.
Remington DL, Thornsberry JM, Matsuoka Y, Wilson LM, Whitt SR, Doebley J, Kresovich S, Goodman MM, Buckler ES: Structure of linkage disequilibrium and phenotypic associations in the maize genome. Proc Natl Acad Sci U S A. 2001, 98: 11479-11484.
Ching A, Caldwell KS, Jung M, Dolan M, Smith OS, Tingey S, Morgante M, Rafalski AJ: SNP frequency, haplotype structure and linkage disequilibrium in elite maize inbred lines. BMC Genet. 2002, 3: 19-
Palaisa KA, Morgante M, Williams M, Rafalski A: Contrasting effects of selection on sequence diversity and linkage disequilibrium at two phytoene synthase loci. Plant Cell. 2003, 15: 1795-1806.
Palaisa K, Morgante M, Tingey S, Rafalski A: Long-range patterns of diversity and linkage disequilibrium surrounding the maize Y1 gene are indicative of an asymmetric selective sweep. Proc Natl Acad Sci U S A. 2004, 101: 9885-9890.
Garris AJ, McCouch SR, Kresovich S: Population structure and its effect on haplotype diversity and linkage disequilibrium surrounding the xa5 locus of rice (Oryza sativa L.). Genetics. 2003, 165: 759-769.
Nordborg M, Borevitz JO, Bergelson J, Berry CC, Chory J, Hangenblag J, Kreitman M, Maloof JN, Noyes T, Oefner PJ, Stahl EA, Weigel D: The extent of linkage disequilibrium in Arabidopsis thaliana. Nat Genet. 2002, 30: 190-193.
Nordborg M, Hu TT, Ishino Y, Jhaveri J, Toomajian C, Zheng H, Bakker E, Calabrese P, Gladstone J, Goyal R, Jakobsson M, Kim S, Morozov Y, Padhukasahasram B, Plagnol V, Rosenberg NA, Shah C, Wall JD, Wang J, Zhao K, Kalbfleisch T, Schulz V, Kreitman M, Bergelson J: The pattern of polymorphism in Arabidopsis thaliana. PLoS Biol. 2005, 3: e196-
Caldwell KS, Russell J, Langridge P, Powell W: Extreme population-dependent linkage disequilibrium detected in an inbreeding plant species, Hordeum vulgare. Genetics. 2006, 172: 557-567.
Massman J, Cooper B, Horsley R, Neate S, Dill-Macky R, Chao S, Dong Y, Schwarz P, Muehlbauer GJ, Smith KP: Genome-wide association mapping of Fusarium head blight resistance in contemporary barley breeding germplasm. Mol Breed. 2011, 27: 439-454.
Flint-Garcia SA, Thornsberry JM, Buckler ES: Structure of linkage disequilibrium in plants. Annual Rev Plant Biology. 2003, 54: 357-374.
Pasam RK, Sharma R, Malosetti M, van Eeuwijk FA, Haseneyer G, Kilian B, Graner A: Genome-wide association studies for agronomical traits in a worldwide spring barley collection. BMC Plant Biol. 2012, 12: 16-
Pouilly N, Delourme R, Alix K, Jenczewski E: Repetitive sequence-derived markers tag centromeres and telomeres and provide insights into chromosome evolution in Brassica napus. Chromosome Res. 2008, 16: 683-700.
Lewontin RC: On measures of gametic disequilibrium. Genetics. 1988, 120: 849-852.
Pichon J-P, Rivière N, Duarte J, Dugas O, Wilmer JA, Gerhardt DJ, Richmond T, Albert TJ, Jeddeloh JA: Rapeseed (B. napus) SNP discovery using a dedicated sequence capture protocol and 454 Sequencing. 2010, San Diego, USA: Proceedings of Plant and Animal Genome Conference, W643-
Faure S, Throude M, Duarte J, Pichon JP, Pauquet J, Rivière N: Wheat and Rapeseed coming to the age of high-throughput SNP discovery. 2011, San Diego, USA: Proceedings of Plant and Animal Genome Conference, P168-
Albert TJ, Molla MN, Muzny DM, Nazareth L, Wheeler D, Song XZ, Richmond TA, Middle CM, Rodesch MJ, Packard CJ, Weinstock GM Gibbs RA: Direct selection of human genomic loci by microarray hybridization. Nat Methods. 2007, 4: 903-905.
Chevreux B, Pfisterer T, Drescher B, Driesel AJ, Müller WE, Wetter T, Suhai S: Using the miraEST assembler for reliable and automated mRNA transcript assembly and SNP detection in sequenced ESTs. Genome Res. 2004, 14: 1147-1159.
Trick M, Cheung F, Drou N, Fraser F, Lobenhofer EK, Hurban P, Magusin A, Town CD, Bancroft I: A newly-developed community microarray resource for transcriptome profiling in Brassica species enables the confirmation of Brassica-specific expressed sequences. BMC Plant Biol. 2009, 50:
Lander ES, Green P, Abrahamson J, Barlow A, Daly MJ, Lincoln SE, Newberg LA: MAPMAKER: an interactive computer package for constructing primary genetic linkage maps of experimental and natural populations. Genomics. 1987, 1: 174-181.
Albini G, Falque M, Joets J: ActionMap: a web-based software that automates loci assignments to framework maps. Nucleic Acids Res. 2003, 31: 3815-3818.
Kosambi DD: The estimation of map distance from recombination values. Ann Eugen. 1944, 12: 172-175.
Arcade A, Labourdette A, Falque M, Mangin B, Chardon F, Charcosset A, Joets J: BioMercator: integrating genetic maps and QTL towards discovery of candidate genes. Bioinformatics. 2004, 20: 2324-2326.
Liu KJ, Muse SV: PowerMarker: an integrated analysis environment for genetic marker analysis. Bioinformatics. 2005, 21: 2128-2129.
Excoffier L, Laval G, Schneider S: Arlequin (version 3.0): An integrated software package for population genetics data analysis. Evol Bioinforma. 2005, 1: 47-50.
Perrier X, Flori A, Bonnot F: Data analysis methods. In: Genetic diversity of cultivated tropical plants. Edited by: Hamon P, Seguin M, Perrier X, Glaszmann JC. 2003, Enfield: Science Publishers, 43-76. Montpellier.
Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MAR, Bender D, Maller J, Sklar P, de Bakker PIW, Daly MJ, Sham PC: PLINK: a tool set for whole-genome association and population-based linkage analyses. American J Human Genet. 2007, 81: 559-575.
Gaut BS, Long AD: The lowdown on linkage disequilibrium. Plant Cell. 2003, 15: 1502-1506.
Shin J-H, Blay S, McNeney B, Graham J: LDheatmap: an R function for graphical display of pairwise linkage disequilibria between single nucleotide polymorphisms. J Stat Soft. 2006, 16: 1-10.
We acknowledge the financial support of SOFIPROTEOL under the FSRSO (Fond de Soutien à la Recherche Semencière Oléagineuse) project “SNP colza” and of INRA (Institut National de la Recherche Agronomique) under a AIP BioRessources grant. We thank Anne Laperche for critical reading of the manuscript. We acknowledge the Genetic Resource Center (BrACySol, UMR IGEPP Ploudaniel, France) for providing the seeds of some B. napus varieties.
The authors declare no competing interests.
RD contributed to the design of the study, performed genetic diversity and linkage disequilibrium analyses and wrote the manuscript. CF carried out genetic mapping and BLAST analyses. BFF contributed to the genetic diversity and linkage disequilibrium studies. GL contributed to genetic mapping and BLAST analyses. JPP, JD, IA and NR selected the SNPs after sequence capture studies. NR carried out genetic mapping. CF, VG, NL, AM, MP, GT acquired the genotyping data. MB and JPM organized the production of genotyping data. PB contributed to the design of the study. JPM coordinated the project. JP helped analyze the results and participated in the coordination of the study. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Table S1: Description of the individual genetic maps obtained on the ‘Darmor-bzh’ x ‘Yudal’ (DYDH), ‘Tapidor’ x Ningyou7’ (TNDH), ‘Aviso’ x ‘Abumasari’ (AADH) and ‘Aviso’ x ‘Montego’ (AMDH) DH populations. Total number of markers and SNPs per LG and on the A, C and whole genomes and ratio of marker and SNP number on the A genome relatively to the C genome (A/C) are indicated. (XLSX 21 KB)
Additional file 2: Table S2: Summary of data obtained for all the markers: assignment to linkage groups, collinearity blocks, centromere, erucic genes and glucosinolate QTL positions, position on the integrated and the individual maps, segregation distortion and PIC values on the different collections. (XLSX 1 MB)
Additional file 3: Figure S1: Alignment of the four individual maps obtained for the TNDH, DYDH, AADH and AMDH populations and the integrated map (on the right). The blocks as defined by Schranz et al.  based on their collinearity with A. thaliana are indicated with the capital letters A to X on the left of each linkage group. (PDF 826 KB)
Additional file 4: Figure S2: Dot-plots obtained for each linkage group between the integrated map and all four individual TNDH, DYDH, AADH and AMDH maps. The marker order on the vertical axis is from the four individual maps and the marker order on the horizontal axis is from the integrated map. Cumulative genetic distance in cM is indicated on each axis. (PDF 471 kb) (PDF 472 KB)
Additional file 5: Table S3: List of the 222 SNPs that were assigned to different linkage groups (LG) according to the mapping populations. (XLSX 59 KB)
Additional file 6: Table S4: List and characteristics of the 1536 public SNPs : Position on the integrated map (or on the individual maps for multi-loci SNPs), Context sequence and SNP position, Arabidopsis and B. rapa hits. (XLSX 289 KB)
Additional file 7: Figure S3: Distribution of PIC values along the linkage groups in the fodder rape (FO), spring (SOSR) and winter (WOSR) oilseed rape and the different seed quality subgroups (“++”,” 0+”,” 00”) within WOSR. PIC was averaged across a sliding window of 10 cM with a step of one cM. (PDF 351 KB)
Additional file 8: Figure S4: Plots of r2 as a function of genetic distance (in cM) between pairs of SNPs on each linkage group in the whole, spring, winter and “00” winter oilseed rape collections. Red curves show the non-linear regressions trend line of r2 versus genetic distance. (PDF 1 MB)
Additional file 9: Figure S5: Linkage disequilibrium heatmaps for each linkage group in the whole, spring (SOSR), winter (WOSR) and “00” WOSR collections. The putative position of the centromeres is indicated by a black arrow. (PDF 5 MB)
About this article
Cite this article
Delourme, R., Falentin, C., Fomeju, B.F. et al. High-density SNP-based genetic map development and linkage disequilibrium assessment in Brassica napusL. BMC Genomics 14, 120 (2013). https://doi.org/10.1186/1471-2164-14-120
- Linkage Disequilibrium
- Linkage Group
- Double Haploid
- Oilseed Rape
- Erucic Acid