The Cryptosporidium Parvum Transcriptome during In Vitro Development

Abstract

Cryptosporidiosis is caused by an obligate intracellular parasite that has eluded global transcriptional or proteomic analysis of the intracellular developmental stages. The transcript abundance for 3,302 genes (87%) of the Cryptosporidium parvum protein coding genome was elucidated over a 72 hr infection within HCT8 cells using Real Time-PCR. The parasite had detectable transcription of all genes in vitro within at least one time point tested, and adjacent genes were not co-regulated. Five genes were not detected within the first 24 hr of infection, one containing two AP2 domains. The fewest genes detected were at 2 hr post infection, while 30% (985) of the genes have their highest expression at 48 and/or 72 hr. Nine expression clusters were formed over the entire 72 hr time course and indicate patterns of transcriptional increases at each of the 7 time points collected except 36 hr, including genes paralleling parasite 18S rRNA transcript levels. Clustering within only the first 24 hr of infection indicates spikes in expression at each of the 4 time points, a group paralleling 18S rRNA transcript levels, and a cluster with peaks at both 6 and 24 hr. All genes were classified into 18 functional categories, which were unequally distributed across clusters. Expression of metabolic, ribosomal and proteasome proteins did not parallel 18S rRNA levels indicating distinct biochemical profiles during developmental stage progression. Proteins involved in translation are over-represented at 6 hr, while structural proteins are over-represented at 12 hr. Standardization methods identified 107 genes with <80% at a single of its total expression at a single time point over 72 hr. This comprehensive transcriptome of the intracellular stages of C. parvum provides insight for understanding its complex development following parasitization of intestinal epithelial cells.

About the Authors:

Mary J. Mauzy

Shinichiro Enomoto

Cheryl A. Lancto

Mitchell S. Abrahamsen

Mark S. Rutherford

* E-mail: [email protected]

Affiliation: Department of Veterinary and Biomedical Sciences, College of Veterinary Medicine, University of Minnesota, St. Paul, Minnesota, United States of America

Introduction

Cryptosporidium species are global contaminants of surface water and are the second leading cause of human gastrointestinal illness in the United States. Reported incidence is highest in children, yet seroprevalence is significant in all age categories [1], [2]. Due to its resistance to standard water chlorine disinfection, Cryptosporidium is a public health concern and a potential water-borne bioterrorism agent due to its low infectious dose (as low as 10 oocysts) and its ability to be stably delivered to the human population en masse [3]. Illness varies from profuse, self-limiting diarrhea to life threatening malabsorption and dehydration depending on immune status. Effective therapeutics have not been formulated because the eukaryotic parasite has a condensed genome lacking many of the traditional drug targets [4]. Most of the remaining genes have remained functionally uncharacterized, thereby limiting pharmacological targets [5], [6].

Apicomplexa are parasitic eukaryotes noted for undergoing both asexual and sexual replicative stages during their life cycle. Cryptosporidium spp. complete their life cycle within a single host utilizing only epithelial cells. Ingestion of an oocyst results in excystation of four sporozoites in the gastrointestinal tract. Following attachment to the host epithelial cell, the parasite resides within an intracellular but extracytoplasmic parasitophorous vacuole derived from the host cell membrane. Therein, sporozoites mature into trophozoites which then progress through asexual replication (3–4 rounds of mitosis) in ∼24 hr to form type 1 meronts that release 6–8 merozoites. These merozoites infect new epithelial cells to either repeat asexual replication generating more type 1 meronts, or through an unknown, presumably environmental trigger progress through sexual development resulting in type 2 meronts. Type 2 meronts release 4 merozoites that develop either into micro- or macro-gamonts that continue through sporogony to produce infectious oocysts that are excreted in high numbers in the feces. Sexual development of Cryptosporidium has been morphologically described in vitro, indicating the environmental switch is present, yet monolayer cultures are unable to maintain continuous propagation [7], [8].

Little is known regarding how Cryptosporidium regulates developmental stage progression or the different cellular requirements required by each stage. The intricate enclosure of the parasite in a parasitophorous vacuole on the host cell surface has limited high-throughput analyses of the transcriptional or proteomic repertoire of Cryptosporidium to the sporozoite stage [9]–[12]. Morphological characterization of the parasite is also limited because many stages are macroscopically identical even though the parasite utilizes both asexual and sexual developmental progression. To gain insight into Cryptosporidium biology and development, we completed a genome-wide transcriptome analysis over a 72 hr in vitro infection of human epithelial cells using the zoonotic species, C. parvum. Real time-PCR (RT-PCR) for 3,302 genes (87% of the protein coding genes) indicated that each gene has detectable transcription in at least one time point assessed. Further characterization of gene expression indicates 9 clusters over the entire 72 hr time course, while the initial asexual replication cycle (2–24 hr) produced 6 clusters with both clustering techniques distinguishing genes with coordinate patterns of expression. Additionally, we identified genes for which mRNA levels spiked at single time points, suggestive of fluctuations in transcript density dependent on the parasitic developmental stage present. Herein we describe the first comprehensive temporal transcriptome analysis of C. parvum during in vitro development in epithelial cells.

Materials and Methods

C. parvum parasites

C. parvum oocysts (Iowa strain) were purchased (Bunch Grass Farm, Drury, ID) and stored in antibiotics at 4°C for less than 2 months prior to use. Before infecting the cells, oocysts were surface sterilized by treatment with a 33% bleach (sodium hypochlorite) solution (1 ml/3×107 oocysts) on ice for 7 minutes, washed profusely with Hank's buffered saline solution (HBSS), and stored in HBSS at 4°C overnight.

C. parvum infection model

Human ileocecal adenocarcinoma (HCT8, ATCC CCL-244; American Type Culture Collection, Rockville, MD.) cells were cultured in a humidified incubator at 37°C in an atmosphere containing 5% CO2 on 22 mm glass coverslips contained in 6 well plates or 10 cm2 dishes until confluency was reached in 6 days. Culture media [RPMI 1640 medium supplemented with 10% fetal bovine serum (FBS), 1 mM sodium pyruvate, 50 U/ml penicillin G, 50 U/ml streptomycin, and 0.25 µg/ml amphotericin B (Fungizone) (pH 7.4)] was changed every 24–48 hours as required. After six days the culture media was replaced with Cp-Up media [RPMI 1640 media containing 10% fetal bovine serum, 15 mM HEPES, 50 mM glucose, 10 µg/ml bovine insulin, 35 µg/ml ascorbic acid, 1.0 µg/ml folic acid, 4.0 µg/ml 4-aminobenzoic acid, 2.0 µg/ml calcium pantothenate, 50 U/ml penicillin G, 50 U/ml streptomycin, and 0.25 µg/ml amphotericin B (Fungizone) (pH 7.4)] and cultured for an additional four days with media refreshment as necessary [13], [14]. Prior to infection, oocysts were warmed to room temperature and inoculated onto culture monolayers at 1.5×106 oocysts/well for coverslips or 2−2.5×107 oocysts/10 cm2 dish as previously described [14]. Following a 2 hr excystation period, the unexcysted oocysts and free sporozoites were washed from monolayers with warm HBSS and cultures were incubated in Cp-Up media for the specified time points at 37°C [7], [14]. Infection rate was 80%–90% at 24 hr depending on the batch and storage period of oocysts. Cells without infection were used as mock controls. At the specified time periods, monolayers were washed once briefly in warm HBSS and the coverslips were fixed in PBS containing 4% paraformaldehyde for 15 minutes followed by four washings in PBS and stored at 4°C. The 10 cm2 dishes were rinsed once in PBS then lysed in TRIzol (Invitrogen) and stored at −80°C until RNA isolation.

Indirect immunofluorescence

C. parvum infected, paraformaldehyde fixed HCT8 cell coverslips were permeabilized by treating with 0.15% Triton-X 100 in PBS for 10 min. Coverslips were washed and nonspecific binding sites were blocked for 40 min by using 2.5% fetal calf serum +2.5% goat serum. Coverslips were incubated for 1 hr with either Cp-65.10, a pan monoclonal antibody that recognizes all C. parvum life-stages, or a secondary control antibody. Following washing, the coverslips were incubated with AlexaFluor 568 (Invitrogen). The coverslips were washed, counterstained with DAPI and mounted to slides. Photomicrographs were captured at 40X using a Nikon microscope equipped with a high resolution Zeiss Axiovert 2000, with an Axiocam attachment.

RNA isolation.

RNA was isolated from samples stored in TRIzol following the manufacturers protocol. In brief, 0.2 ml chloroform was added per 1 ml of TRIzol used, mixed briefly and incubated for 3 minutes at room temperature. The aqueous layer was recovered after separation via centrifugation at 10,000 rpm for 20 minutes. Five hundred µl of isopropanol was added per 1 ml TRIzol used, incubated at room temperature for 10 minutes and removed by centrifugation at 10,000 rpm for 20 minutes. The resulting pellet was washed first in 75% ethanol (1 ml per 1 ml TRIzol used), then 70% ethanol (1 ml per 1 ml TRIzol used), with pelleting of RNA at 10,000 rpm after each wash. After removal of the final wash, RNA was resuspended in molecular grade water at a concentration of 1–2 µg/µl. RNA recovery and integrity was verified on a formaldehyde gel prior to use.

DNase treatment.

DNA contamination was degraded utilizing the Turbo DNA-free kit (Applied Biosystems) following the manufacturers recommendations. 50 µl reactions were constructed with 30 µg of RNA and 5 µl of 10× buffer. 2 µl of Turbo DNase was incubated with the sample at 37°C for 30 minutes, followed by a 2 minute room temperature incubation with 10 µl of inactivating reagent. The inactivating reagent was removed via centrifugation at 10,000× g for 90 seconds. The RNA was collected and quantified for cDNA synthesis.

cDNA synthesis.

To obtain enough cDNA for the whole genome transcriptome and to reduce qRT-PCR variability inherent to cDNA synthesis [15], [16], the cDNAs for each time point and replicate were made in multiple 20 µl reaction volumes, and then replicate cDNA reactions for each time point were individually pooled for each of the four separate time courses. Aliquots were stored at −80°C until used in the qPCR reactions. cDNA synthesis was accomplished using Superscript III cDNA synthesis kit (Invitrogen), with the following modifications from the manufacturers protocol. Two micrograms of DNased RNA was utilized, with 200 ng of random hexamer primers. Denaturing was done at 65°C for 5 minutes, followed by the synthesis reaction with incubations at 25°C for 10 minutes, 50°C for 60 minutes with the reaction being terminated at 85°C for 20 minutes. 1 U of recombinant RNase H (New England BioLabs) was incubated with each sample for 20 minutes at 37°C to remove template RNA. Each cDNA synthesis reaction included a negative control lacking reverse transcriptase to verify proper DNase treatment. cDNA synthesis, removal of genomic DNA contaminants (Figure S1), and DNA degradation was verified using C. parvum 18S ribosomal RNA primers (F-CTCCACCAACTAAGAACGGCC, R-TAGAGATTGGAGGTTGTTCCT).

