Complete mitochondrial genomes of Baylisascaris schroederi, Baylisascaris ailuri and Baylisascaris transfuga from giant panda, red panda and polar bear

Characterisation of the mitochondrial genome and phylogenetic analysis of Toxocara apodemi (Nematoda: Ascarididae)

We first sequenced and characterised the complete mitochondrial genome of Toxocara apodeme, then studied the evolutionary relationship of the species within Toxocaridae. The complete mitochondrial genome was amplified using PCR with 14 specific primers. The mitogenome length was 14303 bp in size, including 12 PCGs (encoding 3,423 amino acids), 22 tRNAs, 2 rRNAs, and 2 NCRs, with 68.38% A+T contents. The mt genomes of T. apodemi had relatively compact structures with 11 intergenic spacers and 5 overlaps. Comparative analyses of the nucleotide sequences of complete mt genomes showed that T. apodemi had higher identities with T. canis than other congeners. A sliding window analysis of 12 PCGs among 5 Toxocara species indicated that nad4 had the highest sequence divergence, and cox1 was the least variable gene. Relative synonymous codon usage showed that UUG, ACU, CCU, CGU, and UCU most frequently occurred in the complete genomes of T. apodemi. The Ka/Ks ratio showed that all Toxocara mt genes were subject to purification selection. The largest genetic distance between T. apodemi and the other 4 congeneric species was found in nad2, and the smallest was found in cox2. Phylogenetic analyses based on the concatenated amino acid sequences of 12 PCGs demonstrated that T. apodemi formed a distinct branch and was always a sister taxon to other congeneric species. The present study determined the complete mt genome sequences of T. apodemi, which provide novel genetic markers for further studies of the taxonomy, population genetics, and systematics of the Toxocaridae nematodes.

Morphology and ASAP analysis of the important zoonotic nematode parasite Baylisascaris procyonis (Stefahski and Zarnowski, 1951), with molecular phylogenetic relationships of Baylisascaris species (Nematoda: Ascaridida)

Species of Baylisascaris (Nematoda: Ascarididae) are of great veterinary and zoonotic significance, owing to cause Baylisascariosis or Baylisascariasis in wildlife, captive animals and humans. However, the phylogenetic relationships of the current 10 Baylisascaris species remain unclear. Moreover, our current knowledge of the detailed morphology and morphometrics of the important zoonotic species B. procyonis is still insufficient. The taxonomical status of B. procyonis and B. columnaris remains under debate. In the present study, the detailed morphology of B. procyonis was studied using light and scanning electron microscopy based on newly collected specimens from the raccoon Procyon lotor (Linnaeus) in China. The results of the ASAP analysis and Bayesian inference (BI) using the 28S, ITS, cox1 and cox2 genetic markers did not support that B. procyonis and B. columnaris represent two distinct species. Integrative morphological and molecular assessment challenged the validity of B. procyonis, and suggested that B. procyonis seems to represent a synonym of B. columnaris. Molecular phylogenetic results indicated that the species of Baylisascaris were grouped into 4 clades according to their host specificity. The present study provided new insights into the taxonomic status of B. procyonis and preliminarily clarified the phylogenetic relationships of Baylisascaris species.

Mitogenomic phylogenies suggest the resurrection of the subfamily Porrocaecinae and provide insights into the systematics of the superfamily Ascaridoidea (Nematoda: Ascaridomorpha), with the description of a new species of Porrocaecum

BACKGROUND

The family Toxocaridae is a group of zooparasitic nematodes of veterinary, medical and economic significance. However, the evolutionary relationship of Porrocaecum and Toxocara, both genera currently classified in Toxocaridae, and the monophyly of the Toxocaridae remain under debate. Moreover, the validity of the subgenus Laymanicaecum in the genus Porrocaecum is open to question. Due to the scarcity of an available genetic database, molecular identification of Porrocaecum nematodes is still in its infancy.

METHODS

A number of Porrocaecum nematodes collected from the Eurasian marsh harrier Circus aeruginosus (Linnaeus) (Falconiformes: Accipitridae) in the Czech Republic were identified using integrated morphological methods (light and scanning electron microscopy) and molecular techniques (sequencing and analyzing the nuclear 18S, 28S and ITS regions). The complete mitochondrial genomes of the collected nematode specimens and of Porrocaecum (Laymanicaecum) reticulatum (Linstow, 1899) were sequenced and annotated for the first time. Phylogenetic analyses of ascaridoid nematodes based on the amino acid sequences of 12 protein-coding genes of mitochondrial genomes were performed using maximum likelihood and Bayesian inference.

RESULTS

A new species of Porrocaecum, named P. moraveci n. sp., is described based on the morphological and genetic evidence. The mitogenomes of P. moraveci n. sp. and P. reticulatum both contain 36 genes and are 14,517 and 14,210 bp in length, respectively. Comparative mitogenomics revealed that P. moraveci n. sp. represents the first known species with three non-coding regions and that P. reticulatum has the lowest overall A + T content in the mitogenomes of ascaridoid nematodes tested to date. Phylogenetic analyses showed the representatives of Toxocara clustered together with species of the family Ascarididae rather than with Porrocaecum and that P. moraveci n. sp. is a sister to P. reticulatum.

CONCLUSIONS

The characterization of the complete mitochondrial genomes of P. moraveci n. sp. and P. reticulatum is reported for the first time. Mitogenomic phylogeny analyses indicated that the family Toxocaridae is non-monophyletic and that the genera Porrocaecum and Toxocara do not have an affinity. The validity of the subgenus Laymanicaecum in Porrocaecum was also rejected. Our results suggest that: (i) Toxocaridae should be degraded to a subfamily of the Ascarididae that includes only the genus Toxocara; and (ii) the subfamily Porrocaecinae should be resurrected to include only the genus Porrocaecum. The present study enriches the database of ascaridoid mitogenomes and provides a new insight into the systematics of the superfamily Ascaridoidea.

The Complete Mitogenome of Toxocara vitulorum: Novel In-Sights into the Phylogenetics in Toxocaridae

Toxocara vitulorum (Ascaridida: Nematoda) is one of the most common intestinal nematodes of cattle and buffalos and, therefore, represents a serious threat to their populations worldwide. Despite its significance in veterinary health the epidemiology, population genetics, and molecular ecology of this nematode remain poorly understood. The mitogenome can yield a foundation for studying these areas and assist in the surveillance and control of T. vitulorum. Herein, the first whole mitogenome of T. vitulorum was sequenced utilizing Illumina technology and characterized with bioinformatic pipeline analyses. The entire genome of T. vitulorum was 15,045 bp in length and contained 12 protein-coding genes (PCGs), 22 transfer RNAs (tRNAs), and two ribosomal RNAs (rRNAs). The gene arrangement (GA) of T. vitulorum was similar to those of other Toxocara species under GA3. The whole genome showed significant levels of AT and GC skew. Comparative mitogenomics including sequence identities, Ka/Ks, and sliding window analysis, indicated a purifying selection of 12 PCGs with cox1 and nad6 having the lowest and highest evolutionary rate, respectively. Whole amino acid sequence-based phylogenetic analysis supported a novel sister-species relationship of T. vitulorum with the congeneric species Toxocara canis, Toxocara cati, and Toxocara malaysiensis in the family Toxocaridae. Further, 12 (PCGs) single gene-based phylogenies suggested that nad4 and nad6 genes shared same topological trees with that of the whole genome, suggesting that these genes were suitable as novel genetic markers for phylogenetic and evolutionary studies of Ascaridida species. This complete mitogenome of T. vitulorum refined phylogenetic relationships in Toxocaridae and provided the resource of markers for population genetics, systematics, and epidemiology of this bovine nematode.

