The erratic mitochondrial clock: variations of mutation rate, not population size, affect mtDNA diversity across birds and mammals

Background  

During the last ten years, major advances have been made in characterizing and understanding the evolution of mitochondrial DNA, the most popular marker of molecular biodiversity. Several important results were recently reported using mammals as model organisms, including (i) the absence of relationship between mitochondrial DNA diversity and life-history or ecological variables, (ii) the absence of prominent adaptive selection, contrary to what was found in invertebrates, and (iii) the unexpectedly large variation in neutral substitution rate among lineages, revealing a possible link with species maximal longevity. We propose to challenge these results thanks to the bird/mammal comparison. Direct estimates of population size are available in birds, and this group presents striking life-history trait differences with mammals (higher mass-specific metabolic rate and longevity). These properties make birds the ideal model to directly test for population size effects, and to discriminate between competing hypotheses about the causes of substitution rate variation.     

Threatened fish species in the Northeast Atlantic are functionally rare

Aim

The criteria used to define the International Union for Conservation of Nature (IUCN) Red List categories are essentially based on demographic parameters at the species level, but they do not integrate species' traits or their roles in ecosystems. Consequently, current IUCN-based protection measures may not be sufficient to conserve ecosystem functioning and services. Some species may have a singular combination of traits associated with unique functions. Such functionally distinct species are increasingly recognized as a key facet of biodiversity since they are, by definition, functionally irreplaceable. The aim of this study is to investigate whether threatened species are also functionally rare and to identify which traits determine extinction risk.

Location

European continental shelf seas.

Time period

1984–2020.

Major taxa studied

Marine fish.

Methods

Using newly compiled trait information of 425 marine fish species in European waters, and more than 30 years of scientific bottom trawl surveys, we estimated the functional distinctiveness, restrictedness and scarcity of each species and cross-referenced it with their IUCN conservation status.

Results

In European continental shelf seas, 38% of the species threatened with extinction (9 out of 24 species) were identified as the most functionally distinct. By mapping extinction risk in the multidimensional species trait space, we showed that species with the greatest risk of extinction are long-lived and of high trophic level. We also identified that the most functionally distinct species are sparsely distributed (4% of the total area on average) and have scarce abundances (<1% of the relative mean abundance of common species).

Main Conclusions

Because a substantial proportion of threatened species are functionally distinct and thus may play unique roles in ecosystem functioning, we stress that species traits—especially functional rarity—should become an indispensable step in the development of conservation management plans.     

Global drivers of population density in terrestrial vertebrates

Aim

Although the effects of life history traits on population density have been investigated widely, how spatial environmental variation influences population density for a large range of organisms and at a broad spatial scale is poorly known. Filling this knowledge gap is crucial for global species management and conservation planning and to understand the potential impact of environmental changes on multiple species.

Location

Global.

Time period

Present.

Major taxa studied

Terrestrial amphibians, reptiles, birds and mammals.

Methods

We collected population density estimates for a range of terrestrial vertebrates, including 364 estimates for amphibians, 850 for reptiles, 5,667 for birds and 7,651 for mammals. We contrasted the importance of life history traits and environmental predictors using mixed models and tested different hypotheses to explain the variation in population density for the four groups. We assessed the predictive accuracy of models through cross‐validation and mapped the partial response of vertebrate population density to environmental variables globally.

Results

Amphibians were more abundant in wet areas with high productivity levels, whereas reptiles showed relatively higher densities in arid areas with low productivity and stable temperatures. The density of birds and mammals was typically high in temperate wet areas with intermediate levels of productivity. The models showed good predictive abilities, with pseudo‐R² ranging between 0.68 (birds) and 0.83 (reptiles).

Main conclusions

Traits determine most of the variation in population density across species, whereas environmental conditions explain the intraspecific variation across populations. Species traits, resource availability and climatic stability have a different influence on the population density of the four groups. These models can be used to predict the average species population density over large areas and be used to explore macroecological patterns and inform conservation analyses.     

Age at the onset of senescence in birds and mammals is predicted by early-life performance

