Phylogeography of a mountain lizard species: an ancient fragmentation process mediated by riverine barriers in the Liolaemus monticola complex (Sauria: Liolaemidae)

Liolaemus monticola is a mountain lizard species, with a widespread distribution from central Chile that displays several highly polymorphic chromosomal races. Our study determined the phylogeographic structuring and relationships among three chromosomal races of L. monticola in Chile. Mitochondrial DNA (mtDNA) sequences of the cytochrome b gene were examined using the following phylogenetic methods: maximum parsimony, maximum likelihood, Bayesian inference and nested clade phylogeographic analyses (NCPAs). These methods revealed two major monophyletic clades (north and south) in the L. monticola species, with non-overlapping geographical locations separated by the Maipo and Yeso rivers (except one hybrid, from a zone of secondary contact). The NCPA showed that a past fragmentation process likely resulted in the separation of the two clades. The southern clade includes all samples of the 'Southern, 2 n  = 34' race; the northern clade is comprised of all remaining derived chromosomal races: the 'Northern, 2 n  = 38–40 and the Multiple Fission, 2 n  = 42–44' races. Our results support the hypothesis of a geographical and genetic split resulting from allopatric processes caused by riparian barriers acting over a long time period. The inferred biogeographical scenario shows that populations have moved from the south to the north using the Andean mountains as the primary corridor for dispersal.     

Phylogeny and phylogeography of the Liolaemus darwinii complex (Squamata: Liolaemidae): evidence for introgression and incomplete lineage sorting

Although mitochondrial DNA markers have several properties that make them suitable for phylogeographic studies, they are not free of difficulties. Phylogeographic inferences within and between closely related species can be mislead by introgression and retention of ancestral polymorphism. Here we combine different phylogenetic, phylogeographic, and population genetic methods to extract the maximum information from the Liolaemus darwinii complex. We estimate the phylogeographic structure of L. darwinii across most of its distributional range, and we then estimate relationships between L. darwinii and the syntopic species L. laurenti and L. grosseorum. Our results suggest that range expansion of these lineages brought them into secondary contact in areas where they are presently in syntopy. Here we present the first evidence for introgression in lizards from temperate South America (of L. danwinii mitochondrial DNA into L. laurenti and L. grosseorum), and for incomplete lineage sorting (between L. darwinii and L. laurenti). We show that a combination of methods can provide additional support for inferences derived from any single method and thus provide more robust interpretations and narrow the range of plausible hypotheses about mechanisms and processes of divergence. Additional studies are needed in this group of lizards and in other codistributed groups to determine if Pleistocene climatic changes could be a general factor influencing the evolutionary history of a regional biota.     

Phylogeography of Caribbean lizard malaria: tracing the history of vector‐borne parasites

The Anolis lizards of the eastern Caribbean islands are parasitized by several species of malaria parasites (Plasmodium). Here I focus on two species of Plasmodium, using molecular data (mitochondrial cytochrome b sequences) to recover the phylogeography of the parasites throughout the Lesser Antilles and Puerto Rico. The two parasites were originally described as a single species, P. azurophilum, which infects both red and white blood cells. Here the two species are termed P. azurophilum Red and P. azurophilum White based on their host cell type. Six haplotypes were found in 100 infections sequenced of P. azurophilum Red and six in 45 infections of P. azurophilum White. Nested clade analysis revealed a significant association of geographical location and clades as well as a pattern of past fragmentation of parasite populations. This is consistent with the hypothesis that vector‐borne parasites such as malaria may be subject to frequent local extinctions and recolonizations. Comparison of the phylogeography of the lizard and parasites shows only weak concordance; that is, the parasites colonized the lizards in the islands, but dispersal events between islands via vectors or failed lizard colonizations were present. The two parasites had different histories, P. azurophilum Red colonized the islands from both the north and south, and P. azurophilum White originated in the central Lesser Antilles, probably from P. azurophilum Red, then moved to both north and south. This is the first study to examine the biogeography of a pair of sibling species of vector‐borne parasites within an island archipelago system.     

Multilocus phylogeography of the Patagonian lizard complex Liolaemus kriegi (Iguania: Liolaemini)

