Skip to main content

The trypanosome transcriptome is remodelled during differentiation but displays limited responsiveness within life stages



Trypanosomatids utilise polycistronic transcription for production of the vast majority of protein-coding mRNAs, which operates in the absence of gene-specific promoters. Resolution of nascent transcripts by polyadenylation and trans-splicing, together with specific rates of mRNA turnover, serve to generate steady state transcript levels that can differ in abundance across several orders of magnitude and can be developmentally regulated. We used a targeted oligonucleotide microarray, representing the strongly developmentally-regulated T. brucei membrane trafficking system and ~10% of the Trypanosoma brucei genome, to investigate both between-stage, or differentiation-dependent, transcriptome changes and within-stage flexibility in response to various challenges.


6% of the gene cohort are developmentally regulated, including several small GTPases, SNAREs, vesicle coat factors and protein kinases both consistent with and extending previous data. Therefore substantial differentiation-dependent remodeling of the trypanosome transcriptome is associated with membrane transport. Both the microarray and qRT-PCR were then used to analyse transcriptome changes resulting from specific gene over-expression, knockdown, altered culture conditions and chemical stress. Firstly, manipulation of Rab5 expression results in co-ordinate changes to clathrin protein expression levels and endocytotic activity, but no detectable changes to steady-state mRNA levels, which indicates that the effect is mediated post-transcriptionally. Secondly, knockdown of clathrin or the variant surface glycoprotein failed to perturb transcription. Thirdly, exposure to dithiothreitol or tunicamycin revealed no evidence for a classical unfolded protein response, mediated in higher eukaryotes by transcriptional changes. Finally, altered serum levels invoked little transcriptome alteration beyond changes to expression of ESAG6/7, the transferrin receptor.


While trypanosomes regulate mRNA abundance to effect the major changes accompanying differentiation, a given differentiated state appears transcriptionally inflexible. The implications of the absence of a transcriptome response in trypanosomes for both virulence and models of life cycle progression are discussed.


Trypanosoma brucei is the causative agent of sleeping sickness in humans and N'gana in cattle, and has a major economic and morbidity impact across much of Africa [1]. The principal mechanism of immune evasion in the mammalian host is antigenic variation, sequential expression of immunologically distinct variant surface glycoproteins (VSGs) at the cell surface [2], but additional mechanisms, including manipulation of the host immune system and antibody clearance from the surface also participate [35]. The full impact of these mechanisms on virulence is not understood at this time. The mammalian part of the life cycle is complex, involving penetration of multiple tissue spaces in the mammalian host, encompassing the bloodstream, lymphatic system and central nervous systems. Despite morphological similarities between mammalian forms, it is unknown if a developmental stimulus is required for tissue tropism or if this results from stochastic events. Further, differentiation to the procyclic insect form on entering the tsetse fly vector is accompanied by massive cellular remodeling, including replacement of the VSG coat with procyclins, a family of acidic glycoproteins, activation of the mitochondrion, changes to cellular and organelle morphology, altered cell cycle checkpoints and attenuation of endocytic activity [610]. Multiple stages have been described during fly infection, with clear changes in morphology and surface antigen expression [7, 11, 12]. How these alterations are controlled at the transcriptional level remains unclear; for example each form could represent a distinct developmental stage, implying stimulus-driven differentiation, or result from transcriptome flexibility, allowing modulation of levels of certain transcripts in response to altered conditions.

Polycistronic transcription in trypanosomes [13] precludes promoter-driven control of transcription, while specific degradation appears to be the major mechanism underpinning regulation of mRNA steady state levels. Microarray hybridization and real time (RT) PCR methods can monitor steady-state RNA levels, and are independent of mechanisms controlling mRNA abundance. An earlier whole genome microarray study [14] identified ~2% of ORFs as developmentally regulated, but several developmentally regulated factors, including the clathrin heavy chain and Rab11 [10, 15, 16] were not detected in that analysis.

Membrane trafficking is a defining characteristic of eukaryotic cells, playing major roles in nutrient uptake, turnover, signalling, immune defence, apoptosis and many other processes. Mechanisms regulating transport remain only partly understood in any organism, but transcriptional, post-transcriptional and post-translational processes are all implicated, with changes in transcription [17], complex kinase integration [18], and formation of multiple protein-protein interactions [19, 20] contributing to control. Further, copious evidence implicates small GTPases in regulating distinct aspects of transport, including Rabs in vesicle fusion [21], ARFs in regulation of membrane coat systems [22] and Rho-related proteins in cytoskeletal interactions [23, 24]. Expression profiling for multiple human and mouse tissues indicates correlations between expression of certain Rabs and the SNAREs, adaptors, and other proteins with which they interact [17] and facilitates construction of potential interaction hubs and prediction of novel pathways.

Understanding of membrane transport is comparatively advanced in trypanosomes [25, 26], and various cellular changes are associated with life cycle progression [6], most notably an order of magnitude increase in endocytic activity in the mammalian infective bloodstream form (BSF) versus the insect procyclic (PCF) stages [10]. Increased endocytic activity in the BSF is likely related to immune evasion [4, 26, 27]. By contrast, knowledge of trypanosome signal transduction is poor [28]. Few Ras-like GTPases or their corresponding regulatory factors are present in trypanosomes, heterotrimeric GTPases are absent and there are no obvious receptor-type tyrosine kinases [2931]. There is a complex predicted trypanosome kinome [31] but limited understanding as to how signaling is mediated via these factors [3234]. There is evidence for phosphatidylinositol-mediated signaling, but this is comparatively unexplored [35], and a large novel family of receptor-coupled adenylate cyclases in kinetoplastids argues for significant environmental sensing via novel mechanisms [36, 37].

We selected the trypanosome trafficking system as a clear example of a developmentally regulated process, where sufficient data are available for meaningful investigation of transcriptome changes. Using microarray analysis and quantitative RT-PCR, we assessed developmental regulation of trypanosome membrane transport and transcriptional flexibility by targeted genetic manipulation and altered environmental conditions designed to modulate trafficking pathways. From work in other organisms, we anticipated the involvement of protein kinases, small GTPases and, in the case of the unfolded protein response (UPR), a system that responds to increased concentrations of non-native polypeptides within the endoplasmic reticulum, clear evidence for a transcriptional mechanism. Our data indicates developmental remodeling of the transcriptome, but essentially no response to altered conditions.


Microarray design

Due to limitations of ORF assignment for specific spots of the previously described whole genome microarray platform [14], we designed a more targeted and fully validated array to address developmental regulation of the trafficking system in T. brucei. Major classes of proteins involved in membrane transport include small GTPases, vesicle coats proteins, kinases and phosphatases, motility factors mediating movement of vesicles along cytoskeletal elements, tethering factors for binding of vesicles to target membranes, SNARE proteins, ESCRTs and numerous other molecules involved in protein sorting and membrane deformation. The trypanosome genome was parsed extensively for such factors [3841] and, based on these and additional searches, a subgenome microarray was constructed. The array also included putative ER chaperones, peptidases bearing an ER signal peptide and potentially involved in degradation within the trafficking system, homologues of kinases implicated in membrane traffic in other systems [18], sphingolipid biosynthesis enzymes with roles in exocytosis of glycosylphosphatidylinositol (GPI)-anchored proteins, as well as trypanosome-specific surface antigens, e.g. BARP, CRAM, and the invariant surface glycoprotein family, or ISGs [4244]. In addition we also included oligonucleotides against ESAG6 and ESAG7, which encode the trypanosome transferrin receptor, a molecule reported to exhibit expression that is modulated by growth conditions [45].

673 ORFs in the T. brucei genome with a potential role in membrane trafficking were identified, representing ~8% of the predicted protein coding potential. Representatives of most components of the major gene families are found in the T. brucei genome, but notable exceptions are the absence of several coat systems [46], a relatively small number of GTPases and associated factors [30], and the absence of BAR and SH2 domains; BAR domains are used extensively for membrane deformation in higher eukaryotes [47], while the absence of SH2 domains may indicate limited or unconventional signalling through protein tyrosine kinase pathways [40, 48]. Apart from clustering of multicopy genes, there is no location bias for a specific chromosome or region within a chromosome for these ORFs [see Additional file 1]. An additional 123 ORFs were also included as a reference. The first array (v1.0), used for comparison of life stages and analysis of Rab5 isoform over-expression, comprised 694 oligonucleotides, 56 of which targeted reference ORFs. A modified array (v1.1), used for all other analyses, included an additional 49 oligonucleotides. The vast majority of oligonucelotides are gene-specific but some are expected to hybridise to multiple transcripts from multicopy genes [see Additional file 2].

Overview of differential expression between trypanosome life-stages

We extracted RNA from exponentially growing cultures of T. brucei bloodstream and procyclic forms, labelled and competitively hybridised the cDNA to the microarray. As expected, the majority of genes were equivalently expressed in the two life stages, but a small cohort did present significant differential expression (Figure 1). The heat-map (Figure 1B) indicates that the data for the most significantly developmentally regulated genes are highly reproducible across all replicates (8 arrays, 4 replicate spots per array).

Figure 1

Overview of developmental transcriptome changes in the membrane trafficking system of trypanosomes. Panel A: Scatter plot of raw data for all 3600 spots on a representative microarray used for developmental expression experiments. Cy5 fluorescence is plotted on the Y-axis (bloodstream, BSF) and Cy3 fluorescence on the X-axis (procyclic, PCF). Spots with a BSF/PCF ratio above two are highlighted in red, while spots with a PCF/BSF ratio above two are highlighted in green. Panel B: Clustering of the data for eight microarray experiments comparing BSF to PCF, representing four biological replicates plus relevant dye-swaps. The scale indicates the colour scheme for the z-score of the data, i.e. how far and in what direction, the ratio for each spot deviates from the mean for each array, expressed in units of standard deviation; bright red indicates significant upregulation in BSF, bright green indicates significant upregulation in PCF, dark colours or black indicate no differential expression between the two developmental stages. For each target gene, the four replicate spots on the array were averaged. White indicates data points rejected because of too much variability in the replicate spots. A small but significant number of ORFs exhibit strong developmental regulation.

To identify which mRNAs are differentially expressed, limma and fspma packages were used to interrogate the data [4951]. For comparison with previous work [14], ratios above 2.5 (+/- 1.325, log2 scale) were considered highly significant. Using this threshold, 8% of the genes represented on the array are upregulated in BSF and 5% in PCF, but excluding the reference set these numbers reduce to 4% in BSF and 2% in PCF. Using a whole genome microarray, Brems and cowrokers found only 2% of trypanosome genes were developmentally regulated, and while different pre-processing steps, in both mRNA preparation and data analysis, make direct comparison difficult, the present data suggest that the trafficking cohort likely exhibits bias in developmental regulation. There is also greater upregulation in BSFs compared to PCFs, consistent with increased activity in the trafficking system in bloodstream stages [10]. Brems et al. found 80 genes upregulated in PCF and 34 in BSF while, using the same threshold, the present analysis identifies 16 upregulated in PCF and 34 in BSF, from a much smaller cohort. These data suggest that significant transcriptional level remodelling of the membrane trafficking system accompanies life-cycle progression in trypanosomes.

Strongly developmentally regulated genes include several reference ORFs with known regulated expression, i.e mRNAs for procyclin, trans-sialidase, CRAM, CAP17, and CAP5.5 are increased in PCF, while ESAG6, ESAG7, ISG65, ISG75, GPI-PLC, and CAP15 are upregulated in BSF (Table 1). These data are fully consistent with previous work and support the validity of the microarray approach for identifying previously uncharacterised developmentally expressed genes. Rab11, also developmentally regulated [15], gave an expression ratio below 2.5 by both limma and fspma analysis. Hence, analysis using a more relaxed threshold to identify developmentally regulated mRNAs was also used. Table 1 lists all ORFs which either gave a ratio above 2.0 (+/- 1.00, log2 scale) by limma or fspma analysis, or which had a B-statistic value greater than 1.0 in limma, indicating statistical significance in differential expression, despite a low relative expression ratio. Excluding the reference set, 101 genes with a possible role in membrane transport are differentially expressed in T. brucei; 64 upregulated in BSF, 37 in PCF, i.e. 8% and 5% respectively.

Table 1 Transcripts differentially expressed between the procyclic and bloodstream forms of T. brucei.

25% of the most significantly developmentally regulated genes were validated by qRT-PCR analysis, using the same mRNA samples as used for the microarray analysis. The vast majority of ORFs tested exhibited congruence between qRT-PCR and microarray analysis, i.e. differential regulation in the same direction (BSF or PCF, Table 1). Significantly, all ORFs with a ratio above 2.5 on the array also gave a significant ratio (i.e. above 2) in qRT-PCR. Two genes with known developmental regulation, the clathrin heavy chain [10, 16] and TbMSP-A [52], did not demonstrate significant differential regulation by microarray analysis, but highly significant ratios were obtained by qRT-PCR (data not shown). Various factors may be responsible for this array artefact, which were not pursued further; however these data specifically highlight examples of false-negative data in microarray analysis.

Differential expression of trafficking factors in trypanosomes

Considering the differentially expressed genes by functional class indicates that a broad remodelling accompanies differentiation, i.e. developmental changes are not restricted to specific protein classes or pathways. Specifically, members of the chaperones, some proteases and non-PI kinases are upregulated in PCF, while the small GTPases, additional proteases, SNARE proteins and PI kinases represent the major factors upregulated in BSFs (Figure 2, Additional file 6). An increased requirement for chaperones may be due to the greater heterogeneity of the procyclic cell surface compared to the bloodstream form, the latter being dominated by VSG. A differential requirement for secretory pathway proteases is consistent with distinct roles for these molecules between the two life stages, while upregulation of both small GTPases and SNARE proteins in BSFs is broadly consistent with increased activity within the trafficking system in the mammalian stage.

Figure 2

Significantly developmentally expressed trypanosome genes grouped by functional class. Genes represented by array oligonucleotides were grouped by function, based on sequence similarity to annotated sequences, domain annotation, GO terms and additional criteria. The proportion in each functional class significantly upregulated in BSF are shown in red and in PCF in green. The remaining fraction of genes in each class are shown in grey. The first bar illustrates the proportion of developmentally regulated transcripts in the entire studied gene cohort. For raw data [see Additional file 6].

Specificity and control factors

