Phylogeographic analysis reveals northerly refugia for the riverine amphibianTriturus dobrogicus(Caudata: Salamandridae)

Detailed DNA barcoding of mayflies in a small European country proved how far we are from having comprehensive barcode reference libraries

Mayflies (Ephemeroptera) are among the crucial water and habitat quality bioindicators. However, despite their intensive long-term use in various studies, more reliable mayfly DNA barcode data have been produced in a negligible number of countries, and only ~40% of European species had been barcoded with less than 50% of families covered. Despite being carried out in a small area, our study presents the second-most species-rich DNA reference library of mayflies from Europe and the first comprehensive view from an important biodiversity hotspot such as the Western Carpathians. Within 1153 sequences, 76 morphologically determined species were recorded and added to the Barcode of Life Data System (BOLD) database. All obtained sequences were assigned to 97 BINs, 11 of which were unique and three represented species never barcoded before. Sequences of 16 species with high intraspecific variability were divided into 40 BINs, confirming the presence of cryptic lineages. Due to the low interspecific divergence and the non-existing barcoding gap, sequences of six species were assigned to three shared BINs. Delimitation analyses resulted in 79 and 107 putative species respectively. Bayesian and maximum-likelihood phylogenies confirmed the monophyly of almost all species and complexes of cryptic taxa and proved that DNA barcoding distinguishes almost all studied mayfly species. We have shown that it is still sufficient to thoroughly investigate the fauna of a small but geographically important area to enrich global databases greatly. In particular, the insights gained here transcend the local context and may have broader implications for advancing barcoding efforts.

Genomic analysis reveals complex population structure within the smooth newt, Lissotriton vulgaris, in Central Europe

Species with wide-range distributions usually display high genetic variation. This variation can be partly explained by historical lineages that were temporally isolated from each other and are back into secondary reproductive contact, and partly by local adaptations. The smooth newt (Lissotriton vulgaris) is one of the most widely distributed amphibians species across Eurasia and forms a species complex with a partially overlapping distribution and morphology. In the present study, we explored the population genomic structure of smooth newt lineages in the Carpathian Basin (CB) relying on single-nucleotide polymorphisms. Our dataset included new and previously published data to study the secondary contact zone between lineages in the CB and also tested for the barrier effect of rivers to gene flow between these lineages. We confirmed the presence of the South L. v. vulgaris Lineage distributed in Transdanubia and we provided new distribution records of L. v. ampelensis inhabiting the eastern territories of the CB. High genetic diversity of smooth newts was observed, especially in the North Hungarian Mountains and at the interfluves of the main rivers in the South with four distinct lineages of L. v. vulgaris and one lineage of L. v. ampelensis showing a low level of admixture with the spatially closest L. v. vulgaris lineage. Moreover, admixture detected at the interfluve of the main rivers (i.e. Danube and Tisza) suggested a secondary contact zone in the area. Finally, we found that the river Danube has a very weak effect on population divergence, while the river Tisza is a geographical barrier limiting gene flow between smooth newt lineages. As the range boundaries of L. v. vulgaris and L. v. ampelensis in the CB coincide with the river Tisza, our study underpins the influence of rivers in lineage diversification.

Strong intraspecific phylogenetic and karyotypic diversification in Isophya modestior (Orthoptera: Tettigoniidae: Phaneropterinae)

Isophya modestior (Orthoptera: Tettigoniidae) is a species distributed in central and south-eastern Europe, where its distribution is largely separated by two large rivers (Sava and Danube). Since previous studies on the song and morphology of the stridulatory file across its complete geographic range showed that the species is separated into two main groups, we decided to use phylogenetic and karyological analyses in order to evaluate the status of the previously analysed populations.

Phylogenetic analyses showed the existence of two major clades within I. modestior with very high bootstrap values and posterior probabilities—Clade A: present on the Balkan Peninsula, Slovenia (Inner Carniola), Italy, Pannonian Serbia (Vršac Mts and Deronje) and Austria (Burgenland and Lower Austria); Clade B: present in Slovenia (Upper Carniola), Croatia and Austria (Carinthia), Pannonian Serbia (Fruška Gora Mt.) and Hungary.

