Limited genetic differentiation among breeding, molting, and wintering groups of the threatened Steller’s eider: the role of historic and contemporary factors

Genetic diversity and shallow genetic differentiation of the endangered scaly-sided merganser Mergus squamatus

Examining patterns of genetic diversity are crucial for conservation planning on endangered species, while inferring the underlying process of recent anthropogenic habitat modifications in the context potential long-term demographic changes remains challenging. The globally endangered scaly-sided merganser (SSME), Mergus squamatus, is endemic to a narrow range in Northeast Asia, and its population has recently been contracted into two main breeding areas. Although low genetic diversity has been suggested in the Russian population, the genetic status and demographic history of these individuals have not been fully elucidated. We therefore examined the genetic diversity and structure of the breeding populations of the SSME and investigated the relative importance of historical and recent demographic changes to the present-day pattern of genetic diversity. Using 10 nuclear microsatellite (SSR) markers and mitochondrial DNA (mtDNA) control region sequences, we found limited female-inherited genetic diversity and a high level of nuclear genetic diversity. In addition, analysis of both markers consistently revealed significant but weak divergence between the breeding populations. Inconsistent demographic history parameters calculated from mtDNA and bottleneck analysis results based on SSR suggested a stable historical effective population size. By applying approximate Bayesian computation, it was estimated that populations started to genetically diverge from each other due to recent fragmentation events caused by anthropogenic effects rather than isolation during Last Glacial Maximum (LGM) and post-LGM recolonization. These results suggest that limited historical population size and shallow evolutionary history may be potential factors contributing to the contemporary genetic diversity pattern of breeding SSME populations. Conservation efforts should focus on protecting the current breeding habitats from further destruction, with priority given to both the Russian and Chinese population, as well as restoring the connected suitable breeding grounds.

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.

Use of genetic mark-recapture to estimate breeding site fidelity and philopatry in a threatened sea duck population, Alaska-breeding Steller’s eiders

The Steller’s eider Polysticta stelleri is a sea duck that breeds in Arctic tundra regions of Russia and Alaska (USA). The Alaska-breeding population is listed as ‘threatened’ under the US Endangered Species Act because of a perceived contraction of the breeding range in North America. Understanding demography of the listed population is critical for evaluating measures that can lead to increased abundance and thus, long-term viability. Specifically, estimates of return rates to breeding areas by adult females and natal areas by juvenile females are needed for planning effective recovery actions. We used a suite of polymorphic loci to genotype individuals and generated genetic profiles of nesting females and female offspring from nest materials collected between 1995 and 2016 in a ~170 km2 study area near Utqiagvik, Alaska. We analyzed capture histories of genetically identified individuals to estimate breeding site fidelity, temporary emigration, and philopatry. From a sample of 365 nests, we found that breeding site fidelity of adult females was high (0.91 ± 0.07 SE), and temporary emigration was also high (0.77 ± 0.06) and annually variable (range 0.34-0.97). From egg shell remains of 124 hatched females, we observed 9 recaptures as nesting adults, suggesting that philopatry was also high (range 0.6-1.0). Given the relatively high rates of adult female breeding site fidelity and female philopatry that we estimated, management actions that reduce mortality of adult females and increase annual productivity are likely to help maintain the population of a few hundred breeding Steller’s eiders on the Arctic Coastal Plain of Alaska.

Coast to coast: High genomic connectivity in North American scoters

Dispersal shapes demographic processes and therefore is fundamental to understanding biological, ecological, and evolutionary processes acting within populations. However, assessing population connectivity in scoters (Melanitta sp.) is challenging as these species have large spatial distributions that span remote landscapes, have varying nesting distributions (disjunct vs. continuous), exhibit unknown levels of dispersal, and vary in the timing of the formation of pair bonds (winter vs. fall/spring migration) that may influence the distribution of genetic diversity. Here, we used double-digest restriction-associated DNA sequence (ddRAD) and microsatellite genotype data to assess population structure within the three North American species of scoter (black scoter, M. americana; white-winged scoter, M. deglandi; surf scoter, M. perspicillata), and between their European congeners (common scoter, M. nigra; velvet scoter, M. fusca). We uncovered no or weak genomic structure (ddRAD Φ ST < 0.019; microsatellite F ST < 0.004) within North America but high levels of structure among European congeners (ddRAD Φ ST > 0.155, microsatellite F ST > 0.086). The pattern of limited genomic structure within North America is shared with other sea duck species and is often attributed to male-biased dispersal. Further, migratory tendencies (east vs. west) of female surf and white-winged scoters in central Canada are known to vary across years, providing additional opportunities for intracontinental dispersal and a mechanism for the maintenance of genomic connectivity across North America. In contrast, the black scoter had relatively elevated levels of divergence between Alaska and Atlantic sites and a second genetic cluster found in Alaska at ddRAD loci was concordant with its disjunct breeding distribution suggestive of a dispersal barrier (behavioral or physical). Although scoter populations appear to be connected through a dispersal network, a small percentage (<4%) of ddRAD loci had elevated divergence which may be useful in linking areas (nesting, molting, staging, and wintering) throughout the annual cycle.

