Differential phenotypic plasticity of subalpine trees predicts trait integration under climate warming.

Understanding limiting factors of phenotypic plasticity is essential given its critical role in shaping biological adaptation and evolution in changing environments. It has been proposed that the pattern of phenotypic correlation could constrain trait plasticity. However, the interplay between phenotypic plasticity and integration has remained contentious. We experimentally simulated climate warming in juveniles of three subalpine tree species by exposing them to three-year in situ open-top chambers (OTCs), and then measured functional plasticity of 72 eco-physiological traits to evaluate whether phenotypic integration constituted an intrinsic constraint to plasticity. We also tested the relationship between the differences in plasticity and maintenance in trait integration. Phenotypic plasticity was positively associated with integration in deciduous tree species under warming. The difference in the plasticity of two paired traits could predict their integration in different environments, where traits displaying more similar plasticity were more likely to be correlated. Our study showed no indication that phenotypic integration constrained plasticity. More importantly, we demonstrated that differential plasticity between traits might result in a notable reorganization of the trait associations, and that warming commonly induced a tighter phenotype. Our study provides new insights into the interplay between phenotypic plasticity and integration in subalpine trees under climate warming.

Trait divergence and opposite above- and below-ground strategies facilitate moso bamboo invasion into subtropical evergreen broadleaf forest.

Understanding the invasion of moso bamboo (Phyllostachys edulis) into adjacent evergreen broadleaf forest based on functional traits is crucial due to its significant influence on ecosystem processes. However, existing research has primarily focused on above- or below-ground traits in isolation, lacking a comprehensive integration of both. In this study, we conducted a trait-based analysis including 23 leaf traits and 11 root traits in three forest types - bamboo forest, mixed bamboo and broadleaf forest, and evergreen broadleaf forest - to investigate trait differences, phenotypic integration, and above- and below-ground resource strategies in bamboo and broadleaf species. Our findings demonstrated significant differences in leaf and root key traits between bamboo and broadleaf species, strongly supporting the "phenotypic divergence hypothesis". Bamboo exhibited stronger trait correlations compared to broadleaf species, indicating higher phenotypic integration. Above- and below-ground strategies were characterized by trade-offs rather than coordination, resulting in a multi-dimensional trait syndrome. Specifically, a unidimensional leaf economics spectrum revealed that bamboo with higher leaf N concentrations (LNC), P concentrations (LPC), and specific leaf area (SLA) adopted a "fast acquisitive" above-ground strategy, while broadleaf species with thicker leaves employed a "slow conservative" above-ground strategy. A two-dimensional root trait syndrome indicated a "conservation" gradient with bamboo adopting a "slow conservative" below-ground strategy associated with higher root tissue density (RTD), and broadleaf species exhibiting a "fast acquisitive" below-ground strategy linked to higher root N concentrations (RNC) and P concentrations (RPC), and a "collaboration" gradient probably ranging from broadleaf species with a "do-it-yourself" strategy characterized by high specific root length (SRL), to bamboo adopting an "outsourcing" strategy with thicker roots. In conclusion, key trait divergence from coexisting broadleaf species, higher phenotypic integration, and multi-dimensional opposite above- and below-ground resource strategies confer competitive advantages to moso bamboo, shedding light on the mechanistic understanding of its invasion into subtropical evergreen broadleaf forest and providing theoretical guidance for maintaining the stability of subtropical forest ecosystem.

Polyploidy and the evolution of phenotypic integration: Network analysis reveals relationships among anatomy, morphology, and physiology.

