General quantitative genetic methods for comparative biology: phylogenies, taxonomies and multi-trait models for continuous and categorical characters

Metabolic loads and the costs of metazoan reproduction

Reproduction includes two energy investments-the energy in the offspring and the energy expended to make them. The former is well understood, whereas the latter is unquantified but often assumed to be small. Without understanding both investments, the true energy costs of reproduction are unknown. We present a framework for estimating the total energy costs of reproduction by combining data on the energy content of offspring (direct costs) and the metabolic load of bearing them (indirect costs). We find that direct costs typically represent the smaller fraction of the energy expended on reproduction. Mammals pay the highest reproductive costs (excluding lactation), ~90% of which are indirect. Ectotherms expend less on reproduction overall, and live-bearing ectotherms pay higher indirect costs compared with egg-layers. We show that the energy demands of reproduction exceed standard assumptions.

Maternal transfer of per- and polyfluoroalkyl substances (PFAS) in wild birds: A systematic review and meta-analysis

Per- and polyfluoroalkyl substances (PFAS) are synthetic chemicals widely used in consumer products. PFAS can accumulate in animal tissues, resulting in biomagnification and adverse effects on wildlife, such as reproductive impairment. In bird species, PFAS are transferred from mothers to eggs along with essential nutrients and may affect embryo development. However, the extent of maternal PFAS transfer across different species and compounds remains poorly understood. Here, we conducted a systematic review and meta-analysis to quantify maternal PFAS transfer in wild birds and investigate potential sources of variation. We tested the moderating effects of compounds' physicochemical properties and biological traits of studied birds. The dataset included 505 measurements of PFAS concentration and 371 effect sizes derived from 13 studies on 16 bird species and 25 compounds. Overall, across all studies and species, we found a 41% higher concentration of PFAS in offspring than in mothers. Specifically, contaminants were concentrated in the yolk, longer and heavier compounds showed preferential transfer, larger clutch size was associated with decreased PFAS transfer and a higher transfer rate was shown in species with piscivorous and opportunistic/diverse diets. A validation assessment showed good robustness of the overall meta-analytic result. Given the crucial role of birds in maintaining ecological balance, this research article has relevant implications for modelling the impacts of PFAS on wildlife, ecosystems, and human health.

Nest complexity correlates with larger brain size but smaller body mass across bird species

Amniotes differ substantially in absolute and relative brain size after controlling for allometry, and numerous hypotheses have been proposed to explain brain size evolution. Brain size is thought to correlate with processing capacity and the brain's ability to support complex manipulation such as nest-building skills. The increased complexity of nest structure is supposed to be a measure of an ability to manipulate nesting material into the required shape. The degree of nest-structure complexity is also supposed to be associated with body mass, partly because small species lose heat faster and delicate and insulated nests are more crucial for temperature control of eggs during incubation by small birds. Here, we conducted comparative analyses to test these hypotheses by investigating whether the complexity of species-typical nest structure can be explained by brain size and body mass (a covariate also to control for allometric effects on brain size) across 1353 bird species from 147 families. Consistent with these hypotheses, our results revealed that avian brain size increases as the complexity of the nest structure increases after controlling for a significant effect of body size, and also that a negative relationship exists between nest complexity and body mass.

New horizons for comparative studies and meta-analyses

Comparative analyses and meta-analyses are key tools to elucidate broad biological principles, yet the two approaches often appear different in purpose. We propose an integrated approach that can generate deeper insights into ecoevolutionary processes. Marrying comparative and meta-analytic approaches will allow for (i) a more accurate investigation of drivers of biological variation, (ii) a greater ability to account for sources of non-independence in experimental data, (iii) more effective control of publication bias, and (iv) improved transparency and reproducibility. Stronger integration of meta-analytic and comparative studies can also broaden the scope from species-centric investigations to community-level responses and function-valued traits (e.g., reaction norms). We illuminate commonalities, differences, and the transformative potential of combining these methodologies for advancing ecology and evolutionary biology.

Small skeletons show size-specific scaling: an exploration of allometry in the mammalian lumbar spine

Studies of vertebrate bone biomechanics often focus on skeletal adaptations at upper extremes of body mass, disregarding the importance of skeletal adaptations at lower extremes. Yet mammals are ancestrally small and most modern species have masses under 5 kg, so the evolution of morphology and function at small size should be prioritized for understanding how mammals subsist. We examined allometric scaling of lumbar vertebrae in the small-bodied Philippine endemic rodents known as cloud rats, which vary in mass across two orders of magnitude (15.5 g-2700 g). External vertebral dimensions scale with isometry or positive allometry, likely relating to body size and nuances in quadrupedal posture. In contrast to most mammalian trabecular bone studies, bone volume fraction and trabecular thickness scale with positive allometry and isometry, respectively. It is physiologically impossible for these trends to continue to the upper extremes of mammalian body size, and we demonstrate a fundamental difference in trabecular bone allometry between large- and small-bodied mammals. These findings have important implications for the biomechanical capabilities of mammalian bone at small body size; for the selective pressures that govern skeletal evolution in small mammals; and for the way we define 'small' and 'large' in the context of vertebrate skeletons.

Significant shifts in latitudinal optima of North American birds

Changes in climate can alter environmental conditions faster than most species can adapt. A prediction under a warming climate is that species will shift their distributions poleward through time. While many studies focus on range shifts, latitudinal shifts in species' optima can occur without detectable changes in their range. We quantified shifts in latitudinal optima for 209 North American bird species over the last 55 y. The latitudinal optimum (m) for each species in each year was estimated using a bespoke flexible non-linear zero-inflated model of abundance vs. latitude, and the annual shift in m through time was quantified. One-third (70) of the bird species showed a significant shift in their optimum. Overall, mean peak abundances of North American birds have shifted northward, on average, at a rate of 1.5 km per year (±0.58 SE), corresponding to a total distance moved of 82.5 km (±31.9 SE) over the last 55 y. Stronger poleward shifts at the continental scale were linked to key species' traits, including thermal optimum, habitat specialization, and territoriality. Shifts in the western region were larger and less variable than in the eastern region, and they were linked to species' thermal optimum, habitat density preference, and habitat specialization. Individual species' latitudinal shifts were most strongly linked to their estimated thermal optimum, clearly indicating a climate-driven response. Displacement of species from their historically optimal realized niches can have dramatic ecological consequences. Effective conservation must consider within-range abundance shifts. Areas currently deemed "optimal" are unlikely to remain so.

Comparing Bayesian models for the genetic evaluation of oocytes and embryo counts in Dairy Gir cattle

Count traits are usually explored in livestock breeding programs, and they usually do not fit into normal distribution, requiring alternatives to adjust the phenotype to estimate accurate genetic parameters and breeding values. Alternatively, distribution such as Poisson can be used to evaluate count traits. This study aimed to compare and discuss the genetic evaluation for oocyte and embryo counts considering Gaussian (untransformed variable - LIN; transformed by logarithm - LOG; transformed by Anscombe - ANS) and Poisson (POI) distributions. The data comprised 11,343 total oocytes (TO), viable oocytes (VO), cleaved embryos (CE), and viable embryo (VE) records of ovum pick-up from 1740 Dairy Gir heifers and cows. The genetic parameters and breeding values were estimated by the MCMCglmm package of the R software. The posterior means of heritability varied from 0.40 (LIN) to 0.49 (POI) for TO, 0.39 (LIN) to 0.49 (POI) for VO, 0.30 (LOG) to 0.41 (POI) for cleaved embryos, and 0.19 (LIN) to 0.32 (POI) for viable embryos. The posterior means of repeatability varied from 0.56 (LIN) to 0.65 (POI) for TO, 0.53 (LOG) to 0.63 (POI) for VO, 0.44 (LOG) to 0.60 (POI) for CE, and 0.36 (LOG) to 0.56 (POI) for VE. Deviance information criterion and mean squared residuals indicated that POI model should be used for the genetic evaluation of embryo and oocyte count traits. Spearman's rank correlation between estimated breeding value (EBV) for embryo and oocyte count traits computed by POI, LOG, and ANS models was high (ranging from 0.77 to 0.99), indicating little reranking among the best animals. The POI model is the most adequate for genetic evaluation, resulting in more reliable EBV of oocyte and embryo count traits for Dairy Gir cattle.

