Mature Erythrocytes of Iguana iguana (Squamata, Iguanidae) Possess Functional Mitochondria

Abstract

Electron microscopy analyses of Iguana iguana blood preparations revealed the presence of mitochondria within erythrocytes with well-structured cristae. Fluorescence microscopy analyses upon incubation with phalloidin-FITC, Hoechst 33342 and mitochondrial transmembrane potential (Δψ_m)-sensitive probe MitoTracker Red indicated that mitochondria i) widely occur in erythrocytes, ii) are polarized, and iii) seem to be preferentially confined at a "perinuclear" region, as confirmed by electron microscopy. The analysis of NADH-dependent oxygen consumption showed that red blood cells retain the capability to consume oxygen, thereby providing compelling evidence that mitochondria of Iguana erythrocytes are functional and capable to perform oxidative phosphorylation.

About the Authors:

Giuseppina Di Giacomo

Affiliation: Dipartimento di Biologia, Università Tor Vergata, Roma, Italia

Silvia Campello

Affiliations Dipartimento di Biologia, Università Tor Vergata, Roma, Italia, IRCCS Fondazione Santa Lucia, Roma, Italia

Mauro Corrado

Affiliations IRCCS Fondazione Santa Lucia, Roma, Italia, Istituto Telethon Dulbecco, Istituto Veneto di Medicina Molecolare, Padova, Italia

Livia Di Giambattista

Affiliation: Dipartimento di Biologia, Università Tor Vergata, Roma, Italia

Claudia Cirotti

Affiliation: Dipartimento di Biologia, Università Tor Vergata, Roma, Italia

Giuseppe Filomeni

* E-mail: [email protected] (GF); [email protected] (GG)

Affiliations Dipartimento di Biologia, Università Tor Vergata, Roma, Italia, Cell Stress and Survival Unit, Danish Cancer Society Research Center, Copenhagen, Denmark

Gabriele Gentile

* E-mail: [email protected] (GF); [email protected] (GG)

Affiliation: Dipartimento di Biologia, Università Tor Vergata, Roma, Italia

Introduction

In vertebrates, erythrocytes originate from multipotent hematopoietic stem cells that start a process of differentiation leading to first burst-forming units-erythroid (BFUe), then colony-forming units-erythroid (CFUe), normoblasts, erythroblasts, reticulocytes, and ultimately mature erythrocytes [1]. In mammals, during erythroblast maturation, erythroblasts normally lose their nucleus and mitochondria. In fact, in mammals, nuclear retention in erythroblasts is associated with specific conditions or pathologies [2]. In a traditional view, mammals’ anucleate-erythrocyte condition would be considered as an evolutionary response to the need to optimize surface/volume for hemoglobin and the need to cope with a circulatory system that presents end-blood-vessels smaller than in other vertebrates [3]. Flexibility and augmented capability of big nucleated erythrocytes to squeeze through small capillaries would be derived from the extrusion of the nucleus and also other cellular organelles, such as endoplasmic reticulum, not required for erythrocytes’ main function as oxygen carriers. However, such a theory has recently been questioned and evidence has been provided showing that regulation of cytoskeleton for deformation into the biconcave shape may not necessarily require nuclear extrusion. Instead, the extrusion of both mitochondria and nucleus may help mammal erythrocytes against oxidative stress caused by high-sugar and high-heme conditions [4]. Moreover, the extrusion of both mitochondria and nucleus could in general be beneficial for erythrocytes as both source and target of damage derived from the generation of reactive oxygen species (ROS) would be eliminated. As mitochondria consume oxygen during oxidative phosphorylation, their extrusion during erythrocyte maturation should also increase the efficiency of oxygen transport in mature erythrocytes. Thus, in this perspective, loss of mitochondria could benefit erythrocytes of all vertebrates. Such a hypothesis has been recently tackled and contrasted by Stier and collaborators [5] who demonstrated the presence of functional mitochondria in zebra finch erythrocytes and found that levels of oxidative stress in the blood of zebra finches are not higher than those observed in a species of mammal of comparable size.

