Diversity and biogeography of South American mud turtles elucidated by multilocus DNA sequencing (Testudines: Kinosternidae).

Kinosternon is the most speciose genus of extant turtles, with 22 currently recognized species, distributed across large parts of the Americas. Most species have small distributions, but K. leucostomum and K. scorpioides range from Mexico to South America. Previous studies have found discordance between mitochondrial and nuclear phylogenies in some kinosternid groups, with the current taxonomy following the nuclear-based results. Herein, based on extended molecular, geographic, and taxonomic sampling, we explore the phylogeographic structure and taxonomic limits for K. leucostomum and the K. scorpioides group and present a fossil-calibrated nuclear time tree for Kinosternon. Our results reveal contrasting differentiation patterns for the K. scorpioides group and K. leucostomum, despite overlapping distributions. Kinosternon leucostomum shows only shallow geographic divergence, whereas the K. scorpioides group is polyphyletic with up to 10 distinct taxa, some of them undescribed. We support the elevation of K. s. albogulare and K. s. cruentatum to species level. Given the deep divergence within the genus Kinosternon, we propose the recognition of three subgenera, Kinosternon, Cryptochelys and Thyrosternum, and the abandonment of the group-based classification, at least for the K. leucostomum and K. scorpioides groups. Our results show an initial split in Kinosternon that gave rise to two main radiations, one Nearctic and one mainly Neotropical. Most speciation events in Kinosternon occurred during the Quaternary and we hypothesize that they were mediated by both climatic and geological events. Additionally, our data imply that at least three South American colonizations occurred, two in the K. leucostomum group, and one in the K. scorpioides group. Additionally, we hypothesize that discordance between mitochondrial and nuclear phylogenetic signal is due to mitochondrial capture from an extinct kinosternine lineage.

Mother knows best: nest-site choice homogenizes embryo thermal environments among populations in a widespread ectotherm.

Species with large geographical ranges provide an excellent model for studying how different populations respond to dissimilar local conditions, particularly with respect to variation in climate. Maternal effects, such as nest-site choice greatly affect offspring phenotypes and survival. Thus, maternal behaviour has the potential to mitigate the effects of divergent climatic conditions across a species' range. We delineated natural nesting areas of six populations of painted turtles (Chrysemys picta) that span a broad latitudinal range and quantified spatial and temporal variation in nest characteristics. To quantify microhabitats available for females to choose, we also identified sites within the nesting area of each location that were representative of available thermal microhabitats. Across the range, females nested non-randomly and targeted microhabitats that generally had less canopy cover and thus higher nest temperatures. Nest microhabitats differed among locations but did not predictably vary with latitude or historic mean air temperature during embryonic development. In conjunction with other studies of these populations, our results suggest that nest-site choice is homogenizing nest environments, which buffers embryos from thermally induced selection and could slow embryonic evolution. Thus, although effective at a macroclimatic scale, nest-site choice is unlikely to compensate for novel stressors that rapidly increase local temperatures. This article is part of the theme issue 'The evolutionary ecology of nests: a cross-taxon approach'.

Sex-specific innate immunity and ageing in long-lived fresh water turtles (Kinosternon flavescens: Kinosternidae).

BACKGROUND: The progressive deregulation of the immune system with age, termed immunosenescence, has been well studied in mammalian systems, but studies of immune function in long-lived, wild, non-mammalian populations are scarce. In this study we leverage a 38-year mark-recapture study to quantify the relationships among age, sex, survival, reproductive output and the innate immune system in a long-lived reptile, yellow mud turtles (Kinosternon flavescens; Testudines; Kinosternidae). METHODS: We estimated rates of survival and age-specific mortality by sex based on mark-recapture data for 1530 adult females and 860 adult males over 38 years of captures. We analyzed bactericidal competence (BC), and two immune responses to foreign red blood cells - natural antibody-mediated haemagglutination (NAbs), and complement-mediated haemolysis ability (Lys) - in 200 adults (102 females; 98 males) that ranged from 7 to 58 years of age captured in May 2018 during their emergence from brumation, and for which reproductive output and long-term mark-recapture data were available. RESULTS: We found that females are smaller and live longer than males in this population, but the rate of accelerating mortality across adulthood is the same for both sexes. In contrast, males exhibited higher innate immunity than females for all three immune variables we measured. All immune responses also varied inversely with age, indicating immunosenescence. For females that reproduced in the preceding reproductive season, egg mass (and therefore total clutch mass) increased with age,. In addition to immunosenescence of bactericidal competence, females that produced smaller clutches also had lower bactericidal competence. CONCLUSIONS: Contrary to the general vertebrate pattern of lower immune responses in males than females (possibly reflecting the suppressive effects of androgens), we found higher levels of all three immune variables in males. In addition, contrary to previous work that found no evidence of immunosenescence in painted turtles or red-eared slider turtles, we found a decrease in bactericidal competence, lysis ability, and natural antibodies with age in yellow mud turtles.

