Oesophagostomum dentatum and Oesophagostomum quadrispinulatum: characterization of the complete mitochondrial genome sequences of the two pig nodule worms

The complete mitochondrial genome of Morishitium polonicum (Trematoda, Cyclocoelidae) and its phylogenetic implications

Trematodes can adversely impact the health and survival of wild animals. The trematode family Cyclocoelidae, which includes large digenean bird parasites, lacks molecular analysis, and reclassifications have not been supported. This study produced the first fully assembled and annotated mitochondrial genome sequence for the trematode Morishitium polonicum. The whole length of the M. polonicum (GenBank accession number: OP930879) mitogenome is 14083 bp, containing 22 transfer ribonucleic acids (tRNAs), 2 ribosomal RNAs (rRNAs, rrnL and rrnS), and a noncoding control section (D-loop) 13777 to 13854 bp in length. The 12 PCG areas have 3269 codons and a total length of 10053 bp, which makes up 71.38% of the mitochondrial genome's overall sequence. Most (10/12) of the PCGs that code for proteins begin with ATG, while the nad4L and nad1 genes have a GTG start codon. Phylogenetic analysis using the concatenated nucleotide sequences of 12 PCGs, and the ML tree analysis results showed that M. polonicum is more closely related to with Echinostomatidae and Fasciolidae, which indicates that the family Cyclocoelidae is more closely associated with Echinochasmidae. This study provides mtDNA information, and analysis of mitogenomic structure and evolution. Moreover, we aimed to understand the phylogenetic relationships of this fluke.

Characterization of the complete mitochondrial genome of Plagiorchis multiglandularis (Digenea, Plagiorchiidae): Comparison with the members of Xiphidiatan species and phylogenetic implications

Plagiorchis multiglandularis Semenov, 1927 is a common fluke of birds and mammals, with significant impacts on animals and also human health. However, the systematics of Plagiorchiidae remain ambiguous. In the present study, the complete mitochondrial (mt) genome of P. multiglandularis cercariae was sequenced and compared with other digeneans in the order Xiphidiata. The complete circular mt genome of P. multiglandularis was 14,228 bp in length. The mitogenome contains 12 protein-coding genes and 22 transfer RNA genes. The 3' end of nad4L overlaps the 5' end of nad4 by 40 bp, while the atp8 gene is absent. Twenty-one transfer RNA genes transcribe products with conventional cloverleaf structures, while one transfer RNA gene has unpaired D-arms. Comparative analysis with related digenean trematodes revealed that A + T content of mt genome of P. multiglandularis was significantly higher among all the xiphidiatan trematodes. Phylogenetic analyses demonstrated that Plagiorchiidae formed a monophyletic branch, in which Plagiorchiidae are more closely related to Paragonimidae than Prosthogonimidae. Our data enhanced the Plagiorchis mt genome database and provides molecular resources for further studies of Plagiorchiidae taxonomy, population genetics and systematics.

Comprehensive Molecular Characterization of the Mitochondrial Genome of the Takin Lungworm Varestrongylus eleguneniensis (Strongylida: Protostrongylidae)

The takin lungworm Varestrongylus eleguneniensis (Strongylida: Protostrongylidae) causes lethal bronchopneumonia and represents severe threats to captive and wild populations. However, until now there has been very limited information available concerning the molecular epidemiology and evolutionary biology of V. eleguneniensis. Mitochondrial genomes (mtDNAs) can provide resources for investigations in these areas and, therefore, can assist with the surveillance and control of this lungworm. Herein, the complete mtDNA of V. eleguneniensis was sequenced and characterized with Illumina pipeline analyses. This circular genome (13,625 bp) encoded twelve protein-coding genes (PCGs), two rRNAs, and twenty-two tRNAs, with notable levels of AT and GC skews. Comparative genomics revealed a purifying selection among PCGs, with cox1 and nad6 having the lowest and the highest evolutionary rate, respectively. Genome-wide phylogenies showed a close relationship between V. eleguneniensis and Protostrongylus rufescens in Strongylida. Single gene (PCGs or rRNAs)-based phylogenies indicated that cox1 and nad5 genes shared the same family-level topology with that inferred from genomic datasets, suggesting that both genes could be suitable genetic markers for evolutionary and phylogenetic studies of Strongylida species. This was the first mtDNA of any member of the genus Varestrongylus, and its comprehensive molecular characterization represents a new resource for systematic, population genetic and evolutionary biological studies of Varestrongylus lungworms in wildlife.

Phylogenetic Relationships of the Strongyloid Nematodes of Australasian Marsupials Based on Mitochondrial Protein Sequences

Simple Summary

Parasitic strongyloid nematodes endemic to the gastrointestinal tracts of Australasian marsupials are one of the most diverse groups of mammalian parasites. These nematodes are currently placed in the family Chabertiidae comprising two subfamilies, namely the Cloacininae and Phascolostrongylinae. Their current classification relies primarily on morphological features and has not been validated using molecular data. This study aimed to determine the phylogenetic relationships of the Cloacininae and Phascolostrongylinae within the family Chabertiidae and their relationship with other groups of strongyloid nematodes from non-marsupial hosts, using mitochondrial protein sequence datasets. The findings supported the recognition of the family Cloacinidae, containing the Cloacininae and Phascolostrongylinae, as a monophyletic group within the Strongyloidea. However, the subfamily Phascolostrongylinae was paraphyletic, and the relationships of individual genera corresponded with their host families. Genera of the Cloacininae and Phascolostrongylinae occurring in macropod hosts were more closely related compared to genera of the Phascolostrongylinae occurring in wombats. This study suggests an alternative hypothesis for the origin of marsupial strongyloid nematodes in vombatid hosts that should be explored further using molecular approaches and more widespread sampling.

