Snail species, abundance and physicochemical water factors by contact site and season.

Source data used in the publication: Berger et al. (2021) - Provisional title: Whole-genome sequencing of Schistosoma mansoni reveals extensive diversity with limited selection despite repeated mass drug administration. These data were used to generate all figures used in the publication and all files are organised and labelled specifically to run with the custom code that uses these data can be found at: DOI:10.5281/zenodo.4975909. (I'm happy to answer any questions just email [email protected], manuscript that goes with these data is about to be resubmitted) File descriptions:
prunedData.eigenvec - Figure 2a&b - PCA eigenvectors prunedData.eigenval - Figure 2a&b- PCA eigenvalues autosomes.mdist - Figure 2c - PLINK distance matrix used to build the neighbour joining phylogeny all.pi.pixy.schools.txt - Figure 2d - Nucleotide diversity results for each school subpopulation autosomes.dxy.5kb.schools.txt - Figure 2e - Autosomal D_XY results between school subpopulations. autosomes.fst.5kb.schools.txt - Figure 2e - Autosomal F_ST results between school subpopulations. admixture_all.txt - Figure 2f - ADMIXTURE results for each sample and population sizes, column 1 represents number of populations (K), columns 3-8 represent admixture values for each population. sfs.csv - Figure 3a, Supplementary figure 10a - Site frequency spectra (allelic proportions at each frequency bin) for each school. TD.all.txt - Figure 3b - Tajima's D values calculated in 5 kb windows for each school subpopulation. ALL.MAYUGE.IHS.ihs.out.100bins.norm.txt.zip - Figure 4a - Normalised iHS scores for the Mayuge district parasite populations (Selscan output). ALL.TORORO.IHS.ihs.out.100bins.norm.txt.zip - Figure 4b - Normalised iHS scores for the Tororo district parasite populations (Selscan output). ALL.MAYUGEvsTORORO.xpehh.xpehh.out.norm.txt.zip - Figure 4c - Normalised XP-EHH scores between Mayuge and Tororo parasite populations. MAYUGE_TORORO_2000.windowed.weir.txt.zip - Figure 4d - F_ST values calculated between Mayuge and Tororo populations in 2kb windows. MAYUGE_PI.windowed.pi.zip - Figure 4e, Supplementary figure 12a&c - Nucleotide diversity values calculated in 2 kb windows for Mayuge populations.
TORORO_PI.windowed.pi.zip - Figure 4e, Supplementary figure 12a&c- Nucleotide diversity values calculated in 2 kb windows for Kocoge populations (Tororo district). all.pi.treat.fix.txt.zip - Figure 5a - Nucleotide diversity results for each treatment subpopulation autosomes.dxy.5kb.treatment.txt - Figure 5b - Autosomal D_XY results between clearance phenotype subpopulations. autosomes.fst.5kb.treatment.txt - Figure 5b - Autosomal F_ST results between clearance phenotype subpopulations. fst.windows.2kb.treatment.txt.zip - Figure 5c - F_ST values for comparisons between different treatment groups (Pre-treatment, post-treatment (good clearers), post-treatment (poor clearers)) assoc_err_binary.txt.zip - Figure 5d - Results of binary trait association between miracidia sampled from hosts with good clearance phenotypes (where treatment appeared to be highly effective) and miracidia isolated post-treatment from hosts with poor clearance phenotypes (where miracidia are potentially derived from parasites that survived treatment. assoc_err_linear.txt.zip - Figure 5e - Results of linear regression genome-wide association study with the ERR estimates for all 198 samples, using the mean of the posterior ERR estimates from Crellen et al. (2016) as a quantitative trait. median.coverage.txt - Supplementary figure 1 - Normalised depth of read coverage (column 4) calculated in 25 kb windows (columns 2&3) across all samples for all chromosomes (column 1).
cohort.genotyped.txt.zip - Supplementary figure 2(a-f) - Variant quality site values (used to inform variant site retention or removal).
hard_filtered.imiss.txt - Supplementary figure 2g - Per sample variant missingness (used to inform quality control)
hard_filtered_filtindv.lmiss.txt.zip - Supplementary figure 2h - Per site missingness (used to inform quality control) cv_scores.txt - Supplementary figure 5 - ADMIXTURE CV scores (column 2) for each population size (1). *_SMC_SE.csv - Supplementary figure 6 -SMC++ results (from 25 subsampled replicates) for each school subpopulation and outgroup samples. smcpp.csv - Supplementary figure 7 - SMC++ results for each school subpopulation and outgroup samples.
pi.per_host.txt.zip - Supplementary figure 8a-d - Nucleotide diversity values for each host infrapopulation.
sexing.csv - Supplementary figure 9- inferred sex (based on differential read coverage over pseudoautosomal and Z-specific regions of the Z chromosome). MAYUGE_TAJIMA_D.Tajima.D.2kb.txt.zip - Supplementary figure 11 - Tajima's D values calculated for the Mayuge population in 2kb windows. Tororo_TAJIMA_D.Tajima.D.2kb.txt.zip - Supplementary figure 11 - Tajima's D values calculated for the Tororo population in 2kb windows. genes.bed - Supplementary figures 13-18 - Coordinates of gene models (S. mansoni v7 annotation). KOCOGE_SITE_PI.sites.pi.txt.zip - Supplementary figures 13-18 - Per site nucleotide diversity values MAYUGE_TORORO_sites.weir.fst.txt.zip - Supplementary figures 13a-18a - Per site F_ST values between Mayuge and Tororo populations. kocoge_median.ld.txt.zip - Supplementary figure 19 - The decay of linkage disequilibrium with genomic distance between all sites within 50 kb for the Kocoge parasite samples. Chromosomes are shown in column 1, distance in column 2, median values in column 3. mayuge_median.ld.txt.zip - Supplementary figure 19 - The decay of linkage disequilibrium with genomic distance between all sites within 50 kb for the Mayuge parasite samples. Chromosomes are shown in column 1, distance in column 2, median values in column 3. median.sample.coverage.txt - Median chromosomal depth of read coverage for each sample. schools.list - List of samples and schools where they were sampled. coverage_5kb.windows.txt.zip - Per sample depth of read coverage in 5 kb windows. Columns 4,5,6 represent the median, mean and sstev of coverage for each 5kb window (columns 2&3) along each chromosome (column 1).

Source data used in the publication: Berger et al. (2021) - Provisional title: 'Whole-genome sequencing of Schistosoma mansoni reveals extensive diversity with limited selection despite mass drug administration'. These data were used to generate all figures used in the publication and all files are organised and labelled specifically to run with the custom code that uses these data can be found at: http://doi.org/10.5281/zenodo.4975908.
File descriptions: SOURCE DATA.zip - All source data for all figures. Figure 1b: supplementary_data_9.txt - Metadata Figure 2a&b: 207_PCA.eigenvec - PCA eigenvectors 207_PCA.eigenval - PCA eigenvalues Figure 2c: autosomes.mdist - PLINK distance matrix used to build the neighbour joining phylogeny Figure 2d: all.pi.pixy.schools.txt - Nucleotide diversity results for each school subpopulation. Figure 2e: autosomes.dxy.5kb.schools.txt - Autosomal D_XY results between school subpopulations. autosomes.fst.5kb.schools.txt - Autosomal F_ST results between school subpopulations. Figure 2f: admixture_all.txt - ADMIXTURE results for each sample and population sizes, column 1 represents number of populations (K), columns 3-8 represent admixture values for each population. Figure 3a, Supplementary figure 10a: sfs.csv - Site frequency spectra (allelic proportions at each frequency bin) for each school. Figure 3b: TD.all.txt - Tajima's D values calculated in 5 kb windows for each school subpopulation. Figure 4a, Supplementary figures 13-18: ALL.MAYUGE.IHS.ihs.out.100bins.norm.txt.zip - Normalised iHS scores for the Mayuge district parasite populations (Selscan output). Figure 4b, Supplementary figures 13-18: ALL.TORORO.IHS.ihs.out.100bins.norm.txt.zip - - Normalised iHS scores for the Tororo district parasite populations (Selscan output). Figure 4c, Supplementary figures 13-18: ALL.MAYUGEvsTORORO.xpehh.xpehh.out.norm.txt.zip - - Normalised XP-EHH scores between Mayuge and Tororo parasite populations. Figure 4d, Supplementary figures 13-18: MAYUGE_TORORO_2000.windowed.weir.txt.zip - F_ST values calculated between Mayuge and Tororo populations in 2kb windows. Figure 4e, Supplementary figures 12a&c: MAYUGE_PI.windowed.pi.zip - Nucleotide diversity values calculated in 2 kb windows for Mayuge populations. TORORO_PI.windowed.pi.zip - Nucleotide diversity values calculated in 2 kb windows for Kocoge populations (Tororo district). Figure 5a: all.pi.treat.fix.txt.zip - Nucleotide diversity results for each treatment subpopulation Figure 5b autosomes.dxy.5kb.treatment.txt - - Autosomal D_XY results between clearance phenotype subpopulations. autosomes.fst.5kb.treatment.txt - Autosomal F_ST results between clearance phenotype subpopulations. Figure 5c: fst.windows.2kb.treatment.txt.zip - F_ST values for comparisons between different treatment groups (Pre-treatment, post-treatment (good clearers), post-treatment (poor clearers)) Figure 5d: assoc_err_binary.txt.zip - Results of binary trait association between miracidia sampled from hosts with good clearance phenotypes (where treatment appeared to be highly effective) and miracidia isolated post-treatment from hosts with poor clearance phenotypes (where miracidia are potentially derived from parasites that survived treatment. Figure 5e: assoc_err_linear.txt.zip - - Results of linear regression genome-wide association study with the ERR estimates for all 198 samples, using the mean of the posterior ERR estimates from Crellen et al. (2016) as a quantitative trait. Supplementary figure 1: median.coverage.txt - Normalised depth of read coverage (column 4) calculated in 25 kb windows (columns 2&3) across all samples for all chromosomes (column 1). Supplementary figure 2a-f: cohort.genotyped.txt.zip - - Variant quality site values (used to inform variant site retention or removal). Supplementary figure 2g: hard_filtered.imiss.txt - Per sample variant missingness (used to inform quality control). Supplementary figure 2h: hard_filtered_filtindv.lmiss.txt.zip - Per site missingness (used to inform quality control). Supplementary figure 3a, 4a, 4b: prunedData.eigenvec - PCA eigenvectors prunedData.eigenval - PCA eigenvalues Supplementary figure 3b: pruned_data.mdist.csv - Distance matrix used as the basis for the neighbour joining phylogeny. Supplementary figure 5: cv_scores.txt - ADMIXTURE coefficient of variation scores (column 2) for each population size (1). Supplementary figure 6: *_SMC_SE.csv - SMC++ results (from 25 subsampled replicates) for each school subpopulation and outgroup samples. Supplementary Figure 7: smcpp.csv - SMC++ results for each school subpopulation and outgroup samples. Supplementary Figure 8a-d pi.per_host.txt.zip - Nucleotide diversity values for each host infrapopulation. Supplementary Figure 9: sexing.csv - inferred sex (based on differential read coverage over pseudoautosomal and Z-specific regions of the Z chromosome). Supplementary Figure 10b: sfs_res.csv - residuals for the SFS analysis in 3a/10a. Supplementary Figure 11: MAYUGE_TAJIMA_D.Tajima.D.2kb.txt.zip - Tajima's D values calculated for the Mayuge population in 2kb windows. Tororo_TAJIMA_D.Tajima.D.2kb.txt.zip - Tajima's D values calculated for the Tororo population in 2kb windows. Supplementary Figures 13-18: genes.bed - Coordinates of gene models (S. mansoni v7 annotation). KOCOGE_SITE_PI.sites.pi.txt.zip - Per site nucleotide diversity values MAYUGE_TORORO_sites.weir.fst.txt.zip - Per site F_ST values between Mayuge and Tororo populations. coverage_5kb.windows.txt.zip - Per sample depth of read coverage in 5 kb windows. Columns 4,5,6 represent the median, mean and sstev of coverage for each 5kb window (columns 2&3) along each chromosome (column 1). median.sample.coverage.txt - Median chromosomal depth of read coverage for each sample. Supplementary Figure 19: kocoge_median.ld.txt.zip - - The decay of linkage disequilibrium with genomic distance between all sites within 50 kb for the Kocoge parasite samples. Chromosomes are shown in column 1, distance in column 2, median values in column 3. mayuge_median.ld.txt.zip - The decay of linkage disequilibrium with genomic distance between all sites within 50 kb for the Mayuge parasite samples. Chromosomes are shown in column 1, distance in column 2, median values in column 3. Misc files: schools.list - List of samples and schools where they were sampled.

Source data used in the publication: Berger et al. (2021) - Provisional title: Whole-genome sequencing of Schistosoma mansoni reveals extensive diversity with limited selection despite repeated mass drug administration. These data were used to generate all figures used in the publication and all files are organised and labelled specifically to run with the custom code that uses these data can be found at: http://doi.org/10.5281/zenodo.4975908. (I'm happy to answer any questions just email [email protected])
File descriptions: Figure 2a&b: 207_PCA.eigenvec - PCA eigenvectors 207_PCA.eigenval - PCA eigenvalues Figure 2c: autosomes.mdist - PLINK distance matrix used to build the neighbour joining phylogeny Figure 2d: all.pi.pixy.schools.txt - Nucleotide diversity results for each school subpopulation. Figure 2e: autosomes.dxy.5kb.schools.txt - Autosomal D_XY results between school subpopulations. autosomes.fst.5kb.schools.txt - Autosomal F_ST results between school subpopulations. Figure 2f: admixture_all.txt - ADMIXTURE results for each sample and population sizes, column 1 represents number of populations (K), columns 3-8 represent admixture values for each population. Figure 3a, Supplementary figure 10a: sfs.csv - Site frequency spectra (allelic proportions at each frequency bin) for each school. Figure 3b: TD.all.txt - Tajima's D values calculated in 5 kb windows for each school subpopulation. Figure 4a, Supplementary figures 13-18: ALL.MAYUGE.IHS.ihs.out.100bins.norm.txt.zip - Normalised iHS scores for the Mayuge district parasite populations (Selscan output). Figure 4b, Supplementary figures 13-18: ALL.TORORO.IHS.ihs.out.100bins.norm.txt.zip - - Normalised iHS scores for the Tororo district parasite populations (Selscan output). Figure 4c, Supplementary figures 13-18: ALL.MAYUGEvsTORORO.xpehh.xpehh.out.norm.txt.zip - - Normalised XP-EHH scores between Mayuge and Tororo parasite populations. Figure 4d, Supplementary figures 13-18: MAYUGE_TORORO_2000.windowed.weir.txt.zip - F_ST values calculated between Mayuge and Tororo populations in 2kb windows. Figure 4e, Supplementary figures 12a&c: MAYUGE_PI.windowed.pi.zip - Nucleotide diversity values calculated in 2 kb windows for Mayuge populations. TORORO_PI.windowed.pi.zip - Nucleotide diversity values calculated in 2 kb windows for Kocoge populations (Tororo district). Figure 5a: all.pi.treat.fix.txt.zip - Nucleotide diversity results for each treatment subpopulation Figure 5b autosomes.dxy.5kb.treatment.txt - - Autosomal D_XY results between clearance phenotype subpopulations. autosomes.fst.5kb.treatment.txt - Autosomal F_ST results between clearance phenotype subpopulations. Figure 5c: fst.windows.2kb.treatment.txt.zip - F_ST values for comparisons between different treatment groups (Pre-treatment, post-treatment (good clearers), post-treatment (poor clearers)) Figure 5d: assoc_err_binary.txt.zip - Results of binary trait association between miracidia sampled from hosts with good clearance phenotypes (where treatment appeared to be highly effective) and miracidia isolated post-treatment from hosts with poor clearance phenotypes (where miracidia are potentially derived from parasites that survived treatment. Figure 5e: assoc_err_linear.txt.zip - - Results of linear regression genome-wide association study with the ERR estimates for all 198 samples, using the mean of the posterior ERR estimates from Crellen et al. (2016) as a quantitative trait. Supplementary figure 1: median.coverage.txt - Normalised depth of read coverage (column 4) calculated in 25 kb windows (columns 2&3) across all samples for all chromosomes (column 1). Supplementary figure 2a-f: cohort.genotyped.txt.zip - - Variant quality site values (used to inform variant site retention or removal). Supplementary figure 2g: hard_filtered.imiss.txt - Per sample variant missingness (used to inform quality control). Supplementary figure 2h: hard_filtered_filtindv.lmiss.txt.zip - Per site missingness (used to inform quality control). Supplementary figure 3a, 4a, 4b: prunedData.eigenvec - PCA eigenvectors prunedData.eigenval - PCA eigenvalues Supplementary figure 3b: pruned_data.mdist.csv - Distance matrix used as the basis for the neighbour joining phylogeny. Supplementary figure 5: cv_scores.txt - ADMIXTURE coefficient of variation scores (column 2) for each population size (1). Supplementary figure 6: *_SMC_SE.csv - SMC++ results (from 25 subsampled replicates) for each school subpopulation and outgroup samples. Supplementary figure 7: smcpp.csv - SMC++ results for each school subpopulation and outgroup samples. Supplementary figure 8a-d pi.per_host.txt.zip - Nucleotide diversity values for each host infrapopulation. Supplementary figure 9: sexing.csv - inferred sex (based on differential read coverage over pseudoautosomal and Z-specific regions of the Z chromosome). Supplementary figure 11: MAYUGE_TAJIMA_D.Tajima.D.2kb.txt.zip - Tajima's D values calculated for the Mayuge population in 2kb windows. Tororo_TAJIMA_D.Tajima.D.2kb.txt.zip - Tajima's D values calculated for the Tororo population in 2kb windows. Supplementary figures 13-18: genes.bed - Coordinates of gene models (S. mansoni v7 annotation). KOCOGE_SITE_PI.sites.pi.txt.zip - Per site nucleotide diversity values MAYUGE_TORORO_sites.weir.fst.txt.zip - Per site F_ST values between Mayuge and Tororo populations. coverage_5kb.windows.txt.zip - Per sample depth of read coverage in 5 kb windows. Columns 4,5,6 represent the median, mean and sstev of coverage for each 5kb window (columns 2&3) along each chromosome (column 1). median.sample.coverage.txt - Median chromosomal depth of read coverage for each sample. Supplementary figure 19: kocoge_median.ld.txt.zip - - The decay of linkage disequilibrium with genomic distance between all sites within 50 kb for the Kocoge parasite samples. Chromosomes are shown in column 1, distance in column 2, median values in column 3. mayuge_median.ld.txt.zip - The decay of linkage disequilibrium with genomic distance between all sites within 50 kb for the Mayuge parasite samples. Chromosomes are shown in column 1, distance in column 2, median values in column 3. Misc files: schools.list - List of samples and schools where they were sampled.

Source data used in the publication: Berger et al. (2021) - Provisional title: 'Population genomics of Schistosoma mansoni reveals extensive diversity with limited selection despite repeated mass drug administration'. These data were used to generate all figures used in the publication and all files are organised and labelled specifically to run with the custom code that uses these data can be found at: http://doi.org/10.5281/zenodo.4975908. (I'm happy to answer any questions just email [email protected])
File descriptions: Figure 1b: supplementary_data_9.txt - Metadata Figure 2a&b: 207_PCA.eigenvec - PCA eigenvectors 207_PCA.eigenval - PCA eigenvalues Figure 2c: autosomes.mdist - PLINK distance matrix used to build the neighbour joining phylogeny Figure 2d: all.pi.pixy.schools.txt - Nucleotide diversity results for each school subpopulation. Figure 2e: autosomes.dxy.5kb.schools.txt - Autosomal D_XY results between school subpopulations. autosomes.fst.5kb.schools.txt - Autosomal F_ST results between school subpopulations. Figure 2f: admixture_all.txt - ADMIXTURE results for each sample and population sizes, column 1 represents number of populations (K), columns 3-8 represent admixture values for each population. Figure 3a, Supplementary figure 10a: sfs.csv - Site frequency spectra (allelic proportions at each frequency bin) for each school. Figure 3b: TD.all.txt - Tajima's D values calculated in 5 kb windows for each school subpopulation. Figure 4a, Supplementary figures 13-18: ALL.MAYUGE.IHS.ihs.out.100bins.norm.txt.zip - Normalised iHS scores for the Mayuge district parasite populations (Selscan output). Figure 4b, Supplementary figures 13-18: ALL.TORORO.IHS.ihs.out.100bins.norm.txt.zip - - Normalised iHS scores for the Tororo district parasite populations (Selscan output). Figure 4c, Supplementary figures 13-18: ALL.MAYUGEvsTORORO.xpehh.xpehh.out.norm.txt.zip - - Normalised XP-EHH scores between Mayuge and Tororo parasite populations. Figure 4d, Supplementary figures 13-18: MAYUGE_TORORO_2000.windowed.weir.txt.zip - F_ST values calculated between Mayuge and Tororo populations in 2kb windows. Figure 4e, Supplementary figures 12a&c: MAYUGE_PI.windowed.pi.zip - Nucleotide diversity values calculated in 2 kb windows for Mayuge populations. TORORO_PI.windowed.pi.zip - Nucleotide diversity values calculated in 2 kb windows for Kocoge populations (Tororo district). Figure 5a: all.pi.treat.fix.txt.zip - Nucleotide diversity results for each treatment subpopulation Figure 5b autosomes.dxy.5kb.treatment.txt - - Autosomal D_XY results between clearance phenotype subpopulations. autosomes.fst.5kb.treatment.txt - Autosomal F_ST results between clearance phenotype subpopulations. Figure 5c: fst.windows.2kb.treatment.txt.zip - F_ST values for comparisons between different treatment groups (Pre-treatment, post-treatment (good clearers), post-treatment (poor clearers)) Figure 5d: assoc_err_binary.txt.zip - Results of binary trait association between miracidia sampled from hosts with good clearance phenotypes (where treatment appeared to be highly effective) and miracidia isolated post-treatment from hosts with poor clearance phenotypes (where miracidia are potentially derived from parasites that survived treatment. Figure 5e: assoc_err_linear.txt.zip - - Results of linear regression genome-wide association study with the ERR estimates for all 198 samples, using the mean of the posterior ERR estimates from Crellen et al. (2016) as a quantitative trait. Supplementary figure 1: median.coverage.txt - Normalised depth of read coverage (column 4) calculated in 25 kb windows (columns 2&3) across all samples for all chromosomes (column 1). Supplementary figure 2a-f: cohort.genotyped.txt.zip - - Variant quality site values (used to inform variant site retention or removal). Supplementary figure 2g: hard_filtered.imiss.txt - Per sample variant missingness (used to inform quality control). Supplementary figure 2h: hard_filtered_filtindv.lmiss.txt.zip - Per site missingness (used to inform quality control). Supplementary figure 3a, 4a, 4b: prunedData.eigenvec - PCA eigenvectors prunedData.eigenval - PCA eigenvalues Supplementary figure 3b: pruned_data.mdist.csv - Distance matrix used as the basis for the neighbour joining phylogeny. Supplementary figure 5: cv_scores.txt - ADMIXTURE coefficient of variation scores (column 2) for each population size (1). Supplementary figure 6: *_SMC_SE.csv - SMC++ results (from 25 subsampled replicates) for each school subpopulation and outgroup samples. Supplementary Figure 7: smcpp.csv - SMC++ results for each school subpopulation and outgroup samples. Supplementary Figure 8a-d pi.per_host.txt.zip - Nucleotide diversity values for each host infrapopulation. Supplementary Figure 9: sexing.csv - inferred sex (based on differential read coverage over pseudoautosomal and Z-specific regions of the Z chromosome). Supplementary Figure 10b: sfs_res.csv - residuals for the SFS analysis in 3a/10a. Supplementary Figure 11: MAYUGE_TAJIMA_D.Tajima.D.2kb.txt.zip - Tajima's D values calculated for the Mayuge population in 2kb windows. Tororo_TAJIMA_D.Tajima.D.2kb.txt.zip - Tajima's D values calculated for the Tororo population in 2kb windows. Supplementary Figures 13-18: genes.bed - Coordinates of gene models (S. mansoni v7 annotation). KOCOGE_SITE_PI.sites.pi.txt.zip - Per site nucleotide diversity values MAYUGE_TORORO_sites.weir.fst.txt.zip - Per site F_ST values between Mayuge and Tororo populations. coverage_5kb.windows.txt.zip - Per sample depth of read coverage in 5 kb windows. Columns 4,5,6 represent the median, mean and sstev of coverage for each 5kb window (columns 2&3) along each chromosome (column 1). median.sample.coverage.txt - Median chromosomal depth of read coverage for each sample. Supplementary Figure 19: kocoge_median.ld.txt.zip - - The decay of linkage disequilibrium with genomic distance between all sites within 50 kb for the Kocoge parasite samples. Chromosomes are shown in column 1, distance in column 2, median values in column 3. mayuge_median.ld.txt.zip - The decay of linkage disequilibrium with genomic distance between all sites within 50 kb for the Mayuge parasite samples. Chromosomes are shown in column 1, distance in column 2, median values in column 3. Misc files: schools.list - List of samples and schools where they were sampled.

Source data used in the publication: Berger et al. (2021) - Provisional title: 'Whole-genome sequencing of Schistosoma mansoni reveals extensive diversity with limited selection despite mass drug administration'. These data were used to generate all figures used in the publication and all files are organised and labelled specifically to run with the custom code that uses these data can be found at: http://doi.org/10.5281/zenodo.4975908. (I'm happy to answer any questions just email [email protected])
File descriptions: source_data.xlsx - All source data for all figures. Figure 1b: supplementary_data_9.txt - Metadata Figure 2a&b: 207_PCA.eigenvec - PCA eigenvectors 207_PCA.eigenval - PCA eigenvalues Figure 2c: autosomes.mdist - PLINK distance matrix used to build the neighbour joining phylogeny Figure 2d: all.pi.pixy.schools.txt - Nucleotide diversity results for each school subpopulation. Figure 2e: autosomes.dxy.5kb.schools.txt - Autosomal D_XY results between school subpopulations. autosomes.fst.5kb.schools.txt - Autosomal F_ST results between school subpopulations. Figure 2f: admixture_all.txt - ADMIXTURE results for each sample and population sizes, column 1 represents number of populations (K), columns 3-8 represent admixture values for each population. Figure 3a, Supplementary figure 10a: sfs.csv - Site frequency spectra (allelic proportions at each frequency bin) for each school. Figure 3b: TD.all.txt - Tajima's D values calculated in 5 kb windows for each school subpopulation. Figure 4a, Supplementary figures 13-18: ALL.MAYUGE.IHS.ihs.out.100bins.norm.txt.zip - Normalised iHS scores for the Mayuge district parasite populations (Selscan output). Figure 4b, Supplementary figures 13-18: ALL.TORORO.IHS.ihs.out.100bins.norm.txt.zip - - Normalised iHS scores for the Tororo district parasite populations (Selscan output). Figure 4c, Supplementary figures 13-18: ALL.MAYUGEvsTORORO.xpehh.xpehh.out.norm.txt.zip - - Normalised XP-EHH scores between Mayuge and Tororo parasite populations. Figure 4d, Supplementary figures 13-18: MAYUGE_TORORO_2000.windowed.weir.txt.zip - F_ST values calculated between Mayuge and Tororo populations in 2kb windows. Figure 4e, Supplementary figures 12a&c: MAYUGE_PI.windowed.pi.zip - Nucleotide diversity values calculated in 2 kb windows for Mayuge populations. TORORO_PI.windowed.pi.zip - Nucleotide diversity values calculated in 2 kb windows for Kocoge populations (Tororo district). Figure 5a: all.pi.treat.fix.txt.zip - Nucleotide diversity results for each treatment subpopulation Figure 5b autosomes.dxy.5kb.treatment.txt - - Autosomal D_XY results between clearance phenotype subpopulations. autosomes.fst.5kb.treatment.txt - Autosomal F_ST results between clearance phenotype subpopulations. Figure 5c: fst.windows.2kb.treatment.txt.zip - F_ST values for comparisons between different treatment groups (Pre-treatment, post-treatment (good clearers), post-treatment (poor clearers)) Figure 5d: assoc_err_binary.txt.zip - Results of binary trait association between miracidia sampled from hosts with good clearance phenotypes (where treatment appeared to be highly effective) and miracidia isolated post-treatment from hosts with poor clearance phenotypes (where miracidia are potentially derived from parasites that survived treatment. Figure 5e: assoc_err_linear.txt.zip - - Results of linear regression genome-wide association study with the ERR estimates for all 198 samples, using the mean of the posterior ERR estimates from Crellen et al. (2016) as a quantitative trait. Supplementary figure 1: median.coverage.txt - Normalised depth of read coverage (column 4) calculated in 25 kb windows (columns 2&3) across all samples for all chromosomes (column 1). Supplementary figure 2a-f: cohort.genotyped.txt.zip - - Variant quality site values (used to inform variant site retention or removal). Supplementary figure 2g: hard_filtered.imiss.txt - Per sample variant missingness (used to inform quality control). Supplementary figure 2h: hard_filtered_filtindv.lmiss.txt.zip - Per site missingness (used to inform quality control). Supplementary figure 3a, 4a, 4b: prunedData.eigenvec - PCA eigenvectors prunedData.eigenval - PCA eigenvalues Supplementary figure 3b: pruned_data.mdist.csv - Distance matrix used as the basis for the neighbour joining phylogeny. Supplementary figure 5: cv_scores.txt - ADMIXTURE coefficient of variation scores (column 2) for each population size (1). Supplementary figure 6: *_SMC_SE.csv - SMC++ results (from 25 subsampled replicates) for each school subpopulation and outgroup samples. Supplementary Figure 7: smcpp.csv - SMC++ results for each school subpopulation and outgroup samples. Supplementary Figure 8a-d pi.per_host.txt.zip - Nucleotide diversity values for each host infrapopulation. Supplementary Figure 9: sexing.csv - inferred sex (based on differential read coverage over pseudoautosomal and Z-specific regions of the Z chromosome). Supplementary Figure 10b: sfs_res.csv - residuals for the SFS analysis in 3a/10a. Supplementary Figure 11: MAYUGE_TAJIMA_D.Tajima.D.2kb.txt.zip - Tajima's D values calculated for the Mayuge population in 2kb windows. Tororo_TAJIMA_D.Tajima.D.2kb.txt.zip - Tajima's D values calculated for the Tororo population in 2kb windows. Supplementary Figures 13-18: genes.bed - Coordinates of gene models (S. mansoni v7 annotation). KOCOGE_SITE_PI.sites.pi.txt.zip - Per site nucleotide diversity values MAYUGE_TORORO_sites.weir.fst.txt.zip - Per site F_ST values between Mayuge and Tororo populations. coverage_5kb.windows.txt.zip - Per sample depth of read coverage in 5 kb windows. Columns 4,5,6 represent the median, mean and sstev of coverage for each 5kb window (columns 2&3) along each chromosome (column 1). median.sample.coverage.txt - Median chromosomal depth of read coverage for each sample. Supplementary Figure 19: kocoge_median.ld.txt.zip - - The decay of linkage disequilibrium with genomic distance between all sites within 50 kb for the Kocoge parasite samples. Chromosomes are shown in column 1, distance in column 2, median values in column 3. mayuge_median.ld.txt.zip - The decay of linkage disequilibrium with genomic distance between all sites within 50 kb for the Mayuge parasite samples. Chromosomes are shown in column 1, distance in column 2, median values in column 3. Misc files: schools.list - List of samples and schools where they were sampled.

Source data used in the publication: Berger et al. (2021) - Provisional title: 'Whole-genome sequencing of Schistosoma mansoni reveals extensive diversity with limited selection despite mass drug administration'. These data were used to generate all figures used in the publication and all files are organised and labelled specifically to run with the custom code that uses these data can be found at: http://doi.org/10.5281/zenodo.4975908.
File descriptions: SOURCE DATA.zip - All source data for all figures. Figure 1b: supplementary_data_9.txt - Metadata Figure 2a&b: 207_PCA.eigenvec - PCA eigenvectors 207_PCA.eigenval - PCA eigenvalues Figure 2c: autosomes.mdist - PLINK distance matrix used to build the neighbour joining phylogeny Figure 2d: all.pi.pixy.schools.txt - Nucleotide diversity results for each school subpopulation. Figure 2e: autosomes.dxy.5kb.schools.txt - Autosomal D_XY results between school subpopulations. autosomes.fst.5kb.schools.txt - Autosomal F_ST results between school subpopulations. Figure 2f: admixture_all.txt - ADMIXTURE results for each sample and population sizes, column 1 represents number of populations (K), columns 3-8 represent admixture values for each population. Figure 3a, Supplementary figure 10a: sfs.csv - Site frequency spectra (allelic proportions at each frequency bin) for each school. Figure 3b: TD.all.txt - Tajima's D values calculated in 5 kb windows for each school subpopulation. Figure 4a, Supplementary figures 13-18: ALL.MAYUGE.IHS.ihs.out.100bins.norm.txt.zip - Normalised iHS scores for the Mayuge district parasite populations (Selscan output). Figure 4b, Supplementary figures 13-18: ALL.TORORO.IHS.ihs.out.100bins.norm.txt.zip - - Normalised iHS scores for the Tororo district parasite populations (Selscan output). Figure 4c, Supplementary figures 13-18: ALL.MAYUGEvsTORORO.xpehh.xpehh.out.norm.txt.zip - - Normalised XP-EHH scores between Mayuge and Tororo parasite populations. Figure 4d, Supplementary figures 13-18: MAYUGE_TORORO_2000.windowed.weir.txt.zip - F_ST values calculated between Mayuge and Tororo populations in 2kb windows. Figure 4e, Supplementary figures 12a&c: MAYUGE_PI.windowed.pi.zip - Nucleotide diversity values calculated in 2 kb windows for Mayuge populations. TORORO_PI.windowed.pi.zip - Nucleotide diversity values calculated in 2 kb windows for Kocoge populations (Tororo district). Figure 5a: all.pi.treat.fix.txt.zip - Nucleotide diversity results for each treatment subpopulation Figure 5b autosomes.dxy.5kb.treatment.txt - - Autosomal D_XY results between clearance phenotype subpopulations. autosomes.fst.5kb.treatment.txt - Autosomal F_ST results between clearance phenotype subpopulations. Figure 5c: fst.windows.2kb.treatment.txt.zip - F_ST values for comparisons between different treatment groups (Pre-treatment, post-treatment (good clearers), post-treatment (poor clearers)) Figure 5d: assoc_err_binary.txt.zip - Results of binary trait association between miracidia sampled from hosts with good clearance phenotypes (where treatment appeared to be highly effective) and miracidia isolated post-treatment from hosts with poor clearance phenotypes (where miracidia are potentially derived from parasites that survived treatment. Figure 5e: assoc_err_linear.txt.zip - - Results of linear regression genome-wide association study with the ERR estimates for all 198 samples, using the mean of the posterior ERR estimates from Crellen et al. (2016) as a quantitative trait. Supplementary figure 1: median.coverage.txt - Normalised depth of read coverage (column 4) calculated in 25 kb windows (columns 2&3) across all samples for all chromosomes (column 1). Supplementary figure 2a-f: cohort.genotyped.txt.zip - - Variant quality site values (used to inform variant site retention or removal). Supplementary figure 2g: hard_filtered.imiss.txt - Per sample variant missingness (used to inform quality control). Supplementary figure 2h: hard_filtered_filtindv.lmiss.txt.zip - Per site missingness (used to inform quality control). Supplementary figure 3a, 4a, 4b: prunedData.eigenvec - PCA eigenvectors prunedData.eigenval - PCA eigenvalues Supplementary figure 3b: pruned_data.mdist.csv - Distance matrix used as the basis for the neighbour joining phylogeny. Supplementary figure 5: cv_scores.txt - ADMIXTURE coefficient of variation scores (column 2) for each population size (1). Supplementary figure 6: *_SMC_SE.csv - SMC++ results (from 25 subsampled replicates) for each school subpopulation and outgroup samples. Supplementary Figure 7: smcpp.csv - SMC++ results for each school subpopulation and outgroup samples. Supplementary Figure 8a-d pi.per_host.txt.zip - Nucleotide diversity values for each host infrapopulation. Supplementary Figure 9: sexing.csv - inferred sex (based on differential read coverage over pseudoautosomal and Z-specific regions of the Z chromosome). Supplementary Figure 10b: sfs_res.csv - residuals for the SFS analysis in 3a/10a. Supplementary Figure 11: MAYUGE_TAJIMA_D.Tajima.D.2kb.txt.zip - Tajima's D values calculated for the Mayuge population in 2kb windows. Tororo_TAJIMA_D.Tajima.D.2kb.txt.zip - Tajima's D values calculated for the Tororo population in 2kb windows. Supplementary Figures 13-18: genes.bed - Coordinates of gene models (S. mansoni v7 annotation). KOCOGE_SITE_PI.sites.pi.txt.zip - Per site nucleotide diversity values MAYUGE_TORORO_sites.weir.fst.txt.zip - Per site F_ST values between Mayuge and Tororo populations. coverage_5kb.windows.txt.zip - Per sample depth of read coverage in 5 kb windows. Columns 4,5,6 represent the median, mean and sstev of coverage for each 5kb window (columns 2&3) along each chromosome (column 1). median.sample.coverage.txt - Median chromosomal depth of read coverage for each sample. Supplementary Figure 19: kocoge_median.ld.txt.zip - - The decay of linkage disequilibrium with genomic distance between all sites within 50 kb for the Kocoge parasite samples. Chromosomes are shown in column 1, distance in column 2, median values in column 3. mayuge_median.ld.txt.zip - The decay of linkage disequilibrium with genomic distance between all sites within 50 kb for the Mayuge parasite samples. Chromosomes are shown in column 1, distance in column 2, median values in column 3. Misc files: schools.list - List of samples and schools where they were sampled.

Raw anonymised Schistosoma mansoni and human host ABO blood group dataset from school-aged children, aged 6-14, from three high endemicity villages in Mayuge District Uganda. Data include age, sex, height, weight, village, date, timepoint pre and post treatment with praziquantel and albendazole, up to three days per timepoint of duplicate Kato-Katz slides with infection intensity data for S. mansoni and soil-transmitted helminth infections, including hookworm, Ascaris lumbricoides, Trichuris trichiura, Hymenolepis nana, and Enterobius vermicularis, recorded, and ABO and Rhesus positive or negative blood group type.

Time to reinfection and microsatellite data for epidemiological and genetic analyses

Dataset to underpin associated article.