Biogeographic barriers,Pleistocene refugia,and climatic gradients in the southeastern Nearctic drive diversification in cornsnakes (Pantherophis guttatus complex)

The southeastern Nearctic is a biodiversity hotspot that is also rich in cryptic species. Numerous hypotheses (e.g., vicariance, local adaptation, and Pleistocene speciation in glacial refugia) have been tested in an attempt to explain diversification and the observed pattern of extant biodiversity. However, previous phylogeographic studies have both supported and refuted these hypotheses. Therefore, while data support one or more of these diversification hypotheses, it is likely that taxa are forming within this region in species‐specific ways. Here, we generate a genomic data set for the cornsnakes (Pantherophis guttatus complex), which are widespread across this region, spanning both biogeographic barriers and climatic gradients. We use phylogeographic model selection combined with hindcast ecological niche models to determine regions of habitat stability through time. This combined approach suggests that numerous drivers of population differentiation explain the current diversity of this group of snakes. The Mississippi River caused initial speciation in this species complex, with more recent divergence events linked to adaptations to ecological heterogeneity and allopatric Pleistocene refugia. Lastly, we discuss the taxonomy of this group and suggest there may be additional cryptic species in need of formal recognition.     

Local and Regional Scale Genetic Variation in the Cape Dune Mole-Rat,Bathyergus suillus

The distribution of genetic variation is determined through the interaction of life history, morphology and habitat specificity of a species in conjunction with landscape structure. While numerous studies have investigated this interplay of factors in species inhabiting aquatic, riverine, terrestrial, arboreal and saxicolous systems, the fossorial system has remained largely unexplored. In this study we attempt to elucidate the impacts of a subterranean lifestyle coupled with a heterogeneous landscape on genetic partitioning by using a subterranean mammal species, the Cape dune mole-rat (Bathyergus suillus), as our model. Bathyergus suillus is one of a few mammal species endemic to the Cape Floristic Region (CFR) of the Western Cape of South Africa. Its distribution is fragmented by rivers and mountains; both geographic phenomena that may act as geographical barriers to gene-flow. Using two mitochondrial fragments (cytochrome b and control region) as well as nine microsatellite loci, we determined the phylogeographic structure and gene-flow patterns at two different spatial scales (local and regional). Furthermore, we investigated genetic differentiation between populations and applied Bayesian clustering and assignment approaches to our data. Nearly every population formed a genetically unique entity with significant genetic structure evident across geographic barriers such as rivers (Berg, Verlorenvlei, Breede and Gourits Rivers), mountains (Piketberg and Hottentots Holland Mountains) and with geographic distance at both spatial scales. Surprisingly, B. suillus was found to be paraphyletic with respect to its sister species, B. janetta–a result largely overlooked by previous studies on these taxa. A systematic revision of the genus Bathyergus is therefore necessary. This study provides a valuable insight into how the biology, life-history and habitat specificity of animals inhabiting a fossorial system may act in concert with the structure of the surrounding landscape to influence genetic distinctiveness and ultimately speciation.     

Rivers,not refugia,drove diversification in arboreal,sub‐Saharan African snakes

The relative roles of rivers versus refugia in shaping the high levels of species diversity in tropical rainforests have been widely debated for decades. Only recently has it become possible to take an integrative approach to test predictions derived from these hypotheses using genomic sequencing and paleo‐species distribution modeling. Herein, we tested the predictions of the classic river, refuge, and river‐refuge hypotheses on diversification in the arboreal sub‐Saharan African snake genus Toxicodryas. We used dated phylogeographic inferences, population clustering analyses, demographic model selection, and paleo‐distribution modeling to conduct a phylogenomic and historical demographic analysis of this genus. Our results revealed significant population genetic structure within both Toxicodryas species, corresponding geographically to river barriers and divergence times from the mid‐Miocene to Pliocene. Our demographic analyses supported the interpretation that rivers are indications of strong barriers to gene flow among populations since their divergence. Additionally, we found no support for a major contraction of suitable habitat during the last glacial maximum, allowing us to reject both the refuge and river‐refuge hypotheses in favor of the river‐barrier hypothesis. Based on conservative interpretations of our species delimitation analyses with the Sanger and ddRAD data sets, two new cryptic species are identified from east‐central Africa. This study highlights the complexity of diversification dynamics in the African tropics and the advantages of integrative approaches to studying speciation in tropical regions.     

