The radiation continuum and the evolution of frog diversity.

Most of life's vast diversity of species and phenotypes is often attributed to adaptive radiation. Yet its contribution to species and phenotypic diversity of a major group has not been examined. Two key questions remain unresolved. First, what proportion of clades show macroevolutionary dynamics similar to adaptive radiations? Second, what proportion of overall species richness and phenotypic diversity do these adaptive-radiation-like clades contain? We address these questions with phylogenetic and morphological data for 1226 frog species across 43 families (which represent >99% of all species). Less than half of frog families resembled adaptive radiations (with rapid diversification and morphological evolution). Yet, these adaptive-radiation-like clades encompassed ~75% of both morphological and species diversity, despite rapid rates in other clades (e.g., non-adaptive radiations). Overall, we support the importance of adaptive-radiation-like evolution for explaining diversity patterns and provide a framework for characterizing macroevolutionary dynamics and diversity patterns in other groups.

Bridging Performance and Adaptive Landscapes to Understand Long-Term Functional Evolution.

AbstractUnderstanding functional adaptation demands an integrative framework that captures the complex interactions between form, function, ecology, and evolutionary processes. In this review, we discuss how to integrate the following two distinct approaches to better understand functional evolution: (1) the adaptive landscape approach (ALA), aimed at finding adaptive peaks for different ecologies, and (2) the performance landscape approach (PLA), aimed at finding performance peaks for different ecologies. We focus on the Ornstein-Uhlenbeck process as the evolutionary model for the ALA and on biomechanical modeling to estimate performance for the PLA. Whereas both the ALA and the PLA have each given insight into functional adaptation, separately they cannot address how much performance contributes to fitness or whether evolutionary constraints have played a role in form-function evolution. We show that merging these approaches leads to a deeper understanding of these issues. By comparing the locations of performance and adaptive peaks, we can infer how much performance contributes to fitness in species' current environments. By testing for the relevance of history on phenotypic variation, we can infer the influence of past selection and constraints on functional adaptation. We apply this merged framework in a case study of turtle shell evolution and explain how to interpret different possible outcomes. Even though such outcomes can be quite complex, they represent the multifaceted relations among function, fitness, and constraints.

A Cautionary Note on "A Cautionary Note on the Use of Ornstein Uhlenbeck Models in Macroevolutionary Studies".

Models based on the Ornstein-Uhlenbeck process have become standard for the comparative study of adaptation. Cooper et al. (2016) have cast doubt on this practice by claiming statistical problems with fitting Ornstein-Uhlenbeck models to comparative data. Specifically, they claim that statistical tests of Brownian motion may have too high Type I error rates and that such error rates are exacerbated by measurement error. In this note, we argue that these results have little relevance to the estimation of adaptation with Ornstein-Uhlenbeck models for three reasons. First, we point out that Cooper et al. (2016) did not consider the detection of distinct optima (e.g. for different environments), and therefore did not evaluate the standard test for adaptation. Second, we show that consideration of parameter estimates, and not just statistical significance, will usually lead to correct inferences about evolutionary dynamics. Third, we show that bias due to measurement error can be corrected for by standard methods. We conclude that Cooper et al. (2016) have not identified any statistical problems specific to Ornstein-Uhlenbeck models, and that their cautions against their use in comparative analyses are unfounded and misleading. [adaptation, Ornstein-Uhlenbeck model, phylogenetic comparative method.].

Ecology, sexual dimorphism, and jumping evolution in anurans.

