Tracking of quiescence in Leishmania by quantifying the expression of GFP in the ribosomal DNA locus

Under stressful conditions some microorganisms adopt a quiescent stage characterized by a reversible non or slow proliferative condition that allows their survival. This adaptation was only recently discovered in Leishmania. We developed an in vitro model and a biosensor to track quiescence at population and single cell levels. The biosensor is a GFP reporter gene integrated within the 18S rDNA locus, which allows monitoring the expression of 18S rRNA (rGFP expression). We showed that rGFP expression decreased significantly and rapidly during the transition from extracellular promastigotes to intracellular amastigotes and that it was coupled in vitro with a decrease in replication as measured by BrdU incorporation. rGFP expression was useful to track the reversibility of quiescence in live cells and showed for the first time the heterogeneity of physiological stages among the population of amastigotes in which shallow and deep quiescent stages may coexist. We also validated the use of rGFP expression as a biosensor in animal models of latent infection. Our models and biosensor should allow further characterization of quiescence at metabolic and molecular level.

prokaryotes as Mycobacterium to eukaryotes as Saccharomyces [14][15][16] . Quiescence complexity was further addressed by a third study that explored the homogeneity of quiescence among L. major amastigotes in vivo 17 . The authors found that in persistent infections of the self-healing C57BL/6J mice model, there are two populations among amastigotes: one in a slow proliferative stage (generation time of 60 hours, semi quiescent) and another with no evidence of proliferation (quiescent). These results suggested that amastigotes can thus adopt different physiological stages during the course of the infection.
Molecular mechanism triggering quiescence in Leishmania are still unknown and specific molecular markers are not yet available. As the amastigotes may present two physiologically different stages, one quiescent but another proliferative, there is a need for a biosensor that would allow to distinguish both sub-populations for further molecular characterization. Considering that transcription of rDNA is (i) regulated by a well-defined promotor in Leishmania -in contrast to other genes-and (ii) down-regulated in quiescent stages, we hypothesized that the expression of a reporter gene integrated in the rDNA locus could serve as a biosensor of quiescence. We used here GFP as a reporter gene and demonstrated in L. braziliensis and L. mexicana reference strains that it can be used as a biosensor of the proliferative stage of parasites and as a tool to distinguish quiescent stages among amastigotes at a single cell level. Using this biosensor, we have characterized in vitro the heterogeneity of physiologically downregulated stages among Leishmania amastigotes and the kinetic of the entrance and exit of these stages. The results of experimental work in vivo and the evaluation of quiescence in a panel of clinical isolates are also shown.

