Integrated phenotypes: understanding trait covariation in plants and animals

Functional trait response to aridity based on leaf trait network analysis in the Hexi Corridor, China

BACKGROUND

Aridity acts as a strong environmental filter for plants and is predicted to intensify in the future, resulting in changes to leaf functional traits. However, few studies explore how interactions of multiple traits result in leaf trait tradeoff strategies along an aridity gradient, and whether trait separation occurs with increasing aridity intensity. This study examines the impact of long-term aridity on 14 plant leaf traits in two arid areas (arid and hyper-arid) in the Hexi Corridor, China. A leaf trait network (LTN) was constructed to study how leaf trait tradeoff strategies differ between the two areas. Structural equation modeling (SEM) was used to identify the direct and indirect effects of aridity and functional diversity (as measured by community weighted means and functional dispersion) on leaf nutrient concentration.

RESULTS

LTN shows trait separation, poor synergy among traits, and low resource utilization. Correlation analyses showed that the mass ratio hypothesis is dominant, and aridity is positively correlated with leaf relative water content (RWC) and leaf phosphorus content, and negatively correlated with leaf nitrogen content (LNC). SEM results indicated that LNC is directly affected by aridity, RWC, leaf carbon content, and plant height. Aridity and functional dispersion directly affects leaf phosphorus content.

CONCLUSIONS

Results indicate that increasing drought weakens plant coordination among specific traits, and the main change in plant trait tradeoff strategies is reflected in the separation of nutrient traits. Exploring the change of the tradeoff among traits along the aridity gradient can better understand the adaptation process of plants to aridity and the process of community function change.

Climate: The dominant factor influencing the spatial distribution pattern of the leaf trait network of Populus euphratica along the main stream of the Tarim River

Leaves are the primary interface through which plants interact with the environment, their functional traits (morphology, anatomy, physiology) directly reflecting ecological strategies that mediate species-environment interactions. These traits link plant performance to ecosystem processes, shaping species distributions and coexistence via their complex relationships with climatic and edaphic factors. Based on previous work, we selected 20 P. euphratica trees along the desert riparian forest of the main stream of the Tarim River for leaf sample collection and habitat survey. We used 27 leaf traits of P. euphratica to visualize the leaf trait network (LTN). Through network structure parameters, such as edge density, diameter, average path length, and average clustering coefficient, the spatial pattern of the LTN and its relationships with 19 climatic factors and 11 soil factors were discussed using principal component analysis and correlation analysis. The results showed that: (1) there were significant differences in the parameters of the leaf trait network of P. euphratica along the main stream of the Tarim River. The variation coefficients of the diameter and average path length were the largest, respectively, whereas that of the average clustering coefficient was the smallest. (2) Among the parameters, only Modularity was significantly correlated with STK and SOM. But Average clustering coefficient was significantly positively correlated with Isothermality, Average path length and Diameter were significantly positively correlated with Min temperature of coldest month and the Average clustering coefficient was significantly negatively correlated with Min temperature of coldest month, the Diameter was significantly positively correlated with Precipitation of wettest month. In general, the correlation between climate factors and LTNs was stronger than soil factors. (3) The explanatory power of climatic factors alone on the leaf traits of P. euphratica was generally higher than that of soil factors, indicating that climatic conditions play a more decisive role in shaping the network structure of leaf traits of P. euphratica. However, the influence of soil conditions on some LTNs parameters cannot be ignored.The spatial variability of leaf trait networks is driven by climate and soil factors, with climate dominating along the Tarim River's main course.

Macroevolution along developmental lines of least resistance in fly wings

Evolutionary change requires genetic variation, and a reigning paradigm in biology is that rates of microevolution can be predicted from estimates of available genetic variation within populations. However, the accuracy of such predictions should decay on longer evolutionary timescales, as the influence of genetic constraints diminishes. Here we show that intrinsic developmental variability and standing genetic variation in wing shape in two distantly related flies, Drosophila melanogaster and Sepsis punctum, are aligned and predict deep divergence in the dipteran phylogeny, spanning >900 taxa and 185 million years. This alignment cannot be easily explained by constraint hypotheses unless most of the quantified standing genetic variation is associated with deleterious side effects and is effectively unusable for evolution. However, phenotyping of 71 genetic lines of S. punctum revealed no covariation between wing shape and fitness, lending no support to this hypothesis. We also find little evidence for genetic constraints on the pace of wing shape evolution along the dipteran phylogeny. Instead, correlational selection related to allometric scaling, simultaneously shaping developmental variability and deep divergence in fly wings, emerges as a potential explanation for the observed alignment. This suggests that pervasive natural selection has the potential to shape developmental architectures of some morphological characters such that their intrinsic variability predicts their long-term evolution.

Modularity and integration of the neural arch and vertebral centrum in primates

The vertebral column consists of multiple homologous elements that have specialized within and between taxa and serve important functions in positional support and as protection for the central nervous system. The study of modularity and integration provides new insights into the evolution of complex structures such as the vertebral column. Patterns of modularity and integration may reflect underlying genetic-developmental patterns and facilitate evolution. Previous studies have identified mixed modularity patterns within and between elements across mammals generally, within primates and carnivorans. Here, we assess modularity within and between elements in the complete post-axial vertebral column in four catarrhine taxa: Macaca (n = 96), Hylobates (n = 77), Pan (n = 92), and Homo (n = 151). We use the Covariance Ratio (CR) to estimate r2 and the standardized eigenvalues (SVE) variance for comparative purposes. Our results show that there is general, widespread integration within the catarrhine vertebral column, both within and between elements. Hominoids tend to display greater modularity than do macaques, but these estimates are rarely significant. Clusters of modularity in the mid-cervical and upper thoracic regions may relate to special nervous system structures in these areas, and locomotor behaviors in general may influence patterns of modularity in primates. In particular, we find that size is a pervasive factor affecting integration among vertebral elements, though its effects on specific structures are variable. Our results generally do not agree with those found across mammals or within carnivorans, and future studies should focus on genus-level assessments across a variety of taxa.

Phylogeny does not predict the outcome of heterospecific pollen-pistil interactions in a species-rich alpine plant community

PREMISE

Co-occurring plant species that share generalist pollinators often exchange pollen. This heterospecific pollen transfer (HPT) impacts male and female reproductive success through pollen loss and reductions in seed set, respectively. The resulting fitness cost of HPT imposes selection on reproductive traits (e.g., floral color and shape), yet we currently lack strong predictors for the post-pollination fate of heterospecific pollen, especially within community and phylogenetic contexts.

METHODS

We investigated the fate of heterospecific pollen at three distinct stages of plant reproduction: (1) pollen germination on the stigma, (2) pollen tube growth in the style, and (3) fertilization of ovules. We experimentally crossed 11 naturally co-flowering species in the subalpine meadows of the Colorado Rocky Mountains, across a spectrum of phylogenetic relatedness. Using generalized linear mixed models and generalized linear models, we evaluated the effect of parental species identity and phylogenetic relatedness on pollen tube growth at each reproductive stage.

RESULTS

We found that heterospecific pollen tubes can germinate and grow within pistils at each reproductive stage, even when parental species are >100 My divergent. There was no significant effect of phylogenetic distance on heterospecific pollen success, and no evidence for a mechanism that suspends heterospecific pollen germination or pollen tube growth within heterospecific stigmas or styles.

