Search tips
Search criteria 


Logo of plosntdPLoS Neglected Tropical DiseasesView this ArticleSubmit to PLoSGet E-mail AlertsContact UsPublic Library of Science (PLoS)
PLoS Negl Trop Dis. 2017 November; 11(11): e0006095.
Published online 2017 November 20. doi:  10.1371/journal.pntd.0006095
PMCID: PMC5714378

Prevalence and risk factors associated with Leishmania infection in Trang Province, southern Thailand

Jipada Manomat, Data curation, Formal analysis, Project administration, Writing – original draft,1 Saovanee Leelayoova, Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Supervision, Validation, Writing – original draft, Writing – review & editing,2 Lertwut Bualert, Data curation, Project administration, Validation,3 Peerapan Tan-ariya, Conceptualization, Supervision, Validation, Writing – review & editing,1 Suradej Siripattanapipong, Formal analysis, Validation, Writing – review & editing,1 Mathirut Mungthin, Conceptualization, Data curation, Methodology, Supervision, Writing – review & editing,2 Tawee Naaglor, Formal analysis, Validation,2 and Phunlerd Piyaraj, Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Writing – original draft2,*
Fabiano Oliveira, Editor



Autochthonous cutaneous and visceral leishmaniasis (VL) caused by Leishmania martiniquensis and Leishmania siamensis have been considered emerging infectious diseases in Thailand. The disease burden is significantly underestimated, especially the prevalence of Leishmania infection among HIV-positive patients.


A cross-sectional study was conducted to determine the prevalence and risk factors associated with Leishmania infection among patients with HIV/AIDS living in Trang province, southern Thailand, between 2015 and 2016. Antibodies against Leishmania infection were assayed using the direct agglutination test (DAT). DNA of Leishmania was detected by ITS1-PCR using the buffy coat. Species of Leishmania were also identified.


Of 724 participants, the prevalence of Leishmania infection was 25.1% (182/724) using either DAT or PCR assays. Seroprevalence of Leishmania infection was 18.5% (134/724), while Leishmania DNA detected by the PCR method was 8.4% (61/724). Of these, 24.9% (180/724) were asymptomatic, whereas 0.3% (2/724) were symptomatic VL and VL/CL (cutaneous leishmaniasis). At least five species were identified: L. siamensis, L. martiniquensis, L. donovani complex, L. lainsoni, and L. major. Multivariate analysis showed that CD4+ levels <500 cells/μL and living in stilt houses were independently associated with Leishmania infection. Those who were PCR positive for Leishmania DNA were significantly associated with a detectable viral load, whereas non-injection drug use (NIDU) and CD4+ levels <500 cells/μL were potential risk factors of Leishmania seropositivity.


A magnitude of the prevalence of underreporting Leishmania infection among Thai patients with HIV was revealed in this study. Effective public health policy to prevent and control disease transmission is urgently needed.

Author summary

Visceral leishmaniasis (VL) in Thailand is caused by two causative agents, Leishmania martiniquensis and Leishmania siamensis. A public health concern brought us to investigate the magnitude of Leishmania infection among individuals with HIV living in an affected area, Trang province, southern Thailand. The results showed a high seroprevalence of Leishmania infection. Using PCR-based technique, DNA detection in the buffy coat was 8.4%, and 1.8% were symptomatic VL. Asymptomatic Leishmania infection could play an important role in disease transmission. Risk factors associated with Leishmania infection among Thai patients with HIV were firstly described. Those who were NIDU and lived in stilt houses were associated with Leishmania infection. Individuals who had lower immunity with detectable viral load were more likely to contract the infection. Interestingly, not only L. martiniquensis and L. siamensis were identified but also L. donovani complex, L. major, and L. lainsoni, were firstly reported among indigenous Thai people. These findings could lead to effective intervention and prevention methods to control leishmaniasis in Thailand. Further studies are needed to investigate the disease development in asymptomatic infections as well as evaluate the prevalence in other regions of Thailand.


Co-infection of leishmaniasis and human immunodeficiency virus (HIV) is a major public health problem globally. Leishmania and HIV each promote the activation of the other, causing host immune impairment. The co-infection results in treatment failure, high relapse, and high mortality rate [1]. Meta-analysis has revealed that the direct agglutination test (DAT) gave high sensitivity and specificity for serodiagnosis of VL when compared to other serological tests [2]. However, low sensitivity of serological tests for VL diagnosis has been shown among these patients due to defective host immunity [3]. To increase the sensitivity of Leishmania DNA detection, the polymerase chain reaction (PCR) using blood samples has been suggested [4]. The internal transcribed spacer 1 (ITS1)-PCR method has been recommended to detect Leishmania DNA [5].

