Fine-scale sampling reveals distinct isolation by distance patterns in chum salmon (Oncorhynchus keta) populations occupying a glacially dynamic environment

Population structure of chum salmon and selection on the markers collected for stock identification

Genetic stock identification (GSI) is a major management tool of Pacific salmon (Oncorhynchus Spp.) that has provided rich genetic baseline data of allozymes, microsatellites, and single-nucleotide polymorphisms (SNPs) across the Pacific Rim. Here, we analyzed published data sets for adult chum salmon (Oncorhynchus keta), namely 10 microsatellites, 53 SNPs, and a mitochondrial DNA locus (mtDNA3, control region, and NADH-3 combined) in samples from 495 locations in the same distribution range (n = 61,813). TreeMix analysis of the microsatellite loci identified the greatest convergence toward Japanese/Korean populations and suggested two admixture events from Japan/Korea to Russia and the Alaskan Peninsula. The SNPs had been purposively collected from rapidly evolving genes to increase the power of GSI. The largest expected heterozygosity was observed in Japanese/Korean populations for microsatellites, whereas it was largest in Western Alaskan populations for SNPs, reflecting the SNP discovery process. A regression of SNP population structures on those of microsatellites indicated the selection of the SNP loci according to deviations from the predicted structures. Specifically, we matched the sampling locations of the SNPs with those of the microsatellites and performed regression analyses of SNP allele frequencies on a 2-dimensional scaling (MDS) of matched locations obtained from microsatellite pairwise F ST values. The MDS first axis indicated a latitudinal cline in American and Russian populations, whereas the second axis showed differentiation of Japanese/Korean populations. The top five outlier SNPs included mtDNA3, U502241 (unknown), GnRH373, ras1362, and TCP178, which were identified by principal component analysis. We summarized the functions of 53 nuclear genes surrounding SNPs and the mtDNA3 locus by referring to a gene database system and propose how they may influence the fitness of chum salmon.

Understanding population structure in an evolutionary context: population-specific FST and pairwise FST

Populations are shaped by their history. It is crucial to interpret population structure in an evolutionary context. Pairwise F_ST measures population structure, whereas population-specific F_ST measures deviation from the ancestral population. To understand the current population structure and a population’s history of range expansion, we propose a representation method that overlays population-specific F_ST estimates on a sampling location map, and on an unrooted neighbor-joining tree and a multi-dimensional scaling plot inferred from a pairwise F_ST distance matrix. We examined the usefulness of our procedure using simulations that mimicked population colonization from an ancestral population and by analyzing published human, Atlantic cod, and wild poplar data. Our results demonstrated that population-specific F_ST values identify the source population and trace the evolutionary history of its derived populations. Conversely, pairwise F_ST values represent the current population structure. By integrating the results of both estimators, we obtained a new picture of the population structure that incorporates evolutionary history. The generalized least squares estimate of genome-wide population-specific F_ST indicated that the wild poplar population expanded its distribution to the north, where daylight hours are long in summer, to coastal areas with abundant rainfall, and to the south where summers are dry. Genomic data highlight the power of the bias-corrected moment estimators of F_ST, whether global, pairwise, or population-specific, that provide unbiased estimates of F_ST. All F_ST moment estimators described in this paper have reasonable processing times and are useful in population genomics studies.

Temporal dynamics of genetic clines of invasive European green crab (Carcinus maenas) in eastern North America

Two genetically distinct lineages of European green crabs (Carcinus maenas) were independently introduced to eastern North America, the first in the early 19th century and the second in the late 20th century. These lineages first came into secondary contact in southeastern Nova Scotia, Canada (NS), where they hybridized, producing latitudinal genetic clines. Previous studies have documented a persistent southward shift in the clines of different marker types, consistent with existing dispersal and recruitment pathways. We evaluated current clinal structure by quantifying the distribution of lineages and fine-scale hybridization patterns across the eastern North American range (25 locations, ~39 to 49°N) using informative single nucleotide polymorphisms (SNPs; n = 96). In addition, temporal changes in the genetic clines were evaluated using mitochondrial DNA and microsatellite loci (n = 9-11) over a 15-year period (2000-2015). Clinal structure was consistent with prior work demonstrating the existence of both northern and southern lineages with a hybrid zone occurring between southern New Brunswick (NB) and southern NS. Extensive later generation hybrids were detected in this region and in southeastern Newfoundland. Temporal genetic analysis confirmed the southward progression of clines over time; however, the rate of this progression was slower than predicted by forecasting models, and current clines for all marker types deviated significantly from these predictions. Our results suggest that neutral and selective processes contribute to cline dynamics, and ultimately, highlight how selection, hybridization, and dispersal can collectively influence invasion success.

