Developmental plasticity of cardiac anoxia-tolerance in juvenile common snapping turtles ( Chelydra serpentina)

Living on the edge: Crayfish as drivers to anoxification of their own shelter microenvironment

Burrowing is a common trait among crayfish thought to help species deal with adverse environmental challenges. However, little is known about the microhabitat ecology of crayfish taxa in relation to their burrows. To fill this knowledge gap, we assessed the availability of oxygen inside the crayfish shelter by series of in-vivo and in-silico modelling experiments. Under modeled condition, we found that, except for the entrance region of the 200 mm, a flooded burrow microenvironment became anoxic within 8 h, on average. Multiple 12-hour day-night cycles, with burrows occupied by crayfish for 12 h and empty for 12 h, were not sufficient for refreshing the burrow microenvironment. We then examined the degree to which crayfish species with different propensities for burrowing are tolerant of self-created anoxia. From these experiments, primary and secondary burrowers showed best and most consistent tolerance-exhibiting ≥ 64% survival to anoxia and 25-91% survival of ≥ 9 h at anoxia, respectively. Tertiary burrowers exhibited little to no tolerance of anoxia with 0-50% survival to anoxia and only one species exhibiting survival (2%) of ≥ 9 h at anoxia. Results suggest that moderate to strongly burrowing crayfish can quickly draw down the dissolved oxygen in burrow water but appear to have conserved a legacy of strong tolerance of anoxia from their monophyletic ancestors-the lobsters-whereas tertiary burrowers have lost (or never evolved) this ability.

Hypoxic incubation at 50% of atmospheric levels shifts the cardiovascular response to acute hypoxia in American alligators, Alligator mississippiensis

We designed a series of studies to investigate whether hypoxia (10% O2) from 20% of incubation to hatching, or from 20 to 50% of incubation, affects cardiovascular function when juvenile American alligators reached an age of 4-5 years compared to juveniles that were incubated in 21% O2. At this age, we measured blood flows in all the major arteries as well as heart rate, blood pressure, and blood gases in animals in normoxia and acute hypoxia (10% O2 and 5% O2). In all three groups, exposure to acute hypoxia of 10% O2 caused a decrease in blood O2 concentration and an increase in heart rate in 4-5-year-old animals, with limited effects on blood flow in the major outflow vessels of the heart. In response to more acute hypoxia (5% O2), where blood O2 concentration decreased even further, we measured increased heart rate and blood flow in the right aorta, subclavian artery, carotid artery, and pulmonary artery; however, blood flow in the left aorta either decreased or did not change. Embryonic exposure to hypoxia increased the threshold for eliciting an increase in heart rate indicative of a decrease in sensitivity. Alligators that had been incubated in hypoxia also had higher arterial PCO2 values in normoxia, suggesting a reduction in ventilation relative to metabolism.

Effects of Developmental Hypoxia on the Vertebrate Cardiovascular System

Developmental hypoxia has profound and persistent effects on the vertebrate cardiovascular system, but the nature, magnitude, and long-term outcome of the hypoxic consequences are species specific. Here we aim to identify common and novel cardiovascular responses among vertebrates that encounter developmental hypoxia, and we discuss the possible medical and ecological implications.

Diverse Response Pattern to Anoxia in Three Freshwater Turtle Species

With increasing water eutrophication and global warming, anoxia and hypoxia are becoming more and more common in water environments. Most vertebrates have a limited tolerance to anoxia of only a few minutes, but some species, such as turtles, can survive for months being exposed to anoxia. Antioxidant defense systems may have a potential role in resisting anoxia stress in freshwater turtles. The three-keeled pond turtle Chinemys reevesii, the snapping turtle Chelydra serpentina and the soft-shelled turtle Pelodiscus sinensis are three popular aquaculture species and share similar habitats in China. While C. reevesii and C. serpentina are hard-shelled turtles with poor skin permeability, P. sinensis is soft-shelled turtle whose skin permeability is good. We examined the antioxidant defense responses in different tissues of the three turtle species under acute anoxia stress for 10 h and subsequently recovered for 24 h in order to reveal the response patterns of the antioxidant defense system of the three turtle species that differed in morphological structure and life history strategy. We found that the antioxidant response patterns to acute anoxia stress were tissue- and species-specific. The soft-shelled turtle was more sensitive to anoxia than the hard-shelled turtles. Under anoxia stress, the three species kept the activities of most antioxidant enzymes stable. C. reevesii and P. sinensis were highly dependent on vitamin C in antioxidant defense, while high activities of structural antioxidant enzymes were found in the tissues of C. serpentina. The above diverse patterns may be related with adaptive evolution of morphological structure and physiological functions of the three turtle species.

