Accessibility / Report Error

Reconstruction of the evolutionary biogeography reveals the origin and diversification of Arisaema (Araceae)

ABSTRACT

Explaining where plants originated and how plants migrated between eastern Asia and other Northern Hemisphere regions is still challenging. The genus Arisaema is considered a good example to illuminate this process, due to its wide distribution in both the tropical and temperate zones. Based on comprehensive taxon sampling, the historical biogeography of Arisaema was conducted using the molecular data from four chloroplast DNA regions. The phylogeny of Arisaema was studied by using ML and BI methods. Divergence times of Arisaema species were estimated by the uncorrelated lognormal Bayesian method. The ancestral area reconstruction of Arisaema was performed with BioGeoBEARS and Bayes-DIVA. Our results indicated that the crown group of Arisaema originated in Southeast Asia ca. 31.28 Ma. After the initial origin, Arisaema was inferred to have become widespread in Southeast Asia. Members in Himalayas, Africa and North China, Japan, the Korean Peninsula region (NJK region) were derived from migrants originating in Southeast Asia. American members were derived from migrants originating in the NJK region. Arisaema originated in Southeast Asia during the Oligocene with the development of tropical rainforest and a warm, moist climate. During the Oligocene to Miocene, Southeast Asia was the original center for the diversification of Arisaema in the NJK region.

Keywords:
Arisaema; biogeography; origin; Eocene; Southeast Asia

Introduction

The genus Arisaema is one of the largest genera in the family Araceae in the World with about 200 species (Ohi-Toma et al. 2016Ohi-Toma T, Wu S, Murata H, Murata J. 2016. An updated genus-wide phylogenetic analysis of Arisaema (Araceae) with reference to sections. Botanical Journal of the Linnean Society 182: 100-114.). It is distributed in the temperate and warm-temperate areas from the East to South-East of North America to northern Mexico, from tropical East and North-East Africa to the Arabian Peninsula, and from temperate East Asia to South and South-East Asia and Malesia (Govaerts et al. 2002Govaerts R, Frodin DG, Bogner J, et al. 2002. Arisaema. World checklist and Bibliography of Araceae (and Acoraceae). Royal Botanic Gardens, Kew, Richmond. p. 181-218.; Gusman & Gusman 2006Gusman G, Gusman L. 2006. The genus Arisaema: a monograph for botanists and nature lovers, second revised and enlarged edition. Ruggell, Ganter Verlag.; Murata 2011Murata J. 2011. The picture book of plant systematics in color: Arisaema in Japan. Hokuryukan, Tokyo. (in Japanese with taxonomic keys in English); Ohi-Toma et al. 2016Ohi-Toma T, Wu S, Murata H, Murata J. 2016. An updated genus-wide phylogenetic analysis of Arisaema (Araceae) with reference to sections. Botanical Journal of the Linnean Society 182: 100-114.). In the traditional medicine of India, China and Vietnam, some species of Arisaema are used as local medicine (Li et al. 2010Li H, Zhu G, Murata J. 2010. Arisaema (Araceae). In: Wu ZY, Raven PH, Hong DY. (eds.) Flora of China. Vol. 23. Beijing and St. Louis, Science Press and Missouri Botanical Garden Press. p. 43-69.). The young leaves and inflorescences of some species are used as vegetable or food plants (Gusman & Gusman 2006Gusman G, Gusman L. 2006. The genus Arisaema: a monograph for botanists and nature lovers, second revised and enlarged edition. Ruggell, Ganter Verlag.).

In the family Araceae, the genus Arisaema was placed in the tribe Arisaematae, subfamily Aroideae (Mayo et al. 1997Mayo SJ, Bogner J, Boyce PC. 1997. The genera of Araceae. Kew, Richmond, Royal Botanic Gardens.; Renner & Zhang 2004Renner SS, Zhang LB. 2004. Biogeography of the Pistia clade (Araceae): based on chloroplast and mitochondrial DNA sequences and Bayesian divergence time inference. Systematic Biology 53: 422-432. ). Arisaema is characterised by being seasonal or evergreen, monosexual or bisexual plants; having tuberous or rhizomatous stem, with or without pseudostem; leaf blades diverse from trifoliolate, pedate to radiate. Inflorescence solitary, bearing spadix monosexual or bisexual with the appendix morphology very diverse; male flowers with 2-5 stamens; female flowers with unilocular ovary containing several basal ovules (Fig. 1).

Figure 1
Diversity of morphological characters in the genus Arisaema. A - B: Diversity of habits and leaf blades; D - I: Diversity of inflorescences; J - O: Diversity of spadixes. A. A. menglaense; B. Arisaema sp.; C. A. rhizomatum; D. A. pierreanum; E. A. lidaense; F. A. averyanovii; G. A. franchetianum; H. A. rhizomatum; I. A. langbianense; J & K. A. menglaense; L. A. pingbianense; M. A. rhizomatum; N. A. rostratum; O. A. pierreanum (Photograph: V. D. Nguyen).

Since being established in 1831, there have been a wide variety of studies published within the genus Arisaema and several systems of classification proposed by Schott (1860Schott HW. 1860. Prodromus systematis aroidearum. Typis congregationis. Vienna, Mechitharisticae.), Engler (1920Engler A. 1920. Araceae-Aroideae, Araceae-Pistioidea. In: Engler A. (ed.) Das Pflanzenreich. Vol. 73. Leipzig, Verlag von Wilhelm Engelmann. p. 149-220.), Nakai (1950Nakai T. 1950. Classes, ordines, familiae, subfamiliae, tribus, genera nova quae attient ad plantas koreanas (Supplementum). Journal Japanese Botany 25: 5-7.), Hara (1971Hara H. 1971. A Revision of the Eastern Himalaya species of the genus Arisaema (Araceae). In: Flora of Eastern Himalaya, Second report. Tokyo, University of Tokyo Press. p. 321-354.), Murata (1984Murata J. 1984. An attempt at an infrageneric classification of genus Arisaema (Araceae). Journal of the Faculty of Science, University of Tokyo, Section III. Botany 13: 431-482., 1991Murata J. 1991. The systematic position of Arisaema nepenthoides and A. wattii (Araceae). Kew Bulletin 46: 119-128.), Gusman & Gusman (2002Gusman G, Gusman L. 2002. The genus Arisaema. A monograph for botanists and nature lovers. Liechtenstein, Ganter Verlag., 2006Gusman G, Gusman L. 2006. The genus Arisaema: a monograph for botanists and nature lovers, second revised and enlarged edition. Ruggell, Ganter Verlag.), Li et al. (2010Li H, Zhu G, Murata J. 2010. Arisaema (Araceae). In: Wu ZY, Raven PH, Hong DY. (eds.) Flora of China. Vol. 23. Beijing and St. Louis, Science Press and Missouri Botanical Garden Press. p. 43-69.), Murata (2011)Murata J. 2011. The picture book of plant systematics in color: Arisaema in Japan. Hokuryukan, Tokyo. (in Japanese with taxonomic keys in English), Murata et al. (2013)Murata J, Nagamasu H, Ohashi H. 2013. A nomenclatural review on the infrageneric classification of Arisaema (Araceae). Journal Japanese Botany 88: 36-45. and Ohi-Toma et al. (2016Ohi-Toma T, Wu S, Murata H, Murata J. 2016. An updated genus-wide phylogenetic analysis of Arisaema (Araceae) with reference to sections. Botanical Journal of the Linnean Society 182: 100-114.). All of the above studies focused on systematics and have provided a foundation for investigating the phylogeny and historical biogeography of the genus Arisaema. Renner et al. (2004Renner SS, Zhang LB, Murata J. 2004. A chloroplast phylogeny of Arisaema (Araceae) illustrates Tertiary floristic links between Asia, North America, and east Africa. American Journal of Botany 91: 881-888.) conducted a phylogenetic study based on 81 individuals of Arisaema. Results of this study supported the monophyly of this genus, but several polytomies were found in the phylogenetic trees of Arisaema. Additionally, the study also estimated the divergence time of the genus, provided the main biogeographical event of Arisaema, and considered floristic links among East Africa, Arabia, the Himalayan region, China, and North America. However, this study used only one calibration point (the fossil similar to A. triphyllum at 18 million years ago) and penalized likelihood. The resulting time estimates are not sufficiently reliable to clarify divergence time of Arisaema, and the biogeography of Arisaema was not discussed in detail.

Ohi-Toma et al. (2016Ohi-Toma T, Wu S, Murata H, Murata J. 2016. An updated genus-wide phylogenetic analysis of Arisaema (Araceae) with reference to sections. Botanical Journal of the Linnean Society 182: 100-114.) introduced a new phylogenetic model using maximum parsimony and maximum likelihood methods to analyze 124 species of Arisaema. This work showed that the genus Arisaema includes 15 sections: Anomala, Arisaema, Attenuata, Clavata, Decipientia, Dochafa, Fimbriata, Flagellarisaema, Frenchetiana, Nepenthoidea, Odorata, Pistillata, Sinarisaema, Tenuipistillata and Tortuosa. Fourteen of these sections were recognized by Murata (2011Murata J. 2011. The picture book of plant systematics in color: Arisaema in Japan. Hokuryukan, Tokyo. (in Japanese with taxonomic keys in English), 2013Murata J, Nagamasu H, Ohashi H. 2013. A nomenclatural review on the infrageneric classification of Arisaema (Araceae). Journal Japanese Botany 88: 36-45.), and section Fimbriata was reconfirmed as distinct from other sections in the genus. The updated phylogeny of Arisaema helped to clarify the relationship between phylogeny and geographical distribution for some sections in the genus, but this study did not focus on the evolutionary history of the genus.

Nauheimer et al. (2012Nauheimer L, Metzler D, Renner SS. 2012. Global history of the ancient monocot family Araceae inferred with models accounting for past continental positions and previous ranges based on fossils. New Phytologist 195: 938-950.) studied the global history of the family Araceae. This work included ten species of Arisaema, well calibrated with several fossils of Araceae, and indicated that Araceae began to diversify in the Early Cretaceous, and that all eight subfamilies existed before the K⁄T boundary. Early lineages persisted in Laurasia, and there were several relatively recent entries from Laurasia into Africa, South America, South-East Asia and Australia. Arisaema was estimated to have appeared during the Oligocene, after splitting from the ancestor of the group along with several other genera at the crown of Araceae. However, due to limitations of Arisaema samples, the results of divergence times and biogeography of Arisaema were not mentioned. Thus, the questions on historical biogeography and distribution hypothesis of this genus have not been clarified.

