Brassica species, including crops such as cabbage, turnip and oilseed, display enormous phenotypic variation. Brassica genomes have all undergone a whole-genome triplication (WGT) event with unknown effects on phenotype diversification. We resequenced 199 Brassica rapa and 119 Brassica oleracea accessions representing various morphotypes and identified signals of selection at the mesohexaploid subgenome level. For cabbage morphotypes with their typical leaf-heading trait, we identified four subgenome loci that show signs of parallel selection among subgenomes within B. rapa, as well as four such loci within B. oleracea. Fifteen subgenome loci are under selection and are shared by these two species. We also detected strong subgenome parallel selection linked to the domestication of the tuberous morphotypes, turnip (B. rapa) and kohlrabi (B. oleracea). Overall, we demonstrated that the mesohexaploidization of the two Brassica genomes contributed to their diversification into heading and tuber-forming morphotypes through convergent subgenome parallel selection of paralogous genes.

BackgroundFlowering time is an important trait in Brassica rapa crops. FLOWERING LOCUS C (FLC) is a MADS-box transcription factor that acts as a potent repressor of flowering. Expression of FLC is silenced when plants are exposed to low temperature, which activates flowering. There are four copies of FLC in B. rapa. Analyses of different segregating populations have suggested that BraA.FLC.a (BrFLC1) and BraA.FLC.b (BrFLC2) play major roles in controlling flowering time in B. rapa.ResultsWe analyzed the BrFLC2 sequence in nine B. rapa accessions, and identified a 57-bp insertion/deletion (InDel) across exon 4 and intron 4 resulting in a non-functional allele. In total, three types of transcripts were identified for this mutated BrFLC2 allele. The InDel was used to develop a PCR-based marker, which was used to screen a collection of 159 B. rapa accessions. The deletion genotype was present only in oil-type B. rapa, including ssp. oleifera and ssp. tricolaris, and not in other subspecies. The deletion genotype was significantly correlated with variation in flowering time. In contrast, the reported splicing site variation in BrFLC1, which also leads to a non-functional locus, was detected but not correlated with variation in flowering time in oil-type B. rapa, although it was correlated with variation in flowering time in vegetable-type B. rapa.ConclusionsOur results suggest that the naturally occurring deletion mutation across exon 4 and intron 4 in BrFLC2 gene contributes greatly to variation in flowering time in oil-type B. rapa. The observed different relationship between BrFLC1 or BrFLC2 and flowering time variation indicates that the control of flowering time has evolved separately between oil-type and vegetable-type B. rapa groups.

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. 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.