Integrative phylogenomics and morphology reveal the evolution and biogeography of Encephalartos (Zamiaceae)

Abstract

Encephalartos, an African endemic genus within the Zamiaceae, comprises 65 extant species whose phylogenetic relationships have remained unresolved due to limited genetic differentiation observed in previous studies. This research reconstructs the evolutionary history of Encephalartos utilizing 3545 single-copy nuclear genes derived from transcriptomes of 64 species. The study estimates divergence times and reconstructs ancestral states for 12 key morphological traits. Phylogenetic analyses definitively resolve eight major clades, supported by both molecular and morphological evidence. Although these clades partially align with previous morphology and geography based classifications, the genomic data provide novel insights, necessitating a revised infrageneric system. Biogeographic reconstructions indicate that Encephalartos originated in southern Africa during the Oligocene (~26.3 Ma), subsequently dispersing into eastern and northern Africa through the Zimbabwe–Mozambique corridor during the Miocene, followed by expansion into Central Africa. Speciation rates decreased markedly during the Pliocene and Pleistocene, potentially due to intensified climatic drying and cooling. Morphological character mapping identified ancestral traits, including aerial stems, green leaves, and red sarcotesta. Specific transitions such as subterranean stems in clade IV and bluish-green leaves in clades II and V, further substantiate clade differentiation. These findings resolve long-standing taxonomic uncertainties and emphasize the Oligocene-Miocene as a crucial period for Encephalartos diversification, influenced by Cenozoic climate change. This research establishes a robust framework for future systematic and conservation studies while demonstrating the effectiveness of transcriptome data in resolving phylogenies of slowly evolving lineages.

Keywords

biogeography / cycads / Encephalartos / phylogenetics / transcriptome data / Zamiaceae

References

[1] Abascal F, Zardoya R, Telford MJ. 2010. TranslatorX: Multiple alignment of nucleotide sequences guided by amino acid translations. Nucleic Acids Research 38: W7–W13.

[2] Balkenhol N, Dudaniec RY, Krutovsky KV, Johnson JS, Cairns DM, Segelbacher G, Selkoe KA, von der Heyden S, Wang IJ, Selmoni O, Joost S. 2019. Landscape genomics: Understanding relationships between environmental heterogeneity and genomic characteristics of populations. In: Rajora OP ed. Population genomics: Concepts, approaches and applications. Cham: Springer International Publishing. 261–322.

[3] Bezuidenhout H. 2020. Assessment results (2015–2018) of re-established poached cycad trees in Addo Elephant National Park, Eastern Cape, South Africa. Koedoe 62: a1626.

[4] Bolger AM, Lohse M, Usadel B. 2014. Trimmomatic: A flexible trimmer for illumina sequence data. Bioinformatics 30: 2114–2120.

[5] Calonje M, Meerow AW, Griffith MP, Salas-Leiva D, Vovides AP, Coiro M, Francisco-Ortega J. 2019. A time-calibrated species tree phylogeny of the New World cycad genus Zamia L. (Zamiaceae, Cycadales). International Journal of Plant Sciences 180: 286–314.

[6] Calonje M, Stevenson DW, Osborne R. 2025. The World List of Cycads (Version 2025.09.04-r1). Coral Gables, FL: Montgomery Botanical Center. https://doi.org/10.5281/zenodo.17059408

[7] Chaiprasongsuk M, Mingmuang M, Thongpan A, Namwongprom K. 2007. Molecular identification of Encephalartos (Zamiaceae) species and their relationships to morphological characters. Agriculture and Natural Resources 41: 43–60.

[8] Clugston JAR, Griffith MP, Kenicer GJ, Husby CE, Calonje MA, Stevenson DW, Little DP. 2016. Zamia (Zamiaceae) phenology in a phylogenetic context: Does in situ reproductive timing correlate with ancestry? Edinburgh Journal of Botany 73: 345–370.

[9] Clugston JAR, Ruhsam M, Kernicher GJ, Henwood M, Milne RI, Nagalingum N. 2022. Conservation genomics of an Australian cycad Cycas calcicola, and the absence of key genotypes in botanic gardens. Conservation Genetics 23: 449–465.