The fitness consequences of wildlife conservation translocations: a meta-analysis

Conservation translocation is a common strategy to offset mounting rates of population declines through the transfer of captive- or wild-origin organisms into areas where conspecific populations are imperilled or completely extirpated. Translocations that supplement existing populations are referred to as reinforcements and can be conducted using captive-origin animals [ex situ reinforcement (ESR)] or wild-origin animals without any captive ancestry [in situ reinforcement (ISR)]. These programs have been criticized for low success rates and husbandry practices that produce individuals with genetic and performance deficits, but the post-release performance of captive-origin or wild-origin translocated groups has not been systematically reviewed to quantify success relative to wild-resident control groups. To assess the disparity in post-release performance of translocated organisms relative to wild-resident conspecifics and examine the association of performance disparity with organismal and methodological factors across studies, we conducted a systematic review and meta-analysis of 821 performance comparisons from 171 studies representing nine animal classes (101 species). We found that translocated organisms have 64% decreased odds of out-performing their wild-resident counterparts, supporting claims of systemic issues hampering conservation translocations. To help identify translocation practices that could maximize program success in the future, we further quantified the impact of broad organismal and methodological factors on the disparity between translocated and wild-resident conspecific performance. Pre-release animal enrichment significantly reduced performance disparities, whereas our results suggest no overall effects of taxonomic group, sex, captive generation time, or the type of fitness surrogate measured. This work is the most comprehensive systematic review to date of animal conservation translocations in which wild conspecifics were used as comparators, thereby facilitating an evaluation of the overall impact of this conservation strategy and identifying specific actions to increase success. Our review highlights the need for conservation managers to include both sympatric and allopatric wild-reference groups to ensure the post-release performance of translocated animals can be evaluated. Further, our analyses identify pre-release animal enrichment as a particular strategy for improving the outcomes of animal conservation translocations, and demonstrate how meta-analysis can be used to identify implementation choices that maximize translocated animal contributions to recipient population growth and viability.

Adaptive responses to living in stressful habitats: Do invasive and native plant populations use different strategies?

Plants inhabit stressful environments characterized by a variety of stressors, including mine sites, mountains, deserts, and high latitudes. Populations from stressful and reference (non-stressful) sites often have performance differences. However, while invasive and native species may respond differently to stressful environments, there is limited understanding of the patterns in reaction norms of populations from these sites. Here, we use phylogenetically controlled meta-analysis to assess the performance of populations under stress and non-stress conditions. We ask whether stress populations of natives and invasives differ in the magnitude of lowered performance under non-stress conditions and if they vary in the degree of performance advantage under stress. We also assessed whether these distinctions differ with stress intensity. Our findings revealed that natives not only have greater adaptive advantages but also more performance reductions than invasives. Populations from very stressful sites had more efficient adaptations, and performance costs increased with stress intensity in natives only. Overall, the results support the notion that adaptation is frequently costless. Reproductive output was most closely associated with adaptive costs and benefits. Our study characterized the adaptive strategies used by invasive and native plants under stressful conditions, thereby providing important insights into the limitations of adaptation to extreme sites.

Dyadic inter-group cooperation in shotgun hunting activities in a Congo Basin village

Understanding the dynamics of inter-group cooperation in human adaptation has been the subject of recent empirical and theoretical studies in evolutionary anthropology, beginning to fill gaps in our knowledge of how interactions across political, economic and social domains can - and often do - lead to stable, large-scale cooperation. Here we investigate dyadic intergroup cooperation in shotgun hunting in the Republic of the Congo. In the Congo Basin, inter-group cooperation between foragers and farmers is at the centre of an exchange system maintained by traditional norms and institutions such as fictive kinship. Here, we focused on what factors predict cooperative shotgun hunting exchanges between BaYaka and Yambe. We conducted structured interviews with 48 BaYaka hunters and 18 Yambe men who organise hunts in a village along the Motaba River. We used Bayesian multilevel regression models to investigate the influence of Yambe and BaYaka attributes on probability of dyadic cooperation. We found that BaYaka men's reputations as skilled hunters and their family size each predicted cooperation in shotgun hunting, whereas there was no effect of Yambe attributes (status, wealth, family size). We discuss the results in terms of evolutionary models of men as hunters and inter-group cooperation, as well as biodiversity conservation implications.

Foraging predicts the evolution of warning coloration and mimicry in snakes

Warning coloration and Batesian mimicry are classic examples of Darwinian evolution, but empirical evolutionary patterns are often paradoxical. We test whether foraging costs predict the evolution of striking coloration by integrating genetic and ecological data for aposematic and mimetic snakes (Elapidae and Dipsadidae). Our phylogenetic comparison on a total of 432 species demonstrated that dramatic changes in coloration were well predicted by foraging strategy. Multiple tests consistently indicated that warning coloration and conspicuous mimicry were more likely to evolve in species where foraging costs of conspicuous appearance were relaxed by poor vision of their prey, concealed habitat, or nocturnal activity. Reversion to crypsis was also well predicted by ecology for elapids but not for dipsadids. In contrast to a theoretical prediction and general trends, snakes' conspicuous coloration was correlated with secretive ecology, suggesting that a selection regime underlies evolutionary patterns. We also found evidence that mimicry of inconspicuous models (pitvipers) may have evolved in association with foraging demand for crypsis. These findings demonstrate that foraging is an important factor necessary to understand the evolution, persistence, and diversity of warning coloration and mimicry of snakes, highlighting the significance of additional selective factors in solving the warning coloration paradox.

Unifying approaches from statistical genetics and phylogenetics for mapping phenotypes in structured populations

In both statistical genetics and phylogenetics, a major goal is to identify correlations between genetic loci or other aspects of the phenotype or environment and a focal trait. In these two fields, there are sophisticated but disparate statistical traditions aimed at these tasks. The disconnect between their respective approaches is becoming untenable as questions in medicine, conservation biology, and evolutionary biology increasingly rely on integrating data from within and among species, and once-clear conceptual divisions are becoming increasingly blurred. To help bridge this divide, we derive a general model describing the covariance between the genetic contributions to the quantitative phenotypes of different individuals. Taking this approach shows that standard models in both statistical genetics (e.g., Genome-Wide Association Studies; GWAS) and phylogenetic comparative biology (e.g., phylogenetic regression) can be interpreted as special cases of this more general quantitative-genetic model. The fact that these models share the same core architecture means that we can build a unified understanding of the strengths and limitations of different methods for controlling for genetic structure when testing for associations. We develop intuition for why and when spurious correlations may occur using analytical theory and conduct population-genetic and phylogenetic simulations of quantitative traits. The structural similarity of problems in statistical genetics and phylogenetics enables us to take methodological advances from one field and apply them in the other. We demonstrate this by showing how a standard GWAS technique-including both the genetic relatedness matrix (GRM) as well as its leading eigenvectors, corresponding to the principal components of the genotype matrix, in a regression model-can mitigate spurious correlations in phylogenetic analyses. As a case study of this, we re-examine an analysis testing for co-evolution of expression levels between genes across a fungal phylogeny, and show that including covariance matrix eigenvectors as covariates decreases the false positive rate while simultaneously increasing the true positive rate. More generally, this work provides a foundation for more integrative approaches for understanding the genetic architecture of phenotypes and how evolutionary processes shape it.

