Male mate preference as an agent of fecundity selection in a polymorphic salamander

Color polymorphisms are associated with variation in other traits which may affect individual fitness, and these color‐trait associations are expected to contribute to nonrandom mating in polymorphic species. The red‐backed salamander (Plethodon cinereus) exhibits a polymorphism in dorsal pattern: striped and unstriped, and previous studies have suggested that they may mate nonrandomly. However, the mechanism(s) contributing to this behavior remain unclear. Here we consider the role that male preference may have in driving mating behavior in P. cinereus. We limit our focus to striped individuals because this morph is most likely to be choosy given their dominant, aggressive behavioral profiles relative to unstriped males. Specifically, we evaluated (a) whether striped males preferentially associate with females with respect to her dorsum color, size, and body condition and (b) if so, whether female traits are evaluated via visual or chemical cues. We also considered whether the frequency of another male social behavior, nose taps, was associated with mate preferences. We found that striped male P. cinereus nose tapped more often to preferred females. However, males only assessed potential mates via chemical cues, preferring larger females overall. Reproductive phenology data on a sample of gravid females drawn from the same population indicated that the color morphs do not differ in reproductive traits, but larger females have greater fecundity. Given our findings, we conclude that female P. cinereus are under fecundity selection, mediated by male preference. In this manner, male mating behavior contributes to observations of nonrandom mate associations in this population of P. cinereus.     

Methodological considerations for detection of terrestrial small‐body salamander eDNA and implications for biodiversity conservation

Environmental DNA (eDNA) can be used as an assessment tool to detect populations of threatened species and provide fine‐scale data required to make management decisions. The objectives of this project were to use quantitative PCR (qPCR) to: (i) detect spiked salamander DNA in soil, (ii) quantify eDNA degradation over time, (iii) determine detectability of salamander eDNA in a terrestrial environment using soil, faeces, and skin swabs, (iv) detect salamander eDNA in a mesocosm experiment. Salamander eDNA was positively detected in 100% of skin swabs and 66% of faecal samples and concentrations did not differ between the two sources. However, eDNA was not detected in soil samples collected from directly underneath wild‐caught living salamanders. Salamander genomic DNA (gDNA) was detected in all qPCR reactions when spiked into soil at 10.0, 5.0, and 1.0 ng/g soil and spike concentration had a significant effect on detected concentrations. Only 33% of samples showed recoverable eDNA when spiked with 0.25 ng/g soil, which was the low end of eDNA detection. To determine the rate of eDNA degradation, gDNA (1 ng/g soil) was spiked into soil and quantified over seven days. Salamander eDNA concentrations decreased across days, but eDNA was still amplifiable at day 7. Salamander eDNA was detected in two of 182 mesocosm soil samples over 12 weeks (n = 52 control samples; n = 65 presence samples; n = 65 eviction samples). The discrepancy in detection success between experiments indicates the potential challenges for this method to be used as a monitoring technique for small‐bodied wild terrestrial salamander populations.     

Questioning the use of an amphibian colour morph as an indicator of climate change

The effects of recent climate changes on earth ecosystems are likely among the most important ecological concerns in human history. Good bioindicators are essential to properly assess the magnitude of these changes. In the last decades, studies have suggested that the morph proportion of the eastern red‐backed salamander (Plethodon cinereus), one of the most widely distributed and abundant vertebrate species in forests of eastern North America, could be used as a proxy for monitoring climate changes. Based on new discoveries in the northern areas of the species' range and on one of the largest compilation ever made for a vertebrate in North America (236 109 observations compiled from 1880 to 2013 in 1148 localities), we demonstrate however that climatic and geographic variables do not influence the colour morph proportions in P. cinereus populations. Consequently, we show that the use of colour morph proportions of this species do not perform as an indicator of climate change. Our findings indicate that bioindicator paradigms can be significantly challenged by new ecological research and more representative databases.     

Ecological resistance surfaces predict fine‐scale genetic differentiation in a terrestrial woodland salamander

