Genetic differentiation among populations of a Hispaniolan trunk anole that exhibit geographical variation in dewlap colour

The genetic basis of conspicuous coloration in the Guadeloupean anole: Evolution by sexual and ecological selection

Understanding how natural selection acts on the genome and contributes to the process of speciation is a primary aim of the study of evolution. Here we used natural variation in two subspecies of the Guadeloupean anole (Anolis marmoratus ssp.), from the island of Guadeloupe in the Lesser Antilles, to explore the genomic basis of adaptation and speciation in Anolis lizards. These subspecies inhabit distinct ecological environments and display marked differences in adult male color and pattern. We sequenced the complete genomes of 20 anoles, 10 from each subspecies, at 1.4× coverage. We used genome-wide scans of population differentiation, allele frequency spectrum, and linkage disequilibrium to characterize the genomic architecture within and between the subspecies. While most of the genome was undifferentiated, we observed five large divergent regions. Within these regions we identified blocks, 5 kb pairs in length, enriched for fixed single nucleotide polymorphisms. These blocks encompass 97 genes, two of which are candidate pigmentation genes. One is melanophilin (mlph), which helps transport melanosomes within melanocytes. The other is a cluster of differentiation 36 (cd36), which regulates carotenoid pigment sequestration. We used high-pressure liquid chromatography to confirm that carotenoid pigments are significantly more abundant in the conspicuous orange-pigmented skin of male A. m. marmoratus suggesting that cd36 may be regulating pigment deposition in this tissue. We identified for the first time a carotenoid gene that is a potential target of divergent sexual selection and may be contributing to the early stages of speciation in Anolis lizards.

Dewlap colour variation in Anolis sagrei is maintained among habitats within islands of the West Indies

Animal signals evolve in an ecological context. Locally adapting animal sexual signals can be especially important for initiating or reinforcing reproductive isolation during the early stages of speciation. Previous studies have demonstrated that dewlap colour in Anolis lizards can be highly variable between populations in relation to both biotic and abiotic adaptive drivers at relatively large geographical scales. Here, we investigated differentiation of dewlap coloration among habitat types at a small spatial scale, within multiple islands of the West Indies, to test the hypothesis that similar local adaptive processes occur over smaller spatial scales. We explored variation in dewlap coloration in the most widespread species of anole, Anolis sagrei, across three characteristic habitats spanning the Bahamas and the Cayman Islands, namely beach scrub, primary coppice forest and mangrove forest. Using reflectance spectrometry paired with supervised machine learning, we found significant differences in spectral properties of the dewlap between habitats within small islands, sometimes over very short distances. Passive divergence in dewlap phenotype associated with isolation-by-distance did not seem to explain our results. On the other hand, these habitat-specific dewlap differences varied in magnitude and direction across islands, and thus, our primary test for adaptation-parallel responses across islands-was not supported. We suggest that neutral processes or selection could be involved in several ways, including sexual selection. Our results shed new light on the scale at which signal colour polymorphism can be maintained in the presence of gene flow, and the relative role of local adaptation and other processes in driving these patterns of dewlap colour variation across islands.

Ecological Opportunity from Innovation, not Islands, Drove the Anole Lizard Adaptive Radiation

Islands are thought to facilitate adaptive radiation by providing release from competition and predation. Anole lizards are considered a classic example of this phenomenon: different ecological specialists ('ecomorphs') evolved in the Caribbean Greater Antilles (Cuba, Hispaniola, Jamaica, and Puerto Rico), resulting in convergent assemblages that are not observed in mainland Latin America. Yet, the role of islands in facilitating adaptive radiation is more often implied than directly tested, leaving uncertain the role of biogeography in stimulating diversification. Here, we assess the proposed "island effect" on anole diversification using Bayesian phylogenetic comparative methods that explicitly incorporate rate heterogeneity across the tree and demonstrate two cases of would-be false positives. We discovered that rates of speciation and morphological evolution of island and mainland anoles are equivalent, implying that islands provide no special context for exceptionally rapid diversification. Likewise, rates of evolution were equivalent between island anoles that arose via in situ versus dispersal-based mechanisms, and we found no evidence for island-specific rates of speciation or morphological evolution. Nonetheless, the origin of Anolis is characterized by a speciation pulse that slowed over time - a classic signature of waning ecological opportunity. Our findings cast doubt on the notion that islands catalyzed the anole adaptive radiation and instead point to a key innovation, adhesive toe pads, which facilitated the exploitation of many arboreal niches sparsely utilized by other iguanian lizards. The selective pressures responsible for arboreal niche diversification differ between islands and the mainland, but the tempo of diversification driven by these discordant processes is indistinguishable.

