Multiple Glacial Refugia of the Low-Dispersal Ground Beetle Carabus irregularis: Molecular Data Support Predictions of Species Distribution Models

Classical glacial refugia such as the southern European peninsulas were important for species survival during glacial periods and acted as sources of post-glacial colonisation processes. Only recently, some studies have provided evidence for glacial refugia north of the southern European peninsulas. In the present study, we combined species distribution models (SDMs) with phylogeographic analyses (using mitochondrial DNA = mtDNA) to investigate if the cold-adapted, stenotopic and flightless ground beetle species, Carabus irregularis, survived the Last Glacial Maximum (LGM) in classical and/or other refugia. SDMs (for both a western European and for a Carpathian subgroup) were calculated with MAXENT on the basis of 645 species records to predict current and past distribution patterns. Two mtDNA loci (CO1 and ND5, concatenated sequence length: 1785 bp) were analyzed from 91 C. irregularis specimens to reconstruct the phylogeography of Central and eastern European populations and to estimate divergence times of the given lineages. Strong intra-specific genetic differentiation (inter-clade Φ_ST values ranged from 0.92 to 0.99) implied long-term isolation of major clades and subsclades. The high divergence between the nominate subspecies and the Carpathian subspecies C. i. montandoni points to two independent species rather than subspecies (K-2P distance 0.042 ± 0.004; supposed divergence of the maternal lineages dated back 1.6 to 2.5 million years BP) differing not only morphologically but also genetically and ecologically from each other. The SDMs also inferred classical as well as other refugia for C. irregularis, especially north of the Alps, in southeastern Europe and in the Carpathians. The coincidences between the results of both methods confirm the assumption of multiple glacial refugia for the studied species and the usefulness of combining methodological approaches for the understanding of the history of low-dispersal insect species.     

Quaternary refugia in southwestern Iran: insights from two sympatric moth species (Insecta, Lepidoptera)

Alternating glacial and interglacial periods led to range shifts (contractions and expansions), persistence in distinct glacial refugia and extinction events in various temperate organisms. Today, the integrative analysis of molecular markers and spatial distribution models conducted for multiple taxa allows the detection of phylogeographical patterns, thus reconstructing major biogeographical events in their shared evolutionary history. In this study, the effects of past climate change on the evolutionary history of two sympatric moth species (Gnopharmia colchidaria s.l. and G. kasrunensis) and their host plants (Prunus scoparia and P. fenzliana) were inferred for the largely neglected biodiversity hot spot Iran. We complementarily analyzed the population structure of both moth species (187 specimens, based on COI) in congruence with batched species distribution models (SDMs) for all four taxa and for the times of the Last Glacial Maximum (21 ky BP), 6 ky BP and today. Coincidence of SDMs and the distribution of haplotype lineages indicated a shared refugium for the southwestern Zagros Mountains and potential species-specific refugial areas in the southern Caucasus and the Kope-Dagh Mountains. Both moth species experienced past population expansion.     

Taxonomy and Distribution of Freshwater Pearl Mussels (Unionoida: Margaritiferidae) of the Russian Far East

The freshwater pearl mussel family Margaritiferidae includes 13 extant species, which are all listed by IUCN as endangered or vulnerable taxa. In this study, an extensive spatial sampling of Margaritifera spp. across the Russian Far East (Amur Basin, Kamchatka Peninsula, Kurile Archipelago and Sakhalin Island) was conducted for a revision of their taxonomy and distribution ranges. Based on their DNA sequences, shell and soft tissue morphology, three valid species were identified: Margaritifera dahurica (Middendorff, 1850), M. laevis (Haas, 1910) and M. middendorffi (Rosén, 1926). M. dahurica ranges across the Amur basin and some of the nearest river systems. M. laevis is distributed in Japan, Sakhalin Island and the Kurile Archipelago. M. middendorffi was previously considered an endemic species of the Kamchatka. However, it is widespread in the rivers of Kamchatka, Sakhalin Island, the Kurile Islands (across the Bussol Strait, which is the most significant biogeographical boundary within the archipelago), and, likely, in Japan. The Japanese species M. togakushiensis Kondo & Kobayashi, 2005 seems to be conspecific with M. middendorffi because of similar morphological patterns, small shell size (<100 mm long) and overlapped ranges, but it is in need of a separate revision. Phylogenetic analysis reveals that two NW Pacific margaritiferid species, M. laevis and M. middendorffi, formed a monophyletic 18S rDNA clade together with the North American species M. marrianae and M. falcata. The patterns that were found in these Margaritifera spp. are similar to those of freshwater fishes, indicating multiple colonizations of Eastern Asia by different mitochondrial lineages, including an ancient Beringian exchange between freshwater faunas across the Pacific.     

