Phytogeographic history of the Tea family inferred through high-resolution phylogeny and fossils

Bayesian tip-dated phylogeny and biogeography of Cissampelideae (Menispermaceae): Mitigating the effects of homoplastic morphological characters

The integration of morphological and molecular data is essential to understand the affinities of fossil taxa and spatio-temporal evolutionary processes of organisms. However, homoplastic morphological characters can mislead the placement of fossil taxa and impact downstream analyses. Here, we provide an example of how to mitigate effectively the effect of morphological homoplasy on the placement of fossil taxa and biogeographic inferences of Cissampelideae. We assembled three data types, morphological data only, morphological data with a molecular scaffold and combined morphological and molecular data. By removing high-level homoplastic morphological data or reweighting the morphological characters, we conducted 15 parsimony, 12 undated Bayesian and four dated Bayesian analyses. Our results show that the 14 selected Cissampelideae fossil taxa are placed poorly when based only on morphological data, but the addition of molecular scaffold and combination of morphological and molecular data greatly improve the resolution of fossil nodes. We raise the monotypic Stephania subg. Botryodiscia to generic status and discover that three fossils previously assigned to Stephania should be members of Diploclisia. The Bayesian tip-dated tree recovered by removing homoplastic morphological characters with a Rescaled Consistency Index <0.25 has the highest stratigraphic fit and consequently generates more reasonable biogeographic reconstruction for Cissampelideae. Cissampelideae began to diversify in Asia in the latest Cretaceous and subsequently dispersed to South America around the Cretaceous-Palaeogene boundary. Two dispersal events from Asia to Africa occurred in the Early Eocene and the Late Eocene-Late Oligocene, respectively. These findings provide guidelines and practical methods for mitigating the effects of homoplastic morphological characters on fossil placements and Bayesian tip-dating, as well as insights into the past tropical floristic exchanges among different continents.

Fossils indicate marine dispersal in osteoglossid fishes, a classic example of continental vicariance

The separation of closely related terrestrial or freshwater species by vast marine barriers represents a biogeographical riddle. Such cases can provide evidence for vicariance, a process whereby ancient geological events like continental rifting divided ancestral geographical ranges. With an evolutionary history extending tens of millions of years, freshwater ecology, and distribution encompassing widely separated southern landmasses, osteoglossid bonytongue fishes are a textbook case of vicariance attributed to Mesozoic fragmentation of the Gondwanan supercontinent. Largely overlooked fossils complicate the clean narrative invoked for extant species by recording occurrences on additional continents and in marine settings. Here, we present a new total-evidence phylogenetic hypothesis for bonytongue fishes combined with quantitative models of range evolution and show that the last common ancestor of extant osteoglossids was likely marine, and that the group colonized freshwater settings at least four times when both extant and extinct lineages are considered. The correspondence between extant osteoglossid relationships and patterns of continental fragmentation therefore represents a striking example of biogeographical pseudocongruence. Contrary to arguments against vicariance hypotheses that rely only on temporal or phylogenetic evidence, these results provide direct palaeontological support for enhanced dispersal ability early in the history of a group with widely separated distributions in the modern day.

Historical climate change and vicariance events contributed to the intercontinental disjunct distribution pattern of ash species (Fraxinus, Oleaceae)

The Northern Hemisphere temperate forests exhibit a disjunct distributional pattern in Europe, North America, and East Asia. Here, to reveal the promoter of intercontinental disjunct distribution, Fraxinus was used as a model organism to integrate abundant fossil evidence with high-resolution phylogenies in a phytogeographic analysis. We constructed a robust phylogenetic tree using genomic data, reconstructed the geographic ancestral areas, and evaluated the effect of incorporating fossil information on the reconstructed biogeographic history. The phylogenetic relationships of Fraxinus were highly resolved and divided into seven clades. Fraxinus originated in western North America during Eocene, and six intercontinental dispersal events and five intercontinental vicariance events were occured. Results suggest that climate change and vicariance contributed to the intercontinental disjunct distribution pattern of Fraxinus. Moreover, results highlight the necessity of integrating phylogenetic relationship and fossil to improve the reliability of inferred biogeographic events and our understanding of the processes underlying disjunct distributions.

