Genome‐wide association analyses reveal complex genetic architecture underlying natural variation for flowering time in canola

Optimum flowering time is the key to maximize canola production in order to meet global demand of vegetable oil, biodiesel and canola‐meal. We reveal extensive variation in flowering time across diverse genotypes of canola under field, glasshouse and controlled environmental conditions. We conduct a genome‐wide association study and identify 69 single nucleotide polymorphism (SNP) markers associated with flowering time, which are repeatedly detected across experiments. Several associated SNPs occur in clusters across the canola genome; seven of them were detected within 20 Kb regions of a priori candidate genes; FLOWERING LOCUS T, FRUITFUL, FLOWERING LOCUS C, CONSTANS, FRIGIDA, PHYTOCHROME B and an additional five SNPs were localized within 14 Kb of a previously identified quantitative trait loci for flowering time. Expression analyses showed that among FLC paralogs, BnFLC.A2 accounts for ~23% of natural variation in diverse accessions. Genome‐wide association analysis for FLC expression levels mapped not only BnFLC.C2 but also other loci that contribute to variation in FLC expression. In addition to revealing the complex genetic architecture of flowering time variation, we demonstrate that the identified SNPs can be modelled to predict flowering time in diverse canola germplasm accurately and hence are suitable for genomic selection of adaptative traits in canola improvement programmes.     

Prediction of flowering time in Brassica oleracea using a quantitative trait loci-based phenology model

Uniformly developing plants with a predictable time to harvest or flowering under unfavourable climate conditions are a major breeding goal in crop species. The main flowering regulators and their response to environmental signals have been identified in Arabidopsis thaliana and homologues of flowering genes have been mapped in many crop species. However, it remains unclear which genes determine within and across genotype flowering time variability in Brassica oleracea and how genetic flowering time regulation is influenced by environmental factors. The goal of this study is model-based prediction of flowering time in a B. oleracea DH-line population using genotype-specific and quantitative trait loci (QTL) model input parameters. A QTL-based phenology model accounting for genotypic differences in temperature responses during vernalisation and non-temperature-sensitive durations from floral transition to flowering was evaluated in two field trials. The model was parameterised using original genotype-specific model input parameters and QTL effects. The genotype-specific model parameterisation showed accurate predictability of flowering time if floral induction was promoted by low temperature (R(2) = 0.81); unfavourably high temperatures reduced predictability (R(2) = 0.65). Replacing original model input parameters by QTL effects reduced the capability of the model to describe across-genotype variability (R(2) = 0.59 and 0.50). Flowering time was highly correlated with a model parameter accounting for vernalisation effects. Within-genotype variability was significantly correlated with the same parameter if temperature during the inductive phase was high. We conclude that flowering time variability across genotypes was largely due to differences in vernalisation response, although it has been shown elsewhere that the candidate FLOWERING LOCUS C (FLC) did not co-segregate with flowering time in the same population. FLC independent vernalisation pathways have been described for several species, but not yet for B. oleracea.     

SNP markers‐based map construction and genome‐wide linkage analysis in Brassica napus

An Illumina Infinium array comprising 5306 single nucleotide polymorphism (SNP) markers was used to genotype 175 individuals of a doubled haploid population derived from a cross between Skipton and Ag‐Spectrum, two Australian cultivars of rapeseed (Brassica napus L.). A genetic linkage map based on 613 SNP and 228 non‐SNP (DArT, SSR, SRAP and candidate gene markers) covering 2514.8 cM was constructed and further utilized to identify loci associated with flowering time and resistance to blackleg, a disease caused by the fungus Leptosphaeria maculans. Comparison between genetic map positions of SNP markers and the sequenced Brassica rapa (A) and Brassica oleracea (C) genome scaffolds showed several genomic rearrangements in the B. napus genome. A major locus controlling resistance to L. maculans was identified at both seedling and adult plant stages on chromosome A07. QTL analyses revealed that up to 40.2% of genetic variation for flowering time was accounted for by loci having quantitative effects. Comparative mapping showed Arabidopsis and Brassica flowering genes such as Phytochrome A/D, Flowering Locus C and agamous‐Like MADS box gene AGL1 map within marker intervals associated with flowering time in a DH population from Skipton/Ag‐Spectrum. Genomic regions associated with flowering time and resistance to L. maculans had several SNP markers mapped within 10 cM. Our results suggest that SNP markers will be suitable for various applications such as trait introgression, comparative mapping and high‐resolution mapping of loci in B. napus.     

