First molecular subtyping and phylogeny of Blastocystis sp. isolated from domestic and synanthropic animals (dogs, cats and brown rats) in southern Iran
Parasites & Vectors volume 13, Article number: 365 (2020)
Blastocystis sp. is a common intestinal protist that infects humans and many animals globally. Thus far, 22 subtypes (STs) have been identified in mammalian and avian hosts. Since various STs are common to humans and animals, it was suggested that some human infections might arise from zoonotic transmission. Therefore, the aim of this study was to assess the presence of Blastocystis sp. in domestic (dogs and cats) and synanthropic animals (rats) of Fars Province, Iran, and to genetically characterize the samples.
A total of 400 fresh faecal samples from 154 dogs, 119 cats, and 127 rats were inspected by direct microscopy, Wheatley’s trichrome staining, in vitro culture, and 18S rRNA gene nested-PCR. Finally, sequencing and phylogenetic analyses were performed.
Out of 400 samples, 47 (11.8%) and 61 (15.3%) samples were detected as positive by direct wet mount and culture, respectively. Molecular analysis detected a larger number of positive samples (n = 70, 17.5%): nested-PCR showed that 29 (18.8%) dogs, 21 (17.7%) cats, and 20 (15.8%) rats were infected by Blastocystis sp. Sequence analysis of positive samples indicated the presence of zoonotic STs in all investigated host species. Specifically, ST2 (allele 9), ST3 (allele 34), ST4 (allele 94), ST7 (allele 99), ST8 (allele 21), and ST10 (allele 152) were detected in dogs; ST1 (allele 2), ST3 (allele 34), ST4 (allele 94), ST10 (allele 152), and ST14 (allele 159) were detected in cats; and ST1 (allele 2), ST3 (allele 34), and ST4 (allele 92) were detected in rats.
Our data suggest that domestic dogs and cats can serve as possible reservoirs for in-contact humans, especially those who handle shelter-resident and client-owned animals. Moreover, rats as synanthropic animals can function as a potential source of human infections. Conversely, humans can act as a source of infections to animals. These results should be reinforced in future molecular epidemiological studies.
Blastocystis sp. is an anaerobic eukaryotic protist belonging to the phylum Stramenopiles. It is isolated from the lumen of humans and animals and has four different morphological forms (central-body, granular, amoeboid and cystic) . Binary fission, budding, and plasmotomy are the forms of reproduction observed in Blastocystis sp. . Blastocystis sp. has a global distribution, with a recently reported prevalence of up to 100% in humans . It is believed that Blastocystis sp. is transmitted via the faecal-oral route, water, food, and direct person-to-person contact [1, 4]. There is supporting evidence that zoonotic transmission of Blastocystis sp. from animals to humans living in a community-based environment may occur via close contact with animal enclosures [5,6,7].
In recent years, several studies have reported that Blastocystis sp. is a causative agent of diarrhoea, abdominal pain, flatulence, irritable bowel syndrome, inflammatory bowel disease, and urticaria [8, 9]. However, its pathogenicity is still a topic of some debate, as asymptomatic carriage is common [10, 11]. Considerable evidence suggests that Blastocystis sp. is capable of producing cysteine proteinases that interpose interleukin-8 release from enterocytes, impelling enterocytes apoptosis and augmenting gut permeability while also potentially evading recognition by Toll-like receptors [12, 13].
Blastocystis sp. has been shown to exhibit extensive genetic diversity [14, 15]. Hitherto, 22 different subtypes (STs) of Blastocystis sp., consisting of ST1 to ST17, ST21, and ST23 to ST26 have been described in humans and a variety of animals based on polymorphisms in 18S rRNA gene sequences . Nevertheless, the prevalence of the different subtypes varies among countries and among regions within the same country. ST1 to ST9 and ST12 have been isolated from humans, with ST3 being the most frequent subtype, and all except ST9 have also been reported in animals [17, 18]. ST5 is prevalently isolated from livestock, ST6 and ST7 from birds, and ST8 from non-human primates (NHPs) [19,20,21,22]. The proof that ST5 to ST8 have only been fortuitously discovered in humans has been interpreted as suggestive of zoonotic transmission [19,20,21,22]. ST10 to ST17, ST21, and ST23 to ST26 have been recorded only in non-human hosts thus far [14,15,16].
A substantial number of animal species have been examined for the identification of Blastocystis sp. [5, 14, 20,21,22,23]. However, more studies are required to characterize the zoonotic transmission of this parasite and to identify whether other subtypes of Blastocystis sp. exist. The analysis of DNA extracted directly from stool samples is considered to be highly sensitive, providing the means for genotyping and subtyping [24, 25]. Diagnostic real-time PCR has been recently introduced; however, so far it has not been used for screening .
Blastocystis sp. has been reported to have prevalence rates of 2‒47% in livestock and 8‒67% in captive animals from zoological gardens [14, 20,21,22, 27]. Presently, an expanding interest in Blastocystis sp. research fortifies the necessity for investigations of this parasite in areas where it has not yet been assessed. In Iran, molecular studies were performed in relation to human infection . Concerning animals, Blastocystis sp. was previously genotyped from hooded crows and pigeons , cattle  and wild boars , while there are no molecular studies regarding the occurrence of Blastocystis sp. in animals living in close contact with humans, such as pets and synanthropic animals. Therefore, the aim of this study was to assess the presence of Blastocystis sp. in dogs, cats and rats of Fars Province, southern Iran, and to genetically characterize the samples by 18S rRNA gene sequencing. To the best of our knowledge, this is the first and largest epidemiological study executed on these animals in Iran.
For the present cross-sectional study, 400 fresh faecal samples from 154 dogs (Canis lupus familiaris), 119 cats (Felis catus domesticus), and 127 brown rats (Rattus norvegicus) were collected from December 2016 to October 2018. The study was performed on dogs and cats regardless of their race, age or sex.
Dogs’ fresh faecal samples were collected individually per-rectum under an anaesthesia regime during a trap, neuter, and release (TNR) sterilization programme to restrain stray dog populations controlled by the Fars Veterinary Administration.
After consulting with the vector control unit of the Fars Veterinary Administration, searching and identifying the active colonies of rodents in different places of Fars Province were carried out and wild rodents were trapped. Cats and brown rats were lured in the same spatiotemporal location using Sherman baited cage-traps with tinned fish and roasted walnuts, respectively. Wire cage-traps were placed at rodents’ burrows entrances in gardens and grasslands around the residential houses in different cities and counties of Fars Province. Traps were set in the evening and were checked the next day in the morning. Faeces of captive animals was collected from their cages to minimize contamination using disposable spoons, stored in tagged polystyrene flasks, and transported as early as possible to the head laboratory. Gross macroscopic examination of faecal samples was done to check for consistency. The consistencies of the collected faecal specimens were formed and soft. Every time a mean of 10 traps was used for screening and processing.
Samples were processed within 2‒4 h after collection. Faecal smears were prepared and first examined by direct microscopy using both a saline and iodine wet mount preparation. In addition, faecal smears were stained with Wheatley’s trichrome stain.
In vitro culture was concurrently carried out for all stool samples using Jones’ medium supplemented with 10% heat-inactivated horse serum (Gibco, Frankfurt, Germany), 10% Pen-Strep (1000 IU/ml and 500 μg/ml) and powdered rice-starch. Nearly 200 mg stool specimens were inoculated into 5 ml screw-cap tubes containing 3 ml of Jones’ medium using a sterile disposable applicator stick. The culture tubes were incubated at 37 °C for 72 h . The presence of any of the 4 morphologies of Blastocystis sp. was observed daily by placing one drop of culture product onto SAF-coated coverslips and staining with Wheatley’s trichrome stain. Positive samples were subcultured on LYSGM medium supplemented with 10% heat-inactivated horse serum (Gibco) for mass cultivation .
To lessen impediments from raw fibres and impurities, the faecal specimens were sieved and washed 3 times with distilled water by centrifugation at 1500×g for 10 min. Total genomic DNA of Blastocystis sp. was extracted using a QIAamp® Fast DNA Stool Mini Kit (Qiagen, Hilden, Germany) according to the manufacturer’s instructions. DNA quality was verified by NanoDrop® (Thermo Fisher Scientific, Carlsbad, CA, USA) measurements. The final DNA was eluted in 100 μl of AE buffer to increase its concentration and stored at − 20 °C until use.
Amplification of a 1100-bp fragment of the Blastocystis sp. 18S rRNA gene was carried out using a two-step nested-PCR. This method was chosen since it provides a large fragment for precise sequencing and phylogenetic analysis, and as it is a pan-Blastocystis sp. technique, it permits the identification of known and unknown subtypes of Blastocystis sp. .
The primers RD3 (5′-GGG ATC CTG ATC CTT CCG CAG GTT CAC CTA C-3′); RD5 (5′-GGA AGC TTA TCT GGT TGA TCC TGC CAG TA-3′); BlF (5′-GGA GGT AGT GAC AAT AAA TC-3′); and BlR (5′-CGT TCA TGA TGA ACA ATT AC-3′) were used for the first and second rounds of the nested-PCR. The first PCR reactions were performed in 25 μl reaction mixtures containing 1× PCR buffer (10 mM Tris-HCl, pH 9, 50 mM KCl), 1.5 mM MgCl2, 0.2 mM of each dNTPs (Roche, Alameda, CA, USA), 1.25 U Platinum®Taq High-Fidelity DNA polymerase (Invitrogen, Groningen, the Netherlands), 40 µg/ml of BSA, 10 pM of each forward and reverse primers, and 2 μl genomic DNA.
The PCR conditions consisted of pre-denaturation at 95 °C for 5 min, 40 cycles of denaturation at 95 °C for 30 s, annealing at 65 °C for 45 s, extension at 72 °C for 1 min and a final extension at 72 °C for 10 min using an Eppendorf Mastercycler Gradient PCR machine (Eppendorf, Hamburg, Germany). The second PCR cycle was carried out with similar reaction systems except using the PCR products of the first cycle as the template, and 54 °C as the annealing temperature. Electrophoresis was performed by adding 2 µl of the PCR products and a 100-bp molecular marker (GenScript, Tokyo, Japan) loaded on 1.5% agarose gels (AddGene, Watertown, MA, USA) in TBE buffer and staining with SYBR® Safe DNA Gel Stain (Thermo Fisher Scientific) for 1 h at 90 V. Bands were observed by UV light and photographed (Uvitec, Cambridge, UK). Positive and negative controls were included in every PCR run.
Sequencing and phylogenetic analysis
Bands of the expected size were excised from the agarose gel and purified using a QIAquick® Gel Extraction Kit (Qiagen) according to the manufacturer’s instructions. Purified products were sequenced in both directions using nest-2 primers at Roche Molecular Diagnostics (Roche, Mannheim, Germany) by capillary electrophoresis using a Big DyeTM Terminator Cycle Sequencing Kit (Applied Biosystems, Foster City, CA, USA) in an ABI PRISM® 3730 automated sequencer.
Raw sequencing data in both the forward and reverse directions were checked using Chromas 2.6.6 (Technelysium, Brisbane, Australia). Special attention was paid to the double peaks (indicative of mixed infections by different subtypes or alleles) and the accuracy of the nucleotides was guaranteed. Multiple sequence alignment was performed using the MUSCLE algorithm of MEGA-X . The consensus sequences were then compared with homologous sequences available in the GenBank database using the nBLAST program. The sequences were assembled and edited with BioEdit (v.7.2.6; https://www.bioedit.com). Subtypes were determined by exact match or an identity ≥ 98% against all known Blastocystis sp. subtypes, with a query coverage of ≥ 98%.
All established sequences were submitted to the multilocus sequence typing (MLST) database (https://pubmlst.org/blastocystis/) for subtype confirmation and relevant allele identification.
Two molecular phylogenetic trees (one for dogs and cats, and another for rats) were constructed with the Neighbour-Joining (NJ) method, and genetic distances were calculated with the Maximum Composite Likelihood model in MEGA-X . Bootstrap analysis (with 1000 replicates) was carried out to define the robustness of the findings. Proteromonas lacertae and Blastocystis lapemi were used as outgroups for the two analyses.
The agreement between microscopy, culture, and PCR results was measured using the Cohen’s kappa coefficient (κ) (GraphPad Prism®, Melbourne, Australia) and results were interpreted as follows: no agreement (κ ≤ 0); slight agreement (0.01 < κ < 0.2); fair agreement (0.21 < κ < 0.4); moderate agreement (0.41 < κ < 0.6); substantial agreement (0.61 < κ < 0.8); and almost perfect agreement (0.81 < κ ≤ 1).
Fisher’s exact test was used to compare the frequency of Blastocystis sp. carriage between dogs, cats, and rats (based on molecular results). The significance level was established at P < 0.05.
Infection rates, subtyping and allele analysis
Overall, 400 faecal samples were collected in the current study from 154 dogs, 119 cats, and 127 rats and investigated for Blastocystis sp. Out of 400 faecal samples, 47 (11.8%) samples were shown to be positive by direct wet mount and Wheatley’s trichrome staining examination. The most common finding was the central-body form of the protozoan (Additional file 1: Figure S1). Among 400 stool samples, 61 (15.3%) were found to be positive by in vitro cultivation in Jones’ medium. The most common forms of Blastocystis sp. in Jones’ medium were the central-body and granular forms. κ-values of 0.461, 0.706 and 0.864 were appraised between microscopy, culture, and PCR (Table 1).
Nested-PCR was carried out for amplification of a 1100-bp band (Additional file 2: Figure S2). Overall, 29 faecal samples from dogs (18.8%), 21 faecal samples from cats (17.7%), and 20 faecal samples from rats (15.8%) were detected as positive by nested-PCR (Table 2). No statistically significant differences in prevalence values were recorded among the three species studied (P > 0.05). All positive samples were sequenced successfully. The partial sequences of the 18S rRNA gene of Blastocystis sp. obtained in this study were deposited in the GenBank database under accession numbers MN264509–MN264522.
Among the 70 Blastocystis-positive samples, 8 subtypes were detected: 6 zoonotic STs (ST1–ST4, ST7, and ST8), ST10, and ST14. Dogs were the hosts to the widest range of STs, followed by cats and rats. Remarkably, 6 subtypes were found in the dogs, including ST2 (allele 9), ST3 (allele 34), ST4 (allele 94), ST7 (allele 99), ST8 (allele 21), and ST10 (allele 152). Meanwhile, 5 subtypes were detected in the cats, comprising ST1 (allele 2), ST3 (allele 34), ST4 (allele 94), ST10 (allele 152), and ST14 (allele 159). In addition, 3 subtypes were found in the rats, including ST4 (allele 92), ST3 (allele 34), and ST1 (allele 2) (Table 2). No mixed infections involving different STs of the protist were identified. No intra-subtype genetic heterogeneity was observed within detected subtypes.
A total of 14 representative sequences were obtained from 70 Blastocystis sp. isolates in the present study (Table 2). The sequences obtained in this survey shared high identity with the Blastocystis sp. sequences registered in GenBank.
When the phylogenetic analysis was executed, in contrast to the 17 reference subtype sequences in GenBank, we procured a precise discrimination of Blastocystis sp. subtypes compatible with that established by BLAST queries. Because our study confirmed that ST4 was the most dominant subtype in southern Iranian rodent population, so two phylogenetic trees were constructed to compare this subtype with other rat-derived ST4. All sequences obtained from rats, cats, and dogs were closely related to animal- or human-derived sequences in GenBank and clustered together (Figs. 1 and 2).
Inter-subtype genetic polymorphism was discerned between detected subtypes and reference subtypes registered in GenBank. According to homology analysis, the two ST1 sequences derived from rats and cats (MN264511 and MN264518, respectively) were the same and were 99.36% and 98.91% similar to those from pigs in France (AY135404) and from humans in Japan (AB023499). The only ST2 sequence obtained from dogs (MN264512) shared 99.36% and 98.27% identity with that from pigs in the Philippines (EU445487) and from monkeys in Japan (AB107969). Among the three ST3 sequences, one rat-derived and one cat-derived Blastocystis sp. isolate produced the same sequences (MN264510 and MN264519, respectively), with these sequences being 99.45%, 99.36%, and 99.36% similar to those in cattle in Japan (AB107965), humans in the Philippines (EU445496), and pigs in Germany (MK801366), respectively. The remaining dog-derived ST3 sequence (MN264513) was 99.27% similar to that from dogs in Malaysia (KX234611).
Among the three ST4 isolates, three different sequences were obtained. The rat-derived ST4 sequence (MN264509) was 99.36% and 99.27% similar to that from rats in France (AY135408) and from rats in Japan (AB071000 and MH127499), respectively. The cat-derived ST4 sequence (MN264520) shared 99.36% identity with that from humans in the USA (JN682513) and from rats in Indonesia (MH127488). The remaining dog-derived ST4 sequence (MN264514) was 99.36% similar to that from dogs in Malaysia (KX234618). The only dog-derived ST7 sequence (MN264516) was 99.36% similar to that from chickens in France (AY135410). In addition, the only dog-derived ST8 sequence (MN264515) was 99.36% similar to that from dogs in Malaysia (KX234612).
Among the two ST10 isolates, two different sequences were obtained. The dog-derived ST10 sequence (MN264517) shared 99.36% identity with that from dogs in Malaysia (KX234616). The cat-derived ST10 sequence (MN264521) was 99.36% similar to that from camels in Libya (KC148207) and from goats in Malaysia (KX234638). The only cat-derived ST14 sequence (MN264522) shared 99.27% identity with that from cattle in England (KC148205).
Blastocystis sp. is an enteric protist with a universal distribution. It is assumed that multiple zoonotic isolates of Blastocystis sp. exist with regular animal-to-human and human-to-animal transmission and as a large potential reservoir in animals for infections in humans [4, 35,36,37]. Recent genomic data revealed putative virulence factors and illustrated deleterious clinical outcomes of Blastocystis sp. on the intestinal barrier, leading to conceivable models of pathogenesis [10,11,12,13, 38].
Blastocystis sp. has been secluded from a large spectrum of animals, such as birds, pigs, cattle, sheep, goats, wild carnivores, NHPs, and less frequently rodents, amphibians, reptiles, and insects [21,22,23, 39,40,41]. Blastocystis sp. has cryptic host specificity and is heeded potential zoonotic protozoa since infections in humans have been connected with contact with various mammals and birds [21, 42, 43]. Furthermore, successful experimental infections of rats after oral administration of ST1, ST3, and ST4 purified cysts isolated from human stool samples confirmed the transmission of this protist between human and animal hosts [44, 45]. A higher risk of Blastocystis sp. infection was found in people with close animal contact, supporting the hypothesis of transmission from animals to humans [5, 6, 46, 47].
In a territory where Blastocystis sp. is largely present in humans, there is a necessity to determine the load of this infection in animals and to identify the reservoirs or sources of infection in the community. In Iran, several related studies have been accomplished on the human population , while few surveys have been performed for domestic and wild animals [29,30,31]. To the best of our knowledge, our study represents the first inspection of and report on the prevalence of Blastocystis sp. subtypes isolated from domestic and synanthropic animals (dogs, cats and brown rats) in Fars Province, Iran. Three diagnostic investigations were conducted in parallel on the 400 faecal samples included in the present study: direct microscopy; in vitro culture; and nested-PCR. Studies have shown that PCR is a valid technique with high sensitivity, specificity, reliability and reproducibility . In the present study, nested-PCR allowed us to identify the largest number of positive samples. Previous studies carried out in animals in Iran showed higher infection rates in hooded crows and pigeons (42.9% and 44.4%, respectively)  and in wild boars (25.0%) , although a lower infection rate was reported in cattle (9.6%) .
Globally, epidemiological studies conducted on dogs and cats to investigate their potential role as natural reservoirs of human Blastocystis sp. infection revealed disparate prevalence values. Surveys targeting sheltered and household canine and feline populations did not demonstrate the presence of Blastocystis sp. in Japan and Spain [48, 49], while infection rates in the range of 9.7‒37% have been documented in shelter dogs in the USA , in stray dogs in India , in dogs housed in Italian rescue shelters , and in symptomatic dogs and cats attending a veterinary clinic in Chile .
Concerning brown rats, the recorded prevalence of 15.8% was in line with that reported from previous surveys on rodents, ranging from 8.4% to 45.9% [52, 53]. The large number of positive brown rats, synanthropic animals that share spaces with humans and contaminate the environment with faeces and urine, strengthens their role as potential reservoir hosts for Blastocystis sp.
Furthermore, this study is the first to assess the subtype diversity of Blastocystis sp. based on the 18S rRNA gene in dogs, cats and rats in Fars Province, Iran. Currently, human Blastocystis sp. isolates are classified as ST1-ST9 and ST12, and all except ST9 have also been identified in other animals [17, 18, 25, 39]. In particular, ST1 to ST4 are among the most shared subtypes in humans and have been associated with human disease, having low host specificity and probable zoonotic connotation [54, 55].
In our study, 8 subtypes were identified: ST1–ST4; ST7; ST8; ST10; and ST14. It is noteworthy that the three investigated species harboured not only the same STs, but also the same alleles of each subtype, indicating a possible common environmental source of infection. ST1, detected in the present study in cats and brown rats, was previously reported in captive wild animals, NHPs, sheep, goats, pigs, water voles, marsupials and birds in various locations [21,22,23, 27, 56,57,58,59,60,61]. ST2, recorded in eight dogs, is considered the dominant subtype in NHPs . The zoonotic potential of NHP isolates was previously investigated, in which ST2 isolates with identical sequences were found in rhesus monkeys and children living in the same area of Nepal . Moreover, ST1 and ST2 have been found in zookeepers and five primate species in Australia . ST3 is the Blastocystis sp. subtype with the highest prevalence in humans worldwide and is probably a human species-specific subtype [35, 54, 62]. The finding of ST3 in the three investigated species supported the potential role of domestic and synanthropic animals in human Blastocystis sp. infections.
In the present study, ST1, ST3, and ST4 were detected in brown rats. ST4 was the most prevalent, and it was also found in cats and dogs. Generally, rodents have been shown to be the chief animal reservoir of ST4, which is an ancestry with recent entry into the human population [63, 64]. ST4 has been the most prevalent subtype found in Danish Blastocystis-positive patients presenting with acute diarrhoea . Interestingly, ST4 is reported at low prevalence rates in asymptomatic, apparently healthy individuals in Spain . Most epidemiological studies describing ST4 in human populations depict allele 42, allele 91, and allele 133 of this subtype [49, 65]. However, the ST4 genetic variant reported in rats here is allele 92. This is an interesting point that should be further investigated. In addition, ST4 has been found in NHPs, giraffes, kangaroos, dogs, pigs and ostriches [20, 33, 57, 61, 63]. Moreover, ST3 and ST4 were identified in cockroaches [41, 61, 66].
In the present survey, ST7 and ST8 were also detected in dogs. ST7 is one of the most frequent subtypes in birds and is principally considered an avian subtype [57, 62]. In addition to birds, ST7 has been reported in pigs, goats and dogs [22, 33, 56, 58]. Furthermore, ST8 has been reported more prevalently in NHP handlers, suggesting a zoonotic scatter from primates to their handlers [61, 63]. Recently, ST8 was identified in pigs and marsupials [20, 61].
ST10, detected in the present study both in dogs and cats, has been reported exclusively in non-human hosts, including wild animals and livestock (cattle, yaks, camels, pigs, sheep and goats) in Libya, UK, China, Malaysia, Japan and the UAE [15, 56, 58, 59, 67,68,69,70]. ST10 was also detected in cats from the USA and dogs from France [7, 71].
Finally, we detected ST14 in two cats. ST14 is similar to ST10 in terms of host range, infecting some common livestock and some artiodactyls [57, 68]. Recently, ST14 was isolated from hooded crows and pigeons in Iran .
In dogs, ST1 to ST3 have been identified in Italy , France , Australia , and the Philippines , suggesting that dogs could be involved in the transmission of Blastocystis sp. to humans. Dogs are proficient in shedding possibly zoonotic subtypes and may act as incidental zoonotic reservoirs for infection. Particularly, direct evidence of zoonotic transmission has been granted by other surveys examining human and canine/feline populations living in the same spatiotemporal setting. In a prior study, eight dogs in close contact with 11 symptomatic family members infected by Blastocystis sp. were all positive for the parasite and harboured at least one subtype shared with each of the corresponding patients, suggesting that the source of infection of the owners was their household dogs . Blastocystis ST2 to ST5 were concurrently found in people and domestic dogs living in an urban community in the Philippines . Likewise, ST5 was identified in a villager and a dog both living in a village in Thailand signifying the probability of dogs as sources of zoonotic transmission of Blastocystis sp. . Moreover, all Blastocystis sp. isolated from dogs in India and Indonesia were characterized as zoonotic subtypes [33, 75]. Contrarily, a recent study demonstrate that pet dogs and cats play a tiny role as natural reservoirs of human Blastocystis sp. infection in Spain .
Specifically, stray dogs, which generally live in areas with poor sanitation and hygiene, were identified to be at greater risk of transferring Blastocystis sp. than domestic dogs . Stray dogs in India carried a diverse range of Blastocystis sp. subtypes, including ST1 and ST4 to ST6, while Australian and Cambodian dogs harboured only ST1 and ST2, respectively . Furthermore, ST1 was found in shelter-resident dogs and cats in the Pacific Northwest of the USA, emphasizing the potential role of stray dogs and cats as the origin of human Blastocystis sp. infection .
People who have close contact with dogs and cats for occupational or leisure reasons are at a high risk of acquiring Blastocystis sp. infection. Faced with the fact that coprophagia is a common instinct in dogs, the greater prevalence and diversity of subtypes found in dogs could be ascribed to their increased environmental exposure to faecal content from humans, cattle, sheep, goats and birds. Meanwhile, Blastocystis sp. in animal faeces can enter creeks and rivers by surface run-off after extreme rainfall, which allows water contamination downstream and broad geographical dispersion of Blastocystis sp. The environmental pollution level, size of cysts, resistance of cysts to water treatment methods such as chlorination, and poor elimination of cysts during the filtration process permit the zoonotic transmission of Blastocystis sp. .
The findings of this study showed that dogs, cats, and rats in Fars Province, Iran are shedding Blastocystis sp. Eight subtypes were characterized, with subtype overlaps between animals. Dogs, cats and rats are presumed to be part of the transmission dynamics of this infection in Fars Province. Although the role of these animals as a possible natural reservoir of Blastocystis sp. remains obscure, it seems that these animals could represent possible reservoirs of zoonotic transmission of Blastocystis sp. These data also accentuate the significance of screening other hosts of Blastocystis sp. in order to fully characterize the epidemiology of this protist. Given that extensive genetic diversity exists within Blastocystis sp. subtypes, future molecular characterization and comparisons of dog, cat, rat, human and other mammalian Blastocystis sp. subtypes using multilocus sequence typing (MLST) performed within communities where Blastocystis sp. is endemic in dogs, cats and humans will ideally shed further light on their role as natural hosts for infection. This epidemiological cohort study also furnishes essential information for taking preventive and control measures that should help lessen the risks of zoonotic transmission of Blastocystis sp. to humans.
Availability of data and materials
All data generated or analysed during this study are included in this published article. Sequences were submitted to the GenBank database under the accession numbers MN264509–MN264522.
trap, neuter, and release
sodium acetate-acetic acid-formalin
liver extract-yeast extract-serum-gastric mucin
bovine serum albumin
Tan KSW. New insights on classification, identification, and clinical relevance of Blastocystis sp. Clin Microbiol Rev. 2008;21:639–65.
Zhang X, Zhang S, Qiao J, Wu X, Zhao L, Liu Y, Fan X. Ultrastructural insights into morphology and reproductive mode of Blastocystis hominis. Parasitol Res. 2012;110:1165–72.
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.
Lee LI, Chye TT, Karmacharya BM, Govind SK. Blastocystis sp.: waterborne zoonotic organism, a possibility? Parasit Vectors. 2012;5:130.
Parkar U, Traub RJ, Vitali S, Elliot A, Levecke B, Robertson I, et al. Molecular characterization of Blastocystis isolates from zoo animals and their animal-keepers. Vet Parasitol. 2010;169:8–17.
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.
Ruaux CG, Stang BV. Prevalence of Blastocystis in shelter-resident and client-owned companion animals in the US Pacific Northwest. PLoS ONE. 2014;9:e107496.
Casero RD, Mongi F, Sánchez A, Ramírez JD. Blastocystis and urticaria: examination of subtypes and morphotypes in an unusual clinical manifestation. Acta Trop. 2015;148:156–61.
Cifre S, Gozalbo M, Ortiz V, Soriano JM, Merino JF, Trelis M. Blastocystis subtypes and their association with irritable bowel syndrome. Med Hypotheses. 2018;116:4–9.
Yason JA, Liang YR, Png CW, Zhang Y, Tan KSW. Interactions between a pathogenic Blastocystis subtype and gut microbiota: in vitro and in vivo studies. Microbiome. 2019;7:30.
Lepczyńska M, Białkowska J, Dzika E, Piskorz-Ogórek K, Korycińska J. Blastocystis: how do specific diets and human gut microbiota affect its development and pathogenicity? Eur J Clin Microbiol Infect Dis. 2017;36:1531–40.
Puthia MK, Lu J, Tan KSW. Blastocystis ratti contains cysteine proteases that mediate interleukin-8 response from human intestinal epithelial cells in an NF-κB-dependent manner. Eukaryot Cell. 2008;7:435–43.
Teo JDW, MacAry PA, Tan KSW. Pleiotropic effects of Blastocystis spp. subtypes 4 and 7 on ligand-specific toll-like receptor signaling and NF-κB activation in a human monocyte cell line. PLoS ONE. 2014;9:e89036.
Stensvold CR, Alfellani MA, Clark CG. Levels of genetic diversity vary dramatically between Blastocystis subtypes. Infect Genet Evol. 2012;12:263–73.
Alfellani MA, Taner-Mulla D, Jacob AS, Imeede CA, Yoshikawa H, Stensvold CR, Clark CG. Genetic diversity of Blastocystis in livestock and zoo animals. Protist. 2013;164:497–509.
Stensvold CR, Clark CG. Pre-empting Pandora’s box: Blastocystis subtypes revisited. Trends Parasitol. 2020;36:229–32.
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.
Ramírez JD, Sánchez A, Hernández C, Flórez C, Bernal MC, Giraldo JC, et al. Geographic distribution of human Blastocystis subtypes in South America. Infect Genet Evol. 2016;41:32–5.
Yan Y, Su S, Ye J, Lai X, Lai R, Liao H, et al. Blastocystis sp. subtype 5: a possibly zoonotic genotype. Parasitol Res. 2007;101:1527–32.
Roberts T, Stark D, Harkness J, Ellis J. Subtype distribution of Blastocystis isolates from a variety of animals from New South Wales, Australia. Vet Parasitol. 2013;196:85–9.
Wang J, Gong B, Liu X, Zhao W, Bu T, Zhang W, et al. Distribution and genetic diversity of Blastocystis subtypes in various mammal and bird species in northeastern China. Parasit Vectors. 2018;11:522.
Wang J, Gong B, Yang F, Zhang W, Zheng Y, Liu A. Subtype distribution and genetic characterizations of Blastocystis in pigs, cattle, sheep and goats in northeastern China’s Heilongjiang Province. Infect Genet Evol. 2018;57:171–6.
Calero-Bernal R, Santín M, Maloney JG, Martín-Pérez M, Habela MA, Fernández-García JL, et al. Blastocystis sp. subtype diversity in wild carnivore species from Spain. J Eukaryot Microbiol. 2020;67:273–8.
Roberts T, Barratt J, Harkness J, Ellis J, Stark D. Comparison of microscopy, culture, and conventional polymerase chain reaction for detection of Blastocystis sp. in clinical stool samples. Am J Trop Med Hyg. 2011;84:308–12.
Ramírez JD, Sánchez LV, Bautista DC, Corredor AF, Flórez AC, Stensvold CR. Blastocystis subtypes detected in humans and animals from Colombia. Infect Genet Evol. 2014;22:223–8.
Stensvold CR, Ahmed UN, Andersen LO, Nielsen HV. Development and evaluation of a genus-specific, probe-based, internal-process-controlled real-time PCR assay for sensitive and specific detection of Blastocystis spp. J Clin Microbiol. 2012;50:1847–51.
Betts EL, Gentekaki E, Thomasz A, Breakell V, Carpenter AI, Tsaousis AD. Genetic diversity of Blastocystis in non-primate animals. Parasitology. 2018;145:1228–34.
Badparva E, Ezatpour B, Mahmoudvand H, Behzadifar M, Behzadifar M, Kheirandish F. Prevalence and genotype analysis of Blastocystis hominis in Iran: a systematic review and meta-analysis. Arch Clin Infect Dis. 2017;12:e36648.
Asghari A, Sadraei J, Pirestani M, Mohammadpour I. First molecular identification and subtype distribution of Blastocystis sp. isolated from hooded crows (Corvus cornix) and pigeons (Columba livia) in Tehran Province. Iran. Comp Immunol Microbiol Infect Dis. 2019;62:25–30.
Badparva E, Sadraee J, Kheirandish F. Genetic diversity of Blastocystis isolated from cattle in Khorramabad, Iran. Jundishapur J Microbiol. 2015;8:e14810.
Solaymani-Mohammadi S, Rezaian M, Hooshyar H, Mowlavi GR, Babaei Z, Anwar MA. Intestinal protozoa in wild boars (Sus scrofa) in western Iran. J Wildl Dis. 2004;40:801–3.
Clark CG, Stensvold CR. Blastocystis: isolation, xenic cultivation, and cryopreservation. Curr Protoc Microbiol. 2016;43:20A.1.1–8.
Wang W, Cuttell L, Bielefeldt-Ohmann H, Inpankaew T, Owen H, Traub RJ. Diversity of Blastocystis subtypes in dogs in different geographical settings. Parasit Vectors. 2013;6:215.
Kumar S, Stecher G, Li M, Knyaz C, Tamura K. MEGA X: molecular evolutionary genetics analysis across computing platforms. Mol Biol Evol. 2018;35:1547–9.
Noël C, Dufernez F, Gerbod 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.
Stensvold CR, Nielsen HV, Mølbak K, Smith HV. Pursuing the clinical significance of Blastocystis-diagnostic limitations. Trends Parasitol. 2009;25:23–9.
Stensvold CR, Lewis HC, Hammerum AM, Porsbo LJ, Nielsen SS, Olsen KE, et al. Blastocystis: unravelling potential risk factors and clinical significance of a common but neglected parasite. Epidemiol Infect. 2009;137:1655–63.
Ajjampur SSR, Tan KSW. Pathogenic mechanisms in Blastocystis spp.-interpreting results from in vitro and in vivo studies. Parasitol Int. 2016;65:772–9.
Clark CG, van der Giezen M, Alfellani MA, Stensvold CR. Recent developments in Blastocystis research. Adv Parasitol. 2013;82:1–32.
Alfellani MA, Jacob AS, Perea NO, Krecek RC, Taner-Mulla D, Verweij JJ, et al. Diversity and distribution of Blastocystis sp. subtypes in non-human primates. Parasitology. 2013;140:966–71.
Yoshikawa H, Koyama Y, Tsuchiya E, Takami K. Blastocystis phylogeny among various isolates from humans to insects. Parasitol Int. 2016;65:750–9.
Yoshikawa H, Abe N, Wu Z. PCR-based identification of zoonotic isolates of Blastocystis from mammals and birds. Microbiology. 2004;150:1147–51.
Vaisusuk K, Saijuntha W, Sedlak S, Thanchomnang T, Pilap W, Suksavate W, et al. Blastocystis subtypes detected in long-tailed macaques in Thailand-further evidence of cryptic host specificity. Acta Trop. 2018;184:78–82.
Růžková J, Květoňová D, Jirků M, Lhotská Z, Stensvold CR, Parfrey LW, et al. Evaluating rodent experimental models for studies of Blastocystis ST1. Exp Parasitol. 2018;191:55–61.
Defaye M, Nourrisson C, Baudu E, Warwzyniak I, Bonnin V, Bonnet M, et al. Efficient and reproducible experimental infections of rats with Blastocystis spp. PLoS ONE. 2018;13:e0207669.
Yoshikawa H, Wu Z, Pandey K, Pandey BD, Sherchand JB, Yanagi T, Kanbara H. Molecular characterization of Blastocystis isolates from children and rhesus monkeys in Kathmandu, Nepal. Vet Parasitol. 2009;160:295–300.
Lee IL, Tan TC, Tan PC, Nanthiney DR, Biraj MK, Surendra KM, Suresh KG. Predominance of Blastocystis sp. subtype 4 in rural communities, Nepal. Parasitol Res. 2012;110:1553–62.
Abe N, Nagoshi M, Takami K, Sawano Y, Yoshikawa H. A survey of Blastocystis sp. in livestock, pets, and zoo animals in Japan. Vet Parasitol. 2002;106:203–12.
Paulos S, Köster PC, de Lucio A, Hernández-de-Mingo M, Cardona GA, Fernández-Crespo JC, et al. 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.
Gazzonis AL, Marangi M, Zanzani SA, Villa L, Giangaspero A, Manfredi MT. Molecular epidemiology of Blastocystis sp. in dogs housed in Italian rescue shelters. Parasitol Res. 2019;118:3011–7.
López J, Abarca K, Paredes P, Inzunza E. Intestinal parasites in dogs and cats with gastrointestinal symptoms in Santiago, Chile. Rev Med Chil. 2006;134:193–200.
Chai Y, Deng L, Liu H, Yao J, Zhong Z, Fu H, et al. First subtyping of Blastocystis sp. from pet rodents in southwestern China. Int J Parasitol Parasites Wildl. 2020;11:143–8.
Farah Haziqah MT, Mohd Zain SN, Chandrawathani P, Premaalatha B, Mohd Khairul Nizam MK, Arutchelvan R, Suresh K. Genetic diversity of rodent Blastocystis sp. from Peninsular Malaysia. Trop Biomed. 2018;35:586–92.
Jones MS, Whipps CM, Ganac RD, Hudson NR, Boorom K. Association of Blastocystis subtype 3 and 1 with patients from an Oregon community presenting with chronic gastrointestinal illness. Parasitol Res. 2009;104:341–5.
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;84:883–5.
Song JK, Yin YL, Yuan YJ, Tang H, Ren GJ, Zhang HJ, et al. First genotyping of Blastocystis sp. in dairy, meat, and cashmere goats in northwestern China. Acta Trop. 2017;176:277–82.
Cian A, El Safadi D, Osman M, Moriniere R, Gantois N, Benamrouz-Vanneste S, et al. Molecular epidemiology of Blastocystis sp. in various animal groups from two French zoos and evaluation of potential zoonotic risk. PLoS ONE. 2017;12:e0169659.
Song JK, Hu RS, Fan XC, Wang SS, Zhang HJ, Zhao GH. Molecular characterization of Blastocystis from pigs in Shaanxi province of China. Acta Trop. 2017;173:130–5.
Zhao GH, Hu XF, Liu TL, Hu RS, Yu ZQ, Yang WB, et al. Molecular characterization of Blastocystis sp. in captive wild animals in Qinling Mountains. Parasitol Res. 2017;116:2327–33.
Li WC, Wang K, Gu Y. Occurrence of Blastocystis sp. and Pentatrichomonas hominis in sheep and goats in China. Parasit Vectors. 2018;11:93.
Valença-Barbosa C, do Bomfim TCB, Teixeira BR, Gentile R, Neto SFdC, Magalhães BS, et al. Molecular epidemiology of Blastocystis isolated from animals in the state of Rio de Janeiro, Brazil. PLoS ONE. 2019;14:e0210740.
Greige S, El Safadi D, Bécu N, Gantois N, Pereira B, Chabé M, et al. Prevalence and subtype distribution of Blastocystis sp isolates from poultry in Lebanon and evidence of zoonotic potential. Parasit Vectors. 2018;11:389.
Stensvold CR, Alfellani MA, Nørskov-Lauritsen S, Prip K, Victory EL, Maddox C, et al. Subtype distribution of Blastocystis isolates from synanthropic and zoo animals and identification of a new subtype. Int J Parasitol. 2009;39:473–9.
Katsumata M, Yoshikawa H, Tokoro M, Mizuno T, Nagamoto T, Hendarto J, et al. Molecular phylogeny of Blastocystis isolates from wild rodents captured in Indonesia and Japan. Parasitol Res. 2018;117:2841–6.
Jiménez PA, Jaimes JE, Ramírez JD. A summary of Blastocystis subtypes in North and South America. Parasit Vectors. 2019;12:376.
Yoshikawa H, Wu Z, Howe J, Hashimoto T, Geok-Choo NG, Tan KSW. Ultrastructural and phylogenetic studies on Blastocystis isolates from cockroaches. J Eukaryot Microbiol. 2007;54:33–7.
Tan TC, Tan PC, Sharma R, Sugnaseelan S, Suresh KG. Genetic diversity of caprine Blastocystis from Peninsular Malaysia. Parasitol Res. 2013;112:85–9.
Masuda A, Sumiyoshi T, Ohtaki T, Matsumoto J. Prevalence and molecular subtyping of Blastocystis from dairy cattle in Kanagawa, Japan. Parasitol Int. 2018;67:702–5.
Ren M, Song JK, Yang F, Zou M, Wang PX, Wang D, et al. First genotyping of Blastocystis in yaks from Qinghai Province, northwestern China. Parasit Vectors. 2019;12:171.
AbuOdeh R, Ezzedine S, Madkour M, Stensvold CR, Samie A, Nasrallah G, et al. Molecular subtyping of Blastocystis from diverse animals in the United Arab Emirates. Protist. 2019;170:125679.
Osman M, Bories J, El Safadi D, Poirel MT, Gantois N, Benamrouz-Vanneste S, et al. Prevalence and genetic diversity of the intestinal parasites Blastocystis sp. and Cryptosporidium spp. in household dogs in France and evaluation of zoonotic transmission risk. Vet Parasitol. 2015;214:167–70.
Nagel R, Cuttell L, Stensvold CR, Mills PC, Bielefeldt-Ohmann H, Traub RJ. Blastocystis subtypes in symptomatic and asymptomatic family members and pets and response to therapy. Intern Med J. 2012;42:1187–95.
Belleza ML, Reyes JC, Tongol-Rivera PN, Rivera WL. Subtype analysis of Blastocystis sp. isolates from human and canine hosts in an urban community in the Philippines. Parasitol Int. 2016;65:291–4.
Parkar U, Traub RJ, Kumar S, Mungthin M, Vitali S, Leelayoova S, et al. Direct characterization of Blastocystis from feces by PCR and evidence of zoonotic potential. Parasitology. 2007;134:359–67.
Yoshikawa H, Tokoro M, Nagamoto T, Arayama S, Asih PBS, Rozi IE, Syafruddin D. Molecular survey of Blastocystis sp. from humans and associated animals in an Indonesian community with poor hygiene. Parasitol Int. 2016;65:780–4.
Suresh K, Smith HV, Tan TC. Viable Blastocystis cysts in Scottish and Malaysian sewage samples. Appl Environ Microbiol. 2005;71:5619–20.
The authors would like to thank and appreciate all the veterinary staff throughout Fars Province, Iran, who contributed to the collection of animal faecal samples.
No funding or sponsorship was received either for the study or article publication.
Ethics approval and consent to participate
The research protocol was reviewed and approved by the Institutional Animal Care Committee (IACC) of Shiraz University of Medical Sciences (No. 93‒2817) for animal ethical standards.
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.
The central-body form of Blastocystis sp. stained with Wheatley’s trichrome stain (1000× magnification).
Electrophoresis of PCR products of Blastocystis sp. DNA extracted from dog, cat, and rat faecal samples based on the 18S rRNA gene. The nine lanes contain the products from the negative control (Lane 9), positive control (Lane 8), positive samples of Blastocystis sp. that shows a specific diagnostic band of the 1100-bp fragment (Lanes 1-7), and a molecular size marker (Lane MM).
About this article
Cite this article
Mohammadpour, I., Bozorg-Ghalati, F., Gazzonis, A.L. et al. First molecular subtyping and phylogeny of Blastocystis sp. isolated from domestic and synanthropic animals (dogs, cats and brown rats) in southern Iran. Parasites Vectors 13, 365 (2020). https://doi.org/10.1186/s13071-020-04225-9
- Phylogenetic analysis
- Canis lupus familiaris
- Felis catus domesticus
- Rattus norvegicus