Premise: Most traits are polygenic and most genes are pleiotropic, resulting in complex, integrated phenotypes. Polyploidy presents an excellent opportunity to explore the evolution of phenotypic integration as entire genomes are duplicated, allowing for new associations among traits and potentially leading to enhanced or reduced phenotypic integration. Despite the multivariate nature of phenotypic evolution, studies often rely on simplistic bivariate correlations that cannot accurately represent complex phenotypes or data reduction techniques that can obscure specific trait relationships. Methods: We apply network modeling, a common gene co-expression analysis, to the study of phenotypic integration to identify multivariate patterns of phenotypic evolution, including anatomy and morphology (structural) and physiology (functional) traits in response to whole genome duplication in the genus Brassica. Results: We identify four key structural traits that are overrepresented in the evolution of phenotypic integration. Seeding networks with key traits allowed us to identify structure-function relationships not apparent from bivariate analyses. In general, allopolyploids exhibited larger, more robust networks indicative of increased phenotypic integration compared to diploids. Discussion: Phenotypic network analysis may provide important insights into the effects of selection on non-target traits, even when they lack direct correlations with the target traits. Network analysis may allow for more nuanced predictions of both natural and artificial selection.

The Paradox of Predictability Provides a Bridge Between Micro- and Macroevolution.

The relationship between the evolutionary dynamics observed in contemporary populations (microevolution) and evolution on timescales of millions of years (macroevolution) has been a topic of considerable debate. Historically, this debate centers on inconsistencies between microevolutionary processes and macroevolutionary patterns. Here, we characterize a striking exception: emerging evidence indicates that standing variation in contemporary populations and macroevolutionary rates of phenotypic divergence are often positively correlated. This apparent consistency between micro- and macroevolution is paradoxical because it contradicts our previous understanding of phenotypic evolution and is so far unexplained. Here, we explore the prospects for bridging evolutionary timescales through an examination of this "paradox of predictability." We begin by explaining why the divergence-variance correlation is a paradox, followed by data analysis to show that the correlation is a general phenomenon across a broad range of temporal scales, from a few generations to tens of millions of years. Then we review complementary approaches from quantitative-genetics, comparative morphology, evo-devo, and paleontology to argue that they can help to address the paradox from the shared vantage point of recent work on evolvability. In conclusion, we recommend a methodological orientation that combines different kinds of short-term and long-term data using multiple analytical frameworks in an interdisciplinary research program. Such a program will increase our general understanding about how evolution works within and across timescales.

Integration of attractive and defensive phytochemicals is unlikely to constrain chemical diversification in a perennial herb.

Diversification of plant chemical phenotypes is typically associated with spatially and temporally variable plant-insect interactions. Floral scent is often assumed to be the target of pollinator-mediated selection, whereas foliar compounds are considered targets of antagonist-mediated selection. However, floral and vegetative phytochemicals can be biosynthetically linked and may thus evolve as integrated phenotypes. Utilizing a common garden of 28 populations of the perennial herb Arabis alpina (Brassicaceae), we investigated integration within and among floral scent compounds and foliar defense compounds (both volatile compounds and tissue-bound glucosinolates). Within floral scent volatiles, foliar volatile compounds, and glucosinolates, phytochemicals were often positively correlated, and correlations were stronger within these groups than between them. Thus, we found no evidence of integration between compound groups indicating that these are free to evolve independently. Relative to self-compatible populations, self-incompatible populations experienced stronger correlations between floral scent compounds, and a trend toward lower integration between floral scent and foliar volatiles. Our study serves as a rare test of integration of multiple, physiologically related plant traits that each are potential targets of insect-mediated selection. Our results suggest that independent evolutionary forces are likely to diversify different axes of plant chemistry without major constraints.

Context-dependent effects of thermal acclimation on physiological correlates of animal personality in Asiatic toads.

