Effects of fragmentation on genetic diversity in island populations of the Aegean wall lizard Podarcis erhardii (Lacertidae, Reptilia)

What Influences the Prevalence and Intensity of Haemoparasites and Ectoparasites in an Insular Lizard?

Island biogeography theories predict that characteristics such as island size, age, and isolation interplay in host-parasite dynamics. In this study, we analyzed haemogregarines of the Aegean wall lizard, Podarcis erhardii, to investigate how island characteristics relate to parasite prevalence and intensity. A previous assessment of 19 Greek island populations suggested that isolation time and host population density were key predictors of haemogregarines. Here, by combining microscopy and genetic techniques, we extend this previous study to four additional islands: Syros, Folegandros, Santorini and Nea Kameni. We also recorded the prevalence of ticks and mites, definitive hosts for these parasites. The genetically identified haemogregarines are part of a clade with parasites from other lizard species, including some considered as Karyolysus, but others assigned to Hepatozoon. The prevalence of these parasites differed significantly between islands, while their intensity did not. The presence of ticks was associated with endoparasite prevalence, and males were more frequently infected by haemogregarines than females. Combining our data with that of the previous study, we found no significant impact of the island age and area on parasite prevalence. We also confirmed the presence of the unrelated parasite genus Schellackia through microscopy and DNA sequencing, which is the first record of this genus in this host species. Our results further highlight the complexity of host-parasite systems.

Eat or be eaten? An observation of Podarcis erhardii consuming Scolopendra cingulata from Andros Island, Cyclades, Greece

Podarcis wall lizards mainly feed on coleopterans, orthopterans, arachnids, and other small invertebrates. However, Aegean wall lizards, Podarcis erhardii, are widely distributed across Aegean islands and are increasingly observed eating non-traditional food items ranging from plant material to conspecific eggs and body parts. Here, we report the first documented case of P. erhardii consuming a large centipede, Scolopendra cingulata. The predator-prey relationship between these species has appeared to go both ways and may intensify on islands.

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.

Low diversity, little genetic structure but no inbreeding in a high-density island endemic pit-viper Gloydius shedaoensis

Island species and their ecosystems play an important role in global biodiversity preservation, and many vulnerable island species are conservation priorities. Although insular habitat likely facilitates the species diversification process, it may also aggravate the fragility of these species with high risk of inbreeding. The Shedao pit-viper Gloydius shedaoensis is an island endemic species with an extremely high population density, which has been categorized as vulnerable in the IUCN (International Union for the Conservation of Nature and Natural Resources) Red List. We collected 13,148 SNP (Single Nucleotide Polymorphism) from across its genome and examined its genetic diversity and demographic history. The Shedao pit-viper has a low genetic diversity but shows no sign of inbreeding. Furthermore, population genetic structure analysis, including the neighbor-joining tree, principal coordinate analysis, clustering, and spatial autocorrelation, revealed a general lack of spatial structure. Only the isolation by distance residues suggested a weak patchiness. Overall, the population is nearly panmictic and gene flow is evenly distributed across the island. A large number of individuals, small size of the island, and the lack of population structure likely all contribute to the lack of inbreeding in this species. We also detected signs of male-biased dispersal, which likely is another inbreeding avoidance strategy. Historical demographic analysis suggested that the historical population size and distribution of the species are much larger than their current ones. The multiple transgressive–regressive events since the Late Pleistocene are likely the main cause of the population size changes. Taken together, our results provide a basic scientific foundation for the conservation of this interesting and important species.

Taxonomic reassessment and phylogenetic placement of Cyrtodactylus phuketensis (Reptilia, Gekkonidae) based on morphological and molecular evidence

