Tracing glacial refugia of Triturus newts based on mitochondrial DNA phylogeography and species distribution modeling

Spatial genetic structure in European marbled newts revealed with target enrichment by sequence capture

European marbled newts come in two species that have abutting ranges. The northern species, Triturus marmoratus, is found in France and the northern part of the Iberian Peninsula, whereas the southern species, T. pygmaeus, is found in the southwestern corner of the Iberian Peninsula. We study the intraspecific genetic differentiation of the group because morphological data show geographical variation and because the Iberian Peninsula is a recognized center of speciation and intraspecific genetic diversity for all kinds of organisms, amphibians included. We use target enrichment by sequence capture to generate c. 7 k nuclear DNA markers. We observe limited genetic exchange between the species, which confirms their distinctiveness. Both species show substantial genetic structuring that is only in part mirrored by morphological variation. Genetically differentiated groups are found in the south (T. marmoratus) and west (T. pygmaeus) of the species ranges. Our observations highlight the position of the Iberian Peninsula as a hotspot for genetic differentiation.

A synoptic review of the Amphibians of Iran: bibliography, taxonomy, synonymy, distribution, conservation status, and identification key to the eggs, larvae, and adults

This study provides an illustrated account, a comprehensive update of the systematics, and a bibliography of the 15 species of anurans in five families, eight genera; and of the six species of urodeles in two families, four genera in Iran. Bufonidae, with eight species, is the most diverse family; Salamandridae has five species and Ranidae has four species. This study also presents updated identification keys for the eggs, larvae, and metamorphosed amphibians of Iran. We designated specimen NMW 19855.1 as neotype of Pelophylax persicus (Schneider, 1799) comb. nov.. Along with distribution maps obtained from all the reliable localities and museum specimens known at this time, the modelled habitat of species, and for the first time, the National Red List of amphibians based on the IUCN red list categories and criteria. Based on our evaluation we propose to categorize Bufo eichwaldi, Paradactylodon persicus, Neurergus derjugini, and N. kaiseri as Vulnerable at National Red List, and to move Bufotes (Calliopersa) luristanicus, B. (C.) surdus, Firouzophrynus olivaceus, and Rana pseudodalmatina from the category of Least Concern (LC) to Near Threatened (NT). The National Red List of amphibians that we propose has significant implications for endangered species management and conservation. Forty-one percent of amphibian species in Iran are endemic to the country, and more than forty percent of the Iranian amphibians are at risk of extinction. Zagros Mountain forest and Hyrcaniain forests have more than 80% (i.e. 18 species) of the diversity of Iranian amphibians. A considerable amount of scientific literature published on Iranian amphibians in Persian language is not easily accessible to researchers outside Iran. This monograph attempts to remedy the situation and provides broader access to international herpetology. We recognize that taxonomy is always in a state of flux, and the names and synonymies used here reflect our current view.

A two-species distribution model for parapatric newts, with inferences on their history of spatial replacement

Related species often engage in abutting or overlapping contact zones with various strengths of interspecific competition. Biotic interactions such as these preclude the registration of the full profile of environmental variables that would describe the otherwise larger species ranges. Here, I advocate to forego full range species distribution modelling and instead focus on the ecography of the contact zone, for example with ‘two-species distribution models’ (TSDMs), in which presence data are contrasted against the background of environmental data. The newts Triturus cristatus and Triturus marmoratus meet in the west of France. A countrywide TSDM suggests that the contact zone of the species is located at a climatic gradient, in line with their north-eastern vs. south-western ranges. The species are also ecologically segregated by elevation and forestation, which is in line with a documented movement of the contact zone caused by hedgerow removal in lowland areas. Hindcasts for the Holocene suggest that the species contact zone was positioned at either the same place as at present or more to the south, depending on the amount of forestation. A forecast under climate warming predicts a fast movement to the north, but this scenario is deemed unrealistic. One reason is that recent habitat loss compromises dispersal and range expansion. Other species pairs to which TSDMs have been applied are listed for comparison.

Population structure and adaptive variation of Helichrysum italicum (Roth) G. Don along eastern Adriatic temperature and precipitation gradient

Immortelle (Helichrysum italicum (Roth) G. Don; Asteraceae) is a perennial plant species native to the Mediterranean region, known for many properties with wide application mainly in perfume and cosmetic industry. A total of 18 wild H. italicum populations systematically sampled along the eastern Adriatic environmental gradient were studied using AFLP markers to determine genetic diversity and structure and to identify loci potentially responsible for adaptive divergence. Results showed higher levels of intrapopulation diversity than interpopulation diversity. Genetic differentiation among populations was significant but low, indicating extensive gene flow between populations. Bayesian analysis of population structure revealed the existence of two genetic clusters. Combining the results of F_ST - outlier analysis (Mcheza and BayeScan) and genome-environment association analysis (Samβada, LFMM) four AFLP loci strongly associated with the bioclimatic variables Bio03 Isothermality, Bio08 Mean temperature of the wettest quarter, Bio15 Precipitation seasonality, and Bio17 Precipitation of driest quarter were found to be the main variables driving potential adaptive genetic variation in H. italicum along the eastern Adriatic environmental gradient. Redundancy analysis revealed that the partitioning of genetic variation was mainly associated with the adaptation to temperature oscillations. The results of the research may contribute to a clearer understanding of the importance of local adaptations for the genetic differentiation of Mediterranean plants and allow the planning of appropriate conservation strategies. However, considering that the identified outlier loci may be linked to genes under selection rather than being the target of natural selection, future studies must aim at their additional analysis.

