Comparative transcriptome analysis of Eimeria necatrix third-generation merozoites and gametocytes reveals genes involved in sexual differentiation and gametocyte development

Transcriptomic Insights into the Developmental Dynamics of Eimeria acervulina: A Comparative Study of a Precocious Line and the Wild Type

Coccidiosis, a parasitic disease caused by single or multiple Eimeria species, leads to significant economic losses in the poultry industry. The Eimeria life cycle includes schizogony, gametogony, and sporogony. To investigate the dynamics of gene expression and regulatory networks during the development of Eimeria acervulina, we employed time-course transcriptomics to rigorously compare the gene expression patterns between a precocious line (PL) and the wild type (WT) of E. acervulina. The results revealed that the PL enters into gametogony 12 h earlier than the WT, and both the PL and WT exhibited distinct clustering patterns during the development phase. A weighted gene co-expression network analysis (WGCNA) identified genes specifically expressed at four distinct developmental stages, schizogony, gametogony, sporulated oocysts, and unsporulated oocysts, clarifying the key biological processes at each stage. This study used global transcriptome profiling to elucidate molecular variations throughout the E. acervulina life cycle, providing critical insights into molecular characterization and valuable resources for investigating other apicomplexan parasites of public health importance.

CRISPR-Cas9-based method for isolating microgametes of Eimeria tenella

Eimeria tenella infections are known to cause severe caecal damage and death of the infected chicken. Gamogony is an essential stage in E. tenella life cycle and in the establishment of coccidiosis. Prior research had extensively explored isolation and separation of the parasite gametes - microgamete (male) and macrogamete (female). However, there is little information on the efficient, highly purified and distinctly separated male and female gametes. In this study, we generated a genome editing line expressing mCherry fluorescent protein fused with GCS1 protein in E. tenella by using Toxoplasma gondii CRISPR-Cas9 system, flow cytometry and fluorescence microscopy. This allowed precise separation of E. tenella male and female gametes in the transgenic parasite population. The separation of male and female gametes would not only build on our understanding of E. tenella transmission, but it would also facilitate development of gametocidal compounds as drug targets for E. tenella infection.

Comparative proteomic analysis of wall-forming bodies and oocyst wall reveals the molecular basis underlying oocyst wall formation in Eimeria necatrix

BACKGROUND

The durable oocyst wall formed from the contents of wall-forming bodies (WFBs) protects Eimeria parasites from harsh conditions and enhances parasite transmission. Comprehending the contents of WFBs and proteins involved in oocyst wall formation is pivotal to understanding the mechanism of the oocyst wall formation and the search for novel targets to disrupt parasite transmission.

METHODS

Total proteins extracted from WFBs and the oocyst wall of Eimeria necatrix were subjected to comparative proteomic analysis using tandem mass tag in conjunction with liquid chromatography tandem-mass spectrometry techniques. After functional clustering analysis of the identified proteins, three proteins, including E. necatrix disulfide isomerase (EnPDI), thioredoxin (EnTrx) and phosphoglycerate kinase (EnPGK), were selected for further study to confirm their potential roles in oocyst wall formation.

RESULTS

A total of 3009 and 2973 proteins were identified from WFBs and the oocyst wall of E. necatrix, respectively. Among these proteins, 1102 were identified as differentially expressed proteins, of which 506 were upregulated and 596 downregulated in the oocyst wall compared to the WFBs. A total of 108 proteins, including compositional proteins of the oocyst wall, proteases, oxidoreductases, proteins involved in glycosylation, proteins involved in synthesis of the acid-fast lipid layer and proteins related to transport, were proposed to be involved in oocyst wall formation. The approximate molecular sizes of native EnPDI, EnTrx and EnPGK proteins were 55, 50 and 45 kDa, respectively. EnPDI was present in both type 1 and type 2 WFBs, EnTrx was present only in type 2 WFB2 and EnPGK was present only in type 1 WFBs, whereas all of them were localized to the outer layer of the oocyst wall, indicating that all of them participate in the formation of the oocyst wall.

