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.     

Postcatastrophe Population Dynamics and Density Dependence of an Endemic Island Duck

ABSTRACT Laysan ducks (Anas laysanensis) are restricted to approximately 9 km² in the Northwestern Hawaiian Islands, USA. To evaluate the importance of density dependence for Laysan ducks, we conducted a Bayesian analysis to estimate the parameters of a Gompertz model and the magnitude of process variation and observation error based on the fluctuations in Laysan duck abundance on Laysan Island from 1994 to 2007. This model described a stationary distribution for the population at carrying capacity that fluctuates around a long-term mean of 456 ducks and is between 316 to 636 ducks 95% of the time. This range of expected variability can be used to identify changes in population size that warn of catastrophic events. Density-dependent population dynamics may explain the recovery of Laysan duck from catastrophic declines and allow managers to identify population monitoring thresholds.     

Persistence of an endangered native duck,feral mallards,and multiple hybrid swarms across the main Hawaiian Islands

Interspecific hybridization is recognized as an important process in the evolutionary dynamics of both speciation and the reversal of speciation. However, our understanding of the spatial and temporal patterns of hybridization that erode versus promote species boundaries is incomplete. The endangered, endemic koloa maoli (or Hawaiian duck, Anas wyvilliana) is thought to be threatened with genetic extinction through ongoing hybridization with an introduced congener, the feral mallard (A. platyrhynchos). We investigated spatial and temporal variation in hybrid prevalence in populations throughout the main Hawaiian Islands, using genomic data to characterize population structure of koloa, quantify the extent of hybridization, and compare hybrid proportions over time. To accomplish this, we genotyped 3,308 double‐digest restriction‐site‐associated DNA (ddRAD) loci in 425 putative koloa, mallards, and hybrids from populations across the main Hawaiian Islands. We found that despite a population decline in the last century, koloa genetic diversity is high. There were few hybrids on the island of Kaua?i, home to the largest population of koloa. By contrast, we report that sampled populations outside of Kaua?i can now be characterized as hybrid swarms, in that all individuals sampled were of mixed koloa × mallard ancestry. Further, there is some evidence that these swarms are stable over time. These findings demonstrate spatial variation in the extent and consequences of interspecific hybridization, and highlight how islands or island‐like systems with small population sizes may be especially prone to genetic extinction when met with a congener that is not reproductively isolated.     

Effects of Asymmetric Nuclear Introgression,Introgressive Mitochondrial Sweep,and Purifying Selection on Phylogenetic Reconstruction and Divergence Estimates in the Pacific Clade of Locustella Warblers

When isolated but reproductively compatible populations expand geographically and meet, simulations predict asymmetric introgression of neutral loci from a local to invading taxon. Genetic introgression may affect phylogenetic reconstruction by obscuring topology and divergence estimates. We combined phylogenetic analysis of sequences from one mtDNA and 12 nuDNA loci with analysis of gene flow among 5 species of Pacific Locustella warblers to test for presence of genetic introgression and its effects on tree topology and divergence estimates. Our data showed that nuDNA introgression was substantial and asymmetrical among all members of superspecies groups whereas mtDNA showed no introgression except a single species pair where the invader''s mtDNA was swept by mtDNA of the local species. This introgressive sweep of mtDNA had the opposite direction of the nuDNA introgression and resulted in the paraphyly of the local species'' mtDNA haplotypes with respect to those of the invader. Тhe multilocus nuDNA species tree resolved all inter- and intraspecific relationships despite substantial introgression. However, the node ages on the species tree may be underestimated as suggested by the differences in node age estimates based on non-introgressing mtDNA and introgressing nuDNA. In turn, the introgressive sweep and strong purifying selection appear to elongate internal branches in the mtDNA gene tree.     

Phylogeography of the mallard Anas platyrhynchos from Eurasia inferred from sequencing of the mtDNA control region

Phylogenetic relationships, demographic history, and geographic distribution of the mtDNA haplotypes of the mallard Anas platyrhynchos were examined in three populations, Indian, Northern European, and Far Eastern. Two divergent halotype groups, A and B, were found in the Far Eastern population, while haplotypes identified in Northern European and Far Eastern populations were exclusively of the A group. The presence of B group haplotypes in the Far Eastern population can be explained either in terms of hybridization of the mallard with spot-billed duck Anas zonorhyncha at the south of the Russian Far East, or by the mtDNA paraphyly in mallards. In general, mallards from Eurasia were characterized by low genetic population differentiation along with slightly expressed phylogeographic structure. The most differentiated was the population from India (??_st = 0.076?0.077), while the difference between Northern European and Far Eastern populations was extremely low (??_st = 0.0029). Differentiation of Anas platyrhynchos Indian population was determined by the fact that a part of the population, inhabiting southern and eastern coasts of the Hindustan Peninsula, was resident.     