Landscape genetics has seen tremendous advances since its introduction, but parameterization and optimization of resistance surfaces still poses significant challenges. Despite increased availability and resolution of spatial data, few studies have integrated empirical data to directly represent ecological processes as genetic resistance surfaces. In our study, we determine the landscape and ecological factors affecting gene flow in the western slimy salamander (Plethodon albagula). We used field data to derive resistance surfaces representing salamander abundance and rate of water loss through combinations of canopy cover, topographic wetness, topographic position, solar exposure and distance from ravine. These ecologically explicit composite surfaces directly represent an ecological process or physiological limitation of our organism. Using generalized linear mixed‐effects models, we optimized resistance surfaces using a nonlinear optimization algorithm to minimize model AIC. We found clear support for the resistance surface representing the rate of water loss experienced by adult salamanders in the summer. Resistance was lowest at intermediate levels of water loss and higher when the rate of water loss was predicted to be low or high. This pattern may arise from the compensatory movement behaviour of salamanders through suboptimal habitat, but also reflects the physiological limitations of salamanders and their sensitivity to extreme environmental conditions. Our study demonstrates that composite representations of ecologically explicit processes can provide novel insight and can better explain genetic differentiation than ecologically implicit landscape resistance surfaces. Additionally, our study underscores the fact that spatial estimates of habitat suitability or abundance may not serve as adequate proxies for describing gene flow, as predicted abundance was a poor predictor of genetic differentiation.     

Demographic consequences of climate variation along an elevational gradient for a montane terrestrial salamander

Climate change represents a significant threat to amphibians. However, for many species, the relationship between demography and climate is unknown, which limits predictive models. Here, we describe the life history variation of Plethodon montanus using capture–recapture data over a period of 4 years, along an elevational gradient to determine how survival and growth vary with climate, and how these relationships vary with elevation. We used a hierarchical model to estimate asymptotic size and growth rate and a spatial Cormack-Jolly-Seber model to estimate the probability of capture and survival and dispersal variance. We found that during the active season, growth and survival rates are both positively affected by precipitation; however, while survival was positively affected by temperature at all elevations, higher temperatures led to a decrease in growth at lower elevations, while at higher elevations the opposite was true. During the inactive season, we found reduced growth rates, whereas survival was lower compared with the active at lower elevations but was higher at higher elevations. Increased inactive season temperatures resulted in decreased survival while we found that temperature, amount of snow, and elevation interacted to influence survival. At low elevations, which were warmer, survival decreased with increasing snowfall but at higher elevations, survival generally increased with increasing snowfall. Our results demonstrate that understanding how the environment can affect salamander demography to develop mechanistic models will require knowledge of the actual environmental conditions experienced by a given population as well as an understanding of the overall differences in climate at a given site.     

Evolutionary response to global change: Climate and land use interact to shape color polymorphism in a woodland salamander

Evolutionary change has been demonstrated to occur rapidly in human‐modified systems, yet understanding how multiple components of global change interact to affect adaptive evolution remains a critical knowledge gap. Climate change is predicted to impose directional selection on traits to reduce thermal stress, but the strength of directional selection may be mediated by changes in the thermal environment driven by land use. We examined how regional climatic conditions and land use interact to affect genetically based color polymorphism in the eastern red‐backed salamander (Plethodon cinereus). P. cinereus is a woodland salamander with two primary discrete color morphs (striped, unstriped) that have been associated with macroclimatic conditions. Striped individuals are most common in colder regions, but morph frequencies can be variable within climate zones. We used path analysis to analyze morph frequencies among 238,591 individual salamanders across 1,170 sites in North America. Frequency of striped individuals was positively related to forest cover in populations occurring in warmer regions (>7°C annually), a relationship that was weak to nonexistent in populations located in colder regions (≤7°C annually). Our results suggest that directional selection imposed by climate warming at a regional scale may be amplified by forest loss and suppressed by forest persistence, with a mediating effect of land use that varies geographically. Our work highlights how the complex interaction of selection pressures imposed by different components of global change may lead to divergent evolutionary trajectories among populations.     

