- Research article
- Open Access
Quantitative trait loci analysis and genome-wide comparison for silique related traits in Brassica napus
- Xiaodong Wang†1, 2,
- Li Chen†1, 3,
- Aina Wang†4,
- Hao Wang4Email author,
- Jianhua Tian4,
- Xiaoping Zhao4,
- Hongbo Chao1,
- Yajun Zhao4,
- Weiguo Zhao4,
- Jun Xiang3,
- Jianping Gan3 and
- Maoteng Li1Email author
© Wang et al. 2016
Received: 13 November 2015
Accepted: 15 March 2016
Published: 22 March 2016
Yield of rapeseed is determined by three components: silique number, seed number per silique and thousand seed weight. Seed number per silique and thousand seed weight are influenced by silique length, seed density, silique breadth, silique thickness and silique volume. Some QTLs for silique traits have been reported in B. napus, however, no studies have focused on the six agronomic traits (seed number per silique, silique length, silique breadth, silique thickness, seed density and silique volume) simultaneously, and the genetic determinism of such complex traits have not been fully elucidated.
In this study, the six silique traits were evaluated using 348 lines of a doubled haploid population, the KN population. The results showed that 2, 4, 1, 1 and 2 QTLs explaining > 10 % of phenotypic variation were obtained for silique length, silique breadth, silique thickness, seed number per silique and silique volume, respectively. Notably, three major effect QTLs (cqSB-C6-1, cqSB-C6-2 and cqSV-C6-3) were identified in at least three environments, and 17 unique QTLs controlling at least two traits were obtained. A high-density consensus map containing 1225 markers was constructed for QTL comparison by combining the KN map with other five published maps. The comparative results revealed that 14, 13 and 11 QTLs for silique breadth, silique thickness and silique volume might be the potential new QTLs because few QTLs for these traits were reported in B. napus. In addition, potential new QTLs for silique length (11), seed number per silique (6) and seed density (5) were also identified. Twenty-five candidate genes underlying 27 QTLs for silique related traits were obtained.
This study constructed QTL analysis in B. napus, and obtained 60 consensus QTLs for six silique related traits. The potential new QTLs will enhance our understanding of the genetic control of silique traits, and the stable QTLs provided the targets for improving seed yield in future. These findings provided comprehensive insights into the genetic network affecting silique traits at QTL level in B. napus.
Brassica napus (AACC, 2n = 38) is cultivated all over the world for the production of vegetable oil, animal feed and biodiesel. At present, developing high-yielding cultivars is one of the most important tasks of rapeseed breeders. Yield of rapeseed is determined by three components: silique number (SN), seed number per silique (SPS) and thousand seed weight (TSW) . SN and SPS determine the total number of seeds per plant. Path analysis identified that TSW has the greatest effect on seed yield, followed by the number of pods per plant . SPS and TSW are influenced by silique length (SL), seed density (SD), silique breadth (SB), silique thickness (ST) and silique volume (SV). Silique related traits showed significant relationships between each other, such as SL has a significant positive relationship with SPS [3, 4]. Therefore, improving seed yield through coordination between these traits is an important breeding goal.
The silique traits are all complex quantitative traits controlled by polygenes and highly influenced by environmental conditions. Application of molecular markers for QTL mapping has proved to be a powerful genetic approach to dissect quantitative traits. In recent years, QTL analysis of silique traits, such as SL, TSW, SPS and SD has been carried out using different populations in B. napus [1, 3–6]. More recently, a major QTL on chromosome A9 simultaneously affects TSW and SL in B. napus was isolated by fine mapping and association analysis . Li et al. (2015) report the cloning and characterization of BnaC9.SMG7b on C9, a major QTL that controls SPS in B. napus . In addition, SB, ST and SV are also important silique characteristics for seed yield. Adamskia et al. (2009) found that seed size correlated with oil and protein contents, and large seeds normally have better adaptability during germination . Larger values of SL, SB and ST lead to larger SV, and this may lead to larger seed size and stronger photosynthesis. However, there has been little research concerning SB, ST and SV so far.
Several consensus maps have been constructed by integrating different linkage maps based on common molecular markers [10–13], and QTL hotspots for many agronomic traits in B. napus have been identified, including flowering time , seed yield , oil content [11, 15], heterosis-related traits  and plant height . Zhou et al. (2014) carried out in silico mapping to integrate 1960 QTLs, which increased the density of targeted QTL-linked markers and validated the stable existence of QTLs across different populations . Comparison of QTL position and validation of conserved QTLs through map integration is of great significance for breeders to utilize QTLs effectively.
Although QTLs for silique traits have been identified [10–13], the candidate genes related to these QTLs have not been fully elucidated. Comparative mapping among related species is a powerful tool for genetic studies of transferring genomic information from the well-studied species to those more genetically complicated . Nowadays, much effort focused on comparative analysis between Brassica and Arabidopsis. For example, Raman et al. (2015) conduct a genome-wide association study for flowering time in B. napus and identify seven single nucleotide polymorphism markers were detected within 20 Kb regions of Arabidopsis candidate genes, including BnFLC.A2 accounts for ~23 % of natural variation in diverse accessions . Numerous of predicted genes for other traits were obtained, such as phosphorus homeostasis  and SPS . To facilitate the utilization of QTLs for the Brassica community, comparative mapping among B. napus, Arabidopsis, B. rapa and B. oleracea genomes is necessary to obtain candidate genes, and the released B. rapa, B. oleracea and B. napus genome information may accelerate map-based cloning [23–25].
A large population, a high-density genetic map and replicated experiments in multiple environments are considered as three key factors for increasing statistical power and precision in detecting QTLs . In our previous study, KN doubled haploid (DH) population was derived from a cross between ‘KenC-8’ and ‘N53-2’, which contained 348 lines . This population also showed high phenotypic variation for silique related traits. Using KN DH population, the aim of this study was to detect QTLs for SL, SB, ST, SPS, SV and SD in multiple environments, and to study the relationships among these traits (including TSW). We also constructed a consensus map to perform QTL comparisons between different populations, and identified candidate genes through comparative genome analysis.
Trait variation and correlation analysis among different silique traits
Phenotypic correlations among silique related traits, seed yield and yield-related traits in KN population
QTL mapping for SL, SB, ST, SPS, SV and SD
Consensus QTL obtained for the six silique related traits
-0.27 ~ -0.22
-0.29 ~ -0.21
0.20 ~ 0.39
-0.15 ~ -0.13
0.11 ~ 0.16
0.20 ~ 0.27
0.20 ~ 0.26
0.08 ~ 0.09
-0.09 ~ -0.08
1.56 ~ 2.02
Effect analysis of QTLs relative higher phenotype variance explanation for the silique related traits
6.75 ± 0.77
6.94 ± 0.83
4.95 ± 0.44
5.05 ± 0.60
3.67 ± 0.07
3.60 ± 0.15
3.58 ± 0.10
3.52 ± 0.10
3.36 ± 0.09
3.31 ± 0.16
3.57 ± 0.07
3.73 ± 0.11
25.04 ± 22.05
21.66 ± 37.39
23.94 ± 31.95
19.68 ± 47.00
20.51 ± 34.25
18.89 ± 28.45
0.79 ± 0.03
0.68 ± 0.02
0.83 ± 0.04
0.71 ± 0.02
0.69 ± 0.03
0.59 ± 0.05
For SB, 20 identified QTLs were obtained with PVE ranging from 3.44 % to 15.06 % (Additional file 3). Finally, these 20 identified QTLs were integrated into 14 consensus QTLs, including each two consensus QTLs were detected in three and two microenvironments (Table 2, Fig. 2). QTLs cqSB-C6-1 and cqSB-C6-2 with PVE >10 % were considered as major QTLs, as they were repeatedly detected in three microenvironments (10DL, 11DL and 12DL). Among the identified QTLs, qSB-C6-3 had the highest PVE and had a relatively high LOD value of 9.61. To identify its effect, genotype data of the DH population for the nearest marker e4m5-260 on C6 was collected, and a similar result for SL was observed (Table 3).
For ST, 16 identified QTLs were obtained with PVE of 3.26 %–10.41 % (Additional file 3). The 16 QTLs were integrated into 13 consensus QTLs, and two QTLs were detected in two microenvironments (Table 2, Fig. 2). QTL qST-A9-1, with the highest PVE of 10.41 %, was closely linked with marker e8m13–248 on A9, and the phenotypic data of DH lines carrying ‘A’ alleles was significantly smaller than lines carrying ‘B’ alleles, suggesting that alleles responsible for increasing ST existed in the male parent ‘KenC-8’ (Table 3).
For SV, 16 identified QTLs with PVE of 3.33 %–10.39 % were obtained (Additional file 3). The 16 identified QTLs were finally integrated into 11 consensus QTLs, including one and three QTLs were detected in three and two microenvironments, respectively (Table 2, Fig. 2). QTL cqSV-C6-3 that integrated from qSV-C6-3 (PVE of 10.14 %), qSV-C6-4 (10.39 %) and qSV-C6-5 (8.26 %) was considered as a major QTL. Further investigation revealed that the phenotypic data of lines carrying ‘A’ and ‘B’ alleles of tightly linked marker Na12A02 on C6 (qSV-C6-4 with the highest PVE) showed significant differences. The average SV was significantly larger in lines carrying ‘A’ alleles than ‘B’ alleles, suggesting that the genes responsible for increasing SV exist in the female parent ‘N53-2’ (Table 3).
Eight QTLs for SPS with PVE within the range of 4.32 %–11.16 % were obtained (Additional file 3), and then were integrated into six consensus QTLs. QTL cqSPS-C6-3 was detected in three environments (Table 2, Fig. 2), which was integrated from qSPS-C6-3, qSPS-C6-4 and qSPS-C6-5 with PVE of 4.32 %, 11.16 % and 5.94 %, respectively. QTL qSPS-C6-4 (with the highest PVE) was closely linked with marker CB10632 on C6, and further analysis revealed that the related phenotypic data of lines carrying ‘A’ alleles had significantly more SPS, on average, than lines carrying ‘B’ alleles. This suggests that the genes increasing SPS existed in the male parent ‘N53-2’ (Table 3).
For SD, six QTLs were identified with PVE of 4.28 %–6.71 % (Additional file 3), none of these were considered as major QTL. The six identified QTLs were integrated into five consensus QTLs. Only cqSD-C6-2 was detected in two microenvironments (Table 2).
Forty-two QTLs for TSW were detected in eight experiments, and 26 consensus QTLs were obtained, and three and six consensus QTLs were detected in four and two microenvironments, respectively (Additional file 3) . The PVE of these QTLs were within the range of 3.15 %–19.62 %. QTLs cqSW-A7-2 and cqSW-C1-1 were two major QTLs with PVE >10 % in two different experiments, respectively.
Detection of pleiotropic unique QTLs for silique related traits
The list of 17 pleiotropic unique QTL for silique related traits in KN population
Consensus map construction and QTL comparison
Five mapping populations (SS, BE, QN, TN and HY) were selected for construction of consensus map and QTL comparison (Additional file 6). For exact comparisons among populations, QTLs collected in each population were used in the first round of meta-analysis to obtain consensus QTLs (Additional file 6). The KN DH map was treated as a reference map, and the markers and QTLs were projected from the five chosen populations onto the KN map. Finally, a consensus map containing 1225 markers was constructed, with the exceptions of A9, C4, C7 and C8 that lacked a common marker (Additional file 6 and 7). A total of 117 consensus QTLs, including 76 on the A genome and 41 on the C genome, were located from the chosen five populations (Additional file 6). Finally, 34 QTLs (26 in the A genome and eight in C genome), including 7, 1, 1 and 25 QTLs controlling SL, SD, SPS and TSW were aligned to the consensus map, respectively (Additional file 8). The number of aligned QTLs on each chromosome ranged from one (A6, A10 and C9) to 11 (A3).
Comparative genome analysis and candidate gene identification
On the basis of the distribution of 24 conserved chromosomal blocks of Arabidopsis described by Schranz et al. , 33 blocks and 127 islands in the KN DH map were identified with respect to the Arabidopsis genome according to 141 markers with known sequence information .
All of the QTLs with CIs containing homologous genes were separately compared to the physical genomic regions of B. rapa (A genome) and B. oleracea (C genome) (Additional file 9). Based on the map alignment between Arabidopsis and B. napus, a total of 88 orthologs of 37 genes that controlled these silique traits in Arabidopsis were mapped onto the synteny blocks and islands (Additional file 10), and 25 genes were located in the CI of 27 QTLs (Additional file 9 and 11). Among the 25 genes, several were shown to have a relationship with the seven silique traits, including EXS, GASA4, ANT, TTG2 and AP2. Especially on chromosome C6, gene ATAGP19 (BnaA07g24460D) and CUL3B (scaffold_213) were associated with unique QTLs controlling SV, SD, ST and SW. The ATAGP19 mutant could lead to fewer siliques, less seed production and several abnormalities in cell size, number, shape and packing . CUL3B is essential for normal embryogenesis . At the same time, sequences of these genes available for B. rapa and B. oleracea were also mapped onto related chromosomes (Additional file 10).
In the present study, the genetic basis for the six silique traits (SPS, SL, SB, ST, SD and SV) at QTL level was performed in the KN DH population containing 348 lines, possibly the largest population yet used to perform QTL mapping for silique traits. SL is known to be positively correlated with SPS [31, 32], and the same result was also obtained in the present study. No significant positive relationship was observed between SL and TSW, which differed from previous results [1, 2]. The most persuasive explanation for the phenomenon was that the population used in the present study was different from other published populations, and each mapping population might carry different alleles for silique related traits and represent its own genetic background. For example, Yang et al. (2012) indicated that SL was significantly correlated with SW, using a RIL derived from a cross between an EMS mutant with extremely long SL and SW, and an inbred line with regular SL and SW . By comparison, the KN population for SL and SW showed high phenotypic variation but the two parents showed no significant difference. To date, little research has been conducted on SB, ST and SV and their relationships with other yield-related factors. In this study, SL, SB and SV all showed significant positive relationships with SPS and SD (Table 1). The present results revealed that SPS was significant negatively correlated with TSW, as also observed by Zhang et al. (2011) . It is noteworthy that all the six silique traits in the present study showed highly and positively correlated with the level of seed yield (Table 1), suggesting that at least some of the underlying genetic determinism or same genes are shared across the formation and development of these traits. These silique traits colud be served as important traits for section of high yield cultivar for breeding purpose. These results suggest that coordinating the balance among these agronomic traits of B. napus was the premise of ensuring the high and stable yield production.
Fine mapping and map-based cloning are efficient ways to dissect the genetic bases of quantitative traits and to identify genes underlying these traits. QTLs with large PVE and can be detected in different environments are more suitable for fine mapping and map-based cloning . Integrating QTL results across different environments can determine which QTLs are less influenced by environmental factors, and can help to isolate environment-specific QTLs . In the present study, 15 stable consensus QTLs expressed in at least two microenvironments, and 45 environment-specific QTLs were obtained (Table 2). The high proportion of environment-specific QTLs suggested a large impact of the natural environment on silique-associated traits. Most of them had minor effects with low LOD scores. It is clear that use such minor QTLs for MAS or clone such QTLs via fine mapping is unlikely. However, these QTLs provide an important genetic resource for further research . The explanation for the existence of stable QTLs was that they could be responsible for main genetic effects with high LOD scores, as QTLs with major effects are more likely to be stable across multiple environments [33-35]. In the present study, three QTLs with PVE > 10 % controlling SB (cqSB-C6-1 and cqSB-C6-2) and SV (cqSV-C6-2) were identified in at least three microenvironments in DL (Table 2). In addition, cqSB-C6-1 and cqSB-C6-2 were both integrated from three identified QTLs detected in DL and with mean additive effects of 0.22 and 0.22, respectively (Additional file 3). These results showed that when alleles from ‘KenC-8’ existed in the two QTLs, the SB could increase about 0.8 mm in DL. These kinds of QTLs might be worthy of attention when doing MAS for developing varieties with special adaptability. In addition, unique QTL uqC6-4 would be the most interesting for breeding. The uqC6-4 could control as much as six traits except for SL, which were all with positive additive effect. Compared with the seed yield and oil content QTLs detected in the same population, uqC6-4 was co-localized with the major QTL for seed yield (cqSY-C6-2), which was repeatedly detected in three environments , and was also co-localized with the QTL (qOC-C6-4) for seed oil content consistently expressed in four environments . Among the two parents of KN DH population, ‘KenC-8’ was the parental of the variety ‘Zayou59’ released in 1996 in China, and ‘N53-2’ was a DH line with seed oil content >50 % . QTLs detected in KN population could be used in breeding, not only for developing high-yielding cultivars, but also for increasing seed oil content.
Although the genetic determinism of these silique traits remain largely unknown in rapeseed, correlation analysis with the same DH population may provide some insights. Most pairs of traits showed significant phenotypic correlations in this study, which was also reflected by the genome location and effects of QTL detected. Many QTL clusters were observed on chromosomes A2, A3, A5, A6, A7, A9, C6 and C9 (Table 4). With the overlapping CIs of QTLs controlling different traits, 17 unique QTLs that controlled at least two traits were obtained and the unique QTL (uqC6-4) could control as much as six traits (Fig. 3). Co-localizations between QTLs for some traits were expected such for SB/ST/TSW due to the correlation (Table 1). This were also in agreement with the pleiotropic unique QTLs, such as uqA9-3 and uqC6-4 control the three traits simultaneously, and uqA3-2 and uqA7-2 have effect on SB and TSW, and uqA3-4 and uqC6-2 effect SB and ST simultaneously, and uqC9-2 have effect on ST and SW (Additional file 5). Previous studies on silique traits also demonstrated the existence of pleiotropic QTLs in B. napus [1, 4]. In addition, the consistency of QTL loci on the chromosome for various traits may facilitate selection efficiency by selecting markers closely associated with these traits. Furthermore, if the direction of the additive effect of the QTLs were the same, selection would be easier and more effective, such as for QTLs clustered in the same region of C6 (Additional file 5).
In general, a consensus map is a powerful tool to validate conserved QTLs across populations. To compare QTLs detected in different populations, a consensus map containing 1225 markers was constructed (Additional file 6). To our knowledge, there are few reports concerning QTL analysis for SB, ST and SV, and the 14, 13 and 11 consensus QTLs in the present study might be new QTLs. Compared with 34 projected QTLs from other five maps, three QTLs for TSW were the same as those identified in previous studies, including cqSW-A3-1, cqSW-A3-3 and cqSW-A4. In addition, 11, 6 and 5 consensus QTLs for SL, SPS and SD observed in the present study were differed from those of the other five populations, and hence might considered as potential new QTLs. The detected QTLs for SL had very high PVE, reaching 53.4 % and 54.4 % on A9 [1, 36], while no QTL for SL was obtained on A9 in the present study, possibly due to only minor differences in SL between the two parents. Comparison among different populations would also facilitate an assessment of intraspecific variation for these QTLs.
The establishment of comparative mapping between Arabidopsis and B. napus is a powerful tool to identify candidate genes underlying QTLs . In this study, 25 candidate genes were assigned to the CIs of 27 QTLs (Additional file 11), and several of these genes were shown to have a relationship with the six silique traits. Studies have revealed that the genes ANT, ARF2 and TOC1 underlying QTLs exist on chromosomes A1–A3 [1, 18, 36], whereas ANT underlying QTLs aligned on A3 and C1 was observed in the present study. Fan et al. (2010) detected QTLs for TSW on A7, which might correspond to TTG2 and MINI3 , whereas TTG2 underlying QTLs was aligned on A4 and A5 in the present results. As well as ANT, TTG2 and MINI3, gene EXS that controlled embryo development was aligned in the CIs of QTLs controlling SPS . GASA4, which positively affected both seed size and total seed yield, was located in the CI of a QTL controlling TSW . AP2 was found underlying QTLs controlling TSW on A2 and SL on C1, and ap2 mutants could produce larger seeds than wild type . On A5, IKU1 controlling seed size was located in the CI of QTLs controlling SPS and ST . Genes FIS2 and MEA were aligned in the CIs of QTLs controlling SPS and ST on A5 and QTLs controlling SL and ST on A6, respectively. These two genes repressed seed development in the absence of pollination . On C5, CUL3A was associated with QTLs controlling SL, which is essential for normal embryogenesis . ATAGP19 and CUL3B were associated with unique QTLs controlling SV, SD, ST and SW; and the AtAGP19 mutant had fewer siliques, less seed production and several abnormalities in cell size, number, shape and packing . Compared with currently published B. napus genome information, related B. napus genome information was also obtained. These findings provided useful resources on candidate genes for seed yield improvement.
We identified 82 individual QTLs for six silique traits (SPS, SL, SB, ST, SD and SV), and then integrated these QTLs into 60 consensus QTLs by meta-analysis. Among them, all of the 14, 13 and 11 consensus QTLs for SB, ST and SV might be the potential new QTLs because few QTLs for these traits were reported in B. napus. After compared with five published populations, 11, 6 and 5 potential new consensus QTLs for SL, SPS and SD were obtained. By in silico mapping analysis, we identified 25 candidate genes underlying 27 QTLs for silique related traits. Our results provided useful information into the regulatory model for the control of silique in B. napus.
Plant material and genetic linkage map
In this study, a segregating DH population, named KN DH, from a cross between male ‘KenC-8’ and female ‘N53-2’ was used to detect QTLs for silique traits. The KN DH population consisted of 348 lines, which was previously used for developing the linkage map (KN map) and QTL evaluation for oil content . The KN genetic linkage map was constructed with 403 molecular markers, including 275 simple sequence repeats, 117 sequence-related amplified polymorphisms, 10 sequence tagged sites and one intron fragment length polymorphism, which covered a total length of 1783.9 cM. The genotype of the population had two kinds of bands, the genotypes of the two parents: ‘KenC-8’ (B) and ‘N53-2’ (A). If the additive effect of QTLs was negative, this suggested the increasing allele originated from ‘KenC-8’, and otherwise it was from ‘N53-2’.
Field experiment and data collection
The field experiments were carried out in two locations: Wuhan in the Hubei Province (a semi winter-type rapeseed growing area, coded WH) and Dali in the Shannxi Province (a winter-type rapeseed growing area, coded DL). The seeds were sown in 2009–2011 in DL and 2012 in WH; and the related experiments were named 10DL, 11DL, 12DL and 13WH, respectively. The experiment location of WH was the experiment base of Huazhong University of Science and Technology, and DL was the experiment base of Hybrid Rapeseed Research Center of Shaanxi Province. No specific permissions were required for the field trials. Year–location combinations were treated as microenvironments, and then these microenvironments were divided into two contrasting macroenvironments: semi-winter and winter. For SL, SPS and SD, phenotypic data were collected in all of the four microexperiments. For SB, ST and SV, phenotypic data were collected in three microenvironments: 10DL, 11DL and 12DL. For TSW, phenotypic data was used to detect its relationships with the above six silique traits, and the QTLs for TSW obtained by Zhao et al. (2016) in eight microenvironments were used for detecting of unique QTLs and QTL comparison in the present study (Additional fine 3) . For each trial, all 348 lines together with their parents were grown in a completely randomized block design with two replications in WH, and three replications in DL. Each line was grown in two rows per plot in all three locations, every row contained about 12 plants, and the distance between rows and between plants was 0.4 and 0.2 m, respectively. The field management followed normal agricultural practice.
All DH lines were open-pollinated in all environments to keep a full record of fertilization rates. At the mature stage, five representational plants in the middle of each plot were selected to analyze the six silique traits. Then, six well-developed siliques, randomly selected from the first branch adjacent to the main inflorescence, were used to investigate phenotypic data. As the rapeseed silique is oval-shaped, the maximum width of a silique in the middle position was treated as “breadth”, and the minimum-data at the same position was treated as “thickness”. Phenotypic data were collected as follows. (1) SL (cm): estimated from the average length (excluding the length of beak) of the six well-developed chosen siliques. (2) SB and ST (both mm): estimated from the average breadth and thickness, respectively, of the six chosen siliques using a Vernier caliper. (3) SPS: estimated from the average value of six chosen siliques. (4) SV (mL): estimated from the average volume of the six chosen siliques using a cylinder containing water, from the volume change following submersion of siliques . (5) SD: measured according to SD = SPS/SL. All phenotype data were collected in three measurements per replication.
QTL analysis and meta-analysis for the silique traits
Basic statistical analysis was implemented using SPSS 18.0 software (SPSS Inc., Chicago, IL, USA). QTL analysis for SL, SB, ST, SPS, SV and SD were performed by the QTL Cartographer 2.5 software with composite interval mapping (CIM) . A walking speed was set to 2 cM and a window size of 10 cM with five background cofactors. The LOD threshold was determined by 1000-permutation test based upon a 5 % error rate , and a LOD of 2.7–3.2 was used to identify the existence of QTLs in each environment, and these QTLs were termed ‘identified QTLs’. The method for QTL nomenclature was as described by Wang et al. (2013) . A number was added if more than one QTL was located on the same chromosome, for instance, qSL-A5-2.
Identified QTLs repeatedly detected in different microenvironments and with overlapping CIs for the same trait were integrated into a consensus QTL by meta-analysis using BioMercator 2.1 software with default parameters. [46, 47], and were then named with initial letters ‘cq’, for instance, cqSL-A5-2. Consensus QTLs occurring in at least one microenvironment with PVE > 20 %, or in at least two microenvironments with PVE > 10 %, were considered as major QTLs . The consensus QTLs with overlapping CIs for different traits were further integrated into one unique QTL and designated with initial letters ‘uq’ (for example, uq-A5-1). If a consensus QTL had no overlapping CI with others, it was also regarded as a unique QTL.
Construction of the consensus map and QTL comparison
Consensus map construction and QTL comparison were carried out with BioMercator 2.1 [46, 47]. The populations chosen for comparison included three DH and two RIL populations: SS , BE , QN , TN  and HY  (Additional file 6). A ‘two-round’ strategy was used for QTL comparison . In the first round, QTLs identified in each population were collected, and those QTLs in one population of the same trait with overlapping CIs were integrated into consensus QTLs using meta-analysis. The consensus QTLs were named with the population abbreviation followed with ‘cq’ and the name of traits, a hyphen (-) and the linkage group. A number was added if more than one QTL was located on same chromosome (e.g. QN-cqSL-A10-1). In the second round, based on the common markers, the markers in homologous chromosomes were projected from other maps onto the reference KN DH map to construct a consensus map. Then the consensus QTLs detected in different populations were aligned to the KN map and QTLs were compared between published QTLs and our results.
Map alignment of B. napus with Arabidopsis, B. rapa and B. oleracea
Among the 403 markers mapped in the KN genetic map, 141 markers with known sequence information were treated as anchor markers for map alignment between B. napus and Arabidopsis . Finally, 33 blocks and 127 islands were aligned between Arabidopsis pseudochromosomes and the consensus KN DH linkage map . The related genes that may control the seven silique traits of Arabidopsis on each block were located according to their physical position in the Arabidopsis genome. Their genetic position on the KN DH map was obtained according to their closest anchor markers in the same synteny block. If the aligned gene(s) was located in the CI of a QTL, then the orthologous candidate genes were assumed to be associated with the target QTL. When compared with B. rapa and B. oleracea, the related genes consisting of each chromosome were obtained based on sequence information of molecular markers. To identify and locate the putative genes in B. napus, the homologous sequences of genes related to silique traits in the databases of B. rapa and B. oleracea were searched with the BLASTn program (E value < 1E–20 when using Arabidopsis genes related to silique traits as queries). The resulting sequences from the search were first mapped onto the chromosomes of B. rapa or B. oleracea and then placed on the B. napus linkage groups, based on the homologous collinear relationships between B. napus and B. rapa/B. oleracea via the Arabidopsis genome.
Availability of data and material
The datasets supporting the conclusions of this article are included within the article and its additional files.
The work was supported by the National Basic Research Program of China (2015CB150205), International Cooperation in Science and Technology Projects (2014DFA32210), Jiangsu Collaborative Innovation Center for Modern Crop Production and New Century Talents Support Program of the Ministry of Education of China (NCET110172).
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Yang P, Shu C, Chen L, Xu J, Wu J, Liu K. Identification of a major QTL for silique length and seed weight in oilseed rape (Brassica napus L.). Theor Appl Genet. 2012;125(2):285–96.View ArticlePubMedGoogle Scholar
- Oezer H, Oral E, Dogru U. Relationships between yield and yield components on currently improved spring rapeseed cultivars. Turk J Agric For. 1999;23:603–8.Google Scholar
- Zhang L, Yang G, Liu P, Hong D, Li S, He Q. Genetic and correlation analysis of silique-traits in Brassica napus L. by quantitative trait locus mapping. Theor Appl Genet. 2011;122(1):21–31.View ArticlePubMedGoogle Scholar
- 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(3):851–61.View ArticlePubMedPubMed CentralGoogle Scholar
- Ding G, Zhao Z, Liao Y, Hu Y, Shi L, Long Y, Xu F: Quantitative trait loci for seed yield and yield-related traits, and their responses to reduced phosphorus supply in Brassica napus. Ann Bot-London. 2012;109:747–59.View ArticleGoogle Scholar
- Chen W, Zhang Y, Liu X, Chen B, Tu J, Fu T. Detection of QTL for six yield-related traits in oilseed rape (Brassica napus) using DH and immortalized F2 populations. Theor Appl Genet. 2007;115(6):849–58.View ArticlePubMedGoogle Scholar
- Liu J, Hua W, Hu Z, Yang H, Zhang L, Li R, Deng L, Sun X, Wang X, Wang H: Natural variation in ARF18 gene simultaneously affects seed weight and silique length in polyploid rapeseed. Proc Natl Acad Sci U S A. 2015;112(37):5123–32.View ArticleGoogle Scholar
- Li S, Chen L, Zhang L, Li X, Liu Y, Wu Z, Dong F, Wan L, Liu K, Hong D et al: BnaC9.SMG7b functions as a positive regulator of number of seeds per silique in rapeseed (Brassica napus L.) by regulating the formation of functional female gametophytes. Plant Physiol. 2015;169(4):2744–60.PubMedGoogle Scholar
- Adamski NM, Anastasiou E, Eriksson S, O'Neill CM, Lenhard M. Local maternal control of seed size by KLUH/CYP78A5-dependent growth signaling. Proc Natl Acad Sci U S A. 2009;106(47):20115–20.View ArticlePubMedPubMed CentralGoogle Scholar
- Raman H, Raman R, Kilian A, Detering F, Long Y, Edwards D, Parkin IA, Sharpe AG, Nelson MN, Larkan N et al: A consensus map of rapeseed (Brassica napus L.) based on diversity array technology markers: applications in genetic dissection of qualitative and quantitative traits. BMC Genomics. 2013;14:277.View ArticlePubMedPubMed CentralGoogle Scholar
- Wang X, Wang H, Long Y, Li D, Yin Y, Tian J, Chen L, Liu L, Zhao W, Zhao Y et al: Identification of QTLs associated with oil content in a high-oil Brassica napus cultivar and construction of a high-density consensus map for QTLs comparison in B. napus. PLoS One. 2013;8(12):e80569.View ArticlePubMedPubMed CentralGoogle Scholar
- Scoles G, Fu F, Liu L, Chai Y, Chen L, Yang T, Jin M, Ma A, Yan X, Zhang Z et al: Localization of QTLs for seed color using recombinant inbred lines of Brassica napus in different environments. Genome. 2007;50(9):840–54.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- Long Y, Shi J, Qiu D, Li R, Zhang C, Wang J, Hou J, Zhao J, Shi L, Park BS et al: Flowering time quantitative trait loci analysis of oilseed Brassica in multiple environments and genomewide alignment with Arabidopsis. Genetics. 2007;177(4):2433–44.PubMedPubMed CentralGoogle Scholar
- Delourme R, Falentin C, Huteau V, Clouet V, Horvais R, Gandon B, Specel S, Hanneton L, Dheu JE, Deschamps M et al: Genetic control of oil content in oilseed rape (Brassica napus L.). Theor Appl Genet. 2006;113(7):1331–45.View ArticlePubMedGoogle Scholar
- Shi J, Li R, Zou J, Long Y, Meng J. A dynamic and complex network regulates the heterosis of yield-correlated traits in rapeseed (Brassica napus L.). PLoS One. 2011;6(7):e21645.View ArticlePubMedPubMed CentralGoogle Scholar
- Wang X, Wang H, Long Y, Liu L, Zhao Y, Tian J, Zhao W, Li B, Chen L, Chao H et al: Dynamic and comparative QTL analysis for plant height in different developmental stages of Brassica napus L. Theor Appl Genet. 2015;128(6):1175–92.View ArticlePubMedGoogle Scholar
- Zhou Q, Fu D, Mason AS, Zeng Y, Zhao C, Huang Y. In silico integration of quantitative trait loci for seed yield and yield-related traits in Brassica napus. Mol Breeding. 2014;33(4):881–94.View ArticleGoogle Scholar
- Ding G, Liao Y, Yang M, Zhao Z, Shi L, Xu F. Development of gene-based markers from functional Arabidopsis thaliana genes involved in phosphorus homeostasis and mapping in Brassica napus. Euphytica. 2011;181(3):305–22.View ArticleGoogle Scholar
- Raman H, Raman R, Coombes N, Song J, Prangnell R, Bandaranayake C, et al. Genome-wide association analyses reveal complex genetic architecture underlying natural variation for flowering time in canola. Cell & Environment: Plant; 2015. doi:10.1111/pce.12644.Google Scholar
- Cai G, Yang Q, Yang Q, Zhao Z, Chen H, Wu J, Fan C, Zhou Y: Identification of candidate genes of QTLs for seed weight in Brassica napus through comparative mapping among Arabidopsis and Brassica species. BMC Genet. 2012;13:105.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhang L, Li S, Chen L, Yang G. Identification and mapping of a major dominant quantitative trait locus controlling seeds per silique as a single Mendelian factor in Brassica napus L. Theor Appl Genet. 2012;125(4):695–705.View ArticlePubMedGoogle Scholar
- Chalhoub B, Denoeud F, Liu S, Parkin IA, Tang H, Wang X, Chiquet J, Belcram H, Tong C, Samans B et al: Early allopolyploid evolution in the post-Neolithic Brassica napus oilseed genome. Science. 2014;950–953.Google Scholar
- Liu S, Liu Y, Yang X, Tong C, Edwards D, Parkin IA, Zhao M, Ma J, Yu J, Huang S et al: The Brassica oleracea genome reveals the asymmetrical evolution of polyploid genomes. Nat Commun. 2014;5:3930.PubMedPubMed CentralGoogle Scholar
- Wang X, Wang H, Wang J, Sun R, Wu J, Liu S, Bai Y, Mun JH, Bancroft I, Cheng F et al: The genome of the mesopolyploid crop species Brassica rapa. Nat Genet. 2011;43(10):1035–9.View ArticlePubMedGoogle Scholar
- Jiang C, Zeng Z. Multiple trait analysis of genetic mapping for quantitative trait loci. Genetics. 1995;140(3):1111–27.PubMedPubMed CentralGoogle Scholar
- Zhao W, Wang X, Wang H, Tian J, Li B, Chen L, Chao H, Long Y, Gan J, Liang W et al: Genome-wide identification of QTL for seed yield and yield-related traits and construction of a high-density consensus map for QTL comparison in Brassica napus. Front Plant Sci. 2016;7:17.PubMedPubMed CentralGoogle Scholar
- 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–42.View ArticlePubMedGoogle Scholar
- Yang J, Sardar HS, McGovern KR, Zhang Y, Showalter AM. A lysine-rich arabinogalactan protein in Arabidopsis is essential for plant growth and development, including cell division and expansion. Plant J. 2007;49(4):629–40.View ArticlePubMedGoogle Scholar
- Thomann A, Brukhin V, Dieterle M, Gheyeselinck J, Vantard M, Grossniklaus U, et al. Arabidopsis CUL3A and CUL3B genes are essential for normal embryogenesis. Plant J. 2005;43(3):437–48.View ArticlePubMedGoogle Scholar
- Aytaç Z, Kınacı G. Genetic variability and association studies of some quantitative characters in winter rapeseed (Brassica napus L). Afr J Biotechnol. 2009;8:3547–54.Google Scholar
- Chay P, Thurling N. Identification of genes controlling pod length in spring rapeseed, Brassica napus L., and their utilization for yield improvement. Plant Breeding. 1989;103(1):54–62.View ArticleGoogle Scholar
- Zhuang J, Lin H, Lu J, Qian H, Hittalmani S, Huang N, Zheng K: Analysis of QTL × environment interaction for yield components and plant height in rice. Theor Appl Genet. 1997;95:799–808.View ArticleGoogle Scholar
- Tanksley SD. Mapping polygenes. Annu Rev Genet. 1993;27:205–33.View ArticlePubMedGoogle Scholar
- Paterson AH, Damon S, Hewitt JD, Zamir D, Rabinowitch HD, Lincoln SE, Lander ES, Tanksley SD: Mendelian factors underlying quantitative traits in tomato: comparison across species, generations, and environments. Genetics. 1991;127(1):181–97.PubMedPubMed CentralGoogle Scholar
- Li N, Shi J, Wang X, Liu G, Wang H. A combined linkage and regional association mapping validation and fine mapping of two major pleiotropic QTLs for seed weight and silique length in rapeseed (Brassica napus L.). BMC Plant Biol. 2014;14:114.View ArticlePubMedPubMed CentralGoogle Scholar
- Fan C, Cai G, Qin J, Li Q, Yang M, Wu J, Fu T, Liu K, Zhou Y: Mapping of quantitative trait loci and development of allele-specific markers for seed weight in Brassica napus. Theor Appl Genet. 2010;121(7):1289–301.View ArticlePubMedGoogle Scholar
- Canales C, Bhatt AM, Scott R, Dickinson H. EXS, a putative LRR receptor kinase, regulates male germline cell number and tapetal identity and promotes seed development in Arabidopsis. Curr Biol. 2002;12(20):1718–27.View ArticlePubMedGoogle Scholar
- Roxrud I, Lid SE, Fletcher JC, Schmidt EDL, Opsahl-Sorteberg HG. GASA4, one of the 14-member Arabidopsis GASA family of small polypeptides, regulates flowering and seed development. Plant Cell Physiol. 2007;48(3):471–83.View ArticlePubMedGoogle Scholar
- Ohto M, Floyd SK, Fischer RL, Goldberg RB, Harada JJ. Effects of APETALA2 on embryo, endosperm, and seed coat development determine seed size in Arabidopsis. Sex Plant Reprod. 2009;22(4):277–89.View ArticlePubMedPubMed CentralGoogle Scholar
- Luo M, Dennis ES, Berger F, Peacock WJ, Chaudhury A. MINISEED3 (MINI3), a WRKY family gene, and HAIKU2 (IKU2), a leucine-rich repeat (LRR) KINASE gene, are regulators of seed size in Arabidopsis. Proc Natl Acad Sci U S A. 2005;102(48):17531–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Luo M, Bilodeau P, Dennis ES, Peacock WJ, Chaudhury A. Expression and parent-of-origin effects for FIS2, MEA, and FIE in the endosperm and embryo of developing Arabidopsis seeds. Proc Natl Acad Sci U S A. 2000;97(19):10637–42.View ArticlePubMedPubMed CentralGoogle Scholar
- Escribano MR, Santalla M, De Ron AM. Genetic diversity in pod and seed quality traits of common bean populations from northwestern Spain. Euphytica. 1997;93(1):71–81.View ArticleGoogle Scholar
- Wang S, Basten CJ, Zeng ZB. Windows QTL cartographer 2.5. Department of statistics. Raleigh, NC: North Carolina State University; 2012.Google Scholar
- Doerge RW, Churchill GA. Permutation tests for multiple loci affecting a quantitative character. Genetics. 1996;142(1):285–94.PubMedPubMed CentralGoogle Scholar
- 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–6.View ArticlePubMedGoogle Scholar
- Goffinet B, Gerber S. Quantitative trait loci: a meta-analysis. Genetics. 2000;155(1):463–73.PubMedPubMed CentralGoogle Scholar