Adaptive processes drive ecomorphological convergent evolution in antwrens (Thamnophilidae)

Phylogenetic niche conservatism (PNC) and convergence are contrasting evolutionary patterns that describe phenotypic similarity across independent lineages. Assessing whether and how adaptive processes give origin to these patterns represent a fundamental step toward understanding phenotypic evolution. Phylogenetic model‐based approaches offer the opportunity not only to distinguish between PNC and convergence, but also to determine the extent that adaptive processes explain phenotypic similarity. The Myrmotherula complex in the Neotropical family Thamnophilidae is a polyphyletic group of sexually dimorphic small insectivorous forest birds that are relatively homogeneous in size and shape. Here, we integrate a comprehensive species‐level molecular phylogeny of the Myrmotherula complex with morphometric and ecological data within a comparative framework to test whether phenotypic similarity is described by a pattern of PNC or convergence, and to identify evolutionary mechanisms underlying body size and shape evolution. We show that antwrens in the Myrmotherula complex represent distantly related clades that exhibit adaptive convergent evolution in body size and divergent evolution in body shape. Phenotypic similarity in the group is primarily driven by their tendency to converge toward smaller body sizes. Differences in body size and shape across lineages are associated to ecological and behavioral factors.     

Phenotypic integration between claw and toepad traits promotes microhabitat specialization in the Anolis adaptive radiation

The performance of an organism in its environment frequently depends more on its composite phenotype than on individual phenotypic traits. Thus, understanding environmental adaptation requires investigating patterns of covariation across functionally related traits. The replicated adaptive radiations of Greater Antillean Anolis lizards are characterized by ecological and morphological convergence, thus, providing an opportunity to examine the role of multiple phenotypes in microhabitat adaptation. Here, we examine integrated claw and toepad morphological evolution in relation to habitat partitioning across the adaptive radiations of Greater Antillean anoles. Based on analysis of 428 specimens from 57 species, we found that different aspects of claw morphology were associated with different perch dimensions, with claw height positively associated with perch diameter and claw curvature positively associated with perch height. Patterns of integration also varied across claw and toepad traits, likely driven by correlative selection for performance on smoother and rougher substrates. Finally, rates of evolution differed between claw and toepad traits, with claw length evolving faster than all other traits despite having no predicted functional importance. Our results highlight the multivariate nature of phenotypic adaptation and suggest that phenotypic integration across Greater Antillean anoles is driven by fine‐scale correlative selection based on structural habitat specialization.     

Discordant patterns of genetic and phenotypic differentiation in five grasshopper species codistributed across a microreserve network

Conservation plans can be greatly improved when information on the evolutionary and demographic consequences of habitat fragmentation is available for several codistributed species. Here, we study spatial patterns of phenotypic and genetic variation among five grasshopper species that are codistributed across a network of microreserves but show remarkable differences in dispersal‐related morphology (body size and wing length), degree of habitat specialization and extent of fragmentation of their respective habitats in the study region. In particular, we tested the hypothesis that species with preferences for highly fragmented microhabitats show stronger genetic and phenotypic structure than codistributed generalist taxa inhabiting a continuous matrix of suitable habitat. We also hypothesized a higher resemblance of spatial patterns of genetic and phenotypic variability among species that have experienced a higher degree of habitat fragmentation due to their more similar responses to the parallel large‐scale destruction of their natural habitats. In partial agreement with our first hypothesis, we found that genetic structure, but not phenotypic differentiation, was higher in species linked to highly fragmented habitats. We did not find support for congruent patterns of phenotypic and genetic variability among any studied species, indicating that they show idiosyncratic evolutionary trajectories and distinctive demographic responses to habitat fragmentation across a common landscape. This suggests that conservation practices in networks of protected areas require detailed ecological and evolutionary information on target species to focus management efforts on those taxa that are more sensitive to the effects of habitat fragmentation.     

Diversity trends and their ontogenetic basis: an exploration of allometric disparity in rodents

