3D printed models are an accurate, cost-effective, and reproducible tool for quantifying terrestrial thermal environments

Predicting ecological responses to rapid environmental change has become one of the greatest challenges of modern biology. One of the major hurdles in forecasting these responses is accurately quantifying the thermal environments that organisms experience. The distribution of temperatures available within an organism's habitat is typically measured using data loggers called operative temperature models (OTMs) that are designed to mimic certain properties of heat exchange in the focal organism. The gold standard for OTM construction in studies of terrestrial ectotherms has been the use of copper electroforming which creates anatomically accurate models that equilibrate quickly to ambient thermal conditions. However, electroformed models require the use of caustic chemicals, are often brittle, and their production is expensive and time intensive. This has resulted in many researchers resorting to the use of simplified OTMs that can yield substantial measurement errors. 3D printing offers the prospect of robust, easily replicated, morphologically accurate, and cost-effective OTMs that capture the benefits but alleviate the problems associated with electroforming. Here, we validate the use of OTMs that were 3D printed using several materials across eight lizard species of different body sizes and living in habitats ranging from deserts to tropical forests. We show that 3D printed OTMs have low thermal inertia and predict the live animal's equilibration temperature with high accuracy across a wide range of body sizes and microhabitats. Finally, we developed a free online repository and database of 3D scans (https://www.3dotm.org/) to increase the accessibility of this tool to researchers around the world and facilitate ease of production of 3D printed models. 3D printing of OTMs is generalizable to taxa beyond lizards. If widely adopted, this approach promises greater accuracy and reproducibility in studies of terrestrial thermal ecology and should lead to improved forecasts of the biological impacts of climate change.

Variable stoichiometric and macronutrient responses to lizard predation in Ozark glade grasshopper communities

The General Stress Paradigm (GSP) predicts that prey body compositions should shift under chronic predation as prey increase body carbon and decrease body nitrogen content through dietary changes, heightened metabolism, reduced dietary efficiency, and the breakdown of nitrogen rich tissues to make labile carbohydrates available. In our study, we explored how the elemental and macronutrient content along with the morphology of three abundant Ozark glade grasshopper species differed between glades with and without predatory collared lizard (Crotaphytus collaris) populations. Our results indicated that lichen grasshoppers (Trimerotropis saxatilis) increased body C:N ratios in response to predators. Scudder's short-wing grasshoppers (Melanoplus scudderi) increased both body %C and %protein content, while the handsome grasshoppers (Syrbula admirabilis) did not significantly respond to the presence of collared lizards. None of the three grasshopper species showed morphological responses to predation. We also found that elemental and macronutrient content of grasshoppers was not always significantly correlated and was not associated with the same environmental factors, indicating a need to incorporate both perspectives in future research and utilize more accurate macromolecular assays. Overall, we found support for some aspects of the GSP in field-active animals and add to the growing body of evidence that predator-induced changes in prey body composition are more complex than predicted by the original GSP.

A chromosome-level genome assembly for the eastern fence lizard (Sceloporus undulatus), a reptile model for physiological and evolutionary ecology

Background

High-quality genomic resources facilitate investigations into behavioral ecology, morphological and physiological adaptations, and the evolution of genomic architecture. Lizards in the genus Sceloporus have a long history as important ecological, evolutionary, and physiological models, making them a valuable target for the development of genomic resources.

Findings

We present a high-quality chromosome-level reference genome assembly, SceUnd1.0 (using 10X Genomics Chromium, HiC, and Pacific Biosciences data), and tissue/developmental stage transcriptomes for the eastern fence lizard, Sceloporus undulatus. We performed synteny analysis with other snake and lizard assemblies to identify broad patterns of chromosome evolution including the fusion of micro- and macrochromosomes. We also used this new assembly to provide improved reference-based genome assemblies for 34 additional Sceloporus species. Finally, we used RNAseq and whole-genome resequencing data to compare 3 assemblies, each representing an increased level of cost and effort: Supernova Assembly with data from 10X Genomics Chromium, HiRise Assembly that added data from HiC, and PBJelly Assembly that added data from Pacific Biosciences sequencing. We found that the Supernova Assembly contained the full genome and was a suitable reference for RNAseq and single-nucleotide polymorphism calling, but the chromosome-level scaffolds provided by the addition of HiC data allowed synteny and whole-genome association mapping analyses. The subsequent addition of PacBio data doubled the contig N50 but provided negligible gains in scaffold length.

Conclusions

These new genomic resources provide valuable tools for advanced molecular analysis of an organism that has become a model in physiology and evolutionary ecology.

Functional and Environmental Constraints on Prey Capture Speed in a Lizard

Movement is an important component of animal behavior and determines how an organism interacts with its environment. The speed at which an animal moves through its environment can be constrained by internal (e.g., physiological state) and external factors (e.g., habitat complexity). When foraging, animals should move at speeds that maximize prey capture while minimizing mistakes (i.e., missing prey, slipping). We used experimental arenas containing obstacles spaced in different arrays to test how variation in habitat complexity influenced attack distance, prey capture speed, and foraging success in the Prairie Lizard. Obstacles spaced uniformly across arenas resulted in 15% slower prey capture speed and 30–38% shorter attack distance compared to arenas with no obstacles or with obstacles clustered in opposite corners of the arena. Prey capture probability was not influenced by arena type or capture speed, but declined with increasing attack distance. Similarly, the probability of prey consumption declined with attack distance across arena types. However, prey consumption probability declined with increasing prey capture speed in more open arenas but not in the cluttered arena. Foraging accuracy declined with increasing speed in more open arenas, and remained relatively constant when obstacles were in closer proximity. Foraging success was primarily constrained by intrinsic properties (speed-maneuverability tradeoff) when ample space was available, but environmental conditions had a greater impact on foraging success in “cluttered” habitats. This empirical test of theoretical predictions about optimal movement speeds in animals provides a step forward in understanding how animals select speeds in nature.

Role of phenotypic plasticity in morphological differentiation between watersnake populations

An individual's morphology is shaped by the environmental pressures it experiences, and the resulting morphological response is the culmination of both genetic factors and environmental (non-genetic) conditions experienced early in its life (i.e. phenotypic plasticity). The role that phenotypic plasticity plays in shaping phenotypes is important, but evidence for its influence is often mixed. We exposed female neonate diamond-backed watersnakes (Nerodia rhombifer) from populations experiencing different prey-size regimes to different feeding treatments to test the influence of phenotypic plasticity in shaping trophic morphology. We found that snakes in a large-prey treatment from a population frequently encountering large prey exhibited a higher growth rate in body size than individuals in a small-prey treatment from the same population. This pattern was not observed in snakes from a population that regularly encounters small prey. We also found that regardless of treatment, snakes from the small-prey population were smaller at birth than snakes from the large-prey population and remained so throughout the study. These results suggest that the ability to plastically respond to environmental pressures may be population-specific. These results also indicate a genetic predisposition towards larger body sizes in a population where large prey items are more common.

Covariation between Thermally Mediated Color and Performance Traits in a Lizard

Physiological changes in response to environmental cues are not uncommon. Temperature has strong, predictable effects on many traits, such that many traits in ectotherms follow stereotyped thermal performance curves in response to increasing temperature. The prairie lizard-an abundant lizard throughout the central United States-has thermally sensitive, blue abdominal and throat patches. Currently, the role of these patches is not well understood. In this study, we set out to investigate whether individual plasticity in patch color paralleled individual plasticity in sprint speed (do they covary), and if the plasticity in these two patches signal redundant or independent information, testing competing hypotheses suggested for the evolution of multiple signals. We found that both abdominal and throat patch hue follow classical thermal performance curves, suggesting that at the species level hue is a good predictor of sprinting ability. At the individual level, we found that color and performance were statistically repeatable, so individuals with relatively high phenotypic values maintain relatively high phenotypic values across all temperatures. Additionally, we found that abdominal and patch hue covary with sprinting speed at the individual level. Together, these results suggest that the bluest individuals are the fastest individuals across temperatures. However, we found that abdominal and throat patch hue do not covary with each other at the individual level, suggesting that these signals may have independent functions. The importance of examining the function of individual variation cannot be overstated, and overall, more work is needed to better understand both the proximate and ultimate mechanisms underlying signal plasticity in this species and others.

Variation in thermal tolerance of North American ants

Changing climates are predicted to alter the distribution of thermal niches. Small ectotherms such as ants may be particularly vulnerable to heat injury and death. We quantified the critical thermal maxima of 92 ant colonies representing 14 common temperate ant species. The mean CTmax for all measured ants was 47.8 °C (±0.27; range=40.2-51.2 °C), and within-colony variation was lower than among-colony variation. Critical thermal maxima differed among species and were negatively correlated with body size. Results of this study illustrate the importance of accounting for mass, among and within colony variation, and interspecific differences in diel activity patterns, which are often neglected in studies of ant thermal physiology.

Individual (co)variation in standard metabolic rate, feeding rate, and exploratory behavior in wild-caught semiaquatic salamanders

Repeatability is an important concept in evolutionary analyses because it provides information regarding the benefit of repeated measurements and, in most cases, a putative upper limit to heritability estimates. Repeatability (R) of different aspects of energy metabolism and behavior has been demonstrated in a variety of organisms over short and long time intervals. Recent research suggests that consistent individual differences in behavior and energy metabolism might covary. Here we present new data on the repeatability of body mass, standard metabolic rate (SMR), voluntary exploratory behavior, and feeding rate in a semiaquatic salamander and ask whether individual variation in behavioral traits is correlated with individual variation in metabolism on a whole-animal basis and after conditioning on body mass. All measured traits were repeatable, but the repeatability estimates ranged from very high for body mass (R = 0.98), to intermediate for SMR (R = 0.39) and food intake (R = 0.58), to low for exploratory behavior (R = 0.25). Moreover, repeatability estimates for all traits except body mass declined over time (i.e., from 3 to 9 wk), although this pattern could be a consequence of the relatively low sample size used in this study. Despite significant repeatability in all traits, we find little evidence that behaviors are correlated with SMR at the phenotypic and among-individual levels when conditioned on body mass. Specifically, the phenotypic correlations between SMR and exploratory behavior were negative in all trials but significantly so in one trial only. Salamanders in this study showed individual variation in how their exploratory behavior changed across trials (but not body mass, SMR, and feed intake), which might have contributed to observed changing correlations across trials.

The effects of temperature and activity on intraspecific scaling of metabolic rates in a lungless salamander

The scaling of metabolic rate with body mass holds substantial predictive power as many biological processes depend on energy. A significant body of theory has been developed based on the assumption that metabolic rate scales with body mass as a power function with an exponent of 0.75, and that this scaling relationship is independent of temperature. Here we test this hypothesis at the intraspecific level in a lungless salamander using data on both standard and maximal metabolic rates (SMR and MMR, respectively). We also address a recently proposed alternative explanation that predicts systematic variation in this mass-scaling exponent, the metabolic level boundaries hypothesis (MLB). Consistent with predictions of the metabolic theory of ecology the mass scaling of SMR and MMR were independent of temperature, however, we find evidence that the mass-scaling exponent for SMR and MMR differ significantly from 0.75. Further, our data do not provide strong support for MLB. Mass-scaling exponents for MMR generally exceed those for SMR, although these differences are rarely statistically significant.

Gambusia quadruncus (Cyprinodontiformes: Poeciliidae): a new species of mosquitofish from east-central México

Gambusia quadruncus n. sp., the llanos mosquitofish, is described from east-central México. The region inhabited by the species represents a hotspot of diversity of Gambusia, and G. quadruncus sometimes coexists with at least three congeners. The species differs from its closest relative, Gambusia affinis, in several characteristics with plausible effects on reproductive isolation, e.g. body size, body and fin morphology, male genital morphology (distal tip of gonopodium) and female anal spot morphology (colouration near the urogenital sinus). Moreover, combined analysis of mitochondrial and nuclear gene sequence data (c. 2158 total base pairs) indicates reciprocal monophyly of G. quadruncus and its sister species G. affinis, with levels of genetic divergence suggesting the two species diverged from one another over a million years ago. The origin of G. quadruncus may reflect a vicariant event associated with Pliocene orogenesis in the Tamaulipas Arch and a frontal section of the Sierra Madre Oriental (Lleran Mesas). Gambusia quadruncus inhabits a variety of freshwater habitats across several river drainages, with its range spanning at least 350 km from north to south, covering over 25 000 km(2). A key to aid identification of the species is provided.

ECOLOGICAL SPECIATION IN GAMBUSIA FISHES

Although theory indicates that natural selection can facilitate speciation as a by-product, demonstrating ongoing speciation via this by-product mechanism in nature has proven difficult. We examined morphological, molecular, and behavioral data to investigate ecology's role in incipient speciation for a post-Pleistocene radiation of Bahamas mosquitofish (Gambusia hubbsi) inhabiting blue holes. We show that adaptation to divergent predator regimes is driving ecological speciation as a by-product. Divergence in body shape, coupled with assortative mating for body shape, produces reproductive isolation that is twice as strong between populations inhabiting different predator regimes than between populations that evolved in similar ecological environments. Gathering analogous data on reproductive isolation at the interspecific level in the genus, we find that this mechanism of speciation may have been historically prevalent in Gambusia. These results suggest that speciation in nature can result as a by-product of divergence in ecologically important traits, producing interspecific patterns that persist long after speciation events have completed.

Population structure and history of a phenotypically variable teiid lizard (Ameiva chrysolaema) from Hispaniola: the influence of a geologically complex island

Ameiva chrysolaema is distributed across the island of Hispaniola in the West Indies. The species is restricted to dry lowlands between major mountain ranges and along the southern and eastern coasts. Phylogenetic and phylogeographic analyses of mtDNA sequence variation from 14 sampling localities identify at least three independent evolutionary lineages, separated from one another by major mountain ranges. Nested clade phylogeographic analysis (NCPA) suggests a complex history of population fragmentation, consistent with geological evidence of seawater incursions into the Azua and Enriquillo basins during the Pliocene/Pleistocene (approximately 1.6 mya). Significantly negative Fu's F(S) values and parameters of mismatch distributions suggest that formerly fragmented populations have recently expanded their ranges. Significantly large average population clade distances (APCD) for two sampling localities in the Azua basin suggest secondary contact at these localities of previously separated populations. The distribution of haplotypes among polymorphic populations of A. chrysolaema suggests that variation in dorsal pattern represents a polymorphism within evolutionary lineages. Ameiva leberi is ecologically indistinguishable from and syntopic with A. chrysolaema. Genetic data suggest that A. leberi is a junior synonym of A. chrysolaema.

Partial island submergence and speciation in an adaptive radiation: a multilocus analysis of the Cuban green anoles

Sympatric speciation is often proposed to account for species-rich adaptive radiations within lakes or islands, where barriers to gene flow or dispersal may be lacking. However, allopatric speciation may also occur in such situations, especially when ranges are fragmented by fluctuating water levels. We test the hypothesis that Miocene fragmentation of Cuba into three palaeo-archipelagos accompanied species-level divergence in the adaptive radiation of West Indian Anolis lizards. Analysis of morphology, mitochondrial DNA (mt DNA) and nuclear DNA in the Cuban green anoles (carolinensis subgroup) strongly supports three pre dictions made by this hypothesis. First, three geographical sets of populations, whose ranges correspond with palaeo-archipelago boundaries, are distinct and warrant recognition as independent evolutionary lineages or species. Coalescence of nuclear sequence fragments sampled from these species and the large divergences observed between their mtDNA haplotypes suggest separation prior to the subsequent unification of Cuba ca. 5 Myr ago. Second, molecular phylogenetic relationships among these species reflect historical geographical relationships rather than morphological similarity. Third, all three species remain distinct despite extensive geographical contact subsequent to island unification, occasional hybridization and introgression of mtDNA haplotypes. Allopatric speciation initiated during partial island submergence may play an important role in speciation during the adaptive radiation of Anolis lizards.