Persistent individual variation in behaviour, or 'personality', is a widespread phenomenon in animals, and understanding the evolution of animal personality is a key task of current biology. Natural selection has been proposed to promote the integration of personality with animal 'intrinsic states', such as metabolic or endocrine traits, and this integration varies with ecological conditions. However, these external ecological modulatory effects have rarely been examined. Here, we investigate the effects of thermal acclimation on between-individual covariations between physiology and behaviour in Asiatic toads (Bufo gargarizans) along an altitudinal gradient. Our results reveal that the thermal modulatory effects on the covariations depend on the altitudinal population. Specifically, at low altitudes, between-individual covariations are highly plastic, with risk-taking behaviour covarying with baseline glucocorticoids (GCs) under warm acclimation, but risk-taking and exploration behaviour covarying with resting metabolic rate (RMR) under cold acclimation. In contrast, between-individual covariations are relatively fixed at high altitudes, with risk-taking behaviour consistently covarying with baseline GCs. Furthermore, at low altitudes, changes in covariations between RMR and personality are associated with adjustment of energy management models. Evidently, animal physiological states that determine or covary with personality can adapt according to the seasonal thermal environment and the thermal evolutionary background of populations. Our findings highlight the importance of a multi-system physiological approach to understand the evolution of animal personality.

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.

Conserved patterns and locomotor-related evolutionary constraints in the hominoid vertebral column.

The evolution of the hominoid lineage is characterized by pervasive homoplasy, notably in regions such as the vertebral column, which plays a central role in body support and locomotion. Few isolated and fewer associated vertebrae are known for most fossil hominoid taxa, but identified specimens indicate potentially high levels of convergence in terms of both form and number. Homoplasy thus complicates attempts to identify the anatomy of the last common ancestor of hominins and other taxa and stymies reconstructions of evolutionary scenarios. One way to clarify the role of homoplasy is by investigating constraints via phenotypic integration, which assesses covariation among traits, shapes evolutionary pathways, and itself evolves in response to selection. We assessed phenotypic integration and evolvability across the subaxial (cervical, thoracic, lumbar, sacral) vertebral column of macaques (n = 96), gibbons (n = 77), chimpanzees (n = 92), and modern humans (n = 151). We found a mid-cervical cluster that may have shifted cranially in hominoids, a persistent thoracic cluster that is most marked in chimpanzees, and an expanded lumbosacral cluster in hominoids that is most expanded in gibbons. Our results highlight the highly conserved nature of the vertebral column. Taxa appear to exploit existing patterns of integration and ontogenetic processes to shift, expand, or reduce cluster boundaries. Gibbons appear to be the most highly derived taxon in our sample, possibly in response to their highly specialized locomotion.

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.

Coevolution with toxic prey produces functional trade-offs in sodium channels of predatory snakes.

Seemingly unrelated traits often share the same underlying molecular mechanisms, potentially generating a pleiotropic relationship whereby selection shaping one trait can simultaneously compromise another. While such functional trade-offs are expected to influence evolutionary outcomes, their actual relevance in nature is masked by obscure links between genotype, phenotype, and fitness. Here, we describe functional trade-offs that likely govern a key adaptation and coevolutionary dynamics in a predator-prey system. Several garter snake (Thamnophis spp.) populations have evolved resistance to tetrodotoxin (TTX), a potent chemical defense in their prey, toxic newts (Taricha spp.). Snakes achieve TTX resistance through mutations occurring at toxin-binding sites in the pore of snake skeletal muscle voltage-gated sodium channels (NaV1.4). We hypothesized that these mutations impair basic NaV functions, producing molecular trade-offs that should ultimately scale up to compromised organismal performance. We investigate biophysical costs in two snake species with unique and independently evolved mutations that confer TTX resistance. We show electrophysiological evidence that skeletal muscle sodium channels encoded by toxin-resistant alleles are functionally compromised. Furthermore, skeletal muscles from snakes with resistance genotypes exhibit reduced mechanical performance. Lastly, modeling the molecular stability of these sodium channel variants partially explains the electrophysiological and muscle impairments. Ultimately, adaptive genetic changes favoring toxin resistance appear to negatively impact sodium channel function, skeletal muscle strength, and organismal performance. These functional trade-offs at the cellular and organ levels appear to underpin locomotor deficits observed in resistant snakes and may explain variation in the population-level success of toxin-resistant alleles across the landscape, ultimately shaping the trajectory of snake-newt coevolution.

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.

Are animal personality, body condition, physiology and structural size integrated? A comparison of species, populations and sexes, and the value of study replication.

