Genomic analyses of arugula (Eruca sativa) accessions reveal domestication history and breeding insights

Abstract

Eruca sativa, or arugula, a leafy vegetable from the Mediterranean region, is gaining popularity worldwide. However, its domestication history and genetic diversity remain poorly understood. This study investigates the genetic background of globally collected arugula accessions, their domestication history, and the implications for future breeding strategies. Among the 202 arugula accessions analysed, four distinct groups were identified. Group 1 (G1) appears to be misidentified as arugula and is more closely related to radish (Raphanus sativus) and Brassica spp. Group 2 (G2), Group 3 (G3), and Group 4 (G4) likely represent arugula populations evolved separately, with G4 showing greater genetic isolation and drifts. Geographically, G2 accessions cluster in Europe, G3 in North Africa, and G4 in South Asia (furthest from the Mediterranean region). Phylogenetic and demographic history analyses suggest that G2 is more ancient than G3 and G4. While G3 and G4 share a more recent common ancestor, small scale introgression between G2 and G3 was detected. Selective sweeps signatures revealed 12 genes under positive selection in G4 compared to G2, and these 12 genes are likely associated with adaptive divergency, suggesting that arugula has not undergone stringent artificial selection. A potential hybrid zone between Turkey, Israel, and Iran was identified. Accessions from this zone may serve as a valuable resource for screening new arugula varieties with superior traits. Although 202 accessions were studied, some exhibited high genetic similarity. To enhance germplasm diversity for breeding, future efforts should prioritize expanding collections in regions where arugula has naturalized.

Keywords

Genetic diversity / Domestication history / Hybrid zone / Adaptive divergence / Crop breeding

References

[1] AgrawalN, GuptaM, BangaSS, Heslop-HarrisonJP. Identification of chromosomes and chromosome rearrangements in crop brassicas and Raphanus sativus: a cytogenetic toolkit using synthesized massive oligonucleotide libraries. Front Plant Sci, 2020, 11. 598039

[2] AlexanderDH, NovembreJ, LangeK. Fast model-based estimation of ancestry in unrelated individuals. Genome Res, 2009, 19(9): 1655-1664.

[3] BehrAA, LiuKZ, Liu-FangG, NakkaP, RamachandranS. pong: fast analysis and visualization of latent clusters in population genetic data. Bioinformatics, 2016, 32(18): 2817-2823.

[4] BellL, WagstaffC. Rocket science: a review of phytochemical & health-related research in Eruca & Diplotaxis species. Food Chemistry: X, 2019, 1. 100002

[5] ChangCC, ChowCC, TellierLC, VattikutiS, PurcellSM, LeeJJ. Second-generation PLINK: rising to the challenge of larger and richer datasets. Gigascience, 2015, 417.

[6] ChenH, PattersonN, ReichD. Population differentiation as a test for selective sweeps. Genome Res, 2010, 20(3): 393-402.

[7] DaiL, LiX, GuanC, ZhongJ. On ovary and embryo culture of hybrid between rapeseeds and Rocketsalad. J Hunan Agric Univ (Nat Sci), 2004, 3: 201-204

[8] DanecekP, AutonA, AbecasisG, AlbersCA, BanksE, DePristoMA, et al.. The variant call format and VCFtools. Bioinformatics, 2011, 27(15): 2156-2158.

[9] DanecekP, BonfieldJK, LiddleJ, MarshallJ, OhanV, PollardMO, et al.. Twelve years of SAMtools and BCFtools. Gigascience, 2021, 102giab008.

[10] DongG, DuL, YangL, LuM, QiuM, LiH, et al.. Dispersal of crop-livestock and geographical-temporal variation of subsistence along the Steppe and Silk Roads across Eurasia in prehistory. Sci China Earth Sci, 2022, 65(7): 1187-1210.

[11] GolkarP, BakhtiariMA. Evaluation of genetic diversity in the world collection of Eruca sativa L. using oil content, fatty acids and molecular markers. Ind Crops Prod, 2020, 148. 112280

[12] GreuterW, BurdetHM, LongGMed-Checklist: a critical inventory of vascular plants of the circum-mediterranean countries, 2008, Geneva. OPTIMA and Conservatoire et Jardin botaniques. 3

[13] Guijarro-RealC, NavarroA, EspositoS, FestaG, MacellaroR, Di CesareC, et al.. Large scale phenotyping and molecular analysis in a germplasm collection of rocket salad (Eruca vesicaria) reveal a differentiation of the gene pool by geographical origin. Euphytica, 2020, 216353.

[14] HallM, JoblingJ, RogersG. Factors affecting growth of perennial wall rocket and annual garden rocket. Int J Veg Sci, 2012, 18(4): 393-411.

[15] Hanelt P, Buttner R, Mansfeld R. Mansfeld's encyclopedia of agricultural and horticultural crops (except ornamentals). Springer; 2001. https://doi.org/10.1007/978-3-540-30442-5.

[16] KangL, QianL, ZhengM, ChenL, ChenH, YangL, et al.. Genomic insights into the origin, domestication and diversification of Brassica juncea. Nat Genet, 2021, 53(9): 1392-1402.

[17] KieferC, WillingEM, JiaoWB, SunH, PiednoelM, HumannU, et al.. Interspecies association mapping links reduced CG to TG substitution rates to the loss of gene-body methylation. Nat Plants, 2019, 5(8): 846-55.