This study presents a detailed phylogeographical analysis of one of the most conspicuous groups of lizards in northwestern Patagonia, the Liolaemus kriegi complex. This region is geographically very complex as a result of Andean orogeny and subsequent volcanism coupled with a long history of glaciations and climatic changes. For 247 individuals we sequenced one mitochondrial gene (cytochrome b) and for a subset we sequenced another mitochondrial gene [12S ribosomal RNA (12S)] and two nuclear fragments [kinesin family member 24 (KIF24) and BA3 ribosomal RNA (BA3)]. We obtained gene trees and mitochondrial and nuclear haploytpe networks, and estimated genetic distances between the main lineages and basic molecular diversity indices. We also performed spatial analysis of molecular variance (SAMOVA) and Bayesian Skyline Plot (BSP) analyses, and concordant patterns from different lines of evidence permitted delimitation of seven lineages: two described species, Liolaemus buergeri and Liolaemus tregenzai; four candidate species, Liolaemus sp. A, Liolaemus sp. B, Liolaemus sp. C, and Liolaemus sp. D; and one lineage that includes all individuals from the geographical range of Liolaemus ceii and L. kriegi, referred to as L. kriegi + L. ceii. We discuss the evolutionary processes that may contribute to the origin of these lineages and their taxonomic and conservation implications. © 2014 The Linnean Society of London, Biological Journal of the Linnean Society, 2014, 113 , 256–269.     

Phylogeography of the mole-shrew (Anourosorex yamashinai) in Taiwan: implications of interglacial refugia in a high-elevation small mammal

To test the Pleistocene interglacial refugia hypothesis with a high-elevation mammal, we studied the phylogeography of the mole-shrew (Anourosorex yamashinai) using partial mitochondrial cytochrome b gene sequences (737 bases). This shrew is endemic to Taiwan. It is mainly distributed in the highlands from 1000 to 2500 m in elevation. We examined 103 specimens from 24 localities in three mountain ranges of Taiwan and found 36 haplotypes. These haplotypes separated into two major phylogroups (Northern and Southern) plus a minor phylogroup (Houhuan) of only one haplotype. This demonstrated strong association with geography. The formation of these three phylogroups may be the result of interglacial refugia during the middle Pleistocene. Distinct sublineages were not found within each major phylogroup, suggesting that the populations (phylogroups) explosively expanded from the interglacial refugia of ancestral founder haplotypes. The present distribution pattern of haplotypes suggests that Mount Houhuan is an effective refugium in central Taiwan. It was not possible to specify the refugia for the Northern and Southern phylogroups.     

Phylogeography of the tree lizard, Urosaurus ornatus: responses of populations to past climate change

Isolation due to both geological barriers and range contractions during the Pleistocene glacial maxima has been an important cause of diversification of arid-adapted species in the North American deserts. Tree lizards, Urosaurus ornatus, are distributed across much of the southwestern arid regions and can tolerate a wide range of environments. Thus, they may have avoided large-scale shifts in distribution caused by Pleistocene climate change and any subsequent evolutionary impacts. Cytochrome b sequences were sampled from U. ornatus across the northern part of their range to test if current structure of these populations resulted from post-Pleistocene range expansion and habitat fragmentation, or prior geological isolation. Phylogenetic analyses found geographical structuring of populations consistent with a model of long-term geographical isolation corresponding to each of the desert regions. The two post-Pleistocene hypotheses were not well supported as estimated times of divergence predated the retreat of the last continental ice sheet. Populations in different regions were impacted by different processes. Southern populations of U. ornatus appear to have remained largely independent of more derived northern and eastern populations during Pleistocene climate change, while populations in regions containing more derived populations showed evidence of more recent range expansion (Colorado Plateau). As populations of U. ornatus attest to, the complex and dynamic history of the southwestern USA has left a deep-rooted and multifaceted imprint on genetic and phylogeographical structure of the species living there.     

Phylogeography of lethal male fighting in a social spider mite