A comparison of chromosomes of 51 specimens revealed discrete differences between their karyotypes. The physical characteristics of the karyotypes included chromosome number (2n), sex chromosome (X) morphology and C-banding patterns. The standard chromosome complement of 50 specimens from different localities is characterized by 2n = 30 + X0 in males. In one male collected in Fruška Gora Mt. (Andrevlje), the chromosome number was reduced to 2n = 28 + neo-XY. Therefore, further cytogenetic studies involving larger samples, especially from Fruška Gora Mt., are needed in order to gain a more comprehensive view of the chromosome evolution in this group of Isophya species.

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.

Caddisfly (Trichoptera, Insecta) fauna and assemblages of the north-eastern part of the Pannonian Lowland (West Ukraine, Transcarpathia)

The caddisfly fauna of the Transcarpathian part of the Pannonian Lowland was poorly studied formerly. Here, we present the results of a six-year survey (2015-2020) carried out in four sampling sites of the Ukrainian part of the Bereg Plain and provide the actualised checklist of this area. Actually, 7346 specimens of 53 caddisfly species were collected. The number of known caddisfly species increased from 13 to 61. Two species Hydropsycheguttata and Parasetodesrespersellus, which formerly were considered extinct in the Pannonian Ecoregion, were detected and another especially rare species (e.g. Cyrnusflavidus) was also recorded. The fauna of the region cover a significant part of both Hungarian and Ukrainian caddisfly fauna. Assemblages of four characteristic habitat types of the region showed significant differences considering their quantitative and qualitative composition, substrate, current, hydrological- and feeding types. The high diversity and natural value of the small lowland watercourses were proven using a new Caddisfly Conservation Index (CCI) calculated, based on vulnerability and rarity of species. The fauna and assemblages showed a unique character mainly independent from large rivers of the region.

Climate, immigration and speciation shape terrestrial and aquatic biodiversity in the European Alps

Quaternary climate fluctuations can affect speciation in regional biodiversity assembly in two non-mutually exclusive ways: a glacial species pump, where isolation in glacial refugia accelerates allopatric speciation, and adaptive radiation in underused adaptive zones during ice-free periods. We detected biogeographic and genetic signatures associated with both mechanisms in the assembly of the biota of the European Alps. Age distributions of endemic and widespread species within aquatic and terrestrial taxa (amphipods, fishes, amphibians, butterflies and flowering plants) revealed that endemic fish evolved only in lakes, are highly sympatric, and mainly of Holocene age, consistent with adaptive radiation. Endemic amphipods are ancient, suggesting preglacial radiation with limited range expansion and local Pleistocene survival, perhaps facilitated by a groundwater-dwelling lifestyle. Terrestrial endemics are mostly of Pleistocene age and are thus more consistent with the glacial species pump. The lack of evidence for Holocene adaptive radiation in the terrestrial biome is consistent with faster recolonization through range expansion of these taxa after glacial retreats. More stable and less seasonal ecological conditions in lakes during the Holocene may also have contributed to Holocene speciation in lakes. The high proportion of young, endemic species makes the Alpine biota vulnerable to climate change, but the mechanisms and consequences of species loss will likely differ between biomes because of their distinct evolutionary histories.

Genetic traces of hybrid zone movement across a fragmented habitat

Theoretical and empirical studies suggest that the structure and position of hybrid zones can change over time. Evidence for moving hybrid zones has been directly inferred by repeated sampling over time, or indirectly through the detection of genetic footprints left by the receding species and the resulting asymmetric patterns of introgression across markers. We here investigate a hybrid zone formed by two subspecies of the Iberian golden-striped salamander, Chioglossa lusitanica, using a panel of 35 nuclear loci (31 SNPs and 4 allozymes) and one mitochondrial locus in a transect in central Portugal. We found concordant and coincident clines for most of the nuclear loci (n=22, 63%), defining a narrow hybrid zone of ca. 6 km wide, with the centre positioned ca. 15 km south of the Mondego river. Asymmetric introgression was observed at another 14 loci. Their clines are displaced towards the north, with positions located either close to the Mondego river (n=6), or further northwards (n=8). We interpret these profiles as genetic traces of the southward displacement of C. lusitanica lusitanica by C. l. longipes over the wider Mondego river valley. We noted the absence of significant linkage disequilibrium and we inferred low levels of effective selection per locus against hybrids, suggesting that introgression in the area of species replacement occurred under a neutral diffusion process. A species distribution model suggests that the C. lusitanica hybrid zone coincides with a narrow corridor of fragmented habitat. From the position of the displaced clines, we infer that patches of locally suitable habitat trapped some genetic variants that became disassociated from the southward moving hybrid zone. This study highlights the influence of habitat availability on hybrid zone movement.

