Cryptic diversity and phylogeographic patterns of Mediodactylus species in the Eastern Mediterranean region

Cryptic diversity poses a great obstacle in our attempts to assess the current biodiversity crisis and may hamper conservation efforts. The gekkonid genus Mediodactylus, a well-known case of hidden species and genetic diversity, has been taxonomically reclassified several times during the last decade. Focusing on the Mediterranean populations, a recent study within the M. kotschyi species complex using classic mtDNA/nuDNA markers suggested the existence of five distinct species, some being endemic and some possibly threatened, yet their relationships have not been fully resolved. Here, we generated genome-wide SNPs (using ddRADseq) and applied molecular species delimitation approaches and population genomic analyses to further disentangle these relationships. Τhe most extensive nuclear dataset, so far, encompassing 2,360 loci and ∼ 699,000 bp from across the genome of Mediodactylus gecko, enabled us to resolve previously obscure phylogenetic relationships among the five, recently elevated, Mediodactylus species and to support the hypothesis that the taxon includes several new, undescribed species. Population genomic analyses within each of the proposed species showed strong genetic structure and high levels of genetic differentiation among populations.

Whole-genome SNP genotyping unveils ancestral and recent introgression in wild and domestic goats

After the domestication of goats around 10,000 years before the present (BP), humans transported goats far beyond the range of their wild ancestor, the bezoar goat. This brought domestic goats into contact with many wild goat species such as ibex and markhor, enabling introgression between domestic and wild goats. To investigate this, while shedding light on the taxonomic status of wild and domestic goats, we analysed genome-wide SNP data of 613 specimens from 14 taxonomic units, including Capra hircus, C. pyrenaica, C. ibex (from Switzerland, Austria, Germany and Slovenia), C. aegagrus aegagrus, C. a. cretica, C. h. dorcas, C. caucasica caucasica, C. c. severtzovi, C. c. cylindricornis, C. falconeri, C. sibirica sibirica, C. s. alaiana and C. nubiana, as well as Oreamnos americanus (mountain goat) as an outgroup. To trace gene flow between domestic and wild goats, we integrated genotype data of local goat breeds from the Alps as well as from countries such as Spain, Greece, Türkiye, Egypt, Sudan, Iran, Russia (Caucasus and Altai) and Pakistan. Our phylogenetic analyses displayed a clear separation between bezoar-type and ibex-type clades with wild goats from the Greek islands of Crete and Youra clustered within domestic goats, confirming their feral origin. Our analyses also revealed gene flow between the lineages of Caucasian tur and domestic goats that most likely occurred before or during early domestication. Within the clade of domestic goats, analyses inferred gene flow between African and Iberian goats. The detected events of introgression were consistent with previous reports and offered interesting insights into the historical relationships among domestic and wild goats.

Genome sequencing of 2000 canids by the Dog10K consortium advances the understanding of demography, genome function and architecture

BACKGROUND

The international Dog10K project aims to sequence and analyze several thousand canine genomes. Incorporating 20 × data from 1987 individuals, including 1611 dogs (321 breeds), 309 village dogs, 63 wolves, and four coyotes, we identify genomic variation across the canid family, setting the stage for detailed studies of domestication, behavior, morphology, disease susceptibility, and genome architecture and function.

RESULTS

We report the analysis of > 48 M single-nucleotide, indel, and structural variants spanning the autosomes, X chromosome, and mitochondria. We discover more than 75% of variation for 239 sampled breeds. Allele sharing analysis indicates that 94.9% of breeds form monophyletic clusters and 25 major clades. German Shepherd Dogs and related breeds show the highest allele sharing with independent breeds from multiple clades. On average, each breed dog differs from the UU_Cfam_GSD_1.0 reference at 26,960 deletions and 14,034 insertions greater than 50 bp, with wolves having 14% more variants. Discovered variants include retrogene insertions from 926 parent genes. To aid functional prioritization, single-nucleotide variants were annotated with SnpEff and Zoonomia phyloP constraint scores. Constrained positions were negatively correlated with allele frequency. Finally, the utility of the Dog10K data as an imputation reference panel is assessed, generating high-confidence calls across varied genotyping platform densities including for breeds not included in the Dog10K collection.

CONCLUSIONS

We have developed a dense dataset of 1987 sequenced canids that reveals patterns of allele sharing, identifies likely functional variants, informs breed structure, and enables accurate imputation. Dog10K data are publicly available.

Infectious disease threats to amphibians in Greece: new localities positive for Batrachochytrium dendrobatidis