[18] Lan M. Dian Nan Ben Cao. 1959.

[19] LiH, DurbinR. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics, 2009, 25(14): 1754-1760.

[20] Lidin OG. Tanegashima-the arrival of Europe in Japan. Routledge; 2003.

[21] LiuH, LiQ, XingY. Genes contributing to domestication of rice seed traits and its global expansion. Genes, 2018.

[22] LiuX, FuYX. Stairway plot 2: demographic history inference with folded SNP frequency spectra. Genome Biol, 2020, 211280.

[23] MalinskyM, MatschinerM, SvardalH. Dsuite-fast D-statistics and related admixture evidence from VCF files. Mol Ecol Resour, 2021, 21(2): 584-595.

[24] ManichaikulA, MychaleckyjJC, RichSS, DalyK, SaleM, ChenWM. Robust relationship inference in genome-wide association studies. Bioinformatics, 2010, 26(22): 2867-2873.

[25] MatsuzawaY, MekiyanonS, KanekoY, BangSW, WakuiK, TakahataY. Male sterility in alloplasmic Brassica rapa L. carrying Eruca sativa cytoplasm. Plant Breed, 1999, 118(1): 82-84

[26] Morales M, Janick J. Arugula: a promising specialty leaf vegetable. In: Janick J, Whipkey A, editors. n.d.

[27] NguyenLT, SchmidtHA, von HaeselerA, MinhBQ. IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol Biol Evol, 2015, 32(1): 268-74.

[28] Nothnagel T, Budahn H, Schrader O, Klocke E, editors. Intergeneric hybridization between Brassica oleracea and Eruca sativa. VI International Symposium on Brassicas and XVIII Crucifer Genetics Workshop 1005; 2012.

[29] Olesen MH, Romme HA, Kampmann HH. Rucola plants with cytoplasmic male sterility (CMS). US8247655B2. Google; 2012.

[30] ParkJS, ParkJH, ParkYD. Construction of pseudomolecule sequences of Brassica rapa ssp. pekinensis inbred line CT001 and analysis of spontaneous mutations derived via sexual propagation. PLoS One, 2019, 149. e0222283

[31] PickrellJK, PritchardJK. Inference of population splits and mixtures from genome-wide allele frequency data. PLoS Genet, 2012, 811. e1002967

[32] Pignone D, Gómez-Campo C. Eruca. Wild crop relatives: genomic and breeding resources. In: Kole C, editor. Springer; 2011. p. 149–60. https://doi.org/10.1007/978-3-642-21102-7.

[33] R Core Team. R: A language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing; 2023.

[34] Renfrew J, Sanderson H.Herbs and vegetables. In: Ghillean P, Mark N, editors. 2012.

[35] RStudio Team. RStudio: Integrated Development for R. PBC, Boston, MA: RStudio; 2020.

[36] SastryE. Taramira (Eruca sativa) and its improvement–a review. Agric Rev, 2003, 24(4): 235-249

[37] SihagRC. Pollination ecology of rocket (Eruca vesicaria (L.) Cav. ssp. sativa (Mill.) Thell) in the semi-arid environments of Northwest India: native bees are the major pollinators. Ecologies, 2023, 4(3): 580-594.

[38] SubramanianB, GaoS, LercherMJ, HuS, ChenWH. Evolview v3: a webserver for visualization, annotation, and management of phylogenetic trees. Nucleic Acids Res, 2019, 47(W1): W270-W275.

[39] SunW, GuanC, MengY, LiuZ, ZhangT, LiX, et al.. Intergeneric crosses between Eruca sativa Mill. and Brassica species. Zuo Wu Xue Bao, 2005, 1: 36-42

[40] SzklarczykD, GableAL, LyonD, JungeA, WyderS, Huerta-CepasJ, et al.. STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res, 2019, 47(D1): D607-D613.

[41] TraboulayDMColumbus and Las Casas: the conquest and Christianization of America, 1492–1566, 1994, Maryland. University Press of America.

[42] WangY, ChuP, YangQ, ChangS, ChenJ, HuM, et al.. Complete mitochondrial genome of Eruca sativa Mill. (Garden rocket). PLoS One, 2014, 98. e105748

[43] WarwickSI, GugelRK, Gómez-CampoC, JamesT. Genetic variation in Eruca vesicaria (L.) Cav. Plant Genet Resour, 2007, 5(3): 142-53.

[44] XuG, ZhangX, ChenW, ZhangR, LiZ, WenW, et al.. Population genomics of Zea species identifies selection signatures during maize domestication and adaptation. BMC Plant Biol, 2022, 22172.

[45] Yaniv Z, editor. Tradition, uses and research on rocket in Israel. In: Rocket: a Mediterranean crop for the world. Report of a workshop. Legnaro (Padova); 1996.

[46] Zafar-PashanezhadM, ShahbaziE, GolkarP, ShiranB. Genetic variation of Eruca sativa L. genotypes revealed by agro-morphological traits and ISSR molecular markers. Ind Crops Prod, 2020, 145. 111992

Funding

Singapore Food Agency(SFS_RND_SUFP_001_04)

RIGHTS & PERMISSIONS

The Author(s)