In Thailand, the first autochthonous VL case was documented in a 3-year-old girl living in a southern province in 1999 [6]. Until 2012, L. siamensis was firstly reported in a patient with HIV in Trang province. Since then, CL and/or VL have been sporadically reported in immunocompetent and immunocompromised patients predominantly in the south and north of Thailand and about 40% were patients with HIV/AIDS. L. martiniquensis was the predominant causative agent while L. siamensis was indigenously reported in only one Thai patient [7].

Information of the true prevalence of Leishmania infection among Thai patients with HIV, a high risk group, is still lacking. Thus, the objectives of this study were to determine the prevalence and the risk factors associated with Leishmania infection among patients with HIV/AIDS in Trang province, southern Thailand.


Study design and population

A cross-sectional study of Leishmania infection was conducted between February 2015 and February 2016. Eligible participants were >18 years old and attending an HIV clinic, Trang Hospital, Trang province. They visited the clinic every 6 months for follow-up testing and to receive antiretroviral therapy (ART). They lived in ten districts of Trang province, other nine provinces located in the south, and three provinces in other regions of Thailand. Clinical information of participants was collected from patients’ medical records.

Ethics statement

Written informed consent was obtained from all participants. All participants were >18 years old. All analyzed data were anonymized. The research protocol was approved by the Ethics Committee of the Royal Thai Army Medical Department and the Ethics Committee of Mahidol University, Thailand.

Blood collection

Eight milliliters of EDTA anti-coagulated blood samples were collected. The whole blood was centrifuged at 900 × g for 10 minutes to separate the plasma and buffy coat and was then kept at −20°C until further use.


Seropositivity of Leishmania infection was defined as detection of antibodies in individuals who were exposed to Leishmania infection and being either symptomatic or asymptomatic.

Asymptomatic Leishmania infection was defined as individuals who experienced no symptoms of VL but presented a positive test by DAT or PCR assays.

Symptomatic VL was defined as individuals having a history of fever lasting at least 2 weeks with splenomegaly. One or more of the following clinical characteristics may be observed: hepatomegaly, weight loss, anemia, leucopenia, thrombocytopenia, and hypergammaglobulinemia. Detection of the parasites must be confirmed under microscopic examination or by PCR assay using any clinical samples (e.g., bone marrow aspirates, lymph node, blood, and/or other biopsy samples).

Detection of Leishmania antibodies

Leishmania antibodies were assayed using the commercial DAT kit (Biomedical Research) according to the manufacturer’s instruction. The positive plasma control was obtained from confirmed VL cases using the PCR method. For the negative control, plasma from healthy individuals was used. The cutoff value of positive DAT titers was ≥1:100 following manufacturer recommendation.

Leishmania DNA detection

DNA was extracted from 200 μL of buffy coat sample using Gen UP gDNA Kit (Biotech). Nested PCR was used to amplify the ITS1 region of the ribosomal DNA (rDNA) gene of Leishmania. In the primary PCR, primers LITSR and L5.8S were used to amplify the 319–348 amplicons [8]. The newly designed secondary primers LITSR2 (CTG-GAT-CAT-TTT-CCG-ATG-ATT) and L5.8S inner (GTT-ATG-TGA-GCC-GTT-ATC-C) generated 230–280 amplicons depending on Leishmania species. PCR reactions were performed using the MJ Mini thermal cycler (BioRad) in volumes of 25 μL, containing 12.5 pmol of each primer, 0.2 mM dNTP, 1.5 mM MgCl2, 1× PCR buffer, 1 U of Taq DNA polymerase, and 4 μL of DNA template. DNA of L. martiniquensis promastigotes (MHOM/MQ/92/MAR1) was used as the positive control. The condition was started by pre-denaturation at 94°C for 3 minutes followed by 35 cycles: denaturation at 94°C for 1 minute, annealing temperature at 54°C for 30 seconds, and extension at 72°C for 30 seconds. Final extension was at 72°C for 5 minutes. PCR products were separated by electrophoresis in 1.5% agarose gel stained with SYBR Safe (Invitrogen). The results were visualized and documented by Molecular Imager Gel Doc XR+ System with Imager Lab 3.0 (BioRad).

Sequence analysis

Positive PCR products were sent to U2Bio Co. Ltd., South Korea for sequencing. Chromatograms were validated using BioEdit version 7.0.1. The sequences were multiple-aligned with reference Leishmania strains retrieved from GenBank. The phylogenetic tree was constructed by using the neighbor-joining (NJ) method using the MEGA program, version 7.0. The reliability was tested by 1,000 bootstrap replications and Tajima-Nei was selected for the DNA substitution model of phylogenetic analysis.


