Determining evolutionary origin and phylogenetic relationships of mallard-like ducks of Oceania, greater Indonesia, and the Philippines with ddRAD-seq data.

AIM: We aim to determine the evolutionary origins and population genetics of mallard-like ducks of Oceania, greater Indonesia, and the Philippines. LOCATION: Oceania, greater Indonesia, and the Philippines. TAXON: Mallard (Anas platyrhynchos), Pacific black duck (A. superciliosa spp.), and Philippine duck (A. luzonica) METHODS: Thousands of nuclear ddRAD-seq loci and the mitochondrial DNA control region were assayed across individuals representative of each species' range. We assessed population structure and phylogenetic relationships, as well as estimated demographic histories to reconstruct the biogeographical history of each species. RESULTS: Philippine and Pacific black ducks represent unique genetic lineages that diverged from the mallard 1-2 million years ago. We find no support for the Philippine duck representing a hybrid species as once posited; however, their low levels of genetic diversity requires further attention. We find a lack of substructure among Philippine ducks. However, we found pronounced differentiation between subspecies of Pacific black ducks, especially between A. s. superciliosa from New Zealand and A. s. rogersi from Australia, Papua New Guinea, and Timor-Leste, Indonesia. Anas superciliosa pelewensis gave mixed results; individuals from the Solomon Islands were differentiated from the other subspecies, but those from the island of Aunu'u, American Samoa, were genetically more similar to A. s. rogersi than A. s. pelewensis samples from the Solomon Islands. Finally, we find limited evidence of interspecific gene flow at evolutionary scales, and mallard introgression among contemporary samples. MAIN CONCLUSIONS: Mallard-like ducks radiated across Oceania, greater Indonesia, and the Philippines within the last 2 million years. Only the Pacific black duck showed unique sub-structuring that largely followed known sub-species ranges, except for A. s. pelewensis. We posit that the high interrelatedness among Solomon Island samples suggests that their genetic distinctiveness may simply be the result of high levels of genetic drift. In contrast, we conclude that mainland Australian Pacific black ducks were the most likely source for the recent colonization of American Samoa. As a result, our findings suggest that either the A. s. pelewensis subspecies designations and/or its geographical range may require re-evaluation. Continued re-evaluation of evolutionary and taxonomic relationships is necessary when attempting to reconstruct and understand biogeographical histories, with important implications towards any attempts to implement conservation strategies.

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

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

Population genomic analyses reveal population structure and major hubs of invasive Anopheles stephensi in the Horn of Africa.

Anopheles stephensi invasion in the Horn of Africa (HoA) poses a substantial risk of increased malaria disease burden in the region. An understanding of the history of introduction(s), establishment(s) and potential A. stephensi sources in the HoA is needed to predict future expansions and establish where they may be effectively controlled. To this end, we take a landscape genomic approach to assess A. stephensi origins and spread throughout the HoA, information essential for vector control. Specifically, we assayed 2070 genome-wide single nucleotide polymorphisms across 214 samples spanning 13 populations of A. stephensi from Ethiopia and Somaliland collected in 2018 and 2020, respectively. Principal component and genetic ancestry analyses revealed clustering that followed an isolation-by-distance pattern, with genetic divergence among the Ethiopian samples significantly correlating with geographical distance. Additionally, genetic relatedness was observed between the northeastern and east central Ethiopian A. stephensi populations and the Somaliland A. stephensi populations. These results reveal population differentiation and genetic connectivity within HoA A. stephensi populations. Furthermore, based on genetic network analysis, we uncovered that Dire Dawa, the site of a spring 2022 malaria outbreak, was one of the major hubs from which sequential founder events occurred in the rest of the eastern Ethiopian region. These findings can be useful for the selection of sites for heightened control to prevent future malaria outbreaks. Finally, we did not detect significant genotype-environmental associations, potentially due to the recency of their colonization and/or other anthropogenic factors leading to the initial spread and establishment of A. stephensi. Our study highlights how coupling genomic data at landscape levels can shed light into even ongoing invasions.

Evolutionary and ecological drivers of local adaptation and speciation in a North American avian species complex.

