European colonization by the spined loach (Cobitis taenia) from Ponto-Caspian refugia based on mitochondrial DNA variation

Whole-genome analysis reveals phylogenetic and demographic history of Eurasian perch

The contemporary diversity and distribution of species are shaped by their evolutionary and ecological history. This can be deciphered with the help of phylogenetic and demographic analysis methods, ideally combining and supplementing information from mitochondrial and nuclear genomes. In this study, we investigated the demographic history of Eurasian perch (Perca fluviatilis), a highly adaptable teleost with a distribution range across Eurasia. We combined whole-genome resequencing data with available genomic resources to analyse the phylogeny, phylogeography, and demographic history of P. fluviatilis populations from Europe and Siberia. We identified five highly diverged evolutionary mtDNA lineages, three of which show a strong signal of admixture in the Baltic Sea region. The estimated mean divergence time between these lineages ranged from 0.24 to 1.42 million years. Based on nuclear genomes, two distinct demographic trajectories were observed in European and Siberian samples reflecting contrasting demographic histories ca. 30,000-100,000 years before the present. A comparison of mtDNA and nuclear DNA evolutionary trees and AMOVA revealed concordances, as well as incongruences, between the two types of data, most likely reflecting recent postglacial colonization and hybridization events. Overall, our findings demonstrate the power and usefulness of genome-wide information for delineating historical processes that have shaped the genome of P. fluviatilis. We also highlight the added value of data-mining existing transcriptomic resources to complement novel sequence data, helping to shed light on putative glacial refugia and postglacial recolonization routes.

Increased spatial resolution of sampling in the Carpathian basin helps to understand the phylogeny of central European stream-dwelling gudgeons

Background

Phylogenetic studies of widespread European fish species often do not completely cover their entire distribution area, and some areas are often excluded from analyses than others. For example, Carpathian stocks are often omitted from these surveys or are under-represented in the samples. However, this area served as an extra-Mediterranean refugia for many species; therefore, it is assumed that fish stocks here may show special phylogenetic features. For this reason, increased spatial resolution of sampling, namely revealing genetic information from unexamined Carpathian areas within the range of doubtful taxa, may help us better understand their phylogenetic features. To test this hypothesis, a phylogenetic investigation using a partial mtCR sequence data was conducted on 56 stream-dwelling freshwater fish (Gobio spp.) individuals collected from 11 rivers of the data-deficient Southeastern Carpathian area. Moreover, we revieved the available phylogenetic data of Middle-Danubian stream-dwelling gudgeon lineages to delineate their distribution in the area.

Results

Seven out of the nine detected haplotypes were newly described, suggesting the studied area hosts distinct and diverse Gobio stocks. Two valid species (G. obtusirostris, G. gobio), and a haplogroup with doubtful phylogenetic position” G. sp. 1" were detected in the area, showing a specific spatial distribution pattern. The distribution of the detected lineages in the Middle-Danubian area correspond with recent and paleo hydrogeographic features, at the same time mainly on their bordering areas show considerable overlap.

Conclusions

Despite the relatively limited geographic range of the study, our results provide important information which contributes to a better understanding of the phylogenetic, taxonomic and distribution features of Central European gudgeons. The genetically confirmed distribution data of the indicated lineages corresponds well with the recent and near-recent hydrogeographic features of the area, and assumes several hybrid zones in the Carpathian Basin. Additionally, the results show that the middle and lower Danubian watershed cannot be excluded from the range of G. gobio. Moreover, the” G. sp. 1", is slightly differentiated but phylogenetically distinct entity, and is the only Gobio taxa thus far detected in the middle and lower Tisza-basin. However, further investigations are necessary to clarify the taxonomic position of this group.

Phylogeography and potential glacial refugia of terrestrial gastropod Faustina faustina (Rossmässler, 1835) (Gastropoda: Eupulmonata: Helicidae) inferred from molecular data and species distribution models

Faustina faustina is a conchologically highly diverse forest gastropod with several morphological forms. It is a Carpathian species, but it also occurs in northern isolated localities, where it was probably introduced. We performed the first phylogeographic analysis of 22 populations, based on three molecular markers: COI, ITS-2, and 28S rRNA. Genetic data were complemented by paleo-distribution models of spatial occupancy during the Last Glacial Maximum to strengthen inferences of refugial areas. We discovered high genetic variability of COI sequences with p-distances between haplotypes ranged from 0.2 to 18.1% (6.3–16.6% between clades). For nuclear markers, a haplotype distribution pattern was revealed. Species distribution models indicated a few potential refugia in the Carpathians, with the most climatically stable and largest areas in the Southern Carpathians. In some climate scenarios, putative microrefugia were also predicted in the Western and Eastern Carpathians, and in the Apuseni Mts. Our results suggest the glacial in situ survival of F. faustina and its Holocene expansion in the Sudetes. Although our genetic data as well as shell phenotypes showed considerable variation within and between studied populations, the molecular species delimitation approaches still imply only one single species. Our study contributes to the understanding of the impact of processes on shaping contemporary population genetic structure and diversity in low-dispersal, forest species.

Cytogenetic Characterization of Seven Novel satDNA Markers in Two Species of Spined Loaches (Cobitis) and Their Clonal Hybrids

