Evolution of the Leishmania braziliensis species complex from amplified fragment length polymorphisms, and clinical implications

Diversity and dissemination of viruses in pathogenic protozoa

Viruses are the most abundant biological entities on Earth and play a significant role in the evolution of many organisms and ecosystems. In pathogenic protozoa, the presence of viruses has been linked to an increased risk of treatment failure and severe clinical outcome. Here, we studied the molecular epidemiology of the zoonotic disease cutaneous leishmaniasis in Peru and Bolivia through a joint evolutionary analysis of Leishmania braziliensis and their dsRNA Leishmania virus 1. We show that parasite populations circulate in tropical rainforests and are associated with single viral lineages that appear in low prevalence. In contrast, groups of hybrid parasites are geographically and ecologically more dispersed and associated with an increased prevalence, diversity and spread of viruses. Our results suggest that parasite gene flow and hybridization increased the frequency of parasite-virus symbioses, a process that may change the epidemiology of leishmaniasis in the region.

Parasite hybridization promotes spreading of endosymbiotic viruses

Viruses are the most abundant biological entities on Earth and play a significant role in the evolution of many organisms and ecosystems. In pathogenic protozoa, the presence of endosymbiotic viruses has been linked to an increased risk of treatment failure and severe clinical outcome. Here, we studied the molecular epidemiology of the zoonotic disease cutaneous leishmaniasis in Peru and Bolivia through a joint evolutionary analysis of Leishmania braziliensis parasites and their endosymbiotic Leishmania RNA virus. We show that parasite populations circulate in isolated pockets of suitable habitat and are associated with single viral lineages that appear in low prevalence. In contrast, groups of hybrid parasites were geographically and ecologically dispersed, and commonly infected from a pool of genetically diverse viruses. Our results suggest that parasite hybridization, likely due to increased human migration and ecological perturbations, increased the frequency of endosymbiotic interactions known to play a key role in disease severity.

Sequencing of hsp70 for discernment of species from the Leishmania (Viannia) guyanensis complex from endemic areas in Colombia

BACKGROUND

Colombia is ranked very high among countries with the highest numbers of endemic Leishmania species (n = 9) causing human disease. Although much effort has been devoted to generating simple and specific tools for Leishmania species identification, challenges remain in the discrimination of species belonging to the Leishmania (Viannia) guyanensis complex: L. (V.) guyanensis and L. (V.) panamensis.

METHODS

A set of seven reference strains of species belonging to the L. (Leishmania) and L. (Viannia) subgenera, clinical strains from human cases of cutaneous leishmaniasis (CL; n = 26) and samples collected from sylvatic mammals and sand flies (n = 7) from endemic areas in Colombia were analyzed in this study. The heat-shock protein 70 gene (hsp70) was amplified by PCR from DNA extracted from logarithmic-phase promastigotes or tissue samples, and the PCR products were sequenced. Sequence alignment was performed against a set of previously published and curated sequences, and phylogenetic analysis based on the maximum-likelihood and Bayesian inference approaches was conducted. Haplotype diversity among strains and species of the L. (V.) guyanensis complex was explored using a median-joining network.

RESULTS

Sequencing of the hsp70 gene for L. (Viannia) spp. typing was comparable to species identification using isoenzyme electrophoresis or monoclonal antibodies. Complete species matching was found, except for one sylvatic sample with an identity yet unsolved. Among the L. (V.) panamensis clinical strains, two distinctive phylogenetic clusters were found to correlate with two different zymodemes: L. (V.) panamensis Z2.2 and Z2.3. Analysis of samples from sylvatic environments identified novel records of naturally infected wild mammal and sand fly species.

CONCLUSIONS

Our results support the adequacy of hsp70 gene sequencing as a single-locus approach for discrimination of L. (Viannia) spp., as well as for exploring the genetic diversity within the L. (V.) guyanensis complex.

Genomic diversity and genetic variation of Leishmania panamensis within its endemic range

Species belonging to the Leishmania (Viannia) subgenus are important causative agents of cutaneous and mucocutaneous leishmaniasis in Central and South America. These parasites possess several distinctive biological features that are influenced by their genetics, population structure, and genome instability. To date, several studies have revealed varying degrees of genetic diversity within Leishmania species. Particularly, in species of the L. (Viannia) subgenus, a generalized high intraspecific genetic diversity has been reported, although, conflicting conclusions have been drawn using different molecular techniques. Despite being the most common Leishmania species circulating in Panama and Colombia, few studies have analyzed clinical samples of Leishmania panamensis using whole-genome sequencing, and their restricted number of samples has limited the information they can provide to understand the population structure of L. panamensis. Here, we used next generation sequencing (NGS) to explore the genetic diversity of L. panamensis within its endemic range, analyzing data from 43 isolates of Colombian and Panamanian origin. Our results show the occurrence of three well-defined geographically correlated groups, and suggests the possible occurrence of additional phylogeographic groups. Furthermore, these results support the existence of a mixed mode of reproduction in L. panamensis, with varying frequencies of events of genetic recombination occurring primarily within subpopulations of closely related strains. This study offers important insights into the population genetics and reproduction mode of L. panamensis, paving the way to better understand their population structure and the emergence and maintenance of key eco-epidemiological traits.

Genetic diversity and population structure of Leishmania (Viannia) braziliensis in the Peruvian jungle

Background

Human cutaneous leishmaniasis caused by Leishmania (Viannia) braziliensis is highly prevalent in the Peruvian jungle, where it affects military forces deployed to fight against drug trafficking and civilian people that migrate from the highland to the lowland jungle for economic activities such as mining, agriculture, construction, and chestnut harvest.

We explored the genetic diversity and population structure of 124 L. (V.) braziliensis isolates collected from the highland (Junín, Cusco, and Ayacucho) and lowland Peruvian jungle (Loreto, Ucayali, and Madre de Dios). All samples were genotyped using Multilocus Microsatellite Typing (MLMT) of ten highly polymorphic markers.

Principal findings