Gene rearrangements in the mitochondrial genome of ten ascaris species and phylogenetic implications for Ascaridoidea and Heterakoidea families

The Ascaridoidea family and Heterakoidea family are the most common and typical representative of large parasites. Although our understanding of these parasites' diversity has expanded by analyses of some mitochondrial genes, there is limited information on these species' evolutionary rates. Here we determined ten complete mitogenome sequences of five subfamilies of Ascaridoidea and one subfamily of Heterakoidea. The phylogenetic tree divided the Ascaridoidea into six monophyletic major clades, and the divergence time of Heterakoidea family and Ascaridoidea family can be placed during the early Carboniferous Period (300-360 Mya). The reconstruction of the ancestral state showed that the gene orders of all species in Ascaridoidea were conserved, and the Heterakoidea had obvious genome rearrangement. The conserved blocks between them were divided into five and the main types are tandem-duplication/random loss (TDRL). These results will help to better understand the gene rearrangements and evolutionary position of ascaris species.

Population genetic structure and morphological diversity of Cruzia tentaculata (Nematoda: Ascaridida), a parasite of marsupials (Didelphinae), along the Atlantic Forest on the eastern coast of South America

Cruzia tentaculata is a helminth parasite of marsupials and has a wide geographic distribution from Mexico to Argentina. The aim of this study was to analyse the genetic population structure of this nematode along the Atlantic Forest biome. Cruzia tentaculata specimens were recovered from Didelphis aurita, Didelphis albiventris and Philander quica in 9 localities. Morphological and morphometric data were investigated for phenotypic diversity among localities and hosts using multivariate discriminant analysis of principal components. Phylogenetic relationships of C. tentaculata were determined using maximum likelihood and Bayesian inference. The population structure was analysed by fixation indices, molecular variance analysis, Tajima's D and Fu's Fs neutrality tests, Mantel tests and Bayesian clustering analysis. A higher significant morphometric difference for males was observed between localities. In the haplogroup networks, 2 groups were recovered, separating locations from the north and from the south/southeast. The morphometric variation in C. tentaculata between different localities was compatible with this north and southeast/south pattern, suggesting adaptation to different ecological conditions. Population genetic analyses suggested a pattern of evolutionary processes driven by Pleistocene glacial refugia in the northeast and southeast of the Atlantic Forest based on the distribution of genetic diversity.

Sequence analyses of mitochondrial gene may support the existence of cryptic species within Ascaridia galli

Ascaridia galli (Nematoda: Ascaridiidae) is the most common intestinal roundworm of chickens and other birds with a worldwide distribution. Although A. galli has been extensively studied, knowledge of the genetic variation of this parasite in detail is still insufficient. The present study examined genetic variation in the mitochondrial cytochrome c oxidase subunit 1 (cox1) gene among A. galli isolates (n = 26) from domestic chickens in Hunan Province, China. A portion of the cox1 (pcox1) gene was amplified by polymerase chain reaction separately from adult A. galli individuals and the amplicons were subjected to sequencing from both directions. The length of the sequences of pcox1 is 441 bp. Although the intra-specific sequence variation within A. galli is 0-7.7%, the inter-specific sequence differences among other members of the infraorder Ascaridomorpha were 11.4-18.9%. Phylogenetic analyses based on the maximum likelihood method using the sequences of pcox1 confirmed that all of the Ascaridia isolates were A. galli, and also resolved three distinct clades. Taken together, the findings suggest that A. galli may represent a complex of cryptic species. Our results provide an additional genetic marker for the management of A. galli in chickens and other birds.

Regulatory effects of a novel cysteine protease inhibitor in Baylisascaris schroederi migratory larvae on mice immune cells

Background

The giant panda (Ailuropoda melanoleuca) is a well-known, rare and endangered species. Baylisascaris schroederi is a pathogenic ascarid. Infection with B. schroederi may cause death in giant pandas. At present, the immune evasion mechanism of B. schroederi is little known. Cysteine protease inhibitors (CPI) play important roles in the regulation of host immune responses against certain nematodes. In this study, we focused on the analysis of the regulation of B. schroederi migratory larvae CPI (rBsCPI-1) on mice immune cells.

Methods

First, the pattern recognition receptors on the surface of peripheral blood mononuclear cells (PBMCs) and the signal pathways that transduce extracellular signals into the nucleus activated by rBsCPI-1 were identified. Then, the regulatory effects of rBsCPI-1 on PBMCs physiological activities were detected. Finally, the effects of rBsCPI-1 on TLR signaling pathway activation and NF-κB phosphorylation in mice immunized with recombinant protein were analysed.

Results

The results suggested that rBsCPI-1 secreted by B. schroederi migratory larvae is mainly recognized by TLR2 and TLR4 on PBMCs. Extracellular signals are transduced into the nucleus through the MAPK and NF-κB signaling pathways, enhancing the phagocytosis, migration, and apoptosis of PBMCs; meanwhile, rBsCPI-1 induces high expression of NO. Thus, rBsCPI-1 plays a role in immune regulation. In addition, the high expression of negative regulatory factors also ensured that TLR activation is maintained at the optimal level.

Conclusions

rBsCPI-1 can transduce regulatory signals into immune cells by activating the TLR2/4-NF-κB/MAPK signaling pathway, having a certain regulatory effect on the physiological activities. Meanwhile, rBsCPI-1 can maintain the immune response in a balance by limiting the over-activation of the TLRs signaling pathway and thus contributes to B. schroederi immune evasion.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-022-05240-8.

Genome of the Giant Panda Roundworm Illuminates Its Host Shift and Parasitic Adaptation

Baylisascaris schroederi, a roundworm (ascaridoid) parasite specific to the bamboo-feeding giant panda (Ailuropoda melanoleuca), represents a leading cause of mortality in wild giant panda populations. Here, we present a 293-megabase chromosome-level genome assembly of B. schroederi to infer its biology, including host adaptations. Comparative genomics revealed an evolutionary trajectory accompanied by host-shift events in ascaridoid parasite lineages after host separations, suggesting their potential for transmission and rapid adaptation to new hosts. Genomic and anatomical lines of evidence, including expansion and positive selection of genes related to the cuticle and basal metabolisms, indicate that B. schroederi undergoes specific adaptations to survive in the sharp-edged bamboo-enriched gut of giant pandas by structurally increasing its cuticle thickness and efficiently utilizing host nutrients through gut parasitism. Additionally, we characterized the secretome of B. schroederi and predicted potential drug and vaccine targets for new control strategies. Overall, this genome resource provides new insights into the host adaptation of B. schroederi to the giant panda as well as the host-shift events in ascaridoid parasite lineages. Our findings on the unique biology of B. schroederi will also aid in the development of prevention and treatment measures to protect giant panda populations from roundworm parasitism.

Chromosome-scale assembly and whole-genome sequencing of 266 giant panda roundworms provide insights into their evolution, adaptation and potential drug targets

