Heat-shock protein 70 gene sequencing for Leishmania species typing in European tropical infectious disease clinics.

We describe Leishmania species determination on clinical samples on the basis of partial sequencing of the heat-shock protein 70 gene (hsp70), without the need for parasite isolation. The method is especially suited for use in non-endemic infectious disease clinics dealing with relatively few cases on an annual basis, for which no fast high throughput diagnostic tests are needed. We show that the results obtained from this gene are in nearly perfect agreement with those from multilocus enzyme electrophoresis, which is still considered by many clinicians and the World Health Organization (WHO) as the gold standard in Leishmania species typing. Currently, 203 sequences are available that cover the entire hsp70 gene region analysed here, originating from a total of 41 leishmaniasis endemic countries, and representing 15 species and sub-species causing human disease. We also provide a detailed laboratory protocol that includes a step-by-step procedure of the typing methodology, to facilitate implementation in diagnostic laboratories.

Universal PCR assays for the differential detection of all Old World Leishmania species.

For the epidemiological monitoring and clinical case management of leishmaniasis, determination of the causative Leishmania species gains importance. Current assays for the Old World often suffer from drawbacks in terms of validation on a geographically representative sample set and the ability to recognize all species complexes. We want to contribute to standardized species typing for Old World leishmaniasis. We determined the ribosomal DNA internal transcribed spacer 1 sequence of 24 strains or isolates, and validated four species-specific polymerase chain reactions (PCRs) amplifying this target. They discriminate L. aethiopica, L. tropica, L. major, and the L. donovani complex, use the same cycling conditions, and include an internal amplification control. Our PCRs amplify 0.1 pg of Leishmania DNA, while being 100% specific for species identification on an extensive panel of geographically representative strains and isolates. Similar results were obtained in an endemic reference laboratory in Kenya. Species could also be identified in clinical specimens. The presented PCRs require only agarose gel detection, and have several other advantages over many existing assays. We outline potential problems, suggest concrete solutions for transferring the technique to other settings, and deliver the proof-of-principle for analyzing clinical samples.

Increased metacyclogenesis of antimony-resistant Leishmania donovani clinical lines.

Mathematical models predict that the future of epidemics of drug-resistant pathogens depends in part on the competitive fitness of drug-resistant strains. Considering metacyclogenesis (differentiation process essential for infectivity) as a major contributor to the fitness of Leishmania donovani, we tested its relationship with pentavalent antimony (SbV) resistance in clinical lines. Different methods for the assessment of metacyclogenesis were cross-validated: gene expression profiling (META1 and SHERP), morphometry (microscopy and FACS), in vitro infectivity to macrophages and resistance to complement lysis. This was done on a model constituted by 2 pairs of reference strains cloned from a SbV-resistant and -sensitive isolate. We selected the most adequate parameter and extended the analysis of metacyclogenesis diversity to a sample of 20 clinical lines with different in vitro susceptibility to the drug. The capacity of metacyclogenesis, as measured by the complement lysis test, was shown to be significantly higher in SbV-resistant clinical lines of L. donovani than in SbV-sensitive lines. Together with other lines of evidence, it is concluded that L. donovani constitutes a unique example and model of drug-resistant pathogens with traits of increased fitness. These findings raise a fundamental question about the potential risks of selecting more virulent pathogens through massive chemotherapeutic interventions.

Comparison of gene expression patterns among Leishmania braziliensis clinical isolates showing a different in vitro susceptibility to pentavalent antimony.

