Evolutionary melting pots: a biodiversity hotspot shaped by ring diversifications around the Black Sea in the Eastern tree frog (Hyla orientalis)

The Extraordinary Diversity of Merodon avidus Complex (Diptera: Syrphidae)-Adding New Areas, New Species and a New Molecular Marker

In this paper, the Merodon avidus (Diptera, Syrphidae) species complex was revised, whereupon we discovered and described four new species for science: Merodon atroavidus Vujić, Radenković et Likov sp. nov., M. magnus Vujić, Kočiš Tubić et Ačanski sp. nov., M. nigroscutum Vujić, Radenković et Likov sp. nov. and M. pseudomoenium Vujić, Kočiš Tubić et Ačanski sp. nov. An integrative taxonomy approach was used to delimit species boundaries. Two molecular markers (the mitochondrial COI gene and nuclear 28S rRNA gene-newly analysed marker for the complex) and geometric morphometry of the wing shape, together with morphological data and distribution, successfully separated all species from the complex. The morphological variability of the analysed species is described and discussed and an illustrated diagnostic key for typical morpho-forms of species from the M. avidus complex is presented. A distribution map of all investigated species from the complex is provided. The level of endemicity of the M. avidus complex was discussed.

Pericentromeric recombination suppression and the 'large X effect' in plants

X chromosome was reported to be a major contributor to isolation between closely related species-the 'large X' effect (LXE). The causes of LXE are not clear, but the leading theory is that it is caused by recessive species incompatibilities exposed in the phenotype due to the hemizygosity of X-linked genes in the heterogametic sex. However, the LXE was also reported in species with relatively recently evolved sex chromosomes where Y chromosome is not completely degenerate and X-linked genes are not hemizygous, such as the plant Silene latifolia. Recent genome sequencing and detailed genetic mapping in this species revealed a massive (> 330 Mb) non- or rarely-recombining pericentromeric region on the X chromosome (Xpr) that comprises ~ 90% of the chromosome and over 13% of the entire genome. If any of the Xpr genes are involved in species incompatibilities, this would oppose interspecific gene flow for other genes tightly linked in the Xpr. Here we test the hypothesis that the previously reported LXE in S. latifolia is caused by the lack of recombination on most of the X chromosome. Based on genome-wide analysis of DNA polymorphism and gene expression in S. latifolia and its close cross-compatible relative S. dioica, we report that the rarely-recombining regions represent a significant barrier for interspecific gene flow. We found little evidence for any additional factors contributing to the LXE, suggesting that extensive pericentromeric recombination suppression on the X-chromosome is the major if not the only cause of the LXE in S. latifolia and S. dioica.

The ecological importance of hybridization

Hybridization as an evolutionary process has been studied in depth over the past few decades. Research has focused on its role in shaping reproductive barriers, its adaptive value, and its genomic consequences. In contrast, our knowledge of ecological dimensions of hybridization is still in its infancy, despite hybridization being an inherently ecological interaction. Using examples from various organisms, we show that hybridization can affect and be affected by non-reproductive interactions, including predation, competition, parasitism, mutualism, commensalism, and organism-environment interactions, with significant implications for community structure and ecosystem functioning. However, since these dimensions of hybridization have mostly been revealed from studies designed to decipher other evolutionary processes, we argue that much of the eco-evolutionary importance of hybridization is yet to be discovered.

Unrecognized diversity, genetic structuring, and phylogeography of the genus Triplophysa (Cypriniformes: Nemacheilidae) sheds light on two opposite colonization routes during Quaternary glaciation that occurred in the Qilian Mountains

In recent years, the influence of historical geological and climatic events on the evolution of flora and fauna in the Tibetan Plateau has been a hot research topic. The Qilian Mountain region is one of the most important sources of biodiversity on the Qinghai-Tibet Plateau. Many species existed in the region during the Pleistocene glacial oscillation, and the complex geographical environment provided suitable conditions for the survival of local species. The shrinkage, expansion, and transfer of the distribution range and population size of species have significant effects on genetic diversity and intraspecific differentiation. To reveal the effects of geological uplift and climate oscillation on the evolution of fish populations in the Qilian Mountains, we investigated the genetic structure, phylogenetic relationship, and phylogeographical characteristics of genus Triplophysa species in the Qilian Mountains using the mitochondrial DNA gene (COI), three nuclear genes (RAG1, sRH, and Myh6) and 11 pairs of nuclear microsatellite markers. We collected 11 species of genus Triplophysa living in the Qilian Mountains, among which Triplophysa hsutschouensis and Triplophysa papillosolabiata are widely distributed in the rivers on the northern slope of the Qilian Mountains. There was a high degree of lineage differentiation among species, and the genetic diversity of endemic species was low. The different geographical groups of T. papillosolabiata presented some allogeneic adaptation and differentiation, which was closely related to the changes in the river system. Except for the population expansion event of T. hsutschouensis during the last glacial period of the Qinghai-Tibet Plateau (0.025 MYA), the population sizes of other plateau loach species remained stable without significant population expansion. Starting from the east and west sides of the Qilian Mountains, T. hsutschouensis, and T. papillosolabiata showed two species colonization routes in opposite directions. The geological events of the uplift of the Qinghai-Tibet Plateau and the climatic oscillation of the Quaternary glaciation had a great influence on the genetic structure of the plateau loach in the Qilian Mountains, which promoted the genetic differentiation of the plateau loach and formed some unique new species. The results of this study have important guiding significance for fish habitat protection in the Qilian Mountains.

