Multilocus genetic diversity and historical biogeography of the endemic wall lizard from Ibiza and Formentera,Podarcis pityusensis(Squamata: Lacertidae)

DNA metabarcoding the diet of Podarcis lizards endemic to the Balearic Islands

Dietary studies are essential to unravel the functioning of ecosystems and ultimately to understand biodiversity. This task, which at first may seem simple, becomes especially complex in those cases of omnivorous species with highly variable diets. In this regard, the emergence of next-generation DNA sequencing methodologies represents a powerful tool to address the problem. Here we implement a high-throughput metabarcoding strategy based on the analysis of four molecular markers aimed at sequencing both mitochondrial (animal prey) and chloroplast (diet plants) genome fragments from fecal samples of two lizard species endemic to the Balearic Archipelago (Podarcis lilfordi and P. pityusensis) obtained through non-invasive methods. The results allowed for the characterization of their diets with a high degree of taxonomic detail and have contributed a large number of new trophic records. The reported diets are based mainly on the consumption of arthropods, mollusks and plants from a diversity of taxonomic orders, as well as carrion and marine subsidies. Our analyses also reveal inter- and intra-specific differences both in terms of seasonality and geographical distribution of the sampled lizard populations. These molecular findings provide new insights into the trophic interactions of these threatened endemic lizards in their unique and isolated ecosystems.

Faecal Microbiota Divergence in Allopatric Populations of Podarcis lilfordi and P. pityusensis, Two Lizard Species Endemic to the Balearic Islands

Gut microbial communities provide essential functions to their hosts and are known to influence both their ecology and evolution. However, our knowledge of these complex associations is still very limited in reptiles. Here we report the 16S rRNA gene faecal microbiota profiles of two lizard species endemic to the Balearic archipelago (Podarcis lilfordi and P. pityusensis), encompassing their allopatric range of distribution through a noninvasive sampling, as an alternative to previous studies that implied killing specimens of these IUCN endangered and near-threatened species, respectively. Both lizard species showed a faecal microbiome composition consistent with their omnivorous trophic ecology, with a high representation of cellulolytic bacteria taxa. We also identified species-specific core microbiota signatures and retrieved lizard species, islet ascription, and seasonality as the main factors in explaining bacterial community composition. The different Balearic Podarcis populations are characterised by harbouring a high proportion of unique bacterial taxa, thus reinforcing their view as unique and divergent evolutionary entities.

Fossil-calibrated time tree of Podarcis wall lizards provides limited support for biogeographic calibration models

Podarcis wall lizards are endemic to the Mediterranean Basin where they represent the predominant reptile group. Despite being extensively used as model organisms in evolutionary and ecological studies their phylogeny and historical biogeography are still incompletely resolved. Moreover, molecular clock calibrations used in wall lizard phylogeography are based on the assumption of vicariant speciation triggered by the abrupt Mediterranean Sea level rise at the end of the Messinian salinity crisis (MSC). However, the validity of this biogeographic calibration remains untested. In this study we inferred a robust time tree based on multilocus data and fossil calibrations using both gene concatenation and species-tree approaches and including models with gene-flow. We found five deeply divergent, geographically coherent, and well-supported clades comprising species from i) Iberian Peninsula and North Africa; ii) Western Mediterranean islands, iii) Sicilian and Maltese islands; and iv-v) Balkan region and Aegean islands. The mitochondrial tree shows some inconsistencies with the species tree that warrant future investigation. Diversification of main clades is estimated in a short time frame during the Middle Miocene and might have been associated with a period of global climate cooling with the establishment of a marked climatic zonation in Europe. Cladogenetic events within the main clades are scattered throughout the time tree, from the Late Miocene to the Early Pleistocene, suggesting that speciation events in wall lizards reflect a complex interplay between regional topography, climate and geological history rather than a shared major climatic or paleogeographic event. Our absolute time estimates, as well as a relative dating approach, demonstrate that the assumption of a causal link between sea-level rise at the end of the MSC and the diversification of many island endemics is not justified. This study reinforces the notion that multiple dispersal and vicariant events, at different time frames, are required to explain current allopatric distributions and to account for the historical assembly of Mediterranean biota, and cautions against the use of biogeographic calibrations based on the assumption of vicariance.

Whole Plastome Sequencing Within Silene Section Psammophilae Reveals Mainland Hybridization and Divergence With the Balearic Island Populations