When males fight for access to females, such conflict rarely escalates into lethal fight because the risks and costs involved, that is, severe injury or death, are too high. The social spider mite, Stigmaeopsis miscanthi, does exhibit lethal male fights, and this male–male aggressiveness varies among populations. To understand the evolution of lethal fighting, we investigated aggressiveness in 42 populations and phylogenetic relationships in 47 populations along the Japanese archipelago. By analysis of the male weapon morph, a proxy for aggressiveness, we confirmed the existence of a mildly aggressive (ML) form, besides the low aggression (LW) and high aggression (HG) forms reported earlier. To evaluate demographic history of these three forms, we employed the approximate Bayesian computation approach using mtCOI sequences and taking into consideration the postlast glacial expansion history of the host plant, Miscanthus sinensis. As results, hierarchical split models are more likely to explain the observed genetic pattern than admixture models, and the ML form in the subtropical region was considered the ancestral group. The inferred demographic history was consistent with the one reconstructed for the host plant in a previous study. The LW form was split from the ML form during the last glacial period (20,000–40,000 years BP), and subsequently, the HG form was split from the ML form at the end of or after the last glacial period (5,494–10,988 years BP). The results also suggest that the mite invaded Japan more than once, resulting in the present parapatric distribution of LW and HG forms in eastern Japan.     

Phylogeography of the red-tailed chipmunk (Tamias ruficaudus), a northern Rocky Mountain endemic

The northern Rocky Mountains have experienced a complex history of geological events and environmental fluctuation, including Pleistocene glaciation. To provide an initial assessment of the genetic impact of this history on the regional biota we estimated phylogenetic relationships within Tamias ruficaudus, a regional endemic, from cytochrome b sequence variation using parsimony, maximum likelihood, and nested clade analysis. Analyses of sequence variation in 187 individuals from 43 localities across the distribution of T. ruficaudus indicate a history of vicariance events and range fluctuation consistent with successive periods of extensive Pleistocene glaciation in the northern Rocky Mountains. Intraspecific divergence levels (c. 4.7% uncorrected) and phylogenetic structure are consistent with a genealogical vicariance initiated prior to the Late Pleistocene, whereas nested clade analyses indicate more recent population history structured by both fragmentation and range expansion. A comparison of sequence variation with bacular morphology indicates that the two genetically and morphologically differentiated entities exhibit a zone of differential character introgression. Sequence data support a multiple refugia hypothesis and provide a phylogeographical case study for the ongoing synthesis of regional biogeography for northern Rocky Mountain endemics.     

Phylogeography of the dusky shrew, Sorex monticolus (Insectivora, Soricidae): insight into deep and shallow history in northwestern North America

Phylogenetic relationships among the dusky shrew (Sorexmonticolus) and eight related species (S. bairdi, S. bendirii, S. neomexicanus, S.ornatus, S. pacificus, S. palustris, S. sonomae and S.vagrans) were assessed using sequences from the mitochondrial cytochrome b gene (801 bp). Analyses using parsimony and maximum likelihood revealed significant molecular variation not reflected in previous morphological studies of these species. Conversely, three morphologically defined species (S.bairdi, S.neomexicanus and S.pacificus) were poorly differentiated. Sorexornatus and S.vagrans represented basal taxa for a more inclusive group that included: (i) a widespread Continental clade containing S.monticolus (Arizona to Alaska, including S. neomexicanus); (ii) a Coastal clade containing S.monticolus (Oregon to south-east Alaska, including S. bairdi and S. pacificus); (iii) the semiaquatic species (S. bendirii and S. palustris); and (iv) S.sonomae. Additional subdivision was observed within the Continental clade corresponding to populations from the northern and southern Rocky Mountains. Average uncorrected sequence divergence between the Coastal and Continental clades was 5.3% (range 4.5-6.2%), which exceeds many interspecific comparisons within this species complex and within the genus Sorex. Lack of resolution of internal nodes within topologies suggests a deep history of rapid diversification within this group. Late Pleistocene/Holocene glacial perturbations are reflected in the shallow phylogeographic structure within these clades in western North America. Our results suggest also that S. monticolus is not monophyletic under current taxonomic nomenclature. This perspective on phylogeographic history was developed within a growing comparative framework for other organisms in western North America.     

Genetic and morphological variation in a Mediterranean glacial refugium: evidence from Italian pygmy shrews,Sorex minutus (Mammalia: Soricomorpha)