Phylotranscriptomic analyses reveal deep gene tree discordance in Camellia (Theaceae)

Gene tree discordance is a significant legacy of biological evolution. Multiple factors can result in incongruence among genes, such as introgression, incomplete lineage sorting (ILS), gene duplication or loss. Resolving the background of gene tree discordance is a critical way to uncover the process of species diversification. Camellia, the largest genus in Theaceae, has controversial taxonomy and systematics due in part to a complex evolutionary history. We used 60 transcriptomes of 55 species, which represented 15 sections of Camellia to investigate its phylogeny and the possible causes of gene tree discordance. We conducted gene tree discordance analysis based on 1,617 orthologous low-copy nuclear genes, primarily using coalescent species trees and polytomy tests to distinguish hard and soft conflict. A selective pressure analysis was also performed to assess the impact of selection on phylogenetic topology reconstruction. Our results detected different levels of gene tree discordance in the backbone of Camellia, and recovered rapid diversification as one of the possible causes of gene tree discordance. Furthermore, we confirmed that none of the currently proposed sections of Camellia was monophyletic. Comparisons among datasets partitioned under different selective pressure regimes showed that integrating all orthologous genes provided the best phylogenetic resolution of the species tree of Camellia. The findings of this study reveal rapid diversification as a major source of gene tree discordance in Camellia and will facilitate future investigation of reticulate relationships at the species level in this important plant genus.

The evolutionary history of bees in time and space

Bees are the most significant pollinators of flowering plants. This partnership began ca. 120 million years ago, but the uncertainty of how and when bees spread across the planet has greatly obscured investigations of this key mutualism. We present a novel analysis of bee biogeography using extensive new genomic and fossil data to demonstrate that bees originated in Western Gondwana (Africa and South America). Bees likely originated in the Early Cretaceous, shortly before the breakup of Western Gondwana, and the early evolution of any major bee lineage is associated with either the South American or African land masses. Subsequently, bees colonized northern continents via a complex history of vicariance and dispersal. The notable early absences from large landmasses, particularly in Australia and India, have important implications for understanding the assembly of local floras and diverse modes of pollination. How bees spread around the world from their hypothesized Southern Hemisphere origin parallels the histories of numerous flowering plant clades, providing an essential step to studying the evolution of angiosperm pollination syndromes in space and time.

Phylogenomics and historical biogeography of Hydrangeeae (Hydrangeaceae) elucidate the effects of geologic and climatic dynamics on diversification

Demonstrating the process of transregional biogeography and mechanisms underlying evolutionary radiations is crucial to understanding biological evolution. Here, we use Hydrangeeae (Hydrangeaceae), a tribe with a unique disjunct distribution and complex trait variations, using a solid phylogenetic framework, to investigate how geographical and climatic factors interact with functional traits to trigger plant evolutionary radiations. We constructed the first highly supported and dated phylogenetic framework using 79 protein-coding genes obtained from 81 plastomes, representing 63 species and all major clades, and found that most extant species originated from asynchronous diversification of two lineages undergoing repeated expansion and retraction, at middle and high latitudes of the Northern Hemisphere between East Asia and North America, during the Eocene to Pleistocene (driven by geologic and climatic dynamics). In accordance with these drivers, interactions of flora between central-eastern China and Japan occurred frequently after the Late Tertiary. We found that resource limitation and range fragmentation probably accelerated the diversification of Hydrangeeae, which supports the resource-use hypothesis. Our study sheds light on the evolutionary radiation and assembly of flora within East Asia, and the East Asian-North American disjunction, through integration of phylogenomic and biogeographic data with functional trait and ecological data.

Five long-distance dispersals shaped the major intercontinental disjunctions in Tectariaceae s.l. (Polypodiales, Polypodiopsida)

