Comparison of Flowering Time Genes in Brassica Rapa, B. Napus and Arabidopsis Thaliana

The major difference between annual and biennial cultivars of oilseed Brassica napus and B. rapa is conferred by genes controlling vernalization-responsive flowering time. These genes were compared between the species by aligning the map positions of flowering time quantitative trait loci (QTLs) detected in a segregating population of each species. The results suggest that two major QTLs identified in B. rapa correspond to two major QTLs identified in B. napus. Since B. rapa is one of the hypothesized diploid parents of the amphidiploid B. napus, the vernalization requirement of B. napus probably originated from B. rapa. Brassica genes also were compared to flowering time genes in Arabidopsis thaliana by mapping RFLP loci with the same probes in both B. napus and Arabidopsis. The region containing one pair of Brassica QTLs was collinear with the top of chromosome 5 in A. thaliana where flowering time genes FLC, FY and CO are located. The region containing the second pair of QTLs showed fractured collinearity with several regions of the Arabidopsis genome, including the top of chromosome 4 where FRI is located. Thus, these Brassica genes may correspond to two genes (FLC and FRI) that regulate flowering time in the latest flowering ecotypes of Arabidopsis.     

Conserved structure and function of the Arabidopsis flowering time gene CONSTANS in Brassica napus

The Arabidopsis thaliana CONSTANS (CO) gene which promotes flowering in long days was recently isolated by chromosome walking. The mapping of QTLs controlling flowering time in Brassica species has identified genomic regions that contain homologues of the CO gene. Four genes homologous to the Arabidopsis CO gene were isolated from a pair of homoeologous loci in each of two doubled-haploid Brassica napus lines displaying different flowering times, N-o-1 and N-o-9. The four genes, BnCOa1, BnCOa9, BnCOb1 and BnCOb9, are located on linkage groups N10 and N19, and are highly similar to each other and to the Arabidopsis CO gene. Two regions of the proteins are particularly well conserved, a N-terminal region with two putative zinc fingers and a C-terminal region which may contain a nuclear localization signal. All four genes appear to be expressed in B. napus. The BnCOa1 allele was shown to complement the co-2 mutation in Arabidopsis in a dosage-dependent manner causing earlier flowering than in wild type under both long- and short-day conditions.     

Arabidopsis thaliana: a source of candidate disease-resistance genes for Brassica napus.

Common structural and amino acid motifs among cloned plant disease-resistance genes (R genes), have made it possible to identify putative disease-resistance sequences based on DNA sequence identity. Mapping of such R-gene homologues will identify candidate disease-resistance loci to expedite map-based cloning strategies in complex crop genomes. Arabidopsis thaliana expressed sequence tags (ESTs) with homology to cloned plant R genes (R-ESTs), were mapped in both A. thaliana and Brassica napus to identify candidate R-gene loci and investigate intergenomic collinearity. Brassica R-gene homologous sequences were also mapped in B. napus. In total, 103 R-EST loci and 36 Brassica R-gene homologous loci were positioned on the N-fo-61-9 B. napus genetic map, and 48 R-EST loci positioned on the Columbia x Landsberg A. thaliana map. The mapped loci identified collinear regions between Arabidopsis and Brassica which had been observed in previous comparative mapping studies; the detection of syntenic genomic regions indicated that there was no apparent rapid divergence of the identified genomic regions housing the R-EST loci.     

Comparative mapping between Arabidopsis thaliana and Brassica nigra indicates that Brassica genomes have evolved through extensive genome replication accompanied by chromosome fusions and frequent rearrangements.

Chromosome organization and evolution in the Brassicaceae family was studied using comparative linkage mapping. A total of 160 mapped Arabidopsis thaliana DNA fragments identified 284 homologous loci covering 751 cM in Brassica nigra. The data support that modern diploid Brassica species are descended from a hexaploid ancestor, and that the A. thaliana genome is similar in structure and complexity to those of each of the hypothetical diploid progenitors of the proposed hexaploid. Thus, the Brassica lineage probably went through a triplication after the divergence of the lineages leading to A. thaliana and B. nigra. These duplications were also accompanied by an exceptionally high rate of chromosomal rearrangements. The average length of conserved segments between A. thaliana and B. nigra was estimated at 8 cM. This estimate corresponds to approximately 90 rearrangements since the divergence of the two species. The estimated rate of chromosomal rearrangements is higher than any previously reported data based on comparative mapping. Despite the large number of rearrangements, fine-scale comparative mapping between model plant A. thaliana and Brassica crops is likely to result in the identification of a large number of genes that affect important traits in Brassica crops.     