Postfledging Survival of Laysan Ducks

ABSTRACT Precise and unbiased estimates of demographic parameters are necessary for effective population monitoring and to parameterize population models (e.g., population viability analyses). This is especially important for endangered species, where recovery planning and managers' decisions can influence species persistence. In this study, we used mark—recapture methods to estimate survival of fledged juveniles (hatch-yr [HY]) and adult (after-hatch-yr [AHY]) Laysan ducks (Anas laysanensis), an endangered anatid restricted to Laysan Island in the northwestern Hawaiian Islands. To better understand population dynamics, we examined how survival varied as a function of Laysan duck density during 1998–2004. Using random effects models, we also quantified process variation in survival, thereby quantifying the appropriate source of variation for future population models. The dataset supported variation in survival that was time (yr), age (AHY vs. HY), and sex specific. Due to small sample sizes, we did not examine time specificity in the survival of HY ducks. Survival of HY ducks was 0.832 (SE = 0.087) for females (n = 21) and 0.999 (SE < 0.001) for males (n = 15) during 1998–2001. Trends in time and density lacked support as sources of variation in the survival of AHY ducks during 1998–2004. After-hatch-year survival ranged from 0.792 (SE = 0.033) to 0.999 (SE < 0.001). Where we modeled survival as a random effect, annual survival for AHY females was 0.881 (SE = 0.017) and process variation (σ_s) was 0.034. For AHY males, annual survival (μ_s) was 0.906 (SE = 0.019) and process variation (σ_s) was 0.040. This information will improve existing population viability analysis models for Laysan ducks. We believe that monitoring the source and translocation populations will be paramount for increasing our understanding of Laysan duck dynamics, recovery planning, and population management.     

Identification of endangered Hawaiian ducks (Anas wyvilliana), introduced North American mallards (A. platyrhynchos) and their hybrids using multilocus genotypes

Hawaiian ducks (Anas wyvilliana), or koloa, are endemic to the Hawaiian Islands and are listed as a federal and state endangered species. Hybridization between koloa and introduced mallards (A. platyrhynchos) is believed to be a primary threat to the recovery of koloa. We evaluated the utility of two sets of nuclear markers (microsatellite loci and amplified fragment length polymorphisms) and a variable portion of the mitochondrial DNA control region to distinguish among koloa, mallards, and hybrids. We show that microsatellite and AFLP markers can be used to distinguish between koloa and mallard-koloa hybrids with a high degree of confidence. For all but one of the putative koloa in our sample, the posterior probability of belonging to the koloa category was >0.90. Similarly all but one of the mallard-koloa hybrids were assigned to the hybrid category with posterior probabilities >0.98. Subsets of markers led to poorer resolution among koloa, mallard and hybrid categories. Among a sample of 61 koloa, hybrids and mallards, we found 25 different mtDNA haplotypes, belonging to two groups of haplotypes (A and B) identified previously in mallards and their relatives. All putative koloa samples exhibited group B haplotypes, of which 65% comprised one haplotype, while the rest were divided among four haplotypes. All Hawai’i mallard samples exhibited haplotypes that belonged to group A. Hybrids and California mallards exhibited haplotypes belonging to both groups, but a majority were of group A, suggesting that hybridization may more commonly involve mating between Hawai’i mallard females and koloa males.     

A new tool for studying waterfowl immune and metabolic responses: Molecular level analysis using kinome profiling

Here, we describe the design of an Anas‐specific kinome peptide array that can be used to study the immunometabolic responses of mallard and American black duck to pathogens, contaminants, and environmental stress. The peptide arrays contain 2,642 unique phosphorylate‐able peptide sequences representing 1,900 proteins. These proteins cover a wide array of metabolic and immunological processes, and 758 Gene Ontology Biological processes are statistically significantly represented on the duck peptide array of those 164 contain the term “metabolic” and 25 “immune.” In addition, we conducted a comparison of mallard to American black duck at a genetic and proteomic level. Our results show a significant genomic and proteomic overlap between these two duck species, so that we have designed a cross‐reactive peptide array capable of studying both species. This is the first reported development of a wildlife species‐specific kinome peptide array.     

