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

Non-vernalization requirement for flowering in Brassica rapa conferred by a dominant allele of FLOWERING LOCUS T

KEY MESSAGE

We identified and characterized a dominant FT allele for flowering without vernalization in Brassica rapa, while demonstrating its potential for deployment in breeding to accelerate flowering in various Brassicaceae crops. Controlling the timing of flowering is key to improving yield and quality of several agricultural crops including the Brassicas. Many Brassicaceae crops possess a conserved flowering mechanism in which FLOWERING LOCUS C (FLC) represses the transcription of flowering activators such as FLOWERING LOCUS T (FT) during vernalization. Here, we employed genetic analysis based on next-generation sequencing to identify a dominant FT allele, BraA.FT.2-C, for flowering in the absence of vernalization in the Brassica rapa cultivar 'CHOY SUM EX CHINA 3'. BraA.FT.2-C harbors two large insertions upstream of its coding region and is expressed without vernalization, despite FLC expression. We show that BraA.FT.2-C offers an opportunity to introduce flowering without vernalization requirement into winter-type brassica crops, including B. napus, which have many functional FLC paralogs. Furthermore, we demonstrated the feasibility of using B. rapa harboring BraA.FT.2-C as rootstock for grafting to induce flowering in radish (Raphanus sativus), which requires vernalization for flowering. We believe that the ability of BraA.FT.2-C to overcome repression by FLC can have significant applications in brassica crops breeding to increase yields by accelerating or delaying flowering.

Rapid, nonparallel genomic evolution of Brassica rapa (field mustard) under experimental drought

While we know that climate change can potentially cause rapid phenotypic evolution, our understanding of the genetic basis and degree of genetic parallelism of rapid evolutionary responses to climate change is limited. In this study, we combined the resurrection approach with an evolve-and-resequence design to examine genome-wide evolutionary changes following drought. We exposed genetically similar replicate populations of the annual plant Brassica rapa derived from a field population in southern California to four generations of experimental drought or watered conditions in a greenhouse. Genome-wide sequencing of ancestral and descendant population pools identified hundreds of SNPs that showed evidence of rapidly evolving in response to drought. Several of these were in stress response genes, and two were identified in a prior study of drought response in this species. However, almost all genetic changes were unique among experimental populations, indicating that the evolutionary changes were largely nonparallel, despite the fact that genetically similar replicates of the same founder population had experienced controlled and consistent selection regimes. This nonparallelism of evolution at the genetic level is potentially because of polygenetic adaptation allowing for multiple different genetic routes to similar phenotypic outcomes. Our findings help to elucidate the relationship between rapid phenotypic and genomic evolution and shed light on the degree of parallelism and predictability of genomic evolution to environmental change.

A functionally divergent SOC1 homolog improves soybean yield and latitudinal adaptation

Soybean (Glycine max) grows in a wide range of latitudes, but it is extremely sensitive to photoperiod, which reduces its yield and ability to adapt to different environments. Therefore, understanding of the genetic basis of soybean adaptation is of great significance for breeding and improvement. Here, we characterized Tof18 (SOC1a) that conditions early flowering and growth habit under both short-day and long-day conditions. Molecular analysis confirmed that the two SOC1 homologs present in soybeans (SOC1a and SOC1b) underwent evolutionary functional divergence, with SOC1a having stronger effects on flowering time and stem node number than SOC1b due to transcriptional differences. soc1a soc1b double mutants showed stronger functional effects than either of the single mutants, perhaps due to the formation of SOC1a and SOC1b homodimers or heterodimers. Additionally, Tof18/SOC1a improves the latitudinal adaptation of cultivated soybeans, highlighting the functional importance of SOC1a. The Tof18G allele facilitates adaptation to high latitudes, whereas Tof18A facilitates adaptation to low latitudes. We demonstrated that SOC1s contribute to floral induction in both leaves and shoot apex through inter-regulation with FTs. The SOC1a-SOC1b-Dt2 complex plays essential roles in stem growth habit by directly binding to the regulatory sequence of Dt1, making the genes encoding these proteins potential targets for genome editing to improve soybean yield via molecular breeding. Since the natural Tof18A allele increases node number, introgressing this allele into modern cultivars could improve yields, which would help optimize land use for food production in the face of population growth and global warming.