While it is known that nucleus is present in erythrocytes of all vertebrates, with a few exceptions including mammals, it is less clear to what extent mature nucleated erythrocytes in vertebrates also possess functional mitochondria. Importantly, Stier and collaborators [5] applied a multiple approach to investigate the presence of functional mitochondria in mature nucleated erythrocytes in birds, as prior studies, based on biochemistry or ultra-morphology, had led to incomplete results [6–10]. Direct and indirect evidence has been cumulated for fish, showing that functional mitochondria occur in erythrocytes of several species [11–15]. In amphibians, mitochondria were found close to the nuclear envelope in mature erythrocytes of the urodele Triturus cristatus [16] and few anuran (Rana) species [7,17]. For the last, indirect evidence of functional mitochondria was also provided [6,18]. For reptiles, Yasuzumi [19] stated that mitochondria were found close to the nuclear envelope in mature erythrocytes of turtles, although no image was provided. Indirect evidence was instead offered to support the existence of functional mitochondria in erythrocytes of a viperid snake [20] and an agamid lizard [21].

Here we used electron and fluorescence microscopy to document the existence of mitochondria in the cytoplasm of mature erythrocytes of Iguana iguana (Squamata, Iguanidae). We also used a Clark-type oxygen electrode to demonstrate electron flow to oxygen as a result of oxidative phosphorylation in Iguana erythrocytes.

Materials and Methods

Blood collection and fractionation

Blood (2 mL) was drawn from the caudal vein of 2 adult green iguana (Iguana iguana) males and 2 adult females hosted at the BIOPARCO of Rome (Italy). Iguanas were not related by kinship. Blood samples were obtained using heparinized syringes and stored in ice after collection. Samples were then immediately transported to the laboratories of the Department of Biology (Tor Vergata University) for subsequent analyses. As soon as in the laboratory, whole blood was fractionated by centrifugation at 1500 g for 12 min at room temperature with gently acceleration and deceleration to avoid mixing different phases. This methodology allows separate the blood into an upper plasma layer, a lower red blood cell layer, and a thin interface containing the white blood cells (buffy coat). Plasma and, a later stage, also leucocytes were aspirated off from the top of the tube, while erythrocytes were collected by the bottom to avoid any even little contamination. Erythrocytes thus prepared were next used for different analyses.

Ethic statement

Animal manipulation and blood sampling were performed according to a protocol that minimized animal stress, in accordance with the European Community guidelines and with the approval of relevant National (Ministry of Health) and local (Institutional Animal Care and Use, Tor Vergata University; BIOPARCO Foundation) ethical committees.

Erythrocyte mitochondrial DNA amplification

In order to demonstrate the presence of mitochondrial DNA in Iguana erythrocytes, we extracted DNA from purified erythrocytes of each of the four iguanas using the QIAamp DNA Mini Kit (Qiagen). We PCR amplified and sequenced a mitochondrial DNA target 903-bp region including 709 bases of the 3’ end of the ND4 subunit of the nicotinamide adenine dinucleotide dehydrogenase gene and the tRNA genes histidine, serine, and leucine (in part). Primers and PCR conditions were as used by Malone and collaborators [22]. Sequences were deposited in Genbank (accessions KT326936, KT326937, KT326938, and KT326939). To safely exclude amplification of nuclear copies of mtDNA (numts) we compared our sequences with the corresponding functional sequence (accession AJ278511) obtained from DNA extracted from heart and liver [23]. Sequences were aligned by using MUSCLE ver. 3.8.31 [24]. We used DNASP [25] to calculate the number of synonymous, non-synonymous, and stop-codon substitutions between sequences.

Electron microscopy evaluation of erythrocyte ultrastructure

Erythrocyte ultrastructure was obtained by transmission electron microscopy (TEM). Red blood cells were separated and fixed with 2.5% glutaraldehyde with 0.1 M sodium phosphate, pH 7.4 for 1 hour at ice temperature. Sample were post-fixed with osmium tetroxide, stained with uranyl acetate, dehydrated in ethanol and embedded in Epon resin. After sectioning, samples then were collected on uncoated nickel grids and observed.

Fluorescence microscopy analyses

