Mitochondrial phylogeography of the moor frog, Rana arvalis

The moor frog Rana arvalis is a lowland species with a broad Eurasiatic distribution, from arctic tundra through forest to the steppe zone. Its present-day range suggests that glacial refugia of this frog were located outside southern European peninsulas. We studied the species-wide phylogeographical pattern using sequence variation in a 682 base pairs fragment of mtDNA cytochrome b gene; 223 individuals from 73 localities were analysed. Two main clades, A and B, differing by c. 3.6% sequence divergence were detected. The A clade is further subdivided into two subclades, AI and AII differing by 1.0%. All three lineages are present in the Carpathian Basin (CB), whereas the rest of the species range, including huge expanses of Eurasian lowlands, are inhabited solely by the AI lineage. We infer that AII and B lineages survived several glacial cycles in the CB but did not expand, at least in the present interglacial, to the north. The geographical distribution and genealogical relationships between haplotypes from the AI lineage indicate that this group had two glacial refugia, one located in the eastern part of the CB and the other probably in southern Russia. Populations from both refugia contributed to the colonization of the western part of the range, whereas the eastern part was colonized from the eastern refugium only. The effective population size as evidenced by theta(ML) is an order of magnitude higher in the AI lineage than in the AII and B lineages. Demographic expansion was detected in all three lineages.     

Using nested clade analysis to assess the history of colonization and the persistence of populations of an Iberian Lizard

The distribution of the lizard Lacerta schreiberi is likely to have been severely affected by the climatic cycles that have influenced the Iberian Peninsula. Information about the species ecology and Iberian physiogeography was used to generate specific hypotheses about episodes of colonization and subsequent population persistence. These hypotheses generated predictions about the distribution of genetic variation, which were tested using nested clade analysis (NCA) supplemented by analysis of molecular variance (amova). Two predictions were confirmed by NCA; that is those that specified multiple and allopatric refugia. However, the remaining three predictions were not corroborated by the analyses. Firstly, a simple analysis of the distribution of genetic variability failed to detect an expected difference in the pattern of colonization between the inland mountain system and the coastal region. Moreover, while NCA did detect the expected genetic pattern in southern coastal populations, it was explained in terms of long-distance migration, which seems implausible because of the extent of unsuitable habitat. A more likely cause of the pattern is population fragmentation and a reduction in population size caused during the Holocene. Finally, NCA also failed to detect a northwestern population expansion, which is supported by other evidence. We conclude that NCA has a limited ability to detect range expansion led by individuals with more ancestral (interior) haplotypes.     

The Impacts of Urbanization on the Distribution and Body Condition of the Rice-paddy Frog(Fejervarya multistriata) and Gold-striped Pond Frog(Pelophylax plancyi) in Shanghai,China

Previous studies have suggested that urbanization presents a major threat to anuran populations. However, very few studies have looked at the relationship between urbanization and anuran body condition. We investigated whether the distribution and body condition of the rice-paddy frog(Fejervarya multistriata) and gold-striped pond frog(Pelophylax plancyi) are influenced by increasing urbanization in Shanghai, China. Four study sites with six indicators of the major land-cover types were scored to indicate their position on an urbanization gradient. We found that both the density and body condition of F. multistriata declined significantly along this gradient. Although we observed a significant difference in body condition of P. plancyi among study sites with different degrees of urbanization, we did not find any corresponding significant differences in population density. Our results indicate that both the densities and body condition of these two anuran species show a negative relationship with increasing urbanization, but that the density of P. plancyi was only slightly affected in Shanghai.     

Geological events play a larger role than Pleistocene climatic fluctuations in driving the genetic structure of Quasipaa boulengeri (Anura: Dicroglossidae)

