Phylogeography of two European newt species--discordance between mtDNA and morphology

Morphological variation, modularity and integration in the scapula and humerus of Lissotriton newts

The modular organization of tetrapod paired limbs and girdles, influenced by the expression of Hox genes is one of the primary driving forces of the evolution of animal locomotion. The increased morphological diversification of the paired limbs is correlated with reduced between-limb covariation, while correlation within the elements is usually higher than between the elements. The tailed amphibians, such as Lissotriton newts, have a biphasic lifestyle with aquatic and terrestrial environments imposing different constraints on limb skeleton. By employing the methods of computerized microtomography and 3D geometric morphometrics, we explored the pattern of morphological variation, disparity, modularity and morphological integration in the proximal parts of the anterior limbs of six species of Eurasian small bodied newts. Although the species significantly differ in limb shape, there is a great overlap in morphology of scapula and humerus, and there are no differences in morphological disparity. For the scapula, the shape differences related to the duration of the aquatic period are in length, depth and curvature. The shape of the humerus is not affected by the length of aquatic period, and shape differences between the species are related to robustness of the body. The length of aquatic period has statistically supported phylogenetic signal. The scapula and humerus are structures of varying modularity. For the humerus, the strongest support on the phylogenetic level was for the capitulum/shaft hypothesis, which can also be interpreted as functional modularity. For the scapula, the greatest support was for the antero-posterior hypothesis of modularity in case of Lissotriton vulgaris, which can be explained by different functional roles and muscle insertion patterns, while there was no phylogenetic modularity. The modularity patterns seem to correspond with the general tetrapod pattern, with modularity being more pronounced in the distal structure. The future research should include more salamandrid taxa with different habitat preferences and both adult and larval stages, in order to explore how size, phylogeny and ecology affect the morphology and covariation patterns of limbs.

Limited evidence for genetic differentiation or adaptation in two amphibian species across replicated rural-urban gradients

Urbanization leads to complex environmental changes and poses multiple challenges to organisms. Amphibians are highly susceptible to the effects of urbanization, with land use conversion, habitat destruction, and degradation ranked as the most significant threats. Consequently, amphibians are declining in urban areas, in both population numbers and abundance, however, the effect of urbanization on population genetic parameters remains unclear. Here, we studied the genomic response to urbanization in two widespread European species, the common toad Bufo bufo (26 localities, 480 individuals), and the smooth newt Lissotriton vulgaris (30 localities, 516 individuals) in three geographic regions: southern and northern Poland and southern Norway. We assessed genome-wide SNP variation using RADseq (ca. 42 and 552 thousand SNPs in toads and newts, respectively) and adaptively relevant major histocompatibility complex (MHC) class I and II genes. The results linked most of the genetic differentiation in both marker types to regional (latitudinal) effects, which also correspond to historical biogeography. Further, we did not find any association between genetic differentiation and level of urbanization at local scales for either species. However, urban smooth newts, but not toads, have lower levels of within-population genome-wide diversity, suggesting higher susceptibility to the negative effects of urbanization. A decreasing level of genetic diversity linked to increasing urbanization was also found for MHC II in smooth newts, while the relationship between MHC class I diversity and urbanization differed between geographic regions. We did not find any effects of urbanization on MHC diversity in the toad populations. Although two genetic environment association analyses of genome-wide data, LFMM and BayPass, revealed numerous (219 in B. bufo and 7040 in L. vulgaris) SNPs statistically associated with urbanization, we found a marked lack of repeatability between geographic regions, suggesting a complex and multifaceted response to natural selection elicited by life in the city.

Ventral and lateral spot patterns differentiation between three smooth newt species (Amphibia: Salamandridae: Lissotriton)

The location, shape and area of dark spots on the belly, throat and sides of the body of three closely related species of smooth newts (Lissotriton kosswigi, L. schmidtleri and L. vulgaris) were analysed. Differences were found between species and sexes in the location and shape of dark spots on the belly, throat and lateral surfaces of the body. Discriminant analysis of these spot characteristics allows to identificate of males of all three species with a high degree of confidence (85–91%). In females, only L. vulgaris was correctly distinguished from the other two species (accuracy 81–94%). Anatolian and Thracian populations of L. schmidtleri have very similar patterns of dark spots, which confirm their conspecificity. The differences in the location and size of dark spots make possible successful identification of the species during field research and study of museum specimens. The developed method could be useful for comparative studies of other animal species which have spotted camouflage or aposematic colouration.

Variation in vertebrae shape across small-bodied newts reveals functional and developmental constraints acting upon the trunk region

The salamander vertebral column is largely undifferentiated with a series of more or less uniform rib-bearing presacral vertebrae traditionally designated as the trunk region. We explored regionalization of the salamander trunk in seven species and two subspecies of the salamander genus Lissotriton by the combination of microcomputed tomography scanning and geometric morphometrics. The detailed information on trunk vertebral shape was subjected to a multidimensional cluster analysis and a phenotypic trajectory analysis. With these complementary approaches, we observed a clear morphological regionalization. Clustering analysis showed that the anterior trunk vertebrae (T1 and T2) have distinct morphologies that are shared by all taxa, whereas the subsequent, more posterior vertebrae show significant disparity between species. The phenotypic trajectory analysis revealed that all taxa share a common pattern and amount of shape change along the trunk region. Altogether, our results support the hypothesis of a conserved anterior-posterior developmental patterning which can be associated with different functional demands, reflecting (sub)species' and, possibly, regional ecological divergences within species.

Phylogenetic analyses of Eurasian lynx (Lynx lynx Linnaeus, 1758) including new mitochondrial DNA sequences from Iran