Data are now available to formulate a hypothesis for the relationships and distribution hypothesis of Arisaema species. Here we report phylogenetic and biogeographical analyses based mainly on published molecular data from chloroplast DNA regions from 127 of the 200 species of Arisaema. Our major aims are to reconstruct the biogeographical history of Arisaema reconciled with present day distributions, fossil data, and geological information. We explored the following questions: (1) How old is Arisaema? (2) What is the ancestral area of Arisaema and when were the major biomes colonized? (3) What is the best explanation for Arisaema’s large number of species and wide distribution?

Materials and methods

Sampling, DNA extraction, amplification and sequencing

Four chloroplast DNA regions (trnL-trnF, rpl20-rps12, psbB-psbH and rps2-rpoC2) were studied. The DNA sequences from Arisaema species were last downloaded on 12 September 2020. Duplicate and uncertain sequences were excluded. Two new sequences were obtained for an undescribed species recently found in Vietnam. In total, our data were obtained for 129 species including 127 identified species of Arisaema for all 15 sections of the genus that were recognized by Ohi-Toma et al. (2016Ohi-Toma T, Wu S, Murata H, Murata J. 2016. An updated genus-wide phylogenetic analysis of Arisaema (Araceae) with reference to sections. Botanical Journal of the Linnean Society 182: 100-114.). Thus, the molecular data in this study represents most geographical and taxonomic diversity in the genus with approximately 84 % of Arisaema species. Voucher information and GenBank accession numbers are presented in Tab. S1.

For the sample from Vietnam, we extracted genomic DNA from silica gel dried leaves using the CTAB procedure (Doyle & Doyle 1987Doyle JJ, Doyle JL. 1987. A rapid DNA isolation procedure for small quantities of fresh leaf tissue. Phytochemical bulletin 19: 11-15.). Polymerase chain reactions and sequencing were performed using the primers presented in Ohi-Toma et al. (2016Ohi-Toma T, Wu S, Murata H, Murata J. 2016. An updated genus-wide phylogenetic analysis of Arisaema (Araceae) with reference to sections. Botanical Journal of the Linnean Society 182: 100-114.). PCR products were separated and visualized using an ABI3730 automated sequencer (Applied Biosystems, USA). All sequences were aligned in Geneious v.8.0.5 (Kearse et al. 2012Kearse M, Moir R, Wilson A, et al. 2012. Geneious basic: an integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics 28: 1647-1649.).

Phylogenetic analyses

The phylogenetic analyses of Arisaema were conducted using maximum likelihood (ML) and Bayesian inference (BI) methods. The GTR + I + G substitution model was determined as a best-fitting model for the combined datasets by the Akaike Information Criterion (AIC) as implemented in jModelTest v.2.1.6 (Darriba et al. 2012Darriba D, Taboada GL, Doallo R, Posada D. 2012. JModelTest 2: more models, new heuristics and high-performance computing. Nature Methods 9: 772.) available in the CIPRES Science Gateway (Miller et al. 2010Miller MA, Pfeiffer W, Schwartz T. 2010. Creating the CIPRES Science Gateway for inference of large phylogenetic trees. In: Proceedings of the gateway computing environments workshop (GCE). New Orleans, Institute of Electrical and Electronics Engineers (IEEE). p. 1-8.).

The ML analyses were performed using RAxML v.8.2.8 (Stamatakis 2014Stamatakis A. 2014. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics 30: 1312–1313.) available in the CIPRES Science Gateway (Miller et al. 2010Miller MA, Pfeiffer W, Schwartz T. 2010. Creating the CIPRES Science Gateway for inference of large phylogenetic trees. In: Proceedings of the gateway computing environments workshop (GCE). New Orleans, Institute of Electrical and Electronics Engineers (IEEE). p. 1-8.), with the GTR + I + G substitution model and 1000 bootstrap replicates for statistical support. The Bayesian analysis was performed in MrBayes v.3.1.2 (Ronquist & Huelsenbeck 2003Ronquist F, Huelsenbeck JP. 2003. MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics 19: 1572-1574.) on the CIPRES Science Gateway (Miller et al. 2010Miller MA, Pfeiffer W, Schwartz T. 2010. Creating the CIPRES Science Gateway for inference of large phylogenetic trees. In: Proceedings of the gateway computing environments workshop (GCE). New Orleans, Institute of Electrical and Electronics Engineers (IEEE). p. 1-8.) using the GTR + I + G substitution model. The Markov chain Monte Carlo (MCMC) algorithm was run for 10 million generations with a total of four chains, starting from a random tree, and trees were sampled every 1000 generations. The program Tracer v.1.6 (Rambaut & Drummond 2007Rambaut A, Drummond AJ. 2007. Tracer. Version 1.4. http://beast.bio.ed.ac.uk/Tracer. 11 Oct. 2007.
http://beast.bio.ed.ac.uk/Tracer...
) was used to check that effective sample sizes (ESSs) were obtained for all relevant parameters assessed (> 200). With the first 25 % of sampled generations discarded as burn-in, the 50 % majority-rule consensus tree and Bayesian posterior probabilities (PP) were obtained using the remaining trees.

Divergence time estimation

The divergence times of Arisaema were estimated by using the uncorrelated lognormal Bayesian method in BEAST v.1.8.0 (Drummond et al. 2012Drummond AJ, Suchard MA, Xie D, Rambaut A. 2012. Bayesian phylogenetics with BEAUti and the BEAST 1.7. Molecular Biology and Evolution 29: 1969-1973.). All dating runs relied on a Yule process tree prior. A lognormal distribution was used for the calibration points (Deng et al. 2014Deng JB, Drew BT, Mavrodiev EV, Gitzendanner MA, Soltis PS, Soltis DE. 2014. Phylogeny, divergence times, and historical biogeography of the angiosperm family Saxifragaceae. Molecular Phylogenetics Evolution 83: 86-98.; Liu et al. 2018Liu B, Le CT, Barrett RL, et al. 2018. Historical biogeography of Loranthaceae (Santalales): Diversification agrees with emergence of tropical forests and radiation of songbirds. Molecular Phylogenetics and Evolution 124: 199-212.), and a normal distribution was used for the secondary calibration points (Lu et al. 2013Lu LM, Wang W, Chen ZD, Wen J. 2013. Phylogeny of the non-monophyletic Cayratia Juss. Vitaceae and implications for character evolution and biogeography. Molecular Phylogenetics and Evolution 68: 502-515.; Gonçalves et al. 2020Gonçalves DJP, Shimizu GH, Ortiz EM, Jansen RK, Simpson BB. 2020. Historical biogeography of Vochysiaceae reveals an unexpected perspective of plant evolution in the Neotropics. American Journal of Botany 107: 1-17 ). We conducted two separate MCMC runs of 100 million generations with samples taken every 10,000 generations. Tracer was used to check if ESSs for all relevant parameters were well above 200 and that stationarity had been reached. A maximum credibility tree was then built by TreeAnnotator v.1.8.0 with the initial 25 % of trees discarded as burn-in. The final result was visualized in Figtree v.1.4.0 (Rambaut 2009Rambaut A. 2009. FigTree v.1.4. http://tree.bio.ed.ac.uk/software/figtree/. 21 Dec. 2009.
http://tree.bio.ed.ac.uk/software/figtre...
).

The fossils of Araceae have been reported from a number of Late Cretaceous localities, including several from the Coniacian of Eastern Asia, the Campanian of Central Europe and North America, the Maastrichtian of India and South America (Herrera et al. 2008Herrera FA, Jaramillo CA, Dilcher DL, Wing SL, Gosmez NC. 2008. Fossil Araceae from a Paleocene neotropical rainforest in Colombia. American Journal of Botany 95: 1569-1583.; Nauheimer et al. 2012Nauheimer L, Metzler D, Renner SS. 2012. Global history of the ancient monocot family Araceae inferred with models accounting for past continental positions and previous ranges based on fossils. New Phytologist 195: 938-950.). However, fossils recognized as Arisaema are limited. Knowlton (1925Knowlton FH. 1925. Arisaema hesperia n. sp. U.S. Geological Survey, Reston, Professional Paper 140: 29-30., 1926Knowlton FH. 1926. Flora of the Latah formation of Spokane, Washington, and Coeur d’Alene, Idaho. Virginia, U.S. Geological Survey, Reston, Professional Paper 140A, 1-81.) recognized Arisaema hesperia as a fossil from the Mid Miocene Latah Formation near Spokane, based on the close match with infructescences of living North American A. triphyllum. Renner et al. (2004Renner SS, Zhang LB, Murata J. 2004. A chloroplast phylogeny of Arisaema (Araceae) illustrates Tertiary floristic links between Asia, North America, and east Africa. American Journal of Botany 91: 881-888.) and Nauheimer et al. (2012)Nauheimer L, Metzler D, Renner SS. 2012. Global history of the ancient monocot family Araceae inferred with models accounting for past continental positions and previous ranges based on fossils. New Phytologist 195: 938-950. used this fossil to estimate the divergence times of Arisaema and Araceae, respectively. We thus used Arisaema hesperia to calibrate A. triphyllum at 18 million years ago (Ma). We used a lognormal distribution with the minimum age constraints 18 Ma taken as zero offset, and both the lognormal mean and the standard deviation set to 1 (95 % HPD: 18.52, 32.08 Ma). No further suitable fossils of Arisaema are currently available inside or outside the USA, so the phylogenetic tree was calibrated using a normally distributed prior based on divergence time (Lu et al. 2013Lu LM, Wang W, Chen ZD, Wen J. 2013. Phylogeny of the non-monophyletic Cayratia Juss. Vitaceae and implications for character evolution and biogeography. Molecular Phylogenetics and Evolution 68: 502-515.; Gonçalves et al. 2020Gonçalves DJP, Shimizu GH, Ortiz EM, Jansen RK, Simpson BB. 2020. Historical biogeography of Vochysiaceae reveals an unexpected perspective of plant evolution in the Neotropics. American Journal of Botany 107: 1-17 ). Nauheimer et al. (2012)Nauheimer L, Metzler D, Renner SS. 2012. Global history of the ancient monocot family Araceae inferred with models accounting for past continental positions and previous ranges based on fossils. New Phytologist 195: 938-950. concluded that Arisaema was splitted from the same ancestor of a group that gave-rise to genera Pinellia, Lazarum, Thenophonum, Typhonium, Sauromatum, Eminium, Helicodiceros, Biarum, Arum, and Dracunculus, at ca. 32 Ma (95 % HPD of 24 - 40 Ma) using BEAST. We used this as the divergence age of Arisaema for our secondary calibration of the stem node for this genus at 32 Ma. We assumed a normal distribution with a mean of 32 Ma, and a standard deviation of 4.9, giving a 95 % HPD of 23.94 - 40.06 Ma.

Ancestral area reconstruction

