Genetic variability among Trichuris ovis isolates from different hosts in Guangdong Province, China revealed by sequences of three mitochondrial genes

Comparison of the first time detected Oesophagostomum asperum with Oesophagostomum columbianum in sheep and goats in Bangladesh based on the trinity: Morphology, morphometry and genetic diversity

Oesophagostomum spp. (Family: Chabertiidae) is keeping a low profile in terms of severity in Bangladesh while maintaining economic loss through disguise within sheep and goats. The study was performed to identify prevalence, confirmation of species through morphology and morphometry followed by phylogeny using ITS2 and COX1 genes. In total 384 slaughterhouse-sourced small and large intestines were pooled from Mymensingh, Kishoreganj, Netrokona, Sherpur and Tangail districts of Mymensingh division. Followed by isolation, O. columbianum and O. asperum were identified following their key morphological features. Notably, O. asperum was first time detected in Bangladesh. The overall prevalence of Oesophagostomum spp. was found 60.93%. The prevalence of O. columbianum (64.95%) was almost double than that of O. asperum (35.04%). Among several characters, only the distance between anus to tail tip showed a significant morphological disparity in female. The Neighbor-joining (NJ) phylogenic trees based on ITS2 and COX1 genes confirmed the study species. The first time identified O. asperum along with morphometry and phylogeny will add value to the fact that nematodes are invisibly present with high prevalence in this country. This study will help to draw specific attention to command a practical control strategy for intervening in economic loss.

Bacteria associated with ovine gut parasites Trichuris ovis and Haemonchus contortus

An associated microbiome of any host helps it in different metabolic processes ranging from the decomposition of food to the maturation of gametes. Organisms with a parasitic mode of life, though present at nutritious sites inside their host, maintain their own microbiome. Nevertheless, the comprehensive characterization and functionality of microbiome in parasitic organisms remain understudied. We selected two nematode parasites of Kashmir Merino sheep viz;Haemonchus contortus and Trichuris ovis based on their higher prevalence, difference in mode of nutrition, habitation site and effect on host. The objective of the study was to explore the bacteria associated with these parasitic nematodes of sheep. We adopted a 16S rRNA metagenomic sequencing approach to estimate and compare the bacterial communities present in these two nematode species. Nematode parasites from Kashmir Merino sheep were identified morphologically and confirmed with DNA characterization. H. contortus was dominated by phylum Proteobacteria (57%), Firmicutes (25%), Bacteroidota (15%) and Actinobacteriota (3%). Conversely, T. ovis showed Proteobacteria (78%) followed by Firmicutes (8%), Bacteroidota (8%), Actinobacteriota (1%), Fusobacteriota (1%) and other phyla (4%). This study provides a comprehensive account of the microbiome composition of H. contortus and T. ovis, both of which are highly prevalent among Kashmir Merino sheep. Additionally, T. ovis exhibited a greater bacterial diversity compared to H. contortus. Notably, these nematodes were found to harbor certain pathogenic bacteria. This study can further be carried forward in gaining insights into the complex relationship between the microbiota of a parasite and its pathogenicity, reproductive potential and host microbiome modification.

Morphological and molecular characterization of Aspiculuris tetraptera (nematoda: Heteroxynematidae) from Mus musculus (rodentia: Muridae) in Saudi Arabia

Aspiculuris tetraptera a pinworm of mice, is an important parasite in institutions with mice colonies for both research and teaching purposes. Infection with this parasite has impact on biomedical research. This is likely due to the availability of the parasite’s eggs in the environment, therefore can easily be transmitted and infection is generally asymptomatic. No information regarding the prevalence, morphology or phylogeny is available on A. tetraptera from Saudi Arabia. A group of 50 laboratory mice were investigated for the presence of A. tetraptera. Worms were described morphologically and molecular characterization was attempted using 18S rRNA and Cytochrome Oxidase Subunit I genes. The prevalence of A. tetraptera infestation in the laboratory mice examined was found to be 46%. Morphological description indicated that the worms belong to A. tetraptera and this was confirmed by molecular characterization. Both regions studied have shown that the worm under investigation grouped with A. tetraptera. 18S rDNA sequences obtained in the present study showed high identity with sequences from A. tetraptera while Cytochrome c Oxidase subunit I gene (COI) sequences showed intraspecific variation resulted into two haplotypes from the isolates in the present study. A. tetraptera was recorded for the first time from Saudi Arabia. Molecular characterization has shown, based on the COI sequences, that the Saudi isolates of A. tetraptera are distinct.