The most abundantly expressed small GTPases at the mRNA level are Ran, ARL3A, both Rab1A and Rab1B, and several members of the Arf subfamily (Figure 3). Ran is a very abundant protein in most cells, while high expression of Arf1, Rab1, and ARL3A possibly suggest a considerable emphasis on flux through the ER and Golgi complex and flagellum biogenesis respectively [5355]. Developmentally expressed small GTPases are restricted to the BSF. This cohort includes markers for the recycling endosome, Rab11 [15], the endosomes/Golgi, TbArf1 [53], and the flagellum, TbArl3A [55]. Identification of Rab11 is consistent with earlier work [15] and our data provide the first evidence for developmental regulation of TbArf1 and TbArl3A. The data also identify Rab1B, RabX3, Ras-like GTPases TbRX3 and TbRHP, and a putative GTPase with an Hsr1 domain (Tb927.4.2380) as developmentally regulated; with no experimental information for these proteins the significance of these changes is not clear [29, 56]. Importantly, we see no strong developmental regulation for any of the putative GTPase-activating proteins (GAPs) and guanine nucleotide exchange factors (GEFs) included in the study. Some developmental regulation of TBC family Rab-GAPs is detectable by qRT-PCR, but at levels unlikely to be detected by the array (C. Gabernet-Castello and MCF, unpublished data).

Figure 3

Relative steady state mRNA levels of selected trypanosome genes. Fluorescence intensity data for eight microarray experiments comparing BSF to PCF, representing four biological replicates plus dye-swaps, was normalised to the total signal intensity of each slide, averaged and plotted as bar-graphs against ORF designation or accession number. Error bars indicate the standard deviation, dark bars are BSF and light bars PCF. Panel A: Small GTPases. ARL and Rab annotations are based on [39], Ras-like annotations based on [29] and additional data at GeneDB. Panel B: Cell surface proteins. Panel C: SNARE family proteins. Annotations based on domain architecture and similarity to sequences described for L. major [59]. Panel D: Putative secretory pathway proteases. ORFs were included for predicted proteases that bear an ER targeting sequence and/or a predicted trans-membrane domain. In each cohort, a limited number of transcripts are highly expressed.

In addition to a number of small GTPases, including Rab subfamily members, five predicted SNARE mRNAs exhibit increased expression in BSFs (Table 1). While no direct experimental information is available for these molecules, Tb09.211.3920, the most abundant Qa SNARE at the mRNA level (Figure 3), is likely involved in post-Golgi transport on account of similarity to trans-Golgi network (TGN) SNAREs SYP42 from Arabidopsis thaliana and syntaxin16/Tlg2p from S. cerevisiae [57, 58], while the L. major orthologue, LmjF35.2720, localises close to the Golgi complex [59]. Tb10.70.7410 and Tb09.160.2420 are similar to the A. thaliana early endosomal SNARE VAMP727 and the ER-located SYP72, respectively [58], while Tb927.7.6440 and Tb927.3.3720 do not display significant similarity to any specific SNARE. By contrast, Tb10.61.0870, is the only upregulated SNARE message in the procyclic stage, with weak similarity to Aspergillus Bos1, suggesting a putative role in anterior transport from the ER; Tb10.61.0870 is also the most highly expressed Qb SNARE (Figure 3). These data are suggestive of specific upregulation of mRNAs encoding endocytic and exocytic SNAREs in BSF, with evidence for limited PCF regulation of the secretory pathway. On average, Qc and R SNAREs are more highly expressed than Qa and Qb SNAREs (Figure 3). The T. brucei N-ethylmaleimide sensitive factor (NSF, Tb927.1.1560), an ATPase responsible for regulation of SNARE complex disassembly [60], and a hypothetical protein similar to γ-soluble NSF attachment protein (γ-SNAP, Tb10.70.4620) are upregulated in BSF, which may indicate increased overall SNARE activity, but a further putative SNAP (Tb927.1.3110) is upregulated in PCF. Interestingly, both Tb927.1.1560 (NSF) and Tb927.1.3110 (SNAP) have severe RNAi knockdown growth phenotypes [61]. Expression of Vps16, a component of the HOPS tethering complex functioning at the late endosome, is increased in the BSF (Table 1). Overall, multiple factors likely to be controlling endocytic pathways, encompassing SNAREs, Rabs and possibly also the HOPS complex, appear more highly expressed in the BSF, consistent with earlier data but also suggesting more extensive remodelling than previously suspected.

Coat proteins

The vast majority of vesicle coat factor mRNAs are expressed at equivalent levels in both major life stages. Of those that are regulated, most notable is the clathrin light chain, Tb10.61.1910, upregulated in BSF and in agreement with increased endocytosis in this life stage [16]. Further, Tb11.01.6880, containing a predicted emp24 domain, suggesting a role as a cargo receptor in COPII-mediated transport [62], is upregulated in BSF, whereas a second emp24-related ORF, Tb927.4.4350, is upregulated in PCF. Should these proteins indeed represent emp24 isoforms, this expression profile indicates potential developmental-stage specific remodelling of COPII cargo selection providing differentiation in ER exit pathways.

A putative epsinR (Tb11.50.0006), is also significantly upregulated in BSF; epsinR is associated with clathrin coated vesicles and AP1 [63]. The AP1 subunit is significantly upregulated in BSF, but the μ-subunit of AP3, is upregulated in PCF (Table 1); both complexes have a role in endosome to TGN transport. Although no further adaptin components exhibit highly significant differential expression, a general trend of upregulation of AP1 component in BSFs and downregulation of AP3 and AP4 was evident on the microarray [see Additional file 3] which suggests developmental regulation. From previous analyses, the AP1 γ-subunit is constitutively expressed [64], but the β-chain is upregulated in BSF [16], suggesting possible subunit-specific regulation. To further confirm the significance of these data polyclonal antibodies were raised against the δ-subunit of AP3. Western blotting confirmed significant developmental regulation of AP3δ, with expression greatly augmented in the PCF [see Additional file 3]. These data are also consistent with evidence suggesting increased emphasis on lysosomal trafficking in the PCF [26, 65].


Several proteases are known components of the trypanosome surface and endocytic system or predicted to be, based on the presence of a signal peptide. The gp63 isoform TbMSP-B (Tb927.8.1610), is upregulated in PCF, while TbMSP-C (Tb10.70.5290) is upregulated in BSF (Table 1), confirming previous analyses [52]. TbMSP-A, previously identified as BSF-specific [52], is only slightly upregulated by microarray, but is validated strongly by qRT-PCR (results not shown). Further, metacaspase MCA3 (Tb927.6.930) showed significant upregulation in BSF, also confirming earlier data [66]. Microarray data indicate that MCA5 (Tb09.211.4760) is upregulated in PCF, but qRT-PCR data failed to validate this, suggesting more equivalent expression, also in agreement with earlier work [66]. Two subtilysin-like serine peptidases (Tb927.3.4230, Tb11.02.1280), a CAAX prenyl protease (Tb09.211.0680), and a thimet oligopeptidase A (Tb927.7.190) are upregulated in BSF. The thimet oligopeptidase Tb927.7.190 was recently suggested to play a role in dysregulated kinin metabolism observed in the plasma of trypanosome-infected hosts and to contribute to vascular lesions observed in African trypanosomiasis [67], consistent with augmented expression in BSF.

The most abundant protease class message is for Tb927.6.960, an L-cathepsin-like cysteine peptidase, and likely a lysosomal protein [68], which appears equivalently expressed in both life stages (Figure 3). Some of the developmentally regulated proteases are also abundant at the mRNA level (Figure 3), namely MSP-B (Tb927.8.1610), MCA5 (Tb09.211.4760), a subtilysin-like serine peptidase (Tb927.3.4230), and thimet oligopeptidase A (Tb927.7.190). Hence the major proteolytic activity likely resides within the cysteine protease class, as previously reported [69], but these new data suggest a significant level of developmental remodelling, as well as potential contributions from additional protease classes.

Chaperones and ER lumenal proteins

Several putative ER chaperones, likely to aid folding of secreted coat proteins, were also studied. Of these, Tb927.7.5790 is homologous to ERp44, a thioredoxin of the protein disulfide isomerase (PDI) family, which may have a role in controlling the redox state of the ER, and is upregulated in BSF [70]. Two additional PDI proteins are also upregulated in the bloodstream stage of T. brucei [71, 72], suggesting remodelling of these factors between life stages.

In the PCF stage, two predicted mitochondrial chaperones, an Hsp70 (Tb927.6.3740) and an Hsp60 (Tb10.70.0280), are upregulated, consistent both with previous data [73, 74] and increased mitochondrial activity [6]. Additionally, an Hsp70 (Tb11.01.3110) with high similarity to human Hsc70, an uncoating ATPase for clathrin coated vesicles, is upregulated in PCF [7577]. Since clathrin-mediated endocytosis is upregulated in BSF but lysosomal targeting upregulated in PCF, these data may indicate a role in lysosomal targeting, but there are multiple Hsp70 ORFs in T. brucei and Hsp70 has been previously reported to be upregulated slightly in BSF [78]. Two novel ER lumenal proteins, Tb11.01.2640 and Tb11.01.8120, with a predicted signal peptide and C-terminal HDEL ER-retrieval motif are developmentally regulated, the former with augmented expression in BSF and the latter in PCF. Both are localised to the ER, and give a severe growth phenotype by RNAi knockdown, suggesting an important function within the ER (T. Sergeenko and MCF, unpublished data). Overall, the expression profile of the chaperones does not provide evidence for substantial differential remodelling of the ER folding environment, and may suggest a comparatively similar set of factors are required for synthesis of the cell surface in PCF and BSF, despite major differences in the structures of the superabundant surface antigens.

Lipid kinases

Two highly conserved lipid kinases with roles in intracellular trafficking are known in trypanosomes: Fab1 (Tb11.47.0002), a PI 5-kinase and TbVps34 (Tb927.8.6210), a PI 3-kinase involved in endocytosis and Golgi segregation [35, 46]. Neither of these factors shows evidence for differential expression, but several other predicted phosphatidylinositol kinases are developmentally regulated (Table 1). A type III phosphatidylinositol 4-kinase, Tb927.3.4020, highly similar to S. cerevisiae Stt4p, is upregulated in PCF, consistent with a recent report of a role for a type III PI 4-kinase in protein trafficking, maintenance of Golgi and cell morphology and cytokinesis in PCFs [79]. A predicted phosphatidylinositol 3-kinase, Tb10.6k15.2060, is also upregulated in PCF and most likely represents an atypical protein kinase of the TOR-related family [35]. Significantly, human FRAP1, a member of the TOR-related kinase family, is required for clathrin-mediated endocytosis [18]. Tb11.01.6980, an atypical PI 3-kinase with a PI(3)P-binding FYVE domain and a putative N-terminal trans-membrane helix is highly upregulated in BSF; on account of the FYVE domain this protein is expected to localise to the endosomal system [35]. Further, Tb10.70.2440, a predicted PI4P 5-kinase, is also upregulated in BSF, while two serine threonine kinases with pleckstrin homology (PH) domains, capable of recognising PI4,5P lipids, are upregulated in BSF. Taken together these data suggest significant remodelling of the PI-based signalling system between life stages, together with evidence for increased expression of several putative endosomal functioning PI-kinases in the BSF.

Protein kinases

A genome-wide analysis of the roles of kinases in clathrin- and caveolin-mediated endocytosis was performed recently in H. sapiens [18]; significantly, caveolin-mediated endocytosis is restricted to the metazoa [46]. We looked for orthologues of the kinases that gave an endocytosis phenotype in the Pelkmans et al . study [18] and included these in the array to analyse their differential expression.

Two metabolic pathway kinases and three putative signalling kinases are upregulated in BSF. These include Tb10.70.5800, a glycosomal hexokinase, expected due to increased glycolytic activity in BSFs [14, 80]; this kinase is similar to H. sapiens hexokinase II, required for clathrin-mediated endocytosis [18]. Pyruvate kinase (Tb10.61.2680) is also upregulated, in agreement with earlier work [81], and interestingly the related H. sapiens kinase, PKLR, inhibits caveolae-mediated endocytosis [18]. The three signalling kinases upregulated in BSF all show similarities to endocytosis kinases [18]. Tb927.6.1780, a putative MAP-kinase similar to human MAPK8, is activated in response to temperature stress [34]. A MAP-kinase kinase kinase, Tb10.61.2490, has similarity to human MAP3K8, and finally a calcium-dependent protein kinase, Tb927.4.3770, is highly similar to CAMK4, which participates in a calcium/calmodulin-dependent protein kinase cascade [82]. These data clearly suggest a role for kinase activity in BSF trafficking.

Eight primarily metabolic pathway kinases are upregulated in PCF. One is a nucleoside diphosphate kinase (Tb11.01.7800) which is involved in the equilibration of cellular pools of nucleotide triphosphates, and which binds telomeres and acts as a Rho-GAP in other systems [8385]. This protein localises predominantly to the nucleus in T. brucei [86] and shows similarity to human NME6, found to inhibit clathrin-mediated endocytosis [18]. A ribokinase (Tb11.03.0090) upregulated in PCF is highly similar to human RBSK, required for clathrin-mediated endocytosis [18]. A putative NEK kinase (Tb927.3.3190) is related to human NEK2, 6, and 7, all of which are required for clathrin-mediated endocytosis [18], and a putative UTP-glucose-1-phosphate uridylyltransferase 2 (Tb10.389.0330) is related to human UGP2, required for clathrin-mediated endocytosis [18]. A glycosomal glycerol kinase (Tb09.211.3540) related to human CARKL, which is required for caveolin/lipid raft-mediated endocytosis [18], was upregulated in PCF, although equivalent levels of specific activity for this enzyme have been reported for both life stages [80, 87]. An arginine kinase (Tb09.160.4560) involved in the management of cellular ATP energy reserves [88] is similar to human CKB, required for clathrin-mediated endocytosis [18]. A putative dephospho-CoA kinase (Tb927.6.710) which catalyzes the final step in CoA biosynthesis, is related to human DPCK, which is required for caveolin/lipid raft-mediated endocytosis [18]. A phosphoglycerate kinase (Tb927.1.720), similar to human PGK1, which is involved in both clathrin-mediated and caveolin/lipid raft-mediated endocytosis [18] was up in PCF, consistent with previous analyses for T. brucei PGKB [89], although the oligonucleotide used does not allow us to differentiate PGK isoforms. A putative Ser/Thr and Tyr dual specificity phosphatase (Tb10.406.0240) possibly involved in MAP kinase signalling, is also upregulated in PCF. Taken together these data suggest multiple kinases are involved in control of trafficking in both life stages, with more signalling kinases upregulated in BSF, and these data provide a number of candidates for further analysis.


Covalent attachment of ubiquitin to surface molecules is an important mechanism for internalization and lysosomal targeting [90, 91], and is present in trypanosomes [92, 93]. Five ORFs likely to participate in the trypanosome ubiquitylation system are upregulated in BSF; Tb09.211.3610, an ubiquitin-activating E1 enzyme, Tb927.4.2710, an ubiquitin-conjugating E2 enzyme, Tb927.6.2370, an ubiqutin ligase E3, Tb11.46.0014, an ubiquitin carboxyl-terminal hydrolase similar to S. cerevisiae Doa4p protein, and Tb10.70.1130, which contains two putative ubiquitin interaction motifs, but otherwise has no predicted function. Upregulation of a near complete pathway for protein ubiquitylation and deubiquitylation in BSFs is suggestive of stage-specific ubiquitylation, but clearly direct evidence is required [94]. RNAi of the E1 Tb09.211.3610 shows a strong growth defect but does not affect the stability of ISG65, a surface protein subject to ubiquitylation [93]. A protein containing a VHS domain (Tb11.01.2540) which is normally found in the Vps27 component of ESCRT 0, and a putative Vps46 component of the ESCRT III-associated complex (Tb11.02.0580) are upregulated in procyclics, further supporting differential regulation of the ubiquitin sorting system.