DNA barcodes evidence the contact zone of eastern and western caddisfly lineages in the Western Carpathians

The region of the Western Carpathians is, among other aspects, very important for survival and diversity of European freshwater fauna due to the presence of a large number of (sub)mountain springs and streams. However, these ecologically and faunistically diversified habitats are still understudied in the context of genetic diversity and population structure of their inhabitants. This study focuses on genetic diversity and distribution patterns of the caddisfly Rhyacophila tristis, common and widespread representative of mountain freshwater fauna. Analysis of the COI mitochondrial marker revealed presence of the western and eastern lineages, with samples from both lineages being grouped in BOLD (Barcode of Life Data System) into separate BINs (Barcode Index Numbers). Our data indicates that eastern lineage (BIN_E) is more closely related to the Balkan populations than to co-occurring western lineage (BIN_W), and that the contact zone of the lineages passes through the W Carpathians. The study revealed phylogeographic and demographic differences between lineages, supporting hypothesis of their evolutionary independence and specific ecological preferences. The obtained genetic data of the R. tristis population from W Carpathians improved our knowledge about population genetics of this aquatic species and can contribute to understanding the state and evolution of biodiversity of freshwater ecosystems in Europe.

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.

Mitochondrial DNA diversity of the alpine newt (Ichthyosaura alpestris) in the Carpathian Basin: evidence for multiple cryptic lineages associated with Pleistocene refugia

The phylogeography and molecular taxonomy of the Alpine newt, Ichthyosaura alpestris, has been intensively studied in the past. However, previous studies did not include a comprehensive sampling from the Carpathian Basin, possibly a key region in the evolution of the species. We used a 1251 bp long fragment of the mitochondrial genome to infer the species’ evolutionary history in central-eastern Europe by assigning isolated Carpathian Basin populations from 6 regions to previously defined mtDNA lineages. We also revised the morphology-based intraspecific taxonomy of the species in the light of new genetic data. Alpine newt populations from the Carpathian Basin represented two different mitochondrial lineages. The Mátra, Bükk and Zemplén Mts populations can be assigned to the Western lineage of the nominotypical subspecies. Bakony and Őrség populations showed high haplotype diversity and formed a separate clade within the Western lineage, suggesting that the Carpathian Basin might have provided cryptic refugia for Alpine newt populations in their cold-continental forest-steppe landscapes during the younger Pleistocene. Newts from Apuseni Mts were related to the Eastern lineage but formed a distinct clade within this lineage. Considering the morphological and genetic differentiation of the Bakony and Őrség populations, consistent with a long independent evolutionary history, we propose these populations be referred to as Ichthyosaura alpestris bakonyiensis (Dely, 1964). We provide a redescription of this poorly known subspecies.

The tale of springs and streams: how different aquatic ecosystems impacted the mtDNA population structure of two riffle beetles in the Western Carpathians

The Western Carpathians are a particularly interesting part of the Carpathian Arc. According to recent molecular data upon aquatic and terrestrial taxa, this mountain area is an important biodiversity hotspot of Europe. Moreover, the W Carpathians include rich systems of karst springs inhabited by specific fauna, where molecular diversity and phylogeographic patterns are yet to be fully explored. Our study aims to compare population genetic structure and molecular diversity of two related and commonly co-occurring riffle beetles, Elmis aenea (PWJ Müller, 1806) and Limnius perrisi (Dufour, 1843) in the springs and streams of the W Carpathians using the mitochondrial DNA barcoding fragment of the cytochrome c oxidase subunit I gene (COI). The relatively stable thermal and chemical conditions of springs throughout unfavourable climatic settings make these highly specific lotic systems potentially ideal for a long-term survival of some aquatic biota. Populations of both elmid species were relatively homogeneous genetically, with a single dominant haplotype. However, we revealed that E. aenea significantly dominated in the springs, while L. perrisi preferred streams. Relative isolation of the springs and their stable conditions were reflected in significantly higher molecular diversity of the E. aenea population in comparison to L. perrisi. The results of Bayesian Skyline Plot analysis also indicated the exceptional position of springs regarding maintaining the population size of E. aenea. On the other hand, it seems that streams in the W Carpathians provide more effective dispersal channels for L. perrisi, whose population expanded much earlier compared to E. aenea. Present study points out that different demographic histories of these two closely related elmid species are manifested by their different habitat preference and molecular diversity.

