Phylogeographic patterns of a lower intertidal isopod in the Gulf of California and the Caribbean and comparison with other intertidal isopods

A growing body of knowledge on the diversity and evolution of intertidal isopods across different regions worldwide has enhanced our understanding on biological diversification at the poorly studied, yet vast, sea–land interface. High genetic divergences among numerous allopatric lineages have been identified within presumed single broadly distributed species. Excirolana mayana is an intertidal isopod that is commonly found in sandy beaches throughout the Gulf of California. Its distribution in the Pacific extends from this basin to Colombia and in the Atlantic from Florida to Venezuela. Despite its broad distribution and ecological importance, its evolutionary history has been largely neglected. Herein, we examined phylogeographic patterns of E. mayana in the Gulf of California and the Caribbean, based on maximum‐likelihood and Bayesian phylogenetic analyses of DNA sequences from four mitochondrial genes (16S rDNA, 12S rDNA, cytochrome oxidase I gene, and cytochrome b gene). We compared the phylogeographic patterns of E. mayana with those of the coastal isopods Ligia and Excirolana braziliensis (Gulf of California and Caribbean) and Tylos (Gulf of California). We found highly divergent lineages in both, the Gulf of California and Caribbean, suggesting the presence of multiple species. We identified two instances of Atlantic–Pacific divergences. Some geographical structuring among the major clades found in the Caribbean is observed. Haplotypes from the Gulf of California form a monophyletic group sister to a lineage found in Venezuela. Phylogeographic patterns of E. mayana in the Gulf of California differ from those observed in Ligia and Tylos in this region. Nonetheless, several clades of E. mayana have similar distributions to clades of these two other isopod taxa. The high levels of cryptic diversity detected in E. mayana also pose challenges for the conservation of this isopod and its fragile environment, the sandy shores.     

Vicariance and Its Impact on the Molecular Ecology of a Chinese Ranid Frog Species-Complex (Odorrana schmackeri,Ranidae)

Paleogeological events and Pleistocene climatic fluctuations have had profound influences on the genetic patterns and phylogeographic structure of species in southern China. In this study, we investigated the population genetic structure and Phylogeography of the Odorrana schmackeri species complex, mountain stream-dwelling odorous frogs, endemic to southern China. We obtained mitochondrial sequences (1,151bp) of the complete ND2 gene and two flanking tRNAs of 511 individuals from 25 sites for phylogeographic analyses. Phylogenetic reconstruction revealed seven divergent evolutionary lineages, with mean pairwise (K2P) sequence distances from 7.8% to 21.1%, except for a closer ND2 distance (3.4%). The complex geological history of southern China drove matrilineal divergence in the O. schmackeri species complex into highly structured geographical units. The first divergence between lineage A+B and other lineages (C-G) had likely been influenced by the uplift of coastal mountains of Southeast China during the Mio-Pliocene period. The subsequent divergences between the lineages C-G may have followed the formation of the Three Gorges and the intensification of the East Asian summer monsoon during the late Pliocene and early Pleistocene. Demographic analyses indicated that major lineages A and C have been experienced recent population expansion (c. 0.045–0.245 Ma) from multiple refugia prior to the Last Glacial Maximum (LGM). Molecular analysis suggest that these seven lineages may represent seven different species, three described species and four cryptic species and should at least be separated into seven management units corresponding to these seven geographic lineages for conservation.     

Phylogeography of two New Zealand lizards: McCann's skink (Oligosoma maccanni) and the brown skink (O. zelandicum)