To determine the risk factors and outcomes of Leishmania infection, standardized questionnaires were used. Enrolled subjects with HIV were interviewed face-to-face covering demographic data, socioeconomic status, clinical symptoms, and associated risk behaviors.

Statistical analysis

The association between potential risk factors and Leishmania infection was assessed by univariate and multivariate logistic regression analysis. Odds ratios and 95% confidence intervals (CI) were calculated and p values <0.05 were considered statistically significant. All analyses were performed using STATA, version SE14 (Stata Corporation, College Station, TX, USA). (


Study population and characteristics

A total of 724 participants with HIV were enrolled in this study. Of these, 643 (88.8%) filled out questionnaires. Living areas of participants were as follows: 570 (88.6%) lived in Trang province and 67 (10.4%) lived in nine other provinces located in the south. Only five (0.8%) were from other regions of Thailand. The mean age was 43.6 ± 8.5 years. The characteristics of the enrolled subjects are shown in Table 1. Their clinical characteristics and risk behaviors during the past one year (S1 and S2 Tables) included history of injection drug users (IDUs) at 16.6% while 13.7% were non-injection drug users (NIDUs). A total of 512 (79.6%) subjects lived in non-stilt houses while 131 (20.4%) lived in stilt houses. Most of the individuals (68.6%) had CD4+ levels more than 500 cells/μL and only 9.6% were less than 200 cells/μL.

Table 1
The characteristics of enrolled patients with HIV (n = 643) were analyzed using three categories: i) patients who were either seropositive by DAT analysis with titers of >100 or positive by PCR assay; ii) patients who were seropositive by DAT ...

Prevalence of Leishmania infection

Three categories of data analysis of the prevalence of Leishmania infection were performed (Table 1). The first group comprised patients who were either seropositive by DAT analysis with titers of ≥100 or positive by PCR assay. The second category involved patients who were seropositive by DAT analysis with titers of ≥100, and the last group comprised patients who were positive only by PCR assay. The prevalence of Leishmania infection using positive results either by DAT or PCR assays was 25.1% (182/724). Seropositive cases comprised 18.5% (134/724) while Leishmania DNA detection by PCR was 8.4% (61/724). Only 1.8% (13/724) were positive using both methods (Table 2). Table 3 shows numbers of positive Leishmania infection using DAT, the titer of DAT, and PCR. Thus, the overall prevalence of asymptomatic Leishmania infection was 24.9 (180/724). Tables Tables44 and and55 show affected areas of Leishmania infections.

Table 2
Number of Leishmania infections by DAT analysis with titers of ≥100 or by PCR assay (n = 724).
Table 3
Numbers of positive DAT analysis together with numbers of positive PCR results at the titers of 1:100 to 1:6400.
Table 4
Patients who were either seropositive by DAT analysis with titers of ≥100 or positive by PCR assay among 13 provinces.
Table 5
Patients who were either seropositive by DAT analysis with titers of ≥100 or positive by PCR assay among 10 districts of Trang province.

Regarding the analysis of those who were either seropositive by DAT or positive by PCR assay, the prevalence of Leishmania infection significantly differed among participants who lived in stilt houses compared to those living in non-stilt houses (p = 0.03), those who developed jaundice (p = 0.02), having opportunistic infection (p = 0.002) especially tuberculosis (p = 0.001), and those having low CD4+ levels <500 cells/μL (p = 0.003) (Fig 1). No significant difference was found among age group, gender, educational level, occupation, working outdoors at night, average income, years of HIV diagnosis, viral load, history of going abroad, drug use (IDUs/NIDUs), pet/animal owner, animal shed nearby the house, plantation nearby the house, and bed net use (S1 and S2 Tables).

Fig 1
The percent of positive Leishmania infections by clinical characteristics and risk behaviors of enrolled patients with HIV (n = 643) was characterized into three categories: i) positivity diagnosed by DAT analysis with titers of ≥100 or by PCR ...

Symptomatic VL

Two patients showed symptomatic VL. The first case was a 39-year-old herdsman living in Phuket province. He had a history of NIDU. He presented with a fever for more than two weeks. Many nodules were observed on the trunk. Laboratory findings revealed CD4+ levels at 173 cells/μL with an undetectable viral load. The DAT titer was 1:6400. The causative agent was L. martiniquensis, which was identified by the nested ITS1-PCR using the buffy coat and skin biopsy. The patient was treated with amphotericin B. However, he died of disease progression one year after initial VL diagnosis. The second case involved a 41-year-old male who originally lived in Trang province. He worked on a rubber plantation. He was both an IDU and NIDU. He developed epistaxis and bleeding gums. Laboratory findings revealed pancytopenia and CD4+ levels of 622 cells/μL with an undetectable viral load. Intracellular amastigotes were observed from the lymph node biopsy with Giemsa stain. The DAT titer was 1:6400 and the causative agent was L. martiniquensis, which was identified using the nested ITS1-PCR of the buffy coat. The patient was treated with amphotericin B but he died within one week after treatment.