Adding Biotic Interactions into Paleodistribution Models: A Host-Cleptoparasite Complex of Neotropical Orchid Bees

Orchid bees compose an exclusive Neotropical pollinators group, with bright body coloration. Several of those species build their own nests, while others are reported as nest cleptoparasites. Here, the objective was to evaluate whether the inclusion of a strong biotic interaction, such as the presence of a host species, improved the ability of species distribution models (SDMs) to predict the geographic range of the cleptoparasite species. The target species were Aglae caerulea and its host species Eulaema nigrita. Additionally, since A. caerulea is more frequently found in the Amazon rather than the Cerrado areas, a secondary objective was to evaluate whether this species is increasing or decreasing its distribution given South American past and current climatic conditions. SDMs methods (Maxent and Bioclim), in addition with current and past South American climatic conditions, as well as the occurrences for A. caerulea and E. nigrita were used to generate the distribution models. The distribution of A. caerulea was generated with and without the inclusion of the distribution of E. nigrita as a predictor variable. The results indicate A. caerulea was barely affected by past climatic conditions and the populations from the Cerrado savanna could be at least 21,000 years old (the last glacial maximum), as well as the Amazonian ones. On the other hand, in this study, the inclusion of the host-cleptoparasite interaction complex did not statistically improve the quality of the produced models, which means that the geographic range of this cleptoparasite species is mainly constrained by climate and not by the presence of the host species. Nonetheless, this could also be caused by unknown complexes of other Euglossini hosts with A. caerulea, which still are still needed to be described by science.     

Allozyme characterisation of four Mediterranean evergreen oak species

Polymorphism at 11 enzyme loci was used to compare the four evergreen oak species of the Mediterranean Basin and to establish their taxonomical relationships. As a comparison, two evergreen oak species from Tibet, located in the primary distribution area of evergreen oaks, were analysed at the same loci. Cladistic analyses led to the separation of a cluster of four species, namely the three Mediterranean Q. ilex, Q. coccifera and Q. alnifolia and the Tibetan Q. aquifolioides. The other Tibetan species, Q. semicarpifolia, and Q. suber, from the Western Mediterranean Basin, were very distinct genetically from the four other species. These results were not consistent with previous taxonomical treatment, based on morphology alone, which classified Q. ilex and Q. coccifera in “ilex group” (=subgenus Sclerophyllodrys O. Schwartz), and the other four species in “cerris group” (subgenus Cerris, according to Schwartz). Allozyme variation has thus provided useful information to clarify the taxonomy of evergreen oaks.     

Using genetic diversity information to establish core collections of Stylosanthes capitata and Stylosanthes macrocephala

Stylosanthes species are important forage legumes in tropical and subtropical areas. S. macrocephala and S. capitata germplasm collections that consist of 134 and 192 accessions, respectively, are maintained at the Brazilian Agricultural Research Corporation Cerrados (Embrapa-Cerrados). Polymorphic microsatellite markers were used to assess genetic diversity and population structure with the aim to assemble a core collection. The mean values of H_O and H_E for S. macrocephala were 0.08 and 0.36, respectively, whereas the means for S. capitata were 0.48 and 0.50, respectively. Roger’s genetic distance varied from 0 to 0.83 for S. macrocephala and from 0 to 0.85 for S. capitata. Analysis with STRUCTURE software distinguished five groups among the S. macrocephala accessions and four groups among those of S. capitata. Nei’s genetic diversity was 27% in S. macrocephala and 11% in S. capitata. Core collections were assembled for both species. For S. macrocephala, all of the allelic diversity was represented by 23 accessions, whereas only 13 accessions were necessary to represent all allelic diversity for S. capitata. The data presented herein evidence the population structure present in the Embrapa-Cerrados germplasm collections of S. macrocephala and S. capitata, which may be useful for breeding programs and germplasm conservation.     