Genetic and morphological data demonstrate hybridization and backcrossing in a pair of salamanders at the far end of the speciation continuum

Deeply diverged yet hybridizing species provide a system to investigate the final stages of the speciation process. We study a hybridizing pair of salamander species-the morphologically and genetically drastically different newts Triturus cristatus and T. marmoratus-with a panel of 32 nuclear and mitochondrial genetic markers. Morphologically identified hybrids are mostly of the F1 generation and mothered by T. cristatus. The sex ratio of the F1 hybrid class is reciprocally skewed, with a preponderance of females in T. cristatus-mothered hybrids and males in T. marmoratus-mothered hybrids. This amounts to the Haldane effect operating in one direction of the cross. Deeper generation hybrids are occasionally produced, possibly including F1 hybrid × backcross hybrid offspring. Interspecific gene flow is low, yet skewed toward T. cristatus. This asymmetry may be caused by hybrid zone movement, with the superseding species being predisposed to introgression. The persisting gene flow between deeply differentiated species supports the notion that full genetic isolation may be selected against. Conversely, published morphological data suggest that introgressive hybridization is detrimental, with digital malformations occurring more frequently in the area of sympatry. Finally, to assist field identification, both within the area of natural range overlap and concerning anthropogenic introductions elsewhere, we document the phenotypical variation of two generations of hybrids compared with both parental species. We suggest that fluctuating range boundaries, ecological segregation, cytonuclear incompatibilities and hybrid breakdown through Bateson-Dobzhansky-Muller incompatibilities all contribute to species integrity, despite incomplete isolation during secondary contact.

Pleistocene allopatric differentiation followed by recent range expansion explains the distribution and molecular diversity of two congeneric crustacean species in the Palaearctic

Pleistocene glaciations had a tremendous impact on the biota across the Palaearctic, resulting in strong phylogeographic signals of range contraction and rapid postglacial recolonization of the deglaciated areas. Here, we explore the diversity patterns and history of two sibling species of passively dispersing taxa typical of temporary ponds, fairy shrimps (Anostraca). We combine mitochondrial (COI) and nuclear (ITS2 and 18S) markers to conduct a range-wide phylogeographic study including 56 populations of Branchinecta ferox and Branchinecta orientalis in the Palaearctic. Specifically, we investigate whether their largely overlapping ranges in Europe resulted from allopatric differentiation in separate glacial refugia followed by a secondary contact and reconstruct their postglacial recolonization from the inhabited refugia. Our results suggest the existence of distinct refugia for the two species, with genetic divergence among intraspecific lineages consistent with late Pleistocene glacial cycles. While B. ferox lineages originated from Mediterranean refugia, the origin of B. orientalis lineages was possibly located on the Pannonian Plain. We showed that most dispersal events predominantly happened within 100 km, coupled with several recent long-distance events (> 1000 km). Hence the regional habitat density of suitable habitats in Central Europe is possibly a key to the co-existence of the two species. Overall, our study illustrates how isolation in combination with stochastic effects linked to glacial periods are important drivers of the allopatric differentiation of Palaearctic taxa.

Batrachochytrium salamandrivorans Threat to the Iberian Urodele Hotspot

The recent introduction of the chytrid fungus Batrachochytrium salamandrivorans into northeastern Spain threatens salamander diversity on the Iberian Peninsula. We assessed the current epidemiological situation with extensive field sampling of urodele populations. We then sought to delineate priority regions and identify conservation units for the Iberian Peninsula by estimating the susceptibility of Iberian urodeles using laboratory experiments, evidence from mortality events in nature and captivity and inference from phylogeny. None of the 1395 field samples, collected between 2015 and 2021 were positive for Bsal and no Bsal-associated mortality events were recorded, in contrast to the confirmed occurrence of Bsal outbreak previously described in 2018. We classified five of eleven Iberian urodele species as highly susceptible, predicting elevated mortality and population declines following potential Bsal emergence in the wild, five species as intermediately susceptible with variable disease outcomes and one species as resistant to disease and mortality. We identified the six conservation units (i.e., species or lineages within species) at highest risk and propose priority areas for active disease surveillance and field biosecurity measures. The magnitude of the disease threat identified here emphasizes the need for region-tailored disease abatement plans that couple active disease surveillance to rapid and drastic actions.

Stability in the Balkans: phylogeography of the endemic Greek stream frog, Rana graeca