A major question in behavioural ecology is why behaviour, physiology and morphology are often integrated into syndromes. In great tits, Parus major, for example, explorative males are larger (vs. smaller) and leaner (vs. heavier) compared to less explorative individuals. Unfortunately, considerable debate exists on whether patterns found in specific studies are replicable. This debate calls for study replication among species, populations and sexes. We measured behavioural (exploration), physiological (breathing rate) and morphological traits (body mass, tarsus length, wing length, bill length) in two species (great vs. blue tits Cyanistes caeruleus), two populations (Forstenrieder Park vs. Starnberg) and two sexes (males vs. females). We then tested whether the same pattern of integration characterized all unique combinations of these three biological categories (hereafter called datasets). We used a multi-year repeated measures set-up to estimate among-individual trait correlation matrices for each dataset. We then used structural equation modelling to test for size-dependent behaviour and physiology, size-corrected (i.e. size-independent) behaviour-physiology correlations and size-corrected body mass-dependent behaviour and physiology. Finally, we used meta-analyses to test which structural paths were generally (vs. conditionally) supported (vs. unsupported). We found general and consistent support for size-dependent physiology and size-corrected body mass-dependent physiology across datasets: faster breathers were smaller but heavier for their size. Unexpectedly, condition-dependent behaviour was not supported: explorative birds were neither leaner, nor was this relationship heterogeneous across datasets. All other hypothesized patterns were dataset-specific: the covariance between size and behaviour, and between behaviour and physiology differed in sign between datasets, and both were, on average, not supported. This heterogeneity was not explained by any of our moderators: species, population or sex. The specific pattern of size- and condition-dependent physiology reported for a unique combination of species, population, and sex, thus generally predicted those in others. Patterns of size- or condition-dependent behaviour (i.e. 'personality'), or behaviour-physiology syndromes reported in specific datasets, by contrast, did not. These findings call for studies revealing the ecological background of this variation and highlight the value of study replication to help understand whether patterns of phenotypic integration reported in one study can be generalized.

Exact expressions and numerical evaluation of average evolvability measures for characterizing and comparing [Formula: see text] matrices.

Theory predicts that the additive genetic covariance ([Formula: see text]) matrix determines a population's short-term (in)ability to respond to directional selection-evolvability in the Hansen-Houle sense-which is typically quantified and compared via certain scalar indices called evolvability measures. Often, interest is in obtaining the averages of these measures across all possible selection gradients, but explicit formulae for most of these average measures have not been known. Previous authors relied either on approximations by the delta method, whose accuracy is generally unknown, or Monte Carlo evaluations (including the random skewers analysis), which necessarily involve random fluctuations. This study presents new, exact expressions for the average conditional evolvability, average autonomy, average respondability, average flexibility, average response difference, and average response correlation, utilizing their mathematical structures as ratios of quadratic forms. The new expressions are infinite series involving top-order zonal and invariant polynomials of matrix arguments, and can be numerically evaluated as their partial sums with, for some measures, known error bounds. Whenever these partial sums numerically converge within reasonable computational time and memory, they will replace the previous approximate methods. In addition, new expressions are derived for the average measures under a general normal distribution for the selection gradient, extending the applicability of these measures into a substantially broader class of selection regimes.

Leaf phenotypic plasticity coupled with integration facilitates the adaptation of plants to enhanced N deposition.