CONCLUSIONS

Our results show that even in communities where HPT is common, pre-zygotic post-pollination mechanisms do not provide strong barriers to interspecific fertilization. HPT can result in the loss of ovules even between highly diverged plant species.

Intraspecific Trait Variation in Seedlings Reveals Independence Between Leaf and Root Traits but a Lack of an Independent "Collaboration Axis" Belowground

Plant functional traits help determine resource acquisition strategies. Global trends at the interspecific scale suggest independence between leaf and root traits described by three functional dimensions: resource acquisition above- and belowground and degree of mycorrhizal collaboration belowground. However, there are ecological and evolutionary reasons to expect different patterns of variation within species, especially within seedlings-the stage at which most tree mortality occurs. Describing the intraspecific patterns of trait variation in seedlings will improve the understanding of tree populations' ability to cope with environmental change. We ask the following questions: (1) How do traits above- and belowground co-vary within species? (2) How do traits relate to soil nutrients and light conditions? We collected root and leaf traits on 131 seedlings from four naturally occurring woody species across eight sites in a temperate, deciduous broadleaf forest in the USA. We measured traits reflecting resource use strategies-specific leaf area, leaf nitrogen, root nitrogen, and root tissue density-and those defining the collaboration axis-specific root length and root diameter. We measured light conditions for each seedling and soil nitrogen and phosphorus to examine the relationship between traits and abiotic conditions using a novel multivariate regression analysis approach. We found that above- and belowground traits segregated into independent functional axes and that the collaboration axis merged with the belowground resource-acquisition axis. We found limited associations between abiotic factors and traits. Our findings suggest that within species, there might be additional constraints to adjust to soil conditions and therefore impact response to environmental change.

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.

Cumulative and interannual effects of reproduction in eastern grey kangaroos

Reproduction can reduce energy allocation to other life-history traits such as survival and growth. Resource constraints give rise to (co)variation in life-history traits and to heterogeneity in energy acquisition and allocation. At each reproductive opportunity, females face a choice between allocation to current reproduction or to maintenance. Many studies compare reproductive trade-offs between two consecutive years, but few account for the cumulative effects of reproduction over multiple years, a crucial factor in understanding life-history evolution in long-lived iteroparous species. We compared short- (interannual) and long-term (cumulative) reproductive trade-offs with a 14-year capture-mark-recapture study of eastern grey kangaroos, where females can have substantial skeletal growth for several years after maturation. We used a multivariate approach to compare how interannual and multi-annual cumulative reproduction affected growth (n = 378 measurements), mass change (n = 376 measurements) and subsequent reproduction (n = 388 measurements), and to quantify (co)variation between these traits among individuals (n = 107) and years (n = 14). Interannually, young females that reproduced experienced decreased skeletal growth compared to young females that did not reproduce. Reproductive females of all ages experienced reduced mass gain and weaning probability in the following year. The cumulative effects of multiple reproductions included decreased skeletal growth, mass gain and weaning probability in the following year. These effects increased with age and reproductive rate. We found positive trait correlations between mass change, leg growth and subsequent reproduction among individuals and years, though weaker at the cumulative than interannual level. Females experience dynamic interannual and cumulative trade-offs. Our analyses of cumulative costs of reproduction revealed long-term trade-offs as well as cumulative costs that were not apparent when estimating interannual costs. Trait correlations suggested heterogeneity in growth and reproduction among females. Years of increased growth were followed by years of increased reproduction, and years of poor growth were followed by years of poor reproduction. Our exploration of both interannual and cumulative costs of reproduction underscores the need to account for long-term reproductive histories to better understand reproductive trade-offs in long-lived iteroparous species.

Adaptive Radiation Without Independent Stages of Trait Evolution in a Group of Caribbean Anoles

Adaptive radiation involves diversification along multiple trait axes, producing phenotypically diverse, species-rich lineages. Theory generally predicts that multi-trait evolution occurs via a "stages" model, with some traits saturating early in a lineage's history, and others diversifying later. Despite its multidimensional nature, however, we know surprisingly little about how different suites of traits evolve during adaptive radiation. Here, we investigated the rate, pattern, and timing of morphological and physiological evolution in the anole lizard adaptive radiation from the Caribbean island of Hispaniola. Rates and patterns of morphological and physiological diversity are largely unaligned, corresponding to independent selective pressures associated with structural and thermal niches. Cold tolerance evolution reflects parapatric divergence across elevation, rather than niche partitioning within communities. Heat tolerance evolution and the preferred temperature evolve more slowly than cold tolerance, reflecting behavioral buffering, particularly in edge-habitat species (a pattern associated with the Bogert effect). In contrast to the nearby island of Puerto Rico, closely related anoles on Hispaniola do not sympatrically partition thermal niche space. Instead, allopatric and parapatric separation across biogeographic and environmental boundaries serves to keep morphologically similar close relatives apart. The phenotypic diversity of this island's adaptive radiation accumulated largely as a by-product of time, with surprisingly few exceptional pulses of trait evolution. A better understanding of the processes that guide multidimensional trait evolution (and nuance therein) will prove key in determining whether the stages model should be considered a common theme of adaptive radiation.

Root and shoot phenology, architecture, and organ properties: an integrated trait network among 44 herbaceous wetland species

Integrating traits across above- and belowground organs offers comprehensive insights into plant ecology, but their various functions also increase model complexity. This study aimed to illuminate the interspecific pattern of whole-plant trait correlations through a network lens, including a detailed analysis of the root system. Using a network algorithm that allows individual traits to belong to multiple modules, we characterize interrelations among 19 traits, spanning both shoot and root phenology, architecture, morphology, and tissue properties of 44 species, mostly herbaceous monocots from Northern Ontario wetlands, grown in a common garden. The resulting trait network shows three distinct yet partially overlapping modules. Two major trait modules indicate constraints of plant size and form, and resource economics, respectively. These modules highlight the interdependence between shoot size, root architecture and porosity, and a shoot-root coordination in phenology and dry-matter content. A third module depicts leaf biomechanical adaptations specific to wetland graminoids. All three modules overlap on shoot height, suggesting multifaceted constraints of plant stature. In the network, individual-level traits showed significantly higher centrality than tissue-level traits do, demonstrating a hierarchical trait integration. The presented whole-plant, integrated network suggests that trait covariation is essentially function-driven rather than organ-specific.

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.

Diversity of Intraspecific Patterns of Brain Region Size Covariation in Fish

Traits often do not evolve in isolation or vary independently of other traits. Instead, they can be affected by covariation, both within and across species. However, the importance of within-species trait covariation and, critically, the degree to which it varies between species has yet to be thoroughly studied. Brain morphology is a trait of great ecological and behavioral importance, with regions that are hypothesized to vary in size based on behavioral and cognitive demands. Sizes of brain regions have also been shown to covary with each other across various taxa. Here, we test the degree to which covariation in brain region sizes within species has been conserved across 10 teleost fish species. These 10 species span five orders, allowing us to examine how phylogenetic proximity influences similarities in intraspecific trait covariation. Our results showed a trend that similar patterns of brain region size covariation occur in more closely related species. Interestingly, there were certain brain region pairs that showed similar levels of covariation across all species regardless of phylogenetic distance, such as the telencephalon and optic tectum, while others, such as the olfactory bulb and the hypothalamus, varied more independently. Ultimately, the patterns of brain region covariation shown here suggest that evolutionary mechanisms or constraints can act on specific brain regions independently, and that these constraints can change over evolutionary time.

