The influence of multiple functional demands on morphological diversification: A test on turtle shells

Organismal parts are often involved in the performance of more than one function. The role of trade‐offs in influencing phenotypic evolution of such parts is well‐studied; less well‐understood is their role in influencing phenotypic diversity. Increases in the number of functions a part is involved in may inhibit subsequent diversification, as the number of trade‐offs increases. Alternately, such an increase might promote phenotypic diversification, by increasing adaptive landscape complexity and promoting specialization for different roles. We compare these predictions by testing whether aquatic turtle shells, which resist loads, act as hydrodynamic elements, facilitate self‐righting, and exchange heat with the environment, differ in phenotypic diversity from those of terrestrial species, which perform all the same functions except for hydrodynamics. We used 53 3D landmarks digitized on 2722 specimens of 274 hard‐shelled turtle species to quantify shell shape variation, and a set of phylogenetic hypotheses to examine evolutionary patterns. Terrestrial turtles consistently had higher phenotypic diversity than aquatic species. Differences are not due to differences in the rates of evolution between the two groups, but rather differences in evolutionary mode. Thus this study supports the traditional view of the role of multiple functions in determining phenotypic diversity.     

Functional performance of turtle humerus shape across an ecological adaptive landscape

The concept of the adaptive landscape has been invaluable to evolutionary biologists for visualizing the dynamics of selection and adaptation, and is increasingly being used to study morpho‐functional data. Here, we construct adaptive landscapes to explore functional trade‐offs associated with variation in humerus morphology among turtles adapted to three different locomotor environments: marine, semiaquatic, and terrestrial. Humerus shape from 40 species of cryptodire turtles was quantified using a pseudolandmark approach. Hypothetical shapes were extracted in a grid across morphospace and four functional traits (strength, stride length, mechanical advantage, and hydrodynamics) measured on those shapes. Quantitative trait modeling was used to construct adaptive landscapes that optimize the functional traits for each of the three locomotor ecologies. Our data show that turtles living in different environments have statistically different humeral shapes. The optimum adaptive landscape for each ecology is defined by a different combination of performance trade‐offs, with turtle species clustering around their respective adaptive peak. Further, species adhere to pareto fronts between marine–semiaquatic and semiaquatic–terrestrial optima, but not between marine–terrestrial. Our study demonstrates the utility of adaptive landscapes in informing the link between form, function, and ecological adaptation, and establishes a framework for reconstructing turtle ecological evolution using isolated humeri from the fossil record.     

Evolution of morphology and locomotor performance in anurans: relationships with microhabitat diversification

The relationships between morphology, performance, behavior and ecology provide evidence for multiple and complex phenotypic adaptations. The anuran body plan, for example, is evolutionarily conserved and shows clear specializations to jumping performance back at least to the early Jurassic. However, there are instances of more recent adaptation to habit diversity in the post‐cranial skeleton, including relative limb length. The present study tested adaptive models of morphological evolution in anurans associated with the diversity of microhabitat use (semi‐aquatic arboreal, fossorial, torrent, and terrestrial) in species of anuran amphibians from Brazil and Australia. We use phylogenetic comparative methods to determine which evolutionary models, including Brownian motion (BM) and Ornstein‐Uhlenbeck (OU) are consistent with morphological variation observed across anuran species. Furthermore, this study investigated the relationship of maximum distance jumped as a function of components of morphological variables and microhabitat use. We found there are multiple optima of limb lengths associated to different microhabitats with a trend of increasing hindlimbs in torrent, arboreal, semi‐aquatic whereas fossorial and terrestrial species evolve toward optima with shorter hindlimbs. Moreover, arboreal, semi‐aquatic and torrent anurans have higher jumping performance and longer hindlimbs, when compared to terrestrial and fossorial species. We corroborate the hypothesis that evolutionary modifications of overall limb morphology have been important in the diversification of locomotor performance along the anuran phylogeny. Such evolutionary changes converged in different phylogenetic groups adapted to similar microhabitat use in two different zoogeographical regions.     

Biomechanics on the half shell: functional performance influences patterns of morphological variation in the emydid turtle carapace

