Secondary contact and changes in coastal habitat availability influence the nonequilibrium population structure of a salmonid (Oncorhynchus keta)

Historical and contemporary processes driving the origin and structure of an admixed population within a contact zone between subspecies of a north temperate diadromous fish

Hybridization between divergent lineages can result in losses of distinct evolutionary taxa. Alternatively, hybridization can lead to increased genetic variability that may fuel local adaptation and the generation of novel traits and/or taxa. Here, we examined single-nucleotide polymorphisms generated using genotyping-by-sequencing in a population of Dolly Varden char (Pisces: Salmonidae) that is highly admixed within a contact zone between two subspecies (Salvelinus malma malma, Northern Dolly Varden [NDV] and S. m. lordi, Southern Dolly Varden [SDV]) in southwestern Alaska to assess the spatial distribution of hybrids and to test hypotheses on the origin of the admixed population. Ancestry analysis revealed that this admixed population is composed of advanced generation hybrids between NDV and SDV or advanced backcrosses to SDV; no F1 hybrids were detected. Coalescent-based demographic modelling supported the origin of this population about 55,000 years ago by secondary contact between NDV and SDV with low levels of contemporary gene flow. Ancestry in NDV and SDV varies within the watershed and ancestry in NDV was positively associated with distance upstream from the sea, contingent on habitat-type sampled, and negatively associated with the number of migrations that individual fish made to the sea. Our results suggest that divergence between subspecies over hundreds of thousands of years may not be associated with significant reproductive isolation, but that elevated diversity owing to hybridization may have contributed to adaptive divergence in habitat use and life history.

Insights from a chum salmon (Oncorhynchus keta) genome assembly regarding whole-genome duplication and nucleotide variation influencing gene function

Chum salmon are ecologically important to Pacific Ocean ecosystems and commercially important to fisheries. To improve the genetic resources available for this species, we sequenced and assembled the genome of a male chum salmon using Oxford Nanopore read technology and the Flye genome assembly software (contig N50: ∼2 Mbp, complete BUSCOs: ∼98.1%). We also resequenced the genomes of 59 chum salmon from hatchery sources to better characterize the genome assembly and the diversity of nucleotide variants impacting phenotype variation. With genomic sequences from a doubled haploid individual, we were able to identify regions of the genome assembly that have been collapsed due to high sequence similarity between homeologous (duplicated) chromosomes. The homeologous chromosomes are relics of an ancient salmonid-specific genome duplication. These regions were enriched with genes whose functions are related to the immune system and responses to toxins. From analyzing nucleotide variant annotations of the resequenced genomes, we were also able to identify genes that have increased levels of variants thought to moderately impact gene function. Genes related to the immune system and the detection of chemical stimuli (olfaction) had increased levels of these variants based on a gene ontology enrichment analysis. The tandem organization of many of the enriched genes raises the question of why they have this organization.

Genetic assessment of the effects of streamscape succession on coho salmon Oncorhynchus kisutch colonization in recently deglaciated streams

Measures of genetic diversity within and among populations and historical geomorphological data on stream landscapes were used in model simulations based on approximate Bayesian computation (ABC) to examine hypotheses of the relative importance of stream features (geomorphology and age) associated with colonization events and gene flow for coho salmon Oncorhynchus kisutch breeding in recently deglaciated streams (50-240 years b.p.) in Glacier Bay National Park (GBNP), Alaska. Population estimates of genetic diversity including heterozygosity and allelic richness declined significantly and monotonically from the oldest and largest to youngest and smallest GBNP streams. Interpopulation variance in allele frequency increased with increasing distance between streams (r = 0·435, P < 0·01) and was inversely related to stream age (r = -0·281, P < 0·01). The most supported model of colonization involved ongoing or recent (<10 generations before sampling) colonization originating from large populations outside Glacier Bay proper into all other GBNP streams sampled. Results here show that sustained gene flow from large source populations is important to recently established O. kisutch metapopulations. Studies that document how genetic and demographic characteristics of newly founded populations vary associated with successional changes in stream habitat are of particular importance to and have significant implications for, restoration of declining or repatriation of extirpated populations in other regions of the species' native range.

Temporal sampling helps unravel the genetic structure of naturally occurring populations of a phytoparasitic nematode. 2. Separating the relative effects of gene flow and genetic drift