Abstract

Australasian marsupials harbour a diverse group of gastrointestinal strongyloid nematodes. These nematodes are currently grouped into two subfamilies, namely the Cloacininae and Phascolostrongylinae. Based on morphological criteria, the Cloacininae and Phascolostrongylinae were defined as monophyletic and placed in the family Cloacinidae, but this has not been supported by molecular data and they are currently placed in the Chabertiidae. Although molecular data (internal transcribed spacers of the nuclear ribosomal RNA genes or mitochondrial protein-coding genes) have been used to verify morphological classifications within the Cloacininae and Phascolostrongylinae, the phylogenetic relationships between the subfamilies have not been rigorously tested. This study determined the phylogenetic relationships of the subfamilies Cloacininae and Phascolostrongylinae using amino acid sequences conceptually translated from the twelve concatenated mitochondrial protein-coding genes. The findings demonstrated that the Cloacininae and Phascolostrongylinae formed a well-supported monophyletic assemblage, consistent with their morphological classification as an independent family, Cloacinidae. Unexpectedly, however, the subfamily Phascolostrongylinae was split into two groups comprising the genera from macropodid hosts (kangaroos and wallabies) and those from vombatid hosts (wombats). Genera of the Cloacininae and Phascolostrongylinae occurring in macropodid hosts were more closely related compared to genera of the Phascolostrongylinae occurring in wombats that formed a sister relationship with the remaining genera from macropods. These findings provide molecular evidence supporting the monophyly of the family Cloacinidae and an alternative hypothesis for the origin of marsupial strongyloid nematodes in vombatid hosts that requires further exploration using molecular approaches and additional samples

Cytochrome c oxidase subunit 1 gene reveals species composition and phylogenetic relationships of Oesophagostomum spp. infecting pigs in northeastern Brazil

Helminths of the genus Oesophagostomum cause enteric diseases and affect domestic animals such as pigs. The aim of this study was to explore the species composition and genetic diversity of Oesophagostomum spp. infecting pigs in close contact with humans in the state of Piauí, Brazil. Eighty-seven fecal samples were collected for parasitological tests and molecular analysis. Through microscopy, the overall positivity rate for strongyliform eggs was 81.6% among the pigs studied. Forty-two strongyliform egg samples were subjected to PCR and six cox1 sequences (637 bp) were identified for the genus Oesophagostomum. The sequences were identified as Oesophagostomum dentatum, O. quadrispinulatum and O. columbianum. In the phylogenetic tree and haplotype network, 89 sequences were separated into seven clusters, which also included reference sequences from GenBank. Oesophagostomum dentatum and O. quadrispinulatum were seen to be closely related species and formed a monophyletic group related to O. aculeatum. Oesophagostomum columbianum showed similarity with sequences from parasites infecting small ruminants and the clade was positioned closer to O. bifurcum. High interspecific diversity was found and intraspecific diversity varied according to the species. This was the first study to characterize Oesophagostomum DNA sequences obtained from pigs in Brazil.

Phylogenetic relationships of the nematode subfamily Phascolostrongylinae from macropodid and vombatid marsupials inferred using mitochondrial protein sequence data

BACKGROUND

The subfamily Phascolostrongylinae (Superfamily Strongyloidea) comprises nematodes that are parasitic in the gastrointestinal tracts of macropodid (Family Macropodidae) and vombatid (Family Vombatidae) marsupials. Currently, nine genera and 20 species have been attributed to the subfamily Phascolostrongylinae. Previous studies using sequence data sets for the internal transcribed spacers (ITS) of nuclear ribosomal DNA showed conflicting topologies between the Phascolostrongylinae and related subfamilies. Therefore, the aim of this study was to validate the phylogenetic relationships within the Phascolostrongylinae and its relationship with the families Chabertiidae and Strongylidae using mitochondrial amino acid sequences.

METHODS

The sequences of all 12 mitochondrial protein-coding genes were obtained by next-generation sequencing of individual adult nematodes (n = 8) representing members of the Phascolostrongylinae. These sequences were conceptually translated and the phylogenetic relationships within the Phascolostrongylinae and its relationship with the families Chabertiidae and Strongylidae were inferred from aligned, concatenated amino acid sequence data sets.

RESULTS

Within the Phascolostrongylinae, the wombat-specific genera grouped separately from the genera occurring in macropods. Two of the phascolostrongyline tribes were monophyletic, including Phascolostrongylinea and Hypodontinea, whereas the tribe Macropostrongyloidinea was paraphyletic. The tribe Phascolostrongylinea occurring in wombats was closely related to Oesophagostomum spp., also from the family Chabertiidae, which formed a sister relationship with the Phascolostrongylinae.