Molecular mapping of qualitative and quantitative loci for resistance to Leptosphaeria maculans causing blackleg disease in canola (Brassica napus L.)

Blackleg, caused by Leptosphaeria maculans, is one of the most important diseases of oilseed and vegetable crucifiers worldwide. The present study describes (1) the construction of a genetic linkage map, comprising 255 markers, based upon simple sequence repeats (SSR), sequence-related amplified polymorphism, sequence tagged sites, and EST-SSRs and (2) the localization of qualitative (race-specific) and quantitative (race non-specific) trait loci controlling blackleg resistance in a doubled-haploid population derived from the Australian canola (Brassica napus L.) cultivars Skipton and Ag-Spectrum using the whole-genome average interval mapping approach. Marker regression analyses revealed that at least 14 genomic regions with LOD ≥ 2.0 were associated with qualitative and quantitative blackleg resistance, explaining 4.6-88.9 % of genotypic variation. A major qualitative locus, designated RlmSkipton (Rlm4), was mapped on chromosome A7, within 0.8 cM of the SSR marker Xbrms075. Alignment of the molecular markers underlying this QTL region with the genome sequence data of B. rapa L. suggests that RlmSkipton is located approximately 80 kb from the Xbrms075 locus. Molecular marker-RlmSkipton linkage was further validated in an F(2) population from Skipton/Ag-Spectrum. Our results show that SSR markers linked to consistent genomic regions are suitable for enrichment of favourable alleles for blackleg resistance in canola breeding programs.     

Quantitative Trait Loci for Thermal Time to Flowering and Photoperiod Responsiveness Discovered in Summer Annual-Type Brassica napus L

Time of flowering is a key adaptive trait in plants and is conditioned by the interaction of genes and environmental cues including length of photoperiod, ambient temperature and vernalisation. Here we investigated the photoperiod responsiveness of summer annual-types of Brassica napus (rapeseed, canola). A population of 131 doubled haploid lines derived from a cross between European and Australian parents was evaluated for days to flowering, thermal time to flowering (measured in degree-days) and the number of leaf nodes at flowering in a compact and efficient glasshouse-based experiment with replicated short and long day treatments. All three traits were under strong genetic control with heritability estimates ranging from 0.85–0.93. There was a very strong photoperiod effect with flowering in the population accelerated by 765 degree-days in the long day versus short day treatments. However, there was a strong genetic correlation of line effects (0.91) between the long and short day treatments and relatively low genotype x treatment interaction indicating that photoperiod had a similar effect across the population. Bivariate analysis of thermal time to flowering in short and long days revealed three main effect quantitative trait loci (QTLs) that accounted for 57.7% of the variation in the population and no significant interaction QTLs. These results provided insight into the contrasting adaptations of Australian and European varieties. Both parents responded to photoperiod and their alleles shifted the population to earlier flowering under long days. In addition, segregation of QTLs in the population caused wide transgressive segregation in thermal time to flowering. Potential candidate flowering time homologues located near QTLs were identified with the aid of the Brassica rapa reference genome sequence. We discuss how these results will help to guide the breeding of summer annual types of B. napus adapted to new and changing environments.     

Comparative mapping of loci controlling winter survival and related traits in oilseed Brassica rapa and B. napus

