Skip to main content

Diverse regulatory factors associate with flowering time and yield responses in winter-type Brassica napus



Flowering time, plant height and seed yield are strongly influenced by climatic and day-length adaptation in crop plants. To investigate these traits under highly diverse field conditions in the important oilseed crop Brassica napus, we performed a genome-wide association study using data from diverse agroecological environments spanning three continents.


A total of 158 European winter-type B.napus inbred lines were genotyped with 21,623 unique, single-locus single-nucleotide polymorphism (SNP) markers using the Brassica 60 K-SNP Illumina® Infinium consortium array. Phenotypic associations were calculated in the panel over the years 2010–2012 for flowering time, plant height and seed yield in 5 highly diverse locations in Germany, China and Chile, adding up to 11 diverse environments in total.


We identified 101 genome regions associating with the onset of flowering, 69 with plant height, 36 with seed yield and 68 cross-trait regions with potential adaptive value. Within these regions, B.napus orthologs for a number of candidate adaptation genes were detected, including central circadian clock components like CIRCADIAN CLOCK- ASSOCIATED 1 (Bna.CCA1) and the important flowering-time regulators FLOWERING LOCUS T (Bna.FT) and FRUITFUL (Bna.FUL).


Gene ontology (GO) enrichment analysis of candidate regions suggested that selection of genes involved in post-transcriptional and epigenetic regulation of flowering time may play a potential role in adaptation of B. napus to highly divergent environments. The classical flowering time regulators Bna.FLC and Bna.CO were not found among the candidate regions, although both show functional variation. Allelic effects were additive for plant height and yield, but not for flowering time. The scarcity of positive minor alleles for yield in this breeding pool points to a lack of diversity for adaptation that could restrict yield gain in the face of environmental change.


Our study provides a valuable framework to further improve the adaptability and yield stability of this recent allopolyploid crop under changing environments. The results suggest that flowering time regulation within an adapted B. napus breeding pool is driven by a high number of small modulating processes rather than major transcription factors like Bna.CO. In contrast, yield regulation appears highly parallel, therefore yield could be increased by pyramiding positively associated haplotypes.


Environmental adaptation is decisive for optimal yield of major crop plants. This includes adaptation of the reproductive system to the prevailing climatic conditions (i.e., season, day length, and temperature) and an appropriate response to biotic and abiotic stresses. Adaptation of the reproductive system is mainly achieved by setting the time for flowering, which is considered to be the transition from the vegetative to the generative phase of the plant lifecycle. This switch determines how much the plant invests into vegetative growth and, consequently, how many nutrients will be available for remobilization at the time of seed filling [1]. This control occurs both directly, depending on the leaf area and/or the root architecture, and indirectly, via synchronization with day length, water availability and ambient temperature. In temperate regions, flowering time also determines the time left for seed maturation before winter. The complex task of finding the optimal time to switch is controlled by the flowering gene network, representing a set of interacting pathways of transcription factors, photosensors, enzymes and miRNAs (reviewed in depth in [28]).

In Arabidopsis thaliana it has been shown before that optimal flowering time is correlated with higher seed yield [9]. A recent study in Brassica rapa identified B.rapa FLOWERING LOCUS C (Bra.FLC) as a candidate for a quantitative trait locus (QTL) contributing to transgenerative fitness, i.e., the influence of flowering time of the mother plant on germination of the offspring [10]. In Brassica napus, earlier QTL studies showed colocalization of QTL for flowering and traits like yield and plant height [1113], suggesting that QTL-associated flowering genes like B.napus FLOWERING LOCUS C (Bna.FLC), FLOWERING LOCUS T (Bna.FT), VERNALISATION INSENSITIVE 3 (Bna.VIN3) and ZEITLUPE (Bna.ZTL) have a major impact on these traits. This is indicative for a more general role of flowering regulators in plant development, making them global regulators of plant adaptation. Such a role has been demonstrated in A.thaliana for Ath.FLC [14], CRYPTOCHROME 2 (Ath.CRY2), FRIGIDA (Ath.FRI) and PHYTOCHROME C (Ath.PHYC) [1].

Brassica napus (rapeseed, oilseed rape, canola) is a recent allotetraploid crop species with high economic impact in Europe, Asia, Australia and North America [15]. Its huge phenotypic variation groups the species into different morphotypes with distinct flowering behavior, vernalisation requirement and winter hardiness. Being closely related to the model crucifer A.thaliana, B.napus serves as appropriate model for effects of allopolyploidy [1619]. It harbors two related subgenomes from a recent interspecific hybridization (A subgenome with 10 chromosomes from B.rapa and C subgenome with 9 chromosomes from B. oleracea). Polyploidy, although widely spread among crop species, poses special challenges for molecular breeding, as the similarity of subgenomes can impede development of locus-specific markers, cloning of candidate genes or knock-out techniques. Recently developed resources like the Brassica 60 k SNP Illumina® Infinium consortium array and the publication of the Brassica napus genomic sequence [20] now allow fine-mapping of genotype-trait associations on a genome-wide scale.

The genetics of flowering, plant development and yield have been studied extensively in Brassica crop species by QTL analysis in biparental populations [1113, 2126] and association genetics [2729]. Biparental QTL studies generally suffer from their low resolution and only assess alleles present in the mapping parents; therefore the results may not be applicable within or across breeding pools. Genome-wide association studies overcome these restrictions, given reasonable marker saturation and an appropriately diverse population. Cai and co-authors [27] used 1191 markers to study 192 inbred lines to map yield-related traits, whereas Li and co-authors [28] increased the number of markers almost 20-fold by using the Brassica 60 k SNP array to study seed weight and quality traits in 472 diverse B.napus lines. High-resolution identification of the well-known FATTY ACID ELONGASE 1 (Bna.FAE1) homologs controlling erucic acid content demonstrated the power of the GWAS approach to map key genes for important traits using densely-spaced genome-wide SNP markers.

Genes for global plant adaptation are highly interesting markers for plant breeding. Biparental QTL studies have demonstrated that variants of the major flowering-time regulatory genes like FLOWERING LOCUS C (Bna.FLC) and FLOWERING LOCUS T (Bna.FT) are the major causes for the broad ecogeographical diversification within B.napus into spring-type, winter-type and semi-winter rapeseed gene pools [12, 22, 26, 30]. Because of the extremely strong phenotypic effects of these few loci, in across-pool populations the genetic effects of loci with regulatory roles within each morphotype tend to be masked, making it difficult to reveal subtle effects of floral regulators on local adaptation. These regulators are often important for breeding because they allow breeders to adjust yield-related traits within adapted gene pools. Association mapping within pools can potentially reveal allelic variation in regulatory loci responsible for local adaptability, for example to environmental variation or climate change.

To reveal if selected members of the flowering network are valuable regulatory factors of plant adaptability, we phenotyped a set of 158 diverse winter-type B.napus inbred lines in extremely different environments in three location in Germany (2010–2012), and one location in Chile (2011–2012) and China (2012). Phenotypic data was collected for onset of flowering, plant height and yield and used to calculate genome-wide associations with genotype data from the Brassica 60 k SNP array. By investigating the gene content and annotations within chromosome regions in linkage disequilibrium (LD) with QTL for these three traits, we aimed to gain insight into the mechanisms of adaptation to fluctuations in photoperiodicity and climatic fluctuation in this important global crop.


Field trials

The beginning of flowering (BOF) varied considerably among years and locations, spanning over a range of 97 days from 157 days after sowing (DAS) in Temuco (Chile), 2012, to 254 DAS in Rauischholzhausen (Germany) 2010. The order of flowering locations for earliest to latest flowering was Temuco (Chile)/Beibei (China)/Gross Gerau/Rauischholzhausen/Giessen (all Germany), whereby the difference was not significant between Giessen and Rauischholzhausen in 2011 (see Fig. 1, Table 1). Although this order was retained during the years 2010–2012, absolute flowering time values of each location also varied significantly between years, reflecting clearly different climatic conditions regarding, temperature and precipitation from 2010 to 2012 and varying day length between the locations (see Fig. 2). Correlations of flowering time between locations show that flowering behaviour was most similar between Gross Gerau and Rauischholzhausen, and most divergent between Giessen and Temuco. This most probably reflects differences in precipitation and day length (Fig. 2). Plant height and seed yield were recorded in all of the German environments, while plant height was also recorded in Temuco. Correlations between flowering time and plant height were always significantly positive (p = 0.05) with the exception of Temuco 2012, which was positive, but not significant. Correlations between flowering time and seed yield were negative (with the exception of Giessen 2011) and significantly so for Gross Gerau (all years) and Rauischholzhausen 2012. Correlations between plant height and yield showed no trend and were not significant (correlations see Additional file 1: Table S1).

Fig. 1

