Additive QTLs on three chromosomes control flowering time in woodland strawberry (Fragaria vesca L.)

Flowering time is an important trait that affects survival, reproduction and yield in both wild and cultivated plants. Therefore, many studies have focused on the identification of flowering time quantitative trait locus (QTLs) in different crops, and molecular control of this trait has been extensively investigated in model species. Here we report the mapping of QTLs for flowering time and vegetative traits in a large woodland strawberry mapping population that was phenotyped both under field conditions and in a greenhouse after flower induction in the field. The greenhouse experiment revealed additive QTLs in three linkage groups (LG), two on both LG4 and LG7, and one on LG6 that explain about half of the flowering time variance in the population. Three of the QTLs were newly identified in this study, and one co-localized with the previously characterized FvTFL1 gene. An additional strong QTL corresponding to previously mapped PFRU was detected in both field and greenhouse experiments indicating that gene(s) in this locus can control the timing of flowering in different environments in addition to the duration of flowering and axillary bud differentiation to runners and branch crowns. Several putative flowering time genes were identified in these QTL regions that await functional validation. Our results indicate that a few major QTLs may control flowering time and axillary bud differentiation in strawberries. We suggest that the identification of causal genes in the diploid strawberry may enable fine tuning of flowering time and vegetative growth in the closely related octoploid cultivated strawberry.


