Gene migration for re-emerging amebiasis in Iran's northwest-Iraq borders: a microevolutionary scale for reflecting epidemiological drift of Entamoeba histolytica metapopulations

Gene migration of giardiasis in Iran; a microevolutionary scale for reflecting transmission patterns of Giardia lamblia assemblages in symptomatic patients

In the microevolutionary scale of Giardia lamblia, the gene migration indicates how G. lamblia assemblages have transmitted between adjacent counties. 33 positive fecal samples were taken from patients suffering gastrointestinal disorders (nausea, bloating, burping constipation and fatty diarrhea) at Tabriz and Ardabil cities, where located in the cold regions of northwest Iran. Following parasitological examinations, DNA samples were extracted, amplified and digested by single-step PCR-RFLP assay, targeting the glutamate dehydrogenase (gdh) locus to distinguish within and between assemblages A and B. PCR products were directly sequenced to reconfirm their heterogeneity traits and phylogenetic analysis. Of the 33 isolates, 81.9% (n: 27), 9% (n: 3) and 9% (n: 3) were successfully identified as assemblages A (genotype AII), B (genotype BIII) and the mixed of genotypes AII and B, respectively. Despite the presence of heterogeneous clinical backgrounds, a low genetic diversity of sub-assemblage AII was identified among symptomatic cases. A low value of pairwise fixation index showed that G. lamblia sub-assemblage AII is not genetically differentiated among northwest regions of Iran. The occurrence of haplotypes TAB-1/ARD-1 between two regional populations indicates that there is a dawn of G. lamblia gene flow due to transfer of alleles through host mobility and/or ecological alterations. To assess the hypothetical evolutionary scenario, further studies are essential for multilocus genotyping of G. lamblia in tropical regions of Iran and neighboring countries.

Differential detection of Entamoeba histolytica, Entamoeba dispar and Entamoeba moshkovskii in faecal samples using nested multiplex PCR in west of Iran

This study aimed to determine the prevalence of Entamoeba histolytica, Entamoeba dispar and Entamoeba moshkovskii (collectively referred to as Entamoeba complex), using microscopic and molecular methods in Kurdistan Province, northwest of Iran. The relationship between positive Entamoeba species and clinical symptoms was also investigated. Eight positive Entamoeba complex, as well as four Entamoeba complex-like isolates, were detected by microscopic stool examination. DNA was extracted from all positive and from 55 randomly selected negative stool samples. PCR was performed using species-specific 18S rRNA primers for the Entamoeba complex. All positive PCR samples were sequenced. In total, 14 (1.01%) out of 1383 isolates, i.e. 12 microscopy-positive and Entamoeba complex-like isolates and two out of 55 microscopy-negative isolates, were identified via PCR and sequencing. Overall, 0.58% (8/1383) of the isolates were E. dispar, 0.14% (2/1383) E. histolytica, 0.07% (1/1383) E. moshkovskii and 0.22% (3/1383) were mixed of E. histolytica and E. dispar. Based on our findings, the prevalence of E. dispar is greater than that of E. histoltyica. On the other hand, a case of E. moshkovskii was reported for the first time in this region. It seems that some gastrointestinal symptoms may be attributed to Entamoeba species.

Amoebiasis in Iran: a systematic review and meta-analysis

A comprehensive meta-analysis study was performed to estimate the reliable national prevalence and molecular epidemiology of amoebiasis in Iran. Nine English and Persian databases were searched to achieve the relevant studies. Pooled estimates were generated and meta-regression was performed. We identified 71 eligible articles involving 330 930 subjects from 25 provinces to be included in the final analysis. Moreover, 17 studies compromising 462 polymerase chain reaction (PCR)-positive isolates performed molecular analysis to inter-species differentiation. The pooled prevalence of Entamoeba infection among Iranian population was about 1% (95% CI 0.8-2.0%). Moreover, regarding Human Development Index (HDI), a higher prevalence was observed in undeveloped provinces. Out of 462 PCR-positive isolates, 83% (95% CI 69-94%) and 12% (95% CI 3-24%) were Entamoeba dispar, Entamoeba histolytica, respectively. In subgroup analysis based on molecular results, in general, population prevalence of Entamoeba dispar and E. histolytica were 91% (95% CI 80-99%) and 7%, (95% CI 0-19%), respectively, while prevalence of these species in patients with gastrointestinal disorders were 75% (95% CI 45-96%) and 18% (95% CI 1-43%), respectively. Our findings indicate the low burden of amoebiasis in Iran. E. dispar, that is mostly non-pathogenic, was identified as most prevalent species. Nevertheless, we suggest more public health interventions in areas with lower HDI.

Molecular evaluation of pvdhfr and pvmdr-1 mutants in Plasmodium vivax isolates after treatment with sulfadoxine/pyrimethamine and chloroquine in Iran during 2001-2016

