Search

Abstract Background Haemaphysalis qinghaiensis is a tick species distributed only in China. Due to its ability to transmit a variety of pathogens, including species of the genera Anaplasma, Rickettsia, Babesia, and Theileria, it seriously endangers livestock husbandry. However, the microbial community of the midgut of H. qinghaiensis females collected from yaks and Tibetan sheep has not yet been characterized using metagenomic sequencing technology. Methods Haemaphysalis qinghaiensis were collected from the skins of yaks and Tibetan sheep in Gansu Province, China. Genomic DNA was extracted from the midguts and midgut contents of fully engorged H. qinghaiensis females collected from the two hosts. Metagenomic sequencing technology was used to analyze the microbial community of the two groups. Results Fifty-seven phyla, 483 genera, and 755 species were identified in the two groups of samples. The ticks from the two hosts harbored common and unique microorganisms. At the phylum level, the dominant common phyla were Proteobacteria, Firmicutes, and Mucoromycota. At the genus level, the dominant common genera were Anaplasma, Ehrlichia, and Pseudomonas. At the species level, bacteria including Anaplasma phagocytophilum, Ehrlichia minasensis, and Pseudomonas aeruginosa along with eukaryotes such as Synchytrium endobioticum and Rhizophagus irregularis, and viruses such as the orf virus, Alphadintovirus mayetiola, and Parasteatoda house spider adintovirus were detected in both groups. In addition, the midgut of H. qinghaiensis collected from yaks had unique microbial taxa including two phyla, eight genera, and 23 species. Unique microorganisms in the midgut of H. qinghaiensis collected from Tibetan sheep included two phyla, 14 genera, and 32 species. Kyoto Encyclopedia of Genes and Genomes enrichment analysis revealed that the functional genes of the microbiome of H. qinghaiensis were annotated to six pathways, and the metabolic pathways included 11 metabolic processes, in which the genes involved in carbohydrate metabolism were the most abundant, followed by the genes involved in lipid metabolism. Conclusions These findings indicate that most of the microbial species in the collected H. qinghaiensis ticks were the same in both hosts, but there were also slight differences. The analytical data from this study have enhanced our understanding of the midgut microbial composition of H. qinghaiensis collected from different hosts. The database of H. qinghaiensis microbe constructed from this study will lay the foundation for predicting tick-borne diseases. Furthermore, a comprehensive understanding of tick microbiomes will be useful for understanding vector competency and interactions with ticks and midgut microorganisms. Graphical abstract