Evolutionary drivers of sexual signal variation in Amazon Slender Anoles

Phenotypic variation among populations, as seen in the signaling traits of many species, provides an opportunity to test whether similar factors generate repeated phenotypic patterns in different parts of a species' range. We investigated whether genetic divergence, abiotic gradients, and sympatry with closely related species explain variation in the dewlap colors of Amazon Slender Anoles, Anolis fuscoauratus. To this aim, we characterized dewlap diversity in the field with respect to population genetic structure and evolutionary relationships, assessed whether dewlap phenotypes are associated with climate or landscape variables, and tested for nonrandom associations in the distributions of A. fuscoauratus phenotypes and sympatric Anolis species. We found that dewlap colors vary among but not within sites in A. fuscoauratus. Regional genetic clusters included multiple phenotypes, while populations with similar dewlaps were often distantly related. Phenotypes did not segregate in environmental space, providing no support for optimized signal transmission at a local scale. Instead, we found a negative association between certain phenotypes and sympatric Anolis species with similar dewlap color attributes, suggesting that interactions with closely related species promoted dewlap divergence among A. fuscoauratus populations. Amazon Slender Anoles emerge as a promising system to address questions about parallel trait evolution and the contribution of signaling traits to speciation.

A morphometric assessment of species boundaries in a widespread anole lizard (Squamata: Dactyloidae)

Cryptic species – genetically distinct species that are morphologically difficult to distinguish – present challenges to systematists. Operationally, cryptic species are very difficult to identify and sole use of genetic data or morphological data can fail to recognize evolutionarily isolated lineages. We use morphometric data to test species boundaries hypothesized with genetic data in the North Caribbean bark anole (Anolis distichus), a suspected species complex. We use univariate and multivariate analyses to test if candidate species based on genetic data can be accurately diagnosed. We also test alternative species delimitation scenarios with a model fitting approach that evaluates normal mixture models capable of identifying morphological clusters. Our analyses reject the hypothesis that the candidate species are diagnosable. Neither uni- nor multivariate morphometric data distinguish candidate species. The best-supported model included two morphological clusters; however, these clusters were uneven and did not align with a plausible species divergence scenario. After removing two related traits driving this result, only one cluster was supported. Despite substantial differentiation revealed by genetic data, we recover no new evidence to delimit species and refrain from taxonomic revision. This study highlights the importance of considering other types of data along with molecular data when delimiting species.

Does breeding season variation affect evolution of a sexual signaling trait in a tropical lizard clade?

Sexually selected traits can be expected to increase in importance when the period of sexual behavior is constrained, such as in seasonally restricted breeders. Anolis lizard male dewlaps are classic examples of multifaceted signaling traits, with demonstrated intraspecific reproductive function reflected in courtship behavior. Fitch and Hillis found a correlation between dewlap size and seasonality in mainland Anolis using traditional statistical methods and suggested that seasonally restricted breeding seasons enhanced the differentiation of this signaling trait. Here, we present two tests of the Fitch-Hillis Hypothesis using new phylogenetic and morphological data sets for 44 species of Mexican Anolis. A significant relationship between dewlap size and seasonality is evident in phylogenetically uncorrected analyses but erodes once phylogeny is accounted for. This loss of strong statistical support for a relationship between a key aspect of dewlap morphology and seasonality also occurs within a species complex (A. sericeus group) that inhabits seasonal and aseasonal environments. Our results fail to support seasonality as a strong driver of evolution of Anolis dewlap size. We discuss the implications of our results and the difficulty of disentangling the strength of single mechanisms on trait evolution when multiple selection pressures are likely at play.

Phylogeographic and phenotypic outcomes of brown anole colonization across the Caribbean provide insight into the beginning stages of an adaptive radiation

