Evidence for two domestication lineages supporting a middle-eastern origin for Brassica oleracea crops from diversified kale populations

Abstract

Brassica oleracea displays enormous phenotypic variation, including vegetables like cabbage, broccoli, cauliflower, kohlrabi, kales etc. Its domestication has not been clarified, despite several genetic studies and investigations of ancient literature. We used 14 152 high-quality SNP markers for population genetic studies and species-tree estimation (treating morphotypes as “species”) using SVD-quartets coalescent-modelling of a collection of 912 globally distributed accessions representing ten morphotypes of B. oleracea, wild B. oleracea accessions and nine related C9 Brassica species. Our genealogical tree provided evidence for two domestication lineages, the “leafy head” lineage (LHL) and the “arrested inflorescence” lineage (AIL). It also showed that kales are polyphyletic with regards to B. oleracea morphotypes, which fits ancient literature describing highly diverse kale types at around 400 BC. The SVD-quartets species tree topology showed that different kale clades are sister to either the LHL or the AIL. Cabbages from the middle-east formed the first-branching cabbage-clade, supporting the hypothesis that cabbage domestication started in the middle-east, which is confirmed by archeological evidence and historic writings. We hypothesize that cabbages and cauliflowers stem from kales introduced from Western Europe to the middle-east, possibly transported with the tin-trade routes in the Bronze age, to be re-introduced later into Europe. Cauliflower is the least diverse morphotype showing strong genetic differentiation with other morphotypes except broccoli, suggesting a strong genetic bottleneck. Genetic diversity reduced from landraces to modern hybrids for almost all morphotypes. This comprehensive Brassica C-group germplasm collection provides valuable genetic resources and a sound basis for B. oleracea breeding.

References

[1] Bonnema G, Del Carpio DP, Zhao JJ . Diversity analysis and molecular taxonomy of Brassica vegetable crops. In: Kole C, Sadowski J, eds, Genetics, genomics and breeding of crop Plants.Enfield, USA: Science Publishers, 2011, 81-124.

[2] Kole C, Henry RJ . Genetics, Genomics and Breeding of Crop Plants. USA: Science Publishers; 2010.

[3] Dias JS . Portuguese perennial kale: a relic leafy vegetable crop. Genet Resour Crop Evol. 2012; 59: 1201-6.

[4] Guo N, Han S, Zong M et al. Identification and differential expression analysis of anthocyanin biosynthetic genes in leaf color variants of ornamental kale. BMC Genomics. 2019; 20: 564.

[5] Crozier AA . The Cauliflower. USA: Ann Arbor, Mich., Register Publishing Company; 1891.

[6] Maggioni L. Domestication of Brassica Oleracea L. [Ph.D. thesis]. Alnarp, Sweden: Swedish University of Agricultural Sciences. 2015.

[7] Tribulato A, Donzella E, Sdouga D et al. Bio-morphological characterization of Mediterranean wild and cultivated Brassica species. In: VII International Symposium on Brassicas Eds ME Cartea, P Soengas and M Francisco. Spain. Acta Horticulturae 1202. 2017, 9-16.

[8] Snogerup S, Gustafsson M, Von Bothmer R . Brassica sect. Brassica (Brassicaceae) I. taxonomy and variation. Willdenowia. 1990; 19: 271-365.

[9] Kianian SF, Quiros CF . Trait inheritance, fertility, and genomic relationships of some n=9 brassica species. Genet Resour Crop Evol. 1992; 39: 165-75.

[10] Bothmer R, Gustafsson M, Snogerup S . Brassica sect. Brassica (Brassicaceae). II. Inter-and intraspecific crosses with cultivars of B. oleracea. Genet Resour Crop Evol. 1995; 42: 165-78.

[11] Mabry ME, Turner-Hissong SD, Gallagher EY et al. The evolutionary history of wild, domesticated, and feral brassica Oleracea (Brassicaceae). Mol Biol Evol. 2021; 38: 4419-34.