High polymorphism and genetic diversity were found in Peruvian isolates of L. (V.) braziliensis. Most markers are not in Hardy-Weinberg equilibrium; this deviation is most likely caused by local inbreeding, as shown by the positive F_IS values. Linkage Disequilibrium in subpopulations was not strong, suggesting the reproduction was not strictly clonal. Likewise, for the first time, two genetic clusters of this parasite were determined, distributed in both areas of the Peruvian jungle, which suggested a possible recent colonization event of the highland jungle from the lowland jungle.

Conclusions

L. (V.) braziliensis exhibits considerable genetic diversity with two different clusters in the Peruvian jungle. Migration analysis suggested a colonization event between geographical areas of distribution. Although no human migration was observed at the time of sampling, earlier displacement of humans, reservoirs, or vectors could have been responsible for the parasite spread in both regions.

L. (V.) braziliensis is widespread in the Peruvian jungle region. In this region, the departments of Cusco and Madre de Dios account for a large number of patients that get infected while working in the virgin forest. For the first time, we described ample genetic diversity among Peruvian L. (V.) braziliensis isolates with new alleles that were not previously reported in South America. In addition, two different genetic clusters or subpopulations of L. (V.) braziliensis in the Amazon Rainforest of Peru were described. This finding reveals the important distribution of parasite populations and suggests a possible colonization event between ecoregions of the highland and lowland Peruvian jungle independently of recent human migration.

Leishmania and the Model of Predominant Clonal Evolution

As it is the case for other pathogenic microorganisms, the respective impact of clonality and genetic exchange on Leishmania natural populations has been the object of lively debates since the early 1980s. The predominant clonal evolution (PCE) model states that genetic exchange in these parasites’ natural populations may have a high relevance on an evolutionary scale, but is not sufficient to erase a persistent phylogenetic signal and the existence of bifurcating trees. Recent data based on high-resolution markers and genomic polymorphisms fully confirm the PCE model down to a microevolutionary level.

Diagnostic Efficacy of Recombinase-Polymerase-Amplification Coupled with Lateral Flow Strip Reading in Patients with Cutaneous Leishmaniasis from the Amazonas Rainforest of Perú

Cutaneous leishmaniasis (CL) is highly prevalent in rural and sylvatic regions of Latin America, with an estimated 55,000 annual cases. Diagnosis in resource-limited areas still relies on microscopy of dermal scrapings, while more sensitive methods like PCR are not attainable due to costs and lack of adequate health infrastructure. Isothermal amplification of Leishmania DNA can be performed without sophisticated equipment and training and may become a point of care (POC) test for health care centers with scarce resources. We evaluated the efficacy of recombinase-polymerase-amplification (RPA-LF) to diagnose CL in 226 patients attending a clinic in Puerto Maldonado within the Peruvian Amazon basin. Conventional PCR targeting kinetoplast DNA (kDNA-PCR) was used as the gold standard. Eight of 226 patients were considered true negatives (microscopy, kDNA-PCR, and RPA-LF negative), while RPA-LF resulted positive in 186 of 204 kDNA-PCR positive patients, yielding 91.2% (confidence interval [CI] = 86.5-94.4%) sensitivity and 93% (CI 88.6-95.8%) positive predictive value. There were 14% (32/226) discrepant samples alternating positive and negative results in similar proportions between both tests. Quantitative PCR used to resolve the discrepancies suggested that they occurred in samples with scarce parasite numbers as determined by high cycle threshold (Ct) values (≥32; cutoff 35.5). Microscopy had the lowest sensitivity of all methods (45.4%). Nested real-time PCR performed in 71 samples determined that Leishmania (Viannia) braziliensis was highly prevalent (69/71), and Leishmania (Viannia) lainsoni was present in only two isolates. Results indicated that RPA-LF has POC potential for CL endemic areas, yet further simplification and optimization coupled with field validation will be necessary to confirm its broad applicability.

Tegumentary leishmaniasis by Leishmania braziliensis complex in Cochabamba, Bolivia including the presence of L. braziliensis outlier: Tegumentary leishmaniasis in Cochabamba, Bolivia

Leishmaniasis is caused by protozoans of the Leishmania genus, which includes more than 20 species capable of infecting humans worldwide. In the Americas, the most widespread specie is L. braziliensis, present in 18 countries including Bolivia. The taxonomic position of the L. braziliensis complex has been a subject of controversy, complicated further by the recent identification of a particular subpopulation named L. braziliensis atypical or outlier. The aim of this study was to carry out a systematic analysis of the L. braziliensis complex in Bolivia and to describe the associated clinical characteristics. Forty-one strains were analyzed by sequencing an amplified 1245 bp fragment of the hsp70 gene, which allowed its identification as: 24 (59%) L. braziliensis, 16 (39%) L. braziliensis outlier, and one (2%) L. peruviana. In a dendrogram constructed, L. braziliensis and L. peruviana are grouped in the same cluster, whilst L. braziliensis outlier appears in a separate branch. Sequence alignment allowed the identification of five non-polymorphic nucleotide positions (288, 297, 642, 993, and 1213) that discriminate L. braziliensis and L. peruviana from L. braziliensis outlier. Moreover, nucleotide positions 51 and 561 enable L. peruviana to be discriminated from the other two taxa. A greater diversity was observed in L. braziliensis outlier than in L. braziliensis-L. peruviana. The 41 strains came from 32 patients with tegumentary leishmaniasis, among which 22 patients (69%) presented cutaneous lesions (11 caused by L. braziliensis and 11 by L. braziliensis outlier) and 10 patients (31%) mucocutaneous lesions (eight caused by L. braziliensis, one by L. braziliensis outlier, and one by L. peruviana). Nine patients (28%) simultaneously provided two isolates, each from a separate lesion, and in each case the same genotype was identified in both. Treatment failure was observed in six patients infected with L. braziliensis and one patient with L. peruviana.

The Maze Pathway of Coevolution: A Critical Review over the Leishmania and Its Endosymbiotic History