Interspecific hybridization is a powerful evolutionary force. However, the investigation of hybrids requires the application of methodologies that provide efficient and indubitable identification of both parental subgenomes in hybrid individuals. Repetitive DNA, and especially the satellite DNA sequences (satDNA), can rapidly diverge even between closely related species, hence providing a useful tool for cytogenetic investigations of hybrids. Recent progress in whole-genome sequencing (WGS) offers unprecedented possibilities for the development of new tools for species determination, including identification of species-specific satDNA markers. In this study, we focused on spined loaches (Cobitis, Teleostei), a group of fishes with frequent interspecific hybridization. Using the WGS of one species, C. elongatoides, we identified seven satDNA markers, which were mapped by fluorescence in situ hybridization on mitotic and lampbrush chromosomes of C. elongatoides, C. taenia and their triploid hybrids (C. elongatoides × 2C. taenia). Two of these markers were chromosome-specific in both species, one had centromeric localization in multiple chromosomes and four had variable patterns between tested species. Our study provided a novel set of cytogenetic markers for Cobitis species and demonstrated that NGS-based development of satDNA cytogenetic markers may provide a very efficient and easy tool for the investigation of hybrid genomes, cell ploidy, and karyotype evolution.

High cryptic diversity of bitterling fish in the southern West Palearctic

South-east Europe, along with the adjacent region of south-west Asia, is an important biodiversity hotspot with high local endemism largely contributed by contemporary continental lineages that retreated to southern refugia during colder Quaternary periods. We investigated the genetic diversity of the European bitterling fish (Rhodeus amarus) species complex (Cyprinidae) across its range in the western Palearctic, but with a particular emphasis in the region of Balkan, Pontic and Caspian refugia. We genotyped 12 polymorphic microsatellite loci and a partial sequence of mitochondrial gene cytochrome b (CYTB) for a set of 1,038 individuals from 60 populations. We used mtDNA sequences to infer phylogenetic relationships and historical demography, and microsatellite markers to describe fine-scale genetic variability and structure. Our mtDNA analysis revealed six well-supported lineages, with limited local co-occurrence. Two lineages are distributed throughout central and western Europe (lineages "A" and "B"), with two zones of secondary contact. Another two lineages were restricted to the Ponto-Aegean region of Greece (lineages "C" and "D") and the final two lineages were restricted south of the Caucasus mountains (lineage "E" from the Black Sea watershed and lineage "F" from the Caspian watershed). A signal of recent expansion was revealed in the two widespread lineages and the Ponto-Aegean lineage "C". The geographic distribution of clusters detected by nuclear microsatellites corresponded well with mitochondrial lineages and demonstrated finely sub-structured populations. A profound population structure suggested a significant role of genetic drift in differentiation among lineages. Lineage divergence in the Ponto-Aegean and Caspian regions are substantial, supporting the validity of two described endemic species (Rhodeus meridionalis as lineage "D" and Rhodeus colchicus as lineage "E") and invite taxonomic evaluation of the other two southern lineages (Thracean "C" and Caspian "F").

Diversity of spined loaches from Asia Minor in a phylogenetic context (Teleostei: Cobitidae)

Accurate determination of species diversity in areas of high endemicity, particularly those lacking comprehensive systematic knowledge, represents a challenge for both taxonomists and conservationists. This need is particularly evident in areas greatly affected by anthropogenic disturbances such as the Eastern Mediterranean and its freshwater environments. To improve our knowledge of Eastern Mediterranean freshwater fishes, we phylogenetically studied Western Palearctic Cobitis species, focusing on those found in Turkey. Overall, our results provide a robust framework to assess the number of species of Cobitis. Phylogenetic reconstructions based on mitochondrial (cyt b) and nuclear (RAG1) sequences show seven major clades (Clades 1-7) grouping all Western Palearctic Cobitis species, except C. melanoleuca. In general, each major clade comprises Cobitis species that inhabit geographically close areas and have similar secondary sexual characters. Multiple divergent lineages were identified in our analyses, some of which were highly divergent such as the ones inhabiting Turkish freshwaters. Moreover, in some analyses, several of the identified lineages were incongruent with a priori defined species. Furthermore, our analyses identified eight potentially new candidate species, six that had been suggested in previous studies and two that are reported here for the first time. Our results reveal Turkey as the area with the greatest diversity of spined loaches in the Mediterranean.

Comparing RADseq and microsatellites to infer complex phylogeographic patterns, an empirical perspective in the Crucian carp, Carassius carassius, L

The conservation of threatened species must be underpinned by phylogeographic knowledge. This need is epitomized by the freshwater fish Carassius carassius, which is in decline across much of its European range. Restriction site-associated DNA sequencing (RADseq) is increasingly used for such applications; however, RADseq is expensive, and limitations on sample number must be weighed against the benefit of large numbers of markers. This trade-off has previously been examined using simulation studies; however, empirical comparisons between these markers, especially in a phylogeographic context, are lacking. Here, we compare the results from microsatellites and RADseq for the phylogeography of C. carassius to test whether it is more advantageous to genotype fewer markers (microsatellites) in many samples, or many markers (SNPs) in fewer samples. These data sets, along with data from the mitochondrial cytochrome b gene, agree on broad phylogeographic patterns, showing the existence of two previously unidentified C. carassius lineages in Europe: one found throughout northern and central-eastern European drainages and a second almost exclusively confined to the Danubian catchment. These lineages have been isolated for approximately 2.15 m years and should be considered separate conservation units. RADseq recovered finer population structure and stronger patterns of IBD than microsatellites, despite including only 17.6% of samples (38% of populations and 52% of samples per population). RADseq was also used along with approximate Bayesian computation to show that the postglacial colonization routes of C. carassius differ from the general patterns of freshwater fish in Europe, likely as a result of their distinctive ecology.