Meristic co-evolution and genomic co-localization of lateral line scales and vertebrae in Central American cichlid fishes

Meristic traits are often treated as distinct phenotypes that can be used to differentiate and delineate recently diverged species. For instance, the number of lateral line scales and vertebrae, two traits that vary substantially among Neotropical Heroine cichlid species, have been previously suggested to co-evolve. These meristic traits could co-evolve due to shared adaptive, developmental, or genetic factors. If they were found to be genetically or developmentally non-independent, this might require a more general re-evaluation of their role in evolutionary or taxonomic studies. We expanded a previous analysis of correlated evolution of meristic traits (lateral line scales and vertebrae counts) in these fishes to include a range of phylogenetic reconstructions as well as the analyses of 13 Nicaraguan Midas cichlid species (Amphilophus spp.). Additionally, we performed qualitative traits locus (QTL) mapping in a F2 laboratory-reared hybrid population from two ecologically divergent Midas cichlid fish species to discover and evaluate whether genomic co-segregation might explain the observed patterns of meristic co-evolution. Meristic values for these traits were found to morphologically differentiate some species of the Midas cichlid adaptive radiation. Our QTL analysis pinpointed several genomic regions associated with divergence in these traits and highlighted the potential for genomic co-segregation of the lateral line and vertebrae numbers on two chromosomes. Further, our phylogenetic comparative analyses consistently recovered a significant positive evolutionary correlation between the counts of lateral line scale and vertebrae numbers in Neotropical cichlids. Hence, the findings of genomic co-segregation could partially explain the co-evolution of these two meristic traits in these species. Continuing to unravel the genetic architecture governing meristic divergence helps to better understand both trait correlations and the utility of meristic traits in taxonomic diagnoses and how traits in phenotypes might be expected to co-evolve.

Variable Craniofacial Shape and Development among Multiple Cave-Adapted Populations of Astyanax mexicanus

Astyanax mexicanus is a freshwater fish species with blind cave morphs and sighted surface morphs. Like other troglodytic species, independently evolved cave-dwelling A. mexicanus populations share several stereotypic phenotypes, including the expansion of certain sensory systems, as well as the loss of eyes and pigmentation. Here, we assess the extent to which there is also parallelism in craniofacial development across cave populations. Since multiple forces may be acting upon variation in the A. mexicanus system, including phylogenetic history, selection, and developmental constraint, several outcomes are possible. For example, eye regression may have triggered a conserved series of compensatory developmental events, in which case we would expect to observe highly similar craniofacial phenotypes across cave populations. Selection for cave-specific foraging may also lead to the evolution of a conserved craniofacial phenotype, especially in regions of the head directly associated with feeding. Alternatively, in the absence of a common axis of selection or strong developmental constraints, craniofacial shape may evolve under neutral processes such as gene flow, drift, and bottlenecking, in which case patterns of variation should reflect the evolutionary history of A. mexicanus. Our results found that cave-adapted populations do share certain anatomical features; however, they generally did not support the hypothesis of a conserved craniofacial phenotype across caves, as nearly every pairwise comparison was statistically significant, with greater effect sizes noted between more distantly related cave populations with little gene flow. A similar pattern was observed for developmental trajectories. We also found that morphological disparity was lower among all three cave populations versus surface fish, suggesting eye loss is not associated with increased variation, which would be consistent with a release of developmental constraint. Instead, this pattern reflects the relatively low genetic diversity within cave populations. Finally, magnitudes of craniofacial integration were found to be similar among all groups, meaning that coordinated development among anatomical units is robust to eye loss in A. mexicanus. We conclude that, in contrast to many conserved phenotypes across cave populations, global craniofacial shape is more variable, and patterns of shape variation are more in line with population structure than developmental architecture or selection.

Sex-specific scaling of leaf phosphorus vs. nitrogen under unequal reproductive requirements in Eurya japonica, a dioecious plant

PREMISE

Previous work searching for sexual dimorphism has largely relied on the comparison of trait mean vectors between sexes in dioecious plants. Whether trait scaling (i.e., the ratio of proportional changes in covarying traits) differs between sexes, along with its functional significance, remains unclear.

METHODS

We measured 10 vegetative traits pertaining to carbon, water, and nutrient economics across 337 individuals (157 males and 180 females) of the diocious species Eurya japonica during the fruiting season in eastern China. Piecewise structural equation modeling was employed to reveal the scaling relationships of multiple interacting traits, and multivariate analysis of (co)variance was conducted to test for intersexual differences.

RESULTS

There was no sexual dimorphism in terms of trait mean vectors across the 10 vegetative traits in E. japonica. Moreover, most relationships for covarying trait pairs (17 out of 19) exhibited common scaling slopes between sexes. However, the scaling slopes for leaf phosphorus (P) vs. nitrogen (N) differed between sexes, with 5.6- and 3.0-fold increases of P coinciding with a 10-fold increase of N in male and female plants, respectively.

CONCLUSIONS

The lower ratio of proportional changes in P vs. N for females likely reflects stronger P limitation for their vegetative growth, as they require greater P investments in fruiting. Therefore, P vs. N scaling can be a key avenue allowing for sex-specific strategic optimization under unequal reproductive requirements. This study uncovers a hidden aspect of secondary sex character in dioecious plants, and highlights the use of trait scaling to understand sex-defined economic strategies.

Many-to-many mapping: A simulation study of how the number of traits and tasks affect the evolution of form and function

Many-to-many mapping of form-to-function posits that multiple morphological and physiological traits affect the performance of multiple tasks in an organism, and that redundancy and multitasking occur simultaneously to shape the evolution of an organism's phenotype. Many-to-many mapping is expected to be ubiquitous in nature, yet little is known about how it influences the evolution of organismal phenotype. The F-matrix is a powerful tool to study these issues because it describes how multiple traits affect multiple tasks. We undertook a simulation study using the F-matrix to test how the number of traits and the number of tasks affect trait integration and evolvability, as well as the relationships among tasks. We found that as the number of traits and/or tasks increases, the relationships between the tasks and the integration between the traits become weaker, and that the evolvability of the traits increases, all resulting in a system that is freer to evolve. We also found that as the number of traits increases, performance tradeoffs tend to become weaker, but only to a point. Our work shows that it is important to consider not only multiple traits, but also the multitude of tasks that those traits carry out when studying form-function relationships. We suggest that evolution of these relationships follows functional lines of least resistance, which are less defined in more complex systems, resulting in a mechanism for diversification.

Increased precipitation alters the effects of nitrogen deposition on soil bacterial and fungal communities in a temperate forest