Reconstructing the phylogenetic relationships within Caryophyllaceae tribe Sileneae has been obscured by hybridization and incomplete lineage sorting. Silene is the largest genus in the Caryophyllaceae, and unraveling its evolutionary history has been particularly challenging. In order to infer the phylogenetic relationships among the five species in Silene section Psammophilae, we have performed a genome skimming approach to acquire the complete plastid genome (cpDNA), nuclear ribosomal cistron (nrDNA), and partial mitochondrial genome (mtDNA). We have included 26 populations, representing the range of each species' distribution. This section includes five morphologically similar species endemic to the Iberian Peninsula and Balearic Islands (Ibiza and Formentera), yet some of them occupy distinct edaphic habitats (e.g. maritime sands, calcareous sandstones). In addition to phylogeographic analyses, genetic structuring using the chloroplast data set was inferred with Discriminant Analysis of Principal Components (DAPC), analyses of molecular variance (AMOVA), and a partial Mantel test. Reference-guided assembly of 50 bp single-end and 250 bp paired-end Illumina reads produced the nearly complete cpDNA genome (154 kbp), partial mtDNA genome (from 81 to 114 kbp), and the nrDNA cistron (6.4 kbp). Selected variable regions of the cpDNA and mtDNA assemblies were confirmed by Sanger sequencing. Phylogenetic analyses of the mainland populations reveal incongruence among the three genomes. None of the three data sets produced relationships consistent with taxonomy or geography. In contrast, Silene cambessedesii, present in the Balearic Islands, is the only species that forms a strongly supported monophyletic clade in the cpDNA genome and is strongly differentiated with respect to the remaining taxa of the Iberian Peninsula. These results contrast with those obtained for mainland populations. Across the entire analysis, only one well-supported mainland clade of Silene littorea and Silene stockenii emerges from the southern region of the Iberian Peninsula. DAPC and AMOVA results suggest the absence of genetic structure among mainland populations of Silene section Psammophilae, whereas partial Mantel test discarded spatial correlation of genetic differentiation. The widespread incongruence between morphology-based taxonomic boundaries and phylogeography suggests a history of interspecific hybridization, in which only a substantial geographic barrier, like isolation by the Mediterranean Sea, was sufficient to create and maintain species boundaries in Silene section Psammophilae.

Spotlight on islands: on the origin and diversification of an ancient lineage of the Italian wall lizard Podarcis siculus in the western Pontine Islands

Groups of proximate continental islands may conceal more tangled phylogeographic patterns than oceanic archipelagos as a consequence of repeated sea level changes, which allow populations to experience gene flow during periods of low sea level stands and isolation by vicariant mechanisms during periods of high sea level stands. Here, we describe for the first time an ancient and diverging lineage of the Italian wall lizard Podarcis siculus from the western Pontine Islands. We used nuclear and mitochondrial DNA sequences of 156 individuals with the aim of unraveling their phylogenetic position, while microsatellite loci were used to test several a priori insular biogeographic models of migration with empirical data. Our results suggest that the western Pontine populations colonized the islands early during their Pliocene volcanic formation, while populations from the eastern Pontine Islands seem to have been introduced recently. The inter-island genetic makeup indicates an important role of historical migration, probably due to glacial land bridges connecting islands followed by a recent vicariant mechanism of isolation. Moreover, the most supported migration model predicted higher gene flow among islands which are geographically arranged in parallel. Considering the threatened status of small insular endemic populations, we suggest this new evolutionarily independent unit be given priority in conservation efforts.

Evolutionary history of Podarcis tiliguerta on Corsica and Sardinia

Background

Podarcis tiliguerta is a wall lizard endemic to the Mediterranean islands of Corsica and Sardinia. Previous findings of high mtDNA and morphological diversity have led to the suggestion that it may represent a species complex. Here, we analysed mitochondrial and nuclear markers (mtDNA, 3110 bp; 6 nDNA loci, 3961 bp) in P. tiliguerta sampled from thirty-two localities across Corsica and Sardinia.

Results

We find much greater intraspecific genetic divergence than between sister species of other Mediterranean island Podarcis, i.e., between P. lilfordi and P. pityusensis. We detected three mtDNA clusters in Corsica (North, South-East and South-West) and either two or three in Sardinia (North vs. South) depending on the clustering method. Only one or two nDNA groups were identified within each main island (again, depending on the method). A Bayesian time-calibrated multispecies coalescent tree was obtained from mtDNA and provided statistical support for a Miocene origin of the species (13.87 Ma, 95% HPD: 18.30–10.77 Ma). The posterior mean divergence time for the Corsican and Sardinian lineages was 12.75 Ma ago (95% HPD: 16.94–9.04 Ma).

Conclusion

The results support the evolutionary distinctiveness of Corsican and Sardinian populations and also indicate a lack of post-divergence migration despite periods of contact being possible. Further to this, species delimitation analyses of Corsican and Sardinian lineages provided statistical support for their recognition as distinct (sister) taxa. Our results provide new insights into the biogeography of the Mediterranean biodiversity hotspot, and contribute important findings relevant to the systematics and evolution of this speciose lizard genus.

Electronic supplementary material

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

An investigation for population maintenance mechanism in a miniature garden: genetic connectivity or independence of small islet populations of the Ryukyu five-lined skink

The Ryukyu five-lined skink (Plestiodon marginatus) is an island lizard that is even found in tiny islets with less than half a hectare of habitat area. We hypothesized that the island populations are maintained under frequent gene flow among the islands or independent of each other. To test our hypotheses, we investigated genetic structure of 21 populations from 11 land-bridge islands that were connected during the latest glacial age, and 4 isolated islands. Analyses using mitochondrial cytochrome b gene sequence (n = 67) and 10 microsatellite loci (n = 235) revealed moderate to high levels of genetic differentiation, existence of many private alleles/haplotypes in most islands, little contemporary migration, a positive correlation between genetic variability and island area, and a negative correlation between relatedness and island area. These evidences suggest a strong effect of independent genetic drift as opposed to gene flow, favoring the isolation hypothesis even in tiny islet populations. Isolation-by-distance effect was demonstrated and it became more prominent when the 4 isolated islands were excluded, suggesting that the pattern is a remnant of the land-bridge age. In a few island populations, however, the possibility of occasional overwater dispersals was partially supported and therefore could not be ruled out.