Is Niagara Falls a barrier to gene flow in riverine fishes? A test using genome-wide SNP data from seven native species
Data files
Jun 14, 2021 version files 8.23 MB
-
Ambloplites.genepop.gen
217.25 KB
-
Ambloplites.usnps.geno
30.03 KB
-
Ambloplites.usnps.geno.xls
54.27 KB
-
Ambloplites.vcf
222.18 KB
-
Ameiurus.genepop.gen
384.69 KB
-
Ameiurus.usnps.geno
53.47 KB
-
Ameiurus.usnps.geno.xls
61.95 KB
-
Ameiurus.vcf
379.83 KB
-
Catostomus.genepop.gen
2.08 MB
-
Catostomus.usnps.geno
289.44 KB
-
Catostomus.usnps.geno.xls
74.24 KB
-
Catostomus.vcf
2.13 MB
-
Micropterus.genepop.gen
229.44 KB
-
Micropterus.usnps.geno
31.70 KB
-
Micropterus.usnps.geno.xls
60.93 KB
-
Micropterus.vcf
230.33 KB
-
MODEL_MSFS-2.obs
337 B
-
MODEL_MSFS-3.obs
399 B
-
MODEL_MSFS-4.obs
268 B
-
MODEL_MSFS-5.obs
362 B
-
MODEL_MSFS-6.obs
348 B
-
MODEL_MSFS.obs
366 B
-
MODEL-2.est
501 B
-
MODEL-2.tpl
808 B
-
MODEL-3.est
620 B
-
MODEL-3.tpl
805 B
-
MODEL-4.est
501 B
-
MODEL-4.tpl
790 B
-
MODEL-5.est
501 B
-
MODEL-5.tpl
793 B
-
MODEL-6.est
620 B
-
MODEL-6.tpl
790 B
-
MODEL-7.est
501 B
-
MODEL.est
382 B
-
MODEL.tpl
805 B
-
Mox_macro.genepop.gen
303.55 KB
-
Mox_macro.usnps.geno
40.90 KB
-
Mox_macro.usnps.geno.xls
46.59 KB
-
Mox_macro.vcf
529.40 KB
-
Mox_valen.genepop.gen
30.87 KB
-
Mox_valenc.usnps.geno
4.06 KB
-
Mox_valenc.usnps.geno.xls
33.28 KB
-
Mox_valenc.vcf
69.60 KB
-
Perca.genepop.gen
259.47 KB
-
Perca.usnps.geno
36.10 KB
-
Perca.usnps.geno.xls
73.22 KB
-
Perca.vcf
246.73 KB
-
Values_used_in_fsc2_analysis.xlsx
9.26 KB
Abstract
Since the early Holocene, fish population genetics in the Laurentian Great Lakes have been shaped by the dual influences of habitat structure and post-glacial dispersal. Riverscape genetics theory predicts that longitudinal habitat corridors and unidirectional downstream water-flow drive the downstream accumulation of genetic diversity, whereas post-glacial dispersal theory predicts that fish genetic diversity should decrease with increasing distance from glacial refugia. This study examines populations of seven native fish species codistributed above and below the 58 m-high Niagara Falls – a hypothesized barrier to gene flow in aquatic species. A better understanding of Niagara Falls’ role as a barrier to gene flow and dispersal is needed to identify drivers of Great Lakes genetic diversity and guide strategies to limit exotic species invasions. We used genome-wide SNPs and coalescent models to test whether populations are: (1) genetically distinct, consistent with the Niagara Falls barrier hypothesis; (2) more genetically diverse upstream, consistent with post-glacial expansion theory, or downstream, consistent with the riverscape habitat theory; and, (3) are migrating either upstream or downstream past Niagara Falls. We found that genetic diversity is consistently greater below Niagara Falls and the falls are an effective barrier to downstream migration, but at least some species have likely dispersed upstream past the falls since the time of their formation yet before construction of the Welland Canal. Models restricting migration to after opening of the Welland Canal were generally rejected. These results help explain how river habitat features affect aquatic species’ genetic diversity and highlight the need to better understand post-glacial dispersal pathways.
Methods
These are both original data and permutations of data generated via a three-enzyme restriction site associated DNA (3RAD) library preparationa and sequencing pipeline applied to seven native fish species that are codistributed above and below Niagara Falls, in the Niagara River between lakes Erie and Ontario at the Canada – United States border. The seven species examined are: Ambloplites rupestris, Ameiurus nebulosus, Catostomus commersoni, Micropterus salmoides, Moxostoma macrolepidotum, Moxostoma valenciennesi, and Perca flavescens. Raw 3RAD sequence reads were clustered and SNPs were called using the program Pyrad. Pyrad outputs provided here are the unlinked SNP matrices and the VCF files generated for each species. Also provided are the input files for our principle coordinate analyses in .xls format, input files for our diveRsity analyses in Genepop format, and input files for our fastsimcoal2 demographic model testing analyses.
Usage notes
Each data set contains data for each of the seven species *except* the fastsimcoal2 datasets, which exclude Moxostoma valenciennesi, for which there were two few individuals to run demographic analyses.