Helminth diseases have long been a threat to the health of humans and animals. Roundworms are important organisms for studying parasitic mechanisms, disease transmission and prevention. The study of parasites in the giant panda is of importance for understanding how roundworms adapt to the host. Here, we report a high‐quality chromosome‐scale genome of Baylisascaris schroederi with a genome size of 253.60 Mb and 19,262 predicted protein‐coding genes. We found that gene families related to epidermal chitin synthesis and environmental information processes in the roundworm genome have expanded significantly. Furthermore, we demonstrated unique genes involved in essential amino acid metabolism in the B. schroederi genome, inferred to be essential for the adaptation to the giant panda‐specific diet. In addition, under different deworming pressures, we found that four resistance‐related genes (glc‐1, nrf‐6, bre‐4 and ced‐7) were under strong positive selection in a captive population. Finally, 23 known drug targets and 47 potential drug target proteins were identified. The genome provides a unique reference for inferring the early evolution of roundworms and their adaptation to the host. Population genetic analysis and drug sensitivity prediction provide insights revealing the impact of deworming history on population genetic structure of importance for disease prevention.

PHYLOGEOGRAPHY OF BAYLISASCARIS PROCYONIS (RACCOON ROUNDWORM) IN NORTH AMERICA

Sequences of the mitochondrial cytochrome c oxidase 1 (COI) gene of 115 Baylisascaris procyonis individuals from 13 U.S. states and 1 Canadian province were obtained from 44 raccoon hosts to assess genetic variation and geographic structure. The maximum genetic distance between individuals was low (1.6%), consistent with a single species. Moderate COI haplotype (h = 0.60) and nucleotide (π = 0.0053) diversity were found overall. Low haplotype diversity was found among samples east of the Mississippi River (h = 0.036), suggesting that historical growth and expansion of raccoon populations in this region could be responsible for high parasite gene flow or a selective sweep of B. procyonis mtDNA. There was low genetic structure (average Φst = 0.07) for samples east of the continental divide, but samples from Colorado showed higher diversity and differentiation from midwestern and eastern samples. There was marked genetic structure between samples from east and west of the continental divide, with no haplotypes shared between these regions. There was no significant isolation by distance among any of these geographic samples. The phylogeographic patterns for B. procyonis are similar to genetic results reported for their raccoon definitive hosts. The phylogeographic divergence of B. procyonis from east and west of the continental divide may involve vicariance resulting from Pleistocene glaciation and associated climate variation.

The mitogenome of Ophidascaris wangi isolated from snakes in China

Different species of the genus Ophidascaris (Baylis, 1921; Nematoda: Ascaridida, Ascaridoidea) are intestinal parasites of various snake species. More than 30 Ophidascaris species have been reported worldwide; however, few molecular genetic studies have been conducted on this genus. We sequenced the complete mitogenome of Ophidascaris wangi parasitizing two snake species of the family Colubridae, i.e., Elaphe carinata (Günther, 1864) and Dinodon rufozonatum. The mitogenome sequence of O. wangi was approximately 14,660 base pairs (bp) long and encoded 36 genes, including 12 protein-coding genes (PCGs), 2 ribosomal RNA (rRNA) genes, and 22 transfer RNA genes. Gene arrangement, genome content, and transcription direction were in line with those in Toxascaris leonina (Linstow, 1902; Ascaridida: Ascarididae). Phylogenetics of O. wangi and other ascaridoids were reconstructed based on the concatenated amino acid sequences of 12 PCGs, and on nucleotide sequences of 12 PCGs and two rRNA genes. Phylogenetic analyses were performed using maximum likelihood and Bayesian inference methods, and the results suggested that O. wangi constitutes a sister clade of Ascaris, Parascaris, Baylisascaris, and Toxascaris within the family Ascarididae, which is a sister clade of Toxocaridae. The mitogenome sequence of O. wangi obtained from the present study will be useful for future identification of the nematode worms in the genus Ophidascaris and will increase the understanding of population genetics, molecular epidemiology, and phylogenetics of ascaridoid nematodes in snakes.

The mitochondrial genome of the sheep roundworm Ascaris ovis (Ascaridida: Nematoda) from Southwest China

The Ascaris roundworms (Ascaridida: Nematoda), one of the commonest soil-transmitted helminths (STHs), can cause ascariasis with significant socioeconomic and public health impact. In this study, the mitochondrial genome of Ascaris ovis, a representative of this genus from the sheep in Southwest China was determined using Illumina sequencing technology. The assembled genome was 14,205 bp in size and encoded 36 genes, including 12 protein-coding genes, 22 transfer RNA genes, and two ribosomal RNA genes. Phylogenetic analysis showed that A. ovis grouped with the congeneric Ascaris lumbricoides of humans, Ascaris spp. of non-human primates and Ascaris suum of pigs and together formed a monophyletic group relationship with either species of Baylisascaris/Toxascaris/Parascaris, species of Toxocara, species of Anisakis/Pseudoterranova or species of Ascaridia/Heterakis in the order Ascaridida, supporting its genetic similarity with A. lumbricoides, A. suum, and other congeneric species. The cumulative mitochondrial DNA data of this genus should contribute to a better understanding of the phylogenetic relationships among these roundworms.

The morphological and molecular characterization of Baylisascaris devosi Sprent, 1952 (Ascaridoidea, Nematoda), collected from Pine marten (Martes martes) in Iran

Background

Baylisascaris devosi is an intestinal nematode found in several carnivores including fisher, wolverine, Beech marten, American marten and sable in different parts of the world, but this nematode has not been reported from Pine marten. Therefore, this study aimed to identify Baylisascaris isolated from a Pine marten in Iran using morphological and molecular approaches.

Methods

Specimens of B. devosi were collected from one road-killed Pine marten in northern Iran. Morphological features were evaluated using scanning electron microscopy, energy dispersive x-ray analysis and ion sectioning. The molecular characterization was carried out using partial Cox1, LSU rDNA and ITS-rDNA genes.

Results

The nematodes isolated from the Pine marten were confirmed to be B. devosi based on the morphological features and the sequence of ribosomal and mitochondrial loci. X-ray scans (EDAX) were completed on gallium cut structures (papillae, eggs, male spike and mouth denticles) of B. devosi using a dual-beam scanning electron microscope. The male spike and mouth denticles had a high level of hardening elements (Ca, P, S), helping to explain the chemical nature and morphology of the worm. Based on these genetic marker analyses, our sequence had the greatest similarity with Russian B. devosi isolated from sable.

Conclusions

In this study, to our knowledge, the occurrence of B. devosi infection in Pine marten is reported for the first time. Molecular analysis showed that these three genes are suitable molecular markers for identification and inferring phylogenetic relationships of Baylisascaris species. Furthermore, the high divergence of Cox1 between Baylisascaris species indicates that Cox1 could be used for their phylogenetic and taxonomic studies.

Molecular characterization of ascaridoid parasites from captive wild carnivores in China using ribosomal and mitochondrial sequences

BACKGROUND

Despite the public health importance of toxocariasis/toxascariasis, only a few species of these ascaridoid parasites from wild canine and feline carnivores have been studied at the molecular level so far. Poor understanding of diversity, host distribution and the potential (zoonotic) transmission of the ascaridoid species among wild animals negatively affects their surveillance and control in natural settings. In this study, we updated previous knowledge by profiling the genetic diversity and phylogenetic relationships of ascaridoid species among eleven wild canine and feline animals on the basis of a combined analysis of the ribosomal internal transcribed spacer region (ITS) gene and the partial mitochondrial cytochrome c oxidase subunit 2 (cox2) and NADH dehydrogenase subunit 1 (nad1) genes.

