High-altitude adaptation is accompanied by strong signatures of purifying selection in the mitochondrial genomes of three Andean waterfowl

Abstract

Evidence from a variety of organisms points to convergent evolution on the mitochondria associated with a physiological response to oxygen deprivation or temperature stress, including mechanisms for high-altitude adaptation. Here, we examine whether demography and/or selection explains standing mitogenome nucleotide diversity in high-altitude adapted populations of three Andean waterfowl species: yellow-billed pintail (Anas georgica), speckled teal (Anas flavirostris), and cinnamon teal (Spatula cyanoptera). We compared a total of 60 mitogenomes from each of these three duck species (n = 20 per species) across low and high altitudes and tested whether part(s) or all of the mitogenome exhibited expected signatures of purifying selection within the high-altitude populations of these species. Historical effective population sizes (N_e) were inferred to be similar between high- and low-altitude populations of each species, suggesting that selection rather than genetic drift best explains the reduced genetic variation found in mitochondrial genes of high-altitude populations compared to low-altitude populations of the same species. Specifically, we provide evidence that establishment of these three Andean waterfowl species in the high-altitude environment, coincided at least in part with a persistent pattern of negative purifying selection acting on oxidative phosphorylation (OXPHOS) function of the mitochondria. Our results further reveal that the extent of gene-specific purifying selection has been greatest in the speckled teal, the species with the longest history of high-altitude occupancy.

Introduction

Oxidative phosphorylation (OXPHOS) refers to the synthesis of adenosine triphosphate (ATP) in the mitochondria, which is driven by the transfer of protons across the inner mitochondrial membrane via the electron transport chain (ETC). Mitochondrial-encoded OXPHOS and other nuclear-encoded genes involved in the ETC serve a critical function in the production of cellular energy fueled by the TCA cycle through substrate oxidation derived from carbohydrates, fatty acids, and proteins [1]. Consequently, mitochondria are constantly exposed to selective pressures as individuals adapt to the energetic requirements of different environments; and thus, are considered an important driver of speciation [2–4] as well as leading indicators of population divergence for studies involving historical phylogeography [5, 6]. To date, across lineages, many studies have identified adaptive mitochondrial variation associated with a physiological response to changes in oxygen (O₂) availability, temperature, or oxidative stress, including differential mitochondrial gene expression and regulation [7–9], increased or decreased density of mitochondria in skeletal muscle [10, 11], examples of increased or decreased OXPHOS capacity and catalytic efficiency [12, 13], as well as positively-selected genetic variation [14–16]. In birds that have adapted to hypoxic, high-altitude conditions there also is ample evidence of selective pressures on mitogenomic protein coding regions [17–19]. Finally, adaptation to high-altitude environments where O₂ availability is severely limited frequently has been shown to follow predictable paths of parallel or convergent evolution [20–23]. However, while this generally has been shown to be true for genes under strong directional selection like hemoglobin [24], this has not always been shown to be the case for studies of mitochondrial function in which much more variable responses generally have been revealed [13, 25–30].

While direct assessments of mitochondrial function are crucial for establishing causation, measuring of standing genetic variation offers another route to infer the process of natural selection. Standing genetic variation generally refers to the amount of genetic diversity present in a population before or after natural selection has occurred, and as such represents the diversity (or lack-thereof) that selection has to act upon. However, disentangling the roles that demography and selection play shaping standing genetic variation can be complicated, as the former can leave genomic signatures resembling the latter or vice versa, especially for a single non-recombining linkage group such as the mitochondrial DNA (mtDNA). This may be particularly important in the case of recent founder events and/or niche invasions, in which subsequent population expansion leads to a shift in the allele frequency spectrum resembling either positive directional selection, or yet another type of selection, purifying selection [31]. Purifying selection refers to the selective removal of deleterious alleles and is expected to play a pervasive role in maintaining the status quo of biochemical function through stabilizing selection. Like directional selection, purifying selection is also characterized by a reduction in standing genetic variation; thus, it is important to assess possible influences of both demography and selection when establishing whether particular genes or genetic regions show decreased standing variation due to any of these processes.

The Andes are the world’s highest mountains outside of Asia and possess some of the largest high plateaus in the world. In South America, the Altiplano is the largest such plateau extending from northwestern Argentina to southern Peru. The Altiplano contains many wetlands, including Lake Titicaca, that host a large number of waterfowl [32, 33]. Most prominent are the common dabbling duck (Anas & Spatula) species whose high-altitude populations are known to have independently colonized habitats in the Andes from closely-related ancestral low-altitude lineages [28]. To date, research on this taxonomic group has identified genetic convergence in both the hypoxia-inducible factor (HIF) pathway, as well as hemoglobin (Hb) genes [20, 22, 34, 35]. However, outside of these well-studied genes, genetic variation associated with other physiological mechanisms for high-altitude living remains relatively unexplored in this group. Here, we seek to determine whether convergent evolution resulting from either directional selection or purifying selection has played out in the mitogenomes of these same species of ducks due to its role in OXPHOS as described above.

We compared the standing genetic diversity of each low- and high-altitude population of three Andean duck species using ten mitogenomes per population. The species compared include the yellow-billed pintail (Anas georgica), speckled teal (Anas flavirostris), and cinnamon teal (Spatula cyanoptera; syn. Anas cyanoptera), which have been the focus of previous physiological studies spanning O₂ transport [20, 22, 28, 36–41]. These paired low- and high-altitude population pairs provide natural replicates to study the genetic attributes linked to demography caused by historical founder events involving colonization of cold, high-altitude, hypoxic environments from ancestral low-altitude environments. These pairwise comparisons also represent contrasting depths of divergence times between the populations. Whereas respective low- and high-altitude populations of yellow-billed pintail and cinnamon teal diverged more recently, these mtDNA data reveal that speckled teal represent much deeper divergence by approximately two orders of magnitude [42–44]. Thus, these species allow us to infer potential differences caused by the effect of demography versus selective forces on standing genetic diversity among the population-pairs, including how historical colonization of the high-altitude environment impacted mitogenomic diversity either through directional selection or purifying selection. Ultimately, we find little evidence that variable demographic processes have shaped mitogenome diversity in these species. Instead, we present evidence consistent with the prevalence of stronger purifying selection in high-altitude populations compared to low-altitude populations of each species, with the strongest evidence for purifying selection occurring in the species that has occupied the high-altitude environment for the longest period of evolutionary time.

Materials and methods

Specimen collection and DNA extraction

