Morphological variation in salamanders and their potential response to climate change

Viviparity is associated with larger female size and higher sexual size dimorphism in a reproductively bimodal lizard

Squamate reptiles are central for studying phenotypic correlates of evolutionary transitions from oviparity to viviparity because these transitions are numerous, with many of them being recent. Several models of life-history theory predict that viviparity is associated with increased female size, and thus more female-biased sexual size dimorphism (SSD). Yet, the corresponding empirical evidence is overall weak and inconsistent. The lizard Zootoca vivipara, which occupies a major part of Northern Eurasia and includes four viviparous and two non-sister oviparous lineages, represents an excellent model for testing these predictions. We analysed how sex-specific body size and SSD is associated with parity mode, using body length data for nearly 14,000 adult individuals from 97 geographically distinct populations, which cover almost the entire species' range and represent all six lineages. Our analyses controlled for lineage identity, climatic seasonality (the strongest predictor of geographic body size variation in previous studies of this species) and several aspects of data heterogeneity. Parity mode, lineage and seasonality are significantly associated with female size and SSD; the first two predictors accounted for 14%-26% of the total variation each, while seasonality explained 5%-7%. Viviparous populations exhibited a larger female size than oviparous populations, with no concomitant differences in male size. The variation of male size was overall low and poorly explained by our predictors. Albeit fully expected from theory, the strong female bias of the body size differences between oviparous and viviparous populations found in Z. vivipara is not evident from available data on three other lizard systems of closely related lineages differing in parity mode. We confront this pattern with the data on female reproductive traits in the considered systems and the frequencies of evolutionary changes of parity mode in the corresponding lizard families and speculate why the life-history correlates of live-bearing in Z. vivipara are distinct. Comparing conspecific populations, our study provides the most direct evidence for the predicted effect of parity mode on adult body size but also demonstrates that the revealed pattern may not be general. This might explain why across squamates, viviparity is only weakly associated with larger size.

Deep ecomorphological and genetic divergence in Steller's Jays (Cyanocitta stelleri, Aves: Corvidae)

The relationship between ecology and morphology is a cornerstone of evolutionary biology, and quantifying variation across environments can shed light on processes that give rise to biodiversity. Three morphotypes of the Steller's Jay (Cyanocitta stelleri) occupy different ecoregions in western North America, which vary in climate and landcover. These morphotypes (Coastal, Interior, Rocky Mountain) differ in size, plumage coloration, and head pattern. We sampled 1080 Steller's Jays from 68 populations (plus 11 outgroups) to address three main questions using data on morphology, plumage, genetics (mtDNA, microsatellites), and ecological niches: (1) How do phenotypic and genetic traits vary within and among populations, morphotypes, and ecoregions? (2) How do population-level differences in Steller's Jays compare with other sister species pairs of North American birds? (3) What can we infer about the population history of Steller's Jays in relation to past climates, paleoecology, and niche evolution? We found substantial morphological, genetic, and ecological differentiation among morphotypes. The greatest genetic divergence separated Coastal and Interior morphotypes from the Rocky Mountain morphotype, which was associated with warmer, drier, and more open habitats. Microsatellites revealed additional structure between Coastal and Interior groups. The deep mtDNA split between Coastal/Interior and Rocky Mountain lineages of Steller's Jay (ND2 ~ 7.8%) is older than most North American avian sister species and dates to approximately 4.3 mya. Interior and Rocky Mountain morphotypes contact across a narrow zone with steep clines in traits and reduced gene flow. The distribution of the three morphotypes coincides with divergent varieties of ponderosa pine and Douglas fir. Species distribution models support multiple glacial refugia for Steller's Jays. Our integrative dataset combined with extensive geographic sampling provides compelling evidence for recognizing at least two species of Steller's Jay.

Morphological Variation and Its Environmental Correlates in the Taihangshan Swelled-Vented Frog across the Qinling Mountains

Simple Summary

Amphibians have weak dispersal abilities and are sensitive to environmental changes, resulting in their disproportionately high risk of extinction, with many species’ populations rapidly declining. Therefore, it is critical for amphibian conservation to understand their adaptive potential by exploring how amphibians respond to environmental changes based on morphological variations. Our results showed that morphological traits of Feirana taihangnica significantly differed among ages. Along with the increase in annual mean temperature, snout-vent length showed an anti-hump trend, indicating no support for Bergmann’s rule. Mean ultraviolet-B of the highest and lowest months were positively and negatively correlated with head width, thigh length and tibia width, respectively. The present study can help understand the effects of environmental changes on morphological variations of this mountain frog species and its adaptive potential, providing important implications for species conservation.