Distribution of flowering time. Values in Days after sowing (DAS) for all 11 environments in 3 years (2010–2012) and 5 locations (Germany: Giessen, Gross Gerau, Rauischholzhausen; Chile: Temuco; China: Beibei, Chongqing). The onset of flowering was attributed when 10 % of the plants in a plot reached BBCH stage 61

Table 1 Phenotypic means for all locations and years
Fig. 2

Climate data for all environments tested. The graphs show mean values (daylength, temperature) and monthly sums (precipitation) for each location (daylength) and each environment (temperature, precipitation). The environments are abbreviated as follows: Germany: Giessen (GI), Gross Gerau (GG), Rauischholzhausen (RH); Chile: Temuco (TE); China: Beibei, Chongqing (BB), with their respective year of harvest. For precipitation, consecutive years are placed consecutively beneath each other

Population structure

PCA plots using 3 clusters for components 1/2 and 2/3 are shown in Additional file 2: Figure S1. The first four principal components explained 7.04, 6.46, 4.33 and 3.85 % of the genetic variation, respectively. This revealed slight population structure indicated by 3 clusters (see Additional file 2: Figure S1). Identity by descent (IBD) values among clusters 1 and 2 were higher than for the total population, while cluster 3 showed slightly lower IBD values than the total population. This cluster includes three winter-type fodder rapes and one exotic accession and shows greater genetic variance than the population average. F ST values among cluster pairs revealed highest relatedness between clusters 2 and 3, whereas cluster 1 was approximately equally distant from clusters 2 and 3. This means that cluster 1 clearly separates from clusters 2 and 3, whereas clusters 2 and 3 separate into a defined subpopulation with higher relatedness (cluster 2) and a less defined subpopulation with more diverse material (cluster 3). Total population F ST was found to be 0.1273, indicating a medium to high relatedness across the overall population. Based on these findings and on further testing of population stratification, GWAS was conducted using a mixed model approach with PC-adjustment.

Linkage disequilibrium

Mean pairwise r 2 for all chromosomes was found to be 0.105. When comparing the decay of LD among the two subgenomes, it is evident that the C subgenome shows a much slower average LD decay than the A subgenome (Fig. 3). Whereas LD in the A subgenome drops to less than 0.2 within 494.8 kbp, in the C subgenome LD does not drop under 0.2 until an average of 3389.5 kbp. When investigating into this phenomenon by chromosome, it becomes evident that chromosomes C01 and C04 are the major causes for the slow LD decay in the C genome, with additional contributions from chromosomes C07 and C08 (Fig. 3). This means that chromosomes C02, C03, C05 and C09 show a behavior similar to the A genome. LD decay on chromosome C08 shows an unexpected shape with a clearly visible peak centered at 5256.0 kbp, which might represent a signature of co-selection. On chromosome C07, there is a similar, but less pronounced peak centered on 5414.1 kbp.

Fig. 3

Patterns of linkage disequilibrium (LD) decay. Left: the total genome compared to the respective A and C subgenomes, Right: each of the 19 B.napus chromosome separately. Chromosomes C01 and C04 are indicated by arrows

Association analysis

Onset of flowering

Association analysis for beginning of flowering (BOF) in the context of absolute flowering time (BOF-ABS) and relative flowering time (BOF-REL, see Methods for definition) defined 62 and 47 regions of interest, respectively, 8 of which overlapped. Examples of Manhattan plots for BOF-ABS are shown in Fig. 4, along with their quantile-quantile (QQ) plots. All regions were compared to the positions of known B. napus gene models which were published in [20]. Gene ontology terms were queried for flowering-related terms like “flower development”, “vegetative to reproductive transition”, “photoperiod” or “vernalisation” (see methods for full list). For BOF-ABS a total of 181 flowering related genes were identified in 49 candidate regions, while 13 regions were left unexplained by flowering related gene ontology. For BOF-REL, 232 flowering related genes were identified in 38 regions of interest, leaving 9 regions unexplained. A total of 39 genes overlapped between both analyses. A full list of all candidate genes with their closest and highest associating SNP markers, including their phenotypic effects, is shown in Additional file 3: Table S2. For BOF-ABS, the list contains 31 transcription factors, 2 translation initiation factors, 1 translation elongation factor and 18 protein kinases. Furthermore, 2 histones and 11 histone interacting factors were among the candidates. For BOF-REL we found 42 transcription factors, 3 translation initiation factors and 19 protein kinases, as well as 25 histone interacting factors.

Fig. 4

Manhattan plots and QQ plots. Describing marker-trait associations for absolute values of the onset of flowering (BOF-ABS) in the three divergent environments Gross Gerau, Rauischholzhausen (both Germany) and Temuco (Chile) in 2011. Associations were calculated using a mixed model with PC adjustment, using 21,623 single-copy, genome-wide SNP markers. The solid line represents the significance threshold (p = 0.1) for associations after Bonferroni correction, the dashed line represents the level considered as significant in this analysis (−log10(p) > 4), and the dotted line represents the threshold for suggestive associations(−log10(p) > 3). The A and C subgenomes are displayed in a slightly different color code (dark: A subgenome, bright: C subgenome)

Among the detected transcription factors within BOF-ABS candidate regions, orthologs to 7 A.thaliana agamous-like (AGL) transcription factors were found, 6 of which are known to play a role in determining floral meristem identity [31]. In the analysis of BOF-REL we detected 8 orthologs to agamous-like transcription factors as well as a copy of the transcriptional co-repressor LEUNIG (BnaC01g05930D) which is known to negatively regulate AGAMOUS (BnaC03g62970D, found for BOF-ABS). Copies of Bna.CO and Bna.FLC were not detected among the candidate regions, although an ortholog of FRIGIDA (BnaA03g13320D), a central vernalisation regulator, was detected in the Chinese environment in BOF-REL, in relative proximity to a copy of FLC. A summary of interesting candidate genes, including the abovementioned, is given in Table 2. Overall, we detected only rare examples of early-acting flowering regulators, while late-acting flowering regulators and floral homeotic genes were more predominant among the associated regions.

Table 2 Examples of candidate genes within LD blocks associated with absolute flowering time (BOF-ABS) and relative flowering time (BOF-REL)

Plant height

For plant height (HEI-ABS), 69 regions of interest were found to be associated with the available phenotype data from the environments in Chile and Germany. Gene ontology terms were searched for terms like “photosynthesis”, “catabolic”, “cell growth” or “development”, resulting in a list of 712 genes covering 65 from 69 regions of interest (see methods for full list of GO terms). Out of the regions of interest for plant height, 17 overlapped with regions associated with BOF-ABS, containing for example copies of the agamous-like transcription factors AGL32 (BnaA09g05410D), FRUITFUL (FUL; BnaA09g05500D, BnaC02g41870D) and FOREVER YOUNG FLOWER (FYF; BnaA09g05900D) (Fig. 5). Overall, we detected 8 AGL transcription factors, including a copy of SEPALLATA 3 (SEP3; BnaA09g28040D) differing from the copy found for BOF-ABS (BnaC05g21200D). The candidate regions also contained orthologs of the flowering genes GIGANTEA (GI; BnaA09g30390D), TEMPRANILLO1 (TEM1; BnaA07g35130D) and FPA (BnaC08g29740D). Moreover, we found orthologs to the squamosa-promotor binding-like proteins SPL1 (BnaA05g00780D), SPL4 (BnaC06g10070D), SPL10 (BnaA09g27950D) and SPL11(BnaA09g27960D). To summarize, a high overlap of associating regions was observed between flowering time and plant height. Moreover, some regions associated with plant height contained flowering-related genes for which associations to flowering time itself were not detected.

Fig. 5

Local representation of marker-trait associations with absolute values for flowering time (BOF-ABS), plant height (HEI-ABS) and seed yield (YIE-ABS) in a region of B.napus chromosome A09. Each environment is represented in a unique colour consistent between representations, as shown by the colour key. The solid line represents the significance threshold 0.1 after Bonferroni correction, the dashed line represents the significance threshold in this analysis (−log10(p) > 4), and the dotted line represents the threshold for suggestive associations (−log10(p) > 3). The positions of 5 candidate genes within the peak interval are indicated by black squared dots. The representation below shows the local conservation of linkage disequilibrium (r 2) among the SNPs in the trait-associated region is represented by a heat map


