Historical contingency and ecological determinism interact to prime speciation in sticklebacks, Gasterosteus

Historical contingency and determinism are often cast as opposing paradigms under which evolutionary diversification operates. It may be, however, that both factors act together to promote evolutionary divergence, although there are few examples of such interaction in nature. We tested phylogenetic predictions of an explicit historical model of divergence (double invasions of freshwater by marine ancestors) in sympatric species of three-spined sticklebacks (Gasterosteus aculeatus) where determinism has been implicated as an important factor driving evolutionary novelty. Microsatellite DNA variation at six loci revealed relatively low genetic variation in freshwater populations, supporting the hypothesis that they were derived by colonization of freshwater by more diverse marine ancestors. Phylogenetic and genetic distance analyses suggested that pairs of sympatric species have evolved multiple times, further implicating determinism as a factor in speciation. Our data also supported predictions based on the hypothesis that the evolution of sympatric species was contingent upon 'double invasions' of postglacial lakes by ancestral marine sticklebacks. Sympatric sticklebacks, therefore, provide an example of adaptive radiation by determinism contingent upon historical conditions promoting unique ecological interactions, and illustrate how contingency and determinism may interact to generate geographical variation in species diversity     

An integrative method for delimiting cohesion species: finding the population-species interface in a group of Californian trapdoor spiders with extreme genetic divergence and geographic structuring

Here we present an objective, repeatable approach to delineating species when populations are divergent and highly structured geographically using the Californian trapdoor spider species complex Aptostichus atomarius Simon as a model system. This system is particularly difficult because under strict criteria of geographical concordance coupled with estimates of genetic divergence, an unrealistic number of population lineages would qualify as species (20 to 60). Our novel phylogeographic approach, which is generally applicable but particularly relevant to highly structured systems, uses genealogical exclusivity to establish a topological framework to examine lineages for genetic and ecological exchangeability in an effort to delimit cohesion species. Both qualitative assessments of habitat and niche-based distribution modeling are employed to evaluate selective regime and ecological interchangeability among genetic lineages; adaptive divergence among populations is weighted more heavily than simple geographical concordance. Based on these analyses we conclude that five cohesion species should be recognized, three of which are new to science.     

Nuclear gene genealogies reveal historical,demographic and selective factors associated with speciation in field crickets

Population genetics theory predicts that genetic drift should eliminate shared polymorphism, leading to monophyly or exclusivity of populations, when the elapsed time between lineage-splitting events is large relative to effective population size. We examined patterns of nucleotide variation in introns at four nuclear loci to relate processes affecting the history of genes to patterns of divergence among natural populations and species. Ancestral polymorphisms were shared among three recognized species, Gryllus firmus, G. pennsylvanicus, and G. ovisopis, and genealogical patterns suggest that successive speciation events occurred recently and rapidly relative to effective population size. High levels of shared polymorphism among these morphologically, behaviorally, and ecologically distinct species indicate that only a small fraction of the genome needs to become differentiated for speciation to occur. Among the four nuclear gene loci there was a 10-fold range in nucleotide diversity, and patterns of polymorphism and divergence suggest that natural selection has acted to maintain or eliminate variation at some loci. While nuclear gene genealogies may have limited applications in phylogeography or other approaches dependent on population monophyly, they provide important insights into the historical, demographic, and selective forces that shape speciation.     

Estimating a geographically explicit model of population divergence