Winter survival is an important characteristic of oilseedBrassica that is seeded in the fall in northern climates,and it may be affected by genetic variation for other cold-regulated traits,such as freezing tolerance and vernalization responsive flowering time. Weanalyzed immortalized populations of oilseed Brassica rapa(recombinant inbred lines) and B. napus (double haploidlines) derived from crosses of annual and biennial types in order to comparethe map positions and effects of quantitative trait loci controlling wintersurvival, nonacclimated and acclimated freezing tolerances, and flowering time.The B. napus population was evaluated in multiple winters,and six of the 16 total significant QTL for winter survival were detected inmore than one winter. Correspondence in the map positions of QTL controllingdifferent traits within species provided evidence that some alleles causinggreater acclimated freezing tolerance and later flowering time also contributedto increased winter survival. Correspondence in the map positions of QTLbetween species provided evidence for allelic variation at homologous loci inB. rapa and B. napus. The potentialrole of some candidate genes in regulating these traits is discussed.     

A Brassica rapa Linkage Map of EST-based SNP Markers for Identification of Candidate Genes Controlling Flowering Time and Leaf Morphological Traits

For identification of genes responsible for varietal differences in flowering time and leaf morphological traits, we constructed a linkage map of Brassica rapa DNA markers including 170 EST-based markers, 12 SSR markers, and 59 BAC sequence-based markers, of which 151 are single nucleotide polymorphism (SNP) markers. By BLASTN, 223 markers were shown to have homologous regions in Arabidopsis thaliana, and these homologous loci covered nearly the whole genome of A. thaliana. Synteny analysis between B. rapa and A. thaliana revealed 33 large syntenic regions. Three quantitative trait loci (QTLs) for flowering time were detected. BrFLC1 and BrFLC2 were linked to the QTLs for bolting time, budding time, and flowering time. Three SNPs in the promoter, which may be the cause of low expression of BrFLC2 in the early-flowering parental line, were identified. For leaf lobe depth and leaf hairiness, one major QTL corresponding to a syntenic region containing GIBBERELLIN 20 OXIDASE 3 and one major QTL containing BrGL1, respectively, were detected. Analysis of nucleotide sequences and expression of these genes suggested possible involvement of these genes in leaf morphological traits.     

Genetic architecture of the circadian clock and flowering time in <Emphasis Type="Italic">Brassica rapa</Emphasis>

The circadian clock serves to coordinate physiology and behavior with the diurnal cycles derived from the daily rotation of the earth. In plants, circadian rhythms contribute to growth and yield and, hence, to both agricultural productivity and evolutionary fitness. Arabidopsis thaliana has served as a tractable model species in which to dissect clock mechanism and function, but it now becomes important to define the extent to which the Arabidopsis model can be extrapolated to other species, including crops. Accordingly, we have extended our studies to the close Arabidopsis relative and crop species, Brassica rapa. We have investigated natural variation in circadian function and flowering time among multiple B. rapa collections. There is wide variation in clock function, based on a robust rhythm in cotyledon movement, within a collection of B. rapa accessions, wild populations and recombinant inbred lines (RILs) derived from a cross between parents from two distinct subspecies, a rapid cycling Chinese cabbage (ssp. pekinensis) and a Yellow Sarson oilseed (ssp. trilocularis). We further analyzed the RILs to identify the quantitative trait loci (QTL) responsible for this natural variation in clock period and temperature compensation, as well as for flowering time under different temperature and day length settings. Most clock and flowering-time QTL mapped to overlapping chromosomal loci. We have exploited micro-synteny between the Arabidopsis and B. rapa genomes to identify candidate genes for these QTL.     

The quantitative-genetic and QTL architecture of trait integration and modularity in Brassica rapa across simulated seasonal settings

