Role of FRIGIDA and FLOWERING LOCUS C in determining variation in flowering time of Arabidopsis

Allopolyploidy expanded gene content but not pangenomic variation in the hexaploid oilseed Camelina sativa

Ancient whole-genome duplications are believed to facilitate novelty and adaptation by providing the raw fuel for new genes. However, it is unclear how recent whole-genome duplications may contribute to evolvability within recent polyploids. Hybridization accompanying some whole-genome duplications may combine divergent gene content among diploid species. Some theory and evidence suggest that polyploids have a greater accumulation and tolerance of gene presence-absence and genomic structural variation, but it is unclear to what extent either is true. To test how recent polyploidy may influence pangenomic variation, we sequenced, assembled, and annotated 12 complete, chromosome-scale genomes of Camelina sativa, an allohexaploid biofuel crop with 3 distinct subgenomes. Using pangenomic comparative analyses, we characterized gene presence-absence and genomic structural variation both within and between the subgenomes. We found over 75% of ortholog gene clusters are core in C. sativa and <10% of sequence space was affected by genomic structural rearrangements. In contrast, 19% of gene clusters were unique to one subgenome, and the majority of these were Camelina specific (no ortholog in Arabidopsis). We identified an inversion that may contribute to vernalization requirements in winter-type Camelina and an enrichment of Camelina-specific genes with enzymatic processes related to seed oil quality and Camelina's unique glucosinolate profile. Genes related to these traits exhibited little presence-absence variation. Our results reveal minimal pangenomic variation in this species and instead show how hybridization accompanied by whole-genome duplication may benefit polyploids by merging diverged gene content of different species.

A systematic review on the implications of concurrent heat and drought stress in modulating floral development in plants

The continuous change in climate, along with irregular rainfall patterns, poses a significant threat to sustainable agricultural productivity worldwide. Both high temperatures and drought stress are key factors limiting crop growth, and with global climate change, the occurrence of combined heat and drought stress is expected to rise. This will further exacerbate the vulnerability of agricultural yield. Simultaneous heat and drought stress is prevalent in field conditions, and while extensive research has been done on the individual effects of heat and drought stress on plants, little is known about the molecular mechanisms underlying plant acclimation to a combination of these stressors. The reproductive stage, especially the flowering phase, has been identified as the most sensitive to both heat and drought stress, leading to sterility in plants. However, our understanding of the combined stress response in commonly used crop plants is still limited. Hence, it is crucial to study and comprehend the effects and interactions between high temperatures and drought stress during the reproductive stages of crops. This review delves into the morpho-physiological changes in reproductive organs of various plant species under combined heat and drought stress and also details the molecular regulation of the mechanism of combined stress tolerance in plants. Notably, the article incorporates expression analyses of candidate genes in rice flowers, emphasizing the utilization of modern biotechnological methods to enhance stress tolerance in plants. Overall, the review provides a comprehensive insight into the regulation of floral development in plants following concurrent heat and drought stress.

C-TERMINAL DOMAIN PHOSPHATASE-LIKE 1 promotes flowering with TAF15b by repressing the floral repressor gene FLOWERING LOCUS C

Arabidopsis TATA-BINDING PROTEIN-ASSOCIATED FACTOR15b (TAF15b) is a plant-specific component of the transcription factor IID complex. TAF15b is involved in the autonomous pathway for flowering and represses the transcription of FLOWERING LOCUS C (FLC), a major floral repressor in Arabidopsis. While components of the autonomous flowering pathway have been extensively studied, scant attention has been directed toward elucidating the direct transcriptional regulators responsible for repressing FLC transcription. Here, we demonstrate that C-TERMINAL DOMAIN PHOSPHATASE-LIKE 1 (CPL1) is a physical and functional partner of TAF15b, playing a role in FLC repression. CPL1 is a protein phosphatase that dephosphorylates the C-terminal domain of RNA polymerase II (Pol II). Through the immunoprecipitation and mass spectrometry technique, we identified CPL1 as an interacting partner of TAF15b. Similar to taf15b, the cpl1 mutant showed a late-flowering phenotype caused by an increase in FLC levels. Additionally, the increase in cpl1 was correlated with the enrichment of phosphorylated Pol II in the FLC chromatin, as expected. We also discovered that CPL1 and TAF15b share additional common target genes through transcriptome analysis. These results suggest that TAF15b and CPL1 cooperatively repress transcription through the dephosphorylation of Pol II, especially at the FLC locus.

Ecological trade-offs drive phenotypic and genetic differentiation of Arabidopsis thaliana in Europe

Plant diversity is shaped by trade-offs between traits related to competitive ability, propagule dispersal, and stress resistance. However, we still lack a clear understanding of how these trade-offs influence species distribution and population dynamics. In Arabidopsis thaliana, recent genetic analyses revealed a group of cosmopolitan genotypes that successfully recolonized Europe from its center after the last glaciation, excluding older (relict) lineages from the distribution except for their north and south margins. Here, we tested the hypothesis that cosmopolitans expanded due to higher colonization ability, while relicts persisted at the margins due to higher tolerance to competition and/or stress. We compared the phenotypic and genetic differentiation between 71 European genotypes originating from the center, and the south and north margins. We showed that a trade-off between plant fecundity and seed mass shapes the differentiation of A. thaliana in Europe, suggesting that the success of the cosmopolitan groups could be explained by their high dispersal ability. However, at both north and south margins, we found evidence of selection for alleles conferring low dispersal but highly competitive and stress-resistance abilities. This study sheds light on the role of ecological trade-offs as evolutionary drivers of the distribution and dynamics of plant populations.

Genome-wide discovery of InDels and validation of PCR-Based InDel markers for earliness in a RIL population and genotypes of lentil (Lens culinaris Medik.)

The systematic identification of insertion/deletion (InDel) length polymorphisms from the entire lentil genome can be used to map the quantitative trait loci (QTL) and also for the marker-assisted selection (MAS) for various linked traits. The InDels were identified by comparing the whole-genome resequencing (WGRS) data of two extreme bulks (early- and late-flowering bulk) and a parental genotype (Globe Mutant) of lentil. The bulks were made by pooling 20 extreme recombinant inbred lines (RILs) each, derived by crossing Globe Mutant (late flowering parent) with L4775 (early flowering parent). Finally, 734,716 novel InDels were identified, which is nearly one InDel per 5,096 bp of lentil genome. Furthermore, 74.94% of InDels were within the intergenic region and 99.45% displayed modifier effects. Of these, 15,732 had insertions or deletions of 20 bp or more, making them amenable to the development of PCR-based markers. An InDel marker I-SP-356.6 (chr. 3; position 356,687,623; positioned 174.5 Kb from the LcFRI gene) was identified as having a phenotypic variance explained (PVE) value of 47.7% for earliness when validated in a RIL population. Thus, I-SP-356.6 marker can be deployed in MAS to facilitate the transfer of the earliness trait to other elite late-maturing cultivars. Two InDel markers viz., I-SP-356.6 and I-SP-383.9 (chr. 3; linked to LcELF3a gene) when tested in 9 lentil genotypes differing for maturity duration, clearly distinguished three early (L4775, ILL7663, Precoz) and four late genotypes (Globe Mutant, MFX, L4602, L830). However, these InDels could not be validated in two genotypes (L4717, L4727), suggesting either absence of polymorphism and/or presence of other loci causing earliness. The identified InDel markers can act as valuable tools for MAS for the development of early maturing lentil varieties.

