Using genetics and otolith microchemistry to investigate larval summer flounder connectivity throughout the species range from 1989-2012: Using multiple natural tags provides evidence for extensive larval dispersal across space and through time in summer flounder by Hoey et al. (2020).

These data, methods and scripts are provided in the interests of open science. If you have questions or find errors, please let me know.

Citable as

Contact: Jennifer Hoey ([email protected])

• Larvae Sampling Database.csv: csv of metadata for all sampled larvae, whether sequenced or not
• masterPADElarvae.txt: Data of 523 larval summer flounder, including date, location, 10 outlier loci, SL and developmental stage
• masterPADElarvae_metadata.txt: Metadata pertaining to masterPADElarve.txt
• otolith_data.txt: raw larval otolith microchemistry data
• otolith_data_nozeros.txt: larval otolith microchemistry data where zeros have been randomly replaced by a number between zero and the smallest non-zero number for that element
• pop.allele.counts10.txt: Adult summer flounder allele counts at outlier loci from Hoey & Pinsky (2018)
• SNP.DP3g95maf05lm75.FIL.recode.vcf: all SNPs output following filtering described in Bioinformatics.txt
• SNP.DP3g95maf05lm75.FIL.recode.firstsnp.vcf: same as SNP.DP3g95maf05lm75.FIL.recode.vcf but only first SNP per contig containing 1904 loci across 293 larvae and 235 adult summer flounder
• structure_input_Mar14_2017_528_1904.str: structure-formatted file of 1904 loci across 293 larvae and 235 adult summer flounder used for many of the genetic analyses

• Bioinformatics.txt: Bioinformatic steps taken to go from sequencing reads to SNP genotypes, as well as steps taken to produce filtered SNP genotypes SNP.DP3g95maf05lm75.FIL.recode.vcf

• assignment_confusion_matrix.R: Script to sample n individuals from each of 10 spawning regions, calculate genotype likelihoods of multilocus genotypes, determine most likely population and calculate % assignment and plot confusion matrices
• conceptual_dispersal_diagram.R: Script to plot conceptual_dispersal_diagram.png
• GAMS_to_generate_likelihoods.R: fits GAMS between allele frequencies and distance along coast, calculates observed individual genotype likelihoods across 10 populations, plots indiv_assignments_10pops.png and calculates observed cluster genotype likelihoods from individual likelihoods
• group assignment GAMS.R: Given larval clustering data and GAM-determined adult allele frequencies calculated from GAMS_to_generate_likelihoods.R, simulates larval genotypes for clusters of size N, then calculates the distribution of genotype likelihoods for each of ten adult populations. Plots observed cluster genotype likelihood on top of emperical distribution and calculates a p-value
• otoliths.R: Using larval otolith microchemistry data for NMDS, ANOVAs, LDAs, clustering, clustering pie charts, etc.
• plot allele freqs.R: Plots allele frequencies of 10 GAM-determined locations to produce outlierloci_lineplot_GAM_10pops.png
• recentIBD.R: Fst calculations for larvae captured 2008-2012
• structurethrutime.R: bulk of larval genomic analyses, including Fst, AMOVA
• larvalstructure.R: spatial and temporal PCAs, generating larval_pca_regionalpops_1904loci.png and larval_pca_temppops_1904loci.png

• cluster197_piecharts_bytimeperiod.png: pie charts of larval clusters with otolith-only data for each time period (Figure S8)
• nmds_byperiod_10elements.png: nmds plots for each time period using all 10 elements (Figure 3)
• lda_barplots.png: LDA barplots for training (data within 68% confidence ellipses), test & combined datasets (Figure S9)
• lda_barplots_3timeperiods.png: LDA barplots for each time period separately (Figure S11)
• lda_barplots_late.png: LDA barplots for training (data within 68% confidence ellipses), test & combined datasets for 2008-2012 larvae only
• otolith_overtimelogged_boxplots.png: box plots of otolith microchemistry by time period (Figure S7)
• otolith_sitelogged_boxplots.png: box plots of otolith microchemistry by ingress site (Figure S6)
• post_group_membership_prob.png: frequency distributions of predicted group membership probabilities from LDA classification results (Figure S10)

• larvae293_1904loci.pdf: STRUCTURE plot of 293 larvae using 1904 loci (Figure S3)
• larval_pca_regionalpops_1904loci.png: regional larval PCA (Figure S2A)
• larval_pca_temppops_1904loci.png: temporal larval PCA (Figure S2B)
• indiv_assignments_10pops.png: barplots of individual assignment by ingress site (Figure S5)
• outlierloci_lineplot_GAM_10pops.png: Adult allele frequencies at 10 outlier loci across 10 GAM-determined locations (Figure S1)

• cluster_piecharts_bytimeperiod.png: pie charts of clusters in each time period, colored by ingress site (Figure S12)
• obs_likelihoods_early_6clusters_10pops.txt: Likelihoods of larval clusters from the early time period (1989-1993) across all 10 GAM-determined adult spawning locations, with Pop 10 being northern-most and Pop 1 being southern-most
• obs_likelihoods_middle_2clusters_10pops.txt: Likelihoods of larval clusters from the middle time period (1998-2002) across all 10 GAM-determined adult spawning locations, with Pop 10 being northern-most and Pop 1 being southern-most
• obs_likelihoods_late_3clusters_10pops.txt: Likelihoods of larval clusters from the late time period (2008-2012) across all 10 GAM-determined adult spawning locations, with Pop 10 being northern-most and Pop 1 being southern-most
• conceptual_dispersal_diagram.png: conceptual diagram depicting where larval clusters originated and where they ingressed to (Figure 4)
• indiv_and_cluster_assignments_hists_10GAMpops.png: likelihood differences at the individual versus cluster levels (Figure S4)
• cluster_likelihoods_10GAMSgroups.png: likelihoods of each cluster plotted across 10 GAM-determined locations for each time period

These make up Figure 2A-E

• assignment_confusion_matrix.png: confusion matrix of individuals
• assignment_confusion_matrix_5individuals.png: confusion matrix of groups of 5 individuals
• assignment_confusion_matrix_10individuals.png: confusion matrix of groups of 10 individuals
• assignment_confusion_matrix_20individuals.png: confusion matrix of groups of 20 individuals
• assignment_confusion_matrix_30individuals.png: confusion matrix of groups of 30 individuals