Search

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, andA. c. carolinensis) with a close relative, the yellow-billed teal (Anas flavirostris) from South America. 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 showedA. c. crecca,A. c. nimia, andA. c. carolinensisclustering together to form one polytomous clade, withA.flavirostrissister 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, andcarolinensis−flavirostris). Given prior work, gene flow was expected among the Holarctic taxa, but gene flow between North Americancarolinensisand South Americanflavirostris(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). Ultraconserved elements are a powerful tool for simultaneously studying systematics and population genomics in systems like this.Graphical Abstract

Birds living in riverine environments may show weak population structure because high dispersal abilities required to track habitat dynamics can result in gene flow over broad spatial scales. Alternatively, the configuration of river networks may result in restricted dispersal within river courses or basins, leading to high genetic structure. Although several bird species are riverine specialists in the Andes, no study has extensively evaluated the population genetic structure of any of them. We examined evidence from genetic and morphological data to address questions about the biogeography and taxonomy of the Torrent Duck (Merganetta armata), a riverine specialist bird with a broad distribution in Andean riverine habitats which certainly comprises different subspecies and may comprise more than one species. We found deep subdivisions of Torrent Duck populations from the northern, central and southern portions of the Andes. These lineages, which partly coincide with subspecies described based on plumage variation and body size, do not share mtDNA haplotypes, have private nuclear alleles and exhibit marked differences in morphometric traits. Some geographic barriers presumably restricting gene flow between groups partially coincide with those associated with major genetic breaks in forest species with similar distributions along the Andes, suggesting that bird assemblages including species occupying different habitats were likely affected by common biogeographical events. The three groups of Torrent Ducks may be considered different species under some species definitions and are distinct evolutionary lineages to be conserved and managed separately. Fil: Gutiérrez Pinto, Natalia. University of Nebraska; Estados Unidos. Universidad de los Andes; Colombia Fil: McCracken, Kevin G.. University of Alaska; Estados Unidos. University of Miami; Estados Unidos. University Of Miami. Rosenstiel School Of Marine Atmospheric Science; Estados Unidos Fil: Tubaro, Pablo Luis. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Parque Centenario. Museo Argentino de Ciencias Naturales "Bernardino Rivadavia"; Argentina Fil: Kopuchian, Cecilia. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Nordeste. Centro de Ecología Aplicada del Litoral. Universidad Nacional del Nordeste. Centro de Ecología Aplicada del Litoral; Argentina. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Parque Centenario. Museo Argentino de Ciencias Naturales "Bernardino Rivadavia"; Argentina Fil: Astié, Andrea Alejandra. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Mendoza. Instituto Argentino de Investigaciones de las Zonas Áridas. Provincia de Mendoza. Instituto Argentino de Investigaciones de las Zonas Áridas. Universidad Nacional de Cuyo. Instituto Argentino de Investigaciones de las Zonas Áridas; Argentina Fil: Cadena, Carlos Daniel. Universidad de los Andes; Colombia

The Eurasian wigeon (Mareca penelope) is one of the most numerous migrant species of waterfowl in the Palearctic. Annually, significant part of the world’s wigeon population makes seasonal flights over distances from tens to thousands or more kilometers. According to different estimates based on banding data, five geographic populations of the species were described in the Palearctic. However, distinct borders between the populations have not been identified. At the same time, no phylogeographic studies have been carried out for the complete native range of wigeon so far. In addition to the fundamental importance of such a study, knowledge of the genetic structure of populations is necessary for the development of measures to increase the number of and preserve this valuable game species. The aim of our work was a phylogeographic analysis of the wigeon across its vast native range in the Palearctic including ducks wintering in North America. We examined genetic diversity and differentiation of wigeon populations identified with banding data, phylogenetic relationships of mtDNA haplotypes and demographic history of populations and species as a whole by sequencing a 661 base-pair 5’-fragment of the mitochondrial control region from 195 individual ducks collected throughout the Palearctic and Nearctic. Genetic diversity was high in all studied populations. A reconstruction of haplotypes phylogeny revealed the absence of geographic structure in the data. Nonetheless, analysis of molecular variance (AMOVA) identified two groups of populations: EuropeanSiberian and East Asian. The former included wigeons from Europe, Siberia and the Atlantic coast of North America, and the latter comprised ducks from Russian Far East, Kamchatka Peninsula, Chukotka Autonomous District, the Aleutian Islands, Alaska, and the Pacific coast of 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.

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-O(2) 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-O(2) 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.

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.

