R. Graham Reynolds

Assortative mating in poison-dart frogs based on an ecologically important trait

Abstract The origin of new species can be influenced by both deterministic and stochastic factors. Mate choice and natural selection may be important deterministic causes of speciation (as opposed to the essentially stochastic factors of geographic isolation and genetic drift). Theoretical models predict that speciation is more likely when mate choice depends on an ecologically important trait that is subject to divergent natural selection, although many authors have considered such mating/ecology pleiotropy, or “magic-traits” to be unlikely. However, phenotypic signals are important in both mate choice and ecological processes such as avoiding predation. In chemically defended species, it may be that the phenotypic characteristics influencing mate choice are the same signals being used to transmit a warning to potential predators, although few studies have demonstrated this in wild populations. We tested for assortative mating between two color morphs of the Strawberry Poison-Dart Frog, Dendrobates pumilio, a group with striking geographic variation in aposematic color patterns. We found that females significantly prefer individuals of their own morph under two different light treatments, indicating strong assortative mating based on multiple coloration cues that are also important ecological signals. This study provides a rare example of one phenotypic trait affecting both ecological viability and nonrandom mating, indicating that mating/ecology pleiotropy is plausible in wild populations, particularly for organisms that are aposematically colored and visually orienting.

What can DNA tell us about biological invasions

It is often hoped that population genetics can answer questions about the demographic and geographic dynamics of recent biological invasions. Conversely, invasions with well-known histories are sometimes billed as opportunities to test methods of population genetic inference. In both cases, underappreciated limitations constrain the usefulness of genetic methods. The most significant is that human-caused invasions have occurred on historical timescales that are orders of magnitude smaller than the timescales of mutation and genetic drift for most multicellular organisms. Analyses based on the neutral theory of molecular evolution cannot resolve such rapid dynamics. Invasion histories cannot be reconstructed in the same way as biogeographic changes occurring over millenia. Analyses assuming equilibrium between mutation, drift, gene flow, and selection will rarely be applicable, and even methods designed for explicitly non-equilibrium questions often require longer timescales than the few generations of most invasions of current concern. We identified only a few population genetic questions that are tractable over such short timescales. These include comparison of alternative hypotheses for the geographic origin of an invasion, testing for bottlenecks, and hybridization between native and invasive species. When proposing population genetic analysis of a biological invasion, we recommend that biologists ask (i) whether the questions to be addressed will materially affect management practice or policy, and (ii) whether the proposed analyses can really be expected to address important population genetic questions. Despite our own enthusiasm for population genetic research, we conclude that genetic analysis of biological invasions is justified only under exceptional circumstances.

Toward a Tree-of-Life for the boas and pythons: Multilocus species-level phylogeny with unprecedented taxon sampling

Snakes in the families Boidae and Pythonidae constitute some of the most spectacular reptiles and comprise an enormous diversity of morphology, behavior, and ecology. While many species of boas and pythons are familiar, taxonomy and evolutionary relationships within these families remain contentious and fluid. A major effort in evolutionary and conservation biology is to assemble a comprehensive Tree-of-Life, or a macro-scale phylogenetic hypothesis, for all known life on Earth. No previously published study has produced a species-level molecular phylogeny for more than 61% of boa species or 65% of python species. Using both novel and previously published sequence data, we have produced a species-level phylogeny for 84.5% of boid species and 82.5% of pythonid species, contextualized within a larger phylogeny of henophidian snakes. We obtained new sequence data for three boid, one pythonid, and two tropidophiid taxa which have never previously been included in a molecular study, in addition to generating novel sequences for seven genes across an additional 12 taxa. We compiled an 11-gene dataset for 127 taxa, consisting of the mitochondrial genes CYTB, 12S, and 16S, and the nuclear genes bdnf, bmp2, c-mos, gpr35, rag1, ntf3, odc, and slc30a1, totaling up to 7561 base pairs per taxon. We analyzed this dataset using both maximum likelihood and Bayesian inference and recovered a well-supported phylogeny for these species. We found significant evidence of discordance between taxonomy and evolutionary relationships in the genera Tropidophis, Morelia, Liasis, and Leiopython, and we found support for elevating two previously suggested boid species. We suggest a revised taxonomy for the boas (13 genera, 58 species) and pythons (8 genera, 40 species), review relationships between our study and the many other molecular phylogenetic studies of henophidian snakes, and present a taxonomic database and alignment which may be easily used and built upon by other researchers.