CONCLUSIONS

To the best of our knowledge, this is the first report on the proteomes of WFBs and the oocyst wall of E. necatrix. The data obtained from this study form a basis for deciphering the molecular mechanisms underlying oocyst wall formation of Eimeria parasites. They also provide valuable resources for future studies on the development of novel therapeutic agents and vaccines aimed at combating coccidian transmission.

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.

Proteomic Analysis of Fractionated Eimeria tenella Sporulated Oocysts Reveals Involvement in Oocyst Wall Formation

Eimeria tenella is the most pathogenic intracellular protozoan parasite of the Eimeria species. Eimeria oocyst wall biogenesis appears to play a central role in oocyst transmission. Proteome profiling offers insights into the mechanisms governing the molecular basis of oocyst wall formation and identifies targets for blocking parasite transmission. Tandem mass tags (TMT)-labeled quantitative proteomics was used to analyze the oocyst wall and sporocysts of E. tenella. A combined total of 2865 E. tenella proteins were identified in the oocyst wall and sporocyst fractions; among these, 401 DEPs were identified, of which 211 were upregulated and 190 were downregulated. The 211 up-regulated DEPs were involved in various biological processes, including DNA replication, fatty acid metabolism and biosynthesis, glutathione metabolism, and propanoate metabolism. Among these proteins, several are of interest for their likely role in oocyst wall formation, including two tyrosine-rich gametocyte proteins (EtGAM56, EtSWP1) and two cysteine-rich proteins (EtOWP2, EtOWP6). Concurrently, 96 uncharacterized proteins may also participate in oocyst wall formation. The present study significantly expands our knowledge of the proteome of the oocyst wall of E. tenella, thereby providing a theoretical basis for further understanding of the biosynthesis and resilience of the E. tenella oocyst wall.

Unravelling the sexual developmental biology of Cystoisospora suis, a model for comparative coccidian parasite studies

Introduction

The apicomplexan parasite Cystoisospora suis has global significance as an enteropathogen of suckling piglets. Its intricate life cycle entails a transition from an asexual phase to sexual development, ultimately leading to the formation of transmissible oocysts.

Methods

To advance our understanding of the parasite's cellular development, we complemented previous transcriptome studies by delving into the proteome profiles at five distinct time points of in vitro cultivation through LC/MS-MS analysis.

Results

A total of 1,324 proteins were identified in the in vitro developmental stages of C. suis, and 1,082 proteins were identified as significantly differentially expressed. Data are available via ProteomeXchange with identifier PXD045050. We performed BLAST, GO enrichment, and KEGG pathway analyses on the up- and downregulated proteins to elucidate correlated events in the C. suis life cycle. Our analyses revealed intriguing metabolic patterns in macromolecule metabolism, DNA- and RNA-related processes, proteins associated with sexual stages, and those involved in cell invasion, reflecting the adaptation of sexual stages to a nutrient-poor and potentially stressful extracellular environment, with a focus on enzymes involved in metabolism and energy production.

Discussion

These findings have important implications for understanding the developmental biology of C. suis as well as other, related coccidian parasites, such as Eimeria spp. and Toxoplasma gondii. They also support the role of C. suis as a new model for the comparative biology of coccidian tissue cyst stages.

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.

What Do We Know about Surface Proteins of Chicken Parasites Eimeria?