A single amino acid residue substitution in BraA04g017190.3C, a histone methyltransferase, results in premature bolting in Chinese cabbage (Brassica rapa L. ssp. Pekinensis)

BACKGROUND

Flowering is an important inflection point in the transformation from vegetative to reproductive growth, and premature bolting severely decreases crop yield and quality.

RESULTS

In this study, a stable early-bolting mutant, ebm3, was identified in an ethyl methanesulfonate (EMS)-mutagenized population of a Chinese cabbage doubled haploid (DH) line 'FT'. Compared with 'FT', ebm3 showed early bolting under natural cultivation in autumn, and curled leaves. Genetic analysis showed that the early-bolting phenotype was controlled by a single recessive nuclear gene. Modified MutMap sequencing, genotyping analyses and allelism test provide strong evidence that BrEBM3 (BraA04g017190.3 C), encoding the histone methyltransferase CURLY LEAF (CLF), was the strongly candidate gene of the emb3. A C to T base substitution in the 14th exon of BrEBM3 resulted in an amino acid change (S to F) and the early-bolting phenotype of emb3. The mutation occurred in the SET domain (Suppressor of protein-effect variegation 3-9, Enhancer-of-zeste, Trithorax), which catalyzes site- and state-specific lysine methylation in histones. Tissue-specific expression analysis showed that BrEBM3 was highly expressed in the flower and bud. Promoter activity assay confirmed that BrEBM3 promoter was active in inflorescences. Subcellular localization analysis revealed that BrEBM3 localized in the nucleus. Transcriptomic studies supported that BrEBM3 mutation might repress H3K27me3 deposition and activate expression of the AGAMOUS (AG) and AGAMOUS-like (AGL) loci, resulting in early flowering.

CONCLUSIONS

Our study revealed that an EMS-induced early-bolting mutant ebm3 in Chinese cabbage was caused by a nonsynonymous mutation in BraA04g017190.3 C, encoding the histone methyltransferase CLF. These results improve our knowledge of the genetic and genomic resources of bolting and flowering, and may be beneficial to the genetic improvement of Chinese cabbage.

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.

Rapid evolutionary changes in gene expression in response to climate fluctuations

There is now abundant evidence of rapid evolution in natural populations, but the genetic mechanisms of these changes remain unclear. One possible route to rapid evolution is through changes in the expression of genes that influence traits under selection. We examined contemporary evolutionary gene expression changes in plant populations responding to environmental fluctuations. We compared genome‐wide gene expression, using RNA‐seq, in two populations of Brassica rapa collected over four time points between 1997 and 2014, during which precipitation in southern California fluctuated dramatically and phenotypic and genotypic changes occurred. By combining transcriptome profiling with the resurrection approach, we directly examined evolutionary changes in gene expression over time. For both populations, we found a substantial number of differentially expressed genes between generations, indicating rapid evolution in the expression of many genes. Using existing gene annotations, we found that many changes occurred in genes involved in regulating stress responses and flowering time. These appeared related to the fluctuations in precipitation and were potentially adaptive. However, the evolutionary changes in gene expression differed across generations within and between populations, indicating largely independent evolutionary trajectories across populations and over time. Our study provides strong evidence for rapid evolution in gene expression, and indicates that changes in gene expression can be one mechanism of rapid evolutionary responses to selection episodes. This study also illustrates that combining resurrection studies with transcriptomics is a powerful approach for investigating evolutionary changes at the gene regulatory level, and will provide new insights into the genetic basis of contemporary evolution.

see also the Perspective by Emily B. Josephs.

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.

BrFLC5: a weak regulator of flowering time in Brassica rapa

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

Highly preserved roles of Brassica MIR172 in polyploid Brassicas: ectopic expression of variants of Brassica MIR172 accelerates floral transition