Within organisms, groups of traits with different functions are frequently modular, such that variation among modules is independent and variation within modules is tightly integrated, or correlated. Here, we investigated patterns of trait integration and modularity in Brassica rapa in response to three simulated seasonal temperature/photoperiod conditions. The goals of this research were to use trait correlations to understand patterns of trait integration and modularity within and among floral, vegetative and phenological traits of B. rapa in each of three treatments, to examine the QTL architecture underlying patterns of trait integration and modularity, and to quantify how variation in temperature and photoperiod affects the correlation structure and QTL architecture of traits. All floral organs of B. rapa were strongly correlated, and contrary to expectations, floral and vegetative traits were also correlated. Extensive QTL co-localization suggests that covariation of these traits is likely due to pleiotropy, although physically linked loci that independently affect individual traits cannot be ruled out. Across treatments, the structure of genotypic and QTL correlations was generally conserved. Any observed variation in genetic architecture arose from genotype × environment interactions (GEIs) and attendant QTL × E in response to temperature but not photoperiod.     

Mapping of days to flower and seed yield in spring oilseed <Emphasis Type="Italic">Brassica napus</Emphasis> carrying genome content introgressed from <Emphasis Type="Italic">Brassica oleracea</Emphasis>

Earliness of flowering and maturity and high seed yield are important objectives of breeding spring Brassica napus canola. Previously, we have introgressed earliness of flowering from Brassica oleracea into spring B. napus canola through interspecific crossing between these two species. In this paper, we report quantitative trait locus (QTL) mapping of days to flower and seed yield by use of publicly available markers and markers designed based on flowering time genes and a doubled haploid population, derived from crossing of the spring canola parent and an early flowering line developed from a B. napus × B. oleracea cross, tested in nine field trials for over 5 years. Five genomic regions associated with days to flower were identified on C1, C2, C3, and C6 of which the single QTL of C1 was detected in all trials; in all cases, the allele introgressed from B. oleracea reduced the number of days to flower. BLASTn search in the Brassica genomes located the physical position of the QTL markers and identified putative flowering time genes in these regions. In the case of seed yield, ten QTL from eight linkage groups were detected; however, none could be consistently detected in all trials. The QTL region of C1 associated with days to flower did not show significant association with seed yield in more than 80% of the field trials; however, in a single trial, the allele introgressed from B. oleracea exerted a negative effect on seed yield. Thus, the genomic regions and molecular markers identified in this research could potentially be used in breeding for the development of early flowering B. napus canola cultivars without affecting seed yield in a majority of the environments.     

Genome-Wide Delineation of Natural Variation for Pod Shatter Resistance in Brassica napus

Resistance to pod shattering (shatter resistance) is a target trait for global rapeseed (canola, Brassica napus L.), improvement programs to minimise grain loss in the mature standing crop, and during windrowing and mechanical harvest. We describe the genetic basis of natural variation for shatter resistance in B. napus and show that several quantitative trait loci (QTL) control this trait. To identify loci underlying shatter resistance, we used a novel genotyping-by-sequencing approach DArT-Seq. QTL analysis detected a total of 12 significant QTL on chromosomes A03, A07, A09, C03, C04, C06, and C08; which jointly account for approximately 57% of the genotypic variation in shatter resistance. Through Genome-Wide Association Studies, we show that a large number of loci, including those that are involved in shattering in Arabidopsis, account for variation in shatter resistance in diverse B. napus germplasm. Our results indicate that genetic diversity for shatter resistance genes in B. napus is limited; many of the genes that might control this trait were not included during the natural creation of this species, or were not retained during the domestication and selection process. We speculate that valuable diversity for this trait was lost during the natural creation of B. napus. To improve shatter resistance, breeders will need to target the introduction of useful alleles especially from genotypes of other related species of Brassica, such as those that we have identified.     

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