Risk factor analysis

Univariate and multivariate analysis of risk factors for acquiring Leishmania infections by DAT titers of ≥100 or positive by PCR assays are shown in Table 6. Univariate analysis showed that participants who lived in stilt houses had higher risk (OR = 1.58, 95% CI = 1.04–2.39) of contracting Leishmania when compared with those living in non-stilt houses. CD4+ levels between 200 and 500 cells/μL were at higher risk than those who had CD4+ levels >500 cells/μL (OR = 1.70, 95% CI = 1.29–2.97). After adjusting for age, gender, NIDUs, history of traveling abroad, pet owners and raising animals in housing areas, bed net use, animal shed and plantation nearby the house, underlying diseases, viral load, and duration of HIV diagnosis, multivariate logistic regression analysis revealed that those living in stilt houses had greater risk (OR = 1.60, 95% CI = 1.04–2.47) of acquiring the infection when compared with those living in non-stilt houses. In addition, those who had CD4+ levels 200–500 cells/μL (OR = 2.13, 95% CI = 1.36–3.32) and <200 cells/μL (OR = 1.98, 95% CI = 1.06–3.73) had higher risk of contracting Leishmania than those who had CD4+ levels >500 cells/μL. In addition, multivariate analysis of the associated risk factors of Leishmania infection using only seropositivity showed that participants who were NIDU had higher risk of presenting Leishmania seropositivity than those who were not (OR = 2.23, 95% CI = 1.27–3.92). In addition, those who had CD4+ levels 200–500 cells/μL were also at higher risk of being seropositive than those who had CD4+ levels >500 cells/μL (OR = 2.09, 95% CI = 1.27–3.44) (S3 Table). Multivariate analysis of the associated risk factors of Leishmania infection using only positive PCR results showed that those who had a detectable viral load >50 copies/mL were at higher risk of acquiring detectable Leishmania DNA in the blood (OR = 2.31, 95% CI = 1.01–5.29) than those who had an undetectable viral load after adjusting for those variables as mentioned above (S4 Table).

Table 6
Univariate and multivariate analysis of associated risk factors of Leishmania infection using seropositive results by DAT analysis with titers of ≥100 or positive by PCR assay.

Phylogenetic analysis and Leishmania species identification

Of 61 samples, nucleotide sequencing was successful for 49. These sequences, together with 16 reference sequences of different Leishmania species, were included to construct the phylogenetic tree using the NJ method (Fig 2). The phylogenetic analyses grouped the sequences into five separated clades. The majority of the samples (20 (40.8%) and 13 (26.5%)) were closely related to L. siamensis and L. martiniquensis, respectively, whereas ten (20.4%) were closely related to L. donovani complex, five (10.2%) were related to L. lainsoni, and one (2.1%) was related to L. major.

Fig 2
The unrooted phylogenetic tree inferred from the ITS1 sequences using the NJ method.


This was the first study providing important information of the prevalence of co-infection of Leishmania among Thai patients with HIV who had been regularly attending the HIV clinic in Trang province. The prevalence of Leishmania infection was approximately one-fourth of the 724 participants determined by either DAT or PCR assays. In the past, Leishmania infections were previously reported in five provinces in the south (Surat Thani [6], Phang-Nga [9], Trang [10], Songkhla [11], and Satun [12]) where most people mainly earn their living in agricultural sectors. The climate and humidity in the south are suitable for the sand fly’s habitat where potential sand fly vectors have been reported [13, 14]. Our study revealed four new affected areas of Leishmania infections in the south (Phuket, Krabi, Nakhon Si Thammarat, and Phatthalung provinces). Thus, at present, overall affected areas cover nine southern provinces.

For VL/HIV co-infection, the sensitivities of DAT to detect antibodies against Leishmania infection were 50–84% [15]. Related studies of VL in immunocompetent subjects in endemic areas in India and Iran used DAT titers at different cutoff values ranging from 1:800 to 1:3200 [1618]. A cutoff value at a titer of 1:100 was previously used to screen VL among HIV-positive patients who developed clinical symptoms in northeast Iran [19]. A low DAT titer of 1:200 was also detected in an immunocompetent VL Thai patient caused by L. martiniquensis [9]. In this study, one patient having a DAT titer of 1:100 also produced a positive PCR result. Therefore, a positive serological test at low titers could have diagnostic value to detect the infection. Thus, the cutoff values of DAT varied when conducted in different study populations as well as areas of study where cross-reactivity of DAT against other blood parasite infections could have occurred. A systematic review revealed that DAT titers detected in symptomatic patients were higher than those of asymptomatic patients [20]. Patients presenting very high DAT titer would have significantly greater disease progression than those presenting low titers [21]. In this study, clinical characteristics of VL were observed in two symptomatic cases that had DAT titers of 1:6400 together with positive PCR results. Thus, a close follow-up for those asymptomatic infections showing positive results of DAT or PCR is needed.