For yield (YIE-ABS), 36 regions of interest were defined by the associations detected with the available phenotype data from the German environments (no data available for Chile and China). Of those, 6 regions overlapped with BOF-ABS and 4 with HEI-ABS, while 1 region on chromosome A09 overlapped between all of them (Fig. 5). Gene ontology terms corresponding to the genes among these regions were searched for terms like “lipid”, “terpenoid”, “biosynthetic process” or “embryo (sac) development” (see methods for full list of GO terms). The 568 detected genes covered 34 from 36 regions of interest. Noteably, 31 % of the candidate genes within these regions had an orthologous copy in another candidate region, whereas for flowering and plant height this was the case for only 8 % (BOF-ABS), 9 % (BOF-REL) and 9 % (HEI-ABS) of the detected genes. This was mainly due to several regions of interest being homoeologous, for example on chromosomes A03/A09/C09, A04/C04 and A05/C05, but also due to adjacent duplicated genes within the same region. Homoeology between associated regions was not observed for BOF-ABS, and only a very short segment on chromosomes A02/C02 was found to be homoeologous between regions associated with BOF-REL (containing the two copies of LUMINIDEPENDENS; BnaA02g20680D and BnaC02g27000D). HEI-ABS showed 3 regions homoeologous to each other on chromosomes A03/C03, A06/C05 and A09/C02, respectively. Examples of duplicate copies of a gene or gene family are given in Table 3.

Table 3 Examples of duplicated candidate genes within LD blocks associated with yield

Some classical flowering genes were also found among those regions, for example two copies of FUL (BnaA03g39830D, BnaA09g05510D) and FLOWERING PROMOTING FACTOR 1 (FPF1; BnaA09g04780D, BnaC09g04250D), and one copy each of TFL1 (BnaA10g26300D) and FT (BnaC04g14850D) along with two copies of PSEUDO RESPONSE REGULATOR 5 (PRR5; BnaA09g04840D, BnaC09g04380D) and one copy of PRR1 (BnaC09g05250D).

Overall, homoeologous chromosome regions were commonly associated to yield, and many of these regions carry copies of flowering-time candidate genes. This indicates a parallel selection among B. napus subgenomes.

Multi-trait associations

For multi-trait comparative analysis (COMP) we defined 68 regions of interest, taking into account the relative phenotypic values only for German environments in which all traits were measured. The regions of interest showed at least 1 suggestive association for one trait and at least 1 supporting association for another trait. We considered genes with GO terms referring to all traits under study as described before. This identified 739 candidate genes covering 67 from 68 regions of interest. Among the candidates were orthologs to well-known components of the circadian clock like CIRCADIAN CLOCK ASSISTED 1 (CCA1; BnaA05g01050D) and TIME FOR COFFEE (TIC; BnaC03g42780D), which both occurred in a region inversely affecting plant height and yield in Gross Gerau in 2011. We moreover detected a copy of FT that was inversely associated with plant height and yield in Rauischholzhausen in 2011 (BnaA07g25310D). The 3 AGL transcription factors Bna.AGL6, Bna.AGL8/FUL, Bna.AGL65 were again detected. AGL6 is a regulator of FLC and FT [32], while AGL65 is necessary for pollen development [33]. AGL6 and FUL (the same copy as in YIE-ABS) are found in regions affecting all traits.

Copies of the he three SPL orthologs SPL1, SPL4 and SPL8 (BnaA10g00110D) were also detected in COMP. For SPL1 and SPL4 the same copies were detected as for HEI-ABS, while SPL8 is known to regulate anther development. Moreover, we detected a copy of LIGHT-DEPENDENT SHORT HYPOCOTYLS (LSH1; BnaC05g47190D), a phytochrome-dependent regulator of seedling development, in a region also associated with contrasting effects on flowering and height in Giessen in 2011 (early/small; late/high).

Central flowering regulators

Except for a homolog of FRI on chromosome A03 (BnaA03g13320D), none of the central flowering regulators (e.g., Bna.FT, Bna.CO or Bna.GI) or central vernalisation genes (e.g., Bna.FLC) were found among the candidate genes for the onset of flowering (BOF-ABS/ BOF-REL). Moreover, the A03 FRI copy was only detected in association with BOF-REL in the sub-tropical environment of Beibei 2012. A paralog of FLC was found in a region with suggestive association, also in Beibei 2012. When checking for closest marker associations of further copies of FLC, some copies (on chromosomes A01, A02 and C02) showed supporting, but no suggestive or significant associations in either BOF-REL or BOF-ABS. On the other hand, we found homologs of FUL (BnaA09g05500D, BnaC02g41870D), GI (BnaA09g30390D) and TEM1 (BnaA07g35130D) for HEI-ABS, TFL1 (BnaA10g26300D), FT (BnaC04g14850D) for YIE-ABS, and FT (BnaA07g25310D) for COMP.

Resequencing all copies of B. napus FT for a subset of 121 genotypes using a sequence capture approach (methods according to [34]) revealed non-synonymous SNP variation for BnaA07g25310D (R21Q, I49L) and BnaC04g14850D (H81Y). Whereas BnaA07g25310D exhibited heterozygous SNP calls, suggesting more than one homologous locus, clear genotyping calls were obtained for BnaC04g14850D. The T allele was associated with later flowering; taller plants and inferior yield compared to the C allele, with significant effects in different environments for flowering, plant height and yield (see Fig. 6).

Fig. 6

Boxplots showing yield values for genotypes carrying the FT-C-allele (blue) and the FT-T-allele (red) for each location and year. The environments are abbreviated as follows: Giessen (GI), Gross Gerau (GG), Rauischholzhausen (RH) with their respective year of harvest. Significant differences tested with the student’s t-test are labelled (*p < 0.05, **p < 0.01, ***p > 0.001)

Resequencing across all copies of Bna.CO revealed 10 nonsynonymous SNPs for BnaA03g27860D (4), BnaA10g18430D (1), BnaC03g32910D (1) and BnaAnng39160D (4). Testing those SNPs for association to flowering time, plant height and yield revealed associations for SNPs on BnaA10g18430D (E203D, associated with flowering time and yield), BnaC03g32910D (E79Q, associated with flowering time) and BnaAnng39160D (H165Q, associated with plant height, Q210H, associated with flowering time and yield). For Bna.SOC1, 4 nonsynonymous SNPs were detected for BnaC04g53290D (1), BnaC05g00840D (1) and BnaCnng45490D (2). Testing for associations revealed associations for SNPs on BnaC05g00840D (E191D, associated with plant height) and BnaCnng45490D (D122E, associated with all traits). We detected supporting associations only for SNPs closest to BnaA03g27860D (COMP, yield) and BnaC03g32910D (COMP, height), whereas no associations were detected for SNPs closest to BnaA10g18430D, BnaC09g46500D, BnaA04g26320D and BnaC05g00840D . LD between the adjacent flanking SNPs was low to very low (r 2 = 0.17 for BnaA03g27860D, 0.002 for BnaA10g18430D, 0.42 for BnaC03g32910D, 0.002 for BnaC09g46500D, and r 2 = 0.295 for BnaC04g53290D).

Allelic effects

Testing for allelic effects in BOF-ABS, HEI ABS and YIE-ABS, we observed significant additive effects for plant height and yield but not for flowering time. Whereas the phenotypic effect was subtle for plant height, it was strong for yield. Generally, the major allele induced an early, short, high-yielding phenotype, whereas the number of major alleles leading to late, high, low-yielding plants was small. This suggests that selection within this population has strongly driven phenotype towards early-flowering, smaller and better yielding plants (see Fig. 7, Additional file 4 Figure S2 for an example).

Fig. 7

Boxplots showing yield values for genotypes carrying the specified number of positive and negative alleles for Gross-Gerau in 2010

Gene ontology enrichment

GO enrichment analysis for all 1733 genes among the detected regions in BOF-ABS revealed significant enrichment of terms referring to post-transcriptional regulation of gene expression, for example RNA interference, miRNA regulation and RNA methylation. Terms related to epigenetic regulation were also significantly enriched, like methylation-dependent chromatin silencing, as well as dephosphorylation processes, proteasome core complex assembly and cell cycle (see Fig. 8). Vegetative phase change and response to stimulus were also among the significantly enriched GO terms. For BOF-REL, the most enriched terms among the 3133 genes were genes connected with response to oomycetes and negative regulation of programmed cell death. Petal development and aging were also enriched (Additional file 5: Figure S3).

Fig. 8

Gene ontology scatterplot constructed with REVIGO in R [39] for all genes in the associated regions for absolute values of flowering time (BOF-ABS). The scatterplot accounts for all GO terms with a count of min. 10., and were constructed using the Arabidopsis thaliana GO term database, using a similarity of 0.5, with SimRel as semantic similarity measure. Colours indicate the p-value of enrichment according to the legend. The size of each bubble reflects the count of each term among the enriched term list. Terms with a –log(p) > 1.5 are ascribed to their bubbles

GO enrichment analysis carried out for 2249 genes among regions of interest for HEI-ABS revealed the importance of small molecule metabolic processes for plant height, including monocarboxylic acid metabolism, cellular ketone metabolic process, carbohydrate metabolic process, alcohol catabolic process and the electron transport chain. The second group of enriched GO terms included processes involved in developmental regulation, like lateral root development, cell fate commitment and leaf senescence. We found also terms referring to virus and insect responses, to regulation of protein metabolic processes, as well as the cellular biogenic amine metabolism (Additional file 5: Figure S3).