It has been hypothesized that most morphological evolution occurs by allometric differentiation. Because rodents encapsulate a phenomenal amount of taxonomic diversity and, among several clades, contrasting levels of morphological diversity, they represent an excellent subject to address the question: how variable are allometric patterns during evolution? We investigated the influence of phylogenetic relations and ecological factors on the results of the first quantification of allometric disparity among rodents by exploring allometric space, a multivariate morphospace here derived from, and encapsulating all, the ontogenetic trajectories of 34 rodent species from two parallel phylogenetic radiations. Disparity was quantified using angles between ontogenetic trajectories for different species and clades. We found an overlapping occupation of allometric space by muroid and hystricognath species, revealing both clades possess similar abilities to evolve in different directions of phenotypic space, and anatomical diversity does not act to constrain the labile nature of allometric patterning. Morphological features to enable efficient processing of food serve to group rodents in allometric space, reflecting the importance of convergent morphology, rather than shared evolutionary history, in the generation of allometric patterns. Our results indicate that the conserved level of morphological integration found among primates cannot simply be extended to all mammals.     

An evolutionary perspective on leaf economics: phylogenetics of leaf mass per area in vascular plants

In plant leaves, resource use follows a trade‐off between rapid resource capture and conservative storage. This “worldwide leaf economics spectrum” consists of a suite of intercorrelated leaf traits, among which leaf mass per area, LMA, is one of the most fundamental as it indicates the cost of leaf construction and light‐interception borne by plants. We conducted a broad‐scale analysis of the evolutionary history of LMA across a large dataset of 5401 vascular plant species. The phylogenetic signal in LMA displayed low but significant conservatism, that is, leaf economics tended to be more similar among close relatives than expected by chance alone. Models of trait evolution indicated that LMA evolved under weak stabilizing selection. Moreover, results suggest that different optimal phenotypes evolved among large clades within which extremes tended to be selected against. Conservatism in LMA was strongly related to growth form, as were selection intensity and phenotypic evolutionary rates: woody plants showed higher conservatism in relation to stronger stabilizing selection and lower evolutionary rates compared to herbaceous taxa. The evolutionary history of LMA thus paints different evolutionary trajectories of vascular plant species across clades, revealing the coordination of leaf trait evolution with growth forms in response to varying selection regimes.     

Quantification and decomposition of environment‐selection relationships

In nature, selection varies across time in most environments, but we lack an understanding of how specific ecological changes drive this variation. Ecological factors can alter phenotypic selection coefficients through changes in trait distributions or individual mean fitness, even when the trait‐absolute fitness relationship remains constant. We apply and extend a regression‐based approach in a population of Soay sheep (Ovis aries) and suggest metrics of environment‐selection relationships that can be compared across studies. We then introduce a novel method that constructs an environmentally structured fitness function. This allows calculation of full (as in existing approaches) and partial (acting separately through the absolute fitness function slope, mean fitness, and phenotype distribution) sensitivities of selection to an ecological variable. Both approaches show positive overall effects of density on viability selection of lamb mass. However, the second approach demonstrates that this relationship is largely driven by effects of density on mean fitness, rather than on the trait‐fitness relationship slope. If such mechanisms of environmental dependence of selection are common, this could have important implications regarding the frequency of fluctuating selection, and how previous selection inferences relate to longer term evolutionary dynamics.     

Phylogenetic convergence and multiple shell shape optima for gliding scallops (Bivalvia: Pectinidae)