Real Time PCR (RT-PCR).

Twenty microliter reactions were constructed using a 4 µl of template from a 1∶100 dilution of synthesized cDNA (2 µg RNA starting concentration), 0.1 µM primer pairs, and 2× AccuQuant SYBR Green SuperMix, with Low Rox (Quanta BioSciences, Inc.). RT-PCR was run on Stratagene Mx3000P or Mx3005P thermocyclers, with cycle conditions of: 1 cycle of initial denaturation (95°C for 1 minute), 45 PCR cycles (95°C for 1 minute, 58°C for 30 sec. and ramping to 60°C for 30 sec.). Amplification was followed by melt curve generation via 1 minute at 95°C ramping down to 57°C in 2 degree intervals followed by ramping back up to 95°C for 30 seconds. All primers (Table S2) were designed using Primer3(http://frodo.wi.mit.edu/primer3/) and tested on 1) C. parvum genomic DNA to confirm amplification and product size, and 2) HCT8 cDNA to disqualify primers amplifying host transcripts indistinguishable from parasite products. The C. parvum genome has few gene paralogues, and few Cryptosporidium genes have introns [12], [17], which allowed primer testing on Cryptosporidium genomic DNA. Successful primer pairs (3,302) were used to query the HCT8/Cryptosporidium cDNA constructed from total RNA isolated at 2, 6, 12, 24, 36, 48 and 72 hours (hr) post addition of oocysts to the cell culture. Confirmation of successful 18S rRNA amplification utilizing both the amplification and disassociation curves for each reaction allowed for inclusion of the data for subsequent analysis. For each gene of interest, the amplification and disassociation curves were also verified for consistency, and wells were individually disqualified from analysis if the Ct was >40, the amplification curve slope was skewed from average, or the product amplified did not have the correct melt temperature. The unprocessed RT-PCR experimental data are available on cryptogenome.umn.edu, on cryptoDB (cryptodb.org), and will also be available in the NCBI repository.

The raw fluorescence data was curve-fitted to calculate the concentration of template in each sample.

The fluorescence (F(c) as a function of the PCR cycle (c) was fitted to a four point logistic curve:where the four parameters represent, Fb; fluorescence base, Fmax; fluorescence maximum, Ch; cycle at midpoint of the rise, and β; slope. The initial fluorescence (I0) was calculated by extrapolating back to cycle  = 0. This results in a value of transcript per 18S rRNA concentration per 2 µg of cDNA.

Hence the relative DNA concentration is obtained by normalizing the I0 of a gene to that of 18S ribosomal RNA (Table S1) as previously validated for RT-PCR [8], [18], [19]. The median gene expression was calculated using 4 independent infection time courses and standardization to the total transcript detected over the entire time course per gene.

Functional Annotation.

As described in Jensen et al [20], all genes were assigned to one of 18 functional categories based on gene annotation comparison between cryptogenome.umn.edu and cryptodb.org as well as literature review. The 18 functional categories used were (categories modified for this analysis are explained in parenthesis): DNA-associated, Interacting, Metabolism, Oxidative Stress, Phosphorylation, Protease related, Protein folding, Protein transport/modification, RNA-associated, Transcription, Translation, Transporter, Other, Hypothetical, Conserved Hypothetical, Apicomplexan conserved (genes with or without known function that are conserved across Aplicomplexa), and Cryptosporidium (genes of interest within Cryptosporidium, such as mucins, COWPs, TSPs, TRAPs and MEDLE family).

Computational analysis.

Cluster analysis was done using the DIvisive ANAlysis Clustering (DIANA) algorithm available in the ‘cluster’ package of R [21], [22]. Data was analyzed as the percent relative expression for each gene's time course as described above using the absolute cosine angle distance metric [23]. Heatmaps were visualized using the R package ‘gplots’ [24].

Results

Our understanding of Cryptosporidium infection has remained quite limited due to the parasite's tight association with the host intestinal epithelial cell. We used Real Time (RT)-PCR to describe the transcriptional profile of 3,302 C. parvum genes during a 72 hr in vitro infection time course using the epithelial cell line, HCT8, resulting in the first comprehensive transcriptome for this obligate intracellular parasite. Four thousand thirty-three primer pairs were tested, and 731 were rejected based on irregularity in amplification efficiency or because they amplified host cell cDNA. Four independent in vitro infections were prepared and assessed by RT- PCR. Gene expression ranged from 10−10 to 0.5 transcripts relative to C. parvum 18S rRNA, while the median expression relative to rRNA across the time course did not change (Table 1). Approximately two-thirds of the transcripts identified were expressed at all time points, albeit at different relative transcript levels.

[Figure omitted. See PDF.]

Table 1. Transcriptome Summary Statistics.

https://doi.org/10.1371/journal.pone.0031715.t001

Morphological description of C. parvum developmental stage progression

In agreement with prior data [25], removal of unexcysted oocysts after 2 hr incubation with the host cell culture layer results in relatively homogenous parasite morphology during the first 24 hr of infection (Figure 1). At 2 hr sporozoites have attached and appear as small, dense structures atop the host cell. By 6 hr the sporozoites have expanded into trophozoites, which are medium sized spheres atop the host cell. Thereafter, immature (12 hr) and mature (24 hr) meront development is observed. These appear as larger staining spheres, with merozoites observed within more mature meronts. However, once a full cycle of asexual replication has been completed (between 24–36 hr p.i.) and merozoite reinfection has begun, the relative homogeneity previously observed is diminished as the various parasitic stages progress at different time scales. Morphological description is also complicated by the occurrence of sexual replication within type 2 meronts, which are not distinguishable from asexually derived (type 1) meronts by light microscopy.

[Figure omitted. See PDF.]

Figure 1. Quantification of Parasite Stage Progression.

HCT-8 coverslips were infected and stained with pan-specific antibody to quantify C. parvum infection. >300 parasites were counted on 15–30 frames for each of four replicates for each time point. The average for each of the 3 stages of parasite determined was converted into a percentage of the total parasite detected per time point. zoite  =  sporozoite/merozoite, trophozoite  =  all trophozoites (size of stained parasite was > zoite but <merozoite), meront  =  all parasites staining > trophozoite.

https://doi.org/10.1371/journal.pone.0031715.g001

Gene Expression is not Correlated to Chromosomal Location

All 3,302 genes analyzed were detectable within at least one time point of the 72 hr time course. The transcriptional density across each time point follows a normal distribution, with an increase in detectable genes (non-zero values) occurring from 6, 12, and 24 hr, coinciding with a decrease in C. parvum 18S rRNA Ct values (Table 1). C. parvum 18S rRNA Ct values increase slightly after 36 hr. However, there is no decline in transcript abundance measured as either total number of detectable genes or total mRNA content at 72 hr (Table 1, Figure S2).

Previous work had suggested a transcriptional cascade of gene regulation [26], but our data suggests C. parvum utilizes a complex transcriptional regulatory system. First, a single chromosome exhibiting prominent transcription is not observed at any time point regardless of transcript quantity (Figure 2). Second, by examining the quantity of transcript, tandem genes and divergent genes do not have correlated expression as a norm (data not shown). Third, genes within a cluster (see below) do not map to the same chromosomal location. Complex transcriptional regulation was further verified by examining the tandem cluster of peptidases found on chromosome 3 (Figure 3). C. parvum has a large number of M16 family proteases and their kin. Four genes contain a full protease catalytic domain and 15 degenerate genes appear to be related to M16 enzymes. The most notable are the tandem Insulin Degrading Enzymes (IDEs) found near the telomere of chromosome 3. Using relaxed criteria to identify the metallohydrolase domain, the cluster contained 12 consecutive protease genes, cgd3_4170 through cgd3_4280, spanning a 57 kb region. This grouping of genes appears to be formed by gene duplication. All of these IDE genes were expressed, but at vastly different levels (Figure 3B) and with different expression patterns (Figure 3C). For example, cgd3_4240 was expressed at low levels early in the infection with a sharp peak in transcript levels at 48 hr, while all others showed greater expression throughout the time course. In contrast, cgd3_4260 was constitutively expressed at much higher levels throughout the examined time course. The 12 IDE genes sorted into five different expression pattern clusters, further confirming that C. parvum gene expression is independent of chromosome location, even for co-localized gene families.

[Figure omitted. See PDF.]

Figure 2. Expression by Chromosome.

For each time point, the number of genes with detectable expression (non-zero) was quantified by chromosomal location.

https://doi.org/10.1371/journal.pone.0031715.g002

[Figure omitted. See PDF.]

Figure 3. Chromosome 3 Peptidase Loci are not Co-expressed.

Twelve genes, including 11 peptidase genes are present in tandem on chromosome 3. (A) The protease motifs were identified and color coded based on the conservation level. (B) RT-PCR for each gene was normalized to C. parvum 18S rRNA and the median for four biological replicates for the 7 time point time course is shown. (C) The normalized expression was then averaged by each gene's total expression for all 7 time points.

https://doi.org/10.1371/journal.pone.0031715.g003

Expression profiling of an in vitro 72 hour time course

RT-PCR provides a quantitative evaluation of each gene transcript present, but to delineate patterns of gene expression we focused on the relative change of gene expression rather than the absolute transcript abundance during the 72 hr time period examined. To determine which genes have similar patterns of expression, we normalized to C. parvum 18S rRNA and computed a relative expression for each gene by dividing its expression at each time point by the sum of gene expression for all 7 time points. The expression profile is a vector of 7 relative percentage values, and all genes were clustered to identify shared patterns of expression using the Divisive Analysis (DIANA) algorithm (Figure 4) [21], [22], [27]. Utilizing the distance cutoff of 100, 9 clusters were defined for genes expressed over the 72 hr infection period. Based on putative functional annotations and protein domain similarities, we classified each gene into 1 of 18 categories [20]. Examination of the distribution of functional categories was assessed between clusters (Table 2, Figure 5A). All functional categories except the 2–24 hr clustering ‘Other’ were unequally distributed between clusters at a significance level <0.05 (Figure 5A). In addition, genes annotated as apicomplexan, Cryptosporidium, hypothetical or conserved hypothetical proteins were randomly distributed across all clusters (data not shown).

[Figure omitted. See PDF.]

Figure 4. C. parvum Transcriptome 2–72 hr Relative Expression Cluster Analysis.

The DIANA algorithm was used to cluster each gene's relative expression using the absolute cosign distance metric. For the heatmap each gene is represented in one row and each column represents the percent expression at one time point. (A) The relative expression (%) was determined for each gene's median normalized value over the 7 time point, 72 hr time course (3302 genes). Nine clusters (1–9) are present with the dendrogram distance cutoff of 100 as indicated with the colored bar to the right of the heatmap. (B) For the 72 hr time course, each clusters' average expression (+/− SD) is graphed, with the number of genes found in each cluster within the parentheses.

https://doi.org/10.1371/journal.pone.0031715.g004

[Figure omitted. See PDF.]

Figure 5. Distribution of Genes Across Clusters.

The percentage of each of 13 functional categories was determined across all clusters. Five categories were omitted: Apicomplexan, Cryptosporidium, Oxidative Stress, Hypothetical, and Conserved Hypothetical. Chi-square testing was used to determine any significance in the distribution of the number of genes across clusters using the null hypothesis that the functional classifications would be equally distributed across clusters. Data was using the 2–72 hr (A) and 2–24 hr (B) data sets. Chi-square significance: * = <0.05, ** = <0.01, *** = <0.001, **** = <0.0001.

https://doi.org/10.1371/journal.pone.0031715.g005

[Figure omitted. See PDF.]

Table 2. Functional Categories Distribution across 2-72 hr.

https://doi.org/10.1371/journal.pone.0031715.t002

Cluster 7 represented an expression pattern following the normalizing gene's expression (C. parvum 18S rRNA). Cluster 7 had a high proportion of genes associated with protein transport or modification, but also contained multiple genes for DNA-associated proteins involved in replication, including 5 termination factors and 4 histone-modifying enzymes. Genes with a high expression within this cluster were ethanolamine-phosphatecytidyl transferase and NADP-dependent alcohol dehydrogenase.

Five clusters contained genes with a defined peak in expression at a single time point. Cluster 6 contained the largest number of genes, and expression peaked at 2 hr followed by a lower expression at the remaining time points (Figure 4B). The largest functional category expressed in this cluster was that of RNA-associated genes. Notably, 22 helicases, including Dbp8p-like helicase, Sgs1, Has1p, Rok1p-family helicase, DRS1, Mtr4p-like, and 8 involved in ribosome biogenesis, and 11 spliceosome proteins, including PRP22, splicing factor 3b, FBP21, and U5snRNP had their highest expression at 2 hr. Many transcription-related genes were found in this cluster, including 11 polymerase-associated proteins. Additionally, 19 ribosome and 13 tRNA-associated transcripts were present. We also observed high expression at 2 hr for 6 H+ transporters, 3 mitosis-associated genes (XPMC2, DIM1 and Ard1) and 8 DNA replication-associated genes (including Separin/esp1, SMC1, UVR8-like protein with RCC1 domains, and Cbf5p). Within this cluster thymidine kinase, phosphatidate cytidylyltransferase, initiation factor 3 subunit 8, and 60S ribosomal subunit L7A had the highest magnitude of expression.

As shown in Figure 4B, genes in cluster 9 had a defined peak in expression at 6 hr with a small rise again at 24 hr. The largest functional category of genes with this expression pattern was translation-associated transcripts that included 56 ribosome and 17 tRNA associated genes. Genes with high expression within this cluster included a mucin, Rad23, 40S ribosomal proteins S3 and S19, 60S ribosomal protein L8, and elongation factor 1-gamma.

Clusters 3 and 8 both contained genes with spikes in expression at 12 hr, with cluster 8 having a higher average percentage of expression at 12 hr and little expression thereafter. In contrast, genes in cluster 3 showed almost equal peaks of expression at 12 and 36 hr. For both clusters, the functional class most highly represented are structural genes, in particular cytoskeletal components that include two unique myosins, four actin, three kinesin, and two tubulin gene transcripts. However, cluster 8 contained many genes associated with metabolism, whereas cluster 3 contained transcripts for a variety of predicted transporter genes.

Although there was not a cluster containing genes with peak expression exclusively at 24 hr, genes in cluster 5 showed a bimodal expression pattern that included a small peak at 6 hr and a larger peak at 24 hr. This cluster contained the fewest genes overall, and the largest functional categories represented were structural protein, transcription-related, and DNA-associated genes. The replication factors minichromosome maintenance complex component (MCM) 6 and MCM 7 were found in this cluster along with the gene for Structural Maintenance of Chromosomes (SMC) 2 and the DNA replication-related chl12 gene, which has one of the highest transcript copy numbers within this cluster.

Three clusters (1, 2, and 4) indicate a pattern of increased expression within the last three time points, which correspond to time periods that sexual replication stages are observed. Cluster 2 represented genes with a tightly defined spike in transcript levels at 48 hr, while cluster 1 contained genes in which transcript levels peaked at 72 hr. Cluster 4 contained genes with little or no expression until 36 hr. Afterwards, expression of these genes increased at 48 hr and remained high with a slight drop off at 72 hr. The largest functional category within cluster 4 was DNA-associated genes, including ten DNA repair genes and eleven replication-related genes. Genes encoding enzymes associated with metabolism were the highest represented functional group within clusters 1 and 2, with cluster 2 containing transcripts for three amino acid transporters and the DNA-associated genes ARD1 and topoisomerase II, as well as four COWP genes. Two genes in this cluster, COWP1 and CpCCP2, exhibited high transcript numbers. A spike in expression for DNA-associated genes including DMC1, SMC1 and SMC3 was observed at 72 hr in cluster 1. A similar spike was observed for three cyclin homologues and the remaining four COWP genes, with COWP8 having a highly detected copy number.

Expression profiling of the asexual replication cycle (2–24 hr)

Morphological homogeneity is highest (Figure 1) in the first 24 hr of infection because sexual replication stages do not appear in culture until after 24 hr [28]. We further focused analyses on the first 24 hr to identify genes expressed by asexual developmental stages. The DIANA algorithm distinguished six clusters (divisive coefficient  = 0.98) using the dendrogram distance cutoff of 1 (Figure 6) [24], [27]. One cluster (cluster E) contained genes mirroring 18S rRNA transcript levels and likely represented housekeeping functions for all asexual stages. Interestingly, 30% of the genes in cluster E were hypothetical proteins (Table 3), suggesting we know little of even the basic biochemistry for this pathogen. Four clusters represented genes with a spike in expression at each of the four time points, while cluster C contained genes with a bimodal expression pattern at 6 and 24 hr. Across the 6 clusters, almost 50% of the genes were assigned between cluster D (6 hr spiked expression pattern) and cluster B (12 hr spiked expression), while cluster A (24 hr maximum expression) had the fewest genes. Five genes were not included in this cluster analysis because they had no detectable expression within the first 24 hr (Table 4), including cgd4_1110, which is a hypothetical gene with a predicted AP2 domain at both the N- and C-termini.

[Figure omitted. See PDF.]

Figure 6. C. parvum Transcriptome 2–24 hr Relative Expression Cluster Analysis.

The DIANA algorithm was used to cluster each gene's relative expression using the absolute cosign distance metric. For the heatmap each gene is represented in one row and each column represents the percent expression at one time point. (A)The relative expression (%) was determined for each gene's median normalized value over the first 24 hr of the time course, which includes 4 time points (3297 genes). Six clusters (A–F) of genes are present with a dendrogram distance cutoff of 1, as depicted with the color bar to the right of the heatmap. (B) The average gene expression per cluster (+/− SD) using the first 4 time points (24 hr) is plotted with the number of genes found in each cluster within the ().

https://doi.org/10.1371/journal.pone.0031715.g006

[Figure omitted. See PDF.]

Table 3. Functional Categories Distribution across 2-24 hr.

https://doi.org/10.1371/journal.pone.0031715.t003

[Figure omitted. See PDF.]

Table 4. Genes with no Detectable Expression in the first 24 hr.

https://doi.org/10.1371/journal.pone.0031715.t004

This focused analysis revealed additional genes that exhibited high expression at 2 hr within the first 24 hr, but showed increased expression at later time points 36, 48 or 72 hr when analyzing the entire time course. Using this approach we found 96 genes differing within the 2 hr clusters, of which 48% are annotated as hypothetical proteins. At 2 hr post-infection, transcripts may represent newly transcribed genes or may have been carried forward from sporozoite packaging within the oocyst. Further, disregarding expression data collected at 36–72 hr and focusing on the first round of asexual development (2–24 hr) resulted in differing gene expression patterns compared to analysis with the entire time course (2–72 hr). With this asexual specific analysis (Figure 5B), the 2 hr time point now showed the highest relative expression of genes CpCCp1, COWP1, COWP2, glutathioneperoxidase, triosephosphate isomerase, Cyclin A, PRP2, and RNA polymerase II B6 and B10 subunits.

Similarly, at 6 hr we found a 153 gene difference compared to the 72 hr clustering, with only 35% being hypothetical. Transcripts for two mucins, the DNA polymerase epsilon subunit, MUS81 endonuclease, CDC45, a chromatin protein similar to BDF1, 13 metabolism related genes, HSP40, RAB11, PRP5, and a FLX1-like transporter all peaked at 6 hr, which corresponds to the time point when a high percentage of trophozoites were present.

An additional 223 genes show an expression peak at 12 hr, half of which were annotated as hypothetical proteins. Three Cryptosporidium specific genes, TRAPC1, TSP1, and a mucin, were all highly expressed. Additionally, DMC1, topoisomerase VIA, MCM10, endonulcease III, histone H4, and 14 additional metabolism related genes had increased expression at this time point.

At 24 hr we observed 44 genes varying from the previous 72 hr cluster scheme, 39% of which were annotated as hypothetical. Genes with a predicted function included TSP6, MCM6 and 7, SMC2 and 4, RNA polymerase II CTD/NL1 interacting protein, and translation initiation factor if-2B beta.

The largest difference found was in the genes within this analysis that have a bimodal expression at 6 and 24 hr. Four hundred-two genes found within cluster C are not found in cluster 5, with 35% being hypothetical. This group of genes, found in cluster C (Figure 5), included eight Apicomplexan genes, including lipin2, ORC subunit 5, and cgd7_4100 (membrane protein), and five Cryptosporidium related genes: TSP10, TSP4, mucin and a COWP. The largest category represented was that of DNA-associated genes, such as repC1, DNA polymerase theta, delta 1 and 2, MCM2 and 3, bub1, ORC1 and 2, PDNA, HDAC2, topoisomerase I, and DNAse 1. In addition, cluster C included CHC1 clathrin heavy chain, as well as the Adaptin AP complex subunits alpha, beta, mu and gamma.

Spikes in Gene Expression

Seventy-four percent (2,424) of the genes detected during the first 24 hr of infection have detectable expression at all four time points (Figure 7). Of the genes with variable detection, 42 were identified at only one time point: one at 2 hr, two at 12 hr, and thirty-nine at 24 hr post-infection (Figure 7). This patterning further suggests that the number of transcripts detected is dependent on the parasite stage present and changes as the parasite develops. To ascertain large shifts in transcription across the entire time course used, we identified transcripts with >80% of its total expression within a single time point over the 72 hr time course, indicative of at least a 4-fold induction of transcription. This analysis resulted in 107 time-point specific genes, 72 of which show >80% of their total expression at either the 48 or 72 hr time points only (Table S3). The residual 20% expression for these 72 genes fell primarily in either the 48 or 72 hr time point, suggesting they are expressed by sexual developmental stages only.

[Figure omitted. See PDF.]

Figure 7. Unique Asexual Stage Gene Expression.

Gene transcription fluctuates during the infection time course. Most genes are detected at all 4 time points within the first 24 hours of infection (asexual replication). Gene lists were constructed containing the gene IDs for all genes with detectable transcription (non-zero) at 2, 6, 12, and 24 hr time points. These lists were compared to show the number of genes with overlapping expression at all combination of the four time points used. Blue = 2 hr, yellow = 6 hr, green = 12 hr, and red = 24 hr [72].

https://doi.org/10.1371/journal.pone.0031715.g007

An unexpectedly large number (59%) of the 107 genes with >80% expression at one time point were annotated as hypothetical genes (Table S3). The highest number of remaining genes were categorized as Cryptosporidium-specific, interacting, and phosphorylation-related genes (Table S3). Genes having a functional annotation included the Rab GDP dissociation inhibitor (cgd7_4740) spiking at 36 hr and phosphoglucomutase and two myosins, which spiked at 12 hr. At 48 hr we observed a spike in expression of seven Cryptosporidium related genes, including five mucins, CpCCP3, and COWP6. Also of note at 48 hr were genes encoding translation initiation factor 3 subunit 10 and a predicted oocyst wall protein (CpPOWP) 4 [29]. Genes represented at 72 hr included COWP7, CpPOWP3 [29], DMC1, pyruvate decarboxylase, and cyclin A.

Discussion

We applied RT-PCR and cluster analyses to describe changes in the C. parvum transcriptome during the first 72 hr in vitro infection of intestinal epithelial cells. The time points assayed (2, 6, 12, 24, 36, 48, and 72 hr) were predicted to be integral transitional states during in vitro parasite development based on our morphological assessment [25] (Figure 1) and other reports [30]. Transcripts were measured using four independent culture replicates. We detected transcripts for all 3,302 genes tested (Table 1). This detection rate was surprising considering that the intensely studied genome of Sachromyces cerevisiae contained 12.3% “dubious” ORFs with no transcripts detected [31]. When transcript abundance for each gene was standardized as a percentage of its total transcript abundance across all time points, we found 107 genes with spikes in expression of >80% at a single time point.

The dataset produced 9 clusters based on relative level of expression in which the 18 functional categories were distributed unequally. The representation of functional categories of genes at each time point suggested a choreographed developmental program for this parasite during infection. At 2 hr post infection, sporozoites expressed proportionately more RNA-associated, transcription and transporter genes. At 6 hr trophozoites expressed more translation-related genes, at 12 hr meronts exhibited more transcripts for structural protein-related genes, transporters, and enzymes required for metabolism, while the 24 hr meronts had an over-representation of DNA-associated genes. At 36 hr merozoites have reinfected new epithelial cells for another round of replication, and DNA-associated genes were again overrepresented. However, by 48 and 72 hr post-infection, metabolic enzymes were the overrepresented functional category. Comprehensively, our transcriptome suggests a cascade of gene expression consistent with unique biologies for each developmental stage.

Sporozoites: preparation for parasitism and transcription

At 2 hr, the sporozoite has attached to the host cell and appears to be ramping up transcription in preparation for cellular parasitization [32]. In our analysis, transcripts for 2,532 genes had detectable expression at 2 hr, representing the smallest population out of the 7 time points examined, as might be predicted based on the assumption that sporozoites are quiescent within the oocyst. Previously published EST data [9]–[12], although not quantitative, indicates <1,000 non-repetitive transcripts which suggests that host attachment triggers a robust gene expression profile as the parasite begins to capture nutrients from the host cell. The freshly excysted sporozoite is a transient state that is primed for translation, based on the high amount of ribosomal proteins found in both EST and proteomic data [9]–[11]. Excystation itself increases expression of ribosomal protein genes, yet our data showed that additional translational machinery is not highly transcribed until after 12 hr [33], even though many transcripts for previously undetected genes appear at 2 hr. This illustrates that Cryptosporidium sporozoites are packaged with all components needed to attach, commence parasitophorous vacuole formation, and initiate transcription following attachment to the host cell. This rapid increase in transcription following attachment is similar to the ring stage of Plasmodium and tachyzoite replication in Toxoplasma [34]–[36]. These events have not been captured in prior Cryptosporidium EST and proteome analyses and represent the power in examining a global transcriptome approach.

The majority (96%) of the genes detected in the first 24 hr were expressed at all four time points (Figure 7), and 16% had elevated expression at 2 hr. However, transcripts for cgd8_730, a possible Apicomplexan protein, were only detected at 2 hr within the first 24 hr. The P. falciparum homologue for cgd8_730, PFD0910w, shows highest expression during gametogenesis which would suggest that cgd8_730 in sporozoites is a remnant transcript that was expressed during sexual development [37]. Toward this end, we observed very low cgd8_730 transcripts levels at 2 hr, but transcripts increased >1,000-fold at 48 and 72 hr, consistent with expression during gametogenesis. The 2 hr transcriptome included many transporters, DNA-associated proteins, and transcription-related genes, which peaked at 2 hr and were turned off or greatly down-regulated in proceeding time points. This is surprising since many of these genes might be expected to show constitutive expression over the course of an infection. The possibility that the observed down-regulation is an artifact due to normalization must be considered. However, normalization using parasite 18S rRNA gene quantity had been previously established under these culture conditions to correlate with parasite density [7], [8]. Even so, we observed that ribosomal components decreased in expression at later time points (72 hr) despite a lack of corresponding decrease in total number of C. parvum transcripts or parasite numbers detected (Table 1). Because individual parasite stages cannot be readily purified away from the parasitized host cell, the measured RNA concentration used to initiate cDNA synthesis is largely derived from the host cell RNA and thus this phenomena may also be related to host cell growth or demise during infection.

The Cryptosporidium genome lacks genes for energy synthesis beyond glycolysis. Therefore, expression of transporters was expected to initiate early following attachment. Fewer than half of the annotated transporters had been detected previously in freshly excysted sporozoites [4], [6], [9], [12]. Our analyses found that almost 20% of the genes annotated as transporters, including most of the sugar and nucleotide-like transporters, had high expression at 2 hr. Additionally, two putative amino acid transporters (cgd3_470 and cgd5_1630) had 50% of their relative expression at this time point. This is consistent with the parasite needing quick access to the host cell contents for nutrients necessary for developmental stage progression.

The Cryptosporidium specific genes TSP8 (cgd6_780) and CpCCP1(cgd7_1730) were expressed at high levels at 2 hr. TSP8 proteins are abundant in sporozoites and type 1 merozoites [38]. Our data indicated that 79% of TSP8's relative expression was at 2 hr, suggesting the protein is made during the first round of asexual development following sporozoite invasion and either was not required for later asexual divisions or represents a very stable protein. In apparent contrast to prior work showing that the CpCCP1 protein was abundant in sporozoites and was embedded throughout the circumference of the parasitophorous vacuole during intracellular development [39], we found CpCCP1 had >90% relative transcript expression in the combined 48 and 72 hr time points. This apparent conflict can be explained if the CpCCP1 protein, as opposed to new gene transcription, was packaged in oocysts and carried forward rather than induced de novo gene expression following excystation. It may also suggest that CpCCP1 is more of a generic protein important for parasitophorous vacuole integrity throughout the entirety of the infection rather than being used only by sporozoites during initial host cell invasion.

Trophozoites: Protein synthesis and degradation

The largest number of genes had peak expression at 6 hr. Trophozoites are the major morphological stage present at this time and our data showed an expression increase of genes required for protein translation. In addition, we observed increased transcription for genes encoding proteins required for folding, modification and transport of the newly translated polypeptides. This was supported with the detection of 15 chaperone proteins, including many heat shock proteins (HSPs). Additionally, the trophozoite is transcribing almost the entire complement of proteasome genes which could be used to recycle either proteins no longer needed by the sporozoite or newly acquired host proteins to facilitate new protein synthesis (Figure S3).

In conjunction with parasite preparation for new protein synthesis, trophozoites are enlarging in preparation for mitosis and division into merozoites through schizogony [30], [40]. Nutrient availability, protein synthesis, and ribosome biogenesis are major factors needed for mitosis in eukaryotic organisms [41]–[43]. Our analyses showed that all but three of the parasite ribosomal and tRNA synthesis genes peaked in expression at 6 hr (Figure S3). The exceptions, 60S ribosome L7A (cgd3_760), tRNA related adenyl cyclase associated protein (cgd5_440), and L35 (cgd3_830), had peaks at 12 hr or bimodally at 12 and 36 hr, respectively. Schizogony involves multiple rounds of DNA division followed by daughter cell budding. Our data clearly supported that the trophozoite stage was preparing for chromosomal duplication and separation with the expression of transcripts for CDC45, all replication licensing factors (MCM), replication factors C2, C5 and repA1, replication protein A1 large subunit, origin recognition complex (ORC)1 and ORC2, and the spindle checkpoint gene bub1. Plasmodium lacks orthologs for CDC45 and Bub1, but ORC1 transcription is highest within the schizont stage, with protein detected in the trophozoite [44]–[46]. ORC proteins have not been detected in the Cryptosporidium sporozoite proteomes, but we observe increased expression of all ORCs at 24 hr (Figure S3). This suggests the merozoite may be prepackaging the transcripts necessary to proceed with DNA-replication following attachment.

Immature (12 hr) Meronts: Dividing nuclei

It has been previously reported in other Apicomplexa that nuclear division occurs within the immature type 1 meront and is followed by cytokinesis, resulting in the formation of the mature type 1 meront [47]. The 12 hr transcriptome supports a similar division process for C. parvum. Transcripts for the H2B, H3 and H4 histones were highly expressed at 12 hr. Acetyltransferase GCN5, which activates transcription of genes required for DNA replication and movement from G1/S based on nutrient availability, MCM10, which regulates DNA replication initiation in fission and budding yeast [48], and PMS1, which functions in DNA repair, also showed maximum expression at 12 hr [49].

The immature C. parvum meront transcriptome also included genes for structural components found to precede zoite formation in other Apicomplexa, such as the apical complex proteins that include components of the microneme and rhoptry, as well as the formation of the Inner Membrane Complex [50]-[52]. These complexes are composed of a variety of membrane skeletal proteins, such as actins, myosins, tubulins, and kinesins. Within this 12 hr cluster were five Cryptosporidium-specific genes, TSP3, TSP11, TSP12, TRAP-C1, and Cp15/60, confirming previous transcription data [14], [53]. TSP3, TSP11, and TRAP-C1 have been documented as microneme-associated proteins, further suggesting the parasite is preparing for the future separation of daughter cells [9]. Cp15/60 (cgd7_4540) and cgd3_1570 are previously reported sporozoite antigens, but their high transcript levels at 12 hr in our culture system suggested they are also components of merozoites.

Metabolic enzymes required for energy conversion and potential storage of sugars were expressed at their highest level at 12 hr in the developing merozoites and may be used to complete development or store energy for the traversal to the next host cell. Overall, 8 of the 9 genes used in glycolysis were present in our data set (Figure 8). Three of these genes were constitutively expressed (cgd6_3800 [hexokinase], cgd2_2130 [pyrophosphate-dependent phosphofructokinase], and cgd1_2040 [pyruvate kinase]), whereas cgd2_3200 [glucose-6-phosphate isomerase] spiked in expression at 6 hr. Transcripts for the remaining enzymes cgd1_3020 [fructose-1 disphosphate aldolase], cgd6_3790 [glyceraldehyde 3-phosphate dehydrogenase], cgd7_4270 [phosphoglycerate mutase], and cgd5_1960 [enolase] all spiked at 12 hr. It is interesting to note that these same four proteins are found in the sporozoite proteome [10], suggesting that the type I meront/merozoite also has a high need for energy production. Non-constitutive expression of these enzymes is not surprising since Toxoplasma also demonstrates stage-specific expression of enolase and pyruvate kinase [54]–[56]. In addition, cgd3_3100, a putative sugar transporter, was highly expressed at this time point providing further evidence for an increase in metabolism in type I meronts. Interestingly, most of the host cell glycolytic pathway genes have been analyzed during a very similar infection time course [57]. Combining data generated in this study indicates that while the parasite is increasing cgd1_3020, cgd6_3790, and cgd5_1960 expression at 12 hr, host cell expression of orthologs is decreasing [58] suggesting that the host cell is decreasing its own glucose usage in response to parasite invasion.

[Figure omitted. See PDF.]

Figure 8. Glycolysis Related Transcripts.

The relative expression using 2-72 hr data is graphed for the annotated enzymes involved in glycolysis.

https://doi.org/10.1371/journal.pone.0031715.g008

Mature Meronts/Merozoites: cell cycle completion and DNA repair

The 24 hr time point contained the highest fraction of mature type I meronts with 6–8 merozoites that have completed cytokineses (Figure 1). Although both sporozoites and merozoites infect epithelial cells within the same host, their different expression profiles and distribution of functional categories between 2 hr and 24 hr indicated that sporozoites are biochemically distinct from merozoites. Most notably DNA-associated genes involved in replication and mitosis were specifically elevated at 24 hr. Surprisingly, transcripts for two core components for the condesin complex, SMC2 and SMC4 [59], [60], topoisomerase II, DNA ligase I, and RAD45 were elevated after cytokinesis. Since merozoite cytokinesis is largely complete by 24 hr, this suggests that either these proteins are produced and packaged for subsequent infection by merozoites or the transcripts have a very long half-life. The pre-packaging of proteins needed in the next round of DNA synthesis within merozoites was illustrated further by the high expression of CDC6, MCM6, and MCM7 at 24 hr. In eukaryotes, the MCM2-7 complex assembles on ORC at early G1 to form a prereplicative complex [61]. Activation of the replicative origin by CDC6 prevents cells from entering S phase when environmental conditions are unfavorable. Here, Cryptosporidium merozoites appear poised to initiate the next cell cycle upon reinfection. Other MCMs were more highly expressed at earlier time points. Taken together, the transcriptome is consistent with a mechanism in which the sporozoites infect with replication complexes that are missing a few subunits, and that DNA replication ensues only after translation of the missing factors. Alternatively, the asynchrony of MCM synthesis could also indicate that Cryptosporidium, like Plasmodium, may use a different helicase formation that is developmentally regulated during asexual development, similar to the expression and use of the MCM2/6/7 complex in Plasmodium [62].

Sexual stage-associated transcripts

The previous discussion was based on the clustering of the first 24 hr because of the relatively synchronous developmental stage progression permits us to look specifically at the transcripts required for one round of asexual replication. When taking the relative expression utilizing the entire 72 hr time course, several replication-associated transcripts spiked in expression at later time points, appearing to mute expression within the first 24 hr. In later time points, additional sexual stages are developing, and certain expression patterns may be diluted by simultaneous and asynchronous progression of both asexual and sexual developmental programs. For example, gametogenesis occurs after 24 hr and requires rapid DNA replication. This would explain the large increase in DNA replication genes within the later time points (Figure S3). Additionally, it is highly plausible that like other Apicomplexans, Cryptosporidium sexual replication may not use schizogony, thus changing the transcriptional requirements for sexual development [50].

At 36 hr reinfection of epithelial cells by merozoites has occurred, which coincides with an overrepresentation of DNA-associated genes. Genes in cluster 3 had similar peaks in expression at 12 and 36 hr, whereas cluster 8 genes are much more highly expressed at 12 hr. The cyclical expression of some genes between 12 hr and 36 hr (cluster 3) suggests a reinfection of asexually replicating merozoites had a similar gene expression pattern to that of the first round of asexual replication. However, having a separate subset of genes that does not follow this cyclical patterning (cluster 8) also suggests that there is a difference in the cultures at these time points. We noted that only cgd7_4740, Rab GDP disassociation inhibitor, has a high (84%) relative expression at 36 hr, suggesting it is specific to a new stage that is present at this time. Interestingly, cgd4_3110, a putative nucleotide-sugar transporter and cgd2_800, a Major facilitator super family transporter, also had the majority of their expression at 36 hr, which may indicate a change in the nutritional requirements of this new developing parasite stage.

The culture at 48 hr has increased morphological complexity, with over half of the organisms appearing to be of meront size or greater. Despite this developmental stage heterogeneity, the appearance of genes that peaked at 48 hr (clusters 2 and 4) indicated the presence of new developmental stages. Metabolic enzymes and transporters were the most overrepresented functional categories expressed at this time point. Metabolic enzymes such as fatty acid synthetase, S-adenosylhomocysteinase, and pyruvate:ferredoxin oxidoreductase showed increased expression as the time course ensues. This is surprising considering that the streamlined metabolic pathways in Cryptosporidium might predict a more constitutive expression pattern rather than the striking peak observed here. Increased expression of these enzymes could reflect a shift in the metabolic need of the sexual stages developing at this time or may be due to nutrient depletion of the host cells forcing the parasite to shift its energy sources. Relative expression of the amino acid transporters cgd3_2730, cgd3_2050, and cgd3_4350 was highest at 48 hr and further emphasizes a shift in metabolic requirements. A large number of genes at this time point appeared to encode extracellular proteins, including 5 mucin-like genes, one of which (cgd3_720), showed no detectable expression until 48 hr. Three protein kinases and one phosphatase were highly expressed and may be key regulators in the life cycle progression of this parasite.

Clearly the 72 hr transcriptome indicates significant changes within the culture, and 37 of the 107 putative developmental stage-specific genes (>80% expression at one time point) peaked at 72 hr. Metabolic enzymes are again the overrepresented category at 72 hr, with 82% of pyruvate decarboxylase total transcripts observed at this time. Relative expression of cgd7_4800 (amino acid transporter) and cgd2_590 (nucleotide-sugar transporter) also peak at this time point, which is consistent with a change in the metabolic requirements of the parasite. There are eighteen AP2 domain containing genes in our dataset out of the nineteen predicted in Cryptosporidium [63]. The hypothetical protein, cgd4_1110, an AP2 domain containing protein, had no detectable expression in the first 24 hr, but has over 95% of its expression between 48 and 72 hr [63]–[65]. In addition, another AP2 domain containing gene (cgd2_3490) had the majority of its expression during the 48 and 72 hr time points, whereas the AP2 domain containing proteins, cgd4_2950 and cgd8_3130 had >60% expression at 12 hr. This expression patterning suggests that like the other Apicomplexa, stage-specific utilization of AP2 domain containing transcription factors by C. parvum controls developmental stage-associated gene expression.

Cryptosporidium sporozoites infect as diploid organisms and meiosis is required to proceed through zygote formation. Genes known to be integral for meiosis in many organisms were highly expressed specifically at 48 and 72 hr in our culture. Dmc1, CDK (cgd5_2510) and Cyclin A (cgd3_4050) all had their highest expression at 72 hr, while Hop1 peaked at 48 hr. ORC4 had over 50% of its expression at 48 hr alone, while ORC2 and ORC5 both demonstrate cyclic expression that peaked at 6, 24, and 48 hr. When analyzing all 7 time points, MCM 2, 3, and 5–7 were more constitutively expressed, although MCM 2, 3 and 10 had higher relative expression at 48 hr. We found ORC1 and CDC6 expression to follow similar patterns, until 72 hr where CDC6 remained high while ORC1 expression plummeted (Figure S3). This suggests that ORC1, CDC6, and perhaps OCR4 may form a replication initiation complex during the sexual stages present at later time points, whereas ORC2 and ORC5 have more universal roles in DNA replication and differs from the model in which ORC1-ORC6 and CDC6 function as a stoichiometric complex for DNA origin initiation [66]. However, recent mapping of ORC binding has identified many sites that are not consistent with origin of replications [67], [68]. It is likely that due to the streamlined genome of Cryptosporidium many proteins possess multiple functions, in this case perhaps an additional chromatin-associated role for ORC1.

Conclusions

We observed transcripts for all 3,200+ genes assessed within at least one time point over 72 hr in vitro infection of HCT-8 cells. This 72 hr time course covered the known length of time parasite can be propagated in vitro without a large decrease in parasite numbers detected through 18S rRNA Ct values. Following the removal of unexcysted oocysts at 2 hr incubation, the first 24 hr of the time course remained highly homogenous despite clear life cycle progression. Thus, transcripts specific to the predominant parasitic stage at specific time points were discovered and transcriptome differences between sporozoites and merozoites were discernible. Over 75% of the genome had detectable transcription throughout the 72 hr time course, yet the identification of genes with large spikes in expression at a single time point was particularly intriguing as they may encode developmental stage-specific proteins [53], [69]–[71]. It is even more interesting that the vast majority of our putative developmental stage-specific genes are of unknown function and lack any recognizable protein motif. These putative stage specific genes may reveal biology that is unique to a specific developmental stage and may be useful for purifying individual stages from the heterogeneous cultures or as chemotherapeutic targets.

Supporting Information

[Figure omitted. See PDF.]

Figure S1.

cDNA Synthesis Control. Lack of genomic contamination was confirmed for the cDNA time courses with RT-PCR of replicate reactions of samples with (black) and without (red) reverse transcription. All the reactions without reverse transcriptase had a Ct value >35 with amplicon melt temperatures ranging from 73–75°C. Due to the nature of ribosomal primers, anything >35 with a melt temperature <83–85°C is deemed as no product detected.

https://doi.org/10.1371/journal.pone.0031715.s001

(TIFF)

Figure S2.

Categorization of Gene Expression. RT-PCR was run on 3302 genes in the C. parvum genome and normalized to C. parvum 18S rRNA expression. Each gene's median expression from four biological replicates was categorized for each time point as being: zero (not detected), high (log10 normalized value >2 standard deviations (SD) from mean), moderate (log10 normalized value within 2 SD from mean), or low (log10 normalized value <2 SD from mean) expression at each time point.

https://doi.org/10.1371/journal.pone.0031715.s002

(TIF)

Figure S3.

Relative Expression for DNA Replication Associated Proteins, Proteasome Subunits, Ribosome related and tRNA synthesis. The relative expression for 2–72 hr data is graphed for the (A) ORC, (B) MCM and (C) Meiosis related genes, (D) tRNA synthesis, (E) ribosomal proteins, (F) alpha and (G) beta proteasome subunits.

https://doi.org/10.1371/journal.pone.0031715.s003

(TIF)

Table S1.

C. parvum 18S rRNA-normalized transcript abundance for each of the four infection time courses.

https://doi.org/10.1371/journal.pone.0031715.s004

(PDF)

Table S2.

Primers used for C. parvum RT-PCR.

https://doi.org/10.1371/journal.pone.0031715.s005

(PDF)

Table S3.

C. parvum genes with >80% of total expression at a single time point.

https://doi.org/10.1371/journal.pone.0031715.s006

(PDF)

Acknowledgments

The authors thank Dr. Cavan Reilly for his statistical assistance and Dr. Michael Murtaugh for critical reading of the manuscript. We are also grateful to Julie Hendrickson, Xinxin “Shirley” Yao, and Karen Jasper for their technical assistance.

Author Contributions

Conceived and designed the experiments: MSA SE CAL. Performed the experiments: MJM CAL. Analyzed the data: MJM SE MSR. Contributed reagents/materials/analysis tools: MSA MSR. Wrote the paper: MJM MSR.

Citation: Mauzy MJ, Enomoto S, Lancto CA, Abrahamsen MS, Rutherford MS (2012) The Cryptosporidium Parvum Transcriptome during In Vitro Development. PLoS ONE7(3): e31715. https://doi.org/10.1371/journal.pone.0031715

References

1. Frost FJ, Muller TB, Calderon RL, Craun GF (2004) Analysis of serological responses to Cryptosporidium antigen among NHANES III participants. Ann Epidemiol 14: 473–478.FJ FrostTB MullerRL CalderonGF Craun2004Analysis of serological responses to Cryptosporidium antigen among NHANES III participants.Ann Epidemiol14473478

2. Centers for Disease Control and Prevention (2010) Cryptosporidiosis surveillance-united states, 2006–2008 and giardiasis surveillance- United States, 2006–2008. MMWR 59(No. SS-6): 1–12.Centers for Disease Control and Prevention2010Cryptosporidiosis surveillance-united states, 2006–2008 and giardiasis surveillance- United States, 2006–2008.MMWR59No. SS-6112

3. Chappell CL, Okhuysen PC, Langer-Curry R, Widmer G, Akiyoshi DE, et al. (2006) Cryptosporidium hominis: Experimental challenge of healthy adults. Am J Trop Med Hyg 75: 851–857.CL ChappellPC OkhuysenR. Langer-CurryG. WidmerDE Akiyoshi2006Cryptosporidium hominis: Experimental challenge of healthy adults.Am J Trop Med Hyg75851857

4. Thompson RC, Olson ME, Zhu G, Enomoto S, Abrahamsen MS, et al. (2005) Cryptosporidium and cryptosporidiosis. Adv Parasitol 59: 77–158.RC ThompsonME OlsonG. ZhuS. EnomotoMS Abrahamsen2005Cryptosporidium and cryptosporidiosis.Adv Parasitol5977158

5. Xu P, Widmer G, Wang Y, Ozaki LS, Alves JM, et al. (2004) The genome of Cryptosporidium hominis. Nature 431: 1107–1112.P. XuG. WidmerY. WangLS OzakiJM Alves2004The genome of Cryptosporidium hominis.Nature43111071112

6. Abrahamsen MS, Templeton TJ, Enomoto S, Abrahante JE, Zhu G, et al. (2004) Complete genome sequence of the apicomplexan, Cryptosporidium parvum. Science 304: 441–445.MS AbrahamsenTJ TempletonS. EnomotoJE AbrahanteG. Zhu2004Complete genome sequence of the apicomplexan, Cryptosporidium parvum.Science304441445

7. Abrahamsen MS, Schroeder AA (1999) Characterization of intracellular Cryptosporidium parvum gene expression. Mol Biochem Parasitol 104: 141–146.MS AbrahamsenAA Schroeder1999Characterization of intracellular Cryptosporidium parvum gene expression.Mol Biochem Parasitol104141146

8. Cai X, Woods KM, Upton SJ, Zhu G (2005) Application of quantitative real-time reverse transcription-PCR in assessing drug efficacy against the intracellular pathogen Cryptosporidium parvum in vitro. Antimicrob Agents Chemother 49: 4437–4442.X. CaiKM WoodsSJ UptonG. Zhu2005Application of quantitative real-time reverse transcription-PCR in assessing drug efficacy against the intracellular pathogen Cryptosporidium parvum in vitro.Antimicrob Agents Chemother4944374442

9. Sanderson SJ, Xia D, Prieto H, Yates J, Heiges M, et al. (2008) Determining the protein repertoire of Cryptosporidium parvum sporozoites. Proteomics 8: 1398–1414.SJ SandersonD. XiaH. PrietoJ. YatesM. Heiges2008Determining the protein repertoire of Cryptosporidium parvum sporozoites.Proteomics813981414

10. Snelling WJ, Lin Q, Moore JE, Millar BC, Tosini F, et al. (2007) Proteomics analysis and protein expression during sporozoite excystation of Cryptosporidium parvum (coccidia, apicomplexa). Mol Cell Proteomics 6: 346–355.WJ SnellingQ. LinJE MooreBC MillarF. Tosini2007Proteomics analysis and protein expression during sporozoite excystation of Cryptosporidium parvum (coccidia, apicomplexa).Mol Cell Proteomics6346355

11. Strong WB, Nelson RG (2000) Gene discovery in Cryptosporidium parvum: Expressed sequence tags and genome survey sequences. Contrib Microbiol 6: 92–115.WB StrongRG Nelson2000Gene discovery in Cryptosporidium parvum: Expressed sequence tags and genome survey sequences.Contrib Microbiol692115

12. Yamagishi J, Wakaguri H, Sugano S, Kawano S, Fujisaki K, et al. (2011) Construction and analysis of full-length cDNA library of Cryptosporidium parvum. Parasitol Int 60: 199–202.J. YamagishiH. WakaguriS. SuganoS. KawanoK. Fujisaki2011Construction and analysis of full-length cDNA library of Cryptosporidium parvum.Parasitol Int60199202

13. Upton SJ, Tilley M, Brillhart DB (1995) Effects of select medium supplements on in vitro development of Cryptosporidium parvum in HCT-8 cells. J Clin Microbiol 33: 371–375.SJ UptonM. TilleyDB Brillhart1995Effects of select medium supplements on in vitro development of Cryptosporidium parvum in HCT-8 cells.J Clin Microbiol33371375

14. Deng M, Templeton TJ, London NR, Bauer C, Schroeder AA, et al. (2002) Cryptosporidium parvum genes containing thrombospondin type 1 domains. Infect Immun 70: 6987–6995.M. DengTJ TempletonNR LondonC. BauerAA Schroeder2002Cryptosporidium parvum genes containing thrombospondin type 1 domains.Infect Immun7069876995

15. Kubista M, Andrade JM, Bengtsson M, Forootan A, Jonak J, et al. (2006) The real-time polymerase chain reaction. Mol Aspects Med 27: 95–125.M. KubistaJM AndradeM. BengtssonA. ForootanJ. Jonak2006The real-time polymerase chain reaction.Mol Aspects Med2795125

16. Stahlberg A, Hakansson J, Xian X, Semb H, Kubista M (2004) Properties of the reverse transcription reaction in mRNA quantification. Clin Chem 50: 509–515.A. StahlbergJ. HakanssonX. XianH. SembM. Kubista2004Properties of the reverse transcription reaction in mRNA quantification.Clin Chem50509515

17. Wakaguri H, Suzuki Y, Sasaki M, Sugano S, Watanabe J (2009) Inconsistencies of genome annotations in apicomplexan parasites revealed by 5′-end-one-pass and full-length sequences of oligo-capped cDNAs. BMC Genomics 10: 312.H. WakaguriY. SuzukiM. SasakiS. SuganoJ. Watanabe2009Inconsistencies of genome annotations in apicomplexan parasites revealed by 5′-end-one-pass and full-length sequences of oligo-capped cDNAs.BMC Genomics10312

18. Schroeder AA, Brown AM, Abrahamsen MS (1998) Identification and cloning of a developmentally regulated Cryptosporidium parvum gene by differential mRNA display PCR. Gene 216: 327–334.AA SchroederAM BrownMS Abrahamsen1998Identification and cloning of a developmentally regulated Cryptosporidium parvum gene by differential mRNA display PCR.Gene216327334

19. Fontaine M, Guillot E (2003) Study of 18S rRNA and rDNA stability by real-time RT-PCR in heat-inactivated Cryptosporidium parvum oocysts. FEMS Microbiol Lett 226: 237–243.M. FontaineE. Guillot2003Study of 18S rRNA and rDNA stability by real-time RT-PCR in heat-inactivated Cryptosporidium parvum oocysts.FEMS Microbiol Lett226237243

20. Jensen BC, Sivam D, Kifer CT, Myler PJ, Parsons M (2009) Widespread variation in transcript abundance within and across developmental stages of Trypanosoma brucei. BMC Genomics 10: 482.BC JensenD. SivamCT KiferPJ MylerM. Parsons2009Widespread variation in transcript abundance within and across developmental stages of Trypanosoma brucei.BMC Genomics10482

21. R Development Core Team (2009) R: A language and environment for statistical computing. 2.10.1. R Development Core Team2009R: A language and environment for statistical computing. 2.10.1.Available: http://www.R-project.org. Accessed July 1, 2011. Available: http://www.R-project.org. Accessed July 1, 2011.

22. Maechler M, Rousseeuw P, Struyf A, Hubert M, Hornick K (2011) Cluster analysis basics and extensions. R package version 1.14.1: M. MaechlerP. RousseeuwA. StruyfM. HubertK. Hornick2011Cluster analysis basics and extensions.R package version1.14.1

23. Eisen MB, Spellman PT, Brown PO, Botstein D (1998) Cluster analysis and display of genome-wide expression patterns. Proc Natl Acad Sci U S A 95: 14863–14868.MB EisenPT SpellmanPO BrownD. Botstein1998Cluster analysis and display of genome-wide expression patterns.Proc Natl Acad Sci U S A951486314868

24. Warnes GR (2009) Includes R source code and/or documentation contributed by (in alphabetical order):. In: Bolker B, Bonebakker L, Gentleman R, Huber W, Liaw A, Lumley T, Maechler M, Magnusson A, Moeller S, Schwartz M, Venables W, editors. GR Warnes2009Includes R source code and/or documentation contributed by (in alphabetical order):B. BolkerL. BonebakkerR. GentlemanW. HuberA. LiawT. LumleyM. MaechlerA. MagnussonS. MoellerM. SchwartzW. Venables(2009) Gplots: Various R programming tools for plotting data. 2.7.4. Available: http://CRAN.R-project.org/package=gplots via the Internet. Accessed July 1, 2011. (2009) Gplots: Various R programming tools for plotting data. 2.7.4. Available: http://CRAN.R-project.org/package=gplots via the Internet. Accessed July 1, 2011.

25. Liu J, Enomoto S, Lancto CA, Abrahamsen MS, Rutherford MS (2008) Inhibition of apoptosis in Cryptosporidium parvum infected intestinal epithelial cells is dependent on survivin. Infect Immun 76: 3784–3792.J. LiuS. EnomotoCA LanctoMS AbrahamsenMS Rutherford2008Inhibition of apoptosis in Cryptosporidium parvum infected intestinal epithelial cells is dependent on survivin.Infect Immun7637843792

26. Templeton TJ, Iyer LM, Anantharaman V, Enomoto S, Abrahante JE, et al. (2004) Comparative analysis of apicomplexa and genomic diversity in eukaryotes. Genome Res 14: 1686–1695.TJ TempletonLM IyerV. AnantharamanS. EnomotoJE Abrahante2004Comparative analysis of apicomplexa and genomic diversity in eukaryotes.Genome Res1416861695

27. Kaufman L, Rousseeuw PJ (1990) Finding groups in data: An introduction to cluster analysis. New York: John Wiley & Sons, Inc. 342 p.L. KaufmanPJ Rousseeuw1990Finding groups in data: An introduction to cluster analysisNew YorkJohn Wiley & Sons, Inc342

28. Borowski H, Thompson RC, Armstrong T, Clode PL (2010) Morphological characterization of Cryptosporidium parvum life-cycle stages in an in vitro model system. Parasitology 137: 13–26.H. BorowskiRC ThompsonT. ArmstrongPL Clode2010Morphological characterization of Cryptosporidium parvum life-cycle stages in an in vitro model system.Parasitology1371326

29. Chatterjee A, Banerjee S, Steffen M, O'Connor RM, Ward HD, et al. (2010) Evidence for mucin-like glycoproteins that tether sporozoites of Cryptosporidium parvum to the inner surface of the oocyst wall. Eukaryot Cell 9: 84–96.A. ChatterjeeS. BanerjeeM. SteffenRM O'ConnorHD Ward2010Evidence for mucin-like glycoproteins that tether sporozoites of Cryptosporidium parvum to the inner surface of the oocyst wall.Eukaryot Cell98496

30. Dubey JP, Speer CA, Fayer R (1990) Cryptosporidiosis of man and animals. Boca Raton, FL: CRC Press. JP DubeyCA SpeerR. Fayer1990Cryptosporidiosis of man and animalsBoca Raton, FLCRC Press

31. SGD Project. Saccharomyces genome database (2010) SGD Project. Saccharomyces genome database2010Available: http://www.yeastgenome.org/cache/genomeSnapshot.html. Accessed January 20, 2011. Available: http://www.yeastgenome.org/cache/genomeSnapshot.html. Accessed January 20, 2011.

32. Joe A, Verdon R, Tzipori S, Keusch GT, Ward HD (1998) Attachment of Cryptosporidium parvum sporozoites to human intestinal epithelial cells. Infect Immun 66: 3429–3432.A. JoeR. VerdonS. TziporiGT KeuschHD Ward1998Attachment of Cryptosporidium parvum sporozoites to human intestinal epithelial cells.Infect Immun6634293432

33. Jenkins MC, O'Brien C, Miska K, Schwarz RS, Karns J, et al. (2011) Gene expression during excystation of Cryptosporidium parvum oocysts. Parasitol Res 109: 509–513.MC JenkinsC. O'BrienK. MiskaRS SchwarzJ. Karns2011Gene expression during excystation of Cryptosporidium parvum oocysts.Parasitol Res109509513

34. Behnke MS, Wootton JC, Lehmann MM, Radke JB, Lucas O, et al. (2010) Coordinated progression through two subtranscriptomes underlies the tachyzoite cycle of Toxoplasma gondii. PLoS One 5: e12354.MS BehnkeJC WoottonMM LehmannJB RadkeO. Lucas2010Coordinated progression through two subtranscriptomes underlies the tachyzoite cycle of Toxoplasma gondii.PLoS One5e12354

35. Gritzmacher CA, Reese RT (1984) Protein and nucleic acid synthesis during synchronized growth of Plasmodium falciparum. J Bacteriol 160: 1165–1167.CA GritzmacherRT Reese1984Protein and nucleic acid synthesis during synchronized growth of Plasmodium falciparum.J Bacteriol16011651167

36. de Rojas MO, Wasserman M (1985) Temporal relationships on macromolecular synthesis during the asexual cell cycle of Plasmodium falciparum. Trans R Soc Trop Med Hyg 79: 792–796.MO de RojasM. Wasserman1985Temporal relationships on macromolecular synthesis during the asexual cell cycle of Plasmodium falciparum.Trans R Soc Trop Med Hyg79792796

37. Le Roch KG, Zhou Y, Blair PL, Grainger M, Moch JK, et al. (2003) Discovery of gene function by expression profiling of the malaria parasite life cycle. Science 301: 1503–1508.KG Le RochY. ZhouPL BlairM. GraingerJK Moch2003Discovery of gene function by expression profiling of the malaria parasite life cycle.Science30115031508

38. Putignani L, Possenti A, Cherchi S, Pozio E, Crisanti A, et al. (2008) The thrombospondin-related protein CpMIC1 (CpTSP8) belongs to the repertoire of micronemal proteins of Cryptosporidium parvum. Mol Biochem Parasitol 157: 98–101.L. PutignaniA. PossentiS. CherchiE. PozioA. Crisanti2008The thrombospondin-related protein CpMIC1 (CpTSP8) belongs to the repertoire of micronemal proteins of Cryptosporidium parvum.Mol Biochem Parasitol15798101

39. Tosini F, Agnoli A, Mele R, Gomez Morales MA, Pozio E (2004) A new modular protein of Cryptosporidium parvum, with ricin B and LCCL domains, expressed in the sporozoite invasive stage. Mol Biochem Parasitol 134: 137–147.F. TosiniA. AgnoliR. MeleMA Gomez MoralesE. Pozio2004A new modular protein of Cryptosporidium parvum, with ricin B and LCCL domains, expressed in the sporozoite invasive stage.Mol Biochem Parasitol134137147

40. Fayer R (1997) Cryptosporidium and cryptosporidiosis. Boca Raton, FL: CRC Press. R. Fayer1997Cryptosporidium and cryptosporidiosisBoca Raton, FLCRC Press

41. Jorgensen P, Rupes I, Sharom JR, Schneper L, Broach JR, et al. (2004) A dynamic transcriptional network communicates growth potential to ribosome synthesis and critical cell size. Genes Dev 18: 2491–2505.P. JorgensenI. RupesJR SharomL. SchneperJR Broach2004A dynamic transcriptional network communicates growth potential to ribosome synthesis and critical cell size.Genes Dev1824912505

42. Bernstein KA, Bleichert F, Bean JM, Cross FR, Baserga SJ (2007) Ribosome biogenesis is sensed at the start cell cycle checkpoint. Mol Biol Cell 18: 953–964.KA BernsteinF. BleichertJM BeanFR CrossSJ Baserga2007Ribosome biogenesis is sensed at the start cell cycle checkpoint.Mol Biol Cell18953964

43. Moore SA (1988) Kinetic evidence for a critical rate of protein synthesis in the saccharomyces cerevisiae yeast cell cycle. J Biol Chem 263: 9674–9681.SA Moore1988Kinetic evidence for a critical rate of protein synthesis in the saccharomyces cerevisiae yeast cell cycle.J Biol Chem26396749681

44. Li JL, Cox LS (2003) Characterisation of a sexual stage-specific gene encoding ORC1 homologue in the human malaria parasite Plasmodium falciparum. Parasitol Int 52(1): 41–52.JL LiLS Cox2003Characterisation of a sexual stage-specific gene encoding ORC1 homologue in the human malaria parasite Plasmodium falciparum.Parasitol Int5214152

45. Mehra P, Biswas AK, Gupta A, Gourinath S, Chitnis CE, et al. (2005) Expression and characterization of human malaria parasite Plasmodium falciparum origin recognition complex subunit 1. Biochem Biophys Res Commun 337: 955–966.P. MehraAK BiswasA. GuptaS. GourinathCE Chitnis2005Expression and characterization of human malaria parasite Plasmodium falciparum origin recognition complex subunit 1.Biochem Biophys Res Commun337955966

46. Chen F, Mackey AJ, Stoeckert CJ Jr, Roos DS (2006) OrthoMCL-DB: Querying a comprehensive multi-species collection of ortholog groups. 34. (Database issue): Nucleic Acids Res. pp. D363–8.F. ChenAJ MackeyCJ Stoeckert JrDS Roos2006OrthoMCL-DB: Querying a comprehensive multi-species collection of ortholog groups34Database issueNucleic Acids ResD3638

47. Agop-Nersesian C, Naissant B, Ben Rached F, Rauch M, Kretzschmar A, et al. (2009) Rab11A-controlled assembly of the inner membrane complex is required for completion of apicomplexan cytokinesis. PLoS Pathog 5(1): e1000270.C. Agop-NersesianB. NaissantF. Ben RachedM. RauchA. Kretzschmar2009Rab11A-controlled assembly of the inner membrane complex is required for completion of apicomplexan cytokinesis.PLoS Pathog51e1000270

48. Maiorano D, Lutzmann M, Mechali M (2006) MCM proteins and DNA replication. Curr Opin Cell Biol 18: 130–136.D. MaioranoM. LutzmannM. Mechali2006MCM proteins and DNA replication.Curr Opin Cell Biol18130136

49. Kolodner RD, Marsischky GT (1999) Eukaryotic DNA mismatch repair. Curr Opin Genet Dev 9: 89–96.RD KolodnerGT Marsischky1999Eukaryotic DNA mismatch repair.Curr Opin Genet Dev98996

50. Striepen B, Jordan CN, Reiff S, van Dooren GG (2007) Building the perfect parasite: Cell division in apicomplexa. PLoS Pathog 3(6): e78.B. StriepenCN JordanS. ReiffGG van Dooren2007Building the perfect parasite: Cell division in apicomplexa.PLoS Pathog36e78

51. Gordon JL, Beatty WL, Sibley LD (2008) A novel actin-related protein is associated with daughter cell formation in Toxoplasma gondii. Eukaryot Cell 7: 1500–1512.JL GordonWL BeattyLD Sibley2008A novel actin-related protein is associated with daughter cell formation in Toxoplasma gondii.Eukaryot Cell715001512

52. Santos JM, Lebrun M, Daher W, Soldati D, Dubremetz JF (2009) Apicomplexan cytoskeleton and motors: Key regulators in morphogenesis, cell division, transport and motility. Int J Parasitol 39: 153–162.JM SantosM. LebrunW. DaherD. SoldatiJF Dubremetz2009Apicomplexan cytoskeleton and motors: Key regulators in morphogenesis, cell division, transport and motility.Int J Parasitol39153162

53. Hall N, Carlton J (2005) Comparative genomics of malaria parasites. Curr Opin Genet Dev 15: 609–613.N. HallJ. Carlton2005Comparative genomics of malaria parasites.Curr Opin Genet Dev15609613

54. Cleary MD, Singh U, Blader IJ, Brewer JL, Boothroyd JC (2002) Toxoplasma gondii asexual development: Identification of developmentally regulated genes and distinct patterns of gene expression. Eukaryot Cell 1: 329–340.MD ClearyU. SinghIJ BladerJL BrewerJC Boothroyd2002Toxoplasma gondii asexual development: Identification of developmentally regulated genes and distinct patterns of gene expression.Eukaryot Cell1329340

55. Denton H, Roberts CW, Alexander J, Thong KW, Coombs GH (1996) Enzymes of energy metabolism in the bradyzoites and tachyzoites of Toxoplasma gondii. FEMS Microbiol Lett 137: 103–108.H. DentonCW RobertsJ. AlexanderKW ThongGH Coombs1996Enzymes of energy metabolism in the bradyzoites and tachyzoites of Toxoplasma gondii.FEMS Microbiol Lett137103108

56. Yang S, Parmley SF (1995) A bradyzoite stage-specifically expressed gene of Toxoplasma gondii encodes a polypeptide homologous to lactate dehydrogenase. Mol Biochem Parasitol 73: 291–294.S. YangSF Parmley1995A bradyzoite stage-specifically expressed gene of Toxoplasma gondii encodes a polypeptide homologous to lactate dehydrogenase.Mol Biochem Parasitol73291294

57. Deng M, Lancto CA, Abrahamsen MS (2004) Cryptosporidium parvum regulation of human epithelial cell gene expression. Int J Parasitol 34: 73–82.M. DengCA LanctoMS Abrahamsen2004Cryptosporidium parvum regulation of human epithelial cell gene expression.Int J Parasitol347382

58. Xu T, Ping J, Yu Y, Yu F, Yu Y, et al. (2010) Revealing parasite influence in metabolic pathways in apicomplexa infected patients. BMC Bioinformatics 11: Suppl 11S13.T. XuJ. PingY. YuF. YuY. Yu2010Revealing parasite influence in metabolic pathways in apicomplexa infected patients.BMC Bioinformatics11Suppl 11S13

59. Hirano T, Mitchison TJ (1994) A heterodimeric coiled-coil protein required for mitotic chromosome condensation in vitro. Cell 79: 449–458.T. HiranoTJ Mitchison1994A heterodimeric coiled-coil protein required for mitotic chromosome condensation in vitro.Cell79449458

60. Sutani T, Yuasa T, Tomonaga T, Dohmae N, Takio K, et al. (1999) Fission yeast condensin complex: Essential roles of non-SMC subunits for condensation and Cdc2 phosphorylation of Cut3/SMC4. Genes Dev 13: 2271–2283.T. SutaniT. YuasaT. TomonagaN. DohmaeK. Takio1999Fission yeast condensin complex: Essential roles of non-SMC subunits for condensation and Cdc2 phosphorylation of Cut3/SMC4.Genes Dev1322712283

61. Bell SP, Dutta A (2002) DNA replication in eukaryotic cells. Annu Rev Biochem 71: 333–374.SP BellA. Dutta2002DNA replication in eukaryotic cells.Annu Rev Biochem71333374

62. Patterson S, Robert C, Whittle C, Chakrabarti R, Doerig C, et al. (2006) Pre-replication complex organization in the atypical DNA replication cycle of Plasmodium falciparum: Characterization of the mini-chromosome maintenance (MCM) complex formation. Mol Biochem Parasitol 145: 50–59.S. PattersonC. RobertC. WhittleR. ChakrabartiC. Doerig2006Pre-replication complex organization in the atypical DNA replication cycle of Plasmodium falciparum: Characterization of the mini-chromosome maintenance (MCM) complex formation.Mol Biochem Parasitol145505910.1016/j.molbiopara.2005.09.006. 10.1016/j.molbiopara.2005.09.006.

63. Balaji S, Babu MM, Iyer LM, Aravind L (2005) Discovery of the principal specific transcription factors of apicomplexa and their implication for the evolution of the AP2-integrase DNA binding domains. Nucleic Acids Res 33: 3994–4006.S. BalajiMM BabuLM IyerL. Aravind2005Discovery of the principal specific transcription factors of apicomplexa and their implication for the evolution of the AP2-integrase DNA binding domains.Nucleic Acids Res3339944006

64. Yuda M, Iwanaga S, Shigenobu S, Kato T, Kaneko I (2010) Transcription factor AP2-sp and its target genes in malarial sporozoites. Mol Microbiol 75: 854–863.M. YudaS. IwanagaS. ShigenobuT. KatoI. Kaneko2010Transcription factor AP2-sp and its target genes in malarial sporozoites.Mol Microbiol7585486310.1111/j.1365-2958.2009.07005.x. 10.1111/j.1365-2958.2009.07005.x.

65. Heiges M, Wang H, Robinson E, Aurrecoechea C, Gao X, et al. (2006) CryptoDB: A Cryptosporidium bioinformatics resource update. 34. (Database issue): Nucleic Acids Res. pp. D419–22.M. HeigesH. WangE. RobinsonC. AurrecoecheaX. Gao2006CryptoDB: A Cryptosporidium bioinformatics resource update34Database issueNucleic Acids ResD41922

66. Duncker BP, Chesnokov IN, McConkey BJ (2009) The origin recognition complex protein family. Genome Biol 10(3): 214.BP DunckerIN ChesnokovBJ McConkey2009The origin recognition complex protein family.Genome Biol103214

67. Dueber EC, Costa A, Corn JE, Bell SD, Berger JM (2011) Molecular determinants of origin discrimination by Orc1 initiators in archaea. Nucleic Acids Res 39: 3621–3631.EC DueberA. CostaJE CornSD BellJM Berger2011Molecular determinants of origin discrimination by Orc1 initiators in archaea.Nucleic Acids Res3936213631

68. Grainge I, Gaudier M, Schuwirth BS, Westcott SL, Sandall J, et al. (2006) Biochemical analysis of a DNA replication origin in the archaeon Aeropyrum pernix. J Mol Biol 363: 355–369. 1.I. GraingeM. GaudierBS SchuwirthSL WestcottJ. Sandall2006Biochemical analysis of a DNA replication origin in the archaeon Aeropyrum pernix.J Mol Biol363355369. 1

69. Bishop R, Shah T, Pelle R, Hoyle D, Pearson T, et al. (2005) Analysis of the transcriptome of the protozoan Theileria parva using MPSS reveals that the majority of genes are transcriptionally active in the schizont stage. Nucleic Acids Res 33: 5503–5511.R. BishopT. ShahR. PelleD. HoyleT. Pearson2005Analysis of the transcriptome of the protozoan Theileria parva using MPSS reveals that the majority of genes are transcriptionally active in the schizont stage.Nucleic Acids Res3355035511

70. Bozdech Z, Llinas M, Pulliam BL, Wong ED, Zhu J, et al. (2003) The transcriptome of the intraerythrocytic developmental cycle of Plasmodium falciparum. PLoS Biol 1(1): E5.Z. BozdechM. LlinasBL PulliamED WongJ. Zhu2003The transcriptome of the intraerythrocytic developmental cycle of Plasmodium falciparum.PLoS Biol11E5

71. Behnke MS, Radke JB, Smith AT, Sullivan WJ Jr, White MW (2008) The transcription of bradyzoite genes in Toxoplasma gondii is controlled by autonomous promoter elements. Mol Microbiol 68: 1502–1518.MS BehnkeJB RadkeAT SmithWJ Sullivan JrMW White2008The transcription of bradyzoite genes in Toxoplasma gondii is controlled by autonomous promoter elements.Mol Microbiol6815021518

72. Oliveros JC (2007) VENNY. JC Oliveros2007VENNY.an interactive tool for comparing lists with venn diagrams. May 1, 2011. an interactive tool for comparing lists with venn diagrams. May 1, 2011.