Genome-wide parallelism underlies contemporary adaptation in urban lizards.

Urbanization drastically transforms landscapes, resulting in fragmentation, degradation, and the loss of local biodiversity. Yet, urban environments also offer opportunities to observe rapid evolutionary change in wild populations that survive and even thrive in these novel habitats. In many ways, cities represent replicated "natural experiments" in which geographically separated populations adaptively respond to similar selection pressures over rapid evolutionary timescales. Little is known, however, about the genetic basis of adaptive phenotypic differentiation in urban populations nor the extent to which phenotypic parallelism is reflected at the genomic level with signatures of parallel selection. Here, we analyzed the genomic underpinnings of parallel urban-associated phenotypic change in Anolis cristatellus, a small-bodied neotropical lizard found abundantly in both urbanized and forested environments. We show that phenotypic parallelism in response to parallel urban environmental change is underlain by genomic parallelism and identify candidate loci across the Anolis genome associated with this adaptive morphological divergence. Our findings point to polygenic selection on standing genetic variation as a key process to effectuate rapid morphological adaptation. Identified candidate loci represent several functions associated with skeletomuscular development, morphology, and human disease. Taken together, these results shed light on the genomic basis of complex morphological adaptations, provide insight into the role of contingency and determinism in adaptation to novel environments, and underscore the value of urban environments to address fundamental evolutionary questions.

Testing for heterogeneous rates of discrete character evolution on phylogenies.

Many hypotheses in the field of phylogenetic comparative biology involve specific changes in the rate or process of trait evolution. This is particularly true of approaches designed to connect macroevolutionary pattern to microevolutionary process. We present a method designed to test whether the rate of evolution of a discrete character has changed in one or more clades, lineages, or time periods. This method differs from other related approaches (such as the 'covarion' model) in that the 'regimes' in which the rate or process is postulated to have changed are specified a priori by the user, rather than inferred from the data. Similarly, it differs from methods designed to model a correlation between two binary traits in that the regimes mapped onto the tree are fixed. We apply our method to investigate the rate of dewlap color and/or caudal vertebra number evolution in Caribbean and mainland clades of the diverse lizard genus Anolis. We find little evidence to support any difference in the evolutionary process between mainland and island evolution for either character. We also examine the statistical properties of the method more generally and show that it has acceptable type I error, parameter estimation, and power. Finally, we discuss some general issues of frequentist hypothesis testing and model adequacy, as well as the relationship of our method to existing models of heterogeneity in the rate of discrete character evolution on phylogenies.

The Difficulty of Predicting Evolutionary Change in Response to Novel Ecological Interactions: A Field Experiment with Anolis Lizards.

AbstractDetermining whether and how evolution is predictable is an important goal, particularly as anthropogenic disturbances lead to novel species interactions that could modify selective pressures. Here, we use a multigeneration field experiment with brown anole lizards (Anolis sagrei) to test hypotheses about the predictability of evolution. We manipulated the presence/absence of predators and competitors of A. sagrei across 16 islands in the Bahamas that had preexisting brown anole populations. Before the experiment and again after roughly five generations, we measured traits related to locomotor performance and habitat use by brown anoles and used double-digest restriction enzyme-associated DNA sequencing to estimate genome-wide changes in allele frequencies. Although previous work showed that predators and competitors had characteristic effects on brown anole behavior, diet, and population sizes, we found that evolutionary change at both phenotypic and genomic levels was difficult to forecast. Phenotypic changes were contingent on sex and habitat use, whereas genetic change was unpredictable and not measurably correlated with phenotypic changes, experimental treatments, or other environmental factors. Our work shows how differences in ecological context can alter evolutionary outcomes over short timescales and underscores the difficulty of forecasting evolutionary responses to multispecies interactions in natural conditions, even in a well-studied system with ample supporting ecological information.

Mesozoic origin of coleoid cephalopods and their abrupt shifts of diversification patterns.