CONCLUSION

The current phylogenetic relationship within the subfamily Phascolostrongylinae supports findings from a previous study based on ITS sequence data. This study contributes also to the understanding of the phylogenetic position of the subfamily Phascolostrongylinae within the Chabertiidae. Future studies investigating the relationships between the Phascolostrongylinae and Cloacininae from macropodid marsupials may advance our knowledge of the phylogeny of strongyloid nematodes in marsupials.

The mitochondrial genome sequence analysis of Ophidascaris baylisi from the Burmese python (Python molurus bivittatus)

Ophidascaris species are parasitic roundworms that inhabit the python gut, resulting in severe granulomatous lesions or even death. However, the classification and nomenclature of these roundworms are still controversial. Our study aims to identify a snake roundworm from the Burmese python (Python molurus bivittatus) and analyze the mitochondrial genome. We identified this roundworm as Ophidascaris baylisi based on the morphology and cytochrome c oxidase subunit I (cox1) sequence. Ophidascaris baylisi complete mitochondrial genome was 14,784 bp in length, consisting of two non-coding regions and 36 mitochondrial genes (12 protein-coding genes, 22 tRNA genes, and two rRNA genes). The protein-coding genes used TTG, ATG, ATT, or TTA as start codons and TAG, TAA, or T as stop codons. All tRNA genes showed a TV-loop structure, except trnS1^AGN and trnS2^UCN revealed a D-loop structure. The mitochondrial large ribosomal subunit 16S (rrnL) and small ribosomal subunit 12S (rrnS) were 956 bp and 700 bp long, respectively. Phylogenetic analysis based on O. baylisi mitochondrial protein-coding genes demonstrated that O. baylisi clustered with the family Ascarididae members and was most closely related to Ophidascaris wangi. These results may enhance the nematode mitochondrial genome database and provide valuable molecular markers for further research on the taxonomy, phylogeny, and genetic relationships of Ophidascaris nematodes.

Comparative analysis of mitochondrial DNA datasets indicates that Cylicostephanus minutus represents a species complex

Cyathostomins are one kind of the most important parasites in equids. Cylicostephanus minutus is a member of the subfamily Cyathostominae. In the present study, we determined the complete mitochondrial (mt) genomes from four Cs. minutus isolates and reconstructed the phylogenetic relationship of Strongylidae to test the hypothesis that Cs. minutus represents a species complex. The complete mt genome sequences of Cs. minutus were 13,772-13,822 bp in length, and contained 36 genes (12 protein coding genes, 22 tRNA genes, two rRNA genes), and two non-coding regions (NCRs). The intraspecific identity of nucleotide sequences and amino acid sequences in Cs. minutus (1-4) were 89.3-97.9% and 97.0-98.8%, respectively. Two operational taxonomic units (OTUs) were determined based on the mt genome sequences, OTU 2 (Csm 1 and Csm 2) and OTU 3 (Csm 3 and Csm 4). Sequence analysis showed the divergence between OTU 2 and OTU 3 was 8.9-10.7%. Pairwise comparisons of 12 protein coding genes between OTU 2 and OTU 3 showed a difference of 3.0-13.3% at the nucleotide level and 0-6.7% at the amino acid level. Phylogenetic analysis showed the separation of Cs. minutus isolates from the same host into different distinct clades based on mt genomes. Comparisons of partial mt cox1, nad5, and cytb and ITS2 sequences from 20 Cs. minutus isolates from the same host and the same geographical location with other Cs. minutus sequences available in GenBank revealed significant nucleotide differences. Phylogenetic analysis showed a separation of Cs. minutus into three distinct clades. Thus, the comparative and phylogenetic analyses of mtDNA datasets indicated that Cs. minutus represents a complex of at least three species. Our results have further confirmed the existence of a cryptic Cs. minutus species, and provides a reference for the taxonomical, population genetics, and systematics studies of other cyathostomin species.

Mitochondrial genome evidence suggests Cooperia sp. from China may represent a distinct species from Cooperia oncophora from Australia

Cooperia spp. are parasitic nematodes parasitizing in small intestine of ruminants with a worldwide distribution. Infection of ruminants with Cooperia species can cause severe enteritis, causing significant socio-economic losses to the livestock industry. However, it is yet to know whether there is genetic diversity in mitochondrial (mt) DNA sequences of Cooperia nematodes from different geographic regions. The objective of the present study was to examine sequence difference in mt genomes between Cooperia sp. from China and other Cooperia species. We determined the sequences of the internal transcribed spacer (ITS-1 and ITS-2) of nuclear ribosomal DNA (rDNA) of 11 Cooperia specimens collected from the small intestine of a Tianzhu White yak in Gansu Province, northwestern China, which had 99% similarity with that of C. oncophora from Brazil (GenBank accession Number: AJ544290) in ITS-1, and 99% similarity with those from Denmark (AB245040), Scotland and Australia (AJ000032) in ITS-2, indicating that specimens used in the present study should at least represent parasites in Cooperia. We then determined the complete mt genome sequences of one representative specimen of Cooperia sp. from China (CspC), compared the mt DNA sequences with that of C. oncophora from Australia (COA, GQ888713), and conducted phylogenetic analysis with selected nematodes using both maximum likelihood (ML) and Bayesian inference (BI) methods based on both concatenated 12 PCGs, rrnL and rrnS sequences and partial cox2 sequences. The complete mt genome sequence of CspC (KY769271) is 13, 583 bp in length, which is 91 bp shorter than that from COA. The sequence difference over the entire mt genome between CspC and COA was 12.2% in nucleotide and 6.3% in inferred amino acids, with nad4L and nad1 being the most variable and the most conserved PCGs, respectively. Phylogenetic analysis indicated that CspC and COA were closely-related but distinct taxa. The determination of mt genome sequences for Cooperia sp. from China also provides novel resources for further studies of taxonomy, systematics and population genetics of Cooperia from different geographical locations.