Studying wild pathogen populations in natural ecosystems offers the opportunity to better understand the evolutionary dynamics of biotic diseases in crops and to enhance pest control strategies. We used simulations and genetic markers to investigate the spatial and temporal population genetic structure of wild populations of the beet cyst nematode Heterodera schachtii on a wild host plant species, the sea beet (Beta vulgaris spp. maritima), the wild ancestor of cultivated beets. Our analysis of the variation of eight microsatellite loci across four study sites showed that (i) wild H. schachtii populations displayed fine-scaled genetic structure with no evidence of substantial levels of gene flow beyond the scale of the host plant, and comparisons with simulations indicated that (ii) genetic drift substantially affected the residual signals of isolation-by-distance processes, leading to departures from migration-drift equilibrium. In contrast to what can be suspected for (crop) field populations, this showed that wild cyst nematodes have very low dispersal capabilities and are strongly disconnected from each other. Our results provide some key elements for designing pest control strategies, such as decreasing passive dispersal events to limit the spread of virulence among field nematode populations.

Nonequilibrium Conditions Explain Spatial Variability in Genetic Structuring of Little Penguin (Eudyptula minor)

Factors responsible for spatial structuring of population genetic variation are varied, and in many instances there may be no obvious explanations for genetic structuring observed, or those invoked may reflect spurious correlations. A study of little penguins (Eudyptula minor) in southeast Australia documented low spatial structuring of genetic variation with the exception of colonies at the western limit of sampling, and this distinction was attributed to an intervening oceanographic feature (Bonney Upwelling), differences in breeding phenology, or sea level change. Here, we conducted sampling across the entire Australian range, employing additional markers (12 microsatellites and mitochondrial DNA, 697 individuals, 17 colonies). The zone of elevated genetic structuring previously observed actually represents the eastern half of a genetic cline, within which structuring exists over much shorter spatial scales than elsewhere. Colonies separated by as little as 27 km in the zone are genetically distinguishable, while outside the zone, homogeneity cannot be rejected at scales of up to 1400 km. Given a lack of additional physical or environmental barriers to gene flow, the zone of elevated genetic structuring may reflect secondary contact of lineages (with or without selection against interbreeding), or recent colonization and expansion from this region. This study highlights the importance of sampling scale to reveal the cause of genetic structuring.

Is isolation by adaptation driving genetic divergence among proximate Dolly Varden char populations?

Numerous studies of population genetics in salmonids and other anadromous fishes have revealed that population structure is generally organized into geographic hierarchies (isolation by distance), but significant structure can exist in proximate populations due to varying selective pressures (isolation by adaptation). In Chignik Lakes, Alaska, anadromous Dolly Varden char (Salvelinus malma) spawn in nearly all accessible streams throughout the watershed, including those draining directly to an estuary, Chignik Lagoon, into larger rivers, and into lakes. Collections of Dolly Varden fry from 13 streams throughout the system revealed low levels of population structure among streams emptying into freshwater. However, much stronger genetic differentiation was detected between streams emptying into freshwater and streams flowing directly into estuarine environments. This fine-scale reproductive isolation without any physical barriers to migration is likely driven by differences in selection pressures across freshwater and estuarine environments. Estuary tributaries had fewer larger, older juveniles, suggesting an alternative life history of smolting and migration to the marine environment at a much smaller size than occurs in the other populations. Therefore, genetic data were consistent with a scenario where isolation by adaptation occurs between populations of Dolly Varden in the study system, and ecological data suggest that this isolation may partially be a result of a novel Dolly Varden life history of seawater tolerance at a smaller size than previously recognized.

Source-sink estimates of genetic introgression show influence of hatchery strays on wild chum salmon populations in Prince William Sound, Alaska

The extent to which stray, hatchery-reared salmon affect wild populations is much debated. Although experiments show that artificial breeding and culture influence the genetics of hatchery salmon, little is known about the interaction between hatchery and wild salmon in a natural setting. Here, we estimated historical and contemporary genetic population structures of chum salmon (Oncorhynchus keta) in Prince William Sound (PWS), Alaska, with 135 single nucleotide polymorphism (SNP) markers. Historical population structure was inferred from the analysis of DNA from fish scales, which had been archived since the late 1960’s for several populations in PWS. Parallel analyses with microsatellites and a test based on Hardy-Weinberg proportions showed that about 50% of the fish-scale DNA was cross-contaminated with DNA from other fish. These samples were removed from the analysis. We used a novel application of the classical source-sink model to compare SNP allele frequencies in these archived fish-scales (1964–1982) with frequencies in contemporary samples (2008–2010) and found a temporal shift toward hatchery allele frequencies in some wild populations. Other populations showed markedly less introgression, despite moderate amounts of hatchery straying. The extent of introgression may reflect similarities in spawning time and life-history traits between hatchery and wild fish, or the degree that hybrids return to a natal spawning area. The source-sink model is a powerful means of detecting low levels of introgression over several generations.