Differential regulation of FLOWERING LOCUS C expression by vernalization in cabbage and Arabidopsis

Evolution of major flowering pathway integrators in Orchidaceae

The Orchidaceae is a mega-diverse plant family with ca. 29,000 species with a large variety of life forms that can colonize transitory habitats. Despite this diversity, little is known about their flowering integrators in response to specific environmental factors. During the reproductive transition in flowering plants a vegetative apical meristem (SAM) transforms into an inflorescence meristem (IM) that forms bracts and flowers. In model grasses, like rice, a flowering genetic regulatory network (FGRN) controlling reproductive transitions has been identified, but little is known in the Orchidaceae. In order to analyze the players of the FRGN in orchids, we performed comprehensive phylogenetic analyses of CONSTANS-like/CONSTANS-like 4 (COL/COL4), FLOWERING LOCUS D (FD), FLOWERING LOCUS C/FRUITFULL (FLC/FUL) and SUPRESSOR OF OVEREXPRESSION OF CONSTANS 1 (SOC1) gene lineages. In addition to PEBP and AGL24/SVP genes previously analyzed, here we identify an increase of orchid homologs belonging to COL4, and FUL gene lineages in comparison with other monocots, including grasses, due to orchid-specific gene lineage duplications. Contrariwise, local duplications in Orchidaceae are less frequent in the COL, FD and SOC1 gene lineages, which points to a retention of key functions under strong purifying selection in essential signaling factors. We also identified changes in the protein sequences after such duplications, variation in the evolutionary rates of resulting paralogous clades and targeted expression of isolated homologs in different orchids. Interestingly, vernalization-response genes like VERNALIZATION1 (VRN1) and FLOWERING LOCUS C (FLC) are completely lacking in orchids, or alternatively are reduced in number, as is the case of VERNALIZATION2/GHD7 (VRN2). Our findings point to non-canonical factors sensing temperature changes in orchids during reproductive transition. Expression data of key factors gathered from Elleanthus auratiacus, a terrestrial orchid in high Andean mountains allow us to characterize which copies are actually active during flowering. Altogether, our data lays down a comprehensive framework to assess gene function of a restricted number of homologs identified more likely playing key roles during the flowering transition, and the changes of the FGRN in neotropical orchids in comparison with temperate grasses.

Molecular Evaluation of the Effects of FLC Homologs and Coordinating Regulators on the Flowering Responses to Vernalization in Cabbage (Brassica oleracea var. capitata) Genotypes

The flowering loci of cabbage must be understood to boost their productivity. In this study, to clarify the flowering mechanisms of cabbage, we examined the three flowering repressors BoFLC1, 2 and 3, and the flowering regulators BoGI, BoCOOLAIR, and BoVIN3 of early (CAB1), middle (CAB3), and late (CAB5) flowering cabbage genotypes. Analysis of allele-specifically amplified genomic DNA and various sequence alignments demonstrated that maximal insertions and deletions influenced cabbage flowering behavior, notably in CAB3 and CAB5. Phylogenetic studies showed that BoFLC1, 2, and 3 in the CAB1, 3, and 5 genotypes had the highest homologies to other Brassica species, with CAB3 and 5 the most similar. Although CAB3 and CAB5 have comparable genetic patterns, flowering repressors and flowering regulators were investigated individually with and without vernalization to determine their minor flowering differences. The expression investigation revealed that vernalized CAB5 downregulated all BoFLC genes compared to CAB3 and, in contrast, CAB3 exhibited upregulated BoCOOLAIR. We hypothesized that the CAB3 BoFLC locus' additional insertions may have led to BoCOOLAIR overexpression and BoFLC downregulation. This study sheds light on cabbage genotypes-particularly those of CAB1 and CAB5-and suggests that structural variations in BoFLC2 and 3 bind flowering regulators, such as COOLAIR, which may affect cabbage flowering time.

Identification of vernalization-related genes and cold memory element (CME) required for vernalization response in radish (Raphanus sativus L.)

Floral transition is accelerated by exposure to long-term cold like winter in plants, which is called as vernalization. Acceleration of floral transition by vernalization is observed in a diversity of biennial and perennial plants including Brassicaceae family plants. Scientific efforts to understand molecular mechanism underlying vernalization-mediated floral transition have been intensively focused in model plant Arabidopsis thaliana. To get a better understanding on floral transition by vernalization in radish (Raphanus sativus L.), we investigated transcriptomic changes taking place during vernalization in radish. Thousands of genes were differentially regulated along time course of vernalization compared to non-vernalization (NV) sample. Twelve major clusters of DEGs were identified based on distinctive expression profiles during vernalization. Radish FLC homologs were shown to exert an inhibition of floral transition when transformed into Arabidopsis plants. In addition, DNA region containing RY motifs located within a Raphanus sativus FLC homolog, RsFLC1 was found to be required for repression of RsFLC1 by vernalization. Transgenic plants harboring disrupted RY motifs were impaired in the enrichment of H3K27me3 on RsFLC1 chromatin, thus resulting in the delayed flowering in Arabidopsis. Taken together, we report transcriptomic profiles of radish during vernalization and demonstrate the requirement of RY motif for vernalization-mediated repression of RsFLC homologs in radish (Raphanus sativus L.).

Genetic dissection of morphological variation between cauliflower and a rapid cycling Brassica oleracea line

To improve resolution to small genomic regions and sensitivity to small-effect loci in the identification of genetic factors conferring the enlarged inflorescence and other traits of cauliflower while also expediting further genetic dissection, 104 near-isogenic introgression lines (NIILs) covering 78.56% of the cauliflower genome, were selected from an advanced backcross population using cauliflower [Brassica oleracea var. botrytis L., mutant for Orange gene (ORG)] as the donor parent and a rapid cycling line (TO1434) as recurrent parent. Subsets of the advanced backcross population and NIILs were planted in the field for 8 seasons, finding 141 marker-trait associations for 15 leaf-, stem-, and flower-traits. Exemplifying the usefulness of these lines, we delineated the previously known flower color gene to a 4.5 MB interval on C3; a gene for small plant size to a 3.4 MB region on C8; and a gene for large plant size and flowering time to a 6.1 MB region on C9. This approach unmasked closely linked QTL alleles with opposing effects (on chr. 8) and revealed both alleles with expected phenotypic effects and effects opposite the parental phenotypes. Selected B. oleracea NIILs with short generation time add new value to widely used research and teaching materials.

Upregulation of tandem duplicated BoFLC1 genes is associated with the non-flowering trait in Brassica oleracea var. capitata

Tandem duplicated BoFLC1 genes (BoFLC1a and BoFLC1b), which were identified as the candidate causal genes for the non-flowering trait in the cabbage mutant 'nfc', were upregulated during winter in 'nfc'. The non-flowering natural cabbage mutant 'nfc' was discovered from the breeding line 'T15' with normal flowering characteristics. In this study, we investigated the molecular basis underlying the non-flowering trait of 'nfc'. First, 'nfc' was induced to flower using the grafting floral induction method, and three F₂ populations were generated. The flowering phenotype of each F₂ population was widely distributed with non-flowering individuals appearing in two populations. QTL-seq analysis detected a genomic region associated with flowering date at approximately 51 Mb on chromosome 9 in two of the three F₂ populations. Subsequent validation and fine mapping of the candidate genomic region using QTL analysis identified the quantitative trait loci (QTL) at 50,177,696-51,474,818 bp on chromosome 9 covering 241 genes. Additionally, RNA-seq analysis in leaves and shoot apices of 'nfc' and 'T15' plants identified 19 and 15 differentially expressed genes related to flowering time, respectively. Based on these results, we identified tandem duplicated BoFLC1 genes, which are homologs of the floral repressor FLOWERING LOCUS C, as the candidate genes responsible for the non-flowering trait of 'nfc'. We designated the tandem duplicated BoFLC1 genes as BoFLC1a and BoFLC1b. Expression analysis revealed that the expression levels of BoFLC1a and BoFLC1b were downregulated during winter in 'T15' but were upregulated and maintained during winter in 'nfc'. Additionally, the expression level of the floral integrator BoFT was upregulated in the spring in 'T15' but hardly upregulated in 'nfc'. These results suggest that the upregulated levels of BoFLC1a and BoFLC1b contributed to the non-flowering trait of 'nfc'.

Brassica rapa orphan gene BR1 delays flowering time in Arabidopsis