In the early 2000s, numerous cases of European amphibian population declines and mass die-offs started to emerge. Investigating those events led to the discovery that wild European amphibians were confronted with grave disease threats caused by introduced pathogens, namely the amphibian and the salamander chytrid fungi Batrachochytrium dendrobatidis (Bd) and B. salamandrivorans (Bsal) and ranaviruses. In Greece, Bd was previously documented among wild amphibian populations in 2 different locations and 3 different species. However, no disease-related mass declines or mortality events have been reported. In this work, we build upon previous findings with new, subsequently obtained data, resulting in a 225-sample dataset of 14 species from 17 different locations throughout Greece, in order to examine the occurrence status of all 3 pathogens responsible for emerging infectious diseases in European amphibians. No positive samples for Bsal or ranavirus were recorded in any location. We confirmed the presence of Bd in 4 more localities and in 4 more species, including 1 urodelan (Macedonian crested newt Triturus macedonicus) and 1 introduced anuran (American bullfrog Lithobates catesbeianus). All insular localities were negative for Bd, except for Crete, where Bd was identified in 2 different locations. Again, no mass declines or die-offs were recorded in any Bd-positive area or elsewhere. However, given the persistence of Bd across Greece over the past ~20 yr, monitoring efforts should continue, and ideally be further expanded.

Endogenous Viral Elements in Shrew Genomes Provide Insights into Pestivirus Ancient History

As viral genomic imprints in host genomes, endogenous viral elements (EVEs) shed light on the deep evolutionary history of viruses, ancestral host ranges, and ancient viral-host interactions. In addition, they may provide crucial information for calibrating viral evolutionary timescales. In this study, we conducted a comprehensive in silico screening of a large data set of available mammalian genomes for EVEs deriving from members of the viral family Flaviviridae, an important group of viruses including well-known human pathogens, such as Zika, dengue, or hepatitis C viruses. We identified two novel pestivirus-like EVEs in the reference genome of the Indochinese shrew (Crocidura indochinensis). Homologs of these novel EVEs were subsequently detected in vivo by molecular detection and sequencing in 27 shrew species, including 26 species representing a wide distribution within the Crocidurinae subfamily and one in the Soricinae subfamily on different continents. Based on this wide distribution, we estimate that the integration event occurred before the last common ancestor of the subfamily, about 10.8 million years ago, attesting to an ancient origin of pestiviruses and Flaviviridae in general. Moreover, we provide the first description of Flaviviridae-derived EVEs in mammals even though the family encompasses numerous mammal-infecting members. This also suggests that shrews were past and perhaps also current natural reservoirs of pestiviruses. Taken together, our results expand the current known Pestivirus host range and provide novel insight into the ancient evolutionary history of pestiviruses and the Flaviviridae family in general.

Ecological niche differentiation in the Anatolian rock lizards (Genus: Anatololacerta) (Reptilia: Lacertidae) of the Anatolian Peninsula and Aegean Islands

The genus Anatololacerta is distributed in the eastern Mediterranean region including Asia Minor and some east Aegean islands. Recent phylogenetic studies suggested that this genus displayed cryptic diversity and was divided into five species: A. anatolica, A. pelasgiana, A. ibrahimi, A. finikensis and A. danfordi. The ecological niche differentiations of these species have not been studied so far. Our aims for this study were to predict the potential suitable habitats for the species nested in genus Anatololacerta, and to examine the niche overlaps and differentiations via identity and background tests. The occurrence data were obtained from literature and our own field surveys. Occurrence records were rarefied and assessed in a 30 arc-second resolution layer, compatible with several bioclimatic and topographic variables. Species distribution analyses were performed using maximum entropy approach and pairwise niche comparisons were evaluated by identity and background tests. Our results demonstrated that the species delimitation among this genus was not only affected by geographic isolation but also that precipitation and temperature influenced the habitat suitability for these species. Predicted suitability usually well matched the actual species distributions. Moreover, the niche overlap (identity test) analyses verified that allopatric Anatololacerta species show clear ecological differentiations. However, a niche overlap between parapatric species A. pelasgiana and A. finikensis, was confirmed by identity and background tests. It has been suggested that these parapatric species could be more affected by microclimatological parameters than the others. The results of our study are in agreement with the latest phylogenetic study within this genus.

The story of a rock-star: multilocus phylogeny and species delimitation in the starred or roughtail rock agama, Laudakia stellio (Reptilia: Agamidae)

Situated at the junction of three continents, Europe, Asia and Africa, the Eastern Mediterranean is an ideal region to study the effects of palaeogeography, ecology and long human presence on animal evolution. Laudakia stellio (Squamata: Agamidae) is found across this region and offers an excellent opportunity for such studies. The high morphological variation across their range suggests that these lizards might represent a species complex. This is the first study exploring their evolutionary history, using molecular markers and individuals from all described subspecies. We employed the latest phylogenetic and species-delimitation methods to identify all distinct evolutionary lineages, their genetic variation and divergence times. The phenotypical diversity of L. stellio matches its genetic differentiation: almost all subspecies correspond to well-supported retrieved subclades and additional distinct lineages representing intermediate morphs have been retrieved. ‘Laudakia stellio’ represents three distinct evolutionary entities that diverged during the Plio-Pleistocene transition, which we propose as distinct species. One includes Greek and Turkish populations, as well as cryptic Anatolian lineages. The second comprises all other Near East populations and the third is endemic to Cyprus. Our results indicate a role of humans in shaping present distribution patterns, and highlight the importance of the Aegean, Anatolia and the Levant as glacial refugia and diversity hotspots.

Genomic consequences of colonisation, migration and genetic drift in barn owl insular populations of the eastern Mediterranean

