Gas exchange, oxygen transport and metabolism in high-altitude waterfowl

High-altitude life poses physiological challenges to all animals due to decreased environmental oxygen (O2) availability (hypoxia) and cold. Supporting high metabolic rates and body temperatures with limited O2 is challenging. Many birds, however, thrive at high altitudes. The O2-transport cascade describes the pathway involved in moving O2 from the environment to the tissues encompassing: (i) ventilation, (ii) pulmonary O2 diffusion, (iii) circulation, (iv) tissue O2 diffusion, and (v) mitochondrial O2 use for ATP production. Shared avian traits such as rigid lungs with cross-current gas exchange and unidirectional airflow aid in O2 acquisition and transport in all birds. Many high-altitude birds, however, have evolved enhancements to some or all steps in the cascade. In this review, we summarize the current literature on gas exchange and O2 transport in high-altitude birds, providing an overview of the O2-transport cascade that principally draws on the literature from high-altitude waterfowl, the most well-studied group of high-altitude birds. We close by discussing two important avenues for future research: distinguishing between the influences of plasticity and evolution and investigating whether the morphological and physiological differences discussed contribute to enhanced locomotor or thermogenic performance, a potential critical link to fitness.This article is part of the theme issue 'The biology of the avian respiratory system'.

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

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 (Ne) 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.

High-altitude adaptation in vertebrates as revealed by mitochondrial genome analyses

The high-altitude environment may drive vertebrate evolution in a certain way, and vertebrates living in different altitude environments might have different energy requirements. We hypothesized that the high-altitude environment might impose different influences on vertebrate mitochondrial genomes (mtDNA). We used selection pressure analyses and PIC (phylogenetic independent contrasts) analysis to detect the evolutionary rate of vertebrate mtDNA protein-coding genes (PCGs) from different altitudes. The results showed that the ratio of nonsynonymous/synonymous substitutions (dN/dS) in the mtDNA PCGs was significantly higher in high-altitude vertebrates than in low-altitude vertebrates. The seven rapidly evolving genes were shared by the high-altitude vertebrates, and only one positive selection gene (ND5 gene) was detected in the high-altitude vertebrates. Our results suggest the mtDNA evolutionary rate in high-altitude vertebrates was higher than in low-altitude vertebrates as their evolution requires more energy in a high-altitude environment. Our study demonstrates the high-altitude environment (low atmospheric O2 levels) drives vertebrate evolution in mtDNA PCGs.

Adaptive introgression of the beta-globin cluster in two Andean waterfowl

Introgression of beneficial alleles has emerged as an important avenue for genetic adaptation in both plant and animal populations. In vertebrates, adaptation to hypoxic high-altitude environments involves the coordination of multiple molecular and cellular mechanisms, including selection on the hypoxia-inducible factor (HIF) pathway and the blood-O2 transport protein hemoglobin (Hb). In two Andean duck species, a striking DNA sequence similarity reflecting identity by descent is present across the ~20 kb β-globin cluster including both embryonic (HBE) and adult (HBB) paralogs, though it was yet untested whether this is due to independent parallel evolution or adaptive introgression. In this study, we find that identical amino acid substitutions in the β-globin cluster that increase Hb-O2 affinity have likely resulted from historical interbreeding between high-altitude populations of two different distantly-related species. We examined the direction of introgression and discovered that the species with a deeper mtDNA divergence that colonized high altitude earlier in history (Anas flavirostris) transferred adaptive genetic variation to the species with a shallower divergence (A. georgica) that likely colonized high altitude more recently possibly following a range shift into a novel environment. As a consequence, the species that received these β-globin variants through hybridization might have adapted to hypoxic conditions in the high-altitude environment more quickly through acquiring beneficial alleles from the standing, hybrid-origin variation, leading to faster evolution.

Interspecific Gene Flow and the Evolution of Specialization in Black and White Rhinoceros

Africa's black (Diceros bicornis) and white (Ceratotherium simum) rhinoceros are closely related sister-taxa that evolved highly divergent obligate browsing and grazing feeding strategies. Although their precursor species Diceros praecox and Ceratotherium mauritanicum appear in the fossil record ∼5.2 Ma, by 4 Ma both were still mixed feeders, and were even spatiotemporally sympatric at several Pliocene sites in what is today Africa's Rift Valley. Here, we ask whether or not D. praecox and C. mauritanicum were reproductively isolated when they came into Pliocene secondary contact. We sequenced and de novo assembled the first annotated black rhinoceros reference genome and compared it with available genomes of other black and white rhinoceros. We show that ancestral gene flow between D. praecox and C. mauritanicum ceased sometime between 3.3 and 4.1 Ma, despite conventional methods for the detection of gene flow from whole genome data returning false positive signatures of recent interspecific migration due to incomplete lineage sorting. We propose that ongoing Pliocene genetic exchange, for up to 2 My after initial divergence, could have potentially hindered the development of obligate feeding strategies until both species were fully reproductively isolated, but that the more severe and shifting paleoclimate of the early Pleistocene was likely the ultimate driver of ecological specialization in African rhinoceros.

Convergent changes in muscle metabolism depend on duration of high-altitude ancestry across Andean waterfowl

High-altitude environments require that animals meet the metabolic O2 demands for locomotion and thermogenesis in O2-thin air, but the degree to which convergent metabolic changes have arisen across independent high-altitude lineages or the speed at which such changes arise is unclear. We examined seven high-altitude waterfowl that have inhabited the Andes (3812-4806 m elevation) over varying evolutionary time scales, to elucidate changes in biochemical pathways of energy metabolism in flight muscle relative to low-altitude sister taxa. Convergent changes across high-altitude taxa included increased hydroxyacyl-coA dehydrogenase and succinate dehydrogenase activities, decreased lactate dehydrogenase, pyruvate kinase, creatine kinase, and cytochrome c oxidase activities, and increased myoglobin content. ATP synthase activity increased in only the longest established high-altitude taxa, whereas hexokinase activity increased in only newly established taxa. Therefore, changes in pathways of lipid oxidation, glycolysis, and mitochondrial oxidative phosphorylation are common strategies to cope with high-altitude hypoxia, but some changes require longer evolutionary time to arise.