Evolutionary basis of mitonuclear discordance between sister species of mole salamanders (Ambystoma sp.)

Distinct genetic markers should show similar patterns of differentiation between species reflecting their common evolutionary histories, yet there are increasing examples of differences in the biogeographic distribution of species‐specific nuclear (nuDNA) and mitochondrial DNA (mtDNA) variants within and between species. Identifying the evolutionary processes that underlie these anomalous patterns of genetic differentiation is an important goal. Here, we analyse the putative mitonuclear discordance observed between sister species of mole salamanders (Ambystoma barbouri and A. texanum) in which A. barbouri‐specific mtDNA is found in animals located within the range of A. texanum. We test three hypotheses for this discordance (undetected range expansion, mtDNA introgression, and hybridization) using nuDNA and mtDNA data analysed with methods that varied in the parameters estimated and the timescales measured. Results from a Bayesian clustering technique (structure ), bidirectional estimates of gene flow (migrate ‐n and IMa2) and phylogeny‐based methods (*beast , buck y) all support the conclusion that the discordance is due to geographically restricted mtDNA introgression from A. barbouri into A. texanum. Limited data on species‐specific tooth morphology match this conclusion. Significant differences in environmental conditions exist between sites where A. texanum with and without A. barbouri‐like mtDNA occur, suggesting a possible role for selection in the process of introgression. Overall, our study provides a general example of the value of using complimentary analyses to make inferences of the directionality, timescale, and source of mtDNA introgression in animals.     

Population Genetics of a Translocated Population of Mottled Ducks and Allies

Translocating species is an important management tool to establish or expand the range of species. Success of translocations requires an understanding of potential consequences, including whether a sufficient number of individuals were used to minimize founder effects and if interspecific hybridization poses a threat. We provide an updated and comprehensive genetic assessment of a 1970s–1980s translocation and now established mottled duck (Anas fulvigula) population in South Carolina, USA. In addition to examining the population genetics of these mottled ducks, we simulated expected genetic assignments for generational hybrids (F1–F10), permitting formal purity assignment across samples to identify true hybrids and establish hybridization rates. In addition to wild mallards (A. platyrhynchos), we tested for presence of hybrids with migrant American black ducks (A. rubripes) and released domestic game-farm mallards (A. p. domesticus). We used wild reference populations of North American mallard-like ducks and sampled game-farm mallards from 2 sites in South Carolina that could potentially interbreed with mottled ducks. Despite 2 different subspecies of mottled duck (Florida [A. f. fulvigula] and the Western Gulf Coast [A. f. maculatlus]) used in original translocations, we determined the gene pool of the Western Gulf Coast mottled duck was overwhelmingly represented in South Carolina's current population. We found no evidence of founder effects or inbreeding and concluded the original translocation of 1,285 mottled ducks was sufficient to maintain current genetic diversity. We identified 7 hybrids, including an F1 and 3 late-staged (i.e., F2–F3 backcrosses) mottled duck × black duck hybrids, 1 F2-mottled duck backcrossed with a wild mallard, and 2 F3-mottled ducks introgressed with game-farm mallard. We estimated a 15% hybridization rate in our mottled duck dataset; however, the general lack of F1 and intermediate hybrids were inconsistent with scenarios of high hybridization rates or presence of a hybrid swarm. Instead, our results suggested a scenario of infrequent interspecific hybridization between South Carolina's mottled ducks and congeners. We concluded that South Carolina's mottled duck population is sufficiently large now to absorb current hybridization rates because 85% of sampled mottled ducks were pure. These results demonstrate the importance in managing and maintaining large parental populations to counter hybridization. As such, future population management of mottled ducks in South Carolina will benefit from increased geographical and continued sampling to monitor hybridization rates with closely related congeners. We also suggest that any future translocations of mottled ducks to coastal South Carolina should originate from the Western Gulf Coast. © 2021 The Wildlife Society.     

Multilocus phylogeography of Australian teals (Anas spp.): a case study of the relationship between vagility and genetic structure

