DOI: 10.5281/zenodo.4940589
title: Data release: Whole-genome sequencing of Schistosoma mansoni reveals extensive diversity with limited selection despite repeated mass drug administration

description: Source data used in the publication: Berger et al. (2021) - Provisional title: Whole-genome sequencing of <em>Schistosoma mansoni</em> 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 db22@sanger.ac.uk, manuscript that goes with these data is about to be resubmitted) <strong>File descriptions:</strong><br> prunedData.eigenvec - <strong>Figure 2a&b </strong>- PCA eigenvectors prunedData.eigenval - <strong>Figure 2a&b- </strong>PCA eigenvalues autosomes.mdist - <strong>Figure 2c</strong> - PLINK distance matrix used to build the neighbour joining phylogeny all.pi.pixy.schools.txt - <strong>Figure 2d</strong> - Nucleotide diversity results for each school subpopulation autosomes.dxy.5kb.schools.txt - <strong>Figure 2e </strong>- Autosomal D<sub>XY</sub> results between school subpopulations. autosomes.fst.5kb.schools.txt - <strong>Figure 2e </strong>- Autosomal F<sub>ST</sub> results between school subpopulations. admixture_all.txt - <strong>Figure 2f </strong>- 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 - <strong>Figure 3a, Supplementary figure 10a</strong> - Site frequency spectra (allelic proportions at each frequency bin) for each school. TD.all.txt - <strong>Figure 3b</strong> - Tajima's D values calculated in 5 kb windows for each school subpopulation. ALL.MAYUGE.IHS.ihs.out.100bins.norm.txt.zip - <strong>Figure 4a</strong> - Normalised iHS scores for the Mayuge district parasite populations (Selscan output). ALL.TORORO.IHS.ihs.out.100bins.norm.txt.zip - <strong>Figure 4b - </strong>Normalised iHS scores for the Tororo district parasite populations (Selscan output). ALL.MAYUGEvsTORORO.xpehh.xpehh.out.norm.txt.zip - <strong>Figure 4c</strong> - Normalised XP-EHH scores between Mayuge and Tororo parasite populations. MAYUGE_TORORO_2000.windowed.weir.txt.zip - <strong>Figure 4d - </strong>F<sub>ST</sub> values calculated between Mayuge and Tororo populations in 2kb windows. MAYUGE_PI.windowed.pi.zip - <strong>Figure 4e, Supplementary figure 12a&c - </strong>Nucleotide diversity values calculated in 2 kb windows for Mayuge populations. <br> TORORO_PI.windowed.pi.zip - <strong>Figure 4e, Supplementary figure 12a&c</strong>- Nucleotide diversity values calculated in 2 kb windows for Kocoge populations (Tororo district). all.pi.treat.fix.txt.zip - <strong>Figure 5a </strong>- Nucleotide diversity results for each treatment subpopulation autosomes.dxy.5kb.treatment.txt - <strong>Figure 5b </strong>- Autosomal D<sub>XY</sub> results between clearance phenotype subpopulations. autosomes.fst.5kb.treatment.txt - <strong>Figure 5b </strong>- Autosomal F<sub>ST</sub> results between clearance phenotype subpopulations. fst.windows.2kb.treatment.txt.zip - <strong>Figure 5c </strong>- F<sub>ST</sub> values for comparisons between different treatment groups (Pre-treatment, post-treatment (good clearers), post-treatment (poor clearers)) assoc_err_binary.txt.zip - <strong>Figure 5d </strong>- 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 - <strong>Figure 5e </strong>- 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 - <strong>Supplementary figure 1</strong> - Normalised depth of read coverage (column 4) calculated in 25 kb windows (columns 2&3) across all samples for all chromosomes (column 1).<br> cohort.genotyped.txt.zip - <strong>Supplementary figure 2(a-f) </strong>- Variant quality site values (used to inform variant site retention or removal). <br> hard_filtered.imiss.txt - <strong>Supplementary figure 2g </strong>- Per sample variant missingness (used to inform quality control)<br> hard_filtered_filtindv.lmiss.txt.zip - <strong>Supplementary figure 2h </strong>- Per site missingness (used to inform quality control) cv_scores.txt - <strong>Supplementary figure 5 </strong>-<strong> </strong>ADMIXTURE CV scores (column 2) for each population size (1). *_SMC_SE.csv - <strong>Supplementary figure 6 -</strong>SMC++ results (from 25 subsampled replicates) for each school subpopulation and outgroup samples. smcpp.csv -<strong> Supplementary figure 7 </strong>- SMC++ results for each school subpopulation and outgroup samples. <br> pi.per_host.txt.zip - <strong>Supplementary figure 8a-d</strong> - Nucleotide diversity values for each host infrapopulation. <br> sexing.csv - <strong>Supplementary figure 9</strong>- 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 - <strong>Supplementary figure 11 </strong>- Tajima's D values calculated for the Mayuge population in 2kb windows. Tororo_TAJIMA_D.Tajima.D.2kb.txt.zip - <strong>Supplementary figure 11 </strong>- Tajima's D values calculated for the Tororo population in 2kb windows. genes.bed - <strong>Supplementary figures 13-18</strong> - Coordinates of gene models (<em>S. mansoni </em>v7 annotation). KOCOGE_SITE_PI.sites.pi.txt.zip - <strong>Supplementary figures 13-18 </strong>- Per site nucleotide diversity values MAYUGE_TORORO_sites.weir.fst.txt.zip - <strong>Supplementary figures 13a-18a </strong>- Per site F<sub>ST</sub> values between Mayuge and Tororo populations. kocoge_median.ld.txt.zip - <strong>Supplementary figure 19 </strong>- 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 - <strong>Supplementary figure 19 </strong>- 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).