Hypothetical open reading frames

A considerable fraction of the analysed ORFs are annotated as hypothetical, and were included on account of sequence features shared with characterized trafficking factors. Several of these exhibit developmental regulation. For example, the C2 domain is a Ca2+-dependent membrane-targeting module [9597]; two ORFs with predicted C2 domains, Tb09.160.3240 and Tb927.7.3550, are upregulated in procyclics, and one, Tb927.4.2080, is upregulated in BSF. Interestingly, Tb927.7.3550 localises between the inner face of the plasma membrane and the sub-pellicular corset of microtubules in T. brucei, and includes a highly charged region characteristic of a tubulin-binding domain [98]. RNAi evidence indicates a clear role for Tb927.7.3550 in cytoskeletal function [98].

Other hypothetical genes that show significant differential regulation at the mRNA level are, in PCF: Tb927.7.2640, which contains part of a Sec34 domain (although it is not a true COG3/Sec34 orthologue [41]), and Tb11.47.0022, which contains a signal peptide and a trans-membrane helix; both ORFs are highly conserved between trypanosomatids but have no identifiable orthologues in other organisms. In BSF: Tb11.01.4460, which contains a cysteine-rich repeat region found in the low-density lipoprotein (LDL) receptor, and Tb927.5.4520, an extended coiled-coil protein are upregulated. The flagellum adhesion glycoprotein fla1 was upregulated in PCF, as has been reported previously [99], but the closely related fla2 [100] was upregulated in BSF, indicating that the two gene products may have similar functions but distinct expression patterns.

Surface antigens

Messages for BSF-specific ESAG6/7 and ISG65/75 surface antigens are highly abundant (Figure 3). By contrast, CRAM, a possible lipoprotein receptor abundant in PCF, and BARP, a putative glycolipid-lipid raft associated protein are not so highly expressed (Figure 3). Interestingly, some of the oligonucleotides for BARP showed upregulation in PCF and some in BSF (Table 1); BSF-specific and epimastigote-specific BARP expression have been reported previously [12, 44]. The differential expression levels for the ISG75 mRNAs (Figure 3) represent mostly the degeneracy of each oligonucleotide [see Additional file 2], but the mRNAs for the ISG65 genes on chromosome 5 (Tb927.5.1390-1410-1430) are clearly more abundant than those for the ISG65 genes on chromosome 2 (Tb927.2.3280-3290-3300-3310). We also included oligonucleotides for several ISG-like proteins, four of which, Tb927.5.310, Tb927.5.420, Tb927.5.430, and Tb927.5.630, are significantly upregulated in BSF, suggesting additional and so far uncharacterised developmental modification of the parasite surface.

Influence of stress and growth conditions on the transcriptome

The detection of a significant cohort of developmentally regulated trypanosome mRNAs by microarray, and their validation by qRT-PCR, provided confidence that our procedure could be exploited for investigation of transcriptome changes associated with differing conditions. In several model organisms, the interactions between gene products has been probed by analysis of coordinated changes to the transcriptome in response to challenges to the cell. Such studies have included several changes to growth conditions, insult with specific chemicals to invoke the unfolded protein response (UPR), and targeted manipulations of gene expression [101104]. Here we used all three approaches to probe the T. brucei transcriptome in order to determine the level of flexibility, and to assess if such data could facilitate construction of interaction and pathway maps based on coordinate expression profiles.

Iron starvation

Bloodstream trypanosomes are sensitive to iron starvation, and upregulate ESAG6/7 transferrin receptor (TfR) mRNA and protein levels when placed in low transferrin or iron-depleted media [45, 105] and adjust TfR expression to compensate for reduced endocytosis [65]. Both observations suggest a specific transcriptional-level response resulting from an iron-sensing mechanism. The transcriptomes of BSFs cultured in different proportions of fetal bovine serum (FBS), in the absence of FBS or in the absence of FBS but supplemented with bovine holo-transferrin were analysed. An increase in ESAG6/7 mRNA was detected under FBS depletion by the array (Table 2), which was validated by qRT-PCR in agreement with previous analyses [45, 105]. mRNAs for procyclin, the gp63 isoform TbMSP-C (Tb10.70.5290) and the major cysteine peptidase (Tb927.6.960) were upregulated in serum-depleted cells (Table 2). In cells grown in 30% FBS, procyclin mRNA, as well as an ISG-like protein (Tb927.5.410) and a putative calcium-dependent protein kinase (Tb927.4.3770) were downregulated (Table 2), which contrasts with a moderate upregulation of procyclin in serum free conditions. These alterations were confirmed by qRT-PCR, and may suggest a suppression of procyclin expression by serum factors. However, overall there were no additional transcriptome changes when cells were grown in 30% FBS or in 0% FBS with or without transferrin. Scatter plots of the entire array data set indicate few transcripts that fall outside of the region considered as equivalent, in marked contrast to the scattergram obtained for developmental changes (compare Figure 4B–D with Figure 4A or Figure 1). These data suggest a very limited and specific ability of the BSF to respond to alterations in levels of nutrients, which may be restricted primarily to ESAG6/7.

Table 2 Expression changes accompanying environmental challenge in T. brucei.
Figure 4

Scatter plots of transcriptome data from trypanosomes subjected to altered culturing conditions and other challenges. In each case, Cy5 fluorescence is plotted on the Y-axis and Cy3 fluorescence on the X-axis. Spots with a Cy5/Cy3 or Cy3/Cy5 ratio above two, i.e. significantly differentially expressed are highlighted in black, spots with a ratio below two are shown in grey. Panel A: BSF (Cy5) vs PCF (Cy3); data are identical to Figure 1 and reproduced for comparison. Panel B: BSF in 30% FBS (Cy5) vs BSF (Cy3). Panel C: BSF (Cy5) vs BSF in 0% FBS, 5 mg/ml BSA (Cy3). Panel D: BSF (Cy5) vs BSF in 0% FBS, 5 mg/ml, BSA 0.3 mg/ml Tfn (Cy3). Panel E: BSF in 1 mM DTT for 1 hr (Cy5) vs BSF (Cy3). Panel F: BSF in 1 mM DTT for 4 hr (Cy5) vs BSF (Cy3). Panel G: BSF in 5 μg/ml tunicamycin for 4 hr (Cy5) vs BSF (Cy3). Panel H: BSF in 5 μg/ml tunicamycin for 24 hr (Cy5) vs BSF (Cy3). Panel I: VSG RNAi in BSF, induced for 24 hr (Cy5) vs uninduced (Cy3). Panel J: VSG RNAi in BSF, induced for 72 hr (Cy5) vs uninduced (Cy3). Panel K: CLH RNAi in BSF, induced for 24 hr (Cy5) vs uninduced (Cy3). Panel L: BSF (Cy5) vs MITat1.1 grown in vivo in rats (Cy3). The scatter plots for most experiments (panels B-K) indicate few transcripts that fall outside of the region considered as constitutive, in marked contrast to the scattergrams obtained for developmental changes (panel A) and in the comparison between in vitro versus in vivo BSF cultures (panel L).

Unfolded protein response (UPR)

Transcriptome responsiveness was further tested by attempts to invoke the UPR, a classic ER-based pathway stimulated by the presence of increased levels of unfolded polypeptides [106]. The UPR is mediated via transcriptional responses in both yeast and metazoan cells [104, 107], and can be reliably activated by addition of tunicamycin or dithiothreitol (DTT). Using tunicamycin and DTT at concentrations that invoke a UPR in mammalian cells, Arabidopsis, yeast and other systems [108112], we found that, in T. brucei, tunicamycin resulted in growth arrest over a period of up to 24 hours, while DTT exposure led to rapid cell death [see Additional file 4]. This concentration of tunicamycin also efficiently inhibits trypanosome N-glycosylation [113]. Therefore, we analysed the transcriptome at 5 μg/ml and 1 mM for tunicamycin and DTT, respectively, under conditions where cells remained viable, as assessed by motility.

Essentially no significant change to the transcriptome was observed by either treatment. Very few ORFs showed significant differential regulation under DTT or tunicamycin treatment, and scatterplots indicate few transcripts falling outside of the equivalence region, in marked contrast to developmental changes (Table 2, Figure 4E–H). Only one hypothetical chaperone was upregulated after 4 hr in 1 mM DTT, while Rab11, ESAG7, and an ISG-like protein (Tb927.5.410) were downregulated. By contrast, a PDI class ORF, Tb10.6k15.2290, was downregulated after 24 hr of tunicamycin treatment, while ESAG6/7 as well as procyclin were upregulated in cells treated with tunicamycin (Table 2). Further, total cell lysates were examined for increased expression of the major ER chaperone, BiP, by Western blotting and qRT-PCR. No alteration in expression of this marker, which is rapidly and strongly induced by the UPR in higher eukaryotes, was observed following DTT or tunicamycin exposure [see Additional file 4]. These data indicate an essential lack of transcriptome responsiveness to altered conditions by the trypanosome, and specifically the complete absence of a classical UPR. Significantly, this likely contributes to the extreme sensitivity of trypanosomes to DTT.

Knockdown of clathrin and VSG

As an additional strategy for investigating transcriptome responsiveness, the mRNAs of two highly important proteins involved in endocytosis and surface architecture, specifically the clathrin heavy chain (CLH) and VSG, were suppressed with RNAi. Both proteins are essential and CLH RNAi results in rapid cell death. For CLH, knockdown leads to a complete block to endocytosis [114], while VSG RNAi results in arrest of cell cycle progression [115]. We hypothesized that if trypanosomes were able to sense alterations in trafficking and respond to these changes, then depletion of these two ORFs by RNAi would be expected to elicit a response.

We analysed the transcriptome of cells 24 hours after induction of the CLH RNAi, as the CLH protein levels are significantly reduced at that stage [see Additional file 5], growth is arrested, and the big eye morphology is clearly visible (data not shown). Despite a clear phenotype, we observed very limited transcriptome changes (Figure 4K, Table 2). Only TbMSP-C and an ISG-like ORF (Tb927.5.410) were significantly upregulated, while only a putative Rab effector, GMP-PDEδ (Tb927.2.4580) and a putative CBL-interacting protein kinase (Tb927.8.870) were significantly downregulated (Table 2). For the VSG RNAi, we analysed cells 24 and 72 hours after induction to examine both early and late effects of VSG depletion. No major transcriptome changes were observed (Figure 4I–J, Table 2) despite growth arrest [see Additional file 5]. Indeed, only one ORF, a putative RabGEF (Tb927.3.2430) was significantly upregulated (Table 2). The absence of a coordinated response within the transcriptome as a result of these severe perturbations of cellular homeostasis suggests that trypanosomes lack mechanisms for response at the mRNA level to altered circumstances.

The transcriptome is unaltered in cells with augmented endocytic activity

As a final test of transcriptome responsiveness we analysed the effects of over-expression of the early endosome small GTPases Rab5A and Rab5B in PCF, which results in significant augmentation of fluid phase and receptor-mediated endocytic activity [116]. The expression level of CLH was increased, providing a mechanistic explanation for augmented endocytic activity (Figure 5A). These data extend earlier work showing loss of CLH protein following Rab5A or 5B knockdown in BSF [117], and demonstrating coordinated expression of clathrin and Rab5.

Figure 5

Limited differential expression in Rab5-overexpressing PCF cells. Panel A: Western blot analysis for the clathrin heavy chain (CLH) and Rab5A, for wildtype BSFs and PCFs, PCFs overexpressing Rab5A (5AWT), Rab5AQL (GTP-locked mutant) and Rab5BQL (GTP-locked mutant) cell lines. BiP was used as a loading control. Panel B: Scatter plots of raw data for representative microarrays for BSF (Cy5) vs PCF (Cy3) experiment (plot as in Figure 1 for comparison), and Rab5-overexpessing lines (Cy5) vs PCF (Cy3) as indicated. Cy5 fluorescence is plotted on the Y-axis and Cy3 fluorescence on the X-axis. Spots with a Cy5/Cy3 or Cy3/Cy5 ratio above two are highlighted in black, spots with a ratio below two are shown in grey. Panel C: Relative expression levels of Rab5A, Rab5B, CLH, RabX3, COPIε, and adenylate kinase 3 (AK3) in wildtype SMBs and PCFs, and in the Rab5-overexpressor PCF lines, as assessed by qRT-PCR.

Very limited differential expression was seen between Rab5 over-expressing and control cells (Figure 5B, Table 2). Analysis by array, qRT-PCR and Western blot confirmed that the cell lines were overexpressing Rab5A or Rab5B but only detected significant upregulation for RabX3, adenylate kinase 3 (Tb10.70.5150), and COPIε mRNAs in the Rab5AWT, Rab5AQL and Rab5BQL overexpressor lines (Table 2, Figure 5C), none of which are expected as specific changes related to early endosome activity. Further, several mRNAs implicated as altered by array analysis failed to be validated by qRT-PCR, suggesting small or insignificant changes only (Table 2). Importantly, we were unable to detect a significant increase in CLH mRNA in any of the cell lines by qRT-PCR, despite significant protein level changes (Figure 5). These data suggest, firstly, that alteration of CLH levels in these cells proceeds via exclusively post-transcriptional mechanisms and, secondly, that despite clear perturbation to endocytic activity, this is not reflected in an obvious or coherent alteration to the underlying transcriptome.

Significant transcriptome alterations between in vivo and in vitro cultured trypanosomes

The vast majority of molecular studies on T. brucei are carried out using in vitro cultures. Recently, whole-genome microarray analysis indicated significant differences between in vivo and in vitro cultures of the malaria parasite Plasmodium falciparum [118]. Therefore, as an additional investigation into the influence of culturing conditions on intracellular transport we compared mRNA from bloodstream form MITat1.1 trypanosomes extracted from infected rats and from SMB cells cultured under standard conditions. A significant number of differentially abundant mRNAs are evident from the scattergram of these arrays (Figure 4L), and limma analysis indicated that 35 ORFs are significantly upregulated in the in vitro culture, while three ORFs are significantly downregulated (Table 3). These include ORFs for an RNA helicase, procyclin, CAP5.5, several small GTPases and their effectors and several kinases (Table 3). This expression pattern did not reveal an obvious functional pattern, in the sense that ORFs from many biological systems appeared affected. To eliminate the possibility that these alterations arise from cold-shock during isolation of MITat1.1 trypomastigotes, we also extracted RNA from SMB cultures incubated at 4°C for 1 hr. We observed very few changes between normal and cold-shocked SMB cultures (Table 3), with only CAP5.5 and a Sec34-like ORF significantly upregulated in the cold treatment, both of which are procyclic markers. Given the rapidity of isolation of the mRNA from MITat1.1, which is achieved in less than one hour, and their incubation at 4°C during the procedure, which is expected to reduce metabolic activity, we consider these findings to be unlikely ascribable to factors associated with preparation of the samples as the differences in several mRNA levels are too great to have accumulated in such a short period. Therefore, these data highlight significant differences between in vivo and in vitro cultured trypanosomes; this issue would be best pursued by whole genome microarray analysis which is beyond the scope of the present study.