Development of a Droplet Digital PCR for Detection of Trichuriasis in Sheep

Trichuriasis is a serious threat to the economic development of animal husbandry. This research aimed to establish a droplet digital PCR (ddPCR) method to detect Trichuris spp. for the early diagnosis and prevention of trichuriasis in sheep. The real-time quantitative PCR (qPCR) and ddPCR methods were used for the detection of nematodes by targeted amplification of the ITS gene. Each means was evaluated to optimize the limit of detection and reproducibility. For a recombinant plasmid, the qPCR results showed that the detection limit was 31.7 copies per reaction. In contrast to qPCR, ddPCR was able to detect concentrations below 3.17 copies per reaction. Both assays exhibited good reproducibility. However, the ddPCR method was more stable for low-copy-number detection. This new assay was specific for Trichuris spp. and did not cross-react with other relevant gastrointestinal nematodes. A total of 98 clinical samples were tested with both assays. The results showed that the positive rate of ddPCR (80.6%) was higher than that of qPCR (72.4%). This method could be used as an efficient molecular biology tool to test for Trichuris spp. and could be a new valuable tool for the clinical diagnosis and prevention of trichuriasis.

Sequencing and phylogenetic analysis of mitochondrial cox1 and nad1 genes in Toxocara canis and Toxascaris leonina isolates from Iran

Toxocara canis and Toxascaris leonina are the most important ascaridoid nematodes of the family Toxocaridae. The present study was aimed to characterisation and analysis of genetic variation within and among T. canis and T. leonina isolates obtained from Iran by sequencing partial mitochondrial cytochrome c oxidase subunit 1 (pcox1) and partial NADH dehydrogenase subunit 1 (pnad1) genes. A total number of 134 adult nematodes belonging to Toxocaridae family were collected from stray dogs in Alborz province, Iran during 2015 and 2016. Polymerase chain reaction (PCR) was performed and products were sequenced. Sequences of two mitochondrial cox1 and nad1 genes were compared with other sequences in the GenBank, while multiple sequences alignment analysis was performed using the Bioedit and MEGA6 software and phylogenetic tree was plotted. For all isolates, amplicons of about 450 and 350 base pairs (bp) were successfully produced by PCR for cox1 and nad1, respectively. All sequences of T. canis isolates from present study were 100% homologous across the nad1 gene but not in the cox1 gene. The results indicate that the PCR method based on sequence of cox1 and nad1 genes is a suitable technique for the differentiation of T. canis and T. leonina species and that mtDNA regions could be used as genetic markers for the identification and differentiation of Toxocara species.

Molecular characterization of Trichuris spp. from captive animals based on mitochondrial markers

Monoxenous parasites may easily infect animals in captivity, and nematodes belonging to the genus Trichuris are commonly reported in zoological gardens worldwide. Infections in captive animals should be accurately monitored and the characterization of pathogens is highly advisable, as a tool to infer possible routes of intra- and interspecific transmission pathways and to assess the related zoonotic potential. Whipworms are usually identified on the basis of few morphological features of adults males and eggs and by an host-affiliation criterion. Given the strong morphological convergence of adaptive traits and the possible occurrence of hybridization and/or cross-infections events, the use of molecular methods is of great utility. Here, we analysed two partial mitochondrial loci, the cytochrome c oxidase I and the cytochrome b regions, in Trichuris spp. infecting four animal species hosted in the Bioparco Zoological Garden of Rome. Results from molecular systematics, compared to previous data, suggested that the five Trichuris taxa recovered were well separated, showing a significant degree of host affiliation (herbivorous, primates/swine and rodents/canids). The screw horn antelopes and the camels were infected with two variants of Trichuris ovis; Trichuris sp. similar to those observed in rodents from South America was infecting the Patagonian maras. Moreover, Trichuris from the ring-tailed lemur showed a great similarity to Trichuris infecting the Japanese macaque previously analysed from the same zoological garden, and clustering together with Trichuris trichiura, posing a potential zoonotic threat for visitors and workers.