Migratory culture, population structure and stock identity in North Pacific beluga whales (Delphinapterus leucas)

The annual return of beluga whales, Delphinapterus leucas, to traditional seasonal locations across the Arctic may involve migratory culture, while the convergence of discrete summering aggregations on common wintering grounds may facilitate outbreeding. Natal philopatry and cultural inheritance, however, has been difficult to assess as earlier studies were of too short a duration, while genetic analyses of breeding patterns, especially across the beluga's Pacific range, have been hampered by inadequate sampling and sparse information on wintering areas. Using a much expanded sample and genetic marker set comprising 1,647 whales, spanning more than two decades and encompassing all major coastal summering aggregations in the Pacific Ocean, we found evolutionary-level divergence among three geographic regions: the Gulf of Alaska, the Bering-Chukchi-Beaufort Seas, and the Sea of Okhotsk (Φst = 0.11-0.32, Rst = 0.09-0.13), and likely demographic independence of (Fst-mtDNA = 0.02-0.66), and in many cases limited gene flow (Fst-nDNA = 0.0-0.02; K = 5-6) among, summering groups within regions. Assignment tests identified few immigrants within summering aggregations, linked migrating groups to specific summering areas, and found that some migratory corridors comprise whales from multiple subpopulations (PBAYES = 0.31:0.69). Further, dispersal is male-biased and substantial numbers of closely related whales congregate together at coastal summering areas. Stable patterns of heterogeneity between areas and consistently high proportions (~20%) of close kin (including parent-offspring) sampled up to 20 years apart within areas (G = 0.2-2.9, p>0.5) is the first direct evidence of natal philopatry to migration destinations in belugas. Using recent satellite telemetry findings on belugas we found that the spatial proximity of winter ranges has a greater influence on the degree of both individual and genetic exchange than summer ranges (rwinter-Fst-mtDNA = 0.9, rsummer-Fst-nDNA = 0.1). These findings indicate widespread natal philopatry to summering aggregation and entire migratory circuits, and provide compelling evidence that migratory culture and kinship helps maintain demographically discrete beluga stocks that can overlap in time and space.

Spatial heterogeneity as a genetic mixing mechanism in highly philopatric colonial seabirds

How genetic diversity is maintained in philopatric colonial systems remains unclear, and understanding the dynamic balance of philopatry and dispersal at all spatial scales is essential to the study of the evolution of coloniality. In the King penguin, Aptenodytes patagonicus, return rates of post-fledging chicks to their natal sub-colony are remarkably high. Empirical studies have shown that adults return year after year to their previous breeding territories within a radius of a few meters. Yet, little reliable data are available on intra- and inter-colonial dispersal in this species. Here, we present the first fine-scale study of the genetic structure in a king penguin colony in the Crozet Archipelago. Samples were collected from individual chicks and analysed at 8 microsatellite loci. Precise geolocation data of hatching sites and selective pressures associated with habitat features were recorded for all sampling locations. We found that despite strong natal and breeding site fidelity, king penguins retain a high degree of panmixia and genetic diversity. Yet, genetic structure appears markedly heterogeneous across the colony, with higher-than-expected inbreeding levels, and local inbreeding and relatedness hotspots that overlap predicted higher-quality nesting locations. This points towards heterogeneous population structure at the sub-colony level, in which fine-scale environmental features drive local philopatric behaviour, while lower-quality patches may act as genetic mixing mechanisms at the colony level. These findings show how a lack of global genetic structuring can emerge from small-scale heterogeneity in ecological parameters, as opposed to the classical model of homogeneous dispersal. Our results also emphasize the importance of sampling design for estimation of population parameters in colonial seabirds, as at high spatial resolution, basic genetic features are shown to be location-dependent. Finally, this study stresses the importance of understanding intra-colonial dispersal and genetic mixing mechanisms in order to better estimate species-wide gene flows and population dynamics.

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.

Multilocus phylogeography (mitochondrial, autosomal and Z-chromosomal loci) and genetic consequence of long-distance male dispersal in Black-throated tits (Aegithalos concinnus)

