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.

Detection of Leishmania RNA Virus 1 in Leishmania (Viannia) panamensis Isolates, Panama.

We detected Leishmania RNA virus 1 (LRV1) in 11 isolates of Leishmania (Viannia) panamensis collected during 2014-2019 from patients from different geographic areas in Panama. The distribution suggested a spread of LRV1 in L. (V.) panamensis parasites. We found no association between LRV1 and an increase in clinical pathology.

Genetic Diversity of Trypanosoma cruzi in Panama Inferred by Multi-locus Sequence Typing of Mitochondrial Genes.

The objective of this study was to provide information on Trypanosoma cruzi genetic diversity among isolates obtained from different biological sources circulating in endemic areas of Panama. Initial discrete typing units (DTUs) assignment was performed evaluating three single locus molecular markers (mini-exon, heat shock protein 60 and glucose-6-phosphate isomerase genes). Further diversity within TcI lineages was explored using a multi-locus sequence typing approach with six maxicircle genes. Haplotype network analysis and evolutionary divergency estimations were conducted to investigate the genetic relatedness between Panamanian TcI isolates and isolates from different endemic regions in the Americas. Our molecular approach validated that TcI is the predominant DTU circulating in Panama across different hosts and vector species, but also confirmed the presence of TcIII and TcVI circulating in the country. The phylogenetic tree topography for most Panamanian TcI isolates displayed a high level of genetic homogeneity between them. The haplotype network analysis inferred a higher genetic diversity within Panamanian TcI isolates, displaying eight different haplotypes circulating in endemic regions of the country, and revealed geographical structuring among TcI from different endemic regions in the Americas. This study adds novelty on the genetic diversity of T. cruzi circulating in Panama and complements regional phylogeographic studies regarding intra-TcI variations.

Surveillance and genotype characterization of zoonotic trypanosomatidae in Didelphis marsupialis in two endemic sites of rural Panama.

Didelphis marsupialis has been reported as a competent reservoir for trypanosomatid parasites infections. The aim of this study was to measure Trypanosoma cruzi, T. rangeli, and Leishmania spp. infection rates and to characterize discrete typing units (DTUs) of T. cruzi in D. marsupialis from two Chagas disease endemic sites in Panama. Blood from 57 wild-caught D. marsupialis were examined from two rural communities, Las Pavas (N = 18) and Trinidad de las Minas (N = 39). Twenty-two (38.60%) opossums were positive for flagellates by general hemoculture. T. cruzi infection was confirmed by positive hemoculture and/or kDNA based PCR performed in 31/57 (54.39%) blood samples from opossums. T. rangeli infection was confirmed by hemoculture and/or TrF/R2-Primer PCR assay applied on 12/57 (21.05%) blood samples. Nine (15.79%) D. marsupialis harbored T. cruzi/T. rangeli coinfections. All opossums tested negative for Leishmania spp. by PCR assays based on kDNA and HSP70 gene amplification. There was a significant association between T. cruzi infection and site (Fisher exact test, p = 0.02), with a higher proportion of T. cruzi infected opossums in Las Pavas (77.78%, n = 14/18) compared to Trinidad de las Minas (43.59%, n = 17/39). A significant association was found between habitat type and T. cruzi infection in opossums across both communities, (X2 = 6.91, p = 0.01, df = 1), with a higher proportion of T. cruzi infection in opossums captured in forest remnants (76%, 19/25) compared to peridomestic areas (37.5%, 12/32). T. rangeli detection, but not T. cruzi detection, may be improved by culture followed by PCR. TcI was the only DTU detected in 22 T. cruzi samples using conventional and real-time PCR. Eight T. rangeli positive samples were characterized as KP1(-)/lineage C. Trypanosome infection data from this common synanthropic mammal provides important information for improved surveillance and management of Chagas disease in endemic regions of Panama.

Plasmodium vivax Genetic Diversity in Panama: Challenges for Malaria Elimination in Mesoamerica.

