An Overview on the Main Assemblages and Sub-assemblages of Giardia intestinalis in the Western Half of Iran

authors:

avatar Zahra Eslamirad ORCID 1 , avatar Reza Hajihossein ORCID 1 , avatar Azam Moslemi ORCID 2 , *

Department of Parasitology and Mycology, School of Medicine, Arak University of Medical Sciences, Arak, Iran
Department of Biostatistics, School of Medicine, Arak University of Medical Sciences, Arak, Iran

how to cite: Eslamirad Z, Hajihossein R, Moslemi A. An Overview on the Main Assemblages and Sub-assemblages of Giardia intestinalis in the Western Half of Iran. Shiraz E-Med J. 2022;23(10):e116536. https://doi.org/10.5812/semj-116536.

Abstract

Context:

Studies have shown that various Assemblages, sub-assemblage (subtypes) and genotypes of Giardia intestinalis have multiple hosts; therefore, their distribution reservoirs as well as pattern of epidemiological distribution are different. The present study was conducted to summarize the results of genetic studies on G. intestinalis in Iran.

Evidence Acquisition:

To investigate the articles on the sub-assemblage of G. intestinalis in Iran, a systematic search was conducted in Persian and English databases. The search process led to the entry of 23 articles into this systematic review. Also, in this study, to estimate the ratio of the dominant sub-assemblage of Giardia parasite in Iran, meta-analysis was used and a significant level of 0.05 was considered.

Results:

The results of this study showed that investigations on assemblages and sub-assemblage of G. intestinalis had been mostly made in the western half of Iran. The overall prevalence of A, B, and the mixed assemblages of Giardia intestinalis in these areas of Iran was estimated to be 0.56, 0.27, and 0.16, respectively. Also, the prevalence of AI and AII sub-assemblage of G. intestinalis in the same areas of Iran was 0.34 and 0.5, respectively. The prevalence of BIII and BIV sub-assemblage of this parasite in mentioned areas was further found to be 0.2 and 0.06, respectively. In addition, the total prevalence of all sub-assemblage (AI, AII, BIII, & BIV) was calculated to be 0.18 (P < 0.001).

Conclusions:

Given that assemblage A of G. intestinalis has the highest prevalence in the human society of western half of Iran, in this area of the country, the main epidemiological pattern of Giardia transmission is zoonotic. However, to identify the source of this parasite spread, more studies are needed.

1. Context

Giardiasis is a global disease caused by a flagellate protozoan called Giardia. The only pathogenic species of this parasite in humans is Giardia intestinalis (also known as G. duodenalis and G. lamblia). Although this disease is often asymptomatic, the most common symptom is diarrhea. The causative agent of the disease is resistant cysts, which are excreted from the host intestine and transmitted by oral-fecal, water, and food (1, 2). Thus, identifying the parasite's reservoir is very important to prevent it.

Studies have shown different hosts and as a result different reservoir for various assemblages of G. intestinalis. Based on its protein and DNA polymorphism, eight major genetic groups or assemblages (A, B, C, D, E, F, G, and H) have been identified for this parasite (3-5). Numerous studies have been also performed to determine the sub-assemblages of this parasite around the world (6). Most of human isolates are in assemblages A and B (3, 4, 6). Because these two assemblages could also infect other mammals, the zoonotic transmission potential of G. intestinalis is epidemiologically important (3). In some parts of Iran, studies have been conducted to identify the sub-assemblages of G. intestinalis in humans. So, it is valuable to summarize these studies to determine the predominant sub-assemblages (subtypes) and epidemiological transmission pattern of this parasite in Iran. On the other hand, such studies could identify information gaps regarding the sub-assemblage of this parasite, thereby opening new avenues for further research. Hence, the present study aimed to provide an overview on the determination of genotypes of G. intestinalis and summarize the studies on this area in Iran.

2. Evidence Acquisition

2.1. Search Process

All English or Persian articles and dissertations performed on Giardia genotype in Iran (from 1980 to 2020) indexed in PubMed, Google scholar, Science Direct, Scopus, Medline, Medlib, Scientific Information Database (SID), IranMedex, IranDoc, and Magiran were collected and reviewed. The keywords were a combination of Giardia, intestinal protozoa, assemblage, genotype, G. intestinalis, G. duodenalis, G. lamblia, Giardia genotype, Iran, Islamic Republic of Iran, genotype of Giardia in Iran. Figure 1 shows the search process.

Steps to select and approve of articles for systematic review and meta-analysis
Steps to select and approve of articles for systematic review and meta-analysis

2.2. Selection criteria

The following inclusion criteria were considered for this study. Full papers and dissertations related to determining Giardia genotypes based on different genes of this parasite, genetic characters’ and human Giardia genotype in Iran published from 1980 to 2020 were included in this work. In contrast, case studies, experimental studies, and animal studies, as well as duplicate articles were excluded.

2.3. Data Extraction

Selected papers were carefully reviewed by two researchers and information, such as first author, year of publication, type of study, location of study, language, subjects, sample size, diagnosis test, assemblage or genotype detection technique, and genotype were extracted and recorded.

2.4. Statistical Analysis

In this study, meta-analysis was used to estimate the proportion of dominant genotypes of Giardia parasite in Iran (7). Confidence intervals were also calculated using the exact binomial distribution. Further, publication bias was investigated using the Egger test (8). Heterogeneity among the studies was considered by benefiting from the Cochran’s Q test with a significant level lower than 0.1 and I2 statistic greater than 60% (9). Due to the observed heterogeneity, we used random effects model to estimate the prevalence of predominant genotype of Giardia parasite in Iran.

