Specimen collection is essential for modern science.

Natural history museums are vital repositories of specimens, samples and data that inform about the natural world; this Formal Comment revisits a Perspective that advocated for the adoption of compassionate collection practices, querying whether it will ever be possible to completely do away with whole animal specimen collection.

A set of principles and practical suggestions for equitable fieldwork in biology.

Field biology is an area of research that involves working directly with living organisms in situ through a practice known as "fieldwork." Conducting fieldwork often requires complex logistical planning within multiregional or multinational teams, interacting with local communities at field sites, and collaborative research led by one or a few of the core team members. However, existing power imbalances stemming from geopolitical history, discrimination, and professional position, among other factors, perpetuate inequities when conducting these research endeavors. After reflecting on our own research programs, we propose four general principles to guide equitable, inclusive, ethical, and safe practices in field biology: be collaborative, be respectful, be legal, and be safe. Although many biologists already structure their field programs around these principles or similar values, executing equitable research practices can prove challenging and requires careful consideration, especially by those in positions with relatively greater privilege. Based on experiences and input from a diverse group of global collaborators, we provide suggestions for action-oriented approaches to make field biology more equitable, with particular attention to how those with greater privilege can contribute. While we acknowledge that not all suggestions will be applicable to every institution or program, we hope that they will generate discussions and provide a baseline for training in proactive, equitable fieldwork practices.

Species Delimitation, Phylogenomics, and Biogeography of Sulawesi Flying Lizards: A Diversification History Complicated by Ancient Hybridization, Cryptic Species, and Arrested Speciation.

The biota of Sulawesi is noted for its high degree of endemism and for its substantial levels of in situ biological diversification. While the island's long period of isolation and dynamic tectonic history have been implicated as drivers of the regional diversification, this has rarely been tested in the context of an explicit geological framework. Here, we provide a tectonically informed biogeographical framework that we use to explore the diversification history of Sulawesi flying lizards (the Draco lineatus Group), a radiation that is endemic to Sulawesi and its surrounding islands. We employ a framework for inferring cryptic speciation that involves phylogeographic and genetic clustering analyses as a means of identifying potential species followed by population demographic assessment of divergence-timing and rates of bi-directional migration as means of confirming lineage independence (and thus species status). Using this approach, phylogenetic and population genetic analyses of mitochondrial sequence data obtained for 613 samples, a 50-SNP data set for 370 samples, and a 1249-locus exon-capture data set for 106 samples indicate that the current taxonomy substantially understates the true number of Sulawesi Draco species, that both cryptic and arrested speciations have taken place, and that ancient hybridization confounds phylogenetic analyses that do not explicitly account for reticulation. The Draco lineatus Group appears to comprise 15 species-9 on Sulawesi proper and 6 on peripheral islands. The common ancestor of this group colonized Sulawesi ~11 Ma when proto-Sulawesi was likely composed of two ancestral islands, and began to radiate ~6 Ma as new islands formed and were colonized via overwater dispersal. The enlargement and amalgamation of many of these proto-islands into modern Sulawesi, especially during the past 3 Ma, set in motion dynamic species interactions as once-isolated lineages came into secondary contact, some of which resulted in lineage merger, and others surviving to the present. [Genomics; Indonesia; introgression; mitochondria; phylogenetics; phylogeography; population genetics; reptiles.].

Highly-multiplexed and efficient long-amplicon PacBio and Nanopore sequencing of hundreds of full mitochondrial genomes.

BACKGROUND: Mitochondrial genome sequences have become critical to the study of biodiversity. Genome skimming and other short-read based methods are the most common approaches, but they are not well-suited to scale up to multiplexing hundreds of samples. Here, we report on a new approach to sequence hundreds to thousands of complete mitochondrial genomes in parallel using long-amplicon sequencing. We amplified the mitochondrial genome of 677 specimens in two partially overlapping amplicons and implemented an asymmetric PCR-based indexing approach to multiplex 1,159 long amplicons together on a single PacBio SMRT Sequel II cell. We also tested this method on Oxford Nanopore Technologies (ONT) MinION R9.4 to assess if this method could be applied to other long-read technologies. We implemented several optimizations that make this method significantly more efficient than alternative mitochondrial genome sequencing methods. RESULTS: With the PacBio sequencing data we recovered at least one of the two fragments for 96% of samples (~ 80-90%) with mean coverage ~ 1,500x. The ONT data recovered less than 50% of input fragments likely due to low throughput and the design of the Barcoded Universal Primers which were optimized for PacBio sequencing. We compared a single mitochondrial gene alignment to half and full mitochondrial genomes and found, as expected, increased tree support with longer alignments, though whole mitochondrial genomes were not significantly better than half mitochondrial genomes. CONCLUSIONS: This method can effectively capture thousands of long amplicons in a single run and be used to build more robust phylogenies quickly and effectively. We provide several recommendations for future users depending on the evolutionary scale of their system. A natural extension of this method is to collect multi-locus datasets consisting of mitochondrial genomes and several long nuclear loci at once.