We still have little knowledge concerning the phylogeography of amphibians and reptiles from the Balkan Peninsula compared with the other two Mediterranean peninsulas. This raises concerns for endemic taxa from these peninsulas, because it might interfere with further conservation efforts. Here we focus on the endemic Greek stream frog (Rana graeca) and reconstruct its biogeography and evolutionary history. Using four genetic markers (Cytb, 16S, COI and BDNF) in > 350 sequences covering the whole distribution range, we conducted phylogenetic, demographic and ecological niche analyses, which revealed the phylogeography of this species. Surprisingly, this examination of R. graeca reveals a very shallow level of intraspecific genetic variability through the Balkans, with two main, statistically supported lineages having a partly sympatric distribution. The most variable marker was Cytb, which showed 19 haplotypes in 123 analysed sequences in the whole species distribution area. Here presented genetic data, together with the environmental niche projection and demographic analyses suggest that R. graeca was probably affected only marginally by climatic oscillations, with the Hellenides as the most suitable area for the occurrence of the species in different geological periods. This is consistent with the observed genetic diversity, which is mostly related to these mountains. Although the species shows a certain level of phenotypic variability and ecological preferences, this might be related to species plasticity affected by the micro-climatic conditions in small areas, which merits further research. Comparing phylogeography of other amphibian and reptile species in the Balkans, we showed that the observed pattern represents a new view on the phylogeography of the Balkan herpetofauna.

Niche differentiation between deeply divergent phylogenetic lineages of an endemic newt: implications for Species Distribution Models

Species distribution models (SDMs) treat species as a single unit, neglecting intraspecific variation. Few studies address the impact of intraspecific variation on SDM performance, and none of them account for the well-known inter-algorithm variability in prediction performance. The endemic Iberian amphibian Lissotriton boscai comprises two geographically highly structured phylogenetic lineages, which allowed us to explore how intraspecific variation affects the overall performance of SDMs and the predicted ecological niche. We built species and lineage distribution models using three different presence-only algorithms. We also tested for niche overlap, niche equivalency and niche similarity, using an ordination technique. We found differences in the predicted potential distribution of the two lineages and the underlying environmental factors. Moreover, intraspecific differences in model predictive capacity existed irrespective of which algorithm was used to build the distribution models. This was coupled with lineages showing a low degree of niche overlap and occurring in relatively different environmental niches spaces. The intraspecific variation observed in L. boscai led to an improved intraspecific predictivity of the lineage level based-distribution models. There was partial spatial agreement between the niche overlap and independently reported secondary contact zones. Thus SDMs built only at the species level may be too naive to predict impacts of global change on species distributions.

Genetic Divergence Across Glacial Refugia Despite Interglacial Gene Flow in a Crested Newt

MtDNA-based phylogeography has illuminated the impact of the Pleistocene Ice Age on species distribution dynamics and the build-up of genetic divergence. The well-known shortcomings of mtDNA in biogeographical inference can be compensated by integrating multilocus data and species distribution modelling into phylogeography. We re-visit the phylogeography of the Italian crested newt (Triturus carnifex), a species distributed in two of Europe’s main glacial refugia, the Balkan and Italian Peninsulas. While a new 51 nuclear DNA marker dataset supports the existence of three lineages previously suggested by mtDNA (Balkan, northern Italy and southern Italy), the nuclear DNA dataset also provides improved resolution where these lineages have obtained secondary contact. We observe geographically restricted admixture at the contact between the Balkan and northern Italy gene pools and identify a potential mtDNA ghost lineage here. At the contact between the northern and southern Italy gene pools we find admixture over a broader area, as well as asymmetric mtDNA introgression. Our species distribution model is in agreement with a distribution restricted to distinct refugia during Pleistocene glacial cycles and postglacial expansion with secondary contact. Our study supports: (1) the relevance of the north-western Balkan Peninsula as a discrete glacial refugium; (2) the importance of north-eastern Italy and the northern Apennine as suture zones; and (3) the applicability of a refugia-within-refugia scenario within the Italian Peninsula.

Reconstruction of past distribution for the Mongolian toad, Strauchbufo raddei (Anura: Bufonidae) using environmental modeling

The use of ecological models enables determining the current distribution of species, but also their past distribution when matching climatic conditions are available. In specific cases, they can also be used to determine the likelihood of fossils to belong to the same species-under the hypothesis that all individuals of a species have the same ecological requirements. Here, using environmental modeling, we reconstructed the distribution of the Mongolian toad, Strauchbufo raddei, since the Last Glacial Maximum and thus covering the time period between the Late Pleistocene and the Holocene. We found the range of the species to have shifted over time, with the LGM population clustered around the current southern range of the species, before expanding east and north during the Pleistocene, and reaching the current range since the mid-Holocene. Finally, we determined that the ecological conditions during the life-time of the mid-Pleistocene fossils attributed to the species in Europe were too different from the one of the extant species or fossils occurring at the same period in Asia to belong to the same species.

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.

Post-glacial phylogeography and variation in innate immunity loci in a sylvatic rodent, bank vole Myodes glareolus