Wind-dispersed seeds blur phylogeographic breaks: The complex evolutionary history of Populus lasiocarpa around the Sichuan Basin

The strength of phylogeographic breaks can vary among species in the same area despite being subject to the same geological and climate history due to differences in biological traits. Several important phylogeographic breaks exist around the Sichuan Basin in Southwest China but few studies have focused on wind-dispersed plants. Here, we investigated the phylogeographic patterns and the evolutionary history of Populus lasiocarpa, a wind-pollinated and wind-dispersed tree species with a circum-Sichuan Basin distribution in southwest China. We sequenced and analyzed three plastid DNA fragments (ptDNA) and eight nuclear microsatellites (nSSRs) of 265 individuals of P. lasiocarpa from 21 populations spanning the entire distribution range. Distribution patterns based on nSSR data revealed that there are three genetic groups in P. lasiocarpa. This is consistent with the three phylogeographic breaks (Sichuan Basin, the Kaiyong Line and the 105°E line), where the Sichuan basin acts as the main barrier to gene flow between western and eastern groups. However, the distribution pattern based on ptDNA haplotypes poorly matched the phylogeographic breaks, and wind-dispersed seeds may be one of the main contributing factors. Species distribution modelling suggested a larger potential distribution in the last glacial maximum with a severe bottleneck during the last interglacial. A DIYABC model also suggested a population contraction and expansion for both western and eastern lineages. These results indicate that biological traits are likely to affect the evolutionary history of plants, and that nuclear molecular markers, which experience higher levels of gene flow, might be better indicators of phylogeographic breaks.

New insights on patterns of genetic admixture and phylogeographic history in Iberian high mountain populations of midwife toads

Multiple Quaternary glacial refugia in the Iberian Peninsula, commonly known as "refugia within refugia", allowed diverging populations to come into contact and admix, potentially boosting substantial mito-nuclear discordances. In this study, we employ a comprehensive set of mitochondrial and nuclear markers to shed light onto the drivers of geographical differentiation in Iberian high mountain populations of the midwife toads Alytes obstetricans and A. almogavarii from the Pyrenees, Picos de Europa and Guadarrama Mountains. In the three analysed mountain regions, we detected evidence of extensive mito-nuclear discordances and/or admixture between taxa. Clustering analyses identified three major divergent lineages in the Pyrenees (corresponding to the eastern, central and central-western Pyrenees), which possibly recurrently expanded and admixed during the succession of glacial-interglacial periods that characterised the Late Pleistocene, and that currently follow a ring-shaped diversification pattern. On the other hand, populations from the Picos de Europa mountains (NW Iberian Peninsula) showed a mitochondrial affinity to central-western Pyrenean populations and a nuclear affinity to populations from the central Iberian Peninsula, suggesting a likely admixed origin for Picos de Europa populations. Finally, populations from the Guadarrama Mountain Range (central Iberian Peninsula) were depleted of genetic diversity, possibly as a consequence of a recent epidemic of chytridiomycosis. This work highlights the complex evolutionary history that shaped the current genetic composition of high mountain populations, and underscores the importance of using a multilocus approach to better infer the dynamics of population divergence.

The introduction of three cryptic tree frog species in the Dutch coastal dunes challenges conservation paradigms

Invasive species are considered one of the main drivers of the sixth mass extinction. Conservation solutions depend on whether a species is also indigenous to the country it invades (i.e., beyond its native range). In the case of invasive cryptic species, genetic tools are required to establish their identity. We illustrate these issues with the human-mediated colonization of the Dutch coastal dunes by Hyla tree frogs. Although previously assumed to concern the indigenous common tree frog H. arborea, European tree frogs comprise a complex of allopatric cryptic species, meaning the taxonomic identity of introduced Dutch populations warrants investigation. We sequence mtDNA for 173 individuals from native and introduced populations across the Netherlands and compare our dataset with hundreds of Hyla haplotypes previously barcoded in the Western Palearctic. Two of the dune populations carry an mtDNA haplotype of the native species H. arborea that occurs naturally elsewhere in the Netherlands. In contrast, mtDNA assigned to the eastern tree frog H. orientalis was detected in all three other dune populations. In one of these populations mtDNA of the Italian tree frog H. intermedia was also found. Not one, but three species of tree frogs have thus been introduced to the Dutch coastal dunes, only one of them being native to the Netherlands. This situation causes a conservation conundrum as some introduced populations are lawfully protected but could pose a threat to local biodiversity. Regarding the ‘true’ exotic tree frog species, all conservation options should be considered.

