article: Genome sequencing of Leishmania species that give rise to a range of disease phenotypes in the host has revealed highly conserved gene content and synteny across the genus. Only a small number of genes are differentially distributed between the three species sequenced to date, L. major, L. infantum and L. braziliensis. It is not yet known how many of these genes are expressed in the disease-promoting intracellular amastigotes of these species or whether genes conserved between the species are differentially expressed in the host. We have used customised oligonucleotide microarrays to confirm that all of the differentially distributed genes identified by genome comparisons are expressed in intracellular amastigotes, with only a few of these subject to regulation at the RNA level. In the first large-scale study of gene expression in L. braziliensis, we show that only ∼9% of the genes analysed are regulated in their RNA expression during the L. braziliensis life cycle, a figure consistent with that observed in other Leishmania species. Comparing amastigote gene expression profiles between species confirms the proposal that Leishmania transcriptomes undergo little regulation but also identifies conserved genes that are regulated differently between species in the host. We have also investigated whether host immune competence influences parasite gene expression, by comparing RNA expression profiles in L. major amastigotes derived from either wild-type (BALB/c) or immunologically compromised (Rag2−/− γc−/−) mice. While parasite dissemination from the site of infection is enhanced in the Rag2−/− γc−/− genetic background, parasite RNA expression profiles are unperturbed. These findings support the hypothesis that Leishmania amastigotes are pre-adapted for intracellular survival and undergo little dynamic modulation of gene expression at the RNA level. Species-specific parasite factors contributing to virulence and pathogenicity in the host may be limited to the products of a small number of differentially distributed genes or the differential regulation of conserved genes, either of which are subject to translational and/or post-translational controls.
  Infection with species of the kinetoplastid parasite, Leishmania, results in a spectrum of diseases in man, termed the leishmaniases [1], [2]. These range from the non-fatal chronic cutaneous lesions arising from L. major infection to mucocutaneous leishmaniasis usually associated with L. braziliensis (classified within the sub-genus L. Viannia) and the often fatal visceralising disease, most commonly associated with L. donovani infection in the Indian sub-continent, L. chagasi in Brazil and L. infantum in the Mediterranean basin. (The last two species are generally considered to be genetically identical [3]). While the species of infecting parasite can play a defining role in disease type, the genetic background and immune response of the host are also major factors in determining clinical outcome [4], [5], [6], [7], [8], [9]. Understanding the relative contribution of these different components may enhance our understanding of pathogenicity in the leishmaniases. Sequencing and comparison of the genomes of representative lab-adapted strains of L. major, L. infantum and L. braziliensis have revealed strong conservation of gene content and synteny, with only a small number of genes identified as differentially distributed between species [1], [10]. This subset of genes, together with sequences preferentially expressed in intracellular amastigotes and/or showing differential expression between species, may be important in facilitating parasite survival and maintenance within the host. The best-characterised example of the former class is the L. donovani complex-specific A2 gene coding for an amastigote protein of as yet unknown function which, when expressed in L. major, leads to increased parasite dissemination to the viscera [11]. Recent expression profiling has identified 3–9% of genes that are modulated at the RNA level between life cycle stages of several Leishmania species [12], [13], [14], [15]. Moreover, comparisons of L. mexicana amastigote parasites grown axenically with those maintained within macrophages, either in vitro or in vivo, have shown that axenic (extracellular) amastigotes are more similar to extracellular promastigotes than to macrophage-derived (intracellular) amastigotes in their RNA profiles [12]. These observations emphasise the importance of using parasites isolated ex vivo to investigate the mechanisms of intracellular survival. To date, no comparative expression profiling has been performed on L. Viannia spp. , despite the relative divergence of the genome of L. braziliensis from that of L. major or L. infantum [1]. A complicating factor in the analysis of Leishmania gene expression is the almost complete absence of defined RNA polymerase II promotors in kinetoplastid species, coupled with the characteristic bidirectional polycistronic transcription units found on individual chromosomes [16], [17], [18]. In these organisms, polycistronic precursor RNAs (which may be expressed constitutively) are processed by coupled trans-splicing and polyadenylation [19] to generate mature mRNA transcripts for translation. Expression of individual genes is regulated post-transcriptionally, a process largely dependent on RNA stabilisation mechanisms [17], [18]. This post-transcriptional regulation, coupled with an as yet unknown extent of regulation at the level of translation, results in variable correlations between gene and protein expression levels [20], [21]. Such factors place greater emphasis on identifying those genes which undergo regulation during the life cycle while also looking for differences in the relative levels of expression of conserved genes. Given the importance of intracellular Leishmania stages to disease in man, this study has focused on amastigote gene expression, comparing RNA expression profiles between the three sequenced Leishmania species to identify any significant differences that may be functionally relevant in these infective parasite stages. To achieve this end, an oligonucleotide array was designed representative of 4 functional classes of genes, together comprising ∼10% of the genome. These targets included (i) all genes identified as differentially-distributed between the 3 sequenced Leishmania species [1]; (ii) all genes containing amino acid repeats within their open reading frames [22]; (iii) all genes encoding proteins predicted (with high probability) to be co-translationally modified by N-myristoylation [23]; (iv) a range of control genes of known RNA expression profiles. The rationale for this choice of target genes is that their expressed products provide a subset for analysis as putative targets for therapeutic intervention. 2–4% of the 3 sequenced Leishmania genomes, including 30% of the differentially-distributed genes, contain amino acid repeats within their protein-coding regions, as identified by the RepSeq web utility [22]. Repeat domains are associated with a range of functions relevant to host survival in other protozoan parasites, including antigenic variation, host-cell receptor binding and intracellular protein-protein interactions in Plasmodium species [24], [25], [26]. Amino acid repeat-containing proteins encoded within the Leishmania genomes include several kinases, cysteine peptidases, putative surface antigen proteins and the infective stage-specific HASPB (formerly named GBP), expressed in amastigotes of all L. Leishmania species analysed but absent from L. (Viannia.) braziliensis [1], [27], [28], [29], [30]. HASPB is also a target for N-myristoylation, a co-translational protein modification catalysed by the enzyme N-myristoyltransferase (NMT) that is a genetically-validated target for Leishmania drug development [31], [32]. Hybridisation of these customised oligonucleotide arrays with amastigote RNA derived from footpad lesions (L. major), spleens (L. infantum) and RAW 264. 7 macrophages (L. braziliensis) revealed that only a small number of the target genes are differentially expressed at the mRNA level between species and that fewer still are regulated during the parasite life cycle. In addition, comparison of expression profiles derived from L. major amastigotes isolated from hosts of differing immune competence (BALB/c vs. Rag2−/− γc−/− mice [33]) demonstrated that host immune pressure has little effect on parasite gene expression at the RNA level. Taken as a whole, the data presented here suggest that parasites do not respond dynamically to host immune pressure and that any influence of varying transcript levels on virulence and pathogenicity of different Leishmania species is likely to result from the differential expression of conserved genes between species and/or the expression of a small number of genes that are differentially distributed between species.
  The three representative Leishmania proteomes (L. major, L. braziliensis and L. infantum) were analysed for the presence of amino acid repeats using default settings at the RepSeq web utility (www. repseq. org; [22]); 256 proteins were identified. A further 62 proteins were predicted to be N-terminally myristoylated with high confidence using the NMT Predictor and Myristoylator programs as described [23]. The genes identified as differentially distributed between the three Leishmania species [1] were reanalysed using BLAST and the latest assemblies of the three genomes available at the time on GeneDB (L. major published assembly [10], L. infantum v2 and L. braziliensis v1). 242 of these genes were included in this study (Table S1); the remainder have been removed during annotation revision. The three sequenced genome reference strains used were L. major MHOM/IL/80/Friedlin; L. infantum clone JPCM5 MCAN/ES/98/LLM-877; L. braziliensis clone M2904 MHOM/BR/75M2904. Procyclic and metacyclic parasites were maintained at 26°C in pH 7. 0 and pH 5. 5 media respectively; for L. major and L. braziliensis, 1×M199 medium was supplemented with 10% Gibco fetal bovine serum (FBS, Invitrogen) and penicillin-streptomycin (Invitrogen) and, in the case of L. braziliensis, 2% human male urine. L. infantum was cultured in Gibco HOMEM (Invitrogen) supplemented with 10% FBS (Biosera) and penicillin-streptomycin. The pH of all media was adjusted using orthophosphoric acid. Metacyclic RNAs were validated by RT-qPCR profiling of the metacyclic specific markers SHERP [34] in L. major and L. infantum and Meta1 [35] in L. braziliensis. For L. major in vivo passage, BALB/c and Rag2−/−γc−/− mice [33], [36] were infected sub-cutaneously in the rear footpad with 4×106 late stage L. major metacyclics purified by Percoll gradient fractionation [37]. Footpad thickness was measured weekly using digital calipers, taking an average of three readings. At four weeks post-infection, mice were euthanized and tissues were sampled for analysis. For L. infantum in vivo passage, infection was initiated in hamsters by intraperitoneal inoculation of 107 stationary phase metacyclics as described [38], with euthanasia and spleen removal 9–12 months later. BALB/c mice (Charles River Laboratories) and Rag2−/−γc−/− mice (bred in the Centre for Immunology and Infection) were housed in pathogen-free conditions at the University of York. All animal work was conducted under UK Home Office Licence requirements and after institutional ethical review. Attempts to purify amastigotes from the footpads of BALB/c, C57BL/6 and Rag2−/−γc−/− mice infected with the L. braziliensis genome strain were unsuccessful (despite establishment of lymph nodes infections with the same parasites), so Percoll-purified metacyclics were used to infect RAW 264. 7 macrophages maintained in DMEM medium supplemented with 10% FBS, l-glutamine and sodium pyruvate (1 mM final concentrations). Amastigotes were obtained from infected macrophages after 3 days culture at 34°C. For all species, infections were initiated with parasites that had been recently passaged in vivo to reduce any impact of long term in vitro culture on parasite gene expression and virulence. Infected tissues were homogenised and fractionated (centrifugation at 100 g, 4°C) to remove larger debris. Host cells in the resulting supernatant were pelleted (2000 g, 4°C) and treated with 0. 05% w/v saponin (Sigma Aldrich) for lysis and release of amastigotes. Host cell debris was removed by layering parasites over Percoll (1. 037 g/ml) and centrifugation at 2000 g, 4°C for 1 hr. Amastigote yields were assessed using a haemocytometer and Beckman Z Series Coulter Counter. Parasite species was verified by restriction fragment length polymorphism analysis of amastigote DNA [39]. For immunoblotting, 5×106 host-derived amastigotes were pelleted, washed in PBS and resuspended in an SDS load mix supplemented with Pepstatin A (50 ug/ml). Treatment with 3% β-mercaptoethanol was followed by a 10 min incubation at 95°C prior to 12% SDS-PAGE, electroblotting on immobilon-P membrane (Millipore) and probing with anti-HASPB [40], [41] at 1∶1250 dilution and anti-NMT [31] at 1∶1000, as previously described. ECL Western Blotting Detection Reagents (GE Healthcare Lifesciences) were used in conjunction with a G∶BOX Chemi Imaging System and GeneSnap V. 7. 0 software (Syngene) to determine relative expression of detected proteins. L. major parasite burdens in mouse footpad lesions were determined by purification and manual counting while those in the liver and spleen tissues were assessed from duplicate samples using quantitative PCR, adapted from the method of [42]. DNA was extracted from tissues using a DNeasy Blood and Tissue kit (Qiagen) following overnight digestion with Proteinase K at 55°C. Quantitation of L. major DNA was performed using primers specific for a 116-bp fragment of the kinetoplastid minicircle DNA: forward, 5′-CCTATTTTACACCAACCCCCAGT-3′ (JW11) and reverse, 5′-GGGTAGGGGCGTTCTGCGAAA-3′ (JW12) [43]. DNA standards were created for each tissue sampled, by spiking naïve tissue with known numbers of L. major promastigotes and extracting the DNA as described above. For RT-qPCR reactions, samples were analysed using Power SYBR Green PCR Master Mix and an ABI 7300 sequence analyzer (Applied Biosystems). Samples were subjected to an initial holding step at 95°C for 10 min, followed by an amplification step (40 cycles of 95°C for 15 sec and 60°C for 1 min) and a single dissociation step of 95°C for 15 sec, 60°C for 20 sec and 95°C for 15 sec. The total parasite burden per organ was calculated by determining the total number of parasites per g of tissue sampled and multiplying by total organ mass. Amastigote pellets were resuspended in TRIzol Reagent (Invitrogen) at a concentration of 1 ml/108 cells and total RNA extracted using the protocol provided. Extracted RNA was purified using a MEGAclear kit (Ambion) and the quantity and quality verified using a NanoDrop ND-1000 Spectrophotometer and Agilent 2100 Bioanalyzer, respectively, prior to storage at −80°C. cDNA synthesis and subsequent labelling was performed using the Amino Allyl cDNA labeling kit (Ambion) and Alexa Flour 555 and 647 dyes (Invitrogen). 70-mer oligonucleotides were designed and synthesised by Operon Biotechnologies using parameters such that a single oligo had >90% identity with the gene target in each of the three species but <70% identity with any other gene on the array. All oligos (Table S2) were BLAST-searched against the NCBI databases for the mouse and hamster genomes, to minimise potential cross hybridisation. The microarrays were printed on Nexterion Slides E (Schott) using a QArrayMini (Genetix) as follows: 5 replicates were printed for each probe in a random pattern totalling 5376 probes per slide (including buffer and negative controls). cDNA from procyclic RNA from each of the three species was hybridised against the arrays for quality control, revealing that all probes were hybridised (data not shown). Prior to hybridisation, each slide was washed at room temperature as follows: 1×5 min in 0. 1% Triton X-100,2×2 min in 1 mM HCL, 1×10 min in 100 mM KCL. Slides were subsequently blocked by incubation in 1×Nexterion Blocking E (Schott) for 15 min and then dried by centrifugation. Labelled cDNA pellets were resuspended in Nexterion Hyb and hybridised on the slides in Corning hybridisation chambers (Fisher) for 16 hrs at 45°C. The slides were then washed in 2×Sodium Chloride Citrate (SSC) /0. 1% sodium dodecyl sulphate (SDS) for 5 min, 2×SSC for a further 5 min, 0. 2×SSC for 5 min and dried by centrifugation at 200 g for 3 min. The slides were scanned using an Axon 4000A scanner and fluorescence signal intensities of the array features and local backgrounds measured by GenePix Pro 5. 0. For each slide, flags were automatically assigned to each probe and manually evaluated, based on the quality of hybridisation, to assist in normalisation within and between arrays. All raw data extracted from GenePix were imported into R 2. 2. 1 for normalization and statistical analysis using the Bioconductor packages (www. bioconductor. org), predominantly Limma. Briefly, gpr files were read and background correction performed using the Kooperberg model. QC analysis, entailing density distribution analysis coupled with correlation distance heatmap plotting, was used to identify substandard arrays, which were removed from further analysis. Array data were normalised within arrays using printtiploess, followed by normalisation between arrays using the “Aquantile” method, which normalises intensity distribution across all slides. Subsequently, the mean intensity was determined for the 5 replicate spots on each array prior to single channel linear modelling and contrast fitting of non-control array features, with eBayes methods to generate moderated t-statistics. Four contrasts were performed: L. braziliensis procyclics vs. metacyclics vs. amastigotes, L. major amastigotes vs. procyclics, L. infantum amastigotes vs. procyclics and L. major amastigotes from BALB/c mice vs. L. major amastigotes from Rag2−/−γc−/− mice. Gene lists were generated, with accompanying statistics, from the eBayes fit data. The decideTests method was applied, with an adjust p value threshold of 0. 01, to classify genes as to their significance in one of more contrasts. A fold change of 1. 7 was then used as threshold for determining differentially expressed genes, similar to other recent studies [13], [14], [20]. In accordance with MIAME (Minimum Information About a Microarray Experiments) regulations [44], all data were deposited into ArrayExpress database [45] at www. ebi. ac. uk under the accession number E-MEXP-2063. For host-derived amastigote analysis, a minimum of 3 biological replicates (each comprising 2–3 technical replicates) per species were used. 2–3 technical replicates were also used for each of the following analyses: culture-derived amastigotes, metacyclics and procyclics RT-qPCR was performed on selected genes to validate the microarray data, using 1–2 biological replicates (independent of those used on the arrays). cDNA was synthesised using the Omniscript RT Kit (Qiagen). Cross-species primers were designed by aligning gene sequences using CLUSTALW (http: //align. genome. jp/) and then selecting primers based on the output of Primaclade [46]. Each biological sample was run with 4–5 technical replicates. Data generated were analysed [47] and normalised using a range of suitable targets (including published controls; [14], [35], [48]). All primers used are shown in Table S3.
  Expression profiles were generated for each L. braziliensis life cycle stage using mRNA from culture-derived procyclics and metacyclics and RAW 264. 7 macrophage-derived amastigotes. The profiles were analysed using scatter plots (Figure 1A–C), revealing that ∼9% (60/678) of the genes probed were significantly regulated (fold change >1. 7, p-value<0. 05) in their expression between life cycle stages. Thus, although the arrays used for this analysis of L. braziliensis gene expression represent only 10% of the genome (and are biased towards specific gene subsets), the proportion of stage-regulated genes (∼9%) is of the same order of magnitude as that observed in whole genome analysis of L. major (7%), L. mexicana (3. 5%) and L. infantum (9. 3%) [12], [14]. Of the 60 regulated genes in L. braziliensis (Table 2, Figure 2), 5 were preferentially expressed in amastigotes (LbrM05_V2. 0380 - microtubule-associated protein, LbrM07_V2. 0360 - ATP-dependent DEAD/H RNA helicase, LbrM09_V2. 0960 - calmodulin, LbrM11_V2. 0560 - 40S ribosomal protein S21 and LbrM31_V2. 2570 - 3′ nucleotidase/nuclease) with a further 35 preferentially expressed in metacyclics and 11 in procyclics. Analysis of the specific gene subsets identified none of the NMT targets as upregulated in amastigotes while 5 of these genes showed preferential expression in either procyclics and/or metacyclics. 9. 8% (25/256) of the genes coding for proteins containing amino acid repeats showed stage specific regulation, with 3 of these encoding proteins of unknown function that are preferentially expressed in the procyclic and metacyclic stages (as compared with amastigotes). Interestingly, all of the genes preferentially expressed in L. braziliensis amastigotes in this analysis, with the exception of the calmodulin array, contain repeat motifs. If these RNA differences correlate with protein levels, this result could indicate that repeat-containing proteins are functionally important in the intracellular stage of the life cycle and provide a useful sub-set of biomarkers with diagnostic and vaccination potential. Pair-wise comparisons of amastigote RNA profiles for each of the three species identified only a small number of conserved genes showing differential expression. Varying numbers of probes were directly comparable between the species (as some genes are present in two rather than three species) including 608 between L. braziliensis and L. major, 615 between L. braziliensis and L. infantum and 686 between L. major and L. infantum. The expression profiles are shown in Figure 1D–F, with the number of differentially expressed genes in each species summarised in Figure 3 (with details in Table S4). These analyses were used to identify genes showing significantly increased amastigote expression in one species compared with the other two (Table 3). Five genes were identified as showing significantly increased expression in L. braziliensis amastigotes: 2 genes with predicted functions (LbrM08_V2. 0810 - cathepsin L-like protease, LbrM25_V2. 1020 and LbrM35_V2. 3370 - ADP-ribosylation factor) and 3 genes of unknown function (LbrM06_V2. 0730, LbrM20_V2. 0550 and LbrM21_V2. 0910). A further 5 genes showed significantly increased expression in L. major amastigotes, 4 of which code for proteins of predicted functions (LmjF06. 0400 - Fructose biphosphate aldolase, LmjF24. 1730 - protein kinase, LmjF31. 0830 - Triacyglycerol lipase, LmjF36. 5370- tyrosine specific protein phosphatase) while the last codes for an unknown protein (LmjF31. 1400). A total of 15 genes showed significantly increased expression in L. infantum amastigotes alone, including 10 coding for proteins of unknown function (LinJ05_V3. 0670, LinJ07_V3. 0270, LinJ07_V3. 0950, LinJ10_V3. 1230, LinJ13_V3. 0770, LinJ14_V3. 0470, LinJ18_V3. 1160, LinJ24_V3. 1530, LinJ29_V3. 0110 and LinJ31_V3. 1470). and 5 for proteins of predicted function (LinJ12_v4. 0671 – surface antigen protein, LinJ04_V3. 0320 - mitochondrial exoribonuclease DSS-1, LinJ11_V3. 0400 - tubulin-tyrosine ligase-like protein, LinJ31_V3. 1490 - surface membrane protein gp46-like protein and the proteophosphoglycan ppg3 array on chromosome 35 (chromosome 34 in L. braziliensis). Interestingly, gp46 has also been identified as an upregulated gene in L. donovani amastigotes derived from patients with PKDL (post-kala azar dermal leishmaniasis; [52]). The ppg3 gene is located within an array of related sequences (still undergoing refinement of annotation due to their repetitive structure) on chromosome 35 of Leishmania species and codes for one of the glycoprotein constituents of filamentous proteophosphoglycan (fPPG), a major component of promastigote secretory gel (PSG). PSG is produced by parasites within the gut of the sandfly vector, where it contributes to a physical blockage that promotes host blood-feeding and subsequent parasite transmission [53], [54]. The Leishmania PPGs belong to a novel class of serine- and threonine-rich proteins possessing large conserved Ala-Pro-Ser repeats that represent major antigenic determinants in all Leishmania species examined to date [53], [55], [56]. These repeats are extensively modified by phosphodiester-linked oligosaccharides and terminal manno-oligosaccharides, while the modified PPGs are either water soluble and secreted or membrane-bound. PPGs are expressed in both promastigotes and amastigotes, localising to the flagellar pocket, endosomes, lysosomes and the external surface of the parasite [57], [58], [59]. Surface PPGs are important in parasite-macrophage interactions: they apparently bind to host cells to modulate macrophage function during early infection, both by inhibiting tumor necrosis factor-α production and synergizing with interferon-γ to stimulate the production of nitric oxide [60]. The significance of increased RNA expression from the PPG gene array in L. infantum amastigotes (as compared with L. major and L. braziliensis amastigotes) is unknown; no studies have yet investigated the L. infantum PPGs or addressed whether one or more of these molecules might perform a species-specific role in establishing and maintaining amastigotes within host macrophages. Calmodulin, the ubiquitous eukaryotic intracellular calcium receptor, plays a role in the regulation of many cellular proteins and transmembrane ion transporters with wide-ranging downstream physiological consequences. Calmodulin is encoded by a 3-gene array on chromosome 9 in Leishmania but has not yet been widely studied in this organism, despite linkage to the modulation of a plasma membrane Ca2+-ATPase in L. donovani [61]. In this present study, both the array analyses and subsequent RT-qPCRs show decreased abundance of calmodulin transcripts in amastigotes of L. major and L. infantum as compared with procyclic and metacyclic parasites. This correlates with the observed changes in calmodulin protein during in vitro differentiation of L. donovani [62]. By contrast, calmodulin transcript levels are upregulated in L. braziliensis amastigotes (Figure 4). To date, it is unclear whether these RNA differences reflect mRNAs derived from one or several of the gene copies, nor how this distribution might vary between the species. It is also unknown whether the varying RNA levels correlate with protein expression in L. braziliensis and how this might influence calmodulin function in different amastigote populations. The small number of differentially distributed genes identified by comparative genomic analysis of L. major, L. infantum and L. braziliensis [1] include several well-characterised stage specific genes that are implicated in vector transmission and virulence in the host (e. g. A2, HASPB, SHERP [11], [34], [63]. The expression profiles for the ∼200 differentially distributed genes were determined from the L. braziliensis life cycle analysis generated here and supplemented with published microarray data generated from life cycle studies of L. major and L. infantum [12], [14] (Table S1). These cumulative results have confirmed that all genes differentially distributed between the three species are expressed in all life cycle stages. Only 34 of these sequences are regulated between stages, however: 19 in L. major, 9 in L. infantum and 6 in L. braziliensis. The majority of these genes encode proteins of unknown function while the remainder are the HASPA, HASPB and SHERP genes from L. Leishmania species as well as an adenine phosphoribosyltransferase (LbrM26_V2. 0120) and a unique gene family in L. braziliensis (LbrM23_V2. 1110 and LbrM23_V2. 1120) and the L. infantum specific LinJ22_V3. 0670. Previous comparisons of the life cycle expression profiles of L. major, L. infantum and L. mexicana have revealed similar numbers of genes that are regulated during the life cycle [12], [13], [14]. The most interesting observation associated with these data is that, while these regulated genes are conserved between all three species, whether they are regulated or not appears to be a species-specific property. As an example of this, only 114/1228 of the conserved regulated genes of L. major and L. infantum undergo modulation in both species [14]. A similar effect has been observed in comparisons between L. major and L. mexicana promastigote RNA expression profiles [12]. In this current study, of the 21 genes identified as being stage regulated in L. braziliensis (Table S5), 14 are only regulated in this species while 6 share regulation with either L. major or L. infantum and one, a 3′-nucleotidase/nuclease, is regulated in all three species and shows upregulation in amastigotes [12], [14]. Leishmania parasites are purine auxotrophs and require mechanisms to salvage these essential nutrients from the sandfly gut and the host [64], [65], [66]. The 3′-nucleotidase/nuclease [67] degrades and dephosphorylates exogenous purine sources to nucleosides that can be transported across the plasma membrane [67]. Independent studies in L. donovani [68] and L. mexicana [69] have shown that 3′-nucleotidase/nuclease mRNA transcripts are upregulated in promastigotes when parasites are starved of purines but not in amastigotes. One interpretation of these conflicting observations could be that there are additional gene copies expressing nucleotidase/nuclease activity during the parasite life cycle. The complete list of all regulated genes identified in L. braziliensis, L. major and L. infantum is presented in Table S5. The influence of host immune pressure on gene expression in Leishmania has not been studied in vivo to our knowledge. In this study, parasite burdens, lesion development, and parasite gene expression profiles were examined in immunocompetent BALB/c and immunocompromised Rag2−/−γc−/− mice which have a mixed genetic background of 129 Ola, BALB/c and C57BL/6 [33]. These latter mice possess a mutation in the common γ-chain receptor (a component of the receptors for IL-2, -4, -7, -9 and -15) and in the recombination activating gene 2 (Rag2; required for V (J) D rearrangement) and are characterised by a complete absence of mature T lymphocytes, B lymphocytes and NK cells [33]. Footpad infections were established by inoculation with 4×106 metacyclic L. major parasites and resultant lesion formation was measured at weekly intervals. The Rag2−/−γc−/− mice displayed no signs of inflammation in the footpad and did not develop cutaneous lesions. In contrast, vigorous lesion development was seen in BALB/c mice with a greater than 2-fold increase in footpad thickness at one month post-infection (Figure 5A). Parasite burdens in the footpads of Rag2−/−γc−/− mice were significantly lower than in the BALB/c mice; however these immunocompromised mice showed 2. 5-fold higher parasite burdens in the liver and spleen (Figure 5B–D). Enhanced parasite burdens in visceral organs may arise from accelerated dissemination of parasites through the lymphatic system from the site of infection, as Rag2−/−γc−/− mice lack peripheral lymph nodes. In addition, survival and proliferation of parasites may be increased in these tissues due to the inability of Rag2−/−γc−/− mice to mount protective T and B cell dependant responses required for parasite clearance. Expression profiling of the amastigotes taken from the footpads of these mouse strains revealed only 3 genes that showed a significant difference in expression (out of 700 genes on the array). Only one of these genes (LmjF06. 0720, coding for a protein of unknown function) was preferentially expressed in amastigotes from BALB/c mice. The other 2 genes showed increased expression in the Rag2−/−γc−/− mice: LmjF24. 1840, a lysophospholipase and LmjF35. 4730, another “unknown” gene (Table 4). Independent analysis of HASPB transcript and protein abundance in these amastigotes (by RT-qPCR and immunoblotting respectively) revealed no significant differences in expression (Figure 6). This lack of amastigote RNA modulation, assuming that it correlates with protein expression levels, confirms that Leishmania parasites are pre-adapted toward intracellular survival, regardless of the state of the host immune system. These data are of relevance to our interpretation of the cross-species amastigote expression results (Table 3) given that the amastigotes were sourced from rodents of different genetic backgrounds. Small microarrays, targeted towards 9% of the genes conserved between the three representative Leishmania species and an additional ∼250 genes differentially distributed between these species, were designed and hybridised with host-derived amastigote RNA isolated from cutaneous lesions (L. major), spleens (L. infantum) and RAW 264. 7 macrophages (L. braziliensis). ∼9% of the probed genes were identified as being regulated during the L. braziliensis life cycle, a figure comparable to that observed using whole genome arrays for L. major, L. mexicana and L. infantum [12], [14]. It is interesting to note however, that while most regulated genes are conserved between all three representative species, the majority of genes regulated in one species are not regulated in the others. Moreover, comparative expression profiles generated for L. major, L. infantum and L. braziliensis amastigotes revealed that species-specific differential regulation of conserved genes was common and this may impact on parasite survival in the host. These species-specific differences require further study that should focus on determining whether relative protein abundances are affected by increased transcript abundances. Attempts to correlate the relative abundance of mRNA and protein for individual genes in L. infantum amastigotes resulted in only modest agreement (∼53% [21]), with these differences attributed to the extensive post-transcriptional, translational and post-translational regulation operating in Leishmania parasites [13], [21]. Finally, the lack of modulation of gene expression profiles in L. major parasites responding to different immune pressures is consistent with the hypothesis that Leishmania parasites are constitutively adapted toward survival in a range of hosts [20].
summary: The single-celled parasite Leishmania, transmitted by sand flies in more than 88 tropical and sub-tropical countries globally, infects man and other mammals, causing a spectrum of diseases called the leishmaniases. Over 12 million people are currently infected worldwide with 2 million new cases reported each year. The type of leishmaniasis that develops in the mammalian host is dependent on the species of infecting parasite and the immune response to infection (that can be influenced by host genetic variation). Our research is focused on identifying parasite factors that contribute to pathogenicity in the host and understanding how these might differ between parasite species that give rise to the different clinical forms of leishmaniasis. Molecules of this type might lead to new therapeutic tools in the longer term. In this paper, we report a comparative analysis of gene expression profiles in three Leishmania species that give rise to different types of disease, focusing on the intracellular stages that reside in mammalian macrophages. Our results show that there are only a small number of differences between these parasite species, with host genetics playing only a minor role in influencing the parasites' response to their intracellular habitat. These small changes may be significant, however, in determining the clinical outcome of infection.
keywords: genetics and genomics/comparative genomics
  genetics and genomics/gene expression
  infectious diseases/neglected tropical diseases
  immunology/immune response
  infectious diseases/skin infections
  infectious diseases/protozoal infections
title: Comparative Expression Profiling of Leishmania: Modulation in Gene Expression between Species and in Different Host Genetic Backgrounds