The New Zealand skink fauna has proven to be an ideal taxonomic group in which to examine the impact of climatic and geological processes on the evolution of the New Zealand biota since the Pliocene. Here we examine the phylogeography of McCann's skink (Oligosoma maccanni) in order to gain insight into the relative contribution of Pliocene and Pleistocene processes on patterns of genetic structure in the South Island biota, and investigate the phylogeography of the brown skink (O. zelandicum) to examine whether Cook Strait landbridges facilitated geneflow between the North and South Islands in the late-Pleistocene. We obtained mitochondrial DNA sequence data (ND2 and ND4; 1282bp) from across the range of both species. We examined the phylogeographic patterns evident in each species using Neighbour-Joining, Maximum Likelihood and Bayesian methods. We found substantial phylogeographic structure within O. maccanni, with seven distinct clades identified. Divergences among clades are estimated to have occurred during the Pliocene. Populations in the Otago/Southland region (south of the Waitaki River valley) formed a well-supported lineage within O. maccanni. A substantial genetic break was evident between populations in east and west Otago, either side of the Nevis-Cardrona fault system, while north-south genetic breaks were evident within the Canterbury region. Within-clade divergences in O. maccanni appear to have occurred during the mid- to late-Pleistocene. Shimodaira-Hasegawa topology tests indicated that the 'Garston' skink is not genetically distinct from O. maccanni. There was only relatively minor phylogeographic structure within O. zelandicum, with divergences among populations occurring during the mid- to late-Pleistocene. Our genetic data supports a single colonisation of the North Island by O. zelandicum from the South Island, with the estimated timing of this event (0.46mya) consistent with the initial formation of Cook Strait.     

Lineage divergence and speciation in the Web‐toed Salamanders (Plethodontidae: Hydromantes) of the Sierra Nevada,California

Peripatric speciation and the importance of founder effects have long been controversial, and multilocus sequence data and coalescent methods now allow hypotheses of peripatric speciation to be tested in a rigorous manner. Using a multilocus phylogeographical data set for two species of salamanders (genus Hydromantes) from the Sierra Nevada of California, hypotheses of recent divergence by peripatric speciation and older, allopatric divergence were tested. Phylogeographical analysis revealed two divergent lineages within Hydromantes platycephalus, which were estimated to have diverged in the Pliocene. By contrast, a low‐elevation species, Hydromantes brunus, diverged from within the northern lineage of H. platycephalus much more recently (mid‐Pleistocene), during a time of major climatic change in the Sierra Nevada. Multilocus species tree estimation and coalescent estimates of divergence time, migration rate, and growth rate reject a scenario of ancient speciation of H. brunus with subsequent gene flow and introgression from H. platycephalus, instead supporting a more recent divergence with population expansion. Although the small, peripheral distribution of H. brunus suggests the possibility of peripatric speciation, the estimated founding population size of the species was too large to have allowed founder effects to be important in its divergence. These results provide evidence for both recent speciation, most likely tied to the climatic changes of the Pleistocene, and older lineage divergence, possibly due to geological events, and add to evidence that Pleistocene glacial cycles were an important driver of diversification in the Sierra Nevada.     

Explosive radiation and spatial expansion across the cold environments of the Old World in an avian family

Our objective was to elucidate the biogeography and speciation patterns in an entire avian family, which shows a complex pattern of overlapping and nonoverlapping geographical distributions, and much variation in plumage, but less in size and structure. We estimated the phylogeny and divergence times for all of the world's species of Prunella based on multiple genetic loci, and analyzed morphometric divergence and biogeographical history. The common ancestor of Prunella was present in the Sino‐Himalayan Mountains or these mountains and Central Asia–Mongolia more than 9 million years ago (mya), but a burst of speciations took place during the mid‐Pliocene to early Pleistocene. The relationships among the six primary lineages resulting from that differentiation are unresolved, probably because of the rapid radiation. A general increase in sympatry with increasing time since divergence is evident. With one exception, species in clades younger than c. 3.7 my are allopatric. Species that are widely sympatric, including the most recently diverged (2.4 mya) sympatric sisters, are generally more divergent in size/structure than allo‐/parapatric close relatives. The distributional pattern and inferred ages suggest divergence in allopatry and substantial waiting time until secondary contact, likely due to competitive exclusion. All sympatrically breeding species are ecologically segregated, as suggested by differences in size/structure and habitat. Colonizations of new areas were facilitated during glacial periods, followed by fragmentation during interglacials—contrary to the usual view that glacial periods resulted mainly in fragmentations.     

