Human-induced morphological shifts in an island lizard

Understanding the evolutionary consequences of anthropogenic change is an emerging topic in evolutionary biology. While highly sensitive species may go extinct in response to anthropogenic habitat alteration, those with broader environmental tolerances may persist and adapt to the changes. Here, we use morphological data from the brown anole (Anolis sagrei), a lizard species that lives in both natural and human-disturbed habitats, to examine the impact of anthropogenic habitat alteration. We find populations inhabiting disturbed habitats were significantly larger in snout-vent length, hindspan, and mass and provide evidence that the observed divergence in hindspan is driven by human-induced changes in habitat structure. Populations were found to be genetically distinct among islands but are not genetically differentiated between habitat types on islands. Thus, the observed pattern of intra-island morphological differences cannot be explained by separate founding populations. Rather, our results are consistent with morphological differences between habitats having arisen in situ on each island. Results underscore the significant impact anthropogenic change may have on evolutionary trajectories of populations that persist in human-altered habitats.


Introduction
Two-thirds of the world's terrestrial area is now devoted directly to supporting human populations, either through agriculture, fisheries, urbanization, or infrastructure (Millenium Ecosystem Assessment 2005). Species that are able to persist in human-altered habitats are often subject to environmental conditions that are very different from those in their natural habitat. Habitat degradation can lead to shifts in community structure that in turn alter predation pressure and competitive regimes (Klein 1989;Eggleston et al. 2005;Foley et al. 2005;Laeser et al. 2005). Shifts in microclimatic conditions such as light availability and thermal characteristics are among the other effects of habitat alteration (Sumner et al. 1999). Divergent natural selection in human-altered habitats may lead to unique phenotypes not found in natural habitats (Hendry et al. 2006) and shift the evolutionary trajectories of species (Smith and Bernatchez 2008). Even if the novel traits are environmentally induced plastic responses, selection pressures acting on individuals can alter the population's evolutionary trajectory via genetic assimilation (West-Eberhard 2003). Human-induced habitat alterations can therefore have long-term effects on populations inhabiting disturbed areas. Escalating human disturbance underscores the importance of understanding the full spectrum of evolutionary impacts that humans are having on natural populations (Smith and Bernatchez 2008;Hendry et al. 2011;Sih et al. 2011).
Anolis lizards provide an exceptional opportunity to investigate the evolutionary impacts of anthropogenic habitat alteration because of the well-documented association between morphology and habitat use in this group (Williams 1983;Losos 1990;Losos et al. 1994Losos et al. , 2001Langerhans et al. 2006). The adaptive radiation of Anolis lizards is one of the most extraordinary vertebrate radiations known (Schluter 2000). Six distinct morphological classes, termed ecomorphs (Williams 1983), have evolved repeatedly and independently, giving rise to the nearly 150 species found in the Caribbean today (Jackman et al. 1999). Ecomorphs are characterized by divergence in body size and limb length, traits that are highly correlated with habitat use (Losos 1990(Losos , 1994Irschick and Losos 1998) and which are highly heritable in the species studied here (data for A. sagrei, body size, as indexed by snout-vent length (SVL): h 2 = 1.02 ± 0.29, hindlimb length: h 2 = 0.78 ± 0.13) (Calsbeek and Bonneaud 2008). Species with relatively long limbs are most often found perching on broad diameter substrates (e.g. tree trunks), whereas species with relatively shorter limbs are more often found on narrow diameter substrates (e.g. twigs). These correlations are thought to have arisen through natural selection favoring fast-moving, long-limbed lizards on broad surfaces (enhanced maximum sprint speed) and slower, more sure-footed lizards on narrow surfaces (enhanced agility) (Losos and Sinervo 1989;Losos 1990). These correlations between morphology and habitat use are considered adaptive and are thought to have played an important role in the adaptive radiation of these lizards (Losos et al. 2001). Recent work has also shown that the patterns of divergence in morphology among species are mirrored at the intraspecific level Marnocha 2008).
Anolis sagrei is a small (40-70 mm, SVL), semi-arboreal lizard with a broad tropical and subtropical distribution. It is a member of the trunk-ground ecomorph class and the most common anole in the Bahamas. On the islands of the Bahamas, A. sagrei occupies both its natural forest habitat and habitat that has been greatly altered by human activities. Here, we use morphological measurements from populations of A. sagrei to test the hypothesis that anthropogenic habitat alteration has lead to morphological divergence between lizards living in disturbed and natural habitats. In particular, we examine whether human-induced shifts in habitat structure have resulted in divergence in limb length, a key character in the diversification of Anolis lizards. Finally, we assess genetic variation among our study populations to test whether the patterns of morphological divergence are congruent with patterns of genetic divergence.

