Flowering Time Gene Variation in Brassica Species Shows Evolutionary Principles

A transcriptomic time-series reveals differing trajectories during pre-floral development in the apex and leaf in winter and spring varieties of Brassica napus

Oilseed rape (Brassica napus) is an important global oil crop, with spring and winter varieties grown commercially. To understand the transcriptomic differences between these varieties, we collected transcriptomes from apex and leaf tissue from a spring variety, Westar, and a winter variety, Tapidor, before, during, and after vernalisation treatment, until the plants flowered. Large transcriptomic differences were noted in both varieties during the vernalisation treatment because of temperature and day length changes. Transcriptomic alignment revealed that the apex transcriptome reflects developmental state, whereas the leaf transcriptome is more closely aligned to the age of the plant. Similar numbers of copies of genes were expressed in both varieties during the time series, although key flowering time genes exhibited expression pattern differences. BnaFLC copies on A2 and A10 are the best candidates for the increased vernalisation requirement of Tapidor. Other BnaFLC copies show tissue-dependent reactivation of expression post-cold, with these dynamics suggesting some copies have retained or acquired a perennial nature. BnaSOC1 genes, also related to the vernalisation pathway, have expression profiles which suggest tissue subfunctionalisation. This understanding may help to breed varieties with more consistent or robust vernalisation responses, of special importance due to the milder winters resulting from climate change.

Auxin-regulated timing of transition from vegetative to reproductive growth in rapeseed (Brassica napus L.) under different nitrogen application rates

Accelerating the differentiation of floral meristem (FM) from shoot apical meristems (SAM) which determines the conversion from vegetative to reproductive growth is of great significance for the production of rapeseed (Brassica napus L.). In this research, the mechanisms of different nitrogen (N) application rates (low N, N1; normal N, N2; and high N, N3) on different FM development stages triggering the regulation of FM differentiation genes through the auxin biosynthetic and signal transduction were investigated. We found that the stage of FM differentiation, which was identified through a stereomicroscope and scanning electron microscope, came 4 and 7 days earlier under high N rate than under normal and low N levels, with the seed yield increased by 11.1 and 22.6%, respectively. Analysis of the auxin and its derivatives contents showed that the main biosynthesis way of auxin was the indole acetaldehyde oxime (IAOx) pathway, with 3-Indole acetonitrile dramatically accumulated during FM differentiation. At the same time, an obvious decrease of IAA contents at each FM differentiation stage was detected, and then gradually rose. Results of the expression of genes involved in auxin biosynthesis, auxin signaling transduction, and FM identification under five FM differentiation stages and three nitrogen application rates showed that genes involved in auxin biosynthesis were regulated before the FM differentiation stage, while the regulation of FM identity genes appeared mainly at the middle and later periods of the five stages, and the regulation level of genes varied under different N rates. Taken together, a high nitrogen rate could accelerate the initiation of FM differentiation, and auxin involved a lot in this regulation.

Refining bulk segregant analyses: ontology-mediated discovery of flowering time genes in Brassica oleracea

BACKGROUND

Bulk segregant analysis (BSA) can help identify quantitative trait loci (QTLs), but this may result in substantial bycatch of functionally irrelevant genes.

RESULTS

Here we develop a Gene Ontology-mediated approach to zoom in on specific genes located inside QTLs identified by BSA as implicated in a continuous trait. We apply this to a novel experimental system: flowering time in the giant woody Jersey kale, which we phenotyped in four bulks of flowering onset. Our inferred QTLs yielded tens of thousands of candidate genes. We reduced this by two orders of magnitude by focusing on genes annotated with terms contained within relevant subgraphs of the Gene Ontology. A pathway enrichment test then led to the circadian rhythm pathway. The genes that enriched this pathway are attested from previous research as regulating flowering time. Within that pathway, the genes CCA1, FT, and TSF were identified as having functionally significant variation compared to Arabidopsis. We validated and confirmed our ontology-mediated results through genome sequencing and homology-based SNP analysis. However, our ontology-mediated approach produced additional genes of putative importance, showing that the approach aids in exploration and discovery.

CONCLUSIONS

Our method is potentially applicable to the study of other complex traits and we therefore make our workflows available as open-source code and a reusable Docker container.

Conserved and distinct roles of H3K27me3 demethylases regulating flowering time in Brassica rapa