A NEW BAYESIAN METHOD FOR FITTING EVOLUTIONARY MODELS TO COMPARATIVE DATA WITH INTRASPECIFIC VARIATION

Phylogenetic comparative methods that incorporate intraspecific variability are relatively new and, so far, not especially widely used in empirical studies. In the present short article we will describe a new Bayesian method for fitting evolutionary models to comparative data that incorporates intraspecific variability. This method differs from an existing likelihood‐based approach in that it requires no a priori inference about species means and variances; rather it takes phenotypic values from individuals and a phylogenetic tree as input, and then samples species means and variances, along with the parameters of the evolutionary model, from their joint posterior probability distribution. One of the most novel and intriguing attributes of this approach is that jointly sampling the species means with the evolutionary model parameters means that the model and tree can influence our estimates of species mean trait values, not just the reverse. In the present implementation, we first apply this method to the most widely used evolutionary model for continuously valued phenotypic trait data (Brownian motion). However, the general approach has broad applicability, which we illustrate by also fitting the λ model, another simple model for quantitative trait evolution on a phylogeny. We test our approach via simulation and by analyzing two empirical datasets obtained from the literature. Finally, we have implemented the methods described herein in a new function for the R statistical computing environment, and this function will be distributed as part of the ‘phytools’ R library.

Phenotypic shifts in urban areas in the tropical lizard Anolis cristatellus

Urbanization is an increasingly important dimension of global change, and urban areas likely impose significant natural selection on the species that reside within them. Although many species of plants and animals can survive in urban areas, so far relatively little research has investigated whether such populations have adapted (in an evolutionary sense) to their newfound milieu. Even less of this work has taken place in tropical regions, many of which have experienced dramatic growth and intensification of urbanization in recent decades. In the present study, we focus on the neotropical lizard, Anolis cristatellus. We tested whether lizard ecology and morphology differ between urban and natural areas in three of the most populous municipalities on the island of Puerto Rico. We found that environmental conditions including temperature, humidity, and substrate availability differ dramatically between neighboring urban and natural areas. We also found that lizards in urban areas use artificial substrates a large proportion of the time, and that these substrates tend to be broader than substrates in natural forest. Finally, our morphological data showed that lizards in urban areas have longer limbs relative to their body size, as well as more subdigital scales called lamellae, when compared to lizards from nearby forested habitats. This shift in phenotype is exactly in the direction predicted based on habitat differences between our urban and natural study sites, combined with our results on how substrates are being used by lizards in these areas. Findings from a common‐garden rearing experiment using individuals from one of our three pairs of populations provide evidence that trait differences between urban and natural sites may be genetically based. Taken together, our data suggest that anoles in urban areas are under significant differential natural selection and may be evolutionarily adapting to their human‐modified environments.

Effects of climatic and geological processes during the pleistocene on the evolutionary history of the northern cavefish, Amblyopsis spelaea (teleostei: amblyopsidae).

Climatic and geological processes associated with glaciation cycles during the Pleistocene have been implicated in influencing patterns of genetic variation and promoting speciation of temperate flora and fauna. However, determining the factors promoting divergence and speciation is often difficult in many groups because of our limited understanding of potential vicariant barriers and connectivity between populations. Pleistocene glacial cycles are thought to have significantly influenced the distribution and diversity of subterranean invertebrates; however, impacts on subterranean aquatic vertebrates are less clear. We employed several hypothesis‐driven approaches to assess the impacts of Pleistocene climatic and geological changes on the Northern Cavefish, Amblyopsis spelaea, whose current distribution occurs near the southern extent of glacial advances in North America. Our results show that the modern Ohio River has been a significant barrier to dispersal and is correlated with patterns of genetic divergence. We infer that populations were isolated in two refugia located north and south of the Ohio River during the most recent two glacial cycles with evidence of demographic expansion in the northern isolate. Finally, we conclude that climatic and geological processes have resulted in the formation of cryptic forms and advocate recognition of two distinct phylogenetic lineages currently recognized as A. spelaea.