Table 3 Expression changes associated with in vivo versus in vitro culture of T. brucei.


Coordinate control of gene expression is critical for life cycle progression in trypanosomes, and extensive remodeling of the parasite cell is well documented. Here we used the intracellular trafficking system as a model for transcriptome changes. We selected this process on account of well documented alterations to endocytosis accompanying developmental progression, a good level of characterization in trypanosomes, together with, by analogy to higher eukaryotes, the expectation that many gene products will collaborate in regulation. Endomembrane transport involves ~10% of the protein coding potential within the eukaryotic genome [119, 120] and position and time-dependent targeting, assembly of macromolecular complexes and cytoskeletal interactions all participate in controlling a system demanding high fidelity in delivery of molecules to specific subcellular destinations. Further, in higher eukaryotes, transcriptional flexibility and coordination are associated with intracellular transport [17, 121]. Due to the polycistronic mode of transcription in trypanosomes, changes in relative levels of mRNAs result primarily from alterations in the efficiency of nascent RNA processing or half-life [13]. Here we detected clear evidence for developmental alterations to the trypanosome transcriptome, but found very little flexibility within a given life stage.

Firstly, ~6% of transport-associated transcripts are developmentally regulated, a greater fraction than reported for genome-wide transcription [14]. These data indicate that developmental remodeling of endomembrane transport is underpinned by alterations to mRNA abundance, as expected. We also observed prominent differential expression of transport-associated mRNAs between trypanosomes cultured in vitro and those isolated from a mammalian host. A greater proportion of mRNAs were upregulated in the bloodstream stage compared to the insect form, correlating with increased endocytic activity. Prominent upregulated factors include Rab GTPases and SNAREs; significantly, simple upregulation of Rab and/or SNARE proteins is sufficient to augment specific transport pathways [116, 122124]. While the cohort of upregulated mammalian stage mRNAs is consistent with increased endocytic activity, there is evidence that the Rab5/Rab11 VSG endocytosis and recycling system [25, 122, 125] is only part of the developmental change. For PCFs, the upregulated cohort is consistent with lysosomal degradative pathways playing a more prominent role [26]. Direct examination of several of the gene products implicated here is clearly required for a more detailed understanding of trypanosome differentiation.

Secondly, and in clear contrast, extremely limited changes to mRNA levels were encountered in response to a wide variety of altered states, including major changes to culturing conditions, insult with agents disrupting protein folding and RNAi of critical mRNAs encoding factors essential for endocytosis and cell surface maintenance. Despite polycistronic transcription, a priori there is no compelling reason to assume that turnover of mRNA is uncoupled from signaling systems and, moreover, the observation that iron-deficiency leads to increased expression of the ESAG6/7 transferrin receptor mRNAs and protein provides a potential example of signal-mediated changes to mRNA abundance [45, 65, 105]. Analysis of a large gene cohort indicates that the transferrin receptor is a special case, and we were unable to find compelling evidence for alterations to levels of most mRNAs. The transferrin receptor genes are located in an RNA Pol I-transcribed subtelomeric array [126] and thus may be regulated differentially from most genes that are transcribed by RNA Pol II [127]. We also observed mRNA level alterations in a further Pol I-derived transcript, procyclin, in response to several challenges, but we note that procyclin expression is unusually sensitive to many factors (for example [128]). Regardless of these exceptions, our data indicate that the vast majority of trypanosome mRNAs are unmodulated, suggesting that, for most messages, transcription and turnover are unaffected by signaling pathways. This implies that the trypanosome has a very limited capacity to respond to altered circumstances in vivo, except via pre-programmed differentiation pathways.

The absence of transcriptome responsiveness to environmental cues has several consequences. Firstly, lack of a classical UPR indicates ER insensitivity to increased concentrations of non-native polypeptides, and implies that changes in biosynthetic output similarly will not be sensed. This is consistent with extreme sensitivity to concentrations of DTT that are well tolerated in higher eukaryotes. Whether the UPR is absent from additional excavates is unknown, but preliminary evidence suggests that Giardia also lacks the pathway [112]. Secondly, failure to modify the transcriptome following gene knockdown or over-expression, but where a major impact on transport is found, suggests that modulation of trafficking, and potentially other systems, within a given life stage is restricted to proteomic changes [117]. Thirdly, and most significantly, the data indicate an inability of the trypanosome endomembrane system to adapt to altering conditions and, coupled with the possibility that such unresponsiveness extends to additional cellular functions, has profound implications.

Limited transcriptional flexibility is clearly a restriction for the trypanosome, but is in agreement with a modest small GTPase signaling network [30], which may indicate a rather limited signal transduction capacity. As each life stage exists within a highly controlled environment the trypanosome's ability to modulate the trafficking system has likely been lost. The environment within the insect host is poorly understood, but the bloodstream pH range in H. sapiens is 7.35 – 7.45, blood glucose is at 4 – 8 mM and transferrin concentration is maintained at ~300 ± 50 μg/dl [129]; major alterations to these parameters are lethal, and hence there is little need for the parasite to retain mechanisms for responsiveness. By contrast, a new transcriptome would be triggered by appropriate signals when the trypanosome enters a new compartment or host. Hence the trypanosome life cycle may be viewed as a set of inflexible, sequential transcriptional profiles that serve to adapt the parasite to each new environment (Figure 6). The in vitro cultured BSF is a convenient model for the mammalian bloodstream trypomastigote but, during infection, trypanosomes also invade tissue spaces, the lymphatic system and, in later disease stages, the cerebral spinal fluid [1]. Each environment differs greatly from the bloodstream in composition and hence transcriptome inflexibility suggests that a differentiation event may be required for adaptation. Further, the major transcriptome differences between cultured and animal-isolated trypanosomes underscores the potential for currently uncharted, but dramatic, adaptations associated with mammalian infection.

Figure 6

Transcriptional flexibility and inflexibility in differentiation and responsiveness. Upper panel: flexible system. Gene cohorts 1 and 3 are developmentally regulated, and either highly expressed or not expressed; examples of these types of gene products are the trypanosome surface antigens, VSG in the bloodstream form (red) and procyclin in the insect stage (green). The vast majority of genes fall into cohort 2, where, for example, either small or large changes to transcription could result from alterations to the environment (light and dark blue), or a more continually altering transcriptional profile is present that may seek to track changing conditions (purple). This behavior may propagate from one life stage to the next (light and dark blue) or be lost (purple) resulting in altered transcriptional flexibility for genes between life stages. Such a profile is found in higher eukaryotes, including humans and yeast, and probably also many protists, including E. gracilis. Lower panel: inflexible system. In this model gene cohorts 1 and 3 behave as before, but transcription of the genes in cohort 2 remains unchanged. The relative levels of mRNAs from the genes in this cohort may remain constant following differentiation (light blue) or be significantly altered (dark blue and purple). Such a profile is observed here for T. brucei and has been reported previously for P. falciparum, and is potentially a result of a parasitic life style where the host is responsible for provision of a homeostatic environment.

A similar transcriptional rigidity is suggested for the intra-erythrocytic stages of the malaria parasite, Plasmodium falciparum, possibly also related to the fact that Plasmodium enjoys a homeostatic host environment [130132]. Trypanosomatids (Excavata) and Plasmodia (Chromalveolata) are highly divergent and separated by many non-parasitic or free-living species, and thus transcriptome unresponsiveness most likely is the result of independent secondary loss accompanying adoption of the parasitic life style. Consistent with this view, environmental stress treatments do induce transcriptional responses in Euglena gracilis [133], a free-living protist closely related to kinetoplastids. Further, as Plasmodium utilises a more conventional promoter-mediated mechanism for transcription, this also rules out a trivial explanation that transcriptome inflexibility is simply the result of polycistronic transcription. Examination of transcriptome regulation in more free-living organisms related to both Plasmodium and trypanosomes would be highly informative.

Finally, these data indicate that knockdown of a specific mRNA does not alter the transcript levels of additional gene products coordinated functionally with the RNAi target. Therefore, in T. brucei RNAi is unlikely to suffer from off target effects as seen in higher eukaryotes [134136], and hence is, in this regard, unusually clean. Secondary effects are likely mediated via proteome changes, and could explain the similar phenotypes observed for many knockdowns [61]. This is particularly prominent for cell cycle/cytokinesis defects that frequently emerge rather late following knockdown. Such effects may be the result of generic loss of the normal proteome, rather than absence of specific factors.


The life cycle of T. brucei includes a dramatic differentiation event as the parasite progresses from the insect to the mammalian host. Differentiation encompasses considerable remodeling of the membrane trafficking system, which is also a vital component of the host-parasite interface. Using a subgenome microarray targeting the T. brucei trafficking systems together with qRT-PCR analysis, we find clear evidence for developmental alteration to ~6% of the trypanosome transport-associated transcriptome. This indicates that developmental remodeling of endomembrane transport is underpinned by alterations to mRNA abundance, and that such changes may be analysed by microarray. Most protein-coding mRNAs in trypanosomatids are produced by polycistronic transcription, with specific rates of decay as the major mechanism controlling steady state transcript levels. Multiple stimuli, including major changes to culturing conditions, insult with agents disrupting protein folding, and RNAi of critical mRNAs encoding factors essential for endocytosis and cell surface maintenance, elicited extremely limited changes to the transcriptome. Besides alterations in expression of the transferrin receptor upon iron stress, which is transcribed by Pol I, the vast majority of trypanosome chromosomal Pol II-derived mRNAs are unmodulated. These findings imply a very limited capacity of trypanosomes to respond to altered circumstances in vivo, except via pre-determined differentiation pathways (Figure 6). Hence mRNA levels are insensitive to external stimuli and most responses are mediated at the proteome level. These data suggest a new view of the trypanosome life cycle as a progression through sequential, inflexible transcriptional profiles that serve to adapt the parasite to each new host environment. As similar transcriptional rigidity is suggested for the malaria parasite, Plasmodium falciparum, but transcriptional flexibility is found in the closely related and free-living protist, Euglena gracilis, it is likely that loss of transcriptional flexibility arose through the parasitic lifestyle.


Open reading frame selection

We selected ORFs encoding components of the trypanosome transport system using the geneDB interface [137]. The genome was parsed with four criteria; (i) annotation keywords, including description information and wildcards, for protein classes involved in endocytosis and membrane trafficking (e.g. clathrin, adaptin, SNARE, Rab, ADP-ribosylation factor, coatomer, Sec, Vps), (ii) Pfam/Interpro domain annotation (e.g. ARF/SAR superfamily, C2, ENTH, PH, SH3, SNF7, VHS, BAR), (iii) "Biological Process" Gene Ontology (GO) annotation (e.g. endocytosis, vesicle mediated transport, GTPase activity) and (iv) sequence similarity by BLAST. The target gene selection was augmented by parsing the literature for candidate factors, i.e. yeast and mammalian proteins that have been functionally demonstrated to be involved in trafficking. We thus identified 673 ORFs for inclusion.

Array construction

Oligonucleotides were designed with OligoArray 2.0 [138] specifying a 67 mer within one kilobase of the 3'-end of each target protein coding region. 694 oligonucleotides were designed, targeting 796 genes (several ORFs are multicopy and hence multiple transcripts are sampled by a single oligonucleotide in such cases). The array encompassed 56 oligonucleotides targeting 123 reference ORFs, whose expression pattern has been established previously. Oligonucleotides were manufactured by Illumina (Invitrogen) and included a 5' amino modification for attachment to Codelink activated slides (GE Healthcare). Arrays were printed using a BioRobotics Microgrid (Genomic Solutions) equipped with 16 NanoPins (Matrix Technology Corporation) in the Centre for Microarray Resources, Department of Pathology, University of Cambridge. Slides were prepared, washed and dried according to the manufacturer's instructions (GE Healthcare). Each oligonucleotide was spotted in quadruplicate, while α- and β-tubulin were spotted at the corners of the grid and reference oligonucleotides and additional spotting controls (empty, water, spotting buffer) were distributed randomly.

Cells, in vitro and in vivo culture conditions

Cultures of T. brucei brucei Lister 427 single marker bloodstream form (SMB) and procyclic form (PCF) strain 427 were maintained in HMI-9 (Invitrogen) and SDM-79 (Sigma) media, respectively [139141]. Rab5A and Rab5B over-expressing PCF lines have been described previously [116]. A VSG knockdown line was generated by expression of a double strand fragment of the 221 VSG ORF (PCR primers: TTTCTGCAGCGGTCACTATG and TGCTTGTTGTTGGCTTTCAG) cloned into Eam 1105I-digested p2T7TABlue [61]. Tetracycline-responsive SMB cells were transfected with Not I-digested p2T7TABlue-VSG by electroporation, selected and maintained with 2.5 μg/ml G418 and 2.5 μg/ml hygromycin. RNAi was induced by adding 1 μg/ml tetracycline. The clathrin RNAi BSF line was generated as described [114] except that the TbCLH gene fragment was subcloned into the p2T7TA vector [114] and an Amaxa nucleofector was used for transfection [142]. Transformants were maintained under 2.5 μg/ml neomycin and 2.5 μg/ml hygromycin selection, and RNAi was induced by adding 1 μg/ml tetracycline. For treatment with dithiothreitol (DTT) and tunicamycin, SMB cells were grown in HMI-9 to 1 × 106 cells/ml, and the media supplemented with either DTT (Sigma) at 1 mM, or tunicamycin (Sigma) at 5 μg/ml. Samples for RNA extraction were taken after 1 hr and 4 hr for DTT addition and 4 hr and 24 hr for tunicamycin. To assess the influence of serum on gene expression, cells were grown in standard media until a concentration of 1 × 106 was reached and harvested by centrifugation at 800 g for 10 minutes. The cell pellet was re-suspended in media lacking fetal bovine serum (FBS), and supplemented with 5 mg/ml BSA with or without 0.3 mg/ml bovine holo-transferrin (Sigma), to simulate physiological concentrations of albumin and transferrin in media supplemented with 10% FBS [129]. For 30% FBS treatment cells were grown in HMI-9 media supplemented with 30% FBS.