The Eurasian lynx (Lynx lynx) is one of the widespread felids in Eurasia; however, relatively little is known about the Asian subspecies, and especially the Iranian populations, which comprise the most southwestern part of its range. The current study aimed to assess the phylogenetic status of Iranian populations relative to other populations of Eurasia, by sequencing a 613 bp fragment of the mitochondrial control region. In total, 44 haplotypes were recorded from 83 sequences throughout Eurasia, two of which were found in Iran. The haplotype (H1) is dominant in all Iranian lynx populations and identical to specimens from SW Russia and central China. The second haplotype (H2) is unique and was recorded only from Ghazvin Province in the central Alborz Mountains. Both haplotypes occur in Ghazvin Province. The phylogenetic tree and a median-joining network identified four clades (i.e., East, West 1, West 2, and South). These results are congruent with previous studies and suggest that Eurasian lynx was restricted to the southern part of its range during the glacial maxima and expanded from there to East Asia and to Europe during several independent re-colonization events. The Caucasus region most like plays an important role as a refugium during glacial cycles.

Reticulate evolution in Conidae: Evidence of nuclear and mitochondrial introgression

Conidae is a hyperdiverse family of marine snails that has many hallmarks of adaptive radiation. Hybridization and introgression may contribute to such instances of rapid diversification by generating novel gene combinations that facilitate exploitation of distinct niches. Here we evaluated whether or not these mechanisms may have contributed to the evolutionary history of a subgenus of Conidae (Virroconus). Several observations hint at evidence of past introgression for members of this group, including incongruence between phylogenetic relationships inferred from mitochondrial gene sequences and morphology and widespread sympatry of many Virroconus species in the Indo-West Pacific. We generated and analyzed transcriptome data of Virroconus species to (i) infer a robust nuclear phylogeny, (ii) assess mitochondrial and nuclear gene tree discordance, and (iii) formally test for introgression of nuclear loci. We identified introgression of mitochondrial genomes and nuclear gene regions between ancestors of one pair of Virroconus species, and mitochondrial introgression between another pair. We also found evidence of adaptive introgression of conotoxin venom loci between a third pair of species. Together, our results demonstrate that hybridization and introgression impacted the evolutionary history of Virroconus and hence may have contributed to the adaptive radiation of Conidae.

Settling taxonomic and nomenclatural problems in brine shrimps, Artemia (Crustacea: Branchiopoda: Anostraca), by integrating mitogenomics, marker discordances and nomenclature rules

High morphological plasticity in populations of brine shrimp subjected to different environmental conditions, mainly salinity, hindered for centuries the identification of the taxonomic entities encompassed within Artemia. In addition, the mismatch between molecular and morphological evolution rates complicates the characterization of evolutionary lineages, generating taxonomic problems. Here, we propose a phylogenetic hypothesis for Artemia based on two new complete mitogenomes, and determine levels of congruence in the definition of evolutionary units using nuclear and mtDNA data. We used a fossil of Artemia to calibrate the molecular clock and discuss divergence times within the genus. The hypothesis proposed herein suggests a more recent time frame for lineage splitting than previously considered. Phylogeographic analyses were performed using GenBank available mitochondrial and nuclear markers. Evidence of gen e flow, identified through discordances between nuclear and mtDNA markers, was used to reconsider the specific status of some taxa. As a result, we consider Artemia to be represented by five evolutionary units: Southern Cone, Mediterranean-South African, New World, Western Asian, and Eastern Asian Lineages. After an exhaustive bibliographical revision, unavailable names for nomenclatural purposes were discarded. The remaining available names have been assigned to their respective evolutionary lineage. The proper names for the evolutionary units in which brine shrimps are structured remain as follows: Artemia persimilis Piccinelli & Prosdocimi, 1968 for the Southern Cone Lineage, Artemia salina (Linnaeus, 1758) for the Mediterranean-SouthAfrican Lineage, Artemia urmiana Günther, 1899 for the Western Asian Lineage, and Artemia sinica Cai, 1989 for the Eastern Asian Lineage. The name Artemia monica Verrill, 1869 has nomenclatural priority over A. franciscana Kellogg, 1906 for naming the New World Lineage. New synonymies are proposed for A. salina (= C. dybowskii Grochowski, 1896 n. syn., and A. tunisiana Bowen & Sterling, 1978 n. syn.), A. monica (= A. franciscana Kellogg, 1906 n. syn., and A. salina var. pacifica Sars, 1904 n. syn.); A. urmiana (= B. milhausenii Fischer de Waldheim, 1834 n. syn., A. koeppeniana Fischer, 1851 n. syn., A. proxima King, 1855 n. syn., A. s. var. biloba Entz, 1886 n. syn., A. s. var. furcata Entz, 1886 n. syn., A. asiatica Walter, 1887 n. syn., A. parthenogenetica Bowen & Sterling, 1978 n. syn., A. ebinurica Qian & Wang, 1992 n. syn., A. murae Naganawa, 2017 n. syn., and A. frameshifta Naganawa & Mura, 2017 n. syn.). Internal deep nuclear structuring within the A. monica and A. salina clades, might suggest the existence of additional evolutionary units within these taxa.

Genetic diversity, phylogenetic position and morphometric analysis of Astacus colchicus (Decapoda, Astacidae): a new insight into Eastern European crayfish fauna