Variation in temperature of peak trait performance constrains adaptation of arthropod populations to climatic warming

The capacity of arthropod populations to adapt to long-term climatic warming is currently uncertain. Here we combine theory and extensive data to show that the rate of their thermal adaptation to climatic warming will be constrained in two fundamental ways. First, the rate of thermal adaptation of an arthropod population is predicted to be limited by changes in the temperatures at which the performance of four key life-history traits can peak, in a specific order of declining importance: juvenile development, adult fecundity, juvenile mortality and adult mortality. Second, directional thermal adaptation is constrained due to differences in the temperature of the peak performance of these four traits, with these differences expected to persist because of energetic allocation and life-history trade-offs. We compile a new global dataset of 61 diverse arthropod species which provides strong empirical evidence to support these predictions, demonstrating that contemporary populations have indeed evolved under these constraints. Our results provide a basis for using relatively feasible trait measurements to predict the adaptive capacity of diverse arthropod populations to geographic temperature gradients, as well as ongoing and future climatic warming.

Plumage and eggshell colouration covary with the level of sex-specific parental contributions to nest building in birds

Interspecific variation in sex-specific contributions to prenatal parental care, including avian nest building, is becoming increasingly better understood as we amass more information on more species. We examined whether sex-specific nest building contributions covary with the colouration of parents and their eggs in 521 species of Western Palearctic birds. Having colourful plumage and laying colourful eggs are costly because of the deposition of pigments in feathers and eggs and/or forming costly nanostructural substrates in feathers, and so it might be expected that those costs covary with the costs of nest building at the level of individuals and/or across species to produce of a suite of codivergent traits. Using a phylogenetically informed approach, we tested the hypothesis that species in which females alone invest energy building nests exhibit less sexual plumage dichromatism. However, we found comparative support for the opposite of this prediction. We then tested that species in which females alone build nests lay more colourful, and costlier, eggs because the dual costs of building nests and laying colourful eggs can only be borne by higher quality individuals. As expected, we found that species in which females build nests alone or together with males are more likely to lay colourfully pigmented eggs relative to species in which only males build nests. Finally, stochastic character mapping provided evidence of the repeated evolution of female-only nest building. Interspecific sex differences in plumage colouration therefore covary in a complex manner with female pre- (nest building) and post-copulatory (egg production) investment in reproduction.

Telomere length is not a useful tool for chronological age estimation in animals

Telomeres are short repetitive DNA sequences capping the ends of chromosomes. Telomere shortening occurs during cell division and may be accelerated by oxidative damage or ameliorated by telomere maintenance mechanisms. Consequently, telomere length changes with age, which was recently confirmed in a large meta-analysis across vertebrates. However, based on the correlation between telomere length and age, it was concluded that telomere length can be used as a tool for chronological age estimation in animals. Correlation should not be confused with predictability, and the current data and studies suggest that telomeres cannot be used to reliably predict individual chronological age. There are biological reasons for why there is large individual variation in telomere dynamics, which is mainly due to high susceptibility to a wide range of environmental, but also genetic factors, rendering telomeres unfeasible as a tool for age estimation. The use of telomeres for chronological age estimation is largely a misguided effort, but its occasional reappearance in the literature raises concerns that it will mislead resources in wildlife conservation.

Taste shaped the use of botanical drugs

The perception of taste and flavour (a combination of taste, smell, and chemesthesis), here also referred to as chemosensation, enables animals to find high-value foods and avoid toxins. Humans have learned to use unpalatable and toxic substances as medicines, yet the importance of chemosensation in this process is poorly understood. Here, we generate tasting-panel data for botanical drugs and apply phylogenetic generalised linear mixed models to test whether intensity and complexity of chemosensory qualities as well as particular tastes and flavours can predict ancient Graeco-Roman drug use. We found chemosensation to be strongly predictive of therapeutic use: botanical drugs with high therapeutic versatility have simple yet intense tastes and flavours, and 21 of 22 chemosensory qualities predicted at least one therapeutic use. In addition to the common notion of bitter tasting medicines, we also found starchy, musky, sweet, and soapy drugs associated with versatility. In ancient Greece and Rome, illness was thought to arise from imbalance in bodily fluids or humours, yet our study suggests that uses of drugs were based on observed physiological effects that are often consistent with modern understanding of chemesthesis and taste receptor pharmacology.

Primate social organization evolved from a flexible pair-living ancestor

Explaining the evolution of primate social organization has been fundamental to understand human sociality and social evolution more broadly. It has often been suggested that the ancestor of all primates was solitary and that other forms of social organization evolved later, with transitions being driven by various life history traits and ecological factors. However, recent research showed that many understudied primate species previously assumed to be solitary actually live in pairs, and intraspecific variation in social organization is common. We built a detailed database from primary field studies quantifying the number of social units expressing different social organizations in each population. We used Bayesian phylogenetic models to infer the probability of each social organization, conditional on several socioecological and life history predictors. Here, we show that when intraspecific variation is accounted for, the ancestral social organization of primates was inferred to be variable, with the most common social organization being pair-living but with approximately 10 to 20% of social units of the ancestral population deviating from this pattern by being solitary living. Body size and activity patterns had large effects on transitions between types of social organizations. As in other mammalian clades, pair-living is closely linked to small body size and likely more common in ancestral species. Our results challenge the assumption that ancestral primates were solitary and that pair-living evolved afterward emphasizing the importance of focusing on field data and accounting for intraspecific variation, providing a flexible statistical framework for doing so.

Insights to HIV-1 coreceptor usage by estimating HLA adaptation with Bayesian generalized linear mixed models

The mechanisms triggering the human immunodeficiency virus type I (HIV-1) to switch the coreceptor usage from CCR5 to CXCR4 during the course of infection are not entirely understood. While low CD4+ T cell counts are associated with CXCR4 usage, a predominance of CXCR4 usage with still high CD4+ T cell counts remains puzzling. Here, we explore the hypothesis that viral adaptation to the human leukocyte antigen (HLA) complex, especially to the HLA class II alleles, contributes to the coreceptor switch. To this end, we sequence the viral gag and env protein with corresponding HLA class I and II alleles of a new cohort of 312 treatment-naive, subtype C, chronically-infected HIV-1 patients from South Africa. To estimate HLA adaptation, we develop a novel computational approach using Bayesian generalized linear mixed models (GLMMs). Our model allows to consider the entire HLA repertoire without restricting the model to pre-learned HLA-polymorphisms. In addition, we correct for phylogenetic relatedness of the viruses within the model itself to account for founder effects. Using our model, we observe that CXCR4-using variants are more adapted than CCR5-using variants (p-value = 1.34e-2). Additionally, adapted CCR5-using variants have a significantly lower predicted false positive rate (FPR) by the geno2pheno[coreceptor] tool compared to the non-adapted CCR5-using variants (p-value = 2.21e-2), where a low FPR is associated with CXCR4 usage. Consequently, estimating HLA adaptation can be an asset in predicting not only coreceptor usage, but also an approaching coreceptor switch in CCR5-using variants. We propose the usage of Bayesian GLMMs for modeling virus-host adaptation in general.

Leukocyte Concentrations Are Isometric in Reptiles Unlike in Endotherms