INTRODUCTION: Evaluation of Leishmania drug susceptibility depends on in vitro Sb(V) susceptibility assays, which are labour-intensive and may give a biased view of the true parasite resistance. Molecular markers are urgently needed to improve and simplify the monitoring of Sb(V)-resistance. We analysed here the gene expression profile of 21 L. braziliensis clinical isolates in vitro defined as Sb(V)-resistant and -sensitive, in order to identify potential resistance markers. METHODS: The differential expression of 13 genes involved in Sb(V) metabolism, oxidative stress or housekeeping functions was analysed during in vitro promastigote growth. RESULTS: Expression profiles were up-regulated for 5 genes only, each time affecting a different set of isolates (mosaic picture of gene expression). Two genes, ODC (ornithine decarboxylase) and TRYR (trypanothione reductase), showed a significantly higher expression rate in the group of Sb(V)-resistant compared to the group of Sb(V)-sensitive parasites (P<0.01). However, analysis of individual isolates showed both markers to explain only partially the drug resistance. DISCUSSION: Our results might be explained by (i) the occurrence of a pleiotropic molecular mechanism leading to the in vitro Sb(V) resistance and/or (ii) the existence of different epi-phenotypes not revealed by the in vitro Sb(V) susceptibility assays, but interfering with the gene expression patterns.

Development and evaluation of different PCR-based typing methods for discrimination of Leishmania donovani isolates from Nepal.

INTRODUCTION: Leishmania donovani, the causative agent of visceral leishmaniasis in the Indian subcontinent, has been reported to be genetically homogeneous. In order to support ongoing initiatives to eliminate the disease, highly discriminative tools are required for documenting the parasite population and dynamics. METHODS: Thirty-four clinical isolates of L. donovani from Nepal were analysed on the basis of size and restriction endonuclease polymorphisms of PCR amplicons from kinetoplast minicircle DNA, 5 nuclear microsatellites, and nuclear loci encoding glycoprotein 63, cysteine proteinase B, and hydrophilic acylated surface protein B. We present and validate a procedure allowing standardized analysis of kDNA fingerprint patterns. RESULTS: Our results show that parasites are best discriminated on the basis of kinetoplast minicircle DNA (14 genotypes) and 1 microsatellite defining 7 genotypes, while the remaining markers discriminated 2 groups or were monomorphic. Combination of all nuclear markers revealed 8 genotypes, while extension with kDNA data yielded 18 genotypes. CONCLUSION: We present tools that allow discrimination of closely related L. donovani strains circulating in the Terai region of Nepal. These can be used to study the micro-epidemiology of parasite populations, determine the geographical origin of infections, distinguish relapses from re-infection, and monitor the spread of particular variants.

Heat-shock protein 70 PCR-RFLP: a universal simple tool for Leishmania species discrimination in the New and Old World.

INTRODUCTION: Species typing in leishmaniasis gains importance in diagnostics, epidemiology, and clinical studies. A restriction fragment length polymorphism (RFLP) assay of PCR amplicons from a partial heat-shock protein 70 gene (hsp70) had been described for the New World, allowing identification of some species. METHODS: Based on an initial in silico analysis of 51 hsp70 sequences, most of which we recently determined in the frame of a phylogenetic study, species-specific restriction sites were identified. These were tested by PCR-RFLP on 139 strains from 14 species, thereby documenting both inter- and intra-species variability. RESULTS: Our assay could identify Leishmania infantum, L. donovani, L. tropica, L. aethiopica, L. major, L. lainsoni, L. naiffi, L. braziliensis, L. peruviana, L. guyanensis, and L. panamensis by applying 2 subsequent digests. L. mexicana, L. amazonensis, and L. garnhami did not generate species-specific restriction fragment patterns. CONCLUSION: Currently no assay is available for global Leishmania species discrimination. We present a universal PCR-RFLP method allowing identification of most medically relevant Old and New World Leishmania species on the basis of a single PCR, obviating the need to perform separate PCRs. The technique is simple to perform and can be implemented in all settings where PCR is available.

Evolutionary history of Trypanosoma cruzi according to antigen genes.

