Next Article in Journal
Transmission of Soil Transmitted Helminthiasis in the Mifi Health District (West Region, Cameroon): Low Endemicity but Still Prevailing Risk
Previous Article in Journal
Environmental Correlates of Prevalence of an Intraerythrocytic Apicomplexan Infecting Caribbean Damselfish
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Molecular Detection and Genotyping of Enteric Protists in Asymptomatic Schoolchildren and Their Legal Guardians in Madrid, Spain

by
Pamela Carolina Köster
1,†,
Alejandro Dashti
1,†,
Lucia Reh
1,2,
Begoña Bailo
1,
Marta Hernández-de-Mingo
1,
Aly Salimo Muadica
1 and
David Carmena
1,*
1
Parasitology Reference and Research Laboratory, National Centre for Microbiology, Health Institute Carlos III, Ctra. Majadahonda-Pozuelo Km 2, Majadahonda, 28220 Madrid, Spain
2
European Program for Public Health Microbiology Training, European Centre for Disease Prevention and Control, 169 73 Solna, Sweden
*
Author to whom correspondence should be addressed.
These two authors contributed equally to this work.
Submission received: 26 March 2021 / Revised: 25 April 2021 / Accepted: 26 April 2021 / Published: 2 May 2021

Abstract

:
Asymptomatic carriage of diarrhoea-causing enteric protist parasites in the general population is poorly understood, particularly in medium- to high-income countries. This molecular epidemiological survey investigates the presence, molecular diversity, and household transmission of Giardia duodenalis, Cryptosporidium spp., Blastocystis sp., and Enterocystozoon bieneusi in schoolchildren aged 2–13 years (n = 74) and their legal guardians (n = 6) in Madrid, Spain. Enteroparasite detection and genotyping was conducted in stool samples by molecular (PCR and Sanger sequencing) methods. Potential associations linked to infections were investigated through epidemiological questionnaires. Giardia duodenalis was the most prevalent enteric parasite found (14%, 95% CI: 7.1–23), followed by Blastocystis sp. (10%, 95% CI: 6.2–22) and Cryptosporidium spp. (3.8%, 95% CI: 0.78–11). None of the participants tested positive for E. bieneusi. Sequence analyses revealed the presence of G. duodenalis assemblage B, sub-assemblage BIV in a single child. The three Cryptosporidium isolates obtained were assigned to C. hominis, two of them belonging to the gp60 subtype IbA10G2. Four Blastocystis subtypes were identified including ST2 (38%, 3/8), ST3 (25%, 2/8), ST4 (25%, 2/8), and ST8 (12%, 1/8). All G. duodenalis and Cryptosporidium isolates were detected in children only. Blastocystis ST3 and ST4 were circulating in members of the same household. Blastocystis carriage rates increased with the age of the participants. Presence of diarrhoea-causing enteric protists was common in apparently healthy children.

1. Introduction

Several enteric protistan species contribute to the burden of diarrhoeal illness globally. Those considered of higher public health relevance include the protozoa Cryptosporidium spp., Entamoeba histolytica, Giardia duodenalis, and Cyclospora cayetanensis, and the microsporidia Enterocytozoon bieneusi [1,2,3,4]. Infections by these well-recognized pathogens affect primarily young children and immunocompromised individuals in poor-resource settings characterised by insufficient access to safe drinking water and sanitation [5]. Other species such as the stramenopile Blastocystis sp. and the trichomonad Dientamoeba fragilis have a less-defined pathogenic potential [6,7]. In contrast with the relative abundance of data existing in endemic areas, a comparatively much lower amount of information is currently available on the epidemiology of diarrhoea-causing protistan pathogens in medium- to high-income countries. In these settings, enteric protist infections are erroneously perceived as less relevant threats to public health because of better health standards, lower prevalence rates, less severe symptoms, and rarely fatal consequences [8]. This is in spite that some of these pathogens are frequent cause of food- and waterborne outbreaks of gastrointestinal disease globally [9,10], resulting in significant economic burden associated to medical and treatment costs [11]. Additionally, in developed countries, G. duodenalis and Cryptosporidium spp. infections have been associated with growth retardation and cognitive impairment in children [12] and with long-term complications including perceived food intolerance, irritable bowel syndrome, and chronic fatigue in the general population [13].
Giardia duodenalis, Cryptosporidium spp., Blastocystis sp., and E. bieneusi exhibit extensive intra-species genetic diversity, allowing the identification of several genotypes/subtypes with variable host specificities and ranges. Some of them can infect mammal species other than humans and are therefore zoonotic. Giardia duodenalis is currently regarded as a multi-species complex comprising eight (A to H) distinct assemblages, of which assemblages A and B cause the bulk of human infections. Infection cases by assemblages C-F have been sporadically reported mainly in children and immunocompromised individuals [14]. Cryptosporidium encompasses at least 45 valid species, with C. hominis and C. parvum accounting for more than 90% of documented human cases of cryptosporidiosis [15]. At least 22 subtypes (ST) have been identified within Blastocystis sp., with ST1–4 being responsible for most of the human cases reported globally; infections by ST5–9 and ST12 are only sporadically documented [16]. Finally, several hundred E. bieneusi genotypes have been defined and clustered into 11 phylogenetic groups, with Group 1 and Group 2 including those of zoonotic potential [17].
In the Autonomous Region of Madrid (ARM, Central Spain), the occurrence and genetic diversity of G. duodenalis and Cryptosporidium spp. was initially assessed in children attending day care centres in Majadahonda [18], and in symptomatic individuals of all ages seeking medical attention at two major hospitals in Madrid city [19,20]. A subsequent large-scale molecular epidemiological survey investigated these pathogens, in addition to Blastocystis sp. and E. bieneusi, in asymptomatic schoolchildren in Leganés [21,22]. These studies provided important insights in the understanding of the epidemiology of enteric protists in Spain including (i) G. duodenalis and Blastocystis (but not Cryptosporidium spp. or E. bieneusi) were a common finding in the stools of apparently healthy children, (ii) transmission was mainly of anthropic nature at local foci (e.g., schools), and iii) the occurrence of gastrointestinal symptoms (mainly diarrhoea) was not linked to a given pathogen species/genotype. The present study attempts to complement and expand previous knowledge in a different schoolchildren population in Madrid city. To assess the potential occurrence of within-household transmission events, legal guardians of participating schoolchildren were also investigated.