AbstractHow do large and small reptiles defend against infections, given the consequences of body mass for physiology and disease transmission? Functionally equivalent mammalian and avian granulocytes increased disproportionately with body mass (i.e., scaled hypermetrically), such that large organisms had higher concentrations than expected by a prediction of proportional protection across sizes. However, as these scaling relationships were derived from endothermic animals, they do not necessarily inform the scaling of leukocyte concentration for ectothermic reptiles that have a different physiology and evolutionary history. Here, we asked whether and how lymphocyte and heterophil concentrations relate to body mass among more than 120 reptile species. We compared these relationships to those found in birds and mammals and to existing scaling frameworks (i.e., protecton, complexity, rate of metabolism, or safety factor hypotheses). Both lymphocyte and heterophil concentrations scaled almost isometrically among reptiles. In contrast, functionally equivalent granulocytes scaled hypermetrically and lymphocytes scaled isometrically in birds and mammals. Life history traits were also poor predictors of variation in reptilian heterophil and lymphocyte concentrations. Our results provide insight into differences in immune protection in birds and mammals relative to that in reptiles through a comparative lens. The shape of scaling relationships differs, which should be considered when modeling disease dynamics among these groups.

Genomic signatures of convergent shifts to plunge-diving behavior in birds

Understanding the genetic basis of convergence at broad phylogenetic scales remains a key challenge in biology. Kingfishers (Aves: Alcedinidae) are a cosmopolitan avian radiation with diverse colors, diets, and feeding behaviors-including the archetypal plunge-dive into water. Given the sensory and locomotor challenges associated with air-water transitions, kingfishers offer a powerful opportunity to explore the effects of convergent behaviors on the evolution of genomes and phenotypes, as well as direct comparisons between continental and island lineages. Here, we use whole-genome sequencing of 30 diverse kingfisher species to identify the genomic signatures associated with convergent feeding behaviors. We show that species with smaller ranges (i.e., on islands) have experienced stronger demographic fluctuations than those on continents, and that these differences have influenced the dynamics of molecular evolution. Comparative genomic analyses reveal positive selection and genomic convergence in brain and dietary genes in plunge-divers. These findings enhance our understanding of the connections between genotype and phenotype in a diverse avian radiation.

Phylogenetic meta-analysis reveals system-specific behavioural type-behavioural predictability correlations

The biological significance of behavioural predictability (environment-independent within-individual behavioural variation) became accepted recently as an important part of an individual's behavioural strategy besides behavioural type (individual mean behaviour). However, we do not know how behavioural type and predictability evolve. Here, we tested different evolutionary scenarios: (i) the two traits evolve independently (lack of correlations) and (ii) the two traits' evolution is constrained (abundant correlations) due to either (ii/a) proximate constraints (direction of correlations is similar) or (ii/b) local adaptations (direction of correlations is variable). We applied a set of phylogenetic meta-analyses based on 93 effect sizes across 44 vertebrate and invertebrate species, focusing on activity and risk-taking. The general correlation between behavioural type and predictability did not differ from zero. Effect sizes for correlations showed considerable heterogeneity, with both negative and positive correlations occurring. The overall absolute (unsigned) effect size was high (Zr = 0.58), and significantly exceeded the null expectation based on randomized data. Our results support the adaptive scenario: correlations between behavioural type and predictability are abundant in nature, but their direction is variable. We suggest that the evolution of these behavioural components might be constrained in a system-specific way.

The functional mechanism behind the latitudinal pattern of liana diversity: Freeze-thaw embolism reduces the ecological performance of liana species

There is a strong decrease in liana diversity along latitudinal and altitudinal gradients at a global scale, and there is a marked difference in liana diversity between tropical and temperate ecosystems. From these observations, it has been proposed that cold temperatures would restrict the ecological patterns of liana because of their vascular system's vulnerability to freeze-thaw embolism. Our objective was to establish the functional mechanism that drives the loss of liana diversity along a latitudinal temperature gradient. We evaluate the ecological performance of liana in 10 different species based on the apical growth rate, as well as functional traits associated with efficiency (maximum hydraulic conductivity and percentage conductivity lost) and safety of water transport (vessel diameter, vessel density, wood density, and root pressure). We found that at the colder (more southern) site within the latitudinal gradient, liana species showed lower performance, with a fivefold decrease in their apical growth rate as compared to the warmer (more northern) sites. We postulate that this lower performance results from a much lower water transport efficiency (26.1-fold decrease as compared to liana species that inhabit warmer sites) that results from higher freeze-thaw (37.5% of PLC) and reduction of vessel diameter (3 times narrower). These results are unmistakable evidence that cold temperature restricts liana performance: in a cold environment, liana species exhibit a strong decrease in performance, low efficiency, and higher safety of water transport. Conversely, at warmer sites, we found that liana species exhibit functional strategies associated with higher performance, higher efficiency, and lower safety of water transport capacity. This trade-off between efficiency and safety of water transport and their effects on performance could explain the latitudinal pattern of liana diversity.

A mixed-model approach for estimating drivers of microbiota community composition and differential taxonomic abundance

Next-generation sequencing (NGS) and metabarcoding approaches are increasingly applied to wild animal populations, but there is a disconnect between the widely applied generalized linear mixed model (GLMM) approaches commonly used to study phenotypic variation and the statistical toolkit from community ecology typically applied to metabarcoding data. Here, we describe the suitability of a novel GLMM-based approach for analyzing the taxon-specific sequence read counts derived from standard metabarcoding data. This approach allows decomposition of the contribution of different drivers to variation in community composition (e.g., age, season, individual) via interaction terms in the model random-effects structure. We provide guidance to implementing this approach and show how these models can identify how responsible specific taxonomic groups are for the effects attributed to different drivers. We applied this approach to two cross-sectional data sets from the Soay sheep population of St. Kilda. GLMMs showed agreement with dissimilarity-based approaches highlighting the substantial contribution of age and minimal contribution of season to microbiota community compositions, and simultaneously estimated the contribution of other technical and biological factors. We further used model predictions to show that age effects were principally due to increases in taxa of the phylum Bacteroidetes and declines in taxa of the phylum Firmicutes. This approach offers a powerful means for understanding the influence of drivers of community structure derived from metabarcoding data. We discuss how our approach could be readily adapted to allow researchers to estimate contributions of additional factors such as host or microbe phylogeny to answer emerging questions surrounding the ecological and evolutionary roles of within-host communities. IMPORTANCE NGS and fecal metabarcoding methods have provided powerful opportunities to study the wild gut microbiome. A wealth of data is, therefore, amassing across wild systems, generating the need for analytical approaches that can appropriately investigate simultaneous factors at the host and environmental scale that determine the composition of these communities. Here, we describe a generalized linear mixed-effects model (GLMM) approach to analyze read count data from metabarcoding of the gut microbiota, allowing us to quantify the contributions of multiple host and environmental factors to within-host community structure. Our approach provides outputs that are familiar to a majority of field ecologists and can be run using any standard mixed-effects modeling packages. We illustrate this approach using two metabarcoding data sets from the Soay sheep population of St. Kilda investigating age and season effects as worked examples.

Diversity, distribution and intrinsic extinction vulnerability of exploited marine bivalves