The mitochondrial transmembrane potential (Δψ_m) of iguana's erythrocytes was analyzed taking advantage of the Δψ_m-sensitive probe MitoTracker Red CMXRos by fluorescence microscopy. In detail, Iguana blood samples were diluted in PBS (1:5 v/v) and stained for 30 min with 50 nM of MitoTracker Red CMXRos. This is a derivative of X-Rosamine that enters the cells due its lipophilic nature and accumulates into functional mitochondria driven by Δψ_m. Once in the organelles, the probe is very well retained also upon fixation with paraformaldehyde and, at variance with other probes (e.g., Mito Tracker Red FM, or Green FM), is recommended by the manufacturer for applications in fixed cells, and for multicolor labeling experiments. Indeed, Mito Tracker Red CMXRos red fluorescence is well resolved from the green fluorescence of other probes. Afterwards, blood was smeared on a d-polylisinated slide, fixed with 4% paraformaldehyde for 10 minutes and permeabilized with 0.2% Triton X-100. Each step, blood smear was gently washed three times with cold PBS to prevent any loss of the cells. Next, FITC-conjugated phalloidin (1:1000) was added to the blood smears and incubated at RT for 40 min to visualize actin fibers, that especially in the erythrocytes, are mostly confined underneath the cell membrane. Hoechst 33342, that we used to stain DNA (and therefore identify the cell region occupied by nuclei), enters equally living and permeabilized cells and strongly bounds dsDNA. Therefore, after washing with PBS, blood smears were incubated with 1mL of Hoechst 33342 (1:10000 in PBS) for 30 min, gently washed three times in PBS and visualized by fluorescence microscopy.

Oxygen consumption measurement

Oxygen consumption in intact erythrocytes was determined at 25°C using a Clark-type oxygen electrode equipped with thermostatic control and magnetic stirring. Purified erythrocytes obtained from 0.25 mL of Iguana blood specimens were re-suspended in the same volume of PBS and added to in 1.25 mL of measurement buffer (125 mM KCl, 10 mM Tris/MOPS, pH 7.4, 10 μM Tris/EGTA, pH 7.4, 1 mM K₂HPO₄) upon incubation with digitonin. 5 mM glutamate, 2.5 mM malate and 1 mM NADH were also added to the experimental buffer to sustain complex I-dependent oxygen consumption. In particular, NADH provides the reducing equivalents for the electron transport chain, while malate and glutamate fuel the malate-aspartate shuttle, which is required for the virtual intra-matrix uptake of NADH. At the end of each measurement, 5 μM of complex I (NADH dehydrogenase) inhibitor rotenone was always added to verify that oxygen consumption was effectively driven by complex I. Measurements were also repeated in the presence of catalase at the concentration of 0.5 U in order to rule out any side oxygen consumption due to its partial reduction to H₂O₂. Indeed, it is well known that H₂O₂ production frequently occurs in erythrocytes starting from one-electron reduction of oxygen by heme-contained iron of hemoglobin. Superoxide produced (spontaneously or by superoxide dismutase-mediated catalysis) dismutates in H₂O₂. Therefore, the addition of catalase into the oxygraphic chamber is a tool commonly used to understand whether oxygen is fully reduced by electron transport chain to H₂O, or only partially to H₂O₂. In the first case catalase addition is basically ineffective and it does not change oxygen consumption rate. By contrast, when hydrogen peroxide is produced by any source, catalase rapidly dismutates H₂O₂ into water and molecular oxygen through the reaction

In accordance with this reaction, for every 2 moles of H₂O₂, catalase re-generates 1 mole of O₂, this increasing oxygen concentration inside the oxygraphic chamber and virtually decreasing oxygen consumption up to half, case in which all the available oxygen is reduced to H₂O₂ and not to H₂O by the mitochondrial respiratory chain.

Statistical analyses