The description of the genus Leishmania as the causative agent of leishmaniasis occurred in the modern age. However, evolutionary studies suggest that the origin of Leishmania can be traced back to the Mesozoic era. Subsequently, during its evolutionary process, it achieved worldwide dispersion predating the breakup of the Gondwana supercontinent. It is assumed that this parasite evolved from monoxenic Trypanosomatidae. Phylogenetic studies locate dixenous Leishmania in a well-supported clade, in the recently named subfamily Leishmaniinae, which also includes monoxenous trypanosomatids. Virus-like particles have been reported in many species of this family. To date, several Leishmania species have been reported to be infected by Leishmania RNA virus (LRV) and Leishbunyavirus (LBV). Since the first descriptions of LRVs decades ago, differences in their genomic structures have been highlighted, leading to the designation of LRV1 in L. (Viannia) species and LRV2 in L. (Leishmania) species. There are strong indications that viruses that infect Leishmania spp. have the ability to enhance parasitic survival in humans as well as in experimental infections, through highly complex and specialized mechanisms. Phylogenetic analyses of these viruses have shown that their genomic differences correlate with the parasite species infected, suggesting a coevolutionary process. Herein, we will explore what has been described in the literature regarding the relationship between Leishmania and endosymbiotic Leishmania viruses and what is known about this association that could contribute to discussions about the worldwide dispersion of Leishmania.

Molecular Identification of Parasites Causing Cutaneous Leishmaniasis in Panama

Isolates from 475 cutaneous leishmaniasis (CL) patients from three endemic regions were studied by three typing techniques. The molecular analysis from lesion scrapings based on hsp70 PCR-RFLP showed that 78.1% (371/475) restriction patterns corresponded to Leishmania (Viannia) panamensis, 19% (90/475) to Leishmania (Viannia) guyanensis, and 3.0% (14/475) to Leishmania (Viannia) braziliensis. Promastigotes isolated by culture from lesions of 228 patients (48.0%, 228/475) were identified by multi-locus enzyme electrophoresis. Of them, 95.2% (217/228) were typified as L. (V.) panamensis, 1.3% (3/228) as L. (V.) guyanensis, 2.2% (5/228) as L. (V.) braziliensis, and 1.3% (3/228) as hybrids (L. [V.] braziliensis/L. [V.] panamensis). However, a partial sequencing analysis of the hsp70 gene from 77 selected samples showed 16.9% (13/77) typified as L. (V.) panamensis, 68.8% (53/77) as Leishmania (V.) sp., 1, 3.9% (3/77) as L. (V.) guyanensis, 1.3% (1/77) as L. (V.) braziliensis outlier, 2.6% (2/77) as Leishmania (Viannia) naiffi, 2.6% as (2/77) Leishmania (V.) sp., and 2 and 3.9% (3/77) hybrid isolates of L. (V.) braziliensis/L. (V.) guyanensis. These results confirm L. (V.) panamensis as the predominant species and cause of CL lesions in Panama and that L. (V.) guyanensis, L. (V.) braziliensis, and L. (V.) naiffi are circulating to a lower degree. Furthermore, the determination of parasite isolates belonging to atypical clusters and hybrid isolates suggests the circulation of genetic variants with important implications for the epidemiology and clinical follow-up of CL in Panama. No evidence of the existence of parasites of the Leishmania (Leishmania) mexicana complex in Panamanian territory was found in this study.

Distribution, treatment outcome and genetic diversity of Leishmania species in military personnel from Colombia with cutaneous leishmaniasis

BACKGROUND

Leishmaniasis is one of the most important infectious diseases affecting the Colombian National Army due to the high number of reported cases and exposure throughout military operations in endemic areas. The main aim of this study was to estimate the geographical distribution along with the genetic diversity and treatment outcome of Leishmania species in Colombian military personnel.

METHODS

Skin lesion samples by smear and aspirate were collected in 136 patients having parasitological cutaneous leishmaniasis (CL) diagnosis. DNA was extracted, the nuclear marker heat shock protein 70 (HSP70) was amplified by PCR and sequenced. Leishmania species were identified by BLASTn. The geo-spatial distribution of the identified parasites was determined according to the possible site of infection. Gene tree was constructed by maximum likelihood (ML), diversity indices (π, h) were estimated and haplotype network was constructed under the Templeton-Crandall-Sing algorithm in order to determine the geographic relationships of the genetic variants of Leishmania species circulating in Colombian military population.

RESULTS

The species were identified in 77.94% of the samples, with a predominance of L. braziliensis (65.09%), followed by L. panamensis (31.13%), L. naiffi by the first time reported in Colombia in two patients (1.89%) as well as L. lindenbergi in a single patient (0.945%) with possible infection in the municipality of Miraflores, Guaviare and L. infantum in a single patient (0.945%) notified with CL in the municipality of Tumaco, Nariño. The phylogenetic analysis was consistent according to bootstrap, showing four strongly differentiated clades.

CONCLUSIONS

The geo-spatial distribution suggested that L. braziliensis has a greater abundance, while L. panamensis has a greater dispersion. The phylogenetic relationships of Leishmania species in Colombian military personnel was estimated with the confirmation of two new species circulating without prior report in the country and a species with no background for CL in the Colombian army. A substantial genetic diversity of Leishmania braziliensis was defined. This study contributes through the understanding of the molecular epidemiology to the CL transmission in Colombia.

Ecological divergence and hybridization of Neotropical Leishmania parasites

Parasites are interesting models for studying speciation processes because they have a high potential for specialization, thanks to the intimate ecological association with their hosts and vectors. Yet little is known about the circumstances under which new parasite lineages emerge. Here we studied the genome diversity of parasites of the Leishmania braziliensis species complex that inhabit both Amazonian and Andean biotas in Peru. We identify three major parasite lineages that occupy particular ecological niches and show that these emerged during forestation changes over the past 150,000 y. We furthermore discovered that meiotic recombination between Amazonian and Andean lineages resulted in full-genome hybrids presenting mixed mitochondrial genomes, providing insights into the genetic consequences of hybridization in parasitic protozoa.