The meta-analysis of proportions was further performed by using “Metaprop” commands on STATA software (Ver.14). The significance level was considered as 0.05 (7). Freeman-Tukey double arcsine transformation was used in the random effects model when the proportion in some studies was very large and close to 100% (7). Additionally, a Forest plot was drawn based on the random effects model for the estimated proportion of assemblages with confidence interval of 95%.

3. Results

After searching the databases, 173 articles were found related to the Giardia exploration in Iran. After excluding the unrelated and duplicated studies, a total number of 23 associated articles to the sub-assemblages of human G. intestinalis in Iran were included. Figure 1 shows the process of designing and searching for articles, and Figure 2 illustrates the geographical areas of the included articles (Figures 1 & 2).

The geographical areas related to the reviewed articles (Provinces in which has been studied the Giardia intestinalis genotype)
The geographical areas related to the reviewed articles (Provinces in which has been studied the Giardia intestinalis genotype)

From 23 studies, one was cohort study and the rest were cross-sectional. These studies had been conducted in 13 provinces of Iran. Four studies in Khuzestan, three studies in Fars, three studies in Isfahan, two studies in East Azerbaijan, two studies in Tehran, two studies in Kerman, and one study in other provinces (7 provinces) (Figure 2). Totally, 1059 human specimens were examined in all the 23 included studies. The PCR-RFLP technique had been used for detecting of Giardia assemblages in 11 studies, while in others, the sequencing technique had been used (Table 1).

Table 1.

Extraction of Information About the Assemblage of Giardia intestinalis from the Selected Articles a

LineFirst AuthorType of StudyPlace of Study (City)Characteristic of the Sample Under StudyNumber of Samples Molecularly TestedDiagnosis MethodDetermination Method AssemblageAssemblageDescription and Conclusion
ABMixNeg
1Babaei et al. (10)Cross sectionalTehranStool samples were reported to be positive microscopically38PCR (Based on the gdh gene)RFLP33 (87)3 (7.8)2 (5.2)-Epidemiologic transmission pattern was anthroponotic
2Fallah et al. (11)Cross sectionalTabrizLaboratory samples of Tabriz Children's Hospital laboratory and reference laboratory samples that were positive microscopic method34PCR (Based on the tim gene)Seq17 (50)13 (38)1 (3)3
3Fallah et al. (12)Cross sectionalTabrizTabriz Children's Hospital and reference laboratory samples34PCR (Based on the gdh gene)RFLP6 (17.6)12 (35.3)-16Epidemiologic transmission pattern was anthroponotic
4Etemadi et al. (13)Cross sectionalKermanKerman Medical Centers30PCR (Based on the gdh gene)RFLP23 (76.6)7 (23.4)--
5Pestehchian et al. (14)Cross sectionalIsfahanIsfahan Health Centers67PCR (Based on the gdh gene)RFLP40 (59.7)25 (37.3)2 (3)-Epidemiologic transmission pattern was anthroponotic and zoonotic.
6Sarkari et al. (15)Cross sectionalFars provinceFars Province Laboratories205PCR (Based on the gdh gene)RFLP128 (62.4)36 (17.6)8 (3.9)33Epidemiologic transmission pattern was anthroponotic
7Manouchehri Naeini et al. (16)Cross sectionalShahreKordShahrekord Medical Diagnosis Laboratories31PCR (Based on the tpi gene)Seq11 (35.5)16 (51.6)4 (12.9)-Association was seen between diarrhea and genotype A. And association was seen between asymptomatic infection and genotype B
8Rafiei et al. (17)Cross sectionalAhwazAhwaz Health Centers100PCR (Based on the gdh gene)RFLP18 (18)28 (28)54 (54)-There is a association between clinical symptoms and genotype A
9Roointan et al. (18)Sectional (September 2011 to July 2012)AhwazChildren referred to health center clinics50PCR (Based on the gdh gene)RFLP5 (10)8 (16)37 (74)-
10Etemadi et al. (19)Cross sectionalKermanThree medical diagnostic laboratories30PCR (Based on the gdh gene)RFLP23 (76.6)7 (23.4)--
11Pestechian et al. (20)Cross sectionalIsfahan provinceMedical diagnostic laboratories67PCR (Based on the tpi gene)RFLP40 (59.7)25 (37.3)2 (3)-
12Rayani et al. (21)Cross sectionalShirazHealth centers and hospitals in Shiraz50PCR (Based on the SSU-rDNA)Seq32 (80)8 (20)-10Epidemiologic transmission pattern was anthroponotic
13Rayani et al. (22)Cross sectionalShirazHealth centers and hospitals in Shiraz40PCR (Based on the gdh gene)Seq32 (80)8 (20)--Epidemiologic transmission pattern was anthroponotic
14Bahrami et al. (23)Cross sectionalKurdistan province14 Medical Diagnosis Laboratories23Nested-PCR (Based on the tpi gene)Seq12 (52)11 (48)--Epidemiologic transmission pattern was zoonotic
15Nasiri Goorabi et al. (24)Cross sectionalBaharestan cityHealth centers of Baharestan city in Tehran province25PCR (Based on the β-giardin gene)RFLP16 (64)9 (36)--Epidemiologic transmission pattern was anthroponotic and anthropozoonotic
16Hooshyar et al. (25)Cross-sectionalKashanMedical Diagnosis Laboratories47PCR (Based on the gdh gene)RFLP24 (51.1)9 (19.1)11 (23.4)3Epidemiologic transmission pattern anthropozoonotic
17Effati et al. (26)Cross-sectionalAlborz ProvinceShahid Bahonar Children's Hospital of Karaj And Medical diagnostic laboratories Fardis, Nazar abad and Eshtehard15Nested-PCR (Based on the tpi gene)Seq2 (13.3)5 (33.3)-8 (53.4)Epidemiologic transmission pattern anthropozoonotic and zoonotic
18Kashinahanji et al. (27)Cross sectionalHamedanClients of health centers23Nested-PCR (Based on the tpi gene)Seq18 (78.2)5 (21.7)--Epidemiologic transmission pattern zoonotic
19Shahnazi et al. (28)Cross sectional)QazvinPeople who come in contact with food20PCR (Based on the gdh gene)Seq11 (55)6 (30)-3
20Mirrezaie et al. (29)Cross sectionalAndimeshkImam Ali Hospital and Health Center in Andimeshk40PCR (Based on the gene gdh and β-giardin)Seq23 (57.5)17 (42.5)--Only the number of samples that were sequencing was listed, but 84 samples were tested for PCR.
21Mahmoudi et al. (30)Cross-sectional aRashtVisitors to 3 educational hospitals of Gilan University of Medical Sciences41PCR (Based on the gdh gene)Seq38 (92.6)3 (7.4)--60 samples were PCR, but 41 samples were sent for sequencing.
22Rafiei et al. (31)Cohort studyShushtarReferrals to health centers24PCR (Based on genes SSU-rRNA )Seq12 (50)12 (50)--Based on SSU-rRNA gene 90 samples was PCR test, which confirmed 82 samples, but PCR based on gdh, Tpi, and bg genes was performed on 24 samples
23Abdi et al. (32)Cross sectionalArakGiardia-infected samples from Arak University of Medical Sciences Medical Diagnosis Laboratory it was gathered25PCR (Based on the betagiardin gene)Seq20 (80)---Assemblage of 5 samples was unknown