Intercontinental disjunct distributions can arise either from vicariance, from long-distance dispersal, or through extinction of an ancestral population with a broader distribution. Tectariaceae s.l., a clade of ferns in Polypodiales with ca. 300 species mainly distributed in the tropics and subtropics, provide an excellent opportunity to investigate global distribution patterns. Here, we assembled a dataset of eight plastid markers and one nuclear marker of 636 (92% increase of the earlier largest sampling) accessions representing ca. 210 species of all eight genera in Tectariaceae s.l. (Arthropteridaceae, Pteridryaceae, and Tectariaceae s.s.) and 35 species of other families of eupolypods Ⅰ. A new phylogeny is reconstructed to study the biogeography and trait-associated diversification. Our major results include: (1) a distinct lineage of Tectaria sister to the rest of the American Tectaria is identified; (2) Tectariaceae s.l., and the three families: Arthropteridaceae (Arthropteris), Pteridryaceae (Draconopteris, Malaifilix, Polydictyum, Pteridrys), and Tectariaceae s.s. (Hypoderris, Tectaria, and Triplophyllum), might have all originated in late Cretaceous; (3) only five intercontinental dispersals occurred in Pteridryaceae and Tectariaceae s.s. giving rise to their current intercontinental disjunction; (4) we provide the second evidence in ferns that a long-distance dispersal between Malesia and Americas during the Paleocene to Eocene led to the establishment/origin of a new genus (Draconopteris); and (5) diversification rate of each state of leaf dissection is different, and the lowest is in the simple-leaved taxa.

Phylogenomic analyses of Camellia support reticulate evolution among major clades

Camellia (Theaceae) is a morphologically highly diverse genus of flowering plants and includes many famous species with high economic value, and the phylogeny of this genus is not fully resolved. We used 95 transcriptomes from 87 Camellia species and identified 1481 low-copy genes to conduct a detailed analysis of the phylogeny of this genus according to various data-screening criteria. The results show that, very different from the two existing classification systems of Camellia, 87 species are grouped into 8 main clades and two independent species, and that all 8 clades except Clade 8 were strongly supported by almost all the coalescent or concatenated trees using different gene subsets. However, the relationships among these clades were weakly supported and different from analyses using different gene subsets; furthermore, they do not agree with the phylogeny from chloroplast genomes of Camellia. Additional analyses support reticulate evolution (probably resulting from introgression or hybridization) among some major Camellia lineages, providing explanation for extensive gene tree conflicts. Furthermore, we inferred that together with the formation of East Asian subtropical evergreen broad-leaved forests, Camellia underwent a radiative divergence of major clades at 23 ∼ 19 Ma in the late Miocene then had a subsequent species burst at 10 ∼ 5 Ma. Principal component and cluster analyses provides new insights into morphological changes underlying the evolution of Camellia and a reference to further clarify subgenus and sections of this genus. The comprehensive study here including a nuclear phylogeny and other analyses reveal the rapid evolutionary history of Camellia.

Out of Madagascar, repeatedly: The phylogenetics and biogeography of Dombeyoideae (Malvaceae s.l.)