Evaluating within‐population variability in behavior and demography for the adaptive potential of a dispersal‐limited species to climate change

Multiple pathways exist for species to respond to changing climates. However, responses of dispersal‐limited species will be more strongly tied to ability to adapt within existing populations as rates of environmental change will likely exceed movement rates. Here, we assess adaptive capacity in Plethodon cinereus, a dispersal‐limited woodland salamander. We quantify plasticity in behavior and variation in demography to observed variation in environmental variables over a 5‐year period. We found strong evidence that temperature and rainfall influence P. cinereus surface presence, indicating changes in climate are likely to affect seasonal activity patterns. We also found that warmer summer temperatures reduced individual growth rates into the autumn, which is likely to have negative demographic consequences. Reduced growth rates may delay reproductive maturity and lead to reductions in size‐specific fecundity, potentially reducing population‐level persistence. To better understand within‐population variability in responses, we examined differences between two common color morphs. Previous evidence suggests that the color polymorphism may be linked to physiological differences in heat and moisture tolerance. We found only moderate support for morph‐specific differences for the relationship between individual growth and temperature. Measuring environmental sensitivity to climatic variability is the first step in predicting species' responses to climate change. Our results suggest phenological shifts and changes in growth rates are likely responses under scenarios where further warming occurs, and we discuss possible adaptive strategies for resulting selective pressures.     

Modeling optimal responses and fitness consequences in a changing Arctic

Animals must balance a series of costs and benefits while trying to maximize their fitness. For example, an individual may need to choose how much energy to allocate to reproduction versus growth, or how much time to spend on vigilance versus foraging. Their decisions depend on complex interactions between environmental conditions, behavioral plasticity, reproductive biology, and energetic demands. As animals respond to novel environmental conditions caused by climate change, the optimal decisions may shift. Stochastic dynamic programming provides a flexible modeling framework with which to explore these trade‐offs, but this method has not yet been used to study possible changes in optimal trade‐offs caused by climate change. We created a stochastic dynamic programming model capturing trade‐off decisions required by an individual adult female polar bear (Ursus maritimus) as well as the fitness consequences of her decisions. We predicted optimal foraging decisions throughout her lifetime as well as the energetic thresholds below which it is optimal for her to abandon a reproductive attempt. To explore the effects of climate change, we shortened the spring feeding period by up to 3 weeks, which led to predictions of riskier foraging behavior and higher reproductive thresholds. The resulting changes in fitness may be interpreted as a best‐case scenario, where bears adapt instantaneously and optimally to new environmental conditions. If the spring feeding period was reduced by 1 week, her expected fitness declined by 15%, and if reduced by 3 weeks, expected fitness declined by 68%. This demonstrates an effective way to explore a species' optimal response to a changing landscape of costs and benefits and highlights the fact that small annual effects can result in large cumulative changes in expected lifetime fitness.     

Recent spatial and temporal changes in body size of terrestrial vertebrates: probable causes and pitfalls

Geographical and temporal variations in body size are common phenomena among organisms and may evolve within a few years. We argue that body size acts much like a barometer, fluctuating in parallel with changes in the relevant key predictor(s), and that geographical and temporal changes in body size are actually manifestations of the same drivers. Frequently, the principal predictors of body size are food availability during the period of growth and ambient temperature, which often affects food availability. Food availability depends on net primary productivity that, in turn, is determined by climate and weather (mainly temperature and precipitation), and these depend mainly on solar radiation and other solar activities. When the above predictors are related to latitude the changes have often been interpreted as conforming to Bergmann's rule, but in many cases such interpretations should be viewed with caution due to the interrelationships among various environmental predictors. Recent temporal changes in body size have often been related to global warming. However, in many cases the above key predictors are not related to either latitude and/or year, and it is the task of the researcher to determine which particular environmental predictor is the one that determines food availability and, in turn, body size. The chance of discerning a significant change in body size depends to a large extent on sample size (specimens/year). The most recent changes in body size are probably phenotypic, but there are some cases in which they are partly genetic.     

