Background: Leishmania (L) are intracellular protozoan parasites that are able to survive and replicate within the harsh and potentially hostile phagolysosomal environment of mammalian mononuclear phagocytes. A complex interplay then takes place between the macrophage (MPhi) striving to eliminate the pathogen and the parasite struggling for its own survival. To investigate this host-parasite conflict at the transcriptional level, in the context of monocyte-derived human MPhis (MDM) infection by L. major metacyclic promastigotes, the quantitative technique of serial analysis of gene expression (SAGE) was used.

Results: After extracting mRNA from resting human MPhis, Leishmania-infected human MPhis and L. major parasites, three SAGE libraries were constructed and sequenced generating up to 28,173; 57,514 and 33,906 tags respectively (corresponding to 12,946; 23,442 and 9,530 unique tags). Using computational data analysis and direct comparison to 357,888 publicly available experimental human tags, the parasite and the host cell transcriptomes were then simultaneously characterized from the mixed cellular extract, confidently discriminating host from parasite transcripts. This procedure led us to reliably assign 3,814 tags to MPhis' and 3,666 tags to L. major parasites transcripts. We focused on these, showing significant changes in their expression that are likely to be relevant to the pathogenesis of parasite infection: (i) human MPhis genes, belonging to key immune response proteins (e.g., IFNgamma pathway, S100 and chemokine families) and (ii) a group of Leishmania genes showing a preferential expression at the parasite's intra-cellular developing stage.

Conclusion: Dual SAGE transcriptome analysis provided a useful, powerful and accurate approach to discriminating genes of human or parasitic origin in Leishmania-infected human MPhis. The findings presented in this work suggest that the Leishmania parasite modulates key transcripts in human MPhis that may be beneficial for its establishment and survival. Furthermore, these results provide an overview of gene expression at two developmental stages of the parasite, namely metacyclic promastigotes and intracellular amastigotes and indicate a broad difference between their transcriptomic profiles. Finally, our reported set of expressed genes will be useful in future rounds of data mining and gene annotation.

Figure 1

Figure 1

Venn diagram comparing the parasite-infected...

Figure 1

Venn diagram comparing the parasite-infected MDM and L. major SAGE libraries with the...

Figure 2

Figure 2

Comparison of gene expression modulation...

Figure 2

Comparison of gene expression modulation in the L. major -infected "MDM+Lm" library to...

Figure 3

Figure 3

Covariance Analysis of Gene expression...

Figure 3

Covariance Analysis of Gene expression profiles raised with the 500 most abundant tags...

Figure 4

Figure 4

Examples of gene transcripts categorized...

Figure 4

Examples of gene transcripts categorized into functional classes involved in defense MΦ programs ...