The results of statistical analysis to estimate the prevalence of A, B, as well as the mixed assemblages of G. intestinalis in the Iranian population are as follows.

Results of the Egger test revealed that in 23 studies, the proportion of assemblage A lacked publication bias (P = 0.72), the proportion of assemblage B, had publication bias (P = 0.086), while the proportion of both assemblage (mix) was deficient in publication bias (P = 0.71). On the other hand, results of the Cochran's Q test and I2 statistic showed a heterogeneity among studies for assemblages A (I2 = 91.06%, Q = 246.02, P-value<0.001), B (I2 = 71.37%, Q = 73.35, P-value < 0.001), and mixture of both assemblages (I2 = 96.34%, Q = 218.55 I2, P-value < 0.001). So, we applied the random effects model in this meta-analysis. Figure 3A & B) displays the forest plot for the estimated proportions of assemblage A and B with confidence interval of 95% in western half of Iran. As shown in this figure, the overall proportion of assemblage A in this region of country is 0.56 with a confidence interval of 95% (0.46, 0.67) (P < 0.001), whilst the overall proportion of assemblage B has been estimated at 0.27 in the same region with a confidence interval of 95% (0.22, 0.33) (P < 0.001). Moreover, Figure 3C presents the forest plot for the mixed proportion of both A and B assemblages. This figure shows that the overall mixed proportion of both cases is 0.16 in mentioned area with a confidence interval of 95% (0.04, 0.35) (P < 0.001).

Forest plot based on the random effects model for assemblages of Giardia intestinalis in the western half of Iran. (A) assemblage A, (B) assemblage B (C) the mixed A and B assemblages
Forest plot based on the random effects model for assemblages of Giardia intestinalis in the western half of Iran. (A) assemblage A, (B) assemblage B (C) the mixed A and B assemblages

As noticed in Table 2, out of 23 selected studies, sub-assemblages of G. intestinalis have been presented in 16 studies.

Table 2.

Extraction of Information About Sub-assemblages of Giardia intestinalis from the Reviewed Articles a

LineFirst AuthorPlace of StudySample SizeAssemblage (Genotype or Sub Assemblage)
ABMixNeg
AIAIIBIIIBIVB Novel
1Babaei et al. (10)Tehran3833 (87)-3 (7.8)--2 (5.2)-
2Fallah et al. (11)Tabriz34-6 (33.3)8 (44.4)4 (22.2)--16
3Etemadi et al. (13)Kerman305 (16.6)18 (60)7 (23.4)----
4Pestehchian et al. (14)Esfahan67-40 (59.7)23 (34.32)2 (2.98)-2 (2.98)-
5Sarkari et al. (15)Fars province205-128 (74.41)30 (17.44)6 (3.49)-8 (4.66)33
6Rafiei et al. (17)Ahwaz100-18 (18)28 (28)--54 (54)-
7Roointan et al. (18)Ahwaz50-5 (10)8 (16)--37 (74)-
8Etemadi et al. (19)Kerman305 (16.6)18 (60)7 (23.4)----
9Rayani et al. (21)Shiraz50-32 (80)4 (10)4 (10)--10
10Rayani et al. (22)Shiraz40-32 (80)4 (10)4 (10)---
11Bahrami et al. (23)Kurdistan234 (17.4)8 (34.8)5 (21.7)4 (17.4)2 (8.7)--
12Hooshyar et al. (25)Kashan47-24 (51.1)7 (14.9)2 (4.2)-11 (23.4)3
13Effati et al. (26)Alborz4-1 (25)2 (50)1 (25)---
14Shahnazi et al. (28)Qazvin2011 (55)6 (30)---3
15Mahmoudi et al. (30)Rasht41-38 (92.6)-3 (7.4)---
16Rafiei et al. (31)Shushtar24-12 (50)12 (50)--

