Using genomic data to revisit an early example of reproductive character displacement in HaitianAnolislizards

Genomic analysis of the rhesus macaque (Macaca mulatta) and the cynomolgus macaque (Macaca fascicularis) uncover polygenic signatures of reinforcement speciation

Speciation can involve phases of divergent adaptation in allopatry and ecological/reproductive character displacement in sympatry or parapatry. Reproductive character displacement can result as a means of preventing hybridization, a process known as reinforcement speciation. In this study, we use whole-genome sequencing (WGS) of two closely related primate species that have experienced introgression in their history, the rhesus (Macaca mulatta) and cynomolgus (M. fascicularis) macaques, to identify genes exhibiting reproductive character displacement and other patterns consistent with reinforcement speciation. Using windowed scans of various population genetic statistics to identify signatures of reinforcement, we find 184 candidate genes associated with a variety of functions, including an overrepresentation of multiple neurological functions and several genes involved in sexual development and gametogenesis. These results are consistent with a variety of genes acting in a reinforcement process between these species. We also find signatures of introgression of the Y-chromosome that confirm previous studies suggesting male-driven introgression of M. mulatta into M. fascicularis populations. This study uses WGS to find evidence of the process of reinforcement in primates that have medical and conservation relevance.

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.

Experimental test of selection against hybridization as a driver of avian signal divergence

Signal divergence may be pivotal in the generation and maintenance of new biodiversity by allowing closely related species to avoid some costs of co-occurrence. In birds, closely related, sympatric species are more divergent in their colour patterns than those that live apart, but the selective pressures driving this pattern remain unclear. Traditionally, signal divergence among sympatric species is thought to result from selection against hybridization, but broad evidence is lacking. Here, we conducted field experiments on naïve birds using spectrometer-matched, painted 3D-printed models to test whether selection against hybridization drives colour pattern divergence in the genus Poecile. To address selection for male colour pattern divergence without the influence of learning or the evolution of female discrimination in sympatry, we simulated secondary contact between Poecile species, and conducted mate choice experiments on naïve, allopatric females. We found that female black-capped chickadees (Poecile atricapillus) are equally likely to perform copulation solicitation displays to sympatric and allopatric heterospecific congeners when they are paired with conspecifics, but exhibit a strong preference for less divergent males when presented with paired heterospecific congeners. These results suggest that increased colour pattern divergence among sympatric species can reduce the likelihood of mixed mating in some contexts, and therefore should be favoured by selection against hybridization.

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.

Experimental tests of selection against heterospecific aggression as a driver of avian colour pattern divergence

Signal divergence is thought to reduce the costs of co-occurrence for closely related species and may thereby be important in the generation and maintenance of new biodiversity. In birds, closely related, sympatric species are more divergent in their colour patterns than those that live apart, but the selective pressures driving sympatric divergence in colour pattern are not well-understood. Here, we conducted field experiments on naïve birds using spectrometer-matched, painted, 3D-printed models to test whether selection against heterospecific aggression might drive colour pattern divergence in the genus Poecile. We found that territorial male black-capped chickadees (P. atricapillus) are equally likely to attack sympatric and allopatric congeners, and wintering flocks are equally likely to visit feeders occupied by sympatric and allopatric congeners, despite sympatric congeners being more divergent in colour pattern. These results suggest that either the concerted evolution of additional traits (e.g. discrimination), or interactions in sympatry that promote learning, is required if colour pattern divergence among sympatric species is to reduce heterospecific aggression. Alternatively, colour pattern divergence among sympatric species may be caused by other selective pressures, such as selection against hybridization or habitat partitioning and secondary signal adaptation.

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.

The evolution of polymorphism in the warning coloration of the Amazonian poison frog Adelphobates galactonotus

While intraspecific variation in aposematic signals can be selected for by different predatory responses, their evolution is also contingent on other processes shaping genetic variation. We evaluate the relative contributions of selection, geographic isolation, and random genetic drift to the evolution of aposematic color polymorphism in the poison frog Adelphobates galactonotus, distributed throughout eastern Brazilian Amazonia. Dorsal coloration was measured for 111 individuals and genetic data were obtained from 220 individuals at two mitochondrial genes (mtDNA) and 7963 Single Nucleotide Polymorphisms (SNPs). Four color categories were described (brown, blue, yellow, orange) and our models of frog and bird visual systems indicated that each color was distinguishable for these taxa. Using outlier and correlative analyses we found no compelling genetic evidence for color being under divergent selection. A time-calibrated mtDNA tree suggests that the present distribution of dorsal coloration resulted from processes occurring during the Pleistocene. Separate phylogenies based on SNPs and mtDNA resolved the same well supported clades, each containing different colored populations. Ancestral character state analysis provided some evidence for evolutionary transitions in color type. Genetic structure was more strongly associated with geographic features, than color category, suggesting that the distribution of color is explained by localized processes. Evidence for geographic isolation together with estimates of low effective population size implicates drift as playing a key role in color diversification. Our results highlight the relevance of considering the neutral processes involved with the evolution of traits with important fitness consequences.

Phylogeography and population genomics of a lotic water beetle across a complex tropical landscape

