Sexual signal evolution and patterns of assortative mating across an intraspecific contact zone.

Contact zones provide important insights into the evolutionary processes that underlie lineage divergence and speciation. Here, we use a contact zone to ascertain speciation potential in the red-eyed treefrog (Agalychnis callidryas), a brightly coloured and polymorphic frog that exhibits unusually high levels of intraspecific variation. Populations of A. callidryas differ in a number of traits, several of which are known sexual signals that mediate premating reproductive isolation in allopatric populations. Along the Caribbean coast of Costa Rica, a ~100 km contact zone, situated between two phenotypically and genetically divergent parent populations, contains multiple colour pattern phenotypes and late-generation hybrids. This contact zone provides the opportunity to examine processes that are important in the earliest stages of lineage divergence. We performed analyses of colour pattern variation in five contact zone sites and six parental sites and found complex, continuous colour variation along the contact zone. We found discordance between the geographic distribution of colour pattern and previously described genomic population structure. We then used a parental site and contact zone site to measure assortative mating and directional selection from naturally-occurring amplectant mating pairs. We found assortative mating in a parental population, but no assortative mating in the contact zone. Furthermore, we uncovered evidence of directional preference towards the adjacent parental phenotype in the contact zone population, but no directional preference in the parent population. Combined, these data provide insights into potential dynamics at the contact zone borders and indicate that incipient speciation between parent populations will be slowed.

Genetic isolation by distance underlies colour pattern divergence in red-eyed treefrogs (Agalychnis callidryas).

Investigating the spatial distribution of genetic and phenotypic variation can provide insights into the evolutionary processes that shape diversity in natural systems. We characterized patterns of genetic and phenotypic diversity to learn about drivers of colour-pattern diversification in red-eyed treefrogs (Agalychnis callidryas) in Costa Rica. Along the Pacific coast, red-eyed treefrogs have conspicuous leg colour patterning that transitions from orange in the north to purple in the south. We measured phenotypic variation of frogs, with increased sampling at sites where the orange-to-purple transition occurs. At the transition zone, we discovered the co-occurrence of multiple colour-pattern morphs. To explore possible causes of this variation, we generated a single nucleotide polymorphism data set to analyse population genetic structure, measure genetic diversity and infer the processes that mediate genotype-phenotype dynamics. We investigated how patterns of genetic relatedness correspond to individual measures of colour pattern along the coast, including testing for the role of hybridization in geographic regions where orange and purple phenotypic groups co-occur. We found no evidence that colour-pattern polymorphism in the transition zone arose through recent hybridization. Instead, a strong pattern of genetic isolation by distance indicates that colour-pattern variation was either retained through other processes such as ancestral colour polymorphisms or ancient secondary contact, or else it was generated by novel mutations. We found that phenotype changes along the Pacific coast more than would be expected based on genetic divergence and geographic distance alone. Combined, our results suggest the possibility of selective pressures acting on colour pattern at a small geographic scale.

La investigación de la distribución espacial de la variación genética y fenotípica puede proporcionar información sobre los procesos evolutivos que dan forma a la diversidad en los sistemas naturales. Caracterizamos patrones de diversidad genética y fenotípica para conocer los impulsores de la diversificación de patrones de color en ranas con ojos rojos (Agalychnis callidryas) en Costa Rica. A lo largo de la costa del Pacífico, las ranas con ojos rojos tienen un patrón de color llamativo en las patas que cambia de naranja en el norte a púrpura en el sur. Medimos la variación fenotípica de las ranas en los sitios del Pacífico, con un mayor muestreo en los sitios donde ocurre la transición de naranja a púrpura. En la zona de transición, descubrimos la co-ocurrencia de múltiples morfos de patrones de color. Para explorar las posibles causas de esta variación, generamos un conjunto de datos SNP con secuenciación RAD para analizar la estructura genética de la población, medir la diversidad genética e inferir los procesos que median la dinámica genotipo-fenotipo. Investigamos cómo los patrones de parentesco genético se corresponden con medidas individuales de patrón de color a lo largo de la costa, incluidas las pruebas del papel de la hibridación en regiones geográficas donde coexisten grupos fenotípicos naranja y morado. No encontramos evidencia de que el polimorfismo del patrón de color en la zona de transición surgiera a través de una hibridación o introgresión reciente. En cambio, un fuerte patrón de aislamiento genético por distancia indica que la variación del patrón de color se retuvo a través de otros procesos, como los polimorfismos de color ancestrales, el contacto secundario antiguo o la generada por mutaciones novedosas. Descubrimos que el fenotipo de color cambia a lo largo de la costa del Pacífico más de lo que se esperaría solo por la divergencia genética y la distancia geográfica. Combinados, nuestros resultados sugieren la posibilidad de que las presiones selectivas actúen sobre el patrón de color a pequeña escala geográfica.