An important question in evolutionary biology is how often, and to what extent, do similar ecologies elicit distantly related taxa to evolve towards the same phenotype? In some scenarios, the repeated evolution of particular phenotypes may be expected, for instance when species are exposed to common selective forces that result from strong functional demands. In bivalved scallops (Pectinidae), some species exhibit a distinct swimming behaviour (gliding), which requires specific biomechanical attributes to generate lift and reduce drag during locomotive events. Further, a phylogenetic analysis revealed that gliding behaviour has independently evolved at least four times, which raises the question as to whether these independent lineages have also converged on a similar phenotype. Here, we test the hypothesis that gliding scallops display shell shape convergence using a combination of geometric morphometrics and phylogenetic comparative methods that evaluate patterns of multivariate trait evolution. Our findings reveal that the gliding species display less morphological disparity and significant evolutionary convergence in morphospace, relative to expectations under a neutral model of Brownian motion for evolutionary phenotypic change. Intriguingly, the phylomorphospace patterns indicate that gliding lineages follow similar evolutionary trajectories to not one, but two regions of morphological space, and subsequent analyses identified significant differences in their biomechanical parameters, suggesting that these two groups of scallops accomplish gliding in different ways. Thus, whereas there is a clear gliding morphotype that has evolved convergently across the phylogeny, functionally distinct morphological subforms are apparent, suggesting that there may be two optima for the gliding phenotype in the Pectinidae.     

Genetic covariance structure of growth in the salamander Ambystoma macrodactylum

The size of an organism at any point during ontogeny often has fitness consequences through either direct selection on size or through selection on size-related morphological, performance, or life history traits. However, the evolutionary response to selection on size across ontogeny (a growth trajectory) may be limited by genetic correlations across ages. Here we characterize the phenotypic and genetic covariance structure of length and mass growth trajectories in a natural population of larval Ambystoma macrodactylum using function-valued quantitative genetic analyses and principal component decomposition. Most of the phenotypic and genetic variation in both growth trajectories appears to be confined to a single principal component describing a pattern of positive covariation among sizes across all ages. Higher order principal components with no significant associated genetic variation were identified for both trajectories, suggesting that evolution towards certain patterns of negative covariation between sizes across ages is constrained. The well-characterized positive relationship between size at metamorphosis and fitness in pond-breeding amphibians predicts that the across-age covariance structure will strongly limit evolution only if there is negative selection on size prior to metamorphosis. The pattern of genetic covariation observed in this study is similar to that observed in other vertebrate taxa, indicating that size may often be highly genetically and phenotypically integrated across ontogeny. Additionally, we find that phenotypic and genetic analyses of growth trajectories can yield qualitatively similar patterns of covariance structure.     

Predicting evolution with generalized models of divergent selection: a case study with poeciliid fish

Over the past century and half since the process of natural selection was first described, one enduring question has captivated many, "how predictable is evolution?" Because natural selection comprises deterministic components, the course of evolution may exhibit some level of predictability across organismal groups. Here, I provide an early appraisal of the utility of one particular approach to understanding the predictability of evolution: generalized models of divergent selection (GMDS). The GMDS approach is meant to provide a unifying framework for the science of evolutionary prediction, offering a means of better understanding the causes and consequences of phenotypic and genetic evolution. I describe and test a GMDS centered on the evolution of body shape, size of the gonopodium (sperm-transfer organ), steady-swimming abilities, fast-start swimming performance, and reproductive isolation between populations in Gambusia fishes (Family Poeciliidae). The GMDS produced some accurate evolutionary predictions in Gambusia, identifying variation in intensity of predation by piscivorous fish as a major factor driving repeatable and predictable phenotypic divergence, and apparently playing a key role in promoting ecological speciation. Moreover, the model's applicability seems quite general, as patterns of differentiation in body shape between predator regimes in many disparate fishes match the model's predictions. The fact that such a simple model could yield accurate evolutionary predictions in distantly related fishes inhabiting different geographic regions and types of habitat, and experiencing different predator species, suggests that the model pinpointed a causal factor underlying major, shared patterns of diversification. The GMDS approach appears to represent a promising method of addressing the predictability of evolution and identifying environmental factors responsible for driving major patterns of replicated evolution.     

Future directions in the ontogeny of plant defence: understanding the evolutionary causes and consequences