Coleoids are the most diverse group of cephalopod mollusks. While their origin is date during the Mesozoic, the diversification pattern is unknown. However, two hypotheses have been proposed. The first suggests an increasing diversification rate after the Cretaceous-Paleogene extinction event (K-Pg) as consequence of empty habitats left by the ammonites and belemnites. The second hypothesis proposes a mid-Cenozoic increase in diversification rate related to distributional changes during ice ages and biotic interactions. To test these hypotheses, we estimated a lineage through time (LTT) and the gamma-statistic along with model-based diversification rates. These analyses were conducted on a dated molecular phylogeny for coleoids that we reconstructed using five molecular markers (cytochrome b, 16S rRNA, cytochrome oxidase I, rhodopsin, and PAX-6). Our divergence time estimation suggests that coleoids originated in the Mesozoic Era (Middle Triassic) and that both main clades (Decapodiformes and Octopodiformes) diverged in the Cretaceous/Jurassic Period. The LTT, gamma statistic, and diversification rates inferred with the Bayesian Analysis of Macro-evolutionary Mixtures (BAMM), indicate an acceleration in diversification rate over time since the origin of coleoids. Additionally, BAMM allowed us to detect abrupt increases in diversification rate before and after the K-Pg boundary. Our results partially support both hypotheses as all analyses indicate that the coleoid diversification rate was increasing during the Cenozoic. However, our results also indicate increasing diversification rates before the K-Pg boundary. We propose that the radiation of coleoids has been shaped by an acceleration in diversification rate over time, including exceptional episodes of abrupt increases before and after the K-Pg boundary.

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.

Correlated evolution of flower size and seed number in flowering plants (monocotyledons).

Background and Aims: Kin selection theory predicts that a parent may minimize deleterious effects of competition among seeds developing within ovaries by increasing the genetic relatedness of seeds within an ovary. Alternatively, the number of developing seeds could be reduced to one or a few. It has also been suggested that single or few seeded fruits may be correlated with small flowers, and multi-ovulate ovaries or many seeded fruits may be associated with large flowers with specialized pollination mechanisms. We examined the correlation between flower size and seed number in 69 families of monocotyledons to assess if correlations are significant and independent of phylogeny. Methods: We first examined the effect of phylogenetic history on the evolution of these two traits, flower size and seed number, and then mapped correlations between them on the latest phylogenetic tree of monocotyledons. Results: The results provide phylogenetically robust evidence of strong correlated evolution between flower size and seed number and show that correlated evolution of traits is not constrained by phylogenetic history of taxa. Moreover, the two character combinations, small flowers and a single or few seeds per fruit, and large flowers and many seeded fruits, have persisted in monocotyledons longer than other trait combinations. Conclusions: The analyses support the suggestion that most angiosperms may fall into two categories, one with large flowers and many seeded fruits and the other with small flowers and single or few seeded fruits, and kin selection within ovaries may explain the observed patterns.

Linking locomotor performance to morphological shifts in urban lizards.

Urban habitats are drastically modified from their natural state, creating unique challenges and selection pressures for organisms that reside in them. We compared locomotor performance of Anolis lizards from urban and forest habitats on tracks differing in angle and substrate, and found that using artificial substrates came at a cost: lizards ran substantially slower and frequently lost traction on man-made surfaces compared to bark. We found that various morphological traits were positively correlated with sprint speed and that these same traits were significantly larger in urban compared to forest lizards. We found that urban lizards ran faster on both man-made and natural surfaces, suggesting similar mechanisms improve locomotor performance on both classes of substrate. Thus, lizards in urban areas may be under selection to run faster on all flat surfaces, while forest lizards face competing demands of running, jumping and clinging to narrow perches. Novel locomotor challenges posed by urban habitats likely have fitness consequences for lizards that cannot effectively use man-made surfaces, providing a mechanistic basis for observed phenotypic shifts in urban populations of this species.

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.

phytools 2.0: an updated R ecosystem for phylogenetic comparative methods (and other things).

Phylogenetic comparative methods comprise the general endeavor of using an estimated phylogenetic tree (or set of trees) to make secondary inferences: about trait evolution, diversification dynamics, biogeography, community ecology, and a wide range of other phenomena or processes. Over the past ten years or so, the phytools R package has grown to become an important research tool for phylogenetic comparative analysis. phytools is a diverse contributed R library now consisting of hundreds of different functions covering a variety of methods and purposes in phylogenetic biology. As of the time of writing, phytools included functionality for fitting models of trait evolution, for reconstructing ancestral states, for studying diversification on trees, and for visualizing phylogenies, comparative data, and fitted models, as well numerous other tasks related to phylogenetic biology. Here, I describe some significant features of and recent updates to phytools, while also illustrating several popular workflows of the phytools computational software.

Divergence in coloration and ecological speciation in the Anolis marmoratus species complex.