The population history of Garra orientalis (Teleostei: Cyprinidae) using mitochondrial DNA and microsatellite data with approximate Bayesian computation

Background

The South China landmass has been characterized by a complex geological history, including mountain lifting, climate changes, and river capture/reversal events. To determine how this complexity has influenced the landmass’s phylogeography, our study examined the phylogeography of Garra orientalis, a cyprinid widely distributed in South China, using sequences from the mitochondrial DNA control region and cytochrome b gene (1887 bp) and polymorphisms of thirteen microsatellite loci.

Results

In total, 157 specimens were collected from eight populations. All 88 mtDNA haplotypes were identified as belonging to three major lineages, and these lineages were almost allopatric in their distributions. The results of a statistical dispersal-vicariance analysis suggested that the ancestral populations of G. orientalis were distributed south of the Yunkai Mountains, including on Hainan Island. The mtDNA data revealed a strong relationship between phylogeny and geography. In the microsatellite analysis, a total of 339 alleles with an average of 26 alleles per locus were observed across thirteen microsatellite loci. A clustering algorithm for microsatellite data revealed an admixture-like genetic structure. Although the mtDNA and microsatellite data sets displayed a discordant population structure, the results of an approximate Bayesian computation approach showed that these two markers revealed congruent historical signals. The population history of G. orientalis reflects vicariance events and dispersal related to the complex geological history of South China.

Conclusion

Our results (i) found that the discordances between mtDNA and microsatellite markers were accounted for by admixtures; (ii) showed that the Wuzhishan and Yinggeling mountain ranges and Qiongzhou Strait were important barriers limiting gene exchange between populations on both sides; (iii) indicated that during glaciation and inter-glacial periods, the strait and continental shelves were exposed and sank, which contributed with the dispersion and differentiation of populations; and (iv) displayed that the admixtures between lineages took place in coastal populations and then colonized the tributaries of the Pearl River.

Genomic population structure of freshwater-resident and anadromous ide (Leuciscus idus) in north-western Europe

Climate change experts largely agree that future climate change and associated rises in oceanic water levels over the upcoming decades, will affect marine salinity levels. The subsequent effects on fish communities in estuarine ecosystems however, are less clear. One species that is likely to become increasingly affected by changes in salinity is the ide (Leuciscus idus). The ide is a stenohaline freshwater fish that primarily inhabits rivers, with frequent anadromous behavior when sea salinity does not exceed 15%. Unlike most other anadromous Baltic Sea fish species, the ide has yet to be subjected to large-scale stocking programs, and thus provides an excellent opportunity for studying the natural population structure across the current salinity gradient in the Danish Belts. To explore this, we used Genotyping-by-Sequencing to determine genomic population structure of both freshwater resident and anadromous ide populations in the western Baltic Sea region, and relate the results to the current salinity gradient and the demographic history of ide in the region. The sample sites separate into four clusters, with all anadromous populations in one cluster and the freshwater resident populations in the remaining three. Results demonstrate high level of differentiation between sites hosting freshwater resident populations, but little differentiation among anadromous populations. Thus ide exhibit the genomic population structure of both a typical freshwater species, and a typical anadromous species. In addition to providing a first insight into the population structure of north-western European ide, our data also (1) provide indications of a single illegal introduction by man; (2) suggest limited genetic effects of heavy pollution in the past; and (3) indicate possible historical anadromous behavior in a now isolated freshwater population.

Phylogeographical analysis of a cold-temperate freshwater fish, the Amur sleeper (Perccottus glenii) in the Amur and Liaohe River basins of Northeast Asia

The Amur sleeper Perccottus glenii (Perciformes, Gobioidei, Odontobutidae) is well known as an invasive fish in the river basins of Eastern and Central Europe, but its genetic background is unavailable across its native habitats in northeast Asia. In this study, we used the mitochondrial cytochrome b gene by sampling 19 populations of P. glenii across its native distributional areas of Liaohe and Amur River basins to explore its evolutionary history. Phylogenetic analyses identified three major clades within P. glenii, among which Clade A and Clade B were co-distributed in the Liaohe and Amur River basins, and Clade C was restricted to the latter. Molecular dating showed that the splits of Clades A, B and C have happened in the late Early-early Middle Pleistocene and the most recent common ancestors of these clades have been presented in the late Middle-early Late Pleistocene. The P. glenii showed very high levels of genetic structure among populations (ΦST = 0.801), probably due to the characters of its life histories with very limited dispersal ability. The admixture of different clades in some populations of P. glenii probably reflects historical secondary contact. These findings indicate that Pleistocene climatic oscillation and river capture were major determinants for genetic variations and evolutionary history of the P. glenii.

Low genetic diversity in wide-spread Eurasian liver fluke Opisthorchis felineus suggests special demographic history of this trematode species

Opisthorchis felineus or Siberian liver fluke is a trematode parasite (Opisthorchiidae) that infects the hepato-biliary system of humans and other mammals. Despite its public health significance, this wide-spread Eurasian species is one of the most poorly studied human liver flukes and nothing is known about its population genetic structure and demographic history. In this paper, we attempt to fill this gap for the first time and to explore the genetic diversity in O. felineus populations from Eastern Europe (Ukraine, European part of Russia), Northern Asia (Siberia) and Central Asia (Northern Kazakhstan). Analysis of marker DNA fragments from O. felineus mitochondrial cytochrome c oxidase subunit 1 and 3 (cox1, cox3) and nuclear rDNA internal transcribed spacer 1 (ITS1) sequences revealed that genetic diversity is very low across the large geographic range of this species. Microevolutionary processes in populations of trematodes may well be influenced by their peculiar biology. Nevertheless, we suggest that lack of population genetics structure observed in O. felineus can be primarily explained by the Pleistocene glacial events and subsequent sudden population growth from a very limited group of founders. Rapid range expansion of O. felineus through Asian and European territories after severe bottleneck points to a high dispersal potential of this trematode species.