Multiple flowering time QTLs within several Brassica species could be the result of duplicated copies of one ancestral gene.

Quantitative trait locus (QTL) analysis was used to study the evolution of genes controlling the timing of flowering in four Brassica genomes that are all extensively replicated. Comparative mapping showed that a chromosomal region from the top of Arabidopsis thaliana chromosome 5 corresponded to three homoeologous copies in each of the diploid species Brassica nigra, B. oleracea, and B. rapa and six copies in the amphidiploid B. juncea. QTLs were detected in two of the three replicated segments in each diploid genome and in three of the six replicated segments in B. juncea. These results indicate that, for the studied trait, multiple QTLs resulting from genome duplication is the rule rather than the exception. Brassica homologues to two candidate genes (CO and FLC) identified from the corresponding A. thaliana region were mapped. CO homologues mapped close to the QTL peaks in eight of nine QTLs, while FLC homologues mapped farther away in those cases where the mapping resolution allowed a comparison. Thus, our data are consistent with the hypothesis that all the major QTLs we detected in the different species of Brassica could be the result of duplicated copies of the same ancestral gene, possibly the ancestor of CO.     

Naturally occurring indel variation in the Brassica nigra COL1 gene is associated with variation in flowering time

Previous QTL mapping identified a Brassica nigra homolog to Arabidopsis thaliana CO as a candidate gene affecting flowering time in B. nigra. Transformation of an A. thaliana co mutant with two different alleles of the B. nigra CO (Bni COa) homolog, one from an early-flowering B. nigra plant and one from a late one, did not show any differential effect of the two alleles on flowering time. The DNA sequence of the coding region of the two alleles was also identical, showing that nucleotide variation influencing flowering time must reside outside the coding region of Bni COa. In contrast, the nucleotide sequence of the B. nigra COL1 (Bni COL1) gene located 3.5 kb upstream of Bni COa was highly diverged between the alleles from early and late plants. One indel polymorphism in the Bni COL1 coding region, present in several natural populations of B. nigra, displayed a significant association with flowering time within a majority of these populations. These data indicate that a quantitative trait nucleotide (QTN) affecting flowering time is located within or close to the Bni COL1 gene. The intergenic sequence between Bni COL1 and Bni COa displayed a prominent peak of divergence 1 kb downstream of the Bni COL1 coding region. This region could contain regulatory elements for the downstream Bni COa gene. Our data suggest that a naturally occurring QTN for flowering time affects the function or expression of either Bni COL1 or Bni COa.     

Flowering time variation in oilseed rape (Brassica napus L.) is associated with allelic variation in the FRIGIDA homologue BnaA.FRI.a

Oilseed rape (Brassica napus L.) is a major oil crop which is grown worldwide. Adaptation to different environments and regional climatic conditions involves variation in the regulation of flowering time. Winter types have a strong vernalization requirement whereas semi-winter and spring types have a low vernalization requirement or flower without exposure to cold, respectively. In Arabidopsis thaliana, FRIGIDA (FRI) is a key regulator which inhibits floral transition through activation of FLOWERING LOCUS C (FLC), a central repressor of flowering which controls vernalization requirement and response. Here, four FRI homologues in B. napus were identified by BAC library screening and PCR-based cloning. While all homologues are expressed, two genes were found to be differentially expressed in aerial plant organs. One of these, BnaA.FRI.a, was mapped to a region on chromosome A03 which co-localizes with a major flowering time quantitative trait locus in multiple environments in a doubled-haploid mapping population. Association analysis of BnaA.FRI.a revealed that six SNPs, including at least one at a putative functional site, and one haplotype block, respectively, are associated with flowering time variation in 248 accessions, with flowering times differing by 13-19 d between extreme haplotypes. The results from both linkage analysis and association mapping indicate that BnaA.FRI.a is a major determinant of flowering time in oilseed rape, and suggest further that this gene also contributes to the differentiation between growth types. The putative functional polymorphisms identified here may facilitate adaptation of this crop to specific environments through marker-assisted breeding.     

Contrasting genome organisation: two regions of the Brassica oleracea genome compared with collinear regions of the Arabidopsis thaliana genome.