Orphan genes are essential to the emergence of species-specific traits and the process of evolution, lacking sequence similarity to any other identified genes. As they lack recognizable domains or functional motifs, however, efforts to characterize these orphan genes are often difficult. Flowering is a key trait in Brassica rapa, as premature bolting can have a pronounced adverse impact on plant quality and yield. Bolting resistance-related orphan genes, however, have yet to be characterized. In this study, an orphan gene designated BOLTING RESISTANCE 1 (BR1) was identified and found through gene structural variation analyses to be more highly conserved in Chinese cabbage than in other available accessions. The expression of BR1 was increased in bolting resistant Chinese cabbage and decreased in bolting non-resistant type, and the expression of some mark genes were consist with bolting resistance phenotype. BR1 is primarily expressed in leaves at the vegetative growth stage, and the highest BR1 expression levels during the flowering stage were observed in the flower buds and silique as compared to other tissue types. The overexpression of BR1 in Arabidopsis was associated with enhanced bolting resistance under long day (LD) conditions, with these transgenic plants exhibiting significant decreases in stem height, rosette radius, and chlorophyll content. Transcriptomic sequencing of WT and BR1OE plants showed the association of BR1 with other bolting resistance genes. Transcriptomic sequencing and qPCR revealed that six flowering integrator genes and one chlorophyll biosynthesis-related gene were downregulated following BR1 overexpression. Six key genes in photoperiodic flowering pathway exhibited downward expression trends in BR1OE plants, while the expression of floral repressor AtFLC gene was upregulated. The transcripts of these key genes were consistent with observed phenotypes in BR1OE plants, and the results indicated that BR1 may function through vernalization and photoperiodic pathway. Instead, the protein encoded by BR1 gene was subsequently found to localize to the nucleus. Taken together, we first propose that orphan gene BR1 functions as a novel regulator of flowering time, and these results suggested that BR1 may represent a promising candidate gene to support the selective breeding of Chinese cabbage cultivars with enhanced bolting resistance.

A 215-bp indel at intron I of BoFLC2 affects flowering time in Brassica oleracea var. capitata during vernalization

In response to cold, a 215-bp deletion at intron I of BoFLC2 slows its silencing activity by feedback to the core genes of the PHD-PRC2 complex, resulting in late flowering in cabbage. Cabbage is a plant-vernalization-responsive flowering type. In response to cold, BoFLC2 is an important transcription factor, which allows cabbage plants to remain in the vegetative phase. However, there have been few reports on the detailed and functional effects of genetic variation in BoFLC2 on flowering time in cabbage. Herein, BoFLC2^E and BoFLC2^L, cloned from extremely early and extremely late flowering cabbages, respectively, exhibited a 215-bp indel at intron I, three non-synonymous SNPs and a 3-bp indel at exon II. BoFLC2^L was found to be related to late flowering, as verified in 40 extremely early/late flowering accessions, a diverse set of cabbage inbred lines and two F2 generations by using indel-FLC2 marker. Among the genetic variation of BoFLC2, the 215-bp deletion at intron I was the main reason for the delayed flowering time, as verified in the transgenic progenies of seed-vernalization-responsive Arabidopsis thaliana (Col) and rapid cycler B. oleracea (TO1000, boflc2). This is the first report to show that the intron I indel of BoFLC2 affects the flowering time of cabbage. Although the intron I 215-bp indel between BoFLC2^E and BoFLC2^L did not cause alternative splicing, it slowed BoFLC2^L silencing during vernalization and feedback to the core genes of the PHD-PRC2 complex, resulting in their lower transcription levels. Our study not only provides an effective molecular marker-assisted selective strategy for identifying bolting-resistant resources and breeding improved varieties in cabbage, but also provides an entry point for exploring the mechanisms of flowering time in plant-vernalization-responsive plants.

Non-vernalization requirement in Chinese kale caused by loss of BoFLC and low expressions of its paralogs

We identified the loss of BoFLC gene as the cause of non-vernalization requirement in B. oleracea. Our developed codominant marker of BoFLC gene can be used for breeding program of B. oleracea crops. Many species of the Brassicaceae family, including some Brassica crops, require vernalization to avoid pre-winter flowering. Vernalization is an unfavorable trait for Chinese kale (Brassica oleracea var. chinensis Lei), a stem vegetable, and therefore it has been lost during its domestication/breeding process. To reveal the genetics of vernalization variation, we constructed an F₂ population through crossing a Chinese kale (a non-vernalization crop) with a kale (a vernalization crop). Using bulked segregant analysis (BSA) and RNA-seq, we identified one major quantitative trait locus (QTL) controlling vernalization and fine-mapped it to a region spanning 80 kb. Synteny analysis and PCR-based sequencing results revealed that compared to that of the kale parent, the candidate region of the Chinese kale parent lost a 9,325-bp fragment containing FLC homolog (BoFLC). In addition to the BoFLC gene, there are four other FLC homologs in the genome of B. oleracea, including Bo3g005470, Bo3g024250, Bo9g173370, and Bo9g173400. The qPCR analysis showed that the BoFLC had the highest expression among the five members of the FLC family. Considering the low expression levels of the four paralogs of BoFLC, we speculate that its paralogs cannot compensate the function of the lost BoFLC, therefore the presence/absence (PA) polymorphism of BoFLC determines the vernalization variation. Based on the PA polymorphism of BoFLC, we designed a codominant marker for the vernalization trait, which can be used for breeding programs of B. oleracea crops.

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.

Gene co-expression network analysis reveals key pathways and hub genes in Chinese cabbage (Brassica rapa L.) during vernalization

BACKGROUND

Vernalization is a type of low temperature stress used to promote rapid bolting and flowering in plants. Although rapid bolting and flowering promote the reproduction of Chinese cabbages (Brassica rapa L. ssp. pekinensis), this process causes their commercial value to decline. Clarifying the mechanisms of vernalization is essential for its further application. We performed RNA sequencing of gradient-vernalization in order to explore the reasons for the different bolting process of two Chinese cabbage accessions during vernalization.

RESULTS

There was considerable variation in gene expression between different-bolting Chinese cabbage accessions during vernalization. Comparative transcriptome analysis and weighted gene co-expression network analysis (WGCNA) were performed for different-bolting Chinese cabbage during different vernalization periods. The biological function analysis and hub gene annotation of highly relevant modules revealed that shoot system morphogenesis and polysaccharide and sugar metabolism caused early-bolting 'XBJ' to bolt and flower faster; chitin, ABA and ethylene-activated signaling pathways were enriched in late-bolting 'JWW'; and leaf senescence and carbohydrate metabolism enrichment were found in the two Chinese cabbage-related modules, indicating that these pathways may be related to bolting and flowering. The high connectivity of hub genes regulated vernalization, including MTHFR2, CPRD49, AAP8, endoglucanase 10, BXLs, GATLs, and WRKYs. Additionally, five genes related to flower development, BBX32 (binds to the FT promoter), SUS1 (increases FT expression), TSF (the closest homologue of FT), PAO and NAC029 (plays a role in leaf senescence), were expressed in the two Chinese cabbage accessions.

CONCLUSION

The present work provides a comprehensive overview of vernalization-related gene networks in two different-bolting Chinese cabbages during vernalization. In addition, the candidate pathways and hub genes related to vernalization identified here will serve as a reference for breeders in the regulation of Chinese cabbage production.

Total FLC transcript dynamics from divergent paralogue expression explains flowering diversity in Brassica napus

Flowering time is a key adaptive and agronomic trait. In Arabidopsis, natural variation in expression levels of the floral repressor FLOWERING LOCUS C (FLC) leads to differences in vernalization. In Brassica napus there are nine copies of FLC. Here, we study how these multiple FLC paralogues determine vernalization requirement as a system. We collected transcriptome time series for Brassica napus spring, winter, semi-winter, and Siberian kale crop types. Modelling was used to link FLC expression dynamics to floral response following vernalization. We show that relaxed selection pressure has allowed expression of FLC paralogues to diverge, resulting in variation of FLC expression during cold treatment between paralogues and accessions. We find that total FLC expression dynamics best explains differences in cold requirement between cultivars, rather than expression of specific FLC paralogues. The combination of multiple FLC paralogues with different expression dynamics leads to rich behaviour in response to cold and a wide range of vernalization requirements in B. napus. We find evidence for different strategies to determine the response to cold in existing winter rapeseed accessions.

Repressive chromatin modification underpins the long-term expression trend of a perennial flowering gene in nature

