Diversity of woodland strawberry inflorescences arises from heterochrony regulated by TERMINAL FLOWER 1 and FLOWERING LOCUS T.

A vast variety of inflorescence architectures have evolved in angiosperms. Here, we analyze the diversity and development of the woodland strawberry (Fragaria vesca) inflorescence. Contrary to historical classifications, we show that it is a closed thyrse: a compound inflorescence with determinate primary monopodial axis and lateral sympodial branches, thus combining features of racemes and cymes. We demonstrate that this architecture is generated by 2 types of inflorescence meristems differing in their geometry. We further show that woodland strawberry homologs of TERMINAL FLOWER 1 (FvTFL1) and FLOWERING LOCUS T (FvFT1) regulate the development of both the racemose and cymose components of the thyrse. Loss of functional FvTFL1 reduces the number of lateral branches of the main axis and iterations in the lateral branches but does not affect their cymose pattern. These changes can be enhanced or compensated by altering FvFT1 expression. We complement our experimental findings with a computational model that captures inflorescence development using a small set of rules. The model highlights the distinct regulation of the fate of the primary and higher-order meristems, and explains the phenotypic diversity among inflorescences in terms of heterochrony resulting from the opposite action of FvTFL1 and FvFT1 within the thyrse framework. Our results represent a detailed analysis of thyrse architecture development at the meristematic and molecular levels.

Diversity of the volatilome and the fruit size and shape in European woodland strawberry (Fragaria vesca).

Woodland strawberry (Fragaria vesca subsp. vesca) is a wild relative of cultivated strawberry (F. × ananassa) producing small and typically conical fruits with an intense flavor and aroma. The wild strawberry species, F. vesca, is a rich resource of genetic and metabolic variability, but its diversity remains largely unexplored and unexploited. In this study, we aim for an in-depth characterization of the fruit complex volatilome by GC-MS as well as the fruit size and shape using a European germplasm collection that represents the continental diversity of the species. We report characteristic volatilome footprints and fruit phenotypes of specific geographical areas. Thus, this study uncovers phenotypic variation linked to geographical distribution that will be valuable for further genetic studies to identify candidate genes or develop markers linked to volatile compounds or fruit shape and size traits.

Multi-model GWAS reveals key loci for horticultural traits in reconstructed garden strawberry.

The cultivated garden strawberry (Fragaria × ananassa) has a rich history, originating from the hybridization of two wild octoploid strawberry species in the 18th century. Two-step reconstruction of Fragaria × ananassa through controlled crossings between pre-improved selections of its parental species is a promising approach for enriching the breeding germplasm of strawberry for wider adaptability. We created a population of reconstructed strawberry by hybridizing elite selections of F. virginiana and F. chiloensis. A replicated field experiment was conducted to evaluate the population's performance for eleven horticulturally important traits, over multiple years. Population structure analyses based on Fana-50 k SNP array data confirmed pedigree-based grouping of the progenies into four distinct groups. As complex traits are often influenced by environmental variables, and population structure can lead to spurious associations, we tested multiple genome-wide association study (GWAS) models. GWAS uncovered 39 quantitative trait loci (QTL) regions for eight traits distributed across twenty chromosomes, including 11 consistent and 28 putative QTLs. Candidate genes for traits including winter survival, flowering time, runnering vigor, and hermaphrodism were identified within the QTL regions. To our knowledge, this study marks the first comprehensive investigation of adaptive and horticultural traits in a large, multi-familial reconstructed strawberry population using SNP markers.

The GATA factor HANABA TARANU promotes runner formation by regulating axillary bud initiation and outgrowth in cultivated strawberry.

A runner, as an elongated branch, develops from the axillary bud (AXB) in the leaf axil and is crucial for the clonal propagation of cultivated strawberry (Fragaria × ananassa Duch.). Runner formation occurs in at least two steps: AXB initiation and AXB outgrowth. HANABA TARANU (HAN ) encodes a GATA transcription factor that affects AXB initiation in Arabidopsis and promotes branching in grass species, but the underlying mechanism is largely unknown. Here, the function of a strawberry HAN homolog FaHAN in runner formation was characterized. FaHAN transcripts can be detected in the leaf axils. Overexpression (OE) of FaHAN increased the number of runners, mainly by enhancing AXB outgrowth, in strawberry. The expression of the strawberry homolog of BRANCHED1 , a key inhibitor of AXB outgrowth in many plant species, was significantly downregulated in the AXBs of FaHAN -OE lines, whereas the expression of the strawberry homolog of SHOOT MERISTEMLESS, a marker gene for AXB initiation in Arabidopsis, was upregulated. Moreover, several genes of gibberellin biosynthesis and cytokinin signaling pathways were activated, whereas the auxin response pathway genes were repressed. Further assays indicated that FaHAN could be directly activated by FaNAC2, the overexpression of which in strawberry also increased the number of runners. The silencing of FaNAC2 or FaHAN inhibited AXB initiation and led to a higher proportion of dormant AXBs, confirming their roles in the control of runner formation. Taken together, our results revealed a FaNAC2-FaHAN pathway in the control of runner formation and have provided a means to enhance the vegetative propagation of cultivated strawberry.