The GO enrichment performed for the 2144 genes among the interesting regions for YIE-ABS showed a strong enrichment for metabolic pathways like triterpenoid metabolism, steroid metabolism or regulation of lipid metabolic processes. Pathways like sulfur, phosphorus and nitrogen compound metabolism were also enriched, as well as positive regulation of developmental processes, cell proliferation and gene expression (Additional file 5: Figure S3).

For 2162 genes detected in COMP, we found a strong enrichment for response to environmental factors like abiotic, biotic or endogenous stimuli, hyperosmotic response or response to heat. Moreover there was enrichment for biosynthesis of both small and large molecules, (e.g., flavonoids), immune system processes and macromolecule methylation (Additional file 5: Figure S3).


The diversity set used to assemble our winter-type B. napus panel was designed to represent the available breeding gene pool in the species [35, 36]. Genotyping the selected 158 winter lines with the Brassica 60 k Illumina® Infinium SNP array allows deep insight into the genetic variation present in this subset, although a large part of genotyping information is lost due to ambiguity in the reference genome. Chequerboard patterns in LD (Fig. 5) also indicate that there still may be issues concerning the reference assembly, as markers linked to each other are separated by ones which are not linked to them. Another explanation for this phenomenon could be rearrangements in the genomes of some lines, a common phenomenon in B. napus [20], which could decrease local LD. Since almost all accessions were adapted winter-type materials, the relatedness within the selected population was medium to high, reflecting the low genetic diversity expected within this gene pool [37]. The narrow genetic diversity was reflected by the relatively low associations detected for all traits. The low effects of associated loci also reflect the nature of the traits under investigation, which are influenced by a large number of different small-effect loci, reducing single contributions. The investigation of allelic effects for absolute values of flowering time, plant height and yield revealed that the population was subjected to strong selection towards an early flowering, small growing and high yielding phenotype. The results suggest that there is still potential to increase yield by pyramiding positive alleles. This might also be possible with regard to plant height, but does not hold true for flowering time. Flowering time genes do not act in parallel, but rather in a complex adaptive network, therefore effects are not stable over different environments. The suggestion that pyramiding of haplotypes may have a positive effect on yield is mirrored by observations of large homoeologous regions associated to yield (see below), which also points to parallel selection. All the same, the number of positive alleles is limited, demonstrating the importance of increasing diversity to ensure stable yield progress in changing environments [38].

Due to the low association levels, we used a significance cutoff based on FDR, to avoid loss of interesting loci, and used different approaches for detection of important regulators. Calculating genome-wide associations with absolute values for the onset of flowering, plant height and yield allowed us to find regions responsible for differences between genotypes (BOF-ABS, HEI-ABS, YIE-ABS), whereas calculating with the difference of the respective absolute value from the genotypic mean allowed us to detect regions responsible for differential population behaviour in different environments (BOF-REL). Comparisons between such relative association analyses for the same set of environments allowed detection of loci with adaptive potential (COMP).

The GO term enrichment analyses enabled us to identify the most important biological processes influencing the traits under study. Although the function of candidate genes cannot be confirmed without further reverse genetic analyses, this method allows dissection of trait regulation into pathways even without knowing the exact genetic players [39]. For BOF-ABS, the GO enrichment analysis revealed a superordinate influence of post-transcriptional gene expression regulation on floral initiation. This underlines the particular importance of epigenetic processes like methylation and RNA interference in determination of flowering time. Epigenetic factors play a huge role in adaptive traits, as for example response to abiotic or biotic stress [40]. This is consistent with our results, as non-optimal conditions like pests, water/nutrient deficiency or temperature stress are expected to play a significantly higher role in field trials than in the controlled conditions under which experiments for gene function are performed. This may in part explain why our analysis of flowering time detected only a few classical flowering time genes that do not explain all detected associations. Another reason may lie in the population, which might have been strongly selected for certain alleles. All the same, detection of some of those genes gives hints to pathways that play a role in flowering control under field variable conditions. The A. thaliana homologs of the detected candidates Bna SEP3, Bna SEP4, Bna FUL, Bna.AGL42, Bna.AGL71 and Bna AG are all regulated directly or indirectly by SOC1 [41, 42], which in turn is regulated by NUCLEAR FACTOR Y. We therefore propose that the SOC1 regulatory pathway has a major impact on flowering time determination of winter-type B.napus under field conditions (Fig. 9). Sequencing for all Bna.SOC1 copies in a subset of 121 genotypes revealed a low functional diversity. This indicates Bna.SOC1 may be highly conserved due to its central role. The only localized Bna.SOC1 copy that shows functional variation (BnaC04g53290D) lay in a region which was not associated to any trait at all. Two copies which could not be assigned to a specific location due to the ambiguity in the reference assembly carried functional variation which was associated to the phenotypic variation.

Fig. 9

Schematic representation of interactions in the Arabidopsis SOC1 regulation module for all respective candidate genes. The scheme includes all interaction partners which were detected as candidate genes for the onset of flowering in this study. For clarity, the FT copy detected in YIE-ABS was also added. Simple arrows represent regulation, double arrows represent interaction

Moreover, a number of candidate genes detected in BOF-REL, like Bna.FWA, Bna.TOE2 and Bna.GRP3 [43, 44], may relate to epigenetic control of flowering. Copies of TFL1 have recently been found to delay flowering time and increase seed yield in B.napus [45].

The GO enrichment for BOF-REL revealed a prevalence of developmental GO terms like negative regulation of programmed cell death, petal or xylem development and aging, indicating that variation at flowering loci impacts many genes downstream of the known regulatory cascades. References to small RNA regulation coincide with the detection of similar AGL transcription factors from the SOC1 regulation module (Bna.SEP4, Bna.AGL42, Bna.AGL71, Bna.LEUNIG) to those identified for BOF-ABS. Interestingly, we also detected two adjacent copies of Bna.AGL1/SHP1, whose A. thaliana homolog was found to be an antagonist of FUL in fruit dehiscence [46], indicating that that this gene might also interact with FUL in flowering (Fig. 9).

The only basic flowering time gene associated with BOF-REL was Bna.FRI, a well-described vernalisation regulator, in the Chinese environment. Bna.FRI.a, corresponding to BnaA03g13320, has already been established as a major flowering time determinator in B.napus [47]. Ath.FRI was also found before to play a role in climatic adaptation [1, 48]. As vernalisation is a quantitative process [49], we hypothesize that, within gene pools adapted to cold winter conditions, functional variation in central vernalisation genes may only play a role when temperatures are barely inductive for vernalisation, but not in environments with long periods of temperatures around and below freezing. Natural variation in vernalisation requirement was previously observed in Arabidopsis [50]. Similarly, we were also unable to detect significant associations with central integrators of the photoperiod pathway like CO or GI.

GO enrichment analysis for HEI-ABS and YIE-ABS revealed the primary importance of metabolic efficiency for plant height and yield. Whereas small molecule catabolism and the electron transport chain were implicated in control of plant height, the metabolism of higher-order molecules like triterpenoids or lipids appeared to be more influential in regulation of seed yield. Both traits enriched GO terms related to nitrogen metabolism, whereas phosphorus and sulfur metabolism were only enriched for yield, putatively in the context of lipid and glucosinolate metabolism. This may reflect the generally more advanced breeding level and higher yields of modern oilseed rape forms which also have high oil and low seed glucosinolate contents.

Although HEI-ABS and YIE-ABS give a similar GO enrichment, the overlap of associated regions was the lowest between all traits, while the highest overlap was observed for BOF-ABS and HEI-ABS. This is consistent with the respective phenotypic correlations, where we could not see a correlation trend between plant height and yield, but a significantly positive correlation between the onset of flowering and plant height was observed. This indicates that the switch to the reproductive phase determined plant height, however plant height was neither decisive nor indicative for seed yield. All the same, there was a negative correlation between the onset of flowering and yield, supporting the assumption that flowering time influences both traits independently of each other. Notably, some of the overlapping regions contained genes from the SOC1 regulation module (Bna.FUL, Bna.FYF).

It was striking that for YIE-ABS, more than 30 % of the candidates were detected with more than one copy, caused by several regions of interest being homoeologous, indicating co-expression of those genes between the subgenomes. The highest and largest regions of interest were found on homologous regions of chromosomes A09 and C09. The region on A09 does not colocalize with QTL for seed weight detected by [28], and nor did a further region on A07, probably because [28] used a population spanning winter-type, semi-winter and spring-type B.napus representing highly divergent gene pools. All the same, the high subgenome homology between associated regions shows that seed yield in B.napus may indeed benefit from fixed heterosis among homoeologous alleles, as proposed by Basunanda and co-workers [51]. This is strengthened by the pyramiding allelic effects for yield (Fig. 7).

