skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Inferring the evolution of reproductive isolation in a lineage of fossil threespine stickleback, Gasterosteus doryssus
Darwin attributed the absence of species transitions in the fossil record to his hypothesis that speciation occurs within isolated habitat patches too geographically restricted to be captured by fossil sequences. Mayr's peripatric speciation model added that such speciation would be rapid, further explaining missing evidence of diversification. Indeed, Eldredge and Gould's original punctuated equilibrium model combined Darwin's conjecture, Mayr's model and 124 years of unsuccessfully sampling the fossil record for transitions. Observing such divergence, however, could illustrate the tempo and mode of evolution during early speciation. Here, we investigate peripatric divergence in a Miocene stickleback fish,Gasterosteus doryssus.This lineage appeared and, over approximately 8000 generations, evolved significant reduction of 12 of 16 traits related to armour, swimming and diet, relative to its ancestral population. This was greater morphological divergence than we observed between reproductively isolated, benthic-limnetic ecotypes of extantGasterosteus aculeatus. Therefore, we infer that reproductive isolation was evolving.However, local extinction ofG. doryssuslineages shows how young, isolated, speciating populations often disappear, supporting Darwin's explanation for missing evidence and revealing a mechanism behind morphological stasis. Extinction may also account for limited sustained divergence within the stickleback species complex and help reconcile speciation rate variation observed across time scales.  more » « less
Award ID(s):
2145830
PAR ID:
10558891
Author(s) / Creator(s):
; ; ; ; ; ; ;
Publisher / Repository:
Royal Society Publishing
Date Published:
Journal Name:
Proceedings of the Royal Society B: Biological Sciences
Volume:
291
Issue:
2021
ISSN:
0962-8452
Subject(s) / Keyword(s):
anagenetic speciation Gasterosteus aculeatus punctuated equilibrium time series ephemeral speciation ecological speciation
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Ryan Folk (Ed.)
    Phylogenetic divergence-time estimation has been revolutionized by two recent developments: 1) total-evidence dating (or "tip-dating") approaches that allow for the incorporation of fossils as tips in the analysis, with their phylogenetic and temporal relationships to the extant taxa inferred from the data and 2) the fossilized birth-death (FBD) class of tree models that capture the processes that produce the tree (speciation, extinction, and fossilization) and thus provide a coherent and biologically interpretable tree prior. To explore the behavior of these methods, we apply them to marattialean ferns, a group that was dominant in Carboniferous landscapes prior to declining to its modest extant diversity of slightly over 100 species. We show that tree models have a dramatic influence on estimates of both divergence times and topological relationships. This influence is driven by the strong, counter-intuitive informativeness of the uniform tree prior, and the inherent nonidentifiability of divergence-time models. In contrast to the strong influence of the tree models, we find minor effects of differing the morphological transition model or the morphological clock model. We compare the performance of a large pool of candidate models using a combination of posterior-predictive simulation and Bayes factors. Notably, an FBD model with epoch-specific speciation and extinction rates was strongly favored by Bayes factors. Our best-fitting model infers stem and crown divergences for the Marattiales in the mid-Devonian and Late Cretaceous, respectively, with elevated speciation rates in the Mississippian and elevated extinction rates in the Cisuralian leading to a peak diversity of ∼2800 species at the end of the Carboniferous, representing the heyday of the Psaroniaceae. This peak is followed by the rapid decline and ultimate extinction of the Psaroniaceae, with their descendants, the Marattiaceae, persisting at approximately stable levels of diversity until the present. This general diversification pattern appears to be insensitive to potential biases in the fossil record; despite the preponderance of available fossils being from Pennsylvanian coal balls, incorporating fossilization-rate variation does not improve model fit. In addition, by incorporating temporal data directly within the model and allowing for the inference of the phylogenetic position of the fossils, our study makes the surprising inference that the clade of extant Marattiales is relatively young, younger than any of the fossils historically thought to be congeneric with extant species. This result is a dramatic demonstration of the dangers of node-based approaches to divergence-time estimation, where the assignment of fossils to particular clades is made a priori (earlier node-based studies that constrained the minimum ages of extant genera based on these fossils resulted in much older age estimates than in our study) and of the utility of explicit models of morphological evolution and lineage diversification. [Bayesian model comparison; Carboniferous; divergence-time estimation; fossil record; fossilized birth–death; lineage diversification; Marattiales; models of morphological evolution; Psaronius; RevBayes.] 
    more » « less
  2. Folk, Ryan (Ed.)
    Abstract Phylogenetic divergence-time estimation has been revolutionized by two recent developments: 1) total-evidence dating (or "tip-dating") approaches that allow for the incorporation of fossils as tips in the analysis, with their phylogenetic and temporal relationships to the extant taxa inferred from the data and 2) the fossilized birth-death (FBD) class of tree models that capture the processes that produce the tree (speciation, extinction, and fossilization) and thus provide a coherent and biologically interpretable tree prior. To explore the behavior of these methods, we apply them to marattialean ferns, a group that was dominant in Carboniferous landscapes prior to declining to its modest extant diversity of slightly over 100 species. We show that tree models have a dramatic influence on estimates of both divergence times and topological relationships. This influence is driven by the strong, counter-intuitive informativeness of the uniform tree prior, and the inherent nonidentifiability of divergence-time models. In contrast to the strong influence of the tree models, we find minor effects of differing the morphological transition model or the morphological clock model. We compare the performance of a large pool of candidate models using a combination of posterior-predictive simulation and Bayes factors. Notably, an FBD model with epoch-specific speciation and extinction rates was strongly favored by Bayes factors. Our best-fitting model infers stem and crown divergences for the Marattiales in the mid-Devonian and Late Cretaceous, respectively, with elevated speciation rates in the Mississippian and elevated extinction rates in the Cisuralian leading to a peak diversity of $${\sim}$$2800 species at the end of the Carboniferous, representing the heyday of the Psaroniaceae. This peak is followed by the rapid decline and ultimate extinction of the Psaroniaceae, with their descendants, the Marattiaceae, persisting at approximately stable levels of diversity until the present. This general diversification pattern appears to be insensitive to potential biases in the fossil record; despite the preponderance of available fossils being from Pennsylvanian coal balls, incorporating fossilization-rate variation does not improve model fit. In addition, by incorporating temporal data directly within the model and allowing for the inference of the phylogenetic position of the fossils, our study makes the surprising inference that the clade of extant Marattiales is relatively young, younger than any of the fossils historically thought to be congeneric with extant species. This result is a dramatic demonstration of the dangers of node-based approaches to divergence-time estimation, where the assignment of fossils to particular clades is made a priori (earlier node-based studies that constrained the minimum ages of extant genera based on these fossils resulted in much older age estimates than in our study) and of the utility of explicit models of morphological evolution and lineage diversification. [Bayesian model comparison; Carboniferous; divergence-time estimation; fossil record; fossilized birth–death; lineage diversification; Marattiales; models of morphological evolution; Psaronius; RevBayes.] 
    more » « less
  3. Abstract Allometric scaling describes the relationship of trait size to body size within and among taxa. The slope of the population‐level regression of trait size against body size (i.e. static allometry) is typically invariant among closely related populations and species. Such invariance is commonly interpreted to reflect a combination of developmental and selective constraints that delimit a phenotypic space into which evolution could proceed most easily. Thus, understanding how allometric relationships do eventually evolve is important to understanding phenotypic diversification. In a lineage of fossil Threespine Stickleback (Gasterosteus doryssus), we investigated the evolvability of static allometric slopes for nine traits (five armour and four non‐armour) that evolved significant trait differences across 10 samples over 8500 years. The armour traits showed weak static allometric relationships and a mismatch between those slopes and observed evolution. This suggests that observed evolution in these traits was not constrained by relationships with body size, perhaps because prior, repeated adaptation to freshwater habitats by Threespine Stickleback had generated strong selection to break constraint. In contrast, for non‐armour traits, we found stronger allometric relationships. Those allometric slopes did evolve on short time scales. However, those changes were small and fluctuating and the slopes remained strong predictors of the evolutionary trajectory of trait means over time (i.e. evolutionary allometry), supporting the hypothesis of allometry as constraint. 
    more » « less
  4. Living turtles are characterized by extraordinarily low species diversity given their age. The clade’s extensive fossil record indicates that climate and biogeography may have played important roles in determining their diversity. We investigated this hypothesis by collecting a molecular dataset for 591 individual turtles that, together, represent 80% of all turtle species, including representatives of all families and 98% of genera, and used it to jointly estimate phylogeny and divergence times. We found that the turtle tree is characterized by relatively constant diversification (speciation minus extinction) punctuated by a single threefold increase. We also found that this shift is temporally and geographically associated with newly emerged continental margins that appeared during the Eocene−Oligocene transition about 30 million years before present. In apparent contrast, the fossil record from this time period contains evidence for a major, but regional, extinction event. These seemingly discordant findings appear to be driven by a common global process: global cooling and drying at the time of the Eocene−Oligocene transition. This climatic shift led to aridification that drove extinctions in important fossil-bearing areas, while simultaneously exposing new continental margin habitat that subsequently allowed for a burst of speciation associated with these newly exploitable ecological opportunities. 
    more » « less
  5. Abstract Species delimitation is an imperative first step toward understanding Earth's biodiversity, yet what constitutes a species and the relative importance of the various processes by which new species arise continue to be debatable. Species delimitation in spiders has traditionally used morphological characters; however, certain mygalomorph spiders exhibit morphological homogeneity despite long periods of population‐level isolation, absence of gene flow, and consequent high degrees of molecular divergence. Studies have shown strong geographic structuring and significant genetic divergence among several species complexes within the trapdoor spider genusAptostichus, most of which are restricted to the California Floristic Province (CAFP) biodiversity hotspot. Specifically, theAptostichus icenogleicomplex, which comprises the three sibling species,A. barackobamai,A. isabella, andA. icenoglei, exhibits evidence of cryptic mitochondrial DNA diversity throughout their ranges in Northern, Central, and Southern California. Our study aimed to explicitly test species hypotheses within this assemblage by implementing a cohesion species‐based approach. We used genomic‐scale data (ultraconserved elements, UCEs) to first evaluate genetic exchangeability and then assessed ecological interchangeability of genetic lineages. Biogeographical analysis was used to assess the likelihood of dispersal versus vicariance events that may have influenced speciation pattern and process across the CAFP's complex geologic and topographic landscape. Considering the lack of congruence across data types and analyses, we take a more conservative approach by retaining species boundaries withinA. icenoglei. 
    more » « less