The taxonomy and phylogeny of the Cyrtodactyluspulchellus complex along the Thai-Malay Peninsular region has been the focus of many recent studies and has resulted in the recognition of 17 species. However, the majority of these studies were focused on Peninsular and insular Malaysia where there were specimens and genetic vouchers. The taxonomic status and phylogenetic relationships of the Thai species in this complex remain unresolved, due to the lack of genetic material of some species, especially C.phuketensis and C.macrotuberculatus from Thai populations. In this study, we investigated the phylogenetic relationship between C.phuketensis and its closely related species C.macrotuberculatus, using both morphometric and molecular data. Phylogenetic analyses of mitochondrial NADH dehydrogenase subunit 2 (ND2) gene revealed that C.phuketensis is embedded within a C.macrotuberculatus clade with 1.45–4.20% (mean 2.63%) uncorrected pairwise sequence divergences. Morphological comparisons showed nearly identical measurements of C.phuketensis and C.macrotuberculatus and overlapping ranges in meristic characters. Based on these data, C.phuketensis is considered to be a variant of C.macrotuberculatus, thus rendering C.phuketensis a junior synonym of C.macrotuberculatus.

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.

Morphological and ecological divergence in two populations of European glass lizard, Pseudopus apodus (Squamata: Anguidae)

The European glass lizard, Pseudopus apodus (Pallas, 1775), is a large, legless lizard with wide distribution across south-eastern Europe and eastern and central Asia. To date, morphological diversification among populations on a geographically small scale has not yet been reported in this lizard. Thus, we investigated the morphological variations and corresponding differences in habitat utilization between two populations of P. apodus inhabiting the same biogeographical zone within a relatively close geographic area. We hypothesized that minor differences in habitat could induce a significant level of morphological differentiation, thus indicating morphological plasticity in this species on a small geographical scale. We sampled 164 individuals (92 from the Croatian mainland and 72 from the island of Cres). Results showed that P. apodus indeed exhibited morphological differences between populations in the same biogeographical zone within a relatively close geographic area, with the Cres Island individuals being generally larger than the individuals from the mainland. Some ecological characteristics were similar in both populations (e.g., soil temperature, distance to hiding place), whereas others were distinct (e.g., air temperature and humidity). In addition, vegetation cover differed between the two sites, with more vegetation present on the mainland than on the island. Furthermore, the Cres Island population showed clear sexual dimorphism, which was absent in the mainland population.

Isolation and predation drive gecko life-history evolution on islands

Insular animals are thought to be under weak predation pressure and increased intraspecific competition compared with those on the mainland. Thus, insular populations are predicted to evolve ‘slow’ life histories characterized by fewer and smaller clutches of larger eggs, a pattern called the ‘island syndrome’. To test this pattern, we collected data on egg volume, clutch size and laying frequency of 31 Aegean Island populations of the closely related geckos of the Mediodactylus kotschyi species complex. We tested how predation pressure, resource abundance, island area and isolation influenced reproductive traits. Isolation and predation were the main drivers of variation in life-history traits. Higher predator richness seemed to promote faster life histories, perhaps owing to predation on adults, whereas the presence of boas promoted slower life histories, perhaps owing to release from predation by rats on the eggs of geckos. Insular geckos followed only some of the predictions of the ‘island syndrome’. Predation pressure seemed to be more complex than expected and drove life histories of species in two opposing directions. Our results highlight the importance of considering the identity of specific predators in ecological studies.

Island biogeography theory explains the genetic diversity of a fragmented rock ptarmigan (Lagopus muta) population

The island biogeography theory is one of the major theories in ecology, and its applicability to natural systems is well documented. The core model of the theory, the equilibrium model of island biogeography, predicts that species diversity on an island is positively related to the size of the island, but negatively related by the island's distance to the mainland. In recent years, ecologists have begun to apply this model when investigating genetic diversity, arguing that genetic and species diversity might be influenced by similar ecological processes. However, most studies have focused on oceanic islands, but knowledge on how the theory applies to islands located on the mainland (e.g., mountain islands, forest islands) is scarce. In this study, we examined how the size and degree of isolation of mountain islands would affect the genetic diversity of an alpine bird, the rock ptarmigan (Lagopus muta). Within our study area, we defined the largest contiguous mountain area as the mainland, while smaller mountains surrounding the mainland were defined as islands. We found that the observed heterozygosity (H o) was significantly higher, and the inbreeding coefficient (F is) significantly lower, on the mainland compared to islands. There was a positive significant relationship between the unbiased expected heterozygosity (H n.b.) and island size (log km2), but a negative significant relationship between H o and the cost distance to the mainland. Our results are consistent with the equilibrium model of island biogeography and show that the model is well suited for investigating genetic diversity among islands, but also on islands located on the mainland.

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.

