<i>Leishmania infantum</i> Genetic Diversity and<i> Lutzomyia longipalpis</i> Mitochondrial Haplotypes in Brazil

Ribolla, Paulo Eduardo Martins; Gushi, Letícia Tsieme; Pires e Cruz, Maria do Socorro; Costa, Carlos Henrique Nery; Costa, Dorcas Lamounier; Lima Júnior, Manoel Sebastião da Costa; Dorval, Maria Elizabeth Moraes Cavalheiros; Gutierrez de Oliveira, Alessandra; da Cunha Santos, Mirella Ferreira; Fonseca Camargo-Neves, Vera Lúcia; Fortaleza, Carlos Magno Castello Branco; Alonso, Diego Peres

doi:https://doi.org/10.1155/2016/9249217

BioMed Research International

On this page

Abstract Introduction Methods Results Discussion Conclusions Acknowledgments References Copyright Related Articles

Research Article | Open Access

Volume 2016 | Article ID 9249217 | https://doi.org/10.1155/2016/9249217

Leishmania infantum Genetic Diversity and Lutzomyia longipalpis Mitochondrial Haplotypes in Brazil

Paulo Eduardo Martins Ribolla,¹Letícia Tsieme Gushi,¹Maria do Socorro Pires e Cruz,²Carlos Henrique Nery Costa,³Dorcas Lamounier Costa,³Manoel Sebastião da Costa Lima Júnior,⁴Maria Elizabeth Moraes Cavalheiros Dorval,⁴Alessandra Gutierrez de Oliveira,⁵Mirella Ferreira da Cunha Santos,⁵Vera Lúcia Fonseca Camargo-Neves,⁶Carlos Magno Castello Branco Fortaleza,⁷and Diego Peres Alonso¹

Academic Editor: Amogh A. Sahasrabuddhe

Received06 Nov 2015

Revised22 Feb 2016

Accepted07 Mar 2016

Published29 Mar 2016

Abstract

Leishmania infantum is the etiological agent of visceral leishmaniasis (VL) in the Americas with domestic dogs being its major reservoir hosts. The main VL vector is the sandfly Lutzomyia longipalpis, while other Lutzomyia species may play a role in disease transmission. Although the genetic structure of L. infantum populations has been widely evaluated, only a few studies have addressed this subject coupled to the genetic structure of the respective sandfly vectors. In this study, we analyzed the population structure of L. infantum in three major VL endemic areas in Brazil and associated it with Lutzomyia longipalpis geographic structure.

1. Introduction

Leishmaniases are parasitic diseases caused by protozoans from the genus Leishmania, which are transmitted by the bite of female sandflies from the family Psychodidae. The clinical manifestations of leishmaniases are particularly diverse and present different characteristics: visceral leishmaniasis (VL), the most severe one; mucocutaneous leishmaniasis, characterized as a mutilating disease; diffuse cutaneous leishmaniasis, caused by a deficient cellular immune response; and cutaneous leishmaniasis, which causes single or multiple lesions on the skin. The epidemiology of leishmaniasis is highly complex: there are 20 known species of Leishmania pathogenic to humans and at least 30 species of sandflies vectors. Furthermore, this disease can be designated as a zoonosis, which involves animals as the reservoir hosts or as an anthroponosis, when humans are the only source of parasites for sandflies. Leishmaniasis is widely spread in 98 countries and 3 territories, from which more than 70% are developing countries and 13 are among the least developed ones [1].

Visceral leishmaniasis can be either an anthroponosis (e.g., in the Indian subcontinent) or a zoonosis (e.g., in the Mediterranean or in the Americas), and it is characterized by chronic evolution and systemic involvement, which if untreated may result in death. In the Americas, Leishmania infantum is the etiological agent of the disease and Brazil accounts for over 90% of the cases in the continent [1, 2]. Domestic dogs are the proven reservoir hosts in rural and urban areas, while the role of naturally infected wild mammals (canids and marsupials) as L. infantum reservoir hosts is still controversial [3]. The main sandfly vector is Lutzomyia longipalpis, but other Lutzomyia species might play a role in disease transmission; for example, in Corumbá, Mato Grosso do Sul, naturally Leishmania-infected Lu. cruzi have been discovered and because there is still no evidence of Lu. longipalpis in this region, that sandfly is considered the main vector [4, 5].