Comparative phylogeography and demographic history of European shads (Alosa alosa and A. fallax) inferred from mitochondrial DNA

Background

Comparative broad-scale phylogeographic studies of aquatic organisms provide insights on biotic responses to the paleohydrological dynamics associated with climatic oscillations. These insights can be used to formulate a framework for understanding the evolutionary history of a species or closely related taxa as well as aid in predictive modeling of further responses to climate change. Anadromous fishes constitute interesting models for understanding the relative importance of environmental versus biological factors in shaping intraspecific genetic substructure on the interface between marine and freshwater realms. European shads, Alosa alosa and A. fallax are anadromous species that have persisted through historical large-scale environmental perturbations and now additionally face an array of anthropogenic challenges. A comprehensive phylogeographic investigation of these species is needed to provide insights on both the historical processes that have shaped their extant genetic structure and diversity, and the prospects for their future management and conservation.

Results

Despite introgressive hybridization, A. alosa and A. fallax are genetically divergent, congruent with previous studies. Three similarly divergent mtDNA clades were recognized within both A. fallax and A. alosa, most likely originating during common periods of isolation during the Pleistocene among the studied oceanographic regions. Periods of basin isolation apparently extended to the Black Sea as additional Alosa clades occur there. The present day geographic distribution of genetic diversity within European Alosa sp. suggests the existence of a strong but permeable barrier between the Atlantic and Mediterranean seas, as shown for a number of other aquatic species. Overall mtDNA diversity is considerably lower for A. alosa compared to A. fallax, suggesting that the former species is more sensitive to climatic as well as anthropogenic changes. For A. fallax, migration from the Mediterranean to the Atlantic was detected but not in the opposite direction, with colonization of the North Atlantic probably occurring after last glacial maximum.

Conclusion

The similar haplotype network topologies between the two species support a common intraspecific history of isolation. Despite these similarities, A. alosa and A. fallax have clearly responded differently to the hydrological dynamics of the Pleistocene, as reflected in their distinct demographic histories. As the species additionally occupy different ecological niches it should not be surprising that they differ in resilience to natural or human-mediated climatic changes. For A. fallax, it is further clear that its demographic response to large-scale hydrological events is not synchronized between the Atlantic and Mediterranean basins. These regional and species-specific differences should be incorporated into future predictive modeling of biological response to climate change as well as current management concepts.

Dynamic formation of asexual diploid and polyploid lineages: multilocus analysis of Cobitis reveals the mechanisms maintaining the diversity of clones

Given the hybrid genomic constitutions and increased ploidy of many asexual animals, the identification of processes governing the origin and maintenance of clonal diversity provides useful information about the evolutionary consequences of interspecific hybridization, asexuality and polyploidy. In order to understand the processes driving observed diversity of biotypes and clones in the Cobitis taenia hybrid complex, we performed fine-scale genetic analysis of Central European hybrid zone between two sexual species using microsatellite genotyping and mtDNA sequencing. We found that the hybrid zone is populated by an assemblage of clonally (gynogenetically) reproducing di-, tri- and tetraploid hybrid lineages and that successful clones, which are able of spatial expansion, recruit from two ploidy levels, i.e. diploid and triploid. We further compared the distribution of observed estimates of clonal ages to theoretical distributions simulated under various assumptions and showed that new clones are most likely continuously recruited from ancestral populations. This suggests that the clonal diversity is maintained by dynamic equilibrium between origination and extinction of clonal lineages. On the other hand, an interclonal selection is implied by nonrandom spatial distribution of individual clones with respect to the coexisting sexual species. Importantly, there was no evidence for sexually reproducing hybrids or clonally reproducing non-hybrid forms. Together with previous successful laboratory synthesis of clonal Cobitis hybrids, our data thus provide the most compelling evidence that 1) the origin of asexuality is causally linked to interspecific hybridization; 2) successful establishment of clones is not restricted to one specific ploidy level and 3) the initiation of clonality and polyploidy may be dynamic and continuous in asexual complexes.

Postglacial colonization of the Qinling Mountains: phylogeography of the swelled vent frog (Feirana quadranus)

BACKGROUND

The influence of Pleistocene climatic fluctuations on intraspecific diversification in the Qinling-Daba Mountains of East Asia remains poorly investigated. We tested hypotheses concerning refugia during the last glacial maximum (LGM) in this region by examining the phylogeography of the swelled vent frog (Feirana quadranus; Dicroglossidae, Anura, Amphibia).

METHODOLOGY/PRINCIPAL FINDINGS

