New insights into flowering date in Prunus: fine mapping of a major QTL in sweet cherry

Abstract

Flowering date is an important trait in Prunus fruit species, especially for their adaptation in a global warming context. Numerous quantitative trait loci (QTLs) have been identified and a major one was previously located on LG4. The objectives of this study were to fine-map this QTL in sweet cherry, to identify robust candidate genes by using the new sweet cherry genome sequence of the cultivar “Regina” and to define markers usable in marker-assisted selection (MAS). We performed QTL analyses on two populations derived from crosses using cultivars “Regina” and “Garnet” as parents. The first one (n = 117) was phenotyped over ten years, while the second one (n = 1386) was evaluated during three years. Kompetitive allele specific PCR (KASP) markers located within the QTL region on LG4 were developed and mapped within this region, consisting in the first fine mapping in sweet cherry. The QTL interval was narrowed from 380 kb to 68 kb and candidate genes were identified by using the genome sequence of “Regina”. Their expression was analyzed from bud dormancy period to flowering in cultivars “Regina” and “Garnet”. Several genes, such as PavBOI-E3, PavSR45a and PavSAUR71, were differentially expressed in these two cultivars and could be then considered as promising candidate genes. Two KASP markers were validated using a population derived from a cross between cultivars “Regina” and “Lapins” and two collections, including landraces and modern cultivars. Thanks to the high synteny within the Prunus genus, these results give new insights into the control of flowering date in Prunus species and pave the way for the development of molecular breeding strategies.

References

[1] Quero-García J, Schuster M, López-Ortega G, Charlot G. Sweet cherry varieties and improvement. In: Quero-García J, Iezzoni A, Pulawska J, Lang G, eds. Cherries: Botany, Production and Uses.CABI (UK), 2017, 60-94.

[2] Radičević S, Cerović R, Nikolić D, Đorđević M . The effect of genotype and temperature on pollen tube growth and fertilization in sweet cherry (Prunus avium L.). Euphytica. 2016; 209: 121-36.

[3] Wenden B, Campoy JA, Jensen M, López-Ortega G . Climatic limiting factors: temperature. In: Quero-García J, Iezzoni A, Pulawska J, Lang G, eds. Cherries: Botany, Production and Uses.CABI (UK), 2017, 166-88.

[4] Dirlewanger E, Quero-García J, Le Dantec et al. Comparison of the genetic determinism of two key phenological traits, flowering and maturity dates, in three Prunus species: peach, apricot and sweet cherry. Heredity (Edind). 2012; 109: 280-92.

[5] Lang GA, Early JD, Martin GC, Darnell RL . Endo-, Para-, and ecodormancy: physiological terminology and classification for dormancy research. HortScience. 1987; 22: 371-7.

[6] Luedeling E. Climate change impacts on winter chill for temperate fruit and nut production: a review. Sci Hortic. 2012; 144: 218-29.

[7] Couvillon GA, Erez A . Influence of prolonged exposure to chilling temperatures on bud break and heat requirement for bloom of several fruit species. J Am Soc Hortic Sci. 1985; 110: 47-50.

[8] Castède S, Campoy JA, García JQ et al. Genetic determinism of phenological traits highly affected by climate change in Prunus avium: flowering date dissected into chilling and heat requirements. New Phytol. 2014; 202: 703-15.

[9] Calle A, Cai L, Iezzoni A, Wünsch A . Genetic dissection of bloom time in low chilling sweet cherry (Prunus avium L.) using a multi-family QTL approach. Front Plant Sci. 2020; 10: 1647.

[10] Sánchez-Pérez R, Del Cueto J, Dicenta F, Martínez-Gómez P . Recent advancements to study flowering time in almond and other Prunus species. Front Plant Sci. 2014; 5: 334.

[11] Dirlewanger E, Graziano E, Joobeur T et al. Comparative mapping and marker-assisted selection in Rosaceae fruit crops. Proc Natl Acad Sci U S A. 2004; 101: 9891-6.

[12] Alburquerque N, García-Montiel F, Carrillo A, Burgos L . Chilling and heat requirements of sweet cherry cultivars and the relationship between altitude and the probability of satisfying the chill requirements. Environ Exp Bot. 2008; 64: 162-70.