The introgression of winter germplasm into spring canola (Brassica napus L.) represents a novel approach to improve seed yield of hybrid spring canola. In this study, quantitative trait loci (QTL) for seed yield and other traits were genetically mapped to determine the effects of genomic regions introgressed from winter germplasm into spring canola. Plant materials used comprised of two populations of doubled haploid (DH) lines having winter germplasm introgression from two related French winter cultivars and their testcrosses with a spring line used in commercial hybrids. These populations were evaluated for 2 years at two locations (Wisconsin, USA and Saskatchewan, Canada). Genetic linkage maps based on RFLP loci were constructed for each DH population. Six QTL were detected in the testcross populations for which the winter alleles increased seed yield. One of these QTL explained 11 and 19% of the phenotypic variation in the two Canadian environments. The winter allele for another QTL that increased seed yield was linked in coupling to a QTL allele for high glucosinolate content, suggesting that the transition of rapeseed into canola could have resulted in the loss of favorable seed yield alleles. Most QTL for which the introgressed allele decreased seed yield of hybrids mapped to genomic regions having homoeologous non-reciprocal transpositions. This suggests that allelic configurations created by these rearrangements might make an important contribution to genetic variation for complex traits in oilseed B. napus and could account for a portion of the heterotic effects in hybrids. Electronic Supplementary Material Supplementary material is available to authorised users in the online version of this article at .     

BrFLC5: a weak regulator of flowering time in <Emphasis Type="Italic">Brassica rapa</Emphasis>

Key message

A splicing site mutation in BrFLC5, a non-syntenic paralogue of FLOWERING LOCUS C, was demonstrated to be related to flowering time variation in Brassica rapa.

Abstract

Flowering time regulation in Brassica rapa is more complex than in Arabidopsis, as there are multiple paralogues of flowering time genes in B. rapa. Brassica rapa contains four FLOWERING LOCUS C (FLC) genes, three of which are syntenic orthologues of AtFLC, while BrFLC5 is not. BrFLC1, BrFLC2, and BrFLC3 have been reported to be involved in flowering time regulation. However, BrFLC5 has thus far been deemed a pseudogene. We detected two alternative splicing patterns of BrFLC5 resulting from a nucleotide mutation (G/A) at the first nucleotide of intron 3 (named as Pi3+1(G/A)). Genotyping of BrFLC5Pi3?+?1(G/A) for 301 B. rapa accessions showed that this single nucleotide polymorphism was significantly related to flowering time variation (p?<?0.001). In the collection, the frequency of the functional G allele (35.2%) was much lower than that of the nonfunctional A allele (59.1%); however, the frequency of the G allele was very high among the turnips (83.6%). An F2 population segregating at this locus was developed to analyze the genetic effect of BrFLC5. The result showed that the G allele individuals began to bolt two days later than the A allele individuals, indicating that BrFLC5 is a weak regulator of flowering time. BrFLC5 was expressed at the lowest level among the three analyzed BrFLCs. The late allele (G allele) was dominant to the early allele (A allele) at the BrFLC5 locus, which was in contrast to that of BrFLC1 and BrFLC2. This characteristic suggests that BrFLC5 would be more efficient for breeding premature bolting resistance in B. rapa.

     

Mapping loci controlling vernalization requirement in Brassica rapa

Brassica cultivars are classified as biennial or annual based on their requirement for a period of cold treatment (vernalization) to induce flowering. Genes controlling the vernalization requirement were identified in a Brassica rapa F₂ population derived from a cross between an annual and a biennial oilseed cultivar by using an RFLP linkage map and quantitative trait locus (QTL) analysis of flowering time in F₃ lines. Two genomic regions were strongly associated with variation for flowering time of unvernalized plants and alleles from the biennial parent in these regions delayed flowering. These QTLs had no significant effect on flowering time after plants were vernalized for 6 weeks, suggesting that they control flowering time through the requirement for vernalization. The two B. rapa linkage groups containing these QTLs had RFLP loci in common with two B. napus linkage groups that were shown previously to contain QTLs for flowering time. An RFLP locus detected by the cold-induced gene COR6.6 cloned from Arabidopsis thaliana mapped very near to one of the B. rapa QTLs for flowering time.     