Comparative analysis of mitochondrial DNA datasets indicates that Toxascaris leonina represents a species complex

Background

Toxascaris leonina is one of the most common intestinal parasites of canids and felids. In this study, we characterised the entire mitochondrial (mt) genome sequence of T. leonina from the cheetah and compared it with that of T. leonina from the dog.

Results

The entire mt genome sequence of T. leonina from the cheetah is 14,685 bp in size, which is 375 bp longer than that from the dog, and it is 408 bp longer than that from the South China tiger. The overall nucleotide sequence (except for the non-coding region) identity was 92.8% between the two mt genomes of T. leonina from the cheetah and the dog. For the 12 protein-coding genes, sequence difference between T. leonina from the cheetah and the dog was 5.0–9.7% at the nucleotide level and 1.0–7.2% at the amino acid level. Moreover, comparison of mt cox1 sequences among T. leonina isolates (n = 23) from different hosts revealed substantial nucleotide differences (10.6%). Phylogenetic analysis showed the separation of T. leonina from canid and felid hosts into three distinct clades.

Conclusions

Taken together, these mtDNA datasets indicate that T. leonina from canid and felid hosts represents a species complex. Our results have implications for further studies of the molecular epidemiology, systematics and population genetics of this nematode.

Electronic supplementary material

The online version of this article (10.1186/s13071-019-3447-2) contains supplementary material, which is available to authorized users.

Characterization of the Complete Mitochondrial Genome of Ostertagia trifurcata of Small Ruminants and its Phylogenetic Associations for the Trichostrongyloidea Superfamily

The complete mitochondrial (mt) genome of Ostertagia trifurcata, a parasitic nematode of small ruminants, has been sequenced and its phylogenetic relationship with selected members from the superfamily Trichostrongyloidea was investigated on the basis of deduced datasets of mt amino acid sequences. The entire mt genome of Ostertagia trifurcata is circular and 14,151 bp in length. It consists of a total of 36 genes comprising 12 genes coding for proteins (PCGs), 2 genes for ribosomal RNA (rRNA), 22 transfer RNA (tRNA) genes and 2 non-coding regions, since all genes are transcribed in the same direction. The phylogenetic analysis based on the concatenated datasets of predicted amino acid sequences of the 12 protein coding genes supported monophylies of the Haemonchidae, Dictyocaulidae and Molineidae families, but rejected monophylies of the Trichostrongylidae family. The complete characterization and provision of the mtDNA sequence of Ostertagia trifurcata provides novel genetic markers for molecular epidemiological investigations, systematics, diagnostics and population genetics of Ostertagia trifurcata and its correspondents.

Characterization of the complete mitochondrial genome of Metastrongylus salmi (M. salmi) derived from Tibetan pigs in Tibet, China

The present study was designed to determine and analyze the mt genomes of Metastrongylus salmi (M. salmi), and reveal the phylogenetic relationships of this parasite using mt DNA sequences. Results showed that the complete mt genome of M. salmi was 13722 bp containing 12 protein-coding genes (cox1-3, nad1-6, nad4L, atp6 and cytb), 22 transfer RNA genes, and 2 ribosomal RNA genes (rrnL and rrnS). The overall A+T content was 73.54% and the nucleotide composition was A (23.52%), C (6.14%), G (19.60%), T (50.02%), and N (UCAG) (0.73%). A total of 4237 amino acids are encoded from the Tibetan isolates of M. salmi mt genomes. The ATA was predicted as the most common starting codon with 41.7% (5/12 protein genes); and 11 of the 12 protein genes were found to have a TAG or TAA translation termination codon. By clustering together the phylogenetic trees of Tibetan M. salmi and Austrian M. salmi, the M. salmi isolated from Tibetan pigs was found to be highly homological with that stemmed from Austrian one. This information provides meaningful insights into the phylogenetic position of the M. salmi China isolate and represents a useful resource for selecting molecular markers for diagnosis and population studies.

According to mitochondrial DNA evidence, Parascaris equorum and Parascaris univalens may represent the same species