Plant defence often varies by orders of magnitude as plants develop from the seedling to juvenile to mature and senescent stages. Ontogenetic trajectories can involve switches among defence traits, leading to complex shifting phenotypes across plant lifetimes. While considerable research has characterised ontogenetic trajectories for now hundreds of plant species, we still lack a clear understanding of the molecular, ecological and evolutionary factors driving these patterns. In this study, we identify several non‐mutually exclusive factors that may have led to the evolution of ontogenetic trajectories in plant defence, including developmental constraints, resource allocation costs, multi‐functionality of defence traits, and herbivore selection pressure. Evidence from recent physiological studies is highlighted to shed light on the underlying molecular mechanisms involved in the regulation and activation of these developmental changes. Overall, our goal is to promote new research avenues that would provide evidence for the factors that have promoted the evolution of this complex lifetime phenotype. Future research focusing on the questions and approaches identified here will advance the field and shed light on why defence traits shift so dramatically across plant ontogeny, a widespread but poorly understood ecological pattern.     

A non‐invasive geometric morphometrics method for exploring variation in dorsal head shape in urodeles: sexual dimorphism and geographic variation in Salamandra salamandra

The study of morphological variation among and within taxa can shed light on the evolution of phenotypic diversification. In the case of urodeles, the dorso‐ventral view of the head captures most of the ontogenetic and evolutionary variation of the entire head, which is a structure with a high potential for being a target of selection due to its relevance in ecological and social functions. Here, we describe a non‐invasive procedure of geometric morphometrics for exploring morphological variation in the external dorso‐ventral view of urodeles' head. To explore the accuracy of the method and its potential for describing morphological patterns we applied it to two populations of Salamandra salamandra gallaica from NW Iberia. Using landmark‐based geometric morphometrics, we detected differences in head shape between populations and sexes, and an allometric relationship between shape and size. We also determined that not all differences in head shape are due to size variation, suggesting intrinsic shape differences across sexes and populations. These morphological patterns had not been previously explored in S. salamandra , despite the high levels of intraspecific diversity within this species. The methodological procedure presented here allows to detect shape variation at a very fine scale, and solves the drawbacks of using cranial samples, thus increasing the possibilities of using collection specimens and alive animals for exploring dorsal head shape variation and its evolutionary and ecological implications in urodeles. J. Morphol. 278:475–485, 2017. © 2017 Wiley Periodicals, Inc.     

The relative roles of adaptation and phylogeny in determination of larval traits in diversifying anuran lineages

I measured phenotypic traits important to the fitness of larval anurans to assess the relative roles of ancestral trait value and selective regime in determining present-day phenotypes. The positions of 14 species from three taxonomic families and three different habitats in a phenotypic space defined by 19 traits provided measures of taxonomic and ecological similarity. The distribution of phenotypic distances among species revealed that neither taxonomy nor habitat overwhelmingly determined phenotype. There appear to be multiple ways in which anurans can exploit pond types. However, the direction of phenotypic movement was not random from one species to the next. Independent contrasts revealed significant correlations in the evolution of traits that were consistent among lineages. These correlations reflected well-known trade-offs that result from functional relationships among the constituent traits. Although there is no simple pattern in the distribution of mean phenotypes across environments and lineages, the pattern of the evolutionary trajectories that created that distribution is consistent with a predictive theory of multivariate evolution.     

Evolution in a changing environment: a case study with great tit fledging mass

Heritable phenotypic traits under significant and consistent directional selection often fail to show the expected evolutionary response. A potential explanation for this contradiction is that because environmental conditions change constantly, environmental change can mask an evolutionary response to selection. We combined an "animal model" analysis with 36 years of data from a long-term study of great tits (Parus major) to explore selection on and evolution of a morphological trait: body mass at fledging. We found significant heritability of this trait, but despite consistent positive directional selection on both the phenotypic and the additive genetic component of body mass, the population mean phenotypic value declined rather than increased over time. However, the mean breeding value for body mass at fledging increased over time, presumably in response to selection. We show that the divergence between the response to selection observed at the levels of genotype and phenotype can be explained by a change in environmental conditions over time, that is, related both to increased spring temperature before breeding and elevated population density. Our results support the suggestion that measuring phenotypes may not always give a reliable impression of evolutionary trajectories and that understanding patterns of phenotypic evolution in nature requires an understanding of how the environment has itself changed.     