For in vivo grown cells, the blood from rats infected with monomorphic bloodstream forms (MITat1.1 at 1 × 109 cells/ml parasitemia) was taken by aortic puncture into heparinized syringe. The buffy coat fraction was prepared by centrifugation of infected blood (600 g, 10 min at 4°C). All subsequent treatments were performed at 4°C. Trypanosomes were purified from the buffy coat using PSG (NaH2PO4 3 mM, Na2HPO4 50 mM, NaCl 44 mM, sucrose 100 mM, adenosine 0.1 mM, pH 8.0). The cells were washed three times in PSG and then resuspended at 109 cells/ml.

RNA extraction

RNA was extracted using the Qiagen RNeasy kit following the manufacturer's instructions. 1 × 108 BSF (SMB) or 5 × 107 PCF cells were used per extraction. RNA quantity and quality was measured using a NanoDrop ND-1000 (NanoDrop Technologies) and denaturing formaldehyde agarose gels [143] prior to further analysis. For the MITat1.1 cells, RNA was prepared using the Absolutely RNA kit (Stratagene) according to the manufacturer's instructions.

Array hybridization

For each experiment, hybridisations were performed with at least three biological replicates with Cy3/Cy5 dye-swaps. Preparation of amplified and labelled cDNA targets were essentially as described [144] except that AmpliTaq DNA polymerase Buffer II (Applied Biosystems) was used for the cDNA amplification. Following hybridisation, slides were washed for 5 min in 2 × SSC, 5 min in 0.1 × SSC/0.1% SDS, 5 min in 0.1 × SSC and rinsed in dH2O. All washes were performed at room temperature. After washing, slides were dried by centrifugation and scanned on an Axon 4100 scanner using GenePix software (Molecular Devices). Fluorescence intensity was normalised so that the total signal for the Cy3 and the Cy5 channels were equivalent. Raw image data were extracted using Bluefuse (BlueGnome).

Array data analysis

To identify significant differential expression, the microarray data were analysed using limma (version 2.4.7) [50] from the Bioconductor open-source project [49] running under R (version 2.2.1), and a custom program written by K. Kelly (Cambridge). Data pre-processing comprised within-array print-tip Loess normalisation and between-array scale normalisation. To identify differentially expressed mRNAs a linear model was fitted to the data, using as weights the square root of the "confidence" value given by Bluefuse for each spot (which depends on the area, circularity and uniformity of the spot, signal to noise ratio, signal saturation, and the consistency of the Cy3/Cy5 ratio in different pixels across the spot). Adjustment for multiple testing was carried out using the false discovery rate (FDR) method [145] in limma. For clustering analysis, data from the Bluefuse analysis were processed using the Gene Expression Pattern Analysis Suite (GEPAS) [146, 147] using the DNMAD and SOTARRAY programs. Replicates were merged after filtering replicates with values > 0.5 – 1.0 above the median of replicates. Hierarchical clustering was performed using UPGMA and Euclidian distance with settings at normal.

Quantitative real time (qRT)-PCR

For cDNA synthesis, 2 μg RNA was diluted to 10 μl with diethylpyrocarbonate (DEPC)-treated water and denatured at 70°C, 5 min. 15 μl of a reaction mix was added (2.5 μl dNTPs (25 mM stock), 5 μl 5× reverse transcription buffer (Invitrogen), 2 μl 100 mM DTT, 0.5 μl RNAseOUT (recombinant ribonuclease inhibitor, 40 U/μl, Invitrogen), 2 μl oligo dT (T30VN, 10 μM stock), 0.5 μl Superscript II Reverse Transcriptase (200 U/μl Invitrogen), and 2.5 μl DEPC-treated water) and incubated at 37°C for 1 hr, heat-inactivated at 90°C, 5 min and finally diluted to 200 μl with DEPC-treated water.

For qRT-PCR, 5 μl of cDNA was used in a 25 μl reaction including IQ SYBR Green Supermix (BioRad) with 0.4 μM gene-specific forward and reverse primers [see Additional file 2]. qRT-PCR reactions were performed in white thin-wall polypropylene multiplate 96-well unskirted PCR plates (BioRad) sealed with microseal 'B' adhesive (BioRad). Reactions were performed in a BioRad MiniOpticon real time PCR detection system and included an initial denaturation at 95°C for 3 min, 40 cycles of 95°C 30 seconds, 58°C 30 sec, 72°C 6 min (with a signal read at the end of each cycle), and a final melting curve to check fidelity from 60 – 95°C, with a signal read every 1°C. Gene-specific 20 bp primers for each gene were designed using Primer3 [148] specified to amplify a ~120 bp fragment (+/- 10 bp) in the last kilobase of the 3' end of the open reading frame and coinciding with the region recognised by the corresponding array oligonucleotide. Primer pairs were validated in silico using Amplify [149] to minimise the probability of mispriming or formation of primer dimers and secondary structure.

Antibody production and Western blot analysis

TbδAd antisera was generated against an expressed gene fragment (residues 222–508) amplified from T. brucei genomic DNA using Herculase DNA polymerase (Stratagene) with the following primers: 5'-GGAATTCCCAGCTTCTTAGGTCTAGCGGT-3' and 5'-CCGCTCGAGCTCAACAATAGTTCGCATGTGC-3'. The 865 bp PCR product was subcloned into pPCR-Script, excised with Eco RI, and inserted into pGEX-2TK. Polyclonal rabbit antibodies were raised against TbδAd-GST recombinant fusion protein, which was SDS-PAGE-purified due to insolubility when expressed in E. coli, mixed with RIBI adjuvant (Sigma) and used to immunize rabbits on four immunisations spread over a period of five months. 0.5 mg recombinant protein was used per immunisation course.