INTRODUCTION
Synchronizing floral development with local climatic conditions is important for crop and wild plants because of its significant impact on reproduction, yield and ultimately survival. The genetic control of flowering time is thus of immense significance. As a result, numerous mapping studies have focused on the identification of flowering time QTLs in different crops, and a detailed molecular understanding of flowering time regulation in model plants has facilitated the identification of causal genes. [1][2][3][4] In the cultivated strawberry, one of the most economically important berry crops throughout the world, flowering time is a major breeding target because of year-round demand for fresh berries. To extend the cropping season, breeders are developing early and late flowering cultivars as well as everbearing (EB; also called dayneutral) cultivars, which display a continuous flowering habit. [5][6][7][8] In seasonal flowering 'Junebearing' strawberries, flower induction occurs in autumn, and the timing of the induction correlates with flowering time in the following season. 9 Both short days (SD) and cool temperatures can induce flowering. Typically, flower induction occurs independently of day-length at temperatures below~13°C, whereas SDs are required at intermediate temperatures, and at temperatures above 20-24°C floral induction is inhibited. [10][11][12][13] There is however significant variation in critical temperature limits between cultivars, 10,14 and some cultivars possess an obligatory SD requirement for flower induction even under cool temperatures. 15 In contrast to seasonal flowering strawberries, EB cultivars flower earlier under long days (LD) than SDs. 16,17 Both seasonal flowering and EB habits are also found in the woodland strawberry. 18 As the cultivated strawberry is a complex allo-octoploid (2n = 8x = 56), the woodland strawberry (F. vesca; 2n = 2x = 14) and other diploid species have been used as a genetically facile surrogate system for genetic investigation. Several genetic linkage maps have been developed for diploid strawberry species. Initially saturated maps were developed using large numbers of microsatellite (SSR) and other PCR-based markers, [19][20][21][22][23][24][25] and the Fragaria SSR-based reference linkage map was used to anchor the F. vesca 'Hawaii-4' genome sequence. 26 More recently, high-throughput genotyping has been employed using various platforms to produce dense SNP-based linkage maps that have enabled more precise sequence scaffold anchoring and orientation, 27 whereas the most comprehensive linkage map for diploid Fragaria to date is that of F. iinumae produced using the Fragaria Axiom SNP genotyping array (Axiom IStraw90, Affymetrix, Santa Clara, CA, USA), supplemented with SNPs scored using genotyping-bysequencing (GBS). 28,29 The linkage maps have been used to map a number of genetic loci in diploid strawberry species, including fruit color, runnering, leaf color, and seasonal flowering, 20,[30][31][32] and represent a powerful tool for the genetic dissection of traits of economic importance in the cultivated strawberry.
Previous genetic mapping studies have shown that the gene causing seasonal flowering habit in woodland strawberry is located on linkage group 6 (LG6). 21,31 This gene has been recently shown to encode a strong repressor of flowering, TERMINAL FLOWER 1 (FvTFL1), and a 2-bp deletion in the first exon of FvTFL1 causes expression of the EB phenotype in woodland strawberry. 32 Furthermore, homologs of FvTFL1 have been found to repress flowering in several rosaceous species suggesting the conservation of the genetic mechanism of flowering time regulation. [33][34][35][36] Only a few other flowering time genes have been functionally characterized in woodland strawberry. Using genetic transformation, recent studies have revealed how homologs of FLOWERING LOCUS T (FvFT1) and SUPPRESSOR OF THE OVEREXPRESSION OF CONSTANS1 (FvSOC1) mediate photoperiod and temperature signals and how they are integrated to control flower induction. 32,[37][38][39] In the seasonal flowering woodland strawberry, cool temperatures below 13°C cause photoperiod-independent downregulation of FvTFL1 by an unknown mechanism leading to flower induction, whereas at higher temperatures, photoperiod controls flowering through FvTFL1. Under LD, leaf-expressed FvFT1 activates FvSOC1 in the shoot apex leading to the upregulation of FvTFL1 and the maintenance of a vegetative stage, whereas under SD, this FvFT1-FvSOC1-FvTFL1 pathway is silenced and flower induction occurs. However, at temperatures above 20°C FvTFL1 is highly upregulated by an unknown activator that functions independently of FvSOC1, and plants remain vegetative. 37,39 This model contrasts with the mechanism in Arabidopsis in which both FT and SOC1 function as floral activators. 1 However, FvFT1 is a strong floral activator in EB woodland strawberries that are lacking functional FvTFL1. 32,38 The genetic control of the EB habit has also been studied in the cultivated strawberry. Weebadde et al. 5 suggested the presence of several QTLs with significant genotype × environment interaction, 5 whereas other studies indicated that a single dominant locus causes continuous flowering. 40,41 The presence of a single QTL was supported by recent genetic mapping studies by two groups, 6,8 who found a major dominant QTL on LG4 in two independent crossing populations. 6,8 Moreover, Honjo et al. 42 demonstrated that the EB trait is controlled by the same gene in the EB parents of these populations. 42 In addition to the EB trait, this locus called PFRU also controls the production of runners, which are long shoots that enable the efficient clonal reproduction of the species. 6 Although studies indicate that different genes cause the EB phenotype in the woodland and cultivated strawberries, a recent study has shown that TFL1 homologs encode major floral repressors in both species, and that the silencing of this gene also causes the EB habit in cultivated strawberry. 43 Moreover, the regulation of FaTFL1 correlates with flower induction in different seasonal flowering cultivars. 43,44 As strawberries form terminal inflorescences from the apical meristems of the crowns, the number of inflorescences, and consequently the yield potential, depends on the number of side shoots called branch crowns that are produced. [45][46][47][48] SD conditions promote the differentiation of strawberry axillary buds to branch crowns, whereas under LDs axillary buds typically differentiate to runners. 10,14,45 Changes in gibberellin biosynthesis and signaling are known to mediate the photoperiodic differentiation of axillary buds. 49 This differentiation is not understood at the molecular level, although runnerless mutants are known in woodland strawberry. 18 Here, we aimed to identify novel genetic loci controlling flowering time and axillary bud differentiation using the large woodland strawberry F2 mapping population produced by Koskela et al. 32 We demonstrate that five QTLs with additive effects explain about half of the flowering time variance found in the population. One of the QTLs co-localized with the previously identified FvTFL1 and another with PFRU, whereas the others were newly identified in this study. In addition, we report the mapping of two QTLs that affect the number of branch crowns and runners in our mapping population.