We used the Wilcoxon-Mann-Whitney test of a difference between two groups to test the null hypotheses that oxygen consumption was not different between: i) solutions containing only buffer (no treatment) and solutions containing buffer and erythrocytes (treatment 1); ii) solutions containing only buffer (no treatment) and solutions containing buffer and erythrocytes + catalase (treatment 2). We also used the Wilcoxon signed-rank test for matched pairs to test the null hypothesis that oxygen consumption was not different between treatments 1 and 2. As statistical power is affected by the size of the effect and the size of the sample used to detect it, we performed a post-hoc analysis of achieved power (1-β) for each test, given α = 0.05, N_sample1 = 4, N_sample2 = 4, and effect size (d and d_z), with correlation (r) between treatments 1 and 2 being equal to 0.85. In a conservative approach, we estimated power of tests using the minimum asymptotic relative efficiency method for two-tailed tests, as implemented in G*Power ver.3.1.9.2 [26].

Remaining statistical analysis was performed by using Past ver. 2.16 [27].

Results and Discussion

The four mitochondrial DNA sequences obtained were between 95% and 96% similar to the sequence used as functional reference. No stop codon substitutions were observed in none of the four coding sequences analyzed, while the maximum number of synonymous and non-synonymous substitutions between sequences was 26 and 4, respectively. The maximum number of synonymous and non-synonymous substitutions between sequences increased to 30 and 5, respectively, when the reference sequence was included in the analysis. This is strongly suggestive that, for all the four iguanas, we amplified a functional copy of the ND4 gene from erythrocyte mtDNA.

In order to assess whether erythrocytes from Iguana iguana contain mitochondria, we first performed electron microscopy analyses of blood preparations to visualize ultrastructure details of red blood cells. Fig 1 shows representative images of nucleated erythrocytes displaying organelles that did not appear widely distributed in the cytosol, but were present in small numbers, mainly localized at the nuclear border. High magnification images unambiguously showed that those organelles were mitochondria with well-structured cristae. However, due to the cut plane, some mitochondria were not included in the section, and only some erythrocytes showed further intra-cytoplasmic structures resembling an undefined pool of internal membranes (e.g., traces of endoplasmic reticulum). Based on these results, which were consistent with prior studies of vertebrate erythrocytes reporting occasional retention of mitochondria and other organelles during the maturation phases [16,28], we could not assess whether mitochondria were regularly or just occasionally present in Iguana erythrocytes. In fact, we were still not able to distinguish if they were in a dismantling phase, reasonably during erythrocytes maturation or conversely, whether they were retaining their functionality.

[Figure omitted. See PDF.]

Fig 1. Ultrastructure of Iguana erythrocytes.

(A-D) Four different fields of blood specimens containing Iguana erythrocytes. Nuclei, mitochondria (mit) and uncharacterized internal membrane pools resembling endoplasmic reticulum (*) are indicated. High magnification images (4 and 10X) are shown on the right of each picture.

https://doi.org/10.1371/journal.pone.0136770.g001

To address these issues, we performed fluorescence microscopy analyses upon incubation with phalloidin-FITC, Hoechst 33342 and the mitochondrial transmembrane potential (Δψ_m)-sensitive probe MitoTracker red CMXRos to highlight cell boundaries, nuclei and mitochondria, respectively. Images shown in Fig 2 unambiguously indicate that mitochondria widely occur in all erythrocytes collected from Iguana and they seem to be preferentially confined at a "perinuclear" region, reinforcing the evidence obtained by electron microscopy.

[Figure omitted. See PDF.]

Fig 2. Fluorescence images of Iguana erythrocytes.

Iguana erythrocytes stained with: phalloidin-FITC (green) to visualize actin fibers below the plasma membrane; the DNA-binding dye Hoechst 33342 (blue) to highlight nuclei; the Δψ_m-sensitive probe MitoTracker Red CMXRos (red) to localize polarized mitochondria. Superimposition (merge) of the three fluorescences is shown at the bottom, along with a 3X magnification detail of one erythrocyte.

https://doi.org/10.1371/journal.pone.0136770.g002

The fluorescence microscopy analyses provided the direct evidence that each iguanas’ erythrocyte has mitochondria, clearly suggesting that, by retaining their Δψ_m, they are virtually able to generate the proton-motive force required to produce adenosine triphosphate (ATP).