[13] Sánchez-Pérez R, Dicenta F, Martínez-Gómez P . Inheritance of chilling and heat requirements for flowering in almond and QTL analysis. Tree Genet Genomes. 2012; 8: 379-89.

[14] Campoy JA, Ruiz D, Allderman L et al. The fulfilment of chilling requirements and the adaptation of apricot (Prunus armeniaca L.) in warm winter climates: an approach in Murcia (Spain) and the Western cape (South Africa). Eur J Agron. 2012; 37: 43-55.

[15] Okie WR, Blackburn B . Increasing chilling reduces heat requirement for floral Budbreak in peach. HortScience. 2011; 46: 245-52.

[16] Cai L, Stegmeir T, Sebolt A et al. Identification of bloom date QTLs and haplotype analysis in tetraploid sour cherry (Prunus cerasus). Tree Genet Genomes. 2018; 14: 22.

[17] Tavaud M, Zanetto A, David JL et al. Genetic relationships between diploid and allotetraploid cherry species (Prunus avium, Prunus × gondouinii and Prunus cerasus). Heredity. 2004; 93: 631-8.

[18] Bielenberg DG, Wang Y(E, Li Z et al. Sequencing and annotation of the evergrowing locus in peach [Prunus persica (L.) Batsch] reveals a cluster of six MADS-box transcription factors as candidate genes for regulation of terminal bud formation. Tree Genet Genomes. 2008; 4: 495-507.

[19] Zhebentyayeva TN, Fan S, Chandra A et al. Dissection of chilling requirement and bloom date QTLs in peach using a whole genome sequencing of sibling trees from an F2 mapping population. Tree Genet Genomes. 2014; 10: 35-51.

[20] Castède S, Campoy JA, Le Dantec et al. Mapping of candidate genes involved in bud dormancy and flowering time in sweet cherry (Prunus avium). PLoS One. 2015; 10: e0143250.

[21] Verde I, Abbott AG, Scalabrin S et al. The high-quality draft genome of peach (Prunus persica) identifies unique patterns of genetic diversity, domestication and genome evolution. Nat Genet. 2013; 45: 487-94.

[22] Peace C, Bassil N, Main D et al. Development and evaluation of a genome-wide 6K SNP Array for diploid sweet cherry and tetraploid sour cherry. PLoS One. 2012; 7: e48305.

[23] Vanderzande S, Zheng P, Cai L et al. The cherry 6+9K SNP array: a cost-effective improvement to the cherry 6K SNP array for genetic studies. Sci Rep. 2020; 10: 7613.

[24] Vimont N, Fouché M, Campoy JA et al. From bud formation to flowering: transcriptomic state defines the cherry developmental phases of sweet cherry bud dormancy. BMC Genomics. 2019; 20: 974.

[25] Pinosio S, Marroni F, Zuccolo A et al. A draft genome of sweet cherry (Prunus avium L.) reveals genome-wide and local effects of domestication. Plant J. 2020; 103: 1420-32.

[26] Calle A, Serradilla MJ, Wünsch A . QTL mapping of phenolic compounds and fruit colour in sweet cherry using a 6+9K SNP array genetic map. Sci Hortic. 2021; 280: 109900.

[27] Fan S, Bielenberg DG, Zhebentyayeva TN et al. Mapping quantitative trait loci associated with chilling requirement, heat requirement and bloom date in peach (Prunus persica). New Phytol. 2010; 185: 917-30.

[28] Calle A, Grimplet J, Le Dantec L, Wünsch A . Identification and characterization of DAMs mutations associated with Early blooming in sweet cherry, and validation of DNA-based markers for selection. Front Plant Sci. 2021; 12: 621491.

[29] Boudehri K, Bendahmane A, Cardinet G et al. Phenotypic and fine genetic characterization of the D locus controlling fruit acidity in peach. BMC Plant Biol. 2009; 9: 59.

[30] Pirona R, Eduardo I, Pacheco I et al. Fine mapping and identification of a candidate gene for a major locus controlling maturity date in peach. BMC Plant Biol. 2013; 13: 166.

[31] Lu Z, Niu L, Chagné D et al. Fine mapping of the temperature-sensitive semi-dwarf (Tssd) locus regulating the internode length in peach (Prunus persica). Mol Breed. 2016; 36: 20.

[32] Marimon N, Luque J, Arús P, Eduardo I . Fine mapping and identification of candidate genes for the peach powdery mildew resistance gene Vr3. Hortic Res. 2020; 7: 1-9.

[33] Lu Z, Pan L, Wei B et al. Fine mapping of the gene controlling the fruit skin hairiness of Prunus persica and its uses for MAS in progenies. Plants (Basel). 2021; 10: 1433.

[34] Cai L, Quero-García J, Barreneche T et al. A fruit firmness QTL identified on linkage group 4 in sweet cherry (Prunus avium L.) is associated with domesticated and bred germplasm. Sci Rep. 2019; 9: 5008.

[35] Calle A, Wünsch A . Multiple-population QTL mapping of maturity and fruit-quality traits reveals LG4 region as a breeding target in sweet cherry (Prunus avium L.). Hortic Res. 2020; 7: 1-13.

[36] Quero-García J, Letourmy P, Campoy JA et al. Multi-year analyses on three populations reveal the first stable QTLs for tolerance to rain-induced fruit cracking in sweet cherry (Prunus avium L.). Hortic Res. 2021; 8: 136.

[37] Bielenberg DG, Rauh B, Fan S et al. Genotyping by sequencing for SNP-based linkage map construction and QTL analysis of chilling requirement and bloom date in peach [ Prunus persica (L.) Batsch] . PLoS One. 2015; 10: e0139406.

[38] Socias R, Felipe AJ, Aparisi JG . A major gene for flowering time in almond. Plant Breed. 1999; 118: 443-8.

[39] Silva C, Garcia-Mas J, Sánchez AM et al. Looking into flowering time in almond (Prunus dulcis (mill) D. A. Webb): the candidate gene approach. Theor Appl Genet. 2005; 110: 959-68.

[40] Nguyen KT, Park J, Park E et al. The Arabidopsis RING domain protein BOI inhibits flowering via CO-dependent and CO-independent mechanisms. Mol Plant. 2015; 8: 1725-36.

[41] Ali GS, Palusa SG, Golovkin M et al. Regulation of plant developmental processes by a novel splicing factor. PLoS One. 2007; 2: e471.

[42] Chen SL, Rooney TJ, Hu AR et al. Quantitative proteomics reveals a role for SERINE/ARGININE-rich 45 in regulating RNA metabolism and modulating transcriptional suppression via the ASAP complex in Arabidopsis thaliana. Front Plant Sci. 2019; 10: 1116.

[43] Ausin I, Greenberg MVC, Li CF, Jacobsen SE . The splicing factor SR45 affects the RNA-directed DNA methylation pathway in Arabidopsis. Epigenetics. 2012; 7: 29.

[44] Carvalho RF, Szakonyi D, Simpson CG et al. The Arabidopsis SR45 splicing factor, a negative regulator of sugar signaling, modulates SNF1-related protein kinase 1 stability. Plant Cell. 2016; 28: 1910-25.

[45] Stortenbeker N, Bemer M . The SAUR gene family: the plant’s toolbox for adaptation of growth and development. J Exp Bot. 2019; 70: 17-27.

[46] Liu J, Sherif SM . Hormonal orchestration of bud dormancy cycle in deciduous Woody perennials. Front Plant Sci. 2019; 10: 1136.

[47] Sun M, Qian X, Chen C et al. Ectopic expression of GsSRK in Medicago sativa reveals its involvement in plant architecture and salt stress responses. Front Plant Sci. 2018; 9: 226.

[48] Naeem ul Hassan M, Zainal Z, Ismail I . Plant kelch containing F-box proteins: structure, evolution and functions. RSC Adv. 2015; 5: 42808-14.

[49] Liu H, Zhou X, Li Q et al. CCT domain-containing genes in cereal crops: flowering time and beyond. Theor Appl Genet. 2020; 133: 1385-96.

[50] Datta S, Hettiarachchi GHCM, Deng X-W, Holm M . Arabidopsis CONSTANS-LIKE3 is a positive regulator of red light signaling and root growth. Plant Cell. 2006; 18: 70-84.

[51] Tripathi P, Carvallo M, Hamilton EE et al. Arabidopsis B-BOX32 interacts with CONSTANS-LIKE3 to regulate flowering. PNAS. 2017; 114: 172-7.

[52] Singh RK, Svystun T, AlDahmash B et al. Photoperiod- and temperature-mediated control of phenology in trees - a molecular perspective. New Phytol. 2017; 213: 511-24.

[53] Bittner F, Oreb M, Mendel RR . ABA3 is a molybdenum cofactor sulfurase required for activation of aldehyde oxidase and xanthine dehydrogenase in Arabidopsis thaliana. J Biol Chem. 2001; 276: 40381-4.

[54] Pereira-Santana A, Alcaraz LD, Castaño E et al. Comparative genomics of NAC transcriptional factors in angiosperms: implications for the adaptation and diversification of flowering plants. PLoS One. 2015; 10: e0141866.

[55] Migicovsky Z, Gardner KM, Money D et al. Genome to phenome mapping in apple using historical data. Plant Genome. 2016; 9: 1-15.

[56] Shahnejat-Bushehri S, Tarkowska D, Sakuraba Y, Balazadeh S . Arabidopsis NAC transcription factor JUB1 regulates GA/BR metabolism and signalling. Nat Plants. 2016; 2: 1-9.

[57] Wu A, Allu AD, Garapati P et al. JUNGBRUNNEN1, a reactive oxygen species-responsive NAC transcription factor, regulates longevity in Arabidopsis. Plant Cell. 2012; 24: 482-506.

[58] Groppi A, Liu S, Cornille A et al. Population genomics of apricots unravels domestication history and adaptive events. Nat Commun. 2021; 12: 3956.

[59] Campoy JA, Lerigoleur-Balsemin E, Christmann H et al. Genetic diversity, linkage disequilibrium, population structure and construction of a core collection of Prunus avium L. landraces and bred cultivars. BMC Plant Biol. 2016; 16: 49.

[60] Fadón E, Herrero M, Rodrigo J . Flower development in sweet cherry framed in the BBCH scale. Sci Hortic. 2015; 192: 141-7.

[61] Bates D, Mächler M, Bolker B, Walker S . Fitting linear mixed-effects models using lme4 . J Stat Softw. 2015; 67: 1-48.

[62] Klagges C, Campoy JA, Quero-García J et al. Construction and comparative analyses of highly dense linkage maps of two sweet cherry intra-specific progenies of commercial cultivars. PLoS One. 2013; 8: e54743.

[63] Grattapaglia D, Sederoff R . Genetic linkage maps of Eucalyptus grandis and Eucalyptus urophylla using a pseudo-testcross: mapping strategy and RAPD markers. Genetics. 1994; 137: 1121-37.

[64] Van Ooijen JW . JoinMap^® 4, software for the calculation of genetic linkage maps in experimental populations . Wageningen: Kyazma BV, 2006.

[65] Voorrips RE . MapChart: software for the graphical presentation of linkage maps and QTLs. J Hered. 2002; 93: 77-8.

[66] Dirlewanger E, Le-Dantec L, Fouché M et al. New sweet cherry genomic tools and their use in marker-assisted breeding. Acta Hortic. 2020; 533-42.

[67] Bernard A, Marrano A, Donkpegan A et al. Association and linkage mapping to unravel genetic architecture of phenological traits and lateral bearing in Persian walnut (Juglans regia L.). BMC Genomics. 2020; 21: 203.

[68] Le Dantec L, Girollet N, Jérôme G et al. Assembly and annotation of ‘Regina’ sweet cherry genome. 2020.

[69] Conesa A, Götz S, García-Gómez JM et al. Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005; 21: 3674-6.

[70] Bolger AM, Lohse M, Usadel B . Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014; 30: 2114-20.

[71] Vimont N, Schwarzenberg A, Domijan M et al. Fine tuning of hormonal signaling is linked to dormancy status in sweet cherry flower buds. Tree Physiol. 2021; 41: 544-61.

[72] Anders S, Pyl PT, Huber W . HTSeq-a python framework to work with high-throughput sequencing data. Bioinformatics. 2015; 31: 166-9.

[73] Wagner GP, Kin K, Lynch VJ . Measurement of mRNA abundance using RNA-seq data: RPKM measure is inconsistent among samples. Theory Biosci. 2012; 131: 281-5.