Panama and all nations within the Mesoamerican region have committed to eliminate malaria within this decade. With more than 90% of the malaria cases in this region caused by Plasmodium vivax, an efficient national/regional elimination plan must include a comprehensive study of this parasite's genetic diversity. Here, we retrospectively analyzed P. vivax genetic diversity in autochthonous and imported field isolates collected in different endemic regions in Panama from 2007 to 2020, using highly polymorphic markers (csp, msp-1, and msp-3α). We did the analysis using molecular techniques that are cost-effective for malaria molecular surveillance within Mesoamerica. Thus, we used molecular analyses that are feasible for malaria molecular surveillance within the region, and that can provide useful information for policy and decision making about malaria elimination. We also evaluated if haplotypes established by combining the genotypes found in these genes were associated with relevant epidemiological variables and showed structure across the transmission foci that have been observed in Panama. Ten different haplotypes were identified, some of them strongly associated with geographical origin, age, and collection year. Phylogenetic analysis of csp (central repeat domain) revealed that both major variant types (vk210 and vk247) were circulating in Panama. Variant vk247 was restricted to the eastern endemic regions, while vk210 was predominant (77.3%) and widespread, displaying higher diversity (14 alleles) and geographically biased alleles. The regional implications of these molecular findings for the control of P. vivax malaria to achieve elimination across Mesoamerica are discussed.

Evaluation of cytochrome b sequence to identify Leishmania species and variants: the case of Panama.

BACKGROUND: The genetic heterogeneity of Leishmania parasites is a major factor responsible for the wide variety of Leishmania-associated manifestations. Consequently, understanding the genetic make-up of Leishmania species using suitable molecular markers is an important component of realising local and regional scale disease risk. The cytochrome b (cytb) is frequently used to type New World Leishmania species. However, its potential to discriminate Leishmania species and variants requires further evaluation. OBJECTIVES: To explore the capacity of cytb gene to identify New World Leishmania species and variants and to develop an approach able to type local Leishmania species and variants. METHODS: We retrieved 360 partial and complete Leishmania cytb gene sequences publicly available in GenBank database to study all single nucleotide polymorphisms (SNPs) across the cytb gene that differentiate New World Leishmania species. This information was used to develop an approach based upon the polymorphisms found in a DNA segment of 948bp. We also compared the typing results found with this technique with the polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) profiling obtained using HSP70 gene as target. One hundred Panamanian isolates were used to both typed Leishmania species and assess local genetic variability. FINDINGS: We found complete agreement between our cytb approach and the PCR-RFLP profiling method based on HSP70 for Leishmania species identification. Ninety-two isolates were identified as L. panamensis, although other Viannia species were found circulating at a lower frequency. Three L. panamensis haplotypes were identified in Panamanian provinces. We also provide an initial report of L. guyanensis haplotypes circulating in Panama. MAIN CONCLUSIONS: Cytb gene sequence encompasses key main SNPs that aid to identify Leishmania species. The cytb approach developed with this information was able to identify and assess genetic variability of local Leishmania species found in this study.

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-restriction fragment length polymorphism (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.

Evaluation of cytochrome b sequence to identify Leishmania species and variants: the case of Panama

BACKGROUND The genetic heterogeneity of Leishmania parasites is a major factor responsible for the wide variety of Leishmania-associated manifestations. Consequently, understanding the genetic make-up of Leishmania species using suitable molecular markers is an important component of realising local and regional scale disease risk. The cytochrome b (cytb) is frequently used to type New World Leishmania species. However, its potential to discriminate Leishmania species and variants requires further evaluation. OBJECTIVES To explore the capacity of cytb gene to identify New World Leishmania species and variants and to develop an approach able to type local Leishmania species and variants. METHODS We retrieved 360 partial and complete Leishmania cytb gene sequences publicly available in GenBank database to study all single nucleotide polymorphisms (SNPs) across the cytb gene that differentiate New World Leishmania species. This information was used to develop an approach based upon the polymorphisms found in a DNA segment of 948bp. We also compared the typing results found with this technique with the polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) profiling obtained using HSP70 gene as target. One hundred Panamanian isolates were used to both typed Leishmania species and assess local genetic variability. FINDINGS We found complete agreement between our cytb approach and the PCR-RFLP profiling method based on HSP70 for Leishmania species identification. Ninety-two isolates were identified as L. panamensis, although other Viannia species were found circulating at a lower frequency. Three L. panamensis haplotypes were identified in Panamanian provinces. We also provide an initial report of L. guyanensis haplotypes circulating in Panama. MAIN CONCLUSIONS Cytb gene sequence encompasses key main SNPs that aid to identify Leishmania species. The cytb approach developed with this information was able to identify and assess genetic variability of local Leishmania species found in this study.