Abstract

The Taihangshan swelled-vented frog (Feirana taihangnica), an endemic species to the Qinling Mountains, central China, has experienced a dramatic population decline over the last few decades. The aim of this work was to quantify morphological variation in F. taihangnica across the Qinling Mountains and examine environmental correlates of this variation of morphological traits. We implemented a hierarchical partitioning to estimate the independent contribution of each environmental variable on morphological variations. Temperature seasonality was the greatest contributor in variations of snout-vent length (SVL) and head width, and ultraviolet-B (UV-B) radiation of the lowest month was the most influential on both thigh length and tibia width. Then, we used generalized additive models to analyze the relationship between each environmental factor and morphological trait variations. Along the increasing of annual mean temperature, SVL decreased firstly and then increased, indicating no support for Bergmann’s rule. Furthermore, SVL was negatively correlated with annual precipitation, while positively with temperature seasonality. The mean UV-B of the highest and lowest months was positively and negatively correlated with head width, thigh length and tibia width, respectively. The results of this study help us to understand adaptive potential of this mountain frog species via morphological variations in the light of environmental changes.

Identifying climate refugia for high-elevation Alpine birds under current climate warming predictions

Identifying climate refugia is key to effective biodiversity conservation under a changing climate, especially for mountain-specialist species adapted to cold conditions and highly threatened by climate warming. We combined species distribution models (SDMs) with climate forecasts to identify climate refugia for high-elevation bird species (Lagopus muta, Anthus spinoletta, Prunella collaris, Montifringilla nivalis) in the European Alps, where the ecological effects of climate changes are particularly evident and predicted to intensify. We considered future (2041-2070) conditions (SSP585 scenario, four climate models) and identified three types of refugia: (1) in-situ refugia potentially suitable under both current and future climate conditions, ex-situ refugia suitable (2) only in the future according to all future conditions, or (3) under at least three out of four future conditions. SDMs were based on a very large, high-resolution occurrence dataset (2901-12,601 independent records for each species) collected by citizen scientists. SDMs were fitted using different algorithms, balancing statistical accuracy, ecological realism and predictive/extrapolation ability. We selected the most reliable ones based on consistency between training and testing data and extrapolation over distant areas. Future predictions revealed that all species (with the partial exception of A. spinoletta) will undergo a range contraction towards higher elevations, losing 17%-59% of their current range (larger losses in L. muta). We identified ~15,000 km² of the Alpine region as in-situ refugia for at least three species, of which 44% are currently designated as protected areas (PAs; 18%-66% among countries). Our findings highlight the usefulness of spatially accurate data collected by citizen scientists, and the importance of model testing by extrapolating over independent areas. Climate refugia, which are only partly included within the current PAs system, should be priority sites for the conservation of Alpine high-elevation species and habitats, where habitat degradation/alteration by human activities should be prevented to ensure future suitability for alpine species.

How to approach the study of syndromes in macroevolution and ecology

Syndromes, wherein multiple traits evolve convergently in response to a shared selective driver, form a central concept in ecology and evolution. Recent work has questioned the existence of some classic syndromes, such as pollination and seed dispersal syndromes. Here, we discuss some of the major issues that have afflicted research into syndromes in macroevolution and ecology. First, correlated evolution of traits and hypothesized selective drivers is often relied on as the only evidence for adaptation of those traits to those hypothesized drivers, without supporting evidence. Second, the selective driver is often inferred from a combination of traits without explicit testing. Third, researchers often measure traits that are easy for humans to observe rather than measuring traits that are suited to testing the hypothesis of adaptation. Finally, species are often chosen for study because of their striking phenotypes, which leads to the illusion of syndromes and divergence. We argue that these issues can be avoided by combining studies of trait variation across entire clades or communities with explicit tests of adaptive hypotheses and that taking this approach will lead to a better understanding of syndrome-like evolution and its drivers.

Physical and ecological isolation contribute to maintain genetic differentiation between fire salamander subspecies