Mapping population and phenotyping
Selected lines from a F2 population derived from the cross F. vesca f. semperflorens 'Hawaii-4' (H4) × F. vesca subsp. vesca (FV), (denoted H4 × FV), 32  DNA extraction and quantification DNA was extracted from young leaves of the parents 'H4' and 'FV', along with the F1 and 335 F2 seedlings of the progeny, using a modified version of the CTAB method. 31 The DNA was quantified using a Nanodrop spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA) and was diluted 1:100 (5-10 ng μl − 1 ) for use in PCR. For genome resequencing an additional purification step was carried out. NaCl was added at a final concentration of 0.2 M, and DNA was precipitated by adding 2/3 volumes of cold isopropanol. After centrifugation, the pellet was washed with 100 μl of 70% ethanol and re-suspended in 50 μl of TE.
Amplification and scoring of genetic markers A total of 73 published SSR primer pairs (Supplementary File 1) labeled on the forward primer with either 6-FAM or HEX fluorescent dyes [19][20][21][22]25,[50][51][52] were screened for polymorphism in the parental and F1 lines using the 'Type-it' PCR mastermix (Qiagen, Hilden, Germany), following the PCR protocol reported in Sargent et al. 53 In addition, 12 SSR primer pairs were tested using the fluorescent labeling method described by Schuelke 24 (Supplementary File 1). Samples were diluted 1:75, separated by capillary electrophoresis (Foster City, CA, USA) on a genetic analyzer (ABI 3100, ABI 3130 XL or ABI 3700, Applied Biosystems), and the resultant data were collected and analyzed using GeneMapper or Peak Scanner 2 (Applied Biosystems).
SNPs were identified in the parental genome sequences and the temperature switch PCR (TSP) protocol 23 was used to develop assays for 49 SNP markers within flowering time QTL regions identified in this investigation. In brief, two sets of primers, a locus specific (LS) set and a nested-locus specific (NLS) set were designed per SNP using Primer3 v.0. 40 (ref. 54) and NetPrimer (http://www.premierbiosoft.com/netprimer/netprlaunch/netprlaunch.html), respectively. The LS primers were designed to have a Tm between 60-65°C (optimal 63°C) and a product size 4400 bp, whereas the NLS primers were designed to have a T m of 43-47°C (optimal 45°C) and were extended using a non-complementary 5′ tail region that increased the overall primer T m to 52-55°C (optimal 53°C). 23 PCR was performed as described previously. 23 The primer concentrations and the numbers of PCR cycles were adjusted for each primer individually to optimize marker amplification for ease of scoring. Products were visualized Mapping of flowering time QTLs in strawberry S Samad et al.
over UV light following electrophoresis through a 1.5% (w/v) TAE agarose gel containing gel red (Applied Bioprobes, Rockville, MD, USA) at a constant voltage of 100 V for 2-3 h, depending on the sizes of the fragments produced. The TSP loci were named UH, followed by the chromosome number and finally a unique primer pair reference number (Supplementary File 1).
Genotyping-by-sequencing DNA samples (30-100 ng μl − 1 ) of 188 seedlings selected to represent the full seedling population genotypically (based on the SSR and TSP data) and phenotypically, were sent to the Cornell University where genotyping-bysequencing (GBS) was carried out according to Elshire et al. 55 GBS libraries were sequenced using the Illumina HiSeq2500 platform (San Diego, CA, USA). GBS tags with a minimum number of three reads and unique alignment to the reference genome (https://www.rosaceae.org/species/fragaria_ vesca/genome_v2.0.a1) 27 were identified, and from them SNPs were called by Genomic Diversity Facility in Cornell using the GBS-pipeline 56 implemented in the TASSEL 3.0 software application. 57 To remove genotyping errors, several filtering steps were performed for each polymorphic site using vcftools v. 0.1.12b. 58 A site was accepted if genotype calling rate was over 90% per site, the minor allele frequency was 40.1, and the proportion of heterozygote genotypes fit to the Hardy-Weinberg expectations (2pq, Po 0.01). Missing data was imputed with BEAGLE4 (ref. 59) and the TASSEL5 was used to produce ABH file with 4672 SNPs. Following an initial round of linkage mapping (described below), a second filtering and manual imputation step was performed using the rationale previously reported by Ward et al. 60 This second imputation step removed suspected false homozygous calls that created multiple double-recombination events, inflating mapping distances and confounding local marker ordering.

Heritability, linkage mapping and QTL analysis
The broad sense heritability for flowering time and numbers of runners and branch crowns were calculated in a population of 188 individuals according the formula H 2 = σ 2 g/(σ 2 g+σ 2 e), where σ 2 g = (MS(genotype) − MS (environment))/r and σ 2 e = MS(environment). MS(genotype) and MS(environment) were obtained from analysis of variance for completely randomized design as the mean sums of squares for genotype and residual error, respectively, and r was the effective number of replicates (three to four per genotype). Estimation of these variance components were performed by REML as implemented in PROC VARCOMP procedure in SAS/STAT software, version 9.4 (Copyright 2013, SAS Institute, Cary, NC, USA). Linkage mapping was performed using JOINMAP 4.1 (Kyazma, Wageningen, The Netherlands) 61 using maximum likelihood and the default mapping criteria. The maps presented were illustrated using MapChart 2.2. 62 The significance of association between individual markers on the H4 × FV linkage map and the phenotypes was calculated with linear regression using MapQTL 6.0 (Kyazma). 63 QTLs were identified employing interval mapping with a step size of 1 cM, considering a maximum of five neighboring markers. The genome-wide log of odds (LOD) threshold was determined over 10 000 permutations and the most significant markers were then used as co-factors for restricted multiple QTL mapping (rMQM) with a step size of 1 cM. Linkage mapping and QTL analyses were performed on 335 individuals using SSR and TSP markers and on 186 individuals using all marked data.
Genome resequencing and the identification of single-nucleotide polymorphisms between the 'H4' and 'FV' genomes  58 and only sites having at least 10-fold coverage and o50-fold coverage were accepted for further analysis. Genomic regions spanning flowering time QTL were identified from the genomic positions of flanking SSR markers mapped in the full mapping progeny, within which polymorphic SNPs were identified for subsequent marker development.
Identification of candidate genes in the QTL regions To identify candidate genes for the QTLs, predicted genes located between the markers flanking the most significant markers were BLAST searched against the RefSeq or UniRef90 protein database using the blastx algorithm in Blast2GO 68 or blast+ (2.4.0). Query sequences longer than 8000 bp were restricted to the first 8000 bp for Blast2GO as this is the maximum length accepted by the software. Homologs of known flowering time genes were selected as candidate genes. Parental genetic variation in the coding regions of candidate genes were retrieved in the Illumina whole genome resequencing data using vcftools, 58 and translated amino acid sequences were compared to reveal non-synonymous variation in the coding regions.  Figure 2). Variation was also observed in the vegetative development of the F2 lines. By the end of the experiment, most F2 lines produced 2-4 branch crowns per plant. However, a significant number of genotypes produced an average of less than two, or more than four branch crowns (Figure 1b). The total number of runners produced per genotype varied between 5 and 22 (Figure 1c). The broad sense heritability (H 2 ) was 0.48 and 0.41 for branch crowns and runners, respectively. The number of runners was positively correlated with the number of branch crowns (r = 0.785; Figure 1d). However, no correlation was found between flowering time and the number of branch crowns or runners. 1 755 780 tags (73.4%) were uniquely aligned on the H4 reference genome. Initial SNP calling revealed 133 930 SNPs, but after Hardy-Weinberg filtering and removal of SNPs that had more than 10% missing data and a minor allele frequency below 0.1, an ABH file of 4672 SNPs was retained. Following initial mapping, imputation and filtering, 2395 loci in 186 F2 lines, comprising 53 SSR/TSP markers and 2342 GBS-derived SNP markers coalesced into a linkage map of the seven expected linkage groups that covered a total genetic distance of 558 cM (Figure 2). Plotting genetic distances against the physical positions of the markers on the Fvb genome sequence assembly demonstrated that the genetic map was predominantly co-linear with the genome sequence ( Figure 2) and thus marker placement on the map could be considered generally reliable.

Segregation
Flowering time QTLs in greenhouse and field experiments QTL analysis was first carried out using greenhouse flowering time data on the progeny of 335 seedlings and segregation data for 60 genetic loci. This analysis revealed four significant QTLs, one on both LG6 and LG7, and two on LG4 (Supplementary Figure 4). QTL analysis was also performed in growth chamber flowering time data using 32 F2 lines that were selected based on their extreme flowering time phenotypes in the greenhouse experiment. Using these lines, the same four QTLs with a LOD value above 3 were detected in both greenhouse and growth chamber flowering time data (Table 1). To increase the mapping resolution, we performed QTL analysis using the additional GBS based markers on selected 186 seedlings including individuals that exhibited informative recombination close to the QTL regions mapped in a larger progeny. This high-resolution mapping identified one additional QTL peak on LG7 and narrowed down the other QTL peaks ( Figure 3). The peak of the LOD of the QTL on LG6 corresponded to the position of the FvTFL1 gene previously identified by Koskela et al. 32 (Figure 3). 32 This QTL explained 17% of the observed variance; plants homozygous for 'FV' alleles of FvTFL1 flowered later than heterozygous lines ( Table 2). The peaks of the LOD for two significant loci on LG4 were associated with markers SFvb4_7853414 and bx083, both explaining about 15% of the observed variance. In both cases, the allele from the 'FV' parent promoted flowering. The peaks of the LOD of the significant QTLs on LG7 were associated with markers SFvb7_16204472 and SFvb7_21710529 that explained over 12% of the observed variance each, with the allele of the 'FV' parent promoting earlier flowering. The analysis of different haplotype combinations indicated that QTLs on all three LGs had additive effects on flowering time (Figure 4; Table 2). Plants with 'FV' markers in both LG4 and LG7 flowered first, whereas the latest plants typically had one or two 'H4' alleles in both LGs. In addition, lines containing two functional FvTFL1 alleles tended to flower later than heterozygous lines, regardless of LG4 and LG7 marker genotypes.
The analysis of field phenotypes for flowering time revealed highly significant QTL only on LG4 ( Figure 5). The peak of the LOD of this QTL was associated with marker bx083 as in the greenhouse dataset, but there was also another equally high peak~1 Mb upstream of bx083. This QTL explained almost 25% of the phenotypic variance. Another minor QTL just above the genome-wide significance threshold of 3 was detected on LG2 with a peak of the LOD associated with the marker CFVCT020 at 11.94 cM. However, no QTLs were found on LG6 or LG7 in the field using phenotypic data of 186 F2 lines.
Runnering and branch crown evaluation QTL were revealed on LG4 and LG5 associated with branch crown formation and runnering. The peak of the LOD of the significant QTL for the number of branch crowns was associated with the marker SFvb4_29399865 at 53.184 cM on LG4 and with the marker SFvb_25913832 at 65.718 cM on LG5 (Figures 6a and b), and QTLs for the number of runners were located in the same regions (Supplementary Figure 5). Both QTLs explained over 20% of the variance in the number of branch crowns, but their effect on runnering was less pronounced (Table 2). Increased numbers of runners and branch crowns were associated with 'FV' alleles on both LG4 and LG5.
Candidate genes for QTLs To identify candidate genes around the QTLs, BLAST searches were carried out on the genomic regions between the flanking markers. In most cases, no obvious candidate genes were detected inside the narrow QTL peaks, but promising candidates were found very close to the highest LOD values. FvFT2 (mrna04680.1-v1.0-hybrid) that was also identified in another QTL mapping study, 69 was found within the~400 kb QTL region in the end of the LG4. However, no candidate genes were searched for the other QTL on the LG4 because the physical position of the QTL was uncertain. On LG7, a homolog of EARLY FLOWERING 6 (ELF6, gene23255-v1.0-hybrid) 70 was located about 600 kb upstream of the first QTL peak, and a gene encoding a homolog of floral repressor TFL1 (FvCENTRORADIALIS1, FvCEN1; gene13304-v1.0-hybrid) 32 was found~650 kb upstream of the second QTL peak. Candidate genes were also identified within the QTL intervals controlling branch crown and runner formation. On LG4, a TCP transcription factor (FvTCP7, gene04759-v1.0-hybrid) 71 was found 150 kb upstream of the QTL peak. Furthermore, genes encoding two closely related MADS transcription factors that are homologous to DORMANCY ASSOCIATED MADS BOX (DAM, mrna12119.1-v1.0-hybrid, mrna12120.1-v1.0-hybrid), 2,72,73 were found inside the QTL region on the LG5. No non-synonymous variation was found in the coding sequence of any candidate gene between the parental genome sequences (data not shown).

DISCUSSION
Extension of the strawberry production season in the open field is of significant economic importance, and two different breeding strategies have been implemented to maximize the length of the season. Several breeding programs are focusing on EB cultivars that can produce berries for several months during a single season, 6,8,42 whereas the breeding of early and late ripening cultivars is an alternative approach to extend the ripening season. 7 Previous QTL mapping studies in strawberry have focused on the perpetual flowering habit, 5,6,8,42 and no QTLs have been reported for earliness. Here, we report the identification of additive QTLs on three LGs that explain over 50% of the observed 18-day variance in flowering time in woodland strawberry mapping population previously reported by Koskela et al. 32 Moreover, we have found two QTLs that affect branch crown formation, the trait that has a direct effect on strawberry yield potential. 45,47,48 Additive QTLs on three linkage groups control flowering time in woodland strawberry QTL mapping using replicated greenhouse flowering time data in our H4 × FV F2 population revealed one new QTL on the LG4, two on the LG7, and two additional QTLs that co-localized with previously mapped QTLs on the LG4 and LG6. Four of these QTLs were also detected in a replicated growth chamber experiment using 32 F2 selected lines that exhibited extreme flowering time phenotypes in a greenhouse experiment, indicating that these QTLs are robust at least in controlled climate. In the non-replicated  field experiment, the analysis using GBS data on 186 lines revealed major QTL only on LG4. However, QTL mapping on a larger population of 335 lines using SSR markers showed LOD values close to 4.0 in both LG6 and LG7 (data not shown), but a replicated field experiment is needed to confirm these QTLs in the field. On the LG6, the QTL mapped on the same region than a previously characterized FvTFL1 gene that encodes a repressor of flowering. 32 In cultivated strawberry, a major QTL for perpetual flowering and runnering (PFRU) was previously identified on the LG4. 6,8,42 A recent fine-mapping study placed this QTL close to marker bx083, 69 and the same marker was located in the QTL peak also in our greenhouse and field experiments. Therefore, our data suggest a new role for PFRU in the control of earliness, in addition to its role in the flowering period, at least in woodland strawberry. The effect of PFRU on flowering time seems to be robust in different environments. However, the other QTL, that is located upstream of PFRU on LG4, was only detected in controlled climate experiments indicating that it may function only in specific environments.
We found other QTLs on LG7 where no flowering related QTLs have been previously detected in strawberries. These QTLs showed high LOD values in the greenhouse experiment with the F2 population of 186 lines (GBS data), but not in the field indicating that they were also influenced by the environment. One of these QTLs matched also with the QTL region detected in our growth chamber experiment (Table 1). This QTL was associated with the marker BFaCT044 at 22 166 857 bp on LG7 in a set of 32 early and late flowering lines in both growth chamber and greenhouse, and about 400 kb upstream (the most significant marker SFvb7_21710529) in the GBS dataset of 186 lines.
The analysis of allele combinations on QTL positions demonstrated their additive effects on flowering time ( Figure 4; Table 2). Clearly, 'H4' specific markers on LG4 and LG7 QTL regions, in both heterozygous and homozygous condition delayed flowering. Furthermore, homozygous FvTFL1 alleles from 'FV' parent caused additional delay compared to heterozygous plants. Plants homozygous for 'H4' alleles of FvTFL1 were excluded from this study because they exhibit perpetual flowering habit. 32 Candidate genes for flowering time QTLs We searched for candidate flowering time genes in the genomic regions containing the markers with the highest LOD if the physical and genetic maps on the region were mostly co-linear ( Figure 2). This was not the case for the newly identified QTL on the LG4, and therefore, a better genome assembly is needed to enable reliable candidate gene search on this region.
Candidate genes for PFRU, that co-localize with the QTL identified in our greenhouse and field studies on the LG4, has been recently reported. 69 These genes include FvFT2 and FvCDF2 that are homologs of known regulators of photoperiodic flowering. In Arabidopsis, CDFs are floral repressors that contribute to photoperiodic flowering by down-regulating both CO and FT mRNA expression especially in the morning. 74 Furthermore, specific alleles of potato CDF gene were selected during domestication allowing the cultivation of this SD crop in LD conditions in northern latitudes. 75 Tissue specific gene expression analyses in woodland strawberry showed that FvFT2 is mainly expressed in flowers and fruits. 32,76 However, both FvFT2 and FvCDF2 may also have important roles in the control of photoperiodic development in strawberries and are thus good candidates for future functional studies.
FvELF6 was identified as a candidate gene for the QTL on LG7. ELF6 was originally identified as a repressor of photoperiodic flowering in Arabidopsis. 70 It is a histone demethylase that control gene expression through chromatin regulation. 77,78 As histone methylation is a common mechanism in the control of flowering time, 79 and a homologous gene to ELF6 co-localizes with a QTL for spring heat requirement for flowering on almond LG2, 2 thus, ELF6 is an interesting candidate for future studies.
The homolog of TFL1/CEN is located in the QTL region for flowering time in peach 72,80 and gene functional studies showed that TFL1 homologs encode major floral repressors in several rosaceous species. [32][33][34][35] Interestingly, TFL1/CEN homologs were identified as candidate genes for two QTL regions in this study. A QTL on the LG6 co-localize with FvTFL1 that encodes a major floral repressor, 32 and closely related FvCEN1 is located near the second QTL peak on LG7 and awaits for functional validation. Our finding that plants containing only one functional FvTFL1 copy flowered   earlier than plants homozygous for functional 'FV' alleles provides interesting opportunities to control flowering time, particularly in polyploid crops such as the octoploid cultivated strawberry. In fact, a recent report showed that FaTFL1 is a major floral repressor also in cultivated strawberry, and that altered regulation of FaTFL1 mRNA expression is associated with different flowering responses in this species. 43 To enable more efficient breeding of flowering time, further studies are needed to uncover whether the variation in FaTFL1 expression between cultivars is caused by allelic variation of FaTFL1 itself or by upstream regulators. Among these regulators, the expression of FT1 and SOC1 homologs correlate with TFL1 mRNA levels only in specific conditions in woodland and cultivated strawberries suggesting that TFL1 alleles or its other regulators should be major targets of future investigations. 39,43 Whether some of the candidate genes identified in this study control flowering time through FvTFL1 is an interesting open question.
Identification of QTLs controlling crown branching and runner formation Although we mainly focused on flowering time, we also identified two QTLs on the LG4 and LG5 controlling the number of branch crowns and runners in woodland strawberry. These traits showed a strong positive correlation because the production of branch crowns increases the number of axillary buds that are able to form runners under LD conditions. Therefore, these QTL primarily control the branching of the leaf rosette in woodland strawberry. On the LG4, the QTL co-localized with previously mapped PFRU locus. 6,69 However, the peak of the LOD was located about 700 kb downstream of the flowering time QTL mapped in this study. This indicates that two closely located QTLs may control flowering and axillary bud differentiation to runners and branch crowns in our population, but additional studies are needed to confirm this hypothesis.
We identified a gene encoding TCP transcription factor close to the most significant marker on the LG4. This TCP, recently named as FvTCP7, 71 shows the highest similarity with Arabidopsis TCP14 and TCP15, which control cell proliferation and internode elongation. 81 As runners are long shoots with elongated internodes this is a promising candidate gene for the strawberry axillary bud differentiation. Initial characterization of FvTCP7 indicated that it is highly expressed in vegetative tissues including runners, but also in flower buds, and the protein is localized to nucleus. 71 However, the function of FvTCP7 in strawberry is unknown. Further studies are needed to elucidate what is the relationship between the branching QTL detected here and PFRU that controls runner formation, in addition to perpetual flowering, in cultivated strawberry. 6,8 Also the newly identified QTL region on the LG5 contained interesting candidate genes for future studies including two genes homologous to DAM and SHORT VEGETATIVE PHASE (SVP) that encode major regulators of dormancy in Rosaceae and floral repressor in Arabidopsis, respectively. 2,72,73,82 Although strawberries do not have true dormancy, their growth vigor is strongly reduced under SDs in autumn including the cessation of runner formation, and a long period of chilling is needed to resume normal vegetative development. 83 Functional studies would reveal whether identified FvDAM transcription factors control dormancy and what is their role in the control of axillary bud differentiation. Furthermore, the possible interaction of the branching/runnering QTLs identified in this study with the FvSOC1 and gibberellin pathway, that were previously shown to control runnering in woodland strawberry, 37,49 should be elucidated.

CONCLUSIONS
Here, we have reported the mapping of three new flowering time QTLs in woodland strawberry, which co-localize with homologs of known flowering time regulators that have not been functionally characterized in this species. In addition, we have shown that PFRU, another QTL that was previously reported to control the duration of flowering and the number of runners in cultivated strawberry, 6,69 affects flowering time, crown branching and runner formation in woodland strawberry. Although we were able to detect four out of five flowering time QTLs in two independent experiments in controlled climate, only one of these showed a LOD value above 3 in the non-replicated field experiment. Therefore, further replicated experiments are needed to reveal the role of these QTLs in different environments. On the basis of the evidence shown here, identified QTLs control flowering time additively along with a previously identified floral repressor FvTFL1 that has quantitative effect on flowering time. As TFL1 homologs likely integrate both light and temperature signals to control flowering time in woodland and cultivated strawberries, 39,43 we suggest that candidate genes on the LG4 and LG7 QTL regions may control flowering time through TFL1. However, further studies are needed to clarify causal genes and to reveal their roles in the control of flowering time and vegetative traits. As cultivated strawberry is an octoploid species with at least two subgenome donors, 27,84 we hypothesize that it may contain significant genetic variation in these regulators that may enable fine tuning of flowering time and shoot architecture in different environments through marker assisted selection or genomic selection strategies that have recently been developed in this species. 85