|Year : 2020 | Volume
| Issue : 11 | Page : 515-520
Genetic variation and phylogenetic relationship of Hypoderaeum conoideum (Bloch, 1782) Dietz, 1909 (Trematoda: Echinostomatidae) inferred from nuclear and mitochondrial DNA sequences
Chairat Tantrawatpan1, Weerachai Saijuntha2
1 Division of Cell Biology, Department of Preclinical Sciences, Faculty of Medicine, Thammasat University, Rangsit Campus, Pathumthani 12120, Thailand
2 Walai Rukhavej Botanical Research Institute, Mahasarakham University, Maha Sarakham 44150, Thailand
|Date of Submission||05-Apr-2020|
|Date of Decision||15-May-2020|
|Date of Acceptance||16-May-2020|
|Date of Web Publication||30-Sep-2020|
Walai Rukhavej Botanical Research Institute, Mahasarakham University, Maha Sarakham 44150
Source of Support: This research was supported by Faculty of Medicine, Thammasat University, Thailand to CT, grant number 2-18/2562, Conflict of Interest: None
Objective: To explore genetic variations of Hypoderaeum conoideum collected from domestic ducks from 12 different localities in Thailand and Lao PDR, as well as their phylogenetic relationship with American and European isolates.
Methods: The nucleotide sequences of their nuclear ribosomal DNA (ITS), mitochondrial cytochrome c oxidase subunit 1 (CO1), and NADH dehydrogenase subunit 1 (ND1) were used to analyze genetic diversity indices.
Results: We found relatively high levels of nucleotide polymorphism in ND1 (4.02%), whereas moderate and low levels were observed in CO1 (2.11%) and ITS (0.96%), respectively. Based on these polymorphisms, the 20 ND1, 12 CO1, and 18 ITS haplotypes were classified, and several common haplotypes were observed in all samples. At least three major lineages, namely American, European and Asian lineages, have been classified by phylogenetic analyses based on ND1 sequences.
Conclusions: Our report demonstrates that the ND1 gene is the most suitable genetic marker to explore genetic variation and phylogenetic relationship of Hypoderaeum conoideum. However, a combination of all loci for ND1, CO1 and ITS would be of great value toward further genetic investigation of this endemic worldwide parasite. Thus, comprehensive molecular genetic analyses of Hypoderaeum conoideum from its worldwide distribution is needed to further understanding of the evolutionary and systematic relationships of this parasite.
Keywords: Echinostomes; Genetic diversity; Genetic differentiation; Nuclear ribosomal DNA; Cytochrome c oxidase subunit 1; NADH dehydrogenase subunit 1
|How to cite this article:|
Tantrawatpan C, Saijuntha W. Genetic variation and phylogenetic relationship of Hypoderaeum conoideum (Bloch, 1782) Dietz, 1909 (Trematoda: Echinostomatidae) inferred from nuclear and mitochondrial DNA sequences. Asian Pac J Trop Med 2020;13:515-20
|How to cite this URL:|
Tantrawatpan C, Saijuntha W. Genetic variation and phylogenetic relationship of Hypoderaeum conoideum (Bloch, 1782) Dietz, 1909 (Trematoda: Echinostomatidae) inferred from nuclear and mitochondrial DNA sequences. Asian Pac J Trop Med [serial online] 2020 [cited 2022 May 18];13:515-20. Available from: https://www.apjtm.org/text.asp?2020/13/11/515/295362
| 1. Introduction|| |
Hypoderaeum (H.) conoideum (Bloch, 1782) Dietz, 1909 is a species of digenetic trematode in the family Echinostomatidae, which is the causative agent of human and animal echinostomiasis,. The life cycle of H. conoideum has been extensively studied experimentally,. A wide variety of freshwater snails, especially the lymnaeid species, are the first intermediate hosts and shed cercariae. Aquatic animals, such as bivalves, fishes, and tadpoles, often act as second intermediate hosts. Eating these raw or partially-cooked aquatic animals has been identified as the primary mode of transmission,. It has been noted that wild rats, chickens, ducks, and geese may be infected and then consumed,,. Heavy infection of H. conoideum was observed in free-grazing ducks in Thailand and Lao PDR, and there have been reports of human echinostomiasis in Thailand,.
To date, at least four medically important echinostome species have been documented in Thailand: Echinostoma (E.) revolutum (Froelich, 1802), Echinostoma ilocanum (Garrison, 1908), Artyfechinostomum (A.) malayanum (Leiper, 1911) Mendheim, 1943, and H. conoideum. However, as the eggs of all these echinostomes are morphologically very similar, there is extreme difficulty in differential diagnosis. Thus, several molecular markers/techniques have been developed to solve this problem,,,. Molecular markers such as allozymes as well as nuclear and mitochondrial DNA sequences have been effective in the genetic study of echinostomes,,,. Nuclear ribosomal DNA, mitochondrial cytochrome c oxidase subunit 1 (CO1) and NADH dehydrogenase subunit 1 (ND1) have been widely used to explore genetic variation of echinostomes. CO1 and ND1 sequences were used to differentiate E. revolutum and E. miyagawai infections in free-grazing ducks. Recently, species-specific reverse primers have been designed to target different ND1 fragment sizes in multiplex polymerase chain reaction (PCR) to differentiate four echinostomes.
Previous reports demonstrate that the genetic variations present in these echinostomes play an important role when developing these molecular markers/techniques. While comprehensive genetic investigation has been undertaken in some species, e.g. the “37 collar-spined” or the “revolutum” group including E. miyagawai and E. revolutum,,, there have been few similar investigations into H. conoideum. To our knowledge, very few published reports of genetic investigation into H. conoideum, e.g. mitochondrial genome sequencing and genetic differentiation from E. revolutum , outlined possible phylogenetic relationships with some species of the family Echinostomatidae using nuclear DNA sequence and karyological analysis,. Therefore, our study aims to use nuclear ribosomal DNA (ITS), mitochondrial CO1, and the ND1 gene to explore H. conoideum genetic variations in domestic ducks from Thailand and Lao PDR.
| 2. Materials and methods|| |
2.1. Sample collection and DNA extraction
Duck intestines were purchased from local slaughterhouses in several localities around Thailand and from one location in Vientiane, Lao PDR [Table 1] and [Figure 1]. The intestines were then opened to retrieve adult H. conoideum. The adult worms were identified by flattening them between two glass slides, held apart by a small piece of filter paper, and then examined under a light microscope. Species were identified based on elongate-shaped testis, presence of small collar and 45-53 collar-spines. Adult worms were then washed several times in normal saline before being fixed in molecular grade ethanol and stored at –20 °C until use. A single worm was placed into a 1.5 mL vial, and the alcohol removed by vacuum centrifuge. Then an E.Z.N.A.® Tissue DNA kit (Omega Bio-tek, USA) was used to extract DNA following the manufacturer’s instructions. The biosafety ethics for this research were approved by the Institute Biosafety Committee of Thammasat University: 173/2561.
|Figure 1: Map of sampling localities. More details of the locality codes are shown in Table 1.|
Click here to view
2.2. PCR and DNA sequencing
To amplify the DNA target regions, PCR was performed using primer pairs of ITS1F (5’-GGTGAACCTGCGGAAGGATC-3’) and ITS2R (5’-AGCGGGTAATCACGACTG-3’) for ITS1-5.8S-ITS2 fragments. The primers for CO1 and ND1 amplification are current design and use in this study, namely Hcon-CO1F (5’-GATGCCGGTGTTGATTGGTG-3’) and Hcon-COR (5’-CACACAAACCCGACGAGGTA-3’) for CO1, Hcon-ND1F (5’A TTCTGCTTACTTAGTTTTGAGGAG-3’) and Hcon-ND1R (5’-TAATCATAACGAACCCGAGGTAAC-3’) for ND1. For all DNA fragment amplification, PCR was performed in 25 μL, including approximately 100 ng of template DNA, 1 μL of each primer (each 10 pmol/μL), 2.5 μL of 10× Ex-Taq buffer (Mg2+ plus), 2 μL dNTPs (each 2.5 mmol/L), 0.125 μL of Ex-Taq DNA polymerase (5 U/μL), and 19.375 μL of deionized water. PCR employed an initial denaturing step at 95 °C for 4 min, 35 cycles of denaturing at 94 °C for 40 s, annealing at 50 °C for 60 s, extending at 72 °C for 60 s, and a final extending step of 72 °C for 8 min. PCR products were separated into 1% agarose gel electrophoresis, and visualized with GelRed™ Nucleic Acid Gel Stain (Biotium, Inc., Hayward, CA). The amplified bands were cut and purified using E.Z.N.A.® Gel Extraction kit (Omega Bio-tek, USA) and subsequently sequenced at 1st BASE Laboratories, Malaysia.
2.3. Data analyses
The sequences were checked and edited using ABI sequence scanner v1.0 and BioEdit v7.2.6 software. Multiple sequence alignment was performed using the ClustalW program. Genetic diversity indices, including segregated sites, haplotype number, haplotype diversity, and nucleotide diversity, were computed and generated by DnaSp v5 software. Median joining haplotype networks were generated using Network v.5.0.11 (http://www. fluxus-engineering.com) based on CO1, ND1, and ITS sequences. Phylogenetic relationship was performed based on ND1 sequences of H. conoideum from different localities in Thailand and Lao PDR, as well as a sequence of H. conoideum from Hubei, China, including Hypoderaeum sp. from America and Europe available in GenBank. The maximum likelihood and neighbor joining trees were performed using MEGA X with nodal support estimated using 1 000 bootstrap re-sampling.
| 3. Results|| |
The new sequences generated in this study were deposited in GenBank with the accession numbers MT127126 to MT127176, MT159471 to MT159526 and MT175378 to 175434 for the partial sequences of ITS1-5.8S-ITS2, CO1 and ND1, respectively. In the case of ND1, the 820 bp fragment was successfully amplified, but only 522 bp were readable and reliable for analyses after DNA sequencing. As the variable sites of ITS1-5.8S-ITS2 were found in the ITS1 and ITS2 regions. Although all samples underwent PCR to amplify all three DNA regions, some were not successfully amplified or did not provide reliable DNA sequencing results in all regions. Thus, we have an unequal number of samples examined in each locus [Table 2].
|Table 2: Molecular diversity indices and neutrality test of ITS, CO1 and ND1 sequences in the Hypoderaeum conoideum populations from Thailand and Lao PDR.|
Click here to view
Molecular diversity index comparisons between three loci found that ND1 had the highest nucleotide polymorphism with variable nucleotide sites at 4.023%, whereas CO1 and ITS sequences had 2.111% and 0.964%, respectively [Table 2]. Of these variable sites, nucleotide transition was most commonly observed in all three DNA regions, whereas nucleotide transversion were only occasionally present in CO1 and ND1 genes. Insertion/deletion (indel) was only visible in the ITS regions at four nucleotide sites [Table 2]. Nucleotide and haplotype (gene) diversity were observed at (0.007 76 ± 0.000 60) and (0.927 ± 0.017) for ND1, (0.004 02 ± 0.000 40) and (0.789 ± 0.038) for CO1, (0.001 10 ± 0.000 15) and (0.678 ± 0.041) for ITS sequences analyses. Neutrality analysis found that negative values were observed in all Tajima’s D, Fu and Li’s D, and Fu and Li’s F calculations [Table 2], revealing an excess of rare alleles due to a recent population expansion. There were 20 haplotypes of ND1, 12 of CO1, and 18 of ITS; of these, 13 ND1, 9 CO1, and 10 ITS haplotypes were unique to a particular locality [Table 2]. Several common haplotypes were networked between different localities, such as haplotypes C1, C2 and C10 for CO1; haplotypes I6, I3, I12, I15, I16 and I17 for ITS; haplotypes N1, N6, N7 and N14 for ND1 [Figure 2].
|Figure 2: Haplotype networks constructed from CO1 (A), ITS (B) and ND1 (C) sequences. Each circle represents a haplotype named with a different number. Different printed patterns in the haplotype networks correspond to their geographical localities separated into 12 different localities.|
Click here to view
Phylogenetic tree constructed based on ND1 sequence [Figure 3] demonstrated that the genetic structure of H. conoideum is related to continent geographical scale and defined into three different lineages supported by high bootstrap values. All samples from Thailand and Lao PDR were closely aligned with a sequence of H. conoideum from Hubei, China and defined as Asian lineage. Whereas the Hypoderaeum sp. ND1 sequences from USA and Canada were monophyletic group as American lineage and closely aligned with two sequences from Spain which was classified as European lineage [Figure 3].
|Figure 3: Phylogenetic tree constructed based on ND1 sequences of Hypoderaeum conoideum from Thailand, Lao PDR and China, including Hypoderaeum sp. form Spain, USA and Canada. Nodal supports are bootstrap values obtained first by neighbor joining and then maximum likelihood analyses.|
Click here to view
| 4. Discussion|| |
To the best of our knowledge, this study is the first to use three genetic markers including two mitochondrial genes, i.e. CO1 and ND1, and nuclear DNA sequences, i.e. ITS1-5.8S-ITS2, to explore genetic variations in H. conoideum recovered from domestic ducks in Thailand and Lao PDR. Nucleotide polymorphism of H. conoideum was detected at all the loci examined. The mitochondrial ND1 sequence revealed the highest polymorphism and was most suitable for use as a genetic marker to explore H. conoideum genetic variations, as also previously reported in the genetic variation investigations of the “37 collar-spined” or the “revolutum” group of Echinostoma,. Some combination of all three loci incorporated within our study could potentially be of great use for further genetic investigation of this parasite, which is endemic worldwide.
We found a high amount of genetic variation in H. conoideum populations from the different geographical localities. We found many common or shared haplotypes of H. conoideum in Thailand and Lao PDR. This indicates there has clearly been some gene flow among the populations between these two regions. This may have occurred as farmers in both areas feed their free-grazing ducks on naturally occurring snails by rotating them among rice paddies in Thailand,. It is likely the same in Lao PDR.
As expected, there was no genetic sub-structuring related to different regions in Thailand and Lao PDR observed in our study. However, there were reports of geographically influenced genetic differentiation within E. revolutum and E. miyagawai at the continental scale, such as in the Americas and Australia; isolates from Europe and Southeast Asia are defined as the Eurasian population,. We found that the lineages of H. conoideum and Hypoderaeum sp. were difference at a continental scale by separated into American, European and Asian lineages supporting the previous reports of genetic structured of the 37 collar-spined echinostomes,.
Genetic variations of H. conoideum populations in Thailand and Lao PDR may potentially be caused by several factors, such as specialization of infecting different species of intermediate and final hosts. Many aquatic animals can act as intermediate hosts, and the final hosts could be birds, rodents or other mammals, including humans. The various species or cryptic species of snail intermediate hosts located in a particular area may contain differently adapted genotypes/haplotypes. As previously reported, a unique haplotype based on the ND1 sequence was noted in A. malayanum recovered from different snails acting as intermediate hosts. Moreover, temporal factors such as seasonal/climatic variations may also have influenced the genetic differentiation and sub-structuring of H. conoideum populations. Previous evidence of the genetic differentiation in E. revolutum populations from the same locality in Thailand related to temporal factors was reported. To facilitate more understanding of the population genetics and structure of this parasite, H. conoideum samples should be collected from different intermediate/final hosts as well as at different times of year in a particular locality and then genetically analyzed.
In Southeast Asia, there is still a high prevalence of echinostome infection in humans, particularly in Lao PDR, Vietnam, Cambodia, and Thailand,,,. Co-endemic and mixed infection of echinostomes, e.g. H. conoideum and E. revolutum in duck hosts or H. conoideum and A. malayanum in rodent hosts (Ribas et al, unpublished) have been observed. All of these previous studies have identified echinostome specimens solely based on morphology; in a recent report, however, multiplex PCR, based on species-specific reverse primers and targeting variable regions of the ND1 gene, was able to differentially diagnose human echinostomes in Thailand. Other mitochondrial genes could also be used to specifically detect E. revolutum metacercariae in their snail hosts,. The basic information of nucleotide sequence variations from our study can be applied in the molecular diagnosis/detection of H. conoideum as well as other medically important echinostomes. It would be worthwhile to attempt the molecular identification of future human echinostomiasis cases in Southeast Asia. Not only would this increase our understanding of the host-parasite life cycle, but it may hopefully provide us with valuable information on how to interrupt its maturation and, thus prevent future infestations in intermediate and human hosts.
Conflict of interest statement
The authors declare they have no conflict of interest.
We would like to thank Dr. Adrian R. Plant for English proofreading.
This research was supported by Faculty of Medicine, Thammasat University, Thailand to CT, grant number 2-18/2562.
CT and WS contributed equally to this article and read and approved the final manuscript.
| References|| |
Toledo R, Esteban JG. An update on human echinostomiasis. Trans R
Soc Trop Med Hyg
Chai JY. Echinostomes. In: Chai JY (ed.) Human intestinal flukes from discovery to treatment and control
. Netherland: Springer Nature BV; 2019, p. 169-344.
Muñoz-Antoli C, Toledo R, Esteban JG. The life cycle and transmission dynamics of the larval stages of Hypoderaeum conoideum. J Helminthol
Azizi H, Farahnak A, Mobedi I, Molaei Rad MB. Experimental life cycle of Hypoderaeum conoideum
(Block, 1872) Diez, 1909 (Trematoda: Echinostomatidae) parasite from the North of Iran. Iran J Parasitol
Toledo R, Muñoz-Antolí C, Pérez M, Esteban JG. Miracidial infectivity of Hypoderaeum conoideum
(Trematoda: Echinostomatidae): Differential susceptibility of two lymnaeid species. Parasitol Res
Chai JY. Echinostomes in humans. In: Fried B, Toledo R (eds.) The biology of echinostomes: From the molecule to the community
. New York: Springer; 2009, p. 147-183.
Rim HJ. Handbook Series in Zoonoses
. Section C: Parasitic zoonoses, Vol. III. Trematode zoonoses. Boca Raton: CRC Press; 1982.
Saijuntha W, Duenngai K, Tantrawatpan C. Zoonotic echinostome infections in free-grazing ducks in Thailand. Korean J Parasitol
Toledo R, Muñoz-Antoli C, Esteban JG. Intestinal trematode infections. Ad Exp Med Biol
Yokogawa M, Harinasuta C, Charoenlarp P. Hypoderaeum conoideum
(Block, 1782) Diez, 1909, a common intestinal fluke of man in the northeast Thailand. Jpn J Parasitol
Harinasuta T, Bunnag D, Radomyos P. Intestinal fluke infections. In: Pawlowski ZS (ed.) Intestinal helminthic infections: Bailliere’s clinical tropical medicine and communicable diseases
. London: Bailliere Tindall; 1987, p. 695-721.
Tantrawatpan C, Saijuntha W. Multiplex PCR development for the differential detection of four medically important echinostomes (Trematoda: Echinostomatidae) in Thailand. Acta Trop
Georgieva S, Selbach C, Falt nková A, Soldánová M, Sures B, Skínisson K, et al. New cryptic species of the ‘revolutum’
group of Echinostoma
(Digenea: Echinostomatidae) revealed by molecular and morphological data. Parasit Vectors
Saijuntha W, Tantrawatpan C, Sithithaworn P, Duenngai K, Agatsuma T, Andrews RH, et al. Genetic characterization and phylogenetic analysis of echinostomes from free-grazing ducks in Thailand using ribosomal DNA sequences. Southeast Asian J Trop Med Public Health
Nagataki M, Tantrawatpan C, Agatsuma T, Sugiura T, Duenngai K, Sithithaworn P, et al. Mitochondrial DNA sequences of 37 collar-spined echinostomes (Digenea: Echinostomatidae) in Thailand and Lao PDR reveals presence of two species: Echinostoma revolutum
and E. miyagawai. Infect Genet Evol
Saijuntha W, Sithithaworn P, Andrews RH. Genetic differentiation of Echinostoma revolutum
and Hypodereaum conoideum
from domestic ducks in Thailand by multilocus enzyme electrophoresis. J Helminthol
Saijuntha W, Tantrawatpan C, Sithithaworn P, Andrews RH, Petney TN. Spatial and temporal genetic variation of Echinostoma revolutum
(Trematoda: Echinostomatidae) from Thailand and the Lao PDR. Acta Trop
Buddhachat K, Chontananarth T. Is species identification of Echinostoma revolutum
using mitochondrial DNA barcoding feasible with high-resolution melting analysis? Parasitol Res
Yang X, Gasser RB, Koehler AV, Wang L, Zhu K, Chen L, et al. Mitochondrial genome of Hypoderaeum conideum
comparison with selected trematodes. Parasit Vectors
Staneviciute G, Stunženas V, Petkeviciute R. Phylogenetic relationships of some species of the family Echiostomatidae Odner, 1910 (Trematoda), inferred from nuclear rDNA sequences and karyological analysis. Comp Cytogenet
Radomyos P, Tangtrongchitr A, Krudsood S, Wilairatana P, Looareesuwan S. Atlas of medical parasitology
[in Thai]. 7th ed. Bangkok: Parbpim press; 2004.
Areekit S, Singhaphan P, Khuchareontaworn S, Kanjanavas P, Sriyaphai T, Pakpitchareon A, et al. Intraspecies variation of Brugia
spp. in cat reservoirs using complete ITS sequences. Parasitol Res
Hall TA. BioEdit: A user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. in Nucleic acids symposium series
. London: Information Retrieval Ltd; 1999, p. 1979-2000.
Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam H, et al. Clustal W and Clustal X version 2.0. Bioinformatics
Librado P, Rozas J. DnaSP v5: A software for comprehensive analysis of DNA polymorphism data. Bioinformatics
Bandelt HJ, Forster P, Rohl A. Median-joining networks for inferring intraspecific phylogenies. Mol Biol Evol
Kumar S, Stecher G, Li M, Knyaz C, Tamura K. MEGA X: Molecular evolutionary genetics analysis across computing platforms. Mol Biol Evol
Tantrawatpan C, Saijuntha W, Sithithaworn P, Andrews RH, Petney TN. Genetic differentiation of Artyfechinostomum malayanum
and A. sufrartyfex
(Trematoda: Echinostomatidae) based on the internal transcribed spacer sequences. Parasitol Res
Chai JY, Sohn WM, Yong TS, Eom KS, Min DY, Hoang EH, et al. Echinostome flukes recovered from humans in Khammouane Province, Lao PDR. Korean J Parasitol
Sohn WM, Chai JY, Yong TS, Eom KS, Yoon CH, Sinuon M, et al. Echinostoma revolutum
infection in children, Pursat Province, Cambodia. Emerg Infect Dis
Anucherngchai S, Chontananarth T, Tejangkura T, Chai JY. The study of Cytochrome B (CYTB): Species-specific detection and phylogenetic relationship of Echinostoma revolutum
, (Froelich, 1802). J Parasit Dis
[Figure 1], [Figure 2], [Figure 3]
[Table 1], [Table 2]
|This article has been cited by|
||Intercontinental distributions, phylogenetic position and life cycles of species of Apharyngostrigea (Digenea, Diplostomoidea) illuminated with morphological, experimental, molecular and genomic data
| ||Sean A. Lockea,Fabiana B. Dragob,Danimar López-Hernándezc,Fred D. Chibwanad,Verónica Núñezb,Alex Van Dama,María Fernanda Achinellye,Pieter T.J. Johnsonf,Jordana Costa Alves de Assisc,Alan Lane de Meloc,Hudson Alves Pintoc |
| ||International Journal for Parasitology. 2021; |
|[Pubmed] | [DOI]|