The study of insular populations was key in the development of evolutionary theory. The successful colonisation of an island depends on the geographic context, and specific characteristics of the organism and the island, but also on stochastic processes. As a result, apparently identical islands may harbour populations with contrasting histories. Here, we use whole genome sequences of 65 barn owls to investigate the patterns of inbreeding and genetic diversity of insular populations in the eastern Mediterranean Sea. We focus on Crete and Cyprus, islands with similar size, climate and distance to mainland, that provide natural replicates for a comparative analysis of the impacts of microevolutionary processes on isolated populations. We show that barn owl populations from each island have a separate origin, Crete being genetically more similar to other Greek islands and mainland Greece, and Cyprus more similar to the Levant. Further, our data show that their respective demographic histories following colonisation were also distinct. On the one hand, Crete harbours a small population and maintains very low levels of gene flow with neighbouring populations. This has resulted in low genetic diversity, strong genetic drift, increased relatedness in the population and remote inbreeding. Cyprus, on the other hand, appears to maintain enough gene flow with the mainland to avoid such an outcome. Our study provides a comparative population genomic analysis of the effects of neutral processes on a classical island‐mainland model system. It provides empirical evidence for the role of stochastic processes in determining the fate of diverging isolated populations.

Population genomics of wall lizards reflects the dynamic history of the Mediterranean Basin

The Mediterranean Basin has experienced extensive change in geology and climate over the past six million years. Yet, the relative importance of key geological events for the distribution and genetic structure of the Mediterranean fauna remains poorly understood. Here, we use population genomic and phylogenomic analyses to establish the evolutionary history and genetic structure of common wall lizards (Podarcis muralis). This species is particularly informative because, in contrast to other Mediterranean lizards, it is widespread across the Iberian, Italian, and Balkan Peninsulas, and in extra-Mediterranean regions. We found strong support for six major lineages within P. muralis, which were largely discordant with the phylogenetic relationship of mitochondrial DNA. The most recent common ancestor of extant P. muralis was likely distributed in the Italian Peninsula, and experienced an “Out-of-Italy” expansion following the Messinian salinity crisis (∼5 Mya), resulting in the differentiation into the extant lineages on the Iberian, Italian, and Balkan Peninsulas. Introgression analysis revealed that both inter- and intraspecific gene flows have been pervasive throughout the evolutionary history of P. muralis. For example, the Southern Italy lineage has a hybrid origin, formed through admixture between the Central Italy lineage and an ancient lineage that was the sister to all other P. muralis. More recent genetic differentiation is associated with the onset of the Quaternary glaciations, which influenced population dynamics and genetic diversity of contemporary lineages. These results demonstrate the pervasive role of Mediterranean geology and climate for the evolutionary history and population genetic structure of extant species.

The distribution and biogeography of slow worms (Anguis, Squamata) across the Western Palearctic, with an emphasis on secondary contact zones

The slow-worm lizards (Anguis) comprise five species occurring throughout most of the Western Palearctic. Although these species are relatively uniform morphologically – with the exception of A. cephallonica, which exhibits a quite unique morphology – they are genetically deeply divergent. Here, we provide detailed distribution maps for each species and discuss their biogeography and conservation based on updated genetic data and a robust distribution database. We pay particular attention to the so called ‘grey zone’, which typically represents secondary contact zones and in some cases confirmed or presumed hybrid zones. Four of the five species live in parapatry, while only two species, A. cephallonica and A. graeca from the southern Balkans occur in partial sympatry. Further research should focus on the eco-evolutionary interactions between species in contact, including their hybridization rates, to reveal deeper details of the slow-worm evolutionary and natural history.

The Island of Extremes: Giants and Dwarfs on a Small Remote Island

Body size evolution on islands is widely studied and hotly debated. Gigantism and dwarfism are thought to evolve under strong natural selection, especially on small remote islands. We report a curious co-occurrence of both dwarf and giant lizards on the same small, remote island (Plakida): the largest Podarcis erhardii (Lacertidae) and smallest Mediodactylus kotschyi sensu lato; Gekkonidae — the two commonest insular reptiles in the Aegean Sea. The geckos of Plakida have a peculiar tail-waving behavior, documented here for the first time in this genus. We suspect that P. erhardii evolved large size to consume geckos and the geckos evolved a unique tail-waving behavior as a defensive mechanism.

Endoparasites of European Wildcats (Felis silvestris) in Greece

The European wildcat (Felis silvestris) is the only wild felid living in Greece. Wildcat populations are declining due to anthropogenic and phenological unfavourable conditions, and parasites may have an additional negative impact. In the present study, the occurrence of endoparasites in wildcats in Greece and the potential threats posed to wildcats, domestic animals, and humans in the study areas has been investigated. In a six-year period, 23 road-killed wildcats and 62 wildcat faecal samples were collected from different areas of the country. Necropsy for the detection of endoparasites and standard parasitological examinations of faecal samples were performed. Parasites were morphologically identified and, in selected cases, molecularly analysed. All necropsied wildcats (100%) were infected by three to 10 different parasite taxa, with the most prevalent being Taenia taeniaeformis (73.9%), Toxocara cati (60.9%), Angiostrongylus chabaudi (56.5%), Ancylostoma tubaeforme (39.1%), Cylicospirura spp. (34.8%), Troglostrongylus brevior (34.8%), and Capillaria aerophila (33.8%). Of the 62 faecal samples examined, 53 (85.5%) were positive for one or more parasite elements (larvae, eggs, or oocysts). The most frequent were T. cati (45.2%), A. chabaudi (29%), C. aerophila (24.2%), and Ancylostomatidae (17.7%). This is the first survey on endoparasites affecting wildcats in Greece. Some of the parasites here found are frequent in domestic and wild felids, while others, i.e., Oslerus rostratus and Cylicospirura petrowi, were described for the first time in the European wildcat. Most of them have a significant pathogenic potential, causing severe to hazardous diseases to infected felids and some, under specific circumstances, can also threaten human health.

