Alterations to cavefish red blood cells provide evidence of adaptation to reduced subterranean oxygen

Scientific Reports volume 12, Article number: 3735 (2022) Cite this article

3212 Accesses

3 Citations

108 Altmetric

Metrics details

Animals inhabiting extreme environments allow the powerful opportunity to examine adaptive evolution in response to diverse pressures. One such pressure is reduced oxygen, commonly present at high-altitude and subterranean environments. Cave-dwelling animals must also deal with darkness and starvation, both of which have been rigorously studied as key forces driving the evolution of cave-associated traits. Interestingly, hypoxia as an environmental pressure has received less attention. Here we examined putatively adaptive phenotypes evolving in a freshwater teleost fish, Astyanax mexicanus, which includes both surface- and cave-dwelling forms. This model system also provides the opportunity to identify convergent responses to hypoxia, owing to the presence of numerous natural and independently-colonised cave populations, alongside closely-related surface conspecifics. The focus of this study is hemoglobin, an essential molecule for oxygen transport and delivery. We found that multiple cave populations harbor a higher concentration of hemoglobin in their blood, which is coincident with an increase in cave morph erythrocyte size compared to surface fish. Interestingly, both cave and surface morphs have comparable numbers of erythrocytes per unit of blood, suggesting elevated hemoglobin is not due to overproduction of red blood cells. Alternatively, owing to an increased cell area of erythrocytes in cavefish, we reason that they contain more hemoglobin per erythrocyte. These findings support the notion that cavefish have adapted to hypoxia in caves through modulation of both hemoglobin production and erythrocyte size. This work reveals an additional adaptive feature of Astyanax cavefish, and demonstrates that coordinated changes between cellular architecture and molecular changes are necessary for organisms evolving under intense environmental pressure.

Our understanding of the complex interactions between extreme environmental conditions and organismal response to these pressures remains incomplete. A powerful natural model to study adaptation in extreme environments is the blind Mexican cavefish, Astyanax mexicanus (Fig. 1a). This species allows for direct comparisons of two extant morphotypes, a river-dwelling ‘surface’ morph and obligate subterranean ‘cave’ morphs. Thirty extant cave populations1 are distributed across networks of subterranean caves within the exposed limestone karst of the Sierra de El Abra region of northeast Mexico (Fig. 1a). Phenotypic changes have accompanied this surface-to-cave transition, over an estimated period of ~ 20 to 500 thousand years2,3, including loss of eyes, pigmentation and enhancement of non-visual senses. In the 8 decades since their discovery4, considerable progress has been made towards understanding the mechanisms underlying regressive features, however the environmental pressures of the El Abra karst system have mostly examined the effects of limited light and nutrition5.

Multiple populations of Astyanax mexicanus cavefish have higher concentrations of hemoglobin compared to surface fish. Astyanax mexicanus consists of cave and surface morphotypes. The cave morphs are distributed across numerous caves of the Sierra de El Abra in northeastern Mexico, while surface morphs reside in nearby rivers, streams and lakes (a). Each of three cavefish populations displayed statistically significantly higher concentrations of hemoglobin than the conspecific surface morph (b, Chica WRS p = 0.005, Tinaja WRS p = 0.0005, Pachón WRS p = 0.0001). Hemoglobin concentrations were measured via spectrophotometric analysis (c). No outliers were detected.