Reliable molecular differentiation of Trichuris ovis and Trichuris discolor from sheep (Ovis orientalis aries) and roe deer (Capreolus capreolus) and morphological characterisation of their females: morphology does not work sufficiently

The main aim of the study was to evaluate associations between morphological variability of Trichuris females from sheep and roe deer and their rDNA polymorphism in whipworm populations from the Czech Republic. The results introduced the use of new molecular markers based on the internal transcribed spacer (ITS)1-5.8S RNA-ITS2 region polymorphisms, as useful tools for the unambiguous differentiation of congeners Trichuris ovis and Trichuris discolor. These markers revealed both parasites in roe deer and in sheep; however, T. ovis females predominated in sheep while T. discolor females occurred mostly in roe deer. Additional analysis of ITS1-5.8 rRNA-ITS2 discovered the genetic uniformity of the analysed T. discolor but high haplotype variation of T. ovis. Simultaneously, molecularly designated female individuals of both species were categorised into four morphotypes (MT) on the basis of morphology of genital pore area. MT1 and MT4 (vulvar opening on everted vaginal appendage/on visible cuticular bulge) occurred only in T. ovis, MT2 (uneverted vagina-vulvar opening without any elevation) was identified only in T. discolor and MT3 (transient type of vulvar opening on a small swelling) was observed in both species. Statistical analysis of biometric data confirmed that morphology of vulva is not a reliable marker for the species determination. On the basis of the ITS1-5.8S RNA-ITS2 region variability, we carried out a phylogenetic analysis (maximum likelihood method, Hasegawa-Kishino-Yano model) which showed that T. ovis haplotypes from the Czech Republic and Ireland and T. discolor haplotypes from the Czech Republic, Spain, Iran and Japan are sister OTUs.

Nuclear and mitochondrial genes for inferring Trichuris phylogeny

Nucleotide sequences of the triose phosphate isomerase (TPI) gene (624 bp) and mitochondrial cytochrome b (cob) gene (520 bp) were obtained by PCR and evaluated for utility in inferring the phylogenetic relationships among Trichuris species. Published sequences of one other nuclear gene (18S or SSU rRNA, 1816-1846 bp) and one additional mitochondrial (mtDNA) gene (cytochrome oxidase 1, cox1, 342 bp) were also analyzed. Maximum likelihood and Bayesian inference methods were used to infer phylogenies for each gene separately but also for the combined mitochondrial data (two genes), the combined nuclear data (two genes), and the total evidence (four gene) dataset. Few Trichuris clades were uniformly resolved across separate analyses of individual genes. For the mtDNA, the cob gene trees had greater phylogenetic resolution and tended to have higher support values than the cox1 analyses. For nuclear genes, the SSU gene trees had slightly greater resolution and support values than the TPI analyses, but TPI was the only gene with reliable support for the deepest nodes in the tree. Combined analyses of genes yielded strongly supported clades in most cases, with the exception of the relationship among Trichuris clades 1, 2, and 3, which showed conflicting results between nuclear and mitochondrial genes. Both the TPI and cob genes proved valuable for inferring Trichuris relationships, with greatest resolution and support values achieved through combined analysis of multiple genes. Based on the phylogeny of the combined analysis of nuclear and mitochondrial genes, parsimony mapping of definitive host utilization depicts artiodactyls as the ancestral hosts for these Trichuris, with host-shifts into primates, rodents, and Carnivora.

Taxonomy and phylogeny of Trichuris globulosa Von Linstow, 1901 from camels. A review of Trichuris species parasitizing herbivorous