Life-history theory predicts that traits involved in maturity, reproduction and survival correlate along a fast–slow continuum of life histories. Evolutionary theories and empirical results indicate that senescence-related traits vary along this continuum, with slow species senescing later and at a slower pace than fast species. Because senescence patterns are typically difficult to estimate from studies in the wild, here we propose to predict the associated trait values in the frame of life-history theory. From a comparative analysis based on 81 free-ranging populations of 72 species of birds and mammals, we find that a nonlinear combination of fecundity, age at first reproduction and survival over the immature stage can account for ca two-thirds of the variance in the age at the onset of actuarial senescence. Our life-history model performs better than a model predicting the onset based on generation time, and it only includes life-history traits during early life as explanatory variables, i.e. parameters that are both theoretically expected to shape senescence and are measurable within relatively short studies. We discuss the good-fit of our life-history model to the available data in the light of current evolutionary theories of senescence. We further use it to evaluate whether studies that provided no evidence for senescence lasted long enough to include the onset of senescence.     

Imputing missing data in plant traits: A guide to improve gap-filling

Aim

Globally distributed plant trait data are increasingly used to understand relationships between biodiversity and ecosystem processes. However, global trait databases are sparse because they are compiled from many, mostly small databases. This sparsity in both trait space completeness and geographical distribution limits the potential for both multivariate and global analyses. Thus, ‘gap-filling’ approaches are often used to impute missing trait data. Recent methods, like Bayesian hierarchical probabilistic matrix factorization (BHPMF), can impute large and sparse data sets using side information. We investigate whether BHPMF imputation leads to biases in trait space and identify aspects influencing bias to provide guidance for its usage.

Innovation

We use a fully observed trait data set from which entries are randomly removed, along with extensive but sparse additional data. We use BHPMF for imputation and evaluate bias by: (1) accuracy (residuals, RMSE, trait means), (2) correlations (bi- and multivariate) and (3) taxonomic and functional clustering (valuewise, uni- and multivariate). BHPMF preserves general patterns of trait distributions but induces taxonomic clustering. Data set–external trait data had little effect on induced taxonomic clustering and stabilized trait–trait correlations.

Main Conclusions

Our study extends the criteria for the evaluation of gap-filling beyond RMSE, providing insight into statistical data structure and allowing better informed use of imputed trait data, with improved practice for imputation. We expect our findings to be valuable beyond applications in plant ecology, for any study using hierarchical side information for imputation.     

Early response of herbaceous vegetation to Rhododendron ponticum subsp. baeticum invasion in European Atlantic forests

Questions

Rhododendron ponticum subsp. baeticum is an invasive shrub of growing concern in continental Europe, but little is known about its impact on native plant communities. Here we ask: do environmental conditions differ between forest stands invaded by it and uninvaded stands? Do these differences correlate with R. ponticum's cover? Are these differences associated with differences in taxonomic and functional diversity of vascular plant species of the herb layer? Can these vegetation changes be explained by the sorting of certain life-history traits by R. ponticum-induced environmental changes?

Location

Several forests invaded by R. ponticum in the French Atlantic domain.

Methods

We recorded vegetation composition and a number of environmental variables in 400-m² plots that were established in 64 paired forest stands (32 invaded vs 32 uninvaded). We compiled traits from existing databases. We computed several metrics of taxonomic and functional diversity. We compared environmental variables and diversity metrics between invaded and uninvaded stands. We used correlation and regression analyses to relate them with R. ponticum's cover. We ran RLQ and fourth-corner analyses to explore the relationships between R. ponticum invasion, environmental variables, species traits, and vegetation composition.

Results

Independent of its abundance, R. ponticum invasion was associated with lower light arrival at the forest floor and increased litter thickness. Concomitantly, species richness and diversity and trait diversity were reduced. The major driver of species assemblages was soil pH, which strongly interacted with the invasion gradient. R. ponticum did not sort species according to traits associated with shade tolerance and thick-litter tolerance. However, tree and shrub saplings were more abundant in invaded than uninvaded stands, at the expense of graminoid and fern species.

Conclusions

As R. ponticum becomes the dominant shrub, it exerts new selection forces on life-history traits of extant species, mostly via reduced light availability, increased litter thickness, and physical competition, thereby reducing taxonomic and functional diversity of the herb layer, without impeding tree and shrub self-regeneration, at least in the short term.     

BioDeepTime: A database of biodiversity time series for modern and fossil assemblages

Motivation

We have little understanding of how communities respond to varying magnitudes and rates of environmental perturbations across temporal scales. BioDeepTime harmonizes assemblage time series of presence and abundance data to help facilitate investigations of community dynamics across timescales and the response of communities to natural and anthropogenic stressors. BioDeepTime includes time series of terrestrial and aquatic assemblages of varying spatial and temporal grain and extent from the present-day to millions of years ago.

Main Types of Variables Included

