DOI: 10.5061/dryad.q83bk3jg2
title: Population collapse in viviparid gastropods of the Lake Victoria ecoregion started before the Last Glacial Maximum

description: For the purpose of reproducibility, we here provide the datasets and R
 script supporting the analyses of the paper “Population collapse in
 viviparid gastropods of the Lake Victoria ecoregion started before the
 Last Glacial Maximum” by Van Bocxlaer et al. This paper has been accepted
 for publication in Molecular Ecology on 31 July 2020. In this study, we
 examine the population structure of the clade of Bellamya gastropods that
 occupies the Lake Victoria ecoregion with the aim to relate past
 environmental change with demography and diversification dynamics. The
 here provided datasets include 1) an alignment of a fragment of the gene
 cytochrome c oxidase subunit 1 for 60 specimens; 2) genotype data for 321
 individuals from 39 localities for 15 microsatellite loci (total dataset);
 3) a regrouped genotype dataset (282 specimens from 21 localities for 15
 microsatellite loci), which was used for some analyses in our study.
 Analyses were performed with various programs as reported in our paper.
 Here we provide input and result files (33 files in total) for these
 analyses, complemented with an R script that readily allows reproducing
 the majority of our inquiries.
description: Methodological information is provided in the material and
 methods section of to “Population collapse in viviparid gastropods of the
 Lake Victoria ecoregion started before the Last Glacial Maximum” by Van
 Bocxlaer et al. and accepted for publication in Molecular
 Ecology.