We obtained complete mitochondrial ND2 gene sequences of 224 individuals from 34 populations of Feirana quadranus for phylogeographic analyses. Additionally, we obtained nuclear tyrosinase gene sequences of 68 F. quadranus, one F. kangxianensis and three F. taihangnica samples to test for mitochondrial introgression among them. Phylogenetic analyses based on all genes revealed no introgression among them. Phylogenetic analyses based on ND2 datasets revealed that F. quadranus was comprised of six lineages which were separated by deep valleys; the sole exception is that the Main Qinling and Micang-Western Qinling lineages overlap in distribution. Analyses of population structure indicated restricted gene flow among lineages. Coalescent simulations and divergence dating indicated that the basal diversification within F. quadranus may be associated with the dramatic uplifts of the Tibetan Plateau during the Pliocene. Coalescent simulations indicated that Wuling, Daba, and Western Qinling-Micang-Longmen Mountains were refugia for F. quadranus during the LGM. Demographic analyses indicated that the Daba lineage experienced population size increase prior to the LGM but the Main Qinling and the Micang-Western Qinling lineages expanded in population size and range after the LGM, and the other lineages almost have stable population size or slight slow population size decline.

CONCLUSIONS/SIGNIFICANCE

The Qinling-Daba Mountains hosted three refugia for F. quadranus during the LGM. Populations that originated in the Daba Mountains colonized the Main Qinling Mountains after the LGM. Recent sharp expansion of the Micang-Western Qinling and Main Qinling lineages probably contribute to their present-day secondary contact.

Nuclear and mitochondrial phylogeography of the European fire-bellied toads Bombina bombina and Bombina variegata supports their independent histories

Exact location and number of glacial refugia still remain unclear for many European cold-blooded terrestrial vertebrates. We performed a fine-scaled multilocus phylogeographic analysis of two Bombina species combining mitochondrial variation of 950 toads from 385 sites and nuclear genes (Rag-1, Ncx-1) from a subset of samples to reconstruct their colonization and contemporary variation patterns. We identified the lowlands northwest of the Black Sea and the Carpathians to be important refugial areas for B. bombina and B. variegata, respectively. This result emphasizes the importance of Central European refugia for ectothermic terrestrial species, far north of the Mediterranean areas regarded as exclusive glacial refugia for the animals. Additional refugia for B. variegata have been located in the southern Apennines and Balkans. In contrast, no evidence for the importance of other east European plains as refugial regions has been found. The distribution of mtDNA and Ncx-1 variation suggests the presence of local refugia near the Black Sea for B. bombina; however, coalescent simulations did not allow to distinguish whether one or two refugia were present in the region. Strong genetic drift apparently accompanied postglacial expansions reducing diversity in the colonization areas. Extended sampling, coupled with the multilocus isolation with migration analysis, revealed a limited and geographically restricted gene flow from the Balkan to Carpathian populations of B. variegata. However, despite proximity of inferred B. bombina and B. variegata refugia, gene exchange between them was not detected.

Population history, phylogeography, and conservation genetics of the last Neotropical mega-herbivore, the lowland tapir (Tapirus terrestris)

BACKGROUND

Understanding the forces that shaped Neotropical diversity is central issue to explain tropical biodiversity and inform conservation action; yet few studies have examined large, widespread species. Lowland tapir (Tapirus terrrestris, Perissodactyla, Tapiridae) is the largest Neotropical herbivore whose ancestors arrived in South America during the Great American Biotic Interchange. A Pleistocene diversification is inferred for the genus Tapirus from the fossil record, but only two species survived the Pleistocene megafauna extinction. Here, we investigate the history of lowland tapir as revealed by variation at the mitochondrial gene Cytochrome b, compare it to the fossil data, and explore mechanisms that could have shaped the observed structure of current populations.

RESULTS

Separate methodological approaches found mutually exclusive divergence times for lowland tapir, either in the late or in the early Pleistocene, although a late Pleistocene divergence is more in tune with the fossil record. Bayesian analysis favored mountain tapir (T. pinchaque) paraphyly in relation to lowland tapir over reciprocal monophyly, corroborating the inferences from the fossil data these species are sister taxa. A coalescent-based analysis rejected a null hypothesis of allopatric divergence, suggesting a complex history. Based on the geographic distribution of haplotypes we propose (i) a central role for western Amazonia in tapir diversification, with a key role of the ecological gradient along the transition between Andean subcloud forests and Amazon lowland forest, and (ii) that the Amazon river acted as an barrier to gene flow. Finally, the branching patterns and estimates based on nucleotide diversity indicate a population expansion after the Last Glacial Maximum.

CONCLUSIONS

This study is the first examining lowland tapir phylogeography. Climatic events at the end of the Pleistocene, parapatric speciation, divergence along the Andean foothill, and role of the Amazon river, have similarly shaped the history of other taxa. Nevertheless further work with additional samples and loci is needed to improve our initial assessment. From a conservation perspective, we did not find a correspondence between genetic structure in lowland tapir and ecogeographic regions proposed to define conservation priorities in the Neotropics. This discrepancy sheds doubt into this scheme's ability to generate effective conservation planning for vagile species.

Phylogenetic relationships and biogeographical patterns in Circum-Mediterranean subfamily Leuciscinae (Teleostei, Cyprinidae) inferred from both mitochondrial and nuclear data

BACKGROUND

Leuciscinae is a subfamily belonging to the Cyprinidae fish family that is widely distributed in Circum-Mediterranean region. Many efforts have been carried out to deciphering the evolutionary history of this group. Thus, different biogeographical scenarios have tried to explain the colonization of Europe and Mediterranean area by cyprinids, such as the "north dispersal" or the "Lago Mare dispersal" models. Most recently, Pleistocene glaciations influenced the distribution of leuciscins, especially in North and Central Europe. Weighing up these biogeographical scenarios, this paper constitutes not only the first attempt at deciphering the mitochondrial and nuclear relationships of Mediterranean leuciscins but also a test of biogeographical hypotheses that could have determined the current distribution of Circum-Mediterranean leuciscins.