To finally answer this question, we measured NADH-dependent oxygen consumption in iguanas’ erythrocytes by a Clark-type oxygraph. In order to rule out that oxygen consumption was due to any possible side partial reduction to H₂O₂ catalyzed by heme iron of hemoglobin (or by mitochondria itself by a futile partial reduction of oxygen), and not driven by mitochondrial electron transport chain, we performed these analyses in the presence or absence of the H₂O₂-scavenging enzyme catalase. Indeed, by catalyzing the disproportionation of H₂O₂ in O₂ and H₂O, catalase regenerates oxygen in the oxygraphic chamber and virtually produces a decrease of the oxygen consumption rate, which will be proportional to the concentration of H₂O₂. Raw data are presented in Table 1.

[Figure omitted. See PDF.]

Table 1. Oxygen consumption of Iguana erythrocytes.

https://doi.org/10.1371/journal.pone.0136770.t001

Values displayed in Fig 3 represent the mean ± SD of single measurements performed on N = 4 individuals. The achieved power was > 0.98 for tests i) and ii), confirming that N_sample1 = 4 and N_sample2 = 4 were adequate for our statistical analyses. For the Wilcoxon signed-rank test between treatments 1 and 2, the achieved power was 0.45.

[Figure omitted. See PDF.]

Fig 3. Oxygen consumption of Iguana erythrocytes.

Purified erythrocytes obtained from 0.25 mL of blood were dissolved in 1.25 mL of measurement buffer containing 2.5 mM malate, 5 mM glutamate and 1 mM NADH to fuel complex I. Upon incubation with digitonin, oxygen consumption was measured for up to 5 minutes in the presence or absence of catalase (0.5 U) and expressed as nmol of oxygen per minute. Five μM rotenone was finally added to the oxygraphic chamber in order to irreversibly inhibit complex I activity and verify that oxygen consumption was dependent on complex I-driven mitochondrial respiration. (Buffer, measurement buffer without red blood cells; RBC, red blood cells; RBC + rot, red blood cells incubated with rotenone; RBC + cat, red blood cells incubated with catalase). *p = 0.030; n.s., not significant (p = 0.176).

https://doi.org/10.1371/journal.pone.0136770.g003

Results show that: i) red blood cells retained the capability to consume oxygen, and ii) catalase did not significantly modulate oxygen concentration, suggesting that this was not (or only moderately, in case we made a Type II error) associated with a side-production of H₂O₂, but reasonably due to mitochondrial activity.

Energy is required by erythrocytes to sustain several physiological functions and, with regard to the energetic metabolism—although with differences among taxa—nucleated mature erythrocytes have been reported to be metabolically more active than mammalian anucleated ones [29]. Besides glycolysis, functional oxidative phosphorylation could guarantee further energy to drive several energy-requiring processes and provide metabolites capable to modulate respiratory physiology [9]. In fact, among the potential roles of mitochondria inside nucleated erythrocytes, recently discussed by Stier and collaborators [5], it has been shown that the rate of ATP production by oxidative phosphorylation modulates hemoglobin affinity to oxygen in snake erythrocytes [20]. In particular, it has been reported that, under hypoxic conditions, the concentration of nucleoside triphosphates (primarily ATP) decreases with a trend inversely correlated with oxygen affinity [30]. Bartlett [31] reviewed the occurrence of possible hemoglobin-affinity modifiers − i.e. red-cell organic phosphates (RCOP) − in different vertebrate species including Iguana iguana, with emphasis on 2,3-diphosphoglycerate (DPG), inositol pentaphosphate (IP₅), ATP and guanosine triphosphate (GTP). Interestingly, out of the four compounds analyzed, ATP resulted highly predominant in Iguana erythrocytes. Additionally, it has been shown that Iguana hemoglobin possesses structural requisites for ATP binding [32]. Admittedly, other factors also contribute to affect oxygen-affinity of hemoglobin in reptiles with adaptation to hypoxia, including temperature, other allosteric interactions with cellular effectors, and gene-based changes in hemoglobin structure and intrinsic oxygen-binding [33–35]. However, the role of oxidative phosphorylation, as long as the cellular mechanisms responsible for decreases in ATP concentration during hypoxia still await an in-depth investigation in reptile nucleated erythrocytes.

Acknowledgments