Long-term stability of sex chromosome gene content allows accurate qPCR-based molecular sexing across birds

Embryos, juveniles, and even adults of many bird species lack pronounced external sexually dimorphic characteristics. Accurate identification of sex is crucial for research (e.g., developmental, population, and evolutionary studies), management of wildlife species, and captive breeding programmes for both conservation and poultry. An accurate molecular sexing method applicable across the entire bird radiation is theoretically possible thanks to the long-term stability of their ZZ/ZW sex chromosomes, but current methods are not applicable in a wide range of bird lineages. Here, we developed a novel molecular sexing method based on the comparison of gene copy number variation by quantitative real-time PCR (qPCR) in conserved Z-specific genes (CHRNA6, DDX4, LPAR1, TMEM161B, VPS13A), i.e. genes linked to Z but absent from W chromosomes. We tested the method across three paleognath and 70 neognath species covering the avian phylogeny. In addition, we designed primers for four Z-specific genes (DOCK8, FUT10, PIGG and PSD3) for qPCR-based molecular sexing in three paleognath species. We have demonstrated that the genes DOCK8, FUT10, PIGG and PSD3 can identify sex in paleognath birds and the genes CHRNA6, DDX4, TMEM161B, and VPS13A can reveal sex in neognath birds. The gene LPAR1 can be used to accurately identify sex in both paleognath and neognath species. Along with outlining a novel method of practical importance for molecular sexing in birds, our study also documents in detail the conservation of sex chromosomes across the avian phylogeny.

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.

The wall lizards of the Balkan peninsula: Tackling questions at the interface of phylogenomics and population genomics

Wall lizards of the genus Podarcis (Sauria, Lacertidae) are the predominant reptile group in southern Europe, including 24 recognized species. Mitochondrial DNA data have shown that, with the exception of P. muralis, the Podarcis species distributed in the Balkan peninsula form a species group that is further sub-divided into two subgroups: the one of "P. tauricus" consisting of P. tauricus, P. milensis, P. gaigeae, and P. melisellensis, and the other of "P. erhardii" comprising P. erhardii, P. levendis, P. cretensis, and P. peloponnesiacus. In an attempt to explore the Balkan Podarcis phylogenomic relationships, assess the levels of genetic structure and to re-evaluate the number of extant species, we employed phylogenomic and admixture approaches on ddRADseq (double digested Restriction site Associated DNA sequencing) genomic data. With this efficient Next Generation Sequencing approach, we were able to obtain a large number of genomic loci randomly distributed throughout the genome and use them to resolve the previously obscure phylogenetic relationships among the different Podarcis species distributed in the Balkans. The obtained phylogenomic relationships support the monophyly of both aforementioned subgroups and revealed several divergent lineages within each subgroup, stressing the need for taxonomic re-evaluation of Podarcis' species in Balkans. The phylogenomic trees and the species delimitation analyses confirmed all recently recognized species (P. levendis, P. cretensis, and P. ionicus) and showed the presence of at least two more species, one in P. erhardii and the other in P. peloponnesiacus.

More evidence of cryptic diversity in Anatololacerta species complex Arnold, Arribas and Carranza, 2007 (Squamata: Lacertidae) and re-evaluation of its current taxonomy

Genetic diversity is not always congruent with phenotypic heterogeneity, resulting in cryptic species complexes which cause a great struggle for scientists trying to define ‘species’ and describe relationships among taxa. Anatololacerta is a lizard genus distributed in southern and western Anatolia and some neighboring Aegean islands. Three morphospecies were recognized in Anatololacerta but a recent molecular study revealed the presence of cryptic diversity within the genus which led to the raise of a subspecies to species level. Currently the genus includes the species A. anatolica, A. danfordi, A. budaki and A. pelasgiana. Using a comprehensive sampling concerning both the number of specimens (218 specimens) and the genetic markers (3 nuclear and 3 mitochondrial), we performed phylogenetic analyses including tree reconstruction, species delimitation and divergence times estimation. The results revealed the occurrence of one more cryptic lineage which should be regarded as a separate species for which the name A. ibrahimi stat. nov. has priority. The existence of five well differentiated species with parapatric distributions in Anatololacerta is strongly supported. There is also evidence of recent and rapid radiation of the genus which probably causes phylogenetic relationships between these species to remain largely unresolved. At last, we proceeded to some nomenclatorial changes: The current name A. budaki was synonymized with A. pelasgiana because specimens of the type-locality of A. budaki are assigned genetically to A. pelasgiana. The genetic lineage including specimens currently assigned to A. budaki was named A. finikensis stat. nov., raising the subspecies A. b. finikensis to species level.