Acclimation of circadian rhythms in woodland strawberries (Fragaria vesca L.) to Arctic and mid-latitude photoperiods.

BACKGROUND: Though many abiotic factors are constantly changing, the photoperiod is a predictable factor that enables plants to time many physiological responses. This timing is regulated by the circadian clock, yet little is known about how the clock adapts to the differences in photoperiod between mid-latitudes and high latitudes. The primary objective of this study was to compare how clock gene expression is modified in four woodland strawberry (Fragaria vesca L.) accessions originating from two different populations in Italy (IT1: Tenno, Italy, 45°N, IT4: Salorno, Italy, 46°N) and two in Northern Norway (NOR2: Alta, Norway, 69°N, NOR13: Indre Nordnes, Norway 69°N) when grown under simulated daylength conditions of an Arctic or mid-latitude photoperiod. The second objective was to investigate whether population origin or the difference in photoperiod influenced phytohormone accumulation. RESULTS: The Arctic photoperiod induced lower expression in IT4 and NOR13 for six clock genes (FvLHY, FvRVE8, FvPRR9, FvPRR7, FvPRR5, and FvLUX), in IT1 for three genes (FvLHY, FvPRR9, and FvPRR5) and in NOR2 for one gene (FvPRR9). Free-running rhythms for FvLHY in IT1 and IT4 were higher after the Arctic photoperiod, while the free-running rhythm for FvLUX in IT4 was higher after the mid-latitude photoperiod. IT1 showed significantly higher expression of FvLHY and FvPRR9 than all other accessions, as well as significantly higher expression of the circadian regulated phytohormone, abscisic acid (ABA), but low levels of salicylic acid (SA). NOR13 had significantly higher expression of FvRVE8, FvTOC1, and FvLUX than all other accessions. NOR2 had extremely low levels of auxin (IAA) and high levels of the jasmonate catabolite, hydroxyjasmonic acid (OH-JA). CONCLUSIONS: Our study shows that circadian rhythms in Fragaria vesca are driven by both the experienced photoperiod and genetic factors, while phytohormone levels are primarily determined by specific accessions' genetic factors rather than the experienced photoperiod.

Two MYB activators of anthocyanin biosynthesis exhibit specialized activities in petiole and fruit of diploid strawberry.

The R2R3-MYB transcription factor FveMYB10 is a major regulator of anthocyanin pigmentation in the red fruits of strawberry. fvemyb10 loss-of-function mutants form yellow fruits but still accumulate purple-colored anthocyanins in the petioles, suggesting that anthocyanin biosynthesis is under distinct regulation in fruits and petioles. From chemical mutagenesis in the diploid wild strawberry Fragaria vesca, we identified a green petioles (gp)-1 mutant that lacks anthocyanins in petioles. Using mapping-by-sequencing and transient functional assays, we confirmed that the causative mutation resides in a FveMYB10-Like (FveMYB10L) gene and that FveMYB10 and FveMYB10L function independently in the fruit and petiole, respectively. In addition to their tissue-specific regulation, FveMYB10 and FveMYB10L respond differently to changes in light quality, produce distinct anthocyanin compositions, and preferentially activate different downstream anthocyanin biosynthesis genes in their respective tissues. This work identifies a new regulator of anthocyanin synthesis and demonstrates that two paralogous MYB genes with specialized functions enable tissue-specific regulation of anthocyanin biosynthesis in fruit and petiole tissues.

Allelic Variation of MYB10 Is the Major Force Controlling Natural Variation in Skin and Flesh Color in Strawberry (Fragaria spp.) Fruit.