Results of statistical analysis are presented below to estimate the proportion of sub-assemblages A, B, and mixed of Giardia intestinalis in the population of the western half of Iran.

Furthermore, results of the Egger test indicated that the proportion of Giardia intestinalis sub-assemblages A (AI, AII) in the mentioned area lacked publication bias (P > 0.05). On the other hand, the Cochran's Q test and I2statistic showed a heterogeneity among the studies for AI sub-assemblages (I2 = 97.13%, Q = 209.02, P-value <0.001), and AII sub-assemblages (I2 = 92.44%, Q = 185.11, P-value <0.001). So, we applied the random effects model in this meta-analysis. As seen in Figure 4A-D, the overall proportion of AI (Figure 4A) was estimated 0.34 with a confidence interval of 95% (0.04, 0.74) (P < 0.05). In addition, the overall proportion of AII (Figure 4B) was estimated to be 0.5 with a confidence interval of 95% (0.36, 0.64) (P < 0.001).

Forest plot based on the random effects model for genotypes (sub assemblages) of Giardia intestinalis in the western half of Iran. (A) AI; (B) AII (C) BIII (D) BІѴ (E) mixed A and B sub assemblages
Forest plot based on the random effects model for genotypes (sub assemblages) of Giardia intestinalis in the western half of Iran. (A) AI; (B) AII (C) BIII (D) BІѴ (E) mixed A and B sub assemblages

Similarly, results of the Egger test demonstrated that the proportion of BIII (Figure 4C) sub-assemblages had no publication bias (P > 0.05), while the proportion of BIV sub-assemblages did (P < 0.05). The Cochran's Q test and I2 statistic for BIII sub-assemblages (I2 = 66.85%, Q = 42.23, P-value < 0.001), and BIV sub-assemblages (I2 = 45.38%, Q = 14.65, P-value <0.1) showed heterogeneity among studies. Therefore, the random effects model was used in the proportion meta-analysis of G. intestinalis sub-assemblages in the western half of Iran. Additionally, based on Figure 4D), the overall proportions of BIII and BIV were estimated to be 0.2 [with a confidence interval of 95% (0.15, 0.26) (P < 0.001)] and 0.06 [with a confidence interval of 95% (0.03, 0.09) (P < 0.001)], respectively.

On the other hand, results of the Egger test indicated that the mixed proportion of different sub-assemblages of G. intestinalis in mentioned region lacked publication bias (P = 0.675). In addition, outcomes of the Cochran's Q test and I2 statistic verified the heterogeneity among the studies (I2 = 97.13%, Q = 209.02, P-value < 0.001). Therefore, the random effects model was used in the proportion meta-analysis for various sub-assemblages of G. intestinalis. Eventually, based on Figure 4E, we can notice that the overall proportion for various G. intestinalis sub-assemblages is 0.18 with a confidence interval of 95% (0.02, 0.43) (P < 0.001).

4. Discussion

In this systematic study, the prevalence of assemblages and sub-assemblages of human G. intestinalis in different geographical regions of Iran was analyzed. The result of this study showed that most of studies on assemblages and sub-assemblages of G. intestinalis have been performed in the western half of Iran. There was little of the same information in the eastern half of this country. The main results of this study are as follows.

The overall prevalence of assemblages A, B, and the mixed assemblage (A and B) in the mentioned area were 0.56, 0.27, and 0.16, respectively. Also, prevalence of AI, AII, BIII and BIV sub-assemblages of were 0.34, 0.5, 0.2 and 0.06, respectively. In addition, the total prevalence of all sub-assemblages (AI, AII and BIII, BIV) was estimated to be 0.18.

Giardiasis is one of the most important health problems worldwide because the parasite that causes this disease, G. intestinalis, is a gastrointestinal protozoan common among human, domestic, and wild animals (3, 5). Hence, the epidemiological study of Giardiasis is important to identify the host spectrum of different species of this parasite as well as its assemblages, sub-assemblages, strains and genotypes. On the other hand, such studies, if performed using molecular methods, will help us improve our understanding about the zoonotic transmission potential of animal’s Giardia species and to determine how many cases of human Giardiasis have animal source (3). Although being sporadic, molecular taxonomic studies can also clarify the relationship between the parasitic and its host genotypes, pathogenesis, and clinical symptoms (33).

Also, different species of Giardia parasite have different hosts. Currently, six species of this parasite have been recognized by researchers. Of these six species, only G. intestinalis (also known as G. lamblia and G. duodenalis) could infect humans and many other mammals. Studies have shown that these species of Giardia have the widest range of hosting and the greatest health importance (34). The results of epidemiological and genotypic studies have also confirmed the possibility of zoonotic transmission of G. intestinalis (33).