Comparative morphometric analysis of lungs of the semifossorial giant pouched rat (Cricetomys gambianus) and the subterranean Nigerian mole rat (Cryptomys foxi)

Lungs of the rodent species, the African giant pouched rat (Cricetomys gambianus) and the Nigerian mole rat (Cryptomys foxi) were investigated. Significant morphometric differences exist between the two species. The volume of the lung per unit body mass was 2.7 times larger; the respiratory surface area 3.4 times greater; the volume of the pulmonary capillary blood 2 times more; the harmonic mean thickness of the blood-gas (tissue) barrier (τht) ~29% thinner and; the total pulmonary morphometric diffusing capacity (DLo2) for O2 2.3 times more in C. foxi. C. gambianus occupies open burrows that are ventilated with air while C. foxi lives in closed burrows. The less morphometrically specialized lungs of C. gambianus may be attributed to its much larger body mass (~6 times more) and possibly lower metabolic rate and its semifossorial life whereas the 'superior' lungs of C. foxi may largely be ascribed to the subterranean hypoxic and hypercapnic environment it occupies. Compared to other rodents species that have been investigated hitherto, the τht was mostly smaller in the lungs of the subterranean species and C. foxi has the highest mass-specific DLo2. The fossorial- and the subterranean rodents have acquired various pulmonary structural specializations that relate to habitats occupied.

Cardiovascular responses to progressive hypoxia in ducks native to high altitude in the Andes

The cardiovascular system is critical for delivering O2 to tissues. Here we examine the cardiovascular responses to progressive hypoxia in four high-altitude Andean duck species compared to four related low-altitude populations in North America, tested at their native altitude. Ducks were exposed to stepwise decreases in inspired partial pressure of O2 while we monitored heart rate, O2 consumption rate, blood O2 saturation, haematocrit (Hct), and blood haemoglobin concentration [Hb]. We calculated O2 pulse (the product of stroke volume and the arterial-venous O2 content difference), blood O2 concentration, and heart rate variability. Regardless of altitude, all eight populations maintained O2 consumption rate with minimal change in heart rate or O2 pulse, indicating that O2 consumption was maintained by either a constant arterial-venous O2 content difference (an increase in the relative O2 extracted from arterial blood) or by a combination of changes in stroke volume and the arterial-venous O2 content difference. Three high-altitude taxa (yellow-billed pintails, cinnamon teal, and speckled teal) had higher Hct and [Hb], increasing the O2 content of arterial blood, and potentially providing a greater reserve for enhancing O2 delivery during hypoxia. Hct and [Hb] between low- and high-altitude populations of ruddy duck were similar, representing a potential adaptation to diving life. Heart rate variability was generally lower in high-altitude ducks, concurrent with similar or lower heart rates than low-altitude ducks, suggesting a reduction in vagal and sympathetic tone. These unique features of the Andean ducks differ from previous observations in both Andean geese and bar-headed geese, neither of which exhibit significant elevations in Hct or [Hb] compared to their low-altitude relatives, revealing yet another avian strategy for coping with high altitude.

Rapid phenotypic evolution with shallow genomic differentiation during early stages of high elevation adaptation in Eurasian Tree Sparrows

Known as the ‘third polar region’, the Qinghai-Tibet Plateau represents one of the harshest highland environments in the world and yet a number of organisms thrive there. Previous studies of birds, animals and humans have focused on well-differentiated populations in later stages of phenotypic divergence. The adaptive processes during the initial phase of highland adaptation remain poorly understood. We studied a human commensal, the Eurasian Tree Sparrow, which has followed human agriculture to the Qinghai-Tibet Plateau. Despite strong phenotypic differentiation at multiple levels, in particular in muscle-related phenotypes, highland and lowland populations show shallow genomic divergence and the colonization event occurred within the past few thousand years. In a one-month acclimation experiment investigating phenotypic plasticity, we exposed adult lowland tree sparrows to a hypoxic environment and did not observe muscle changes. Through population genetic analyses, we identified a signature of polygenic adaptation, whereby shifts in allele frequencies are spread across multiple loci, many of which are associated with muscle-related processes. Our results reveal a case of positive selection in which polygenic adaptation appears to drive rapid phenotypic evolution, shedding light on early stages of adaptive evolution to a novel environment.

Control of breathing and respiratory gas exchange in high-altitude ducks native to the Andes

We examined the control of breathing and respiratory gas exchange in six species of high-altitude duck that independently colonized the high Andes. We compared ducks from high-altitude populations in Peru (Lake Titicaca at ∼3800 m above sea level; Chancay River at ∼3000-4100 m) with closely related populations or species from low altitude. Hypoxic ventilatory responses were measured shortly after capture at the native altitude. In general, ducks responded to acute hypoxia with robust increases in total ventilation and pulmonary O₂ extraction. O₂ consumption rates were maintained or increased slightly in acute hypoxia, despite ∼1-2°C reductions in body temperature in most species. Two high-altitude taxa - yellow-billed pintail and torrent duck - exhibited higher total ventilation than their low-altitude counterparts, and yellow-billed pintail exhibited greater increases in pulmonary O₂ extraction in severe hypoxia. In contrast, three other high-altitude taxa - Andean ruddy duck, Andean cinnamon teal and speckled teal - had similar or slightly reduced total ventilation and pulmonary O₂ extraction compared with low-altitude relatives. Arterial O₂ saturation (S _aO2 ) was elevated in yellow-billed pintails at moderate levels of hypoxia, but there were no differences in S _aO2  in other high-altitude taxa compared with their close relatives. This finding suggests that improvements in S _aO2  in hypoxia can require increases in both breathing and haemoglobin-O₂ affinity, because the yellow-billed pintail was the only high-altitude duck with concurrent increases in both traits compared with its low-altitude relative. Overall, our results suggest that distinct physiological strategies for coping with hypoxia can exist across different high-altitude lineages, even among those inhabiting very similar high-altitude habitats.