Patterns of genetic variation can provide valuable insights for deciphering the relative roles of different evolutionary processes in species differentiation. However, population-genetic models for studying divergence in geographically structured species are generally lacking. Since these are the biogeographic settings where genetic drift is expected to predominate, not only are population-genetic tests of hypotheses in geographically structured species constrained, but generalizations about the evolutionary processes that promote species divergence may also be potentially biased. Here we estimate a population-divergence model in montane grasshoppers from the sky islands of the Rocky Mountains. Because this region was directly impacted by Pleistocene glaciation, both the displacement into glacial refugia and recolonization of montane habitats may contribute to differentiation. Building on the tradition of using information from the genealogical relationships of alleles to infer the geography of divergence, here the additional consideration of the process of gene-lineage sorting is used to obtain a quantitative estimate of population relationships and historical associations (i.e., a population tree) from the gene trees of five anonymous nuclear loci and one mitochondrial locus in the broadly distributed species Melanoplus oregonensis. Three different approaches are used to estimate a model of population divergence; this comparison allows us to evaluate specific methodological assumptions that influence the estimated history of divergence. A model of population divergence was identified that significantly fits the data better compared to the other approaches, based on per-site likelihood scores of the multiple loci, and that provides clues about how divergence proceeded in M. oregonensis during the dynamic Pleistocene. Unlike the approaches that either considered only the most recent coalescence (i.e., information from a single individual per population) or did not consider the pattern of coalescence in the gene genealogies, the population-divergence model that best fits the data was estimated by considering the pattern of gene lineage coalescence across multiple individuals, as well as loci. These results indicate that sampling of multiple individuals per population is critical to obtaining an accurate estimate of the history of divergence so that the signal of common ancestry can be separated from the confounding influence of gene flow-even though estimates suggest that gene flow is not a predominant factor structuring patterns of genetic variation across these sky island populations. They also suggest that the gene genealogies contain information about population relationships, despite the lack of complete sorting of gene lineages. What emerges from the analyses is a model of population divergence that incorporates both contemporary distributions and historical associations, and shows a latitudinal and regional structuring of populations reminiscent of population displacements into multiple glacial refugia. Because the population-divergence model itself is built upon the specific events shaping the history of M. oregonensis, it provides a framework for estimating additional population-genetic parameters relevant to understanding the processes governing differentiation in geographically structured species and avoids the problems of relying on overly simplified and inaccurate divergence models. The utility of these approaches, as well as the caveats and future improvements, for estimating population relationships and historical associations relevant to genetic analyses of geographically structured species are discussed.     

Molecular phylogenetics at the population/species interface in cave spiders of the southern Appalachians (Araneae:Nesticidae:Nesticus)

This paper focuses on the relationship between population genetic structure and speciation mechanisms in a monophyletic species group of Appalachian cave spiders (Nesticus). Using mtDNA sequence data gathered from 256 individuals, I analyzed patterns of genetic variation within and between populations for three pairs of closely related sister species. Each sister-pair comparison involves taxa with differing distributional and ecological attributes; if these ecological attributes are reflected in basic demographic differences, then speciation might proceed differently across these sister taxa comparisons. Both frequency-based and gene tree analyses reveal that the genetic structure of the Nesticus species studied is characterized by similar and essentially complete population subdivision, regardless of differences in general ecology. These findings contrast with results of prior genetic studies of cave-dwelling arthropods that have typically revealed variation in population structure corresponding to differences in general ecology. Species fragmentation through both extrinsic and intrinsic evolutionary forces has resulted in discrete, perhaps independent, populations within morphologically defined species. Large sequence divergence values observed between populations suggest that this independence may extend well into the past. These patterns of mtDNA genealogical structure and divergence imply that species as morphological lineages are currently more inclusive than basal evolutionary or phylogenetic units, a suggestion that has important implications for the study of speciation mechanisms.      

Geographical and Ecological Drivers of Mitonuclear Genetic Divergence in a Mediterranean Grasshopper

The study of the neutral and/or selective processes driving genetic variation in natural populations is central to determine the evolutionary history of species and lineages and understand how they interact with different historical and contemporary components of landscape heterogeneity. Here, we combine nuclear and mitochondrial data to study the processes shaping genetic divergence in the Mediterranean esparto grasshopper (Ramburiella hispanica). Our analyses revealed the presence of three main lineages, two in Europe that split in the Early-Middle Pleistocene and one in North Africa that diverged from the two European ones after the Messinian. Lineage-specific potential distribution models and tests of environmental niche differentiation suggest that the phylogeographic structure of the species was driven by allopatric divergence due to the re-opening of the Gibraltar strait at the end of the Messinian (Europe–Africa split) and population fragmentation in geographically isolated Pleistocene climatic refugia (European split). Although we found no evidence for environment as an important driver of genetic divergence at the onset of lineage formation, our analyses considering the spatial distribution of populations and different aspects of landscape composition suggest that genetic differentiation at mitochondrial loci was largely explained by environmental dissimilarity, whereas resistance-based estimates of geographical distance were the only predictors of genetic differentiation at nuclear markers. Overall, our study shows that although historical factors have largely shaped concordant range-wide patterns of mitonuclear genetic structure in the esparto grasshopper, different contemporary processes (neutral gene flow vs. environmental-based selection) seem to be governing the spatial distribution of genetic variation in the two genomes.     