Parascarosis is caused mainly by parasitic infections with Parascaris equorum and Parascaris univalens, the most common ascarid nematodes, in the small intestine of equines. Parascarosis often causes severe illness and even death in foals and yearlings. In this study, we obtained the complete sequence of the P. equorum mitochondrial (mt) genome and compared its organization and structure with that of P. equorum Japan isolate (nearly complete), and the complete mtDNA sequences of P. univalens Switzerland and USA isolates. The complete mtDNA genome of P. equorum China isolate is 13,899 base pairs (bp), making it the smallest of the four genomes. All four Parascaris mt genomes are circular, and all genes are transcribed in the same direction. The P. equorum mtDNA genome consists of 12 protein-coding genes, two ribosomal RNA genes, 22 transfer (t) RNA genes and one non-coding region, which is consistent with P. equorum Japan isolate and P. univalens Switzerland isolate but distinct from P. univalens USA isolate, which has 20 tRNA genes. Differences in nucleotide sequences of the four entire mt genomes range from 0.1–0.9%, and differences in total amino acid sequences of protein-coding genes are 0.2–2.1%. Phylogenetic analyses showed that the four Parascaris species clustered in a clade, indicating that P. equorum and P. univalens are very closely related. These mt genome datasets provide genetic evidence that P. equorum and P. univalens may represent the same species, which will be of use in further studies of the taxonomy, systematics and population genetics of ascarids and other nematodes.

Characterization of the complete mitochondrial genome of Marshallagia marshalli and phylogenetic implications for the superfamily Trichostrongyloidea

Marshallagia marshalli (Nematoda: Trichostrongylidae) infection can lead to serious parasitic gastroenteritis in sheep, goat, and wild ruminant, causing significant socioeconomic losses worldwide. Up to now, the study concerning the molecular biology of M. marshalli is limited. Herein, we sequenced the complete mitochondrial (mt) genome of M. marshalli and examined its phylogenetic relationship with selected members of the superfamily Trichostrongyloidea using Bayesian inference (BI) based on concatenated mt amino acid sequence datasets. The complete mt genome sequence of M. marshalli is 13,891 bp, including 12 protein-coding genes, 22 transfer RNA genes, and 2 ribosomal RNA genes. All protein-coding genes are transcribed in the same direction. Phylogenetic analyses based on concatenated amino acid sequences of the 12 protein-coding genes supported the monophylies of the families Haemonchidae, Molineidae, and Dictyocaulidae with strong statistical support, but rejected the monophyly of the family Trichostrongylidae. The determination of the complete mt genome sequence of M. marshalli provides novel genetic markers for studying the systematics, population genetics, and molecular epidemiology of M. marshalli and its congeners.

Prevalence of three Oesophagostomum spp. from Tibetan Pigs analyzed by Genetic Markers of nad1, cox3 and ITS1

Parasitic nematodes of Oesophagostomum spp., commonly known as 'nodular worms' are one of the most widely distributed and prevalent emerging zoonotic nematodes. However, little is known about the prevalence and gene characteristics of those parasites in Tibetan pigs. Therefore, a study was carried out to investigate the prevalence, isolation and identification of Oesophagostomum spp from Tibetan pigs by genetic markers of nad1,cox3 and ITS1 for the first time. The results revealed that the infection rate of O. dentatum and O. quadrispinulatum by genetic markers of nad1 was 81.13%; 35 (66.04%); the O. dentatum and O. quadrispinulatum by genetic markers cox3 was 66.04%, and O. dentatum and O. stephanostomum by genetic markers ITS1 was found to be 77.36%. Interestingly, the O. stephanostomum specie was identified and isolated from 90.48% stomach and 69.23% colon samples by genetic markers of ITS1. The present study, for the first time has described the presence and genetic characterization of Oesophagostomun spp of O. dentatum, O. quadrispinulatum and especially O. stephanostomum in Tibetan pigs from the high and remote Tibetan plateau. A public concern should be raised in terms of economical losses and severe public health problem.

The mitochondrial genome of Ancylostoma tubaeforme from cats in China

Ancylostoma tubaeforme may infect canids, felids and humans, and pose a potential risk to public health. Polymerase chain reaction (PCR) techniques were used to amplify the complete mitochondrial (mt) genome sequence of A. tubaeforme from cats and to analyse its sequence characteristics after molecular identification based on the internal transcribed spacer ITS1+ sequence. The results show that the complete mt genome sequence (GenBank accession number KY070315) of A. tubaeforme from cats was 13,730 bp in length, including 12 protein-coding genes, 22 transfer RNA (tRNA) genes, two ribosomal RNA (rRNA) genes, two non-coding regions and an AT-rich region. The nucleotide content of A and T was 77.93%, biased toward A and T. Twelve protein-coding genes used ATT, TTG and GTG as initiation codons, and TAA, TAG, TA and T as termination codons. The length of the 22 tRNA genes ranged from 52 to 62 bp, their predicted secondary structures were D loops and V loops. The lengths of the two rRNAs were 958 and 697 bp. Phylogenetic analyses showed that A. tubaeforme from cats was in the lineage of Ancylostoma, having a close phylogenetic relationship with A. caninum. This study reports for the first time the mt genome of A. tubaeforme from cats in China, which could enhance the mt genome database of Ancylostomatidae nematodes, and it offers the scientific basis for further studies in the genetic diversity of hookworms among different hosts.

Complete mitochondrial genomes of Gnathostoma nipponicum and Gnathostoma sp., and their comparison with other Gnathostoma species