Sexual dimorphism (SD) is a common feature of animals, and selection for sexually dimorphic traits may affect both functional morphological traits and organismal performance. Trait evolution through natural selection can also vary across environments. However, whether the evolution of organismal performance is distinct between the sexes is rarely tested in a phylogenetic comparative context. Anurans commonly exhibit sexual size dimorphism, which may affect jumping performance given the effects of body size on locomotion. They also live in a wide variety of microhabitats. Yet the relationships among dimorphism, performance, and ecology remain underexamined in anurans. Here, we explore relationships between microhabitat use, body size, and jumping performance in males and females to determine the drivers of dimorphic patterns in jumping performance. Using methods for predicting jumping performance through anatomical measurements, we describe how fecundity selection and natural selection associated with body size and microhabitat have likely shaped female jumping performance. We found that the magnitude of sexual size dimorphism (where females are about 14% larger than males) was much lower than dimorphism in muscle volume, where females had 42% more muscle than males (after accounting for body size). Despite these sometimes-large averages, phylogenetic t-tests failed to show the statistical significance of SD for any variable, indicating sexually dimorphic species tend to be closely related. While SD of jumping performance did not vary among microhabitats, we found female jumping velocity and energy differed across microhabitats. Overall, our findings indicate that differences in sex-specific reproductive roles, size, jumping-related morphology, and performance are all important determinants in how selection has led to the incredible ecophenotypic diversity of anurans.

Redefining Possible: Combining Phylogenomic and Supersparse Data in Frogs.

The data available for reconstructing molecular phylogenies have become wildly disparate. Phylogenomic studies can generate data for thousands of genetic markers for dozens of species, but for hundreds of other taxa, data may be available from only a few genes. Can these two types of data be integrated to combine the advantages of both, addressing the relationships of hundreds of species with thousands of genes? Here, we show that this is possible, using data from frogs. We generated a phylogenomic data set for 138 ingroup species and 3,784 nuclear markers (ultraconserved elements [UCEs]), including new UCE data from 70 species. We also assembled a supermatrix data set, including data from 97% of frog genera (441 total), with 1-307 genes per taxon. We then produced a combined phylogenomic-supermatrix data set (a "gigamatrix") containing 441 ingroup taxa and 4,091 markers but with 86% missing data overall. Likelihood analysis of the gigamatrix yielded a generally well-supported tree among families, largely consistent with trees from the phylogenomic data alone. All terminal taxa were placed in the expected families, even though 42.5% of these taxa each had >99.5% missing data and 70.2% had >90% missing data. Our results show that missing data need not be an impediment to successfully combining very large phylogenomic and supermatrix data sets, and they open the door to new studies that simultaneously maximize sampling of genes and taxa.

Phylogenetic analysis of adaptation in comparative physiology and biomechanics: overview and a case study of thermal physiology in treefrogs.

Comparative phylogenetic studies of adaptation are uncommon in biomechanics and physiology. Such studies require data collection from many species, a challenge when this is experimentally intensive. Moreover, researchers struggle to employ the most biologically appropriate phylogenetic tools for identifying adaptive evolution. Here, we detail an established but greatly underutilized phylogenetic comparative framework - the Ornstein-Uhlenbeck process - that explicitly models long-term adaptation. We discuss challenges in implementing and interpreting the model, and we outline potential solutions. We demonstrate use of the model through studying the evolution of thermal physiology in treefrogs. Frogs of the family Hylidae have twice colonized the temperate zone from the tropics, and such colonization likely involved a fundamental change in physiology due to colder and more seasonal temperatures. However, which traits changed to allow colonization is unclear. We measured cold tolerance and characterized thermal performance curves in jumping for 12 species of treefrogs distributed from the Neotropics to temperate North America. We then conducted phylogenetic comparative analyses to examine how tolerances and performance curves evolved and to test whether that evolution was adaptive. We found that tolerance to low temperatures increased with the transition to the temperate zone. In contrast, jumping well at colder temperatures was unrelated to biogeography and thus did not adapt during dispersal. Overall, our study shows how comparative phylogenetic methods can be leveraged in biomechanics and physiology to test the evolutionary drivers of variation among species.

Improving inference and avoiding overinterpretation of hidden-state diversification models: Specialized plant breeding has no effect on diversification in frogs.

The hidden-state speciation and extinction (HiSSE) model helps avoid spurious results when testing whether a character affects diversification rates. However, care must be taken to optimally analyze models and interpret results. Recently, Tonini et al. (TEA hereafter) studied anuran (frog and toad) diversification with HiSSE methods. They concluded that their focal state, breeding in phytotelmata, increases net diversification rates. Yet this conclusion is counterintuitive, because the state that purportedly increases net diversification rates is 14 times rarer among species than the alternative. Herein, I revisit TEA's analyses and demonstrate problems with inferring model likelihoods, conducting post hoc tests, and interpreting results. I also reevaluate their top models and find that diverse strategies are necessary to reach the parameter values that maximize each model's likelihood. In contrast to TEA, I find no support for an effect of phytotelm breeding on net diversification rates in Neotropical anurans. In particular, even though the most highly supported models include the focal character, averaging parameter estimates over hidden states shows that the focal character does not influence diversification rates. Finally, I suggest ways to better analyze and interpret complex diversification models-both state-dependent and beyond-for future studies in other organisms.