In Brazil, VL typically occurred in rural settings, but since 1980 its incidence has been changing due to widespread urban outbreaks. The first major VL urban epidemic in the country happened in Teresina, Piauí State. Since then, epidemics occurred in Natal (Rio Grande do Norte) and São Luís (Maranhão), and the disease subsequently spread to other regions of the country. Autochthonous cases were recently described for the first time in the southernmost State of Rio Grande do Sul. The current epidemiological scenario of VL leaves no doubt regarding the severity of the situation and the unchecked geographic spread of the disease. In the 1990s, only 10% of the cases occurred outside the Northeast Region, but in 2007 the proportion reached 50% of cases. From 2006 to 2008, autochthonous transmission of VL was reported in more than 1,200 municipalities in 21 states [6].

The broad spectrum of leishmaniasis-associated symptoms, coupled with the wide diversity of vertebrate and invertebrate host species, suggests that both parasites’ and hosts’ genetic backgrounds determine the patterns of the disease [7]. On the other hand, clonal diversity and genetic heterogeneity, which can cause variability in parasite virulence, are quite common in Leishmania [8].

Several studies showed that genetic variability of L. infantum in Brazil is low, with restricted diversity and limited population clustering. In a recent work assessing parasite populations distributed over 18 states, three major clustered populations could be inferred using microsatellite typing. When the analysis is performed in parasites from closely related geographic regions, the overall diversity is even lower [9–11].

When we look at sandfly genetic variability, there is compelling evidence that the Lutzomyia population structure in Brazil is complex, with different genotypes identified depending on the geographic region assessed and also the species involved in parasite transmission [12–14].

Based on these studies, it is logical to hypothesize that the interactions of L. infantum genotypic variants with different hosts and vector populations may ultimately influence the transmission dynamics and severity of eventual outbreaks. Hence, assessing the genetic structure of both vector populations and parasites may help us to understand the dynamics of vector-parasite interactions and the epidemiological aspects of American visceral leishmaniasis. Here, we used PCR-RFLP of kinetoplast minicircle DNA (kDNA) to identify L. infantum genotypic variants from three VL endemic areas in Brazil: Teresina in Piauí State, Campo Grande in Mato Grosso do Sul State, and Bauru in São Paulo State. kDNA-RFLP analysis when compared to microsatellite genotyping has proven to be more sensitive to examine genetic data of closely related sympatric L. infantum strains [15]. In addition, in order to identify different haplotypes of Lu. longipalpis and Lu. cruzi sandflies from those three VL endemic areas, we used mitochondrial 12S rDNA sequencing. As a maternal inheritance, rapidly evolving, nonrecombinant and haploid molecular marker, 12S rDNA is suitable to trace genealogy and evolutionary history. To our knowledge, this is the first study that seeks to compare genetic variability of Leishmania infantum parasites to the genetic structure of its vectors in Brazil.

2. Methods

2.1. Ethics Statement

For insect collections in Mato Grosso do Sul State, we obtained a permanent license for collecting and transporting zoological material N° 25592-1 on behalf of Dr. Alessandra Gutierrez de Oliveira, issued by the System of Authorization and Information on Biodiversity of the Brazilian Institute of Environment and Renewable Natural Resources (Sisbio/IBAMA). For insect collections in São Paulo State and Piauí State, no specific permissions were required since the specimens were kindly provided by the Center for the Control of Endemic Diseases (SUCEN) and Federal Piauí State University, respectively. The collections were performed at private residences, whose owners personally granted permission to enter their backyards to capture the sandflies. All of these residences were located in urban areas and no endangered or protected species were collected in this study.

2.2. Sandfly Collections

Sandflies were captured by both manual collection and electric traps. Manual collection was performed with electric aspirators, restricting the use of a Castro catcher to locations where aspiration could not be used. The selected collection points were preassessed in order to establish the best capturing location in the peridomicile. At each selected point, modified CDC light traps were installed from 6 p.m. to 6 a.m.

The collections took place in different areas in Brazil and were performed by the respective local teams: São Paulo (SP) State, performed by Center for the Control of Endemic Diseases (SUCEN); Piauí (PI) State, performed by Piauí Federal University; and Mato Grosso do Sul (MS) State, carried out by Mato Grosso do Sul Federal University.

Lu. cruzi was collected in Corumbá (MS) and Lu. longipalpis in all other places: Campo Grande (MS), Teresina (PI), Andradina (SP), Araçatuba (SP), and Birigui (SP). All identified insects were kept in 70% ethanol until use.

2.3. Sandfly Genomic DNA Isolation

