Simultaneous genotyping of snails and infecting trematode parasites using high-throughput amplicon sequencing.
Data files
Jul 20, 2021 version files 5.28 MB
-
BLAST_expected.csv
10.48 KB
-
BLASTED.csv
10.74 KB
-
COI1_snail_alignment.nex
13.56 KB
-
COI1_trematode_alignment.nex
42.04 KB
-
COI2_snails_alignment.nex
7.32 KB
-
COI2_trematode_alignment.nex
36.74 KB
-
cytb_trematode_alignment.nex
64.96 KB
-
dada2_script_nano.R
5.14 KB
-
ID_sample.csv
44.84 KB
-
infection.csv
5.76 KB
-
ITS1_snail_alignment.nex
27.27 KB
-
ITS2_trematode_alignment.nex
304.29 KB
-
NAD1_snails_alignment.nex
11.04 KB
-
NAD1_trematode_alignment.nex
64.84 KB
-
RD_diagnostic.csv
4.71 KB
-
README.txt
5.57 KB
-
seqtab_nochim_long.csv
4.43 MB
-
unique_ASVs_dada2.fasta
193.72 KB
Abstract
Several methodological issues currently hamper the study of entire trematode communities within populations of their intermediate snail hosts. Here we develop a new workflow using high-throughput amplicon sequencing to simultaneously genotype snail hosts and their infecting trematode parasites. We designed primers to amplify 4 snail and 5 trematode markers in a single multiplex PCR. While also applicable to other genera, we focused on medically and economically important snail genera within the Superorder Hygrophila and targeted a broad taxonomic range of parasites within the Class Trematoda. We tested the workflow using 417 Biomphalaria glabrata specimens experimentally infected with Schistosoma rodhaini, two strains of Schistosoma mansoni, and combinations thereof. We evaluated the reliability of infection diagnostics, the robustness of the workflow, its specificity related to host and parasite identification, and the sensitivity to detect co-infections, immature infections, and changes of parasite biomass during the infection process. Finally, we investigated its applicability in wild-caught snails of other genera naturally infected with diverse trematode assemblages. After stringent quality control the workflow allows the identification of snails to species level, and of trematodes to taxonomic levels ranging from family to strain. It is sensitive to detect immature infections and changes in parasite biomass described in previous experimental studies. Co-infections were successfully identified, opening the possibility to examine parasite-parasite interactions such as interspecific competition. Altogether, these results demonstrate that our workflow provides a powerful tool to analyze the processes shaping trematode communities within natural snail populations.
The dataset includes bash and R script used ot process the output of our high-throughput amplicon sequencing experiments from raw reads to amplicon sequence variants (ASVs). We also included the R script used to further analyze the ASVs profile of each sample, as well as the R environment of that script to spare users the need to rerun it. The input files required to run the scripts are also included, except for the raw reads which were uploaded as fastq files on NCBI SRA: PRJNA743924. Finally, the abbreviations and their meaning are listed in the README.