Population genetics of two chromatic morphs of the Chagas disease vector Rhodnius pallescens Barber, 1932 in Panamá.

Rhodnius pallescens is the principal vector of Chagas disease in Panama. Recently a dark chromatic morph has been discovered in the highlands of Veraguas Province. Limited genetic studies have been conducted with regards to the population structure and dispersal potential of Triatominae vectors, particularly in R. pallescens. Next generation sequencing methods such as RADseq and complete mitochondrial DNA (mtDNA) genome sequencing have great potential for examining vector biology across space and time. Here we utilize a RADseq method (3RAD), along with complete mtDNA sequencing, to examine the population structure of the two chromatic morpho types of R. pallescens in Panama. We sequenced 105 R. pallescens samples from five localities in Panama. We generated a 2216 SNP dataset and 6 complete mtDNA genomes. RADseq showed significant differentiation among the five localities (FCT = 0.695; P = .004), but most of this was between localities with the dark vs. light chromatic morphs (Veraguas vs. Panama Oeste). The mtDNA genomes showed a 97-98% similarity between dark and light chromatic morphs across all genes and a 502 bp insert in light morphs. Thus, both the RADseq and mtDNA data showed highly differentiated clades with essentially no gene flow between the dark and light chromatic morphs from Veraguas and central Panama respectively. We discuss the growing evidence showing clear distinctions between these two morpho types with the possibility that these are separate species, an area of research that requires further investigation. Finally, we discuss the cost-effectiveness of 3RAD which is a third of the cost compared to other RADseq methods used recently in Chagas disease vector research.

Minimum free energy predicted base pairing in the 39 nt spliced leader and 5' UTR of calmodulin mRNA from Trypanosoma cruzi: influence of the multiple trans-splicing sites.

We analyzed the compositional changes and the stable base pairs in the predicted secondary structure of the 5' UTR calmodulin mRNA in T. cruzi. The three copies of calmodulin in T. cruzi genome display variable position of the trans splicing sites and give rise to several mRNA that differs slightly on 5' UTR composition in the epimastigote stage. We show that the pattern of high probability base pairs in the minimum free energy predicted secondary structures of the calmodulin 5' UTR remains unchanged despite the nucleotide composition variation. However, the 39 nt spliced leader (mini-exon, the 5' exon sequence transferred to trypanosome mRNAs by the mechanism of trans splicing) shows a variable pattern of high and low probability base pairing as consequence of the altered composition of the 5' UTR.

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.

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.

Cell disruption using a different methodology for proteomics analysis of Trypanosoma cruzi strains.

We have developed a cell disruption method to produce a protein extract using Trypanosoma cruzi cells based on a straightforward hypoosmotic lysis protocol. The procedure consists of three steps: incubation of the cells in a hypoosmotic lysis buffer, sonication in a water bath, and centrifugation. The final protein extract was designated TcS12. The stages of cell disruption at different incubation times were monitored by differential interference contrast microscopy. After 30min of incubation in lysis buffer at 4°C, the T. cruzi epimastigote forms changed from slender to round-shaped parasites. Nevertheless, cell disruption took place following sonication of the sample for 30min. The efficiency of the methodology was also validated by flow cytometry, which resulted in 72% of propidium iodide (PI)-labeled cells. To estimate the protein extraction yield and the differential protein expression, the proteomics profile of four T. cruzi strains (CL-Brener, Dm28c, Y, and 4167) were analyzed by liquid chromatography tandem mass spectrometry (LCMS/MS) on a SYNAPT HDMS system using the label-free MS(E) approach. ProteinLynx Global Server (version 2.5) with Expression(E) analysis identified a total of 1153 proteins and revealed 428 differentially expressed proteins among the strains. Gene ontology analysis showed that not only cytosolic proteins but also nuclear and organellar ones were present in the extract.