Flowering time: From physiology, through genetics to mechanism

Plant species have evolved different requirements for environmental/endogenous cues to induce flowering. Originally, these varying requirements were thought to reflect the action of different molecular mechanisms. Thinking changed when genetic and molecular analysis in Arabidopsis thaliana revealed that a network of environmental and endogenous signaling input pathways converge to regulate a common set of "floral pathway integrators." Variation in the predominance of the different input pathways within a network can generate the diversity of requirements observed in different species. Many genes identified by flowering time mutants were found to encode general developmental and gene regulators, with their targets having a specific flowering function. Studies of natural variation in flowering were more successful at identifying genes acting as nodes in the network central to adaptation and domestication. Attention has now turned to mechanistic dissection of flowering time gene function and how that has changed during adaptation. This will inform breeding strategies for climate-proof crops and help define which genes act as critical flowering nodes in many other species.

Genome-wide association study identifies a novel BMI1A QTL allele that confers FLC expression diversity in Arabidopsis thaliana

Identification and understanding of the genetic basis of natural variations in plants are essential for comprehending their phenotypic adaptation. Here, we report a genome-wide association study (GWAS) of FLOWERING LOCUS C (FLC) expression in 727 Arabidopsis accessions. We identified B LYMPHOMA MOLONEY MURINE LEUKEMIA VIRUS INSERTION REGION 1 HOMOLOG 1A (BMI1A) as a causal gene for one of the FLC expression quantitative trait loci (QTLs). Loss of function in BMI1A increases FLC expression and delays flowering time at 16 °C significantly compared with the wild type (Col-0). BMI1A activity is required for histone H3 lysine 27 trimethylation (H3K27me3) accumulation at the FLC, MADS AFFECTING FLOWERING 4 (MAF4), and MAF5 loci at low ambient temperature. We further uncovered two BMI1A haplotypes associated with the natural variation in FLC expression and flowering time at 16 °C, and demonstrated that polymorphisms in the BMI1A promoter region are the main contributor. Different BMI1A haplotypes are strongly associated with geographical distribution, and the low ambient temperature-sensitive BMI1A variants are associated with a lower mean temperature of the driest quarter of their collection sites compared with the temperature-non-responsive variants, indicating that the natural variations in BMI1A have adaptive functions in FLC expression and flowering time regulation. Therefore, our results provide new insights into the natural variations in FLC expression and flowering time diversity in plants.

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.

An overview of floral regulatory genes in annual and perennial plants

The floral initiation in angiosperms is a complex process influenced by endogenous and exogenous signals. With this approach, we aim to provide a comprehensive review to integrate this complex floral regulatory process and summarize the regulatory genes and their functions in annuals and perennials. Seven primary paths leading to flowering have been discovered in Arabidopsis under several growth condition that include; photoperiod, ambient temperature, vernalization, gibberellins, autonomous, aging and carbohydrates. These pathways involve a series of interlinked signaling pathways that respond to both internal and external signals, such as light, temperature, hormones, and developmental cues, to coordinate the expression of genes that are involved in flower development. Among them, the photoperiodic pathway was the most important and conserved as some of the fundamental loci and mechanisms are shared even by closely related plant species. The activation of floral regulatory genes such as FLC, FT, LFY, and SOC1 that determine floral meristem identity and the transition to the flowering stage result from the merging of these pathways. Recent studies confirmed that alternative splicing, antisense RNA and epigenetic modification play crucial roles by regulating the expression of genes related to blooming. In this review, we documented recent progress in the floral transition time in annuals and perennials, with emphasis on the specific regulatory mechanisms along with the application of various molecular approaches including overexpression studies, RNA interference and Virus-induced flowering. Furthermore, the similarities and differences between annual and perennial flowering will aid significant contributions to the field by elucidating the mechanisms of perennial plant development and floral initiation regulation.

Responsiveness to long days for flowering is reduced in Arabidopsis by yearly variation in growing season temperatures

Conservative flowering behaviours, such as flowering during long days in summer or late flowering at a high leaf number, are often proposed to protect against variable winter and spring temperatures which lead to frost damage if premature flowering occurs. Yet, due the many factors in natural environments relative to the number of individuals compared, assessing which climate characteristics drive these flowering traits has been difficult. We applied a multidisciplinary approach to 10 winter-annual Arabidopsis thaliana populations from a wide climactic gradient in Norway. We used a variable reduction strategy to assess which of 100 climate descriptors from their home sites correlated most to their flowering behaviours when tested for responsiveness to photoperiod after saturation of vernalization; then, assessed sequence variation of 19 known environmental-response flowering genes. Photoperiod responsiveness inversely correlated with interannual variation in timing of growing season onset. Time to flowering appeared driven by growing season length, curtailed by cold fall temperatures. The distribution of FLM, TFL2 and HOS1 haplotypes, genes involved in ambient temperature response, correlated with growing-season climate. We show that long-day responsiveness and late flowering may be driven not by risk of spring frosts, but by growing season temperature and length, perhaps to opportunistically maximize growth.

Causal role of a promoter polymorphism in natural variation of the Arabidopsis floral repressor gene FLC

Noncoding polymorphism frequently associates with phenotypic variation, but causation and mechanism are rarely established. Noncoding single-nucleotide polymorphisms (SNPs) characterize the major haplotypes of the Arabidopsis thaliana floral repressor gene FLOWERING LOCUS C (FLC). This noncoding polymorphism generates a range of FLC expression levels, determining the requirement for and the response to winter cold. The major adaptive determinant of these FLC haplotypes was shown to be the autumnal levels of FLC expression. Here, we investigate how noncoding SNPs influence FLC transcriptional output. We identify an upstream transcription start site (uTSS) cluster at FLC, whose usage is increased by an A variant at the promoter SNP-230. This variant is present in relatively few Arabidopsis accessions, with the majority containing G at this site. We demonstrate a causal role for the A variant at -230 in reduced FLC transcriptional output. The G variant upregulates FLC expression redundantly with the major transcriptional activator FRIGIDA (FRI). We demonstrate an additive interaction of SNP-230 with an intronic SNP+259, which also differentially influences uTSS usage. Combinatorial interactions between noncoding SNPs and transcriptional activators thus generate quantitative variation in FLC transcription that has facilitated the adaptation of Arabidopsis accessions to distinct climates.

Plant environmental memory: implications, mechanisms and opportunities for plant scientists and beyond