Paleoclimatic and paleogeological events have been identified as being the two main drivers of genetic structuring in extant organisms. We used a montane stream‐dwelling frog, Quasipaa boulengeri, to explore the relative roles played by these drivers on species in southern China, a region needing thorough studies. We detected four major matrilines, and no broadly distributed haplotypes occurred. The complex orogenesis of south‐western China drove matrilineal divergence in Q. boulengeri into highly structured geographical units. These matrilines subsequently persisted in situ with stable populations rather than undergoing expansions during glacial cycling. The unification of the upper and middle Yangtze River in the Three Gorges mountain region mediated downstream colonization of this frog. Analyses identified geological events as playing a larger role than climatic fluctuations in driving the population history of Q. boulengeri. Nuclear allele analyses indicated gene flow; this maintained genetic cohesion of the species. South‐eastern Sichuan Basin was identified as the area of secondary contact for several matrilines, and this area deserves further study and special protection.     

Comparative phylogeography of the two pink salmon broodlines: an analysis based on a mitochondrial DNA genealogy

Over most of their natural northern Pacific Ocean range, pink salmon (Oncorhynchus gorbuscha) spawn in a habitat that was repeatedly and profoundly affected by Pleistocene glacial advances. A strictly two-year life cycle of pink salmon has resulted in two reproductively isolated broodlines, which spawn in alternating years and evolved as temporal replicates of the same species. To study the influence of historical events on phylogeographical and population genetic structure of the two broodlines, we first reconstructed a fine-scale mtDNA haplotype genealogy from a sample of 80 individuals and then determined the geographical distribution of the major genealogical assemblages for 718 individuals sampled from nine Alaskan and eastern Asian even- and nine odd-year pink salmon populations. Analysis of restriction site states in seven polymerase chain reaction (PCR)-amplified mtDNA regions (comprising 97% of the mitochondrial genome) using 13 endonucleases resolved 38 haplotypes, which clustered into five genealogical lineages that differed from 0.065 to 0.225% in net sequence divergence. The lineage sorting between broodlines was incomplete, which suggests a recent common ancestry. Within each lineage, haplotypes exhibited star-like genealogies indicating recent population growth. The depth of the haplotype genealogy is shallow ( approximately 0.5% of nucleotide sequence divergence) and probably reflects repeated decreases in population size due to Pleistocene glacial advances. Nested clade analysis (NCA) of geographical distances showed that the geographical distribution observed for mitochondrial DNA (mtDNA) haplotypes resulted from alternating influences of historical range expansions and episodes of restricted dispersal. Analyses of molecular variance showed weak geographical structuring of mtDNA variation, except for the strong subdivision between Asian and Alaskan populations within the even-year broodline. The genetic similarities observed among and within geographical regions probably originated from postglacial recolonizations from common sources rather than extensive gene flow. The phylogeographical and population genetic structures differ substantally between broodlines. This can be explained by stochastic lineage sorting in glacial refugia and perhaps different recolonization routes in even- and odd-year broodlines.     

Palaeodistribution modelling and genetic evidence highlight differential post‐glacial range shifts of a rain forest conifer distributed across a latitudinal gradient

