The evolution of multiple colour mechanisms is correlated with diversification in sunbirds (Nectariniidae).

How and why certain groups become speciose is a key question in evolutionary biology. Novel traits that enable diversification by opening new ecological niches are likely important mechanisms. However, ornamental traits can also promote diversification by opening up novel sensory niches and thereby creating novel inter-specific interactions. More specifically, ornamental colours may enable more precise and/or easier species recognition, and may act as key innovations by increasing the number of species-specific patterns and promoting diversification. While the influence of colouration on diversification is well-studied, the influence of the mechanisms that produce those colours (e.g. pigmentary, nanostructural) is less so, even though the ontogeny and evolution of these mechanisms differ. We estimated a new phylogenetic tree for 121 sunbird species and combined colour data of 106 species with a range of phylogenetic tools to test the hypothesis that the evolution of novel colour mechanisms increases diversification in sunbirds, one of the most colourful bird clades. Results suggest that (1) the evolution of novel colour mechanisms expands the visual sensory niche, increasing the number of achievable colours. (2) Structural colouration diverges more readily across the body than pigment-based colouration, enabling an increase in colour complexity. (3) Novel colour mechanisms might minimize trade-offs between natural and sexual selection such that colour can function both as camouflage and conspicuous signal. (4) Despite structural colours being more colourful and mobile, only melanin-based colouration is positively correlated with net diversification. Together, these findings explain why colour distances increase with increasing number of sympatric species, even though packing of colour space would predict otherwise.

Flagellar export apparatus and ATP synthetase: Homology evidenced by synteny predating the Last Universal Common Ancestor.

We report evidence further supporting homology between proteins in the F1 FO -ATP synthetase and the bacterial flagellar motor (BFM). BFM proteins FliH, FliI, and FliJ have been hypothesized to be homologous to FO -b + F1 -δ, F1 -α/ß, and F1 -γ, with similar structure and interactions. We conduct a further test by constructing a gene order dataset, examining the order of fliH, fliI, and fliJ genes across the phylogenetic breadth of flagellar and nonflagellar type 3 secretion systems, and comparing this to published surveys of gene order in the F1 FO -ATP synthetase, its N-ATPase relatives, and the bacterial/archaeal V- and A-type ATPases. Strikingly, the fliHIJ gene order was deeply conserved, with the few exceptions appearing derived, and exactly matching the widely conserved F-ATPase gene order atpFHAG, coding for subunits b-δ-α-γ. The V/A-type ATPases have a similar conserved gene order. Our results confirm homology between these systems, and suggest a rare case of synteny conserved over billions of years, predating the Last Universal Common Ancestor (LUCA).

A new method for determining ribosomal DNA copy number shows differences between Saccharomyces cerevisiae populations.

Ribosomal DNA genes (rDNA) encode the major ribosomal RNAs and in eukaryotes typically form tandem repeat arrays. Species have characteristic rDNA copy numbers, but there is substantial intra-species variation in copy number that results from frequent rDNA recombination. Copy number differences can have phenotypic consequences, however difficulties in quantifying copy number mean we lack a comprehensive understanding of how copy number evolves and the consequences. Here we present a genomic sequence read approach to estimate rDNA copy number based on modal coverage to help overcome limitations with existing mean coverage-based approaches. We validated our method using Saccharomyces cerevisiae strains with known rDNA copy numbers. Application of our pipeline to a global sample of S. cerevisiae isolates showed that different populations have different rDNA copy numbers. Our results demonstrate the utility of the modal coverage method, and highlight the high level of rDNA copy number variation within and between populations.

Novel Integrative Modeling of Molecules and Morphology across Evolutionary Timescales.

Evolutionary models account for either population- or species-level processes but usually not both. We introduce a new model, the FBD-MSC, which makes it possible for the first time to integrate both the genealogical and fossilization phenomena, by means of the multispecies coalescent (MSC) and the fossilized birth-death (FBD) processes. Using this model, we reconstruct the phylogeny representing all extant and many fossil Caninae, recovering both the relative and absolute time of speciation events. We quantify known inaccuracy issues with divergence time estimates using the popular strategy of concatenating molecular alignments and show that the FBD-MSC solves them. Our new integrative method and empirical results advance the paradigm and practice of probabilistic total evidence analyses in evolutionary biology.[Caninae; fossilized birth-death; molecular clock; multispecies coalescent; phylogenetics; species trees.].