Plants are extremely plastic organisms. They continuously receive and integrate environmental information and adjust their growth and development to favour fitness and survival. When this integration of information affects subsequent life stages or the development of subsequent generations, it can be considered an environmental memory. Thus, plant memory is a relevant mechanism by which plants respond adaptively to different environments. If the cost of maintaining the response is offset by its benefits, it may influence evolutionary trajectories. As such, plant memory has a sophisticated underlying molecular mechanism with multiple components and layers. Nonetheless, when mathematical modelling is combined with knowledge of ecological, physiological, and developmental effects as well as molecular mechanisms as a tool for understanding plant memory, the combined potential becomes unfathomable for the management of plant communities in natural and agricultural ecosystems. In this review, we summarize recent advances in the understanding of plant memory, discuss the ecological requirements for its evolution, outline the multilayered molecular network and mechanisms required for accurate and fail-proof plant responses to variable environments, point out the direct involvement of the plant metabolism and discuss the tremendous potential of various types of models to further our understanding of the plant's environmental memory. Throughout, we emphasize the use of plant memory as a tool to unlock the secrets of the natural world.

Pan-European study of genotypes and phenotypes in the Arabidopsis relative Cardamine hirsuta reveals how adaptation, demography, and development shape diversity patterns

We study natural DNA polymorphisms and associated phenotypes in the Arabidopsis relative Cardamine hirsuta. We observed strong genetic differentiation among several ancestry groups and broader distribution of Iberian relict strains in European C. hirsuta compared to Arabidopsis. We found synchronization between vegetative and reproductive development and a pervasive role for heterochronic pathways in shaping C. hirsuta natural variation. A single, fast-cycling ChFRIGIDA allele evolved adaptively allowing range expansion from glacial refugia, unlike Arabidopsis where multiple FRIGIDA haplotypes were involved. The Azores islands, where Arabidopsis is scarce, are a hotspot for C. hirsuta diversity. We identified a quantitative trait locus (QTL) in the heterochronic SPL9 transcription factor as a determinant of an Azorean morphotype. This QTL shows evidence for positive selection, and its distribution mirrors a climate gradient that broadly shaped the Azorean flora. Overall, we establish a framework to explore how the interplay of adaptation, demography, and development shaped diversity patterns of 2 related plant species.

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

Molecular mechanisms of adaptive evolution in wild animals and plants

Wild animals and plants have developed a variety of adaptive traits driven by adaptive evolution, an important strategy for species survival and persistence. Uncovering the molecular mechanisms of adaptive evolution is the key to understanding species diversification, phenotypic convergence, and inter-species interaction. As the genome sequences of more and more non-model organisms are becoming available, the focus of studies on molecular mechanisms of adaptive evolution has shifted from the candidate gene method to genetic mapping based on genome-wide scanning. In this study, we reviewed the latest research advances in wild animals and plants, focusing on adaptive traits, convergent evolution, and coevolution. Firstly, we focused on the adaptive evolution of morphological, behavioral, and physiological traits. Secondly, we reviewed the phenotypic convergences of life history traits and responding to environmental pressures, and the underlying molecular convergence mechanisms. Thirdly, we summarized the advances of coevolution, including the four main types: mutualism, parasitism, predation and competition. Overall, these latest advances greatly increase our understanding of the underlying molecular mechanisms for diverse adaptive traits and species interaction, demonstrating that the development of evolutionary biology has been greatly accelerated by multi-omics technologies. Finally, we highlighted the emerging trends and future prospects around the above three aspects of adaptive evolution.

AGO1 and HSP90 buffer different genetic variants in Arabidopsis thaliana

Argonaute 1 (AGO1), the principal protein component of microRNA-mediated regulation, plays a key role in plant growth and development. AGO1 physically interacts with the chaperone HSP90, which buffers cryptic genetic variation in plants and animals. We sought to determine whether genetic perturbation of AGO1 in Arabidopsis thaliana would also reveal cryptic genetic variation, and if so, whether AGO1-dependent loci overlap with those dependent on HSP90. To address these questions, we introgressed a hypomorphic mutant allele of AGO1 into a set of mapping lines derived from the commonly used Arabidopsis strains Col-0 and Ler. Although we identified several cases in which AGO1 buffered genetic variation, none of the AGO1-dependent loci overlapped with those buffered by HSP90 for the traits assayed. We focused on 1 buffered locus where AGO1 perturbation uncoupled the traits days to flowering and rosette leaf number, which are otherwise closely correlated. Using a bulk segregant approach, we identified a nonfunctional Ler hua2 mutant allele as the causal AGO1-buffered polymorphism. Introduction of a nonfunctional hua2 allele into a Col-0 ago1 mutant background recapitulated the Ler-dependent ago1 phenotype, implying that coupling of these traits involves different molecular players in these closely related strains. Taken together, our findings demonstrate that even though AGO1 and HSP90 buffer genetic variation in the same traits, these robustness regulators interact epistatically with different genetic loci, suggesting that higher-order epistasis is uncommon. Plain Language Summary Argonaute 1 (AGO1), a key player in plant development, interacts with the chaperone HSP90, which buffers environmental and genetic variation. We found that AGO1 buffers environmental and genetic variation in the same traits; however, AGO1-dependent and HSP90-dependent loci do not overlap. Detailed analysis of a buffered locus found that a nonfunctional HUA2 allele decouples days to flowering and rosette leaf number in an AGO1-dependent manner, suggesting that the AGO1-dependent buffering acts at the network level.

Vernalization-triggered expression of the antisense transcript COOLAIR is mediated by CBF genes

To synchronize flowering time with spring, many plants undergo vernalization, a floral-promotion process triggered by exposure to long-term winter cold. In Arabidopsis thaliana, this is achieved through cold-mediated epigenetic silencing of the floral repressor, FLOWERING LOCUS C (FLC). COOLAIR, a cold-induced antisense RNA transcribed from the FLC locus, has been proposed to facilitate FLC silencing. Here, we show that C-repeat (CRT)/dehydration-responsive elements (DREs) at the 3'-end of FLC and CRT/DRE-binding factors (CBFs) are required for cold-mediated expression of COOLAIR. CBFs bind to CRT/DREs at the 3'-end of FLC, both in vitro and in vivo, and CBF levels increase gradually during vernalization. Cold-induced COOLAIR expression is severely impaired in cbfs mutants in which all CBF genes are knocked-out. Conversely, CBF-overexpressing plants show increased COOLAIR levels even at warm temperatures. We show that COOLAIR is induced by CBFs during early stages of vernalization but COOLAIR levels decrease in later phases as FLC chromatin transitions to an inactive state to which CBFs can no longer bind. We also demonstrate that cbfs and FLC_ΔCOOLAIR mutants exhibit a normal vernalization response despite their inability to activate COOLAIR expression during cold, revealing that COOLAIR is not required for the vernalization process.

Epistatic interactions among multiple copies of FLC genes with naturally occurring insertions correlate with flowering time variation in radish