Results
rGfp expression in Leishmania promastigotes and amastigotes. The decreased transcription of the rDNA locus is a trait of quiescence in several microorganisms 15,18,19 . Therefore, we postulated that if the GFP gene was integrated within the rDNA locus, the subsequent measurement of its expression (rGFP expression) would serve as a quiescence biosensor and would allow to discriminate quiescent populations among alive parasites. In L. braziliensis, the analysis at whole population level showed that rGFP expression in alive cells (as indicated by their integrity of membrane with NucRed Dead staining, Fig. S1) was the highest in the highly proliferative Pro Log , with an average value of 1995.0 ± 5.7 (mean RFU ± SD). The rGFP expression decreased progressively and significantly to 1361.0 ± 6.2, 217.3 ± 46.6 and 23.37 ± 2.6 in Pro Sta , Amas Axe and Amas Int respectively (One way ANOVA, P = 0.0002) ( Fig. 1A-C). In L. mexicana the same sequential and significant decrease in rGFP expression was observed between the same life stages. (Fig. S2).
We further analyzed the flow cytometry results at single cell level. We observed that for both species, Ama Axe constitute a more heterogeneous population than promastigotes, as evidenced by the thickness of the rGFP expression histogram (Fig. 1B). The analysis of the variances of each population in L. braziliensis showed that Ama Axe have statistically significant higher variability (coefficient of variation (CV) = 90,2) in comparison to Pro Log (CV = 36,0) or Pro Sta (CV = 36,9) (F test < 0.0001). The acquisition of new images with settings to optimize the detection of GFP confirmed the heterogeneity of downregulated stages in Ama Int even among the ones residing in the same macrophages (Fig. 1D). We then calculated in each life stage the proportion of cells that have quantifiable levels of rGFP expression. Using a very stringent gate which allows only 0.04% of false positive in the negative control (non rGFP Leishmania strain), we found that in L. braziliensis the percentage of cells with quantifiable rGFP expression was 99.6 ± 0.1, 99.4 ± 0.1, 52.3 ± 12.8 and 0.56 ± 0.2 (mean ± SD) in Pro Log , Pro Sta , Ama Axe and Ama Int respectively ( Fig. 2A). Very similar results were shown in L. mexicana: in Ama Int of that species, 4.0 ± 0.7% (mean ± SD) reached quantifiable but low rGFP expression (Fig. 2B). This population can be observed in the Fig. S2B as a small peak at the right side of the corresponding graph (arrow). Overall, the results show that the downregulation of rGFP expression is an inherent trait of amastigotes and that it develops rapidly as it is already seen in Ama Axe .
Evaluation of rGFP expression and its relationship with cell proliferation. The deeply downregulated levels of rGFP expression in amastigotes indicated that amastigotes are in a quiescent stage which is characterized by a low or non-proliferative condition. In order to establish the relationship between rGFP expression and the proliferative condition of parasites we evaluated simultaneously both parameter by immunofluorescence. We quantified (i) the proportion of proliferative cells with the use of the BrdU assay which is the standard molecular method to measure DNA replication, together with (ii) the immunofluorescence of rGFP. Images were acquired with a confocal microscope and the quantification of both parameters was performed by flow cytometry (Fig. 3A). Fast incorporation of BrdU and high rGFP expression (~80000 RFU) was observed in Pro Log under conditions where they were replicating with a doubling time of ~8 hrs. After 4, 8 and 24 hrs of incubation with BrdU, there was 47, 80 and 90% of BrdU + Pro Log cells respectively. As expected, with the decreased rGFP in Pro Sta (~30000 RFU) most of the cells were BrdU − with only a small proportion of the population (~10%.) being BrdU + for all the times of exposition (Fig. 3B). The homogeneity of the Pro log population was evinced by the pattern of the scatter plots for the codetection of GFP and incorporation of BrdU by flow cytometry. After 8 hrs of incubation with BrdU most of the population was in the quadrant UR which contains BrdU + and GFP + cells. Negligible fraction of cells were found in the other three quadrants (Fig. 3C, scatter plots P_a to P_c). In both promastigote stages of L. braziliensis the mean of rGFP expression on the population BrdU + was significantly higher in comparison to the one of BrdU − (Fig. 3D).
The population of Ama Axe did not show cellular duplication over 24 hrs, which is in agreement with their even more decreased levels of rGFP expression (~16000 RFU) and slow although progressive increase in their BrdU + population, evincing their diminished rate of proliferation. After 4, 8 and 24 hrs of incubation with BrdU, there was 6, 11 and 29% of BrdU + cells respectively (Fig. 3B). Despite of the homogeneity in the scatter plots for the detection of GFP and BrdU in Pro log an heterogenous pattern was shown in Ama axe in the evaluated times of exposure to BrdU. The bulk of the BrdU + population had low but quantifiable levels of GFP (quadrant UR), while a negligible fraction of BrdU + cells were GFP − (quadrant LR); the GFP − population was also BrdU − and they were found in the quadrant LL of each scatter plot (please see Fig. 3C, scatter plots P_i to P_l). Taking into consideration the incorporation of BrdU after 4 hrs of incubation, we calculated that the replication time in L. braziliensis Pro Log is ~7.6 hrs while it would be ~62.4 hrs in Ama Axe if we consider them as a homogeneous population. However, considering the results above, the occurrence of two physiologically different sub populations in Ama Axe is likely: one with delayed incorporation of BrdU and low but positive rGFP expression and a second one with negligible incorporation of BrdU and unquantifiable rGFP expression (Fig. 3C). In this case the adjustment of the generation time for the rGFP + population gives a replication time of 48.8 hrs. In the rGFP − cells, generation time (if there is any replication) cannot be estimated. Similar results were obtained in L. mexicana where the decreased level of rGFP expression was also associated with a slow rate of proliferation as shown by the diminished percentage of cells incorporating BrdU (Fig. S3).
Evaluation of reversibility of the quiescent stage. We evaluated the kinetics and reversibility of quiescence by the quantification of rGFP expression. The arrested or slowly growing population should be able to revert to its faster rate of growth once parasites are in an environment without nutritional limitations or stressful conditions. With this purpose we analyzed the rGFP expression in L. braziliensis and L. mexicana daily during the process of differentiation from promastigotes to Ama Axe (in response to the acidic pH and high temperatures). In L. braziliensis the levels of rGFP expression during the process of amastigogenesis in comparison to Pro Log decreased significantly to 64 ± 0.2% (mean ± SD), 43.1 ± 0.8%, 28.2 ± 1.4% and 19.8 ± 0.1% on day1, day 2, day 3 and day 4 respectively (One way ANOVA, P < 0.0001) (Fig. 4A). Similar results were found during the amastigogenesis of L. mexicana: rGFP expression decreased to 88.5 ± 2.2 (mean ± SD), 38.9 ± 1.0, 18.7 ± 0.6 and 13.5 ± 0.9% on day 1, day 2, day 3 and day 4 respectively (One way ANOVA P = 0.003) (Fig. 4A). We also evaluated the reversibility of rGFP expression when the parasites differentiate back from Ama Axe to promastigote. The rGFP expression on the whole population of L. braziliensis increased progressively from day 1 (32.1%) to day 3 (86.3%) and growth was resumed after day 1. More heterogeneity was observed in L. mexicana: on day 1 we saw that although most of the population initiated the increase of their levels of rGFP expression, a fraction of them maintained a low rGFP expression as evidenced by the tail in the left of the histogram (Fig. 4B).