Poultry is the first source of animal protein for human consumption. In a changing world, this sector is facing new challenges, such as a projected increase in demand, higher standards of food quality and safety, and reduction of environmental impact. Chicken coccidiosis is a highly widespread enteric disease caused by Eimeria spp. which causes significant economic losses to the poultry industry worldwide; however, the impact on family poultry holders or backyard production-which plays a key role in food security in small communities and involves mainly rural women-has been little explored. Coccidiosis disease is controlled by good husbandry measures, chemoprophylaxis, and/or live vaccination. The first live vaccines against chicken coccidiosis were developed in the 1950s; however, after more than seven decades, none has reached the market. Current limitations on their use have led to research in next-generation vaccines based on recombinant or live-vectored vaccines. Next-generation vaccines are required to control this complex parasitic disease, and for this purpose, protective antigens need to be identified. In this review, we have scrutinised surface proteins identified so far in Eimeria spp. affecting chickens. Most of these surface proteins are anchored to the parasite membrane by a glycosylphosphatidylinositol (GPI) molecule. The biosynthesis of GPIs, as well as the role of currently identified surface proteins and interest as vaccine candidates has been summarised. The potential role of surface proteins in drug resistance and immune escape and how these could limit the efficacy of control strategies was also discussed.

Examination of gametocyte protein 22 localization and oocyst wall formation in Eimeria necatrix using laser confocal microscopy and scanning electron microscopy

BACKGROUND

Eimeria parasite infection occurs via ingestion of oocysts. The robust, bilayer oocyst wall is formed from the contents of wall-forming bodies (WFBs), WFB1 and WFB2, located exclusively in macrogametocytes. Eimeria necatrix gametocyte proteins 22 and 59 (EnGAM22 and EnGAM59) have been found to localize to WFBs and the oocyst wall. However, the exact localization of these two proteins is not clear.

METHODS

WFBs of E. necatrix were extracted from purified gametocytes using a cutoff filter and the extracts of purified WFBs and gametocytes were analyzed using sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) and immunoblotting. Then, the localization of EnGAM22 and EnGAM59 proteins was determined using an indirect immunofluorescence assay. Finally, the development of macrogametocytes and the oocyst wall of E. necatrix was analyzed using laser confocal microscopy and scanning electron microscopy.

RESULTS

Purified WFBs had the same shape and size as those observed in macrogametocytes. EnGAM22 protein localized to WFB1, whereas EnGAM59 protein localized to WFB2. Both EnGAM22 and EnGAM59 native proteins were detected in the extracts of WFBs and gametocytes. The outer layer of the oocyst wall was formed by the release of the contents of WFB1 at the surface of the macrogametocyte to form an anti-EnGAM22 positive layer. WFB2 then appeared to give rise to the inner layer, which was anti-EnGAM59 positive.

CONCLUSIONS

EnGAM22 and EnGAM59 proteins localized to WFB1 and WFB2 and were involved in the formation of the outer and inner layers of the oocyst wall of E. necatrix, respectively. The processes of macrogametogenesis and oocyst wall formation of E. necatrix are similar to other Eimeria parasites. The anti-EnGAM22 antibody could be used as a tool to track the transport and secretion of proteins in WFB1 during oocyst development.

Evidence of high-efficiency cross fertilization in Eimeria acervulina revealed using two lines of transgenic parasites

Eimeria species are apicomplexan parasites with a direct life cycle consisting of a replicative phase involving multiple rounds of asexual replication in the intestine or other organs including kidneys, liver, and gallbladder, depending on the species, followed by a sexual phase or gamogony involving the development and fertilization of gametes, an essential process for Eimeria transmission. Recent advances in the genetic manipulation of these parasites made it possible to conduct genetic crosses combined with genomic approaches to elucidate the genetic determinants of Eimeria development, virulence, drug resistance, and immune evasion. Here, we employed genetic techniques to generate two transgenic Eimeria acervulina lines, EaGAM56 and EaHAP2, each expressing two unique fluorescent proteins, with one controlled by a constitutive promotor for cross-efficiency analysis and the other by a male or female gametocyte stage-specific promoter to observe sexual development. The expression of fluorescent proteins in the transgenic lines was analyzed in different developmental stages of the E. acervulina life cycle by immunoblotting and by examination of frozen sections using fluorescence microscopy. The effect of infective doses on cross-fertilization was further investigated by conducting several genetic crosses between the two transgenic lines at different doses and ratios. Two transgenic lines expressing constitutive and gametocyte-specific fluorescence proteins were generated and characterized. These transgenic parasites display synchronous development in chickens, comparable with that of the wild type. Genetic crosses between the two transgenic parasites showed that a high rate of oocysts co-expressing the two reporters could be achieved following inoculation with high doses of infective oocysts. We further showed that the proportion of co-transfected oocysts can be modulated by altering the ratio of the transgenic parental lines. Higher infective doses and similar numbers of functional gametocytes from the parents increase the rate of cross-fertilization. Our data highlight the usefulness of genetic manipulation and fluorescently-labeled transgenic gametocytes as tools to study Eimeria development and to elucidate the factors that modulate sexual development. This work sets the stage for the implementation of novel approaches to investigate other aspects of Eimeria pathogenesis, virulence, and drug susceptibility and resistance.

