- Open Access
American tegumentary leishmaniasis caused by Leishmania (Viannia) braziliensis: assessment of parasite genetic variability at intra- and inter-patient levels
Parasites & Vectorsvolume 6, Article number: 189 (2013)
The genetic variability of Leishmania (Viannia) braziliensis was assessed at intra and interpatient levels of individuals with different clinical manifestations of American tegumentary leishmaniasis (ATL).
Fifty-two samples, of which 13 originated from cutaneous lesions and 39 from mucosal lesions, provided by 35 patients, were examined by low-stringency single-specific-primer PCR (LSSP-PCR) and phenetic analysis. Genetic variability of L. (V.) braziliensis, in kinetoplast DNA (kDNA) signatures, was compared both from different patients and from different lesions of the same patient. Phenetic analysis was performed to evaluate the degree of heterogeneity of the kDNA minicircles. In order to evaluate inter and intrapatient L. (V.) braziliensis genetic variability, the percentage of shared bands and analysis of the coefficients of similarity were analyzed.
Different genetic profiles, representing kDNA signatures of the parasite, were obtained by LSSP-PCR analysis of each sample. Phenetic analysis grouped genetic profiles of different levels of differentiation from more similar to most divergent. The percentage of shared bands at the inter and intrapatient levels was 77% and 89%, respectively. Comparison of the average inter and intrapatient coefficients of similarity and their standard deviations were statistically significant (p < 0.001).
Genetic variability at the intrapatient level was less pronounced than that between different patients. A conceptual model was proposed to better understand the complexity at both levels.
Leishmaniasis is a disease ranging from potential self-healing cutaneous lesions to mucosal disorders to more severe visceral affliction possibly resulting in death, all depending upon the species of Leishmania, the protozooan etiological agent.
American tegumentary leishmaniasis (ATL) caused by Leishmania (Viannia) braziliensis is characterized by chronicity, latency and metastatic tendencies. In some cases, after the appearance of an initial cutaneous lesion, multiple cutaneous lesions (disseminated leishmaniasis, DL) and/or mucosal lesions (mucosal leishmaniasis, ML) can arise as a consequence of the dissemination of the parasite through the blood and lymphatic systems [1–3]. Cutaneous lesions can remain active for many years and can coexist with mucosal lesions, as exhibited in the mucocutaneous leishmaniasis (MCL) form . Less than 5% of cutaneous leishmaniasis cases develop destructive lesions in the upper respiratory tract, primarily in the nasal mucosa , and more than 25% of the DL cases present damaged mucosa simultaneously with cutaneous manifestations . However, approximately 16% of the ML cases do not present any previous history of an initial cutaneous lesion .
Factors such as the immunological response of the patient and the genetic constitution of the parasite could influence the clinical manifestations and prognosis of this disease. It is also quite possible that these factors play a role in parasite persistence. Parasite persistence and reactivation have been reported to be contributing factors towards recurring episodes of ATL [8–10]. However, the mechanism of the dissemination and emergence of L. (V.) braziliensis to new cutaneous and mucosal tissue, the pathogenic process as well as the natural history of the recurrent leishmaniasis is not yet fully understood.
Human infection by L. (V.) braziliensis produces a broad spectrum of clinical manifestations that can be attributed to intraspecific variability of the parasite. The contribution of the parasite to the diversity of ATL patient clinical conditions has been examined in several studies aiming to investigate whether certain species and/or clonal populations of Leishmania are correlated with any clinical manifestations of the disease. Intraspecific heterogeneity among members of the subgenera Viannia and Leishmania was detected by genetic and molecular studies, which have afforded a better understanding of their epidemiology in diverse endemic areas [11–14].
Different polymorphic molecular markers based on the polymerase chain reaction (PCR) have been used in several studies of the L. (V.) braziliensis genetic variability and molecular epidemiology. The importance of this particular species was demonstrated by the existence of various intraspecific genetic polymorphisms and the variability of the clinical manifestations observed in ATL. Previous studies have established a correlation between genetic polymorphisms and the eco-epidemiological data of isolates [15–17]. Such studies reinforce the clonal population structure of Leishmania, as proposed by Tibayrenc et al. . Leishmania species naturally circulate as a group of heterogeneous subpopulations, and demonstration of the polyclonality of the initial inoculum of multiple strains has been reported . These mechanisms could influence the clinical prognosis of the disease and the efficiency of therapeutic strategies.
Our group recently introduced the application of low-stringency single-specific-primer PCR (LSSP-PCR) towards the investigation of intraspecific polymorphisms in the variable region of the kinetoplast DNA (kDNA) minicircles of L. (V.) braziliensis in ATL patients [20, 21]. The purpose of the present study is to further investigate parasite subpopulations, with a particular focus on the genetic variability of L. (V.) braziliensis at intra and interpatient levels, using the same molecular marker associated with phenetic analyses. Towards this end, samples obtained from cutaneous and mucosal lesions of patients with different clinical manifestations of ATL from Rio de Janeiro, Brazil were analyzed.
This study included 35 patients from Rio de Janeiro, Brazil, presenting ATL and mucosal involvement, each attended at the Laboratório de Vigilância em Leishmanioses/Lab. VigiLeish, Instituto de Pesquisa Clínica Evandro Chagas – IPEC/Fiocruz, Rio de Janeiro, from January 2005 to December 2009. Patients were grouped according to their clinical manifestations in the following three categories: i) mucocutaneous leishmaniasis (MCL) for patients that presented cutaneous and mucosal lesions in contiguous and concomitant forms, ii) mucosal leishmaniasis (ML) for patients that presented only mucosal lesions (nasal and/or oral) and iii) disseminated leishmaniasis (DL) for patients presenting more than 10 cutaneous lesions and mucosal involvement. All patients were subjected to dermatological and otolaryngological examination, including direct inspection of the upper respiratory tract using rigid optic Karl Storz 70° and 30°, to assess, respectively, the presence or absence of larynx and pharyngeal or nasal lesions. This study was approved by the Ethical Committee of Research at IPEC/FIOCRUZ (process 0016.0.009-02), and all of the patients signed consent forms.
Molecular diagnosis by PCR and low-stringency single-specific-primer PCR (LSSP-PCR) reactions
Patient biopsy samples from cutaneous and/or mucosal lesions were processed for molecular analysis. Extraction of DNA and molecular diagnostic PCR-based assays were performed as previously described [21, 22]. The 750 bp fragments generated by specific PCR were purified with a Wizard PCR Prep System (Promega, Madison, WI) according to the manufacturer’s instructions. The purified products were then submitted to a second round of amplification (LSSP-PCR) with only one specific primer under low stringency conditions. Reaction conditions and visualization of the products were performed as described by Oliveira et al. .
Bands varying in size from 350–750 bp were selected for phenetic analysis. Analysis of the profiles generated by LSSP-PCR was performed by calculating the levels of similarity using the Simple Matching Coefficient (Sm). The matrix of similarity was transformed into a dendrogram through the Unweighted Pair Group Method Arithmetical Average (UPGMA) algorithm with the NTSYS version 2.1 program (Exeter Software).
Comparison of inter and intrapatient genetic variability
Interpatient genetic variability was evaluated by: i) obtaining the coefficients of similarity between the patient sample and all of the other samples (from all of the patients), ii) calculating the arithmetic averages and their respective standard deviations of those interpatient coefficients of similarity and iii) calculating the average of those mean values, with the aim of estimating the average interpatient similarity (i.e. the variability existing among parasites recovered from lesions of different patients). Intrapatient genetic variability was estimated by calculating the average of the mean coefficient of similarity of each patient, taking into account the cutaneous and mucosal lesions present in the same patient. The averages and standard deviations (SD) of the inter and intrapatient coefficients of similarity were compared applying the stud ent’s t-test.
Comparison of inter- and intra-patient similarities was also performed by analyzing the percentages of shared characters (bands), defined by genetic profiles generated by LSSP-PCR. The sharing frequency of each character between interpatient samples was calculated after first calculating the average of each sharing frequency. The same analysis was performed to compare different lesions from the same patient.
Clinical-epidemiological characteristics of the patients
The clinical forms of ATL exhibited by the patients were: 16 with MCL, 17 with ML and 2 with DL. Of these patients, 25 (71.4%) were male and 10 (28.6%) were female. Patient age varied from 12 to 77 years. The presence of lesions in the nasal mucosa was observed in 77% of the patients. The two patients with DL (patients 34 and 35) exhibited concomitant cutaneous and mucosal lesions. Only one patient with MCL (patient 5) displayed concomitant nasal and oral mucosal lesions. Of the 17 patients with ML, 12 reported the absence of an initial cutaneous lesion or scar suggestive of leishmaniasis. Data are shown in Table 1.
Molecular diagnostics and analysis of genetic variability
From the 35 patients, 52 tissue samples were analyzed, of which 13 originated from cutaneous lesions and 39 from mucosal lesions. Samples from cutaneous and mucosal lesions were only available from 9 of the 16 patients with MCL. From 2 of the patients with ML (patients 31 and 32), fragments of the mucosal nasal lesions were collected at different times over a period of one year. From patient 34 with DL, three samples of cutaneous lesions were analyzed, two of which were collected at the same time as the sample of the mucosal lesion and one that was collected 6 months after treatment.
Internal controls for PCR inhibition were performed by amplifying a 110 bp fragment of the human β-globin gene. The expected 110 bp fragment was obtained in each sample tested (data not shown), confirming the absence of internal inhibitors. The 750 bp diagnostic band was evident in all samples analyzed (data not shown).
Different genetic profiles, representing kDNA signatures of the parasite, were obtained by LSSP-PCR analysis of each sample. Figure 1-A shows selected representative results. Phenetic analysis grouped the LSSP-PCR profiles according to their coefficients of similarity (Sm), which diverged from 0.53 to 0.94 (Figure 1-B). The dendrogram illustrates that in two patients with MCL (patients 12 and 11) samples of cutaneous and mucosal lesions were genetically similar to each other, sharing a coefficient of similarity of 0.94. In two other patients, one with MCL (patient 5) and one with DL (patient 35), the nasal and oral mucosal lesions also presented similar genetic profiles (Sm = 0.94), although somewhat differing genetic profiles (Sm = 0.78 and Sm = 0.82, respectively) were apparent in the respective cutaneous lesions. Patients 31 and 27, with the ML clinical condition, also showed similar genetic profiles in the samples of their mucosal nasal lesions. On the other hand, distinct genetic profiles were detected in samples from the cutaneous and mucosal lesions of patients 14, 3, 8 and in the nasal mucosal lesions of patients 31 and 32, which were collected at different times over a period of one year. Interestingly, very divergent profiles were observed in the cutaneous and mucosal lesion samples of patient 34 and also in the mucosal lesion samples of patients 31, 17 and 32.
In order to evaluate interpatient genetic variability, the percentage of shared bands was analyzed. This approach was validated by first demonstrating that, on average, each character (band) was shared by 77% of the patients. Analysis of interpatient coefficients of similarity revealed a similar average of 0.69 (SD = 0.05). On the other hand, the average characters shared at the intrapatient level was 89%, and the average of the average coefficients of similarity was 0.80 (SD = 0.12). Comparison of the average inter and intrapatient coefficients of similarity and their standard deviations were statistically significant (p < 0.001). Table 2 represents the inter and intrapatient values of the coefficient of similarity for selected representative samples.
The conceptual model
We have proposed a conceptual model (Figure 2) as a way to explain the dynamics of the circulating populations of L. (V.) braziliensis at both intra and interpatient levels in a specific endemic area.
Our group has recently been applying different molecular tools towards the study of specific diagnostics and the genotypic behaviour of L. (V.) braziliensis subpopulations in patients with mucocutaneous leishmaniasis [20, 21]. These studies served as a basis for the present investigation, which aimed to explore the relationship between genetic diversity of L. (V.) braziliensis at inter and intrapatient levels presenting distinct clinical manifestations of ATL. At the same time, this study also aimed to evaluate such genetic diversity following two hypothetical scenarios: (i) the individual, by him/herself, acting as a biological filter whereby highly adapted populations are selected and (ii) eco-epidemiological interactions occurring in the endemic area as reflecting on the dynamics of the circulating populations of L. (V.) braziliensis. According to our conceptual model, the situation at intrapatient level takes into account the participation of both parasite virulence factors and the immunological system of the host and that the reactivation of the disease would reflect in the emergence of different parasite subpopulations. However, at the interpatient level circumstance is that heterogeneous subpopulations circulate among man and vectors as well as domestic and sylvatic animals in simple and complex epidemiological interactions. In such specific cases individuals would be exposed to distinct selective pressures, to the polyclonality of the initial inoculum and also to multiple independent infections. If the biological entity that evolves is the clonal lineage , one can imagine different clones being reshuffled and propagated in the natural circulating population, sharing different factors of virulence or competing for selective advantages such as distinct growth rates and differential tropisms [23, 24].
Genetic variability of L. (V.) braziliensis, demonstrated here by polymorphisms of kDNA signatures, was identified among the patients and among different lesions of the same patient. Phenetic analysis was adopted to evaluate the degree of heterogeneity of the kDNA minicircles and group them into genetic profiles of different levels of differentiation, from more similar to more divergent. Although this study did not uncover a relationship between parasite genotype and clinical manifestations of ATL, considerably divergent genetic profiles have been detected in nasal mucosal lesions of patients with ML. Schriefer et al.  noted that a higher frequency of a determined genotype of L. (V.) braziliensis in certain isolates of ML originated from the same endemic area in the state of Bahia.
This idea of multiple independent infections may be an explanation corroborated by data from patient 34, with DL, whose results display significant genetic diversity in the minicircles from distinct cutaneous and mucosal lesions collected both at the same time and also at different moments over a short period of time.
Despite the existence of heterogeneity, genotypic similarities between the samples from cutaneous and mucosal lesions were evident on the intrapatient level. The average of the coefficients of similarities (0.80) was also of major significance, indicating that the samples are genetically related. Therefore, the intrapatient genetic variability, less than that between different patients (77%), must be a result that is also in accordance with the percentage of shared intrapatient characters (89%).
Elegant studies have demonstrated that the host immune system plays an important role in the selection of parasite populations [26–28]. Considering this, it is possible that man may act as a biological filter towards the selection of parasite subpopulations most adapted to develop in the host organism, and in this case, subpopulations can be responsible for either cutaneous or mucosal pathologies as mentioned in our conceptual model.
From the epidemiological point of view, the complex cycle of transmissions in which the individuals are exposed is a very important factor for the possibility of reinfections through the accumulation of multiple independent infections in patients residing in endemic areas [8, 29]. The complexity of transmission cycles, involving sylvatic and domestic mammals together with the phlebotamine vectors surely contributes to the large diversity of the natural population of L. (V.) braziliensis[30, 31]. An understanding of the genetic diversity of circulating parasite populations in a specific endemic area, as well as an investigation of the clonality of the initial lesions, is fundamentally important and assists us in interpreting the inter and intrapatient variability.
This study addresses the genetic variability of L. (V.) braziliensis, demonstrated by polymorphisms of kDNA signatures, not only among the patients but also among different lesions of the same patient. Despite the existence of heterogeneity, genotypic similarities between the samples from cutaneous and mucosal lesions were evident at the intrapatient level, the genetic variability being less pronounced than the variability between different patients. A conceptual model was proposed to better understand the complexity at both levels.
Costa JM, Marsden PD, Llanos-Cuentas EA, Netto EM, Carvalho EM, Barral A, Rosa AC, Cuba CC, Magalhães AV, Barreto AC: Disseminated cutaneous leishmaniasis in a field clinic in Bahia, Brazil: a report of eight cases. J Trop Med Hyg. 1986, 89: 319-323.
Martinez JE, Arias AL, Escobar MA, Saravia NG: Haemoculture of Leishmania (Viannia) braziliensis from two cases of mucosal leishmaniasis: re-examination of haematogenous dissemination. Trans R Soc Trop Med Hyg. 1992, 86: 392-394. 10.1016/0035-9203(92)90233-3.
Barral A, Guerreiro J, Bomfim G, Correia D, Barral-Netto M, Carvalho EM: Lymphadenopathy as the first sign of human cutaneous infection by Leishmania braziliensis. AmJTrop Med Hyg. 1995, 53: 256-259.
Marsden PD, Llanos-Cuentas EA, Lago EL, Cuba CC, Barreto AC, Costa JM: Human mucocutaneous leishmaniasis in Três Braços, Bahia-Brazil. An area of Leishmania braziliensis braziliensis transmission. III Mucosal disease presentation and initial evolution. Rev Soc Bras Med Trop. 1984, 17: 179-186. 10.1590/S0037-86821984000400004.
Marsden PD: Mucosal leishmaniasis (“espundia” Escomel, 1911). Trans R Soc Trop Med Hyg. 1986, 80: 859-876. 10.1016/0035-9203(86)90243-9.
Turetz ML, Machado PR, Ko AI, Alves F, Bittencourt A, Almeida RP, Mobashery N, Johnson WD, Carvalho EM: Disseminated leishmaniasis: a new and emerging form of leishmaniasis observed in northeastern Brazil. J Infect Dis. 2002, 186: 1829-1834. 10.1086/345772.
Weigle K, Saravia NG: Natural history, clinical evolution, and the host-parasite interaction in New World cutaneous leishmaniasis. Clin Dermatol. 1996, 14: 433-450. 10.1016/0738-081X(96)00036-3.
Saravia NG, Weigle K, Segura I, Gianninim SH, Pacheco RS, Labrada LA, Goncalves A: Recurrent lesions in human Leishmania braziliensis infection – reactivation or reinfection?. Lancet. 1990, 336: 398-402. 10.1016/0140-6736(90)91945-7.
Schubach AO, Haddad F, Oliveira-Neto MP, Degrave W, Pirmez C, Grimaldi G, Fernandes O: Detection of Leishmania DNA by polymerase chain reaction in scars of treated human patients. J Infect Dis. 1998, 178: 911-914. 10.1086/515355.
Mendonça MG, Brito MEF, Rodrigues EHG, Bandeira V, Jardim ML, Abath FGC: Persitence of Leishmania parasites in scars after clinical cure of American cutaneous leishmaniasis: is there a sterile cure?. J Infect Dis. 2004, 189: 1018-1023. 10.1086/382135.
Pacheco RS, Lopes UG, Morel CM, Grimaldi G, Momen H:Schizodeme analysis ofLeishmaniaisolates and comparison with some phenotypic techniques. Leishmania. Taxonomy and Phyogeny. Edited by: Rioux JA. 1986, Montpellier: IMEEE, 57-65.
Tomaz-Soccol V, Lanotte G, Rioux J, Pratlong F, Martini-Dumas A, Serres E: Monophyletic origin of the genus Leishmania Ross, 1903. Ann Parasitol Hum Comp. 1993, 68: 107-108.
Cupolillo E, Grimaldi G, Momen H: A general classification of New World Leishmania using numerical zymotaxonomy. AmJTrop Med Hyg. 1994, 50: 296-311.
Pacheco RS, Fernandes O, Salinas G, Segura I, Momen H, Degrave W, Saravia NG, Campbell DA: Intraspecific heterogeneity in the mini-exon gene localization of Leishmania (Viannia) panamensis and Leishmania (Viannia) guyanensis from Colombia. J Parasitol. 2000, 86: 1250-1253.
Gomes RF, Macedo AM, Pena SD, Melo MN:Leishmania(Viannia)braziliensis: genetic relationships between strains isolated from different areas of Brazil as revealed by DNA fingerprinting and RAPD.Exp Parasitol. 1995, 80: 681-687. 10.1006/expr.1995.1084.
Ishikawa EA, Silveira FT, Magalhães AL, RB G j, Melo MN, Gomes R, Silveira TG, Shaw JJ: Genetic variation in populations of Leishmania species in Brazil. Trans R Soc Trop Med Hyg. 2002, 96: S111-S121.
Cupolillo E, Brahim LR, Toaldo CB, Oliveira-Neto MP, Brito ME, Falqueto A, Naiff MF, Grimaldi G: Genetic polymorphism and molecular epidemiology of Leishmania (Viannia) braziliensis from different hosts and geographic areas in Brazil. J Clin Microbiol. 2003, 41: 3126-3132. 10.1128/JCM.41.7.3126-3132.2003.
Tibayrenc M, Kjellberg F, Ayala FJ: A clonal theory of parasitic protozoa: the population structures of Entamoeba, Giardia, Leishmania, Naegleria, Plasmodium, Trichomonas, and Trypanosoma and their medical and taxonomical consequences. Proc Natl Acad Sci U S A. 1990, 87: 2414-2418. 10.1073/pnas.87.7.2414.
Pacheco RS, Grimaldi G, Momen H, Morel CM: Population heterogeneity among clones of New World Leishmania species. Parasitology. 1990, 100: 393-398. 10.1017/S0031182000078677.
Baptista C, Schubach AO, Madeira MF, Leal CA, Pires MQ, Oliveira FS, Conceição-Silva F, Valete-Rosalino CM, Salgueiro MM, Pacheco RS: Leishmania (Viannia) braziliensis genotypes identified in lesions of patients with atypical or typical manifestations of tegumentary leishmaniasis: evaluation by two molecular markers. Exp Parasitol. 2009, 121: 317-322. 10.1016/j.exppara.2008.12.006.
Oliveira FS, Valete-Rosalino CM, Schubach AO, Pacheco RS: kDNA minicircle signatures of Leishmania (Viannia) braziliensis in oral and nasal mucosa from mucosal leishmaniasis patients. Diagn Microbiol Infect Dis. 2010, 66: 361-365. 10.1016/j.diagmicrobio.2009.11.004.
Oliveira FS, Pirmez C, Pires MQ, Brazil RP, Pacheco RS: PCR-based diagnosis for detection of Leishmania in skin and blood of rodents from an endemic area of cutaneous and visceral leishmaniasis in Brazil. Vet Parasitol. 2005, 15: 219-227.
Ayala FJ: Trypanosoma and Leishmania have clonal population structures of epidemiological significance. Biol Res. 1993, 26: 47-63.
Pacheco RS, Brito CM:Reflections on the population dynamics ofTrypanosoma cruzi: heterogeneity versus plasticity.Mem Inst Oswaldo Cruz. 1999, 94: 199-201.
Schriefer A, Schriefer AL, Góes-Neto A, Guimarães LH, Carvalho LP, Almeida RP, Machado PR, Lessa HA, Jesus AR, Riley LW: Multiclonal Leishmania braziliensis population structure and its clinical implication in a region of endemicity for American tegumentary leishmaniasis. Infect Immun. 2004, 72: 508-514. 10.1128/IAI.72.1.508-514.2004.
Jiménez M, Ferrer-Dufol M, Cañavate C, Gutiérrez-Solar B, Molina R, Laguna F, López-Vélez R, Cercenado E, Daudén E, Blázquez J: Variability of Leishmania (Leishmania) infantum among stocks from immunocompromised, immunocompetent patients and dogs in Spain. FEMS Microbiol Lett. 1995, 131: 197-204.
Macedo AM, Pena SDJ: Genetic variability of trypanosoma cruzi: implications for the pathogenesis of Chagas’ disease. Parasitol Today. 1998, 14: 119-124. 10.1016/S0169-4758(97)01179-4.
Brito CMM, Pires MQ, Pacheco RS: Chagas disease and HIV co-infection: genetic analyses of two Trypanosoma cruzi strains under experimental immunosuppression. Kinetoplastid Biol Dis. 2003, 2: 17-23. 10.1186/1475-9292-2-17.
Gontijo CM, Silva ES, Fuccio MB, Sousa MC, Pacheco RS, Dias ES, Andrade Filho JD, Brazil RP, Melo MN: Epidemiological studies of an outbreak of cutaneous leishmaniasis in the Rio Jequitinhonha Valley, Minas Gerais, Brazil. Acta Trop. 2002, 81: 143-150. 10.1016/S0001-706X(01)00205-4.
Lainson R, Shaw JJ:New World Leishmaniasis – The neotropicalLeishmaniaspecies. Topley & Wilson´s Microbiology & Microbial Infections. Edited by: Cox FEG, Kreier JP, Wakelin D. 1998, Arnold, London, United Kingdom, 241-266.
Brito ME, Andrade MS, Mendonça MG, Silva CJ, Almeida EL, Lima BS, Félix SM, Abath FG, Graça GC, Porrozzi R: Species diversity of Leishmania (Viannia) parasites circulating in an endemic area for cutaneous leishmaniasis located in the Atlantic rainforest region of northeastern Brazil. Trop Med Int Health. 2009, 14: 1278-1286. 10.1111/j.1365-3156.2009.02361.x.
The authors acknowledge the team of dermatologists and otolaryngologists at the Laboratório de Vigilância em Leishmanioses – IPEC/Fiocruz for collecting patient biopsies. This work received financial support from Instituto Kinder do Brasil (IKB), Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq), and Fundação Carlos Chagas de Apoio a Pesquisa do Estado do Rio de Janeiro (FAPERJ). R.S. Pacheco and A.O. Schubach are investigators of CNPq.
English review and revision by Mitchell Raymond Lishon, native of Chicago, Illinois, USA – UCLA, 1969.
The author declared that they have no competing interest.
FSO and RSP conceived the study; FSO and RSP designed the study; FSO performed the PCR assays and the genetic characterization of the parasites; CMVR was responsible for othorinolaryngologic evaluation and collection of the biological samples; SJBP and FACC contributed design and statistical analysis; RSP, CMVR and AOS analysed and interpreted the data; FSO drafted the manuscript. RSP critically revised the manuscript for intellectual content. All authors read and approved the final manuscript. FSO and RSP are guarantors of the paper. RSP and AOS are investigators of CNPq, Brazil.