Trx4, a novel thioredoxin protein, is important for Toxoplasma gondii fitness

BACKGROUND

To successfully replicate within the host cell, Toxoplasma gondii employs several mechanisms to overcome the host cell defenses and mitigate the harmful effects of the free radicals resulting from its own metabolic processes using effectors such as thioredoxin proteins. In this study, we characterize the location and functions of a newly identified thioredoxin in T. gondii, which was named Trx4.

METHODS

We characterized the functional role of Trx4 in T. gondii Type I RH and Type II Pru strains by gene knockout and studied its subcellular localization by endogenous protein HA tagging using CRISPR-Cas9 gene editing. The enzyme-catalyzed proximity labeling technique, the TurboID system, was employed to identify the proteins in proximity to Trx4.

RESULTS

Trx4 was identified as a dense granule protein of T. gondii predominantly expressed in the parasitophorous vacuole (PV) and was partially co-localized with GRA1 and GRA5. Functional analysis showed that deletion of trx4 markedly influenced the parasite lytic cycle, resulting in impaired host cell invasion capacity in both RH and Pru strains. Mutation of Trx domains in Trx4 in RH strain revealed that two Trx domains were important for the parasite invasion. By utilizing the TurboID system to biotinylate proteins in proximity to Trx4, we identified a substantial number of proteins, some of which are novel, and others are previously characterized, predominantly distributed in the dense granules. In addition, we uncovered three novel proteins co-localized with Trx4. Intriguingly, deletion of trx4 did not affect the localization of these three proteins. Finally, a virulence assay demonstrated that knockout of trx4 resulted in a significant attenuation of virulence and a significant reduction in brain cyst loads in mice.

CONCLUSIONS

Trx4 plays an important role in T. gondii invasion and virulence in Type I RH strain and Type II Pru strain. Combining the TurboID system with CRISPR-Cas9 technique revealed many PV-localized proximity proteins associated with Trx4. These findings suggest a versatile role of Trx4 in mediating the processes that occur in this distinctive intracellular membrane-bound vacuolar compartment.

The Toxoplasma monocarboxylate transporters are involved in the metabolism within the apicoplast and are linked to parasite survival

The apicoplast is a four-membrane plastid found in the apicomplexans, which harbors biosynthesis and organelle housekeeping activities in the matrix. However, the mechanism driving the flux of metabolites, in and out, remains unknown. Here, we used TurboID and genome engineering to identify apicoplast transporters in Toxoplasma gondii. Among the many novel transporters, we show that one pair of apicomplexan monocarboxylate transporters (AMTs) appears to have evolved from a putative host cell that engulfed a red alga. Protein depletion showed that AMT1 and AMT2 are critical for parasite growth. Metabolite analyses supported the notion that AMT1 and AMT2 are associated with biosynthesis of isoprenoids and fatty acids. However, stronger phenotypic defects were observed for AMT2, including in the inability to establish T. gondii parasite virulence in mice. This study clarifies, significantly, the mystery of apicoplast transporter composition and reveals the importance of the pair of AMTs in maintaining the apicoplast activity in apicomplexans.

Research Note: Preliminary functional analysis of EGF-like domains of Eimeria tenella microneme protein 7 identified in sporozoites and merozoites

Eimeria tenella, an obligate intracellular apicomplexan parasite, is the major causative agent of chicken coccidiosis. Some epidermal growth factor (EGF)-like domain-containing proteins of other members of apicomplexan parasites have been reported to contribute to parasite survival. To date, however, EGF-like domain-containing proteins of E. tenella are not well studied. In this study, a gene fragment that encodes 4 EGF-like domains of E. tenella microneme protein 7 (EGF-EtMIC7) was amplified and expressed using an Escherichia coli expression system. Following generation of polyclonal antibodies that recognize recombinant EGF-EtMIC7 (rEGF-EtMIC7), the expression of EtMIC7 in sporozoites and merozoites was examined. Moreover, its roles in cellular regulation were investigated. The native EtMIC7 in E. tenella sporozoites and merozoites was detected by using Western blot and indirect immunofluorescence assays. rEGF-EtMIC7 could activate Akt, whereas blockade of EGF receptor (EGFR) failed to induce Akt phosphorylation. Compared with the control group, LMH cells treated with rEGF-EtMIC7 showed increased cell proliferation and expressed higher levels of B cell leukemia/lymphoma 2 (BCL-2). These findings contribute to the better understanding of parasite-host interactions at the molecular level during E. tenella infection.

Expression profiles of host miRNAs and circRNAs and ceRNA network during Toxoplasma gondii lytic cycle

Toxoplasma gondii is an opportunistic protozoan parasite that is highly prevalent in the human population and can lead to adverse health consequences in immunocompromised patients and pregnant women. Noncoding RNAs, such as microRNAs (miRNAs) and circular RNAs (circRNAs), play important regulatory roles in the pathogenesis of many infections. However, the differentially expressed (DE) miRNAs and circRNAs implicated in the host cell response during the lytic cycle of T. gondii are unknown. In this study, we profiled the expression of miRNAs and circRNAs in human foreskin fibroblasts (HFFs) at different time points after T. gondii infection using RNA sequencing (RNA-seq). We identified a total of 7, 7, 27, 45, 70, 148, 203, and 217 DEmiRNAs and 276, 355, 782, 1863, 1738, 6336, 1229, and 1680 DEcircRNAs at 1.5, 3, 6, 9, 12, 24, 36, and 48 h post infection (hpi), respectively. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses revealed that the DE transcripts were enriched in immune response, apoptosis, signal transduction, and metabolism-related pathways. These findings provide new insight into the involvement of miRNAs and circRNAs in the host response to T. gondii infection.

A newly characterized dense granule protein (GRA76) is important for the growth and virulence of Toxoplasma gondii

Pathogenicity of the zoonotic pathogen Toxoplasma gondii largely depends on the secretion of effector proteins into the extracellular milieu and host cell cytosol, including the dense granule proteins (GRAs). The protein-encoding gene TGME49_299780 was previously identified as a contributor to parasite fitness. However, its involvement in parasite growth, virulence and infectivity in vitro and in vivo remains unknown. Here, we comprehensively examined the role of this new protein, termed GRA76, in parasite pathogenicity. Subcellular localization revealed high expression of GRA76 in tachyzoites inside the parasitophorous vacuole (PV). However, its expression was significantly decreased in bradyzoites. A CRISPR-Cas9 approach was used to knock out the gra76 gene in the T. gondii type I RH strain and type II Pru strain. The in vitro plaque assays and intracellular replication showed the involvement of GRA76 in replication of RH and Pru strains. Deletion of the gra76 gene significantly decreased parasite virulence, and reduced the brain cyst burden in mice. Using RNA sequencing, we detected a significant increase in the expression of bradyzoite-associated genes such as BAG1 and LDH2 in the PruΔgra76 strain compared with the wild-type Pru strain. Using an in vitro bradyzoite differentiation assay, we showed that loss of GRA76 significantly increased the propensity for parasites to form bradyzoites. Immunization with PruΔgra76 conferred partial protection against acute and chronic infection in mice. These findings show the important role of GRA76 in the pathogenesis of T. gondii and highlight the potential of PruΔgra76 as a candidate for a live-attenuated vaccine.

Immunization with Live-Attenuated RHΔhad2a Strain Confers Partial Protective Immunity against Acute and Chronic Infection of Toxoplasma gondii in Mice