Validation of rGFP expression as biosensor of quiescence in animal models with latent infections.
We further wanted to validate the use of rGFP expression as a marker of quiescence in vivo. We monitored amastigotes up to 4 months post infection in animal models where there is not progress to disease but instead latent infections (data not shown). As we expected that recovered quiescent amastigotes would show deeply downregulated levels of rGFP expression, we used CellTracker Deep Red pre-stained promastigotes of rGFP lines to infect the animals (Fig. S4). This staining allowed us tracking the quiescent amastigotes, even if the rGFP was downregulated below the detection limit. In addition, the presence of this staining in the cells also confirmed the quiescent condition of amastigotes. Indeed, proliferative cells are expected to dilute progressively the CellTracker Deep Red until it becomes non-detectable: this was verified in proliferative promastigotes, in which the dye became non detectable on day 6 ( Fig. 5A,B). In contrast, in both hamsters infected with L. braziliensis and mice infected with L. mexicana, intracellular amastigotes positives for CellTracker Deep Red were visualized 15, 60 and 120 days post infection. This confirms a very slow rate of proliferation. Also, the retention of CellTracker Deep Red correlated with the non-detectable (most of the events) or negligible rGFP expression in each visualized amastigote; amastigotes with high rGFP expression and positive for the CellTracker Deep Red were never detected (Fig. 5C). In order to confirm that the amastigotes with negligible rGFP expression and positive for the cell tracker do represent viable cells, an aliquot of the harvested cells was placed in Tobie's medium to allow the differentiation of amastigotes back to promastigotes and their further propagation. For both species, promastigotes were recovered in 100% of the animals 60 days post infection. At 120 days post infection, promastigotes were recovered in 75% and 100% of the animals in L. braziliensis and L. mexicana respectively (Fig. 5D). In order to confirm that the negligible expression of rGFP of in vivo amastigotes is due to their quiescent stage and not due to the loss of GFP within the rDNA locus, we measured the rGFP expression in the isolates recovered 60 days and 120 days post infection. All promastigotes differentiated from the mice-derived amastigotes recovered a pattern of rGFP expression comparable to the original parasites before the infection. For instance, the parasites isolated from 4 mice, 60 days post infection, generated populations with 99.2 ± 0.4% (mean ± SD) of GFP positive cells and with 2077 ± 341 RFU (Fig. 5E). Similar results were obtained with L. mexicana.

Quiescence of L. braziliensis amastigotes in clinical isolates.
In order to extend the relevance of our study to (recent) clinical isolates, we derived rGFP cloned lines from each of them and evaluated their capability to adopt a quiescent stage by measuring the rGFP expression (Fig. 6). We found that as observed in L. braziliensis M2904 rGFP Cl2, Ama Axe of the 8 lines showed downregulation of rGFP expression in Ama Axe , with levels ranging from ~26% to only ~6% of the levels of Pro Log (Fig. 6B). These results show all the clinical isolates can adopt a quiescent stage, and that there is some variability in the extend of downregulation among these.