Epigenetic regulation is necessary for optimal organism development and preservation of gene expression profiles in the cell. In plants, the trimethylation of histone H3 lysine 27 (H3K27me3) is a silencing epigenetic mark relevant for developmental transitions like flowering. The floral transition is a key agronomic trait; however, the epigenetic mechanisms of flowering time regulation in crops remain poorly understood. Here we study the Jumonji H3K27me3 demethylases BraA.REF6 and BraA.ELF6 in Brassica rapa. Phenotypic characterization of novel mutant lines and genome-wide H3K27me3 chromatin immunoprecipitation and transcriptomic analyses indicated that BraA.REF6 plays a greater role than BraA.ELF6 in fine-tuning H3K27me3 levels. In addition, we found that braA.elf6 mutants were early flowering due to high H3K27me3 levels at B. rapa homologs of the floral repressor FLC. Unlike mutations in Arabidopsis thaliana, braA.ref6 mutants were late flowering without altering the expression of B. rapa FLC genes. Remarkably, we found that BraA.REF6 regulated a number of gibberellic acid (GA) biosynthetic genes, including a homolog of GA1, and that GA-treatment complemented the late flowering mutant phenotype. This study increases our understanding of the epigenetic regulation of flowering time in B. rapa, highlighting conserved and distinct regulatory mechanisms between model and crop species.

A premature stop codon in BrFLC2 transcript results in early flowering in oilseed-type Brassica rapa plants

