A divergent haplotype with a large deletion at the berry color locus causes a white-skinned phenotype in grapevine

Abstract

The current genetic model explaining berry skin color in Vitis vinifera is incomplete and fails to predict berry skin color phenotypes for one allele of VvMybA1, referred to as VvMybA1_SUB. Our study focuses on this specific allele, revealing that the haplotype containing VvMybA1_SUB (haplotype F) represents an ancient lineage of the berry color locus. Within haplotype F, we identified two functional subhaplotypes, HapF1 and HapF2, associated with black-skinned phenotype, and one non-functional subhaplotype, HapF^DEL, responsible for white-skinned phenotype. HapF1 likely originated from wild populations domesticated in the Near East and subsequently spread globally with the expansion of viticulture. In contrast, HapF2 has a more restricted distribution and may have emerged from hybridization events between cultivated grapevines and local wild populations as viticulture migrated to the Italian peninsula. Furthermore, we found that in white-skinned berry cultivar, HapF has undergone a large deletion at the berry color locus, removing the majority of the VvMybA genes. Previous works suggested a single common origin for white-skinned varieties during grapevine domestication. Our results challenge this notion, instead proposing that white-skinned grape cultivars arose at least twice during grapevine domestication history. Alongside the major haplotype A, some white-skinned cultivars, such as cv. ‘Sultanina’ harbor HapF^DEL. Since HapF^DEL is present only in table grape varieties, we suggest that it likely arose from a recent mutational event and dispersed along the ancient Silk Road into East Asia. These findings enhance our understanding of the genetic diversity and evolutionary trajectory of grapevine cultivars, offering insights into their domestication and spread across different geographical regions.

Acknowledgements

The authors thank Cécile Marchal, director of INRAE grapevine collection of Vassal-Montpellier for the plant material and documentation, Dr Olivier Schumpp, head of the virology group at Agroscope for providing support and facilities, Dr Agnes Dienesné Nagy for analyzing wine anthocyanin content and Dr Franco Röckel for sharing plant material from the Geilweilerhof (Germany) germplasm and Zecca Odoardo from the Institut Agricole in Aosta (Italy) for exchanging plant material.

Author contributions

J.S.R. designed the experiment, managed the project, and wrote the manuscript. J.L.S., V.Z., T.L. and D.C. provided support, facilities, and financial resources. J.B. conducted the laboratory work for genomic sequencing and PCR analysis. T.L. and J.S.R. collected biological samples. D.R., J.B., and J.S.R. performed the ONT sequencing and transcriptomic analysis. D.R. and J.S.R performed the data analysis. T.V. conducted the microvinifications. K.A. was responsible for genetic mapping and QTL analysis. E.D. performed the microsatellite analysis. D.C. conducted the WGS experiment and provided guidance on genome evolution. All authors reviewed, revised, and approved the final manuscript.

Data availability

DNA sequencing and RNA-Seq datasets were deposited in the Sequence Read Archive under Bioproject PRJNA1166439.

Conflict of interest statement

The authors declare no conflicts of interest.

Supplementary Data

Supplementary data is available at Horticulture Research online.

References

[1] Vogt T. Phenylpropanoid biosynthesis. Mol Plant. 2010; 3:2-20

[2] Albert NW, Lafferty DJ, Moss SMA. et al. Flavonoids—flowers, fruit, forage and the future. JRSocNZ. 2023; 53:304-31

[3] Slinkard K, Singleton V. Phenol content of grape skins and the loss of ability to make anthocyanins by mutation. Vitis. 1984; 23: 173-5

[4] Guasch-Jané MR, Andrés-Lacueva C, Jáuregui O. et al. First evi-dence of white wine in ancient Egypt from Tutankhamun’s tomb. JArchaeolSci. 2006; 33:1075-80

[5] Levadoux L. La sélection et l’hybridation Chez la Vigne. Montpellier: Imprimerie Charles Déhan; 1951:

[6] Sparvoli F, Martin C, Scienza A. et al. Cloning and molecular analysis of structural genes involved in flavonoid and stilbene biosynthesis in grape (Vitis vinifera L.). Plant Mol Biol. 1994; 24: 743-55

[7] Boss PK, Davies C, Robinson SP. Anthocyanin composition and anthocyanin pathway gene expression in grapevine sports dif-fering in berry skin colour. Aust J Grape Wine Res. 1996; 2:163-70

[8] Kobayashi S, Ishimaru M, Hiraoka K. et al. Myb-related genes of theKyoho grape(Vitis labruscana) regulate anthocyanin biosyn-thesis. Planta. 2002; 215:924-33