The transcriptome from asexual to sexual in vitro development of Cystoisospora suis (Apicomplexa: Coccidia)

The apicomplexan parasite Cystoisospora suis is an enteropathogen of suckling piglets with woldwide distribution. As with all coccidian parasites, its lifecycle is characterized by asexual multiplication followed by sexual development with two morphologically distinct cell types that presumably fuse to form a zygote from which the oocyst arises. However, knowledge of the sexual development of C. suis is still limited. To complement previous in vitro studies, we analysed transcriptional profiles at three different time points of development (corresponding to asexual, immature and mature sexual stages) in vitro via RNASeq. Overall, transcription of genes encoding proteins with important roles in gametes biology, oocyst wall biosynthesis, DNA replication and axonema formation as well as proteins with important roles in merozoite biology was identified. A homologue of an oocyst wall tyrosine rich protein of Toxoplasma gondii was expressed in macrogametes and oocysts of C. suis. We evaluated inhibition of sexual development in a host-free culture for C. suis by antiserum specific to this protein to evaluate whether it could be exploited as a candidate for control strategies against C. suis. Based on these data, targets can be defined for future strategies to interrupt parasite transmission during sexual development.

Identification and Protective Efficacy of Eimeria tenella Rhoptry Kinase Family Protein 17

Simple Summary

Approximately 8000 genes of Eimeria tenella have been found by genome sequencing, whereas very few data are currently available regarding E. tenella rhoptry kinase family proteins. In this study, the coding sequence of the rhoptry kinase family protein 17 of E. tenella (EtROP17) was cloned and expressed in Escherichia coli, and then the protective efficacy of the recombinant EtROP17 (rEtROP17) was assessed in chickens. Sequence analysis showed that a single base difference at position 1901 of the ROP17 of the SD-01 strain was observed compared with that of the Houghton strain. EtROP17 was expressed in the merozoite stage of E. tenella and may be a potential vaccine candidate against coccidiosis.

Abstract

Eimeria tenella encodes a genome of approximately 8000 genes. To date, however, very few data are available regarding E. tenella rhoptry kinase family proteins. In the present study, the gene fragment encoding the mature peptide of the rhoptry kinase family protein 17 of E. tenella (EtROP17) was amplified by PCR and expressed in E. coli. Then, we generated polyclonal antibodies that recognize EtROP17 and investigated the expression of EtROP17 in the merozoite stage of E. tenella by immunofluorescent staining and Western blot analysis. Meanwhile, the protective efficacy of rEtROP17 against E. tenella was evaluated in chickens. Sequencing analysis showed that a single base difference at sequence position 1901 was observed between the SD-01 strain and the Houghton strain. EtROP17 was expressed in the merozoite stage of E. tenella. The results of the animal challenge experiments demonstrated that vaccination with rEtROP17 significantly reduced cecal lesions and oocyst outputs compared with the challenged control group. Our findings indicate that EtROP17 could serve as a potential candidate for developing a new vaccine against E. tenella.

Eimeria proteins: order amidst disorder