The Taxonomic Status of Mazama bricenii and the Significance of the Táchira Depression for Mammalian Endemism in the Cordillera de Mérida,Venezuela

We studied the taxonomy and biogeography of Mazama bricenii, a brocket deer classified as Vulnerable by the IUCN, drawing on qualitative and quantitative morphology and sequences of the mitochondrial cytochrome-b gene. We used Ecological Niche Modeling (ENM) to evaluate the hypothesis that M. bricenii of the Venezuelan Cordillera de Mérida (CM) might have become isolated from populations of its putative sister species, Mazama rufina, in the Colombian Cordillera Oriental (CO). This hypothesis assumes that warm, dry climatic conditions in the Táchira Depression were unsuitable for the species. Our analyses did not reveal morphological differences between specimens geographically attributable to M. bricenii and M. rufina, and phylogenetic analyses of molecular data recovered M. bricenii nested within the diversity of M. rufina. These results indicate that M. bricenii should be regarded as a junior synonym of M. rufina. ENM analyses revealed the existence of suitable climatic conditions for M. rufina in the Táchira Depression during the last glacial maximum and even at present, suggesting that gene flow between populations in the CO and CM may have occurred until at least the beginning of the current interglacial period and may continue today. Because this pattern might characterize other mammals currently considered endemic to the CM, we examined which of these species match two criteria that we propose herein to estimate if they can be regarded as endemic to the CM with confidence: (1) that morphological or molecular evidence exists indicating that the putative endemic taxon is distinctive from congeneric populations in the CO; and (2) that the putative endemic taxon is restricted to either cloud forest or páramo, or both. Only Aepeomys reigi, Cryptotis meridensis, and Nasuella meridensis matched both criteria; hence, additional research is necessary to assess the true taxonomic status and distribution of the remaining species thought to be CM endemics.     

Evolutionary history of Caribbean species of Myotis,with evidence of a third Lesser Antillean endemic

Currently, four species of Myotis are known from the islands of the Caribbean (Myotis dominicensis, M. martiniquensis, M. nesopolus, and M. nigricans). Myotis dominicensis and M. martiniquensis are endemic to the Lesser Antilles, whereas M. nesopolus and M. nigricans are considered conspecific with mainland populations. Recent phylogenetic and phylogeographic studies provided hypotheses regarding the origin and diversification of M. dominicensis and M. martiniquensis. However, these studies focused primarily on convergent morphology or distribution patterns of this genus and not on the evolutionary history of Caribbean Myotis. Here, we explore variation across multiple datasets generated from Caribbean Myotis. We present morphologic and genetic (mitochondrial and nuclear) data from an extensive sample of Caribbean Myotis species, including the previously unsampled taxa M. martiniquensis nyctor and M. nesopolus. Our data indicate that the historically recognized subspecies M. m. nyctor is genetically and morphologically distinct from M. martiniquensis, warranting recognition of a third Caribbean endemic—Myotis nyctor. Moreover, we provide evidence of unrecognized species-level variation in Caribbean and northern South American populations of Myotis.     

Demography of the columnar cactus Neobuxbaumia macrocephala: a comparative approach using population projection matrices

Neobuxbaumia macrocephala is a long-lived columnar cactus endemic to the Tehuacán-Cuicatlán Valley in south-central Mexico. This plant has a very restricted distribution and few recruitment events have been detected in its populations. In this study, we analyze the N. macrocephala demographic pattern using a projection matrix in order to determine the main limiting factors of this species. To accomplish this goal, we compare our results with those obtained for another species of the same genus, N. tetetzo. Considering that both species inhabit the same valley, we believe that this comparative study will offer insights into the main demographic limitations of N. macrocephala. Results showed that these species of columnar cacti have similar demographic patterns in which survival is the process with the highest relative contribution to λ, followed by growth and reproduction. Of all the life cycle stages, seeds and seedlings have the lowest survival probabilities due to a high mortality caused by seed predation and effects of direct solar radiation on germinated seeds. The estimated growth rates indicate that populations of these species of Neobuxbaumia are in a numerical equilibrium. With respect to reproduction, N. macrocephala produce a lower number of seeds per plant than N. tetetzo. This low level of sexual reproduction may decrease the probability of establishment of new individuals in N. macrocephala populations. It is suggested that pollen limitation and pre-dispersal seed predation could be some factors that limit the distribution and abundance of this columnar cactus. This revised version was published online in June 2006 with corrections to the Cover Date.     