Consecutive Virgin Births in the New World Boid Snake, the Colombian Rainbow Boa, Epicrates maurus

Until recently, facultative automictic parthenogenesis within the squamate reptiles exhibiting ZZ:ZW genetic sex determination has resulted in single reproductive events producing male (ZZ) or female (ZW) offspring. With the recent discovery of viable parthenogenetically produced female (WW) Boa constrictors, the existence of further parthenogenetic events resulting in WW females was questioned. Here, we provide genetic evidence for consecutive virgin births by a female Colombian rainbow boa (Epicrates maurus), resulting in the production of WW females likely through terminal fusion automixis. Samples were screened at 22 microsatellite loci with 12 amplifying unambiguous products. Of these, maternal heterozygosity was observed in 4, with the offspring differentially homozygous at each locus. This study documents the first record of parthenogenesis within the genus Epicrates, a second within the serpent lineage Boidae, and the third genetically confirmed case of consecutive virgin births of viable offspring within any vertebrate lineage. Unlike the recent record in Boa constrictors, the female described here was isolated from conspecifics from birth, demonstrating that males are not required to stimulate parthenogenetic reproduction in this species and possibly other Boas.

Molecular phylogeny and historical biogeography of West Indian boid snakes (Chilabothrus)

The evolutionary and biogeographic history of West Indian boid snakes (Epicrates), a group of nine species and 14 subspecies, was once thought to be well understood; however, new research has indicated that we are missing a clear understanding of the evolutionary relationships of this group. Here, we present the first multilocus, species-tree based analyses of the evolutionary relationships, divergence times, and historical biogeography of this clade with data from 10 genes and 6256 bp. We find evidence for a single colonization of the Caribbean from mainland South America in the Oligocene or early Miocene, followed by a radiation throughout the Greater Antilles and Bahamas. These findings support the previous suggestion that Epicrates sensu lato Wagler is paraphyletic with respect to the anacondas (Eunectes Wagler), and hence we restrict Epicrates to the mainland clade and use the available name Chilabothrus Duméril and Bibron for the West Indian clade. Our results suggest some diversification occurred within island banks, though most species divergence events seem to have occurred in allopatry. We also find evidence for a remarkable diversification within the Bahamian archipelago suggesting that the recognition of another Bahamian endemic species C. strigilatus is warranted.

A taxonomic revision of boas (Serpentes: Boidae).

Large molecular datasets including many species and loci have greatly improved our knowledge of snake phylogeny, particularly within the group including boas (Table 1). Recent taxonomic revisions using molecular phylogenies have clarified some of the previously contentious nomenclature of the group (Wilcox et al. 2002; Lawson et al. 2004; Burbrink 2005; Noonan & Chippindale 2006), resulting in a robust taxonomy that is mostly concordant with the phylogeny as currently known, which includes ~85% of described, extant species (Pyron et al. 2013; Reynolds et al. 2014). However, a few unresolved issues remain, related primarily to the rules of the International Code of Zoological Nomenclature ( the Code hereafter) and the application of Linnaean ranks (International Commission on Zoological Nomenclature et al. 1999).

Genetic analysis of a novel invasion of Puerto Rico by an exotic constricting snake

The tropical island Puerto Rico is potentially vulnerable to invasion by some species of exotic snakes; however, until now no established populations had been reported. Here we report and genetically characterize the nascent invasion of Puerto Rico by an exotic constricting snake of the family Boidae (Boa constrictor) using mtDNA and microsatellite data. Over 150 individual B. constrictor have been removed from Mayagüez municipality since May 2011, and our results from the genetic analysis of 32 individuals suggest that this population was recently founded by individuals of one subspecies from a genetic lineage common to zoo and breeding collections, but that the potential propagule pool consists of two subspecies. We also suggest that anthropogenic long-distance dispersal within the island of Puerto Rico may be occurring from the established population, with implications for further establishment across the island. This study represents the first report of the naturalization of an invasive species of boid snake in Puerto Rico and will be important in determining mitigation strategies for this invasion as well as providing a basis for comparison to other on-going studies of invasive snakes.