Brassicaceae crops, which underwent whole-genome triplication during their evolution, have multiple copies of flowering-related genes. Interactions among multiple gene copies may be involved in flowering time regulation; however, this mechanism is poorly understood. In this study, we performed comprehensive, high-throughput RNA sequencing analysis to identify candidate genes involved in the extremely late-bolting (LB) trait in radish. Then, we examined the regulatory roles and interactions of radish FLOWERING LOCUS C (RsFLC) paralogs, the main flowering repressor candidates. Seven flowering integrator genes, five vernalization genes, nine photoperiodic/circadian clock genes and eight genes from other flowering pathways were differentially expressed in the early-bolting (EB) cultivar 'Aokubinagafuto' and LB radish cultivar 'Tokinashi' under different vernalization conditions. In the LB cultivar, RsFLC1 and RsFLC2 expression levels were maintained after 40 days of cold exposure. Bolting time was significantly correlated with the expression rates of RsFLC1 and RsFLC2. Using the EB × LB F2 population, we performed association analyses of genotypes with or without 1910- and 1627-bp insertions in the first introns of RsFLC1 and RsFLC2, respectively. The insertion alleles prevented the repression of their respective FLC genes under cold conditions. Interestingly, genotypes homozygous for RsFLC2 insertion alleles maintained high RsFLC1 and RsFLC3 expression levels under cold conditions, and two-way analysis of variance revealed that RsFLC1 and RsFLC3 expression was influenced by the RsFLC2 genotype. Our results indicate that insertions in the first introns of RsFLC1 and RsFLC2 contribute to the late-flowering trait in radish via different mechanisms. The RsFLC2 insertion allele conferred a strong delay in bolting by inhibiting the repression of all three RsFLC genes, suggesting that radish flowering time is determined by epistatic interactions among multiple FLC gene copies.

Transcriptome analysis of critical genes related to flowering in Mikania micrantha at different altitudes provides insights for a potential control

BACKGROUND

Mikania micrantha is a vine with strong invasion ability, and its strong sexual reproduction ability is not only the main factor of harm, but also a serious obstacle to control. M. micrantha spreads mainly through seed production. Therefore, inhibiting the flowering and seed production of M. micrantha is an effective strategy to prevent from continuing to spread.

RESULT

The flowering number of M. micrantha is different at different altitudes. A total of 67.01 Gb of clean data were obtained from nine cDNA libraries, and more than 83.47% of the clean reads were mapped to the reference genome. In total, 5878 and 7686 significantly differentially expressed genes (DEGs) were found in E2 vs. E9 and E13 vs. E9, respectively. Based on the background annotation and gene expression, some candidate genes related to the flowering pathway were initially screened, and their expression levels in the three different altitudes in flower bud differentiation showed the same trend. That is, at an altitude of 1300 m, the flower integration gene and flower meristem gene were downregulated (such as SOC1 and AP1), and the flowering inhibition gene was upregulated (such as FRI and SVP). Additionally, the results showed that there were many DEGs involved in the hormone signal transduction pathway in the flower bud differentiation of M. micrantha at different altitudes.

CONCLUSIONS

Our results provide abundant sequence resources for clarifying the underlying mechanisms of flower bud differentiation and mining the key factors inhibiting the flowering and seed production of M. micrantha to provide technical support for the discovery of an efficient control method.

Burning questions for a warming and changing world: 15 unknowns in plant abiotic stress

We present unresolved questions in plant abiotic stress biology as posed by 15 research groups with expertise spanning eco-physiology to cell and molecular biology. Common themes of these questions include the need to better understand how plants detect water availability, temperature, salinity, and rising carbon dioxide (CO2) levels; how environmental signals interface with endogenous signaling and development (e.g. circadian clock and flowering time); and how this integrated signaling controls downstream responses (e.g. stomatal regulation, proline metabolism, and growth versus defense balance). The plasma membrane comes up frequently as a site of key signaling and transport events (e.g. mechanosensing and lipid-derived signaling, aquaporins). Adaptation to water extremes and rising CO2 affects hydraulic architecture and transpiration, as well as root and shoot growth and morphology, in ways not fully understood. Environmental adaptation involves tradeoffs that limit ecological distribution and crop resilience in the face of changing and increasingly unpredictable environments. Exploration of plant diversity within and among species can help us know which of these tradeoffs represent fundamental limits and which ones can be circumvented by bringing new trait combinations together. Better defining what constitutes beneficial stress resistance in different contexts and making connections between genes and phenotypes, and between laboratory and field observations, are overarching challenges.

Molecular epigenetic mechanisms for the memory of temperature stresses in plants

The sessile plants encounter various stresses; some are prolonged, whereas some others are recurrent. Temperature is crucial for plant growth and development, and plants often encounter adverse high temperature fluctuations (heat stresses) as well as prolonged cold exposure such as seasonal temperature drops in winter when grown in temperate regions. Many plants can remember past temperature stresses to get adapted to adverse local temperature changes to ensure survival and/or reproductive success. Here, we summarize chromatin-based mechanisms underlying acquired thermotolerance or thermomemory in plants and review recent progresses on molecular epigenetic understanding of 'remembering of prolonged cold in winter' or vernalization, a process critical for various over-wintering plants to acquire competence to flower in the coming spring. In addition, perspectives on future study in temperature stress memories of economically-important crops are discussed.

Genome-Wide Identification of Maize Protein Arginine Methyltransferase Genes and Functional Analysis of ZmPRMT1 Reveal Essential Roles in Arabidopsis Flowering Regulation and Abiotic Stress Tolerance

Histone methylation, as one of the important epigenetic regulatory mechanisms, plays a significant role in growth and developmental processes and stress responses of plants, via altering the methylation status or ratio of arginine and lysine residues of histone tails, which can affect the regulation of gene expression. Protein arginine methyltransferases (PRMTs) have been revealed to be responsible for histone methylation of specific arginine residues in plants, which is important for maintaining pleiotropic development and adaptation to abiotic stresses in plants. Here, for the first time, a total of eight PRMT genes in maize have been identified and characterized in this study, named as ZmPRMT1-8. According to comparative analyses of phylogenetic relationship and structural characteristics among PRMT gene family members from several representative species, all maize 8 PRMT proteins were categorized into three distinct subfamilies. Further, schematic structure and chromosome location analyses displayed evolutionarily conserved structure features and an unevenly distribution on maize chromosomes of ZmPRMT genes, respectively. The expression patterns of ZmPRMT genes in different tissues and under various abiotic stresses (heat, drought, and salt) were determined. The expression patterns of ZmPRMT genes indicated that they play a role in regulating growth and development and responses to abiotic stress. Eventually, to verify the biological roles of ZmPRMT genes, the transgenic Arabidopsis plants overexpressing ZmPRMT1 gene was constructed as a typical representative. The results demonstrated that overexpression of ZmPRMT1 can promote earlier flowering time and confer enhanced heat tolerance in transgenic Arabidopsis. Taken together, our results are the first to report the roles of ZmPRMT1 gene in regulating flowering time and resisting heat stress response in plants and will provide a vital theoretical basis for further unraveling the functional roles and epigenetic regulatory mechanism of ZmPRMT genes in maize growth, development and responses to abiotic stresses.

Mechanisms of Vernalization-Induced Flowering in Legumes