Abstract Background Toxoplasma gondii and Neospora caninum are closely related protozoan parasites that are considered important causes of abortion in livestock, causing huge economic losses. Hunan Province ranks 12th in the production of beef and mutton in China. However, limited data are available on the seroprevalence, risk factors and molecular characterization of T. gondii and N. caninum in beef cattle and goats in Hunan province, China. Methods Sera of 985 beef cattle and 1147 goats were examined for the presence of specific antibodies against T. gondii using indirect hemagglutination test (IHAT) and anti-N. caninum IgG using competitive-inhibition enzyme-linked immunoassay assay (cELISA). Statistical analysis of possible risk factors was performed using PASW Statistics. Muscle samples of 160 beef cattle and 160 goats were examined for the presence of T. gondii DNA (B1 gene) and N. caninum DNA (Nc-5 gene) by nested PCR. The B1 gene-positive samples were genotyped at 10 genetic markers using the multilocus nested PCR-RFLP (Mn-PCR-RFLP). Results Specific IgG against T. gondii were detected in 8.3% (82/985) and 13.3% (153/1147) and against N. caninum in 2.1% (21/985) and 2.0% (23/1147) of the beef cattle and goats, respectively. Based on statistical analysis, the presence of cats, semi-intensive management mode and gender were identified as significant risk factors for T. gondii infection in beef cattle. Age was a significant risk factor for T. gondii infection in goats (P 3 years was a significant risk factor for N. caninum infection in beef cattle (P

Abstract Background Increasing evidence has shown that non-coding RNA (ncRNA) molecules play fundamental roles in cells, and many are stable in body fluids as circulating RNAs. Study on these ncRNAs will provide insights into toxoplasmosis pathophysiology and/or help reveal diagnostic biomarkers. Methods We performed a high-throughput RNA-Seq study to comprehensively profile the microRNAs (miRNAs) and PIWI-interacting RNAs (piRNAs) in rabbit serum and urine after infection with Toxoplasma gondii oocysts during the whole infection process. Results Total RNA extracted from serum and urine samples of acutely infected [8 days post-infection (DPI)], chronically infected (70 DPI) and uninfected rabbits were subjected to genome-wide small RNA sequencing. We identified 2089 miRNAs and 2224 novel piRNAs from the rabbit sera associated with T. gondii infection. Meanwhile, a total of 518 miRNAs and 4182 novel piRNAs were identified in the rabbit urine associated with T. gondii infection. Of these identified small ncRNAs, 1178 and 1317 serum miRNAs and 311 and 294 urine miRNAs were identified as differentially expressed (DE) miRNAs in the acute and chronic stages of infections, respectively. A total of 1748 and 1814 serum piRNAs and 597 and 708 urine piRNAs were found in the acute and chronic infection stages, respectively. Of these dysregulated ncRNAs, a total of 88 common DE miRNAs and 120 DE novel piRNAs were found in both serum and urine samples of infected rabbits. Conclusions These findings provide valuable data for revealing the physiology of herbivore toxoplasmosis caused by oocyst infection. Circulating ncRNAs identified in this study are potential novel diagnostic biomarkers for the detection/diagnosis of toxoplasmosis in herbivorous animals. Graphical Abstract

Abstract Background The family Capillariidae is a group of thread-like nematodes of 27 genera and over 300 species that infect a great variety of hosts including humans. Among these, some taxa such as the genus Aonchotheca have remained controversial regarding their systematic status for decades. The aim of the current study was to verify Aonchotheca’s systemic status and to further determine whether it is a distinct genus from Capillaria using molecular and phylogenetic analyses. Results We sequenced the mitochondrial (mt) genome and nuclear small subunit (18S) rRNA gene of Aonchotheca putorii, a representative species of the genus, and investigated its systematic status in Trichinellida using maximum likelihood and Bayesian inference. The differences in amino acid sequences of 13 protein-coding genes were 12.69–67.35% among Aonchotheca, Capillaria, Eucoleus, and Pseudocapillaria with cox1 (12.69%) and atp8 (67.35%) as the most and the least conserved gene, respectively, and the difference of two mt rRNAs was 18.61–34.15%. Phylogenetic analyses of the complete mt genome and 18S rRNAs unequivocally showed that Aonchotheca was a distinct genus from Capillaria. Conclusions Large difference exists among Aonchotheca, Capillaria, Eucoleus, and Pseudocapillarias. Aonchotheca putorii is the first species in the genus Aonchotheca for which a complete mitogenome has been sequenced. These data are useful for phylogenetics, systematics and the evolution of Capillariidae. Graphical Abstract

Abstract Background Haemaphysalis flava is a hematophagous ectoparasite that acquires the nutrition needed for development and reproduction by sucking blood and digesting the blood meal. During blood-sucking and blood-meal digestion, the prevention of blood coagulation is important for this tick. Previous studies have shown that heat shock cognate 70 (HSC70) protein has certain anticoagulant activities, but its immunogenicity remains unclear. Also, whether the mutation of individual bases of the TKD-like peptide of HSC70 through the overlap extension method can change its anticoagulant activities and immunogenicity remains to be investigated. Methods The gene encoding the HSC70 protein was cloned from a complementary DNA library synthesized from H. flava. The coding gene of the TKD-like peptide of HSC70 was mutated into a TKD peptide coding gene (HSC70TKD) using the overlap extension method. Escherichia coli prokaryotic expression plasmids were constructed to obtain the recombinant proteins of HSC70 (rHSC70) and HSC70TKD (rHSC70TKD). The purified rHSC70 and rHSC70TKD were evaluated at different concentrations for anticoagulant activities using four in vitro clotting assays. Emulsifying recombinant proteins with complete and incomplete Freund’s adjuvants were subcutaneously immunized in Sprague Dawley rats. The serum antibody titers and serum concentrations of interferon-gamma (IFN-γ) and interleukin-4 (IL-4) were detected using an indirect enzyme-linked immunosorbent assay to assess the immunogenicity of rHSC70 and rHSC70TKD. Results The open reading frame of HSC70 was successfully amplified and found to have a length of 1958 bp. The gene encoding the TKD-like peptide of HSC70 was artificially mutated, with the 1373-position adenine (A) of the original sequence mutated into guanine (G), the 1385-position cytosine (C) mutated into G and the 1386-position G mutated into C. rHSC70 and rHSC70TKD that fused with His-tag were obtained using the expression plasmids pET-28a-HSC70 and pET-28a-HSC70TKD, respectively. rHSC70 and rHSC70TKD prolonged the thrombin time (TT) and reduced the fibrinogen (FIB) content in the plasma, but did not affect the prothrombin time (PT) or activated partial thromboplastin time (APTT) when compared to the negative control. Interestingly, the ability of rHSC70TKD to prolong the TT and reduce the FIB content in the plasma was better than that of rHSC70. The specific antibody titers of both rHSC70 and rHSC70TKD in rat serum reached 1:124,000 14 days after the third immunization. The serum concentration of IFN-γ in the rHSC70TKD group was higher than that in the rHSC70 group. The rHSC70 group has the highest serum concentration of IL-4, and the serum concentration of IL-4 in the rHSC70TKD group was higher than that in the negative group. Conclusions rHSC70 and rHSC70TKD exhibited anticoagulant activities by prolonging the TT and reducing the FIB content in vitro. rHSC70TKD had better anticoagulant activities than rHSC70. Both rHSC70 and rHSC70TKD had good immunogenicity and induced humoral and cellular immunity.

Abstract Background Fleas are the most economically significant blood-feeding ectoparasites worldwide. Ctenocephalides felis and Pulex irritans can parasitize various animals closely related to humans and are of high veterinary significance. Methods In this study, 82 samples were collected from 7 provinces of China. Through studying the nuclear genes ITS1 and EF-1α and two different mitochondrial genes cox1 and cox2, the population genetics and genetic variation of C. felis and P. irritans in China were further investigated. Results The intraspecies differences between C. felis and P. irritans ranged from 0 to 3.9%. The interspecific variance in the EF-1α, cox1, and cox2 sequences was 8.2–18.3%, while the ITS1 sequence was 50.1–52.2%. High genetic diversity was observed in both C. felis and P. irritans, and the nucleotide diversity of cox1 was higher than that of cox2. Moderate gene flow was detected in the C. felis and P. irritans populations. Both species possessed many haplotypes, but the haplotype distribution was uneven. Fu's Fs and Tajima's D tests showed that C. felis and P. irritans experienced a bottleneck effect in Guangxi Zhuang Autonomous Region and Henan province. Evolutionary analysis suggested that C. felis may have two geographical lineages in China, while no multiple lineages of P.irritans were found. Conclusions Using sequence comparison and the construction of phylogenetic trees, we found a moderate amount of gene flow in the C. felis and P. irritans populations. Both species possessed many haplotypes, but the distribution of haplotypes varied among the provinces. Fu’s Fs and Tajima’s D tests indicated that both species had experienced a bottleneck effect in Guangxi and Henan provinces. Evolutionary analysis suggested that C. felis may have two geographical lineages in China, while no multiple lineages of P.irritans were found. This study will help better understand fleas' population genetics and evolutionary biology. Graphical Abstract

Abstract Background Fleas (Insecta: Siphonaptera) are obligatory hematophagous ectoparasites of humans and animals and serve as vectors of many disease-causing agents. Despite past and current research efforts on fleas due to their medical and veterinary importance, correct identification and robust phylogenetic analysis of these ectoparasites have often proved challenging. Methods We decoded the complete mitochondrial (mt) genome of the human flea Pulex irritans and nearly complete mt genome of the dog flea Ctenocephalides canis, and subsequently used this information to reconstruct the phylogeny of fleas among Endopterygota insects. Results The complete mt genome of P. irritans was 20,337 bp, whereas the clearly sequenced coding region of the C. canis mt genome was 15,609 bp. Both mt genomes were found to contain 37 genes, including 13 protein-coding genes, 22 transfer RNA genes and two ribosomal RNA genes. The coding region of the C. canis mt genome was only 93.5% identical to that of the cat flea C. felis, unequivocally confirming that they are distinct species. Our phylogenomic analyses of the mt genomes showed a sister relationship between the order Siphonaptera and orders Diptera + Mecoptera + Megaloptera + Neuroptera and positively support the hypothesis that the fleas in the order Siphonaptera are monophyletic. Conclusions Our results demonstrate that the mt genomes of P. irritans and C. canis are different. The phylogenetic tree shows that fleas are monophyletic and strongly support an order-level objective. These mt genomes provide novel molecular markers for studying the taxonomy and phylogeny of fleas in the future.

Abstract Background Ticks, which are ectoparasites of animals, may carry multiple pathogens. The cattle tick Rhipicephalus microplus is an important bovine parasite in China. However, the midgut microbiome of R. microplus from China has not been characterized via metagenomic methods. Methods Rhipicephalus microplus were collected from cattle in the city of Changsha in Hunan province, China. The DNA of the midgut contents was extracted from fully engorged adult female R. microplus. A DNA library was constructed and sequenced using an Illumina HiSeq sequencing platform. SOAPdenovo software was used to assemble and analyze the clean data. The latent class analysis algorithm applied to system classification by MEGAN software was used to annotate the information on the species’ sequences. DIAMOND software was used to compare unigenes with the Kyoto Encyclopedia of Genes and Genomes (KEGG) database, and functional annotation was carried out based on the results of the comparison. Results The dominant phyla in the five samples were Firmicutes, Proteobacteria, and Actinobacteria. Streptococcus, Mycobacterium, Anaplasma, Enterococcus, Shigella, Lactobacillus, Brachyspira, Pseudomonas, Enterobacter, Bacillus, and Lactococcus were the dominant genera in the five samples. The endosymbiotic bacterium Wolbachia was also detected in all of the samples. Mycobacterium malmesburyense, Streptococcus pneumoniae, Anaplasma phagocytophilum, Enterococcus faecium, Shigella sonnei, Enterococcus faecalis, Lactobacillus casei, Brachyspira hampsonii, Pseudomonas syringae, Enterobacter cloacae, and Lactococcus garvieae were the dominant species in the five samples. In addition to these bacterial species, we also detected some eukaryotes, such as Rhizophagus irregularis, Enterospora canceri, Smittium culicis, Zancudomyces culisetae, Trachipleistophora hominis, and viruses such as orf virus, human endogenous retrovirus type W, enzootic nasal tumor virus of goats, bovine retrovirus CH15, and galidia endogenous retrovirus in all of the samples at the species level. The results of the annotated KEGG pathway predictions for the gene functions of the midgut microflora of R. microplus indicated genes involved in lipid and amino acid metabolism, infectious diseases (e.g., Streptococcus pneumonia infection, human granulocytic anaplasmosis, Shigella sonnei infection, Salmonella enterica infection, and pathogenic Escherichia coli infection), and cancer. Conclusions Our study revealed that the midgut microbiome of R. microplus is not only composed of a large number of bacteria, but that a portion also comprises eukaryotes and viruses. The data presented here enhance our understanding of this tick’s midgut microbiome and provide fundamental information for the control of ticks and tick-borne diseases. Graphical Abstract

Abstract Background Fragmented mitochondrial (mt) genomes and extensive mt gene rearrangements have been frequently reported from parasitic lice (Insecta: Phthiraptera). However, relatively little is known about the mt genomes from the family Philopteridae, the most species-rich family within the suborder Ischnocera. Methods Herein, we use next-generation sequencing to decode the mt genome of Falcolipeurus suturalis and compare it with the mt genome of F. quadripustulatus. Phylogenetic relationships within the family Philopteridae were inferred from the concatenated 13 protein-coding genes of the two Falcolipeurus lice and members of the family Philopteridae using Bayesian inference (BI) and maximum likelihood (ML) methods. Results The complete mt genome of F. suturalis is a circular, double-stranded DNA molecule 16,659 bp in size that contains 13 protein-coding genes, 22 transfer RNA genes, two ribosomal RNA genes, and three non-coding regions. The gene order of the F. suturalis mt genome is rearranged relative to that of F. quadripustulatus, and is radically different from both other louse species and the putative ancestral insect. Phylogenetic analyses revealed clear genetic distinctiveness between F. suturalis and F. quadripustulatus (Bayesian posterior probabilities = 1.0 and bootstrapping frequencies = 100), and that the genus Falcolipeurus is sister to the genus Ibidoecus (Bayesian posterior probabilities = 1.0 and bootstrapping frequencies = 100). Conclusions These datasets help to better understand gene rearrangements in lice and the phylogenetic position of Falcolipeurus and provide useful genetic markers for systematic studies of bird lice. Graphic abstract

Abstract Background Human gnathostomiasis is a food-borne zoonosis. Its etiological agents are the third-stage larvae of Gnathostoma spp. Human gnathostomiasis is often reported in developing countries, but it is also an emerging disease in developed countries in non-endemic areas. The recent surge in cases of human gnathostomiasis is mainly due to the increasing consumption of raw freshwater fish, amphibians, and reptiles. Methods This article reviews the literature on Gnathostoma spp. and the disease that these parasites cause in humans. We review the literature on the life cycle and pathogenesis of these parasites, the clinical features, epidemiology, diagnosis, treatment, control, and new molecular findings on human gnathostomiasis, and social-ecological factors related to the transmission of this disease. Conclusions The information presented provides an impetus for studying the parasite biology and host immunity. It is urgently needed to develop a quick and sensitive diagnosis and to develop an effective regimen for the management and control of human gnathostomiasis.

Abstract Background The family Hoplopleuridae contains at least 183 species of blood-sucking lice, which widely parasitize both mice and rats. Fragmented mitochondrial (mt) genomes have been reported in two rat lice (Hoplopleura kitti and H. akanezumi) from this family, but some minichromosomes were unidentified in their mt genomes. Methods We sequenced the mt genome of the rat louse Hoplopleura sp. with an Illumina platform and compared its mt genome organization with H. kitti and H. akanezumi. Results Fragmented mt genome of the rat louse Hoplopleura sp. contains 37 genes which are on 12 circular mt minichromosomes. Each mt minichromosome is 1.8–2.7 kb long and contains 1–5 genes and one large non-coding region. The gene content and arrangement of mt minichromosomes of Hoplopleura sp. (n = 3) and H. kitti (n = 3) are different from those in H. akanezumi (n = 3). Phylogenetic analyses based on the deduced amino acid sequences of the eight protein-coding genes showed that the Hoplopleura sp. was more closely related to H. akanezumi than to H. kitti, and then they formed a monophyletic group. Conclusions Comparison among the three rat lice revealed variation in the composition of mt minichromosomes within the genus Hoplopleura. Hoplopleura sp. is the first species from the family Hoplopleuridae for which a complete fragmented mt genome has been sequenced. The new data provide useful genetic markers for studying the population genetics, molecular systematics and phylogenetics of blood-sucking lice.

Abstract Background Cryptosporidium viatorum is a minor Cryptosporidium pathogen in humans. Currently, there is limited information regarding the prevalence and genotypes of C. viatorum in animals in China. Methods In this study, 228 faecal samples were collected from two wild rat species (Leopoldamys edwardsi and Berylmys bowersi) in Chongqing Municipality and Guangdong Province, China. These specimens were analyzed for C. viatorum and then subtyped it using PCR and sequence analysis of the small subunit ribosomal RNA (SSU rRNA) and 60-kilodalton glycoprotein (gp60) genes, respectively. Results A total of 25 (11.0%) faecal samples were tested positive for C. viatorum by SSU rRNA assay. Of these samples, 4 (3.6%) came from L. edwardsi and 21 (18.0%) from B. bowersi. Of the 25 C. viatorum-positive samples, 17 were successfully amplified at the gp60 gene locus, which represented four subtypes belonging to two subtype families, including XVa (XVaA6, XVaA3g, XVaA3h) and XVc (XVcA2G1). Phylogenetic analysis based on the gp60 amino acid sequences indicated that all of the C. viatorum isolates grouped together, supporting the conclusion that C. viatorum from the wild rats represent two subtype families. Conclusions These results indicate an occurrence of C. viatorum XVa subtype family from rats which is genetically identical to those found in humans. Our findings suggest that wild rats may be a potential source of human cryptosporidiosis.

Abstract Background Toxocara canis, a globally distributed roundworm, can cause debilitating disease in dogs and humans; however, little is known about the metabolomic response of the hosts to T. canis infection. There is an increasing need to understand the metabolic mechanisms underlying the pathogenesis of T. canis infection in dogs. Here, we examined the metabolomic changes in Beagle dogsʼ serum following T. canis infection using LC-MS/MS. Results The metabolic profiles of Beagle dogsʼ serum were determined at 12 h, 24 h, 10 d and 36 d after oral infection with 300 infectious T. canis eggs by LC-MS/MS. We tested whether the T. canis-associated differentially abundant metabolites could distinguish the serum of infected dogs from controls, as measured by the area under the receiver operating characteristic (ROC) curve (AUC). The differentially expressed metabolites were further evaluated by principal components analysis and pathway enrichment analysis. A total of 5756 and 5299 ions were detected in ESI+ and ESI− mode, respectively. ROC curve analysis revealed nine and five metabolite markers, at 12 hpi and 24 hpi to 36 dpi, respectively, with potential diagnostic value for toxocariasis. The levels of taurocholate, estradiol, prostaglandins and leukotriene were significantly changed. Primary bile acid biosynthesis pathway, steroid hormone biosynthesis pathway and biosynthesis of unsaturated fatty acids pathway were significantly altered by T. canis infection. Conclusions These findings show that T. canis infection can induce several changes in the dog serum metabolome and that the metabolic signature associated with T. canis infection in dogs has potential for toxocariasis diagnosis.

Abstract 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.

Abstract Background Ticks and tick-borne diseases are of major public health concern. Currently, development of vaccines against ticks is considered crucial for their control. A critical step in this process is the screening of viable antigens. Faeces are byproducts of digestion and blood meal utilization, and partly reflect the vitality and vector potential of ticks. However, an integrated analysis of proteins in tick faeces is lacking. The present study explored the protein components in the faeces of the tick Haemaphysalis flava, by liquid chromatography–tandem mass spectrometry (LC/MS-MS) to identify potential protein antigens for vaccine development against ticks. Methods Faeces from adult H. flava engorged females were collected. Proteins were extracted from faeces, and the trypsin-digested peptides were analyzed by LC/MS-MS. High confidence proteins were identified based on unique peptides revealed by MS. Potential faecal protein genes, as well as their sources, were also characterized by searching previous transcriptome datasets from the salivary glands and midgut of H. flava. Results In total, 21 were recognized with confidence. Amongst these, 18 were of likely tick origin, while three proteins (serum albumin, haemoglobin α and β subunits) were likely from hosts. Seventeen unigenes corresponding to these proteins were retrieved by searching our previous H. flava salivary glands and midgut transcriptomic datasets. Some proteins were reported to prevent blood clotting, play a role in immunity and antibiosis, and formation of musculature. The functions of the remaining proteins are unknown. Conclusions Identifying antigens for tick vaccine development is feasible by analyzing the faecal proteome as well as the transcriptomes of salivary glands and midguts. The vast number of proteins detected in tick faeces highlights the complexity of blood digestion in ticks, a field that needs more investigation.

Abstract Background Melophagus ovinus, one of the most common haematophagous ectoparasites of sheep, can cause anaemia and reductions in weight gain, wool growth and hide value. However, no information is available about the microfloral structure of the midgut of this ectoparasite. In the present study, we investigated the microbial community structure of the midgut contents of fully engorged female and male M. ovinus using Illumina HiSeq. Results The phylum showing the highest abundance was Proteobacteria (99.9%). The dominant bacterial genera in females and males were Bartonella, Arsenophonus and Wolbachia. Some less abundant bacterial genera were also detected, including Enterobacter, Acinetobacter, Halomonas, Shewanella, Bacillus and Staphylococcus. Conclusions Bartonella, Arsenophonus and Wolbachia were the dominant bacterial genera in the midgut of female and male M. ovinus. Although detected, Enterobacter, Acinetobacter, Halomonas, Shewanella, Bacillus and Staphylococcus showed low abundances. Importantly, this is the first report of the presence of Arsenophonus, Wolbachia, Enterobacter, Halomonas, Shewanella, Bacillus and Staphylococcus in the midgut of M. ovinus.

Abstract Background Rumen flukes parasitize the rumen and reticulum of ruminants, causing paramphistomiasis. Over the years, there has been considerable debate as to whether Paramphistomum leydeni and Paramphistomum cervi are the same or distant species. Methods In the present study, the complete mitochondrial (mt) genome of P. leydeni was amplified using PCR-based sequencing and compared with that of P. cervi. The second internal transcribed spacer (ITS-2) of nuclear ribosomal DNA (rDNA) of P. leydeni specimens (n = 6) and P. cervi specimens (n = 8) was amplified and then sequenced. Phylogenetic relationship of the concatenated amino acid sequence data for 12 protein-coding genes of the two rumen flukes and selected members of Trematoda was evaluated using Bayesian inference (BI). Results The complete mt genome of P. leydeni was 14,050 bp in size. Significant nucleotide difference between the P. leydeni mt genome and that of P. cervi (14.7%) was observed. For genetic divergence in ITS-2, sequence difference between P. leydeni and P. cervi was 3.1%, while no sequence variation was detected within each of them. Phylogenetic analysis indicated that P. leydeni and P. cervi are closely-related but distinct rumen flukes. Conclusions Results of the present study support the proposal that P. leydeni and P. cervi represent two distinct valid species. The mt genome sequences of P. leydeni provide plentiful resources of mitochondrial markers, which can be combined with nuclear markers, for further comparative studies of the biology of P. leydeni and its congeners from China and other countries.

Abstract Background Gongylonema pulchrum (Nematoda: Gongylonematidae), a thread-like spirurid gullet worm, infects a range of mammalian definitive hosts, including cattle, pigs, equines, goats, primates and humans, and can cause gongylonemiasis. Methods In the present study, the complete mitochondrial (mt) genome of G. pulchrum was obtained using Long-range PCR and subsequent primer walking. The phylogenetic position of G. pulchrum within the Spiruromorpha was established using Bayesian analyses of the protein-coding genes at the amino acid level. Results The length of this AT-rich (75.94%) mt genome is 13,798 bp. It contains 12 protein-coding genes, two ribosomal RNA genes, 22 transfer RNA genes and one non-coding region. The gene arrangement is the same as those of Thelazia callipaeda (Thelaziidae) and Setaria digitata (Onchocercidae), but distinct from that of Heliconema longissimum (Physalopteridae). Phylogenetic analyses, based on the concatenated amino acid sequence data for all 12 protein-coding genes using Bayesian inference (BI) method, showed that G. pulchrum (Gongylonematidae) was more closely related to Spirocerca lupi (Spiruroidea) than other members of the infraorder Spiruromorpha. Conclusions The present study represents the first mt genome sequence for the family Gongylonematidae, which provides the opportunity to develop novel genetic markers for studies of epidemiology, population genetics and systematics of this nematode of human and animal health significance.

Abstract 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.

Abstract Background Limited available sequence information has greatly impeded population genetics, phylogenetics and systematics studies in the subclass Acari (mites and ticks). Mitochondrial (mt) DNA is well known to provide genetic markers for investigations in these areas, but complete mt genomic data have been lacking for many Acari species. Herein, we present the complete mt genome of the scab mite Psoroptes cuniculi. Methods P. cuniculi was collected from a naturally infected New Zealand white rabbit from China and identified by morphological criteria. The complete mt genome of P. cuniculi was amplified by PCR and then sequenced. The relationships of this scab mite with selected members of the Acari were assessed by phylogenetic analysis of concatenated amino acid sequence datasets by Bayesian inference (BI), maximum likelihood (ML) and maximum parsimony (MP). Results This mt genome (14,247 bp) is circular and consists of 37 genes, including 13 genes for proteins, 22 genes for tRNA, 2 genes for rRNA. The gene arrangement in mt genome of P. cuniculi is the same as those of Dermatophagoides farinae (Pyroglyphidae) and Aleuroglyphus ovatus (Acaridae), but distinct from those of Steganacarus magnus (Steganacaridae) and Panonychus citri (Tetranychidae). Phylogenetic analyses using concatenated amino acid sequences of 12 protein-coding genes, with three different computational algorithms (BI, ML and MP), showed the division of subclass Acari into two superorders, supported the monophylies of the both superorders Parasitiformes and Acariformes; and the three orders Ixodida and Mesostigmata and Astigmata, but rejected the monophyly of the order Prostigmata. Conclusions The mt genome of P. cuniculi represents the first mt genome of any member of the family Psoroptidae. Analysis of mt genome sequences in the present study has provided new insights into the phylogenetic relationships among several major lineages of Acari species.

Abstract Background Fascioliasis is an important and neglected disease of humans and other mammals, caused by trematodes of the genus Fasciola. Fasciola hepatica and F. gigantica are valid species that infect humans and animals, but the specific status of Fasciola sp. (‘intermediate form’) is unclear. Methods Single specimens inferred to represent Fasciola sp. (‘intermediate form’; Heilongjiang) and F. gigantica (Guangxi) from China were genetically identified and characterized using PCR-based sequencing of the first and second internal transcribed spacer regions of nuclear ribosomal DNA. The complete mitochondrial (mt) genomes of these representative specimens were then sequenced. The relationships of these specimens with selected members of the Trematoda were assessed by phylogenetic analysis of concatenated amino acid sequence datasets by Bayesian inference (BI). Results The complete mt genomes of representatives of Fasciola sp. and F. gigantica were 14,453 bp and 14,478 bp in size, respectively. Both mt genomes contain 12 protein-coding genes, 22 transfer RNA genes and two ribosomal RNA genes, but lack an atp 8 gene. All protein-coding genes are transcribed in the same direction, and the gene order in both mt genomes is the same as that published for F. hepatica. Phylogenetic analysis of the concatenated amino acid sequence data for all 12 protein-coding genes showed that the specimen of Fasciola sp. was more closely related to F. gigantica than to F. hepatica. Conclusions The mt genomes characterized here provide a rich source of markers, which can be used in combination with nuclear markers and imaging techniques, for future comparative studies of the biology of Fasciola sp. from China and other countries.

Abstract 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 (rrn S) and the least conserved nad 2 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.

Background: Haemaphysalis flava is a hematophagous ectoparasite that acquires the nutrition needed for development and reproduction by sucking blood and digesting the blood meal. During blood-sucking and blood-meal digestion, the prevention of blood coagulation is important for this tick. Previous studies have shown that heat shock cognate 70 (HSC70) protein has certain anticoagulant activities, but its immunogenicity remains unclear. Also, whether the mutation of individual bases of the TKD-like peptide of HSC70 through the overlap extension method can change its anticoagulant activities and immunogenicity remains to be investigated. Methods: The gene encoding the HSC70 protein was cloned from a complementary DNA library synthesized from H. flava. The coding gene of the TKD-like peptide of HSC70 was mutated into a TKD peptide coding gene (HSC70TKD) using the overlap extension method. Escherichia coli prokaryotic expression plasmids were constructed to obtain the recombinant proteins of HSC70 (rHSC70) and HSC70TKD (rHSC70TKD). The purified rHSC70 and rHSC70TKD were evaluated at different concentrations for anticoagulant activities using four in vitro clotting assays. Emulsifying recombinant proteins with complete and incomplete Freund's adjuvants were subcutaneously immunized in Sprague Dawley rats. The serum antibody titers and serum concentrations of interferon-gamma (IFN-γ) and interleukin-4 (IL-4) were detected using an indirect enzyme-linked immunosorbent assay to assess the immunogenicity of rHSC70 and rHSC70TKD. Results: The open reading frame of HSC70 was successfully amplified and found to have a length of 1958 bp. The gene encoding the TKD-like peptide of HSC70 was artificially mutated, with the 1373-position adenine (A) of the original sequence mutated into guanine (G), the 1385-position cytosine (C) mutated into G and the 1386-position G mutated into C. rHSC70 and rHSC70TKD that fused with His-tag were obtained using the expression plasmids pET-28a-HSC70 and pET-28a-HSC70TKD, respectively. rHSC70 and rHSC70TKD prolonged the thrombin time (TT) and reduced the fibrinogen (FIB) content in the plasma, but did not affect the prothrombin time (PT) or activated partial thromboplastin time (APTT) when compared to the negative control. Interestingly, the ability of rHSC70TKD to prolong the TT and reduce the FIB content in the plasma was better than that of rHSC70. The specific antibody titers of both rHSC70 and rHSC70TKD in rat serum reached 1:124,000 14 days after the third immunization. The serum concentration of IFN-γ in the rHSC70TKD group was higher than that in the rHSC70 group. The rHSC70 group has the highest serum concentration of IL-4, and the serum concentration of IL-4 in the rHSC70TKD group was higher than that in the negative group. Conclusions: rHSC70 and rHSC70TKD exhibited anticoagulant activities by prolonging the TT and reducing the FIB content in vitro. rHSC70TKD had better anticoagulant activities than rHSC70. Both rHSC70 and rHSC70TKD had good immunogenicity and induced humoral and cellular immunity. [ABSTRACT FROM AUTHOR]

Abstract Background Toxoplasma gondii is a highly prevalent protozoan parasite infecting a wide range of animals and humans. The epidemiological and biological diversity of T. gondii has resulted in a high genetic variation and unusual population structure in this parasite. This study examined sequence diversity in dense granule 5 (GRA5) gene among T. gondii isolates from different hosts and geographical regions. Methods The entire genome region of the GRA5 gene was amplified and sequenced from 14 T. gondii isolates, and phylogenetic relationship among these T. gondii isolates was reconstructed using Bayesian inference (BI) and maximum parsimony (MP) based on the GRA5 sequences. Results The complete sequence of the GRA5 gene was 1614 bp in length for strains TgCatBr5 and MAS, but 1617 bp for the other 12 strains. Sequence analysis identified 41 (0–1.7%) variable nucleotide positions among all isolates, with 18 variations of these being in the coding region. Variable positions in the coding region resulted in 11 amino acid substitutions, and a deletion of 3 bp in the strains TgCatBr5 and MAS leading to the deletion of one amino acid. Sequence variations resulted in the existence of polymorphic restriction sites for endonucleases AatII and MluI, allowing the differentiation of the three major clonal lineage types I, II and III by PCR-RFLP. Phylogenetic analyses using BI and MP supported the clear differentiation of the examined T. gondii strains into their respective genotypes. Conclusions This study demonstrated the existence of sequence variability in the GRA5 gene sequence among T. gondii isolates from different hosts and geographical regions, which allowed the differentiation of the examined T. gondii strains into their respective genotypes, suggesting that this highly polymorphic GRA5 locus may provide a new genetic marker for population genetic studies of T. gondii isolates.