2. Results

2.1. Occurrence of Protist Enteroparasites

A total of 80 stool samples and associated questionnaires were collected from asymptomatic schoolchildren (n = 74) and their legal guardians (n = 6). Of these, 44 (55%) were male and 34 (43%) were female. The sex of two participants (2.5%) was unknown. The age of the participating schoolchildren ranged from 2 to 13 years with a median age of 7 years. The age of the participating legal guardians ranged from 38 to 55 years with a median age of 44 years. The age of four participants was unknown. The number of participants by age group was as follow: 0–4 years, n = 14; 5–9 years, n = 46; 10–14 years, n = 10; >35 years, n = 6. Of the 80 participants, four (5.0%) were born in countries other than Spain including Bolivia (1), China (1), Cuba (1), and Ecuador (1). This variable was unknown for five (6.3%) participants. Only those PCR amplicons confirmed by Sanger sequencing were considered as positive for any of the enteric protist species investigated here. The full dataset used to determine the prevalence rates and molecular diversity of G. duodenalis, Cryptosporidium spp., Blastocystis sp., and E. bieneusi and the sociodemographic and clinical variables used to estimate risk infections in the present study are presented as a spreadsheet (Supplementary Table S1).
Overall, 20 (25%, 95% CI: 16–36) of the volunteering participants carried at least one enteric protist species; 11 individuals (14%, 95% CI: 7.1–23) were infected with G. duodenalis, three (3.8, 95% CI: 0.78–11) with Cryptosporidium spp., and eight (10%, 95% CI: 6.2–22) with Blastocystis sp. None of the participants tested positive for E. bieneusi. Co-infections by G. duodenalis and Blastocystis sp. were identified in two participants (2.5%).
The distribution of protist enteroparasites according to the sex and age group of the volunteering participants is shown in Figure 1. Giardia duodenalis (18% vs. 8.8%) and Blastocystis sp. (14% vs. 5.9%) were more prevalently found in males than in females, whereas the opposite was true for Cryptosporidium spp. (2.3% vs. 5.9%). Giardia duodenalis was detected in schoolchildren only, particularly in those aged 0–4 years. A similar trend was observed for Cryptosporidium spp. In contrast, Blastocystis sp. followed a marked age-related pattern, with individuals of older age bearing higher carriage rates.

2.2. Molecular Characterization of G. duodenalis Isolates

The 11 DNA isolates that yielded a positive result for G. duodenalis by real-time PCR (qPCR) generated cycle threshold (Ct) values ranging from 18 to 37 (median: 33). Of these, 55% (6/11) had qPCR Ct values >30. Genotyping data were only available for a single isolate (from a nine years-old male) with a Ct = 18 that was successfully amplified at the glutamate dehydrogenase (gdh), ß-giardin (bg), and triose phosphate isomerase (tpi) genes of the parasite. Sequence analyses identified this isolate as assemblage B and sub-assemblage BIV (Table 1).

2.3. Molecular Characterization of Cryptosporidium spp. Isolates

All three Cryptosporidium-positive isolates detected by PCR amplification of the small subunit ribosomal RNA (ssu rRNA) gene were identified as C. hominis, and their partial sequences had 100% identity with reference sequence AF108865. Two of them yielded amplicons of the expected size at the PCR targeting the 60 kDa glycoprotein (gp60) gene of the parasite. Sequence analyses revealed the presence of gp60 subtype IbA10G2 in these samples, in both cases identical to reference sequence AY262031.

2.4. Molecular Characterization of Blastocystis sp. Isolates

Sequence analyses of the eight isolates identified as Blastocystis sp. allowed the identification of four subtypes (ST) including ST2 (38%, 3/8), ST3 (25%, 2/8), ST4 (25%, 2/8), and ST8 (12%, 1/8) (Table 2). Alleles 12 and 11+12 were identified within ST2, allele 34 within ST3, allele 42 within ST4, and allele 95 within ST8. A 55 years-old father and his 9-years old daughter were infected by ST4 and ST3 in the same household, respectively.

2.5. Risk Association Analysis

The results of the univariable analysis of potential risk associations between sociodemographic, epidemiological, clinical, and behavioural factors and a higher likelihood of infection/carriage by G. duodenalis, Cryptosporidium spp., or Blastocystis sp. are shown in Table 3. For Cryptosporidium spp., the only risk association determined was contact with family members with diarrhoea (p = 0.01). Despite the marked age-related pattern observed for Blastocystis sp. in Figure 1c (with individuals of older age bearing higher prevalence rates of the protist), this distribution was not statistically significant (p = 0.1). The same result was obtained when the four original age groups (0‒4, 5‒9, 10‒14, and >35 years) were considered (X-squared = 1.9365, df = 3, p-value = 0.6). No risk associations could be demonstrated for G. duodenalis and Blastocystis sp. infections. Because of the low number of associations found, multivariable analysis was not attempted.

3. Discussion