This study uses the carapace of emydid turtles to address hypothesized differences between terrestrial and aquatic species. Geometric morphometrics are used to quantify shell shape, and performance is estimated for two shell functions: shell strength and hydrodynamics. Aquatic turtle shells differ in shape from terrestrial turtle shells and are characterized by lower frontal areas and presumably lower drag. Terrestrial turtle shells are stronger than those of aquatic turtles; many-to-one mapping of morphology to function does not entirely mitigate a functional trade-off between mechanical strength and hydrodynamic performance. Furthermore, areas of morphospace characterized by exceptionally poor performance in either of the functions are not occupied by any emydid species. Though aquatic and terrestrial species show no significant differences in the rate of morphological evolution, aquatic species show a higher lineage density, indicative of a greater amount of convergence in their evolutionary history. The techniques employed in this study, including the modeling of theoretical shapes to assess performance in unoccupied areas of morphospace, suggest a framework for future studies of morphological variation.     

Evidence for a multi‐peak adaptive landscape in the evolution of trophic morphology in terapontid fishes

Using the Australian marine‐freshwater terapontid fishes as a model system, we examined the role of dietary phenotypic optima in an adaptive macro‐evolutionary landscape. Comparative modelling relying on both a priori and data‐driven identification of selective regimes suggested multi‐peak models as best describing much of the dietary phenotypic landscape of terapontids. Both approaches identified common phenotypic optima for different lineages of marine and freshwater herbivores, and minimal differentiation between carnivores and omnivores, irrespective of their phylogenetic relationships, as the model best describing morphological evolution. Significant correlations also existed between these phenotypic axes and proportions of non‐animal dietary items in species’ diets. While simulation results provided evidence for a multi‐peak adaptive landscape in the evolution of trophic morphology in terapontids, they could not rule out chance convergence in these adaptive peaks. However, they do provide scope for identifying areas for more detailed, functionally specific study of phenotypic convergence in herbivorous terapontid trophic habits. © 2014 The Linnean Society of London, Biological Journal of the Linnean Society, 2014, 113 , 623–634.     

Cope's rule and the adaptive landscape of dinosaur body size evolution

The largest known dinosaurs weighed at least 20 million times as much as the smallest, indicating exceptional phenotypic divergence. Previous studies have focused on extreme giant sizes, tests of Cope's rule, and miniaturization on the line leading to birds. We use non‐uniform macroevolutionary models based on Ornstein–Uhlenbeck and trend processes to unify these observations, asking: what patterns of evolutionary rates, directionality and constraint explain the diversification of dinosaur body mass? We find that dinosaur evolution is constrained by attraction to discrete body size optima that undergo rare, but abrupt, evolutionary shifts. This model explains both the rarity of multi‐lineage directional trends, and the occurrence of abrupt directional excursions during the origins of groups such as tiny pygostylian birds and giant sauropods. Most expansion of trait space results from rare, constraint‐breaking innovations in just a small number of lineages. These lineages shifted rapidly into novel regions of trait space, occasionally to small sizes, but most often to large or giant sizes. As with Cenozoic mammals, intermediate body sizes were typically attained only transiently by lineages on a trajectory from small to large size. This demonstrates that bimodality in the macroevolutionary adaptive landscape for land vertebrates has existed for more than 200 million years.     

Time and environment explain the current richness distribution of non‐marine turtles worldwide

Ecological, historical, and evolutionary hypotheses are important to explain geographical diversity gradients in many clades, but few studies have combined them into a single analysis allowing a comparison of their relative importance. This study aimed to evaluate the relative importance of ecological, historical, and evolutionary hypotheses in explaining the current global distribution of non‐marine turtles, a group whose distribution patterns are still poorly explored. We used data from distribution range maps of 336 species of non‐marine turtles, environmental layers, and phylogeny to obtain richness estimates of these animals in 2° × 2° cells and predictors related to ecological, evolutionary and historical hypotheses driving richness patterns. Then we used a path analysis to evaluate direct and indirect effects of the predictors on turtle richness. Ancestral area reconstruction was also performed in order to evaluate the influence of time‐for‐speciation in the current diversity of the group. We found that environmental variables had the highest direct effects on non‐marine turtle richness, whereas diversification rates and area available in the last 55 million yr minimally influenced turtle distributions. We found evidence for the time‐for‐speciation effect, since regions colonized early were generally richer than recently colonized regions. In addition, regions with a high number of colonization events had a higher number of turtle species. Our results suggested that ecological processes may influence non‐marine turtle richness independent of diversification rates, but they are probably related to dispersal abilities. However, colonization time was also an important component that must be taken into account. Finally, our study provided additional support for the importance of ecological (climate and productivity) and historical (time‐for‐speciation and dispersal) processes in shaping current biodiversity patterns.     

