On the origin of the recent herpetofauna of Sicily: Comparative phylogeography using homologous mitochondrial and nuclear genes

DNA barcoding reveals cryptic species in the sea slater Ligia italica (Crustacea, Isopoda) from Tunisia

Barcoding studies have provided significant insights into phylogenetic relationships among species belonging to the genus Ligia (Crustacea, Isopoda). Herein the diversity of the Italian sea slater Ligia italica from Tunisia is studied for the first time. Samples were collected from 18 localities in Tunisia, and the analysis included previously published sequences from Italy and Greece available in GenBank. Bayesian and Maximum Likelihood phylogenetic analyses were carried out using a fragment of the mitochondrial COI gene. Putative cryptic species were explored using the 'barcode gap' approach in the software ASAP. A genetic landscape shape analysis was carried out using the program Alleles in Space. The analyses revealed highly divergent and well-supported clades of L. italica dispersed across Tunisia (Clades A1 and A2), Greece (Clade B) and Italy (Clades C1 and C2). High genetic dissimilarity among clades suggested that L. italica constitute a cryptic species complex. Divergence among different L. italica lineages (Clades A, B and C) occurred around 7-4.5 Ma. The detected high genetic distances among clades did not result from atypical mitochondrial DNAs or intracellular infection by Wolbachia bacteria. The complex history of the Mediterranean Sea appears to have played a significant role in shaping the phylogeographic pattern of Ligia italica. Additional morphological and molecular studies are needed to confirm the existence of cryptic species in Ligia italica in Mediterranean.

Colonization routes uncovered in a widely introduced Mediterranean gecko, Tarentola mauritanica

In this study, we aimed to understand the contemporary and ancient colonization routes of the Moorish gecko, Tarentola mauritanica, using simple sequence repeats. By analyzing the genetic diversity of populations in different regions, we found that Morocco is the genetic diversity hotspot for the species, followed by the Iberian Peninsula. However, historical gene flow estimates identified the Iberian Peninsula, not Morocco, as the primary contributor of colonizing individuals, along with continental Italy to a lesser extent. Currently, mainland Italy is the main source of introduced individuals, likely due to the plant nursery trade. The study suggests that human-facilitated introductions from various geographical origins, with numerous regions colonized through continental Italy during two distinct periods, are responsible for the recurrent entry of individuals belonging to the European lineage of T. mauritanica into the Mediterranean and Macaronesia. These findings can inform better monitoring surveys and conservation programs by identifying putative current colonization routes of alien species.

Aliens Coming by Ships: Distribution and Origins of the Ocellated Skink Populations in Peninsular Italy

Simple Summary

Commercial routes are reported as the main cause of biological invasions. Particularly, naval trade may accidentally bring several species to new areas where they are not native. This is particularly evident for coastal areas, where most biological invasions occur. In our work, we reported, for the first time, the presence of the ocellated skink, native to the largest Italian islands (Sardinia, Sicily and surrounding islets in a port area of continental Central Italy). We collected several individuals of this alien population and we sampled them for molecular analyses, comparing them with those naturally occurring in Sardinia, Sicily and the Mediterranean basin, including individuals accidentally introduced to peninsular Southern Italy. Differently from what previously suggested, the nucleus in Portici (Southern Italy) may have originated from Sardinia. The intense cork trade and touristic traffic between Sardinia and Southern Tuscany may have been responsible for the introduction of this lizard also to Central Italy.

Abstract

The ocellated skink (Chalcides ocellatus) is a widespread lizard, naturally distributed between the Maghreb and coastal Pakistan, with few insular populations in the Mediterranean coastal area. Some populations of this species have also been recorded in peninsular Italy, Campania and Southern Tuscany due to accidental introductions via touristic and commercial routes. In this work, we conducted genetic analyses on mitochondrial DNA COXI, cytb and 16S mtDNA genes on a sample of Italian insular and peninsular populations. Differently from what previously suggested, the nucleus in Portici (Southern Italy) may have originated from Sardinia. The intense trade and touristic traffic between Sardinia and Southern Tuscany may have been responsible for the introduction of this lizard also to Central Italy.

Dispersal patterns of Squamata in the Mediterranean: An evolutionary perspective

