Consistent changes in muscle phenotype and mitochondrial abundance underlie dive performance across multiple lineages of diving ducks.

Diving animals must sustain high muscle activity with finite oxygen (O2) to forage underwater. Studies have shown that some diving mammals exhibit changes in the metabolic phenotype of locomotory muscles compared to non-divers, but the pervasiveness of such changes across diving animals is unclear, particularly among diving birds. Here, we examine whether changes in muscle phenotype and mitochondrial abundance are associated with dive capacity across 17 species of ducks from three distinct evolutionary clades (tribes) in the subfamily Anatinae - the longest diving sea ducks, the mid-tier diving pochards, and the non-diving dabblers. In the gastrocnemius (the primary swimming and diving muscle), mitochondrial volume density in both oxidative and glycolytic fiber types were 70% and 30% higher in sea ducks compared to dabblers, respectively. These differences were associated with preferential proliferation of the subsarcolemmal subfraction, the mitochondria adjacent to the cell membrane and nearest to capillaries, relative to the intermyofibrillar subfraction. Capillary density and capillary-to-fiber ratio were positively correlated with mitochondrial volume density, with no variation in the density of oxidative fiber types across tribes. In the pectoralis, sea ducks had greater abundance of oxidative fiber types than dabblers, whereas pochards were intermediate between the two. These data suggest that skeletal muscles of sea ducks have a heightened capacity for aerobic metabolism and an enhanced ability to utilize O2 stores in the blood and muscle while diving.

Consistent changes in muscle metabolism underlie dive performance across multiple lineages of diving ducks.

Diving animals must sustain high activity with limited O2 stores to successfully capture prey. Studies suggest that increasing body O2 stores supports breath-hold diving, but less is known about metabolic specializations that underlie underwater locomotion. We measured maximal activities of 10 key enzymes in locomotory muscles (gastrocnemius and pectoralis) to identify biochemical changes associated with diving in pathways of oxidative and substrate-level phosphorylation and compared them across three groups of ducks-the longest diving sea ducks (eight spp.), the mid-tier diving pochards (three spp.) and the non-diving dabblers (five spp.). Relative to dabblers, both diving groups had increased activities of succinate dehydrogenase and cytochrome c oxidase, and sea ducks further showed increases in citrate synthase (CS) and hydroxyacyl-CoA dehydrogenase (HOAD). Both diving groups had relative decreases in capacity for anaerobic metabolism (lower ratio of lactate dehydrogenase to CS), with sea ducks also showing a greater capacity for oxidative phosphorylation and lipid oxidation (lower ratio of pyruvate kinase to CS, higher ratio of HOAD to hexokinase). These data suggest that the locomotory muscles of diving ducks are specialized for sustaining high rates of aerobic metabolism, emphasizing the importance of body O2 stores for dive performance in these species.

Population genomics indicate three different modes of divergence and speciation with gene flow in the green-winged teal duck complex.

The processes leading to divergence and speciation can differ broadly among taxa with different life histories. We examine these processes in a small clade of ducks with historically uncertain relationships and species limits. The green-winged teal (Anas crecca) complex is a Holarctic species of dabbling duck currently categorized as three subspecies (Anas crecca crecca, A. c. nimia, and A. c. carolinensis) with a close relative, the yellow-billed teal (Anas flavirostris) from South America. A. c. crecca and A. c. carolinensis are seasonal migrants, while the other taxa are sedentary. We examined divergence and speciation patterns in this group, determining their phylogenetic relationships and the presence and levels of gene flow among lineages using both mitochondrial and genome-wide nuclear DNA obtained from 1,393 ultraconserved element (UCE) loci. Phylogenetic relationships using nuclear DNA among these taxa showed A. c. crecca, A. c. nimia, and A. c. carolinensis clustering together to form one polytomous clade, with A. flavirostris sister to this clade. This relationship can be summarized as (crecca, nimia, carolinensis)(flavirostris). However, whole mitogenomes revealed a different phylogeny: (crecca, nimia)(carolinensis, flavirostris). The best demographic model for key pairwise comparisons supported divergence with gene flow as the probable speciation mechanism in all three contrasts (crecca-nimia, crecca-carolinensis, and carolinensis-flavirostris). Given prior work, gene flow was expected among the Holarctic taxa, but gene flow between North American carolinensis and South American flavirostris (M âˆ¼0.1-0.4 individuals/generation), albeit low, was not expected. Three geographically oriented modes of divergence are likely involved in the diversification of this complex: heteropatric (crecca-nimia), parapatric (crecca-carolinensis), and (mostly) allopatric (carolinensis-flavirostris). Our study shows that ultraconserved elements are a powerful tool for simultaneously studying systematics and population genomics in systems with historically uncertain relationships and species limits.

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.