Natural and anthropogenic landscape factors shape functional connectivity of an ecological specialist in urban Southern California.

Identifying how natural (i.e., unaltered by human activity) and anthropogenic landscape variables influence contemporary functional connectivity in terrestrial organisms can elucidate the genetic consequences of environmental change. We examine population genetic structure and functional connectivity among populations of a declining species, the Blainville's horned lizard (Phrynosoma blainvillii), in the urbanized landscape of the Greater Los Angeles Area in Southern California, USA. Using single nucleotide polymorphism data, we assessed genetic structure among populations occurring at the interface of two abutting evolutionary lineages, and at a fine scale among habitat fragments within the heavily urbanized area. Based on the ecology of P. blainvillii, we predicted which environmental variables influence population structure and gene flow and used gravity models to distinguish among hypotheses to best explain population connectivity. Our results show evidence of admixture between two evolutionary lineages and strong population genetic structure across small habitat fragments. We also show that topography, microclimate, and soil and vegetation types are important predictors of functional connectivity, and that anthropogenic disturbance, including recent fire history and urban development, are key factors impacting contemporary population dynamics. Examining how natural and anthropogenic sources of landscape variation affect contemporary population genetics is critical to understanding how to best manage sensitive species in a rapidly changing landscape.

Genetic and phenotypic evidence of a contact zone between divergent colour morphs of the iconic red-eyed treefrog.

Hybrid zones act as natural laboratories where divergent genomes interact, providing powerful systems for examining the evolutionary processes underlying biological diversity. In this study, we characterized patterns of genomic and phenotypic variation resulting from hybridization between divergent intraspecific lineages of the Neotropical red-eyed treefrog (Agalychnis callidryas). We found genetic evidence of a newly discovered contact zone and phenotypic novelty in leg colour-a trait suspected to play a role in mediating assortative mating in this species. Analysis of hybrid ancestry revealed an abundance of later-generation Fn individuals, suggesting persistence of hybrids in the contact zone. Hybrids are predominantly of southern ancestry but are phenotypically more similar to northern populations. Genome-wide association mapping revealed QTL with measurable effects on leg-colour variation, but further work is required to dissect the architecture of this trait and establish causal links. Further, genomic cline analyses indicated substantial variation in patterns of introgression across the genome. Directional introgression of loci associated with different aspects of leg colour are inherited from each parental lineage, creating a distinct hybrid colour pattern. We show that hybridization can generate new phenotypes, revealing the evolutionary processes that potentially underlie patterns of phenotypic diversity in this iconic polytypic frog. Our study is consistent with a role of hybridization and sexual selection in lineage diversification, evolutionary processes that have been implicated in accelerating divergence in the most phenotypically diverse species.

Regional variation in drivers of connectivity for two frog species (Rana pretiosa and R. luteiventris) from the U.S. Pacific Northwest.

Comparative landscape genetics has uncovered high levels of variability in which landscape factors affect connectivity among species and regions. However, the relative importance of species traits versus environmental variation for predicting landscape patterns of connectivity is unresolved. We provide evidence from a landscape genetics study of two sister taxa of frogs, the Oregon spotted frog (Rana pretiosa) and the Columbia spotted frog (Rana luteiventris) in Oregon and Idaho, USA. Rana pretiosa is relatively more dependent on moisture for dispersal than R. luteiventris, so if species traits influence connectivity, we predicted that connectivity among R. pretiosa populations would be more positively associated with moisture than R. luteiventris. However, if environmental differences are important drivers of gene flow, we predicted that connectivity would be more positively related to moisture in arid regions. We tested these predictions using eight microsatellite loci and gravity models in two R. pretiosa regions and four R. luteiventris regions (n = 1,168 frogs). In R. pretiosa, but not R. luteiventris, connectivity was positively related to mean annual precipitation, supporting our first prediction. In contrast, connectivity was not more positively related to moisture in more arid regions. Various temperature metrics were important predictors for both species and in all regions, but the directionality of their effects varied. Therefore, the pattern of variation in drivers of connectivity was consistent with predictions based on species traits rather than on environmental variation.