At the present work, we carried out a morph-biometrical and molecular study of Trichuris species isolated from Camelus dromedarius from Iran and from Ovis aries from South Africa comparatively with other species of Trichuris from different herbivorous hosts and geographical regions. The population from camels from Iran was identified as Trichuris globulosa. Two different morphometrically populations of Trichuris sp. from sheep from South Africa were identified: Trichuris ovis and Trichuris skrjabini. Ribosomal data did not reveal significate differences in the ITS2 sequences between T. ovis and T. globulosa to assess a specific determination. The mitochondrial data suggest that T. globulosa constitute a different genetic lineage to T. ovis. Cytochrome c-oxidase and cytochrome b partial gene sequences corroborated the existence of a different genetic lineage of T. ovis from sheep of South Africa that would be closely related to the populations of T. globulosa from camels from Iran. The cytochrome c-oxidase and cytochrome b partial gene sequences of T. globulosa have been reported for the first time.

Absence of genetic structure in Baylisascaris schroederi populations, a giant panda parasite, determined by mitochondrial sequencing

Background

Infection with the parasitic nematode, Baylisascaris schroederi (Ascaridida: Nematoda), is one of the most important causes of death in giant pandas, and was responsible for half of deaths between 2001 and 2005. Mitochondrial (mt) DNA sequences of parasites can unveil their genetic diversity and depict their likely dynamic evolution and therefore may provide insights into parasite survival and responses to host changes, as well as parasite control.

Methods

Based on previous studies, the present study further annotated the genetic variability and structure of B. schroederi populations by combining two different mtDNA markers, ATPase subunit 6 (atp6) and cytochrome c oxidase subunit I (cox1). Both sequences were completely amplified and genetically analyzed among 57 B. schroederi isolates, which were individually collected from ten geographical regions located in three important giant panda habitats in China (Minshan, Qionglai and Qinling mountain ranges).

Results

For the DNA dataset, we identified 20 haplotypes of atp6, 24 haplotypes of cox1, and 39 haplotypes of atp6 + cox1. Further haplotype network and phylogenetic analyses demonstrated that B. schroederi populations were predominantly driven by three common haplotypes, atp6 A1, cox1 C10, and atp6 + cox1 H11. However, due to low rates of gene differentiation between the three populations, both the atp6 and cox1 genes appeared not to be significantly associated with geographical divisions. In addition, high gene flow was detected among the B. schroederi populations, consistent with previous studies, suggesting that this parasite may be essentially homogenous across endemic areas. Finally, neutrality tests and mismatch analysis indicated that B. schroederi had undergone earlier demographic expansion.

Conclusions

These results confirmed that B. schroederi populations do not follow a pattern of isolation by distance, further revealing the possible existence of physical connections before geographic separation. This study should also contribute to an improved understanding of the population genetics and evolutionary biology of B. schroederi and assist in the control of baylisascariasis in giant pandas.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-014-0606-3) contains supplementary material, which is available to authorized users.

Sequence variation in two mitochondrial DNA regions and internal transcribed spacer among isolates of the nematode Oesophagostomum asperum originating from goats in Hunan Province, China