Landscape Connectivity Limits the Predicted Impact of Fungal Pathogen Invasion

Infectious diseases are major drivers of biodiversity loss. The risk of fungal diseases to the survival of threatened animals in nature is determined by a complex interplay between host, pathogen and environment. We here predict the risk of invasion of populations of threatened Mediterranean salamanders of the genus Lyciasalamandra by the pathogenic chytrid fungus Batrachochytrium salamandrivorans by combining field sampling and lab trials. In 494 samples across all seven species of Lyciasalamandra, B. salamandrivorans was found to be absent. Single exposure to a low (1000) number of fungal zoospores resulted in fast buildup of lethal infections in three L. helverseni. Thermal preference of the salamanders was well within the thermal envelope of the pathogen and body temperatures never exceeded the fungus' thermal critical maximum, limiting the salamanders' defense opportunities. The relatively low thermal host preference largely invalidates macroclimatic based habitat suitability predictions and, combined with current pathogen absence and high host densities, suggests a high probability of local salamander population declines upon invasion by B. salamandrivorans. However, the unfavorable landscape that shaped intraspecific host genetic diversity, lack of known alternative hosts and rapid host mortality after infection present barriers to further, natural pathogen dispersal between populations and thus species extinction. The risk of anthropogenic spread stresses the importance of biosecurity in amphibian habitats.

Species list of the European herpetofauna – 2020 update by the Taxonomic Committee of the Societas Europaea Herpetologica

The last species list of the European herpetofauna was published by Speybroeck, Beukema and Crochet (2010). In the meantime, ongoing research led to numerous taxonomic changes, including the discovery of new species-level lineages as well as reclassifications at genus level, requiring significant changes to this list. As of 2019, a new Taxonomic Committee was established as an official entity within the European Herpetological Society, Societas Europaea Herpetologica (SEH). Twelve members from nine European countries reviewed, discussed and voted on recent taxonomic research on a case-by-case basis. Accepted changes led to critical compilation of a new species list, which is hereby presented and discussed. According to our list, 301 species (95 amphibians, 15 chelonians, including six species of sea turtles, and 191 squamates) occur within our expanded geographical definition of Europe. The list includes 14 non-native species (three amphibians, one chelonian, and ten squamates).

Morphometrics and genetics highlight the complex history of Eastern Mediterranean spiny mice

Spiny mice of the Acomys cahirinus group display a complex geographical structure in the Eastern Mediterranean area, as shown by previous genetic and chromosomal studies. To better elucidate the evolutionary relationships between insular populations from Crete and Cyprus and continental populations from North Africa and Cilicia in Turkey, genetic and morphometric variations were investigated, based on mitochondrial D-loop sequences, and the size and shape of the first upper molar. The Cypriot and the Cilician populations show idiosyncratic divergence in molar size and shape, while Cretan populations present a geographical structure with at least three differentiated subpopulations, as shown by congruent distributions of haplogroups, Robertsonian fusions and morphometric variation. A complex history of multiple introductions is probably responsible for this structure, and insular isolation coupled with habitat shift should have further promoted a pronounced and rapid morphological evolution in molar size and shape on Crete and Cyprus.

Impacts of the synthetic androgen Trenbolone on gonad differentiation and development - comparisons between three deeply diverged anuran families

Using a recently developed approach for testing endocrine disruptive chemicals (EDCs) in amphibians, comprising synchronized tadpole exposure plus genetic and histological sexing of metamorphs in a flow-through-system, we tested the effects of 17β-Trenbolone (Tb), a widely used growth promoter in cattle farming, in three deeply diverged anuran families: the amphibian model species Xenopus laevis (Pipidae) and the non-models Bufo(tes) viridis (Bufonidae) and Hyla arborea (Hylidae). Trenbolone was applied in three environmentally and/or physiologically relevant concentrations (0.027 µg/L (10-10 M), 0.27 µg/L (10-9 M), 2.7 µg/L (10-8 M)). In none of the species, Tb caused sex reversals or masculinization of gonads but had negative species-specific impacts on gonad morphology and differentiation after the completion of metamorphosis, independently of genetic sex. In H. arborea and B. viridis, mounting Tb-concentration correlated positively with anatomical abnormalities at 27 µg/L (10-9 M) and 2.7 µg/L (10-8 M), occurring in X. laevis only at the highest Tb concentration. Despite anatomical aberrations, histologically all gonadal tissues differentiated seemingly normally when examined at the histological level but at various rates. Tb-concentration caused various species-specific mortalities (low in Xenopus, uncertain in Bufo). Our data suggest that deep phylogenetic divergence modifies EDC-vulnerability, as previously demonstrated for Bisphenol A (BPA) and Ethinylestradiol (EE2).

Phylogeography of a cryptic speciation continuum in Eurasian spadefoot toads (Pelobates)