Brassica crop species are of worldwide importance and are closely related to the model plant Arabidopsis thaliana for which the complete genome sequence has recently been established. We investigated collinearity of marker order by comparing two contrasting regions of the Brassica oleracea genome with homologous regions of A. thaliana. Although there is widespread replication of marker loci in both A. thaliana and B. oleracea, we found that a combination of genetic markers mapped in B. oleracea, including RFLPs, CAPS, and SSRs allowed comparison and interpretation of medium-scale chromosomal organisation and rearrangements. The interpretation of data was facilitated by hybridising probes onto the whole A. thaliana genome, as represented by BAC contigs. Twenty marker loci were sampled from the whole length of the shortest B. oleracea linkage group, 06, and 21 from a 30.4-cM section of the longest linkage group, 03. There is evidence of locus duplication on linkage group 06. Locus order is well conserved between a putative duplicated region of 10.5 cM and a discrete region comprising 25 cM of A. thaliana chromosome I. This was supported by evidence from seven paralogous loci, three of which were duplicated in a 30.6-cM region of linkage group 06. The pattern of locus order for the remainder of linkage group 06 and the sampled section of linkage group 03 was more complex when compared with the A. thaliana genome. Although there was some conservation of locus order between markers on linkage group 03 and approximately 9 cM of A. thaliana chromosome I, this was superimposed upon a complex pattern of additional loci that were replicated in both A. thaliana and B. oleracea. The results are discussed in the context of the ability to use collinear information to assist map-based cloning.     

Comparative physical mapping of segments of the genome of Brassica oleracea var. alboglabra that are homoeologous to sequenced regions of chromosomes 4 and 5 of Arabidopsis thaliana

Due to their relatedness to Arabidopsis thaliana (Arabidopsis), the cultivated Brassica species represent the first group of crops with which to evaluate comparative genomics approaches to understanding biological processes and manipulating traits. We have constructed a high-quality binary BAC library (JBo) from genomic DNA of Brassica oleracea var. alboglabra, in order to underpin such investigations. Using the Arabidopsis genome sequence and clones from the JBo library, we have analysed aspects of gene conservation and microsynteny between a 222 kb region of the genome of Arabidopsis and homoeologous segments of the genome of B. oleracea. All 19 predicted genes tested were found to hybridize to clones in the JBo library, indicating a high level of gene conservation. Further analyses and physical mapping with the BAC clones identified allowed us to construct clone contig maps and analyse in detail the gene content and organization in the set of paralogous segments identified in the genome of B. oleracea. Extensive divergence of gene content was observed, both between the B. oleracea paralogous segments and between them and their homoeologous segment within the genome of Arabidopsis. However, the genes present show highly conserved collinearity with their orthologues in the genome of Arabidopsis. We have identified one example of a Brassica gene in a non-collinear position and one rearrangement. Some of the genes not present in the discernible homoeologous regions appear to be located elsewhere in the B. oleracea genome. The implications of our findings for comparative map-based cloning of genes from crop species are discussed.     

Physical organization of the major duplication on Brassica oleracea chromosome O6 revealed through fluorescence in situ hybridization with Arabidopsis and Brassica BAC probes.

The close relationship between Brassica oleracea and Arabidopsis thaliana has been used to explore the genetic and physical collinearity of the two species, focusing on an inverted segmental chromosome duplication within linkage group O6 of B. oleracea. Genetic evidence suggests that these segments share a common origin with a region of Arabidopsis chromosome 1. Brassica oleracea and Arabidopsis bacterial artificial chromosome probes have been used for fluorescence in situ hybridization analysis of B. oleracea pachytene chromosomes to further characterize the inverted duplication. This has been highly effective in increasing the local resolution of the cytogenetic map. We have shown that the physical order of corresponding genetic markers is highly conserved between the duplicated regions in B. oleracea and the physical lengths of the regions at pachytene are similar, while the genetic distances are considerably different. The physical marker order is also well conserved between Arabidopsis and B. oleracea, with only one short inversion identified. Furthermore, the relative physical distances between the markers in one segment of B. oleracea and Arabidopsis have stayed approximately the same. The efficacy of using fluorescence in situ hybridization, together with other forms of physical and genetic mapping, for elucidating such issues relating to synteny is discussed.     

In silico integration of quantitative trait loci for seed yield and yield-related traits in Brassica napus

Mapping quantitative trait loci (QTLs) is a foundation for molecular marker-assisted selection and map-based gene cloning. During the past decade, numerous QTLs for seed yield (SY) and yield-related traits in Brassica napus L. have been identified. However, integration of these results in order to compare QTLs from different mapping populations has not been undertaken, due to the lack of common molecular markers between studies. Using previously reported Brassica rapa and Brassica oleracea genome sequences, we carried out in silico integration of 1,960 QTLs associated with 13 SY and yield-related traits from 15 B. napus mapping experiments over the last decade. A total of 736 SY and yield-related QTLs were mapped onto 283 loci in the A and C genomes of B. napus. These QTLs were unevenly distributed across the 19 B. napus chromosomes, with the most on chromosome A3 and the least on chromosome C6. Our integrated QTL map identified 142 loci where the conserved QTLs were detected and 25 multifunctional loci, mostly for the traits of flowering time (FT), plant height, 1,000-seed weight, maturity time and SY. These conserved QTLs and multifunctional loci may result from pleiotropism or clustered genes. At the same time, a total of 146 genes underlying the QTLs for FT and other yield-related traits were identified by comparative mapping with the Arabidopsis genome. These results facilitate the retrieval of B. napus SY and yield-related QTLs for research communities, increase the density of targeted QTL-linked markers, validate the existence of QTLs across different populations, and advance the fine mapping of genes.     