Some of the most important insights into the ecological and evolutionary processes of diversification and speciation have come from studies of island adaptive radiations, yet relatively little research has examined how these radiations initiate. We suggest that Anolis sagrei is a candidate for understanding the origins of the Caribbean Anolis adaptive radiation and how a colonizing anole species begins to undergo allopatric diversification, phenotypic divergence and, potentially, speciation. We undertook a genomic and morphological analysis of representative populations across the entire native range of A. sagrei, finding that the species originated in the early Pliocene, with the deepest divergence occurring between western and eastern Cuba. Lineages from these two regions subsequently colonized the northern Caribbean. We find that at the broadest scale, populations colonizing areas with fewer closely related competitors tend to evolve larger body size and more lamellae on their toepads. This trend follows expectations for post-colonization divergence from progenitors and convergence in allopatry, whereby populations freed from competition with close relatives evolve towards common morphological and ecological optima. Taken together, our results show a complex history of ancient and recent Cuban diaspora with populations on competitor-poor islands evolving away from their ancestral Cuban populations regardless of their phylogenetic relationships, thus providing insight into the original diversification of colonist anoles at the beginning of the radiation. Our research also supplies an evolutionary framework for the many studies of this increasingly important species in ecological and evolutionary research.

Biogeography, Systematics, and Ecomorphology of Pacific Island Anoles

: Anoles are regarded as important models for understanding dynamic processes in ecology and evolution. Most work on this group has focused on species in the Caribbean Sea, and recently in mainland South and Central America. However, the Eastern Tropical Pacific (ETP) is home to seven species of anoles from three unique islands (Islas Cocos, Gorgona, and Malpelo) that have been largely overlooked. Four of these species are endemic to single islands (Norops townsendi on Isla Cocos, Dactyloa agassizi on Isla Malpelo, D. gorgonae and N. medemi on Isla Gorgona). Herein, we present a phylogenetic analysis of anoles from these islands in light of the greater anole phylogeny to estimate the timing of divergence from mainland lineages for each species. We find that two species of solitary anoles (D. agassizi and N. townsendi) diverged from mainland ancestors prior to the emergence of their respective islands. We also present population-wide morphological data suggesting that both display sexual size dimorphism, similar to single-island endemics in the Caribbean. All lineages on Isla Gorgona likely arose during past connections with South America, and ecologically partition their habitat. Finally, we highlight the importance of conservation of these species and island fauna in general.

Cenozoic aridization in Central Eurasia shaped diversification of toad-headed agamas (Phrynocephalus; Agamidae, Reptilia)

We hypothesize the phylogenetic relationships of the agamid genus Phrynocephalus to assess how past environmental changes shaped the evolutionary and biogeographic history of these lizards and especially the impact of paleogeography and climatic factors. Phrynocephalus is one of the most diverse and taxonomically confusing lizard genera. As a key element of Palearctic deserts, it serves as a promising model for studies of historical biogeography and formation of arid habitats in Eurasia. We used 51 samples representing 33 of 40 recognized species of Phrynocephalus covering all major areas of the genus. Molecular data included four mtDNA (COI, ND2, ND4, Cytb; 2,703 bp) and four nuDNA protein-coding genes (RAG1, BDNF, AKAP9, NKTR; 4,188 bp). AU-tests were implemented to test for significant differences between mtDNA- and nuDNA-based topologies. A time-calibrated phylogeny was estimated using a Bayesian relaxed molecular clock with nine fossil calibrations. We reconstructed the ancestral area of origin, biogeographic scenarios, body size, and the evolution of habitat preference. Phylogenetic analyses of nuDNA genes recovered a well-resolved and supported topology. Analyses detected significant discordance with the less-supported mtDNA genealogy. The position of Phrynocephalus mystaceus conflicted greatly between the two datasets. MtDNA introgression due to ancient hybridization best explained this result. Monophyletic Phrynocephalus contained three main clades: (I) oviparous species from south-western and Middle Asia; (II) viviparous species of Qinghai–Tibetan Plateau (QTP); and (III) oviparous species of the Caspian Basin, Middle and Central Asia. Phrynocephalus originated in late Oligocene (26.9 Ma) and modern species diversified during the middle Miocene (14.8–13.5 Ma). The reconstruction of ancestral areas indicated that Phrynocephalus originated in Middle East–southern Middle Asia. Body size miniaturization likely occurred early in the history of Phrynocephalus. The common ancestor of Phrynocephalus probably preferred sandy substrates with the inclusion of clay or gravel. The time of Agaminae radiation and origin of Phrynocephalus in the late Oligocene significantly precedes the landbridge between Afro-Arabia and Eurasia in the Early Miocene. Diversification of Phrynocephalus coincides well with the mid-Miocene climatic transition when a rapid cooling of climate drove progressing aridification and the Paratethys salinity crisis. These factors likely triggered the spreading of desert habitats in Central Eurasia, which Phrynocephalus occupied. The origin of the viviparous Tibetan clade has been associated traditionally with uplifting of the QTP; however, further studies are needed to confirm this. Progressing late Miocene aridification, the decrease of the Paratethys Basin, orogenesis, and Plio–Pleistocene climate oscillations likely promoted further diversification within Phrynocephalus. We discuss Phrynocephalus taxonomy in scope of the new analyses.