In the northern hemisphere, the spatial structure of many taxa has been shaped by migration patterns after the last glaciation, and phylogeography based on mtDNA variation may reflect the post-glacial demography. The mtDNA lineages are expected to differ in their adaptations to local conditions but little is known about the impact of these conditions on functional genetic variation. Here, we answer this question through an analysis of geographic variation and selection patterns in seven innate immunity genes in free-living bank voles Myodes glareolus from 10 localities across species range assigned to different lineages based on mtDNA. We found clear discrepancies between population structure in mtDNA and each of the studied innate immunity genes. There was no uniform pattern of spatial variation at immunity loci, they differed in the levels of polymorphism, and the results of neutrality tests were not consistent over loci. Each locus comprised a few common haplotypes shared between mitochondrial lineages and studied locations, plus numerous haplotypes unique for each studied site. Our results suggest that the diversity of innate immunity genes cannot be explained solely in terms of demographic processes, and that the observed polymorphism may be attributed to local selection. The strength and direction of selection differed between loci, even within the same gene family, which underlines how crucial it is to take a complex approach while studying the selection patterns acting on immune-related genes.

Phylogeography of Potamon ibericum (Brachyura: Potamidae) identifies Quaternary glacial refugia within the Caucasus biodiversity hot spot

Refugia are critical for the maintenance of biodiversity during the periods of Quaternary climatic oscillations. The long-term persistence of refugial populations in a large continuous refugium has resulted in a homogenous pattern of genetic structure among populations, while highly structured evolutionary lineages characterize the restriction of refugial populations to smaller subrefugia. These mechanisms have resulted in the identification of hot spots of biodiversity within putative glacial refugia. We studied phylogeography of Potamon ibericum (Brachyura: Potamidae) in the drainages of the western Caucasus biodiversity hot spot (i.e., Colchis and the Caucasus) to infer spatial genetic structure and potential refugia for a freshwater crab in this region. These areas have traditionally considered as a refugium due to the presence of Tertiary relict species. We integrated population genetic data and historical demographic analysis from cytochrome oxidase subunit I sequences and paleoclimatic data from species distribution modeling (SDM). The results revealed the lack of phylogeographic structure and provided evidence for demographic expansion. The SDM presented a rather homogenous and large refugium that extended from northeast Turkey to Colchis during the last glacial period. In contrast to these findings, previous phylogeographic study on P. ibericum of the eastern Caucasus biodiversity hot spot (i.e., Hyrcania) identified multiple independent refugia. By combining these results, we explain the significance of this important western Palearctic hot spot of biological diversity in shaping the geographic distribution of intraspecific genetic diversity in a freshwater taxon.

Non-native freshwater fauna in Portugal: A review

We present the most updated list of non-native freshwater fauna established in Portugal, including the Azores and Madeira archipelagos. This list includes 67 species at national level but corresponds to 84 species records, of which 53 are in the mainland, 23 in the Azores and 8 in Madeira archipelagos. We also discuss the progression of the cumulative number of introductions since 1800 and identify the most probable vectors of introduction, main taxonomic groups and their regions of origin. Furthermore, we review the existing knowledge about ecological and economic impacts, invasion risk and potential distribution of invaders, under present and future climatic conditions, and the applied management actions, including the production of legislation. Along the 20th century the number of successful introductions increased at an approximate rate of two new species per decade until the beginning of 1970s. Since then, this rate increased to about 14 new species per decade. These introductions were mainly a result of fisheries, as contaminants or for ornamental purposes. Fish and mollusks are the taxonomic groups with more established species, representing more than half of the total. Most species (>70%) are native from other regions of Europe and North America. Studies about ecological or socioeconomic impacts are more common for fish, crustaceans and mollusks. Impacts for most amphibians, reptiles and mammals are not thoroughly studied. A few studies on the impacts and management actions of health-threatening mosquitoes are also available. The potential distribution in the Portuguese territory was modelled for 26 species. Only a minority of these models provides projections of distributions under scenarios of future climate change. A comparison of the Portuguese and EU legislation shows large discrepancies in the invasive species lists. Using the EU list and a ranking procedure for the national context, we identify freshwater species of high national concern for which actions are urgently needed.

Landscape genetics reveal broad and fine-scale population structure due to landscape features and climate history in the northern leopard frog (Rana pipiens) in North Dakota

Prehistoric climate and landscape features play large roles structuring wildlife populations. The amphibians of the northern Great Plains of North America present an opportunity to investigate how these factors affect colonization, migration, and current population genetic structure. This study used 11 microsatellite loci to genotype 1,230 northern leopard frogs (Rana pipiens) from 41 wetlands (30 samples/wetland) across North Dakota. Genetic structure of the sampled frogs was evaluated using Bayesian and multivariate clustering methods. All analyses produced concordant results, identifying a major east-west split between two R. pipiens population clusters separated by the Missouri River. Substructuring within the two major identified population clusters was also found. Spatial principal component analysis (sPCA) and variance partitioning analysis identified distance, river basins, and the Missouri River as the most important landscape factors differentiating R. pipiens populations across the state. Bayesian reconstruction of coalescence times suggested the major east-west split occurred ~13-18 kya during a period of glacial retreat in the northern Great Plains and substructuring largely occurred ~5-11 kya during a period of extreme drought cycles. A range-wide species distribution model (SDM) for R. pipiens was developed and applied to prehistoric climate conditions during the Last Glacial Maximum (21 kya) and the mid-Holocene (6 kya) from the CCSM4 climate model to identify potential refugia. The SDM indicated potential refugia existed in South Dakota or further south in Nebraska. The ancestral populations of R. pipiens in North Dakota may have inhabited these refugia, but more sampling outside the state is needed to reconstruct the route of colonization. Using microsatellite genotype data, this study determined that colonization from glacial refugia, drought dynamics in the northern Great Plains, and major rivers acting as barriers to gene flow were the defining forces shaping the regional population structure of R. pipiens in North Dakota.