Sex,competition and mimicry: an eco-evolutionary model reveals unexpected impacts of ecological interactions on the evolution of phenotypes in sympatry

Phenotypic evolution in sympatric species can be strongly impacted by species interactions, either mutualistic or antagonistic. Heterospecific reproductive behaviours between sympatric species have been shown to favour phenotypic divergence of traits used as sexual cues. Those traits may also be involved in local adaptation or in other types of species interactions and, as a result, undergo complex evolutions across sympatric species. Here we focus on mimicry and study how reproductive interference may impair phenotypic convergence between species with various levels of defence. We use a deterministic model assuming two sympatric species where individuals can display two different warning colour patterns. This eco-evolutionary model explores how ecological interactions shape phenotypic evolution within sympatric species. We investigate the effect of 1) the opposing density-dependent selections exerted on colour patterns by predation and reproductive behaviour and 2) the impact of relative species and phenotype abundances on the fitness costs faced by each individual depending on their species and phenotype. Our model shows that reproductive interference may limit the convergent effect of mimetic interactions and may promote phenotypic divergence between Müllerian mimics. The divergent and convergent evolution of traits also strongly depends on the relative species and phenotype abundances and levels of trophic competition, highlighting how the eco-evolutionary feedbacks between phenotypic evolution and species abundances may result in strikingly different evolutionary routes.     

Patterns of univariate and multivariate plasticity to elevated carbon dioxide in six European populations of Arabidopsis thaliana

The impact of elevated carbon dioxide on plants is a growing concern in evolutionary ecology and global change biology. Characterizing patterns of phenotypic integration and multivariate plasticity to elevated carbon dioxide can provide insights into ecological and evolutionary dynamics in future human‐altered environments. Here, we examined univariate and multivariate responses to carbon enrichment in six functional traits among six European accessions of Arabidopsis thaliana. We detected phenotypic plasticity in both univariate and multivariate phenotypes, but did not find significant variation in plasticity (genotype by environment interactions) within or among accessions. Eigenvector, eigenvalue variance, and common principal components analyses showed that elevated carbon dioxide altered patterns of trait covariance, reduced the strength of phenotypic integration, and decreased population‐level differentiation in the multivariate phenotype. Our data suggest that future carbon dioxide conditions may influence evolutionary dynamics in natural populations of A. thaliana.     

Variation,selection and evolution of function-valued traits

We describe an emerging framework for understanding variation, selection and evolution of phenotypic traits that are mathematical functions. We use one specific empirical example – thermal performance curves (TPCs) for growth rates of caterpillars – to demonstrate how models for function-valued traits are natural extensions of more familiar, multivariate models for correlated, quantitative traits. We emphasize three main points. First, because function-valued traits are continuous functions, there are important constraints on their patterns of variation that are not captured by multivariate models. Phenotypic and genetic variation in function-valued traits can be quantified in terms of variance-covariance functions and their associated eigenfunctions: we illustrate how these are estimated as well as their biological interpretations for TPCs. Second, selection on a function-valued trait is itself a function, defined in terms of selection gradient functions. For TPCs, the selection gradient describes how the relationship between an organism's performance and its fitness varies as a function of its temperature. We show how the form of the selection gradient function for TPCs relates to the frequency distribution of environmental states (caterpillar temperatures) during selection. Third, we can predict evolutionary responses of function-valued traits in terms of the genetic variance-covariance and the selection gradient functions. We illustrate how non-linear evolutionary responses of TPCs may occur even when the mean phenotype and the selection gradient are themselves linear functions of temperature. Finally, we discuss some of the methodological and empirical challenges for future studies of the evolution of function-valued traits.     