Study sites
During June and July of 2003, we captured A. sagrei from four pairs of sites (natural and disturbed) with a pair on each of the following islands: New Providence (25°04¢N, 77°22¢W), Harbour Island (25°30¢N, 76°38¢W), Eleuthera (25°15¢N, 76°18.5¢W), and Great Exuma (23°30¢N, 75°46¢W) (Fig. 1). On each island, the natural and disturbed study sites were within 5 km of each other. Natural sites consisted of densely wooded areas (termed blackland coppice) containing a great diversity of plant species and were under little or no human influence. These areas have a dense canopy typically 30-40 feet high and are composed of broadleaf hardwood trees such as Caribbean mahogany (Sweitenia mahagoni) and gumbolimbo (Bursera simaruba). The shady, humid understory supports numerous plant species, including bromeliads (Tillandsia spp.), orchids (Spiranthes spp. and Habenaria spp.), and a wide variety of shrubby plants. Disturbed habitat, by striking contrast, is highly impacted by humans and consists of largely open areas characterized by the presence of buildings and other man-made structures such as roads, sidewalks, and fences. Disturbed habitats have virtually no canopy cover and the vegetation consists primarily of palm trees (Arecaceae family), low shrubbery, and manicured lawns. Disturbed study sites are all located in towns that have been occupied by humans for at least 200-350 years (Craton 2007), the equivalent of hundreds of generations for Anolis lizards.
To quantitatively characterize habitat at each of our eight sites, we used satellite remote sensing data derived from the vegetation continuous field (Hansen et al. 2002) product from the MODIS satellite (https://lpdaac.usgs.gov/lpdaac/products/modis_overview) at 500-m resolution as a measure of the percentage of tree cover. We extracted the percent tree cover at our sampling sites using ArcGIS 9.3 (ESRI 2004) and visually verified the extracted data by means of optical satellite imagery. Although the environment in the Bahamas varies over small scales, the resolution of the data used is sufficiently high to provide an accurate characterization of the habitat available at each of the study sites.

Field methods
At each of the eight study sites, we captured approximately 20 adult male A. sagrei by hand or using a silk noose. We made the following measurements on each individual: snout-vent length (SVL; from the tip of the snout to the anterior side of the cloaca), hindspan (from one femoral-tibial joint to the other), forespan (from one humero-radio-ulnar joint to the other), gape width (width of head measured at insertion of jaw), head depth (depth of head measured at ear openings), toe pad width (across the widest part of the fourth toe of the right hindlimb), and mass ( Fig. S1). SVL was measured to the nearest millimeter with a metal ruler, mass was measured to the nearest 0.1 g with a Pesola 10.0 g spring scale, and all other measurements were taken with dial calipers to the nearest 0.1 mm. A 2-mm piece of tissue from the distal tip of the tail was collected as a source of DNA for genetic analyses. Each individual was released to its original point of capture after being marked temporarily with a small spot of paint to prevent recapture. We also recorded the diameter of the structure that each individual was perched on (perch diameter) just prior to capture. As A. sagrei are largely sedentary and often use the same perch repeatedly (Calsbeek 2009a), perch diameter at first sighting is a relatively good estimate of habitat use. Lizards that were moving when first sighted were not scored for perch diameter.

Morphological analyses
To assess the effect of habitat type on morphological characters, we used a linear mixed model with two habitats per island (habitat crossed with island) and island as a random effect to account for variation owing to island. Significance values were adjusted using a standard Bonferroni correction. This analysis was performed using the procedure xtmixed in Stata statistical software (StataCorp 2007).
To allow for visualization of size-adjusted morphological traits, we calculated residuals values from the linear regressions of the traits on SVL (all traits are highly correlated with body size in these lizards) (SAS 2006).

Analysis of perch data
Because the presence of flat perches prevented us from treating perch diameter as a continuous variable, we assigned the perch diameters to five categories. Perch diameters were categorized as follows: category 1: 0.1-2.9 cm, category 2: 3.0-5.9 cm, category 3: 6-8.9 cm, category 4: 9-40 cm, and category 5: flat surfaces. We then performed a contingency analysis to test for differences in the frequency distributions of perch diameters between natural and disturbed habitats (SAS 2006).

Analysis of genetic data
We extracted genomic DNA from tail tissue using the QIAGEN QIAamp DNA Mini. Six microsatellite loci were amplified using PCR, and length polymorphism (18-20 individuals per site) was assessed on an automated DNA sequencer (ABI 3700) using fluorescent labeling of one primer in each pair.
We tested for linkage disequilibrium among loci and for Hardy-Weinberg equilibrium using GENEPOP (Raymond and Rousset 1995). We also used GENEPOP to estimate F ST values between study sites. We tested the statistical significance of pairwise F ST values using the G test in the program FSTAT (Goudet 1995). To control for type I error, we applied a sequential Bonferonni correction (Rice 1989). Because relative measures of differentiation, such as estimates of F ST , can be difficult to compare (Hedrick 1999), we also estimated Nei's standard genetic distance (D S ) using FSTAT. We generated a neighborjoining network tree using the allele-sharing distance D AS (Jin and Chakraboty 1993) and the program POPULA-TIONS v.1.2.28 (Langella 2002).

Results
Results from the linear mixed model indicate a significant difference in morphology between lizards occupying disturbed and natural habitats. Of the seven morphological traits, SVL, hindspan, and mass significantly differed between disturbed and natural habitats (Table 1). Individuals in disturbed habitats were larger in SVL, hindspan, and mass than individuals in natural habitats (Fig. 2, Table 1). None of the other morphological characters differed significantly between habitat types after Bonferroni correction (Table 1). Mean and variance values for each of the morphological traits are listed in Table S1.
The percentage of tree cover as determined from satellite remote sensing differed between habitat types, with significantly greater cover in natural sites (t-test, n = 8, F 1,7 = 30.66, P = 0.002). Tree cover in natural sites was much greater than disturbed sites ranging from 56% to 85% and from 2% to 38%, respectively (Fig. 3). These results mirrored observations from fieldwork on the ground that showed relatively few trees in disturbed sites compared with natural sites. Finally, we found a significant negative relationship between percentage tree cover and hind limb length (r = )0.79, P = 0.01).
The frequency distributions of perch diameters significantly differed between habitat types (contingency analysis, n = 160, v 2 = 35.24, P < 0.001). We measured a high frequency of broad perches in disturbed habitat where Model formula for all other traits: Trait = b 0i + b 1 · habitat + b 2 · SVL + e ij b 0i = s 00 + l i . d.f. = 3. *Significant at a Bonferroni-corrected a of 0.007.  percentage canopy cover was lowest. By contrast, narrow perches, but not broad perches, were abundant in natural habitat (Fig. 4). We detected no significant linkage disequilibrium between pairs of microsatellite loci and therefore consider the six loci to be independent markers. Some of the six loci were out of Hardy-Weinberg equilibrium in some of the populations. However, none of the populations were out of equilibrium at all loci, nor was any locus out of equilibrium in all populations. We therefore conclude that our results are not adversely affected by the presence of null alleles.
We detected no significant intra-island population structure between habitats, indicating that the intra-island populations are not derived from separate populations. Pairwise population comparisons of F ST values were not significantly different from zero between pairs of natural and disturbed sites on each of the four islands (Table 2). We did, however, detect genetic divergence between populations on different islands. Pairwise F ST values were significantly different from zero in 18 of the 24 inter-island population comparisons (Table 2). Nei's genetic distance values indicated a similar pattern, with genetic divergence among islands (D S range: 0.312-1.069) far exceeding genetic divergence within islands (D S range: 0.196-0.371) ( Table 2). The eight populations were strongly clustered into four island groups based on neighbor-joining analysis of allele-sharing genetic distances (Fig. 5).