Unusual evolution of tree frog populations in the Chernobyl exclusion zone

Despite the ubiquity of pollutants in the environment, their long‐term ecological consequences are not always clear and still poorly studied. This is the case concerning the radioactive contamination of the environment following the major nuclear accident at the Chernobyl nuclear power plant. Notwithstanding the implications of evolutionary processes on the population status, few studies concern the evolution of organisms chronically exposed to ionizing radiation in the Chernobyl exclusion zone. Here, we examined genetic markers for 19 populations of Eastern tree frog (Hyla orientalis) sampled in the Chernobyl region about thirty years after the nuclear power plant accident to investigate microevolutionary processes ongoing in local populations. Genetic diversity estimated from nuclear and mitochondrial markers showed an absence of genetic erosion and higher mitochondrial diversity in tree frogs from the Chernobyl exclusion zone compared to other European populations. Moreover, the study of haplotype network permitted us to decipher the presence of an independent recent evolutionary history of Chernobyl exclusion zone's Eastern tree frogs caused by an elevated mutation rate compared to other European populations. By fitting to our data a model of haplotype network evolution, we suspected that Eastern tree frog populations in the Chernobyl exclusion zone have a high mitochondrial mutation rate and small effective population sizes. These data suggest that Eastern tree frog populations might offset the impact of deleterious mutations because of their large clutch size, but also question the long‐term impact of ionizing radiation on the status of other species living in the Chernobyl exclusion zone.

Morphometry of two cryptic tree frog species at their hybrid zone reveals neither intermediate nor transgressive morphotypes

Under incomplete reproductive isolation, secondary contact of diverged allopatric lineages may lead to the formation of hybrid zones that allow to study recombinants over several generations as excellent systems of genomic interactions resulting from the evolutionary forces acting on certain genes and phenotypes. Hybrid phenotypes may either exhibit intermediacy or, alternatively, transgressive traits, which exceed the extremes of their parents due to epistasis and segregation of complementary alleles. While transgressive morphotypes have been examined in fish, reptiles, birds, and mammals, studies in amphibians are rare. Here, we associate microsatellite-based genotypes with morphometrics-based morphotypes of two tree frog species of the Hyla arborea group, sampled across a hybrid zone in Poland, to understand whether the genetically differentiated parental species also differ in morphology between each other and their hybrids and whether secondary contact leads to the evolution of intermediate or transgressive morphotypes. Using univariate approaches, explorative multivariate methods (principal component analyses) as well as techniques with prior grouping (discriminant function analyses), we find that morphotypes of both parental species and hybrids differ from each other. Importantly, hybrid morphotypes are neither intermediate nor transgressive but found to be more similar to H. orientalis than to H. arborea.

Morphological diversification of Mediterranean anurans: the roles of evolutionary history and climate

Investigation of the ecological and evolutionary mechanisms governing the origin and diversification of species requires integrative approaches that often have to accommodate strong discordance among datasets. A common source of conflict is the combination of morphological and molecular characters with different evolutionary rates. Resolution of these discordances is crucial to assess the relative roles of different processes in generating and maintaining biodiversity. Anuran amphibians provide many examples of morphologically similar, genetically divergent lineages, posing questions about the relative roles of phylogeny and ecological factors in phenotypic evolution. We focused on three circum-Mediterranean anuran genera (Hyla, Alytes and Discoglossus), characterizing morphological and environmental disparity and comparing diversity patterns across biological levels of organization. Using a comparative phylogenetic framework, we tested how shared ancestry and climatic factors come together to shape phenotypic diversity. We found higher morphological differentiation within Hyla and Alytes than in Discoglossus. Body size and limb morphology contributed most to inter- and intraspecific morphological variation in Hyla and Alytes, but there was no strong phylogenetic signal, indicating that shared ancestry does not predict patterns of phenotypic divergence. In contrast, we uncovered a significant association between morphology and climatic descriptors, supporting the hypothesis that morphological disparity between species results from adaptive evolution.

Diversity, distribution and molecular species delimitation in frogs and toads from the Eastern Palaearctic