Dispersals have been shown to be critical to the evolution of the long isolated but megadiverse flora of Madagascar and the surrounding islands of the western Indian Ocean, but we are just beginning to understand the directionality of these dispersals. With more than half of its species occurring in the western Indian Ocean region (WIOR), the paleotropical subfamily Dombeyoideae provides a particularly useful case study through which to better understand the biogeography of the WIOR, and yet its biogeography is poorly understood. Here we sampled six molecular markers from all 20 genera in the Dombeyoideae to reconstruct the most complete phylogeny to date for the subfamily. From this, divergence times, calibrated with three fossils (two dombeyoid, one malvoid), and ancestral range estimations were hypothesized. Biogeographic stochastic mapping (BSM) analyses on the maximum clade credibility tree were completed and compared to BSM analyses on 1,000 trees randomly sampled from the posterior distribution of trees resulting from the dating analysis. We found the Dombeyoideae crown node diverged ca. 53 million years ago out of a broad ancestral range involving all three major areas of its distribution: Madagascar, Africa, and Asia. The majority of diversification and dispersals in the subfamily occurred within the last ca. 10 million years, mostly from the Pliocene onwards. There were roughly five dispersals from Madagascar to Africa (and only one in reverse), at least six from Madagascar to surrounding islands of the WIOR (Mascarenes and Comoros), and one dispersal from Madagascar to Asia (and ca. 1 in reverse). Other long-distance dispersals included one from Africa to St. Helena and one from Africa to Australasia, both from within the most widespread clade, the Cheirolaena & allies clade, and one dispersal from Asia to Africa. Critically, the Dombeyoideae provide strong evidence for considering the island of Madagascar as a source for the colonization of continents, as well as the surrounding islands of the WIOR. Furthermore, narrow sympatry was a key process in the evolution of the subfamily, particularly in Madagascar and the Mascarenes.

Phylogenomic insights into the origin and evolutionary history of evergreen broadleaved forests in East Asia under Cenozoic climate change

The evergreen versus deciduous leaf habit is an important functional trait for adaptation of forest trees and has been hypothesized to be related to the evolutionary processes of the component species under paleoclimatic change, and potentially reflected in the dynamic history of evergreen broadleaved forests (EBLFs) in East Asia. However, knowledge about the shift of evergreen versus deciduous leaf with the impact of paleoclimatic change using genomic data remains rare. Here, we focus on the Litsea complex (Lauraceae), a key lineage with dominant species of EBLFs, to gain insights into how evergreen versus deciduous trait shifted, providing insights into the origin and historical dynamics of EBLFs in East Asia under Cenozoic climate change. We reconstructed a robust phylogeny of the Litsea complex using genome-wide single-nucleotide variants (SNVs) with eight clades resolved. Fossil-calibrated analyses, diversification rate shifts, ancestral habit, ecological niche modelling and climate niche reconstruction were employed to estimate its origin and diversification pattern. Taking into account studies on other plant lineages dominating EBLFs of East Asia, it was revealed that the prototype of EBLFs in East Asia probably emerged in the Early Eocene (55-50 million years ago [Ma]), facilitated by the greenhouse warming. As a response to the cooling and drying climate in the Middle to Late Eocene (48-38 Ma), deciduous habits were evolved in the dominant lineages of the EBLFs in East Asia. Up to the Early Miocene (23 Ma), the prevailing of East Asian monsoon increased the extreme seasonal precipitation and accelerated the emergence of evergreen habits of the dominant lineages, and ultimately shaped the vegetation resembling that of today.

Phylogenomics and the flowering plant tree of life

The advances accelerated by next-generation sequencing and long-read sequencing technologies continue to provide an impetus for plant phylogenetic study. In the past decade, a large number of phylogenetic studies adopting hundreds to thousands of genes across a wealth of clades have emerged and ushered plant phylogenetics and evolution into a new era. In the meantime, a roadmap for researchers when making decisions across different approaches for their phylogenomic research design is imminent. This review focuses on the utility of genomic data (from organelle genomes, to both reduced representation sequencing and whole-genome sequencing) in phylogenetic and evolutionary investigations, describes the baseline methodology of experimental and analytical procedures, and summarizes recent progress in flowering plant phylogenomics at the ordinal, familial, tribal, and lower levels. We also discuss the challenges, such as the adverse impact on orthology inference and phylogenetic reconstruction raised from systematic errors, and underlying biological factors, such as whole-genome duplication, hybridization/introgression, and incomplete lineage sorting, together suggesting that a bifurcating tree may not be the best model for the tree of life. Finally, we discuss promising avenues for future plant phylogenomic studies.

Pseudopollen in Camellia oleifera and its implications for pollination ecology and taxonomy

BACKGROUND AND AIMS