In this study, serological and molecular diagnosis among individuals with HIV was not in concordance with other studies [22]. Our results showed that PCR positivity was low when compared with numbers of DAT positivity. The potential reasons among DAT positive individuals who might become PCR negative could include degradation and clearance of Leishmania DNA after infection, which corresponds to development of protective immunity due to the use of antiviral drugs for HIV. A positive PCR test among DAT negative individuals could occur when the individual was bitten by a Leishmania infected sand fly, but either immunity has not yet developed or antibody levels are too low to be detected by the methods employed, especially in HIV-positive populations [23].

Before using HAART, asymptomatic infection in HIV-Leishmania co-infection in Europe was 4–33% [24]. However, the incidence has been reduced to 20% after ART drugs had been given to all individuals with HIV [3]. In this study, all patients with HIV regularly received HAART treatment that could restore TH1 cytokine and antibody production [25]. HAART-treated HIV patients demonstrated a better ability to control Leishmania infection [26]. Many patients with subclinical VL did not develop clinical symptoms after taking HAART medications while some developed the disease [3]. Other risk factors (e.g., stage of HIV infection, parasite virulence, drug resistance, nutritional status, age, and gender) may be involved in disease progression [1, 27]. Prospective studies are needed to determine disease progression in this population.

In this study, most enrolled participants (88.6%) were not randomly selected and originally lived in Trang province. Thus, our results do not represent the prevalence of co-infection of Leishmania and patients with HIV of the country, and they do not represent the prevalence in each district of Trang province. However, this study showed the magnitude of the seroprevalence and significant numbers of subclinical results of Leishmania DNA detection circulating in the blood in Thai individuals with HIV.

In southern Europe, IDUs were the most important risk factor accounting for more than 90% of all cases [28]. Our results showed no significant difference in prevalence among IDUs, while seroprevalence was associated with NIDUs. Leishmaniasis was also associated with socioeconomic status. In India, housing materials such as mud, plants, and earthen floors were risk factors for VL [20, 29]. Our study showed that living in stilt houses was an independent associated risk factor for Leishmania infection. The presence of stilts provided an open area under the house that might have increased chances of sand fly bites to humans as well as providing resting sites for sand flies. Sand flies frequently bite at dusk [30] when most people spend their time at the open area of the house. In addition, CD4+ levels play an important role to protect the host from opportunistic infections. Related studies showed that the first episode of symptomatic VL diagnosis involved more than 80% of patients with HIV who had low CD4+ levels [26, 31]. Our results also confirmed that low CD4+ levels (<200 cells/μL) as well as 200 to 500 cells/μL significantly increased the risk of Leishmania infection.

In this study, Leishmania DNA detection by PCR assay was associated with detectable viral load. No correlation has been reported between PCR positivity and CD4+ levels, whereas the correlation between HIV viral load and parasitemia was observed among asymptomatic patients [22]. Clinical progression of HIV/AIDS was simultaneously promoted by VL. HIV infection enhances parasite growth by modulating significant cytokine response to Leishmania while the parasite upregulates viral expression [3].

Public health awareness of Leishmania infection in Thailand started when one autochthonous VL was reported in 1996. However, at that time, the disease was uncommon as well as unfamiliar to most physicians, which could have led to a lot of underreporting of leishmaniasis cases in the past 20 years. Using the molecular method, this cross-sectional study was the first to systematically estimate the prevalence of Leishmania infection among patients with HIV, a high risk group, revealing not only L. siamensis and L. martiniquensis infection but also infections of other species (e.g., L. donovani complex, L. lainsoni, and L. major) that already existed in the affected area. From the phylogenetic analysis, the ITS1 region is one of the more powerful targets used to discriminate the Leishmania species [32]. Our previous study showed that the ITS1 region had the highest sensitivity to detect L. martiniquensis and L. siamensis compared with the other genes: hsp70, cyt b, and SSU-rRNA [5]. Additionally, we used hsp70-PCR and kDNA-PCR to amplify DNA samples of L. major, L. donovani, L. lainsoni, L. siamensis, and L. martiniquenensis. Unfortunately, negative PCR results were obtained due to lower sensitivities of PCR amplification using these genes.

