Population genetics and phylogeographic history of the insular lizard Podarcis lilfordi (Gunther, 1874) from the Balearic Islands based on genome-wide polymorphic data

Islands provide a great system to explore the processes that maintain genetic diversity and promote local adaptation. We explored the genomic diversity of the Balearic lizard Podarcis lilfordi, an endemic species characterized by numerous small insular populations with large phenotypic diversity. Using the newly available genome for this species, we characterized more than 300,000 SNPs, merging genotyping-by-sequencing (GBS) data with previously published restriction site-associated DNA sequencing (RAD-Seq) data, providing a dataset of 16 island populations (191 individuals) across the range of species distribution (Menorca, Mallorca, and Cabrera). Results indicate that each islet hosts a well-differentiated population (F ST = 0.247 ± 0.09), with no recent immigration/translocation events. Contrary to expectations, most populations harbor a considerable genetic diversity (mean nucleotide diversity, P i = 0.144 ± 0.021), characterized by overall low inbreeding values (F IS < 0.1). While the genetic diversity significantly decreased with decreasing islet surface, maintenance of substantial genetic diversity even in tiny islets suggests variable selection or other mechanisms that buffer genetic drift. Maximum-likelihood tree based on concatenated SNP data confirmed the existence of the two major independent lineages of Menorca and Mallorca/Cabrera. Multiple lines of evidence, including admixture and root testing, robustly placed the origin of the species in the Mallorca Island, rather than in Menorca. Outlier analysis mainly retrieved a strong signature of genome differentiation between the two major archipelagos, especially in the sexual chromosome Z. A set of proteins were target of multiple outliers and primarily associated with binding and catalytic activity, providing interesting candidates for future selection studies. This study provides the framework to explore crucial aspects of the genetic basis of phenotypic divergence and insular adaptation.

Life span, growth, senescence and island syndrome: Accounting for imperfect detection and continuous growth

Small vertebrates on islands are expected to attain a larger body size, and a greater survival than their mainland counterparts. Comparative studies have questioned whether lizards exhibit this set of adaptations, referred to as the 'island syndrome'. We collected data on 730 individuals the endemic Lilford's lizard Podarcis lilfordi throughout a 10-year period on a small island of the Balearic archipelago (Spain). We coupled a growth function with a capture-mark-recapture model to simultaneously estimate size- and sex-dependent growth rate and survival. To put our results into a wider context, we conducted a systematic review of growth, life span and age at maturity in different Podarcis species comparing insular and mainland populations. We found a low average growth coefficient (0.56 and 0.41 year^-1 for males and females to reach an asymptotic size of 72.3 and 65.6 mm respectively), a high annual survival probability of 0.81 and 0.79 in males and females, and a large variability between individuals in growth parameters. Survival probability decreased with body size in both sexes, indicating a senescence pattern typical of long-lived species or in populations with a low extrinsic mortality. Assuming a constant survival after sexual maturity, at about 2 years old, the average life span was 6.18 years in males and 8.99 in females. The oldest animal was a male last captured at an estimated age of ≥13 years and still alive at the end of the study. Our results agree with the predictions of the 'island syndrome' for survival, life span and growth parameters. A comparative analysis of these values across 29 populations of 16 different species of Podarcis indicated that insular lizards grow slower and live longer than their mainland counterparts. However, our data differed from other island populations of the same species, suggesting that island-specific characteristics play an additional role to isolation. Within this study we developed an analytical approach to study the body size-dependent survival of small reptiles. We discuss its applicability to contrast hypotheses on senescence in different sexes of this species, and provide the code used to integrate the growth and capture-mark-recapture models.

Agent-based models reveal limits of mark-release-recapture estimates for the rare butterfly, Bhutanitis thaidina (Lepidoptera: Papilionidae)

Insect diversity and abundance are in drastic decline worldwide, but quantifying insect populations to better conserve them is a difficult task. Mark-release-recapture (MRR) is widely used as an ecological indicator for insect populations, but the accuracy of MRR estimates can vary with factors such as spatial scale, sampling effort and models of inference. We conducted a 3-year MRR study of B. thaidina in Yanzigou valley, Mt. Gongga but failed to obtain sufficient data for a robust population estimate. This prompted us to integrate B. thaidina life history information to parameterize agent-based models and evaluate the conditions under which successful MRR studies could be conducted. We evaluated: (1) the performance of MRR models under different landscape types, and (2) the influence of experimental design on the accuracy and variance of MRR-based estimates. Our simulations revealed systematic underestimates of true population parameters by MRR models when sampling effort was insufficient. In a total of 2772 simulations, subjective decisions in sampling protocol (e.g., frequency, number of sampling locations, use of spatially explicit models, type of estimands) accounted for nearly half of the variation in estimates. We conclude that MRR-based estimates could be improved with the addition of more field-specific parameters.