Study on allozymes has further shown that all human isolates of G. intestinalis belong to two genetic assemblages (A and B) each of which contains at least four genetic clusters or sub-assemblages (I to IV). On the other hand, both of these assemblages are zoonotic whose examination can help us track the source of infection (35). Further, in Brazil, the findings of two separate studies showed that in communities, where people care for pets, humans and animals (especially cats) could become infected with Giardia assemblage A. In such communities, there seems to be the possibility of cross-transmission of this parasite between animals and humans (36, 37). Another study in Poland verified assemblage A in Giardia infected cats (38). Research works in Mexico, Spain, and Jamaica approved the same assemblage of G. intestinalis (A) in dogs that had close connection with human (39-41). Furthermore, a study conducted to determine the sub-assemblage of G. duodenalis in livestock of Urmia, Iran showed that the predominant assemblage among the livestock of this city is type E, which is not zoonotic (42). Similar results were obtained in other studies conducted on Arabian horses, and ruminants in Ahvaz and Yazd province, Iran, respectively (43, 44). As well as studies performed in Turkey on horse and cattle (45, 46).

As the results of this systematic study showed, prevalence of assemblage A of G. intestinalis with estimated of 0.56 is the dominant assemblage in the human society of western half of Iran. By comparing Giardia sub-assemblages in animals and humans, it can be concluded that in mentioned region of Iran zoonotic transmission of this parasite from livestock and pets is very scarce. In contrast, a general comparison between the prevalence of Giardia in this region of Iran, Brazil, Mexico, Spain, and Jamaica suggested that the source of contamination in these countries may be stray domestic animals, including cats and dogs (36, 39-41). Currently, the large number of stray cats and dogs are the most important urban problems in Iran. These animals can be involved in the transmission and spread of a number of zoonotic parasitic infections including Giardia. Works on the prevalence and genotype determination of Giardia in Iranian domestic animals is so little that we had to report results of the only study in this field in Iran, which showed that the prevalence of G. duodenalis in cats in Ahvaz was 2% and 3.33% by microscopic method and immunochromatography assay, respectively (47).

Moreover, studies in other countries showed that assemblage A of G. intestinalis had the most possibility of zoonotic transmission since it is the most common species of this parasite in animals and it also lacks a specific host. It should be also noted that humans are usually infected with genotype AII, while animals are often infected with type AI (34). Results of three studies conducted in Iran to determine the genotypes of the G. duodenalis in animals also showed that although the predominant assemblage of the parasite in these animals is of category E, a number of these animals were infected with assemblage A and genotype AI (42-44). Thus, it would be likely that the source of zoonotic transmission of Giardia is not livestock (cattle) and horses and there might be other animals responsible for this transmission in Iran.

The results of determining assemblage of Giardia in two neighboring countries of Iran showed that the dominant assemblage in Saudi Arabia and Turkey were A and B, respectively (4).

In addition, the association between the genotypes of the G. intestinalis and clinical symptoms was investigated in a human experimental study. The results showed that the clinical signs of Giardiasis appear only in people infected with the GS / M isolate of the parasite. This insolate belongs to Assemblage B of Giardia (3, 33). Further, results of this study led to a theory stating that pathogenesis of parasite changes with variation of parasite strain. Results of animal experimental studies also reinforced this theory (3, 33). Moreover, a study on the infected Egyptian school children confirmed the symptoms of Giardiasis appeared in patients with Assemblage B of Giardia (48). A similar study in Cuba also confirmed the same assemble (B) (48, 49).

In Iran, several studies have been performed to investigate the relationship between Giardia genotype and the presence or absence of clinical signs. For example, Rafiei et al. showed that the predominant genotype of Giardia in Ahvaz is AII and BIII, but no significant difference was observed between the presence/absence of symptoms and the incidence of these genotypes (17). In another study by Etemadi et al., patients with assemblage A suffered from milder symptoms (mild diarrhea) than those with assemblage B (13).

On the other hand, according to the findings of the current study, it seems that in most studies that have been performed to determine the genotype of Giardia parasite in Iran, in fact, sub-assemblages the parasite were identified and reported.

To achieve more accurate results in Iran and have the possibility for better comparisons with other countries’ studies, it is necessary to conduct more comprehensive studies in different parts of Iran regarding the prevalence and genetic characteristics (assemblage, sub-assemblage and genotype) of G. intestinalis in humans, livestock, and pets.

4.1. Conclusions

Taking into account all the above-discussed studies, it is obvious that the epidemiological transmission pattern of G. intestinalis in Iran is zoonotic transmission. However, it seems important to determine the genetic characteristics of the parasite in order to discover the source of contamination in human societies. Therefore, to reduce the risk of this disease transmission from animal to human (zoonotic), revision of treatment, prevention and control guidelines of this disease in cats and other pets is a priority. On the other hand, the results of the current study shed new light on the prospects and opportunities for future research in the field of epidemiology and genetics of this parasite. Hence, we recommend researchers to use the results of this systematic study to identify areas of Iran that need research in the field of genetic and epidemiological studies of Giardia.

Acknowledgements