Evidences for a shared history for spectacled salamanders, haplotypes and climate

The so-called glacial refugia, formed during the Pleistocene climatic oscillations, played a major role in shaping the distribution of European species, triggering migrations or isolating populations. Many of these events were recently investigated by genetic data, mainly for the European Last Glacial stage, in the Iberic, Italian and Greek-Balkan peninsulas. The amphibian genus Salamandrina, the most ancient living salamandrid lineage, was widespread in Europe until the climatic oscillations of Miocene probably forced it to shelter in the only suitable territory at that time, the Apennines. Nowadays this genus is endemic of peninsular Italy with two parapatric species, S. perspicillata and S. terdigitata, sharing an area of secondary contact formed after the Last Glacial Maximum. Climate is generally identified as the key factor for the interpretation of genetic data. In this research, we directly measure climate influences on the two Salamandrina known species through Ensemble Modelling techniques and post-modelling GIS analyses, integrating updated genetic data in this process. Our results confirm the hypotheses of southwards (and subsequent northwards) shifts, identify glacial refugia and corridors used for the post-glacial re-colonization. Finally, we map a contact zone deserving more sampling effort to disentangle the introgression and hybridization observed.

Oxidative cost of interspecific hybridization: a case study of two Triturus species and their hybrids

Oxidative stress has most recently been suggested as one of the possible mechanisms responsible for reduced fitness of hybrids. To explore possible oxidative cost of hybridization, we examined anti-oxidant defence system parameters (superoxide dismutase, catalase, glutathione peroxidase, glutathione S-transferase, glutathione reductase, glutathione, SH groups), their interconnectedness (index of integration) and levels of oxidative damage [concentrations of lipid peroxides, TBARS (thiobarbituric acid reactive substances)] in laboratory-reared newt species, Triturus macedonicus and T. ivanbureschi, and their hybrid. Our results showed that parental species differed in anti-oxidant defence system parameters, but not in the levels of integration of the whole system and oxidative damage. Individuals of T. ivanbureschi had higher activities of superoxide dismutase, glutathione S-transferase and concentrations of glutathione. Hybrid individuals of crested newts displayed higher levels of the anti-oxidant defence system (higher superoxide dismutase, catalase, glutathione peroxidase activities and concentrations of SH groups), and a lower overall correlation of anti-oxidant system (lower index of integration) in comparison with both parental species, suggesting that they may possess a less efficient anti-oxidant defence system and a higher investment in maintaining oxidative balance. The higher investment in the anti-oxidant system could divert limited resources away from other functions and affect further hybrid fitness. The presented findings contribute to a better understanding of the anti-oxidant defence system of crested newts and their interspecies differences, and support the hypothesis that oxidative stress is one of the costs of interspecific hybridization.

Examining the interglacial high-elevation refugia scenario in East Asian subtropical mountain systems with the frog species Leptobrachium liui

The effects of Quaternary climatic oscillations on the distributions of organisms in different parts of the world are not equally well understood, limiting the ability to understand the determinants of biodiversity. Compared with the mountain regions in southern Europe and southwestern North America, such effects on high-elevation species in the East Asian subtropical mountain systems located in southern and southeastern China have seldom been addressed. In this study, using Leptobrachium liui (Megophryidae), we made one of the earliest attempts to examine the interglacial high-elevation refugia scenario in these Asian mountains. Based on our current understanding of the study system, we formulated a hypothesis that these frogs of western origin were distributed more widely and continuously during glacial phases, allowing eastward dispersal, and that they are currently isolated in interglacial refugia at higher elevations. Microsatellite data and mitochondrial and nuclear sequence data were obtained with extensive sampling followed by the synthesis of phylogeographic and population genetic analyses and modeling of the species distribution. The analyses revealed a sequential eastward divergence of microsatellite clusters and gene lineages accompanied by a decline in genetic diversity. Molecular dating estimates revealed divergence events in the Pleistocene, and a reduction in local populations was inferred to have occurred at a time comparable to the end of the last glacial. Strong genetic isolation by distance reflecting a more continuous historical distribution was detected. Furthermore, environmental niche models inferred a wide planar distribution during the last glacial maximum, providing further support for the hypothesis.

Absence of heterosis in hybrid crested newts

Relationships between phylogenetic relatedness, hybrid zone spatial structure, the amount of interspecific gene flow and population demography were investigated, with the newt genus Triturus as a model system. In earlier work, a bimodal hybrid zone of two distantly related species combined low interspecific gene flow with hybrid sterility and heterosis was documented. Apart from that, a suite of unimodal hybrid zones in closely related Triturus showed more or less extensive introgressive hybridization with no evidence for heterosis. We here report on population demography and interspecific gene flow in two Triturus species (T. macedonicus and T. ivanbureschi in Serbia). These are two that are moderately related, engage in a heterogeneous uni-/bimodal hybrid zone and hence represent an intermediate situation. This study used 13 diagnostic nuclear genetic markers in a population at the species contact zone. This showed that all individuals were hybrids, with no parentals detected. Age, size and longevity and the estimated growth curves are not exceeding that of the parental species, so that we conclude the absence of heterosis in T. macedonicus-T. ivanbureschi. Observations across the genus support the hypothesis that fertile hybrids allocate resources to reproduction and infertile hybrids allocate resources to growth. Several Triturus species hybrid zones not yet studied allow the testing of this hypothesis.