BioDeepTime currently contains 7,437,847 taxon records from 10,062 assemblage time series, each with a minimum of 10 time steps. Age constraints, sampling method, environment and taxonomic scope are provided for each time series.

Spatial Location and Grain

The database includes 8752 unique sampling locations from freshwater, marine and terrestrial ecosystems. Spatial grain represented by individual samples varies from quadrats on the order of several cm² to grid cells of ~100 km².

Time Period and Grain

BioDeepTime in aggregate currently spans the last 451 million years, with the 10,062 modern and fossil assemblage time series ranging in extent from years to millions of years. The median extent of modern time series is 18.7 years and for fossil series is 54,872 years. Temporal grain, the time encompassed by individual samples, ranges from days to tens of thousands of years.

Major Taxa and Level of Measurement

The database contains information on 28,777 unique taxa with 4,769,789 records at the species level and another 271,218 records known to the genus level, including time series of benthic and planktonic foraminifera, coccolithophores, diatoms, ostracods, plants (pollen), radiolarians and other invertebrates and vertebrates. There are to date 7012 modern and 3050 fossil time series in BioDeepTime.

Software Format

SQLite, Comma-separated values.     

Environmental variation is a major predictor of global trait turnover in mammals

Aim

To evaluate how environment and evolutionary history interact to influence global patterns of mammal trait diversity (a combination of 14 morphological and life‐history traits).

Location

The global terrestrial environment.

Taxon

Terrestrial mammals.

Methods

We calculated patterns of spatial turnover for mammalian traits and phylogenetic lineages using the mean nearest taxon distance. We then used a variance partitioning approach to establish the relative contribution of trait conservatism, ecological adaptation and clade specific ecological preferences on global trait turnover.

Results

We provide a global scale analysis of trait turnover across mammalian terrestrial assemblages, which demonstrates that phylogenetic turnover by itself does not predict trait turnover better than random expectations. Conversely, trait turnover is consistently more strongly associated with environmental variation than predicted by our null models. The influence of clade‐specific ecological preferences, reflected by the shared component of phylogenetic turnover and environmental variation, was considerably higher than expectations. Although global patterns of trait turnover are dependent on the trait under consideration, there is a consistent association between trait turnover and environmental predictive variables, regardless of the trait considered.

Main conclusions

Our results suggest that changes in phylogenetic composition are not always coupled with changes in trait composition on a global scale and that environmental conditions are strongly associated with patterns of trait composition across species assemblages, both within and across phylogenetic clades.     

Pre-adaptation and adaptation shape trait-environment matching in the Neotropics

Aim

Functional traits shape the distribution of taxa across environments. However, it remains unclear whether trait and environmental niche evolution are correlated, and what happened first: trait change facilitating environment shifts (pre-adaptation) or environmental change leading to trait change (adaptation). We focus on a species-rich Neotropical legume radiation to shed light on this enigma.

Location

Neotropics.

Time Period

Cenozoic.

Major Taxa Studied

Fabaceae: Papilionoidae: Swartzia.

Methods

We assembled leaflet, fruit and petal size data from monographs and herbarium collections for 86 to 96% of the c. 180 Swartzia species, inferred a dated Swartzia phylogenetic tree from existing DNA sequences covering 38% of the species and integrated these with distribution, soil and climate data. We used phylogenetic linear regression to quantify trait–environment relationships and applied comparative methods to evaluate modes of correlated evolution between traits and environments.

Results

Leaflet and petal size were strongly linked to climate, while fruit size was not associated with climate or soil characteristics. Evolutionary transitions to relatively low rainfall and low temperature environments were conditional on the evolution of small leaflets, whereas transitions to wet and warm environments were preceded by the evolution of larger leaflets. In contrast, transitions to the warmest or coldest environments were followed, rather than preceded, by petal loss.

Main Conclusion

Our results show that the macroevolution of functional traits has influenced the broad-scale distribution of Swartzia across Neotropical rainforest, seasonally dry, montane and inundated habitats. We suggest that trait evolution is conditional on environmental change but both pre-adaptive and adaptive processes may occur. These processes are important to understand the distribution of diversity at both regional (e.g. Amazonia) and global biogeographical scales.     

Hung out to dry: choice of priority ecoregions for conserving threatened neotropical anurans depends on life-history traits

Background

In the Neotropics, nearly 35% of amphibian species are threatened by habitat loss, habitat fragmentation, and habitat split; anuran species with different developmental modes respond to habitat disturbance in different ways. This entails broad-scale strategies for conserving biodiversity and advocates for the identification of high conservation-value regions that are significant in a global or continental context and that could underpin more detailed conservation assessments towards such areas.