A total of 60 individuals composed of 10 low- and 10 high-altitude samples across three species of ducks were used [22] (Fig 1). For the cinnamon teal, individuals from low-altitude populations are the S. c. cyanoptera subspecies (n = 10; elevation range 7–23 m) and from high-altitude are the S. c. orinoma subspecies (n = 10; elevation range 3533–3,871 m) [43]. For the speckled teal, individuals from low-altitude populations are the A. f. flavirostris subspecies (n = 10; elevation range 77–860 m) and from high-altitude are the A. f. oxyptera subspecies (n = 10; elevation range 3,211–4,405 m) [44]. For the yellow-billed pintail, individuals from both populations are taxonomically identified as the same subspecies Anas georgica spinicauda. A total of 20 yellow-billed pintails were collected from low- (n = 10; elevation range 292–914 m) and high-altitude (n = 10; elevation range 3,332–4,070 m) [42]. Genomic DNA was extracted from tissue using a DNeasy Tissue Kit (Qiagen, Valencia, California, USA) following manufacturers protocols.

[Figure omitted. See PDF.]

Species distribution map layers were obtained from [93].

Target-enrichment sequencing

We utilized in-solution target capture to selectively enrich libraries for the complete mitochondrial genome prior to NGS sequencing [45]. All steps of the process were performed by MYcroarray (Ann Arbor, MI). In short, a custom MYbaits® biotinylated ssRNA target capture baitset was designed with 1,359 120mer probes at 2x tiling density from five previously published Anas mitochondrial genomes (A. acuta, NC_024631.1; A. crecca, KC771255.1; A. formosa, NC_015482.1, A. platyrhynchos NC_009684/EU009397.1; and A. poecilrhyncha NC_022418.1). DNA samples were then subjected to capture library synthesis and preparation. The prepared libraries were hybridized with the custom mitochondrial baits. Following hybridization, target regions were purified on magnetic beads followed by post-hybridization amplification to add indexing sequences. Ultimately, 250–300 bp insert sizes were targeted with sequencing of a single multiplexed library performed on an Illumina Hi-Seq platform using paired-end (100 bp) chemistry. Illumina reads were deposited in NCBI’s Sequence Read Archive (SRA; http://www.ncbi.nlm.nih.gov/sra; SRA data PRJNA508951).

Mitochondrial genome assembly and annotation

All sequences were de-multiplexed by individual with adapters trimmed and quality filtered (Q < 30) post-sequencing using the FASTX-Toolkit v. 0.0.13 package [46]. Trimming was performed utilizing fastq-clipper (AGATCGGAAGAGC), and remaining sequences were then filtered by length and quality using fastq-quality filter (reads <20 bp, and Q < 30). The mitochondrial genome for each individual was assembled with MITObim v. 1.7 [47] using as a reference the full mitochondrial genome of the mallard (Anas platyrynchnos; EU009397.1), which is approximately equidistant from all study species to avoid a potential source of bias. For genome reconstruction, the “quick” and “trim” options were used, with an iteration limit set to 30. Mapping results were aligned with MAFFT [48] and annotated in Geneious v. 8.1.6 [49] based on 90% threshold similarity to the mallard reference. Protein coding, tRNA, and rRNA genes were extracted from these annotations and then subjected to individual analyses.

Population genetic and selection analyses

All comparisons were made between low- and high-altitude populations within species. For each species, we calculated pairwise F_ST [50] for all variable sites for the full mitochondrial genomes as well as individual gene regions—protein coding, and non-protein coding (D-Loop, tRNAs, rRNA). This was done for each of the three species using Arlequin v. 3.5 [51] and MEGA v. 7 [52]. We also calculated Tajima’s D [53, 54], which compares the mean number of pairwise differences and number of segregating sites (e.g., Tajima’s θ vs. Watterson’s θ), scaled as such so that each are expected to be the same in a neutrally evolving population of constant size. Where this is not the case, D < 0 reveals a significant excess of high frequency polymorphisms and can be indicative of a recent selective sweep caused by directional selection or purifying selection (selection) and/or a founder event followed by population expansion (demography). D > 0 by contrast reveals a significant excess of both low- and high-frequency polymorphism and generally is indicative of either balancing selection (selection) or population contraction (demography).

Next, we utilized codon-based Z-tests to compare the relative abundance of synonymous and nonsynonymous substitutions in these mitochondrial sequences, where Z = (d_N−d_S) / √ (Var(d_S) + Var(d_N)). The null hypothesis (H₀) is that the ratio is equal (d_N = d_S), whereas d_N > d_S indicates directional (positive) selection and d_N < d_S indicates purifying (negative) selection. These analyses were performed estimating variance via a bootstrap method (500 replications) as implemented in MEGA v. 7 [54, 55].

To examine evidence for selection for specific codons, we used the SLAC (Single Likelihood Ancestor Counting) method, MEME (Mixed Effects Model of Evolution), and FUBAR (A Fast, Unconstrained Bayesian AppRoximation for Inferring Selection) through the program HyPhy v. 2.5 [56]. SLAC also assesses selection using the nonsynonymous/synonymous ratio (dN/dS), which is defined as ω, with values of ω < 1, ω = 1, and ω > 1 also signifying negative/purifying selection, neutral evolution, and positive/directional selection, respectively. However, these modules estimate ω at every codon in the alignments and report which codons show evidence of positive or negative selection, based on the expected and observed numbers of synonymous and nonsynonymous substitutions as inferred using maximum-likelihood [57]. Unlike SLAC, MEME allows ω to vary across codons as well as across branches of the phylogeny, allowing it to detect a small proportion of branches that are evolving under positive directional selection or under episodic selection [58]. To avoid a high false-positive rate, due to the reduced number of sequences, sites with P values < 0.1 for both models were considered significant [57]. Finally, FUBAR uses a Bayesian approach to infer nonsynoymous (d_N) and synonymous (d_S) substitution rates on a per-site basis for a given coding alignment and corresponding phylogeny, and assumes that the selection pressure for each site is constant along the entire phylogeny [59]. To avoid a high false-positive rate, due to the reduced number of sequences, sites with posterior probabilities > 0.9 for both models were considered significant.

Demographic analyses ‐ population size fluctuation

Changes in effective population size (N_e) across time were inferred using Bayesian Skyline plots as implemented in BEAST v. 1.8 [60]. We ran linear skyline models with a strict molecular clock, and optimum base-pair substitution models determined for each dataset based on Bayesian Information Criterion (BIC) scores as estimated in MEGA v. 10 [52]. Multiple chains were run for 100 million steps yielding effective sample sizes (ESS) of at least 200. The first 10% of runs were discarded as “burn-in”. All operators (transitional kernals or ‘moves’) in the Markov chain Monte Carlo were optimized automatically. Results of the analyses were illustrated using Tracer v. 1.5 [61].

Divergence between low- and high-altitude populations

Finally, we calculated pairwise genetic divergence between each pair of low- and high-altitude populations. This was done for the full mtDNA and each protein coding gene listed above using the Tamura-Nei [62] model of substitution in MEGA v. 10 [52]. The purpose of this analysis was to compare relative divergence times between low- and high-altitude populations among the three species and therefore provide some perspective on the relative duration of occupancy of the high-altitude environment as recently examined in the context of biochemical changes in pathways of energy metabolism [28], in which species that had occupied the highlands for longer periods of evolutionary time had gradually shifted activities of key enzymes related to OXPHOS capacity relative to newcomers to the highlands.

Results

Assembly features and mitogenome organization

We attained an average sequencing depth of 232x across mitogenomes for all species. The final consensus genome sizes per species were 16,601bp (cinnamon teal), 16,602 bp (speckled teal), and 16,616 bp (yellow-billed pintail), respectively (S1 Table). Across these species, the light strand includes NADH6 and 8 tRNAs, whereas the heavy strand includes the remaining 12 protein coding genes, 14 tRNAs and the 12S and 16S rRNAs. The non-coding control region (D-loop) in these species is located between ND6/tRNA-Glu and the tRNA-Phe/12S ribosomal RNA genes similar to other waterfowl species [63]. Overall composition and arrangement of the mitochondrial genome is typical of birds and matches other Anseriformes species ranging from 16,608–16,594 bp in size [64–67].

A maximum-likelihood tree of all 60 individuals including the 6 outgroups (Fig 2) was constructed using GARLI [68, 69], in order to show general phylogenetic relationships between the species/populations. Outgroups for this analysis included A. acuta (NC_024631.1), A. clypeata (NC_028346.1), A. crecca (KC771255.1), A. formosa (NC_015482.1), A. platyrhynchos (EU009397.1), and A. poecilrhyncha (NC_022418.1).

[Figure omitted. See PDF.]

Numbers at nodes represent percent consensus support; “high” = individuals from high-altitude populations, “low” = individuals from low-altitude populations.

Demographic analyses ‐ population size fluctuation

Plotting the median and 95% confidence intervals revealed largely stable estimates of N_e across time for each analyzed low- and high-altitude population of speckled teal and cinnamon teal (Fig 3A and 3B). However, whereas speckled and cinnamon teal show almost constant effective population size through time, yellow-billed pintails appear to have more recently experienced exponentially increasing numbers, particularly in the case of the low-altitude yellow-billed pintail population (Fig 3C). The optimum base-pair substitution models for these analyses were determined to be either HKY85 or the similar TrN93 model.

[Figure omitted. See PDF.]

Bayesian skyline plots for low- and high-altitude populations of (A) speckled teal, (B) cinnamon teal, and (C) yellow-billed pintail.

Population genetic analyses

Comparisons within species between low- and high-altitude populations of cinnamon teal and yellow-billed pintail revealed extremely low mitochondrial differentiation, whereas comparisons within speckled teal revealed much higher differentiation for most mitochondrial gene regions, corresponding to approximately an order of magnitude deeper genetic divergence (Fig 4 and S2 Table), therefore reflecting a longer period of high-altitude occupancy by speckled teal.

[Figure omitted. See PDF.]

For the cinnamon teal, the overall F_ST between low- and high-altitude was 0.02, with 64 segregating sites. For the yellow-billed pintail, F_ST between low- and high -altitude was 0.00, with 139 segregating sites. In contrast, F_ST between low- and high -altitude speckled teal was 0.71, with 206 segregating sites. In speckled teal, there were 119 polymorphisms that were fixed between low- and high-altitude populations, whereas no fixed differences occurred in the other two species. There were few instances of non-synonymous changes between high and low-altitude populations of cinnamon teal or yellow-billed pintail, and none with significant F_ST values. In cinnamon teal, there was a non-synomymous change in ND1 (Met160Val), ND4 (Ala163Thr), and none in yellow-billed pintail. Furthermore, SNPs that were significantly different in speckled teal (P < 0.05) were present in both noncoding (D-Loop, l-rRNA, tRNA-Ala, tRNA-Thr) and a large majority of protein coding genes; specifically, there were 5 instances of various non-synonymous changes resulting a high degree of fixation (F_ST > 0.75) in the speckled teal, specifically in ATP6 (Tyr13His), ND1 (Met325Thr), ND2 (Ile231Leu), ND3 (Met107Val), and ND5 (Ile128Val) (S3 Table). There were additional nonsynonymous changes with a moderate degree of fixation in speckled teal, including COX1 (Leu401Phe), ND2 (Met245Ile), ND3 (Pro76Leu, Val84Ile, Val107Met), and ND5 (Arg476Thr) (S3 Table).

Tests for selection across the mitochondrial genome

Selection was also assayed at the level of specific codons using HyPhy [56], which uses ω (d_N/d_S ratio). The intra-species SLAC analysis identified no credible instances of positive selection. The MEME analysis showed no evidence of episodic diversifying selection across any of the three species, whereas FUBAR showed extensive evidence of purifying selection in all three species across all genes (S4 Table). This same pattern arises when all three species data was combined, with widespread evidence of purifying selection, and no sites significantly under positive selection. Overall, the results from the three analyses (SLAC, MEME, FUBAR) identified many more examples of purifying selection, as compared to positive selection in all three species (S4 Table).

Applying standard neutrality tests to the protein coding genes, both high- and low-altitude populations across species had significantly negative Tajima’s D estimates (Figs 5 and 6 and S5–S7 Tables). Specifically, cinnamon teal and yellow-billed pintail showed an excess of low-frequency SNPs, as well as paucity of high-frequency SNPs. By contrast, speckled teal showed an excess of low- and high-frequency SNPs segregating between the low- and high-altitude populations. Thus, the speckled teal showed a strikingly different site-frequency spectrum characteristic of a more deeply diverged pair of populations.

[Figure omitted. See PDF.]

[Figure omitted. See PDF.]

Finally, the codon-based Z-test was used to test hypotheses about the type of selection (neutral, positive, purifying) predominating between low- and high-altitude populations. For all three species, the Z-test for individual protein coding genes showed evidence for either neutrality or purifying selection (Table 1). None of the analyses in any of the three species yielded significant results for positive selection. However, the low-altitude speckled teal population showed significantly fewer genes under purifying selection, relative to the high-altitude speckled teal population (10; X² = 7.539, P < 0.01). This pattern was not seen in either the cinnamon teal or the yellow-billed pintail (P > 0.05). However, comparing the incidence of sites under purifying selection across all low-altitude populations to all high-altitude populations showed evidence of more widespread purifying selection (X² = 6.477, P < 0.01). In addition, COX2 and ND1 genes appear to be under purifying selection in the high-altitude populations across all three species, whereas the ND4-6 genes appear to be under purifying selection in both low- and high-altitude populations. Overall, neutrality tests with different sensitivities to demographic effects showed significant purifying selection in the protein coding OXPHOS genes of all high-altitude populations (Figs 5 and 6 and S5–S7 Tables).

[Figure omitted. See PDF.]

Because tRNA and rRNA genes do not produce an amino acid sequence, quantifying non-synonymous or synonymous mutations is not possible; thus, selection could only be assessed using Tajima’s D. Both tRNA and rRNA regions showed evidence of purifying selection and/or population expansion as compared to both the non-coding D-loop and protein-coding regions (S7 Table). Otherwise, each high-altitude population of each species showed patterns consistent with purifying selection.

Discussion

In this study, we present a comparative analysis of mitogenome evolution in low- and high-altitude populations of three Andean waterfowl species that colonized the high Andes at different times in history. We performed various analyses aimed to disentangle the relative contributions of demographic and selective processes on the extant standing genetic diversity of each respective population’s mitogenomes. First, we found few or no observable changes in effective population size through time when comparing each respective high- versus low-altitude population within each of the three duck species (Fig 3). All populations within species were similar in this regard. Importantly, the lack of significant variation in effective population size between low- and high-altitude populations through time for each of the six lineages suggests that (1) all three species of these waterfowl did not experience strong reductions in genetic diversity following historical founder events in which the Andes were colonized from low-altitude ancestral populations, and (2) high-altitude populations continue to maintain their effective population size similar to the ancestral size. Given the lack of observable changes in effective population size between populations within species, we conclude that reductions in mitogenome variation of high-altitude populations have more likely been the result of stronger and more frequent instances of purifying selection than experienced by their low-altitude counterparts. In particular, we find evidence of the strongest purifying selection on the OXPHOS genes of the high-altitude populations, likely to be the result of selection for more efficient metabolism in high-altitude habitats [28, 70–72].

Patterns of selection in the three Andean duck species

For the cinnamon teal and yellow-billed pintail, recent divergence has resulted in little or no difference between the high- and low-altitude populations in polymorphism across the entirety of the mitochondrial genome, which is consistent with their more recent demographic history in the highlands (Fig 4). However, the mitochondrial genomes of both species show evidence for purifying selection overall (Table 1 and Figs 5 and 6), with the strongest apparent negative selection being exerted on the coding regions for cinnamon teal, and noncoding regions (tRNAs, rRNAs) for the yellow-billed pintail. For the cinnamon teal, the high-altitude population is considered a separate subspecies, due to striking differences in plumage patterns and body size [73, 74], whereas low- and high-altitude populations of the yellow-billed pintail are classified as the same subspecies [42].

In contrast, the speckled teal showed the strongest evidence for purifying selection on standing mitochondrial genetic diversity. Despite a deeper divergence [44] between the high- and low-altitude populations, mitochondrial genome divergence was extreme (F_ST = 0.71) compared to nuclear DNA (F_ST = 0.06). In this species, our results suggest a pattern of non-neutral processes playing a role in divergence, which mirrors previous research on the evolution of Hb-O₂ affinity across each of these species [20, 75, 76], as well as patterns of variation consistent with strong selection on the HIF pathway [22].

Altogether, our results also show a significant role for purifying selection in both low- and high-altitude populations of all three species, but especially in protein-coding regions of high-altitude populations (Figs 5 and 6 and Table 1, X² = 6.477, P < 0.01). This implies that adaptation to hypoxic environments in these duck species has been demarcated by differences in the strength of purifying selection—“stronger” at high altitude (i.e., more significant instances) vs. “weaker” at low altitude (i.e., fewer significant instances). This contrasting and predictable pattern of selection suggests that the mitochondrial genome has played a key role in adaptation to high-altitude, hypoxic environments. Adaptation of nuclear and mitochondrial proteins at the level of respiratory function is thought to be stringent because of how tightly intertwined the functions of the mitochondrial and nuclear genomes are [77–80]. The compatibility between OXPHOS mitochondrial-nuclear subunits is particularly apparent in prior documented examples of population-level divergence, where it has effectively created isolating barriers leading to population structuring and speciation [81–84]. Thus, given the functional importance of the genes encoded by the mitochondrial genome, selection for optimal catalytic capacity/regulatory efficiently is expected to lead to stronger patterns of purifying selection.

Mitochondrial efficiency is rooted in the relationship of the mitochondrially-encoded subunits to their nuclear-encoded subunits, and specifically the coupling of ATP production by chemiosmosis to the electrochemical proton gradient established by the electron transport system [85]. This mito-nuclear association is thought to be coevolved with both nuclear and mitochondrial genomes evolving complementary changes in the other to ensure correct mitochondrial functioning [79, 86]. Mitochondrial efficiency has been shown to underlie adaptations to extreme environments, like at high-altitude [12, 71], as well as strongly locomotive species such as the waterfowl studied here [87–89]. However, OXPHOS is not perfectly coupled because protons can leak across the inner membrane independent of ATP synthesis [90]. Uncoupling reduces the membrane potential, leading to waste and a reduction in ATP production efficiency [91]. However, uncoupling is also important for thermogenesis, and it also has importance for reactive oxygen species (ROS) production. Since ROS production is positively correlated with the steepness of the proton gradient, uncoupling caused by proton leak can minimize ROS production and therefore prevent oxidative damage [90], which also may be a key factor influencing evolution at high altitude [25, 30].

Recent biochemical assessment of ATP production in high-altitude duck species showed respiratory capacity in high-altitude ducks was significantly associated with elevated activities of mitochondrial enzymes, and oxidative capacity [13, 28], whereupon it was suggested that high-altitude ducks may take advantage of O₂-costly fuels (i.e., lipids), provided that sufficient tissue O₂ supply is maintained by evolved or plastic changes throughout the O₂ transport pathway [20, 92]. Thus, evolved changes in mitochondrial catalytic efficiency could be the reason behind such pervasive signature of purifying selection in high-altitude populations.

Conclusions

The mitochondria is essential for producing energy in the form of ATP through the electron transport chain and its corresponding OXPHOS related complexes and enzymatic pathways, which also involve thermogenesis and ROS production. With these roles, the mitochondrial genome is frequently associated with adaptation to various environments including extreme cold environments and hypoxic environments. We assessed the role of standing genetic variation in the mitochondrial genome as it potentially facilitated adaptation to hypoxic, high-altitude environments in three Andean duck species. Based on the observed relative contributions of demography related to changes in effective population size in these species, we argue that purifying selection has played a strong role at high altitude given the consistency of results across multiple taxa with varying degrees of similarity and difference in demographic history.

Specifically, the results from the three study species suggest that there is selection on the mitogenome across species, but that the same biochemical mechanisms do not always result from the same environmental pressures, as illustrated by two of the same species, in which previous work showed evidence of both short- and long-term adaptation through the hypoxia-inducible factor pathway [22] and hemoglobin complexes [20, 75, 76]. Thus, species invade a new habitat potentially using other mechanisms but then adapt to it—i.e., speckled teal has had enough time for selection to show signs of distinct biochemical pathways, whereas more recently arrived species have not [28].

These three Andean waterfowl species successfully colonized and, are thriving in the same high-altitude environment, with purifying selection acting to maintain the OXPHOS unit’s ability to operate optimally. The signature of purifying selection on the mitogenome is especially apparent for the more deeply diverged speckled teal with the longer history of occupancy at high altitude. Specifically, a substantial increase in purifying selection resulting in decreased nucleotide diversity in this species, is likely due not only to the invasion of high-altitude niches but also due to elapsed time at altitude. These results open up additional questions associated with how purifying selection of this type influences mito-nuclear interactions, plus a general need for a greater understanding of the physiology and biochemistry associated with mitochondria function in this extreme environment. Finally, it should be noted that achieving a better understanding of the unique biochemistry and physiology of these species in this extreme environment should also be beneficial to the long-term conservation of these high Andean species, especially as factors related to mitochondrial function relate to capacity for movement and thermogenesis in a rapidly changing climate where suitable wetlands may be shifting.

Supporting information

S1 Table. Sample and assembly statistics associated with the mitochondrial genome assemblies of the three Andean duck species.

https://doi.org/10.1371/journal.pone.0294842.s001

(XLSX)

S2 Table. Tamura-Nei genetic distance between low- and high-altitude populations of the three Andean duck species.

https://doi.org/10.1371/journal.pone.0294842.s002

(XLSX)

S3 Table. Significant F_ST values in the mitochondrial genome of speckled teal.

https://doi.org/10.1371/journal.pone.0294842.s003

(XLSX)

S4 Table. Significant instances of selection (P > 0.9) across sites using FUBAR (Fast, Unconstrained Bayesian AppRoximation) for each of the three species: Cinnamon Teal (CT), Speckled Teal (ST), Yellow Billed Pintail (YBP). Prob [P] is the posterior probability.

https://doi.org/10.1371/journal.pone.0294842.s004

(XLSX)

S5 Table. Tajima’s D statistics calculated in MEGA for mitochondrial gene regions of three Andean duck species.

https://doi.org/10.1371/journal.pone.0294842.s005

(XLSX)

S6 Table. Tajima’s D statistics for each OXPHOS complex within the mitochondrial genome.

https://doi.org/10.1371/journal.pone.0294842.s006

(XLSX)

S7 Table. Tajima’s D statistics for the non-coding elements (tRNA, rRNA and dLoop) within the mitochondrial genome.

https://doi.org/10.1371/journal.pone.0294842.s007

(XLSX)

S1 Fig.

https://doi.org/10.1371/journal.pone.0294842.s008

(JPG)

Acknowledgments

We thank the many people and provincial and federal governments in Argentina, Peru, and Bolivia who assisted us with fieldwork for many years. Funding was provided by Alaska EPSCoR (NSF EPS-0092040, EPS-0346770), the National Science Foundation (DEB-0444748 and IOS-0949439), Frank M. Chapman Fund at the American Museum of Natural History, and the Kushlan Endowment for Waterbird Biology and Conservation at the University of Miami. AMG was supported by an UM College of Arts and Sciences Maytag Fellowship and Graduate Summer Research Fellowship. This work was supported in part by the high-performance computing and data storage resources operated by the Research Computing Systems Group at the University of Alaska Fairbanks, Geophysical Institute. AMG and KGM conceived the idea/design/experiment, AMG, PL and REW performed the experiments and analyzed the data, AMG, PL, REW and KGM wrote the paper. Specimen collections were carried out under University of Alaska Fairbanks IACUC protocols #02–01 and #05–05.

Citation: Graham AM, Lavretsky P, Wilson RE, McCracken KG (2024) High-altitude adaptation is accompanied by strong signatures of purifying selection in the mitochondrial genomes of three Andean waterfowl. PLoS ONE 19(1): e0294842. https://doi.org/10.1371/journal.pone.0294842

About the Authors:

Allie M. Graham

Roles: Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Validation, Writing – original draft, Writing – review & editing

E-mail: [email protected]

Affiliation: Eccles Institute for Human Genetics, University of Utah, Salt Lake City, UT, United States of America

Philip Lavretsky

Roles: Formal analysis, Methodology, Software, Validation

Affiliation: Department of Biological Sciences, University of Texas at El Paso, El Paso, TX, United States of America

Robert E. Wilson

Roles: Data curation, Formal analysis, Investigation, Methodology, Resources, Writing – review & editing

Affiliation: School of Natural Resources and Nebraska State Museum, University of Nebraska–Lincoln, Lincoln, NE, United States of America

ORICD: https://orcid.org/0000-0003-1800-0183

Kevin G. McCracken

Roles: Conceptualization, Data curation, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Writing – original draft, Writing – review & editing

Affiliations: Department of Biology, University of Miami, Coral Gables, FL, United States of America, Department of Marine Biology and Ecology, Rosenstiel School of Marine, Atmospheric, and Earth Science, University of Miami, Miami, FL, United States of America, Human Genetics and Genomics, University of Miami Miller School of Medicine, Miami, FL, United States of America, University of Alaska Museum, University of Alaska Fairbanks, Fairbanks, AK, United States of America

ORICD: https://orcid.org/0000-0002-1477-8839

[/RAW_REF_TEXT]

References

1. Ballard J, Melvin R. Linking the mitochondrial genotype to the organismal phenotype. Mol Ecol. 2010;19(8):1523–39. pmid:20345689

2. Gershoni M, Templeton AR, Mishmar D. Mitochondrial bioenergetics as a major motive force of speciation. BioEssays. 2009;31(6):642–50. pmid:19408245

3. Das J. The role of mitochondrial respiration in physiological and evolutionary adaptation. BioEssays. 2006;28(9):890–901. pmid:16937356

4. Breton S, Milani L, Ghiselli F, Guerra D, Stewart DT, Passamonti M. A resourceful genome: updating the functional repertoire and evolutionary role of animal mitochondrial DNAs. Trends Genet. 2014;30(12):555–64. pmid:25263762

5. UniProt: the universal protein knowledgebase. Nucleic Acids Res. 2017;45(D1):D158–D69. pmid:27899622

6. Zink RM, Barrowclough GF. Mitochondrial DNA under siege in avian phylogeography. Mol Ecol. 2008;17(9):2107–21. pmid:18397219

7. Gracey AY, Fraser EJ, Li W, Fang Y, Taylor RR, Rogers J, et al. Coping with cold: an integrative, multitissue analysis of the transcriptome of a poikilothermic vertebrate. Proc Natl Acad Sci U S A. 2004;101(48):16970–5. pmid:15550548

8. Whitehead A, Crawford DL. Neutral and adaptive variation in gene expression. Proc Natl Acad Sci U S A. 2006;103(14):5425–30. pmid:16567645

9. Scott GR, Elogio TS, Lui MA, Storz JF, Cheviron ZA. Adaptive modifications of muscle phenotype in high-altitude deer mice are associated with evolved changes in gene regulation. Mol Biol Evol. 2015:msv076. pmid:25851956

10. Hoppeler H, Vogt M, Weibel ER, Flück M. Response of skeletal muscle mitochondria to hypoxia. Exp Physiol. 2003;88(1):109–19. pmid:12525860

11. Scott GR, Egginton S, Richards JG, Milsom WK. Evolution of muscle phenotype for extreme high altitude flight in the bar-headed goose. Proc. R. Soc. B: Biol. Sci. 2009;276(1673):3645–53.

12. Zhang Z-Y, Chen B, Zhao D-J, Kang L. Functional modulation of mitochondrial cytochrome c oxidase underlies adaptation to high-altitude hypoxia in a Tibetan migratory locust. Proc. R. Soc. B: Biol. Sci. 2013;280(1756):20122758.

13. Dawson NJ, Ivy CM, Alza L, Cheek R, York JM, Chua B, et al. Mitochondrial physiology in the skeletal and cardiac muscles is altered in torrent ducks, Merganetta armata, from high altitudes in the Andes. J Exp Biol. 2016;219(23):3719–28. pmid:27618861

14. Hassanin A, Ropiquet A, Couloux A, Cruaud C. Evolution of the mitochondrial genome in mammals living at high altitude: new insights from a study of the tribe Caprini (Bovidae, Antilopinae). J Mol Evol. 2009;68(4):293–310. pmid:19294454

15. Yu L, Wang X, Ting N, Zhang Y. Mitogenomic analysis of Chinese snub-nosed monkeys: evidence of positive selection in NADH dehydrogenase genes in high-altitude adaptation. Mitochondrion. 2011;11(3):497–503. pmid:21292038

16. Melo-Ferreira J, Vilela J, Fonseca MM, da Fonseca RR, Boursot P, Alves PC. The elusive nature of adaptive mitochondrial DNA evolution of an arctic lineage prone to frequent introgression. Genome Biol Evol. 2014;6(4):886–96. pmid:24696399

17. Cai Q, Qian X, Lang Y, Luo Y, Xu J, Pan S, et al. Genome sequence of ground tit Pseudopodoces humilis and its adaptation to high altitude. Genome Biol. 2013;14(3):1. pmid:23537097

18. Toews DP, Mandic M, Richards JG, Irwin DE. Migration, mitochondria, and the yellow‐rumped warbler. Evolution. 2014;68(1):241–55. pmid:24102562

19. Scott GR, Schulte PM, Egginton S, Scott AL, Richards JG, Milsom WK. Molecular evolution of cytochrome c oxidase underlies high-altitude adaptation in the bar-headed goose. Mol Biol Evol. 2011;28(1):351–63. pmid:20685719

20. Natarajan C, Projecto-Garcia J, Moriyama H, Weber RE, Muñoz-Fuentes V, Green AJ, et al. Convergent evolution of hemoglobin function in high-altitude Andean waterfowl involves limited parallelism at the molecular sequence level. PLoS Genet. 2015;11(12):e1005681. pmid:26637114

21. Zhu X, Guan Y, Signore AV, Natarajan C, DuBay SG, Cheng Y, et al. Divergent and parallel routes of biochemical adaptation in high-altitude passerine birds from the Qinghai-Tibet Plateau. Proc Natl Acad Sci U S A. 2018;115(8):1865–70. pmid:29432191

22. Graham AM, McCracken KG. Convergent evolution on the hypoxia-inducible factor (HIF) pathway genes EGLN1 and EPAS1 in high-altitude ducks. Heredity. 2019;122(6):819–32. pmid:30631144

23. Projecto-Garcia J, Natarajan C, Moriyama H, Weber RE, Fago A, Cheviron ZA, et al. Repeated elevational transitions in hemoglobin function during the evolution of Andean hummingbirds. Proc Natl Acad Sci U S A. 2013;110(51):20669–74. pmid:24297909

24. Storz JF. Hemoglobin–oxygen affinity in high-altitude vertebrates: is there evidence for an adaptive trend? J Exp Biol. 2016;219(20):3190–203. pmid:27802149

25. Murray AJ, Horscroft JA. Mitochondrial function at extreme high altitude. J Physiol. 2016;594(5):1137–49. pmid:26033622

26. Horscroft JA, Kotwica AO, Laner V, West JA, Hennis PJ, Levett DZ, et al. Metabolic basis to Sherpa altitude adaptation. Proc Natl Acad Sci U S A. 2017;114(24):6382–7. pmid:28533386

27. Dawson NJ, Lyons SA, Henry DA, Scott GR. Effects of chronic hypoxia on diaphragm function in deer mice native to high altitude. Acta Physiol. 2018;223(1):e13030. pmid:29316265

28. Dawson NJ, Alza L, Nandal G, Scott GR, McCracken KG. Convergent changes in muscle metabolism depend on duration of high-altitude ancestry across Andean waterfowl. Elife. 2020;9:e56259. pmid:32729830

29. Mahalingam S, McClelland GB, Scott GR. Evolved changes in the intracellular distribution and physiology of muscle mitochondria in high‐altitude native deer mice. J Physiol. 2017;595(14):4785–801. pmid:28418073

30. Scott GR, Guo KH, Dawson NJ. The mitochondrial basis for adaptive variation in aerobic performance in high-altitude deer mice. Integr Comp Biol. 2018;58(3):506–18. pmid:29873740

31. Cvijović I, Good BH, Desai MM. The effect of strong purifying selection on genetic diversity. Genetics. 2018;209(4):1235–78. pmid:29844134

32. Fjeldså J. Origin, evolution, and status of the avifauna of Andean wetlands. Ornithol Monogr. 1985:85–112.

33. Caziani SM, Derlindati EJ, Tálamo A, Sureda AL, Trucco CE, Nicolossi G. Waterbird richness in Altiplano wetlands of northwestern Argentina. Waterbirds. 2001;103–17.

34. McCracken KG, Barger CP, Sorenson MD. Phylogenetic and structural analysis of the HbA (aA/bA) and HbD (aD/bA) hemoglobin genes in two high-altitude waterfowl from the Himalayas and the Andes: Bar-headed goose (Anser indicus) and Andean goose (Chloephaga melanoptera). Mol Phylogenet Evol. 2010;56(2):649–58.

35. McCracken KG, Barger CP, Bulgarella M, Johnson KP, Kuhner MK, Moore AV, et al. Signatures of high-altitude adaptation in the major hemoglobin of five species of Andean dabbling ducks. Am Nat. 2009;174(5):631–50. pmid:19788356

36. Bakkeren C, Smith E, York JM, Chua B, McCracken KG, Milsom WK. A morphometric analysis of the lungs of high‐altitude ducks and geese. J Anat. 2020;237(1):188–96. pmid:32173858

37. Laguë SL, Ivy CM, York JM, Chua BA, Alza L, Cheek R, et al. Cardiovascular responses to progressive hypoxia in ducks native to high altitude in the Andes. J Exp Biol. 2020;223(5):jeb211250.

38. York JM, Chua BA, Ivy CM, Alza L, Cheek R, Scott GR, et al. Respiratory mechanics of eleven avian species resident at high and low altitude. J Exp Biol. 2017;220(6):1079–89. pmid:28298466

39. Ivy CM, Lague SL, York JM, Chua BA, Alza L, Cheek R, et al. Control of breathing and respiratory gas exchange in high-altitude ducks native to the Andes. J Exp Biol. 2019;222(7):jeb198622.

40. Lague SL, Chua B, Alza L, Scott GR, Frappell PB, Zhong Y, et al. Divergent respiratory and cardiovascular responses to hypoxia in bar-headed geese and Andean birds. J Exp Biol. 2017;220(22):4186–94. pmid:29141880

41. Milsom W, Scott G, Frappell P, McKracken K. Different strategies for convective O2 transport in high altitude birds: a graphical analysis. Comp Biochem Physiol, Part A Mol Integr Physiol. 2020:110871. pmid:33321176

42. McCracken KG, Bulgarella M, Johnson KP, Kuhner MK, Trucco J, Valqui TH, et al. Gene flow in the face of countervailing selection: adaptation to high-altitude hypoxia in the bA hemoglobin subunit of yellow-billed pintails in the Andes. Mol Biol Evol. 2009;26(4):815–27.

43. Wilson RE, Peters JL, McCracken KG. Genetic and phenotypic divergence between low- and high- altitude populations of two recently diverged cinnamon teal subspecies. Evolution. 2013;67(1):170–84. pmid:23289570

44. Graham AM, Lavretsky P, Muñoz-Fuentes V, Green AJ, Wilson RE, McCracken KG. Migration-selection balance drives genetic differentiation in genes associated with high-altitude function in the speckled teal (Anas flavirostris) in the Andes. Genome Biol Evol. 2018;10(1):14–32. pmid:29211852

45. Gnirke A, Melnikov A, Maguire J, Rogov P, LeProust EM, Brockman W, et al. Solution hybrid selection with ultra-long oligonucleotides for massively parallel targeted sequencing. Nat Biotechnol. 2009;27(2):182–9. pmid:19182786

46. Gordon A, Hannon G. Fastx-toolkit. Computer program distributed by the author, website http://hannonlabcshledu/fastx_toolkit/indexhtml [accessed 2014–2015]. 2010.

47. Hahn C, Bachmann L, Chevreux B. Reconstructing mitochondrial genomes directly from genomic next-generation sequencing reads—a baiting and iterative mapping approach. Nucleic Acids Res. 2013;41(13):e129–e. pmid:23661685

48. Katoh K, Misawa K, Kuma Ki, Miyata T. MAFFT: a novel method for rapid multiple sequence alignment based on fast Fourier transform. Nucleic Acids Res. 2002;30(14):3059–66. pmid:12136088

49. Kearse M, Moir R, Wilson A, Stones-Havas S, Cheung M, Sturrock S, et al. Geneious Basic: an integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics. 2012;28(12):1647–9. pmid:22543367

50. Weir BS, Cockerham CC. Estimating F-statistics for the analysis of population structure. evolution. 1984:1358–70. pmid:28563791

51. Excoffier L, Lischer HE. Arlequin suite v. 3.5: a new series of programs to perform population genetics analyses under Linux and Windows. Mol Ecol Resour. 2010;10(3):564–7.

52. Kumar S, Stecher G, Tamura K. MEGA7: Molecular Evolutionary Genetics Analysis v. 7.0 for bigger datasets. Mol Biol Evol. 2016;33(7)1870–1874.

53. Tajima F. Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics. 1989;123(3):585–95. pmid:2513255

54. Nei M, Kumar S. Molecular evolution and phylogenetics: Oxford university press; 2000.

55. Nei M, Gojobori T. Simple methods for estimating the numbers of synonymous and nonsynonymous nucleotide substitutions. Mol Biol Evol. 1986;3(5):418–26. pmid:3444411

56. Pond SL, Frost SD, Muse SV. HyPhy: hypothesis testing using phylogenies. Bioinformatics. 2005;21(5):676–9. pmid:15509596

57. Pond SLK, Frost SD. Not so different after all: a comparison of methods for detecting amino acid sites under selection. Mol Biol Evol. 2005;22(5):1208–22. pmid:15703242

58. Murrell B, Wertheim JO, Moola S, Weighill T, Scheffler K, Pond SLK. Detecting individual sites subject to episodic diversifying selection. PLoS Genet. 2012;8(7):e1002764. pmid:22807683

59. Murrell B, Moola S, Mabona A, Weighill T, Sheward D, Kosakovsky Pond SL, et al. FUBAR: a fast, unconstrained bayesian approximation for inferring selection. Mol Biol Evol. 2013;30(5):1196–205. pmid:23420840

60. Drummond AJ, Rambaut A. BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol Biol. 2007;7(1):214. pmid:17996036

61. Rambaut A, Drummond AJ, Xie D, Baele G, Suchard MA. Posterior summarization in Bayesian phylogenetics using Tracer 1.7. Syst Biol. 2018;67(5):901. pmid:29718447

62. Tamura K, Nei M. Estimation of the number of nucleotide substitutions in the control region of mitochondrial DNA in humans and chimpanzees. Mol Biol Evol. 1993;10(3):512–26. pmid:8336541

63. Liu G, Zhou L, Zhang L, Luo Z, Xu W. The complete mitochondrial genome of Bean goose (Anser fabalis) and implications for Anseriformes taxonomy. PLoS One. 2013;8(5):e63334. pmid:23717412

64. Yan L, Zhang C, Pan T, Zhou W, Hu C, Chang Q, et al. Mitochondrial genome of the Anas acuta (Anatidae: Anas). Mitochondrial DNA. 2015;26(2):297–8. pmid:24047167

65. Zhou W, Zhang C, Pan T, Yan L, Hu C, Xue C, et al. The complete mitochondrial genome of Anas poecilorhyncha (Anatidae: Anas). Mitochondrial DNA. 2015;26(2):265–6. pmid:24021013

66. Pan T, Ren L, Wang H, Chen J, Zhang B. Mitochondrial genome of the Anas falcata (Anatidae: Anas). Mitochondrial DNA. 2014;25(2):111–2. pmid:23586338

67. Hu C, Chang Q, Zhou W, Yan L, Pan T, Xue C, et al. Mitochondrial genome of the Anas crecca (Anatidae: Anas). Mitochondrial DNA. 2015;26(4):625–6. pmid:24438263

68. Bazinet AL, Zwickl DJ, Cummings MP. A gateway for phylogenetic analysis powered by grid computing featuring GARLI 2.0. Syst Biol. 2014;63(5):812–8. pmid:24789072

69. Zwickl D. GARLI: genetic algorithm for rapid likelihood inference. See http://wwwbioutexasedu/faculty/antisense/garli/Garlihtml. 2006.

70. Ge R-L, Simonson TS, Cooksey RC, Tanna U, Qin G, Huff CD, et al. Metabolic insight into mechanisms of high-altitude adaptation in Tibetans. Mol Genet Metab. 2012;106(2):244–7. pmid:22503288

71. Zhao D, Zhang Z, Cease A, Harrison J, Kang L. Efficient utilization of aerobic metabolism helps Tibetan locusts conquer hypoxia. BMC Genomics. 2013;14(1):1. pmid:24047108

72. Mahalingam S, Cheviron ZA, Storz JF, McClelland GB, Scott GR. Chronic cold exposure induces mitochondrial plasticity in deer mice native to high altitudes. J Physiol. 2020;598(23):5411–26. pmid:32886797

73. Wilson RE, Eaton M, McCracken KG, y Sudamérica N. color divergence among cinnamon Teal (Anas cyanoptera) subspecies from North America and South America. Ornitol Neotrop. 2008;19:307–14.

74. Wilson RE, Valqui TH, McCracken KG. Ecogeographic variation in Cinnamon Teal (Anas cyanoptera) along elevational and latitudinal gradients. Ornithol Monogr. 2010;67(1):141–61.

75. McCracken KG, Barger CP, Bulgarella M, Johnson KP, Sonsthagen SA, Trucco J, et al. Parallel evolution in the major haemoglobin genes of eight species of Andean waterfowl. Mol Ecol. 2009;18(19):3992–4005. pmid:19754505

76. Graham AM, Peters JL, Wilson RE, Muñoz-Fuentes V, Green AJ, Dorfsman DA, et al. Adaptive introgression of the beta-globin cluster in two Andean waterfowl. Heredity. 2021:1–17. pmid:33903741

77. Wolff JN, Ladoukakis ED, Enríquez JA, Dowling DK. Mitonuclear interactions: evolutionary consequences over multiple biological scales. Philos. Trans. R. Soc. Lond., B, Biol. Sci.. 2014;369(1646):20130443. pmid:24864313

78. Dowling DK, Friberg U, Lindell J. Evolutionary implications of non-neutral mitochondrial genetic variation. Trends Ecol Evol. 2008;23(10):546–54. pmid:18722688

79. Rand DM, Haney RA, Fry AJ. Cytonuclear coevolution: the genomics of cooperation. Trends Ecol Evol. 2004;19(12):645–53. pmid:16701327

80. Blier PU, Dufresne F, Burton RS. Natural selection and the evolution of mtDNA-encoded peptides: evidence for intergenomic co-adaptation. Trends Genet. 2001;17(7):400–6. pmid:11418221

81. Ellison CK, Burton RS. Disruption of mitochondrial function in interpopulation hybrids of Tigriopus californicus. Evolution. 2006;60(7):1382–91. pmid:16929655

82. Meiklejohn CD, Holmbeck MA, Siddiq MA, Abt DN, Rand DM, Montooth KL. An incompatibility between a mitochondrial tRNA and its nuclear-encoded tRNA synthetase compromises development and fitness in Drosophila. PLoS Genet. 2013;9(1):e1003238. pmid:23382693

83. Niehuis O, Judson AK, Gadau J. Cytonuclear genic incompatibilities cause increased mortality in male F2 hybrids of Nasonia giraulti and N. vitripennis. Genetics. 2008;178(1):413–26. pmid:18202384

84. Bolnick DI, Turelli M, Lopez-Fernández H, Wainwright PC, Near TJ. Accelerated mitochondrial evolution and “Darwin’s corollary”: asymmetric viability of reciprocal F1 hybrids in Centrarchid fishes. Genetics. 2008;178(2):1037–48. pmid:18245356

85. Mitchell P. Coupling of phosphorylation to electron and hydrogen transfer by a chemi-osmotic type of mechanism. Nature. 1961;191(4784):144–8. pmid:13771349

86. Havird JC, Weaver RJ, Milani L, Ghiselli F, Greenway R, Ramsey AJ, et al. Beyond the powerhouse: integrating mitonuclear evolution, physiology, and theory in comparative biology. Integr Comp Biol. 2019;59(4):856–63. pmid:31504533

87. Shen Y-Y, Shi P, Sun Y-B, Zhang Y-P. Relaxation of selective constraints on avian mitochondrial DNA following the degeneration of flight ability. Genome Res. 2009;19(10):1760–5. pmid:19617397

88. Shen Y-Y, Liang L, Zhu Z-H, Zhou W-P, Irwin DM, Zhang Y-P. Adaptive evolution of energy metabolism genes and the origin of flight in bats. Proc Natl Acad Sci U S A. 2010;107(19):8666–71. pmid:20421465

89. Mitterboeck TF, Adamowicz SJ. Flight loss linked to faster molecular evolution in insects. Proc. R. Soc. B: Biol. Sci. 2013;280(1767):20131128. pmid:23884090

90. Divakaruni AS, Brand MD. The regulation and physiology of mitochondrial proton leak. Physiol. 2011;26(3):192–205. pmid:21670165

91. Demine S, Renard P, Arnould T. Mitochondrial uncoupling: a key controller of biological processes in physiology and diseases. Cells. 2019;8(8):795. pmid:31366145

92. Maina JN, McCracken KG, Chua B, York JM, Milsom WK. Morphological and morphometric specializations of the lung of the Andean goose, Chloephaga melanoptera: a lifelong high-altitude resident. PLoS One. 2017;12(3):e0174395. pmid:28339478

93. Ridgely R, Allnutt T, Brooks T, McNicol D, Mehlman D, Young B, et al. Digital distribution maps of the birds of the Western Hemisphere. 2003.