Blood chemistry and biliverdin differ according to reproduction and tourism in a free-living lizard.

Changes in the physiological health of species are an essential indicator of changing conditions and environmental challenges. Reponses to environmental challenges can often induce stress, influence physiology, and change metabolism in organisms. Here we tested blood chemistry parameters indicative of stress and metabolic activity using an i-STAT point-of-care blood analyzer in seven populations of free-ranging rock iguanas exposed to varying levels of tourism and supplemental feeding. We found significant differences in blood chemistry (glucose, oxygen, carbon dioxide, hematocrit, hemoglobin, calcium, potassium, and biliverdin levels) among populations exposed to varying levels of tourism, and some variation between sexes and reproductive states. However, different variables are not directly related to one another, suggesting that the causal physiological pathways driving tourism-induced differences are influenced by mechanisms that are not detected by common analyses of blood chemistry. Future work should investigate upstream regulators of these factors affected by tourism. Regardless, these blood metrics are known to be both stress sensitive and related to metabolic activity, suggesting that exposure to tourism and associated supplemental feeding by tourists are generally driven by stress-related changes in blood chemistry, biliverdin, and metabolism.

Phylogenomics and historical biogeography of West Indian Rock Iguanas (genus Cyclura).

The genus Cyclura includes nine extant species and six subspecies of West Indian Rock Iguanas and is one of the most imperiled genera of squamate reptiles globally. An understanding of species diversity, evolutionary relationships, diversification, and historical biogeography in this group is crucial for implementing sound long-term conservation strategies. We collected DNA samples from 1 to 10 individuals per taxon from all Cyclura taxa (n = 70 ingroup individuals), focusing where possible on incorporating individuals from different populations of each species. We also collected 1-2 individuals from each of seven outgroup species of iguanas (Iguana delicatissima; five Ctenosaura species) and Anolis sagrei (n = 12 total outgroup individuals). We used targeted genomic sequence capture to isolate and to sequence 1,872 loci comprising of 687,308 base pairs (bp) from each of the 82 individuals from across the nuclear genome. We extracted mitochondrial reads and assembled and annotated mitogenomes for all Cyclura taxa plus outgroup species. We present well-supported phylogenomic gene tree/species tree analyses for all extant species of Cyclura using ASTRAL-III, SVDQuartets, and StarBEAST2 methods, and discuss the taxonomic, biogeographic, and conservation implications of these data. We find a most recent common ancestor of the genus 9.91 million years ago. The earliest divergence within Cyclura separates C. pinguis from a clade comprising all other Cyclura. Within the latter group, a clade comprising C. carinata from the southern Lucayan Islands and C. ricordii from Hispaniola is the sister taxon to a clade comprising the other Cyclura. Among the other Cyclura, the species C. cornuta and C. stejnegeri (from Hispaniola and Isla Mona) form the sister taxon to a clade of species from Jamaica (C. collei), Cuba and Cayman Islands (C. nubila and C. lewisi), and the eastern (C. rileyi) and western (C. cychlura) Lucayan Islands. Cyclura cychlura and C. rileyi form a clade whose sister taxa are C. nubila and C. lewisi. Cyclura collei is the sister taxon to these four species combined.

Diverse aging rates in ectothermic tetrapods provide insights for the evolution of aging and longevity.