The ancestral area reconstruction of Arisaema was conducted using (1) BioGeoBEARS (Matzke 2013Matzke NJ. 2013. Probabilistic historical biogeography: new models for founderevent speciation, imperfect detection, and fossils allow improved accuracy and model-testing. Frontiers of Biogeography 5: 242-248.) implemented in R (R Development Core Team 2016R Development Core Team. 2016. R: A language and environment for statistical computing. Vienna, Austria, R foundation for statistical computing.), and (2) a Bayesian approach to dispersal-vicariance analysis (Bayes-DIVA; Nylander et al. 2008Nylander JA, Olsson U, Alstrom P, Sanmartín I. 2008. Accounting for phylogenetic uncertainty in biogeography: a Bayesian approach to dispersal vicariance analysis of the thrushes (Aves: Turdus). Systematic Biologists 57: 257-268.) implemented in RASP v.3.2 (Yu et al. 2011Yu Y, Harris AJ, He X. 2011. RASP (Reconstruct Ancestral State in Phylogenies) 2.0 Beta. http://mnh.scu.edu.cn/soft/blog/RASP. 19 Mar. 2012.
http://mnh.scu.edu.cn/soft/blog/RASP...
) using datasets from the BEAST analysis.

For the BioGeoBEARS, two likelihood-based models were implemented: Dispersal-Extinction-Cladogenesis (DEC) (Ree & Smith 2008Ree RH, Smith SA. 2008. Maximum likelihood inference of geographic range evolution by dispersal, local extinction, and cladogenesis. Systematic Biology 57: 4-14.) and the likelihood version of dispersal-vicariance (DIVA) (Ronquist 1997Ronquist F. 1997. Dispersal-vicariance analysis: a new approach to the quantification of historical biogeography. Systematic Biology 46: 195-203. ; herein DIVALIKE). In the two models, an additional ‘‘j” parameter (founder event/jump speciation) was added, which allows descendant lineages to have a different area from the direct ancestor (Matzke 2013Matzke NJ. 2013. Probabilistic historical biogeography: new models for founderevent speciation, imperfect detection, and fossils allow improved accuracy and model-testing. Frontiers of Biogeography 5: 242-248.; Vasconcelos et al. 2017Vasconcelos TNC, Proença CEB, Ahmad B, et al. 2017. Myrteae phylogeny, calibration, biogeography and diversification patterns: increased understanding in the most species rich tribe of Myrtaceae. Molecular Phylogenetics Evolution 109: 113-137. ; Liu et al. 2018Liu B, Le CT, Barrett RL, et al. 2018. Historical biogeography of Loranthaceae (Santalales): Diversification agrees with emergence of tropical forests and radiation of songbirds. Molecular Phylogenetics and Evolution 124: 199-212.). In total, four models resulted in the BioGeoBEARS analysis. In addition, we estimated the number and type of biogeographical events by biogeographical stochastic mapping (BSM) implemented in ‘BioGeoBEARS’ (Matzke 2015Matzke NJ. 2015. Stochastic mapping under biogeographical models. http://phylo.wikidot.com/biogeobears#stochastic_mapping. 1 Jun. 2015.
http://phylo.wikidot.com/biogeobears#sto...
) under the best-fit biogeographical model (see results). According to Dupin et al. (2017Dupin J, Matzke NJ, Särkinen T, et al. 2017. Bayesian estimation of the global biogeographical history of the Solanaceae. Journal Biogeography 44: 887-899.), we divided the biogeographical events under the models including within-area speciation, vicariance, and dispersal events (range expansions and founder events). Frequency of events was estimated from the mean and standard deviation of event counts from 100 BSMs.

The Bayes-DIVA method allows minimal phylogenetic uncertainty by utilizing the prior distribution of trees resulting from BEAST and generating credibility support values for alternative phylogenetic relationships (Nylander et al. 2008Nylander JA, Olsson U, Alstrom P, Sanmartín I. 2008. Accounting for phylogenetic uncertainty in biogeography: a Bayesian approach to dispersal vicariance analysis of the thrushes (Aves: Turdus). Systematic Biologists 57: 257-268.; Yu et al. 2011Yu Y, Harris AJ, He X. 2011. RASP (Reconstruct Ancestral State in Phylogenies) 2.0 Beta. http://mnh.scu.edu.cn/soft/blog/RASP. 19 Mar. 2012.
http://mnh.scu.edu.cn/soft/blog/RASP...
). For the Bayes-DIVA analysis, we performed based on 10,000 trees from the BEAST analysis and computed a condensed tree as the final representative tree with the first 2500 trees discarded as burn-in.

The biogeographic data for species within Arisaema were compiled from the distribution information that describes the extant species of the genus and their relatives in the literatures (Hara 1971Hara H. 1971. A Revision of the Eastern Himalaya species of the genus Arisaema (Araceae). In: Flora of Eastern Himalaya, Second report. Tokyo, University of Tokyo Press. p. 321-354.; Gusman 2006Gusman G. 2006. Arisaema siangense (Araceae), a new species from Arunachal Pradesh, India. Systematics and Geography of Plants 76: 229–233.; Li et al. 2010Li H, Zhu G, Murata J. 2010. Arisaema (Araceae). In: Wu ZY, Raven PH, Hong DY. (eds.) Flora of China. Vol. 23. Beijing and St. Louis, Science Press and Missouri Botanical Garden Press. p. 43-69.), and herbarium specimens. Five biogeographical areas were defined based on the geographic barriers and floristic divisions of Arisaema (Hara 1971Hara H. 1971. A Revision of the Eastern Himalaya species of the genus Arisaema (Araceae). In: Flora of Eastern Himalaya, Second report. Tokyo, University of Tokyo Press. p. 321-354.; McLoughlin 2001McLoughlin S. 2001. The breakup history of Gondwana and its impact on pre-Cenozoic floristic provincialism. Australian Journal Botany 49: 271-300.; Gusman 2006Gusman G. 2006. Arisaema siangense (Araceae), a new species from Arunachal Pradesh, India. Systematics and Geography of Plants 76: 229–233.; Li et al. 2010Li H, Zhu G, Murata J. 2010. Arisaema (Araceae). In: Wu ZY, Raven PH, Hong DY. (eds.) Flora of China. Vol. 23. Beijing and St. Louis, Science Press and Missouri Botanical Garden Press. p. 43-69.) A: Southeast (SE) Asia (including mainland of Southeast Asia, Indochina and Malesia, but excluding New Guinea); B: Africa (including the coastal area of the Arabian Peninsula and East Africa); C: America (including North America, Mexico); D: North China, Japan, the Korean Peninsula region (NJK region); E: Himalayas (including Pakistan, India, Himalaya region, and SW China).

Arisaema is distributed in different areas in three continents (Asia, Africa, and America). There are multiple differences in climatic and geographic barriers among these continents. The morphology and habitat of Arisaema species in Southeast Asia and NJK region or Himalayas are different. Murata (1984Murata J. 1984. An attempt at an infrageneric classification of genus Arisaema (Araceae). Journal of the Faculty of Science, University of Tokyo, Section III. Botany 13: 431-482.) recognized that in Arisaema, tuberous plants with the spilodistichous leaf arrangement are always deciduous and are better adapted to winter dormancy in cool climates. Those characters mostly belong to species distributed in NJK region, America and Himalayas. Murata (1990)Murata J. 1990. Present status of Arisaema systematics. The Botanical Magazine, Tokyo 103: 371-382. also noted that plants with rhizomatous stem and quincuncial leaf arrangement are evergreen plants, better adapted to non-seasonal growth habits under tropical climate. Those characters are found in members mainly distributed in low elevations of Southern China to SE Asia. Therefore, we defined Southeast Asia, NJK region (North China, Japan, and the Korean Peninsula region) and Himalayas as three separate biogeographical areas.

Results

Phylogenetic relationships

The lengths of the four plastid data sets: trnLF, rpl20-rps12, psbB-psbH and rps2-poC2 were 472, 1079, 736, 566 bps respectively. The combined dataset resulted in a matrix with 2853 characters. The topology from ML and BI analyses of the combined dataset were highly congruent, only differing in low resolution nodes. We thus present the BI tree with BS and PP values in Fig. 2.

Figure 2
Majority rule consensus tree of Arisaema based on the combined datasets of four DNA regions (trnL-trnF, rpl20-rps12, psbB-psbH and rps2-rpoC2). ML bootstrap values and posterior probabilities (PP) of the BI analysis are presented above the branches. “-” indicates the support values less than 50%. The circumscription of sections follows the classification of Ohi-Toma et al. (2016Ohi-Toma T, Wu S, Murata H, Murata J. 2016. An updated genus-wide phylogenetic analysis of Arisaema (Araceae) with reference to sections. Botanical Journal of the Linnean Society 182: 100-114.). (A). Arisaema elephas; (B). Arisaema averyanovii; (C). Arisaema yunnanense; (D). Arisaema balansae; (E). Arisaema rhizomatum; (F). Arisaema consanguineum; (G). Arisaema erubescens; (H). Arisaema amurense. Photograph: V. D. Nguyen (B, D, E, F), B. Liu (A, C, G, H).

The phylogenetic results reported here are more advanced compared to other previous analyses (Renner et al. 2004Renner SS, Zhang LB, Murata J. 2004. A chloroplast phylogeny of Arisaema (Araceae) illustrates Tertiary floristic links between Asia, North America, and east Africa. American Journal of Botany 91: 881-888.; Ohi-Toma et al. 2016Ohi-Toma T, Wu S, Murata H, Murata J. 2016. An updated genus-wide phylogenetic analysis of Arisaema (Araceae) with reference to sections. Botanical Journal of the Linnean Society 182: 100-114.) based on the robust dataset. Our phylogenetic analyses supported Arisaema as monophyletic, with 15 sections being recognized within the genus (Fig. 2). Polytomies are distributed in the phylogenetic tree, however, placement and support values for most clades are improved. The basal of the phylogenetic trees includes eight sections Tenuipistillia, Dochafa, Tortuosa, Arisaema, Fimbriata, Attenuata, Odorata, and Anomala. However, the positions of the seven sections are not stable. The three sections Tenuipistillia, Dochafa, and Tortuosa weakly formed a clade. Section Arisaema was well supported as a single clade (BS: 98 %, PP: 1.0) (Fig. 2). While, four sections Fimbriata, Attenuata, Odorata, and Anomala weakly formed a clade, with Fimbriata poorly supported as sister to the three rest sections (Fig. 2). The middle portion of the Arisaema phylogenetic tree includes the remaining seven sections Clavata, Flagellarisaema, Nepenthoidea, Decipientia, Pistillata, Franchetiana and Sinarisaema. The sections Nepenthoidea and Pistillata are not supported as monophyletic. This situation was also seen in Ohi-Toma et al. (2016)Ohi-Toma T, Wu S, Murata H, Murata J. 2016. An updated genus-wide phylogenetic analysis of Arisaema (Araceae) with reference to sections. Botanical Journal of the Linnean Society 182: 100-114..