[10] Coetzee J. 1993. African flora since the terminal Jurassic. Biological relationships between Africa and South America 37: 37–61.

[11] Coiro M, Allio R, Mazet N, Seyfullah LJ, Condamine FL. 2023. Reconciling fossils with phylogenies reveals the origin and macroevolutionary processes explaining the global cycad biodiversity. New Phytologist 240: 1616–1635.

[12] Donaldson JS, Hill KD, Stevenson DW. 2003. Cycads of the World: An Overview. In: JS Donaldson ed. Cycads: Status survey and conservation action plan. Gland, Switzerland and Cambridge, UK: IUCN. 3–8.

[13] Emms DM, Kelly S. 2019. Orthofinder: Phylogenetic orthology inference for comparative genomics. Genome Biology 20: 1–14.

[14] Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, Adiconis X, Fan L, Raychowdhury R, Zeng Q. 2011. Full-length transcriptome assembly from RNA-seq data without a reference genome. Nature Biotechnology 29: 644–652.

[15] Griffith MP, Clase T, Toribio P, Encarnación Piñeyro Y, Jiménez F, Gratacos X, Sanchez V, Meerow A, Meyer A, Kramer A, Fant J, Havens K, Magellan TM, Dosmann M, Hoban S. 2020. Can a botanic garden metacollection better conserve wild plant diversity? A case study comparing pooled collections with an ideal sampling model. International Journal of Plant Sciences 181: 485–496.

[16] Grobbelaar N. 2004. Cycads with special reference to the Southern African species. South Africa: N. Grobbelaar, Pretoria.

[17] Gutiérrez-Ortega JS, Salinas-Rodríguez MM, Martínez JF, Molina-Freaner F, Pérez-Farrera MA, Vovides AP, Matsuki Y, Suyama Y, Ohsawa TA, Watano Y, Kajita T. 2018. The phylogeography of the cycad genus Dioon (Zamiaceae) clarifies its Cenozoic expansion and diversification in the Mexican transition zone. Annals of Botany 121: 535–548.

[18] Habib S, Gong Y, Dong S, Lindstrom A, Stevenson DW, Wu H, Zhang S. 2023. Phylotranscriptomics shed light on intrageneric relationships and historical biogeography of Ceratozamia (Cycadales). Plants 12: 478.

[19] Habib S, Gong Y, Dong S, Lindstrom A, William Stevenson D, Liu Y, Wu H, Zhang S. 2022. Phylotranscriptomics reveal the spatio-temporal distribution and morphological evolution of Macrozamia, an australian endemic genus of Cycadales. Annals of Botany 130: 671–685.

[20] Helfrich P, Rieb E, Abrami G, Lücking A, Mehler A. 2018. TreeAnnotator: Versatile Visual Annotation of Hierarchical Text Relations. In: N Calzolari, K Choukri, C Cieri, T Declerck, S Goggi, K Hasida, Isahara H, Maegaard B, Mariani J, Mazo H, Moreno A, Odijk J, Piperidis S, Tokunaga T eds. Proceedings of the Eleventh International Conference on Language Resources and Evaluation (LREC 2018). Japan: European Language Resources Association (ELRA).

[21] Hurter P, Glen H. 1996. Encephalartos hirsutus (Zamiaceae): A newly described species from South Africa. South African Journal of Botany 62: 46–48.

[22] Huson DH, Bryant D. 2006. Application of phylogenetic networks in evolutionary studies. Molecular Biology and Evolution 23: 254–267.

[23] IUCN. 2025. The IUCN red list of threatened species. Version 2025-1 [online]. Available from https://www.iucnredlist.org. [accessed 05 May 2025].

[24] Jacobs BF. 2004. Palaeobotanical studies from Tropical Africa: Relevance to the evolution of forest, woodland and savannah biomes. Philosophical Transactions of the Royal Society of London. Series B: Biological Sciences 359: 1573–1583.