Ecological separation in a polymorphic terrestrial salamander

1. When studying speciation, researchers commonly examine reproductive isolation in recently diverged populations. Polymorphic species provide an opportunity to examine the role of reproductive isolation in populations that may be in the process of divergence. 2. We examined a polymorphic population of Plethodon cinereus (red-backed salamanders) for evidence of sympatric ecological separation by colour morphology. Recent studies have correlated temperature and climate with colour morphology in this species, but no studies have looked at differences in diet or mate choice between colour morphs. We used artificial cover objects to assess salamander diet, mating preference and surface activity over a 2-year period at a field site in north-eastern Ohio. 3. We detected differences in diet between two colour morphs, striped and unstriped. The diets of striped individuals were significantly more diverse and were made up of more profitable prey than the diets of unstriped salamanders. 4. Opposite sex pairs were made up of individuals of the same colour morph and striped males were found more often with larger females than were unstriped males. 5. We corroborate findings of earlier studies suggesting that the unstriped form is adapted to warmer conditions. Unstriped individuals were the first to withdraw from the forest floor as temperatures fell in the late fall. We found no evidence that the colour morphs responded differently to abiotic factors such as soil moisture and relative humidity, and responses to surface temperatures were also equivocal. 6. We conclude that the two colour morphs exhibit some degree of ecological separation and tend to mate assortatively, but are unlikely to be undergoing divergence given the observed frequency of intermorph pairings.     

Genetic variation in Plethodon cinereus and Plethodon hubrichti from in and around a contact zone

Climate change poses several challenges to biological communities including changes in the frequency of encounters between closely related congeners as a result of range shifts. When climate change leads to increased hybridization, hybrid dysfunction or genetic swamping may increase extinction risk—particularly in range‐restricted species with low vagility. The Peaks of Otter Salamander, Plethodon hubrichti, is a fully terrestrial woodland salamander that is restricted to ~18 km of ridgeline in the mountains of southwestern Virginia, and its range is surrounded by the abundant and widespread Eastern Red‐backed Salamander, Plethodon cinereus. In order to determine whether these two species are hybridizing and how their range limits may be shifting, we assessed variation at eight microsatellite loci and a 1,008 bp region of Cytochrome B in both species at allopatric reference sites and within a contact zone. Our results show that hybridization between P. hubrichti and P. cinereus either does not occur or is very rare. However, we find that diversity and differentiation are substantially higher in the mountaintop endemic P. hubrichti than in the widespread P. cinereus, despite similar movement ability for the two species as assessed by a homing experiment. Furthermore, estimation of divergence times between reference and contact zone populations via approximate Bayesian computation is consistent with the idea that P. cinereus has expanded into the range of P. hubrichti. Given the apparent recent colonization of the contact zone by P. cinereus, future monitoring of P. cinereus range limits should be a priority for the management of P. hubrichti populations.     

Bison body size and climate change

The relationship between body size and temperature of mammals is poorly resolved, especially for large keystone species such as bison (Bison bison). Bison are well represented in the fossil record across North America, which provides an opportunity to relate body size to climate within a species. We measured the length of a leg bone (calcaneal tuber, DstL) in 849 specimens from 60 localities that were dated by stratigraphy and ¹⁴C decay. We estimated body mass (M) as M = (DstL/11.49)³. Average annual temperature was estimated from δ¹⁸O values in the ice cores from Greenland. Calcaneal tuber length of Bison declined over the last 40,000 years, that is, average body mass was 37% larger (910 ± 50 kg) than today (665 ± 21 kg). Average annual temperature has warmed by 6°C since the Last Glacial Maximum (~24–18 kya) and is predicted to further increase by 4°C by the end of the 21st century. If body size continues to linearly respond to global temperature, Bison body mass will likely decline by an additional 46%, to 357 ± 54 kg, with an increase of 4°C globally. The rate of mass loss is 41 ± 10 kg per°C increase in global temperature. Changes in body size of Bison may be a result of migration, disease, or human harvest but those effects are likely to be local and short‐term and not likely to persist over the long time scale of the fossil record. The strong correspondence between body size of bison and air temperature is more likely the result of persistent effects on the ability to grow and the consequences of sustaining a large body mass in a warming environment. Continuing rises in global temperature will likely depress body sizes of bison, and perhaps other large grazers, without human intervention.     