Phylogeographic structure and deep lineage diversification of the red alga Chondrus ocellatus Holmes in the Northwest Pacific

A major goal of phylogeographic analysis using molecular markers is to understand the ecological and historical variables that influence genetic diversity within a species. Here, we used sequences of the mitochondrial Cox1 gene and nuclear internal transcribed spacer to reconstruct its phylogeography and demographic history of the intertidal red seaweed Chondrus ocellatus over most of its geographical range in the Northwest Pacific. We found three deeply separated lineages A, B and C, which diverged from one another in the early Pliocene–late Miocene (c. 4.5–7.7 Ma). The remarkably deep divergences, both within and between lineages, appear to have resulted from ancient isolations, accelerated by random drift and limited genetic exchange between regions. The disjunct distributions of lineages A and C along the coasts of Japan may reflect divergence during isolation in scattered refugia. The distribution of lineage B, from the South China Sea to the Korean Peninsula, appears to reflect postglacial recolonizations of coastal habitats. These three lineages do not coincide with the three documented morphological formae in C. ocellatus, suggesting that additional cryptic species may exist in this taxon. Our study illustrates the interaction of environmental variability and demographic processes in producing lineage diversification in an intertidal seaweed and highlights the importance of phylogeographic approaches for discovering cryptic marine biodiversity.     

A geographic mosaic of evolutionary lineages within the insular endemic newt Euproctus montanus

Islands are hotspots of biodiversity, with a disproportionately high fraction of endemic lineages, often of ancient origin. Nevertheless, intra‐island phylogeographies are surprisingly scarce, leading to a scanty knowledge about the microevolutionary processes induced on island populations by Plio‐Pleistocene climatic oscillations, and the manner in which these processes contributed to shape their current genetic diversity. We investigated the phylogeography, historical demography and species distribution models of the Corsican endemic newt Euproctus montanus (north‐western Mediterranean). As for many island endemics, the continuous distribution of E. montanus throughout its range has hitherto been considered as evidence for a single large population, a belief that also guided the species' categorization for conservation purposes. Instead, we found a geographic mosaic of ancient evolutionary lineages, with five main clades of likely Pliocene origin (2.6–5.8 My), all but one restricted to northern Corsica. Moreover, the copresence between main lineages in the same population was limited to a single case. As also suggested by growing literature on intra‐island phylogeographic variation, it seems that the extensive use of simplifying assumption on the population structure and historical demography of island populations—both in theoretical and applicative studies—should be carefully reconsidered, a claim that is well exemplified by the case presented here.     

Multilocus phylogeographic assessment of the California Mountain Kingsnake (Lampropeltis zonata) suggests alternative patterns of diversification for the California Floristic Province

Phylogeographic inference can determine the timing of population divergence, historical demographic processes, patterns of migration, and when extended to multiple species, the history of communities. Single‐locus analyses can mislead interpretations of the evolutionary history of taxa and comparative analyses. It is therefore important to revisit previous single‐locus phylogeographic studies, particularly those that have been used to propose general patterns for regional biotas and the processes responsible for generating inferred patterns. Here, we employ a multilocus statistical approach to re‐examine the phylogeography of Lampropeltis zonata. Using nonparametic and Bayesian species delimitation, we determined that there are two well‐supported species within L. zonata. Ecological niche modelling supports the delimitation of these taxa, suggesting that the two species inhabit distinct climatic environments. Gene flow between the two taxa is low and appears to occur unidirectionally. Further, our data suggest that gene flow was mediated by females, a rare pattern in snakes. In contrast to previous analyses, we determined that the divergence between the two lineages occurred in the late Pliocene (c. 2.07 Ma). Spatially and temporally, the divergence of these lineages is associated with the inundation of central California by the Monterey Bay. The effective population sizes of the two species appear to have been unaffected by Pleistocene glaciation. Our increased sampling of loci for L. zonata, combined with previously published multilocus analyses of other sympatric species, suggests that previous conclusions reached by comparative phylogeographic studies conducted within the California Floristic Province should be reassessed.     