The phylogeny of European crayfish fauna, especially with respect to Eastern European species, is still far from being completely resolved. To fill this gap, we analyzed most of the European crayfish species focusing on the phylogenetic position of the endemic crayfish Astacus colchicus, inhabiting Georgia. Three mitochondrial and one nuclear marker were used to study evolutionary relationships among European crayfish species, resulting in the unique phylogenetic position of A. colchicus indicating independent species status to A. astacus. Phylogenetic analyses revealed a deep molecular divergence of A. colchicus in comparison to A. astacus (6.5-10.9% in mtDNA and 1.1% in nDNA) as well as to Pontastacus leptodactylus and P. pachypus (5.5-10.0% in mtDNA and 1.4-2.4% in nDNA). Absent ventral process on second male pleopod and abdominal somites II and III with pleura rounded lacking prominent spines clearly indicate taxonomic assignment to the genus Astacus; however, the species is distributed almost in the middle of Ponto-Caspian area typical by occurrence of the genus Pontastacus. Several morphological indices linked to head length, carapace, and total body length and width were found to demonstrate apparent differences between A. colchicus and A. astacus. Although this study provides a novel insight into European crayfish phylogeography, we also point out the gaps in comprehensive study of the P. leptodactylus species complex, which could reveal details about the potential species status of particular species and subspecies within this genus.

The impact of Miocene orogeny for the diversification of Caucasian Epeorus (Caucasiron) mayflies (Ephemeroptera: Heptageniidae)

A common hypothesis for the high biodiversity of mountains is the diversification driven by orogeny creating conditions for rapid in situ speciation of resident lineages. The Caucasus is a young mountain system considered as a biodiversity hotspot; however, the origin and evolution of its diversity remain poorly understood. This study focuses on mayflies of the subgenus Caucasiron, one of the most diversified stenotopic mayflies inhabiting various types of streams throughout the Caucasus. Using the time-calibrated phylogeny based on two mitochondrial (COI, 16S) and three nuclear (EF-1α, wg, 28S) gene fragments, we tested the role of Caucasian orogeny in biogeography, diversification patterns, and altitudinal diversification of Caucasiron mayflies. We found that orogeny promoted the lineage diversification of Caucasiron in the Miocene. The highest diversification rate corresponding with the uplift of mountains was followed by a significant slowdown towards the present suggesting minor influence of Pleistocene climatic oscillations on the speciation. The Caucasiron lineages cluster into three principal clades originating in the Upper Miocene. We found a strong support that one of the three clades diversified via allopatric speciation in the Greater Caucasus isolated in the Parathetys Sea. The other two clades originating most likely outside the Greater Caucasus diversified towards high and low altitude, respectively, indicating possible role of climatic factors and/or passive uplift on their differentiation. Current high Caucasiron diversity in the Greater Caucasus is a result of in situ speciation and later immigration from adjacent mountain ranges after the Parathetys Sea retreat. Our phylogeny supported the monophyly of Rhithrogeninae, Epeorus s.l., Caucasiron, and Iron. Epeorus subgenus Ironopsis was found paraphyletic, with its European representatives more closely related to Epeorus s.str. than to Iron. Therefore, we re-arranged taxa treated within Ironopsis to comply with the phylogeny recovered herein.

Mitogenomics of historical type specimens of Australasian turtles: clarification of taxonomic confusion and old mitochondrial introgression

Diagnosability is central to taxonomy as are type specimens which define taxa. New advances in technologies and the discovery of new informative traits must be matched with previous taxonomic decisions based on name-bearing type specimens. Consequently, the challenge of sequencing highly degraded DNA from historical types becomes an inevitability to resolve the very many taxonomic issues arising from, by modern standards, poor historical species descriptions leading to difficulties to assign names to genetic clusters identified from fresh material. Here we apply high-throughput parallel sequencing and sequence baiting to reconstruct the mitogenomes from 18 type specimens of Australasian side-necked turtles (Chelidae). We resolve a number of important issues that have confused the taxonomy of this family, and analyse the mitogenomes of the types and those of fresh material to improve our understanding of the phylogenetic relationships of this morphologically conservative group. Together with previously published nuclear genomic data, our study provides evidence for multiple old mitochondrial introgressions.

Herps without borders: a new newt case and a review of transalpine alien introductions in western Europe

Biogeographic processes have led to different evolutionary taxa occurring in the northern and southern edges of the Alpine Mountains in Western Europe. The integrity of this diversity is being challenged by frequent human-mediated trans-alpine translocations, sometimes leading to biological invasions. Several alien terrestrial vertebrates of south Alpine origins (Italy, Swiss Ticino) are causing damages to native north Alpine fauna. In this paper, we used molecular tools to characterize the understudied case of the Mediterranean smooth newt (Lissotriton vulgaris meridionalis) expanding in the outskirts of Geneva since its introduction before 1975. Mitochondrial and nuclear DNA sequencing suggest that these exotic populations are a mixture between two diverged L. v. meridionalis lineages from central Italy, and traces of potential hybridization with the native L. v. vulgaris was detected. This situation echoes many other trans-alpine alien introductions. We review all comparable cases of southern to northern Alps introductions in vertebrates, including seven reptiles and four amphibians. The majority of south alpine alien lineages were presumably imported voluntarily by enthusiasts and appear to perform better in the disturbed habitats found in the anthropogenic landscapes of Western Europe compared to their native north Alpine counterparts. Most pose serious threats to related species of similar ecology, through direct competition, predation and introgressive hybridization. Difficulties to detect alien species on time lead to significant conservation costs. Better education together with more appropriate and reactive management plans will be necessary to limit the impact of future alien introductions.

Comprehensive molecular phylogeny of barn owls and relatives (Family: Tytonidae), and their six major Pleistocene radiations

