Global Environmental Factors Impact the Evolution of Adult Hemoglobins in Squamata Reptiles (Lizards and Snakes) and Terrestrial Turtles

Convergent evolution of oxygen transport mechanisms arises from respiratory proteins adapting to similar environmental pressures. We examined this relationship between adult hemoglobin subunits (Hbs: HBA1, HBAD, HBB1, and HBB2) found in land reptiles (lizards, snakes, and turtles) with their global distribution variables: Altitude, latitude, ambient temperature, and biomass production. We found that biomass was positively associated with the synonymous substitution rate (dS) of HBAD, while it showed the opposite trend for HBB2 in snakes. Additionally, latitude was negatively related to the dS of HBB2 in snakes, but nonsignificant with other Hbs. Altitude was negatively associated with ω = dN/dS of HBA1 and HBAD, whereas temperature showed a similar negative trend with the ω of HBAD across reptiles and in HBB2 of snakes. At amino acid sites, we found most were conserved except for 11 (two near the heme-binding pocket) across Hbs. These fast-changing sites shifted from polar to nonpolar residues, showing a pattern seen in high-altitude mammals. Our results highlight that in reptiles (i) Hbs are diversifying at individual amino acid sites while generally some subunits exhibiting lower ω rates at higher altitudes and hotter temperatures, with the later and higher biomass ecosystems also linked to increases in dS; (ii) HBBs are the most conserved of the Hbs; (iii) latitudinal gradients only show a significant association with the dS of HBB2 in snakes; and (iv) gene conversion events occurred across HBBs in reptiles, which confound their homology assignation, except for snakes that evidenced a single major duplication in their HBBs.