FISH-mapping of rDNAs and Arabidopsis BACs on pachytene complements of selected Brassicas.

To improve resolution of physical mapping on Brassica chromosomes, we have chosen the pachytene stage of meiosis where incompletely condensed bivalents are much longer than their counterparts at mitotic metaphase. Mapping with 5S and 45S rDNA sequences demonstrated the advantage of pachytene chromosomes in efficient physical mapping and confirmed the presence of a novel 5S rDNA locus in Brassica oleracea, initially identified by genetic mapping using restriction fragment length polymorphism (RFLP). Fluorescence in situ hybridization (FISH) analysis visualized the presence of the third 5S rDNA locus on the long arm of chromosome C2 and confirmed the earlier reports of two 45S rDNA loci in the B. oleracea genome. FISH mapping of low-copy sequences from the Arabidopsis thaliana bacterial artificial chromosome (BAC) clones on the B. oleracea chromosomes confirmed the expectation of efficient and precise physical mapping of meiotic bivalents based on data available from A. thaliana and indicated conserved organization of these two BAC sequences on two B. oleracea chromosomes. Based on the heterologous in situ hybridization with BACs and their mapping applied to long pachytene bivalents, a new approach in comparative analysis of Brassica and A. thaliana genomes is discussed.     

A sequence-tagged linkage map of Brassica rapa

A detailed genetic linkage map of Brassica rapa has been constructed containing 545 sequence-tagged loci covering 1287 cM, with an average mapping interval of 2.4 cM. The loci were identified using a combination of 520 RFLP and 25 PCR-based markers. RFLP probes were derived from 359 B. rapa EST clones and amplification products of 11 B. rapa and 26 Arabidopsis. Including 21 SSR markers provided anchors to previously published linkage maps for B. rapa and B. napus and is followed as the referenced mapping of R1-R10. The sequence-tagged markers allowed interpretation of the pattern of chromosome duplications within the B. rapa genome and comparison with Arabidopsis. A total of 62 EST markers showing a single RFLP band were mapped through 10 linkage groups, indicating that these can be valuable anchoring markers for chromosome-based genome sequencing of B. rapa. Other RFLP probes gave rise to 2-5 loci, inferring that B. rapa genome duplication is a general phenomenon through 10 chromosomes. The map includes five loci of FLC paralogues, which represent the previously reported BrFLC-1, -2, -3, and -5 and additionally identified BrFLC3 paralogues derived from local segmental duplication on R3.     

The association of flowering time quantitative trait loci with duplicated regions and candidate loci in Brassica oleracea.

A population of 150 doubled haploid lines of rapid cycling Brassica oleracea, derived from an F1 from a var. alboglabra x var. italica cross, was scored for flowering time in two trials. Using information on 82 mapped molecular markers, spread evenly across the nine linkage groups, QTL were identified at six locations; one each on linkage groups O2 and O3 and two each on linkage groups O5 and O9. In total, these QTL explained 58 and 93% of the genetical variation in the two trials. Three of these QTL, on linkage groups O2, O3, and O9, were situated in regions showing considerable homology both with each other and with chromosome regions of B. nigra that have been shown to affect flowering time. These same regions are all homologous to a single tract of Arabidopsis chromosome 5, which contains a number of the flowering-related genes, one or more of which may be candidates for the QTL found in Brassica.     

Recent degeneration of an old duplicated flowering time gene in Brassica nigra