Interspecific phylogeography of the Stator limbatus species complex: the geographic context of speciation and specialization

Diversification in phytophagous insects is often attributed to a propensity toward specialization and to a tendency for speciation to be associated with host-shifts. Phylogenetic analysis revealed a sister relationship between the generalist Stator limbatus and the specialist host-shifted Stator beali, providing a system to examine the genealogical and geographic origins of the main processes involved in this diversification: host-shifts, specialization, and reproductive isolation. We examine the interspecific phylogeographic relationships between these species using mitochondrial DNA sequence data. S. beali is derived within S. limbatus, rendering the latter paraphyletic and suggesting a budding process of speciation. The inherent polarity in this genealogical pattern indicates that the specialist habit, clumping oviposition behavior, and distinct genitalia of S. beali are all derived from the ancestral S. limbatus. The phylogeography of S. limbatus also shows strong geographic structure with divergences corresponding to known biogeographic boundaries, indicating that this evolutionary signal has not been erased by the vagaries of history. However, the derivation of S. beali and the evolution of reproductive isolation between the two species does not correspond to these known biogeographic boundaries, as S. beali and its sister clade of S. limbatus are restricted to the same geographic province. The geographic proximity of diversification combined with a divergence time estimated at the beginning of the Pleistocene indicates that speciation likely occurred very rapidly, although further genetic and ecological work is necessary to examine the mode of speciation. This study provides the historical context for ongoing evolutionary, ecological, and quantitative genetic research on the divergence in diet breadth between these species.     

Geographical versus ecological isolation of closely related black flies (Diptera: Simuliidae) inferred from phylogeny, geography, and ecology

To investigate patterns of geographical and ecological separation among morphologically similar, closely related species of black flies, we integrated ecological, geographical, and phylogenetic information, based on multiple gene sequences, for 12 species in the subgenus Gomphostilbia in Thailand. Molecular characters supported the monophyly of the Simulium ceylonicum species group, but not of the Simulium batoense species group, suggesting that revisionary work is needed for the latter. Both ecological and geographical isolation of similar taxa were revealed. Stream velocity and altitude were among the principal ecological factors differing between closely related species. Most closely related species in the subgenus Gomphostilbia overlap geographically, suggesting the possibility of sympatric speciation driven by ecological divergence. Geographical isolation via dispersal also might have contributed to species divergence, while Pleistocene climate changes possibly influenced population genetic structure, demographic history, and speciation of some members of the subgenus.     

Hierarchical comparative analysis of genetic and genitalic geographical structure: testing patterns of male and female genital evolution in the scarab beetle Phyllophaga hirticula (Coleoptera: Scarabaeidae)

It is generally accepted that genitalia are among the fastest evolving characters in insects and that selection on these structures may increase speciation rates in groups with polygamous mating systems. If selection is causing genitalic divergence between or among populations of a species, one prediction is that geographical structure of genitalic morphology would be in place before genetic structure of a rapidly evolving neutral marker. The current study tests this hypothesis in the geographically widespread scarab beetle Phyllophaga hirticula by evaluating whether standing variation in male and female genitalia is more or less geographically structured than a mitochondrial genetic marker. Geographical structure of mitochondrial (mt)DNA and male and female genitalic shape were analysed using analysis of variance, multivariate analysis of variance, Mantel tests, and tests of spatial autocorrelation. The results show that, although female genitalia are more geographically structured than mtDNA, male genitalia are not. This pattern suggests that selection on female genitalic variation may be causing divergence of these structures among populations.  © 2009 The Linnean Society of London, Biological Journal of the Linnean Society , 2009, 96 , 135–149.     