Wolbachia Infections Mimic Cryptic Speciation in Two Parasitic Butterfly Species,Phengaris teleius and P. nausithous (Lepidoptera: Lycaenidae)

Deep mitochondrial divergence within species may result from cryptic speciation, from phylogeographic isolation or from endosymbiotic bacteria like Wolbachia that manipulate host reproduction. Phengaris butterflies are social parasites that spend most of their life in close relationship with ants. Previously, cryptic speciation has been hypothesised for two Phengaris species based on divergent mtDNA sequences. Since Phengaris species are highly endangered, the existence of cryptic species would have drastic consequences for conservation and management. We tested for cryptic speciation and alternative scenarios in P. teleius and P. nausithous based on a comprehensive sample across their Palaearctic ranges using COI gene sequences, nuclear microsatellites and tests for Wolbachia. In both species a deep mitochondrial split occurring 0.65–1.97 myrs ago was observed that did not correspond with microsatellite data but was concordant with Wolbachia infection. Haplotypes previously attributed to cryptic species were part of the Wolbachia-infected clades. In both species remaining phylogeographic structure was largely consistent between mitochondrial and nuclear genomes. In P. teleius several mitochondrial and nuclear groups were observed in East Asia while a single haplogroup and nuclear cluster prevailed across continental Eurasia. Neutrality tests suggested rapid demographic expansion into that area. In contrast, P. nausithous had several mitochondrial and nuclear groups in Europe, suggesting a complex phylogeographic history in the western part of the species range. We conclude that deep intraspecific divergences found in DNA barcode studies do not necessarily need to represent cryptic speciation but instead can be due to both infection by Wolbachia and phylogeographic structure.     

Examining the Role of Effective Population Size on Mitochondrial and Multilocus Divergence Time Discordance in a Songbird

Estimates of speciation times are subject to a number of potential errors. One source of bias is that effective population size (N_e) has been shown to influence substitution rates. This issue is of particular interest for phylogeographic studies because population sizes can vary dramatically among genetically structured populations across species’ ranges. In this study, we used multilocus data to examine temporal phylogeographic patterns in a widespread North American songbird, the Northern Cardinal (Cardinalis cardinalis). Species tree estimation indicated that the phylogeographic structure of C. cardinalis was comprised of four well-supported mainland lineages with large population sizes (large N_e) and two island lineages comprised of much smaller populations (small N_e). We inferred speciation times from mtDNA and multilocus data and found there was discordance between events that represented island-mainland divergences, whereas both estimates were similar for divergences among mainland lineages. We performed coalescent simulations and found that the difference in speciation times could be attributed to stochasticity for a recently diverged island lineage. However, the magnitude of the change between speciation times estimated from mtDNA and multilocus data of an older island lineage was substantially greater than predicted by coalescent simulations. For this divergence, we found the discordance in time estimates was due to a substantial increase in the mtDNA substitution rate in the small island population. These findings indicate that in phylogeographic studies the relative tempo of evolution between mtDNA and nuclear DNA can become highly discordant in small populations.     

Comparative phylogeography of trans‐Andean freshwater fishes based on genome‐wide nuclear and mitochondrial markers