Trypanosoma cruzi, the agent of Chagas disease is associated with a very high clinical and epidemiological pleomorphism. This might be better understood through studies on the evolutionary history of the parasite. We explored here the value of antigen genes for the understanding of the evolution within T. cruzi. We selected 11 genes and 12 loci associated with different functions and considered to be involved in host-parasite interaction (cell adhesion, infection, molecular mimicry). The polymorphism of the respective genes in a sample representative of the diversity of T. cruzi was screened by PCR-RFLP and evolutionary relationships were inferred by phenetic analysis. Our results support the classification of T. cruzi in 2 major lineages and 6 discrete typing units (DTUs). The topology of the PCR-RFLP tree was the one that better fitted with the epidemiological features of the different DTUs: (i) lineage I, being encountered in sylvatic as well as domestic transmission cycles, (ii) IIa/c being associated with a sylvatic transmission cycle and (iii) IIb/d/e being associated with a domestic transmission cycle. Our study also supported the hypothesis that the evolutionary history of T. cruzi has been shaped by a series of hybridization events in the framework of a predominant clonal evolution pattern.

Evaluation of an in vitro and in vivo model for experimental infection with Leishmania (Viannia) braziliensis and L. (V.) peruviana.

Leishmania (Viannia) braziliensis and L. (V.) peruviana are two parasite species characterized by a very different pathogenicity in humans despite a high genetic similarity. We hypothesized previously that L. (V.) peruviana would descend from L. (V.) braziliensis and would have acquired its 'peruviana' character during the southward colonization and adaptation of the transmission cycle in the Peruvian Andes. In order to have a first appreciation of the differences in virulence between both species, we evaluated an in vitro and in vivo model for experimental infection. A procedure was adapted to enrich culture forms in infective stages and the purified metacyclics were used to infect macrophage cell lines and golden hamsters. The models were tested with 2 representative strains of L. (V.) braziliensis from cutaneous and mucosal origin respectively and 2 representative strains of L. (V.) peruviana from Northern and Southern Peru respectively. Our models were reproducible and sensitive enough to detect phenotypic differences among strains. We showed in vitro as well as in vivo that the L. (V.) braziliensis was more infective than L. (V.) peruviana. Furthermore, we found that in vitro infectivity patterns of the 4 strains analysed, were in agreement with the geographical structuring of parasite populations demonstrated in our previous studies. Further work is needed to confirm our results with more strains of different geographical origin and their specific clinical outcome. However, our data open new perspectives for understanding the process of speciation in Leishmania and its implications in terms of pathogenicity.

American tegumentary leishmaniasis: antigen-gene polymorphism, taxonomy and clinical pleomorphism.

Multi-locus enzyme electrophoresis is the current gold standard for the genetic characterisation of Leishmania. However, this method is time-consuming and, more importantly, cannot be directly applied to parasites present in host tissue. PCR-based methods represent an ideal alternative but, to date, a multi-locus analysis has not been applied to the same sample. This has now been achieved with a sample of 55 neotropical isolates (Leishmania (Viannia) braziliensis, L. (V.) peruviana, L. (V.) guyanensis, L. (V.) lainsoni and L. (L.) amazonensis), using five different genes as targets, four of which encoded major Leishmania antigens (gp63, Hsp70, H2B and Cpb). Our multi-locus approach strongly supports the current taxonomy and demonstrates a highly robust method of distinguishing different strains. Within L. (V.) braziliensis, we did not encounter so far specific genetic differences between parasites isolated from cutaneous and mucosal lesions. Interestingly, results provided by each of the different antigen-genes in the species considered, were different, suggesting different selective pressures. Our work emphasises the need for a multi-disciplinary approach to study the clinical pleomorphism of leishmaniasis.

Relation between variation in copy number of ribosomal RNA encoding genes and size of harbouring chromosomes in Leishmania of subgenus Viannia.