Anthropogenic nitrogen (N) deposition and increased precipitation are known to alter soil microbial communities. However, the combined effects of elevated N deposition and increased precipitation on soil microbial community dynamics and co-occurrence networks in temperate forests remain elusive. In this study, we conducted a field manipulation experiment by applying N solution and water to the forest canopy to simulate natural N deposition and increased precipitation in a temperate forest. We collected samples in the litter layer, organic soil layer, and mineral soil layer in 2018-2019 after 6-7 years of N and water treatments, and explored how elevated N deposition and increased precipitation regulate soil microbial diversity, community composition, and co-occurrence networks in different soil layers and at different sampling times. We found that the effects of N deposition and increased precipitation on soil microbial communities varied with soil layers and sampling times. Compared to the ambient environment, single canopy N addition (CN) or single canopy water addition (CW) did not affect bacterial Shannon diversity in the mineral soil layer in 2018, but the combined canopy N and water additions (CNW) decreased it in this layer at this time. CN increased fungal OTU richness in the organic and mineral soil layers in 2018; however, CW and CNW did not have an effect on it in the same layer at the same time. CW and CNW, but not CN, significantly affected bacterial and fungal community compositions in the litter layer in 2018 and in the organic soil layer in 2019. In contrast, CN, but not CW or CNW, significantly affected fungal community composition in the litter layer in 2019. CNW exhibited higher complexities of bacterial and fungal co-occurrence networks than CN and the ambient environment, indicating increased precipitation can strengthen the effect of N deposition on the complexity of bacterial and fungal co-occurrence networks. Our findings suggest that increased precipitation alters the effects of atmospheric N deposition on soil bacterial and fungal communities in this temperate forest, depending on soil layer and sampling time. Moreover, both bacterial and fungal community compositions are sensitive to increased precipitation, but the bacterial community composition is more sensitive to N deposition than the fungal community composition in the organic and mineral soil layers in this forest.

Adaptive potential of maritime pine under contrasting environments

BACKGROUND

Predicting the adaptability of forest tree populations under future climates requires a better knowledge of both the adaptive significance and evolvability of measurable key traits. Phenotypic plasticity, standing genetic variation and degree of phenotypic integration shape the actual and future population genetic structure, but empirical estimations in forest tree species are still extremely scarce. We analysed 11 maritime pine populations covering the distribution range of the species (119 families and 8 trees/family, ca. 1300 trees) in a common garden experiment planted at two sites with contrasting productivity. We used plant height as a surrogate of fitness and measured five traits (mean and plasticity of carbon isotope discrimination, specific leaf area, needle biomass, Phenology growth index) related to four different strategies (acquisitive economics, photosynthetic organ size, growth allocation and avoidance of water stress).

RESULTS

Estimated values of additive genetic variation would allow adaptation of the populations to future environmental conditions. Overall phenotypic integration and selection gradients were higher at the high productivity site, while phenotypic integration within populations was higher at the low productivity site. Response to selection was related mainly to photosynthetic organ size and drought-avoidance mechanisms rather than to water use efficiency. Phenotypic plasticity of water use efficiency could be maladaptive, resulting from selection for height growth.

CONCLUSIONS

Contrary to the expectations in a drought tolerant species, our study suggests that variation in traits related to photosynthetic organ size and acquisitive investment of resources drive phenotypic selection across and within maritime pine populations. Both genetic variation and evolvability of key adaptive traits were considerably high, including plasticity of water use efficiency. These characteristics would enable a relatively fast micro-evolution of populations in response to the ongoing climate changes. Moreover, differentiation among populations in the studied traits would increase under the expected more productive future Atlantic conditions.

Evaluating modularity in the hominine skull related to feeding biomechanics

OBJECTIVES

Modular architecture of traits in complex organisms can be important for morphological evolution at micro- and sometimes macroevolutionary scales as it may influence the tempo and direction of changes to groups of traits that are essential for particular functions, including food acquisition and processing. We tested several distinct hypotheses about craniofacial modularity in the hominine skull in relation to feeding biomechanics.

MATERIALS AND METHODS

First, we formulated hypothesized functional modules for craniofacial traits reflecting specific demands of feeding biomechanics (e.g., masseter leverage/gape or tooth crown mechanics) in Homo sapiens, Pan troglodytes, and Gorilla gorilla. Then, the pattern and strength of modular signal was quantified by the covariance ratio coefficient and compared across groups using covariance ratio effect size. Hierarchical clustering analysis was then conducted to examine whether a priori-defined functional modules correspond to empirically recovered clusters.

RESULTS

There was statistical support for most a priori-defined functional modules in the cranium and half of the functional modules in the mandible. Modularity signal was similar in the cranium and mandible, and across the three taxa. Despite a similar strength of modularity, the empirically recovered clusters do not map perfectly onto our priori functional modules, indicating that further work is needed to refine our hypothesized functional modules.

CONCLUSION

The results suggest that modular structure of traits in association with feeding biomechanics were mostly shared with humans and the two African apes. Thus, conserved patterns of functional modularity may have facilitated evolutionary changes to the skull during human evolution.

Spatial and temporal heterogeneity alter the cost of plasticity in Pristionchus pacificus

Phenotypic plasticity, the ability of a single genotype to produce distinct phenotypes under different environmental conditions, has become a leading concept in ecology and evolutionary biology, with the most extreme examples being the formation of alternative phenotypes (polyphenisms). However, several aspects associated with phenotypic plasticity remain controversial, such as the existence of associated costs. While already predicted by some of the pioneers of plasticity research, i.e. Schmalhausen and Bradshaw, experimental and theoretical approaches have provided limited support for the costs of plasticity. In experimental studies, one common restriction is the measurement of all relevant parameters over long time periods. Similarly, theoretical studies rarely use modelling approaches that incorporate specific experimentally-derived fitness parameters. Therefore, the existence of the costs of plasticity remains disputed. Here, we provide an integrative approach to understand the cost of adaptive plasticity and its ecological ramifications, by combining laboratory data from the nematode plasticity model system Pristionchus pacificus with a stage-structured population model. Taking advantage of measurements of two isogenic strains grown on two distinct diets, we illustrate how spatial and temporal heterogeneity with regard to the distribution of resources on a metapopulation can alter the outcome of the competition and alleviate the realized cost of plasticity.

Morphological evolution of bird wings follows a mechanical sensitivity gradient determined by the aerodynamics of flapping flight

The physical principles that govern the function of biological structures also mediate their evolution, but the evolutionary drivers of morphological traits within complex structures can be difficult to predict. Here, we use morphological traits measured from 1096 3-dimensional bird wing scans from 178 species to test the interaction of two frameworks for relating morphology to evolution. We examine whether the evolutionary rate (σ2) and mode is dominated by the modular organization of the wing into handwing and armwing regions, and/or the relationship between trait morphology and functional output (i.e. mechanical sensitivity, driven here by flapping flight aerodynamics). Our results support discretization of the armwing and handwing as morphological modules, but morphological disparity and σ2 varied continuously with the mechanical sensitivity gradient and were not modular. Thus, mechanical sensitivity should be considered an independent and fundamental driver of evolutionary dynamics in biomechanical traits, distinct from morphological modularity.

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.

A data-driven framework to model the organism-environment system