Natural environments require organisms to possess robust mechanisms allowing responses to seasonal trends. In Arabidopsis halleri, the flowering regulator AhgFLC shows upregulation and downregulation phases along with long-term past temperature, but the underlying machinery remains elusive. Here, we investigate the seasonal dynamics of histone modifications, H3K27me3 and H3K4me3, at AhgFLC in a natural population. Our advanced modelling and transplant experiments reveal that H3K27me3-mediated chromatin regulation at AhgFLC provides two essential properties. One is the ability to respond to the long-term temperature trends via bidirectional interactions between H3K27me3 and H3K4me3; the other is the ratchet-like character of the AhgFLC system, i.e. reversible in the entire perennial life cycle but irreversible during the upregulation phase. Furthermore, we show that the long-term temperature trends are locally indexed at AhgFLC in the form of histone modifications. Our study provides a more comprehensive understanding of H3K27me3 function at AhgFLC in a complex natural environment.

The flowering regulator FLC shows upregulation and downregulation phases along with long-term past temperature in Arabidopsishalleri. Here, the authors reveal that H3K27me3-mediated chromatin regulation at AhgFLC provides the ability to respond to both the seasonal temperature trends and the perennial life cycle.

The histone modification H3 lysine 27 tri-methylation has conserved gene regulatory roles in the triplicated genome of Brassica rapa L

Brassica rapa L. is an important vegetable and oilseed crop. We investigated the distribution of the histone mark tri-methylation of H3K27 (H3K27me3) in B. rapa and its role in the control of gene expression at two stages of development (2-day cotyledons and 14-day leaves) and among paralogs in the triplicated genome. H3K27me3 has a similar distribution in two inbred lines, while there was variation of H3K27me3 sites between tissues. Sites that are specific to 2-day cotyledons have increased transcriptional activity, and low levels of H3K27me3 in the gene body region. In 14-day leaves, levels of H3K27me3 were associated with decreased gene expression. In the triplicated genome, H3K27me3 is associated with paralogs that have tissue-specific expression. Even though B. rapa and Arabidopsis thaliana are not closely related within the Brassicaceae, there is conservation of H3K27me3-marked sites in the two species. Both B. rapa and A. thaliana require vernalization for floral initiation with FLC being the major controlling locus. In all four BrFLC paralogs, low-temperature treatment increases H3K27me3 at the proximal nucleation site reducing BrFLC expression. Following return to normal temperature growth conditions, H3K27me3 spreads along all four BrFLC paralogs providing stable repression of the gene.

Genome-wide identification of lncRNAs during hickory (Carya cathayensis) flowering

Non-coding RNAs with lengths greater than 200 bp are known as long non-coding RNAs (lncRNAs), and these RNAs play important role in gene regulation and plant development. However, to date, little is known regarding the role played by lncRNAs during flowering in hickory (Carya cathayensis). Here, we performed whole transcriptome RNA-sequencing of samples from hickory female and male floral buds, in which three samples (H0311PF, H0318PF, and H0402PF) represent pre-flowering, flowering, and post-flowering, respectively, while eight male samples collected from May 8th to June 13th as this time course are the key stage for male floral bud differentiation. We identified 2163 lncRNAs in hickory during flowering, containing 213 intronic, 1488 intergenic, and 462 antisense lncRNAs. We noticed that 510 and 648 lncRNAs were differentially expressed corresponding to female and male floral buds, respectively. And some of the lncRNAs were in a tightly tissue-specific or stage-specific manner. To further understand the roles of the lncRNAs, we predicted the function of the lncRNAs in cis- and trans-acting modes. The results showed that 924 lncRNAs were cis-correlated with 1536 protein-coding genes, while 1207 lncRNAs co-expressed (trans-acting) with 7432 protein-coding genes (R > 0.95 or R < - 0.95). These lncRNAs were all enriched in flower development-associated biological processes, i.e., circadian rhythm, vernalization response, response to gibberellin, inflorescence development, floral organ development, etc. To further understand the relationships between lncRNAs and floral-core genes, we build a co-expressing lncRNA-mRNA flowering network. We classified these floral genes into different pathway (photoperiod, vernalization, gibberellin, autonomous, and sucrose pathway) according to their particular functions. We found a set of lncRNAs that preferentially expressed in these pathways. The network showed that some lncRNAs (i.e., XLOC_038669, XLOC_017938) functioned in a particular flowering time pathway, while others (i.e., XLOC_011251, XLOC_04110) were involved in multiple pathway. Furthermore, some lncRNAs (i.e., XLOC_038669, XLOC_009597, and XLOC_049539) played roles in single or multiple pathways via interaction with each other. This study provides a genome-wide survey of hickory flower-related lncRNAs and will contribute to further understanding of the molecular mechanism underpinning flowering in hickory.

Genetic dissection of flowering time in Brassica rapa responses to temperature and photoperiod

The Brassica rapa (B. rapa) species displays enormous phenotypic diversity, with leafy vegetables, storage root vegetables and oil crops. These different crops all have different flowering time, which determine their growing season and cultivation area. Little is known about the effects of diverse temperature and day-lengths on flowering time QTL associated with FLC paralogues. We phenotyped the flowering time of a doubled haploid population, established from a cross between Yellow sarson and Pak choi under diverse environmental conditions. We identified flowering-time QTL (fQTL) in different photoperiod and temperature regimes in the greenhouse, and studied their colocation with known flowering time genes. As several fQTL colocalized with FLC paralogues, we studied the expression patterns of four FLC paralogues during the course of vernalization in parental lines. Under all environmental conditions tested the major fQTL that mapped to the BrFLC2_A02 locus was detected, however its effect decreased when plants were grown at low temperatures. Another fQTL that mapped to the FLC paralogue, BrFLC5_A03 was also identified under all tested environments, while no fQTL colocated with BrFLC1_A10 or BrFLC3_A03. Furthermore, the vernalization treatment decreased expression of all BrFLC paralogues in the parental lines, and showed the lowest transcript level after 28 days of vernalization. Transcript abundance stayed low after returning the plants for seven days to normal growth temperature. Interestingly, transcript abundance of BrFLC3_A03 and BrFLC5_A03 was repressed much stronger and already reached lowest levels after 14d in the early-flowering type YS-143. This study improves understanding of the effects of daylength and vernalization on flowering time in B. rapa and the role of the different BrFLC paralogues therein.

Non-vernalization Flowering and Seed Set of Cabbage Induced by Grafting Onto Radish Rootstocks

Cabbage (Brassica oleracea var. capitata) requires a long-term low-temperature exposure for floral induction, causing a delay in the breeding cycle. The objective of this study is to develop a method to induce flowering in cabbage without low-temperature treatment, using a grafting method. We conducted grafting experiments using two flower-induced Chinese kale cultivars (B. oleracea var. alboglabra) and seven radish cultivars/accessions as rootstocks and investigated the flowering response of grafted cabbage scions without low-temperature treatment. "Watanabe-seiko No.1" cabbage, when grafted onto the two Chinese kale cultivars, did not formed flower buds. Flowering was successfully induced in "Watanabe-seiko No.1" by grafting onto three out of the seven tested radish cultivars, and in "Kinkei No.201" and "Red cabbage" by grafting onto one tested radish cultivar. In "Watanabe-seiko No.1," the earliest flower bud appearance was observed at 29 days after grafting. Seeds were also obtained from the three cabbage cultivars that flowered by grafting. Gene expression analysis of "Watanabe-seiko No.1" cabbage scions which formed flower buds by grafting, revealed high expression of the homolog of the floral integrator, SUPPRESSOR OF OVEREXPRESSION OF CONSTANS 1 (BoSOC1), at the time of flower bud appearance. However, in the same leaf samples, we observed low expression of two homologs of florigen, FLOWERING LOCUS T (BoFT.C2 and BoFT.C6). In addition, two homologs of the floral repressor FLOWERING LOCUS C (BoFLC3 and BoFLC4), which are known to be down-regulated before flower bud differentiation in the vernalization pathway, were highly expressed, indicating that grafting onto radish induces cabbage flowering independently of the vernalization pathway. The expression level of the radish FT homolog (RsFT) in "Rat's tail-G2," which had highly induced flowering in the grafted cabbage scion, was higher than in the other radish cultivars. However, although "Rat's tail-CH" effectively induced flowering in the cabbage scion, the expression of RsFT was low in this cultivar. In this study, floral induction of non-vernalized cabbage cannot be explained by the expression levels of RsFT in rootstock plants, alone. The flowering of non-vernalized cabbage would be induced by transmissible agents from rootstocks and not by the expression of cabbage FT, BoFT, from the scion itself.