Gnathostomiasis is a foodborne zoonotic parasitosis caused by Gnathostoma nematodes. It has caused significant public problems worldwide, but its molecular biology is limited. The purpose of this study was to decode the complete mitochondrial (mt) genomes of Gnathostoma nipponicum and Gnathostoma sp., and compare their mt sequences with other Gnathostoma species. The complete mt genome sequences were amplified by long-range PCR and determined by subsequent primer walking. The complete mt genomes of G. nipponicum and Gnathostoma sp. were 14,093bp and 14,391bp, respectively. Both of the two mt genomes contain 12 protein-coding genes (PCGs), 2 ribosomal RNA genes and 22 transfer RNA genes. The gene order and transcription direction are the same as G. spinigerum and G. doloresi. The sequence difference across the entire mt genomes varied from 14.4% to 18.2% between G. nipponicum, Gnathostoma sp., G. spinigerum and G. doloresi of Japan and China isolates. Phylogenetic analyses by Bayesian inference (BI) using concatenated amino acid sequences of 12 PCGs showed that G. nipponicum and Gnathostoma sp. are two distinctive species of Gnathostoma, and G. nipponicum are more closely related to Gnathostoma sp. than to G. spinigerum. The mtDNA datasets provide abundant resources of novel markers, which can be used for the studies of molecular epidemiology and diagnosis of Gnathostoma spp.

The complete mitochondrial genomes of Gnathostoma doloresi from China and Japan

Gnathostoma doloresi is one of the neglected pathogens causing gnathostomiasis. Although this zoonotic parasite leads to significant socioeconomic concerns globally, little is known of its genetics and systematics. In the present study, we sequenced and characterized the complete mitochondrial (mt) genomes of G. doloresi isolates from China and Japan. The lengths of the mt genomes of the G. doloresi China and Japan isolates are 13,809 and 13,812 bp, respectively. Both mt genomes encode 36 genes, including 12 protein-coding genes (PCGs), 2 ribosomal RNA genes, and 22 transfer RNA genes. The gene order, transcription direction, and genome content are identical with its congener G. spinigerum. Phylogenetic analyses based on concatenated amino acid sequences of 12 PCGs by Bayesian inference (BI) indicated that G. doloresi are closely related to G. spinigerum. Our data provide an invaluable resource for studying the molecular epidemiology, phylogenetics, and population genetics of Gnathostoma spp. and should have implications for further studies of the diagnosis, prevention, and control of gnathostomiasis in humans and animals.

Complete Mitochondrial genome of an equine intestinal parasite, Triodontophorus brevicauda (Chromadorea: Strongylidae): the first characterization within the genus

The complete mitochondrial (mt) genome sequence of Triodontophorus brevicauda, an intestinal equine nematode parasite was determined for the first time. The circular T. brevicauda mt genome is 14,305 bp in length and contains 36 genes, of which 12 code for protein, 22 for transfer RNA, and two for ribosomal RNA, and lacks atp8 mtDNA gene. Phylogenetic analysis based on the concatenated amino acid sequence of the 12 protein-coding genes was performed using three different tree-building methods. The Strongyloidea cluster divides into two large branches, and each nematode family included in our study forms an independent clade, though paraphyly confounds the issue at some nodes. T. brevicauda clusters together with Cylicocyclus insignis with high statistical support. The mtDNA data in this study not only provide a new mtDNA resource for phylogeny, but also become a novel and useful genetic marker for further studies on the identification, population genetics, and molecular epidemiology of the genus Triodontophorus in equine.

Mitochondrial genome plasticity among species of the nematode genus Meloidogyne (Nematoda: Tylenchina)

The mitochondrial (mt) genomes of the plant-parasitic root-knot nematodes Meloidogyne arenaria, Meloidogyne enterolobii and Meloidogyne javanica were sequenced and compared with those of three other root-knot nematode species in order to explore the mt genome plasticity within Meloidogyne. The mt genomes of M. arenaria, M. enterolobii and M. javanica are circular, with an estimated size of 18.8, 18.9 and 19.6 kb, respectively. Compared to other nematodes these mt genomes are larger, due to the presence of large non-coding regions. The mt genome architecture within the genus Meloidogyne varied in the position of trn genes and in the position, length and nucleotide composition of non-coding regions. These variations were observed independent of the species' natural environments or reproductive modes. M. enterolobii showed three main non-coding regions whereas Meloidogyne chitwoodi, Meloidogyne incognita, M. javanica and M. arenaria had two non-coding regions, and Meloidogyne graminicola had a unique large non-coding region interrupted by two trn genes. trn genes were positioned in different regions of the mt genomes in M. chitwoodi, M. enterolobii and M. graminicola, whereas the trn gene order was identical between M. arenaria, M. incognita and M. javanica. Importantly, M. graminicola had extra copies of trnV and trnS2. High divergence levels between the two copies of each trn might indicate duplication events followed by random loss and mutations in the anticodon. Tree-based methods based on amino acid sequences of 12 mt protein-coding genes support the monophyly for the tropical and mitotic parthenogenetic species, M. arenaria, M. enterolobii, M. incognita and M. javanica and for a clade that includes the meiotic parthenogenetic species, M. chitwoodi and M. graminicola. A comparison of the mt genome architecture in plant-parasitic nematodes and phylogenetic analyses support that Pratylenchus is the most recent ancestor of root-knot nematodes.

