Accumulation of mutations in genes associated with sexual reproduction contributed to the domestication of a vegetatively propagated staple crop, enset

Abstract

Enset (Ensete ventricosum (Welw.) Cheesman) is a drought tolerant, vegetatively propagated crop that was domesticated in Ethiopia. It is a staple food for more than 20 million people in Ethiopia. Despite its current importance and immense potential, enset is among the most genetically understudied and underexploited food crops. We collected 230 enset wild and cultivated accessions across the main enset producing regions in Ethiopia and applied amplified fragment length polymorphism (AFLP) and genotype by sequencing (GBS) analyses to these accessions. Wild and cultivated accessions were clearly separated from each other, with 89 genes found to harbour SNPs that separated wild from cultivated accessions. Among these, 17 genes are thought to be involved in flower initiation and seed development. Among cultivated accessions, differentiation was mostly associated with geographical location and with proximity to wild populations. Our results indicate that vegetative propagation of elite clones has favoured capacity for vegetative growth at the expense of capacity for sexual reproduction. This is consistent with previous reports that cultivated enset tends to produce non-viable seeds and flowers less frequently than wild enset.

References

[1] Martínez-Ainsworth, N. E. & Tenaillon, M. I. Superheroes and masterminds of plant domestication. C. R. Biol. 339, 268-273 (2016).

[2] Denham, T. et al. The domestication syndrome in vegetatively propagated field crops. Ann. Bot. 125, 581-597 (2020).

[3] Cheesman, E. Classification of the Bananas: The Genus Musa L. Kew Bull. 2, 106-117 (1947).

[4] Tamrat, S. et al. Germination ecology of wild and domesticated Ensete ventricosum: evidence for maintenance of sexual reproductive capacity in a vegetatively propagated perennial crop. Preprint at https://doi.org/10.1101/2020.04.30.055582 (2020).

[5] Borrell, J. S. et al. Enset in Ethiopia: a poorly characterized but resilient starch staple. Ann. Bot. 123, 747-766 (2019).

[6] Eshetae, M. A., Hailu, B. T. & Demissew, S. Spatial characterization and distribution modelling of Ensete ventricosum (wild and cultivated) in Ethiopia. Geocarto Int. 1-16 (2019).

[7] Guzzon, F. & Müller, J. V. Current availability of seed material of enset (Ensete ventricosum, Musaceae) and its Sub-Saharan wild relatives. Genet. Resour. Crop Evol. 63, 185-191 (2016).

[8] Harrison, J. et al. A draft genome sequence for Ensete ventricosum, the drought-tolerant “tree against hunger”. Agronomy 4, 13-33 (2014).

[9] Brandt, S. A. et al. The tree against hunger. Enset-based agricultural systems in Ethiopia American Association for the Advancement of Science, Washington DC (1997).

[10] Quinlan, M. B., Quinlan, R. J. & Dira, S. Sidama agro-pastoralism and ethnobiological classification of its primary plant, Enset (Ensete ventricosum). Ethnobiol. Lett. 5, 116-125 (2014).

[11] Negash, A., Tsegaye, A., van Treuren, R. & Visser, B. AFLP analysis of enset clonal diversity in south and southwestern Ethiopia for conservation. Crop Sci. 42, 1105-1111 (2002).

[12] Birmeta, G., Nybom, H. & Bekele, E. Distinction between wild and cultivated enset (Ensete ventricosum) gene pools in Ethiopia using RAPD markers. Hereditas 140, 139-148 (2004).

[13] Tobiaw, D. C. & Bekele, E. Analysis of genetic diversity among cultivated enset (Ensete ventricosum) populations from Essera and Kefficho, southwestern part of Ethiopia using inter simple sequence repeats (ISSRs) marker. Afr. J. Biotechnol. 10, 15697-15709 (2013).

[14] Gerura, F. N. et al. Genetic diversity and population structure of enset (Ensete ventricosum Welw Cheesman) landraces of Gurage zone, Ethiopia. Genet. Resour. Crop Evol. 66, 1813-1824 (2019).

[15] McKey, D., Elias, M., Pujol, B. & Duputié, A. The evolutionary ecology of clonally propagated domesticated plants. New Phytol. 186, 318-332 (2010).

[16] Olango, T. M., Tesfaye, B., Pagnotta, M. A., Pè, M. E. & Catellani, M. Development of SSR markers and genetic diversity analysis in enset (Ensete ventricosum (Welw.) Cheesman), an orphan food security crop from Southern Ethiopia. BMC Genet. 16, 98 (2015).

[17] Kamanda, I. et al. Genetic diversity of provitamin-A cassava (Manihot esculenta Crantz) in Sierra Leone. Genet. Resour. Crop Evol. 67, 1-16 (2020).

[18] Mandel, J., Dechaine, J., Marek, L. & Burke, J. Genetic diversity and population structure in cultivated sunflower and a comparison to its wild progenitor, Helianthus annuus L. Theor. Appl. Genet. 123, 693-704 (2011).

[19] Becerra, V., Paredes, M., Ferreira, M. E., Gutiérrez, E. & Díaz, L. M. Assessment of the genetic diversity and population structure in temperate japonica rice germplasm used in breeding in Chile, with SSR markers. Chil. J. Agric. Res. 77, 15-26 (2017).

[20] Yang, H. et al. Genetic divergence between Camellia sinensis and its wild relatives revealed via genome-wide SNPs from RAD sequencing. PLoS ONE 11, e0151424 (2016).

[21] Zhao, D.-W., Yang, J.-B., Yang, S.-X., Kato, K. & Luo, J.-P. Genetic diversity and domestication origin of tea plant Camellia taliensis (Theaceae) as revealed by microsatellite markers. BMC Plant Biol. 14, 14 (2014).