The fruits of diploid and octoploid strawberry (Fragaria spp) show substantial natural variation in color due to distinct anthocyanin accumulation and distribution patterns. Anthocyanin biosynthesis is controlled by a clade of R2R3 MYB transcription factors, among which MYB10 is the main activator in strawberry fruit. Here, we show that mutations in MYB10 cause most of the variation in anthocyanin accumulation and distribution observed in diploid woodland strawberry (F. vesca) and octoploid cultivated strawberry (F ×ananassa). Using a mapping-by-sequencing approach, we identified a gypsy-transposon in MYB10 that truncates the protein and knocks out anthocyanin biosynthesis in a white-fruited F. vesca ecotype. Two additional loss-of-function mutations in MYB10 were identified among geographically diverse white-fruited F. vesca ecotypes. Genetic and transcriptomic analyses of octoploid Fragaria spp revealed that FaMYB10-2, one of three MYB10 homoeologs identified, regulates anthocyanin biosynthesis in developing fruit. Furthermore, independent mutations in MYB10-2 are the underlying cause of natural variation in fruit skin and flesh color in octoploid strawberry. We identified a CACTA-like transposon (FaEnSpm-2) insertion in the MYB10-2 promoter of red-fleshed accessions that was associated with enhanced expression. Our findings suggest that cis-regulatory elements in FaEnSpm-2 are responsible for enhanced MYB10-2 expression and anthocyanin biosynthesis in strawberry fruit flesh.

Methylome, transcriptome, and phenotype changes induced by temperature conditions experienced during sexual reproduction in Fragaria vesca.

Temperature conditions experienced during embryogenesis and seed development may induce epigenetic changes that increase phenotypic variation in plants. Here we investigate if embryogenesis and seed development at two different temperatures (28 vs. 18°C) result in lasting phenotypic effects and DNA methylation changes in woodland strawberry (Fragaria vesca). Using five European ecotypes from Spain (ES12), Iceland (ICE2), Italy (IT4), and Norway (NOR2 and NOR29), we found statistically significant differences between plants from seeds produced at 18 or 28°C in three of four phenotypic features investigated under common garden conditions. This indicates the establishment of a temperature-induced epigenetic memory-like response during embryogenesis and seed development. The memory effect was significant in two ecotypes: in NOR2 flowering time, number of growth points and petiole length were affected, and in ES12 number of growth points was affected. This indicates that genetic differences between ecotypes in their epigenetic machinery, or other allelic differences, impact this type of plasticity. We observed statistically significant differences between ecotypes in DNA methylation marks in repetitive elements, pseudogenes, and genic elements. Leaf transcriptomes were also affected by embryonic temperature in an ecotype-specific manner. Although we observed significant and lasting phenotypic change in at least some ecotypes, there was considerable variation in DNA methylation between individual plants within each temperature treatment. This within-treatment variability in DNA methylation marks in F. vesca progeny may partly be a result of allelic redistribution from recombination during meiosis and subsequent epigenetic reprogramming during embryogenesis.

Signatures of polygenic adaptation align with genome-wide methylation patterns in wild strawberry plants.

Epigenetic inheritance can drive adaptive evolution independently of DNA sequence variation. However, to what extent epigenetic variation represents an autonomous evolutionary force remains largely elusive. Through gene ontology and comparative analyses of genomic and epigenomic variation of wild strawberry plants raised in distinct drought settings, we characterised genome-wide covariation between single nucleotide polymorphisms (SNPs) and differentially methylated cytosines (DMCs). Covariation between SNPs and DMCs was independent of genomic proximity, but instead associated with fitness-related processes such as stress responses, genome regulation and reproduction. We expected this functional SNP-DMC covariation to be driven by adaptive evolution canalising SNP and DMC variation, but instead observed significantly lower covariation with DMCs for adaptive rather than for neutral SNPs. Drought-induced DMCs frequently co-varied with tens of SNPs, suggesting high genomic redundancy as a broad potential basis for polygenic adaptation of gene expression. Our findings suggest that stress-responsive DMCs initially co-vary with many SNPs under increased environmental stress, and that natural selection acting upon several of these SNPs subsequently reduces standing covariation with stress-responsive DMCs. Our study supports DNA methylation profiles that represent complex quantitative traits rather than autonomous evolutionary forces. We provide a conceptual framework for polygenic regulation and adaptation shaping genome-wide methylation patterns in plants.

Woodland strawberry axillary bud fate is dictated by a crosstalk of environmental and endogenous factors.