Comparative phylogeography of three skink species (Oligosoma moco, O. smithi, O. suteri; Reptilia: Scincidae) in northeastern New Zealand

Sea-level fluctuations during the Pliocene and Pleistocene have shaped the landscape of the Northland region of New Zealand. We examined the comparative phylogeography of three skink species (Oligosoma moco, O. smithi, O. suteri) in northeastern New Zealand in order to investigate the impact of the historical processes that have prevailed since the Pliocene on the Northland fauna. O. moco, O. smithi and O. suteri have similar distributions across northeastern New Zealand, frequently occurring in sympatry. We obtained sequence data from across the entire range of each species, targeting the ND2 mitochondrial gene. Using Neighbor-Joining, Maximum likelihood and Bayesian methods, our analysis revealed contrasting phylogeographic patterns in each species. We found substantial phylogeographic structure within O. moco, with three distinct clades identified. Similarly, deep phylogeographic divergence was evident within O. smithi, with three distinct clades present. Clade 1 included O. smithi populations from the Three Kings Islands and the western coastline of Northland, while Clade 2 encompassed the remainder of the range. However, since Clade 3 corresponded to a described species (O. microlepis), O. smithi might represent a species complex. In both O. moco and O. smithi, divergences among clades are estimated to have occurred in the Pliocene, with divergences within clades occurring during the Pleistocene. In contrast, genetic divergence among O. suteri populations was extremely limited and indicative of more recent divergences during the Pleistocene. The lack of phylogeographic structure in O. suteri might be a consequence of its oviparous reproductive mode, which restricted its distribution to warm northern refugia during glacial maxima. Differences in the ecology and biology of each species might have produced contrasting responses to the same historical processes, and ultimately diverse phylogeographic patterns. Our study reveals an absence of consistent and concordant phylogeographic patterns in the Northland biota, even within the same taxonomic group.     

Isolation‐driven divergence: speciation in a widespread North American songbird (Aves: Certhiidae)

Lineage, or true ‘species’, trees may differ from gene trees because of stochastic processes in molecular evolution leading to gene‐tree heterogeneity. Problems with inferring species trees because of excessive incomplete lineage sorting may be exacerbated in lineages with rapid diversification or recent divergences necessitating the use of multiple loci and individuals. Many recent multilocus studies that investigate divergence times identify lineage splitting to be more recent than single‐locus studies, forcing the revision of biogeographic scenarios driving divergence. Here, we use 21 nuclear loci from regional populations to re‐evaluate hypotheses identified in an mtDNA phylogeographic study of the Brown Creeper (Certhia americana), as well as identify processes driving divergence. Nuclear phylogeographic analyses identified hierarchical genetic structure, supporting a basal split at approximately 32°N latitude, splitting northern and southern populations, with mixed patterns of genealogical concordance and discordance between data sets within the major lineages. Coalescent‐based analyses identify isolation, with little to no gene flow, as the primary driver of divergence between lineages. Recent isolation appears to have caused genetic bottlenecks in populations in the Sierra Madre Oriental and coastal mountain ranges of California, which may be targets for conservation concerns.     

Multilocus phylogeography of a widespread savanna–woodland‐adapted rodent reveals the influence of Pleistocene geomorphology and climate change in Africa's Zambezi region

