The steady-state transcriptome of the four major life-cycle stages of Trypanosoma cruzi
© Minning et al; licensee BioMed Central Ltd. 2009
Received: 09 April 2009
Accepted: 07 August 2009
Published: 07 August 2009
Chronic chagasic cardiomyopathy is a debilitating and frequently fatal outcome of human infection with the protozoan parasite, Trypanosoma cruzi. Microarray analysis of gene expression during the T. cruzi life-cycle could be a valuable means of identifying drug and vaccine targets based on their appropriate expression patterns, but results from previous microarray studies in T. cruzi and related kinetoplastid parasites have suggested that the transcript abundances of most genes in these organisms do not vary significantly between life-cycle stages.
In this study, we used whole genome, oligonucleotide microarrays to globally determine the extent to which T. cruzi regulates mRNA relative abundances over the course of its complete life-cycle. In contrast to previous microarray studies in kinetoplastids, we observed that relative transcript abundances for over 50% of the genes detected on the T. cruzi microarrays were significantly regulated during the T. cruzi life-cycle. The significant regulation of 25 of these genes was confirmed by quantitative reverse-transcriptase PCR (qRT-PCR). The T. cruzi transcriptome also mirrored published protein expression data for several functional groups. Among the differentially regulated genes were members of paralog clusters, nearly 10% of which showed divergent expression patterns between cluster members.
Taken together, these data support the conclusion that transcript abundance is an important level of gene expression regulation in T. cruzi. Thus, microarray analysis is a valuable screening tool for identifying stage-regulated T. cruzi genes and metabolic pathways.
Trypanosoma cruzi, the etiologic agent of Chagas disease in humans, is a protozoan parasite which assumes four morphological stages during its cycle in insect and mammalian hosts. In the reduviid insect vector, T. cruzi epimastigotes replicate extracellularly in the lumen of the gut. As the parasite stages reach the posterior end of the gut they attach to the wall of the rectum and convert to non-replicative, infective metacyclic trypomastigotes that are released in the feces when the bug takes a blood meal. The metacyclic trypomastigotes enter the mammalian host when they are rubbed into the bite wound or other open skin or mucosal surfaces and invade host cells. The intracytoplasmic parasites convert to the replicative, amastigote stage and undergo many rounds of division before transforming into elongated, motile trypomastigote stages which are released when the host cell ruptures. The trypomastigotes are disseminated in the blood and lymph where they may infect virtually any nucleated cell or be taken up by the insect vector to complete the life cycle. A detailed description of the T. cruzi life-cycle with accompanying diagram can be accessed at the web site for the Centers for Disease Control and Prevention http://www.dpd.cdc.gov/dpdx/html/trypanosomiasisamerican.htm.
There are no vaccines for Chagas disease, and specific chemotherapy for T. cruzi infection is hindered by adverse side-effects and questionable efficacy. In order to broadly screen for genes with relevant expression patterns, and to learn more about the biology of T. cruzi in general, we and others have conducted transcriptomic and proteomic analyses of T. cruzi life-cycle stages [1–11]. Although these analyses have highlighted potential drug and vaccine targets, they have been somewhat limited in scope, largely non-quantitative, and rarely correlated transcript and protein abundances. The latter issue is of particular importance for T. cruzi and other kinetoplastids, due to the generally accepted view that regulation of gene expression in the kinetoplastids is almost entirely post-transcriptional (reviewed in [12, 13]). In fact, only one kinetoplastid RNA pol II promoter has been described despite repeated attempts in many laboratories (reviewed i ). Other studies have implicated mRNA processing , translational repression [16–18], polysome recruitment , and codon adaptation  in the regulation of gene expression in the kinetoplastids, all processes that would be predicted to mitigate the role of mRNA abundance regulation in determining protein expression levels. Indeed, previous microarray studies in the kinetoplastids have revealed relatively modest numbers of genes whose transcript abundances are significantly regulated (reviewed in ). There is, however, plentiful evidence that transcript levels in these parasites are controlled by mRNA decay, involving 3'-UTRs and RNA-binding proteins [19, 22–28]. To determine the extent of mRNA abundance regulation in T. cruzi globally, we performed whole-genome, DNA microarray analysis of the four life-cycle stages of T. cruzi. The results of these analyses were compared with existing protein expression data for T. cruzi to determine the correlation, if any, between transcript and protein relative abundances in this human pathogen.
T. cruzi displays significant stage-regulation of relative transcript abundances for thousands of its genes
With or without the contribution of large gene family members, greater than 50% of the Trypanosoma cruzi genes detected on whole-genome, oligonucleotide microarrays were significantly stage-regulated at the RNA level.
# spots (percent of total)
# genes (percent of total)
large gene families3
minus large gene families
Of the genes upregulated in two stages, 76% (1083/1423), were co-upregulated in stages occurring in the same host (mammalian hosts – amastigotes and trypomastigotes; insect hosts – epimastigotes and metacyclic trypomastigotes), whereas approximately 20% (282/1423) were co-upregulated in stages with similar biological functions. For example, genes co-upregulated in amastigotes and epimastigotes (dividing stages) were enriched in those involved in DNA repair, including proliferative cell nuclear antigen (PCNA). Likewise, genes co-upregulated in trypomastigotes and metacyclic trypomastigotes (non-dividing, infective stages) included never in mitosis (NIMA) related kinase and were enriched for trans-sialidases (TS), which are involved in invasion (reviewed in ).
The validity of T. cruzi microarray data is supported by qRT-PCR and by comparison to existing proteomic data
The accuracy of the microarray results was confirmed by quantitative reverse-transcriptase PCR (qRT-PCR) for a subset of the significantly regulated genes.
The T. cruzi microarray data provide expanded coverage of important functional groups
The concordance between the presently reported transcriptomic data and published proteomic data suggests that the mRNA relative abundances observed here will be useful indicators of the protein expression levels of genes for which protein expression data has not been reported, including many genes annotated as 'hypothetical.' Current proteomic methods often fail to detect or quantify proteins that are expressed at low levels, (e.g. protein kinases and phosphatases; reviewed in ) or proteins that have physical properties making them less amenable to proteomic analyses, (e.g. membrane proteins; reviewed in ). For example, this transcriptomic analysis detected 30 protein phosphatases and 128 protein kinases, which were significantly regulated at the transcript level (Additional files 1 and 2), and few of which have published protein expression data. Likewise, our transcriptome analysis detected 166 significantly regulated genes annotated as 'integral-to-membrane' (Additional file 3).
Members of many T. cruzi paralog clusters exhibit significant stage-regulated divergence in mRNA relative abundance patterns
In addition to increasing transcript abundance, another presumed function of gene duplication is to allow for the development of new protein functions . Such neofunctionalization may have occurred in the T. cruzi 60S ribosomal L18 protein paralog group (Figure 6). Tc00.1047053503395.40, which was upregulated in trypomastigotes, has a 130 amino acid N-terminal extension not found in any other publicly-available 60S ribosomal L18 sequence, suggesting the possibility that this paralog has a novel function in T. cruzi trypomastigotes.
Transcriptomic analyses have not been as aggressively pursued in kinetoplastid organisms as they have been in other protozoan parasites, such as Plasmodium. This is likely due to the view that the predominantly post-transcriptional nature of gene expression regulation in kinetoplastids makes microarray studies in these organisms generally less informative than microarray studies in organisms for which transcription initiation plays a larger role in gene expression regulation. In support of this, the ranges of observed mRNA ratios between life-cycle stages in kinetoplastid organisms have been relatively narrow compared to the ranges observed in organisms with canonical RNA pol II promoters. For example, 1,000-fold induction in relative transcript abundance for CD70 was observed in a microarray study comparing Human T-lymphotropic virus type 1 (HTLV-1) carrying T-cell lines versus HTLV-1-negative T-cell lines , whereas in kinetoplastid microarray studies, fold inductions above ~8-fold are rarely observed [3, 42–46].
Further contributing to the debate over the utility of microarray studies in kinetoplastids has been the general observation that relatively few genes have exhibited significant stage regulation of mRNA relative abundances in previous microarray studies. For instance, Diehl et al., using arrays of 21,024 PCR-amplified genome shot-gun library inserts from T. brucei to compare relative transcript abundances between in vitro cultured human, long slender forms and procyclic parasites, observed that, although 75% of the array elements detected transcripts, only 2% displayed significant differences between the two life-cycle stages . The observation that only 2% of the detected transcripts differed significantly in relative abundance between two very different life-cycle stages supported the concept that transcript abundances in T. brucei are largely constitutive, especially considering the two life-cycle stages were compared directly on the microarrays. This picture of predominantly constitutive genome expression at the level of transcript abundance in kinetoplastids was reinforced manifold by subsequent studies in T. brucei  and Leishmania spp. (reviewed in ).
In light of the findings from T. brucei and Leishmania spp. microarray studies as well as findings from our previous, limited microarray study in T. cruzi , we sought to determine if transcript levels are globally variable between life-cycle stages of T. cruzi using whole-genome oligonucleotide microarrays. When cDNAs from each life-cycle stage were co-hybridized with a reference cDNA sample comprised of all four life-cycle stages on oligonucleotide, whole genome microarrays, we observed that the relative transcript abundances for over 50% of the genes detected on the arrays were significantly regulated between the life-cycle stages.
The microarray data for selected genes were validated by qRT-PCR. There was generally good agreement between the two platforms in terms of direction of regulation, however there were some quantitative differences, especially in terms of the extent of up or downregulation estimated by these two techniques. These quantitative differences likely were the result of sequence-specific effects and differences in dynamic range for the two platforms [49, 50]. In addition, the RNA samples used for qRT-PCR were not the identical samples used for the microarrays.
Recent improvements in microarray design, such as the use of oligonucleotide probes designed for more uniform hybridization kinetics and with lower likelihood of cross-hybridization than amplicon-derived probes and the use of microarrays with greater genome coverage could in part account for why our study identified a higher percentage of significantly stage-regulated genes than previous microarray studies in other kinetoplastids [1–4, 39, 42]. However, even a recent study by Rochette et al. that used Leishmania whole genome oligonucleotide microarrays to directly compare amastigotes and procyclics reported that only 7% of the genes in Leishmania infantum and 9.3% of the genes in L. major were developmentally regulated . It is also possible that the use of different statistical methods for analysis of microarray data my yield widely different estimates of the frequency of regulated genes. However, using less rigorous statistical criteria for determining significance, such as fold-change, results in higher false discovery rates  and thus would be predicted to overestimate the number of significantly regulated genes. Our reanalysis of the Papadopoulou data from L. infantum using our statistical methods found 715 genes (7% of the genes fulfilling our criteria for detection) to be significantly regulated between amastigotes and promastigotes with an FDR of 1.55% (an FDR nearly 4 times greater than what we report for T. cruzi). Thus, different statistical analysis methods cannot account for the differential estimates of the frequency of regulated genes in Leishmania and T. cruzi.
By incorporating all four life cycle stages of T. cruzi in our analysis, in contrast to the 2 stage comparisons done with other kinetoplastids (e.g. procyclic to amastigote  or procyclic to bloodstream forms ) we have increased the possibility of identifying a greater number of genes regulated in expression in at least one of the 4 life cycle stages. However, the contribution of increased number of stages studied to the overall number of significantly regulated genes was greatly mitigated by the fact that we used a reference design in which each life-cycle stage was compared to a mixture of all four life-cycle stages. This design minimized the number of microarrays required to compare the four life-cycle stages (each stage vs. reference = 24 hybridizations in the present study, each stage vs. each other stage = 36 hybridizations for equivalent replication), and also simplified the analysis, but had the disadvantage of buffering the ratios observed. By contrast, previous microarray studies in the kinetoplastids involved directly comparing one life-cycle stage with another on the microarrays, an experimental design which should maximize the number of significant ratios. Thus our estimate that >50% of genes in T. cruzi have significantly stage-regulated transcript levels is much more likely to be an underestimate than an overestimate.
The number of significantly stage-regulated transcripts reported here for T. cruzi is also likely to be underestimated because >2,000 of the oligonucleotide probes detected transcripts in two or fewer stages and were removed from further analysis for statistical reasons (significance determination for such genes would have been dubious, because SAM analysis of four groups (life-cycle stages) was unreliable for genes present in less than three groups). Also, almost certainly some of the genes that we reported as non-significantly regulated are significantly regulated in transitional stages of the life-cycle. For instance, Saxena et al., using microarrays of ~8300 PCR-amplified genome survey sequencing (GSS) clone inserts, identified 344 protein coding genes significantly regulated during axenic promastigote-to-amastigote differentiation in Leishmania donovani . Of these, 136 genes, nearly 40% of the total number of differentially regulated genes identified in the study, displayed transient up or downregulation that would have been missed by only comparing fully differentiated parasites. Thus our analysis of only fully differentiated stages likely underestimated the percentage of genes we found significantly regulated at the transcript level during the T. cruzi life cycle, suggesting that time course studies looking at transitional stages in T. cruzi would be fruitful.
The disparity in the degree of transcript abundance regulation between T. cruzi and the other sequenced kinetoplastids suggests there may be differences in the number and kind of RNA binding proteins in their genomes, because almost all regulation of mRNA abundance in these organisms is at the level of mRNA stability. However, the genomes of T. cruzi, T. brucei, and L. major all have remarkably similar complements of RNA binding proteins (reviewed in ). Of the 77 proteins with RNA Recognition Motifs (RRMs) in T. cruzi, T. brucei, and L. major, only two were unique to T. cruzi, RBP4 (Tc00.1047053508901.20) and DRBD8 (Tc00.1047053503709.10, Tc00.1047053509581.50), both of which were upregulated in metacyclic trypomastigotes.
Present knowledge of the kinetoplastids does not elicit an obvious biological explanation for why gene expression regulation in T. cruz i might be so different from its closest relatives. Perhaps the most striking difference between T. cruzi and T. brucei and L. major is that only T. cruzi has both intracellular and extracellular life-cycle stages in the mammalian host. Amastigotes must replicate intracellularly and, although sequestered from the actions of antibodies and complement, are subjected to attack by CD8+ cytotoxic T lymphocytes (CTL)  and to the effects of host cytokines  and antimicrobials, including reactive oxygen species and nitric oxide . Trypomastigotes are extracellular and thus must evade host defenses including complement, antibodies, and phagocytes. An enhanced capacity for regulating transcript abundances may have provided T. cruzi with the agility necessary to rapidly shift between such different environments.
The finding that T. cruzi divergently regulates mRNA relative abundances for members of paralog clusters was unexpected, given the post-transcriptional mode of gene expression regulation in this parasite. Divergent expression of duplicate genes is believed to occur by the accumulation of mutations in non-coding cis regulatory sequences, such as promoters . Since kinetoplastids appear to be largely devoid of canonical promoters, it will be interesting to study the 3' untranslated regions of paralogs with divergent expression patterns to identify the sequences responsible for the differential mRNA abundances.
The differentially regulated T. cruzi duplicate genes identified here will be useful to study the retention and neofunctionalization of duplicate genes in an ancient eukaryote lacking transcriptional control of gene expression. In such an environment, the contributions of mRNA stability control and sub-cellular localization to the evolution of new genes can be studied in isolation from the effects of transcriptional control. Gene knockout studies are currently underway to determine if the array-identified, divergently-expressed paralogs, such as the trypomastigote upregulated 60S ribosomal L18 protein, have novel functions.
T. cruzi displays stage-regulated control of mRNA abundances for a significant proportion of its genes despite the fact that regulation of gene expression in the kinetoplastids is primarily post-transcriptional and despite the fact that previous genomic studies of the kinetoplastids have found that only a low proportion of their genomes are stage-regulated. Moreover, T. cruzi oligonucleotide microarray measurements largely agree with known protein expression data for key functional groups. Taken together, T. cruzi oligonucleotide microarray analysis is a useful screening tool to determine stage-regulated gene expression in this human pathogen.
Brazil strain T. cruzi trypomastigotes were grown in monolayers of Vero cells (ATCC no. CCL-81) in RPMI supplemented with 5% horse serum as previously described . Emergent trypomastigotes were harvested daily and examined by light microscopy to determine the percentages of amastigotes and trypomastigotes. Only preparations containing > 95% trypomastigotes were used in the subsequent studies. Amastigotes were prepared from axenically-induced trypomastigotes as described previously . Briefly, emergent trypomastigote samples were centrifuged at 3,000 × g for 15 m at room temperature, brought to density of 5 × 106/ml in Protein-Free Hybridoma Medium (PFHM-II, Gibco, Bethesda, MD) at pH 5.0, and incubated at 37°C until greater than 95% of the parasites were fully converted to the amastigote stage as determined by microscopic examination. T. cruzi epimastigotes were grown in Liver Infusion Tryptose media (LIT) as previously described . Cultures were harvested during mid-log phase by centrifugation at 3,000 × g for 10 m at room temperature. Metacyclic trypomastigotes were obtained from epimastigotes by axenic induction as previously described . Briefly, log phase epimastigote cultures were centrifuged at 3,000 × g for 15 m at room temperature, brought to a density of 5 × 106/ml in Complete Grace's Insect Medium (Sigma no. G8142, Sigma, St. Louis, MO) supplemented with 10% fetal bovine serum, pH 6.6, and incubated at 29°C for 10–14 d. The percentages of metacyclics were determined by microscopic examination of parasites stained with Dif-Quick (Baxter Diagnostics, McGaw Park, IL).
Preparation of labeled cDNAs
Total RNA was isolated from parasites with the High Pure Total RNA Isolation Kit (Roche) per the manufacturer's instructions, including on-column DNase I digestion. Labeled first strand cDNA was synthesized from the total RNA samples using Superscript™ II reverse transcriptase (Life Technologies, Grand Island, NY) as previously described . Briefly, 10 μg RNA were combined with 2 μg oligo-d(T) (Life Technologies), brought to 10 μl with water, denatured at 70°C for 10 m, and cooled on ice. Final 30 ul reactions were assembled with final concentrations of 1× First Strand Buffer (Life Technologies), 10 mM DTT (Life Technologies), 0.5 mM dATP, dCTP, and dGTP (Amersham Pharmacia Biotech, Piscataway, NJ), 0.1 mM dTTP (Amersham Pharmacia Biotech), 0.1 mM Cy3- or Cy5-dUTP (Amersham Pharmacia Biotech), and 13.3 U/μl Superscript II reverse transcriptase (Life Technologies). The reaction mixtures were incubated at 42°C for 2 h and stopped with the addition of 1.5 μL 20 mM EDTA. To degrade the RNA template, 1.5 μL 500 mM NaOH were added and the reactions were heated at 70°C for 10 m. The reactions were neutralized by adding 1.5 μL 500 mM HCl, and unincorporated fluorescent nucleotides were removed using GFX columns per the manufacturer's instructions (Amersham Pharmacia Biotech). The purified products were eluted using 50 μL TE, pH 8.0, and the labeled cDNA was completely dried in a SpeedVac (Savant Instruments, Holbrook, NY, USA) and resuspended in 10 μL water.
The T. cruzi microarrays were obtained from the Pathogen Functional Genomics Resource Center (PFGRC). The microarray description is available at http://pfgrc.tigr.org. The 12,288 unique array oligonucleotides were sense-strand 70-mers designed against open reading frames in the annotated CL Brener reference genome sequence and were printed in duplicate. The microarrays contained an additional 500 control oligonucleotides designed from Arabidopsis sequences, also printed in duplicate.
To assure accuracy in the spot-to-gene annotation, the sequences of the oligonucleotides used on the TIGR T. cruzi microarrays were remapped onto the T. cruzi genome and assigned to genes. The mapping was complicated by the fact that T. cruzi is diploid and the CL-Brener strain used for genome sequencing is a hybrid strain . Thus, the two alleles for genes in the T. cruz i genome frequently differ by 1–2% in their coding sequences , and the terms 'allele' and 'gene' are conflated when describing the T. cruzi genome. Moreover, complete chromosomes were never assembled for the T. cruzi genome , and triploidy for some loci has been confirmed . As a result, the total number of 'genes' identified by the oligonucleotides on the microarrays was greater than the number of unique oligonucleotides. Thus, oligonucleotides mapping with greater than 80% homology to three or fewer 'genes' were used in the subsequent analyses, and oligonucleotides that mapped to more than three genes were not included in the enumeration of significantly regulated genes. Additional file 4 shows the results of our remapping of the oligonucleotides to the T. cruzi annotated genome, the original TIGR mapping, and the results of the cross-hybridization prediction. Additional file 5 shows the representation of large gene families (cluster designations per ) on the microarrays. The mucin associated surface protein (MASP) family in particular was poorly represented on the microarrays.
The microarrays were prehybridized immediately preceding hybridization as described on the PFGRC website. Briefly, arrays are incubated in prehybridization solution (4× SSC, 40% formamide, 0.1% SDS, 0.5% BSA, and 25 mM Tris-HCl, pH 8.0) at 42°C for 1 h. Slides are then washed at room temperature in 0.2% SDS, followed by 3 washes with water, and dried by centrifugation.
Six hybridizations were performed for each life-cycle stage. The hybridizations consisted of three dye-swap experiments from three independent samples (biological replicates). In each case, the experimental sample was from a single life-cycle stage and the control sample was an equal mixture of all four life-cycle stages. Hybridizations were performed as previously described . Briefly, equal volumes of control and experimental labeled cDNAs were combined with mouse CoT1-DNA (Life Technologies), poly(A)-DNA (Amersham Pharmacia Biotech), and hybridization buffer to yield final concentrations of 4× SSC, 40% formamide, 0.1% SDS, 0.5 μg/μl CoT1-DNA, 0.5 μg/μl poly(A)-DNA, and 25 mM Tris-Cl, pH 8.0). Samples were heat denatured at 95°C for 3 m, centrifuged at 12,000 × g at room temperature for 1 m, applied to the prehybridized array in a MAUI mixer chamber (BioMicro Systems, Inc., Salt Lake City, Utah), placed in a MAUI hybridization apparatus and hybridized at 42°C for 16–20 h with constant mixing. Following hybridization, the array, with attached mixer, was removed from the MAUI apparatus, the mixer was removed from the array, and the array was washed once for 2 m at room temperature in 300 ml PBS with 0.05% SDS, followed by 3 washes for 2 m each at room temperature in 300 ml PBS, and one wash for 1 m at room temperature in 300 ml 0.2× PBS. The slide was then dried by centrifugation.
Scanning and data analysis
Immediately following hybridization and washing, arrays were scanned on a ScanArray 4000 (GSI Lumonics, Wilmington, MA) at the Integrated Biotech Laboratories at the University of Georgia. Microarrays were scanned at multiple laser and photomultiplier tube (PMT) settings in each channel to obtain Tagged Image File Format (TIFF) images of matching sensitivity in the two channels. Microarray scans were quantified using The Institute for Genomic Research (TIGR) SpotFinder module of TM4 http://www.tm4.org. Individual sub-grids within each scan-specific grid were aligned manually. Background correction was used to flag spots with signal intensities less than local background intensity plus one standard deviation. The Otsu algorithm  was used to identify spots. The signal intensity files generated in SpotFinder were imported into the TIGR Microarray Data Analysis System (MIDAS) module of TM4 http://www.tm4.org using default parameters: e.g. MIDAS filtered spots with signal intensity in either channel less than 1, spots with signal intensities less than 2 × local background in either channel, and spots with integrated signal intensities less than 10,000. Since spots typically contained approximately 100 pixels, integrated intensities of 10,000 corresponded to mean signal intensities of 100. Normalization was performed by locally weighted least squares regression (LOWESS) within each sub-grid using the default smoothing parameter of 0.33. Within-slide replicate spots analysis was used to convert the signal intensities of one replicate spot in a pair of replicate spots to the geometric mean of it and its replicate and to set the signal intensities for the other spot to zero, thus compressing the data. The within-slide replicate spots analysis also had the effect of setting the signal intensities to zero for both spots in a replicate pair to zero if the signal intensity in either channel for either spot was zero. The normalized signal intensity files thus generated in MIDAS were paired by dye-swap and checked for dye-swap consistency, again using the default parameters. Spots with greater than 2 standard deviations difference in log2 ratios between dye-swap replicates were assigned zero intensities. Spots that were within threshold consistency were averaged (geometric mean). Thus, each set of two input dye-swap files resulted in one output file. The above described operations, including setting minimum intensity cut-offs, setting signal intensities to zero for spots with signal intensities not significantly above background, averaging of within-slide replicate spots, and setting signal intensities to zero for spots with significant dye bias were standard data manipulations for array data QC and the utilization of technical replication (replicate spots and dye-swap hybridizations) to eliminate the influence of low-quality spots and/or dye bias from the downstream data analyses , thus leaving only biological variation to be evaluated, which was the intent of this study. Example data analysis reports for the MIDAS analyses are in Additional file 6. These analyses generated 12 result files, one from each biological replicate for each of the four life-cycle stages. Each row in the final result files corresponded to eight independent measurements (2 channels × 2 replicate spots × 2 dye-swap hybridizations). Using the TIGR MultiExperiment Viewer (MeV) software the result files were filtered to remove genes with zero intensity in more than 3 result files (percentage cut-off of 75%), leaving 10,256 spots. The genes and experiments were median centered and normalized using the default method for this software. To find spots with significant, repeatable, between-stage variation we analyzed the result files by multi-class SAM for 4 groups (amastigotes, trypomastigotes, epimastigotes, metacyclic trypomastigotes – 3 files each). The settings for the SAM analysis were: number of permutations = 1,000 (10 times the default value), select S0 using Tusher et al. method  (default), calculate q values? = yes (not default – slow), imputation engine = K nearest neighbors (default), number of neighbors = 10 (default), and construct hierarchical trees? = yes (not default). Significant spots were selected with a median false discovery rate (FDR) of 0.06071% and a 90th percentile FDR of 0.40582%. The microarray data were deposited in the Gene Expression Omnibus http://www.ncbi.nlm.nih.gov/geo/ under the accession GSE14641.
Genes from functional groups of interest were selected for qRT-PCR analysis in a given life-cycle stage based on the gene being significantly regulated in the microarray analysis, as determined by SAM.
Total RNA was prepared from 3 independent samples from each life-cycle stage as described above. RNA samples were treated with DNase I (Promega) at a final concentration of 1 U/μg RNA at 37°C for 45 m. RNA was purified from the DNase digestions using RNeasy RNA purification columns (Qiagen). First strand cDNA reactions from 5 μg total RNA were primed with a mixture of oligo d(T) and a cDNA primer (Tc-18S cDNA: AAGAAATATCGGTGAACTTTCG) specific to the 3' end of T. cruzi 18S rRNA. Briefly, 0.5 μg oligo d(T)V and 5 pmoles T. cruzi 18S rRNA primer were added to the RNA samples and the primer/template mixtures were denatured at 65°C for 5 m. The samples were cooled to 42°C and final 20 μl reactions were assembled containing 1× first strand buffer, 10 mM DTT, 0.5 mM dNTP, 40 U RNasin, and 200 U SuperScript II reverse transcriptase. The cDNA reactions were incubated at 42°C for 1.5 h. For each RNA sample, duplicate control cDNA reactions were prepared in which the reverse transcriptase was omitted. Following first strand synthesis, the reactions were heat inactivated at 70°C for 10 m, then treated with 2 U RNase H (Invitrogen) at 37°C for 45 m.
Quantifications of selected T. cruzi genes were performed for stage and reference cDNA samples in duplicate on an iCyler (Bio-Rad Laboratories, Hercules, CA) with an iQ5 Multicolor Real-Time PCR Detection System (Bio-Rad). The primer sequences used in the qRT-PCR analyses are available in Additional file 7. Reactions were prepared containing 5 pmoles forward and reverse primers, 1× iQ SYBR® Green Supermix (Bio-Rad), and 2 μl template DNA. Standard curves were prepared for each run using known quantities of T. cruzi genomic DNA (ten-fold dilutions beginning at 15 ng/μl) and primers for the gene being quantified. The raw quantifications were calculated using the iQ5 Optical Detection System software and normalized to the 18S rRNA values for each sample. The final stage/reference ratios were the averages, for each gene, of the 9 possible normalized stage/reference comparisons (3 stage samples × 3 reference samples).
The oligonucleotide microarrays for this project were provided by the PFGRC. We are grateful to L. Cheng, C. Boehlke, and G. Cooley for excellent technical assistance. This work was supported by NIH grants R03AI064622 and P01AI044979.
- Krieger M, Freund A, Avila A, Muniz B, Probst C, Pavoni D, Sunaga D, Marchini F, Picchi G, Lenzi K, et al: Analysis of the gene expression program during the cellular differentiation of Trypanosoma cruzi (metacyclogenesis) through microarray hybridization. MPM XIII: 2002. 2002, Woods Hole, MA: Marine Biological Laboratory, 232B-Google Scholar
- Avila AR, Dallagiovanna B, Yamada-Ogatta SF, Monteiro-Goes V, Fragoso SP, Krieger MA, Goldenberg S: Stage-specific gene expression during Trypanosoma cruzi metacyclogenesis. Genet Mol Res. 2003, 2 (1): 159-168.PubMedGoogle Scholar
- Minning TA, Bua J, Garcia GA, McGraw RA, Tarleton RL: Microarray profiling of gene expression during trypomastigote to amastigote transition in Trypanosoma cruzi. Mol Biochem Parasitol. 2003, 131 (1): 55-64. 10.1016/S0166-6851(03)00189-0.View ArticlePubMedGoogle Scholar
- Baptista CS, Vencio RZ, Abdala S, Valadares MP, Martins C, de Braganca Pereira CA, Zingales B: DNA microarrays for comparative genomics and analysis of gene expression in Trypanosoma cruzi. Mol Biochem Parasitol. 2004, 138 (2): 183-194. 10.1016/j.molbiopara.2004.06.017.View ArticlePubMedGoogle Scholar
- Paba J, Ricart CA, Fontes W, Santana JM, Teixeira AR, Marchese J, Williamson B, Hunt T, Karger BL, Sousa MV: Proteomic analysis of Trypanosoma cruzi developmental stages using isotope-coded affinity tag reagents. Journal of proteome research. 2004, 3 (3): 517-524. 10.1021/pr034075o.View ArticlePubMedGoogle Scholar
- Paba J, Santana JM, Teixeira AR, Fontes W, Sousa MV, Ricart CA: Proteomic analysis of the human pathogen Trypanosoma cruzi. Proteomics. 2004, 4 (4): 1052-1059. 10.1002/pmic.200300637.View ArticlePubMedGoogle Scholar
- Parodi-Talice A, Duran R, Arrambide N, Prieto V, Pineyro MD, Pritsch O, Cayota A, Cervenansky C, Robello C: Proteome analysis of the causative agent of Chagas disease: Trypanosoma cruzi. Int J Parasitol. 2004, 34 (8): 881-886. 10.1016/j.ijpara.2004.05.002.View ArticlePubMedGoogle Scholar
- Atwood JA, Weatherly DB, Minning TA, Bundy B, Cavola C, Opperdoes FR, Orlando R, Tarleton RL: The Trypanosoma cruzi proteome. Science. 2005, 309 (5733): 473-476. 10.1126/science.1110289.View ArticlePubMedGoogle Scholar
- Atwood JA, Minning T, Ludolf F, Nuccio A, Weatherly DB, Alvarez-Manilla G, Tarleton R, Orlando R: Glycoproteomics of Trypanosoma cruzi trypomastigotes using subcellular fractionation, lectin affinity, and stable isotope labeling. Journal of proteome research. 2006, 5 (12): 3376-3384. 10.1021/pr060364b.View ArticlePubMedGoogle Scholar
- Parodi-Talice A, Monteiro-Goes V, Arrambide N, Avila AR, Duran R, Correa A, Dallagiovanna B, Cayota A, Krieger M, Goldenberg S, et al: Proteomic analysis of metacyclic trypomastigotes undergoing Trypanosoma cruzi metacyclogenesis. J Mass Spectrom. 2007, 42 (11): 1422-1432. 10.1002/jms.1267.View ArticlePubMedGoogle Scholar
- Koumandou VL, Natesan SK, Sergeenko T, Field MC: The trypanosome transcriptome is remodelled during differentiation but displays limited responsiveness within life stages. BMC genomics. 2008, 9: 298-10.1186/1471-2164-9-298.PubMed CentralView ArticlePubMedGoogle Scholar
- Clayton CE: Life without transcriptional control? From fly to man and back again. Embo J. 2002, 21 (8): 1881-1888. 10.1093/emboj/21.8.1881.PubMed CentralView ArticlePubMedGoogle Scholar
- Haile S, Papadopoulou B: Developmental regulation of gene expression in trypanosomatid parasitic protozoa. Curr Opin Microbiol. 2007, 10 (6): 569-577. 10.1016/j.mib.2007.10.001.View ArticlePubMedGoogle Scholar
- Campbell DA, Thomas S, Sturm NR: Transcription in kinetoplastid protozoa: why be normal?. Microbes Infect. 2003, 5 (13): 1231-1240. 10.1016/j.micinf.2003.09.005.View ArticlePubMedGoogle Scholar
- Jager AV, De Gaudenzi JG, Cassola A, D'Orso I, Frasch AC: mRNA maturation by two-step trans-splicing/polyadenylation processing in trypanosomes. Proc Natl Acad Sci USA. 2007, 104 (7): 2035-2042. 10.1073/pnas.0611125104.PubMed CentralView ArticlePubMedGoogle Scholar
- Nardelli SC, Avila AR, Freund A, Motta MC, Manhaes L, de Jesus TC, Schenkman S, Fragoso SP, Krieger MA, Goldenberg S, et al: Small-subunit rRNA processome proteins are translationally regulated during differentiation of Trypanosoma cruzi. Eukaryot Cell. 2007, 6 (2): 337-345. 10.1128/EC.00279-06.PubMed CentralView ArticlePubMedGoogle Scholar
- Holetz FB, Correa A, Avila AR, Nakamura CV, Krieger MA, Goldenberg S: Evidence of P-body-like structures in Trypanosoma cruzi. Biochem Biophys Res Commun. 2007, 356 (4): 1062-1067. 10.1016/j.bbrc.2007.03.104.View ArticlePubMedGoogle Scholar
- Cassola A, De Gaudenzi JG, Frasch AC: Recruitment of mRNAs to cytoplasmic ribonucleoprotein granules in trypanosomes. Mol Microbiol. 2007, 65 (3): 655-670. 10.1111/j.1365-2958.2007.05833.x.View ArticlePubMedGoogle Scholar
- McNicoll F, Muller M, Cloutier S, Boilard N, Rochette A, Dube M, Papadopoulou B: Distinct 3'-untranslated region elements regulate stage-specific mRNA accumulation and translation in Leishmania. J Biol Chem. 2005, 280 (42): 35238-35246. 10.1074/jbc.M507511200.View ArticlePubMedGoogle Scholar
- Horn D: Codon usage suggests that translational selection has a major impact on protein expression in trypanosomatids. BMC genomics. 2008, 9 (1): 2-10.1186/1471-2164-9-2.PubMed CentralView ArticlePubMedGoogle Scholar
- Duncan R: DNA microarray analysis of protozoan parasite gene expression: outcomes correlate with mechanisms of regulation. Trends Parasitol. 2004, 20 (5): 211-215. 10.1016/j.pt.2004.02.008.View ArticlePubMedGoogle Scholar
- Bringaud F, Muller M, Cerqueira GC, Smith M, Rochette A, El-Sayed NM, Papadopoulou B, Ghedin E: Members of a large retroposon family are determinants of post-transcriptional gene expression in Leishmania. PLoS pathogens. 2007, 3 (9): 1291-1307. 10.1371/journal.ppat.0030136.View ArticlePubMedGoogle Scholar
- Dallagiovanna B, Correa A, Probst CM, Holetz F, Smircich P, de Aguiar AM, Mansur F, da Silva CV, Mortara RA, Garat B, et al: Functional genomic characterization of mRNAs associated with TcPUF6, a pumilio-like protein from Trypanosoma cruzi. J Biol Chem. 2008, 283 (13): 8266-8273. 10.1074/jbc.M703097200.PubMed CentralView ArticlePubMedGoogle Scholar
- D'Orso I, Frasch AC: Functionally different AU- and G-rich cis-elements confer developmentally regulated mRNA stability in Trypanosoma cruzi by interaction with specific RNA-binding proteins. J Biol Chem. 2001, 276 (19): 15783-15793. 10.1074/jbc.M010959200.View ArticlePubMedGoogle Scholar
- D'Orso I, Frasch AC: TcUBP-1, a developmentally regulated U-rich RNA-binding protein involved in selective mRNA destabilization in trypanosomes. J Biol Chem. 2001, 276 (37): 34801-34809. 10.1074/jbc.M102120200.View ArticlePubMedGoogle Scholar
- Di Noia JM, D'Orso I, Sanchez DO, Frasch AC: AU-rich elements in the 3'-untranslated region of a new mucin-type gene family of Trypanosoma cruzi confers mRNA instability and modulates translation efficiency. J Biol Chem. 2000, 275 (14): 10218-10227. 10.1074/jbc.275.14.10218.View ArticlePubMedGoogle Scholar
- Robles A, Clayton C: Regulation of an amino acid transporter mRNA in Trypanosoma brucei. Mol Biochem Parasitol. 2008, 157 (1): 102-106. 10.1016/j.molbiopara.2007.09.005.View ArticlePubMedGoogle Scholar
- Noe G, De Gaudenzi JG, Frasch AC: Functionally related transcripts have common RNA motifs for specific RNA-binding proteins in trypanosomes. BMC molecular biology. 2008, 9 (1): 107-10.1186/1471-2199-9-107.PubMed CentralView ArticlePubMedGoogle Scholar
- Tusher VG, Tibshirani R, Chu G: Significance analysis of microarrays applied to the ionizing radiation response. Proc Natl Acad Sci USA. 2001, 98 (9): 5116-5121. 10.1073/pnas.091062498.PubMed CentralView ArticlePubMedGoogle Scholar
- Yoshida N: Molecular basis of mammalian cell invasion by Trypanosoma cruzi. An Acad Bras Cienc. 2006, 78 (1): 87-111.View ArticlePubMedGoogle Scholar
- Huang da W, Sherman BT, Lempicki RA: Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nature protocols. 2009, 4 (1): 44-57. 10.1038/nprot.2008.211.View ArticlePubMedGoogle Scholar
- Dennis G, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, Lempicki RA: DAVID: Database for Annotation, Visualization, and Integrated Discovery. Genome biology. 2003, 4 (5): P3-10.1186/gb-2003-4-5-p3.View ArticlePubMedGoogle Scholar
- Brandao A, Urmenyi T, Rondinelli E, Gonzalez A, de Miranda AB, Degrave W: Identification of transcribed sequences (ESTs) in the Trypanosoma cruzi genome project. Mem Inst Oswaldo Cruz. 1997, 92 (6): 863-866. 10.1590/S0074-02761997000600024.View ArticlePubMedGoogle Scholar
- Buscaglia CA, Campo VA, Frasch AC, Di Noia JM: Trypanosoma cruzi surface mucins: host-dependent coat diversity. Nature reviews. 2006, 4 (3): 229-236. 10.1038/nrmicro1351.PubMedGoogle Scholar
- Aebersold R, Mann M: Mass spectrometry-based proteomics. Nature. 2003, 422 (6928): 198-207. 10.1038/nature01511.View ArticlePubMedGoogle Scholar
- Josic D, Clifton JG: Mammalian plasma membrane proteomics. Proteomics. 2007, 7 (16): 3010-3029. 10.1002/pmic.200700139.View ArticlePubMedGoogle Scholar
- Teixeira SM, Russell DG, Kirchhoff LV, Donelson JE: A differentially expressed gene family encoding "amastin," a surface protein of Trypanosoma cruzi amastigotes. J Biol Chem. 1994, 269 (32): 20509-20516.PubMedGoogle Scholar
- Ohno S: Evolution by gene duplication. 1970, New York,: SpringerView ArticleGoogle Scholar
- Diehl S, Diehl F, El-Sayed NM, Clayton C, Hoheisel JD: Analysis of stage-specific gene expression in the bloodstream and the procyclic form of Trypanosoma brucei using a genomic DNA-microarray. Mol Biochem Parasitol. 2002, 123 (2): 115-123. 10.1016/S0166-6851(02)00138-X.View ArticlePubMedGoogle Scholar
- Luu VD, Brems S, Hoheisel JD, Burchmore R, Guilbride DL, Clayton C: Functional analysis of Trypanosoma brucei PUF1. Mol Biochem Parasitol. 2006, 150 (2): 340-349. 10.1016/j.molbiopara.2006.09.007.View ArticlePubMedGoogle Scholar
- Cohen-Freue G, Holzer TR, Forney JD, McMaster WR: Global gene expression in Leishmania. Int J Parasitol. 2007, 37 (10): 1077-1086. 10.1016/j.ijpara.2007.04.011.View ArticlePubMedGoogle Scholar
- Saxena A, Worthey EA, Yan S, Leland A, Stuart KD, Myler PJ: Evaluation of differential gene expression in Leishmania major Friedlin procyclics and metacyclics using DNA microarray analysis. Mol Biochem Parasitol. 2003, 129 (1): 103-114. 10.1016/S0166-6851(03)00100-2.View ArticlePubMedGoogle Scholar
- Rochette A, Raymond F, Ubeda JM, Smith M, Messier N, Boisvert S, Rigault P, Corbeil J, Ouellette M, Papadopoulou B: Genome-wide gene expression profiling analysis of Leishmania major and Leishmania infantum developmental stages reveals substantial differences between the two species. BMC genomics. 2008, 9: 255-10.1186/1471-2164-9-255.PubMed CentralView ArticlePubMedGoogle Scholar
- De Gaudenzi J, Frasch AC, Clayton C: RNA-binding domain proteins in Kinetoplastids: a comparative analysis. Eukaryot Cell. 2005, 4 (12): 2106-2114. 10.1128/EC.4.12.2106-2114.2005.PubMed CentralView ArticlePubMedGoogle Scholar
- Tarleton RL, Koller BH, Latour A, Postan M: Susceptibility of beta 2-microglobulin-deficient mice to Trypanosoma cruzi infection. Nature. 1992, 356 (6367): 338-340. 10.1038/356338a0.View ArticlePubMedGoogle Scholar
- Abrahamsohn IA: Cytokines in innate and acquired immunity to Trypanosoma cruzi infection. Braz J Med Biol Res. 1998, 31 (1): 117-121. 10.1590/S0100-879X1998000100015.View ArticlePubMedGoogle Scholar
- Fichera LE, Albareda MC, Laucella SA, Postan M: Intracellular growth of Trypanosoma cruzi in cardiac myocytes is inhibited by cytokine-induced nitric oxide release. Infect Immun. 2004, 72 (1): 359-363. 10.1128/IAI.72.1.359-363.2004.PubMed CentralView ArticlePubMedGoogle Scholar
- Li WH, Yang J, Gu X: Expression divergence between duplicate genes. Trends Genet. 2005, 21 (11): 602-607. 10.1016/j.tig.2005.08.006.View ArticlePubMedGoogle Scholar
- Piras R, Piras MM, Henriquez D: The effect of inhibitors of macromolecular biosynthesis on the in vitro infectivity and morphology of Trypanosoma cruzi trypomastigotes. Mol Biochem Parasitol. 1982, 6 (2): 83-92. 10.1016/0166-6851(82)90067-6.View ArticlePubMedGoogle Scholar
- Tomlinson S, Vandekerckhove F, Frevert U, Nussenzweig V: The induction of Trypanosoma cruzi trypomastigote to amastigote transformation by low pH. Parasitology. 1995, 110: 547-554. 10.1017/S0031182000065264.View ArticlePubMedGoogle Scholar
- Rondinelli E, Silva R, Carvalho JF, de Almeida Soares CM, de Carvalho EF, de Castro FT: Trypanosoma cruzi: an in vitro cycle of cell differentiation in axenic culture. Exp Parasitol. 1988, 66 (2): 197-204. 10.1016/0014-4894(88)90091-4.View ArticlePubMedGoogle Scholar
- Isola EL, Lammel EM, Gonzalez Cappa SM: Trypanosoma cruzi: differentiation after interaction of epimastigotes and Triatoma infestans intestinal homogenate. Exp Parasitol. 1986, 62 (3): 329-335. 10.1016/0014-4894(86)90039-1.View ArticlePubMedGoogle Scholar
- Sturm NR, Vargas NS, Westenberger SJ, Zingales B, Campbell DA: Evidence for multiple hybrid groups in Trypanosoma cruzi. Int J Parasitol. 2003, 33 (3): 269-279. 10.1016/S0020-7519(02)00264-3.View ArticlePubMedGoogle Scholar
- El-Sayed NM, Myler PJ, Bartholomeu DC, Nilsson D, Aggarwal G, Tran AN, Ghedin E, Worthey EA, Delcher AL, Blandin G, et al: The genome sequence of Trypanosoma cruzi, etiologic agent of Chagas disease. Science. 2005, 309 (5733): 409-415. 10.1126/science.1112631.View ArticlePubMedGoogle Scholar
- Obado SO, Taylor MC, Wilkinson SR, Bromley EV, Kelly JM: Functional mapping of a trypanosome centromere by chromosome fragmentation identifies a 16-kb GC-rich transcriptional "strand-switch" domain as a major feature. Genome research. 2005, 15 (1): 36-43. 10.1101/gr.2895105.PubMed CentralView ArticlePubMedGoogle Scholar
- Liao P, Chen T, Chung P: A fast algorithm for multilevel thresholding. Journal of Information Science and Engineering. 2001, 17: 713-727.Google Scholar
- Quackenbush J: Microarray data normalization and transformation. Nature genetics. 2002, 32 (Suppl): 496-501. 10.1038/ng1032.View ArticlePubMedGoogle Scholar
- Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215 (3): 403-410.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.