Hummingbird blood traits track oxygen availability across space and time.

Predictable trait variation across environments suggests shared adaptive responses via repeated genetic evolution, phenotypic plasticity or both. Matching of trait-environment associations at phylogenetic and individual scales implies consistency between these processes. Alternatively, mismatch implies that evolutionary divergence has changed the rules of trait-environment covariation. Here we tested whether species adaptation alters elevational variation in blood traits. We measured blood for 1217 Andean hummingbirds of 77 species across a 4600-m elevational gradient. Unexpectedly, elevational variation in haemoglobin concentration ([Hb]) was scale independent, suggesting that physics of gas exchange, rather than species differences, determines responses to changing oxygen pressure. However, mechanisms of [Hb] adjustment did show signals of species adaptation: Species at either low or high elevations adjusted cell size, whereas species at mid-elevations adjusted cell number. This elevational variation in red blood cell number versus size suggests that genetic adaptation to high altitude has changed how these traits respond to shifts in oxygen availability.

Their fates intertwined: diversification patterns of the Asian gliding vertebrates may have been forged by dipterocarp trees.

The repeated evolution of gliding in diverse Asian vertebrate lineages is hypothesized to have been triggered by the dominance of tall dipterocarp trees in the tropical forests of Southeast Asia. These dipterocarp forests have acted as both centres of diversification and climatic refugia for gliding vertebrates, and support most of their extant diversity. We predict similarities in the diversification patterns of dipterocarp trees and gliding vertebrates, and specifically test whether episodic diversification events such as rate shifts and/or mass extinctions were temporally congruent in these groups. We analysed diversification patterns in reconstructed timetrees of Asian dipterocarps, the most speciose gliding vertebrates from different classes (Draco lizards, gliding frogs and Pteromyini squirrels) and compared them with similar-sized clades of non-gliding relatives (Diploderma lizards, Philautus frogs and Callosciurinae squirrels) from Southeast Asia. We found significant declines in net-diversification rates of dipterocarps and the gliding vertebrates during the Pliocene-Pleistocene, but not in the non-gliding groups. We conclude that the homogeneity and temporal coincidence of these rate declines point to a viable ecological correlation between dipterocarps and the gliding vertebrates. Further, we suggest that while the diversification decay in dipterocarps was precipitated by post-Miocene aridification of Asia, the crises in the gliding vertebrates were induced by both events concomitantly.

Bewildering biogeography: Waves of dispersal and diversification across southern Wallacea by bent-toed geckos (genus: Cyrtodactylus).

Bent-toed Geckos, genus Cyrtodactylus, are one of the most diverse terrestrial vertebrate groups, and their range extends from South Asia into Australo-Papua and adjacent Pacific islands. Given the generally high faunal endemism on Wallacean islands, it is rather paradoxical that the diversity in these geckos appears to be so low (21 species in Wallacea, 15 in the Philippines) compared with continental shelf assemblages (>300 species on Sunda + Sahul Shelves + adjacent islands). To determine whether this shortfall was real or an artifact of historical undersampling, we analyzed mitochondrial DNA sequences of hundreds of southern Wallacean samples (Lesser Sundas + southern Maluku). After screening to guide sample selection for target capture data collection, we obtained a 1150-locus genomic dataset (1,476,505 bp) for 119 samples of southern Wallacean and closely related lineages. The results suggest that species diversity of Cyrtodactylus in southern Wallacea is vastly underestimated, with phylogenomic and clustering analyses suggesting as many as 25 candidate species, in contrast to the 8 currently described. Gene exchange between adjacent candidate species is absent or minimal across the archipelago with only one case of > 0.5 migrants per generation. Biogeographical analysis suggests that the hitherto unrecognized diversity is the result of at least three independent dispersals from Sulawesi or its offshore islands into southern Wallacea between 6 and 14 Ma, with one invasion producing small-bodied geckos and the other two or three producing larger-bodied geckos. The smaller-bodied laevigatus group appears to be able to coexist with members of either larger-bodied clade, but we have yet to find members of the two larger-bodied clades occurring in sympatry, suggesting that ecological partitioning or competitive exclusion may be shaping individual island assemblages.