Apicomplexans are important pathogens that cause severe infections in humans and animals. The biology and pathogeneses of these parasites have shown that proteins are intrinsically modulated during developmental transitions, physiological processes and disease progression. Also, proteins are integral components of parasite structural elements and organelles. Among apicomplexan parasites, Eimeria species are an important disease aetiology for economically important animals wherein identification and characterisation of proteins have been long-winded. Nonetheless, this review seeks to give a comprehensive overview of constitutively expressed Eimeria proteins. These molecules are discussed across developmental stages, organelles and sub-cellular components vis-à-vis their biological functions. In addition, hindsight and suggestions are offered with intention to summarise the existing trend of eimerian protein characterisation and to provide a baseline for future studies.

Sexual Development in Non-Human Parasitic Apicomplexa: Just Biology or Targets for Control?

Simple Summary

Cellular reproduction is a key part of the apicomplexan life cycle, and both mitotic (asexual) and meiotic (sexual) cell divisions produce new individual cells. Sexual reproduction in most eukaryotic taxa indicates that it has had considerable success during evolution, and it must confer profound benefits, considering its significant costs. The phylum Apicomplexa consists of almost exclusively parasitic single-celled eukaryotic organisms that can affect a wide host range of animals from invertebrates to mammals. Their development is characterized by complex steps in which asexual and sexual replication alternate and the fertilization of a macrogamete by a microgamete results in the formation of a zygote that undergoes meiosis, thus forming a new generation of asexual stages. In apicomplexans, sex is assumed to be induced by the (stressful) condition of having to leave the host, and either gametes or zygotes (or stages arising from it) are transmitted to a new host. Therefore, sex and meiosis are linked to parasite transmission, and consequently dissemination, which are key to the parasitic lifestyle. We hypothesize that improved knowledge of the sexual biology of the Apicomplexa will be essential to design and implement effective transmission-blocking strategies for the control of the major parasites of this group.

Abstract

The phylum Apicomplexa is a major group of protozoan parasites including gregarines, coccidia, haemogregarines, haemosporidia and piroplasms, with more than 6000 named species. Three of these subgroups, the coccidia, hemosporidia, and piroplasms, contain parasites that cause important diseases of humans and animals worldwide. All of them have complex life cycles involving a switch between asexual and sexual reproduction, which is key to their development. Fertilization (i.e., fusion of female and male cells) results in the formation of a zygote that undergoes meiosis, forming a new generation of asexual stages. In eukaryotes, sexual reproduction is the predominant mode of recombination and segregation of DNA. Sex is well documented in many protist groups, and together with meiosis, is frequently linked with transmission to new hosts. Apicomplexan sexual stages constitute a bottleneck in the life cycle of these parasites, as they are obligatory for the development of new transmissible stages. Consequently, the sexual stages represent attractive targets for vaccination. Detailed understanding of apicomplexan sexual biology will pave the way for the design and implementation of effective transmission-blocking strategies for parasite control. This article reviews the current knowledge on the sexual development of Apicomplexa and the progress in transmission-blocking vaccines for their control, their advantages and limitations and outstanding questions for the future.

Dynamically expressed genes provide candidate viability biomarkers in a model coccidian

Eimeria parasites cause enteric disease in livestock and the closely related Cyclosporacayetanensis causes human disease. Oocysts of these coccidian parasites undergo maturation (sporulation) before becoming infectious. Here, we assessed transcription in maturing oocysts of Eimeria acervulina, a widespread chicken parasite, predicted gene functions, and determined which of these genes also occur in C. cayetanensis. RNA-Sequencing yielded ~2 billion paired-end reads, 92% of which mapped to the E. acervulina genome. The ~6,900 annotated genes underwent temporally-coordinated patterns of gene expression. Fifty-three genes each contributed >1,000 transcripts per million (TPM) throughout the study interval, including cation-transporting ATPases, an oocyst wall protein, a palmitoyltransferase, membrane proteins, and hypothetical proteins. These genes were enriched for 285 gene ontology (GO) terms and 13 genes were ascribed to 17 KEGG pathways, defining housekeeping processes and functions important throughout sporulation. Expression differed in mature and immature oocysts for 40% (2,928) of all genes; of these, nearly two-thirds (1,843) increased their expression over time. Eight genes expressed most in immature oocysts, encoding proteins promoting oocyst maturation and development, were assigned to 37 GO terms and 5 KEGG pathways. Fifty-six genes underwent significant upregulation in mature oocysts, each contributing at least 1,000 TPM. Of these, 40 were annotated by 215 GO assignments and 9 were associated with 18 KEGG pathways, encoding products involved in respiration, carbon fixation, energy utilization, invasion, motility, and stress and detoxification responses. Sporulation orchestrates coordinated changes in the expression of many genes, most especially those governing metabolic activity. Establishing the long-term fate of these transcripts in sporulated oocysts and in senescent and deceased oocysts will further elucidate the biology of coccidian development, and may provide tools to assay infectiousness of parasite cohorts. Moreover, because many of these genes have homologues in C. cayetanensis, they may prove useful as biomarkers for risk.

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