Biodiversity analyses can greatly benefit from coherent species delimitation schemes and up-to-date distribution data. In this article, we have made the daring attempt to delimit and map described and undescribed lineages of anuran amphibians in the Eastern Palaearctic (EP) region in its broad sense. Through a literature review, we have evaluated the species status considering reproductive isolation and genetic divergence, combined with an extensive occurrence dataset (nearly 85k localities). Altogether 274 native species from 46 genera and ten families were retrieved, plus eight additional species introduced from other realms. Independent hotspots of species richness were concentrated in southern Tibet (Medog County), the circum-Sichuan Basin region, Taiwan, the Korean Peninsula and the main Japanese islands. Phylogeographic breaks responsible for recent in situ speciation events were shared around the Sichuan Mountains, across Honshu and between the Ryukyu Island groups, but not across shallow water bodies like the Yellow Sea and the Taiwan Strait. Anuran compositions suggested to restrict the zoogeographical limits of the EP to East Asia. In a rapidly evolving field, our study provides a checkpoint to appreciate patterns of species diversity in the EP under a single, spatially explicit, species delimitation framework that integrates phylogeographic data in taxonomic research.

Lack of impact of radiation on blood physiology biomarkers of Chernobyl tree frogs

BACKGROUND

Human actions have altered natural ecosystems worldwide. Among the many pollutants released to the environment, ionizing radiation can cause severe damage at different molecular and functional levels. The accident in the Chernobyl Nuclear Power Plant (1986) caused the largest release of ionizing radiation to the environment in human history. Here, we examined the impact of the current exposure to ionizing radiation on blood physiology biomarkers of adult males of the Eastern tree frog (Hyla orientalis) inhabiting within and outside the Chernobyl Exclusion Zone. We measured the levels of eight blood parameters (sodium, potassium, chloride, ionized calcium, total carbon dioxide, glucose, urea nitrogen, and anion gap), physiological markers of homeostasis, as well as of liver and kidney function.

RESULTS

Levels of blood physiology biomarkers did not vary in function of the current exposure of tree frogs to ionizing radiation within the Chernobyl Exclusion Zone. Physiological blood levels were similar in frogs inhabiting Chernobyl (both in areas with medium-high or low radiation) than in tree frogs living outside Chernobyl exposed only to background radiation levels.

CONCLUSIONS

The observed lack of effects of current radiation levels on blood biomarkers can be a consequence of the low levels of radiation currently experienced by Chernobyl tree frogs, but also to the fact that our sampling was restricted to active breeding males, i.e. potentially healthy adult individuals. Despite the clear absence of effects of current radiation levels on physiological blood parameters in tree frogs, more research covering different life stages and ecological scenarios is still needed to clarify the impact of ionizing radiation on the physiology, ecology, and dynamics of wildlife inhabiting radioactive-contaminated areas.

Pleistocene climate fluctuations as the major driver of genetic diversity and distribution patterns of the Caspian green lizard, Lacerta strigata Eichwald, 1831

Green lizards of the genus Lacerta have served as excellent models for studying the impact of Pleistocene climatic oscillations on genetic structures. The Caspian green lizard, Lacerta strigata, occupies various habitats across the Caucasus and the South Caspian Sea, with the Hyrcanian Forests and north of the Alborz Mountains forming the core of the range. This study aimed to re-examine the phylogenetic relationships of L. strigata with other congeneric members and to assess the genetic structure and historical demography of the species. Furthermore, Species Distribution Models (SDMs) were performed to infer the species' potential habitat suitability and were then projected on climate scenarios reflecting current and past (6 ky and 21 ky before present) conditions. A total of 39 individuals collected from most of the distribution range, together with additional lacertid species sequences from the GenBank database, were examined using mtDNA (Cyt b and 12S ribosomal RNA) and nuclear (C-mos and β-fibrinogen) sequence data. Based on the phylogenetic analyses, L. strigata was found to be a sister taxon to all other members of the genus. The species included two main clades (regional western and eastern) that diverged in a period between the Early and Middle Pleistocene. Based on the BBM and S-Diva analyses, both dispersal and vicariance events explained the phylogeographic structure of the species in the Hyrcanian Forests. The historical demographic analyses using Bayesian skyline plots showed a mild increase in the effective population size from about 120 Kya for the western regional clade. According to phylogeographic structures and SDMs evidence, as in other species within the region, it appears that the south of the Caspian Sea (Hyrcanian Forests), and the Alborz Mountains acted as multiple refugia during cold periods and promoted expansion outwards amid the warm periods. Overall, the results provided evidence that the genetic structure of the species has been influenced by the Pleistocene climatic fluctuations.

Ring distribution patterns-diversification or speciation? Comparative phylogeography of two small mammals in the mountains surrounding the Sichuan Basin