Chromosomal size polymorphism in Leishmania of subgenus Viannia has been correlated with eco-geography. The sizes of chromosomes bearing rDNA genes were determined in 69 isolates. A considerable size-variation was observed, ranging from 1100 to 1500 kb. Chromosomes of L.(V.). braziliensis, L.(V.)guyanensis and L.(V.) peruviana from northern Peru were significantly larger (200 kb) than those of L.(V.) peruviana from southern Peru. In addition, 31 out of 69 isolates presented each two different-sized homologues of the rDNA chromosome. Long range restriction mapping of three different-sized rDNA chromosomes from L.(V.)braziliensis M2903 and L.(V.)peruviana HB31 (north) and LC106 (south) each revealed three fragments delimited by PmeI restriction sites: two constant in size (the centre and one extremity of the chromosome) and one variable (the other extremity, containing a single cluster of rDNA genes). Further analysis of the M2903 rDNA chromosome allowed the localization of its 140 kb rDNA cluster at 85 kb from the telomeric end. Two arguments indicated that size-variation of the rDNA chromosome is partially due to amplification/deletion of the clustered rDNA genes: (i) size-variation of the cluster-containing fragment was proportional to the size-variation of the whole chromosome, and (ii) hybridization signal intensity of the rDNA chromosome with a small subunit rDNA probe strongly correlated with chromosomal size. Nevertheless, DNA sequences present between the rDNA cluster and the telomere might also play a role in chromosomal size polymorphism. In addition, our data suggest that rDNA gene copy number (20-40 copies cell(-1) under a diploid hypothesis) in subgenus Viannia is lower than reported previously.

Chromosome rearrangement in Leishmania mexicana M379.

Circular extrachromosomal elements were observed in a variety of Leishmania species. We show here that two lines originating from the same isolate have been found to contain a circular DNA molecule of 26.6 kb and a linear chromosome of about 250 kb, respectively, which share a homology of more than 20 kb. The circular DNA molecule and its related region on the linear chromosome were cloned and their restriction maps compared. This investigation reveals information about chromosome rearrangement in L. mexicana M379. Further examination will enable us to understand the nature of chromosome rearrangement such as circularization or linearization.

Size-polymorphism of mini-exon gene-bearing chromosomes among natural populations of Leishmania, subgenus Viannia.

In order to explore genomic plasticity at the level of the mini-exon gene-bearing chromosome in natural populations of Leishmania, the molecular karyotype of 84 Leishmania stocks belonging to subgenus Viannia, originating mostly from Peru and Bolivia, and differing according to eco-geographical and clinical parameters, was resolved and hybridised with a mini-exon probe. The results suggest that size variation of the mini-exon gene-bearing chromosome is frequent and important (up to 245-kb size-difference), and partially involves variation (up to 50%) in copy number of mini-exon genes. There is no significant size-difference between mini-exon-bearing chromosomes of Peruvian and Bolivian populations of cutaneous and mucosal isolates of Leishmania (Viannia) braziliensis, but there is between eco-geographical populations of Leishmania (Viannia) peruviana. Leishmania (V.) peruviana presented a significantly smaller mini-exon-bearing chromosome than the other species of subgenus Viannia. The contrast between the general chromosome size heterogeneity and the homogeneity observed in some Peruvian Andean areas is discussed in terms of selective pressure.

Molecular epidemiology and diagnosis of Leishmania: what have we learnt from genome structure, dynamics and function?

