1
|
Kitada S, Kishino H. Population structure of chum salmon and selection on the markers collected for stock identification. Ecol Evol 2021; 11:13972-13985. [PMID: 34707832 PMCID: PMC8525185 DOI: 10.1002/ece3.8102] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/10/2021] [Revised: 08/25/2021] [Accepted: 08/27/2021] [Indexed: 11/25/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Shuichi Kitada
- Tokyo University of Marine Science and TechnologyTokyoJapan
| | - Hirohisa Kishino
- Graduate School of Agriculture and Life SciencesThe University of TokyoTokyoJapan
- Present address:
The Research Institute of Evolutionary BiologyTokyoJapan
| |
Collapse
|
2
|
Kitada S, Nakamichi R, Kishino H. Understanding population structure in an evolutionary context: population-specific FST and pairwise FST. G3-GENES GENOMES GENETICS 2021; 11:6364900. [PMID: 34549777 PMCID: PMC8527463 DOI: 10.1093/g3journal/jkab316] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 01/20/2021] [Accepted: 08/27/2021] [Indexed: 01/04/2023]
Abstract
Populations are shaped by their history. It is crucial to interpret population structure in an evolutionary context. Pairwise FST measures population structure, whereas population-specific FST 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 FST estimates on a sampling location map, and on an unrooted neighbor-joining tree and a multi-dimensional scaling plot inferred from a pairwise FST 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 FST values identify the source population and trace the evolutionary history of its derived populations. Conversely, pairwise FST 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 FST 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 FST, whether global, pairwise, or population-specific, that provide unbiased estimates of FST. All FST moment estimators described in this paper have reasonable processing times and are useful in population genomics studies.
Collapse
Affiliation(s)
- Shuichi Kitada
- Tokyo University of Marine Science and Technology, Tokyo 108-8477, Japan
| | | | - Hirohisa Kishino
- Graduate School of Agriculture and Life Sciences, The University of Tokyo, Tokyo 113-8657, Japan.,The Research Institute of Evolutionary Biology, Tokyo 138-0098, Japan
| |
Collapse
|
3
|
Lehnert SJ, DiBacco C, Jeffery NW, Blakeslee AMH, Isaksson J, Roman J, Wringe BF, Stanley RRE, Matheson K, McKenzie CH, Hamilton LC, Bradbury IR. Temporal dynamics of genetic clines of invasive European green crab ( Carcinus maenas) in eastern North America. Evol Appl 2018; 11:1656-1670. [PMID: 30344634 PMCID: PMC6183463 DOI: 10.1111/eva.12657] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/28/2017] [Revised: 05/28/2018] [Accepted: 06/03/2018] [Indexed: 12/12/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Sarah J. Lehnert
- Northwest Atlantic Fisheries CentreFisheries and Oceans CanadaSt. John'sNewfoundlandCanada
| | - Claudio DiBacco
- Bedford Institute of OceanographyFisheries and Oceans CanadaDartmouthNova ScotiaCanada
| | - Nicholas W. Jeffery
- Bedford Institute of OceanographyFisheries and Oceans CanadaDartmouthNova ScotiaCanada
| | | | - Jonatan Isaksson
- Gund Institute for EnvironmentUniversity of VermontBurlingtonVermont
| | - Joe Roman
- Gund Institute for EnvironmentUniversity of VermontBurlingtonVermont
| | - Brendan F. Wringe
- Bedford Institute of OceanographyFisheries and Oceans CanadaDartmouthNova ScotiaCanada
| | - Ryan R. E. Stanley
- Bedford Institute of OceanographyFisheries and Oceans CanadaDartmouthNova ScotiaCanada
| | - Kyle Matheson
- Northwest Atlantic Fisheries CentreFisheries and Oceans CanadaSt. John'sNewfoundlandCanada
| | - Cynthia H. McKenzie
- Northwest Atlantic Fisheries CentreFisheries and Oceans CanadaSt. John'sNewfoundlandCanada
| | - Lorraine C. Hamilton
- Aquatic Biotechnology LaboratoryBedford Institute of OceanographyDartmouthNova ScotiaCanada
| | - Ian R. Bradbury
- Northwest Atlantic Fisheries CentreFisheries and Oceans CanadaSt. John'sNewfoundlandCanada
| |
Collapse
|
4
|
Fine-scale temperature-associated genetic structure between inshore and offshore populations of sea scallop (Placopecten magellanicus). Heredity (Edinb) 2018; 122:69-80. [PMID: 29773897 DOI: 10.1038/s41437-018-0087-9] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2017] [Revised: 04/16/2018] [Accepted: 04/16/2018] [Indexed: 11/08/2022] Open
Abstract
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.
Collapse
|
5
|
Waples RK, Seeb JE, Seeb LW. Congruent population structure across paralogous and nonparalogous loci in Salish Sea chum salmon (Oncorhynchus keta). Mol Ecol 2017; 26:4131-4144. [PMID: 28452089 DOI: 10.1111/mec.14163] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/14/2016] [Revised: 04/01/2017] [Accepted: 04/10/2017] [Indexed: 12/19/2022]
Abstract
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.
Collapse
Affiliation(s)
- R K Waples
- School of Aquatic and Fishery Sciences, University of Washington, Seattle, WA, USA
| | - J E Seeb
- School of Aquatic and Fishery Sciences, University of Washington, Seattle, WA, USA
| | - L W Seeb
- School of Aquatic and Fishery Sciences, University of Washington, Seattle, WA, USA
| |
Collapse
|
7
|
Vera-Escalona I, Habit E, Ruzzante DE. Echoes of a distant time: effects of historical processes on contemporary genetic patterns in Galaxias platei in Patagonia. Mol Ecol 2015; 24:4112-28. [PMID: 26147523 DOI: 10.1111/mec.13303] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/20/2015] [Revised: 06/21/2015] [Accepted: 06/24/2015] [Indexed: 12/11/2022]
Abstract
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.
Collapse
Affiliation(s)
- Iván Vera-Escalona
- Department of Biology, Dalhousie University, 1355 Oxford St., Halifax, NS, Canada
| | - Evelyn Habit
- Departamento de Sistemas Acuáticos, Facultad de Ciencias Ambientales y Centro EULA-Chile, Universidad de Concepción, Barrio Universitario s/n, Concepción, Chile
| | - Daniel E Ruzzante
- Department of Biology, Dalhousie University, 1355 Oxford St., Halifax, NS, Canada
| |
Collapse
|
8
|
Burridge CP, Peucker AJ, Valautham SK, Styan CA, Dann P. Nonequilibrium Conditions Explain Spatial Variability in Genetic Structuring of Little Penguin (Eudyptula minor). J Hered 2015; 106:228-37. [PMID: 25833231 PMCID: PMC4406270 DOI: 10.1093/jhered/esv009] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/24/2014] [Accepted: 02/09/2015] [Indexed: 11/28/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Christopher P Burridge
- From the School of Biological Sciences, University of Tasmania, Private Bag 55, Hobart, Tasmania 7001, Australia (Burridge and Valautham); the School of Life and Environmental Sciences, Deakin University, Warrnambool, Victoria 3280, Australia (Peucker and Styan); the School of Energy and Resources, UCL Australia, Adelaide, South Australia 5000, Australia (Styan); and the Research Department, Phillip Island Nature Parks, Cowes, Victoria 3922, Australia (Dann).
| | - Amanda J Peucker
- From the School of Biological Sciences, University of Tasmania, Private Bag 55, Hobart, Tasmania 7001, Australia (Burridge and Valautham); the School of Life and Environmental Sciences, Deakin University, Warrnambool, Victoria 3280, Australia (Peucker and Styan); the School of Energy and Resources, UCL Australia, Adelaide, South Australia 5000, Australia (Styan); and the Research Department, Phillip Island Nature Parks, Cowes, Victoria 3922, Australia (Dann)
| | - Sureen K Valautham
- From the School of Biological Sciences, University of Tasmania, Private Bag 55, Hobart, Tasmania 7001, Australia (Burridge and Valautham); the School of Life and Environmental Sciences, Deakin University, Warrnambool, Victoria 3280, Australia (Peucker and Styan); the School of Energy and Resources, UCL Australia, Adelaide, South Australia 5000, Australia (Styan); and the Research Department, Phillip Island Nature Parks, Cowes, Victoria 3922, Australia (Dann)
| | - Craig A Styan
- From the School of Biological Sciences, University of Tasmania, Private Bag 55, Hobart, Tasmania 7001, Australia (Burridge and Valautham); the School of Life and Environmental Sciences, Deakin University, Warrnambool, Victoria 3280, Australia (Peucker and Styan); the School of Energy and Resources, UCL Australia, Adelaide, South Australia 5000, Australia (Styan); and the Research Department, Phillip Island Nature Parks, Cowes, Victoria 3922, Australia (Dann)
| | - Peter Dann
- From the School of Biological Sciences, University of Tasmania, Private Bag 55, Hobart, Tasmania 7001, Australia (Burridge and Valautham); the School of Life and Environmental Sciences, Deakin University, Warrnambool, Victoria 3280, Australia (Peucker and Styan); the School of Energy and Resources, UCL Australia, Adelaide, South Australia 5000, Australia (Styan); and the Research Department, Phillip Island Nature Parks, Cowes, Victoria 3922, Australia (Dann)
| |
Collapse
|