Ecological diversification associated with the benthic‐to‐pelagic transition by North American minnows

Ecological opportunity is often regarded as a key factor that explains why diversity is unevenly distributed across life. Colonization of novel environments or adaptive zones may promote diversification. North American minnows exhibit an ancestral benthic‐to‐pelagic habitat shift that coincided with a burst in diversification. Here, we evaluate the phenotypic and ecological implications of this habitat shift by assessing craniofacial and dietary traits among 34 species and testing for morphology–diet covariation, convergence and adaptive optima. There were several instances of morphology–diet covariation such as correlations between mouth angle and the consumption of terrestrial insects and between relative gut length and the consumption of algae. After accounting for size and phylogenetic nonindependence, benthic species had longer heads, longer snouts, eyes positioned higher on their head, smaller mouth angles and longer digestive tracts than pelagic minnows. Benthic minnows also consumed more algae but less terrestrial insects, by volume, than pelagic minnows. Lastly, there were three distinct evolutionary regimes and more convergence in morphology and dietary characteristics than expected under a Brownian motion model of evolution. These findings indicate that colonization of the pelagic zone by minnows involved myriad phenotypic and dietary changes associated with exploitation of terrestrial subsidies. Thus, minnows exhibit phenotype–dietary covariation, an expansion of ecological roles and a burst in diversification rates in response to the ecological opportunity afforded by the colonization of a novel habitat.     

Disruptive selection as a driver of evolutionary branching and caste evolution in social insects

Theory suggests that evolutionary branching via disruptive selection may be a relatively common and powerful force driving phenotypic divergence. Here, we extend this theory to social insects, which have novel social axes of phenotypic diversification. Our model, built around turtle ant (Cephalotes) biology, is used to explore whether disruptive selection can drive the evolutionary branching of divergent colony phenotypes that include a novel soldier caste. Soldier evolution is a recurrent theme in social insect diversification that is exemplified in the turtle ants. We show that phenotypic mutants can gain competitive advantages that induce disruptive selection and subsequent branching. A soldier caste does not generally appear before branching, but can evolve from subsequent competition. The soldier caste then evolves in association with specialized resource preferences that maximize defensive performance. Overall, our model indicates that resource specialization may occur in the absence of morphological specialization, but that when morphological specialization evolves, it is always in association with resource specialization. This evolutionary coupling of ecological and morphological specialization is consistent with recent empirical evidence, but contrary to predictions of classical caste theory. Our model provides a new theoretical understanding of the ecology of caste evolution that explicitly considers the process of adaptive phenotypic divergence and diversification.     

Diversification on an ecologically constrained adaptive landscape

We used phylogenetic analysis of body-size ecomorphs in a crustacean species complex to gain insight into how spatial complexity of ecological processes generates and maintains biological diversity. Studies of geographically widespread species of Hyalella amphipods show that phenotypic evolution is tightly constrained in a manner consistent with adaptive responses to alternative predation regimes. A molecular phylogeny indicates that evolution of Hyalella ecomorphs is characterized by parallel evolution and by phenotypic stasis despite substantial levels of underlying molecular change. The phylogeny suggests that species diversification sometimes occurs by niche shifts, and sometimes occurs without a change in niche. Moreover, diversification in the Hyalella ecomorphs has involved the repeated evolution of similar phenotypic forms that exist in similar ecological settings, a hallmark of adaptive evolution. The evolutionary stasis observed in clades separated by substantial genetic divergence, but existing in similar habitats, is also suggestive of stabilizing natural selection acting to constrain phenotypic evolution within narrow bounds. We interpret the observed decoupling of genetic and phenotypic diversification in terms of adaptive radiation on an ecologically constrained adaptive landscape, and suggest that ecological constraints, perhaps acting together with genetic and functional constraints, may explain the parallel evolution and evolutionary stasis inferred by the phylogeny.     

The stationary distribution of a continuously varying strategy in a class-structured population under mutation-selection-drift balance