Fire Does Not Strongly Affect Genetic Diversity or Structure of a Common Treefrog in the Endangered Florida Scrub.

Fire regimes influence natural populations of organisms in diverse ways, via direct effects on population dynamics as well as indirect effects on habitat and ecosystem processes. Although many amphibian species have evolved to persist in fire-dependent ecosystems, the effects of fire on the genetic diversity of amphibian populations remain relatively unexplored. We examined how different aspects of fire history relate to population genetic diversity and structure of an abundant anuran, Hyla femoralis, in a large, intact area of Florida scrub containing hundreds of seasonally inundated ponds. Specifically, we assessed the overall population genetic structure and examined whether variation in time since fire, fire intensity, or historical fire frequency at breeding sites explained spatial variation in genetic diversity. Based on our sampling of 17 breeding aggregations within the 2,100-ha study area, neither recent nor frequent fire reduce genetic diversity or restrict connectivity among ponds for H. femoralis. Overall, mean effective population sizes were large (average range = 68-572). We detected a positive trend between effective population size (Ne) and average intensity of the most-recent fire, with this factor explaining 42% of the variation in Ne. Our results contrast with previous studies that consistently demonstrate strong relationships between fire history and population genetic structure of scrub-associated lizard species, suggesting that H. femoralis is resilient to a wide range of fire regimes. More generally, our study contributes to understanding the roles of life-history characteristics and environmental unpredictability in shaping organisms' responses to fire.

Variation in male spermiation response to exogenous hormones among divergent populations of Red-eyed Treefrogs.

BACKGROUND: The non-lethal collection of sperm from live males is an important component for multiple captive-breeding techniques, including assisted reproductive technology (ART) protocols, sperm cryopreservation and in vitro fertilization. However, in amphibians, the type and amount of hormone necessary to induce spermiation can be highly variable, even among closely related species. We are unaware of any studies that have examined the spermiation response to exogenous hormones across highly differentiated populations within a species. METHODS: We examined variation in sperm viability and production in response to the hormone LHRH among four divergent populations of the Red-eyed Treefrog (Agalychnis callidryas). We hypothesized that these highly differentiated populations would show variability in sperm count and viability in response to two dosages, 2 µg/g and 4 µg/g, of the hormone LHRH. We collected spermic urine 3 h post injection (PI). We then examined variation in spermiation at 3, 7, 12, and 24 h PI of LHRH for two allopatric populations that previously showed evidence of premating behavioral isolation. RESULTS: One population of Red-eyed Treefrog exhibited reduced sperm viability but not count in response to the hormone LHRH compared to all other populations. In addition, we found peak viability at 3 h PI for the allopatric population comparison. There was no difference in sperm production within or between populations at 3, 7, 12, or 24 h PI. For both studies, intrapopulation variation was high. CONCLUSION: ART often focuses on threatened species with small, isolated populations, which could evolve localized differences due to the evolutionary process of drift and selection. The high variation in response and the population-level differences in sperm viability we observed demonstrate that practitioners of ART should consider the possibility of divergent responses to hormones which may affect study design and animal receptivity to ART protocols.

Minimal prezygotic isolation between ecologically divergent sibling species.

Divergence in mating signals typically accompanies speciation. We examine two ecologically divergent sibling species of crickets to assess the degree and timing of the evolution of prezygotic reproductive isolation. Gryllus saxatilis occurs in rocky habitats throughout western North America with long-winged individuals capable of long-distance dispersal; Gryllus navajo is endemic to red-rock sandstone areas of south-eastern Utah and north-eastern Arizona and has short-winged individuals only capable of limited dispersal. Previous genetic work suggested some degree of introgression and/or incomplete lineage sorting is likely. Here we: (1) use restriction site associated DNA sequencing (RAD-seq) genetic data to describe the degree of genetic divergence among species and populations; (2) examine the strength of prezygotic isolation by (i) quantifying differences among male mating songs, and (ii) testing whether females prefer G. saxatilis or G. navajo calling songs. Our results show that genetically distinct "pure" species populations and genetically intermediate populations exist. Male mating songs are statistically distinguishable, but the absolute differences are small. In playback experiments, females from pure populations had no preference based on song; however, females from a genetically intermediate population preferred G. navajo song. Together these results suggest that prezygotic isolation is minimal, and mediated by female behaviour in admixed populations.