[25] Kalyaanamoorthy S, Minh BQ, Wong TK, Von Haeseler A, Jermiin LS. 2017. ModelFinder: Fast model selection for accurate phylogenetic estimates. Nature Methods 14: 587–589.

[26] Katoh K, Standley DM. 2013. MAFFT multiple sequence alignment software version 7: Improvements in performance and usability. Molecular Biology and Evolution 30: 772–780.

[27] Keppel G. 2009. Morphological variation, an expanded description and ethnobotanical evaluation of Cycas seemannii A.Braun (Cycadaceae). The South Pacific Journal of Natural Science 27: 20–27.

[28] Konings KM. 2006. Life history traits of South African Encephalartos spp. (Zamiaceae) and their implications for understanding population structure, responses to threats and effective conservation action. Masters Dissertation. South Africa: University of Cape Town.

[29] Kück P, Longo GC. 2014. FASconCAT-G: Extensive functions for multiple sequence alignment preparations concerning phylogenetic studies. Frontiers in Zoology 11: 81.

[30] Laetsch DR, Blaxter ML. 2017. Kinfin: Software for taxon-aware analysis of clustered protein sequences. G3: Genes, Genomes, Genetics 7: 3349–3357.

[31] Landis MJ, Matzke NJ, Moore BR, Huelsenbeck JP. 2013. Bayesian analysis of biogeography when the number of areas is large. Systematic Biology 62: 789–804.

[32] Lemmon EM, Lemmon AR. 2013. High-throughput genomic data in systematics and phylogenetics. Annual Review of Ecology, Evolution, and Systematics 44: 99–121.

[33] Lindstrom A, Habib S, Dong S, Gong Y, Liu J, Calonje M, Stevenson D, Zhang S. 2024. Transcriptome sequencing data provide a solid base to understand the phylogenetic relationships, biogeography and reticulated evolution of the genus Zamia L. (Cycadales: Zamiaceae). Annals of Botany 134: 747–768.

[34] Liu J, Lindstrom A, Gong Y, Dong S, Liu Y, Zhang S, Xun G. 2024. Eco-evolutionary evidence for the global diversity pattern of Cycas (Cycadaceae). Journal of Integrative Plant Biology 66: 1170–1191.

[35] Liu J, Zhang S, Nagalingum N, Chiang YC, Lindstrom A, Xun G. 2018. Phylogeny of the gymnosperm genus Cycas L. (Cycadaceae) as inferred from plastid and nuclear loci based on a large-scale sampling: Evolutionary relationships and taxonomical implications. Molecular Phylogenetics and Evolution 127: 87–97.

[36] Liu Y, Wang S, Li L, Yang T, Dong S, Wei T, Wu S, Liu Y, Gong Y, Feng X, Ma J, Chang G, Huang J, Yang Y, Wang H, Liu M, Xu Y, Liang H, Yu J, Cai Y, Zhang Z, Fan Y, Mu W, Sahu SK, Liu S, Lang X, Yang L, Li N, Habib S, Yang Y, Lindstrom AJ, Liang P, Goffinet B, Zaman S, Wegrzyn JL, Li D, Liu J, Cui J, Sonnenschein EC, Wang X, Ruan J, Xue J-Y, Shao Z-Q, Song C, Fan G, Li Z, Zhang L, Liu J, Liu Z-J, Jiao Y, Wang X-Q, Wu H, Wang E, Lisby M, Yang H, Wang J, Liu X, Xu X, Li N, Soltis PS, Van de Peer Y, Soltis DE, Gong X, Liu H, Zhang S. 2022. The Cycas genome and the early evolution of seed plants. Nature Plants 8: 389–401.

[37] Mabunda MA. 2007. Species-level phylogenetic reconstruction of the African cycad genusEncephalartos (Zamiaceae). Masters Dissertation. South Africa: University of the Western Cape.

[38] Maddison WP, Maddison DR. 2023. Mesquite: A modular system for evolutionary analysis. Version 3.81 [online]. Available from http://www.mesquiteproject.org [accessed 15 September 2024].