Reference genome of an iconic lizard in western North America, Blainville's horned lizard Phrynosoma blainvillii.

Genome assemblies are increasingly being used to identify adaptive genetic variation that can help prioritize the population management of protected species. This approach may be particularly relevant to species like Blainville's horned lizard, Phrynosoma blainvillii, due to its specialized diet on noxious harvester ants, numerous adaptative traits for avoiding predation (e.g. cranial horns, dorsoventrally compressed body, cryptic coloration, and blood squirting from the orbital sinuses), and status as Species of Special Concern in California. Rangewide decline since the early 20th century, the basis of its conservation status, has been driven mainly by habitat conversion, over-collecting, and invasion of a non-native ant that displaces its native ant prey base. Here, we report on a scaffold-level genome assembly for P. blainvillii as part of the California Conservation Genomics Project (CCGP), produced using Pacific Biosciences HiFi long reads and Hi-C chromatin-proximity sequencing technology. The de novo assembly has 78 scaffolds, a total length of ~2.21 Gb, a scaffold N50 length of ~352 Mb, and BUSCO score of 97.4%. This is the second species of Phrynosoma for which a reference genome has been assembled and represents a considerable improvement in terms of contiguity and completeness. Combined with the landscape genomics data being compiled by the CCGP, this assembly will help strategize efforts to maintain and/or restore local genetic diversity, where interventions like genetic rescue, translocation, and strategic land preservation may be the only means by which P. blainvillii and other low-vagility species can survive in the fragmented habitats of California.

Phylogenomic Analysis Reveals Dispersal-Driven Speciation and Divergence with Gene Flow in Lesser Sunda Flying Lizards (Genus Draco).

The Lesser Sunda Archipelago offers exceptional potential as a model system for studying the dynamics of dispersal-driven diversification. The geographic proximity of the islands suggests the possibility for successful dispersal, but this is countered by the permanence of the marine barriers and extreme intervening currents that are expected to hinder gene flow. Phylogenetic and species delimitation analyses of flying lizards (genus Draco) using single mitochondrial genes, complete mitochondrial genomes, and exome-capture data sets identified 9-11 deeply divergent lineages including single-island endemics, lineages that span multiple islands, and parapatrically distributed nonsister lineages on the larger islands. Population clustering and PCA confirmed these genetic boundaries with isolation-by-distance playing a role in some islands or island sets. While gdi estimates place most candidate species comparisons in the ambiguous zone, migration estimates suggest 9 or 10 species exist with nuclear introgression detected across some intra-island contact zones. Initial entry of Draco into the archipelago occurred at 5.5-7.5 Ma, with most inter-island colonization events having occurred between 1-3 Ma. Biogeographical model testing favors scenarios integrating geographic distance and historical island connectivity, including an initial stepping-stone dispersal process from the Greater Sunda Shelf through the Sunda Arc as far eastward as Lembata Island. However, rather than reaching the adjacent island of Pantar by dispersing over the 15-km wide Alor Strait, Draco ultimately reached Pantar (and much of the rest of the archipelago) by way of a circuitous route involving at least five overwater dispersal events. These findings suggest that historical geological and oceanographic conditions heavily influenced dispersal pathways and gene flow, which in turn drove species formation and shaped species boundaries. [Biogeography; genomics, Indonesia; lizards; phylogeography; reptiles].

Rivers of Indochina as potential drivers of lineage diversification in the spotted flying lizard (Draco maculatus) species complex.

Southeast Asia hosts a rich concentration of biodiversity within multiple biodiversity hotspots. Indochina, a region with remarkably high levels of in situ diversification, possesses five major rivers (Ayeyarwady, Chiang Mai, Mekong, Red, and Salween), several of which coincide with phylogenetic breaks of terrestrial taxa. Draco maculatus possesses a range that stretches across Indochina, which widespread geographic distribution along with potential discrete variation within subspecies alludes to the possibility of this taxon constituting multiple divergent lineages. Using sequence data from three mitochondrial (12S, 16S, and ND2) and three nuclear (BDNF, CMOS, and PNN) genes, we provide the first estimated phylogeny of this hypothesized species complex and examine its phylogeographic architecture with maximum likelihood and Bayes factor delimitation (BFD) approaches. Our results support multiple divergent lineages with phylogenetic breaks coincident with rivers, indicating that river barriers may be contributing to the elevated levels of in situ diversification of Indochina.