Understanding historical influences of climate and physiographic barriers in shaping patterns of biodiversity remains limited for many regions of the world. For mammals of continental Africa, phylogeographic studies, particularly for West African lineages, implicate both geographic barriers and climate oscillations in shaping small mammal diversity. In contrast, studies for southern African species have revealed conflicting phylogenetic patterns for how mammalian lineages respond to both climate change and geologic events such as river formation, especially during the Pleistocene. However, these studies were often biased by limited geographic sampling or exclusively focused on large‐bodied taxa. We exploited the broad southern African distribution of a savanna–woodland‐adapted African rodent, Gerbilliscus leucogaster (bushveld gerbil) and generated mitochondrial, autosomal and sex chromosome data to quantify regional signatures of climatic and vicariant biogeographic phenomena. Results indicate the most recent common ancestor for all G. leucogaster lineages occurred during the early Pleistocene. We documented six divergent mitochondrial lineages that diverged ~0.270–0.100 mya, each of which was geographically isolated during periods characterized by alterations to the course of the Zambezi River and its tributaries as well as regional ‘megadroughts’. Results demonstrate the presence of a widespread lineage exhibiting demographic expansion ~0.065–0.035 mya, a time that coincides with savanna–woodland expansion across southern Africa. A multilocus autosomal perspective revealed the influence of the Kafue River as a current barrier to gene flow and regions of secondary contact among divergent mitochondrial lineages. Our results demonstrate the importance of both climatic fluctuations and physiographic vicariance in shaping the distribution of southern African biodiversity.     

High genetic diversity and stable Pleistocene distributional ranges in the widespread Mexican red oak Quercus castanea Née (1801) (Fagaceae)

The Mexican highlands are areas of high biological complexity where taxa of Nearctic and Neotropical origin and different population histories are found. To gain a more detailed view of the evolution of the biota in these regions, it is necessary to evaluate the effects of historical tectonic and climate events on species. Here, we analyzed the phylogeographic structure, historical demographic processes, and the contemporary period, Last Glacial Maximum (LGM) and Last Interglacial (LIG) ecological niche models of Quercus castanea, to infer the historical population dynamics of this oak distributed in the Mexican highlands. A total of 36 populations of Q. castanea were genotyped with seven chloroplast microsatellite loci in four recognized biogeographic provinces of Mexico: the Sierra Madre Occidental (western mountain range), the Central Plateau, the Trans‐Mexican Volcanic Belt (TMVB, mountain range crossing central Mexico from west to east) and the Sierra Madre del Sur (SMS, southern mountain range). We obtained standard statistics of genetic diversity and structure and tested for signals of historical demographic expansions. A total of 90 haplotypes were identified, and 29 of these haplotypes were restricted to single populations. The within‐population genetic diversity was high (mean h_S = 0.72), and among‐population genetic differentiation showed a strong phylogeographic structure (N_ST = 0.630 > G_ST = 0.266; p < .001). Signals of demographic expansion were identified in the TMVB and the SMS. The ecological niche models suggested a considerable percentage of stable distribution area for the species during the LGM and connectivity between the TMVB and the SMS. High genetic diversity, strong phylogeographic structure, and ecological niche models suggest in situ permanence of Q. castanea populations with large effective population sizes. The complex geological and climatic histories of the TMVB help to explain the origin and maintenance of a large proportion of the genetic diversity in this oak species.     

Phylogeography of northern Dolly Varden Salvelinus malma malma based on analysis of mitochondrial DNA

The northern Dolly Varden, Salvelinus malma malma, is a typical representative of arctic fauna distributed in northeastern Asia and northwestern North America. Because its spawning habitats were affected by Pleistocene glacial advances over most of its natural range, S. m. malma is among the most interesting objects of phylogeographic and microevolutionary studies. We reconstructed the genealogy of mtDNA haplotypes from 27 Alaskan and Asian populations to study the influence of glacial and geological vicariance events on the contemporary population genetic structure, phylogeographic subdivision and distribution of the northern Dolly Varden. Analysis of restriction site states in three PCR‐amplified mtDNA regions (ND1/ND2, ND5/ND6, Cytb/D‐loop; 47% of the mitochondrial genome) resolved 75 haplotypes in 436 fish. Similar patterns of subspecific variation apparent from hierarchical diversity and nested clade analyses of mtDNA haplotypes identify weak spatial differentiation and low levels of divergence. Our results suggest that (1) demographic history has been influenced by historical range expansions and recent isolation by distance, (2) present populations from Asia and North America were colonized from one main Beringian Refugium, and (3) that this taxon's ancestral population probably experienced a bottleneck in the Beringian Refugium during the late Pleistocene (Wisconsin) glacial period.     