In this comprehensive survey of microbiomes of >900 species, including 315 mammals and 491 birds, we find a striking convergence of the microbiomes of birds and animals that fly. In nonflying mammals, diet and short-term evolutionary relatedness drive the microbiome, and many microbial species are specific to a particular kind of mammal, but flying mammals and birds break this pattern with many microbes shared across different species, with little correlation either with diet or with relatedness of the hosts. This finding suggests that adaptation to flight breaks long-held relationships between hosts and their microbes., Diet and host phylogeny drive the taxonomic and functional contents of the gut microbiome in mammals, yet it is unknown whether these patterns hold across all vertebrate lineages. Here, we assessed gut microbiomes from ∼900 vertebrate species, including 315 mammals and 491 birds, assessing contributions of diet, phylogeny, and physiology to structuring gut microbiomes. In most nonflying mammals, strong correlations exist between microbial community similarity, host diet, and host phylogenetic distance up to the host order level. In birds, by contrast, gut microbiomes are only very weakly correlated to diet or host phylogeny. Furthermore, while most microbes resident in mammalian guts are present in only a restricted taxonomic range of hosts, most microbes recovered from birds show little evidence of host specificity. Notably, among the mammals, bats host especially bird-like gut microbiomes, with little evidence for correlation to host diet or phylogeny. This suggests that host-gut microbiome phylosymbiosis depends on factors convergently absent in birds and bats, potentially associated with physiological adaptations to flight. Our findings expose major variations in the behavior of these important symbioses in endothermic vertebrates and may signal fundamental evolutionary shifts in the cost/benefit framework of the gut microbiome.

Flightlessness in birds is the product of changes in suites of characters-including increased body size and reduced anterior limbs-that have evolved repeatedly and independently under similar ecological conditions (generally insularity). It remains unknown whether this phenotypic convergence extends to the genomic level, partially because many losses of flight occurred long ago (such as in penguins or ratites), thus complicating the study of the genetic pathways to flightlessness. Here, we use genome sequencing to study the evolution of flightlessness in a group of ducks that are current and dynamic exemplars of this major functional transition. These recently diverged Tachyeres steamer ducks differ in their ability to fly: one species is predominantly flighted and three are mainly flightless. Through a genome-wide association analysis, we identify two narrow candidate genomic regions implicated in the morphological changes that led to flightlessness, and reconstruct the number of times flightlessness has evolved in Tachyeres. The strongest association is with DYRK1A, a gene that when knocked out in mice leads to alterations in growth and bone morphogenesis. These findings, together with phylogenetic and demographic analyses, imply that the genomic changes leading to flightlessness in Tachyeres may have evolved once, and that this trait remains functionally polymorphic in two species.

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 systemsthat have recently colonized extremeenvironments. 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 geneflowacross 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 DNAdamage, 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

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.

Under drift-mutation equilibrium, genetic diversity is expected to be correlated with effective population size (Ne). Changes in population size and gene flow are two important processes that can cause populations to deviate from this expected relationship. In this study, we used DNA sequences from six independent loci to examine the influence of these processes on standing genetic diversity in endemic mottled ducks (Anas fulvigula) and geographically widespread mallards (A. platyrhynchos), two species known to hybridize. Mottled ducks have an estimated census size that is about two orders-of-magnitude smaller than that of mallards, yet these two species have similar levels of genetic diversity, especially at nuclear DNA. Coalescent analyses suggest that a population expansion in the mallard at least partly explains this discrepancy, but the mottled duck harbors higher genetic diversity and apparent Ne than expected for its census size even after accounting for a population decline. Incorporating gene flow into the model, however, reduced the estimated Ne of mottled ducks to 33 % of the equilibrium Ne and yielded an estimated Ne consistent with census size. We also examined the utility of these loci to distinguish among mallards, mottled ducks, and their hybrids. Most putatively pure individuals were correctly assigned to species, but the power for detecting hybrids was low. Although hybridization with mallards potentially poses a conservation threat to mottled ducks by creating a risk of extinction by hybridization, introgression of mallard alleles has helped maintain high genetic diversity in mottled ducks and might be important for the adaptability and survival of this species.

Heteropatric differentiation is a mode of speciation with gene flow in which divergence occurs between lineages that are in sympatry and allopatry at different times during cyclic spatial movements. Empirical evidence suggests that heteropatric differentiation may prove to be common among seasonally migratory organisms. We examined genetic differentiation between the sedentary Aleutian Islands population of green-winged teal (Anas crecca-nimia) and its close migratory relative, the Eurasian, or Old World (OW), Anas c. crecca population, a portion of which passes through the range of nimia during its seasonal migrations. We also examined its relationship with the parapatric North American, New World (NW), A. c. carolinensis population. Sequence data from eight nuclear introns and the mtDNA control region showed that the nimia-crecca divergence occurred much more recently than the deeper crecca-carolinensis split (~83 000 years vs. ~1.1 Myr). Despite considerable spatial overlap between crecca and nimia during seasonal migration, three key predictions of heteropatric differentiation are supported: significant genetic divergence (overall mean Φst = 0.07), low gene flow (2Ne m ~ 1.8), and an effective population size in nimia that is not especially low (Ne ~ 80 000 individuals). Similar levels of gene flow have come into nimia from carolinensis, but no detectable nuclear gene flow has gone out of nimia into either OW (crecca) or NW (carolinensis) populations. We infer that adaptations of these populations to local optima in different places (e.g. each matching their reproductive effort to different resource blooms) promote genetic isolation and divergence despite periods of sympatry between them, as the heteropatric model predicts.

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.