The field-derived sandflies were grinded with the help of a plastic pestle in 1.5 mL tubes containing 300 μL of 5% Chelex® (Bio-Rad). The solution was then vortexed for 15 s, centrifuged at 11,000 g for 20 s, and incubated at 80°C for 30 min, after which the procedure was repeated. The supernatant was finally removed, transferred to another 1.5 mL microcentrifuge tube, and stored at −20°C. We had an average of 45 ng of DNA per sandfly measured with NanoDrop™ 1000 (Thermo Scientific).

2.4. Sandfly Mitochondrial 12S rDNA Amplification and Sequencing

PCR amplification of the Lutzomyia sp. 12S rDNA mitochondrial region was performed with the primers T1B (5′-AAACTAGGATTAGATACCT-3′) and T2B (5′-AATGAGAGCGACGGGCGATG-3′), according to Beati et al. [16]. Reactions of 25 μL were set up as follows: 13.7 μL of ultrapure water, 2.5 μL of 10x Platinum buffer (Life Technologies), 1.0 μL MgCl2 (50 mM), 0.5 μL dNTPs (0.1 mM), 1.0 μL of each oligonucleotide (10 pmol/μL), 0.3 μL of Platinum Taq, (Life Technologies; 5 U/μL), and 10 ng of genomic DNA. The reaction was carried out in a thermal cycler as follows: 5 cycles of 94°C for 15 s, 51°C for 30 s, and 68°C for 30 s, followed by 25 cycles of 94°C for 15 s, 53°C for 30 s, and 70°C for 30 s, and a final extension step of 70°C for 5 min. The amplified DNA fragments were UV visualized after electrophoresis on 1% agarose gel stained with ethidium bromide.

The resulting DNA fragments were purified with ExoSAP-IT kit (GE Healthcare), according to the manufacturer’s protocol. The 20 μL sequencing reactions consisted of 2 μL of BigDye Terminator (Life Technologies), 6.0 μL of BigDye Terminator 5x Sequencing Buffer (Life Technologies), 3.2 μL of the primers (1 pmol/μL), 4.8 μL of ultrapure water, and 200 ng of DNA measured with NanoDrop 1000 (Thermo Scientific). All reactions were carried out in a thermal cycler, with 35 cycles of 95°C for 20 s, 50°C for 15 s, and 60°C for 2 min. The amplified DNA was precipitated with 80 μL of 65% isopropanol, washed with 200 μL of 70% ethanol, and air-dried for 5 min. Before injection, samples were resuspended in 10 μL of HI-DI formamide (Life Technologies) and heated at 95°C for 3 min for DNA denaturation and immediately cooled on ice. Sample processing occurred in an ABI377 automatic sequencer.

2.5. Sequencing Analysis