Biogeographic barriers potentially restrict gene flow but variation in dispersal or vagility can influence the effectiveness of these barriers among different species and produce characteristic patterns of population genetic structure. The objective of this study was to investigate interspecific and intraspecific genetic structure in two closely related species that differ in several life‐history characteristics. The grey teal Anas gracilis is geographically widespread throughout Australia with a distribution that crosses several recognized biogeographic barriers. This species has high vagility as its extensive movements track broad‐scale patterns in rainfall. In contrast, the closely related chestnut teal A. castanea is endemic to the mesic southeastern and southwestern regions of Australia and is more sedentary. We hypothesized that these differences in life‐history characteristics would result in more pronounced population structuring in the chestnut teal. We sequenced five nuclear loci (nuDNA) for 49 grey teal and 23 chestnut teal and compared results to published mitochondrial DNA (mtDNA) sequences. We used analysis of molecular variance to examine population structure, and applied coalescent based approaches to estimate demographic parameters. As predicted, chestnut teal were more strongly structured at both mtDNA and nuDNA (Φ_ST= 0.163 and 0.054, respectively) than were grey teal (Φ_ST < 0.0001 for both sets of loci). Surprisingly, a greater proportion of the total genetic variation was partitioned among populations within species (Φ_SC= 0.014 and 0.047 for nuDNA and mtDNA, respectively) than between the two species (Φ_CT < 0.0001 for both loci). The ‘Isolation with Migration’ coalescent model suggested a late Pleistocene divergence between the taxa, but remarkably, a deeper divergence between the southeastern and southwestern populations of chestnut teal. We conclude that dispersal potential played a prominent role in the structuring of populations within these species and that divergent selection associated with ecology and life history traits likely contributed to rapid and recent speciation in this pair.     

Assessing the fidelity of ancient DNA sequences amplified from nuclear genes

To date, the field of ancient DNA has relied almost exclusively on mitochondrial DNA (mtDNA) sequences. However, a number of recent studies have reported the successful recovery of ancient nuclear DNA (nuDNA) sequences, thereby allowing the characterization of genetic loci directly involved in phenotypic traits of extinct taxa. It is well documented that postmortem damage in ancient mtDNA can lead to the generation of artifactual sequences. However, as yet no one has thoroughly investigated the damage spectrum in ancient nuDNA. By comparing clone sequences from 23 fossil specimens, recovered from environments ranging from permafrost to desert, we demonstrate the presence of miscoding lesion damage in both the mtDNA and nuDNA, resulting in insertion of erroneous bases during amplification. Interestingly, no significant differences in the frequency of miscoding lesion damage are recorded between mtDNA and nuDNA despite great differences in cellular copy numbers. For both mtDNA and nuDNA, we find significant positive correlations between total sequence heterogeneity and the rates of type 1 transitions (adenine --> guanine and thymine --> cytosine) and type 2 transitions (cytosine --> thymine and guanine --> adenine), respectively. Type 2 transitions are by far the most dominant and increase relative to those of type 1 with damage load. The results suggest that the deamination of cytosine (and 5-methyl cytosine) to uracil (and thymine) is the main cause of miscoding lesions in both ancient mtDNA and nuDNA sequences. We argue that the problems presented by postmortem damage, as well as problems with contamination from exogenous sources of conserved nuclear genes, allelic variation, and the reliance on single nucleotide polymorphisms, call for great caution in studies relying on ancient nuDNA sequences.     

Mallard–Black Duck Hybridization and Population Genetic Structure in North Carolina

North Carolina, USA, represents the southern extent of the American black duck's (Anas rubripes) breeding range. Mallards (A. platyrhynchos) are present on the breeding grounds of the American black duck and hybridization is observed between these species; therefore, we assessed the genetic integrity, hybridization rates, and population structure of this local breeding population. We extracted genomic and mitochondrial DNA from chorioallantoic membranes and contour feathers from monitored black duck nests. We then prepared the extracted DNA for analysis using high-throughput DNA sequencing methods (ddRAD-seq). First, we assessed nuclear and mitochondrial population structure, genetic diversity, and differentiation across samples from North Carolina, and compared them against 199 genetically vetted mallards, black ducks, and mallard × black duck hybrids that served as genetic references. Next, we tested for parentage and sibling relationship and overall relatedness of black ducks in North Carolina. We recovered strong population structure and high co-ancestry across genetic markers due to interrelatedness among sampled nests in North Carolina and concluded that black ducks have been locally breeding in this area for a prolonged period of time. Despite a high level of interrelatedness among our samples, nucleotide diversity was similar to the reference continental black duck population, suggesting little effect of genetic drift, including inbreeding. Additionally, we conclude that molecular diversity of black ducks in North Carolina is maintained at reference population levels through the influx of genetic material from unrelated, migrating male black ducks. Finally, we report a hybridization level of 47.5%, covering 3 filial generations. Of identified hybrids, 54.7% and 53% were the direct result of interbreeding between black ducks and captive-reared or wild mallards, respectively. We conclude that because of high rates of interspecific hybridization and successive backcrossing events, introgression from wild and feral mallards is occurring into this population of breeding black ducks and requires careful consideration in future management efforts. © 2021 The Wildlife Society.     