RESULTS

In total, three genetically distinct ascaridoid lineages were determined to be present among these wild carnivores sampled, including Toxocara canis in Alopex lagopus and Vulpes vulpes, Toxocara cati in Felis chaus, Prionailurus bengalensis and Catopuma temmincki and Toxascaris leonina in Canis lupus, Panthera tigris altaica, Panthera tigris amoyensis, Panthera tigris tigris, Panthera leo and Lynx lynx. Furthermore, it was evident that T. leonina lineage split into three well-supported subclades depending on their host species, i.e. wild felids, dogs and wolves and foxes, based on integrated genetic and phylogenetic evidence, supporting that a complex of T. leonina other than one species infecting these hosts.

CONCLUSIONS

These results provide new molecular insights into classification, phylogenetic relationships and epidemiological importance of ascaridoids from wild canids and felids and also highlight the complex of the taxonomy and genetics of Toxascaris in their wild and domestic carnivorous hosts.

Review on parasites of wild and captive giant pandas (Ailuropoda melanoleuca): Diversity, disease and conservation impact

The giant panda (Ailuropoda melanoleuca) is a rare species with a small global population size, and lives in the wild in only a few fragmented mountain ranges of Southwest China. Parasitic infections are among the important causes of death of giant pandas that hamper their group development. We reviewed the parasitic infections prevailing in giant pandas, and the parasitic diversity, diseases and their impact on conservation of this animal. A total of 35 parasitic species were documented in giant pandas, belonging to nematode (n = 6), trematode (n = 1), cestode (n = 2), protozoa (n = 9), and ectozoa (n = 17 (tick = 13, mite = 2, and flea = 2)). Among them, Baylisascaris schroederi had the highest prevalence and was the leading cause of death for giant pandas. Some parasites caused asymptomatic infections in giant pandas, and their health implications for the pandas remain unknown. As a whole, parasites are reported to be an important threat to the conservation of the giant pandas. Regular deworming and environmental disinfection appear to be effective ways to prevent captive giant pandas from parasitoses. In wild panda populations, parasitic control measures are suggested to include detailed examination of the ecology of the host-parasite assembly, with particular attention to density-dependent transmission. The parasitic pathogenesis and detection methods together with their biology, epidemiology, treatment, prevention and control need to be further studied for better protection of giant pandas from parasitoses.

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.

Molecular phylogenetics and species-level systematics of Baylisascaris

Nucleotide sequences representing nine genes and five presumptive genetic loci were used to infer phylogenetic relationships among seven Baylisascaris species, including one species with no previously available molecular data. These genes were used to test the species status of B. procyonis and B. columnaris using a coalescent approach. Phylogenetic analysis based on combined analysis of sequence data strongly supported monophyly of the genus and separated the species into two main clades. Clade 1 included B. procyonis, B. columnaris, and B. devosi, species hosted by musteloid carnivores. Clade 2 included B. transfuga and B. schroederi from ursids, B. ailuri, a species from the red panda (a musteloid), and B. tasmaniensis from a marsupial. Within clade 2, geographic isolates of B. transfuga, B. schroederi (from giant panda), and B. ailuri formed a strongly supported clade. In certain analyses (e.g., some single genes), B. tasmaniensis was sister to all other Baylisascaris species rather than sister to the species from ursids and red panda. Using one combination of priors corresponding to moderate population size and shallow genetic divergence, the multispecies coalescent analysis of B. procyonis and B. columnaris yielded moderate support (posterior probability 0.91) for these taxa as separate species. However, other prior combinations yielded weak or no support for delimiting these taxa as separate species. Similarly, tree topologies constrained to represent reciprocal monophyly of B. columnaris and B. procyonis individuals (topologies consistent with separate species) were significantly worse in some cases, but not others, depending on the dataset analyzed. An expanded analysis of SNPs and other genetic markers that were previously suggested to distinguish between individuals of B. procyonis and B. columnaris was made by characterization of additional individual nematodes. The results suggest that many of these SNPs do not represent fixed differences between nematodes derived from raccoon and skunk hosts.

•

A phylogenetic hypothesis for Baylisascaris species was produced using nine genes.

•

Genetic data was generated for two new species- B. devosi and B. tasmaniensis.

•

Baylisascaris devosi and B. tasmaniensis were part of a monophyletic Baylisascaris.

•

B. procyonis (raccoon) and B. columnaris (skunk) could not be reliably distinguished.

•

Established SNPs may not be diagnostic for Baylisascaris from raccoons and skunks.

Genetic characterisation and phylogenetic status of whipworms (Trichuris spp.) from captive non-human primates in China, determined by nuclear and mitochondrial sequencing

Background

Whipworms (Nematoda: Trichuridae), among the most common soil-transmitted helminths (STHs), can cause the socioeconomically important disease trichuriasis in various mammalian hosts including humans and non-human primates. For many years, Trichuris from non-human primates has been assigned to the same species as the one infecting humans Trichuris trichiura. More recently, several molecular reports challenged this assumption following recognition of a Trichuris species complex observed in humans and non-human primates. A refined concept for species limits within Trichuris contributes to an understanding of diversity and the potential (zoonotic) transmission among humans and non-human primates. In this study, we expanded previous investigations by exploring the diversity of Trichuris among eight primates including three Asian autochthonous species (i.e. Rhinopithecus roxellana, Rhinopithecus bieti and Nomascus leucogenys). Species-level identification, whether novel or assignable to known lineages of Trichuris, was based on analyses of nuclear internal transcribed spacers (ITS) and mitochondrial cytochrome c oxidase subunit 1 (cox1) genes.

Results

In total, seven genetically distinct subgroups of whipworms were determined to be present among the primates sampled. Most Trichuris lineages, including Subgroups 1, 1’, 3, 5 and 6, showed a broad host range and were not restricted to particular primate species; in addition to T. trichiura, a complex of Trichuris species was shown infecting primates. Furthermore, it was assumed that Trichuris spp. from either N. leucogenys and P. hamadryas or R. roxellana and R. bieti, respectively, were conspecific. Each pair was indicated to be a discrete lineage of Trichuris, designated, respectively, as Subgroups 1 or 1’ and 2, based on integrated genetic and phylogenetic evidence.

Conclusion

These results emphasise that the taxonomy and genetic variations of Trichuris are more complicated than previously acknowledged. These cumulative molecular and phylogenetic data provide a better understanding of the taxonomy, genetics and evolutionary biology of the whipworms.

Electronic supplementary material

The online version of this article (10.1186/s13071-018-3100-5) contains supplementary material, which is available to authorized users.

A frequent roundworm Baylisascaris transfuga in overpopulated brown bears (Ursus arctos) in Slovakia: a problem worthy of attention

The genus Baylisascaris (order Ascaridida) includes numerous relatively host-specific nematodes, which are common in intestines of wild mammals. Some of them may have impact on veterinary and public health, as their larvae have the potential to cause visceral, ocular, and/or neural larva migrans in a wide range of mammals, birds, and humans. Baylisascaris transfuga is a parasite occurring in a range of bear species throughout the world. We present the current data on B. transfuga occurrence in brown bears from a relatively restricted territory of the Poľana Protected Landscape Area in Central Slovakia, obtained by traditional methods (faecal examination, morphology). Species affiliation was confirmed by employing molecular markers generating nuclear 28S and mitochondrial cox1 sequences in adult worms. Based on 17 examined samples (15 excrements and two intestines of young bear females), the occurrence of B. transfuga in the surveyed area was assessed as 52.9%. Both bear females were infected with adult and juvenile worms. Due to the high density of bears in the locality, the high infection rate with ascarids, and the huge number of eggs produced by the parasites, it is apparent that the respective environment, including the inhabited areas, might be markedly contaminated by Baylisascaris eggs. The ability of B. transfuga to serve as a zoonotic agent has not been unambiguously proved; however, this attribute should be considered and subjected to further research.