In the section Anomala, our taxon which was collected in Tuyen Quang province in northern Vietnam (Arisaema sp. VN25) was placed closely relative to A. petelotii, A. hainanense and A. omkoiense, and was strongly supported by the molecular data (Fig. 2). However, our species from Tuyen Quang province is close to A. grapsospadix and A. balansae morphologically. It can be distinguished from A. grapsospadix by having no pseudostem, a spathe without a white area at the base of spathe lamina and from A. balansae by the spathe lacking an auricle at the tube mouth and appendix covered by short neuters. Thus, it should be described as a new species based on the evidence here for future studies.

Divergence times

The divergence time estimations for Arisaema are presented in Fig. 3. The stem age of Arisaema was estimated to be 33.27 Ma (95 % HPD: 29.82, 36.58 Ma). The group including sections Fimbriata, Odorata, Attenuata, and Anomala diverged early at the crown node of Arisaema 27.86 Ma (95 % HPD: 22.86, 32.99 Ma). The African member A. chimperianum diverged 8.83 Ma (95 % HPD: 2.52, 15.51 Ma) at node 6 during late Miocene (Fig. 3), while the American Arisaema are not monophyletic and diverged since the Mid Miocene to Pliocene (Fig. 3). A comparison dating from this study, Nauheimer et al. (2012Nauheimer L, Metzler D, Renner SS. 2012. Global history of the ancient monocot family Araceae inferred with models accounting for past continental positions and previous ranges based on fossils. New Phytologist 195: 938-950.) and Renner et al. (2004Renner SS, Zhang LB, Murata J. 2004. A chloroplast phylogeny of Arisaema (Araceae) illustrates Tertiary floristic links between Asia, North America, and east Africa. American Journal of Botany 91: 881-888.) is presented in Tab. S2.

Figure 3
Maximum clade credibility tree inferred from BEAST based on the combined datasets of four DNA regions. The bars around node ages indicate 95 % highest posterior density intervals. Node constraints are indicated with stars. Nodes of interests were marked as 0-18.

Ancestral area reconstruction

The results of ancestral area reconstruction from BioGeoBEARS and Bayes-DIVA are highly congruent. However, the results from BioGeoBEARS showed higher resolution for more nodes than Bayes-DIVA. Among the four models of BioGeoBEARS, the results indicated higher log likelihood values for three parameters in comparison to two parameters (Tab. 1), thus indicating jump speciation (i.e., dispersal between non-adjacent areas) as an important pattern in range variation of Arisaema. Moreover, the BioGeoBEARS analyses showed DEC+j as the best-fit biogeographical model. Thus, we only present the reconstruction of BioGeoBEARS under the DEC+j model (Fig. 4). The result of ancestral area reconstructions from Bayes-DIVA is shown in Fig. S1. The number of nodes in Fig. 3 and Fig. 4 are consistent, and a summary of the divergence age estimations and ancestral area reconstruction is presented in Tab. 2. Our biogeographical stochastic mapping (BSM) results under the best-fit biogeographical model DEC+j revealed that most biogeographical events comprise within-area speciation (74 %) and dispersals (23 %), with very few (3 %) vicariant events (Tab. S3).

Table 1
Comparison of the fit of different models of biogeographical range evolution and model specific estimates for different parameters (d = dispersal, e = extinction, j = weight of jump dispersal (founder speciation)).

Table 2
Summary of supported clades, divergence age estimation (Ma) by BEAST and ancestral area reconstructions by BioGeoBEARS and Bayes-DIVA. Node numbers refer to Figs. 3, 4. Area abbreviations are as follows: A: Southeast (SE) Asia (including mainland of Southeast Asia, Indochina and Malesia, but excluding New Guinea); B: Africa (including the coastal area of the Arabian Peninsula and East Africa); C: America (including North America, Mexico); D: North China, Japan, the Korean Peninsula region (NJK region); E: Himalayas (including E Pakistan, N India, Himalaya region, and SW China).

Figure 4
Ancestral area reconstruction of Arisaema by BioGeoBEARS (j = 0.000517, LnL = −115.90). Geologic time scale is shown at the bottom. Area abbreviations are as follows: A: Southeast (SE) Asia (including mainland of Southeast Asia, Indochina and Malesia, but excluding New Guinea); B: Africa (including the coastal area of the Arabian Peninsula and East Africa); C: America (including North America, Mexico); D: North China, Japan, the Korean Peninsula region (NJK region); E: Himalayas (including E Pakistan, N India, Himalaya region, and SW China).

Our results suggest that Arisaema most likely originated in Southeast Asia (area A) (node 0, Fig. 4), and subsequently diversified in Southeast Asia during the Oligocene (node 1, Fig. 4). The group including Southeast Asian and Himalayas members of four sections Fimbriata, Anomala, Attenuata and Odorata originated in Southeast Asia at node 1 during the Mid Oligocene. In this clade, A. album, the endemic species of India, diverged in the late Miocene at node 3.

Similarly, all remaining sections likely originated in Southeast Asia and then spread to other areas during the Oligocene to Miocene (Fig. 4). The dispersal events from Southeast Asia to Himalayas, Southeast Asia to NJK region or NJK region to America were presented multiple times within Arisaema (Fig. 4).

Discussion

Origin of Arisaema in Southeast Asia

The biogeographical analyses and divergence time estimations suggest that Arisaema originated in Southeast Asia ca. 33.27 Ma (95 % HPD: 29.82, 36.58 Ma) during the early Oligocene (node 0, Figs. 3, 4; Tab. 2), and soon after the initial origin Arisaema started to diverge at 31.28 Ma (95 % HPD: 27.08, 35.45 Ma; node 1, Figs. 3, 4; Tab. 2). Forests (including tropical forests and temperate forests or mountainous forests at high elevation) were present throughout the Northern Hemisphere during the early Cenozoic, which allowed for intercontinental dispersal of plants (Wei et al. 2015Wei R, Xiang QP, Schneider H, et al. 2015. Eurasian origin, boreotropical migration and transoceanic dispersal in the pantropical fern genus Diplazium (Athyriaceae). Journal Biogeography 42, 1809-1819.; West et al. 2015West CK, Greenwood DR, Basinger JF. 2015. Was the Arctic Eocene ‘rainforest’ monsoonal? Estimates of seasonal precipitation from early Eocene megafloras from Ellesmere Island, Nunavut. Earth Planetary Science Letters 427: 18-30.). In fact, both the tropical forests and temperate forests (including mountainous forests at high elevation) can harbor Arisaema, with many species also occurring in the deciduous forests and subalpine forests in East Asia and East North America. Among the sections of Arisaema, only the sect. Anomala (ca. 25 spp.), sect. Attenuata (ca. 20 spp.) and sect. Fimbriata (ca. 2 spp.) are distributed in tropical forests, others (ca. 150 spp.) in temperate forests, or relatively higher elevation forests of tropical regions (i.e., the mountainous forests in Himalaya, 1000-3000 m) (Murata 1984Murata J. 1984. An attempt at an infrageneric classification of genus Arisaema (Araceae). Journal of the Faculty of Science, University of Tokyo, Section III. Botany 13: 431-482., 1990Murata J. 1990. Present status of Arisaema systematics. The Botanical Magazine, Tokyo 103: 371-382.; Mayo et al. 1997Mayo SJ, Bogner J, Boyce PC. 1997. The genera of Araceae. Kew, Richmond, Royal Botanic Gardens.; Gusman & Gusman 2002Gusman G, Gusman L. 2002. The genus Arisaema. A monograph for botanists and nature lovers. Liechtenstein, Ganter Verlag.; Govaerts et al. 2002Govaerts R, Frodin DG, Bogner J, et al. 2002. Arisaema. World checklist and Bibliography of Araceae (and Acoraceae). Royal Botanic Gardens, Kew, Richmond. p. 181-218.; Murata et al. 2013Murata J, Nagamasu H, Ohashi H. 2013. A nomenclatural review on the infrageneric classification of Arisaema (Araceae). Journal Japanese Botany 88: 36-45.; Ohi-Toma et al. 2016Ohi-Toma T, Wu S, Murata H, Murata J. 2016. An updated genus-wide phylogenetic analysis of Arisaema (Araceae) with reference to sections. Botanical Journal of the Linnean Society 182: 100-114.).

Based on our observation in the field and suggestion from Mayo et al. (1998Mayo SJ, Bogner J, Boyce PC. 1998. Araceae. In: Kubitzki K. (ed.) The families and genera of vascular plants. New York, Springer-Verlag Berlin Heidelberg GmbH. Vol. 4. p. 26-74.), Arisaema and several genera of Araceae prefer wet to swampy conditions, or grow near open water. They prefer growing in woodlands or in shady places under trees in deciduous forests and rainforests. Tropical forests during the Paleocene/Eocene were located at middle paleolatitudes where the climate was more favorable (Morley 2000Morley RJ. 2000. Origin and evolution of tropical rain forests. New York, Wiley.; 2003Morley RJ. 2003. Interplate dispersal paths for megathermal angiosperms. Perspectives in Plant Ecology, Evolution and Systematics 6: 5-20.; Jablonski 2004Jablonski NG. 2004. Primate homeland: forests and the evolution of primates during the Tertiary and Quaternary in Asia. Anthropological Science 113: 117-122.; Schuettpelz & Pryer 2009Schuettpelz E, Pryer KM. 2009. Evidence for a Cenozoic radiation of ferns in an angiosperm-dominated canopy. Proceedings of the National Academy of Sciences USA 106: 11200-11205.; Couvreur et al. 2011Couvreur TLP, Forest F, Baker WJ. 2011. Origin and global diversification patterns of tropical rain forests: inferences from a complete genus-level phylogeny of palms. BMC Biology 9: 44.). By late Eocene to early Oligocene, they retreated to lower latitudes (Contreras et al. 2013Contreras L, Pross J, Bijl PK, et al. 2013. Early to Middle Eocene vegetation dynamics at the Wilkes Land Margin (Antarctica). Review of Palaeobotany and Palynology 197: 119-142.). Modern tropical forests are dominated by angiosperms in taxa such as Fabaceae (Lavin et al. 2005Lavin M., Herendeen PS, Wojciechowski MF. 2005. Evolutionary rates analysis of Leguminosae implicates a rapid diversification of lineages during the Tertiary. Systermatic Biology 54: 575-594. ), Sapindales (Wang et al. 2009Wang H, Moore MJ, Soltis PS, et al. 2009. Rosid radiation and the rapid rise of angiosperm-dominated forests. Proceedings of the National Academy of Sciences USA 106: 3853-3858.), Rubiaceae (Antonelli et al. 2009Antonelli A, Nylander JAA, Persson C, et al. 2009. Tracing the impact of the Andean uplift on Neotropical plant evolution. Proceedings of the National Academy of Sciences USA 106: 9749-9754.) and Fagaceae (Xiang et al. 2014Xiang XG, Wang W, Li RQ, et al. 2014. Large-scale phylogenetic analyses reveal fagalean diversification promoted by the interplay of diaspores and environments in the Paleogene. Perspectives in Plant Ecology, Evolution and Systematics 16: 101-110. ). The woody families above appear to have diverged or radiated before the time of origin of Arisaema and attained a level of ecological prominence comparable to that of today (Crane & Herendeen 2006Crane PR, Herendeen PS. 2006. Cretaceous floras containing angiosperm flowers and fruits from eastern North America. Review of Palaeobotany and Palynology 90: 319-337. ). During the late Eocene and early Oligocene, the climate was warm and moist in most of Southeast Asia (Zachos et al. 2001Zachos JC, Pagani M, Sloan L, Thomas E, Billups K. 2001. Trends, rhythms, and aberrations in global climate 65 Ma to present. Science 292: 686-693.). Thus, the climate condition and forest biomes in this period could have provided habitats for the origin, survival, and radiation differentiation of Arisaema, as it likes to grow in the understory, along forest margins, and in light gaps in forest.