The rising use of sulfadoxine/pyrimethamine (SP) in the treatment of chloroquine (CQ)-resistant Plasmodium falciparum has resulted in increased exposure to P. vivax isolates in Iran, where both species are being circulated. In this investigation, the frequency of pvdhfr and pvmdr-1 mutants was assessed in P. vivax strains during 2001-2016 after the introduction of SP/CQ in malarious areas of Iran. The P. vivax isolates (n, 52) were obtained from autochthonous samples in Southeast Iran during 2015-2016. The genomic DNA was extracted and examined using nested polymerase chain reaction-(PCR) and sequencing. Mutations were detected in pvdhfr codons P33L (21.2%), T61 M (25%), S93H (3.9%), and S117 T (1.9%) and 5 isolates showed double mutations (33 L/61 M, 7.7%; 33 L/117 T, 1.9%). No mutation was identified in pvdhfr codons F57 and S58. The pvmdr-1 1076 L mutation was detected in 93.3% of P. vivax isolates. The findings indicated that the frequency of three codons of pvdhfr F57/S58/S117 has decreased from 2001 (1.05%/7.0%/16.9%) to 2016 (0%/0%/1.9%). Genomic analysis of pvmdr-1 showed that the frequency of 1076 L has gradually increased from 2013 (93%) to 2016 (93.3%) (P > .05). The results demonstrated that P. vivax isolates are probably being exited under SP pressure, which reflects the appropriate level of training for field microscopists, as established by Iranian policymakers. Emergent pvdhfr codons 33L, 61M, and 93H should be noticed in plausible drug tolerance and treatment plans. The high prevalence of pvmdr-1 1076L mutation shows that efficacy of CQ combination with primaquine may be in danger of being compromised, however further investigations are needed to evaluate the clinical importance of CQ-resistant P. vivax isolates.

Genetic variability and transcontinental sharing of Giardia duodenalis infrapopulations determined by glutamate dehydrogenase gene

Microevolutionary data of Giardia duodenalis sub-assemblages is a prerequisite for determining the invasion zoonotic patterns of the parasite. To infer transmission patterns that could not be differentiated by the phenotypic features, a population genetic investigation is crucial for the elucidation of the genetic structure of G. duodenalis among the continents. Forty G. duodenalis positive fecal samples were collected from different foci of Northwest Iran. The specimens were subjected to Trichrome staining and sucrose gradient flotation. DNA samples were extracted, amplified, and sequenced by targeting glutamate dehydrogenase (gdh) gene. The global gdh sequences of sub-assemblages AII and BIV retrieved from NCBI GenBank were analyzed to estimate diversity indices, neutrality indices, and gene migration tests. Sequencing analyses indicated various levels of genetic variability of sub-assemblages AII and BIV among the five continents. Sub-assemblage BIV had greater genetic variability (haplotype diversity: 0.975; nucleotide diversity: 0.04246) than sub-assemblage AII. The statistical Fst value demonstrated that the genetic structure of sub-assemblages AII and BIV are moderately differentiated between European-American populations (Fst: 0.05352-0.15182), whereas a significant differentiation was not seen among other geographical population pairs. We conclude that a high gene flow of G. duodenalis sub-assemblages AII and BIV is unequivocally sharing among the continents. The current findings strengthen our knowledge to assess the evolutionary patterns of G. duodenalis in endemic foci of the world and it will become the basis of public health policy to control human giardiasis.

Phylogeography, genetic variability and structure of Acanthamoeba metapopulations in Iran inferred by 18S ribosomal RNA sequences: A systematic review and meta-analysis

OBJECTIVE

To verify phylogeography and genetic structure of Acanthamoeba populations among the Iranian clinical isolates and natural/artificial environments distributed in various regions of the country.

METHODS

We searched electronic databases including Medline, PubMed, Science Direct, Scopus and Google Scholar from 2005 to 2016. To explore the genetic variability of Acanthamoeba sp, 205 sequences were retrieved from keratitis patients, immunosuppressed cases and environmental sources as of various geographies of Iran.

RESULTS

T4 genotype was the predominant strain in Iran, and the rare genotypes belonged to T2, T3, T5 (Acanthamoeba lenticulata), T6, T9, T11, T13 and T15 (Acanthamoeba jacobsi). A total of 47 unique haplotypes of T4 were identified. A parsimonious network of the sequence haplotypes demonstrated star-like feature containing haplogroups IR6 (34.1%) and IR7 (31.2%) as the most common haplotypes. In accordance with the analysis of molecular variance, the high value of haplotype diversity (0.612-0.848) of Acanthamoeba T4 represented genetic variability within populations. Neutrality indices of the 18S ribosomal RNA demonstrated negative values in all populations which represented a considerable divergence from neutrality. The majority of genetic diversity belonged to the infected contact lens and dust samples in immunodeficiency and ophthalmology wards, which indicated potential routes for exposure to a pathogenic Acanthamoeba sp. in at-risk individuals. A pairwise fixation index (F_ST) was from low to high values (0.02433-0.41892). The statistically F_ST points out that T4 is genetically differentiated between north-west, north-south and central-south metapopulations, but not differentiated between west-central, west-south, central-south, and north-central isolates.

CONCLUSIONS

An occurrence of IR6 and IR7 displays that possibly a gene flow of Acanthamoeba T4 occurred after the founder effect or bottleneck experience through ecological changes or host mobility. This is the first systematic review and meta-analysis providing new approaches into gene migration and transmission patterns of Acanthamoeba sp, and targeting at the high-risk individuals/sources among the various regions of Iran.