Analysis of ambient temperature-responsive transcriptome in shoot apical meristem of heat-tolerant and heat-sensitive broccoli inbred lines during floral head formation

BACKGROUND

Head formation of broccoli (Brassica oleracea var. italica) is greatly reduced under high temperature (22 °C and 27 °C). Broccoli inbred lines that are capable of producing heads at high temperatures in summer are varieties that are unique to Taiwan. However, knowledge of the early-activated pathways of broccoli head formation under high temperature is limited.

RESULTS

We compared heat-tolerant (HT) and heat-sensitive (HS) transcriptome of broccoli under different temperatures. Weighted gene correlation network analysis (WGCNA) revealed that genes involved in calcium signaling pathways, mitogen-activated protein kinase (MAPK) cascades, leucine-rich repeat receptor-like kinases (LRR-RLKs), and genes coding for heat-shock proteins and reactive oxygen species homeostasis shared a similar expression pattern to BoFLC1, which was highly expressed at high temperature (27 °C). Of note, these genes were less expressed in HT than HS broccoli at 22 °C. Co-expression analysis identified a model for LRR-RLKs in survival-reproduction tradeoffs by modulating MAPK- versus phytohormones-signaling during head formation. The difference in head-forming ability in response to heat stress between HT and HS broccoli may result from their differential transcriptome profiles of LRR-RLK genes. High temperature induced JA- as well as suppressed auxin- and cytokinin-related pathways may facilitate a balancing act to ensure fitness at 27 °C. BoFLC1 was less expressed in HT than HS at 22 °C, whereas other FLC homologues were not. Promoter analysis of BoFLC1 showed fewer AT dinucleotide repeats in HT broccoli. These results provide insight into the early activation of stress- or development-related pathways during head formation in broccoli. The identification of the BoFLC1 DNA biomarker may facilitate breeding of HT broccoli.

CONCLUSIONS

In this study, HT and HS broccoli genotypes were used to determine the effect of temperature on head formation by transcriptome profiling. On the basis of the expression pattern of high temperature-associated signaling genes, the HS transcriptome may be involved in stress defense instead of transition to the reproductive phase in response to heat stress. Transcriptome profiling of HT and HS broccoli helps in understanding the molecular mechanisms underlying head-forming capacity and in promoting functional marker-assisted breeding.

Complex Horticultural Quality Traits in Broccoli Are Illuminated by Evaluation of the Immortal BolTBDH Mapping Population

Improving horticultural quality in regionally adapted broccoli (Brassica oleracea var. italica) and other B. oleracea crops is challenging due to complex genetic control of traits affecting morphology, development, and yield. Mapping horticultural quality traits to genomic loci is an essential step in these improvement efforts. Understanding the mechanisms underlying horticultural quality enables multi-trait marker-assisted selection for improved, resilient, and regionally adapted B. oleracea germplasm. The publicly-available biparental double-haploid BolTBDH mapping population (Chinese kale × broccoli; N = 175) was evaluated for 25 horticultural traits in six trait classes (architecture, biomass, phenology, leaf morphology, floral morphology, and head quality) by multiple quantitative trait loci mapping using 1,881 genotype-by-sequencing derived single nucleotide polymorphisms. The physical locations of 56 single and 41 epistatic quantitative trait locus (QTL) were identified. Four head quality QTL (OQ_C03@57.0, OQ_C04@33.3, OQ_CC08@25.5, and OQ_C09@49.7) explain a cumulative 81.9% of phenotypic variance in the broccoli heading phenotype, contain the FLOWERING LOCUS C (FLC) homologs Bo9g173400 and Bo9g173370, and exhibit epistatic effects. Three key genomic hotspots associated with pleiotropic control of the broccoli heading phenotype were identified. One phenology hotspot reduces days to flowering by 7.0 days and includes an additional FLC homolog Bo3g024250 that does not exhibit epistatic effects with the three horticultural quality hotspots. Strong candidates for other horticultural traits were identified: BoLMI1 (Bo3g002560) associated with serrated leaf margins and leaf apex shape, BoCCD4 (Bo3g158650) implicated in flower color, and BoAP2 (Bo1g004960) implicated in the hooked sepal horticultural trait. The BolTBDH population provides a framework for B. oleracea improvement by targeting key genomic loci contributing to high horticultural quality broccoli and enabling de novo mapping of currently unexplored traits.

Genetic dissection of the mechanism of flowering time based on an environmentally stable and specific QTL in Brassica napus

Flowering time is an important agronomic trait that is highly influenced by the environment. To elucidate the genetic mechanism of flowering time in rapeseed (Brassica napus L.), a genome-wide QTL analysis was performed in a doubled haploid population grown in winter, semi-winter and spring ecological conditions. Fifty-five consensus QTLs were identified after combining phenotype and genomic data, including 12 environment-stable QTLs and 43 environment-specific QTLs. Importantly, six major QTLs for flowering time were identified, of which two were considered environment-specific QTLs in spring ecological condition and four were considered environment-stable QTLs in winter and semi-winter ecological conditions. Through QTL comparison, 18 QTLs were colocalized with QTLs from six other published studies. Combining the candidate genes with their functional annotation, in 49 of 55 consensus QTLs, 151 candidate genes in B. napus corresponding to 95 homologous genes in Arabidopsis thaliana related to flowering were identified, including BnaC03g32910D (CO), BnaA02g12130D (FT) and BnaA03g13630D (FLC). Most of the candidate genes were involved in different flowering regulatory pathways. Based on re-sequencing and differences in sequence annotation between the two parents, we found that regions containing some candidate genes have numerous non-frameshift InDels and many non- synonymous mutations, which might directly lead to gene functional variation. Flowering time was negativly correlated with seed yield and thousand seed weight based on a QTL comparison of flowering time and seed yield traits, which has implications in breeding new early-maturing varieties of B. napus. Moreover, a putative flowering regulatory network was constructed, including the photoperiod, circadian clock, vernalization, autonomous and gibberellin pathways. Multiple copies of genes led to functional difference among the different copies of homologous genes, which also increased the complexity of the flowering regulatory networks. Taken together, the present results not only provide new insights into the genetic regulatory network underlying the control of flowering time but also improve our understanding of flowering time regulatory pathways in rapeseed.

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.

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.

Whole-transcriptome analysis reveals genetic factors underlying flowering time regulation in rapeseed (Brassica napus L.)

Rapeseed (Brassica napus L.), one of the most important sources of vegetable oil and protein-rich meals worldwide, is adapted to different geographical regions by modification of flowering time. Rapeseed cultivars have different day length and vernalization requirements, which categorize them into winter, spring, and semiwinter ecotypes. To gain a deeper insight into genetic factors controlling floral transition in B. napus, we performed RNA sequencing (RNA-seq) in the semiwinter doubled haploid line, Ningyou7, at different developmental stages and temperature regimes. The expression profiles of more than 54,000 gene models were compared between different treatments and developmental stages, and the differentially expressed genes were considered as targets for association analysis and genetic mapping to confirm their role in floral transition. Consequently, 36 genes with association to flowering time, seed yield, or both were identified. We found novel indications for neofunctionalization in homologs of known flowering time regulators like VIN3 and FUL. Our study proved the potential of RNA-seq along with association analysis and genetic mapping to identify candidate genes for floral transition in rapeseed. The candidate genes identified in this study could be subjected to genetic modification or targeted mutagenesis and genotype building to breed rapeseed adapted to certain environments.

The production and characterization of a BoFLC2 introgressed Brassica rapa by repeated backcrossing to an F1

Flowering time is an important agronomic trait for Brassica rapa crops, and previous breeding work in Brassica has successfully transmitted other important agronomic traits from donor species. However, there has been no previous attempts to produce hybrids replacing the original Brassica FLC alleles with alien FLC alleles. In this paper, we introduce the creation of a chromosome substitution line (CSSL) containing a homozygous introgression of Flowering Locus C from Brassica oleracea (BoFLC2) into a B. rapa genomic background, and characterize the CSSL line with respect to the parental cultivars. The preferential transmission of alien chromosome inheritance and the pattern of transmission observed during the production of the CSSLs are also discussed.

Arabidopsis inositol phosphate kinases IPK1 and ITPK1 constitute a metabolic pathway in maintaining phosphate homeostasis