The tropical Andes are an important natural laboratory to understand speciation in many taxa. Here we examined the evolutionary history of parasites of the Leishmania braziliensis species complex based on whole-genome sequencing of 67 isolates from 47 localities in Peru. We first show the origin of Andean Leishmania as a clade of near-clonal lineages that diverged from admixed Amazonian ancestors, accompanied by a significant reduction in genome diversity and large structural variations implicated in host–parasite interactions. Within the Andean species, patterns of population structure were strongly associated with biogeographical origin. Molecular clock and ecological niche modeling suggested that the history of diversification of the Andean lineages is limited to the Late Pleistocene and intimately associated with habitat contractions driven by climate change. These results suggest that changes in forestation over the past 150,000 y have influenced speciation and diversity of these Neotropical parasites. Second, genome-scale analyses provided evidence of meiotic-like recombination between Andean and Amazonian Leishmania species, resulting in full-genome hybrids. The mitochondrial genome of these hybrids consisted of homogeneous uniparental maxicircles, but minicircles originated from both parental species. We further show that mitochondrial minicircles—but not maxicircles—show a similar evolutionary pattern to the nuclear genome, suggesting that compatibility between nuclear-encoded mitochondrial genes and minicircle-encoded guide RNA genes is essential to maintain efficient respiration. By comparing full nuclear and mitochondrial genome ancestries, our data expand our appreciation on the genetic consequences of diversification and hybridization in parasitic protozoa.

Leishmania naiffi and lainsoni in French Guiana: Clinical features and phylogenetic variability

In French Guiana, five species are associated with Cutaneous Leishmaniasis (CL). Though infections with Leishmania guyanensis, L. (V.) braziliensis and L. (L.) amazonensis have been extensively described, there are few available clinical and genetic data on L. (V.) lainsoni and L. (V.) naiffi. We determined the clinical and epidemiological features of all cases of CL due to L. (V.) naiffi and L. (V.) lainsoni diagnosed in French Guiana between 2003 and 2019. Phylogenetic analysis was performed by sequencing a portion of HSP70 and cyt b genes. Five cases of L. naiffi and 25 cases of L. lainsoni were reported. Patients infected by L. (V.) lainsoni were usually infected on gold camps, mostly along the Maroni river (60%), while L. naiffi was observed in French patients infected on the coast (100%). A high number of pediatric cases (n = 5; 20%) was observed for L. (V.) lainsoni. A mild clinical course was observed for all cases of L. (V.) naiffi. HSP70 and cyt b partial nucleotide sequence analysis revealed different geographical clusters within L. (V.) naiffi and L. (V.) lainsoni but no association were found between phylogenetic and clinical features. Our data suggest distinct socio-epidemiological features for these two Leishmania species. Patients seem to get infected with L. (V.) naiffi during leisure activities in anthropized coastal areas, while L. (V.) lainsoni shares common features with L. (V.) guyanensis and braziliensis and seems to be acquired during professional activities in primary forest regions. Phylogenetic analysis has provided information on the intraspecific genetic variability of L. (V.) naiffi and L. (V.) lainsoni and how these genotypes are distributed at the geographic level.

Leishmania (Viannia) braziliensis type 2 as probable etiological agent of canine cutaneous leishmaniasis in Brazilian Amazon

Canine cutaneous leishmaniasis (CCL) is a zoonosis of public health interest, and in the Americas, Leishmania (Viannia) braziliensis has been identified as the main etiological agent. The present study sought to investigate Leishmania spp. infection in domestic dogs from a rural area of the Xapuri municipality, Acre state, Brazilian Amazonia. For this purpose, visits were carried out to domiciles where the human cases of American cutaneous leishmaniasis (ACL) occurred, followed by the clinical evaluation of the animals in search of clinical signs suggestive of CCL. Blood samples were collected from 40 dogs, 13 of which had lesions suggestive of CCL, and biopsies of these lesions were performed. The methods used were Neal, Novy, and Nicolle’s (NNN) medium cultures and direct parasitological examination. Further, to detect and characterize Leishmania DNA some molecular techniques were performed such as conventional polymerase chain reaction (PCR) and sequencing targeting SSU rDNA and ITS1, restriction fragment length polymorphism (RFLP) and high resolution melting (HRM) analysis targeting hsp70. The investigation revealed that the results obtained from the parasitological methods were negative. In PCR by ITS1 and network topology sequences, six strains from dogs, isolated from the Peruvian Andes, appeared identical to Leishmania (Viannia) braziliensis type 2 (99–100%). By other molecular methods these samples turned out to be positive to Leishmania (Viannia) sp.. The diagnosis of Leishmania in domestic dogs from Acre state showed a high proportion of infected animals, and the occurrence of L. braziliensis type 2 in Brazil for the first time. This new report suggests that L. braziliensis type 2 is both trans- and cis-Andean. However, more studies are needed regarding the clinical and diagnostic aspects of this species of Leishmania.

Phylogenetic Studies

Phylogenetics is an important component of the systems biology approach. Knowledge about evolution of the genus Leishmania is essential to understand various aspects of basic biology of these parasites, such as parasite-host or parasite-vector relationships, biogeography, or epidemiology. Here, we present a comprehensive guideline for performing phylogenetic studies based on DNA sequence data, but with principles that can be adapted to protein sequences or other molecular markers. It is presented as a compilation of the most commonly used genetic targets for phylogenetic studies of Leishmania, including their respective primers for amplification and references, as well as details of PCR assays. Guidelines are, then, presented to choose the best targets in relation to the types of samples under study. Finally, and importantly, instructions are given to obtain optimal sequences, alignments, and datasets for the subsequent data analysis and phylogenetic inference. Different bioinformatics methods and software for phylogenetic inference are presented and explained. This chapter aims to provide a compilation of methods and generic guidelines to conduct phylogenetics of Leishmania for nonspecialists.

Analysis of genetic polymorphisms and tropism in East African Leishmania donovani by Amplified Fragment Length Polymorphism and kDNA minicircle sequencing