Murata et al. (1993Murata J, Ren C, Murata H, et al. 1993. Arisaema flavum, a honey-producing Araceae. Abstract 1203, 15th International Botanical Congress, Yokohama, Japan.), Wada & Uemura (1994Wada N, Uemura S. 1994. Seed dispersal and predation by small rodents on the herbaceous understory plantSymplocarpus renifolius. American Midland Naturalist 132: 320-327.), Vogel & Martens (2000Vogel S, Martens J. 2000. A survey of the function of the lethal kettle traps of Arisaema (Araceae), with records of pollinating fungus gnats from Nepal. Botanical Journal of the Linnean Society 133: 61-100.), Renner et al. (2004Renner SS, Zhang LB, Murata J. 2004. A chloroplast phylogeny of Arisaema (Araceae) illustrates Tertiary floristic links between Asia, North America, and east Africa. American Journal of Botany 91: 881-888.), Thomason (2020Thomason K. 2020. Rodent seed preference trials of ecologically important plants in Kentucky. Thesis projects. Bowling Green Kentucky, Mahurin Honors College Capstone Experience at Western Kentucky University.) noted that fruits of Arisaema are few-seed berries with bright colors (orange, yellow and red) and represent striking adaptations for birds and small rodents dispersal. Seeds were generally resistant to separate digestive processes in birds and small rodents. Therefore, seeds can be dispersed far to other areas (Kleyheeg et al. 2018Kleyheeg E, Claessens M, Soons MB. 2018. Interactions between seed traits and digestive processes determine the germinability of bird-dispersed seeds. PLOS ONE 13(4): e0195026. ). Jetz et al. (2012Jetz W, Thomas GH, Joy JB, Hartnann K, Mooers AO. 2012. The global diversity of birds in space and time. Nature 491: 444-448.) inferred that birds underwent a major increase in diversification rate since ca. 50 Ma. Therefore, it is likely that the dispersal and diversification of Arisaema was closely related to birds. Additionally, the islands of Southeast Asia appeared during the late Eocene to early Oligocene (Holloway & Hall 1998Holloway DJ, Hall R. 1998. SE Asian geology and biogeography: an introduction. In: Biogeography and Geological Evolution of SE Asia. Leiden, Backhuys Publishers, p. 1-23.; Meimberg et al. 2001Meimberg H, Wistuba A, Dittrich P, Heubl G. 2001. Molecular phylogeny of Nepenthaceae based on cladistic analysis of plastid trnK intron sequence data. Plant Biology 3: 164-175. ; Morley 2003Morley RJ. 2003. Interplate dispersal paths for megathermal angiosperms. Perspectives in Plant Ecology, Evolution and Systematics 6: 5-20.; Clayton et al. 2009Clayton JW, Soltis PS, Soltis DE. 2009. Recent long-distance dispersal overshadows ancient biogeographical patterns in a pantropical angiosperm family (Simaroubaceae, Sapindales). Systematic Biology 58: 395-410.; Hall 2012Hall R. 2012. Late Jurassic-Cenozoic reconstructions of the Indonesian region and the Indian Ocean. Tectonophysics 570-571: 1-41. ; Zahirovic et al. 2014Zahirovic S, Seton M, Müller RD. 2014. The Cretaceous and Cenozoic tectonic evolution of Southeast Asia. Solid Earth 5: 227-273. ). Birds with the ability to travel at least several hundred kilometers may have spread the seeds of Arisaema widely.

In addition, Arisaema species have scents that attract mammals and especially rodents. Some rodent species act as pollinators and seed dispersers for Arisaema (Gusman & Gusman 2006Gusman G, Gusman L. 2006. The genus Arisaema: a monograph for botanists and nature lovers, second revised and enlarged edition. Ruggell, Ganter Verlag.), while Fabre et al. (2012Fabre P, Hautier L, Dimitrov D, Douzery EJP. 2012. A glimpse on the pattern of rodent diversification: a phylogenetic approach. BMC Evolutionary Biology 12: 88.) suggested that many extant families of rodent originated during the Paleogene before the end of the Oligocene (ca. 31 Ma). Oliver et al. (2017Oliver A, Sanisidro O, Baatarjav B, Niiden I, Daxner-Hock G. 2017. Diversification rates in Ctenodactylidae (Rodentia, Mammalia) from Mongolia. Palaeobiodiversity and Palaeoenvironments 97: 51-65.) suggested that Ctenodactylidae (Rodentia) from Mongolia has an initial burst in diversification in the early Oligocene. To settle in the evolving forest, birds and rodents of SE Asia supported the spread and diversification of Arisaema in Southeast Asia (Fig. 4), producing many endemic species in Southeast Asia, and generated high diversity in Asia, especially NJK region as suggested by Govaerts et al. (2002Govaerts R, Frodin DG, Bogner J, et al. 2002. Arisaema. World checklist and Bibliography of Araceae (and Acoraceae). Royal Botanic Gardens, Kew, Richmond. p. 181-218.).

According to results from BSM (Tab. S3), within-area speciation is the most frequent type of biogeographical event in Arisaema. Within-area speciation events appear preponderant in most of the large clades of Arisaema that are endemic to single areas as defined in our study. Additionally, we found that dispersal without “range contractions” was the principal driver of range evolution which occurred about nine times more often than vicariance (Tab. S3). Dupin et al. (2017Dupin J, Matzke NJ, Särkinen T, et al. 2017. Bayesian estimation of the global biogeographical history of the Solanaceae. Journal Biogeography 44: 887-899.) suggested that dispersal events seem to be the most common factor for the distribution of plants, whereas, vicariant events are important factors for range evolution in many animals. Thus dispersal events between both proximal and distant regions appear to have been frequent in the historical biogeography of Arisaema and may have been facilitated by the presence of forests.

On morphology, Hara (1971Hara H. 1971. A Revision of the Eastern Himalaya species of the genus Arisaema (Araceae). In: Flora of Eastern Himalaya, Second report. Tokyo, University of Tokyo Press. p. 321-354.) documented that a primitive type is a narrow cylindrical appendix character. This character occurs mainly in the sections Anomala and Fimbriatum. Both the two sections are distributed mainly in lowland SE Asia and Southern China. Also following Hara (1971)Hara H. 1971. A Revision of the Eastern Himalaya species of the genus Arisaema (Araceae). In: Flora of Eastern Himalaya, Second report. Tokyo, University of Tokyo Press. p. 321-354., this character developed two evolution trends. One trend is appendix becoming more stout and clavate or even capitates forms at the upper portion and truncate at the base with a stipe. Those characters of the appendix are found in the sections Pistillata, Sinnarisaema, Franchetiana, Decipientia, Nepenthoidea, and Dochafa. It means that Arisaema species possibly dispersed from SE Asia following two directions. One direction is going up to NJK region and to America and the other is from SE Asia up to high mountains of South-East Asia, Himalayas and Africa. The other trend is the appendix being elongated and often curved, flagellate tapering to a thin thread. This style of appendix was found in the sections Flagellarisaema, Tenuipistillata, Tortuosa, and Arisaema. All of them are distributed in the Himalayas. The evolutionary trends above can be developed parallel in different regions and different groups.

Migration to Himalayas and Africa from Southeast Asia

In the middle of the tree, the ancestral area reconstruction inferred that Arisaema migrated from Southeast Asia to Himalayas several times (Fig. 4). In addition, we found that Arisaema migrated from Southeast Asia to Himalayas and with no recognized reversed migration from Himalayas even though the two areas are very close. The Indian subcontinent began rifting from Australia-Antarctica ca. 136 Ma (Gibbons et al. 2013Gibbons AD, Whittaker JM, Müller RD. 2013. The breakup of East Gondwana: Assimilating constraints from Cretaceous ocean basins around India into a best-fit tectonic model. Journal Geophysical Research 118: 808-822.), and then collided with Eurasia at 43-50 Ma (McLoughlin 2001McLoughlin S. 2001. The breakup history of Gondwana and its impact on pre-Cenozoic floristic provincialism. Australian Journal Botany 49: 271-300.; Morley 2003Morley RJ. 2003. Interplate dispersal paths for megathermal angiosperms. Perspectives in Plant Ecology, Evolution and Systematics 6: 5-20.; Zahirovic et al. 2014Zahirovic S, Seton M, Müller RD. 2014. The Cretaceous and Cenozoic tectonic evolution of Southeast Asia. Solid Earth 5: 227-273. ). Our results indicated that all the Himalayas samples diverged after 21.28 Ma: A. jacquemontii ca. 21.28 Ma (95 % HPD: 13.88, 28.91 Ma), A. album ca. 8.06 Ma (95 % HPD: 4.22, 12.83 Ma), A. sahyadricum and A. murrayi ca. 1.17 Ma (95 % HPD: 0.06, 3.68 Ma), A. auriculatum and A. shimienense ca. 0.67 Ma (95 % HPD: 0, 3.51 Ma), A. exappendiculatum ca. 3.07 Ma (95 % HPD: 0.78, 6.31 Ma), A. erubescens ca. 1.13 Ma (95 % HPD: 0.03, 3.20 Ma), A. ciliatum ca. 2.76 Ma, much later than the collision between India and Eurasia. Moreover, Nauheimer et al. (2012Nauheimer L, Metzler D, Renner SS. 2012. Global history of the ancient monocot family Araceae inferred with models accounting for past continental positions and previous ranges based on fossils. New Phytologist 195: 938-950.) indicated that all genera of Pistia clade in Aroideae originated in Eurasia (excluding India) except Pistia and Protarum. Thus, the Himalayas members of Arisaema could only be derive from Southeast Asian ancestor by bird and small rodent dispersal (Fig. 4), and the uplift of the Himalayas may have hindered the reverse migration to Asia by the Himalayan members, thus forming more species endemic in that area. For the case of A. album, the endemic species from NE India, our molecular results strongly supported its placement within Attenuata section. During the colonization of Arisaema in Southeast Asia, A. album likely originated by within-area speciation events at 8 Ma after the collision between India and Eurasia when the two areas were continuous and the colonization took place near India.