This paper reviews our exploration of the dynamics of the Leishmania genome and its contribution to epidemiology and diagnosis. We used as a model Peruvian populations of L. (Viannia) braziliensis and L. (V.) peruviana, 2 species very close phylogenetically, but phenotypically very different in biotope and pathology. We initially focused on karyotype analysis. Our data showed that chromosomes were subject to a fast rate of evolution, and were sensitive indicators of genetic drift. Therefore, molecular karyotyping appeared an adequate tool for monitoring (i) emergence of close species, (ii) ecogeographical differentiation at the intraspecific level, and (iii) strain 'fingerprinting'. Chromosome size variation was mostly due to the number of tandemly repeated genes (rDNA, mini-exon, gp63, and cysteine proteinase genes), and could involve the deletion of unique genes (L. (V.) braziliensis-specific gp63 families). Considering the importance of these genes in parasitism, their rearrangement might have functional implications: adaptation to different environments and pleomorphic pathogenicity. Our knowledge of genome structure and dynamics was used to develop new polymerase chain reaction (PCR) techniques. Amplification of gp63 genes followed by cleavage with restriction enzymes and study of restriction fragment length polymorphism (gp63 PCR-RFLP) allowed the discrimination of all species tested, even directly in biopsies with 95% sensitivity (compared with PCR amplification of kinetoplast deoxyribonucleic acid). At the intra-specific level, RFLP was also observed and corresponded to mutations in major immunogen domains of gp63. These seem to be under strong selection pressure, and the technique should facilitate addressing how the host's immune pressure may modulate parasite population structure. Altogether, gp63 PCR-RFLP represents a significant operational improvement over the other techniques for molecular epidemiology and diagnosis: it combines sensitivity, discriminatory power and prognostic value.

Identification of Lebanese dermotropic putative Leishmania archibaldi isolates by gp63 PCR-RFLP.

Leishmania stocks isolated from cutaneous lesions in Lebanon were characterized by PCR methods. The stocks were typed as putative L. (L.) archibaldi (gp63 PCR-RFLP), belonging to 2 different genotypes (PCR-based schizodeme analysis). This constitutes the first report on the presence of L. (L.) archibaldi in the Middle East.

In vitro promastigote fitness of putative Leishmania (Viannia) braziliensis/Leishmania (Viannia) peruviana hybrids.

In order to initiate studies on the phenotypic properties of hybrids vs. their putative parents, the in vitro growth behaviour of promastigotes was compared for 15 stocks characterised as Leishmania (Viannia) braziliensis, Leishmania (Viannia) peruviana and putative hybrids (isolated from the Eastern Andean valley of Huanuco, Peru). Five sets of three stocks, each set including a L.(V.)braziliensis, a L.(V.)peruviana and a putative hybrid, were constituted randomly and counted daily close to isolation from man (ten to 18 subcultures). Hybrids and L.(V.)peruviana presented similar growth characteristics, and they displayed a growth capacity (growth rate and cell density at stationary phase) significantly lower than the one of L.(V.)braziliensis. Following prolonged in vitro maintenance of one of the sets, the hybrid kept its lower growth capacity. The contrast between the difficulty to grow in vitro these putative hybrids, and their high isolation rate from natural populations is discussed.

A new focus of cutaneous leishmaniasis due to Leishmania amazonensis in a Sub Andean region of Bolivia.

We detected a new outbreak focus with high incidence of cutaneous leishmaniasis in the Sub Andean region of La Paz. This area was never considered previously as an endemic zone of leishmaniasis. Leishmania stocks from human lesions were isolated: three stocks were explored by pulse field gradient electrophoresis, showing evidence for their affiliation to the L. mexicana complex. Eight stocks were submitted to isoenzyme electrophoresis and compared with five reference strains: L. amazonensis, L. braziliensis, L. chagasi, L. mexicana and L. pifanoi. Close genetic proximity was evidenced between newly isolated parasites and the reference stock of L. amazonensis, whereas high divergence was observed between them and either the L. pifanoi, L. mexicana, L. braziliensis and L. chagasi reference strains.

Gene expression profiling of Leishmania (Leishmania) donovani: overcoming technical variation and exploiting biological variation.

Gene expression profiling is increasingly used in the field of infectious diseases for characterization of host, pathogen and the nature of their interaction. The purpose of this study was to develop a robust, standardized method for comparative expression profiling and molecular characterization of Leishmania donovani clinical isolates. The limitations and possibilities associated with expression profiling in intracellular amastigotes and promastigotes were assessed through a series of comparative experiments in which technical and biological parameters were scrutinized. On a technical level, our results show that it is essential to use parasite harvesting procedures that involve minimal disturbance of the parasite's environment in order to 'freeze' gene expression levels instantly; this is particularly a delicate task for intracellular amastigotes and for specific 'sensory' genes. On the biological level, we demonstrate that gene expression levels fluctuate during in vitro development of both intracellular amastigotes and promastigotes. We chose to use expression-curves rather than single, specific, time-point measurements to capture this biological variation. Intracellular amastigote protocols need further refinement, but we describe a first generation tool for high-throughput comparative molecular characterization of patients' isolates, based on the changing expression profiles of promastigotes during in vitro differentiation.