Testing paleolimnological predictions with molecular data: the origins of Holarctic Eubosmina

Zooplankton of the family Bosminidae have a unique paleolimnological record in many Holarctic lakes that provides a near continuous record of morphological change for thousands of years. If this morphological change could be interpreted reliably, then a rarely achieved direct observation of macroevolution would be feasible. We tested paleolimnological predictions derived from morphological variation found in the genus Eubosmina using mtDNA and nuclear DNA sequence variation from geographically distant Holarctic sites. The mtDNA and nDNA trees were congruent but genetic divergence was inversely associated with morphological divergence. The three most genetically divergent groups belonged to Eubosmina longispina, whose phylogeography and genetic divergence was consistent with glacial vicariance. The genetic evidence also supported the hypothesis that at least two Nearctic species were recent European introductions. Finally, the genetic evidence was consistent with paleolimnology in the finding of several proposed species undergoing rapid morphological evolution and being post-glacially derived from European E. longispina. The results suggested that lacustrine bosminids are susceptible to geographic speciation processes, and that morphological interpretation of diversity in paleolimnology can be markedly improved by genetic studies.     

Diversification in the northern neotropics: mitochondrial and nuclear DNA phylogeography of the iguana Ctenosaura pectinata and related species

While Quaternary climatic changes are considered by some to have been a major factor promoting speciation within the neotropics, others suggest that much of the neotropical species diversity originated before the Pleistocene. Using mitochondrial and nuclear sequence data, we evaluate the relative importance of Pleistocene and pre-Pleistocene events within the evolutionary history of the Mexican iguana Ctenosaura pectinata , and related species. Results support the existence of cryptic lineages with strong mitochondrial divergence (> 4%) among them. Some of these lineages form zones of secondary contact, with one of them hybridizing with C. hemilopha . Evolutionary network analyses reveal the oldest populations of C. pectinata to be those of the northern and southern Mexican coastal regions. Inland and mid-latitudinal coastal populations are younger in age as a consequence of a history of local extinction within these regions followed by re-colonization. Estimated divergence times suggest that C. pectinata originated during the Pliocene, whereas geographically distinct mitochondrial DNA lineages first started to diverge during the Pliocene, with subsequent divergence continuing through the Pleistocene. Our results highlight the influence of both Pliocene and Pleistocene events in shaping the geographical distribution of genetic variation within neotropical lowland organisms. Areas of high genetic diversity in southern Mexico were detected, this finding plus the high levels of genetic diversity within C. pectinata , have implications for the conservation of this threatened species.     

Evolutionary and demographic history of the Californian scrub white oak species complex: an integrative approach

Understanding the factors promoting species formation is a major task in evolutionary research. Here, we employ an integrative approach to study the evolutionary history of the Californian scrub white oak species complex (genus Quercus). To infer the relative importance of geographical isolation and ecological divergence in driving the speciation process, we (i) analysed inter‐ and intraspecific patterns of genetic differentiation and employed an approximate Bayesian computation (ABC) framework to evaluate different plausible scenarios of species divergence. In a second step, we (ii) linked the inferred divergence pathways with current and past species distribution models (SDMs) and (iii) tested for niche differentiation and phylogenetic niche conservatism across taxa. ABC analyses showed that the most plausible scenario is the one considering the divergence of two main lineages followed by a more recent pulse of speciation. Genotypic data in conjunction with SDMs and niche differentiation analyses support that different factors (geography vs. environment) and modes of speciation (parapatry, allopatry and maybe sympatry) have played a role in the divergence process within this complex. We found no significant relationship between genetic differentiation and niche overlap, which probably reflects niche lability and/or that multiple factors, have contributed to speciation. Our study shows that different mechanisms can drive divergence even among closely related taxa representing early stages of species formation and exemplifies the importance of adopting integrative approaches to get a better understanding of the speciation process.     

Cryptic diversity,diversification and vicariance in two species complexes of Tomocerus (Collembola,Tomoceridae) from China