Speciation is a process in which genetic drift and selection cause divergence over time. However, there is no rule dictating the time required for speciation, and even low levels of gene flow hinder divergence, so that taxa may be poised at the threshold of speciation for long periods of evolutionary time. We sequenced mitochondrial DNA (mtDNA) and eight nuclear introns (nuDNA) to estimate genomic levels of differentiation and gene flow between the Eurasian common teal (Anas crecca crecca) and the North American green-winged teal (Anas crecca carolinensis). These ducks come into contact in Beringia (north-eastern Asia and north-western North America) and have probably done so, perhaps cyclically, since the Pliocene-Pleistocene transition, ~2.6 Ma, when they apparently began diverging. They have diagnosable differences in male plumage and are 6.9% divergent in the mtDNA control region, with only 1 of 58 crecca and 2 of 86 carolinensis having haplotypes grouping with the other. Two nuclear loci were likewise strongly structured between these teal (Φ(st) ≥ 0.35), but six loci were undifferentiated or only weakly structured (Φ(st) = 0.0-0.06). Gene flow between crecca and carolinensis was ~1 individual per generation in both directions in mtDNA, but was asymmetrical in nuDNA, with ~1 and ~20 individuals per generation immigrating into crecca and carolinensis, respectively. This study illustrates that species delimitation using a single marker oversimplifies the complexity of the speciation process, and it suggests that even with divergent selection, moderate levels of gene flow may stall the speciation process short of completion.

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.

Two species of waterfowl living at high altitude provide a prominent example of parallel adaptation at the molecular level. The bar-headed goose (Anser indicus) breeds at high elevations in central Asia and migrates across the Himalayas, where the partial pressure of oxygen (O(2)) is one-third of sea level. In South America, the distantly related Andean goose (Chloephaga melanoptera) is endemic to the high Andes. Both species exhibit increased blood-O(2) affinity, which has been attributed to the effects of single amino acid substitutions in the major hemoglobin. Here we present phylogenetic analyses of the swans and geese (Anserinae) and South American sheldgeese (Anatinae) using the three genes that encode the major (HbA) and minor (HbD) hemoglobin isoforms. We sought to determine whether two amino acid substitutions that have been the focus of extensive biochemical analysis (Ala-alpha(A)119 and Ser-beta(A)55) are uniquely derived in bar-headed goose and Andean goose, respectively, and to examine evidence of molecular adaptation at other positions in hemoglobin genes by comparing these two high-altitude taxa to their closest relatives. Bayesian analysis of the alpha(A)-, alpha(D)-, and beta(A)-subunit genes produced well-resolved phylogenies, with high posterior probabilities and bootstrap values for most genera. The bar-headed goose is likely sister to all other Anser species. Andean goose, the sole highland representative of the South American sheldgeese is either sister to the other Chloephaga species or sister to Neochen. In the bar-headed goose, four derived substitutions were observed in HbA (alpha(A)12, 18, 63, 119) and two in HbD (alpha(D)2, 47). Four derived substitutions in Andean goose include three in HbA (alpha(A)8, 77; beta(A)86) and two in HbD (alpha(D)9; beta(A)86). Considering both highland species, four substitutions (Ala-alpha(A)8, Ala-alpha(A)12, Ser-alpha(A)18, Leu-alpha(D)9) were located at adjacent positions on the A helix (or AB corner) of the alpha-chains, three others (Thr-alpha(A)77, Ser-beta(A)86, Ser-alpha(D)2) were in close proximity to inositolpentaphosphate (IP(5)) binding sites, and Ala-alpha(A)119 occurred at an alphabeta intersubunit contact. Ser-beta(A)55, which is involved in the same alphabeta intersubunit contact and was previously shown to increase Hb-O(2) affinity, is not unique to Andean goose, but is a synapomorphy of the South American sheldgeese, a clade of predominantly lowland waterfowl. Our findings illustrate the importance of understanding phylogenetic relationships and polarity of character-state changes when making inferences about adaptive evolution.