Discussion
Human activities are dramatically altering ecosystems across the planet, but the evolutionary consequences of these changes remain poorly understood (Smith and Bernatchez 2008). Studies of human-driven evolutionary change have documented rapid evolution in the behavior, morphology, and life history of native insects in response to introduced host plants (Carroll and Boyd 1992; Perch diameters are categorized as follows: category 1: perch diameters of 0.1-2.9 cm, category 2: 3.0-5.9 cm, category 3: 6-8.9 cm, category 4: 9-40 cm, and category 5: flat surfaces. Frequency distributions of perch diameters significantly differed between natural and disturbed sites, with a higher frequency of broad perches in disturbed habitat.

Human-induced morphological shifts
Marnocha et al. Singer et al. 1993;Filchak et al. 2000). A number of studies have also shown that increased turbidity because of human activities may disrupt mate recognition in fish and thereby cause extinction through hybridization (Seehausen et al. 1997(Seehausen et al. , 2008Gow et al. 2006;Taylor et al. 2006). A recent surge in work on this topic has led to discoveries of human-induced evolutionary change in a diversity of organisms (Bradshaw and Holzapfel 2001;Parmesan 2006;Slabbekoorn and Ripmeester 2008;Waples et al. 2008), but much remains unknown. Here, we contribute to this growing body of knowledge by examining anthropogenic impacts in Anolis lizards, a model system for studies of adaptive radiation. Our results demonstrate morphological divergence, but not genetic divergence, between populations of A. sagrei inhabiting disturbed and natural habitats on four islands of the Bahamas. Lizards in disturbed habitats were significantly larger in SVL, hindspan, and mass than those in natural habitats. Our results indicate that populations are genetically distinct between each island but are not   genetically differentiated between habitat types on each island. Thus, the observed pattern of intra-island morphological differences cannot be explained by separate founding populations for each habitat nor differential immigration by distinct genetic types. Rather, our results are consistent with morphological differences between habitats having arisen in situ on each island. Lizards in disturbed habitats had longer hindspans than those in natural habitats, likely reflecting differences in habitat. Estimates of percent tree cover from satellite remote sensing across islands show that disturbed habitats have roughly 30-70% less canopy cover than natural sites. Furthermore, hind limb length was negatively correlated with percentage canopy cover, consistent with the assertion that the loss of trees in human-altered habitats leads to morphological divergence. Not only are the disturbed sites more open than the natural sites, they also provide very different perches for the lizards that occupy them. The frequency distributions of perch diameters differed between habitat types, with a higher frequency of broad perches in disturbed habitats. These results are consistent with considerable evidence from previous studies that natural selection shapes the association between limb length and perch diameter in anoles (Williams 1983;Losos 1990;Calsbeek and Irschick 2007). In addition, our results suggest that anthropogenic habitat alteration can have direct impacts on the limb length of these lizards. Many of the broad perch surfaces in the disturbed habitats were novel structures such as walls, boards, and fence posts. The presence of these artificial perches appears to have altered the morphology of lizards inhabiting the disturbed sites.
Differences in habitat structure may also explain the divergence in body size between lizards occupying disturbed and natural habitats. Among ecomorphs, species that perch on large diameter surfaces tend to be larger than those that perch on smaller diameter surfaces (Williams 1972(Williams , 1983. This correlation has also been shown at the intraspecific level . While the heritability of the traits examined here is high (Calsbeek and Bonneaud 2008), suggesting they may evolve by natural selection, the observed divergence in body size may reflect either a plastic response to environmental variation or evolutionary shifts or some combination of both. Future work will focus on common garden breeding experiments from these populations to resolve this question.
Given the differences in tree cover between habitat types, thermal environment likely differs between natural and disturbed environments. This variation in addition to differences in factors such as food availability may explain the observed differences in body size between natural and disturbed study sites. Recent results from a laboratory study that reared lizards under varying food and temperature conditions showed that phenotypic plasticity varied greatly depending on the trait but was not apparent in hindlimb (Marnocha 2008). Hatchlings were placed in treatments that varied in temperature and food availability and were measured every 14 days until they reached adult size (defined as 196 days of age). Male SVL, as well as other indices of body size, differed between treatments in both the growth trajectories and at adult size. In contrast to overall body size, there was no evidence of plasticity in hindlimb length across the varying temperature and food conditions. Although further work is needed, these results suggests that hindlimb differences in the wild are more likely to result from selection pressures than plastic responses to food and temperature. The high heritability of this trait, the lack of evidence of plasticity in the previously described study, as well as the multitude of studies showing limb length shifts in response to natural selection Marnocha 2008;Calsbeek 2009b) suggests that the differences in limb length found here are unlikely to be due solely to environmental effects and may have an evolutionary origin.
Our results demonstrate that anthropogenic habitat alteration may carry evolutionary consequences for species that persist in degraded habitats. We show divergence in several morphological characters between A. sagrei occupying natural and disturbed habitats, despite genetic evidence that the populations occupying the two habitat types are genetically undifferentiated and likely derived from the same founding population. The observed morphological shifts may be a fortuitously adaptive plastic response, the result of an adaptive norm of reaction that has evolved in response to habitat variation, or a genetic response to habitat-specific selection pressures. Future research will be needed to fully address this question. Nevertheless, the influence of human activities on morphology may have important evolutionary implications for Anolis lizards. The repeated evolution of ecomorphs on different islands highlights the importance of ecology in driving diversification in anoles Jackman et al. 1999). Here, we show that human activities, by impacting the ecology, may also be impacting the course of evolution and the patterns of diversification in this group of lizards. comments on previous drafts. We are grateful to the Bahamas Department of Agriculture for permission to conduct research. Funding for this study was provided by grants to EM from the United States Environmental Protection Agency, the American Museum of Natural History, and The Explorer's Club; a National Geographic Society grant and a National Science Foundation grant to TBS.

Literature cited Supporting Information
Additional Supporting Information may be found in the online version of this article: Figure S1. Photographs demonstrating how the following morphological characters were measured: (A) forespan, (B) snout-vent length, (C) gape width, (D) hindspan, (E) head depth, and (F) toe pad width. Table S1. Morphological traits (raw values) in each of the habitat types. Mass is shown in grams; all other measurements are in millimeters. Statistical analyses took SVL into account as a covariate, so that all other trait values were size-adjusted.
Please note: Wiley-Blackwell are not responsible for the content or functionality of any supporting materials supplied by the authors. Any queries (other than missing material) should be directed to the corresponding author for the article.