Vernalization is the requirement for exposure to low temperatures to trigger flowering. The best knowledge about the mechanisms of vernalization response has been accumulated for Arabidopsis and cereals. In Arabidopsis thaliana, vernalization involves an epigenetic silencing of the MADS-box gene FLOWERING LOCUS C (FLC), which is a flowering repressor. FLC silencing releases the expression of the main flowering inductor FLOWERING LOCUS T (FT), resulting in a floral transition. Remarkably, no FLC homologues have been identified in the vernalization-responsive legumes, and the mechanisms of cold-mediated transition to flowering in these species remain elusive. Nevertheless, legume FT genes have been shown to retain the function of the main vernalization signal integrators. Unlike Arabidopsis, legumes have three subclades of FT genes, which demonstrate distinct patterns of regulation with respect to environmental cues and tissue specificity. This implies complex mechanisms of vernalization signal propagation in the flowering network, that remain largely elusive. Here, for the first time, we summarize the available information on the genetic basis of cold-induced flowering in legumes with a special focus on the role of FT genes.

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.

2b-RAD Genotyping of the Seagrass Cymodocea nodosa Along a Latitudinal Cline Identifies Candidate Genes for Environmental Adaptation

Plant populations distributed along broad latitudinal gradients often show patterns of clinal variation in genotype and phenotype. Differences in photoperiod and temperature cues across latitudes influence major phenological events, such as timing of flowering or seed dormancy. Here, we used an array of 4,941 SNPs derived from 2b-RAD genotyping to characterize population differentiation and levels of genetic and genotypic diversity of three populations of the seagrass Cymodocea nodosa along a latitudinal gradient extending across the Atlantic-Mediterranean boundary (i.e., Gran Canaria—Canary Islands, Faro—Portugal, and Ebro Delta—Spain). Our main goal was to search for potential outlier loci that could underlie adaptive differentiation of populations across the latitudinal distribution of the species. We hypothesized that such polymorphisms could be related to variation in photoperiod-temperature regime occurring across latitudes. The three populations were clearly differentiated and exhibited diverse levels of clonality and genetic diversity. Cymodocea nodosa from the Mediterranean displayed the highest genotypic richness, while the Portuguese population had the highest clonality values. Gran Canaria exhibited the lowest genetic diversity (as observed heterozygosity). Nine SNPs were reliably identified as outliers across the three sites by two different methods (i.e., BayeScan and pcadapt), and three SNPs could be associated to specific protein-coding genes by screening available C. nodosa transcriptomes. Two SNPs-carrying contigs encoded for transcription factors, while the other one encoded for an enzyme specifically involved in the regulation of flowering time, namely Lysine-specific histone demethylase 1 homolog 2. When analyzing biological processes enriched within the whole dataset of outlier SNPs identified by at least one method, “regulation of transcription” and “signalling” were among the most represented. Our results highlight the fundamental importance signal integration and gene-regulatory networks, as well as epigenetic regulation via DNA (de)methylation, could have for enabling adaptation of seagrass populations along environmental gradients.

The genetic basis of natural variation in the timing of vegetative phase change in Arabidopsis thaliana

The juvenile-to-adult transition in plants is known as vegetative phase change and is marked by changes in the expression of leaf traits in response to a decrease in the level of miR156 and miR157. To determine whether this is the only mechanism of vegetative phase change, we measured the appearance of phase-specific leaf traits in 70 natural accessions of Arabidopsis thaliana. We found that leaf shape was poorly correlated with abaxial trichome production (two adult traits), that variation in these traits was not necessarily correlated with the level of miR156, and that there was little to no correlation between the appearance of adult-specific vegetative traits and flowering time. We identified eight quantitative trait loci controlling phase-specific vegetative traits from a cross between the Columbia (Col-0) and Shakdara (Sha) accessions. Only one of these quantitative trait loci includes genes known to regulate vegetative phase change (MIR156A and TOE1), which were expressed at levels consistent with the precocious phenotype of Sha. Our results suggest that vegetative phase change is regulated both by the miR156/SPL module and by genes specific to different vegetative traits, and that natural variation in vegetative phase change can arise from either source.

Identification and Characterization of the MIKC-Type MADS-Box Gene Family in Brassica napus and Its Role in Floral Transition

Increasing rapeseed yield has always been a primary goal of rapeseed research and breeding. However, flowering time is a prerequisite for stable rapeseed yield and determines its adaptability to ecological regions. MIKC-type MADS-box (MICK) genes are a class of transcription factors that are involved in various physiological and developmental processes in plants. To understand their role in floral transition-related pathways, a genome-wide screening was conducted with Brassica napus (B. napus), which revealed 172 members. Using previous data from a genome-wide association analysis of flowering traits, BnaSVP and BnaSEP1 were identified as candidate flowering genes. Therefore, we used the CRISPR/Cas9 system to verify the function of BnaSVP and BnaSEP1 in B. napus. T0 plants were edited efficiently at the BnaSVP and BnaSEP1 target sites to generate homozygous and heterozygous mutants with most mutations stably inherited by the next generation. Notably, the mutant only showed the early flowering phenotype when all homologous copies of BnaSVP were edited, indicating functional redundancy between homologous copies. However, no changes in flowering were observed in the BnaSEP1 mutant. Quantitative analysis of the pathway-related genes in the BnaSVP mutant revealed the upregulation of SUPPRESSOR OF OVEREXPRESSION OF CONSTANS 1 (SOC1) and FLOWERING LOCUS T (FT) genes, which promoted early flowering in the mutant. In summary, our study created early flowering mutants, which provided valuable resources for early maturing breeding, and provided a new method for improving polyploid crops.

Natural variation and improved genome annotation of the emerging biofuel crop field pennycress (Thlaspi arvense)

The Brassicaceae family comprises more than 3,700 species with a diversity of phenotypic characteristics, including seed oil content and composition. Recently, the global interest in Thlaspi arvense L. (pennycress) has grown as the seed oil composition makes it a suitable source for biodiesel and aviation fuel production. However, many wild traits of this species need to be domesticated to make pennycress ideal for cultivation. Molecular breeding and engineering efforts require the availability of an accurate genome sequence of the species. Here, we describe pennycress genome annotation improvements, using a combination of long- and short-read transcriptome data obtained from RNA derived from embryos of 22 accessions, in addition to public genome and gene expression information. Our analysis identified 27,213 protein-coding genes, as well as on average 6,188 biallelic SNPs. In addition, we used the identified SNPs to evaluate the population structure of our accessions. The data from this analysis support that the accession Ames 32872, originally from Armenia, is highly divergent from the other accessions, while the accessions originating from Canada and the United States cluster together. When we evaluated the likely signatures of natural selection from alternative SNPs, we found 7 candidate genes under likely recent positive selection. These genes are enriched with functions related to amino acid metabolism and lipid biosynthesis and highlight possible future targets for crop improvement efforts in pennycress.

Adaptive significance of flowering time variation across natural seasonal environments in Arabidopsis thaliana

The relevance of flowering time variation and plasticity to climate adaptation requires a comprehensive empirical assessment. We investigated natural selection and the genetic architecture of flowering time in Arabidopsis through field experiments in Europe across multiple sites and seasons. We estimated selection for flowering time, plasticity and canalization. Loci associated with flowering time, plasticity and canalization by genome-wide association studies were tested for a geographic signature of climate adaptation. Selection favored early flowering and increased canalization, except at the northernmost site, but was rarely detected for plasticity. Genome-wide association studies revealed significant associations with flowering traits and supported a substantial polygenic inheritance. Alleles associated with late flowering, including functional FRIGIDA variants, were more common in regions experiencing high annual temperature variation. Flowering time plasticity to fall vs spring and summer environments was associated with GIGANTEA SUPPRESSOR 5, which promotes early flowering under decreasing day length and temperature. The finding that late flowering genotypes and alleles are associated with climate is evidence for past adaptation. Real-time phenotypic selection analysis, however, reveals pervasive contemporary selection for rapid flowering in agricultural settings across most of the species range. The response to this selection may involve genetic shifts in environmental cuing compared to the ancestral state.