Visceral leishmaniasis (VL), the most severe form of leishmaniasis, is caused by Leishmania donovani. In addition to fatal VL, these parasites also cause skin diseases in immune-competent and -suppressed people, post-kala azar dermal leishmaniasis (PKDL) and HIV/VL co-infections, respectively. Genetic polymorphism in 36 Ethiopian and Sudanese L. donovani strains from VL, PKDL and HIV/VL patients was examined using Amplified Fragment Length Polymorphism (AFLP), kDNA minicircle sequencing and Southern blotting. Strains were isolated from different patient tissues: in VL from lymph node, spleen or bone marrow; and in HIV/VL from skin, spleen or bone marrow. When VL and PKDL strains from the same region in Sudan were examined by Southern blotting using a DNA probe to the L. donovani 28S rRNA gene only minor differences were observed. kDNA sequence analysis distributed the strains in no particular order among four clusters (A - D), while AFLP analysis grouped the strains according to geographical origin into two major clades, Southern Ethiopia (SE) and Sudan/Northern Ethiopia (SD/NE). Strains in the latter clade were further divided into subpopulations by zymodeme, geography and year of isolation, but not by clinical symptoms. However, skin isolates showed significantly (p < 0.0001) fewer polymorphic AFLP fragments (average 10 strains = 348.6 ± 8.1) than VL strains (average 26 strains = 383.5 ± 3.8).

Andean cutaneous leishmaniasis (Andean-CL, uta) in Peru and Ecuador: the causative Leishmania parasites and clinico-epidemiological features

This study provides comprehensive information on the past and current status of the Andean cutaneous leishmaniasis (Andean-CL, uta) in Peru and Ecuador, mainly focusing on the causative Leishmania parasites and clinico-epidemiological features. Available information and data including our unpublished works were analyzed thoroughly. Endemic regions of the Andean-CL (uta) in Peru run from the north Piura/Cajamarca to the south Ayacucho at a wide range of the Pacific watersheds of the Andes through several departments, while in Ecuador those exist at limited and spotted areas in the country's mid-southwestern two provinces, Azuay and Chimborazo. The principal species of the genus Leishmania are completely different at subgenus level, L. (Viannia) peruviana in Peru, and L. (Leishmania) mexicana and L. (L.) major-like (infrequent occurrence) in Ecuador. The Peruvian uta is now prevalent in different age and sex groups, being not clearly defined as found in the past. The precise reasons are not known and should be elucidated further, though probable factors, such as emergence of other Leishmania parasites, non-immune peoples' migration into the areas, etc., were discussed briefly in the text. The Andean-CL cases in Ecuador are more rural than before, probably because of a rapid development of the Leishmania-positive communities and towns, and the change of life-styles of the inhabitants, including newly constructed houses and roads in the endemic areas. Such information is helpful for future management of the disease, not only for Leishmania-endemic areas in the Andes but also for other endemic areas.

Identification of Leishmania by Matrix-Assisted Laser Desorption Ionization-Time of Flight (MALDI-TOF) Mass Spectrometry Using a Free Web-Based Application and a Dedicated Mass-Spectral Library

Human leishmaniases are widespread diseases with different clinical forms caused by about 20 species within the Leishmania genus. Leishmania species identification is relevant for therapeutic management and prognosis, especially for cutaneous and mucocutaneous forms. Several methods are available to identify Leishmania species from culture, but they have not been standardized for the majority of the currently described species, with the exception of multilocus enzyme electrophoresis. Moreover, these techniques are expensive, time-consuming, and not available in all laboratories. Within the last decade, mass spectrometry (MS) has been adapted for the identification of microorganisms, including Leishmania However, no commercial reference mass-spectral database is available. In this study, a reference mass-spectral library (MSL) for Leishmania isolates, accessible through a free Web-based application (mass-spectral identification [MSI]), was constructed and tested. It includes mass-spectral data for 33 different Leishmania species, including species that infect humans, animals, and phlebotomine vectors. Four laboratories on two continents evaluated the performance of MSI using 268 samples, 231 of which were Leishmania strains. All Leishmania strains, but one, were correctly identified at least to the complex level. A risk of species misidentification within the Leishmania donovani, L. guyanensis, and L. braziliensis complexes was observed, as previously reported for other techniques. The tested application was reliable, with identification results being comparable to those obtained with reference methods but with a more favorable cost-efficiency ratio. This free online identification system relies on a scalable database and can be implemented directly in users' computers.

Use of AFLP for the study of eukaryotic pathogens affecting humans

Amplified fragment length polymorphism (AFLP) is a genotyping technique based on PCR amplification of specific restriction fragments from a particular genome. The methodology has been extensively used in plant biology to solve a variety of scientific questions, including taxonomy, molecular epidemiology, systematics, population genetics, among many others. The AFLP share advantages and disadvantages with other types of molecular markers, being particularly useful in organisms with no previous DNA sequence knowledge. In eukaryotic pathogens, the technique has not been extensively used, although it has the potential to solve many important issues as it allows the simultaneous examination of hundreds or even thousands of polymorphic sites in the genome of the organism. Here we describe the main applications published on the use of AFLP in eukaryotic pathogens, with emphasis in species of the groups fungi, protozoa and helminths, and discuss the role of this methodology in the context of new techniques derived from the advances of the next generation sequencing.

Comparison of Leishmania typing results obtained from 16 European clinical laboratories in 2014

Leishmaniasis is endemic in southern Europe, and in other European countries cases are diagnosed in travellers who have visited affected areas both within the continent and beyond. Prompt and accurate diagnosis poses a challenge in clinical practice in Europe. Different methods exist for identification of the infecting Leishmania species. Sixteen clinical laboratories in 10 European countries, plus Israel and Turkey, conducted a study to assess their genotyping performance. DNA from 21 promastigote cultures of 13 species was analysed blindly by the routinely used typing method. Five different molecular targets were used, which were analysed with PCR-based methods. Different levels of identification were achieved, and either the Leishmania subgenus, species complex, or actual species were reported. The overall error rate of strains placed in the wrong complex or species was 8.5%. Various reasons for incorrect typing were identified. The study shows there is considerable room for improvement and standardisation of Leishmania typing. The use of well validated standard operating procedures is recommended, covering testing, interpretation, and reporting guidelines. Application of the internal transcribed spacer 1 of the rDNA array should be restricted to Old World samples, while the heat-shock protein 70 gene and the mini-exon can be applied globally.

Is Predominant Clonal Evolution a Common Evolutionary Adaptation to Parasitism in Pathogenic Parasitic Protozoa, Fungi, Bacteria, and Viruses?