Contrasting effects of climate on juvenile body size in a Southern Hemisphere passerine bird

Despite extensive research on the topic, it has been difficult to reach general conclusions as to the effects of climate change on morphology in wild animals: in particular, the effects of warming temperatures have been associated with increases, decreases or stasis in body size in different populations. Here, we use a fine‐scale analysis of associations between weather and offspring body size in a long‐term study of a wild passerine bird, the cooperatively breeding superb fairy‐wren, in south‐eastern Australia to show that such variation in the direction of associations occurs even within a population. Over the past 26 years, our study population has experienced increased temperatures, increased frequency of heatwaves and reduced rainfall – but the mean body mass of chicks has not changed. Despite the apparent stasis, mass was associated with weather across the previous year, but in multiple counteracting ways. Firstly, (i) chick mass was negatively associated with extremely recent heatwaves, but there also positive associations with (ii) higher maximum temperatures and (iii) higher rainfall, both occurring in a period prior to and during the nesting period, and finally (iv) a longer‐term negative association with higher maximum temperatures following the previous breeding season. Our results illustrate how a morphological trait may be affected by both short‐ and long‐term effects of the same weather variable at multiple times of the year and that these effects may act in different directions. We also show that climate within the relevant time windows may not be changing in the same way, such that overall long‐term temporal trends in body size may be minimal. Such complexity means that analytical approaches that search for a single ‘best’ window for one particular weather variable may miss other relevant information, and is also likely to make analyses of phenotypic plasticity and prediction of longer‐term population dynamics difficult.     

Sequence capture and next‐generation sequencing of ultraconserved elements in a large‐genome salamander

Amidst the rapid advancement in next‐generation sequencing (NGS) technology over the last few years, salamanders have been left behind. Salamanders have enormous genomes—up to 40 times the size of the human genome—and this poses challenges to generating NGS data sets of quality and quantity similar to those of other vertebrates. However, optimization of laboratory protocols is time‐consuming and often cost prohibitive, and continued omission of salamanders from novel phylogeographic research is detrimental to species facing decline. Here, we use a salamander endemic to the southeastern United States, Plethodon serratus, to test the utility of an established protocol for sequence capture of ultraconserved elements (UCEs) in resolving intraspecific phylogeographic relationships and delimiting cryptic species. Without modifying the standard laboratory protocol, we generated a data set consisting of over 600 million reads for 85 P. serratus samples. Species delimitation analyses support recognition of seven species within P. serratus sensu lato, and all phylogenetic relationships among the seven species are fully resolved under a coalescent model. Results also corroborate previous data suggesting nonmonophyly of the Ouachita and Louisiana regions. Our results demonstrate that established UCE protocols can successfully be used in phylogeographic studies of salamander species, providing a powerful tool for future research on evolutionary history of amphibians and other organisms with large genomes.     

Biology, not environment, drives major patterns in maximum tetrapod body size through time

Abiotic and biological factors have been hypothesized as controlling maximum body size of tetrapods and other animals through geological time. We analyse the effects of three abiotic factors--oxygen, temperature and land area--on maximum size of Permian-Jurassic archosauromorphs and therapsids, and Cenozoic mammals, using time series generalized least-squares regression models. We also examine maximum size growth curves for the Permian-Jurassic data by comparing fits of Gompertz and logistic models. When serial correlation is removed, we find no robust correlations, indicating that these environmental factors did not consistently control tetrapod maximum size. Gompertz models--i.e. exponentially decreasing rate of size increase at larger sizes--fit maximum size curves far better than logistic models. This suggests that biological limits such as reduced fecundity and niche space availability become increasingly limiting as larger sizes are reached. Environmental factors analysed may still have imposed an upper limit on tetrapod body size, but any environmentally imposed limit did not vary substantially during the intervals examined despite variation in these environmental factors.     