Marine bivalves are important components of ecosystems and exploited by humans for food across the world, but the intrinsic vulnerability of exploited bivalve species to global changes is poorly known. Here, we expand the list of shallow-marine bivalves known to be exploited worldwide, with 720 exploited bivalve species added beyond the 81 in the United Nations FAO Production Database, and investigate their diversity, distribution and extinction vulnerability using a metric based on ecological traits and evolutionary history. The added species shift the richness hotspot of exploited species from the northeast Atlantic to the west Pacific, with 55% of bivalve families being exploited, concentrated mostly in two major clades but all major body plans. We find that exploited species tend to be larger in size, occur in shallower waters, and have larger geographic and thermal ranges-the last two traits are known to confer extinction-resistance in marine bivalves. However, exploited bivalve species in certain regions such as the tropical east Atlantic and the temperate northeast and southeast Pacific, are among those with high intrinsic vulnerability and are a large fraction of regional faunal diversity. Our results pinpoint regional faunas and specific taxa of likely concern for management and conservation.

Human and animal dominance hierarchies show a pyramidal structure guiding adult and infant social inferences

This study investigates the structure of social hierarchies. We hypothesized that if social dominance relations serve to regulate conflicts over resources, then hierarchies should converge towards pyramidal shapes. Structural analyses and simulations confirmed this hypothesis, revealing a triadic-pyramidal motif across human and non-human hierarchies (114 species). Phylogenetic analyses showed that this pyramidal motif is widespread, with little influence of group size or phylogeny. Furthermore, nine experiments conducted in France found that human adults (N = 120) and infants (N = 120) draw inferences about dominance relations that are consistent with hierarchies' pyramidal motif. By contrast, human participants do not draw equivalent inferences based on a tree-shaped pattern with a similar complexity to pyramids. In short, social hierarchies exhibit a pyramidal motif across a wide range of species and environments. From infancy, humans exploit this regularity to draw systematic inferences about unobserved dominance relations, using processes akin to formal reasoning.

Multi-trait and multi-environment genomic prediction for flowering traits in maize: a deep learning approach

Maize (Zea mays L.), the third most widely cultivated cereal crop in the world, plays a critical role in global food security. To improve the efficiency of selecting superior genotypes in breeding programs, researchers have aimed to identify key genomic regions that impact agronomic traits. In this study, the performance of multi-trait, multi-environment deep learning models was compared to that of Bayesian models (Markov Chain Monte Carlo generalized linear mixed models (MCMCglmm), Bayesian Genomic Genotype-Environment Interaction (BGGE), and Bayesian Multi-Trait and Multi-Environment (BMTME)) in terms of the prediction accuracy of flowering-related traits (Anthesis-Silking Interval: ASI, Female Flowering: FF, and Male Flowering: MF). A tropical maize panel of 258 inbred lines from Brazil was evaluated in three sites (Cambira-2018, Sabaudia-2018, and Iguatemi-2020 and 2021) using approximately 290,000 single nucleotide polymorphisms (SNPs). The results demonstrated a 14.4% increase in prediction accuracy when employing multi-trait models compared to the use of a single trait in a single environment approach. The accuracy of predictions also improved by 6.4% when using a single trait in a multi-environment scheme compared to using multi-trait analysis. Additionally, deep learning models consistently outperformed Bayesian models in both single and multiple trait and environment approaches. A complementary genome-wide association study identified associations with 26 candidate genes related to flowering time traits, and 31 marker-trait associations were identified, accounting for 37%, 37%, and 22% of the phenotypic variation of ASI, FF and MF, respectively. In conclusion, our findings suggest that deep learning models have the potential to significantly improve the accuracy of predictions, regardless of the approach used and provide support for the efficacy of this method in genomic selection for flowering-related traits in tropical maize.

Experimental estimates of germline mutation rate in eukaryotes: a phylogenetic meta-analysis

Mutation is the ultimate source of all genetic variation, and over the last 10 years the ready availability of whole-genome sequencing has permitted direct estimation of mutation rate for many non-model species across the tree of life. In this meta-analysis, we make a comprehensive search of the literature for mutation rate estimates in eukaryotes, identifying 140 mutation accumulation (MA) and parent-offspring (PO) sequencing studies covering 134 species. Based on these data, we revisit differences in the single-nucleotide mutation (SNM) rate between different phylogenetic lineages and update the known relationships between mutation rate and generation time, genome size, and nucleotide diversity-while accounting for phylogenetic nonindependence. We do not find a significant difference between MA and PO in estimated mutation rates, but we confirm that mammal and plant lineages have higher mutation rates than arthropods and that unicellular eukaryotes have the lowest mutation rates. We find that mutation rates are higher in species with longer generation times and larger genome sizes, even when accounting for phylogenetic relationships. Moreover, although nucleotide diversity is positively correlated with mutation rate, the gradient of the relationship is significantly less than one (on a logarithmic scale), consistent with higher mutation rates in populations with smaller effective size. For the 29 species for which data are available, we find that indel mutation rates are positively correlated with nucleotide mutation rates and that short deletions are generally more common than short insertions. Nevertheless, despite recent progress, no estimates of either SNM or indel mutation rates are available for the majority of deeply branching eukaryotic lineages-or even for most animal phyla. Even among charismatic megafauna, experimental mutation rate estimates remain unknown for amphibia and scarce for reptiles and fish.

Symbioses shape feeding niches and diversification across insects

For over 300 million years, insects have relied on symbiotic microbes for nutrition and defence. However, it is unclear whether specific ecological conditions have repeatedly favoured the evolution of symbioses, and how this has influenced insect diversification. Here, using data on 1,850 microbe-insect symbioses across 402 insect families, we found that symbionts have allowed insects to specialize on a range of nutrient-imbalanced diets, including phloem, blood and wood. Across diets, the only limiting nutrient consistently associated with the evolution of obligate symbiosis was B vitamins. The shift to new diets, facilitated by symbionts, had mixed consequences for insect diversification. In some cases, such as herbivory, it resulted in spectacular species proliferation. In other niches, such as strict blood feeding, diversification has been severely constrained. Symbioses therefore appear to solve widespread nutrient deficiencies for insects, but the consequences for insect diversification depend on the feeding niche that is invaded.

Masting is uncommon in trees that depend on mutualist dispersers in the context of global climate and fertility gradients

The benefits of masting (volatile, quasi-synchronous seed production at lagged intervals) include satiation of seed predators, but these benefits come with a cost to mutualist pollen and seed dispersers. If the evolution of masting represents a balance between these benefits and costs, we expect mast avoidance in species that are heavily reliant on mutualist dispersers. These effects play out in the context of variable climate and site fertility among species that vary widely in nutrient demand. Meta-analyses of published data have focused on variation at the population scale, thus omitting periodicity within trees and synchronicity between trees. From raw data on 12 million tree-years worldwide, we quantified three components of masting that have not previously been analysed together: (i) volatility, defined as the frequency-weighted year-to-year variation; (ii) periodicity, representing the lag between high-seed years; and (iii) synchronicity, indicating the tree-to-tree correlation. Results show that mast avoidance (low volatility and low synchronicity) by species dependent on mutualist dispersers explains more variation than any other effect. Nutrient-demanding species have low volatility, and species that are most common on nutrient-rich and warm/wet sites exhibit short periods. The prevalence of masting in cold/dry sites coincides with climatic conditions where dependence on vertebrate dispersers is less common than in the wet tropics. Mutualist dispersers neutralize the benefits of masting for predator satiation, further balancing the effects of climate, site fertility and nutrient demands.

Systematic review of avian hatching failure and implications for conservation