Many traits and/or strategies expressed by organisms are quantitative phenotypes. Because populations are of finite size and genomes are subject to mutations, these continuously varying phenotypes are under the joint pressure of mutation, natural selection and random genetic drift. This article derives the stationary distribution for such a phenotype under a mutation-selection-drift balance in a class-structured population allowing for demographically varying class sizes and/or changing environmental conditions. The salient feature of the stationary distribution is that it can be entirely characterized in terms of the average size of the gene pool and Hamilton's inclusive fitness effect. The exploration of the phenotypic space varies exponentially with the cumulative inclusive fitness effect over state space, which determines an adaptive landscape. The peaks of the landscapes are those phenotypes that are candidate evolutionary stable strategies and can be determined by standard phenotypic selection gradient methods (e.g. evolutionary game theory, kin selection theory, adaptive dynamics). The curvature of the stationary distribution provides a measure of the stability by convergence of candidate evolutionary stable strategies, and it is evaluated explicitly for two biological scenarios: first, a coordination game, which illustrates that, for a multipeaked adaptive landscape, stochastically stable strategies can be singled out by letting the size of the gene pool grow large; second, a sex-allocation game for diploids and haplo-diploids, which suggests that the equilibrium sex ratio follows a Beta distribution with parameters depending on the features of the genetic system.     

Adaptive landscape and functional diversity of Neotropical cichlids: implications for the ecology and evolution of Cichlinae (Cichlidae; Cichliformes)

Morphological, lineage and ecological diversity can vary substantially even among closely related lineages. Factors that influence morphological diversification, especially in functionally relevant traits, can help to explain the modern distribution of disparity across phylogenies and communities. Multivariate axes of feeding functional morphology from 75 species of Neotropical cichlid and a stepwise‐AIC algorithm were used to estimate the adaptive landscape of functional morphospace in Cichlinae. Adaptive landscape complexity and convergence, as well as the functional diversity of Cichlinae, were compared with expectations under null evolutionary models. Neotropical cichlid feeding function varied primarily between traits associated with ram feeding vs. suction feeding/biting and secondarily with oral jaw muscle size and pharyngeal crushing capacity. The number of changes in selective regimes and the amount of convergence between lineages was higher than expected under a null model of evolution, but convergence was not higher than expected under a similarly complex adaptive landscape. Functional disparity was compatible with an adaptive landscape model, whereas the distribution of evolutionary change through morphospace corresponded with a process of evolution towards a single adaptive peak. The continentally distributed Neotropical cichlids have evolved relatively rapidly towards a number of adaptive peaks in functional trait space. Selection in Cichlinae functional morphospace is more complex than expected under null evolutionary models. The complexity of selective constraints in feeding morphology has likely been a significant contributor to the diversity of feeding ecology in this clade.     

Tightly congruent bursts of lineage and phenotypic diversification identified in a continental ant radiation

Adaptive diversification is thought to be shaped by ecological opportunity. A prediction of this ecological process of diversification is that it should result in congruent bursts of lineage and phenotypic diversification, but few studies have found this expected association. Here, we study the relationship between rates of lineage diversification and body size evolution in the turtle ants, a diverse Neotropical clade. Using a near complete, time‐calibrated phylogeny we investigated lineage diversification dynamics and body size disparity through model fitting analyses and estimation of per‐lineage rates of cladogenesis and phenotypic evolution. We identify an exceptionally high degree of congruence between the high rates of lineage and body size diversification in a young clade undergoing renewed diversification in the ecologically distinct Chacoan biogeographical region of South America. It is likely that the region presented turtle ants with novel ecological opportunity, which facilitated a nested burst of diversification and phenotypic evolution within the group. Our results provide a compelling quantitative example of tight congruence between rates of lineage and phenotypic diversification, meeting the key predicted pattern of adaptive diversification shaped by ecological opportunity.     

Possible human impacts on adaptive radiation: beak size bimodality in Darwin's finches

Adaptive radiation is facilitated by a rugged adaptive landscape, where fitness peaks correspond to trait values that enhance the use of distinct resources. Different species are thought to occupy the different peaks, with hybrids falling into low-fitness valleys between them. We hypothesize that human activities can smooth adaptive landscapes, increase hybrid fitness and hamper evolutionary diversification. We investigated this possibility by analysing beak size data for 1755 Geospiza fortis measured between 1964 and 2005 on the island of Santa Cruz, Galápagos. Some populations of this species can display a resource-based bimodality in beak size, which mirrors the greater beak size differences among species. We first show that an historically bimodal population at one site, Academy Bay, has lost this property in concert with a marked increase in local human population density. We next show that a nearby site with lower human impacts, El Garrapatero, currently manifests strong bimodality. This comparison suggests that bimodality can persist when human densities are low (Academy Bay in the past, El Garrapatero in the present), but not when they are high (Academy Bay in the present). Human activities may negatively impact diversification in 'young' adaptive radiations, perhaps by altering adaptive landscapes.     