Studying the genetic differentiation in a unique geographical area contributes to understanding the process of speciation. Here, we explore the spatial genetic structure and underlying formation mechanism of two congeneric small mammal species (Apodemus draco and A. chevrieri), which are mainly distributed in the mountains surrounding the lowland Sichuan Basin, southwest China. We applied a set of comparative phylogeographical analyses to determine their genetic diversification patterns, combining mitochondrial (Cytb and COI) and nuclear (microsatellite loci) markers, with dense sampling throughout the range (411 A. draco from 21 sites and 191 A. chevrieri from 22 sites). Moreover, we performed three complementary statistical methods to investigate the correlation between genotype and geographical and environmental components, and predicted the potential suitable distributional range under the present and historical climate conditions. Our results suggest that both species have experienced allopatric differentiation and admixture in historical periods, resulting in a ring-shape diversification, under the barrier effect of the Sichuan Basin. We infer that the tectonic events of the Qinghai-Tibetan Plateau and climatic oscillations during the Quaternary played an important role on the genetic divergence of the two species by providing environmental heterogeneity and geographical variation. Our study reveals a case of two sympatric small mammals following a ring-shaped diversification pattern and provides insight into the process of differentiation.

Reconstructing hotspots of genetic diversity from glacial refugia and subsequent dispersal in Italian common toads (Bufo bufo)

Genetic diversity feeds the evolutionary process and allows populations to adapt to environmental changes. However, we still lack a thorough understanding of why hotspots of genetic diversity are so 'hot'. Here, we analysed the relative contribution of bioclimatic stability and genetic admixture between divergent lineages in shaping spatial patterns of genetic diversity in the common toad Bufo bufo along the Italian peninsula. We combined population genetic, phylogeographic and species distribution modelling (SDM) approaches to map ancestral areas, glacial refugia, and secondary contact zones. We consistently identified three phylogeographic lineages, distributed in northern, central and southern Italy. These lineages expanded from their ancestral areas and established secondary contact zones, before the last interglacial. SDM identified widespread glacial refugia in peninsular Italy, sometimes located under the present-day sea-level. Generalized linear models indicated genetic admixture as the only significant predictor of the levels of population genetic diversity. Our results show that glacial refugia contributed to preserving both levels and patterns of genetic diversity across glacial-interglacial cycles, but not to their formation, and highlight a general principle emerging in Mediterranean species: higher levels of genetic diversity mark populations with substantial contributions from multiple genetic lineages, irrespective of the location of glacial refugia.

A reassessment of the diversity of green toads (Bufotes) in the circum-Aegean region

The dynamic biogeography of glacial refugia may cause complex patterns of genetic admixture between parapatric taxa, which in turn can mislead their systematics, diversity, and distributions. We investigated this issue for green toads (Bufotes) inhabiting the circum-Aegean region, a biodiversity hotspot of the Eastern Mediterranean. A previous phylogeographic study based on mitochondrial and microsatellite loci identified the hybrid zone between the European (viridis) and Anatolian (sitibundus) lineages of B. viridis all over the Balkan Peninsula, but subsequent range-wide genomic analyses (>1000 SNPs) located this transition in Turkey, a thousand kilometers eastwards. In order to clarify the diversity and taxonomy of the circum-Aegean populations, we reconciled these conflicting findings by integrating previous data with pure sitibundus individuals. Our results confirmed that the viridis/sitibundus hybrid zone extends in Western Anatolia, but that southeastern European populations feature cytonuclear discordances and a high and structured microsatellite diversity. This remarkable situation may stem from a massive geographic displacement of the hybrid zone during the last glacial fluctuations, an underappreciated yet seemingly common feature among the herpetofauna of the region. Our study thus contributes to the rising view that mitochondrial DNA can be a poor predictor of current phylogeographic structure, hence the need for genomic data, especially for narrowly distributed taxa. Finally, the analyses unambiguously support the distinction of a micro-endemic clade of green toads unique to some Cyclades islands, for which we provide a formal taxonomic description.

Extensive cytonuclear discordance in a crested newt from the Balkan Peninsula glacial refugium

Integration of multilocus data and species distribution modelling into phylogeography allows mitochondrial DNA (mtDNA)-based scenarios to be fine-tuned. We address the question of whether extensive mtDNA substructuring in the crested newt Triturus macedonicus from the Balkan Peninsula is matched in the nuclear genome. We determine the intraspecific population structure based on 52 nuclear DNA markers and project a species distribution model on climate layers for the Last Glacial Maximum. We show that T. macedonicus accumulated nuclear DNA population structure in an area predicted to have been climatically stable during the Pleistocene, with four nuclear DNA groups in the western part of the species range. The distribution of these nuclear DNA groups shows little agreement with that of mtDNA structuring, which shows three highly distinct species-specific clades and a fourth one introgressed from another crested newt species. This cytonuclear discordance conveys that historical biogeographical scenarios based on mtDNA exclusively should be interpreted with caution. Our findings further highlight the important role the Balkan Peninsula has played in the evolution and preservation of European biodiversity.