Plant architecture is defined by fates and positions of meristematic tissues and has direct consequences on yield potential and environmental adaptation of the plant. In strawberries (Fragaria vesca L. and F. × ananassa Duch.), shoot apical meristems can remain vegetative or differentiate into a terminal inflorescence meristem. Strawberry axillary buds (AXBs) are located in leaf axils and can either remain dormant or follow one of the two possible developmental fates. AXBs can either develop into stolons needed for clonal reproduction or into branch crowns (BCs) that can bear their own terminal inflorescences under favorable conditions. Although AXB fate has direct consequences on yield potential and vegetative propagation of strawberries, the regulation of AXB fate has so far remained obscure. We subjected a number of woodland strawberry (F. vesca L.) natural accessions and transgenic genotypes to different environmental conditions and growth regulator treatments to demonstrate that strawberry AXB fate is regulated either by environmental or endogenous factors, depending on the AXB position on the plant. We confirm that the F. vesca GIBBERELLIN20-oxidase4 (FvGA20ox4) gene is indispensable for stolon development and under tight environmental regulation. Moreover, our data show that apical dominance inhibits the outgrowth of the youngest AXB as BCs, although the effect of apical dominance can be overrun by the activity of FvGA20ox4. Finally, we demonstrate that the FvGA20ox4 is photoperiodically regulated via FvSOC1 (F. vesca SUPPRESSOR OF OVEREXPRESSION OF CONSTANS1) at 18°C, but at higher temperature of 22°C an unidentified FvSOC1-independent pathway promotes stolon development.

Influence of Arctic light conditions on crop production and quality.

The natural light conditions above the Arctic Circle are unique in terms of annual variation creating special growth conditions for crop production. These include low solar elevations, very long daily photosynthetic light periods, midnight sun/absence of dark nights, and altered spectral distribution depending on solar elevation. All these factors are known to affect the growth and the metabolism of plants, although their influence on northern crop plants has not yet been reviewed. The ongoing global warming is especially affecting the temperature × light interactions in the Arctic, and understanding the impact on crop production and plant metabolism will be important for an Arctic contribution to global food production. Arctic light conditions have a strong influence on the timing of plant development, which together with temperature limits the number of cultivars suitable for Arctic agriculture. This review compiles information from the reports about the effects of light conditions at high latitudes on growth, biomass production, flowering and quality of the crop plants and discusses the gained knowledge and the key gaps to be addressed.

Altered regulation of TERMINAL FLOWER 1 causes the unique vernalisation response in an arctic woodland strawberry accession.

Vernalisation requirement is an agriculturally important trait that postpones the development of cold-sensitive floral organs until the spring. The family Rosaceae includes many agriculturally important fruit and berry crops that suffer from crop losses caused by frost injury to overwintering flower buds. Recently, a vernalisation-requiring accession of the Rosaceae model woodland strawberry (Fragaria vesca) has been identified in northern Norway. Understanding the molecular basis of the vernalisation requirement in this accession would advance the development of strawberry cultivars better adapted to temperate climate. We use gene silencing, gene expression analysis, genetic mapping and population genomics to study the genetic basis of the vernalisation requirement in woodland strawberry. Our results indicate that the woodland strawberry vernalisation requirement is endemic to northern Norwegian population, and mapping data suggest the orthologue of TERMINAL FLOWER1 (FvTFL1) as the causal floral repressor. We demonstrate that exceptionally low temperatures are needed to downregulate FvTFL1 and to make these plants competent to induce flowering at low postvernalisation temperatures in the spring. We show that altered regulation of FvTFL1 in the northern Norwegian woodland strawberry accession postpones flower induction until the spring, allowing plants to avoid winter injuries of flower buds that commonly occur in temperate regions.

Fragaria vesca CONSTANS controls photoperiodic flowering and vegetative development.

According to the external coincidence model, photoperiodic flowering occurs when CONSTANS (CO) mRNA expression coincides with light in the afternoon of long days (LDs), leading to the activation of FLOWERING LOCUS T (FT). CO has evolved in Brassicaceae from other Group Ia CO-like (COL) proteins which do not control photoperiodic flowering in Arabidopsis. COLs in other species have evolved different functions as floral activators or even as repressors. To understand photoperiodic development in the perennial rosaceous model species woodland strawberry, we functionally characterized FvCO, the only Group Ia COL in its genome. We demonstrate that FvCO has a major role in the photoperiodic control of flowering and vegetative reproduction through runners. FvCO is needed to generate a bimodal rhythm of FvFT1 which encodes a floral activator in the LD accession Hawaii-4: a sharp FvCO expression peak at dawn is followed by the FvFT1 morning peak in LDs indicating possible direct regulation, but additional factors that may include FvGI and FvFKF1 are probably needed to schedule the second FvFT1 peak around dusk. These results demonstrate that although FvCO and FvFT1 play major roles in photoperiodic development, the CO-based external coincidence around dusk is not fully applicable to the woodland strawberry.