Organisms modify their development and function in response to the environment. At the same time, the environment is modified by the activities of the organism. Despite the ubiquity of such dynamical interactions in nature, it remains challenging to develop models that accurately represent them, and that can be fitted using data. These features are desirable when modeling phenomena such as phenotypic plasticity, to generate quantitative predictions of how the system will respond to environmental signals of different magnitude or at different times, for example, during ontogeny. Here, we explain a modeling framework that represents the organism and environment as a single coupled dynamical system in terms of inputs and outputs. Inputs are external signals, and outputs are measurements of the system in time. The framework uses time-series data of inputs and outputs to fit a nonlinear black-box model that allows to predict how the system will respond to novel input signals. The framework has three key properties: it captures the dynamical nature of the organism-environment system, it can be fitted with data, and it can be applied without detailed knowledge of the system. We study phenotypic plasticity using in silico experiments and demonstrate that the framework predicts the response to novel environmental signals. The framework allows us to model plasticity as a dynamical property that changes in time during ontogeny, reflecting the well-known fact that organisms are more or less plastic at different developmental stages.

Accounting for trait variability and coordination in predictions of drought-induced range shifts in woody plants

Functional traits offer a promising avenue to improve predictions of species range shifts under climate change, which will entail warmer and often drier conditions. Although the conceptual foundation linking traits with plant performance and range shifts appears solid, the predictive ability of individual traits remains generally low. In this review, we address this apparent paradox, emphasizing examples of woody plants and traits associated with drought responses at the species' rear edge. Low predictive ability reflects the fact not only that range dynamics tend to be complex and multifactorial, as well as uncertainty in the identification of relevant traits and limited data availability, but also that trait effects are scale- and context-dependent. The latter results from the complex interactions among traits (e.g. compensatory effects) and between them and the environment (e.g. exposure), which ultimately determine persistence and colonization capacity. To confront this complexity, a more balanced coverage of the main functional dimensions involved (stress tolerance, resource use, regeneration and dispersal) is needed, and modelling approaches must be developed that explicitly account for: trait coordination in a hierarchical context; trait variability in space and time and its relationship with exposure; and the effect of biotic interactions in an ecological community context.

What do molecular laws of life mean for species: absolute restrictions or mere suggestions?

Evolutionary biologists are interested in finding universal patterns of covariation between macroscopic and molecular traits. Knowledge of such laws of life can be essential for understanding the course of evolutionary processes. Molecular parameters are presumably close to fundamental limits set to all organisms by laws of physics and chemistry. Thus, laws of life that include such parameters are hypothesized to be similar at both wide interspecific levels of variation and narrower levels of intraspecific and intraindividual variation in different species. In this Commentary, I discuss examples where the significance or direction of such molecular laws of life can be compared at different levels of biological variation: (1) the membrane pacemaker theory of metabolism, (2) the correlation between variation in metabolic rate and mitochondrial efficiency and (3) the allometric scaling of metabolism. All three examples reveal that covariations within species or individuals that include molecular parameters do not always follow patterns observed between species. I conclude that limits set by molecular laws of life can be circumvented (at least to some degree) by changes in other traits, and thus, they usually do not impose strict limitations on minor within-species evolutionary changes (i.e. microevolution). I also briefly discuss some of the most promising perspectives for future studies on the universality of molecular laws of life.

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.

CYCLOIDEA-like genes control floral symmetry, floral orientation, and nectar guide patterning

Actinomorphic flowers usually orient vertically (relative to the horizon) and possess symmetric nectar guides, while zygomorphic flowers often face horizontally and have asymmetric nectar guides, indicating that floral symmetry, floral orientation, and nectar guide patterning are correlated. The origin of floral zygomorphy is dependent on the dorsoventrally asymmetric expression of CYCLOIDEA (CYC)-like genes. However, how horizontal orientation and asymmetric nectar guides are achieved remains poorly understood. Here, we selected Chirita pumila (Gesneriaceae) as a model plant to explore the molecular bases for these traits. By analyzing gene expression patterns, protein-DNA and protein-protein interactions, and encoded protein functions, we identified multiple roles and functional divergence of 2 CYC-like genes, i.e. CpCYC1 and CpCYC2, in controlling floral symmetry, floral orientation, and nectar guide patterning. CpCYC1 positively regulates its own expression, whereas CpCYC2 does not regulate itself. In addition, CpCYC2 upregulates CpCYC1, while CpCYC1 downregulates CpCYC2. This asymmetric auto-regulation and cross-regulation mechanism might explain the high expression levels of only 1 of these genes. We show that CpCYC1 and CpCYC2 determine asymmetric nectar guide formation, likely by directly repressing the flavonoid synthesis-related gene CpF3'5'H. We further suggest that CYC-like genes play multiple conserved roles in Gesneriaceae. These findings shed light on the repeated origins of zygomorphic flowers in angiosperms.

Embodied cognitive morphogenesis as a route to intelligent systems

The embryological view of development is that coordinated gene expression, cellular physics and migration provides the basis for phenotypic complexity. This stands in contrast with the prevailing view of embodied cognition, which claims that informational feedback between organisms and their environment is key to the emergence of intelligent behaviours. We aim to unite these two perspectives as embodied cognitive morphogenesis, in which morphogenetic symmetry breaking produces specialized organismal subsystems which serve as a substrate for the emergence of autonomous behaviours. As embodied cognitive morphogenesis produces fluctuating phenotypic asymmetry and the emergence of information processing subsystems, we observe three distinct properties: acquisition, generativity and transformation. Using a generic organismal agent, such properties are captured through models such as tensegrity networks, differentiation trees and embodied hypernetworks, providing a means to identify the context of various symmetry-breaking events in developmental time. Related concepts that help us define this phenotype further include concepts such as modularity, homeostasis and 4E (embodied, enactive, embedded and extended) cognition. We conclude by considering these autonomous developmental systems as a process called connectogenesis, connecting various parts of the emerged phenotype into an approach useful for the analysis of organisms and the design of bioinspired computational agents.

Plant trait networks reveal adaptation strategies in the drylands of China

BACKGROUND

Plants accomplish multiple functions by the interrelationships between functional traits. Clarifying the complex relationships between plant traits would enable us to better understand how plants employ different strategies to adapt to the environment. Although increasing attention is being paid to plant traits, few studies focused on the adaptation to aridity through the relationship among multiple traits. We established plant trait networks (PTNs) to explore the interdependence of sixteen plant traits across drylands.

RESULTS

Our results revealed significant differences in PTNs among different plant life-forms and different levels of aridity. Trait relationships for woody plants were weaker, but were more modularized than for herbs. Woody plants were more connected in economic traits, whereas herbs were more connected in structural traits to reduce damage caused by drought. Furthermore, the correlations between traits were tighter with higher edge density in semi-arid than in arid regions, suggesting that resource sharing and trait coordination are more advantageous under low drought conditions. Importantly, our results demonstrated that stem phosphorus concentration (SPC) was a hub trait correlated with other traits across drylands.

CONCLUSIONS

The results demonstrate that plants exhibited adaptations to the arid environment by adjusting trait modules through alternative strategies. PTNs provide a new insight into understanding the adaptation strategies of plants to drought stress based on the interdependence among plant functional traits.

Protein-lipid Association in Lizard Chemical Signals