Climatic niche evolution is faster in sympatric than allopatric lineages of the butterfly genus Pyrgus

Understanding how speciation relates to ecological divergence has long fascinated biologists. It is assumed that ecological divergence is essential to sympatric speciation, as a mechanism to avoid competition and eventually lead to reproductive isolation, while divergence in allopatry is not necessarily associated with niche differentiation. The impact of the spatial context of divergence on the evolutionary rates of abiotic dimensions of the ecological niche has rarely been explored for an entire clade. Here, we compare the magnitude of climatic niche shifts between sympatric versus allopatric divergence of lineages in butterflies. By combining next-generation sequencing, parametric biogeography and ecological niche analyses applied to a genus-wide phylogeny of Palaearctic Pyrgus butterflies, we compare evolutionary rates along eight climatic dimensions across sister lineages that diverged in large-scale sympatry versus allopatry. In order to examine the possible effects of the spatial scale at which sympatry is defined, we considered three sets of biogeographic assignments, ranging from narrow to broad definition. Our findings suggest higher rates of niche evolution along all climatic dimensions for sister lineages that diverge in sympatry, when using a narrow delineation of biogeographic areas. This result contrasts with significantly lower rates of climatic niche evolution found in cases of allopatric speciation, despite the biogeographic regions defined here being characterized by significantly different climates. Higher rates in allopatry are retrieved when biogeographic areas are too widely defined-in such a case allopatric events may be recorded as sympatric. Our results reveal the macro-evolutionary significance of abiotic niche differentiation involved in speciation processes within biogeographic regions, and illustrate the importance of the spatial scale chosen to define areas when applying parametric biogeographic analyses.

Condition dependence of shared traits differs between sympatric Anolis lizards

In many species, sexually selected combat and display traits can confer dramatic fitness benefits to males by aiding in mate acquisition, so individuals maximally invest energy into their growth and maintenance. Such traits are deemed condition-dependent, as the energy that is available for investment depends on the health and condition of the individual. Condition dependence is present in a wide range of traits across many taxa, but the extent to which condition dependence varies among shared traits in closely related species is poorly understood. We tested for condition dependence in two Anolis lizards, the ground anole (Anolis humilis) and the slender anole (Anolis limifrons). Specifically, we measured dewlap area, jaw length, jaw width, and sprint speed, and then regressed one of two indices of body condition on each of these traits for both species. We found that dewlap area and jaw width exhibited condition dependence in ground anoles, but not in slender anoles. Sprint speed and jaw length were not condition-dependent in either species. The presence of condition dependence in ground anoles, but not slender anoles, implies evolutionary liability in the condition dependence of shared traits. Additionally, the fact that condition dependence was only detected in ground anoles, which have a greater relative dewlap size, suggests a potential role for signal evolution or strength of sexual selection in the evolution of condition dependence. Finally, our research suggests that variation in condition dependence of sexually selected traits could have implications for the evolution of dewlap diversity among Anolis species.

Molecular ecology studies of species radiations: current research gaps, opportunities and challenges

Understanding the drivers and limits of species radiations is a crucial goal of evolutionary genetics and molecular ecology, yet research on this topic has been hampered by the notorious difficulty of connecting micro- and macroevolutionary approaches to studying the drivers of diversification. To chart the current research gaps, opportunities and challenges of molecular ecology approaches to studying radiations, we examine the literature in the journal Molecular Ecology and revisit recent high-profile examples of evolutionary genomic research on radiations. We find that available studies of radiations are highly unevenly distributed among taxa, with many ecologically important and species-rich organismal groups remaining severely understudied, including arthropods, plants and fungi. Most studies employed molecular methods suitable over either short or long evolutionary time scales, such as microsatellites or restriction site-associated DNA sequencing (RAD-seq) in the former case and conventional amplicon sequencing of organellar DNA in the latter. The potential of molecular ecology studies to address and resolve patterns and processes around the species level in radiating groups of taxa is currently limited primarily by sample size and a dearth of information on radiating nuclear genomes as opposed to organellar ones. Based on our literature survey and personal experience, we suggest possible ways forward in the coming years. We touch on the potential and current limitations of whole-genome sequencing (WGS) in studies of radiations. We suggest that WGS and targeted ('capture') resequencing emerge as the methods of choice for scaling up the sampling of populations, species and genomes, including currently understudied organismal groups and the genes or regulatory elements expected to matter most to species radiations.