RESULTS

A total of 4439 characters (mitochondrial + nuclear) from 321 individuals of 176 leuciscine species rendered a well-supported phylogeny, showing fourteen main lineages. Analyses of independent mitochondrial and nuclear markers supported the same main lineages, but basal relationships were not concordant. Moreover, some incongruence was found among independent mitochondrial and nuclear phylogenies. The monophyly of some poorly known genera such as Pseudophoxinus and Petroleuciscus was rejected. Representatives of both genera belong to different evolutionary lineages. Timing of cladogenetic events among the main leuciscine lineages was gained using mitochondrial and all genes data set.

CONCLUSIONS

Adaptations to a predatory lifestyle or miniaturization have superimposed the morphology of some species. These species have been separated into different genera, which are not supported by a phylogenetic framework. Such is the case of the genera Pseudophoxinus and Petroleuciscus, which real taxonomy is not well known. The diversification of leuciscine lineages has been determined by intense vicariant events following the paleoclimatological and hydrogeological history of Mediterranean region. We propose different colonization models of Mediterranean region during the early Oligocene. Later vicariance events promoted Leuciscinae diversification during Oligocene and Miocene periods. Our data corroborate the presence of leuciscins in North Africa before the Messinian salinity crisis. Indeed, Messinian period appears as a stage of gradually Leuciscinae diversification. The rise of humidity at the beginning of the Pliocene promoted the colonization and posterior isolation of newly established freshwater populations. Finally, Pleistocene glaciations determined the current European distribution of some leuciscine species.

Multiple factors have shaped the phylogeography of Chinese spiny loach Cobitis sinensis in Taiwan as inferred from mitochondrial DNA variation

Mitochondrial DNA cytochrome b sequences (1140 bp) in 61 specimens of Chinese spiny loach Cobitis sinensis from 12 drainages in Taiwan were identified as two major clades, exhibiting a southern and a northern distribution, north of TzengWen and south of TzengWen (including TzengWen), respectively. The divergence time between these two phylogroups was estimated at 7.34-9.06 million years before present (B.P.), but these two phylogroups were formed c. 3.41-4.23 and 2.22-2.75 M B.P., respectively. Moreover, geological events have been recalculated that Taiwan Island emerged above sea level at an estimate of c. 4-5 M B.P., and quickly became its present shape at c. 2 M B.P. through mountain building. These results suggest that these two major clades of C. sinensis in Taiwan might originate from two different continental populations, since the island's initial isolation in the Pliocene. Within southern Taiwan, the initial colonization was hypothesized to be in KaoPing River, followed by its northward dispersal. The high divergence between KaoPing and TzengWen was influenced by glaciations and landforms. Within north Taiwan, the colonization was from the Miaoli Plateau through western Taiwan to north-eastern and northern Taiwan. This dispersal pattern is concordant with the previously proposed hypothesis. Apparently, both geological and phylogeographic evidence suggested that river capture of the upper Takia River by the LanYang River promoted range expansion in freshwater fishes and also indicated that the Central Range within Taiwan did not act as a barrier to the dispersal of C. sinensis.

Shallow phylogeographic structuring of Vimba vimba across Europe suggests two distinct refugia during the last glaciation

Genetic variation and geographical structuring of vimba Vimba vimba were analysed across 26 sites (80 individuals) by means of mtDNA sequences (cyt b gene, mitochondrial control region) to localize hypothesized glacial refugia and to reconstruct postglacial recoloniation routes. Although genetic diversity among sequenced individuals was low, a combined analysis of the two sequenced fragments revealed a western (central and northern Europe: Danube, Elbe and lakes of Sweden) and an eastern clade (eastern Europe: Dnieper-South Bug, Don, Neman). Furthermore, a number of divergent ancestral haplotypes distributed around the Black and Caspian Seas became apparent. Mismatch analyses supported a sudden expansion model for the populations of the western clade between 50 and 10 000 bp. Overall, the study provides strong evidence for a northward and westward expansion of V. vimba from two refugial regions located in the Danubian drainage and the northern Pontic regions respectively.

Matrilinear phylogeography and demographical patterns of Rutilus rutilus: implications for taxonomy and conservation

A phylogeographic analysis of mitochondrial DNA sequence variation was carried out to infer the geographical distribution of the genealogical lineages and the historical demography of roach Rutilus rutilus (L.). A total of 265 individuals from 52 sites covering most of the Eurasian distribution range were sequenced for a 475 bp fragment of the mitochondrial cytochrome b gene. The monophyletic roach contained two deep clades that dated back to the Pliocene. The Ponto-Caspian clade comprised populations from Greece to Siberia with a likely palaeorefugium at the west coast of the Caspian Sea. This clade largely corresponds to individuals with morphological features described as Rutilus heckelii. The west European clade included individuals from central and western Europe with the Danube and Dniester basins as possible palaeorefugia. This clade largely corresponds to individuals with morphological features described as R. rutilus. A suture-zone of the two main lineages was observed along the coastal region of the Black Sea. The neutrality tests and the mismatch distributions indicated a demographic expansion during the Middle-Pleistocene for both clades.

Unexpected relationships of substructured populations in Chinese Locusta migratoria

Background