Testing for adaptive radiation: A new approach applied to Madagascar frogs.

Adaptive radiation is a key topic at the intersection of ecology and evolutionary biology. Yet the definition and identification of adaptive radiation both remain contentious. Here, we introduce a new approach for identifying adaptive radiations that combines key aspects of two widely used definitions. Our approach compares evolutionary rates in morphology, performance, and diversification between the candidate radiation and other clades. We then apply this approach to a putative adaptive radiation of frogs from Madagascar (Mantellidae). We present new data on morphology and performance from mantellid frogs, then compare rates of diversification and multivariate evolution of size, shape, and performance between mantellids and other frogs. We find that mantellids potentially pass our test for accelerated rates of evolution for shape, but not for size, performance, or diversification. Our results demonstrate that clades can have accelerated phenotypic evolution without rapid diversification (dubbed "adaptive non-radiation"). We also highlight general issues in testing for adaptive radiation, including taxon sampling and the problem of including another adaptive radiation among the comparison clades. Finally, we suggest that similar tests should be conducted on other putative adaptive radiations on Madagascar, comparing their evolutionary rates to those of related clades outside Madagascar. Based on our results, we speculate that older Madagascar clades may show evolutionary patterns more similar to those on a continent than an island.

Disturbed habitats locally reduce the signal of deep evolutionary history in functional traits of plants.

The functioning of present ecosystems reflects deep evolutionary history of locally cooccurring species if their functional traits show high phylogenetic signal (PS). However, we do not understand what drives local PS. We hypothesize that local PS is high in undisturbed and stressful habitats, either due to ongoing local assembly of species that maintained ancestral traits, or to past evolutionary maintenance of ancestral traits within habitat species-pools, or to both. We quantified PS and diversity of 10 traits within 6704 local plant communities across 38 Dutch habitat types differing in disturbance or stress. Mean local PS varied 50-fold among habitat types, often independently of phylogenetic or trait diversity. Mean local PS decreased with disturbance but showed no consistent relationship to stress. Mean local PS exceeded species-pool PS, reflecting nonrandom subsampling from the pool. Disturbance or stress related more strongly to mean local than to species-pool PS. Disturbed habitats harbour species with evolutionary divergent trait values, probably driven by ongoing, local assembly of species: environmental fluctuations might maintain different trait values within lineages through an evolutionary storage effect. If functional traits do not reflect phylogeny, ecosystem functioning might not be contingent on the presence of particular lineages, and lineages might establish evolutionarily novel interactions.

Interspecific variation in bristle number on forewings of tiny insects does not influence clap-and-fling aerodynamics.

Miniature insects must overcome significant viscous resistance in order to fly. They typically possess wings with long bristles on the fringes and use a clap-and-fling mechanism to augment lift. These unique solutions to the extreme conditions of flight at tiny sizes (<2 mm body length) suggest that natural selection has optimized wing design for better aerodynamic performance. However, species vary in wingspan, number of bristles (n) and bristle gap (G) to diameter (D) ratio (G/D). How this variation relates to body length (BL) and its effects on aerodynamics remain unknown. We measured forewing images of 38 species of thrips and 21 species of fairyflies. Our phylogenetic comparative analyses showed that n and wingspan scaled positively and similarly with BL across both groups, whereas G/D decreased with BL, with a sharper decline in thrips. We next measured aerodynamic forces and visualized flow on physical models of bristled wings performing clap-and-fling kinematics at a chord-based Reynolds number of 10 using a dynamically scaled robotic platform. We examined the effects of dimensional (G, D, wingspan) and non-dimensional (n, G/D) geometric variables on dimensionless lift and drag. We found that: (1) increasing G reduced drag more than decreasing D; (2) changing n had minimal impact on lift generation; and (3) varying G/D minimally affected aerodynamic forces. These aerodynamic results suggest little pressure to functionally optimize n and G/D. Combined with the scaling relationships between wing variables and BL, much wing variation in tiny flying insects might be best explained by underlying shared growth factors.