Aim We examine the range expansion/contraction dynamics during the last glacial cycle of the late‐successional tropical rain forest conifer Podocarpus elatus using a combination of modelling and molecular marker analyses. Specifically, we test whether distributional changes predicted by environmental niche modelling are in agreement with (1) the glacial maximum contractions inferred from the southern fossil record, and (2) population genetic‐based estimates of range disjunctions and demographic dynamics. In addition, we test whether northern and southern ranges are likely to have experienced similar expansion/contraction dynamics. Location Eastern Australian tropical and subtropical rain forests. Methods Environmental niche modelling was completed for three time periods during the last glacial cycle and was interpreted in light of the known palynology. We collected 109 samples from 32 populations across the entire range of P. elatus. Six microsatellite loci and Bayesian coalescence analysis were used to infer population expansion/contraction dynamics, and five sequenced loci (one plastid and four nuclear) were used to quantify genetic structure/diversity. Results Environmental niche modelling suggested that the northern and southern ranges of P. elatus experienced different expansion/contraction dynamics. In the northern range, the habitat suitable for P. elatus persisted in a small refugial area during the Last Glacial Maximum (LGM, 21 ka) and then expanded during the post‐glacial period. Conversely, in the south suitable habitat was widespread during the LGM but subsequently contracted. These differential dynamics were supported by Bayesian analyses of the population genetic data (northern dispersal) and are consistent with the greater genetic diversity in the south compared with the north. A contact zone between the two genetically divergent groups (corresponding to the Macleay Overlap Zone) was supported by environmental niche modelling and molecular analyses. Main conclusions The climatic fluctuations of the Quaternary have differentially impacted the northern and southern ranges of a broadly distributed rain forest tree in Australia. Recurrent contraction/expansion cycles contributed to the genetic distinction between northern and southern distributions of P. elatus. By combining molecular and environmental niche modelling evidence, this unique study undermines the general assumption that broadly distributed species respond in a uniform way to climate change.     

Mitochondrial phylogeography and genetic diversity of Tibetan gazelle (Procapra picticaudata): implications for conservation

The Tibetan gazelle (Procapra picticaudata) is a threatened species and distributed on the Qinghai-Tibet Plateau of China (Qinghai Province, Tibet Autonomous Region and the adjacent Gansu Province, Sichuan Province, and Xinjiang Uigur Autonomous Region). Small peripheral populations of Tibetan gazelle were once found in northern Sikkim and Ladakh, but now these are close to extinction. To describe the evolutionary history and to assess the genetic diversity within this monotypic species and population structure among different geographic locations in China, we sequenced mitochondrial DNA from the control region (CR) and cytochrome (cyt) b gene for 46 individuals from 12 geographic localities in Qinghai, Tibet, Xinjiang, Gansu, and Sichuan. A total of 25 CR haplotypes and 16 cyt b haplotypes were identified from these gazelle samples. CR haplotype diversity (0.98+/-0.01) and nucleotide diversity (0.08+/-0.009) were both high. Phylogenetic trees indicate that the Tibetan gazelle in China can be divided into three main clades: Tibet, Sichuan (SCH) and Qinghai-Arjin Shan-Kekexili (QH-ARJ-KKXL). Analysis of molecular variance (AMOVA) and network analysis consistently support this geographic structure in both datasets. Significant differentiation between populations argues for the presence of management units (MUs). Such differentiation may reflect a geographic separation resulting from the uplift of the Qinghai-Tibet Plateau during the Late Pliocene and Pleistocene. Mismatch distribution analysis implies that Tibetan gazelle has undergone complex population changes. We suggest that the present population structure has resulted from habitat fragmentation during the recent glacial period on the Qinghai-Tibet Plateau and population expansion from glacial refugia after the glacial period. It is likely that the present populations of Tibetan gazelle exhibit a pattern reminiscent of several bottlenecks and expansions in the recent past.     

Phylogeography of Cavernularia hultenii: evidence of slow genetic drift in a widely disjunct lichen

Population structure and history is poorly known in most lichenized ascomycetes. Many species display large-scale infraspecific disjunctions, which have been explained alternately by range fragmentation in species of high age and widespread long-distance dispersal. Using the lichen Cavernularia hultenii, which is widely disjunct across North America and Europe, Pleistocene and Holocene population history was inferred. The internal transcribed spacer (ITS) and part of the the intergenic spacer (IGS) region of the nuclear ribosomal DNA were sequenced in 300 individuals representing 62 populations across the range of the species. While four ancestral haplotypes are found in all areas, none of the observed tip haplotypes is present in more than one of the three part ranges. Although this is evidence for a past fragmentation event, nested clade analysis (NCA) remains equivocal in the choice between allopatric fragmentation and long-distance dispersal. Mismatch distributions indicate exponential population growth, probably during postglacial invasion of C. hultenii into formerly glaciated areas of western North America. The presence of one southern and at least one northern glacial refugium in South Central Alaska is inferred. Evidence for another refugium in the Queen Charlotte Islands or Alexander Archipelago is inconclusive because of sparse sampling. However, a range expansion was not confirmed unambiguously by NCA. The limited power of NCA to infer past range fragmentations and expansions is due apparently to the shallow haplotype network and widespread ancestral haplotypes. This can be explained by slow genetic drift causing incomplete removal of ancestral haplotypes from the postfragmentation and postexpansion areas.     