Mottled ducks,feral mallards,and their hybrids in Florida

The nonmigratory and endemic Florida mottled duck (Anas fulvigula fulvigula) is facing conservation threats from the combined effects of urbanization and introgressive hybridization with feral mallards (Anas platyrhynchos) and mallard x mottled duck hybrids. In the past, the status of the Florida mottled duck population was assessed during annual aerial surveys and most brown ducks (mottled ducks, mallards, and hybrids of them) detected during the survey would have been mottled ducks. But the release of domesticated mallards for aesthetic purposes has led to increases in the prevalence of mallards-hybrids (mallards or mallard x mottled duck hybrids) throughout peninsular Florida, USA, and because it is impossible to differentiate among mottled ducks, female mallards, and hybrids during aerial surveys, helicopter surveys were halted in 2009 until state researchers could conduct a range-wide study to determine what proportion of brown ducks are mottled ducks versus mallards-hybrids. We used plumage keys and high-resolution photography to categorize brown ducks from 557 wetland grid points as either mottled ducks or mallards-hybrids. Of the 5,179 brown ducks categorized, 40.1% were mottled ducks and 59.9% were mallards-hybrids. We used logistic regression analysis to model the interactive effect of a site's latitude and level of urbanization (urban gradient value within a 2-km buffer) to generate a predictive raster surface (1-km resolution) of the study area with values corresponding to the probability that a brown duck observed within a cell is a pure mottled duck. Predicted values will be used as correction factors when estimating final mottled duck population abundance from brown-duck survey data. Additionally, the predictive raster surface will be used to identify wetlands where mottled ducks remain predominant so that these sites can be targeted for preservation. Overall, mallards-hybrids outnumbered mottled ducks throughout most of peninsular Florida, especially in more urbanized regions, and their current prevalence rate presents a serious conservation threat, via hybridization, to extant mottled duck populations.     

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.     

SEQUENTIAL RADIATIONS AND PATTERNS OF SPECIATION IN THE HAWAIIAN CRICKET GENUS LAUPALA INFERRED FROM DNA SEQUENCES

The tremendous diversity of endemic Hawaiian crickets is thought to have originated primarily through intraisland radiations, in contrast to an interisland mode of diversification in the native Hawaiian Drosophila. The Hawaiian cricket genus Laupala (family Gryllidae) is one of several native genera of flightless crickets found in rain-forest habitat across the Hawaiian archipelago. I examined the phylogenetic relationships among mitochondrial DNA (mtDNA) sequences sampled from 17 species of Laupala, including the 12S ribosomal RNA (rRNA), transfer RNA (RNA)^val and 16S rRNA regions. The distribution of mtDNA variants suggests that species within Laupala are endemic to single islands. The phylogenetic estimate produced from both maximum likelihood and maximum parsimony supports the hypothesis that speciation in Laupala occurred mainly within islands. The inferred biogeographical history suggests that diversification in Laupala began on Kauai, the oldest rain-forested Hawaiian island. Subsequently, colonization to younger islands in the archipelago resulted in a radiation of considerable phylogenetic diversity. Phylogenetic patterns in mtDNA are not congruent with prior systematic or taxonomic hypotheses. Hypotheses that may explain the conflict between the phylogenetic patterns of mtDNA variation and the species taxonomy are discussed.     

Speciation genomics and a role for the Z chromosome in the early stages of divergence between Mexican ducks and mallards