Palaeoclimatic models help to understand current distribution of Caucasian forest species

Spatial and temporal constraints on dispersal explain the absence of species from areas with potentially suitable conditions. Previous studies have shown that post‐glacial recolonization has shaped the current ranges of many species, yet it is not completely clear to what extent interspecific differences in range size depend on different dispersal rates. The inferred boundaries of glacial refugia are difficult to validate, and may bias spatial distribution models (SDMs) that consider post‐glacial dispersal constraints. We predicted the current distribution of 12 Caucasian forest plants and animals, factoring in the effective geographical distance from inferred glacial refugia as an additional predictor. To infer glacial refugia, we tested the transferability of the current SDMs based on the distribution of climatic variables, and projected the most transferable ones onto two climate scenarios simulated for the Last Glacial Maximum (LGM). We then calculated least‐cost distances from the inferred refugia, using elevation as a friction surface, and recalculated the current SDMs incorporating the distances as an additional variable. We compared the predictive powers of the initial with the final SDMs. The palaeoclimatic simulation that best matched the distribution of species was assumed to represent the closest fit to the true palaeoclimate. SDMs incorporating refugial distance performed significantly better for all but one studied species, and the Model for Interdisciplinary Research on Climate (MIROC) climatic simulation provided a more convincing pattern of the LGM climate than the Community Climate System Model (CCSM) simulation. Our results suggest that the projection of suitable habitat models onto past climatic conditions may yield realistic boundaries of glacial refugia, and that the current distribution of forest species in the study region is strongly associated with locations of former refugia. We inferred six major forest refugia throughout western Asia: (1) Colchis; (2) western Anatolia; (3) western Taurus; (4) the upper reaches of the Tigris River; (5) the Levant; and (6) the southern Caspian basin. The boundaries of the modelled refugia were substantially broader than the refugia boundaries inferred solely from pollen records. Thus, our method could be used to: (1) improve models of current species distributions by considering the dispersal histories of the species; and (2) validate alternative reconstructions of palaeoclimate with current distribution data. © 2011 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, 105 , 231–248.     

Identification of New Genomospecies in the Mycobacterium terrae Complex

Members of the Mycobacterium terrae complex are slow-growing, non-chromogenic acid-fast bacilli found in the natural environment and occasionally in clinical material. These genetically closely-related members are difficult to differentiate by conventional phenotypic and molecular tests. In this paper we describe the use of whole genome data for the identification of four strains genetically similar to Mycobacterium sp. JDM601, a newly identified member of the M. terrae complex. Phylogenetic information from the alignment of genome-wide orthologous genes and single nucleotide polymorphisms show consistent clustering of the four strains together with M. sp. JDM601 into a distinct clade separate from other rapid and slow growing mycobacterial species. More detailed inter-strain comparisons using average nucleotide identity, tetra-nucleotide frequencies and analysis of synteny indicate that our strains are closely related to but not of the same species as M. sp. JDM601. Besides the 16S rRNA signature described previously for the M. terrae complex, five more hypothetical proteins were found that are potentially useful for the rapid identification of mycobacterial species belonging to the M. terrae complex. This paper illustrates the versatile utilization of whole genome data for the delineation of new bacterial species and introduces four new genomospecies to add to current members in the M. terrae complex.     

Past,present, and future predictions on the suitable habitat of the Slender racer (Orientocoluber spinalis) using species distribution models

Species distribution models (SDMs) across past, present, and future timelines provide insights into the current distribution of these species and their reaction to climate change. Specifically, if a species is threatened or not well‐known, the information may be critical to understand that species. In this study, we computed SDMs for Orientocoluber spinalis, a monotypic snake genus found in central and northeast Asia, across the past (last interglacial, last glacial maximum, and mid‐Holocene), present, and future (2070s). The goal of the study was to understand the shifts in distribution across time, and the climatic factors primarily affecting the distribution of the species. We found the suitable habitat of O. spinalis to be persistently located in cold‐dry winter and hot summer climatic areas where annual mean temperature, isothermality, and annual mean precipitation were important for suitable habitat conditions. Since the last glacial maximum, the suitable habitat of the species has consistently shifted northward. Despite the increase in suitable habitat, the rapid alterations in weather regimes because of climate change in the near future are likely to greatly threaten the southern populations of O. spinalis, especially in South Korea and China. To cope with such potential future threats, understanding the ecological requirements of the species and developing conservation plans are urgently needed.     