Additionally, our ancestral area reconstruction analyses indicated that Arisaema migrated from Southeast Asia to Africa only once at node 6 ca. 8.83 Ma (95 % HPD: 2.52, 15.51 Ma) during the late Miocene. The Asia-Africa disjunction could be explained by three major competing hypotheses: (1) migration via the Eocene boreotropical forests (Davis et al. 2002Davis CC, Bell CD, Fritsch PW, Mathews S. 2002. Phylogeny of Acridocarpus, Brachylophon (Malpighiaceae): implications for tertiary tropical floras and Afroasian biogeography. Evolution 56: 2395-2405.; Weeks et al. 2005Weeks A, Daly DC, Simpson BB. 2005. The phylogenetic history and biogeography of the frankincense and myrrh family (Burseraceae) based on nuclear and chloroplast sequence data. Molecular Phylogenetics Evolution 35: 85-101.; Muellner et al. 2006Muellner AN, Savolainen V, Samuel R, Chase MW. 2006. The mahogany family “out-of-Africa”: divergence time estimation, global biogeographic patterns inferred from plastid rbcL DNA sequences, extant, and fossil distribution of diversity. Molecular Phylogenetics and Evolution 40: 236-250.); (2) overland migration via land bridges (Zhou et al. 2012Zhou LL, Su YCF, Thomas DC, Saunders RMK. 2012. ‘Out-of-Africa’ dispersal of tropical floras during the Miocene climatic optimum: evidence from Uvaria (Annonaceae). Journal Biogeography 39: 322-335.); and (3) recent transoceanic long-distance dispersal (Renner 2004Renner SS. 2004. Plant dispersal across the tropical Atlantic by wind and sea currents. International Journal of Plant Sciences 165: S23-S33.; Yuan et al. 2005Yuan YM, Wohlhauser S, Moller M, Klackenberg J, Callmander MW, Kupfer P. 2005. Phylogeny and biogeography of Exacum (Gentianaceae): a disjunctive distribution in the Indian ocean basin resulted from long distance dispersal and extensive radiation. Systematics Biology 54: 21-34.; Clayton et al. 2009Clayton JW, Soltis PS, Soltis DE. 2009. Recent long-distance dispersal overshadows ancient biogeographical patterns in a pantropical angiosperm family (Simaroubaceae, Sapindales). Systematic Biology 58: 395-410.). The first hypothesis relied on the existence of extensive boreotropical forests from Europe across South Asia associated with the Eocene thermal maximum ca. 50-52 Ma (Zachos et al. 2001Zachos JC, Pagani M, Sloan L, Thomas E, Billups K. 2001. Trends, rhythms, and aberrations in global climate 65 Ma to present. Science 292: 686-693.; 2008Zachos JC, Gerald RD, Richard EZ. 2008. An early Cenozoic perspective on greenhouse: warming and carbon-cycle dynamics. Nature 451: 279-283.). During the optimum climatic period, most regions of northern Africa, including areas now occupied by the Sahara, were dominated by rainforests (Raven & Axelrod 1974Raven PH, Axelrod DI. 1974. Angiosperm biogeography and past continental movements. Annals Missouri Botany Garden 61: 539-673.; Axelrod & Raven 1978Axelrod DI, Raven PH. 1978. Late Cretaceous and Tertiary vegetation history of Africa. In: Werger MJA. (ed.) Biogeography and ecology of southern Africa. The Hague, Dr. W. Junk b.v. Publishers. p. 77-130.). These tropical forests could have extended much further north than extant forests, enabling floristic exchanges between Africa and Asia in both directions (Wolfe 1975Wolfe JA. 1975. Some aspects of plant geography of the Northern Hemisphere during the Late Cretaceous and Tertiary. Annals Missouri Botany Garden 62: 264-279.; Morley 2000Morley RJ. 2000. Origin and evolution of tropical rain forests. New York, Wiley.). This hypothesis has been used to explain the disjunction between Asia and Africa in some angiosperm families (Davis et al. 2002Davis CC, Bell CD, Fritsch PW, Mathews S. 2002. Phylogeny of Acridocarpus, Brachylophon (Malpighiaceae): implications for tertiary tropical floras and Afroasian biogeography. Evolution 56: 2395-2405., 2005Davis CC, Webb CO, Wurdack KJ, Jaramillo CA, Donoghue MJ. 2005. Explosive radiation of Malpighiales supports a Mid-Cretaceous origin of modern tropical rain forests. American Naturalist 165: E36-E65.; Weeks et al. 2005Weeks A, Daly DC, Simpson BB. 2005. The phylogenetic history and biogeography of the frankincense and myrrh family (Burseraceae) based on nuclear and chloroplast sequence data. Molecular Phylogenetics Evolution 35: 85-101.; Lu et al. 2013Lu LM, Wang W, Chen ZD, Wen J. 2013. Phylogeny of the non-monophyletic Cayratia Juss. Vitaceae and implications for character evolution and biogeography. Molecular Phylogenetics and Evolution 68: 502-515.; Liu et al. 2018Liu B, Le CT, Barrett RL, et al. 2018. Historical biogeography of Loranthaceae (Santalales): Diversification agrees with emergence of tropical forests and radiation of songbirds. Molecular Phylogenetics and Evolution 124: 199-212.; Zhou et al. 2019Zhou Z, Hu JJ, Wen J, Sun H. 2019. Morphometric, phylogenetic and biogeographic analyses of Pyrularia (Santalales), a parasitic disjunct lineage between eastern Asia and eastern North America. Taxon 64: 47-71.). However, the time of this expansion of boreotropical forests did not coincide with the Asian - African divergence time of Arisaema during the late Miocene. The second hypothesis “overland migration via land bridges” indicated that in the early Miocene ca. 20 Ma, the collision between the Africa continent with the Eurasia continent and the tentative closure of the Tethys sea enabled land connections and land bridge formation between Africa and Southwest Asia. These geological events, coinciding with a warm phase in the middle Miocene climatic optimum, facilitated extensive tropical biotic exchange between Asia and Africa across the Arabian Peninsula. However, the divergence time of Arisaema is much later than the time of the hypothesis. The last hypothesis “recent transoceanic long-distance dispersal” has been used to explain the disjunction between Asia and Africa in several genera of angiosperms: Adansonia (Malvaceae) (Baum et al. 1998Baum DA, Small RL, Wendel JF. 1998. Biogeography and floral evolution of baobabs (Adansonia, Bombacaceae) as inferred from multiple data sets. Systematic Biology 47: 181-207.); Osbeckia (Melastomataceae) (Renner & Meyer 2001Renner SS, Meyer K. 2001. Melastomataceae come full circle: biogeographic reconstruction and molecular clock dating. Evolution 55: 1315-1324.; Renner 2004Renner SS. 2004. Plant dispersal across the tropical Atlantic by wind and sea currents. International Journal of Plant Sciences 165: S23-S33.); Exacum (Gentianaceae) (Yuan et al. 2005Yuan YM, Wohlhauser S, Moller M, Klackenberg J, Callmander MW, Kupfer P. 2005. Phylogeny and biogeography of Exacum (Gentianaceae): a disjunctive distribution in the Indian ocean basin resulted from long distance dispersal and extensive radiation. Systematics Biology 54: 21-34.); Eurycoma, Brucea and Soulamea (Simaroubaceae) (Clayton et al. 2009Clayton JW, Soltis PS, Soltis DE. 2009. Recent long-distance dispersal overshadows ancient biogeographical patterns in a pantropical angiosperm family (Simaroubaceae, Sapindales). Systematic Biology 58: 395-410.); Loranthaceae (Liu et al. 2018Liu B, Le CT, Barrett RL, et al. 2018. Historical biogeography of Loranthaceae (Santalales): Diversification agrees with emergence of tropical forests and radiation of songbirds. Molecular Phylogenetics and Evolution 124: 199-212.). This hypothesis may be more feasible for Arisaema, with a long-distance dispersal from Southeast Asia to Africa at ca. 8.83 Ma (95 % HPD: 2.52, 15.51 Ma) in the late Miocene (node 6 in Tab. 2 and Fig. 4). The taxa involved in these migrations possibly evolved into the Asian endemic species, African endemic species and species shared by both areas.

In addition, our results of the divergence between African A. schimperianum and Himalayan A. costatum at 8.83 Ma (95 % HPD: 2.52, 15.51 Ma) (node 6, Figs. 3, 4) is younger than Renner et al. (2004Renner SS, Zhang LB, Murata J. 2004. A chloroplast phylogeny of Arisaema (Araceae) illustrates Tertiary floristic links between Asia, North America, and east Africa. American Journal of Botany 91: 881-888.) (19 to 31 Ma). Renner et al. (2004)Renner SS, Zhang LB, Murata J. 2004. A chloroplast phylogeny of Arisaema (Araceae) illustrates Tertiary floristic links between Asia, North America, and east Africa. American Journal of Botany 91: 881-888. noted that their age estimates between Africa and Himalayas should be improved in future study. Thus, our result is more reliable than Renner et al. (2004)Renner SS, Zhang LB, Murata J. 2004. A chloroplast phylogeny of Arisaema (Araceae) illustrates Tertiary floristic links between Asia, North America, and east Africa. American Journal of Botany 91: 881-888. with the advances of fossil constraint, analysis methods and intensive sampling.

Dispersal from Southeast Asia to North China, Japan and the Korean Peninsula and vice versa

Our ancestral area reconstruction inferred that Arisaema dispersal from Southeast Asia to NJK region and vice versa multiple times since the early Miocene at node 9 ca. 21.27 Ma (Fig. 4). The global climate has been warm since the late Oligocene and early Miocene before it became optimum in the Mid Miocene; after that period, the climate started cooling (Zachos et al. 2001Zachos JC, Pagani M, Sloan L, Thomas E, Billups K. 2001. Trends, rhythms, and aberrations in global climate 65 Ma to present. Science 292: 686-693.). The forests have been present since the early Cenozoic in the Northern Hemisphere. Moreover, the two areas are very close, and always contiguous. Thus we hypothesize that Southeast Asian Arisaema was dispersed to the NJK region during a warm climate and via moist forest areas.