Landscape features shape patterns of gene flow among populations, ultimately determining where taxa lay along the continuum between panmixia to complete reproductive isolation. Gene flow can be restricted, leading to population differentiation in two non-exclusive ways: "physical isolation", in which geographic distance in combination with the landscape features restricts movement of individuals promoting genetic drift, and "ecological isolation", in which adaptive mechanisms constrain gene flow between different environments via divergent natural selection. In central Iberia, two fire salamander subspecies occur in parapatry across elevation gradients along the Iberian Central System mountains, while in the adjacent Montes de Toledo Region only one of them occurs. By integrating population and landscape genetic analyses, we show a ubiquitous role of physical isolation between and within mountain ranges, with unsuitable landscapes increasing differentiation between populations. However, across the Iberian Central System, we found strong support for a significant contribution of ecological isolation, with low genetic differentiation in environmentally homogeneous areas, but high differentiation across sharp transitions in precipitation seasonality. These patterns are consistent with a significant contribution of ecological isolation in restricting gene flow among subspecies. Overall, our results suggest that ecological divergence contributes to reduce genetic admixture, creating an opportunity for lineages to follow distinct evolutionary trajectories.

Capturing response differences of species distribution to climate and human pressures by incorporating local adaptation: Implications for the conservation of a critically endangered species

Considering local adaptation has been increasingly involved in forecasting species distributions under climate change and the management of species conservation. Herein, we take the critically endangered Chinese giant salamander (Andrias davidianus) that has both a low dispersal ability and distinct population divergence in different regions as an example. Basin-scale models that represent different populations in the Huanghe River Basin (HRB), the Yangtze River Basin (YRB), and the Pearl River Basin (PRB) were established using ensemble species distribution models. The species ranges under the future human population density (HPD) and climate change were predicted, and the range loss was evaluated for local basins in 2050 and 2070. Our results showed that the predominant factors affecting species distributions differed among basins, and the responses of the species occurrence to HPD and climate factors were distinctly different from northern to southern basins. Future HPD changes would be the most influential factor that engenders negative impacts on the species distribution in all three basins, especially in the HRB. Climate change will likely be less prominent in decreasing the species range, excluding in the YRB and PRB under the highest-emissions scenario in 2050. Overall, the high-emissions scenario would more significantly aggravate the negative impacts produced by HPD change in both 2050 and 2070, with maximum losses of species ranges in the HRB, YRB, and PRB of 83.4%, 60.0%, and 53.5%, respectively, under the scenarios of the combined impacts of HPD and climate changes. We proposed adapted conservation policies to effectively protect the habitat of this critically endangered animal in different basins based on the outcomes. Our research addresses the importance of incorporating local adaptation into species distribution modeling to inform conservation and management decisions.

Photographic database of the European cave salamanders, genus Hydromantes

European Hydromantes are a group of eight salamander species often occurring in subterranean habitats, which are a difficult environment to explore. All Hydromantes are strictly protected species and thus, low-impact methodologies to study these salamanders are strongly needed. Here we used a photographic technique to produce a large dataset of European Hydromantes, providing standardised pictures of 1,052 individuals belonging to the eight species, including hybrids as well. With our methodology we were able to reduce the handling time of individuals, and produce high quality pictures useful to investigate multiple life traits of these endangered species. Furthermore, the standardised photos provided here can be used for future comparisons of individuals from the surveyed populations.

Variation in body size and sexual size dimorphism in the most widely ranging lizard: testing the effects of reproductive mode and climate

Reproductive mode, ancestry, and climate are hypothesized to determine body size variation in reptiles but their effects have rarely been estimated simultaneously, especially at the intraspecific level. The common lizard (Zootoca vivipara) occupies almost the entire Northern Eurasia and includes viviparous and oviparous lineages, thus representing an excellent model for such studies. Using body length data for >10,000 individuals from 72 geographically distinct populations over the species' range, we analyzed how sex-specific adult body size and sexual size dimorphism (SSD) is associated with reproductive mode, lineage identity, and several climatic variables. Variation in male size was low and poorly explained by our predictors. In contrast, female size and SSD varied considerably, demonstrating significant effects of reproductive mode and particularly seasonality. Populations of the western oviparous lineage (northern Spain, south-western France) exhibited a smaller female size and less female-biased SSD than those of the western viviparous (France to Eastern Europe) and the eastern viviparous (Eastern Europe to Far East) lineages; this pattern persisted even after controlling for climatic effects. The phenotypic response to seasonality was complex: across the lineages, as well as within the eastern viviparous lineage, female size and SSD increase with increasing seasonality, whereas the western viviparous lineage followed the opposing trends. Altogether, viviparous populations seem to follow a saw-tooth geographic cline, which might reflect the nonmonotonic relationship of body size at maturity in females with the length of activity season. This relationship is predicted to arise in perennial ectotherms as a response to environmental constraints caused by seasonality of growth and reproduction. The SSD allometry followed the converse of Rensch's rule, a rare pattern for amniotes. Our results provide the first evidence of opposing body size-climate relationships in intraspecific units.