The present study examined sequence variability in two mitochondrial DNA (mtDNA) regions, namely cytochrome c oxidase subunit 1 (cox1) and NADH dehydrogenase subunit 1 (nad1), and internal transcribed spacer (ITS) of nuclear ribosomal DNA (rDNA) among Oesophagostomum asperum isolates from goats in Hunan Province, China. A portion of the cox1 (pcox1), nad1 (pnad1) genes and the ITS (ITS1+5.8S rDNA+ITS2) rDNA were amplified by polymerase chain reaction (PCR) separately from adult O. asperum individuals and the representative amplicons were subjected to sequencing from both directions. The lengths of pcox1, pnad1 and ITS rDNA were 366 bp, 681 bp and 785 bp, respectively. The A+T contents of gene sequences were 71.5–72% for pcox1, 73.7–74.2% for pnad1 and 58–58.8% for ITS rDNA. Intra-specific sequence variations within O. asperum were 0–1.6% for pcox1, 0–1.9% for pnad1 and 0–1.7% for ITS rDNA, while inter-specific sequence differences among members of the genus Oesophagostomum were significantly higher, being 11.1–12.5%, 13.3–17.7% and 8.5–18.6% for pcox1, pnad1 and ITS rDNA, respectively. Phylogenetic analyses using combined sequences of pcox1 and pnad1, with three different computational algorithms (Bayesian inference, maximum likelihood and maximum parsimony), revealed distinct groups with high statistical support. These findings demonstrated the existence of intra-specific variation in mtDNA and rDNA sequences among O. asperum isolates from goats in Hunan Province, China, and have implications for studying molecular epidemiology and population genetics of O. asperum.

Sequence variability in four mitochondrial genes among pinworm Aspicularis tetraptera isolates from laboratory mice in four provinces, China

This research aimed at exploring sequence variability in four mitochondrial (mt) genes, namely, cytochrome c oxidase subunit 1 (cox1), cytochrome b (cytb) and NADH dehydrogenase subunits 1 and 5 (nad1 and nad5), among pinworm Aspicularis tetraptera isolates from laboratory mice in four different provinces, China. A part of the cox1 (pcox1), cytb (pcytb), nad1 and nad5 genes (pnad1 and pnad5) were amplified separately from individual pinworms by polymerase chain reaction (PCR) and sequenced to determine sequence variations and examine their phylogenetic relationships. Herein, the intra-specific sequence variations within A. tetraptera were 0-0.5% for pcox1, 0-1.4% for pcytb, 0-1.8% for pnad1 and 0-1.7% for pnad5, respectively. In contrast, the inter-specific sequence differences among members of the Oxyuridae were significantly higher, being 13.7-17.0% for pcox1, 24.5-34.7% for pcytb, 26.6-29.6% for pnad1 and 24.4-25.5% for pnad5, respectively. Three methods, namely, Bayesian inference (BI), maximum likelihood (ML) and maximum parsimony (MP), were used for phylogenetic analyses based on the combined sequences of the four mt gene sequences, and the results indicated that all A. tetraptera samples form monophyletic groups, but samples from the same geographical origin did not always cluster together. This study demonstrated the existence of low-level intra-specific variation in four mtDNA sequences among A. tetraptera isolates from laboratory mice in different geographic regions in China, indicating no obvious geographical distinction among A. tetraptera isolates in China. These findings have important implications for studying systematics, molecular epidemiology and population genetics of A. tetraptera.

Genetic differences between Chabertia ovina and C. erschowi revealed by sequence analysis of four mitochondrial genes

This study examined sequence differences in mitochondrial cytochrome c oxidase subunit 1 (cox1), large subunit ribosomal RNA (rrnL) and NADH dehydrogenase subunits 1 and 4 (nad1 and nad4) between Chabertia ovina and C. erschowi from yaks in Qinghai and goats in Shaanxi provinces, China. A part of the cox1 (pcox1), rrnL (prrnL), nad1 and nad4 genes (pnad1 and pnad4) were amplified separately from individual nematodes by PCR and sequenced. The length of the sequences of pcox1, prrnL, pnad1 and pnad4 was 441 bp, 450 bp, 526 bp and 914 bp for C. ovina, and 441 bp, 451 bp, 517 bp and 810 bp for C. erschowi, respectively. The intra-specific sequence variations within C. ovina were 0.2-2.9% for pcox1, 0-0.9% for prrnL, 0.6-2.3% for pnad1, and 0.4-2.0% for pnad4, and were 0.5-1.6% for pcox1, 0-1.1% for prrnL, 0.2-1.7% for pnad1, and 0.4-1.1% for pnad4 within C. erschowi. Whereas, the inter-specific sequence differences between the two species were obviously higher, being 11.6-12.9% for pcox1, 9.8-11.1% for prrnL, 14.4-15.9% for pnad1, and 16.4-17.7% for pnad4. Phylogenetic analyses using Bayesian inference (BI), based on combined sequences of four genes, indicated that the C. ovina and C. erschowi represent distinct species. These results demonstrate that these mt gene sequences provide novel genetic markers for the identificaiton and differentiation C. ovina and C. erschowi, and have implications for studying the population genetics and molecular epidemiology of Chabertia spp.