This molecular-based survey provides novel data on the occurrence and molecular diversity of enteric protist species in Spanish paediatric populations, contributing to improve our knowledge on their epidemiology and transmission. Obtained results complement and expand those generated by our research group in the ARM. The main findings of the study are the confirmation that the protist enteroparasites G. duodenalis, Cryptosporidium spp., and Blastocystis sp. are a common finding in the stools of apparently healthy schoolchildren, and to a lesser extent, in their legal guardians. Evidence supporting the occurrence of within-household transmission was also provided for Blastocystis sp.
Giardia duodenalis was the most prevalent protist enteroparasite found in the present study (14%). Compared with previous PCR-based data available in the ARM, that infection rate was slightly lower than those reported in young children attending day care centres in Majadahonda (16%, 15/90) and in a similar schoolchildren population in Leganés (17%, 263/1,512) [18,22]. In the rest of the country, paediatric infections by G. duodenalis have been identified in the range of 3–36% by light and immunofluorescence microscopy, ELISA for the detection of copro-antigens or rapid diagnostic tests [23,24]. Regarding genotyping data, sequence analyses revealed sub-assemblage BIV in the only G. duodenalis isolate successfully genotyped in the present study. In the ARM, BIV was also the only G. duodenalis sub-assemblage detected in young children attending day care centres [18], and the most prevalent genetic variant of the parasite (79%) in a large school population [22]. The fact that BIV was also the predominant sub-assemblage (62%) in outpatients with clinical manifestations in this very same Spanish region strongly suggests that the genotype of G. duodenalis does not determine per se the outcome of the infection and the occurrence of symptoms [19].
The observed overall prevalence of Cryptosporidium infection (4%) was higher than those (1–3%) reported by PCR-based methods in similar studies conducted in the ARM [18,22]. Infections rates ranging from 1% to 10% have been previously documented in asymptomatic paediatric populations by non-molecular methods in other Spanish regions [24,25,26]. Not surprisingly, C. hominis was the only Cryptosporidium species found in the present study, affecting children only. Indeed, C. hominis has been reported as the most common (65–88%) Cryptosporidium species infecting symptomatic Spanish populations [12,27,28], including those (82%) surveyed in the ARM [20]. In addition, C. hominis caused 71% of the Cryptosporidium infections detected in asymptomatic schoolchildren in the ARM [22]. As in the case of the present survey, most of the above-mentioned C. hominis infections were attributed to the gp60 subtype IbA10G2. Taken together, all these molecular data indicate that the frequency of C. hominis infection is similar in all the Spanish populations investigated to date irrespectively of their clinical status.
In this survey 10% of the investigated children and adults carried Blastocystis sp., a prevalence rate lower than that (13%) identified also by PCR in an equivalent (but much larger) asymptomatic schoolchildren population in the ARM [22]. Both prevalence rates were in the low range of those (8–23%) previously reported by conventional microscopy in other Spanish populations including asymptomatic preschool- and schoolchildren [29], symptomatic outpatients [30], and HIV-positive children [31]. Our sequence analysis confirmed the presence of four Blastocystis STs, being ST2 the most prevalent (38%), followed by ST3 (25%), ST4 (25%), and ST8 (12%). Obtained results were remarkably similar to those (ST2: 36%, ST1: 23%, ST3: 22%, ST4: 19%, ST8: 1%) documented previously in schoolchildren in the ARM [22], with the exception that in the present survey, ST1 was not detected. Interestingly, ST2 (62%) was also the most prevalent Blastocystis ST identified in a community survey conducted in northern Spain [32]. An important contribution of this study is the demonstration that the distribution of Blastocystis sp. followed a clear age-related pattern, with individuals of older age bearing the highest carriage rates. This finding may be the direct consequence of the persistent enteric colonization of the protist, suggesting that Blastocystis sp. should be primarily regarded as a non-pathogenic commensal that forms part of the gut microbiota of the host. Similar age-related patterns for Blastocystis sp. carriage have been observed in healthy children in Spain [21], in orphan children and their caregivers in Thailand [33], and in symptomatic and asymptomatic schoolchildren in Mozambique [34]. In addition, the finding of ST3 and ST4 in a father and his daughter sharing the same household highlights the relevance of asymptomatic carriers in the spreading of faecally–orally microorganisms including Blastocystis sp.
Confirming previous results in apparently healthy schoolchildren populations [22], the microsporidia E. bieneusi was not detected in the surveyed paediatric and adult population. In Spain, E. bieneusi infections have been identified mainly in immunocompromised individuals including HIV/AIDS patients [31,35], transplant recipients [36], and elderly people [37], indicative of the opportunistic nature of this pathogen.
In this study, no animal-adapted species/genetic variants of G. duodenalis and Cryptosporidium spp. were identified, suggesting that the main source of infections by these pathogens should be anthropic in nature, most likely through direct human-to-human contact. However, the finding of Blastocystis ST8 in a 9-year-old male with no obvious risk factors for parasite infections was interesting. This very same Blastocystis ST was previously detected in a 10-year-old female in the ARM [22]. ST8 carriage has been documented primarily in captive and wild non-human primates, but also in swine, several avian species, and rodents [38]. Additionally, zoonotic transmission of ST8 has been documented between captive primates and their handlers in a zoological garden in the United Kingdom [39]. In our study, the source of infection by ST8 was unknown, but a zoonotic transmission event was highly suspected.
The results obtained in the present study and some of the conclusions reached may be biased by some design and methodological limitations. The most important was related to the reduced sample size (particularly those from legal guardians) that have affected the statistical power of the analyses conducted. Low G. duodenalis burdens (evidenced by the elevated qPCR Ct values observed) compromised the success rate of the PCR protocols used for genotyping purposes, as the amplified markers (gdh, bg, and tpi) were all single-copy genes with limited diagnostic sensitivity. In addition, our qPCR protocol lacked amplification internal control. Therefore, it is possible that an unknown number of true Giardia-positive samples were undetected due to the inhibition of the amplification reactions.

4. Materials and Methods

4.1. Study Area and Stool Sample Collection

A cross-sectional molecular epidemiological study of diarrhoea-causing enteric parasites including the protozoan G. duodenalis and Cryptosporidium spp., the stramenopile Blastocystis sp., and the microsporidia E. bieneusi was carried out in voluntary asymptomatic schoolchildren (2–13 years) attending a public school in Madrid (central Spain) in October 2017. An informative meeting was held for interested families and parents/legal guardians were also invited to participate; participating schoolchildren and parents/legal guardians were provided with sampling kits (sterile polystyrene plastic flask with spatula and a unique identification number) to obtain individual stool samples. Parents/legal guardians assisted in the collection of stool samples from consenting schoolchildren and brought the samples to school. Samples were collected by members of the research team at scheduled times, transported to the Spanish National Centre for Microbiology, and stored at 4 °C. All stool samples were processed within 3 days of collection.