Throughout the speciation process, genomic divergence can be differentially impacted by selective pressures, as well as gene flow and genetic drift. Disentangling the effects of these evolutionary mechanisms remains challenging, especially for nonmodel organisms. Accounting for complex evolutionary histories and contemporary population structure often requires sufficient sample sizes, for which the expense of full genomes remains prohibitive. Here, we demonstrate the utility of partial-genome sequence data for range-wide samples to shed light into the divergence process of two closely related ducks, the Mexican duck (Anas diazi) and mallard (A. platyrhynchos). We determine the role of selective and neutral processes during speciation of Mexican ducks by integrating evolutionary and demographic modelling with genotype-environment and genotype-phenotype association testing. First, evolutionary models and demographic analyses support the hypothesis that Mexican ducks originally diverged ~300,000 years ago in climate refugia arising during a glacial period in southwest North America, and that subsequent environmental selective pressures played a key role in divergence. Mexican ducks then showed cyclical demographic patterns that probably reflected repeated range expansions and contractions, along with bouts of gene flow with mallards during glacial cycles. Finally, we provide evidence that sexual selection acted on several phenotypic traits as a co-evolutionary process, facilitating the development of reproductive barriers that initially arose due to strong ecological selection. More broadly, this work reveals that the genomic and phenotypic patterns observed across species complexes are the result of myriad factors that contribute in dynamic ways to the evolutionary trajectories of a lineage.

Phylogenomics reveals ancient and contemporary gene flow contributing to the evolutionary history of sea ducks (Tribe Mergini).

Insight into complex evolutionary histories continues to build through broad comparative phylogenomic and population genomic studies. In particular, there is a need to understand the extent and scale that gene flow contributes to standing genomic diversity and the role introgression has played in evolutionary processes such as hybrid speciation. Here, we investigate the evolutionary history of the Mergini tribe (sea ducks) by coupling multi-species comparisons with phylogenomic analyses of thousands of nuclear ddRAD-seq loci, including Z-sex chromosome and autosomal linked loci, and the mitogenome assayed across all extant sea duck species in North America. All sea duck species are strongly structured across all sampled marker types (pair-wise species ΦST > 0.2), with clear genetic assignments of individuals to their respective species, and phylogenetic relationships recapitulate known relationships. Despite strong species integrity, we identify at least 18 putative hybrids; with all but one being late generational backcrosses. Most interesting, we provide the first evidence that an ancestral gene flow event between long-tailed ducks (Clangula hyemalis) and true Eiders (Somateria spp.) not only moved genetic material into the former species, but likely generated a novel species - the Steller's eider (Polysticta stelleri) - via hybrid speciation. Despite generally low contemporary levels of gene flow, we conclude that hybridization has and continues to be an important process that shifts novel genetic variation between species within the tribe Mergini. Finally, we outline methods that permit researchers to contrast genomic patterns of contemporary versus ancestral gene flow when attempting to reconstruct potentially complex evolutionary histories.

Population genetics and conservation of recently discovered springsnails in Arizona.

Establishing baseline geographical distributions of extant genetic diversity is increasingly important for future conservation efforts of freshwater species. We analyse the mitochondrial cytochrome c oxidase subunit I (COI) gene to taxonomically characterize 233 samples from recently discovered springsnail populations throughout 17 sites in Arizona, USA. A total of 28 unique COI haplotypes were recovered, with the number of haplotypes ranging from 1 to 4 by population in Arizona. Phylogenetic analyses resulted in haplotypes from 13 of 17 locations in Arizona being successfully identified to species, with five described and three undescribed species in the genus Pyrgulopsis (Hydrobiidae). Future work will require in-depth morphological work to clarify the taxonomic status of these putatively novel species. Importantly, among recovered species, we identified haplotype diversity of the critically endangered Three Forks springsnail, Pyrgulopsis trivialis, in the eastern Gila watershed, which will inform wildlife managers in deciding which source populations to use in reintroduction efforts. We discuss possible causes for observed population structure of Arizona's springsnail populations, with suggestions for the future sampling schemes necessary for the conservation of this uniquely important freshwater mollusc.

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.

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.

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.

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

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

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

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

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

simRestore: A decision-making tool for adaptive management of the native genetic status of wild populations.

Anthropogenic hybridization, or higher and non-natural rates of gene flow directly and indirectly induced by human activities, is considered a significant threat to biodiversity. The primary concern for conservation is the potential for genomic extinction and loss of adaptiveness for native species due to the extensive introgression of non-native genes. To alleviate or reverse trends for such scenarios requires the direct integration of genomic data within a model framework for effective management. Towards this end, we developed the simRestore R program as a decision-making tool that integrates ecological and genomic information to simulate ancestry outcomes from optimized conservation strategies. In short, the program optimizes supplementation and removal strategies across generations until a set native genetic threshold is reached within the studied population. Importantly, in addition to helping with initial decision-making, simulations can be updated with the outcomes of ongoing efforts, allowing for the adaptive management of populations. After demonstrating functionality, we apply and optimize among actionable management strategies for the endangered Hawaiian duck for which the current primary threat is genetic extinction through ongoing anthropogenic hybridization with feral mallards. Simulations demonstrate that supplemental and removal efforts can be strategically tailored to move the genetic ancestry of Hawaii's hybrid populations towards Hawaiian duck without the need to completely start over. Further, we discuss ecological parameter sensitivity, including which factors are most important to ensure genetic outcomes (i.e. number of offspring). Finally, to facilitate use, the program is also available online as a Shiny Web application.

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