At the Last Glacial Maximum (LGM), the southern European peninsulas were important refugia for temperate species. Current genetic subdivision of species within these peninsulas may reflect past population subdivision at the LGM, as in ‘refugia within refugia’, and/or at other time periods. In the present study, we assess whether pygmy shrew populations from different regions within Italy are genetically and morphologically distinct. One maternally and two paternally inherited molecular markers (cytochrome b and Y‐chromosome introns, respectively) were analysed using several phylogenetic methods. A geometric morphometric analysis was performed on mandibles to evaluate size and shape variability between populations. Mandible shape was also explored with a functional approach that considered the mandible as a first‐order lever affecting bite force. We found genetically and morphologically distinct European, Italian, and southern Italian groups. Mandible size increased with decreasing latitude and southern Italian pygmy shrews exhibited mandibles with the strongest bite force. It is not clear whether or not the southern Italian and Italian groups of pygmy shrews occupied different refugia within the Italian peninsula at the LGM. It is likely, however, that geographic isolation earlier than the LGM on islands at the site of present‐day Calabria was important in generating the distinctive southern Italian group of pygmy shrews, and also the genetic groups in other small vertebrates that we review here. Calabria is an important hotspot for genetic diversity, and is worthy of conservation attention. © 2010 The Linnean Society of London, Biological Journal of the Linnean Society, 2010, 100 , 774–787.     

Phylogeography of Southern Water Vole (Arvicola sapidus): evidence for refugia within the Iberian glacial refugium?

The role of Southern European peninsulas as glacial refugia for temperate species has been widely established, but phylogeographic patterns within refugia are being only recently addressed. Here we describe the phylogeographic patterns for Southern water vole ( Arvicola sapidus ) in its whole distribution across Iberia and France. Control region and cytochrome b sequences were obtained for 228 samples from 130 localities across Iberia and France. Eighty-five haplotypes were found in total yielding a high overall mitochondrial diversity (π = 0.027; H  = 0.974). Phylogeographic structure was relatively shallow (3.1% average intraspecific divergence) with few supported clades and 95% and 90% maximum parsimony unconnected networks, but significant, as reflected in increased pairwise nucleotide divergences with distance ( r  = 0.197, P  = 0.03) and significant autocorrelation up to ∼500 km. Spatial analysis of molecular variance analysis detected seven geographical groups explaining 43.73% of the total mitochondrial variation. We detected demographic expansions in three of these groups. A recent colonization of France from Iberia was suggested and estimated around 62 000 years bp by an isolation-with-migration model. Our results suggest the contribution of episodes of isolation in glacial subrefugia in Iberia, but seem to exclude a long-term isolation over successive glacial cycles. Phylogeographic divergence was probably tempered by relatively large population sizes and rapid and extensive mixing among subrefugia during interglacials, that might have eroded the phylogeographic structure accumulated at glacial peaks. Phenotypic differences in A. sapidus do not delineate historically isolated intraspecific divisions and do not warrant subspecific delimitations. Our results do support the existence of subrefugia within Iberia and their role in promoting intraspecific divergences.     

Phylogeography of Peromyscus furvus (Rodentia; muridae) based on cytochrome b sequence data

We conducted phylogenetic analyses of cytochrome b sequence data to assess genetic variation within and among the three allopatric segments of the distribution of Peromyscus furvus from the Sierra Madre Oriental in eastern Mexico. We identified 24 unique haplotypes among the 54 individuals examined and genetic distances ranged up to 0.078 substitutions per site. Populations from the central portion of the range formed a monophyletic unit, whereas samples from the southern distributional unit were polyphyletic. Furthermore, the southernmost population sampled may represent a distinct species. This high degree of genetic differentiation among populations, currently recognized as conspecific, mirrors the result of other genetic studies of highland rodents in Mesoamerica. Together these studies indicate that the region, already considered hyperdiverse on the basis of species diversity and endemism, may contain considerably greater diversity than is currently appreciated.     

Phylogeography of the banded killifish (Fundulus diaphanus): glacial races and secondary contact

Spatial patterns of morphological and genetic polymorphisms between the two banded killifish sub-species ( Fundulus diaphanus diaphanus and F. d. menona ) provide evidence that these sub-species represent distinct evolutionary lineages corresponding to glacial races of Atlantic and Mississippian origins. Individuals with intermediate phenotypes in Lake Ontario and the upper St Lawrence River probably result from secondary contacts between these races. The hybrid zone unexpectedly extends to Lake Erie and involves the introgression of the F. d. diaphanus mitochondrial genome into fish with morphology like that of F. d. menona .     