No evidence for pericardial restraint in the snapping turtle (Chelydra serpentina) following pharmacologically-induced bradycardia at rest or during exercise

Most animals elevate cardiac output during exercise through a rise in heart rate (f_H), whilst stroke volume (V_S) remains relatively unchanged. Cardiac pacing reveals that elevating f_H alone does not alter cardiac output, which is instead largely regulated by the peripheral vasculature. In terms of myocardial oxygen demand, an increase in f_H is more costly than that which would incur if V_S instead were to increase. We hypothesized that f_H must increase because any substantial rise in V_S would be constrained by the pericardium. To investigate this hypothesis, we explored the effects of pharmacologically-induced bradycardia, with ivabradine treatment, on V_S at rest and during exercise in the common snapping turtle (Chelydra serpentina) with intact or opened pericardium. We first showed that, in isolated myocardial preparations, ivabradine exerted a pronounced positive inotropic effect on atrial tissue, but only minor effects on ventricle. Ivabradine reduced f_H in vivo, such that exercise tachycardia was attenuated. Pulmonary and systemic V_S rose in response to ivabradine. The rise in pulmonary V_S largely compensated for the bradycardia at rest, leaving total pulmonary flow unchanged by ivabradine, although ivabradine reduced pulmonary blood flow during swimming (exercise x ivabradine interaction, P<0.05). Although systemic V_S increased, systemic blood flow was reduced by ivabradine both at rest and during exercise, in spite of ivabradine's potential to increase cardiac contractility. Opening the pericardium had no effect on f_H, V_S or blood flows before or after ivabradine, indicating that the pericardium does not constrain V_S in turtles, even during pharmacologically-induced bradycardia.

Changes in hemoglobin function and isoform expression during embryonic development in the American alligator, Alligator mississippiensis

In the developing embryos of egg-laying vertebrates, O₂ flux takes place across a fixed surface area of the eggshell and the chorioallantoic membrane. In the case of crocodilians, the developing embryo may experience a decrease in O₂ flux when the nest becomes hypoxic, which may cause compensatory adjustments in blood O₂ transport. However, whether the switch from embryonic to adult hemoglobin isoforms (isoHbs) plays some role in these adjustments is unknown. Here, we provide a detailed characterization of the developmental switch of isoHb synthesis in the American alligator, Alligator mississippiensis. We examined the in vitro functional properties and subunit composition of purified alligator isoHbs expressed during embryonic developmental stages in normoxia and hypoxia (10% O₂). We found distinct patterns of isoHb expression in alligator embryos at different stages of development, but these patterns were not affected by hypoxia. Specifically, alligator embryos expressed two main isoHbs: HbI, prevalent at early developmental stages, with a high O₂ affinity and high ATP sensitivity, and HbII, prevalent at later stages and identical to the adult protein, with a low O₂ affinity and high CO₂ sensitivity. These results indicate that whole blood O₂ affinity is mainly regulated by ATP in the early embryo and by CO₂ and bicarbonate from the late embryo until adult life, but the developmental regulation of isoHb expression is not affected by hypoxia exposure.

Developmental programming of DNA methylation and gene expression patterns is associated with extreme cardiovascular tolerance to anoxia in the common snapping turtle

Background