Comparative studies of mortality in the wild are necessary to understand the evolution of aging; yet, ectothermic tetrapods are underrepresented in this comparative landscape, despite their suitability for testing evolutionary hypotheses. We present a study of aging rates and longevity across wild tetrapod ectotherms, using data from 107 populations (77 species) of nonavian reptiles and amphibians. We test hypotheses of how thermoregulatory mode, environmental temperature, protective phenotypes, and pace of life history contribute to demographic aging. Controlling for phylogeny and body size, ectotherms display a higher diversity of aging rates compared with endotherms and include phylogenetically widespread evidence of negligible aging. Protective phenotypes and life-history strategies further explain macroevolutionary patterns of aging. Analyzing ectothermic tetrapods in a comparative context enhances our understanding of the evolution of aging.

Complex tourism and season interactions contribute to disparate physiologies in an endangered rock iguana.

To promote survival and fitness, organisms use a suite of physiological systems to respond to both predictable and unpredictable changes in the environment. These physiological responses are also influenced by changes in life history state. The continued activation of physiological systems stemming from persistent environmental perturbations enable animals to cope with these challenges but may over time lead to significant effects on the health of wildlife. In the present study, we tested how varying environmental perturbations driven by tourism and associated supplemental feeding affects the energetics, corticosterone and immunity of six discrete populations of the northern Bahamian rock iguana (Cyclura cychlura inornata and Cyclura cychlura figginsi). We studied populations within and outside the reproductive season and quantified tourist numbers during sample collection. Specifically, we measured clutch size, body condition, plasma energy metabolites, reactive oxygen species, baseline corticosterone concentrations and immune function of male and female iguanas from each population to address whether (i) disparate physiologies are emerging across a gradient of tourism and feeding, (ii) both subspecies respond similarly and (iii) responses vary with season/reproductive condition. We found significant effects of tourism level, season and their interaction on the physiology of both C. c. inornata and C. c. figginsi, supporting the idea that tourism is leading to the divergence of phenotypes. Specifically, we found elevated plasma energy metabolites, oxidative stress and a measure of innate immunity (bactericidal ability), but reduced corticosterone concentrations with increasing tourism in both subspecies of rock iguanas. These physiological metrics differ according to the level of tourism in both subspecies and persist across seasons despite variation with natural seasonal and reproductive changes. These findings suggest that anthropogenic disturbance results in disparate physiologies in northern Bahamian rock iguanas.

Glucose tolerance of iguanas is affected by high-sugar diets in the lab and supplemental feeding by ecotourists in the wild.

There is great interspecific variation in the nutritional composition of natural diets, and the varied nutritional content is physiologically tolerated because of evolutionarily based balances between diet composition and processing ability. However, as a result of landscape change and human exposure, unnatural diets are becoming widespread among wildlife without the necessary time for evolutionary matching between the diet and its processing. We tested how a controlled, unnatural high glucose diet affects glucose tolerance using captive green iguanas, and we performed similar glucose tolerance tests on wild Northern Bahamian rock iguanas that are either frequently fed grapes by tourists or experience no such supplementation. We evaluated both short and longer-term blood glucose responses and corticosterone (CORT) concentrations as changes have been associated with altered diets. Experimental glucose supplementation in the laboratory and tourist feeding in the wild both significantly affected glucose metabolism. When iguanas received a glucose-rich diet, we found greater acute increases in blood glucose following a glucose challenge. Relative to unfed iguanas, tourist-fed iguanas had significantly lower baseline CORT, higher baseline blood glucose, and slower returns to baseline glucose levels following a glucose challenge. Therefore, unnatural consumption of high amounts of glucose alters glucose metabolism in laboratory iguanas with short-term glucose treatment and free-living iguanas exposed to long-term feeding by tourists. Based on these results and the increasing prevalence of anthropogenically altered wildlife diets, the consequences of dietary changes on glucose metabolism should be further investigated across species, as such changes in glucose metabolism have health consequences in humans (e.g. diabetes).

Climate effects on nesting phenology in Nebraska turtles.