A Unique Mitochondrial Gene Block Inversion in Antarctic Trematomin Fishes: A Cautionary Tale.

Many Antarctic notothenioid fishes have major rearrangements in their mitochondrial (mt) genomes. Here, we report the complete mt genomes of 3 trematomin notothenioids: the bald notothen (Trematomus (Pagothenia) borchgrevinki), the spotted notothen (T. nicolai), and the emerald notothen (T. bernacchii). The 3 mt genomes were sequenced using next-generation Illumina technology, and the assemblies verified by Sanger sequencing. When compared with the canonical mt gene order of the Antarctic silverfish (Pleuragramma antarctica), we found a large gene inversion in the 3 trematomin mt genomes that included tRNAIle, ND1, tRNALeu2, 16S, tRNAVal, 12S, tRNAPhe, and the control region. The trematomin mt genomes contained 3 intergenic spacers, which are thought to be the remnants of previous gene and control region duplications. All control regions included the characteristic conserved regulatory sequence motifs. Although short-read next-generation DNA sequencing technology has allowed the rapid and cost-effective sequencing of a large number of complete mt genomes, it is essential in all cases to verify the assembly in order to prevent the publication and use of erroneous data.

Trait-dependent dispersal in rails (Aves: Rallidae): Historical biogeography of a cosmopolitan bird clade.

The ability of lineages to disperse over evolutionary timescales may be influenced by the gain or loss of traits after adaptation to new ecological conditions. For example, rails (Aves: Rallidae) have many cases of flightless insular endemic species that presumably evolved after flying ancestors dispersed over large ocean barriers and became isolated. Nonetheless, the details of how flying and its loss have influenced the clade's historical biogeography are unknown, as is the importance of other predictors of dispersal such as the geographic distance between regions. Here, we used a dated phylogeny of 158 species of rails to compare trait-dependent and trait-independent biogeography models in BioGeoBEARS. We evaluated a probabilistic historical biogeographical model that allows geographic range and flight to co-evolve and influence dispersal ability on a phylogeny. The best-fitting dispersal model was a trait-dependent dispersal (DEC + j + x + t21 + m1) that accrued 85.2% of the corrected Akaike Information Criterion (AICc) model weight. The distance-dependence parameter, x was estimated at -0.54, ranging from -0.49 to -0.65 across models, suggesting that a doubling of dispersal distance results in an approximately 31% decrease in dispersal rate (2-0.54 = 0.69). The estimated rate of loss of flight (t21) was similar across all models (~0.029 loss events per lineage per million years). The multiplier on dispersal rate when a lineage is non-flying, m1, is estimated to be 0.38 under this model. Surprisingly, the estimate of m1 was not 0.0, probably because the loss of flight is so common in the rails that entire clades of flightless species are found in the data, forcing the model to attribute some dispersal to flightless lineages. These results indicate that long-distance dispersal over macroevolutionary timespans can be modelled, rather than simply attributed to chance, allowing support for different hypotheses to be quantified and limitations to be identified. Overall, by combining new analytical methods with a comprehensive phylogeny, we use a quantitative framework to show how traits influence dispersal capacity and eventually shape geographical distributions at a macroevolutionary scale.

Statistical Comparison of Trait-Dependent Biogeographical Models Indicates That Podocarpaceae Dispersal Is Influenced by Both Seed Cone Traits and Geographical Distance.

The ability of lineages to disperse long distances over evolutionary timescales may be influenced by the gain or loss of traits adapted to enhance local, ecological dispersal. For example, some species in the southern conifer family Podocarpaceae have fleshy cones that encourage bird dispersal, but it is unknown how this trait has influenced the clade's historical biogeography, or its importance compared with other predictors of dispersal such as the geographic distance between regions. We answer these questions quantitatively by using a dated phylogeny of 197 species of southern conifers (Podocarpaceae and their sister family Araucariaceae) to statistically compare standard, trait-independent biogeography models with new BioGeoBEARS models where an evolving trait can influence dispersal probability, and trait history, biogeographical history, and model parameters are jointly inferred. We validate the method with simulation-inference experiments. Comparing all models, those that include trait-dependent dispersal accrue 87.5% of the corrected Akaike Information Criterion (AICc) model weight. Averaged across all models, lineages with nonfleshy cones had a dispersal probability multiplier of 0.49 compared with lineages with fleshy cones. Distance is included as a predictor of dispersal in all credible models (100% model weight). However, models with changing geography earned only 22.0% of the model weight, and models submerging New Caledonia/New Zealand earned only 0.01%. The importance of traits and distance suggests that long-distance dispersal over macroevolutionary timespans should not be thought of as a highly unpredictable chance event. Instead, long-distance dispersal can be modeled, allowing statistical model comparison to quantify support for different hypotheses.

