Full-length transcriptome analysis and identification of transcript structures in Eimeria necatrix from different developmental stages by single-molecule real-time sequencing

Comparative proteomic analysis across the developmental stages of the Eimeria tenella

Eimeria tenella is the main pathogen responsible for coccidiosis in chickens. The life cycle of E. tenella is, arguably, the least complex of all Coccidia, with only one host. However, it presents different developmental stages, either in the environment or in the host and either intracellular or extracellular. Its signaling and metabolic pathways change with its different developmental stages. Until now, little is known about the developmental regulation and transformation mechanisms of its life cycle. In this study, protein profiles from the five developmental stages, including unsporulated oocysts (USO), partially sporulated (7 h) oocysts (SO7h), sporulated oocysts (SO), sporozoites (S) and second-generation merozoites (M2), were harvested using the label-free quantitative proteomics approach. Then the differentially expressed proteins (DEPs) for these stages were identified. A total of 314, 432, 689, and 665 DEPs were identified from the comparison of SO7h vs USO, SO vs SO7h, S vs SO, and M2 vs S, respectively. By conducting weighted gene coexpression network analysis (WGCNA), six modules were dissected. Proteins in blue and brown modules were calculated to be significantly positively correlated with the E. tenella developmental stages of sporozoites (S) and second-generation merozoites (M2), respectively. In addition, hub proteins with high intra-module degree were identified. Gene Ontology (GO) and Kyoto Encyclopedia of Gene and Genomes (KEGG) pathway enrichment analyses revealed that hub proteins in blue modules were involved in electron transport chain and oxidative phosphorylation. Hub proteins in the brown module were involved in RNA splicing. These findings provide new clues and ideas to enhance our fundamental understanding of the molecular mechanisms underlying parasite development.

Localization in vivo and in vitro confirms EnApiAP2 protein encoded by ENH_00027130 as a nuclear protein in Eimeria necatrix

Introduction

Apicomplexan AP2 family of proteins (ApiAP2) are transcription factors (TFs) that regulate parasite growth and development, but little is known about the ApiAP2 TFs in Eimeria spp. ENH_00027130 sequence is predicted to encode a Eimeria necatrix ApiAP2 protein (EnApiAP2).

Methods

The cDNAs encoding full-length and truncated EnApiAP2 protein were cloned and sequenced, respectively. Then, the two cDNAs were cloned into the pET28a(+) expression vector and expressed expressed in Escherichia coli BL21. The mouse polyclonal antibody (pAb) and monoclonal antibody (mAb) against recombinant EnApiAP2 (rEnApiAP2) and EnApiAP2tr (rEnApiAP2tr) were prepared and used to localize the native EnApiAP2 protein in E. necatrix, respectively. Finally, the recombinant pEGFP-C1-ΔNLS-EnApiAP2s (knockout of a nuclear localization sequence, NLS) and pEGFP-C1-EnApiAP2 plasmid were constructed and transfected into DF-1 cells, respectively, to further observe subcellular localization of EnApiAP2 protein.

Results

The EnApiAP2 gene had a size of 5019 bp and encoded 1672 amino acids, containing a conserved AP2 domain with a secondary structure consisting of an α-helix and three antiparallel β-strands. The rEnApiAP2 and rEnApiAP2tr were predominantly expressed in the form of inclusion bodies, and could be recognized by the 6×His tag mAb and the serum of convalescent chickens after infection with E. necatrix, respectively. The native EnApiAP2 protein was detected in sporozoites (SZ) and second generation merozoites (MZ-2) extracts, with a size of approximately 210 kDa. A quantitative real-time PCR (qPCR) analysis showed that the transcription level of EnApiAP2 was significantly higher in SZ than in MZ-2, third generation merozoites (MZ-3) and gametocytes (P<0.01). EnApiAP2 protein was localized in the nuclei of SZ, MZ-2 and MZ-3 of E. necatrix. The protein of EnApiAP2 was localized in the nucleus of the DF-1 cells, whereas the ΔNLS-EnApiAP2 was expressed in the cytoplasm, which further confirmed that EnApiAP2 is nucleoprotein.

Discussion

EnApiAP2 protein encoded by ENH_00027130 sequence was localized in the nucleus of E. necatrix parasites, and relied on the NLS for migration to DF-1 cell nucleus. The function of EnApiAP2 need further study.

Expression of a SAG protein with a CAP domain from Eimeria necatrix and its role in invasion and immunoprotection