Convergent evolution on the hypoxia-inducible factor (HIF) pathway genes EGLN1 and EPAS1 in high-altitude ducks

During periods of reduced O₂ supply, the most profound changes in gene expression are mediated by hypoxia-inducible factor (HIF) transcription factors that play a key role in cellular responses to low-O₂ tension. Using target-enrichment sequencing, we tested whether variation in 26 genes in the HIF signaling pathway was associated with high altitude and therefore corresponding O₂ availability in three duck species that colonized the Andes from ancestral low-altitude habitats in South America. We found strong support for convergent evolution in the case of two of the three duck species with the same genes (EGLN1, EPAS1), and even the same exons (exon 12, EPAS1), exhibiting extreme outliers with a high probability of directional selection in the high-altitude populations. These results mirror patterns of adaptation seen in human populations, which showed mutations in EPAS1, and transcriptional regulation differences in EGLN1, causing changes in downstream target transactivation, associated with a blunted hypoxic response.

Genetic variation in PTPN1 contributes to metabolic adaptation to high-altitude hypoxia in Tibetan migratory locusts

Animal and human highlanders have evolved distinct traits to enhance tissue oxygen delivery and utilization. Unlike vertebrates, insects use their tracheal system for efficient oxygen delivery. However, the genetic basis of insect adaptation to high-altitude hypoxia remains unexplored. Here, we report a potential mechanism of metabolic adaptation of migratory locusts in the Tibetan Plateau, through whole-genome resequencing and functional investigation. A genome-wide scan revealed that the positively selected genes in Tibetan locusts are predominantly involved in carbon and energy metabolism. We observed a notable signal of natural selection in the gene PTPN1, which encodes PTP1B, an inhibitor of insulin signaling pathway. We show that a PTPN1 coding mutation regulates the metabolism of Tibetan locusts by mediating insulin signaling activity in response to hypoxia. Overall, our findings provide evidence for the high-altitude hypoxia adaptation of insects at the genomic level and explore a potential regulatory mechanism underlying the evolved metabolic homeostasis.

Vertebrate adaptation to high-altitude life has been extensively investigated, while invertebrates are less well-studied. Here, the authors find signals of adaptive evolution in genomes of migratory locusts from the Tibetan Plateau, and implicate a PTPN1 coding mutation in their hypoxia response.

Elevational divergence in the great tit complex revealed by major hemoglobin genes

Gene flow and demographic history can play important roles in the adaptive genetic differentiation of species, which is rarely understood in the high-altitude adaptive evolution of birds. To elucidate genetic divergence of populations in the great tit complex (Parus major, P. minor and P. cinereus) at different elevations, we compared the genetic structure and gene flow in hemoglobin genes with neutral loci. Our results revealed the elevationally divergent structure of α^A-globin gene, distinctive from that of the β^A-globin gene and neutral loci. We further investigated gene flow patterns among the populations in the central-northern (> 1,000 m a.s.l.), south-eastern (< 1,000 m a.s.l.) regions and the Southwest Mountains (> 2,000 m a.s.l.) in China. The high-altitude (> 1,000 m a.s.l.) diverged α^A-globin genetic structure coincided with higher α^A-globin gene flow between highland populations, in contrast to restricted neutral gene flow concordant with the phylogeny. The higher α^A-globin gene flow suggests the possibility of adaptive evolution during population divergence, contrary to the lower α^A-globin gene flow homogenized by neutral loci during population expansion. In concordance with patterns of historical gene flow, genotypic and allelic profiles provide distinctive patterns of fixation in different high-altitude populations. The fixation of alleles at contrasting elevations may primarily due to highland standing variants α^A49Asn/72Asn/108Ala originating from the south-western population. Our findings demonstrate a pattern of genetic divergence with gene flow in major hemoglobin genes depending on population demographic history.

Neutral and functionally important genes shed light on phylogeography and the history of high-altitude colonization in a widespread New World duck

Phylogeographic studies often infer historical demographic processes underlying species distributions based on patterns of neutral genetic variation, but spatial variation in functionally important genes can provide additional insights about biogeographic history allowing for inferences about the potential role of adaptation in geographic range evolution. Integrating data from neutral markers and genes involved in oxygen (O2)-transport physiology, we test historical hypotheses about colonization and gene flow across low- and high-altitude regions in the Ruddy Duck (Oxyura jamaicensis), a widely distributed species in the New World. Using multilocus analyses that for the first time include populations from the Colombian Andes, we also examined the hypothesis that Ruddy Duck populations from northern South America are of hybrid origin. We found that neutral and functional genes appear to have moved into the Colombian Andes from both North America and southern South America, and that high-altitude Colombian populations do not exhibit evidence of adaptation to hypoxia in hemoglobin genes. Therefore, the biogeographic history of Ruddy Ducks is likely more complex than previously inferred. Our new data raise questions about the hypothesis that adaptation via natural selection to high-altitude conditions through amino acid replacements in the hemoglobin protein allowed Ruddy Ducks to disperse south along the high Andes into southern South America. The existence of shared genetic variation with populations from both North America and southern South America as well as private alleles suggests that the Colombian population of Ruddy Ducks may be of old hybrid origin. This study illustrates the breadth of inferences one can make by combining data from nuclear and functionally important loci in phylogeography, and underscores the importance of complete range-wide sampling to study species history in complex landscapes.

Migration-Selection Balance Drives Genetic Differentiation in Genes Associated with High-Altitude Function in the Speckled Teal (Anas flavirostris) in the Andes