Methodology/Principal Findings

We identified key ecoregion sets for anuran conservation using an algorithm that favors complementarity (beta-diversity) among ecoregions. Using the WWF''s Wildfinder database, which encompasses 700 threatened anuran species in 119 Neotropical ecoregions, we separated species into those with aquatic larvae (AL) or terrestrial development (TD), as this life-history trait affects their response to habitat disturbance. The conservation target of 100% of species representation was attained with a set of 66 ecoregions. Among these, 30 were classified as priority both for species with AL and TD, 26 were priority exclusively for species with AL, and 10 for species with TD only. Priority ecoregions for both developmental modes are concentrated in the Andes and in Mesoamerica. Ecoregions important for conserving species with AL are widely distributed across the Neotropics. When anuran life histories were ignored, species with AL were always underrepresented in priority sets.

Conclusions/Significance

The inclusion of anuran developmental modes in prioritization analyses resulted in more comprehensive coverage of priority ecoregions–especially those essential for species that require an aquatic habitat for their reproduction–when compared to usual analyses that do not consider this life-history trait. This is the first appraisal of the most important regions for conservation of threatened Neotropical anurans. It is also a first endeavor including anuran life-history traits in priority area-selection for conservation, with a clear gain in comprehensiveness of the selection process.     

Energetics and life-history of bats in comparison to small mammals

Mammals can be aligned along a slow-fast life-history continuum and a low–high metabolic rate continuum based on their traits. Small non-volant mammals occupy the fast/high end in both continua with high reproductive rates and short life spans linked with high mass-specific metabolic rates. Bats occupy the high end of the metabolic continuum, but the slow end of the life-history continuum with low reproductive rates and long life spans. Typically, both continua are linked, and similar life-history traits of species are reflected in more similar metabolic rates. We therefore hypothesized that metabolic rates are similar in species with similar life-history traits. Resting metabolic rates (RMR) were measured for three ecologically and morphologically similar sympatric bat species (Myotis nattereri, M. bechsteinii, and Plecotus auritus; Vespertilionidae) and compared to data from other similar-sized, temperate insectivorous mammals with other life-history strategies. The bat species share similar life-histories and RMRs, both of which differ from the remaining mammals and therefore supporting our hypothesis. To verify that bats are similar in RMR, two energetically contrasting periods were compared. RMRs in post-lactating females did not differ between bat species. It was, however, positively correlated with parasite load in both Myotis species. However, RMRs differed during energy-demanding pregnancy where M. nattereri had the significantly lowest RMR, suggesting metabolic compensation as an energy-saving strategy. We conclude that the energy requirements of bat species with similar life-history traits resemble each other during periods of low energetic demands and are more similar to each other than to other small temperate mammals.     

Insular shifts and trade-offs in life-history traits in pond frogs in the Zhoushan Archipelago, China

Island and mainland populations of animal species often differ strikingly in life-history traits such as clutch size, egg size, total reproductive effort and body size. However, despite widespread recognition of insular shifts in these life-history traits in birds, mammals and reptiles, there have been no reports of such life-history shifts in amphibians. Furthermore, most studies have focused on one specific life-history trait without explicit consideration of coordinated evolution among these intimately linked life-history traits, and thus the relationships among these traits are poorly studied. Here we provide the first evidence of insular shifts and trade-offs in a coordinated suite of life-history traits for an amphibian species, the pond frog Rana nigromaculata . Life-history data were collected from eight islands in the Zhoushan Archipelago and neighboring mainland China. We found consistent, significant shifts in all life-history traits between mainland and island populations. Island populations had smaller clutch sizes, larger egg sizes, larger female body size and invested less in total reproductive effort than mainland populations. Significant negative relationships were found between egg size and clutch size and between egg size and total reproductive effort among frog populations after controlling for the effects of body size. Therefore, decreased reproductive effort and clutch size, larger egg size and body size in pond frogs on islands were selected through trade-offs as an overall life-history strategy. Our findings contribute to the formation of a broad, repeatable ecological generality for insular shifts in life-history traits across a range of terrestrial vertebrate taxa.     

Phylogenetic prediction of the maximum per capita rate of population growth