Discussion
Many pathogens can remain in their host and cause a subclinical infection after chemotherapy or self-healing. These latent infections can extend during the lifespan of the host and constitute a reservoir for reactivation of the disease in an unpredictable time. This is the case in leishmaniasis, where relapse is common in the absence of drug resistance 20,21 , and complications like muco-cutaneous leishmaniasis or post kala-azar dermal leishmaniasis occur months to years after healing of primary cutaneous or visceral lesions, respectively 22,23 . In other pathogens, in vivo models of subclinical infections have shown evidence that hosts harbor the pathogen in a quiescent stage with no or slow replication 8,9,11,24 . However, the physiology of these stages is barely characterized in vivo because of their scarce number and the difficulty to detect them. We aimed here to develop and validate a simple bio-marker of quiescence and an in vitro model for further studies of this phenomenon in Leishmania. In this study, we inserted the GFP gene into the rDNA locus to evaluate if the expression of this reporter gene driven by the ribosomal promoter may serve as a biosensor of the metabolic activity of the parasite. We showed that (i) rGFP expression decreased significantly and rapidly during the in vitro transition from extracellular promastigotes to intracellular amastigotes, (ii) the decrease in rGFP expression was coupled in vitro with a decrease in replication as measured by BrdU incorporation, (iii) modulation of rGFP expression both in vitro and in vivo was reversible, the signal of the reporter gene increasing again during the transition from amastigote to promastigote, (iv) the use of rGFP expression allowed to identify heterogeneity among amastigotes and (v) our rGFP bio-sensor was applicable in vivo.
A common feature of quiescent populations in a broad range of microbes is the reduction in the abundance of ribosomal RNAs and corresponding proteins 15,25,26 . Down-regulation of biosynthetic machinery, including ribosomes, fits with cells that do not need to duplicate their DNA and synthetize all the proteins and other macromolecules to generate daughter cells. This was clearly demonstrated in Mycobacterium bovis where the fast proliferative stages had ~3880 ribosomes per cell while the slow proliferative stages had 687 ribosomes 27 . In L. infantum Ama Axe showed a decreased rate of synthesis of proteins and amount of polyribosomes 28 . Similarly, our previous study in L. braziliensis showed that Ama Axe downregulate their levels of rRNA to 15% of the levels found in Pro Log 13 . In present study, we showed that Ama Axe of both L. mexicana and L. braziliensis had an overall downregulated rGFP expression equivalent to 10% of the levels found in Pro Log . This reduction becomes even more stringent in Ama Int as the overall expression of the rGFP dropped to less than 1% in L. mexicana and 3% in L. braziliensis in comparison to the levels found in Pro Log . This higher drop can be an adaptative mechanism that allows the amastigotes to survive in the phagolysosome where they will face the limitation of micronutrients and the macrophage microbicidal processes [29][30][31] . We further extended the documentation of the quiescence phenomenon to clinical isolates. Considerable variability in the downregulation of rGFP expression suggests that clinical isolates could have different propensity to adopt the quiescent stage.
The inverse relationship between proliferation and biogenesis of ribosomes offers a set of potential molecular markers for identification of quiescent cells 32 . In the present study, the results of the BrdU assay for proliferation confirm that rGFP expression can be used as a marker for the identification of quiescent populations of amastigotes as in Ama Axe only the population with measurable but low rGFP expression showed some incorporation of BrdU while the population with no rGFP expression did not show any detectable incorporation of BrdU. The deep down regulation of biosynthesis (here measured by rGFP expression) and the negligible incorporation of BrdU in Ama Int are in agreement with the constant infection index in peritoneal macrophages infected with L. braziliensis and L. mexicana over 72 hrs post infection: this was an early indicator of the lower rate of proliferation of Ama Int in comparison to the ~7 hrs required by promastigotes [33][34][35] . The use of rGFP as a marker of quiescence will facilitate further molecular characterization of quiescent subpopulations.
In pathogens reversibility of the quiescent stage to a proliferative condition is imperative to enable transmission and survival of the population 36 . This was demonstrated here, as both axenic and rodent ex vivo amastigotes were also capable to rapidly resume their active biosynthetic level and proliferative condition when they were placed in conditions optimal to differentiate back to the proliferative promastigote stage. Future studies, monitoring rGFP expression of quiescent cells in real time with technologies as time lapse microscopy and microfluidics would allow to establish the cut off levels of rGFP expression indicative of parasites resuming the proliferative condition. These approaches may facilitate the selection/or sorting of these cells for further characterization of the molecular mechanism needed for both the entry and exit of quiescence.
Leishmania amastigotes may show a range of physiologic heterogeneity. It has been shown that 4 months after the infection with L. major, the healed foot pad of C57BL/6J mice harbors two populations of amastigotes: one that replicates with a time of generation of 60 hrs and a second one that do not show evidence of replication 17 . Monitoring rGFP expression in L. mexicana and L. braziliensis amastigotes supported this concept of heterogeneity among amastigotes through the co-existence of amastigotes with varying downregulated physiological stages. Similarly, in L. amazonensis it has been shown that clones derived of parasites freshly taken from lesion have three phenotypes: fast, slow and intermediate progression to pathology in BALB/c mice 37 . These phenotypes were not strictly genetically based as a fast clone generated also sub clones with the three phenotypes 37 . Our results in Leishmania cloned populations support the growing evidence that quiescence in pathogens is a dynamic phenomenon where subpopulations of shallow to deep quiescent stages may form the whole population 9,38,39 . It www.nature.com/scientificreports www.nature.com/scientificreports/ is postulated that this heterogeneity allows the survival of the population rather than individual cells in different environmental conditions 36,37 .
In vivo Leishmania infections have different outcome depending on the strain, model and the number of parasites inoculated (38). In the present study we used Leishmania strains which upon inoculation generated asymptomatic infection during the 4 months of the experimental work. Although no lesions were developed, we found amastigotes in the tissues. The quiescent condition of parasites in these latent infections was confirmed by two indicators: (i) the negligible expression of rGFP expression and (ii) the retention of the CellTracker Deep Red up to 120 days post infection. Our results show that the inoculation of a moderate dose of stationary parasites in hamsters and mice for L. braziliensis and L. mexicana respectively (and the strains here evaluated) represents a good model for the further molecular characterization of quiescent population in latent infections. Although, the integration of fluorescent or luminescent reporters within the rDNA locus is useful to monitor quiescence, its application to monitor infections in vivo could be limited as the downregulation in the transcription of the rDNA locus may bring the expression of the reporter under the limit of detection of the in vivo imaging systems. Therefore, the lack of signal in animal tissues does not necessarily indicate the absence of infection. In order to monitor infections a second reporter gene should be introduced in another locus with constitutive expression in amastigote.
By using rGFP as biosensor of quiescence, we showed that L. braziliensis and L. mexicana strains produce a heterogeneous population of quiescent amastigotes that are characterized by a decreased rate of proliferation and a lower capacity of biosynthesis, both reversible. The use of transgenic rGFP parasites and single cell approaches will offer valuable tools for the further molecular in vitro studies of quiescent subpopulation among amastigotes. The relevance of quiescence in subclinical infections and its potential role in non-heritable drug tolerance in Leishmania urge additional exploration as it has been shown that it plays an important role in other pathogens [40][41][42] .