Adaptive divergence in coloration is expected to produce reproductive isolation in species that use colourful signals in mate choice and species recognition. Indeed, many adaptive radiations are characterized by differentiation in colourful signals, suggesting that divergent selection acting on coloration may be an important component of speciation. Populations in the Anolis marmoratus species complex from the Caribbean island of Guadeloupe display striking divergence in the colour and pattern of adult males that occurs over small geographic distances, suggesting strong divergent selection. Here we test the hypothesis that divergence in coloration results in reduced gene flow among populations. We quantify variation in adult male coloration across a habitat gradient between mesic and xeric habitats, use a multilocus coalescent approach to infer historical demographic parameters of divergence, and examine gene flow and population structure using microsatellite variation. We find that colour variation evolved without geographic isolation and in the face of gene flow, consistent with strong divergent selection and that both ecological and sexual selection are implicated. However, we find no significant differentiation at microsatellite loci across populations, suggesting little reproductive isolation and high levels of contemporary gene exchange. Strong divergent selection on loci affecting coloration probably maintains clinal phenotypic variation despite high gene flow at neutral loci, supporting the notion of a porous genome in which adaptive portions of the genome remain fixed whereas neutral portions are homogenized by gene flow and recombination. We discuss the impact of these findings for studies of colour evolution and ecological speciation.

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.

Repeated modification of early limb morphogenesis programmes underlies the convergence of relative limb length in Anolis lizards.

The independent evolution of similar morphologies has long been a subject of considerable interest to biologists. Does phenotypic convergence reflect the primacy of natural selection, or does development set the course of evolution by channelling variation in certain directions? Here, we examine the ontogenetic origins of relative limb length variation among Anolis lizard habitat specialists to address whether convergent phenotypes have arisen through convergent developmental trajectories. Despite the numerous developmental processes that could potentially contribute to variation in adult limb length, our analyses reveal that, in Anolis lizards, such variation is repeatedly the result of changes occurring very early in development, prior to formation of the cartilaginous long bone anlagen.

A simple hierarchical model for heterogeneity in the evolutionary correlation on a phylogenetic tree.

Numerous questions in phylogenetic comparative biology revolve around the correlated evolution of two or more phenotypic traits on a phylogeny. In many cases, it may be sufficient to assume a constant value for the evolutionary correlation between characters across all the clades and branches of the tree. Under other circumstances, however, it is desirable or necessary to account for the possibility that the evolutionary correlation differs through time or in different sections of the phylogeny. Here, we present a method designed to fit a hierarchical series of models for heterogeneity in the evolutionary rates and correlation of two quantitative traits on a phylogenetic tree. We apply the method to two datasets: one for different attributes of the buccal morphology in sunfishes (Centrarchidae); and a second for overall body length and relative body depth in rock- and non-rock-dwelling South American iguanian lizards. We also examine the performance of the method for parameter estimation and model selection using a small set of numerical simulations.

Testing for genetic assimilation with phylogenetic comparative analysis: Conceptual, methodological, and statistical considerations.

Genetic assimilation is a process that leads to reduced phenotypic plasticity during adaptation to novel conditions, a potentially important phenomenon under global environmental change. Null expectations when testing for genetic assimilation, however, are not always clear. For instance, the statistical artifact of regression to the mean could bias us toward detecting genetic assimilation when it has not occurred. Likewise, the specific mechanism underlying plasticity expression may affect null expectations under neutral evolution. We used macroevolutionary numerical simulations to examine both of these important issues and their interaction, varying whether plasticity evolves, the evolutionary mechanism, trait measurement error, and experimental design. We also modified an existing reaction norm correction method to account for phylogenetic nonindependence. We found (1) regression to the mean is pervasive and can generate spurious support for genetic assimilation; (2) experimental design and post hoc correction can minimize this spurious effect; and (3) neutral evolution can produce patterns consistent with genetic assimilation without constraint or selection, depending on the mechanism of plasticity expression. Additionally, we reanalyzed published macroevolutionary data supporting genetic assimilation, and found that support was reduced after proper correction. Considerable caution is thus required whenever investigating genetic assimilation and reaction norm evolution at macroevolutionary scales.

A variable-rate quantitative trait evolution model using penalized-likelihood.

In recent years it has become increasingly popular to use phylogenetic comparative methods to investigate heterogeneity in the rate or process of quantitative trait evolution across the branches or clades of a phylogenetic tree. Here, I present a new method for modeling variability in the rate of evolution of a continuously-valued character trait on a reconstructed phylogeny. The underlying model of evolution is stochastic diffusion (Brownian motion), but in which the instantaneous diffusion rate (σ2) also evolves by Brownian motion on a logarithmic scale. Unfortunately, it's not possible to simultaneously estimate the rates of evolution along each edge of the tree and the rate of evolution of σ2 itself using Maximum Likelihood. As such, I propose a penalized-likelihood method in which the penalty term is equal to the log-transformed probability density of the rates under a Brownian model, multiplied by a 'smoothing' coefficient, λ, selected by the user. λ determines the magnitude of penalty that's applied to rate variation between edges. Lower values of λ penalize rate variation relatively little; whereas larger λ values result in minimal rate variation among edges of the tree in the fitted model, eventually converging on a single value of σ2 for all of the branches of the tree. In addition to presenting this model here, I have also implemented it as part of my phytools R package in the function multirateBM. Using different values of the penalty coefficient, λ, I fit the model to simulated data with: Brownian rate variation among edges (the model assumption); uncorrelated rate variation; rate changes that occur in discrete places on the tree; and no rate variation at all among the branches of the phylogeny. I then compare the estimated values of σ2 to their known true values. In addition, I use the method to analyze a simple empirical dataset of body mass evolution in mammals. Finally, I discuss the relationship between the method of this article and other models from the phylogenetic comparative methods and finance literature, as well as some applications and limitations of the approach.