Coupling GIS spatial analysis and Ensemble Niche Modelling to investigate climate change-related threats to the Sicilian pond turtle Emys trinacris, an endangered species from the Mediterranean

The pond turtle Emys trinacris is an endangered endemic species of Sicily showing a fragmented distribution throughout the main island. In this study, we applied "Ensemble Niche Modelling", combining more classical statistical techniques as Generalized Linear Models and Multivariate Adaptive Regression Splines with machine-learning approaches as Boosted Regression Trees and Maxent, to model the potential distribution of the species under current and future climatic conditions. Moreover, a "gap analysis" performed on both the species' presence sites and the predictions from the Ensemble Models is proposed to integrate outputs from these models, in order to assess the conservation status of this threatened species in the context of biodiversity management. For this aim, four "Representative Concentration Pathways", corresponding to different greenhouse gases emissions trajectories were considered to project the obtained models to both 2050 and 2070. Areas lost, gained or remaining stable for the target species in the projected models were calculated. E. trinacris' potential distribution resulted to be significantly dependent upon precipitation-linked variables, mainly precipitation of wettest and coldest quarter. Future negative effects for the conservation of this species, because of more unstable precipitation patterns and extreme meteorological events, emerged from our analyses. Further, the sites currently inhabited by E. trinacris are, for more than a half, out of the Protected Areas network, highlighting an inadequate management of the species by the authorities responsible for its protection. Our results, therefore, suggest that in the next future the Sicilian pond turtle will need the utmost attention by the scientific community to avoid the imminent risk of extinction. Finally, the gap analysis performed in GIS environment resulted to be a very informative post-modeling technique, potentially applicable to the management of species at risk and to Protected Areas' planning in many contexts.

The study of larval tail morphology reveals differentiation between two Triturus species and their hybrids

In amphibians, morphological differentiation and disparity at the larval and post-metamorphic ontogenetic stages can diverge, owing to various contrasting environments and different selective pressures. In the monophyletic clade of nineTriturusnewt species, five different morphotypes can be recognized, but information on larval morphology is limited. Here we explore divergence of larval morphology inTriturus ivanbureschi,T. macedonicus, and their F1 hybrids. These two genetically and morphologically distinct crested newt species hybridize in nature and form a relatively wide hybrid zone in the central part of the Balkan Peninsula. Using a geometric morphometric approach and multivariate statistics, we evaluated differences of tail size and shape, colouration pattern, and the presence of a tail filament at the mid-larval stage in larvae reared under controlled laboratory conditions. We chose the tail as the main propulsive organ crucial for locomotion, feeding, and escaping predators. We found thatTriturus ivanbureschiandT. macedonicuslarvae differ in tail shape, but not in tail size. Two groups of F1 hybrid larvae (obtained from reciprocal crossing) were similar to each other, but differed from the parental species in size and shape of the tail, colouration pattern, and the presence of a tail filament. Our results indicate that, like adults, larvae diverge morphologically and hybrid larvae do not exhibit intermediate morphology of the parental species.

A signature of dynamic biogeography: enclaves indicate past species replacement

Understanding how species have replaced each other in the past is important to predicting future species turnover. While past species replacement is difficult to detect after the fact, the process may be inferred from present-day distribution patterns. Species with abutting ranges sometimes show a characteristic distribution pattern, where a section of one species range is enveloped by that of the other. Such an enclave could indicate past species replacement: when a species is partly supplanted by a competitor, but a population endures locally while the invading species moves around and past it, an enclave forms. If the two species hybridize and backcross, the receding species is predicted to leave genetic traces within the expanding one under a scenario of species replacement. By screening dozens of genes in hybridizing crested newts, we uncover genetic remnants of the ancestral species, now inhabiting an enclave, in the range of the surrounding invading species. This independent genetic evidence supports the past distribution dynamics we predicted from the enclave. We suggest that enclaves provide a valuable tool in understanding historical species replacement, which is important because a major conservation concern arising from anthropogenic climate change is increased species replacement in the future.

Unraveling climate influences on the distribution of the parapatric newts Lissotriton vulgaris meridionalis and L. italicus

Background

Climate is often considered as a key ecological factor limiting the capability of expansion of most species and the extent of suitable habitats. In this contribution, we implement Species Distribution Models (SDMs) to study two parapatric amphibians, Lissotriton vulgaris meridionalis and L. italicus, investigating if and how climate has influenced their present and past (Last Glacial Maximum and Holocene) distributions. A database of 901 GPS presence records was generated for the two newts. SDMs were built through Boosted Regression Trees and Maxent, using the Worldclim bioclimatic variables as predictors.

Results

