BACKGROUND
all organisms are capable of adapting to their environment, which is usually achieved by adjusting gene expression levels, often at transcription initiation.
it has been the goal of many studies to determine the changes in transcription in response to varying environmental conditions, such as when the cells are under stress, drug pressure, in different environments or subject to immune responses.
traditional genome-wide analysis of gene expression of cells under different conditions or, in the case of parasites, at different life cycle stages, has mainly been carried out by hybridization-based methods such as microarrays  <cit> .
such hybridization-based approaches are subject to non-specific hybridization, cross hybridization and nonlinear and saturable hybridization kinetics, providing relative rather than direct quantitative expression levels that are not likely to be comparable between experiments <cit> .
two alternative approaches to gene expression analysis are sequence based and have become increasingly popular due to recent developments in sequencing technologies <cit> .
the first is the direct sequencing of cdna, termed rnaseq  <cit> .
the second approach is based on sequencing serial analysis of gene expression  libraries ), a method that generates a digital output proportional to the number of transcripts per mrna <cit> .
these methods are limited only by the depth of coverage of the sequencing undertaken and both approaches have the benefit of not requiring presynthesised oligonucleotide probes , allowing the direct enumeration of transcript molecules, i.e.
digital quantification, which is directly comparable across different experiments.
the sage/dge approach  is based on the addition of specific adapters to poly cdna, which has been digested with a restriction enzyme, typically nlaiii.
further digestion of the cdna with an enzyme that recognises within the adapter but cleaves  <dig> bp downstream generates a  <dig> base tag from that transcript and unlike traditional sage, where the tags are ligated together and then sequenced , the tags are directly sequenced using next generation sequencing technology, allowing direct quantification of the number of transcripts in a sample.
this process will generate a sequence tag from transcripts that contain nlaiii sites, preferentially from the 3' end of the transcript.
this is typically either in the 3' untranslated regions  or the 3' end of the open reading frame .
transcription profiling using the sage and more recently the dge method has been used extensively to compare expression profiles in a range of tissues  <cit> .
the advantage of dge is the larger dynamic range obtained per experiment.
however, there are a number of disadvantages with the dge approach.
it does not generate the added value of determining 5' and 3' utrs or alternative splice variants, which is an important consideration for the analysis of most organisms, .
dge also requires a reference transcriptome for comparison and is confined to the analysis of genes containing the four-base restriction site  used in the library construction.
if the 5' and 3' utrs of genes are not included in the reference transcriptome, transcripts with nlaiii sites in the 3'utrs or within  <dig> bp of the stop codon will not be represented.
despite only obtaining a partial transcriptome of the organism under investigation, dge is a useful tool in the analysis of gene expression at an unprecedented level of sensitivity, particularly when comparing two very similar samples.
this sensitive and cost-effective approach can be applied to the study of gene expression in cells under a spectrum of experimental or natural conditions.
one group of organisms that is subjected to dramatic environmental changes throughout their life cycle, including larges changes in temperature, nutrients and host immune defences is the parasitic protozoa.
the evolutionarily ancient single-celled parasite trypanosoma brucei gambiense has a particularly complicated life cycle, passing between the tsetse vector and different mammalian hosts species including humans, where it causes a fatal disease in humans, african sleeping sickness <cit> .
throughout its life cycle, the parasite has to adapt to extreme changes in its environment, in response to which it differentiates into several morphologically distinct forms, involving organelle repositioning and changes to its metabolism <cit> .
to what extent are these gross morphological and metabolic changes controlled by alterations in the transcriptome of the parasites between different life cycle stages?
studies to date investigating the transcriptome of african trypanosomes have focused mainly on the non-human infective sub-species, t. b. brucei, revealing that transcription in trypanosomes is unusual.
the majority of trypanosome genes are transcribed by rna polymerase ii in long polycistronic units, which are then processed by trans-splicing that adds a 39-nucleotide spliced leader to every mrna to produce monocistronic mature mrnas  <cit> .
however evidence to date suggests that, unlike bacterial operons, most genes in polycistronic units are not functionally related and are not co-regulated.
the steady state levels of mrna appear to be determined predominantly post-transcriptionally by mrna degradation controlled by the 3'utr of transcripts <cit> .
there are notable exceptions to this pattern of gene regulation in that some gene arrays that encode the parasite's variant surface glycoproteins  and procyclin expressions sites in bloodstream forms and procyclic forms, respectively, are co-regulated and co-transcribed by rna polymerase i <cit> .
this raises the question, are there other co-regulated polycistronic units within the genome?
if so, when are they expressed and what do they encode?
are co-regulated polycistronic units transcribed by rna polymerase i alone?
while many genes in the genome appear to be arranged in the same orientation and are presumably transcribed as polycistronic units, several small open reading frames , typically between 150- <dig> bp in size, with the opposite orientation have been annotated in the genome as being 'unlikely' <cit> .
do these small orfs represent a novel class of trypanosome genes that are transcribed in the opposite orientation from other genes in the polycistron or are they an artefact of gene prediction software?
in order to address these questions an analysis of the transcript abundance of these orfs is required.
microarray analysis has been carried out to examine changes in transcript abundance between the mammalian bloodstream forms and the insect procyclic forms of the related parasite t. b. brucei and the results suggest that relatively few parasite genes are differentially regulated between the two life cycle stages <cit> .
in this study we add to the transcriptome analysis of trypanosomes, demonstrating the use of dge to examine the transcriptome of two life cycle stages of the clinically relevant sub-species, t. b. gambiense.
RESULTS
library production and dge tag annotation
in this study we determined and compared the partial transcriptome of two different life cycle stages of the pleomorphic t. b. gambiense strain, stib 386; the procyclic  stage, which is not human-infective, and the potentially human-infective bloodstream stage.
rna from parasites of both life cycle stages was prepared from exponentially growing cells.
three dge libraries were made from replicate cultures for each life cycle stage and sequenced using solexa  technology.
the three procyclic form libraries  and three bloodstream form libraries  generated  <dig>  ×  <dig> and  <dig>  ×  <dig> tags, respectively, for the two classes of samples, after filtering for poor quality sequencing scores .
the normalized tag counts are presented as supplementary data  and all raw data is available from gene expression omnibus  accession number gse <dig> and will be made available on tritrypdb <cit> .
the number of tags generated using this methodology is at least 300-fold greater than the number of tags obtained for a t. b. brucei strain, using sage with traditional sequencing methods <cit> .
the tags were aligned to the fully annotated reference transcriptome of t. b. brucei, treu  <dig> <cit>  allowing for a two base pair mismatch for any polymorphisms between the reference genome and the genome of stib  <dig>  as the 5' and 3' utr of most t. b. brucei genes have not be defined, the annotated transcriptome contained only orfs, limiting this analysis to tags that align uniquely to the orf of genes.
the reference transcriptome represented the chromosome internal genes and did not extend past the sub-telomeric regions of the chromosomes or include the minichromosomes.
the number of tags that were of good quality and aligned uniquely to the reference genome for each library is presented in table  <dig>  providing expression data for  <dig>  genes.
however there were a large number of tags that either aligned to multiple sites or no sites within the annotated reference transcriptome, in part reflecting the incomplete nature of the in silico transcriptome and the distribution of nlaiii sites in the genome.
it is likely that many of the tags that did not align to the transcriptome actually aligned to the 3'utr of transcripts as genes in the t. brucei genome have relatively long 3'utrs with a median size of  <dig> nt <cit> .
the dynamic range of dge spanned four orders of magnitude.
the most abundant transcript in all procyclic form libraries was that for the heat shock gene, hsp <dig> , with a maximum tag count in procyclic form library a of  <dig>  tags accounting for  <dig> % of uniquely aligned tags for that sample.
however, the tag counts for the majority of genes were low in both the bloodstream and procyclic forms.
reproducibility
digital gene expression has been shown to have very low technical variability <cit> .
in order to examine the ability of dge to detect subtle differences between biological replicates, dge libraries generated from rna extracted from three identical cultures grown under the same conditions for each life cycle stage were compared.
the number of tags generated for each library was similar .
figure  <dig> shows a series of scatter plots comparing the normalised tag counts  for the replicate libraries in each pairwise combination.
a high correlation was observed between replicates, indicating a high degree of reproducibility of technical and biological replicates, with an average pearson's correlation coefficient of r =  <dig>  for the procyclic form replicates and r =  <dig>  for the bloodstream form replicates following log transformation .
bloodstream form parasites when grown at high density are capable of differentiating into the next life cycle stage, the short stumpy form, which is cell cycle arrested and pre-adapted to life in the tsetse fly <cit> .
the inability to synchronise pleomorphic cells could result in the presence in some cultures of short stumpy differentiated forms, resulting in reduced reproducibility in bloodstream form cultures compared to procyclic forms.
although no short stumpy forms were observed in any of the bloodstream form cultures, and the overall correlation coefficient for bloodstream form replicates was reasonably high, similar to that for procyclic forms, a small proportion of trypanosomes could have begun to differentiate into stumpy forms.
indeed, pad <dig>  a gene previously found to be enriched in stumpy forms over slender bloodstream forms and procyclics, was up-regulated 5-fold in the bloodstream forms in this analysis, although not statistically significant at the  <dig>  false discovery rate  .
this indicated that some of the bloodstream form cells used in this analysis may have begun the differentiation process to stumpy forms.
the variation between replicate cultures, although small, clearly demonstrated the value of performing biological replicates.
differentially expressed genes
the expressed tags that uniquely aligned to the reference transcriptome generated expression data for  <dig>  genes, approximately 81% of the total number of genes  in the annotated genome <cit>  .
using a threshold of  <dig>  average fold change in expression for comparison with previous studies using microarrays,  <dig>  genes were up-regulated in bloodstream forms and only  <dig> genes up-regulated in procyclic forms, a total of 28% differentially expressed genes , compared to 2% using microarray analysis  <cit> .
however, this difference was less an indication of biological significance than a product of the arbitrary nature of using a fold-change cut-off, which was insensitive to differences in factors affecting statistical power such as technical sensitivity, sample size and variability among biological replicates.
in this analysis we have employed a statistically valid method of identifying differentially expressed genes taking into account biological replicates, termed rank products <cit> , which ranked genes by degree of differential expression and estimated both gene-specific p-values and a fdr that allowed many hypotheses to be tested simultaneously.
intuitively, setting a low, stringent fdr would define a smaller but better supported set of genes.
using a highly stringent fdr of  <dig> , where 10% of identified genes would be expected to be false positives, a total of  <dig> genes were up-regulated in bloodstream forms compared to procyclics and  <dig> genes were up-regulated in procyclic forms compared to bloodstream forms.
at a slightly less stringent fdr of  <dig> ,  <dig> genes were up-regulated in bloodstream forms and  <dig> genes were up-regulated in procyclic forms .
the overall difference in regulation is illustrated by a quantile-quantile  plot  where the observed -log <dig> for each gene was plotted against the expected value under the null hypothesis of no difference in gene expression between life cycle stages.
the observed values for both bloodstream and procyclic stages rose above the y = x line at low p-values  > 1) indicating that there were many more low p-values than expected by chance, i.e.
many genes were strongly differentially expressed, with more genes differentially expressed in bloodstream stages than in procyclic stages.
throughout this analysis all genes with differential gene expression are referred to as being up-regulated either in bloodstream or procyclics.
however this does not imply a regulatory mechanism merely a comparative description that could have been generated by a down-regulation in one of the life-cycle stages.
remarkably the fold change of many of these differentially expressed genes was extremely high , ranging from  <dig> to 634-fold up-regulated in bloodstream forms and  <dig>  to  <dig>  up-regulated in procyclic forms.
the distribution of the detected fold change is illustrated in the volcano plot , where the statistical significance of each gene was plotted against fold change.
the plot shows that the tags with the highest average differences between life cycle stages  also showed the greatest degree of significance.
another feature of the plot was the tendency of procyclics to show more significant differences than bloodstream forms at fold changes close to zero, reflecting the greater reproducibility within procyclics leading to higher statistical power relative to the bloodstream group.
differentially expressed genes, up-regulated in bloodstream forms
at the  <dig>  fdr,  <dig> genes have been identified as being up-regulated in bloodstream forms.
thirty two of these were hypothetical genes of unknown function , four of which  have been previously identified in t. b. brucei as being up-regulated in bloodstream forms in a proteome study <cit>  and one  by microarray analysis <cit> .
the identification of these possible coding sequences as being up-regulated in bloodstream forms indicates that they have been correctly classed as genes in the t. b. brucei genome confirming the gene model used for their identification.
as expected, a large proportion of genes  that were differentially regulated at the  <dig>  fdr level were associated with antigenic variation; the vsg genes, and related expression site associated genes , both of which were expressed at very high levels in bloodstream forms  <cit> .
the vsgs are the major surface coat proteins unique to bloodstream forms, which are encoded by members of a large gene family, although only one vsg gene is expressed at any one time.
the vsg gene repertoire is extremely diverse with very few vsg gene sequences being shared between strains <cit> .
in this analysis, there were three predominant vsg genes that together represent 66% of the vsg tags that uniquely aligned to the reference transcriptome.
they aligned to tb <dig> .
<dig>  tb <dig> .
<dig> and tb09v <dig>  pseudogenes in the treu  <dig> reference transcriptome.
one hundred and two other vsg genes that uniquely aligned the reference transcriptome were also expressed within this cell population at low levels.
it is likely that there were many other vsg gene tags present in the library that did not align to the reference transcriptome, due to strain specific nature of repertoire of vsg genes.
according to current dogma, only one expression site, and therefore only one vsg gene, is active in any one trypanosome cell <cit> , this suggests that there were multiple populations of cells expressing different vsg genes within the bloodstream form cultures analysed.
an alternative explanation would be that there were low levels of multiple vsg mrnas being expressed at the one time in one cell.
the difference in metabolism between the two different life-cycle stages due to the change in the main energy source, from glucose in the bloodstream to proline in procyclic forms <cit> , was reflected in the up-regulation in bloodstream forms of the tht1- gene encoding a glucose transporter, , a gene encoding a putative glycerol uptake protein  and the ald gene encoding fructose-bisphosphate aldolase , consistent with previous microarray or proteome studies <cit> .
other genes that were up-regulated in bloodstream forms include those involved in intracellular protein transport or modification, i.e.
the p <dig> gene that encodes a lysosomal/endosomal membrane protein , nsf gene that encodes the vesicular-fusion protein  and the gpi-plc gene that encodes glycosylphosphatidylinositol-specific phospholipase c .
the gene tb <dig> .
<dig> encoding a putative iron/ascorbate oxidoreductase family protein has not been previously identified as being up-regulated in bloodstream forms using microarray analysis <cit> , although it did appear to be up-regulated in bloodstream forms at the protein level <cit> .
although this gene was detected only at a low level in bloodstream forms  and not detected at all in any of the procyclic form libraries, it was still possible to identify a highly significant difference in transcript abundance, illustrating the sensitivity of digital tag technology to detect differences even for low-abundance transcripts.
similarly, three genes encoding leucine-rich repeat proteins that were up-regulated in bloodstream forms were either not detected at all in the procyclic form libraries  or at extremely low levels .
other genes that were up-regulated in bloodstream forms at the  <dig>  fdr include genes that encode dynein heavy chain proteins putatively involved in intraflagellar transport , the chaperone protein, dnaj , sumo, involved in post-translational modification , glutathionylspermidine synthetase  and the cysteine peptidase, calpain , none of which had previously been shown to be up-regulated in bloodstream forms in t. brucei by microarray analysis <cit> .
by setting a stringent fdr, a small number of well-supported genes were identified in this study, however this would undoubtedly exclude a number of genuinely differentially expressed genes, some of which were previously identified in other studies as being differentially expressed.
for comparison, we have taken these previously identified differentially expressed genes described in the microarray studies of koumandou et al.
<cit>  and brems et al.
<cit>  and compared them, where possible, with the dge data in this analysis .
this shows that the vast majority of genes that were previously identified as being up-regulated in bloodstream forms in t. b. brucei, were also found to be up-regulated in the related parasite, t. b. gambiense, albeit with approximately half of these genes having a p-value greater than  <dig> .
a differentially expressed cluster of genes
of the  <dig> genes that were up-regulated in bloodstream forms above the  <dig>  fdr threshold, five contiguous genes  were up-regulated between  <dig> and 90-fold .
these genes, located at a convergent strand switch region at an internal site on chromosome  <dig>  were three esags  and two procyclin associated genes .
although most esags are expressed from telomeric expression sites  <cit> , there is evidence that some esags can also be expressed from chromosome internal regions of the genome <cit> .
the identification of two pag genes that were highly expressed in bloodstream forms compared to procyclic forms was somewhat surprising, as related genes in this multi-gene family were more highly expressed in the procyclic forms in t. brucei  <cit> .
indeed, copies of pag <dig> and pag <dig> found on chromosome  <dig> have been shown to be expressed from a polycictronic unit containing the procyclic surface coat genes, ep <dig> and ep2procyclin <cit> , which were highly expressed in procyclic forms .
the functions of the pag proteins are unknown but there is strong amino acid similarity to the transferrin binding proteins, esag <dig> and esag <dig> <cit> .
in order to verify that these genes were up-regulated in bloodstream forms, rt-pcr was performed for each gene in turn, using primers designed to unique regions of each gene.
the rt-pcr clearly confirmed the differential expression of all five genes .
rnaseq data, released prior to publication, for the related sub-species t. b. brucei also clearly showed the steady-state rna levels of these genes was far greater in bloodstream forms than procyclics <cit> .
differentially expressed gene clusters were defined as two or more contiguous genes, which were significantly differentially expressed at the  <dig>  fdr.
the gene accession identification number, protein description and average fold change is given for each gene.
neighbouring genes within a polycistronic transcription unit are usually regulated independently due to variable mrna stability, largely regulated by sequences in the 3'utr of the mrna.
such post-transcriptional regulation of gene expression has become a distinguishing feature of t. b. brucei and the basis of much research <cit> .
here five consecutive genes were co-regulated and highly expressed in bloodstream form trypanosomes.
these genes are situated next to a region where neighbouring polycistronic units meet, termed a strand switch region.
interestingly, divergent and convergent strand switch regions being associated with transcription start and stop sites, respectively <cit> .
for this gene cluster, the strand switch region is convergent and likely to contain transcription stop sites.
there are remarkably short intergenic regions between each of the five orf in this transcription unit, ranging in size from  <dig> bp to  <dig> bp, which could indicate extremely short 3'utrs of these genes and perhaps the lack of any signals for mrna degradation.
such highly co-regulated genes are reminiscent of the bloodstream form specific telomeric expression sites, which are unusual in that they are transcribed by rna polymerase i <cit> .
other rna polymerase i transcription units include the chromosome internal rrna loci and the ep procyclin and gpeet loci, which are up-regulated in procyclic forms in t. b. brucei and which also contain members of the pag multi-gene family <cit> .
however, there was no obvious sequences homology between regions upstream of this polycistronic unit and the rna polymerase i promoters of the procyclin and vsg loci <cit> .
other differentially expressed gene clusters up-regulated in bloodstream forms
systematic analysis of three or more contiguous differentially expressed genes  revealed four additional differentially expressed gene clusters  that were up-regulated in bloodstream forms and one cluster up-regulated in procyclic forms .
gene cluster  <dig> contains four genes that were up-regulated  <dig> to 122- fold on chromosome  <dig>  this gene cluster is present in a polycistronic unit that is positioned at a convergent strand switch region and encodes a member of the pag gene family , in a similar manner to the differentially expressed gene cluster on chromosome  <dig>  two other differentially expressed gene clusters also contain esag <dig> and esag <dig> and, in all cases, the differentially expressed genes appear not to be syntenic with other kinetoplastids and most contain hypothetical genes unique to t. brucei.
all gene clusters described also show differential mrna levels in the rnaseq data of t. b. brucei <cit> , except for cluster  <dig>  which was possibly a sub-species differentially expressed gene cluster.
what role these genes play in the parasite's life cycle clearly warrants further investigation.
differentially expressed genes, up-regulated in procyclics
at the  <dig>  fdr level,  <dig> genes were identified as being up-regulated in procyclic forms  fewer than were up-regulated in bloodstream forms.
the most up-regulated gene identified in procyclic form trypanosomes encoded the procyclic stage specific surface antigen pssa- <dig>   <cit> .
the genes encoding the major procyclic surface coat protein procyclin, ep <dig> and ep <dig> <cit> , do not contain nlaiii sites within the orf, and so these genes did not appear in the in silico derived transcriptome to which the tags were aligned.
however a tag matching the 3'utr of ep <dig> was identified in the dataset of unaligned tags and was expressed at extremely high level in procyclic forms, up-regulated 250-fold.
a large proportion of the genes up-regulated in procyclic forms were annotated as hypothetical  providing some insight into their possible role in trypanosome development.
interestingly, several small open reading frames annotated as being unlikely to encode proteins <cit>  were identified as being differentially expressed and particularly up-regulated in procyclic forms , suggesting the gene annotation is correct.
known developmentally regulated mrnas involved in mitochondrial biogenesis and metabolism in procyclic forms that were found to be differentially expressed at the protein level were identified in the present analysis, gmdh  and ppdk  <cit> .
while many of the cytochrome oxidase subunit genes  were identified as being up-regulated in procyclic forms <cit>  at the  <dig>  fdr, they fall just below the stringent  <dig>  fdr threshold .
other genes up-regulated in procyclic forms include tbts , an amino acid transporter gene , a gene encoding a calpain-like protein  and the surface protein gene, cram , which have all been previously identified as being up-regulated in procyclic form trypanosomes  <cit> .
several other genes identified in the present study as being up-regulated in procyclic forms were not identified in previous studies at the mrna level, including a putative udp-gal or udp-glcnac-dependent glycosytransferase gene , although this was found to be up-regulated in procyclics in a comparative proteome study <cit> .
interestingly there are  <dig> tandemly repeated copies of this gene, but it is only the last gene in the array that was up-regulated in procyclic forms, the rest being up-regulated approximately three-fold in bloodstream forms .
intriguingly, while the other genes in the tandem array share near-identical 3'utr sequences, tb <dig> .
<dig> has a very different 3'utr, consistent with the hypothesis that the 3'utr contains regulatory elements that determine differential gene expression <cit> .
previous analysis of the tandemly repeated gene cluster of phosphoglycerol kinase gene copies , have indicated that each gene is regulated differently in the polycistronic unit, with pgka being constitutively expressed, pgkb being procyclic form specific and pgkc being bloodstream form specific  <cit> .
analysis of the expression patterns of these genes using dge indicated that very few unique tags aligned to the orf of these genes, partly due to the similarity in gene sequence between all three pgk gene copies reducing the number of informative tags.
in addition, the majority of tags for pgkc aligned to the 3'utr, which was not represented in the reference transcriptome.
the few unique tags that were generated for these genes did not show any significant difference in tag numbers between life cycle stages.
assuming the expression pattern for these genes is the same in t.b.
gambiense and t.b.
brucei, the analysis of this gene family illustrates that this technique only provided relative expression differences for a portion of the differentially expressed transcripts and did not detect all differentially expressed genes.
the most highly expressed gene in procyclic forms was hsp <dig>  although this gene was also among the top four genes expressed in all bloodstream form libraries, it appeared to be up-regulated in procyclic forms  <dig> -fold .
interestingly, hsp <dig> has been previously identified as being up-regulated in procyclic forms of the pleomorphic strain, t. b. brucei  <dig>  but not the monomorphic strain t. b. brucei lister  <dig> <cit> , indicating its differential expression was strain specific.
whole genome sequencing of the stib  <dig> strain revealed an estimated  <dig> copies of this gene in the genome , which could explain the high expression levels.
comparing the previously identified genes that have been described as being up-regulated in procyclic forms in the microarray studies of koumandou et al.
<cit>  and brems et al <cit>  for t. b. brucei with the dge data in this analysis , many genes appear to be up-regulated in both types of analysis, confirming previous analysis and suggesting similarities in differentiation pathways between the different sub-species.
however there were several genes, which gave contradictory results, i.e.
they appeared to be down-regulated  in procyclic forms, albeit not statistically significantly .
CONCLUSIONS
in this study, the partial transcriptomes of bloodstream and procyclic cells of t. b. gambiense were compared using dge, providing expression data for  <dig> genes, constituting 81% of the genome.
however this still leaves 19% the genome with no expression information.
some of those genes could be specific to the other life cycle stages not analysed here, such as epimastigotes, while others could be unrepresented due to the limitations of the dge technique.
these limitations include only sampling genes with nlaiii sites within their open reading frames, an underrepresentation of genes with long 3'utrs and the inability to align many tags uniquely to the transcriptome.
while dge allows transcription profiling of a large proportion of the genes in the genome, a more comprehensive analysis of the transcriptome could be have been obtained by rnaseq analysis.
in addition, allowing for a two base pair mismatch between the tag and the reference transcriptome to take account of sequence differences between the t. b. brucei and t. b. gambiense genomes could result in misalignment of the tag to the reference transcriptome and could also introduce a bias in some of the analysis.
however, this approach is particularly useful when comparing paired samples, in this case different life cycle stages of the same strain and the concordance between the data generated by dge and that of microarrays, albeit for a different sub-species, indicates that this approach is valid.
dge has revealed that a larger proportion of genes in t. b. gambiense are stage-regulated in contrast to two previous microarray studies on t. b. brucei.
a recent third microarray study <cit>  examining gene expression in a time course of parasite differentiation revealed many of the differentially expressed genes identified here were common to all studies.
where this study differs from microarray analysis is in the far greater sensitivity of dge compared to hybridization-based techniques, which coupled with its high reproducibility means that a larger number of genes that had not previously been identified as being stage-regulated could be identified as being differentially expressed with a high degree of statistical power.
the second reason for the difference in gene expression between the dge analysis and the previous microarray studies is that different strains/sub-species were used.
strains of t. brucei can vary naturally in a range of different phenotypes, for example, in terms of virulence, growth, drug resistance and resistance to human serum.
undoubtedly such phenotypic differences will be as a result of different gene expression patterns, although strain-specific differences have been relatively under investigated in this species with the majority of research being focused on a single laboratory-adapted line.
the increased power to detect strain-specific differences means that dge could be a useful tool to examine naturally occurring variable phenotypes, or strains that have been selected for particular phenotypes, and to understand the pathways involved in these processes.
one particular phenomenon revealed by dge, which had not previously been observed from microarray studies, but which was evident in the pre-publication of rnaseq data of t. brucei <cit>  is the presence of several differentially expressed clusters of genes present on what appears to be the same polycistronic units.
these clusters are reminiscent of the rna polymerase i transcribed polycistronic units involved in the expression of vsg and procyclin genes, which encode the parasite's surface coat.
many of these differentially expressed gene clusters are located in regions of the genome that are unique to t. brucei, with no synteny with other kinetoplastids.
these regions are also associated with strand switch sites.
the strand switch regions of differentially expressed gene clusters identified here are convergent and do not obviously correspond to the possible rna polymerase ii transcription start sites previously identified <cit> .
it is tempting to speculate that these islands of differentially expressed genes are transcribed by the rna polymerase that controls transcription of other differentially expressed genes, i.e.
rna polymerase i. indeed four of the gene clusters contain either pags or esags, families of genes that are often transcribed by rna polymerase i in either procyclin or bloodstream expression sites, respectively.
however at present it is unclear which rna polymerase transcribes these life cycle specific gene clusters, although understanding the basis of the control of the differentially expressed gene clusters clearly warrants further investigation to determine if this phenomenon is due to transcript stability of mrna or transcription rates.
in conclusion we have shown that dge analysis generates a wealth of data, revealing a large number of hitherto unknown differentially regulated genes and providing new insights into this unusual pathogen.

 all organisms are capable of adapting to their environment, which is usually achieved by adjusting gene expression levels, often at transcription initiation. it has been the goal of many studies to determine the changes in transcription in response to varying environmental conditions, such as when the cells are under stress, drug pressure, in different environments or subject to immune responses. traditional genome-wide analysis of gene expression of cells under different conditions or, in the case of parasites, at different life cycle stages, has mainly been carried out by hybridization-based methods such as microarrays. such hybridization-based approaches are subject to non-specific hybridization, cross hybridization and nonlinear and saturable hybridization kinetics, providing relative rather than direct quantitative expression levels that are not likely to be comparable between experiments. two alternative approaches to gene expression analysis are sequence based and have become increasingly popular due to recent developments in sequencing technologies. the first is the direct sequencing of cdna, termed rnaseq. the second approach is based on sequencing serial analysis of gene expression libraries ), a method that generates a digital output proportional to the number of transcripts per mrna. these methods are limited only by the depth of coverage of the sequencing undertaken and both approaches have the benefit of not requiring presynthesised oligonucleotide probes, allowing the direct enumeration of transcript molecules, i.e. digital quantification, which is directly comparable across different experiments. the sage/dge approach is based on the addition of specific adapters to poly cdna, which has been digested with a restriction enzyme, typically nlaiii. further digestion of the cdna with an enzyme that recognises within the adapter but cleaves bp downstream generates a base tag from that transcript and unlike traditional sage, where the tags are ligated together and then sequenced, the tags are directly sequenced using next generation sequencing technology, allowing direct quantification of the number of transcripts in a sample. this process will generate a sequence tag from transcripts that contain nlaiii sites, preferentially from the 3' end of the transcript. this is typically either in the 3' untranslated regions or the 3' end of the open reading frame. transcription profiling using the sage and more recently the dge method has been used extensively to compare expression profiles in a range of tissues. the advantage of dge is the larger dynamic range obtained per experiment. however, there are a number of disadvantages with the dge approach. it does not generate the added value of determining 5' and 3' utrs or alternative splice variants, which is an important consideration for the analysis of most organisms,. dge also requires a reference transcriptome for comparison and is confined to the analysis of genes containing the four-base restriction site used in the library construction. if the 5' and 3' utrs of genes are not included in the reference transcriptome, transcripts with nlaiii sites in the 3'utrs or within bp of the stop codon will not be represented. despite only obtaining a partial transcriptome of the organism under investigation, dge is a useful tool in the analysis of gene expression at an unprecedented level of sensitivity, particularly when comparing two very similar samples. this sensitive and cost-effective approach can be applied to the study of gene expression in cells under a spectrum of experimental or natural conditions. one group of organisms that is subjected to dramatic environmental changes throughout their life cycle, including larges changes in temperature, nutrients and host immune defences is the parasitic protozoa. the evolutionarily ancient single-celled parasite trypanosoma brucei gambiense has a particularly complicated life cycle, passing between the tsetse vector and different mammalian hosts species including humans, where it causes a fatal disease in humans, african sleeping sickness. throughout its life cycle, the parasite has to adapt to extreme changes in its environment, in response to which it differentiates into several morphologically distinct forms, involving organelle repositioning and changes to its metabolism. to what extent are these gross morphological and metabolic changes controlled by alterations in the transcriptome of the parasites between different life cycle stages? studies to date investigating the transcriptome of african trypanosomes have focused mainly on the non-human infective sub-species, t. b. brucei, revealing that transcription in trypanosomes is unusual. the majority of trypanosome genes are transcribed by rna polymerase ii in long polycistronic units, which are then processed by trans-splicing that adds a 39-nucleotide spliced leader to every mrna to produce monocistronic mature mrnas. however evidence to date suggests that, unlike bacterial operons, most genes in polycistronic units are not functionally related and are not co-regulated. the steady state levels of mrna appear to be determined predominantly post-transcriptionally by mrna degradation controlled by the 3'utr of transcripts. there are notable exceptions to this pattern of gene regulation in that some gene arrays that encode the parasite's variant surface glycoproteins and procyclin expressions sites in bloodstream forms and procyclic forms, respectively, are co-regulated and co-transcribed by rna polymerase i. this raises the question, are there other co-regulated polycistronic units within the genome? if so, when are they expressed and what do they encode? are co-regulated polycistronic units transcribed by rna polymerase i alone? while many genes in the genome appear to be arranged in the same orientation and are presumably transcribed as polycistronic units, several small open reading frames, typically between 150- bp in size, with the opposite orientation have been annotated in the genome as being'unlikely'. do these small orfs represent a novel class of trypanosome genes that are transcribed in the opposite orientation from other genes in the polycistron or are they an artefact of gene prediction software? in order to address these questions an analysis of the transcript abundance of these orfs is required. microarray analysis has been carried out to examine changes in transcript abundance between the mammalian bloodstream forms and the insect procyclic forms of the related parasite t. b. brucei and the results suggest that relatively few parasite genes are differentially regulated between the two life cycle stages. in this study we add to the transcriptome analysis of trypanosomes, demonstrating the use of dge to examine the transcriptome of two life cycle stages of the clinically relevant sub-species, t. b. gambiense.  library production and dge tag annotation in this study we determined and compared the partial transcriptome of two different life cycle stages of the pleomorphic t. b. gambiense strain, stib 386; the procyclic stage, which is not human-infective, and the potentially human-infective bloodstream stage. rna from parasites of both life cycle stages was prepared from exponentially growing cells. three dge libraries were made from replicate cultures for each life cycle stage and sequenced using solexa technology. the three procyclic form libraries and three bloodstream form libraries generated × and × tags, respectively, for the two classes of samples, after filtering for poor quality sequencing scores. the normalized tag counts are presented as supplementary data and all raw data is available from gene expression omnibus accession number gse and will be made available on tritrypdb. the number of tags generated using this methodology is at least 300-fold greater than the number of tags obtained for a t. b. brucei strain, using sage with traditional sequencing methods. the tags were aligned to the fully annotated reference transcriptome of t. b. brucei, treu allowing for a two base pair mismatch for any polymorphisms between the reference genome and the genome of stib as the 5' and 3' utr of most t. b. brucei genes have not be defined, the annotated transcriptome contained only orfs, limiting this analysis to tags that align uniquely to the orf of genes. the reference transcriptome represented the chromosome internal genes and did not extend past the sub-telomeric regions of the chromosomes or include the minichromosomes. the number of tags that were of good quality and aligned uniquely to the reference genome for each library is presented in table providing expression data for genes. however there were a large number of tags that either aligned to multiple sites or no sites within the annotated reference transcriptome, in part reflecting the incomplete nature of the in silico transcriptome and the distribution of nlaiii sites in the genome. it is likely that many of the tags that did not align to the transcriptome actually aligned to the 3'utr of transcripts as genes in the t. brucei genome have relatively long 3'utrs with a median size of nt. the dynamic range of dge spanned four orders of magnitude. the most abundant transcript in all procyclic form libraries was that for the heat shock gene, hsp, with a maximum tag count in procyclic form library a of tags accounting for % of uniquely aligned tags for that sample. however, the tag counts for the majority of genes were low in both the bloodstream and procyclic forms. reproducibility digital gene expression has been shown to have very low technical variability. in order to examine the ability of dge to detect subtle differences between biological replicates, dge libraries generated from rna extracted from three identical cultures grown under the same conditions for each life cycle stage were compared. the number of tags generated for each library was similar. figure shows a series of scatter plots comparing the normalised tag counts for the replicate libraries in each pairwise combination. a high correlation was observed between replicates, indicating a high degree of reproducibility of technical and biological replicates, with an average pearson's correlation coefficient of r = for the procyclic form replicates and r = for the bloodstream form replicates following log transformation. bloodstream form parasites when grown at high density are capable of differentiating into the next life cycle stage, the short stumpy form, which is cell cycle arrested and pre-adapted to life in the tsetse fly. the inability to synchronise pleomorphic cells could result in the presence in some cultures of short stumpy differentiated forms, resulting in reduced reproducibility in bloodstream form cultures compared to procyclic forms. although no short stumpy forms were observed in any of the bloodstream form cultures, and the overall correlation coefficient for bloodstream form replicates was reasonably high, similar to that for procyclic forms, a small proportion of trypanosomes could have begun to differentiate into stumpy forms. indeed, pad a gene previously found to be enriched in stumpy forms over slender bloodstream forms and procyclics, was up-regulated 5-fold in the bloodstream forms in this analysis, although not statistically significant at the false discovery rate. this indicated that some of the bloodstream form cells used in this analysis may have begun the differentiation process to stumpy forms. the variation between replicate cultures, although small, clearly demonstrated the value of performing biological replicates. differentially expressed genes the expressed tags that uniquely aligned to the reference transcriptome generated expression data for genes, approximately 81% of the total number of genes in the annotated genome. using a threshold of average fold change in expression for comparison with previous studies using microarrays, genes were up-regulated in bloodstream forms and only genes up-regulated in procyclic forms, a total of 28% differentially expressed genes, compared to 2% using microarray analysis. however, this difference was less an indication of biological significance than a product of the arbitrary nature of using a fold-change cut-off, which was insensitive to differences in factors affecting statistical power such as technical sensitivity, sample size and variability among biological replicates. in this analysis we have employed a statistically valid method of identifying differentially expressed genes taking into account biological replicates, termed rank products, which ranked genes by degree of differential expression and estimated both gene-specific p-values and a fdr that allowed many hypotheses to be tested simultaneously. intuitively, setting a low, stringent fdr would define a smaller but better supported set of genes. using a highly stringent fdr of, where 10% of identified genes would be expected to be false positives, a total of genes were up-regulated in bloodstream forms compared to procyclics and genes were up-regulated in procyclic forms compared to bloodstream forms. at a slightly less stringent fdr of, genes were up-regulated in bloodstream forms and genes were up-regulated in procyclic forms. the overall difference in regulation is illustrated by a quantile-quantile plot where the observed -log for each gene was plotted against the expected value under the null hypothesis of no difference in gene expression between life cycle stages. the observed values for both bloodstream and procyclic stages rose above the y = x line at low p-values > 1) indicating that there were many more low p-values than expected by chance, i.e. many genes were strongly differentially expressed, with more genes differentially expressed in bloodstream stages than in procyclic stages. throughout this analysis all genes with differential gene expression are referred to as being up-regulated either in bloodstream or procyclics. however this does not imply a regulatory mechanism merely a comparative description that could have been generated by a down-regulation in one of the life-cycle stages. remarkably the fold change of many of these differentially expressed genes was extremely high, ranging from to 634-fold up-regulated in bloodstream forms and to up-regulated in procyclic forms. the distribution of the detected fold change is illustrated in the volcano plot, where the statistical significance of each gene was plotted against fold change. the plot shows that the tags with the highest average differences between life cycle stages also showed the greatest degree of significance. another feature of the plot was the tendency of procyclics to show more significant differences than bloodstream forms at fold changes close to zero, reflecting the greater reproducibility within procyclics leading to higher statistical power relative to the bloodstream group. differentially expressed genes, up-regulated in bloodstream forms at the fdr, genes have been identified as being up-regulated in bloodstream forms. thirty two of these were hypothetical genes of unknown function, four of which have been previously identified in t. b. brucei as being up-regulated in bloodstream forms in a proteome study and one by microarray analysis. the identification of these possible coding sequences as being up-regulated in bloodstream forms indicates that they have been correctly classed as genes in the t. b. brucei genome confirming the gene model used for their identification. as expected, a large proportion of genes that were differentially regulated at the fdr level were associated with antigenic variation; the vsg genes, and related expression site associated genes, both of which were expressed at very high levels in bloodstream forms. the vsgs are the major surface coat proteins unique to bloodstream forms, which are encoded by members of a large gene family, although only one vsg gene is expressed at any one time. the vsg gene repertoire is extremely diverse with very few vsg gene sequences being shared between strains. in this analysis, there were three predominant vsg genes that together represent 66% of the vsg tags that uniquely aligned to the reference transcriptome. they aligned to tb. tb. and tb09v pseudogenes in the treu reference transcriptome. one hundred and two other vsg genes that uniquely aligned the reference transcriptome were also expressed within this cell population at low levels. it is likely that there were many other vsg gene tags present in the library that did not align to the reference transcriptome, due to strain specific nature of repertoire of vsg genes. according to current dogma, only one expression site, and therefore only one vsg gene, is active in any one trypanosome cell, this suggests that there were multiple populations of cells expressing different vsg genes within the bloodstream form cultures analysed. an alternative explanation would be that there were low levels of multiple vsg mrnas being expressed at the one time in one cell. the difference in metabolism between the two different life-cycle stages due to the change in the main energy source, from glucose in the bloodstream to proline in procyclic forms, was reflected in the up-regulation in bloodstream forms of the tht1- gene encoding a glucose transporter,, a gene encoding a putative glycerol uptake protein and the ald gene encoding fructose-bisphosphate aldolase, consistent with previous microarray or proteome studies. other genes that were up-regulated in bloodstream forms include those involved in intracellular protein transport or modification, i.e. the p gene that encodes a lysosomal/endosomal membrane protein, nsf gene that encodes the vesicular-fusion protein and the gpi-plc gene that encodes glycosylphosphatidylinositol-specific phospholipase c. the gene tb. encoding a putative iron/ascorbate oxidoreductase family protein has not been previously identified as being up-regulated in bloodstream forms using microarray analysis, although it did appear to be up-regulated in bloodstream forms at the protein level. although this gene was detected only at a low level in bloodstream forms and not detected at all in any of the procyclic form libraries, it was still possible to identify a highly significant difference in transcript abundance, illustrating the sensitivity of digital tag technology to detect differences even for low-abundance transcripts. similarly, three genes encoding leucine-rich repeat proteins that were up-regulated in bloodstream forms were either not detected at all in the procyclic form libraries or at extremely low levels. other genes that were up-regulated in bloodstream forms at the fdr include genes that encode dynein heavy chain proteins putatively involved in intraflagellar transport, the chaperone protein, dnaj, sumo, involved in post-translational modification, glutathionylspermidine synthetase and the cysteine peptidase, calpain, none of which had previously been shown to be up-regulated in bloodstream forms in t. brucei by microarray analysis. by setting a stringent fdr, a small number of well-supported genes were identified in this study, however this would undoubtedly exclude a number of genuinely differentially expressed genes, some of which were previously identified in other studies as being differentially expressed. for comparison, we have taken these previously identified differentially expressed genes described in the microarray studies of koumandou et al. and brems et al. and compared them, where possible, with the dge data in this analysis. this shows that the vast majority of genes that were previously identified as being up-regulated in bloodstream forms in t. b. brucei, were also found to be up-regulated in the related parasite, t. b. gambiense, albeit with approximately half of these genes having a p-value greater than. a differentially expressed cluster of genes of the genes that were up-regulated in bloodstream forms above the fdr threshold, five contiguous genes were up-regulated between and 90-fold. these genes, located at a convergent strand switch region at an internal site on chromosome were three esags and two procyclin associated genes. although most esags are expressed from telomeric expression sites, there is evidence that some esags can also be expressed from chromosome internal regions of the genome. the identification of two pag genes that were highly expressed in bloodstream forms compared to procyclic forms was somewhat surprising, as related genes in this multi-gene family were more highly expressed in the procyclic forms in t. brucei. indeed, copies of pag and pag found on chromosome have been shown to be expressed from a polycictronic unit containing the procyclic surface coat genes, ep and ep2procyclin, which were highly expressed in procyclic forms. the functions of the pag proteins are unknown but there is strong amino acid similarity to the transferrin binding proteins, esag and esag. in order to verify that these genes were up-regulated in bloodstream forms, rt-pcr was performed for each gene in turn, using primers designed to unique regions of each gene. the rt-pcr clearly confirmed the differential expression of all five genes. rnaseq data, released prior to publication, for the related sub-species t. b. brucei also clearly showed the steady-state rna levels of these genes was far greater in bloodstream forms than procyclics. differentially expressed gene clusters were defined as two or more contiguous genes, which were significantly differentially expressed at the fdr. the gene accession identification number, protein description and average fold change is given for each gene. neighbouring genes within a polycistronic transcription unit are usually regulated independently due to variable mrna stability, largely regulated by sequences in the 3'utr of the mrna. such post-transcriptional regulation of gene expression has become a distinguishing feature of t. b. brucei and the basis of much research. here five consecutive genes were co-regulated and highly expressed in bloodstream form trypanosomes. these genes are situated next to a region where neighbouring polycistronic units meet, termed a strand switch region. interestingly, divergent and convergent strand switch regions being associated with transcription start and stop sites, respectively. for this gene cluster, the strand switch region is convergent and likely to contain transcription stop sites. there are remarkably short intergenic regions between each of the five orf in this transcription unit, ranging in size from bp to bp, which could indicate extremely short 3'utrs of these genes and perhaps the lack of any signals for mrna degradation. such highly co-regulated genes are reminiscent of the bloodstream form specific telomeric expression sites, which are unusual in that they are transcribed by rna polymerase i. other rna polymerase i transcription units include the chromosome internal rrna loci and the ep procyclin and gpeet loci, which are up-regulated in procyclic forms in t. b. brucei and which also contain members of the pag multi-gene family. however, there was no obvious sequences homology between regions upstream of this polycistronic unit and the rna polymerase i promoters of the procyclin and vsg loci. other differentially expressed gene clusters up-regulated in bloodstream forms systematic analysis of three or more contiguous differentially expressed genes revealed four additional differentially expressed gene clusters that were up-regulated in bloodstream forms and one cluster up-regulated in procyclic forms. gene cluster contains four genes that were up-regulated to 122- fold on chromosome this gene cluster is present in a polycistronic unit that is positioned at a convergent strand switch region and encodes a member of the pag gene family, in a similar manner to the differentially expressed gene cluster on chromosome two other differentially expressed gene clusters also contain esag and esag and, in all cases, the differentially expressed genes appear not to be syntenic with other kinetoplastids and most contain hypothetical genes unique to t. brucei. all gene clusters described also show differential mrna levels in the rnaseq data of t. b. brucei, except for cluster which was possibly a sub-species differentially expressed gene cluster. what role these genes play in the parasite's life cycle clearly warrants further investigation. differentially expressed genes, up-regulated in procyclics at the fdr level, genes were identified as being up-regulated in procyclic forms fewer than were up-regulated in bloodstream forms. the most up-regulated gene identified in procyclic form trypanosomes encoded the procyclic stage specific surface antigen pssa-. the genes encoding the major procyclic surface coat protein procyclin, ep and ep, do not contain nlaiii sites within the orf, and so these genes did not appear in the in silico derived transcriptome to which the tags were aligned. however a tag matching the 3'utr of ep was identified in the dataset of unaligned tags and was expressed at extremely high level in procyclic forms, up-regulated 250-fold. a large proportion of the genes up-regulated in procyclic forms were annotated as hypothetical providing some insight into their possible role in trypanosome development. interestingly, several small open reading frames annotated as being unlikely to encode proteins were identified as being differentially expressed and particularly up-regulated in procyclic forms, suggesting the gene annotation is correct. known developmentally regulated mrnas involved in mitochondrial biogenesis and metabolism in procyclic forms that were found to be differentially expressed at the protein level were identified in the present analysis, gmdh and ppdk. while many of the cytochrome oxidase subunit genes were identified as being up-regulated in procyclic forms at the fdr, they fall just below the stringent fdr threshold. other genes up-regulated in procyclic forms include tbts, an amino acid transporter gene, a gene encoding a calpain-like protein and the surface protein gene, cram, which have all been previously identified as being up-regulated in procyclic form trypanosomes. several other genes identified in the present study as being up-regulated in procyclic forms were not identified in previous studies at the mrna level, including a putative udp-gal or udp-glcnac-dependent glycosytransferase gene, although this was found to be up-regulated in procyclics in a comparative proteome study. interestingly there are tandemly repeated copies of this gene, but it is only the last gene in the array that was up-regulated in procyclic forms, the rest being up-regulated approximately three-fold in bloodstream forms. intriguingly, while the other genes in the tandem array share near-identical 3'utr sequences, tb. has a very different 3'utr, consistent with the hypothesis that the 3'utr contains regulatory elements that determine differential gene expression. previous analysis of the tandemly repeated gene cluster of phosphoglycerol kinase gene copies, have indicated that each gene is regulated differently in the polycistronic unit, with pgka being constitutively expressed, pgkb being procyclic form specific and pgkc being bloodstream form specific. analysis of the expression patterns of these genes using dge indicated that very few unique tags aligned to the orf of these genes, partly due to the similarity in gene sequence between all three pgk gene copies reducing the number of informative tags. in addition, the majority of tags for pgkc aligned to the 3'utr, which was not represented in the reference transcriptome. the few unique tags that were generated for these genes did not show any significant difference in tag numbers between life cycle stages. assuming the expression pattern for these genes is the same in t.b. gambiense and t.b. brucei, the analysis of this gene family illustrates that this technique only provided relative expression differences for a portion of the differentially expressed transcripts and did not detect all differentially expressed genes. the most highly expressed gene in procyclic forms was hsp although this gene was also among the top four genes expressed in all bloodstream form libraries, it appeared to be up-regulated in procyclic forms -fold. interestingly, hsp has been previously identified as being up-regulated in procyclic forms of the pleomorphic strain, t. b. brucei but not the monomorphic strain t. b. brucei lister, indicating its differential expression was strain specific. whole genome sequencing of the stib strain revealed an estimated copies of this gene in the genome, which could explain the high expression levels. comparing the previously identified genes that have been described as being up-regulated in procyclic forms in the microarray studies of koumandou et al. and brems et al for t. b. brucei with the dge data in this analysis, many genes appear to be up-regulated in both types of analysis, confirming previous analysis and suggesting similarities in differentiation pathways between the different sub-species. however there were several genes, which gave contradictory results, i.e. they appeared to be down-regulated in procyclic forms, albeit not statistically significantly.  in this study, the partial transcriptomes of bloodstream and procyclic cells of t. b. gambiense were compared using dge, providing expression data for genes, constituting 81% of the genome. however this still leaves 19% the genome with no expression information. some of those genes could be specific to the other life cycle stages not analysed here, such as epimastigotes, while others could be unrepresented due to the limitations of the dge technique. these limitations include only sampling genes with nlaiii sites within their open reading frames, an underrepresentation of genes with long 3'utrs and the inability to align many tags uniquely to the transcriptome. while dge allows transcription profiling of a large proportion of the genes in the genome, a more comprehensive analysis of the transcriptome could be have been obtained by rnaseq analysis. in addition, allowing for a two base pair mismatch between the tag and the reference transcriptome to take account of sequence differences between the t. b. brucei and t. b. gambiense genomes could result in misalignment of the tag to the reference transcriptome and could also introduce a bias in some of the analysis. however, this approach is particularly useful when comparing paired samples, in this case different life cycle stages of the same strain and the concordance between the data generated by dge and that of microarrays, albeit for a different sub-species, indicates that this approach is valid. dge has revealed that a larger proportion of genes in t. b. gambiense are stage-regulated in contrast to two previous microarray studies on t. b. brucei. a recent third microarray study examining gene expression in a time course of parasite differentiation revealed many of the differentially expressed genes identified here were common to all studies. where this study differs from microarray analysis is in the far greater sensitivity of dge compared to hybridization-based techniques, which coupled with its high reproducibility means that a larger number of genes that had not previously been identified as being stage-regulated could be identified as being differentially expressed with a high degree of statistical power. the second reason for the difference in gene expression between the dge analysis and the previous microarray studies is that different strains/sub-species were used. strains of t. brucei can vary naturally in a range of different phenotypes, for example, in terms of virulence, growth, drug resistance and resistance to human serum. undoubtedly such phenotypic differences will be as a result of different gene expression patterns, although strain-specific differences have been relatively under investigated in this species with the majority of research being focused on a single laboratory-adapted line. the increased power to detect strain-specific differences means that dge could be a useful tool to examine naturally occurring variable phenotypes, or strains that have been selected for particular phenotypes, and to understand the pathways involved in these processes. one particular phenomenon revealed by dge, which had not previously been observed from microarray studies, but which was evident in the pre-publication of rnaseq data of t. brucei is the presence of several differentially expressed clusters of genes present on what appears to be the same polycistronic units. these clusters are reminiscent of the rna polymerase i transcribed polycistronic units involved in the expression of vsg and procyclin genes, which encode the parasite's surface coat. many of these differentially expressed gene clusters are located in regions of the genome that are unique to t. brucei, with no synteny with other kinetoplastids. these regions are also associated with strand switch sites. the strand switch regions of differentially expressed gene clusters identified here are convergent and do not obviously correspond to the possible rna polymerase ii transcription start sites previously identified. it is tempting to speculate that these islands of differentially expressed genes are transcribed by the rna polymerase that controls transcription of other differentially expressed genes, i.e. rna polymerase i. indeed four of the gene clusters contain either pags or esags, families of genes that are often transcribed by rna polymerase i in either procyclin or bloodstream expression sites, respectively. however at present it is unclear which rna polymerase transcribes these life cycle specific gene clusters, although understanding the basis of the control of the differentially expressed gene clusters clearly warrants further investigation to determine if this phenomenon is due to transcript stability of mrna or transcription rates. in conclusion we have shown that dge analysis generates a wealth of data, revealing a large number of hitherto unknown differentially regulated genes and providing new insights into this unusual pathogen. 