Phylogeography of the Alpine salamander, Salamandra atra (Salamandridae) and the influence of the Pleistocene climatic oscillations on population divergence

Fifty individuals of the endemic Alpine salamander, Salamandra atra, representing 13 populations throughout the range of the two currently recognized subspecies, atra and aurorae, were examined for sequence variation in a large portion (1050 bp) of the mitochondrial cytochrome b gene. We revealed a large number of mitochondrial DNA (mtDNA) haplotypes (10). Interpopulation sequence divergence was very low, ranging from 0 to 3.1%. The relationships among haplotypes were poorly resolved. The divergence time estimate between several mtDNA haplotypes suggested a pre-Pleistocene differentiation approximately 3 million years ago. Moreover, the impact of the Pleistocene glaciations on the phylogeographical patterns appears to have been secondary, although a somewhat reduced genetic variability was found in populations living in areas that were directly affected by the glaciation.     

Phylogeography and population structure of thornback rays (Raja clavata L., Rajidae)

The phylogeography of thornback rays (Raja clavata) was assessed from European waters, using five nuclear microsatellite loci and mitochondrial cytochome b sequences. Strong regional differentiation was found between the Mediterranean basin, the Azores and the European continental shelf. Allelic and haplotype diversities were high in Portuguese populations, consistent with the existence of a refugium along the Iberian Peninsula. Unexpectedly, high diversity was also found in the English Channel/North Sea area. The lowest genetic diversity was found in the Black Sea. Populations sampled from the Mediterranean, Adriatic and Black Seas were characterized by a single mitochondrial haplotype. This haplotype was also the most ancestral and widespread outside of the Mediterranean basin except for the Azores. Populations from the Azores were dominated by a second ancestral haplotype which was shared with British populations. Results from multidimensional scaling, amova and nested clade analysis indicate that British waters are a secondary contact zone recolonized from at least two refugia--one around the Iberian Peninsula and one possibly in the Azores. Links to a potential refugium known as the Hurd Deep, between Cornwall and Brittany, are discussed. Finally, a historical demographic analysis indicates that thornback ray populations started to expand between 580,000 and 362,000 years ago, which suggests that the Last Glacial Maximum (20,000 years ago) had mainly affected the distribution of populations rather than population size.     

Phylogeography of the common vole (Microtus arvalis) with particular emphasis on the colonization of the Orkney archipelago

To investigate the human introduction of the common vole Microtus arvalis onto the Orkney islands, the complete cytochrome b gene was sequenced in 41 specimens from both Orkney (four localities) and elsewhere in their range (26 localities). Orkney voles belonged to the same phylogenetic lineage, 'Western', as individuals from France and Spain indicating southwestern Europe as the most likely source area for the islands. This result is of interest with respect to the movement and trading links of the Neolithic people who likely transported the voles. As well as the Western lineage, our phylogenetic trees revealed three other purely European lineages: the 'Italian' (single specimen from N. Italy), the 'Central' (Germany, Netherlands, Denmark) and the 'Eastern' (Hungary, Slovakia, Poland, Ukraine, Finland, European Russia). Individuals from European Russia, W. Siberia, Georgia, Ukraine and Armenia formed a fifth distinct lineage coinciding with the distribution of the 'obscurus' chromosomal form of M. arvalis. These phylogeographical data suggest that M. arvalis occupied multiple refugia during the last glaciation.     

Phylogeographic history of Patagonian lizards of the Liolaemus elongatus complex (Iguania: Liolaemini) based on mitochondrial and nuclear DNA sequences

In this study, we present a phylogeographic analysis of a group of lizards distributed in north‐western Patagonia, the Liolaemus elongatus complex. We sequenced 581 individuals for one mitochondrial gene (cytochrome‐b), and for a subset, we sequenced another mitochondrial gene (12S rRNA) and two nuclear genes: kinesin family member 24 (KIF24) and the anonymous nuclear locus LDAB1D. We estimated gene trees, mitochondrial and nuclear haploytpe networks, standard molecular diversity indices, genetic distances between lineages and Bayesian skyline plots. Our results provide evidence for recognition of seven species previously described within the L. elongatus complex: Liolaemus antumalguen, Liolaemus chillanensis, Liolaemus carlosgarini, Liolaemus burmeisteri, Liolaemus smaug, Liolaemus elongatus and Liolaemus crandalli, but we did not find sufficient evidence to support Liolaemus choique, Liolaemus shitan or Liolaemus sp. 6 as distinct species. We identified four candidate species (Liolaemus sp. 1, Liolaemus sp. 2, Liolaemus sp. 3 and Liolaemus sp. 7), and we discuss evolutionary processes that may have contributed to the origin of these lineages and their taxonomic and conservation implications.     