Avian hatching failure is a widespread phenomenon, affecting around 10% of all eggs that are laid and not lost to predation, damage, or desertion. Our understanding of hatching failure is limited in terms of both its underpinning mechanisms and its occurrence across different populations. It is widely acknowledged that rates of hatching failure are higher in threatened species and in populations maintained in captivity compared to wild, non-threatened species, but these differences have rarely been quantified and any broader patterns remain unexplored. To examine the associations between threat status, management interventions, and hatching failure across populations we conducted a phylogenetically controlled multilevel meta-analysis across 231 studies and 241 species of birds. Our data set included both threatened (Critically Endangered, Endangered, and Vulnerable) and non-threatened (Near Threatened and Least Concern) species across wild and captive populations, as well as 'wild managed' ('free-living') populations. We found the mean overall rate of hatching failure across all populations to be 16.79%, with the hatching failure rate of wild, non-threatened species being 12.40%. We found that populations of threatened species experienced significantly higher mean hatching failure than populations of non-threatened species. Different levels of management were also associated with different rates of hatching failure, with wild populations experiencing the lowest rate of hatching failure, followed by wild managed populations, and populations in captivity experiencing the highest rate. Similarly, populations that were subject to the specific management interventions of artificial incubation, supplementary feeding, and artificial nest provision displayed significantly higher rates of hatching failure than populations without these interventions. The driver of this correlation between hatching failure and management remains unclear, but could be an indirect result of threatened species being more likely to have lower hatching success and also being more likely to be subject to management, indicating that conservation efforts are fittingly being focused towards the species potentially most at risk from extinction. This is the most comprehensive comparative analysis of avian hatching failure that has been conducted to date, and the first to quantify explicitly how threat status and management are associated with the rate of hatching failure in a population. We discuss the implications of our results, focusing on their potential applications to conservation. Although we identified several factors clearly associated with variation in hatching failure, a significant amount of heterogeneity was not explained by our meta-analytical model, indicating that other factors influencing hatching failure were not included here. We discuss what these factors might be and suggest avenues for further research. Finally, we discuss the inconsistency in how hatching failure is defined and reported within the literature, and propose a standardised definition to be used in future studies which will enable better comparison across populations and ensure that the most accurate information is used to support management decisions.

Ornaments indicate parasite load only if they are dynamic or parasites are contagious

Choosing to mate with an infected partner has several potential fitness costs, including disease transmission and infection-induced reductions in fecundity and parental care. By instead choosing a mate with no, or few, parasites, animals avoid these costs and may also obtain resistance genes for offspring. Within a population, then, the quality of sexually selected ornaments on which mate choice is based should correlate negatively with the number of parasites with which a host is infected ("parasite load"). However, the hundreds of tests of this prediction yield positive, negative, or no correlation between parasite load and ornament quality. Here, we use phylogenetically controlled meta-analysis of 424 correlations from 142 studies on a wide range of host and parasite taxa to evaluate explanations for this ambiguity. We found that ornament quality is weakly negatively correlated with parasite load overall, but the relationship is more strongly negative among ornaments that can dynamically change in quality, such as behavioral displays and skin pigmentation, and thus can accurately reflect current parasite load. The relationship was also more strongly negative among parasites that can transmit during sex. Thus, the direct benefit of avoiding parasite transmission may be a key driver of parasite-mediated sexual selection. No other moderators, including methodological details and whether males exhibit parental care, explained the substantial heterogeneity in our data set. We hope to stimulate research that more inclusively considers the many and varied ways in which parasites, sexual selection, and epidemiology intersect.

Phylogenomic comparative methods: Accurate evolutionary inferences in the presence of gene tree discordance

Phylogenetic comparative methods have long been a mainstay of evolutionary biology, allowing for the study of trait evolution across species while accounting for their common ancestry. These analyses typically assume a single, bifurcating phylogenetic tree describing the shared history among species. However, modern phylogenomic analyses have shown that genomes are often composed of mosaic histories that can disagree both with the species tree and with each other-so-called discordant gene trees. These gene trees describe shared histories that are not captured by the species tree, and therefore that are unaccounted for in classic comparative approaches. The application of standard comparative methods to species histories containing discordance leads to incorrect inferences about the timing, direction, and rate of evolution. Here, we develop two approaches for incorporating gene tree histories into comparative methods: one that constructs an updated phylogenetic variance-covariance matrix from gene trees, and another that applies Felsenstein's pruning algorithm over a set of gene trees to calculate trait histories and likelihoods. Using simulation, we demonstrate that our approaches generate much more accurate estimates of tree-wide rates of trait evolution than standard methods. We apply our methods to two clades of the wild tomato genus Solanum with varying rates of discordance, demonstrating the contribution of gene tree discordance to variation in a set of floral traits. Our approaches have the potential to be applied to a broad range of classic inference problems in phylogenetics, including ancestral state reconstruction and the inference of lineage-specific rate shifts.

Eggshell composition and surface properties of avian brood-parasitic species compared with non-parasitic species

The eggs of avian obligate brood-parasitic species have multiple adaptations to deceive hosts and optimize development in host nests. While the structure and composition of the eggshell in all birds is essential for embryo growth and protection from external threats, parasitic eggs may face specific challenges such as high microbial loads, rapid laying and ejection by the host parents. We set out to assess whether eggshells of avian brood-parasitic species have either (i) specialized structural properties, to meet the demands of a brood-parasitic strategy or (ii) similar structural properties to eggs of their hosts, due to the similar nest environment. We measured the surface topography (roughness), wettability (how well surfaces repel water) and calcium content of eggshells of a phylogenetically and geographically diverse range of brood-parasitic species (representing four of the seven independent lineages of avian brood-parasitic species), their hosts and close relatives of the parasites. These components of the eggshell structure have been demonstrated previously to influence such factors as the risk of microbial infection and overall shell strength. Within a phylogenetically controlled framework, we found no overall significant differences in eggshell roughness, wettability and calcium content between (i) parasitic and non-parasitic species, or (ii) parasitic species and their hosts. Both the wettability and calcium content of the eggs from brood-parasitic species were not more similar to those of their hosts' eggs than expected by chance. By contrast, the mean surface roughness of the eggs of brood-parasitic species was more similar to that of their hosts' eggs than expected by chance, suggesting brood-parasitic species may have evolved to lay eggs that match the host nest environment for this trait. The lack of significant overall differences between parasitic and non-parasitic species, including hosts, in the traits we measured, suggests that phylogenetic signal, as well as general adaptations to the nest environment and for embryo development, outweigh any influence of a parasitic lifestyle on these eggshell properties.

Measuring the shape of mortality across animals and plants: Alternatives to H entropy metrics reveal hidden type IV survivorship curves and associations with parental care at macro-ecological scales

The shape of mortality, or how mortality is spread across an organism's life course, is fundamental to a range of biological processes, with attempts to quantify it rooted in ecology, evolution, and demography. One approach to quantify the distribution of mortality over an organism's life is the use of entropy metrics whose values are interpreted within the classical framework of survivorship curves ranging from type I distributions, with mortality concentrated in late life stages, to type III survivorship curves associated with high early stage mortality. However, entropy metrics were originally developed using restricted taxonomic groups and the behavior of entropy metrics over larger scales of variation may make them unsuitable for wider-ranging contemporary comparative studies. Here, we revisit the classic survivorship framework and, using a combination of simulations and comparative analysis of demography data spanning the animal and plant kingdoms, we show that commonly used entropy metrics cannot distinguish between the most extreme survivorship curves, which in turn can mask important macroecological patterns. We show how using H entropy masks a macroecological pattern of how parental care is associated with type I and type II species and for macroecological studies recommend the use of metrics, such as measures of area under the curve. Using frameworks and metrics that capture the full range of variation of survivorship curves will aid in our understanding of the links between the shape of mortality, population dynamics, and life history traits.

Insights into mammalian TE diversity through the curation of 248 genome assemblies