What Determines the Distinct Morphology of Species with a Particular Ecology? The Roles of Many-to-One Mapping and Trade-Offs in the Evolution of Frog Ecomorphology and Performance.

Organisms inhabiting a specific environment often have distinct morphology, but the factors that affect this fit are unclear when multiple morphological traits affect performance in multiple behaviors. Does the realized morphology of a species reflect a compromise in performance among behaviors (i.e., trade-offs)? Or does many-to-one mapping result in morphological distinctness without compromising performance across behaviors? The importance of these principles in organismal design has rarely been compared at the macroevolutionary scale. Here I study 191 species of frogs from around the world that inhabit different microhabitats, using models of phenotypic evolution to examine how form-function relationships may explain the fit between ecology and morphology. I found three key results. First, despite being distinct in leg morphology, ecomorphs were similar in jumping performance. Second, ecomorphs that regularly swim showed higher swimming performance, which paralleled the higher leg muscle mass in these taxa. Third, many-to-one mapping of form onto function occurred at all but the highest levels of both jumping and swimming performance. The seemingly contradictory first two results were explained by the third: when one behavior occurs in all species while another is restricted to a subset, many-to-one mapping allows species with distinct ecologies to have distinct body forms that reflect their specialized behavior while maintaining similar performance in a more general shared behavior.

The relationship between habitat use and body shape in geckos.

Geckos are a highly diverse group of lizards, with more than 1,700 species that exhibit a wide range of behaviors, ecologies, and sizes. However, no study has examined links between habitat use and body shape in pad-bearing geckos. We set out to answer a basic question using a data set of pad-bearing geckos (112 species, 103 pad-bearing, 9 padless, 42 genera): Do geckos that occur in different habitats also differ in body shape? Overall, we found that body shape was surprisingly similar among our sample of pad-bearing species, with the exception of the genus Uroplatus, which was clearly distinct from other geckos due to its depressed body and long limbs. However, the padless geckos differed in body shape from the pad-bearing geckos by having longer arms and legs and less rotund bodies. We found that about half of the pad-bearing species primarily inhabit trees, with the other half, divided approximately equally among rocks, the ground, and mixed habitats. We found no significant links between habitat use and body shape, nor any propensity for larger species to occupy different habitats than smaller species. Padless species tend to inhabit rock and ground substrates. Our results indicate that pad-bearing geckos have a relatively uniform body form, which contrasts with to their diversity in color, size, and behavior. Indeed, our data show that the general gecko body shape is suitable for a wide range of habitats, ranging from arboreal to terrestrial. This pattern is a departure from other ecomorphological studies and suggests that geckos may not easily fit into the mold of adaptive radiation, as suggested by prior studies.

Diversification rates are more strongly related to microhabitat than climate in squamate reptiles (lizards and snakes).

Patterns of species richness among clades can be directly explained by the ages of clades or their rates of diversification. The factors that most strongly influence diversification rates remain highly uncertain, since most studies typically consider only a single predictor variable. Here, we explore the relative impacts of macroclimate (i.e., occurring in tropical vs. temperate regions) and microhabitat use (i.e., terrestrial, fossorial, arboreal, aquatic) on diversification rates of squamate reptile clades (lizards and snakes). We obtained data on microhabitat, macroclimatic distribution, and phylogeny for >4000 species. We estimated diversification rates of squamate clades (mostly families) from a time-calibrated tree, and used phylogenetic methods to test relationships between diversification rates and microhabitat and macroclimate. Across 72 squamate clades, the best-fitting model included microhabitat but not climatic distribution. Microhabitat explained ∼37% of the variation in diversification rates among clades, with a generally positive impact of arboreal microhabitat use on diversification, and negative impacts of fossorial and aquatic microhabitat use. Overall, our results show that the impacts of microhabitat on diversification rates can be more important than those of climate, despite much greater emphasis on climate in previous studies.