In 1997, Tsou described the special differentiation of the connective tissues of some species of Theaceae to produce single-celled powders with unique patterns called pseudopollen. The purpose of this study was to investigate the morphological structure of the pseudopollen of Camellia oleifera (Theaceae) and to study the morphology of pseudopollen in seven other Camellia species.

METHODS

Scanning electron microscopy, paraffin section, light microscopy, transmission electron microscopy, histochemistry.

KEY RESULT

C. oleifera pseudopollen was similar to normal pollen in macroscopic morphology but different microscopically. The normal pollen was starch-rich and yellow, with mostly reticulate exine ornamentation. In contrast, the pseudopollen was a white powder, single-celled and rich in protein, with parallel unbranched ridge lines on the outer wall, and originated from the parenchyma of the connective tissues. There are also differences in the micro-characteristics of normal and pseudopollen among different species in Camellia.

CONCLUSION

There are great differences in morphological structure between C. oleifera and other species in Camellia normal pollen and pseudopollen; these results may indicate that the pseudopollen can be used as a taxonomic basis for Camellia, and the macroscopic similarity between pseudopollen and pollen and histochemical characteristics of pseudopollen can be a pollination strategy.

Phylogeographic structure of Heteroplexis (Asteraceae), an endangered endemic genus in the limestone karst regions of southern China

Though the karst regions in south and southwest China are plant diversity hotspots, our understanding of the phylogeography and evolutionary history of the plants there remains limited. The genus Heteroplexis (Asteraceae) is one of the typical representative plants isolated by karst habitat islands, and is also an endangered and endemic plant to China. In this study, species-level phylogeographic analysis of the genus Heteroplexis was conducted using restriction site-associated DNA sequencing (RADseq). The genetic structure showed a clear phylogeographic structure consistent with the current species boundaries in the H. microcephala, H. incana, H. vernonioides, H. sericophylla, and H. impressinervia. The significant global (R = 0.37, P < 0.01) and regional (R = 0.650.95, P < 0.05) isolation by distance (IBD) signals among species indicate strong geographic isolation in the karst mountains, which may result in chronically restricted gene flow and increased genetic drift and differentiation. Furthermore, the phylogeographic structure of Heteroplexis suggested a southward migration since the last glacial period. Demographic analysis revealed the karst mountains as a refuge for Heteroplexis species. Finally, both Treemix and ABBA-BABA statistic detected significant historical gene flow between species. Significant historical gene flow and long-term stability of effective population size (Ne) together explain the high genome-wide genetic diversity among species (π = 0.05370.0838). However, the recent collapse of Ne, widespread inbreeding within populations, and restricted contemporary gene flow suggest that Heteroplexis species are probably facing a high risk of genetic diversity loss. Our results help to understand the evolutionary history of karst plants and guide conservation.

Phylogenomics Resolves the Phylogeny of Theaceae by Using Low-Copy and Multi-Copy Nuclear Gene Makers and Uncovers a Fast Radiation Event Contributing to Tea Plants Diversity

Simple Summary

The Theaceae includes more than 300 species of great morphological diversity and has immense economic, cultural, and ornamental values. However, the evolutionary history of this family remains elusive. We integrated 91 genomes and transcriptome datasets of Theaceae and successfully resolved the phylogeny of Theaceae including relatives of cultivated tea plants from both extensive low-copy and multi-copy nuclear gene markers. Bayes-based molecular dating revealed that the ancestor of the tea family originated slightly earlier than the K-Pg boundary (Mass extinction events including the extinction of dinosaurs) with early diversification of three tribes associated with the Early Eocene Climatic Optimum. Further speciation analysis suggested a sole significant diversification shift rate in the common ancestor of Camellia associated with the Mid-Miocene Climatic Optimum. Collectively, polyploidy events, and key morphological innovation characters, such as pericarp with seed coat hardening, could possibly contribute to the Theaceae species diversity.

Abstract