We examined transposable element (TE) content of 248 placental mammal genome assemblies, the largest de novo TE curation effort in eukaryotes to date. We found that although mammals resemble one another in total TE content and diversity, they show substantial differences with regard to recent TE accumulation. This includes multiple recent expansion and quiescence events across the mammalian tree. Young TEs, particularly long interspersed elements, drive increases in genome size, whereas DNA transposons are associated with smaller genomes. Mammals tend to accumulate only a few types of TEs at any given time, with one TE type dominating. We also found association between dietary habit and the presence of DNA transposon invasions. These detailed annotations will serve as a benchmark for future comparative TE analyses among placental mammals.

Variation in ontogenetic trajectories of limb dimensions in humans is attributable to both climatic effects and neutral evolution

Previous studies showed that there is variation in ontogenetic trajectories of human limb dimensions and proportions. However, little is known about the evolutionary significance of this variation. This study used a global sample of modern human immature long bone measurements and a multivariate linear mixed-effects model to study 1) whether the variation in ontogenetic trajectories of limb dimensions is consistent with ecogeographic predictions and 2) the effects of different evolutionary forces on the variation in ontogenetic trajectories. We found that genetic relatedness arising from neutral (nonselective) evolution, allometric variation associated with the change in size, and directional effects from climate all contributed to the variation in ontogenetic trajectories of all major long bone dimensions in modern humans. After accounting for the effects of neutral evolution and holding other effects considered in the current study constant, extreme temperatures have weak, positive associations with diaphyseal length and breadth measurements, while mean temperature shows negative associations with diaphyseal dimensions. The association with extreme temperatures fits the expectations of ecogeographic rules, while the association with mean temperature may explain the observed among-group variation in intralimb indices. The association with climate is present throughout ontogeny, suggesting an explanation of adaptation by natural selection as the most likely cause. On the other hand, genetic relatedness among groups, as structured by neutral evolutionary factors, is an important consideration when interpreting skeletal morphology, even for nonadult individuals.

Complex plumages spur rapid color diversification in kingfishers (Aves: Alcedinidae)

Colorful signals in nature provide some of the most stunning examples of rapid phenotypic evolution. Yet, studying color pattern evolution has been historically difficult owing to differences in perceptual ability of humans and analytical challenges with studying how complex color patterns evolve. Island systems provide a natural laboratory for testing hypotheses about the direction and magnitude of phenotypic change. A recent study found that plumage colors of island species are darker and less complex than continental species. Whether such shifts in plumage complexity are associated with increased rates of color evolution remains unknown. Here, we use geometric morphometric techniques to test the hypothesis that plumage complexity and insularity interact to influence color diversity in a species-rich clade of colorful birds-kingfishers (Aves: Alcedinidae). In particular, we test three predictions: (1) plumage complexity enhances interspecific rates of color evolution, (2) plumage complexity is lower on islands, and (3) rates of plumage color evolution are higher on islands. Our results show that more complex plumages result in more diverse colors among species and that island species have higher rates of color evolution. Importantly, we found that island species did not have more complex plumages than their continental relatives. Thus, complexity may be a key innovation that facilitates evolutionary response of individual color patches to distinct selection pressures on islands, rather than being a direct target of selection itself. This study demonstrates how a truly multivariate treatment of color data can reveal evolutionary patterns that might otherwise go unnoticed.

Innovation across 13 ungulate species: problem solvers are less integrated in the social group and less neophobic

Innovation is the ability to solve new problems or find novel solutions to familiar problems, and it is known to provide animals with crucial fitness benefits. Although this ability has been extensively studied in some taxa, the factors that predict innovation within and across species are still largely unclear. In this study, we used a novel foraging task to test 111 individuals belonging to 13 ungulate species-a still understudied taxon. To solve the task, individuals had to open transparent and opaque cups with food rewards, by removing their cover. We assessed whether individual factors (neophobia, social integration, sex, age, rank) and socio-ecological factors (dietary breadth, fission-fusion dynamics, domestication, group size) predicted participation and performance in the task. Using a phylogenetic approach, we showed that success was higher for less neophobic and socially less integrated individuals. Moreover, less neophobic individuals, individuals of domesticated species and having higher fission-fusion dynamics were more likely to participate in the task. These results are in line with recent literature suggesting a central role of sociality and personality traits to successfully deal with novel challenges, and confirm ungulates as a promising taxon to test evolutionary theories with a comparative approach.

Advanced methods and implementations for the meta-analyses of animal models: Current practices and future recommendations

Meta-analytic techniques have been widely used to synthesize data from animal models of human diseases and conditions, but these analyses often face two statistical challenges due to complex nature of animal data (e.g., multiple effect sizes and multiple species): statistical dependency and confounding heterogeneity. These challenges can lead to unreliable and less informative evidence, which hinders the translation of findings from animal to human studies. We present a literature survey of meta-analysis using animal models (animal meta-analysis), showing that these issues are not adequately addressed in current practice. To address these challenges, we propose a meta-analytic framework based on multilevel (linear mixed-effects) models. Through conceptualization, formulations, and worked examples, we illustrate how this framework can appropriately address these issues while allowing for testing new questions. Additionally, we introduce other advanced techniques such as multivariate models, robust variance estimation, and meta-analysis of emergent effect sizes, which can deliver robust inferences and novel biological insights. We also provide a tutorial with annotated R code to demonstrate the implementation of these techniques.

Evolutionary history of hominin brain size and phylogenetic comparative methods

An absolutely and relatively large brain has traditionally been viewed as a distinctive characteristic of the Homo genus, with anatomically modern humans presented at the apex of a long line of progressive increases in encephalization. Many studies continue to focus attention on increasing brain size in the Homo genus, while excluding measures of absolute and relative brain size of more geologically recent, smaller brained, hominins such as Homo floresiensis, and Homo naledi and smaller brained Homo erectus specimens. This review discusses the benefits of using phylogenetic comparative methods to trace the diverse changes in hominin brain evolution and the drawbacks of not doing so.

Extended parental provisioning and variation in vertebrate brain sizes

Large brains provide adaptive cognitive benefits but require unusually high, near-constant energy inputs and become fully functional well after their growth is completed. Consequently, young of most larger-brained endotherms should not be able to independently support the growth and development of their own brains. This paradox is solved if the evolution of extended parental provisioning facilitated brain size evolution. Comparative studies indeed show that extended parental provisioning coevolved with brain size and that it may improve immature survival. The major role of extended parental provisioning supports the idea that the ability to sustain the costs of brains limited brain size evolution.

Correlated evolution of social organization and lifespan in mammals

Discerning the relationship between sociality and longevity would permit a deeper understanding of how animal life history evolved. Here, we perform a phylogenetic comparative analysis of ~1000 mammalian species on three states of social organization (solitary, pair-living, and group-living) and longevity. We show that group-living species generally live longer than solitary species, and that the transition rate from a short-lived state to a long-lived state is higher in group-living than non-group-living species, altogether supporting the correlated evolution of social organization and longevity. The comparative brain transcriptomes of 94 mammalian species identify 31 genes, hormones and immunity-related pathways broadly involved in the association between social organization and longevity. Further selection features reveal twenty overlapping pathways under selection for both social organization and longevity. These results underscore a molecular basis for the influence of the social organization on longevity.

Human activities favour prolific life histories in both traded and introduced vertebrates