Parasites of the Giant Panda: A Risk Factor in the Conservation of a Species

The giant panda, with an estimated population size of 2239 in the world (in 2015), is a global symbol of wildlife conservation that is threatened by habitat loss, poor reproduction and limited resistance to some infectious diseases. Of these factors, some diseases caused by parasites are considered as the foremost threat to its conservation. However, there is surprisingly little published information on the parasites of the giant panda, most of which has been disseminated in the Chinese literature. Herein, we review all peer-reviewed publications (in English or Chinese language) and governmental documents for information on parasites of the giant pandas, with an emphasis on the intestinal nematode Baylisascaris schroederi (McIntosh, 1939) as it dominates published literature. The purpose of this chapter is to: (i) review the parasites recorded in the giant panda and describe what is known about their biology; (ii) discuss key aspects of the pathogenesis, diagnosis, treatment and control of key parasites that are reported to cause clinical problems and (iii) conclude by making some suggestions for future research. This chapter shows that we are only just 'scratching the surface' when it comes to parasites and parasitological research of the giant panda. Clearly, there needs to be a concerted research effort to support the conservation of this iconic species.

Status of gastrointestinal parasites in Red Panda of Nepal

Red pandas are known to be highly susceptible to endoparasites, which can have a prominent impact on the population dynamics of this endangered species. There are very limited published reports on prevalence and risk of parasites in wild populations of red panda, especially localized reports. This study attempts to provide an in-depth insight of the status of endoparasites in red pandas, which is critical for strengthening conservation efforts. A total of 272 fecal samples were collected through systematic sampling across the red panda distribution range in Nepal and coprological examination was completed using standard techniques. It was followed by an estimation of prevalence and mean intensity of parasites, as well as statistical analysis, which was carried out using R statistical software. Parasite prevalence was documented in 90.80% (n = 247) out of 272 samples examined which includes seven different species along with three genera of parasites belonging to Protozoans (3 species), Cestodes (1 genus, 1 species) and Nematodes (2 genera, 3 species). Nematodes predominated in all infected samples (87.62%). Prevalence of Ancyclostoma duodenale (n = 227, 70.06%), having a mean intensity of 3.45 ± 2.88 individuals per sample, was observed, followed by Ascaris lumbricoides (n = 19, 5.86%) and Entamoeba histolytica (n = 24, 7.41%). Eight variables for assessing the determinants of infestation were tested: protected areas; non-protected areas; aspect; elevation; slope; and distance to water sources, herding stations, and settlements. Only the settlement displayed significant association (β = −1534e−04, t =  − 2.192, p = 0.0293) though each parasite species displayed dissimilar association with different variables. This study indicates the urgent need of improving existing herding practice through habitat zonation, rotational grazing, medication of livestock, and prohibition of open defecation within and around red panda habitat.

Beyond the raccoon roundworm: The natural history of non-raccoon Baylisascaris species in the New World

A total of 10 species of Baylisascaris, a genus of ascaridoid nematodes, occur worldwide and 6 of them occur in the New World. Most of the Baylisascaris species have a similar life cycle with carnivorous mammals or marsupials serving as definitive hosts and a smaller prey host serving as paratenic (or intermediate) hosts. However, one species in rodents is unique in that it only has one host. Considerable research has been conducted on B. procyonis, the raccoon roundworm, as it is a well-known cause of severe to fatal neurologic disease in humans and many wildlife species. However, other Baylisascaris species could cause larva migrans but research on them is limited in comparison. In addition to concerns related to the potential impacts of larva migrans on potential paratenic hosts, there are many questions about the geographic ranges, definitive and paratenic host diversity, and general ecology of these non-raccoon Baylisascaris species. Here, we provide a comprehensive review of the current knowledge of New World Baylisascaris species, including B. columnaris of skunks, B. transfuga and B. venezuelensis of bears, B. laevis of sciurids, B. devosi of gulonids, B. melis of badgers, and B. potosis of kinkajou. Discussed are what is known regarding the morphology, host range, geographic distribution, ecoepidemiology, infection dynamics in definitive and paratenic hosts, treatment, and control of these under-studied species. Also, we discuss the currently used molecular tools used to investigate this group of parasites. Because of morphologic similarities among larval stages of sympatric Baylisascaris species, these molecular tools should provide critical insight into these poorly-understood areas, especially paratenic and definitive host diversity and the possible risk these parasites pose to the health to the former group. This, paired with traditional experimental infections, morphological analysis, and field surveys will lead to a greater understanding of this interesting and important nematode genus.

Fatty-binding protein and galectin of Baylisascaris schroederi: Prokaryotic expression and preliminary evaluation of serodiagnostic potential

Baylisascaris schroederi is a common parasite of captive giant pandas. The diagnosis of this ascariasis is normally carried out by a sedimentation-floatation method or PCR to detect eggs in feces, but neither method is suitable for early diagnosis. Fatty acid-binding protein (FABP) and galectin (GAL) exist in various animals and participate in important biology of parasites. Because of their good immunogenicity, they are seen as potential antigens for the diagnosis of parasitic diseases. In this study, we cloned and expressed recombinant FABP and GAL from B. schroederi (rBs-FABP and rBs-GAL) and developed indirect enzyme-linked immunosorbent assays (ELISAs) to evaluate their potential for diagnosing ascariasis in giant pandas. Immunolocalization showed that Bs-FABP and Bs-GAL were widely distributed in adult worms. The ELISA based on rBs-FABP showed sensitivity of 95.8% (23/24) and specificity of 100% (12/12), and that based on rBs-GAL had sensitivity of 91.7% (22/24) and specificity of 100% (12/12).

Characterization of the complete mitochondrial genome of Ortleppascaris sinensis (Nematoda: Heterocheilidae) and comparative mitogenomic analysis of eighteen Ascaridida nematodes

Ortleppascaris sinensis (Nematoda: Ascaridida) is a dominant intestinal nematode of the captive Chinese alligator. However, the epidemiology, molecular ecology and population genetics of this parasite remain largely unexplored. In this study, the complete mitochondrial (mt) genome sequence of O. sinensis was first determined using a polymerase chain reaction (PCR)-based primer-walking strategy, and this is also the first sequencing of the complete mitochondrial genome of a member of the genus Ortleppascaris. The circular mitochondrial genome (13,828 bp) of O. sinensis contained 12 protein-coding, 22 transfer RNA and 2 ribosomal RNA genes, but lacked the ATP synthetase subunit 8 gene. Finally, phylogenetic analysis of mtDNAs indicated that the genus Ortleppascaris should be attributed to the family Heterocheilidae. It is necessary to sequence more mtNDAs of Ortleppascaris nematodes in the future to test and confirm our conclusion. The complete mitochondrial genome sequence of O. sinensis reported here should contribute to molecular diagnosis, epidemiological investigations and ecological studies of O. sinensis and other related Ascaridida nematodes.

Molecular characterization of Baylisascaris devosi Sprent, 1952 (Ascaridoidea, Nematoda) from Kamchatka sables