The maximum per capita rate of population growth, r, is a central measure of population biology. However, researchers can only directly calculate r when adequate time series, life tables and similar datasets are available. We instead view r as an evolvable, synthetic life-history trait and use comparative phylogenetic approaches to predict r for poorly known species. Combining molecular phylogenies, life-history trait data and stochastic macroevolutionary models, we predicted r for mammals of the Caniformia and Cervidae. Cross-validation analyses demonstrated that, even with sparse life-history data, comparative methods estimated r well and outperformed models based on body mass. Values of r predicted via comparative methods were in strong rank agreement with observed values and reduced mean prediction errors by approximately 68 per cent compared with two null models. We demonstrate the utility of our method by estimating r for 102 extant species in these mammal groups with unknown life-history traits.     

Strength of density feedback in census data increases from slow to fast life histories

Life-history theory predicts an increasing rate of population growth among species arranged along a continuum from slow to fast life histories. We examine the effects of this continuum on density-feedback strength estimated using long-term census data from >700 vertebrates, invertebrates, and plants. Four life-history traits (Age at first reproduction, Body size, Fertility, Longevity) were related statistically to Gompertz strength of density feedback using generalized linear mixed-effects models and multi-model inference. Life-history traits alone explained 10 to 30% of the variation in strength across species (after controlling for time-series length and phylogenetic nonindependence). Effect sizes were largest for body size in mammals and longevity in birds, and density feedback was consistently stronger for smaller-bodied and shorter-lived species. Overcompensatory density feedback (strength <-1) occurred in 20% of species, predominantly at the fast end of the life-history continuum, implying relatively high population variability. These results support the idea that life history leaves an evolutionary signal in long-term population trends as inferred from census data. Where there is a lack of detailed demographic data, broad life-history information can inform management and conservation decisions about rebound capacity from low numbers, and propensity to fluctuate, of arrays of species in areas planned for development, harvesting, protection, and population recovery.     

Evolution of embryonic developmental period in the marine bird families Alcidae and Spheniscidae: roles for nutrition and predation?

Background  

Nutrition and predation have been considered two primary agents of selection important in the evolution of avian life history traits. The relative importance of these natural selective forces in the evolution of avian embryonic developmental period (EDP) remain poorly resolved, perhaps in part because research has tended to focus on a single, high taxonomic-level group of birds: Order Passeriformes. The marine bird families Alcidae (auks) and Spheniscidae (penguins) exhibit marked variation in EDP, as well as behavioural and ecological traits ultimately linked to EDP. Therefore, auks and penguins provide a unique opportunity to assess the natural selective basis of variation in a key life-history trait at a low taxonomic-level. We used phylogenetic comparative methods to investigate the relative importance of behavioural and ecological factors related to nutrition and predation in the evolution of avian EDP.     

Unveiling the environmental drivers of intraspecific body size variation in terrestrial vertebrates

Aim

Whether intraspecific spatial patterns in body size are generalizable across species remains contentious, as well as the mechanisms underlying these patterns. Here we test several hypotheses explaining within-species body size variation in terrestrial vertebrates including the heat balance, seasonality, resource availability and water conservation hypotheses for ectotherms, and the heat conservation, heat dissipation, starvation resistance and resource availability hypotheses for endotherms.

Location

Global.

Time period

1970–2016.

Major taxa studied

Amphibians, reptiles, birds and mammals.

Methods

We collected 235,905 body size records for 2,229 species (amphibians = 36; reptiles = 81; birds = 1,545; mammals = 567) and performed a phylogenetic meta-analysis of intraspecific correlations between body size and environmental variables. We further tested whether correlations differ between migratory and non-migratory bird and mammal species, and between thermoregulating and thermoconforming ectotherms.

Results

For bird species, smaller intraspecific body size was associated with higher mean and maximum temperatures and lower resource seasonality. Size–environment relationships followed a similar pattern in resident and migratory birds, but the effect of resource availability on body size was slightly positive only for non-migratory birds. For mammals, we found that intraspecific body size was smaller with lower resource availability and seasonality, with this pattern being more evident in sedentary than migratory species. No clear size–environment relationships were found for reptiles and amphibians.

Main conclusions