the evolutionarily ancient parasite, trypanosoma brucei, is unusual in that the majority of its genes are regulated post-transcriptionally, leading to the suggestion that transcript abundance of most genes does not vary significantly between different life cycle stages despite the fact that the parasite undergoes substantial cellular remodelling and metabolic changes throughout its complex life cycle. to investigate this in the clinically relevant sub-species, trypanosoma brucei gambiense, which is the causative agent of the fatal human disease african sleeping sickness, we have compared the transcriptome of two different life cycle stages, the potentially human-infective bloodstream forms with the non-human-infective procyclic stage using digital gene expression analysis. the transcriptome analysis revealed the existence of several differentially expressed gene clusters within the genome, indicating that contiguous genes, presumably from the same polycistronic unit, are co-regulated either at the level of transcription or transcript stability. dge analysis is extremely sensitive for detecting gene expression differences, revealing firstly that a far greater number of genes are stage-regulated than had previously been identified and secondly and more importantly, this analysis has revealed the existence of several differentially expressed clusters of genes present on what appears to be the same polycistronic units, a phenomenon which had not previously been observed in microarray studies. these differentially regulated clusters of genes are in addition to the previously identified rna polymerase i polycistronic units of variant surface glycoproteins and procyclin expression sites, which encode the major surface proteins of the parasite. this raises a number of questions regarding the function and regulation of the gene clusters that clearly warrant further study. over eleven million unique tags were generated, producing expression data for genes, covering 81% of the genes in the genome. compared to microarray analysis of the related t. b. brucei parasite, approximately times more genes with a -fold change in expression levels were detected.
 the evolutionarily ancient parasite, trypanosoma brucei, is unusual in that the majority of its genes are regulated post-transcriptionally, leading to the suggestion that transcript abundance of most genes does not vary significantly between different life cycle stages despite the fact that the parasite undergoes substantial cellular remodelling and metabolic changes throughout its complex life cycle. to investigate this in the clinically relevant sub-species, trypanosoma brucei gambiense, which is the causative agent of the fatal human disease african sleeping sickness, we have compared the transcriptome of two different life cycle stages, the potentially human-infective bloodstream forms with the non-human-infective procyclic stage using digital gene expression analysis.  over eleven million unique tags were generated, producing expression data for genes, covering 81% of the genes in the genome. compared to microarray analysis of the related t. b. brucei parasite, approximately times more genes with a -fold change in expression levels were detected. the transcriptome analysis revealed the existence of several differentially expressed gene clusters within the genome, indicating that contiguous genes, presumably from the same polycistronic unit, are co-regulated either at the level of transcription or transcript stability.  dge analysis is extremely sensitive for detecting gene expression differences, revealing firstly that a far greater number of genes are stage-regulated than had previously been identified and secondly and more importantly, this analysis has revealed the existence of several differentially expressed clusters of genes present on what appears to be the same polycistronic units, a phenomenon which had not previously been observed in microarray studies. these differentially regulated clusters of genes are in addition to the previously identified rna polymerase i polycistronic units of variant surface glycoproteins and procyclin expression sites, which encode the major surface proteins of the parasite. this raises a number of questions regarding the function and regulation of the gene clusters that clearly warrant further study. 