Musk Ducks (Biziura lobata) are endemic to Australia and occur as two geographically isolated populations separated by the Nullarbor Plain, a vast arid region in southern Australia. We studied genetic variation in Musk Duck populations at coarse (eastern versus western Australia) and fine scales (four sites within eastern Australia). We found significant genetic structure between eastern and western Australia in the mtDNA control region (ΦST = 0.747), one nuclear intron (ΦST = 0.193) and eight microsatellite loci (FST = 0.035). In contrast, there was little genetic structure between Kangaroo Island and adjacent mainland regions within eastern Australia. One small population of Musk Ducks in Victoria (Lake Wendouree) differed from both Kangaroo Island and the remainder of mainland eastern Australia, possibly due to genetic drift exacerbated by inbreeding and small population size. The observed low pairwise distance between the eastern and western mtDNA lineages (0.36%) suggests that they diverged near the end of the Pleistocene, a period characterised by frequent shifts between wet and arid conditions in central Australia. Our genetic results corroborate the display call divergence and Mathews’ (Austral Avian Record 2:83–107, 1914) subspecies classification, and confirm that eastern and western populations of Musk Duck are currently isolated from each other.

We investigated local genetic associations among female Pacific common eiders (Somateria mollissima v-nigrum) nesting in a stochastic Arctic environment within two groups of barrier islands (Simpson Lagoon and Mikkelsen Bay) in the Beaufort Sea, Alaska. Nonrandom genetic associations were observed among nesting females using regional spatial autocorrelation analyses for distance classes up to 1000 m in Simpson Lagoon. Nearest-neighbour analyses identified clusters of genetically related females with positive lr values observed for 0-13% and 0-7% of the comparisons in Simpson Lagoon and Mikkelsen Bay, respectively, across years. These results indicate that a proportion of females are nesting in close proximity to more genetically related individuals, albeit at low frequency. Such kin groupings may form through active association between relatives or through natal philopatry and breeding site fidelity. Eiders nest in close association with driftwood, which is redistributed annually by seasonal storms. Yet, genetic associations were still observed. Microgeographic structure may thus be more attributable to kin association than natal philopatry and site fidelity. However, habitat availability may also influence the level of structure observed. Regional structure was present only within Simpson Lagoon and this island group includes at least three islands with sufficient driftwood for colonies, whereas only one island at Mikkelsen Bay has these features. A long-term demographic study is needed to understand more fully the mechanisms that lead to fine-scale genetic structure observed in common eiders breeding in the Beaufort Sea.

Hemoglobins play a key role in oxygen transport. High-oxygen-affinity hemoglobins are adaptive in hypoxic environments. To better understand adaptation to high-altitude hypoxia, we extracted RNA from the bone marrow of six Crested Ducks (Lophonetta specularioides) inhabiting the central high Andes of Peru (4,218–4,605 m elevation) and sequenced >2,000 expressed sequence tags (EST) from a non-normalized complementary-DNA (cDNA) library. Overall, we identified 1,692 ESTs in the expression profile representing 462 different genes. Among those, the ESTs that occurred at the highest frequency were the αA (major) hemoglobin subunit (HBA2; 22.7%), leukocyte cell-derived chemotaxin 2 (LECT2; 10%), αD (minor) hemoglobin subunit (HBA1; 9.6%), beta defensins (DEFB; 6.3%), and the βA hemoglobin subunit (HBB; 3.7%). These results provide the first quantitative identification of gene expression in bone marrow of individuals inhabiting high-altitude regions and are in agreement with the known hemopoietic and i...

Moose (Alces alces) are highly mobile mammals that occur across arboreal regions of North America, Europe, and Asia. Alaskan moose (Alces alces gigas) range across much of Alaska and are primary herbivore consumers, exerting a prominent influence on ecosystem structure and functioning. Increased knowledge gained from population genetics provides insights into their population dynamics, history, and dispersal of these unique large herbivores and can aid in conservation efforts. We examined the genetic diversity and population structure of moose (n = 141) with 8 polymorphic microsatellites from 6 regions spanning much of Alaska. Expected heterozygosity was moderate (H(E) = 0.483-0.612), and private alleles ranged from 0 to 6. Both F(ST) and R(ST) indicated significant population structure (P < 0.001) with F(ST) < 0.109 and R(ST) < 0.125. Results of analyses from STRUCTURE indicated 2 prominent population groups, a mix of moose from the Yakutat and Tetlin regions versus all other moose, with slight substructure observed among the second population. Estimates of dispersal differed between analytical approaches, indicating a high level of historical or current gene flow. Mantel tests indicated that isolation-by-distance partially explained observed structure among moose populations (R(2) = 0.45, P < 0.01). Finally, there was no evidence of bottlenecks either at the population level or overall. We conclude that weak population structure occurs among moose in Alaska with population expansion from interior Alaska westward toward the coast.

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.