A genomic assessment of species boundaries and hybridization in a group of highly polymorphic anoles (distichus species complex)

Delimiting young species is one of the great challenges of systematic biology, particularly when the species in question exhibit little morphological divergence. Anolis distichus, a trunk anole with more than a dozen subspecies that are defined primarily by dewlap color, may actually represent several independent evolutionary lineages. To test this, we utilized amplified fragment length polymorphisms (AFLP) genome scans and genetic clustering analyses in conjunction with a coalescent‐based species delimitation method. We examined a geographically widespread set of samples and two heavily sampled hybrid zones. We find that genetic divergence is associated with a major biogeographic barrier, the Hispaniolan paleo‐island boundary, but not with dewlap color. Additionally, we find support for hypotheses regarding colonization of two Hispaniolan satellite islands and the Bahamas from mainland Hispaniola. Our results show that A. distichus is composed of seven distinct evolutionary lineages still experiencing a limited degree of gene flow. We suggest that A. distichus merits taxonomic revision, but that dewlap color cannot be relied upon as the primary diagnostic character.

Diversification in wild populations of the model organism Anolis carolinensis: A genome-wide phylogeographic investigation

The green anole (Anolis carolinensis) is a lizard widespread throughout the southeastern United States and is a model organism for the study of reproductive behavior, physiology, neural biology, and genomics. Previous phylogeographic studies of A. carolinensis using mitochondrial DNA and small numbers of nuclear loci identified conflicting and poorly supported relationships among geographically structured clades; these inconsistencies preclude confident use of A. carolinensis evolutionary history in association with morphological, physiological, or reproductive biology studies among sampling localities and necessitate increased effort to resolve evolutionary relationships among natural populations. Here, we used anchored hybrid enrichment of hundreds of genetic markers across the genome of A. carolinensis and identified five strongly supported phylogeographic groups. Using multiple analyses, we produced a fully resolved species tree, investigated relative support for each lineage across all gene trees, and identified mito‐nuclear discordance when comparing our results to previous studies. We found fixed differences in only one clade—southern Florida restricted to the Everglades region—while most polymorphisms were shared between lineages. The southern Florida group likely diverged from other populations during the Pliocene, with all other diversification during the Pleistocene. Multiple lines of support, including phylogenetic relationships, a latitudinal gradient in genetic diversity, and relatively more stable long‐term population sizes in southern phylogeographic groups, indicate that diversification in A. carolinensis occurred northward from southern Florida.

Adaptive signal coloration maintained in the face of gene flow in a Hispaniolan Anolis Lizard

Background

Studies of geographic variation can provide insight into the evolutionary processes involved in the early stages of biological diversification. In particular, multiple, replicated cases of geographic trait divergence present a powerful approach to study how patterns of introgression and adaptive divergence can vary with geographic space and time. In this study, we conduct replicated, fine-scaled molecular genetic analyses of striking geographic dewlap color variation of a Hispaniolan Anolis lizard, Anolis distichus, to investigate whether adaptive trait divergence is consistently associated with speciation, whereby genetic divergence is observed with neutral markers, or whether locally adapted traits are maintained in the face of continued gene flow.

Results

We find instances where shifts in adaptive dewlap coloration across short geographic distances are associated with reproductive isolation as well as maintained in the face of gene flow, suggesting the importance of both processes in maintaining geographic dewlap variation.

Conclusion

Our study suggests that adaptive dewlap color differences are maintained under strong divergent natural selection, but this divergence does not necessarily lead to anole speciation.

Electronic supplementary material

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

Ecological distributions, phenological isolation, and genetic structure in sympatric and parapatric populations of the Larrea tridentata polyploid complex

PREMISE OF THE STUDY

Polyploidy is widely recognized as a mechanism of diversification. Contributions of polyploidy to specific pre- and postzygotic barriers-and classifications of polyploid speciation as "ecological" vs. "non-ecological"-are more contentious. Evaluation of these issues requires comprehensive studies that test ecological characteristics of cytotypes as well as the coincidence of genetic structure with cytotype distributions.