Sexual Dichromatism Drives Diversification within a Major Radiation of African Amphibians.

Theory predicts that sexually dimorphic traits under strong sexual selection, particularly those involved with intersexual signaling, can accelerate speciation and produce bursts of diversification. Sexual dichromatism (sexual dimorphism in color) is widely used as a proxy for sexual selection and is associated with rapid diversification in several animal groups, yet studies using phylogenetic comparative methods to explicitly test for an association between sexual dichromatism and diversification have produced conflicting results. Sexual dichromatism is rare in frogs, but it is both striking and prevalent in African reed frogs, a major component of the diverse frog radiation termed Afrobatrachia. In contrast to most other vertebrates, reed frogs display female-biased dichromatism in which females undergo color transformation, often resulting in more ornate coloration in females than in males. We produce a robust phylogeny of Afrobatrachia to investigate the evolutionary origins of sexual dichromatism in this radiation and examine whether the presence of dichromatism is associated with increased rates of net diversification. We find that sexual dichromatism evolved once within hyperoliids and was followed by numerous independent reversals to monochromatism. We detect significant diversification rate heterogeneity in Afrobatrachia and find that sexually dichromatic lineages have double the average net diversification rate of monochromatic lineages. By conducting trait simulations on our empirical phylogeny, we demonstrate that our inference of trait-dependent diversification is robust. Although sexual dichromatism in hyperoliid frogs is linked to their rapid diversification and supports macroevolutionary predictions of speciation by sexual selection, the function of dichromatism in reed frogs remains unclear. We propose that reed frogs are a compelling system for studying the roles of natural and sexual selection on the evolution of sexual dichromatism across micro- and macroevolutionary timescales.

Macroevolutionary Patterns of Sexual Size Dimorphism Among African Tree Frogs (Family: Hyperoliidae).

Sexual size dimorphism (SSD) is shaped by multiple selective forces that drive the evolution of sex-specific body size, resulting in male or female-biased SSD. Stronger selection on one sex can result in an allometric body size scaling relationship consistent with Rensch's rule or its converse. Anurans (frogs and toads) generally display female-biased SSD, but there is variation across clades and the mechanisms driving the evolution of SSD remain poorly understood. We investigated these topics in a diverse family of African treefrogs (Hyperoliidae). Hyperoliids display traits considered rare among amphibians, including sexual dichromatism and protogynous sex change. Using phylogenetic comparative methods, we tested if adult ecology, sexual dichromatism, and sex change were predictors of body size or SSD. We also tested whether hyperoliids displayed allometric interspecific body size scaling relationships. We found a majority of hyperoliid taxa display female-biased SSD, but that adult ecology and sexual dichromatism are poor predictors of sex-specific body size and SSD. Regardless of the groupings analyzed (partitioned by clades or traits), we found support for isometric body size scaling. However, we found that sex change is a significant predictor of SSD variation. Species in the Hyperolius viridiflavus complex, which putatively display this trait, show a significant reduction in SSD and are frequently sexually monomorphic in size. Although protogynous sex change needs to be validated for several of these species, we tentatively propose this trait is a novel mechanism influencing anuran body size evolution. Beyond this association, additional factors that shape the evolution of anuran body size and SSD remain elusive.

Venom Composition in a Phenotypically Variable Pit Viper ( Trimeresurus insularis) across the Lesser Sunda Archipelago.