[12] Maggioni L, von Bothmer R, Poulsen G et al. Genetic diversity and population structure of leafy kale and Brassica rupestris Raf. In South Italy. Hereditas. 2014; 151: 145-58.

[13] Maggioni L, von Bothmer R, Poulsen G, Lipman E . Domestication, diversity and use of Brassica oleracea L., based on ancient Greek and Latin texts. Genet Resour Crop Evol. 2018; 65: 137-59.

[14] Gómez-Campo C, Prakash S. 2 Origin and Domestication in Developments in Plant Genetics and Breeding Vol. 4. Edited by Gómez-Campo C . Amsterdam: Elsevier; 1999: 33-58.

[15] Berger D, Soles JS, Giumlia-Mair AR et al. Isotope systematics and chemical composition of tin ingots from Mochlos (Crete) and other late bronze age sites in the eastern Mediterranean Sea: An ultimate key to tin provenance? PLoS One. 2019; 14: e0218326.

[16] Penhallurick RD. Tin in Antiquity: Its Mining and Trade throughout the Ancient World with Particular Reference to Cornwall. Edited by Penhallurick RD . London: Institute of Metal; 2008.

[17] Couvreur TL, Franzke A, Al-Shehbaz IA et al. Molecular phylogenetics, temporal diversification, and principles of evolution in the mustard family (Brassicaceae). Mol Biol Evol. 2010; 27: 55-71.

[18] Guo N, Wang S, Gao L et al. Genome sequencing sheds light on the contribution of structural variants to Brassica oleracea diversification. BMC Biol. 2021; 19: 1-15.

[19] Zeven A . Sixteenth to eighteenth century depictions of cole crops, (Brassica oleracea L.), turnips (B. Rapa L. cultivar group vegetable turnip) and radish (Raphanus sativus L.) from Flan Ders and the present-day Netherlands. In: ISHS Brassica Symposium-IX Crucifer Genetics Workshop 407. (Ed. J.S. Dias, I. Crute, A.A. Monteiro). Acta Hort. 1996, 29-34.

[20] El-Esawi MA, Germaine K, Bourke P, Malone R . Genetic diversity and population structure of Brassica oleracea germplasm in Ireland using SSR markers. C R Biol. 2016; 339: 133-40.

[21] Izzah NK, Lee J, Perumal S et al., eds. Microsatellite-based analysis of genetic diversity in 91 commercial Brassica oleracea L. cultivars belonging to six varietal groups. Genet Resour Crop Evol. 2013; 60: 1967-86.

[22] van Hintum TJ, van de Wiel CC, Visser DL et al. The distribution of genetic diversity in a Brassica oleracea gene bank collection related to the effects on diversity of regeneration, as measured with AFLPs. Theor Appl Genet. 2007; 114: 777-86.

[23] Stansell Z, Hyma K, Fresnedo-Ramírez J et al. Genotyping-by-sequencing of Brassica oleracea vegetables reveals unique phylogenetic patterns, population structure and domestication footprints. Hortic Res. 2018; 5: 1-10.

[24] Cheng F, Sun R, Hou X et al. Subgenome parallel selection is associated with morphotype diversification and convergent crop domestication in Brassica rapa and Brassica oleracea. Nat Genet. 2016; 48: 1218-24.

[25] Pelc SE, Couillard DM, Stansell ZJ, Farnham MW . Genetic diversity and population structure of collard landraces and their relationship to other Brassica oleracea crops. Plant Genome. 2015; 8: eplantgenome2015.04.0023.

[26] Farnham M, Davis E, Morgan J, Smith J . Neglected landraces of collard (Brassica oleracea L. var. viridis) from the Carolinas (USA). Genet Resour Crop Evol. 2008; 55: 797-801.

[27] Warwick SI, Sauder CA . Phylogeny of tribe Brassiceae (Brassicaceae) based on chloroplast restriction site polymorphisms and nuclear ribosomal internal transcribed spacer and chloroplast trn L intron sequences. Can J Bot. 2005; 83: 467- 83.