The effects of historical fragmentation on major histocompatibility complex class II β and microsatellite variation in the Aegean island reptile, Podarcis erhardii

The major histocompatibility complex (MHC) plays a key role in disease resistance and is the most polymorphic gene region in vertebrates. Although habitat fragmentation is predicted to lead to a loss in MHC variation through drift, the impact of other evolutionary forces may counter this effect. Here we assess the impact of selection, drift, migration, and recombination on MHC class II and microsatellite variability in 14 island populations of the Aegean wall lizard Podarcis erhardii. Lizards were sampled from islands within the Cyclades (Greece) formed by rising sea levels as the last glacial maximum approximately 20,000 before present. Bathymetric data were used to determine the area and age of each island, allowing us to infer the corresponding magnitude and timing of genetic bottlenecks associated with island formation. Both MHC and microsatellite variation were positively associated with island area, supporting the hypothesis that drift governs neutral and adaptive variation in this system. However, MHC but not microsatellite variability declined significantly with island age. This discrepancy is likely due to the fact that microsatellites attain mutation‐drift equilibrium more rapidly than MHC. Although we detected signals of balancing selection, recombination and migration, the effects of these evolutionary processes appeared negligible relative to drift. This study demonstrates how land bridge islands can provide novel insights into the impact of historical fragmentation on genetic diversity as well as help disentangle the effects of different evolutionary forces on neutral and adaptive diversity.

Importance of dispersal routes that minimize open-ocean movement to the genetic structure of island populations

Islands present a unique scenario in conservation biology, offering refuge yet imposing limitations on insular populations. The Kimberley region of northwestern Australia has more than 2500 islands that have recently come into focus as substantial conservation resources. It is therefore of great interest for managers to understand the driving forces of genetic structure of species within these island archipelagos. We used the ubiquitous bar-shouldered skink (Ctenotus inornatus) as a model species to represent the influence of landscape factors on genetic structure across the Kimberley islands. On 41 islands and 4 mainland locations in a remote area of Australia, we genotyped individuals across 18 nuclear (microsatellite) markers. Measures of genetic differentiation and diversity were used in two complementary analyses. We used circuit theory and Mantel tests to examine the influence of the landscape matrix on population connectivity and linear regression and model selection based on Akaike's information criterion to investigate landscape controls on genetic diversity. Genetic differentiation between islands was best predicted with circuit-theory models that accounted for the large difference in resistance to dispersal between land and ocean. In contrast, straight-line distances were unrelated to either resistance distances or genetic differentiation. Instead, connectivity was determined by island-hopping routes that allow organisms to minimize the distance of difficult ocean passages. Island populations of C. inornatus retained varying degrees of genetic diversity (NA = 1.83 - 7.39), but it was greatest on islands closer to the mainland, in terms of resistance-distance units. In contrast, genetic diversity was unrelated to island size. Our results highlight the potential for islands to contribute to both theoretical and applied conservation, provide strong evidence of the driving forces of population structure within undisturbed landscapes, and identify the islands most valuable for conservation based on their contributions to gene flow and genetic diversity.

Intraspecific Colour Variation among Lizards in Distinct Island Environments Enhances Local Camouflage

Within-species colour variation is widespread among animals. Understanding how this arises can elucidate evolutionary mechanisms, such as those underlying reproductive isolation and speciation. Here, we investigated whether five island populations of Aegean wall lizards (Podarcis erhardii) have more effective camouflage against their own (local) island substrates than against other (non-local) island substrates to avian predators, and whether this was linked to island differences in substrate appearance. We also investigated whether degree of local substrate matching varied among island populations and between sexes. In most populations, both sexes were better matched against local backgrounds than against non-local backgrounds, particularly in terms of luminance (perceived lightness), which usually occurred when local and non-local backgrounds were different in appearance. This was found even between island populations that historically had a land connection and in populations that have been isolated relatively recently, suggesting that isolation in these distinct island environments has been sufficient to cause enhanced local background matching, sometimes on a rapid evolutionary time-scale. However, heightened local matching was poorer in populations inhabiting more variable and unstable environments with a prolonged history of volcanic activity. Overall, these results show that lizard coloration is tuned to provide camouflage in local environments, either due to genetic adaptation or changes during development. Yet, the occurrence and extent of selection for local matching may depend on specific conditions associated with local ecology and biogeographic history. These results emphasize how anti-predator adaptations to different environments can drive divergence within a species, which may contribute to reproductive isolation among populations and lead to ecological speciation.