L. martiniquensis and L. siamensis were the predominant species detected in this study. Additionally, L. infantum (which causes VL) was previously reported in an HIV negative individual living in Bangkok [33]. The L. donovani complex species are the major causative agents of VL worldwide. The distribution of the complex species might have been introduced by travelers or workers from VL-endemic areas to Thailand [34]. L. major, a causative agent for CL in the Old World, causes zoonotic transmission, especially Afghanistan and India [31]. VL, caused by L. major, has occasionally been reported among patients with HIV [3]. L. donovani complex species and L. major have been reported in China, Bangladesh, India, and Nepal [15]. Thus, distribution of these two species in Thailand could be possible. L. lainsoni infection causes localized CL and was reported in South America in Bolivia, Peru, Suriname, French Guiana, and Brazil [15, 35]. This is the first report of L. major and L. lainsoni infection among individuals with HIV in Thailand. In addition, no report of L. lainsoni has been documented in the Old World, especially among patients with HIV.

To prevent and control VL, understanding disease epidemiology is extremely important. A cohort study conducted in this population as well as studies of potential vectors and animal reservoirs are needed. Moreover, large scale molecular epidemiological studies in other high morbidity areas are required for this emerging disease.

Supporting information

S1 Checklist

STROBE checklist for cross-sectional studies.


S1 Table

Clinical characteristics in the past year of enrolled patients with HIV (n = 643) were analyzed using three categories: i) patients who were either seropositive by DAT analysis with titers of >100 or positive by PCR assay, ii) patients who were seropositive by DAT analysis with titers of >100, and iii) patients who were positive only by PCR assay.


S2 Table

Characteristics of risk behaviors of enrolled patients with HIV (n = 643) were analyzed using three categories: i) patients who were either seropositive by DAT analysis with titers of >100 or positive by PCR assay, ii) patients who were seropositive by DAT analysis with titers of >100, and iii) patients who were positive only by PCR assay.


S3 Table

Univariate and multivariate analysis of associated risk factors of Leishmania infection using seropositive results by DAT analysis with titers of >100.


S4 Table

Univariate and multivariate analysis of associated risk factors of Leishmania infection using seropositive results by PCR assay.



We express our appreciation to all participants for their enthusiastic support of this study.

Funding Statement

This work was supported by the National Science and Technology Development Agency, Science and Achievement Scholarship of Thailand, and the Phramongkutklao College of Medicine Research Fund. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Data Availability

Data Availability

The datasets supporting the finding of this article are included in Fig 2. The data are available through the NCBI database (GenBank) under the accession no. KY982629-KY982677.