A frequent response of organisms to climate change is altering the timing of reproduction, and advancement of reproductive timing has been a common reaction to warming temperatures in temperate regions. We tested whether this pattern applied to two common North American turtle species over the past three decades in Nebraska, USA. The timing of nesting (either first date or average date) of the Common Snapping Turtle (Chelydra serpentina) was negatively correlated with mean December maximum temperatures of the preceding year and mean May minimum and maximum temperatures in the nesting year and positively correlated with precipitation in July of the previous year. Increased temperatures during the late winter and spring likely permit earlier emergence from hibernation, increased metabolic rates and feeding opportunities, and accelerated vitellogenesis, ovulation, and egg shelling, all of which could drive earlier nesting. However, for the Painted Turtle (Chrysemys picta), the timing of nesting was positively correlated with mean minimum temperatures in September, October, December of the previous year, February of the nesting year, and April precipitation. These results suggest warmer fall, and winter temperature may impose an increased metabolic cost to painted turtles that impedes fall vitellogenesis, and April rains may slow the completion of vitellogenesis through decreased basking opportunities. For both species, nest deposition was highly correlated with body size, and larger females nested earlier in the season. Although average annual ambient temperatures have increased over the last four decades of our overall fieldwork at our study site, spring temperatures have not yet increased, and hence, nesting phenology has not advanced at our site for Chelydra. While Chrysemys exhibited a weak trend toward later nesting, this response was likely due to increased recruitment of smaller females into the population due to nest protection and predator control (Procyon lotor) in the early 2000s. Should climate change result in an increase in spring temperatures, nesting phenology would presumably respond accordingly, conditional on body size variation within these populations.

Turtles and Tortoises Are in Trouble.

Turtles and tortoises (chelonians) have been integral components of global ecosystems for about 220 million years and have played important roles in human culture for at least 400,000 years. The chelonian shell is a remarkable evolutionary adaptation, facilitating success in terrestrial, freshwater and marine ecosystems. Today, more than half of the 360 living species and 482 total taxa (species and subspecies combined) are threatened with extinction. This places chelonians among the groups with the highest extinction risk of any sizeable vertebrate group. Turtle populations are declining rapidly due to habitat loss, consumption by humans for food and traditional medicines and collection for the international pet trade. Many taxa could become extinct in this century. Here, we examine survival threats to turtles and tortoises and discuss the interventions that will be needed to prevent widespread extinction in this group in coming decades.

Understanding reproductive allometry in turtles: A slippery "slope".

Measures of reproductive output in turtles are generally positively correlated with female body size. However, a full understanding of reproductive allometry in turtles requires logarithmic transformation of reproductive and body size variables prior to regression analyses. This allows for slope comparisons with expected linear or cubic relationships for linear to linear and linear to volumetric variables, respectively. We compiled scaling data using this approach from published and unpublished turtle studies (46 populations of 25 species from eight families) to quantify patterns among taxa. Our results suggest that for log-log comparisons of clutch size, egg width, egg mass, clutch mass, and pelvic aperture width to shell length, all scale hypoallometrically despite theoretical predictions of isometry. Clutch size generally scaled at ~1.7 to 2.0 (compared to an isometric expectation of 3.0), egg width at ~0.5 (compared to an expectation of 1.0), egg mass at ~1.1 to 1.3 (3.0), clutch mass at ~2.5 to 2.8 (3.0), and pelvic aperture width at 0.8-0.9 (1.0). We also found preliminary evidence that scaling may differ across years and clutches even in the same population, as well as across populations of the same species. Future investigators should aspire to collect data on all these reproductive parameters and to report log-log allometric analyses to test our preliminary conclusions regarding reproductive allometry in turtles.

Geographic variation in thermal sensitivity of early life traits in a widespread reptile.

Taxa with large geographic distributions generally encompass diverse macroclimatic conditions, potentially requiring local adaptation and/or phenotypic plasticity to match their phenotypes to differing environments. These eco-evolutionary processes are of particular interest in organisms with traits that are directly affected by temperature, such as embryonic development in oviparous ectotherms. Here we examine the spatial distribution of fitness-related early life phenotypes across the range of a widespread vertebrate, the painted turtle (Chrysemys picta). We quantified embryonic and hatchling traits from seven locations (in Idaho, Minnesota, Oregon, Illinois, Nebraska, Kansas, and New Mexico) after incubating eggs under constant conditions across a series of environmentally relevant temperatures. Thermal reaction norms for incubation duration and hatchling mass varied among locations under this common-garden experiment, indicating genetic differentiation or pre-ovulatory maternal effects. However, latitude, a commonly used proxy for geographic variation, was not a strong predictor of these geographic differences. Our findings suggest that this macroclimatic proxy may be an unreliable surrogate for microclimatic conditions experienced locally in nests. Instead, complex interactions between abiotic and biotic factors likely drive among-population phenotypic variation in this system. Understanding spatial variation in key life-history traits provides an important perspective on adaptation to contemporary and future climatic conditions.