Strawberry homologue of terminal flower1 integrates photoperiod and temperature signals to inhibit flowering.

Photoperiod and temperature are major environmental signals affecting flowering in plants. Although molecular pathways mediating these signals have been well characterized in the annual model plant Arabidopsis, much less information is known in perennials. Many perennials including the woodland strawberry (Fragaria vesca L.) are induced to flower in response to decreasing photoperiod and temperature in autumn and they flower following spring. We showed earlier that, in contrast with Arabidopsis, the photoperiodic induction of flowering in strawberry occurs in short days (SD) when the decrease in FvFT1 (flowering locus T) and FvSOC1 (suppressor of the overexpression of constans1) expression leads to lower mRNA levels of the floral repressor, FvTFL1 (terminal flower1). By using transgenic lines and gene expression analyses, we show evidence that the temperature-mediated changes in the FvTFL1 mRNA expression set critical temperature limits for the photoperiodic flowering in strawberry. At temperatures below 13 °C, low expression level of FvTFL1 in both SD and long days (LD) allows flower induction to occur independently of the photoperiod. Rising temperature gradually increases FvTFL1 mRNA levels under LD, and at temperatures above 13 °C, SD is required for the flower induction that depends on the deactivation of FvSOC1 and FvTFL1. However, an unknown transcriptional activator, which functions independently of FvSOC1, enhances the expression of FvTFL1 at 23 °C preventing photoperiodic flowering. We suggest that the observed effect of the photoperiod × temperature interaction on FvTFL1 mRNA expression may allow strawberry to induce flowers in correct time in different climates.

TERMINAL FLOWER1 is a breeding target for a novel everbearing trait and tailored flowering responses in cultivated strawberry (Fragaria × ananassa Duch.).

The effects of daylength and temperature on flowering of the cultivated octoploid strawberry (Fragaria × ananassa Duch.) have been studied extensively at the physiological level, but information on the molecular pathways controlling flowering in the species is scarce. The flowering pathway has been studied at the molecular level in the diploid short-day woodland strawberry (F. vesca L.), in which the FLOWERING LOCUS T1 (FvFT1)-SUPPRESSOR OF OVEREXPRESSION OF CONSTANS1 (FvSOC1)-TERMINAL FLOWER1 (FvTFL1) pathway is essential for the correct timing of flowering. In this work, we show by transgenic approach that the silencing of the floral repressor FaTFL1 in the octoploid short-day cultivar 'Elsanta' is sufficient to induce perpetual flowering under long days without direct changes in vegetative reproduction. We also demonstrate that although the genes FaFT1 and FaSOC1 show similar expression patterns in different cultivars, the regulation of FaTFL1 varies widely from cultivar to cultivar and is correlated with floral induction, indicating that the transcription of FaTFL1 occurs at least partially independently of the FaFT1-FaSOC1 module. Our results indicate that changing the expression patterns of FaTFL1 through biotechnological or conventional breeding approaches could result in strawberries with specific flowering and runnering characteristics including new types of everbearing cultivars.

The Fragaria vesca homolog of suppressor of overexpression of constans1 represses flowering and promotes vegetative growth.

In the annual long-day plant Arabidopsis thaliana, suppressor of overexpression of constans1 (SOC1) integrates endogenous and environmental signals to promote flowering. We analyzed the function and regulation of the SOC1 homolog (Fragaria vesca [Fv] SOC1) in the perennial short-day plant woodland strawberry (Fragaria vesca). We found that Fv SOC1 overexpression represses flower initiation under inductive short days, whereas its silencing causes continuous flowering in both short days and noninductive long days, similar to mutants in the floral repressor Fv terminal flower1 (Fv TFL1). Molecular analysis of these transgenic lines revealed that Fv SOC1 activates Fv TFL1 in the shoot apex, leading to the repression of flowering in strawberry. In parallel, Fv SOC1 regulates the differentiation of axillary buds to runners or axillary leaf rosettes, probably through the activation of gibberellin biosynthetic genes. We also demonstrated that Fv SOC1 is regulated by photoperiod and Fv flowering locus T1, suggesting that it plays a central role in the photoperiodic control of both generative and vegetative growth in strawberry. In conclusion, we propose that Fv SOC1 is a signaling hub that regulates yearly cycles of vegetative and generative development through separate genetic pathways.