For Western blotting, 1 × 107 trypanosome cells were harvested by centrifuging at 800 g, 10 min, 4°C, washed in phosphate-buffered saline (PBS) and resuspended in 37.5 μl PBS and 12.5 μl 4× Laemmli SDS loading buffer [143]. Lysates were denatured at 95°C for 10 min and analysed by 12.5% SDS-PAGE [143]. Equivalence of protein loading was verified by Coomassie Blue staining of a duplicate gel. Proteins were electrophoretically transferred to Immobilon-P membrane (Millipore) using a wet transfer tank (Hoefer Instruments). Nonspecific binding was blocked for 1 hr with Tris-buffered saline, pH 7.4, 0.2% Tween (TBST) supplemented with 5% freeze-dried milk. The membrane was then incubated for 1 hr with primary antibody, also diluted in TBST-milk. Rabbit polyclonal antibody against trypanosome BiP (a kind gift of James Bangs, University of Wisconsin) or the trypanosome clathrin heavy chain [16] were used at 1:1,000 dilution, and rabbit anti-TbδAd (see above) at 1:1,000. After washing 3 times for 5 min in TBST, a commercial secondary anti-IgG rabbit horse-radish-peroxidase conjugate (Sigma) was used at 1:10,000 in TBST for 1 hr, and bound-antibodies were detected on Biomax MR-1 films (Kodak) using H2O2-activated Luminol (Sigma) in 100 mM Tris pH 8.5.


  1. 1.

    Barrett MP, Burchmore RJ, Stich A, Lazzari JO, Frasch AC, Cazzulo JJ, Krishna S: The trypanosomiases. Lancet. 2003, 362 (9394): 1469-1480. 10.1016/S0140-6736(03)14694-6.

    Google Scholar 

  2. 2.

    Horn D, Barry JD: The central roles of telomeres and subtelomeres in antigenic variation in African trypanosomes. Chromosome Res. 2005, 13 (5): 525-533. 10.1007/s10577-005-0991-8.

    CAS  Google Scholar 

  3. 3.

    Barry JD: Capping of variable antigen on Trypanosoma brucei, and its immunological and biological significance. J Cell Sci. 1979, 37: 287-302.

    CAS  Google Scholar 

  4. 4.

    O'Beirne C, Lowry CM, Voorheis HP: Both IgM and IgG anti-VSG antibodies initiate a cycle of aggregation-disaggregation of bloodstream forms of Trypanosoma brucei without damage to the parasite. Mol Biochem Parasitol. 1998, 91 (1): 165-193. 10.1016/S0166-6851(97)00191-6.

    Google Scholar 

  5. 5.

    Sternberg JM: Human African trypanosomiasis: clinical presentation and immune response. Parasite Immunol. 2004, 26 (11–12): 469-476. 10.1111/j.0141-9838.2004.00731.x.

    CAS  Google Scholar 

  6. 6.

    Matthews KR: The developmental cell biology of Trypanosoma brucei. J Cell Sci. 2005, 118 (Pt 2): 283-290. 10.1242/jcs.01649.

    CAS  PubMed Central  Google Scholar 

  7. 7.

    Vassella E, Acosta-Serrano A, Studer E, Lee SH, Englund PT, Roditi I: Multiple procyclin isoforms are expressed differentially during the development of insect forms of Trypanosoma brucei. J Mol Biol. 2001, 312 (4): 597-607. 10.1006/jmbi.2001.5004.

    CAS  Google Scholar 

  8. 8.

    Rout MP, Field MC: Isolation and characterization of subnuclear compartments from Trypanosoma brucei. Identification of a major repetitive nuclear lamina component. J Biol Chem. 2001, 276 (41): 38261-38271.

    CAS  Google Scholar 

  9. 9.

    Hammarton TC: Cell cycle regulation in Trypanosoma brucei. Mol Biochem Parasitol. 2007, 153 (1): 1-8. 10.1016/j.molbiopara.2007.01.017.

    CAS  PubMed Central  Google Scholar 

  10. 10.

    Natesan SK, Peacock L, Matthews K, Gibson W, Field MC: Activation of endocytosis as an adaptation to the mammalian host by trypanosomes. Eukaryot Cell. 2007, 6 (11): 2029-2037. 10.1128/EC.00213-07.

    CAS  PubMed Central  Google Scholar 

  11. 11.

    Sharma R, Peacock L, Gluenz E, Gull K, Gibson W, Carrington M: Asymmetric cell division as a route to reduction in cell length and change in cell morphology in trypanosomes. Protist. 2008, 159 (1): 137-151. 10.1016/j.protis.2007.07.004.

    Google Scholar 

  12. 12.

    Urwyler S, Studer E, Renggli CK, Roditi I: A family of stage-specific alanine-rich proteins on the surface of epimastigote forms of Trypanosoma brucei. Mol Microbiol. 2007, 63 (1): 218-228. 10.1111/j.1365-2958.2006.05492.x.

    CAS  Google Scholar 

  13. 13.

    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.

    CAS  PubMed Central  Google Scholar 

  14. 14.

    Brems S, Guilbride DL, Gundlesdodjir-Planck D, Busold C, Luu V-D, Schanne M, Hoheisel J, Clayton C: The transcriptomes of Trypanosoma brucei Lister 427 and TREU927 bloodstream and procyclic trypomastigotes. Mol Biochem Parasitol. 2005, 139 (2): 163-172. 10.1016/j.molbiopara.2004.11.004.

    CAS  Google Scholar 

  15. 15.

    Jeffries TR, Morgan GW, Field MC: A developmentally regulated rab11 homologue in Trypanosoma brucei is involved in recycling processes. J Cell Sci. 2001, 114 (Pt 14): 2617-2626.

    CAS  Google Scholar 

  16. 16.

    Morgan GW, Allen CL, Jeffries TR, Hollinshead M, Field MC: Developmental and morphological regulation of clathrin-mediated endocytosis in Trypanosoma brucei. J Cell Sci. 2001, 114 (Pt 14): 2605-2615.

    CAS  Google Scholar 

  17. 17.

    Gurkan C, Lapp H, Alory C, Su AI, Hogenesch JB, Balch WE: Large-scale profiling of Rab GTPase trafficking networks: the membrome. Mol Biol Cell. 2005, 16 (8): 3847-3864. 10.1091/mbc.E05-01-0062.

    CAS  PubMed Central  Google Scholar 

  18. 18.

    Pelkmans L, Fava E, Grabner H, Hannus M, Habermann B, Krausz E, Zerial M: Genome-wide analysis of human kinases in clathrin- and caveolae/raft-mediated endocytosis. Nature. 2005, 436 (7047): 78-86. 10.1038/nature03571.

    CAS  Google Scholar 

  19. 19.

    Christoforidis S, McBride HM, Burgoyne RD, Zerial M: The Rab5 effector EEA1 is a core component of endosome docking. Nature. 1999, 397 (6720): 621-625. 10.1038/17618.

    CAS  Google Scholar 

  20. 20.

    Zerial M, McBride H: Rab proteins as membrane organizers. Nat Rev Mol Cell Biol. 2001, 2 (2): 107-117. 10.1038/35052055.

    CAS  Google Scholar 

  21. 21.

    McBride HM, Rybin V, Murphy C, Giner A, Teasdale R, Zerial M: Oligomeric complexes link Rab5 effectors with NSF and drive membrane fusion via interactions between EEA1 and syntaxin 13. Cell. 1999, 98 (3): 377-386. 10.1016/S0092-8674(00)81966-2.

    CAS  Google Scholar 

  22. 22.

    Liu W, Duden R, Phair RD, Lippincott-Schwartz J: ArfGAP1 dynamics and its role in COPI coat assembly on Golgi membranes of living cells. J Cell Biol. 2005, 168 (7): 1053-1063. 10.1083/jcb.200410142.

    CAS  PubMed Central  Google Scholar 

  23. 23.

    Stamnes M: Regulating the actin cytoskeleton during vesicular transport. Curr Opin Cell Biol. 2002, 14 (4): 428-433. 10.1016/S0955-0674(02)00349-6.

    CAS  Google Scholar 

  24. 24.

    Randazzo PA: RhoD, Src, and hDia2C in endosome motility. Dev Cell. 2003, 4 (3): 287-288. 10.1016/S1534-5807(03)00067-4.

    CAS  Google Scholar 

  25. 25.

    Field MC, Carrington M: Intracellular membrane transport systems in Trypanosoma brucei. Traffic. 2004, 5 (12): 905-913. 10.1111/j.1600-0854.2004.00234.x.

    CAS  Google Scholar 

  26. 26.

    Field MC, Natesan SK, Gabernet-Castello C, Koumandou VL: Intracellular trafficking in the trypanosomatids. Traffic. 2007, 8 (6): 629-639. 10.1111/j.1600-0854.2007.00558.x.

    CAS  Google Scholar 

  27. 27.

    Engstler M, Pfohl T, Herminghaus S, Boshart M, Wiegertjes G, Heddergott N, Overath P: Hydrodynamic flow-mediated protein sorting on the cell surface of trypanosomes. Cell. 2007, 131 (3): 505-515. 10.1016/j.cell.2007.08.046.

    CAS  Google Scholar 

  28. 28.

    Parsons M, Ruben L: Pathways involved in environmental sensing in trypanosomatids. Parasitol Today. 2000, 16 (2): 56-62. 10.1016/S0169-4758(99)01590-2.

    CAS  Google Scholar 

  29. 29.

    Field MC: Signalling the genome: the Ras-like small GTPase family of trypanosomatids. Trends Parasitol. 2005, 21 (10): 447-450. 10.1016/

    CAS  Google Scholar 

  30. 30.

    Field MC, O'Reilly AJ: How complex is GTPase signaling in trypanosomes?. Trends Parasitol. 2008, 24 (6): 253-7. 10.1016/

    CAS  Google Scholar 

  31. 31.

    Parsons M, Worthey EA, Ward PN, Mottram JC: Comparative analysis of the kinomes of three pathogenic trypanosomatids: Leishmania major, Trypanosoma brucei and Trypanosoma cruzi. BMC Genomics. 2005, 6: 127-10.1186/1471-2164-6-127.

    PubMed Central  Google Scholar 

  32. 32.

    Muller IB, Domenicali-Pfister D, Roditi I, Vassella E: Stage-specific requirement of a mitogen-activated protein kinase by Trypanosoma brucei. Mol Biol Cell. 2002, 13 (11): 3787-3799. 10.1091/mbc.E02-02-0093.

    CAS  PubMed Central  Google Scholar 

  33. 33.

    Domenicali Pfister D, Burkard G, Morand S, Renggli CK, Roditi I, Vassella E: A mitogen-activated protein kinase controls differentiation of bloodstream forms of Trypanosoma brucei. Eukaryot Cell. 2006, 5 (7): 1126-1135. 10.1128/EC.00094-06.

    PubMed Central  Google Scholar 

  34. 34.

    Guttinger A, Schwab C, Morand S, Roditi I, Vassella E: A mitogen-activated protein kinase of Trypanosoma brucei confers resistance to temperature stress. Mol Biochem Parasitol. 2007, 153 (2): 203-206. 10.1016/j.molbiopara.2007.02.001.

    Google Scholar 

  35. 35.

    Hall BS, Gabernet-Castello C, Voak A, Goulding D, Natesan SK, Field MC: TbVps34, the trypanosome orthologue of Vps34, is required for Golgi complex segregation. J Biol Chem. 2006, 281 (37): 27600-27612. 10.1074/jbc.M602183200.

    CAS  Google Scholar 

  36. 36.

    Alexandre S, Paindavoine P, Hanocq-Quertier J, Paturiaux-Hanocq F, Tebabi P, Pays E: Families of adenylate cyclase genes in Trypanosoma brucei. Mol Biochem Parasitol. 1996, 77 (2): 173-182. 10.1016/0166-6851(96)02591-1.

    CAS  Google Scholar 

  37. 37.

    Naula C, Seebeck T: Cyclic AMP signaling in trypanosomatids. Parasitol Today. 2000, 16 (1): 35-38. 10.1016/S0169-4758(99)01582-3.

    CAS  Google Scholar 

  38. 38.

    Dacks JB, Field MC: Eukaryotic cell evolution from a comparative genomic perspective: the endomembrane system. Organelles, genomes and eukaryote phylogeny: an evolutionary synthesis in the age of genomics. Edited by: Hirt R, Horner D. 2004, GRC Press, 309-334.

    Google Scholar 

  39. 39.

    Berriman M, Ghedin E, Hertz-Fowler C, Blandin G, Renauld H, Bartholomeu DC, Lennard NJ, Caler E, Hamlin NE, Haas B, Bohme U, Hannick L, Aslett MA, Shallom J, Marcello L, Hou L, Wickstead B, Alsmark UC, Arrowsmith C, Atkin RJ, Barron AJ, Bringaud F, Brooks K, Carrington M, Cherevach I, Chillingworth TJ, Churcher C, Clark LN, Corton CH, Cronin A, Davies RM, Doggett J, Djikeng A, Feldblyum T, Field MC, Fraser A, Goodhead I, Hance Z, Harper D, Harris BR, Hauser H, Hostetler J, Ivens A, Jagels K, Johnson D, Johnson J, Jones K, Kerhornou AX, Koo H, Larke N, Landfear S, Larkin C, Leech V, Line A, Lord A, Macleod A, Mooney PJ, Moule S, Martin DM, Morgan GW, Mungall K, Norbertczak H, Ormond D, Pai G, Peacock CS, Peterson J, Quail MA, Rabbinowitsch E, Rajandream MA, Reitter C, Salzberg SL, Sanders M, Schobel S, Sharp S, Simmonds M, Simpson AJ, Tallon L, Turner CM, Tait A, Tivey AR, Van Aken S, Walker D, Wanless D, Wang S, White B, White O, Whitehead S, Woodward J, Wortman J, Adams MD, Embley TM, Gull K, Ullu E, Barry JD, Fairlamb AH, Opperdoes F, Barrell BG, Donelson JE, Hall N, Fraser CM, Melville SE, El-Sayed NM: The genome of the African trypanosome Trypanosoma brucei. Science. 2005, 309 (5733): 416-422. 10.1126/science.1112642.

    CAS  Google Scholar 

  40. 40.

    Engstler M, Bangs JD, Field MC: Intracellular transport systems in trypanosomes: function, evolution and virulence. Trypanosomes – after the genome. Edited by: Barry JD, Mottram JC, McCulloch R, Acosta-Serrano A. 2007, Norfolk: Horizon Scientific Press

    Google Scholar 

  41. 41.

    Koumandou VL, Dacks JB, Coulson RM, Field MC: Control systems for membrane fusion in the ancestral eukaryote; evolution of tethering complexes and SM proteins. BMC Evol Biol. 2007, 7: 29-10.1186/1471-2148-7-29.

    PubMed Central  Google Scholar 

  42. 42.

    Lee MG, Bihain BE, Russell DG, Deckelbaum RJ, Ploeg Van der LH: Characterization of a cDNA encoding a cysteine-rich cell surface protein located in the flagellar pocket of the protozoan Trypanosoma brucei. Mol Cell Biol. 1990, 10 (9): 4506-4517.

    CAS  PubMed Central  Google Scholar 

  43. 43.

    Ziegelbauer K, Multhaup G, Overath P: Molecular characterization of two invariant surface glycoproteins specific for the bloodstream stage of Trypanosoma brucei. J Biol Chem. 1992, 267 (15): 10797-10803.

    CAS  Google Scholar 

  44. 44.

    Nolan DP, Jackson DG, Biggs MJ, Brabazon ED, Pays A, Van Laethem F, Paturiaux-Hanocq F, Elliott JF, Voorheis HP, Pays E: Characterization of a novel alanine-rich protein located in surface microdomains in Trypanosoma brucei. J Biol Chem. 2000, 275 (6): 4072-4080. 10.1074/jbc.275.6.4072.

    CAS  Google Scholar 

  45. 45.

    Mussmann R, Engstler M, Gerrits H, Kieft R, Toaldo CB, Onderwater J, Koerten H, van Luenen HGAM, Borst P: Factors affecting the level and localization of the transferrin receptor in Trypanosoma brucei. J Biol Chem. 2004, 279 (39): 40690-40698. 10.1074/jbc.M404697200.

    CAS  Google Scholar 

  46. 46.

    Field MC, Gabernet-Castello C, Dacks JB: Reconstructing the evolution of the endocytic system: insights from genomics and molecular cell biology. Adv Exp Med Biol. 2007, 607: 84-96.

    Google Scholar 

  47. 47.

    Gallop JL, McMahon HT: BAR domains and membrane curvature: bringing your curves to the BAR. Biochem Soc Symp. 2005, 223-231. 72

  48. 48.

    Machida K, Mayer BJ: The SH2 domain: versatile signaling module and pharmaceutical target. Biochim Biophys Acta. 2005, 1747 (1): 1-25.

    CAS  Google Scholar 

  49. 49.

    Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, Ellis B, Gautier L, Ge Y, Gentry J, Hornik K, Hothorn T, Huber W, Iacus S, Irizarry R, Leisch F, Li C, Maechler M, Rossini AJ, Sawitzki G, Smith C, Smyth G, Tierney L, Yang JY, Zhang J: Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004, 5 (10): R80-10.1186/gb-2004-5-10-r80.

    PubMed Central  Google Scholar 

  50. 50.

    Smyth GK: Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004, 3: Article3-

    Google Scholar 

  51. 51.

    Sykacek P, Furlong RA, Micklem G: A friendly statistics package for microarray analysis. Bioinformatics. 2005, 21 (21): 4069-4070. 10.1093/bioinformatics/bti663.

    CAS  Google Scholar 

  52. 52.

    LaCount DJ, Gruszynski AE, Grandgenett PM, Bangs JD, Donelson JE: Expression and function of the Trypanosoma brucei major surface protease (GP63) genes. J Biol Chem. 2003, 278 (27): 24658-24664. 10.1074/jbc.M301451200.

    CAS  Google Scholar 

  53. 53.

    Price HP, Stark M, Smith DF: Trypanosoma brucei ARF1 plays a central role in endocytosis and golgi-lysosome trafficking. Mol Biol Cell. 2007, 18 (3): 864-873. 10.1091/mbc.E06-08-0736.

    CAS  PubMed Central  Google Scholar 

  54. 54.

    Dhir V, Goulding D, Field MC: TbRAB1 and TbRAB2 mediate trafficking through the early secretory pathway of Trypanosoma brucei. Mol Biochem Parasitol. 2004, 137 (2): 253-265. 10.1016/j.molbiopara.2004.06.007.

    CAS  Google Scholar 

  55. 55.

    Sahin A, Lemercier G, Tetaud E, Espiau B, Myler P, Stuart K, Bakalara N, Merlin G: Trypanosomatid flagellum biogenesis: ARL-3A is involved in several species. Exp Parasitol. 2004, 108 (3–4): 126-133. 10.1016/j.exppara.2004.07.010.

    CAS  Google Scholar 

  56. 56.

    Ackers JP, Dhir V, Field MC: A bioinformatic analysis of the RAB genes of Trypanosoma brucei. Mol Biochem Parasitol. 2005, 141 (1): 89-97. 10.1016/j.molbiopara.2005.01.017.

    CAS  Google Scholar 

  57. 57.

    Abeliovich H, Grote E, Novick P, Ferro-Novick S: Tlg2p, a yeast syntaxin homolog that resides on the Golgi and endocytic structures. J Biol Chem. 1998, 273 (19): 11719-11727. 10.1074/jbc.273.19.11719.

    CAS  Google Scholar 

  58. 58.

    Uemura T, Ueda T, Ohniwa RL, Nakano A, Takeyasu K, Sato MH: Systematic analysis of SNARE molecules in Arabidopsis : dissection of the post-Golgi network in plant cells. Cell Struct Funct. 2004, 29 (2): 49-65. 10.1247/csf.29.49.

    CAS  Google Scholar 

  59. 59.

    Besteiro S, Coombs GH, Mottram JC: The SNARE protein family of Leishmania major. BMC Genomics. 2006, 7: 250-10.1186/1471-2164-7-250.

    PubMed Central  Google Scholar 

  60. 60.

    Jahn R, Scheller RH: SNAREs – engines for membrane fusion. Nat Rev Mol Cell Biol. 2006, 7 (9): 631-643. 10.1038/nrm2002.

    CAS  Google Scholar 

  61. 61.

    Subramaniam C, Veazey P, Redmond S, Hayes-Sinclair J, Chambers E, Carrington M, Gull K, Matthews K, Horn D, Field MC: Chromosome-wide analysis of gene function by RNA interference in the african trypanosome. Eukaryot Cell. 2006, 5 (9): 1539-1549. 10.1128/EC.00141-06.

    CAS  PubMed Central  Google Scholar 

  62. 62.

    Belden WJ, Barlowe C: Distinct roles for the cytoplasmic tail sequences of Emp24p and Erv25p in transport between the endoplasmic reticulum and Golgi complex. J Biol Chem. 2001, 276 (46): 43040-43048. 10.1074/jbc.M108113200.

    CAS  Google Scholar 

  63. 63.

    Saint-Pol A, Yelamos B, Amessou M, Mills IG, Dugast M, Tenza D, Schu P, Antony C, McMahon HT, Lamaze C, Johannes L: Clathrin adaptor epsinR is required for retrograde sorting on early endosomal membranes. Dev Cell. 2004, 6 (4): 525-538. 10.1016/S1534-5807(04)00100-5.

    CAS  Google Scholar 

  64. 64.

    Allen CL, Liao D, Chung WL, Field MC: Dileucine signal-dependent and AP-1-independent targeting of a lysosomal glycoprotein in Trypanosoma brucei. Mol Biochem Parasitol. 2007, 156 (2): 175-190. 10.1016/j.molbiopara.2007.07.020.

    CAS  Google Scholar 

  65. 65.

    Hall BS, Smith E, Langer W, Jacobs LA, Goulding D, Field MC: Developmental variation in Rab11-dependent trafficking in Trypanosoma brucei. Eukaryot Cell. 2005, 4 (5): 971-980. 10.1128/EC.4.5.971-980.2005.

    CAS  PubMed Central  Google Scholar 

  66. 66.

    Helms MJ, Ambit A, Appleton P, Tetley L, Coombs GH, Mottram JC: Bloodstream form Trypanosoma brucei depend upon multiple metacaspases associated with RAB11-positive endosomes. J Cell Sci. 2006, 119 (Pt 6): 1105-1117. 10.1242/jcs.02809.

    CAS  Google Scholar 

  67. 67.

    Morty RE, Vadasz I, Bulau P, Dive V, Oliveira V, Seeger W, Juliano L: Tropolysin, a new oligopeptidase from African trypanosomes. Biochemistry. 2005, 44 (44): 14658-14669. 10.1021/bi051035k.

    CAS  Google Scholar 

  68. 68.

    Caffrey CR, Hansell E, Lucas KD, Brinen LS, Alvarez Hernandez A, Cheng J, Gwaltney SL, Roush WR, Stierhof YD, Bogyo M, Steverding D, McKerrow JH: Active site mapping, biochemical properties and subcellular localization of rhodesain, the major cysteine protease of Trypanosoma brucei rhodesiense. Mol Biochem Parasitol. 2001, 118 (1): 61-73. 10.1016/S0166-6851(01)00368-1.

    CAS  Google Scholar 

  69. 69.

    Mottram JC, Helms MJ, Coombs GH, Sajid M: Clan CD cysteine peptidases of parasitic protozoa. Trends Parasitol. 2003, 19 (4): 182-187. 10.1016/S1471-4922(03)00038-2.

    CAS  Google Scholar 

  70. 70.

    Anelli T, Alessio M, Mezghrani A, Simmen T, Talamo F, Bachi A, Sitia R: ERp44, a novel endoplasmic reticulum folding assistant of the thioredoxin family. Embo J. 2002, 21 (4): 835-844. 10.1093/emboj/21.4.835.

    CAS  PubMed Central  Google Scholar 

  71. 71.

    Hsu MP, Muhich ML, Boothroyd JC: A developmentally regulated gene of trypanosomes encodes a homologue of rat protein-disulfide isomerase and phosphoinositol-phospholipase C. Biochemistry. 1989, 28 (15): 6440-6446. 10.1021/bi00441a042.

    CAS  Google Scholar 

  72. 72.

    Rubotham J, Woods K, Garcia-Salcedo JA, Pays E, Nolan DP: Characterization of two protein disulfide isomerases from the endocytic pathway of bloodstream forms of Trypanosoma brucei. J Biol Chem. 2005, 280 (11): 10410-10418. 10.1074/jbc.M409375200.

    CAS  Google Scholar 

  73. 73.

    Bringaud F, Peyruchaud S, Baltz D, Giroud C, Simpson L, Baltz T: Molecular characterization of the mitochondrial heat shock protein 60 gene from Trypanosoma brucei. Mol Biochem Parasitol. 1995, 74 (1): 119-123. 10.1016/0166-6851(95)02486-7.

    CAS  Google Scholar 

  74. 74.

    Tyler KM, Matthews KR, Gull K: The bloodstream differentiation-division of Trypanosoma brucei studied using mitochondrial markers. Proc Biol Sci. 1997, 264 (1387): 1481-1490. 10.1098/rspb.1997.0205.

    CAS  PubMed Central  Google Scholar 

  75. 75.

    Chiang HL, Terlecky SR, Plant CP, Dice JF: A role for a 70-kilodalton heat shock protein in lysosomal degradation of intracellular proteins. Science. 1989, 246 (4928): 382-385. 10.1126/science.2799391.

    CAS  Google Scholar 

  76. 76.

    Gething MJ, Sambrook J: Protein folding in the cell. Nature. 1992, 355 (6355): 33-45. 10.1038/355033a0.

    CAS  Google Scholar 

  77. 77.

    Eisenberg E, Greene LE: Multiple roles of auxilin and hsc70 in clathrin-mediated endocytosis. Traffic. 2007, 8 (6): 640-646. 10.1111/j.1600-0854.2007.00568.x.

    CAS  Google Scholar 

  78. 78.

    Lee MG, Ploeg Van der LH: Transcription of the heat shock 70 locus in Trypanosoma brucei. Mol Biochem Parasitol. 1990, 41 (2): 221-231. 10.1016/0166-6851(90)90185-O.

    CAS  Google Scholar 

  79. 79.

    Rodgers MJ, Albanesi JP, Phillips MA: Phosphatidylinositol 4-kinase III-beta is required for golgi maintenance and cytokinesis in Trypanosoma brucei. Eukaryot Cell. 2007, 6 (7): 1108-1118. 10.1128/EC.00107-07.

    CAS  PubMed Central  Google Scholar 

  80. 80.

    Hart DT, Misset O, Edwards SW, Opperdoes FR: A comparison of the glycosomes (microbodies) isolated from Trypanosoma brucei bloodstream form and cultured procyclic trypomastigotes. Mol Biochem Parasitol. 1984, 12 (1): 25-35. 10.1016/0166-6851(84)90041-0.

    CAS  Google Scholar 

  81. 81.

    Barnard JP, Pedersen PL: Alteration of pyruvate metabolism in African trypanosomes during differentiation from bloodstream into insect forms. Arch Biochem Biophys. 1994, 313 (1): 77-82. 10.1006/abbi.1994.1361.

    CAS  Google Scholar 

  82. 82.

    Soderling TR: The Ca-calmodulin-dependent protein kinase cascade. Trends Biochem Sci. 1999, 24 (6): 232-236. 10.1016/S0968-0004(99)01383-3.

    CAS  Google Scholar 

  83. 83.

    Postel EH, Weiss VH, Beneken J, Kirtane A: Mutational analysis of NM23-H2/NDP kinase identifies the structural domains critical to recognition of a c-myc regulatory element. Proc Natl Acad Sci USA. 1996, 93 (14): 6892-6897. 10.1073/pnas.93.14.6892.

    CAS  PubMed Central  Google Scholar 

  84. 84.

    Nosaka K, Kawahara M, Masuda M, Satomi Y, Nishino H: Association of nucleoside diphosphate kinase nm23-H2 with human telomeres. Biochem Biophys Res Commun. 1998, 243 (2): 342-348. 10.1006/bbrc.1997.8097.

    CAS  Google Scholar 

  85. 85.

    Chopra P, Koduri H, Singh R, Koul A, Ghildiyal M, Sharma K, Tyagi AK, Singh Y: Nucleoside diphosphate kinase of Mycobacterium tuberculosis acts as GTPase-activating protein for Rho-GTPases. FEBS Lett. 2004, 571 (1–3): 212-216. 10.1016/j.febslet.2004.06.073.

    CAS  Google Scholar 

  86. 86.

    Hunger-Glaser I, Hemphill A, Shalaby T, Hanni M, Seebeck T: Nucleoside diphosphate kinase of Trypanosoma brucei. Gene. 2000, 257 (2): 251-257. 10.1016/S0378-1119(00)00401-7.

    CAS  Google Scholar 

  87. 87.

    Steinborn K, Szallies A, Mecke D, Duszenko M: Cloning, heterologous expression and kinetic analysis of glycerol kinase (TbGLK1) from Trypanosoma brucei. Biol Chem. 2000, 381 (11): 1071-1077. 10.1515/BC.2000.132.

    CAS  Google Scholar 

  88. 88.

    Pereira CA, Alonso GD, Torres HN, Flawia MM: Arginine kinase: a common feature for management of energy reserves in African and American flagellated trypanosomatids. J Eukaryot Microbiol. 2002, 49 (1): 82-85. 10.1111/j.1550-7408.2002.tb00346.x.

    CAS  Google Scholar 

  89. 89.

    Haanstra JR, Stewart M, Luu VD, van Tuijl A, Westerhoff HV, Clayton C, Bakker BM: Control and Regulation of Gene Expression: quantitative analysis of the expression of phosphoglycerate kinase in bloodstream form Trypanosoma brucei. J Biol Chem. 2008, 283 (5): 2495-2507. 10.1074/jbc.M705782200.

    CAS  Google Scholar 

  90. 90.

    Katzmann DJ, Odorizzi G, Emr SD: Receptor downregulation and multivesicular-body sorting. Nat Rev Mol Cell Biol. 2002, 3 (12): 893-905. 10.1038/nrm973.

    CAS  Google Scholar 

  91. 91.

    Williams RL, Urbe S: The emerging shape of the ESCRT machinery. Nat Rev Mol Cell Biol. 2007, 8 (5): 355-368. 10.1038/nrm2162.

    CAS  Google Scholar 

  92. 92.

    Chung WL, Carrington M, Field MC: Cytoplasmic targeting signals in transmembrane invariant surface glycoproteins of trypanosomes. J Biol Chem. 2004, 279 (52): 54887-54895. 10.1074/jbc.M409311200.

    CAS  Google Scholar 

  93. 93.

    Chung W-L, Leung KF, Natesan SKA, Carrington M, Field MC: Ubiquitylation is required for degradation of trans-membrane surface proteins in trypanosomes. Traffic. 2008,

    Google Scholar 

  94. 94.

    Leung KF, Dacks JB, Field MC: Evolution of the multivesicular body ESCRT machinery; retention across the eukaryotic lineage. Traffic. 2008,

    Google Scholar 

  95. 95.

    Davletov BA, Sudhof TC: A single C2 domain from synaptotagmin I is sufficient for high affinity Ca2+/phospholipid binding. J Biol Chem. 1993, 268 (35): 26386-26390.

    CAS  Google Scholar 

  96. 96.

    Gawler DJ: Points of convergence between Ca2+ and Ras signalling pathways. Biochim Biophys Acta. 1998, 1448 (2): 171-182. 10.1016/S0167-4889(98)00141-4.

    CAS  Google Scholar 

  97. 97.

    Katan M, Allen VL: Modular PH and C2 domains in membrane attachment and other functions. FEBS Lett. 1999, 452 (1–2): 36-40. 10.1016/S0014-5793(99)00531-1.

    CAS  Google Scholar 

  98. 98.

    Baines A, Gull K: WCB is a C2 domain protein defining the plasma membrane – sub-pellicular microtubule corset of kinetoplastid parasites. Protist. 2008, 159 (1): 115-125. 10.1016/j.protis.2007.09.001.

    CAS  Google Scholar 

  99. 99.

    Nozaki T, Haynes PA, Cross GA: Characterization of the Trypanosoma brucei homologue of a Trypanosoma cruzi flagellum-adhesion glycoprotein. Mol Biochem Parasitol. 1996, 82 (2): 245-255. 10.1016/0166-6851(96)02741-7.

    CAS  Google Scholar 

  100. 100.

    LaCount DJ, Barrett B, Donelson JE: Trypanosoma brucei FLA1 is required for flagellum attachment and cytokinesis. J Biol Chem. 2002, 277 (20): 17580-17588. 10.1074/jbc.M200873200.

    CAS  Google Scholar 

  101. 101.

    DeRisi JL, Iyer VR, Brown PO: Exploring the metabolic and genetic control of gene expression on a genomic scale. Science. 1997, 278 (5338): 680-686. 10.1126/science.278.5338.680.

    CAS  Google Scholar 

  102. 102.

    Eisen MB, Spellman PT, Brown PO, Botstein D: Cluster analysis and display of genome-wide expression patterns. Proc Natl Acad Sci USA. 1998, 95 (25): 14863-14868. 10.1073/pnas.95.25.14863.

    CAS  PubMed Central  Google Scholar 

  103. 103.

    Hughes TR, Marton MJ, Jones AR, Roberts CJ, Stoughton R, Armour CD, Bennett HA, Coffey E, Dai H, He YD, Kidd MJ, King AM, Meyer MR, Slade D, Lum PY, Stepaniants SB, Shoemaker DD, Gachotte D, Chakraburtty K, Simon J, Bard M, Friend SH: Functional discovery via a compendium of expression profiles. Cell. 2000, 102 (1): 109-126. 10.1016/S0092-8674(00)00015-5.

    CAS  Google Scholar 

  104. 104.

    Travers KJ, Patil CK, Wodicka L, Lockhart DJ, Weissman JS, Walter P: Functional and genomic analyses reveal an essential coordination between the unfolded protein response and ER-associated degradation. Cell. 2000, 101 (3): 249-258. 10.1016/S0092-8674(00)80835-1.

    CAS  Google Scholar 

  105. 105.

    Fast B, Kremp K, Boshart M, Steverding D: Iron-dependent regulation of transferrin receptor expression in Trypanosoma brucei. Biochem J. 1999, 342 (Pt 3): 691-696. 10.1042/0264-6021:3420691.

    CAS  PubMed Central  Google Scholar 

  106. 106.

    Patil C, Walter P: Intracellular signaling from the endoplasmic reticulum to the nucleus: the unfolded protein response in yeast and mammals. Curr Opin Cell Biol. 2001, 13 (3): 349-355. 10.1016/S0955-0674(00)00219-2.

    CAS  Google Scholar 

  107. 107.

    Schroder M, Kaufman RJ: The mammalian unfolded protein response. Annu Rev Biochem. 2005, 74: 739-789. 10.1146/annurev.biochem.73.011303.074134.

    Google Scholar 

  108. 108.

    Arvas M, Pakula T, Lanthaler K, Saloheimo M, Valkonen M, Suortti T, Robson G, Penttila M: Common features and interesting differences in transcriptional responses to secretion stress in the fungi Trichoderma reesei and Saccharomyces cerevisiae. BMC Genomics. 2006, 7: 32-10.1186/1471-2164-7-32.

    PubMed Central  Google Scholar 

  109. 109.

    Benedetti C, Fabbri M, Sitia R, Cabibbo A: Aspects of gene regulation during the UPR in human cells. Biochem Biophys Res Commun. 2000, 278 (3): 530-536. 10.1006/bbrc.2000.3838.

    CAS  Google Scholar 

  110. 110.

    Liang SH, Zhang W, McGrath BC, Zhang P, Cavener DR: PERK (eIF2alpha kinase) is required to activate the stress-activated MAPKs and induce the expression of immediate-early genes upon disruption of ER calcium homoeostasis. Biochem J. 2006, 393 (Pt 1): 201-209.

    CAS  PubMed Central  Google Scholar 

  111. 111.

    Martinez IM, Chrispeels MJ: Genomic analysis of the unfolded protein response in Arabidopsis shows its connection to important cellular processes. Plant Cell. 2003, 15 (2): 561-576. 10.1105/tpc.007609.

    CAS  PubMed Central  Google Scholar 

  112. 112.

    Reiner DS, McCaffery JM, Gillin FD: Reversible interruption of Giardia lamblia cyst wall protein transport in a novel regulated secretory pathway. Cell Microbiol. 2001, 3 (7): 459-472. 10.1046/j.1462-5822.2001.00129.x.

    CAS  Google Scholar 

  113. 113.

    Ferguson MA, Duszenko M, Lamont GS, Overath P, Cross GA: Biosynthesis of Trypanosoma brucei variant surface glycoproteins. N-glycosylation and addition of a phosphatidylinositol membrane anchor. J Biol Chem. 1986, 261 (1): 356-362.

    CAS  Google Scholar 

  114. 114.

    Allen CL, Goulding D, Field MC: Clathrin-mediated endocytosis is essential in Trypanosoma brucei. Embo J. 2003, 22 (19): 4991-5002. 10.1093/emboj/cdg481.

    CAS  PubMed Central  Google Scholar 

  115. 115.

    Sheader K, Vaughan S, Minchin J, Hughes K, Gull K, Rudenko G: Variant surface glycoprotein RNA interference triggers a precytokinesis cell cycle arrest in African trypanosomes. Proc Natl Acad Sci USA. 2005, 102 (24): 8716-8721. 10.1073/pnas.0501886102.

    CAS  PubMed Central  Google Scholar 

  116. 116.

    Pal A, Hall BS, Nesbeth DN, Field HI, Field MC: Differential endocytic functions of Trypanosoma brucei Rab5 isoforms reveal a glycosylphosphatidylinositol-specific endosomal pathway. J Biol Chem. 2002, 277 (11): 9529-9539. 10.1074/jbc.M110055200.

    CAS  Google Scholar 

  117. 117.

    Hall B, Allen CL, Goulding D, Field MC: Both of the Rab5 subfamily small GTPases of Trypanosoma brucei are essential and required for endocytosis. Mol Biochem Parasitol. 2004, 138 (1): 67-77. 10.1016/j.molbiopara.2004.07.007.

    CAS  Google Scholar 

  118. 118.

    Daily JP, Le Roch KG, Sarr O, Ndiaye D, Lukens A, Zhou Y, Ndir O, Mboup S, Sultan A, Winzeler EA, Wirth DF: In vivo transcriptome of Plasmodium falciparum reveals overexpression of transcripts that encode surface proteins. J Infect Dis. 2005, 191 (7): 1196-1203. 10.1086/428289.

    CAS  PubMed Central  Google Scholar 

  119. 119.

    Bonangelino CJ, Chavez EM, Bonifacino JS: Genomic screen for vacuolar protein sorting genes in Saccharomyces cerevisiae. Mol Biol Cell. 2002, 13 (7): 2486-2501. 10.1091/mbc.02-01-0005.

    CAS  PubMed Central  Google Scholar 

  120. 120.

    Quenneville NR, Conibear E: Toward the systems biology of vesicle transport. Traffic. 2006, 7 (7): 761-768. 10.1111/j.1600-0854.2006.00436.x.

    CAS  Google Scholar 

  121. 121.

    Schuldiner M, Collins SR, Thompson NJ, Denic V, Bhamidipati A, Punna T, Ihmels J, Andrews B, Boone C, Greenblatt JF, Weissman JS, Krogan NJ: Exploration of the function and organization of the yeast early secretory pathway through an epistatic miniarray profile. Cell. 2005, 123 (3): 507-519. 10.1016/j.cell.2005.08.031.

    CAS  Google Scholar 

  122. 122.

    Pal A, Hall BS, Jeffries TR, Field MC: Rab5 and Rab11 mediate transferrin and anti-variant surface glycoprotein antibody recycling in Trypanosoma brucei. Biochem J. 2003, 374 (Pt 2): 443-451. 10.1042/BJ20030469.

    CAS  PubMed Central  Google Scholar 

  123. 123.

    Wick P, Gansel X, Oulevey C, Page V, Studer I, Durst M, Sticher L: The expression of the t-SNARE AtSNAP33 is induced by pathogens and mechanical stimulation. Plant Physiol. 2003, 132 (1): 343-351. 10.1104/pp.102.012633.

    CAS  PubMed Central  Google Scholar 

  124. 124.

    Spurlin BA, Park SY, Nevins AK, Kim JK, Thurmond DC: Syntaxin 4 transgenic mice exhibit enhanced insulin-mediated glucose uptake in skeletal muscle. Diabetes. 2004, 53 (9): 2223-2231. 10.2337/diabetes.53.9.2223.

    CAS  Google Scholar 

  125. 125.

    Grunfelder CG, Engstler M, Weise F, Schwarz H, Stierhof YD, Morgan GW, Field MC, Overath P: Endocytosis of a glycosylphosphatidylinositol-anchored protein via clathrin-coated vesicles, sorting by default in endosomes, and exocytosis via RAB11-positive carriers. Mol Biol Cell. 2003, 14 (5): 2029-2040. 10.1091/mbc.E02-10-0640.

    CAS  PubMed Central  Google Scholar 

  126. 126.

    Gunzl A, Bruderer T, Laufer G, Schimanski B, Tu LC, Chung HM, Lee PT, Lee MG: RNA polymerase I transcribes procyclin genes and variant surface glycoprotein gene expression sites in Trypanosoma brucei. Eukaryot Cell. 2003, 2 (3): 542-551. 10.1128/EC.2.3.542-551.2003.

    PubMed Central  Google Scholar 

  127. 127.

    Palenchar JB, Bellofatto V: Gene transcription in trypanosomes. Mol Biochem Parasitol. 2006, 146 (2): 135-141. 10.1016/j.molbiopara.2005.12.008.

    CAS  Google Scholar 

  128. 128.

    Engstler M, Boshart M: Cold shock and regulation of surface protein trafficking convey sensitization to inducers of stage differentiation in Trypanosoma brucei. Genes Dev. 2004, 18: 2798-2811. 10.1101/gad.323404.

    CAS  PubMed Central  Google Scholar 

  129. 129.

    Putnam FW: The plasma proteins. 1984, London: Academic Press, IV: 2

    Google Scholar 

  130. 130.

    Wirth DF: Biological revelations. Nature. 2002, 419 (6906): 495-496. 10.1038/419495a.

    CAS  Google Scholar 

  131. 131.

    Coulson RM, Hall N, Ouzounis CA: Comparative genomics of transcriptional control in the human malaria parasite Plasmodium falciparum. Genome Res. 2004, 14 (8): 1548-1554. 10.1101/gr.2218604.

    CAS  PubMed Central  Google Scholar 

  132. 132.

    Llinas M, Bozdech Z, Wong ED, Adai AT, DeRisi JL: Comparative whole genome transcriptome analysis of three Plasmodium falciparum strains. Nucleic Acids Res. 2006, 34 (4): 1166-1173. 10.1093/nar/gkj517.

    CAS  PubMed Central  Google Scholar 

  133. 133.

    Dos Santos Ferreira V, Rocchetta I, Conforti V, Bench S, Feldman R, Levin MJ: Gene expression patterns in Euglena gracilis : insights into the cellular response to environmental stress. Gene. 2007, 389 (2): 136-145. 10.1016/j.gene.2006.10.023.

    CAS  Google Scholar 

  134. 134.

    Jackson AL, Bartz SR, Schelter J, Kobayashi SV, Burchard J, Mao M, Li B, Cavet G, Linsley PS: Expression profiling reveals off-target gene regulation by RNAi. Nat Biotechnol. 2003, 21 (6): 635-637. 10.1038/nbt831.

    CAS  Google Scholar 

  135. 135.

    Aagaard L, Rossi JJ: RNAi therapeutics: principles, prospects and challenges. Adv Drug Deliv Rev. 2007, 59 (2–3): 75-86. 10.1016/j.addr.2007.03.005.

    CAS  PubMed Central  Google Scholar 

  136. 136.

    Sarov M, Stewart AF: The best control for the specificity of RNAi. Trends Biotechnol. 2005, 23 (9): 446-448.

    CAS  Google Scholar 

  137. 137.

    GeneDB database. []

  138. 138.

    Rouillard J-M, Zuker M, Gulari E: OligoArray 2.0: design of oligonucleotide probes for DNA microarrays using a thermodynamic approach. Nucleic Acids Res. 2003, 31 (12): 3057-3062. 10.1093/nar/gkg426.

    CAS  PubMed Central  Google Scholar 

  139. 139.

    Brun R, Schonenberger : Cultivation and in vitro cloning of procyclic culture forms of Trypanosoma brucei in a semi-defined medium. Short communication. Acta Trop. 1979, 36 (3): 289-292.

    CAS  Google Scholar 

  140. 140.

    Hirumi H, Hirumi K: Axenic culture of African trypanosome bloodstream forms. Parasitol Today. 1994, 10 (2): 80-84. 10.1016/0169-4758(94)90402-2.

    CAS  Google Scholar 

  141. 141.

    Wirtz E, Leal S, Ochatt C, Cross GA: A tightly regulated inducible expression system for conditional gene knock-outs and dominant-negative genetics in Trypanosoma brucei. Mol Biochem Parasitol. 1999, 99 (1): 89-101. 10.1016/S0166-6851(99)00002-X.

    CAS  Google Scholar 

  142. 142.

    Burkard G, Fragoso CM, Roditi I: Highly efficient stable transformation of bloodstream forms of Trypanosoma brucei. Mol Biochem Parasitol. 2007, 153 (2): 220-223. 10.1016/j.molbiopara.2007.02.008.

    CAS  Google Scholar 

  143. 143.

    Sambrook J, Fritsch EF, Maniatis T: Molecular cloning: A laboratory manual. 1989, Cold Spring Harbor Press, 2

    Google Scholar 

  144. 144.

    Petalidis L, Bhattacharyya S, Morris GA, Collins VP, Freeman TC, Lyons PA: Global amplification of mRNA by template-switching PCR: linearity and application to microarray analysis. Nucleic Acids Res. 2003, 31 (22): e142-10.1093/nar/gng142.

    CAS  PubMed Central  Google Scholar 

  145. 145.

    Reiner A, Yekutieli D, Benjamini Y: Identifying differentially expressed genes using false discovery rate controlling procedures. Bioinformatics. 2003, 19 (3): 368-375. 10.1093/bioinformatics/btf877.

    CAS  Google Scholar 

  146. 146.

    Gene Expression Pattern Analysis Suite. []

  147. 147.

    Herrero J, Al-Shahrour F, Diaz-Uriarte R, Mateos A, Vaquerizas JM, Santoyo J, Dopazo J: GEPAS: A web-based resource for microarray gene expression data analysis. Nucleic Acids Res. 2003, 31 (13): 3461-3467. 10.1093/nar/gkg591.

    CAS  PubMed Central  Google Scholar 

  148. 148.

    Primer3 primer design program. []

  149. 149.

    Amplify. []