Evaluation of three different DNA extraction methods from blood samples collected in dried filter paper in Plasmodium subpatent infections from the Amazon region in Brazil.

Asymptomatic Plasmodium infection is a new challenge for public health in the American region. The polymerase chain reaction (PCR) is the best method for diagnosing subpatent parasitemias. In endemic areas, blood collection is hampered by geographical distances and deficient transport and storage conditions of the samples. Because DNA extraction from blood collected on filter paper is an efficient method for molecular studies in high parasitemic individuals, we investigated whether the technique could be an alternative for Plasmodium diagnosis among asymptomatic and pauciparasitemic subjects. In this report we compared three different methods (Chelex®-saponin, methanol and TRIS-EDTA) of DNA extraction from blood collected on filter paper from asymptomatic Plasmodium-infected individuals. Polymerase chain reaction assays for detection of Plasmodium species showed the best results when the Chelex®-saponin method was used. Even though the sensitivity of detection was approximately 66% and 31% for P. falciparum and P. vivax, respectively, this method did not show the effectiveness in DNA extraction required for molecular diagnosis of Plasmodium. The development of better methods for extracting DNA from blood collected on filter paper is important for the diagnosis of subpatent malarial infections in remote areas and would contribute to establishing the epidemiology of this form of infection.

Evaluation of three different DNA extraction methods from blood samples collected in dried filter paper in Plasmodium subpatent infections from the Amazon region in Brazil / Avaliação de três métodos de extração de DNA obtidos de amostras de sangue coletadas em papel de filtro em infecções subpatentes de Plasmodium da região Amazônica no Brasil

Asymptomatic Plasmodium infection is a new challenge for public health in the American region. The polymerase chain reaction (PCR) is the best method for diagnosing subpatent parasitemias. In endemic areas, blood collection is hampered by geographical distances and deficient transport and storage conditions of the samples. Because DNA extraction from blood collected on filter paper is an efficient method for molecular studies in high parasitemic individuals, we investigated whether the technique could be an alternative for Plasmodium diagnosis among asymptomatic and pauciparasitemic subjects. In this report we compared three different methods (Chelex®-saponin, methanol and TRIS-EDTA) of DNA extraction from blood collected on filter paper from asymptomatic Plasmodium-infected individuals. Polymerase chain reaction assays for detection of Plasmodium species showed the best results when the Chelex®-saponin method was used. Even though the sensitivity of detection was approximately 66% and 31% for P. falciparum and P. vivax, respectively, this method did not show the effectiveness in DNA extraction required for molecular diagnosis of Plasmodium. The development of better methods for extracting DNA from blood collected on filter paper is important for the diagnosis of subpatent malarial infections in remote areas and would contribute to establishing the epidemiology of this form of infection.

Infecção assintomática por Plasmodium é um novo desafio para a saúde pública no Brasil. A reação em cadeia da polimerase (PCR) é o melhor método para detectar baixas parasitemias presentes em pacientes com infecção assintomática. Nas áreas endêmicas, a coleta de sangue total é dificultada pela distancia geográfica, transporte e adequada armazenagem das amostras. A coleta de sangue em papel de filtro pode ser uma alternativa nessas áreas de difícil acesso. Neste estudo foram comparados três diferentes métodos de extração de ADN a partir de papel de filtro usando como controle extração a partir de sangue total. O protocolo Chelex®-Saponina foi o que obteve o melhor resultado quando comparado com os outros três protocolos. No entanto a sensibilidade foi de 66,7% para o P. falciparum e 31,6% para o P. vivax. Conclui-se que em caso de infecção assintomática o papel de filtro não é ainda uma boa alternativa para coleta de amostras.

Echinococcus granulosus genotypes circulating in alpacas (Lama pacos) and pigs (Sus scrofa) from an endemic region in Peru.

The identification of the genotypes of Echinococcus granulosus present in livestock and wild animals within regions endemic for cystic echinococcosis (CE) is epidemiologically important. Individual strains display different biological characteristics that contribute to outbreaks of CE and that must be taken into account in the design of intervention programs. In this study, samples of hydatid cysts due to E. granulosus were collected from alpacas (4) in Puno and pigs (8) in Ayacucho in Peru, an endemic region for CE. Polymerase chain reaction amplification and DNA sequencing of specific regions of the mitochondrial cytochrome C oxidase subunit 1 and NADH dehydrogenase subunit 1 genes confirmed the presence of a strain common to sheep, the G1 genotype, in alpacas. Two different strains of E. granulosus were identified in pigs: the G1 and the G7 genotypes. This is the first report of the G1 genotype of E. granulosus in alpacas in endemic regions of CE in Peru.