[9] Klempnauer K-H, Gonda TJ, Michael Bishop J. Nucleotide sequence of the retroviral leukemia gene v-myb and its cellular progenitor c-myb: the architecture of a transduced oncogene. Cell. 1982; 31:453-63

[10] Paz-Ares J, Ghosal D, Wienand U. et al. The regulatory c1 locus of Zea mays encodes a protein with homology to myb proto-oncogene products and with structural similarities to transcrip-tional activators. EMBO J. 1987; 6:3553-8

[11] Doligez A, Bouquet A, Danglot Y. et al. Genetic mapping of grapevine (Vitis vinifera L.) applied to the detection of QTLs for seedlessness and berry weight. Theor Appl Genet. 2002; 105: 780-95

[12] Walker AR, Lee E, Bogs J. et al. White grapes arose through the mutation of two similar and adjacent regulatory genes. Plant J. 2007; 49:772-85

[13] Fournier-Level A, et al.Lacombe T, le Cunff L. Evolution of the VvMybA gene family, the major determinant of berry colour in cultivated grapevine (Vitis vinifera L.). Heredity. 2010; 104:351-62

[14] Azuma A. Genetic and environmental impacts on the biosyn-thesis of anthocyanins in grapes. The Horticulture Journal. 2018; 87: 1-17

[15] Kobayashi S, Goto-Yamamoto N, Hirochika H. Retrotransposon-induced mutations in grape skin color. Science. 2004; 304:982

[16] Vezzulli S, Leonardelli L, Malossini U. et al. Pinot blanc and Pinot gris arose as independent somatic mutations of Pinot noir. JExp Bot. 2012; 63:6359-69

[17] Walker AR, Lee E, Robinson SP. Two new grape cultivars, bud sports of Cabernet Sauvignon bearing pale-coloured berries, are the result of deletion of two regulatory genes of the berry colour locus. Plant Mol Biol. 2006; 62:623-35

[18] Carbonell-Bejerano P, Royo C, Torres-Pérez R. et al. Catastrophic unbalanced genome rearrangements cause somatic loss of berry color in grapevine. Plant Physiol. 2017; 175:786-801

[19] Lijavetzky D, Ruiz-García L, Cabezas JA. et al. Molecular genet-ics of berry colour variation in table grape. Mol Gen Genomics. 2006; 276:427-35

[20] Azuma A, Kobayashi S, Mitani N. et al. Genomic and genetic analysis of Myb-related genes that regulate anthocyanin biosyn-thesis in grape berry skin. Theor Appl Genet. 2008; 117:1009-19

[21] Vouillamoz J, Maigre D, Meredith CP. Microsatellite analysis of ancient alpine grape cultivars: pedigree reconstruction of Vitis vinifera L. ’Cornalin du Valais’. Theor Appl Genet. 2003; 107: 448-54

[22] Raimondi S, Tumino G, Ruffa P. et al. DNA-based genealogy reconstruction of Nebbiolo, Barbera and other ancient grapevine cultivars from northwestern Italy. Sci Rep. 2020; 10:15782

[23] D’Onofrio C, Tumino G, Gardiman M. et al. Parentage atlas of Italian grapevine varieties as inferred from SNP genotyping. Front. Plant Sci. 2021; 11:11

[24] Lacombe T, Boursiquot JM, Laucou V. et al. Large-scale parentage analysis in an extended set of grapevine cultivars (Vitis vinifera L.). Theor Appl Genet. 2013; 126:401-14

[25] Cipriani G, Spadotto A, Jurman I. et al. The SSR-based molecular profile of 1005 grapevine (Vitis vinifera L.) accessions uncovers new synonymy and parentages, and reveals a large admixture amongst varieties of different geographic origin. Theor Appl Genet. 2010; 121:1569-85

[26] Balic I, Vizoso P, Nilo-Poyanco R. et al. Transcriptome analysis during ripening of table grape berry cv. Thompson seedless. PLoS One. 2018; 13:e0190087

[27] Massonnet M, Cochetel N, Minio A. et al. The genetic basis of sex determination in grapes. Nat Commun. 2020; 11:2902

[28] Cohen P, Bacilieri R, Ramos-Madrigal J. et al. Ancient DNA from a lost Negev Highlands desert grape reveals a late antiquity wine lineage. Proc Natl Acad Sci. 2023; 120:e2213563120

[29] Cochetel N, Minio A, Massonnet M. et al. Diploid chromosome-scale assembly of the Muscadinia rotundifolia genome supports chromosome fusion and disease resistance gene expansion dur-ing Vitis and Muscadinia divergence. G3 (Bethesda). 2021;11