We gratefully thank the Foundation BIOPARCO di Roma and in particular Fulvio Fraticelli, Klaus G. Friedrich, Stefano Micarelli, Daniele Macale, Giuliana Galetto, Diego Reggianti, and Massimiliano Pellicciotta for their assistance and participation. We thank an anonymous reviewer for the constructive criticism expressed on a previous version of this work.

Author Contributions

Conceived and designed the experiments: GF GG. Performed the experiments: GF GDG SC CC MC LDG. Analyzed the data: GF GDG SC CC GG. Contributed reagents/materials/analysis tools: GF GG. Wrote the paper: GG GF. Blood drawing and preparation for the experiments: GG. Obtained permission to draw and use blood: GG. Contributed to the critical review of the manuscript: GDG SC CC.

Citation: Di Giacomo G, Campello S, Corrado M, Di Giambattista L, Cirotti C, Filomeni G, et al. (2015) Mature Erythrocytes of Iguana iguana (Squamata, Iguanidae) Possess Functional Mitochondria. PLoS ONE 10(9): e0136770. https://doi.org/10.1371/journal.pone.0136770

References

1. Moritz KM, Lim GB, Wintour EM (1997) Developmental regulation of erythropoietin and erythropoiesis. Am J Physiol Regul Integr Comp Physiol 273: R1829–R1844.

2. Thrall MA. Erythrocyte Morphology. In: Thrall MA, Weiser G, Allison R, Campbell TW, editors. Veterinary hematology and clinical chemistry. Oxford: Wiley-Blackwell; 2012. pp. 61–74.

3. Snyder GK, Sheafor BA (1999) Red blood cells: centerpiece in the evolution of the vertebrate circulatory system. Am Zool 39: 189–198.

4. Zhang ZW, Cheng J, Xu F, Chen YE, Du JB, Yuan M, et al. (2011) Red blood cell extrudes nucleus and mitochondria against oxidative stress. Iubmb Life 63: 560–565. pmid:21698761

5. Stier A, Bize P, Schull Q, Zoll J, Singh F, Geny B, et al. (2013) Avian erythrocytes have functional mitochondria, opening novel perspectives for birds as animal models in the study of ageing. Front Zool 10: 33. pmid:23758841

6. Hunter AS, Hunter FR (1957) A comparative study of erythrocyte metabolism. J Cell Comp Physiol 49: 479–502.

7. Davies HG (1961) Structure in nucleated erythrocytes. J Biophys Biochem Cytol 9: 671–687. pmid:13720098

8. Harris JR, Brown JN (1971) Fractionation of the avian erythrocyte: an ultrastructural study. J Ultrastruct Res 36: 8–23. pmid:4105480

9. Beam KG, Alper SL, Palade GE, Greengard P (1979) Hormonally regulated phosphoprotein of turkey erythrocytes: localization to plasma membrane. J Cell Biol 83: 1–15. pmid:229109

10. Brasch K, Adams GH, Neelin JM (1974) Evidence for erythrocyte-specific histone modification and structural changes in chromatin during goose erythrocyte maturation. J Cell Sci 15: 659–677. pmid:4426929

11. Boutilier RG, Ferguson RA (1989) Nucleated red-cell function: metabolism and pH regulation. Can J Zool 67: 2986–2993.

12. Eddy FB (1977) Oxygen uptake by rainbow trout blood, Salmo gairdneri. J Fish Biol 10: 87–90.

13. Phillips MCL, Moyes CD, Tufts BL (2000) The effects of cell ageing on metabolism in rainbow trout (Oncorhynchus mykiss) red blood cells. J Exp Biol 203: 1039–1045. pmid:10683163

14. Pica A, Scacco S, Papa F, De Nitto E, Papa S (2001) Morphological and biochemical characterization of mitochondria in Torpedo red blood cells. Comp Biochem Physiol B Biochem Mol Biol 128: 213–219. pmid:11207435

15. Tiihonen K, Nikinmaa M (1991) Substrate utilization by carp (Cyprinus Carpio) erythrocytes. J Exp Biol 161: 509–514.