The Neotropical region represents one of the greatest biodiversity hot spots on earth. Despite its unparalleled biodiversity, regional comparative phylogeographic studies are still scarce, with most focusing on model clades (e.g. birds) and typically examining a handful of loci. Here, we apply a genome‐wide comparative phylogeographic approach to test hypotheses of codiversification of freshwater fishes in the trans‐Andean region. Using target capture methods, we examined exon data for over 1,000 loci combined with complete mitochondrial genomes to study the phylogeographic history of five primary fish species (>150 individuals) collected from eight major river basins in Northwestern South America and Lower Central America. To assess their patterns of genetic structure, we inferred genealogical concordance taking into account all major aspects of phylogeography (within loci, across multiple genes, across species and among biogeographic provinces). Based on phylogeographic concordance factors, we tested four a priori biogeographic hypotheses, finding support for three of them and uncovering a novel, unexpected pattern of codiversification. The four emerging inter‐riverine patterns are as follows: (a) Tuira + Atrato, (b) Ranchería + Catatumbo, (c) Magdalena system and (d) Sinú + Atrato. These patterns are interpreted as shared responses to the complex uplifting and orogenic processes that modified or sundered watersheds, allowing codiversification and speciation over geological time. We also find evidence of cryptic speciation in one of the species examined and instances of mitochondrial introgression in others. These results help further our knowledge of the historical geographic factors shaping the outstanding biodiversity of the Neotropics.     

Phylogeography of a stream-dwelling frog (Pseudacris cadaverina) in southern California

Recent phylogeographic studies of animal taxa in California have revealed common geographic patterns of evolutionary divergence and genetic diversity that are generally attributable to landscape influences. However, there remains a paucity of knowledge on the evolution of freshwater taxa in southern California. Here, we investigate phylogeographic patterns in a stream-dwelling frog (Pseudacris cadaverina). Two hundred and twenty-one individuals were collected from 46 populations across the species’ range in southern California. Using 1100 bp of sequence data from cytochrome b and tRNA-Glu, we conducted phylogenetic analyses, analysis of molecular variance, and nested clade phylogeographic analysis to gain insight into the factors contributing to the distribution of genetic diversity in P. cadaverina. We tested for evidence of two putative phylogeographic breaks and tested hypotheses that genetic diversity in this species is partitioned into (1) major watersheds, (2) mountain ranges, and (3) coastal and desert regions. Our results suggest that the eastern Transverse Ranges are the center of origin for extant P. cadaverina lineages and that the observed genetic structure in this species was established during the Pleistocene Epoch. There is strong support for three major haplotype groups and a Transverse Range break in P. cadaverina that is concordant with breaks found in numerous other taxa. The distribution of genetic diversity in P. cadaverina is due in large part to the separation of populations into different major watersheds and mountain ranges. Gene flow appears to be generally limited among disjunct populations throughout the region and some desert populations have been isolated by historical habitat fragmentation.     

Distribution modelling and statistical phylogeography: an integrative framework for generating and testing alternative biogeographical hypotheses

Statistical phylogeographic studies contribute to our understanding of the factors that influence population divergence and speciation, and that ultimately generate biogeographical patterns. The use of coalescent modelling for analyses of genetic data provides a framework for statistically testing alternative hypotheses about the timing and pattern of divergence. However, the extent to which such approaches contribute to our understanding of biogeography depends on how well the alternative hypotheses chosen capture relevant aspects of species histories. New modelling techniques, which explicitly incorporate spatio-geographic data external to the gene trees themselves, provide a means for generating realistic phylogeographic hypotheses, even for taxa without a detailed fossil record. Here we illustrate how two such techniques – species distribution modelling and its historical extension, palaeodistribution modelling – in conjunction with coalescent simulations can be used to generate and test alternative hypotheses. In doing so, we highlight a few key studies that have creatively integrated both historical geographic and genetic data and argue for the wider incorporation of such explicit integrations in biogeographical studies.     

As old as the hills: Pliocene palaeogeographical processes influence patterns of genetic structure in the widespread,common shrub Banksia sessilis