The owl family Tytonidae comprises two genera: Phodilus, limited to the forests of central Africa and South-East Asia, and the ubiquitous Tyto. The genus Tyto is majorly represented by the cosmopolitan Common Barn Owl group, with more than 30 subspecies worldwide. Discrete differences in body size and plumage colouration have led to the classification of this family into many species and subspecies, but the taxonomic status and phylogenetic relationships between taxa remain unclear, and in some groups controversial. Although several previous studies attempted to resolve this problem, they have been limited in their taxonomic and geographical coverage, or have relied on restricted molecular evidence and low sample sizes. Based on the most comprehensive sampling to date (16 out of 17 Tyto species, and one out of three Phodilus species), a multi-locus approach using seven mitochondrial and two nuclear markers, and taking advantage of field data and museum collections available worldwide, our main questions in this study were: (1) what are the phylogenetic relationships and classification status of the whole family; (2) when and where did the most important speciation events occur? We confirm that the Common Barn Owl, Tyto alba is divided into three main evolutionary units: the American Barn Owl, T. furcata; the Western Barn Owl, T. alba; and the Eastern Barn Owl, T. javanica, and suggest a Late Miocene (ca. 6 mya) Australasian and African origin of the group. Our results are supported by fossil age information, given that the most recent common ancestor between the Tytonidae genera Phodilus and Tyto was probably from the Oligocene (ca. 28 mya) of Australasia. We finally reveal six major Pleistocene radiations of Tyto, all resulting in wide-range distributions.

The distributions of the six species constituting the smooth newt species complex (Lissotriton vulgaris sensu lato and L. montandoni) – an addition to the New Atlas of Amphibians and Reptiles of Europe

The ‘smooth newt’, the taxon traditionally referred to as Lissotriton vulgaris, consists of multiple morphologically distinct taxa. Given the uncertainty concerning the validity and rank of these taxa, L. vulgaris sensu lato has often been treated as a single, polytypic species. A recent study, driven by genetic data, proposed to recognize five species, L. graecus, L. kosswigi, L. lantzi, L. schmidtleri and a more restricted L. vulgaris. The Carpathian newt L. montandoni was confirmed to be a closely related sister species. We propose to refer to this collective of six Lissotriton species as the smooth newt or Lissotriton vulgaris species complex. Guided by comprehensive genomic data from throughout the range of the smooth newt species complex we 1) delineate the distribution ranges, 2) provide a distribution database, and 3) produce distribution maps according to the format of the New Atlas of Amphibians and Reptiles of Europe, for the six constituent species. This allows us to 4) highlight regions where more research is needed to determine the position of contact zones.

Geographic patterns of phenotypic diversity in incipient species of North American blister beetles (Coleoptera: Meloidae) are not determined by species niches, but driven by demography along the speciation process

The Epicauta stigmata complex is a group of blister beetles composed of three parapatric or sympatric species that occur in central Mexico to southern USA: E. stigmata, E. uniforma and E. melanochroa. These species are morphologically very similar, and are mainly distinguished by body colour differences. Here we assessed whether phenotypic divergence in coloration patterns define evolutionary units within the complex. We studied the phylogenetic relationships, demographic history and concordances between morphological and ecological traits in the group. The complex apparently had a demographic history of recent population expansion during the last glaciation period 75000 to 9500 years ago. The three species show no reciprocal monophyly, and thus their allospecificity was not confirmed. The current distribution of haplotypes and the genetic divergences in these taxa can be explained by either recent mitochondrial introgression events caused by hybridisation or by incomplete lineage sorting. Colour pattern differences in the complex are not likely a product of local selection acting over a common genetic background. We suggest that phenotypic divergence in colour patterns during an incipient speciation process might be seen as an enhancing factor of cohesion within each of the three evolutionary units.

Isolation and gene flow in a speciation continuum in newts

Because reproductive isolation often evolves gradually, differentiating lineages may retain the potential for genetic exchange for prolonged periods, providing an opportunity to quantify and to understand the fundamental role of gene flow during speciation. Here we delimit evolutionary lineages, reconstruct the phylogeny and infer gene flow in newts of the Lissotriton vulgaris species complex based on 74 nuclear markers sampled from 127 localities. We demonstrate that distinct lineages along the speciation continuum in newts exchange nontrivial amounts of genes, affecting their evolutionary trajectories. By integrating a wide array of methods, we delimit nine evolutionary lineages and show that two principal factors have driven their genetic differentiation: time since the last common ancestor determining levels of shared ancestral polymorphism, and shifts in geographic distributions determining the extent of secondary contact. Post-divergence gene flow, indicative of evolutionary non-independence, has been most extensive in Central Europe, while four southern European lineages have acquired the population-genetic hallmarks of independent species (L. graecus, L. kosswigi, L. lantzi, L. schmidtleri). We obtained strong statistical support for widespread mtDNA introgression following secondary contact, previously suggested by discordance between mtDNA phylogeny and morphology. Our study reveals long-term evolutionary persistence of evolutionary lineages that may periodically exchange genes with one another: although some of these lineages may become extinct or fuse, others will acquire complete reproductive isolation and will carry signatures of this complex history in their genomes.

Linkage Map of Lissotriton Newts Provides Insight into the Genetic Basis of Reproductive Isolation

Linkage maps are widely used to investigate structure, function, and evolution of genomes. In speciation research, maps facilitate the study of the genetic architecture of reproductive isolation by allowing identification of genomic regions underlying reduced fitness of hybrids. Here we present a linkage map for European newts of the Lissotriton vulgaris species complex, constructed using two families of F2 L. montandoni × L. vulgaris hybrids. The map consists of 1146 protein-coding genes on 12 linkage groups, equal to the haploid chromosome number, with a total length of 1484 cM (1.29 cM per marker). It is notably shorter than two other maps available for salamanders, but the differences in map length are consistent with cytogenetic estimates of the number of chiasmata per chromosomal arm. Thus, large salamander genomes do not necessarily translate into long linkage maps, as previously suggested. Consequently, salamanders are an excellent model to study evolutionary consequences of recombination rate variation in taxa with large genomes and a similar number of chromosomes. A complex pattern of transmission ratio distortion (TRD) was detected: TRD occurred mostly in one family, in one breeding season, and was clustered in two genomic segments. This is consistent with environment-dependent mortality of individuals carrying L. montandoni alleles in these two segments and suggests a role of TRD blocks in reproductive isolation. The reported linkage map will empower studies on the genomic architecture of divergence and interactions between the genomes of hybridizing newts.

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.