16. Tooze J, Davies HG (1967) Light- and electron- microscope studies on the spleen of the newt Triturus cristatus: the fine structure of erythropoietic cells. J Cell Sci 2: 617–640. pmid:6080549

17. Lessler MA, Herrera FM (1962) Electron-microscope studies of x-ray damage to frog blood cells. Radiat Res 17: 111–117. pmid:14464354

18. Bratosin D, Estaquier J, Slomianny C, Tissier J-P, Quatannens B, Bulai T, et al. (2004) On the evolution of erythrocyte programmed cell death: apoptosis of Rana esculenta nucleated red blood cells involves cysteine proteinase activation and mitochondrion permeabilization. Biochimie 86: 183–192. pmid:15134833

19. Yasuzumi G (1960) Licht- und elektronenmikroskopische Studien an kernhaltingen Erythrocyten. Z Zellforsch 51: 325–335.

20. Ogo SH, Bernardes CF, Glass ML, Torsoni MA, Vercesi AE (1993) Functional mitochondria in snake Bothrops alternatus erythrocytes and modulation of HbO2 affinity by mitochondrial ATP. J Comp Physiol B 163: 614–619.

21. Olsson M, Wilson M, Uller T, Mott B, Isaksson C, Healey M, et al. (2008) Free radicals run in lizard families. Biol Lett 4: 186–188. pmid:18211861

22. Malone CL, Wheeler T, Taylor JF, Davis SK (2000) Phylogeography of the Caribbean Rock Iguana (Cyclura): implications for conservation and insights on the biogeographic history of the West Indies. Mol Phylogen Evol 17: 269–279.

23. Janke A, Erpenbeck D, Nilsson M, Arnason U (2001) The mitochondrial genomes of the iguana (Iguana iguana) and the caiman (Caiman crocodylus): implications for amniote phylogeny. Proc R Soc Lond B Biol Sci 268: 623–631.

24. Edgar RC (2004) MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucl Acids Res 32: 1792–1797. pmid:15034147

25. Librado P, Rozas J (2009) DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics. 251451–1452.

26. Faul F, Erdfelder E, Lang AG, Buchner A (2007) G*Power 3: A flexible statistical power analysis program for the social, behavioral, and biomedical sciences. Behav Res Methods 39: 175–191. pmid:17695343

27. Hammer Ø, Harper DAT, Ryan PD (2001) PAST: Paleontological Statistics Software Package for Education and Data Analysis. Palaeontol Electron 4: 1–9.

28. Zentgraf H, Deumling B, Jarasch ED, Franke WW (1971) Nuclear membranes and plasma membranes from hen erythrocytes. I. Isolation, characterization, and comparison. J Biol Chem 246: 2986–2995. pmid:4251996

29. Nikinmaa M. Vertebrate red blood cells: adaptations of function to respiratory requirements. Berlin: Springer-Verlag; 1990.

30. Herman JK, Ingermann RL (1996) Effects of hypoxia and hyperoxia on oxygen-transfer properties of the blood of a viviparous snake. J Exp Biol 199: 2061–2070. pmid:9319971

31. Bartlett GR (1980) Phosphate compounds in vertebrate red blood cells. Am Zool 20: 103–114.

32. Rucknagel KP, Braunitzer G, Wiesner H (1988) Hemoglobins of reptiles. The primary structures of the alpha I- and beta I-chains of common iguana (Iguana iguana) hemoglobin. Biol Chem Hoppe Seyler 369: 1143–1150. pmid:3242545

33. Grigg GC, Wells RMG, Beard LA (1993) Allosteric control of oxygen binding by haemoglobin during embryonic development in the crocodile Crocodylus porosus: the role of red cell organic phosphates and carbon dioxide. J Exp Biol 175: 15–32.

34. Weber RE (2007) High-altitude adaptations in vertebrate hemoglobins. Respir Physiol Neurobiol 158: 132–142. pmid:17561448

35. Damsgaard C, Storz JF, Hoffmann FG, Fago A (2013) Hemoglobin isoform differentiation and allosteric regulation of oxygen binding in the turtle, Trachemys scripta. Am J Physiol Regul Integr Comp Physiol 305: R961–R967. pmid:23986362