Previous Article | Next Article ![]()
Journal of Clinical Microbiology, July 2003, p. 3126-3132, Vol. 41, No. 7
0095-1137/03/$08.00+0 DOI: 10.1128/JCM.41.7.3126-3132.2003
Copyright © 2003, American Society for Microbiology. All Rights Reserved.
Laboratório de Pesquisas em Leishmaniose-Instituto Oswaldo Cruz,1 Hospital Evandro Chagas,Fundação Oswaldo Cruz Rio de JaneiroBrazil,2 Centro de Pesquisas Ageu Magalhães, Fundação Oswaldo CruzRecifeBrazil,3 Universidade Federal do Espírito SantoVitoriaBrazil,4 Instituto Nacional de Pesquisas da AmazôniaManausBrazil5
Received 18 November 2002/ Returned for modification 31 January 2003/ Accepted 31 March 2003
|
|
|---|
|
|
|---|
Most of the environmental factors affecting the epidemiology of the various leishmaniases are still poorly understood. Wild mammals serve as reservoirs for most of the New World Leishmania spp. (21), but there is increasing evidence that some of the human pathogenic Leishmania strains can be maintained in both sylvan and urban cycles. In the case of L. (V.) braziliensis, the principal vertebrate hosts in the sylvan cycle have not been identified, but there is evidence that dogs, horses, and donkeys may serve as reservoir hosts of this parasite (15). The existence of an urban cycle involving peridomestic sandfly species for L. (V.) braziliensis reflects the ability of these parasites and their vectors to adapt to changes in their original forested habitats with important public health implications. Studies using molecular techniques to characterize L. (V.) braziliensis populations from different regions have shown a relationship between level of similarity among the parasite populations (12, 24) and their geographic range, but recent data have also indicated that the considerable variability detected among these parasites is more probably related to the sandfly vector(s) and/or animal reservoir(s) involved in the transmission cycles (18).
Pathogens that produce many different genetic variants are more prone to infect multiple hosts (37). Although several studies have discussed the polymorphism observed in natural populations of different Leishmania species (8, 25, 34), until now there has been little information available about the genetic variability of the parasites and the correlation with ecoepidemiological features of the disease (16). The risk factors for clinical leishmaniasis are still poorly understood and probably are influenced by host and parasite features.
Considering the public health importance of leishmaniasis caused by L. (V.) braziliensis in Brazil and the role of the genetic polymorphism of the parasites in the epidemiology of the disease, we applied multilocus enzyme electrophoresis (MLEE) and the analysis of the restriction fragment length polymorphism (RFLP) of the internal transcribed spacers (ITS) of the rRNA genes as typing methods to determine the level of genetic variation in natural population of L. (V.) braziliensis derived from different endemic areas. MLEE are rapidly evolving markers and have been the most universally accepted means of characterizing and differentiating the many diverse types of leishmanial parasite (5, 30). The electrophoretic mobilities and banding patterns of several different enzymes are usually checked and compared to produce enzyme profiles that designate zymodemes. The relationships among the zymodemes are often verified by numerical taxonomic methods (5, 22, 30, 36). Recently, the ITS of the rRNA genes were reported to be useful markers to type Leishmania parasites (6, 33, 34). These regions are flanking the ribosomal genes, which are found in all eukaryotic organisms, and are clustered in tandem arrays on different chromosomes. The ITS regions of the rRNA genes evolve more rapidly than the functional domains that they flank, and phylogenetic studies of several organisms have demonstrated that they are useful to infer about the relationship among closely related populations of organisms (17).
The observed genetic variability and the relationships among the isolates were the basis for investigating whether the genetic differences in L. (V.) braziliensis reflect distinct ecoepidemiological features of the infection since these strains are endemic in the distinct areas studied, thus bringing to light fundamental information for future control programs of the disease.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Results of molecular typing and provenance details of 101 L. (V.) braziliensis strains collected over a 23-year period that were studied for their genetic relatedness
|
Intergenic region typing (IRT). Genomic RFLP analysis of the ITS between the small and large subunits of the ribosomal DNA (rDNA) locus (PCR-RFLP ITSrDNA) of leishmanial parasites used in the present study (Table 1) was done, as described elsewhere (6). The ITS were PCR amplified and digested with restriction enzymes (AluI, BstUI, EcoRI, FspI, HaeIII, HhaI, RsaI, and TaqI). Digestion products were separated by high-resolution electrophoresis in a 12% polyacrylamide gel by using the Genephor apparatus (Amersham Pharmacia Biotech). The fragment sizes were estimated relative to molecular weight markers and the banding patterns were used to group the isolates in genotypes (each genotype grouped isolates with the same banding pattern for all of the restriction enzyme used) and were also subjected to numerical analysis.
Numerical analysis of MLEE or PCR-RFLP ITSrDNA data. For this numerical analysis, a matrix with the presence or absence of the bands was built and analyzed by using the NTSYS software program (version 2.1; Exeter Software). The similarity matrices constructed by using the simple matching coefficient (Ssm = m/n, where m is the number of matches and n is the total number of characters) were transformed into phenograms by the UPGMA (unweighted pair-group method with arithmetic averages) algorithm (35). The confidence of the groups was tested applying a bootstrap analysis with 1,000 replicates using FreeTree software (version 0.9.1.50). The average and standard deviation of the similarities observed among individuals from each area (considering the clusters observed by the phenetic analysis) where the infection is endemic were calculated as a parameter of genetic diversity.
|
|
|---|
Affinities between zymodemes were calculated by using the simple matching similarity coefficient and were transformed into a phenogram by using mean distances between groups (Fig. 1A). The 15 L. (V.) braziliensis zymodemes could be clustered into two major groups, with high level of statistic support as demonstrated by the bootstrap value obtained for each group (Fig. 1A). Each group clustered strains collected on either the Atlantic coast or the Amazon region. One group (comprising isolates from CL or ML patients) was less diverse compared to the other group (comprising isolates from humans and various wild animals and sylvan sandflies). Furthermore, the later group contained two parasites (Z-53 and Z-69) that clustered independently of other zymodemes. These two zymodemes were assigned to the second group based on their geographic origin and because they presented a higher level of similarity with the zymodemes of this group than to the ones of the other group.
![]() View larger version (29K): [in a new window] |
FIG. 1. (A) MLEE phenogram showing the similarity between groups of L. (V.) braziliensis zymodemes. Unless indicated otherwise by the number in parentheses, IOC/Z contains a unique isolate. (B) Molecular tree showing the clustering of L. (V.) braziliensis genotypes classified in this work (PCR-RFLP ITSrDNA). The bootstrap values (italic numbers) are shown for the principal clusters. Details about strains are given in Table 1. The three main clusters represent parasites collected from different hosts and geographic areas in Brazil (i.e., the Atlantic coast [A and B] or the Amazon region [C]). Isolates characterized as IOC/Z-27 (outlined with continuous and dashed lines) presented genotypes in either group A and B. A, Strains collected in rural localities (Espírito Santo and Pernambuco) where different sandfly and animal species are involved in the transmission cycle; B, strains from urban areas (Espírito Santo and Rio de Janeiro) where transmission is related with a cycle involving domestic animals and preidomestic sandfly species; and C, strains from the Amazon region (in this case, the parasites are maintained in an enzootic cycle involving various wild animals and sylvan sand flies).
|
![]() View larger version (57K): [in a new window] |
FIG. 2. Acrylamide (12%) gel electrophoresis comparison of ITS fragment patterns generated with the restriction enzymes BstUI (A) and HhaI (B), among selected strains of L. (V.) braziliensis from Brazil. The codes in each lane represent the genotypes and the respective zymodemes (IOC/Z). The genotypic groups correspond to those described in Fig. 1. (A) Differences among groups A, B, and C; (B) differences among isolates from the same zymodeme.
|
Geographic distribution and genetic diversity of the genotypes. A geographic structuring was observed for some genotypes but not for others that was related to genotypes from more distant areas where the organisms were endemic (Fig. 1 and 3). The genetic similarity average based on the MLEE and IRT data was calculated for each locality, and it was assumed as a parameter to evaluate the level of genetic diversity. We observed that the isolates from the Amazonian region presented the highest level of genetic diversity and that the isolates from Rio de Janeiro were more homogeneous (Fig. 3). The isolates from Espírito Santo were considered as being involved in two distinct transmission cycles. The Espírito Santo isolates from urban areas grouped together with the isolates from Rio de Janeiro (group B) and those from rural areas clustered with isolates from Pernambuco (group A) (Fig. 3).
![]() View larger version (53K): [in a new window] |
FIG. 3. Map of Brazil showing the areas from which the L. (V.) braziliensis isolates originated. A schematic and summarized dendrogram showing the relationship among the genotypic groups (A, B, and C as defined in Fig. 1) was superimposed. The level of genetic diversity found among the isolates from each area where the parasites were endemic was estimated by using the average of similarities as follows (averages ± the standard deviations of MLEE and IRT, respectively): Amazonia (0.75 ± 0.082 and 0.76 ± 0.095), Pernambuco (0.92 ± 0.01 and 0.83 ±0.082), Espírito Santo A (1.0 and 0.89 ± 0.05), Espírito Santo B (1.0 and 0.83 ± 0.11), and Rio de Janeiro (1.0 and 0.93 ± 0.035). The Brazilian states where the isolates were collected included Amazonia, Rondônia, Pernambuco, Espírito Santo, and Rio de Janeiro (see Table 1).
|
|
|
|---|
The present systematic molecular study revealed L. (V.) braziliensis strains to be an extremely diverse population. The multilocus data on a diverse collection of strains from different hosts and ecologic regions in Brazil indicate a clonal population structure, with some zymodemes (e.g., IOC/Z-27) widely distributed and others seemingly unique and localized to a particular endemic focus (e.g., IOC/Z-35 and IOC/Z-45). The IRT analysis resulted in a very high degree of discrimination among isolates from diverse regions, confirming the ability of this technique to discriminate between closely related parasites (6). The three distinct genotypic groups described herein may represent important genetic variation within L. (V.) braziliensis, which could explain the plasticity of these parasites and of their ability to adapt to changing ecological conditions. Changes in the ITSrDNA locus, which is a noncoding region, should not reflect the ability to adapt to ecological changes. However, because of the high level of polymorphism in the composition of the ITSrDNA locus, several molecular epidemiological studies have been conducted by using this tool. Analysis of the ITS regions of the nuclear rRNA operon for different parasites has demonstrated the utility of this marker in detect genetic diversity and the association of this with several epidemiological aspects (9, 26)
The present study also indicates that most of the genotypes are specific to some geographic areas (e.g., genotype B1 contains isolates from one geographic location, in Espírito Santo, that were collected over a distinct period of time [between 1982 and 1995]), showing the importance of clonal dissemination versus sexual reproduction for this organism in nature, at least occurring in some regions. The stability of this genotype adapting to some environmental modifications through years in the same area is a clear proof of the aforementioned statement.
L. (V.) braziliensis can infect vertebrate hosts of different species and orders (21) and has been associated with a number of different sandfly species (15). Our analyses indicated that the molecular diversity found in parasites from the Amazon Basin is apparently related to the great number of sandfly vector(s) and/or animal reservoir(s) involved in the transmission cycles (13, 14). In contrast, the L. (V.) braziliensis populations circulating in the Brazilian Atlantic coast (in the present study, representing 81.2% of the isolates collected from human beings and dogs in old established communities in nonforested areas) showed a lower level of heterogeneity than did the Amazonian strains and have peridomestic sandfly species such as Lutzomyia intermedia and Lutzomyia whitmani as the principal suspected vectors (1, 3, 10). The presence of the same genotypes in humans and dogs in the same area suggests that specific transmission cycles could be defined by using the proposed methodologies.
A total of 49% of the parasites collected in areas (Rio de Janeiro or Espírito Santo) where sylvan animals are probably not involved in the transmission cycle were clustered together into the same zymodeme (IOC/Z-27) of the L. (V.) braziliensis reference strain (MHOM/BR/1966/M2903), which was isolated from a patient with American CL in the Serra dos Carajás, Pará State. In this case, a sylvan sandfly, Psychodopygus wellcomei, has been determined to be the vector (20, 31). In contrast, isolates from Pernambuco, where grass mice (Bolomys lasiurus) and black rats (Rattus rattus) have been found to be parasite reservoirs (S. P. Brandão-Filho, M. E. F. Brito, F. G. Carvalho, E. A. Ishikawa, E. Cupolillo, and J. J. Shaw, Abstr. Wordleish 2000, abstr. 100, 2000), were grouped into five distinct zymodemes. Factors that predispose pathogens to infect multiple hosts include high levels of genetic diversity and opportunities for cross-species transmission (37). The extent of host adaptation is therefore linked to and limited by the genetic variability. Pathogens that contain more accumulating mutations should produce more genetic variants and are more likely to be generalists. An alternative explanation to the extent of host adaptation of some organisms is to consider that the infection of a given host is only related to opportunity; therefore, parasites that infect several different hosts are more likely to be isolated in different transmission cycles.
The RFLP-ITSrDNA data show that parasites from urban areas (Espírito Santo and Rio de Janeiro), where transmission is associated with Lutzomyia intermedia, cluster together. Another cluster comprises parasites from rural localities in Espírito Santo and Pernambuco, where Lutzomyia whitmani is the principal vector species. Other studies using different molecular markers have demonstrated a geographic structuring of Leishmania parasites (2, 16, 23, 33, 34), and populations of L. (V.) braziliensis isolates collected from geographically closed areas seem to be more genetically similar (12, 24). However, the present study indicates that the geographic structuring of L. (V.) braziliensis reflects the plasticity of these parasites to adapt to different vector species (or populations of this vectors) involved in the transmission cycle. Furthermore, our results indicate an association between the diversity of vectors involved in the transmission of L. (V.) braziliensis and the genetic diversity of the parasites. Evaluation of the phosphoglycan composition of these leishmanial populations is fundamental to an understanding of the vector-parasite interaction (32). Another important approach is to study the population structure of the sandfly vectors in areas where the genetic diversity of L. (V.) braziliensis has been determined.
It was recently suggested that the dispersion of Leishmania clones is associated with the vector and reservoir movements (18). However, what is important to note is the plasticity of leishmanial parasites, their ability to adapt to changes in ecologic conditions, and a possible evolutionary relationship between them and their vectors as observed in the Old World (28).
The clinical expression of leishmanial infection is dependent on a number of different factors (15). The identification of virulence factors that contribute to the pathogenesis of L. (V.) braziliensis infections are poorly understood, mainly because of the lack of good in vivo models. A better understanding of the genetic basis of pathogenesis in ML is needed, as are studies on the role of host factors in immune susceptibility to this form of the disease. Our genotyping studies did not identify a pathogenic clone specifically associated with the development of ML. However, in other analyses we were able to detect distinct genotypic types of a single L. (V.) braziliensis isolate from a patient with ML (unpublished data). Competitive interactions among coinfecting clones could affect host immune responses and play a role in determining the immunopathology of the progressive and destructive form of ML. Some underlying causes of its hyperergic response might include the following factors: (i) the resistance of some parasite clones to elimination (27) and the persistence of "allergic" antigen that evokes hypergic hypersensitivity inflammatory responses and (ii) autoimmune phenomena related to antigens cross-reactive between leishmanial parasites and host tissues (11, 29).
In conclusion, the combination of MLEE and IRT-ITSrDNA is useful for studies of the interface of population genetic and epidemiology. The application of these molecular techniques could allow us to better determine the genetic diversity of various Leishmania parasites circulating in nature. This information is expected to yield new approaches to mitigate parasite transmission and virulence in humans and eventually to improve parasite control.
This work was supported by Fundação Oswaldo Cruz, Fundação de Amparo a Pesquisa do Estado do Rio de Janeiro FAPERJ, Conselho Nacional de Desenvolvimento Científico e Tecnológico CNPq, and PRONEX III/CNPq/MCT (Brazil).
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»