We propose that predominant clonal evolution (PCE) in microbial pathogens be defined as restrained recombination on an evolutionary scale, with genetic exchange scarce enough to not break the prevalent pattern of clonal population structure. The main features of PCE are (1) strong linkage disequilibrium, (2) the widespread occurrence of stable genetic clusters blurred by occasional bouts of genetic exchange ('near-clades'), (3) the existence of a "clonality threshold", beyond which recombination is efficiently countered by PCE, and near-clades irreversibly diverge. We hypothesize that the PCE features are not mainly due to natural selection but also chiefly originate from in-built genetic properties of pathogens. We show that the PCE model obtains even in microbes that have been considered as 'highly recombining', such as Neisseria meningitidis, and that some clonality features are observed even in Plasmodium, which has been long described as panmictic. Lastly, we provide evidence that PCE features are also observed in viruses, taking into account their extremely fast genetic turnover. The PCE model provides a convenient population genetic framework for any kind of micropathogen. It makes it possible to describe convenient units of analysis (clones and near-clades) for all applied studies. Due to PCE features, these units of analysis are stable in space and time, and clearly delimited. The PCE model opens up the possibility of revisiting the problem of species definition in these organisms. We hypothesize that PCE constitutes a major evolutionary strategy for protozoa, fungi, bacteria, and viruses to adapt to parasitism.

Calmodulin Polymerase Chain Reaction-Restriction Fragment Length Polymorphism for Leishmania Identification and Typing

A precise identification of Leishmania species involved in human infections has epidemiological and clinical importance. Herein, we describe a preliminary validation of a restriction fragment length polymorphism assay, based on the calmodulin intergenic spacer region, as a tool for detecting and typing Leishmania species. After calmodulin amplification, the enzyme HaeIII yielded a clear distinction between reference strains of Leishmania mexicana, Leishmania amazonensis, Leishmania infantum, Leishmania lainsoni, and the rest of the Viannia reference species analyzed. The closely related Viannia species: Leishmania braziliensis, Leishmania panamensis, and Leishmania guyanensis, are separated in a subsequent digestion step with different restriction enzymes. We have developed a more accessible molecular protocol for Leishmania identification/typing based on the exploitation of part of the calmodulin gene. This methodology has the potential to become an additional tool for Leishmania species characterization and taxonomy.

Molecular identification of Leishmania spp. clinical isolates from Colombia based on hsp70 gene

INTRODUCTION

Leishmaniasis is highly prevalent in Colombia, where at least six different species can cause disease of varying clinical presentations in humans. The identification of the infecting species is quite important for prognosis, therapeutics and epidemiology. Different techniques with variable discriminatory power have been used for the identification. 

OBJECTIVE

To carry out the molecular identification of Leishmania species through the amplification of a fragment of the hsp70 gene. 

MATERIALS AND METHODS

Molecular amplification of the hsp70 gene fragment (PCR-hsp70) followed by restriction fragment length polymorphism analysis (RFLP) was done for identification purposes using DNA from 81 clinical isolates of Leishmania. 

RESULTS

A single amplicon was obtained for all samples analyzed. The enzymatic restrictions of the 81 PCR products identified 70 with a banding pattern corresponding to L. braziliensis with two different patterns (62 and eight isolates, respectively), nine isolates compatible with L. panamensis and two with L. guyanensis. The geographical origin of the isolates is consistent with previous reports about the distribution of the corresponding species in Colombia. 

CONCLUSIONS

The PCR-hsp70/RFLP technique used is a valid tool for the identification of Leishmania species isolated from clinical samples of patients in Colombia, which may also be applicable to the study of strains obtained from vectors and reservoirs with epidemiological significance.

Comparative genomic analysis of Leishmania (Viannia) peruviana and Leishmania (Viannia) braziliensis

BACKGROUND

The Leishmania (Viannia) braziliensis complex is responsible for most cases of New World tegumentary leishmaniasis. This complex includes two closely related species but with different geographic distribution and disease phenotypes, L. (V.) peruviana and L. (V.) braziliensis. However, the genetic basis of these differences is not well understood and the status of L. (V.) peruviana as distinct species has been questioned by some. Here we sequenced the genomes of two L. (V.) peruviana isolates (LEM1537 and PAB-4377) using Illumina high throughput sequencing and performed comparative analyses against the L. (V.) braziliensis M2904 reference genome. Comparisons were focused on the detection of Single Nucleotide Polymorphisms (SNPs), insertions and deletions (INDELs), aneuploidy and gene copy number variations.

RESULTS

We found 94,070 variants shared by both L. (V.) peruviana isolates (144,079 in PAB-4377 and 136,946 in LEM1537) against the L. (V.) braziliensis M2904 reference genome while only 26,853 variants separated both L. (V.) peruviana genomes. Analysis in coding sequences detected 26,750 SNPs and 1,513 indels shared by both L. (V.) peruviana isolates against L. (V.) braziliensis M2904 and revealed two L. (V.) braziliensis pseudogenes that are likely to have coding potential in L. (V.) peruviana. Chromosomal read density and allele frequency profiling showed a heterogeneous pattern of aneuploidy with an overall disomic tendency in both L. (V.) peruviana isolates, in contrast with a trisomic pattern in the L. (V.) braziliensis M2904 reference. Read depth analysis allowed us to detect more than 368 gene expansions and 14 expanded gene arrays in L. (V.) peruviana, and the likely absence of expanded amastin gene arrays.

CONCLUSIONS

The greater numbers of interspecific SNP/indel differences between L. (V.) peruviana and L. (V.) braziliensis and the presence of different gene and chromosome copy number variations support the classification of both organisms as closely related but distinct species. The extensive nucleotide polymorphisms and differences in gene and chromosome copy numbers in L. (V.) peruviana suggests the possibility that these may contribute to some of the unique features of its biology, including a lower pathology and lack of mucosal development.

Association of the Endobiont Double-Stranded RNA Virus LRV1 With Treatment Failure for Human Leishmaniasis Caused by Leishmania braziliensis in Peru and Bolivia