The impact of Quaternary glaciation on the development of phylogeographic structure in plant species is well documented. In unglaciated landscapes, phylogeographic patterns tend to reflect processes relating to persistence and stochasticity, yet other factors, associated with the palaeogeographical history of the landscape, including geomorphological events, can also have a significant influence. The unglaciated landscape of south‐western Western Australia is an ideal location to observe these ancient drivers of lineage diversification, with tectonic activity associated with the Darling Fault in the late Pliocene attributed to patterns of deep phylogeographic divergence in a widespread tree from this region. Interestingly, other species within this region have not shown this pattern and this palaeogeographical boundary therefore presents an opportunity to examine age and historical distribution of plant species endemic to this region. In this study, we assess patterns of genetic diversity and structure across 28 populations of the widespread shrub Banksia sessilis using three cpDNA markers and nine nuclear microsatellite markers. Sixteen cpDNA haplotypes were identified, comprising two major chloroplast DNA lineages that are estimated to have diverged in the Pliocene, approximately 3.3 million years ago. This timing coincides with major geomorphological processes in the landscape, including the separation of the Darling Plateau from the adjacent Swan Coastal Plain, as well as eustatic changes on the Swan Coastal Plain that are likely to have resulted in the physical isolation of historical plant lineages. Chloroplast lineages were broadly aligned with populations associated with older lateritic soils of the Darling Plateau and Geraldton sandplains or the younger sandy soils associated with the Swan Coastal Plain and Southern Coastline. This structural pattern of lateritic versus non‐lateritic division was not observed in the nuclear microsatellite data that identified three genetic clades that roughly corresponded to populations in the North, South, and Central portions of the distributions.     

Phylogeography and Post-Glacial Recolonization in Wolverines (Gulo gulo) from across Their Circumpolar Distribution

Interglacial-glacial cycles of the Quaternary are widely recognized in shaping phylogeographic structure. Patterns from cold adapted species can be especially informative - in particular, uncovering additional glacial refugia, identifying likely recolonization patterns, and increasing our understanding of species’ responses to climate change. We investigated phylogenetic structure of the wolverine, a wide-ranging cold adapted carnivore, using a 318 bp of the mitochondrial DNA control region for 983 wolverines (n = 209 this study, n = 774 from GenBank) from across their full Holarctic distribution. Bayesian phylogenetic tree reconstruction and the distribution of observed pairwise haplotype differences (mismatch distribution) provided evidence of a single rapid population expansion across the wolverine’s Holarctic range. Even though molecular evidence corroborated a single refugium, significant subdivisions of population genetic structure (0.01< Φ_ST <0.99, P<0.05) were detected. Pairwise Φ_ST estimates separated Scandinavia from Russia and Mongolia, and identified five main divisions within North America - the Central Arctic, a western region, an eastern region consisting of Ontario and Quebec/Labrador, Manitoba, and California. These data are in contrast to the nearly panmictic structure observed in northwestern North America using nuclear microsatellites, but largely support the nuclear DNA separation of contemporary Manitoba and Ontario wolverines from northern populations. Historic samples (c. 1900) from the functionally extirpated eastern population of Quebec/Labrador displayed genetic similarities to contemporary Ontario wolverines. To understand these divergence patterns, four hypotheses were tested using Approximate Bayesian Computation (ABC). The most supported hypothesis was a single Beringia incursion during the last glacial maximum that established the northwestern population, followed by a west-to-east colonization during the Holocene. This pattern is suggestive of colonization occurring in accordance with glacial retreat, and supports expansion from a single refugium. These data are significant relative to current discussions on the conservation status of this species across its range.     

Testing comparative phylogeographic models of marine vicariance and dispersal using a hierarchical Bayesian approach

Background  

Marine allopatric speciation is an enigma because pelagic larval dispersal can potentially connect disjunct populations thereby preventing reproductive and morphological divergence. Here we present a new hierarchical approximate Bayesian computation model (HABC) that tests two hypotheses of marine allopatric speciation: 1.) "soft vicariance", where a speciation involves fragmentation of a large widespread ancestral species range that was previously connected by long distance gene flow; and 2.) peripatric colonization, where speciations in peripheral archipelagos emerge from sweepstakes colonizations from central source regions. The HABC approach analyzes all the phylogeographic datasets at once in order to make across taxon-pair inferences about biogeographic processes while explicitly allowing for uncertainty in the demographic differences within each taxon-pair. Our method uses comparative phylogeographic data that consists of single locus mtDNA sequences from multiple co-distributed taxa containing pairs of central and peripheral populations. We use the method on two comparative phylogeographic data sets consisting of cowrie gastropod endemics co-distributed in the Hawaiian (11 taxon-pairs) and Marquesan archipelagos (7 taxon-pairs).     

