Search

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.

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.

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.

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.

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.

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.

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.

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.