Springtails (Collembola) are a group of arthropods that are found in terrestrial ecosystems throughout the world. Two species complexes, Tomocerus ocreatus and T. nigrus complexes, are widely distributed in the southern and northern parts of eastern China, respectively. There is a poor understanding of the species diversity within these complexes and of the factors affecting their diversification and dispersal. Species delimitation using a general mixed Yule coalescent model and a Bayesian multilocus approach recognized 22 DNA‐based species. This supports the presence of extensive cryptic diversity in species that are geographically widespread. In addition to genetic differences, we discovered corresponding morphological differences in jumping organs among the major clades. Analyses of divergence times and historical biogeographical processes revealed that ocreatus and nigrus complexes originated in southern and northern China, respectively. We estimated their divergence at 27.8–44.9 Mya during the Eocene–Oligocene, at the time when the transmeridional Qinling–Dabie Mountains uplifted and formed the north–south geographical boundary of eastern China. Diversification analyses suggest that the subsequent orogenesis of the Qinghai–Tibetan Plateau in western China had little impact on divergences within the two species complexes so that they maintained their geographical patterns from the Paleogene to the present day. Our findings also point to a potentially important influence of the Qinling–Dabie Mountains on patterns of animal speciation and distribution in China.     

Evolutionary patterns and genetic structure in localized and widespread species in the Stylidium caricifolium complex (Stylidiaceae)

The Stylidium caricifolium (Stylidiaceae) complex consists of seven currently recognized species and a taxon of putative hybrid origin endemic to southwest Western Australia. These taxa vary in geographical distribution from widespread, extending over a range of 500 km, to extremely localized, covering a range of only 0.5 km. Patterns of allozyme variation were investigated in 61 populations covering all taxa and two closely related species. Measures of genetic diversity were consistently lower and in some cases significantly lower in four rare and geographically restricted taxa compared with their widespread relatives. In contrast, genetic diversity in two other localized taxa was comparable or higher than in the widespread taxa. The level of divergence among populations was moderate to high, with a significant trend of higher F(ST) values for the widespread species to lower values for the geographically restricted and rare taxa. Phylogenetic relationships and levels of divergence indicate that most taxa are probably relictual rather than recently evolved. Geographical localization and rarity in this complex can be attributed to a range of factors associated with habitat specificity, historical and ecological processes that characterize the southwest region, and mode of origin.     

Idiosyncratic responses to drivers of genetic differentiation in the complex landscapes of Isthmian Central America

Isthmian Central America (ICA) is one of the most biodiverse regions in the world, hosting an exceptionally high number of species per unit area. ICA was formed <25 million years ago and, consequently, its biotic assemblage is relatively young and derived from both colonization and in situ diversification. Despite intensive taxonomic work on the local fauna, the potential forces driving genetic divergences and ultimately speciation in ICA remain poorly studied. Here, we used a landscape genetics approach to test whether isolation by distance, topography, habitat suitability, or environment drive the genetic diversity of the regional frog assemblage. To this end, we combined data on landscape features and mitochondrial DNA sequence variation for nine codistributed amphibian species with disparate life histories. In five species, we found that at least one of the factors tested explained patterns of genetic divergence. However, rather than finding a general pattern, our results revealed idiosyncratic responses to historical and ecological processes, indicating that intrinsic life-history characteristics may determine the effect of different drivers of isolation on genetic divergence in ICA. Our work also suggests that the convergence of several factors promoting isolation among populations over a heterogeneous landscape might maximize genetic differentiation, despite short geographical distances. In conclusion, abiotic factors and geographical features have differentially affected the genetic diversity across the regional frog assemblage. Much more complex models (i.e., considering multiple drivers), beyond simple vicariance of Caribbean and Pacific lineages, are needed to better understand the evolutionary history of ICA’s diverse biotas.Subject terms: Evolutionary genetics, Ecological genetics, Biogeography     

Phylogeography and divergence date estimates of a lichen species complex with a disjunct distribution pattern