Emerging studies have suggested that there is a close link between inositol phosphate (InsP) metabolism and cellular phosphate (P_i ) homeostasis in eukaryotes; however, whether a common InsP species is deployed as an evolutionarily conserved metabolic messenger to mediate P_i signaling remains unknown. Here, using genetics and InsP profiling combined with P_i -starvation response (PSR) analysis in Arabidopsis thaliana, we showed that the kinase activity of inositol pentakisphosphate 2-kinase (IPK1), an enzyme required for phytate (inositol hexakisphosphate; InsP₆ ) synthesis, is indispensable for maintaining P_i homeostasis under P_i -replete conditions, and inositol 1,3,4-trisphosphate 5/6-kinase 1 (ITPK1) plays an equivalent role. Although both ipk1-1 and itpk1 mutants exhibited decreased levels of InsP₆ and diphosphoinositol pentakisphosphate (PP-InsP₅ ; InsP₇ ), disruption of another ITPK family enzyme, ITPK4, which correspondingly caused depletion of InsP₆ and InsP₇ , did not display similar P_i -related phenotypes, which precludes these InsP species from being effectors. Notably, the level of d/l-Ins(3,4,5,6)P₄ was concurrently elevated in both ipk1-1 and itpk1 mutants, which showed a specific correlation with the misregulated P_i phenotypes. However, the level of d/l-Ins(3,4,5,6)P₄ is not responsive to P_i starvation that instead manifests a shoot-specific increase in the InsP₇ level. This study demonstrates a more nuanced picture of the intersection of InsP metabolism and P_i homeostasis and PSRs than has previously been elaborated, and additionally establishes intermediate steps to phytate biosynthesis in plant vegetative tissues.

Arabidopsis inositol phosphate kinases IPK1 and ITPK1 constitute a metabolic pathway in maintaining phosphate homeostasis

Emerging studies have suggested that there is a close link between inositol phosphate (InsP) metabolism and cellular phosphate (Pi ) homeostasis in eukaryotes; however, whether a common InsP species is deployed as an evolutionarily conserved metabolic messenger to mediate Pi signaling remains unknown. Here, using genetics and InsP profiling combined with Pi -starvation response (PSR) analysis in Arabidopsis thaliana, we showed that the kinase activity of inositol pentakisphosphate 2-kinase (IPK1), an enzyme required for phytate (inositol hexakisphosphate; InsP6 ) synthesis, is indispensable for maintaining Pi homeostasis under Pi -replete conditions, and inositol 1,3,4-trisphosphate 5/6-kinase 1 (ITPK1) plays an equivalent role. Although both ipk1-1 and itpk1 mutants exhibited decreased levels of InsP6 and diphosphoinositol pentakisphosphate (PP-InsP5 ; InsP7 ), disruption of another ITPK family enzyme, ITPK4, which correspondingly caused depletion of InsP6 and InsP7 , did not display similar Pi -related phenotypes, which precludes these InsP species from being effectors. Notably, the level of d/l-Ins(3,4,5,6)P4 was concurrently elevated in both ipk1-1 and itpk1 mutants, which showed a specific correlation with the misregulated Pi phenotypes. However, the level of d/l-Ins(3,4,5,6)P4 is not responsive to Pi starvation that instead manifests a shoot-specific increase in the InsP7 level. This study demonstrates a more nuanced picture of the intersection of InsP metabolism and Pi homeostasis and PSRs than has previously been elaborated, and additionally establishes intermediate steps to phytate biosynthesis in plant vegetative tissues.

Identification of Morus notabilis MADS-box genes and elucidation of the roles of MnMADS33 during endodormancy

The MADS-box genes encode transcriptional regulators with various functions especially during floral development. A total of 54 putative Morus notabilis MADS-box genes (MnMADSs) were identified and phylogenetically classified as either type I (17 genes) or type II (37 genes). The detected genes included three FLOWERING LOCUS C-like (MnFLC-like) genes, MnMADS33, MnMADS50, and MnMADS7. MnFLC-like proteins could directly or indirectly repress promoter activity of the mulberry FLOWERING LOCUS T-like (MnFT) gene. Transgenic Arabidopsis thaliana overexpressing MnFLC-like genes exhibited delayed flowering and down-regulated expression of FT and SUPPRESSOR OF OVEREXPRESSION OF CONSTANS1 (SOC1). The gene expression analyses in floral bud indicated that MnMADS33 expression increased, while MnFT expression decreased during the induction of dormancy in response to cold conditions. Dormancy release resulted in the down-regulation of MnMADS33 expression and the up-regulation of MnFT expression. Furthermore, abscisic acid promoted the transcription of MnMADS33 and MnFT, although the expression level of MnFT gradually decreased. These results are consistent with the hypothesis that MnMADS33 negatively regulated the expression of MnFT to repress dormancy release and flowering in mulberry. This study may be relevant for future investigations regarding the effects of MnMADS genes on mulberry flowering development.

Epigenetic regulation of agronomical traits in Brassicaceae

Epigenetic regulation, covalent modification of DNA and changes in histone proteins are closely linked to plant development and stress response through flexibly altering the chromatin structure to regulate gene expression. In this review, we will illustrate the importance of epigenetic influences by discussing three agriculturally important traits of Brassicaceae. (1) Vernalization, an acceleration of flowering by prolonged cold exposure regulated through epigenetic silencing of a central floral repressor, FLOWERING LOCUS C. This is associated with cold-dependent repressive histone mark accumulation, which confers competency of consequence vegetative-to-reproductive phase transition. (2) Hybrid vigor, in which an F₁ hybrid shows superior performance to the parental lines. Combination of distinct epigenomes with different DNA methylation states between parental lines is important for increase in growth rate in a hybrid progeny. This is independent of siRNA-directed DNA methylation but dependent on the chromatin remodeler DDM1. (3) Self-incompatibility, a reproductive mating system to prevent self-fertilization. This is controlled by the S-locus consisting of SP11 and SRK which are responsible for self/non-self recognition. Because self-incompatibility in Brassicaceae is sporophytically controlled, there are dominance relationships between S haplotypes in the stigma and pollen. The dominance relationships in the pollen rely on de novo DNA methylation at the promoter region of a recessive allele, which is triggered by siRNA production from a flanking region of a dominant allele.

Comparative analysis of molecular and physiological traits between perennial Arabis alpina Pajares and annual Arabidopsis thaliana Sy-0

Annual plants complete life cycle in a year while perennial plants maintain growth for several years. Arabis alpina, a polycarpic perennial, is a close relative of monocarpic annual Arabidopsis. Pajares is an accession of A. alpina requiring vernalization, a long-term cold for flowering. Arabidopsis shows holistic flowering whereas Pajares shows idiographic flowering, producing axillary branches under variable developmental phases from juvenile, adult vegetative to reproductive phases. To understand the molecular mechanism behind diverse phases of axillary branches, we analyzed the levels of primary miR156 expressions because miR156-SPL module is a key regulator for developmental phase transition. We found that in Pajares, miR156 levels were highly variable among the axillary branches, which causes differential sensitivity to vernalization. Thus, the axillary branches expressing high levels of miR156 remain in juvenile phase even after vernalization, whereas the axillary branches expressing low levels of miR156 produce flowers after vernalization. In contrast, every axillary branches of Arabidopsis winter annual Sy-0 expressed similar levels of miR156 and synchronously responded to vernalization, which causes holistic flowering. Therefore, we suggest that variable miR156 expression levels and the resulting differential response to vernalization among axillary branches are distinctive features determining polycarpic perenniality of A. alpina Pajares.

Flowering Time Gene Variation in Brassica Species Shows Evolutionary Principles

Flowering time genes have a strong influence on successful reproduction and life cycle adaptation. However, their regulation is highly complex and only well understood in diploid model systems. For crops with a polyploid background from the genus Brassica, data on flowering time gene variation are scarce, although indispensable for modern breeding techniques like marker-assisted breeding. We have deep-sequenced all paralogs of 35 Arabidopsis thaliana flowering regulators using Sequence Capture followed by Illumina sequencing in two selected accessions of the vegetable species Brassica rapa and Brassica oleracea, respectively. Using these data, we were able to call SNPs, InDels and copy number variations (CNVs) for genes from the total flowering time network including central flowering regulators, but also genes from the vernalisation pathway, the photoperiod pathway, temperature regulation, the circadian clock and the downstream effectors. Comparing the results to a complementary data set from the allotetraploid species Brassica napus, we detected rearrangements in B. napus which probably occurred early after the allopolyploidisation event. Those data are both a valuable resource for flowering time research in those vegetable species, as well as a contribution to speciation genetics.