Historical biogeography of lowland species of toads (Bufo) across the Trans-Mexican Neovolcanic Belt and the Isthmus of Tehuantepec

Aim In this study, we investigate phylogeographic structure in two different species groups of lowland toads. First, we further investigate strict parapatry of the Pliocene‐vicariant Bufo valliceps/B. nebulifer species pair. Secondly, we test for similar phylogeographic structure in the distantly related toad B. marinus, a species we hypothesize will show a Pleistocene dispersal across the same area. Location The eastern extension of the Trans‐Mexican Neovolcanic Belt (TMNB) contacts the Atlantic Coast in central Veracruz, Mexico. Although it is not a massive structure at this eastern terminus, the TMNB has nonetheless effected vicariance and subsequent speciation in several groups of animals. The Isthmus of Tehuantepec unites the North American continent with Nuclear Central America and is also known to be a biogeographic barrier for many taxa. Methods We use sequence data from two mitochondrial DNA genes (c. 550 base‐pairs (bp) of 16S and c. 420 bp of cyt b) from 58 individuals of the B. valliceps/nebulifer complex, collected from 24 localities. We also present homologous sequence data from 23 individuals of B. marinus, collected from 12 localities. We conduct maximum‐parsimony, maximum‐likelihood and Bayesian analyses to investigate phylogeographic structure. We then use parsimony‐ and likelihood‐based topology tests to assess alternative phylogenetic hypotheses and use a previously calibrated molecular rate of evolution to estimate dates of divergence. Results Our results further define the parapatric contact zone across the TMNB between the Pliocene‐vicariant sister species B. valliceps and B. nebulifer. In contrast, phylogenetic structure among populations of B. marinus across the TMNB is much shallower, suggesting a more recent Pleistocene dispersal in this species. In addition, we found phylogeographic structure associated with the Isthmus of Tehuantepec in both species groups. Main conclusions The existence of a Pliocene–Pleistocene seaway across the Isthmus of Tehuantepec has been controversial. Our data depict clades on either side of the isthmus within two distinct species (B. valliceps and B. marinus), although none of the clades associated with the isthmus, for either species, are reciprocally monophyletic. In the B. valliceps/B. nebulifer complex, the TMNB separation appears to predate the isthmian break, whereas in B. marinus dispersal across the TMNB has occurred subsequent to the presence of a barrier at the Isthmus of Tehuantepec.     

Molecular Markers Reveal Limited Population Genetic Structure in a North American Corvid,Clark’s Nutcracker (Nucifraga columbiana)

The genetic impact of barriers and Pleistocene glaciations on high latitude resident species has not been widely investigated. The Clark’s nutcracker is an endemic North American corvid closely associated with Pinus-dominated forests. The nutcracker’s encompasses known barriers to dispersal for other species, and glaciated and unglaciated areas. Clark’s nutcrackers also irruptively disperse long distances in search of pine seed crops, creating the potential for gene flow among populations. Using the highly variable mitochondrial DNA control region, seven microsatellite loci, and species distribution modeling, we examined the effects of glaciations and dispersal barriers on population genetic patterns and population structure of nutcrackers. We sequenced 900 bp of mitochondrial control region for 169 individuals from 15 populations and analysed seven polymorphic microsatellite loci for 13 populations across the Clark’s nutcracker range. We used species distribution modeling and a range of phylogeographic analyses to examine evolutionary history. Clark’s nutcracker populations are not highly differentiated throughout their range, suggesting high levels of gene flow among populations, though we did find some evidence of isolation by distance and peripheral isolation. Our analyses suggested expansion from a single refugium after the last glacial maximum, but patterns of genetic diversity and paleodistribution modeling of suitable habitat were inconclusive as to the location of this refugium. Potential barriers to dispersal (e.g. mountain ranges) do not appear to restrict gene flow in Clark’s nutcracker, and postglacial expansion likely occurred quickly from a single refugium located south of the ice sheets.     