Despite the low overlap of candidate regions for yield with flowering and plant height, we detected one region on A09 which was conserved for all traits, containing AGL transcription factor Bna.FUL and histone methylation factor Bna.GRP3. While the former was detected in 2 copies for BOF-REL, the latter was detected in 3 copies for YIE-ABS. FUL was found to regulate not only flowering meristems, but all meristems [52] and has been shown to be regulated by NTT, leading to altered patterning of the silique [53].

Although not annotated as such, we also consider some flowering and clock genes found among the associated regions to be possible adaptation loci influencing yield. These include two copies each of Bna.FUL, Bna.FPF1 and Bna.PRR5, and one copy each of Bna.TFL1, Bna.FT and Bna.PRR1/TOC1. Some of those genes have already been found to have an impact on seed yield: Ath.FPF1 was proposed to have a further role in shade avoidance [54]. Ath.TOC1 and Ath.PRR5 are components of the circadian clock, which has been shown to influence overall plant fitness in allopolyploids [55]. Bna.TFL1 and Bna.FT have already been shown to influence seed yield [45], whereas this was also the case for Ath.FUL [53]. Further evidence is provided by sequence data for all copies of Bna.FT in a subset of 121genotypes (Fig. 6), showing a significant yield reduction for the minor allele in two environments (Giessen and Gross-Gerau 2012), whereas the trend was similar when all environments were considered.

To find further candidates for adaptation, we performed a comparative analysis for relative values in all German environments. An additional copy of Bna.FUL was implicated in changes in all three traits, further suggesting that this gene may represent a general adaptation factor. Other genes which could be confirmed by COMP were Bna.SPL1 and Bna.SPL4.

Promising candidates for adaption genes are the clock components Bna.CCA1 and Bna.TIC. Both were detected in Gross-Gerau in 2011 in a region inversely affecting height and yield. As discussed before, a general influence of the circadian clock on plant fitness (e.g., through regulation of chlorophyll content and starch metabolism) was already shown for allopolyploids [55]. It was also interesting to note that a copy of FT (BnaA07g25310D) was associated with both plant height and yield in Rauischholzhausen in 2011. This corresponds to the findings for YIE-ABS, where we detected a further copy of FT (BnaC04g14850D) in a region associating with yield, underlined by independent sequence data, and with previous data suggesting that B. napus FT influences seed yield [45].

Interestingly, sequencing also revealed considerable functional variance for Bna.CO. Three SNPs obtained from resequencing a subset of 121 lines for CO copies were associated with flowering time. Bna.CO are strong candidates for flowering time determination [56, 57] and have been considered as a candidates in many QTL studies [11, 2426], but were not detected in any region associated to flowering time in the present study. All the same, the very low LD between the flanking markers of two copies onA10 and C09 in particular suggests that strong selection might have broken down LD in this area. This highlights the importance of an improved reference sequence assembly, which will then allow anchoring of all important candidate genes to associated regions. Meanwhile, the presented data represent a valuable base for genetic studies and may assist in developing new breeding strategies for this important oil crop.


We identified 101 B.napus genome regions associating with the onset of flowering, 69 with plant height and 36 with seed yield in highly diverse field environments. A total of 68 cross-trait regions were delineated. Generally, the pattern of associating regions was very environment specific, especially for flowering time, indicating a high influence of local field conditions. Several potential candidate genes for plant adaptation were found in chromosome regions associating with two or more traits. These include components of the circadian clock like Bna.CCA1 and Bna.TIC, the flowering regulators Bna.FT and Bna.FUL, and developmental factors like Bna.SPL4. Using GO enrichment we found that trait-associated regions are enriched for genes involved in post-transcriptional regulation of gene expression, indicating that processes like RNA interference represent an important regulatory mechanism for flowering time determination in winter-type B.napus under field conditions. Whereas homologs of the basic flowering regulators Bna.CO and Bna.FLC where not detected in regions associated with flowering time, we found orthologs to several genes interacting either directly or indirectly with SOC1, an important signal integrator in the flowering pathway. Interestingly, AGL transcription factors were particularly strongly represented. Several flowering time regulators, including an homolog of TFL1, were detected in regions associated with seed yield. Our data suggest that numerous genes directly or indirectly involved in flowering regulation play pleiotropic roles in adaptation to environmental stresses, and thus represent important genetic factors to be considered for adaptive breeding. The interaction of post-transcriptional and epigenetic networks with adaptation genes and stress networks in contrast to classical transcription factors represents an important area for future research in this regard,. Regarding yield, pyramiding positive haplotypes detected in this study might assist in increasing yield further, but introgression of new genetic diversity into the winter oilseed rape breeding pool is necessary to stabilize yield gains in the face of environmental change.


Plant Material

The plant material used in the present study was selected from the ERANET-ASSYST B.napus diversity set, a panel of over 500 homozygous inbred lines (S5 or greater) described in [35, 36]. A panel of 158 adapted winter-type accessions were grown in different locations in Germany (Rauischholzhausen, Giessen, Gross Gerau) from 2009 until 2013, in southwest China (Beibei, Chongqing) from 2012 to 2013, and in central Chile (Temuco) in 2011 and 2012, as described in [34]. In all cases the plants were grown in autumn-sown trials and harvested the following spring or summer. In Germany, where winters generally have prolonged periods with temperatures well below freezing, the trials were sown in late August or early September, with harvest the following July. In southern-hemisphere Temuco, where the winter is comparatively short and temperatures are seldom below zero, the trials were sown in mid-April and harvested in January, while in the sub-tropical climate of Beibei winter-type rapeseed must be planted late September or early October to prevent flowering before winter. Although temperatures are then sufficiently cold to vernalise winter rapeseed, the cold period is very short so that flowering occurs much quicker than under northern European conditions.

Field trials

Plots for phenotyping and DNA extraction were sown in a completely randomized block design with a harvest plot size of 3 x 1.25 m in 1 replicate (containing around 200 plants) except for Gross Gerau 2012, where 2 replicates were sown. Positional effects were controlled using an early flowering and a late flowering control genotype in 4 replicates, including them in the block design. When both genotypes were significantly separated using Student’s t-test on a p < 0.05 base, we concluded that the spatial impact was negligible. This was the case for all locations studied with the exception of GI2011 (p = 0.251). Single replication trials are standard procedure for determination of flowering time due to the high heritability. Plant protection and fertilisation was performed according to normal practice at the respective locations, as necessary for the prevailing conditions.

DNA isolation

Leaf material for genomic DNA extraction was harvested from each accession in spring 2012 from the field trial in Giessen, Germany. Pooled leaf samples were taken from at least 5 different plants per genotype, immediately shock-frozen in liquid nitrogen and kept at −20 °C until extraction. Leaf material was ground in liquid nitrogen with a mortar and pestle. DNA was extracted using a common CTAB protocol modified from Doyle and Doyle (1990) as described earlier [34]. DNA concentration was determined using a Qubit fluorometer and the Qubit dsDNA BR assay kit (Life Technologies, Darmstadt, Germany) according to the manufacturer’s protocol. DNA quantity and purity was further checked on 0.5 % agarose gel (3 V/cm, 0.5xTBE, 120 min).


Flowering times were recorded for 158 accessions over a total of 11 environments from 2010–2012 at the 5 locations Giessen, Germany (50° 35′ N, 8° 40′ E), Gross Gerau, Germany (49° 55′ N, 8° 29′ E), Rauischholzhausen, Germany (50° 46′ N, 8° 53′ E), Beibei, Chongqing, China (29° 33′ N, 106° 34′ E) and Temuco, Chile (38° 44′ S, 72° 36′ W)). The onset of flowering was defined as the day where 10 % of the plants in a plot (3.75 m2) reached BBCH 61. Plant height in cm was measured at the onset of flowering. Plots were harvested in July (Germany) and January (Chile) and seed yield (corrected to dt/ha) was determined.

SNP genotyping and data preprocessing

The 158 accessions were genotyped using the Brassica 60 K Illumina® Infinium SNP array by TraitGenetics GmbH (Gatersleben, Germany). Flanking sequences of 52,160 called SNPs were localized by BLAST to the B.napus Darmor-bzh reference genome assembly v4.1 [20], allowing no less than 50 bp overlap, 95 % sequence identity, and no gaps. SNPs with ambiguous position and SNPs on chr_random were excluded from the analysis, leaving 28,698 unique, locus-specific SNPs with a defined genome position. Heterozygous calls were treated as missing values. After preprocessing the marker set for non-missing marker values >0.9, minor allele frequency > 0.01 and non-missing individual markers >0.8 we retained 21,623 unique SNP markers and a population of 140 individuals for genome-wide association analysis. This resulted in a mean marker density of 1138 markers per chromosome, covering approximately 32.8 kbp/marker. Data preprocessing was performed with R (version 3.1.0) using the package GenABEL [58].

Association analysis