Nucleotide polymorphism affecting FLC expression underpins heading date variation in horticultural brassicas

Variation in flowering time and response to overwintering has been exploited to breed brassica vegetables that can be harvested year-round. Our knowledge of flowering time control now enables the investigation of the molecular basis of this important variation. Here, we show that a major determinant of heading date variation in Brassica oleracea is from variation in vernalization response through allelic variation at FLOWERING LOCUS C.C2 (BoFLC4). We characterize two alleles of BoFLC.C2 that are both functional and confer a requirement for vernalization, but they show distinct expression dynamics in response to cold. Complementation experiments in Arabidopsis thaliana revealed that the allelic variation results from cis polymorphism at BoFLC.C2, which quantitatively influences the degree of cold-induced epigenetic silencing. This results in one allelic variant conferring consistently later heading under both glasshouse and field conditions through reduced environmental sensitivity. Our results suggest that breeding of brassica varieties for commercially valuable variation in heading date has been achieved through the selection of cis polymorphism at FLC, similar to that underpinning natural variation in A. thaliana. This understanding will allow for the selection of alleles with distinct sensitivities to cold and robust heading dates under variable climatic conditions, and will facilitate the breeding of varieties more resistant to climate change.

Quantitative trait loci controlling leaf appearance and curd initiation of cauliflower in relation to temperature

QTL regions on chromosomes C06 and C09 are involved in temperature dependent time to curd induction in cauliflower. Temperature is the main environmental factor influencing curding time of cauliflower (Brassica oleracea var. botrytis). Temperatures above 20-22 °C inhibit development towards curding even in many summer cultivars. To identify quantitative trait loci (QTL) controlling curding time and its related traits in a wide range of different temperature regimes from 12 to 27 °C, a doubled haploid (DH) mapping population segregating for curding time was developed and days to curd initiation (DCI), leaf appearance rate (LAR), and final leaf number (FLN) were measured. The population was genotyped with 176 single nucleotide polymorphism (SNP) markers. Composite interval mapping (CIM) revealed repeatedly detected QTL for DCI on C06 and C09. The estimated additive effect increased at high temperatures. Significant QTL × environment interactions (Q × E) for FLN and DCI on C06 and C09 suggest that these hotspot regions have major influences on temperature mediated curd induction. 25 % of the DH lines did not induce curds at temperatures higher than 22 °C. Applying a binary model revealed a QTL with LOD >15 on C06. Nearly all lines carrying the allele of the reliable early maturing parental line (PL) on that locus induced curds at high temperatures while only half of the DH lines carrying the allele of the unreliable PL reached the generative phase during the experiment. Large variation in LAR was observed. QTL for LAR were detected repeatedly in several environments on C01, C04 and C06. Negative correlations between LAR and DCI and QTL co-localizations on C04 and C06 suggest that LAR has also effects on development towards curd induction.

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.

Lineage isolation in the face of active gene flow in the coastal plant wild radish is reinforced by differentiated vernalisation responses

Background

The respective role and relative importance of natural selection and gene flow in the process of population divergence has been a central theme in the speciation literature. A previous study presented conclusive evidence that wild radish on Japanese islands comprises two genetically isolated lineages: the southern and northern groups. However, a general understanding of the lineage isolation with frequent seed flow of the coastal plant species is still unclear. We surveyed nucleotide polymorphisms over 14 nuclear loci in 72 individuals across the Japan–Ryukyu Islands Arc to address the demographic history of wild radish utilising the isolation-with-migration (IM) model. In addition, we investigated the flowering times of individuals in different wild radish lineages, with and without cold exposure, to assess their respective vernalisation responses.

Results

Coalescent simulations suggested that divergence between the southern and northern lineages of wild radish began ~18,000 years ago, initially during the Last Glacial Maximum (LGM) period. The gene flow from the southern to northern groups was considerably higher than that in the opposite direction, indicating effective dispersal of viable seeds via the northward Kuroshio Current. Our greenhouse experiments indicated that cold exposure was not required for flowering in the southern group, but could advance the date of flowering, suggesting that vernalisation would be facultative in the southern group. In contrast, the northern group was either unable to flower or flowered later without prior cold exposure, and thus had an obligate requirement for cold treatment.

Conclusions

The south–north lineage divergence in wild radish could be triggered by a directional change in the sea current during the ice age, despite gene flow due to the high dispersability and longevity of seeds. We also found that temperature profoundly affected the vernalisation responses of wild radish, which may repress reproductive success and ultimately drive and reinforce intra-specific differentiation between the two lineages of wild radish. This study provides new insights into the maintenance of lineage differentiation with on-going gene flow in coastal plants.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0655-7) contains supplementary material, which is available to authorized users.

Identification of plant vacuolar transporters mediating phosphate storage

Plant vacuoles serve as the primary intracellular compartments for inorganic phosphate (Pi) storage. Passage of Pi across vacuolar membranes plays a critical role in buffering the cytoplasmic Pi level against fluctuations of external Pi and metabolic activities. Here we demonstrate that the SPX-MFS proteins, designated as PHOSPHATE TRANSPORTER 5 family (PHT5), also named Vacuolar Phosphate Transporter (VPT), function as vacuolar Pi transporters. Based on ³¹P-magnetic resonance spectroscopy analysis, Arabidopsis pht5;1 loss-of-function mutants accumulate less Pi and exhibit a lower vacuolar-to-cytoplasmic Pi ratio than controls. Conversely, overexpression of PHT5 leads to massive Pi sequestration into vacuoles and altered regulation of Pi starvation-responsive genes. Furthermore, we show that heterologous expression of the rice homologue OsSPX-MFS1 mediates Pi influx to yeast vacuoles. Our findings show that a group of Pi transporters in vacuolar membranes regulate cytoplasmic Pi homeostasis and are required for fitness and plant growth.

The plant vacuole acts as a storage compartment for inorganic phosphate and buffers cytoplasmic phosphate concentration. Here, Liu et al. identify a group of vacuolar phosphate transporters in Arabidopsis that are required for plant growth in response to fluctuating availability of phosphate.

Tuning growth cycles of Brassica crops via natural antisense transcripts of BrFLC

Several oilseed and vegetable crops of Brassica are biennials that require a prolonged winter cold for flowering, a process called vernalization. FLOWERING LOCUS C (FLC) is a central repressor of flowering. Here, we report that the overexpression of natural antisense transcripts (NATs) of Brassica rapa FLC (BrFLC) greatly shortens plant growth cycles. In rapid-, medium- and slow-cycling crop types, there are four copies of the BrFLC genes, which show extensive variation in sequences and expression levels. In Bre, a biennial crop type that requires vernalization, five NATs derived from the BrFLC2 locus are rapidly induced under cold conditions, while all four BrFLC genes are gradually down-regulated. The transgenic Bre lines overexpressing a long NAT of BrFLC2 do not require vernalization, resulting in a gradient of shortened growth cycles. Among them, a subset of lines both flower and set seeds as early as Yellow sarson, an annual crop type in which all four BrFLC genes have non-sense mutations and are nonfunctional in flowering repression. Our results demonstrate that the growth cycles of biennial crops of Brassica can be altered by changing the expression levels of BrFLC2 NATs. Thus, BrFLC2 NATs and their transgenic lines are useful for the genetic manipulation of crop growth cycles.

Variation in the flowering time orthologs BrFLC and BrSOC1 in a natural population of Brassica rapa

Understanding the genetic basis of natural phenotypic variation is of great importance, particularly since selection can act on this variation to cause evolution. We examined expression and allelic variation in candidate flowering time loci in Brassica rapa plants derived from a natural population and showing a broad range in the timing of first flowering. The loci of interest were orthologs of the Arabidopsis genes FLC and SOC1 (BrFLC and BrSOC1, respectively), which in Arabidopsis play a central role in the flowering time regulatory network, with FLC repressing and SOC1 promoting flowering. In B. rapa, there are four copies of FLC and three of SOC1. Plants were grown in controlled conditions in the lab. Comparisons were made between plants that flowered the earliest and latest, with the difference in average flowering time between these groups ∼30 days. As expected, we found that total expression of BrSOC1 paralogs was significantly greater in early than in late flowering plants. Paralog-specific primers showed that expression was greater in early flowering plants in the BrSOC1 paralogs Br004928, Br00393 and Br009324, although the difference was not significant in Br009324. Thus expression of at least 2 of the 3 BrSOC1 orthologs is consistent with their predicted role in flowering time in this natural population. Sequences of the promoter regions of the BrSOC1 orthologs were variable, but there was no association between allelic variation at these loci and flowering time variation. For the BrFLC orthologs, expression varied over time, but did not differ between the early and late flowering plants. The coding regions, promoter regions and introns of these genes were generally invariant. Thus the BrFLC orthologs do not appear to influence flowering time in this population. Overall, the results suggest that even for a trait like flowering time that is controlled by a very well described genetic regulatory network, understanding the underlying genetic basis of natural variation in such a quantitative trait is challenging.