Limited phylogeographic and genetic connectivity in Acacia species of low stature in an arid landscape

Widespread plant species are expected to maintain genetic diversity and gene flow via pollen and seed dispersal. Stature is a key life history trait that affects seed and potentially pollen dispersal, with limited stature associated with limited dispersal and greater genetic differentiation. We sampled Hill’s tabletop wattle (Acacia hilliana) and curry wattle (Acacia spondylophylla), two co‐distributed, widespread, Acacia shrubs of low stature, across the arid Pilbara region of north‐western Australia. Using chloroplast sequence and nuclear microsatellite data we evaluated patterns of population genetic and phylogeographic diversity and structure, demographic signals, ratios of pollen to seed dispersal, evidence for historical refugia, and association between elevation and diversity. Results showed strong phylogeographic (chloroplast, G _ST = 0.831 and 0.898 for A. hilliana and A. spondylophylla, respectively) and contemporary (nuclear, F _ST = 0.260 and 0.349 for A. hilliana and A. spondylophylla, respectively) genetic structure in both species. This indicates limited genetic connectivity via seed and pollen dispersal associated with Acacia species of small stature compared to taller tree and shrub acacias across the Pilbara bioregion. This effect of stature on genetic structure is superimposed on moderate levels of genetic diversity that were expected based on widespread ranges (haplotype diversity h = 25 and 12; nuclear diversity He = 0.60 and 0.47 for A. hilliana and A. spondylophylla, respectively). Contemporary genetic structure was congruent at the greater landscape scale, especially in terms of strong genetic differentiation among geographically disjunct populations in less elevated areas. Measures of diversity and connectivity were associated with traits of greater geographic population proximity, population density, population size, and greater individual longevity, and some evidence for range expansion in A. hilliana. Results illustrate that low stature is associated with limited dispersal and greater patterns of genetic differentiation for congenerics in a common landscape and highlight the complex influence of taxon‐specific life history and ecological traits to seed and pollen dispersal.     

Paleodistributions and comparative molecular phylogeography of leafcutter ants (Atta spp.) provide new insight into the origins of Amazonian diversity

The evolutionary basis for high species diversity in tropical regions of the world remains unresolved. Much research has focused on the biogeography of speciation in the Amazon Basin, which harbors the greatest diversity of terrestrial life. The leading hypotheses on allopatric diversification of Amazonian taxa are the Pleistocene refugia, marine incursion, and riverine barrier hypotheses. Recent advances in the fields of phylogeography and species-distribution modeling permit a modern re-evaluation of these hypotheses. Our approach combines comparative, molecular phylogeographic analyses using mitochondrial DNA sequence data with paleodistribution modeling of species ranges at the last glacial maximum (LGM) to test these hypotheses for three co-distributed species of leafcutter ants (Atta spp.). The cumulative results of all tests reject every prediction of the riverine barrier hypothesis, but are unable to reject several predictions of the Pleistocene refugia and marine incursion hypotheses. Coalescent dating analyses suggest that population structure formed recently (Pleistocene-Pliocene), but are unable to reject the possibility that Miocene events may be responsible for structuring populations in two of the three species examined. The available data therefore suggest that either marine incursions in the Miocene or climate changes during the Pleistocene--or both--have shaped the population structure of the three species examined. Our results also reconceptualize the traditional Pleistocene refugia hypothesis, and offer a novel framework for future research into the area.     

Comparative phylogeography of the western Indian Ocean reef fauna