We evaluated the effects of selection on the molecular evolution of the functional domains of the mammalian cytochrome b gene as it relates to physicochemical properties shown to correlate with rates of amino acid replacement. Two groups of mammals were considered: pocket gophers of the rodent family Geomyidae, and cetaceans and ungulates of the monophyletic taxon Cetartopdactyla. Several characteristics of cytochrome b evolution were common to both mammal groups. The evolution of the matrix domain reflected the region's relative lack of function. Goodness of fit to neutral expectations indicated that external influences have had very little effect on the evolution of the matrix, although in some cases conservative and moderate changes have been favored. Although rates of synonymous nucleotide substitution have been relatively high, the transmembrane domain exhibited poor goodness of fit to neutral expectations. However, the evolution of the transmembrane domain has been constrained by negative selection, allowing a preponderance of conservative and moderate amino acid replacements. We hypothesize that a high rate of substitution is maintained in spite of negative selection because the codons of the transmembrane coding region are predisposed to conservative changes in all amino acid properties. The evolutionary patterns of the intermembrane domain in pocket gophers and cetartiodactyls, however, were very different. Changes inferred from the pocket gopher phylogenetic tree exhibited a significant fit to neutral expectations for each of the amino acid properties. Changes inferred from the cetartiodactyl tree exhibited significant fit to neutral expectations for polarity and isoelectric point, but not for composition, molecular volume, polar requirement, or hydropathy. In each case, lack of fit was due to selection that promoted conservative or moderate change, with the noteworthy exception of polar requirement. We detected an unexpectedly large change in polar requirement (from aspartic acid to threonine) in two separate lineages (Camelus bactrianus and all cetaceans) at amino acid position 159. This inferred change occurred in a region of the cyt-b protein that directly interacts with external surface proteins of the cytochrome bc(1) complex and resulted in a reversion to a more common character state in vertebrates.

The mottled duck (Anas fulvigula) is a year-round endemic resident of the Gulf Coast and one of two non-migratory dabbling ducks that inhabit North America. To investigate population genetic structure of allopatric mottled duck populations, we collected 5’ control region sequences (bp 78–774) from the mitochondria of 219 mottled ducks sampled at 11 widely spaced geographic localities in Texas, Louisiana, and Florida and compared them to each other and to homologous sequences from 4 Mexican ducks (A. diazi), 13 American black ducks (A. rubripes), and 10 mallards (A. platyrhynchos). We identified 57 unique haplotypes composed of 665 or 666 nucleotides in the 246 control region sequences. Of the 665 homologous positions, 8.3% (n = 55) vary among haplotypes, and 98.2% (n = 54) of these occur within the first 351 nucleotides from the 5’ end of the outgroup sequence. Neighbor-joining analysis shows a large distal clade (52.5% of mottled ducks sampled in our study) composed of two reciprocally monophyletic clades of mottled duck haplotypes, one of which is endemic to Texas and Louisiana and the other endemic to Florida. No mottled ducks sampled in Florida occur in the clade composed of mottled ducks from Texas and Louisiana or vice versa, suggesting that (1) an enduring geographic split has existed for many years between east and west, and (2) gene flow currently is non-existent (or at least undetectable) across the central Gulf Coast. The remaining 47.5% of mottled ducks sampled in our study branch basally from this derived clade, show substantially less hierarchical structure, and fall into various lineage groups of mixed species composition with no geographic or species-specific pattern. Pairwise FST values corroborate the pattern of strong differentiation observed between Texas/Louisiana and Florida. Our findings are consistent with a pattern of partial lineage sorting from a polymorphic ancestral gene pool reshuffled by hybridizing mallards. Control region data and patterns of divergence in mallard-like species worldwide, furthermore, suggest that mottled ducks are close relatives of Mexican ducks, and in turn nested within black ducks. Genetic similarities to nominate mallards are less likely to be the product of common ancestry, but the result of past hybridization with a dichromatic mallard ancestor that invaded North America from Asia many generations ago. Our findings have several important consequences for the conservation biology of mottled ducks across the Gulf Coast and our understanding of the phylogeography of mallard-like species worldwide.

