Genome-wide discovery and development of polymorphic microsatellites from Leishmania panamensis parasites circulating in central Panama

Insights into the Genetic Diversity of Leishmania (Viannia) panamensis in Panama, Inferred via Multilocus Sequence Typing (MLST)

Leishmaniasis is a disease caused by parasites of the genus Leishmania and transmitted by sand fly vectors. Tegumentary leishmaniasis is the most prevalent clinical outcome in Latin America, afflicting people from 18 countries. In Panama, the annual incidence rate of leishmaniasis is as high as 3000 cases, representing a major public health problem. In endemic regions, L. panamensis is responsible for almost eighty percent of human cases that present different clinical outcomes. These differences in disease outcomes could be the result of the local interplay between L. panamensis variants and human hosts with different genetic backgrounds. The genetic diversity of L. panamensis in Panama has only been partially explored, and the variability reported for this species is based on few studies restricted to small populations and/or with poor resolutive markers at low taxonomic levels. Accordingly, in this study, we explored the genetic diversity of sixty-nine L. panamensis isolates from different endemic regions of Panama, using an MLST approach based on four housekeeping genes (Aconitase, ALAT, GPI and HSP70). Two to seven haplotypes per locus were identified, and regional differences in the genetic diversity of L. panamensis were observed. A genotype analysis evidenced the circulation of thirteen L. panamensis genotypes, a fact that might have important implications for the local control of the disease.

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.

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.

Genomic Diversification, Structural Plasticity, and Hybridization in Leishmania (Viannia) braziliensis

Leishmania (Viannia) braziliensis is an important Leishmania species circulating in several Central and South American countries. Among Leishmania species circulating in Brazil, Argentina and Colombia, L. braziliensis has the highest genomic variability. However, genomic variability at the whole genome level has been only studied in Brazilian and Peruvian isolates; to date, no Colombian isolates have been studied. Considering that in Colombia, L. braziliensis is a species with great clinical and therapeutic relevance, as well as the role of genetic variability in the epidemiology of leishmaniasis, we analyzed and evaluated intraspecific genomic variability of L. braziliensis from Colombian and Bolivian isolates and compared them with Brazilian isolates. Twenty-one genomes were analyzed, six from Colombian patients, one from a Bolivian patient, and 14 Brazilian isolates downloaded from public databases. The results obtained of Phylogenomic analysis showed the existence of four well-supported clades, which evidenced intraspecific variability. The whole-genome analysis revealed structural variations in the somy, mainly in the Brazilian genomes (clade 1 and clade 3), low copy number variations, and a moderate number of single-nucleotide polymorphisms (SNPs) in all genomes analyzed. Interestingly, the genomes belonging to clades 2 and 3 from Colombia and Brazil, respectively, were characterized by low heterozygosity (~90% of SNP loci were homozygous) and regions suggestive of loss of heterozygosity (LOH). Additionally, we observed the drastic whole genome loss of heterozygosity and possible hybridization events in one genome belonging to clade 4. Unique/shared SNPs between and within the four clades were identified, revealing the importance of some of them in biological processes of L. braziliensis. Our analyses demonstrate high genomic variability of L. braziliensis in different regions of South America, mainly in Colombia and suggest that this species exhibits striking genomic diversity and a capacity of genomic hybridization; additionally, this is the first study to report whole-genome sequences of Colombian L. braziliensis isolates.

Intraspecific Genomic Divergence and Minor Structural Variations in Leishmania (Viannia) panamensis

Leishmania (Viannia) panamensis is one of the most important Leishmania species associated with cutaneous leishmaniasis (CL) in Latin America. Despite its wide geographic distribution and pathogenic potential in humans and animals, the genomic variability of this species is low compared with other Leishmania species circulating in the same geographical area. No studies have reported a detailed analysis of the whole genome of L. panamensis from clinical isolates using DNA high-throughput sequencing to clarify its intraspecific genomic variability or plausible divergence. Therefore, this study aimed to evaluate the intraspecific genomic variability of L. panamensis from Colombia and Panama. A total of 22 genomes were analyzed, 19 from Colombian patients with CL and three genomes from Panama obtained from public databases. The phylogenomic analysis revealed the potential existence of three well-supported clades as evidence of intraspecific divergence. Additionally, the whole-genome analysis showed low structural variations in terms of ploidy, copy number variations, and single-nucleotide polymorphisms (SNPs). SNPs shared among all clades were identified, revealing their importance in different biological processes of L. panamensis. The findings not only expand our knowledge of intraspecific genomic variability of one of the most important Leishmania species in South America but also highlights the possible existence of different clades/lineages/subpopulations across a geographic scale.