Mutation in TERMINAL FLOWER1 reverses the photoperiodic requirement for flowering in the wild strawberry Fragaria vesca.

Photoperiodic flowering has been extensively studied in the annual short-day and long-day plants rice (Oryza sativa) and Arabidopsis (Arabidopsis thaliana), whereas less is known about the control of flowering in perennials. In the perennial wild strawberry, Fragaria vesca (Rosaceae), short-day and perpetual flowering long-day accessions occur. Genetic analyses showed that differences in their flowering responses are caused by a single gene, SEASONAL FLOWERING LOCUS, which may encode the F. vesca homolog of TERMINAL FLOWER1 (FvTFL1). We show through high-resolution mapping and transgenic approaches that FvTFL1 is the basis of this change in flowering behavior and demonstrate that FvTFL1 acts as a photoperiodically regulated repressor. In short-day F. vesca, long photoperiods activate FvTFL1 mRNA expression and short days suppress it, promoting flower induction. These seasonal cycles in FvTFL1 mRNA level confer seasonal cycling of vegetative and reproductive development. Mutations in FvTFL1 prevent long-day suppression of flowering, and the early flowering that then occurs under long days is dependent on the F. vesca homolog of FLOWERING LOCUS T. This photoperiodic response mechanism differs from those described in model annual plants. We suggest that this mechanism controls flowering within the perennial growth cycle in F. vesca and demonstrate that a change in a single gene reverses the photoperiodic requirements for flowering.

The regulation of seasonal flowering in the Rosaceae.

Molecular mechanisms regulating the flowering process have been extensively studied in model annual plants; in perennials, however, understanding of the molecular mechanisms controlling flowering has just started to emerge. Here we review the current state of flowering research in perennial plants of the rose family (Rosaceae), which is one of the most economically important families of horticultural plants. Strawberry (Fragaria spp.), raspberry (Rubus spp.), rose (Rosa spp.), and apple (Malus spp.) are used to illustrate how photoperiod and temperature control seasonal flowering in rosaceous crops. We highlight recent molecular studies which have revealed homologues of terminal flower1 (TFL1) to be major regulators of both the juvenile to adult, and the vegetative to reproductive transitions in various rosaceous species. Additionally, recent advances in understanding of the regulation of TFL1 are discussed.

Major transcriptomic differences are induced by warmer temperature conditions experienced during asexual and sexual reproduction in Fragaria vesca ecotypes.

A major challenge for plants in a rapidly changing climate is to adapt to rising temperatures. Some plants adapt to temperature conditions by generating an epigenetic memory that can be transmitted both meiotically and mitotically. Such epigenetic memories may increase phenotypic variation to global warming and provide time for adaptation to occur through classical genetic selection. The goal of this study was to understand how warmer temperature conditions experienced during sexual and asexual reproduction affect the transcriptomes of different strawberry (Fragaria vesca) ecotypes. We let four European F. vesca ecotypes reproduce at two contrasting temperatures (18 and 28°C), either asexually through stolon formation for several generations, or sexually by seeds (achenes). We then analyzed the transcriptome of unfolding leaves, with emphasis on differential expression of genes belonging to the epigenetic machinery. For asexually reproduced plants we found a general transcriptomic response to temperature conditions but for sexually reproduced plants we found less significant responses. We predicted several splicing isoforms for important genes (e.g. a SOC1, LHY, and SVP homolog), and found significantly more differentially presented splicing event variants following asexual vs. sexual reproduction. This difference could be due to the stochastic character of recombination during meiosis or to differential creation or erasure of epigenetic marks during embryogenesis and seed development. Strikingly, very few differentially expressed genes were shared between ecotypes, perhaps because ecotypes differ greatly both genetically and epigenetically. Genes related to the epigenetic machinery were predominantly upregulated at 28°C during asexual reproduction but downregulated after sexual reproduction, indicating that temperature-induced change affects the epigenetic machinery differently during the two types of reproduction.