References

  • 1.

    Farthing M. Giardiasis. Protozoan diseases. New York, USA: Arnold Publications; 2010. p. 562-78.

  • 2.

    Adam RD. Biology of Giardia lamblia. Clin Microbiol Rev. 2001;14(3):447-75. [PubMed ID: 11432808]. [PubMed Central ID: PMC88984]. https://doi.org/10.1128/CMR.14.3.447-475.2001.

  • 3.

    Feng Y, Xiao L. Zoonotic potential and molecular epidemiology of Giardia species and giardiasis. Clin Microbiol Rev. 2011;24(1):110-40. [PubMed ID: 21233509]. [PubMed Central ID: PMC3021202]. https://doi.org/10.1128/CMR.00033-10.

  • 4.

    Ryan U, Caccio SM. Zoonotic potential of Giardia. Int J Parasitol. 2013;43(12-13):943-56. [PubMed ID: 23856595]. https://doi.org/10.1016/j.ijpara.2013.06.001.

  • 5.

    Ryan UM, Feng Y, Fayer R, Xiao L. Taxonomy and molecular epidemiology of Cryptosporidium and Giardia - a 50 year perspective (1971-2021). Int J Parasitol. 2021;51(13-14):1099-119. [PubMed ID: 34715087]. https://doi.org/10.1016/j.ijpara.2021.08.007.

  • 6.

    Heyworth MF. Giardia duodenalis genetic assemblages and hosts. Parasite. 2016;23:13. [PubMed ID: 26984116]. [PubMed Central ID: PMC4794627]. https://doi.org/10.1051/parasite/2016013.

  • 7.

    Nyaga VN, Arbyn M, Aerts M. Metaprop: a Stata command to perform meta-analysis of binomial data. Arch Public Health. 2014;72(1):39. [PubMed ID: 25810908]. [PubMed Central ID: PMC4373114]. https://doi.org/10.1186/2049-3258-72-39.

  • 8.

    Egger M, Davey Smith G, Schneider M, Minder C. Bias in meta-analysis detected by a simple, graphical test. BMJ. 1997;315(7109):629-34. [PubMed ID: 9310563]. [PubMed Central ID: PMC2127453]. https://doi.org/10.1136/bmj.315.7109.629.

  • 9.

    Guyatt GH, Oxman AD, Kunz R, Woodcock J, Brozek J, Helfand M, et al. GRADE guidelines: 7. Rating the quality of evidence--inconsistency. J Clin Epidemiol. 2011;64(12):1294-302. [PubMed ID: 21803546]. https://doi.org/10.1016/j.jclinepi.2011.03.017.

  • 10.

    Babaei Z, Ourmazdi H, Akhlaghi L, Rezaei S, Razmjou E, Soltani AS, et al. Molecular characterization of the Iranian isolates of Giardia lamblia: application of the glutamate dehydrogenase gene. Iran J Public Health. 2008;37(2):75-82.

  • 11.

    Fallah E, Nahavandi KH, Jamali R, Mahdavi B, Asgharzade M. Genetic Characterization of Giardia intestinalis Strains from Patients Having Sporadic Giardiasis by Using PCR Assay. J Med Sci. 2008;8(3):310-5. https://doi.org/10.3923/jms.2008.310.315.

  • 12.

    Fallah E, Nahavandi KH, Jamali R, Poor BM, Asgharzade M. Molecular Identification of Giardia duodenalis Isolates from Human and Animal Reservoirs by PCR-RFLP. J Biol Sci. 2008;8(5):896-901. https://doi.org/10.3923/jbs.2008.896.901.

  • 13.

    Etemadi S, Ziyaali N, Babaei Z, Fasihi M, Ziyaali A, Salari Z, et al. [The correlation between clinical signs and genotypes of Giardia duodenalis isolated from patients with giardiasis in Kerman city]. J Kerman Univ Med Sci. 2011;18(4):330-9. Persian.

  • 14.

    Pestehchian N, Rasekh H, Babaei Z, Yousefi HA, Eskandarian AA, Kazemi M, et al. Identification of genotypes of Giardia duodenalis human isolates in Isfahan, Iran, using polymerase chain reaction - Restriction Fragment Length polymorphism. Adv Biomed Res. 2012;1:84. [PubMed ID: 23946932]. [PubMed Central ID: PMC3724326]. https://doi.org/10.4103/2277-9175.105166.

  • 15.

    Sarkari B, Ashrafmansori A, Hatam GR, Motazedian MH, Asgari Q, Mohammadpour I. Genotyping of Giardia lamblia isolates from human in southern Iran. Trop Biomed. 2012;29(3):366-71. [PubMed ID: 23018499].

  • 16.

    Manouchehri Naeini K, Hosseini SA, Gholipour A, Babaei Z, Taghipoor S. [Genotyping of Giardia duodenalis isolates in individuals with and without chronic diarrhea using Polymerase Chain Reaction]. J Mazandaran Univ Med Sci. 2012;22(95):39-46. Persian.

  • 17.

    Rafiei A, Roointan ES, Samarbafzadeh AR, Shayesteh AA, Shamsizadeh A, Pourmahdi Borujeni M. Investigation of Possible Correlation between Giardia duodenalis Genotypes and Clinical Symptoms in Southwest of Iran. Iran J Parasitol. 2013;8(3):389-95. [PubMed ID: 24454431]. [PubMed Central ID: PMC3887239].

  • 18.

    Roointan ES, Rafiei A, Samarbaf-Zadeh AR, Shayesteh AA, Shamsizadeh A, Pourmahdi Borujeni M. Genotype Analysis of Giardia lamblia Isolated From Children in Ahvaz, Southwest of Iran. Jundishapur J Microbiol. 2013. https://doi.org/10.5812/jjm.6443.

  • 19.

    Etemadi S, Sarani Ali Abadi P, Babaei Z. Genotype of Giardia duodenalis isolates from humans using glultamate dehydrogenase (gdh) gens by PCRRFLP. Life Sci J. 2013;10(3s):373-7.

  • 20.

    Pestechian N, Rasekh H, Rostami-Nejad M, Yousofi HA, Hosseini-Safa A. Molecular identification of Giardia lamblia; is there any correlation between diarrhea and genotyping in Iranian population? Gastroenterol Hepatol Bed Bench. 2014;7(3):168-72.

  • 21.

    Rayani M, Unyah NZ, Hatam G. Molecular identification of Giardia duodenalis isolates from Fars province, Iran. Iran J Parasitol. 2014;9(1):70.

  • 22.

    Rayani M, Hatam G, Unyah NZ, Ashrafmansori A, Abdullah WO, Hamat RA. Phylogenetic Analysis of Giardia lamblia Human Genotypes in Fars Province, Southern Iran. Iran J Parasitol. 2017;12(4):522-33. [PubMed ID: 29317877]. [PubMed Central ID: PMC5756302].

  • 23.

    Bahrami F, Zamini GH, Haghighi A, Khademerfan MB. Detection and molecular identification of human Giardia isolates in the West of Iran. Biomed Res (India). 2017;28(13):5687-92.

  • 24.

    Nasiri Goorabi L, Pirestani M, Sadraei J. [Genotyping of Giardia duodenalis by β-giardin gene in asymptomatic patients]. J Mazandaran Univ Med Sci. 2017;27(150):27-34. Persian.

  • 25.

    Hooshyar H, Ghafarinasab S, Arbabi M, Delavari M, Rasti S. Genetic Variation of Giardia lamblia Isolates from Food-handlers in Kashan, Central Iran. Iran J Parasitol. 2017;12(1):83-9. [PubMed ID: 28761464]. [PubMed Central ID: PMC5522702].

  • 26.

    Effati F, Dalimi AH, Pirestani M. [A survey on Giardia and Cryptosporidium infection and genotyping common Giardia in children in Alborz Province]. Pathol Res. 2018;21(3):133-9. Persian.

  • 27.

    Kashinahanji M, Haghighi A, Bahrami F, Fallah M, Saidijam M, Matini M, et al. Giardia lamblia assemblages A and B isolated from symptomatic and asymptomatic persons in Hamadan, west of Iran. J Parasit Dis. 2019;43(4):616-23. [PubMed ID: 31749533]. [PubMed Central ID: PMC6841826]. https://doi.org/10.1007/s12639-019-01139-x.

  • 28.

    Shahnazi M, Naghizadeh F, Hajialilo E, Alizadeh SA, Saraei M, Alipour M. Investigation of Giardia intestinalis Genotypes among the Food Handlers of Qazvin, Iran. Iran J Parasitol. 2019;14(4):534-41. [PubMed ID: 32099556]. [PubMed Central ID: PMC7028238].

  • 29.

    Mirrezaie E, Beiromvand M, Tavalla M, Teimoori A, Mirzavand S. Molecular Genotyping of Giardia duodenalis in Humans in the Andimeshk County, Southwestern Iran. Acta Parasitol. 2019;64(2):376-83. [PubMed ID: 30968348]. https://doi.org/10.2478/s11686-019-00051-x.

  • 30.

    Mahmoudi MR, Mahdavi F, Ashrafi K, Forghanparast K, Rahmati B, Mirzaei A, et al. Report of Giardia assemblages and giardiasis in residents of Guilan province-Iran. Parasitol Res. 2020;119(3):1083-91. [PubMed ID: 31965303]. https://doi.org/10.1007/s00436-019-06595-1.

  • 31.

    Rafiei A, Baghlaninezhad R, Koster PC, Bailo B, Hernandez de Mingo M, Carmena D, et al. Multilocus genotyping of Giardia duodenalis in Southwestern Iran. A community survey. PLoS One. 2020;15(2). e0228317. [PubMed ID: 32027684]. [PubMed Central ID: PMC7004373]. https://doi.org/10.1371/journal.pone.0228317.

  • 32.

    Abdi MJ, Yosefi Chigan P, Eslamirad Z, Hajihossein R, Ahmadi A. [Investigating of Human Giardia intestinalis Assemblage in Markazi province, Iran]. Pars Jahrom Univ Med Sci. 2020;18(3):43-7. Persian. https://doi.org/10.52547/jmj.18.3.43.

  • 33.

    Muhsen K, Levine MM. A systematic review and meta-analysis of the association between Giardia lamblia and endemic pediatric diarrhea in developing countries. Clin Infect Dis. 2012;55 Suppl 4:S271-93. [PubMed ID: 23169940]. [PubMed Central ID: PMC3502312]. https://doi.org/10.1093/cid/cis762.

  • 34.

    Monis PT, Caccio SM, Thompson RC. Variation in Giardia: towards a taxonomic revision of the genus. Trends Parasitol. 2009;25(2):93-100. [PubMed ID: 19135417]. https://doi.org/10.1016/j.pt.2008.11.006.

  • 35.

    Monis PT, Thompson RC. Cryptosporidium and Giardia-zoonoses: fact or fiction? Infect Genet Evol. 2003;3(4):233-44. [PubMed ID: 14636685]. https://doi.org/10.1016/j.meegid.2003.08.003.

  • 36.

    Ramirez-Ocampo S, Cotte-Alzate JD, Escobedo AA, Rodriguez-Morales AJ. Prevalence of zoonotic and non-zoonotic genotypes of Giardia intestinalis in cats: a systematic review and meta-analysis. Infez Med. 2017;25(4):326-38. [PubMed ID: 29286011].

  • 37.

    Souza SL, Gennari SM, Richtzenhain LJ, Pena HF, Funada MR, Cortez A, et al. Molecular identification of Giardia duodenalis isolates from humans, dogs, cats and cattle from the state of Sao Paulo, Brazil, by sequence analysis of fragments of glutamate dehydrogenase (gdh) coding gene. Vet Parasitol. 2007;149(3-4):258-64. [PubMed ID: 17900812]. https://doi.org/10.1016/j.vetpar.2007.08.019.

  • 38.

    Piekara-Stepinska A, Piekarska J, Gorczykowski M, Bania J. Genotypes of Giardia duodenalis in Household Dogs and Cats from Poland. Acta Parasitol. 2021;66(2):428-35. [PubMed ID: 33040220]. [PubMed Central ID: PMC8166709]. https://doi.org/10.1007/s11686-020-00292-1.

  • 39.

    Godinez-Galaz EM, Veyna-Salazar NP, Olvera-Ramirez AM, Milian-Suazo F, Perea-Razo CA, Bernal-Reynaga R, et al. Prevalence and Zoonotic Potential of Giardia intestinalis in Dogs of the Central Region of Mexico. Animals (Basel). 2019;9(6):325. [PubMed ID: 31174344]. [PubMed Central ID: PMC6617256]. https://doi.org/10.3390/ani9060325.

  • 40.

    Dado D, Montoya A, Blanco MA, Miro G, Saugar JM, Bailo B, et al. Prevalence and genotypes of Giardia duodenalis from dogs in Spain: possible zoonotic transmission and public health importance. Parasitol Res. 2012;111(6):2419-22. [PubMed ID: 22983168]. https://doi.org/10.1007/s00436-012-3100-x.

  • 41.

    Lee MF, Cadogan P, Eytle S, Copeland S, Walochnik J, Lindo JF. Molecular epidemiology and multilocus sequence analysis of potentially zoonotic Giardia spp. from humans and dogs in Jamaica. Parasitol Res. 2017;116(1):409-14. [PubMed ID: 27885466]. https://doi.org/10.1007/s00436-016-5304-y.

  • 42.

    Malekifard F, Ahmadpour M. Molecular detection and identification of Giardia duodenalis in cattle of Urmia, northwest of Iran. Vet Res Forum. 2018;9(1):81-5. [PubMed ID: 29719668]. [PubMed Central ID: PMC5913565].

  • 43.

    Jafari H, Razi Jalali MH, Seyfi Abad Shapouri M, Haji Hajikolaii MR. Prevalence and genotyping of Giardia duodenalis among Arabian horses in Ahvaz, southwest of Iran. Arch Razi Inst. 2016;71(3):177-81.

  • 44.

    Kiani-Salmi N, Fattahi-Bafghi A, Astani A, Sazmand A, Zahedi A, Firoozi Z, et al. Molecular typing of Giardia duodenalis in cattle, sheep and goats in an arid area of central Iran. Infect Genet Evol. 2019;75:104021. [PubMed ID: 31494270]. https://doi.org/10.1016/j.meegid.2019.104021.

  • 45.

    Demircan K, Onder Z, Duzlu O, Yildirim A, Okur M, Ciloglu A, et al. First Molecular Detection and Phylogenetic Analyses of Zoonotic Giardia intestinalis in Horses in Turkey. J Equine Vet Sci. 2019;80:56-60. [PubMed ID: 31443835]. https://doi.org/10.1016/j.jevs.2019.06.017.

  • 46.

    Onder Z, Simsek E, Duzlu O, Yetismis G, Ciloglu A, Okur M, et al. Molecular prevalence and genotyping of Giardia duodenalis in cattle in Central Anatolia Region of Turkey. Parasitol Res. 2020;119(9):2927-34. [PubMed ID: 32562064]. https://doi.org/10.1007/s00436-020-06771-8.

  • 47.

    Mosallanejad B, Avizeh R, Jalali MR, Alborzi A. Prevalence of Giardia duodenalis Infection in Household Cats of Ahvaz District, South-West of Iran. Iran J Parasitol. 2010;5(3):27-34. [PubMed ID: 22347252]. [PubMed Central ID: PMC3279840].

  • 48.

    Hussein EM, Ismail OA, Mokhtar AB, Mohamed SE, Saad RM. Nested PCR targeting intergenic spacer (IGS) in genotyping of Giardia duodenalis isolated from symptomatic and asymptomatic infected Egyptian school children. Parasitol Res. 2017;116(2):763-71. [PubMed ID: 27975120]. https://doi.org/10.1007/s00436-016-5347-0.

  • 49.

    Jerez Puebla LE, Nunez FA, Santos LP, Rivero LR, Silva IM, Valdes LA, et al. Molecular analysis of Giardia duodenalis isolates from symptomatic and asymptomatic children from La Habana, Cuba. Parasite Epidemiol Control. 2017;2(3):105-13. [PubMed ID: 29774288]. [PubMed Central ID: PMC5952668]. https://doi.org/10.1016/j.parepi.2017.05.003.