Assessing patterns of connectivity at the community and population levels is relevant to marine resource management and conservation. The present study reviews this issue with a focus on the western Indian Ocean (WIO) biogeographic province. This part of the Indian Ocean holds more species than expected from current models of global reef fish species richness. In this study, checklists of reef fish species were examined to determine levels of endemism in each of 10 biogeographic provinces of the Indian Ocean. Results showed that the number of endemic species was higher in the WIO than in any other region of the Indian Ocean. Endemic species from the WIO on the average had a larger body size than elsewhere in the tropical Indian Ocean. This suggests an effect of peripheral speciation, as previously documented in the Hawaiian reef fish fauna, relative to other sites in the tropical western Pacific. To explore evolutionary dynamics of species across biogeographic provinces and infer mechanisms of speciation, we present and compare the results of phylogeographic surveys based on compilations of published and unpublished mitochondrial DNA sequences for 19 Indo-Pacific reef-associated fishes (rainbow grouper Cephalopholis argus, scrawled butterflyfish Chaetodon meyeri, bluespot mullet Crenimugil sp. A, humbug damselfish Dascyllus abudafur/Dascyllus aruanus, areolate grouper Epinephelus areolatus, blacktip grouper Epinephelus fasciatus, honeycomb grouper Epinephelus merra, bluespotted cornetfish Fistularia commersonii, cleaner wrasse Labroides sp. 1, longface emperor Lethrinus sp. A, bluestripe snapper Lutjanus kasmira, unicornfishes Naso brevirosris, Naso unicornis and Naso vlamingii, blue-spotted maskray Neotrygon kuhlii, largescale mullet Planiliza macrolepis, common parrotfish Scarus psicattus, crescent grunter Terapon jarbua, whitetip reef shark Triaenodon obesus) and three coastal Indo-West Pacific invertebrates (blue seastar Linckia laevigata, spiny lobster Panulirus homarus, small giant clam Tridacna maxima). Heterogeneous and often unbalanced sampling design, paucity of data in a number of cases, and among-species discrepancy in phylogeographic structure precluded any generalization regarding phylogeographic patterns. Nevertheless, the WIO might have been a source of haplotypes in some cases and it also harboured an endemic clade in at least one case. The present survey also highlighted likely cryptic species. This may eventually affect the accuracy of the current checklists of species, which form the basis of some of the recent advances in Indo-West Pacific marine ecology and biogeography.     

Integrating coalescent and ecological niche modeling in comparative phylogeography

Understanding the factors that contribute to the formation of population genetic structure is a central goal of phylogeographic research, but achieving this goal can be complicated by the stochastic variance inherent to genetic processes. Statistical approaches to testing phylogeographic hypotheses accommodate this stochasticity by evaluating competing models of putative historical population structure, often by simulating null distributions of the expected variance. The effectiveness of these tests depends on the biological realism of the models. Information from the fossil record can aid in reconstructing the historical distributions of some taxa. However, for the majority of taxa, which lack sufficient fossils, paleodistributional modeling can provide valuable spatial-geographic data concerning ancestral distributions. Paleodistributional models are generated by projecting ecological niche models, which predict the current distribution of each species, onto a model of past climatic conditions. Here, we generate paleodistributional models describing the suitable habitat during the last glacial maximum for lineages from the mesic forests of the Pacific Northwest of North America, and use these models to generate alternative phylogeographic hypotheses. Coalescent simulations are then used to test these hypotheses to improve our understanding of the historical events that promoted the formation of population genetic structure in this ecosystem. Results from Pacific Northwest mesic forest organisms demonstrate the utility of these combined approaches. Paleodistribution models and population genetic structure are congruent across three amphibian lineages, suggesting that they have responded in a concerted manner to environmental change. Two other species, a willow and a water vole, despite being currently codistributed and having similar population genetic structure, were predicted by the paleodistributional model to have had markedly different distributions during the last glacial maximum. This suggests that congruent phylogeographic patterns can arise from incongruent ancestral distributions. Paleodistributional models introduce a much-needed spatial-geographic perspective to statistical phylogeography. In conjunction with coalescent models of population genetic structure, they have the potential to improve our understanding of the factors that promote population divergence and ultimately produce regional patterns of biodiversity.     