Highly migratory species are usually expected to have minimal population substructure because strong gene flow has the effect of homogenizing genetic variation over geographical populations, counteracting random drift, selection and mutation. The migratory locust Locusta migratoria belongs to a monotypic genus, and is an infamous pest insect with exceptional migratory ability – with dispersal documented over a thousand kilometers. Its distributional area is greater than that of any other locust or grasshopper, occurring in practically all the temperate and tropical regions of the eastern hemisphere. Consequently, minimal population substructuring is expected. However, in marked contrast to its high dispersal ability, three geographical subspecies have been distinguished in China, with more than nine being biologically and morphologically identified in the world. Such subspecies status has been under considerable debate.

Results

By multilocus microsatellite genotyping analysis, we provide ample genetic evidence for strong population substructure in this highly migratory insect that conforms to geography. More importantly, our genetic data identified an unexpected cryptic subdivision and demonstrated a strong affiliation of the East China locusts to those in Northwest/Northern China. The migratory locusts in China formed three distinct groups, viz. (1) the Tibetan group, comprising locusts from Tibet and nearby West China high mountain regions; this is congruent with the previously recognized Tibetan subspecies, L. m. tibetensis; (2) the South China group, containing locusts from the Hainan islands; this corresponds to the Southeast Asia oriental tropical subspecies L. m. manilensis; (3) the North China group, including locusts from the Northwest and Northern China (the Asiatic subspecies L. m. migratoria), Central China and Eastern China regions. Therefore, the traditional concept on Locusta subspecies status established from Uvarov in 1930s needs to be revised. The three groups of locusts probably have separate evolutionary histories that were most likely linked to Quaternary glaciations events, and derived from different ancestral refugial populations following postglacial expansions.

Conclusion

The migratory locust populations in China have differentiated into three genetically distinct groups despite high dispersal capability. While this clarified long-standing suspicions on the subspecific diversification of this species in China, it also revealed that the locusts in the vast area of East China are not the oriental subspecies but the Asiatic subspecies, an unexpected substructuring pattern. The distribution pattern of the three locust groups in China may be primarily defined by adaptive differentiation coupled to Quaternary glaciations events. Our results are of general significance both for locust research and for phylogeographical study of flora and fauna in China, illustrating the potential importance of phylogeographical history in shaping the divergence and distribution patterns of widespread species with strong dispersal ability.

Phylogeography of the greater horseshoe bat, Rhinolophus ferrumequinum: contrasting results from mitochondrial and microsatellite data

Phylogeographical studies are typically based on haplotype data, occasionally on nuclear markers such as microsatellites, but rarely combine both. This is unfortunate because the use of markers with contrasting modes of inheritance and rates of evolution might provide a more accurate and comprehensive understanding of a species' history. Here we present a detailed study of the phylogeography of the greater horseshoe bat, Rhinolophus ferrumequinum, using 1098 bp of the mitochondrial ND2 gene from 45 localities from across its Palaearctic range to infer population history. In addition, we re-analysed a large microsatellite data set available for this species and compared the results of both markers to infer population relationships and the historical processes influencing them. We show that mtDNA, the most popular marker in phylogeography studies, yielded a misleading result, and would have led us to conclude erroneously that a single expansion had taken place in Europe. Only by combining the mitochondrial and microsatellite data sets are we able to reconstruct the species' history and show two colonization events in Europe, one before the Last Glacial Maximum (LGM) and one after it. Combining markers also revealed the importance of Asia Minor as an ancient refugium for this species and a source population for the expansion of the greater horseshoe bat into Europe before the LGM.

Making it on their own: sperm-dependent hybrid fishes (Cobitis) switch the sexual hosts and expand beyond the ranges of their original sperm donors

Interspecific hybridization may result in asexual hybrid lineages that reproduce via parthenogenesis. Contrary to true parthenogens, sperm-dependent asexuals (gynogens and hybridogens) are restricted to the range of bisexual species, generally the parental taxa, by their need for a sperm donor. It has been documented that asexual lineages may rarely use sperm from a non-parental species or even switch a host. The available literature reports do not allow distinguishing, between whether such host switches arise by the expansion of asexuals out of their parental's range (and into that of another's) or by the local extinction of a parental population followed by a host switch. The present study combines new and previously collected data on the distribution and history of gynogenetic spined loaches (Cobitis) of hybrid origin. We identified at least three clonal lineages that have independently switched their sperm dependency to different non-parental Cobitis species, and in cases incorporated their genomes. Our current knowledge of European Cobitis species and their hybrids suggests that this pattern most probably results from the expansion of gynogenetic lineages into new areas. Such expansion was independent of the original parental species. This suggests that sperm dependence is not as restrictive to geographical expansion when compared with true parthenogenesis as previously thought.

Development of nine polymorphic microsatellite loci in the spined loach, Cobitis taenia, and cross-species amplification in the related species C. elongatoides, C. taurica and C. tanaitica

Nine polymorphic microsatellite loci were developed for the spined loach, Cobitis taenia (Teleostei: Cobitidae). The loci were validated using 50 individuals from a population in Belgium. Moderate to high levels of polymorphism were detected (two to 11 alleles). In addition, most markers amplified successfully in three closely related taxa that are known to hybridize with C. taenia: C. elongatoides, C. taurica and C. tanaitica. Some of the loci are most likely diagnostic among species. These markers will be valuable for the study of the historical and contemporary interactions within C. taenia and the Cobitis species complex.

Ancient divisions, recent expansions: phylogeography and population genetics of the round goby Apollonia melanostoma