Description, biology and conservation of a new species of Australian tree frog (Amphibia: Anura: Hylidae: Litoria) and an assessment of the remaining populations of Litoria genimaculata Horst, 1883: systematic and conservation implications of an unusual speciation event

The Australian populations of the green-eyed tree frog Litoria genimaculata consist of a northern and southern genetic lineage that meet in a mosaic contact zone comprising two independent areas of contact: one where the main ranges of the lineages overlap, and the second where a population of the southern lineage is isolated within the range of the northern lineage. A recent study failed to find significant reproductive isolation between the main ranges of the two lineages, despite deep genetic divergence, partial postzygotic isolation, and call differences. The study did, however, demonstrate rapid phenotypic divergence and speciation of the isolated population of the southern lineage from both the parapatric northern lineage and from the allopatric, but genetically similar, main range of the southern lineage. Herein, the isolated population of the southern lineage is described as a distinct species, Litoria myola sp. nov. , whereas the remainder of the southern lineage and the northern lineage are retained as a single, paraphyletic species, Litoria genimaculata . Resolving this unusual systematic situation demonstrates the value of using multiple lines of evidence in delimiting species. Litoria myola sp. nov. has a very small distribution and population size and warrants a Critically Endangered listing (B1, 2) under IUCN criteria. Threats and management recommendations are outlined, and the conservation of hybrid zones as areas of evolutionary novelty is discussed.  © 2007 The Linnean Society of London, Biological Journal of the Linnean Society , 2007, 91 , 549–563.     

Northern glacial refugia for the pygmy shrew Sorex minutus in Europe revealed by phylogeographic analyses and species distribution modelling

The southern European peninsulas (Iberian, Italian and Balkan) are traditionally recognized as glacial refugia from where many species colonized central and northern Europe after the Last Glacial Maximum (LGM). However, evidence that some species had more northerly refugia is accumulating from phylogeographic, palaeontological and palynological studies, and more recently from species distribution modelling (SDM), but further studies are needed to test the idea of northern refugia in Europe. Here, we take a rarely implemented multidisciplinary approach to assess if the pygmy shrew Sorex minutus, a widespread Eurasian mammal species, had northern refugia during the LGM, and if these influenced its postglacial geographic distribution. First, we evaluated the phylogeographic and population expansion patterns using mtDNA sequence data from 123 pygmy shrews. Then, we used SDM to predict present and past (LGM) potential distributions using two different training data sets, two different algorithms (Maxent and GARP) and climate reconstructions for the LGM with two different general circulation models. An LGM distribution in the southern peninsulas was predicted by the SDM approaches, in line with the occurrence of lineages of S. minutus in these areas. The phylogeographic analyses also indicated a widespread and strictly northern‐central European lineage, not derived from southern peninsulas, and with a postglacial population expansion signature. This was consistent with the SDM predictions of suitable LGM conditions for S. minutus occurring across central and eastern Europe, from unglaciated parts of the British Isles to much of the eastern European Plain. Hence, S. minutus likely persisted in parts of central and eastern Europe during the LGM, from where it colonized other northern areas during the late‐glacial and postglacial periods. Our results provide new insights into the glacial and postglacial colonization history of the European mammal fauna, notably supporting glacial refugia further north than traditionally recognized.     

Geography best explains global patterns of genetic diversity and postglacial co‐expansion in marine turtles