Extra-Mediterranean glacial refugia in a Mediterranean faunal element: the phylogeography of the chalk-hill blue Polyommatus coridon (Lepidoptera, Lycaenidae)

Most warm-adapted Central European species are thought to have survived ice ages exclusively in Mediterranean refugia. During recent years, this point of view has been questioned. Therefore, we tested the hypothesis that extra-Mediterranean refugia also played a role in warm-adapted insect species and selected the chalk-hill blue, Polyommatus coridon. We sequenced two mitochondrial loci (COI, CR) in 150 individuals from 30 populations covering nearly the complete range. Minimum spanning networks and other statistical analyses concordantly revealed four genetic lineages with strong phylogeographic signal: a western group in Italy, France and western/central Germany, an eastern lineage in the Balkan Peninsula, the Carpathian Basin and eastern Central Europe, an Alpine group with populations in the Alps and southern Germany and a Pyrenean group. Our results are generally consistent with previous analyses for P. coridon based on allozymes and DNA sequences, but provide additional insights. We propose that these four lineages have evolved during allopatry in different glacial refugia, two in typical Mediterranean refugia (Apennines and Balkan Peninsulas), but two in extra-Mediterranean areas south of the Alps and Pyrenees. This supports survival of warm-adapted organisms in these regions in close geographic proximity to the refugia of high mountain species.

Selective Landscapes in newt Immune Genes Inferred from Patterns of Nucleotide Variation

Host–pathogen interactions may result in either directional selection or in pressure for the maintenance of polymorphism at the molecular level. Hence signatures of both positive and balancing selection are expected in immune genes. Because both overall selective pressure and specific targets may differ between species, large-scale population genomic studies are useful in detecting functionally important immune genes and comparing selective landscapes between taxa. Such studies are of particular interest in amphibians, a group threatened worldwide by emerging infectious diseases. Here, we present an analysis of polymorphism and divergence of 634 immune genes in two lineages of Lissotriton newts: L. montandoni and L. vulgaris graecus. Variation in newt immune genes has been shaped predominantly by widespread purifying selection and strong evolutionary constraint, implying long-term importance of these genes for functioning of the immune system. The two evolutionary lineages differ in the overall strength of purifying selection which can partially be explained by demographic history but may also signal differences in long-term pathogen pressure. The prevalent constraint notwithstanding, 23 putative targets of positive selection and 11 putative targets of balancing selection were identified. The latter were detected by composite tests involving the demographic model and further validated in independent population samples. Putative targets of balancing selection encode proteins which may interact closely with pathogens but include also regulators of immune response. The identified candidates will be useful for testing whether genes affected by balancing selection are more prone to interspecific introgression than other genes in the genome.

The importance of molecular markers and primer design when characterizing biodiversity from environmental DNA

Environmental DNA (eDNA) comprises DNA fragments that have been shed into the environment by organisms, and which can be extracted from environmental samples such as water or soil. Characterization of eDNA can allow researchers to infer the presence or absence of species from a particular site without the need to locate and identify individuals, and therefore may provide an extremely valuable tool for quantifying biodiversity. However, as is often the case with relatively new protocols, methodological challenges remain. A number of earlier reviews have discussed these challenges, but none have provided extensive treatment of the critical decisions surrounding molecular markers and primer development for use in eDNA assays. This review discusses a number of options and approaches that can be used when determining which primers and gene regions are most appropriate for either targeted species detection or metabarcoding macro-organisms from eDNA. The latter represents a new field that is growing rapidly, and which has the potential to revolutionize future assessments of community and ecosystem diversity.

DNA barcoding of crickets, katydids and grasshoppers (Orthoptera) from Central Europe with focus on Austria, Germany and Switzerland

We present a DNA barcoding study on the insect order Orthoptera that was generated in collaboration between four barcoding projects in three countries, viz. Barcoding Fauna Bavarica (Germany), German Barcode of Life, Austrian Barcode of Life and Swiss Barcode of Life. Our data set includes 748 COI sequences from 127 of the 162 taxa (78.4%) recorded in the three countries involved. Ninety-three of these 122 species (76.2%, including all Ensifera) can be reliably identified using DNA barcodes. The remaining 26 caeliferan species (families Acrididae and Tetrigidae) form ten clusters that share barcodes among up to five species, in three cases even across different genera, and in six cases even sharing individual barcodes. We discuss incomplete lineage sorting and hybridization as most likely causes of this phenomenon, as the species concerned are phylogenetically young and hybridization has been previously observed. We also highlight the problem of nuclear mitochondrial pseudogenes (numts), a known problem in the barcoding of orthopteran species, and the possibility of Wolbachia infections. Finally, we discuss the possible taxonomic implications of our barcoding results and point out future research directions.

North African hybrid sparrows (Passer domesticus, P. hispaniolensis) back from oblivion - ecological segregation and asymmetric mitochondrial introgression between parental species