Based on our observation in the field, Arisaema species in subtropical and tropical areas are evergreen leaves, whereas members in temperature or mountainous forests at high elevation have seasonal leaves. Our results recognized that in the NJK region, Arisaema diverged during periods of global cooling of Miocene and Pliocene (Fig. 4). In fact, based on our observations in the field, Arisaema in the NJK region or in mountainous forests at high elevation always has seasonal leaves. This could be a mechanism of adaptation to living conditions when dispersal from tropical region to temperate region or mountainous forests at high elevation regions. In addition, the mechanism of adaptation to living conditions can be made the diversification of Arisaema in the NJK region (Figs. 4, S1).

Dispersal to America from North China, Japan and the Korean Peninsula during Miocene

The biogeographical analyses show that the colonization to America of Arisaema has occurred multiple times (Fig. 4). The disjunction of Arisaema between the NJK region and America can be explained by the hypothesis of the Bering land bridge. The Bering land bridge was mostly exposed during the Neogene for trans-Pacific floristic changes until about 3.5 Ma (Wen 1999Wen J. 1999. Evolution of eastern Asian and eastern North American disjunct distributions in flowering plants. Annual Review Ecology Evolution Systematics 30: 421-455.; Ickert-Bond et al. 2009Ickert-Bond SM, Murray DF, DeChaine E. 2009. Contrasting patterns of plant distribution in Beringia. Alaska Park Science 8: 26-32.). Our ancestral area reconstruction indicated that Arisaema dispersal to America from NJK region ca. 18 Ma and 5 Ma (Figs. 3, 4). The sharp global temperature was decreased in the Paleogene (Zachos et al. 2001Zachos JC, Pagani M, Sloan L, Thomas E, Billups K. 2001. Trends, rhythms, and aberrations in global climate 65 Ma to present. Science 292: 686-693.) and subsequent events, such as continental-wide glaciations, drying of the mid-continental region of North America and orogenies in western North America and the western part of eastern Asia. However, global climate optimum in the mid Miocene (ca. 17-15 Ma) could support the dispersal of Arisaema. Thus we hypothesized that Arisaema was dispersed to America from the NJK region by bird via the Bering land bridge during the Miocene. These results highlight that active migrations during favorable periods may have resulted in the widespread distributions of many lineages. Our results also indicate the importance of the Miocene in the Northern Hemisphere intercontinental disjunction of temperate and tropical lineages.

Conclusions

This study represents a comprehensive phylogenetic and biogeographic reconstruction of the genus Arisaema, providing an understanding of the evolutionary history of the genus through time and space. Arisaema originated in Southeast Asia at ca. 31 Ma during the Oligocene with the presence and development of the tropical rainforest and suitable climatic conditions. Important events in the historical biogeography of Arisaema occurred from Oligocene to Miocene, with Southeast Asia as the center of origin and diversification of the genus in North China, Japan and the Korean Peninsula. Finally, this study demonstrates the usefulness of using large secondary molecular datasets and biogeography to determine the evolutionary origin and diversification of plant taxa that occur in Asia.

Acknowledgments

This research is funded by Graduate University of Science and Technology under grant number GUST.STS.ĐT2020-ST01. We also thank Dr. Lars Nauheimer of Australian Tropical Herbarium, Australia, and Dr. Peter C. Boyce from Ludwig-Maximilians-Universität München, Germany for their help in initiating this study and critical reading of the manuscript.