Echinococcus granulosus genotypes circulating in alpacas (Lama pacos) and pigs (Sus scrofa) from an endemic region in Peru

The identification of the genotypes of Echinococcus granulosus present in livestock and wild animals within regions endemic for cystic echinococcosis (CE) is epidemiologically important. Individual strains display different biological characteristics that contribute to outbreaks of CE and that must be taken into account in the design of intervention programs. In this study, samples of hydatid cysts due to E. granulosus were collected from alpacas (4) in Puno and pigs (8) in Ayacucho in Peru, an endemic region for CE. Polymerase chain reaction amplification and DNA sequencing of specific regions of the mitochondrial cytochrome C oxidase subunit 1 and NADH dehydrogenase subunit 1 genes confirmed the presence of a strain common to sheep, the G1 genotype, in alpacas. Two different strains of E. granulosus were identified in pigs: the G1 and the G7 genotypes. This is the first report of the G1 genotype of E. granulosus in alpacas in endemic regions of CE in Peru.

Rocky Mountain spotted fever in Panama: a cluster description.

Rocky Mountain spotted fever (RMSF) is a tick-borne infection caused by Rickettsia rickettsii. We report a cluster of fatal cases of RMSF in 2007 in Panama, involving a pregnant woman and two children from the same family.  The woman presented with a fever followed by respiratory distress, maculopapular rash, and an eschar at the site from which a tick had been removed.  She died four days after disease onset.  This is the second published report of an eschar in a patient confirmed by PCR to be infected with R. rickettsii.  One month later, the children presented within days of one another with fever and rash and died three and four days after disease onset. The diagnosis was confirmed by immunohistochemistry, PCR and sequencing of the genes of R. rickettsii in tissues obtained at autopsy. 

Human trypanosomiasis in the eastern region of the Panama Province: new endemic areas for Chagas disease.

The epidemiology of Chagas disease was studied in five rural communities located in the eastern region of the Panama Province. Serological tests for Trypanosoma cruzi infection revealed a prevalence of 5.88% (12/204). Hemocultures coupled with polymerase chain reaction (PCR) analysis showed a Trypanosoma rangeli infection rate of 5.88% (12/204). An overall trypanosome infection index of 11.76% (24/204) was detected in this population. A total of 121 triatomine specimens were collected in domestic and peridomestic habitats. Rhodnius pallescens was confirmed as the predominant species. Molecular analysis showed that 17.8% (13/73) of the examined insects were positive for T. cruzi, 17.8% (13/73) for T. rangeli, and 35.6% (26/73) presented mixed infections. Among 73 R. pallescens evaluated, 16.4% (12/73) contained opossum blood meals. The epidemiological implications of these findings are discussed.

Molecular epidemiology of American tegumentary leishmaniasis in Panama.

American tegumentary leishmaniasis is an increasing public health problem in Panama. This study describes the clinical characteristics and the molecular epidemiology of leishmaniasis in Panama over a 5-year period (2004-2008). Additionally, we applied a polymerase chain reaction (PCR)-restriction fragment length polymorphism (RFLP)-based assay to identify Leishmania species in clinical isolates, skin scrapings, and sandflies specimens. Whereas 60.3% of cases were detected with conventional parasitologic techniques (smear or in vitro culture), the PCR detected 72% positive patients. Our clinical-epidemiologic data corroborate the high incidence of L. (Viannia) panamensis and provide evidence of peridomestic and/or domestic transmission. Mucosal involvement was observed in 4.2% of the patients. The overall natural infection rate with Leishmania in 103 pools of sandflies was 0.46%. Lutzomyia gomezi and Lutzomya panamensis were the prevalent species incriminated as vectors at the capture sites in central Panama. This study contributes to a better knowledge of the current epidemiology of tegumentary leishmaniasis in Panama.