Finite element modeling of shell shape in the freshwater turtle Pseudemys concinna reveals a trade‐off between mechanical strength and hydrodynamic efficiency

Aquatic species can experience different selective pressures on morphology in different flow regimes. Species inhabiting lotic regimes often adapt to these conditions by evolving low‐drag (i.e., streamlined) morphologies that reduce the likelihood of dislodgment or displacement. However, hydrodynamic factors are not the only selective pressures influencing organismal morphology and shapes well suited to flow conditions may compromise performance in other roles. We investigated the possibility of morphological trade‐offs in the turtle Pseudemys concinna. Individuals living in lotic environments have flatter, more streamlined shells than those living in lentic environments; however, this flatter shape may also make the shells less capable of resisting predator‐induced loads. We tested the idea that “lotic” shell shapes are weaker than “lentic” shell shapes, concomitantly examining effects of sex. Geometric morphometric data were used to transform an existing finite element shell model into a series of models corresponding to the shapes of individual turtles. Models were assigned identical material properties and loaded under identical conditions, and the stresses produced by a series of eight loads were extracted to describe the strength of the shells. “Lotic” shell shapes produced significantly higher stresses than “lentic” shell shapes, indicating that the former is weaker than the latter. Females had significantly stronger shell shapes than males, although these differences were less consistent than differences between flow regimes. We conclude that, despite the potential for many‐to‐one mapping of shell shape onto strength, P. concinna experiences a trade‐off in shell shape between hydrodynamic and mechanical performance. This trade‐off may be evident in many other turtle species or any other aquatic species that also depend on a shell for defense. However, evolution of body size may provide an avenue of escape from this trade‐off in some cases, as changes in size can drastically affect mechanical performance while having little effect on hydrodynamic performance. J. Morphol. 2011. © 2011 Wiley‐Liss, Inc.     

Morphological convergence of shell shape in distantly related scallop species (Mollusca: Pectinidae)

Morphological convergence is a central concept in evolutionary biology, but convergent patterns remain under‐studied in nonvertebrate organisms. Some scallop species exhibit long‐distance swimming, a behaviour whose biomechanical requirements probably generate similar selective regimes. We tested the hypothesis that shell shape similarity in long‐distance swimming species is a result of convergent evolution. Using landmark‐based geometric morphometrics, we quantified shell shape in seven species representing major behavioural habits. All species displayed distinct shell shapes, with the exception of the two long‐distance swimmers, whose shells were indistinguishable. These species also displayed reduced morphological variance relative to other taxa. Finally, a phylogenetic simulation revealed that these species were more similar in their shell shape than was expected under Brownian motion, the model of character evolution that best described changes in shell shape. Together, these findings reveal that convergent evolution of shell shape occurs in scallops, and suggest that selection for shell shape and behaviour may be important in the diversification of the group. © 2011 The Linnean Society of London, Zoological Journal of the Linnean Society, 2011, 163 , 571–584.     

SELECTION FOR MECHANICAL ADVANTAGE UNDERLIES MULTIPLE CRANIAL OPTIMA IN NEW WORLD LEAF‐NOSED BATS

Selection for divergent performance optima has been proposed as a central mechanism underlying adaptive radiation. Uncovering multiple optima requires identifying forms associated with different adaptive zones and linking those forms to performance. However, testing and modeling the performance of complex morphologies like the cranium is challenging. We introduce a three‐dimensional finite‐element (FE) model of the cranium that can be morphed into different shapes by varying simple parameters to investigate the relationship between two engineering‐based measures of performance, mechanical advantage and von Mises stress, and four divergent adaptive zones occupied by New World Leaf‐nosed bats. To investigate these relationships, we tested the fit of Brownian motion and Ornstein–Uhlenbeck models of evolution in mechanical advantage and von Mises stress using dated multilocus phylogenies. The analyses revealed three performance optima for mechanical advantage among species from three adaptive zones: bats that eat nectar; generalized insectivores, omnivores and some frugivores; and bats that specialize on hard canopy fruits. Only two optima, one corresponding to nectar feeding, were consistently uncovered for von Mises stress. These results suggest that mechanical advantage played a larger role than von Mises stress in the radiation of New World Leaf‐nosed bats into divergent adaptive zones.     

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.