Genetic variability in the three mitochondrial genes among Oesophagostomum asperum isolates from different regions in Shaanxi and Hunan Provinces, China

The present study examined sequence variations in three mitochondrial DNA (mtDNA) regions, namely, NADH dehydrogenase subunit 5 (nad5), adenosine triphosphate subunit 6 (atp6) and cytochrome c oxidase subunit 3 (cox3), among Oesophagostomum asperum isolates from different regions in Shaanxi and Hunan provinces, China. The lengths for partial sequences of nad5 (pnad5), atp6 (patp6) and cox3 (pcox3) were 427 bp, 381 bp and 337 bp, respectively. The intra-specific sequence variations among all O. asperum samples were 0-2.11%, 0-1.84% and 0-1.48% for pnad5, patp6 and pcox3, respectively, while the inter-specific sequence differences among Oesophagostomum species in pig and small ruminants were 18-21.3% for pnad5, 18.3-24.5% for patp6 and 10.6-13.7% for pcox3. A phylogenetic analysis based on combined sequences of three mtDNA fragments indicated that all O. asperum isolates were grouped in one solid clade, and the Oesophagostomum spp. from pig were located in another clade. However, these mtDNA fragments could not reveal genetic relationships of geographical isolates of O. asperum in China. These results provided valuable information for studying population genetics of Oesophagostomum spp., and for controlling Oesophagostomum infection in animals as well as humans.

Genetic variability of Baylisascaris schroederi from the Qinling subspecies of the giant panda in China revealed by sequences of three mitochondrial genes

The present study examined the variations in three mitochondrial (mt) DNA sequences, namely cytochrome b (cytb), cytochrome c oxidase subunit 3 (cox3) and NADH dehydrogenase subunit 5 (nad5), among Baylisascaris schroederi isolates from the Qinling subspecies of the giant panda in Shaanxi province, northwestern China. No differences in length were detected in the three mt fragments from different isolates. The intra-specific sequence variations within all B. schroederi samples were 0-2.6% for pcytb, 0-1.8% for pcox3 and 0-2.1% for pnad5, while the inter-specific sequence differences among members of the genus Baylisascaris were 8.2-15.2%, 6.2-15.9% and 8.4-16.0% for pcytb, pcox3, pnad5, respectively. A phylogenetic analysis of the combined sequences of pcytb, pcox3 and pnad 5 showed that all B. schroederi samples in the present study were located in two large clusters, with one cluster containing samples from giant pandas in Sichuan province. These findings provide basic information for further study of molecular epidemiology and control of B. schroederi infection in the Qinling subspecies of the giant panda and throughout China.

Genetic variability among Dicrocoelium dendriticum isolates from different regions in Shaanxi Province, China revealed by sequences of three mitochondrial genes

The genetic variations in three mitochondrial DNA (mtDNA) regions, namely portion of cytochrome c oxidase subunit 1 (pcox1), NADH dehydrogenase subunit 1 (pnad1) and cytochrome b (pcytb), were examined for Dicrocoelium dendriticum samples isolated from different origins in Shaanxi Province, northwestern China. The intra-specific sequence differences within D. dendriticum samples were 0-0.52% for pcox1, 0-0.73% for pnad1 and 0-0.58% for pcytb. Phylogenetic analyses based on combined sequences of three mtDNA showed that all D. dendriticum samples were clustered together in same clade of Paragonimus westermani. But the phylogenetic trees could not reveal geographically genetic relationships of D. dendriticum isolates in this province. These findings will provide basic information for further study of molecular epidemiology and control of D. dendriticum infection in this province as well as in China.