description: Here we briefly describe
 the 33 provided files. In what follows we refer to “Population collapse in
 viviparid gastropods of the Lake Victoria ecoregion started before the
 Last Glacial Maximum” by Van Bocxlaer et al. and accepted for publication
 in Molecular Ecology as “our paper”. <i>VanBocxlaer_et_al_Bellamya_LakeVictoriaEcoregion_Rscripts.R</i> This file contains all R scripts used to analyze the viviparid data with respect to the aims of our paper. It also contains full session information to maximize reproducibility. <i>Bellamya_COI_alignment.fas</i> This dataset contains the alignment for 60 specimens of a 593-basepair fragment of the gene cytochrome <i>c</i> oxidase subunit 1 that was used for network analysis as reported in our paper. The sequence codes provided in this file correspond to those reported in Table S1 of our paper. This table also links to NCBI GenBank accession numbers. <i>Bellamya_COI_labels.txt</i> This file links the sequence codes used in the abovementioned fasta file with population codes and information on waterbodies/genetic clusters as reported in our paper. <i>Bellamya_TotalDataset_RawGenotypeTable.txt</i> This file contains our total microsatellite dataset, i.e. genotype information for 321 individuals from 39 localities genotyped for 15 microsatellite loci, in simple tabular format. <i>Bellamya_TotalDataset_GENEPOP.txt</i> This file contains our total microsatellite dataset, i.e. genotype information for 321 individuals from 39 localities genotyped for 15 microsatellite loci, in GenePop format, which allows easy conversion to other formats. <i>Bellamya_RegroupedDataset_GENEPOP.txt</i> This file contains our regrouped microsatellite dataset, i.e. genotype information for 282 individuals from 21 localities genotyped for 15 microsatellite loci, in GenePop format, which allows easy conversion to other formats. <i>Bellamya_TotalData_GenePop_input.gen</i> Most analyses reported in our paper have been performed with 13 instead of 15 microsatellite loci, because two loci displayed substantial evidence for null alleles. This file contains our total microsatellite dataset excluding these two microsatellite loci (Bel_L50 and Bel_L45). The data are provided in GenePop format, which allows easy conversion to other formats. <i>Bellamya_TotalData_GenePop_input_II.gen</i> This file is identical to “<i>Bellamya_TotalData_GenePop_input.gen</i>”, however, with null alleles coded as ‘999’ instead of ‘000’. <i>Bellamya_TotalData_MorphoID.txt</i> This file contains morphological identifications for all specimens in our total microsatellite dataset (with rows in the same order as in <i>Bellamya_TotalData_GenePop_input.gen</i>). A more complete overview of these data is provided in Table S2 of our paper. <i>Bellamya_RegroupedData_GenePop_input.gen </i> Most analyses reported in our paper have been performed with 13 instead of 15 microsatellite loci, because two loci displayed substantial evidence for null alleles. This file contains our regrouped microsatellite dataset excluding these two microsatellite loci (Bel_L50 and Bel_L45). The data are provided in GenePop format, which allows easy conversion to other formats. <i>Bellamya_RegroupedData_GenePop_input_II.gen</i> This file is identical to “<i>Bellamya_ RegroupedData_GenePop_input.gen</i>”, however, with null alleles coded as ‘999’ instead of ‘000’. <i>Bellamya_RegroupedData_3POP_GenePop_input.gen</i> This file contains our regrouped microsatellite dataset (321 individuals from 39 localities genotyped for 13 microsatellite loci), now organized in 3 gene pools following our results of analyses of spatial patterns of genetic structure. <i>Bellamya_All_Loci_LD.txt</i> This file contains our data on pairwise linkage disequilibrium among loci. The associated <i>p</i>-values are not corrected yet for multiple tests (such a correction is performed in our R script). <i>Bellaymya_Fstat_Fst_13loci.txt</i> This file contains pairwise population <i>F<sub>ST</sub></i> values for our regrouped dataset (282 individuals from 21 localities genotyped for 13 microsatellite loci) in tabular format. <i>Bellamya_Pairwise_FST_pvalues.txt</i> This file contains pairwise population <i>F<sub>ST</sub></i> values as a list, with the associated <i>p</i>-values. The associated <i>p</i>-values are not corrected yet for multiple tests (such a correction is performed in our R script). Bellamya_Pairwise_FST_pvalues_geodist.txt This file contains pairwise population <i>F<sub>ST</sub></i> values as a list, with the associated <i>p</i>-values (corrected for the false discovery rate) and the geographic distance between the localities of each pairwise comparison. <i>Bellamya_LocalityData.txt</i> This file contains Table S5 of our manuscript together with some additional information for the purpose of preparing figures in R. <i>Bellamya_LVG_Fi_data.txt</i> This file contains data on individual inbreeding coefficients (F<sub>i</sub>) for localities within the Lake Victoria gene pool with associated information on how many morphospecies were found at each locality (‘M’ for multiple, ‘S’ for single). This information has been extracted from Table S5 of our paper. <i>STRUCTURE_39pop_321ind_13loci_InputData.txt</i> This file contains our total microsatellite dataset in the format required by the software STRUCTURE. <i>Structure3Groups_assignementFile_Indivs_FINAL.txt</i> This file contains the assignment probabilities for each individual in our total microsatellite dataset, as obtained from STRUCTURE analysis with 3 genetic clusters. <i>SPAGeDi-RegroupedDataset_AllPops_GeneticInputData.txt</i> This file contains our regrouped microsatellite dataset (282 individuals from 21 localities, genotyped for 13 microsatellite loci) in the format required for analyses with the software SPAGeDi. <i>SPAGeDi-RegroupedDataset_AllPops_GeoDistances_indiv.txt</i> This file contains geographical distances for each pairwise comparison between two individuals in our regrouped microsatellite dataset (282 individuals from 21 localities, genotyped for 13 microsatellite loci), used by the software SPAGeDi. <i>SPAGeDi-RegroupedDataset_LVG_GeneticInputData.txt</i> This file contains the genotypic input dataset for analyses with the software SPAGeDi on the populations of the Lake Victoria gene pool (212 individuals from 16 localities, genotyped for 13 microsatellite loci). <i>SPAGeDi-RegroupedDataset_LVG_GeoDistances_indiv.txt</i> This file contains geographical distances for each pairwise comparison between two individuals in the Lake Victoria gene pool from our regrouped dataset (212 individuals from 16 localities, genotyped for 13 microsatellite loci), used by the software SPAGeDi. <i>Bellamya_SPAGeDi_Total_Output_GeoDist_13loci.txt</i> Output summary of all our analyses with the software SPAGeDi. <i>Bellamya_InputBayesass_migration_13loci.txt</i> This file contains the genetic input data for conversion to the input format required by the software BayesAss. <i>Bayesass_InputFile.txt</i> This file contains the genetic input data in the format required by the software BayesAss. Results of the BayesAss analyses are reported in Table S7 of our paper. <i>LVictoria_pop1_Final.mss</i> Input dataset for <i>Bellamya</i> gastropods of the Lake Victoria gene pool for analysis with the software DIYABC. Given that DIYABC treats null alleles as missing data all 15 microsatellite loci were included. This dataset is a subset compiled from the regrouped dataset. <i>LVictoria_pop1_Final_noADM.mss</i> Input dataset for <i>Bellamya</i> gastropods of the Lake Victoria gene pool for analysis with the software DIYABC, excluding individuals that show signs of admixture. Given that DIYABC treats null alleles as missing data all 15 microsatellite loci were included. This dataset is a subset compiled from the regrouped dataset. <i>VictoriaNile_pop1_Final.mss</i> Input dataset for <i>Bellamya</i> gastropods of the Lake Kyoga gene pool for analysis with the software DIYABC. Given that DIYABC treats null alleles as missing data all 15 microsatellite loci were included. This dataset is a subset compiled from the regrouped dataset. <i>LAlbertWhiteNile_pop1_Final.mss</i> Input dataset for <i>Bellamya gastropods</i> of the Lake Albert gene pool for analysis with the software DIYABC. Given that DIYABC treats null alleles as missing data all 15 microsatellite loci were included. This dataset is a subset compiled from the regrouped dataset. <i>LAlbertWhiteNile_pop1_Final_noADM.mss</i> Input dataset for <i>Bellamya gastropods</i> of the Lake Albert gene pool for analysis with the software DIYABC, excluding individuals that show signs of admixture. Given that DIYABC treats null alleles as missing data all 15 microsatellite loci were included. This dataset is a subset compiled from the regrouped dataset. <i>ABC_robustness_migration.txt</i> This file contains the results of our analysis as to how migration affects approximate Bayesian computations. Each row represents the results for a single simulated dataset (33,000 datasets in total), which are summarized in Figure S7 of our paper.