Continental-scale determinants of population trends in European amphibians and reptiles

The continuous decline of biodiversity is determined by the complex and joint effects of multiple environmental drivers. Still, a large part of past global change studies reporting and explaining biodiversity trends have focused on a single driver. Therefore, we are often unable to attribute biodiversity changes to different drivers, since a multivariable design is required to disentangle joint effects and interactions. In this work, we used a meta-regression within a Bayesian framework to analyze 843 time series of population abundance from 17 European amphibian and reptile species over the last 45 years. We investigated the relative effects of climate change, alien species, habitat availability, and habitat change in driving trends of population abundance over time, and evaluated how the importance of these factors differs across species. A large number of populations (54%) declined, but differences between species were strong, with some species showing positive trends. Populations declined more often in areas with a high number of alien species, and in areas where climate change has caused loss of suitability. Habitat features showed small variation over the last 25 years, with an average loss of suitable habitat of 0.1%/year per population. Still, a strong interaction between habitat availability and the richness of alien species indicated that the negative impact of alien species was particularly strong for populations living in landscapes with less suitable habitat. Furthermore, when excluding the two commonest species, habitat loss was the main correlate of negative population trends for the remaining species. By analyzing trends for multiple species across a broad spatial scale, we identify alien species, climate change, and habitat changes as the major drivers of European amphibian and reptile decline.

Rapid phenotypic evolution following shifts in life cycle complexity

Life cycle strategies have evolved extensively throughout the history of metazoans. The expression of disparate life stages within a single ontogeny can present conflicts to trait evolution, and therefore may have played a major role in shaping metazoan forms. However, few studies have examined the consequences of adding or subtracting life stages on patterns of trait evolution. By analysing trait evolution in a clade of closely related salamander lineages we show that shifts in the number of life cycle stages are associated with rapid phenotypic evolution. Specifically, salamanders with an aquatic-only (paedomorphic) life cycle have frequently added vertebrae to their trunk skeleton compared with closely related lineages with a complex aquatic-to-terrestrial (biphasic) life cycle. The rate of vertebral column evolution is also substantially lower in biphasic lineages, which may reflect the functional compromise of a complex cycle. This study demonstrates that the consequences of life cycle evolution can be detected at very fine scales of divergence. Rapid evolutionary responses can result from shifts in selective regimes following changes in life cycle complexity.

Little effect of climate change on body size of herbivorous beetles

Ongoing climate change affects various aspects of an animal's life, with important effects on distribution range and phenology. The relationship between global warming and body size changes in mammals and birds has been widely studied, with most findings indicating a decline in body size over time. Nevertheless, little data exist on similar size change patterns of invertebrates in general and insects in particular, and it is unclear whether insects should decrease in size or not with climate warming. We measured over 4000 beetle specimens, belonging to 29 beetle species in 8 families, collected in Israel during the last 100 years. The sampled species are all herbivorous. We examined whether beetle body size had changed over the years, while also investigating the relationships between body size and annual temperature, precipitation, net primary productivity (NPP) at the collection site and collection month. None of the environmental variables, including the collection year, was correlated with the size of most of the studied beetle species, while there were strong interactions of all variables with species. Our results, though mostly negative, suggest that the effect of climate change on insect body size is species-specific and by no means a general macro-ecological rule. They also suggest that the intrapopulation variance in body size of insects collected as adults in the field is large enough to conceal intersite environmental effects on body size, such as the effect of temperature and NPP.

Evidence for complex life cycle constraints on salamander body form diversification

Metazoans display a tremendous diversity of developmental patterns, including complex life cycles composed of morphologically disparate stages. In this regard, the evolution of life cycle complexity promotes phenotypic diversity. However, correlations between life cycle stages can constrain the evolution of some structures and functions. Despite the potential macroevolutionary consequences, few studies have tested the impacts of life cycle evolution on broad-scale patterns of trait diversification. Here we show that larval and adult salamanders with a simple, aquatic-only (paedomorphic) life cycle had an increased rate of vertebral column and body form diversification compared to lineages with a complex, aquatic-terrestrial (biphasic) life cycle. These differences in life cycle complexity explain the variations in vertebral number and adult body form better than larval ecology. In addition, we found that lineages with a simple terrestrial-only (direct developing) life cycle also had a higher rate of adult body form evolution than biphasic lineages, but still 10-fold lower than aquatic-only lineages. Our analyses demonstrate that prominent shifts in phenotypic evolution can follow long-term transitions in life cycle complexity, which may reflect underlying stage-dependent constraints.