One environmental pressure commonly present in subterranean caves, reduced oxygen levels, has received far less attention. Subterranean environments across the globe frequently have reduced oxygen levels compared to surrounding terrestrial environments due to organismal respiration6, absence of photosynthetic input6, and limited air mixing with the surface atmosphere7. Empirical measurements of dissolved oxygen have been conducted in two El Abra caves, Pachón and Tinaja, which demonstrate significantly lower oxygen levels in cave pools compared to surrounding surface waters (DO = 2.97 mg/L in the Pachón cave compared to 8.20 mg/L in the surface environment (Rascón); 59% saturation in the Tinaja cave compared to 80% saturation in the surface environment (Nacimiento del Río Choy)8,9. Although this environmental feature has most likely impacted the evolution of Astyanax cavefish, few studies have addressed the putative adaptive features mitigating this pressure10,11,12.

Owing to the global effects of climate change13,14,15, and the consequential impact on natural dissolved oxygen levels, hypoxia tolerance has been studied in several teleost species. Complex responses to hypoxia include behavioral changes (decreased predator avoidance in grey mullets16), morphological changes (gill remodeling in crucian carp17) and molecular changes (increased expression of the oxygen sensing gene HIF-1 in zebrafish18).

For sufficient oxygen to reach peripheral tissues, vertebrates utilize hemoglobin proteins as oxygen carriers. Hemoglobin molecules include four tetramers, each of which binds a single oxygen molecule, loading four oxygen molecules per hemoglobin protein. Natural oxygen concentrations can vary considerably depending on the environment, and many animals require precise oxygen levels to maintain their standard metabolic rates. Therefore, hemoglobin genes and proteins have likely been subject to evolutionary changes across vertebrate taxa19. These include numerical variation in the genomic organization of hemoglobin genes20, increasing hemoglobin concentration in tambaqui21, or increased hemoglobin binding affinity in red drum22. Here, we examined hemoglobin concentration since it is the primary molecular transport for oxygen delivery23.

We identified changes in cavefish impacting hemoglobin concentration and red blood cell size, which are likely adaptive in the low oxygen cave environment. These alterations, which are absent from closely-related surface morphs, likely confer hypoxia tolerance in Astyanax cavefish given that they were reared (and persist) under the same captive conditions for generations. Thus, this model provides the unique opportunity to identify the genetic basis for changes to hemoglobin transport and hypoxia tolerance in the wild.

Hemoglobin is the primary oxygen transport molecule in nearly all vertebrates19. We first reasoned that hemoglobin levels may be higher in cave morphs as an adaptation to reduced oxygen in the cave environment (compared to overlying terrestrial environments6) based on direct measurements from the El Abra caves where Astyanax cavefish are found8,9. One adaptive mechanism for low oxygen is increasing hemoglobin concentration24,25 to augment the available oxygen for systemic transport.

Using a spectrophotometric approach ("Material and methods"), we found the mean hemoglobin concentration in surface fish (n = 26) was 8.25 ± 0.23 g/dL (Fig. 1b). All three cavefish populations examined (Chica, Tinaja and Pachón cavefish) had significantly higher mean hemoglobin concentrations than surface fish (Chica n = 8, mean = 9.86 ± 0.43 g/dL, Wilcoxon Rank-Sum (WRS) p = 0.005; Tinaja n = 8, mean = 10.60 ± 0.38 g/dL, WRS p = 0.0005; Pachón n = 30, mean = 10.49 ± 0.39 g/dL, WRS p = 0.0001; Fig. 1b), but did not differ from one another (Fig. 1b). Presumably, higher hemoglobin concentrations allow cave morphs to cope with limited dissolved oxygen in the cave. However, diverse cellular features could conceivably contribute to increased hemoglobin concentration, so we next examined cellular features of erythrocytes (i.e., red blood cells).

Hematocrit is the portion of a given blood volume occupied by erythrocytes (Fig. 2a). Certain clinical phenotypes in humans, e.g. polycythemia vera26, are diagnosed based on elevated hematocrit values27,28. Therefore, higher hemoglobin concentrations in cavefish may be due to presence of more erythrocytes compared to surface fish. To determine if elevated hemoglobin is mediated by overproduction of red blood cells, we compared hematocrit between surface and cave morphs.

Hematocrit trends higher in cavefish than surface fish. The fraction of red blood cells to total blood volume (hematocrit) was determined via centrifugation of microcapillary tubes. The erythrocytic phase collects at the distal end of the capillary tube while the buffy coat and plasma collect more proximally (a). Although statistical significance was not achieved when comparing populations individually, (compared to surface; Chica WRS p = 0.115, Tinaja WRS p = 0.018, Pachón WRS p = 0.161), a clear trend of higher hematocrit was noted for each cavefish population compared to surface (b). When comparing by morphotype (surface against all cavefish), however, cavefish do have a statistically significant higher mean hematocrit (surface n = 8, mean = 28.51 ± 0.03%, Cave n = 24, mean = 36.14 ± 0.01%, WRS p = 0.024) than surface fish. One outlier (the lowest Chica value) was detected although statistical significance was unchanged (Table S6).

Hematocrit measures trended higher in every cave population we examined compared to surface fish (Fig. 2b). We first performed a comparison (Wilcoxon Rank Sum test) between surface morphs and cavefish morphs (irrespective of individual cavefish localities) which revealed significant differences between these groups (surface n = 8, mean = 28.51 ± 0.03%, cave n = 24, mean = 36.14 ± 0.01%, WRS p = 0.024). We then performed a post-hoc analysis comparing each of the four populations to one another. These results clearly demonstrated the same trend, but did not achieve statistical significance following Bonferroni correction (compared to surface; Chica WRS p = 0.115, Tinaja WRS p = 0.018, Pachón WRS p = 0.161, owing to the Bonferroni correction, p values must be < 0.0083 to achieve statistical significance). The uniformly elevated hematocrit in cavefish, however, was not surprising given that cavefish have higher hemoglobin concentrations. Consistent with hemoglobin measures, the Tinaja cave population harbored the highest hematocrit value (mean hematocrit = 39.23 ± 0.03%, Fig. 2b), followed by the Pachón cave (mean hematocrit = 35.56 ± 0.03%, Fig. 2b), with the Chica cave population demonstrating the lowest hematocrit (mean hematocrit = 33.63 ± 0.01%, Fig. 2b) of the cave populations. Given these results we next examined numerical variation in red blood cells for each population.

We anticipated cavefish populations may have a higher density of erythrocytes than surface fish reflecting higher hemoglobin concentrations and hematocrit values. Unexpectedly, surface fish had the highest density of erythrocytes (n = 8, mean density = 2.1 × 106 ± 1.4 × 105/mm3, Fig. 3c-c′) followed closely by Pachón cavefish (n = 8, mean density = 2.0 × 106 ± 8.2 × 105/mm3, Fig. 3f-f′). Interestingly, Tinaja and Chica cavefish had slightly lower mean densities (n = 8, 1.6 × 106 ± 1.8 × 105/mm3 and n = 8, 1.5 × 106 ± 6.7 × 104/mm3, Fig. 3d–e′ respectively). Although Chica was statistically different from Pachón (WRS p = 0.001, Fig. 3b), these individuals were obtained commercially as adults ("Material and methods") and therefore their precise age is unknown. Overall, erythrocyte densities were comparable across surface and cave morphs, which was not anticipated given the higher hemoglobin concentrations and hematocrit values in the cave populations. Since erythrocyte density was generally reduced in cavefish, we considered whether differences in erythrocyte cell size could explain the higher hemoglobin concentration in cave morphs.

Erythrocyte density is comparable across Astyanax morphotypes. We determined the density of erythrocytes in a standard blood volume using a hemocytometer (a). Although statistical difference was detected for the Chica population (p < 0.0023), mean erythrocyte density was comparable across morphotypes (b). Erythrocytes were stained with Acridine Orange, visualized under light microscopy (c–f) and counted using ImageJ29 (v2.0.0) (c′–f′, scale bar 130 µm). Two outliers were detected (the lowest value in surface and highest in Chica). Removal of these outliers resulted in a statistically significant difference between Chica and Surface (Table S6).

We measured two-dimensional erythrocyte area in standard blood volumes to determine if red blood cell sizes vary across cave and surface populations. Blood smears were collected, stained and analyzed using light microscopy. Erythrocytes were unambiguously identified ("Material and methods"), manually scored in pixels units (ImageJ29 v2.0.0) and converted to µm for direct comparison. Since cavefish have elevated hematocrit and hemoglobin but comparable erythrocyte densities, we anticipated each cave population may show larger erythrocyte areas. Indeed, surface fish erythrocytes were the smallest across all populations (n = 8, mean area = 71.36 ± 1.43 µm2, Fig. 4a), while erythrocytes from the Tinaja population were on average the largest of the cave populations (n = 8, mean area = 86.05 ± 0.96 µm2, Fig. 4c), followed by Chica (mean area = 85.70 ± 1.38 µm2, Fig. 4b) and Pachón cavefish (mean area = 84.06 ± 1.52 µm2, Fig. 4d). Erythrocytes for all three cave populations were significantly larger than surface fish (Chica WRS p = 0.0002, Tinaja WRS p = 0.0002, Pachón WRS p = 0.0006, Fig. 4e). This cell size increase was not accompanied by an increase in nuclear area. For each population, including surface fish, the two-dimensional area of the nucleus did not differ by more than 0.5 µm2 (14.9 ± 0.24 µm2, Fig. 4f) suggesting erythrocyte area differences are attributable to increased cell cytoplasm, wherein hemoglobin is found.

Two-dimensional surface area of cave erythrocytes is significantly larger than surface erythrocytes. Erythrocytes from each population were stained with Wright–Giemsa and visualized under light microscopy (a–d, scale bar 10 µm). The total (a–d, yellow ring) and nuclear (a–d, dark blue ring) two-dimensional area was calculated for each erythrocyte in ImageJ29 (v2.0.0). The total surface area was significantly higher in all cave populations compared to surface (p < 0.00003) but did not differ from each other (e). Although the total two-dimensional area was higher in cavefish, there were no differences in the two-dimensional area of the nuclei (f), suggesting the differences in total two-dimensional area are due to differences in the cytoplasm, where hemoglobin is found. Two outliers were detected (the highest Chica value and the lowest Pachón value). Removal of these outliers did not change statistical significance (Table S6).

If differences in hematocrit between populations are not governed by increased numbers of erythrocytes, we reasoned that hematocrit variation may be correlated with differences in erythrocyte size. Qualitative comparisons seemed to demonstrate this trend, so to test this quantitatively, we performed correlations for both mean erythrocyte size and mean erythrocyte density, relative to mean hematocrit value. Erythrocyte density measures showed poor correlation to hematocrit (r = − 0.433, Fig. 5a). In contrast, erythrocyte area demonstrated a strong positive correlation with hematocrit (r = 0.869, Fig. 5b), suggesting that variation in elevated hematocrit across Astyanax cavefish populations is a function of erythrocyte size.

Erythrocyte size is positively correlated with hematocrit values in Astyanax cave and surface fish. We performed Pearson's correlations to determine which erythrocytic feature (density or size) is correlated with hematocrit. Mean erythrocyte density was not correlated with hematocrit (r = − 0.4329, a). However, population erythrocyte two-dimensional surface area provided a strong, positive correlation with hematocrit value (r = 0.8688, b).

We evaluated this further by calculating mean corpuscular hemoglobin (MCH, "Material and methods"), a metric that integrates hemoglobin concentration with erythrocyte density to estimate the mass of hemoglobin in a single red blood cell. Consistent with other findings, the surface population had the lowest mean MCH value (mean MCH = 39.41 pg, Table 1) while the Chica population had the highest of the three cave populations, followed by Tinaja, and Pachón (mean MCH = 66.01 pg, 65.02 pg, and 52.20 pg, respectively; Table 1). Thus, increased erythrocyte size (linked to elevated hematocrit) likely contributes to increased hemoglobin levels in natural cavefish populations.

In humans, hematocrit scores provide a relative measure of the percent erythrocyte contribution to a given volume of blood, owing to the fact that red blood cell size is stable across the human population30. Cavefish have higher levels of hemoglobin compared to surface fish. Hematocrit levels are elevated in cavefish, but this is likely a function of larger (rather than more) erythrocytes.

Here, we report that Astyanax cave morphs (from three distinct cave localities: Chica, Pachón and Tinaja) have convergently evolved higher concentrations of blood hemoglobin compared to closely-related surface fish. This elevated hemoglobin could enable cavefish to deliver more oxygen to systemic tissues per unit time. This may allow aerobic respiration to be maintained during periods of reduced oxygen to support sustained foraging activity31,32,33,34. Consistent with increased hemoglobin, cave morphs have higher hematocrit scores. Although this typically signals a greater density of erythrocytes27,28, cavefish did not have more red blood cells compared to surface fish. Rather, erythrocyte sizes varied across cave and surface populations, with cavefish uniformly harboring increased red blood cell sizes which explain, in part, their higher levels of hemoglobin.

Among the challenging environmental pressures encountered by cavefish, such as perpetual darkness and limited nutrition, include lower oxygen compared to overlying surface waters6. Limited oxygen has been shown to affect behavior, physiology, and certain morphologies in other teleost species16,17,18,19,20,21,22,23,24,25,35. Future studies in our system provide the opportunity to directly compare these features between the conspecific ‘surrogate ancestral’ morph—living in well-oxygenated environments—to derived cave morphs that have thrived in low-oxygen waters for tens to hundreds of thousands of years2,3.

The discrepancy in erythrocyte size across Astyanax populations provides at least one explanation for increased hematocrit and hemoglobin concentration in cave morphs. Interestingly, erythrocyte size varies greatly across teleost species. Snyder and Sheafor36 suggested that variation in erythrocyte size maintains efficient oxygen diffusion, but as a function of varying capillary diameters. Although both cell and capillary sizes can vary across taxa, the relationship between the two remain constant. On average, erythrocyte width is 25% larger than capillary diameter36. Interestingly, they found no relationship between erythrocyte size and oxygen carrying capacity (including hemoglobin concentration), thus concluding that erythrocyte size varies to maintain efficient flow and diffusion.

Lay and Baldwin37 investigated whether nucleus size or oxygen delivery better explains erythrocyte size variation in 52 teleost species, and found an inverse relationship between erythrocyte size and hemoglobin concentration. Interestingly, they also reported an inverse relationship between erythrocyte size and aerobic swimming ability28,38. Specifically, fish that swim with short bursts of rapid movement tend to have smaller erythrocytes, and increased aerobic swimming was related to higher numbers of erythrocytes39.

One explanation for this relationship may be that a larger surface area-to-volume ratio in erythrocytes reduces oxygen diffusion distance, allowing faster loading of oxygen when needed37. Since volume increases faster than surface area, and cavefish erythrocytes are larger than those of the surface fish, cavefish erythrocytes likely have a lower surface area-to-volume ratio. Cave morph activity has been characterized by a higher mean velocity compared to surface fish34, however surface fish utilize quick, short, ‘burst’ movements40. Since cavefish possess fewer, but larger, erythrocytes this could reflect their differing natural activity levels. Similar results were discussed by Tang et al.12 in the context of differences in heart size and shape between morphotypes. Cavefish possess smaller, less triangular but therefore ‘spongier’ hearts than surface fish. Spongier hearts have a smaller ratio of compact wall tissue to trabeculae which is associated with a sedentary activity profile41. Further, hearts with more trabeculae have a larger surface area exposed to circulating blood, allowing for increased gas exchange which therefore may offer an advantage for cavefish in hypoxic environments12.

Here we measured several blood phenotypes that likely provide an adaptive advantage for Astyanax cavefish in hypoxic caves. Much is still unknown, however, including the underlying mechanisms that enable hypoxia resistance in this species. Future studies designed to characterize overall gene expression changes in response to hypoxia captivity42,43,44, may reveal key differentially expressed genes impacting hypoxia tolerance. For instance, certain populations of cichlids native to Lake Victoria45,46 demonstrate differential expression of hemoglobin based on variable exposure to hypoxia by ancestral populations. Additionally, coding sequence changes in hemoglobin genes may impact function of the encoded protein, possibly improving oxygen binding capacity in variable ambient oxygen levels over evolutionary time20.

A considerable body of research has identified dozens of adaptive responses to hypoxia in teleosts47. Other adaptations in Astyanax may include changes to gill anatomy and function. For instance, Moran et al.11 noted that Pachón cavefish gills are heavier compared to surface morphs, and the mean gill mass of Pachón/surface F2 hybrids was greater than the parental surface population. This may indicate positive selection occurring as a consequence of hypoxia in caves11. Crucian carp can ‘remodel’ their gills to improve oxygen exchange and limit unnecessary, energetically costly ion exchange17,48. Conceivably, similar alterations to gill morphology may provide more efficient oxygen uptake and energy conservation in hypoxic caves.

Finally, in addition to the evolution of regressive features (e.g., eye loss), many constructive non-visual sensory traits have evolved in cavefish, such as expansions of the lateral line49,50 and taste systems51,52. One unexamined sensory modality in Astyanax is oxygen sensing. Neuroepithelial cells (NECs) are the primary oxygen sensing cells in fish and are distributed throughout a thin epithelial layer within gill filaments and lamellae, covering the efferent aspect of the gill53. Owing to the positioning of these cells, they may play a role in sensing both external (ambient water) and internal (arteriole) changes in the partial pressure of oxygen53. Further, these cells are homologous to type I chemoreceptors located in the carotid body of mammals54. Key differences in ion channel expression of NECs contribute to oxygen tolerance in anoxia-resistant goldfish55. Given that other sensory systems are augmented as a consequence of life in the extreme cave environment, similar alterations in oxygen sensing may have similarly evolved to cope with reduced oxygen conditions.

All animals used in this study were maintained in a satellite aquatic facility at the University of Cincinnati and housed in either 5- or 10-gallon tanks in a custom, reverse-osmosis husbandry unit with multiple filters (carbon, micron, UV and dense particulate filters; Aquaneering, San Diego, CA). This system conditions water to a temperature of 24 ° (± 2 °C), pH of 7.4 (± 0.2), and conductivity of 750 µS/cm (± 50 µS/cm) using real-time dose monitoring. Both Astyanax cave populations (Pachón, Tinaja and Chica) and surface morphs were housed under 12 h:12 h light:dark cycle and fed once daily (~ 9:00am) with a mixture of dry flake food (TetraMin Pro) and system water. Surface fish, Tinaja and Pachón cavefish used in this study were derived from breeding adult Astyanax mexicanus specimens originally provided by Dr. Richard Borowsky. Chica cavefish were commercially acquired, and therefore the precise age of these individuals is unknown. Certain blood parameters, including hematocrit, can change with age in humans, so we minimized biological variation in this study by excluding fish < 40 mm in standard length. Standard length measurements were not recorded for every individual, although every individual did exceed 40 mm standard length.

Animals were anesthetized using ice-cold system water and blood was collected via the caudal vein using a 31G syringe (BD Ultra-Fine™, BD Biosciences, San Jose, CA). We used fish of comparable age and standard length across all populations and found that neither factor significantly impacted any phenotypic score in this study (Supplemental Tables S1–S4). Collections occurred between 12:00 and 5:00 pm to minimize any potential diel effects. All housing conditions and collection methods were approved by the University of Cincinnati Institutional Animal Care and Use Committee (IACUC protocol # 10-01-21-01).

Hemoglobin concentrations were measured with an Aimstrip Hemoglobin Meter (Germaine Laboratories, San Antonio, TX, USA). The meter was calibrated and optically verified before each use per the manufacturer's instructions. 10 µL of whole blood was placed on a test strip (Fig. 1c) and inserted into the Hemoglobin Meter. Concentrations were automatically calculated via spectrophotometry and displayed in g/dL.

Hematocrit values were measured via centrifugation using an LW Scientific ZipCombo centrifuge (LW Scientific, Lawrenceville, GA) outfitted with a 12-place microhematocrit rotor. Whole blood was collected into a capillary tube and sealed at the opposite end with clay. Each capillary tube was spun for 3 min at 12,000 rpm, per manufacturer's recommendation. Blood components were separated into phases: the far more numerous and larger erythrocytes collect in the bottom phase, while the buffy coat (mostly white blood cells and platelets) and blood plasma layer on top of the erythrocytes (Fig. 2a). The length of the erythrocytic phase and the length of the total blood volume were measured. The hematocrit fraction was determined by dividing the length of the erythrocytic phase by the total blood volume.

To determine the density of erythrocytes for each population, an Improved Neubauer Hemocytometer (Hausser Scientific, Horsham, PA) was used (Fig. 3a). 1 µL whole blood was diluted into a staining solution of 500 µL Phosphate Buffered Saline (1×) and 20 µL of a working solution of the nuclear stain Acridine Orange (A1301, Invitrogen by Thermo Fisher Scientific, Waltham, MA, USA). This dilution was incubated at room temperature for 20 min. Following incubation, 10 µL were placed onto the hemocytometer and visualized using a Leica DM2000 LED compound microscope (Leica Microsystems, Wetzlar, Germany). Images were captured with a Leica DMC4500 camera mounted on the microscope and Leica Application Suite (LAS) X software (version 3.0.1.15878).

The hemocytometer is equipped with intricate grid etchings (Fig. 3a) to enable a precise method of counting cell number within a precise area. The largest visible square area using our compound scope was 0.75 mm × 0.75 mm or 0.5625 mm2. Three different 0.5625 mm2 areas were counted and averaged for each individual.

To estimate the total volume, the hemocytometer was designed so that shoulders on either side of the platform support a cover slip exactly 0.1 mm above the platform. Thus, the total volume in this section is 0.05625 mm3. Cells counts were performed using ImageJ29 (v2.0.0) by first using the ‘Color Threshold’ tool. The threshold was adjusted manually for each image to eliminate optical artifacts. The square area was then selected and counting was completed using the ‘Analyze Particles’ tool. This tools permits size exclusion to further ensure artifacts are not counted. Erythrocytes are much larger than white blood cells, therefore we excluded features less than 30 µm2. We calculated the number of cells per 1 mm3 (i.e., cell density) by dividing the raw cell count by the volume (0.05625 mm3), and multiplying by the dilution factor (520).

7.5 µL whole blood was applied to a microscope slide, smeared, air-dried and stained using Wright-Giemsa, Modified stain according to the manufacturers instruction (WG16, Millipore Sigma, Burlington, MA, USA). Nine cells were randomly selected and manually measured and averaged for each individual. The orientation of each selected cell was parallel to the slide (the short and long cell diameter at maximum length), not touching another cell, and free from obvious deformity. All cells were imaged and measured using the ‘Polygon Selection’ and ‘Measure’ tools in ImageJ29. Pixels were converted to µm for each image using the imaged scale bar (LAS X) during microscopy and the ‘Straight Line’ and ‘Set Scale’ tools in ImageJ29. The nuclear area of each erythrocyte was measured using the identical protocol employed to measure total two-dimensional area, as described above.

All statistical analyses were completed in R56 (version 3.6.1) using the base ‘stats’56 and ‘Rmisc’57 packages. At least eight individuals from each population were used for each assay as this was the greatest number of eligible individuals from the smallest in-lab population. Blinding was incorporated when possible by one investigator collecting data (imaging cells), creating an identification key, and removing individual nomenclature before a second investigator quantified and returned the data to investigator one. Wilcoxon rank-sum (WRS) tests were used to detect statistical significance for hemoglobin concentrations, hematocrit, erythrocyte density and erythrocyte size. All p values for these assays were Bonferroni corrected for multiple comparisons. Statistical significance was achieved if p < 0.0083. Pearson's correlations were calculated in R56 using the base ‘stats’56 package. Normality was tested for each population in each assay using the Shapiro–Wilk normality test. All populations followed normal distribution with the exception of Pachón 2d nuclear area measurements (Table S5). An additional analysis was conducted to test the presence of outliers. We tested all populations in each assay using the Hampel filter58. This test did detect outliers, however, removing the outliers affected only a single statistical significance (Fig. 3b, Table S6.). All collected data (including outliers) is presented in this manuscript. Both male and female fish were included in all assays. Sex does not appear to be a biological factor in the blood parameters we examined, based on Wilcoxon rank-sum test analyses (n = 42 for hemoglobin concentration and n = 32 for hematocrit, erythrocyte density, and 2d area measurements, each morph included) for each assay (Supplemental Tables S1–S4). Mean corpuscular hemoglobin (MCH) estimates the mass of hemoglobin protein per erythrocyte. This value is calculated by dividing the number of grams of hemoglobin per liter by the number of red blood cells in millions per milliliter providing picograms of hemoglobin per red blood cell59. Scatter plots were created in R56 using the ‘ggpubr’60 package. Boxplots were created in R56 using the ‘ggplot2’61 package. In these plots, whiskers extend to the smallest or largest value not surpassing 1.5xIQR and horizontal lines represent first and third quartiles and the median. All images were edited in Microsoft PowerPoint (version 16.23).

The University of Cincinnati Institutional Animal Care and Use Committee (IACUC protocol # 10-01-21-01) approved all housing conditions and collection methods for the use of Astyanax mexicanus in this study. Additionally, this study conformed to the principles outlined in the Basel Declaration and the ARRIVE guidelines.

The datasets generated during and/or analysed during the current study are available in Supplementary Information.

Gross, J. B. The complex origin of Astyanax cavefish. BMC Evol. Biol. 12, 105 (2012).

Article PubMed PubMed Central Google Scholar

Fumey, J. et al. Evidence for late Pleistocene origin of Astyanax mexicanus cavefish. BMC Evol. Biol. 18, 43 (2018).

Article PubMed PubMed Central Google Scholar

Herman, A. et al. The role of gene flow in rapid and repeated evolution of cave-related traits in Mexican tetra, Astyanax mexicanus. Mol. Ecol. 27, 4397–4416 (2018).

Article CAS PubMed PubMed Central Google Scholar

Hubbs, C. L. & Innes, W. T. The first known blind fish of the family Characidae: A new genus from Mexico. Occasional Papers of the Museum of Zoology. (1936).

Elliott, W. R. The Astyanax Caves of Mexico: Cavefishes of Tamaulipas, San Luis Potosí, and Guerrero (Association for Mexican Cave Studies, 2018).

Google Scholar

Hervant, F. & Malard, F. In Encyclopedia of Caves (Third Edition) (eds White, W. B. et al.) 8–15 (Academic Press, 2019).

Chapter Google Scholar

Elliott, W. R. Bad air in caves. Am. Cav. Accidents NSS News 55, 396–397 (1997).

Google Scholar

Ornelas-García, P., Pajares, S., Sosa-Jiménez, V. M., Rétaux, S. & Miranda-Gamboa, R. A. Microbiome differences between river-dwelling and cave-adapted populations of the fish Astyanax mexicanus (De Filippi, 1853). PeerJ (San Francisco, CA) 6, e5906 (2018).

Google Scholar

Rohner, N. et al. Cryptic variation in morphological evolution: HSP90 as a capacitor for loss of eyes in cavefish. Science 342, 1372–1375 (2013).

Article ADS CAS PubMed PubMed Central Google Scholar

Hüppop, K. Oxygen consumption of Astyanax fasciatus (Characidae, Pisces): A comparison of epigean and hypogean populations. Environ. Biol. Fishes 17, 299–308 (1986).

Article Google Scholar

Moran, D., Softley, R. & Warrant, E. J. The energetic cost of vision and the evolution of eyeless Mexican cavefish. Sci. Adv. 1, e1500363 (2015).

Article ADS PubMed PubMed Central Google Scholar

Tang, J. L. Y. et al. The developmental origin of heart size and shape differences in Astyanax mexicanus populations. Dev. Biol. 441, 272–284 (2018).

Article CAS PubMed PubMed Central Google Scholar

Jenny, J. et al. Global spread of hypoxia in freshwater ecosystems during the last three centuries is caused by rising local human pressure. Glob. Change Biol. 22, 1481–1489 (2016).

Article ADS Google Scholar

Díaz, R. J. & Breitburg, D. L. In Fish Physiology (eds Richards, J. G. et al.) 1–23 (Academic Press, 2009).

Google Scholar

Diaz, R. J., Nestlerode, J. & Diaz, M. L. A global perspective on the effects of eutrophication and hypoxia on aquatic biota and water quality. (2004).

Marras, S., Killen, S. S., Claireaux, G., Domenici, P. & McKenzie, D. J. Behavioural and kinematic components of the fast-start escape response in fish: Individual variation and temporal repeatability. J. Exp. Biol. 214, 3102–3110 (2011).

Article PubMed Google Scholar

Sollid, J., De Angelis, P., Gundersen, K. & Nilsson, G. E. Hypoxia induces adaptive and reversible gross morphological changes in crucian carp gills. J. Exp. Biol. 206, 3667–3673 (2003).

Article PubMed Google Scholar

Kajimura, S., Aida, K. & Duan, C. Understanding hypoxia-induced gene expression in early development: In vitro and in vivo analysis of hypoxia-inducible factor 1-regulated zebra fish insulin-like growth factor binding protein 1 gene expression. Mol. Cell. Biol. 26, 1142–1155 (2006).

Article CAS PubMed PubMed Central Google Scholar

Storz, J. F. Hemoglobin: Insights Into Protein Structure, Function, and Evolution (Oxford University Press, 2018).

Book Google Scholar

Baalsrud, H. T. et al. Evolution of hemoglobin genes in codfishes influenced by ocean depth. Sci. Rep. 7, 1–10 (2017).

Article CAS Google Scholar

Affonso, E. G. et al. Blood parameters and metabolites in the teleost fish Colossoma macropomum exposed to sulfide or hypoxia. Comp. Biochem. Phys. C. 133, 375–382 (2002).

CAS Google Scholar

Pan, Y. K., Ern, R., Morrison, P. R., Brauner, C. J. & Esbaugh, A. J. Acclimation to prolonged hypoxia alters hemoglobin isoform expression and increases hemoglobin oxygen affinity and aerobic performance in a marine fish. Sci. Rep. 7, 7834 (2017).

Article ADS PubMed PubMed Central Google Scholar

Hardison, R. C. Evolution of hemoglobin and its genes. CSH. Perspect. Med. 2, a011627 (2012).

ADS Google Scholar

Frey, B. J., Weber, R. E., van Aardt, W. J. & Fago, A. The haemoglobin system of the mudfish, Labeo capensis: Adaptations to temperature and hypoxia. Comp. Biochem. Phys. B. 120, 735–742 (1998).

Article Google Scholar

Sun, H., Li, J., Tang, L. & Yang, Z. Responses of crucian carp Carassius auratus to long-term exposure to nitrite and low dissolved oxygen levels. Biochem. Syst. Ecol. 44, 224–232 (2012).

Article CAS Google Scholar

Silver, R. T. & Gjoni, S. The hematocrit value in polycythemia vera: Caveat utilitor. Leuk. Lymphoma 56, 1540–1541 (2015).

Article PubMed Google Scholar

Reinhart, W. H. The optimum hematocrit. Clin. Hemorheol. Microcirc. 64, 575–585 (2016).

Article CAS PubMed Google Scholar

Filho, D. W. et al. Comparative hematology in marine fish. Comp. Biochem. Phys. A. 102, 311–321 (1992).

Article Google Scholar

Schneider, C. A., Rasband, W. S. & Eliceiri, K. W. NIH Image to ImageJ: 25 years of image analysis. Nat. Methods 9, 671–675 (2012).

Article CAS PubMed PubMed Central Google Scholar

McLaren, C. E., Brittenham, G. M. & Hasselblad, V. Statistical and graphical evaluation of erythrocyte volume distributions. Am. J. Physiol. Heart. C. 252, H857–H866 (1987).

Article CAS Google Scholar

Erckens, W. & Martin, W. Exogenous and endogenous control of swimming activity in Astyanax mexicanus (Characidae, Pisces) by direct light response and by a circadian oscillator II. Features of time-controlled behaviour of a cave population and their comparison to a epigean ancestral form. Z. Nat. C 37, 1266–1273 (1982).

Google Scholar

Duboué, E. R., Keene, A. C. & Borowsky, R. L. Evolutionary convergence on sleep loss in cavefish populations. Curr. Biol. 21, 671–676 (2011).

Article PubMed Google Scholar

Yoshizawa, M. et al. Distinct genetic architecture underlies the emergence of sleep loss and prey-seeking behavior in the Mexican cavefish. BMC Biol. 13, 15 (2015).

Article PubMed PubMed Central Google Scholar

Carlson, B. M. & Gross, J. B. Characterization and comparison of activity profiles exhibited by the cave and surface morphotypes of the blind Mexican tetra, Astyanax mexicanus. Comput. Biochem. Phys. C. 208, 114–129 (2018).

CAS Google Scholar

Scott, G. R. et al. Air breathing and aquatic gas exchange during hypoxia in armoured catfish. J. Comp. Physiol. B. 187, 117–133 (2017).

Article CAS PubMed Google Scholar

Snyder, G. K. & Sheafor, B. A. Red blood cells: Centerpiece in the evolution of the vertebrate circulatory system. Am. Zool. 39, 189–198 (1999).

Article Google Scholar

Lay, P. A. & Baldwin, J. What determines the size of teleost erythrocytes? Correlations with oxygen transport and nuclear volume. Fish Physiol. Biochem. 20, 31–35 (1999).

Article CAS Google Scholar

Wells, R. M. G. & Baldwin, J. Oxygen transport potential in tropical reef fish with special reference to blood viscosity and haematocrit. J. Exp. Mar. Biol. Ecol. 141, 131–143 (1990).

Article Google Scholar

Witeska, M. Erythrocytes in teleost fishes: A review. Zool. Ecol. 23, 275–281 (2013).

Article Google Scholar

Beale, A. et al. Circadian rhythms in Mexican blind cavefish Astyanax mexicanus in the lab and in the field. Nat. Commun. 4, 2769 (2013).

Article ADS PubMed Google Scholar

Santer, R. M. & Walker, M. G. Morphological studies on the ventricle of teleost and elasmobranch hearts. J. Zool. 190, 259–272 (1980).

Article Google Scholar

van der Meer David, L. M. et al. Gene expression profiling of the long-term adaptive response to hypoxia in the gills of adult zebrafish. Am. J. Physiol. Reg. I(289), R1512–R1519 (2005).

Google Scholar

Ju, Z., Wells, M. C., Heater, S. J. & Walter, R. B. Multiple tissue gene expression analyses in Japanese medaka (Oryzias latipes) exposed to hypoxia. Comput. Biochem. Phys. C. 145, 134–144 (2007).

Google Scholar

Tiedke, J., Thiel, R. & Burmester, T. Molecular response of estuarine fish to hypoxia: A comparative study with ruffe and flounder from field and laboratory. PLoS One 9, e90778 (2014).

Article ADS PubMed PubMed Central Google Scholar

van den Thillart, G., Wilms, I., Nieveen, M., Weber, R. E. & Witte, F. Hypoxia-induced changes in hemoglobins of Lake Victoria cichlids. J. Exp. Biol. 221, 25 (2018).

Google Scholar

Rutjes, H. A., Nieveen, M. C., Weber, R. E., Witte, F. & Van den Thillart, G. E. E. J. M. Multiple strategies of Lake Victoria cichlids to cope with lifelong hypoxia include hemoglobin switching. Am. J. Physiol-Reg. I. 293, R1376–R1383 (2007).

CAS Google Scholar

Richards, J. G., Farrell, A. P. & Brauner, C. J. Fish Physiology: Hypoxia (Academic Press, 2009).

Google Scholar

Nilsson, G. E. Gill remodeling in fish—a new fashion or an ancient secret?. J. Exp. Biol. 210, 2403–2409 (2007).

Article PubMed Google Scholar

Teyke, T. Morphological differences in neuromasts of the blind cave fish Astyanax hubbsi and the sighted river fish Astyanax mexicanus. Brain Behav. Evol. 35, 23–30 (1990).

Article CAS PubMed Google Scholar

Yoshizawa, M., Gorički, Š, Soares, D. & Jeffery, W. R. Evolution of a behavioral shift mediated by superficial neuromasts helps cavefish find food in darkness. Curr. Biol. 20, 1631–1636 (2010).

Article CAS PubMed PubMed Central Google Scholar

Schemmel, C. Vergleichende untersuchungen an den hautsinnesorganen ober-und unterirdisch lebender Astyanax–Formen. Z. Morphol. Tiere 61, 255–316 (1967).

Article Google Scholar

Varatharasan, N., Croll, R. P. & Franz-Odendaal, T. Taste bud development and patterning in sighted and blind morphs of Astyanax mexicanus. Dev. Dyn. 238, 3056–3064 (2009).

Article PubMed Google Scholar

Perry, S. F., Jonz, M. G. & Gilmour, K. M. Fish Physiology 193–253 (Elsevier, 2009).

Google Scholar

Zachar, P. C. & Jonz, M. G. Neuroepithelial cells of the gill and their role in oxygen sensing. Resp. Physiol. Neurobiol. 184, 301–308 (2012).

Article CAS Google Scholar

Zachar, P. C., Pan, W. & Jonz, M. G. Characterization of ion channels and O2 sensitivity in gill neuroepithelial cells of the anoxia-tolerant goldfish (Carassius auratus). J. Neurophysiol. 118, 3014–3023 (2017).

Article PubMed PubMed Central Google Scholar

R Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. https://www.R-project.org/. (2019).

Hope, R. Rmisc: Rmisc: Ryan Miscellaneous. https://CRAN.R-project.org/package=Rmisc. (2013).

Hampel, F. R. The breakdown points of the mean combined with some rejection rules. Technometrics 27, 95–107 (1985).

Article MathSciNet MATH Google Scholar

Sarma, P. R. In Clinical Methods: The History, Physical, and Laboratory Examinations (eds Walker, H. K. et al.) (Butterworths, 1990).

Google Scholar

Kassambara, A. ggpubr: 'ggplot2' Based Publication Ready Plots. https://CRAN.R-project.org/package=ggpubr. (2019).

Wickham, H. ggplot2. Wiley Interdiscip. Rev. Comput. Stat. 3, 180–185 (2011).

Article Google Scholar

Download references

The authors wish to thank members of the Gross lab, current and former, for helpful discussions concerning this project and assistance with fish husbandry, particularly Daniel Berning and Alyssa Hamm.

Department of Biological Sciences, University of Cincinnati, Cincinnati, OH, 45221, USA

Tyler E. Boggs, Jessica S. Friedman & Joshua B. Gross

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

T.E.B. and J.B.G. designed experiments, T.E.B. and J.S.F. performed experiments, T.E.B. and J.B.G. analyzed the data and T.E.B. and J.B.G. wrote the main text. All authors reviewed the final manuscript.

Correspondence to Joshua B. Gross.

The authors declare no competing interests.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

Boggs, T.E., Friedman, J.S. & Gross, J.B. Alterations to cavefish red blood cells provide evidence of adaptation to reduced subterranean oxygen. Sci Rep 12, 3735 (2022). https://doi.org/10.1038/s41598-022-07619-0

Download citation

Received: 19 November 2021

Accepted: 15 February 2022

Published: 08 March 2022

DOI: https://doi.org/10.1038/s41598-022-07619-0

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.