Resource-based trade-offs and the adaptive significance of seasonal plasticity in butterfly wing melanism

Phenotypic plasticity is the ability of an organism to alter its phenotype in response to environmental cues. This can be adaptive if the cues are reliable predictors of impending conditions and the alterations enhance the organism's ability to capitalize on those conditions. However, since traits do not exist in isolation but as part of larger interdependent systems of traits (phenotypic integration), trade-offs between correlated plastic traits can make phenotypic plasticity non- or maladaptive. We examine this problem in the seasonally plastic wing melanism of a pierid (Order Lepidoptera, Family Pieridae) butterfly, Pieris rapae L. Several wing pattern traits are more melanized in colder than in warmer seasons, resulting in effective thermoregulation through solar absorption. However, other wing pattern traits, the spots, are less melanized during colder seasons than in warmer seasons. Although spot plasticity may be adaptive, reduced melanism of these spots could also be explained by resource-based trade-offs. Theory predicts that traits involved in resource-based trade-offs will be positively correlated when variation among individuals in resource acquisition is greater than variation among individuals in resource allocation strategies, and negatively correlated when variation in allocation is greater than variation in acquisition. Using data from both field studies and laboratory studies that manipulate dietary tyrosine, a melanin precursor, we show that when allocation to thermoregulatory melanism (ventral hindwing, and basal dorsal fore- and hindwing "shading") varies substantially this trait is negatively correlated with spot melanism. However, when there is less variation in allocation to thermoregulatory melanism we find these traits to be positively correlated; these findings are consistent with the resource-based trade-off hypothesis, which may provide a non- or maladaptive hypothesis to explain spot plasticity. We also show that increased dietary tyrosine results in increased spot melanism under some conditions, supporting the more general idea that melanism may involve resource-based costs.

Distribution theories for genetic line of least resistance and evolvability measures

Quantitative genetic theory on multivariate character evolution predicts that a population's response to directional selection is biased toward the major axis of the genetic covariance matrix G-the so-called genetic line of least resistance. Inferences on the genetic constraints in this sense have traditionally been made by measuring the angle of deviation of evolutionary trajectories from the major axis, or more recently by calculating the amount of genetic variance-the Hansen-Houle evolvability-available along the trajectories. However, there have not been clear practical guidelines on how these quantities can be interpreted, especially in a high-dimensional space. This study summarizes pertinent distribution theories for relevant quantities, pointing out that they can be written as ratios of quadratic forms in evolutionary trajectory vectors by taking G as a parameter. For example, a beta distribution with appropriate parameters can be used as a null distribution for squared cosine of the angle of deviation from a major axis or subspace. More general cases can be handled with the probability distribution of ratios of quadratic forms in normal variables. Apart from its use in hypothesis-testing, this latter approach could potentially be used as a heuristic tool for looking into various selection scenarios like directional and/or correlated selection as parameterized with mean and covariance of selection gradients.

Climate change and its effects on body size and shape: the role of endocrine mechanisms

In many organisms, rapidly changing environmental conditions are inducing dramatic shifts in diverse phenotypic traits with consequences for fitness and population viability. However, the mechanisms that underlie these responses remain poorly understood. Endocrine signalling systems often influence suites of traits and are sensitive to changes in environmental conditions; they are thus ideal candidates for uncovering both plastic and evolved consequences of climate change. Here, we use body size and shape, a set of integrated traits predicted to shift in response to rising temperatures with effects on fitness, and insulin-like growth factor-1 as a case study to explore these ideas. We review what is known about changes in body size and shape in response to rising temperatures and then illustrate why endocrine signalling systems are likely to be critical in mediating these effects. Lastly, we discuss research approaches that will advance understanding of the processes that underlie rapid responses to climate change and the role endocrine systems will have. Knowledge of the mechanisms involved in phenotypic responses to climate change will be essential for predicting both the ecological and the long-term evolutionary consequences of a warming climate. This article is part of the theme issue 'Endocrine responses to environmental variation: conceptual approaches and recent developments'.

Why incorporate plant architecture into trait-based ecology?

Trait-based ecology has improved our understanding of the functioning of organisms, communities, ecosystems, and beyond. However, its predictive ability remains limited as long as phenotypic integration and temporal dynamics are not considered. We highlight how the morphogenetic processes that shape the 3D development of a plant during its lifetime affect its performance. We show that the diversity of architectural traits allows us to go beyond organ-level traits in capturing the temporal and spatial dimensions of ecological niches and informing community assembly processes. Overall, we argue that consideration of multilevel topological, geometrical, and ontogenetic features provides a dynamic view of the whole-plant phenotype and a relevant framework for investigating phenotypic integration, plant adaptation and performance, and community structure and dynamics.

Within-individual leaf trait variation increases with phenotypic integration in a subtropical tree diversity experiment

Covariation of plant functional traits, that is, phenotypic integration, might constrain their variability. This was observed for inter- and intraspecific variation, but there is no evidence of a relationship between phenotypic integration and the functional variation within single plants (within-individual trait variation; WTV), which could be key to understand the extent of WTV in contexts like plant-plant interactions. We studied the relationship between WTV and phenotypic integration in c. 500 trees of 21 species in planted forest patches varying in species richness in subtropical China. Using visible and near-infrared spectroscopy (Vis-NIRS), we measured nine leaf morphological and chemical traits. For each tree, we assessed metrics of single and multitrait variation to assess WTV, and we used plant trait network properties based on trait correlations to quantify phenotypic integration. Against expectations, strong phenotypic integration within a tree led to greater variation across leaves. Not only this was true for single traits, but also the dispersion in a tree's multitrait hypervolume was positively associated with tree's phenotypic integration. Surprisingly, we only detected weak influence of the surrounding tree-species diversity on these relationships. Our study suggests that integrated phenotypes allow the variability of leaf phenotypes within the organism and supports that phenotypic integration prevents maladaptive variation.

Covariation of brain and skull shapes as a model to understand the role of crosstalk in development and evolution

Covariation among discrete phenotypes can arise due to selection for shared functions, and/or shared genetic and developmental underpinnings. The consequences of such phenotypic integration are far-reaching and can act to either facilitate or limit morphological variation. The vertebrate brain is known to act as an "organizer" of craniofacial development, secreting morphogens that can affect the shape of the growing neurocranium, consistent with roles for pleiotropy in brain-neurocranium covariation. Here, we test this hypothesis in cichlid fishes by first examining the degree of shape integration between the brain and the neurocranium using three-dimensional geometric morphometrics in an F₅ hybrid population, and then genetically mapping trait covariation using quantitative trait loci (QTL) analysis. We observe shape associations between the brain and the neurocranium, a pattern that holds even when we assess associations between the brain and constituent parts of the neurocranium: the rostrum and braincase. We also recover robust genetic signals for both hard- and soft-tissue traits and identify a genomic region where QTL for the brain and braincase overlap, implicating a role for pleiotropy in patterning trait covariation. Fine mapping of the overlapping genomic region identifies a candidate gene, notch1a, which is known to be involved in patterning skeletal and neural tissues during development. Taken together, these data offer a genetic hypothesis for brain-neurocranium covariation, as well as a potential mechanism by which behavioral shifts may simultaneously drive rapid change in neuroanatomy and craniofacial morphology.

Broad diversity in monoterpene-sesquiterpene balance across wild sunflowers: Implications of leaf and floral volatiles for biotic interactions

PREMISE

As plant lineages diversify across environmental gradients, species are predicted to encounter divergent biotic pressures. This study investigated the evolution of volatile secondary metabolism across species of Helianthus.

METHODS

Leaves and petals of 40 species of wild Helianthus were analyzed via gas chromatography-mass spectrometry to determine volatile secondary metabolite profiles.

RESULTS

Across all species, 500 compounds were identified; 40% were sesquiterpenes, 18% monoterpenes, 3% diterpenes, 4% fatty acid derivatives, and 35% other compounds such as phenolics and small organic molecules. Qualitatively, annuals and species from more arid western climates had leaf compositions with a higher proportion of total monoterpenes, while erect perennials and species from more mesic eastern habitats contained a higher proportion of total sesquiterpenes. Among species, mass-based leaf monoterpene and sesquiterpene abundance were identified as largely orthogonal axes of variation by principal component analysis. Profiles for leaves were not strongly correlated with those of petals.

CONCLUSIONS

Volatile metabolites were highly diverse among wild Helianthus, indicating the value of this genus as a model system and rich genetic resource. The independence of leaf and petal volatile profiles indicates a low level of phenotypic integration between vegetative and reproductive structures, implying vegetative defense and reproductive defense or pollinator attraction functions mediated by terpene profiles in these two organs can evolve without major trade-offs. The major biosynthetic pathways for the major terpenes in wild Helianthus are already well described, providing a road map to deeper inquiry into the drivers of this diversity.

