DOI: https://doi.org/10.21203/rs.3.rs-134785/v2
The objective of the present study was to identify Trichostrongylus species by molecular analysis and also phylogenetic relationships of Trichostrongylus species by mitochondrial Cox1 gene in Guilan province, northern Iran. Abomasum and duodenum contents of 144 livestock were collected from sheep, goats, and cattle in Guilan province. Morphological survey was performed for initial screening. Total DNA was extracted, and the partial region of Cox1 gene was amplified and sequenced. Genetic diversity was calculated and phylogenetic analysis of the data on nucleotide sequence was conducted by MEGA7 software. Three species of Trichostrongylus including T. colubriformis, T. vitrinus, and T. axei were identified by morphological characteristics. The genetic divergence within the species in the present study was observed for T. axei (0-2.5%), T. colubriformis (0.77%), and T. vitrinus (0%). The mean inter-species difference between the three species of Trichostrongylus obtained in this study was 14.4–15.4%. The Cox1 sequences of the members of Trichostrongylus spp. were highly variable and this could be used as a valuable measure to achieve a proper assessment on biodiversity. Sequence data generation from other species of Trichostrongylus will be needed to reconstruct the phylogenetic relationships of this genus of nematodes.
Trichostrongylus nematodes are highly prevalent and considered as gastrointestinal parasitic pathogens among ruminants with worldwide distribution (Roberts and Janovy, 2009, Sharifdini et al., 2017b). Clinical symptoms of humans are mild although in some patients gastrointestinal signs and eosinophilia may occur (Ghanbarzadeh et al., 2018, Wall et al., 2011). These nematodes are major health challenges, causing reduced animal products or even death among the infected animals in severe cases (da Rocha et al., 2020, McLeod, 1995). Several species of the parasite have been reported from herbivores with approximately 12 species identified in humans (Phosuk et al., 2013, Sharifdini et al., 2017b). Also, the frequency of Trichostrongylus spp. in human and animal hosts has been repeatedly reported in Iran (Ghadirian, 1977, Ghadirian et al., 1974, Ghasemikhah et al., 2011, Shahbazi et al., 2012). Ruminant infection was reported from Isfahan, Khuzestan, Mazandaran, Kermanshah, Hormozgan, and West Azerbaijan provinces, with human infections found in Khuzestan, Isfahan, Tehran, Hormozgan, Kermanshah, Mazandaran, Guilan, Sistan & Baluchestan, and West Azerbaijan provinces (Ashrafi et al., 2020, Ashrafi et al., 2015, Ghadirian and Arfaa, 1975, Sharifdini et al., 2020, Sharifdini et al., 2017b).
According to the morphological features reported in previous studies from Iran, several species of nematodes have been identified in human including T. orientalis, T. vitrinus, T. axei, T. colubriformis, T. probolurus, T. skrjabini, T. capricola, and T. lerouxi (Ghadirian, 1977, Ghadirian and Arfaa, 1975, Ghadirian et al., 1974). In recent years, some studies clarified the human infections with T. vitrinus, T. axei, T. colubriformis, and T. longispicularis species in endemic areas of northern Iran with T. colubriformis considered as the predominant species (Sharifdini et al., 2017b, Sharifdini et al., 2017a, Gholami et al., 2015, Ashrafi et al., 2020). Infection with various species of Trichostrongylus including T. colubriformis (Anvari-Tafti et al., 2013, Borji et al., 2010, Ghadirian and Arfaa, 1975, Ghasemikhah et al., 2011, Shahbazi et al., 2012), T. vitrinus (Anvari-Tafti et al., 2013, Borji et al., 2010, Ghadirian and Arfaa, 1975, Ghasemikhah et al., 2011, Shahbazi et al., 2012), T. axei (Ghadirian and Arfaa, 1975), T. capricola (Ghadirian and Arfaa, 1975, Ghasemikhah et al., 2011), T. probolurus (Anvari-Tafti et al., 2013, Borji et al., 2010, Ghadirian and Arfaa, 1975, Ghasemikhah et al., 2011, Shahbazi et al., 2012), T. longispicularis (Ghasemikhah et al., 2011), T. orientalis (Ghadirian and Arfaa, 1975), T. lerouxi (Biocca et al., 1974), T. skrjabini (Ghadirian and Arfaa, 1975), and T. hamatus (Anvari-Tafti et al., 2013) were reported in different herbivores such as sheep (Ghadirian and Arfaa, 1975, Ghasemikhah et al., 2011, Shahbazi et al., 2012), goats (Ghadirian and Arfaa, 1975, Ghasemikhah et al., 2011), cattle (Ghadirian and Arfaa, 1975, Ghasemikhah et al., 2011), buffalos (Ghadirian and Arfaa, 1975, Ghasemikhah et al., 2011), and camels (Anvari-Tafti et al., 2013, Borji et al., 2010) in most parts of Iran. The predominant species of Trichostrongylus among different herbivores are T. colubriformis, T. vitrinus, and T. axei found in most parts of the country (Ghadirian and Arfaa, 1975).
There is a tremendous diversity of the nematodes in the country (Ghasemikhah et al., 2011, Ghasemikhah et al., 2012) however, the molecular approaches, currently available and easily applicable, could accurately identify these species. Molecular studies based on ITS and 28S regions of ribosomal DNA were applied for genetic variation and phylogenetic analysis of Trichostrongylina (de Bellocq et al., 2001, Hoberg et al., 1999, Sharifdini et al., 2017a, von Samson-Himmelstjerna et al., 2002, Pandi et al., 2021). Although, numerous number of studies have focused on ITS2 for analysis of the Trichostrongylidae family in genetic variation, species detection, and phylogenetic relationships (Ghasemikhah et al., 2012, Sharifdini et al., 2017a, Sharifdini et al., 2017b), yet mitochondrial (mt) genomes have the potential to present valuable information. Mt genomes are conserved and present large amounts of sequence data in the organisms, therefore mtDNA are used for evolutionary analyses, taxonomy, population genetics, and systematics studies (dos Santos et al., 2017, Hu et al., 2004, Saccone et al., 1999). There are few studies that have investigated the mitochondrial gene of the Trichostrongylidae family, in which the mtDNA of Marshallagia marshalli, Haemonchus placei, Haemonchus contortus, T. vitrinus, T. axei, Ostertagia trifucata, and Teladorsagia circumcincta species were evaluated for phylogenetic relationship and species identification (Ahmad et al., 2019, Archie and Ezenwa, 2011, dos Santos et al., 2017, Jex et al., 2010, Kuchboev et al., 2020, Palevich et al., 2020, Sun et al., 2018). Taxonomy studies of the nematodes based on sequences of coding mitochondrial genes are more accurate than non-coding ribosomal genes. While mitochondrial genomes are considered as suitable markers for population evolution studies (Kuchboev et al., 2020, Palevich et al., 2020), the studies targeting the mtDNA for identification of Trichostrongylidae family are very limited worldwide with even no single report on mitochondrial gene of the nematodes from Iran. Therefore, the present study focused on molecular phylogenetic analysis based on cytochrome c oxidase subunit I (Cox1) from mitochondrial gene of Trichostrongylus species in northern Iran.
In current study, a total of 144 abomasum and duodenum specimens from livestock, including 72 cattle, 59 sheep, and 13 goats were collected from the abattoir of Talesh district in Guilan province, northern Iran during July to September 2018 (Hosseinnezhad et al., 2021) (Figure. 1).
The Trichostrongylidae family members were isolated by washing the abomasum and duodenum contents followed by passing through the 20, 40, and 100 mesh screens. The helminths captured on mesh screens were examined under stereomicroscope. Morphological features were evaluated after cleaning the worms with normal saline and lactophenol. The samples were preserved in 70% ethanol at room temperature until used (Barghandan et al., 2019).
Male parasites were isolated for DNA extraction. Total genomic DNA was extracted from one male worm of each species of trichostrongyloid nematodes collected from all study animals, using a commercial DNA extraction kit (Yekta Tajhiz Azma, Tehran, Iran) according to the manufacturer’s instructions. The partial region of the cox1 gene with approximately 700bp was amplified using the LCO1490 (5′-GGTCAACAAATCATAAAGATATTGG-3′) and HCO2198 (5′-TAAACTTCAGGGTGACCAAAAAATCA-3′) sequences as cox1 gene forward and reverse primers (Folmer, Black et al. 1994). The thermal PCR profiles included an initial denaturation step at 95°C for 6 minutes followed by 35 cycles of denaturation at 95°C for 45 seconds, annealing at 50°C for 45 seconds, an initial extension step at 72°C for 60 seconds, and a final extension step at 72 ºC for 10 minutes.
The PCR products were sequenced using an ABI 3130xl platform (Applied Biosystems, Foster City, California, USA). The sequences identified by the ABI system were edited and analyzed by BioEdit software (http://www.mbio.ncsu.edu/bioedit/bioedit.html). The sequences were compared with the sequences deposited in the GenBank database by BLAST program (http://www.ncbi.nlm.nih.gov/). The sequences of the three species of genus Trichostrongylus derived from the domestic ruminants and deposited in the GenBank database were marked by the following Accession Numbers: MW051252-MW051254 for T. axei; MW051250 and MW051251 for T. colubriformis, and MW051255 and MW051256 for T. vitrinus).
Multiple sequence alignments were conducted by ClustalW incorporated in the BioEdit software. Phylogenetic tree was constructed by the MEGA7 software (Molecular and Evolution Genetic Analysis v7). The maximum likelihood method based on the Tamura 3-parameter model and maximum–likelihood algorithm was applied. Bootstrap value was done based on 1000 replications in the topology of the tree while Dictyocaulus capreolus considered as out group (Figure. 2). Pairwise distance comparisons clarified the presence of the seven sequences of the three species of Trichostrongylus, isolated in the present study as well as other species using BioEdit software.
All of the study male worms were identified based on the morphological characteristics of male copulatory spicules and gubernaculum (Figure. 1). Trichostrongylus axei was isolated from the cattle, sheep, and goats while T. colubriformis and T. vitrinus were only detected among the sheep and goats. The isolates were successfully amplified for Cox1 gene with specific band. The sequence results confirmed three species of T. colubriformis, T. vitrinus, and T. axei among the specimens. A dendrogram, based on the phylogenetic analysis, showed that the species were placed along, with the same species obtained from the GenBank database, into a distinct cluster of the tree (Figure- 2). The genetic divergence within the species of T. axei, T. colubriformis, and T. vitrinus obtained in this study were 0-2.5%, 0.77%, and 0%, respectively. Two species of Trichostrongylus including T. axei and T. vitrinus isolated from the sheep and goats were quite similar. The intra-species distance rate within the specimens of T. axei, T. colubriformis, and T. vitrinus found in the present study and those available in the GenBank database amounted to 0.95–3.1% (1.9%), 0.19–4.08% (2.4%), and 0-2.32% (1.5%), respectively.
In this study, the mean inter-species differences between our T. axei specimens, compared with T. colubriformis and T. vitrinus isolates, were 14.4% and 14.6%, respectively. Also, the mean genetic difference between the T. colubriformis specimens was 15.4% when compared with T. vitrinus.
Based on our sequences and those deposited in the GenBank, the mean inter-species distance rates between the isolates of T. axei and those of T. colubriformis and T. vitrinus were 13.5% and 14.5%, respectively. Also, the mean genetic diversity between the isolates of T. colubriformis and those of T. vitrinus was 14.9%.
The three species of the Trichostrongylus including T. colubriformis, T. vitrinus, and T. axei identified in the present study, along with the data already reported from Iran confirm that the predominant species in herbivorous animals (Ghadirian and Arfaa, 1975). Iran is one of the most important foci for Trichostrongylus infection among human and animal hosts (Alemi and Arfaa, 1978, Ashrafi et al., 2020, Massoud et al., 1980, Sharifdini et al., 2017a). Proper conditions such as humidity and climate in the northern parts of the country including Mazandaran and Guilan provinces lead to permanent establishment of the life cycle process of soil transmitted helminthes in the regions (Alemi and Arfaa, 1978, Sharifdini et al., 2017a, Sharifdini et al., 2017b).
In the present study the authors used the sequence analysis protocol for detecting the mitochondrial Cox1 gene, whereas several other studies, reported from Iran, employed ITS-rDNA gene specific for the phylogenetic analysis of Trichostrongylus species (Ashrafi et al., 2020, Sharifdini et al., 2017a, Sharifdini et al., 2017b). The nuclear ribosomal gene is widely applied to the studies of deep and shallow phylogenetic relationship in the phylum Nematoda (Holterman et al., 2006, Kiewnick et al., 2014, Sharifdini et al., 2017a, Sharifdini et al., 2017b). Recent studies illustrated that the mitochondrial genes to be the proper options for phylogenomic approach and specifically for the CoxI gene that has mainly been used in population genetic surveys for various nematode parasites of the vertebrates (Kiewnick et al., 2014, Nadler and Hudspeth, 2000, Otranto et al., 2005).
Several studies have demonstrated that the sequence differences between the members of the Trichostrongylus are not noticeable when the detection protocol is based on the ITS2 gene (Ashrafi et al., 2020, Sharifdini et al., 2017a, Sharifdini et al., 2017b). Ashrafi et al. (2020) reported a mean inter-species distance rate of 2.6% within different species of Trichostrongylus while in the current study the mean inter-species variation within our specimens and those available in the GenBank was 13.5–14.9%. Due to the high level divergence in the Cox1 gene, it could be considered as a valuable genetic tool for phylogenetic and taxonomic studies on the members of the Trichostrongylus genus (Ghasemikhah et al., 2012). The phylogenetic tree constructed in our study represented that the three species of T. colubriformis, T. vitrinus, and T. axei were separated in distinct cluster along with the same species obtained from other studies in different countries (Figure. 2). The results of genetic diversity within the species showed that the intra-species distance rate among the present isolates was so close, indicating high proximity of the sequences in the region.
Little information on mitochondrial genes of Trichostrongyloidea superfamily is available. Palevich et al. (2020) in New Zealand investigated the complete mitochondrial genomes of H. contortus and T. circumcincta by phylogenetic analysis (Palevich et al., 2020). Another study, reported from Uzbekistan, was based on ribosomal (ITS2) and mitochondrial (Cox1) of Marshallagia sp. and concluded that the ITS2 sequences has little variation and is not a suitable gene for diagnosing different species, while Cox1 gene shows more diversities (Kuchboev et al., 2020). Ostertagia trifurcata and Marshallagia marshalli were evaluated by phylogenetic analysis of the complete mitochondrial genes in China and the findings introduced complete mt genome sequence of the nematodes as a novel genetic marker for population genetic and molecular epidemiology (Ahmad et al., 2019, Sun et al., 2018). Two other studies, reported from Brazil and Australia, evaluated the complete mitochondrial genes of H. placei, T. circumcincta, T. vitrinus, and T. axei and suggested that the phylogenomics approach of mtDNA could be applied as a new genetic marker in phylogenetic analysis and geographic relationships among different isolates in population genetic studies (dos Santos et al., 2017, Jex et al., 2010). Moreover, the Cox1 and nad4 genes of T. axei were also analyzed for population genetic structure of the nematode in USA (Archie and Ezenwa, 2011). However, additional sequence studies, especially the analysis of both nuclear and mitochondrial genes, are needed to provide a comprehensive understanding of the genetic variations of Trichostrongylus spp. in endemic areas and other parts of Iran.
In the present study three species of T. colubriformis, T. vitrinus, and T. axei were observed among the specimens of Guilan province, northern Iran. This study concluded the genetic diversity of the Cox1 gene is notable and the gene is suitable for analyzing the gene diversity of intra-species distance among helminthes. The scarcity of molecular data on Cox1 gene within Trichostrongylus spp. in various geographical regions and hosts makes it necessary to produce sufficient data on diversities of this gene which eventually leads to reconstruct the total phylogenetic relationships of this group of nematode. Thus, the findings of the present study suggest that the analysis of complete mitochondrial genome to be the focus of further experiments in the future research.
The authors declare that there is no conflict of interest. We would like to appreciate the assistance offered by the colleagues at the Department of Parasitology and Mycology, Medical School, Guilan University of Medical Sciences. This study was financially supported by Guilan University of Medical Sciences.The authors are grateful to the abattoir personnel in Talesh, Guilan Province of Iran. This study was approved by the Medical Ethics Committee of Guilan University of Medical Sciences (Approval code: IR.GUMS.REC.1397.176).