Parallel reduction in flowering time from de novo mutations enable evolutionary rescue in colonizing lineages

Understanding how populations adapt to abrupt environmental change is necessary to predict responses to future challenges, but identifying specific adaptive variants, quantifying their responses to selection and reconstructing their detailed histories is challenging in natural populations. Here, we use Arabidopsis from the Cape Verde Islands as a model to investigate the mechanisms of adaptation after a sudden shift to a more arid climate. We find genome-wide evidence of adaptation after a multivariate change in selection pressures. In particular, time to flowering is reduced in parallel across islands, substantially increasing fitness. This change is mediated by convergent de novo loss of function of two core flowering time genes: FRI on one island and FLC on the other. Evolutionary reconstructions reveal a case where expansion of the new populations coincided with the emergence and proliferation of these variants, consistent with models of rapid adaptation and evolutionary rescue.

Detailing how populations adapted to environmental change is needed to predict future responses, but identifying adaptive variants and detailing their fitness effects is rare. Here, the authors show that parallel loss of FRI and FLC function reduces time to flowering and drives adaptation in a drought prone environment.

A cross-scale approach to unravel the molecular basis of plant phenology in temperate and tropical climates

Plants have evolved to time their leafing, flowering and fruiting in appropriate seasons for growth, reproduction and resting. As a consequence of their adaptation to geographically different environments, there is a rich diversity in plant phenology from temperate and tropical climates. Recent progress in genetic and molecular studies will provide numerous opportunities to study the genetic basis of phenological traits and the history of adaptation of phenological traits to seasonal and aseasonal environments. Integrating molecular data with long-term phenology and climate data into predictive models will be a powerful tool to forecast future phenological changes in the face of global environmental change. Here, we review the cross-scale approach from genes to plant communities from three aspects: the latitudinal gradient of plant phenology at the community level, the environmental and genetic factors underlying the diversity of plant phenology, and an integrated approach to forecast future plant phenology based on genetically informed knowledge. Synthesizing the latest knowledge about plant phenology from molecular, ecological and mathematical perspectives will help us understand how natural selection can lead to the further evolution of the gene regulatory mechanisms in phenological traits in future forest ecosystems.

Functional homoeologous alleles of CONSTANS contribute to seasonal crop type in rapeseed

Two CO paralogs in Brassica napus were confirmed and shown distinct expression pattern and function in promoting flowering and allelic variation s within BnaCO.A10 were found closely associated with ecotype divergence. CONSTANS (CO) is a key gene that responds to photoperiod and in Arabidopsis can promote flowering under long-day (LD) conditions. Brassica napus L. is a major oil crop and close relative of Arabidopsis, and arose via allopolyploidization from the diploids B. rapa (A genome) and B. oleracea (C genome). In this study, we confirmed that B. napus has two CO genes located on the A10 (BnaCO.A10) and C9 (BnaCO.C9) chromosomes. Significant differences in level and temporal pattern of transcription, as well as in protein function, of these homoeologous may have resulted from sequence variation in the promoter as well as in the coding region. Apart from two insertions of 527 bp and 2002 bp in the promoter of BnaCO.C9 that function as transcriptional enhancers, this gene is otherwise highly conserved in both promoter and coding region. However, BnaCO.A10 was classified into two haplotypes and transgene analysis in Arabidopsis and backcross analysis in rapeseed indicated that the winter-type haplotype had a greater effect in promoting flowering than the spring type. We discuss the contribution of CO alleles to species evolution, and for eco-geographic radiation following crop domestication, alongside scope for managing this locus in future breeding.

Association between RsFT, RsFLC and RsCOL5 (A&B) expression and flowering regulation in Japanese wild radish

Flowering is an important step in the life cycle of plants and indicates adaptability to external climatic cues such as temperature and photoperiod. We investigated the expression patterns of core genes related to flowering-time regulation in Japanese wild radish (Raphanus sativus var. raphanistroides) with different vernalization requirements (obligate and facultative) and further identified climatic cues that may act as natural selective forces. Specifically, we analysed flowering-time variation under different cold and photoperiod treatments in Japanese wild radish collected from the Hokkaido (northern lineage) and Okinawa (southern lineage) islands, which experience contrasting climatic cues. The cultivation experiment verified the obligate and facultative vernalization requirements of the northern and southern wild radish accessions, respectively. The expression of major genes involved in flowering time indicated that RsFLC and RsCOL5 (A&B) may interact to regulate flowering time. Notably, floral initiation in the northern lineage was strongly correlated with RsFLC expression, whereas flowering in the southern linage was correlated with induction of RsCOL5-A expression, despite high RsFLC transcript levels. These results suggested that the northern accessions are more sensitive to prolonged cold exposure, whereas the southern accessions are more sensitive to photoperiod. These different mechanisms ultimately confer an optimal flowering time in natural populations in response to locally contrasting climatic cues. This study provides new insights into the variant mechanisms underlying floral pathways in Japanese wild radish from different geographic locations.

Exploring PIF4 's contribution to early flowering in plants under daily variable temperature and its tissue-specific flowering gene network

Molecular mechanisms of how constant temperatures affect flowering time have been largely characterized in the model plant Arabidopsis thaliana; however, the effect of natural daily variable temperature outside laboratories is only partly explored. Several flowering genes have been shown to play important roles in temperature responses, including PHYTOCHROME-INTERACTING FACTOR 4 (PIF4) and FLOWERING LOCUS C (FLC), the two genes encoding for the transcription factors (TFs) that act antagonistically to regulate flowering time by activating and repressing floral integrator FLOWERING LOCUS T (FT), respectively. In this study, we have taken a multidisciplinary approach to explore the contribution of PIF4 to the early flowering observed in the daily variable temperature (VAR) and to broaden its transcriptional network using publicly available transcriptomic data. We observed early flowering in the natural accessions Col-0, C24 and their late flowering hybrid C24xCol grown under VAR, as compared with a constant temperature (CON). The loss-of-function mutation of PIF4 exhibits later flowering in VAR in both the Col-0 parent and the C24xCol hybrid, suggesting that PIF4, at least in part, contributes to acceleration of flowering in the VAR condition. To investigate the interplay between PIF4 and its flowering regulator counterparts, FLC and FT, we performed transcriptional analyses and found that VAR increased PIF4 transcription at the end of the day when temperature peaked at 32°C, when FT transcription was also elevated. On the other hand, we observed a decrease in FLC transcription in the 4-week-old plants grown in VAR, as well as in the plants with PIF4 overexpression grown in CON. These results raise a possibility that PIF4 might also regulate FT indirectly through the repression of FLC, in addition to the well-characterized direct control of PIF4 over FT. To further expand our view on the PIF4-orientated flowering gene network in response to temperature changes, we have constructed a coexpression-transcriptional regulatory network by combining publicly available transcriptomic data and gene regulatory interactions of PIF4 and its closely related flowering genes, PIF5, FLC, and ELF3. The network model reveals conserved and tissue-specific regulatory functions, which are useful for confirming as well as predicting the functions and regulatory interactions between these key flowering genes.

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.