Local adaptation frequently occurs across populations as a result of migration-selection balance between divergent selective pressures and gene flow associated with life in heterogeneous landscapes. Studying the effects of selection and gene flow on the adaptation process can be achieved in systems that have recently colonized extreme environments. This study utilizes an endemic South American duck species, the speckled teal (Anas flavirostris), which has both high- and low-altitude populations. High-altitude speckled teal (A. f. oxyptera) are locally adapted to the Andean environment and mostly allopatric from low-altitude birds (A. f. flavirostris); however, there is occasional gene flow across altitudinal gradients. In this study, we used next-generation sequencing to explore genetic patterns associated with high-altitude adaptation in speckled teal populations, as well as the extent to which the balance between selection and migration have affected genetic architecture. We identified a set of loci with allele frequencies strongly correlated with altitude, including those involved in the insulin-like signaling pathway, bone morphogenesis, oxidative phosphorylation, responders to hypoxia-induced DNA damage, and feedback loops to the hypoxia-inducible factor pathway. These same outlier loci were found to have depressed gene flow estimates, as well as being highly concentrated on the Z-chromosome. Our results suggest a multifactorial response to life at high altitudes through an array of interconnected pathways that are likely under positive selection and whose genetic components seem to be providing an effective genomic barrier to interbreeding, potentially functioning as an avenue for population divergence and speciation.

Respiratory mechanics of eleven avian species resident at high and low altitude

The metabolic cost of breathing at rest has never been successfully measured in birds, but has been hypothesized to be higher than in mammals of a similar size because of the rocking motion of the avian sternum being encumbered by the pectoral flight muscles. To measure the cost and work of breathing, and to investigate whether species resident at high altitude exhibit morphological or mechanical changes that alter the work of breathing, we studied 11 species of waterfowl: five from high altitudes (&gt;3000 m) in Perú, and six from low altitudes in Oregon, USA. Birds were anesthetized and mechanically ventilated in sternal recumbency with known tidal volumes and breathing frequencies. The work done by the ventilator was measured, and these values were applied to the combinations of tidal volumes and breathing frequencies used by the birds to breathe at rest. We found the respiratory system of high-altitude species to be of a similar size, but consistently more compliant than that of low-altitude sister taxa, although this did not translate to a significantly reduced work of breathing. The metabolic cost of breathing was estimated to be between 1 and 3% of basal metabolic rate, as low or lower than estimates for other groups of tetrapods.

Absence of population structure across elevational gradients despite large phenotypic variation in mountain chickadees (Poecile gambeli)