4.2. DNA Extraction and Purification

Genomic DNA was isolated from about 200 mg of each faecal specimen using the QIAamp DNA Stool Mini Kit (Qiagen, Hilden, Germany) according to the manufacturer’s instructions, except that samples mixed with InhibitEX buffer were incubated for 10 min at 95 °C. Extracted and purified DNA samples (200 μL) were kept at 4 °C until further molecular analysis.

4.3. Molecular Detection and Characterization of Giardia duodenalis

A qPCR protocol amplifying a 62 bp-fragment of the ssu rRNA gene of Giardia duodenalis was used as initial screening test for the presence of the parasite [40]. 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 (Supplementary Table S2), and 12.5 μL TaqMan® Gene Expression Master Mix (Applied Biosystems, Foster City, 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 with a qPCR-positive result were re-assessed by sequence-based multi-locus genotyping of the single-copy genes gdh, bg, and tpi of the parasite. A semi-nested PCR was used to amplify a 432-bp fragment of the gdh gene [41]. 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 (Supplementary Table S2). 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. A nested PCR was used to amplify a 511 bp-fragment of the bg gene [42]. 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 (Supplementary Table S2). The primary PCR reaction was carried out with the following amplification conditions: one 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. Finally, a nested PCR was used to amplify a 530 bp-fragment of the tpi gene [43]. PCR reaction mixtures (50 μL) included 2–2.5 μL of template DNA and 0.2 μM of the primer pairs AL3543/AL3546 in the primary reaction and AL3544/AL3545 in the secondary reaction (Supplementary Table S2). Both amplification protocols consisted of an initial denaturation step at 94 °C for 5 min, followed by 35 cycles of 94 °C for 45 s, 50 °C for 45 s, and 72 °C for 1 min, with a final extension of 72 °C for 10 min.

4.4. Molecular Detection and Characterization of Cryptosporidium spp.

The presence of Cryptosporidium spp. was assessed using a nested PCR to amplify a 587-bp fragment of the ssu rRNA gene of the parasite [44]. 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 (Supplementary Table S2). 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. Sub-typing of the isolates identified as C. hominis was attempted at the gp60 gene. Briefly, a nested PCR was conducted to amplify an 870-bp fragment of the gp60 locus [45]. PCR reaction mixtures (50 μL) included 2–3 μL of template DNA and 0.3 μM of the primer pairs AL-3531/AL-3535 in the primary reaction and AL-3532/AL-3534 in the secondary reaction (Supplementary Table S2). The primary PCR reaction consisted of an initial denaturation step of 94 °C for 5 min, followed by 35 cycles of 94 °C for 45 s, 59 °C for 45 s, and 72 °C for 1 min with a final extension of 72 °C for 10 min. The conditions for the secondary PCR were identical to the primary PCR except that the annealing temperature was 50 °C.

4.5. Molecular Detection and Characterization of Blastocystis sp.

Identification of Blastocystis sp. was achieved by a direct PCR targeting the ssu rRNA gene of the parasite [46]. This protocol uses the pan-Blastocystis, barcode primers RD5 and BhRDr (Supplementary Table S2) 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.

4.6. Molecular Detection and Characterization of Enterocytozoon bieneusi

Enterocytozoon bieneusi DNA was detected using a nested PCR to amplify a 390 bp-fragment of the entire internal transcribed spacer as well as portions of the flanking large and small subunits of the rRNA gene of the parasite [47]. Reaction mixtures (50 μL) contained one μL of DNA template and 0.2 μM of the primer pairs EBITS3/EBITS4 in the primary reaction and EBITS1/EBITS2.4 in the secondary reaction (Supplementary Table S2). After denaturation at 94 °C for 3 min, amplification conditions of the primary PCR included 35 cycles of amplification (denaturation at 94 °C for 30 s, annealing at 57 °C for 30 s, and elongation at 72 °C for 40 s), followed by a final extension at 72 °C for 10 min. Conditions for the secondary PCR were identical to the primary PCR except only 30 cycles were carried out with an annealing temperature of 55 °C.
All the direct, semi-nested, and nested PCR protocols described above were conducted on a 2720 thermal cycler (Applied Biosystems). Reaction mixes 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 ABI PRISM 3130 automated DNA sequencer.

4.7. Sequence Analyses

Raw sequencing data in both forward and reverse directions were viewed using the Chromas Lite version 2.1 (Technelysium Pty Ltd., South Brisbane, Australia) 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 version 6 software [48] 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. The sequences obtained in this study have been deposited in GenBank under accession numbers MW810324–MW810326 (G. duodenalis), MW789158 and MW810327 (Cryptosporidium spp.), and MW789153‒MW789157 (Blastocystis sp.).

4.8. Epidemiological Questionnaire

A standardised questionnaire was provided as part of the sampling kit to be completed by the children’s legal guardians. Questions included: (i) demographic characteristics, e.g., age, sex, country of birth, and number of siblings; (ii) behavioural habits, e.g., hand and fruit/vegetable washing and whether there have been any occurrence of diarrhoea in the participant, their family members, their schoolmates, and/or pets; and (iii) additional questions on potential risk factors, e.g., types of drinking water, whether they had swum in pools or natural waters in the 2 weeks prior to sample collection, had any contact with pets and any recent travel abroad. Completed questionnaires, signed informed consents, and individual stool samples were returned for collection by each participating schoolchild as described above.

4.9. Statistical Analyses

Data entry from the epidemiological questionnaires was completed using EpiData version 4.2.0 (EpiData Association, Odense, Denmark) and analysed using R software (R Foundation for Statistical Computing, Vienna, Austria). Prevalence of any G. duodenalis, Cryptosporidium spp., and Blastocystis sp. infections (single or multiple) in the study population was calculated. The chi-squared and/or Fisher’s exact test were used to compare G. duodenalis, Cryptosporidium spp., and Blastocystis sp. infection rates by sex and categorical risk factors (including the other protists as a possible risk factor). A p value < 0.05 was considered evidence of statistical significance.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/parasitologia1020010/s1, Table S1: Full dataset used to determine the prevalence rates and molecular diversity of G. duodenalis, Cryptosporidium spp., Blastocystis sp., and E. bieneusi and the sociodemographic and clinical variables used to estimate risk infections in the present study; Table S2: Oligonucleotides used for the molecular identification and/or characterization of the intestinal parasitic and commensal protist species investigated in the present study.

Author Contributions

Conceptualization, P.C.K. and D.C.; methodology, L.R. and D.C.; validation, D.C.; formal analysis, L.R., A.D., P.C.K., and D.C.; investigation, L.R., A.D., P.C.K., B.B., M.H.-d.-M., and A.S.M.; resources, D.C..; data curation, D.C.; writing—original draft preparation, P.C.K., A.D., and D.C.; writing—review and editing, P.C.K., A.D., and D.C.; visualization, P.C.K. and D.C.; supervision, D.C.; project administration, D.C.; funding acquisition, D.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Health Institute Carlos III (ISCIII), Ministry of Science, Innovation and Universities, Spain, grant number PI16CIII/00024.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki and approved by the Ethics Committee of the Health Institute Carlos III on 17 December 2018 under the reference number CEI PI 90_2018-v2.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

All relevant data are within the paper and its Supplementary Materials. The sequences obtained in this study have been deposited in GenBank under accession numbers MW810324–MW810326 (G. duodenalis), MW789158 and MW810327 (Cryptosporidium spp.), and MW789153‒MW789157 (Blastocystis sp.).

Acknowledgments

The authors are grateful to the participating children and their legal guardians, and to the school principal and teachers for their logistic assistance.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Squire, S.A.; Ryan, U. Cryptosporidium and Giardia in Africa: Current and future challenges. Parasites Vectors 2017, 10, 195. [Google Scholar] [CrossRef] [Green Version]
  2. Carrero, J.C.; Reyes-López, M.; Serrano-Luna, J.; Shibayama, M.; Unzueta, J.; León-Sicairos, N.; de la Garza, M. Intestinal amoebiasis: 160 years of its first detection and still remains as a health problem in developing countries. Int. J. Med. Microbiol. 2020, 310, 151358. [Google Scholar] [CrossRef]
  3. Stentiford, G.D.; Becnel, J.; Weiss, L.M.; Keeling, P.J.; Didier, E.S.; Williams, B.P.; Bjornson, S.; Kent, M.L.; Freeman, M.A.; Brown, M.J.F.; et al. Microsporidia-emergent pathogens in the global food chain. Trends Parasitol. 2016, 32, 336–348. [Google Scholar] [CrossRef] [Green Version]
  4. Giangaspero, A.; Gasser, R.B. Human cyclosporiasis. Lancet Infect. Dis. 2019, 19, e226–e236. [Google Scholar] [CrossRef]
  5. Black, R.E.; Morris, S.S.; Bryce, J. Where and why are 10 million children dying every year? Lancet 2003, 361, 2226–2234. [Google Scholar] [CrossRef]
  6. Stensvold, C.R.; Clark, C.G. Current status of Blastocystis: A personal view. Parasitol. Int. 2016, 65, 763–771. [Google Scholar] [CrossRef]
  7. de Boer, M.D.; Schuurs, T.A.; Vermeer, M.; Ruijs, G.J.H.M.; van der Zanden, A.G.M.; Weel, J.F.; Bruijnesteijn van Coppenraet, L.E.S. Distribution and relevance of Dientamoeba fragilis and Blastocystis species in gastroenteritis: Results from a case-control study. Eur. J. Clin. Microbiol. Infect. Dis. 2020, 39, 197–203. [Google Scholar] [CrossRef]
  8. Fletcher, S.M.; Stark, D.; Harkness, J.; Ellis, J. Enteric protozoa in the developed world: A public health perspective. Clin. Microbiol. Rev. 2012, 25, 420–449. [Google Scholar] [CrossRef] [Green Version]
  9. Efstratiou, A.; Ongerth, J.E.; Karanis, P. Waterborne transmission of protozoan parasites: Review of worldwide outbreaks-An update 2011–2016. Water Res. 2017, 114, 14–22. [Google Scholar] [CrossRef]
  10. Zahedi, A.; Ryan, U. Cryptosporidium-An update with an emphasis on foodborne and waterborne transmission. Res. Vet. Sci. 2020, 132, 500–512. [Google Scholar] [CrossRef]
  11. Torgerson, P.R.; Macpherson, C.N.L. The socioeconomic burden of parasitic zoonoses: Global trends. Vet. Parasitol. 2011, 182, 79–95. [Google Scholar] [CrossRef]
  12. Azcona-Gutiérrez, J.M.; de Lucio, A.; Hernández-de-Mingo, M.; García-García, C.; Soria-Blanco, L.M.; Morales, L.; Aguilera, M.; Fuentes, I.; Carmena, D. 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. [Google Scholar] [CrossRef] [Green Version]
  13. Hanevik, K.; Wensaas, K.A.; Rortveit, G.; Eide, G.E.; Mørch, K.; Langeland, N. Irritable bowel syndrome and chronic fatigue 6 years after Giardia infection: A controlled prospective cohort study. Clin. Infect. Dis. 2014, 59, 1394–1400. [Google Scholar] [CrossRef]
  14. Ryan, U.; Cacciò, S.M. Zoonotic potential of Giardia. Int. J. Parasitol. 2013, 43, 943–956. [Google Scholar] [CrossRef]
  15. Yang, X.; Guo, Y.; Xiao, L.; Feng, Y. Molecular epidemiology of human cryptosporidiosis in low- and middle-income countries. Clin. Microbiol. Rev. 2021, 34, e00087–e00119. [Google Scholar] [CrossRef]
  16. Stensvold, C.R.; Clark, C.G. Pre-empting Pandora’s box: Blastocystis subtypes revisited. Trends Parasitol. 2020, 36, 229–232. [Google Scholar] [CrossRef]
  17. Li, W.; Feng, Y.; Santin, M. Host specificity of Enterocytozoon bieneusi and public health implications. Trends Parasitol. 2019, 35, 436–451. [Google Scholar] [CrossRef]
  18. Mateo, M.; Mateo, M.; Montoya, A.; Bailo, B.; Saugar, J.M.; Aguilera, M.; Fuentes, I.; Carmena, D. Detection and molecular characterization of Giardia duodenalis in children attending day care centers in Majadahonda, Madrid, Central Spain. Medicine 2014, 93, e75. [Google Scholar] [CrossRef]
  19. de Lucio, A.; Martínez-Ruiz, R.; Merino, F.J.; Bailo, B.; Aguilera, M.; Fuentes, I.; Carmena, D. Molecular genotyping of Giardia duodenalis isolates from symptomatic individuals attending two major public hospitals in Madrid, Spain. PLoS ONE 2015, 10, e0143981. [Google Scholar] [CrossRef] [Green Version]
  20. de Lucio, A.; Merino, F.J.; Martínez-Ruiz, R.; Bailo, B.; Aguilera, M.; Fuentes, I.; Carmena, D. Molecular genotyping and sub-genotyping of Cryptosporidium spp. isolates from symptomatic individuals attending two major public hospitals in Madrid, Spain. Infect. Genet. Evol. 2016, 37, 49–56. [Google Scholar] [CrossRef]
  21. Reh, L.; Muadica, A.S.; Köster, P.C.; Balasegaram, S.; Verlander, N.Q.; Chércoles, E.R.; Carmena, D. Substantial prevalence of enteroparasites Cryptosporidium spp., Giardia duodenalis and Blastocystis sp. in asymptomatic schoolchildren in Madrid, Spain, November 2017 to June 2018. Eur. Surveill. 2019, 24, 1900241. [Google Scholar] [CrossRef]
  22. Muadica, A.S.; Köster, P.C.; Dashti, A.; Bailo, B.; Hernández-de-Mingo, M.; Reh, L.; Balasegaram, S.; Verlander, N.Q.; Ruiz Chércoles, E.; Carmena, D. Molecular diversity of Giardia duodenalis, Cryptosporidium spp. and Blastocystis sp. in asymptomatic school children in Leganés, Madrid (Spain). Microorganisms 2020, 8, 466. [Google Scholar] [CrossRef] [Green Version]
  23. Carmena, D.; Cardona, G.A.; Sánchez-Serrano, L.P. Current situation of Giardia infection in Spain: Implications for public health. World J. Clin. Infect. Dis. 2012, 2, 1–12. [Google Scholar] [CrossRef]
  24. Cardona, G.A.; Carabin, H.; Goñi, P.; Arriola, L.; Robinson, G.; Fernández-Crespo, J.C.; Clavel, A.; Chalmers, R.M.; Carmena, D. Identification and molecular characterization of Cryptosporidium and Giardia in children and cattle populations from the province of Álava, North of Spain. Sci. Total Environ. 2011, 412–413, 101–108. [Google Scholar] [CrossRef]
  25. de Lucio, A.; Bailo, B.; Aguilera, M.; Cardona, G.A.; Fernández-Crespo, J.C.; Carmena, D. No molecular epidemiological evidence supporting household transmission of zoonotic Giardia duodenalis and Cryptosporidium spp. from pet dogs and cats in the province of Álava, Northern Spain. Acta Trop. 2017, 170, 48–56. [Google Scholar] [CrossRef]
  26. Rodríguez-Hernández, J.; Canut-Blasco, A.; Martín-Sánchez, A.M. Seasonal prevalences of Cryptosporidium and Giardia infections in children attending day care centres in Salamanca (Spain) studied for a period of 15 months. Eur. J. Epidemiol. 1996, 12, 291–295. [Google Scholar] [CrossRef]
  27. Segura, R.; Prim, N.; Montemayor, M.; Valls, M.E.; Muñoz, C. Predominant virulent IbA10G2 subtype of Cryptosporidium hominis in human isolates in Barcelona: A five-year study. PLoS ONE 2015, 10, e0121753. [Google Scholar] [CrossRef] [Green Version]
  28. Abal-Fabeiro, J.L.; Maside, X.; Llovo, J.; Bartolomé, C. Aetiology and epidemiology of human cryptosporidiosis cases in Galicia (NW Spain), 2000–2008. Epidemiol. Infect. 2015, 143, 3022–3035. [Google Scholar] [CrossRef] [Green Version]
  29. Martín-Sánchez, A.M.; Canut-Blasco, A.; Rodríguez-Hernández, J.; Montes-Martínez, I.; García-Rodríguez, J.A. Epidemiology and clinical significance of Blastocystis hominis in different population groups in Salamanca (Spain). Eur. J. Epidemiol. 1992, 8, 553–559. [Google Scholar]
  30. González-Moreno, O.; Domingo, L.; Teixidor, J.; Gracenea, M. Prevalence and associated factors of intestinal parasitisation: A cross-sectional study among outpatients with gastrointestinal symptoms in Catalonia, Spain. Parasitol. Res. 2011, 108, 87–93. [Google Scholar] [CrossRef]
  31. del Águila, C.; Navajas, R.; Gurbindo, D.; Ramos, J.T.; Mellado, M.J.; Fenoy, S.; Muñoz Fernandez, M.A.; Subirats, M.; Ruiz, J.; Pieniazek, N.J. Microsporidiosis in HIV-positive children in Madrid (Spain). J. Eukaryot. Microbiol. 1997, 44, 84S–85S. [Google Scholar] [CrossRef] [PubMed]
  32. Paulos, S.; Köster, P.C.; de Lucio, A.; Hernández-de-Mingo, M.; Cardona, G.A.; Fernández-Crespo, J.C.; Stensvold, C.R.; Carmena, D. Occurrence and subtype distribution of Blastocystis sp. in humans, dogs and cats sharing household in northern Spain and assessment of zoonotic transmission risk. Zoonoses Public Health 2018, 65, 993–1002. [Google Scholar] [CrossRef]
  33. Pipatsatitpong, D.; Rangsin, R.; Leelayoova, S.; Naaglor, T.; Mungthin, M. Incidence and risk factors of Blastocystis infection in an orphanage in Bangkok, Thailand. Parasites Vectors 2012, 5, 37. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Muadica, A.S.; Balasegaram, S.; Beebeejaun, K.; Köster, P.C.; Bailo, B.; Hernández-de-Mingo, M.; Dashti, A.; Dacal, E.; Saugar, J.M.; Fuentes, I.; et al. Risk associations for intestinal parasites in symptomatic and asymptomatic schoolchildren in central Mozambique. Clin. Microbiol. Infect. 2021, in press. [Google Scholar] [CrossRef] [PubMed]
  35. del Águila, C.; López-Velez, R.; Fenoy, S.; Turrientes, C.; Cobo, J.; Navajas, R.; Visvesvara, G.S.; Croppo, G.P.; Da Silva, A.J.; Pieniazek, N.J. Identification of Enterocytozoon bieneusi spores in respiratory samples from an AIDS patient with a 2-year history of intestinal microsporidiosis. J. Clin. Microbiol. 1997, 35, 1862–1866. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Galván, A.L.; Sánchez, A.M.; Valentín, M.A.; Henriques-Gil, N.; Izquierdo, F.; Fenoy, S.; del Aguila, C. First cases of microsporidiosis in transplant recipients in Spain and review of the literature. J. Clin. Microbiol. 2011, 49, 1301–1306. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Lores, B.; López-Miragaya, I.; Arias, C.; Fenoy, S.; Torres, J.; del Aguila, C. Intestinal microsporidiosis due to Enterocytozoon bieneusi in elderly human immunodeficiency virus–negative patients from Vigo, Spain. Clin. Infect. Dis. 2002, 34, 918–921. [Google Scholar] [CrossRef] [Green Version]
  38. Hublin, J.S.Y.; Maloney, J.G.; Santin, M. Blastocystis in domesticated and wild mammals and birds. Res. Vet. Sci. 2021. [Google Scholar] [CrossRef]
  39. Stensvold, C.R.; Alfellani, M.A.; Nørskov-Lauritsen, S.; Prip, K.; Victory, E.L.; Maddox, C.; Nielsen, H.V.; Clark, C.G. Subtype distribution of Blastocystis isolates from synanthropic and zoo animals and identification of a new subtype. Int. J. Parasitol. 2009, 39, 473–479. [Google Scholar] [CrossRef] [Green Version]
  40. Verweij, J.J.; Schinkel, J.; Laeijendecker, D.; van Rooyen, M.A.; van Lieshout, L.; Polderman, A.M. Real-time PCR for the detection of Giardia lamblia. Mol. Cell. Probes 2003, 17, 223–225. [Google Scholar] [CrossRef]
  41. Read, C.M.; Monis, P.T.; Thompson, R.C. Discrimination of all genotypes of Giardia duodenalis at the glutamate dehydrogenase locus using PCR-RFLP. Infect. Genet. Evol. 2004, 4, 125–130. [Google Scholar] [CrossRef] [PubMed]
  42. Lalle, M.; Pozio, E.; Capelli, G.; Bruschi, F.; Crotti, D.; Cacciò, S.M. 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–213. [Google Scholar] [CrossRef]
  43. Sulaiman, I.M.; Fayer, R.; Bern, C.; Gilman, R.H.; Trout, J.M.; Schantz, P.M.; Das, P.; Lal, A.A.; Xiao, L. Triose phosphate isomerase gene characterization and potential zoonotic transmission of Giardia duodenalis. Emerg. Infect. Dis. 2003, 9, 1444–1452. [Google Scholar] [CrossRef] [PubMed]
  44. Tiangtip, R.; Jongwutiwes, S. Molecular analysis of Cryptosporidium species isolated from HIV-infected patients in Thailand. Trop. Med. Int. Health 2002, 7, 357–364. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Feltus, D.C.; Giddings, C.W.; Schneck, B.L.; Monson, T.; Warshauer, D.; McEvoy, J.M. Evidence supporting zoonotic transmission of Cryptosporidium spp. in Wisconsin. J. Clin. Microbiol. 2006, 44, 4303–4308. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Scicluna, S.M.; Tawari, B.; Clark, C.G. DNA barcoding of Blastocystis. Protist 2006, 157, 77–85. [Google Scholar] [CrossRef]
  47. Buckholt, M.A.; Lee, J.H.; Tzipori, S. Prevalence of Enterocytozoon bieneusi in swine: An 18-month survey at a slaughterhouse in Massachusetts. Appl. Environ. Microbiol. 2002, 68, 2595–2599. [Google Scholar] [CrossRef] [Green Version]
  48. Tamura, K.; Stecher, G.; Peterson, D.; Filipski, A.; Kumar, S. MEGA6: Molecular evolutionary genetics analysis version 6.0. Mol. Biol. Evol. 2013, 30, 2725–2729. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Prevalence rates of the protist enteroparasites investigated in the present study according to the sex and age group of the volunteering participants: (a) Giardia duodenalis, (b) Cryptosporidium spp., and (c) Blastocystis sp.
Figure 1. Prevalence rates of the protist enteroparasites investigated in the present study according to the sex and age group of the volunteering participants: (a) Giardia duodenalis, (b) Cryptosporidium spp., and (c) Blastocystis sp.
Parasitologia 01 00010 g001
Table 1. Diversity, frequency, and molecular features of Giardia duodenalis sequences at the glutamate dehydrogenase (gdh), ß-giardin (bg), and triose phosphate isomerase (tpi) loci generated in the present study. GenBank accession numbers are provided.
Table 1. Diversity, frequency, and molecular features of Giardia duodenalis sequences at the glutamate dehydrogenase (gdh), ß-giardin (bg), and triose phosphate isomerase (tpi) loci generated in the present study. GenBank accession numbers are provided.
LocusAssemblageSub-AssemblageIsolatesReference SequenceStretchSingle Nucleotide PolymorphismsGenBank ID
gdhBBIV1L4050876–491NoneMW810324
bgB1AY072727100–593G159A, C165T, C309T, C324T, C393T, T471CMW810325
tpiBBIV1AF0695601–479NoneMW810326
Table 2. Diversity, frequency, and molecular features of Blastocystis sp. sequences at the small subunit ribosomal RNA (ssu rRNA) locus generated in the present study. GenBank accession numbers are provided.
Table 2. Diversity, frequency, and molecular features of Blastocystis sp. sequences at the small subunit ribosomal RNA (ssu rRNA) locus generated in the present study. GenBank accession numbers are provided.
SubtypeAlleleIsolatesReference SequenceStretchSingle Nucleotide PolymorphismsGenBank ID
ST2122AB0709978–599G185A, T267A, C324T, A477GMW789153
11+121AB07099719–607G185R, T267A, C324T, A477C, T478CMW789154
ST3342AB10796513–602NoneMW789155
ST4422AB07100014–606G78A, G221T, A234GMW789156
ST8951AB10797055–614C128A, A177C, G179A, G190C, T228C, A256TMW789157
Table 3. Univariable analysis. Crude association between Giardia duodenalis, Cryptosporidium spp., and Blastocystis sp. infections and variables of interest. P-values marked in bold indicate numbers that are significant on the 95% confidence limit (CI).
Table 3. Univariable analysis. Crude association between Giardia duodenalis, Cryptosporidium spp., and Blastocystis sp. infections and variables of interest. P-values marked in bold indicate numbers that are significant on the 95% confidence limit (CI).
VariableCategorynGiardia duodenalisCryptosporidium spp.Blastocystis sp.
CasesOR95% CIpCasesOR95% CIpCasesOR95% CIp
GenderMale448 1 6
Female3432.30.56–9.40.3320.370.0–4.30.6022.50.48–130.50
Age group≤7 years407 2 1
>7 years3641.70.45–6.30.4311.80.16–21>0.9950.160.02–1.40.10
Contact with infantsExposed90 0 1
Unexposed61100.0NA0.3330.00.0–NA>0.9951.40.14–140.60
Diarrhoea in familyExposed162 3 1
Unexposed6080.920.09–5.4>0.990Inf1.7–Inf0.0160.600.07–5.4>0.99
Diarrhoea in schoolmatesExposed145 0 0
Unexposed30Inf0.16–Inf0.510NANA>0.9910.00.0–8.40.20
Recent travel abroadExposed252 1 1
(EU countries)Unexposed5180.470.09–2.40.5021.50.13–18>0.9960.310.04–2.70.40
Recent travel abroadExposed81 0 1
(non-EU countries)Unexposed6890.930.10–8.5>0.9930.0NA>0.9961.50.15–140.60
Contact with dogsExposed50 0 0
Unexposed72100.0NA0.4030.0NA>0.9980.00.0–NA>0.99
Contact with catsExposed80 0 2
Unexposed69100.0NA0.5930.0NA>0.9963.50.57–21.30.20
Source of drinking waterTap7711 3 8
Bottled2521.90.40–9.30.5010.970.10–9.8>0.9930.850.21–3.5>0.99
Springs610.830.09–7.8>0.990Inf0.03–Inf>0.9920.230.04–1.50.20
Others410.500.05–5.30.500Inf0.02–Inf>0.9910.350.03–3.80.40
Swimming in poolsExposed332 2 2
Unexposed4170.310.06–1.60.3012.60.22–300.6060.380.07–2.00.30
HandwashingAlways192 1 2
Habitually4951.00.18–5.9>0.9921.30.11–15>0.9960.840.15–4.6>0.99
Never930.240.03–1.80.300Inf0.01–Inf>0.990Inf0.09–Inf0.30
Vegetable washingAlways609 3 5
Habitually1712.80.33–240.440Inf0.10–Inf>0.9930.420.09–1.90.40
Never000.0NA>0.9900.0NA>0.9900.0NA>0.99
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Köster, P.C.; Dashti, A.; Reh, L.; Bailo, B.; Hernández-de-Mingo, M.; Muadica, A.S.; Carmena, D. Molecular Detection and Genotyping of Enteric Protists in Asymptomatic Schoolchildren and Their Legal Guardians in Madrid, Spain. Parasitologia 2021, 1, 83-94. https://0-doi-org.brum.beds.ac.uk/10.3390/parasitologia1020010

AMA Style

Köster PC, Dashti A, Reh L, Bailo B, Hernández-de-Mingo M, Muadica AS, Carmena D. Molecular Detection and Genotyping of Enteric Protists in Asymptomatic Schoolchildren and Their Legal Guardians in Madrid, Spain. Parasitologia. 2021; 1(2):83-94. https://0-doi-org.brum.beds.ac.uk/10.3390/parasitologia1020010

Chicago/Turabian Style

Köster, Pamela Carolina, Alejandro Dashti, Lucia Reh, Begoña Bailo, Marta Hernández-de-Mingo, Aly Salimo Muadica, and David Carmena. 2021. "Molecular Detection and Genotyping of Enteric Protists in Asymptomatic Schoolchildren and Their Legal Guardians in Madrid, Spain" Parasitologia 1, no. 2: 83-94. https://0-doi-org.brum.beds.ac.uk/10.3390/parasitologia1020010

Article Metrics

Back to TopTop