The role of BoFLC2 in cauliflower (Brassica oleracea var. botrytis L.) reproductive development

In agricultural species that are sexually propagated or whose marketable organ is a reproductive structure, management of the flowering process is critical. Inflorescence development in cauliflower is particularly complex, presenting unique challenges for those seeking to predict and manage flowering time. In this study, an integrated physiological and molecular approach was used to clarify the environmental control of cauliflower reproductive development at the molecular level. A functional allele of BoFLC2 was identified for the first time in an annual brassica, along with an allele disrupted by a frameshift mutation (boflc2). In a segregating F₂ population derived from a cross between late-flowering (BoFLC2) and early-flowering (boflc2) lines, this gene behaved in a dosage-dependent manner and accounted for up to 65% of flowering time variation. Transcription of BoFLC genes was reduced by vernalization, with the floral integrator BoFT responding inversely. Overall expression of BoFT was significantly higher in early-flowering boflc2 lines, supporting the idea that BoFLC2 plays a key role in maintaining the vegetative state. A homologue of Arabidopsis VIN3 was isolated for the first time in a brassica crop species and was up-regulated by two days of vernalization, in contrast to findings in Arabidopsis where prolonged exposure to cold was required to elicit up-regulation. The correlations observed between gene expression and flowering time in controlled-environment experiments were validated with gene expression analyses of cauliflowers grown outdoors under 'natural' vernalizing conditions, indicating potential for transcript levels of flowering genes to form the basis of predictive assays for curd initiation and flowering time.

Gene regulatory variation mediates flowering responses to vernalization along an altitudinal gradient in Arabidopsis

Steep environmental gradients provide ideal settings for studies of potentially adaptive phenotypic and genetic variation in plants. The accurate timing of flowering is crucial for reproductive success and is regulated by several pathways, including the vernalization pathway. Among the numerous genes known to enable flowering in response to vernalization, the most prominent is FLOWERING LOCUS C (FLC). FLC and other genes of the vernalization pathway vary extensively among natural populations and are thus candidates for the adaptation of flowering time to environmental gradients such as altitude. We used 15 natural Arabidopsis (Arabidopsis thaliana) genotypes originating from an altitudinal gradient (800-2,700 m above sea level) in the Swiss Alps to test whether flowering time correlated with altitude under different vernalization scenarios. Additionally, we measured the expression of 12 genes of the vernalization pathway and its downstream targets. Flowering time correlated with altitude in a nonlinear manner for vernalized plants. Flowering time could be explained by the expression and regulation of the vernalization pathway, most notably by AGAMOUS LIKE19 (AGL19), FLOWERING LOCUS T (FT), and FLC. The expression of AGL19, FT, and VERNALIZATION INSENSITIVE3 was associated with altitude, and the regulation of MADS AFFECTING FLOWERING2 (MAF2) and MAF3 differed between low- and high-altitude genotypes. In conclusion, we found clinal variation across an altitudinal gradient both in flowering time and the expression and regulation of genes in the flowering time control network, often independent of FLC, suggesting that the timing of flowering may contribute to altitudinal adaptation.

Arabidopsis inositol pentakisphosphate 2-kinase, AtIPK1, is required for growth and modulates phosphate homeostasis at the transcriptional level

Inositol hexakisphosphate (IP6 ) provides a phosphorous reservoir in plant seeds; in addition, along with its biosynthesis intermediates and derivatives, IP6 also plays important roles in diverse developmental and physiological processes. Disruption of the Arabidopsis inositol pentakisphosphate 2-kinase coding gene AtIPK1 was previously shown to reduce IP6 content in vegetative tissues and affect phosphate (Pi) sensing. Here we show that AtIPK1 is required for sustaining plant growth, as null mutants are non-viable. An incomplete loss-of-function mutant, atipk1-1, exhibited disturbed Pi homeostasis and overaccumulated Pi as a consequence of increased Pi uptake activity and root-to-shoot Pi translocation. The atipk1-1 mutants also showed a Pi deficiency-like root system architecture with reduced primary root and enhanced lateral root growth. Transcriptome analysis indicated that a subset of Pi starvation-responsive genes was transcriptionally perturbed in the atipk1-1 mutants and the expression of multiple genes involved in Pi uptake, allocation, and remobilization was increased. Genetic and transcriptional analyses suggest that disturbance of Pi homeostasis caused by atipk1 mutation involved components in addition to PHR1(-like) transcription factors. Notably, the transcriptional increase of a number of Pi starvation-responsive genes in the atipk1-1 mutants is correlated with the reduction of histone variant H2A.Z occupation in chromatin. The myo-inositol-1-phosphate synthase mutants, atmips1 and atmips2 with comparable reduction in vegetative IP6 to that in the atipk1-1 mutants did not overaccumulate Pi, suggesting that Pi homeostasis modulated by AtIPK1 is not solely attributable to IP6 level. This study reveals that AtIPK1 has important roles in growth and Pi homeostasis.

Evolutionary conservation of cold-induced antisense RNAs of FLOWERING LOCUS C in Arabidopsis thaliana perennial relatives

Antisense RNA (asRNA) COOLAIR is expressed at A. thaliana FLOWERING LOCUS C (FLC) in response to winter temperatures. Its contribution to cold-induced silencing of FLC was proposed but its functional and evolutionary significance remain unclear. Here we identify a highly conserved block containing the COOLAIR first exon and core promoter at the 3′ end of several FLC orthologues. Furthermore, asRNAs related to COOLAIR are expressed at FLC loci in the perennials A. alpina and A. lyrata, although some splicing variants differ from A. thaliana. Study of the A. alpina orthologue, PERPETUAL FLOWERING 1 (PEP1), demonstrates that AaCOOLAIR is induced each winter of the perennial life cycle. Introduction of PEP1 into A. thaliana reveals that AaCOOLAIR cis-elements confer cold-inducibility in this heterologous species while the difference between PEP1 and FLC mRNA patterns depends on both cis-elements and species-specific trans-acting factors. Thus, expression of COOLAIR is highly conserved, supporting its importance in FLC regulation.

FLOWERING LOCUS C (FLC) is thought to control the flowering time of A. thaliana in response to winter temperatures, in a process known as vernalization. Here, the authors suggest that the COOLAIR antisense RNA, which is conserved across plant species, acts to repress the expression of FLC during vernalization.

An allotetraploid Brassica napus early-flowering mutant has BnaFLC2-regulated flowering

BACKGROUND

Flowering time is an important agronomic trait, and wide variation in flowering time exists among Brassica napus accessions. GX50 early-flowering mutant, induced from Brassica napus by Ethyl Methane Sulfonate (EMS), exhibits a remarkable early transition from vegetative to reproductive growth.

RESULTS

GX50 plants flowered about 60 days earlier than the control wild-type plant B. napus XY15 under greenhouse conditions. Cytological examination revealed that the GX50 plants form inflorescences as early as from 5 weeks old, flower primordium from 6 weeks old, and siliques from 10 weeks old, whereas 10-week-old XY15 plants are still at vegetative growth stage. To unravel the molecular mechanisms underlying the GX50 flowering phenotype, we analyzed the expression of several key regulatory genes. Expressions of all five BnaFLCs (BnaFLC1 to BnaFLC5), BnaFT and BnaSOC1 were detected. Interestingly, BnaFLCs expression levels were lower in GX50 than those in XY15. Among the five BnaFLCs, only the expression pattern of BnaFLC2 corresponded to the timing of floral organ differentiation in GX50. In agreement with previous knowledge that BnaFLCs repress expression of BnaFT and BnaSOC1, increased levels of BnaFT and BnaSOC1 were observed in GX50 compared with XY15.

CONCLUSION

BnaFLC2, but not the other BnaFLC genes, plays an important role in B. napus GX50 floral transition.

Nitrogen limitation adaptation, a target of microRNA827, mediates degradation of plasma membrane-localized phosphate transporters to maintain phosphate homeostasis in Arabidopsis