Comparative morpho-functional analysis of the humerus and ulna in three Western European moles species of the genus Talpa, including the newly described T. aquitania

The forelimb is involved in many behaviours including locomotion. Notably, the humero-ulnar articulation, implicated in the elbow joint, is of particular importance for both mobility and stability. Functional constraints, induced in part by environmental plasticity, are thought to drive an important part of the bone shape as bone directly responds and remodels in response to both muscle and external forces. In this context, the study of subterranean moles is of particular interest. These moles occupy a hard and heavy medium in comparison with air or water, requiring a powerful body structure to shear and shift the soil. Their general morphology is therefore adapted to digging and to their subterranean lifestyle. The various morpho-functional patterns, which drive diverse abilities according to the environment, are likely targets of natural selection and it is, therefore, useful to understand the relationships between the bone shape and their function. Here, we quantify, through 3D geometric morphometric methods, the interspecific variability in the morphology of the ulna and humerus of three Talpa species, including the new species Talpa aquitania, to infer their potential consequence in species digging performance. We also quantify shape covariation and morphological integration between the humerus and the ulna to test whether these bones evolve as a uniform functional unit or as more or less independent modules. Our results show that interspecific anatomical differences in the humerus and ulna exist among the three species. Shape changes are mostly located at the level of joints and muscle attachments. As the species tend to live in allopatry and the fossorial lifestyle induces strong ecological constraints, interspecific variations could be explained by the properties of the environment in which they live, such as the compactness of the soil. Our results also show that the humerus and ulna are highly integrated. The covariation between the humerus and ulna in moles is dominated by variation in the attachment areas and particularly of the attachment areas of shoulder muscles concerning the humerus, which affect the mechanical force deployed during locomotion and digging. This study also highlights that in the new species, T. aquitania, variations in anatomical structure (general shape and joints) exist and are related to the locality of collect of the individuals.

Environment-driven shifts in inter-individual variation and phenotypic integration within subnetworks of the mussel transcriptome and proteome

The environment can alter the magnitude of phenotypic variation among individuals, potentially influencing evolutionary trajectories. However, environmental influences on variation are complex and remain understudied. Populations in heterogeneous environments might exhibit more variation, the amount of variation could differ between benign and stressful conditions, and/or variation might manifest in different ways among stages of the gene‐to‐protein expression cascade or among physiological functions. Here, we explore these three issues by quantifying patterns of inter‐individual variation in both transcript and protein expression levels among California mussels, Mytilus californianus Conrad. Mussels were exposed to five ecologically relevant treatments that varied in the mean and interindividual heterogeneity of body temperature. To target a diverse set of physiological functions, we assessed variation within 19 expression subnetworks, including canonical stress‐response pathways and empirically derived coexpression clusters that represent a diffuse set of cellular processes. Variation in expression was particularly pronounced in the treatments with high mean and heterogeneous body temperatures. However, with few exceptions, environment‐dependent shifts of variation in the transcriptome were not reflected in the proteome. A metric of phenotypic integration provided evidence for a greater degree of constraint on relative expression levels (i.e., stronger correlation) within expression subnetworks in benign, homogeneous environments. Our results suggest that environments that are more stressful on average – and which also tend to be more heterogeneous – can relax these expression constraints and reduce phenotypic integration within biochemical subnetworks. Context‐dependent “unmasking” of functional variation may contribute to interindividual differences in physiological phenotype and performance in stressful environments.

The role of genus and life span in predicting seed and vegetative trait variation and correlation in Lathyrus, Phaseolus, and Vicia

PREMISE

Annual and perennial life history transitions are abundant among angiosperms, and understanding the phenotypic variation underlying life span shifts is a key endeavor of plant evolutionary biology. Comparative analyses of trait variation and correlation networks among annual and perennial plants is increasingly important as new herbaceous perennial crops are being developed for edible seed. However, it remains unclear how seed to vegetative growth trait relationships correlate with life span.

METHODS

To assess the relative roles of genus and life span in predicting phenotypic variation and trait correlations, we measured seed size and shape, germination proportion, and early-life-stage plant height and leaf growth over 3 mo in 29 annual and perennial, herbaceous congeneric species from three legume genera (Lathyrus, Phaseolus, and Vicia).

RESULTS

Genus was the strongest predictor of seed size and shape variation, and life span consistently predicted plant height and leaf number at single time points. Correlation networks revealed that annual species had significant associations between seed traits and vegetative traits, whereas perennials had no significant seed-vegetative associations. Each genus also differed in the extent of integration between seed and vegetative traits, as well as within-vegetative-trait correlation patterns.

CONCLUSIONS

Genus and life span were important for predicting aspects of early-life-stage phenotypic variation and trait relationships. Differences in phenotypic correlation may indicate that selection on seed size traits will impact vegetative growth differently depending on life span, which has important implications for nascent perennial breeding programs.

On Information Rank Deficiency in Phenotypic Covariance Matrices

This paper investigates a form of rank deficiency in phenotypic covariance matrices derived from geometric morphometric data, and its impact on measures of phenotypic integration. We first define a type of rank deficiency based on information theory, then demonstrate that this deficiency impairs the performance of phenotypic integration metrics in a model system. Lastly, we propose methods to treat for this information rank deficiency. Our first goal is to establish how the rank of a typical geometric morphometric covariance matrix relates to the information entropy of its eigenvalue spectrum. This requires clear definitions of matrix rank, of which we define three: the full matrix rank (equal to the number of input variables), the mathematical rank (the number of non-zero eigenvalues), and the information rank or 'effective rank' (equal to the number of non-redundant eigenvalues). We demonstrate that effective rank deficiency arises from a combination of methodological factors - Generalized Procrustes analysis, use of the correlation matrix, and insufficient sample size - as well as phenotypic covariance. Secondly, we use dire wolf jaws to document how differences in effective rank deficiency bias two metrics used to measure phenotypic integration. The eigenvalue variance characterizes the integration change incorrectly, and the standardized generalized variance lacks the sensitivity needed to detect subtle changes in integration. Both metrics are impacted by the inclusion of many small, but non-zero, eigenvalues arising from a lack of information in the covariance matrix, a problem that usually becomes more pronounced as the number of landmarks increases. We propose a new metric for phenotypic integration that combines the standardized generalized variance with information entropy. This metric is equivalent to the standardized generalized variance, but calculated only from those eigenvalues that carry non-redundant information. It is the standardized generalized variance scaled to the effective rank of the eigenvalue spectrum. We demonstrate that this metric successfully detects the shift of integration in our dire wolf sample. Our third goal is to generalize the new metric to compare data sets with different sample sizes and numbers of variables. We develop a standardization for matrix information based on data permutation, then demonstrate that Smilodon jaws are more integrated than dire wolf jaws. Finally, we describe how our information entropy-based measure allows phenotypic integration to be compared in dense semilandmark data sets without bias, allowing characterization of the information content of any given shape, a quantity we term 'latent dispersion'.

Statistics of eigenvalue dispersion indices: Quantifying the magnitude of phenotypic integration

Analysis of trait covariation plays a pivotal role in the study of phenotypic evolution. The magnitude of covariation is often quantified with statistics based on dispersion of eigenvalues of a covariance or correlation matrix-eigenvalue dispersion indices. This study clarifies the statistical justifications of these statistics and elaborates on their sampling properties. The relative eigenvalue variance of a covariance matrix is known in the statistical literature a test statistic for sphericity, and thus is an appropriate measure of eccentricity of variation. The same of a correlation matrix is equal to the average squared correlation, which has a straightforward interpretation as a measure of integration. Here, expressions for the mean and variance of these statistics are analytically derived under multivariate normality, clarifying the effects of sample size N, number of variables p, and parameters on sampling bias and error. Simulations confirm that approximations involved are reasonably accurate with a moderate sample size (N ≥ 16-64). Importantly, sampling properties of these indices are not adversely affected by a high p:N ratio, promising their utility in high-dimensional phenotypic analyses. They can furthermore be applied to shape variables and phylogenetically structured data with appropriate modifications.

Integration between the cranial boundaries of the nasopharynx and the upper cervical vertebrae in Homo and Pan

