Original Research ARTICLE
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
- 1Hybrid Rapeseed Research Center of Shaanxi Province, Shaanxi Rapeseed Branch of National Centre for Oil Crops Genetic Improvement, Yangling, China
- 2Department of Biotechnology, College of Life Science and Technology, Huazhong University of Science and Technology, Wuhan, China
- 3Key Laboratory of Cotton and Rapeseed, Ministry of Agriculture, Institute of Industrial Crops, Jiangsu Academy of Agricultural Sciences, Nanjing, China
- 4Institute of Biotechnology, Chinese Academy of Agricultural Sciences, Beijing, China
- 5Hubei Collaborative Innovation Center for the Characteristic Resources Exploitation of Dabie Mountains, Huanggang Normal University, Huanggang, China
- 6Department of Applied Biological Science, College of Agriculture and Biotechnology, Zhejiang University, Hangzhou, China
Seed yield (SY) is the most important trait in rapeseed, is determined by multiple seed yield-related traits (SYRTs) and is also easily subject to environmental influence. Many quantitative trait loci (QTLs) for SY and SYRTs have been reported in Brassica napus; however, no studies have focused on seven agronomic traits simultaneously affecting SY. Genome-wide QTL analysis for SY and seven SYRTs in eight environments was conducted in a doubled haploid population containing 348 lines. Totally, 18 and 208 QTLs for SY and SYRTs were observed, respectively, and then these QTLs were integrated into 144 consensus QTLs using a meta-analysis. Three major QTLs for SY were observed, including cqSY-C6-2 and cqSY-C6-3 that were expressed stably in winter cultivation area for 3 years and cqSY-A2-2 only expressed in spring rapeseed area. Trait-by-trait meta-analysis revealed that the 144 consensus QTLs were integrated into 72 pleiotropic unique QTLs. Among them, all the unique QTLs affected SY, except for uq.A6-1, including uq.A2-3, uq.C1-2, uq.C1-3, uq.C6-1, uq.C6-5, and uq.C6-6 could also affect more than two SYRTs. According to the constructed high-density consensus map and QTL comparison from literatures, 36 QTLs from five populations were co-localized with QTLs identified in this study. In addition, 13 orthologous genes were observed, including five each gene for SY and thousand seed weight, and one gene each for biomass yield, branch height, and plant height. The genomic information of these QTLs will be valuable in hybrid cultivar breeding and in analyzing QTL expression in different environments.
Brassica napus (AACC, 2n = 38) originated from hybridization between Brassica rapa (AA, 2n = 20) and Brassica oleracea (CC, 2n = 18; UN, 1935), and is the second most important oilseed crop after soybean (Basunanda et al., 2010). As the global requirements for rapeseed oil and protein are growing rapidly, increasing seed yield (SY) is the main breeding aim at present. SY is directly determined by yield component traits, including thousand seed weight (SW), pod number per plant and seed number per pod (Qzer et al., 1999; Quarrie et al., 2006). In addition, SY is also indirectly influenced by other seed yield related traits (SYRTs), such as biomass yield (BY), plant height (PH), first effective branch height (BH), first effective branch number (FBN), length of main inflorescence (LMI), and pod number of main inflorescence (PMI) in B. napus (Qiu et al., 2006; Li et al., 2007; Shi et al., 2009). Interactions between SY, SW, PH, BH, FBN, LMI, and PMI were observed in previous studies (Yu, 1998; Zhang et al., 2006).
SY and SYRTs are all complex quantitative traits controlled by multiple genes (Kearsey and Pooni, 1998). QTL analysis has proved a powerful genetic approach to dissect complex traits (Paran and Zamir, 2003). Many QTLs for SY and SYRTs have been reported in B. napus, such as QTLs for SY being mainly located on A10, C3, and C6 (Quijada et al., 2006; Udall et al., 2006; Maccaferri et al., 2008). In addition, studies related to QTLs for SY and/or several SYRTs have also been performed (Chen et al., 2007, 2010; Li et al., 2007; Maccaferri et al., 2008; Shi et al., 2009; Basunanda et al., 2010; Ding et al., 2012; Cai et al., 2014). As the genetic backgrounds of different mapping populations for B. napus vary considerably, the number and location of QTLs detected in different populations also differ, thus is very necessary to contrast the QTLs for SY and SYRTs and select the common QTLs in different populations. Although many QTLs for SY and SYRTs have been reported, studies that simultaneously focused on the eight agronomic traits (SY, BY, SW, PH, BH, FBN, LMI, and PMI) are rare. Moreover, the candidate genes for these QTLs have rarely been mentioned. Comparative mapping among the model plant Arabidopsis thaliana with related species is a powerful tool to identify candidate genes. For example, Long et al. (2007) obtained the candidate gene BnFLC10 underlying QTL qFT10-4 and identified the key gene controlling differentiation of winter or spring type rapeseed based on comparative mapping analysis. Shi et al. (2009) and Ding et al. (2011) also obtained the candidate genes controlling flower time and seed phosphorus concentration, respectively, by comparative mapping with the Arabidopsis genome. Comparative mapping among B. napus, Arabidopsis, B. rapa, and B. oleracea genomes is necessary to obtain candidate genes in the confidence intervals (CIs) of QTLs for SY and SYRTs.
In order to increase statistical power and precision of obtaining QTLs, a high-density genetic linkage map is considered as a key factor (Jiang and Zeng, 1995). Several high-density genetic maps for B. napus have been constructed by integrating different linkage maps based on common molecular markers from different populations (Lombard and Delourme, 2001; Scoles et al., 2007; Raman et al., 2013). For example, Lombard and Delourme (2001) constructed a consensus map covering a total length of 2429.0 cM by integrating three individual linkage maps, and Wang et al. (2013) constructed a high-density consensus map with 1335 markers covering 2395.2 cM of the total genome length by merging eight individual linkage maps from different populations. Zhou et al. (2014) used 15 published articles concerning B. napus mapping experiments over the last decade and carried out in silico integration of 1960 QTLs with 13 SY and SYRTs, a total of 736 QTLs were mapped onto 283 loci in the A and C genomes of B. napus.
In the present study, a large doubled haploid (DH) population containing 348 lines was used to investigate the QTLs for SY and SYRTs in multiple environments, and then a consensus map was constructed for QTL comparison between the KN (the population used in this study) and five other published populations. These results provide abundant useful information to further understanding of the genetic mechanisms of SY and SYRTs, and could be used in marker assisted selection for improving SY in B. napus.
Materials and Methods
Plant Material and Field Experiments
A DH population, named KN and containing 348 lines derived from KenC-8 and N53-2, was used in this study (Wang et al., 2013). 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 KN population and its parents were grown in eight environments, including a winter rapeseed area, Dali of Shannxi Province (coded DL), in northwest China for five successive years (September–May of 2008–2009, 2009–2010, 2010–2011, 2011–2012, and 2012–2013); a spring rapeseed area, Sunan of Gansu Province (coded GS), in northwest China for three successive years (April–September of 2010, 2011, and 2012). Year-location combinations were treated as micro-environments, for example, 09DL means that the experiment was carried out in September–May of 2009–2010 at DL. Meanwhile, each year-location combination was treated as a trial.
The field experiments followed a randomized complete block design. The KN population, together with the two parents, was planted in DL and GS with three and two replications, respectively. Each field trial consisted of 348 lines. Each line was grown in a two-row plot with 40 cm between rows and 20 cm between individuals, and row length of 250 cm in all environments.
Measurement of Phenotypic Data of SY and SYRTs
Phenotypic data for SY (g/plant) were recorded with five representative plants in the middle of each plot. These five plants were also used for measurement of other SYRTs: BY (g/plant), SW (g), PH (cm), BH (cm), FBN, LMI (cm), and PMI. Because of a strong requirement for vernalization, the N53-2 and some DH lines did not flower or fully mature in the spring area (10GS, 11GS, and 12GS), and so SY of these DH lines were treated as missing data. SY was the average dry weight of seeds of the five representative individuals. BY was measured as the average total above-ground dry weight of the five plants (excluding the seeds). SW was the average dry weight of 1000 well-filled seeds from the five samples. PH was the average height of the five individuals, measured from the base of the stem to the tip of the main inflorescence. FBN was the number of branches arising from the main stem of each harvested individual. LMI was measured from the bottom to the top of the main inflorescence. PMI was effective pod number from the bottom to the top of the main inflorescence.
Statistical Analysis, QTL Mapping and Meta-Analysis
Estimates of means and variances for the SY and SYRTs were implemented using SPSS 18.0 software (SPSS Inc., Chicago, IL, USA). The QTL information for PH was according to Wang et al. (2015) in 10DL, 11DL, 12DL, 11GS, 12GS, and 12WH. QTL detection for other traits was conducted by composite interval mapping with Windows QTL Cartographer 2.5 software (Wang et al., 2007). The estimated additive effect and phenotypic variation explained (PVE) by each putative QTL were obtained using composite interval mapping model. Significance levels for the log odds score methods (LOD) were determined by 1000-permutation test corresponding to P = 0.05, and LOD of 2.8–3.1 was used to, respectively, identify significant QTLs in each environment, and these QTLs were termed “identified QTL.” QTLs that mapped to the same region with overlapping CIs were assumed to be the same, and BioMercator 2.1 software was used to integrate these QTLs into consensus QTLs using the meta-analysis method (Arcade et al., 2004). If a consensus QTL had at least one environment with PVE ≥ 20% or at least two environments with PVE ≥ 10%, the QTL was defined as a major QTL; the remaining QTLs were defined as minor QTLs (Maccaferri et al., 2008).
The identified QTLs for SY and SYRTs were named according to Wang et al. (2013); for example, the QTL abbreviation “qSY” (q, QTL; SY, seed yield) suffixed with the linkage group (A1–A10, C1–C9), a hyphen (-), and finally the serial number of QTLs in the linkage group (e.g., qSY-A2-1). The QTL integrations were adopted by meta-analysis; for example, the identified QTLs were integrated into consensus QTLs trait-by-trait, and the consensus QTLs for SY and SYRTs with overlapping CIs were integrated into pleiotropic unique QTLs using BioMercator 2.1 (Arcade et al., 2004). The name of consensus QTLs and unique QTLs referred to the name of identified QTLs. For each unique QTL, one or more consensus QTLs for SYRTs were chosen as indicator QTLs, which were defined as potential genetic determinants of the co-localized QTL for SY.
QTL Projection from Other Populations onto the KN Map, and QTL Comparison among the Different Populations
The map projection package of BioMercator 2.1 software was used for QTL projection of SY and SYRTs in five previously reported populations onto the KN genetic map (Table 1), including the QN (Quantum × No. 2127-17; Chen et al., 2007), SE (SI-1300 × Eagle; Li et al., 2007), ER (Express617 × R53; Radoev et al., 2008), TN (Tapidor × Ningyou7; Shi et al., 2009), and BE (B104-2 × Eyou Changjia; Ding et al., 2012) populations. The method for projection of QTLs from different linkage groups was according to Arcade et al. (2004). The method for QTL comparison from different linkage groups was a “two-round” strategy (Shi et al., 2013), and the detailed methods for QTL comparison from different populations were found in Ding et al. (2012) and Jiang et al. (2014). The consensus QTLs from different populations were named with the population abbreviation followed by the consensus QTL names for QTL comparison (e.g., KNcqSY-A2-1).
Candidate Gene Observations by Comparative Mapping among Arabidopsis, B. rapa, B. oleracea, and B. napus
Among the 403 markers mapped in the KN genetic map (Table S1), 141 markers with known sequence information were used for the sequence comparisons of the Arabidopsis genome database with other Brassica species (http://www.arabidopsis.org/). The databases of B. oleracea (Liu et al., 2014), B. rapa (http://brassicadb.org/brad/), and B. napus (http://www.genoscope.fr/brassicanapus/) were used for confirmation of homologous genes on the genomes of Arabidopsis, B. rapa, B. oleracea, and B. napus. Firstly, the 141 markers with known sequence information were used as anchored markers to carry out map alignment between B. napus and Arabidopsis according to the method of Long et al. (2007). If three or more sequence informative markers in the KN population were closely linked within one conserved block of Arabidopsis (Schranz et al., 2006), a synteny block was considered to exist. If there were only one or two sequence informative marker(s), this was recognized as an insertion segment. Secondly, if a synteny block or insertion segment were co-localized with the CI of a QTL, the genes underlying the synteny block or insertion segment were considered as candidate genes for the QTL. Thirdly, the genes of Arabidopsis were used to identify homologous genes in B. rapa, B. oleracea, and B. napus. The detailed methods are found in Long et al. (2007) and Shi et al. (2009).
The authors declare that the experiments comply with the current laws of the country in which they were performed.
Phenotypic Analysis and Genetic Correlation Between SY and SYRTs
The SY and SYRTs of the two parents and the KN population showed differences in most micro-environments (Table 2). There was a wide segregation range of SY, with a continuous normal distribution and transgressive segregation in all trials (Figure 1), suggesting that SY was a quantitative trait with polygenic control. Seven other SYRTs (BY, SW, PH, BH, FBN, LMI, and PMI) also showed a wide segregation range in all trials with normal or near-normal distributions.
Table 2. Mean values and phenotypic variation for SY and SYRTs of KN population evaluated in different micro-environments.
Figure 1. The frequency distribution of SY and SYRTs in multiple environments. The units of the x-axis are the phenotypic values, and the units of the y-axis are the number of lines. SY and SYRTs in different experiments was discriminated using different colored boxes. A unit of measurement: seed yield (g), biomass yield (g), thousand seed weight (g), plant height (cm), first effective branch height (cm), first effective branch number, length of main inflorescence (cm), pod number of main inflorescence.
The correlations between SY and SYRTs showed large differences (Table 3). The results indicated that SY was highly and positively correlated with SYRTs except for FBN, and especially for BY with a correlation coefficient of 0.83. LMI was significantly positively correlated with SY and SYRTs except for FBN. SW was significantly positively correlated with SY (0.31), BY (0.53), PH (0.42), and LMI (0.37). The high correlations among SY and SYRTs indicated that these traits might be controlled by the same kinds of genes in some cases.
Genome-Wide QTL Detection for SY and SYRTs
In total, 226 identified QTLs were observed for SY and SYRTs: 18 for SY and 208 for SYRTs (Table S2). These 226 QTLs were integrated into 144 consensus QTLs, which were located on 18 linkage groups with the exception of A8 (Figure 2, Figure S1).
Figure 2. Genetic linkage map and the location of QTLs for SY and SYRTs in the KN linkage map. The 144 consensus QTLs for SY and SYRTs were distributed on 18 linkage groups with the exception of A8, A1–A10 were represented by the A genome and C1–C9 were represented by the C genome in B. napus. The loci names were listed on the right of the linkage groups, while position of loci were showed on the left side of linkage groups. The consensus QTLs associated with SY and SYRTs were indicated by bars with various backgrounds on the left of each linkage group (Red bar, seed yield; Yellow bar, biomass yield; Green bar, thousand seed weight; Blue bar, plant height; Purple bar, first effective branch height; Brown bar, first effective branch number; Black bar, length of main inflorescence; Dark blue bar, pod number of main inflorescence).
For SY, there were nine consensus QTLs obtained, mainly located on A2, A6, C1, and C6 (Figure 2, Table S2). Four of these QTLs were repeatedly detected in different experiments, including cqSY-C1-2 and cqSY-C6-2, detected in four successive years in DL (09DL, 10DL, 11DL, and 12DL), and cqSY-C6-3 detected in three experiments (09DL, 10DL, and 11DL). In addition, both cqSY-C6-2 and cqSY-C6-3 were assumed to be major QTLs with PVE > 10% in two environments. Meanwhile, cqSY-A2-2, which was only observed in 10GS, was also a major QTL with PVE = 20.91% (Table 4).
For SYRTs, the number of QTLs for different traits clearly differed. For BY, 13 consensus QTLs were obtained, mainly located on A7, C1, and C6 with 4.34–19.96% of PVE (Figure 2, Table S2). Five of these QTLs were repeatedly detected in different experiments, for example, cqBY-C1-2, cqBY-C6-1, and cqBY-C6 were detected in three experiments.
For SW, 25 consensus QTLs were detected, distributed on 11 linkage groups (Figure 2, Table S2). Among them, nine consensus QTLs were repeatedly detected in different environments (Table S3). QTLs cqSW-A7-2, cqSW-C1-1, and cqSW-C1-2 were repeatedly detected in four experiments in winter area—with cqSW-A7-2 and cqSW-C1-1 regarded as two major QTLs with PVE > 10% in two environments (Table 4). The QTL cqSW-C9-1 that appeared in both winter and spring areas was an insensitive QTL in terms of response to environments.
For PH, 18 consensus QTLs were obtained at the mature stage in eight environments (Wang et al., 2015), and mainly located on A3, C6, and C9 (Figure 2, Table S2). Seven QTLs were repeatedly detected in different environments (Table S3), including cqPH-A3-3 detected in six environments, and cqPH-C6-2 and cqPH-C9-5 detected in three environments (Wang et al., 2015).
For BH, 27 consensus QTLs were integrated from 41 identified QTLs, and mainly located on A2, A3, A10, and C9 (Table S2). Ten QTLs were repeatedly detected in different environments (Table S3), for example, cqBH-A3-2 and cqBH-A10 were repeatedly observed both in winter and spring areas. Additionally, QTL cqBH-A2 was regarded as a major QTL with PVE > 10% in 11GS and 12GS. Because BH was highly positively correlated with PH, the major QTL cqBH-A2 for BH might also regulate PH.
For FBN, 25 consensus QTLs were obtained and were located on 12 chromosomes (Table S2). Three QTLs were repeatedly detected in both winter and spring areas: cqFBN-A3-1, cqFBN-A3-2, and cqFBN-C3-2. In addition, cqFBN-A2 and cqFBN-C2 were repeatedly detected in two and three experiments, respectively. One important QTL, cqFBN-C6-1 with PVE > 10% in two environments, was considered as a major QTL (Table 4).
For LMI, 24 identified QTLs were detected and integrated into 17 consensus QTLs, of which four consensus QTLs were repeatedly detected in different experiments (Table S2). For example, cqLMI-A3-1 and cqLMI-A3-2 were repeatedly detected in three and four experiments, respectively. However, no major QTL was observed for LMI. For PMI, 10 consensus QTLs were obtained and only two were repeatedly detected in different environments: cqPMI-A10 and cqPMI-C6 in spring (11GS and 12GS) and winter areas (09DL and 11DL), respectively. No QTLs for PMI reached the standard of a major QTL.
In conclusion, with 22, linkage group A3 had the largest number of consensus QTLs, followed by C6 and C9, both with 17. More than half of the consensus QTLs (95 of 144) for SY and SYRTs were detected in one micro-environment (Table S3). Otherwise, 28, 13, and 7 consensus QTLs were identified in two, three and four micro-environments, respectively (Figure 3). In total, 102 and 34 consensus QTLs were detected in winter (DL) and spring areas (GS), respectively, and only eight appeared in both areas (Figure 3, Table S3). These results indicated that the majority of consensus QTLs were expressed principally in response to a specific environment.
Figure 3. Expression response of 144 consensus QTLs in nature environments. (A) Number of consensus QTLs appeared in one to eight micro-environments. (B) Number of consensus QTLs appeared in winter, spring or both macro-environments.
The Unique QTL Analysis for SY and SYRTs
Of the 144 consensus QTLs for SY and SYRTs, 112 QTLs with overlapping CIs were integrated into 40 unique QTLs, and the remaining 32 consensus QTLs were only detected for one trait (Table S4). Altogether, 72 unique QTLs were obtained, in which 39 unique QTLs, respectively, affected 2–6 different traits. These unique QTLs were considered as pleiotropic (Table S4), such as uq.A2-3 which was integrated from two QTLs for SY and five for SYRTs. Notably, uq.A2-3, uq.C6-5, and uq.C6-6 were integrated from 3 to 6 QTLs of different traits and contained three major QTLs for SY (Table 5).
There were 21 unique QTLs observed, which were integrated from two consensus QTLs that controlled different traits (Table S4). For example, uq.A3-1, uq.A3-2, uq.A3-6, and uq.C5-2 all included QTLs for SW and BH. Likewise, uq.A2-1 and uq.A2-2 were composed of QTLs for SW and PMI. Although uq.A10-3 was a pleiotropic QTL and was integrated from three consensus QTLs, it was only closely linked to FBN and LMI. In addition, uq.A4-1 controlled BY and uq.C6-1 controlled SY, which were both closely linked to the QTL for SW. These results also explained why the QTL for SW was closely linked to the QTL for SY or BY.
Some unique QTLs were integrated from QTLs for more than two SYRTs (Table S5). QTL uq.A3-4 controlled BY, uq.A3-9, and uq.A3-11 controlled LMI, and uq.A7-2 controlled SW, which were all, respectively, tightly linked to the BH and PH. Uq.A7-3 controlled FBN and uq.C1-2 controlled SY, which were both respectively closely linked to BY, SW, and LMI. It is noteworthy that uq.A10-2 and uq.C1-3 were closely related to five SYRTs; uq.A10-2 controlled BH and FBN, and uq.C1-3 controlled SY and BY, which were both tightly linked to SW, PH, and LMI.
Consensus Map Construction and QTL Comparison for SY and SYRTs among Different Mapping Populations
In the present study, five published populations (QN, SE, ER, TN, and BE) for QTL analysis of SY and SYRTs were used for consensus map construction and QTL comparison (Table 1, Table S6). QTLs collected in each population were first integrated into consensus QTLs using BioMercator 2.1 software. A high-density consensus map with 907 molecular markers was constructed (Figure S3). A total of 480 consensus QTLs for SY and SYRTs from five populations were obtained and 166 consensus QTLs were successfully projected onto the consensus map (Figure 4, Figures S2, S3).
Figure 4. The consensus map and QTLs for SY and SYRTs in different populations (A6 and C3). Markers with blue color indicated these makers were projected from other maps on the KN map based on common markers by BioMercator 2.1 software. QTLs for SY and SYRTs detected in different populations were discriminated with different color bars on the left of each linkage group. Red bar, SY (seed yield); Orange bar, BY (biomass yield); Cambridge blue bar, SW (thousand seed weight); Purple bar, PH (plant height); Claybank bar, BH (first effective branch height); Green bar, FBN (first effective branch number); Breen bar, LMI (length of main inflorescence); Blue-green bar, PMI (pod number of main inflorescence). The KN population and five populations were indicated by disks with various backgrounds on the bars of each QTL. Red disk, KN population; Blue disk, TN population; Light green disk, SE population; Light purple disk, ER population; Brown disk, BE population; Yellow disk, QN population.
A total of 34 QTLs for SY were projected onto the KN consensus map, half of which were located on A2, A5, A6, C2, and C6 (Figure 5, Figure S3, Table S7). This revealed that KNcqSY-A2-2, KNcqSY-C6-1, and five QTLs of the TN population (TNqSY-A2-2, TNqSY-A2-3, TNqSY-A2-4, TNqSY-C6-1, and TNqSY-C6-3) were co-localized on A2 and C6, respectively. The results indicated that these QTLs might have some important genes for SY and could be expressed stably in different genetic backgrounds. In addition, five ortholougs genes (GASA4, ATCLH1, RBCS1A, LQY1, and ATGGH1) for SY were aligned in the CIs of these QTLs (Table 6). The two QTLs (KNcqSY-C1-1 and KNcqSY-C1-2) on C1 for SY were not observed in other genetic linkage groups, and these may be two new QTLs in the KN population.
Figure 5. The number and distribution of projection QTLs for SY and SYRTs from five population in B. napus. The x-axis of abscissas indicated 19 linkage groups, y-axis indicated the number of projection QTLs and z-axis indicated SY and SYRTs.
Table 6. Positions of homologous genes of Arabidopsis mapped on the KN linkage map and orthologous genes between B. napus and B. rapa/B. oleracea.
For SW, 40 QTLs were projected onto the KN consensus map (Figure 5, Figure S3, Table S7). KNcqSW-A2-4 and TNqSW-A2-3 were co-localized on A2. KNcqSW-A3-3 and five QTLs of different populations were co-localized on A3, including two QTLs of TN (TNqSW-A3-8 and TNqSW-A3-9) and three of BE (BEcqSW-A3-2, BEcqSW-A3-3, and BEcqSW-A3-4). The co-localized QTLs on KN and TN (KNcqSW-C6-2, TNqSW-C6-3, KNcqSW-C6-3, and TNqSW-C6-2) were also observed on the C6 linkage group (Table S8). A total of five orthologous genes for SW were observed: PTH2, AP2, LCR64, LCR65, and PDF1 (Table 6). These results indicated that the QTLs for SW were reliable and reproducible.
For PH, 36 QTLs were projected onto the KN consensus map, including eight and nine onto A2 and A3, respectively (Figure 5, Figure S3, Table S7). KNcqPH-A2 and two QTLs of TN (TNqPH-A2-1 and TNqPH-A2-2) were co-localized on A2. KNcqPH-A3-2 and three QTLs (TNqPH-A3-1, BEqPH-A3-4, and BEqPH-A3-3) were co-localized on A3. Likewise, KNcqPH-A3-3 and TNqPH-A3-4, KNcqPH-C3 and QNqHPH-C3 were also co-localized on A3 and C3, respectively. In addition, KNcqPH-C6-1 and TNqPH-C6-3, KNcqPH-C6-2 and TNqPH-C6-2 were co-localized in the same CI of C6, respectively. The QTLs on A7 (KNcqPH-A7-1, KNcqPH-A7-2, and KNcqSW-A7-3), A10 (KNcqPH-A10), C1 (KNcqPH-C1), C4 (KNcqPH-C4), and C9 (KNcqPH-C9-1, KNcqPH-C9-2, KNcqPH-C9-3, KNcqPH-C9-4, and KNcqPH-C9-5) were not observed in other genetic linkage groups. One orthologous gene TGH (Bra035958) on A6 was observed (Table 6), but no gene was aligned in the CI of QTLs for PH.
There were 34 QTLs for FBN projected onto the KN consensus map (Figure 5, Figure S3, Table S7), of which 10 were projected onto A3, and only two (KNcqFBN-A3-2 and BEqBN-A3-4) were co-located (Table S8). Two QTLs were projected onto A6, of which KNcqFBN-A6-3 and TNqBN-A6, and KNcqFBN-A6-5 and SEfb6.1 were co-localized. Three QTLs of KN (KNcqFBN-C3-1, KNcqFBN-C3-2, KNcqFBN-C3-3) and four QTLs of QN (QNqHFB-C3-2, QNqHFB-C3-1, QNqFB-C3-1, and QNqFB-C3-2) were co-localized. In addition, KNcqFBN-C5 and QNqFB-C5, and KNcqFBN-C6-1 and TNqBN-C6-1 were co-localized on C5 and C6, respectively. Some potential new QTLs of the KN population were obtained on A2 (KNcqFBN-A2), A7 (KNcqFBN-A7-1), A10 (KNcqFBN-A10-1 and KNcqFBN-A10-2), C2 (KNcqFBN-C2), C4 (KNcqFBN-C4), C7 (KNcqFBN-C7), and C9 (KNcqFBN-C9-1 and KNcqFBN-C9-2).
The traits BY, BH, LMI, and PMI have rarely been studied in other populations, thus, no QTL for PMI was projected onto the KN consensus map. Due to lack of common markers between KN and other maps, only a few QTLs for BY, BH, LMI, and PMI from five populations were projected onto the KN consensus map, including six QTLs for BH and eight each for BY and LMI (Figure S3, Table S8), respectively. However, none of these QTLs were co-located with the QTLs in KN population. In addition, two orthologous genes were observed, including HARDY (Bra017235) for BY and ATPAD4 (Bra006922) for BH (Table 6).
SY and SYRTs for B. napus are complex quantitative traits and easily affected by the environment (Quarrie et al., 2006). In previous studies, Chen et al. (2007) obtained 52 QTLs for six SYRTs and Fan et al. (2010) identified nine QTLs for SW. Chen et al. (2010) obtained 18 QTLs for SY and 22 QTLs for flowering time, while Butruille et al. (1999) found two QTLs for SY. These previous studies did not examine the relationship between SY and SYRTs, this was considered in the present study, showing correlation coefficients in the range of 0.31–0.83, except for FBN.
Several studies have revealed that only a few QTLs for SYRTs were stable in different environments (Li et al., 2007; Shi et al., 2009). The majority of consensus QTLs were also environment specific in the present study, including 70.8 and 23.6% of QTLs there were only detected in winter and spring areas, respectively. Only a few QTLs (5.6%) were detected in both areas, consistent with previous reports. These findings will be helpful for breeders to develop varieties with special adaptability.
Trait-by-trait meta-analysis revealed that 144 consensus QTLs for SY and SYRTs were integrated into 72 pleiotropic unique QTLs. Among them, six of seven unique QTLs for SY were co-located with two to five QTLs for SYRTs. On average, one unique QTL for SY involved 2.5 QTLs for SYRTs (Table 5). For example, uq.A2-3 affected SY as well as five SYRTs: SW, BH, PH, FBN, and LMI. These findings were consistent with the strong correlations between SY and SYRTs. Shi et al. (2009) demonstrated that the QTLs for SY were pleiotropic and synthesized, and numerous SYRTs were potential contributor to tightly link with the QTLs for SY. Li et al. (2007) revealed that QTLs for SY and SYRTs usually had overlapping regions. Similar results were also reported for spring barley (Li et al., 2005) and red clover (Herrmann et al., 2006). The genes identified as having a function for SY and SYRTs also had pleiotropic effects for at least one trait (Ashikari et al., 2005; Lim et al., 2014) or multiple traits (Li et al., 1997; Hall et al., 2005; Quarrie et al., 2006; Burgess-Herbert et al., 2008; Hu et al., 2008). In other words, QTLs for SY and SYRTs might have resulted from pleiotropic QTLs that controlled multiple traits by containing multiple, closely linked, trait-specific genes (Hall et al., 2005).
Indicator QTLs have been successfully used in identifying genes with pleiotropic effects for SY (Shi et al., 2009). Indicator QTLs for SY must be stably expressed, easily measured, and identifiable candidate genes than the co-localized QTLs for SY. Using the similar methods, the candidate genes were successfully cloned by the QTL for BY in wheat (Quarrie et al., 2006), and the QTL for flowering time could be regarded as an indicator QTL for natural variation in rice (Xue et al., 2008). In the present study, four QTLs for SY closely co-existed with QTLs for SW (Table 5). SW can be more easily and precisely measured than SY, and is also less influenced by environmental factors than other SYRTs (Shi et al., 2009; Ding et al., 2012). Therefore, QTLs for SW could be regarded as indicator QTLs for SY. In future research, cloning the indicator QTL KNcqSW-C6-3 for SY is more feasible than cloning KNcqSY-C6-3.
In the present study, 166 QTLs for SY and SYRTs from five populations were projected onto the constructed high-density consensus map. However, only 36 QTLs were co-located with QTLs identified in the KN population, including six QTLs for SY from the TN population and 11 for SW from the TN and BE population. Meanwhile, the co-located QTLs for SY were located on A2 and C6, while the co-located QTLs for SW were mainly located on A3 and C6. Thus, C6 was an important linkage group and included important genes for SY and SW. The reason that only a few QTLs were co-localized might be the lack of sufficient common markers between the different maps (Zhou et al., 2014) and also that there has been little research for some traits, such as PMI and LMI.
The known Arabidopsis genome sequence has been exploited as a tool for comparative analysis between Arabidopsis and Brassicaceae genomes, and the conserved genomic block has been identified in different Brassicaceae species (Arcade et al., 2004; Boivin et al., 2004; Schranz et al., 2006). This provided a method to align candidate genes underlying QTLs controlling important agronomic traits (Long et al., 2007). For SY, the orthologous gene LQY1 was mapped on the C6 chromosome, which encodes a small zinc-finger-containing thylakoid membrane protein of Arabidopsis (Lu, 2011; Jin et al., 2014). Gene LQY1 resulted in a lower quantum yield of photosystem II (PSII) photochemistry and reduced PSII electron transport rate following high-light treatment (Lu et al., 2011). For SW, gene AP2 (Bra011741 and BnaA03g53830D) and PDF1 (Bol039878 and BnaC06g22110D) were underlying the CIs of KNcqSW-A3-3 and KNcqSW-C6-3, respectively. Gene AP2 was involved in the specification of floral organ identity, establishment of floral meristem identity, ovule and seed coat development, and also had a role in controlling seed mass (Ohto et al., 2009). Gene PDF1 is known as plant defense type 1 gene, and conferred high capacities to tolerate and hyperaccumulate zinc and cadmium (Mirouze et al., 2006). In this study, gene PDF1 located on C6 controlled seed size and weight; however, the genetic mechanisms of controlling seed size and weight remained ambiguous. These genes are speculated to be candidate genes for SY or SYRTs in B. napus.
The genetic mechanisms underlying SY and SYRTs were analyzed through QTL analysis in B. napus. A total of 226 QTLs were identified, and were integrated into 144 consensus QTLs. Seven major QTLs were obtained, including three QTLs for SY, two for SW and one each for BH and FBN, respectively. Trait-by-trait meta-analysis revealed that one unique QTL for SY involved 2.5 QTLs for SYRTs. Meanwhile, QTL projection from five different genetic linkage maps onto the KN consensus map showed that 36 QTLs were co-localized with QTLs identified in the KN population. In addition, candidate genes for SY and SYRTs were observed, including five each for SY (GASA4, ATCLH1, RBCS1A, LQY1, and ATGGH1) and SW (PTH2, AP2, LCR64, LCR65, and PDF1), and one each for BY (HARDY), BH (ATPAD4), and PH (TGH). The obtained candidate genes for SY and SYRTs were conducive to fine mapping and key gene cloning. These findings will be valuable in hybrid cultivar breeding and in analyzing QTL expression in different environments.
WZ and XW carried out the QTL analysis and wrote the manuscript. JT, BL, LC, and HC participated in the field experiment. YL, JX, JG, and WL made helpful suggestions to the manuscript. HW and ML designed, led and coordinated the overall study.
This work was supported financially by the National Science Foundation of China (31471532), the National Basic Research Program of China (2015CB150205), International Cooperation in Science and Technology Projects (2014DFA32210), the New Century Talents Support Program of the Ministry of Education of China (NCET110172), the Natural Science Foundation Research Key Projects of Shaanxi Province of China (2013JZ005) and International Cooperation Key Projects of Shaanxi Province of China (2012KW-15).
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fpls.2016.00017
B. napus, Brassica napus; B. rapa, Brassica rapa; B. oleracea, Brassica oleracea; QTL, Quantitative trait locus; DH, Doubled haploid; cM, CentiMorgan; CIs, Confidence intervals; PVE, Phenotypic variation explained; SY, Seed yield; SYRTs, Seed yield related traits; BY, Biomass yield; SW, Thousand seed weight; BH, First effective branch height; PH, Plant height; FBN, First effective branch number; LMI, Length of main inflorescence; PMI, Pod number of main inflorescence; LOD, Log odds score methods.
Arcade, A., Labourdette, A., Falque, M., Mangin, B., Chardon, F., Charcosset, A., et al. (2004). BioMercator: integrating genetic maps and QTL towards discovery of candidate genes. Bioinformatics 20, 2324–2326. doi: 10.1093/bioinformatics/bth230
Basunanda, P., Radoev, M., Ecke, W., Friedt, W., Becker, H. C., and Snowdon, R. J. (2010). Comparative mapping of quantitative trait loci involved in heterosis for seedling and yield traits in oilseed rape (Brassica napus L.). Theor. Appl. Genet. 120, 271–281. doi: 10.1007/s00122-009-1133-z
Boivin, K., Acarkan, A., Mbulu, R. S., Clarenz, O., and Schmidt, R. (2004). The Arabidopsis genome sequence as a tool for genome analysis in Brassicaceae. A comparison of the Arabidopsis and capsella rubella genomes. Plant Physiol. 135, 735–744. doi: 10.1104/pp.104.040030
Burgess-Herbert, S., Cox, A., Tsaih, S. W., and Paigen, B. (2008). Practical applications of the bioinformatics tool box for narrowing quantitative trait loci. Genetics 180, 2227–2235. doi: 10.1534/genetics.108.090175
Butruille, D. V., Guries, R. P., and Osborn, T. C. (1999). Linkage analysis of molecular markers and quantitative trait loci in populations of inbred backcross lines of Brassica napus L. Genetics 153, 949–964.
Cai, D., Xiao, Y., Yang, W., Ye, W., Wang, B., Younas, M., et al. (2014). Association mapping of six yield-related traits in rapeseed (Brassica napus L.). Theor. Appl. Genet. 127, 85–96. doi: 10.1007/s00122-013-2203-9
Chen, G., Geng, J., Rahman, M., Liu, X., Tu, J., Fu, T., et al. (2010). Identification of QTL for oil content, seed yield, and flowering time in oilseed rape (Brassica napus). Euphytica 175, 161–174. doi: 10.1007/s10681-010-0144-9
Chen, W., Zhang, Y., Liu, X., Chen, B., Tu, J., Fu, T., et al. (2007). Detection of QTL for six yield-related traits in oilseed rape (Brassica napus) using DH and immortalized F2 populations. Theor. Appl. Genet. 115, 849–858. doi: 10.1007/s00122-007-0613-2
Ding, G., Liao, Y., Yang, M., Zhao, Z., Shi, L., and Xu, F. (2011). Development of gene-based markers from Arabidopsis thaliana functional genes involved in phosphorpus homeostasis and mapping in Brassica napus. Euphytica 181, 305–322. doi: 10.1007/s10681-011-0428-8
Ding, G., Zhao, Z., Liao, Y., Hu, Y., Shi, L., Long, Y., et al. (2012). Quantitative trait loci for seed yield and yield-related traits, and their responses to reduced phosphorus supply in Brassica napus. Ann. Bot. 109, 747–759. doi: 10.1093/aob/mcr323
Fan, C., Cai, G., Qin, J., Li, Q., Yang, M., Wu, J., et al. (2010). Mapping of quantitative trait loci and development of allele-specific markers for seed weight in Brassica napus. Theor. Appl. Genet. 121, 1289–1301. doi: 10.1007/s00122-010-1388-4
Hall, M. C., Basten, C. J., and Willis, J. H. (2005). Pleiotropic quantitative trait loci contribute to population divergence in traits associated with life-history variation in mimulus guttatus. Genetics 172, 1829–1844. doi: 10.1534/genetics.105.051227
Herrmann, D., Boller, B., Studer, B., Widmer, F., and Kolliker, R. (2006). QTL analysis of seed yield components in red clover (Trifolium pratense L.). Theor. Appl. Genet. 112, 536–545. doi: 10.1007/s00122-005-0158-1
Hu, K., Qiu, D., Shen, X., Li, X., and Wang, S. (2008). Isolation and manipulation of quantitative trait loci for disease resistance in rice using a candidate gene approach. Mol. Plant 1, 786–793. doi: 10.1093/mp/ssn039
Jiang, C., Shi, J., Li, R., Long, Y., Wang, H., Li, D., et al. (2014).Quantitative trait loci that control the oil content variation of rapeseed (Brassica napus L.). Theor. Appl. Genet. 127, 957–968. doi: 10.1007/s00122-014-2271-5
Jin, H., Liu, B., Luo, L., Feng, D., Wang, P., Liu, J., et al. (2014). Hypersensitive to high light1 interacts with low quantum yield of photosystem II1 and functions in protection of photosystem II from photodamage in Arabidopsis. Plant Cell 26, 1213–1229. doi: 10.1105/tpc.113.122424
Li, J., Huang, X., Heinrichs, F., Ganal, M. W., and Roder, M. S. (2005). Analysis of QTLs for yield, yield components, and malting quality in a BC3-DH population of spring barley. Theor. Appl. Genet. 110, 356–363. doi: 10.1007/s00122-004-1847-x
Li, Y., Shen, J., Wang, T., Chen, Q., Zhang, X., Fu, T., et al. (2007). QTL analysis of yield-related traits and their association with functional markers in Brassica napus L. Aust. J. Agr. Res. 58, 759. doi: 10.1071/AR06350
Lim, J., Yang, H., Jung, K., Yoo, S. C., and Paek, N. C. (2014). Quantitative trait locus mapping and candidate gene analysis for plant architecture traits using whole genome re-sequencing in rice. Mol. Cells 37, 149–160. doi: 10.14348/molcells.2014.2336
Liu, S., Liu, Y., Yang, X., Tong, C., Edwards, D., Parkin, I. A., et al. (2014). The Brassica oleracea genome reveals the asymmetrical evolution of polyploid genomes. Nat. Commun. 5:3930. doi: 10.1038/ncomms4930
Lombard, V., and Delourme, R. (2001). A consensus linkage map for rapeseed (Brassica napus L.): construction and integration of three individual maps from DH populations. Theor. Appl. Genet. 103, 491–507. doi: 10.1007/s001220100560
Long, Y., Shi, J., Qiu, D., Li, R., Zhang, C., Wang, J., et al. (2007). Flowering time quantitative trait loci analysis of oilseed Brassica in multiple environments and genomewide alignment with Arabidopsis. Genetics 177, 2433–2444. doi: 10.1534/genetics.107.080705
Lu, Y., Hall, D. A., and Last, R. L. (2011). A small Zinc finger thylakoid protein plays a role in maintenance of Photosystem II in Arabidopsis thaliana. Plant Cell 23, 1861–1875. doi: 10.1105/tpc.111.085456
Maccaferri, M., Sanguineti, M. C., Corneti, S., Ortega, J. L. A., Salem, M. B., Bort, J., et al. (2008). Quantitative trait loci for grain yield and adaptation of durum wheat (Triticum durum Desf.) across a wide range of water availability. Genetics 178, 489–511. doi: 10.1534/genetics.107.077297
Mirouze, M., Sels, J., Richard, O., Czernic, P., Loubet, S., Jacquier, A., et al. (2006). A putative novel role for plant defensins: a defensin from the zinc hyper-accumulating plant, Arabidopsis halleri, confers zinc tolerance. Plant J. 47, 329–342. doi: 10.1111/j.1365-313X.2006.02788.x
Ohto, M. A., Floyd, S. K., Fischer, R. L., Goldberg, R. B., and Harada, J. J. (2009). Effects of APETALA2 on embryo, endosperm, and seed coat development determine seed size in Arabidopsis. Sex. Plant Reprod. 22, 277–289. doi: 10.1007/s00497-009-0116-1
Qiu, D., Morgan, C., Shi, J., Long, Y., Liu, J., Li, R., et al. (2006). A comparative linkage map of oilseed rape and its use for QTL analysis of seed oil and erucic acid content. Theor. Appl. Genet. 114, 67–80. doi: 10.1007/s00122-006-0411-2
Quarrie, S., Pekic, Q. S., Radosevic, R., Kaminska, A., Barnes, J. D., Leverington, M., et al. (2006). Dissecting a wheat QTL for yield present in a range of environments: from the QTL to candidate genes. J. Exp. Bot. 57, 2627–2637. doi: 10.1093/jxb/erl026
Quijada, P. A., Udall, J. A., Lambert, B., and Osborn, T. C. (2006). 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. 113, 549–561. doi: 10.1007/s00122-006-0323-1
Radoev, M., Becker, H. C., and Ecke, W. (2008). Genetic analysis of heterosis for yield and yield components in Rapeseed (Brassica napus L.) by quantitative trait locus mapping. Genetics 179, 1547–1558. doi: 10.1534/genetics.108.089680
Raman, H., Raman, R., Eckermann, P., Coombes, N., Manoli, S., Zou, X., et al. (2013). Genetic and physical mapping of flowering time loci in canola (Brassica napus L.). Theor. Appl. Genet. 126, 119–132. doi: 10.1007/s00122-012-1966-8
Schranz, M. E., Lysak, M. A., and Mitchell-Olds, T. (2006). The ABC's of comparative genomics in the Brassicaceae: building blocks of crucifer genomes. Trends Plant Sci. 11, 535–542. doi: 10.1016/j.tplants.2006.09.002
Scoles, G., Fu, F., Liu, L., Chai, Y., Chen, L., Yang, T., et al. (2007). Localization of QTLs for seed color using recombinant inbred lines of Brassica napus in different environments. Genome 50, 840–854. doi: 10.1139/G07-068
Shi, J., Li, R., Qiu, D., Jiang, C., Long, Y., Morgan, C., et al. (2009). Unraveling the complex trait of crop yield with quantitative trait loci mapping in Brassica napus. Genetics 182, 851–861. doi: 10.1534/genetics.109.101642
Shi, L., Shi, T., Broadley, M. R., White, P. J., Long, Y., Meng, J., et al. (2013). High-throughput root phenotyping screens identify genetic loci associated with root architectural traits in Brassica napus under contrasting phosphate availabilities. Ann. Bot. 112, 381–389. doi: 10.1093/aob/mcs245
Udall, J. A., Quijada, P. A., Lambert, B., and Osborn, T. C. (2006). Quantitative trait analysis of seed yield and other complex traits in hybrid spring rapeseed (Brassica napus L.): 2. Identification of alleles from unadapted germplasm. Theor. Appl. Genet. 113, 597–609. doi: 10.1007/s00122-006-0324-0
Wang, X., Wang, H., Long, Y., Li, D., Yin, Y., Tian, J., et al. (2013). 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 8:e80569. doi: 10.1371/journal.pone.0080569
Wang, X., Wang, H., Long, Y., Liu, L., Zhao, Y., Tian, J., et al. (2015). Dynamic and comparative QTL analysis for plant height in different developmental stages of Brassica napus L. Theor. Appl. Genet. 128, 1175–1192. doi: 10.1007/s00122-015-2498-9
Xue, W., Xing, Y., Weng, X., Zhao, Y., Tang, W., Wang, L., et al. (2008). Natural variation in Ghd7 is an important regulator of heading date and yield potential in rice. Nat. Genet. 40, 761–767. doi: 10.1038/ng.143
Zhang, G., He, Y., Xu, L., Tang, G., and Zhou, W. (2006). Genetic analyses of agronomic and seed quality traits of doubled haploid population in Brassica napus through microspore culture. Euphytica 149, 169–177. doi: 10.1007/s10681-005-9064-5
Keywords: Brassica napus, seed yield, seed yield-related traits, quantitative trait loci, map comparsion, candidate genes
Citation: Zhao W, Wang X, Wang H, Tian J, Li B, Chen L, Chao H, Long Y, Xiang J, Gan J, Liang W and Li M (2016) 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. 7:17. doi: 10.3389/fpls.2016.00017
Received: 25 August 2015; Accepted: 08 January 2016;
Published: 28 January 2016.
Edited by:Juan Francisco Jimenez Bremont, Instituto Potosino de Investigacion Cientifica y Tecnologica, Mexico
Reviewed by:Yinglong Chen, The University of Western Australia, Australia
Daniela Marone, Cereal Research Centre - CRA-CER, Italy
Copyright © 2016 Zhao, Wang, Wang, Tian, Li, Chen, Chao, Long, Xiang, Gan, Liang and Li. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
†These authors have contributed equally to this work.