Abstract
U.S. watermelon (Citrullus lanatus) production is worth ≈$0.5 billion annually to growers and nearly all of them are dependent on reliable synchronized flowering time of triploid cultivars and diploid pollenizers in their production fields. One aspect of this synchronization is time to flowering, the change from the vegetative to reproductive phase of a plant. Flowering time has emerged as one of the key traits in horticultural and agronomic crops to breed for escape from biotic and abiotic stresses. However, very little is known about the control of flowering time in watermelon. The number of genes involved, mode of inheritance, heritability, and the possible candidate genes are all unknown. In this study, quantitative trait loci (QTL) associated with days to first male flower (DMF), days to first female flower (DFF), and the female-male flower interval (FMI) were identified in a ‘Klondike Black Seeded’ × ‘New Hampshire Midget’ recombinant inbred line population over 2 years. Heritability for DMF, DFF, and FMI were 0.43, 0.23, and 0.10, respectively. Control of flowering time was oligogenic with a major, stable, colocalized QTL on chromosome 3 responsible for ≈50% of the phenotypic variation observed for DMF and DFF. This region of the draft genome sequence contains 172 genes, including homologs of the flowering locus T (Cla009504) and tempranillo 1 (Cla000855) genes associated with flowering time in other species. Cla009504 and Cla000855 represent excellent candidate genes toward the development of a functional marker for marker-assisted selection of flowering time in watermelon. In addition to the major QTL on chromosome 3, two other QTL were identified for DMF (chromosomes 2 and 3) and DFF (chromosomes 3 and 11) and one for FMI on chromosome 2. Understanding the genes involved in this trait and the ability to select efficiently for flowering time phenotypes is expected to accelerate the development of new watermelon cultivars in changing environmental conditions.
The change from the vegetative to the reproductive phase of a plant is an important aspect in most commercial agronomic and horticultural crops. A shorter production cycle often means lower input costs for the producer and less environmental impact through energy-saving and reduced pesticide use. A short production cycle also allows for escape from pathogens by harvesting crops before environmental conditions are favorable for high disease pressure (Poland et al., 2009). Flowering time is thus often used to avoid biotic and abiotic stresses. The time to first male and female flower and the interval between the appearance of male and female flowers in monoecious plants have also emerged as key traits used to breed for increased yield under drought conditions (Bolaños and Edmeades, 1993; Chapman and Edmeades, 1999; Richards, 2006; Siddique et al., 1990). The association of drought quantitative trait loci and flowering time loci (Ducrocq et al., 2008) and the realization that selection for specific flowering phenotypes increases yield under drought conditions have contributed to flowering time traits emerging as key breeding priorities for the future (Jung and Müller, 2009).
Flowering traits are complex and it is anticipated that climate change will greatly influence flowering (Craufurd and Wheeler, 2009). As can be expected for such an important trait, a large amount of research has been carried out and the genetic architecture and molecular pathways for flowering time are well described for model systems such as Arabidopsis thaliana and are becoming clearer for important agronomic crops (for recent reviews, see Buckler et al., 2009; Jackson, 2009; Matsoukas et al., 2012; Turck et al., 2008). The picture that is emerging from A. thaliana shows that the molecular pathways that control flowering include a large number of genes representing photoperiod (light signaling and circadian clock), vernalization, autonomous signals, and gibberellin biosynthesis. Very briefly, in the leaf constans [CO (Valverde et al., 2004)] transcription is up-regulated in response to signals received from the photoperiodic induction pathway. The accumulation of CO protein activates flowering locus T [FT (Kardailsky et al., 1999; Kobayashi et al., 1999; Valverde et al., 2004)], which is a floral inducer. FT can also be induced through a CO- independent pathway where FT suppressors are down-regulated. After moving to the meristem, the FT protein promotes other pathway integrator genes (Lee et al., 2008; Liu et al., 2008; Yoo et al., 2005), which affect floral meristem identity genes, leading to flowering (Abe et al., 2005; Lee et al., 2008; Melzer et al., 2008; Wigge et al., 2005).
Watermelon production is responsible for ≈7% of world vegetable production acreage (Food and Agriculture Organization of the United Nations, 2011). However, very little is known about the control of flowering time in this monoecious crop. Commercial cultivars are day-neutral (George, 2009) and male flowers usually open first, followed by female flowers at a ratio of roughly 7:1 (Wehner, 2008). The interval between the opening of the first male flower and first female flower is dependent on the cultivar and environmental conditions.
The percentage of watermelon production in the United States devoted to seedless fruit has increased dramatically from 51% in 2003 to reach 85% in 2009 (U.S. Department of Agriculture, 2011). The seedless fruit is produced on triploid plants, but because pollination is required for fruit set and these plants produce negligible amounts of viable pollen, diploid pollen sources (pollenizers) are required (Boyhan et al., 2000; Maynard, 1992; Maynard and Elmstrom, 1992). As a result of this local preference for seedless fruit, U.S. watermelon production is dependent on synchronized flowering of diploid pollenizers and triploid watermelon cultivars for fruit production. However, the number of genes involved, mode of inheritance, heritability, and the possible candidate genes are all unknown for this crop. These factors are critical impediments to efficient breeding of watermelon cultivars in a changing environment, and especially for the United States where fruit production is dependent on synchronized flowering between two cultivars (a triploid and a pollenizer).
Recently, molecular tools have become available for watermelon that can now be used to advance breeding efforts. We recently produced the first single nucleotide polymorphism (SNP) maps for watermelon (Sandlin et al., 2012) and mapped the first QTL in the species (Prothro et al., 2012a, 2012b, 2013; Sandlin et al., 2012). In 2012, the draft watermelon genome sequence (Guo et al., 2013) became available and here we describe the use of these resources to elucidate the control of flowering time in watermelon.
Materials and Methods
Plant materials and genetic map.
The ‘Klondike Black Seeded’ (KBS) × ‘New Hampshire Midget’ (NHM) recombinant inbred line (RIL) population developed previously to map QTL associated with fruit and seed traits (Prothro et al., 2012a; Sandlin et al., 2012) was advanced in the greenhouse (Jan. to Apr. 2012) through single seed descent (SSD) from F7 to F8 to generate seed for the current experiment.
Trait phenotyping.
During Summer 2012 and 2013, 150 RILs and parental cultivars were sown in seedling trays and transplanted ≈2 weeks later on 23 May 2012 and 22 May 2013 at the Durham Horticulture Farm, in Watkinsville, GA. One plant per RIL/parental cultivar was planted per block in a completely randomized bock design with 10 blocks/replications. Plants were grown according to University of Georgia Cooperative Extension Service recommendations on plastic mulch with between-row spacing of 1.83 m and in-row spacing of 0.9 m. In 2012, the vines were turned to make data collection easier, but, in 2013, high rainfall limited entry into the field to only days essential for data collection.
Data were collected three times per week for days to the opening of the first male flower, days to the opening of the first female flower, and subsequently calculated for female–male interval. Bisexual flowers were scored as female. Pearson correlations were calculated using JMP 9.0.2 (SAS Institute, Cary, NC).
Heritability on family-mean base was calculated as h2 = σ2G /{[σ2G + (σ2G*Y/e) + (σ2e/er)]} for the combined environments (Holland et al., 2003; Nyquist and Baker, 1991), where σ2G equaled genetic variance among the entries, σ2G*Y the variance of genetic × years interaction, σ2e the variance of experimental error, e the number of years, and r the number of replications. Because the data were unbalanced, e and r were computed as harmonic mean of years and harmonic mean of total replications across all years, respectively (Holland et al., 2003).
QTL detection.
Composite interval mapping was used to identify QTL associate with the three traits using Model 6, a walk speed of 2 cM, and a window size of 5 cM (Zeng, 1994; Zeng et al., 1999) with WinQTL Cartographer Version 2.5 (Wang et al., 2011). Data were analyzed separately for different years as well as jointly (averaged over years). Significance of the QTL was determined using 1000 permutations (α = 0.05) (Churchill and Doerge, 1994; Doerge and Churchill, 1996)
Candidate genes.
The sequences of the SNP markers (Sandlin et al., 2012) closest to the QTL (Table 2) were used to determine the approximate location of the stable QTL on chromosome 3 of the draft genome sequence (Guo et al., 2013). The predicted genes in that region were then compared with genes known be involved in flowering time in other crops (Blackman et al., 2011; Ehrenreich et al., 2009; Jung et al., 2012; Matsoukas et al., 2012).
Results and Discussion
The average DMF for NHM and KBS was 9.4 and 21 d, respectively, in 2012 and 11.7 and 25.8 d in 2013 (Fig. 1A). In 2012, the average DFF for NHM and KBS was 18.3 and 29.1 d and in 2013 it was 26.5 and 34.6 d (Fig. 1B). The FMI for the 2 years was 9.5 and 14.7 d for NHM and 7.1 and 8.5 d for KBS, respectively (Fig. 1C) with all traits showing some degree of transgressive segregation in the RIL population. There was a significant positive correlation between DMF and DFF within and across years and a significant negative correlation between DMF and MFI (Table 1). It appears that a larger MFI is associated with earlier male flowers rather than later female flowers, although there was also a significant negative correlation between MFI2013 and DFF2012 (but not DFF2013).
Pearson correlations for days to first male flower (DMF), days to first female flower (DFF) and male–female interval (FMI) for 2012 and 2013 in the ‘Klondike Black Seeded’ × ‘New Hampshire Midget’ watermelon recombinant inbred line population.
Generally, plants flowered later in 2013 (DMFavg = 19.2 d, DFFavg = 30.3 d) than 2012 (DMFavg = 16.1 d, DFFavg = 23.9 d) (Fig. 1A–B), probably as a result of the unusually high amount of rainfall in 2013. The total rainfall for the duration of this experiment was 98.3 mm in 2012 (10 rain days) and 292.9 mm in 2013 (21 rain days), whereas the average for this period from 2002 to 2011 was 137.1 mm with 14.5 rain days. DFF was more affected than DMF, leading to a higher FMI in 2013 (11.1 d) than 2012 (7.7 d), probably as a result of the timing of the high rainfall in the period from 11 to 27 d after transplanting (69.3 mm in 2012 and 276.1 mm in 2013). Heritability for DMF, DFF, and MFI was 0.43 ± 0.035, 0.23 ± 0.030, and 0.10 ± 0.022, respectively. This is much lower than what has been observed for days to male flowering (h2 = 0.84), days to female flowering (h2 = 0.83), and anthesis-silking interval (h2 = 0.68) in the highly replicated maize (Zea maize) populations described by Buckler et al. (2009). In melon (Cucumis melo), heritability of days to anthesis (DA) of 0.64 (Zalapa et al., 2008) has been reported. However, it should be noted that in the latter study DA was defined as the time from transplanting to the time when 50% of plants were flowering rather than the individual plant basis of the current study.
Three QTL were associated with DMF in 2012 and one in 2013, whereas two QTL were identified for DFF and one for MFI in the 2 years (Table 2; Fig. 2). A major QTL for DMF (Qdmf3-1) and DFF (Qdff3-1) located on chromosome 3 explained ≈50% of the phenotypic variance observed in the population and was stable across the 2 years (Table 2; Fig. 2). In 2013 a QTL for FMI (Qfmi3) was also detected in this region. Additional QTL for DMF (Qdmf3-2) and DFF (Qdff3-2) were identified on chromosome 3 in 2012 but were not colocalized. QTL for DMF (Qdmf2) and MFI (Qmfi2) in 2012 were identified at a similar position on chromosome 2, but were not detected in 2013, whereas a QTL for DFF (Qdff11) was detected on chromosome 11 in 2013 but not in 2012. Joint analysis (averaged over years) detected the colocalized QTL on chromosomes 2 (Qdmf2 and Qfmi2) and 3 (Qdmf3-1, Qdff3-1, and Qfmi3).
Genomic regions associated with quantitative trait loci (QTL) for days to first male flower (DMF), days to first female flower (DFF), and male–female interval (FMI) for 2012 (n = 145) and 2013 (n = 144) in the ‘Klondike Black Seeded’ × ‘New Hampshire Midget’ watermelon recombinant inbred line population.
Co-localization for QTL for flowering traits is common and has been observed in other species (Buckler et al., 2009). QTL mapping in other crops suggested that a small number of genes with large effects are associated with flowering time in self-pollinating crops (Cockram et al., 2007; Izawa et al., 2003; Wills and Burke, 2007), whereas in out-crossing crops like maize, flowering time is controlled by a large number of QTL with small additive effects (Buckler et al., 2009). Despite being an out-crossing crop, our results indicate oligogenic control of flowering time in domesticated watermelon. However, it should be noted that watermelon has high natural self-pollination rates, as reflected by the lack of inbreeding depression in the species (Kumar et al., 2013; Kumar and Wehner, 2010; Wehner, 2008). Major colocalized QTL were also found to be associated with days to male and female flower development in Cucurbita pepo (Esteras et al., 2012). It remains to be seen whether oligogenic control of flowering time holds true for wild relatives of these cultivated crops.
Favorable alleles (early flowering) at the major QTL (Qdmf3-1, Qdff3-1) are contributed by the early flowering parent NHM (Table 2). However, the favorable allele at Qdmf2 is contributed by the late flowering KBS parent. This is not unexpected, because antagonistic additive effects are expected for traits with transgressive segregation (deVicente and Tanksley, 1993; Rieseberg et al., 1999).
The LOD-1.5 support interval [≈95% confidence interval (Silva et al., 2012)] for the major QTL on chromosome 3 stretches from 8.2 to 12 cM, which corresponds to approximately the 12.0 to 17 Mbp region on chromosome 3 of the watermelon draft genome (Guo et al., 2013). This region contains 172 predicted genes, including Cla009504, a homolog of FT in squash [Cucurbita moschata, Cm_FTL2 (Lin et al., 2007)] and rice [Oryza sativa, Hd3a (Kojima et al., 2002)] and Cla000855, an AP2/ERF and B3 domain-containing transcription factor with homology to TEM1 in A. thaliana (Castillejo and Pelaz, 2008). Expression of Cm_FTL2 and Hd3a in A. thaliana lead to earlier flowering (Kojima et al., 2002; Lin et al., 2007), whereas TEM1 is a repressor of FT and thus up-regulation of TEM1 leads to delayed flowering (Castillejo and Pelaz, 2008; Osnato et al., 2012). TEM1 also regulates genes associated with gibberellin biosynthesis and thus connects the latter to the photoperiod pathway (Osnato et al., 2012). Cla009504 and Cla000855 represent excellent candidate genes toward the development of functional markers (Andersen and Lübberstedt, 2003) for marker-assisted selection (MAS) for flowering time in watermelon. We are currently sequencing the Cla009504 and Cla000855 alleles in this population and examining gene expression in early and late flowering genotypes. Flowering time traits are the target for MAS in several economically important crops including maize (Ducrocq et al., 2009), canola [Brassica napus (Raman et al., 2013)], soybean [Glycine max (Zhang et al., 2013)], rice (Yano et al., 2001), and wheat [Triticum aestivum (Yan et al., 2006)]. Functional markers, where the genotypic sequence used for selection is the cause of the phenotype, is the preferred marker type because there is no recombination between the marker and the trait gene (Andersen and Lübberstedt, 2003).
To ensure the broad applicability of any developed markers, future research needs to validate the stability of the QTL in multiple environments. During the RIL population seed increase through SSD in the greenhouse, data were also collected for DMF and DFF. Although this was under artificial light (14 h light/10 h dark) in the greenhouse and each line was represented by only a single plant, the major QTL on chromosome 3 (Qdmf3-1 and Qdff3-1) were associated with flowering time [DMFGH-LOD = 26.1, DMFGH-R2 = 50%, DFFGH-LOD = 14.0, DFFGH-R2 = 28% (data not shown)]. This suggests that this QTL will have broad applicability in MAS for flowering time in watermelon.
Other aspects that need to be addressed in future research are the duration and consistency of flowering. These are critical traits in breeding for synchronized flowering where consistent pollen availability from pollenizers at the time when female triploid flowers are receptive, is essential.
Seedless watermelon production has unique challenges that require innovative solutions that will move the science of watermelon breeding forward. We identified a major, stable QTL associated with flowering time in watermelon, which gives us insight into the genetic architecture of this economically important trait. In addition, we identified candidate genes for further study toward elucidation of flowering time pathways and MAS that will reduce the breeding cycle of watermelon cultivars.
Literature Cited
Abe, M., Kobayashi, Y., Yamamoto, S., Daimon, Y., Yamaguchi, A. & Ikeda, Y. 2005 FD, a bZIP protein mediating signals from the floral pathway integrator FT at the shoot apex Science 309 1052 1056
Andersen, J.R. & Lübberstedt, T. 2003 Functional markers in plants Trends Plant Sci. 8 554 560
Blackman, B.K., Rasmussen, D.A., Strasburg, J.L., Raduski, A.R., Burke, J.M., Knapp, S.J., Michaels, S.D. & Rieseberg, L.H. 2011 Contributions of flowering time genes to sunflower domestication and improvement Genetics 187 271 287
Bolaños, J. & Edmeades, G.O. 1993 Eight cycles of selection for drought tolerance in lowland tropical maize. I. Responses in grain yield, biomass, and radiation utilization Field Crops Res. 31 233 252
Boyhan, G.E., Granberry, D.M. & Kelley, T.W. 2000 Commercial watermelon production. Univ. Georgia Coop. Ext. Serv. Bul. 996
Buckler, E.S., Holland, J.B., Bradbury, P.J., Acharya, C.B., Brown, P.J., Browne, C., Ersoz, E., Flint-Garcia, S., Garcia, A., Glaubitz, J.C., Goodman, M.M., Harjes, C., Guill, K., Kroon, D.E., Larsson, S., Lepak, N.K., Li, H., Mitchell, S.E., Pressoir, G., Peiffer, J.A., Rosas, M.O., Rocheford, T.R., Romay, M.C., Romero, S., Salvo, S., Villeda, H.S., Sofia da Silva, H., Sun, Q., Tian, F., Upadyayula, N., Ware, D., Yates, H., Yu, J., Zhang, Z., Kresovich, S. & McMullen, M.D. 2009 The genetic architecture of maize flowering time Science 325 714 718
Castillejo, C. & Pelaz, S. 2008 The balance between CONSTANS and TEMPRANILLO activities determines FT expression to trigger flowering Curr. Biol. 18 1338 1343
Chapman, S.C. & Edmeades, G.O. 1999 Selection improves drought tolerance in tropical maize populations. II. Direct and correlated changes among secondary traits Crop Sci. 39 1315 1324
Churchill, G.A. & Doerge, R.W. 1994 Empirical threshold values for quantitative trait mapping Genetics 138 963 971
Cockram, J., Jones, H., Leigh, F.J., O'Sullivan, D., Powell, W., Laurie, D.A. & Greenland, A.J. 2007 Control of flowering time in temperate cereals: Genes, domestication, and sustainable productivity J. Expt. Bot. 58 1231 1244
Craufurd, P.Q. & Wheeler, T.R. 2009 Climate change and the flowering time of annual crops J. Expt. Bot. 60 2529 2539
deVicente, M.C. & Tanksley, S.D. 1993 QTL analysis of transgressive segregation in an interspecific tomato cross Genetics 134 585 596
Doerge, R.W. & Churchill, G.A. 1996 Permutation tests for multiple loci affecting a quantitative character Genet. Mol. Biol. 142 285 294
Ducrocq, S., Giauffret, C., Madur, D., Combes, V., Dumas, F., Jouanne, S., Coubriche, D., Jamin, P., Moreau, L. & Charcosset, A. 2009 Fine mapping and haplotype structure analysis of a major flowering time quantitative trait locus on maize chromosome 10 Genetics 183 1555 1563
Ducrocq, S., Madur, D., Veyrieras, J.-B., Camus-Kulandaivelu, L., Kloiber-Maitz, M., Presterl, T., Ouzunova, M., Manicacci, D. & Charcosset, A. 2008 Key impact of Vgt1 on flowering time adaptation in maize: Evidence from association mapping and ecogeographical information Genetics 178 2433 2437
Ehrenreich, I.M., Hanzawa, Y., Chou, L., Roe, J.L., Kover, P.X. & Purugganan, M.D. 2009 Candidate gene association mapping of Arabidopsis flowering time Genetics 183 325 335
Esteras, C., Gomez, P., Monforte, A.J., Blanca, J., Vicente-Dolera, N., Roig, C., Nuez, F. & Pico, B. 2012 High-throughput SNP genotyping in Cucurbita pepo for map construction and quantitative trait loci mapping BMC Genomics 13 80
Food and Agriculture Organization of the United Nations 2011 Crop production 2010. 14 Jan. 2013. <http://faostat.fao.org/site/567/default.aspx#ancor>
George, R.A.T. 2009 Vegetable seed production. CABI, Cambridge, MA
Guo, S., Zhang, J., Sun, H., Salse, J., Lucas, W.J., Zhang, H., Zheng, Y., Mao, L., Ren, Y., Wang, Z., Min, J., Guo, X., Murat, F., Ham, B.-K., Zhang, Z., Gao, S., Huang, M., Xu, Y., Zhong, S., Bombarely, A., Mueller, L.A., Zhao, H., He, H., Zhang, Y., Zhang, Z., Huang, S., Tan, T., Pang, E., Lin, K., Hu, Q., Kuang, H., Ni, P., Wang, B., Liu, J., Kou, Q., Hou, W., Zou, X., Jiang, J., Gong, G., Klee, K., Schoof, H., Huang, Y., Hu, X., Dong, S., Liang, D., Wang, J., Wu, K., Xia, Y., Zhao, X., Zheng, Z., Xing, M., Liang, X., Huang, B., Lv, T., Wang, J., Yin, Y., Yi, H., Li, R., Wu, M., Levi, A., Zhang, X., Giovannoni, J.J., Wang, J., Li, Y., Fei, Z. & Xu, Y. 2013 The draft genome of watermelon (Citrullus lanatus) and resequencing of 20 diverse accessions Nat. Genet. 45 51 58
Holland, J.B., Nyquist, W.E. & Cervantes-Martinez, C.T. 2003 Estimating and interpreting heritability for plant breeding: An update Plant Breed. Rev. 22 9 111
Izawa, T., Takahashi, Y. & Yano, M. 2003 Comparative biology comes into bloom: Genomic and genetic comparison of flowering pathways in rice and Arabidopsis Curr. Opin. Plant Biol. 6 113 120
Jackson, S.D. 2009 Plant responses to photoperiod New Phytol. 181 517 531
Jung, C. & Müller, A.E. 2009 Flowering time control and applications in plant breeding Trends Plant Sci. 14 563 573
Jung, C.-H., Wong, C.E., Singh, M.B. & Bhalla, P.L. 2012 Comparative genomic analysis of soybean flowering genes PLoS One 7 e38250
Kardailsky, I., Shukla, V.K., Ahn, J.H., Dagenais, N., Christensen, S.K., Nguyen, J.T., Chory, J., Harrison, M.J. & Weigel, D. 1999 Activation tagging of the floral inducer FT Science 286 1962 1965
Kobayashi, Y., Kaya, H., Goto, K., Iwabuchi, M. & Araki, T. 1999 A pair of related genes with antagonistic roles in mediating flowering signals Science 286 1960 1962
Kojima, S., Takahashi, Y., Kobayashi, Y., Monna, L., Sasaki, T., Araki, T. & Yano, M. 2002 Hd3a, a rice ortholog of the Arabidopsis FT gene, promotes transition to flowering downstream of Hd1 under short-day conditions Plant Cell Physiol. 43 1096 1105
Kumar, R., Dia, M. & Wehner, T.C. 2013 Implications of mating behavior in watermelon breeding HortScience 48 960 964
Kumar, R. & Wehner, T. 2010 Natural outcrossing in watermelon—A review Cucurbit Genet. Coop. Rpt. 33/34 42 43
Lee, J., Oh, M., Park, H. & Lee, I. 2008 SOC1 translocated to the nucleus by interaction with AGL24 directly regulates LEAFY Plant J. 55 832 843
Lin, M.-K., Belanger, H., Lee, Y.-J., Varkonyi-Gasic, E., Taoka, K.-I., Miura, E., Xoconostle-Cázares, B., Gendler, K., Jorgensen, R.A., Phinney, B., Lough, T.J. & Lucas, W.J. 2007 FLOWERING LOCUS T protein may act as the long-distance florigenic signal in the cucurbits Plant Cell 19 1488 1506
Liu, C., Chen, H., Er, H.L., Soo, H.M., Kumar, P.P. & Han, J.H.E.A. 2008 Direct interaction of AGL24 and SOC1 integrates flowering signals in Arabidopsis Dev. Change 135 1481 1491
Matsoukas, I.G., Massiah, A.J. & Thomas, B. 2012 Florigenic and antiflorigenic signaling in plants Plant Cell Physiol. 53 1827 1842
Maynard, D.N. 1992 Growing seedless watermelon. Univ. Florida. Coop. Ext. Serv. Fact Sheet HS687
Maynard, D.N. & Elmstrom, G.W. 1992 Triploid watermelon production practices and cultivars Acta Hort. 318 169 173
Melzer, S., Lens, F., Gennen, J., Vanneste, S., Rohde, A. & Beeckman, T. 2008 Flowering-time genes modulate meristem determinacy and growth form in Arabidopsis thaliana Nat. Genet. 40 1489 1492
Nyquist, W.E. & Baker, R. 1991 Estimation of heritability and prediction of selection response in plant populations Crit. Rev. Plant Sci. 10 235 322
Osnato, M., Castillejo, C., Matías-Hernández, L. & Pelaz, S. 2012 TEMPRANILLO genes link photoperiod and gibberellin pathways to control flowering in Arabidopsis Nature Commun. 3 808
Poland, J.A., Balint-Kurti, P.J., Wisser, R.J., Pratt, R.C. & Nelson, R.J. 2009 Shades of gray: The world of quantitative disease resistance Trends Plant Sci. 14 21 29
Prothro, J., Abdel-Haleem, H., Bachlava, E., White, C., Knapp, S. & McGregor, C. 2013 Quantitative trait loci associated with sex expression in an inter-subspecific watermelon population J. Amer. Soc. Hort. Sci. 38 125 130
Prothro, J., Sandlin, K., Abdel-Haleem, H., Bachlava, E., White, W., Knapp, S. & McGregor, C. 2012a Main and epistatic quantitative trait loci associated with seed size in watermelon J. Amer. Soc. Hort. Sci. 137 452 457
Prothro, J., Sandlin, K., Gill, R., Bachlava, E., White, V., Knapp, S. & McGregor, C. 2012b Mapping of the egusi seed trait locus (eg) and quantitative trait loci associated with seed oil percentage in watermelon J. Amer. Soc. Hort. Sci. 137 311 315
Raman, H., Raman, R., Eckermann, P., Coombes, N., Manoli, S., Zou, X., Edwards, D., Meng, J., Prangnell, R., Stiller, J., Batley, J., Luckett, D., Wratten, N. & Dennis, E. 2013 Genetic and physical mapping of flowering time loci in canola (Brassica napus L.) Theor. Appl. Genet. 126 119 132
Ren, Y., Zhao, H., Kou, Q., Jiang, J., Guo, S., Zhang, H., Hou, W., Zou, X., Sun, H., Gong, G., Levi, A. & Xu, Y. 2012 A high resolution genetic map anchoring scaffolds of the sequenced watermelon genome PLoS One 7 e29453
Richards, R.A. 2006 Physiological traits used in the breeding of new cultivars for water-scarce environments Agr. Water Mgt. 80 197 211
Rieseberg, L.H., Archer, M.A. & Wayne, R.K. 1999 Transgressive segregation, adaptation and speciation Heredity 83 363 372
Sandlin, K.C., Prothro, J.M., Heesacker, A.F., Khalilian, N., Okashah, R., Xiang, W., Bachlava, E., Caldwell, D., Seymour, D., White, V., Chan, E., Tolla, G., White, C., Safran, D., Graham, E., Knapp, S.J. & McGregor, C.E. 2012 Comparative mapping in watermelon [Citrullus lanatus (Thunb.) Matsum. et Nakai] Theor. Appl. Genet. 125 1603 1618
Siddique, K.H.M., Tennant, D., Perry, M.W. & Belford, R.K. 1990 Water use and water use efficiency of old and modern wheat cultivars in a Mediterranean-type environment Aust. J. Agr. Res. 41 431 447
Silva, L.D., Wang, S. & Zeng, Z.-B. 2012 Composite interval mapping and multiple interval mapping: Procedures and guidelines for using windows QTL Cartographer. Quantitative trait loci (QTL) Methods Protocols 871 75 119
Turck, F., Fornara, F. & Coupland, G. 2008 Regulation and identity of florigen: FLOWERING LOCUS T moves center stage Annu. Rev. Plant Biol. 59 573 594
U.S. Department of Agriculture 2011 U.S. Department of Agriculture - National Agricultural Statistics Service. Agricultural statistics. 16 Jan. 2013. <http://www.nass.usda.gov/Publications/Ag_Statistics/index.asp>
Valverde, F., Mouradov, A., Soppe, W., Ravenscroft, D., Samach, A. & Coupland, G. 2004 Photoreceptor regulation of CONSTANS protein in photoperiodic flowering Science 303 1003 1006
Voorrips, R.E. 2002 MapChart: Software for the graphical presentation of linkage maps and QTLs J. Hered. 93 77 78
Wang, S., Basten, C.J. & Zeng, Z.B. 2011 Windows QTL Cartographer 2.5. Dept. Statistics, North Carolina State Univ., Raleigh, NC
Wehner, T. 2008 Watermelon, p. 381–418. In: Prohens, J. and F. Nuez (eds.). Vegetables I: Asteraceae, Brassicaceae, Chenopodicaceae, and Cucurbitaceae. Springer, New York, NY
Wigge, P.A., Kim, M.C., Jaeger, K.E., Busch, W., Schmid, M. & Lohmann, J.U. 2005 Integration of spatial and temporal information during floral induction in Arabidopsis Science 309 1056 1059
Wills, D.M. & Burke, J.M. 2007 Quantitative trait locus analysis of the early domestication of sunflower Genetics 176 2589 2599
Yan, L., Fu, D., Li, C., Blechl, A., Tranquilli, G., Bonafede, M., Sanchez, A., Valarik, M., Yasuda, S. & Dubcovsky, J. 2006 The wheat and barley vernalization gene VRN3 is an orthologue of FT Proc. Natl. Acad. Sci. USA 103 19581 19586
Yano, M., Kojima, S., Takahashi, Y., Lin, H. & Sasaki, T. 2001 Genetic control of flowering time in rice, a short-day plant Plant Physiol. 127 1425 1429
Yoo, S.K., Chung, K.S., Kim, J., Lee, J.H. & Hong, S.M. 2005 CONSTANS activates SUPPRESSOR OF OVEREXPRESSION OF CONSTANS 1 through FLOWERING LOCUS T to promote flowering in Arabidopsis Plant Physiol. 139 770 778
Zalapa, J., Staub, J. & McCreight, J.D. 2008 Variance component analysis of plant architectural traits and fruit yield in melon Euphytica 162 129 143
Zeng, Z.B. 1994 Precision mapping of quantitative trait loci Genetics 136 1457 1468
Zeng, Z.B., Kao, C.H. & Basten, C.J. 1999 Estimating the genetic architecture of quantitative traits Genet. Res. 74 279 289
Zhang, D., Cheng, H., Hu, Z., Wang, H., Kan, G., Liu, C. & Yu, D. 2013 Fine mapping of a major flowering time QTL on soybean chromosome 6 combining linkage and association analysis Euphytica 191 23 33