A stabilized hybrid form of the house sparrow (Passer domesticus) and the Spanish sparrow (P. hispaniolensis) is known as Passer italiae from the Italian Peninsula and a few Mediterranean islands. The growing attention for the Italian hybrid sparrow and increasing knowledge on its biology and genetic constitution greatly contrast the complete lack of knowledge of the long‐known phenotypical hybrid sparrow populations from North Africa. Our study provides new data on the breeding biology and variation of mitochondrial DNA in three Algerian populations of house sparrows, Spanish sparrows, and phenotypical hybrids. In two field seasons, the two species occupied different breeding habitats: Spanish sparrows were only found in rural areas outside the cities and bred in open‐cup nests built in large jujube bushes. In contrast, house sparrows bred only in the town centers and occupied nesting holes in walls of buildings. Phenotypical hybrids were always associated with house sparrow populations. House sparrows and phenotypical hybrids started breeding mid of March, and most pairs had three successive clutches, whereas Spanish sparrows started breeding almost one month later and had only two successive clutches. Mitochondrial introgression is strongly asymmetric because about 75% of the rural Spanish sparrow population carried house sparrow haplotypes. In contrast, populations of the Italian hybrid form, P. italiae, were genetically least diverse among all study populations and showed a near‐fixation of house sparrow haplotypes that elsewhere were extremely rare or that were even unique for the Italian Peninsula. Such differences between mitochondrial gene pools of Italian and North African hybrid sparrow populations provide first evidence that different demographic histories have shaped the extant genetic diversity observed on both continents.

Divergence history of the Carpathian and smooth newts modelled in space and time

Information about demographic history is essential for the understanding of the processes of divergence and speciation. Patterns of genetic variation within and between closely related species provide insights into the history of their interactions. Here, we investigated historical demography and genetic exchange between the Carpathian (Lissotriton montandoni, Lm) and smooth (L. vulgaris, Lv) newts. We combine an extensive geographical sampling and multilocus nuclear sequence data with the approximate Bayesian computation framework to test alternative scenarios of divergence and reconstruct the temporal and spatial pattern of gene flow between species. A model of recent (last glacial period) interspecific gene flow was favoured over alternative models. Thus, despite the relatively old divergence (4-6 mya) and presumably long periods of isolation, the species have retained the ability to exchange genes. Nevertheless, the low migration rates (ca. 10(-6) per gene copy per generation) are consistent with strong reproductive isolation between the species. Models allowing demographic changes were favoured, suggesting that the effective population sizes of both species at least doubled as divergence reaching the current ca. 0.2 million in Lm and 1 million in Lv. We found asymmetry in rates of interspecific gene flow between Lm and one evolutionary lineage of Lv. We suggest that intraspecific polymorphism for hybrid incompatibilities segregating within Lv could explain this pattern and propose further tests to distinguish between alternative explanations. Our study highlights the importance of incorporating intraspecific genetic structure into the models investigating the history of divergence.

Genomic heterogeneity of historical gene flow between two species of newts inferred from transcriptome data

The role of gene flow in species formation is a major unresolved issue in speciation biology. Progress in this area requires information on the long-term patterns of gene flow between diverging species. Here, we used thousands of single-nucleotide polymorphisms derived from transcriptome resequencing and a method modeling the joint frequency spectrum of these polymorphisms to reconstruct patterns of historical gene flow between two Lissotriton newts: L. vulgaris (Lv) and L. montandoni (Lm). We tested several models of divergence including complete isolation and various scenarios of historical gene flow. The model of secondary contact received the highest support. According to this model, the species split from their common ancestor ca. 5.5 million years (MY) ago, evolved in isolation for ca. 2 MY, and have been exchanging genes for the last 3.5 MY Demographic changes have been inferred in both species, with the current effective population size of ca. 0.7 million in Lv and 0.2 million in Lm. The postdivergence gene flow resulted in two-directional introgression which affected the genomes of both species, but was more pronounced from Lv to Lm. Interestingly, we found evidence for genomic heterogeneity of interspecific gene flow. This study demonstrates the complexity of long-term gene flow between distinct but incompletely reproductively isolated taxa which divergence was initiated millions of years ago.

Morphological variation in salamanders and their potential response to climate change

Despite the recognition that some species might quickly adapt to new conditions under climate change, demonstrating and predicting such a fundamental response is challenging. Morphological variations in response to climate may be caused by evolutionary changes or phenotypic plasticity, or both, but teasing apart these processes is difficult. Here, we built on the number of thoracic vertebrae (NTV) in ectothermic vertebrates, a known genetically based feature, to establish a link with body size and evaluate how climate change might affect the future morphological response of this group of species. First, we show that in old-world salamanders, NTV variation is strongly related to changes in body size. Secondly, using 22 salamander species as a case study, we found support for relationships between the spatial variation in selected bioclimatic variables and NTV for most of species. For 44% of species, precipitation and aridity were the predominant drivers of geographical variation of the NTV. Temperature features were dominant for 31% of species, while for 19% temperature and precipitation played a comparable role. This two-step analysis demonstrates that ectothermic vertebrates may evolve in response to climate change by modifying the number of thoracic vertebrae. These findings allow to develop scenarios for potential morphological evolution under future climate change and to identify areas and species in which the most marked evolutionary responses are expected. Resistance to climate change estimated from species distribution models was positively related to present-day species morphological response, suggesting that the ability of morphological evolution may play a role for species' persistence under climate change. The possibility that present-day capacity for local adaptation might help the resistance response to climate change can be integrated into analyses of the impact of global changes and should also be considered when planning management actions favouring species persistence.

Contrasting evolutionary histories of the legless lizards slow worms (Anguis) shaped by the topography of the Balkan Peninsula

Background

Genetic architecture of a species is a result of historical changes in population size and extent of distribution related to climatic and environmental factors and contemporary processes of dispersal and gene flow. Population-size and range contractions, expansions and shifts have a substantial effect on genetic diversity and intraspecific divergence, which is further shaped by gene-flow limiting barriers. The Balkans, as one of the most important sources of European biodiversity, is a region where many temperate species persisted during the Pleistocene glaciations and where high topographic heterogeneity offers suitable conditions for local adaptations of populations. In this study, we investigated the phylogeographical patterns and demographic histories of four species of semifossorial slow-worm lizards (genus Anguis) present in the Balkan Peninsula, and tested the relationship between genetic diversity and topographic heterogeneity of the inhabited ranges.

Results