Clonal propagation and the fast generation of karyotype diversity: An in vitro Leishmania model.

In the present work we studied the karyotype stability during long-term in vitro maintenance in 3 cloned strains of Leishmania (Viannia) peruviana, Leishmania (Viannia) braziliensis and a hybrid between both species. Only the L. (V.) peruviana strain showed an unstable karyotype, even after subcloning. Four chromosomes were studied in detail, each of them characterized by homologous chromosomes of different size (heteromorphy). Variations in chromosome patterns during in vitro maintenance were rapid and discrete, involving loss of heteromorphy or appearance of additional chromosome size variants. The resulting pattern was not the same according to experimental conditions (subinoculation rate or incubation temperature), and interestingly, this was associated with differences in growth behaviour of the respective parasites. No change in total ploidy of the cells was observed by flow cytometry. We discuss several mechanisms that might account for this variation of chromosome patterns, but we favour the occurrence of aneuploidy, caused by aberrant chromosome segregation during mitosis. Our results provide insight into the generation of karyotype diversity in natural conditions and highlight the relativity of the clone concept in parasitology.

Complexity of the major surface protease (msp) gene organization in Leishmania (Viannia) braziliensis: evolutionary and functional implications.

The major surface protease (msp or gp63) of Leishmania plays a major role in the host-parasite interaction. We analysed here the structure of the msp gene locus in Leishmania (Viannia) braziliensis and compared it to results obtained in other species. Physical mapping of cosmid contigs revealed a minimum of 37 genes per haploid genome and at least 8 different msp gene families. Within the same organism, these genes showed a nucleotide sequence varying in certain stretches from 3 to 34%, and a mosaic structure. From an evolutionary point of view, major differences were observed between subgenera Viannia and Leishmania, both in terms of msp gene number and sequence. Within subgenus Viannia, phenetic analysis revealed three clusters in which sequence variants of L. (Viannia) braziliensis and L. (Viannia) guyanensis were interspersed. Functional implications of our results were explored from predicted L. (Viannia) braziliensis protein sequences: regions encoding the msp catalytic site showed a conserved sequence, while regions encoding surface domains possibly involved in the host-parasite interaction (macrophage adhesion sites and immunodominant B-cell and T-cell epitopes) were variable. We speculate that this would be an adaptive strategy of the parasite.

Plasticity of gp63 gene organization in Leishmania (Viannia) braziliensis and Leishmania (Viannia) peruviana.

The genomic organization of gp63 genes in 4 and 7 isolates of Leishmania braziliensis and L. peruviana, respectively was studied by RFLP analysis with 3 restriction enzymes (Bgl I, Sal I and Apa I). Our results showed a marked polymorphism among isolates. Some characters were specific to L. braziliensis or to L. peruviana, and others specific to the respective biogeographical populations of L. peruviana. The average minimum copy number of gp63 genes was found to be higher in L. braziliensis (71) than in L. peruviana (46), suggesting that deletion of gp63 genes might be partially involved in the size decrease of the chromosome bearing gp63 genes, observed between those 2 species (from 700 to 610 kb). Our results may suggest the existence of at least 2 arrays of heterologous gp63 repeats, varying in relative copy number between L. braziliensis and L. peruviana, and among isolates of the latter species. Rearrangement of the gp63 genes was observed during long-term in vitro maintenance of a reference strain of L. braziliensis. These observations document the existence of a dynamic gp63 gene organization in Leishmania of the braziliensis complex.