Pioneering polyploids: the impact of whole-genome duplication on biome shifting in New Zealand Coprosma (Rubiaceae) and Veronica (Plantaginaceae).

The role of whole-genome duplication (WGD) in facilitating shifts into novel biomes remains unknown. Focusing on two diverse woody plant groups in New Zealand, Coprosma (Rubiaceae) and Veronica (Plantaginaceae), we investigate how biome occupancy varies with ploidy level, and test the hypothesis that WGD increases the rate of biome shifting. Ploidy levels and biome occupancy (forest, open and alpine) were determined for indigenous species in both clades. The distribution of low-ploidy (Coprosma: 2x, Veronica: 6x) versus high-ploidy (Coprosma: 4-10x, Veronica: 12-18x) species across biomes was tested statistically. Estimation of the phylogenetic history of biome occupancy and WGD was performed using time-calibrated phylogenies and the R package BioGeoBEARS. Trait-dependent dispersal models were implemented to determine support for an increased rate of biome shifting among high-ploidy lineages. We find support for a greater than random portion of high-ploidy species occupying multiple biomes. We also find strong support for high-ploidy lineages showing a three- to eightfold increase in the rate of biome shifts. These results suggest that WGD promotes ecological expansion into new biomes.

Sodium-powered stators of the bacterial flagellar motor can generate torque in the presence of phenamil with mutations near the peptidoglycan-binding region.

The bacterial flagellar motor powers the rotation that propels the swimming bacteria. Rotational torque is generated by harnessing the flow of ions through ion channels known as stators which couple the energy from the ion gradient across the inner membrane to rotation of the rotor. Here, we used error-prone PCR to introduce single point mutations into the sodium-powered Vibrio alginolyticus/Escherichia coli chimeric stator PotB and selected for motors that exhibited motility in the presence of the sodium-channel inhibitor phenamil. We found single mutations that enable motility under phenamil occurred at two sites: (i) the transmembrane domain of PotB, corresponding to the TM region of the PomB stator from V. alginolyticus and (ii) near the peptidoglycan binding region that corresponds to the C-terminal region of the MotB stator from E. coli. Single cell rotation assays confirmed that individual flagellar motors could rotate in up to 100 µM phenamil. Using phylogenetic logistic regression, we found correlation between natural residue variation and ion source at positions corresponding to PotB F22Y, but not at other sites. Our results demonstrate that it is not only the pore region of the stator that moderates motility in the presence of ion-channel blockers.

Habitat preference modulates trans-oceanic dispersal in a terrestrial vertebrate.

The importance of long-distance dispersal (LDD) in shaping geographical distributions has been debated since the nineteenth century. In terrestrial vertebrates, LDD events across large water bodies are considered highly improbable, but organismal traits affecting dispersal capacity are generally not taken into account. Here, we focus on a recent lizard radiation and combine a summary-coalescent species tree based on 1225 exons with a probabilistic model that links dispersal capacity to an evolving trait, to investigate whether ecological specialization has influenced the probability of trans-oceanic dispersal. Cryptoblepharus species that occur in coastal habitats have on average dispersed 13 to 14 times more frequently than non-coastal species and coastal specialization has, therefore, led to an extraordinarily widespread distribution that includes multiple continents and distant island archipelagoes. Furthermore, their presence across the Pacific substantially predates the age of human colonization and we can explicitly reject the possibility that these patterns are solely shaped by human-mediated dispersal. Overall, by combining new analytical methods with a comprehensive phylogenomic dataset, we use a quantitative framework to show how coastal specialization can influence dispersal capacity and eventually shape geographical distributions at a macroevolutionary scale.

Dating the Species Network: Allopolyploidy and Repetitive DNA Evolution in American Daisies (Melampodium sect. Melampodium, Asteraceae).