[39] Maley J. 1996. The African rain forest–main characteristics of changes in vegetation and climate from the upper cretaceous to the quaternary. Proceedings of the Royal Society of Edinburgh, Section B: Biological Sciences 104: 31–73.

[40] Mankga LT, Yessoufou K, Chitakira M. 2020. On the origin and diversification history of the African genus Encephalartos. South African Journal of Botany 130: 231–239.

[41] Matzke N. 2018. BioGeoBEARS: BioGeography with bayesian (and likelihood) evolutionary analysis with R scripts. version 1.1. 1. San Francisco, CA: GitHub. Available online at: https://github.com/nmatzke/BioGeoBEARS.

[42] McCormack JE, Hird SM, Zellmer AJ, Carstens BC, Brumfield RT. 2013. Applications of next-generation sequencing to phylogeography and phylogenetics. Molecular Phylogenetics and Evolution 66: 526–538.

[43] Minh BQ, Schmidt HA, Chernomor O, Schrempf D, Woodhams MD, Von Haeseler A, Lanfear R. 2020. IQ-TREE 2: New models and efficient methods for phylogenetic inference in the genomic era. Molecular Biology and Evolution 37: 1530–1534.

[44] Medina-Villarreal A, González-Astorga J, de Los Monteros AE. 2019. Evolution of Ceratozamia cycads: A proximate-ultimate approach. Molecular Phylogenetics and Evolution 139: 106530.

[45] Mo YK, Lanfear R, Hahn MW, Minh BQ. 2022. Updated site concordance factors minimize effects of homoplasy and taxon sampling. Bioinformatics 39: btac741.

[46] Morley RJ. 2001. Why are there so many primitive angiosperms in the rain forests of the far East? In: I Metcalfe, JMB Smith, M Morwood, I Davidson eds. Faunal and floral migrations and evolution in SE Asia-Australia. The Netherlands: A. A. Balkema, Rotterdam. 185–199.

[47] Oliver JC. 2013. Microevolutionary processes generate phylogenomic discordance at ancient divergences. Evolution 67: 1823–1830.

[48] Osborne R, Grobbelaar N, Vincent P. 1993. A numerical phenetic study of the genus Encephalartos Lehm. The biology, structure and systematics of the Cycadales. In: DW Stevenson, KJ Norstog eds. Proceedings of CYCAD 90, The Second International Conference on Cycad Biology. Milton: Palm and Cycad Societies of Australia. 279–280.

[49] Plana V. 2004. Mechanisms and tempo of evolution in the African Guineo–Congolian rainforest. Philosophical Transactions of the Royal Society of London. Series B: Biological Sciences 359: 1585–1594.

[50] Prakash S, Grobbelaar N, Van Staden J. 2008. Diversity in Encephalartos woodii collections based on Random Amplified DNA markers (RAPD's) and Inter-Specific Sequence Repeats (ISSR's). South African Journal of Botany 74: 341–344.

[51] Ree RH, Smith SA. 2008. Maximum likelihood inference of geographic range evolution by dispersal, local extinction, and cladogenesis. Systematic Biology 57: 4–14.

[52] Ronquist F. 1997. Dispersal-vicariance analysis: A new approach to the quantification of historical biogeography. Systematic Biology 46: 195–203.

[53] Rousseau P. 2012. A molecular systematic study of the African endemic cycads. Masters Dissertation. South Africa University of Johannesburg.

[54] Sanderson MJ. 2002. Estimating absolute rates of molecular evolution and divergence times: A penalized likelihood approach. Molecular Biology and Evolution 19: 101–109.

[55] Shen W, Le S, Li Y, Hu F. 2016. SeqKit: a cross-platform and ultrafast toolkit for FASTA/Q file manipulation. PLoS One 11: e0163962.

[56] Smith SA, Brown JW, Walker JF. 2018. So many genes, so little time: A practical approach to divergence-time estimation in the genomic era. PLoS One 13: e0197433.

[57] Smith SA, O'Meara BC. 2012. TreePL: Divergence time estimation using penalized likelihood for large phylogenies. Bioinformatics 28: 2689–2690.