Environmental and genetic interactions reveal FLOWERING LOCUS C as a modulator of the natural variation for the plasticity of flowering in Arabidopsis

The timing of flowering initiation depends strongly on the environment, a property termed as the plasticity of flowering. Such plasticity determines the adaptive potential of plants because it provides phenotypic buffer against environmental changes, and its natural variation contributes to evolutionary adaptation. We addressed the genetic mechanisms of the natural variation for this plasticity in Arabidopsis thaliana by analysing a population of recombinant inbred lines derived from Don‐0 and Ler accessions collected from distinct climates. Quantitative trait locus (QTL) mapping in four environmental conditions differing in photoperiod, vernalization treatment and ambient temperature detected the folllowing: (i) FLOWERING LOCUS C (FLC) as a large effect QTL affecting flowering time differentially in all environments; (ii) numerous QTL displaying smaller effects specifically in some conditions; and (iii) significant genetic interactions between FLC and other loci. Hence, the variation for the plasticity of flowering is determined by a combination of environmentally sensitive and specific QTL, and epistasis. Analysis of FLC from Don identified a new and more active allele likely caused by a cis‐regulatory deletion covering the non‐coding RNA COLDAIR. Further characterization of four FLC natural alleles showed different environmental and genetic interactions. Thus, FLC appears as a major modulator of the natural variation for the plasticity of flowering to multiple environmental factors.     

Identification of quantitative trait loci controlling developmental characteristics of Brassica oleracea L

A segregating population of F₁-derived doubled haploid (DH) lines of Brassica oleracea was used to detect and locate QTLs controlling 27 morphological and developmental traits, including leaf, flowering, axillary bud and stem characters. The population resulted from a cross between two very different B. oleracea crop types, an annual cauliflower and a biennial Brussels sprout. A principal component analysis (PCA), based on line means, allowed all the traits to be grouped into distinct categories according to the first five Principal Components. These were: leaf traits (PC1), flowering traits (PC2), axillary bud traits (PC3 and 5) and stem traits (PC4). Between zero and four putative QTL were located per trait, which individually explained between 6% and 43% of the additive genetic variation, using the multiple-marker regression approach to QTL mapping. For lamina width, bare petiole length and stem length two QTL with opposite effects were detected on the same linkage groups. Intra- and inter-specific comparative mapping using RFLP markers identified a QTL on linkage group O8 accounting for variation in vernalisation, which is probably synonymous with a QTL detected on linkage group N19 of Brassica napus. In addition, a QTL for petiole length detected on O3 of this study appeared to be homologous to a QTL detected on another B. oleracea genetic map (Camargo et al. 1995). Received: 28 March 2001 / Accepted: 25 June 2001     

Mapping quantitative trait loci for tissue culture response in VCS3M-DH population of Brassica rapa

Key message

Quantitative trait loci (QTL) controlling callus induction and plant regeneration were identified in the VCS3M-DH population of Brassica rapa.

Abstract

Quantitative trait loci (QTL) controlling callus induction and plant regeneration were identified in the VCS3M-DH population of Brassica rapa. The VCS3M-DH population showed wide and continuous variation in callus induction and shoot regeneration. Significant coefficient correlations were detected between these two parameters. Broad-sense heritability (h ²) for the two traits was around 0.7, indicating genetic regulation of regeneration ability in this population. In the composite interval mapping analysis, two QTLs for callus induction ability, qCi2 and qCi7, were mapped on chromosome A02 and A07, explaining 28.6 % of phenotypic variation. For plant regeneration, four QTLs, qPr6-1 qPr6-2, qPr7, and qPr9 were identified on chromosome A06, A07, and A09, which in total explained 50.1 % of phenotypic variation. Furthermore, 15 putative candidate genes were found on the interval of the six QTLs, which were related to various plant hormones, MADS-box genes, and serine/threonine related genes. These results provide important information to identify genes related to tissue culture ability in B. rapa.