Prevalence and molecular characterization of Strongyloides stercoralis, Giardia duodenalis, Cryptosporidium spp., and Blastocystis spp. isolates in school children in Cubal, Western Angola
Parasites & Vectors volume 11, Article number: 67 (2018)
Human infections by the gastrointestinal helminth Strongyloides stercoralis and the enteric protozoans Giardia duodenalis, Cryptosporidium spp. and Blastocystis spp. are not formally included in the list of 20 neglected tropical diseases prioritised by the World Health Organization. Although largely underdiagnosed and considered of lower public health relevance, these infections have been increasingly demonstrated to cause significant morbidity and even mortality globally, particularly among children living in resource-poor settings.
In this cross-sectional survey the prevalence, frequency and molecular diversity of S. stercoralis, G. duodenalis, Cryptosporidium spp. and Blastocystis spp. were investigated in a school children population in the province of Benguela (Angola). A total of 351 stool samples were collected during January to June 2015. The presence of S. stercoralis and G. duodenalis was confirmed by qPCR methods. Giardia duodenalis assemblages and sub-assemblages were determined by multilocus sequence-based genotyping of the glutamate dehydrogenase and β-giardin genes of the parasite. Detection and identification of Cryptosporidium and Blastocystis species and subtypes was carried out by amplification and sequencing of a partial fragment of the small-subunit ribosomal RNA gene of both protozoan. Analyses of risk factors potentially associated with the transmission of these pathogens were also conducted.
Prevalences of S. stercoralis, G. duodenalis, Cryptosporidium spp., and Blastocystis spp. were estimated at 21.4% (95% CI: 17.1–25.7%), 37.9% (95% CI: 32.8–43.0%), 2.9% (95% CI: 1.1–4.5%) and 25.6% (95% CI: 21.18–30.2%), respectively. Overall, 64.1% (225/351) of the children were infected by at least one of the pathogens investigated. Sequence analyses of the 28 G. duodenalis isolates that were successfully genotyped allowed the identification of sub-assemblages AI (14.3%), AII (14.3%), BIII (7.1%) and BIV (25.0%). Discordant typing results AII/AIII and BIII/BIV were identified in 7.1% and 14.3% of the isolates, respectively. A total of five additional isolates (17.9%) were identified as assemblage B. Three Cryptosporidium species including C. hominis (70%), C. parvum (20%) and C. canis (10%) were found circulating in the children population under study. A total of 75 Blastocystis isolates were assigned to the subtypes ST1 (30.7%), ST2 (30.7%), ST3 (36.0%), ST5 (1.3%) and ST7 (1.3%), respectively. Children younger than seven years of age had significantly higher risk of infections by protozoan enteropathogens (PRR: 1.35, P < 0.01), whereas being underweight seemed to have a protective effect against these infections (PRR: 0.74, P = 0.005).
The burden of disease attributable to human strongyloidiasis, giardiosis, cryptosporidiosis and blastocystosis in Angola is considerably higher than initially estimated in previous surveys. Surveillance and control of these infections should be jointly tackled with formally considered neglected tropical diseases in order to maximize effort and available resources. Our data also demonstrate the added value of using molecular diagnostic methods in high transmission areas.
Gastrointestinal helminthic and protozoan infections (GHPI) are significant contributors to the global burden of disease, particularly in low- and middle-income countries in tropical regions [1, 2]. GHPI primarily affect children in disfavoured settings with limited or no access to safe drinking water, inadequate sanitation, and substandard housing where poor hygiene practices are common . Because enteric nematodes, cestodes, trematodes and protozoans often occur sympatrically in endemic areas, polyparasitism has been suggested to have a synergistic effect in exacerbating detrimental health outcomes in infected individuals . Thus, GHPI have been linked to malabsorption, malnutrition, stunting, chronic anemia, cognitive impairment and failure to thrive [5,6,7,8]. Despite their unquestionable socio-economic and public health impact , the epidemiology of GHPI is still poorly understood in many regions of the world and only partially addressed in the Global Burden of Disease (GBD) studies .
Soil-transmitted helminths (STHs) including Ancylostoma duodenale, Necator americanus, Ascaris lumbricoides, Trichuris trichiura and Strongyloides stercoralis are estimated to infect almost one-sixth of the global population , with up to 370 million people being infected only with S. stercoralis . Out of the approximately 50 species that comprise the genus Strongyloides, only S. stercoralis and S. fuelleborni are infective to humans, although the latter is clinically less important and has a restricted geographical distribution [12, 13]. As other STHs, S. stercoralis is primary transmitted through contact with soil that is contaminated with the nematode. Humans acquire the infection when filariform larvae come in contact with the skin, penetrate it, and migrate through the body. After reaching the small intestine, larvae mature into adult worms and lay their parthenogenetic eggs. Still within the intestinal lumen, the eggs hatch into non-infective rhabditiform larvae, which are excreted with stool into the environment. The larvae molt twice and then develop into infective filariform larvae re-initiating the cycle. Strongyloidiasis is an exception among helminthic infections because the parasite can reproduce within a human host (autoinfection) resulting in long-term infections and causing hyper-infection in debilitated and immunocompromised patients that can be fatal . On the other hand, enteric protozoans Cryptosporidium spp., G. intestinalis and Blastocystis spp., together with Entamoeba histolytica, are regarded as the most common and important causes of protozoan-diarrhoeal disease in humans globally. These parasites are transmitted via the faecal-oral route either indirectly through ingestion of contaminated water or food or directly by contact with infected persons or animals. Only in developing countries, about 200 million people are estimated to have symptomatic giardiasis , whereas cryptosporidiosis is a leading cause of diarrhoeal death in children younger than 5 years globally, only second after rotaviral enteritis . Regarding Blastocystis spp. it is currently estimated that up to 1 billion humans across the world would be colonized/infected with these protozoan species . Although the clinical significance of Blastocystis spp. is still the focus of intense debate, there is mounting in vitro and in vivo evidence linking the presence of these protozoan species with intestinal (e.g. nausea, anorexia, flatulence, acute or chronic diarrhoea, irritable bowel syndrome) and extra-intestinal (e.g. urticarial) disorders [17, 18]. Moreover, invasive and inflammatory potential of the parasites has also been reported .
Furthermore, Cryptosporidium encompasses at least 30 valid species, of which C. hominis and C. parvum are known to cause more than 90% of documented human infections . The only Giardia species that is pathogenic to humans, G. intestinalis, is currently regarded as a multi-species complex divided into eight (A to H) distinct genetic variants (assemblages) with marked differences in host range and specificity. Of these, zoonotic assemblages A and B are commonly reported to infect humans . Blastocystis has also shown to exhibit extensive genetic diversity, allowing the differentiation of at least 17 subtypes (ST), of which STs 1–9 infect humans . STs 10–17 are predominantly found in non-human species but some of them have been occasionally documented as zoonotic agents .
The current epidemiological situation of GHPI in Angola (southern Africa) remains largely unknown. Only a limited number of mainly microscopy-based studies targeting paediatric [8, 24,25,26] and women  populations have been conducted to date in the country. Obtained data revealed the presence of A. lumbricoides (4–17%), hookworms (0.6–14.0%), T. trichiura (0.3–14.0%) and S. stercoralis (0.3–3.5%) in the surveyed populations. Similarly, G. duodenalis and Cryptosporidium spp. were also identified at prevalence rates of 13–22% and 30% of the investigated individuals, respectively, although typing and sub-typing analyses were not carried out. In an attempt to overcome this situation, we present here PCR-based prevalence data of S. stercoralis, G. duodenalis, Cryptosporidium spp. and Blastocystis spp. in school children in Central Angola. Molecular data on the diversity and frequency of the protozoan species/genotypes investigated are also shown.
Cubal is a municipality in the province of Benguela, in the highlands of the central plateau of Angola. The town is approximately 900 m above sea level and located 146 km from the city of Benguela. It extends over an area of 1452 km2 and has an estimated population of near 320,000 inhabitants, of whom 47% are children aged between 5 and 14. The municipality is administratively divided into four urban communes, namely Cubal, Tumbulo, Capupa and Yambala (Fig. 1). The dominant biome in Angola is the savannah. Like the rest of tropical Africa, Angola presents pronounced alternate wet-dry seasons of six-months. The interior highlands have a mild climate with a rainy season (November-April) followed by a dry season (May-October) [27, 28].
Study design and stool sample collection
A cross-sectional study was conducted in primary schools in the four communes of Cubal from January to June 2015. The sample size was estimated at 378 children less than 15 year-old based on a 30% geohelminths prevalence previously reported in Angola [24, 29], a marginal error of 2.5%, and a non-response rate of 15%. Total population figures per commune were taken into account in the calculations. A two-step conglomerate sampling was carried out. First, the four most populated schools in each commune were identified and invited to participate in the present survey. Informative meetings were held with school representatives and parents to explain the goals of the project and the methodological aspects involved. Secondly, children were randomly selected to fill the numbers required in each sampled school. Children who received any antiparasitic drug treatment within the three previous months were excluded. Two consecutive stool specimens were obtained per child, individually collected in sterile polystyrene plastic flasks without any preservative, and suitably labelled. Socio-demographic (gender, age, commune of residence, school) and clinical (body mass index, haemoglobin levels) parameters were also compiled from each participant at the time of sample collection. Body mass index (BMI) values were adjusted for differences in gender and age using growth reference data for children and adolescents, 5–19 year- old . Children whose BMI was less than the 5th percentile were considered underweight . Following World Health Organization (WHO) guidelines, anaemia was defined as haemoglobin level of less than 11 g/dl .
Additionally, standardised data collection spreadsheets were used to identify risk factors potentially associated to the transmission of GHPI. Variables analysed included walking barefoot, household size, water-use practices (source of bath and drinking water, washing raw fruits and vegetables before eating), and contact with domestic animals. Samples were transported to the Microbiology Laboratory of the Nossa Senhora da Paz Hospital (Cubal) for parasite identification by coprological methods on the same day of collection. Aliquots of stool specimens were kept at -20 °C for further molecular analyses.
DNA extraction and purification
Prior to DNA extraction, 1 g of each individual stool sample was re-suspended in 0.9% saline solution and concentrated using the Bioparaprep MINI® device (Leti Diagnostics, Barcelona, Spain), according to the protocol described elsewhere . Total DNA was extracted from ~ 200 mg of concentrated faecal material using the QIAamp® DNA Stool Mini Kit (Qiagen, Hilden, Germany) following the manufacturer’s instructions. Purified DNA samples (200 μl) were stored at -20 °C and shipped to the Parasitology Service at the National Centre for Microbiology (Majadahonda, Spain) for downstream PCR-based diagnostic and genotyping analyses. In order to guarantee the absence of contamination of all the reagents used during the DNA extraction and purification procedure, a water extraction control was routinely included in each sample batch processed.
Molecular detection of Strongyloides stercoralis
Genus-specific primers targeting the small subunit ribosomal RNA (SSU rRNA) gene of Strongyloides spp.  were used to detect the presence of the parasite in a qualitative real-time PCR (qPCR) assay using SybrGreen reagents (Invitrogen, San Diego CA, USA) as described elsewhere . qPCR reactions were performed in a final volume of 25 μl containing 1× Quantimix EasyMaster Mix (Biotools B&M Laboratories, Madrid, Spain), 0.2 μM of each specific primer (Additional file 1: Table S1), 0.5 μl of 50× SybrGreen (Invitrogen, San Diego, CA, USA), and 10 μl of total DNA extracted from stool samples. Purified genomic DNA from Strongyloides venezuelensis L3 was used as positive control. All DNA isolates were assayed in duplicate. A third replicate including 10 ng of S. venezuelensis DNA was also included as internal inhibition control. Negative and no template controls were included in each run. The amplification program consisted of 15 min at 95 °C followed by 50 cycles of 10 s at 95 °C, 10 s at 60 °C and 30 s at 72 °C. DNA amplification and detection of fluorescence at the end of each amplification cycle was performed on a Corbett Rotor Gene™ 6000 real-time PCR system (Qiagen, Hilden, Germany). Data were analysed with Rotor Gene 6000 Series software version 1.7.
Molecular detection and characterization of Giardia duodenalis
Initial detection of G. duodenalis DNA was achieved using a qPCR method targeting a 62 bp region of the SSU rRNA gene of the parasite . Amplification reactions (25 μl) consisted of 3 μl of template DNA, 0.5 μM of primers Gd-80F and Gd-127R, 0.4 μM of probe (Additional file 1: Table S1), and 12.5 μl TaqMan® Gene Expression Master Mix (Applied Biosystems, CA, USA). Detection of parasitic DNA was performed on a Corbett Rotor Gene™ 6000 real-time PCR system (Qiagen) using an amplification protocol consisting on an initial hold step of 2 min at 55 °C and 15 min at 95 °C followed by 45 cycles of 15 s at 95 °C and 1 min at 60 °C. Water (no-template) and genomic DNA (positive) controls were included in each PCR run.
Giardia duodenalis isolates that tested positive by qPCR were subsequently assessed by sequence-based multi-locus genotyping of the genes encoding for the glutamate dehydrogenase (GDH) and ß-giardin (BG) proteins of the parasite. A semi-nested-PCR protocol proposed elsewhere  with minor modifications was used to amplify a ~ 432 bp fragment of the gdh gene. PCR reaction mixtures (25 μl) included 5 μl of template DNA and 0.5 μM of the primer pairs GDHeF/GDHiR in the primary reaction and GDHiF/GDHiR in the secondary reaction (Additional file 1: Table S1). Both amplification protocols consisted of an initial denaturation step at 95 °C for 3 min, followed by 35 cycles of 95 °C for 30 s, 55 °C for 30 s and 72 °C for 1 min, with a final extension of 72 °C for 7 min.
Similarly, a ~ 511 bp fragment of the bg gene of G. duodenalis was amplified using the nested-PCR protocol described by Lalle et al. . PCR reaction mixtures (25 μl) consisted of 3 μl of template DNA and 0.4 μM of the primers sets G7_F/G759_R in the primary reaction and G99_F/G609_R in the secondary reaction (Additional file 1: Table S1). The primary PCR reaction was carried out with the following amplification conditions: 1 step of 95 °C for 7 min, followed by 35 cycles of 95 °C for 30 s, 65 °C for 30 s, and 72 °C for 1 min with a final extension of 72 °C for 7 min. The conditions for the secondary PCR were identical to the primary PCR except that the annealing temperature was 55 °C.
Molecular detection and characterization of Cryptosporidium spp. isolates
The presence of Cryptosporidium spp. was assessed using a nested-PCR protocol to amplify a 587 bp fragment of the SSU rRNA gene of the parasite . Amplification reactions (50 μl) included 3 μl of DNA sample and 0.3 μM of the primer pairs CR-P1/CR-P2 in the primary reaction and CR-P3/CPB-DIAGR in the secondary reaction (Additional file 1: Table S1). Both PCR reactions were carried out as follows: one step of 94 °C for 3 min, followed by 35 cycles of 94 °C for 40 s, 50 °C for 40 s and 72 °C for 1 min, concluding with a final extension of 72 °C for 10 min. Additionally, sub-typing of the isolates identified as C. hominis or C. parvum was attempted at the 60 kDa glycoprotein (gp60) locus following the nested-PCR protocol proposed by Feltus et al. .
Molecular detection and characterization of Blastocystis spp. isolates
Identification of Blastocystis spp. was achieved by a direct PCR protocol targeting the SSU rRNA gene of the parasite . This method uses the pan-Blastocystis, barcode primers RD5 and BhRDr (Additional file 1: Table S1) to amplify a PCR product of ~ 600 bp. Amplification reactions (25 μl) included 5 μl of template DNA and 0.5 μM of the primer set RD5/BhRDr. Amplification conditions consisted of one step of 95 °C for 3 min, followed by 30 cycles of 1 min each at 94, 59 and 72 °C, with an additional 2 min final extension at 72 °C.
All the direct, semi-nested, and nested PCR protocols described above were conducted on a 2720 thermal cycler (Applied Biosystems). Reaction mixes always included 2.5 units of MyTAQ™ DNA polymerase (Bioline GmbH, Luckenwalde, Germany), and 5× MyTAQ™ Reaction Buffer containing 5 mM dNTPs and 15 mM MgCl2. Laboratory-confirmed positive and negative DNA isolates for each parasitic species investigated were routinely used as controls and included in each round of PCR. PCR amplicons were visualized on 2% D5 agarose gels (Conda, Madrid, Spain) stained with Pronasafe nucleic acid staining solution (Conda). Positive-PCR products were directly sequenced in both directions using the internal primer set described above. DNA sequencing was conducted by capillary electrophoresis using the BigDye® Terminator chemistry (Applied Biosystems) on an on ABI PRISM 3130 automated DNA sequencer.
The Chi-square test was used to compare S. stercoralis, G. duodenalis, Cryptosporidium spp., and Blastocystis spp. infection rates in the surveyed human population according to gender, age group, and place of residence. A probability (P) value < 0.05 was considered evidence of statistical significance. Prevalence risk ratios (PRR) and their 95% confidence intervals (CI) were calculated by univariate analyses to assess the association between potential risk factors considered in the individual data collection spreadsheets and infections with these helminth and protozoan pathogens. Data were analysed with the free software RStudio Version 1.0.44 (https://www.rstudio.com/) using the Epitools library.
Sequence and phylogenetic analyses
Raw sequencing data in both forward and reverse directions were viewed using the Chromas Lite version 2.1 sequence analysis program (https://technelysium.com.au/wp/chromas/). The BLAST tool (http://blast.ncbi.nlm.nih.gov/Blast.cgi) was used to compare nucleotide sequences with sequences retrieved from the NCBI GenBank database. Generated DNA consensus sequences were aligned to appropriate reference sequences using the MEGA 6 software  to identify Giardia species and assemblages/sub-assemblages and Cryptosporidium species. Blastocystis sequences were submitted at the Blastocystis 18S database (http://pubmlst.org/blastocystis/) for sub-type confirmation and allele identification.
For the estimation of the phylogenetic inferences among the identified Giardia-positive samples, a phylogenetic tree was inferred using the Neighbor-Joining method in MEGA 6. Only representative, unambiguous (homozygous) sequences were used in the analyses. The evolutionary distances were computed using the Kimura 2-parameter method, and modelled with a gamma distribution. The reliability of the phylogenetic analyses at each branch node was estimated by the bootstrap method using 1000 replications. Representative sequences of the different G. duodenalis assemblages and sub-assemblages were retrieved from the GenBank database and included in the phylogenetic analysis for reference and comparative purposes. The sequences obtained in this study have been deposited in GenBank under accession numbers MF581531–MF581558 (G. duodenalis), MF581559–MF581564 (Cryptosporidium spp.) and MF581565– MF581576 (Blastocystis spp.).
Full sets of stool samples, epidemiological and clinical data, and signed informed consents were obtained from a total of 351 children. In 24.5% (86/351) of the participants only a single stool sample was collected due to field logistic problems. Additional file 2: Table S2 summarizes the main demographic features of the surveyed population. Briefly, 40.5% (142/351) of the children were recruited in Cubal, 23.9% (84/351) in Tumbulo, 14.0% (49/351) in Capupa, and 21.6% (76/351) in Yambala. The male/female ratio was 0.74. Overall, 28.8% (101/351) of the children fell in the age group 4–7 year-old, 51.3% (180/351) in the age group 8–11 year-old, whereas the remaining 19.9% (70/351) were aged between 12 and 15 years.
Prevalence and risk factors for Strongyloides stercoralis infections
The prevalence of human strongyloidiasis, as determined by qPCR, was estimated at 21.4% [95% confidence interval (CI): 17.1–25.7%]. qPCR-positive results had cycle threshold (Ct) values ranging from 22.7 to 41.9 (median: 32.0). Males and females of all age groups were similarly affected, although those aged between 8 and 15 years harboured higher infection rates than younger children (Table 1). Attending to the frequency of detection of Strongyloides spp. by the commune of origin, the only statistically significant difference found was between children from Cubal and Capupa (P < 0.01). None of the variables considered in the study (gender, age group, walking barefoot, household size, bath place, contact with domestic animals, being moderately to severely underweight, and having anaemia) were associated with an increased risk of having strongyloidiasis (Table 2). Water-related practices including type of drinking water source or washing raw vegetables/fruits were not included in the analyses as these variables are not relevant for the transmission of the disease.
Prevalence and risk factors for enteric protozoan infections
PCR-based results identified G. duodenalis as the most prevalent (37.9%; 95% CI: 32.8–43.0%) enteric protozoan parasite identified in the paediatric population under study, followed by Blastocystis spp. (25.6%; 95% CI: 21.18–30.2%), and Cryptosporidium spp. (2.9%; 95% CI: 1.1–4.5%) (Table 1). Gender was not significantly associated to blastocystosis, cryptosporidiosis, or giardiosis, although the later was found more often in male children than in female children. Individuals in the age group of 4–7 years harboured the greatest bulk of G. duodenalis and Cryptosporidium spp. infections. Indeed, children in that particular age group were found significantly more infected by G. duodenalis than children aged between eight and 11 years (P < 0.01). Blastocystis spp. affected children of all ages. Giardiosis or cryptosporidiosis were not statistically more prevalent in any given commune investigated, although differences in occurrence rates of blastocystosis were observed between children from Cubal and Tumbulo (P < 0.01), Cubal and Yambala (P < 0.05), and Tumbulo and Capupa (P < 0.05) (Table 1). Overall, children younger than seven years of age (PRR: 1.35, P < 0.01) were identified as the sub-population more exposed to infections by protozoan enteropathogens (Table 3). Interestingly, being underweight was found to be a protective factor against these infections (PRR: 0.74, P = 0.005). None of the remaining variables considered in the analysis (gender, household size, water-use practices, contact with domestic animals, and having anaemia) were linked to a higher risk of infection (Table 3).
Single and multiple gastrointestinal helminthic and protozoan infections
The proportions of single or multiple infections involving gastrointestinal helminth or protozoan parasites in the children population are shown in Additional file 3: Table S3. Briefly, 64.1% (225/351) of the participants were infected by at least one of the pathogens investigated. G. duodenalis (30.7%; 69/225) followed by Strongyloides spp. (16.9%; 38/225) and Blastocystis spp. (16.9%; 38/225) were the parasite species most frequently found causing single infections. Similarly, co-infections by G. duodenalis + Blastocystis spp. (15.6%; 35/225) and Strongyloides spp. + G. duodenalis (9.3%; 21/225) were the dual combinations more commonly identified. Additionally, a total of six individuals (3%) harboured triple poliparasitisms by Strongyloides spp. + G. duodenalis + Blastocystis spp., and a single individual (0.4%) was found simultaneously infected by the four pathogen species.
Molecular characterization of G. duodenalis isolates
A total of 133 DNA isolates tested positive for G. duodenalis by qPCR, providing Ct values within the range of 19.0–42.2 (median: 33.2). Of these, 17.3% (23/133) and 9.0% (12/133) were successfully amplified at the gdh and bg loci, respectively. A total of 28 isolates were genotyped and/or sub-genotyped by any of the two markers. Multi-locus genotyping data were available for 25.0% (7/28) of the isolates characterised. The low amplification rates obtained for the gdh and bg markers were found to be highly dependent of qPCR Ct values. Only seven gdh and four bg PCR amplicons were obtained from G. duodenalis isolates with qPCR Ct values > 30, which represented 72.9% (97/133) of the total.
Table 4 shows the diversity, frequency and main features of the G. duodenalis sequences generated at the gdh and bg loci in this survey. Sequence analyses revealed the presence of assemblages A (35.7%, 10/28) and B (64.3%, 18/28). No canine (C, D), feline (F), or ruminant (E) assemblages were detected. Out of the 10 assemblage A sequences, 14.3% (4/28) were assigned to the sub-assemblage AI and 14.3% (4/28) to the sub-assemblage AII. Ambiguous AII/AIII results were obtained for two isolates (7.1%; 2/28) at both the gdh/bg markers. Little genetic diversity was observed among AI-AIII isolates both at the gdh and bg loci, with most sequences (particularly within AII) being identical to their corresponding reference sequences. Single-nucleotide polymorphisms (SNPs) were detected at low rates in a limited number of sequences, some of them associated to amino acid substitutions and representing novel genotypes. Similarly, sequence analyses of the 18 isolates ascribed to the assemblage B allowed the identification of sub-assemblages BIII (7.1%; 2/28) and BIV (25.0%; 7/28). Discordant BIII/BIV typing results were determined in 14.3% (4/28) of the isolates at the gdh marker, whereas the remaining isolates (17.9%; 5/28) were allocated into the assemblage B at the bg (but not the gdh) locus. No inter-assemblage mixed infections were detected. Regardless of the genetic marker considered, virtually all B, BIII and BIV sequences differed among them by 4–8 SNPs, including a number of heterozygous positions in the form of double peaks in the electropherograms. Some of the generated sequences corresponded to genotypes not previously deposited in public repository databases.
Figure 2 shows the phylogenetic relationships among the unambiguous (homozygous) sequences at the gdh locus obtained in the present study and reference sequences from NCBI. Sequences of human origin from developed (Spain) and developing (Ethiopia, Mozambique) countries generated hitherto in our laboratory were also included for comparison [42,43,44,45,46]. Sequences belonging to the assemblages A and B grouped together in well-defined clusters, even at the sub-assemblage (AI/AII or BIII/BIV) level. The exception was a BIV sequence (GenBank: MF581542) that did not comfortably fit within the BIII or BIV clades. Differences in branch length were noticeable in both BIII and BIV sequences, reflecting their comparatively elevated rate of nucleotide substitutions per site. Interestingly, most Spanish BIV sequences tended to form a sub-clade separated from BIV sequences from Angola, Ethiopia and Mozambique, an indication of potential variation in genotype prevalence between different geographical areas.
Molecular characterization of Cryptosporidium spp. isolates
Analyses of SSU rDNA sequences demonstrated that human cryptosporidiosis in Cubal is primarily caused by C. hominis (70.0%; 7/10), with C. parvum (20.0%; 2/10) and C. canis (10.0%; 1/10) being much less frequently detected (Table 5). Among the seven sequences assigned to C. hominis, five showed 100% identity with the reference sequence AF108865, with the remaining two sequences corresponding to novel genotypes of the parasite with a variable number of SNPs including substitutions and the deletion of a single nucleotide. Alignment analyses of the two Cryptosporidium isolates identified as C. parvum allowed the identification of 366–489 bp fragments equivalent to positions 538–913 and 537–1025, respectively, of the reference sequence AF112571. Both sequences also represented genetic variants of C. parvum not previously described. Unfortunately, attempts to amplify the C. hominis and C. parvum isolates at the gp60 locus failed repeatedly, so the sub-genotypes of the parasite involved in the children infections remain unknown. Finally, a single isolate was unambiguously assigned to the Cryptosporidium canine-specific species C. canis. Sequence analyses revealed that this isolate was identical to a stretch of sequence of 492 bp comprising positions 526–1017 of the reference sequence AF112576.
The phylogenetic tree depicted in Fig. 3 illustrates the inferred evolutionary history of the SSU rRNA gene sequences generated in the present study and those gathered from the NCBI database for reference and comparison purposes. As expected, sequences assigned to a specific Cryptosporidium species clustered together in well-defined clades demonstrating the robustness of our analyses.
Molecular characterization of Blastocystis spp. isolates
Out of the 90 isolates that tested positive for Blastocystis spp. by PCR, 83.3% (75/90) were successfully subtyped by sequence analyses at the SSU rDNA (barcode region) gene. BLAST searches allowed identification of five Blastocystis subtypes including ST1 (30.7%; 23/75), ST2 (30.7%; 23/75) and ST3 (36.0%; 27/75). Two additional isolates were recognized as ST5 and ST7 (1.3% each), respectively (Fig. 4). Neither mixed infection involving different STs of the parasite nor infections caused by animal-specific ST10-ST17 were recorded. Allele calling using the Blastocystis SSU database allowed the identification of alleles 4 and 81 within ST1, alleles 12, 15 and 71 within ST2, and alleles 34, 36, 37, 38, 39, 51 and 52 within ST3. Alleles 4 (22.7%; 17/75), 12 (18.7%; 14/75) and 36 (9.3%; 7/75) were found as the most represented in the human population under study. A number of isolates (five in ST1, three in ST2, and one in each ST3, ST5 and ST7) could not be analysed at the allele level due to inaccurate or incomplete sequencing data.
Neglected tropical diseases (NTDs) have been undeniably demonstrated as significant contributors to the global burden of disease, particularly in endemic areas from developing countries . Since 2010 NTDs have received increasing attention and research funding from national and international agencies and institutions including WHO . In contrast, other debilitating rather than fatal parasitic diseases have been often deemed less important in comparison to NTDs associated to higher morbidity/mortality rates. Epidemiological data from recent surveys, including the GBD 2010 initiative, have started to challenge this notion . Thus, a number of helminthic and protozoan infections have emerged as pathogens with a much greater public health impact than initially anticipated. This is the case of S. stercoralis, a nematode not even strictly considered an STH , that have been proven to cause significant morbidity  and preventable mortality [14, 51] in numerous resource-poor tropical and subtropical countries . Consequently, an expert recommendation has been made for the inclusion of S. stercoralis in prevalence studies targeting STH globally . Similarly, the disease burden of cryptosporidiosis among young children has been estimated at almost 48 million disability-adjusted life years (DALYs), a figure comparable to that for tuberculosis (49 million) and more than half of the global burden of malaria (83 million) and HIV/AIDS (82 million) . Although not formally an NTD, cryptosporidiosis (together with giardiosis) was included in the Neglected Diseases Initiative launched by WHO in 2004 . In line with this global trend, molecular epidemiological data presented here clearly indicate that Strongyloides spp., G. duodenalis, Cryptosporidium spp., and (to a still uncertain extent), Blastocystis spp. have a substantial public health impact in Angola.
The overall prevalence of human strongyloidiasis obtained in our study (21.4%) was comparable to that (20.7%) previously documented in Ethiopia using a very similar diagnostic approach , but almost doubled the infection rate (12.8%) recently identified in the same Angolan children population by conventional and Baermann methods . Of note, Strongyloides spp. has been undetected or reported at low (0.3–4%) infection rates in a number of epidemiological surveys conducted in Angola , Cameroon , Ethiopia [56, 57] and Kenya . Taken together, these results highlight the convenience of using specific techniques (e.g. Baermann) and/or highly sensitive methods (e.g. PCR) for the accurate detection of the parasite in endemic areas. Interestingly, a PCR-based prevalence of 48% have been recently found in Mozambican general population , indicating that the prevalence of Strongyloides spp. increases with age, probably as a consequence of autoinfection and long-term (chronic) infection events. Our results seem to corroborate this age-related pattern, as most strongyloidiasis cases were found in children older than seven years.
Reported infection rates of giardiosis (0.1–62%) and cryptosporidiosis (0.1–72%) vary widely among African countries depending on the diagnostic method(s) of choice, the targeted sub-population (primarily children aged 0–16 years), and the immune status of the surveyed individuals . The PCR-based prevalence (38%) for G. duodenalis found in the present study was substantially higher than those previously reported in other Angolan paediatric populations, including children presenting with diarrhoea in the Province of Bengo (21.6%)  and school children in the Province of Huíla (20.1%) . These findings very likely reflect the superior diagnostic sensitivity of PCR over non-molecular methods. In contrast, a 30% Cryptosporidium infection rate was found in the later province using a commercial immunochromatographic rapid test, a 10-fold increase compared with the 2.9% prevalence for cryptosporidiosis identified in the present study by PCR. Variations in the epidemiology and transmission of the parasite or differences in the diagnostic specificity of the detection methods used may account for these discrepancies . Although the epidemiology of blastocystosis in Africa is still poorly understood, a number of PCR-based diagnostic surveys seem to suggest that Blastocystis spp. infections are widespread. Prevalence rates of 61–100% have been reported in Nigeria , Senegal  and Tanzania . Based on the same methodology, a much lower (but still considerable) infection rate of 25.6% was found in the children population under study, this being the first attempt to accurately estimate the prevalence of human blastocystosis in Angola.
A number of epidemiological, clinical, environmental, and behavioural factors have been linked to a higher risk of infections by GHPI. These include walking barefoot [56, 65, 66], having a low socio-economic status [67, 68], living in rural areas [68, 69], having contact with livestock [69,70,71], drinking untreated water [69, 72], eating unwashed/raw fruit , belonging to a given age group [25, 74], having diarrhoea , having a poor nutritional status [74, 75], and having anaemia [8, 73], among others. A number of studies have demonstrated that wearing shoes considerably reduce the risk of infection by STHs [56, 65, 66]. In our study none of the variables potentially involved in the transmission of Strongyloides spp. showed a significant association with the presence of the parasite, although the infection was more prevalent in people walking barefoot, reporting contact with domestic animals, or being part of households with six or more individuals. Similar results have been previously described in this very same population in a previous parasitological survey based on conventional techniques for the detection of the parasite . Clearly, more rigorous epidemiological studies are needed to improve our understanding of the actual risk factors involved in the transmission of human strongyloidiasis. However, significantly different PCR-based prevalence rates of Strongyloides spp. were found between children in the communes of Cubal and Capupa, a finding not previously observed when the detection of the parasite was performed by conventional parasitological techniques . These discrepancies may be attributed to differences in population densities leading to poorer hygienic habits (e.g. defecating in open spaces) or insufficient sanitation facilities (e.g. latrines, sewage treatment) in the most populated communes.
Children aged 4–7 years were found more vulnerable to G. duodenalis and Cryptosporidium spp. infections, with the former pathogen reaching statistical significance. Similar findings have been recurrently documented in other studies worldwide, being associated to the poor personal hygiene habits and immature cellular immunity of individuals in that particular age group. Interestingly, when all three protozoan parasites were collectively considered, infections were significantly less frequent in underweight children. This is in contrast with previous findings highlighting the important role of G. duodenalis, Cryptosporidium spp., and, to some extent, also Blastocystis spp. infections, as leading causes of poor absorption of nutrients, inadequate dietary intake, stunting and its associated cognitive and immunologic sequelae , features commonly seen in children from disfavoured settings in developing countries [5,6,7,8]. At present we do not have a clear explanation for this unexpected outcome. In this regard, iron deficiency linked with malnutrition has been associated with reduced odds of G. duodenalis chronic infection  and other infectious diseases including malaria [78, 79] in highly endemic areas. However, in the present study anaemia was not significantly associated with protection against enteric protozoan parasites. Clearly, more research should be conducted to elucidate this controversial point.
A high proportion (64%) of the children population studied was found infected by at least one enteric pathogen species. This percentage is very similar to that (67%) reported in north-western Angola  but considerably higher than that (44%) reported in the southern part of the country . Overall, PCR-based prevalence data presented here clearly reveal that gastrointestinal helminthic and protozoan infections and co-infections are widespread in Cubal, Benguela. Based on their improved diagnostic performance and shorter turnaround time, our data also support the added benefit of adopting, when possible, molecular methods in high transmission areas, as previously demonstrated by other research groups .
Perhaps the major contribution of this study is the detailed account of the genetic diversity of G. duodenalis, Cryptosporidium spp. and Blastocystis spp. isolates from human origin in Angola, a task that had not been conducted to date. Sequence analyses of G. duodenalis isolates at both the gdh and bg loci revealed exciting data. Assemblage B (64%) had a higher prevalence than assemblage A (36%), proportions well in line with those documented in other African countries [21, 43, 44] and world regions [21, 80]. Confirming the results generated by our research group and others, a high level of genotypic diversity (evidenced by the presence of overlapping nucleotide peaks at specific positions in the corresponding sequencing profiles) was observed within assemblage B and, although to a much lesser extent, also within assemblage A. Consequently, virtually all sequences assigned to sub-assemblages AI, BIII and BIV represented distinct genotypic variants of G. duodenalis. In this regard, we have previously reported a much restricted degree of genotypic diversity in a large clinical cohort of patients with giardiosis in Spain, a developed country with a lower endemicity of the infection. In that survey most of the generated BIII and BIV sequences grouped together in a limited number of G. duodenalis genotypes . We, therefore, envisage an epidemiological scenario for giardiosis in Cubal characterised by high infection pressure and transmission intensity of the parasite. Two independent mechanisms have been proposed to explain the large proportion of heterogeneous sequencing profiles commonly observed in assemblage B sequences: (i) the occurrence of true mixed infections, in which PCR amplification bias lead to the generation of sequences belonging to different assemblages/sub-assemblages of the parasite; or (ii) the occurrence of genetic recombination (e.g. intragenic recombination, exchange of alleles, nuclear fusion within cysts) leading to allelic sequence heterozygosity [80,81,82]. Because no intra-assemblage mixed infections were detected in the present study, we favoured recombination as the driving force behind the genetic variation observed in the assemblage B sequences from Angolan human isolates. Remarkably, four assemblage A sequences were confirmed as sub-assemblage AI at the gdh locus. In Africa, this particular genetic variant has only been identified in a limited number of human isolates from Algeria  and Egypt . Within assemblage A, sub-assemblage AI is known to be responsible for 25% of human infections globally, this percentage being 62%–86% in livestock, 73% in domestic dogs and 44% in wildlife species . Overall, it seems reasonable to think that some (if not all) of the AI infections found in this survey may have a zoonotic origin. Unfortunately, a lack of molecular data from animal sources precluded us from confirming this hypothesis.
Cryptosporidium hominis was confirmed as the most prevalent (70%) Cryptosporidium species causing human infections in Angola, followed by C. parvum (20%). These frequencies were in line with those previously reported by most epidemiological surveys in other African countries . Of interest was the identification of novel genetic variants within both C. hominis and C. parvum at the SSU rRNA locus, a fact that may improve comparative sequence analyses that help understand the population dynamics of Cryptosporidium in the country. Available genotyping data seem to indicate that human cryptosporidiosis in Africa is primarily of anthroponotic origin. This assumption is based on two lines of evidence: (i) the typically anthroponotic C. hominis is the dominant Cryptosporidium species reported in African human populations; and (ii) most human infections with C. parvum are caused by the mainly anthroponotically transmitted IIc subtype of the parasite . Our molecular results seem to support this notion, although the relatively low number of isolates analysed and the failure to characterize them at the gp60 locus made advisable the confirmation of these findings in future studies. In contrast, solid evidence of zoonotic transmission was provided by the identification of the canine-specific C. canis infecting a child. Reflecting the opportunistic nature of these infections, C. canis has been sporadically documented in immunocompromised patients and young children in Ethiopia , Kenya  and Nigeria .
Notably, this survey provides novel data on the molecular epidemiology of Blastocystis spp. in Angola. The large set of genotyping results analysed here confirmed the presence at similar proportions of ST1, ST2 and ST3 in the children population investigated, with ST5 and ST7 causing marginal infections only. Very similar subtype frequencies have been documented in hospital outpatients in Tanzania . In contrast, ST1 was the most prevalent Blastocystis subtype circulating among Nigerian school children  and Libyan outpatients , and ST3 in both asymptomatic and clinical populations in Egypt [88, 89], Senegal  and Tunisia . Noticeably, ST4 was not detected in the investigated paediatric population. This Blastocystis subtype is commonly seen in European countries such as Denmark , France , Sweden  and the UK , but rare or absent in African countries [63, 64, 90, 94], evidencing marked differences in subtype distribution by geographical region and (very likely) between groups of individuals. Although not yet conclusively proven, there is mounting evidence suggesting that Blastocystis pathogenicity could be subtype-related. For instance, ST1 has been linked to the aetiology of irritable bowel syndrome , whereas ST3  and ST4 [91, 97] have been predominantly found in patients with gastrointestinal disorders, mainly diarrhoea. Definitively, more research should be conducted in this field to ascertain the potential correlation between ST and virulence. Regarding transmission, recent phylogenetic analyses have demonstrated that most human ST3 infections fell into specific clades and are the consequence of human-to-human transmission . On the other hand, ST5 is common in pigs  and ST7 in ground-dwelling birds , but considered rare in humans . These findings suggest that the ST5 and ST7 infections detected here may be the result of zoonotic transmission events, although the extent of this possibility should be further confirmed in studies targeting both domestic and wildlife animal species in Angola.
A number of technical limitations may have hampered the accuracy of the diagnostic and genotyping results presented here. For instance, the high Ct values observed in a large proportion of the samples that tested positive to Strongyloides spp. and G. duodenalis by qPCR were somehow unexpected in a theoretically endemic area such as Cubal. This fact could be associated to sub-optimal extraction and purification of genomic DNA from stool samples, or inadequate conservation of DNA isolates prior to PCR testing. If true, reported infection rates may be an underestimation of the actual prevalence of GHPI in this geographical area. This could also explain, at least partially, the relatively low percentage of G. duodenalis isolates genotyped at the sub-assemblage level, or the failure to determine C. hominis/C. parvum subtypes at the gp60 marker.
This is the most thorough study assessing the current epidemiological situation of human strongyloidiasis, giardiosis, cryptosporidiosis, and blastocystosis conducted in Angola to date. Our molecular data clearly demonstrate that the burden of illness associated to these neglected parasitic diseases is higher than initially estimated in the country. Taking into consideration that these GHPI share several common features regarding transmission, surveillance, diagnosis, and treatment with other STHs formally considered as NTDs (case of S. stercoralis) and diarrhoea-causing pathogens (case of G. duodenalis, Cryptosporidium spp. and Blastocystis spp.) it seems wise and advisable to adopt a more holistic and integrated approach to tackle the control and eradication of these diseases jointly, particularly in deprived areas in developing countries. This strategy will undoubtedly contribute to maximize effort and resources devoted to the task.
disability-adjusted life year
global burden of disease
gastrointestinal helminthic and protozoan infections
neglected tropical disease
prevalence risk ratio
real-time polymerase chain reaction
- SSU rRNA:
small subunit ribosomal RNA
World Health Organization
Torgerson PR, Devleesschauwer B, Praet N, Speybroeck N, Willingham AL, Kasuga F, et al. World Health Organization estimates of the global and regional disease burden of 11 food-borne parasitic diseases, 2010: a data synthesis. PLoS Med. 2015 Dec 3;12(12):e1001920.
Houweling TA, Karim-Kos HE, Kulik MC, Stolk WA, Haagsma JA, Lenk EJ, et al. Socioeconomic inequalities in neglected tropical diseases: a systematic review. PLoS Negl Trop Dis. 2016;10:e0004546.
Harhay MO, Horton J, Olliaro PL. Epidemiology and control of human gastrointestinal parasites in children. Expert Rev Anti-Infect Ther. 2010;8:219–34.
Petney TN, Andrews RH. Multiparasite communities in animals and humans: frequency, structure and pathogenic significance. Int J Parasitol. 1998;28:377–93.
Berkman DS, Lescano AG, Gilman RH, Lopez SL, Black MM. Effects of stunting, diarrhoeal disease, and parasitic infection during infancy on cognition in late childhood: a follow-up study. Lancet. 2002;359:564–71.
Ijaz MK, Rubino JR. Impact of infectious diseases on cognitive development in childhood and beyond: potential mitigational role of hygiene. Open Infect Dis J. 2012;6:65–70.
Halliez MC, Buret AG. Extra-intestinal and long term consequences of Giardia duodenalis infections. World J Gastroenterol. 2013;19:8974–85.
Oliveira D, Ferreira FS, Atouguia J, Fortes F, Guerra A, Centeno-Lima S. Infection by intestinal parasites, stunting and anemia in school-aged children from southern Angola. PLoS One. 2015;10:e0137327.
Hotez PJ, Fenwick A, Savioli L, Molyneux DH. Rescuing the bottom billion through control of neglected tropical diseases. Lancet. 2009;373:1570–5.
GBD 2015 Mortality and Causes of Death Collaborators. Global, regional, and national life expectancy, all-cause mortality, and cause-specific mortality for 249 causes of death, 1980–2015: a systematic analysis for the Global Burden of Disease Study 2015. Lancet. 2016;388:1459–544.
Buonfrate D, Mena MA, Angheben A, Requena-Mendez A, Muñoz J, Gobbi F, et al. Prevalence of strongyloidiasis in Latin America: a systematic review of the literature. Epidemiol Infect. 2015;143:452–60.
Schär F, Trostdorf U, Giardina F, Khieu V, Muth S, Marti H, et al. Strongyloides stercoralis: global distribution and risk factors. PLoS Negl Trop Dis. 2013;7:e2288.
Ashford RW, Barnish G, Viney ME. Strongyloides fuelleborni kellyi: infection and disease in Papua New Guinea. Parasitol Today. 1992;8:314–8.
Olsen A, van Lieshout L, Marti H, Polderman T, Polman K, Steinmann P, et al. Strongyloidiasis the most neglected of the neglected tropical diseases? Trans R Soc Trop Med Hyg. 2009;103:967–72.
WHO. The World Health Report. Fighting disease fostering development. Geneva: Switzerland; 1996. http://www.who.int/whr/1996/en/. Accessed 15 Jul 2017.
Parija SC, Jeremiah S. Blastocystis: taxonomy, biology and virulence. Trop Parasitol. 2013;3:17–25.
Roberts T, Stark D, Harkness J, Ellis J. Update on the pathogenic potential and treatment options for Blastocystis sp. Gut Pathog. 2014;6:17.
Ajjampur SS, Tan KS. Pathogenic mechanisms in Blastocystis spp. - interpreting results from in vitro and in vivo studies. Parasitol Int. 2016;65:772–9.
Frealle E, El Safadi D, Cian A, Aubry E, Certad G, Osman M, et al. Acute Blastocystis-associated appendicular peritonitis in a child, Casablanca, Morocco. Emerg Infect Dis. 2015;21:91–4.
Ryan U, Fayer R, Xiao L. Cryptosporidium species in humans and animals: current understanding and research needs. Parasitology. 2014;141:1667–85.
Ryan U, Cacciò SM. Zoonotic potential of Giardia. Int J Parasitol. 2013;43:943–56.
Stensvold CR, Suresh GK, Tan KS, Thompson RC, Traub RJ, Viscogliosi E, et al. Terminology for Blastocystis subtypes - a consensus. Trends Parasitol. 2007;23:93–6.
Clark CG, van der Giezen M, Alfellani MA, Stensvold CR. Recent developments in Blastocystis research. Adv Parasitol. 2013;82:1–32.
Sousa-Figueiredo JC, Gamboa D, Pedro JM, Fançony C, Langa AJ, Magalhães RJ, et al. Epidemiology of malaria, schistosomiasis, geohelminths, anemia and malnutrition in the context of a demographic surveillance system in northern Angola. PLoS One. 2012;7:e33189.
Gasparinho C, Mirante MC, Centeno-Lima S, Istrate C, Mayer AC, Tavira L, et al. Etiology of diarrhea in children younger than 5 years attending the Bengo General Hospital in Angola. Pediatr Infect Dis J. 2016;35:e28–34.
de Alegría MLAR, Colmenares K, Espasa M, Amor A, Lopez I, Ninda A, et al. Prevalence of Strongyloides stercoralis and other intestinal parasite infections in school children in a rural area of Angola: a cross-sectional study. Am J Trop Med Hyg. 2017;97:1226–31.
Perfil do Município do Cubal, Província de Benguela 2009. Administraçao Municipal do Cubal: Ediçoes de Angola Lda; 2009.
Bocanegra C, Gallego S, Mendioroz J, Moreno M, Sulleiro E, Salvador F, et al. Epidemiology of schistosomiasis and usefulness of indirect diagnostic tests in school-age children in Cubal. Central Angola PLoS Negl Trop Dis. 2015;9:e0004055.
MINSA, WFP, WHO, UNICEF. Baseline survey for helminth control in school-aged children in Angola Luanda. Ministerio de Saude; 2005.
de Onis M, Onyango AW, Borghi E, Siyam A, Nishida C, Siekmann J. Development of a WHO growth reference for school-aged children and adolescents. Bull World Health Organ. 2007;85:660–7.
Silva AB, Capingana DP, Magalhães P, Gonçalves MA, Molina MD, Rodrigues SL, et al. Cardiovascular risk factors in pre-pubertal schoolchildren in Angola. Cardiovasc J Afr. 2016;27:315–21.
WHO Haemoglobin concentrations for the diagnosis of anaemia and assessment of severity Vitamin and Mineral Nutrition Information System. Geneva: World Health Organization; 2011. http://www.who.int/vmnis/indicators/haemoglobin.pdf. Accessed 15 Jul 2017.
Saugar JM, Merino FJ, Martín-Rabadán P, Fernández-Soto P, Ortega S, Gárate T, et al. Application of real-time PCR for the detection of Strongyloides spp. in clinical samples in a reference center in Spain. Acta Trop. 2015;42:20–5.
Verweij JJ, Canales M, Polman K, Ziem J, Brienen EA, Polderman AM, et al. Molecular diagnosis of Strongyloides stercoralis in faecal samples using real-time PCR. Trans R Soc Trop Med Hyg. 2009;103:342–6.
Verweij JJ, Schinkel J, Laeijendecker D, van Rooyen MA, van Lieshout L, Polderman AM. Real-time PCR for the detection of Giardia lamblia. Mol Cell Probes. 2003;17:223–5.
Read CM, Monis PT, Thompson RC. Discrimination of all genotypes of Giardia duodenalis at the glutamate dehydrogenase locus using PCR-RFLP. Infect Genet Evol. 2004;4:125–30.
Lalle M, Pozio E, Capelli G, Bruschi F, Crotti D, Cacciò SM. Genetic heterogeneity at the beta-giardin locus among human and animal isolates of Giardia duodenalis and identification of potentially zoonotic subgenotypes. Int J Parasitol. 2005;35:207–13.
Tiangtip R, Jongwutiwes S. Molecular analysis of Cryptosporidium species isolated from HIV-infected patients in Thailand. Tropical Med Int Health. 2002;7:357–64.
Feltus DC, Giddings CW, Schneck BL, Monson T, Warshauer D, McEvoy JM. Evidence supporting zoonotic transmission of Cryptosporidium spp. in Wisconsin. J Clin Microbiol. 2006;44:4303–8.
Scicluna SM, Tawari B, Clark CG. DNA barcoding of Blastocystis. Protist. 2006;157:77–85.
Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: Molecular Evolutionary Genetics Analysis version 6.0. Mol Biol Evol. 2013;30:2725–9.
de Lucio A, Martínez-Ruiz R, Merino FJ, Bailo B, Aguilera M, Fuentes I, et al. Molecular genotyping of Giardia duodenalis isolates from symptomatic individuals attending two major public hospitals in Madrid. Spain PLoS One. 2015;10:e0143981.
de Lucio A, Amor-Aramendía A, Bailo B, Saugar JM, Anegagrie M, Arroyo A, et al. Prevalence and genetic diversity of Giardia duodenalis and Cryptosporidium spp. among school children in a rural area of the Amhara region, north-West Ethiopia. PLoS One. 2016;11:e0159992.
Flecha MJ, Benavides CM, Tissiano G, Tesfamariam A, Cuadros J, de Lucio A, et al. Detection and molecular characterisation of Giardia duodenalis, Cryptosporidium spp. and Entamoeba spp. among patients with gastrointestinal symptoms in Gambo Hospital, Oromia Region, southern Ethiopia. Tropical Med Int Health. 2015;20:1213–22.
Azcona-Gutiérrez JM, de Lucio A, Hernández-de-Mingo M, García-García C, Soria-Blanco LM, Morales L, et al. Molecular diversity and frequency of the diarrheagenic enteric protozoan Giardia duodenalis and Cryptosporidium spp. in a hospital setting in northern Spain. PLoS One. 2017;12:e0178575.
Irisarri-Gutiérrez MJ, Hernández-de Mingo M, de Lucio A, Gil H, Morales L, Seguí R, et al. Association between enteric protozoan parasites and gastrointestinal illness among HIV- and tuberculosis-infected individuals in the Chowke district, southern Mozambique. Acta Trop. 2017;170:197–203.
Hotez PJ, Alvarado M, Basáñez MG, Bolliger I, Bourne R, Boussinesq M, et al. The global burden of disease study 2010: interpretation and implications for the neglected tropical diseases. PLoS Negl Trop Dis. 2014;8:e2865.
Hotez PJ. NTDs V.2.0: "blue marble health" - neglected tropical disease control and elimination in a shifting health policy landscape. PLoS Negl Trop Dis. 2013;7:e2570.
Murray CJ, Vos T, Lozano R, Naghavi M, Flaxman AD, Michaud C, et al. Disability-adjusted life years (DALYs) for 291 diseases and injuries in 21 regions, 1990–2010: a systematic analysis for the Global Burden of Disease Study 2010. Lancet. 2012;380:2197–223.
Pullan RL, Smith JL, Jasrasaria R, Brooker SJ. Global numbers of infection and disease burden of soil-transmitted helminth infections in 2010. Parasit Vectors. 2014;7:37.
Geri G, Rabbat A, Mayaux J, Zafrani L, Chalumeau-Lemoine L, Guidet B, et al. Strongyloides stercoralis hyperinfection syndrome: a case series and a review of the literature. Infection. 2015;43:691–8.
Bisoffi Z, Buonfrate D, Montresor A, Requena-Méndez A, Muñoz J, Krolewiecki AJ, et al. Strongyloides stercoralis: a plea for action. PLoS Negl Trop Dis. 2013;7:e2214.
Savioli L, Smith H, Thompson A. Giardia and Cryptosporidium join the 'Neglected Diseases Initiative'. Trends Parasitol. 2006;22:203–8.
Amor A, Rodriguez E, Saugar JM, Arroyo A, López-Quintana B, Abera B, et al. High prevalence of Strongyloides stercoralis in school-age children in a rural highland of north-western Ehiopia: the role of intensive diagnostic work-up. Parasit Vectors. 2016;9:617.
Campbell SJ, Stothard JR, O'Halloran F, Sankey D, Durant T, Ombede DE, et al. Urogenital schistosomiasis and soil-transmitted helminthiasis (STH) in Cameroon: an epidemiological update at Barombi Mbo and Barombi Kotto crater lakes assessing prospects for intensified control interventions. Infect Dis Poverty. 2017;6:49.
Shiferaw MB, Mengistu AD. Helminthiasis: hookworm infection remains a public health problem in Dera district, South Gondar, Ethiopia. PLoS One. 2015;10:e0144588.
Tadege B, Shimelis T. Infections with Schistosoma mansoni and geohelminths among school children dwelling along the shore of the Lake Hawassa, southern Ethiopia. PLoS One. 2017;12:e0181547.
Easton AV, Oliveira RG, O'Connell EM, Kepha S, Mwandawiro CS, Njenga SM, et al. Multi-parallel qPCR provides increased sensitivity and diagnostic breadth for gastrointestinal parasites of humans: field-based inferences on the impact of mass deworming. Parasit Vectors. 2016;9:38.
Meurs L, Polderman AM, Vinkeles Melchers NV, Brienen EA, Verweij JJ, Groosjohan B, et al. Diagnosing polyparasitism in a high-prevalence setting in Beira, Mozambique: Detection of intestinal parasites in fecal samples by microscopy and real-time PCR. PLoS Negl Trop Dis. 2017;11:e0005310.
Squire SA, Ryan U. Cryptosporidium and Giardia in Africa: current and future challenges. Parasit Vectors. 2017;10:195.
Garcia LS, Shimizu RY, Novak S, Carroll M, Chan F. Commercial assay for detection of Giardia lamblia and Cryptosporidium parvum antigens in human fecal specimens by rapid solid-phase qualitative immunochromatography. J Clin Microbiol. 2003;41:209–12.
Poulsen CS, Efunshile AM, Nelson JA, Stensvold CR. Epidemiological aspects of Blastocystis colonization in children in Ilero, Nigeria. Am J Trop Med Hyg. 2016;95:175–9.
El Safadi D, Gaayeb L, Meloni D, Cian A, Poirier P, Wawrzyniak I, et al. Children of Senegal River Basin show the highest prevalence of Blastocystis sp. ever observed worldwide. BMC Infect Dis. 2014;14:164.
Forsell J, Granlund M, Samuelsson L, Koskiniemi S, Edebro H, Evengård B. High occurrence of Blastsocystis sp. subtypes 1-3 and Giardia intestinalis assemblage B among patients in Zanzibar, Tanzania. Parasit Vectors. 2016;9:370.
Bird C, Ame S, Albonico M, Bickle Q. Do shoes reduce hookworm infection in school-aged children on Pemba Island, Zanzibar? A pragmatic trial. Trans R Soc Trop Med Hyg. 2014;108:297–304.
Paige SB, Friant S, Clech L, Malavé C, Kemigabo C, Obeti R, et al. Combining footwear with public health iconography to prevent soil-transmitted helminth infections. Am J Trop Med Hyg. 2017;96:205–13.
Taye B, Desta K, Ejigu S, Dori GU. The magnitude and risk factors of intestinal parasitic infection in relation to human immunodeficiency virus infection and immune status, at ALERT hospital, Addis Ababa, Ethiopia. Parasitol Int. 2014;63:550–6.
Derso A, Nibret E, Munshea A. Prevalence of intestinal parasitic infections and associated risk factors among pregnant women attending antenatal care center at Felege Hiwot Referral Hospital, northwest Ethiopia. BMC Infect Dis. 2016;16:530.
Helmy YA, Von Samson-Himmelstjerna G, Nöckler K, Zessin KH. Frequencies and spatial distributions of Cryptosporidium in livestock animals and children in the Ismailia province of Egypt. Epidemiol Infect. 2015;143:1208–18.
Wegayehu T, Adamu H, Petros B. Prevalence of Giardia duodenalis and Cryptosporidium species infections among children and cattle in north Shewa Zone, Ethiopia. BMC Infect Dis. 2013;13:419.
Adamu H, Petros B, Zhang G, Kassa H, Amer S, Ye J, et al. Distribution and clinical manifestations of Cryptosporidium species and subtypes in HIV/AIDS patients in Ethiopia. PLoS Negl Trop Dis. 2014;8:e2831.
Gyang VP, Chuang TW, Liao CW, Lee YL, Akinwale OP, Orok A, et al. Intestinal parasitic infections: current status and associated risk factors among school aged children in an archetypal African urban slum in Nigeria. J Microbiol Immunol Infect. 2017; https://doi.org/10.1016/j.jmii.2016.09.005 (In press).
Mengist HM, Taye B, Tsegaye A. Intestinal parasitosis in relation to CD4+T cells levels and anemia among HAART initiated and HAART naive pediatric HIV patients in a model ART center in Addis Ababa. Ethiopia. PLoS One. 2015;10:e0117715.
Tellevik MG, Moyo SJ, Blomberg B, Hjøllo T, Maselle SY, Langeland N, et al. Prevalence of Cryptosporidium parvum/hominis, Entamoeba histolytica and Giardia lamblia among young children with and without diarrhea in Dar es Salaam, Tanzania. PLoS Negl Trop Dis. 2015;9:e0004125.
Nguyen NL, Gelaye B, Aboset N, Kumie A, Williams MA, Berhane Y. Intestinal parasitic infection and nutritional status among school children in Angolela, Ethiopia. J Prev Med Hyg. 2012;53:157–64.
Prendergast AJ, Humphrey JH. The stunting syndrome in developing countries. Paediatr Int Child Health. 2014;34:250–65.
Danquah I, Gahutu JB, Ignatius R, Musemakweri A, Mockenhaupt FP. Reduced prevalence of Giardia duodenalis in iron-deficient Rwandan children. Tropical Med Int Health. 2014;19:563–7.
Oppenheimer SJ. Iron and its relation to immunity and infectious disease. J Nutr. 2001;131:616S–33S.
Sazawal S, Black RE, Ramsan M, Chwaya HM, Stoltzfus RJ, Dutta A, et al. Effects of routine prophylactic supplementation with iron and folic acid on admission to hospital and mortality in preschool children in a high malaria transmission setting: community-based, randomised, placebo-controlled trial. Lancet. 2006;367:133–43.
Sprong H, Cacciò SM, van der Giessen JW. ZOOPNET network and partners. Identification of zoonotic genotypes of Giardia duodenalis. PLoS Negl Trop Dis. 2009;3:e558.
Cooper MA, Adam RD, Worobey M, Sterling CR. Population genetics provides evidence for recombination in Giardia. Curr Biol. 2007;17:1984–8.
Ankarklev J, Jerlström-Hultqvist J, Ringqvist E, Troell K, Svärd SG. Behind the smile: cell biology and disease mechanisms of Giardia species. Nat Rev Microbiol. 2010;8:413–22.
Lalle M, Bruschi F, Castagna B, Campa M, Pozio E, Cacciò SM. High genetic polymorphism among Giardia duodenalis isolates from Sahrawi children. Trans R Soc Trop Med Hyg. 2009;103:834–8.
Helmy MM, Abdel-Fattah HS, Rashed L. Real-time PCR/RFLP assay to detect Giardia intestinalis genotypes in human isolates with diarrhea in Egypt. J Parasitol. 2009;95:1000–4.
Gatei W, Wamae CN, Mbae C, Waruru A, Mulinge E, Waithera T, et al. Cryptosporidiosis: prevalence, genotype analysis, and symptoms associated with infections in children in Kenya. Am J Trop Med Hyg. 2006;75:78–82.
Molloy SF, Smith HV, Kirwan P, Nichols RA, Asaolu SO, Connelly L, et al. Identification of a high diversity of Cryptosporidium species genotypes and subtypes in a pediatric population in Nigeria. Am J Trop Med Hyg. 2010;82:608–13.
Abdulsalam AM, Ithoi I, Al-Mekhlafi HM, Al-Mekhlafi AM, Ahmed A, Surin J. Subtype distribution of Blastocystis isolates in Sebha, Libya. PLoS One. 2013;8:e84372.
Hussein EM, Hussein AM, Eida MM, Atwa MM. Pathophysiological variability of different genotypes of human Blastocystis hominis Egyptian isolates in experimentally infected rats. Parasitol Res. 2008;102:853–60.
Souppart L, Moussa H, Cian A, Sanciu G, Poirier P, El Alaoui H, et al. Subtype analysis of Blastocystis isolates from symptomatic patients in Egypt. Parasitol Res. 2010;106:505–11.
Ben Abda I, Maatoug N, Ben Romdhane R, Bouhelmi N, Zallegua N, Aoun K, et al. Prevalence and subtype identification of Blastocystis sp. in healthy individuals in the Tunis area, Tunisia. Am J Trop Med Hyg. 2017;96:202–4.
Stensvold CR, Christiansen DB, Olsen KEP, Nielsen HV. Blastocystis sp. subtype 4 is common in Danish Blastocystis-positive patients presenting with acute diarrhea. Am J Trop Med Hyg. 2011;14:883–5.
El Safadi D, Cian A, Nourrisson C, Pereira B, Morelle C, Bastien P, et al. Prevalence, risk factors for infection and subtype distribution of the intestinal parasite Blastocystis sp. from a large-scale multi-center study in France. BMC Infect Dis. 2016;16:451.
Forsell J, Granlund M, Stensvold CR, Clark CG, Evengard B. Subtype analysis of Blastocystis isolates in Swedish patients. Eur J Clin Microbiol Infect Dis. 2012;14:1689–96.
Alfellani MA, Stensvold CR, Vidal-Lapiedra A, Onuoha ES, Fagbenro-Beyioku AF, Clark CG. Variable geographic distribution of Blastocystis subtypes and its potential implications. Acta Trop. 2013;126:11–8.
Fouad SA, Basyoni MM, Fahmy RA, Kobaisi MH. The pathogenic role of different Blastocystis hominis genotypes isolated from patients with irritable bowel syndrome. Arab J Gastroenterol. 2011;12:194–200.
Tan TC, Suresh KG, Smith HV. Phenotypic and genotypic characterisation of Blastocystis hominis isolates implicates subtype 3 as a subtype with pathogenic potential. Parasitol Res. 2008;104:85–93.
Domínguez-Márquez MV, Guna R, Muñoz C, Gómez-Muñoz MT, Borrás R. High prevalence of subtype 4 among isolates of Blastocystis hominis from symptomatic patients of a health district of Valencia (Spain). Parasitol Res. 2009;105:949–55.
Stensvold CR, Alfellani M, Clark CG. Levels of genetic diversity vary dramatically between Blastocystis subtypes. Infect Genet Evol. 2012;12:263–73.
Wang W, Owen H, Traub RJ, Cuttell L, Inpankaew T, Bielefeldt-Ohmann H. Molecular epidemiology of Blastocystis in pigs and their in-contact humans in Southeast Queensland, Australia, and Cambodia. Vet Parasitol. 2014;203:264–9.
Noël C, Dufernez F, Gerbad D, Edgcomb VP, Delgado-Viscogliosi P, Ho LC, et al. Molecular phylogenies of Blastocystis isolates from different hosts: implications for genetic diversity, identification of species and zoonosis. J Clin Microbiol. 2005;43:348–55.
The authors are grateful to local authorities from Cubal and school principals of the participating schools for their logistic assistance.
This study was funded by the Health Institute Carlos III (ISCIII), Ministry of Economy and Competitiveness under projects CP12/03081 and PI15CIII/00043 (MPY1148/16). Additional funding was also provided by the ISCIII - Collaborative Research Network on Tropical Diseases (RICET) and the European Regional Development Fund (ERDF) under projects RD16CIII/0003/0004 and RD16/0027/0003.
Availability of data and materials
All relevant data are within the article and its additional files. The sequences data were submitted to the GenBank database under the accession numbers MF581531–MF581576.
Ethics approval and consent to participate
This study has been approved by Ethics Committee of the Vall d’Hebron University Hospital and by the Health and Education regional institutions in Angola. Informative meetings were carried out with village chiefs, school teachers, and community leaders to explain the purpose of the study and the procedures involved. Written informed consent was obtained from parents or legal guardians of all children participating in the survey. Individuals diagnosed with intestinal parasites were referred to general practitioners in order to initiate appropriate treatments.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Oligonucleotides used for the molecular identification and/or characterization of Strongyloides stercoralis, Giardia duodenalis, Cryptosporidium spp. and Blastocystis spp. in this study. (DOCX 14 kb)
Main demographic features of the children population recruited in the present study (n = 351) in Cubal, Angola, 2015. (DOCX 14 kb)
Frequency of single and multiple infections by enteric helminthic and protozoan parasites in Cubal, Angola, 2015. (DOCX 14 kb)
About this article
Cite this article
Dacal, E., Saugar, J.M., de Lucio, A. et al. Prevalence and molecular characterization of Strongyloides stercoralis, Giardia duodenalis, Cryptosporidium spp., and Blastocystis spp. isolates in school children in Cubal, Western Angola. Parasites Vectors 11, 67 (2018). https://doi.org/10.1186/s13071-018-2640-z