Chemical communication in terrestrial vertebrates is often built on complex blends, where semiochemical and structural compounds may form an integrated functional unit. In lizards, many species have specialized epidermal glands whose secretions are waxy, homogeneous blends of lipids and proteins, both active in communication. The intimate co-occurrence of such compounds allows us to hypothesize that they should undergo a certain degree of covariation, considering both their semiochemical role and the support-to-lipid function hypothesized for the protein fraction. In order to assess the occurrence and level of protein-lipid covariation, we compared the composition and complexity of the two fractions in the femoral gland secretions of 36 lizard species, combining phylogenetically-informed analysis with tandem mass spectrometry. We found the composition and complexity of the two fractions to be strongly correlated. The composition of the protein fraction was mostly influenced by the relative proportion of cholestanol, provitamin D₃, stigmasterol, and tocopherol, while the complexity of the protein pattern increased with that of lipids. Additionally, two identified proteins (carbonic anhydrase and protein disulfide isomerase) increased their concentration as provitamin D₃ became more abundant. Although our approach does not allow us to decrypt the functional relations between the proteinaceous and lipid components, nor under the semiochemical or structural hypothesis, the finding that the proteins involved in this association were enzymes opens up to new perspectives about protein role: They may confer dynamic properties to the blend, making it able to compensate predictable variation of the environmental conditions. This may expand the view about proteins in the support-to-lipid hypothesis, from being a passive and inert component of the secretions to become an active and dynamic one, thus providing cues for future research.

Pollinator Proboscis Length Plays a Key Role in Floral Integration of Honeysuckle Flowers (Lonicera spp.)

Pollinator-mediated selection is supposed to influence floral integration. However, the potential pathway through which pollinators drive floral integration needs further investigations. We propose that pollinator proboscis length may play a key role in the evolution of floral integration. We first assessed the divergence of floral traits in 11 Lonicera species. Further, we detected the influence of pollinator proboscis length and eight floral traits on floral integration. We then used phylogenetic structural equation models (PSEMs) to illustrate the pathway through which pollinators drive the divergence of floral integration. Results of PCA indicated that species significantly differed in floral traits. Floral integration increased along with corolla tube length, stigma height, lip length, and the main pollinators' proboscis length. PSEMs revealed a potential pathway by which pollinator proboscis length directly selected on corolla tube length and stigma height, while lip length co-varied with stigma height. Compared to species with short corolla tubes, long-tube flowers may experience more intense pollinator-mediated selection due to more specialized pollination systems and thus reduce variation in the floral traits. Along elongation of corolla tube and stigma height, the covariation of other relevant traits might help to maintain pollination success. The direct and indirect pollinator-mediation selection collectively enhances floral integration.

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.

Pleiotropy Structures Plant Height and Seed Weight Scaling in Barley despite Long History of Domestication and Breeding Selection

Size scaling describes the relative growth rates of different body parts of an organism following a positive correlation. Domestication and crop breeding often target the scaling traits in the opposite directions. The genetic mechanism of the size scaling influencing the pattern of size scaling remains unexplored. Here, we revisited a diverse barley (Hordeum vulgare L.) panel with genome-wide single-nucleotide polymorphisms (SNPs) profile and the measurement of their plant height and seed weight to explore the possible genetic mechanisms that may lead to a correlation of the two traits and the influence of domestication and breeding selection on the size scaling. Plant height and seed weight are heritable and remain positively correlated in domesticated barley regardless of growth type and habit. Genomic structural equation modeling systematically evaluated the pleiotropic effect of individual SNP on the plant height and seed weight within a trait correlation network. We discovered seventeen novel SNPs (quantitative trait locus) conferring pleiotropic effect on plant height and seed weight, involving genes with function in diverse traits related to plant growth and development. Linkage disequilibrium decay analysis revealed that a considerable proportion of genetic markers associated with either plant height or seed weight are closely linked in the chromosome. We conclude that pleiotropy and genetic linkage likely form the genetic bases of plant height and seed weight scaling in barley. Our findings contribute to understanding the heritability and genetic basis of size scaling and open a new venue for seeking the underlying mechanism of allometric scaling in plants.

The importance of distinguishing individual differences in 'social impact' versus 'social responsiveness' when quantifying indirect genetic effects on the evolution of social plasticity

Social evolution and the dynamics of social interactions have previously been studied under the frameworks of quantitative genetics and behavioural ecology. In quantitative genetics, indirect genetic effects of social partners on the socially plastic phenotypes of focal individuals typically lack crucial detail already included in treatments of social plasticity in behavioural ecology. Specifically, whilst focal individuals (e.g. receivers) may show variation in their 'responsiveness' to the social environment, individual social partners (e.g. signallers) may have a differential 'impact' on focal phenotypes. Here we propose an integrative framework, that highlights the distinction between responsiveness versus impact in indirect genetic effects for a range of behavioural traits. We describe impact and responsiveness using a reaction norm approach and provide statistical models for the assessment of these effects of focal and social partner identity in different types of social interactions. By providing such a framework, we hope to stimulate future quantitative research investigating the causes and consequences of social interactions on phenotypic evolution.

Towards evolutionary predictions: Current promises and challenges

Evolution has traditionally been a historical and descriptive science, and predicting future evolutionary processes has long been considered impossible. However, evolutionary predictions are increasingly being developed and used in medicine, agriculture, biotechnology and conservation biology. Evolutionary predictions may be used for different purposes, such as to prepare for the future, to try and change the course of evolution or to determine how well we understand evolutionary processes. Similarly, the exact aspect of the evolved population that we want to predict may also differ. For example, we could try to predict which genotype will dominate, the fitness of the population or the extinction probability of a population. In addition, there are many uses of evolutionary predictions that may not always be recognized as such. The main goal of this review is to increase awareness of methods and data in different research fields by showing the breadth of situations in which evolutionary predictions are made. We describe how diverse evolutionary predictions share a common structure described by the predictive scope, time scale and precision. Then, by using examples ranging from SARS-CoV2 and influenza to CRISPR-based gene drives and sustainable product formation in biotechnology, we discuss the methods for predicting evolution, the factors that affect predictability and how predictions can be used to prevent evolution in undesirable directions or to promote beneficial evolution (i.e. evolutionary control). We hope that this review will stimulate collaboration between fields by establishing a common language for evolutionary predictions.

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.

Anatomical Correlates of Cursoriality are Compromised by Body Size and Propensity to Burrow in a Group of Small Mammals (Lagomorpha)

Highly cursorial animals are specialised for fast, sustained running via specific morphological adaptations, notably including changes in limb segment length and mechanical advantage. Members of the order Lagomorpha (hares, rabbits and pikas) vary in cursorial ability; hares are generally highly cursorial, rabbits more frequently saltate, and pikas predominantly trot. Previous investigations of lagomorphs have identified anatomical trends correlated with this ‘cursoriality gradient’, however, the phylogenetic sampling of such investigations has been limited to three American species, namely the American pika (Ochotona princeps), brush rabbit (Sylvilagus bachmani), and black-tailed jackrabbit (Lepus californicus). Here, we expand the phylogenetic sample and body size range by including novel data from Australian samples of the European rabbit (Oryctolagus cuniculus) and European hare (L. europaeus), alongside unpublished data on the Eastern cottontail (S. floridanus). X-ray Computed Tomography and digital landmarking were used to capture proportions within the appendicular skeleton of ~ 40 specimens of each European species. In doubling the number of species studied, we find the previously-identified morphological gradients associated with cursorial behaviour are complicated when evaluated in the larger sample. The relative length and joint velocity of limbs was found to be lower than predicted in European rabbits and hares. Furthermore, we present a novel assessment of morphological integration in the lagomorph appendicular skeleton, finding between-limb covariation patterns that are generally similar to those of other mammals. Broadly, these results suggest cursoriality is only one of many selective forces driving lagomorph skeletal evolution, with variations in body size and fossoriality potentially having measurable impacts.