Eimeria necatrix is a high pathogenic pathogen, which seriously endangers the poultry industry. The surface antigens (SAGs) of Apicomplexa are a kind of membrane protein anchored on the surface of the parasites through its carboxyl terminal glycosylphosphatidylinositol (GPI) structure. However, little is known about GPI-linked surface proteins in E. necatrix. In the present work, the E. necatrix sag gene (Ensag-CAP) was amplified and cloned for expression of the recombinant protein (rEnSAG-CAP). The full length Ensag-CAP gene was 813 bp, coding 270 amino acids with a predicated molecular weight of 28.86 kDa and contained a CAP domain with four sequence motifs CAP1, CAP2, CAP3 and CAP4. The rEnSAG-CAP was about 32 kDa and mainly expressed in a soluble form. Western blot analysis indicated that the rEnSAG-CAP could be recognized by anti-rEnSAG-CAP monoclonal antibody (anti-rEnSAG-CAP McAb) and the convalescent serum of chicken infected with E. necatrix. Native protein of EnSAG-CAP was detected in second-generation merozoites (MZ-2) using anti-rEnSAG-CAP polyclonal antibody (anti-rEnSAG-CAP pAb). The findings from the indirect immunofluorescence assay and enzyme digestion utilizing Bacillus cereus phosphoinositide-specific phospholipase C (PI-PLC) revealed that EnSAG-CAP predominantly localized at the surfaces of SZ and MZ-2 via a GPI anchor. It was observed that EnSAG-CAP can be cleaved from MZ-2 by PI-PLC. Real-time quantitative PCR (qPCR) analysis showed that transcript levels of Ensag-CAP in MZ-2 was significantly higher than that in SZ (P < 0.05). The anti-rEnSAG-CAP McAb in vitro could significantly inhibit the sporozoite invasion into MDBK cells (P < 0.01), which suggests that the protein might participate in sporozoite invasion into MDBK cells. rEnSAG-CAP afforded an immune protection against E. necatrix. The ACI value was 164.99 in the chickens immunized with 200 µg rEnSAG-CAP. Chickens immunized with rEnSAG-CAP had a significantly higher antigen-specific serum IgY response (P < 0.0001). The data indicates that EnSAG-CAP could serve as a potential candidate antigen for the development of a recombinant coccidiosis vaccine.

Comparative transcriptome profiling of Eimeria tenella in various developmental stages and functional analysis of an ApiAP2 transcription factor exclusively expressed during sporogony

BACKGROUND

The apicomplexan parasites Eimeria spp. are the causative agents of coccidiosis, a disease with a significant global impact on the poultry industry. The complex life cycle of Eimeria spp. involves exogenous (sporogony) and endogenous (schizogony and gametogony) stages. Unfortunately, the genetic regulation of these highly dynamic processes, particularly for genes involved in specific developmental phases, is not well understood.

METHODS

In this study, we used RNA sequencing (RNA-Seq) analysis to identify expressed genes and differentially expressed genes (DEGs) at seven time points representing different developmental stages of Eimeria tenella. We then performed K-means clustering along with co-expression analysis to identify functionally enriched gene clusters. Additionally, we predicted apicomplexan AP2 transcription factors in E. tenella using bioinformatics methods. Finally, we generated overexpression and knockout strains of ETH2_0411800 to observe its impact on E. tenella development.

RESULTS

In total, we identified 7329 genes that are expressed during various developmental stages, with 3342 genes exhibiting differential expression during development. Using K-means clustering along with co-expression analysis, we identified clusters functionally enriched for oocyte meiosis, cell cycle, and signaling pathway. Among the 53 predicted ApiAP2 transcription factors, ETH2_0411800 was found to be exclusively expressed during sporogony. The ETH2_0411800 overexpression and knockout strains did not exhibit significant differences in oocyst size or output compared to the parental strain, while the resulting ETH2_0411800 knockout parasite showed a relatively small oocyst output.

CONCLUSIONS

The findings of our research suggest that ETH2_0411800 is not essential for the growth and development of E. tenella. Our study provides insights into the gene expression dynamics and is a valuable resource for exploring the roles of transcription factor genes in regulating the development of Eimeria parasites.

Cloning and expression of a cDNA coding for Eimeria acervulina 25 kDa protein associated with oocyst and sporocyst walls

The purpose of this study was to characterize a gene named EAH 00033530 identified by RNAseq analysis of sporulating Eimeria acervulina oocysts and its encoded protein. Quantitative RT-PCR analysis revealed peak expression of EAH 00033530 mRNA early (3-6 h) in sporulation followed by downregulation at 12-24 h. The gene for EAH 00033530 was expressed in Escherichia coli as a 70 kDa polyHis fusion protein (rEAH 00033530). Antisera prepared against rEAH 00033530 protein identified in immunoblotting a native 25 kDa E. acervulina protein (Ea25) that was present in oocyst-sporocyst extracts after treatment with the reducing agent DTT. Immunofluorescence staining using anti-rEa25 localized the protein to both E. acervulina oocyst and sporocyst walls, but not to sporozoites. The protein may be produced during in vivo oocyst development because immunostaining of duodenal tissue from E. acervulina-infected chickens revealed oocyst wall expression. As observed by ELISA, rEa25 protein appears to elicit a humoral immune response in chickens infected with non-irradiated or radiation-attenuated E. acervulina oocysts.