Species' life histories determine population demographics and thus the probability that introduced populations establish and spread. Life histories also influence which species are most likely to be introduced, but how such 'introduction biases' arise remains unclear. Here, we investigate how life histories affect the probability of trade and introduction in phylogenetic comparative analyses across three vertebrate classes: mammals, reptiles and amphibians. We find that traded species have relatively high reproductive rates and long reproductive lifespans. Within traded species, introduced species have a more extreme version of this same life history profile. Species in the pet trade also have long reproductive lifespans but lack 'fast' traits, likely reflecting demand for rare species which tend to have slow life histories. We identify multiple species not yet traded or introduced but with life histories indicative of high risk of future trade, introduction and potentially invasion. Our findings suggest that species with high invasion potential are favoured in the wildlife trade and therefore that trade regulation is crucial for preventing future invasions.

Sex differences in cooperative coalitions: a mammalian perspective

In group-living species, cooperative tactics can offset asymmetries in resource-holding potential between individuals and alter the outcome of intragroup conflicts. Differences in the kinds of competitive pressures that males and females face might influence the benefits they gain from forming intragroup coalitions. We predicted that there would be a female bias in intragroup coalitions because females (1) are more like to live with kin than males are, and (2) compete over resources that are more readily shared than resources males compete over. We tested this main prediction using information about coalition formation across mammalian species and phylogenetic comparative analyses. We found that for nearly all species in which intragroup coalitions occur, members of both sexes participate, making this the typical mammalian pattern. The presence and frequency of female or male coalitions were not strongly associated with key socio-ecological factors like resource defensibility, sexual dimorphism or philopatry. This suggests that once the ability to form intragroup coalitions emerges in one sex, it is likely to emerge in the other sex as well and that there is no strong phylogenetic legacy of sex differences in this form of cooperation. This article is part of the theme issue 'Cooperation among women: evolutionary and cross-cultural perspectives'.

A Phylogenetic Framework to Simulate Synthetic Interspecies RNA-Seq Data

Interspecies RNA-Seq datasets are increasingly common, and have the potential to answer new questions about the evolution of gene expression. Single-species differential expression analysis is now a well-studied problem that benefits from sound statistical methods. Extensive reviews on biological or synthetic datasets have provided the community with a clear picture on the relative performances of the available methods in various settings. However, synthetic dataset simulation tools are still missing in the interspecies gene expression context. In this work, we develop and implement a new simulation framework. This tool builds on both the RNA-Seq and the phylogenetic comparative methods literatures to generate realistic count datasets, while taking into account the phylogenetic relationships between the samples. We illustrate the usefulness of this new framework through a targeted simulation study, that reproduces the features of a recently published dataset, containing gene expression data in adult eye tissue across blind and sighted freshwater crayfish species. Using our simulated datasets, we perform a fair comparison of several approaches used for differential expression analysis. This benchmark reveals some of the strengths and weaknesses of both the classical and phylogenetic approaches for interspecies differential expression analysis, and allows for a reanalysis of the crayfish dataset. The tool has been integrated in the R package compcodeR, freely available on Bioconductor.

The evolution of sexually dimorphic traits in ecological gradients: an interplay between natural and sexual selection in hummingbirds

Traits that exhibit differences between the sexes have been of special interest in the study of phenotypic evolution. Classic hypotheses explain sexually dimorphic traits via intra-sexual competition and mate selection, yet natural selection may also act differentially on the sexes to produce dimorphism. Natural selection can act either through physiological and ecological constraints on one of the sexes, or by modulating the strength of sexual/social selection. This predicts an association between the degree of dimorphism and variation in ecological environments. Here, we characterize the variation in hummingbird dimorphism across ecological gradients using rich databases of morphology, colouration and song. We show that morphological dimorphism decreases with elevation in the understorey and increases with elevation in mixed habitats, that dichromatism increases at high altitudes in open and mixed habitats, and that song is less complex in mixed habitats. Our results are consistent with flight constraints, lower predation pressure at high elevations and with habitat effects on song transmission. We also show that dichromatism and song complexity are positively associated, while tail dimorphism and song complexity are negatively associated. Our results suggest that key ecological factors shape sexually dimorphic traits, and that different communication modalities do not always evolve in tandem.

Leaf habit affects the distribution of drought sensitivity but not water transport efficiency in the tropics

Considering the global intensification of aridity in tropical biomes due to climate change, we need to understand what shapes the distribution of drought sensitivity in tropical plants. We conducted a pantropical data synthesis representing 1117 species to test whether xylem-specific hydraulic conductivity (K_S ), water potential at leaf turgor loss (Ψ_TLP ) and water potential at 50% loss of K_S (Ψ_P50 ) varied along climate gradients. The Ψ_TLP and Ψ_P50 increased with climatic moisture only for evergreen species, but K_S did not. Species with high Ψ_TLP and Ψ_P50 values were associated with both dry and wet environments. However, drought-deciduous species showed high Ψ_TLP and Ψ_P50 values regardless of water availability, whereas evergreen species only in wet environments. All three traits showed a weak phylogenetic signal and a short half-life. These results suggest strong environmental controls on trait variance, which in turn is modulated by leaf habit along climatic moisture gradients in the tropics.

The evolution of sex similarities in social signals: Climatic seasonality is associated with lower sexual dimorphism and greater elaboration of female and male signals in antbirds (Thamnophilidae)

Selection on signals that mediate social competition varies with resource availability. Climate regulates resource availability, which may affect the strength of competition and selection on signals. Traditionally, this meant that more seasonal, colder, or dryer-overall harsher-environments should favor the elaboration of male signals under stronger male-male competition, increasing sexual dimorphism. However, females also use signals to compete; thus, harsher environments could strengthen competition and favor elaboration of signals in both sexes, decreasing sexual dimorphism. Alternatively, harsher environments could decrease sexual dimorphism due to scarcer resources to invest in signal elaboration in both sexes. We evaluated these contrasting hypotheses in antbirds, a family of Neotropical passerines that varies in female and male signals and occurs across diverse climatic regimes. We tested the association of sexual dimorphism of plumage coloration and songs with temperature, precipitation, and their seasonality. We found that greater seasonality is associated with lower sexual dimorphism in plumage coloration and greater elaboration of visual signals in both sexes, but not acoustic signals. Our results suggest that greater seasonality may be associated with convergent elaboration of female and male visual signals, highlighting the role of signals of both sexes in the evolution of sexual dimorphism.

Within-individual repeatability in telomere length: A meta-analysis in nonmammalian vertebrates

Telomere length is increasingly used as a biomarker of long-term somatic state and future survival prospects. While most studies have overlooked this aspect, biological interpretations based on a given telomere length will benefit from considering the level of within-individual repeatability of telomere length through time. Therefore, we conducted a meta-analysis on 74 longitudinal studies in nonmammalian vertebrates, with the aim to establish the current pattern of within-individual repeatability in telomere length and to identify the methodological (e.g., qPCR/TRF) and biological factors (e.g., age class, phylogeny) that may affect it. While the median within-individual repeatability of telomere length was moderate to high (R = 0.55; 95% CI: 0.05-0.95; N = 82), marked heterogeneity between studies was evident. Measurement method affected the repeatability estimate strongly, with TRF-based studies exhibiting high repeatability (R = 0.80; 95% CI: 0.34-0.96; N = 25), while repeatability of qPCR-based studies was markedly lower and more variable (R = 0.46; 95% CI: 0.04-0.82; N = 57). While phylogeny explained some variance in repeatability, phylogenetic signal was not significant (λ = 0.32; 95% CI: 0.00-0.83). None of the biological factors investigated here significantly explained variation in the repeatability of telomere length, being potentially obscured by methodological differences. Our meta-analysis highlights the high variability in within-individual repeatability estimates between studies and the need to put more effort into separating technical and biological explanations. This is important to better understand to what extent biological factors can affect the repeatability of telomere length and thus the interpretation of telomere length data.