METHODS

We investigated a classical example of autopolyploid speciation, Larrea tridentata, at multiple areas of cytotype co-occurrence. Habitat and phenological differences were compared between diploid, tetraploid, and hexaploid populations on the basis of edaphic, community composition, and flowering time surveys. Frequency of hybridization between diploids and tetraploids was investigated using a diploid-specific chloroplast DNA (cpDNA) marker; genetic structure for all cytotypes was assessed using amplified fragment length polymorphisms (AFLPs).

KEY RESULTS

Across contact zones, we found cytotypes in habitats distinguished by soil and vegetation. We observed modest differences in timing and production of flowers, indicating a degree of assortative mating that was asymmetric between cytotypes. Nonetheless, cpDNA analyses in diploid-tetraploid contact zones suggested that ∼5% of tetraploid plants had hybrid origins involving unilateral sexual polyploidization. Genetic structure of AFLPs largely coincided with cytotype distributions in diploid-tetraploid contact zones. In contrast, there was little structure in areas of contact between tetraploids and hexaploids, suggesting intercytotype gene flow or recurrent hexaploid formation.

CONCLUSIONS

Diploid, tetraploid, and hexaploid cytotypes of L. tridentata are segregated by environmental distributions and flowering phenology in contact zones, with diploid and tetraploid populations having corresponding differences in genetic structure.

Biogeography and evolution of a widespread Central American lizard species complex: Norops humilis, (Squamata: Dactyloidae)

Background

Caribbean anole lizards (Dactyloidae) have frequently been used as models to study questions regarding biogeography and adaptive radiations, but the evolutionary history of Central American anoles (particularly those of the genus Norops) has not been well studied. Previous work has hypothesized a north-to-south dispersal pattern of Central American Norops, but no studies have examined dispersal within any Norops lineages. Here we test two major hypotheses for the dispersal of the N. humilis/quaggulus complex (defined herein, forming a subset within Savage and Guyer’s N. humilis group).

Results

Specimens of the N. humilis group were collected in Central America, from eastern Mexico to the Canal Zone of Panama. Major nodes were dated for comparison to the geologic history of Central America, and ancestral ranges were estimated for the N. humilis/quaggulus complex to test hypothesized dispersal patterns. These lineages displayed a northward dispersal pattern. We also demonstrate that the N. humilis/quaggulus complex consists of a series of highly differentiated mitochondrial lineages, with more conserved nuclear evolution. The paraphyly of the N. humilis species group is confirmed. A spatial analysis of molecular variance suggests that current populations are genetically distinct from one another, with limited mitochondrial gene flow occurring among sites.

Conclusions

The observed south-to-north colonization route within the Norops humilis/quaggulus complex represents the first evidence of a Norops lineage colonizing in a south-to-north pattern, (opposite to the previously held hypothesis for mainland Norops). One previously described taxon (N. quaggulus) was nested within N. humilis, demonstrating the paraphyly of this species; while our analyses also reject the monophyly of the Norops humilis species group (sensu Savage and Guyer), with N. tropidonotus, N. uniformis, and N. marsupialis being distantly related to/highly divergent from the N. humilis/quaggulus complex. Our work sheds light on mainland anole biogeography and past dispersal events, providing a pattern to test against other groups of mainland anoles.

Electronic supplementary material

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

Multilocus phylogenetic analyses of Hispaniolan and Bahamian trunk anoles (distichus species group)

The distichus species group includes six species and 21 subspecies of trunk ecomorph anoles distributed across Hispaniola and its satellite islands as well as the northern Bahamas. Although this group has long served as a model system for studies of reproductive character displacement, adaptation, behavior and speciation, it has never been the subject of a comprehensive phylogenetic analysis. Our goal here is to generate a multilocus phylogenetic dataset (one mitochondrial and seven nuclear loci) and to use this dataset to infer phylogenetic relationships among the majority of the taxa assigned to the distichus species group. We use these phylogenetic trees to address three topics about the group's evolution. First, we consider longstanding taxonomic controversies about the status of several species and subspecies assigned to the distichus species group. Second, we investigate the biogeographic history of the group and specifically test the hypotheses that historical division of Hispaniola into two paleo-islands contributed to the group's diversification and that Bahamian and Hispaniolan satellite island populations are derived from colonists from the main Hispaniolan landmass. Finally, third, we use comparative phylogenetic analyses to test the hypothesis that divergence between pale yellow and darkly pigmented orange or red dewlap coloration has occurred repeatedly across the distichus species group.