Characterization of the intergenic spacer rDNAs of two pig nodule worms, Oesophagostomum dentatum and O. quadrispinulatum

The characteristics of the intergenic spacer rDNAs (IGS rDNAs) of Oesophagostomum dentatum and O. quadrispinulatum isolated from pigs in different geographical locations in Mainland China were determined, and the phylogenetic relationships of the two species were reconstructed using the IGS rDNA sequences. The organization of the IGS rDNA sequences was similar to their organization in other eukaryotes. The 28S-18S IGS rDNA sequences of both O. dentatum and O. quadrispinulatum were found to have variable lengths, that is, 759-762 bp and 937-1128 bp, respectively. All of the sequences contained direct repeats and inverted repeats. The length polymorphisms were related to the different numbers and organization of repetitive elements. Different types and numbers of repeats were found between the two pig nodule species, and two IGS structures were found within O. quadrispinulatum. Phylogenetic analysis showed that all O. dentatum isolates were clustered into one clade, but O. quadrispinulatum isolates from different origins were grouped into two distinct clusters. These results suggested independent species and the existence of genotypes or subspecies within pig nodule worms. Different types and numbers of repeats and IGS rDNA structures could serve as potential markers for differentiating these two species of pig nodule worms.

Characterization of the complete mitochondrial genomes of Nematodirus oiratianus and Nematodirus spathiger of small ruminants

BACKGROUND

Nematodirus spp. are among the most common nematodes of ruminants worldwide. N. oiratianus and N. spathiger are distributed worldwide as highly prevalent gastrointestinal nematodes, which cause emerging health problems and economic losses. Accurate identification of Nematodirus species is essential to develop effective control strategies for Nematodirus infection in ruminants. Mitochondrial DNA (mtDNA) could provide powerful genetic markers for identifying these closely related species and resolving phylogenetic relationships at different taxonomic levels.

METHODS

In the present study, the complete mitochondrial (mt) genomes of N. oiratianus and N. spathiger from small ruminants in China were obtained using Long-range PCR and sequencing.

RESULTS

The complete mt genomes of N. oiratianus and N. spathiger were 13,765 bp and 13,519 bp in length, respectively. Both mt genomes were circular and consisted of 36 genes, including 12 genes encoding proteins, 2 genes encoding rRNA, and 22 genes encoding tRNA. Phylogenetic analyses based on the concatenated amino acid sequence data of all 12 protein-coding genes by Bayesian inference (BI), Maximum likelihood (ML) and Maximum parsimony (MP) showed that the two Nematodirus species (Molineidae) were closely related to Dictyocaulidae.

CONCLUSIONS

The availability of the complete mtDNA sequences of N. oiratianus and N. spathiger not only provides new mtDNA sources for a better understanding of nematode mt genomics and phylogeny, but also provides novel and useful genetic markers for studying diagnosis, population genetics and molecular epidemiology of Nematodirus spp. in small ruminants.

Comparative analyses of the complete mitochondrial genomes of the two ruminant hookworms Bunostomum trigonocephalum and Bunostomum phlebotomum

Bunostomum trigonocephalum and Bunostomum phlebotomum are blood-feeding hookworms of sheep and cattle, causing considerable economic losses to the live stock industries. Studying genetic variability within and among hookworm populations is critical to addressing epidemiological and ecological questions. Mitochondrial (mt) DNA is known to provide useful markers for investigations of population genetics of hookworms, but mt genome sequence data are scant. In the present study, the complete mitochondrial DNA (mtDNA) sequences of the sheep and goat hookworm B. trigonocephalum were determined for the first time, and the mt genome of B. phlebotomum from yak in China was also sequenced for comparative analyses of their gene contents and genome organizations. The lengths of mt DNA sequences of B. trigonocephalum sheep isolate, B. trigonocephalum goat isolate and B. phlebotomum China yak isolate were 13,764bp, 13,771bp and 13,803bp in size, respectively. The identity of the mt genomes was 99.7% between B. trigonocephalum sheep isolate and B. trigonocephalum goat isolate. The identity of B. phlebotomum China yak isolate mt genomes was 85.3% with B. trigonocephalum sheep isolate, and 85.2% with B. trigonocephalum goat isolate. All the mt genes of the two hookworms were transcribed in the same direction and gene arrangements were consistent with those of the GA3 type, including 12 protein-coding genes, 2 rRNA genes and 22 tRNA genes, but lacking ATP synthetase subunit 8 gene. The mt genomes of B. trigonocephalum and B. phlebotomum were similar to prefer bases A and T, the contents of A+T are 76.5% (sheep isolate), 76.4% (goat isolate) and 76.9% (China yak isolate), respectively. Phylogenetic relationships reconstructed using concatenated amino acid sequences of 12 protein-coding genes with three methods (maximum likelihood, Bayesian inference and neighbor joining) revealed that the B. trigonocephalum and B. phlebotomum represent distinct but closely-related species. These data provide novel and useful genetic markers for studying the systematics, and population genetics of the two ruminant hookworms.

Chabertia erschowi (Nematoda) is a distinct species based on nuclear ribosomal DNA sequences and mitochondrial DNA sequences

Background