Presence of morphological integration and modularity of the forcipular apparatus in Lithobius melanops (Chilopoda: Lithobiomorpha: Lithobiidae)

The presence of morphological integration and modularity of the forcipular apparatus, despite its evolutionary significance, has not been analyzed in centipedes. This morphological structure has a crucial role in feeding and defense, thanks to its poisonous part (forcipules), which is important for catching the prey. The aims of our study were: i) to test the hypothesis of modularity of the forcipular apparatus in centipede Lithobius melanops; and ii) to investigate the influence of allometry on overall morphological integration in the aforementioned species using a geometric morphometric approach. The presence of fluctuating asymmetry was obtained by Procrustes ANOVA. Allometry was significant only for the symmetric component of the forcipular apparatus. The modularity hypothesis was not accepted, because the covariance coefficients for symmetric and asymmetric components were lower than 89.5% and 72.1% (respectively) of other RV coefficients obtained by a random contiguous partition of the forcipular apparatus. Results of the present study indicate that allometry does increase the level of morphological integration in the forcipular apparatus. According to our results, the forcipular coxosternite and forcipules could not be considered as separate modules; namely, they probably share similar developmental pathways and function in different forms of behavior and survival in L. melanops.

Evidence for an evo-devo-derived hypothesis on three-dimensional flower shape modularity in a tropical orchid clade

Covarying suites of phenotypic traits, or modules, are increasingly recognized to promote morphological evolution. However, information on how modularity influences flower diversity is rare and lacking for Orchidaceae. Here, we combine high-resolution X-ray computed tomography scanning with three-dimensional geometric morphometrics and phylogenetic comparative methods to test various hypotheses about three-dimensional patterns of flower evolutionary modularity in Malagasy Bulbophyllum orchids and examine rates and modes of module evolution. Based on the four evolutionary modules identified (i.e., sepals, lateral petals, labellum + column-foot, and column-part), our data support the hypothesis that both genetic-developmental and functional adaptive factors shaped evolutionary flower trait covariation in these tropical orchids. In line with "evo-devo" studies, we also find that the labellum evolved independently from the rest of the petal whorl. Finally, we show that modules evolved with different rates, and either in a neutral fashion (only column-part) or under selective constraints, as likely imposed by pollinators. Overall, this study supports current views that modular units can enhance the range and rate of morphological evolution.

Morphological integration of the hominoid postcranium

Previous research has suggested that magnitudes of integration may be distinct in the postcranium of hominoids when compared to other primate species. To test this hypothesis, we estimated and compared magnitudes of integration of eight postcranial bones from three-dimensional surface scans for 57 Hylobates lar, 58 Gorilla gorilla, 60 Pan troglodytes, 60 Homo sapiens, 60 Chlorocebus pygerythrus, and 60 Macaca fascicularis. We tested the hypotheses that 1) magnitudes of integration would be distinct in the postcranium of hominoids compared to cercopithecoids, with the explicit prediction that magnitudes of integration would be lower in hominoids than in cercopithecoids, and 2) girdle elements (scapula, os coxa) would have lower magnitudes of integration across all taxa. Integration was quantified using the integration coefficient of variation from interlandmark distances reflecting anatomical and developmental modules defined according to a priori criteria. A resampling protocol was employed to generate distributions of integration values that were then compared statistically using Mann-Whitney U tests with Bonferroni adjustment. Support for hypothesis 1 was mixed: with the exception of Gorilla, hominoid taxa were less integrated than the cercopithecoids for all anatomical modules. However, Homo, Gorilla, and, to a lesser extent, Pan showed higher integration than Hylobates and the cercopithecoids for homologous limb elements, with magnitudes of integration for both modules being lowest for Hylobates. These results generally support the hypothesis of distinct patterns of magnitudes of integration in the hominoid postcranium. The high integration of Gorilla may be explained by the effects of overall body size. The results supported the predictions of the second hypothesis. Regardless of taxon, the os coxa and scapula were generally the least integrated skeletal elements, while the femur and radius were the most integrated. The lower integration of the girdle elements suggests that the geometric complexities of particular elements may significantly influence study outcomes.

Morphological modularity in the anthropoid axial skeleton

Previous research has found that hominoids have stronger modularity between limb elements than other anthropoids, suggesting that there is less constraint on morphological diversification (e.g., limb proportions) in hominoids in terms of evolutionary independence. However, degrees of modularity in the axial skeleton have not been investigated across a broad range of anthropoid taxa. Thus, it is unknown whether hominoids also have stronger modularity in the axial skeleton than other anthropoids, which has implications for the evolution of diverse torso morphologies in Miocene apes as well as the evolution of novel characteristics in the skull and vertebrae of fossil hominins. In this study, 12 anthropoid genera were sampled to examine degrees of modularity between axial skeletal elements (i.e., cranium, mandible, vertebrae, and sacrum). Covariance ratio coefficients were calculated using variance/covariance matrices of interlandmark distances for each axial skeletal element to evaluate degrees of modularity. The results showed that Alouatta, Hylobates, Gorilla, Pan, and Homo showed generally stronger modularity than other anthropoid taxa when considering all axial skeletal elements. When only considering the vertebral elements (i.e., vertebrae and sacrum), Alouatta, Hylobates, Gorilla, and Pan showed generally stronger modularity than other anthropoid taxa. Humans showed stronger modularity between the skull and vertebrae than other hominoids. Thus, the evolution of novel characteristics in the skull and vertebral column may have been less constrained in fossil hominins due to the dissociation of trait covariation between axial skeletal elements in hominoid ancestors, thus fostering more evolutionary independence between the skull and vertebral column.

Modular regulation of floral traits by a PRE1 homolog in Mimulus verbenaceus: implications for the role of pleiotropy in floral integration

Floral traits often show correlated variation within and among species. For species with fused petals, strong correlations among corolla tube, stamen, and pistil length are particularly prevalent, and these three traits are considered an intra-floral functional module. Pleiotropy has long been implicated in such modular integration of floral traits, but empirical evidence based on actual gene function is scarce. We tested the role of pleiotropy in the expression of intra-floral modularity in the monkeyflower species Mimulus verbenaceus by transgenic manipulation of a homolog of Arabidopsis PRE1. Downregulation of MvPRE1 by RNA interference resulted in simultaneous decreases in the lengths of corolla tube, petal lobe, stamen, and pistil, but little change in calyx and leaf lengths or organ width. Overexpression of MvPRE1 caused increased corolla tube and stamen lengths, with little effect on other floral traits. Our results suggest that genes like MvPRE1 can indeed regulate multiple floral traits in a functional module but meanwhile have little effect on other modules, and that pleiotropic effects of these genes may have played an important role in the evolution of floral integration and intra-floral modularity.

Testing candidate genes linked to corolla shape variation of a pollinator shift in Rhytidophyllum (Gesneriaceae)

