DOI: 10.5061/dryad.qv9s4mwsq
title: Historical, abiotic, and biotic drivers influence contemporary lacustrine fish community composition in the glacial Lake Agassiz basin

description: 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.
description: # Data and code from: Historical, abiotic, and biotic drivers influence
 contemporary lacustrine fish community composition in the glacial Lake
 Agassiz basin 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