Abstract
The determinants of biodiversity patterns can be understood using macroevolutionary analyses. The integration of fossils into phylogenies offers a deeper understanding of processes underlying biodiversity patterns in deep time. Cycadales are considered a relict of a once more diverse and globally distributed group but are restricted to low latitudes today. We still know little about their origin and geographic range evolution. Combining molecular data for extant species and leaf morphological data for extant and fossil species, we study the origin of cycad global biodiversity patterns through Bayesian totalevidence dating analyses. We assess the ancestral geographic origin and trace the historical biogeography of cycads with a time-stratified process-based model. Cycads originated in the Carboniferous on the Laurasian landmass and expanded in Gondwana in the Jurassic. Through now-vanished continental connections, Antarctica and Greenland were crucial biogeographic crossroads for cycad biogeography. Vicariance is an essential speciation mode in the deep and recent past. Their latitudinal span increased in the Jurassic and restrained toward subtropical latitudes in the Neogene in line with biogeographic inferences of high-latitude extirpations. We show the benefits of integrating fossils into phylogenies to estimate ancestral areas of origin and to study evolutionary processes explaining the global distribution of present-day relict groups.