The unwitting inclusion of convergent characters in phylogenetic estimates poses a serious problem for efforts to recover phylogeny. Convergence is not inscrutable, however, particularly when one group of characters tracks phylogeny and another set tracks adaptive history. In such cases, convergent characters may be correlated with one or a few functional anatomical units and readily identifiable by using comparative methods. Stifftail ducks (Oxyurinae) offer one such opportunity to study correlated character evolution and function in the context of phylogenetic reconstruction. Morphological analyses place stifftail ducks as part of a large clade of diving ducks that includes the sea ducks (Mergini), Hymenolaimus, Merganetta, and Tachyeres, and possibly the pochards (Aythyini). Molecular analyses, on the other hand, place stifftails far from other diving ducks and suggest, moreover, that stifftails are polyphyletic. Mitochondrial cytochrome b gene sequences of eight stifftail species traditionally supposed to form a clade were compared with each other and with sequences from 50 other anseriform and galliform species. Stifftail ducks are not the sister group of sea ducks but lie outside the typical ducks (Anatinae). Of the four traditional stifftail genera, monophyly of Oxyura and its sister group relationship with Nomonyx are strongly supported. Heteronetta probably is the sister group of that clade, but support is weak. Biziura is not a true stifftail. Within Oxyura, Old World species (O. australis, O. leucocephala, O. maccoa) appear to form a clade, with New World species (O. jamaicensis, O. vittata) branching basally. Incongruence between molecules and morphology is interpreted to be the result of adaptive specialization and functional convergence in the hind limbs of Biziura and true stifftails. When morphological characters are divided into classes, only hind-limb characters are significantly in conflict with the molecular tree. Likewise, null models of synonymous and nonsynonymous substitution based on patterns of codon-degeneracy and chemical dissimilarity indicate that the nucleotide and amino acid changes postulated by the molecular tree are more plausible than those postulated by the morphological tree. These findings teach general lessons about the utility of highly adaptive characters (in particular those related to foraging ecology) and underscore the problems that convergence can pose for attempts to recover phylogeny. They also demonstrate how the concept of natural data partitions and simple models of evolution (e.g., parsimony, likelihood, neutrality) can be used to test the accuracy of independent phylogenetic estimates and provide arguments in favor of one tree topology over another.

Five species of migratory thrushes (Turdidae) occupy a transcontinental distribution across northern North America. They have largely overlapping breeding ranges, relatively similar ecological niches, and mutualistic relationships with northern woodland communities as insectivores and seed-dispersing frugivores. As an assemblage of ecologically similar species, and given other vertebrate studies, we predicted a shared pattern of genetic divergence among these species between their eastern and western populations, and also that the timing of the coalescent events might be similar and coincident with historical glacial events. To determine how these five lineages effectively established transcontinental distributions, we used mitochondrial cytochrome b sequences to assess genetic structure and lineage coalescence from populations on each side of the continent. Two general patterns occur. Hermit and Swainson’s thrushes (Catharus guttatus and C. ustulatus) have relatively deep divergences between eastern and western phylogroups, probably reflecting shared historic vicariance. The Veery (C. fuscescens), Gray-cheeked Thrush (C. minimus), and American Robin (Turdus migratorius) have relatively shallow divergences between eastern and western populations. However, coalescent and approximate Bayesian computational analyses indicated that among all species as many as five transcontinental divergence events occurred. Divergence within both Hermit and Swainson’s thrushes resembled the divergence between Gray-cheeked Thrushes and Veeries and probably occurred during a similar time period. Despite these species’ ecological similarities, the assemblage exhibits heterogeneity at the species level in how they came to occupy transcontinental northern North America but two general continental patterns at an among-species organizational level, likely related to lineage age.

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. © 2013 Blackwell Publishing Ltd.

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.

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

Hypoxia is one of the most important factors affecting survival at high altitude, and the major hemoglobin protein is a likely target of selection. We compared population genetic structure in the alphaA and betaA hemoglobin subunits (HBA2 and HBB) of five paired lowland and highland populations of Andean dabbling ducks to unlinked reference loci. In the hemoglobin genes, parallel amino acid replacements were overrepresented in highland lineages, and one to five derived substitutions occurred at external solvent-accessible positions on the alpha and beta subunits, at alpha(1)beta(1) intersubunit contacts, or in close proximity to inositolpentaphosphate (IPP) binding sites. Coalescent analyses incorporating the stochasticity of drift and mutation indicated that hemoglobin alleles were less likely to be transferred between highland and lowland populations than unlinked alleles at five other loci. Amino acid replacements that were overrepresented in the highlands were rarely found within lowland populations, suggesting that alleles segregating at high frequency in the highlands may be maladaptive in the lowlands and vice versa. Most highland populations are probably nonmigratory and locally adapted to the Altiplano, but gene flow for several species may be sufficiently high to retard divergence at unlinked loci. Heterozygosity was elevated in the alphaA or betaA subunits of highland populations exhibiting high gene flow between the southern lowlands and the highlands and in highland species that disperse seasonally downslope to midelevation environments from the central Andean plateau. However, elevated heterozygosity occurred more frequently in the alphaA subunit but not simultaneously in both subunits, suggesting that selection may be more constrained by epistasis in the betaA subunit. Concordant patterns among multiple species with different evolutionary histories and depths of historical divergence and gene flow suggest that the major hemoglobin genes of these five dabbling duck species have evolved adaptively in response to high-altitude hypoxia in the Andes.