Invasive snake causes massive reduction of all endemic herpetofauna on Gran Canaria

Invasive snakes represent a serious threat to island biodiversity, being responsible for far-reaching impacts that are noticeably understudied, particularly regarding native reptiles. We analysed the impact of the invasive California kingsnake, Lampropeltis californiae-recently introduced in the Canary Islands-on the abundance of all endemic herpetofauna of the island of Gran Canaria. We quantified the density in invaded and uninvaded sites for the Gran Canaria giant lizard, Gallotia stehlini, the Gran Canaria skink, Chalcides sexlineatus, and Boettger's wall gecko, Tarentola boettgeri. We used spatially explicit capture-recapture and distance-sampling methods for G. stehlini and active searches under rocks for the abundance of the other two reptiles. The abundance of all species was lower in invaded sites, with a reduction in the number of individuals greater than 90% for G. stehlini, greater than 80% for C. sexlineatus and greater than 50% for T. boettgeri in invaded sites. Our results illustrate the severe impact of L. californiae on the endemic herpetofauna of Gran Canaria and highlight the need for strengthened measures to manage this invasion. We also provide further evidence of the negative consequences of invasive snakes on island reptiles and emphasize the need for further research on this matter on islands worldwide.

Differential habitat use by sympatric species of mouse lemurs across a mangrove–dry forest habitat gradient

Understanding the mechanisms by which similar species coexist in sympatry is a major driver of ecological research. Niche partitioning and ecological plasticity can facilitate spatial and habitat use overlap between generalist and specialist species. Mouse lemurs (Microcebus spp.) are a highly speciose group of small primates that are endemic to the forests of Madagascar. In northwestern Madagascar, the relatively widespread M. murinus occurs sympatrically with the microendemic M. ravelobensis. We investigated spatial distributions and densities of these two species across a mangrove–dry forest habitat gradient in Mariarano commune. We used capture-mark-recapture techniques and nocturnal line transect surveys along six transects during June and July 2017. Spatial capture-recapture and distance sampling models were used to estimate lemur densities across habitat types. The congeners displayed differential patterns of spatial distribution and densities. Microcebus murinus was found in similar densities across all habitat types, while M. ravelobensis was found at much higher densities in dry forests compared with mangroves. This suggests that the generalist M. murinus uses a wider array of habitats more evenly than the specialist M. ravelobensis. Our study provides empirical evidence of how cryptic lemur species differ in their habitat use and distribution across an environmental gradient and provides new insights into their ecology in an understudied habitat. Lemurs are one of the most threatened groups of mammals in the world, and understanding how these species are distributed across different forest types is crucial for planning and implementing conservation measures to protect lemur habitat.

Contrasting size-dependent life history strategies of an insular lizard

In many species with continuous growth, body size is an important driver of life-history tactics and its relative importance is thought to reflect the spatio-temporal variability of selective pressures. We developed a deterministic size-dependent integral projection model for 3 insular neighboring lizard populations with contrasting adult body sizes to investigate how size-related selective pressures can influence lizard life-history tactics. For each population, we broke down differences in population growth rates into contributions from size-dependent body growth, survival, and fecundity. A life table response experiment (LTRE) was used to compare the population dynamics of the 3 populations and quantify the contributions of intrinsic demographic coefficients of each population to the population growth rate (λ). Perturbation analyses revealed that the largest adults contributed the most to the population growth rate, but this was not true in the population with the smallest adults and size-independent fertility. Although we were not able to identify a single factor responsible for this difference, the combination of the demographic model on a continuous trait coupled with an LTRE analysis revealed how individuals from sister populations of the same species follow different life strategies and showed different compensatory mechanisms among survival, individual body growth, and fertility. Our results indicate that body size can play a contrasting role even in closely-related and closely-spaced populations.

Diet composition of the lizard Podarcis lilfordi (Lacertidae) on 2 small islands: an individual-resource network approach