Assessing changes in genomic divergence following a century of human-mediated secondary contact among wild and captive-bred ducks.

Along with manipulating habitat, the direct release of domesticated individuals into the wild is a practice used worldwide to augment wildlife populations. We test between possible outcomes of human-mediated secondary contact using genomic techniques at both historical and contemporary timescales for two iconic duck species. First, we sequence several thousand ddRAD-seq loci for contemporary mallards (Anas platyrhynchos) throughout North America and two domestic mallard types (i.e., known game-farm mallards and feral Khaki Campbell's). We show that North American mallards may well be becoming a hybrid swarm due to interbreeding with domesticated game-farm mallards released for hunting. Next, to attain a historical perspective, we applied a bait-capture array targeting thousands of loci in century-old (1842-1915) and contemporary (2009-2010) mallard and American black duck (Anas rubripes) specimens. We conclude that American black ducks and mallards have always been closely related, with a divergence time of ~600,000 years before present, and likely evolved through prolonged isolation followed by limited bouts of gene flow (i.e., secondary contact). They continue to maintain genetic separation, a finding that overturns decades of prior research and speculation suggesting the genetic extinction of the American black duck due to contemporary interbreeding with mallards. Thus, despite having high rates of hybridization, actual gene flow is limited between mallards and American black ducks. Conversely, our historical and contemporary data confirm that the intensive stocking of game-farm mallards during the last ~100 years has fundamentally changed the genetic integrity of North America's wild mallard population, especially in the east. It thus becomes of great interest to ask whether the iconic North American mallard is declining in the wild due to introgression of maladaptive traits from domesticated forms. Moreover, we hypothesize that differential gene flow from domestic game-farm mallards into the wild mallard population may explain the overall temporal increase in differentiation between wild black ducks and mallards, as well as the uncoupling of genetic diversity and effective population size estimates across time in our results. Finally, our findings highlight how genomic methods can recover complex population histories by capturing DNA preserved in traditional museum specimens.

A morphometric analysis of the lungs of high-altitude ducks and geese.

We examined the morphology of the lungs of five species of high-altitude resident ducks from Lake Titicaca in the Peruvian Andes (yellow-billed pintail [Anas georgica], cinnamon teal [Anas cyanoptera orinomus], puna teal [Anas puna], speckled teal [Anas flavirostris oxyptera], and ruddy duck [Oxyura jamaicensis ferruginea]) and compared them with those of the high-altitude migratory bar-headed goose (Anser indicus) and the low-altitude migratory barnacle goose (Branta leucopsis). We then determined the relationship between mass-specific lung volume, the volume densities of the component parts of the lung, and previously reported hypoxia-induced increases in pulmonary O2 extraction. We found that the mass-specific lung volumes and the mass-specific volume of the exchange tissue were larger in the lungs of high-altitude resident birds. The bar-headed goose had a mass-specific lung volume that fell between those of the low-altitude species and the high-altitude residents, but a mass-specific volume of exchange tissue that was not significantly different than that of the high-altitude residents. The data suggest that the mass-specific volume of the lung may increase with evolutionary time spent at altitude. We found an inverse relationship between the percentage increase in pulmonary O2 extraction and the percentage increase in ventilation across species that was independent of the volume density of the exchange tissue, at least for the resident Andean birds.

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.

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 examined the cardiovascular responses to progressive hypoxia in four high-altitude Andean duck species compared with 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 (Hb) concentration. 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 concentration, increasing the O2 content of arterial blood, and potentially providing a greater reserve for enhancing O2 delivery during hypoxia. Hct and Hb concentration 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 concentration compared with their low-altitude relatives, revealing yet another avian strategy for coping with high altitude.

ddRAD-seq data reveal significant genome-wide population structure and divergent genomic regions that distinguish the mallard and close relatives in North America.

Recently evolved species typically share genetic variation across their genomes due to incomplete lineage sorting and/or ongoing gene flow. Given only subtle allele frequency differences at most loci and the expectation that divergent selection may affect only a tiny fraction of the genome, distinguishing closely related species based on multi-locus data requires substantial genomic coverage. In this study, we used ddRAD-seq to sample the genomes of five recently diverged, New World "mallards" (Anas spp.), a group of dabbling duck species characterized by diagnosable phenotypic differences but minimal genetic differentiation. With increased genomic sampling, we aimed to characterize population structure within this group and identify genomic regions that may have experienced divergent selection during speciation. We analyzed 3,017 autosomal ddRAD-seq loci and 177 loci from the Z-chromosome. In contrast to previous studies, the ddRAD-seq data were sufficient to assign individuals to their respective species or subspecies and to generate estimates of gene flow in a phylogenetic framework. We find limited evidence of contemporary gene flow between the dichromatic mallard and several monochromatic taxa, but find evidence for historical gene flow between some monochromatic species pairs. We conclude that the overall genetic similarity of these taxa likely reflects retained ancestral polymorphism rather than recent and extensive gene flow. Thus, despite recurring cases of hybridization in this group, our results challenge the current dogma predicting the genetic extinction of the New World monochromatic dabbling ducks via introgressive hybridization with mallards. Moreover, ddRAD-seq data were sufficient to identify previously unknown outlier regions across the Z-chromosome and several autosomal chromosomes that may have been involved in the diversification of species in this recent radiation.

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

During periods of reduced O2 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-O2 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 O2 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.

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 O2 extraction. O2 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 O2 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 O2 extraction compared with low-altitude relatives. Arterial O2 saturation (SaO2 ) was elevated in yellow-billed pintails at moderate levels of hypoxia, but there were no differences in SaO2  in other high-altitude taxa compared with their close relatives. This finding suggests that improvements in SaO2  in hypoxia can require increases in both breathing and haemoglobin-O2 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.

Respiratory mechanics and morphology of Tibetan and Andean high-altitude geese with divergent life histories.

High-altitude bar-headed geese (Anser indicus) and Andean geese (Chloephaga melanoptera) have been shown to preferentially increase tidal volume over breathing frequency when increasing ventilation during exposure to hypoxia. Increasing tidal volume is a more effective breathing strategy but is also thought to be more mechanically and metabolically expensive. We asked whether there might be differences in the mechanics or morphology of the respiratory systems of high-altitude transient bar-headed geese and high-altitude resident Andean geese that could minimize the cost of breathing more deeply. We compared these two species with a low-altitude migratory species, the barnacle goose (Branta leucopsis). We ventilated anesthetized birds to measure mechanical properties of the respiratory system and used CT scans to quantify respiratory morphology. We found that the respiratory system of Andean geese was disproportionately larger than that of the other two species, allowing use of a deeper breathing strategy for the same energetic cost. The relative size of the respiratory system, especially the caudal air sacs, of bar-headed geese was also larger than that of barnacle geese. However, when normalized to respiratory system size, the mechanical cost of breathing did not differ significantly among these three species, indicating that deeper breathing is enabled by morphological but not mechanical differences between species. The metabolic cost of breathing was estimated to be <1% of basal metabolic rate at rest in normoxia. Because of differences in the magnitude of the ventilatory response, the cost of breathing was estimated to increase 7- to 10-fold in bar-headed and barnacle geese in severe hypoxia, but less than 1-fold in Andean geese exposed to the same low atmospheric PO2.

Convergent Evolution of Hemoglobin Function in High-Altitude Andean Waterfowl Involves Limited Parallelism at the Molecular Sequence Level.

A fundamental question in evolutionary genetics concerns the extent to which adaptive phenotypic convergence is attributable to convergent or parallel changes at the molecular sequence level. Here we report a comparative analysis of hemoglobin (Hb) function in eight phylogenetically replicated pairs of high- and low-altitude waterfowl taxa to test for convergence in the oxygenation properties of Hb, and to assess the extent to which convergence in biochemical phenotype is attributable to repeated amino acid replacements. Functional experiments on native Hb variants and protein engineering experiments based on site-directed mutagenesis revealed the phenotypic effects of specific amino acid replacements that were responsible for convergent increases in Hb-O2 affinity in multiple high-altitude taxa. In six of the eight taxon pairs, high-altitude taxa evolved derived increases in Hb-O2 affinity that were caused by a combination of unique replacements, parallel replacements (involving identical-by-state variants with independent mutational origins in different lineages), and collateral replacements (involving shared, identical-by-descent variants derived via introgressive hybridization). In genome scans of nucleotide differentiation involving high- and low-altitude populations of three separate species, function-altering amino acid polymorphisms in the globin genes emerged as highly significant outliers, providing independent evidence for adaptive divergence in Hb function. The experimental results demonstrate that convergent changes in protein function can occur through multiple historical paths, and can involve multiple possible mutations. Most cases of convergence in Hb function did not involve parallel substitutions and most parallel substitutions did not affect Hb-O2 affinity, indicating that the repeatability of phenotypic evolution does not require parallelism at the molecular level.

Divergent respiratory and cardiovascular responses to hypoxia in bar-headed geese and Andean birds.

Many high-altitude vertebrates have evolved increased capacities in their oxygen transport cascade (ventilation, pulmonary diffusion, circulation and tissue diffusion), enhancing oxygen transfer from the atmosphere to mitochondria. However, the extent of interspecies variation in the control processes that dictate hypoxia responses remains largely unknown. We compared the metabolic, cardiovascular and respiratory responses to progressive decreases in inspired oxygen levels of bar-headed geese (Anser indicus), birds that biannually migrate across the Himalayan mountains, with those of Andean geese (Chloephaga melanoptera) and crested ducks (Lophonetta specularioides), lifelong residents of the high Andes. We show that Andean geese and crested ducks have evolved fundamentally different mechanisms for maintaining oxygen supply during low oxygen (hypoxia) from those of bar-headed geese. Bar-headed geese respond to hypoxia with robust increases in ventilation and heart rate, whereas Andean species increase lung oxygen extraction and cardiac stroke volume. We propose that transient high-altitude performance has favoured the evolution of robust convective oxygen transport recruitment in hypoxia, whereas life-long high-altitude residency has favoured the evolution of structural enhancements to the lungs and heart that increase lung diffusion and stroke volume.

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

Becoming pure: identifying generational classes of admixed individuals within lesser and greater scaup populations.

Estimating the frequency of hybridization is important to understand its evolutionary consequences and its effects on conservation efforts. In this study, we examined the extent of hybridization in two sister species of ducks that hybridize. We used mitochondrial control region sequences and 3589 double-digest restriction-associated DNA sequences (ddRADseq) to identify admixture between wild lesser scaup (Aythya affinis) and greater scaup (A. marila). Among 111 individuals, we found one introgressed mitochondrial DNA haplotype in lesser scaup and four in greater scaup. Likewise, based on the site-frequency spectrum from autosomal DNA, gene flow was asymmetrical, with higher rates from lesser into greater scaup. However, using ddRADseq nuclear DNA, all individuals were assigned to their respective species with >0.95 posterior assignment probability. To examine the power for detecting admixture, we simulated a breeding experiment in which empirical data were used to create F1 hybrids and nine generations (F2-F10) of backcrossing. F1 hybrids and F2, F3 and most F4 backcrosses were clearly distinguishable from pure individuals, but evidence of admixed histories was effectively lost after the fourth generation. Thus, we conclude that low interspecific assignment probabilities (0.011-0.043) for two lesser and nineteen greater scaup were consistent with admixed histories beyond the F3 generation. These results indicate that the propensity of these species to hybridize in the wild is low and largely asymmetric. When applied to species-specific cases, our approach offers powerful utility for examining concerns of hybridization in conservation efforts, especially for determining the generational time until admixed histories are effectively lost through backcrossing.

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.

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.

Phylogenetics of a recent radiation in the mallards and allies (Aves: Anas): inferences from a genomic transect and the multispecies coalescent.

Reconstructing species trees by incorporating information from many independent gene trees reduces the confounding influence of stochastic lineage sorting. Such analyses are particularly important for taxa that share polymorphisms due to incomplete lineage sorting or introgressive hybridization. We investigated phylogenetic relationships among 14 closely related taxa from the mallard (Anas spp.) complex using the multispecies coalescent and 20 nuclear loci sampled from a genomic transect. We also examined how treating recombining loci and hybridizing species influences results by partitioning the data using various protocols. In general, topologies were similar among the various species trees, with major clades consistently composed of the same taxa. However, relationships among these clades and among taxa within clades changed among partitioned data sets. Posterior support generally decreased when filtering for recombination, whereas excluding mallards (Anas platyrhynchos) increased posterior support for taxa known to hybridize with them. Furthermore, branch lengths decreased substantially for recombination-filtered data. Finally, concordance between nuclear and morphometric topologies conflicted with those in the mitochondrial tree, particularly with regard to the placement of the Hawaiian duck (A. wyvilliana), Philippine duck (A. luzonica), and two spot-billed ducks (A. zonorhyncha and A. poecilorhyncha). These results demonstrate the importance of maximizing sequence length and taxon sampling when inferring taxonomic relationships that are confounded by extensive allele sharing.