Correction
1 Sep 2020: Alzaylaee H, Collins RA, Rinaldi G, Shechonge A, Ngatunga B, et al. (2020) Correction: Schistosoma species detection by environmental DNA assays in African freshwaters. PLOS Neglected Tropical Diseases 14(9): e0008721. https://doi.org/10.1371/journal.pntd.0008721 View correction
Figures
Abstract
Background
Schistosomiasis is a neglected tropical parasitic disease associated with severe pathology, mortality and economic loss worldwide. Programs for disease control may benefit from specific and sensitive diagnostic methods to detect Schistosoma trematodes in aquatic environments. Here we report the development of novel environmental DNA (eDNA) qPCR assays for the presence of the human-infecting species Schistosoma mansoni, S. haematobium and S. japonicum.
Methodology/Principal findings
We first tested the specificity of the assays across the three species using genomic DNA preparations which showed successful amplification of target sequences with no cross amplification between the three focal species. In addition, we evaluated the specificity of the assays using synthetic DNA of multiple Schistosoma species, and demonstrated a high overall specificity; however, S. japonicum and S. haematobium assays showed cross-species amplification with very closely-related species. We next tested the effectiveness of the S. mansoni assay using eDNA samples from aquaria containing infected host gastropods, with the target species revealed as present in all infected aquaria. Finally, we evaluated the effectiveness of the S. mansoni and S. haematobium assays using eDNA samples from eight discrete natural freshwater sites in Tanzania, and demonstrated strong correspondence between infection status established using eDNA and conventional assays of parasite prevalence in host snails.
Conclusions/Significance
Collectively, our results suggest that eDNA monitoring is able to detect schistosomes in freshwater bodies, but refinement of the field sampling, storage and assay methods are likely to optimise its performance. We anticipate that environmental DNA-based approaches will help to inform epidemiological studies and contribute to efforts to control and eliminate schistosomiasis in endemic areas.
Author summary
Schistosomiasis, otherwise known as bilharzia or snail fever, is a prevalent human disease found across tropical regions of the world and is a major cause of disability. The disease is acquired from exposure to the schistosome infectious larvae released by infected host snails in freshwaters. Programs to restrict the transmission of schistosomiasis would benefit from rapid and reliable diagnostic methods to detect schistosomes. Here we report a study that has developed new diagnostic tools to identify the DNA from three human-infecting Schistosoma species within water samples. This “environmental DNA” (eDNA) approach requires the filtering and laboratory analyses of water samples, and avoids the requirements to locate, identify and individually test the infectious status of host snails. Our results showed that eDNA methods detect the presence of the parasite in freshwater bodies. However, there is need for further refinement for sampling and laboratory techniques to improve the performance of the assays. We anticipate that eDNA approaches will provide information on the distribution and abundance the water-borne parasites, and potentially contribute to the control and elimination of schistosomiasis.
Citation: Alzaylaee H, Collins RA, Rinaldi G, Shechonge A, Ngatunga B, Morgan ER, et al. (2020) Schistosoma species detection by environmental DNA assays in African freshwaters. PLoS Negl Trop Dis 14(3): e0008129. https://doi.org/10.1371/journal.pntd.0008129
Editor: Robert M. Greenberg, University of Pennsylvania, UNITED STATES
Received: September 11, 2019; Accepted: February 10, 2020; Published: March 23, 2020
Copyright: © 2020 Alzaylaee et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: All relevant data are within the manuscript and its Supporting Information files.
Funding: HA was funded by the Embassy of the Kingdom of Saudi Arabia and Prince Nourah Bin Abdulrahman University. MJG and BPN were supported by the Royal Society-Leverhulme Trust Africa Award AA130107. MJG and RAC were supported by Global Challenges Research Fund pump priming project “Aquatic environmental DNA for enhanced health and food security in Tanzania.” The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Schistosomiasis, a neglected tropical disease caused by trematodes of the genus Schistosoma, is the second most prevalent parasitic disease in humans, after malaria [1, 2]. Studies have indicated that upwards of 207 million people may be infected globally, with a further 779 million individuals at risk of infection [2]. Around 93% of infected people live in sub-Saharan Africa [2, 3], where widespread poverty and the absence of clean water and sanitary facilities contribute to the incidence and propagation of the disease [4, 5]. Schistosomiasis is associated with health problems including malnutrition, anaemia, and impaired childhood development. It is also a major impediment to education, health nutrition and socio-economic development in the infected populations [1, 6, 7, 8].
There are three main Schistosoma species that infect humans. S. haematobium causes urogenital schistosomiasis and is present across Africa and the Middle East. S. mansoni, the agent of intestinal schistosomiasis, is prevalent across Africa, the Middle East and South America. S. japonicum causes hepatic schistosomiasis and is mainly found in the Philippines, Indonesia and China. As trematodes, Schistosoma species have a complex life cycle that requires an intermediate freshwater snail host and a final vertebrate host [1, 8]. Snails that are responsible for schistosomiasis in humans belong to one of three genera: Biomphalaria (for S. mansoni), Bulinus (for S. haematobium) and Oncomelania (for S. japonicum) [8]. Transmission to humans occurs when the free-swimming larval form (cercaria) is shed from infected freshwater snails and penetrates the skin of the definitive host [8, 9]. Once mature, the male and female schistosomes mate, reproduce, and produce eggs that are excreted via urine or faeces, depending on the species. Eggs excreted into freshwater environments hatch and release mobile miracidia which infect snails to continue the cycle [8].
Awareness of safe water sources by human populations in endemic regions is critical for prevention of schistosomiasis [10, 11]. However, assessment and therefore mitigation of the risk posed by bodies of freshwater is not straightforward. Despite large integrated control programs during the last few decades, disease transmission continues due to the absence of rapid and reliable diagnostic tools to detect the schistosomes in endemic areas [12]. Historically the presence of schistosomes has been established by snail-based surveys coupled with microscopic and/or molecular detection of parasitic larvae within the snails [13–18]. These surveys are labour intensive and demand specific training and expertise. Moreover, the sensitivity of this method is generally low given that snail infection rates can be as low as 1–2% of the population, even where there is a high prevalence of infected humans [19, 20].
Over recent years conservation biologists and ecologists have widely adopted the use of environmental DNA (eDNA) methods to establish the presence of species that are difficult to study using other approaches [21–27]. These methods can be particularly useful for studying parasites because they can be challenging to locate and reliably identify. From a parasitological perspective, the term eDNA can be defined as “DNA extracted from an environmental or organismal matrix” [28]. Specifically, this definition includes target DNA from whole microscopic organisms present in the sample, and is not restricted to DNA in solution or cellular debris [28]. Successful detections of trematode parasites using aquatic eDNA include Ribeiroia ondatrae in North America [29], Opisthorchis viverrini in southeast Asia [30], and Fasciola hepatica and Calicophoron daubneyi in Europe [31]. Previous studies have also shown the feasibility of recovering Schistosoma mansoni DNA from both field and laboratory environments [32, 33]. They have also demonstrated that Schistosoma DNA present in laboratory water decayed below levels of detection eight days after the removal of source snails [33]. Collectively, this evidence is highly supportive of the potential of eDNA methods for improved surveillance of all harmful schistosome species [34], but also indicates that improvements are needed to adapt these methods to large-scale surveillance, and to identify species other than S. mansoni.
In this study, we built on previous studies by developing novel eDNA methods to detect human Schistosoma in freshwater bodies. Specifically, we report the development of primers and probes for the mitochondrial 16S rRNA region for S. mansoni, S. haematobium and S. japonicum, and undertake thorough tests of assay specificity. We then validated sensitivity of the S. mansoni assay by testing eDNA from water collected from aquaria holding infected host snails. Finally, we evaluated the applicability and sensitivity of the new assays to detect S. mansoni and S. haematobium eDNA in Tanzanian freshwaters where both species are co-endemic.
Materials and methods
Primer and probe design
Species-specific sets of primers and probes were developed to target short (< 150 basepair) fragments of the mitochondrial DNA 16S rRNA gene in S. japonicum, S. mansoni and S. haematobium. The following GenBank sequences were used as templates: S. mansoni HE601612, S. haematobium DQ157222 and S. japonicum JQ781206. Sequences were aligned to available mitochondrial 16S rRNA sequences of eight other Schistosoma species using ClustalW in Bioedit v.7.2.6 [35], namely Schistosoma curassoni AP017708, Schistosoma mekongi AF217449, Schistosoma margrebowiei AP017709, Schistosoma spindale DQ157223, Schistosoma bovis QMKO01004774, Schistosoma rodhaini LL973454, Schistosoma indicum EF534284 and Schistosoma incognitum EF534285. We used PrimerBlast [36] and PrimerQuest (Integrated DNA Technologies, Coralville IA; https://eu.idtdna.com/Primerquest/Home/Index) to design primers and probes. In both tools we selected default settings, except PCR product size was restricted to between 70 to 170 bp. To avoid cross-species amplification, primers were selected that had at least two mismatches with non-target species. Probes were designed with a dual labelled 5(6)-carboxy-fluorescein (FAM) fluorescent tag at the 5’ end and with a Black Hole Quencher 1 at the 3’ end (Integrated DNA Technologies, Coralville, IA). Final expected lengths of PCR products for S. mansoni, S. haematobium and S. japonicum were 104, 143 and 87 bp, respectively (Table 1).
Specificity of S. mansoni, S. haematobium and S. japonicum assays
Specificity of the primer/probe sets was established using quantitative PCR (qPCR) tests on both genomic and synthetic DNA of known origin. Genomic DNA from S. mansoni, S. haematobium and S. japonicum was sourced from BEI Resources Repository of the National Institute of Allergy and Infectious Diseases (NIAID, Manassas, VA, USA), and diluted to 0.01ng/μl for qPCR tests. The synthetic DNA was generated from GenBank published sequences for S. mansoni, S. haematobium, S. japonicum, S. curassoni, S. margrebowiei, S. spindale, S. bovis, S. rodhaini, S. indicum, S. incognitum and S. mekongi (Table 2). Synthesis was carried out by Eurofins Genomics (Ebersberg, Germany) within a pEX-A128 E. coli vector, and Invitrogen by Thermo-Fisher Scientific (Waltham, MA, USA) within a pMA-T E. coli vector or pMA-RQ (AmpR) E. coli vector. Accession numbers, fragment lengths and total synthesised DNA concentration are shown in S1 Table. To visualise the evolutionary relationships of these mitochondrial sequences, they were aligned using Clustalw in DAMBE7 [37] and a maximum likelihood phylogeny was generated using the HKY+Γ model in Mega-X [38] with branch support estimated from 100 bootstrap replicates (Fig 1).
Numbers above branches represent the branch support (proportion of 100 bootstrap replicates.
The template DNA concentration (copies of synthetic template DNA) for each qPCR reaction was calculated using a copy number and dilution calculator (Thermo-Fisher, Waltham, MA, USA; https://bit.ly/2JPZchd). Each total reaction volume of 5 μl per reaction included: 1 μl DNA extract, 2.5 μl Master Mix (PrimeTime Gene Expression Master Mix; IDT), 1.25 μl molecular grade water and 0.25 μl of primer/probe mix. The primer/probe mix included 4 μl of each primer (100 μM, IDT) and 2 μl of probe (100 μM, IDT) and 40 μl of molecular grade water. Thermocycling conditions were 95°C for 3 min, followed by 45 cycles of 95°C for 0.05 seconds and 60°C for 30 seconds. In addition to test samples, negative no-template control samples were included for quality assurance. The reactions were run on an Eco48 real time qPCR system (PCRmax, Staffordshire, UK) in 48-well plates with ROX normalisation. DNA detection was expressed by cycle threshold (Ct) values.
Sensitivity of S. mansoni and S. haematobium assays
We tested the sensitivity of S. mansoni and S. haematobium assays, but not the S. japonicum assay. To make a quantified stock solution for qPCR, the target 16S rRNA region was first amplified from genomic DNA extracts of S. mansoni and S. haematobium sourced from BEI Resources Repository of the NIAID using the species-specific primers described above, and a conventional end-point PCR. Each PCR comprised 5μl polymerase buffer, 3μl MgCl2 of 25 mM, 0.2μl GoTaq DNA polymerase (5U/ml) (Promega, Madison, WI, USA), 0.5μl dNTPs (10mM each), 12.3 μl molecular grade water, 1μl of each 10 μM primer (IDT) and 2μl template DNA (total volume of 25μl per reaction). The PCR conditions were: 95°C for 2 min of initial denaturation, followed by 40 cycles of 95°C for 30 sec, 55°C for 30 sec, and 72°C for 30 sec, and a final extension at 72°C for 5 min. A negative no-template control was included. The presence of a single PCR product of expected size was confirmed on an 1.5% agarose gel stained with gel red nucleic acid gel stain, 10,000X (Biotium, Fremont, CA) and a PCRSizer 100bp DNA Ladder (Norgen Biotek, Thorold, Canada).
To make the serial dilutions to measure assay sensitivity, each PCR product was then purified using the QIAquick PCR purification kit (Qiagen) following the manufacturer’s protocol. Then, the concentrations of DNA were measured using a Qubit fluorometer (Invitrogen, Waltham, MA). The number of target copies was calculated using the Thermo-Fisher DNA copy number and dilution calculator (above). Then, the stock solution was stored at -20°C for qPCR tests. Before each qPCR test, a dilution series of known concentrations ranging from 1,000,000 copies/μl decreasing ten-fold down to 1 copy/μl was prepared for each assay (Table 3). These serial dilutions were used within two days to limit potential effects of DNA degradation.
Three qPCR technical replicates were run for each ten-fold standard dilution, including negative no-template controls, using the qPCR protocol described above. The limit of detection within a single qPCR reaction (LODI) was defined as the lowest concentration where there is a 95% chance of amplification success in any one individual PCR reaction. We also calculated the lowest concentration with a 95% chance of amplification success in any one of three technical replicate qPCR reactions of the same sample (LODIII), which we considered a useful metric given that qPCR reactions of eDNA are typically conducted in triplicate. These concentration limits were derived from all standards run and calculated by fitting logistic models [39] using CurveExpert Basic 2.1.0 (Hyams Development). The limit of quantification (LOQ) was defined as the lowest concentration at which 90% of all standards run were able to be amplified, following protocols described in [40].
Environmental DNA from aquarium water samples
To validate the accuracy of the S. mansoni qPCR assay, water samples were collected from four aquaria housing Biomphalaria glabrata host snails infected with S. mansoni, at the Wellcome Sanger Institute (WSI), Cambridgeshire, UK. The complete life cycle of Schistosoma mansoni (NMRI strain) is maintained at the WSI by breeding and infecting B. glabrata snails and mice. Mouse infections are performed under the Home Office Project Licence No. P77E8A062 held by GR. We did not specifically quantify the number or biomass of infected snails in aquaria, but the total number of snails in each aquarium was recorded. Additionally, water samples were collected from two aquaria holding non-infected B. glabrata, and a sample of sterile water was collected at the site for use as a negative control. Sampled water was filtered through a Sterivex filter with a pore size 0.22 μm and a polyethersulfone membrane (EMD Millipore corporation, UK), using a peristaltic pump unit. Volumes of water filtered in a single replicate varied between 60–500 ml (Table 4). After filtration, absolute ethanol was pushed through the filter with a sterile 50 ml syringe to preserve the samples, and each filter was kept individually in a labelled 118 ml capacity Whirl-Pak bag (Sigma-Aldrich, UK). Filters were placed in a polystyrene box or a cooler with ice during transportation to the laboratory, and kept at -20°C until being processed. Extraction of eDNA from the individual filter was performed using the DNeasy Power Water Kit (Qiagen, UK) following the manufacturer’s protocol, but eluting into 50 μl of EB buffer. The S. mansoni DNA was quantified using qPCR procedures as described above for the specificity assays. Three PCR replicates were carried out for each sample, alongside three replicates of a negative no-template PCR control, and three replicates of each concentration in a standard curve serial dilution of control positive PCR-derived DNA of S. mansoni (ranging from 1,000,000 copies/μl to 1 copy/μl).
Environmental DNA from the natural environment
To test the S. mansoni and S. haematobium assays on samples from natural habitats, water samples were collected from the surface of freshwaters at eight locations in Tanzania during September 2018 (Table 5; Fig 2). These included locations where both S. mansoni and S. haematobium were plausibly co-endemic due to the presence of host snails and close proximity of human settlements. It also included locations where both S. mansoni and S. haematobium were plausibly absent due to the absence of host snails and the greater distance from settlements. At each location coordinates were obtained using a handheld etrex GPS (Garmin, Olathe, KS), and the pH, conductivity, total dissolved solids and temperature of the water were measured using a multimeter (Hach, Loveland, CO) (Table 5). At each location MJG conducted a 10 minute search for aquatic snails in shallow water benthic habitats including rock surfaces, soft sediments and marginal vegetation. The collected snails from each site were preserved in absolute ethanol and stored in 118 ml capacity whirlpak bags. Additional samples were also collected at sites A and B allowing more thorough testing of the frequency of infected snails at these sites (Table 5).
In total eight sites were surveyed, labelled A-H. The key results from eDNA survey (helix symbol), host snail survey (snail symbol) and test for the presence of schistosomes in snail tissue (cercariae symbol) are given for S. haematobium (H) and S. mansoni (M) assays, with + indicating the site was positive, while a dot indicates no detection was made. Map drawn using the following open source software and data: DIVA-GIS7.5 (https://www.diva-gis.org); Africa river network (15 sec resolution), African drainage basin (15sec resolution) and African Digital Elevation Model data (30 sec resolution) from HydroSHEDS (https://www.hydrosheds.org); African Water Bodies shapefile from RCMRD GeoPortal (http://geoportal.rcmrd.org).
At each of the eight sites, three replicate eDNA samples were collected. Each collection involved the pumping of water though a sterile 0.22 μm Sterivex filter using a sterile 50 ml single-use syringe. Volumes of water filtered varied between 100 to 400 ml (Table 6). After filtration, absolute ethanol was pushed through the filter with a sterile 50 ml syringe to preserve the samples, and filters were kept in labelled 118 ml capacity whirlpak bags. Bottled drinking water was filtered alongside the samples in the field and used as a negative sampling control. DNA was isolated from the filters using the DNeasy Power Water Kit as described above. Each qPCR was performed in triplicate on each eDNA extract, alongside the negative sampling control, the negative non-template PCR control and each of the seven ten-fold serial dilutions of control positive DNA. Analyses of samples for S. mansoni and S. haematobium were conducted separately.
To quantify the infection status of individual molluscs, a small sample of tissue (no more 20mg) was dissected, and DNA extracted using the DNeasy Blood & Tissue Kit (Qiagen, UK) according to the manufacturer’s protocol. Analyses were conducted using the qPCR approach as described above, although only presence or absence of an amplification is reported herein. We chose a PCR-based screening approach as it allowed us the opportunity to collect the samples in the field and screen the samples at a later date accurately for the focal species; however, the method can amplify individuals with prepatent infections that are not currently shedding cercariae [41]. To confirm the identity of PCR products amplified during these tissue assays, the PCR products from three Biomphalaria pfeifferi samples and three Bulinus globosus samples were Sanger sequenced. PCR products were purified using DNA Clean & Concentrator-5 (Zymo Research, Irvine, CA, USA) according to the manufacturer’s protocol, and sequenced by Eurofins Genomics (Ebersberg, Germany) using the forward PCR primer. The identity of derived sequences was confirmed using a default BLAST search against the NCBI GenBank nucleotide database.
Results
Environmental DNA assays exhibit high levels of species specificity
In qPCR experiments testing the specificity of assays against template genomic DNA, no amplifications from negative controls were observed, and we observed no cross-species amplifications between S. mansoni, S. haematobium, S. japonicum.
When synthetic DNA that included species within the same species group (Fig 1) as the target species was employed, the results indicated some cross-species amplification with the S. haematobium and S. japonicum assays albeit with lower levels of qPCR amplification intensity for the non-target species (Table 2). Specifically, the assay for S. japonicum not only amplified synthetic S. japonicum DNA (average Ct score 20.19), but also amplified synthetic DNA of closely-related S. mekongi (average Ct score 28.47). Likewise, the assay for S. haematobium amplified synthetic S. haematobium DNA (average Ct score 27.93), but also amplified synthetic DNA of the closely-related S. bovis (average Ct score 35.67) and S. curassoni (average Ct score 30.59) (Table 2).
High sensitivity of S. mansoni and S. haematobium eDNA assays
Efficiencies of the S. mansoni qPCR assay across 17 dilution series ranged from 91.55 to 110.86%, while R2 ranged from 0.97 to 0.99. The S. mansoni assay amplified all standards between 1,000,000 copies/μl and 100 copies/μl. Amplification success became inconsistent at 10 copies/μl (32/51 PCR reactions) and was rare at the lowest concentration of 1 copy/ μl (2/51 amplifications). The resolved LOQ for the S. mansoni assay was therefore 100 copies/μl. The resolved limit of detection where there was 95% probability of successful individual PCR amplification (LODI) was 41.68 copies/μl, while the 95% probability of successful amplification in at least one of the triplicate qPCRs (LODIII) was 10.84 copies/μl (Table 3; Fig 3).
a) S. mansoni assay–probability of amplification in any one qPCR across all standard templates tested. b) S. mansoni assay–probability of amplification in any one standard template that is subject to triplicate qPCR, c) S. haematobium assay–probability of amplification in any one qPCR across all standard templates tested, d) S. haematobium assay–probability of amplification in any one standard template that is subject to triplicate qPCR. Lines represent logistic models of the form y = a / (1 + b*e(-cx)). a) a = 1.002, b = 27.303, c = 3.826. b) a = 1.000, b = 7.501, c = 4.786. c) a = 1.004, b = 1.009, c = 1.415. b) a = 1.004, b = 0.342, c = 1.399. All models r > 0.995. Grey shading indicates 95% confidence intervals.
Efficiencies of the S. haematobium qPCR assay across 12 dilution series ranged from 94.87 to 105.45%, while R2 ranged from 0.98 to 0.99. The S. haematobium assay amplified all standards between 1,000,000 copies/μl and 1,000 copies/μl. Amplification was mostly successful at 100 copies/μl (34/36 PCR reactions) and 10 copies/μl (29/36 PCR reactions), but was inconsistent at 1 copy/μl (Table 3). The resolved LOQ for the S. haematobium assay was therefore 100 copies/μl. The resolved limit of detection where there was 95% probability of successful individual PCR amplification (LODI) was 108.39 copies/μl, while the 95% probability of successful amplification in at least one of the triplicate qPCRs (LODIII) was 19.19 copies/μl (Table 3; Fig 3).
Environmental DNA assays detected schistosomes from aquarium water samples
The qPCR assay detected S. mansoni DNA in all infected aquaria with mean Ct values ranging from 24.64 to 27.77 (Table 4). There were no positive amplifications from non-infected aquaria, or sterilised water controls. Average qPCR efficiency was 103% across the two dilution series assays (range 99.38–106.95) with mean R2 value of 0.99. All positively amplified samples were above LOQ, LODI and LODIII values for the species assay, and the mean number of copies of DNA detected in each sample ranged from 7,565 to 48,264 copies/μl (Table 4)
Environmental DNA detected the presence of schistosomes from the natural environment
Four out of the eight screened sites were positive for S. mansoni eDNA (sites we named A, B, C and D; Tables 5 & 6, Fig 2). At two of these eDNA positive locations, infected B. globosus gastropods were present (sites A and B), while at two of the locations B. globosus was not found (sites C and D; Table 5). Of the four locations that were eDNA negative, only one had B. globosus present, and these snails were determined to be infected through PCR analysis of tissue DNA (Table 5, Fig 2). Thus, the eDNA assay was congruent with the PCR tests for infected snails in 5/8 (62.5%) of locations. Notably, the estimated numbers of S. mansoni copies resolved in most samples were below the defined LOQ of 100 copies/μl, the LODI of 41.68 copies/μl and the LODIII of 10.84 copies/μl.
For S. haematobium, the target species was detected in eDNA at a single site (Site B, Table 6). At this site all replicates successfully amplified, and infected Bulinus snails were abundant (Table 5). Infected Bulinus were only detected at one other location (Site A), and this site was resolved as negative for eDNA. Thus, the eDNA assay was congruent with the PCR tests for infected snails in 7/8 (87.5%) of locations. The estimated numbers of S. haematobium copies were below the defined LOQ of 100 copies/μl, below the LODI of 108.39 copies/μl, and mainly below the LODIII of 19.19 copies/μl.
Sequencing confirmed the source of PCR product from B. globosus tissue as S. haematobium, (100% BLAST match to GenBank accession EU567132). Similarly, the source of PCR product from B. pfeifferi tissue was confirmed as S. mansoni (100% BLAST match to GenBank accessions AF130787 and LR214937).
Discussion
In this study PCR-based assays have been developed to detect the three main human-infecting Schistosoma species from tropical freshwaters; S. mansoni, S. haematobium and S. japonicum. Probes and primer pairs were designed to amplify the 16S rRNA region of the mitochondrial DNA, with product size <150 bp which makes them well suited to amplify small fragments of degraded DNA. Our study provided newly characterised primers that reliably amplify and distinguish each of the three main species that infect humans, unlike previous investigations of Schistosoma eDNA [32, 33]. In addition, we have exhaustively tested the sensitivity and species-specificity of the assays in-vitro. The latter is critical for the evaluation of Schistosoma eDNA assays, not only because multiple species within the Schistosomatidae family may be present within the same area, potentially leading to false positive results, but also because in-silico testing of primers has frequently been shown to yield inaccurate results when compared to those from in-vitro tests [42].
In addition to testing genomic DNA from the focal species, we chose to use synthetic DNA in the in-vitro species specificity tests. This approach was selected in order to ensure confidence in the species that the DNA represented, without risk of using natural samples potentially contaminated during their collection and/or curation. Additionally, we were able to test our assays against DNA from a broad variety of relevant species, natural samples of which would be difficult to obtain. Our assays consistently amplified both genomic and synthetic DNA of the target species. However, they also amplified closely-related species. Specifically, the S. haematobium assay amplified S. bovis and S. curassoni from the ‘haematobium group’, while the S. japonicum assay amplified S. mekongi from the ‘japonicum group’ (Fig 1). Such cross-species amplification could lead to false-positive amplifications, and in practice the impact of any cross-species amplification will depend on the geographic context where the assay is used. For example, application of the S. haematobium assay at locations in north-west Africa where it occurs with S. curassoni and S. bovis may be affected by false positives. Therefore, during deployment of the assays it may be beneficial to undertake screening to establish Schistosoma species present using alternative morphological or genetic methods. Another consideration is that there is potential for cross-species amplification where introgressive hybridization is present among closely-related species, for example between S. haematobium and S. bovis, between S. haematobium and S. curassoni, and between S. haematobium and S. mattheei [43–45]. In such cases, hybridization can lead to sharing of mitochondrial DNA between parental species, and an environmental DNA assay alone would not be able to distinguish between the species, or identify the presence of hybrids.
We showed it was possible to reliably amplify schistosome eDNA from filters collecting material from natural water bodies and experimental aquaria. None of our negative control samples amplified, and thus we have confidence that the PCRs were amplifying only target schistosome material. It is possible that the filters contained environmental DNA in solution, in cellular debris, and also as whole individual cercariae. From a practical perspective for testing for the presence or absence of the schistosome parasites, it may be of little consequence whether the sample is obtained from partial or whole organisms. However, if the method is to be used for evaluating the relative biomass of schistosomes in the environment, then sampling of whole organisms may be problematic. It is possible that the use of pre-filtration steps may be required to exclude whole organisms, potentially providing more standardised quantitative spatial and temporal comparisons [33]. Alternatively, it would be possible to modify the environmental DNA methods used here to focus sampling exclusively cercariae [46,47]. Such “qPCR cercariometry” requires relatively large volumes of water (e.g. 25 litres) that are filtered through a 20μm mesh zooplankton net, before residual material is collected and concentrated allowing qPCR analyses [46,47]. This cercariometry method would allow the collection of large volumes of genetic material from the environment; however such zooplankton nets are expensive and their re-use requires decontamination between sampling events if they are to remain effective.
Our eDNA assays showed agreement with tests of direct qPCR-based tests of the infection status of snail hosts. Specifically, eDNA showed 62.5% consistency with tests for S. mansoni in host gastropods, and 87.5% consistency with tests for S. haematobium in host gastropods. Thus, our assays are similar in performance to those reported by Sengupta et al. [33] for S. mansoni in Kenya with a 71% agreement between eDNA and conventional survey methods. Importantly, both our study and Sengupta et al. [33] detected cases where the conventional survey failed to detect evidence of S. mansoni that was present in the eDNA assay. Such results could have been because schistosomes were present locally, but host snails were rare and therefore difficult to sample. Alternatively, given that the sites at which we detected eDNA but no host snails had flowing water, our analyses could have detected allochthonous schistosome DNA from upstream locations. In either case, the presence of the schistosome eDNA would be indicative of a risk of infection, and perhaps eDNA surveillance could help to identify locations where risk is present but may otherwise go unnoticed. It is also possible that we sampled material from non-transmissive life stages (eggs, miracidia) derived from local contamination sources. Notably, there were also instances where eDNA assays tested negative, but host infected snails were present. In such circumstances it is possible that cercarial production is limited at those locations, or more likely that either eDNA concentrations were below detection limits, or PCR inhibitors were present. Low eDNA concentration is the most likely explanation, however, as the PowerWater extraction kit used has dedicated inhibitor removal steps, and is designed for these kinds of heterogeneous samples. Further simultaneous sampling of both eDNA and host snails is recommended to more reliably estimate the sampling effort required to have high confidence in negative results from these assays.
Multiple studies have described molecular-based diagnostic techniques for Schistosoma species [48–51]. These methods, alongside technical improvements in the reliability of PCR, have resulted in an increasing number of research centres in low income countries having access to the real-time PCR diagnostic technology [52]. Widespread uptake of eDNA methods for Schistosoma surveillance is therefore feasible in principle. However, the extent to which eDNA methods can replace existing survey protocols will depend on the reliability of the assays for detecting low concentrations of eDNA, as this study identified sites that were eDNA negative but positive by alternative sampling methods. This is most likely because eDNA levels were below detection limits. Specifically, our tests demonstrated that extracted concentrations of schistosome eDNA were at most ~40 copies/ul (which corresponds with ~16 copies/ml of sampled water), which is below our measured limits of detection where individual PCRs becomes 95% reliable. Expectedly, the S. mansoni assay was most successful at locations A and B where the Biomphalaria colonies were large and infection prevalence was extremely high (40% and 88% of snails, respectively). Moreover, the S. haematobium assay was most successful at the location B where the Bulinus colony was large and infection prevalence was as high as 94%. At location A, where the Bulinus colony was small yet infection prevalence was 60%, the assay was negative. Thus, the success of the assay is likely to depend both on snail abundance and infection prevalence. Importantly, all the locations where we identified schistosome-positive snails had much higher rates of infection than is typically observed, which can be around 2% [19], although we recognize that the PCR-based tests of snail infection may have revealed cases of prepatent infections in non-shedding snails [41].
Further work is needed to develop improved methods of collecting, preserving and extracting aqueous environmental DNA in tropical environments for schistosome surveillance, including the number and volume of field replicates, the DNA preservation reagents used, and the extraction protocols. Recently Sengupta et al. [33] estimated that up to seven one litre water samples may be enough to provide 95% confidence in the presence or absence of S. mansoni. In this study we used smaller volumes, with the upper limits determined by filter clogging. It may be that more filters are needed to collect the higher volumes of material required for reliable eDNA assays in turbid environments, which are often characteristic of schistosome transmission sites. The characteristics of the environment may also influence the persistence of eDNA and chances of detection [32,53]. For example, suspended sediment can adsorb DNA [54], and the presence of PCR inhibitors such as humic acids, algae and siliceous sediments can determine the success of qPCR-based eDNA assays [55,56]. Therefore further sampling and tests will be required to optimize the eDNA detection power in light of environmental variation. In particular, tests for the presence of PCR inhibition within sampled sites would also be valuable.
To conclude, rapid and inexpensive diagnosis and surveillance of schistosome prevalence in freshwaters will provide timely advice to stakeholders who could promote interventions to interrupt transmission of schistosomiasis in developing countries. As social and environmental change proceeds ever more rapidly, strategies for control of schistosomiasis must increasingly be adapted to local context, and tools for rapid evaluation of parasite presence at fine spatial scales are therefore urgently needed [57]. Our results, together with those of recent parallel studies [32, 33] provide promising indications that eDNA approaches could be considered as a tool for standard monitoring of parasite presence. However, further research is needed to resolve a sampling design that is reliable across habitat types and robust to idiosyncrasies of tropical fieldwork. Moreover, further work is needed to address the shortcomings of the methods we used in terms of specificity and sensitivity. If such limitations could be overcome, the tool could be invaluable for surveillance studies and provide direct support to control strategies for this neglected tropical disease. However, it has been highlighted that the utility of eDNA assays is dependent on the nature of the schistosome life stages present in the environment [58]. More specifically, eDNA assays that yield positive results are unable to distinguish between a site where transmission is active (production of cercariae) or simply contaminated by eggs and miracidia from urine and faeces. Hence, eDNA assays alone may be unable to reliably distinguish locations where transmission has been successfully interrupted by intervention, and in such cases it will need to be coupled with DNA screening of collected snails [58].
Supporting information
S1 Table. Sequences used to develop synthetic DNA for testing species specificity of probes.
https://doi.org/10.1371/journal.pntd.0008129.s001
(DOCX)
Acknowledgments
We thank staff from the Tanzania Fisheries Research Institute for contributions to fieldwork. The Tanzania Commission for Science and Technology (COSTECH) provided fieldwork permits. We thank Fiona Allan for providing museum samples we used in preliminary tests of the assays. We thank BEI Resources Repository of the National Institute of Allergy and Infectious Diseases for the provision of schistosome material. We thank Richard Wall for his supervisory support, and Michelle Thieme for permission to use HydroSHEDS data for mapping. We thank Marcello Otake Sato and Russell Stothard for valuable comments on the manuscript.
References
- 1. Gryseels B, Polman K, Clerinx J, Kestens L. Human schistosomiasis. Lancet. 2006;368(9541):1106–1118. pmid:16997665
- 2. Steinmann P, Keiser J, Bos R, Tanner M, Utzinger J. Schistosomiasis and water resources development: systematic review, meta-analysis, and estimates of people at risk. Lancet Infect Dis. 2006;6(7):411–425. pmid:16790382
- 3. Hotez PJ, Kamath A. Neglected tropical diseases in sub-Saharan Africa: review of their prevalence, distribution, and disease burden. PLoS Negl Trop Dis. 2009;3(8):e412. pmid:19707588
- 4. Chitsulo L, Engels D, Montresor A, Savioli L. The global status of schistosomiasis and its control. Acta Trop. 2000;77(1):41–51. pmid:10996119
- 5. Adenowo AF, Oyinloye BE, Ogunyinka BI, Kappo AP. Impact of human schistosomiasis in sub-Saharan Africa. Brazi J Infect Dis. 2015;19(2):196–205. pmid:25636189
- 6. King CH, Dickman K, Tisch DJ. Reassessment of the cost of chronic helmintic infection: a meta-analysis of disability-related outcomes in endemic schistosomiasis. Lancet. 2005;365(9470):1561–1569. pmid:15866310
- 7. King CH, Dangerfield-Cha M. The unacknowledged impact of chronic schistosomiasis. Chronic Illn. 2008;4(1):65–79. pmid:18322031
- 8. Colley DG, Bustinduy AL, Secor WE, King CH. Human schistosomiasis. Lancet. 2014;383(9936): 2253–2264. pmid:24698483
- 9. Machacek T, Turjanicova L, Bulantova J, Hrdy J, Horak P, Mikes L. Cercarial dermatitis: a systematic follow-up study of human cases with implications for diagnostics. Parasitol Res. 2018;117(12):3881–3895. pmid:30302587
- 10. Fenwick A, Jourdan P. Schistosomiasis elimination by 2020 or 2030?. Int J Parasitol. 2016;46(7):385–388. pmid:26907938
- 11. Rollinson D, Knopp S, Levitz S, Stothard JR, Tchuente LA, Garba A, et al. Time to set the agenda for schistosomiasis elimination. Acta Trop. 2013;128(2):423–440. pmid:22580511
- 12. Weerakoon KG, Gobert GN, Cai P, McManus DP. Advances in the diagnosis of human schistosomiasis. Clin Microbiol Rev. 2015;28(4):939–967. pmid:26224883
- 13. Wolmarans CT, De Kock KN, Strauss HD, Bornman M. Daily emergence of Schistosoma mansoni and S. haematobium cercariae from naturally infected snails under field conditions. J Helminthol. 2002;76(3):273–277. pmid:12363382
- 14. Kariuki HC, Clennon JA, Brady MS, Kitron U, Sturrock RF, Ouma JH, et al. Distribution patterns and cercarial shedding of Bulinus nasutus and other snails in the Msambweni area, Coast Province, Kenya. Am J Trop Med Hyg. 2004;70(4):449–456. pmid:15100463
- 15. Melo FL, Gomes AL, Barbosa CS, Werkhauser RP, Abath FG. Development of molecular approaches for the identification of transmission sites of schistosomiasis. Trans R Soc Trop Med Hyg. 2006;100(11):1049–1055. pmid:16620892
- 16. Kane RA, Stothard JR, Rollinson D, Leclipteux T, Evraerts J, Standley CJ, et al. Detection and quantification of schistosome DNA in freshwater snails using either fluorescent probes in real-time PCR or oligochromatographic dipstick assays targeting the ribosomal intergenic spacer. Acta Trop. 2013;128(2):241–249. pmid:22100540
- 17. Farghaly A, Saleh AA, Mahdy S, El-Khalik DA, El-Aal NF, Abdel-Rahman SA, et al. Molecular approach for detecting early prepatent Schistosoma mansoni infection in Biomphalaria alexandrina snail host. J Parasit Dis. 2016;40(3):805–812. pmid:27605788
- 18. Bakuza JS, Gillespie R, Nkwengulila G, Adam A, Kilbride E, Mable BK. Assessing S. mansoni prevalence in Biomphalaria snails in the Gombe ecosystem of western Tanzania: the importance of DNA sequence data for clarifying species identification. Parasite Vector. 2017; 10(1): 584. pmid:29169384
- 19. Angelo T, Shahada F, Kassuku A, Mazigo H, Kariuki C, Gouvras A, et al. Population abundance and disease transmission potential of snail intermediate hosts of human schistosomiasis in fishing communities of Mwanza region, north-western, Tanzania. Int J Sci Res. 2014;3(8): 1230–1236.
- 20. Weerakoon K, Gordon C, McManus D. DNA diagnostics for schistosomiasis control. Trop Med Infect Dis. 2018;3(3):81. pmid:30274477
- 21. Brinkman NE, Haugland RA, Wymer LJ, Byappanahalli M, Whitman RL, Vesper SJ. Evaluation of a rapid, quantitative real-time PCR method for enumeration of pathogenic Candida cells in water. Appl Environ Microbiol. 2003;69(3):1775–1782. pmid:12620869
- 22. Walker SF, Salas MB, Jenkins D, Garner TW, Cunningham AA, Hyatt AD, et al. Environmental detection of Batrachochytrium dendrobatidis in a temperate climate. Dis Aquat Org. 2007;77(2):105–112. pmid:17972751
- 23. Ficetola GF, Miaud C, Pompanon F, Taberlet P. Species detection using environmental DNA from water samples. Biol Lett. 2008;4(4):423–425. pmid:18400683
- 24. Jerde CL, Mahon AR, Chadderton WL, Lodge DM. “Sight‐unseen” detection of rare aquatic species using environmental DNA. Conserv Lett. 2011;4(2):150–157.
- 25. Goldberg CS, Pilliod DS, Arkle RS, Waits LP. Molecular detection of vertebrates in stream water: a demonstration using Rocky Mountain tailed frogs and Idaho giant salamanders. PLoS ONE. 2011;6(7): e22746. pmid:21818382
- 26. Thomsen PF, Kielgast JO, Iversen LL, Wiuf C, Rasmussen M, Gilbert MT, et al. Monitoring endangered freshwater biodiversity using environmental DNA. Mol Ecol. 2012; 21(11):2565–2573. pmid:22151771
- 27. Schmidt BR, Kery M, Ursenbacher S, Hyman OJ, Collins JP. Site occupancy models in the analysis of environmental DNA presence/absence surveys: a case study of an emerging amphibian pathogen. Methods Ecol Evol. 2013;4(7):646–653.
- 28. Bass D, Stentiford GD, Littlewood DTJ, Hartikainen H. Diverse applications of environmental DNA methods in parasitology. Trends Parasitol. 2015;31(10):499–513. pmid:26433253
- 29. Huver JR, Koprivnikar J, Johnson PT, Whyard S. Development and application of an eDNA method to detect and quantify a pathogenic parasite in aquatic ecosystems. Ecol Appl. 2015;25(4):991–1002. pmid:26380540
- 30. Hashizume H, Sato M, Sato MO, Ikeda S, Yoonuan T, Sanguankiat S, et al. Application of environmental DNA analysis for the detection of Opisthorchis viverrini DNA in water samples. Acta Trop. 2017;169:1–7. pmid:28108370
- 31. Jones RA, Brophy PM, Davis CN, Davies TE, Emberson H, Stevens PR, et al. Detection of Galba truncatula, Fasciola hepatica and Calicophoron daubneyi environmental DNA within water sources on pasture land, a future tool for fluke control?. Parasite Vector. 2018;11(1):342. pmid:29884202
- 32. Sato MO, Rafalimanantsoa A, Ramarokoto C, Rahetilahy AM, Ravoniarimbinina P, Kawai S, et al. Usefulness of environmental DNA for detecting Schistosoma mansoni occurrence sites in Madagascar. Int J Infect Dis. 2018;76:130–136. pmid:30201503
- 33. Sengupta ME, Hellstrom M, Kariuki HC, Olsen A, Thomsen PF, Mejer H, et al. Environmental DNA for improved detection and environmental surveillance of schistosomiasis. PNAS. 2019;116(18):8931–8940. pmid:30975758
- 34. McManus DP, Gordon C, Weerakoon KG. Testing of water samples for environmental DNA as a surveillance tool to assess the risk of schistosome infection in a locality. Int J Infect Dis. 2018;76:128–129. pmid:30266632
- 35. Hall TA. BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp Ser. 1999; 41: 95–98.
- 36. Ye J, Coulouris G, Zaretskaya I, Cutcutache I, Rozen S, Madden TL. Primer-BLAST: a tool to design target-specific primers for polymerase chain reaction. BMC Bioinformatics. 2012;13(1):134. pmid:22708584
- 37. Xia X. DAMBE7: New and improved tools for data analysis in molecular biology and evolution. Mol Biol Evol. 2018;35(6):1550–1552. pmid:29669107
- 38. Kumar S, Stecher G, Li M, Knyaz C, Tamura K. MEGA X: molecular evolutionary genetics analysis across computing platforms. Mol Biol Evol. 2018;35(6):1547–1549. pmid:29722887
- 39. Forootan A, Sjoback R, Bjorkman J, Sjogreen B, Linz L, Kubista M. Methods to determine limit of detection and limit of quantification in quantitative real-time PCR (qPCR). Biomol Detect Quantif. 2017;12:1–6. pmid:28702366
- 40. Klymus KE, Richter CA, Chapman DC, Paukert C. Quantification of eDNA shedding rates from invasive bighead carp Hypophthalmichthys nobilis and silver carp Hypophthalmichthys molitrix. Biol Cons. 2015;183:77–84.
- 41. Farghaly A, Saleh AA, Mahdy S, El-Khalik DA, El-Aal NFA, Abdel-Rahman SA, Salama MA. Molecular approach for detecting early prepatent Schistosoma mansoni infection in Biomphalaria alexandrina snail host. J Parasit Dis. 2016; 40(3): 805–812. pmid:27605788
- 42. Henriques A, Cereija T, Machado A, Cerca N. In silico vs in vitro analysis of primer specificity for the detection of Gardnerella vaginalis, Atopobium vaginae and Lactobacillus spp. BMC Res Notes. 2012;5(1):637. pmid:23153093
- 43. Huyse T, Webster BL, Geldof S, Stothard JR, Diaw OT, Polman K, et al. Bidirectional introgressive hybridization between a cattle and human schistosome species. PLoS Pathog. 2009;5(9):e1000571. pmid:19730700
- 44. Webster BL, Diaw OT, Seye MM, Webster JP, Rollinson D. Introgressive hybridization of Schistosoma haematobium group species in Senegal: species barrier break down between ruminant and human schistosomes. PLoS Negl Trop Dis. 2013;7(4):e2110. pmid:23593513
- 45. Webster BL, Alharbi MH, Kayuni S, Makaula P, Halstead F, Christiansen R, Juziwelo L, Stanton MC, LaCourse EJ, Rollinson D, Kalua K. Schistosome interactions within the Schistosoma haematobium group, Malawi. Emerg Infect Dis. 2019;25(6):1245–1247. pmid:31107237
- 46. Rudko SP, Reimink RL, Froelich K, Gordy MA, Blankespoor CL, Hanington PC. Use of qPCR-based cercariometry to assess swimmer’s itch in recreational lakes. EcoHealth 2018;15(4):827–839. pmid:30120669
- 47. Rudko SP, Turnbull A, Reimink RL, Froelich K, Hanington PC. Species-specific qPCR assays allow for high-resolution population assessment of four species avian schistosome that cause swimmer's itch in recreational lakes. Int J Parasitol Parasites Wildl. 2019;9: 122–129. pmid:31061794
- 48. Kato-Hayashi N, Kirinoki M, Iwamura Y, Kanazawa T, Kitikoon V, Matsuda H, et al. Identification and differentiation of human schistosomes by polymerase chain reaction. Exp Parasitol. 2010;124(3):325–329. pmid:19948171
- 49. Webster BL, Rollinson D, Stothard JR, Huyse T. Rapid diagnostic multiplex PCR (RD-PCR) to discriminate Schistosoma haematobium and S. bovis. J Helminthol. 2010;84(1):107–114. pmid:19646307
- 50. Webster BL, Emery AM, Webster JP, Gouvras A, Garba A, Diaw O, et al. Genetic diversity within Schistosoma haematobium: DNA barcoding reveals two distinct groups. PLoS Negl Trop Dis. 2012;6(10):e1882. pmid:23145200
- 51. Fernandez-Soto P, Arahuetes JG, Hernandez AS, Aban JL, Santiago BV, Muro A. A loop-mediated isothermal amplification (LAMP) assay for early detection of Schistosoma mansoni in stool samples: a diagnostic approach in a murine model. PLoS Negl Trop Dis. 2014;8(9):e3126. pmid:25187956
- 52. Ten-Hove RJ, Verweij JJ, Vereecken K, Polman K, Dieye L, van-Lieshout L. Multiplex real-time PCR for the detection and quantification of Schistosoma mansoni and S. haematobium infection in stool samples collected in northern Senegal. Trans R Soc Trop Med Hyg. 2008;102(2):179–185. pmid:18177680
- 53. Barnes MA, Turner CR, Jerde CL, Renshaw MA, Chadderton WL, Lodge DM Environmental conditions influence eDNA persistence in aquatic systems. Environ Sci Technol. 2014;48(3):1819–1827. pmid:24422450
- 54. Lorenz MG, Wackernagel W. Adsorption of DNA to sand and variable degradation rates of adsorbed DNA. Appl Environ Microbiol. 1987;53(12):2948–52. pmid:3435148
- 55. McKee AM, Spear SF, Pierson TW. The effect of dilution and the use of a post-extraction nucleic acid purification column on the accuracy, precision, and inhibition of environmental DNA samples. Biol Conserv. 2015;183:70–76.
- 56. Stoeckle BC, Beggel S, Cerwenka AF, Motivans E, Kuehn R, Geist J. A systematic approach to evaluate the influence of environmental conditions on eDNA detection success in aquatic ecosystems. PloS ONE 2017;12(12):e0189119. pmid:29220394
- 57. Booth M, Clements A. Neglected tropical disease control: the case for adaptive, location-specific solutions. Trends Parasitol. 2018;34(4):272–282. pmid:29500033
- 58. Stothard JR, Campbell SJ, Osei-Atweneboana MY, Durant T, Stanton MC, Biritwum NK, Rollinson D, Ombede DRE, Tchuem-Tchuenté LA. Towards interruption of schistosomiasis transmission in sub-Saharan Africa: developing an appropriate environmental surveillance framework to guide and to support ‘end game’ interventions. Infect Dis Poverty 2017;6:10. pmid:28088239