There is no author summary for this article yet. Authors can add summaries to their articles on ScienceOpen to make them more accessible to a non-specialist audience.
Abstract
This hookworm, uncommonly found in humans, has a higher cure rate than that for
Necator americanus hookworm.
Abstract
During 2012–2015, US-bound refugees living in Myanmar–Thailand border camps (n = 1,839)
were surveyed for hookworm infection and treatment response by using quantitative
PCR. Samples were collected at 3 time points: after each of 2 treatments with albendazole
and after resettlement in the United States. Baseline prevalence of
Necator americanus hookworm was 25.4%,
Ancylostoma duodenale 0%, and
Ancylostoma ceylanicum (a zoonosis) 5.4%. Compared with
N. americanus prevalence,
A. ceylanicum hookworm prevalence peaked in younger age groups, and blood eosinophil concentrations
during
A. ceylanicum infection were higher than those for
N. americanus infection. Female sex was associated with a lower risk for either hookworm infection.
Cure rates after 1 dose of albendazole were greater for
A. ceylanicum (93.3%) than
N. americanus (65.9%) hookworm (p<0.001). Lower
N. americanus hookworm cure rates were unrelated to β-tubulin single-nucleotide polymorphisms at
codons 200 or 167.
A. ceylanicum hookworm infection might be more common in humans than previously recognized.
Introduction Ascaris lumbricoides (roundworm), Trichuris trichiura (whipworm) and Necator americanus/Ancylostoma duodenale (hookworms) are the most common species of soil-transmitted helminths infecting humans worldwide. More than a billion people are infected with at least one species and 300 million are estimated to have severe infections with more than one of these parasites [1]. These infections are endemic in tropical and sub-tropical regions of the developing world and are associated with poverty, lack of clean water, and poor sanitation [2]. School-age children are the most at risk of infection with STHs and early childhood infections contribute significantly to debilitation [3]. Infected children can be malnourished and experience stunting growth and intellectual retardation, with cognitive and educational deficits [1]. Because of all these characteristics and according to estimates by the World Health Organization (WHO) [4], STHs are included in the group of the so-called neglected tropical diseases (NTDs). The main intervention to control STH infections at a community level is based on periodic mass drug administration (MDA) of the benzimidazole (BZ) anthelmintic (AH) drugs, albendazole (ABZ) or mebendazole (MBZ) [5], that reduce the prevalence and intensity of infections [6]. However large-scale chemotherapy programmes with these drugs have the potential to exert selection pressures on the causing parasites, which may favour the development of drug resistance. Recently, expansion of MDA programmes for STHs have highlighted the need to monitor for the possibility that resistance may develop [7], [8]. This development would have important adverse consequences on the benefits provided by deworming programmes [9], [10]. BZ resistance in other parasitic nematodes is caused by a single nucleotide polymorphism (SNP) in the β-tubulin gene at codon positions 200 (T→A), 167 (T→A), or 198 (A→C) [11]–[13]. The frequency and relative importance of these different SNPs varies among the nematode species studied to date [13]–[15]. Molecular markers have been developed to identify SNPs in A. lumbricoides, T. trichiura and hookworms and wild-type and mutant-type control plasmids have been constructed to obtain the genotype profiles of “susceptible-type” parasites (which do not have mutations at position 200, 167 and 198) and “mutant-type” parasites (that contain mutations in one of these three positions) [16]. The codon 200 polymorphism has been identified in T. trichiura populations collected from untreated subjects in Kenya and from treated subjects in Panama [17], and in hookworms collected in Haiti in an area periodically treated with ABZ [16]. In the present report, our aim was to investigate the field efficacy of ABZ against STH infections in countries where polymorphism in β-tubulin had previously been identified, and to assess the frequency of each SNP prior to, and after ABZ treatment within each country and for each of the three main nematode species. Materials and Methods Study Design and Study Areas in Haiti, Kenya and Panama Three cross-sectional studies were carried out between 2008 and 2010 in three countries located in different geographical areas, including Haiti in the Caribbean, Kenya in East Africa and Panama in Central America. These studies were integrated into the following control and evaluation programmes of ABZ use: a national health programme to fight lymphatic filariasis (LF) and intestinal worms in Haiti (in collaboration with the Hôpital Ste. Croix and the Centers for Disease Control (CDC)); school health and nutrition programmes in Kenya, and health and nutrition programmes in Panama. Thus, each study included a treatment with ABZ, stool examinations, and genotyping of the β-tubulin gene in eggs collected before and after ABZ treatment to assess the drug efficacy against STHs and to examine at each time point the frequency of possible SNPs associated with ABZ resistance in nematodes of veterinary importance. The study designs in Kenya, Panama and Haiti (Figure 1) were different because each study was part of a separate national and local MDA programme. Thus, different protocols were applied rather than a purposely designed multi-study site single protocol. 10.1371/journal.pntd.0002247.g001 Figure 1 Flow chart of the study designs in Haiti, Kenya and Panama. *The number of eggs genotyped is different from the number of eggs collected. This is due to failures in DNA extraction, PCR amplification or Pyrosequencing. x̄ represents the mean number of eggs sampled per host; SD = standard deviation. The dashed lines mean that the connection between the two boxes is not direct. Study in Haiti The study in Haiti included stool collections prior to, and after a drug treatment given in 2009. In Haiti, the study sites were located in the West and Southeast Haitian departments known to be naive for MDA with ABZ. Individuals from five endemic communities who reached the inclusion criterion (older than two years old) were randomly selected. All potential participants were informed by a trained community leader of the purpose and methods of the study and gave their oral consent (in the case of children, consent was obtained from a parent). A total of 353 stool samples were collected and analyzed prior to ABZ treatment. Treatment was then distributed to all people of each community (400 mg ABZ and 6 mg/kg diethylcarbamazine (DEC)). ABZ was supplied as a donation from GlaxoSmithKline (GSK). Participants in the study were not observed during the treatment administration. However, the compliance to treatment was evaluated post-treatment by a questionnaire and relied on self-reported information of each participant, and on pre and post-treatment egg counts. Samples from participants who were not treated with ABZ were not included in the analysis. Follow-up faecal samples (n = 317) were collected two weeks after the drug treatment. Study in Kenya School-aged children from two schools located in Kwale District were enrolled in the study. Although the schools were selected because at the time of the study they were not involved in any other STH deworming programmes, pupils had received ABZ and DEC six month prior to our study for LF treatment. Stool samples were collected prior to ABZ treatment in both schools (n = 128) and then all children received 400 mg ABZ under the supervision of the school teachers. The drug was provided by the Kenyan Ministry of Health and was manufactured by GSK. Seven days later, a follow-up collection was undertaken (n = 92). In one of the schools, the stool collections before and after ABZ treatment were done over two consecutive days and the results averaged. Study in Panama The field study was conducted in the region of Comarca Ngäbe-Buglé [18]. Stool samples were collected from pre-school children at different time points over a period of 16 months from July 2008 to October 2009. Two treatments with ABZ, as a suspension (200 mg: one to two years of age; 400 mg: three to five years of age) were distributed, once in 2008 and the second. nine months later in 2009 [19]. The drug was provided by the Panamanian Ministry of Health. In the present study, we only analyzed the faecal samples collected prior to treatment and three weeks after treatment in 2009 since the baseline prevalence for STHs was very low in the 2008 samples. Prior to this second ABZ treatment, stool samples were collected from 270 pre-school aged children and follow-up samples were available from 222 children. Children were observed while received the ABZ treatment, and only samples from treated children were included in the analysis. Ethics Statement In the three studies, all instances of consent were informed. Ethical approval for Haiti was obtained from the Institutional Review Board of the Centers for Diseases Control and Prevention, Atlanta, Georgia, US (Dr. Patrick Lammie), and the Ethics Committee of the Hôpital Ste Croix, Haiti and included the collection of stool samples, examination of stool samples for helminth eggs and DNA analysis of helminth eggs. Informed consent was obtained from all adult participants and from parents or legal guardians of all minors. A parent or legal guardian gave consent in every case of child participation. Based on past experience, it was likely that some people in the communities would not be able to read. A waiver of written informed consent on the basis that the research presented no more than minimal risk of harm to the subjects and involved no procedures for which written consent is normally required outside the research context in this setting, was requested and approved. The use of oral consent was previously approved by the IRB. Subjects were offered a written copy of the IRB approved consent form. The contents of the approved consent form were explained to each household, following which verbal consent was obtained. The reader of the consent form and a witness signed a copy of the form to indicate the subject's agreement. The study in Kenya was approved by the Kenya Medical Research Institute (KEMRI) Ethical Review Committee. The informed consent process, which was approved by the KEMRI IRB specifically stated that the study was included in the national programme on surveillance of disease control and followed established government procedure. The school heads organized meetings with their parent teacher association to obtain agreement for the project. An information sheet was provided. It was emphasized that the participation in the study by children was voluntary and that they may refuse to participate. All children were treated with ABZ and praziquantel in the study. Finally, all samples were anonymised. As part of a national control programme, informed consent was not documented for each individual. The content of the study was verbally explained in detail to parents/guardians, teachers and children from each school. It was done orally as it was necessary to have this explained in local languages. The use of oral consent was previously approved by the IRB. A parent or legal guardian gave oral informed consent in every case of child participation. The study in Panama was approved by the McGill University Review Board in Canada, the Instituto Conmemorativo de Gorgas and local indigenous leaders in Panama. Written informed consent was obtained from primary caregivers for their own participation as well as that of their children. They were provided with an explanation of the study, its significance, and of participant requirements and rights. They were given an opportunity to ask questions in Spanish and in the local language. A parent or legal guardian gave consent in every case of child participation. Stool Sample Procedure and Examination In all studies and for each intervention (before and after ABZ treatment), labelled containers were distributed to each participant and collected the following morning from the community leader in Haiti, schools in Kenya, and each participant's home in Panama. Three diagnostic techniques: McMaster, Kato-Katz, and FLOTAC were used to identify STH eggs and determine the number of egg per gram (epg) in the faecal material collected before and after ABZ treatment. In Haiti we used a modified McMaster technique on each collected sample. One gram of faeces was suspended with water and the solution was stirred until it was completely broken apart. The mixture was poured through surgical gauze into a centrifuge tube. After centrifugation for 10 min at 15,000 rpm, the supernatant was poured off and the tube containing the sediment was filled with saturated sucrose solution and then gently stirred. After 10 min, an aliquot of the flotation fluid from the upper surface of the solution was transferred into each compartment of a McMaster chamber. The eggs were counted in both chambers using a low power objective (×10). The number of epg of faeces was obtained by multiplying the total number of eggs counted in the two chambers by 50 [20]. In Kenya, the McMaster (as described above) and Kato-Katz techniques [21] were used on samples from one school and Kato-Katz alone was used in the other school. For Kato-Katz, the number of eggs counted was multiplied by 24 to obtain the epg [21]. In Panama, the identification of eggs and the assessment of epg were performed using Kato-Katz [21] and FLOTAC [22], [23], as previously described [19]. One of the primary objectives of the study was to assess the frequency of SNPs associated with resistance to ABZ (seen in the veterinary nematode Haemonchus contortus) in the β-tubulin gene of STHs collected from untreated and treated subjects. Thus, from the three studies, eggs from positive subjects were recovered using a saturated sucrose solution with centrifugation, and recovered eggs were preserved in 70% alcohol until use for molecular analysis. Examination of the β-tubulin Gene and Assessment of SNPs DNA Isolation from Individual Eggs of STHs Eggs previously preserved in 70% alcohol were washed in distilled water and separated by species under a dissecting microscope. Individual eggs were isolated with a 10 µl pipette and then placed into a PCR tube. Genomic DNA of each individual egg was then extracted according to a protocol elaborated by Lake and colleagues [24] and adapted to STHs. One hundred millilitres of lysis buffer (50 mM KCl, 10 mM Tris pH 8.3, 2.5 mM MgCl2, 0.45% Nodidet P-40, 0.45% Tween 20, 0.01% gelatine) was previously prepared, and an aliquot of 1 ml of this lysis buffer was taken to which 10 µl proteinase K (10 mg/ml) and 10 µl β-mercaptoethanol (Sigma-Aldrich, St Louis, MO, USA) were added. The solution of lysis buffer with proteinase K and β-mercaptoethanol was placed at −20°C for one week. Individual eggs were placed into a labelled PCR tube and 15 µl of the mixture was added. The tube was subsequently frozen at −80°C for 30 min and then incubated at 60°C for 2 h. This procedure of freezing and heating the egg solution helps the digestion of the eggshell to facilitate DNA extraction. The proteinase K was inactivated by heating to 94°C for 15 min. Genotyping of β-tubulin Gene Positions 167, 198 and 200 in A. lumbricoides, T. trichiura and Hookworms After DNA extraction, all samples were subjected to PCR to amplify small fragments surrounding the SNP at codon positions 167, 198, and 200 in the β-tubulin genes of A. lumbricoides (FJ501301.1, GenBank acc. number), T. trichiura (AF034219.1 GenBank acc. number) and hookworms (EF392851.1 GenBank acc. number). All primers used for genotyping with the Pyrosequencer were designed with the PyroMark Assay Design Software (Qiagen, version 2.0) to amplify single egg DNA. To allow for nested PCRs on single eggs, some primers were different from those used in our previous study [17] in which the DNA concentrations were higher. For PCR, 2 µl of lysate was used as template in a 20 µl reaction containing 2 µl 10×PCR buffer, 1 µl [50 mM] MgSO4, 1 unit Platinum Taq DNA Polymerase High Fidelity, 1 µl of sense and antisense primers [10 µM] (details in Supplementary Table S1) 1 µl dNTP mix [10 µM], and distilled water to 20 µl. PCR reactions were performed with the following cycling parameters: 94°C for 3 min followed by 30 cycles of 94°C, 59°C for 45 s and 68°C for 1 min with a final extension step at 68°C for 5 min. A second PCR reaction using the earlier PCR product as template was necessary to visualize the product on an agarose gel. The reaction contained 5 µl 10×PCR buffer, 2 µl [50 mM] MgSO4, 1 unit Platinum Taq DNA Polymerase High Fidelity, 1 µl of sense and antisense primers [10 µM], 1 µl dNTP mix [10 µM], and distilled water to give a final volume of 50 µl. The same primers were used for each species amplified except for the antisense primer for codon position 198–200 of A. lumbricoides (5′ CAGATGTCGTACAAAGCCTCATT 3′, position). For genotyping, all antisense primers were biotinylated at their 5′end as described [17]. The amplification conditions were 94°C for 3 min followed by 40 cycles of 94°C, 59°C for 45 s and 68°C for 1 min with a final extension step at 68°C for 5 min. Sequencing primers for SNP analysis were used to genotype all PCR products at the codon positions 167, 198 and 200 in a Pyrosequencer (Biotage AB, Charlottesville, VA, USA) (Supplementary Table S1). Parasites genotyped did not necessary come from the same subjects sampled before and after treatment. Eggs were extracted from as many hosts as possible, though finding sufficient number post-treatment was restricted by low egg density. Therefore all eggs, from the same community or school, were pooled, before analysis of genotypes of individual eggs. The number of eggs isolated and genotyped are presented in the flow-chart (Figure 1). Statistical Analysis The treatment efficacy on A. lumbricoides, T. trichiura and hookworms was evaluated by the egg reduction rate (ERR) for each diagnostic method applied in each country: McMaster (in Haiti, and Kenya), Kato-Katz (In Kenya and Panama) and FLOTAC (in Panama). The ERR was calculated, at the group level, as the ratio of the difference between the arithmetic mean of the pre- and post-treatment faecal egg count (FEC) to the pre-treatment arithmetic mean, expressed as a percentage, i.e. ignoring individual variability [25]. The negative individuals at baseline were still sampled at the post-treatment collection and the resulting data were included for the calculation of ERR. Uninfected subjects were included in the mean of the FEC. Confidence intervals of each ERR estimate were determined using a bootstrap resampling method (with replacement) over 10,000 replicates in R (version 2.15.0, Vienna Austria, http://www.R-project.org). Genotype frequencies of SNPs at positions 167, 198 and 200 in A. lumbricoides, T. trichiura and hookworms obtained at the pre-treatment collections were compared with the genotype frequencies of the same SNPs obtained at the post-treatment collections using Fisher's exact test within GraphPad Prism (GraphPad software, San Diego, CA, USA). Deviation from Hardy-Weinberg equilibrium (HWE) was analyzed for the β-tubulin gene at position 200 in T. trichiura using Arlequin version 3.1 software [26], where the p-value was calculated based on the Markov-chain method [27]. Deviations from the HWE were not determined for A. lumbricoides or hookworm for reasons explained below. When a departure from HWE was observed in T. trichiura, we estimated the maximum likelihood frequency of a null allele at position 200. This estimate was calculated using an Expectation-Maximization Algorithm of Dempster and colleagues [28] (EM Algorithm, http://132.206.161.123/em.html). Results Genotype Frequencies for the β-tubulin Gene in T. trichiura In T. trichiura, the three codon positions 167, 198, and 200 were found to be polymorphic in samples collected from untreated and treated subjects in Haiti, Kenya, and Panama. In Haiti we analyzed 65 individual T. trichiura eggs from 30 untreated subjects and 38 from 14 treated subjects. We recorded 11% and 47% experimental failure (includes DNA extraction, PCR amplification and Pyrosequencing) in pre- and post-treatment samples, respectively. Before treatment the T→A SNP at codon position 200 (SNP200) was identified at low frequency; 3.1% of individual eggs genotyped were homozygous resistance-type (AA) and 23.1% were heterozygous (TA). After treatment, there was a significant increase in the frequency of the homozygous resistance-type, from 3.1% to 55.3% (p<0.001), and a statistically significant decrease of the homozygous susceptible-type, from 75.4% to 21.1% (p<0.0001) (Figure 2A). The A→C SNP at codon position 198 was also found and showed a statistically significant change in frequency following treatment. However, the changes at codon 198 were less pronounced than at codon 200. After treatment, there was a significant increase of homozygous resistant-type from 3.1% to 13.2% (p<0.001), and a significant decrease of homozygous susceptible-type from 73.8% to 63.2% (Figure 3). 10.1371/journal.pntd.0002247.g002 Figure 2 Genotype frequencies of the β-tubulin gene position 200 in T. trichiura from Haiti and Kenya. Genotype frequencies of T. trichura collected in Haiti (A) and in Kenya (B); Number of individual T. trichiura eggs genotyped according to the available material, in Haiti was 65 in the untreated group (pre-Tx) and 38 in the treated group (post-Tx), in Kenya was 40 in the untreated group and 90 in the treated group. Sequences were diploid, TT indicates the homozygous susceptible-type TTC/TTC, TA the heterozygous TTC/TAC and AA, the homozygous resistance-type TAC/TAC; Tx = treatment, *Indicates a significant difference (p<0.001) in genotype frequency between the pre- and post treatment groups. P-values were obtained by Fisher's exact test. 10.1371/journal.pntd.0002247.g003 Figure 3 Genotype frequencies of the β-tubulin gene position 198 in T. trichiura from Haiti before and after ABZ treatment. Genotype frequencies of T. trichura collected in Haiti. The number of individual T. trichiura eggs genotyped, according to the available material, was 65 in the untreated group (pre-Tx) and 38 in the treated group (post-Tx). Sequences were diploid, AA indicates the homozygous susceptible-type GAA/GAA, AC the heterozygous GAA/GCA and CC, the homozygous resistance-type GCA/GCA; Tx = treatment, *Indicates a significant difference (p<0.001) in genotype frequency between the pre- and post treatment groups. P-values were obtained by Fisher's exact test. In Kenya, 40 individual T. trichiura eggs from 20 untreated subjects and 90 eggs from 31 treated subjects were genotyped. We recorded 30% and 5% experimental failure in pre- and post-treatment samples, respectively. Only the codon 200 SNP was polymorphic. The same trend observed in Haiti was seen with a statistically significant increase, from 51.3% to 68.5% (p = 0.019) of homozygous resistance-type, and a significant decrease, from 48.6% to 21.4% (p = 0.019) of homozygous susceptible-type after treatment (Figure 2B). In Panama we genotyped 19 “pre-treatment” T. trichiura eggs from 10 subjects, which were collected 9 months after an earlier treatment, and 49 post-treatment eggs from 21 subjects collected 3 weeks after this second ABZ treatment. We observed homozygous resistance-type (78.9%) and homozygous susceptible-type (21.1%), at codon 167, in the pre-treatment collection. For codons 198 and 200, 84% of the pre-treatment T. trichiura egg samples failed the PCR amplification and therefore it was not possible to assess the genotype frequencies for these positions. At the post-treatment collection, the codon 167 polymorphism was still present in the treatment survivors as homozygote resistance-type in 16.3% of individual eggs genotyped. For codons 198 and 200, we observed a predominance of homozygous susceptible-type genotypes, 97.6% and 88.1% respectively; and a low frequency of heterozygous, 2.4% (for each position) and homozygous resistance-type, 9.5% (for codon 200). Genotype Frequencies of β-tubulin Gene Positions 167, 198 and 200 in A. lumbricoides In A. lumbricoides, the codon position 167 of the β-tubulin gene was polymorphic in parasites collected from untreated and treated subjects in Haiti, Kenya and Panama, whereas the codon positions 198 and 200 were monomorphic. The SNP at position 167 identified was confirmed by conventional Sanger sequencing and real time PCR (data not shown). In Haiti, we genotyped 37 individual A. lumbricoides eggs from 13 untreated subjects and five eggs from three treated subjects. We recorded 49% and 28% experimental failure in pre- and post-treatment samples. Prior to treatment, homozygous resistance-type (40%) and heterozygous (60%), at codon 167, were present in the population; however, only heterozygotes were detected after treatment. In Kenya, 22 individual A. lumbricoides eggs from 6 untreated subjects and 19 eggs from 4 treated subjects were genotyped. We recorded 15% and 17% experimental failure in pre- and post-treatment samples, respectively. The predominant genotype frequency identified was the homozygous resistance-type (72.7%) at codon 167. After treatment, there was no significant difference in the homozygous resistance-type frequency; however, as observed in A. lumbricoides collected in Haiti, at codon 167 there was a statistically significant increase of heterozygotes from 4.5% to 21.1% (p<0.001), and also a significant decrease of homozygous susceptible-type from 22.7% to 5.3% (p<0.001). In Panama, 53 individual eggs were genotyped from 28 untreated subjects and 70 eggs from 20 treated subjects. We recorded 57% of experimental failure only in pre-treatment samples. As seen previously in Haiti and Kenya, the most abundant genotype at codon 167 was homozygous resistance-type (97.7%). In this population the lowest genotype frequency was identified as homozygous susceptible-type (2.3%). After treatment, the percentage of both genotypes did not change significantly (96.7% and 3.2%, respectively) (Table 1). 10.1371/journal.pntd.0002247.t001 Table 1 Genotype frequencies of β-tubulin position 167 in A. lumbricoides before and after ABZ treatment. Frequency (%) Haiti Kenya Panama Genotypes Pre-Tx4 (n = 37) Post-Tx (n = 5) Pre-Tx (n = 22) Post-Tx (n = 19) Pre-Tx (n = 53) Post-Tx (n = 70) TT 1 22.7a 5.3a TA 2 60 100.0 4.5b 21.1b 2.3 3.2 AA 3 40 72.7 73.7 97.7 96.8 1 TT = homozygous susceptible-type TTC/TTC. 2 TA = heterozygous TTC/TAC. 3 AA = homozygous resistance-type type TAC/TAC. 4 Tx = treatment. The numbers in parentheses indicate the number of individual eggs genotyped. The letters (a,b) indicate significant difference (p<0.001) between the genotypes of the pre- and post treatment groups. P-values were obtained from Fisher's exact test. Genotype Frequencies of β-tubulin Gene Positions 167, 198 and 200 in Hookworms In hookworms, the codons 167 and 198 of the β-tubulin gene were monomorphic in all samples genotyped from Haiti, Kenya and Panama. Codon 200 polymorphism (TAC) was detected in 2 eggs collected in Kenya. In Kenya, 86 individual eggs from 28 untreated subjects and 127 eggs from 34 treated subjects were genotyped. We recorded 19% and 4% experimental failure in pre- and post-treatment samples, respectively. In the pre-treatment collection we identified homozygous resistance-type at low frequency (2.3%) and a predominance of homozygous susceptible-type (97.7%). After treatment the frequencies did not change significantly. In Haiti, we examined 84 hookworm eggs from 31 untreated subjects and 14 from five treated subjects and we did not identify any polymorphism at the SNP sites of interest. We recorded 19% and 63% of experimental failure for the samples from the pre- and post-treatment, respectively. In Panama, all 23 eggs analyzed from nine untreated subjects and 59 from 29 treated subjects were homozygous susceptible-type for all three positions. We recorded 71% and 24% of experimental failure for the samples from the pre- and post-treatment, respectively. Drug Efficacy The arithmetic means of the faecal egg count (FEC) per gram are presented in Table 2. The standard error of the mean obtained shows a high variability of the FEC in the treated and untreated populations. 10.1371/journal.pntd.0002247.t002 Table 2 Faecal egg count reduction rates for A. lumbricoides, T. trichiura and hookworms. Arithmetic mean of FEC1 at pre-treatment (95% CI2) Arithmetic mean of FEC at post treatment (95% CI) ERR3 (%) (95% CI) Country and Parasite McMaster Kato-Katz FLOTAC McMaster Kato-Katz FLOTAC McMaster Kato-Katz FLOTAC Haiti, n 353 4NA NA 317 NA NA Ascaris 370.0 (50.1–928.8) NA NA 1.6 (0–5.0) NA NA 99.9 (99.5–100) NA NA Trichuris 48.4 (25.7–76.3) NA NA 24.5 (5.9–51.1) NA NA 49.7 (0–88.4) NA NA HW5 209.1 (107.7–338.3) NA NA 2.9 (0.8–6.2) NA NA 98.6 (96.1–99.7) NA NA Kenya, n 24 104 NA 42 92 NA Ascaris 2737.5 (12.5–8166.7) 330.5 (12.3–932.2) NA 73.8 (0–221.4) 73.6 (0–195.1) NA 97.3 (0–100) 80.3 (0–100) NA Trichuris 1650.0 (62.5–4750.0) 161.5 (69.2–282.1) NA 217.1 (32.9–537.8) 145.2 (39.2–285.3) NA 86.8 (0–98.7) 10.1 (0–78.1) NA HW 854.2 (64.6–2139.6) 279.6 (171.7–414.8) NA 3.6 (0–9.5) 8.8 (2.3–17.5) NA 89.9 (0–96.1) 96.8 (92.5–99.2) NA Panama, n NA 92 108 NA 65 223 Ascaris NA 8346.2 (4914.9–12414.7) 1900.8 (1060.4–2840.2) NA 3590.0 (451.3–7913.4) 171.9 (50.8–322.9) NA 60.0 (0–95.1) 89.8 (75.8–97.3) Trichuris NA 280.8±170.1 40.2 (11.1–78.1) NA 1.4 (0–3.5) 14.1 (5.2–25.3) NA 99.5 (89.7–100) 65.1 (0–89.1) HW NA 192.5 (1.8–478.1) 72.6 (31.1–125.8) NA 0.9 (0–2.3) 37.9 (8.3–90.4) NA 99.8 (96.6–100) 47.8 (0–89.9) 1 FEC, faecal egg count, (egg per gram of stool). 2 CI, confidence interval. The confidence intervals of the egg reduction rate were calculated using a bootstrap resampling method [45]. 3 ERR, egg reduction rate. The ERR estimates were obtained by dividing the difference between the arithmetic mean of the pre- and post-treatment FEC (at group level, ignoring individual variability) by the mean of the pre-treatment mean count (see text). 4 NA = Not applicable. 5 HW = Hookworm. The ERR estimates calculated for A. lumbricoides, T. trichiura and hookworms for each diagnostic test applied in Haiti, Kenya and Panama are summarized in Table 2. The ERR from the same country estimated using two different diagnostic methods are not directly comparable as the number of samples tested (and therefore the hosts making up the ERR) were different. A direct comparison of the different diagnostic methods per se is beyond the scope of this paper and has been previously discussed elsewhere [29], [30]. In Haiti where only McMaster was applied, the ERR for A. lumbricoides was the highest (99.9% (95% CI 99.5–100.0)) followed by hookworms (98.6% (95% CI 96.1–99.7)) and T. trichiura (49.7% (95% CI 0.0–88.4)). In Kenya the ERR estimates differed between the Kato-Katz and McMaster methods. The highest ERR was obtained for hookworms, with Kato-Katz at 89.95% (95% CI, 0.0–96.1%) and McMaster at 96.8% (95% CI, 92.5–99.2%). This was followed by A. lumbricoides with 97.3% (95% CI, 0.0–100.0%) and 80.3% (95% CI, 0.0–100.0%), respectively, and finally by T. trichiura with 86.8% (95% CI, 0.0–98.7%) and 10.1% (95% CI, 0.0–78%). The ERR estimate based on Kato-Katz was the one considered as the sample sizes available were much higher (n = 104 and n = 92 at pre- and post-treatment, respectively) than those for McMaster (n = 24 and n = 42 at pre- and post-treatment, respectively). In Panama, the ERR estimates for FLOTAC were the ones considered as the number of samples analyzed was greater than those counted by Kato-Katz [19]. In Panama, the highest ERR estimate was for A. lumbricoides, at 89.8% (95% CI, 75.8–97.3%), followed by T. trichiura with 65.1% (95% CI, 0–89.1%), and finally by hookworms with 47.8% (95% CI, 0–89.9%). Hardy-Weinberg Equilibrium Guo's Exact Hardy-Weinberg test [27] showed that there was a significant departure from Hardy-Weinberg expectations recorded in T. trichiura collected in Haiti (p = 0.0366) after treatment, and in Kenya before and after treatment (p<0.0001) for position 200 of the β-tubulin. This disequilibrium was characterized by a deficiency in the number of heterozygotes (Table 3). The estimated frequency of a null allele at position 200 showed no evidence (χ2 = 3.69) that a null allele was responsible of the paucity of heterozygotes. For positions, 167 and 198, respectively, polymorphism was either not found, or the differences between the frequencies of heterozygotes pre- and post-treatment were not significant. 10.1371/journal.pntd.0002247.t003 Table 3 Observed (HO) and expected (HE) heterozygosity tested under HWE for SNP 200 in T. trichiura. Pre- treatment Post- treatment n1 HO HE p-value n HO HE p-value Haiti 65 0.25 0.24 1.00 38 0.26 0.43 0.0366* Kenya 40 0.15 0.50 <0·0001* 90 0.8 0.38 <0·0001* 1 n = sample size. * Indicates a significant p-value; p-values were obtained from Guo's exact test. Discussion One of the main findings of this study was the identification of a SNP at position 200 in T. trichiura samples collected in Haiti and Kenya at the pre-treatment collections. Thus, these findings suggest that the resistance-type allele (TAC200) already existed in these populations prior to ABZ treatment. This result was consistent with a previous report in which the codon 200 polymorphism was identified in adult T. trichiura from Kenya in an ABZ-naïve population [17]. This SNP was also present in human filarial nematodes from Burkina Faso, in samples obtained from pre-treatment patients at a moderate allele frequency (26.2%) [31]. However, it must be taken into consideration that ABZ is a commonly used AH in endemic countries, and yearly school-based helminth control programmes relying on ABZ and MBZ are also common in countries where our studies were performed [32]–[34]. Indeed, in Kenya, six months prior to our study, children from the selected schools had been treated with ABZ in the context of the LF control programme. In Panama, at the baseline collection (nine months earlier), the codon 200 polymorphism was identified (data not shown), and in a previous study carried out in the same region, the same SNP had been detected in T. trichiura [17]. Finally, in Haiti, in the context of the LF control programme, ABZ is widely distributed yearly in combination with DEC at the community level, in surrounding communities [32]. So, it may be possible that some participants in Panama or Haiti from the studied populations had travelled to treated areas or could have been infected with eggs from family members who had been previously treated with ABZ. Thus, it may be possible that parasites have been previously exposed to the drug and that some selection for the SNPs at codon position 200 had occurred prior to the beginning of MDA for STHs, or that in Kenya, selection at codon 200 may have developed following the previous rounds of ABZ treatment in the population. Other findings were the identification of SNPs in T. trichiura at codon 167 in samples from Panama and at codon 198 in samples from Haiti and Kenya, prior to treatment. The identification of the same resistance-type allele at codon 198, in populations that are geographically separated suggests that this allele could be common in T. trichiura populations, even prior to ABZ treatment. However, the frequency of the resistance-type allele at codon 198, and the identification of the resistance-type allele at codon 167 in only one country suggest that these polymorphisms are less common than the codon 200 polymorphism. In the veterinary parasitic nematode Haemonchus contortus, it has been confirmed that the codon 200 polymorphism, which causes BZ resistance, is the predominant SNP associated with resistance, compared to the codon 167 polymorphism [13], or the codon 198 polymorphism [12], [15]. In addition, it was also interesting to note that the codon 167 polymorphism was always found alone and never associated with the SNPs at positions 198 and 200. This is consistent with observations made in H. contortus suggesting that the polymorphism at position 167 (TAC) does not occur in the same allele as the polymorphism at position 200 (TAC) or 198 (GCA) [12], [15]. The post-treatment examination of the resistance-type allele frequencies of T. trichiura collected in the three different geographical areas suggests that ABZ may be selectively eliminating worms carrying susceptible-type alleles and allowing worms with the resistance-type alleles to survive treatment. Indeed, in Haiti and Kenya one round of treatment significantly increased the frequency of the homozygous resistance-type at codon position 200 (TAC/TAC200) and some heterozygote eggs (TTC/TAC200) were also found after treatment. It was also found that the codon 198 polymorphism persisted after one treatment but at lower frequency compared to the codon 200 polymorphism. Thus, our results suggest that ABZ treatment could select for parasites harbouring the codon 200 (TAC) polymorphism. We tested whether the genotype frequencies of the β-tubulin gene at codon 200 obtained before and after treatment were in HWE. In Haiti, the increase of the homozygous resistance-type after treatment, for codon position 200, correlated with a departure from HWE illustrated by a paucity of heterozygotes and an excess of homozygotes. In Kenya, the β-tubulin at position 200 in T. trichiura was not in HWE before treatment as there was an excess of homozygotes. Different hypotheses may explain the HWE imbalance in genotype frequencies; one possible cause could be selection acting on the parasite populations. The above data strongly suggest that the SNP at β-tubulin codon 200 may be responsible for the intermediate efficacy of ABZ against T. trichiura, where the ERR ranged from 10.1% (Kenya) to 65% (Panama). A poor efficacy of ABZ against T. trichiura has been repeatedly observed across a number of different studies [35]. In addition, it is not uncommon to find considerable variation of ERR estimates (ranging from 0 to 90%) in the literature [36], which, in addition to being explained by the use of different diagnostic methods (see above), is in agreement with the high resistance-type allele frequencies observed at loci 200 and 198 in Haiti and Kenya. If drug susceptibility is determined (in part) by these alleles, then AH resistance might develop relatively quickly in locations such as these. Resistance to BZ drugs is usually recessive [11] so high allele frequencies are required to generate the phenotypic diversity that leads to the poor drug efficacy that has been recorded. Another finding of the study was the identification of the resistance-type SNP (based on comparison with studies in some veterinary nematodes), at codon position 167 in A. lumbricoides before and after treatment. In A. lumbricoides collected in Haiti, Kenya and Panama, the allele containing TAC at codon 167 was identified as homozygous at high frequency prior to ABZ treatment. The post-treatment genotype frequencies indicated that ABZ treatment did not change the frequency of the homozygous, resistance-type, codon 167 in Kenya and Panama, where sample sizes were similar, and heterozygotes remained in the population surviving treatment. It was also interesting to note that similar findings have been published for the veterinary parasitic nematode, Teladorsagia circumcincta [13]. In this case, it was demonstrated that BZ treatment did not change the frequency of the homozygous, resistance-associated codon 167 or that of heterozygotes. In a recent study, it was also found that H. contortus heterozygous at codon 167 were associated with susceptible phenotypes [15]. In contrast, in Haiti where the sample sizes were highly different between the pre- and post-treatment groups, we found that the resistance-type genotype frequency of A. lumbricoides at codon position 167 significantly increased after treatment, but the small sample size at the post-treatment collection made it difficult to draw clear conclusions. The genetic data from A. lumbricoides suggest that the polymorphism at codon position 167 is common and may be naturally present. However, ERR estimates indicate that treatment was successful according to WHO standardized thresholds [37]. This implies that in A. lumbricoides, the mutation at codon position 167 may not have an impact on drug efficacy. In contrast, in intestinal trichostrongylid nematodes of livestock (e.g., H. contortus), the SNP 167 is rarely encountered but the homozygous allele containing TAC does confer BZ resistance in some species [38], [39]. It is also common in these species to have SNP associations between codons 167 and 200 that could confer high levels of BZ resistance [15]. It was interesting to find that in A. lumbricoides, all eggs genotyped were susceptible-types at codon positions 198 and 200. It is possible that should polymorphisms be found at codons 200 or 198 in A. lumbricoides, these polymorphisms may change the sensitivity of the parasite to the drug as has been described in veterinary parasites [12], [15]. It will be informative to investigate the β-tubulin gene of A. lumbricoides in areas where MDA has been implemented to determine whether the position 200 is ever polymorphic and whether it correlates with drug response. Finally, the β-tubulin analysis in hookworms, at codon 200, 167 and 198, revealed the presence of a SNP at codon 200 in samples collected in Kenya at pre- and post-treatment. However, the frequency was very low and the ERR indicated a value consistent with ABZ treatment being successful. In other studies in which a reduced drug efficacy was detected polymorphism at codons 200 and 167 was not detected [40], [41]. The main limitation of the present study was the non homogeneity of the study designs in Haiti, Kenya and Panama. Indeed, the samples and data on response to treatment were obtained as part of different collaborations with ongoing treatment protocols rather than being designed specifically for our study. As such it was not possible to standardize the methods used for egg count estimations and exact time samples were taken after treatment. Furthermore, the studied populations were different in the three studies. In Panama and Kenya it was restricted to pre-school and school age children whereas in Haiti, adults and children were sampled. Furthermore, three diagnostic methods have been applied to quantify STHs. Resulting drug efficacies were variable within the same country, depending on the techniques and sample sizes. In another study, the McMaster and Kato-Katz gave similar results for the ERR [25]. It is important to mention that the sample sizes were different for the two diagnostic tests when performed in the same country. Our goal was not to compare test sensitivity for quantifying egg counts. The results reported here should be viewed as those arising from three different field treatment studies in which we used procedures to assess genotype and have compared the genotypes with assessments of ERR for each of the three separate treatment studies, observed with different methods and sampling schedules. As such the ERR results between the studies should not be compared; what was relevant for our investigations was the assessment of ERR by one or more conventional parasitological methods with the genotype determination for samples from the same communities. The other limitation of the study was the low sample size of individual eggs genotyped due to a high percentage of PCR failure. This was caused by the difficulty in amplifying the DNA of individual eggs, possibly due to the presence of inhibitors from faeces that could prevent the DNA extraction [42]. Also, the low DNA concentration may have affected the success of the DNA amplification and pyrosequencing. In T. trichiura eggs from Panama, we recorded more DNA amplification failure for positions 198–200 than for position 167. One of the primers used to amplify the sequence around the codons 198 and 200 was in part in an intron. Polymorphism in this intron may have resulted in poor amplification in some samples. However, the samples that were successfully genotyped allowed identification of SNPs in T. trichiura which may be associated with the ability of worms to survive ABZ treatment. These results stress the importance of implementing, as integral part of MDA programmes, regular monitoring of drug efficacy and associated parasite genetics for prompt detection of drug resistance. In this study, we estimated the ERR at the group level; it is a more accurate estimation of overall ERR because it can include negative individual values. The estimation of the ERR on a group basis by bootstrapping gave a range of [0–89%] in hookworms, whereas when estimated by Halpenny and colleagues [19] on individual basis the ERR was 89%. Likewise, there is an urgent need to consider control strategies that will maintain a low resistance-type SNP prevalence and a high drug efficacy. Strategies have been described to maintain a high level of STH control while also delaying the possible development of drug resistance [7]. Drug combination has been used to combat the problem of AH resistance in veterinary parasites and human parasites [43], [44]. The pros and cons of various combinations should be carefully considered so that the chosen combination will have both a positive impact on drug efficacy and on a reduction of genetic selection for resistance. Conclusion and Future Directions Comparable data on pre- and post-treatment SNP frequencies were obtained at the different study sites, despite differences in the study designs between sites. However, in the context of monitoring and surveillance of STH control programmes, it will be important to conduct multiple studies according to standardized multi-centre protocols in order to allow drug efficacies to be compared. Supporting Information Table S1 Primers used for amplification of STHs β-tubulin SNP regions and for pyrosequencing. (DOCX) Click here for additional data file.
Human infections with Necator americanus and Ancylostoma duodenale hookworms continue to be recognized as a leading cause of iron deficiency anemia and protein malnutrition in developing countries ( 1 ). On the basis of parasitologic surveys of fecal samples, hookworms are estimated to infect 576–740 million persons globally, and over half of the infections occur in Asia and the Pacific regions ( 2 ). Recent molecular-based epidemiologic surveys have shown Ancylostoma ceylanicum to be the second most common hookworm species infecting humans in Asia. In Thailand, Laos, and Malaysia, 6%–23% of persons positive for hookworm eggs were infected with A. ceylanicum helminths ( 3 – 6 ). There are an estimated 19–73 million A. ceylanicum hookworm–infected persons in regions where this zoonotic helminth is known to be endemic ( 7 ). Dogs and cats act as natural reservoirs for hookworm transmission to humans, and the prevalence of A. ceylanicum hookworms in these animals ranges from 24% to 92% in the Asia-Pacific region ( 6 , 8 – 10 ). Much like anthroponotic helminths, A. ceylanicum hookworms have the potential to produce clinical symptoms of ground itch (a pruritic papular hypersensitivity response caused by the entry of helminths into the skin), epigastric pain, diarrhea, and anemia in humans ( 11 – 15 ). However, despite these reports, relatively little is known about the clinical significance and infection dynamics of this zoonotic hookworm in humans, dogs, and cats. Differentiation of the genus of hookworms infecting humans is imperative because each species varies in its biology, life cycle, pathophysiology, and epidemiology, and these differences have key implications when assessing hookworm-associated illnesses and establishing control measures. The internal transcribed spacer (ITS) –1 and –2 regions and the 5.8S region have been used to detect and characterize hookworm infections directly from eggs in human and animal feces ( 6 , 10 , 16 , 17 ). In addition, sequencing of the cytochrome c oxidase subunit 1 (cox1) gene has been successfully used to establish intraspecies genetic differences of many strongylid nematodes, including hookworms ( 18 – 21 ). The aim of our study was to determine the prevalence, associated risk factors, and infection dynamics of hookworm species infection in humans and dogs living in a rural Cambodian village. To carry out this investigation, we used a combination of conventional parasitologic and molecular epidemiologic approaches. Materials and Methods Study Site and Sample Collection The study was conducted in May 2012 in Dong, a rural village in Rovieng District, Preah Vihear Province, Cambodia. Preah Vihear Province is located in northern Cambodia, bordering Thailand and Laos (13°47′N 104°58′E). The climate is tropical; temperatures are warm and hot all year round, and seasons alternate between dry and wet. Subsistence farming (rice, vegetables, and fish) constitutes the primary source of income for the community. Drinking water is sourced from wells, well pumps, and rain water tanks, and just over half of the households own a latrine. All household electricity is battery or generator powered. Approximately half the households feed semidomesticated, free-roaming community dogs. These dogs are allowed to defecate indiscriminately within the village or outside the homes of their owners. The study protocol was approved by the Ethics Committee of the Canton of Basel and Baselland, Switzerland, and the National Ethics Committee Health Research, Ministry of Health, Cambodia. Dong, the village selected for study, had previously been categorized as having endemic soil-transmitted helminths ( 22 ). According to the treatment guidelines of the Cambodian helminths control program, all children attending primary school in the village were administered albendazole (400 mg) and mebendazole (500 mg) twice a year. At completion of the study, all participants who were found positive for Strongyloides spp. were treated with ivermectin (200 μg/kg body weight), and participants infected with other soil-transmitted helminths were treated with albendazole (400 mg). A cross-section of 67 households was randomly selected from a list provided by the Dong village authority. A total of 218 persons from those households were enrolled in the study. Of the 218 persons, 99 (45.4%) were male. The average age of participants was 30.0 years (range 2–84); female participants were marginally older, on average, than male participants (30.3 vs. 29.8 years of age). On the first day of the study, informed consent was obtained from the enrolled participants, and questionnaires were administered during interviews. Interviews with children (i.e., participants 2–17 years of age) were conducted with the assistance of a parent or legal guardian. All study participants responded to a questionnaire covering demographics, dietary habits, personal hygiene, and level of household income and assets. Prelabeled stool containers were distributed to the 218 study participants for collection of feces on the second morning of the study. Fecal samples were collected from participants’ dogs (N = 94), when applicable. Samples (≈3–5 g) from dogs were collected directly from the rectum at time of the participant interview and placed into a sterile plastic container. If insufficient stool was obtained from a dog, the animal was confined within the owner’s property and resampled on the second morning of the study. All fecal samples were chilled immediately in a cool box and transported to a laboratory in Rovieng Health Center (Rovieng District, Preah Vihear Province) within 2 h after collection. After fecal samples arrived at the laboratory, a minimum of 2 g of each sample was placed into a 15-mL centrifuge tube containing 10% formaldehyde for parasitologic analysis, and 1–2 g of each sample was placed into a 15-mL centrifuge tube containing 2.5% potassium dichromate for molecular analysis. These samples were then shipped at room temperature to the School of Veterinary Science, University of Queensland, Gatton, Queensland, Australia, for further analysis. Parasitologic Procedures All fecal samples were examined by microscope. The relative intensity of hookworm infection, in eggs per gram, was determined by floatation, using a sodium nitrate solution (specific gravity 1.20) ( 23 ). DNA Extraction Genomic DNA was extracted directly from fecal samples by using the PowerSoil DNA Isolation Kit (Mo Bio, Carlsbad, CA, USA) according to the manufacturer’s instructions, with the exception that fecal samples were subjected to a 5-min disruption by using 0.5-mm Zirconia/Silica beads (BioSpec Products, Inc., Bartlesville, OK, USA) instead of the beads provided by the manufacturer. Final elution of DNA was made in 100 μL of elution buffer. The extracted DNA was stored at −20°C until required for PCR amplification. Molecular Characterization of Hookworm species in Humans PCR was conducted by using primers RTHW1F and RTHW1R ( 10 ) in 25-μL volumes; each final reaction contained 1× CoralLoad PCR Buffer (QIAGEN Pty Ltd, Hilden, Germany), 12.5 pmol of each primer, 0.5 U of HotStar Taq DNA Polymerase (QIAGEN), and 2 μL of DNA. The cycling conditions were the same as the published protocol ( 10 ) except for an initial denaturation of 5 min at 95°C. A positive control of N. americanus and A. ceylanicum hookworms and negative controls of distilled water were included in each run. PCR amplicons that were ≈380 bp in size, corresponding to Ancylostoma spp. hookworms, were purified by using the PureLink Quick PCR Purification Kit (Life Technologies, Carlsbad, CA, USA) and submitted to the University of Queensland Animal Genetics Laboratory, Gatton, for bidirectional DNA sequencing. Molecular Characterization of Hookworms species in Dogs PCR–restriction fragment length polymorphism (RFLP) characterization of hookworms from dogs was carried out as described ( 17 , 24 ). In brief, RTGHFI and RTABCR1 primers were used to amplify a 545-bp region of ITS-1, 5.8S, and ITS-2 of A. caninum, A. ceylanicum, and Uncinaria stenocephala hookworms. In a separate PCR, a 673-bp region of an A. braziliense hookworm was amplified by using RTGHF1 and a specific reverse primer, RTAYR1. Both PCR reactions consisted of 1× CoralLoad PCR Buffer (QIAGEN), 12.5 pmol of each primer, 0.2 μL of 20 mg/mL bovine serum albumin, 2 μL DNA, and 1 U of HotStar Taq Polymerase (QIAGEN) in a 25-μL reaction. The cycling conditions were as published ( 17 , 24 ), except for an initial denaturation time of 5 min at 95°C. Amplified PCR product (10 μL; RTGHF1/RTABCR1) was digested with HinFI and RsaI endonucleases in separate reactions at 37°C for 3 h. The RFLP patterns generated by each sample were then compared to the expected RFLP profiles for each hookworm species. PCR and DNA Sequencing of cox1 of A. ceylanicum Hookworm Samples from dogs and humans that were positive for A. ceylanicum hookworms were further characterized to a haplotype level by analysis of the mitochondrial gene (cox1). AceyCOX1F (5′-GCTTTTGGTATTGTA-AGACAG-3′) and AceyCOX1R (5′- CTAACAACATAATAAG-TATCATG-3′) were specifically designed to amplify a 377-bp region of the cox1 gene of A. ceylanicum hookworm. The PCR was carried out in 25-μL volumes, with each reaction containing 1× CoralLoad PCR Buffer, 12.5 pmol of each primer, 0.5 U of HotStar Taq DNA polymerase, and 2 μL of DNA. The cycling conditions were 95°C for 5 min, followed by 50 cycles at 94°C for 30 s, 58°C for 30 s, 72°C for 30 s, and a final extension at 72°C for 7 min. A positive control of A. ceylanicum hookworm and a negative control of distilled water were included in the run. PCR-positive samples were purified by using the PureLink Quick PCR Purification Kit according to the manufacturer’s protocol. Bidirectional DNA sequencing was performed by the University of Queensland Animal Genetics Laboratory. Phylogenetic Analyses DNA sequences were analyzed by using the Finch TV version 1.4.0 trace viewer (Geospiza, Inc., Seattle, WA, USA) and aligned by using BioEdit version 7.2.0 (www.mbio.ncsu.edu/BioEdit/bioedit.html) together with the cox1 gene sequences from the following hookworm species: A. ceylanicum Malaysia isolates (GenBank accession nos. KC247727– KC247745, Pos Iskandar [Human] and Sg Bumbun [Human]); A. caninum and A. duodenale (GenBank accession nos. NC012309 and NC003415, respectively); and A. ceylanicum Thailand genotype (GenBank accession no. KF896595). Neighbor-joining analyses were conducted by using Tamura-Nei parameter distance estimates, and the tree was constructed by using Mega4.1 (www.megasoftware.net). Bootstrap analyses were conducted using 1,000 replicates. Statistical Analyses We used STATA version 12 (StataCorp LP, College Station, TX, USA) for data entry and statistical analyses. The prevalence of hookworm infection was calculated by using descriptive statistics for microscopy and molecular results. A univariate model was used to assess potential risk factors associated with hookworm infection; odds ratios and 95% CI were reported. The level of statistical significance was set at p<0.05. Factors that were significant in univariate analysis were evaluated by multivariate analysis, when applicable. Results Prevalence of Hookworm Infection The prevalence of hookworm infection among the 218 persons tested in Dong village was 26.6% (58/218) as determined by microscopic examination and 57.4% (124/218) as determined by PCR based on amplification of the partial ITS gene. Among dogs, 80.9% (76/94) were positive for hookworms by microscopic examination, and 95.7% (90/94) were positive by PCR based on amplification of the partial ITS gene. Molecular Characterization of Hookworm Species Of the 124 persons with positive samples, 64 (51.6%) harbored A. ceylanicum hookworms; 57 (89.0%) of these infections were single infections. An equal percentage of persons, 64 (51.6%), were infected with N. americanus hookworms, mostly as single infections (59/64 [92.2%]), and 4 (3.2%) persons were infected with A. duodenale hookworms (Table). Table Hookworm species found in humans and dogs, Dong village, Rovieng District, Preah Vihear Province, Cambodia, 2012* Infected host, hookworm species No. (%) positive Humans Necator americanus 59 (47.6) Ancylostoma ceylanicum 57 (46.0) A. duodenale 1 (0.8) N. americanus and A. ceylanicum 4 (3.2) A. ceylanicum and A. duodenale 2 (1.6) N. americanus and A. ceylanicum and A. duodenale 1 (0.8) N. americanus and A. duodenale 0 Total 124 (100.) Dogs A. ceylanicum 81 (90.0) A. caninum 5 (5.6) A. ceylanicum and A. caninum 3 (3.3) A. ceylanicum and N. americanus 1 (1.1) Total 90 (100.0) *The presence of hookworms was determined by PCR amplification of the internal transcribed spacer–1 and –2 regions and the 5.8S region and by DNA sequencing. Of the 90 dogs with positive samples, 85 (94.4%) were infected with A. ceylanicum hookworms, mostly (81/85 [95.3%]) as single infections, and 8 (8.9%) were infected with A. caninum hookworms. One dog was found to be shedding N. americanus eggs (Table). Phylogenetic Analysis of cox1 Gene of A. ceylanicum Of 68 human and 82 dog samples positive for hookworms, 28 (41.2%) and 65 (79.3%), respectively, were successfully amplified at the cox1 gene. Of these, 21 PCR-positive amplicons from human samples and 27 PCR-positive amplicons from dog samples were randomly selected for DNA sequencing and subsequent phylogenetic analysis. The phylogenetic tree distinctly separated into 3 clusters; the A. ceylanicum hookworm isolates grouped together and were genetically distinct from A. caninum hookworm isolates (GenBank accession nos. NC012309 and FJ483518) and A. duodenale hookworm isolates (GenBank accession nos. NC003415 and AJ417718). Within A. ceylanicum hookworm isolates, there was strong bootstrap support (100%) for the division of isolates from various geographic locations into 2 clades. The first clade comprised 4 human isolates, 1 from the current study in Cambodia (Human 19) and 3 previously reported human isolates from Malaysia (GenBank accession no. KC772445; Human [Sg Bumbun]; Human [Gurney] and Human [Pos Iskandar]). The second clade comprised a mix of isolates from humans (n = 20) and dogs (n = 27) from villages in Cambodia; humans (n = 5), dogs (n = 11), and cats (n = 2) from Malaysia ( 21 ); and 1 dog in Thailand (GenBank accession no. KF896595). For human- and dog-derived A. ceylanicum hookworms, representative DNA sequences at each cox1 haplotype were submitted to GenBank under accession nos. KF896596–KF896605 (see sequences marked with asterisks in the Technical Appendix Figure). Age-related Prevalence and Intensity of N. americanus and A. ceylanicum Hookworm Infections The prevalence of N. americanus hookworms peaked in persons 31–50 years of age, whereas the prevalence of A. ceylanicum hookworms peaked in persons 15–20 years age and again in persons 31–50 years of age (Figure). The highest egg intensities for single infections attributed to N. americanus and A. ceylanicum hookworms occurred in persons 21–30 years of age (Figure). Figure Prevalence and intensity (eggs per gram) of Necator americanus and Ancylostoma ceylanicum hookworm infections in humans of different ages in rural Dong village, Rovieng District, Preah Vihear Province, Cambodia, 2012. Risk Factors Associated with Hookworm Infection of Humans and Dogs The results of regression analysis showed an increased risk for hookworm infection in persons who did not wear shoes while defecating (odds ratio 6.0, 95% CI 1.1–28.6; p = 0.038). No significant associations were found between the prevalence and intensity of hookworms by age group, sex, household income, or dietary practices. No risk factors of significance were associated with hookworm infection in dogs. Discussion In this study, zoonotic ancylostomiasis caused by A. ceylanicum hookworms was found to be highly endemic among humans in Dong village, Preah Vihear Province, Cambodia; community dogs were the likely zoonotic reservoir. This finding is in stark contrast to the consistent finding by other molecular-based prevalence studies in the region that N. americanus is the predominant hookworm species in humans, followed by A. ceylanicum and A. duodenale hookworms ( 3 , 5 , 6 , 10 ). PCR proved a superior alternative to microscopy-based techniques for the detection of hookworms in fecal samples ( 25 , 26 ). In Dong village, the prevalence of A. ceylanicum hookworms matched that of their anthroponotic counterpart, N. americanus hookworms, and infections with A. ceylanicum hookworms substantially out-numbered those with A. duodenale hookworms. In addition, most infected persons harbored single-species hookworm infections; just over 10% of hookworm-positive persons had mixed-species infections. These results raise questions about the potential infection dynamics between hookworm species within individual hosts. Our study supports an earlier hypothesis ( 7 ) that anthroponotic hookworms may have a cross-protective role in expelling and preventing the subsequent establishment of A. ceylanicum hookworms via a T helper 2 cell response ( 27 ). The major immunologic action against incoming L3 larvae (infective filariform larvae) and L4 larvae (final larval stage within the intestine) is regulated by the infection itself ( 28 ). Thus, the presence of a stable and long-lived (3–6 years) infection with anthroponotic species ( 29 ) may play a role in providing an unsuitable environment for the establishment of incoming larvae of another closely related (albeit potentially shorter-lived and suboptimally host adapted) species—in this case, A. ceylanicum hookworms. Reduced burdens of anthroponotic hookworm species may also have the added advantage of easing density-dependent intraspecific competition for limited resources within the intestinal niche ( 30 ), leading to an opportunistic establishment of A. ceylanicum hookworms. Although data on the natural life span of A. ceylanicum hookworms in humans do not exist, infections in Dutch servicemen 5 months after their return from New Guinea ( 31 ) suggest that chronic infections with this hookworm may occur. The initiating or causal factor for the emergence of highly endemic levels of monospecific infections with A. ceylanicum hookworms in Dong village remains unclear. In this study, potential causal factors for human infection are likely related to the high levels of A. ceylanicum hookworm infections in community dogs. In rural Malaysia, close contact with community dogs and cats was shown often to be associated with human infection with A. ceylanicum hookworms ( 6 ). In Dong village, dogs were reported to defecate indiscriminately in environments shared with humans, leading to widespread environmental contamination with infective A. ceylanicum hookworm larvae. For humans, defecating while bare foot was shown to be the most significant risk factor for infection with both species of hookworms. Whether these factors, coupled with the administration of preventative chemotherapy, led to an increased opportunity for the A. ceylanicum hookworm to replace the niche of its anthroponotic competitors remains unanswered. Either way, integrated control programs aimed at combining chemotherapeutic interventions with improvements in community hygiene and animal health programs will aid in curbing this potentially opportunistic zoonosis. Molecular epidemiologic data gathered from characterization of the cox1 gene of A. ceylanicum hookworms strongly support previous findings ( 21 ) that A. ceylanicum hookworm isolates from humans and animals formed 2 genetically distinct groups, 1 comprising isolates specific to humans and the other comprising isolates from humans, dogs, and cats. Most A. ceylanicum hookworm isolates from humans in Dong village clustered within the zoonotic haplotype, confirming that transmission from dogs to humans has occurred. Genetic groups inferred by the cox1 gene of A. ceylanicum hookworms were found to be independent of geographic source. Whether the 2 primary haplotypes differ in biologic, epidemiologic, and pathophysiologic characteristic warrants further investigation. The transmission dynamics of A. ceylanicum hookworms in humans of different ages largely paralleled that of N. americanus hookworms: persons 21–30 years of age excreted the highest number of eggs. This highly unexpected finding has key implications. First, this finding suggests that the previous classification of A. ceylanicum as an abnormal and minor hookworm of humans ( 32 ) no longer stands. Second, monospecific infections of humans with <100 A. ceylanicum worms have been reported to cause anemia, even in well-nourished persons ( 14 , 33 ). Thus, attention must be directed to A. ceylanicum infection as a major cause of human illness in areas where this zoonosis is endemic. The zoonotic helminth A. ceylanicum can no longer be classified as an abnormal hookworm of humans. Although previous studies have reported this hookworm’s emergence as the second most common human hookworm species in Southeast Asia, our study demonstrates its ability to infect humans at prevalence and intensity levels at par with that of its anthroponotic competitor, the N. americanus hookworm. We hypothesize that expansion of preventative chemotherapy in the absence of concurrent hygiene and animal health programs is a potential causal factor for the emergence of this zoonosis. Attention must be directed to the effects of A. ceylanicum hookworm infections on human health, and a One Health approach should be adopted for the control of this zoonosis. Technical Appendix Phylogenetic tree of Ancylostoma ceylanicum hookworms from 21 humans and 27 dogs in Cambodia together with reference isolates from Malaysia and Thailand.
Publisher:
Centers for Disease Control and Prevention
ISSN
(Print):
1080-6040
ISSN
(Electronic):
1080-6059
Publication date
(Print):
August
2018
Volume: 24
Issue: 8
Pages: 1472-1481
Affiliations
[1]National Institutes of Health, Bethesda, Maryland, USA (E.M. O’Connell, G. Oduro-Boateng,
S. Harrison, T.B. Nutman);
[2]Centers for Disease Control and Prevention, Atlanta, Georgia, USA (T. Mitchell, D.
Lee, M. Weinberg, W.M. Stauffer);
[3]Smith College, Northampton, Massachusetts, USA (M. Papaiakovou, N. Pilotte, S.A. Williams);
[4]University of Massachusetts, Amherst, Masschusetts, USA (N. Pilotte, S.A. Williams);
[5]International Organization for Migration, Mae Sot, Thailand (P. Sakulrak, D. Tongsukh);
[6]University of Minnesota Medical School, Minneapolis, Minnesota, USA (W.M. Stauffer)
Author notes
Address for correspondence: Elise M. O’Connell, National Institute of Allergy and
Infectious Diseases, 4 Center Dr, B105, Bethesda, MD 20892-6612, USA; email:
oconnellem@
123456niaid.nih.gov
Article
Publisher ID:
18-0280
DOI: 10.3201/eid2408.180280
PMC ID: 6056130
PubMed ID: 30014834
SO-VID: 6e787d3d-0f69-4928-9f1c-8ea27f8baeb3
History
Categories
Subject:
Research
Subject:
Research
Subject:
Ancylostoma ceylanicum Hookworm in Myanmar Refugees, Thailand, 2012–2015
scite shows how a scientific paper has been cited by providing the context of the citation, a classification describing whether it supports, mentions, or contrasts the cited claim, and a label indicating in which section the citation was made.