The nematodes of the genus Baylisascaris are common intestinal parasites of sables (Martes (M.) zibellina kamtschadalica Birula, 1916) on the entire territory of Kamchatka peninsula. Partial sequences of Cox I mitochondrial gene were used for molecular characterization of these nematodes, which confirmed the identification based on morphological data as B. devosi Sprent, 1952. Phylogenetic relationships of this Baylisascaris species were also inferred from the ITS rDNA and LSU rDNA sequences. SEM images were provided for taxonomically important morphological features.

Partial molecular characterization of the mitochondrial genome of Baylisascaris columnaris and prevalence of infection in a wild population of Striped skunks

Members of the genus Baylisascaris utilize omnivores or carnivores as their definitive hosts. The best known member of this genus is Baylisascaris procyonis, which is an intestinal parasite of raccoons. The closest relative of B. procyonis is B. columnaris, which utilizes the common skunk as its definitive host. Although B. procyonis has been extensively studied, relatively little is known of B. columnaris. For example, the mitochondrial genome of B. procyonis has been sequenced in its entirety. Conversely, the mitochondrial genome of B. columnaris remains largely unexplored. Likewise, the prevalence of this parasite in its wild host has not been documented. In this study, we collected parasites from a wild population of skunks in the state of Utah, United States. The cytochrome c oxidase subunit 1 and 2 genes, NADH dehydrogenase 2 and several tRNA genes were sequenced from the mitochondrial genome of these parasites. We also determined the prevalence of B. columnaris in a wild population of skunks. In this work we identify several novel polymorphic genetic loci between B. procyonis and B. columnaris. These findings provide additional molecular targets for the differentiation of Baylisascaris species through clarification of genetic differences between B. columnaris and B. procyonis.

Sequence Analysis of Mitochondrial Genome of Toxascaris leonina from a South China Tiger

Toxascaris leonina is a common parasitic nematode of wild mammals and has significant impacts on the protection of rare wild animals. To analyze population genetic characteristics of T. leonina from South China tiger, its mitochondrial (mt) genome was sequenced. Its complete circular mt genome was 14,277 bp in length, including 12 protein-coding genes, 22 tRNA genes, 2 rRNA genes, and 2 non-coding regions. The nucleotide composition was biased toward A and T. The most common start codon and stop codon were TTG and TAG, and 4 genes ended with an incomplete stop codon. There were 13 intergenic regions ranging 1 to 10 bp in size. Phylogenetically, T. leonina from a South China tiger was close to canine T. leonina. This study reports for the first time a complete mt genome sequence of T. leonina from the South China tiger, and provides a scientific basis for studying the genetic diversity of nematodes between 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.

Update on Baylisascariasis, a Highly Pathogenic Zoonotic Infection

Baylisascaris procyonis, the raccoon roundworm, infects a wide range of vertebrate animals, including humans, in which it causes a particularly severe type of larva migrans. It is an important cause of severe neurologic disease (neural larva migrans [NLM]) but also causes ocular disease (OLM; diffuse unilateral subacute neuroretinitis [DUSN]), visceral larva migrans (VLM), and covert/asymptomatic infections. B. procyonis is common and widespread in raccoons, and there is increasing recognition of human disease, making a clinical consideration of baylisascariasis important. This review provides an update for this disease, especially its clinical relevance and diagnosis, and summarizes the clinical cases of human NLM and VLM known to date. Most diagnosed patients have been young children less than 2 years of age, although the number of older patients diagnosed in recent years has been increasing. The recent development of recombinant antigen-based serodiagnostic assays has aided greatly in the early diagnosis of this infection. Patients recovering with fewer severe sequelae have been reported in recent years, reinforcing the current recommendation that early treatment with albendazole and corticosteroids should be initiated at the earliest suspicion of baylisascariasis. Considering the seriousness of this zoonotic infection, greater public and medical awareness is critical for the prevention and early treatment of human cases.

Comparative analyses of the complete mitochondrial genomes of the two murine pinworms Aspiculuris tetraptera and Syphacia obvelata

Pinworms Aspiculuris tetraptera and Syphacia obvelata are important parasitic nematodes of laboratory mice, rat and other rodents. However, the mitochondrial (mt) genome of these parasites have not been known yet. In the present study, the complete mt genomes of A. tetraptera and S. obvelata were sequenced, which were 13,669 bp and 14,235 bp in size, respectively. Both genomes included 12 protein-coding genes, two rRNA genes, 22 tRNA genes and one non-coding region. The mt genomes of A. tetraptera and S. obvelata preferred bases A and T, with the highest for T and the lowest for C. The mt gene arrangements of the two pinworms were the same as that of the GA8 type. Phylogenetic analysis using mtDNA data revealed that the Bayesian inference (BI) trees contained two big branches: species from Oxyuridomorpha, Rhabditomorpha and Ascaridomorpha formed one branch, and those from Spiruromorpha formed another branch with high statistical support. The two murine pinworms A. tetraptera and S. obvelata have closer relationship than to other pinworms. This study provides a foundation for studying the population genetics, systematics and molecular phylogeny of pinworms.

Mitochondrial genomes of Heterakis gallinae and Heterakis beramporia support that they belong to the infraorder Ascaridomorpha

Heterakis gallinae and Heterakis beramporia are the most prevalent nematode infecting native chicken breed, causing major economic losses. In the present study, the complete mitochondrial genomes (mt) of H. gallinae and H. beramporia were amplified by long-PCR and then sequenced. The complete mt genomes of H. gallinae and H. beramporia were 13,973bp and 14,012bp in size, respectively. Both mt genomes contain 12 protein-coding genes, 22 transfer RNA genes and 2 ribosomal RNA genes. All genes are transcribed in the same direction and the gene arrangement is identical to Ascaridia spp. Phylogenetic analysis based on the 12 protein-coding genes revealed that the family Heterakidae (represented by H. gallinae and H. beramporia) was more closely related to the infraorder Ascaridomorpha than it was to the infraorder Oxyuridomorpha. The present study determined the complete mt genome sequences for two Heterakis species, providing useful markers for studying the systematics, population genetics, and molecular epidemiology of these Heterakis parasites.

An annotated catalogue of the ascaridoid nematode parasites of Chinese vertebrates

A catalogue, based on both examined specimens and the published literature, of all the ascaridoid nematodes recorded in China is presented. A total of 95 recognised species, representing 26 genera in five families, are reported. Detailed information on the type-host, type-locality, original reference, synonyms, annotated subsequent references of taxonomic importance, other host records, site of infection, location of type-specimens and distribution are listed for each recognised species. Additional comments on the taxonomic status of some species are also given. Moreover, some nomenclatural changes are proposed: (i) Toxascaris selenarctis Wang, 1965 and T. ailuri Wu, He & Hu, 1987 are placed in synonymy with Baylisascaris transfuga (Rudolphi, 1819); (ii) Raphidascaris lophii Wang & Wu, 1991 is a secondary homonym of R. lophii (Wu, 1949) and a replacement name, R. wangi nom. nov., is proposed for the former species; (iii) Aliascaris aetoplatea Luo, 2001 is transferred to Terranova Leiper & Atkinson, 1914, as T. aetoplatea (Luo, 2001) n. comb., and should be considered a species inquirenda; (iv) Ophidascaris orientalis (Wang, 1965) is resurrected as a valid species; (v) Phocascaris longispiculum Wang & Wu, 1991 and Ophidascaris agkistrodontis Wang, 1979 are treated as incertae sedis; and (vi) Hysterothylacium sauridae Li, Xu & Zhang, 2008 is listed as a nomen nudum.