[28] Cheng F, Jian W, Cai C et al. Genome resequencing and comparative variome analysis in a Brassica rapa and Brassica oleracea collection. Scientific data. 2016; 3: 1-9.

[29] Pritchard JK, Stephens M, Donnelly P . Inference of population structure using multilocus genotype data. Genetics. 2000; 155: 945-59.

[30] Chifman J, Kubatko L . Quartet inference from SNP data under the coalescent model. Bioinformatics. 2014; 30: 3317-24.

[31] Leaché AD, Oaks JR . The utility of single nucleotide polymorphism (SNP) data in phylogenetics. Annu Rev Ecol Evol Syst. 2017; 48: 69-84.

[32] Zhao J, Wang X, Deng B et al. Genetic relationships within Brassica rapa as inferred from AFLP fingerprints. Theor Appl Genet. 2005; 110: 1301-14.

[33] Chen DH, Ronald P. A rapid DNA minipreparation method suitable for AFLP and other PCR applications. Plant Mol Biol Report. 1999; 17: 53-7.

[34] Truong HT, Ramos AM, Yalcin F et al. Sequence-based genotyping for marker discovery and co-dominant scoring in germplasm and populations. PLoS One. 2012; 7(5): e37565.

[35] Kamvar ZN, Tabima JF, Grünwald NJ . Poppr: an R package for genetic analysis of populations with clonal, partially clonal, and/or sexual reproduction. PeerJ. 2014; 2: e281.

[36] Nguyen LT, Schmidt HA, Von Haeseler A, Minh BQ . IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol Biol Evol. 2015; 32: 268-74.

[37] Letunic I, Bork P . Interactive tree of life (iTOL) v3: an online tool for the display and annotation of phylogenetic and other trees. Nucleic Acids Res. 2016; 44: W242-5.

[38] Price AL, Patterson NJ, Plenge RM et al. Principal components analysis corrects for stratification in genome-wide association studies. Nat Genet. 2006; 38: 904-9.

[39] Patterson N, Price AL, Reich D . Population structure and eigenanalysis. PLoS Genet. 2006; 2: e190.

[40] Kopelman NM, Mayzel J, Jakobsson M et al. Clumpak: a program for identifying clustering modes and packaging population structure inferences across K. Mol Ecol Resour. 2015; 15: 1179-91.

[41] Evanno G, Regnaut S, Goudet J . Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol Ecol. 2005; 14: 2611-20.

[42] Maddison W, Maddison D . Mesquite: a modular system for evolutionary analysis Version 3.61. 2019.

[43] Nei M, Li WH . Mathematical model for studying genetic variation in terms of restriction endonucleases. Proc Nat Acad Sci U S A. 1979; 76: 5269-73.

[44] Danecek P, Auton A, Abecasis G et al. The variant call format and VCFtools. Bioinformatics. 2011; 27: 2156-8.

[45] Xu X, Liu X, Ge S et al. Resequencing 50 accessions of cultivated and wild rice yields markers for identifying agronomically important genes. Nat Biotechnol. 2012; 30: 105-11.

[46] Weir BS, Cockerham CC . Estimating F-statistics for the analysis of population structure. Evolution. 1984; 38: 1358-70.

[47] Del Carpio DP, Basnet RK, De Vos RCH et al. The patterns of population differentiation in a Brassica rapa core collection. Theor Appl Genet. 2011; 122: 1105-18.

[48] Bird KA, An H, Gazave E et al. Population structure and phylogenetic relationships in a diverse panel of Brassica rapa L. Front Plant Sci. 2017; 8: 321.

[49] Hartl D. Principles of Population Genetics 488 pp. Eds Hartl D. and Clark A Sunderland, MA: Sinauer Associates; 1980.

[50] Mohammadi SA, Prasanna B . Analysis of genetic diversity in crop plants-salient statistical tools and considerations. Crop Sci. 2003; 43: 1235-48.