Tea is one of the three most popular nonalcoholic beverages globally and has extremely high economic and cultural value. Currently, the classification, taxonomy, and evolutionary history of the tea family are largely elusive, including phylogeny, divergence, speciation, and diversity. For understanding the evolutionary history and dynamics of species diversity in Theaceae, a robust phylogenetic framework based on 1785 low-copy and 79,103 multi-copy nuclear genes from 91 tea plant genomes and transcriptome datasets had been reconstructed. Our results maximumly supported that the tribes Stewartieae and Gordonieae are successive sister groups to the tribe Theeae from both coalescent and super matrix ML tree analyses. Moreover, in the most evolved tribe, Theeae, the monophyletic genera Pyrenaria, Apterosperma, and Polyspora are the successive sister groups of Camellia. We also yield a well-resolved relationship of Camellia, which contains the vast majority of Theaceae species richness. Molecular dating suggests that Theaceae originated in the late L-Cretaceous, with subsequent early radiation under the Early Eocene Climatic Optimal (EECO) for the three tribes. A diversification rate shift was detected in the common ancestors of Camellia with subsequent acceleration in speciation rate under the climate optimum in the early Miocene. These results provide a phylogenetic framework and new insights into factors that likely have contributed to the survival of Theaceae, especially a successful radiation event of genus Camellia members to subtropic/tropic regions. These novel findings will facilitate the efficient conservation and utilization of germplasm resources for breeding cultivated tea and oil-tea. Collectively, these results provide a foundation for further morphological and functional evolutionary analyses across Theaceae.

Phylotranscriptomics of Theaceae: generic-level relationships, reticulation and whole-genome duplication

BACKGROUND AND AIMS

Theaceae, with three tribes, nine genera and more than 200 species, are of great economic and ecological importance. Recent phylogenetic analyses based on plastomic data resolved the relationships among the three tribes and the intergeneric relationships within two of those tribes. However, generic-level relationships within the largest tribe, Theeae, were not fully resolved. The role of putative whole-genome duplication (WGD) events in the family and possible hybridization events among genera within Theeae also remain to be tested further.

METHODS

Transcriptomes or low-depth whole-genome sequencing of 57 species of Theaceae, as well as additional plastome sequence data, were generated. Using a dataset of low-copy nuclear genes, we reconstructed phylogenetic relationships using concatenated, species tree and phylogenetic network approaches. We further conducted molecular dating analyses and inferred possible WGD events by examining the distribution of the number of synonymous substitutions per synonymous site (Ks) for paralogues in each species. For plastid protein-coding sequences , phylogenies were reconstructed for comparison with the results obtained from analysis of the nuclear dataset.

RESULTS

Based on the 610 low-copy nuclear genes (858 606 bp in length) investigated, Stewartieae was resolved as sister to the other two tribes. Within Theeae, the Apterosperma-Laplacea clade grouped with Pyrenaria, leaving Camellia and Polyspora as sister. The estimated ages within Theaceae were largely consistent with previous studies based mainly on plastome data. Two reticulation events within Camellia and one between the common ancestor of Gordonia and Schima were found. All members of the tea family shared two WGD events, an older At-γ and a recent Ad-β; both events were also shared with the outgroups (Diapensiaceae, Pentaphylacaceae, Styracaceae and Symplocaceae).

CONCLUSIONS

Our analyses using low-copy nuclear genes improved understanding of phylogenetic relationships at the tribal and generic levels previously proposed based on plastome data, but the phylogenetic position of the Apterosperma-Laplacea clade needs more attention. There is no evidence for extensive intergeneric hybridization within Theeae or for a Theaceae-specific WGD event. Land bridges (e.g. the Bering land bridge) during the Late Oligocene may have permitted the intercontinental plant movements that facilitated the putative ancient introgression between the common ancestor of Gordonia and Schima.

Phylogenomics and biogeography of Catalpa (Bignoniaceae) reveal incomplete lineage sorting and three dispersal events