The forward and reverse 12S rDNA sequences were manually checked for quality and the polymorphisms confirmed and then matched using the online EMBOSS GUI tool package (http://imed.med.ucm.es/cgi-bin/emboss.pl?_action=input&_app=merger). The obtained consensus sequences were aligned using Clustal X2 software. Polymorphisms in each sequence were identified and a haplotypic diversity test (Table 1) was performed with the DnaSP 5.10 software. Haplotype diagram generation was performed by TCS: phylogenetic network using statistical estimation parsimony software.

2.6. Parasite Samples and DNA Isolation

Parasites used in this study were collected between 2007 and 2009 (Table 2). The DNA from the promastigotes (from all cultured samples used and for the two parasite samples obtained from sandflies) was isolated with Chelex (Bio-Rad). Briefly, 1 mL aliquots of the cultures were transferred to 1.5 mL centrifuge tubes and spun down for 1 min at 10,000 g. The supernatant was discarded and the pellet resuspended in 300 μL of 10% Chelex (w/v). Following, the samples were incubated for 15 min at 95°C and then centrifuged again for 1 min at 10,000 g. The supernatant containing the DNA was then carefully recovered and stored in a new tube at −20°C. For the two parasite samples isolated from sandflies, the whole insect was grinded in 300 μL of 10% Chelex with the help of a motorized tissue grinder, following the same steps above. We had an average of 200 ng of DNA per culture sampled and 20 ng per sample for the two sandfly-derived parasites measured with NanoDrop 1000 (Thermo Scientific).

The DNA of L. infantum amastigotes was extracted following two different approaches. For dog bone marrow aspirates we used the Illustra Blood GenomicPrep Mini Spin kit (GE Healthcare) according to the manufacturer’s recommendations. For slide-fixed human bone marrow aspirates we used the same protocol after scraping the contents of each slide into a 1.5 mL tube, as previously described [17]. We had an average of 100 ng per dog bone marrow sample and 25 ng of DNA per slide measured with NanoDrop 1000 (Thermo Scientific).

2.7. PCR-RFLP of Kinetoplast DNA (kDNA) and RFLP Analysis

We had initially started our analysis using a panel of 7 microsatellite markers (Li22-35, Li23-41, Li45-24, Li71-33, Lm2TG, Lm4TA, and TubCA) [10]. However, only one marker (Li45-24) was polymorphic and, due to its low variability, only two alleles could be identified. For this reason, we decided to perform only PCR-RFLP of kinetoplast DNA (kDNA) and RFLP analysis.

For the analysis of the kinetoplast minicircle DNA, 157 L. infantum isolates were used (Table 2): 98 cultured samples initially isolated from human patients by sternal bone marrow aspiration (44 from Teresina and 54 from Campo Grande), 42 samples from dog bone marrow aspirates from Teresina, 2 samples from sandflies blood-fed on L. infantum-infected dogs from this same study in Teresina, and 15 slide-derived samples originated from bone marrow aspirates of human patients in Bauru, São Paulo State. PCR reactions were performed with primers LINR4 and LIN19 [18] and generated a 720 bp amplicon, which covers almost the entire minicircle. The 50 μL reactions contained 1 mM MgCl2, 10 mM Tris-HCl (pH 8.3), 0.3 pmol of each oligonucleotide, 0.1 mM dNTPs, 1 unit of Taq polymerase (GE Healthcare), and 5 μL of sample DNA. The amplification conditions were as follows: 3 min at 94°C, 33 cycles of 30 s at 95°C, 30 s at 58°C, and 1 min at 72°C, followed by a final extension step of 10 minutes at 72°C. The PCR products were then precipitated with ethanol, resuspended in water and digested with the restriction enzymes RsaI and HpaII (Promega) as previously described [15]. Approximately 1 μg of each PCR product was used per digestion in order to ensure that all reactions had the same initial amount of DNA. Since the products smaller than 100 bp can be confused with primer dimers and the ones larger than 700 bp can be misidentified as undigested products, only the fragments within this range were used in our RFLP analysis.

Data analysis was performed using R software environment. A binary matrix was constructed based on the profile of fragments generated by each digestion, where 1 represents the presence of a fragment and 0 represents its absence. This matrix was converted into a similarity matrix using the package “proxy” and used for cluster analysis. After, -means partitioning method was used to infer the number of clusters using the package “-means” and Agglomerative Hierarchical Clustering dendrogram was built using the binary distance method and ward cluster method with the package “hclust”.

3. Results

3.1. Sandflies Genetic Analysis

DNA was extracted from a total of 140 individuals as follows: 30 individuals from Andradina (SP), Araçatuba (SP), and Birigui (SP); 29 individuals from Teresina (PI); 14 individuals from Campo Grande (MS); 7 individuals from Corumbá (MS), classified morphologically as Lu. cruzi (Figure 1). PCR reactions generated a mitochondrial 12S ribosomal DNA fragment of approximately 360 bp, as previously described [19], which was then partially sequenced (263 bp). Sequences were screened for significant polymorphisms, and 10 variable sites were found (Table 3). When polymorphisms were assessed with DnaSP 5.10 program, 13 haplotypes were generated: six haplotypes (H8, H9, H10, H11, H12, and H13) containing only individuals from Teresina (PI); five haplotypes (H3, H4, H5, H6, and H7) containing only Araçatuba individuals (SP); one haplotype (H1) containing one individual from Corumbá and one individual from Campo Grande (MS); and one haplotype (H2) covering most of the sequences (111 individuals). Data are represented in a diagram of haplotypes (Figure 2). The haplotypic diversity test showed that Teresina presented the highest diversity (0.672), followed by Araçatuba (0.545), Corumbá (0.286), and Campo Grande (0.143). Andradina and Birigui presented no haplotypic diversity at all (Table 1).

3.2. Parasites RFLP Analysis

The kDNA fragments of interest were successfully amplified from the LinR4 and Lin19 oligos used in this study. RFLP analysis of kinetoplast minicircles DNA was also efficient in detecting restriction patterns between different samples. From the 157 tested samples, we could observe 55 unique genotypes in the cluster analysis dendrogram illustrated in Figure 3. -means partitioning identified 6 major clusters; there was a clear distinction between samples from Teresina, which grouped in two almost exclusive clusters, and all other samples; an exclusive Bauru cluster was also found. Two clusters presented with Teresina and Campo Grande samples, and one cluster presented with Bauru and Campo Grande samples. It is noteworthy that Campo Grande is distributed over 3 major clusters, one that groups together with one Teresina major cluster and the other two that are closer to Bauru clusters in a separate branch of the dendrogram. There was no clustering differentiation related to the years of collection.

4. Discussion

During the past 20 years, the epidemiology of VL has been constantly changing due to a continuous urbanization process, an increasing incidence of HIV/Leishmania coinfections, and syringe sharing by intravenous drug users [20] and the identification of novel L. infantum mammalian hosts/reservoirs [21]. This highlights the necessity of molecularly tracking the geographic distribution of different parasite and vector populations in order to enhance the knowledge on basic epidemiological aspects of the disease, such as its natural history and transmission.

Several molecular approaches have been used in the characterization of genetic variants in the genus Leishmania: amplified polymorphic DNA (RAPD) markers [22], analysis by size polymorphism of restriction fragments (RFLP) of the ITS regions ribosomal DNA [23], and kinetoplast DNA [24]; analysis confirmed sequence amplified regions [25]; and analysis of regions of DNA with microsatellite markers [19, 26–29]. We then decided to proceed with PCR-RFLP analysis of minicircle kDNA because it has a higher resolving power when applied to population genetics studies involving either genetically or geographically closely related strains [24, 30, 31]. Our data revealed a clear distinction between samples from Teresina, which grouped in two almost exclusive clusters, and all other samples; an exclusive Bauru cluster was also found. Two clusters presented with Teresina and Campo Grande samples, and one cluster presented with Bauru and Campo Grande samples. These results allowed us to draw a relationship between genetic distance and geographic origin. Interestingly, geographic origin related to diverse genetic background was also found for L. infantum parasites in Brazil in the study performed by Segatto et al. [10].

Our data is partially in accordance with a previous microsatellite based genotyping study performed with parasite populations from all 5 Brazilian regions. In the study, three well-defined populations could be identified; one that was present mostly in Northeast region, (including Piauí State that was sampled in our study) and the other two present in Midwest region (including Mato Grosso and Mato Grosso do Sul States that were sampled in our study). On the other hand, parasites typed in Southeast region (including São Paulo State that was sampled in our study) are closely related to northeastern parasites while in our study they are closely related to Midwestern parasites [9]. Our findings corroborate the use of this technique in Leishmania genotyping studies and reinforce the idea that in some cases, especially when analyzing strains of very close geographical origin, it is the only molecular marker capable of producing detectable patterns of polymorphism [24, 32].

All these genotyping studies on L. infantum suggest that the nuclear genomic variability of this species is likely to be low. Our hypothesis is that the kinetoplast genome can serve as a source of genetic variability for these parasites. The kDNA minicircles are essential for the function of the trypanosomatid’s mitochondrial genes, as minicircles code for guide RNAs, which play an essential role in editing messenger RNA (mRNA) from the maxicircles that contain genes for essential mitochondrial proteins [33]. Therefore, this DNA is more prone to a rapid response to diverse ambient conditions and stress situations, and parasite fitness conferring different selective advantages might depend on which minicircle classes prevail in different Leishmania strains.

A similar phenomenon, known as transkinetoplastidy, has been described in Leishmania and is responsible for changes in minicircles classes when the parasites are challenged with increasing concentrations of drugs that are normally lethal. This will in turn cause a dramatic change in the abundance of certain minicircles classes, which during transkinetoplastidy will be increased or reduced and replaced by a previously less frequent class [34].

When we look at sandfly genetic analysis we can clearly observe a main haplotype (H2) comprising all individuals from Andradina and Birigui, 13 out of 14 individuals from Campo Grande, 6 out of 7 individuals from Corumbá, 20 out of 30 individuals from Araçatuba, and 11 out 30 individuals from Teresina. There is also a major haplotype (H8) comprising only individuals from Teresina (13 out of 29) and minor haplotypes from Araçatuba. From the 12S rDNA sequencing data, it was not possible to differentiate Lu. longipalpis from Lu. cruzi (Corumbá) since there was no haplotype clustering among Corumbá sandflies. This may suggest that the process of speciation is recent or still occurring. A microsatellite based study assessing the genetic variability of Lu. longipalpis and Lu. cruzi populations in Mato Grosso do Sul State showed evidence of introgression and limited gene flow between the two species, corroborating our findings [12].

In general, we can summarize the data obtained from haplotyping as follows: a major haplotype composed of 111 individuals (comprising 89% of SP, 90% of MS, and 38% of PI individuals); a main haplotype composed of 13 individuals exclusively from Teresina and giving rise to other 4 Teresina exclusive haplotypes (62% of individuals from Teresina with exclusive haplotypes); minor haplotypes comprising only individuals from SP (11% total) and from the same locality (Araçatuba).

When we compare data from parasite genotyping with sandfly 12S rDNA sequencing, the correlation of the two datasets is remarkable. Both show most samples from PI clearly separated from the MS and SP ones which are in turn much more related to each other when compared to PI that presented the highest haplotypic diversity (Table 1). The exception comes from the minor vector haplotypes only found in Araçatuba samples. Araçatuba represents an important landmark in the natural history of VL in SP given the fact that the first VL outbreak registered in the state occurred in this location [35, 36]. This could be a possible explanation to its greater number of unique haplotypes as one can assume that coevolution between parasites and vectors happens for a longer time in this area; this is supported by the high haplotypic diversity found for this population (Table 1). Taken together, these data corroborate that the sandfly vector probably plays an important role in shaping the genetic structure of L. infantum in Brazil as described by Ferreira et al. [9].

This work presents new insights towards the understanding of the population structure of L. infantum and Lu. longipalpis from VL endemic areas in Brazil. Further analyses will be needed to elucidate how different vector populations shape the genetic variability of L. infantum.

5. Conclusions

Taken together, our data indicate that the sandfly vector might play a role in selecting specific parasite strains at a regional level and therefore contributing to the genetic structure of L. infantum in Brazil. Assessing the genetic structure of both vector and parasite populations may help us to understand the evolution process surrounding vector-parasite interactions and shed light on a fundamental aspect of the ecoepidemiology of American visceral leishmaniasis.

Competing Interests

The authors declare that they have no competing interests.

Acknowledgments

The authors would like to acknowledge the Fundação de Amparo à Pesquisa do Estado de São Paulo (Research Grant 2009/10030-9 to Paulo Eduardo Martins Ribolla and Ph.D. Fellowship 2006/61151-2 to Diego Peres Alonso).

References

P. Desjeux, “Leishmaniasis: current situation and new perspectives,” Comparative Immunology, Microbiology and Infectious Diseases, vol. 27, no. 5, pp. 305–318, 2004.
View at: Publisher Site | Google Scholar
J. Alvar, I. D. Vélez, C. Bern et al., “Leishmaniasis worldwide and global estimates of its incidence,” PLoS ONE, vol. 7, no. 5, Article ID e35671, 2012.
View at: Publisher Site | Google Scholar
M. E. C. Dorval, E. T. Oshiro, E. Cupollilo, A. C. C. Castro, and T. P. Alves, “Ocorrência de leishmaniose tegumentar americana no Estado do Mato Grosso do Sul associada à infecção por Leishmania (Leishmania) amazonensis,” Revista da Sociedade Brasileira de Medicina Tropical, vol. 39, no. 1, pp. 43–46, 2006.
View at: Publisher Site | Google Scholar
O. Mangabeira Filho, “Sobre a sistemática e biologia dos Phlebotomus do Ceará,” Revista Brasileira de Malariologia e Doenças Tropicais, vol. 21, pp. 3–26, 1969.
View at: Google Scholar
R. D. Ward, A. L. Ribeiro, P. D. Ready, and A. Murtagh, “Reproductive isolation between different forms of Lutzomyia longipalpis (Lutz & Neiva), (Diptera: Psychodidae), the vector of Leishmania donovani chagasi Cunha & Chagas and its significance to kala-azar distribution in South America,” Memórias do Instituto Oswaldo Cruz, vol. 78, no. 3, pp. 269–280, 1983.
View at: Publisher Site | Google Scholar
G. L. Werneck, “Geographic spread of visceral leishmaniasis in Brazil,” Cadernos de Saúde Pública, vol. 26, no. 4, pp. 644–645, 2010.
View at: Publisher Site | Google Scholar
M. E. J. Woolhouse, L. H. Taylor, and D. T. Haydon, “Population biology of multihost pathogens,” Science, vol. 292, no. 5519, pp. 1109–1112, 2001.
View at: Publisher Site | Google Scholar
Y. J.-F. Garin, A. Sulahian, F. Pratlong et al., “Virulence of Leishmania infantum is expressed as a clonal and dominant phenotype in experimental infections,” Infection and Immunity, vol. 69, no. 12, pp. 7365–7373, 2001.
View at: Publisher Site | Google Scholar
G. E. M. Ferreira, B. N. dos Santos, M. E. C. Dorval et al., “The genetic structure of leishmania infantum populations in Brazil and its possible association with the transmission cycle of visceral leishmaniasis,” PLoS ONE, vol. 7, no. 5, Article ID e36242, 2012.
View at: Publisher Site | Google Scholar
M. Segatto, L. S. Ribeiro, D. L. Costa et al., “Genetic diversity of Leishmania infantum field populations from Brazil,” Memórias do Instituto Oswaldo Cruz, vol. 107, no. 1, pp. 39–47, 2012.
View at: Publisher Site | Google Scholar
L. F. D. S. Batista, M. Segatto, C. E. S. Guedes et al., “An assessment of the genetic diversity of Leishmania infantum isolates from infected dogs in Brazil,” American Journal of Tropical Medicine and Hygiene, vol. 86, no. 5, pp. 799–806, 2012.
View at: Publisher Site | Google Scholar
M. F. C. Santos, P. E. M. Ribolla, D. P. Alonso et al., “Genetic structure of lutzomyia longipalpis populations in Mato Grosso do Sul, Brazil, based on microsatellite markers,” PLoS ONE, vol. 8, no. 9, Article ID e74268, 2013.
View at: Publisher Site | Google Scholar
A. S. Araki, F. M. Vigoder, L. G. Bauzer et al., “Molecular and behavioral differentiation among Brazilian populations of Lutzomyia longipalpis (Diptera: Psychodidae: Phlebotominae),” PLoS Neglected Tropical Diseases, vol. 3, no. 1, article e365, 2009.
View at: Publisher Site | Google Scholar
V. de Queiroz Balbino, I. V. Coutinho-Abreu, I. V. Sonoda et al., “Genetic structure of natural populations of the sand fly Lutzomyia longipalpis (Diptera: Psychodidae) from the Brazilian northeastern region,” Acta Tropica, vol. 98, no. 1, pp. 15–24, 2006.
View at: Publisher Site | Google Scholar
D. P. Alonso, D. L. Costa, I. L. De Mendonça, C. H. N. Costa, and P. E. M. Ribolla, “Short report: heterogeneity of Leishmania infantum chagasi kinetoplast DNA in Teresina (Brazil),” American Journal of Tropical Medicine and Hygiene, vol. 82, no. 5, pp. 819–821, 2010.
View at: Publisher Site | Google Scholar
L. Beati, A. G. Cáceres, J. A. Lee, and L. E. Munstermann, “Systematic relationships among Lutzomyia sand flies (Diptera: Psychodidae) of Peru and Colombia based on the analysis of 12S and 28S ribosomal DNA sequences,” International Journal for Parasitology, vol. 34, no. 2, pp. 225–234, 2004.
View at: Publisher Site | Google Scholar
H. Motazedian, M. Karamian, H. A. Noyes, and S. Ardehali, “DNA extraction and amplification of Leishmania from archived, Giemsa-stained slides, for the diagnosis of cutaneous leishmaniasis by PCR,” Annals of Tropical Medicine and Parasitology, vol. 96, no. 1, pp. 31–34, 2002.
View at: Publisher Site | Google Scholar
A. M. Aransay, E. Scoulica, and Y. Tselentis, “Detection and identification of Leishmania DNA within naturally infected sand flies by seminested PCR on minicircle kinetoplastic DNA,” Applied and Environmental Microbiology, vol. 66, no. 5, pp. 1933–1938, 2000.
View at: Publisher Site | Google Scholar
N. Kebede, S. Oghumu, A. Worku, A. Hailu, S. Varikuti, and A. R. Satoskar, “Multilocus microsatellite signature and identification of specific molecular markers for Leishmania aethiopica,” Parasites and Vectors, vol. 6, article 160, 2013.
View at: Publisher Site | Google Scholar
J. Alvar, P. Aparicio, A. Aseffa et al., “The relationship between leishmaniasis and AIDS: the second 10 years,” Clinical Microbiology Reviews, vol. 21, no. 2, pp. 334–359, 2008.
View at: Publisher Site | Google Scholar
M. Gramiccia and L. Gradoni, “The current status of zoonotic leishmaniases and approaches to disease control,” International Journal for Parasitology, vol. 35, no. 11-12, pp. 1169–1180, 2005.
View at: Publisher Site | Google Scholar
A. Toledo, J. Martín-Sánchez, B. Pesson, C. Sanchiz-Marín, and F. Morillas-Márquez, “Genetic variability within the species Leishmania infantum by RAPD. A lack of correlation with zymodeme structure,” Molecular and Biochemical Parasitology, vol. 119, no. 2, pp. 257–264, 2002.
View at: Publisher Site | Google Scholar
E. Cupolillo, L. R. Brahim, C. B. Toaldo et al., “Genetic polymorphism and molecular epidemiology of Leishmania (Viannia) braziliensis from different hosts and geographic areas in Brazil,” Journal of Clinical Microbiology, vol. 41, no. 7, pp. 3126–3132, 2003.
View at: Publisher Site | Google Scholar
T. Laurent, S. Rijal, V. Yardley et al., “Epidemiological dynamics of antimonial resistance in Leishmania donovani: genotyping reveals a polyclonal population structure among naturally-resistant clinical isolates from Nepal,” Infection, Genetics and Evolution, vol. 7, no. 2, pp. 206–212, 2007.
View at: Publisher Site | Google Scholar
J.-P. Gangneux, J. Menotti, F. Lorenzo et al., “Prospective value of PCR amplification and sequencing for diagnosis and typing of Old World Leishmania infections in an area of nonendemicity,” Journal of Clinical Microbiology, vol. 41, no. 4, pp. 1419–1422, 2003.
View at: Publisher Site | Google Scholar
M. B. Jamjoom, R. W. Ashford, P. A. Bates, S. J. Kemp, and H. A. Noyes, “Towards a standard battery of microsatellite markers for the analysis of the Leishmania donovani complex,” Annals of Tropical Medicine and Parasitology, vol. 96, no. 3, pp. 265–270, 2002.
View at: Publisher Site | Google Scholar
B. Bulle, L. Millon, J.-M. Bart et al., “Practical approach for typing strains of Leishmania infantum by microsatellite analysis,” Journal of Clinical Microbiology, vol. 40, no. 9, pp. 3391–3397, 2002.
View at: Publisher Site | Google Scholar
S. Ochsenreither, K. Kuhls, M. Schaar, W. Presber, and G. Schönian, “Multilocus microsatellite typing as a new tool for discrimination of Leishmania infantum MON-1 strains,” Journal of Clinical Microbiology, vol. 44, no. 2, pp. 495–503, 2006.
View at: Publisher Site | Google Scholar
L. Montoya, M. Gállego, B. Gavignet et al., “Application of microsatellite genotyping to the study of a restricted Leishmania infantum focus: different genotype compositions in isolates from dogs and sand flies,” The American Journal of Tropical Medicine and Hygiene, vol. 76, no. 5, pp. 888–895, 2007.
View at: Google Scholar
A. Nasereddin, K. Azmi, C. L. Jaffe et al., “Kinetoplast DNA heterogeneity among Leishmania infantum strains in central Israel and Palestine,” Veterinary Parasitology, vol. 161, no. 1-2, pp. 126–130, 2009.
View at: Publisher Site | Google Scholar
Y. Botilde, T. Laurent, W. Q. Tintaya et al., “Comparison of molecular markers for strain typing of Leishmania infantum,” Infection, Genetics and Evolution, vol. 6, no. 6, pp. 440–446, 2006.
View at: Publisher Site | Google Scholar
G. Van der Auwera, N. R. Bhattarai, S. Odiwuor, and M. Vuylsteke, “Remarks on identification of amplified fragment length polymorphisms linked to SAG resistance in Leishmania,” Acta Tropica, vol. 113, no. 1, pp. 92–93, 2010.
View at: Publisher Site | Google Scholar
B. Liu, Y. Liu, S. A. Motyka, E. E. C. Agbo, and P. T. Englund, “Fellowship of the rings: the replication of kinetoplast DNA,” Trends in Parasitology, vol. 21, no. 8, pp. 363–369, 2005.
View at: Publisher Site | Google Scholar
T. A. Shapiro and P. T. Englund, “The structure and replication of kinetoplast DNA,” Annual Review of Microbiology, vol. 49, pp. 117–143, 1995.
View at: Publisher Site | Google Scholar
V. L. F. de Camargo-Neves and G. Katz, “Leishmaniose visceral americana no Estado de São Paulo,” Revista da Sociedade Brasileira de Medicina Tropical, vol. 32, supplement 2, pp. 63–64, 1999.
View at: Google Scholar
M. Z. Galimberti, G. Katz, V. L. F. de Camargo-Neves et al., “Leishmaniose visceral americana no Estado de São Paulo,” Revista da Sociedade Brasileira de Medicina Tropical, vol. 32, supplement 1, pp. 217–218, 1999.
View at: Google Scholar

Copyright

Copyright © 2016 Paulo Eduardo Martins Ribolla et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

1894

Downloads

1012

Citations