[58] Stamatakis A. 2014. RAxML version 8: A tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics 30: 1312–1313.

[59] Stewart R, Clugston JAR, Williamson J, Niemann HJ, Little DP, van der Bank M. 2023. Species relationships and phylogenetic diversity of the African genus Encephalartos Lehm. (Zamiaceae). South African Journal of Botany 152: 165–173.

[60] Talavera G, Castresana J. 2007. Improvement of phylogenies after removing divergent and ambiguously aligned blocks from protein sequence alignments. Systematic Biology 56: 564–577.

[61] Tao Y, Chen B, Kang M, Liu Y, Wang J. 2021. Genome-wide evidence for complex hybridization and demographic history in a group of Cycas from China. Frontiers in Genetics 12: 717200.

[62] Than C, Ruths D, Nakhleh L. 2008. PhyloNet: A software package for analyzing and reconstructing reticulate evolutionary relationships. BMC Bioinformatics 9: 1–16.

[63] Treutlein J, Vorster P, Wink M. 2005. Molecular relationships in Encephalartos (Zamiaceae, Cycadales) based on nucleotide sequences of nuclear ITS 1&2, rbcl, and genomic ISSR fingerprinting. Plant Biology 7: 79–90.

[64] Van der Bank H, Wink M, Vorster P, Treutlein J, Brand L, van der Bank M, Hurter J. 2001. Allozyme and DNA sequence comparisons of nine species of Encephalartos (Zamiaceae). Biochemical Systematics and Ecology 29: 241–266.

[65] Viljoen CD, van Staden J. 2006. The genetic relationship between Encephalartos natalensis and E. woodii determined using RAPD fingerprinting. South African Journal of Botany 72: 642–645.

[66] Vorster P. 1993. Taxonomy of Encephalartos (Zamiaceae): Taxonomically useful external characteristics. In: DW Stevenson, KJ Norstog eds. Proceedings of CYCAD 90, The Second International Conference on Cycad Biology. Queensland: Palm and Cycad Societies of Australia. 294–299.

[67] Vorster P. 2004. Classification concepts in Encephalartos (Zamiaceae). In: T Walters, R Osborne eds. Cycad classification: Concepts and recommendations. UK: CABI Publishing Wallingford. 69–83.

[68] Vorster P, Van Der Bank F, Van der Bank M, Wink M. 2004. Phylogeny of Encephalartos: Some Eastern Cape species. The Botanical Review 70: 250–259.

[69] Walters T, Osborne R. 2004. Cycad classification: Concepts and recommendations. UK: CABI Publishing, Wallingford.

[70] Whitfield JB, Lockhart PJ. 2007. Deciphering ancient rapid radiations. Trends in Ecology & Evolution 22: 258–265.

[71] Wichura H, Bousquet R, Oberhänsli R, Strecker MR, Trauth MH. 2010. Evidence for Middle Miocene uplift of the East African Plateau. Geology 38: 543–546.

[72] Yessoufou K, Bamigboye SO, Daru BH, van der Bank M. 2014. Evidence of constant diversification punctuated by a mass extinction in the African cycads. Ecology and Evolution 4: 50–58.

[73] Yu Y, Blair C, He X. 2020. RASP 4: Ancestral state reconstruction tool for multiple genes and characters. Molecular Biology and Evolution 37: 604–606.

[74] Zachos J, Pagani M, Sloan L, Thomas E, Billups K. 2001. Trends, rhythms, and aberrations in global climate 65 MA to present. Science 292: 686–693.

[75] Zachos JC, Dickens GR, Zeebe RE. 2008. An early cenozoic perspective on greenhouse warming and carbon-cycle dynamics. Nature 451: 279–283.

[76] Zhang C, Mirarab S. 2022. Weighting by gene tree uncertainty improves accuracy of quartet-based species trees. Molecular Biology and Evolution 39: msac215.

RIGHTS & PERMISSIONS

2025 The Author(s). Journal of Systematics and Evolution published by John Wiley & Sons Australia, Ltd on behalf of Institute of Botany, Chinese Academy of Sciences.