Disjunct species distributions may result from a combination of geologic events and long-distance dispersal. The foliose lichen species complex Leptogium furfuraceum-L. pseudofurfuraceum has an intercontinental disjunction pattern. Populations of this species complex are found in western North America, southern South America, Africa, and southern Europe. We conducted a phylogenetic study to reconstruct the biogeographic history of this species complex using two ribosomal genes (ITS and LSU) and a protein-coding gene (partial RPB2). Results indicated that the complex comprises four geographically restricted genetic lineages. A sister relationship was found between populations from the same hemispheres, incongruent with previous data derived from morphological characteristics and geographical classification schemes. Incorporating Bayesian ancestral area reconstruction and Bayesian divergence time estimation, we proposed an evolutionary hypothesis for the species complex. The results suggested that processes of biotic expansion via transoceanic dispersal were responsible for the species divergence and distribution patterns observed today. This study also expands the view that cryptic speciation is not a rare phenomenon among fungi and lichens.     

Phenotypic population divergence in terrestrial vertebrates at macro scales

Phenotypic divergence between populations, i.e. how much phenotypes within a species vary geographically, is critical to many aspects of ecology and evolution, including eco-geographical trends, speciation and coexistence. Yet, the variation of divergence across species with different ecologies and distributions and the relative role of adaptive causes remains little understood. We predict that genetic control vs. phenotypic plasticity of traits, geographical distance and (assuming adaptation) environmental differences should explain much of the phenotypic variability between populations. We tested these predictions with body sizes of 1447 populations in 98 terrestrial vertebrate species. Population phenotypic variability differs strongly across species, and divergence increases with increasing levels of clade-typical phenotypic plasticity, the area covered by populations and body size. Geographical distance and environmental dissimilarity are similarly important predictors of divergence within species, highlighting a potential role for biotic and environmental conditions. Increased availability of phylogeographical and ecological data should facilitate further understanding of population divergence drivers at broad scales.     

Population differentiation and restricted gene flow in Spanish crossbills: not isolation-by-distance but isolation-by-ecology

Divergent selection stemming from environmental variation may induce local adaptation and ecological speciation whereas gene flow might have a homogenizing effect. Gene flow among populations using different environments can be reduced by geographical distance (isolation-by-distance) or by divergent selection stemming from resource use (isolation-by-ecology). We tested for and encountered phenotypic and genetic divergence among Spanish crossbills utilizing different species of co-occurring pine trees as their food resource. Morphological, vocal and mtDNA divergence were not correlated with geographical distance, but they were correlated with differences in resource use. Resource diversity has now been found to repeatedly predict crossbill diversity. However, when resource use is not 100% differentiated, additional characters (morphological, vocal, genetic) must be used to uncover and validate hidden population structure. In general, this confirms that ecology drives adaptive divergence and limits neutral gene flow as the first steps towards ecological speciation, unprevented by a high potential for gene flow.     

Influence of gene flow on divergence dating – implications for the speciation history of Takydromus grass lizards

Dating the time of divergence and understanding speciation processes are central to the study of the evolutionary history of organisms but are notoriously difficult. The difficulty is largely rooted in variations in the ancestral population size or in the genealogy variation across loci. To depict the speciation processes and divergence histories of three monophyletic Takydromus species endemic to Taiwan, we sequenced 20 nuclear loci and combined with one mitochondrial locus published in GenBank. They were analysed by a multispecies coalescent approach within a Bayesian framework. Divergence dating based on the gene tree approach showed high variation among loci, and the divergence was estimated at an earlier date than when derived by the species‐tree approach. To test whether variations in the ancestral population size accounted for the majority of this variation, we conducted computer inferences using isolation‐with‐migration (IM) and approximate Bayesian computation (ABC) frameworks. The results revealed that gene flow during the early stage of speciation was strongly favoured over the isolation model, and the initiation of the speciation process was far earlier than the dates estimated by gene‐ and species‐based divergence dating. Due to their limited dispersal ability, it is suggested that geographical isolation may have played a major role in the divergence of these Takydromus species. Nevertheless, this study reveals a more complex situation and demonstrates that gene flow during the speciation process cannot be overlooked and may have a great impact on divergence dating. By using multilocus data and incorporating Bayesian coalescence approaches, we provide a more biologically realistic framework for delineating the divergence history of Takydromus.