Download references


We are grateful to the following for assistance with ORF selection: Amanda O'Reilly, Carme Gabernet-Castello, Sara Melville, Mark Carrington and Jeremy Mottram. We also thank Sergio Callejas and Anthony Brown for assistance with array design, Krystyna Kelly and Rob Furlong for advice with data analysis and Karl Hoffmann and Jennifer Fitzpatrick for help with qRT-PCR. Special thanks to Clare Allen for production of the antibody to TbδAd and to Derek Nolan for providing mRNA from MITat1.1 trypanosomes. We also thank Joel Dacks and Mark Carrington for comments on the manuscript. This work was supported by program and project grants from the Wellcome Trust (to MCF). The microarray designs used in this study as well as the relevant experimental data have been deposited in ArrayExpress under accession numbers: A-MEXP-963, A-MEXP-1068, E-MEXP-1528, E-MEXP-1529, E-MEXP-1530, E-MEXP-1533, E-MEXP-1539, E-MEXP-1540.

Author information



Corresponding author

Correspondence to Mark C Field.

Additional information

Authors' contributions

VLK and MCF designed the microarray and the experimental approach, VLK performed the microarray analysis and data interpretation, TS generated and validated the VSG RNAi line and SKAN performed immunochemical analysis of the Rab5 overexpressor lines. VLK and MCF drafted the manuscript and figures; all authors read and approved the final manuscript.

