Historical, abiotic, and biotic drivers influence contemporary lacustrine fish community composition in the glacial Lake Agassiz basin
Data files
Nov 07, 2025 version files 400.67 KB
-
Illyes_and_Mandrak_CJFAS_Lake_Agassiz_code.Rmd
87.70 KB
-
Illyes_and_Mandrak_CJFAS_Lake_Agassiz_data.xlsx
310.10 KB
-
README.md
2.87 KB
Abstract
Historical determinants and local ecological filters tend to drive community assembly at large and smaller spatiotemporal scales, respectively. We examined the effects of differential contact to the Wisconsinan deglaciation dispersal corridor, Lake Agassiz, and habitat heterogeneity and species interactions on contemporary lacustrine fish community composition in 261 northwestern Ontario lakes. Species composition varied significantly across lakes grouped by their most recent glacial lake coverage, driven by turnover, and the number of connections (not-covered, single, and secondary contact) to the dispersal corridor over time, explained by nestedness. Maximum lake depth, elevation, and most recent stage of glacial lake coverage were dominant drivers of community composition patterns; and most of the significant species co-occurrence relationships were positive. While the dominant environmental variables reflect the effects of habitat heterogeneity and isolation from the dispersal corridor, the species co-occurrence analysis results reflect the influence of environmental niche filtering and biotic interactions in shaping lacustrine fish communities. The spatiotemporal evolution of Lake Agassiz played a role as a historical determinant, and habitat heterogeneity and local species interactions were ecological determinants.
Notes for accompanying Illyes and Mandrak_CJFAS_Lake Agassiz_data.xlsx file
Raw fish data with species at risk removed for their protection. The full set of fish data (including species at risk) are available for non-commercial purposes from the Ministry of Natural Resources in Ontario (https://geohub.lio.gov.on.ca) through their standard data-sharing agreements.
Export sheets from Illyes_and_Mandrak_CJFAS_Lake_Agassiz_data.xlsx as .csv to load them in the R script
Sheet 1:
- Lake_Name = name of lake
- LID = unique numerical identifier of lake
- Surf_Area = lake surface area in [ha]
- Max_Depth in [m]
- Z = lake elevation in [m.a.s.l.]
- Species_ri = species richness
- First_stg = first stage of Lake Agassiz coverage
- Last_stg_merge = last stage of Lake Agassiz coverage
- Sec_con_merge = contact type
- Time_fr = time start of Lake Agassiz coverage
- Time_to = time end of Lake Agassiz coverage
- Cover_leng = length of Lake Agassiz coverage in [14C years]
Sheet 2:
- species_code = unique numerical identifier of species
- species_english_name = species common name
- species_latin_name = species scientific name
- species_scientific_name_short = species scientific name abbreviation
- last_stage_merge = last stage of Lake Agassiz coverage
- length_cover = length of Lake Agassiz coverage in [14C years]
- contact_type = contact type
- elevation_masl = lake elevation in [m.a.s.l.]
- area_ha = lake surface area in [ha]
- lake_id = unique numerical identifier of lake
- lake_name = name of lake
Sheet 3: lake group pairwise Jaccard dissimilarity
- Last stagemerge pair 1 = lake group one
- Last stagemerge pair 2 = lake group two
- Last stage merge pair 1 time to = end of Lake Agassiz coverage of lake group one in [14C years]
- Last stage merge pair 2 time to = end of Lake Agassiz coverage of lake group two in [14C years]
- Time difference C14 y = time difference in end of Lake Agassiz coverage between lake group one and two in [14C years]
- Jaccard dissimilarity = dissimilarity between group one and two
- Jaccard dissimilarity nestedness = nestedness component of dissimilarity between group one and two
- Jaccard dissimilarity turnover = turnover component of dissimilarity between group one and two
Sheet 4:
- rows = species by their scientific name abbreviation
- columns = lakes by their LID numbers
See more details on the data in the methods of the manuscript
See Illyes and Illyes_and_Mandrak_CJFAS_Lake_Agassiz_code.Rmd file for R code
Code was created using R 4.2.2 in RStudio V. 2024.12.1+563
Packages used:
betapart
bipartite
cooccur
dplyr
forcats
ggforce
ggplot2
ggpubr
ggrepel
polynom
RColorBrewer
rdacca.hp
tidyverse
vegan