A Review of Infectious Agents in Polar Bears (Ursus maritimus) and Their Long-Term Ecological Relevance

Disease was a listing criterion for the polar bear (Ursus maritimus) as threatened under the Endangered Species Act in 2008; it is therefore important to evaluate the current state of knowledge and identify any information gaps pertaining to diseases in polar bears. We conducted a systematic literature review focused on infectious agents and associated health impacts identified in polar bears. Overall, the majority of reports in free-ranging bears concerned serosurveys or fecal examinations with little to no information on associated health effects. In contrast, most reports documenting illness or pathology referenced captive animals and diseases caused by etiologic agents not representative of exposure opportunities in wild bears. As such, most of the available infectious disease literature has limited utility as a basis for development of future health assessment and management plans. Given that ecological change is a considerable risk facing polar bear populations, future work should focus on cumulative effects of multiple stressors that could impact polar bear population dynamics.

The complete mitochondrial genome of Pseudoterranova azarasi and comparative analysis with other anisakid nematodes

Anisakiasis/anisakidosis caused by anisakid nematodes is an emerging infectious disease that can cause a wide range of clinical syndromes and are difficult to diagnose and treat in humans. In spite of their significance as pathogens, the systematics, genetics, epidemiology and biology of these parasites remain poorly understood. In the present study, we sequenced the complete mitochondrial (mt) genome of Pseudoterranova azarasi, which is one of the most important zoonotic anisakid parasites. The circular mt genome is 13,954 bp in size and encodes of 36 genes, including 12 protein-coding, 2 ribosomal RNA and 22 transfer RNA genes. The mt gene order of P. azarasi is the same as those of Ascaris spp. (Ascarididae), Toxocara spp. (Toxocaridae) and Anisakis simplex (Anisakidae), but distinct from those of Ascaridia spp. (Ascaridiidae) and Cucullanus robustus (Cucullanidae). Phylogenetic analyses based on concatenated amino acid sequences of 12 protein-coding genes by Bayesian inference (BI) showed that Pseudoterranova were more closely related to Anisakis than they were to Contracaecum with strong a posterior probability support. This mt genome provides a novel genetic markers for exploring cryptic/sibling species and host affiliations, and should have implications for the diagnosis, prevention and control of anisakidosis in humans.

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.

Cloning and characterization of a novel sigma-like glutathione S-transferase from the giant panda parasitic nematode, Baylisascaris schroederi

Background

Baylisascaris schroederi, an intestinal nematode of the giant panda, is the cause of the often fatal disease, baylisascariasis. Glutathione S-transferases (GSTs) are versatile enzymes that can affect parasite survival and parasite-host interactions and, are therefore, potential targets for the development of diagnostic tests and vaccines.

Methods

In this study, we identified a full-length cDNA that encoded a novel, secretory sigma-like GST (Bsc-GSTσ) from a B. schroederi-omic dataset. Following cloning and sequencing, sequence and structural analyses and comparative modeling were performed using online-bioinformatics and proteomics tools. The recombinant Bsc-GSTσ (rBsc-GSTσ) protein was prokaryotically expressed and then used to detect antigenicity and reactivity using immunoblotting assays. In addition, the native protein in female adult B. schroederi was located via immunofluorescence techniques, while the preliminary ELISA-based serodiagnostic potential of rBsc-GSTσ was assessed in native and infected mouse sera.

Results

Bsc-GSTσ contained a 621-bp open reading frame that encoded a polypeptide of 206 amino acids with two typical sigma GST domain profiles, including a GST_N_Sigma_like at the N-terminus and a GST_C_Sigma_like at the C-terminus. The presence of an N-terminal signal sequence indicated that Bsc-GSTσ was a secretory protein. Sequence alignment and phylogenetic analyses showed that Bsc-GSTσ was a nematode-specific member of the Sigma class GSTs and shared the closest genetic distance with its homologue in Ascaris suum. Further comparative structure analyses indicated that Bsc-GSTσ possessed the essential structural motifs (e.g., βαβαββα) and the consensus secondary or tertiary structure that is typical for other characterized GSTσs. Immunolocalization revealed strong distributions of native Bsc-GSTσ in the body hypodermis, lateral chords, gut epithelium, gut microvilli, oviduct epithelium, and ovaries of adult female worms, similar to its homologue in A. suum. Building on good immunogenic properties, rBsc-GSTσ-based ELISA exhibited a sensitivity of 79.1% and a specificity of 82.0% to detect anti-B. schroederi IgG antibodies in the sera of experimentally infected mice.

Conclusion

This study presents a comprehensive demonstration of sequence and structural-based analysis of a new, secretory sigma-like GST from a nematode, and its good serodiagnostic performance suggests that rBsc-GSTσ has the potential to detect B. schroederi and, therefore, could be used to develop an ELISA-based serological test to diagnose baylisascariasis in giant pandas.

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.

The mitochondrial genome of Parascaris univalens--implications for a "forgotten" parasite

Background

Parascaris univalens is an ascaridoid nematode of equids. Little is known about its epidemiology and population genetics in domestic and wild horse populations. PCR-based methods are suited to support studies in these areas, provided that reliable genetic markers are used. Recent studies have shown that mitochondrial (mt) genomic markers are applicable in such methods, but no such markers have been defined for P. univalens.

Methods

Mt genome regions were amplified from total genomic DNA isolated from P. univalens eggs by long-PCR and sequenced using Illumina technology. The mt genome was assembled and annotated using an established bioinformatic pipeline. Amino acid sequences inferred from all protein-encoding genes of the mt genomes were compared with those from other ascaridoid nematodes, and concatenated sequences were subjected to phylogenetic analysis by Bayesian inference.

Results

The circular mt genome was 13,920 bp in length and contained two ribosomal RNA, 12 protein-coding and 22 transfer RNA genes, consistent with those of other ascaridoids. Phylogenetic analysis of the concatenated amino acid sequence data for the 12 mt proteins showed that P. univalens was most closely related to Ascaris lumbricoides and A. suum, to the exclusion of other ascaridoids.

Conclusions

This mt genome representing P. univalens now provides a rich source of genetic markers for future studies of the genetics and epidemiology of this parasite and its congener, P. equorum. This focus is significant, given that there is no published information on the specific prevalence and distribution of P. univalens infection in domestic and wild horse populations.

Electronic supplementary material

The online version of this article (doi:10.1186/1756-3305-7-428) contains supplementary material, which is available to authorized users.

A description of the complete mitochondrial genomes of Amphiporus formidabilis, Prosadenoporus spectaculum and Nipponnemertes punctatula (Nemertea: Hoplonemertea: Monostilifera)

We sequenced the complete mitochondrial genomes (mitogenomes) of three Hoplonemertea species, Amphiporus formidabilis, Prosadenoporus spectaculum and Nipponnemertes punctatula, which are 14,616, 14,655 and 15,354 bp in length, respectively. Each of the three circular mitogenomes consists of 37 typical genes and some non-coding regions. The nucleotide composition of the coding strand is biased toward T, almost a half of total nucleotides in these mitogenomes. There are many poly-T tracts across these mitogenomes, which exhibit T-number variation within different clones of protein-coding genes, mainly resulting from false PCR amplification. The major non-coding regions have tandem repeat motifs and hairpin-like structures that may be associated with the initiation of replication or transcription. Data published to date for nemerteans show that Palaeonemertea species usually bear the largest mitogenomes, while representatives in the more recently derived Distromatonemertea clade bear the smallest ones; and that the gene arrangement of mitogenomes seems to be variable within the phylum Nemertea, but stable within either of Heteronemertea and Hoplonemertea.