1. Lindoso JAL, Cunha MA, Queiroz IT, Moreira CHV. Leishmaniasis–HIV coinfection: current challenges. HIV AIDS (Auckl). 2016;8:147–56. [PMC free article] [PubMed]
2. Maia Z, Lirio M, Mistro S, Mendes CM, Mehta SR, Badaro R. Comparative study of rK39 Leishmania antigen for serodiagnosis of visceral leishmaniasis: systematic review with meta-analysis. PLoS Negl Trop Dis. 2012;6(1):e1484 doi: 10.1371/journal.pntd.0001484 [PMC free article] [PubMed]
3. Alvar J, Aparicio P, Aseffa A, Den Boer M, Canavate C, Dedet JP, et al. The relationship between leishmaniasis and AIDS: the second 10 years. Clin Microbiol Rev. 2008;21(2):334–59, table of contents. doi: 10.1128/CMR.00061-07 [PMC free article] [PubMed]
4. de Ruiter CM, van der Veer C, Leeflang MM, Deborggraeve S, Lucas C, Adams ER. Molecular tools for diagnosis of visceral leishmaniasis: systematic review and meta-analysis of diagnostic test accuracy. J Clin Microbiol. 2014;52(9):3147–55. doi: 10.1128/JCM.00372-14 [PMC free article] [PubMed]
5. Hitakarun A, Tan-ariya P, Siripattanapipong S, Mungthin M, Piyaraj P, Naaglor T, et al. Comparison of PCR methods for detection of Leishmania siamensis infection. Parasit Vectors. 2014;7. [PMC free article] [PubMed]
6. Thisyakorn U, Jongwutiwes S, Vanichsetakul P, Lertsapcharoen P. Visceral leishmaniasis: the first indigenous case report in Thailand. Trans R Soc Trop Med Hyg. 1999;93(1):23–4. [PubMed]
7. Leelayoova S, Siripattanapipong S, Manomat J, Piyaraj P, Tan-Ariya P, Bualert L, et al. Leishmaniasis in Thailand: A Review of Causative Agents and Situations. Am J Trop Med Hyg. 2017;96(3):534–42. doi: 10.4269/ajtmh.16-0604 [PubMed]
8. El Tai NO, El Fari M, Mauricio I, Miles MA, Oskam L, El Safi SH, et al. Leishmania donovani: intraspecific polymorphisms of Sudanese isolates revealed by PCR-based analyses and DNA sequencing. Exp Parasitol. 2001;97(1):35–44. doi: 10.1006/expr.2001.4592 [PubMed]
9. Sukmee T, Siripattanapipong S, Mungthin M, Worapong J, Rangsin R, Samung Y, et al. A suspected new species of Leishmania, the causative agent of visceral leishmaniasis in a Thai patient. Int J Parasitol. 2008;38(6):617–22. doi: 10.1016/j.ijpara.2007.12.003 [PubMed]
10. Bualert L, Charungkiattikul W, Thongsuksai P, Mungthin M, Siripattanapipong S, Khositnithikul R, et al. Autochthonous disseminated dermal and visceral leishmaniasis in an AIDS patient, southern Thailand, caused by Leishmania siamensis. Am J Trop Med Hyg. 2012;86(5):821–4. doi: 10.4269/ajtmh.2012.11-0707 [PMC free article] [PubMed]
11. Chusri S, Hortiwakul T, Silpapojakul K, Siriyasatien P. Consecutive cutaneous and visceral leishmaniasis manifestations involving a novel Leishmania species in two HIV patients in Thailand. Am J Trop Med Hyg. 2012;87(1):76–80. doi: 10.4269/ajtmh.2012.11-0749 [PMC free article] [PubMed]
12. Osatakul S, Mungthin M, Siripattanapipong S, Hitakarun A, Kositnitikul R, Naaglor T, et al. Recurrences of visceral leishmaniasis caused by Leishmania siamensis after treatment with amphotericin B in a seronegative child. Am J Trop Med Hyg. 2014;90(1):40–2. doi: 10.4269/ajtmh.13-0398 [PMC free article] [PubMed]
13. Kanjanopas K, Siripattanapipong S, Ninsaeng U, Hitakarun A, Jitkaew S, Kaewtaphaya P, et al. Sergentomyia (Neophlebotomus) gemmea, a potential vector of Leishmania siamensis in southern Thailand. BMC Infect Dis. 2013;13:333 doi: 10.1186/1471-2334-13-333 [PMC free article] [PubMed]
14. Chusri S, Thammapalo S, Chusri S, Thammapalo S, Silpapojakul K, Siriyasatien P. Animal reservoirs and potential vectors of Leishmania siamensis in southern Thailand. Southeast Asian J Trop Med Public Health. 2014;45(1):13–9. [PubMed]
15. Aronson N, Herwaldt BL, Libman M, Pearson R, Lopez-Velez R, Weina P, et al. Diagnosis and Treatment of Leishmaniasis: Clinical Practice Guidelines by the Infectious Diseases Society of America (IDSA) and the American Society of Tropical Medicine and Hygiene (ASTMH). Am J Trop Med Hyg. 2017;96(1):24–45. doi: 10.4269/ajtmh.16-84256 [PMC free article] [PubMed]
16. Topno RK, Das VNR, Ranjan A, Pandey K, Singh D, Kumar N, et al. Asymptomatic Infection with Visceral Leishmaniasis in a Disease-Endemic Area in Bihar, India. Am J Trop Med Hyg. 2010;83(3):502–6. doi: 10.4269/ajtmh.2010.09-0345 [PMC free article] [PubMed]
17. Ostyn B, Gidwani K, Khanal B, Picado A, Chappuis F, Singh SP, et al. Incidence of symptomatic and asymptomatic Leishmania donovani infections in high-endemic foci in India and Nepal: a prospective study. PLoS Negl Trop Dis. 2011;5(10):e1284 doi: 10.1371/journal.pntd.0001284 [PMC free article] [PubMed]
18. Mohebali M. Visceral leishmaniasis in Iran: Review of the Epidemiological and Clinical Features. Iran J Parasitol. 2013;8(3):348–58. [PMC free article] [PubMed]
19. Shafiei R, Mohebali M, Akhoundi B, Galian MS, Kalantar F, Ashkan S, et al. Emergence of co-infection of visceral leishmaniasis in HIV-positive patients in northeast Iran: a preliminary study. Travel Med Infect Dis. 2014;12(2):173–8. doi: 10.1016/j.tmaid.2013.09.001 [PubMed]
20. Hirve S, Boelaert M, Matlashewski G, Mondal D, Arana B, Kroeger A, et al. Transmission Dynamics of Visceral Leishmaniasis in the Indian Subcontinent—A Systematic Literature Review. PLoS Negl Trop Dis. 2016;10(8):e0004896 doi: 10.1371/journal.pntd.0004896 [PMC free article] [PubMed]
21. Hasker E, Malaviya P, Gidwani K, Picado A, Ostyn B, Kansal S, et al. Strong association between serological status and probability of progression to clinical visceral leishmaniasis in prospective cohort studies in India and Nepal. PLoS Negl Trop Dis. 2014;8(1):e2657 doi: 10.1371/journal.pntd.0002657 [PMC free article] [PubMed]
22. Colomba C, Saporito L, Vitale F, Reale S, Vitale G, Casuccio A, et al. Cryptic Leishmania infantum infection in Italian HIV infected patients. BMC Infect Dis. 2009;9:199 doi: 10.1186/1471-2334-9-199 [PMC free article] [PubMed]
23. Sudarshan M, Singh T, Singh AK, Chourasia A, Singh B, et al. Quantitative PCR in Epidemiology for Early Detection of Visceral Leishmaniasis Cases in India. PLoS Negl Trop Dis. 2014;8(12):e3366 doi: 10.1371/journal.pntd.0003366 [PMC free article] [PubMed]
24. Alvar J, Gutierrez-Solar B, Molina R, Lopez-Velez R, Garcia-Camacho A, Martinez P, et al. Prevalence of Leishmania infection among AIDS patients. Lancet. 1992;339(8806):1427. [PubMed]
25. Orsini M, Canela JR, Disch J, Maciel F, Greco D, Toledo A Jr., et al. High frequency of asymptomatic Leishmania spp. infection among HIV-infected patients living in endemic areas for visceral leishmaniasis in Brazil. Trans R Soc Trop Med Hyg. 2012;106(5):283–8. doi: 10.1016/j.trstmh.2012.01.008 [PubMed]
26. Cruz I, Nieto J, Moreno J, Canavate C, Desjeux P, Alvar J. Leishmania/HIV co-infections in the second decade. Indian J Med Res. 2006;123(3):357–88. [PubMed]
27. van Griensven J, Carrillo E, Lopez-Velez R, Lynen L, Moreno J. Leishmaniasis in immunosuppressed individuals. Clin Microbiol Infect. 2014;20(4):286–99. doi: 10.1111/1469-0691.12556 [PubMed]
28. Herrador Z, Gherasim A, Jimenez BC, Granados M, San Martin JV, Aparicio P. Epidemiological changes in leishmaniasis in Spain according to hospitalization-based records, 1997–2011: raising awareness towards leishmaniasis in non-HIV patients. PLoS Negl Trop Dis. 2015;9(3):e0003594 doi: 10.1371/journal.pntd.0003594 [PMC free article] [PubMed]
29. Kumar R, Kumar P, Chowdhary RK, Pai K, Mishra CP, Kumar K, et al. Kala-azar epidemic in Varanasi district, India. Bull World Health Organ. 1999;77(5):371–4. [PubMed]
30. Maroli M, Feliciangeli MD, Bichaud L, Charrel RN, Gradoni L. Phlebotomine sandflies and the spreading of leishmaniases and other diseases of public health concern. Med Vet Entomol. 2013;27(2):123–47. doi: 10.1111/j.1365-2915.2012.01034.x [PubMed]
31. World Health Organization. Control of the leishmaniases: report of a meeting of the WHO Expert Commitee on the Control of Leishmaniases, Geneva, 22–26 March 2010: Geneva: World Health Organization; 2010. 199 p.
32. Leelayoova S, Siripattanapipong S, Hitakarun A, Kato H, Tan-ariya P, Siriyasatien P, et al. Multilocus characterization and phylogenetic analysis of Leishmania siamensis isolated from autochthonous visceral leishmaniasis cases, southern Thailand. BMC Microbiol. 2013;13:60 doi: 10.1186/1471-2180-13-60 [PMC free article] [PubMed]
33. Maharom P, Siripattanapipong S, Mungthin M, Naaglor T, Sukkawee R, Pudkorn R, et al. Visceral leishmaniasis caused by Leishmania infantum in Thailand. Southeast Asian J Trop Med Public Health. 2008;39(6):988–90. [PubMed]
34. Viriyavejakul P, Viravan C, Riganti M, Punpoowong B. Imported cutaneous leishmaniasis in Thailand. Southeast Asian J Trop Med Public Health. 1997;28(3):558–62. [PubMed]
35. Correa JR, Brazil RP, Soares MJ. Leishmania (Viannia) lainsoni (Kinetoplastida: Trypanosomatidae), a divergent Leishmania of the Viannia subgenus—a mini review. Mem Inst Oswaldo Cruz. 2005;100(6):587–92. doi: S0074-02762005000600014 [PubMed]

Articles from PLoS Neglected Tropical Diseases are provided here courtesy of Public Library of Science