Microhabitat and Climatic Niche Change Explain Patterns of Diversification among Frog Families.

A major goal of ecology and evolutionary biology is to explain patterns of species richness among clades. Differences in rates of net diversification (speciation minus extinction over time) may often explain these patterns, but the factors that drive variation in diversification rates remain uncertain. Three important candidates are climatic niche position (e.g., whether clades are primarily temperate or tropical), rates of climatic niche change among species within clades, and microhabitat (e.g., aquatic, terrestrial, arboreal). The first two factors have been tested separately in several studies, but the relative importance of all three is largely unknown. Here we explore the correlates of diversification among families of frogs, which collectively represent ∼88% of amphibian species. We assemble and analyze data on phylogeny, climate, and microhabitat for thousands of species. We find that the best-fitting phylogenetic multiple regression model includes all three types of variables: microhabitat, rates of climatic niche change, and climatic niche position. This model explains 67% of the variation in diversification rates among frog families, with arboreal microhabitat explaining ∼31%, niche rates ∼25%, and climatic niche position ∼11%. Surprisingly, we show that microhabitat can have a much stronger influence on diversification than climatic niche position or rates of climatic niche change.

Testing Convergence Versus History: Convergence Dominates Phenotypic Evolution for over 150 Million Years in Frogs.

Striking evolutionary convergence can lead to similar sets of species in different locations, such as in cichlid fishes and Anolis lizards, and suggests that evolution can be repeatable and predictable across clades. Yet, most examples of convergence involve relatively small temporal and/or spatial scales. Some authors have speculated that at larger scales (e.g., across continents), differing evolutionary histories will prevent convergence. However, few studies have compared the contrasting roles of convergence and history, and none have done so at large scales. Here we develop a two-part approach to test the scale over which convergence can occur, comparing the relative importance of convergence and history in macroevolution using phylogenetic models of adaptive evolution. We apply this approach to data from morphology, ecology, and phylogeny from 167 species of anuran amphibians (frogs) from 10 local sites across the world, spanning ~160 myr of evolution. Mapping ecology on the phylogeny revealed that similar microhabitat specialists (e.g., aquatic, arboreal) have evolved repeatedly across clades and regions, producing many evolutionary replicates for testing for morphological convergence. By comparing morphological optima for clades and microhabitat types (our first test), we find that convergence associated with microhabitat use dominates frog morphological evolution, producing recurrent ecomorphs that together encompass all sampled species in each community in each region. However, our second test, which examines whether and how much species differ from their inferred optima, shows that convergence is incomplete: that is, phenotypes of most species are still somewhat distant from the estimated optimum for each microhabitat, seemingly because of insufficient time for more complete adaptation (an effect of history). Yet, these effects of history are related to past ecologies, and not clade membership. Overall, our study elucidates the dominant drivers of morphological evolution across a major vertebrate clade and shows that evolution can be repeatable at much greater temporal and spatial scales than commonly thought. It also provides an analytical framework for testing other potential examples of large-scale convergence.

Evolutionary conservatism and convergence both lead to striking similarity in ecology, morphology and performance across continents in frogs.

Many clades contain ecologically and phenotypically similar species across continents, yet the processes generating this similarity are largely unstudied, leaving fundamental questions unanswered. Is similarity in morphology and performance across assemblages caused by evolutionary convergence or by biogeographic dispersal of evolutionarily conserved ecotypes? Does convergence to new ecological conditions erase evidence of past adaptation? Here, we analyse ecology, morphology and performance in frog assemblages from three continents (Asia, Australia and South America), assessing the importance of dispersal and convergent evolution in explaining similarity across regions. We find three striking results. First, species using the same microhabitat type are highly similar in morphology and performance across both clades and continents. Second, some species on different continents owe their similarity to dispersal and evolutionary conservatism (rather than evolutionary convergence), even over vast temporal and spatial scales. Third, in one case, an ecologically specialized ancestor radiated into diverse ecotypes that have converged with those on other continents, largely erasing traces of past adaptation to their ancestral ecology. Overall, our study highlights the roles of both evolutionary conservatism and convergence in explaining similarity in species traits over large spatial and temporal scales and demonstrates a statistical framework for addressing these questions in other systems.