Within-species body size variation across endotherms is explained by disparate underlying mechanisms for birds and mammals. Heat conservation (Bergmann's rule) and heat dissipation are the dominant processes explaining biogeographic intraspecific body size variation in birds, whereas in mammals, body size clines are mostly explained by the starvation resistance and resource availability hypotheses. Our findings contribute to a better understanding of the mechanisms behind species adaptations to the environment across their geographic distributions.     

Effect of species, root branching order and season on the root traits of 13 perennial grass species

Aims

Inter-specific comparisons of plant traits may vary depending on intra-specific variation. Here we examine the impact of root branching order and season on key functional root traits for grass species. We also compare root traits among co-existing grass species as a step towards defining root trait syndromes.

Methods

Monocultures of 13 grass species, grown under field conditions and subjected to intensive management, were used to record root trait values for coarse roots (1st order, >0.3?mm), fine roots (2nd and 3rd orders, <0.2?mm) and mixed root samples over three growing seasons.

Results

Branching order and species had a significant effect on root trait values, whereas season showed a marginal effect. The diameter of coarse roots was more variable than that of fine roots and, as expected, coarse roots had higher tissue density and lower specific root length values than fine roots. Principal component analysis run on eight root traits provided evidence for two trait syndromes related to resource acquisition and conservation strategies across grass species.

Conclusions

Our data show that root branching order is the main determinant of root trait variation among species. This highlights the necessity to include the proportion of fine vs coarse roots when measuring traits of mixed root samples.     

Greater variation in boreal plant community composition and community-level traits at local- than regional-scale

Questions

Rapid climate change in northern latitudes is expected to influence plant functional traits of the whole community (community-level traits) through species compositional changes and/or trait plasticity, limiting our ability to anticipate climate warming impacts on northern plant communities. We explored differences in plant community composition and community-level traits within and among four boreal peatland sites and determined whether intra- or interspecific variation drives community-level traits.

Location

Boreal biome of western North America.

Methods

We collected plant community composition and functional trait data along dominant topoedaphic and/or hydrologic gradients at four peatland sites spanning the latitudinal extent of the boreal biome of western North America. We characterized variability in community composition and community-level traits of understorey vascular and moss species both within (local-scale) and among sites (regional-scale).

Results

Against expectations, community-level traits of vascular plant and moss species were generally consistent among sites. Furthermore, interspecific variation was more important in explaining community-level trait variation than intraspecific variation. Within-site variation in both community-level traits and community composition was greater than among-site variation, suggesting that local environmental gradients (canopy density, organic layer thickness, etc.) may be more influential in determining plant community processes than regional-scale gradients.

Conclusions

Given the importance of interspecific variation to within-site shifts in community-level traits and greater variation of community composition within than among sites, we conclude that climate-induced shifts in understorey community composition may not have a strong influence on community-level traits in boreal peatlands unless local-scale environmental gradients are substantially altered.     

ProbABEL package for genome-wide association analysis of imputed data

Background  

Over the last few years, genome-wide association (GWA) studies became a tool of choice for the identification of loci associated with complex traits. Currently, imputed single nucleotide polymorphisms (SNP) data are frequently used in GWA analyzes. Correct analysis of imputed data calls for the implementation of specific methods which take genotype imputation uncertainty into account.     

Range contractions and reduced body mass predicted for endemic skinks of the Cameroon Volcanic Line with future warming

Aim

The vulnerability of montane species to environmental change has been increasingly recognized in recent years. However, most of these species are regionally endemic with restricted distributions, limiting dispersal necessary for avoiding extinction. The outcome of threats posed for montane lizards is further complicated in species exhibiting mass–temperature relationships where body size increases with cooling temperatures, and thus with altitude, causing intraspecific physiological and behavioural differences. We aimed to identify areas suitable for montane endemic skinks of the Cameroon Volcanic Line (CVL) under current and future climates to reveal patterns of persistence and vulnerability based on an intersection of climate and body mass.

Location

Cameroon Volcanic Line

Methods

We recorded occurrences and measured body mass in the field for two CVL-endemic skink species. We supplemented occurrences with online repository records. We projected current and future habitat suitability in the region by implementing bioclimatic species distribution models-based on occurrences. We tested for elevational variations in body mass, and integrated both occurrence and body mass information in a trait-based model to estimate current and future body mass.

Results

Projected currently suitable habitat for both species was limited to higher elevation regions, which are inhabited by numerous other threatened herpetofauna. We additionally detected Bergmann clines in body mass for both species. Given this variation in body mass, trait model projections covered slightly larger geographical ranges than bioclimatic estimates. Under future warming, both models project substantial contractions in suitable areas, potentially constraining species to mountain tops. Through the trait-based approach, we further detected potential warming-induced body mass reductions in projected suitable areas.

Main Conclusions

We demonstrate how combining occurrence records with species trait information in ecological modelling can reveal complementary trends for comprehensive warming impact assessments. Overall, challenges toward the persistence of CVL-endemic skinks should prompt urgent responses in national conservation management and local community engagement.