Floral adaptations to specific pollinators like corolla shape variation often result in reproductive isolation and thus speciation. But despite their ecological importance, the genetic bases of corolla shape transitions are still poorly understood, especially outside model species. Hence, our goal was to identify candidate genes potentially involved in corolla shape variation between two closely related species of the Rhytidophyllum genus (Gesneriaceae family) from the Antilles with contrasting pollination strategies. Rhytidophyllum rupincola has a tubular corolla and is strictly pollinated by hummingbirds, whereas R. auriculatum has more open flowers and is pollinated by hummingbirds, bats, and insects. We surveyed the literature and used a comparative transcriptome sequence analysis of synonymous and non-synonymous nucleotide substitutions to obtain a list of genes that could explain floral variation between R. auriculatum and R. rupincola. We then tested their association with corolla shape variation using QTL mapping in a F2 hybrid population. Out of 28 genes tested, three were found to be good candidates because of a strong association with corolla shape: RADIALIS, GLOBOSA, and JAGGED. Although the role of these genes in Rhytidophyllum corolla shape variation remains to be confirmed, these findings are a first step towards identifying the genes that have been under selection by pollinators and thus involved in reproductive isolation and speciation in this genus.

Leaf Economic and Hydraulic Traits Signal Disparate Climate Adaptation Patterns in Two Co-Occurring Woodland Eucalypts

With climate change impacting trees worldwide, enhancing adaptation capacity has become an important goal of provenance translocation strategies for forestry, ecological renovation, and biodiversity conservation. Given that not every species can be studied in detail, it is important to understand the extent to which climate adaptation patterns can be generalised across species, in terms of the selective agents and traits involved. We here compare patterns of genetic-based population (co)variation in leaf economic and hydraulic traits, climate–trait associations, and genomic differentiation of two widespread tree species (Eucalyptus pauciflora and E. ovata). We studied 2-year-old trees growing in a common-garden trial established with progeny from populations of both species, pair-sampled from 22 localities across their overlapping native distribution in Tasmania, Australia. Despite originating from the same climatic gradients, the species differed in their levels of population variance and trait covariance, patterns of population variation within each species were uncorrelated, and the species had different climate–trait associations. Further, the pattern of genomic differentiation among populations was uncorrelated between species, and population differentiation in leaf traits was mostly uncorrelated with genomic differentiation. We discuss hypotheses to explain this decoupling of patterns and propose that the choice of seed provenances for climate-based plantings needs to account for multiple dimensions of climate change unless species-specific information is available.

Multi-Stemmed Habit in Trees Contributes Climate Resilience in Tropical Dry Forest

Understanding how environmental adaptations mediate plant and ecosystem responses becomes increasingly important under accelerating global environmental change. Multi-stemmed trees, for example, differ in form and function from single-stemmed trees and may possess physiological advantages that allow for persistence during stressful climatic events such as extended drought. Following the worst drought in Hawaii in a century, we examined patterns of stem abundance and turnover in a Hawaiian lowland dry forest (LDF) and a montane wet forest (MWF) to investigate how multi-stemmed trees might influence site persistence, and how stem abundance and turnover relate to key functional traits. We found stem abundance and multi-stemmed trees to be an important component for climate resilience within the LDF. The LDF had higher relative abundance of multi-stemmed trees, stem abundance, and mean stem abundance compared to a reference MWF. Within the LDF, multi-stemmed trees had higher relative stem abundance (i.e., percent composition of stems to the total number of stems in the LDF) and higher estimated aboveground carbon than single-stemmed trees. Stem abundance varied among species and tree size classes. Stem turnover (i.e., change in stem abundance between five-year censuses) varied among species and tree size classes and species mean stem turnover was correlated with mean species stem abundance per tree. At the plot level, stem abundance per tree is also a predictor of survival, though mortality did not differ between multiple- and single-stemmed trees. Lastly, species with higher mean stem abundance per tree tended to have traits associated with a higher light-saturated photosynthetic rate, suggesting greater productivity in periods with higher water supply. Identifying the traits that allow species and forest communities to persist in dry environments or respond to disturbance is useful for forecasting ecological climate resilience or potential for restoration in tropical dry forests.

Experimental evolution of a pheromone signal

Sexual signals are important in speciation, but understanding their evolution is complex as these signals are often composed of multiple, genetically interdependent components. To understand how signals evolve, we thus need to consider selection responses in multiple components and account for the genetic correlations among components. One intriguing possibility is that selection changes the genetic covariance structure of a multicomponent signal in a way that facilitates a response to selection. However, this hypothesis remains largely untested empirically. In this study, we investigate the evolutionary response of the multicomponent female sex pheromone blend of the moth Heliothis subflexa to 10 generations of artificial selection. We observed a selection response of about three-quarters of a phenotypic standard deviation in the components under selection. Interestingly, other pheromone components that are biochemically and genetically linked to the components under selection did not change. We also found that after the onset of selection, the genetic covariance structure diverged, resulting in the disassociation of components under selection and components not under selection across the first two genetic principle components. Our findings provide rare empirical support for an intriguing mechanism by which a sexual signal can respond to selection without possible constraints from indirect selection responses.

Leaf trait network architecture shifts with species-richness and climate across forests at continental scale

Variation in the architecture of trait networks among ecosystems has been rarely quantified, but can provide high resolution of the contrasting adaptation of the whole phenotype. We constructed leaf trait networks (LTNs) from 35 structural, anatomical and compositional leaf traits for 394 tree species in nine forests from tropical to cold-temperate zones in China. Our analyses supported the hypothesis that LTNs would increase in modular complexity across forests in parallel with species-richness and climatic warmth and moisture, due to reduced phenotypic constraints and greater opportunities for niche differentiation. Additionally, we found that within LTNs, leaf economics traits including leaf thickness would have central importance, acting as hub traits with high connectivity due to their contributions to multiple functions. Across the continent, the greater species richness and trait diversity observed in forests under resource-rich climates enable greater complexity in whole phenotype structure and function as indicated by the trait network architecture.

Flower morphology and plant-bee pollinator interactions are related to stamen dimorphism in Melastomataceae

Approximately 20,000 species of flowering plant offer mainly pollen to their pollinators, generally bees. Stamen dimorphism, a floral trait commonly present in some pollen flowers, is thought to be associated with exclusive pollen provision for highly effective bee pollinators. Notwithstanding, little is known about how stamen dimorphism is related to other floral morphological traits and, consequently, plant-pollinator interactions at the community scale. Here we investigated the relationship between stamen dimorphism and other floral morphological traits, as well as the interactions with pollinators in plants of Melastomataceae. We characterized each plant species as stamen dimorphic or stamen isomorphic according to differences in size and shape between stamen sets. Data on interactions between the plants and their bee pollinators were analysed as quantitative bipartite networks. We found that petal and style size and shape were correlated to stamen dimorphism. Stamen dimorphic species present larger flowers and less variable style shapes than stamen isomorphic species. Furthermore, stamen dimorphism is associated with higher richness of visiting bees, i.e. higher ecological generalization. During the evolutionary history of Melastomataceae, the dependence on pollinators for fruit set has possibly favoured the evolution of larger flowers with dimorphic stamens, which in turn are able to make use of a larger spectrum of pollen-collecting bees, leading to ecological generalization.