Toxoplasmosis caused by Toxoplasma gondii is an important zoonosis of human and animal health significance. Current chemical therapeutics have side effects, and no commercially available vaccine is licensed for the prevention of toxoplasmosis in humans and most animals. Developing a safe and effective vaccine with long-term protection against T. gondii infection is necessary to control toxoplasmosis. HAD2a is a key member of the haloacid dehalogenase (HAD) phosphatase family, which is essential for T. gondii daughter budding. However, the role of HAD2a in T. gondii virulence remains unknown. In this study, we successfully constructed the had2a gene knockout strain in the T. gondii-type I RH strain (RHΔhad2a) and determined its role in virulence and vaccination. These results demonstrate that HAD2a played an important role in parasite daughter budding and in vitro replication. Knockout of the had2a gene attenuated the virulence of the T. gondii-type I RH strain. Vaccination with RHΔhad2a tachyzoites induced a Th1-biased immune response, provided partial protection against acute T. gondii infection in mice by highly virulent tachyzoites of RH and PYS (ToxoDB#9, Chinese I) strains, and conferred strong protection against challenge infection by cysts and oocysts of the less virulent type II Pru strain. These results demonstrate that T. gondii had2a is important for its in vitro proliferation and virulence in mice and that RHΔhad2a may be used as a candidate strain to generate a multiple gene knockout live-attenuated strain or be collaboratively applied with other live-attenuated strains to confer more effective protection against T. gondii infection.

The transcription factor AP2XI-2 is a key negative regulator of Toxoplasma gondii merogony

Sexual development in Toxoplasma gondii is a multistep process that culminates in the production of oocysts, constituting approximately 50% of human infections. However, the molecular mechanisms governing sexual commitment in this parasite remain poorly understood. Here, we demonstrate that the transcription factors AP2XI-2 and AP2XII-1 act as negative regulators, suppressing merozoite-primed pre-sexual commitment during asexual development. Depletion of AP2XI-2 in type II Pru strain induces merogony and production of mature merozoites in an alkaline medium but not in a neutral medium. In contrast, AP2XII-1-depleted Pru strain undergoes several rounds of merogony and produces merozoites in a neutral medium, with more pronounced effects observed under alkaline conditions. Additionally, we identified two additional AP2XI-2-interacting proteins involved in repressing merozoite programming. These findings underscore the intricate regulation of pre-sexual commitment by a network of factors and suggest that AP2XI-2 or AP2XII-1-depleted Pru parasites can serve as a model for studying merogony in vitro.

Prevalence and genotype/subtype distribution of Enterocytozoon bieneusi and Blastocystis in donkeys in Shanxi Province, north China

Enterocytozoon bieneusi and Blastocystis may cause diarrhea in humans and various animals. However, little information is available regarding the prevalence and genetic diversity of E. bieneusi and Blastocystis in donkeys. To fill this gap, we molecularly assessed E. bieneusi and Blastocystis in fecal samples from donkeys (n = 815) in Shanxi Province, north China. The overall prevalence of E. bieneusi and Blastocystis in donkeys was 8.1% and 0.2%, respectively. Region and age were risk factors associated with E. bieneusi infection in donkeys. Three internal transcribed spacer (ITS) genotypes of E. bieneusi were identified in the current study, including two previously described genotypes (D and Henan-IV) and one novel genotype (named SXD1). Of which, genotype D was found to be the most prevalent. Phylogenetic analysis demonstrated that the three genotypes belonged to group 1, implying a potential of zoonotic transmission. Multilocus sequence typing showed that 19, 15, 13, and 22 types were identified at the loci MS1, MS3, MS4, and MS7, respectively, forming six multilocus genotypes (MLGs) distributed in the genotype D. One Blastocystis subtype (ST33) was identified, which has previously been reported only in horses. This is the first molecular-based description of E. bieneusi and Blastocystis infections in donkeys in Shanxi Province, north China, contributing to a better understanding of transmission dynamics and molecular epidemiological characteristics of the two intestinal protozoa.

Prevalence and Multilocus Genotyping of Giardia duodenalis in Donkeys in Shanxi Province, North China

Giardia duodenalis is a ubiquitous flagellated protozoan, causing significant economic losses to animal husbandry and posing threats to public health. China ranks the world's sixth largest major producer of donkeys, rearing approximately 2.6 million donkeys in 2019, but limited investigation of G. duodenalis prevalence has been conducted in the past, and it is yet to be known whether donkeys in Shanxi Province are infected with G. duodenalis. In the present study, a total of 815 fecal samples collected from donkeys in representative regions of Shanxi Province, North China, were examined for G. duodenalis using nested PCR. Then, the assemblages and multilocus genotypes (MLGs) were examined based on three established loci: namely, β-giardin (bg), triosephosphate isomerase (tpi), and glutamate dehydrogenase (gdh). The overall prevalence of G. duodenalis in donkeys in Shanxi Province was 16.81% (137/815). The region was identified as the main risk factor for the observed difference in G. duodenalis prevalence in donkeys among the three study areas (χ2 = 21.611, p < 0.001). Assemblages A, E, and B were identified, with the latter as the predominant assemblage. Three MLGs (MLG-novel-1 to 3) were formed based on sequence variation among the three loci. The present study reveals the presence of G. duodenalis in donkeys in Shanxi Province, North China, for the first time, which not only enriches the data on the distribution of G. duodenalis in donkeys in China but also provides useful baseline data for planning control strategies against G. duodenalis infection in the sampled areas.

Live attenuated anticoccidial vaccines for chickens

Chicken coccidiosis, caused by infection with single or multiple Eimeria species, results in significant economic losses to the global poultry industry. Over the past decades, considerable efforts have been made to generate attenuated Eimeria strains, and the use of live attenuated anticoccidial vaccines for disease prevention has achieved tremendous success. In this review, we evaluate the advantages and limitations of the methods of attenuation as well as attenuated Eimeria strains in a historical perspective. Also, we summarize the recent exciting research advances in transient/stable transfection systems and clustered regularly interspaced short palindromic repeats (CRISPR)-based genome editing developed for Eimeria parasites, and discuss trends and challenges of developing live attenuated anticoccidial vaccines based on transgenesis and genome editing.

The Toxoplasma protein phosphatase 6 catalytic subunit (TgPP6C) is essential for cell cycle progression and virulence

Protein phosphatases are post-translational regulators of Toxoplasma gondii proliferation, tachyzoite-bradyzoite differentiation and pathogenesis. Here, we identify the putative protein phosphatase 6 (TgPP6) subunits of T. gondii and elucidate their role in the parasite lytic cycle. The putative catalytic subunit TgPP6C and regulatory subunit TgPP6R likely form a complex whereas the predicted structural subunit TgPP6S, with low homology to the human PP6 structural subunit, does not coassemble with TgPP6C and TgPP6R. Functional studies showed that TgPP6C and TgPP6R are essential for parasite growth and replication. The ablation of TgPP6C significantly reduced the synchronous division of the parasite's daughter cells during endodyogeny, resulting in disordered rosettes. Moreover, the six conserved motifs of TgPP6C were required for efficient endodyogeny. Phosphoproteomic analysis revealed that ablation of TgPP6C predominately altered the phosphorylation status of proteins involved in the regulation of the parasite cell cycle. Deletion of TgPP6C significantly attenuated the parasite virulence in mice. Immunization of mice with TgPP6C-deficient type I RH strain induced protective immunity against challenge with a lethal dose of RH or PYS tachyzoites and Pru cysts. Taken together, the results show that TgPP6C contributes to the cell division, replication and pathogenicity in T. gondii.

Functional Characterization of Eight Zinc Finger Motif-Containing Proteins in Toxoplasma gondii Type I RH Strain Using the CRISPR-Cas9 System

The Zinc finger protein (ZFP) family is widely distributed in eukaryotes and interacts with DNA, RNA, and various proteins to participate in many molecular processes. In the present study, the biological functions of eight ZFP genes in the lytic cycle and the pathogenicity of Toxoplasma gondii were examined using the CRISPR-Cas9 system. Immunofluorescence showed that four ZFPs (RH248270-HA, RH255310-HA, RH309200-HA, and RH236640-HA) were localized in the cytoplasm, and one ZFP (RH273150-HA) was located in the nucleus, while the expression level of RH285190-HA, RH260870-HA, and RH248450-HA was undetectable. No significant differences were detected between seven RHΔzfp strains (RHΔ285190, RHΔ248270, RHΔ260870, RHΔ255310, RHΔ309200, RHΔ248450, and RHΔ236640) and the wild-type (WT) strain in the T. gondii lytic cycle, including plaque formation, invasion, intracellular replication, and egress, as well as in vitro virulence (p > 0.05). However, the RHΔ273150 strain exhibited significantly lower replication efficiency compared to the other seven RHΔzfp strains and the WT strain, while in vivo virulence in mice was not significantly affected. Comparative expression analysis of the eight zfp genes indicates that certain genes may have essential functions in the sexual reproductive stage of T. gondii. Taken together, these findings expand our current understanding of the roles of ZFPs in T. gondii.

Transcriptomic, 16S ribosomal ribonucleic acid and network pharmacology analyses shed light on the anticoccidial mechanism of green tea polyphenols against Eimeria tenella infection in Wuliangshan black-boned chickens

BACKGROUND

Eimeria tenella is an obligate intracellular parasitic protozoan that invades the chicken cecum and causes coccidiosis, which induces acute lesions and weight loss. Elucidating the anticoccidial mechanism of action of green tea polyphenols could aid the development of anticoccidial drugs and resolve the problem of drug resistance in E. tenella.

METHODS

We constructed a model of E. tenella infection in Wuliangshan black-boned chickens, an indigenous breed of Yunnan Province, China, to study the efficacy of green tea polyphenols against the infection. Alterations in gene expression and in the microbial flora in the cecum were analyzed by ribonucleic acid (RNA) sequencing and 16S ribosomal RNA (rRNA) sequencing. Quantitative real-time polymerase chain reaction was used to verify the host gene expression data obtained by RNA sequencing. Network pharmacology and molecular docking were used to clarify the interactions between the component green tea polyphenols and the targeted proteins; potential anticoccidial herbs were also analyzed.

RESULTS

Treatment with the green tea polyphenols led to a reduction in the lesion score and weight loss of the chickens induced by E. tenella infection. The expression of matrix metalloproteinase 7 (MMP7), MMP1, nitric oxide synthase 2 and ephrin type-A receptor 2 was significantly altered in the E. tenella infection plus green tea polyphenol-treated group and in the E. tenella infection group compared with the control group; these genes were also predicted targets of tea polyphenols. Furthermore, the tea polyphenol (-)-epigallocatechin gallate acted on most of the targets, and the molecular docking analysis showed that it has good affinity with interferon induced with helicase C domain 1 protein. 16S ribosomal RNA sequencing showed that the green tea polyphenols had a regulatory effect on changes in the fecal microbiota induced by E. tenella infection. In total, 171 herbs were predicted to act on two or three targets in MMP7, MMP1, nitric oxide synthase 2 and ephrin type-A receptor 2.

CONCLUSIONS

Green tea polyphenols can directly or indirectly regulate host gene expression and alter the growth of microbiota. The results presented here shed light on the mechanism of action of green tea polyphenols against E. tenella infection in chickens, and have implications for the development of novel anticoccidial products.

Molecular Identification and Genotyping of Cryptosporidium spp. and Blastocystis sp. in Cattle in Representative Areas of Shanxi Province, North China

Both Cryptosporidium spp. and Blastocystis sp. are common intestinal protozoa, which can cause zoonotic diseases and economic losses to livestock industry. To evaluate the prevalence and genetic population structure of Cryptosporidium spp. and Blastocystis sp. in beef and dairy cattle in Shanxi Province, north China, a total of 795 fecal samples were collected from beef and dairy cattle in three representative counties in Shanxi Province, and these fecal samples were examined using molecular approaches based on 18S small-subunit ribosomal RNA (SSU rRNA) of Cryptosporidium spp. and Blastocystis sp., respectively. Among 795 cattle fecal samples, 23 were detected as Cryptosporidium-positive and 103 were detected as Blastocystis-positive, and the overall prevalence of Cryptosporidium spp. and Blastocystis sp. in cattle in Shanxi Province was 2.9% and 13.0%, respectively. For Cryptosporidium spp., DNA sequence analysis indicated that all 23 positive samples were identified as C. andersoni. Furthermore, five known subtypes (ST1, ST10, ST14, ST21 and ST26) and three unknown subtypes of Blastocystis sp. were detected among 103 positive samples using DNA sequence analysis. This study reported the occurrence and prevalence of Cryptosporidium spp. and Blastocystis sp. in cattle in Shanxi Province for the first time, which extends the geographical distribution of these two zoonotic parasites and provides baseline data for the prevention and control of these two important zoonotic parasites in cattle in Shanxi Province.

Prevalence and Genetic Characterization of Blastocystis in Sheep and Pigs in Shanxi Province, North China: From a Public Health Perspective

Blastocystis is a common zoonotic intestinal protozoan and causes a series of gastrointestinal symptoms in humans and animals via the fecal-oral route, causing economic losses and posing public health problems. At present, the prevalence and genetic structure of Blastocystis in sheep and pigs in Shanxi province remains unknown. Thus, the present study collected 492 sheep fecal samples and 362 pig fecal samples from three representative counties in northern, central and southern Shanxi province for the detection of Blastocystis based on its SSU rRNA gene. The results showed that the overall prevalence of Blastocystis in the examined sheep and pigs were 16.26% and 14.09%, respectively. Sequences analyses showed that four known subtypes (ST5, ST10, ST14 and ST30) in sheep and two subtypes (ST1 and ST5) in pigs were detected in this study, with ST5 being the predominate subtype among the study areas. Phylogenetic analysis showed that the same subtypes were clustered into the same branch. This study reveals that sheep and pigs in Shanxi province are hosts for multiple Blastocystis subtypes, including the zoonotic subtypes (ST1 and ST5), posing a risk to public health. Baseline epidemiological data are provided that help in improving our understanding of the role of zoonotic subtypes in Blastocystis transmission.

Functional screening reveals Toxoplasma prenylated proteins required for endocytic trafficking and rhoptry protein sorting

In the apicomplexans, endocytosed cargos (e.g., hemoglobin) are trafficked to a specialized organelle for digestion. This follows a unique endocytotic process at the micropore/cytostome in these parasites. However, the mechanism underlying endocytic trafficking remains elusive, due to the repurposing of classical endocytic proteins for the biogenesis of apical organelles. To resolve this issue, we have exploited the genetic tractability of the model apicomplexan Toxoplasma gondii, which ingests host cytosolic materials (e.g., green fluorescent protein[GFP]). We determined an association between protein prenylation and endocytic trafficking, and using an alkyne-labeled click chemistry approach, the prenylated proteome was characterized. Genome editing, using clustered regularly interspaced short palindromic repaet/CRISPR-associated nuclease 9 (CRISPR/Cas9), was efficiently utilized to generate genetically modified lines for the functional screening of 23 prenylated candidates. This identified four of these proteins that regulate the trafficking of endocytosed GFP vesicles. Among these proteins, Rab1B and YKT6.1 are highly conserved but are non-classical endocytic proteins in eukaryotes. Confocal imaging analysis showed that Rab1B and Ras are substantially localized to both the trans-Golgi network and the endosome-like compartments in the parasite. Conditional knockdown of Rab1B caused a rapid defect in secretory trafficking to the rhoptry bulb, suggesting a trafficking intersection role for the key regulator Rab1B. Further experiments confirmed a critical role for protein prenylation in regulating the stability/activity of these proteins (i.e., Rab1B and YKT6.1) in the parasite. Our findings define the molecular basis of endocytic trafficking and reveal a potential intersection function of Rab1B on membrane trafficking in T. gondii. This might extend to other related protists, including the malarial parasites. IMPORTANCE The protozoan Toxoplasma gondii establishes a permissive niche, in host cells, that allows parasites to acquire large molecules such as proteins. Numerous studies have demonstrated that the parasite repurposes the classical endocytic components for secretory sorting to the apical organelles, leaving the question of endocytic transport to the lysosome-like compartment unclear. Recent studies indicated that endocytic trafficking is likely to associate with protein prenylation in malarial parasites. This information promoted us to examine this association in the model apicomplexan T. gondii and to identify the key components of the prenylated proteome that are involved. By exploiting the genetic tractability of T. gondii and a host GFP acquisition assay, we reveal four non-classical endocytic proteins that regulate the transport of endocytosed cargos (e.g., GFP) in T. gondii. Thus, we extend the principle that protein prenylation regulates endocytic trafficking and elucidate the process of non-classical endocytosis in T. gondii and potentially in other related protists.

Seroprevalence of Toxoplasma gondii infection in sheep and cattle in Shanxi Province, North China

Toxoplasmosis is a worldwide zoonotic disease caused by infection with the intracellular protozoan parasite Toxoplasma gondii, posing significant economic losses to the livestock industry. As a major livestock province, little is known of the prevalence of T. gondii infection in sheep and cattle in Shanxi Province, North China. In this study, a total of 1962 blood samples from cattle (n = 978) and sheep (n = 984), collected from 11 administrative cities in Shanxi Province, were examined for antibodies against T. gondii by using the indirect enzyme linked immunosorbent assay (ELISA) kits commercially available. The results showed that antibodies to T. gondii were detected in 306 of the 978 cattle serum samples (31.29%, 95% CI 28.38-34.19), ranging from 12.64% to 60.00% among the different cities. The overall seroprevalence of T. gondii in sheep was 17.78% (175/984, 95% CI 15.40-20.17), ranging from 2.22% to 41.11% among the different administrative cities. The T. gondii seroprevalence was associated with the management mode and geographical location. This is the first report of T. gondii seroprevalence in cattle and sheep in Shanxi Province, North China, which provides baseline data to plan future control strategies for T. gondii infection in this province.

Dynamic RNA profiles in the small intestinal epithelia of cats after Toxoplasma gondii infection

BACKGROUND

Felids are the only definitive hosts of Toxoplasma gondii. However, the biological features of the feline small intestine following T. gondii infection are poorly understood. We investigated the changes in the expression of RNAs (including mRNAs, long non-coding RNAs and circular RNAs) in the small intestinal epithelia of cats following T. gondii infection to improve our understanding of the life cycle of T. gondii and cat responses to T. gondii infection.

METHODS

Fifteen cats were randomly assigned to five groups, and the infection groups were inoculated with 600 tissue cysts of the T. gondii Pru strain by gavage. The small intestinal epithelia of cats were collected at 6, 10, 14, and 30 days post infection (DPI). Using high-throughput RNA sequencing (RNA-seq), we investigated the changes in RNA expression. The expression levels of differentially expressed (DE) genes and non-coding RNAs (ncRNAs) identified by RNA-seq were validated by quantitative reverse transcription PCR (qRT-PCR). Differential expression was determined using the DESeq R package.

RESULTS

In total, 207 annotated lncRNAs, 20,552 novel lncRNAs, 3342 novel circRNAs and 19,409 mRNAs were identified. Among these, 70 to 344 DE mRNAs, lncRNAs and circRNAs were detected, and the post-cleavage binding sites between 725 ncRNAs and 2082 miRNAs were predicted. Using the co-location method, we predicted that a total of 235 lncRNAs target 1044 protein-coding genes, while the results of co-expression analysis revealed that 174 lncRNAs target 2097 mRNAs. Pathway enrichment analyses of the genes targeted by ncRNAs suggested that most ncRNAs were significantly enriched in immune or diseases-related pathways. NcRNA regulatory networks revealed that a single ncRNA could be directly or indirectly regulated by multiple genes or ncRNAs that could influence the immune response of cats. Co-expression analysis showed that 242 circRNAs, mainly involved in immune responses, were significantly associated with T. gondii infection. In contrast, 1352 protein coding RNAs, mainly involved in nucleic acid process/repair pathways or oocyte development pathways, were negatively associated with T. gondii infection.

CONCLUSIONS

This study is the first to reveal the expression profiles of circRNAs, lncRNAs and mRNAs in the cat small intestine following T. gondii infection and will facilitate the elucidation of the role of ncRNAs in the pathogenesis of T. gondii infection in its definitive host, thereby facilitating the development of novel intervention strategies against T. gondii infection in humans and animals.

Interaction between tissue-dwelling helminth and the gut microbiota drives mucosal immunoregulation

Tissue-dwelling helminths affect billions of people around the world. They are potent manipulators of the host immune system, prominently by promoting regulatory T cells (Tregs) and are generally associated with a modified host gut microbiome. However, the role of the gut microbiota in the immunomodulatory processes for these non-intestinal parasites is still unclear. In the present study, we used an extra-intestinal cestode helminth model-larval Echinococcus multilocularis to explore the tripartite partnership (host-helminth-bacteria) in the context of regulating colonic Tregs in Balb/c mice. We showed that larval E. multilocularis infection in the peritoneal cavity attenuated colitis in Balb/c mice and induced a significant expansion of colonic Foxp3⁺ Treg populations. Fecal microbiota depletion and transplantation experiments showed that the gut microbiota contributed to increasing Tregs after the helminth infection. Shotgun metagenomic and metabolic analyses revealed that the gut microbiome structure after infection was significantly shifted with a remarkable increase of Lactobacillus reuteri and that the microbial metabolic capability was reprogrammed to produce more Treg cell regulator-short-chain fatty acids in feces. Furthermore, we also prove that the L. reuteri strain elevated in infected mice was sufficient to promote the colonic Treg frequency and its growth was potentially associated with T cell-dependent immunity in larval E. multilocularis infection. Collectively, these findings indicate that the extraintestinal helminth drives expansions of host colonic Tregs through the gut microbes. This study suggests that the gut microbiome serves as a critical component of anti-inflammation effects even for a therapy based on an extraintestinal helminth.

The antioxidant protein glutaredoxin 1 is essential for oxidative stress response and pathogenicity of Toxoplasma gondii

Glutaredoxins (Grxs) are ubiquitous antioxidant proteins involved in many molecular processes to protect cells against oxidative damage. Here, we study the roles of Grxs in the pathogenicity of Toxoplasma gondii. We show that Grxs are localized in the mitochondria (Grx1), cytoplasm (Grx2), and apicoplast (Grx3, Grx4), while Grx5 had an undetectable level of expression. We generated Δgrx1-5 mutants of T. gondii type I RH and type II Pru strains using CRISPR-Cas9 system. No significant differences in the infectivity were detected between four Δgrx (grx2-grx5) strains and their respective wild-type (WT) strains in vitro or in vivo. Additionally, no differences were detected in the production of reactive oxygen species, total antioxidant capacity, superoxide dismutase activity, and sensitivity to external oxidative stimuli. Interestingly, RHΔgrx1 or PruΔgrx1 exhibited significant differences in all the investigated aspects compared to the other grx2-grx5 mutant and WT strains. Transcriptome analysis suggests that deletion of grx1 altered the expression of genes involved in transport and metabolic pathways, signal transduction, translation, and obsolete oxidation-reduction process. The data support the conclusion that grx1 supports T. gondii resistance to oxidative killing and is essential for the parasite growth in cultured cells and pathogenicity in mice and that the active site CGFS motif was necessary for Grx1 activity.

Modulation of long noncoding RNA (lncRNA) and messenger RNA (mRNA) expression in the liver of Beagle dogs by Toxocara canis infection

BACKGROUND

Long non-coding RNAs (lncRNAs) and messenger RNAs (mRNAs) play crucial roles in regulating various physiological and pathological processes. However, the role of lncRNAs and mRNAs in mediating the liver response during Toxocara canis infection remains incompletely understood.

METHODS

In the present study, the expression profile of lncRNAs and mRNAs was investigated in the liver of Beagle dogs infected by T. canis using high-throughput RNA sequencing.

RESULTS

Compared with the control groups, 876 differentially expressed (DE) lncRNAs and 288 DEmRNAs were identified at 12 h post-infection (hpi), 906 DElncRNAs and 261 DEmRNAs were identified at 24 hpi, and 876 DElncRNAs and 302 DEmRNAs were identified at 36 days post-infection (dpi). A total of 16 DEmRNAs (e.g. dpp4, crp and gnas) were commonly identified at the three infection stages. Enrichment and co-localization analyses identified several pathways involved in immune and inflammatory responses during T. canis infection. Some novel DElncRNAs, such as LNC_015756, LNC_011050 and LNC_011052, were also associated with immune and inflammatory responses. Also, LNC_005105 and LNC_005401 were associated with the secretion of anti-inflammatory cytokines, which may play a role in the healing of liver pathology at the late stage of infection.

CONCLUSIONS

Our data provided new insight into the regulatory roles of lncRNAs and mRNAs in the pathogenesis of T. canis and improved our understanding of the contribution of lncRNAs and mRNAs to the immune and inflammatory response of the liver during T. canis infection.

The Toxoplasma micropore mediates endocytosis for selective nutrient salvage from host cell compartments

Apicomplexan parasite growth and replication relies on nutrient acquisition from host cells, in which intracellular multiplication occurs, yet the mechanisms that underlie the nutrient salvage remain elusive. Numerous ultrastructural studies have documented a plasma membrane invagination with a dense neck, termed the micropore, on the surface of intracellular parasites. However, the function of this structure remains unknown. Here we validate the micropore as an essential organelle for endocytosis of nutrients from the host cell cytosol and Golgi in the model apicomplexan Toxoplasma gondii. Detailed analyses demonstrated that Kelch13 is localized at the dense neck of the organelle and functions as a protein hub at the micropore for endocytic uptake. Intriguingly, maximal activity of the micropore requires the ceramide de novo synthesis pathway in the parasite. Thus, this study provides insights into the machinery underlying acquisition of host cell-derived nutrients by apicomplexan parasites that are otherwise sequestered from host cell compartments.

Genome-wide exploration reveals distinctive northern and southern variants of Clonorchis sinensis in the Far East

Clonorchis sinensis is a carcinogenic liver fluke that causes clonorchiasis-a neglected tropical disease (NTD) affecting ~35 million people worldwide. No vaccine is available, and chemotherapy relies on one anthelmintic, praziquantel. This parasite has a complex life history and is known to infect a range of species of intermediate (freshwater snails and fish) and definitive (piscivorous) hosts. Despite this biological complexity and the impact of this biocarcinogenic pathogen, there has been no previous study of molecular variation in this parasite on a genome-wide scale. Here, we conducted the first extensive nuclear genomic exploration of C. sinensis individuals (n = 152) representing five distinct populations from mainland China, and one from Far East Russia, and revealed marked genetic variation within this species between "northern" and "southern" geographical regions. The discovery of this variation indicates the existence of biologically distinct variants within C. sinensis, which may have distinct epidemiology, pathogenicity and/or chemotherapic responsiveness. The detection of high heterozygosity within C. sinensis specimens suggests that this parasite has developed mechanisms to readily adapt to changing environments and/or host species during its life history/evolution. From an applied perspective, the identification of invariable genes could assist in finding new intervention targets in this parasite, given the major clinical relevance of clonorchiasis. From a technical perspective, the genomic-informatic workflow established herein will be readily applicable to a wide range of other parasites that cause NTDs.

Effect of deleting four Toxoplasma gondii calcium-binding EGF domain-containing proteins on parasite replication and virulence

Several calcium-binding proteins including calcium-dependent protein kinases play important roles in several facets of the intracellular infection cycle of the apicomplexan protozoan parasite Toxoplasma gondii. However, the role of the calcium-binding epidermal growth factor (EGF) domain-containing proteins (CBDPs) remains poorly understood. In this study, we examined the functions of four CBDP genes in T. gondii RH strain of type I by generating knock-out strains using CRISPR-Cas9 system. We investigated the ability of mutant strains deficient in CBDP1, CBDP2, CBDP3, or CBDP4 to form plaques, replicate intracellularly, and egress from the host cells. The results showed that no definite differences between any of these four CBDP mutant strains and the wild-type strain in terms of their ability to form plaques, intracellular replication, and egress. Additionally, CBDP mutants did not exhibit any significant attenuated virulence compared to the wild-type strain in mice. The expression profiles of CBDP2-4 genes were conserved among T. gondii strains of different genotypes, life cycle stages, and developmental forms. Whether other CBDP genes play any roles in the pathogenicity of T. gondii strains of different genotypes remains to be elucidated.

Seroprevalence and Risk Factors of Chlamydia Infection in Cattle in Shanxi Province, North China

The information on Chlamydia infection in cattle is limited in Shanxi Province, north China. This study aimed to investigate the seroprevalence and risk factors of Chlamydia and Chlamydia abortus infection in cattle in Shanxi Province. In November 2020, a large-scale investigation of Chlamydia seroprevalence was conducted on 981 cattle serum samples collected from 40 cattle farms in 11 cities of Shanxi Province. The seroprevalence of Chlamydia and C. abortus was examined by indirect hemagglutination assay (IHA) and enzyme-linked immunosorbent assay (ELISA), respectively. The seroprevalence of Chlamydia and C. abortus was 52.29% (513/981) and 2.96% (29/981), respectively, in cattle in Shanxi Province. Location was identified as a risk factor for Chlamydia and C. abortus infection (p < 0.05). Under different management patterns, the seroprevalence of Chlamydia and C. abortus in large-scale animal farming companies was higher than that in household animal farms and animal farming cooperatives, and only the seroprevalence of Chlamydia was significantly different in different management patterns (p < 0.01). The results showed that there was higher seroprevalence of Chlamydia in cattle in Shanxi Province, while C. abortus was not the dominant species. This study provided baseline information on Chlamydia infection in cattle in Shanxi Province, which constitutes valuable data for monitoring livestock health and preventing potential zoonoses.

A combined miRNA-piRNA signature in the serum and urine of rabbits infected with Toxoplasma gondii oocysts

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.

Differential miRNA expression profiles in the bone marrow of Beagle dogs at different stages of Toxocara canis infection

BACKGROUND

Toxocara canis is distributed worldwide, posing a serious threat to both human and dog health; however, the pathogenesis of T. canis infection in dogs remains unclear. In this study, the changes in microRNA (miRNA) expression profiles in the bone marrow of Beagle dogs were investigated by RNA-seq and bioinformatics analysis.

RESULTS

Thirty-nine differentially expressed (DE) miRNAs (DEmiRNAs) were identified in this study. Among these, four DEmiRNAs were identified at 24 h post-infection (hpi) and all were up-regulated; eight DEmiRNAs were identified with two up-regulated miRNAs and six down-regulated miRNAs at 96 hpi; 27 DEmiRNAs were identified with 13 up-regulated miRNAs and 14 down-regulated miRNAs at 36 days post-infection (dpi). Among these DEmiRNAs, cfa-miR-193b participates in the immune response by regulating the target gene cd22 at 24 hpi. The novel_328 could participate in the inflammatory and immune responses through regulating the target genes tgfb1 and tespa1, enhancing the immune response of the host and inhibiting the infection of T. canis at 96 hpi. In addition, cfa-miR-331 and novel_129 were associated with immune response and self-protection mechanisms at 36 dpi. 20 pathways were significantly enriched by KEGG pathway analysis, most of which were related to inflammatory response, immune response and cell differentiation, such as Cell adhesion molecules (CAMs), ECM-receptor interaction and Focal adhesion.

CONCLUSIONS

These findings suggested that miRNAs of Beagle dog bone marrow play important roles in the pathogenesis of T. canis infection in dogs and provided useful resources to better understand the interaction between T. canis and the hosts.

The protein phosphatase 2A holoenzyme is a key regulator of starch metabolism and bradyzoite differentiation in Toxoplasma gondii

Phenotypic switching between tachyzoite and bradyzoite is the fundamental mechanism underpinning the pathogenicity and adaptability of the protozoan parasite Toxoplasma gondii. Although accumulation of cytoplasmic starch granules is a hallmark of the quiescent bradyzoite stage, the regulatory factors and mechanisms contributing to amylopectin storage in bradyzoites are incompletely known. Here, we show that T. gondii protein phosphatase 2A (PP2A) holoenzyme is composed of a catalytic subunit PP2A-C, a scaffold subunit PP2A-A and a regulatory subunit PP2A-B. Disruption of any of these subunits increased starch accumulation and blocked the tachyzoite-to-bradyzoite differentiation. PP2A contributes to the regulation of amylopectin metabolism via dephosphorylation of calcium-dependent protein kinase 2 at S679. Phosphoproteomics identified several putative PP2A holoenzyme substrates that are involved in bradyzoite differentiation. Our findings provide novel insight into the role of PP2A as a key regulator of starch metabolism and bradyzoite differentiation in T. gondii.

Seroprevalence of Neospora caninum infection and associated risk factors in cattle in Shanxi Province, north China

Neospora caninum is an obligate intracellular parasitic protozoan that can cause abortions in cattle and pose considerable economic losses to the cattle industry. As a major livestock province, little is known of N. caninum infection in cattle in Shanxi Province, north China. In order to investigate the seroprevalence of N. caninum in cattle in Shanxi Province, 978 cattle serum samples were collected from 11 cities in three representative geographical locations in Shanxi Province, and the N. caninum-specific IgG antibodies were examined using an indirect enzyme linked immunosorbent assay (ELISA) kit commercially available. The results showed that 133 of the 978 examined cattle serum samples (13.60%, 95% CI = 11.45-15.75) were positive for N. caninum antibodies, and the seroprevalence in different cities ranged from 0 to 78.89%. The geographical location and management mode were the risk factors associated with N. caninum infection in cattle herds in Shanxi Province. Cattle in Northern and Central Shanxi Province as well as cattle whose management mode is that of large-scale cattle farming companies are more susceptible to N. caninum infection. This was the first large-scale survey of N. caninum seroprevalence and assessment of associated risk factors in cattle in Shanxi Province, which provided baseline information for the prevention and control of N. caninum infection in cattle in Shanxi Province, north China.

The Seroprevalence of Chlamydia Infection in Sheep in Shanxi Province, China

Chlamydia, an obligate intracellular bacterium, can cause chlamydiosis in humans and animals worldwide and also leads to serious economic losses to the sheep industry. However, the information on Chlamydia infection in sheep was limited in Shanxi Province, northern China. In the present study, a total of 984 serum samples of sheep were collected from 11 regions in Shanxi Province, northern China in the autumn of 2020. The antibodies against Chlamydia and Chlamydia abortus were examined by the indirect hemagglutination assay (IHA) and indirect enzyme-linked immunosorbent assay (ELISA), respectively. The result showed that 351 (35.67%, 95% CI 32.68-38.66) of 984 serum samples were positive for Chlamydia, and the seroprevalence ranged from 6.67% to 70.79% among the different regions. In addition, antibodies to C. abortus infection were detected in 78 (7.93%, 95% CI 6.24-9.61) of 984 serum samples, and the seroprevalence ranged from 6.24% to 14.81% among the different regions. This is the first report on the seroprevalence of Chlamydia and C. abortus in sheep in Shanxi province, northern China. The findings provide baseline information for preventing and controlling Chlamydia infection in sheep in Shanxi Province, China.

Transcriptomic analysis of LMH cells in response to the overexpression of a protein of Eimeria tenella encoded by the locus ETH_00028350

A protein of Eimeria tenella (encoded by the locus ETH_00028350) homologous to Toxoplasma gondii dense granule protein 9, designated as EtHGRA9 hereafter, was reported to be expressed in all life cycle stages of E. tenella. However, no data are currently available regarding its functional properties. In the present study, a recombinant vector harboring a 741 bp gene segment encoding the mature form of EtHGRA9 was constructed and transfected into leghorn male hepatoma (LMH) cells. Then, transcriptomic analysis of the transfected LMH cells was carried out by using a high-throughput RNA-seq technology. The LMH cells overexpressing EtHGRA9 was validated by means of Western blotting as well as indirect immunofluorescence staining. The results demonstrated that the expression of 547 genes (275 upregulated genes and 272 downregulated genes) was altered by EtHGRA9. The quantitative real-time polymerase chain reaction (qRT-PCR) validation of the ten genes with differential expression between the two groups was consistent with the transcriptome analysis. According to pathway enrichment analysis for the obtained differentially expressed genes, seven pathways were significantly affected by EtHGRA9, such as cytokine-cytokine receptor interaction, MAPK signaling pathway, and protein processing in endoplasmic reticulum. Our data reveal several possible roles of EtHGRA9 in immune or inflammatory responses, which paves the way for a better understanding of the molecular interplay between E. tenella and its host.

H11-induced immunoprotection is predominantly linked to N-glycan moieties during Haemonchus contortus infection

Nematodes are one of the largest groups of animals on the planet. Many of them are major pathogens of humans, animals and plants, and cause destructive diseases and socioeconomic losses worldwide. Despite their adverse impacts on human health and agriculture, nematodes can be challenging to control, because anthelmintic treatments do not prevent re-infection, and excessive treatment has led to widespread drug resistance in nematode populations. Indeed, many nematode species of livestock animals have become resistant to almost all classes of anthelmintics used. Most efforts to develop commercial anti-nematode vaccines (native or recombinant) for use in animals and humans have not succeeded, although one effective (dead) vaccine (Barbervax) has been developed to protect animals against one of the most pathogenic parasites of livestock animals – Haemonchus contortus (the barber’s pole worm). This vaccine contains native molecules, called H11 and H-Gal-GP, derived from the intestine of this blood-feeding worm. In its native form, H11 alone consistently induces high levels (75-95%) of immunoprotection in animals against disease (haemonchosis), but recombinant forms thereof do not. Here, to test the hypothesis that post-translational modification (glycosylation) of H11 plays a crucial role in achieving such high immunoprotection, we explored the N-glycoproteome and N-glycome of H11 using the high-resolution mass spectrometry and assessed the roles of N-glycosylation in protective immunity against H. contortus. Our results showed conclusively that N-glycan moieties on H11 are the dominant immunogens, which induce high IgG serum antibody levels in immunised animals, and that anti-H11 IgG antibodies can confer specific, passive immunity in naïve animals. This work provides the first detailed account of the relevance and role of protein glycosylation in protective immunity against a parasitic nematode, with important implications for the design of vaccines against metazoan parasites.

Seroprevalence and risk factors of Toxoplasma gondii and Neospora caninum infection in black goats in Yunnan Province, Southwestern China

Toxoplasma gondii and Neospora caninum are two obligate intracellular protozoan parasites that can cause reproductive failure and production losses. To date, there is no data of T. gondii and N. caninum seroprevalence in black goats in Yunnan Province, southwestern China. In the present study, a total of 734 serum samples were collected from black goats in four different counties of Yunnan Province. 734 and 590 serum samples were examined for antibodies against T. gondii and N. caninum by using MAT and indirect ELISA, respectively. A total of 123 and 76 samples were T. gondii-positive and N. caninum-positive, respectively. The overall seroprevalence of T. gondii in black goats was 16.76% (123/734, 95% CI: 14.06-19.46) with the titer ranged from 1:25 to 1:3200. The seroprevalence of N. caninum was 12.88% (76/590, 95% CI: 10.18-15.58). There was significant difference in seroprevalence of N. caninum in different regions (P < 0.01, χ² = 30.63) and age groups (P < 0.05, χ² = 11.85). Significant differences in seroprevalence of T. gondii were observed in different regions (P < 0.05, χ² = 9.21) and different gender groups (P < 0.01, χ² = 12.29). Results of seroprevalence of T. gondii and N. caninum indicated that T. gondii and N. caninum were prevalent parasites in black goats in Yunnan Province. This is the first report of seroprevalence of T. gondii and N. caninum in black goats in Yunnan Province. The results of this study indicated that some measures should be taken to control these two parasites and to reduce economic losses to the livestock industry in Yunnan Province.

Global proteomic profiling of multiple organs of cats (Felis catus) and proteome-transcriptome correlation during acute Toxoplasma gondii infection

Background

Toxoplasma gondii is a protozoan parasite which can infect almost all warm-blooded animals and humans. Understanding the differential expression of proteins and transcripts associated with T. gondii infection in its definitive host (cat) may improve our knowledge of how the parasite manipulates the molecular microenvironment of its definitive host. The aim of this study was to explore the global proteomic alterations in the major organs of cats during acute T. gondii infection.

Methods

iTRAQ-based quantitative proteomic profiling was performed on six organs (brain, liver, lung, spleen, heart and small intestine) of cats on day 7 post-infection by cysts of T. gondii PRU strain (Genotype II). Mascot software was used to conduct the student’s t-test. Proteins with P values < 0.05 and fold change > 1.2 or < 0.83 were considered as differentially expressed proteins (DEPs).

Results

A total of 32,657 proteins were identified in the six organs, including 2556 DEPs; of which 1325 were up-regulated and 1231 were down-regulated. The brain, liver, lung, spleen, heart and small intestine exhibited 125 DEPs, 463 DEPs, 255 DEPs, 283 DEPs, 855 DEPs and 575 DEPs, respectively. Gene Ontology (GO) annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses of all proteins and DEPs in all organs showed that many proteins were enriched in binding, cell part, cell growth and death, signal transduction, translation, sorting and degradation, extracellular matrix remodeling, tryptophan catabolism, and immune system. Correlations between differentially expressed proteins and transcripts were detected in the liver (n = 19), small intestine (n = 17), heart (n = 9), lung (n = 9) and spleen (n = 3).

Conclusions

The present study identified 2556 DEPs in six cat tissues on day 7 after infection by T. gondii PRU strain, and functional enrichment analyses showed that these DEPs were associated with various cellular and metabolic processes. These findings provide a solid base for further in-depth investigation of the complex proteotranscriptomic reprogramming that mediates the dynamic interplays between T. gondii and the different feline tissues.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40249-022-01022-7.

Lipidomic changes in the liver of beagle dogs associated with Toxocara canis infection

A global lipidomic analysis using liquid chromatography–tandem mass spectrometry was performed on the liver of beagle dogs infected with Toxocara canis to profile hepatic lipid species at 12 h post-infection (hpi), 24 hpi, and 36 days post-infection (dpi). This analysis identified six categories and 42 subclasses of lipids, including 173, 64, and 116 differentially abundant lipid species at 12 hpi, 24 hpi, and 36 dpi, respectively. Many of the identified lysophospholipids, such as lysophosphatidylglycerol, lysophosphatidylserine, and lysophosphatidylcholine, may contribute to the migration and development of T. canis during the early infection stage. Pathway analysis revealed significant alterations of several immune-inflammatory pathways, such as the B-cell receptor signaling pathway, the NF-kappa B signaling pathway, and the C-type lectin receptor signaling pathway at 12 and 24 hpi. These findings demonstrate the value of lipidomic profiling in revealing the extent of changes in the composition and abundance of hepatic lipidome caused by T. canis infection and their relevance to the pathophysiology of toxocariasis in beagle dogs.

RNA sequencing reveals dynamic expression of lncRNAs and mRNAs in caprine endometrial epithelial cells induced by Neospora caninum infection

Background

The effective transmission mode of Neospora caninum, with infection leading to reproductive failure in ruminants, is vertical transmission. The uterus is an important reproductive organ that forms the maternal–fetal interface. Neospora caninum can successfully invade and proliferate in the uterus, but the molecular mechanisms underlying epithelial-pathogen interactions remain unclear. Accumulating evidence suggests that host long noncoding RNAs (lncRNAs) play important roles in cellular molecular regulatory networks, with reports that these RNA molecules are closely related to the pathogenesis of apicomplexan parasites. However, the expression profiles of host lncRNAs during N. caninum infection has not been reported.

Methods

RNA sequencing (RNA-seq) analysis was used to investigate the expression profiles of messenger RNAs (mRNAs) and lncRNAs in caprine endometrial epithelial cells (EECs) infected with N. caninum for 24 h (TZ_24h) and 48 h (TZ_48 h), and the potential functions of differentially expressed (DE) lncRNAs were predicted by using Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis of their mRNA targets.

Results

RNA-seq analysis identified 1280.15 M clean reads in 12 RNA samples, including six samples infected with N. caninum for 24 h (TZ1_24h-TZ3_24h) and 48 h (TZ1_48h-TZ3_48h), and six corresponding control samples (C1_24h-C3_24h and C1_48h-C3_48h). Within the categories TZ_24h-vs-C_24h, TZ_48h-vs-C_48h and TZ_48h-vs-TZ_24h, there were 934 (665 upregulated and 269 downregulated), 1238 (785 upregulated and 453 downregulated) and 489 (252 upregulated and 237 downregulated) DEmRNAs, respectively. GO enrichment and KEGG analysis revealed that these DEmRNAs were mainly involved in the regulation of host immune response (e.g. TNF signaling pathway, MAPK signaling pathway, transforming growth factor beta signaling pathway, AMPK signaling pathway, Toll-like receptor signaling pathway, NOD-like receptor signaling pathway), signaling molecules and interaction (e.g. cytokine-cytokine receptor interaction, cell adhesion molecules and ECM-receptor interaction). A total of 88 (59 upregulated and 29 downregulated), 129 (80 upregulated and 49 downregulated) and 32 (20 upregulated and 12 downregulated) DElncRNAs were found within the categories TZ_24h-vs-C_24h, TZ_48h-vs-C_48h and TZ_48h-vs-TZ_24h, respectively. Functional prediction indicated that these DElncRNAs would be involved in signal transduction (e.g. MAPK signaling pathway, PPAR signaling pathway, ErbB signaling pathway, calcium signaling pathway), neural transmission (e.g. GABAergic synapse, serotonergic synapse, cholinergic synapse), metabolism processes (e.g. glycosphingolipid biosynthesis-lacto and neolacto series, glycosaminoglycan biosynthesis-heparan sulfate/heparin) and signaling molecules and interaction (e.g. cytokine-cytokine receptor interaction, cell adhesion molecules and ECM-receptor interaction).

Conclusions

This is the first investigation of global gene expression profiles of lncRNAs during N. caninum infection. The results provide valuable information for further studies of the roles of lncRNAs during N. caninum infection.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-022-05405-5.

Transcriptional changes in LMH cells induced by Eimeria tenella rhoptry kinase family protein 17

Though a number of Eimeria tenella rhoptry kinase family proteins have been identified, little is known about their molecular functions. In the present study, the gene fragment encoding the matured peptide of E. tenella rhoptry kinase family protein 17 (EtROP17) was used to construct a recombinant vector, followed by transfection into leghorn male hepatoma (LMH) cells. Then, the transcriptional changes in the transfected cells were determined by RNA-seq. The expression of EtROP17 in LMH cells was validated by both Western blot and indirect immunofluorescence analysis. Our analysis showed that EtROP17 altered the expression of 309 genes (114 downregulated genes and 195 upregulated genes) in LMH cells. The quantitative real-time polymerase chain reaction (qRT-PCR) results of the selected differentially expressed genes (DEGs) were consistent with the RNA-seq data. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that DEGs were significantly enriched in nine pathways, such as toll-like receptor signaling pathway, ECM-receptor interaction, intestinal immune network for IgA production and focal adhesion. These findings reveal several potential roles of EtROP17, which contribute to understanding the molecular mechanisms underlying the host-parasite interplay.

First report of haemosporidia and associated risk factors in red junglefowl (Gallus gallus) in China

Background

Avian haemosporidia infect both domestic and wild birds, causing anemia, acute tissue degeneration, and depopulation in wild birds. Poultry and wild birds have been reported as common reservoirs of haemosporidia, but limited information is available for red junglefowl (Gallus gallus) in China. The present study investigated the prevalence and molecular characterization of haemosporidia in red junglefowl.

Methods

Blood samples were collected from 234 red junglefowl from Jinghong City of Yunnan Province, and genomic DNA was extracted from these samples. The prevalence of haemosporidia was determined by nested PCR targeting the mitochondrial cytochrome b (cytb) gene. Molecular characterization was investigated based on phylogenetic analysis of cytb sequences, and associated risk factors were analyzed using the Chi-square (χ²) test.

Results

The overall prevalence of haemosporidia was 74.8% (175/234), and three species were identified, namely Haemoproteus enucleator, Leucocytozoon californicus, and Plasmodium juxtanucleare. The prevalence of haemosporidia in adult fowl (81.1%, 107/132) was significantly higher (χ² = 6.32, df = 1, P = 0.012) than that in juveniles (66.7%, 68/102). Three novel haemosporidian lineages were revealed.

Conclusions

This study examined the prevalence and identified species of avian haemosporidians in red junglefowl, providing new information on the molecular epidemiology and geographical distribution of haemosporidian parasites. Our results indicated high prevalence and diverse species distribution of these haemosporidians in red junglefowl. To the best of our knowledge, this is the first record of haemosporidian infection in red junglefowl in China.

Graphical Abstract

Metabolomics study of cat small intestine during the early stage of Toxoplasma gondii oocyst formation identifies potential biomarkers

Toxoplasma gondii (T. gondii) is a zoonotic intracellular protozoan parasite that can invade, replicate and survive in almost all cells of warm-blooded animals. T. gondii infection threatens the life of the fetus or can cause morbidity in the infant. As the only definitive host of T. gondii, felids spread the pathogen mainly by forming oocysts in the small intestines and discharging the oocysts into the ambient environment, consequently polluting water, vegetables, and meat products. In this study, we used untargeted metabolomics technology to study the changes in metabolites that occurred during the early stage of oocyst formation in the cat small intestine following T. gondii infection and attempted to identify metabolic biomarkers that could potentially be used as diagnostic molecular markers in the future. Domestic cats (Felis catus) were infected with T. gondii Pru tissue cysts, and samples of their small intestinal epithelium were collected at 2 and 4 days post-infection (DPI) for metabolic analysis. LC-MS/MS and multivariate statistical analysis were employed to detect metabolomic signatures that discriminated between the infected and control groups. A total of 1673 ions and 1201 ions were obtained in the positive and negative modes, respectively. Of these ions, 175 were up-regulated and 127 were down-regulated in the positive ion mode; whereas, 123 were up-regulated and 81 were down-regulated in the negative ion mode. Three commonly altered ions (0.74_313.0414 m/z, 8.82_615.2621 m/z and 8.16_325.2362 m/z) were determined to have potential research value. Seventy common metabolic pathways were enriched at two time points, with arginine biosynthesis, pyrimidine metabolism, pantothenate and CoA biosynthesis being the three most significant pathways related to T. gondii. The area under the curve (AUC) of differential metabolites combined with relevant literature analysis showed that N-Methylpelletierine and 3,3-Difluoro-17-methyl-5alpha-androstan-17beta-ol have higher predictability and better potential application value than other metabolites. Our analysis of metabolic markers during the early stage of T. gondii oocyst formation in the small intestine of the definitive host (cat) provided novel insight for understanding oocyst development and a theoretical basis for the application of potential biomarkers.

Molecular Detection and Genotyping of Enterocytozoon bieneusi in Pigs in Shanxi Province, North China

Enterocytozoon bieneusi is a common opportunistic intestinal pathogen that can cause acute diarrhea in immunosuppressed humans and animals. Though E. bieneusi has been widely detected in pigs around the world, little is known of its prevalence and genotype distribution in pigs in Shanxi province, north China. In this study, a total of 362 fecal samples were collected from pigs in three representative counties in north, south, and central Shanxi province, China. The prevalence and genotypes of E. bieneusi were investigated by nested PCR amplification of the ribosomal internal transcribed spacer (ITS) region of the ribosomal RNA (rRNA) gene. Overall, the prevalence of E. bieneusi in pigs in Shanxi province was 54.70% (198/362). Statistical analysis showed the difference in prevalence was statistically significant between regions (χ2 = 41.94, df = 2, P &lt; 0.001) and ages (χ2 = 80.37, df = 1, P &lt; 0.001). In addition, 16 genotypes of E. bieneusi were identified in this study by sequence analysis of the ITS region, including 15 known genotypes (EbpC, EbpA, EbpB, pigEb4, PigEBITS5, I, Henan-I, G, WildBoar 7, SH10, EbpD, CHC5, PigSpEb1, PigSpEb2, and CHG19) and one novel genotype (designated as PigSX-1). Phylogenetic analysis revealed that 14 known genotypes and the novel genotype were clustered into Group 1, whereas genotype I belonged to Group 2. To the best of our knowledge, this is the first report on the prevalence and genotypes of E. bieneusi in pigs in Shanxi province. These findings enrich the genetic diversity of E. bieneusi and provide the baseline data for the prevention and control of E. bieneusi in pigs in the study regions.

Population genomics of ancient and modern Trichuris trichiura

The neglected tropical disease trichuriasis is caused by the whipworm Trichuris trichiura, a soil-transmitted helminth that has infected humans for millennia. Today, T. trichiura infects as many as 500 million people, predominantly in communities with poor sanitary infrastructure enabling sustained faecal-oral transmission. Using whole-genome sequencing of geographically distributed worms collected from human and other primate hosts, together with ancient samples preserved in archaeologically-defined latrines and deposits dated up to one thousand years old, we present the first population genomics study of T. trichiura. We describe the continent-scale genetic structure between whipworms infecting humans and baboons relative to those infecting other primates. Admixture and population demographic analyses support a stepwise distribution of genetic variation that is highest in Uganda, consistent with an African origin and subsequent translocation with human migration. Finally, genome-wide analyses between human samples and between human and non-human primate samples reveal local regions of genetic differentiation between geographically distinct populations. These data provide insight into zoonotic reservoirs of human-infective T. trichiura and will support future efforts toward the implementation of genomic epidemiology of this globally important helminth.

The whipworm Trichuris trichiura is a soil-transmitted helminth that causes the neglected tropical disease trichuriasis in humans. Here, the authors produce whole genome sequences of modern and ancient samples from humans and non-human primates to characterise the genomic diversity and evolution of this pathogen.

[Predictive value of the product of plasma colchicine concentration and poisoning time for the prognosis of colchicine poisoning patients]

Objective: To investigate the predictive value of the product of first plasmacolchicine concentration and poisoning time for the prognosis of colchicine poisoning patients, and to provide a basis for early prognosis assessment. Methods: October 2021, patients with colchicine poisoning admitted in the First Affiliated Hospitol of Wenzhou Medical University from January 2017 to September 2021 were collected, including general information such as patient gender, age, oral colchicine dose, poisoning time, the first laboratory test index andplasma colchicine concentration after admission. The patients were divided into survival group and death group according to their prognosis. The differences in clinical indicators such as admission plasma colchicine concentration, blood routine, blood biochemistry, coagulation function, and blood gas analysis were compared between the two groups, and their predictive value for the prognosis of patients were analyzed. Results: A total of 23 patients with colchicine poisoning, aged 20-85 years, were included in this study, of which 15 cases (65.22%) survived and 8 cases (34.78%) died. The first plasma colchicine concentration at admision were 0.42-53.61 ng/ml. The plasma colchicine concentration and the concentration-time product were 10.08-2147.04 h·ng/ml.Compared with the survival group, the plasma colchicine concentration and the concentration-time product in the death group were significantly increased, and the differences were statistically significant (P<0.05). Univariate logistic regression analysis showed that first plasma concentration and poisoning time>132.48 h·ng/ml, high C-reactive protein, high D-dimer, high absolute value of BE were the risk factors for the prognosis of patients with colchicine poisoning (OR=12.000, 95%CI: 1.1181-128.836; OR=1.053, 95%CI: 1.009-1.098; OR=1.219, 95%CI: 1.039-1.429; OR=1.360, 95%CI: 1.1.044-1.773; P<0.05). High prothrombin time activity was protective factor affecting the prognosis of colchicine poisoning patients (OR=0.941, 95%CI: 0.892~0.993; P<0.05). ROC curve analysis showed that the areas under the curves of first plasma concentration and poisoning time, C-reactive protein, absolute value of BE, D-dimer for predicting the prognosis of patients with colchicine poisoning were 0.918, 0.888, 0.867, 0.837, respectively, and the areas under the curves of prothrombin time activityfor predicting the prognosis of patients with colchicine poisoning was 0.788 (P<0.05) . Conclusion: The product of the first plasma colchicine concentration at admission and poisoning time is closely related to the prognosis of patients with colchicine poisoning, it can be used as a predictor for early evaluation of the prognosis of poisoned patients.

Human pediculosis, a global public health problem

Background

Human pediculosis is caused by hematophagous lice, which are transmitted between individuals via direct and/or indirect contact. Despite the public health importance of louse infestation, information concerning the global burden of pediculosis and the epidemiological landscape of louse-borne diseases is limited. The aim of this review was to summarize the biology, epidemiology, diagnosis, and control of lice infestation in humans. We also discussed the latest advances in molecular taxonomy and molecular genetics of lice.

Methods

We searched five electronic bibliographic databases (PubMed, ScienceDirect, CNKI, VIP Chinese Journal Database, and Wanfang Data) and followed a standard approach for conducting scoping reviews to identify studies on various aspects of human lice. Relevant information reported in the identified studies were collated, categorized, and summarized.

Results

A total of 282 studies were eligible for the final review. Human pediculosis remains a public health issue affecting millions of people worldwide. Emerging evidence suggests that head lice and body lice should be considered conspecific, with different genotypes and ecotypes. Phylogenetic analysis based on mitochondrial (mt) cytb gene sequences identified six distinct clades of lice worldwide. In addition to the direct effect on human health, lice can serve as vectors of disease-causing pathogens. The use of insecticides plays a crucial role in the treatment and prevention of louse infestation. Genome sequencing has advanced our knowledge of the genetic structure and evolutionary biology of human lice.

Conclusions

Human pediculosis is a public health problem affecting millions of people worldwide, particularly in developing countries. More progress can be made if emphasis is placed on the use of emerging omics technologies to elucidate the mechanisms that underpin the physiological, ecological, and evolutionary aspects of lice.

Graphic Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s40249-022-00986-w.

Echinococcosis Is Associated with the Increased Prevalence of Intestinal Blastocystis Infection in Tibetans and Host Susceptibility to the Blastocystis in Mice

Blastocystis is a common human intestinal protozoan parasite. Little is known about its prevalence in echinococcosis. This study tested whether Echinococcus multilocularis infection would increase host susceptibility to Blastocystis. A total of 114 fecal samples (68 hydatid disease patients and 46 healthy people) were collected from Tibetans in the Qinghai province in China. The presence of Blastocystis was identified by sequencing of the small subunit (SSU) rRNA gene. Balb/c mice were co-infected with Blastocystis and E. multilocularis and tested for host susceptibility to Blastocystis. The overall Blastocystis prevalence was 12.3%; 16.2% in the patients and 4.4% in healthy people (p < 0.05). Sequence analysis identified three known Blastocystis genotypes, including ST1, ST2, and ST3, and one unknown genotype. Experimental dual infection significantly reduced mouse survival rate (20%), induced more severe signs, and increased intestinal damages with a higher intestinal colonization level of Blastocystis. The mouse model showed that E. multilocularis infection increases host susceptibility to Blastocystis. Our study shows a significantly higher prevalence of Blastocystis in patients with liver echinococcosis and reveals that non-intestinal E. multilocularis infection increases host susceptibility to the Blastocystis. Our results highlight that E. multilocularis infection is associated with Blastocystis. These findings remind us that more attention should be paid to the gut health of the patients with a helminth infection during clinical patient care.

Global profiling of protein lysine malonylation in Toxoplasma gondii strains of different virulence and genetic backgrounds

Lysine malonylation is a post-translational modification (PTM), which regulates many cellular processes. Limited information is available about the level of lysine malonylation variations between Toxoplasma gondii strains of distinct genetic lineages. Yet, insights into such variations are needed to understand the extent to which lysine malonylation contributes to the differences in the virulence and repertoire of virulence factors between T. gondii genotypes. In this study, we profiled lysine malonylation in T. gondii using quantitative liquid chromatography-tandem mass spectrometry (LC-MS/MS) and immuno-affinity purification. This analysis was performed on three T. gondii strains with distinctive pathogenicity in mice, including RH strain (type I), PRU strain (type II), and VEG strain (type III). In total, 111 differentially malonylated proteins and 152 sites were upregulated, and 17 proteins and 17 sites were downregulated in RH strain versus PRU strain; 50 proteins and 59 sites were upregulated, 50 proteins and 53 sites were downregulated in RH strain versus VEG strain; and 72 proteins and 90 sites were upregulated, and 7 proteins and 8 sites were downregulated in VEG strain versus PRU strain. Differentially malonylated proteins were involved in key processes, such as those mediating the regulation of protein metabolism, stress response, glycolysis, and actin cytoskeleton. These results reveal an association between lysine malonylation and intra-species virulence differences in T. gondii and offer a new resource for elucidating the contribution of lysine malonylation to energy metabolism and virulence in T. gondii.

Lysine malonylation has been shown to play important roles in various biological processes in Toxoplasma gondii. Here, we used quantitative liquid chromatography-tandem mass spectrometry (LC-MS/MS) and immuno-affinity purification to test the hypothesis that lysine malonylation underpins the inter-genotype differences in the virulence of T. gondii. Several up-regulated and down- regulated malonylated proteins were identified in the tachyzoites of RH (type I) strain, PRU (type II) strain, and VEG (type III) strain. Differentially regulated malonylated proteins were enriched in many biological and metabolic pathways, and were found to contribute T. gondii energy metabolism, stress response, and infectivity, suggesting the role of lysine malonylation in the regulation of T. gondii virulence. These findings expand our knowledge of lysine malonylation in T. gondii and provide more insight into the mechanisms mediating the virulence differences between T. gondii strains of different genotypes.

Transcriptomic landscape of hepatic lymph nodes, peripheral blood lymphocytes and spleen of swamp buffaloes infected with the tropical liver fluke Fasciola gigantica

The tropical liver fluke Fasciola gigantica is a parasitic helminth that has been frequently reported to infect mammals, typically involving water buffaloes. In this study, we characterized the tissue transcriptional landscape of buffaloes following infection by F. gigantica. RNAs were isolated from hepatic lymph nodes (hLNs), peripheral blood lymphocytes (pBLs), and spleen at 3-, 42- and 70-days post-infection (dpi), and all samples were subjected to RNA sequencing analyses. At 3 dpi, 2603, 460, and 162 differentially expressed transcripts (DETs) were detected in hLNs, pBLs, and spleen, respectively. At 42 dpi, 322, 937, and 196 DETs were detected in hLNs, pBLs, and spleen, respectively. At 70 dpi, 376, 334, and 165 DETs were detected in hLNs, pBLs, and spleen, respectively. Functional enrichment analysis identified upregulated immune-related pathways in the infected tissues involved in innate and adaptive immune responses, especially in hLNs at 42 and 70 dpi, and pBLs at 3 and 42 dpi. The upregulated transcripts in spleen were not enriched in any immune-related pathway. Co-expression network analysis further identified transcriptional changes associated with immune response to F. gigantica infection. Receiver operating characteristic (ROC) curve analysis showed that 107 genes in hLNs, 32 genes in pBLs, and 36 genes in spleen correlated with F. gigantica load. These findings provide new insight into molecular mechanisms and signaling pathways associated with F. gigantica infection in buffaloes.

Fasciola gigantica is a socioeconomically important tropical liver fluke of mammals, causing fascioliasis–a neglected tropical disease. In the present study, RNA sequencing and bioinformatic approach were employed to explore the global transcriptional changes of hepatic lymph nodes (hLNs), peripheral blood lymphocytes (pBLs), and spleen of water buffaloes during F. gigantica infection at 3-, 42-, and 70-days post-infection (dpi). The results revealed significant transcriptional upregulation of genes associated with innate and adaptive immune responses in infected hLNs (42 and 70 dpi) and pBLs (3 and 42 dpi). However, downregulation of transcripts involved in immune response was detected in pBLs at 70 dpi. The downregulated transcripts were enriched in metabolic pathways, such as drug metabolism-cytochrome P450 in infected hLNs at 3 dpi. These findings provide new insight into the pathogenesis of F. gigantica in its natural mammalian host.

The Detection of Toxoplasma gondii in Wild Rats (Rattus norvegicus) on Mink Farms in Shandong Province, Eastern China

Toxoplasma gondii is a worldwide distributed zoonotic pathogen that threatens public health. However, there have been limited data for T. gondii infection in wild rats (Rattus norvegicus) in China. In the present study, a total of 227 wild rats were captured from three mink farms to investigate T. gondii infection in Shandong Province, eastern China. The DNA was extracted from 25 mg rats' brain tissues and subjected to a PCR amplification by targeting to the T. gondii B1. In 227 wild rat samples, 18 samples (7.93%) were positive for T. gondii. Then, the positive samples were further genotyped based on eight genetic markers, including eight nuclear loci (SAG1, 5'-SAG2, and 3'-SAG2, alternative SAG2, SAG3, GRA6, c29-2, and L358) and an apicoplast locus (Apico) by using the multilocus PCR-restriction fragment length polymorphism technology. Of these samples, eight were genotyped at nine nuclear loci, and two were genotyped at eight nuclear loci, forming three known genotypes (ToxoDB no. 43, ToxoDB no. 91, and ToxoDB no. 189) and two new genotypes. The closest ToxoDB genotypes were observed in wild rats, suggesting the differences in the population structure of the T. gondii between breed farm animals and wild rats. These data revealed the genetic variability of T. gondii in wild rats on mink farms in Shandong Province, with possible implication for public health.

Prevalence and genotypes/subtypes of Enterocytozoon bieneusi and Blastocystis sp. in different breeds of cattle in Jiangxi Province, southeastern China

Enterocytozoon bieneusi and Blastocystis sp. are common zoonotic pathogens that parasitize in the small intestine of humans and animals, posing a threat to public health. However, little information is available on the prevalence and genotypes/subtypes of E. bieneusi and Blastocystis sp. in cattle in Jiangxi Province, southeastern China. In the present study, 556 fecal samples of cattle were collected from Nanchang city, Gao'an city, Xinyu city, and Ji'an city in Jiangxi Province. All samples were examined for the presence of E. bieneusi by nested PCR analysis of the ribosomal internal transcribed spacer (ITS) and Blastocystis sp. using PCR targeting the SSU rRNA gene. The overall prevalence of E. bieneusi and Blastocystis sp. was 5.4% (30/556) and 54.9% (305/556), respectively. The prevalence of E. bieneusi in dairy cattle, beef cattle, and buffaloes was 7.9% (13/165), 3.9% (11/283), and 5.6% (6/108), respectively. Eleven E. bieneusi genotypes were identified in this study, including six known genotypes, D (n = 10), I (n = 5), J (n = 4), IV (n = 4), N (n = 1), and BEB4 (n = 1), and five novel genotypes, JX-I to JX-V (n = 1), with genotype D as the predominant genotype in cattle. Phylogenetic analysis showed that six genotypes of E. bieneusi, D, IV, and JX-II to JX-V, were clustered into zoonotic group 1, whereas the remaining five genotypes belonged to group 2. Moreover, seven, seven, four, and five types were identified by multilocus sequence typing (MLST) at the MS1, MS3, MS4, and MS7 loci, respectively, forming three distinct multilocus genotypes (MLGs). In addition, the prevalence of Blastocystis sp. was 42.4% (70/165), 59.4% (168/283), and 62.0% (67/108) in dairy cattle, beef cattle, and buffaloes, respectively. Sequence analysis revealed that ST1, ST5, ST10, and ST14 of Blastocystis sp. were identified in these cattle, with ST10 being the major subtype. ST1 and ST5 are potential zoonotic subtypes. These findings have important implications for the control of E. bieneusi and Blastocystis sp. in cattle in Jiangxi Province.