Phylogenetic relationships of Amazonian anole lizards (Dactyloa): taxonomic implications, new insights about phenotypic evolution and the timing of diversification

The ecology and evolution of Caribbean anoles are well described, yet little is known about mainland anole species. Lack of phylogenetic information limits our knowledge about species boundaries, morphological evolution, and the biogeography of anoles in South America. To help fill this gap, we provide an updated molecular phylogeny of the Dactyloa (Dactyloidae), with emphasis on the punctata species group. By sampling understudied Amazonian taxa, we (i) assess the phylogenetic placement of the 'odd anole', D. dissimilis; (ii) infer the relationships of the proboscis-bearing D. phyllorhina, testing the hypothesis of independent nasal appendage evolution within the anole radiation; and (iii) examine genetic and dewlap color variation in D. punctata and D. philopunctata. Combining multiple nuclear loci with a review of the fossil record, we also (iv) estimate divergence times within the pleurodont iguanian clade of lizards, including Amazonian representatives of Dactyloa and Norops (Dactyloidae) and of Polychrus (Polychrotidae). We recover the five Dactyloa clades previously referred to as the aequatorialis, heteroderma, latifrons, punctata and roquet species groups, as well as a sixth clade composed of D. dissimilis and the non-Amazonian D. neblinina and D. calimae. We find D. phyllorhina to be nested within the punctata group, suggesting independent evolution of the anole proboscis. We consistently recover D. philopunctata nested within D. punctata, and report limited genetic divergence between distinct dewlap phenotypes. The most recent common ancestor of Dactyloa, Anolis and Norops dates back to the Eocene. Most Amazonian taxa within both Dactyloa and Norops diverged in the Miocene, but some diversification events were as old as the late Eocene and late Oligocene. Amazonian Polychrus diverged in the Pliocene. Our findings have broad implications for anole biogeography, disputing recent suggestions that modern dactyloid genera were present in the Caribbean region during the Cretaceous.

Inferred vs realized patterns of gene flow: an analysis of population structure in the Andros Island Rock Iguana

Ecological data, the primary source of information on patterns and rates of migration, can be integrated with genetic data to more accurately describe the realized connectivity between geographically isolated demes. In this paper we implement this approach and discuss its implications for managing populations of the endangered Andros Island Rock Iguana, Cyclura cychlura cychlura. This iguana is endemic to Andros, a highly fragmented landmass of large islands and smaller cays. Field observations suggest that geographically isolated demes were panmictic due to high, inferred rates of gene flow. We expand on these observations using 16 polymorphic microsatellites to investigate the genetic structure and rates of gene flow from 188 Andros Iguanas collected across 23 island sites. Bayesian clustering of specimens assigned individuals to three distinct genotypic clusters. An analysis of molecular variance (AMOVA) indicates that allele frequency differences are responsible for a significant portion of the genetic variance across the three defined clusters (F_st =  0.117, p0.01). These clusters are associated with larger islands and satellite cays isolated by broad water channels with strong currents. These findings imply that broad water channels present greater obstacles to gene flow than was inferred from field observation alone. Additionally, rates of gene flow were indirectly estimated using BAYESASS 3.0. The proportion of individuals originating from within each identified cluster varied from 94.5 to 98.7%, providing further support for local isolation. Our assessment reveals a major disparity between inferred and realized gene flow. We discuss our results in a conservation perspective for species inhabiting highly fragmented landscapes.

Past hybridization between two East Asian long-tailed tits (Aegithalos bonvaloti and A. fuliginosus)

Introduction

Incomplete lineage sorting and hybridization are two major nonexclusive causes of haplotype sharing between species. Distinguishing between these two processes is notoriously difficult as they can generate similar genetic signatures. Previous studies revealed that the mitochondrial DNA (mtDNA) differentiation between two East Asian long-tailed tits (Aegithalos bonvaloti and A. fuliginosus) was extremely low, even lower than intraspecific differentiation in some other long-tailed tits. Using a combination of multilocus and coalescent analyses, we explored the causes of the anomalous lack of mtDNA differentiation between the two species.

Results