Marine subsidies likely cause gigantism of iguanas in the Bahamas.

We utilized natural experiment opportunities presented by differential conditions (presence/absence of seabirds and invasive species) on cays in the Bahamas to study whether interisland variations in food resources contributed to gigantism in Allen Cays Rock Iguanas (Cyclura cychlura inornata). We analyzed the stable carbon (δ13C) and nitrogen (δ15N) isotope values from iguana tissues and resources from each island food web to test the predictions that (1) food webs on islands with seabirds exhibit the influence of marine subsidies from seabird guano, whereas those from non-seabird islands do not, and (2) size differences in iguanas among cays were due to either (a) supplemental food availability from mice and/or seabird carcasses killed by barn owls (Tyto alba) and/or (b) access to more nutrient-rich vegetation fertilized by seabird guano. Food web components from the seabird island (Allen Cay) had 5-9‰ higher δ15N values than those on the other cays and Allen Cay plants contained nearly two times more nitrogen. Bayesian stable isotope mixing models indicated that C3 plants dominated iguana diets on all islands and showed no evidence for consumption of mice or shearwaters. The iguanas on Allen Cay were ~ 2 times longer (48.3 ± 11.6 cm) and ~ 6 times heavier (5499 ± 2847 g) than iguanas on other cays and this was likely from marine-derived subsidies from seabird guano which caused an increase in nitrogen concentration in the plants and a resultant increase in the δ15N values across the entire food web relative to non-seabird islands.

Extreme thermal fluctuations from climate change unexpectedly accelerate demographic collapse of vertebrates with temperature-dependent sex determination.

Global climate is warming rapidly, threatening vertebrates with temperature-dependent sex determination (TSD) by disrupting sex ratios and other traits. Less understood are the effects of increased thermal fluctuations predicted to accompany climate change. Greater fluctuations could accelerate feminization of species that produce females under warmer conditions (further endangering TSD animals), or counter it (reducing extinction risk). Here we use novel experiments exposing eggs of Painted Turtles (Chrysemys picta) to replicated profiles recorded in field nests plus mathematically-modified profiles of similar shape but wider oscillations, and develop a new mathematical model for analysis. We show that broadening fluctuations around naturally male-producing (cooler) profiles feminizes developing embryos, whereas embryos from warmer profiles remain female or die. This occurs presumably because wider oscillations around cooler profiles expose embryos to very low temperatures that inhibit development, and to feminizing temperatures where most embryogenesis accrues. Likewise, embryos incubated under broader fluctuations around warmer profiles experience mostly feminizing temperatures, some dangerously high (which increase mortality), and fewer colder values that are insufficient to induce male development. Therefore, as thermal fluctuations escalate with global warming, the feminization of TSD turtle populations could accelerate, facilitating extinction by demographic collapse. Aggressive global CO2 mitigation scenarios (RCP2.6) could prevent these risks, while intermediate actions (RCP4.5 and RCP6.0 scenarios) yield moderate feminization, highlighting the peril that insufficient reductions of greenhouse gas emissions pose for TSD taxa. If our findings are generalizable, TSD squamates, tuatara, and crocodilians that produce males at warmer temperatures could suffer accelerated masculinization, underscoring the broad taxonomic threats of climate change.

Climatic niche evolution in turtles is characterized by phylogenetic conservatism for both aquatic and terrestrial species.

Understanding how the climatic niche of species evolved has been a topic of high interest in current theoretical and applied macroecological studies. However, little is known regarding how species traits might influence climatic niche evolution. Here, we evaluated patterns of climatic niche evolution in turtles (tortoises and freshwater turtles) and whether species habitat (terrestrial or aquatic) influences these patterns. We used phylogenetic, climatic and distribution data for 261 species to estimate their climatic niches. Then, we compared whether niche overlap between sister species was higher than between random species pairs and evaluated whether niche optima and rates varied between aquatic and terrestrial species. Sister species had higher values of niche overlap than random species pairs, suggesting phylogenetic climatic niche conservatism in turtles. The climatic niche evolution of the group followed an Ornstein-Uhlenbeck model with different optimum values for aquatic and terrestrial species, but we did not find consistent evidence of differences in their rates of climatic niche evolution. We conclude that phylogenetic climatic niche conservatism occurs among turtle species. Furthermore, terrestrial and aquatic species occupy different climatic niches but these seem to have evolved at similar evolutionary rates, reinforcing the importance of habitat in understanding species climatic niches and their evolution.