The effect of phylogeographic history on species boundaries: a comparative framework in Hyla tree frogs

Because it is indicative of reproductive isolation, the amount of genetic introgression across secondary contact zones is increasingly considered in species delimitation. However, patterns of admixture at range margins can be skewed by the regional dynamics of hybrid zones. In this context, we posit an important role for phylogeographic history: hybrid zones located within glacial refugia (putatively formed during the Late-Pleistocene) should be better defined than those located in post-glacial or introduced ranges (putatively formed during the Holocene and the Anthropocene). We test this hypothesis in a speciation continuum of tree frogs from the Western Palearctic (Hyla), featuring ten identified contacts between species spanning Plio-Pleistocene to Miocene divergences. We review the rich phylogeographic literature of this group and examine the overlooked transition between H. arborea and H. molleri in Western France using a multilocus dataset. Our comparative analysis supports a trend that contacts zones resulting from post-glacial expansions and human translocations feature more extensive introgression than those established within refugial areas. Integrating the biogeographic history of incipient species, i.e. their age since first contact together with their genetic divergence, thus appears timely to draw sound evolutionary and taxonomic inferences from patterns of introgression across hybrid zones.

Are glacial refugia hotspots of speciation and cytonuclear discordances? Answers from the genomic phylogeography of Spanish common frogs

Subdivided Pleistocene glacial refugia, best known as "refugia within refugia", provided opportunities for diverging populations to evolve into incipient species and/or to hybridize and merge following range shifts tracking the climatic fluctuations, potentially promoting extensive cytonuclear discordances and "ghost" mtDNA lineages. Here, we tested which of these opposing evolutionary outcomes prevails in northern Iberian areas hosting multiple historical refugia of common frogs (Rana cf. temporaria), based on a genomic phylogeography approach (mtDNA barcoding and RAD-sequencing). We found evidence for both incipient speciation events and massive cytonuclear discordances. On the one hand, populations from northwestern Spain (Galicia and Asturias, assigned to the regional endemic R. parvipalmata), are deeply-diverged at mitochondrial and nuclear genomes (~4 My of independent evolution), and barely admix with northeastern populations (assigned to R. temporaria sensu stricto) across a narrow hybrid zone (~25 km) located in the Cantabrian Mountains, suggesting that they represent distinct species. On the other hand, the most divergent mtDNA clade, widespread in Cantabria and the Basque country, shares its nuclear genome with other R. temporaria s. s. lineages. Patterns of population expansions and isolation-by-distance among these populations are consistent with past mitochondrial capture and/or drift in generating and maintaining this ghost mitochondrial lineage. This remarkable case study emphasizes the complex evolutionary history that shaped the present genetic diversity of refugial populations, and stresses the need to revisit their phylogeography by genomic approaches, in order to make informed taxonomic inferences.

Phylogeographic analyses point to long-term survival on the spot in micro-endemic Lycian salamanders

Lycian salamanders (genus Lyciasalamandra) constitute an exceptional case of micro-endemism of an amphibian species on the Asian Minor mainland. These viviparous salamanders are confined to karstic limestone formations along the southern Anatolian coast and some islands. We here study the genetic differentiation within and among 118 populations of all seven Lyciasalamandra species across the entire genus’ distribution. Based on circa 900 base pairs of fragments of the mitochondrial 16SrDNA and ATPase genes, we analysed the spatial haplotype distribution as well as the genetic structure and demographic history of populations. We used 253 geo-referenced populations and CHELSA climate data to infer species distribution models which we projected on climatic conditions of the Last Glacial Maximum (LGM). Within all but one species, distinct phyloclades were identified, which only in parts matched current taxonomy. Most haplotypes (78%) were private to single populations. Sometimes population genetic parameters showed contradicting results, although in several cases they indicated recent population expansion of phyloclades. Climatic suitability of localities currently inhabited by salamanders was significantly lower during the LGM compared to recent climate. All data indicated a strong degree of isolation among Lyciasalamandra populations, even within phyloclades. Given the sometimes high degree of haplotype differentiation between adjacent populations, they must have survived periods of deteriorated climates during the Quaternary on the spot. However, the alternative explanation of male biased dispersal combined with a pronounced female philopatry can only be excluded if independent nuclear data confirm this result.

A tree of tree frogs around the Black Sea