Environmental fluctuation during embryonic and fetal development can permanently alter an organism’s morphology, physiology, and behaviour. This phenomenon, known as developmental plasticity, is particularly relevant to reptiles that develop in subterranean nests with variable oxygen tensions. Previous work has shown hypoxia permanently alters the cardiovascular system of snapping turtles and may improve cardiac anoxia tolerance later in life. The mechanisms driving this process are unknown but may involve epigenetic regulation of gene expression via DNA methylation. To test this hypothesis, we assessed in situ cardiac performance during 2 h of acute anoxia in juvenile turtles previously exposed to normoxia (21% oxygen) or hypoxia (10% oxygen) during embryogenesis. Next, we analysed DNA methylation and gene expression patterns in turtles from the same cohorts using whole genome bisulfite sequencing, which represents the first high-resolution investigation of DNA methylation patterns in any reptilian species.

Results

Genome-wide correlations between CpG and CpG island methylation and gene expression patterns in the snapping turtle were consistent with patterns observed in mammals. As hypothesized, developmental hypoxia increased juvenile turtle cardiac anoxia tolerance and programmed DNA methylation and gene expression patterns. Programmed differences in expression of genes such as SCN5A may account for differences in heart rate, while genes such as TNNT2 and TPM3 may underlie differences in calcium sensitivity and contractility of cardiomyocytes and cardiac inotropy. Finally, we identified putative transcription factor-binding sites in promoters and in differentially methylated CpG islands that suggest a model linking programming of DNA methylation during embryogenesis to differential gene expression and cardiovascular physiology later in life. Binding sites for hypoxia inducible factors (HIF1A, ARNT, and EPAS1) and key transcription factors activated by MAPK and BMP signaling (RREB1 and SMAD4) are implicated.

Conclusions

Our data strongly suggests that DNA methylation plays a conserved role in the regulation of gene expression in reptiles. We also show that embryonic hypoxia programs DNA methylation and gene expression patterns and that these changes are associated with enhanced cardiac anoxia tolerance later in life. Programming of cardiac anoxia tolerance has major ecological implications for snapping turtles, because these animals regularly exploit anoxic environments throughout their lifespan.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13072-021-00414-7.

Indirect evidence that anoxia exposure and cold acclimation alter transarcolemmal Ca2+ flux in the cardiac pacemaker, right atrium and ventricle of the red-eared slider turtle (Trachemys scripta)

We indirectly assessed if altered transarcolemmal Ca²⁺ flux accompanies the decreased cardiac activity displayed by Trachemys scripta with anoxia exposure and cold acclimation. Turtles were first acclimated to 21 °C or 5 °C and held under normoxic (21N; 5N) or anoxic conditions (21A; 5A). We then compared the response of intrinsic heart rate (f_H) and maximal developed force of spontaneously contracting right atria (F_max,RA), and maximal developed force of isometrically-contracting ventricular strips (F_max,V), to Ni²⁺ (0.1-10 mM), which respectively blocks T-type Ca²⁺ channels, L-type Ca²⁺ channels and the Na⁺-Ca²⁺-exchanger at the low, intermediate and high concentrations employed. Dose-response curves were established in simulated in vivo normoxic (Sim Norm) or simulated in vivo anoxic extracellular conditions (Sim Anx; 21A and 5A preparations). Ni²⁺ decreased intrinsic f_H, F_max,RA and F_max,V of 21N tissues in a concentration-dependent manner, but the responses were blunted in 21A tissues in Sim Norm. Similarly, dose-response curves for F_max,RA and F_max,V of 5N tissues were right-shifted, whereas anoxia exposure at 5 °C did not further alter the responses. The influence of Sim Anx was acclimation temperature-, cardiac chamber- and contractile parameter-dependent. Combined, the findings suggest that: (1) reduced transarcolemmal Ca²⁺ flux in the cardiac pacemaker is a potential mechanism underlying the slowed intrinsic f_H of anoxic turtles at 21 °C, but not 5 °C, (2) a downregulation of transarcolemmal Ca²⁺ flux may aid cardiac anoxia survival at 21 °C and prime the turtle myocardium for winter anoxia and (3) confirm that altered extracellular conditions with anoxia exposure can modify turtle cardiac transarcolemmal Ca²⁺ flux.

Prenatal hypoxia affects scaling of blood pressure and arterial wall mechanics in the common snapping turtle, Chelydra serpentina