The mtDNA divergence between the two species was shallow, while the nuclear DNA (nuDNA) divergence was considerably deeper. The IMa analyses based on the mtDNA dataset suggested relatively high gene flow from A. fuliginosus to A. bonvaloti, while negligible gene flow in the opposite direction. In contrast to mtDNA, the migration rates at autosomal and Z-linked nuDNA loci were negligible or much lower. The NEWHYBRIDS analysis assigned all individuals except one to pure parental species with high posterior probability. The Bayesian skyline plot showed that both species underwent population expansions during the Last Glacial Maximum (LGM), and the ecological niche modelling suggested that their ranges overlapped more during the LGM than at present.

Conclusions

We suggest that historical hybridization, in combination with selective sweep and/or genetic drift might be the main causes of the extremely low mtDNA differentiation between the two species. The hybridization probably occurred mainly between A. fuliginosus females and A. bonvaloti males. The LGM distribution expansion might have facilitated hybridization, while the post-LGM distribution contraction could have facilitated some mtDNA sorting. Ongoing hybridization between the two species might be very limited, but further studies with more samples from the contact zone are needed to test this conclusion.

The role of heritable and dietary factors in the sexual signal of a Hispaniolan Anolis lizard, Anolis distichus

The diversity of sexual signals is astounding, and divergence in these traits is believed to be associated with the early stages of speciation. An increasing number of studies also suggest a role for natural selection in driving signal divergence for effective transmission in heterogeneous environments. Both speciation and adaptive divergence, however, are contingent on the sexual signal being heritable, yet this often remains assumed and untested. It is particularly critical that the heritability of carotenoid-based sexual signals is investigated because such traits may instead be phenotypically plastic indicators of an individual's quality that exhibit no or little heritable variation. We present the first study to investigate the relative contribution of genetic and environmental factors to the striking diversity of dewlap color and pattern in Anolis lizards. Using a breeding experiment with Anolis distichus populations exhibiting different dewlap phenotypes, we raise F1 offspring in a common garden experiment to assess whether dewlap color is inherited. We follow this with carotenoid supplementation to investigate the influence of dietary pigments to dewlap color variation. We find significant differences in several aspects of dewlap color and pattern to persist to the F1 generation (fathers: N = 19; F1 males: N = 50; P < 0.01) with no change in dewlap phenotype with carotenoid supplementation (N = 52; P > 0.05). These results strongly support that genetic differences underlie dewlap color variation, thereby satisfying a key requirement of natural selection. Our findings provide an important stepping-stone to understanding the evolution of an incredibly diverse signal important for sexual selection and species recognition.

Correlation between Anolis lizard dewlap phenotype and environmental variation indicates adaptive divergence of a signal important to sexual selection and species recognition

Although the importance of signals involved in species recognition and sexual selection to speciation is widely recognized, the processes that underlie signal divergence are still a matter of debate. Several possible processes have been hypothesized, including genetic drift, arbitrary sexual selection, and adaptation to local signaling environments. We use comparative analyses to investigate whether the remarkable geographic variation of dewlap phenotype in a Hispaniolan trunk Anolis lizard (A. distichus) is a result of adaptive signal divergence to heterogeneous environments. We recover a repeated pattern of divergence in A. distichus dewlap color, pattern, and size with environmental variation across Hispaniola. These results are aligned with ecological models of signal divergence and provide strong evidence for dewlap adaptation to local signaling environments. We also find that A. distichus dewlaps vary with the environment in a different manner to other previously studied anoles, thus expanding upon previous predictions on the direction dewlaps will diverge in perceptual color space in response to the environment.

Synergy between Allopatry and Ecology in Population Differentiation and Speciation

The general diversity pattern of the Caribbean anole radiation has been described in detail; however, the actual mechanisms at the origin of their diversification remain controversial. In particular, the role of ecological speciation, and the relative importance of divergence in allopatry and in parapatry, is debated. We describe the genetic structure of anole populations across lineage contact zones and ecotones to investigate the effect of allopatric divergence, natural selection, and the combination of both factors on population differentiation. Allopatric divergence had no significant impact on differentiation across the lineage boundary, while a clear bimodality in genetic and morphological characters was observed across an ecotone within a single lineage. Critically, the strongest differentiation was observed when allopatry and ecology act together, leading to a sharp reduction in gene flow between two lineages inhabiting different habitats. We suggest that, for Caribbean anoles to reach full speciation, a synergistic combination of several historical and ecological factors may be requisite.