Complex phylogeography and historical hybridization between sister taxa of freshwater sculpin (Cottus)

Species ranges that span different geographic landscapes frequently contain cryptic species‐ or population‐level structure. Identifying these possible diversification factors can often be accomplished under a comparative phylogeographic framework. However, comparisons suffer if previous studies are limited to a particular group or habitat type. In California, a complex landscape has led to several phylogeographic breaks, primarily in terrestrial species. However, two sister taxa of freshwater fish, riffle sculpin (Cottus gulosus) and Pit sculpin (Cottus pitensis), display ranges based on morphological identifications that do not coincide with these breaks. Using a comprehensive sampling and nuclear, mitochondrial and microsatellite markers, we hypothesized that proposed species ranges are erroneous based on potential hybridization/gene flow between species. Results identified a phylogeographic signature consistent with this hypothesis, with breaks at the Coast Range Mountains and Sacramento/San Joaquin River confluence. Coastal locations of C. gulosus represent a unique lineage, and ‘true’ C. gulosus were limited to the San Joaquin basin, both regions under strong anthropogenic influence and potential conservation targets. C. pitensis limits extended historically throughout the Sacramento/Pit River basin but currently are restricted to the Pit River. Interestingly, locations in the Sacramento River contained low levels of ancestral hybridization and gene flow from C. gulosus but now appear to be a distinct population. The remaining population structure was strongly correlated with Sierra Nevada presence (high) or absence (low). This study stresses the importance of testing phylogeographic breaks across multiple taxa/habitats before conservation decisions are made, but also the potential impact of different geographic landscapes on evolutionary diversification.     

Speciation in the White-breasted Nuthatch (Sitta carolinensis): a multilocus perspective

Inferring the evolutionary and ecological processes that have shaped contemporary species distributions using the geographic distribution of gene lineages is the principal goal of phylogeographic research. Researchers in the field have recognized that inferences made from a single gene, often mitochondrial, can be informative regarding the pattern of diversification but lack conclusive information regarding the evolutionary mechanisms that led to the observed patterns. Here, we use a multilocus (20 loci) data set to explore the evolutionary history of the White‐breasted Nuthatch (Sitta carolinensis). A previous single‐locus study found S. carolinensis is comprised of four reciprocally monophyletic clades geographically restricted to the pine and oak forests of: (i) eastern North America, (ii) southern Rocky Mountain and Mexican Mountain ranges, (iii) Eastern Sierra Nevada and Northern Rocky Mountains and (iv) Pacific slope of North America. The diversification of the clades was attributed to the fragmentation of North American pine and oak woodlands in the Pliocene with subsequent divergences owing to the Pleistocene glacial cycles. Principal component, clustering and species tree analyses of the multilocus data resolved the same four groups or lineages found in the single‐locus study. Coalescent analyses and hypothesis testing of nested isolation and migration models indicate that isolation and not gene flow has been the major evolutionary mechanism responsible for shaping genetic variation, and all the divergence events within S. carolinensis have occurred in response to the Pleistocene glacial cycles.     

Biogeography and host‐related factors trump parasite life history: limited congruence among the genetic structures of specific ectoparasitic lice and their rodent hosts