The habitat template concept applied to a freshwater system indicates that lotic species, or those which occupy permanent habitats along stream courses, are less dispersive than lentic species, or those that occur in more ephemeral aquatic habitats. Thus, populations of lotic species will be more structured than those of lentic species. Stream courses include both flowing water and small, stagnant microhabitats that can provide refuge when streams are low. Many species occur in these microhabitats but remain poorly studied. Here, we present population genetic data for one such species, the tropical diving beetle Exocelina manokwariensis (Dytiscidae), sampled from six localities along a ~300 km transect across the Birds Head Peninsula of New Guinea. Molecular data from both mitochondrial (CO1 sequences) and nuclear (ddRAD loci) regions document fine-scale population structure across populations that are ~45 km apart. Our results are concordant with previous phylogenetic and macroecological studies that applied the habitat template concept to aquatic systems. This study also illustrates that these diverse but mostly overlooked microhabitats are promising study systems in freshwater ecology and evolutionary biology. With the advent of next-generation sequencing, fine-scale population genomic studies are feasible for small nonmodel organisms to help illuminate the effect of habitat stability on species' natural history, population structure and geographic distribution.

Dealing with AFLP genotyping errors to reveal genetic structure in Plukenetia volubilis (Euphorbiaceae) in the Peruvian Amazon

An analysis of the population structure and genetic diversity for any organism often depends on one or more molecular marker techniques. Nonetheless, these techniques are not absolutely reliable because of various sources of errors arising during the genotyping process. Thus, a complex analysis of genotyping error was carried out with the AFLP method in 169 samples of the oil seed plant Plukenetia volubilis L. from small isolated subpopulations in the Peruvian Amazon. Samples were collected in nine localities from the region of San Martin. Analysis was done in eight datasets with a genotyping error from 0 to 5%. Using eleven primer combinations, 102 to 275 markers were obtained according to the dataset. It was found that it is only possible to obtain the most reliable and robust results through a multiple-level filtering process. Genotyping error and software set up influence both the estimation of population structure and genetic diversity, where in our case population number (K) varied between 2–9 depending on the dataset and statistical method used. Surprisingly, discrepancies in K number were caused more by statistical approaches than by genotyping errors themselves. However, for estimation of genetic diversity, the degree of genotyping error was critical because descriptive parameters (He, F_ST, PLP 5%) varied substantially (by at least 25%). Due to low gene flow, P. volubilis mostly consists of small isolated subpopulations (Φ_PT = 0.252–0.323) with some degree of admixture given by socio-economic connectivity among the sites; a direct link between the genetic and geographic distances was not confirmed. The study illustrates the successful application of AFLP to infer genetic structure in non-model plants.

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.

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.

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.

Selection on male sex pheromone composition contributes to butterfly reproductive isolation

Selection can facilitate diversification by inducing character displacement in mate choice traits that reduce the probability of maladaptive mating between lineages. Although reproductive character displacement (RCD) has been demonstrated in two-taxa case studies, the frequency of this process in nature is still debated. Moreover, studies have focused primarily on visual and acoustic traits, despite the fact that chemical communication is probably the most common means of species recognition. Here, we showed in a large, mostly sympatric, butterfly genus, a strong pattern of recurrent RCD for predicted male sex pheromone composition, but not for visual mate choice traits. Our results suggest that RCD is not anecdotal, and that selection for divergence in male sex pheromone composition contributed to reproductive isolation within the Bicyclus genus. We propose that selection may target olfactory mate choice traits as a more common sensory modality to ensure reproductive isolation among diverging lineages than previously envisaged.

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.

No predatory bias with respect to colour familiarity for the aposematic Adelphobates galactonotus (Anura: Dendrobatidae)

Aposematic colouration deters visually oriented predators because conspicuous signals are easier to detect and associate with unpalatability. Consequently, brightly coloured prey that are novel are predicted to be preyed on more than those with bright but typical colours. Here we evaluated whether predatory bias is associated with the colour differences observed at two different localities for a large, conspicuously coloured and poisonous Amazonian frog, Adelphobates galactonotus. At each locality predation experiments were carried out using frog models of two naturally occurring colours of the study species (blue and orange) and a control (brown). We found no evidence that novel colours were more vulnerable to predation than local colours. These results do not therefore support our hypothesis that predatory bias explains the geographic variation of colour in A. galactonotus.

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.

Coevolution and the diversification of life

Coevolution, reciprocal adaptation between two or more taxa, is commonly invoked as a primary mechanism responsible for generating much of Earth's biodiversity. This conceptually appealing hypothesis is incredibly broad in evolutionary scope, encompassing diverse patterns and processes operating over timescales ranging from microbial generations to geological eras. However, we have surprisingly little evidence that large-scale associations between coevolution and diversity reflect a causal relationship at smaller timescales, in which coevolutionary selection is directly responsible for the formation of new species. In this synthesis, we critically evaluate evidence for the often-invoked hypothesis that coevolution is an important process promoting biological diversification. We conclude that the lack of widespread evidence for coevolutionary diversification may be best explained by the fact that coevolution's importance in diversification varies depending on the type of interaction and the scale of the diversification under consideration.

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.