Members of the Arabidopsis thaliana phosphate transporter1 (PHT1) family are key players in acquisition of Pi from the rhizosphere, and their regulation is indispensable for the maintenance of cellular Pi homeostasis. Here, we reveal posttranslational regulation of Pi transport through modulation of degradation of PHT1 proteins by the RING-type ubiquitin E3 ligase, nitrogen limitation adaptation (NLA). Loss of function of NLA caused high Pi accumulation resulting from increases in the levels of several PHT1s at the protein rather than the transcript level. Evidence of decreased endocytosis and ubiquitination of PHT1s in nla mutants and interaction between NLA and PHT1s in the plasma membranes suggests that NLA directs the ubiquitination of plasma membrane-localized PHT1s, which triggers clathrin-dependent endocytosis followed by endosomal sorting to vacuoles. Furthermore, different subcellular localization of NLA and phosphate2 (pho2; a ubiquitin E2 conjugase) and the synergistic effect of the accumulation of PHT1s and Pi in nla pho2 mutants suggest that they function independently but cooperatively to regulate PHT1 protein amounts. Intriguingly, NLA and PHO2 are the targets of two Pi starvation-induced microRNAs, miR827 and miR399, respectively. Therefore, our findings uncover modulation of Pi transport activity in response to Pi availability through the integration of a microRNA-mediated posttranscriptional pathway and a ubiquitin-mediated posttranslational regulatory pathway.

A root chicory MADS box sequence and the Arabidopsis flowering repressor FLC share common features that suggest conserved function in vernalization and de-vernalization responses

Root chicory (Cichorium intybus var. sativum) is a biennial crop, but is harvested to obtain root inulin at the end of the first growing season before flowering. However, cold temperatures may vernalize seeds or plantlets, leading to incidental early flowering, and hence understanding the molecular basis of vernalization is important. A MADS box sequence was isolated by RT-PCR and named FLC-LIKE1 (CiFL1) because of its phylogenetic positioning within the same clade as the floral repressor Arabidopsis FLOWERING LOCUS C (AtFLC). Moreover, over-expression of CiFL1 in Arabidopsis caused late flowering and prevented up-regulation of the AtFLC target FLOWERING LOCUS T by photoperiod, suggesting functional conservation between root chicory and Arabidopsis. Like AtFLC in Arabidopsis, CiFL1 was repressed during vernalization of seeds or plantlets of chicory, but repression of CiFL1 was unstable when the post-vernalization temperature was favorable to flowering and when it de-vernalized the plants. This instability of CiFL1 repression may be linked to the bienniality of root chicory compared with the annual lifecycle of Arabidopsis. However, re-activation of AtFLC was also observed in Arabidopsis when a high temperature treatment was used straight after seed vernalization, eliminating the promotive effect of cold on flowering. Cold-induced down-regulation of a MADS box floral repressor and its re-activation by high temperature thus appear to be conserved features of the vernalization and de-vernalization responses in distant species.

Adaptation to seasonality and the winter freeze

Flowering plants initially diversified during the Mesozoic era at least 140 million years ago in regions of the world where temperate seasonal environments were not encountered. Since then several cooling events resulted in the contraction of warm and wet environments and the establishment of novel temperate zones in both hemispheres. In response, less than half of modern angiosperm families have members that evolved specific adaptations to cold seasonal climates, including cold acclimation, freezing tolerance, endodormancy, and vernalization responsiveness. Despite compelling evidence for multiple independent origins, the level of genetic constraint on the evolution of adaptations to seasonal cold is not well understood. However, the recent increase in molecular genetic studies examining the response of model and crop species to seasonal cold offers new insight into the evolutionary lability of these traits. This insight has major implications for our understanding of complex trait evolution, and the potential role of local adaptation in response to past and future climate change. In this review, we discuss the biochemical, morphological, and developmental basis of adaptations to seasonal cold, and synthesize recent literature on the genetic basis of these traits in a phylogenomic context. We find evidence for multiple genetic links between distinct physiological responses to cold, possibly reinforcing the coordinated expression of these traits. Furthermore, repeated recruitment of the same or similar ancestral pathways suggests that land plants might be somewhat pre-adapted to dealing with temperature stress, perhaps making inducible cold traits relatively easy to evolve.

Genetic and physical mapping of flowering time loci in canola (Brassica napus L.)

We identified quantitative trait loci (QTL) underlying variation for flowering time in a doubled haploid (DH) population of vernalisation-responsive canola (Brassica napus L.) cultivars Skipton and Ag-Spectrum and aligned them with physical map positions of predicted flowering genes from the Brassica rapa genome. Significant genetic variation in flowering time and response to vernalisation were observed among the DH lines from Skipton/Ag-Spectrum. A molecular linkage map was generated comprising 674 simple sequence repeat, sequence-related amplified polymorphism, sequence characterised amplified region, Diversity Array Technology, and candidate gene based markers loci. QTL analysis indicated that flowering time is a complex trait and is controlled by at least 20 loci, localised on ten different chromosomes. These loci each accounted for between 2.4 and 28.6% of the total genotypic variation for first flowering and response to vernalisation. However, identification of consistent QTL was found to be dependant upon growing environments. We compared the locations of QTL with the physical positions of predicted flowering time genes located on the sequenced genome of B. rapa. Some QTL associated with flowering time on A02, A03, A07, and C06 may represent homologues of known flowering time genes in Arabidopsis; VERNALISATION INSENSITIVE 3, APETALA1, CAULIFLOWER, FLOWERING LOCUS C, FLOWERING LOCUS T, CURLY LEAF, SHORT VEGETATIVE PHASE, GA3 OXIDASE, and LEAFY. Identification of the chromosomal location and effect of the genes influencing flowering time may hasten the development of canola varieties having an optimal time for flowering in target environments such as for low rainfall areas, via marker-assisted selection.

Genetic and physical mapping of flowering time loci in canola (Brassica napus L.)

We identified quantitative trait loci (QTL) underlying variation for flowering time in a doubled haploid (DH) population of vernalisation-responsive canola (Brassica napus L.) cultivars Skipton and Ag-Spectrum and aligned them with physical map positions of predicted flowering genes from the Brassica rapa genome. Significant genetic variation in flowering time and response to vernalisation were observed among the DH lines from Skipton/Ag-Spectrum. A molecular linkage map was generated comprising 674 simple sequence repeat, sequence-related amplified polymorphism, sequence characterised amplified region, Diversity Array Technology, and candidate gene based markers loci. QTL analysis indicated that flowering time is a complex trait and is controlled by at least 20 loci, localised on ten different chromosomes. These loci each accounted for between 2.4 and 28.6% of the total genotypic variation for first flowering and response to vernalisation. However, identification of consistent QTL was found to be dependant upon growing environments. We compared the locations of QTL with the physical positions of predicted flowering time genes located on the sequenced genome of B. rapa. Some QTL associated with flowering time on A02, A03, A07, and C06 may represent homologues of known flowering time genes in Arabidopsis; VERNALISATION INSENSITIVE 3, APETALA1, CAULIFLOWER, FLOWERING LOCUS C, FLOWERING LOCUS T, CURLY LEAF, SHORT VEGETATIVE PHASE, GA3 OXIDASE, and LEAFY. Identification of the chromosomal location and effect of the genes influencing flowering time may hasten the development of canola varieties having an optimal time for flowering in target environments such as for low rainfall areas, via marker-assisted selection.

Role of vernalization and of duplicated FLOWERING LOCUS C in the perennial Arabidopsis lyrata

FLOWERING LOCUS C (FLC) is one of the main genes influencing the vernalization requirement and natural flowering time variation in the annual Arabidopsis thaliana. Here we studied the effects of vernalization on flowering and its genetic basis in the perennial Arabidopsis lyrata. Two tandemly duplicated FLC genes (FLC1 and FLC2) were compared with respect to expression and DNA sequence. The effect of vernalization on flowering and on the expression of FLC1 was studied in three European populations. The genetic basis of the FLC1 expression difference between two of the populations was further studied by expression quantitative trait locus (eQTL) mapping and sequence analysis. FLC1 was shown to have a likely role in the vernalization requirement for flowering in A. lyrata. Vernalization decreased its expression and the northern study populations showed higher FLC1 expression than the southern one. eQTL mapping between two of the populations revealed one eQTL affecting FLC1 expression in the genomic region containing the FLC genes. Most FLC1 sequence differences between the study populations were found in the promoter region and in the first intron. Variation in the FLC1 sequence may cause differences in FLC1 expression between late- and early-flowering A. lyrata populations.