Background

Eimeria necatrix is one of the most pathogenic parasites, causing high mortality in chickens. Although its genome sequence has been published, the sequences and complete structures of its mRNA transcripts remain unclear, limiting exploration of novel biomarkers, drug targets and genetic functions in E. necatrix.

Methods

Second-generation merozoites (MZ-2) of E. necatrix were collected using Percoll density gradients, and high-quality RNA was extracted from them. Single-molecule real-time (SMRT) sequencing and Illumina sequencing were combined to generate the transcripts of MZ-2. Combined with the SMRT sequencing data of sporozoites (SZ) collected in our previous study, the transcriptome and transcript structures of E. necatrix were studied.

Results

SMRT sequencing yielded 21,923 consensus isoforms in MZ-2. A total of 17,151 novel isoforms of known genes and 3918 isoforms of novel genes were successfully identified. We also identified 2752 (SZ) and 3255 (MZ-2) alternative splicing (AS) events, 1705 (SZ) and 1874 (MZ-2) genes with alternative polyadenylation (APA) sites, 4019 (SZ) and 2588 (MZ-2) fusion transcripts, 159 (SZ) and 84 (MZ-2) putative transcription factors (TFs) and 3581 (SZ) and 2039 (MZ-2) long non-coding RNAs (lncRNAs). To validate fusion transcripts, reverse transcription-PCR was performed on 16 candidates, with an accuracy reaching up to 87.5%. Sanger sequencing of the PCR products further confirmed the authenticity of chimeric transcripts. Comparative analysis of transcript structures revealed a total of 3710 consensus isoforms, 815 AS events, 1139 genes with APA sites, 20 putative TFs and 352 lncRNAs in both SZ and MZ-2.

Conclusions

We obtained many long-read isoforms in E. necatrix SZ and MZ-2, from which a series of lncRNAs, AS events, APA events and fusion transcripts were identified. Information on TFs will improve understanding of transcriptional regulation, and fusion event data will greatly improve draft versions of gene models in E. necatrix. This information offers insights into the mechanisms governing the development of E. necatrix and will aid in the development of novel strategies for coccidiosis control.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-021-05015-7.

Recent achievements and doors opened for coccidian parasite research and development through transcriptomics of enteric sexual stages

The Coccidia is the largest group of parasites within the Apicomplexa, a phylum of unicellular, obligate parasites characterized by the possession of an apical complex of organelles and structures in the asexual stages of their life cycles, as well as by a sexual reproductive phase that occurs enterically in host animals. Coccidian sexual reproduction involves morphologically distinct microgametes and macrogametes that combine to form a diploid zygote and, ultimately, following meiosis and mitosis, haploid, infectious sporozoites, inside sporocysts within an oocyst. Recent transcriptomic analyses have identified genes involved in coccidian sexual stage development and reproduction, including genes encoding for microgamete- and macrogamete-specific proteins with roles in gamete motility, fusion and fertilization, and in the formation of the resilient oocyst wall that allows coccidians to persist for long periods in the environment. Transcriptomics has also provided important clues about the regulation of gene expression in the transformation of parasites from one developmental stage to the next, a complex sequence of events that may involve transcription factors such as the apicomplexan Apetala2 (ApiAP2) family, alternative splicing, regulatory RNAs and MORC (a microrchida homologue and regulator of sexual stage development in Toxoplasma gondii). The molecular dissection of coccidian sexual development and reproduction by transcriptomic analyses may lead to the development of novel transmission-blocking strategies.