Cutaneous and mucosal leishmaniasis, caused in South America by Leishmania braziliensis, is difficult to cure by chemotherapy (primarily pentavalent antimonials [Sb(V)]). Treatment failure does not correlate well with resistance in vitro, and the factors responsible for treatment failure in patients are not well understood. Many isolates of L. braziliensis (>25%) contain a double-stranded RNA virus named Leishmaniavirus 1 (LRV1), which has also been reported in Leishmania guyanensis, for which an association with increased pathology, metastasis, and parasite replication was found in murine models. Here we probed the relationship of LRV1 to drug treatment success and disease in 97 L. braziliensis-infected patients from Peru and Bolivia. In vitro cultures were established, parasites were typed as L. braziliensis, and the presence of LRV1 was determined by reverse transcription-polymerase chain reaction, followed by sequence analysis. LRV1 was associated significantly with an increased risk of treatment failure (odds ratio, 3.99; P = .04). There was no significant association with intrinsic Sb(V) resistance among parasites, suggesting that treatment failure arises from LRV1-mediated effects on host metabolism and/or parasite survival. The association of LRV1 with clinical drug treatment failure could serve to guide more-effective treatment of tegumentary disease caused by L. braziliensis.

Species typing in dermal leishmaniasis

Leishmania is an infectious protozoan parasite related to African and American trypanosomes. All Leishmania species that are pathogenic to humans can cause dermal disease. When one is confronted with cutaneous leishmaniasis, identification of the causative species is relevant in both clinical and epidemiological studies, case management, and control. This review gives an overview of the currently existing and most used assays for species discrimination, with a critical appraisal of the limitations of each technique. The consensus taxonomy for the genus is outlined, including debatable species designations. Finally, a numerical literature analysis is presented that describes which methods are most used in various countries and regions in the world, and for which purposes.

A primer for Leishmania population genetic studies

Leishmaniases remain a major public health problem. Despite the development of elaborate experimental techniques and sophisticated statistical tools, how these parasites evolve, adapt themselves to new environmental compartments and hosts, and develop resistance to new drugs remains unclear. Leishmania parasites constitute a complex model from a biological, ecological, and epidemiological point of view but also with respect to their genetics and phylogenetics. With this in view, we seek to outline the criteria, caveats, and confounding factors to be considered for Leishmania population genetic studies. We examine how the taxonomic complexity, heterozygosity, intraspecific and interspecific recombination, aneuploidy, and ameiotic recombination of Leishmania intersect with population genetic studies of this parasite.

Evaluation of four single-locus markers for Leishmania species discrimination by sequencing

Several genetic markers have been described for discriminating Leishmania species. In most reported cases, one or a few polymorphisms are the basis of species identification, and the methods were validated on a limited number of strains from a particular geographical region. Therefore, most techniques may underestimate the global intraspecies variability and are applicable only in certain areas. In addition, interlaboratory standardization is mostly absent, complicating comparisons among different studies. Here, we compared species typing results from all sequence polymorphisms found in four popular markers that can be applied directly on clinical samples: the miniexon or spliced leader, the internal transcribed spacer of the ribosomal DNA array, the 7SL RNA gene, and the heat shock protein 70 gene. Clustering was evaluated among 74 Leishmania strains, selected to represent a wide geographic distribution and genetic variability of the medically relevant species of the genus. Results were compared with a multilocus sequence typing (MLST) approach using 7 single-copy household genes and with multilocus enzyme electrophoresis (MLEE), still considered the gold standard by some. We show that strain groupings are highly congruent across the four different single-locus markers, MLST, and MLEE. Overall, the heat shock protein 70 gene and the miniexon presented the best resolutions for separating medically relevant species. As gene sequence analysis is validated here on a global scale, it is advocated as the method of choice for use in genetic, clinical, and epidemiological studies and for managing patients with unknown origins of infection, especially in Western infectious disease clinics dealing with imported leishmaniasis.

The calmodulin intergenic spacer as molecular target for characterization of Leishmania species

Background

Human leishmaniasis is a neglected disease caused by parasites of the genus Leishmania. Clinical aspects of this disease can vary significantly, reflecting the wide range of parasites in the genus Leishmania. Knowing accurately the Leishmania species infecting humans is important for clinical case management and evaluation of epidemiological risk. Calmodulin is an essential gene in trypanosomatids that modulates the calcium metabolism in various cellular activities. Despite its strong conservation in trypanosomatids, it has been recently observed that its untranslated regions (UTR) diverge among species.

Methods

In this study we analyzed the sequences and the absolute dinucleotide frequency of the intergenic spacer of the calmodulin gene (containing both, 3′ and 5′UTR) in nine reference Leishmania species and ten clinical isolates obtained from patients with cutaneous leishmaniasis.

Results

We show that the short calmodulin intergenic spacers exhibit features that make them interesting for applications in molecular characterization and phylogenetic studies of Leishmania. Dendrograms based on sequence alignments and on the dinucleotide frequency indicate that this particular region of calmodulin gene might be useful for species typing between the Leishmania and Viannia subgenera.

Conclusions

Mutations and composition of the calmodulin intergenic spacer from Leishmania species might have taxonomic value as parameters to define if an isolate is identical to a certain species or belongs to one of the two current subgenera.

Population structure and evidence for both clonality and recombination among Brazilian strains of the subgenus Leishmania (Viannia)

Background/Objectives: Parasites of the subgenus Leishmania (Viannia) cause varying clinical symptoms ranging from cutaneous leishmaniases (CL) with single or few lesions, disseminated CL (DL) with multiple lesions to disfiguring forms of mucocutaneous leishmaniasis (MCL). In this population genetics study, 37 strains of L. (V.) guyanensis, 63 of L. (V.) braziliensis, four of L. (V.) shawi, six of L. (V.) lainsoni, seven of L. (V.) naiffi, one each of L. (V.) utingensis and L. (V.) lindenbergi, and one L. (V.) lainsoni/L. naiffi hybrid from different endemic foci in Brazil were examined for variation at 15 hyper-variable microsatellite markers.