Allopolyploidy has played an important role in the evolution of the flowering plants. Genome mergers are often accompanied by significant and rapid alterations of genome size and structure via chromosomal rearrangements and altered dynamics of tandem and dispersed repetitive DNA families. Recent developments in sequencing technologies and bioinformatic methods allow for a comprehensive investigation of the repetitive component of plant genomes. Interpretation of evolutionary dynamics following allopolyploidization requires both the knowledge of parentage and the age of origin of an allopolyploid. Whereas parentage is typically inferred from cytogenetic and phylogenetic data, age inference is hampered by the reticulate nature of the phylogenetic relationships. Treating subgenomes of allopolyploids as if they belonged to different species (i.e., no recombination among subgenomes) and applying cross-bracing (i.e., putting a constraint on the age difference of nodes pertaining to the same event), we can infer the age of allopolyploids within the framework of the multispecies coalescent within BEAST2. Together with a comprehensive characterization of the repetitive DNA fraction using the RepeatExplorer pipeline, we apply the dating approach in a group of closely related allopolyploids and their progenitor species in the plant genus Melampodium (Asteraceae). We dated the origin of both the allotetraploid, Melampodium strigosum, and its two allohexaploid derivatives, Melampodium pringlei and Melampodium sericeum, which share both parentage and the direction of the cross, to the Pleistocene ($<$1.4 Ma). Thus, Pleistocene climatic fluctuations may have triggered formation of allopolyploids possibly in short intervals, contributing to difficulties in inferring the precise temporal order of allopolyploid species divergence of M. sericeum and M. pringlei. The relatively recent origin of the allopolyploids likely played a role in the near-absence of major changes in the repetitive fraction of the polyploids' genomes. The repetitive elements most affected by the postpolyploidization changes represented retrotransposons of the Ty1-copia lineage Maximus and, to a lesser extent, also Athila elements of Ty3-gypsy family.

Taxon cycle predictions supported by model-based inference in Indo-Pacific trap-jaw ants (Hymenoptera: Formicidae: Odontomachus).

Nonequilibrium dynamics and non-neutral processes, such as trait-dependent dispersal, are often missing from quantitative island biogeography models despite their potential explanatory value. One of the most influential nonequilibrium models is the taxon cycle, but it has been difficult to test its validity as a general biogeographical framework. Here, we test predictions of the taxon cycle model using six expected phylogenetic patterns and a time-calibrated phylogeny of Indo-Pacific Odontomachus (Hymenoptera: Formicidae: Ponerinae), one of the ant genera that E.O. Wilson used when first proposing the hypothesis. We used model-based inference and a newly developed trait-dependent dispersal model to jointly estimate ancestral biogeography, ecology (habitat preferences for forest interiors, vs. "marginal" habitats, such as savannahs, shorelines, disturbed areas) and the linkage between ecology and dispersal rates. We found strong evidence that habitat shifts from forest interior to open and disturbed habitats increased macroevolutionary dispersal rate. In addition, lineages occupying open and disturbed habitats can give rise to both island endemics re-occupying only forest interiors and taxa that re-expand geographical ranges. The phylogenetic predictions outlined in this study can be used in future work to evaluate the relative weights of neutral (e.g., geographical distance and area) and non-neutral (e.g., trait-dependent dispersal) processes in historical biogeography and community ecology.

Recent origin and rapid speciation of Neotropical orchids in the world's richest plant biodiversity hotspot.

The Andean mountains of South America are the most species-rich biodiversity hotspot worldwide with c. 15% of the world's plant species, in only 1% of the world's land surface. Orchids are a key element of the Andean flora, and one of the most prominent components of the Neotropical epiphyte diversity, yet very little is known about their origin and diversification. We address this knowledge gap by inferring the biogeographical history and diversification dynamics of the two largest Neotropical orchid groups (Cymbidieae and Pleurothallidinae), using two unparalleled, densely sampled orchid phylogenies (including more than 400 newly generated DNA sequences), comparative phylogenetic methods, geological and biological datasets. We find that the majority of Andean orchid lineages only originated in the last 20-15 million yr. Andean lineages are derived from lowland Amazonian ancestors, with additional contributions from Central America and the Antilles. Species diversification is correlated with Andean orogeny, and multiple migrations and recolonizations across the Andes indicate that mountains do not constrain orchid dispersal over long timescales. Our study sheds new light on the timing and geography of a major Neotropical diversification, and suggests that mountain uplift promotes species diversification across all elevational zones.

Spatiotemporal Diversification of the True Frogs (Genus Rana): A Historical Framework for a Widely Studied Group of Model Organisms.