Proximate causes of altitudinal differences in body size in an agamid lizard

Body size is directly linked to key life history traits such as growth, fecundity, and survivorship. Identifying the causes of body size variation is a critical task in ecological and evolutionary research. Body size variation along altitudinal gradients has received considerable attention; however, the underlying mechanisms are poorly understood. Here, we compared the growth rate and age structure of toad‐headed lizards (Phrynocephalus vlangalii) from two populations found at different elevations in the Qinghai‐Tibetan Plateau. We used mark‐recapture and skeletochronological analysis to identify the potential proximate causes of altitudinal variation in body size. Lizards from the high‐elevation site had higher growth rates and attained slightly larger adult body sizes than lizards from the low‐elevation site. However, newborns produced by high‐elevation females were smaller than those by low‐elevation females. Von Bertalanffy growth estimates predicted high‐elevation individuals would reach sexual maturity at an earlier age and have a lower mean age than low‐elevation individuals. Relatively lower mean age for the high‐elevation population was confirmed using the skeletochronological analysis. These results support the prediction that a larger adult body size of high‐elevation P. vlangalii results from higher growth rates, associated with higher resource availability.     

Genomic data detect corresponding signatures of population size change on an ecological time scale in two salamander species

Understanding the demography of species over recent history (e.g. <100 years) is critical in studies of ecology and evolution, but records of population history are rarely available. Surveying genetic variation is a potential alternative to census‐based estimates of population size, and can yield insight into the demography of a population. However, to assess the performance of genetic methods, it is important to compare their estimates of population history to known demography. Here, we leveraged the exceptional resources from a wetland with 37 years of amphibian mark–recapture data to study the utility of genetically based demographic inference on salamander species with documented population declines (Ambystoma talpoideum) and expansions (A. opacum), patterns that have been shown to be correlated with changes in wetland hydroperiod. We generated ddRAD data from two temporally sampled populations of A. opacum (1993, 2013) and A. talpoideum (1984, 2011) and used coalescent‐based demographic inference to compare alternate evolutionary models. For both species, demographic model inference supported population size changes that corroborated mark–recapture data. Parameter estimation in A. talpoideum was robust to our variations in analytical approach, while estimates for A. opacum were highly inconsistent, tempering our confidence in detecting a demographic trend in this species. Overall, our robust results in A. talpoideum suggest that genome‐based demographic inference has utility on an ecological scale, but researchers should also be cognizant that these methods may not work in all systems and evolutionary scenarios. Demographic inference may be an important tool for population monitoring and conservation management planning.     

Disentangling the effects of spring anomalies in climate and net primary production on body size of temperate songbirds

Body size is implicated in individual fitness and population dynamics. Mounting interest is being given to the effects of environmental change on body size, but the underlying mechanisms are poorly understood. We tested whether body size and body condition are related to ambient temperature (heat maintenance hypothesis), or/and explained by variations in primary production (food availability hypothesis) during the period of body growth in songbirds. We also explored whether annual population‐level variations of mean body size are due to changes of juvenile growth and/or size‐dependent mortality during the first year. For 41 species, from 257 sites across France, we tested for relationships between wing length (n = 107 193) or body condition (n = 82 022) and local anomalies in temperature, precipitation and net primary production (NDVI) during the breeding period, for juveniles and adults separately. Juvenile body size was best explained by primary production: wings were longer in years with locally high NDVI, but not shorter in years with low NDVI. Temperature showed a slightly positive effect. Body condition and adult wing length did not covary with any of the other tested variables. We found no evidence of climate‐driven size‐dependent mortality for the breeding season. In our temperate system, local climatic anomalies explained little of the body size variation. A large part of wing length variance was site‐specific, suggesting that avian size was more dependent on local drivers than global ones. Net primary production influenced juvenile size the most through effects on body growth. We suggest that, during the breeding season in temperate systems, thermoregulatory mechanisms are less involved in juvenile growth than food assimilation.