Methodology/Principal findings: The multilocus microsatellite profiles obtained for the 120 strains were analysed using both model- and distance-based methods. Significant genetic diversity was observed for all L. (Viannia) strains studied. The two cluster analysis approaches identified two principal genetic groups or populations, one consisting of strains of L. (V.) guyanensis from the Amazon region and the other of strains of L. (V.) braziliensis isolated along the Atlantic coast of Brazil. A third group comprised a heterogeneous assembly of species, including other strains of L. braziliensis isolated from the north of Brazil, which were extremely polymorphic. The latter strains seemed to be more closely related to those of L. (V.) shawi, L. (V.) naiffi, and L. (V.) lainsoni, also isolated in northern Brazilian foci. The MLMT approach identified an epidemic clone consisting of 13 strains of L. braziliensis from Minas Gerais, but evidence for recombination was obtained for the populations of L. (V.) braziliensis from the Atlantic coast and for L. (V.) guyanensis.

Conclusions/Significance: Different levels of recombination versus clonality seem to occur within the subgenus L. (Viannia). Though clearly departing from panmixia, sporadic, but long-term sustained recombination might explain the tremendous genetic diversity and limited population structure found for such L. (Viannia) strains.

Cutaneous leishmaniasis (CL) constitutes a significant public health problem in all federal states of Brazil. Most cases are caused by parasites of the subgenus Leishmania (Viannia) which can cause a variety of clinical symptoms ranging from single or few lesions, disseminated CL with multiple lesions, to disfiguring forms of mucocutaneous leishmaniasis. This study has used a multilocus microsatellite typing approach for exploring the genetic diversity and population structure among 120 strains representing different subgenus L. (Viannia) species and different Brazilian CL foci. Genetic diversity within the subgenus was much higher than expected, especially within L. (V.) braziliensis, L. (V.) shawi, L. (V.) naiffi, and L. (V.) lainsoni which were all from the north of Brazil. These strains could not be assigned to well-defined populations, but presented a rather loosely associated group. Strains of L. (V.) braziliensis isolated along the Atlantic coast of Brazil and strains of L. (V.) guyanensis formed, however, two clearly separated populations exhibiting remarkable levels of sexual exchange. The latter finding is in contrast to previous studies suggesting clonal modes of propagation or inbreeding for natural populations of Leishmania parasites and might explain the genetic heterogeneity and limited population structure for Brazilian strains of subgenus L. (Viannia) observed in this study.

HindII and SduI digests of heat-shock protein 70 PCR for Leishmania typing

Restriction fragment length polymorphisms of the heat-shock protein 70 gene have been used for discriminating Leishmania species. Here, we validated HindII as a much cheaper alternative to EcoRII and SduI for discriminating Leishmania (Viannia) braziliensis from Leishmania (Viannia) naiffi and an atypical Leishmania (V.) braziliensis group, which was previously not possible.

AFLP polymorphisms allow high resolution genetic analysis of American Tegumentary Leishmaniasis agents circulating in Panama and other members of the Leishmania genus

American Tegumentary Leishmaniasis is caused by parasites of the genus Leishmania, and causes significant health problems throughout the Americas. In Panama, Leishmania parasites are endemic, causing thousands of new cases every year, mostly of the cutaneous form. In the last years, the burden of the disease has increased, coincident with increasing disturbances in its natural sylvatic environments. The study of genetic variation in parasites is important for a better understanding of the biology, population genetics, and ultimately the evolution and epidemiology of these organisms. Very few attempts have been made to characterize genetic polymorphisms of parasites isolated from Panamanian patients of cutaneous leishmaniasis. Here we present data on the genetic variability of local isolates of Leishmania, as well as specimens from several other species, by means of Amplified Fragment Length Polymorphisms (AFLP), a technique seldom used to study genetic makeup of parasites. We demonstrate that this technique allows detection of very high levels of genetic variability in local isolates of Leishmania panamensis in a highly reproducible manner. The analysis of AFLP fingerprints generated by unique selective primer combinations in L. panamensis suggests a predominant clonal mode of reproduction. Using fluorescently labeled primers, many taxon-specific fragments were identified which may show potential as species diagnostic fragments. The AFLP permitted a high resolution genetic analysis of the Leishmania genus, clearly separating certain groups among L. panamensis specimens and highly related species such as L. panamensis and L. guyanensis. The phylogenetic networks reconstructed from our AFLP data are congruent with established taxonomy for the genus Leishmania, even when using single selective primer combinations. Results of this study demonstrate that AFLP polymorphisms can be informative for genetic characterization in Leishmania parasites, at both intra and inter-specific levels.

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.

Evolution and species discrimination according to the Leishmania heat-shock protein 20 gene

The Leishmania genus comprises up to 35 species, of which 20 are responsible for human disease. However, the taxonomic status for many of them is under discussion. The small Heat Shock Proteins (sHSPs) are physiologically relevant, protecting cellular proteins from aggregation and maintaining cellular viability under intensive stress conditions. In Leishmania, a protein of this class was previously described, the 20-kDa heat-shock protein (HSP20), which is encoded by a single gene. In the present study, we used this target, alone or in combination with hsp70 gene, to investigate the phylogenetic relationships among Leishmania species. Using a pair of degenerate primers it was possible amplifying a 370bp fragment of the hsp20 coding region in 39 strains of very different geographic origins, representing in total 16 Leishmania species (14 if L. chagasi and L. archibaldi are considered synonymous names of L. infantum and L. donovani, respectively). Nucleotide sequences were readily obtained by direct sequencing of the amplification products. Both phylogenetic trees and networks based on either hsp20 sequences or combined datasets of hsp20 and hsp70 sequences were constructed. These phylogenic analyses supported the division of the Leishmania genus into nine species: L. (L.) donovani, L. (L.) major, L. (L.) tropica, L. (L.) aethiopica, L. (L.) mexicana, L. (V.) lainsoni, L. (V.) naiffi, L. (V.) guyanensis and L. (V.) braziliensis. Additionally, by network analysis, the subspecies L. (L.) donovani infantum and L. (V.) braziliensis peruviana were recognized within the L. (L.) donovani and L. (V.) braziliensis species, respectively. Therefore, hsp20 gene was found to be a suitable molecular marker for Leishmania typing and classification purposes. In addition, this study represents a solid contribution to the objective of establishing a more reliable taxonomy for the genus Leishmania.