We inferred phylogenetic relationships, compared genetic structure and historical demography of slow worms using nucleotide sequence variation of mitochondrial DNA. Four Anguis species with mostly parapatric distributions occur in the Balkan Peninsula. They show different levels of genetic diversity. A signature of population growth was detected in all four species but with various courses in particular populations. We found a strong correlation between genetic diversity of slow-worm populations and topographic ruggedness of the ranges (mountain systems) they inhabit. Areas with more rugged terrain harbour higher genetic diversity.

Conclusions

Phylogeographical pattern of the genus Anguis in the Balkans is concordant with the refugia-within-refugia model previously proposed for both several other taxa in the region and other main European Peninsulas. While slow-worm populations from the southern refugia mostly have restricted distributions and have not dispersed much from their refugial areas, populations from the extra-Mediterranean refugia in northern parts of the Balkans have colonized vast areas of eastern, central, and western Europe. Besides climatic historical events, the heterogeneous topography of the Balkans has also played an important role in shaping genetic diversity of slow worms.

Electronic supplementary material

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

Deeply divergent sympatric mitochondrial lineages of the earthworm Lumbricus rubellus are not reproductively isolated

Background

The accurate delimitation of species is essential to numerous areas of biological research. An unbiased assessment of the diversity, including the cryptic diversity, is of particular importance for the below ground fauna, a major component of global biodiversity. On the British Isles, the epigeic earthworm Lumbricus rubellus, which is a sentinel species in soil ecotoxicology, consists of two cryptic taxa that are differentiated in both the nuclear and the mitochondrial (mtDNA) genomes. Recently, several deeply divergent mtDNA lineages were detected in mainland Europe, but whether these earthworms also constitute cryptic species remains unclear. This information is important from an evolutionary perspective, but it is also essential for the interpretation and the design of ecotoxicological projects. In this study, we used genome-wide RADseq data to assess the reproductive isolation of the divergent mitochondrial lineages of L. rubellus that occur in sympatry in multiple localities in Central Europe.

Results

We identified five divergent (up to 16 % net p-distance) mitochondrial lineages of L. rubellus in sympatry. Because the clustering of the RADseq data was according to the population of origin and not the mtDNA lineage, reproductive isolation among the mtDNA lineages was not likely. Although each population contained multiple mtDNA lineages, subdivisions within the populations were not observed for the nuclear genome. The lack of fixed differences and sharing of the overwhelming majority of nuclear polymorphisms between localities, indicated that the populations did not constitute allopatric species. The nucleotide diversity within the populations was high, 0.7–0.8 %.

Conclusions

The deeply divergent mtDNA sympatric lineages of L. rubellus in Central Europe were not reproductively isolated groups. The earthworm L. rubellus, which is represented by several mtDNA lineages in continental Europe, apparently is a single highly polymorphic species rather than a complex of several cryptic species. This study demonstrated the critical importance of the use of multilocus nuclear data for the unbiased assessment of cryptic diversity and for the delimitation of species in soil invertebrates.

Electronic supplementary material

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

First results of the German Barcode of Life (GBOL) - Myriapoda project: Cryptic lineages in German Stenotaenialinearis (Koch, 1835) (Chilopoda, Geophilomorpha)

As part of the German Barcode of Life (GBOL) Myriapoda program, which aims to sequence the COI barcoding fragment for 2000 specimens of Germany's 200 myriapod species in the near future, 44 sequences of the centipede order Geophilomorpha are analyzed. The analyses are limited to the genera Geophilus Leach, 1814 and Stenotaenia Koch, 1847 and include a total of six species. A special focus is Stenotaenia, of which 19 specimens from southern, western and eastern Germany could be successfully sequenced. The Stenotaenia data shows the presence of three to four vastly different (13.7-16.7% p-distance) lineages of the genus in Germany. At least two of the three lineages show a wide distribution across Germany, only the lineage including topotypes of Stenotaenialinearis shows a more restricted distribution in southern Germany. In a maximum likelihood phylogenetic analysis the Italian species Stenotaenia 'sorrentina' (Attems, 1903) groups with the different German Stenotaenialinearis clades. The strongly different Stenotaenialinearis lineages within Germany, independent of geography, are a strong hint for the presence of additional, cryptic Stenotaenia species in Germany.

The distribution and taxonomy of Lissotriton newts in Turkey (Amphibia, Salamandridae)

Two and perhaps three taxa of Lissotriton newt occur in Turkey. Their species status is controversial. The distribution of these taxa and the taxonomic status of each are reviewed and discussed. A database of 128 Turkish Lissotriton localities was compiled and species distribution models were constructed. We reiterate that the presence of Lissotriton (vulgaris) lantzi in Turkey is disputed and needs confirmation. The range of Lissotriton (vulgaris) kosswigi is restricted to north-western Anatolia - given the small global range of this Turkey endemic, a closer look at its conservation status is warranted. The distribution of Lissotritonvulgarisschmidtleri covers western Asiatic and European Turkey. The findings support an allopatric distribution of the Turkish Lissotriton species. We reflect on the biological significance of previously reported morphological intermediates between Lissotriton (vulgaris) kosswigi and Lissotritonvulgarisschmidtleri in the light of the recent proposal to recognize kosswigi at the species level. The available data are in line with species status for Lissotriton (vulgaris) lantzi and Lissotriton (vulgaris) kosswigi. Although Lissotritonvulgarisschmidtleri is a genetically diverged taxon as well, the extent of gene flow with parapatric European Lissotriton taxa is as yet unknown.

Constraint and adaptation in newt toll-like receptor genes