Materials and Methods
Ethics statement. Mice  In the next paragraphs the abbreviation Pro Log and Pro Sta will be used for logarithmic and stationary promastigotes, respectively.
In order to obtain axenic amastigotes, promastigotes (1.0 × 10 6 parasites/mL) were maintained in complete M199. On the fourth day, 1 mL of promastigotes was centrifuged (1500 g) and the pellet was re-suspended in 5 mL of complete M199 at pH 5.5 and incubated at 34 °C. After four days, the resulting axenic amastigotes (Ama Axe ) were harvested or used to obtain in vitro intracellular amastigotes (Ama Int ) by the infections of BALB/c peritoneal macrophages at a ratio of 4 Ama axe to 1 macrophage as previously described elsewhere 43 . Three biological replicates per each stage and condition were evaluated in all the tests in vitro.

Development of Gfp
Leishmania lines. GFP was integrated within the 18S ribosomal DNA locus with the use of the pLEXSY-neo2 system (Jena Bioscience) as previously reported elsewhere 44 . After the generation of the transgenic rGFP parasites, they were cloned by the 'micro-drop' method. Briefly, an appropriate dilution of a 2 days old promastigote culture was prepared from which micro-drops of 1 µl were placed in the middle of the well (96 wells plate). In each well the presence of a drop with a single promastigote was verified microscopically and upon confirmation 150 µl of complete M199 medium was added. After incubation at 26 °C for one week, growth was microscopically examined and established clones were transferred to a biphasic Tobie's medium for further routine culture. The number of integrations of GFP within the rDNA loci was confirmed by whole genome sequencing of a PCR free library with a HiSeq. 4000 platform (BGI). The sequences were assembled having as reference the genome of the strain M2904 as previously described elsewhere 45 . The sequence analysis showed the presence of 1 and 2 copies of GFP per genome in the lines of L. braziliensis and L. mexicana respectively.
Monitoring rGFP expression in living parasites by flow cytometry. Promastigotes and amastigotes of rGFP lines were stained for the evaluation of their cell viability by using the NucRed dead 647 (Thermo fisher scientific): labeling marks dead cells. The samples were analyzed with a calibrated flow cytometer BD FACS Verse TM in the medium flow rate mode. A wild type (non rGFP) line, not stained with NucRed was included in each experiment as negative control. In order to compare the fluorescence among samples the relative fluorescence units (RFU) were acquired with the same settings during all the experiments. The FCS files were analyzed with the FCS 5 express Plus research edition. Briefly, single cells were selected by pulse geometry gate; a first gate was selected with SSC-W and SSC-H, a second gate was created with FCS-H and FCS-A. Among single cells, a third gate was created by plotting the rGFP expression and the NucRed staining; the levels of rGFP expression was evaluated in cells negatives for the NucRed dead staining. (2019) 9:18951 | https://doi.org/10.1038/s41598-019-55486-z www.nature.com/scientificreports www.nature.com/scientificreports/ Double immunodetection of BrdU and rGFP for confocal microscopy and flow cytometry. To test BrdU labeling in promastigotes and amastigotes the cells were cultured in complete M199 media containing 0.05 mM BrdU (Sigma) for 4, 8 or 24 hrs, after which they were washed and fixed in 4% (w/vol) paraformaldehyde in PBS for 20 min and stained as described below. For all BrdU immunodetection experiments, paraformaldehyde-fixed samples were permeabilized with block perm buffer (PBS with 5% normal goat sera and 0.1% Triton-X-100) for 15 min, washed in distilled water, and then immersed in 2 M HCl for 40 min. After extensive washing with PBS, the samples were incubated in a block perm buffer for 15 min. Subsequently, BrdU and rGFP were labelled simultaneously overnight at 4 °C by using 1/100 dilution of the rat monoclonal anti BrdU (Abcam) and 1/100 dilution of the rabbit monoclonal anti GFP (Abcam) as primary antibodies. The next day the samples were washed two times with PBS, after which they were treated with block perm buffer for 5 min. Samples were subsequently incubated for 1.5 hrs with Alexa 555 goat polyclonal anti rat IgG (1/1000) and Alexa 488 goat polyclonal anti rabbit IgG (1/1000) (Thermo fisher scientific). After two washes with PBS the samples were mounted with Prolong Gold antifade (Thermo fisher scientific) for confocal microscopy or diluted in PBS for the analysis with the flow cytometer. Microscopy was performed on a Zeiss LSM 700 laser scanning microscope. Image analysis was performed using ImageJ software. The samples for flow cytometry were analyzed with the BD FACS Verse TM in the medium flow rate mode. The time of generation was estimated by extrapolating the time needed to have 100% of the population BrdU + from the percentage of BrdU + cells after 4 hours of incubation.
Animal models for latent infection with L. braziliensis and L. mexicana GFP lines and screening of living amastigotes ex vivo. The infection with Leishmania in animal models causes several outcomes depending on the line of the rodents, the infecting species, the stage of the parasite and the amount of the inoculum 46 . Here we use a hamster and mice model where there is not development of lesions but instead latent infection. Briefly, twelve mice (BALB/c) or hamsters (LVG Golden Syrian Hamster) provided by Charles River were anesthetized by intraperitoneal injection of Ketamine (1 mg/10 g body weight) and Xylazine (0.08 mg/10 g body weight). Anesthetized animals were infected in both hind footpads by sub-cutaneous injection of 5 × 10 5 stationary phase parasites pre-stained with CellTracker Deep Red in 20 µl of saline solution (0.3 mL syringe and a 30 G needle gauge). Four mice and hamsters were sacrificed at 2 weeks, 8 weeks and 16 weeks post infection.
The feet of animals were collected and minced with scissors. The pieces were placed in 1.0 ml of complete M199 containing 200 µg/ml of Liberase TL (Roche) and incubated at 34 °C for 30 min. Subsequently, the samples were filtered with a 70 μm cell strainer and washed with 10 ml of PBS. After centrifugation (1500 g) the pellets were resuspended in 1 mL of complete M199 and aliquot was placed in a 96 well Glass bottom microplates. The screening of living amastigotes was performed with the confocal microscope with settings for the codetection of GFP and the Celltracker Deep Red.