Speciation, the process by which one species evolves into two or more, is a major focus of ongoing debate, particularly regarding the geographic context in which it occurs. Geographic models of speciation tend to fall into discrete categories, typically referred to as allopatric, parapatric and sympatric speciation, according to whether two groups evolve reproductive isolation while geographically isolated, differentiated but connected by gene flow, or completely co-occurring. Yet molecular studies indicate that full development of reproductive isolation can take very long compared with the timescale at which climatic oscillations occur, such that the geographic context of differentiating forms might change often during the long process to full species. Studies of genetic relationships across the ranges of organisms with low-dispersal distances have the potential to reveal these complex histories. In a particularly elegant example in this issue, Dufresnes et al. () use genetic variation and ecological niche modelling to show that a ring of populations of the eastern tree frog (Hyla orientalis) surrounding the Black Sea had a complex history of geographic differentiation. Alternating phases of geographic fragmentation and phases of gene flow between neighbouring populations have produced a pattern of gradual genetic change connecting the western, southern and eastern sides of the ring, with the northwestern and northeastern forms being most differentiated. In the north, a population in Crimea appears to have been produced through mixture of the two extreme forms. The overall genetic relationships are reminiscent of those found in ring species, which have been used as prime demonstrations of the process of speciation. The difference, however, is that the terminal forms appear to have mixed rather than be reproductively isolated, although more research is needed to infer whether there might be some reproductive isolation on the northern side of the ring.

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").

Mitochondrial phylogeography of the Iberian endemic frog Rana iberica, with implications for its conservation

Genetic characterization of species using phylogeographic approaches represents a basic reference to understand their evolutionary history as well as to identify conservation priorities to protect areas of particular interest regarding evolutionary potential. Even in well-studied regions such information is lacking for the majority of species, including many endemic species with reduced distribution ranges. We investigate the phylogeographic pattern of the Iberian frog Rana iberica, an endemic amphibian restricted to Central and North-Western Iberian Peninsula. Using mitochondrial sequences, we reconstruct the phylogeographic history of the species to test the effect of Quaternary climate changes on the evolutionary diversification of lineages, that is, the differentiation of mitochondrial lineages and the formation of genetic diversity melting pots, and integrate phylogeographic evidence for future conservation planning. Our results indicate the existence of 3 main mitochondrial lineages differentiated during the Upper Pleistocene. Both historical demographic analyses and climatic niche modeling show a strong effect of glacial climate changes, suggesting recurrent range contractions and expansions. Under such circumstances, differentiation took place most likely by isolation in allopatric interglacial refugia. Secondary lineage admixture in northern Portugal generated a broad mixed zone with highest nucleotide diversity. Given its particular evolutionary potential, its reduced distribution and eventual threats under current climate change scenario, conservation priorities should focus on the isolated lineage from Sierra de Guadalupe.

Shifting barriers and phenotypic diversification by hybridisation

The establishment of hybrid taxa relies on reproductive isolation from the parental forms, typically achieved by ecological differentiation. Here, we present an alternative mechanism, in which shifts in the strength and location of dispersal barriers facilitate diversification by hybridisation. Our case study concerns the highly diverse, stenotopic rock-dwelling cichlids of the African Great Lakes, many of which display geographic colour pattern variation. The littoral habitat of these fish has repeatedly been restructured in the course of ancient lake level fluctuations. Genetic data and an experimental cross support the hybrid origin of a distinct yellow-coloured variant of Tropheus moorii from ancient admixture between two allopatric, red and bluish variants. Deficient assortative mating preferences imply that reproductive isolation continues to be contingent on geographic separation. Linking paleolimnological data with the establishment of the hybrid variant, we sketch a selectively neutral diversification process governed solely by rearrangements of dispersal barriers.

One species or four? Yes!...and, no. Or, arbitrary assignment of lineages to species obscures the diversification processes of Neotropical fishes

Species are fundamental units in many biological disciplines, but there is continuing disagreement as to what species are, how to define them, and even whether the concept is useful. While some of this debate can be attributed to inadequate data and insufficient statistical frameworks in alpha taxonomy, an equal part results from the ambiguity over what species are expected to represent by the many who use them. Here, mtDNA data, microsatellite data, and sequence data from 17 nuclear loci are used in an integrated and quantitative manner to resolve the presence of evolutionary lineages, their contemporary and historical structure, and their correspondence to species, in a species complex of Amazonian peacock "bass" cichlids (Cichla pinima sensu lato). Results suggest that the historical narrative for these populations is more complex than can be portrayed by recognizing them as one, two, or four species: their history and contemporary dynamics cannot be unambiguously rendered as discrete units (taxa) at any level without both choosing the supremacy of one delimitation criterion and obscuring the very information that provides insight into the diversification process. This calls into question the utility of species as a rank, term, or concept, and suggests that while biologists may have a reasonable grasp of the structure of evolution, our methods of conveying these insights need updating. The lack of correspondence between evolutionary phenomena and discrete species should serve as a null hypothesis, and researchers should focus on quantifying the diversity in nature at whatever hierarchical level it occurs.

Population genomic analyses reveal a history of range expansion and trait evolution across the native and invaded range of yellow starthistle (Centaurea solstitialis)

