|Year : 2022 | Volume
| Issue : 1 | Page : 48-53
Internal transcribed spacer region 1 as a promising target for detection of intra-specific polymorphisms for Strongyloides stercoralis
Yasmin Sultana1, Fanrong Kong2, Mandira Mukutmoni3, Laila Fahria3, Aleya Begum3, Rogan Lee1
1 Centre for Infectious Diseases and Microbiological Services, ICPMR, Westmead Hospital, Westmead NSW Australia; Westmead Clinical School, Faculty of Medicine and Health, University of Sydney, NSW, Australia
2 Centre for Infectious Diseases and Microbiological Services, ICPMR, Westmead Hospital, Westmead NSW Australia, Australia
3 Department of Zoology, Faculty of Biological Sciences, University of Dhaka, Dhaka, Bangladesh
|Date of Submission||09-Mar-2021|
|Date of Decision||18-Aug-2021|
|Date of Acceptance||13-Sep-2021|
|Date of Web Publication||25-Jun-2022|
Centre for Infectious Diseases and Microbiological Services, ICPMR, Westmead Hospital
Source of Support: None, Conflict of Interest: None
| Abstract|| |
Background: Strongyloides stercoralis, the causative agent of strongyloidiasis, is a parasitic worm that has larvae capable of reinfecting the same host. This nematode infection is therefore difficult to treat and to achieve total cure. Information about genetic variation and differences in drug susceptibility between strains is needed to improve treatment outcomes.
Aim: To develop a polymerase chain reaction (PCR) to identify the intra-species variation among 13 S. stercoralis isolates collected from Bangladesh, USA and Australia.
Material & Methods: PCR assays were designed by using primers targeting S. stercoralis internal transcribed spacer (ITS) regions 1 and 2. Sequence data generated by these PCR products were compared to the existing ITS1/2, 18S and 28S rRNA gene sequences in GenBank for phylogenetic analysis.
Results: Intra-species single nucleotide polymorphisms (SNPs) were identified in ITS1 and in the 5.8S rRNA gene. The generated phylogram grouped the 13 isolates into dog, Orangutan and human clusters.
Conclusion: This method could be used as an epidemiological tool to study strain differences in larger collections of S. stercoralis isolates. The study forms the basis for further development of an ITS-based assay for S. stercoralis molecular epidemiological studies
Keywords: Bangladesh, different countries, internal transcribed spacer regions 1 and 2 (ITS1/2), single-nucleotide polymorphism, Strongyloides stercoralis
|How to cite this article:|
Sultana Y, Kong F, Mukutmoni M, Fahria L, Begum A, Lee R. Internal transcribed spacer region 1 as a promising target for detection of intra-specific polymorphisms for Strongyloides stercoralis. Trop Parasitol 2022;12:48-53
|How to cite this URL:|
Sultana Y, Kong F, Mukutmoni M, Fahria L, Begum A, Lee R. Internal transcribed spacer region 1 as a promising target for detection of intra-specific polymorphisms for Strongyloides stercoralis. Trop Parasitol [serial online] 2022 [cited 2022 Aug 18];12:48-53. Available from: https://www.tropicalparasitology.org/text.asp?2022/12/1/48/348298
| Introduction|| |
It is estimated that up to 100 million people, worldwide, are infected with Strongyloides stercoralis. Another subspecies, Strongyloides. fuelleborni kellyi, from nonhuman primates, also infects humans in Africa, Southeast Asia, and Oceania.
Infection is sustained in individuals by repeated episodes of autoinfection, which makes it difficult to eradicate. Chronically infected individuals usually become asymptomatic, but lifelong carriers can develop fatal symptomatic disease if they become immunocompromised. Successful treatment is dependent on 100% clearance of the parasite. Little is known about the genetic differences of S. stercoralis isolates and whether there is any variation in susceptibility to ivermectin of these isolates. Furthermore, the transmission of S. stercoralis in the community would be better understood if molecular typing were available so that sources of infection and person-to-person transmission could be traced.
The internal transcribed spacer (ITS) region has been used as a target for species identification and investigation of intra-specific strain variation of parasitic helminths.,,,, The genetic variation in a cestode parasite, Echinococcus granulosus, showing distinct genotypes within isolates using ITS1 and ITS2, has been documented. However, there is limited information about the genotypes of S. stercoralis. More studies are required to improve the knowledge of parasite transmission and epidemiology. A study was carried out on intraspecies diversity based on the Strongyloides 18S ribosomal RNA gene; it was found to be highly conserved and not suitable for species identification. RFLP analysis of ITS1 and ITS2 polymerase chain reaction (PCR) products was used as the basis for species identification and strain typing of Strongyloides species in human and dog isolates. A study of S. stercoralis in humans from three endemic provinces in Iran also mentioned intraspecies variations based on ITS region characterization, but details were not provided.
In this study, PCR primers targeting ITS regions 1 and 2 were designed based on the existing S. stercoralis ITS1/2 and 18S and 28S rRNA gene sequences in GenBank. The ITS region was analyzed to identify intraspecies polymorphisms in isolates collected from three countries.
| Subject and Methods|| |
Origin of parasites/stool specimens
S. stercoralis larvae from 11 culture-positive specimens were collected from a high-risk group living in a slum community of Dhaka and one isolate from an Australian case was used in this study. S. stercoralis isolated by Harada–Mori and/or agar plate culture methods was identified by morphological characteristics on microscopy. The late Professor G. Schad (University of Pennsylvania, United States) provided S. stercoralis larvae, originally isolated from a patient in West Virginia, USA, and laboratory cultured in dogs. All isolated parasites (n = 13) were stored at −80o C until used.
Isolation of genomic DNA
Total genomic DNA was isolated from the parasites using a SIGMA GenElute Mammalian Genomic DNA Miniprep Kit (Sigma-Aldrich Pty. Ltd, NSW, AUSTRALIA; Catalogue number G 1N350) according to manufacturer's instructions.
Polymerase chain reaction
Three sets of primer pairs were used in this PCR assay, two for ITS1 and one for the ITS2 region [Table 1]. Conserved sequences suitable for primers were inferred from multiple alignments of all the available S. stercoralis and Strongyloides robustus ITS1, ITS2, 28S-like, and 5.8S rDNA sequences in GenBank. Primers ITS2Fm and ITS1F2 bind to the end of the 18S-like rRNA gene. Primer pairs ITS4Rm and ITS3R were designed to bind to the 5'-end of the 28S-like rRNA gene [Figure 1].
|Table 1: Primers for amplification of Strongyloides stercoralis internal transcribed spacer|
Click here to view
|Figure 1: The alignment of Strongyloides stercoralis and Strongyloides robustus genes for 18S rRNA, ITS1, 5.8S rRNA, ITS2, 28S rRNA sequences. Upper lines sequences were generated by integrated GenBank sequences of JF699148 and EF653265. Lower lines is S. robustus GenBank AB272232.1 sequence. 18S rRNA (green), 5.8S rRNA (pink), and 28S rRNA gene (light blue) sequences are shown according to S. robustus GenBank sequence AB272232.1. (Sequences for S. stercoralis in this region were not available in GenBank). The location of listed Primers is shown by yellow highlight and the given primers names are in bold text near the 3'-end of primer sequences|
Click here to view
Single and nested polymerase chain reaction for internal transcribed spacer 1 and single polymerase chain reaction for internal transcribed spacer 2
Primers pair ITS1F2 and ITS3R was used to amplify ITS1. Some samples failed to amplify in the single PCR so a nested PCR was developed to increase the sensitivity of the assay. In the nested PCR, the outer primer pair ITS2Fm and ITS4Rm was used for the first round reaction, followed by a second-round reaction using inner primer pair ITS1F2 and ITS3R. The ITS2 region primers were 5.8SF1 and 28SR1, which bind to the 5.8S rRNA gene and 28S-like rRNA gene, respectively.
The 25-μl PCR mixture for the amplification of both ITS1 and 2 was prepared as follows: 10 μl template DNA and 0.5 μM forward and reverse primers (Sigma-Aldrich, Sydney, Australia) in a final volume of 25 μL of PCR buffer (HostarTaq master mix, Qiagen). The PCR was performed as follows: 95°C for 15 min, 40 cycles of 95°C for 15 s, 60°C for 30 s, and 72°C for 1 min followed by a further extension at 72°C for 5 min. PCR products were separated on a 2% agarose gel and visualized using SYBR safe stain (Invitrogen, Victoria, Australia).
Purification of polymerase chain reaction product, sequence confirmation
PCR products were purified using the GFXTM PCR DNA and Gel Band Purification Kit (GE Healthcare, Buckinghamshire, UK). PCR products were separated on 2% agarose gel and sequenced by ARGF (http://www.agrf.org.au/sequencing.html) in Westmead Millennium Institute (WMI), Westmead, Sydney, Australia. The correct product was confirmed for every amplified sample by sequence analysis and NCBI BLAST alignment.
Internal transcribed spacer sequence analysis
All the sequences (n = 13) were manually checked against the chromatograms for sequencing quality. Sequences were aligned using Mega 5.05 and Bioedit Sequence Alignment Editor Program, version 184.108.40.206 and further aligned with published Strongyloides reference sequences (GenBank reference number EF545004, EF653264, EF653265, EF653266, JF699148, JF699149 U43576, U43577, U43578, and U43579). Phylogenetic analyses were performed using the program PAUP, version 4.0b10.
| Results|| |
Internal transcribed spacer sequence polymorphism
ITS regions of all 13 clinical isolates were successfully amplified by one or both PCRs. Both ITS1 and ITS2 regions of 11 isolates were amplified. However, in two of our isolates (nos. 63 and 169), reliable sequences were obtained only from the ITS1 regions [Table 2]. ITS sequence alignment showed intraspecies single-nucleotide polymorphisms (SNPs) in ITS1 and in the 5.8S rRNA gene. Polymorphism in the ITS2 region could not be checked due to length variations in our sequence products.
|Table 2: Allele types found from the internal transcribed spacer 1/2 region of the RNA gene|
Click here to view
Subtyping and phylogram analysis
Based on informative SNPs (presence/absence combination) in ITS1 and 5.8S region sequences and integrating GenBank references with our own collected isolates, a total of 12 subtypes/allele types were categorized in our study [Table 2].
A dendrogram was generated using PAUP version 4.0.b10 including 13 of our isolates with 10 GenBank sequences for ITS1 and 5.8S regions [Table 2]. The 13 isolates were divided into five groups; all human isolates were separated into three closely related groups. The branch lengths of the dog isolate and isolates collected from Orangutans (GenBank isolates), from the common node, were much longer than human isolates collected from Bangladesh and Australia. This suggests that either important groups are missing from the analyses or that significant variation occurs among isolates from different hosts, and thus, the pattern of the phylogram in terms of strain differences corresponded to host diversity rather than the geographical distribution. All sequences from isolates in this study have been deposited in GenBank with accession numbers JX489140 to JX489153.
| Discussion|| |
The ITS region has been reported to be a reliable diagnostic target for species identification of bacteria, fungi, and both human and animal parasites.,,, An Iranian study used ITS as a species-specific target for the detection of S. stercoralis in stool, but intraspecies diversity was not discussed. Furthermore, there are only ten S. stercoralis ITS sequences in GenBank, of which only four contain both ITS1 and ITS2. More sequence data, from different geographic regions, are needed to define intra-specific variation in S. stercoralis, to which this study makes a useful contribution. To our knowledge, this is the first comparison of S. stercoralis ITS1 region in isolates collected from three countries, Bangladesh, Australia, and the USA.
The SNPs found in our study in ITS 1 gave some differentiation ability, 12 allele types/subtypes were identified in both ITS1 and 5.8S gene regions (comparing our isolates and GenBank sequences based on ITS1 only) [Figure 2]. The diversity of the isolates observed in the phylogram suggested that S. stercoralis strain heterogeneity remains regardless of the geographic location even when a small number of isolates were studied. The branch lengths among the human isolates show there is little genetic divergence among the human isolates. In contrast, the branch lengths to isolates from orangutans and dogs are longer, and indicative of greater genetic divergence between the human isolates and those collected from other hosts. The evolution and dispersion of human strongyloidiasis was illustrated by the phylogenetic analysis of mitochondrial cytochrome c-oxidase subunit 1 gene (cox 1) gene; this study also highlighted the diversity of Strongyloides spp. strains in different hosts from different geographical regions. Dispersal of S. fuelleborni and S. stercoralis has been addressed with a target to find out the causative agent of strongyloidiasis in a Japanese mammalogist. In the earlier study, interpretation of the constructed phylogenetic tree indicates that distribution of S. stercoralis happened with human migration while spread and diversity of S. fuelleborni occurred with movement of Old World primates through Africa and Asia. Furthermore, maximum-likelihood analysis of S. stercoralis grouped the 10 isolates into dog and primate clades which supports findings in the current study.
|Figure 2: Phylogram based on the internal transcribed spacer 1 region and 5.8S partial sequence of 23 strains of Strongyloides stercoralis. Numbers of mutated positions are shown on top of the branches. Bootstrap values above 50 are shown in bold below the branches. Maximum parsimony analysis. PAUP 4.0b10 (Altivec)|
Click here to view
Variation in response to treatment as the result of gene mutation was suggested long ago. Gene involvement and association between genetic variation and anthelmintic resistance was also discussed in Caenorhabditis elegans. No data are available on drug efficacy and the genetic differences in Strongyloides spp. However, a study on closely related species was carried out where 96%–99% nuclear diversity in Trichostrongylus and in Haemonchus contortus isolates was explained as drug resistant. Further evidence of resistance to anthelmintic in different strains of Teladorsagia circumcinata isolated from goats shows that drug efficacy was different in different strains of this helminth parasite.
| Conclusions|| |
Comparison of SNP in the ITS1 region of S. stercoralis isolates can provide genetic information on the worm isolate and help to identify whether the infection source is from animal or human. We can identify differences of S. stercoralis based on ITS sequencing even with this limited number of isolates. However, broader studies are required to determine whether sufficient genetic data can be used to identify the source of infection (i.e., human or zoonotic) and whether there is any genetic association to anthelminthic efficacy.
The study was conducted in collaboration with the Institute of Epidemiology, Disease Control, and Research (IEDCR), Dhaka, Bangladesh. The late Professor G. Schad (University of Pennsylvania, United States) provided the USA S. stercoralis which was originally isolated from a patient in West Virginia and then laboratory cultured in dogs. Dr. Be-Nazir Ahmed (IEDCR, Bangladesh) helped to organize the collection of stool specimens from Bangladesh. We gratefully acknowledge the assistance of Carolina Firacative for data analysis and suggestions of Dr. Deborah Holt, Menzies School of Health Research, Northern Territory, Australia.
Financial support and sponsorship
Financial support was received from NSW CIDM Public Health Funds and the corresponding author was supported by an Australian AusAID scholarship.
Conflicts of interest
There are no conflicts of interest.
| References|| |
Bethony J, Brooker S, Albonico M, Geiger SM, Loukas A, Diemert D, et al.
Soil-transmitted helminth infections: Ascariasis, trichuriasis, and hookworm. Lancet 2006;367:1521-32.
Viney ME. The biology and genomics of Strongyloides
. Med Microbiol Immunol 2006;195:49-54.
Campbell AJ, Gasser RB, Chilton NB. Differences in a ribosomal DNA sequence of Strongylus
species allows identification of single eggs. Int J Parasitol 1995;25:359-65.
Barber KE, Mkoji GM, Loker ES. PCR-RFLP analysis of the ITS2 region to identify Schistosoma haematobium
and S. bovis from Kenya. Am J Trop Med Hyg 2000;62:434-40.
Hung GC, Chilton NB, Beveridge I, Gasser RB. A molecular systematic framework for equine strongyles based on ribosomal DNA sequence data. Int J Parasitol 2000;30:95-103.
Nuchprayoon S, Junpee A, Poovorawan Y, Scott AL. Detection and differentiation of filarial parasites by universal primers and polymerase chain reaction-restriction fragment length polymorphism analysis. Am J Trop Med Hyg 2005;73:895-900.
Labes EM, Nurcahyo W, Deplazes P, Mathis A. Genetic characterization of Strongyloides
spp. From captive, semi-captive and wild Bornean orangutans (Pongo pygmaeus
) in Central and East Kalimantan, Borneo, Indonesia. Parasitology 2011;138:1417-22.
Rinder H, Rausch RL, Takahashi K, Kopp H, Thomschke A, Löscher T. Limited range of genetic variation in Echinococcus multilocularis. J Parasitol 1997;83:1045-50.
Hasegawa H, Sato H, Fujita S, Nguema PP, Nobusue K, Miyagi K, et al.
Molecular identification of the causative agent of human strongyloidiasis acquired in Tanzania: Dispersal and diversity of Strongyloides
spp. and their hosts. Parasitol Int 2010;59:407-13.
Ramachandran S, Gam AA, Neva FA. Molecular differences between several species of Strongyloides
and comparison of selected isolates of S. stercoralis
using a polymerase chain reaction-linked restriction fragment length polymorphism approach. Am J Trop Med Hyg 1997;56:61-5.
Nilforoushan MR, Mirhendi H., Rezaie S, Meamar AR, Kia EB. A DNA-based identification of Strongyloides stercoralis
isolates from Iran. Iran J Public Health 2007;36:16-20.
Sultana Y, Gilbert GL, Ahmed BN, Lee R. Strongyloidiasis in a high risk community of Dhaka, Bangladesh. Trans R Soc Trop Med Hyg 2012;106:756-62.
Garcia LS. Diagnostic Medical Parasitology. Washington, DC. ASM Press; 2001.
Barry T, Colleran G, Glennon M, Dunican LK, Gannon F. The 16s/23s ribosomal spacer region as a target for DNA probes to identify eubacteria. PCR Methods Appl 1991;1:51-6.
Bridge PD, Schlitt T, Cannon PF, Buddie AG, Baker M, Borman AM. Domain II hairpin structure in ITS1 sequences as an aid in differentiating recently evolved animal and plant pathogenic fungi. Mycopathologia 2008;166:1-16.
Prichard RK. Anthelmintic resistance in nematodes: Extent, recent understanding and future directions for control and research. Int J Parasitol 1990;20:515-23.
Blouin MS, Dame JB, Tarrant CA, Courtney CH. Unusual population genetics of a parasitic nematode: Mtdna variation within and among populations. Evolution 1992;46:470-6.
Jackson F, Jackson E, Coop RL. Evidence of multiple anthelmintic resistance in a strain of Teladorsagia circumcincta (Ostertagia circumcincta) isolated from goats in Scotland. Res Vet Sci 1992;53:371-4.
[Figure 1], [Figure 2]
[Table 1], [Table 2]