For many species, climate oscillations drove cycles of population contraction during cool glacial periods followed by expansion during interglacials. Some groups, however, show evidence of uniform and synchronous expansion, while others display differences in the timing and extent of demographic change. We compared demographic histories inferred from genetic data across marine turtle species to identify responses to postglacial warming shared across taxa and to examine drivers of past demographic change at the global scale. Using coalescent simulations and approximate Bayesian computation (ABC), we estimated demographic parameters, including the likelihood of past population expansion, from a mitochondrial data set encompassing 23 previously identified lineages from all seven marine turtle species. For lineages with a high posterior probability of expansion, we conducted a hierarchical ABC analysis to estimate the proportion of lineages expanding synchronously and the timing of synchronous expansion. We used Bayesian model averaging to identify variables associated with expansion and genetic diversity. Approximately 60% of extant marine turtle lineages showed evidence of expansion, with the rest mainly exhibiting patterns of genetic diversity most consistent with population stability. For lineages showing expansion, there was a strong signal of synchronous expansion after the Last Glacial Maximum. Expansion and genetic diversity were best explained by ocean basin and the degree of endemism for a given lineage. Geographic differences in sensitivity to climate change have implications for prioritizing conservation actions in marine turtles as well as for identifying areas of past demographic stability and potential resilience to future climate change for broadly distributed taxa.     

Range-wide phylogeography of a temperate lizard, the five-lined skink (Eumeces fasciatus)

We used mitochondrial DNA and microsatellite loci to examine the phylogeographic patterns of the most broadly distributed lizard in eastern North America, the five-lined skink (Eumeces fasciatus). We infer that longitudinal phylogeographic patterns in E. fasciatus are consistent with fragmentation due to refugial and post-glacial dynamics, but that deep divergences within the species imply historical fragmentation that predates the Pleistocene. The effect of multiple refugia is implied from our nested clade analyses, including a northern refugium in Wisconsin. Analysis of population structure using nuclear microsatellite data within the species suggests the importance of glacial dynamics in shaping more recent genetic structuring within one widely distributed lineage that ranges from the Mississippi River to the Atlantic Ocean in longitude and from southern Ontario to the Gulf of Mexico in latitude. Results shed light on the historical processes that have influenced current population structure of a temperate lizard, support the striking similarity of longitudinal phylogeographic structure across many herpetofaunal species in eastern North America, and illustrate the utility of employing multiple markers in phylogeographic studies.     

Reconstruction of the colonization route from glacial refugium to the northern distribution range of the European butterfly Polyommatus coridon (Lepidoptera: Lycaenidae)

Genetic analyses are useful tools for reconstructing glacial distribution patterns and postglacial expansion corridors. However, little information is available about the latter. We reconstruct the expansion corridors of the butterfly Polyommatus coridon from its glacial refugium to the northern edge of its current distribution by comparing populations from southern Lower Saxony (central Germany) to other existing genetic data sets. The populations from Lower Saxony clearly belonged to a western lineage that expanded postglacially from the Adriatic‐Mediterranean region. They form part of a southern German group passing through the Burgundian Gap. In the southern German group, populations belong to a western subgroup. Therefore, expansion followed the Rhine valley and through Hesse, finally reaching southern Lower Saxony and western Thuringia in central Germany. Thus, we present a complete colonization route from the glacial refugium to the northern distribution range of P. coridon. Such data are useful for understanding the biogeographic structure and migration corridors for other mobile Mediterranean species.     

Phylogeography of western Pacific Leucetta 'chagosensis' (Porifera: Calcarea) from ribosomal DNA sequences: implications for population history and conservation of the Great Barrier Reef World Heritage Area (Australia)