Mediterranean islands have a high diversity of squamates, although they are unevenly distributed. This variability in the composition of the reptile assemblages across islands may have been influenced by differences in the colonization abilities of these species. To evaluate the dispersal capacities of squamate species, we modeled their sea routes using cost surface models. We estimated the effects of some life-history traits and the phylogenetic signal in the characteristics of the modeled dispersal paths. We hypothesized that a significant phylogenetic signal should be present if the dispersal ability is enhanced by traits shared among evolutionarily related species. The results showed that no phylogenetic signal was present in the characteristics of the dispersal paths (i.e., in the distance traveled/bypassed sea depth). Thus, no superior island-colonizer lineages were detected in Mediterranean Squamata. However, our analyses also revealed that small-sized lizards were superior to other groups of squamates at dispersing over long distances on the sea.

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

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

Stand out from the Crowd: Small-Scale Genetic Structuring in the Endemic Sicilian Pond Turtle

The geographical pattern of genetic diversity was investigated in the endemic Sicilian pond turtle Emys trinacris across its entire distribution range, using 16 microsatellite loci. Overall, 245 specimens of E. trinacris were studied, showing high polymorphic microsatellite loci, with allele numbers ranging from 7 to 30. STRUCTURE and GENELAND analyses showed a noteworthy, geographically based structuring of the studied populations in five well-characterized clusters, supported by a moderate degree of genetic diversity (FST values between 0.075 and 0.160). Possible explanations for the genetic fragmentation observed are provided, where both natural and human-mediated habitat fragmentation of the Sicilian wetlands played a major role in this process. Finally, some conservation and management suggestions aimed at preventing the loss of genetic variability of the species are briefly reported, stressing the importance of considering the five detected clusters as independent Management Units.

Integrating hybrid zone analyses in species delimitation: lessons from two anuran radiations of the Western Mediterranean

Molecular ecologists often rely on phylogenetic evidence for assessing the species-level systematics of newly discovered lineages. Alternatively, the extent of introgression at phylogeographic transitions can provide a more direct test to assign candidate taxa into subspecies or species categories. Here, we compared phylogenetic versus hybrid zone approaches of species delimitation in two groups of frogs from the Western Mediterranean region (Discoglossus and Pelodytes), by using genomic data (ddRAD). In both genera, coalescent analyses recovered almost all nominal taxa as "species". However, the least-diverged pairs D. g. galganoi/jeanneae and P. punctatus/hespericus admix over hundreds of kilometers, suggesting that they have not yet developed strong reproductive isolation and should be treated as conspecifics. In contrast, the comparatively older D. scovazzi/pictus and P. atlanticus/ibericus form narrow contact zones, consistent with species distinctiveness. Due to their complementarity, we recommend taxonomists to combine phylogenomics with hybrid zone analyses to scale the gray zone of speciation, i.e., the evolutionary window separating widely admixing lineages versus nascent reproductively isolated species. The radically different transitions documented here conform to the view that genetic incompatibilities accumulating with divergence generate a weak barrier to gene flow for long periods of time, until their effects multiply and the speciation process then advances rapidly. Given the variability of the gray zone among taxonomic groups, at least from our current abilities to measure it, we recommend to customize divergence thresholds within radiations to categorize lineages for which no direct test of speciation is possible.

A combination of long term fragmentation and glacial persistence drove the evolutionary history of the Italian wall lizard Podarcis siculus

BACKGROUND

The current distribution of genetic diversity is the result of a vast array of microevolutionary processes, including short-term demographic and ecological mechanisms and long-term allopatric isolation in response to Quaternary climatic fluctuations. We investigated past processes that drove the population differentiation and spatial genetic distribution of the Italian wall lizard Podarcis siculus by means of sequences of mitochondrial cytb (n = 277 from 115 localities) and nuclear mc1r and β-fibint7genes (n = 262 and n = 91, respectively) from all its distribution range. The pattern emerging from the genetic data was compared with current and past (last glacial maximum) species distribution modeling (SDM).

RESULTS

We identified seven deeply divergent parapatric clades which presumably remained isolated in different refugia scattered mainly throughout the Tyrrhenian coast. Conversely, the Adriatic coast showed only two haplogroups with low genetic variability. These results appear to agree with the SDM prediction at the last glacial maximum (LGM) indicating a narrow area of habitat suitability along the Tyrrhenian coast and much lower suitability along the Adriatic one. However, the considerable land exposure of the Adriatic coastline favored a glacial colonization of the Balkan Peninsula.

CONCLUSIONS

Our population-level historical demography showed a common trend consistent with glacial expansions and regional persistence during the last glacial maximum. This complex genetic signature appears to be inconsistent with the expectation of the expansion-contraction model and post-LGM (re)colonizations from southern refugia. Hence it is one of an increasing number of cases in which these assumptions are not met, indicating that long-term fragmentation and pre-LGM events such as glacial persistence were more prominent in shaping genetic variation in this temperate species.