Cryptic phylogeographic diversifications provide unique models to examine the role of phylogenetic divergence on the evolution of reproductive isolation, without extrinsic factors such as ecological and behavioural differentiation. Yet, to date very few comparative studies have been attempted within such radiations. Here, we characterize a new speciation continuum in a group of widespread Eurasian amphibians, the Pelobates spadefoot toads, by conducting multilocus (restriction site associated DNA sequencing and mitochondrial DNA) phylogenetic, phylogeographic and hybrid zone analyses. Within the P. syriacus complex, we discovered species-level cryptic divergences (>5 million years ago [My]) between populations distributed in the Near-East (hereafter P. syriacus sensu stricto [s.s.]) and southeastern Europe (hereafter P. balcanicus), each featuring deep intraspecific lineages. Altogether, we could scale hybridizability to divergence time along six different stages, spanning from sympatry without gene flow (P. fuscus and P. balcanicus, >10 My), parapatry with highly restricted hybridization (P. balcanicus and P. syriacus s.s., >5 My), narrow hybrid zones (~15 km) consistent with partial reproductive isolation (P. fuscus and P. vespertinus, ~3 My), to extensive admixture between Pleistocene and refugial lineages (≤2 My). This full spectrum empirically supports a gradual build up of reproductive barriers through time, reversible up until a threshold that we estimate at ~3 My. Hence, cryptic phylogeographic lineages may fade away or become reproductively isolated species simply depending on the time they persist in allopatry, and without definite ecomorphological divergence.

Genomic and mitochondrial evidence of ancient isolations and extreme introgression in the four-lined snake

Comparing mitochondrial and genomic phylogenies is an essential tool for investigating speciation processes, because each genome carries different inheritance properties and evolutionary characteristics. Furthermore, mitonuclear discordance may arise from ecological adaptation, historic isolation, population size changes, and sex-biased dispersal. Closely related taxa are expected to experience gene flow; however, this may not be true for insular populations or populations isolated in refugia. The four-lined snake Elaphe quatuorlineata has a fragmented distribution, separating populations of the Italian and Balkan Peninsulas, whereas several insular Aegean populations of significantly smaller body size (Cyclades island group and Skyros Island, Greece) are currently considered distinct subspecies. We constructed the species-tree phylogeny of this species utilizing genome-wide single nucleotide polymorphisms and a gene-tree based on complete cytochrome b sequences, aiming to detect convergence and discrepancies between biparentally and maternally inherited genomes. Population structuring, phylogenetic patterns and migration events among geographically defined lineages supported our hypothesis of isolation in multiple sub-refugia. Where biogeographical barriers did not restrict subsequent dispersal, extensive genetic exchange occurred between mainland Balkan populations. This process has led to the mitochondrial sweep of an ancestral mitolineage that survived only in peripheral (East Greece) and insular populations (North Cyclades and Skyros). The Central Cyclades represent an ancient lineage for both molecular markers that emerged almost 3.3 Mya. Considering their distinct morphology, insular E. quatuorlineata populations should be the future focus of an extensive sampling, especially since the mitonuclear discordance observed in this species could be related to ecological adaptations, such as the island-dwarfism phenomenon.

Selection, drift, and introgression shape MHC polymorphism in lizards

The major histocompatibility complex (MHC) has long served as a model for the evolution of adaptive genetic diversity in wild populations. Pathogen-mediated selection is thought to be a main driver of MHC diversity, but it remains elusive to what degree selection shapes MHC diversity in complex biogeographical scenarios where other evolutionary processes (e.g. genetic drift and introgression) may also be acting. Here we focus on two closely related green lizard species, Lacerta trilineata and L. viridis, to address the evolutionary forces acting on MHC diversity in populations with different biogeographic structure. We characterized MHC class I exon 2 and exon 3, and neutral diversity (microsatellites), to study the relative importance of selection, drift, and introgression in shaping MHC diversity. As expected, positive selection was a significant force shaping the high diversity of MHC genes in both species. Moreover, introgression significantly increased MHC diversity in mainland populations, with a primary direction of gene flow from L. viridis to L. trilineata. Finally, we found significantly fewer MHC alleles in island populations, but maintained MHC sequence and functional diversity, suggesting that positive selection counteracted the effect of drift. Overall, our data support that different evolutionary processes govern MHC diversity in different biogeographical scenarios: positive selection occurs broadly while introgression acts in sympatry and drift when the population sizes decrease.

Phylogeography of Aegean green toads (Bufo viridis subgroup): continental hybrid swarm vs. insular diversification with discovery of a new island endemic

Background

Debated aspects in speciation research concern the amount of gene flow between incipient species under secondary contact and the modes by which post-zygotic isolation accumulates. Secondary contact zones of allopatric lineages, involving varying levels of divergence, provide natural settings for comparative studies, for which the Aegean (Eastern Mediterranean) geography offers unique scenarios. In Palearctic green toads (Bufo viridis subgroup or Bufotes), Plio-Pleistocene (~ 2.6 Mya) diverged species show a sharp transition without contemporary gene flow, while younger lineages, diverged in the Lower-Pleistocene (~ 1.9 Mya), admix over tens of kilometers. Here, we conducted a fine-scale multilocus phylogeographic analysis of continental and insular green toads from the Aegean, where a third pair of taxa, involving Mid-Pleistocene diverged (~ 1.5 Mya) mitochondrial lineages, earlier tentatively named viridis and variabilis, (co-)occurs.