Functional characterization of regulatory genes governing flowering time is a research priority for breeding earliness in crop Brassicas. Highly polyploid genomes of Brassicas pose challenges in unraveling homeolog gene function. In Arabidopsis, five MIR172 paralogs control flowering time and floral organ identity by down-regulating AP2 and AP2-like genes. The impact of homeolog diversification on MIR172 loci, however, needs to be examined in morphologically diverse Brassicas. Herein, we analyze fractionation status and phylogeny of MIR172 and target AP2 from Brassicas and compare functionality of MIR172 variants representing distinct sub-genomes and progenitor genomes. Copy number analysis revealed higher retention of MIR172 loci relative to AP2 in diploid and amphi-diploid Brassica species. Dendrogram of 87 MIR172 sequences from Brassicaceae showed five major clusters corresponding to MIR172a-MIR172e which further separated into sub-genome and progenitor genome specific clades. Similar groupings were observed in the phylogeny of 11 Brassica AP2 and AP2-like genes. Over-expression of a pair of natural variants for each of MIR172b, MIR172d and MIR172e representing sub-genomes, progenitor genomes and species of Brassicas displayed floral acceleration in all transgenic lines indicating a strong selection pressure on MIR172. All gain-of-function lines, except 35S::MIR172e and 35S::MIR172e' displayed floral organ defects implying altered target spectrum of MIR172e relative to MIR172b and MIR172d. Expression of MIR172e caused marginal earliness in flowering time in B. juncea. In conclusion, this study demonstrates tightly preserved role of homeologs and natural variants of MIR172 family in mediating flowering in Brassicas and suggests their deployment for introgression of early flowering trait.

FLC expression is down-regulated by cold treatment in Diplotaxis tenuifolia (wild rocket), but flowering time is unaffected

Wild rocket (Diplotaxis tenuifolia) has become a very popular salad leaf due to its peppery taste. It is part of the Brassicaceae family and thus has a high level of homology at the DNA level to other Brassica species including Arabidopsis thaliana. The vernalization and photoperiodic requirements of wild rocket have not been reported to date. Photoperiodic experiments described here demonstrate that rocket is a facultative long day plant. To investigate the vernalization requirement, both seed and young plants were given vernalization treatments at 4°C for different lengths of time. A rocket homologue of FLOWERING LOCUS C (DtFLC) was isolated and shown to functionally complement the Arabidopsis FRI⁺flc3 null mutant. Whilst the expression of DtFLC was significantly reduced after just one week of cold treatment, cold treatments of two to eight weeks had no significant effect on bolting time of wild rocket indicating that rocket does not have a vernalization requirement. These findings illustrate that important fundamental differences can exist between model and crop plant species, such as in this case where down-regulation of DtFLC expression does not enable earlier flowering in wild rocket as it does in Arabidopsis and many other Brassica species.

Genome-Wide Characterization of the MADS-Box Gene Family in Radish (Raphanus sativus L.) and Assessment of Its Roles in Flowering and Floral Organogenesis

The MADS-box gene family is an important transcription factor (TF) family that is involved in various aspects of plant growth and development, especially flowering time and floral organogenesis. Although it has been reported in many plant species, the systematic identification and characterization of MADS-box TF family is still limited in radish (Raphanus sativus L.). In the present study, a comprehensive analysis of MADS-box genes was performed, and a total of 144 MADS-box family members were identified from the whole radish genome. Meanwhile, a detailed list of MADS-box genes from other 28 plant species was also investigated. Through the phylogenetic analysis between radish and Arabidopsis thaliana, all the RsMADS genes were classified into two groups including 68 type I (31 Mα, 12 Mβ and 25Mγ) and 76 type II (70 MIKC^C and 6 MIKC^∗). Among them, 41 (28.47%) RsMADS genes were located in nine linkage groups of radish from R1 to R9. Moreover, the homologous MADS-box gene pairs were identified among radish, A. thaliana, Chinese cabbage and rice. Additionally, the expression profiles of RsMADS genes were systematically investigated in different tissues and growth stages. Furthermore, quantitative real-time PCR analysis was employed to validate expression patterns of some crucial RsMADS genes. These results could provide a valuable resource to explore the potential functions of RsMADS genes in radish, and facilitate dissecting MADS-box gene-mediated molecular mechanisms underlying flowering and floral organogenesis in root vegetable crops.