Evidence from a variety of organisms points to convergent evolution on the mitochondria associated with a physiological response to oxygen deprivation or temperature stress, including mechanisms for high-altitude adaptation. Here, we examine whether demography and/or selection explains standing mitogenome nucleotide diversity in high-altitude adapted populations of three Andean waterfowl species: yellow-billed pintail (Anas georgica), speckled teal (Anas flavirostris), and cinnamon teal (Spatula cyanoptera). We compared a total of 60 mitogenomes from each of these three duck species (n = 20 per species) across low and high altitudes and tested whether part(s) or all of the mitogenome exhibited expected signatures of purifying selection within the high-altitude populations of these species. Historical effective population sizes (Ne) were inferred to be similar between high- and low-altitude populations of each species, suggesting that selection rather than genetic drift best explains the reduced genetic variation found in mitochondrial genes of high-altitude populations compared to low-altitude populations of the same species. Specifically, we provide evidence that establishment of these three Andean waterfowl species in the high-altitude environment, coincided at least in part with a persistent pattern of negative purifying selection acting on oxidative phosphorylation (OXPHOS) function of the mitochondria. Our results further reveal that the extent of gene-specific purifying selection has been greatest in the speckled teal, the species with the longest history of high-altitude occupancy.

Human-Induced Range Expansions Result in a Recent Hybrid Zone between Sister Species of Ducks.

Landscapes are consistently under pressure from human-induced ecological change, often resulting in shifting species distributions. For some species, changing the geographical breadth of their niche space results in matching range shifts to regions other than those in which they are formally found. In this study, we employ a population genomics approach to assess potential conservation issues arising from purported range expansions into the south Texas Brush Country of two sister species of ducks: mottled (Anas fulvigula) and Mexican (Anas diazi) ducks. Specifically, despite being non-migratory, both species are increasingly being recorded outside their formal ranges, with the northeastward and westward expansions of Mexican and mottled ducks, respectively, perhaps resulting in secondary contact today. We assessed genetic ancestry using thousands of autosomal loci across the ranges of both species, as well as sampled Mexican- and mottled-like ducks from across overlapping regions of south Texas. First, we confirm that both species are indeed expanding their ranges, with genetically pure Western Gulf Coast mottled ducks confirmed as far west as La Salle county, Texas, while Mexican ducks recorded across Texas counties near the USA-Mexico border. Importantly, the first confirmed Mexican × mottled duck hybrids were found in between these regions, which likely represents a recently established contact zone that is, on average, ~100 km wide. We posit that climate- and land use-associated changes, including coastal habitat degradation coupled with increases in artificial habitats in the interior regions of Texas, are facilitating these range expansions. Consequently, continued monitoring of this recent contact event can serve to understand species' responses in the Anthropocene, but it can also be used to revise operational survey areas for mottled ducks.

Genomic Analysis of Aedes aegypti in the Northern Chihuahuan Desert of Texas and Mexico.

Background: Aedes aegypti, is the primary vector of dengue, Chikungunya, Zika, and yellow fever viruses. Both natural and human-impacted landscapes have selective pressures on Ae. aegypti, resulting in strong genomic structure even within close geographical distances. Materials and Methods: We assess the genetic structure of this medically important mosquito species at the northern leading edge of their distribution in Southwestern USA. Ae. aegypti were collected during 2017 in the urban communities of El Paso and Sparks, Texas (USA) and in the city of Ciudad Juárez, Mexico. Results: Thousands of nuclear loci were sequenced across 260 captured Ae. aegypti. First, we recovered the genetic structure of Ae. aegypti following geography, with all four major collection communities being genetically distinct. Importantly, we found population structure and genetic diversity that suggest rapid expansion through active-short distance dispersals, with Anapra being the likely source for the others. Next, tests of selection recovered eight functional genes across six outliers: calmodulin with olfactory receptor function; the protein superfamily C-type lectin with function in mosquito immune system and development; and TATA box binding protein with function in gene regulation. Conclusion: Despite these populations being documented in the early 2000s, we find that selective pressures on specific genes have already occurred and likely facilitate Ae. aegypti range expansion.