Limited long-distance dispersal success in a Western European fairy shrimp evidenced by nuclear and mitochondrial lineage structuring

Anostraca are known by their ability for long-distance dispersal, but the existence in several species of deep, geographically structured mtDNA lineages suggests their populations are subjected to allopatric differentiation, isolation, and prevalence of local scale dispersion. Tanymastix stagnalis is one of the most widespread species of Anostraca and previous studies revealed an unclear geographical pattern of mtDNA genetic diversity. Here, we analyze populations from the Iberian and Italian Peninsulas, Central Europe, and Scandinavia, with the aim to characterize the patterns of genetic diversity in a spatio-temporal framework using mtDNA and nuclear markers to test gene flow among close populations. For these aims we built a time-calibrated phylogeny and carried out Bayesian phylogeographic analyses using a continuous diffusion model. Our results indicated that T. stagnalis presents a deeply structured genetic diversity, including 7 ancient lineages, some of them even predating the Pleistocene. The Iberian Peninsula harbors high diversity of lineages, with strong isolation and recent absence of gene flow between populations. Dispersal at local scale seems to be the prevailing dispersal mode of T. stagnalis, which exhibits a pattern of isolation-by-distance in the Iberian Peninsula. We remark the vulnerability of most of these lineages, given the limited known geographic distribution of some of them, and the high risk of losing important evolutionary potential for the species.

Species-specific habitat preferences do not shape the structure of a crested newt hybrid zone (Triturus cristatus x T. carnifex)

Reproductive isolation barriers maintain the integrity of species by preventing interspecific gene flow. They involve temporal, habitat or behavioral isolation acting before fertilization, and postzygotic isolation manifested as hybrid mortality or sterility. One of the approaches of how to study reproductive isolation barriers is through the analysis of hybrid zones. In this paper, we describe the structure of a hybrid zone between two crested newt species (Triturus cristatus and T. carnifex) in the southern part of the Czech Republic using morphological, microsatellite, and mitochondrial (mtDNA) markers. Specifically, we tested the hypothesis that the structure of the hybrid zone is maintained by species-specific habitat preferences. Comparing the genetic structure of populations with geographical and ecological parameters, we found that the hybrid zone was structured primarily geographically, with T. cristatus-like populations occurring in the northeast and T. carnifex-like populations in the southwest. Despite T. cristatus tending to occur in deeper ponds and T. carnifex on localities with more shading, the effect of both ecological parameters on the structure of the zone was minimal. Next, we corroborated that T. carnifex individuals and some hybrids possess mtDNA of T. dobrogicus, whose nuclear background was not detected in the studied hybrid zone. Hybridization between T. carnifex and T. dobrogicus (resulting in unidirectional mtDNA introgression) had to predate subsequent formation of the hybrid zone between T. cristatus and T. carnifex. Populations of crested newts in the southern part of the Czech Republic thus represent a genetic mosaic of nuclear and mitochondrial genomes of three species.

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.

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.

Hybrid zone formation and contrasting outcomes of secondary contact over transects in common toads

Much progress in speciation research stems from documenting patterns of morphological and genetic variation in hybrid zones. Contrasting patterns of marker introgression in different sections of the contact can provide valuable insights on the relative importance of various evolutionary mechanisms maintaining species differences in the face of hybridization and gene flow and on hybrid zone temporal and spatial dynamics. We studied species interactions in the common toads Bufo bufo and B. spinosus in France and northwestern Italy using morphological and molecular data from the mitochondrial and nuclear genomes in an extensive survey, including two independent transects west and east of the Alps. At both, we found sharp, coincident and concordant nuclear genetic transitions. However, morphological clines were wider or absent and mtDNA introgression was asymmetric. We discuss alternative, nonexclusive hypotheses about evolutionary processes generating these patterns, including drift, selection, long-distance dispersal and spatial shifts in hybrid zone location and structure. The distribution of intraspecific mtDNA lineages supports a scenario in which B. bufo held a local refugium during the last glacial maximum. Present-day genetic profiles are best explained by an advance of B. spinosus from a nearby Iberian refugium, largely superseding the local B. bufo population, followed by an advance of B. bufo from the Balkans, with prongs north and south of the Alps, driving B. spinosus southwards. A pendulum moving hybrid zone, first northwards and then southwards, explains the wide areas of introgression at either side of the current position of the contact zones.

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.