Multilocus data from the different genomes are essential to understand the phylogeographic history of species, particularly when a species has the male-biased dispersal pattern. Although Black-throated tits (Aegithalos concinnus) are socially monogamous and cooperatively breeding birds, limited observational data suggested that males may have the ability of long-distance dispersal. We have previously detected three highly supported mitochondrial populations within two subspecies of Black-throated tits (A. c. concinnus and A. c. talifuensis). Here, we used several genetic markers with different inheritance patterns to gain insights about their phylogeographic history. Phylogenetic and individual-based Bayesian analysis showed weak geographical structure amongst nuclear sequences (autosomal and Z-chromosomal loci). Coalescent analysis revealed high levels of gene flow among mitochondrial populations, even between allopatric populations. These results strongly suggested that male-biased gene flow was responsible for the discordant cytonuclear phylogeographic patterns. Consistent with expectation on the genetic consequence of long-distance male dispersal, mantel tests revealed a significant pattern of isolation by distance for mitochondrial sequences, but failed to provide a similar pattern for nuclear genes within a continuous population; female Black-throated tits showed a stronger but not significantly different relationship of isolation by distance than males when using mitochondrial DNA alone. We discussed the contribution of male juveniles with delayed dispersal to the non-significantly different IBD patterns between sexes. Our results using multilocus genetic data revealed aspects of the complex evolutionary history of Black-throated tits and the important role of long-distance male dispersal in the population structuring.

Do common eiders nest in kin groups? Microgeographic genetic structure in a philopatric sea duck

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.

Spring migration routes and chronology of surf scoters (Melanitta perspicillata): a synthesis of Pacific coast studies

Understanding interconnectivity among wintering, stopover, and breeding areas of migratory birds is pivotal to discerning how events occurring in each might have a cross-seasonal effect on another. Such information can guide the location and timing of conservation efforts. Thus, we examined spring migration routes, chronology, and stopover use of 85 surf scoters ( Melanitta perspicillata (L., 1758)) marked with satellite transmitters at four Pacific Flyway wintering sites: San Quintin Bay, Baja California; San Francisco Bay, California; Puget Sound, Washington; and Strait of Georgia, British Columbia. Eighty-three percent of marked scoters followed two main routes to the breeding area: a Southern Inland route involving staging in Puget Sound and Strait of Georgia and protracted inland migration, or a Northern Coastal route characterized by short movements along the Pacific coast of British Columbia and southeast Alaska with inland migration initiating from Lynn Canal and surrounding areas. Route choice was related to nesting site latitude in the Canadian Northern Boreal Forest. Data from birds tracked over 2 years indicated strong migration route fidelity, but altered chronology and stopover locations between years. Departure date varied by wintering site, but arrival and apparent settling dates were synchronous, suggesting individuals adjusted migration timing to meet an optimized reproductive schedule.

Population genetic study of mitochondrial DNA in Roseate spoonbill (Aves; Platalea ajaja) breeding colonies from the Pantanal wetlands, Brazil

Five breeding colonies of the Roseate spoonbill (Aves: Platalea ajaja) from two Brazilian wetland areas (Pantanal and Taim marshes) were sampled, and domain I of the mitochondrial DNA control region (483 bp) was sequenced in 50 birds. The average haplotype diversity (h = 0.75, s = 0.071) and average nucleotide diversity (pi = 0.004, s = 0.003) were evaluated, and nonsignificant differences were found among the colonies studied. The lack of differentiation among breeding colonies revealed by AMOVA analysis was explained either as a consequence of high gene flow or recent expansion. The significantly negative results of the neutrality tests (Fu's F ( s ) = -23.271, P < 0.01; Tajima's D = -1.941, P < 0.01) associated with the star shape of the haplotype tree and mismatch distribution data are evidence supporting the idea that these populations underwent a recent demographic expansion in the Pantanal region. The average time since the expansion is estimated to be 25,773 years, and this agrees with a period of increased moisture that occurred during the last glacial period.

Multilocus phylogeography of a holarctic duck: colonization of north america from eurasia by gadwall (Anas strepera)

More than 100 species of birds have Holarctic distributions extending across Eurasia and North America, and many of them likely achieved these distributions by recently colonizing one continent from the other. Mitochondrial DNA (mtDNA) and five nuclear introns were sequenced to test the direction and timing of colonization for a Holarctic duck, the gadwall (Anas strepera). Three lines of evidence suggest gadwalls colonized North America from Eurasia. First, New World (NW) gadwalls had fewer alleles at every locus and 61% of the allelic richness found in Old World (OW) gadwalls. Second, NW gadwalls had lower mtDNA allelic richness than other NW ducks. Third, coalescent analysis suggested that less than 5% of the ancestral population contributed to NW gadwalls at the time of divergence. Gadwalls likely colonized North America during the Late Pleistocene (approximately 81,000 years ago), but the confidence interval on that estimate was large (8500-450,000 years ago). Intercontinental gene flow and selection also likely contributed to genetic diversity in gadwalls. This study illustrates the use of multiple loci and coalescent analyses for critically testing a priori hypotheses regarding dispersal and colonization and provides an independent datapoint supporting an OW to NW bias in the direction of colonization.