Developmental Dynamics and G-Matrices: Can Morphometric Spaces be Used to Model Phenotypic Evolution?

Modern morphometrics, especially geometric morphometrics, is a powerful tool for modeling the evolution and development of the phenotype. Complicated morphological transformations can be simulated by using standard evolutionary genetic equations for processes such as selection and drift in the same morphospaces that are used for empirical morphometric studies. Such applications appear to be consistent with the theory of quantitative evolution of the phenotype. Nevertheless, concerns exist whether simulations of phenotypic changes directly in morphospaces is realistic because trajectories traced in such spaces describe continuous gradations in the phenotype and because the gain and loss of structures is often impossible because morphospaces are necessarily constructed from variables shared in common by all the phenotypes being considered. Competing models of phenotypic change emphasize morphological discontinuity and novelty. Recently developed models of phenotypic evolution that introduce a “phenotypic landscape” between evolutionary genetic constructs like the adaptive landscape and morphospace may correct this shortcoming.     

Selection experiments and the study of phenotypic plasticity1

Abstract Laboratory selection experiments are powerful tools for establishing evolutionary potentials. Such experiments provide two types of information, knowledge about genetic architecture and insight into evolutionary dynamics. They can be roughly classified into two types: (1) artificial selection in which the experimenter selects on a focal trait or trait index, and (2) quasi‐natural selection in which the experimenter establishes a set of environmental conditions and then allows the population to evolve. Both approaches have been used in the study of phenotypic plasticity. Artificial selection experiments have taken various forms including: selection directly on a reaction norm, selection on a trait in multiple environments, and selection on a trait in a single environment. In the latter experiments, evolution of phenotypic plasticity is investigated as a correlated response. Quasi‐natural selection experiments have examined the effects of both spatial and temporal variation. I describe how to carry out such experiments, summarize past efforts, and suggest further avenues of research.     

Reaction norms of Arabidopsis. V. Flowering time controls phenotypic architecture in response to nutrient stress

The evolutionary and environmental stability of character correlations has increasingly been the focus of ecological and quantitative genetic studies. Although the genetic stability of character correlations is a central assumption of quantitative genetic models of phenotypic evolution, theoretical considerations suggest that both the genetic and the phenotypic architecture should change in response to selection and to environmental heterogeneity. We investigate genetic (population) differences and plasticity to nutrient availability of the phenotypic architecture describing the whole-plant phenotype of Arabidopsis thaliana (Brassicaceae). We found significant genetic differences among early and late flowering ecotypes in the relationships between several traits, when a path-analytical model was used to estimate character correlations. Furthermore, we found significant plasticity of several path coefficients when nutrient levels were altered. A whole-plant analysis considering all paths in the model simultaneously confirmed that populations of A. thaliana are characterized by distinct phenotypic architectures, and that these are altered in different ways by environmental changes. We discuss the implications of these findings for our understanding of selective pressure on and response by multivariate phenotypes.     

Phenotypic variation in growth trajectories in the Arctic charr Salvelinus alpinus

Animals with determinate growth have shown little variation in individual growth patterns, but similar analyses for animals with indeterminate growth have been lacking. We analysed the amount of phenotypic variation in growth patterns across ages among individuals of a hatchery-based population of Arctic charr, Salvelinus alpinus, Salmonidae, using the infinite-dimensional model and including the effects of group size structure. There was little phenotypic variation in growth trajectories: individuals that were small (in relation to the mean) early in life were among the smallest 2.5 years later. If the genetic variation reflects phenotypic variation, not much evolutionary change can be expected. Our results show that there are ecological conditions that determine the strong covariation of size across ages, most likely size-related dominance behaviour, which can mask the true variation of growth patterns. Thus, social interactions can have strong evolutionary effects on traits not directly involved in the behavioural interactions.