A New Catalog of Structural Variants in 1,301 A. thaliana Lines from Africa, Eurasia, and North America Reveals a Signature of Balancing Selection at Defense Response Genes

Genomic variation in the model plant Arabidopsis thaliana has been extensively used to understand evolutionary processes in natural populations, mainly focusing on single-nucleotide polymorphisms. Conversely, structural variation has been largely ignored in spite of its potential to dramatically affect phenotype. Here, we identify 155,440 indels and structural variants ranging in size from 1 bp to 10 kb, including presence/absence variants (PAVs), inversions, and tandem duplications in 1,301 A. thaliana natural accessions from Morocco, Madeira, Europe, Asia, and North America. We show evidence for strong purifying selection on PAVs in genes, in particular for housekeeping genes and homeobox genes, and we find that PAVs are concentrated in defense-related genes (R-genes, secondary metabolites) and F-box genes. This implies the presence of a "core" genome underlying basic cellular processes and a "flexible" genome that includes genes that may be important in spatially or temporally varying selection. Further, we find an excess of intermediate frequency PAVs in defense response genes in nearly all populations studied, consistent with a history of balancing selection on this class of genes. Finally, we find that PAVs in genes involved in the cold requirement for flowering (vernalization) and drought response are strongly associated with temperature at the sites of origin.

Quantitative Control of Early Flowering in White Lupin (Lupinus albus L.)

White lupin (Lupinus albus L.) is a pulse annual plant cultivated from the tropics to temperate regions for its high-protein grain as well as a cover crop or green manure. Wild populations are typically late flowering and have high vernalization requirements. Nevertheless, some early flowering and thermoneutral accessions were found in the Mediterranean basin. Recently, quantitative trait loci (QTLs) explaining flowering time variance were identified in bi-parental population mapping, however, phenotypic and genotypic diversity in the world collection has not been addressed yet. In this study, a diverse set of white lupin accessions (n = 160) was phenotyped for time to flowering in a controlled environment and genotyped with PCR-based markers (n = 50) tagging major QTLs and selected homologs of photoperiod and vernalization pathway genes. This survey highlighted quantitative control of flowering time in white lupin, providing statistically significant associations for all major QTLs and numerous regulatory genes, including white lupin homologs of CONSTANS, FLOWERING LOCUS T, FY, MOTHER OF FT AND TFL1, PHYTOCHROME INTERACTING FACTOR 4, SKI-INTERACTING PROTEIN 1, and VERNALIZATION INDEPENDENCE 3. This revealed the complexity of flowering control in white lupin, dispersed among numerous loci localized on several chromosomes, provided economic justification for future genome-wide association studies or genomic selection rather than relying on simple marker-assisted selection.

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.

Naturally occurring circadian rhythm variation associated with clock gene loci in Swedish Arabidopsis accessions

Circadian clocks have evolved to resonate with external day and night cycles. However, these entrainment signals are not consistent everywhere and vary with latitude, climate and seasonality. This leads to divergent selection for clocks which are locally adapted. To investigate the genetic basis for this circadian variation, we used a delayed fluorescence imaging assay to screen 191 naturally occurring Swedish Arabidopsis accessions for their circadian phenotypes. We demonstrate that the period length co-varies with both geography and population sub-structure. Several candidate loci linked to period, phase and relative amplitude error (RAE) were revealed by genome-wide association mapping and candidate genes were investigated using TDNA mutants. We show that natural variation in a single non-synonymous substitution within COR28 is associated with a long-period and late-flowering phenotype similar to that seen in TDNA knock-out mutants. COR28 is a known coordinator of flowering time, freezing tolerance and the circadian clock; all of which may form selective pressure gradients across Sweden. We demonstrate the effect of the COR28-58S SNP in increasing period length through a co-segregation analysis. Finally, we show that period phenotypic tails remain diverged under lower temperatures and follow a distinctive "arrow-shaped" trend indicative of selection for a cold-biased temperature compensation response.

Polygenic adaptation of rosette growth in Arabidopsis thaliana

The rate at which plants grow is a major functional trait in plant ecology. However, little is known about its evolution in natural populations. Here, we investigate evolutionary and environmental factors shaping variation in the growth rate of Arabidopsis thaliana. We used plant diameter as a proxy to monitor plant growth over time in environments that mimicked latitudinal differences in the intensity of natural light radiation, across a set of 278 genotypes sampled within four broad regions, including an outgroup set of genotypes from China. A field experiment conducted under natural conditions confirmed the ecological relevance of the observed variation. All genotypes markedly expanded their rosette diameter when the light supply was decreased, demonstrating that environmental plasticity is a predominant source of variation to adapt plant size to prevailing light conditions. Yet, we detected significant levels of genetic variation both in growth rate and growth plasticity. Genome-wide association studies revealed that only 2 single nucleotide polymorphisms associate with genetic variation for growth above Bonferroni confidence levels. However, marginally associated variants were significantly enriched among genes with an annotated role in growth and stress reactions. Polygenic scores computed from marginally associated variants confirmed the polygenic basis of growth variation. For both light regimes, phenotypic divergence between the most distantly related population (China) and the various regions in Europe is smaller than the variation observed within Europe, indicating that the evolution of growth rate is likely to be constrained by stabilizing selection. We observed that Spanish genotypes, however, reach a significantly larger size than Northern European genotypes. Tests of adaptive divergence and analysis of the individual burden of deleterious mutations reveal that adaptive processes have played a more important role in shaping regional differences in rosette growth than maladaptive evolution.

The rate at which plants grow is a major functional trait in plant ecology. However, little is known about its genetic variation in natural populations. Here, we investigate genetic and environmental factors shaping variation in the growth rate of Arabidopsis thaliana and ask whether genetic variation in plant growth contributes to adaptation to local environmental conditions. We grew plants under two light regimes that mimic latitudinal differences in the intensity of natural light radiation, and measured plant diameter as it grew over time. When the light supply was decreased, plant diameter grew more slowly but reached a markedly larger final size, confirming that plants can adjust their growth to prevailing light conditions. Yet, we also detected significant levels of genetic variation both in growth rate and in how the growth dynamics is adjusted to the light conditions. We show that this variation is encoded by many loci of small effect that are hard to locate in the genome but overall significantly enriched among genes associated with growth and stress reactions. We further observe that Spanish genotypes tended to reach, on average, a significantly larger rosette size than Northern European genotypes. Tests of adaptive divergence indicate that these differences may reflect adaptation to local environmental conditions.

Gene regulatory networks controlled by FLOWERING LOCUS C that confer variation in seasonal flowering and life history