Characterization of Cystoisospora suis sexual stages in vitro

BACKGROUND

The porcine coccidium Cystoisospora suis is characterized by a complex life-cycle during which asexual multiplication is followed by sexual development with two morphologically distinct cell types, the micro- and macrogametes. Genes related to the sexual stages and cell cycle progression were previously identified in related Apicomplexa. Dynein light chain type 1 and male gamete fusion factor HAP2 are restricted to microgametes. Tyrosine-rich proteins and oocyst wall proteins are a part of the oocyst wall. The Rad51/Dmc1-like protein and Nima-related protein kinases are associated with the cell cycle and fertilization process. Here, the sexual stages of C. suis were characterized in vitro morphologically and for temporal expression changes of the mentioned genes to gain insight into this poorly known phase of coccidian development.

METHODS

Sexual stages of C. suis developing in vitro in porcine intestinal epithelial cells were examined by light and electron microscopy. The transcriptional levels of genes related to merozoite multiplication and sexual development were evaluated by quantitative real-time PCR at different time points of cultivation. Transcription levels were compared for parasites in culture supernatants at 6-9 days of cultivation (doc) and intracellular parasites at 6-15 doc.

RESULTS

Sexual stage of C. suis was detected during 8-11 doc in vitro. Microgamonts (16.8 ± 0.9 µm) and macrogamonts (16.6 ± 1.1 µm) are very similar in shape and size. Microgametes had a round body (3.5 ± 0.5 µm) and two flagella (11.2 ± 0.5 µm). Macrogametes were spherical with a diameter of 12.1 ± 0.5 µm. Merozoite gene transcription peaked on 10 doc and then declined. Genes related to the sexual stages and cell cycle showed an upregulation with a peak on 13 doc, after which they declined.

CONCLUSIONS

The present study linked gene expression changes to the detailed morphological description of C. suis sexual development in vitro, including fertilization, meiosis and oocyst formation in this unique model for coccidian parasites. Following this process at the cellular and molecular level will elucidate details on potential bottlenecks of C. suis development (applicable for coccidian parasites in general) which could be exploited as a novel target for control.

Further confirmation of second- and third-generation Eimeria necatrix merozoite DEGs using suppression subtractive hybridization

In our previous study, we obtained a large number of differentially expressed genes (DEGs) between second-generation merozoites (MZ-2) and third-generation merozoites (MZ-3) of Eimeria necatrix using RNA sequencing (RNA-seq). Here, we report two subtractive cDNA libraries for MZ₂ (forward library) and MZ₃ (reverse library) that were constructed using suppression subtractive hybridization (SSH). PCR amplification revealed that the MZ₂ and MZ₃ libraries contained approximately 96.7% and 95% recombinant clones, respectively, and the length of the inserted fragments ranged from 0.5 to 1.5 kb. A total of 106 and 111 unique sequences were obtained from the MZ₂ and MZ₃ libraries, respectively, and were assembled into 13 specific consensus sequences (contigs or genes) (5 from MZ₂ and 8 from MZ₃). The qRT-PCR results revealed that 11 out of 13 genes were differentially expressed between MZ-2 and MZ-3. Of 13 genes, 11 genes were found in both SSH and our RNA-seq data and displayed a similar expression trend between SSH and RNA-seq data, and the remaining 2 genes have not been reported in both E. necatrix genome and our RNA-seq data. Among the 11 genes, the expression trends of 8 genes were highly consistent between SSH and our RNA-seq data. These DEGs may provide specialized functions related to the life-cycle transitions of Eimeria species.