Population structure analysis and visualization were performed in R (version 3.1.0) using the packages SelectionTools (, ape and scatterplot3d that apply Principal Component Analysis based on genetic distances calculated according to the euclidean modified Rodger’s distance method. The most likely number of population subclusters was estimated by plotting the within-cluster sum of squares against the possible number of clusters, ranging from 1 to 15 (see also Additional file 2: Figure S1). K-means clustering was then performed in R using SelectionTools. Genome-wide association testing and calculation of genome-wide allele identity-by-descent were performed with GenABEL. F ST -values were calculated with the program Genepop (version 4.2.2, [59]. Linkage disequilibrium (LD) was calculated in terms of coefficient of determination (r 2) in R, separately for each chromosome. Decay of LD was displayed using smoothed scatterplot curves with a span of 0.1.

Association analysis for flowering time was carried out using a mixed model with PC-adjustment using 21,623 unique SNPs for 11 environments, of which one was replicated. Using Bonferroni correction, only associations with p < 4.6*10-6 (giving a –log10(p) = 5.33) would prove significant on a p < 0.1 level. Predominantly weak associations are expected for highly quantitative traits like flowering time within an adapted gene pool. Because only 3 of the tested environments had SNPs meeting this condition, associations were considered as significant if –log10(p) > 4, as suggestive if –log10(p) > 3 and as supporting if –log10(p) > 2, based on FDR cutoff calculations performed with the R package fdrtool, proposing a cutoff of –log10(p) > 2. SNPs were considered for further analysis if either (1) a significant association was detected in at least one environment, or (2) at least one suggestive and one supporting association were detected in respective environments. Regions of interest were defined using the LD decay around the most significantly associated SNP markers, extending until r2 decayed to less than 0.5, meaning that regions extended from the left unlinked marker to the right unlinked marker. For plant height and yield, the analysis was performed likewise, with 10 and 8 environments, respectively. In a second approach, phenotypes were calculated as deviation from mean for every genotype in every location (see below), and association analysis was performed likewise.

$$ BOF{(rel)}_{envx}= mean\left(BOF(abs)\right)-BOF{(abs)}_{envx} $$

For comparative analysis in German environments, SNPs were considered as associated if they were at least suggestive for one trait and supporting for another.

Association analysis for the beginning of flowering (BOF) was performed in two ways: once using the absolute flowering time measured as days after sowing (considered as BOF-ABS), and once calculated as BOF-REL. This on the one hand allows us to find loci responsible for early or late flowering in the respective environment (BOF-ABS), and on the other hand to detect loci which are responsible for differential flowering behaviour between environments (BOF-REL). For plant height and yield, association analysis was performed exclusively using absolute values (HEI-ABS, YIE-ABS). Candidate genes were selected based on GO terms, considering the terms listed below.

flower photosynthesis lipid
vernalization mitosis terpenoid
photoperiod mitotic biosynthetic process
circadian cell growth catabolic process
floral development gluconeogenesis
vegetative to reproductive biosynth glycolysis
vegetative phase change catabolic amino acid transport
pollen metabolic embryo development
carpel tricarboxylic acid cycle embryo sac development
sepal glycolysis  

Furthermore, we carried out a comparative analysis between all three traits for German environments (COMP). For this, relative values were calculated only considering German environments, to find regions important for differential behaviour of two or more traits, indicating adaptive potential.

Allelic effects

To assess the allelic effects for all associated regions for BOF-ABS, HEI-ABS and YIE-ABS, we determined the allelic state for each highest associating marker and genotype. We assigned a value to each marker, “0” for the reference allele and the value given by association analysis for the alternative allele. We counted the number of positive/early and negative/late alleles and tested for allelic effects using least significant distance tests in R using the package “agricolae”.

GO enrichment analysis

GO terms of all genes detected among the defined regions of interest were extracted in R, using the package GOstats ( ), separately for the tags BP, MF and CC. Each GO term has a unique identifier with a defined frequency in the Brassica napus genome annotations [20]. The enrichment of each GO term is then determined using a hypergeometrical test using a p-value cut-off of 0.05. Results were filtered for GO terms with a count of at least 10, and visualized using the web-based program REVIGO ([39]), using the following settings: database: Arabidopsis thaliana, similarity: 0.5 (small), semantic similarity measure: SimRel.

Availability of supporting data

SNP data are deposited in dbsnp ( and are available as Supporting data (Additional file 6: Table S3).





Biologische Bundesanstalt, Bundessortenamt und CHemische Industrie


Mode of analysis: flowering time, absolute values


Mode of analysis: flowering time, relative values


Biological process


Base pairs


Cellular component

CCA1 :

Circadian clock- associated 1


Mode of analysis: comparative

CRY2 :

Cryptochrome 2




Days after sowing


Desoxyribonucleic acid


False discovery rate


Flowering locus C

FPF1 :

Flowering Promoting Factor 1





FT :

Flowering locus T




Flowering Wageningen


Forever young flower

GA2OX2 :

Gibberellin 2-oxidase 2

GA20OX1 :

Gibberellin 20-oxidase 1

GI :


GRP3 :

Glycine-rich RNA-binding protein 3


Gene ontology


Genome-wide association study


Mode of analysis: plant height, absolute values


Identity by descent


Linkage disequilibrium

LSH1 :

Light-dependent short hypocotyls


Molecular function




Non-coding RNA


Principal Component Analysis


Pseudo response regulator


Phytochrome C


Quantitative trait locus

RLP19 :

Receptor like protein 19


Ribonucleic acid


5th selfing generation


Single nucleotide polymorphism

SOC1 :

Suppressor of Overexpression of Co 1


Squamosa-promotor binding-like proteins

SRR1 :

Sensitivity to red light reduced 1


Tris-borate - EDTA




Time for coffee

TOE2 :

Target of eat 2

VIN3 :

Vernalisation insensitive 3


Mode of analysis: seed yield, absolute values




  1. 1.

    Mendez-Vigo B, Pico FX, Ramiro M, Martinez-Zapater JM, Alonso-Blanco C. Altitudinal and Climatic Adaptation Is Mediated by Flowering Traits and FRI, FLC, and PHYC Genes in Arabidopsis. Plant Physiology. 2011;157:1942–55.

    CAS  Article  PubMed Central  PubMed  Google Scholar 

  2. 2.

    Jaeger KE, Graf A, Wigge PA. The control of flowering in time and space. J Exp Bot. 2006;57:3415–8.

    CAS  Article  PubMed  Google Scholar 

  3. 3.

    Jung C, Muller AE. Flowering time control and applications in plant breeding. Trends Plant Sci. 2009;14:563–73.

    CAS  Article  PubMed  Google Scholar 

  4. 4.

    Srikanth A, Schmid M. Regulation of flowering time: all roads lead to Rome. Cell Mol Life Sci. 2011;68:2013–37.

    CAS  Article  PubMed  Google Scholar 

  5. 5.

    Wigge PA. Ambient temperature signalling in plants. Curr Opin Plant Biol. 2013;16:661–6.

    CAS  Article  PubMed  Google Scholar 

  6. 6.

    Pak H, Guo Y, Chen M, Chen K, Li Y, Hua S, et al. The effect of exogenous methyl jasmonate on the flowering time, floral organ morphology, and transcript levels of a group of genes implicated in the development of oilseed rape flowers (Brassica napus L.). Planta. 2009;231:79–91.

    CAS  Article  PubMed  Google Scholar 

  7. 7.

    de Montaigu A, Toth R, Coupland G. Plant development goes like clockwork. Trends Genet. 2010;26:296–306.

    Article  PubMed  Google Scholar 

  8. 8.

    Khan MRG, Ai X, Zhang J. Genetic regulation of flowering time in annual and perennial plants. WIREs RNA. 2014;5:347–59.

    CAS  Article  PubMed  Google Scholar 

  9. 9.

    Li Y, Huang Y, Bergelson J, Nordborg M, Borevitz JO. Association mapping of local climate-sensitive quantitative trait loci in Arabidopsis thaliana. Proc Natl Acad Sci. 2010;107:21199–204.

    CAS  Article  PubMed Central  PubMed  Google Scholar 

  10. 10.

    Dechaine JM, Brock TM, Weinig C. QTL architecture of reproductive fitness characters in Brassica rapa. BMC Plant Biol. 2014;14:66.

    Article  PubMed Central  PubMed  Google Scholar 

  11. 11.

    Quijada PA, Udall JA, Lambert B, Osborn TC. Quantitative trait analysis of seed yield and other complex traits in hybrid spring rapeseed (Brassica napus L.): 1. Identification of genomic regions from winter germplasm. Theor Appl Genet. 2006;113:549–61.

    CAS  Article  PubMed  Google Scholar 

  12. 12.

    Udall JA, Quijada PA, Lambert B, Osborn TC. 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. 2006;113:597–609.

    CAS  Article  PubMed  Google Scholar 

  13. 13.

    Shi J, Li R, Qiu D, Jiang C, Long Y, Morgan C, et al. Unraveling the Complex Trait of Crop Yield With Quantitative Trait Loci Mapping in Brassica napus. Genetics. 2009;182:851–61.

    CAS  Article  PubMed Central  PubMed  Google Scholar 

  14. 14.

    Deng W, Ying H, Helliwell CA, Taylor JM, Peacock WJ, Dennis ES. FLOWERING LOCUS C (FLC) regulates development pathways throughout the life cycle of Arabidopsis. Proc Natl Acad Sci. 2011;108:6680–5.

    CAS  Article  PubMed Central  PubMed  Google Scholar 

  15. 15.

    Friedt W, Snowdon R. J.: Oilseed rape. In Handbook of Plant Breeding, Vol. 4: Oil crops breeding. Edited by Vollmann J IR. Berlin/ Heidelberg; 2010.

  16. 16.

    Paterson AH, Lan T, Amasino R, Osborn TC, Quiros C. Brassica genomics: a complement to, and early beneficiary of, the Arabidopsis sequence. Genome Biol. 2001;2(3):REVIEWS1011.

    CAS  Article  PubMed Central  PubMed  Google Scholar 

  17. 17.

    Pires JC, Zhao JW, Schranz ME, Leon EJ, Quijada PA, Lukens LN, et al. Flowering time divergence and genomic rearrangements in resynthesized Brassica polyploids (Brassicaceae). Biol J Linn Soc. 2004;82:675–88.

    Article  Google Scholar 

  18. 18.

    Rana D, Boogaart T, O’Neill CM, Hynes L, Bent E, Macpherson L, et al. Conservation of the microstructure of genome segments in Brassica napus and its diploid relatives. Plant J. 2004;40:725–33.

    CAS  Article  PubMed  Google Scholar 

  19. 19.

    Lysak MA. Chromosome triplication found across the tribe Brassiceae. Genome Res. 2005;15:516–25.

    CAS  Article  PubMed Central  PubMed  Google Scholar 

  20. 20.

    Chalhoub B, Denoeud F, Liu S, Parkin IAP, Tang H, Wang X, et al. Early allopolyploid evolution in the post-Neolithic Brassica napus oilseed genome. Science. 2014;345:950–3.

    CAS  Article  PubMed  Google Scholar 

  21. 21.

    Chen W, Zhang Y, Liu X, Chen B, Tu J, Tingdong F. Detection of QTL for six yield-related traits in oilseed rape (Brassica napus) using DH and immortalized F2 populations. Theor Appl Genet. 2007;115:849–58.

    CAS  Article  PubMed  Google Scholar 

  22. 22.

    Raman H, Raman R, Eckermann P, Coombes N, Manoli S, Zou X, et al. Genetic and physical mapping of flowering time loci in canola (Brassica napus L.). Theor Appl Genet. 2013;126:119–32.

    CAS  Article  PubMed  Google Scholar 

  23. 23.

    Li X, Ramchiary N, Dhandapani V, Choi SR, Hur Y, Nou I, et al. Quantitative Trait Loci Mapping in Brassica rapa Revealed the Structural and Functional Conservation of Genetic Loci Governing Morphological and Yield Component Traits in the A, B, and C Subgenomes of Brassica Species. DNA Res. 2013;20:1–16.

    Article  PubMed Central  PubMed  Google Scholar 

  24. 24.

    Fletcher RS, Mullen JL, Heiliger A, McKay JK. QTL analysis of root morphology, flowering time, and yield reveals trade-offs in response to drought in Brassica napus. J Exp Bot. 2014;66(1):245–56.

    Article  PubMed Central  PubMed  Google Scholar 

  25. 25.

    Luo YX, Luo CY, Du DZ, Fu Z, Yao YM, Xu CC, et al. Quantitative trait analysis of flowering time in spring rapeseed (B. napus L.). Euphytica. 2014;200:321–35.

    Article  Google Scholar 

  26. 26.

    Nelson MN, Rajasekaran R, Smith A, Chen S, Beeck CP, Siddique KHM, et al. Quantitative Trait Loci for Thermal Time to Flowering and Photoperiod Responsiveness Discovered in Summer Annual-Type Brassica napus L. PLoS One. 2014;9:e102611.

    Article  PubMed Central  PubMed  Google Scholar 

  27. 27.

    Cai D, Xiao Y, Yang W, Ye W, Wang B, Younas M, et al. Association mapping of six yield-related traits in rapeseed (Brassica napus L.). Theor Appl Genet. 2014;127:85–96.

    CAS  Article  PubMed  Google Scholar 

  28. 28.

    Li F, Chen B, Xu K, Wu J, Song W, Bancroft I, et al. Genome-Wide Association Study Dissects the Genetic Architecture of Seed Weight and Seed Quality in Rapeseed (Brassica napus L.). DNA Res. 2014;21(4):355–67.

    Article  PubMed Central  PubMed  Google Scholar 

  29. 29.

    Hatzig SV, Frisch M, Breuer F, Nesi N, Ducournau S, Wagner M, et al. Genome-wide association mapping unravels the genetic control of seed germination and vigour in Brassica napus. Front Plant Sci. 2015;6:221.

    Article  PubMed Central  PubMed  Google Scholar 

  30. 30.

    Tadege M, Sheldon C, Helliwell CA, Stoutjesdijk P, Dennis ES, Peacock WJ. Control of flowering time by FLC orthologues in Brassica napus. Plant J. 2001;28(5):545–53.

    CAS  Article  PubMed  Google Scholar 

  31. 31.

    Dreni L, Kater MM. MADS reloaded: evolution of the AGAMOUS subfamily genes. New Phytol. 2014;201:717–32.

    CAS  Article  PubMed  Google Scholar 

  32. 32.

    Yoo SK, Wu X, Lee JS, Ahn JH. AGAMOUS-LIKE 6 is a floral promoter that negatively regulates the FLC/MAF clade genes and positively regulates FT in Arabidopsis. Plant J. 2011;65:62–76.

    CAS  Article  PubMed  Google Scholar 

  33. 33.

    Adamczyk BJ, Fernandez DE. MIKC* MADS domain heterodimers are required for pollen maturation and tube growth in Arabidopsis. Plant Physiology. 2009;149:1713–23.

    CAS  Article  PubMed Central  PubMed  Google Scholar 

  34. 34.

    Schiessl S, Samans B, Hüttel B, Reinhardt R, Snowdon RJ. Capturing sequence variation among flowering-time regulatory gene homologs in the allopolyploid crop species Brassica napus. Front Plant Sci. 2014;5:404.

    Article  PubMed Central  PubMed  Google Scholar 

  35. 35.

    Bus A, Hecht J, Huettel B, Reinhardt R, Stich B. High-throughput polymorphism detection and genotyping in Brassica napus using next-generation RAD sequencing. BMC Genomics. 2012;13:281.

    CAS  Article  PubMed Central  PubMed  Google Scholar 

  36. 36.

    Körber N, Wittkop B, Bus A, Friedt W, Snowdon RJ, Stich B. Seedling development in a Brassica napus diversity set and its relationship to agronomic performance. Theor Appl Genet. 2012;125:1275–87.

    Article  PubMed  Google Scholar 

  37. 37.

    Hasan M, Seyis F, Badani AG, Pons-Kühnemann J, Friedt W, Lühs W, et al. Analysis of Genetic Diversity in the Brassica napus L. Gene Pool Using SSR Markers. Genet Resour Crop Evol. 2006;53:793–802.

    CAS  Article  Google Scholar 

  38. 38.

    Snowdon RJ, Iniguez Luy FL. Potential to improve oilseed rape and canola breeding in the genomics era. Plant Breed. 2012;131:351–60.

    CAS  Article  Google Scholar 

  39. 39.

    Supek F, Bošnjak M, Škunca N, Šmuc T, Gibas C. REVIGO Summarizes and Visualizes Long Lists of Gene Ontology Terms. PLoS One. 2011;6:e21800.

    CAS  Article  PubMed Central  PubMed  Google Scholar 

  40. 40.

    Sahu PP, Pandey G, Sharma N, Puranik S, Muthamilarasan M, Prasad M. Epigenetic mechanisms of plant stress responses and adaptation. Plant Cell Rep. 2013;32:1151–9.

    CAS  Article  PubMed  Google Scholar 

  41. 41.

    Lee J, Lee I. Regulation and function of SOC1, a flowering pathway integrator. J. Exp. Bot. 2010;61:2247–54.

    CAS  Article  PubMed  Google Scholar 

  42. 42.

    Dorca-Fornell C, Gregis V, Grandi V, Coupland G, Colombo L, Kater MM. The Arabidopsis SOC1-like genes AGL42, AGL71 and AGL72 promote flowering in the shoot apical and axillary meristems. The Plant Journal. 2011;67:1006–17.

    CAS  Article  PubMed  Google Scholar 

  43. 43.

    Zhang B, Pan X, Cannon CH, Cobb GP, Anderson TA. Conservation and divergence of plant microRNA genes. The Plant Journal. 2006;46:243–59.

    CAS  Article  PubMed  Google Scholar 

  44. 44.

    Causier B, Ashworth M, Guo W, Davies B. The TOPLESS Interactome: A Framework for Gene Repression in Arabidopsis. Plant Physiology. 2012;158:423–38.

    CAS  Article  PubMed Central  PubMed  Google Scholar 

  45. 45.

    Guo Y, Hans H, Christian J, Molina C. Mutations in single FT- and TFL1-paralogs of rapeseed (Brassica napus L.) and their impact on flowering time and yield components. Front. Plant Sci. 2014;5:282

  46. 46.

    Pabon-Mora N, Wong GK, Ambrose BA. Evolution of fruit development genes in flowering plants. Front. Plant Sci. 2014;5:300.

  47. 47.

    Wang N, Qian W, Suppanz I, Wei L, Mao B, Long Y, et al. Flowering time variation in oilseed rape (Brassica napus L.) is associated with allelic variation in the FRIGIDA homologue BnaA.FRI.a. J. Exp. Bot. 2011;62:5641-58.

    CAS  Article  PubMed Central  PubMed  Google Scholar 

  48. 48.

    Choi K, Kim J, Hwang H, Kim S, Park C, Kim SY, et al. The FRIGIDA complex activates transcription of FLC, a strong flowering repressor in Arabidopsis, by recruiting chromatin modification factors. The Plant Cell. 2011;23:289-303.

    CAS  Article  PubMed Central  PubMed  Google Scholar 

  49. 49.

    Kim D, Doyle MR, Sung S, Amasino RM. Vernalization: winter and the timing of flowering in plants. Annu. Rev. Cell Dev. Biol. 2009;25:277-99.

    CAS  Article  PubMed  Google Scholar 

  50. 50.

    Wollenberg AC, Amasino RM. Natural variation in the temperature range permissive for vernalization in accessions of Arabidopsis thaliana. Plant, Cell & Environment. 2012;35:2181-91.

    Article  Google Scholar 

  51. 51.

    Basunanda P, Radoev M, Ecke W, Friedt W, Becker HC, Snowdon RJ: Comparative mapping of quantitative trait loci involved in heterosis for seedling and yield traits in oilseed rape (Brassica napus L.). Theor Appl Genet. 2010;120:271–81.

    CAS  Article  PubMed Central  PubMed  Google Scholar 

  52. 52.

    Melzer S, Lens F, Gennen J, Vanneste S, Rohde A, Beeckman T. Flowering-time genes modulate meristem determinacy and growth form in Arabidopsis thaliana. Nat Genet. 2008;40:1489–92.

    CAS  Article  PubMed  Google Scholar 

  53. 53.

    Chung KS, Lee JH, Lee JS, Ahn JH. Fruit indehiscence caused by enhanced expression of NO TRANSMITTING TRACT in Arabidopsis thaliana. Mol Cells. 2013;35:519–25.

    CAS  Article  PubMed Central  PubMed  Google Scholar 

  54. 54.

    Wang X, Fan S, Song M, Pang C, Wei H, Yu J, et al. Upland cotton gene GhFPF1 confers promotion of flowering time and shade-avoidance responses in Arabidopsis thaliana. PLoS ONE. 2014;9:e91869.

    Article  PubMed Central  PubMed  Google Scholar 

  55. 55.

    Ni Z, Kim E, Ha M, Lackey E, Liu J, Zhang Y, et al. Altered circadian rhythms regulate growth vigour in hybrids and allopolyploids. Nature. 2008;457:327–31.

    Article  PubMed Central  PubMed  Google Scholar 

  56. 56.

    Robert LS, Robson F, Sharpe A, Lydiate D, Coupland G: Conserved structure and function of the Arabidopsis flowering time gene CONSTANS in Brassica napus. Plant Molecular Biology. 1998;37(5):763–72.

    CAS  Article  PubMed  Google Scholar 

  57. 57.

    Osborn TC, Kole C, Parkin IAP, Sharpe AG, Kuiper M, Lydiate DJ, et al. Comparison of flowering time genes in Brassica rapa, B. napus and Arabidopsis thaliana. Genetics. 1997;146(3):1123-9.

    CAS  PubMed Central  PubMed  Google Scholar 

  58. 58.

    Aulchenko YS, Ripke S, Isaacs A, van Duijn CM. GenABEL: an R library for genome-wide association analysis. Bioinformatics. 2007;23:1294-6.

    CAS  Article  PubMed  Google Scholar 

  59. 59.

    Rousset F. genepop’007: a complete re-implementation of the genepop software for Windows and Linux. Mol Ecol Resour. 2008;8:103–6.

    Article  PubMed  Google Scholar 

Download references


The authors thank Mechthild Schwarte (Rauischholzhausen), Malte Lüh (Giessen), Mario Tolksdorf (Gross Gerau) and Fernando Westermeyer (Temuco) for valuable assistance with the field phenotyping, and Sarah Hatzig for bioinformatics support and helpful discussions.


Funding was provided by the German Research Foundation (DFG) grants 14/16-1 and 14/16-2 to RS.

Author information



Corresponding author

Correspondence to Sarah Schiessl.

Additional information

Competing interests

The authors declare to have no competing financial interests.

Authors’ contributions

SS coordinated the phenotyping, carried out the bioinformatic analysis and wrote the manuscript. WQ and FI planned and implemented the field trials in China and Chile, respectively. RS planned and implemented the field trials in Germany and revised the manuscript. All authors read and approved the final manuscript.

Authors’ information

Not applicable.

Availability of data and materials

Not applicable.

Additional files

Additional file 1: Table S1.

Pearson’s correlation coefficient (above diagonal) and p-values (below diagonal) for all phenotypic traits used in this study. Grey fields indicate correlations for the same environment. The environments are abbreviated as follows: Germany: Giessen (GI), Gross Gerau (GG), Rauischholzhausen (RH); Chile: Temuco (TE); China: Beibei, Chongqing (BB), with their respective year of harvest. (XLSX 26 kb)

Additional file 2: Figure S1.

Results of population structure analysis in the association panel of 158 winter-type accessions, displaying 3 clusters as determined by k-means clustering. (A) and (B) show two PCA plots (component 1 and 2; component 2 and 3, respectively). (C) The best number of clusters was estimated by plotting the within-sum of squares against the possible number of clusters. (D) shows the distribution of the genome-wide allele identity-by-descent (IBD) for the total set, whereas (F) to (H) show the same for each cluster. (E) shows the F ST -values between clusters (white fields) and the mean IBD values for each cluster (grey fields). (PNG 132 kb)

Additional file 3: Table S2.

Candidate genes for associating regions for BOF-ABS, BOF-REL, HEI-ABS, YIE-ABS and COMP. Each gene has a unique identifier according to [20] and a short sequence description. The table names the highest associating marker in the defined region in its LD block with r 2 > 0.5 as well as the SNP marker most closely linked to the respective gene with its r 2 value to the highest associating SNP marker, and the closest markers distance from the gene. It also gives the value of the highest association among the respective region. The next columns show the strength of association in the respective environment; * stands for –log10(p) > 2, ** for -log10(p) > 3 and *** for –log10(p) > 4, and the last columns show the absolute effect of the highest associating marker on the respective trait. (XLSX 553 kb)

Additional file 4: Figure S2.

Boxplots showing (A) flowering time values and (B) plant height for genotypes carrying the specified number of positive and negative alleles for Gross-Gerau in 2010. (ZIP 213 kb)

Additional file 5: Figure S3.

Gene ontology scatterplot constructed with REVIGO in R [35] for all genes in the associated regions for the following traits. (A) relative values of flowering time in respect to genotypic mean (BOF-REL) (B) absolute values for plant height (HEI-ABS) (C) absolute values for seed yield (YIE-ABS) (D) relative values for flowering time, plant height and seed yield (COMP). The treemaps account for all GO terms with a count of min. 10., and were constructed using the Arabidopsis thaliana GO term database, using a similarity of 0.5, with SimRel as semantic similarity measure. Colours indicate the p-value of enrichment according to the legend. The size of each bubble reflects the count of each term among the enriched term list. Terms with a –log(p) > 1.5 are ascribed to their bubbles. (ZIP 870 kb)

Additional file 6: Table S3.

SNP positions and allelic states for all 21,623 SNPs used to calculate associations. The positions refer to the reference sequence, version v4.1, published by [20], which can be found on (XLSX 851 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, 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 ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Schiessl, S., Iniguez-Luy, F., Qian, W. et al. Diverse regulatory factors associate with flowering time and yield responses in winter-type Brassica napus . BMC Genomics 16, 737 (2015).

Download citation


  • GWAS
  • Rapeseed
  • Population structure
  • Post-transcriptional regulation
  • SNP marker
  • Gene ontology enrichment