Results

We discovered a new lineage, endemic to Naxos (Central Cyclades), while coastal islands and Crete feature weak genetic differentiation from the continent. In continental Greece, both lineages, viridis and variabilis, form a hybrid swarm, involving massive mitochondrial and nuclear admixture over hundreds of kilometers, without obvious selection against hybrids.

Conclusions

The genetic signatures of insular Aegean toads appear governed by bathymetry and Quaternary sea level changes, resulting in long-term isolation (Central Cyclades: Naxos) and recent land-bridges (coastal islands). Conversely, Crete has been isolated since the end of the Messinian salinity crisis (5.3 My) and Cretan populations thus likely result from human-mediated colonization, at least since Antiquity, from Peloponnese and Anatolia. Comparisons of green toad hybrid zones support the idea that post-zygotic hybrid incompatibilities accumulate gradually over the genome. In this radiation, only one million years of divergence separate a scenario of complete reproductive isolation, from a secondary contact resulting in near panmixia.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1179-0) contains supplementary material, which is available to authorized users.

Glacial allopatry vs. postglacial parapatry and peripatry: the case of hedgehogs

Although hedgehogs are well-known examples of postglacial recolonisation, the specific processes that shape their population structures have not been examined by detailed sampling and fast-evolving genetic markers in combination with model based clustering methods. This study aims to analyse the impacts of isolation within glacial refugia and of postglacial expansion on the population structure of the Northern White-breasted hedgehog (Erinaceus roumanicus). It also discusses the role of the processes at edges of species distribution in its evolutionary history. The maternally inherited mitochondrial control region and the bi-parentally inherited nuclear microsatellites were used to examine samples within the Central Europe, Balkan Peninsula and adjacent islands. Bayesian coalescent inference and neutrality tests proposed a recent increase in the population size. The most pronounced pattern of population structure involved differentiation of the insular populations in the Mediterranean Sea and the population within the contact zone with E. europaeus in Central Europe. An interspecies hybrid was detected for the first time in Central Europe. A low genetic diversity was observed in Crete, while the highest genetic distances among individuals were found in Romania. The recent population in the post-refugial area related to the Balkan Peninsula shows a complex pattern with pronounced subpopulations located mainly in the Pannonian Basin and at the Adriatic and Pontic coasts. Detailed analyses indicate that parapatry and peripatry may not be the only factors that limit range expansion, but also strong microevolutionary forces that may change the genetic structure of the species. Here we present evidence showing that population differentiation may occur not only during the glacial restriction of the range into the refugia, but also during the interglacial range expansion. Population differentiation at the Balkan Peninsula and adjacent regions could be ascribed to diversification in steppe/forest biomes and complicated geomorphology, including pronounced geographic barriers as Carpathians.

Evolutionary stability of sex chromosomes in snakes

Amniote vertebrates possess various mechanisms of sex determination, but their variability is not equally distributed. The large evolutionary stability of sex chromosomes in viviparous mammals and birds was believed to be connected with their endothermy. However, some ectotherm lineages seem to be comparably conserved in sex determination, but previously there was a lack of molecular evidence to confirm this. Here, we document a stability of sex chromosomes in advanced snakes based on the testing of Z-specificity of genes using quantitative PCR (qPCR) across 37 snake species (our qPCR technique is suitable for molecular sexing in potentially all advanced snakes). We discovered that at least part of sex chromosomes is homologous across all families of caenophidian snakes (Acrochordidae, Xenodermatidae, Pareatidae, Viperidae, Homalopsidae, Colubridae, Elapidae and Lamprophiidae). The emergence of differentiated sex chromosomes can be dated back to about 60 Ma and preceded the extensive diversification of advanced snakes, the group with more than 3000 species. The Z-specific genes of caenophidian snakes are (pseudo)autosomal in the members of the snake families Pythonidae, Xenopeltidae, Boidae, Erycidae and Sanziniidae, as well as in outgroups with differentiated sex chromosomes such as monitor lizards, iguanas and chameleons. Along with iguanas, advanced snakes are therefore another example of ectothermic amniotes with a long-term stability of sex chromosomes comparable with endotherms.

Impaired gonadal and somatic development corroborate vulnerability differences to the synthetic estrogen ethinylestradiol among deeply diverged anuran lineages