True frogs of the genus Rana are widely used as model organisms in studies of development, genetics, physiology, ecology, behavior, and evolution. Comparative studies among the more than 100 species of Rana rely on an understanding of the evolutionary history and patterns of diversification of the group. We estimate a well-resolved, time-calibrated phylogeny from sequences of six nuclear and three mitochondrial loci sampled from most species of Rana, and use that phylogeny to clarify the group's diversification and global biogeography. Our analyses consistently support an "Out of Asia" pattern with two independent dispersals of Rana from East Asia to North America via Beringian land bridges. The more species-rich lineage of New World Rana appears to have experienced a rapid radiation following its colonization of the New World, especially with its expansion into montane and tropical areas of Mexico, Central America, and South America. In contrast, Old World Rana exhibit different trajectories of diversification; diversification in the Old World began very slowly and later underwent a distinct increase in speciation rate around 29-18 Ma. Net diversification is associated with environmental changes and especially intensive tectonic movements along the Asian margin from the Oligocene to early Miocene. Our phylogeny further suggests that previous classifications were misled by morphological homoplasy and plesiomorphic color patterns, as well as a reliance primarily on mitochondrial genes. We provide a phylogenetic taxonomy based on analyses of multiple nuclear and mitochondrial gene loci. [Amphibians; biogeography; diversification rate; Holarctic; transcontinental dispersal.

Inferring node dates from tip dates in fossil Canidae: the importance of tree priors.

Tip-dating methods are becoming popular alternatives to traditional node calibration approaches for building time-scaled phylogenetic trees, but questions remain about their application to empirical datasets. We compared the performance of the most popular methods against a dated tree of fossil Canidae derived from previously published monographs. Using a canid morphology dataset, we performed tip-dating using BEAST v. 2.1.3 and MrBayes v. 3.2.5. We find that for key nodes (Canis, approx. 3.2 Ma, Caninae approx. 11.7 Ma) a non-mechanistic model using a uniform tree prior produces estimates that are unrealistically old (27.5, 38.9 Ma). Mechanistic models (incorporating lineage birth, death and sampling rates) estimate ages that are closely in line with prior research. We provide a discussion of these two families of models (mechanistic versus non-mechanistic) and their applicability to fossil datasets.

Primary endosymbiosis events date to the later Proterozoic with cross-calibrated phylogenetic dating of duplicated ATPase proteins.

Chloroplasts and mitochondria descended from bacterial ancestors, but the dating of these primary endosymbiosis events remains very uncertain, despite their importance for our understanding of the evolution of both bacteria and eukaryotes. All phylogenetic dating in the Proterozoic and before is difficult: Significant debates surround potential fossil calibration points based on the interpretation of the Precambrian microbial fossil record, and strict molecular clock methods cannot be expected to yield accurate dates over such vast timescales because of strong heterogeneity in rates. Even with more sophisticated relaxed-clock analyses, nodes that are distant from fossil calibrations will have a very high uncertainty in dating. However, endosymbiosis events and gene duplications provide some additional information that has never been exploited in dating; namely, that certain nodes on a gene tree must represent the same events, and thus must have the same or very similar dates, even if the exact date is uncertain. We devised techniques to exploit this information: cross-calibration, in which node date calibrations are reused across a phylogeny, and cross-bracing, in which node date calibrations are formally linked in a hierarchical Bayesian model. We apply these methods to proteins with ancient duplications that have remained associated and originated from plastid and mitochondrial endosymbionts: the α and ß subunits of ATP synthase and its relatives, and the elongation factor thermo unstable. The methods yield reductions in dating uncertainty of 14-26% while only using date calibrations derived from phylogenetically unambiguous Phanerozoic fossils of multicellular plants and animals. Our results suggest that primary plastid endosymbiosis occurred ∼900 Mya and mitochondrial endosymbiosis occurred ∼1,200 Mya.

Bayesian analysis of a morphological supermatrix sheds light on controversial fossil hominin relationships.

The phylogenetic relationships of several hominin species remain controversial. Two methodological issues contribute to the uncertainty-use of partial, inconsistent datasets and reliance on phylogenetic methods that are ill-suited to testing competing hypotheses. Here, we report a study designed to overcome these issues. We first compiled a supermatrix of craniodental characters for all widely accepted hominin species. We then took advantage of recently developed Bayesian methods for building trees of serially sampled tips to test among hypotheses that have been put forward in three of the most important current debates in hominin phylogenetics--the relationship between Australopithecus sediba and Homo, the taxonomic status of the Dmanisi hominins, and the place of the so-called hobbit fossils from Flores, Indonesia, in the hominin tree. Based on our results, several published hypotheses can be statistically rejected. For example, the data do not support the claim that Dmanisi hominins and all other early Homo specimens represent a single species, nor that the hobbit fossils are the remains of small-bodied modern humans, one of whom had Down syndrome. More broadly, our study provides a new baseline dataset for future work on hominin phylogeny and illustrates the promise of Bayesian approaches for understanding hominin phylogenetic relationships.