Fine-scale temperature-associated genetic structure between inshore and offshore populations of sea scallop (Placopecten magellanicus)

In the northwest Atlantic Ocean, sea scallop (Placopecten magellanicus) has been characterized by a latitudinal genetic cline with a breakpoint between northern and southern genetic clusters occurring at ~45°N along eastern Nova Scotia, Canada. Using 96 diagnostic single-nucleotide polymorphisms (SNPs) capable of discriminating between northern and southern clusters, we examined fine-scale genetic structure of scallops among 27 sample locations, spanning the largest geographic range evaluated in this species to date (~37-51°N). Here, we confirmed previous observations of northern and southern groups, but we show that the boundary between northern and southern clusters is not a discrete latitudinal break. Instead, at latitudes near the previously described boundary, we found unexpected patterns of fine-scale genetic structure occurring between inshore and offshore sites. Scallops from offshore sites, including St. Pierre Bank and the eastern Scotian Shelf, clustered with southern stocks, whereas inshore sites at similar latitudes clustered with northern stocks. Our analyses revealed significant genetic divergence across small spatial scales (i.e., 129-221 km distances), and that spatial structure over large and fine scales was strongly associated with temperature during seasonal periods of thermal minima. Clear temperature differences between inshore and offshore locations may explain the fine-scale structuring observed, such as why southern lineages of scallop occur at higher latitudes in deeper, warmer offshore waters. Our study supports growing evidence that fine-scale population structure in marine species is common, often environmentally associated, and that consideration of environmental and genomic data can significantly enhance the identification of marine diversity and management units.

Congruent population structure across paralogous and nonparalogous loci in Salish Sea chum salmon (Oncorhynchus keta)

Whole-genome duplications are major evolutionary events with a lasting impact on genome structure. Duplication events complicate genetic analyses as paralogous sequences are difficult to distinguish; consequently, paralogs are often excluded from studies. The effects of an ancient whole-genome duplication (approximately 88 MYA) are still evident in salmonids through the persistence of numerous paralogous gene sequences and partial tetrasomic inheritance. We use restriction site-associated DNA sequencing on 10 collections of chum salmon from the Salish Sea in the USA and Canada to investigate genetic diversity and population structure in both tetrasomic and rediploidized regions of the genome. We use a pedigree and high-density linkage map to identify paralogous loci and to investigate genetic variation across the genome. By applying multivariate statistical methods, we show that it is possible to characterize paralogous loci and that they display similar patterns of population structure as the diploidized portion of the genome. We find genetic associations with the adaptively important trait of run-timing in both sets of loci. By including paralogous loci in genome scans, we can observe evolutionary signals in genomic regions that have routinely been excluded from population genetic studies in other polyploid-derived species.

Echoes of a distant time: effects of historical processes on contemporary genetic patterns in Galaxias platei in Patagonia

Interpreting the genetic structure of a metapopulation as the outcome of gene flow over a variety of timescales is essential for the proper understanding of how changes in landscape affect biological connectivity. Here we contrast historical and contemporary connectivity in two metapopulations of the freshwater fish Galaxias platei in northern and southernmost Patagonia where paleolakes existed during the Holocene and Pleistocene, respectively. Contemporary gene flow was mostly high and asymmetrical in the northern system while extremely reduced in the southernmost system. Historical migration patterns were high and symmetric in the northern system and high and largely asymmetric in the southern system. Both systems showed a moderate structure with a clear pattern of isolation by distance (IBD). Effective population sizes were smaller in populations with low contemporary gene flow. An approximate Bayesian computation (ABC) approach suggests a late Holocene colonization of the lakes in the northern system and recent divergence of the populations from refugial populations from east and west of the Andes. For the southern system, the ABC approach reveals that some of the extant G. platei populations most likely derive from an ancestral population inhabiting a large Pleistocene paleolake while the rest derive from a higher-altitude lake. Our results suggest that neither historical nor contemporary processes individually fully explain the observed structure and geneflow patterns and both are necessary for a proper understanding of the factors that affect diversity and its distribution. Our study highlights the importance of a temporal perspective on connectivity to analyse the diversity of spatially complex metapopulations.

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.