Despite it is widely accepted that intrapopulation variation is fundamental to ecological and evolutionary processes, this level of information has only recently been included into network analysis of species/population interactions. When done, it has revealed non-random patterns in the distribution of trophic resources. Nestedness in resource use among individuals is the most recurrent observed pattern, often accompanied by an absence of modularity, but no previous studies examine bipartite modularity. We use network analysis to describe the diet composition of the Balearic endemic lizard Podarcis lilfordi in 2 islets at population and individual levels, based on the occurrence of food items in fecal samples. Our objectives are to 1) compare niche structure at both levels, 2) characterize niche partition using nestedness and modularity, and 3) assess how size, sex, season, and spatial location influence niche structure. At population-level niche width was wide, but narrow at the level of the individual. Both islet networks were nested, indicating similar ranking of the food preferences among individuals, but also modular, which was partially explained by seasonality. Sex and body size did not notably affect diet composition. Large niche overlap and therefore possibly relaxed competition were observed among females in one of the islets and during spring on both islets. Likewise, higher modularity in autumn suggests that higher competition could lead to specialization in both populations, because resources are usually scarce in this season. The absence of spatial location influence on niche might respond to fine-grained spatio-temporally distribution of food resources. Behavioral traits, not included in this study, could also influence resource partitioning.

Hierarchical distance sampling to estimate population sizes of common lizards across a desert ecoregion

Multispecies wildlife monitoring across large geographical regions is important for effective conservation planning in response to expected impacts from climate change and land use. Unlike many species of birds, mammals, and amphibians which can be efficiently sampled using automated sensors including cameras and sound recorders, reptiles are often much more challenging to detect, in part because of their typically cryptic behavior and generally small body sizes. Although many lizard species are more active during the day which makes them easier to detect using visual encounter surveys, they may be unavailable for sampling during certain periods of the day or year due to their sensitivity to temperature. In recognition of these sampling challenges, we demonstrate application of a recent innovation in distance sampling that adjusts for temporary emigration between repeat survey visits. We used transect surveys to survey lizards at 229 sites throughout the Mojave Desert in California, USA, 2016. We estimated a total population size of 82 million (90% CI: 65-99 million) for the three most common species of lizards across this 66,830 km2 ecoregion. We mapped how density at the 1-km2 scale was predicted to vary with vegetation cover and human development. We validated these results against independent surveys from the southern portion of our study area. Our methods and results demonstrate how multispecies monitoring programs spanning arid ecoregions can better incorporate information about reptiles.

Wind farms have cascading impacts on ecosystems across trophic levels

Wind farms are a cleaner alternative to fossil fuels for mitigating the effects of climate change, but they also have complex ecological consequences. In the biodiversity hotspot of the Western Ghats in India, we find that wind farms reduce the abundance and activity of predatory birds (for example, Buteo, Butastur and Elanus species), which consequently increases the density of lizards, Sarada superba. The cascading effects of wind turbines on lizards include changes in behaviour, physiology and morphology that reflect a combination of predator release and density-dependent competition. By adding an effective trophic level to the top of food webs, we find that wind farms have emerging impacts that are greatly underestimated. There is thus a strong need for an ecosystem-wide view when aligning green-energy goals with environment protection.

Relative role of population density and climatic factors in shaping the body growth rate of Lilford’s Wall Lizard (Podarcis lilfordi)

The body growth rate in small reptiles is modulated by per-capita food resources and recent evidences suggested that this constraint is the mechanism underlying differences between cohorts. Per-capita food resources depend on population size and climatic factors, but their relative role in explaining the variations in growth rate is unclear. We used morphological data collected over 6 years to model the body growth of an insular lizard (Lilford’s Wall Lizard, Podarcis lilfordi (Günther, 1874)). We used nonlinear equations to describe the appropriate length-at-age relationship. For each sex, seasonal growth was observed and the oscillatory von Bertalanffy curve was fitted to the data. Three age classes were recognized, and we investigated the relative role of density, spring rainfall, and temperature in explaining the variation of the growth rate in each sex-by-age class. Results showed that the relative role of density and climatic factors varied according to the age considered. While population size and temperature had a negative effect on the growth rate of juveniles, rainfall had a positive influence on the growth of subadults. Adult growth was near zero and constant over time. The different role of density-dependent and climatic factors in explaining age-dependent growth rate provides an important insight in understanding lizard population dynamics and life-history tactics.