The response of leaf functional traits can provide vital insight into the adaptive strategies of plants under global change. However, empirical knowledge on the acclimation of functional coordination between phenotypic plasticity and integration to increased nitrogen (N) deposition is still scarce. The variation of leaf functional traits of two dominant seedling species, Machilus gamblei and Neolitsea polycarpa, across four N deposition rates (0, 3, 6, and 12 kg N ha-1yr-1), along with the relationship between leaf phenotypic plasticity and integration were investigated in a subtropical montane forest. We found that enhanced N deposition promoted the development of seedling traits toward the direction of resource acquisition, including improved leaf N content, specific leaf area and photosynthetic performance. Appropriate N deposition (≤6 kg N ha-1 yr-1) might induce the optimization of leaf functional traits to promote the capability and efficiency of nutrient use and photosynthesis in seedlings. However, excessive N deposition (12 kg N ha-1 yr-1) would result in detrimental effects on leaf morphological and physiological traits, thus inhibiting the efficiency in resource acquisition. A positive relationship occurred between leaf phenotypic plasticity and integration in both seedling species, implied that higher plasticity of leaf functional traits likely led to better integration with other traits under N deposition. Overall, our study emphasized that leaf functional traits could rapidly respond to changes in N resource, while the coordination between leaf phenotypic plasticity and integration can facilitate the adaptation of tree seedlings in coping with enhanced N deposition. Further studies are still needed on the role of leaf phenotypic plasticity and integration in plant fitness for predicting ecosystem functioning and forest dynamics, especially in the context of future high N deposition.

Frank Beach Award Winner: The centrality of the hypothalamic-pituitary-adrenal axis in dealing with environmental change across temporal scales.

Understanding if and how individuals and populations cope with environmental change is an enduring question in evolutionary ecology that has renewed importance given the pace of change in the Anthropocene. Two evolutionary strategies of coping with environmental change may be particularly important in rapidly changing environments: adaptive phenotypic plasticity and/or bet hedging. Adaptive plasticity could enable individuals to match their phenotypes to the expected environment if there is an accurate cue predicting the selective environment. Diversifying bet hedging involves the production of seemingly random phenotypes in an unpredictable environment, some of which may be adaptive. Here, I review the central role of the hypothalamic-pituitary-adrenal (HPA) axis and glucocorticoids (GCs) in enabling vertebrates to cope with environmental change through adaptive plasticity and bet hedging. I first describe how the HPA axis mediates three types of adaptive plasticity to cope with environmental change (evasion, tolerance, recovery) over short timescales (e.g., 1-3 generations) before discussing how the implications of GCs on phenotype integration may depend upon the timescale under consideration. GCs can promote adaptive phenotypic integration, but their effects on phenotypic co-variation could also limit the dimensions of phenotypic space explored by animals over longer timescales. Finally, I discuss how organismal responses to environmental stressors can act as a bet hedging mechanism and therefore enhance evolvability by increasing genetic or phenotypic variability or reducing patterns of genetic and phenotypic co-variance. Together, this emphasizes the crucial role of the HPA axis in understanding fundamental questions in evolutionary ecology.

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 F5 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.

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.

Univariate and multivariate plasticity in response to incubation temperature in an Australian lizard.

Environments, particularly developmental environments, can generate a considerable amount of phenotypic variation through phenotypic plasticity. Plasticity in response to incubation temperature is well characterised in egg-laying reptiles. However, traits do not always vary independently of one another, and studies encompassing a broad range of traits spanning multiple categories are relatively rare but crucial to better understand whole-organism responses to environmental change, particularly if covariation among traits may constrain plasticity. In this study, we investigated multivariate plasticity in response to incubation across three temperatures in the delicate skink, Lampropholis delicata, and whether this was affected by covariation among traits. At approximately 1 month of age, a suite of growth, locomotor performance, thermal physiology and behavioural traits were measured. Plasticity in the multivariate phenotype of delicate skinks was distinct for different incubation temperatures. Cool temperatures drove shifts in growth, locomotor performance and thermal physiology, while hot temperatures primarily caused changes in locomotor performance and behaviour. These differences are likely due to variation in thermal reaction norms, as there was little evidence that covariation among traits or phenotypic integration influenced plasticity, and there was no effect of incubation temperature on the direction or strength of covariation. While there were broad themes in terms of which trait categories were affected by different incubation treatments, traits appeared to be affected independently by developmental temperature. Comparing reaction norms of a greater range of traits and temperatures will enable better insight into these patterns among trait categories, as well as the impacts of environmental change.