Continental Refugium in the Mongolian Plateau during Quaternary Glacial Oscillations: Phylogeography and Niche Modelling of the Endemic Desert Hamster,Phodopus roborovskii

The Mongolian Plateau (MP), which is situated in the interior of Asia and possesses a typical continental climate, experienced harsh climatic conditions during the Quaternary glacial fluctuations. Although these events likely had huge impacts on the local animal populations, the current effects have hardly been explored. To investigate whether the MP supported a refugium along an oceanic-continental gradient (ROCG), and whether this refugium was glacial or interglacial, we investigated the demographic and phylogeographic history of an endemic mammal species, the desert hamster Phodopus roborovskii. We reconstructed the demographic variation, the phylogeographic diffusion, and modelled the potential habitat during historical periods. The genetic diversity in the MP was the highest among all the localities, and the MP was a suitable habitat throughout the modelled historical periods. A phylogeographic diffusion analysis emphasized the importance of the MP as the centre of origin, preservation and spread for P. roborovskii. The homogeneous landscape provided the opportunity for a wide gene flow, which resulted in low resolution of the phylogenetic relationships. Moreover, P. roborovskii was favoured by the interglacial condition, with both its demographical and geographical ranges expanded within the interglacial periods. The range variation from the Last Glacial Maximum to the current condition reflects a distinct longitudinal shift, while both ranges largely contracted from that of the Last Interglacial. Our results support that the MP served as a refugium and spread centre for P. roborovskii during the Quaternary climate fluctuations. The interglacial expansion and the longitudinal shifts highlighted the important effects of precipitations on the distribution range of species adapted to arid and semi-arid during glacial oscillations.     

Conservation Genetics of Kangaroo Mice, Genus Microdipodops

The two currently recognized species of kangaroo mice, Microdipodops megacephalus and M. pallidus, inhabit sandy soils of the Great Basin Desert in western North America. Given their habitat specificity and the fluctuating climate throughout the Pleistocene, kangaroo mice likely endured a turbulent biogeographic history that resulted in disjunct distributions and isolation of genetic lineages. Recent phylogenetic investigations using mitochondrial data have revealed several mitochondrial clades within this genus that may represent cryptic species. These mitochondrial clades are genetically unique, occupy relatively small distributions, and, as such, may be at an increased risk of extinction due to climate change and extensive recent habitat alteration. Herein, we apply haplotype network, population genetic, and historical demographic analyses to mitochondrial data of each Micropdipodops species and mitochondrial clade to assess conservation genetics within kangaroo mice. Results indicate that each mitochondrial clade is a distinct lineage with little to no gene flow occurring among clades. Additionally, historical demographic analyses support past population expansions and identify locations of past refugium for each distinct lineage. Although mitochondrial data indicate that the clades appear to be in approximate genetic equilibrium and have not suffered any extreme bottlenecks over time, there is still concern for the survival of smaller and more vulnerable Microdipodops subpopulations due to impending habitat threats in the Great Basin Desert.     

Comparative and Evolutionary Analyses of Meloidogyne spp. Based on Mitochondrial Genome Sequences

Molecular taxonomy and evolution of nematodes have been recently the focus of several studies. Mitochondrial sequences were proposed as an alternative for precise identification of Meloidogyne species, to study intraspecific variability and to follow maternal lineages. We characterized the mitochondrial genomes (mtDNAs) of the root knot nematodes M. floridensis, M. hapla and M. incognita. These were AT rich (81–83%) and highly compact, encoding 12 proteins, 2 rRNAs, and 22 tRNAs. Comparisons with published mtDNAs of M. chitwoodi, M. incognita (another strain) and M. graminicola revealed that they share protein and rRNA gene order but differ in the order of tRNAs. The mtDNAs of M. floridensis and M. incognita were strikingly similar (97–100% identity for all coding regions). In contrast, M. floridensis, M. chitwoodi, M. hapla and M. graminicola showed 65–84% nucleotide identity for coding regions. Variable mitochondrial sequences are potentially useful for evolutionary and taxonomic studies. We developed a molecular taxonomic marker by sequencing a highly-variable ~2 kb mitochondrial region, nad5-cox1, from 36 populations of root-knot nematodes to elucidate relationships within the genus Meloidogyne. Isolates of five species formed monophyletic groups and showed little intraspecific variability. We also present a thorough analysis of the mitochondrial region cox2-rrnS. Phylogenies based on either mitochondrial region had good discrimination power but could not discriminate between M. arenaria, M. incognita and M. floridensis.     