Baylisascaris potosis n. sp., a new ascarid nematode isolated from captive kinkajou, Potos flavus, from the Cooperative Republic of Guyana

We describe a new nematode species, Baylisascaris potosis n. sp., isolated from captive kinkajou, Potos flavus, from the Cooperative Republic of Guyana. The nematode was found in fecal specimens, identified morphologically, and confirmed genetically. The new species is similar to Baylisascaris procyonis, Baylisascaris columnaris, and other Baylisascaris species, but is distinguished by the position of the male phasmidial pole. Nuclear and mitochondrial DNA sequence analyses confirmed that the new species is phylogenetically distinct from all the members of the genus Baylisascaris, and groups with B. procyonis and B. columnaris. This nematode is the 10th species assigned to the genus Baylisascaris.

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.

Complete mitochondrial genomes of chimpanzee- and gibbon-derived Ascaris isolated from a zoological garden in southwest China

Roundworms (Ascaridida: Nematoda), one of the most common soil-transmitted helminths (STHs), can cause ascariasis in various hosts worldwide, ranging from wild to domestic animals and humans. Despite the veterinary and health importance of the Ascaridida species, little or no attention has been paid to roundworms infecting wild animals including non-human primates due to the current taxon sampling and survey bias in this order. Importantly, there has been considerable controversy over the years as to whether Ascaris species infecting non-human primates are the same as or distinct from Ascaris lumbricoides infecting humans. Herein, we first characterized the complete mitochondrial genomes of two representative Ascaris isolates derived from two non-human primates, namely, chimpanzees (Pan troglodytes) and gibbons (Hylobates hoolock), in a zoological garden of southwest China and compared them with those of A. lumbricoides and the congeneric Ascaris suum as well as other related species in the same order, and then used comparative mitogenomics, genome-wide nucleotide sequence identity analysis, and phylogeny to determine whether the parasites from chimpanzees and gibbons represent a single species and share genetic similarity with A. lumbricoides. Taken together, our results yielded strong statistical support for the hypothesis that the chimpanzee- and gibbon-derived Ascaris represent a single species that is genetically similar to A. lumbricoides, consistent with the results of previous morphological and molecular studies. Our finding should enhance public alertness to roundworms originating from chimpanzees and gibbons and the mtDNA data presented here also serves to enrich the resource of markers that can be used in molecular diagnostic, systematic, population genetic, and evolutionary biological studies of parasitic nematodes from either wild or domestic hosts.

The complete mitochondrial genome of Toxascaris leonina: Comparison with other closely related species and phylogenetic implications

Adults of Toxascaris leonina (Nematoda: Ascarididae) live in the gastrointestinal tract of both dogs and cats, and cause significant economic losses and potential public health problem worldwide. Although many studies have given insights into this significant pathogen, to date, the complete mitochondrial (mt) genome sequence is still not available for T. leonina. Here, we sequenced the complete mt genome of T. leonina. This AT-rich (71.53%) mt genome (14,310bp) is circular and consists of 36 genes, including 12 genes for proteins, 2 genes for rRNA and 22 genes for tRNA. All mt genes of T. leonina are transcribed in the same direction. The gene order is the same as those of Ascaris spp. (Ascarididae), Toxocara spp. (Toxocaridae), Anisakis simplex and Contracaecum rudolphii B (Anisakidae), but distinct from that of Ascaridia spp. (Ascaridiidae). Phylogenetic analyses using concatenated amino acid sequences of 12 protein-coding genes by Bayesian inference (BI) showed distinct groups with high statistical support, and our data confirm that T. leonina is a member of the Ascarididae, and that this family is more closely related to the Toxocaridae rather than the Anisakidae within the Ascaridoidea. The determination of mt genome sequences of T. leonina provides novel genetic markers for studies into the systematics, population genetics and epidemiology of this parasite.

Mitochondrial genomes of Anisakis simplex and Contracaecum osculatum (sensu stricto)--comparisons with selected nematodes

Anisakid nematodes parasitize mainly fish, marine mammals and/or fish-eating birds, and can be transmitted to a range of fish-eating mammals, including humans, where they can cause gastrointestinal disease linked to larval infection or allergic responses. In spite of the animal and human health significance of these parasites, there are still gaps in our understanding of the systematics, biology, epidemiology and ecology of anisakids. Mitochondrial (mt) DNA provides useful genetic markers for investigations in these areas, but complete mt genomic data have been lacking for most anisakids. In the present study, the mt genomes of Anisakis simplex sensu stricto and Contracaecum osculatum sensu stricto were amplified from genomic DNA by long-range polymerase chain reaction and sequenced using 454 technology. The circular mt genomes of these species were 13,926 and 13,823 bp, respectively, and each of them contained 12 protein-coding, 22 transfer RNA, and 2 ribosomal RNA genes consistent for members of the Ascaridida, Oxyurida, Spirurida, Rhabditida and Strongylida. These mt genomes provide a stepping-stone for future comparative analyses of a range of anisakids and a basis for reinvestigating their genetic relationships. In addition, these markers might be used in prospecting for cryptic species and exploring host affiliations.

Potential of recombinant inorganic pyrophosphatase antigen as a new vaccine candidate against Baylisascaris schroederi in mice

The intestinal nematode Baylisascaris schroederi is an important cause of death for wild and captive giant pandas. Inorganic pyrophosphatases (PPases) are critical for development and molting in nematode parasites and represent potential targets for vaccination. Here, a new PPase homologue, Bsc-PYP-1, from B. schroederi was identified and characterized, and its potential as a vaccine candidate was evaluated in a mouse challenge model. Sequence alignment of PPases from nematode parasites and other organisms show that Bsc-PYP-1 is a nematode-specific member of the family I soluble PPases. Immunohistochemistry revealed strong localization of native Bsc-PYP-1 to the body wall, gut epithelium, ovary and uterus of adult female worms. Additionally, Bsc-PYP-1 homologues were found in roundworms infecting humans (Ascaris lumbricoides), swine (Ascaris suum) and dogs (Toxocara canis). In two vaccine trials, recombinant Bsc-PYP-1 (rBsc-PYP-1) formulated with Freund complete adjuvant induced significantly high antigen-specific immunoglobulin (Ig)G but no IgE or IgM responses. Analysis of IgG-subclass profiles revealed a greater increase of IgG1 than IgG2a. Splenocytes from rBsc-PYP-1/FCA-immunized mice secreted low levels of T helper (Th)1-type cytokines, interferon-γ and interleukin (IL)-2, while producing significantly high levels of IL-10 and significantly elevated levels of IL-4 (Th2 cytokines) after stimulation with rBsc-PYP-1 in vitro. Finally, vaccinated mice had 69.02–71.15% reductions (in 2 experiments) in larval recovery 7 days post-challenge (dpc) and 80% survival at 80 dpc. These results suggest that Th2-mediated immunity elicited by rBsc-PYP-1 provides protection against B. schroederi, and the findings should contribute to further development of Bsc-PYP-1 as a candidate vaccine against baylisascariasis.