The genus Trimeresurus comprises a group of venomous pitvipers endemic to Southeast Asia and the Pacific Islands. Of these, Trimeresurus insularis, the White-lipped Island Pitviper, is a nocturnal, arboreal species that occurs on nearly every major island of the Lesser Sunda archipelago. In the current study, venom phenotypic characteristics of T. insularis sampled from eight Lesser Sunda Islands (Flores, Lembata, Lombok, Pantar, Sumba, Sumbawa, Timor, and Wetar) were evaluated via SDS-PAGE, enzymatic activity assays, fibrinogenolytic assays, gelatin zymography, and RP-HPLC, and the Sumbawa sample was characterized by venomic analysis. For additional comparative analyses, venoms were also examined from several species in the Trimeresurus complex, including T. borneensis, T. gramineus, T. puniceus, T. purpureomaculatus, T. stejnegeri, and Protobothrops flavoviridis. Despite the geographical isolation, T. insularis venoms from all eight islands demonstrated remarkable similarities in gel electrophoretic profiles and RP-HPLC patterns, and all populations had protein bands in the mass ranges of phosphodiesterases (PDE), l-amino acid oxidases (LAAO), P-III snake venom metalloproteinases (SVMP), serine proteases, cysteine-rich secretory proteins (CRISP), phospholipases A2 (PLA2), and C-type lectins. An exception was observed in the Lombok sample, which lacked protein bands in the mass range of serine protease and CRISP. Venomic analysis of the Sumbawa venom also identified these protein families, in addition to several proteins of lesser abundance (<1%), including glutaminyl cyclase, aminopeptidase, PLA2 inhibitor, phospholipase B, cobra venom factor, 5'-nucleotidase, vascular endothelial growth factor, and hyaluronidase. All T. insularis venoms exhibited similarities in thrombin-like and PDE activities, while significant differences were observed for LAAO, SVMP, and kallikrein-like activities, though these differences were only observed for a few islands. Slight but noticeable differences were also observed with fibrinogen and gelatin digestion activities. Trimeresurus insularis venoms exhibited overall similarity to the other Trimeresurus complex species examined, with the exception of P. flavoviridis venom, which showed the greatest overall differentiation. Western blot analysis revealed that all major T. insularis venom proteins were recognized by Green Pitviper ( T. albolabris) antivenom, and reactivity was also seen with most venom proteins of the other Trimeresurus species, but incomplete antivenom-venom recognition was observed against P. flavoviridis venom proteins. These results demonstrate significant conservation in the venom composition of T. insularis across the Lesser Sunda archipelago relative to the other Trimeresurus species examined.

Habitat preference modulates trans-oceanic dispersal in a terrestrial vertebrate.

The importance of long-distance dispersal (LDD) in shaping geographical distributions has been debated since the nineteenth century. In terrestrial vertebrates, LDD events across large water bodies are considered highly improbable, but organismal traits affecting dispersal capacity are generally not taken into account. Here, we focus on a recent lizard radiation and combine a summary-coalescent species tree based on 1225 exons with a probabilistic model that links dispersal capacity to an evolving trait, to investigate whether ecological specialization has influenced the probability of trans-oceanic dispersal. Cryptoblepharus species that occur in coastal habitats have on average dispersed 13 to 14 times more frequently than non-coastal species and coastal specialization has, therefore, led to an extraordinarily widespread distribution that includes multiple continents and distant island archipelagoes. Furthermore, their presence across the Pacific substantially predates the age of human colonization and we can explicitly reject the possibility that these patterns are solely shaped by human-mediated dispersal. Overall, by combining new analytical methods with a comprehensive phylogenomic dataset, we use a quantitative framework to show how coastal specialization can influence dispersal capacity and eventually shape geographical distributions at a macroevolutionary scale.

Trait evolution, resource specialization and vulnerability to plant extinctions among Antillean hummingbirds.

Species traits are thought to predict feeding specialization and the vulnerability of a species to extinctions of interaction partners, but the context in which a species evolved and currently inhabits may also matter. Notably, the predictive power of traits may require that traits evolved to fit interaction partners. Furthermore, local abiotic and biotic conditions may be important. On islands, for instance, specialized and vulnerable species are predicted to be found mainly in mountains, whereas species in lowlands should be generalized and less vulnerable. We evaluated these predictions for hummingbirds and their nectar-food plants on Antillean islands. Our results suggest that the rates of hummingbird trait divergence were higher among ancestral mainland forms before the colonization of the Antilles. In correspondence with the limited trait evolution that occurred within the Antilles, local abiotic and biotic conditions-not species traits-correlate with hummingbird resource specialization and the vulnerability of hummingbirds to extinctions of their floral resources. Specifically, hummingbirds were more specialized and vulnerable in conditions with high topographical complexity, high rainfall, low temperatures and high floral resource richness, which characterize the Antillean Mountains. These findings show that resource specialization and species vulnerability to extinctions of interaction partners are highly context-dependent.

Synchronous diversification of Sulawesi's iconic artiodactyls driven by recent geological events.

The high degree of endemism on Sulawesi has previously been suggested to have vicariant origins, dating back to 40 Ma. Recent studies, however, suggest that much of Sulawesi's fauna assembled over the last 15 Myr. Here, we test the hypothesis that more recent uplift of previously submerged portions of land on Sulawesi promoted diversification and that much of its faunal assemblage is much younger than the island itself. To do so, we combined palaeogeographical reconstructions with genetic and morphometric datasets derived from Sulawesi's three largest mammals: the babirusa, anoa and Sulawesi warty pig. Our results indicate that although these species most likely colonized the area that is now Sulawesi at different times (14 Ma to 2-3 Ma), they experienced an almost synchronous expansion from the central part of the island. Geological reconstructions indicate that this area was above sea level for most of the last 4 Myr, unlike most parts of the island. We conclude that emergence of land on Sulawesi (approx. 1-2 Myr) may have allowed species to expand synchronously. Altogether, our results indicate that the establishment of the highly endemic faunal assemblage on Sulawesi was driven by geological events over the last few million years.