Gastrointestinal nematodes of livestock have major socio-economic importance worldwide. In small ruminants, Chabertia spp. are responsible for economic losses to the livestock industries globally. Although much attention has given us insights into epidemiology, diagnosis, treatment and control of this parasite, over the years, only one species (C. ovina) has been accepted to infect small ruminants, and it is not clear whether C. erschowi is valid as a separate species.

Methods

The first and second internal transcribed spacers (ITS-1 and ITS-2) regions of nuclear ribosomal DNA (rDNA) and the complete mitochondrial (mt) genomes of C. ovina and C. erschowi were amplified and then sequenced. Phylogenetic re-construction of 15 Strongylida species (including C. erschowi) was carried out using Bayesian inference (BI) based on concatenated amino acid sequence datasets.

Results

The ITS rDNA sequences of C. ovina China isolates and C. erschowi samples were 852–854 bp and 862 -866 bp in length, respectively. The mt genome sequence of C. erschowi was 13,705 bp in length, which is 12 bp shorter than that of C. ovina China isolate. The sequence difference between the entire mt genome of C. ovina China isolate and that of C. erschowi was 15.33%. In addition, sequence comparison of the most conserved mt small subunit ribosomal (rrnS) and the least conserved nad2 genes among multiple individual nematodes revealed substantial nucleotide differences between these two species but limited sequence variation within each species.

Conclusions

The mtDNA and rDNA datasets provide robust genetic evidence that C. erschowi is a valid strongylid nematode species. The mtDNA and rDNA datasets presented in the present study provide useful novel markers for further studies of the taxonomy and systematics of the Chabertia species from different hosts and geographical regions.

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.

Characterization of the complete mitochondrial genome of Spirocerca lupi: sequence, gene organization and phylogenetic implications

Background

Spirocerca lupi is a life-threating parasitic nematode of dogs that has a cosmopolitan distribution but is most prevalent in tropical and subtropical countries. Despite its veterinary importance in canids, the epidemiology, molecular ecology and population genetics of this parasite still remain unexplored.

Methods

The complete mitochondrial (mt) genome of S. lupi was amplified in four overlapping long fragments using primers designed based on partial cox1, rrnS, cox2 and nad2 sequences. Phylogenetic re-construction of 13 spirurid species (including S. lupi) was carried out using Bayesian inference (BI) based on concatenated amino acid sequence datasets.

Results

The complete mt genome sequence of S. lupi is 13,780 bp in length, including 12 protein-coding genes, 22 transfer RNA genes and two ribosomal RNA genes, but lacks the atp8 gene. The gene arrangement is identical to that of Thelazia callipaeda (Thelaziidae) and Setaria digitata (Onchocercidae), but distinct from that of Dracunculus medinensis (Dracunculidae) and Heliconema longissimum (Physalopteridae). All genes are transcribed in the same direction and have a nucleotide composition high in A and T. The content of A + T is 73.73% for S. lupi, in accordance with mt genomes of other spirurid nematodes sequenced to date. Phylogenetic analyses using concatenated amino acid sequences of the 12 protein-coding genes by BI showed that the S. lupi (Thelaziidae) is closely related to the families Setariidae and Onchocercidae.

Conclusions

The present study determined the complete mt genome sequence of S. lupi. These new mt genome dataset should provide novel mtDNA markers for studying the molecular epidemiology and population genetics of this parasite, and should have implications for the molecular diagnosis, prevention and control of spirocercosis in dogs and other canids.

Identification and characterization of new major sperm protein genes from Oesophagostomum dentatum and Oesophagostomum quadrispinulatum from pigs in China

The present study identified and characterized new major sperm protein (MSP) genes from the two nodule worms Oesophagostomum dentatum and Oesophagostomum quadrispinulatum collected from pigs in China. Total genomic DNA was extracted individually from 10 male nematode samples representing O. dentatum, and 4 male nematode samples representing O. quadrispinulatum. A pair of primers (OMSP1F/MSP1R) was designed based on the MSP gene sequences of Ascaris suum and O. dentatum available in GenBank, and used to amplify the MSP genes from the two porcine nodule worms. The PCR products were purified and subsequently cloned into pGEM-T Easy vector. Recombinants were identified by PCR and sequenced. Sequence analysis revealed that there were two different types of MSP sequences in O. dentatum and O. quadrispinulatum, one contained intron, and the other did not. The lengths of the MSP sequences containing introns were 433 bp or 439 bp in O. dentatum, and 436 bp, 439 bp or 446 bp in O. quadrispinulatum, containing 1 or 2 introns. Five and three new members of the MSP multigene family were identified in O. dentatum and O. quadrispinulatum in this study, respectively. The MSP sequences without introns were 381 bp in length, and can be deduced into 126 amino acids. The sequences of MSP genes containing introns seem to be more conserved than those without introns. The identities of deduced amino acid sequences of the MSP genes containing introns were 96.0-100% within and between the two nodule worms, and were 81.1-93.7% compared with other published MSP sequences of the representative nematodes. The present study identified new MSP genes with introns from O. dentatum and O. quadrispinulatum for the first time. The identification and characterization of newly described MSP genes from O. dentatum and O. quadrispinulatum have implications for further studies of molecular biology and reproduction control of Oesophagostomum spp.