Gene and genome duplications play a major role in the evolution of plant species. The Brassica nigra genome is highly replicated as a result of ancient polyploidization events. Two copies of the flowering time gene CONSTANS (COa and COb) have been identified in B. nigra, and previous studies showed that COa is functional. In the present study, the polymorphism of 92 COb alleles sampled in seven populations was analyzed. Both polymorphism and recombination levels were elevated and varied strongly among populations and 8% of COb alleles exhibit apparently disabling mutations. Sequence data, however, do not provide unambiguous support for the presence of relaxed selective constraint on COb as compared to known functional CO genes. On the one hand, some of the disabling mutations reached high-frequency arguing for a loss of function but, on the other hand, the ratio of nonsynonymous to synonymous nucleotide polymorphism and diversity is low and similar to that observed in other B. nigra CO and CO-like genes, supporting the conservation of some function. We also showed that COb is still transcribed. Finally, the flowering time of Arabidopsis thaliana co mutant plants transformed with COb alleles with and without apparent disabling mutations was similar. We propose that COb was retained for a long period after duplication, but a recent fixation of a detrimental mutation, possibly as an effect of a bottleneck, resulted in its nonfunctionalization. We also speculate as to the presence of subsequent selection for rapid degeneration of the gene.     

Chromosomal mapping of <Emphasis Type="Italic"> Brassica oleracea </Emphasis>based on ESTs from <Emphasis Type="Italic"> Arabidopsis thaliana</Emphasis>: complexity of the comparative map

Expressed sequence tags (ESTs) from the Arabidopsis thaliana sequencing project were used to construct a genetic RFLP map for Brassica oleracea. Of the 110 A. thaliana ESTs tested, 95 were found to be informative RFLP probes in map construction. In total, 212 new loci corresponding to the 95 ESTs were added to the existing genetic map of B. oleracea. The enriched map covers all nine basic linkage groups and confirms that the chromosomes of B. oleracea and A. thaliana are similar in linear organization. However, varying levels of sequence conservation between the chromosomes of B. oleracea and A. thaliana were detected in different regions of the genomes. Long conserved regions encompassing entire chromosome arms in both genomes were identified; these are probably shared by descent. On the other hand, extensive rearrangements were observed in numerous chromosome regions, producing a mosaic of A. thaliana-like segments in the genome of Brassica. The presence of extensive chromosome duplication in A. thaliana was taken into consideration in the construction of the comparative maps of B. oleracea and A. thaliana.     

CONSTANS activates SUPPRESSOR OF OVEREXPRESSION OF CONSTANS 1 through FLOWERING LOCUS T to promote flowering in Arabidopsis

CONSTANS (CO) regulates flowering time by positively regulating expression of two floral integrators, FLOWERING LOCUS T (FT) and SUPPRESSOR OF OVEREXPRESSION OF CO 1 (SOC1), in Arabidopsis (Arabidopsis thaliana). FT and SOC1 have been proposed to act in parallel pathways downstream of CO based on genetic analysis using weak ft alleles, since ft soc1 double mutants showed an additive effect in suppressing the early flowering of CO overexpressor plants. However, this genetic analysis was inconsistent with the sequential induction pattern of FT and SOC1 found in inducible CO overexpressor plants. Hence, to identify genetic interactions of CO, FT, and SOC1, we carried out genetic and expression analyses with a newly isolated T-DNA allele of FT, ft-10. We found that ft-10 almost completely suppressed the early flowering phenotype of CO overexpressor plants, whereas soc1-2 partially suppressed the phenotype, suggesting that FT is the major output of CO. Expression of SOC1 was altered in gain- or loss-of-function mutants of FT, whereas expression of FT remained unchanged in gain- or loss-of-function mutants of SOC1, suggesting that FT positively regulates SOC1 to promote flowering. In addition, inactivation of FT caused down-regulation of SOC1 even in plants overexpressing CO, indicating that FT is required for SOC1 induction by CO. Taken together, these data suggest that CO activates SOC1 through FT to promote flowering in Arabidopsis.     

Assessing the level of collinearity between Arabidopsis thaliana and Brassica napus for A. thaliana chromosome 5.

This study describes a comprehensive comparison of chromosome 5 of the model crucifer Arabidopsis with the genome of its amphidiploid crop relative Brassica napus and introduces the use of in silico sequence homology to identify conserved loci between the two species. A region of chromosome 5, spanning 8 Mb, was found in six highly conserved copies in the B. napus genome. A single inversion appeared to be the predominant rearrangement that had separated the two lineages leading to the formation of Arabidopsis chromosome 5 and its homologues in B. napus. The observed results could be explained by the fusion of three ancestral genomes with strong similarities to modern-day Arabidopsis to generate the constituent diploid genomes of B. napus. This supports the hypothesis that the diploid Brassica genomes evolved from a common hexaploid ancestor. Alignment of the genetic linkage map of B. napus with the genomic sequence of Arabidopsis indicated that for specific regions a genetic distance of 1 cM in B. napus was equivalent to 285 Kb of Arabidopsis DNA sequence. This analysis strongly supports the application of Arabidopsis as a tool in marker development, map-based gene cloning, and candidate gene identification for the larger genomes of Brassica crop species.