Identifying sources of genetic variation and reconstructing invasion routes for non-native introduced species is central to understanding the circumstances under which they may evolve increased invasiveness. In this study, we used genome-wide single nucleotide polymorphisms to study the colonization history of Centaurea solstitialis in its native range in Eurasia and invasions into the Americas. We leveraged this information to pinpoint key evolutionary shifts in plant size, a focal trait associated with invasiveness in this species. Our analyses revealed clear population genomic structure of potential source populations in Eurasia, including deep differentiation of a lineage found in the southern Apennine and Balkan Peninsulas and divergence among populations in Asia, eastern Europe and western Europe. We found strongest support for an evolutionary scenario in which western European populations were derived from an ancient admixture event between populations from eastern Europe and Asia, and subsequently served as the main genetic 'bridgehead' for introductions to the Americas. Introductions to California appear to be from a single source region, and multiple, independent introductions of divergent genotypes likely occurred into the Pacific Northwest. Plant size has evolved significantly at three points during range expansion, including a large size increase in the lineage responsible for the aggressive invasion of the California interior. These results reveal a long history of colonization, admixture and trait evolution in C. solstitialis, and suggest routes for improving evidence-based management decisions for one of the most ecologically and economically damaging invasive species in the western United States.

Phylogeography and postglacial colonization of Central Europe by Anguis fragilis and Anguis colchica

Five slow-worm species are distributed in the Palearctic region. Two species, Anguis fragilis and A. colchica, are native to Central Europe, where only limited information on the phylogeography of the genus exists. Here, we examined the situation in Poland and surrounding regions, where a mitochondrial contact zone between the species is expected. We used new mitochondrial DNA sequences and available published data from Central Europe and the northern Balkans. Haplotypes of both species were recorded in the study area. Anguis fragilis is represented by a single haplogroup, while A. colchica by three haplogroups. This suggests four independent sources/refugia for postglacial colonization of northern Central Europe. The mitochondrial contact zone seems to mirror the borders between lowlands of the North European Plain and East European Plain, and the south-eastern Poland uplands, while the Vistula River does not constitute a barrier. The presence of both species, A. fragilis and A. colchica, in Poland should also be considered for protection by the Polish conservation legislation.

Phylogeography reveals an ancient cryptic radiation in East-Asian tree frogs (Hyla japonica group) and complex relationships between continental and island lineages

Background

In contrast to the Western Palearctic and Nearctic biogeographic regions, the phylogeography of Eastern-Palearctic terrestrial vertebrates has received relatively little attention. In East Asia, tectonic events, along with Pleistocene climatic conditions, likely affected species distribution and diversity, especially through their impact on sea levels and the consequent opening and closing of land-bridges between Eurasia and the Japanese Archipelago. To better understand these effects, we sequenced mitochondrial and nuclear markers to determine phylogeographic patterns in East-Asian tree frogs, with a particular focus on the widespread H. japonica.

Results

We document several cryptic lineages within the currently recognized H. japonica populations, including two main clades of Late Miocene divergence (~5 Mya). One occurs on the northeastern Japanese Archipelago (Honshu and Hokkaido) and the Russian Far-East islands (Kunashir and Sakhalin), and the second one inhabits the remaining range, comprising southwestern Japan, the Korean Peninsula, Transiberian China, Russia and Mongolia. Each clade further features strong allopatric Plio-Pleistocene subdivisions (~2-3 Mya), especially among continental and southwestern Japanese tree frog populations. Combined with paleo-climate-based distribution models, the molecular data allowed the identification of Pleistocene glacial refugia and continental routes of postglacial recolonization. Phylogenetic reconstructions further supported genetic homogeneity between the Korean H. suweonensis and Chinese H. immaculata, suggesting the former to be a relic population of the latter that arose when the Yellow Sea formed, at the end of the last glaciation.

Conclusions

Patterns of divergence and diversity were likely triggered by Miocene tectonic activities and Quaternary climatic fluctuations (including glaciations), causing the formation and disappearance of land-bridges between the Japanese islands and the continent. Overall, this resulted in a ring-like diversification of H. japonica around the Sea of Japan. Our findings urge for important taxonomic revisions in East-Asian tree frogs. First, they support the synonymy of H. suweonensis (Kuramoto, 1980) and H. immaculata (Boettger, 1888). Second, the nominal H. japonica (Günther, 1859) represents at least two species: an eastern (new taxon A) on the northern Japanese and Russian Far East islands, and a southwestern species (n. t. B) on southern Japanese islands and possibly also forming continental populations. Third, these continental tree frogs may also represent an additional entity, previously described as H. stepheni Boulenger, 1888 (senior synonym of H. ussuriensis Nikolskii, 1918). A complete revision of this group requires further taxonomic and nomenclatural analyses, especially since it remains unclear to which taxon the species-epitheton japonica corresponds to.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0814-x) contains supplementary material, which is available to authorized users.