Influence of gene flow on divergence dating - implications for the speciation history of Takydromus grass lizards

Dating the time of divergence and understanding speciation processes are central to the study of the evolutionary history of organisms but are notoriously difficult. The difficulty is largely rooted in variations in the ancestral population size or in the genealogy variation across loci. To depict the speciation processes and divergence histories of three monophyletic Takydromus species endemic to Taiwan, we sequenced 20 nuclear loci and combined with one mitochondrial locus published in GenBank. They were analysed by a multispecies coalescent approach within a Bayesian framework. Divergence dating based on the gene tree approach showed high variation among loci, and the divergence was estimated at an earlier date than when derived by the species-tree approach. To test whether variations in the ancestral population size accounted for the majority of this variation, we conducted computer inferences using isolation-with-migration (IM) and approximate Bayesian computation (ABC) frameworks. The results revealed that gene flow during the early stage of speciation was strongly favoured over the isolation model, and the initiation of the speciation process was far earlier than the dates estimated by gene- and species-based divergence dating. Due to their limited dispersal ability, it is suggested that geographical isolation may have played a major role in the divergence of these Takydromus species. Nevertheless, this study reveals a more complex situation and demonstrates that gene flow during the speciation process cannot be overlooked and may have a great impact on divergence dating. By using multilocus data and incorporating Bayesian coalescence approaches, we provide a more biologically realistic framework for delineating the divergence history of Takydromus.

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.

Landscape models for nuclear genetic diversity and genetic structure in white-footed mice (Peromyscus leucopus)

Dramatic changes in the North American landscape over the last 12 000 years have shaped the genomes of the small mammals, such as the white-footed mouse (Peromyscus leucopus), which currently inhabit the region. However, very recent interactions of populations with each other and the environment are expected to leave the most pronounced signature on rapidly evolving nuclear microsatellite loci. We analyzed landscape characteristics and microsatellite markers of P. leucopus populations along a transect from southern Ohio to northern Michigan, in order to evaluate hypotheses about the spatial distribution of genetic heterogeneity. Genetic diversity increased to the north and was best approximated by a single-variable model based on habitat availability within a 0.5-km radius of trapping sites. Interpopulation differentiation measured by clustering analysis was highly variable and not significantly related to latitude or habitat availability. Interpopulation differentiation measured as FST values and chord distance was correlated with the proportion of habitat intervening, but was best explained by agricultural distance and by latitude. The observed gradients in diversity and interpopulation differentiation were consistent with recent habitat availability being the major constraint on effective population size in this system, and contradicted the predictions of both the postglacial expansion and core-periphery hypotheses.

Genetic effects of recent habitat fragmentation in the Thousand-Island Lake region of southeast China on the distylous herb Hedyotis chrysotricha (Rubiaceae)

PREMISE OF THE STUDY

Known-age artificial-lake islands provide ideal model systems to elucidate the genetic and evolutionary consequences of anthropogenic habitat fragmentation on very recent time scales. Here, we studied a distylous herb, Hedyotis chrysotricha (Rubiaceae), in the artificially created Thousand-Island Lake (TIL) region of southeast China to explore the genetic consequences of islanding for this species. •

METHODS

Seven microsatellite loci were used to genotype 384 individuals of H. chrysotricha from 18 populations to estimate genetic diversity, population structure, and demographic parameters. •

KEY RESULTS