Acute die-offs of amphibian populations worldwide have been linked to the emergence of viral and fungal diseases. Inter and intraspecific immunogenetic differences may influence the outcome of infection. Toll-like receptors (TLRs) are an essential component of innate immunity and also prime acquired defenses. We report the first comprehensive assessment of TLR gene variation for urodele amphibians. The Lissotriton newt TLR repertoire includes representatives of 13 families and is compositionally most similar to that of the anuran Xenopus. Both ancient and recent gene duplications have occurred in urodeles, bringing the total number of TLR genes to at least 21. Purifying selection has predominated the evolution of newt TLRs in both long (∼70 Ma) and medium (∼18 Ma) timescales. However, we find evidence for both purifying and positive selection acting on TLRs in two recently diverged (2-5 Ma) allopatric evolutionary lineages (Lissotriton montandoni and L. vulgaris graecus). Overall, both forms of selection have been stronger in L. v. graecus, while constraint on most TLR genes in L. montandoni appears relaxed. The differences in selection regimes are unlikely to be biased by demographic effects because these were controlled by means of a historical demographic model derived from an independent data set of 62 loci. We infer that TLR genes undergo distinct trajectories of adaptive evolution in closely related amphibian lineages, highlight the potential of TLRs to capture the signatures of different assemblages of pathogenic microorganisms, and suggest differences between lineages in the relative roles of innate and acquired immunity.

Testing the efficiency of nested barriers to dispersal in the Mediterranean high mountain plant Edraianthus graminifolius (Campanulaceae)

Due to strong spatial heterogeneity and limited Pleistocene glaciation, the Balkan Peninsula is a major European biodiversity hot spot. Surprisingly little, however, is known about patterns and processes of intraspecific diversification of its biota in general and of high-altitude species in particular. A well-suited system to test hypotheses with respect to various isolating factors acting at different geographic scales and to explore full-range phylogeographic patterns on the Balkan Peninsula is Edraianthus graminifolius (Campanulaceae), distributed in the western Balkan mountain systems, the southwestern Carpathians and the Apennine Peninsula. To this end, we used a dense population sampling and employed amplified fragment length polymorphism (AFLP) markers and plastid DNA sequences supplemented by ecological niche modelling. The strongest splits were inferred to separate southern and northern Balkan populations from the central ones, from where range extension occurred to the Carpathians and, in more recent times, once or twice to the Apennine Peninsula. The three genetic groups in the western Balkan Peninsula were remarkably congruent among molecular markers, suggesting that the barriers to gene flow acted over long time periods facilitating allopatric differentiation. Each main group of Balkan populations contained genetically and geographically distinct subgroups, which likely are the result of local refugia during warmer periods. Evidently, the topographically highly complex and during the Last Glacial Maximum only locally glaciated Balkan Peninsula is a hot spot of species richness and endemism as well as a sanctuary of intraspecific genetic diversity, even if the underlying causes remain insufficiently understood.

Single nucleotide polymorphisms reveal genetic structuring of the carpathian newt and provide evidence of interspecific gene flow in the nuclear genome

Genetic variation within species is commonly structured in a hierarchical manner which may result from superimposition of processes acting at different spatial and temporal scales. In organisms of limited dispersal ability, signatures of past subdivision are detectable for a long time. Studies of contemporary genetic structure in such taxa inform about the history of isolation, range changes and local admixture resulting from geographically restricted hybridization with related species. Here we use a set of 139 transcriptome-derived, unlinked nuclear single nucleotide polymorphisms (SNP) to assess the genetic structure of the Carpathian newt (Lissotriton montandoni, Lm) and introgression from its congener, the smooth newt (L. vulgaris, Lv). Two substantially differentiated groups of Lm populations likely originated from separate refugia, both located in the Eastern Carpathians. The colonization of the present range in north-western and south-western directions was accompanied by a modest loss of variation; admixture between the two groups has occurred in the middle of the Eastern Carpathians. Local, apparently recent introgression of Lv alleles into several Lm populations was detected, demonstrating increased power for admixture detection in comparison to a previous study based on a limited number of microsatellite markers. The level of introgression was higher in Lm populations classified as admixed than in syntopic populations. We discuss the possible causes and propose further tests to distinguish between alternatives. Several outlier loci were identified in tests of interspecific differentiation, suggesting genomic heterogeneity of gene flow between species.

When a clonal genome finds its way back to a sexual species: evidence from ongoing but rare introgression in the hybridogenetic water frog complex

Besides several exceptions, asexual metazoans are usually viewed as ephemeral sinks for genomes, which become 'frozen' in clonal lineages after their emergence from ancestral sexual species. Here, we investigated whether and at what rate the asexuals are able to introgress their genomes back into the parental sexual population, thus more or less importantly affecting the gene pools of sexual species. We focused on hybridogenetic hybrids of western Palaearctic water frogs (Pelophylax esculentus), which originate through hybridization between P. ridibundus and P. lessonae, but transmit only clonal ridibundus genome into their gametes. Although usually mating with P. lessonae, P. esculentus may upon mating with P. ridibundus or another hybrid produce sexually reproducing P. ridibundus offspring with the introgressed ex-clonal genome. We compared the rate of nuclear amplified fragment length polymorphism (AFLP) and mitochondrial introgression in two types of populations, that is, those where P. ridibundus occurs in isolation and those where it lives with the hybridogens. Although significant differentiation (Φpt) between sexual and clonal ridibundus genomes suggested limited gene flow between sexuals and hybridogens, a non-negligible (~5%) proportion of P. ridibundus bore introgressed mtDNA and AFLP markers. Whereas transfer of mtDNA was exclusively unidirectional, introgression of nuclear markers was bidirectional. The proportion of introgressed P. ridibundus was highest in syntopic populations with P. esculentus, proving an ongoing and site-specific interspecific genetic transfer mediated by hybridogenetic hybrids. It turns out that asexual hybrids are not just a sink for genes of sexual species, but may significantly influence the genetic architecture of their sexual counterparts.