Theory predicts that parallel evolution should be common when the number of beneficial mutations is limited by selective constraints on protein structure. However, confirmation is scarce in natural populations. Here we studied the major haemoglobin genes of eight Andean duck lineages and compared them to 115 other waterfowl species, including the bar-headed goose (Anser indicus) and Abyssinian blue-winged goose (Cyanochen cyanopterus), two additional species living at high altitude. One to five amino acid replacements were significantly overrepresented or derived in each highland population, and parallel substitutions were more common than in simulated sequences evolved under a neutral model. Two substitutions evolved in parallel in the alpha A subunit of two (Ala-alpha 8) and five (Thr-alpha 77) taxa, and five identical beta A subunit substitutions were observed in two (Ser-beta 4, Glu-beta 94, Met-beta 133) or three (Ser-beta 13, Ser-beta 116) taxa. Substitutions at adjacent sites within the same functional protein region were also observed. Five such replacements were in exterior, solvent-accessible positions on the A helix and AB corner of the alpha A subunit. Five others were in close proximity to inositolpentaphosphate binding sites, and two pairs of independent replacements occurred at two different alpha(1)beta(1) intersubunit contacts. More than half of the substitutions in highland lineages resulted in the acquisition of serine or threonine (18 gains vs. 2 losses), both of which possess a hydroxyl group that can hydrogen bond to a variety of polar substrates. The patterns of parallel evolution observed in these waterfowl suggest that adaptation to high-altitude hypoxia has resulted from selection on unique but overlapping sets of one to five amino acid substitutions in each lineage.

When populations become locally adapted to contrasting environments, alleles that have high fitness in only one environment may be quickly eliminated in populations adapted to other environments, such that gene flow is partly restricted. The stronger the selection, the more rapidly immigrant alleles of lower fitness will be eliminated from the population. However, gene flow may continue to occur at unlinked loci, and adaptive divergence can proceed in the face of countervailing gene flow if selection is strong relative to migration (s > m). We studied the population genetics of the major hemoglobin genes in yellow-billed pintails (Anas georgica) experiencing contrasting partial pressures of oxygen in the Andes of South America. High gene flow and weak population subdivision were evident at seven putatively neutral loci in different chromosomal linkage groups. In contrast, amino acid replacements (Ser-beta13, Ser-beta116, and Met-beta133) in the betaA hemoglobin subunit segregated by elevation between lowland and highland populations with significantly elevated F(ST). Migration rates for the betaA subunit alleles were approximately 17-24 times smaller than for five unlinked reference loci, the alphaA hemoglobin subunit (which lacks amino acid replacements) and the mitochondrial DNA control region. The betaA subunit alleles of yellow-billed pintails were half as likely to be transferred downslope, from the highlands to the lowlands, than in the opposite direction upslope. We hypothesize that migration between the lowlands and highlands is restricted by local adaptation, and the betaA hemoglobin subunit is a likely target of selection related to high-altitude hypoxia; however, gene flow may be sufficiently high to retard divergence at most unlinked loci. Individuals homozygous for lowland alleles may have relatively little difficulty dispersing to the highlands initially but may experience long-term fitness reduction. Individuals homozygous for highland genotypes, in contrast, would be predicted to have difficulty dispersing to the lowlands if their hemoglobin alleles confer high oxygen affinity, predicted to result in chronic erythrocytosis at low elevation. Heterozygous individuals may have a dispersal advantage if their hemoglobin has a wider range of function due to the presence of multiple protein isoforms with a mixture of different oxygen affinities.