The leaf turtle population of Phnom Kulen National Park (northwestern Cambodia) has genetic and morphological signatures of hybridization

Cambodia is known to harbour three distinct species of Southeast Asian leaf turtles (Cyclemys spp.), which are heavily traded and common in seizures of wildlife. Confiscated leaf turtles are often released to natural habitats. Thus, an exact knowledge of the distribution of the individual species is of great importance to avoid the introduction of non‐native turtles, posing the risk of competition and hybridization. In this study, we examine a recently discovered leaf turtle population from Phnom Kulen National Park using external morphology, 17 unlinked microsatellite loci and the mitochondrial cytochrome b gene. Leaf turtles from the Phnom Kulen National Park morphologically resemble C. oldhamii, but harbour mitochondrial haplotypes of C. atripons. With respect to microsatellite loci, the turtles are distinct from C. atripons. Unfortunately no material of C. oldhamii was available from Cambodia. We propose that the Phnom Kulen population represents either a natural hybrid swarm of C. atripons and C. oldhamii or a distinct undescribed species with introgressed mitochondria of C. atripons. This underlines that genetic differentiation of wild leaf turtle populations in Cambodia is complex and suggests that this differentiation pattern becomes increasingly threatened by translocations of confiscated individuals. For drawing a definite conclusion about the taxonomic status of the Phnom Kulen population, denser sampling of other Cambodian leaf turtle populations would be required, in particular of C. oldhamii.     

Phyllachora species infecting maize and other grass species in the Americas represents a complex of closely related species

The genus Phyllachora contains numerous obligate fungal parasites that produce raised, melanized structures called stromata on their plant hosts referred to as tar spot. Members of this genus are known to infect many grass species but generally do not cause significant damage or defoliation, with the exception of P. maydis which has emerged as an important pathogen of maize throughout the Americas, but the origin of this pathogen remains unknown. To date, species designations for Phyllachora have been based on host associations and morphology, and most species are assumed to be host specific. We assessed the sequence diversity of 186 single stroma isolates collected from 16 hosts representing 15 countries. Samples included both herbarium and contemporary strains that covered a temporal range from 1905 to 2019. These 186 isolates were grouped into five distinct species with strong bootstrap support. We found three closely related, but genetically distinct groups of Phyllachora are capable of infecting maize in the United States, we refer to these as the P. maydis species complex. Based on herbarium specimens, we hypothesize that these three groups in the P. maydis species complex originated from Central America, Mexico, and the Caribbean. Although two of these groups were only found on maize, the third and largest group contained contemporary strains found on maize and other grass hosts, as well as herbarium specimens from maize and other grasses that include 10 species of Phyllachora. The herbarium specimens were previously identified based on morphology and host association. This work represents the first attempt at molecular characterization of Phyllachora species infecting grass hosts and indicates some Phyllachora species can infect a broad range of host species and there may be significant synonymy in the Phyllachora genus.     

Mitochondrial DNA variation and Quaternary range dynamics in the endangered Yellow Spotted Mountain Newt,Neurergus derjugini (Caudata,Salamandridae)