Parasites and hosts interact across both micro‐ and macroevolutionary scales where congruence among their phylogeographic and phylogenetic structures may be observed. Within southern Africa, the four‐striped mouse genus, Rhabdomys, is parasitized by the ectoparasitic sucking louse, Polyplax arvicanthis. Molecular data recently suggested the presence of two cryptic species within P. arvicanthis that are sympatrically distributed across the distributions of four putative Rhabdomys species. We tested the hypotheses of phylogeographic congruence and cophylogeny among the two parasite lineages and the four host taxa, utilizing mitochondrial and nuclear sequence data. Despite the documented host‐specificity of P. arvicanthis, limited phylogeographic correspondence and nonsignificant cophylogeny was observed. Instead, the parasite–host evolutionary history is characterized by limited codivergence and several duplication, sorting and host‐switching events. Despite the elevated mutational rates found for P. arvicanthis, the spatial genetic structure was not more pronounced in the parasite lineages compared with the hosts. These findings may be partly attributed to larger effective population sizes of the parasite lineages, the vagility and social behaviour of Rhabdomys, and the lack of host‐specificity observed in areas of host sympatry. Further, the patterns of genetic divergence within parasite and host lineages may also be largely attributed to historical biogeographic changes (expansion‐contraction cycles). It is thus evident that the association between P. arvicanthis and Rhabdomys has been shaped by the synergistic effects of parasite traits, host‐related factors and biogeography over evolutionary time.     

Influence of late Quaternary climate change on present patterns of genetic variation in valley oak,Quercus lobata Née

Phylogeography and ecological niche models (ENMs) suggest that late Quaternary glacial cycles have played a prominent role in shaping present population genetic structure and diversity, but have not applied quantitative methods to dissect the relative contribution of past and present climate vs. other forces. We integrate multilocus phylogeography, climate‐based ENMs and multivariate statistical approaches to infer the effects of late Quaternary climate change on contemporary genetic variation of valley oak (Quercus lobata Née). ENMs indicated that valley oak maintained a stable distribution with local migration from the last interglacial period (~120 ka) to the Last Glacial Maximum (~21 ka, LGM) to the present compared with large‐scale range shifts for an eastern North American white oak (Quercus alba L.). Coast Range and Sierra Nevada foothill populations diverged in the late Pleistocene before the LGM [104 ka (28–1622)] and have occupied somewhat distinct climate niches, according to ENMs and coalescent analyses of divergence time. In accordance with neutral expectations for stable populations, nuclear microsatellite diversity positively correlated with niche stability from the LGM to present. Most strikingly, nuclear and chloroplast microsatellite variation significantly correlated with LGM climate, even after controlling for associations with geographic location and present climate using partial redundancy analyses. Variance partitioning showed that LGM climate uniquely explains a similar proportion of genetic variance as present climate (16% vs. 11–18%), and together, past and present climate explains more than geography (19%). Climate can influence local expansion–contraction dynamics, flowering phenology and thus gene flow, and/or impose selective pressures. These results highlight the lingering effect of past climate on genetic variation in species with stable distributions.     

How lizards survived blizzards: phylogeography of the Liolaemus lineomaculatus group (Liolaemidae) reveals multiple breaks and refugia in southern Patagonia and their concordance with other codistributed taxa

Patagonia was shaped by a complex geological history, including the Miocene uplift of the Andes, followed by volcanism, marine introgressions, and extreme climatic oscillations during Pliocene–Pleistocene glaciation–deglaciation cycles. The distributional patterns and phylogenetic relationships of southern patagonian animals and plants were affected in different ways, and those imprints are reflected in the seven phylogeographic breaks and eight refugia that have been previously proposed. In this study, we estimated time‐calibrated phylogenetic/phylogeographic patterns in lizards of the Liolaemus lineomaculatus group and relate them to historical Miocene‐to‐Pleistocene events of Patagonia and the previously proposed phylogeographic patterns. Individuals from 51 localities were sequenced for the mitochondrial marker (cyt‐b) and a subsample of individuals from each mitochondrial lineage was sequenced for one nuclear (LDA12D) and one slow evolving mitochondrial gene (12S). Our analyses revealed strong phylogeographic structure among lineages and, in most cases, no signal of demographic changes through time. The lineomaculatus group is composed of three strongly supported clades (lineomaculatus, hatcheri and kolengh + silvanae), and divergence estimates suggested their origins associated with the oldest known Patagonian glaciation (7–5 Ma); subsequent diversification within the lineomaculatus clade coincided with the large Pliocene glaciations (~3.5 Ma). The lineomaculatus clade includes nine strongly genetically and geographically structured lineages, five of which are interpreted as candidate species. Our findings suggest that some Liolaemus lineages have persisted in situ, each of them in a different refugium, through several glaciation–deglaciation cycles without demographic fluctuations. We also summarize and update qualitative evidence of some shared phylogeographic breaks and refugia among plants, rodents and lizards.     