Island populations had significantly lower mean genetic diversity than those from the western/eastern mainland (e.g., H(E) = 0.381 vs. 0.461) and also displayed higher mean subdivision (F(ST) = 0.12 vs. 0.042/0.051). BayesAss analyses indicated moderate levels of migration rates among most populations, whereas Bottleneck did not provide strong evidence for such effects. In consequence, 2MOD strongly favored a gene flow-drift model over a pure drift model in the study area, but concomitantly revealed a relatively greater influence of drift in the island populations as evidenced by their significantly higher probabilities of allelic coancestry (F = 0.184 vs. 0.085). •

CONCLUSIONS

The observed genetic patterns in H. chrysotricha indicate that recent anthropogenic habitat fragmentation in the TIL region can lead to significant loss of genetic diversity in isolated fragments (islands) due to ongoing drift. By contrast, patterns of random mating, gene flow, and population connectivity have not greatly been modified yet, possibly owing to the species' fruit (seed) dispersal capabilities providing resilience in the face of habitat fragmentation.

Vicariance divergence and gene flow among islet populations of an endemic lizard

Allopatry and allopatric speciation can arise through two different mechanisms: vicariance or colonization through dispersal. Distinguishing between these different allopatric mechanisms is difficult and one of the major challenges in biogeographical research. Here, we address whether allopatric isolation in an endemic island lizard is the result of vicariance or dispersal. We estimated the amount and direction of gene flow during the divergence of isolated islet populations and subspecies of the endemic Skyros wall lizard Podarcis gaigeae, a phenotypically variable species that inhabits a major island and small islets in the Greek archipelago. We applied isolation-with-migration models to estimate population divergence times, population sizes and gene flow between islet-mainland population pairs. Divergence times were significantly correlated with independently estimated geological divergence times. This correlation strongly supports a vicariance scenario where islet populations have sequentially become isolated from the major island. We did not find evidence for significant gene flow within P. g. gaigeae. However, gene-flow estimates from the islet to the mainland populations were positively affected by islet area and negatively by distance between the islet and mainland. We also found evidence for gene flow from one subspecies (P. g. weigandi) into another (P. g. gaigeae), but not in the other direction. Ongoing gene flow between the subspecies suggests that even in this geographically allopatric scenario with the sea posing a strong barrier to dispersal, divergence with some gene flow is still feasible.

Population divergence in chemical signals and the potential for premating isolation between islet- and mainland populations of the Skyros wall lizard (Podarcis gaigeae)

When sexually selected traits diverge because of different local selective environments, premating isolation might arise as a correlated response. However, sexually selected traits might also diverge by stochastic forces. Here, we show that odour-based mate preferences and scent composition have diverged between islet- and mainland populations of Skyros wall lizard, Podarcis gaigeae. We quantified the degree of scent-mediated premating isolation between populations. Islet lizards preferred scent from islet lizards, whereas the mainland populations were less discriminatory. The pheromone compositions differed more between islets than between islet- and mainland populations and did not differ significantly between mainland populations. There was a tendency for population divergence in pheromones to be positively correlated with neutral genetic divergence. This might indicate a role for genetic drift in evolutionary change in these signals and partial decoupling between signals and preferences. Our results suggest that chemical signals and associated mate preferences can diverge through stochastic and selective forces and influence premating isolation.

Climate change and elevated extinction rates of reptiles from Mediterranean Islands

Recent climate change has caused the distributions of many species to shift poleward, yet few empirical studies have addressed which species will be vulnerable to longer-term climate changes. To investigate past consequences of climate change, we calculated the population extinction rates of 35 reptile species from 87 Greek land-bridge islands in the Mediterranean that occurred over the past 16,000 years. Population extinction rates were higher for those species that today have more northern distributions. We further found that northern species requiring cool, mesic habitats had less available suitable habitat among islands, implicating loss of suitable habitat in their elevated extinction rates. These extinctions occurred in the context of increasing habitat fragmentation, with islands shrinking and separating as sea levels rose. Thus, the circumstances faced by reptiles on the islands are similar to challenges for numerous species today that must cope with a changing climate while living in an increasingly human-fragmented landscape. Our island-biogeographical approach to investigating historical population extinctions gives insight into the long-term patterns of species responses to climate change.