Speciation is a continuous and dynamic process, and studying organisms during the early stages of this process can aid in identifying speciation mechanisms. The mallard (Anas platyrhynchos) and Mexican duck (A. [p.] diazi) are two recently diverged taxa with a history of hybridization and controversial taxonomy. To understand their evolutionary history, we conducted genomic scans to characterize patterns of genetic diversity and divergence across the mitochondrial DNA (mtDNA) control region, 3523 autosomal loci and 172 Z‐linked sex chromosome loci. Between the two taxa, Z‐linked loci (Φ_ST = 0.088) were 5.2 times more differentiated than autosomal DNA (Φ_ST = 0.017) but comparable to mtDNA (Φ_ST = 0.092). This elevated Z differentiation deviated from neutral expectations inferred from simulated data that incorporated demographic history and differences in effective population sizes between marker types. Furthermore, 3% of Z‐linked loci, compared to <0.1% of autosomal loci, were detected as outlier loci under divergent selection with elevated relative (Φ_ST) and absolute (d_XY) estimates of divergence. In contrast, the ratio of Z‐linked and autosomal differentiation among the seven Mexican duck sampling locations was close to 1:1 (Φ_ST = 0.018 for both markers). We conclude that between mallards and Mexican ducks, divergence at autosomal markers is largely neutral, whereas greater divergence on the Z chromosome (or some portions thereof) is likely the product of selection that has been important in speciation. Our results contribute to a growing body of literature indicating elevated divergence on the Z chromosome and its likely importance in avian speciation.     

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.     

Nuclear markers confirm taxonomic status and relationships among highly endangered and closely related right whale species

Right whales (genus: Eubalaena) are among the most endangered mammals, yet their taxonomy and phylogeny have been questioned. A phylogenetic hypothesis based on mitochondrial DNA (mtDNA) variation recently prompted a taxonomic revision, increasing the number of right whale species to three. We critically evaluated this hypothesis using sequence data from 13 nuclear DNA (nuDNA) loci as well as the mtDNA control region. Fixed diagnostic characters among the nuclear markers strongly support the hypothesis of three genetically distinct species, despite lack of any diagnostic morphological characters. A phylogenetics analysis of all data produced a strict consensus cladogram with strong support at nodes that define each right whale species as well as relationships among species. Results showed very little conflict among the individual partitions as well as congruence between the mtDNA and nuDNA datasets. These data clearly demonstrate the strength of using numerous independent genetic markers during a phylogenetics analysis of closely related species. In evaluating phylogenetic support contributed by individual loci, 11 of the 14 loci provided support for at least one of the nodes of interest to this study. Only a single marker (mtDNA control region) provided support at all four nodes. A study using any single nuclear marker would have failed to support the proposed phylogeny, and a strong phylogenetic hypothesis was only revealed by the simultaneous analysis of many nuclear loci. In addition, nu DNA and mtDNA data provided complementary levels of support at nodes of different evolutionary depth indicating that the combined use of mtDNA and nuDNA data is both practical and desirable.     

Cryptic and extensive hybridization between ancient lineages of American crows

Most species and therefore most hybrid zones have historically been defined using phenotypic characters. However, both speciation and hybridization can occur with negligible morphological differentiation. Recently developed genomic tools provide the means to better understand cryptic speciation and hybridization. The Northwestern Crow (Corvus caurinus) and American Crow (Corvus brachyrhynchos) are continuously distributed sister taxa that lack reliable traditional characters for identification. In this first population genomic study of Northwestern and American crows, we use genomic SNPs (nuDNA) and mtDNA to investigate the degree of genetic differentiation between these crows and the extent to which they may hybridize. Our results indicate that American and Northwestern crows have distinct evolutionary histories, supported by two nuDNA ancestry clusters and two 1.1%‐divergent mtDNA clades dating to the late Pleistocene, when glacial advances may have isolated crow populations in separate refugia. We document extensive hybridization, with geographic overlap of mtDNA clades and admixture of nuDNA across >900 km of western Washington and western British Columbia. This broad hybrid zone consists of late‐generation hybrids and backcrosses, but not recent (e.g., F1) hybrids. Nuclear DNA and mtDNA clines had concordant widths and were both centred in southwestern British Columbia, farther north than previously postulated. Overall, our results suggest a history of reticulate evolution in American and Northwestern crows, perhaps due to recurring neutral expansion(s) from Pleistocene glacial refugia followed by lineage fusion(s). However, we do not rule out a contributing role for more recent potential drivers of hybridization, such as expansion into human‐modified habitats.