Precipitation-linked variables and the temperature annual range strongly influence the current occurrence patterns of the two Lissotriton species analyzed. The two newts show opposite responses to the most contributing variables, such as BIO7 (temperature annual range), BIO12 (annual precipitation), BIO17 (precipitation of the driest quarter) and BIO19 (precipitation of the coldest quarter). The hypothesis of climate influencing the distributions of these species is also supported by the fact that the co-occurrences within the sympatric area fall in localities characterized by intermediate values of these predictors. Projections to the Last Glacial Maximum and Holocene scenarios provided a coherent representation of climate influences on the past distributions of the target species. Computation of pairwise variables interactions and the discriminant analysis allowed a deeper interpretation of SDMs’ outputs. Further, we propose a multivariate environmental dissimilarity index (MEDI), derived through a transformation of the multivariate environmental similarity surface (MESS), to deal with extrapolation-linked uncertainties in model projections to past climate. Finally, the niche equivalency and niche similarity tests confirmed the link between SDMs outputs and actual differences in the ecological niches of the two species.

Conclusions

The different responses of the two species to climatic factors have significantly contributed to shape their current distribution, through contractions, expansions and shifts over time, allowing to maintain two wide allopatric areas with an area of sympatry in Central Italy. Moreover, our SDMs hindcasting shows many concordances with previous phylogeographic studies carried out on the same species, thus corroborating the scenarios of potential distribution during the Last Glacial Maximum and the Holocene emerging from the models obtained.

Electronic supplementary material

The online version of this article (10.1186/s12983-017-0239-4) contains supplementary material, which is available to authorized users.

Genetic diversity of Diaphorina citri and its endosymbionts across east and south-east Asia

BACKGROUND

Diaphorina citri is the vector of 'Candidatus Liberibacter asiaticus', the most widespread pathogen associated huanglongbing, the most serious disease of citrus. To enhance our understanding of the distribution and origin of the psyllid, we investigated the genetic diversity and population structures of 24 populations in Asia and one from Florida based on the mtCOI gene. Simultaneously, genetic diversity and population structures of the primary endosymbiont (P-endosymbiont) 'Candidatus Carsonella ruddii' and secondary endosymbiont (S-endosymbiont) 'Candidatus Profftella armatura' of D. citri were determined with the housekeeping genes.

RESULT

AMOVA analysis indicated that populations of D. citri and its endosymbionts in east and south-east Asia were genetically distinct from populations in Pakistan and Florida. Furthermore, P-endosymbiont populations displayed a strong geographical structure across east and south-east Asia, while low genetic diversity indicated the absence of genetic structure among the populations of D. citri and its S-endosymbiont across these regions.

CONCLUSION

The 'Ca. C. ruddii' is more diverse and structured than the D. citri and the 'Ca. P. armatura' across east and south-east Asia. Multiple introductions of the psyllid have occurred in China. Management application for controlling the pest is proposed based on the genetic information of D. citri and its endosymbionts. © 2017 Society of Chemical Industry.

A genomic footprint of hybrid zone movement in crested newts

Speciation typically involves a stage in which species can still exchange genetic material. Interspecific gene flow is facilitated by the hybrid zones that such species establish upon secondary contact. If one member of a hybridizing species pair displaces the other, their hybrid zone would move across the landscape. Although theory predicts that moving hybrid zones quickly stagnate, hybrid zones tracked over one or a few decades do not always follow such a limitation. This suggests that hybrid zones have the potential to traverse considerable distances over extended periods of time. When hybrid zones move, introgression is predicted to result in biased gene flow of selectively neutral alleles, from the receding species into the advancing species. We test for such a genomic footprint of hybrid zone movement in a pair of crested newt species (genus Triturus) for which we have a priori support for westward hybrid zone movement. We perform a multilocus phylogeographical survey and conduct Bayesian clustering analysis, estimation of ancestry and heterozygosity, and geographical cline analysis. In a 600 km wide area east of the present day hybrid zone a genomic footprint constitutes empirical evidence consistent with westward hybrid zone movement. The crested newt case suggests that hybrid zone movement can occur over an extensive span of time and space. Inferring hybrid zone movement provides fundamental insight into historical biogeography and the speciation process, and we anticipate that hybrid zones will prove to be far more mobile than currently appreciated.

Postglacial recolonization shaped the genetic diversity of the winter moth (Operophtera brumata) in Europe

Changes in climate conditions, particularly during the Quaternary climatic oscillations, have long been recognized to be important for shaping patterns of species diversity. For species residing in the western Palearctic, two commonly observed genetic patterns resulting from these cycles are as follows: (1) that the numbers and distributions of genetic lineages correspond with the use of geographically distinct glacial refugia and (2) that southern populations are generally more diverse than northern populations (the “southern richness, northern purity” paradigm). To determine whether these patterns hold true for the widespread pest species the winter moth (Operophtera brumata), we genotyped 699 individual winter moths collected from 15 Eurasian countries with 24 polymorphic microsatellite loci. We find strong evidence for the presence of two major genetic clusters that diverged ~18 to ~22 ka, with evidence that secondary contact (i.e., hybridization) resumed ~ 5 ka along a well‐established hybrid zone in Central Europe. This pattern supports the hypothesis that contemporary populations descend from populations that resided in distinct glacial refugia. However, unlike many previous studies of postglacial recolonization, we found no evidence for the “southern richness, northern purity” paradigm. We also find evidence for ongoing gene flow between populations in adjacent Eurasian countries, suggesting that long‐distance dispersal plays an important part in shaping winter moth genetic diversity. In addition, we find that this gene flow is predominantly in a west‐to‐east direction, suggesting that recently debated reports of cyclical outbreaks of winter moth spreading from east to west across Europe are not the result of dispersal.