[30] Park M, Darwish AG, Elhag RI. et al. A multi-locus genome-wide association study reveals the genetics underlying mus-cadine antioxidant in berry skin. Front Plant Sci. 2022; 13: 969301

[31] Lacombe T., Contribution à l’étude de l’histoire évolutive de la vigne cultivée (Vitis vinifera L.) par l’analyse de la diversité génétique neutre et de gènes d’intérêt. 2012, Centre International d’Etudes Supérieures en Sciences Agronomiques: Montpellier p. 328.

[32] Margaryan K, De Lorenzis G, Aroutiounian R. et al. Genetic analysis of the VvMybA1 gene in Armenian grapevines (V. vinifera L.). Vitis. 2015; 54:49-50

[33] Alejandra Milla T, Cabezas JA, Cabello F. et al. Determining the Spanish origin of representative ancient American grapevine varieties. Am J Enol Vitic. 2007; 58:242-51

[34] Lalla Hasna Z, Haddioui A, Rodríguez V. et al. Identification by SNP analysis of a major role for Cayetana Blanca in the genetic network of Iberian Peninsula grapevine varieties. Am J Enol Vitic. 2012; 63:121-6

[35] Berget A. Etude ampélographique des vignobles du Léman, du Valais et du val d’Aoste. Revue du viticulture. 1903;502

[36] Vouillamoz J, Monaco A, Costantini L. et al. The parentage of ’Sangiovese’, the most important Italian wine grape. Vitis. 2007; 46:19-22

[37] Bergamini C, Caputo AR, Gasparro M. et al. Evidences for an alternative genealogy of ’Sangiovese’. Mol Biotechnol. 2013; 53: 278-88

[38] Crespan M, Mercati F, De Lorenzis. et al. Two main distinct evolutionary stories describe the Italian grapevine assortment. Vitis. 2023; 62:30-6

[39] Jiu S, Guan L, Leng X. et al. The role of VvMYBA2r and VvMYBA2w alleles of the MYBA2 locus in the regulation of anthocyanin biosynthesis for molecular breeding of grape (Vitis spp.) skin coloration. Plant Biotechnol J. 2021; 19:1216-39

[40] Azuma A, Udo Y, Sato A. et al. Haplotype composition at the color locus is a major genetic determinant of skin color variation in Vitis × labruscana grapes. Theor Appl Genet. 2011; 122:1427-38

[41] Shimazaki M, Fujita K, Kobayashi H. et al. Pink-colored grape berry is the result of short insertion in intron of color regulatory gene. PLoS One. 2011; 6:e21308

[42] Butelli E, Garcia-Lor A, Licciardello C. et al. Changes in antho-cyanin production during domestication of citrus. Plant Physiol. 2017; 173:2225-42

[43] Castillejo C, Waurich V, Wagner H. et al. Allelic variation of MYB10 is the major force controlling natural variation in skin and flesh color in strawberry (Fragaria spp.) fruit. Plant Cell. 2020; 32:3723-49

[44] Liu C, Qi X, Song L. et al. Large-fragment deletion encompasses the R2R3 MYB transcription factor, PavMYB10.1, causes yellow fruits in sweet cherry ( Prunus avium L.). Sci Hortic. 2023; 309: 111648

[45] Zhou Y, Minio A, Massonnet M. et al. The population genetics of structural variants in grapevine domestication. Nature Plants. 2019; 5:965-79

[46] Röckel F, Moock C, Schwander F. et al. A 69 kbp deletion at the berry color locus is responsible for berry color recovery in Vitis vinifera L. cultivar ‘Riesling Rot’. Int J Mol Sci. 2022; 23:23

[47] Röckel F. Färberreben (Teinturiers) Sowie Rote Farbmutanten weißer Qualitätsrebsorten Entstanden Durch VvmybA-Bedingte Mutationen am Beerenfarblokus. Julius Kühn-Institut: Quedlinburg (Deutsch-land); 2017:128

[48] Maigre D. Humagne gris: première observation d’une forme grise de l’Humagne rouge en Valais. Revue suisse de viticulture, arboriculture, horticulture. 2005; 37:153-5

[49] Pelsy F. Molecular and cellular mechanisms of diversity within grapevine varieties. Heredity. 2010; 104:331-40

[50] Gatta LF. Saggio Sulle Viti e Sui Vini Della Valle d’Aosta. R. Società Agraria; Torino, 1838:

[51] Zhou Y, Massonnet M, Sanjak JS. et al. Evolutionary genomics of grape (Vitis vinifera ssp. vinifera) domestication. Proc Natl Acad Sci. 2017; 114:11715-20