We evaluated the potential effects of homoplasy, ancestral polymorphism, and hybridization as obstacles to resolving phylogenetic relationships within Nomonyx-Oxyura stiff-tailed ducks (Oxyurinae; subtribe Oxyurina). Mitochondrial DNA (mtDNA) control region sequences from 94 individuals supported monophyly of mtDNA haplotypes for each of the six species and provided no evidence of extant incomplete lineage sorting or inter-specific hybridization. The ruddy ducks (O. j. jamaicensis,O. j. andina, O. j. ferruginea) are each others' closest relatives, but the lack of shared haplotypes between O. j. jamaicensis and O. j. ferruginea suggests long-standing historical isolation. In contrast, O. j. andina shares haplotypes with O. j. jamaicensis and O. j. ferruginea, which supports Todd's (1979) and Fjeldsa's (1986) hypothesis that O. j. andina is an intergrade or hybrid subspecies of O. j. jamaicensis and O. j. ferruginea. Control region data and a much larger data set composed of approximately 8800 base pairs of mitochondrial and nuclear sequence for each species indicate that the two New World species, O. vittata and O. jamaicensis, branch basally within Oxyura. A clade of three Old World species (O. australis, O. maccoa, O. leucocephala) is well supported, but different loci and also different characters within the mtDNA data support three different resolutions of the Old World clade, yielding an essentially unresolved trichotomy. Fundamentally different factors limited the resolution of the mtDNA and nuclear gene trees. Gene trees for most nuclear loci were unresolved due to slow rates of mutation and a lack of informative variation, whereas uncertain resolution of the mtDNA gene tree was due to homoplasy. Within the mtDNA, approximately equal numbers of characters supported each of three possible resolutions. Parametric and nonparametric bootstrap analyses suggest that resolution of the mtDNA tree based on ~4300 bp per taxon is uncertain but that complete mtDNA sequences would yield a fully resolved gene tree. A short internode separating O. leucocephala from (O. australis, O. maccoa) in the best mtDNA tree combined with long terminal branches and substantial rate variation among nucleotide sites allowed the small number of changes occurring on the internode to be obscured by homoplasy in a significant portion of simulated data sets. Although most nuclear loci were uninformative, two loci supported a resolution of the Old World clade (O. maccoa, O. leucocephala) that is incongruent with the best mtDNA tree. Thus, incongruence between nuclear and mtDNA trees may be due to random sorting of ancestral lineages during the short internode, homoplasy in the mtDNA data, or both. The Oxyura trichotomy represents a difficult though likely common problem in molecular systematics. Given a short internode, the mtDNA tree has a greater chance of being congruent with the history of speciation because its effective population size (N(e)) is one-quarter that of any nuclear locus, but its resolution is more likely to be obscured by homoplasy. In contrast, gene trees for more slowly evolving nuclear loci will be difficult to resolve due to a lack of substitutions during the internode, and when resolved are more likely to be incongruent with the species history due to the stochastic effects of lineage sorting. We suggest that researchers consider first whether independent gene trees are adequately resolved and then whether those trees are congruent with the species history. In the case of Oxyura, the answer to both questions may be no. Complete mtDNA sequences combined with data from a very large number of nuclear loci may be the only way to resolve such trichotomies.

Species-level DNA phylogenies frequently suffer from two shortcomings--gene trees usually are constructed from a single locus, and often species are represented by only one individual. To evaluate the effect of these two shortcomings, we tested phylogenetic hypotheses within the wigeons and allies, a clade of Anas ducks (Anatidae) composed of five species. We sequenced two nuclear introns from the Z-chromosome-linked chromo-helicase binding protein gene (CHD1Zb and CHD1Za) and the mitochondrial DNA (mtDNA) control region for multiple individuals sampled from widespread geographic locations. We compared these phylogenies to previously published phylogenies constructed from morphology and protein coding regions of mtDNA. Relative to other nuclear introns, CHD showed remarkable phylogenetic utility. Of the 26 CHD1Zb alleles identified, only one was shared between two species, and the combined CHD datasets revealed that four of the five species were consistent with monophyly. Several species shared mtDNA haplotypes, which probably was a result of interspecific hybridization. Overall, the nuclear CHD tree and the mtDNA tree were more congruent with coding regions of mtDNA than they were with morphology.

Mitochondrial cytochrome b sequence data from 15 species of herons (Aves: Ardeidae), representing 13 genera, were compared with DNA hybridization data of single-copy nuclear DNA (scnDNA) from the same species in a taxonomic congruence assessment of heron phylogeny. The two data sets produced a partially resolved, completely congruent estimate of phylogeny with the following basic structure: (Tigrisoma, Cochlearius, (((Zebrilus, (Ixobrychus, Botaurus)), (((Ardea, Casmerodius), Bubulcus), ((Egretta thula, Egretta caerulea, Egretta tricolor), Syrigma), Butorides, Nycticorax, Nyctanassa)))). Because congruence indicated similar phylogenetic information in the two data sets, we used the relatively unsaturated DNA hybridization distances as surrogates of time to examine graphically the patterns and rates of change in cytochrome b distances. Cytochrome b distances were computed either from whole sequences or from partitioned sequences consisting of transitions, transversions, specific codon site positions, or specific protein-coding regions. These graphical comparisons indicated that unpartitioned cytochrome b has evolved at 5-10 times the rate of scnDNA. Third-position transversions appeared to offer the most useful sequence partition for phylogenetic analysis because of their relatively fast rate of substitution (two times that of scnDNA) and negligible saturation. We also examined lineage-based rates of evolution by comparing branch length patterns between the nuclear and cytochrome b trees. The degree of correlation in corresponding branch lengths between cytochrome b and DNA hybridization trees depended on DNA sequence partitioning. When cytochrome b sequences were not partitioned, branch lengths in the cytochrome b and DNA hybridization trees were not correlated. However, when cytochrome b sequences were reduced to third-position transversions (i.e., unsaturated, relatively fast changing data), branch lengths were correlated. This finding suggests that lineage-based rates of DNA evolution in nuclear and mitochondrial genomes are influenced by common causes.

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.