Stabilization of a salamander moving hybrid zone

When related species meet upon postglacial range expansion, hybrid zones are frequently formed. Theory predicts that such zones may move over the landscape until equilibrium conditions are reached. One hybrid zone observed to be moving in historical times (1950–1979) is that of the pond‐breeding salamanders Triturus cristatus and Triturus marmoratus in western France. We identified the ecological correlates of the species hybrid zone as elevation, forestation, and hedgerows favoring the more terrestrial T. marmoratus and pond density favoring the more aquatic T. cristatus. The past movement of the zone of ca. 30 km over three decades has probably been driven by the drastic postwar reduction of the “bocage” hedgerow landscape, favoring T. cristatus over T. marmoratus. No further hybrid zone movement was observed from 1979 to the present. To explain the changing dynamics of the hybrid zone, we propose that it stalled, either because an equilibrium was found at an altitude of ca. 140 m a.s.l. or due to pond loss and decreased population densities. While we cannot rule out the former explanation, we found support for the latter. Under agricultural intensification, ponds in the study area are lost at an unprecedented rate of 5.5% per year, so that remaining Triturus populations are increasingly isolated, hampering dispersal and further hybrid zone movement.

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

Hotspots of intraspecific genetic diversity, which are of primary importance for the conservation of species, have been associated with glacial refugia, that is areas where species survived the Quaternary climatic oscillations. However, the proximate mechanisms generating these hotspots remain an open issue. Hotspots may reflect the long-term persistence of large refugial populations; alternatively, they may result from allopatric differentiation between small and isolated populations, that later admixed. Here, we test these two scenarios in a widely distributed species of tree frog, Hyla orientalis, which inhabits Asia Minor and southeastern Europe. We apply a fine-scale phylogeographic survey, combining fast-evolving mitochondrial and nuclear markers, with a dense sampling throughout the range, as well as ecological niche modelling, to understand what shaped the genetic variation of this species. We documented an important diversity centre around the Black Sea, composed of multiple allopatric and/or parapatric diversifications, likely driven by a combination of Pleistocene climatic fluctuations and complex regional topography. Remarkably, this diversification forms a ring around the Black Sea, from the Caucasus through Anatolia and eastern Europe, with terminal forms coming into contact and partially admixing in Crimea. Our results support the view that glacial refugia generate rather than host genetic diversity and can also function as evolutionary melting pots of biodiversity. Moreover, we report a new case of ring diversification, triggered by a large, yet cohesive dispersal barrier, a very rare situation in nature. Finally, we emphasize the Black Sea region as an important centre of intraspecific diversity in the Palearctic with implications for conservation.

Is the Danube crested newt Triturus dobrogicus polytypic? A review and new nuclear DNA data

The Danube crested newtTriturus dobrogicushas been proposed to comprise two subspecies:T. d. dobrogicusandT. d. macrosoma. Uncertainty exists in the literature over their distribution and diagnosability. We conduct a multilocus phylogeographical survey and review published data to determine whether a two taxon treatment is warranted. Newly produced and published nuclear DNA data suggest intraspecific variation in the Pannonian Plain part of the range, but with extensive genetic admixture, whereas mitochondrial DNA data shows a lack of geographical structuring inT. dobrogicusaltogether. None of the studied morphological characters suggest the presence of two geographical groups inT. dobrogicusunequivocally. Although Danube Delta newts do have relatively short bodies compared to the remainder of the range (the Pannonian and Lower Danube Plains and the Dnepr Delta), we argue that this finding can be explained by phenotypic plasticity – particularly in light of the incongruent evolutionary scenario suggested by genetic data. We conclude that the total body of evidence does not support the two subspecies hypothesis and recommend thatT. dobrogicusis treated as a monotypic species.

Effect of biogeographic history on population vulnerability in European amphibians

The genetic diversity of populations, which contributes greatly to their adaptive potential, is negatively affected by anthropogenic habitat fragmentation and destruction. However, continental-scale losses of genetic diversity also resulted from the population expansions that followed the end of the last glaciation, an element that is rarely considered in a conservation context. We addressed this issue in a meta-analysis in which we compared the spatial patterns of vulnerability of 18 widespread European amphibians in light of phylogeographic histories (glacial refugia and postglacial routes) and anthropogenic disturbances. Conservation statuses significantly worsened with distances from refugia, particularly in the context of industrial agriculture; human population density also had a negative effect. These findings suggest that features associated with the loss of genetic diversity in post-glacial amphibian populations (such as enhanced fixation load or depressed adaptive potential) may increase their susceptibility to current threats (e.g., habitat fragmentation and pesticide use). We propose that the phylogeographic status of populations (i.e., refugial vs. post-glacial) should be considered in conservation assessments for regional and national red lists.