Phylogeography of a widespread lizard complex reflects patterns of both geographic and ecological isolation

A primary challenge for modern phylogeography is understanding how ecology and geography, both contemporary and historical, shape the spatial distribution and evolutionary histories of species. Phylogeographic patterns are the result of many factors, including geology, climate, habitat, colonization history and lineage‐specific constraints. Assessing the relative influences of these factors is difficult because few species, regions and environments are sampled in enough detail to compare competing hypotheses rigorously and because a particular phylogeographic pattern can potentially result from different evolutionary scenarios. The silky anoles (Anolis sericeus complex) of Central America and Mexico are abundant and found in all types of lowland terrestrial habitat, offering an excellent opportunity to test the relative influences of the factors affecting diversification. Here, we performed a range‐wide statistical phylogeographic analysis on restriction site‐associated DNA (RAD) markers from silky anoles and compared the phylogeographic patterns we recovered to historical and contemporary environmental and topographic data. We constructed niche models to compare niche overlap between sister lineages and conducted coalescent simulations to characterize how the major lineages of silky anoles have diverged. Our results revealed that the mode of divergence for major lineage diversification events was geographic isolation, resulting in ecological divergence between lineages, followed by secondary contact. Moreover, comparisons of parapatric sister lineages suggest that ecological niche divergence contributed to isolation by environment in this system, reflecting the natural history differences among populations in divergent environments.     

Phylogeography of field voles (Microtus agrestis) in Eurasia inferred from mitochondrial DNA sequences

In a distribution-wide phylogeographic survey of the field vole (Microtus agrestis), 75 specimens from 56 localities across Eurasia were examined for DNA sequence variation along the whole 1140 base pair (bp) mitochondrial (mt) cytochrome b gene. The species is subdivided into three main mtDNA phylogeographic groups - western, eastern and southern - with largely allopatric distributions. The western phylogeographical group is found in west and central Europe and spread most probably from a glacial refugium in the Carpathians. The eastern group covers a large range from Lithuania to central Asia, and probably originated from a southeast European source (e.g. the southern Urals or the Caucasus). The southern group occupies an area from Portugal to Hungary, with division into two distinct mtDNA sublineages that presumably derive from separate glacial refugia in the Iberian Peninsula. Molecular clock estimates suggest that the western and eastern field vole populations separated during the last glaciation, whereas the southern population dates back 0.5-0.9 Myr. High levels of mtDNA variation indicate relatively large population sizes and subdivisions within phylogeographic groups during the last glaciation. We report a possible new suture zone in east Europe.     

Phylogeography of a post-glacial colonizer: Microtus longicaudus (Rodentia: muridae)

The molecular phylogeography of Microtus longicaudus was investigated with DNA sequences of the mitochondrial cytochrome b gene. We used phylogenetic and pairwise distance methods to reconstruct the history of the species with particular emphasis on the Pacific Northwest. Genetic variation across the species was consistent with vicariant events during the Pleistocene and subsequent northern postglacial expansion following the receding Laurentide and Cordilleran ice sheets. The largest break (> 6% uncorrected sequence divergence) was found to exist between populations found southeast of the Colorado River (eastern Arizona, Colorado, Wyoming and New Mexico) and all other western populations. Other well-supported subclades were composed of samples from: (i) the islands and north coast of southeast Alaska; (ii) eastern Alaska, British Columbia, Washington and Oregon; and (iii) northern California, Idaho and Montana. Within subclades, divergence was low. Our results suggest that the close relationships among haplotypes within northern subclades are a result of recent colonization, whereas higher among-subclade divergence is caused by genetic differentiation during prolonged periods of isolation, possibly as a result of mid-Pleistocene climatic events.