Discordant patterns of evolutionary differentiation in two Neotropical treefrogs.

Comparative studies of codistributed taxa test the degree to which historical processes have shaped contemporary population structure. Discordant patterns of lineage divergence among taxa indicate that species differ in their response to common historical processes. The complex geologic landscape of the Isthmus of Central America provides an ideal setting to test the effects of vicariance and other biogeographic factors on population history. We compared divergence patterns between two codistributed Neotropical frogs (Dendropsophus ebraccatus and Agalychnis callidryas) that exhibit colour pattern polymorphisms among populations, and found significant differences between them in phenotypic and genetic divergence among populations. Colour pattern in D. ebraccatus did not vary with genetic or geographic distance, while colour pattern co-varied with patterns of gene flow in A. callidryas. In addition, we detected significant species differences in the phylogenetic history of populations, gene flow among them, and the extent to which historical diversification and recent gene flow have been restricted by five biogeographic barriers in Costa Rica and Panama. We inferred that alternate microevolutionary processes explain the unique patterns of diversification in each taxon. Our study underscores how differences in selective regimes and species-typical ecological and life-history traits maintain spatial patterns of diversification.

Genetic diversification, vicariance, and selection in a polytypic frog.

Spatial patterns of heritable phenotypic diversity reflect the relative roles of gene flow and selection in determining geographic variation within a species. We quantified color differentiation and genetic divergence among 20 populations of the red-eyed tree frog (Agalychnis callidryas) in lower Central America. Phylogenetic analyses revealed 5 well-supported mitochondrial DNA clades, and we infer from our phylogeny that geographic barriers have played a large role in structuring populations. Two phenotypic characters varied independently among isolated population groups: Flank coloration distinguished Caribbean from Pacific individuals, whereas leg coloration exhibited a more complex geographic pattern. We detected 3 generalized spatial patterns of genetic and phenotypic diversity: 1) phenotypic differentiation in the presence of historical connectivity, 2) phenotypic uniformity across genetically differentiated regions, and 3) codistribution of genetic and phenotypic characters. These patterns indicate that phenotypic diversification is highly regionalized and can result from spatial variation in localized adaptations, geographic isolation, genetic drift, and/or evolutionary stasis. Although the mode of selection underlying color variation was not the focal objective of this study, we discuss the possible roles of natural and sexual selection in mediating population differentiation. Our study underscores the fact that selection gradients vary across relatively small spatial scales, even in species that occupy relatively homogeneous environments.

Fine scale gene flow and individual movements among subpopulations of Centrolene prosoblepon (Anura: Centrolenidae).

Dispersal capabilities determine and maintain local gene flow, and this has implications for population persistence and/or recolonization following environmental perturbations (natural or anthropogenic), disease outbreaks, or other demographic collapses. To predict recolonization and understand dispersal capacity in a stream-breeding frog, we examined individual movement patterns and gene flow among four subpopulations of the Neotropical glassfrog, Centrolene prosoblepon, at a mid-elevation cloud forest site at El Copé, Panama. We measured male movement directly during a two year mark-recapture study, and indirectly with gene flow estimates from mitochondrial DNA sequences (mtDNA). Individuals of this species showed strong site fidelity: over two years, male frogs in all four headwater streams moved very little (mean = 2.33 m; mode = 0 m). Nine individuals changed streams within one or two years, moving 675-1,108 m. For those males moving more than 10 m, movement was biased upstream (p < 0.001). Using mtDNA ND1 gene sequences, we quantified gene flow within and among headwater streams at two spatial scales: among headwater streams within two adjacent watersheds (2.5 km2) and among streams within a longitudinal gradient covering 5.0 km2. We found high gene flow among headwater streams (phi(ST) = 0.007, p = 0.325) but gene flow was more limited across greater distances (phi(CT) = 0.322, p = 0.065), even within the same drainage network. Lowland populations of C. prosoblepon potentially act as an important source of colonists for upland populations in this watershed.