Pleistocene climatic fluctuations drive isolation and secondary contact in the red diamond rattlesnake (Crotalus ruber) in Baja California

Aim

Many studies have investigated the phylogeographic history of species on the Baja California Peninsula, and they often show one or more genetic breaks that are spatially concordant among many taxa. These phylogeographic breaks are commonly attributed to vicariance as a result of geological or climatic changes, followed by secondary contact when barriers are no longer present. We use restriction‐site associated DNA sequence data and a phylogeographic model selection approach to explicitly test the secondary contact hypothesis in the red diamond rattlesnake, Crotalus ruber.

Location

Baja California and Southern California.

Methods

We used phylogenetic and population clustering approaches to identify population structure. We then used coalescent methods to simultaneously estimate population parameters and test the fit of phylogeographic models to the data. We used ecological niche models to infer suitable habitat for C. ruber at the Last Glacial Maximum (LGM).

Results

Crotalus ruber is composed of distinct northern and southern populations with a boundary near the town of Loreto in Baja California Sur. A model of isolation followed by secondary contact provides the best fit to the data, with both divergence and contact occurring in the Pleistocene. We also identify a genomic signature of northern range expansion in the northern population, consistent with LGM niche models showing that the northern‐most portion of the range of C. ruber was not suitable habitat during the LGM.

Main conclusions

We provide the first explicitly model‐based test of the secondary contact model in Baja California and show that populations of C. ruber were isolated before coming back into contact near Loreto, a region that shows phylogeographic breaks for other taxa. Given the timing of divergence and contact, we suggest that climatic fluctuations have driven the observed phylogeographic structure observed in C. ruber and that they may have driven similar patterns in other taxa.     

Gene flow,population growth and a novel substitution rate estimate in a subtidal rock specialist,the black‐faced blenny Tripterygion delaisi (Perciformes,Blennioidei, Tripterygiidae) from the Adriatic Sea

Population histories depend on the interplay between exogeneous and endogeneous factors. In marine species, phylogeographic and demographic patterns are often shaped by sea level fluctuations, water currents and dispersal ability. Using mitochondrial control region sequences (n = 120), we infer phylogeographic structure and historic population size changes of a common littoral fish species, the black‐faced blenny Tripterygion delaisi (Perciformes, Blennioidei, Tripterygiidae) from the north‐eastern Adriatic Sea. We find that Adriatic T. delaisi are differentiated from conspecific populations in the remaining Mediterranean, but display little phylogeographic structure within the Adriatic basin. The pattern is consistent with passive dispersal of planktonic larvae along cyclonic currents within the Adriatic Sea, but limited active dispersal of adults. Demographic reconstructions are consistent with recent population expansion, probably triggered by rising sea levels after the last glacial maximum (LGM). Placing the onset of population growth between the LGM and the warming of surface waters (18 000–13 000 years BP) and employing a novel expansion dating approach, we inferred a substitution rate of 2.61–3.61% per site per MY. Our study is one of only few existing investigations of the genetic structure of animals within the Adriatic basin and is the first to provide an estimate for mitochondrial control region substitution rates in blennioid fishes.