During the past two decades, the round goby Apollonia melanostoma (=Neogobius melanostomus) has expanded its range via shipping transport and canals, extending north and west from the Ponto-Caspian region of Eurasia and to the North American Great Lakes. Exotic populations of the round goby have been very successful in the Baltic Sea and the Great Lakes regions, exerting significant ecological changes. Our study evaluates the population genetic and biogeographical structure of the round goby across its native and nonindigenous ranges, in light of geological history and its expansion pathways. We analyzed seven new nuclear microsatellite loci and mitochondrial DNA cytochrome b gene sequences from 432 individuals in 22 locations. Population structure was tested using F(ST)-analogs, phylogenetic trees, clustering diagrams, Bayesian assignment tests and nested clade analyses. Results show that native populations in the Black vs. the Caspian Sea basins diverge by 1.4% and c. 350,000 years, corresponding to closure of their prior connections and supporting the taxonomic separation of the Black Sea A. m. melanostoma from the Caspian Sea A. m. affinis. Their within-basin populations diverge by approximately 0.4% and 100,000 years. Nonindigenous populations in the Baltic Sea and Danube and Dnieper Rivers trace to separate northern Black Sea origins, whereas the upper Volga River system houses mixed populations of A. m. melanostoma and A. m. affinis. Native populations average twice the genetic diversity of most exotic sites; however, sites in the Volga River system have high diversity due to mixing of the two taxa. Our results highlight how vicariance and anthropogenic disturbances have shaped a rapidly expanding species' genetic heritage.

Longitudinal range expansion and cryptic eastern species in the western Palaearctic oak gallwasp, Andricus coriarius

The oak gallwasp Andricus coriarius is distributed across the Western Palaearctic from Morocco to Iran. It belongs to a clade of host-alternating Andricus species that requires host oaks in two sections of Quercus subgenus Quercus to complete its lifecycle, a requirement that has restricted the historic distribution and dispersal of members of this clade. Here we present nuclear and mitochondrial sequence evidence from the entire geographic range of A. coriarius to investigate the genetic legacy of longitudinal range expansion. We show A. coriarius as currently understood to be para- or polyphyletic, with three evolutionarily independent (but partially sympatric) lineages that diverged c. 10 million years ago (mya). The similarities in gall structure that have justified recognition of single species to date thus represent either strong conservation of an ancestral state or striking convergence. All three lineages originated in areas to the east of Europe, underlining the significance of Turkey, Iran and the Levant as 'cradles' of gallwasp evolution. One of the three lineages gave rise to all European populations, and range expansion from a putative Eastern origin to the present distribution is predicted to have occurred around 1.6 mya.

A phylogeographic analysis of populations of the wood turtle (Glyptemys insculpta) throughout its range

We investigated the phylogeography of wood turtles, Glyptemys insculpta, in North America using 750 bp of the mitochondrial control region from 117 individuals sampled at 29 localities across the species' range. A total of 21 haplotypes were identified and little genetic variation was found. The highest pairwise difference was 2%. From nested clade analysis (NCA), one main postglacial dispersal route was inferred along the east coast, with subsequent westward dispersal. NCA further revealed some patterns of restricted gene flow/dispersal. We propose that wood turtles experienced a combined effect of bottleneck during the Pleistocene as well as of selective sweep, which produced the low level of variation observed. Following the Pleistocene, wood turtles would have undergone a rapid northward expansion from a common southern refugium as glaciers retreated. These findings shed light on where to direct conservation priorities, on conservation strategies needed, and on the potential effects of interpopulation transfers for this vulnerable species.

Molecular evidence of Pleistocene bidirectional faunal exchange between Europe and the Near East: the case of the bicoloured shrew (Crocidura leucodon, Soricidae)

We sequenced 1077 bp of the mitochondrial cytochrome b gene and 511 bp of the nuclear Apolipoprotein B gene in bicoloured shrew (Crocidura leucodon, Soricidae) populations ranging from France to Georgia. The aims of the study were to identify the main genetic clades within this species and the influence of Pleistocene climatic variations on the respective clades. The mitochondrial analyses revealed a European clade distributed from France eastwards to north-western Turkey and a Near East clade distributed from Georgia to Romania; the two clades separated during the Middle Pleistocene. We clearly identified a population expansion after a bottleneck for the European clade based on mitochondrial and nuclear sequencing data; this expansion was not observed for the eastern clade. We hypothesize that the western population was confined to a small Italo-Balkanic refugium, whereas the eastern population subsisted in several refugia along the southern coast of the Black Sea.

Geological Dates and Molecular Rates: Rapid Divergence of Rivers and Their Biotas

We highlight a novel molecular clock calibration system based on geologically dated river reversal and river capture events. Changes in drainage pattern may effect vicariant isolation of freshwater taxa, and thus provide a predictive framework for associated phylogeographic study. As a case in point, New Zealand's Pelorus and Kaituna rivers became geologically isolated from the larger Wairau River system 70 to 130 kyr BP. We conducted mitochondrial DNA phylogeographic analyses of two unrelated freshwater-limited fish taxa native to these river systems (Gobiomorphus breviceps, n = 63; Galaxias divergens, n = 95). Phylogenetic analysis of combined control region and cytochrome b sequences yielded reciprocally monophyletic clades of Pelorus-Kaituna and Wairau haplotypes for each species. Calibrated rates of molecular change based on this freshwater vicariant event are substantially faster than traditionally accepted rates for fishes but consistent with other recent inferences based on geologically young calibration points. A survey of freshwater phylogeographic literature reveals numerous examples in which the ages of recent evolutionary events may have been substantially overestimated through the use of "accepted" calibrations. We recommend that--wherever possible--biologists should start to reassess the conclusions of such studies by using more appropriate molecular calibrations derived from recent geological events.