Catalpa Scop. (Bignoniaceae) is a small genus (8 spp.) of trees that is disjunctly distributed among eastern Asia, eastern United States, and the West Indies. Catalpa bears beautiful inflorescences and have been cultivated as important ornamental trees for landscaping, gardening, and timber. However, the phylogenetic relationships and biogeographic history of the genus have remained unresolved. In this study, we used a large genomic dataset that includes data from the chloroplast (plastomes), and nuclear genomes (ITS and 5,759 single-copy nuclear genes) to reconstruct phylogenetic relationship within Catalpa, test interspecific gene flow events within the genus, and infer its biogeographic history. Our phylogenetic results indicate that Catalpa is monophyletic containing two main clades, section Catalpa and section Macrocatalpa. Section Catalpa is further divided into three subclades. While most relationships are congruent between the chloroplast and nuclear datasets, the position of C. ovata differs, likely due to incomplete lineage sorting. Interspecific gene flow events include C. bungei s.s. with vectors of inheritance from C. duclouxii and C. fargesii, supporting a combination of these three species and recognizing a broadly circumscribed C. bungei s.l. Our biogeographic study suggests three main dispersal events, two of which occurred during the Oligocene. The first dispersal event occurred from southwestern North America and Mexico into the Greater Antilles giving rise to the ancestor of the section of Macrocatalpa. The second dispersal event also occurred from southwestern North America and Mexico, but led to central and northern North America, subsequently reaching China through the Bering land bridge, and also reaching Europe through the North Atlantic land bridge. The third dispersal event took place in the Miocene from China to North America and gave rise to a clade composed of C. bignonioides and C. speciosa. This study uses a phylogenomic approach and biogeographical methods to infer the evolutionary history of Catalpa, highlighting issues associated with gene tree discordance, and suggesting that incomplete lineage sorting likely played an important role in the evolutionary history of Catalpa.

Incorporating Topological and Age Uncertainty into Event-Based Biogeography of Sand Spiders Supports Paleo-Islands in Galapagos and Ancient Connections among Neotropical Dry Forests

Event-based biogeographic methods, such as dispersal-extinction-cladogenesis, have become increasingly popular for attempting to reconstruct the biogeographic history of organisms. Such methods employ distributional data of sampled species and a dated phylogenetic tree to estimate ancestral distribution ranges. Because the input tree is often a single consensus tree, uncertainty in topology and age estimates are rarely accounted for, even when they may affect the outcome of biogeographic estimates. Even when such uncertainties are taken into account for estimates of ancestral ranges, they are usually ignored when researchers compare competing biogeographic hypotheses. We explore the effect of incorporating this uncertainty in a biogeographic analysis of the 21 species of sand spiders (Sicariidae: Sicarius) from Neotropical xeric biomes, based on a total-evidence phylogeny including a complete sampling of the genus. Using a custom R script, we account for uncertainty in ages and topology by estimating ancestral ranges over a sample of trees from the posterior distribution of a Bayesian analysis, and for uncertainty in biogeographic estimates by using stochastic maps. This approach allows for counting biogeographic events such as dispersal among areas, counting lineages through time per area, and testing biogeographic hypotheses, while not overestimating the confidence in a single topology. Including uncertainty in ages indicates that Sicarius dispersed to the Galapagos Islands when the archipelago was formed by paleo-islands that are now submerged; model comparison strongly favors a scenario where dispersal took place before the current islands emerged. We also investigated past connections among currently disjunct Neotropical dry forests; failing to account for topological uncertainty underestimates possible connections among the Caatinga and Andean dry forests in favor of connections among Caatinga and Caribbean + Mesoamerican dry forests. Additionally, we find that biogeographic models including a founder-event speciation parameter (“+J”) are more prone to suffer from the overconfidence effects of estimating ancestral ranges using a single topology. This effect is alleviated by incorporating topological and age uncertainty while estimating stochastic maps, increasing the similarity in the inference of biogeographic events between models with or without a founder-event speciation parameter. We argue that incorporating phylogenetic uncertainty in biogeographic hypothesis-testing is valuable and should be a commonplace approach in the presence of rogue taxa or wide confidence intervals in age estimates, and especially when using models including founder-event speciation.