In temperate regions of the Earth Pleistocene, climatic fluctuations significantly influenced distribution of species. However, little is known on how glacial and interglacial cycles affected range dynamics of the species occupying lower latitudes. In this study, we investigated mitochondrial DNA (mtDNA) variation and reconstructed the potential current and past (during the mid‐Holocene, 6 ka BP, and the Last Glacial Maximum, LGM, 21 ka BP) distribution of Neurergus derjugini, an endangered amphibian species endemic to the mid‐Zagros Mountains in Iran and Iraq. Six haplotypes identified in the control region (D‐loop) form a well‐supported monophyletic clade, distinct from other Neurergus species and revealing a sister relationship to Neurergus kaiseri. Nucleotide diversity quantifying mean divergence between the sequences is low and does not support the recognition of distinct evolutionary lineages in Neurergus derjugini. The landscape connectivity analysis and the haplotype parsimony network reveal higher gene flow rate between the breeding streams in the southern part of the range, while the northern populations are more isolated. The potential distribution of Neurergus derjugini is restricted to valleys close to mountain tops, wherein very high elevations and dry habitats appear to be unsuitable. During the mid‐Holocene and LGM conditions, the range of the species may have been more extended and shifted to lower elevations. These findings show retraction of the Neurergus derjugini range during the Quaternary and indicate that range dynamics of the species occupying lower latitudes may not follow a scenario of glacial retraction and postglacial expansion.     

Gene flow across secondary contact zones of the Emys orbicularis complex in the Western Mediterranean and evidence for extinction and re‐introduction of pond turtles on Corsica and Sardinia (Testudines: Emydidae)

European pond turtles represent a phylogeographically deeply structured complex of distinct taxa. Here, we use mitochondrial DNA sequences (cytochrome b gene) and eight polymorphic microsatellite loci to investigate genetic differentiation and gene flow of Sicilian, Corsican and Sardinian pond turtles and of subspecies involved in two secondary contact zones in the Pyrenean region and Southern Italy. Mitochondrial and microsatellite differentiation is largely concordant in populations from the core regions of the distribution ranges of the studied taxa. Both marker systems provide no evidence for gene flow between Sicilian pond turtles (Emys trinacris) and Southern Italian subspecies of E. orbicularis. By contrast, in the contact zones limited gene flow occurs between distinct subspecies of E. orbicularis. Although the Southern Italian contact zone is significantly older than the Pyrenean contact zone of Holocene age, patterns of asymmetric introgression are similar. Introgressive hybridization leads to the exchange of mitochondria, but microsatellite data indicate only a few individuals with mixed ancestry. This suggests that incipient isolating mechanisms maintain largely discrete nuclear genomic gene pools. Furthermore, this implies that Southern Italy acted as a hotspot rather than as a melting pot of genetic diversity during the last glacial. Pond turtles from Corsica and Sardinia are not differentiated from continental populations of the subspecies E. o. galloitalica, neither in the mitochondrial nor in the quickly evolving microsatellite markers. As the fossil record argues for a continuous presence of pond turtles on both islands since the Middle Pleistocene, this suggests that the native island populations became extinct and the extant turtles were later introduced by prehistoric settlers. The lack of genetic differentiation of pond turtles from Corsica and Sardinia supports the view that the subspecies described from these islands are not valid.     

Back from the Brink: The Holocene History of the Carpathian Barbel Barbus carpathicus

As a result of specific adaptations and habitat preferences strongly rheophilic fish species may show high levels of endemism. Many temperate rheophilic fish species were subjected to a series of range contractions during the Pleistocene, and then successfully expanded during the Holocene, colonising previously abandoned areas. The Carpathian barbel (Barbus carpathicus Kotlík, Tsigenopoulos, Ráb et Berrebi 2002) occurs in the montane streams in three basins of the main Central European rivers in the northern part of the Carpathian range. We used genetic variation within 3 mitochondrial and 9 microsatellite loci to determine a pattern of postglacial expansion in B. carpathicus. We found that overall genetic variation within the species is relatively low. Estimate of time to the most recent common ancestor (tMRCA) of mitochondrial sequences falls within the Holocene. The highest levels of genetic variation found in upper reaches of the Tisa river in the Danube basin suggest that glacial refugia were located in the south-eastern part of the species range. Our data suggest that the species crossed different watersheds at least six times as three genetically distinct groups (probably established in different expansion episodes) were found in northern part of the species range. Clines of genetic variation were observed in both the Danube and Vistula basins, which probably resulted from subsequent bottlenecks while colonizing successive habitats (south eastern populations) or due to the admixture of genetically diverse individuals to a previously uniform population (Vistula basin). Therefore, B. carpathicus underwent both demographic breakdowns and expansions during the Holocene, showing its distribution and demography are sensitive to environmental change. Our findings are important in the light of the current human-induced habitats alterations.