Amphibians are undergoing a global decline. One poorly investigated reason could be the pollution of aquatic habitats by endocrine disrupting compounds (EDCs). We tested the susceptibility to the synthetically stabilized estrogen 17α-ethinylestradiol (EE2) in three deeply diverged anuran species, differing in sex determination systems, types of gonadogenesis and larval ecologies. To understand whether data from the amphibian model Xenopus laevis (Pipidae) are analogous and applicable to only distantly related non-model amphibians, tadpoles of X. laevis, Hyla arborea (Hylidae) and Bufo viridis (Bufonidae) were simultaneously exposed to 50, 500 and 5000ng/L EE2 from hatching until completion of metamorphosis, using a flow-through-system under identical experimental conditions. Comparing molecularly established genetic with histologically assessed phenotypic sex in all species, we have recently shown that EE2 provoked numerous genetic-male-to-phenotypic-female sex reversals and mixed sex individuals, confirming overall its expected feminizing effect. In the present study, we focus on the influence of EE2 on gonadal and somatic development. Anatomy and histology revealed several species-specific effects. In both non-model species, H. arborea and B. viridis, high numbers of anatomically impaired gonads were observed. In H. arborea, exposed to 5000ng/L EE2, numerous underdeveloped gonads were detected. Whereas EE2 did not alter snout-to-vent length and body weight of X. laevis metamorphs, H. arborea showed a treatment-dependent decrease, while B. viridis exhibited an increase in body weight and snout-to-vent length. Apart from a concentration-dependent occurrence of yellowish skin color in several H. arborea, no organ-specific effects were detected. Since EE2 ubiquitously occurs in many aquatic ecosystems and affects sexual and somatic development, among EDCs, it may indeed contribute to amphibian decline. The inter-species variation in developmental EE2-effects corroborates species-specific vulnerability differences towards EDCs between deeply diverged amphibian groups.

A reinvestigation of phylogeny and divergence times of the Ablepharus kitaibelii species complex (Sauria, Scincidae) based on mtDNA and nuDNA genes

Morphological and DNA data support that the East Mediterranean snake-eyed skink Ablepharus kitaibelii represents a species complex that includes four species A. kitaibelii, A. budaki, A. chernovi, and A. rueppellii, highlighting the need of its taxonomic reevaluation. Here, we used Bayesian and Maximum Likelihood methods to estimate the phylogenetic relationships of all members of the complex based on two mitochondrial (cyt b, 16S rRNA) and two nuclear markers (MC1R, and NKTR) and using Chalcides, Eumeces, and Eutropis as outgroups. The biogeographic history of the complex was also investigated through the application of several phylogeographic (BEAST) and biogeographic (BBM) analyses. Paleogeographic and paleoclimatic data were used to support the inferred phylogeographic patterns. The A. kitaibelli species complex exhibits high genetic diversity, revealing cases of hidden diversity and cases of non-monophyletic species such as A. kitaibelii and A. budaki. Our results indicate that A. pannonicus branches off first and a group that comprises specimens of A. kitaibelli and A. budaki from Kastelorizo Island group (southeast Greece) and southwest Turkey, respectively is differentiated from the rest A. kitaibelli and A. budaki populations and may represent a new species. The estimated divergence times place the origin of the complex in the Middle Miocene (∼16Mya) and the divergence of most currently recognized species in the Late Miocene. The inferred ancestral distribution suggests that the complex originated in Anatolia, supposing that several vicariance and dispersal events that are related with the formation of the Mid-Aegean Trench, the Anatolian Diagonal and the orogenesis of the mountain chains in southern and eastern Anatolia have led to current distribution pattern of A. kitaibelii species complex in the Balkans and Middle East.

Sex reversal assessments reveal different vulnerability to endocrine disruption between deeply diverged anuran lineages

Multiple anthropogenic stressors cause worldwide amphibian declines. Among several poorly investigated causes is global pollution of aquatic ecosystems with endocrine disrupting compounds (EDCs). These substances interfere with the endocrine system and can affect the sexual development of vertebrates including amphibians. We test the susceptibility to an environmentally relevant contraceptive, the artificial estrogen 17α-ethinylestradiol (EE2), simultaneously in three deeply divergent systematic anuran families, a model-species, Xenopus laevis (Pipidae), and two non-models, Hyla arborea (Hylidae) and Bufo viridis (Bufonidae). Our new approach combines synchronized tadpole exposure to three EE2-concentrations (50, 500, 5,000 ng/L) in a flow-through-system and pioneers genetic and histological sexing of metamorphs in non-model anurans for EDC-studies. This novel methodology reveals striking quantitative differences in genetic-male-to-phenotypic-female sex reversal in non-model vs. model species. Our findings qualify molecular sexing in EDC-analyses as requirement to identify sex reversals and state-of-the-art approaches as mandatory to detect species-specific vulnerabilities to EDCs in amphibians.

Complete mitochondrial genome of the endemic legless lizard Anguis cephallonica Werner, 1894 and its comparison with mitogenome of Anguis fragilis Linnaeus, 1758

Complete mitochondrial genome of the Peloponnese endemic lizard species Anguis cephallonica is presented in this study. The complete sequence is 17 208 bp long and consists of 13 protein-coding genes, 22 tRNA genes, two rRNA genes and one control region. The gene order is same as in the relative species Anguis fragilis. Length of the 22 tRNA genes varies from 64 bp to 73 bp. The Anguis cephallonica mitogenome base composition is: A (30.5%), T (24.2%), C (30.5%), G (14.8%), with an A + T bias (54.7%). Six protein coding genes have incomplete stop codons. This is the first complete mitogenome described in this species as well as in any endemic Peloponnese lizard. Presented complete mitochondrial genome will form a basis for future comparative analysis within the genus Anguis.