The nasopharynx is an important anatomical structure involved in respiration. Its bony boundaries, including the basicranium and upper cervical vertebrae, may be subject to selective pressures and constraints related to respiratory function. Here, we investigate phenotypic integration, or covariation, between the face, the basicranial boundaries of the nasopharynx, and the atlas and axis to understand constraints affecting these structures. We collected three-dimensional coordinate data from a sample of 80 humans and 44 chimpanzees, and used two-block partial least squares to assess RV (a multivariate generalization of Pearson's r² ), r_PLS , the covariance ratio, and effect size for integration among structures. We find that integration is significant among some of these structures, and that integration between the basicranial nasopharynx and vertebrae and between the face and vertebrae is likely independent. We also find divergences in the pattern of integration between humans and chimpanzees suggesting greater constraints among the human face and nasopharynx, which we suggest are linked to divergent developmental trajectories in the two taxa. Evolutionary changes in human basicranial anatomy, coupled with human-like developmental trajectories, may have required that the face grow to compensate any variation in nasopharyngeal structure. However, we were unable to determine whether the nasopharynx or the face is more strongly integrated with the vertebrae, and therefore whether respiration or biomechanical considerations related to positional behavior may be more strongly tied to vertebral evolution. Future work should focus on greater sample sizes, soft tissue structures, and more diverse taxa to further clarify these findings.

Morphological canalization, integration, and plasticity in response to population density in Abutilon theophrasti: Influences of soil conditions and growth stages

Phenotypic integration and developmental canalization have been hypothesized to constrain the degree of phenotypic plasticity, but little evidence exists, probably due to the lack of studies on the relationships among the three processes, especially for plants under different environments. We conducted a field experiment by subjecting plants of Abutilon theophrasti to three densities, under infertile and fertile soil conditions, and analyzing correlations among canalization, integration, and plasticity in a variety of measured morphological traits after 50 and 70 days, to investigate the relationships among the three variables in response to density and how these responses vary with soil conditions and growth stages. Results showed trait canalization decreased and phenotypic integration and the degree of plasticity (absolute plasticity) in traits increased with density. Phenotypic integration often positively correlated with absolute plasticity, whereas correlations between trait canalization and plasticity were insignificant in most cases, with a few positive ones between canalization and absolute plasticity at low and medium densities. As plants grew, these correlations intensified in infertile soil and attenuated in fertile soil. Our findings suggested the complexity of the relationship between canalization and plasticity: Decreased canalization is more likely to facilitate active plastic responses under more favorable conditions, whereas increased level of integration should mainly be an outcome of plastic responses. Soil conditions and growth stage may affect responses of these correlations to density via modifying plant size, competition strength, and plastic responses in traits. We also predicted that decreased canalization can be advantageous or disadvantageous, and the lack of response to stress may demonstrate a stronger ability of adaptation than passive response, thus should be adaptive plasticity as active response.

Phenotypic integration does not constrain phenotypic plasticity: differential plasticity of traits is associated to their integration across environments

Understanding constraints to phenotypic plasticity is key given its role on the response of organisms to environmental change. It has been suggested that phenotypic integration, the structure of trait covariation, could limit trait plasticity. However, the relationship between plasticity and integration is far from resolved. Using a database of functional plasticity to drought of a Mediterranean shrub that included 20 ecophysiological traits, we assessed environmentally-induced changes in phenotypic integration and whether integration constrained the expression of plasticity, accounting for the within-environment phenotypic variation of traits. Furthermore, we provide the first test of the association between differential trait plasticity and trait integration across an optimum and a stressful environment. Phenotypic plasticity was positively associated with phenotypic integration in both environments, but this relationship was lost when phenotypic variation was considered. The similarity in the plastic response of two traits predicted their integration across environments, with integrated traits having more similar plasticity. Such variation in the plasticity of traits partly explained the lower phenotypic integration found in the stressful environment. We found no evidence that integration may constitute an internal constraint to plasticity. Rather, we present the first empirical demonstration that differences in plastic responses may involve a major reorganization of the relationships among traits, and challenge the notion that stress generally induces a tighter phenotype.

Novel neuroanatomical integration and scaling define avian brain shape evolution and development

How do large and unique brains evolve? Historically, comparative neuroanatomical studies have attributed the evolutionary genesis of highly encephalized brains to deviations along, as well as from, conserved scaling relationships among brain regions. However, the relative contributions of these concerted (integrated) and mosaic (modular) processes as drivers of brain evolution remain unclear, especially in non-mammalian groups. While proportional brain sizes have been the predominant metric used to characterize brain morphology to date, we perform a high-density geometric morphometric analysis on the encephalized brains of crown birds (Neornithes or Aves) compared to their stem taxa—the non-avialan coelurosaurian dinosaurs and Archaeopteryx. When analyzed together with developmental neuroanatomical data of model archosaurs (Gallus, Alligator), crown birds exhibit a distinct allometric relationship that dictates their brain evolution and development. Furthermore, analyses by neuroanatomical regions reveal that the acquisition of this derived shape-to-size scaling relationship occurred in a mosaic pattern, where the avian-grade optic lobe and cerebellum evolved first among non-avialan dinosaurs, followed by major changes to the evolutionary and developmental dynamics of cerebrum shape after the origin of Avialae. Notably, the brain of crown birds is a more integrated structure than non-avialan archosaurs, implying that diversification of brain morphologies within Neornithes proceeded in a more coordinated manner, perhaps due to spatial constraints and abbreviated growth period. Collectively, these patterns demonstrate a plurality in evolutionary processes that generate encephalized brains in archosaurs and across vertebrates.

Complex interactions underlie the correlated evolution of floral traits and their association with pollinators in a clade with diverse pollination systems

Natural selection by pollinators is an important factor in the morphological diversity and adaptive radiation of flowering plants. Selection by similar pollinators in unrelated plants leads to convergence in floral morphology, or "floral syndromes." Previous investigations into floral syndromes have mostly studied relatively small and/or simple systems, emphasizing vertebrate pollination. Despite the importance of multiple floral traits in plant-pollinator interactions, these studies have examined few quantitative traits, so their co-variation and phenotypic integration have been underexplored. To gain better insights into pollinator-trait dynamics, we investigate the model system of the phlox family (Polemoniaceae), a clade of ∼400 species pollinated by a diversity of vectors. Using a comprehensive phylogeny and large dataset of traits and observations of pollinators, we reconstruct ancestral pollination system, accounting for the temporal history of pollinators. We conduct phylogenetically controlled analyses of trait co-variation and association with pollinators, integrating many analyses over phylogenetic uncertainty. Pollinator shifts are more heterogeneous than previously hypothesized. The evolution of floral traits is partially constrained by phylogenetic history and trait co-variation, but traits are convergent and differences are associated with different pollinators. Trait shifts are usually gradual, rather than rapid, suggesting complex genetic and ecological interactions of flowers at macroevolutionary scales.

The Relationship between Hormones, Glucose, and Oxidative Damage Is Condition and Stress Dependent in a Free-Living Passerine Bird

AbstractPhysiological state is an emergent property of the interactions among physiological systems within an intricate network. Understanding the connections within this network is one of the goals in physiological ecology. Here, we studied the relationship between body condition, two neuroendocrine hormones (corticosterone and insulin-like growth factor 1 [IGF-1]) as physiological regulators, and two physiological systems related to resource metabolism (glucose) and oxidative balance (malondialdehyde). We measured these traits under baseline and stress-induced conditions in free-living house sparrows (Passer domesticus). We used path analysis to analyze different scenarios about the structure of the physiological network. Our data were most consistent with a model in which corticosterone was the major regulator under baseline conditions. This model shows that individuals in better condition have lower corticosterone levels; corticosterone and IGF-1 levels are positively associated; and oxidative damage is higher when levels of corticosterone, IGF-1, and glucose are elevated. After exposure to acute stress, these relationships were considerably reorganized. In response to acute stress, birds increased their corticosterone and glucose levels and decreased their IGF-1 levels. However, individuals in better condition increased their corticosterone levels more and better maintained their IGF-1 levels in response to acute stress. The acute stress-induced changes in corticosterone and IGF-1 levels were associated with an increase in glucose levels, which in turn was associated with a decrease in oxidative damage. We urge ecophysiologists to focus more on physiological networks, as the relationships between physiological traits are complex and dynamic during the organismal stress response.

Weak genetic signal for phenotypic integration implicates developmental processes as major regulators of trait covariation

Phenotypic integration is an important metric that describes the degree of covariation among traits in a population, and is hypothesized to arise due to selection for shared functional processes. Our ability to identify the genetic and/or developmental underpinnings of integration is marred by temporally overlapping cell-, tissue- and structure-level processes that serve to continually 'overwrite' the structure of covariation among traits through ontogeny. Here, we examine whether traits that are integrated at the phenotypic level also exhibit a shared genetic basis (e.g. pleiotropy). We micro-CT scanned two hard tissue traits, and two soft tissue traits (mandible, pectoral girdle, atrium and ventricle, respectively) from an F5 hybrid population of Lake Malawi cichlids, and used geometric morphometrics to extract 3D shape information from each trait. Given the large degree of asymmetric variation that may reflect developmental instability, we separated symmetric from asymmetric components of shape variation. We then performed quantitative trait loci (QTL) analysis to determine the degree of genetic overlap between shapes. While we found ubiquitous associations among traits at the phenotypic level, except for a handful of notable exceptions, our QTL analysis revealed few overlapping genetic regions. Taken together, this indicates developmental interactions can play a large role in determining the degree of phenotypic integration among traits, and likely obfuscate the genotype to phenotype map, limiting our ability to gain a comprehensive picture of the genetic contributors responsible for phenotypic divergence.

Context-dependent trait covariances: how plasticity shapes behavioral syndromes

The study of behavioral syndromes aims to understand among-individual correlations of behavior, yielding insights into the ecological factors and proximate constraints that shape behavior. In parallel, interest has been growing in behavioral plasticity, with results commonly showing that animals vary in their behavioral response to environmental change. These two phenomena are inextricably linked—behavioral syndromes describe cross-trait or cross-context correlations, while variation in behavioral plasticity describes variation in response to changing context. However, they are often discussed separately, with plasticity analyses typically considering a single trait (univariate) across environments, while behavioral trait correlations are studied as multiple traits (multivariate) under one environmental context. Here, we argue that such separation represents a missed opportunity to integrate these concepts. Through observations of multiple traits while manipulating environmental conditions, we can quantify how the environment shapes behavioral correlations, thus quantifying how phenotypes are differentially constrained or integrated under different environmental conditions. Two analytical options exist which enable us to evaluate the context dependence of behavioral syndromes—multivariate reaction norms and character state models. These models are largely two sides of the same coin, but through careful interpretation we can use either to shift our focus to test how the contextual environment shapes trait covariances.

Genetic and environmental canalization are not associated among altitudinally varying populations of Drosophila melanogaster

Organisms are exposed to environmental and mutational effects influencing both mean and variance of phenotypes. Potentially deleterious effects arising from this variation can be reduced by the evolution of buffering (canalizing) mechanisms, ultimately reducing phenotypic variability. There has been interest regarding the conditions enabling the evolution of canalization. Under some models, the circumstances under which genetic canalization evolves are limited despite apparent empirical evidence for it. It has been argued that genetic canalization evolves as a correlated response to environmental canalization (congruence model). Yet, empirical evidence has not consistently supported predictions of a correlation between genetic and environmental canalization. In a recent study, a population of Drosophila adapted to high altitude showed evidence of genetic decanalization relative to those from low altitudes. Using strains derived from these populations, we tested if they varied for multiple aspects of environmental canalization We observed the expected differences in wing size, shape, cell (trichome) density and mutational defects between high- and low-altitude populations. However, we observed little evidence for a relationship between measures of environmental canalization with population or with defect frequency. Our results do not support the predicted association between genetic and environmental canalization.

Divergence in Plant Traits and Increased Modularity Underlie Repeated Transitions Between Low and High Elevations in the Andean Genus Leucheria

Understanding why some plant lineages move from one climatic region to another is a mayor goal of evolutionary biology. In the southern Andes plant lineages that have migrated along mountain ranges tracking cold-humid climates coexist with lineages that have shifted repeatedly between warm-arid at low elevations and cold habitats at high elevations. Transitions between habitats might be facilitated by the acquisition of common traits favoring a resource-conservative strategy that copes with drought resulting from either low precipitation or extreme cold. Alternatively, transitions might be accompanied by phenotypic divergence and accelerated evolution of plant traits, which in turn may depend on the level of coordination among them. Reduced integration and evolution of traits in modules are expected to increase evolutionary rates of traits, allowing diversification in contrasting climates. To examine these hypotheses, we conducted a comparative study in the herbaceous genus Leucheria. We reconstructed ancestral habitat states using Maximum Likelihood and a previously published phylogeny. We performed a Phylogenetic Principal Components Analysis on traits, and then we tested the relationship between PC axes, habitat and climate using Phylogenetic Generalized Least Squares (PGLS). Finally, we compared the evolutionary rates of traits, and the levels of modularity among the three main Clades of Leucheria. Our results suggest that the genus originated at high elevations and later repeatedly colonized arid-semiarid shrublands and humid-forest at lower elevations. PGLS analysis suggested that transitions between habitats were accompanied by shifts in plant strategies: cold habitats at high elevations favored the evolution of traits related to a conservative-resource strategy (thicker and dissected leaves, with high mass per area, and high biomass allocation to roots), whereas warm-arid habitats at lower elevations favored traits related to an acquisitive-resource strategy. As expected, we detected higher levels of modularity in the clades that switched repeatedly between habitats, but higher modularity was not associated with accelerated rates of trait evolution.

Relationships between population traits, nonstructural carbohydrates, and elevation in alpine stands of Vaccinium myrtillus

PREMISE

Despite great attention given to the relationship between plant growth and carbon balance in alpine tree species, little is known about shrubs at the treeline. We hypothesized that the pattern of main nonstructural carbohydrates (NSCs) across elevations depends on the interplay between phenotypic trait plasticity, plant-plant interaction, and elevation.

METHODS

We studied the pattern of NSCs (i.e., glucose, fructose, sucrose, and starch) in alpine stands of Vaccinium myrtillus (above treeline) across an elevational gradient. In the same plots, we measured key growth traits (i.e., anatomical stem features) and shrub cover, evaluating putative relationships with NSCs.

RESULTS

Glucose content was positively related with altitude, but negatively related with shrub cover. Sucrose decreased at high altitude and in older populations and increased with higher percentage of vascular tissue. Starch content increased at middle and high elevations and in stands with high shrub cover. Moreover, starch content was negatively related with the number of xylem rings and the percentage of phloem tissue, but positively correlated with the percentage of xylem tissue.

CONCLUSIONS

We found that the increase in carbon reserves across elevations was uncoupled from plant growth, supporting the growth limitation hypothesis, which postulates NSCs accumulate at high elevation as a consequence of low temperature. Moreover, the response of NSC content to the environmental stress caused by elevation was buffered by phenotypic plasticity of plant traits, suggesting that, under climate warming conditions, shrub expansion due to enhanced plant growth would be pronounced in old but sparse stands.

Variation in mouse pelvic morphology maps to locations enriched in Sox9 Class II and Pitx1 regulatory features

Variation in pelvic morphology has a complex genetic basis and its patterning and specification is governed by conserved developmental pathways. Whether the mechanisms underlying the differentiation and specification of the pelvis also produce the morphological covariation on which natural selection may act, is still an open question in evolutionary developmental biology. We use high-resolution quantitative trait locus (QTL) mapping in the F₃₄ generation of an advanced intercross experiment (LG,SM-G₃₄ ) to characterize the genetic architecture of the mouse pelvis. We test the prediction that genomic features linked to developmental patterning and differentiation of the hind limb and pelvis and the regulation of chondrogenesis are overrepresented in QTL. We find 31 single QTL trait associations at the genome- or chromosome-wise significance level coalescing to 27 pleiotropic loci. We recover further QTL at a more relaxed significance threshold replicating locations found in a previous experiment in an earlier generation of the same population. QTL were more likely than chance to harbor Pitx1 and Sox9 Class II chromatin immunoprecipitation-seq features active during development of skeletal features. There was weak or no support for the enrichment of seven more categories of developmental features drawn from the literature. Our results suggest that genotypic variation is channeled through a subset of developmental processes involved in the generation of phenotypic variation in the pelvis. This finding indicates that the evolvability of complex traits may be subject to biases not evident from patterns of covariance among morphological features or developmental patterning when either is considered in isolation.

Enhanced locomotor performance on familiar surfaces is uncoupled from morphological plasticity in Anolis lizards

The radiation of Anolis lizards in the Caribbean is associated with a diversification of the functional match between morphology, habitat use, and locomotor performance. It has been hypothesized that the microhabitat a lizard is reared in can achieve a similar fit of form and function within a species. This predicts that plasticity in the locomotor apparatus is accompanied by changes in perching behavior or improved locomotor performance. To test this, we raised juveniles of two species (Anolis sagrei and Anolis carolinensis) on either broad or narrow surfaces and examined perching behavior and locomotor performance as well as the shape of the pectoral and pelvic girdles, limb length, and thickness of the long bones. Perching behavior was not affected by the habitat surface experienced during ontogeny. However, individuals raised on broad surfaces showed better locomotor performance on broad surfaces, and the magnitude of the effect was as large as the difference between the two species. Both species showed modifications of pectoral and pelvic shape, but only A. carolinensis developed longer limbs on broad surfaces. However, these morphological adjustments induced by physical activity did not explain why lizards raised on broad surfaces performed better. Thus, it appears that early-life experiences can affect both the morphology of the locomotor apparatus and locomotor performance in Anolis lizards, without the two being functionally connected.

Natural selection acting on integrated phenotypes: covariance among functional leaf traits increases plant fitness

Plant functional strategies are usually accomplished through the simultaneous expression of different traits, and hence their correlations should be promoted by natural selection. The adaptive value of correlations among leaf functional traits, however, has not been assessed in natural populations. We estimated intraspecific variation in leaf functional traits related to the primary metabolism and anti-herbivore defence in a population of Turnera velutina. We analysed whether natural selection favoured the expression of individual traits, particular combinations of traits or leaf phenotypic integration. Patterns of covariation among traits were related to water and nitrogen economy, and were similar among genotypes, but the magnitude of their phenotypic integration differed by 10-fold. Although families did not differ in the mean values of leaf functional traits, directional selection favoured low nitrogen content and low chemical defence, high content of chlorophyll, sugar in extrafloral nectar and trichome density. Families with higher phenotypic integration among leaf traits grew faster and produced more flowers. We suggest that the coordinated expression of leaf traits has an adaptive value, probably related to optimisation in the expression of traits related to water conservation and nitrogen acquisition.

Phenotypic integration and modularity drives skull shape divergence in the Arctic fox (Vulpes lagopus) from the Commander Islands

Phenotypic integration and modularity influence morphological disparity and evolvability. However, studies addressing how morphological integration and modularity change for long periods of genetic isolation are scarce. Here, we investigate patterns of phenotypic integration and modularity in the skull of phenotypically and genetically distinct populations of the Artic fox (Vulpes lagopus) from the Commander Islands of the Aleutian belt (i.e. Bering and Mednyi) that were isolated ca 10 000 years by ice-free waters of the Bering sea. We use three-dimensional geometric morphometrics to quantify the strength of modularity and integration from inter-individual variation (static) and from fluctuating asymmetry (random developmental variation) in both island populations compared to the mainland population (i.e. Chukotka) and we investigated how changes in morphological integration and modularity affect disparity and the directionality of trait divergence. Our results indicate a decrease in morphological integration concomitant to an increase in disparity at a developmental level, from mainland to the smallest and farthest population of Mednyi. However, phenotypic integration is higher in both island populations accompanied by a reduction in disparity compared to the population of mainland at a static level. This higher integration may have favoured morphological adaptive changes towards specific feeding behaviours related to the extreme environmental settings of islands. Our study demonstrates how shifts in phenotypic integration and modularity can facilitate phenotypic evolvability at the intraspecific level that may lead to lineage divergence at macroevolutioanry scales.

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.

Asymmetry in genitalia is in sync with lateralized mating behavior but not with the lateralization of other behaviors

Asymmetries in bilateral organisms attract a lot of curiosity given that they are conspicuous departures from the norm. They allow the investigation of the integration at different levels of biological organization. Here we study whether and how behavioral and asymmetrical anatomical traits co-evolved and work together. We ask if asymmetry is determined locally for each trait or at a whole individual level in a species bearing conspicuous asymmetrical genitalia. Asymmetric genitalia evolved in many species; however, in most cases the direction of asymmetry is fixed. Therefore, it has been rarely determined if there is an association between the direction of asymmetry in genitalia and other traits. In onesided livebearer fish of the genus Jenynsia (Cyprinodontiformes, Anablepidae), the anal fin of males is modified into a gonopodium, an intromittent organ that serves to inseminate females. The gonopodium shows a conspicuous asymmetry, with its tip bending either to the left or the right. By surveying 13 natural populations of Jenynsia lineata, we found that both genital morphs are equally common in wild populations. In a series of experiments in a laboratory population, we discovered asymmetry and lateralization for multiple other traits; yet, the degree of integration varied highly among them. Lateralization in exploratory behavior in response to different stimuli was not associated with genital morphology. Interestingly, the direction of genital asymmetry was positively correlated with sidedness of mating preference and the number of neuromasts in the lateral line. This suggests integration of functionally linked asymmetric traits; however, there is no evidence that asymmetry is determined at the whole individual level in our study species.

Life-history evolution under fluctuating density-dependent selection and the adaptive alignment of pace-of-life syndromes

We present a novel perspective on life-history evolution that combines recent theoretical advances in fluctuating density-dependent selection with the notion of pace-of-life syndromes (POLSs) in behavioural ecology. These ideas posit phenotypic co-variation in life-history, physiological, morphological and behavioural traits as a continuum from the highly fecund, short-lived, bold, aggressive and highly dispersive 'fast' types at one end of the POLS to the less fecund, long-lived, cautious, shy, plastic and socially responsive 'slow' types at the other. We propose that such variation in life histories and the associated individual differences in behaviour can be explained through their eco-evolutionary dynamics with population density - a single and ubiquitous selective factor that is present in all biological systems. Contrasting regimes of environmental stochasticity are expected to affect population density in time and space and create differing patterns of fluctuating density-dependent selection, which generates variation in fast versus slow life histories within and among populations. We therefore predict that a major axis of phenotypic co-variation in life-history, physiological, morphological and behavioural traits (i.e. the POLS) should align with these stochastic fluctuations in the multivariate fitness landscape created by variation in density-dependent selection. Phenotypic plasticity and/or genetic (co-)variation oriented along this major POLS axis are thus expected to facilitate rapid and adaptively integrated changes in various aspects of life histories within and among populations and/or species. The fluctuating density-dependent selection POLS framework presented here therefore provides a series of clear testable predictions, the investigation of which should further our fundamental understanding of life-history evolution and thus our ability to predict natural population dynamics.

Detecting canalization and intra-floral modularity in triggerplant (Stylidium) flowers: correlations are only part of the story

BACKGROUND AND AIMS

The Berg hypothesis proposes that specialized-flower traits experience stronger stabilizing selection than non-floral structures and predicts that variation in specialized-flower traits will be mostly uncorrelated with variation in non-floral traits. Similarly, adaptive-accuracy theory predicts lower variation (as a proportion of the mean) in floral traits than in non-floral ones. Both hypotheses can be extended to comparisons between floral traits, where different parts of the flower can be expected to experience different strengths of stabilizing selection, resulting in contrasting patterns of variation. The present study tests these ideas by analysing variation/covariation in those floral traits influencing the location of pollen placement on, and stigma contact with, pollinators ('pollination-mechanics traits', PMTs) in relation to variation/covariation in non-floral traits and floral traits not directly involved in the mechanics of pollination. The prediction was that PMTs are canalized (buffered against genetic and environmental variation) relative to attraction traits, as manifested in lower variances and modular independence.

METHODS

Floral and inflorescence structures of ten species of triggerplants (Stylidium, Stylidiaceae) in south-western Australia were measured; the data were analysed using multivariate and bivariate approaches to detect modular structure of floral and non-floral traits and assess evidence for canalization of PMTs.

KEY RESULTS

Only six of the ten species had PMTs with smaller correlation coefficients than attraction traits, in contrast to the Berg expectation. However, allometric and variance patterns were generally consistent with the predictions of an extended Berg hypothesis and adaptive accuracy. There was modular separation of most floral traits from non-floral traits and clear intra-floral modular structure. PMTs showed lower proportional variation and shallower allometric slopes than pollinator-attraction traits in nine and eight, respectively, of ten species.

CONCLUSIONS

This study demonstrates the value of allometric and variance analyses (in addition to correlation) in assessing the evolutionary significance of floral-trait stability and plasticity.

Covariation and phenotypic integration in chemical communication displays: biosynthetic constraints and eco-evolutionary implications

Chemical communication is ubiquitous. The identification of conserved structural elements in visual and acoustic communication is well established, but comparable information on chemical communication displays (CCDs) is lacking. We assessed the phenotypic integration of CCDs in a meta-analysis to characterize patterns of covariation in CCDs and identified functional or biosynthetically constrained modules. Poorly integrated plant CCDs (i.e. low covariation between scent compounds) support the notion that plants often utilize one or few key compounds to repel antagonists or to attract pollinators and enemies of herbivores. Animal CCDs (mostly insect pheromones) were usually more integrated than those of plants (i.e. stronger covariation), suggesting that animals communicate via fixed proportions among compounds. Both plant and animal CCDs were composed of modules, which are groups of strongly covarying compounds. Biosynthetic similarity of compounds revealed biosynthetic constraints in the covariation patterns of plant CCDs. We provide a novel perspective on chemical communication and a basis for future investigations on structural properties of CCDs. This will facilitate identifying modules and biosynthetic constraints that may affect the outcome of selection and thus provide a predictive framework for evolutionary trajectories of CCDs in plants and animals.

A fly in a tube: Macroevolutionary expectations for integrated phenotypes

Phenotypic integration and modularity are ubiquitous features of complex organisms, describing the magnitude and pattern of relationships among biological traits. A key prediction is that these relationships, reflecting genetic, developmental, and functional interactions, shape evolutionary processes by governing evolvability and constraint. Over the last 60 years, a rich literature of research has quantified patterns of integration and modularity across a variety of clades and systems. Only recently has it become possible to contextualize these findings in a phylogenetic framework to understand how trait integration interacts with evolutionary tempo and mode. Here, we review the state of macroevolutionary studies of integration and modularity, synthesizing empirical and theoretical work into a conceptual framework for predicting the effects of integration on evolutionary rate and disparity: a fly in a tube. While magnitude of integration is expected to influence the potential for phenotypic variation to be produced and maintained, thus defining the shape and size of a tube in morphospace, evolutionary rate, or the speed at which a fly moves around the tube, is not necessarily controlled by trait interactions. Finally, we demonstrate this reduced disparity relative to the Brownian expectation for a given rate of evolution with an empirical example: the avian cranium.

On the role of sex differences for evolution in heterogeneous and changing fitness landscapes: insights from pygmy grasshoppers

Much research has been devoted to study evolution of local adaptations by natural selection, and to explore the roles of neutral processes and developmental plasticity for patterns of diversity among individuals, populations and species. Some aspects, such as evolution of adaptive variation in phenotypic traits in stable environments, and the role of plasticity in predictable changing environments, are well understood. Other aspects, such as the role of sex differences for evolution in spatially heterogeneous and temporally changing environments and dynamic fitness landscapes, remain elusive. An increased understanding of evolution requires that sex differences in development, physiology, morphology, life-history and behaviours are more broadly considered. Studies of selection should take into consideration that the relationships linking phenotypes to fitness may vary not only according to environmental conditions but also differ between males and females. Such opposing selection, sex-by-environment interaction effects of selection and sex-specific developmental plasticity can have consequences for population differentiation, local adaptations and for the dynamics of polymorphisms. Integrating sex differences in analytical frameworks and population comparisons can therefore illuminate neglected evolutionary drivers and reconcile unexpected patterns. Here, I illustrate these issues using empirical examples from over 20 years of research on colour polymorphic Tetrix subulata and Tetrix undulata pygmy grasshoppers, and summarize findings from observational field studies, manipulation experiments, common garden breeding experiments and population genetics studies.This article is part of the theme issue 'Linking local adaptation with the evolution of sex differences'.

Regressions Fit for Purpose: Models of Locust Phase State Must Not Conflate Morphology With Behavior

Phenotypic plasticity often entails coordinated changes in multiple traits. The effects of two alternative environments on multiple phenotypic traits can be analyzed by multivariable binary logistic regression (LR). Locusts are grasshopper species (family Acrididae) with a capacity to transform between two distinct integrated phenotypes or "phases" in response to changes in population density: a solitarious phase, which occurs when densities are low, and a gregarious phase, which arises as a consequence of crowding and can form very large and economically damaging swarms. The two phases differ in behavior, physiology and morphology. A large body of work on the mechanistic basis of behavioral phase transitions has relied on LR models to estimate the probability of behavioral gregariousness from multiple behavioral variables. Mart́ın-Blázquez and Bakkali (2017; [10.1111/eea.12564]10.1111/eea.12564) have recently proposed standardized LR models for estimating an overall "gregariousness level" from a combination of behavioral and, unusually, morphometric variables. Here I develop a detailed argument to demonstrate that the premise of such an overall "gregariousness level" is fundamentally flawed, since locust phase transformations entail a decoupling of behavior and morphology. LR models that combine phenotypic traits with markedly different response times to environmental change are of very limited value for analyses of phase change in locusts, and of environmentally induced phenotypic transitions in general. I furthermore show why behavioral variables should not be adjusted by measures of body size that themselves differ between the two phases. I discuss the models fitted by Mart́ın-Blázquez and Bakkali (2017) to highlight potential pitfalls in statistical methodology that must be avoided when analysing associations between complex phenotypes and alternative environments. Finally, I reject the idea that "standardized models" provide a valid shortcut to estimating phase state across different developmental stages, strains or species. The points addressed here are pertinent to any research on transitions between complex phenotypes and behavioral syndromes.

Toward Understanding the Repeated Occurrence of Associations between Melanin-Based Coloration and Multiple Phenotypes

Melanin is the most widespread pigment in organisms. Melanin-based coloration has been repeatedly observed to be associated with the same traits and in the same direction in different vertebrate and insect species. However, whether any factors that are common to different taxa account for the repeated evolution of melanin-phenotype associations remains unclear. We propose to approach this question from the perspective of convergent and parallel evolution to clarify to what extent different species have evolved the same associations owing to a shared genetic basis and being subjected to similar selective pressures. Our current understanding of the genetic basis of melanin-phenotype associations allows for both convergent and parallel evolution, but this understanding is still limited. Further research is needed to clarify the generality and interdependencies of the different proposed mechanisms (supergenes, pleiotropy based on hormones, or neural crest cells). The general ecological scenarios whereby melanin-based coloration is under selection-protection from ultraviolet radiation, thermoregulation in cold environments, or as a signal of social status-offer a good opportunity to study how melanin-phenotype associations evolve. Reviewing these scenarios shows that some traits associated with melanin-based coloration might be selected together with coloration by also favoring adaptation but that other associated traits might impede adaptation, which may be indicative of genetic constraints. We therefore encourage further research on the relative roles that selection and genetic constraints play in shaping multiple melanin-phenotype associations. Placed into a phylogenetic context, this will help clarify to what extent these associations result from convergent or parallel evolutionary processes and why melanin-phenotype associations are so common across the tree of life.

Phenotypic integration in an extended phenotype: among-individual variation in nest-building traits of the alfalfa leafcutting bee (Megachile rotundata)

Structures such as nests and burrows are an essential component of many organisms' life-cycle and require a complex sequence of behaviours. Because behaviours can vary consistently among individuals and be correlated with one another, we hypothesized that these structures would (1) show evidence of among-individual variation, (2) be organized into distinct functional modules and (3) show evidence of trade-offs among functional modules due to limits on energy budgets. We tested these hypotheses using the alfalfa leafcutting bee, Megachile rotundata, a solitary bee and important crop pollinator. Megachile rotundata constructs complex nests by gathering leaf materials to form a linear series of cells in pre-existing cavities. In this study, we examined variation in the following nest construction traits: reproduction (number of cells per nest and nest length), nest protection (cap length and number of leaves per cap), cell construction (cell size and number of leaves per cell) and cell provisioning (cell mass) from 60 nests. We found a general decline in investment in cell construction and provisioning with each new cell built. In addition, we found evidence for both repeatability and plasticity in cell provisioning with little evidence for trade-offs among traits. Instead, most traits were positively, albeit weakly, correlated (r ~ 0.15), and traits were loosely organized into covarying modules. Our results show that individual differences in nest construction are detectable at a level similar to that of other behavioural traits and that these traits are only weakly integrated. This suggests that nest components are capable of independent evolutionary trajectories.

Phenotypic integration and the evolution of signal repertoires: A case study of treefrog acoustic communication

Animal signals are inherently complex phenotypes with many interacting parts combining to elicit responses from receivers. The pattern of interrelationships between signal components reflects the extent to which each component is expressed, and responds to selection, either in concert with or independently of others. Furthermore, many species have complex repertoires consisting of multiple signal types used in different contexts, and common morphological and physiological constraints may result in interrelationships extending across the multiple signals in species' repertoires. The evolutionary significance of interrelationships between signal traits can be explored within the framework of phenotypic integration, which offers a suite of quantitative techniques to characterize complex phenotypes. In particular, these techniques allow for the assessment of modularity and integration, which describe, respectively, the extent to which sets of traits covary either independently or jointly. Although signal and repertoire complexity are thought to be major drivers of diversification and social evolution, few studies have explicitly measured the phenotypic integration of signals to investigate the evolution of diverse communication systems. We applied methods from phenotypic integration studies to quantify integration in the two primary vocalization types (advertisement and aggressive calls) in the treefrogs Hyla versicolor, Hyla cinerea, and Dendropsophus ebraccatus. We recorded male calls and calculated standardized phenotypic variance-covariance (P) matrices for characteristics within and across call types. We found significant integration across call types, but the strength of integration varied by species and corresponded with the acoustic similarity of the call types within each species. H. versicolor had the most modular advertisement and aggressive calls and the least acoustically similar call types. Additionally, P was robust to changing social competition levels in H. versicolor. Our findings suggest new directions in animal communication research in which the complex relationships among the traits of multiple signals are a key consideration for understanding signal evolution.

Edaphic history over seedling characters predicts integration and plasticity of integration across geologically variable populations of Arabidopsis thaliana

PREMISE OF THE STUDY

Studies on phenotypic plasticity and plasticity of integration have uncovered functionally linked modules of aboveground traits and seedlings of Arabidopsis thaliana, but we lack details about belowground variation in adult plants. Functional modules can be comprised of additional suites of traits that respond to environmental variation. We assessed whether shoot and root responses to nutrient environments in adult A. thaliana were predictable from seedling traits or population-specific geologic soil characteristics at the site of origin.

METHODS

We compared 17 natural accessions from across the native range of A. thaliana using 14-day-old seedlings grown on agar or sand and plants grown to maturity across nutrient treatments in sand. We measured aboveground size, reproduction, timing traits, root length, and root diameter. Edaphic characteristics were obtained from a global-scale dataset and related to field data.

KEY RESULTS

We detected significant among-population variation in root traits of seedlings and adults and in plasticity in aboveground and belowground traits of adult plants. Phenotypic integration of roots and shoots varied by population and environment. Relative integration was greater in roots than in shoots, and integration was predicted by edaphic soil history, particularly organic carbon content, whereas seedling traits did not predict later ontogenetic stages.

CONCLUSIONS

Soil environment of origin has significant effects on phenotypic plasticity in response to nutrients, and on phenotypic integration of root modules and shoot modules. Root traits varied among populations in reproductively mature individuals, indicating potential for adaptive and integrated functional responses of root systems in annuals.

Beyond mean allelic effects: A locus at the major color gene MC1R associates also with differing levels of phenotypic and genetic (co)variance for coloration in barn owls

The mean phenotypic effects of a discovered variant help to predict major aspects of the evolution and inheritance of a phenotype. However, differences in the phenotypic variance associated to distinct genotypes are often overlooked despite being suggestive of processes that largely influence phenotypic evolution, such as interactions between the genotypes with the environment or the genetic background. We present empirical evidence for a mutation at the melanocortin-1-receptor gene, a major vertebrate coloration gene, affecting phenotypic variance in the barn owl, Tyto alba. The white MC1R allele, which associates with whiter plumage coloration, also associates with a pronounced phenotypic and additive genetic variance for distinct color traits. Contrarily, the rufous allele, associated with a rufous coloration, relates to a lower phenotypic and additive genetic variance, suggesting that this allele may be epistatic over other color loci. Variance differences between genotypes entailed differences in the strength of phenotypic and genetic associations between color traits, suggesting that differences in variance also alter the level of integration between traits. This study highlights that addressing variance differences of genotypes in wild populations provides interesting new insights into the evolutionary mechanisms and the genetic architecture underlying the phenotype.

Non-decoupled morphological evolution of the fore- and hindlimb of sabretooth predators

Specialized organisms are useful for exploring the combined effects of selection of functional traits and developmental constraints on patterns of phenotypic integration. Sabretooth predators are one of the most interesting examples of specialization among mammals. Their hypertrophied, sabre-shaped upper canines and their powerfully built forelimbs have been interpreted as adaptations to a highly specialized predatory behaviour. Given that the elongated and laterally compressed canines of sabretooths were more vulnerable to fracture than the shorter canines of conical-tooth cats, it has been long hypothesized that the heavily muscled forelimbs of sabretooths were used for immobilizing prey before developing a quick and precise killing bite. However, the effect of this unique adaptation on the covariation between the fore- and the hindlimb has not been explored in a quantitative fashion. In this paper, we investigate if the specialization of sabretooth predators decoupled the morphological variation of their forelimb with respect to their hindlimb or, in contrast, both limbs vary in the same fashion as in conical-tooth cats, which do not show such extreme adaptations in their forelimb. We use 3D geometric morphometrics and different morphological indices to compare the fore- and hindlimb of conical- and sabretooth predators. Our results indicate that the limb bones of sabretooth predators covary following the same trend of conical-tooth cats. Therefore, we show that the predatory specialization of sabretooth predators did not result in a decoupling of the morphological evolution of their fore- and hindlimbs. The role of developmental constraints and natural selection on this coordinate variation between the fore- and the hindlimb is discussed in the light of this new evidence.

Genetic and environmental integration of the hawkmoth pollination syndrome in Ruellia humilis (Acanthaceae)

BACKGROUND AND AIMS

The serial homology of floral structures has made it difficult to assess the relative contributions of selection and constraint to floral integration. The interpretation of floral integration may also be clouded by the tacit, but largely untested, assumption that genetic and environmental perturbations affect trait correlations in similar ways. In this study, estimates of both the genetic and environmental correlations between components of the hawkmoth pollination syndrome are presented for chasmogamous flowers of Ruellia humilis , including two levels of control for serial homology.

METHODS

A greenhouse population for quantitative genetic analysis was generated by a partial diallel cross between field-collected plants. An average of 634 chasmogamous flowers were measured for each of eight floral traits that contribute to the hawkmoth syndrome. Genetic correlations (across parents) and environmental correlations (across replicate flowers) were estimated by restricted maximum likelihood.

KEY RESULTS

Stigma height, anther height and floral tube length were very tightly integrated in their responses to both genetic and environmental perturbations. The inclusion of floral disc width as a control for serial homology suggests this integration is an adaptive response to correlational selection imposed by pollinators. In contrast, integration of non-homologous traits was low. Furthermore, when comparisons between the dimensions of serially homologous structures were excluded, the genetic and environmental correlation matrices showed little congruence.

CONCLUSIONS

The results suggest that hawkmoths have imposed strong correlational selection on floral traits involved in the deposition and removal of pollen, and that this is a consequence of stabilizing selection on the relative positions of stigmas and anthers in the face of substantial flower size variation. Low integration of other floral traits, and conflicting patterns of genetic and environmental correlations among these traits, suggest weak or no correlational selection within the range of variability expressed within a population.

Using network analysis to study behavioural phenotypes: an example using domestic dogs

Phenotypic integration describes the complex interrelationships between organismal traits, traditionally focusing on morphology. Recently, research has sought to represent behavioural phenotypes as composed of quasi-independent latent traits. Concurrently, psychologists have opposed latent variable interpretations of human behaviour, proposing instead a network perspective envisaging interrelationships between behaviours as emerging from causal dependencies. Network analysis could also be applied to understand integrated behavioural phenotypes in animals. Here, we assimilate this cross-disciplinary progression of ideas by demonstrating the use of network analysis on survey data collected on behavioural and motivational characteristics of police patrol and detection dogs (Canis lupus familiaris). Networks of conditional independence relationships illustrated a number of functional connections between descriptors, which varied between dog types. The most central descriptors denoted desirable characteristics in both patrol and detection dog networks, with 'Playful' being widely correlated and possessing mediating relationships between descriptors. Bootstrap analyses revealed the stability of network results. We discuss the results in relation to previous research on dog personality, and benefits of using network analysis to study behavioural phenotypes. We conclude that a network perspective offers widespread opportunities for advancing the understanding of phenotypic integration in animal behaviour.

Using network analysis to study behavioural phenotypes: an example using domestic dogs

Phenotypic integration describes the complex interrelationships between organismal traits, traditionally focusing on morphology. Recently, research has sought to represent behavioural phenotypes as composed of quasi-independent latent traits. Concurrently, psychologists have opposed latent variable interpretations of human behaviour, proposing instead a network perspective envisaging interrelationships between behaviours as emerging from causal dependencies. Network analysis could also be applied to understand integrated behavioural phenotypes in animals. Here, we assimilate this cross-disciplinary progression of ideas by demonstrating the use of network analysis on survey data collected on behavioural and motivational characteristics of police patrol and detection dogs (Canis lupus familiaris). Networks of conditional independence relationships illustrated a number of functional connections between descriptors, which varied between dog types. The most central descriptors denoted desirable characteristics in both patrol and detection dog networks, with 'Playful' being widely correlated and possessing mediating relationships between descriptors. Bootstrap analyses revealed the stability of network results. We discuss the results in relation to previous research on dog personality, and benefits of using network analysis to study behavioural phenotypes. We conclude that a network perspective offers widespread opportunities for advancing the understanding of phenotypic integration in animal behaviour.

Experimental assessment of factors mediating the naturalization of a globally invasive tree on sandy coastal plains: a case study from Brazil

As all naturalized species are potential invaders, it is important to better understand the determinants of naturalization of alien plants. This study sought to identify traits that enable the alien tree Casuarina equisetifolia to overcome barriers to survival and reproductive and to become naturalized on sandy coastal plains. Restinga vegetation in Brazil was used as a model system to conceptualize and quantify key stressors (high temperature, solar radiation, drought and salinity) which can limit the initial establishment of the plants. Experiments were conducted to evaluate the effects of these environmental factors on seed persistence in the soil (field), germination (laboratory), survival, growth, phenotypic plasticity and phenotypic integration (greenhouse). Results show that the expected viability of the seeds in the soil was 50 months. Seeds germinated in a similar way in constant and alternating temperatures (20-40 °C), except at 40 °C. Low light, and water and salt stresses reduced germination, but seeds recovered germination when stress diminished. Young plants did not tolerate water stress (<2 % of soil moisture) or deep shade. Growth was greater in sunny than in shady conditions. Although a low degree of phenotypic plasticity is important in habitats with multiple stress factors, this species exhibited high germination plasticity, although young plants showed low plasticity. The positive effect of phenotypic integration on plastic expression in the shade shows that in stressful environments traits that show greater phenotypic plasticity values may have significant phenotypic correlations with other characters, which is an important factor in the evolutionary ecology of this invasive species. Long-term seed persistence in the soil, broad germination requirements (temperature and light conditions) and the capacity to survive in a wide range of light intensity favours its naturalization. However, C. equisetifolia did not tolerate water stress and deep shade, which limit its potential to become naturalized on sandy coastal plain.

EMMLi: A maximum likelihood approach to the analysis of modularity

Identification of phenotypic modules, semiautonomous sets of highly correlated traits, can be accomplished through exploratory (e.g., cluster analysis) or confirmatory approaches (e.g., RV coefficient analysis). Although statistically more robust, confirmatory approaches are generally unable to compare across different model structures. For example, RV coefficient analysis finds support for both two- and six-module models for the therian mammalian skull. Here, we present a maximum likelihood approach that takes into account model parameterization. We compare model log-likelihoods of trait correlation matrices using the finite-sample corrected Akaike Information Criterion, allowing for comparison of hypotheses across different model structures. Simulations varying model complexity and within- and between-module contrast demonstrate that this method correctly identifies model structure and parameters across a wide range of conditions. We further analyzed a dataset of 3-D data, consisting of 61 landmarks from 181 macaque (Macaca fuscata) skulls, distributed among five age categories, testing 31 models, including no modularity among the landmarks and various partitions of two, three, six, and eight modules. Our results clearly support a complex six-module model, with separate within- and intermodule correlations. Furthermore, this model was selected for all five age categories, demonstrating that this complex pattern of integration in the macaque skull appears early and is highly conserved throughout postnatal ontogeny. Subsampling analyses demonstrate that this method is robust to relatively low sample sizes, as is commonly encountered in rare or extinct taxa. This new approach allows for the direct comparison of models with different parameterizations, providing an important tool for the analysis of modularity across diverse systems.

The role of pollinators in the evolution of corolla shape variation, disparity and integration in a highly diversified plant family with a conserved floral bauplan

BACKGROUND AND AIMS

Brassicaceae is one of the most diversified families in the angiosperms. However, most species from this family exhibit a very similar floral bauplan. In this study, we explore the Brassicaceae floral morphospace, examining how corolla shape variation (an estimation of developmental robustness), integration and disparity vary among phylogenetically related species. Our aim is to check whether these floral attributes have evolved in this family despite its apparent morphological conservation, and to test the role of pollinators in driving this evolution.

METHODS

Using geometric morphometric tools, we calculated the phenotypic variation, disparity and integration of the corolla shape of 111 Brassicaceae taxa. We subsequently inferred the phylogenetic relationships of these taxa and explored the evolutionary lability of corolla shape. Finally, we sampled the pollinator assemblages of every taxon included in this study, and determined their pollination niches using a modularity algorithm. We explore the relationship between pollination niche and the attributes of corolla shape.

KEY RESULTS

Phylogenetic signal was weak for all corolla shape attributes. All taxa had generalized pollination systems. Nevertheless, they belong to different pollination niches. There were significant differences in corolla shape among pollination niches even after controlling for the phylogenetic relationship of the plant taxa. Corolla shape variation and disparity was significantly higher in those taxa visited mostly by nocturnal moths, indicating that this pollination niche is associated with a lack of developmental robustness. Corolla integration was higher in those taxa visited mostly by hovering long-tongued flies and long-tongued large bees.

CONCLUSIONS

Corolla variation, integration and disparity were evolutionarily labile and evolved very recently in the evolutionary history of the Brassicaceae. These floral attributes were strongly related to the pollination niche. Even in a plant clade having a very generalized pollination system and exhibiting a conserved floral bauplan, pollinators can drive the evolution of important developmental attributes of corolla shape.

Discovering phenotypic causal structure from nonexperimental data

The evolutionary potential of organisms depends on how their parts are structured into a cohesive whole. A major obstacle for empirical studies of phenotypic organization is that observed associations among characters usually confound different causal pathways such as pleiotropic modules, interphenotypic causal relationships and environmental effects. The present article proposes causal search algorithms as a new tool to distinguish these different modes of phenotypic integration. Without assuming an a priori structure, the algorithms seek a class of causal hypotheses consistent with independence relationships holding in observational data. The technique can be applied to discover causal relationships among a set of measured traits and to distinguish genuine selection from spurious correlations. The former application is illustrated with a biological data set of rat morphological measurements previously analysed by Cheverud et al. (Evolution 1983, 37, 895).

A comparison of floral integration between selfing and outcrossing species: a meta-analysis

BACKGROUND AND AIMS

Floral integration is thought to be an adaptation to promote cross-fertilization, and it is often assumed that it increases morphological matching between flowers and pollinators, increasing the efficiency of pollen transfer. However, the evidence for this role of floral integration is limited, and recent studies have suggested a possible positive association between floral integration and selfing. Although a number of explanations exist to account for this inconsistency, to date there has been no attempt to examine the existence of an association between floral integration and mating system. This study hypothesized that if pollinator-mediated pollen movement among plants (outcrossing) is the main factor promoting floral integration, species with a predominantly outcrossing mating system should present higher levels of floral integration than those with a predominantly selfing mating system.

METHODS

A phylogenetically informed meta-analysis of published data was performed in order to evaluate whether mating system (outcrossing vs. selfing) accounts for the variation in floral integration among 64 species of flowering plants. Morphometric floral information was used to compare intra-floral integration among traits describing sexual organs (androecium and gynoecium) and those corresponding to the perianth (calix and corolla).

KEY RESULTS

The analysis showed that outcrossing species have lower floral integration than selfing species. This pattern was caused by significantly higher integration of sexual traits than perianth traits, as integration of the latter group remained unchanged across mating categories.

CONCLUSIONS

The results suggest that the evolution of selfing is associated with concomitant changes in intra-floral integration. Thus, floral integration of sexual traits should be considered as a critical component of the selfing syndrome.

Shifts and disruptions in resource-use trait syndromes during the evolution of herbaceous crops

Trait-based ecology predicts that evolution in high-resource agricultural environments should select for suites of traits that enable fast resource acquisition and rapid canopy closure. However, crop breeding targets specific agronomic attributes rather than broad trait syndromes. Breeding for specific traits, together with evolution in high-resource environments, might lead to reduced phenotypic integration, according to predictions from the ecological literature. We provide the first comprehensive test of these hypotheses, based on a trait-screening programme of 30 herbaceous crops and their wild progenitors. During crop evolution plants became larger, which enabled them to compete more effectively for light, but they had poorly integrated phenotypes. In a subset of six herbaceous crop species investigated in greater depth, competitiveness for light increased during early plant domestication, whereas diminished phenotypic integration occurred later during crop improvement. Mass-specific leaf and root traits relevant to resource-use strategies (e.g. specific leaf area or tissue density of fine roots) changed during crop evolution, but in diverse and contrasting directions and magnitudes, depending on the crop species. Reductions in phenotypic integration and overinvestment in traits involved in competition for light may affect the chances of upgrading modern herbaceous crops to face current climatic and food security challenges.