Phylogenetic origins of local-scale diversity patterns and the causes of Amazonian megadiversity.

What explains the striking variation in local species richness across the globe and the remarkable diversity of rainforest sites in Amazonia? Here, we apply a novel phylogenetic approach to these questions, using treefrogs (Hylidae) as a model system. Hylids show dramatic variation in local richness globally and incredible local diversity in Amazonia. We find that variation in local richness is not explained primarily by climatic factors, rates of diversification (speciation and extinction) nor morphological variation. Instead, local richness patterns are explained predominantly by the timing of colonization of each region, and Amazonian megadiversity is linked to the long-term sympatry of multiple clades in that region. Our results also suggest intriguing interactions between clade diversification, trait evolution and the accumulation of local richness. Specifically, sympatry between clades seems to slow diversification and trait evolution, but prevents neither the accumulation of local richness over time nor the co-occurrence of similar species.

An expanded phylogeny of treefrogs (Hylidae) based on nuclear and mitochondrial sequence data.

The treefrogs (Hylidae) make up one of the most species-rich families of amphibians. With 885 species currently described, they contain >13% of all amphibian species. In recent years, there has been considerable progress in resolving hylid phylogeny. However, the most comprehensive phylogeny to date (Wiens et al., 2006) included only 292 species, was based only on parsimony, provided only poor support for most higher-level relationships, and conflicted with previous hypotheses in several parts (including the monophyly and relationships of major clades of Hylinae). Here, we present an expanded phylogeny for hylid frogs, including data for 362 hylid taxa for up to 11 genes (4 mitochondrial, 7 nuclear), including 70 additional taxa and >270 sequences not included in the previously most comprehensive analysis. The new tree from maximum likelihood analysis is more well-resolved, strongly supported, and concordant with previous hypotheses, and provides a framework for future systematic, biogeographic, ecological, and evolutionary studies.

Community assembly through evolutionary diversification and dispersal in Middle American treefrogs.

How are ecologically diverse organisms added to local assemblages to create the community structure we see today? In general, within a given region or community, a given trait (character state) may either evolve in situ or be added through dispersal after having evolved elsewhere. Here, we develop simple metrics to quantify the relative importance of these processes and then apply them to a case study in Middle American treefrogs. We examined two ecologically important characters (larval habitat and body size) among 39 communities, using phylogenetic and ecological information from 278 species both inside and outside the region. For each character, variation among communities reflects complex patterns of evolution and dispersal. Our results support several general hypotheses about community assembly, which may apply to many other systems: (1) elevation can play an important role in creating patterns of community structure within a region, (2) contrary to expectations, species can invade communities in which species with similar ecological traits are already present, (3) dispersal events tend to occur between areas with similar climatic regimes, and (4) the first lineage to invade a region diversifies the most ecologically, whereas later invasions show limited change.

Phylogenetic evidence for competitively driven divergence: body-size evolution in Caribbean treefrogs (Hylidae: Osteopilus).

Understanding the role of competition in explaining phenotypic diversity is a challenging problem, given that the most divergent species may no longer compete today. However, convergent evolution of extreme body sizes across communities may offer evidence of past competition. For example, many treefrog assemblages around the world have convergently evolved species with very large and small body sizes. To better understand this global pattern, we studied body-size diversification within the small, endemic radiation of Caribbean treefrogs (Osteopilus). We introduce a suite of analyses designed to help reveal the signature of past competition. Diet analyses show that Osteopilus are generalist predators and that prey size is strongly associated with body size, suggesting that body-size divergence facilitates resource partitioning. Community assembly models indicate that treefrog body-size distributions in Jamaica and Hispaniola are consistent with expectations from competition. Phylogenetic analyses show that similar body-size extremes in Jamaica and Hispaniola have originated through parallel evolution on each island, and the rate of body-size evolution in Osteopilus is accelerated relative to mainland treefrogs. Together, these results suggest that competition may have driven the rapid diversification of body sizes in Caribbean treefrogs to the extremes seen in treefrog communities around the world.