Responses to environmental cues synchronize reproduction of higher plants to the changing seasons. The genetic basis of these responses has been intensively studied in the Brassicaceae. The MADS-domain transcription factor FLOWERING LOCUS C (FLC) plays a central role in the regulatory network that controls flowering of Arabidopsis thaliana in response to seasonal cues. FLC blocks flowering until its transcription is stably repressed by extended exposure to low temperatures in autumn or winter and, therefore, FLC activity is assumed to limit flowering to spring. Recent reviews describe the complex epigenetic mechanisms responsible for FLC repression in cold. We focus on the gene regulatory networks controlled by FLC and how they influence floral transition. Genome-wide approaches determined the in vivo target genes of FLC and identified those whose transcription changes during vernalization or in flc mutants. We describe how studying FLC targets such as FLOWERING LOCUS T, SQUAMOSA PROMOTER BINDING PROTEIN-LIKE 15, and TARGET OF FLC AND SVP 1 can explain different flowering behaviours in response to vernalization and other environmental cues, and help define mechanisms by which FLC represses gene transcription. Elucidating the gene regulatory networks controlled by FLC provides access to the developmental and physiological mechanisms that regulate floral transition.

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

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

Less Is More, Natural Loss-of-Function Mutation Is a Strategy for Adaptation

Gene gain and loss are crucial factors that shape the evolutionary success of diverse organisms. In the past two decades, more attention has been paid to the significance of gene gain through gene duplication or de novo genes. However, gene loss through natural loss-of-function (LoF) mutations, which is prevalent in the genomes of diverse organisms, has been largely ignored. With the development of sequencing techniques, many genomes have been sequenced across diverse species and can be used to study the evolutionary patterns of gene loss. In this review, we summarize recent advances in research on various aspects of LoF mutations, including their identification, evolutionary dynamics in natural populations, and functional effects. In particular, we discuss how LoF mutations can provide insights into the minimum gene set (or the essential gene set) of an organism. Furthermore, we emphasize their potential impact on adaptation. At the genome level, although most LoF mutations are neutral or deleterious, at least some of them are under positive selection and may contribute to biodiversity and adaptation. Overall, we highlight the importance of natural LoF mutations as a robust framework for understanding biological questions in general.

Multiple Loci Control Variation in Plasticity to Foliar Shade Throughout Development in Arabidopsis thaliana

The shade avoidance response is a set of developmental changes exhibited by plants to avoid shading by competitors, and is an important model of adaptive plant plasticity. While the mechanisms of sensing shading by other plants are well-known and appear conserved across plants, less is known about the developmental mechanisms that result in the diverse array of morphological and phenological responses to shading. This is particularly true for traits that appear later in plant development. Here we use a nested association mapping (NAM) population of Arabidopsis thaliana to decipher the genetic architecture of the shade avoidance response in late-vegetative and reproductive plants. We focused on four traits: bolting time, rosette size, inflorescence growth rate, and inflorescence size, found plasticity in each trait in response to shade, and detected 17 total QTL; at least one of which is a novel locus not previously identified for shade responses in Arabidopsis. Using path analysis, we dissected each colocalizing QTL into direct effects on each trait and indirect effects transmitted through direct effects on earlier developmental traits. Doing this separately for each of the seven NAM populations in each environment, we discovered considerable heterogeneity among the QTL effects across populations, suggesting allelic series at multiple QTL or interactions between QTL and the genetic background or the environment. Our results provide insight into the development and variation in shade avoidance responses in Arabidopsis, and emphasize the value of directly modeling the relationships among traits when studying the genetics of complex developmental syndromes.

Global insights into duplicated gene expression and alternative splicing in polyploid Brassica napus under heat, cold, and drought stress

Polyploidy has been a prevalent process during plant evolution and it has made a major impact on the structure and evolution of plant genomes. Many important crop plants are polyploid. There is considerable interest in expression patterns of duplicated genes in polyploids. Alternative splicing (AS) is a fundamental aspect of gene expression that produces multiple final transcript types from a single type of mRNAs. The effects of abiotic stress conditions on AS in polyploids has received little attention. We conducted a global transcriptome analysis of Brassica napus, an allotetraploid derived from B. rapa (A_T ) and B. oleracea (C_T ), by RNA-Seq of plants subjected to cold, heat, and drought stress treatments. Analyses of 27,360 pairs of duplicated genes revealed overall A_T subgenome biases in gene expression and C_T subgenome biases in the extent of alternative splicing under all three stress treatments. More genes increased in expression than decreased in response to the stresses. Negative correlations were found between expression levels and AS frequency for each type of AS. Cold stress produced the greatest changes in gene expression and AS. Cold-induced AS changes were more likely to be shared with those generated by drought than by heat stress. We used homeologs of FLC and CCA1 as case studies to show the dynamics of how duplicates in a polyploid respond to cold stress. Our results suggest that divergence in gene expression and AS patterns between duplicated genes may increase the flexibility of polyploids when responding to abiotic stressors.

The molecular mechanism of vernalization in Arabidopsis and cereals: role of Flowering Locus C and its homologs

Winter varieties of plants can flower only after exposure to prolonged cold. This phenomenon is known as vernalization and has been widely studied in the model plant Arabidopsis thaliana as well as in monocots. Through the repression of floral activator genes, vernalization prevents flowering in winter. In Arabidopsis, FLOWERING LOCUS C or FLC is the key repressor during vernalization, while in monocots vernalization is regulated through VRN1, VRN2 and VRN3 (or FLOWERING LOCUS T). Interestingly, VRN genes are not homologous to FLC but FLC homologs are found to have a significant role in vernalization response in cereals. The presence of FLC homologs in monocots opens new dimensions to understand, compare and retrace the evolution of vernalization pathways between monocots and dicots. In this review, we discuss the molecular mechanism of vernalization-induced flowering along with epigenetic regulations in Arabidopsis and temperate cereals. A better understanding of cold-induced flowering will be helpful in crop breeding strategies to modify the vernalization requirement of economically important temperate cereals.

Maternal transmission of the epigenetic 'memory of winter cold' in Arabidopsis

Some plants can 'remember' past environmental experience to become adapted to a given environment. For instance, after experiencing prolonged low-temperature exposure in winter (winter cold), vernalization-responsive plants remember past cold experience when temperature rises in spring, to acquire competence to flower at a later season favourable for seed production^1,2. In Arabidopsis thaliana, prolonged cold induces silencing of the potent floral repressor FLOWERING LOCUS C (FLC) by Polycomb group (PcG) chromatin modifiers. This Polycomb-repressed chromatin state is epigenetically maintained and thus 'memorized' in subsequent growth and development upon return to warmth^1,3. 'Memory of winter cold' has been viewed as being mitotically stable but meiotically unstable^3-5, and thus not to be transmitted intergenerationally. In general, whether and how chromatin-mediated environmental memories are transmitted across generations are unknown in plants. Here, we show that the cold-induced Polycomb-repressed chromatin state at FLC or memory of winter cold is maintained in the egg cell, that is meiotically stable in the process of female gamete formation, and provide evidence that this Polycomb-mediated memory is not maintained in the sperm cell. Moreover, we show that this cold memory is inherited maternally but not paternally to the zygote and early embryos. Our study demonstrates and further provides mechanistic insights into intergenerational transmission of chromatin state-mediated environmental memories in plants.