The high Andes, gene flow and a stable hybrid zone shape the genetic structure of a wide-ranging South American parrot

Background

While the gene flow in some organisms is strongly affected by physical barriers and geographical distance, other highly mobile species are able to overcome such constraints. In southern South America, the Andes (here up to 6,900 m) may constitute a formidable barrier to dispersal. In addition, this region was affected by cycles of intercalating arid/moist periods during the Upper/Late Pleistocene and Holocene. These factors may have been crucial in driving the phylogeographic structure of the vertebrate fauna of the region. Here we test these hypotheses in the burrowing parrot Cyanoliseus patagonus (Aves, Psittaciformes) across its wide distributional range in Chile and Argentina.

Results

Our data show a Chilean origin for this species, with a single migration event across the Andes during the Upper/Late Pleistocene, which gave rise to all extant Argentinean mitochondrial lineages. Analyses suggest a complex population structure for burrowing parrots in Argentina, which includes a hybrid zone that has remained stable for several thousand years. Within this zone, introgression by expanding haplotypes has resulted in the evolution of an intermediate phenotype. Multivariate regressions show that present day climatic variables have a strong influence on the distribution of genetic heterogeneity, accounting for almost half of the variation in the data.

Conclusions

Here we show how huge barriers like the Andes and the regional environmental conditions imposed constraints on the ability of a parrot species to colonise new habitats, affecting the way in which populations diverged and thus, genetic structure. When contact between divergent populations was re-established, a stable hybrid zone was formed, functioning as a channel for genetic exchange between populations.     

Population genetics and geometric morphometrics of the <Emphasis Type="Italic">Bombus ephippiatus</Emphasis> species complex with implications for its use as a commercial pollinator

Mexico and Central America are among the most biodiverse regions on Earth, harboring many species with high levels of interpopulation morphological and genetic diversity. The mountainous topography of this region contains isolated sky island habitats that have the potential to promote speciation. This has been studied in vertebrates, yet few studies have examined the phylogeographic and genetic structure of insect species encompassing this region. Here we investigate geographic patterns of genetic and morphological divergence and speciation among widespread populations of the highly polymorphic bumble bee Bombus ephippiatus and its closest relative B. wilmattae. We used DNA sequences from a fragment of cytochrome oxidase I (COI), genotypes for twelve microsatellite markers, and morphometric data from wings to construct a well-supported inference of the divergences among these taxa. We have found complex patterns of genetic isolation and morphological divergence within B. ephippiatus across its geographic range and present evidence that B. ephippiatus comprises multiple independent evolutionary lineages. The pattern of their diversification corresponds to geographic and environmental isolating mechanisms, including the Mexican highlands, the lowlands of the Isthmus of Tehuantepec in southern Mexico, the Nicaraguan Depression, the patchily distributed volcanic ranges in Nuclear Central America and Pleistocene glacial cycles. These results have important implications for the development and distribution of B. ephippiatus as a commercial pollinator in Mexico and Central America.     

Differentiation between populations of Japanese grenadier anchovy (Coilia nasus) in Northwestern Pacific based on ISSR markers: Implications for biogeography

Northwestern Pacific provides unique scenarios for studying the roles of geography and ecology in driving population divergence and speciation. To elucidate the phylogeographic pattern of Coilia nasus in Northwestern Pacific, we examined ninety individuals from five localities along the coastal regions of China and the Ariake Bay of Japan by using seven ISSR markers. Analyses of molecular variance (AMOVA) showed that genetic differentiation among groups is relatively high (F_CT = 0.1904; P = 0.000). Bayesian analysis of ISSR data also revealed significant population structuring between Chinese and Japanese locations. Phylogenetic reconstructions show reciprocal monophyly in populations between China and the Ariake Bay of Japan. We conclude that the present-day phylogeographic pattern is the result of genetic isolation between Japanese and Chinese populations in the Northwestern Pacific following the glacial retreat, and that life-history traits and ecology may play a pivotal role in shaping the realized geographical distribution pattern of this species.