Montane habitats are characterized by predictably rapid heterogeneity along elevational gradients and are useful for investigating the consequences of environmental heterogeneity for local adaptation and population genetic structure. Food-caching mountain chickadees inhabit a continuous elevation gradient in the Sierra Nevada, and birds living at harsher, high elevations have better spatial memory ability and exhibit differences in male song structure and female mate preference compared to birds inhabiting milder, low elevations. While high elevation birds breed, on average, two weeks later than low elevation birds, the extent of gene flow between elevations is unknown. Despite phenotypic variation and indirect evidence for local adaptation, population genetic analyses based on 18 073 single nucleotide polymorphisms across three transects of high and low elevation populations provided no evidence for genetic differentiation. Analyses based on individual genotypes revealed no patterns of clustering, pairwise estimates of genetic differentiation (F_ST, Nei's D) were very low, and AMOVA revealed no evidence for genetic variation structured by transect or by low and high elevation sites within transects. In addition, we found no consistent evidence for strong parallel allele frequency divergence between low and high elevation sites within the three transects. Large elevation-related phenotypic variation may be maintained by strong selection despite gene flow and future work should focus on the mechanisms underlying such variation.

Distinct patterns of natural selection in Na+/H+ antiporter genes in Populus euphratica and Populus pruinosa

Salt tolerance genes constitute an important class of loci in plant genomes. Little is known about the extent to which natural selection in saline environments has acted upon these loci, and what types of nucleotide diversity such selection has given rise to. Here, we surveyed genetic diversity in three types of Na⁺/H⁺ antiporter gene (SOS, NhaD, and NHX, belonging to the cation/proton antiporter 1 family), which have well‐characterized essential roles in plant salt tolerance. Ten Na⁺/H⁺ antiporter genes and 16 neutral loci randomly selected as controls were sequenced from 17 accessions of two closely related members of the genus Populus, Populus euphratica and Populus pruinosa, section Turanga, which are native to northwest China. The results show that salt tolerance genes are common targets of natural selection in P. euphratica and P. pruinosa. Moreover, the patterns of nucleotide variation across the three types of Na⁺/H⁺ antiporter gene are distinctly different in these two closely related Populus species, and gene flow from P. pruinosa to P. euphratica is highly restricted. Our results suggest that natural selection played an important role in shaping the current distinct patterns of Na⁺/H⁺ antiporter genes, resulting in adaptive evolution in P. euphratica and P. pruinosa.

Mitochondrial physiology in the skeletal and cardiac muscles is altered in torrent ducks, Merganetta armata, from high altitudes in the Andes

Torrent ducks inhabit fast-flowing rivers in the Andes from sea level to altitudes up to 4500 m. We examined the mitochondrial physiology that facilitates performance over this altitudinal cline by comparing the respiratory capacities of permeabilized fibers, the activities of 16 key metabolic enzymes and the myoglobin content in muscles between high- and low-altitude populations of this species. Mitochondrial respiratory capacities (assessed using substrates of mitochondrial complexes I, II and/or IV) were higher in highland ducks in the gastrocnemius muscle - the primary muscle used to support swimming and diving - but were similar between populations in the pectoralis muscle and the left ventricle. The heightened respiratory capacity in the gastrocnemius of highland ducks was associated with elevated activities of cytochrome oxidase, phosphofructokinase, pyruvate kinase and malate dehydrogenase (MDH). Although respiratory capacities were similar between populations in the other muscles, highland ducks had elevated activities of ATP synthase, lactate dehydrogenase, MDH, hydroxyacyl CoA dehydrogenase and creatine kinase in the left ventricle, and elevated MDH activity and myoglobin content in the pectoralis. Thus, although there was a significant increase in the oxidative capacity of the gastrocnemius in highland ducks, which correlates with improved performance at high altitudes, the variation in metabolic enzyme activities in other muscles not correlated to respiratory capacity, such as the consistent upregulation of MDH activity, may serve other functions that contribute to success at high altitudes.

Mountain chickadees from different elevations sing different songs: acoustic adaptation, temporal drift or signal of local adaptation?

Song in songbirds is widely thought to function in mate choice and male-male competition. Song is also phenotypically plastic and typically learned from local adults; therefore, it varies across geographical space and can serve as a cue for an individual's location of origin, with females commonly preferring males from their respective location. Geographical variation in song dialect may reflect acoustic adaptation to different environments and/or serve as a signal of local adaptation. In montane environments, environmental differences can occur over an elevation gradient, favouring local adaptations across small spatial scales. We tested whether food caching mountain chickadees, known to exhibit elevation-related differences in food caching intensity, spatial memory and the hippocampus, also sing different dialects despite continuous distribution and close proximity. Male songs were collected from high and low elevations at two different mountains (separated by 35 km) to test whether song differs between elevations and/or between adjacent populations at each mountain. Song structure varied significantly between high and low elevation adjacent populations from the same mountain and between populations from different mountains at the same elevations, despite a continuous distribution across each mountain slope. These results suggest that elevation-related differences in song structure in chickadees might serve as a signal for local adaptation.

Mito-nuclear discord in six congeneric lineages of Holarctic ducks (genus Anas)

Many species have Holarctic distributions that extend across Europe, Asia and North America. Most genetics research on these species has examined only mitochondrial (mt) DNA, which has revealed wide variance in divergence between Old World (OW) and New World (NW) populations, ranging from shallow, unstructured genealogies to deeply divergent lineages. In this study, we sequenced 20 nuclear introns to test for concordant patterns of OW-NW differentiation between mtDNA and nuclear (nu) DNA for six lineages of Holarctic ducks (genus Anas). Genetic differentiation for both marker types varied widely among these lineages (idiosyncratic population histories), but mtDNA and nuDNA divergence within lineages was not significantly correlated. Moreover, compared with the association between mtDNA and nuDNA divergence observed among different species, OW-NW nuDNA differentiation was generally lower than mtDNA divergence, at least for lineages with deeply divergent mtDNA. Furthermore, coalescent estimates indicated significantly higher rates of gene flow for nuDNA than mtDNA for four of the six lineages. Thus, Holarctic ducks show prominent mito-nuclear discord between OW and NW populations, and we reject differences in sorting rates as the sole cause of the within-species discord. Male-mediated intercontinental gene flow is likely a leading contributor to this discord, although selection could also cause increased mtDNA divergence relative to weak nuDNA differentiation. The population genetics of these ducks contribute to growing evidence that mtDNA can be an unreliable indicator of stage of speciation and that more holistic approaches are needed for species delimitation.

Fin whale MDH-1 and MPI allozyme variation is not reflected in the corresponding DNA sequences

The appeal of genetic inference methods to assess population genetic structure and guide management efforts is grounded in the correlation between the genetic similarity and gene flow among populations. Effects of such gene flow are typically genomewide; however, some loci may appear as outliers, displaying above or below average genetic divergence relative to the genomewide level. Above average population, genetic divergence may be due to divergent selection as a result of local adaptation. Consequently, substantial efforts have been directed toward such outlying loci in order to identify traits subject to local adaptation. Here, we report the results of an investigation into the molecular basis of the substantial degree of genetic divergence previously reported at allozyme loci among North Atlantic fin whale (Balaenoptera physalus) populations. We sequenced the exons encoding for the two most divergent allozyme loci (MDH-1 and MPI) and failed to detect any nonsynonymous substitutions. Following extensive error checking and analysis of additional bioinformatic and morphological data, we hypothesize that the observed allozyme polymorphisms may reflect phenotypic plasticity at the cellular level, perhaps as a response to nutritional stress. While such plasticity is intriguing in itself, and of fundamental evolutionary interest, our key finding is that the observed allozyme variation does not appear to be a result of genetic drift, migration, or selection on the MDH-1 and MPI exons themselves, stressing the importance of interpreting allozyme data with caution. As for North Atlantic fin whale population structure, our findings support the low levels of differentiation found in previous analyses of DNA nucleotide loci.

Gene flow and hybridization between numerically imbalanced populations of two duck species on the subantarctic island of South Georgia

Hybridization is common between species of animals, particularly in waterfowl (Anatidae). One factor shown to promote hybridization is restricted mate choice, which can occur when 2 species occur in sympatry but one is rare. According to the Hubbs principle, or "desperation hypothesis," the rarer species is more likely to mate with heterospecifics. We report the second of 2 independent examples of hybridization between 2 species of ducks inhabiting island ecosystems in the Subantarctic and South Atlantic Ocean. Yellow-billed pintails (Anas georgica) and speckled teal (Anas flavirostris) are abundant in continental South America, where they are sympatric and coexist in mixed flocks. But on South Georgia, an isolated island in the Subantarctic, the pintail population of approximately 6000 pairs outnumbers a small breeding population of speckled teal 300∶1. Using 6 genetic loci (mtDNA and 5 nuclear introns) and Bayesian assignment tests coupled with coalescent analyses, we identified hybrid-origin speckled teal alleles in 2 pintails on South Georgia. While it is unclear whether introgression has also occurred into the speckled teal population, our data suggest that this hybridization was not a recent event, but occurred some time ago. We also failed to identify unequivocal evidence of introgression in a much larger sample of pintails and speckled teal from Argentina using a 3-population "Isolation-with-Migration" coalescent analysis. Combined with parallel findings of hybridization between these same 2 duck species in the Falkland Islands, where population ratios are reversed and pintails are outnumbered by speckled teal 1:10, our results provide further support for the desperation hypothesis, which predicts that scarcity in one population and abundance of another will often lead to hybridization. While the South Georgia pintail population appears to be thriving, it's possible that low density of conspecific mates and inverse density dependence (Allee effect) may be one factor limiting the reproductive output of the speckled teal population, and this situation may persist unless speckled teal increase in abundance on South Georgia.

A gene with major phenotypic effects as a target for selection vs. homogenizing gene flow

Genes with major phenotypic effects facilitate quantifying the contribution of genetic vs. plastic effects to adaptive divergence. A classical example is Ectodysplasin (Eda), the major gene controlling lateral plate phenotype in three-spined stickleback. Completely plated marine stickleback populations evolved repeatedly towards low-plated freshwater populations, representing a prime example of parallel evolution by natural selection. However, many populations remain polymorphic for lateral plate number. Possible explanations for this polymorphism include relaxation of selection, disruptive selection or a balance between divergent selection and gene flow. We investigated 15 polymorphic stickleback populations from brackish and freshwater habitats in coastal North-western Europe. At each site, we tracked changes in allele frequency at the Eda gene between subadults in fall, adults in spring and juveniles in summer. Eda genotypes were also compared for body size and reproductive investment. We observed a fitness advantage for the Eda allele for the low morph in freshwater and for the allele for the complete morph in brackish water. Despite these results, the differentiation at the Eda gene was poorly correlated with habitat characteristics. Neutral population structure was the best predictor of spatial variation in lateral plate number, suggestive of a substantial effect of gene flow. A meta-analysis revealed that the signature of selection at Eda was weak compared to similar studies in stickleback. We conclude that a balance between divergent selection and gene flow can maintain stickleback populations polymorphic for lateral plate number and that ecologically relevant genes may not always contribute much to local adaptation, even when targeted by selection.

Stepwise colonization of the Andes by ruddy ducks and the evolution of novel β-globin variants

Andean uplift played a key role in Neotropical bird diversification, yet past dispersal and genetic adaptation to high-altitude environments remain little understood. Here we use multilocus population genetics to study population history and historical demographic processes in the ruddy duck (Oxyura jamaicensis), a stiff-tailed diving duck comprising three subspecies distributed from Canada to Tierra del Fuego and inhabiting wetlands from sea level to 4500 m in the Andes. We sequenced the mitochondrial DNA, four autosomal introns and three haemoglobin genes (α(A), α(D), β(A)) and used isolation-with-migration (IM) models to study gene flow between North America and South America, and between the tropical and southern Andes. Our analyses indicated that ruddy ducks dispersed first from North America to the tropical Andes, then from the tropical Andes to the southern Andes. While no nonsynonymous substitutions were found in either α globin gene, three amino acid substitutions were observed in the β(A) globin. Based on phylogenetic reconstruction and power analysis, the first β(A) substitution, found in all Andean individuals, was acquired when ruddy ducks dispersed from low altitude in North America to high altitude in the tropical Andes, whereas the two additional substitutions occurred more recently, when ruddy ducks dispersed from high altitude in the tropical Andes to low altitude in the southern Andes. This stepwise colonization pattern accompanied by polarized β(A) globin amino acid replacements suggest that ruddy ducks first acclimatized or adapted to the Andean highlands and then again to the lowlands. In addition, ruddy ducks colonized the Andean highlands via a less common route as compared to other waterbird species that colonized the Andes northwards from the southern cone of South America.

Ecological speciation along an elevational gradient in a tropical passerine bird?

Local adaptation of populations along elevational gradients is well known, but conclusive evidence that such divergence has resulted in the origin of distinct species in parapatry remains lacking. We integrated morphological, vocal, genetic and behavioural data to test predictions pertaining to the hypothesis of parapatric ecological speciation associated with elevation in populations of a tropical montane songbird, the Grey-breasted Wood-wren (Henicorhina leucophrys: Troglodytidae), from the Sierra Nevada de Santa Marta, Colombia. We confirmed that two distinct populations exist along the elevational gradient. Phylogenetic analyses tentatively indicate that the two populations are not sister taxa, suggesting they did not differentiate from a single ancestor along the gradient, but rather resulted from separate colonization events. The populations showed marked divergence in morphometrics, vocalizations and genetic variation in mitochondrial and nuclear loci, and little to no evidence of hybridization. Individuals of both populations responded more strongly to their own local songs than to songs from another elevation. Although the two forms do not appear to have differentiated locally in parapatry, morphological and vocal divergence along the elevational gradient is consistent with adaptation, suggesting a possible link between adaptive evolution in morphology and songs and the origin of reproductive isolation via a behavioural barrier to gene flow. The adaptive value of phenotypic differences between populations requires additional study.

Gene flow networks among American Aedes aegypti populations

The mosquito Aedes aegypti, the dengue virus vector, has spread throughout the tropics in historical times. While this suggests man-mediated dispersal, estimating contemporary connectivity among populations has remained elusive. Here, we use a large mtDNA dataset and a Bayesian coalescent framework to test a set of hypotheses about gene flow among American Ae. aegypti populations. We assessed gene flow patterns at the continental and subregional (Amazon basin) scales. For the Americas, our data favor a stepping-stone model in which gene flow is higher among adjacent populations but in which, at the same time, North American and southeastern Brazilian populations are directly connected, likely via sea trade. Within Amazonia, the model with highest support suggests extensive gene flow among major cities; Manaus, located at the center of the subregional transport network, emerges as a potentially important connecting hub. Our results suggest substantial connectivity across Ae. aegypti populations in the Americas. As long-distance active dispersal has not been observed in this species, our data support man-mediated dispersal as a major determinant of the genetic structure of American Ae. aegypti populations. The inferred topology of interpopulation connectivity can inform network models of Ae. aegypti and dengue spread.

Genetic and phenotypic divergence between low- and high-altitude populations of two recently diverged cinnamon teal subspecies

Spatial variation in the environment can lead to divergent selection between populations occupying different parts of a species' range, and ultimately lead to population divergence. The colonization of new areas can thus facilitate divergence in beneficial traits, yet with little differentiation at neutral genetic markers. We investigated genetic and phenotypic patterns of divergence between low- and high-altitude populations of cinnamon teal inhabiting normoxic and hypoxic regions in the Andes and adjacent lowlands of South America. Cinnamon teal showed strong divergence in body size (PC1; P(ST) = 0.56) and exhibited significant frequency differences in a single nonsynonymous α-hemoglobin amino acid polymorphism (Asn/Ser-α9; F(ST) = 0.60) between environmental extremes, despite considerable admixture of mtDNA and intron loci (F(ST) = 0.004-0.168). Inferences of strong population segregation were further supported by the observation of few mismatched individuals in either environmental extreme. Coalescent analyses indicated that the highlands were most likely colonized from lowland regions but following divergence, gene flow has been asymmetric from the highlands into the lowlands. Multiple selection pressures associated with high-altitude habitats, including cold and hypoxia, have likely shaped morphological and genetic divergence within South American cinnamon teal populations.

Heterogeneity in genetic diversity among non-coding loci fails to fit neutral coalescent models of population history

Inferring aspects of the population histories of species using coalescent analyses of non-coding nuclear DNA has grown in popularity. These inferences, such as divergence, gene flow, and changes in population size, assume that genetic data reflect simple population histories and neutral evolutionary processes. However, violating model assumptions can result in a poor fit between empirical data and the models. We sampled 22 nuclear intron sequences from at least 19 different chromosomes (a genomic transect) to test for deviations from selective neutrality in the gadwall (Anas strepera), a Holarctic duck. Nucleotide diversity among these loci varied by nearly two orders of magnitude (from 0.0004 to 0.029), and this heterogeneity could not be explained by differences in substitution rates alone. Using two different coalescent methods to infer models of population history and then simulating neutral genetic diversity under these models, we found that the observed among-locus heterogeneity in nucleotide diversity was significantly higher than expected for these simple models. Defining more complex models of population history demonstrated that a pre-divergence bottleneck was also unlikely to explain this heterogeneity. However, both selection and interspecific hybridization could account for the heterogeneity observed among loci. Regardless of the cause of the deviation, our results illustrate that violating key assumptions of coalescent models can mislead inferences of population history.

Multilocus coalescent analysis of haemoglobin differentiation between low- and high-altitude populations of crested ducks (Lophonetta specularioides)

Hypoxia is a key factor determining survival, and haemoglobins are targets of selection in species native to high-altitude regions. We studied population genetic structure and evaluated evidence for local adaptation in the crested duck (Lophonetta specularioides). Differentiation, gene flow and time since divergence between highland and lowland populations were assessed for three haemoglobin genes (α(A) , α(D) , β(A) ) and compared to seven reference loci (six autosomal introns and mtDNA). Four derived amino acid replacements were found in the globin genes that had elevated Φ(ST) values between the Andean highlands and Patagonian lowlands. A single β(A) -globin polymorphism at a site known to influence O(2) affinity was fixed for different alleles in the two populations, whereas three α(A) - and α(D) -globin polymorphisms exhibited high heterozygosity in the highlands but not in the lowlands. Coalescent analyses supported restricted gene flow for haemoglobin alleles and mitochondrial DNA but nonzero gene flow for the introns. Simulating genetic data under a drift-migration model of selective neutrality, the β(A) -globin fell outside the 95% confidence limit of simulated data, suggesting that directional selection is maintaining different variants in the contrasting elevational environments, thereby restricting migration of β(A) -globin alleles. The α(A) - and α(D) -globins, by contrast, did not differ from the simulated values, suggesting that variants in these genes are either selectively neutral, or that the effects of selection could not be differentiated from background levels of population structure and linkage disequilibrium. This study illustrates the combined effects of selection and population history on inferring levels of population divergence for a species distributed across an altitudinal gradient in which selection for hypoxia resistance has likely played an important role.

Elevated performance: the unique physiology of birds that fly at high altitudes

Birds that fly at high altitudes must support vigorous exercise in oxygen-thin environments. Here I discuss the characteristics that help high fliers sustain the high rates of metabolism needed for flight at elevation. Many traits in the O(2) transport pathway distinguish birds in general from other vertebrates. These include enhanced gas-exchange efficiency in the lungs, maintenance of O(2) delivery and oxygenation in the brain during hypoxia, augmented O(2) diffusion capacity in peripheral tissues and a high aerobic capacity. These traits are not high-altitude adaptations, because they are also characteristic of lowland birds, but are nonetheless important for hypoxia tolerance and exercise capacity. However, unique specializations also appear to have arisen, presumably by high-altitude adaptation, at every step in the O(2) pathway of highland species. The distinctive features of high fliers include an enhanced hypoxic ventilatory response, an effective breathing pattern, larger lungs, haemoglobin with a higher O(2) affinity, further augmentation of O(2) diffusion capacity in the periphery and multiple alterations in the metabolic properties of cardiac and skeletal muscle. These unique specializations improve the uptake, circulation and efficient utilization of O(2) during high-altitude hypoxia. High-altitude birds also have larger wings than their lowland relatives to reduce the metabolic costs of staying aloft in low-density air. High fliers are therefore unique in many ways, but the relative roles of adaptation and plasticity (acclimatization) in high-altitude flight are still unclear. Disentangling these roles will be instrumental if we are to understand the physiological basis of altitudinal range limits and how they might shift in response to climate change.

Gene flow and hybridization between numerically imbalanced populations of two duck species in the Falkland Islands

Interspecific hybridization is common in plants and animals, particularly in waterfowl (Anatidae). One factor shown to contribute to hybridization is restricted mate choice, which can occur when two species occur in sympatry but one is rare. The Hubbs principle, or “desperation hypothesis,” states that under such circumstances the rarer species is more likely to mate with heterospecifics. Here we report interspecific hybridization between two waterfowl species that coexist in broad sympatry and mixed flocks throughout southern South America. Speckled teal (Anas flavirostris) and yellow-billed pintails (Anas georgica) are abundant in continental South America, but in the Falkland Islands speckled teal outnumber yellow-billed pintails approximately ten to one. Using eight genetic loci (mtDNA and 7 nuclear introns) coupled with Bayesian assignment tests and relatedness analysis, we identified a speckled teal x yellow-billed pintail F₁ hybrid female and her duckling sired by a male speckled teal. Although our sample in the Falkland Islands was small, we failed to identify unequivocal evidence of hybridization or introgression in a much larger sample from Argentina using a three-population “isolation with migration” coalescent analysis. While additional data are needed to determine if this event in the Falkland Islands was a rare singular occurrence, our results provide further support for the “desperation hypothesis,” which states that scarcity in one population and abundance of another will often lead to hybridization.

A molecular phylogeny of the Sierra-Finches (Phrygilus, Passeriformes): extreme polyphyly in a group of Andean specialists

The unparalleled avian diversity of the Neotropics has long been argued to be in large part the evolutionary consequence of the incredible habitat diversity and rugged topography of the Andes mountains. Various scenarios have been proposed to explain how the Andean context could have generated lineage diversification (e.g. vicariant speciation or parapatric speciation across vertical ecological gradients), yet further study on Andean taxa is needed to reveal the relative importance of the different processes. Here we use mitochondrial and nuclear DNA sequences to derive the first phylogenetic hypothesis for Phrygilus (Sierra-Finches), one of the most species-rich genera of mainly Andean passerines. We find strong evidence that the genus is polyphyletic, comprising four distantly related clades with at least nine other genera interspersed between them (Acanthidops, Catamenia, Diglossa, Haplospiza, Idiopsar, Melanodera, Rowettia, Sicalis and Xenodacnis). These four Phrygilus clades coincide with groups previously established mainly on the basis of plumage characters, suggesting single evolutionary origins for each of these. We consider the history of diversification of each clade, analyzing the timing of splitting events, ancestral reconstruction of altitudinal ranges and current geographical distributions. Phrygilus species origins date mainly to the Pleistocene, with representatives diversifying within, out of, and into the Andes. Finally, we explored whether Phrygilus species, especially those with broad altitudinal and latitudinal Andean distributions, showed phylogeographic structure. Our best-sampled taxon (Phrygilus fruticeti) exhibited no clear pattern; however, we found deep genetic splits within other surveyed species, with Phrygilus unicolor being the most extreme case and deserving of further research.

Adaptive traits are maintained on steep selective gradients despite gene flow and hybridization in the intertidal zone

Gene flow among hybridizing species with incomplete reproductive barriers blurs species boundaries, while selection under heterogeneous local ecological conditions or along strong gradients may counteract this tendency. Congeneric, externally-fertilizing fucoid brown algae occur as distinct morphotypes along intertidal exposure gradients despite gene flow. Combining analyses of genetic and phenotypic traits, we investigate the potential for physiological resilience to emersion stressors to act as an isolating mechanism in the face of gene flow. Along vertical exposure gradients in the intertidal zone of Northern Portugal and Northwest France, the mid-low shore species Fucus vesiculosus, the upper shore species Fucus spiralis, and an intermediate distinctive morphotype of F. spiralis var. platycarpus were morphologically characterized. Two diagnostic microsatellite loci recovered 3 genetic clusters consistent with prior morphological assignment. Phylogenetic analysis based on single nucleotide polymorphisms in 14 protein coding regions unambiguously resolved 3 clades; sympatric F. vesiculosus, F. spiralis, and the allopatric (in southern Iberia) population of F. spiralis var. platycarpus. In contrast, the sympatric F. spiralis var. platycarpus (from Northern Portugal) was distributed across the 3 clades, strongly suggesting hybridization/introgression with both other entities. Common garden experiments showed that physiological resilience following exposure to desiccation/heat stress differed significantly between the 3 sympatric genetic taxa; consistent with their respective vertical distribution on steep environmental clines in exposure time. Phylogenetic analyses indicate that F. spiralis var. platycarpus is a distinct entity in allopatry, but that extensive gene flow occurs with both higher and lower shore species in sympatry. Experimental results suggest that strong selection on physiological traits across steep intertidal exposure gradients acts to maintain the 3 distinct genetic and morphological taxa within their preferred vertical distribution ranges. On the strength of distributional, genetic, physiological and morphological differences, we propose elevation of F. spiralis var. platycarpus from variety to species level, as F. guiryi.

Enhanced AFLP genome scans detect local adaptation in high-altitude populations of a small rodent (Microtus arvalis)

Adaptation to adverse environmental conditions such as high altitude requires physiological and/or morphological changes. Genome scans provide a means to identify the genetic basis of such adaptations without previous knowledge about the particular genetic variants or traits under selection. In this study, we scanned 3027 amplified fragment length polymorphisms (AFLP) in four populations of the common vole Microtus arvalis for loci associated with local adaptation and high altitude. We investigated voles from two populations at high elevation (~2000 m a.s.l.) representing the upper limit of the altitudinal distribution of the species and two geographically close low-altitude populations (<600 m a.s.l.). Statistical analysis incorporated a new Bayesian F(ST) outlier approach specifically developed for AFLP markers, which considers the intensity of AFLP bands instead of mere presence/absence and allows to derive population-based estimates of allele frequencies and F(IS) values. Computer simulations showed that this approach increases the statistical power of the detection of AFLP markers under selection almost to the power of single nucleotide polymorphism (SNP) data without compromising specificity. Our enhanced genome scan resulted in 20 prime candidate markers for positive selection, which show mostly extremely high allele frequency differences between the low- and high-altitude populations. The comparison of global- and pairwise-enhanced genome scans demonstrated further that very strong selective signatures may also be associated with single populations suggesting the importance of local adaptation in alpine populations of common voles.