[22] Birmeta, G. Genetic Variability and Biotechnological Studies for the Conservation and Improvement of Ensete ventricosum, Vol. 502 (Sveriges lantbruksuniv., Acta Universitatis Agriculturae Sueciae, 2004).

[23] Simko, I., Haynes, K. G. & Jones, R. W. Assessment of linkage disequilibrium in potato genome with single nucleotide polymorphism markers. Genetics 173, 2237-2245 (2006).

[24] Nascimento, G., Savegnago, R., Ventura, R., Ledur, M. & Munari, D. In Embrapa Suínos e Aves-Artigo em anais de congresso (ALICE) (ed. Goddard, M.) (In: World Congress of Genetics Applied to Livestock Production, 10., American Society of Animal Science, 2014 …).

[25] Birmeta, G. Biotechnological studies on enset (Ensete ventricosum (Welw.) Cheesman)), a food security staple food crop of Ethiopia. Ethiop. J. Niol. Sci. 17, 75-101 (2018).

[26] Burgos-Hernández, M., Hernández, D. G. & Castillo-Campos, G. Genetic diversity and population genetic structure of wild banana Musa ornata (Musaceae) in Mexico. Plant Syst. Evol. 299, 1899-1910 (2013).

[27] Hildebrand, E. in Biodiversity Research in the Horn of Africa Region, Vol. 54 (eds Friis, I. & Ryding, O) 287-309 (Royal Danish Academy of Sciences and Letters, 2001).

[28] Lotterhos, K. E. & Whitlock, M. C. Evaluation of demographic history and neutral parameterization on the performance of FST outlier tests. Mol. Ecol. 23, 2178-2192 (2014).

[29] Hildebrand, E. in Biodiversity Research in the Horn of Africa Region: Proceedings of the Third International Symposium on the Flora of Ethiopia and Eritrea at the Carlsberg Academy, Copenhagen, August 25-27, 1999, 287 (Kgl. Danske Videnskabernes Selskab).

[30] Borrel, J. S. et al. Enset in Ethiopia: a poorly characterised but resilient starch staple. Ann. Bot. 123, 747-766, https://doi.org/10.1093/aob/mcy214 (2019).

[31] Vos, P. et al. AFLP: a new technique for DNA fingerprinting. Nucleic Acids Res. 23, 4407-4414 (1995).

[32] López, C. M. R. et al. Detection and quantification of tissue of origin in salmon and veal products using methylation sensitive AFLPs. Food Chem. 131, 1493-1498 (2012).

[33] Peakall, R. & Smouse, P. E. GenAlEx 6.5: genetic analysis in Excel. Population genetic software for teaching and research-an update. Bioinformatics 28, 2537-2539 (2012).

[34] Xie, H. et al. Global DNA methylation patterns can play a role in defining terroir in grapevine (Vitis vinifera cv. Shiraz). Front. Plant Sci. 8, 1860 (2017).

[35] Konate, M. et al. Salt stress induces non-CG methylation in coding regions of barley seedlings (Hordeum vulgare). Epigenomes 2, 12 (2018).

[36] Lu, F. et al. Switchgrass genomic diversity, ploidy, and evolution: novel insights from a network-based SNP discovery protocol. PLoS Genet. 9, e1003215 (2013).

[37] Torkamaneh, D., Laroche, J. & Belzile, F. Genome-wide SNP calling from genotyping by sequencing (GBS) data: a comparison of seven pipelines and two sequencing technologies. PloS ONE 11, e0161333 (2016).

[38] D’hont, A. et al. The banana (Musa acuminata) genome and the evolution of monocotyledonous plants. Nature 488, 213 (2012).

[39] Pfeifer, B., Wittelsbürger, U., Ramos-Onsins, S. E. & Lercher, M. J. PopGenome: an efficient Swiss army knife for population genomic analyses in R. Mol. Biol. Evol. 31, 1929-1936 (2014).

[40] Liu, K. & Muse, S. V. PowerMarker: an integrated analysis environment for genetic marker analysis. Bioinformatics 21, 2128-2129 (2005).

[41] Bradbury, P. J. et al. TASSEL: software for association mapping of complex traits in diverse samples. Bioinformatics 23, 2633-2635 (2007).

[42] Parks, D. H. et al. GenGIS: A geospatial information system for genomic data. Genome Res. 19, gr. 095612.095109 (2009).

[43] Jombart, T., Devillard, S. & Balloux, F. Discriminant analysis of principal components: a new method for the analysis of genetically structured populations. BMC Genet. 11, 94 (2010).

[44] Pritchard, J. K., Stephens, M. & Donnelly, P. Inference of population structure using multilocus genotype data. Genetics 155, 945-959 (2000).

[45] Purcell, S. et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am. J. Hum. Genet. 81, 559-575 (2007).

[46] Krzywinski, M. I. et al. Circos: an information aesthetic for comparative genomics. Genome Res. 19, (2009).

[47] Beaumont, M. A. & Nichols, R. A. Evaluating loci for use in the genetic analysis of population structure. Proc. R. Soc. Lond. Ser. B 263, 1619-1626 (1996).

[48] Antao, T., Lopes, A., Lopes, R. J., Beja-Pereira, A. & Luikart, G. LOSITAN: a workbench to detect molecular adaptation based on a F ST-outlier method. BMC Bioinforma. 9, 323 (2008).

[49] Bache, S. M. & Wickham, H. magrittr: A Forward-pipe Operator for R. R package version 1.5, https://CRAN.R-project.org/package=magrittr (2014).

[50] Luu, K., Bazin, E. & Blum, M. G. pcadapt: an R package to perform genome scans for selection based on principal component analysis. Mol. Ecol. Resour. 17, 67-77 (2017).