covid19.Explorer: a web application and R package to explore United States COVID-19 data.

Appearing at the end of 2019, a novel virus (later identified as SARS-CoV-2) was characterized in the city of Wuhan in Hubei Province, China. As of the time of writing, the disease caused by this virus (known as COVID-19) has already resulted in over three million deaths worldwide. SARS-CoV-2 infections and deaths, however, have been highly unevenly distributed among age groups, sexes, countries, and jurisdictions over the course of the pandemic. Herein, I present a tool (the covid19.Explorer R package and web application) that has been designed to explore and analyze publicly available United States COVID-19 infection and death data from the 2020/21 U.S. SARS-CoV-2 pandemic. The analyses and visualizations that this R package and web application facilitate can help users better comprehend the geographic progress of the pandemic, the effectiveness of non-pharmaceutical interventions (such as lockdowns and other measures, which have varied widely among U.S. states), and the relative risks posed by COVID-19 to different age groups within the U.S. population. The end result is an interactive tool that will help its users develop an improved understanding of the temporal and geographic dynamics of the SARS-CoV-2 pandemic, accessible to lay people and scientists alike.

Early giant reveals faster evolution of large body size in ichthyosaurs than in cetaceans.

Body sizes of marine amniotes span six orders of magnitude, yet the factors that governed the evolution of this diversity are largely unknown. High primary production of modern oceans is considered a prerequisite for the emergence of cetacean giants, but that condition cannot explain gigantism in Triassic ichthyosaurs. We describe the new giant ichthyosaur Cymbospondylus youngorum sp. nov. with a 2-meter-long skull from the Middle Triassic Fossil Hill Fauna of Nevada, USA, underscoring rapid size evolution despite the absence of many modern primary producers. Surprisingly, the Fossil Hill Fauna rivaled the composition of modern marine mammal faunas in terms of size range, and energy-flux models suggest that Middle Triassic marine food webs were able to support several large-bodied ichthyosaurs at high trophic levels, shortly after ichthyosaur origins.

Phylogenetic signal and evolutionary correlates of urban tolerance in a widespread neotropical lizard clade.

Urbanization is intensifying worldwide, and while some species tolerate and even exploit urban environments, many others are excluded entirely from this new habitat. Understanding the factors that underlie tolerance of urbanization is thus of rapidly growing importance. Here, we examine urban tolerance across a diverse group of lizards: Caribbean members of the neotropical genus Anolis. Our analyses reveal that urban tolerance has strong phylogenetic signal, suggesting that closely related species tend to respond similarly to urban environments. We propose that this characteristic of urban tolerance in anoles may be used to forecast the possible responses of species to increasing urbanization. In addition, we identified several key ecological and morphological traits that tend to be associated with tolerance in Anolis. Specifically, species experiencing hot and dry conditions in their natural environment and those that maintain higher body temperatures tend to have greater tolerance of urban habitats. We also found that tolerance of urbanization is positively associated with toepad lamella number and negatively associated with ventral scale density and relative hindlimb length. The identification of factors that predispose a species to be more or less urban tolerant can provide a starting point for conservation and sustainable development in our increasingly urbanized world.

PCCA: a program for phylogenetic canonical correlation analysis.

UNLABELLED: PCCA (phylogenetic canonical correlation analysis) is a new program for canonical correlation analysis of multivariate, continuously valued data from biological species. Canonical correlation analysis is a technique in which derived variables are obtained from two sets of original variables whereby the correlations between corresponding derived variables are maximized. It is a very useful multivariate statistical method for the calculation and analysis of correlations between character sets. The program controls for species non-independence due to phylogenetic history and computes canonical coefficients, correlations and scores; and conducts hypothesis tests on the canonical correlations. It can also compute a multivariate version of Pagel's lambda, which can then be used in the phylogenetic transformation. AVAILABILITY: PCCA is distributed as DOS/Windows, Mac OS X and Linux/Unix executables with a detailed program manual and is freely available on the World Wide Web at: http://anolis.oeb.harvard.edu/~liam/programs/.