Nonsense-mediated mRNA decay (NMD)-mediated degradation of BrFLC2 transcripts is the main cause of rapid flowering of oilseed-type B. rapa 'LP08' plants. Many Brassica species require vernalization (long-term winter-like cooling) for transition to the reproductive stage. In the past several decades, scientific efforts have been made to discern the molecular mechanisms underlying vernalization in many species. Thus, to identify the key regulators required for vernalization in Brassica rapa L., we constructed a linkage map composed of 7833 single nucleotide polymorphism markers using the late-flowering Chinese cabbage (B. rapa L. ssp. pekinensis) inbred line 'Chiifu' and the early-flowering yellow sarson (B. rapa L. ssp. trilocularis) line 'LP08' and identified a single major QTL on the upper-arm of the chromosome A02. In addition, we compared the transcriptomes of the lines 'Chiifu' and 'LP08' at five vernalization time points, including both non-vernalized and post-vernalization conditions. We observed that BrFLC2 was significantly downregulated in the early flowering 'LP08' and had two deletion sites (one at 4th exon and the other at 3' downstream region) around the BrFLC2 genomic region compared with the BrFLC2 genomic region in 'Chiifu'. Large deletion at 3' downstream region did not significantly affect transcription of both sense BrFLC2 transcript and antisense transcript, BrFLC2as along vernalization time course. However, the other deletion at 4th exon of BrFLC2 resulted in the generation of premature stop codon in BrFLC2 transcript in LP08 line. Cycloheximide treatment of LP08 line showed the de-repressed level of BrFLC2 in LP08, suggesting that low transcript level of BrFLC2 in LP08 might be caused by nonsense-mediated mRNA decay removing the nonsense transcript of BrFLC2. Collectively, this study provides a better understanding of the molecular mechanisms underlying floral transition in B. rapa.

Haplotype variations of major flowering time genes in quinoa unveil their role in the adaptation to different environmental conditions

Response to photoperiod is of major importance in crop production. It defines the adaptation of plants to local environments. Quinoa is a short-day plant which had been domesticated in the Andeans regions. We wanted to understand the adaptation to long-day conditions by studying orthologues of two major flowering time regulators of Arabidopsis, FLOWERING LOCUS T (FT) and CONSTANS (CO) in quinoa accessions with contrasting photoperiod response. By searching the quinoa reference genome sequence, we identified 24 FT and six CO homologs. CqFT genes displayed remarkably different expression patterns between long- and short-day conditions, whereas the influence of the photoperiod on CqCOL expressions was moderate. Cultivation of 276 quinoa accessions under short- and long-day conditions revealed great differences in photoperiod sensitivity. After sequencing their genomes, we identified large sequence variations in 12 flowering time genes. We found non-random distribution of haplotypes across accessions from different geographical origins, highlighting the role of CqFT and CqCOL genes in the adaptation to different day-length conditions. We identified five haplotypes causing early flowering under long days. This study provides assets for quinoa breeding because superior haplotypes can be assembled in a predictive breeding approach to produce well-adapted early flowering lines under long-day photoperiods.

Genome-Wide Analysis of PEBP Genes in Dendrobium huoshanense: Unveiling the Antagonistic Functions of FT/TFL1 in Flowering Time

Dendrobium is a semi-shade epiphytic Orchidaceae herb with important ornamental and medicinal value. Parts of the cultivation of Dendrobium germplasm resources, as well as the identification of medicinal components, are more studied, but the functional characterization of the flowering regulation in Dendrobium plants is less reported. Here, six PEBP family genes (DhFT3, DhFT1, DhMFT, DhTFL1b, DhFT2, and DhTFL1a) were identified from the Dendrobium huoshanense genome. The chromosome-level mapping showed that these genes were sequentially distributed on chromosomes 6, 9, 15, and 17. The paralogous gene DhTFL1b corresponded to DhTFL1a, which was determined through tandem duplication. The gene structure and conserved motif of DhPEBP indicated five PEBP genes apart from DhMFT contained four exons and three introns entirely. The phylogeny analysis showed that the PEBP gene family in A. thaliana, O. sativa, Z. mays, S. lycopersicum, and P. equestris were classified into three subclades, FT, TFL, and MFT, which maintained a high homology with D. huoshanense. The conserved domain of the amino acid demonstrated that two highly conserved short motifs (DPDXP and GXHR) embed in DhPEBPs, which may contribute to the conformation of the ligand binding bag. The 86th position of DhFTs was tyrosine (Y), while the 83th and 87th of DhTFL1s belonged to histidine (H), suggesting they should have distinct functions in flowering regulation. The promoter of six DhPEBPs contained several cis-elements related to hormone induction, light response, and abiotic stress, which indicated they could be regulated by the environmental stress and endogenous signaling pathways. The qRT-PCR analysis of DhPEBPs in short-term days induced by GA indicated the gene expressions of all DhFTs were gradually increased, whereas the expression of DhTFL1 was decreased. The results implied that DhPEBPs have various regulatory functions in modulating flowering, which will provide a scientific reference for the flowering regulation of Dendrobium plants.

Nitrogen Signaling Genes and SOC1 Determine the Flowering Time in a Reciprocal Negative Feedback Loop in Chinese Cabbage (Brassica rapa L.) Based on CRISPR/Cas9-Mediated Mutagenesis of Multiple BrSOC1 Homologs

Precise flowering timing is critical for the plant life cycle. Here, we examined the molecular mechanisms and regulatory network associated with flowering in Chinese cabbage (Brassica rapa L.) by comparative transcriptome profiling of two Chinese cabbage inbred lines, “4004” (early bolting) and “50” (late bolting). RNA-Seq and quantitative reverse transcription PCR (qPCR) analyses showed that two positive nitric oxide (NO) signaling regulator genes, nitrite reductase (BrNIR) and nitrate reductase (BrNIA), were up-regulated in line “50” with or without vernalization. In agreement with the transcription analysis, the shoots in line “50” had substantially higher nitrogen levels than those in “4004”. Upon vernalization, the flowering repressor gene Circadian 1 (BrCIR1) was significantly up-regulated in line “50”, whereas the flowering enhancer genes named SUPPRESSOR OF OVEREXPRESSION OF CONSTANCE 1 homologs (BrSOC1s) were substantially up-regulated in line “4004”. CRISPR/Cas9-mediated mutagenesis in Chinese cabbage demonstrated that the BrSOC1-1/1-2/1-3 genes were involved in late flowering, and their expression was mutually exclusive with that of the nitrogen signaling genes. Thus, we identified two flowering mechanisms in Chinese cabbage: a reciprocal negative feedback loop between nitrogen signaling genes (BrNIA1 and BrNIR1) and BrSOC1s to control flowering time and positive feedback control of the expression of BrSOC1s.

Comparative transcriptomics reveals desynchronisation of gene expression during the floral transition between Arabidopsis and Brassica rapa cultivars

Comparative transcriptomics can be used to translate an understanding of gene regulatory networks from model systems to less studied species. Here, we use RNA-Seq to determine and compare gene expression dynamics through the floral transition in the model species Arabidopsis thaliana and the closely related crop Brassica rapa. We find that different curve registration functions are required for different genes, indicating that there is no single common 'developmental time' between Arabidopsis and B. rapa. A detailed comparison between Arabidopsis and B. rapa and between two B. rapa accessions reveals different modes of regulation of the key floral integrator SOC1, and that the floral transition in the B. rapa accessions is triggered by different pathways. Our study adds to the mechanistic understanding of the regulatory network of flowering time in rapid cycling B. rapa and highlights the importance of registration methods for the comparison of developmental gene expression data.

A novel deletion in FLOWERING LOCUS T modulates flowering time in winter oilseed rape

A novel structural variant was discovered in the FLOWERING LOCUS T orthologue BnaFT.A02 by long-read sequencing. Nested association mapping in an elite winter oilseed rape population revealed that this 288 bp deletion associates with early flowering, putatively by modification of binding-sites for important flowering regulation genes. Perfect timing of flowering is crucial for optimal pollination and high seed yield. Extensive previous studies of flowering behavior in Brassica napus (canola, rapeseed) identified mutations in key flowering regulators which differentiate winter, semi-winter and spring ecotypes. However, because these are generally fixed in locally adapted genotypes, they have only limited relevance for fine adjustment of flowering time in elite cultivar gene pools. In crosses between ecotypes, the ecotype-specific major-effect mutations mask minor-effect loci of interest for breeding. Here, we investigated flowering time in a multiparental mapping population derived from seven elite winter oilseed rape cultivars which are fixed for major-effect mutations separating winter-type rapeseed from other ecotypes. Association mapping revealed eight genomic regions on chromosomes A02, C02 and C03 associating with fine modulation of flowering time. Long-read genomic resequencing of the seven parental lines identified seven structural variants coinciding with candidate genes for flowering time within chromosome regions associated with flowering time. Segregation patterns for these variants in the elite multiparental population and a diversity set of winter types using locus-specific assays revealed significant associations with flowering time for three deletions on chromosome A02. One of these was a previously undescribed 288 bp deletion within the second intron of FLOWERING LOCUS T on chromosome A02, emphasizing the advantage of long-read sequencing for detection of structural variants in this size range. Detailed analysis revealed the impact of this specific deletion on flowering-time modulation under extreme environments and varying day lengths in elite, winter-type oilseed rape.

Comparative Transcriptome Analysis of Early- and Late-Bolting Traits in Chinese Cabbage (Brassica rapa)

Chinese cabbage is one of the most important and widely consumed vegetables in China. The developmental transition from the vegetative to reproductive phase is a crucial process in the life cycle of flowering plants. In spring-sown Chinese cabbage, late bolting is desirable over early bolting. In this study, we analyzed double haploid (DH) lines of late bolting (“Y410-1” and “SY2004”) heading Chinese cabbage (Brassica rapa var. pekinensis) and early-bolting Chinese cabbage (“CX14-1”) (B. rapa ssp. chinensis var. parachinensis) by comparative transcriptome profiling using the Illumina RNA-seq platform. We assembled 721.49 million clean high-quality paired-end reads into 47,363 transcripts and 47,363 genes, including 3,144 novel unigenes. There were 12,932, 4,732, and 4,732 differentially expressed genes (DEGs) in pairwise comparisons of Y410-1 vs. CX14-1, SY2004 vs. CX14-1, and Y410-1 vs. SY2004, respectively. The RNA-seq results were confirmed by reverse transcription quantitative real-time PCR (RT-qPCR). A Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis of DEGs revealed significant enrichment for plant hormone and signal transduction as well as starch and sucrose metabolism pathways. Among DEGs related to plant hormone and signal transduction, six unigenes encoding the indole-3-acetic acid-induced protein ARG7 (BraA02g009130), auxin-responsive protein SAUR41 (BraA09g058230), serine/threonine-protein kinase BSK11 (BraA07g032960), auxin-induced protein 15A (BraA10g019860), and abscisic acid receptor PYR1 (BraA08g012630 and BraA01g009450), were upregulated in both late bolting Chinese cabbage lines (Y410-1 and SY2004) and were identified as putative candidates for the trait. These results improve our understanding of the molecular mechanisms underlying flowering in Chinese cabbage and provide a foundation for studies of this key trait in related species.

Genetics of days to flowering, maturity and plant height in natural and derived forms of Brassica rapa L

Genome wide association studies enabled prediction of many candidate genes for flowering, maturity and plant height under differing day-length conditions. Some genes were envisaged only from derived B. rapa. Flowering and plant height are the key life history traits. These are crucial for adaptation and productivity. Current investigations aimed to examine genotypic differences governing days to flowering, maturity and plant height under contrasting day-length conditions; and identify genomic regions governing the observed phenotypic variations. An association panel comprising 195 inbred lines, representing natural (NR) and derived (DR) forms of Brassica rapa (AA; 2n = 20), was evaluated at two sowing dates and two locations, representing different day-length regimes. Derived B. rapa is a unique pre-breeding material extracted from B. juncea (AABB; 2n = 36). Population structure analysis, using DArT genotypes established derived B. rapa as a genetic resource distinct from natural B. rapa. Genome wide association studies facilitated detection of many trait associated SNPs. Chromosomes A03, A05 and A09 harboured majority of these. Functional annotation of the associated SNPs and surrounding genome space(s) helped to predict 43 candidate genes. Many of these were predicted under specific day-length conditions. Important among these were the genes encoding floral meristem identity (SPL3, SPL15, AP3, BAM2), photoperiodic responses (COL2, AGL18, SPT, NF-YC4), gibberellic acid biosynthesis (GA1) and regulation of flowering (EBS). Some of the predicted genes were detected for DR subpanel alone. Genes controlling hormones, auxins and gibberellins appeared important for the regulation of plant height. Many of the significant SNPs were located on chromosomes harbouring previously reported QTLs and candidate genes. The identified loci may be used for marker-assisted selection after due validation.

QTL-seq identifies BnaFT.A02 and BnaFLC.A02 as candidates for variation in vernalization requirement and response in winter oilseed rape (Brassica napus)

Winter, spring and biennial varieties of Brassica napus that vary in vernalization requirement are grown for vegetable and oil production. Here, we show that the obligate or facultative nature of the vernalization requirement in European winter oilseed rape is determined by allelic variation at a 10 Mbp region on chromosome A02. This region includes orthologues of the key floral regulators FLOWERING LOCUS C (BnaFLC.A02) and FLOWERING LOCUS T (BnaFT.A02). Polymorphism at BnaFLC.A02 and BnaFT.A02, mostly in cis-regulatory regions, results in distinct gene expression dynamics in response to vernalization treatment. Our data suggest allelic variation at BnaFT.A02 is associated with flowering time in the absence of vernalization, while variation at BnaFLC.A02 is associated with flowering time under vernalizing conditions. We hypothesize selection for BnaFLC.A02 and BnaFT.A02 gene expression variation has facilitated the generation of European winter oilseed rape varieties that are adapted to different winter climates. This knowledge will allow for the selection of alleles of flowering time regulators that alter the vernalization requirement of oilseed rape, informing the generation of new varieties with adapted flowering times and improved yields.

Alternative splicing of flowering time gene FT is associated with halving of time to flowering in coconut

Coconut palm has two distinct types-"tall" and "dwarf"-which differ morphologically. Tall coconut varieties need 8-10 years to start flowering, while dwarf coconut varieties only require 3-5 years. We compared seedling and reproductive stage transcriptomes for both coconut types to determine potential molecular mechanisms underlying control of flowering time in coconut. Several key genes in the photoperiod pathway were differentially expressed between seedling and reproductive leaf samples in both tall and dwarf coconut. These genes included suppressor of overexpression of constans (SOC1), flowering locus T (FT), and Apetala 1 (AP1). Alternative splicing analysis of genes in the photoperiod pathway further revealed that the FT gene produces different transcripts in tall compared to dwarf coconut. The shorter alternative splice variant of FT [which included a 6 bp deletion, alternative 3' splicing sites (A3SS)] was found to be exclusively present in dwarf coconut varieties but absent in most tall coconut varieties. Our results provide a valuable information resource as well as suggesting a probable mechanism for differentiation of flowering time onset in coconut, providing a target for future breeding work in accelerating time to flowering in this crop species.

The role of FRIGIDA and FLOWERING LOCUS C genes in flowering time of Brassica rapa leafy vegetables

There is a wide variation of flowering time among lines of Brassica rapa L. Most B. rapa leafy (Chinese cabbage etc.) or root (turnip) vegetables require prolonged cold exposure for flowering, known as vernalization. Premature bolting caused by low temperature leads to a reduction in the yield/quality of these B. rapa vegetables. Therefore, high bolting resistance is an important breeding trait, and understanding the molecular mechanism of vernalization is necessary to achieve this goal. In this study, we demonstrated that BrFRIb functions as an activator of BrFLC in B. rapa. We showed a positive correlation between the steady state expression levels of the sum of the BrFLC paralogs and the days to flowering after four weeks of cold treatment, suggesting that this is an indicator of the vernalization requirement. We indicate that BrFLCs are repressed by the accumulation of H3K27me3 and that the spreading of H3K27me3 promotes stable FLC repression. However, there was no clear relationship between the level of H3K27me3 in the BrFLC and the vernalization requirement. We also showed that if there was a high vernalization requirement, the rate of repression of BrFLC1 expression following prolonged cold treatments was lower.

Joint QTL mapping and transcriptome sequencing analysis reveal candidate flowering time genes in Brassica napus L

Background

Optimum flowering time is a key agronomic trait in Brassica napus. To investigate the genetic architecture and genetic regulation of flowering time in this important crop, we conducted quantitative trait loci (QTL) analysis of flowering time in a recombinant inbred line (RIL) population, including lines with extreme differences in flowering time, in six environments, along with RNA-Seq analysis.

Results

We detected 27 QTLs distributed on eight chromosomes among six environments, including one major QTL on chromosome C02 that explained 11–25% of the phenotypic variation and was stably detected in all six environments. RNA-Seq analysis revealed 105 flowering time-related differentially expressed genes (DEGs) that play roles in the circadian clock/photoperiod, autonomous pathway, and hormone and vernalization pathways. We focused on DEGs related to the regulation of flowering time, especially DEGs in QTL regions.

Conclusions

We identified 45 flowering time-related genes in these QTL regions, eight of which are DEGs, including key flowering time genes PSEUDO RESPONSE REGULATOR 7 (PRR7) and FY (located in a major QTL region on C02). These findings provide insights into the genetic architecture of flowering time in B. napus.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5356-8) contains supplementary material, which is available to authorized users.

Whole-Genome Resequencing of a Worldwide Collection of Rapeseed Accessions Reveals the Genetic Basis of Ecotype Divergence

Rapeseed (Brassica napus), an important oilseed crop, has adapted to diverse climate zones and latitudes by forming three main ecotype groups, namely winter, semi-winter, and spring types. However, genetic variations underlying the divergence of these ecotypes are largely unknown. Here, we report the global pattern of genetic polymorphisms in rapeseed determined by resequencing a worldwide collection of 991 germplasm accessions. A total of 5.56 and 5.53 million single-nucleotide polymorphisms (SNPs) as well as 1.86 and 1.92 million InDels were identified by mapping reads to the reference genomes of "Darmor-bzh" and "Tapidor," respectively. We generated a map of allelic drift paths that shows splits and mixtures of the main populations, and revealed an asymmetric evolution of the two subgenomes of B. napus by calculating the genetic diversity and linkage disequilibrium parameters. Selective-sweep analysis revealed genetic changes in genes orthologous to those regulating various aspects of plant development and response to stresses. A genome-wide association study identified SNPs in the promoter regions of FLOWERING LOCUS T and FLOWERING LOCUS C orthologs that corresponded to the different rapeseed ecotype groups. Our study provides important insights into the genomic footprints of rapeseed evolution and flowering-time divergence among three ecotype groups, and will facilitate screening of molecular markers for accelerating rapeseed breeding.

Translating Flowering Time From Arabidopsis thaliana to Brassicaceae and Asteraceae Crop Species

Flowering and seed set are essential for plant species to survive, hence plants need to adapt to highly variable environments to flower in the most favorable conditions. Endogenous cues such as plant age and hormones coordinate with the environmental cues like temperature and day length to determine optimal time for the transition from vegetative to reproductive growth. In a breeding context, controlling flowering time would help to speed up the production of new hybrids and produce high yield throughout the year. The flowering time genetic network is extensively studied in the plant model species Arabidopsis thaliana, however this knowledge is still limited in most crops. This article reviews evidence of conservation and divergence of flowering time regulation in A. thaliana with its related crop species in the Brassicaceae and with more distant vegetable crops within the Asteraceae family. Despite the overall conservation of most flowering time pathways in these families, many genes controlling this trait remain elusive, and the function of most Arabidopsis homologs in these crops are yet to be determined. However, the knowledge gathered so far in both model and crop species can be already exploited in vegetable crop breeding for flowering time control.

Genome-Wide Identification of Flowering-Time Genes in Brassica Species and Reveals a Correlation between Selective Pressure and Expression Patterns of Vernalization-Pathway Genes in Brassica napus

Flowering time is a key agronomic trait, directly influencing crop yield and quality. Many flowering-time genes have been identified and characterized in the model plant Arabidopsis thaliana; however, these genes remain uncharacterized in many agronomically important Brassica crops. In this study, we identified 1064, 510, and 524 putative orthologs of A. thaliana flowering-time genes from Brassica napus, Brassica rapa, and Brassica oleracea, respectively, and found that genes involved in the aging and ambient temperature pathways were fewer than those in other flowering pathways. Flowering-time genes were distributed mostly on chromosome C03 in B. napus and B. oleracea, and on chromosome A09 in B. rapa. Calculation of non-synonymous (Ka)/synonymous substitution (Ks) ratios suggested that flowering-time genes in vernalization pathways experienced higher selection pressure than those in other pathways. Expression analysis showed that most vernalization-pathway genes were expressed in flowering organs. Approximately 40% of these genes were highly expressed in the anther, whereas flowering-time integrator genes were expressed in a highly organ-specific manner. Evolutionary selection pressures were negatively correlated with the breadth and expression levels of vernalization-pathway genes. These findings provide an integrated framework of flowering-time genes in these three Brassica crops and provide a foundation for deciphering the relationship between gene expression patterns and their evolutionary selection pressures in Brassica napus.

A Genomic Variation Map Provides Insights into the Genetic Basis of Spring Chinese Cabbage (Brassica rapa ssp. pekinensis) Selection

Chinese cabbage is the most consumed leafy crop in East Asian countries. However, premature bolting induced by continuous low temperatures severely decreases the yield and quality of the Chinese cabbage, and therefore restricts its planting season and geographic distribution. In the past 40 years, spring Chinese cabbage with strong winterness has been selected to meet the market demand. Here, we report a genome variation map of Chinese cabbage generated from the resequencing data of 194 geographically diverse accessions of three ecotypes. In-depth analyses of the selection sweeps and genome-wide patterns revealed that spring Chinese cabbage was selected from a specific population of autumn Chinese cabbage around the area of Shandong peninsula in northern China. We identified 23 genomic loci that underwent intensive selection, and further demonstrated by gene expression and haplotype analyses that the incorporation of elite alleles of VERNALISATION INSENTIVE 3.1 (BrVIN3.1) and FLOWER LOCUS C 1 (BrFLC1) is a determinant genetic source of variation during selection. Moreover, we showed that the quantitative response of BrVIN3.1 to cold due to the sequence variations in the cis elements of the BrVIN3.1 promoter significantly contributes to bolting-time variation in Chinese cabbage. Collectively, our study provides valuable insights into the genetic basis of spring Chinese cabbage selection and will facilitate the breeding of bolting-resistant varieties by molecular-marker-assisted selection, transgenic or gene editing approaches.

Subtropical adaptation of a temperate plant (Brassica oleracea var. italica) utilizes non-vernalization-responsive QTLs

While many tropical plants have been adapted to temperate cultivation, few temperate plants have been adapted to the tropics. Originating in Western Europe, Brassica oleracea vernalization requires a period of low temperature and BoFLC2 regulates the transition to floral development. In B. oleracea germplasm selected in Taiwan, a non-vernalization pathway involving BoFLC3 rather than BoFLC2 regulates curd induction. In 112 subtropical breeding lines, specific haplotype combinations of BoFLC3 and PAN (involved in floral organ identity and a positional candidate for additional curd induction variation) adapt B. oleracea to high ambient temperature and short daylength. Duplicated genes permitted evolution of alternative pathways for control of flowering in temperate and tropical environments, a principle that might be utilized via natural or engineered approaches in other plants. New insight into regulation of Brassica flowering exemplifies translational agriculture, tapping knowledge of botanical models to improve food security under projected climate change scenarios.