In reptiles, exposure to hypoxia during embryonic development affects several cardiovascular parameters. These modifications may impose different mechanical stress to the arterial system, and we speculated that the arterial wall of major outflow vessels would be modified accordingly. Since non-crocodilian reptiles possess a partially divided ventricle, ensuing similar systemic and pulmonary systolic pressures, we investigated how morphological and mechanical properties of segments from the left aortic arch (LAo) and the proximal and distal segments of the left pulmonary artery (LPAp and LPAd, respectively) change as body mass (M_b) increases. Eggs from common snapping turtles, Chelydra serpentina, were incubated under normoxia (21% O₂; N21) or hypoxia (10% O₂; H10), hatched and maintained in normoxia thereafter. Turtles (0.11-6.85 kg) were cannulated to measure arterial pressures, and an injection of adrenaline was used to increase pressures. Portions of the LAo, LPAp and LPAd were fixed under physiological hydrostatic pressures for histology and mechanical assessment. Arterial pressures increased with M_b for N21 but not for H10. Although mechanical and functional characteristics from the LPAp and LPAd were similar between N21 and H10, wall thickness from LAo did not change with M_b in the H10 group, thus wall stress increased in larger turtles. This indicates that larger H10 turtles probably experience an elevated probability of arterial wall rupture without concomitant changes in the cardiovascular system to prevent it. Finally, collagen content of the LPAp and LAo was smaller than in LPAd, suggesting a more distensible arterial wall could attenuate higher pressures from larger turtles.

The Long-Term Effects of Developmental Hypoxia on Cardiac Mitochondrial Function in Snapping Turtles

It is well established that adult vertebrates acclimatizing to hypoxic environments undergo mitochondrial remodeling to enhance oxygen delivery, maintain ATP, and limit oxidative stress. However, many vertebrates also encounter oxygen deprivation during embryonic development. The effects of developmental hypoxia on mitochondrial function are likely to be more profound, because environmental stress during early life can permanently alter cellular physiology and morphology. To this end, we investigated the long-term effects of developmental hypoxia on mitochondrial function in a species that regularly encounters hypoxia during development-the common snapping turtle (Chelydra serpentina). Turtle eggs were incubated in 21% or 10% oxygen from 20% of embryonic development until hatching, and both cohorts were subsequently reared in 21% oxygen for 8 months. Ventricular mitochondria were isolated, and mitochondrial respiration and reactive oxygen species (ROS) production were measured with a microrespirometer. Compared to normoxic controls, juvenile turtles from hypoxic incubations had lower Leak respiration, higher P:O ratios, and reduced rates of ROS production. Interestingly, these same attributes occur in adult vertebrates that acclimatize to hypoxia. We speculate that these adjustments might improve mitochondrial hypoxia tolerance, which would be beneficial for turtles during breath-hold diving and overwintering in anoxic environments.

Draft Genome of the Common Snapping Turtle, Chelydra serpentina, a Model for Phenotypic Plasticity in Reptiles

Turtles are iconic reptiles that inhabit a range of ecosystems from oceans to deserts and climates from the tropics to northern temperate regions. Yet, we have little understanding of the genetic adaptations that allow turtles to survive and reproduce in such diverse environments. Common snapping turtles, Chelydra serpentina, are an ideal model species for studying adaptation to climate because they are widely distributed from tropical to northern temperate zones in North America. They are also easy to maintain and breed in captivity and produce large clutch sizes, which makes them amenable to quantitative genetic and molecular genetic studies of traits like temperature-dependent sex determination. We therefore established a captive breeding colony and sequenced DNA from one female using both short and long reads. After trimming and filtering, we had 209.51Gb of Illumina reads, 25.72Gb of PacBio reads, and 21.72 Gb of Nanopore reads. The assembled genome was 2.258 Gb in size and had 13,224 scaffolds with an N50 of 5.59Mb. The longest scaffold was 27.24Mb. BUSCO analysis revealed 97.4% of core vertebrate genes in the genome. We identified 3.27 million SNPs in the reference turtle, which indicates a relatively high level of individual heterozygosity. We assembled the transcriptome using RNA-Seq data and used gene prediction software to produce 22,812 models of protein coding genes. The quality and contiguity of the snapping turtle genome is similar to or better than most published reptile genomes. The genome and genetic variants identified here provide a foundation for future studies of adaptation to climate.