Toxic toad invasion of Wallacea: A biodiversity hotspot characterized by extraordinary endemism.

(a) A map of Wallacea showing islands invaded by Duttaphrynus melanostictus in red, islands inhabited by Varanus komodoensis in blue, and localities of genetic samples in yellow points. (b) A D. melanostictus from Lombok Island. (c) Environmental niche model for the Sunda Islands clade of D. melanostictus projected into Wallacea. Green color indicates very high suitability, yellow color indicates high suitability, and orange color indicates moderate suitability.

Process-Based Species Pools Reveal the Hidden Signature of Biotic Interactions Amid the Influence of Temperature Filtering.

A persistent challenge in ecology is to tease apart the influence of multiple processes acting simultaneously and interacting in complex ways to shape the structure of species assemblages. We implement a heuristic approach that relies on explicitly defining species pools and permits assessment of the relative influence of the main processes thought to shape assemblage structure: environmental filtering, dispersal limitations, and biotic interactions. We illustrate our approach using data on the assemblage composition and geographic distribution of hummingbirds, a comprehensive phylogeny and morphological traits. The implementation of several process-based species pool definitions in null models suggests that temperature-but not precipitation or dispersal limitation-acts as the main regional filter of assemblage structure. Incorporating this environmental filter directly into the definition of assemblage-specific species pools revealed an otherwise hidden pattern of phylogenetic evenness, indicating that biotic interactions might further influence hummingbird assemblage structure. Such hidden patterns of assemblage structure call for a reexamination of a multitude of phylogenetic- and trait-based studies that did not explicitly consider potentially important processes in their definition of the species pool. Our heuristic approach provides a transparent way to explore patterns and refine interpretations of the underlying causes of assemblage structure.

Deep genetic structure and ecological divergence in a widespread human commensal toad.

The Asian common toad (Duttaphrynus melanostictus) is a human commensal species that occupies a wide variety of habitats across tropical Southeast Asia. We test the hypothesis that genetic variation in D. melanostictus is weakly associated with geography owing to natural and human-mediated dispersal facilitated by its commensal nature. Phylogenetic and population genetic analyses of mitochondrial and nuclear DNA sequence variation, and predictive species distribution modelling, unexpectedly recovered three distinct evolutionary lineages that differ genetically and ecologically, corresponding to the Asian mainland, coastal Myanmar and the Sundaic islands. The persistence of these three divergent lineages, despite ample opportunities for recent human-mediated and geological dispersal, suggests that D. melanostictus actually consists of multiple species, each having narrower geographical ranges and ecological niches, and higher conservation value, than is currently recognized. These findings also have implications for the invasion potential of this human commensal elsewhere, such as in its recently introduced ranges on the islands of Borneo, Sulawesi, Seram and Madagascar.

Repeated elevational transitions in hemoglobin function during the evolution of Andean hummingbirds.

Animals that sustain high levels of aerobic activity under hypoxic conditions (e.g., birds that fly at high altitude) face the physiological challenge of jointly optimizing blood-O2 affinity for O2 loading in the pulmonary circulation and O2 unloading in the systemic circulation. At high altitude, this challenge is especially acute for small endotherms like hummingbirds that have exceedingly high mass-specific metabolic rates. Here we report an experimental analysis of hemoglobin (Hb) function in South American hummingbirds that revealed a positive correlation between Hb-O2 affinity and native elevation. Protein engineering experiments and ancestral-state reconstructions revealed that this correlation is attributable to derived increases in Hb-O2 affinity in highland lineages, as well as derived reductions in Hb-O2 affinity in lowland lineages. Site-directed mutagenesis experiments demonstrated that repeated evolutionary transitions in biochemical phenotype are mainly attributable to repeated amino acid replacements at two epistatically interacting sites that alter the allosteric regulation of Hb-O2 affinity. These results demonstrate that repeated changes in biochemical phenotype involve parallelism at the molecular level, and that mutations with indirect, second-order effects on Hb allostery play key roles in biochemical adaptation.