Chromosomal-level reference genome of a wild North American mallard (Anas platyrhynchos).

The mallard (Anas platyrhynchos) is one of the most common, economically, and socially important birds around the world. Mallards were not only an important food source for early humans but eventually becoming intimately linked with people as they were domesticated over the last 2,000 years. To date, mallard genomes are largely reconstructed from samples of domestic or unknown genetic heritage. Here, we report the first high-quality genome assembly and annotation of a genetically vetted wild mallard from North America (NAwild_v1.0). The genome was assembled using a combination of shotgun libraries, proximity ligation Chicago, and Dovetail Hi-C libraries. The final assembly is ∼1.04 Gb in size, with 98.3% of the sequence located in 30 full or nearly full chromosome-level scaffolds, and with a N50/L50 of 79.1 Mb/4 scaffolds. We used a combination of gene prediction and similarity approaches to annotate a total of 23,584 functional genes, of which 19,242 were associated to GO terms. The genome assembly and the set of annotated genes yielded a 95.4% completeness score when compared with the BUSCO aves_odb10 dataset. Next, we aligned 3 previously published mallard genomes to ours, and demonstrate how runs of homozygosity and nucleotide diversity are substantially higher and lower, respectively, to ours and how these artificially changed genomes resulted in profoundly different and biased demographic histories. Our wild mallard assembly not only provides a valuable resource to shed light onto genome evolution, speciation, and other adaptive processes, but also helping with identifying functional genes that have been significantly altered during the domestication process.

Frequency and types of alternative breeding strategies employed by nesting American black ducks in North Carolina.

Although most birds are considered to be at least partially monogamous, molecular evidence continues to uncover that many species can have multiple sexual mates. Many species of Waterfowl (Order Anseriformes) consistently deploy alternative breeding strategies, and although cavity nesting species have been well studied, few attempts to understand rates of alternative breeding strategies exist in the Anatini tribe. Here, we assay mitochondrial DNA and thousands of nuclear markers across 20 broods of American black ducks (Anas rubripes; "black duck") that included 19 females and 172 offspring to study population structure as well as types and rates of secondary breeding strategies in coastal North Carolina. First, we report high levels of relatedness among nesting black ducks and offspring and while 17 (of 19) females were of pure black duck descent, three were found to be black duck x mallard (A. platyrhynchos) hybrids. Next, we evaluated for mismatched mitochondrial DNA and paternity identities across each female's clutch to determine types and frequency of alternative or secondary breeding strategies. Although we report that nest parasitism occurred in two nests, 37% (7 of 19) of the sampled nests were multi-paternal as a result of extra-pair copulation. In addition to being part of a mix of strategies used to increase fecundity by successfully breeding females, we posit nest densities providing easier alternative mate access for males also explains high rates of extra-pair copulation among our sampled black ducks. Ultimately, however, while some proportion of females of many species engage in forms of secondary breeding strategies, we conclude that the decision to do so appears to be seasonally flexible for each individual.

Population genetics and geographic origins of mallards harvested in northwestern Ohio.

The genetic composition of mallards in eastern North America has been changed by release of domestically-raised, game-farm mallards to supplement wild populations for hunting. We sampled 296 hatch-year mallards harvested in northwestern Ohio, October-December 2019. The aim was to determine their genetic ancestry and geographic origin to understand the geographic extent of game-farm mallard introgression into wild populations in more westward regions of North America. We used molecular analysis to detect that 35% of samples were pure wild mallard, 12% were early generation hybrids between wild and game-farm mallards (i.e., F1-F3), and the remaining 53% of samples were assigned as part of a hybrid swarm. Percentage of individuals in our study with some form of hybridization with game-farm mallard (65%) was greater than previously detected farther south in the mid-continent (~4%), but less than the Atlantic coast of North America (~ 92%). Stable isotope analysis using δ2Hf suggested that pure wild mallards originated from areas farther north and west than hybrid mallards. More specifically, 17% of all Ohio samples had δ2Hf consistent with more western origins in the prairies, parkland, or boreal regions of the mid-continent of North America, with 55%, 35%, and 10% of these being genetically wild, hybrid swarm, and F3, respectively. We conclude that continued game-farm introgression into wild mallards is not isolated to the eastern population of mallards in North America, and may be increasing and more widespread than previously detected. Mallards in our study had greater incidence of game-farm hybridization than other locales in the mid-continent but less than eastern North American regions suggesting further need to understand game-farm mallard genetic variation and movement across the continent.