[52] Jiang H-E, Zhang YB, Li X. et al. Evidence for early viticulture in China: proof of a grapevine (Vitis vinifera L., Vitaceae) in the Yanghai Tombs, Xinjiang. JArchaeolSci. 2009; 36:1458-65

[53] Dong Y, Duan S, Xia Q. et al. Dual domestications and origin of traits in grapevine evolution. Science. 2023; 379:892-901

[54] Goto-Yamamoto N, Sawler J, Myles S. Genetic analysis of East Asian grape cultivars suggests hybridization with wild Vitis. PLoS One. 2015; 10:e0140841

[55] Myles S, Boyko AR, Owens CL. et al. Genetic structure and domestication history of the grape. Proc Natl Acad Sci. 2011; 108: 3530-5

[56] Magris G, Jurman I, Fornasiero A. et al. The genomes of 204 Vitis vinifera accessions reveal the origin of European wine grapes. Nat Commun. 2021; 12:7240

[57] Xiao H, Liu Z, Wang N. et al. Adaptive and maladaptive introgres-sion in grapevine domestication. Proc Natl Acad Sci USA. 2023; 120:e2222041120

[58] Bouby L, Chabal L, Bonhomme V. et al. The Holocene history of grapevine (Vitis vinifera) and viticulture in France retraced from a large-scale archaeobotanical dataset. Palaeogeogr Palaeoclimatol Palaeoecol. 2023; 625:111655

[59] D’Onofrio C. Introgression among cultivated and wild grapevine in Tuscany. Frontiers. Plant Sci. 2020; 11:11

[60] Freitas S, Gazda MA, Rebelo MÂ. et al. Pervasive hybridization with local wild relatives in Western European grapevine vari-eties. Sci Adv. 7:eabi8584

[61] Vivas N, Vivas N, Nonier M-F. Estimation and quantification of wine phenolic compounds. Bull Org Int Vigne Vin. 2003; 76:281-303

[62] Gambino G, Perrone I, Gribaudo I. A rapid and effective method for RNA extraction from different tissues of grapevine and other woody plants. Phytochem Anal. 2008; 19:520-5

[63] Debonneville C, Mandelli L, Brodard J. et al. The complete genome of the “Flavescence Dorée” Phytoplasma reveals char-acteristics of low genome plasticity. Biology (Basel). 2022;11:

[64] Kolmogorov M, Yuan J, Lin Y. et al. Assembly of long, error-prone reads using repeat graphs. Nat Biotechnol. 2019; 37:540-6

[65] Kearse M, Moir R, Wilson A. et al. Geneious basic: an integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics. 2012; 28:1647-9

[66] Su W., Ou S., Hufford M.B., Peterson T., A tutorial of EDTA: extensive De novo TE annotator. Methods Mol Biol, 2021. 2250:p. 55-67, New York, NY, Springer US.

[67] SanMiguel P, Gaut BS, Tikhonov A. et al. The paleontology of intergene retrotransposons of maize. Nat Genet. 1998; 20:43-5

[68] Cabanettes F, Klopp C. D-GENIES: dot plot large genomes in an interactive, efficient and simple way. PeerJ. 2018; 6:e4958

[69] Krumsiek J, Arnold R, Rattei T. Gepard: a rapid and sensitive tool for creating dotplots on genome scale. Bioinformatics. 2007; 23: 1026-8

[70] McKenna A, Hanna M, Banks E. et al. The genome analysis toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010; 20: 1297-303

[71] Chang CC, Chow CC, Tellier LC. et al. Second-generation PLINK: rising to the challenge of larger and richer datasets. Gigascience. 2015; 4:7

[72] Canaguier A, Grimplet J, di Gaspero G. et al. A new version of the grapevine reference genome assembly (12X.v2) and of its annotation (VCost.v3). Genom Data. 2017; 14:56-62

[73] Patro R, Duggal G, Love MI. et al. Salmon provides fast and bias-aware quantification of transcript expression. Nat Methods. 2017; 14:417-9

[74] Elshire RJ, Glaubitz JC, Sun Q. et al.A robust, simple genotyping-by-sequencing (GBS) approach for high diversity species. PLoS One. 2011; 6:e19379

[75] Rastas P. Lep-MAP3: robust linkage mapping even for low-coverage whole genome sequencing data. Bioinformatics. 2017; 33: 3726-32

[76] Broman KW, Wu H, Sen Ś. et al. R/QTL: QTL mapping in experi-mental crosses. Bioinformatics. 2003; 19:889-90