References

  • Antonelli A, Nylander JAA, Persson C, et al 2009. Tracing the impact of the Andean uplift on Neotropical plant evolution. Proceedings of the National Academy of Sciences USA 106: 9749-9754.
  • Axelrod DI, Raven PH. 1978. Late Cretaceous and Tertiary vegetation history of Africa. In: Werger MJA. (ed.) Biogeography and ecology of southern Africa. The Hague, Dr. W. Junk b.v. Publishers. p. 77-130.
  • Baum DA, Small RL, Wendel JF. 1998. Biogeography and floral evolution of baobabs (Adansonia, Bombacaceae) as inferred from multiple data sets. Systematic Biology 47: 181-207.
  • Clayton JW, Soltis PS, Soltis DE. 2009. Recent long-distance dispersal overshadows ancient biogeographical patterns in a pantropical angiosperm family (Simaroubaceae, Sapindales). Systematic Biology 58: 395-410.
  • Contreras L, Pross J, Bijl PK, et al 2013. Early to Middle Eocene vegetation dynamics at the Wilkes Land Margin (Antarctica). Review of Palaeobotany and Palynology 197: 119-142.
  • Couvreur TLP, Forest F, Baker WJ. 2011. Origin and global diversification patterns of tropical rain forests: inferences from a complete genus-level phylogeny of palms. BMC Biology 9: 44.
  • Crane PR, Herendeen PS. 2006. Cretaceous floras containing angiosperm flowers and fruits from eastern North America. Review of Palaeobotany and Palynology 90: 319-337.
  • Darriba D, Taboada GL, Doallo R, Posada D. 2012. JModelTest 2: more models, new heuristics and high-performance computing. Nature Methods 9: 772.
  • Davis CC, Bell CD, Fritsch PW, Mathews S. 2002. Phylogeny of Acridocarpus, Brachylophon (Malpighiaceae): implications for tertiary tropical floras and Afroasian biogeography. Evolution 56: 2395-2405.
  • Davis CC, Webb CO, Wurdack KJ, Jaramillo CA, Donoghue MJ. 2005. Explosive radiation of Malpighiales supports a Mid-Cretaceous origin of modern tropical rain forests. American Naturalist 165: E36-E65.
  • Deng JB, Drew BT, Mavrodiev EV, Gitzendanner MA, Soltis PS, Soltis DE. 2014. Phylogeny, divergence times, and historical biogeography of the angiosperm family Saxifragaceae. Molecular Phylogenetics Evolution 83: 86-98.
  • Doyle JJ, Doyle JL. 1987. A rapid DNA isolation procedure for small quantities of fresh leaf tissue. Phytochemical bulletin 19: 11-15.
  • Drummond AJ, Suchard MA, Xie D, Rambaut A. 2012. Bayesian phylogenetics with BEAUti and the BEAST 1.7. Molecular Biology and Evolution 29: 1969-1973.
  • Dupin J, Matzke NJ, Särkinen T, et al 2017. Bayesian estimation of the global biogeographical history of the Solanaceae. Journal Biogeography 44: 887-899.
  • Engler A. 1920. Araceae-Aroideae, Araceae-Pistioidea. In: Engler A. (ed.) Das Pflanzenreich. Vol. 73. Leipzig, Verlag von Wilhelm Engelmann. p. 149-220.
  • Fabre P, Hautier L, Dimitrov D, Douzery EJP. 2012. A glimpse on the pattern of rodent diversification: a phylogenetic approach. BMC Evolutionary Biology 12: 88.
  • Gibbons AD, Whittaker JM, Müller RD. 2013. The breakup of East Gondwana: Assimilating constraints from Cretaceous ocean basins around India into a best-fit tectonic model. Journal Geophysical Research 118: 808-822.
  • Gonçalves DJP, Shimizu GH, Ortiz EM, Jansen RK, Simpson BB. 2020. Historical biogeography of Vochysiaceae reveals an unexpected perspective of plant evolution in the Neotropics. American Journal of Botany 107: 1-17
  • Govaerts R, Frodin DG, Bogner J, et al 2002. Arisaema World checklist and Bibliography of Araceae (and Acoraceae). Royal Botanic Gardens, Kew, Richmond. p. 181-218.
  • Gusman G. 2006. Arisaema siangense (Araceae), a new species from Arunachal Pradesh, India. Systematics and Geography of Plants 76: 229–233.
  • Gusman G, Gusman L. 2002. The genus Arisaema A monograph for botanists and nature lovers. Liechtenstein, Ganter Verlag.
  • Gusman G, Gusman L. 2006. The genus Arisaema: a monograph for botanists and nature lovers, second revised and enlarged edition. Ruggell, Ganter Verlag.
  • Hall R. 2012. Late Jurassic-Cenozoic reconstructions of the Indonesian region and the Indian Ocean. Tectonophysics 570-571: 1-41.
  • Hara H. 1971. A Revision of the Eastern Himalaya species of the genus Arisaema (Araceae). In: Flora of Eastern Himalaya, Second report. Tokyo, University of Tokyo Press. p. 321-354.
  • Herrera FA, Jaramillo CA, Dilcher DL, Wing SL, Gosmez NC. 2008. Fossil Araceae from a Paleocene neotropical rainforest in Colombia. American Journal of Botany 95: 1569-1583.
  • Holloway DJ, Hall R. 1998. SE Asian geology and biogeography: an introduction. In: Biogeography and Geological Evolution of SE Asia. Leiden, Backhuys Publishers, p. 1-23.
  • Ickert-Bond SM, Murray DF, DeChaine E. 2009. Contrasting patterns of plant distribution in Beringia. Alaska Park Science 8: 26-32.
  • Jablonski NG. 2004. Primate homeland: forests and the evolution of primates during the Tertiary and Quaternary in Asia. Anthropological Science 113: 117-122.
  • Jetz W, Thomas GH, Joy JB, Hartnann K, Mooers AO. 2012. The global diversity of birds in space and time. Nature 491: 444-448.
  • Kearse M, Moir R, Wilson A, et al 2012. Geneious basic: an integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics 28: 1647-1649.
  • Kleyheeg E, Claessens M, Soons MB. 2018. Interactions between seed traits and digestive processes determine the germinability of bird-dispersed seeds. PLOS ONE 13(4): e0195026.
  • Knowlton FH. 1925. Arisaema hesperia n. sp U.S. Geological Survey, Reston, Professional Paper 140: 29-30.
  • Knowlton FH. 1926. Flora of the Latah formation of Spokane, Washington, and Coeur d’Alene, Idaho. Virginia, U.S. Geological Survey, Reston, Professional Paper 140A, 1-81.
  • Lavin M., Herendeen PS, Wojciechowski MF. 2005. Evolutionary rates analysis of Leguminosae implicates a rapid diversification of lineages during the Tertiary. Systermatic Biology 54: 575-594.
  • Li H, Zhu G, Murata J. 2010. Arisaema (Araceae). In: Wu ZY, Raven PH, Hong DY. (eds.) Flora of China. Vol. 23. Beijing and St. Louis, Science Press and Missouri Botanical Garden Press. p. 43-69.
  • Liu B, Le CT, Barrett RL, et al 2018. Historical biogeography of Loranthaceae (Santalales): Diversification agrees with emergence of tropical forests and radiation of songbirds. Molecular Phylogenetics and Evolution 124: 199-212.
  • Lu LM, Wang W, Chen ZD, Wen J. 2013. Phylogeny of the non-monophyletic Cayratia Juss. Vitaceae and implications for character evolution and biogeography. Molecular Phylogenetics and Evolution 68: 502-515.
  • Matzke NJ. 2013. Probabilistic historical biogeography: new models for founderevent speciation, imperfect detection, and fossils allow improved accuracy and model-testing. Frontiers of Biogeography 5: 242-248.
  • Matzke NJ. 2015. Stochastic mapping under biogeographical models. http://phylo.wikidot.com/biogeobears#stochastic_mapping 1 Jun. 2015.
    » http://phylo.wikidot.com/biogeobears#stochastic_mapping
  • Mayo SJ, Bogner J, Boyce PC. 1997. The genera of Araceae. Kew, Richmond, Royal Botanic Gardens.
  • Mayo SJ, Bogner J, Boyce PC. 1998. Araceae. In: Kubitzki K. (ed.) The families and genera of vascular plants. New York, Springer-Verlag Berlin Heidelberg GmbH. Vol. 4. p. 26-74.
  • McLoughlin S. 2001. The breakup history of Gondwana and its impact on pre-Cenozoic floristic provincialism. Australian Journal Botany 49: 271-300.
  • Meimberg H, Wistuba A, Dittrich P, Heubl G. 2001. Molecular phylogeny of Nepenthaceae based on cladistic analysis of plastid trnK intron sequence data. Plant Biology 3: 164-175.
  • Miller MA, Pfeiffer W, Schwartz T. 2010. Creating the CIPRES Science Gateway for inference of large phylogenetic trees. In: Proceedings of the gateway computing environments workshop (GCE). New Orleans, Institute of Electrical and Electronics Engineers (IEEE). p. 1-8.
  • Morley RJ. 2000. Origin and evolution of tropical rain forests. New York, Wiley.
  • Morley RJ. 2003. Interplate dispersal paths for megathermal angiosperms. Perspectives in Plant Ecology, Evolution and Systematics 6: 5-20.
  • Muellner AN, Savolainen V, Samuel R, Chase MW. 2006. The mahogany family “out-of-Africa”: divergence time estimation, global biogeographic patterns inferred from plastid rbcL DNA sequences, extant, and fossil distribution of diversity. Molecular Phylogenetics and Evolution 40: 236-250.
  • Murata J. 1984. An attempt at an infrageneric classification of genus Arisaema (Araceae). Journal of the Faculty of Science, University of Tokyo, Section III. Botany 13: 431-482.
  • Murata J. 1990. Present status of Arisaema systematics. The Botanical Magazine, Tokyo 103: 371-382.
  • Murata J. 1991. The systematic position of Arisaema nepenthoides and A. wattii (Araceae). Kew Bulletin 46: 119-128.
  • Murata J. 2011. The picture book of plant systematics in color: Arisaema in Japan. Hokuryukan, Tokyo. (in Japanese with taxonomic keys in English)
  • Murata J, Ren C, Murata H, et al 1993. Arisaema flavum, a honey-producing Araceae. Abstract 1203, 15th International Botanical Congress, Yokohama, Japan.
  • Murata J, Nagamasu H, Ohashi H. 2013. A nomenclatural review on the infrageneric classification of Arisaema (Araceae). Journal Japanese Botany 88: 36-45.
  • Nakai T. 1950. Classes, ordines, familiae, subfamiliae, tribus, genera nova quae attient ad plantas koreanas (Supplementum). Journal Japanese Botany 25: 5-7.
  • Nauheimer L, Metzler D, Renner SS. 2012. Global history of the ancient monocot family Araceae inferred with models accounting for past continental positions and previous ranges based on fossils. New Phytologist 195: 938-950.
  • Nylander JA, Olsson U, Alstrom P, Sanmartín I. 2008. Accounting for phylogenetic uncertainty in biogeography: a Bayesian approach to dispersal vicariance analysis of the thrushes (Aves: Turdus). Systematic Biologists 57: 257-268.
  • Ohi-Toma T, Wu S, Murata H, Murata J. 2016. An updated genus-wide phylogenetic analysis of Arisaema (Araceae) with reference to sections. Botanical Journal of the Linnean Society 182: 100-114.
  • Oliver A, Sanisidro O, Baatarjav B, Niiden I, Daxner-Hock G. 2017. Diversification rates in Ctenodactylidae (Rodentia, Mammalia) from Mongolia. Palaeobiodiversity and Palaeoenvironments 97: 51-65.
  • R Development Core Team. 2016. R: A language and environment for statistical computing. Vienna, Austria, R foundation for statistical computing.
  • Rambaut A. 2009. FigTree v.1.4. http://tree.bio.ed.ac.uk/software/figtree/ 21 Dec. 2009.
    » http://tree.bio.ed.ac.uk/software/figtree/
  • Rambaut A, Drummond AJ. 2007. Tracer. Version 1.4. http://beast.bio.ed.ac.uk/Tracer 11 Oct. 2007.
    » http://beast.bio.ed.ac.uk/Tracer
  • Raven PH, Axelrod DI. 1974. Angiosperm biogeography and past continental movements. Annals Missouri Botany Garden 61: 539-673.
  • Ree RH, Smith SA. 2008. Maximum likelihood inference of geographic range evolution by dispersal, local extinction, and cladogenesis. Systematic Biology 57: 4-14.
  • Renner SS. 2004. Plant dispersal across the tropical Atlantic by wind and sea currents. International Journal of Plant Sciences 165: S23-S33.
  • Renner SS, Meyer K. 2001. Melastomataceae come full circle: biogeographic reconstruction and molecular clock dating. Evolution 55: 1315-1324.
  • Renner SS, Zhang LB. 2004. Biogeography of the Pistia clade (Araceae): based on chloroplast and mitochondrial DNA sequences and Bayesian divergence time inference. Systematic Biology 53: 422-432.
  • Renner SS, Zhang LB, Murata J. 2004. A chloroplast phylogeny of Arisaema (Araceae) illustrates Tertiary floristic links between Asia, North America, and east Africa. American Journal of Botany 91: 881-888.
  • Ronquist F. 1997. Dispersal-vicariance analysis: a new approach to the quantification of historical biogeography. Systematic Biology 46: 195-203.
  • Ronquist F, Huelsenbeck JP. 2003. MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics 19: 1572-1574.
  • Schott HW. 1860. Prodromus systematis aroidearum. Typis congregationis. Vienna, Mechitharisticae.
  • Schuettpelz E, Pryer KM. 2009. Evidence for a Cenozoic radiation of ferns in an angiosperm-dominated canopy. Proceedings of the National Academy of Sciences USA 106: 11200-11205.
  • Stamatakis A. 2014. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics 30: 1312–1313.
  • Thomason K. 2020. Rodent seed preference trials of ecologically important plants in Kentucky. Thesis projects. Bowling Green Kentucky, Mahurin Honors College Capstone Experience at Western Kentucky University.
  • Vasconcelos TNC, Proença CEB, Ahmad B, et al 2017. Myrteae phylogeny, calibration, biogeography and diversification patterns: increased understanding in the most species rich tribe of Myrtaceae. Molecular Phylogenetics Evolution 109: 113-137.
  • Vogel S, Martens J. 2000. A survey of the function of the lethal kettle traps of Arisaema (Araceae), with records of pollinating fungus gnats from Nepal. Botanical Journal of the Linnean Society 133: 61-100.
  • Wada N, Uemura S. 1994. Seed dispersal and predation by small rodents on the herbaceous understory plantSymplocarpus renifolius American Midland Naturalist 132: 320-327.
  • Wang H, Moore MJ, Soltis PS, et al 2009. Rosid radiation and the rapid rise of angiosperm-dominated forests. Proceedings of the National Academy of Sciences USA 106: 3853-3858.
  • Wei R, Xiang QP, Schneider H, et al 2015. Eurasian origin, boreotropical migration and transoceanic dispersal in the pantropical fern genus Diplazium (Athyriaceae). Journal Biogeography 42, 1809-1819.
  • Weeks A, Daly DC, Simpson BB. 2005. The phylogenetic history and biogeography of the frankincense and myrrh family (Burseraceae) based on nuclear and chloroplast sequence data. Molecular Phylogenetics Evolution 35: 85-101.
  • Wen J. 1999. Evolution of eastern Asian and eastern North American disjunct distributions in flowering plants. Annual Review Ecology Evolution Systematics 30: 421-455.
  • West CK, Greenwood DR, Basinger JF. 2015. Was the Arctic Eocene ‘rainforest’ monsoonal? Estimates of seasonal precipitation from early Eocene megafloras from Ellesmere Island, Nunavut. Earth Planetary Science Letters 427: 18-30.
  • Wolfe JA. 1975. Some aspects of plant geography of the Northern Hemisphere during the Late Cretaceous and Tertiary. Annals Missouri Botany Garden 62: 264-279.
  • Xiang XG, Wang W, Li RQ, et al 2014. Large-scale phylogenetic analyses reveal fagalean diversification promoted by the interplay of diaspores and environments in the Paleogene. Perspectives in Plant Ecology, Evolution and Systematics 16: 101-110.
  • Yu Y, Harris AJ, He X. 2011. RASP (Reconstruct Ancestral State in Phylogenies) 2.0 Beta. http://mnh.scu.edu.cn/soft/blog/RASP 19 Mar. 2012.
    » http://mnh.scu.edu.cn/soft/blog/RASP
  • Yuan YM, Wohlhauser S, Moller M, Klackenberg J, Callmander MW, Kupfer P. 2005. Phylogeny and biogeography of Exacum (Gentianaceae): a disjunctive distribution in the Indian ocean basin resulted from long distance dispersal and extensive radiation. Systematics Biology 54: 21-34.
  • Zachos JC, Pagani M, Sloan L, Thomas E, Billups K. 2001. Trends, rhythms, and aberrations in global climate 65 Ma to present. Science 292: 686-693.
  • Zachos JC, Gerald RD, Richard EZ. 2008. An early Cenozoic perspective on greenhouse: warming and carbon-cycle dynamics. Nature 451: 279-283.
  • Zahirovic S, Seton M, Müller RD. 2014. The Cretaceous and Cenozoic tectonic evolution of Southeast Asia. Solid Earth 5: 227-273.
  • Zhou LL, Su YCF, Thomas DC, Saunders RMK. 2012. ‘Out-of-Africa’ dispersal of tropical floras during the Miocene climatic optimum: evidence from Uvaria (Annonaceae). Journal Biogeography 39: 322-335.
  • Zhou Z, Hu JJ, Wen J, Sun H. 2019. Morphometric, phylogenetic and biogeographic analyses of Pyrularia (Santalales), a parasitic disjunct lineage between eastern Asia and eastern North America. Taxon 64: 47-71.

Publication Dates

  • Publication in this collection
    02 May 2022
  • Date of issue
    2022

History

  • Received
    05 Aug 2021
  • Accepted
    21 Jan 2022
Sociedade Botânica do Brasil SCLN 307 - Bloco B - Sala 218 - Ed. Constrol Center Asa Norte CEP: 70746-520 Brasília/DF. - Alta Floresta - MT - Brazil
E-mail: acta@botanica.org.br