Effects of changing climate on aquatic habitat and connectivity for remnant populations of a wide-ranging frog species in an arid landscape.

Amphibian species persisting in isolated streams and wetlands in desert environments can be susceptible to low connectivity, genetic isolation, and climate changes. We evaluated the past (1900-1930), recent (1981-2010), and future (2071-2100) climate suitability of the arid Great Basin (USA) for the Columbia spotted frog (Rana luteiventris) and assessed whether changes in surface water may affect connectivity for remaining populations. We developed a predictive model of current climate suitability and used it to predict the historic and future distribution of suitable climates. We then modeled changes in surface water availability at each time period. Finally, we quantified connectivity among existing populations on the basis of hydrology and correlated it with interpopulation genetic distance. We found that the area of the Great Basin with suitable climate conditions has declined by approximately 49% over the last century and will likely continue to decline under future climate scenarios. Climate conditions at currently occupied locations have been relatively stable over the last century, which may explain persistence at these sites. However, future climates at these currently occupied locations are predicted to become warmer throughout the year and drier during the frog's activity period (May - September). Fall and winter precipitation may increase, but as rain instead of snow. Earlier runoff and lower summer base flows may reduce connectivity between neighboring populations, which is already limited. Many of these changes could have negative effects on remaining populations over the next 50-80 years, but milder winters, longer growing seasons, and wetter falls might positively affect survival and dispersal. Collectively, however, seasonal shifts in temperature, precipitation, and stream flow patterns could reduce habitat suitability and connectivity for frogs and possibly other aquatic species inhabiting streams in this arid region.

Fine scale gene flow and individual movements among subpopulations of Centrolene prosoblepon (Anura: Centrolenidae)

Dispersal capabilities determine and maintain local gene flow, and this has implications for population persistence and/or recolonization following environmental perturbations (natural or anthropogenic), disease outbreaks, or other demographic collapses. To predict recolonization and understand dispersal capacity in a stream-breeding frog, we examined individual movement patterns and gene flow among four subpopulations of the Neotropical glassfrog, Centrolene prosoblepon, at a mid-elevation cloud forest site at El Copé, Panama. We measured male movement directly during a two year mark-recapture study, and indirectly with gene flow estimates from mitochondrial DNA sequences (mtDNA). Individuals of this species showed strong site fidelity: over two years, male frogs in all four headwater streams moved very little (mean = 2.33 m; mode = 0 m). Nine individuals changed streams within one or two years, moving 675-1,108 m. For those males moving more than 10 m, movement was biased upstream (p < 0.001). Using mtDNA ND1 gene sequences, we quantified gene flow within and among headwater streams at two spatial scales: among headwater streams within two adjacent watersheds (2.5 km2) and among streams within a longitudinal gradient covering 5.0 km2. We found high gene flow among headwater streams (phi(ST) = 0.007, p = 0.325) but gene flow was more limited across greater distances (phi(CT) = 0.322, p = 0.065), even within the same drainage network. Lowland populations of C. prosoblepon potentially act as an important source of colonists for upland populations in this watershed.

La capacidad de dispersión determina y mantiene el flujo genético local, y esto tiene implicaciones para la persistencia poblacional y/o la recolonización que sigue a perturbaciones ambientales. Examinamos patrones individuales de movimiento y flujo genético entre subpoblaciones de Centrolene prosoblepon (Anura: Centrolenidae) en un sitio de elevación media en El Copé, Panamá. Medimos directamente el movimiento de los machos durante un estudio de marcado-recaptura, e indirectamente con estimaciones de flujo genético a partir de secuencias de ADN mitocondrial (mtDNA). Los individuos mostraron fuerte fidelidad a su lugar: por más de dos años, las ranas macho de los cuatro arroyos al inicio del río se movieron muy poco (promedio = 2.33 m; moda = 0 m). Nueve individuos cambiaron de corriente de agua en uno o dos años, moviéndose 675-1 108 m. Usando la secuencia genética ND1 del ADN mitocondrial, medimos el flujo genético en dos escalas espaciales: entre arroyos que originan el río (2.5 km2) y entre arroyos con un gradiente longitudinal en 5.0 km2. Encontramos un flujo genético alto entre los arroyos al inicio del río (f = 0.007, p = 0.325 y otro más limitado en distancias mayores (f = 0.322, p = 0.065).