Model selection in historical biogeography reveals that founder-event speciation is a crucial process in Island Clades.

Founder-event speciation, where a rare jump dispersal event founds a new genetically isolated lineage, has long been considered crucial by many historical biogeographers, but its importance is disputed within the vicariance school. Probabilistic modeling of geographic range evolution creates the potential to test different biogeographical models against data using standard statistical model choice procedures, as long as multiple models are available. I re-implement the Dispersal-Extinction-Cladogenesis (DEC) model of LAGRANGE in the R package BioGeoBEARS, and modify it to create a new model, DEC + J, which adds founder-event speciation, the importance of which is governed by a new free parameter, [Formula: see text]. The identifiability of DEC and DEC + J is tested on data sets simulated under a wide range of macroevolutionary models where geography evolves jointly with lineage birth/death events. The results confirm that DEC and DEC + J are identifiable even though these models ignore the fact that molecular phylogenies are missing many cladogenesis and extinction events. The simulations also indicate that DEC will have substantially increased errors in ancestral range estimation and parameter inference when the true model includes + J. DEC and DEC + J are compared on 13 empirical data sets drawn from studies of island clades. Likelihood-ratio tests indicate that all clades reject DEC, and AICc model weights show large to overwhelming support for DEC + J, for the first time verifying the importance of founder-event speciation in island clades via statistical model choice. Under DEC + J, ancestral nodes are usually estimated to have ranges occupying only one island, rather than the widespread ancestors often favored by DEC. These results indicate that the assumptions of historical biogeography models can have large impacts on inference and require testing and comparison with statistical methods.

Climate refugia: joint inference from fossil records, species distribution models and phylogeography.

Climate refugia, locations where taxa survive periods of regionally adverse climate, are thought to be critical for maintaining biodiversity through the glacial-interglacial climate changes of the Quaternary. A critical research need is to better integrate and reconcile the three major lines of evidence used to infer the existence of past refugia - fossil records, species distribution models and phylogeographic surveys - in order to characterize the complex spatiotemporal trajectories of species and populations in and out of refugia. Here we review the complementary strengths, limitations and new advances for these three approaches. We provide case studies to illustrate their combined application, and point the way towards new opportunities for synthesizing these disparate lines of evidence. Case studies with European beech, Qinghai spruce and Douglas-fir illustrate how the combination of these three approaches successfully resolves complex species histories not attainable from any one approach. Promising new statistical techniques can capitalize on the strengths of each method and provide a robust quantitative reconstruction of species history. Studying past refugia can help identify contemporary refugia and clarify their conservation significance, in particular by elucidating the fine-scale processes and the particular geographic locations that buffer species against rapidly changing climate.

Treating fossils as terminal taxa in divergence time estimation reveals ancient vicariance patterns in the palpimanoid spiders.

Incorporation of fossils into biogeographic studies can have a profound effect on the conclusions that result, particularly when fossil ranges are nonoverlapping with extant ranges. This is the case in archaeid spiders, where there are known fossils from the Northern Hemisphere, yet all living members are restricted to the Southern Hemisphere. To better understand the biogeographic patterns of archaeid spiders and their palpimanoid relatives, we estimate a dated phylogeny using a relaxed clock on a combined molecular and morphological data set. Dating information is compared with treating the archaeid fossil taxa as both node calibrations and as noncontemporaneous terminal tips, both with and without additional calibration points. Estimation of ancestral biogeographic ranges is then performed, using likelihood and Bayesian methods to take into account uncertainty in phylogeny and in dating. We find that treating the fossils as terminal tips within a Bayesian framework, as opposed to dating the phylogeny based only on molecular data with the dates coming from node calibrations, removes the subjectivity involved in assigning priors, which has not been possible with previous methods. Our analyses suggest that the diversification of the northern and southern archaeid lineages was congruent with the breakup of Pangaea into Laurasia and Gondwanaland. This analysis provides a rare example, and perhaps the most strongly supported, where a dated phylogeny confirms a biogeographical hypothesis based on vicariance due to the breakup of the ancient continental plates.