Electronic supplementary material


Additional file 1: Absence of clustering of membrane-trafficking genes in the T. brucei genome. Panel A: The number of genes identified on each chromosome predicted to have a role in membrane traffic based on sequence similarity and/or domain architecture are plotted against chromosome size (in Mb). Multicopy genes are represented here as a single gene. Panel B: Raw data on which plot A is based. Approximate chromosome sizes are from GeneDB. (PDF 22 KB)

Details of open reading frame selection, microarray oligonucleotides, qRT-PCR primers, database annotation and prior work

Additional file 2: . Functional groups assigned based on likely orthologs and/or domain annotation from multiple sources. Accession numbers are for T. brucei 927 at GeneDB, except ESAG6 and ESAG7 which are for T. brucei 427 at GeneDB, and VSG oligonucleotides which are based on sequences deposited at NCBI. All oligonucleotide sequences are given in 5'-3' orientation. Array version 1.0 was used for the bloodstream v procyclic comparison and analysis of the Rab5 overexpressor lines; array version 1.1 used in all other experiments. The "Cross-hybridising ORFs" column gives the GeneDB accession numbers of ORFS likely to cross-hybridise to each oligonucleotide, when applicable, as predicted by the OligoArray 2.0 program. qRT-PCR oligonucleotide sequences are given for all ORFs analysed by qRT-PCR. (XLS 176 KB)

Developmental regulation of components of the adaptin complexes AP1, AP3, and AP4 in

Additional file 3: T. brucei. Panel A: Heatmap of bloodstream versus procyclic form signal ratios for the ORFs corresponding to the adaptin complexes, based on eight microarray experiments. For each target gene, the four replicate spots on the array were averaged, following removal of inconsistent replicates. The scale shows the colour scheme for the z-score of the data, indicating how far and in what direction, the ratio for each spot deviates from the mean for each array, expressed in units of standard deviation; bright red indicates significant upregulation in BSF, bright green indicates significant upregulation in PCF, dark colours or black indicate no differential expression between the two developmental stages. AP1σ is highly upregulated in BSF, and AP3μ in PCF (also see Table 1) with a general trend of BSF upregulation for AP1 components and PCF upregulation (or BSF downregulation) for AP3 and AP4 components. Panel B: Western blot analysis of whole cell lysates from BSF (lane B) and PCF (lane P) (1 × 107 cell equivalents) probed with rabbit antisera raised against recombinant TbδAd adaptin subunit (AP3). The scale to the left represents relative molecular mass in kDa. The predicted molecular weight for the AP3δ protein is 125 kDa. Loading equivalence was monitored by Ponceau Red staining of the membrane after transfer (not shown). (PDF 124 KB)

Response of

Additional file 4: T. brucei cells to dithiothreitol and tunicamycin treatment. Panel A: Growth curves for BSF cultures after the addition of dithiothreitol (DTT, 1–10 mM final concentrations) or tunicamycin (5–10 μg/ml final concentrations). Cell numbers diminish rapidly after the addition of DTT (within 4 hours), whereas cell growth is arrested after addition of tunicamycin and cell numbers remain stable for up to 24 hours. Panel B: Western blot analysis of whole cell lysates (1 × 107 cell equivalents) for BiP in a control culture, as well as cultures supplemented with 1 mM DTT and 5 μg/ml tunicamycin. A culture of BSF cells was divided; one flask was retained as control and DTT or tunicamycin was added to the other subcultures. Samples were taken prior to DTT or tunicamycin addition (C0), 1 hr and 4 hr after DTT treatment (D1 and D4, respectively), 1 hr, 4 hr and 24 hr after tunicamycin treatment (T1, T4 and T24), and at 1 hr, 4 hr and 24 hr from the control (C1, C4, C24). Panel C: BiP protein levels normalised to total protein, based on the Western blot in panel B and two additional replicates. The levels of BiP antigen are normalized to 100% for the control culture at time zero (C0). Panel D: Relative levels of BiP mRNA as measured by qRT-PCR after DTT (black bars) and tunicamycin (grey bars) treatment. The levels of BiP mRNA are normalized to 100% for the control culture at time zero (C0). (PDF 223 KB)


Additional file 5: Verification of phenotypes for VSG and CLH RNAi cell lines. Panel A: Growth curves of induced versus uninduced VSG RNAi cell lines, indicating growth arrest in the induced cell line over four days post induction. Panel B: CLH protein levels in induced (+Tet) versus uninduced (-Tet) CLH RNAi cell lines were examined by Western blotting. BiP was used as a loading control. Data from two biological replicates are shown. (PDF 75 KB)

Significantly developmentally expressed trypanosome genes grouped by functional class

Additional file 6: . Numbers of transcripts on which the graph in Figure 2 is based. The number of genes for each functional class correspond to the ORFs featured in Table 1 and discussed in the text. Note that several of the oligonucleotides on the array target multicopy genes (e.g, tubulin, histones) and thus the total number of ORFs subject to differential regulation is higher than the number of ORFs/oligonucleotides given in Table 1. (PDF 22 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and Permissions

About this article

Cite this article

Koumandou, V.L., Natesan, S.K.A., Sergeenko, T. et al. The trypanosome transcriptome is remodelled during differentiation but displays limited responsiveness within life stages. BMC Genomics 9, 298 (2008).

Download citation


  • Unfolded Protein Response
  • Tunicamycin
  • Bloodstream Form
  • Variant Surface Glycoprotein
  • Clathrin Heavy Chain