Leucetta 'chagosensis' is a widespread calcareous sponge, occurring in shaded habitats of Indo-Pacific coral reefs. In this study we explore relationships among 19 ribosomal DNA sequence types (the ITS1-5.8S-ITS2 region plus flanking gene sequences) found among 54 individuals from 28 locations throughout the western Pacific, with focus on the Great Barrier Reef (GBR). Maximum parsimony analysis revealed phylogeographical structuring into four major clades (although not highly supported by bootstrap analysis) corresponding to the northern/central GBR with Guam and Taiwan, the southern GBR and subtropical regions south to Brisbane, Vanuatu and Indonesia. Subsequent nested clade analysis (NCA) confirmed this structure with a probability of > 95%. After NCA of geographical distances, a pattern of range expansion from the internal Indonesian clade was inferred at the total cladogram level, as the Indonesian clade was found to be the internal and therefore oldest clade. Two distinct clades were found on the GBR, which narrowly overlap geographically in a line approximately from the Whitsunday Islands to the northern Swain Reefs. At various clade levels, NCA inferred that the northern GBR clade was influenced by past fragmentation and contiguous range expansion events, presumably during/after sea level low stands in the Pleistocene, after which the northern GBR might have been recolonized from the Queensland Plateau in the Coral Sea. The southern GBR clade is most closely related to subtropical L. 'chagosensis', and we infer that the southern GBR probably was recolonized from there after sea level low stands, based on our NCA results and supported by oceanographic data. Our results have important implications for conservation and management of the GBR, as they highlight the importance of marginal transition zones in the generation and maintenance of species rich zones, such as the Great Barrier Reef World Heritage Area.     

Post-Pleistocene range expansion of the recently imperiled eastern little brown bat (Myotis lucifugus lucifugus) from a single southern refugium

Myotis lucifugus, once among the most widespread and common bats in North America, has been forecast to be extirpated east of the Rockies in as few as 16 years by the spread of white-nose syndrome. Recent genetic research has demonstrated that this species is paraphyletic and part of a broader species complex; however, only one lineage (Myotis lucifugus lucifugus [M. l. lucifugus]) is present in eastern North America. I used molecular tools and niche modeling to validate this and investigate the role that historical biogeography has played in the phylogenetic and population genetic structure of this species to determine if the eastern subspecies represents an evolutionarily distinct population.To establish the genetic structure within M. l. lucifugus, I densely sampled maternity colonies in Minnesota and sequenced 182 individuals for a portion of cytochrome b. Phylogenetic reconstruction and a haplotype network were used to infer the relationships among mitochondrial haplotypes. Population growth statistics were calculated to determine if there was evidence of significant expansion, and an environmental niche model (ENM) was constructed based on conditions during the last glacial maximum (LGM) to illustrate potential glacial refugia. All individuals derived from a single mitochondrial lineage. Genetic evidence points to population growth starting approximately 18 kya. ENM results show that there was likely a single large southern refugium extending across the southeastern United States and possibly several isolated refugia in western North America. Myotis lucifugus lucifugus likely maintained both a large range and a large population during the peaks of the glacial cycles, and its population appears to have expanded following the retreat of the Laurentide ice sheet. This imperiled lineage likely diverged in isolation from other members of the M. lucifugus/western long-eared Myotis during the Pleistocene.     

Amphibian phylogeography in the Antipodes: Refugia and postglacial colonization explain mitochondrial haplotype distribution in the Patagonian frog Eupsophus calcaratus (Cycloramphidae)

Climatic oscillations, heterogeneity in elevation, topographical position, and isolation time in southwestern Patagonia have been important in promoting diversification of the biota. Geological studies have shown that this region had wide ice-free areas during periods of the last glacial maximum and provided forested refugia for the biota during Pleistocene glaciations. In this study, we sampled the endemic frog Eupsophus calcaratus from 20 localities, covering most of its distribution and including glaciated and non-glaciated regions. We collected DNA sequences for three mitochondrial regions (D-loop, cyt b, 16S), and describe patterns of variation consistent with a history of both the displacement to glacial refugia and recent recolonization to extensively glaciated regions. The inferred demographic history and divergence times of the lineages of E. calcaratus suggest that the Pleistocene had profound effects on the genetic patterns within this taxon in which some populations were able to survive in refugia within colder regions followed by demographic increases but without evidence of significant range expansion. The mtDNA gene tree recovers six major haploclades of E. calcaratus, which we consider diagnostic of species lineages. These results contribute to our understanding of how geological events, predominately glacial oscillations, have influenced current population structure of a broad-ranging, ectothermic vertebrate in the Valdivian Forest region of southern South America.     