Energetic investment associated with vitellogenesis induces an oxidative cost of reproduction.

Oxidative stress is a potential cost of reproduction, but conclusive evidence for this relationship is lacking. The goal of this study was to serially assess across a seasonal gradient the relationship between reproduction, circulating plasma energy metabolites and oxidative state. Here, we examine a study animal ideally suited to test for the oxidative costs of reproduction: the Allen Cays Rock Iguana. Female rock iguanas reproduce at varying frequencies, often skipping years, allowing for a comparison between reproductive and non-reproductive females during the same narrow, annual breeding season. This feature of iguana life history enabled us to address not just sex and seasonal differences in physiology, but also potential oxidative costs of reproduction in females. Male and female iguanas were sampled during the early (vitellogenic), late (gravid) and post-reproductive seasons. Ultrasound examinations were performed on females to quantify reproductive investment, and blood samples were collected for physiology assays, which included reactive oxygen metabolites (d-ROMs), antioxidants, triglycerides, free glycerol and glucose. The early reproductive season was characterized by significant increases in reproductive female's triglycerides, free glycerol and oxidative stress compared to late and post-reproductive periods and non-reproductive females and males during all sampling periods. Antioxidants were significantly elevated during the early reproductive season for reproductive females, non-reproductive females and males when compared to late and post-season. Follicle number in early reproductive females was positively related to d-ROMs, triglycerides and free glycerol, negatively related to antioxidants and showed no relationship with glucose. Measures of oxidative stress, d-ROMs and oxidative index were positively correlated with circulating levels of the lipid metabolite free glycerol during the early reproductive period, but this relationship weakened in the late season and disappeared in the post-season. Broadly, this study supports the hypothesis that the relationship between reproduction and oxidative stress is driven by energy investment, being greatest during early reproduction when vitellogenesis is occurring.

Altered spring phenology of North American freshwater turtles and the importance of representative populations.

Globally, populations of diverse taxa have altered phenology in response to climate change. However, most research has focused on a single population of a given taxon, which may be unrepresentative for comparative analyses, and few long-term studies of phenology in ectothermic amniotes have been published. We test for climate-altered phenology using long-term studies (10-36 years) of nesting behavior in 14 populations representing six genera of freshwater turtles (Chelydra, Chrysemys, Kinosternon, Malaclemys, Sternotherus, and Trachemys). Nesting season initiation occurs earlier in more recent years, with 11 of the populations advancing phenology. The onset of nesting for nearly all populations correlated well with temperatures during the month preceding nesting. Still, certain populations of some species have not advanced phenology as might be expected from global patterns of climate change. This collection of findings suggests a proximate link between local climate and reproduction that is potentially caused by variation in spring emergence from hibernation, ability to process food, and thermoregulatory opportunities prior to nesting. However, even though all species had populations with at least some evidence of phenological advancement, geographic variation in phenology within and among turtle species underscores the critical importance of representative data for accurate comprehensive assessments of the biotic impacts of climate change.

An evaluation of the use of pentosidine as a biomarker for ageing turtles.

Concentrations of the biomarker pentosidine have been shown to be useful measures of age for a number of avian and mammalian species. However, no study has examined its usefulness as an age marker in a long-lived ectotherm despite the fact that such a marker could prove useful in understanding age distributions of populations subject to conservation programmes. Therefore, we evaluated pentosidine concentrations in the interdigital webbing of 117 female yellow mud turtles (Kinosternon flavescens) at a 35 year study site in western Nebraska where nearly all turtles are of known age. Pentosidine concentrations were extraordinarily low and positively correlated with age in this turtle, but concentrations were too variable to permit precise estimates of age for turtles of unknown age. These results may reflect the remarkable physiological adaptations of this turtle to low temperatures and oxygen deprivation in a highly seasonal environment requiring prolonged hibernation. Whether pentosidine concentrations in other ectotherms occupying less seasonal environments would be more highly correlated with age remains to be determined. However, our results suggest that patterns of accumulation of pentosidine in ectotherms may be fundamentally different from those in endotherms.