Population genetic structure and approximate Bayesian computation analyses reveal the southern origin and northward dispersal of the oriental fruit moth Grapholita molesta (Lepidoptera: Tortricidae) in its native range

The oriental fruit moth (OFM) Grapholita molesta is one of the most destructive orchard pests. Assumed to be native to China, the moth is now distributed throughout the world. However, the evolutionary history of this moth in its native range remains unknown. In this study, we explored the population genetic structure, dispersal routes and demographic history of the OFM in China and South Korea based on mitochondrial genes and microsatellite loci. The Mantel test indicated a significant correlation between genetic distance and geographical distance in the populations. Bayesian analysis of population genetic structure (baps ) identified four nested clusters, while the geneland analysis inferred five genetic groups with spatial discontinuities. Based on the approximate Bayesian computation approach, we found that the OFM was originated from southern China near the Shilin area of Yunnan Province. The early divergence and dispersal of this moth was dated to the Penultimate glaciation of Pleistocene. Further dispersal from southern to northern region of China occurred before the last glacial maximum, while the expansion of population size in the derived populations in northern region of China occurred after the last glacial maximum. Our results indicated that the current distribution and structure of the OFM were complicatedly influenced by climatic and geological events and human activities of cultivation and wide dissemination of peach in ancient China. We provide an example on revealing the origin and dispersal history of an agricultural pest insect in its native range as well as the underlying factors.     

A highly divergent mitochondrial DNA lineage of Microtus agrestis in southern Europe

The Mediterranean peninsulas constitute important areas for endemism and intraspecific variation, and are likely places for cryptic biodiversity. We assessed the phylogeographic pattern of field voles (Microtus agrestis) in southern and central Europe by sequence analysis of a 385-bp fragment of the mitochondrial cytochrome b gene in 74 specimens from 44 localities. The majority of samples consisted of skulls collected from owl pellets. The data revealed a highly distinct cytochrome b lineage in an area ranging from Portugal to Hungary. This southern field vole phylogroup differed by a sequence divergence of 5.6-7.1% from the remaining haplotypes, a level of divergence comparable to that found between known Microtus sibling species. However, this ancient phylogeographic break that dates back many glacial cycles has not been recognised previously by either morphology or karyotype. The southern cytochrome b lineage was further divided into two well-defined sublineages that appear to have derived from different glacial refugia in the Iberian Peninsula.     

Colonization history of north European field voles (Microtus agrestis) revealed by mitochondrial DNA

The genetic structure of field vole (Microtus agrestis) populations from northern Europe was examined by restriction fragment length polymorphisms of mitochondrial DNA (mtDNA) in 150 individuals from 67 localities. A total of 83 haplotypes was observed, most of which were rare and highly localized geographically. Overall nucleotide diversity was high (134%), but showed a tendency to decrease with higher latitude. Two major mtDNA lineages differing by 2% in nucleotide sequence were identified. A southern mtDNA lineage was observed in field voles from Britain, Denmark and southern and central Sweden, whereas voles from Finland and northern Sweden belonged to a northern lineage. The strict phylogeographic pattern suggests that the present population generic structure in field voles reflects glacial history: the two groups are derived from different glacial refugia, and recolonized Fennoscandia from two directions. A 150–200-km-wide secondary contact zone between the two mtDNA groups was found in northern Sweden. Distinct phylogeographic substructuring was observed within both major mtDNA groups.