Skip to main content
Advertisement
  • Loading metrics

High levels of infectiousness of asymptomatic Leishmania (Viannia) braziliensis infections in wild rodents highlights their importance in the epidemiology of American Tegumentary Leishmaniasis in Brazil

Abstract

Background

The epidemiological significance of wildlife infections with aetiological agents causing human infectious diseases is largely determined by their infection status, contact potential with humans (via vectors for vector-borne diseases), and their infectiousness to maintain onward transmission. This study quantified these parameters in wild and synanthropic naturally infected rodent populations in an endemic region of tegumentary leishmaniasis in northeast Brazil.

Methods

Capture-mark-recapture (CMR) of rodents was conducted over 27 months in domestic/peri domestic environs, household plantations and nearby Atlantic Forest (9,920 single trap nights). Rodent clinical samples (blood and ear tissue) were tested for infection by conventional PCR and quantitative PCR (qPCR) for Leishmania (Viannia) braziliensis, and xenodiagnosis to measure infectiousness to the local sand fly vector.

Results

A total 603 individuals of 8 rodent species were (re)captured on 1,051 occasions. The most abundant species were Nectomys squamipes (245 individuals, 41% of the total catch), Rattus rattus (148, 25%), and Necromys lasiurus (83, 14%). All species were captured in greater relative frequencies in plantations; R. rattus was the only species captured in all three habitats including in and around houses. Four species, comprising 22.6% of individuals captured at least twice, were geolocated in more than one habitat type; 78.6% were infected with L. (V.) braziliensis, facilitating inter-species and inter-habitat transmission. Species specific period prevalence ranged between 0%-62% being significantly higher in N. squamipes (54–62%) and Hollochillus sciureus (43–47%). Xenodiagnosis was performed on 41 occasions exposing 1,879 Nyssomyia whitmani sand flies to five rodent species (37 individuals). Similar mean levels of infectiousness amongst the more common rodent species were observed. Longitudinal xenodiagnosis of the N. squamipes population revealed a persistent level of infectiousness over 13 months follow-up, infecting a median 48% (IQR: 30.1%-64.2%) of exposed blood-fed vectors. The proportion of exposed flies infected was greater in the low compared to in the high seasonal period of vector abundance. L. (V.) braziliensis parasite loads in rodent blood quantified by qPCR were similar across rodent species but did not represent a reliable quantitative marker of infectiousness to sand flies. The standardised risk of rodent infection in plantations was 70.3% relative to 11.3% and 18.4% in peri domestic and forest habitats respectively. R. rattus was the only exception to this trend indicating greatest risk in the peri domestic environment.

Conclusions

The results support the view that a collective assemblage of wild and synanthropic rodent species is an important wild reservoir of L. (V.) braziliensis in this region, with N. squamipes and R. rattus probably playing a key role in transmission within and between habitat types and rodent species. Rodents, and by implication humans, are at risk of infection in all sampled habitats, but more so in homestead plantations. These conclusions are based on one of the longest CMR study of small rodents in an American Tegumentary Leishmaniasis (ATL) foci.

Author summary

Our previous studies on the prevalence of leishmanial infections in wild caught rodents suggested that they represent potential reservoirs of Leishmania (Viannia) braziliensis in northeast Brazil, supported by showing that laboratory bred rodents that were experimentally infected with L. (V.) braziliensis were infectious to Lutzomyia longipalpis. The present field study reports the abundance of rodent species in different ecotopes, the prevalence of Leishmania infections in a cohort of naturally infected rodent species and the transmission potential of the most abundant species by xenodiagnoses using the primary vector, Nyssomyia whitmani. The results incriminate the local guild of rodents as a collective wildlife reservoir in the zoonotic transmission cycle of L. (V.) braziliensis. The greatest infection risk to rodents, and likely also to humans, is in peri domestic plantations, where N. squamipes, Ne. lasiurus and R. rattus, in particular, occur at high densities; the former species supported the highest infection prevalences and were infectious for extended periods. R. rattus was the only species captured supporting infections in and around households in addition to the other habitats. More than 20% of rodents captured at least twice moved between habitat types. There was a qualitative but not quantitative relationship between L. (V.) braziliensis parasite loads in rodent blood samples and xenopositive sand flies.

Introduction

The leishmaniases are a collection of protozoal diseases of humans and animals associated with phylogenetically diverse parasites that belong to the Leishmaniinae subfamily [1]. The genus Leishmania comprises of more than 20 species, transmitted between mammalian hosts by principally blood-seeking female Phlebotomine sand flies, and most involving one or more non-human animal reservoir in their zoonotic transmission cycle. Depending on the infecting Leishmania species, the disease in humans causes cutaneous, mucocutaneous, or visceral clinical manifestations ranging from self-limiting dermatological lesions to life-threatening visceral involvement requiring complex diagnosis and treatment regimes. Currently endemic in up to 92 countries worldwide, prevention or curtailing of transmission and geographical expansion [2,3] requires a full understanding of the aetiological agents, arthropod vectors and animal hosts maintaining transmission [4].

One such challenge is the complex eco-epidemiology of Leishmania (Viannia) braziliensis that causes American Tegumentary Leishmaniasis (ATL). ATL is a human disease of skin and/or mucosal tissues caused by a number of Leishmania parasites throughout the Americas [57]. Subclinical L. (V.) braziliensis infections have been detected in a multitude of wild, synanthropic, and domesticated animal hosts, but their epidemiological significance are unknown. The uncertainties emanate in part from misinterpretation of host infection data which alone are insufficient to differentiate reservoir hosts (sources of Leishmania) from non-reservoir hosts (Leishmania sink hosts) as infection per se does not necessarily equate to their transmission potential to sand fly vectors [810].

Advances in our understanding are further hindered by the lack of experimental studies due to difficulties to establish L. (V.) braziliensis in culture [11], to maintain detectable infections in laboratory animal models [6,12], and standardisation issues associated with the diversity of genetic and molecular profiles of the parasite, coupled with the non-uniformity in the clinical behaviour in laboratory animals [13,14]. Consequently, compared to the other Leishmania species, there are few experimental studies of L. (V.) braziliensis. Moreover, in the absence of longitudinal studies of naturally infected host populations, it is not clear if the course of laboratory infections indeed mimic those in nature.

Almost nothing is known about wildlife host infectiousness to sand fly vectors, and how it relates to host parasite loads and tissue tropism, and if quantitative PCR (qPCR) methods to detect Leishmania genome equivalents in host tissues offers a reliable marker of host infectiousness by which to measure transmission potential at the population-level; xenodiagnosis, the standard, is difficult to perform at scale.

Here, we present the results of a longitudinal capture-mark-recapture (CMR) study of wild and synanthropic rodents naturally infected with L. (V.) braziliensis in northeast Brazil. The specific objectives were to (i) document rodent species abundance in different ecotypes associated with human activity; (ii) quantify their infection prevalence; (iii) measure their infectiousness to the sand fly vector; (iv) determine the relative infection risk by ecotope; (v) investigate the potential of tissue parasite loads as a surrogate marker of host infectiousness.

Methods

Ethical considerations

Study protocols were approved by the Ethics Committee on Animal Use (CEUA / IAM) (No. 017/2011), and permission from IBAMA (Brazilian Institute of Environment and Renewable Natural Resources) through the Authorization for Activities with Scientific Purpose (No. 12749).

Study area

Fieldwork was undertaken in two rural endemic communities, Engenho Raiz de Dentro (8° 42.7944’S, 35° 44.7078’W) and Engenho Refrigério (8° 42.7963’S, 35° 49.4742’W) located 8.73 km apart, in Amaraji municipality in the State of Pernambuco, north-east Brazil (S1 Fig). The municipality is located south of Zona da Mata about 100 km from Recife, 385m above sea level, and belongs to the north-eastern coastal region known as the Coastal Plateaus (Tabuleiros). The municipality population is approximately 22,829 habitants of which 27% live in rural areas [15,16].

Local temperatures range from 18°C to 29°C, with an average annual rainfall of at least 1,000mm, and a high degree of humidity throughout the year (Köppen–Geiger climate classification: Aw—Tropical wet and dry). The municipality of Amaraji is situated in the domains of the Serinhaém River Basin fed by small streams and water sources. Originally the land was divided into extensive sugar cane plantations known as Engenhos (Mills), of which some were expropriated by the government and divided into plots (˜Parcels˜) to be occupied by local farmers who became their owners. This led to an increase in the rural population living in scattered small homesteads (S2 and S3 Figs) and attending their plantations principally of bananas, pineapples and vegetable cash crops. Pockets of sub-deciduous Atlantic forest typically occur on the steep-sided hill tops of oxisols and podzolic soils, while the alluvial soils in the valleys with small rivers and floodplains are preferred for agriculture. Residents regularly work the homestead plantations whereas the commercial sugarcane plantations only receive seasonal workers between the months of November to April for harvesting. Occasionally the plantations are sprayed with insecticides.

In northeast Brazil, the known sand fly vectors of ATL are Nyssomyia whitmani and Migonemyia migonei [17,18]. In this study site only Ny. whitmani is present representing >97% of the total sand flies captured monthly over two years; seasonal low abundance variably occurs during 3–5 months from May to July/October [19], which roughly corresponds to the coldest (August and July) and the driest (September to December) months.

Between 2011 and 2014, 1,082 ATL cases were reported in Pernambuco including 14 in the Amaraji municipality. Leishmania skin test surveys in the region have shown that the actual number of infections is substantially higher, and that the numbers of reported clinical cases does not reflect the levels of transmission [20]. The apparent reductions in case numbers in recent years is partially explained by the high prevalence of acquired immunity against clinical infection.

Small mammal trapping regime

Between May 2012 and August 2014, 20 rounds of CMR were conducted each of two to 11 days duration (S1 Table). Most of the trapping effort (89.5%) was in Engenho Raiz de Dentro due to seasonally limited access to Engenho Refrigério. 15 homesteads located 200-500m apart, as typical in this region, were selected as central reference points for sampling, chosen due their proximity to peri domestic plantations and Atlantic forest patches in order to facilitate repeat trapping. Eighty numbered Tomahawk traps (45cm x 21cm x 21cm) were set per night, 5 traps in domestic habitats (inside houses and outside within the household boundary- labelled here as peri domestic sites), 55 in adjacent plantations, and 20 in Atlantic Forest remnants. Traps were placed in a single line 2m apart in forest and plantations, and placed in the kitchen, bedroom and/or living room inside houses, and in the yard, terrace and/or near animal shelters outside houses. Traps were placed in the same reference points on each CMR round. Traps were baited with pieces of banana and sugar cane.

Animal processing

Captured animals were identified to species using the nomenclature following previous published references to Akodon arviculoides, Oryzomys eliurus, Oryzomys subflavus and Oxymycterus angularis and herein referred to as Akodon cursor, Oligoryzomys nigripes, Cerradomys subflavus and Oxymycterus dasytrichus, respectively [21]. Animals were sexed and assigned to an age-class at first capture (YY- juvenile; JA- young adult; AA-mature adult) based on their body size by experienced field technicians. Rodents were anesthetized with recommended doses of pre-anaesthetic Xylazine 2%, and aesthetic Ketamine 10%. A blood sample of 0.5 to 1mL was obtained from the retrobulbar venous sinus with a sterile Pasteur pipette. A skin fragment was collected from the right ear of each animal by biopsy punch to provide a 25mg skin sample. Any ectoparasites were removed. Animals were then marked with an individual nano microchip (1.25mm x 7.0mm) (Micro-ID, Burgess Hill, West Sussex, UK) read using a Micro-ID hand scanner (Model MIDO1C). Animals were then placed in micro-isolators in the field station for recovery. The next morning when fully alert the rodents were released at the site of capture. Animals that were subsequently re-captured within the same sample round were identified, recorded and released without further diagnostic sampling. Recaptured individual animals were diagnostically sampled at a median interval of 50 (IQR: 35–91) days. During the study, 32 marsupials were also captured, comprising 21 Didelphis albiventris (1 inside a house; 15 in plantations; 5 in forest), 6 Marmosa spp., (3 in plantations; 3 in forest) and 5 Monodelphis domestica domestica (5 in plantations). As marsupials were not the focus of the study they were not sampled and released immediately.

Xenodiagnosis

Xenodiagnoses using wild-caught Ny. whitmani sand flies.

Xenodiagnosis was performed using wild Ny. whitmani sand flies which were trapped by CDC light traps (set 6pm-6am) in animal shelters, stables and chicken coops in the nearby locality of Engenho Raiz Nova where Ny. whitmani is the dominant sand fly species [22]. After separating the blood-fed and non-blood-fed females, the latter were placed in a Barraud cage into which anaesthetised rodent(s) were introduced and flies allowed to feed for 1 hour in darkness. The rodent(s) submitted for xenodiagnosis on any single day depended on the number of rodents and sand flies captured on the night previous. Post exposure, blood-fed and un-fed female sand flies were removed from the cage and separated individually into a tube containing 70% ethanol for subsequent DNA extraction and qPCR testing. The aim was to expose each captured rodent to a minimum of 20 female Ny. whitmani. In the final analysis, a median 34 (95% C.L. 19–48) blood-fed female flies were tested from each xenodiagnosis experiment.

For analyses of the xenodiagnosis data using wild-caught Ny. whitmani, the observed proportion of sand flies infected with L. (V.) braziliensis was adjusted to account for an expected background level of Leishmania natural infection estimated from an independent sample of 250 Ny. whitmani captured during contemporary trapping rounds in the same area. Of these, 6% (15/250) proved infected, hence the value of the upper 95% confidence limit (C.L.) of the mean infection rate in this control sample was subtracted from the observed fly infection rates by xenodiagnosis. Similarly, the upper 95% C.L. of the L. (V.) braziliensis parasite loads estimated by qPCR in the naturally infected sand flies was subtracted from the values detected in xenopositive Ny. whitmani.

Xenodiagnoses using laboratory reared Lutzomyia longipalpis.

For comparison with a likely secondary sand fly vector species, Lutzomyia longipalpis females were obtained from a laboratory colony held at IAM/FIOCRUZ established from flies captured in Passira municipality, located in the Agreste zone of Pernambuco State, and reared as described [23]. Rodents were anesthetized and placed in cages as above with Lu. longipalpis unfed females for 1 hour. After removal, groups of 25 blood-fed female sand flies were maintained in separate Barraud cages supplied with sugar solution for five to 7 days to allow parasite development. They were then stored individually in 70% ethanol for subsequent DNA extraction and qPCR testing. A mean 81 (95% C.L. 42, 121) Lu. longipalpis females were tested from each xenodiagnosis experiment. As Lu. longipalpis sand flies were laboratory bred, adjustments to infection rates as above were not required.

Molecular tests

DNA was extracted from 200μL blood samples using QIAamp DNA Mini Kit Blood (QIAGEN)® and from 25mg skin biopsy samples using QIAamp DNA Mini Kit (QIAGEN)®. The conventional PCR was performed as described previously [24] using the B1 and B2 kDNA primers (B1: 5′-GGGGTTGGTGTAATATAGTGG-3′ / B2: 5′-CTAATTGTGCACGGGGAGG-3′) as previously described [25]. The qPCR protocol’s final volume was 50 ml containing primers [kDNAf-ATGCCTCTGGGTAGGGGCGTTC/kDNAr1-GGGAGCGCGGCCCACTATATT] (5 pmol/ml) and 2X SYBRGreen Master Mix (Applied Biosystems) [26,27]. For qPCR, a standard curve was included on each plate using L. (V.) braziliensis reference strain LB2903 at seven serial dilutions (1ng, 100pg, 10pg, 1pg, 100fg, 10fg, to 1fg). The test sample OD values were compared to the standard curve and converted into arbitrary units as reported hereafter. Ct values <33 were considered valid based on prior validation of qPCR on each sample type. Test samples were run in duplicate using standard conditions (40 cycles of denaturation at 95°C for 15s and annealing at 60°C for 1min), in ABI Prism 7500 Sequence Detection System (Applied Biosystems) with 7500 Software (version 2.0.5) [27].

Statistical analyses

Estimates of infection and infectious status were computed by fitting (0/1) data to mixed-effects binomial complimentary log–log models including rodent ID and time (days) from first capture as random effects, and reported as risk ratios (RR). In cases where models failed to converge, a cluster term (animal ID) was included to control for likely autocorrelation. Models included variables characterising time (days) since first capture, month, and year of trapping (3 levels), age-class (3 levels), sex (2 levels), ecotype (3 levels) and rodent species (8 levels), as appropriate. Similarly, the proportion of sand flies infected by xenodiagnosis was modelled as binomial data using Generalise Linear Models.

Complimentary log–log model fits were achieved by Gauss–Hermite numerical adaptive quadrature of the random-effects estimators (quadchk routine in STATA), validated using 16 integration points in model run comparisons to confirm quadrature fitting accuracy; model runs showed 0.01% variation in resulting estimates, and were thus considered to be reliable. Potential explanatory variables were evaluated in a preliminary step by log–likelihood ratio test (LRT) of nested models.

Data were analysed in STATA v.15 (StataCorp LP, College Station, TX).

Infection risk

To standardise the infection risk relative to trapping effort and trapping success between ecotypes where captured, crude capture rates were adjusted by: (1) where given

a = ∑ number of successful capture nights of rodent species i within ecotype j

b = ∑ number of attempted trap nights within ecotype j

i = the rodent species

j = the ecotype

p = the number of rodent capture events that showed evidence of infection

n = ∑ number of rodent capture events where clinical samples were diagnostically tested

Results

Sampling

Between May 2012 and August 2014, a total 9,920 single trap nights were accomplished in 20 trapping rounds over 27 months. Of these 620 trap nights were in the peri domestic environment, 6,820 in plantations and 2,480 in forest patches, a ratio of 1: 11: 4 respectively (Tables 1 and S1).

Trapping success

A total 603 rodent individuals of eight species were (re)captured on 1,051 occasions, of which the most abundant species, indicated by the capture records, were Nectomys squamipes (n = 245 individuals representing 41% of captured individuals), followed by Rattus rattus (n = 148, 25%), and Necromys lasiurus (n = 83, 14%) (Table 1).

New individuals were captured and recaptured over much of the trapping period, with a noticeable fall-off in trapping success in the latter four trap rounds in part associated with the reduced trapping effort (Fig 1 and S1 Table). With exception of Ol. nigripes and C. subflavus for which only 2–4 individuals were captured (none recaptured), the median number of recapture events per individual rodent was generally similar across species, as was the median recapture interval (45 to 69 days) (S2 Table). N. squamipes was the exception and most common species recaptured up to 11 times each in plantations and forest habitats (S4A and S4C Fig), which was more frequent than the other species, each being recaptured up to 4 times (S4A and S4B Fig). R. rattus was recaptured on relatively fewer occasions and showed a greater affinity to peri domestic and plantation environments (S4D Fig), and was the only species captured in all three sampled habitat types (Table 1). The abundance and affinity with a particular habitat varied between rodent species, however the proportion of total (re)capture events was generally highest in plantations when adjusted for trapping effort (Table 1 and S4C and S4D Fig). Standardised for ecotype and rodent species, trapping success was 33.7%, 47.6% and 18.7% in the peri domestic environment, plantations and forest habitats, respectively (Table 1). Trapping success was significantly higher in plantations than in the peri domestic environment (χ21 = 7.28, P = 0.007), and higher in peri domestic than forest habitat (χ21 = 14.75, P = 0.0001).

thumbnail
Fig 1. Number of rodents captured in each trapping round (A) individuals at first capture; (B) all individuals recaptured.

https://doi.org/10.1371/journal.pntd.0010996.g001

thumbnail
Table 1. Rodent trapping success: number and proportion of total (re)capture events across sampled ecotypes.

https://doi.org/10.1371/journal.pntd.0010996.t001

Infection prevalence

Rodents were defined as infected at each (re)capture if positive by one or more diagnostic test including conventional PCR (blood samples), qPCR (blood or ear skin tissue samples), and xenodiagnosis. Of the 602 animals and 999 sample test results, 42.5% and 42.1% respectively were considered infected (Table 2). Omitting Ol. nigripes and C. subflavus (2 and 4 diagnostic samples only respectively), all species showed a lower period prevalence of infection compared to N. squamipes (z>-3.26, p<0.001), and Holochillus sciureus (z>-2.01, p<0.045), whereas those of N. squamipes and H. sciureus were not dissimilar (z = -0.77, p = 0.442) (Table 2).

The proportion of blood samples that tested positive by qPCR (38.9% [382/983]) was significantly greater than positive by conventional PCR (5.8% [58/999]) (Tables 2 and S3). QPCR sensitivity also differed between clinical samples: only 1.8% of skin tissue samples tested positive compared to 38.9% of blood samples (Table 2). Of the 14 animals with positive skin samples, 11 were negative by blood sample. Based on conventional PCR alone, the infection prevalence of H. sciureus was greater than that of all other species including N. squamipes (z>2.35, p<0.019) (S3 Table).

thumbnail
Table 2. Number and proportion of rodents and clinical samples classified as infected1 and positive by qPCR for L. (V.) braziliensis.

https://doi.org/10.1371/journal.pntd.0010996.t002

Infection risk through time

L. (V.) braziliensis infection risk in N. squamipes significantly increased from time of first capture (z = 2.47, p = 0.014), illustrated by the rise in infection prevalence amongst captured animals at sequential trapping rounds (Fig 2A). Similarly, initial increases were observed in R. rattus (Fig 2A) and in H. sciureus, Ne. lasiurus, and O. dasytrichus (Fig 2B) (z>3.56, p<0.0001). In the case of R. rattus, H. sciureus and Ne. lasiurus this rise was followed by a significant decrease in infection risk by the end of the study (test of days2 quadratic term: z>-2.62, p<0.009); no such decline was observed for N. squamipes or O. dasytrichus. The risk of rodent infection with time from first capture did not significantly vary between ecotypes for any of the rodent species with adequate data (test of ecotype × days interaction term: z<1.36, p>0.175 in each case).

thumbnail
Fig 2. The proportion of individuals infected from time of first capture.

(A) N. squamipes and R. rattus; (B) H. sciureus, Ne. lasiurus, and O. dasytrichus. Data were aggregated over the time frame to achieve adequate numbers.

https://doi.org/10.1371/journal.pntd.0010996.g002

Seasonal variation in infection prevalence by calendar month was indicated for N. squamipes and for the other species collectively, being highest in April to July, and again in December-January (Fig 3). This pattern did not correspond to the previously reported seasonal variation in Ny. whitmani abundance [28], as also shown in Fig 3.

thumbnail
Fig 3. Proportion of rodent populations infection by calendar month.

Prevalence calculated for N. squamipes, and other rodent species by aggregating data. Superimposed is the mean number of Ny. whitmani sand flies per month collected by CDC light trap in this and nearby foci (published data from Brandão-Filho et al 1994 [19]).

https://doi.org/10.1371/journal.pntd.0010996.g003

L. (V.) braziliensis parasite loads in rodent clinical samples

L. (V.) braziliensis parasite loads in infected clinical samples were highly variable and the frequency distributions over-dispersed (Table 3 and S5A and S5B Fig). There was poor correlation between log10 transformed parasite loads in skin and blood samples from the same animal (rs2 = 0.087, p = 0.575). Median parasite loads in blood samples tended to be lower in H. sciureus, N. squamipes, O. dasytrichus and R. rattus compared to in Ne. lasiurus (S6 Fig), though the differences failed to reach statistical significance (z<1.95, p>0.051). There were no differences in blood parasite loads between any of the other species (z<-1.12, p>0.261). Considering the rodent populations as a whole, there was a rise in mean log10 transformed parasite loads in blood samples from time of 1st observed infection (z = 5.35, p<0.000), followed by a significant decline (z = -4.52, p<0.001).

thumbnail
Table 3. Geometric number of L. (V). braziliensis parasites in rodent clinical samples estimated by qPCR.

https://doi.org/10.1371/journal.pntd.0010996.t003

Infection and demography

Controlling for study design variables and rodent species and individual IDs, the generalised probability of infection across rodent species was not dissimilar between sexes (males: 101/327 vs females: 75/275; z = 0.62, p = 0.533), nor between age-classes categorised at first sample: juvenile (JJ) 47.1% (72/153); young adult (JA) 48.8% (137/281); mature adult (AA) 37.5% (212/565) (z<1.74, p>0.082). For individual species with adequate samples, there were no differences between age-classes or sex in N. squamipes, R. rattus (z<1.87, p>0.06), or H. sciureus (Pearson’s χ21<0.80, p>0.36), whereas the risk of infection in Ne. lasiurus and O. dasytrichus at first sample were similarly greater in both young adult and mature adult age-classes compared to in juveniles (z>10.3, p<0.0001).

Infection by ecotope

The three ecotopes appeared to support similar levels of infection (Table 4). However, the period prevalence varied between rodent species and between ecotopes: whereas R. rattus showed a similar infection prevalence across the three ecotypes, potentially higher prevalences were seen in A. cursor in forest compared to in plantations, and the opposite for O. dasytrichus and Ne. lasiurus, though no statistical differences were detected (Table 4).

thumbnail
Table 4. Number (n) and proportion of (re)capture events classified as infected1 by rodent species and ecotype.

https://doi.org/10.1371/journal.pntd.0010996.t004

Applying Eq 1 to account for capture success and trapping effort, a standardised proportional infection risk was calculated for each rodent species by ecotype (Table 5). Plantations represented the highest percentage (70.3%) of the total estimated infection risk compared to 11.3% in the peri domestic environment and 18.4% in forests. For individual species, plantations represented a higher risk of infection than in forests; R. rattus was the only species that differed from this pattern (Table 5).

thumbnail
Table 5. Proportional risk1 of infection in rodent species in sampled habitats.

https://doi.org/10.1371/journal.pntd.0010996.t005

Ecotype loyalty

Of the 186 individuals captured on at least two independent trap rounds, 42 (22.6%) were captured in more than one ecotype, including 21/106 (19.8%) of N. squamipes and 19/30 (63.3%) of R. rattus (Table 6). Of these, 33/42 (78.6%) were classified as infected, and at least 2 N. squamipes and 2 R. rattus tested were xenopositive.

thumbnail
Table 6. Number (proportion) of individual rodents captured in more than one ecotype.

https://doi.org/10.1371/journal.pntd.0010996.t006

Infectiousness to sand flies

Xenodiagnosis was performed 44 times by exposing 5 rodent species to wild-caught Ny. whitmani sand flies, targeting N. squamipes for repeat sampling (n = 39) over an approximate 13-month period (Table 7). Thirty-nine of the 44 xenodiagnosis trials were conducted on plantation-captured animals which excluded the possibility to test for differences in infectiousness between habitat types.

thumbnail
Table 7. Infectiousness of free-ranging wild rodents to wild caught Ny. whitmani sand flies, or to laboratory-bred Lu. longipalpis, measured by xenodiagnosis.

https://doi.org/10.1371/journal.pntd.0010996.t007

Of the 41 xenopositive trials conducted using Ny. whitmani, 34 trials were xenopositive at the time when the animals tested positive by qPCR/PCR, and 7 trials xenopositive (all on N. squamipes) when the animal tested qPCR/PCR negative. Two animals exposed when qPCR/PCR positive were xenonegative. Details for individual rodents, qPCR/PCR results and the number of sand flies tested, are available in supplementary material (S4 Table).

Five xenodiagnostic trials were conducted using colony reared Lu. longipalpis exposed to three N. squamipes and two R. rattus individuals, all of which were xenopositive (Tables 7 and S5). All these animals except one N. squamipes were PCR/qPCR positive at the time of xenodiagnosis (Tables 7 and S5).

N. squamipes infectiousness to Ny. whitmani

Over the 13-month xenodiagnosis study period, members of the N. squamipes population were infectious (Fig 3). The proportion of sand flies that individuals infected varied from 0.0 to 0.657. Aggregating individual data from the time of 1st capture (n = 4–12 xenodiagnosed individuals in each of 5 periods over 388 days), neither the population mean probability of being infectious, nor the proportion of total sand flies that they infected, significantly vary with increasing time (z<-0.970, p>0.146). The median proportion infected was 0.48 (IQR: 0.301–0.642); range: 0.00–0.657, N = 39), illustrated in Fig 3). Aggregating the xenodiagnosis data according to the periods of vector seasonal abundance (low season: variably May-July/October; high season: August/November-May), the proportion of exposed flies infected was higher in the low compared to in the high season (z>-9.14, p<0.001, varying the delimiting months defining the season). The respective median proportions infected were 0.64 of 716 flies (IQR: 0.600–0.657) in 23 trials, and 0.08 of 645 flies (IQR: 0.024–0.307) in 16 trials; the respective proportions were similar for the alternative definition of season.

thumbnail
Fig 4. The association between the log10 transformed L. (V.) braziliensis parasite loads estimated by qPCR in rodent blood samples and in Ny. whitmani sand flies infected during xenodiagnosis.

https://doi.org/10.1371/journal.pntd.0010996.g004

Male N. squamipes appeared to be more infectious than females, infecting a median 65% (IQR: 62.4, 65.7, n = 12) of blood-fed sand flies compared to 43% (IQR: 5.1, 59.9, n = 13) infected by females (z = -4.16, p<0.001). No differences were detected between age-classes (z<-0.74, p>0.456). Aggregated data for the rodent species other than N. squamipes indicated similar trends: a non-significant difference between age-classes (z<-1.02, p>0.308), and a tendency for males to be more infectious to sand flies than females: 52% (IQR: 44.7, 65.7, n = 3) versus 55% (IQR: 24.8, 63.4, n = 4) (z = -2.16, p = 0.031).

Clinical L. (V.) braziliensis loads in sand flies and infectious rodents

The frequency distributions of the log10 transformed parasite loads in sand flies post xenodiagnosis were over-dispersed (S5C Fig). There was no correlation between the parasites loads in rodent blood and in the xenoinfected sand flies (Fig 4), nor in the contemporary paired samples for individuals (rs2 = 0.076, p = 0.602). Likewise, rodent blood parasite loads were not correlated with the proportion of exposed sand flies infected at point xenodiagnosis (rs2 = 0.234, p = 0.106). Few skin biopsy samples proved positive by qPCR for comparative analyses.

Discussion

Previous field studies in the Zona da Mata in Pernambuco, Brazil, detected Leishmania infection in a mosaic of wild and synanthropic rodents and domesticated animals, indicating that they may play a role in zoonotic transmission of L. (V.) braziliensis [7]. The present CMR study in this region focussed on rodents as they represented the largest biomass [7,24]. L. (V.) braziliensis infections were confirmed in 603 individuals of six of the eight rodent species captured. Period prevalence ranged from 12.5% to 61.5% being statistically higher in N. squamipes and H. sciureus populations (Table 2). No infections were detected in the two Ol. nigripes nor four C. subflavus individuals captured.

The epidemiological significance of rodent infections was further investigated by xenodiagnosis allowing the known vector Ny. whitmani to freely blood-feed on N. squamipes, Ne. lasiurus, R. rattus, A. cursor and O. dasytrichus. All individuals of the five species examined proved infectious except two individual N. squamipes in 3 xenodiagnosis experiments (Table 7). The median proportion of blood-fed flies infected by each species was relatively high (>45%) with exception of O. dasytrichus (1%), suggesting that at least the more abundant rodent species (N. squamipes 41% of total individuals captured, R. rattus 25%, and Ne. lasiurus 14%), significantly contribute to transmission. As the blood-fed Ny. whitmani from these experiments were killed and preserved for molecular testing immediately after blood-feeding, the results reflect potential infectiousness rather than transmissibility per se as the latter would require detection of metacyclic promastigotes post extrinsic incubation. Xenodiagnosis using the known permissive species, Lu. longipalpis [29], a possible secondary vector of L. (V.) braziliensis [10], demonstrated the infectiousness of all 3 N. squamipes and 2 R. rattus individuals upon dissection of flies 5–7 days after feeding. In this case, the results more likely reflect transmissibility.

Xenodiagnosis data for naturally infected wildlife hosts of Leishmania are few, and for L. (V.) braziliensis are limited to a single experimental infection study of N. squamipes, Ne. lasiurus and R. rattus using rodent colonies established with wild-caught animals from the current study region [12]. In that study, 10/18, 18/18 and 6/18 individuals respectively were infectious to Lu. longipalpis, results which help support the findings of the current study. The proportions of sand flies infected by these naturally and experimentally infected rodents were substantially higher than reported from xenodiagnosis experiments on regionally suspected rodent species experimentally infected with L. major strains, L. donovani, L. tropica or L.(V.) panamensis [3032], or lagomorphs and wild canids naturally infected with L. infantum [8,28,33].Some of the variation is attributed to the different detection techniques and timing; molecular methods are generally more sensitive than direct microscopic examination [7]. Notwithstanding, the proportion of flies infected here approach those reported for European and Latin American sand fly vectors infected with L. infantum by domestic dogs, the proven reservoir of L. infantum [8, 9].

Additional criteria to incriminate a reservoir host is that it can maintain infection through the non-transmission period i.e. during periods of low or zero vector abundance [9]. The recaptured rodents in this study showed chronic parasitic infections in the absence of any apparent clinical signs. The average duration of detected infection in individuals and in the whole rodent population (Fig 3), exceeded the period of seasonally low abundance of Ny. whitmani, which is variably 3–5 months from May to July/October in this region [19]. Targeting the N. squamipes population for repeat xenodiagnosis, we demonstrate persistence transmission potential over the 13-month observation period, infecting a median 48% (IQR: 30.1%-64.2%) of exposed Ny. whitmani, and observed significantly greater levels of infectiousness during the reported periods of low vector abundance compared to during the periods of high abundance (Fig 3). Since the rodent species occurred in sympatry, we were unable to demonstrate that L. (V.) braziliensis can be maintained by any single species, but comparative longitudinal xenodiagnoses studies would be useful to partition the basic case reproduction number R0 to determine their relative epidemiological contribution to transmission [8].

In a related study and complimentary analyses of the current data, we demonstrate a high force of infection λ = 0.22/month (SD: 0.183) in the rodent population as a whole, with 50% of the population infected as juveniles (2–3 months old), and 80% infected by 8 months old, associated with a short average prepatent period (time to detectable infection) of 1–3 months [34]. The longitudinal infection incidence appeared highest in Ne. lasiurus and lowest in R. rattus but the variation was broad excluding detection of statistical differences between species [34].

Infected rodents were (re)captured in all the major habitat types under study and where humans frequent, suggesting that humans are also at risk of infection with L. (V.) braziliensis. The household plantations, typically comprising bananas, pineapples and vegetable cash crops, proved to be the habitat of greatest rodent abundance, infection prevalences, and relative infection risk (Tables 4 and 5). Vector presence in each of these ecological habitats is well documented [19]. We did not sample the commercial monoculture sugarcane plantations which are extensive in the region, as human presence is limited to seasonal planting and harvest periods by a seasonal labour force. Sand flies are absent in this habitat but on rare occasions rodents have been captured.

Crude infection prevalences were similar or possibly higher in Atlantic forest patches for some rodent species, notably A. cursor, N. squamipes and R. rattus. The two habitats yielded not dissimilar rodent species compositions, but importantly we show that 22.6% of individuals captured at least twice were captured in more than one habitat type including 19.8% of N. squamipes and 66.3% of R. rattus, indicative of the likely connectivity between these habitats (Table 6). Of these animals, 78.6% were infected at some point during the study. R. rattus linked the peri domestic environment to both plantation and forest. Thus, infected rodent dispersions between habitats occurs and we speculate that individual animal home ranges within the different habitats overlap that of sympatric species, thereby providing opportunities for inter-species and inter-habitat transmission. Many of the species in this study have been recorded in a diversity of habitats, with or without L. (V.) braziliensis infections, indicating their ecological adaptability [24,3542].

It is also conceivable that there are concurrent domestic and peri domestic/sylvatic transmission cycles in this region each involving one (or more) primary host species. In this case, the cycles would be linked by the vector Ny. whitmani which is regularly trapped in all three habitats [19]. The evidence to date suggests that N. squamipes and R. rattus may represent the primary hosts in these putative cycles respectively but this cannot be confirmed by the current data. Only R. rattus was captured in and around houses, and in all three sampled ecotypes under study; it is a synanthropic species frequently observed in the walls and ceilings of houses and has consistently shown high Leishmania infection rates in this focus [19] and elsewhere [24,43]. Additional domestic hosts may include dogs of which 20% were previously shown PCR positive [7] with higher prevalences reported elsewhere [43], though xenodiagnosis experiments have yet to be conducted. N. squamipes is physiologically adapted to a semi-aquatic life, though recognised as a generalist using aquatic and terrestrial environments alike [44]. It was the largest rodent species in this study (body weight range: >170–400 gm) and the most abundant species, which could be important for disproportionately attracting blood-seeking sand flies. Five species of N. squamipes occur over an extensive area of the Americas [45,46] suggesting that they may play a primary role in L. (V.) braziliensis transmission beyond northeast Brazil. Seasonally stable populations dominate the small mammal assemblage in ecologically diverse habitats elsewhere [44,47,48], though not exclusively so; species compositions clearly vary between ecological foci [24,35,37,43] and possibly through time. For example, in a different ATL focus in Pernambuco, 52% of 46 small mammals captured in 2003 were marsupials (33.3% infected) compared to 0% of 118 animals captured in 2006 using the same trapping methods [24].

Xenodiagnosis is the gold standard to measure host infectiousness, but it is logistically difficult to perform at scale. Quantification of L. (V.) braziliensis parasite loads in rodent blood by qPCR was found to be imprecise as a potential surrogate marker of host infectiousness. Indeed, there was generally poor correlation between corresponding parasitaemia levels in rodent blood and skin samples, and between skin samples and xenopositive sand flies. On the other hand, 83% of rodents with blood positive for L (V.) braziliensis were xenopositive suggesting that qPCR positivity may be a useful qualitative indicator of infectiousness.

In contrast to the surprisingly few ear tissue samples being qPCR positive, the L. (V.) braziliensis experimental study identified parasite loads in rodent ear tissue to be the best indicator of infectiousness to Lu. longipalpis [12]. Such an association has been shown also for cohorts of dogs naturally infected with L. (L.) infantum [49]. Despite reports of greater PCR sensitivity testing rodent blood compared to rodent skin [41,50], the current results may be due to unidentified technical reasons. We infer this from the success of molecular detection of L. (V.) braziliensis in ear skin similar to in spleen and liver samples from N. squamipes, Ne. lasiurus and R. rattus experimental animals [12], in naturally infected rodent populations [43], and in longitudinal follow-up of rodents experimentally infected with other Leishmania species, where the ear pinnae inoculation site remained a significant source of parasites [12,51,30]. Furthermore, finding significant L. (V.) braziliensis parasite loads in both skin and internal organs of experimental [12] and naturally infected [37,43] wild rodents in Brazil, and the high parasite loads in the blood of the current study animals, suggests that metastasis may be an intrinsic character of L. (V.) braziliensis in rodents, as is the case in humans and other animals [5254]. Parasite tropism and their non-continuous distributions across the skin [55] will clearly influence the accessibility to blood-feeding sand flies, and thus the hosts’ infectious potential.

The chronic and asymptomatic nature of Leishmania infection in wild rodents observed in this population is not atypical [7,43,56,57]. Even following experimental inoculation of N. squamipes, Ne. lasiurus and R. rattus with high dose (5.5 x 106/ml) L. (V.) braziliensis, the skin lesions that developed in 13/26 animals self-resolved within one month of onset [12]. In that study, all animals were infectious to Lu. longipalpis following spontaneous clinical recovery. Similarly, suspect rodent reservoir species experimentally infected with L. major or L. donovani presented largely asymptomatic infections for 20–30 weeks pre-necropsy; all three L. major infected rodent species were infectious when asymptomatic [30,58].

Conclusions

The collective data on infection prevalence, infectiousness and rodent spatial abundance support the view that the rodent assemblage in this focus represents an important wildlife source of L. (V.) braziliensis for onward transmission. N. squamipes at least plays a significant role, and probably also R. rattus, Ne. lasiurus, and A. cursor. Peri domestic plantations appear to be locations of highest infection risk for rodents, and by extrapolation, humans, forming one part of a wider eco-epidemiological network linking domestic and forest habitats. L. (V.) braziliensis infections are documented in a multitude of other wild animal species in this region and elsewhere (reviewed [59]) but infectiousness studies remain to be conducted. A complete understanding of the geographical complexity of enzootic transmission needs to evaluate the relevance of reported mixed Leishmania spp. host infections [24, 43] and regional genetic diversity of L. (V.) braziliensis [60]. In Pernambuco three genetic groups have so far been identified [61], one and two, that have high heterozygosity, occur in areas where Ny. whitmani is the vector, such as Amaraji. The third variant is associated with forests where other sand fly species are suspected vectors. Transmission pathways ultimately depend on the prevailing conditions of the ecosystem and climatic drivers, undoubtedly fuelled by environmental modifications by humans.

Supporting information

S1 Fig. Satellite image showing the location of the capture areas Engenho Raiz de Dentro (8° 42.7944’S, 35° 44.7078’W) and Engenho Refrigério (8° 42.7963’S, 35° 49.4742’W) in relation to Amaraji town, Amaraji municipality, Pernambuco State, Brazil.

https://doi.org/10.1371/journal.pntd.0010996.s001

(DOCX)

S2 Fig. A homestead in the Engenho Raiz de Dentro area.

The open garage on the left side of the house was used as a field laboratory.

https://doi.org/10.1371/journal.pntd.0010996.s002

(DOCX)

S3 Fig. A homestead in the Engenho Refrigério area.

Note the proximity of the banana plantation to the house and on the other side of the lane a fenced field for livestock.

https://doi.org/10.1371/journal.pntd.0010996.s003

(DOCX)

S4 Fig. Frequencies of recaptures of rodent species over the 20 trapping sessions.

(A) Ne. lasiurus, N. squamipes and R. rattus; (B) A. cursor, H. sciureus and O. dasytrichus; (C) N. squamipes in forest and plantations; (D) R. rattus in plantations, forest, and the domestic environment.

https://doi.org/10.1371/journal.pntd.0010996.s004

(DOCX)

S5 Fig.

Frequency distributions of L. (V.) braziliensis parasite loads detected by qPCR in (A) rodent blood (per 200μL); (B) ear skin biopsies (per 25mg); (C) Ny. whitmani sand flies (per individual fly). In (A) N. squamipes alone (right), all other rodent species (left). Note the differences in the Y-axis scales.

https://doi.org/10.1371/journal.pntd.0010996.s005

(DOCX)

S6 Fig. Boxplot of the median (IQR) L. (V.) braziliensis parasite loads detected by qPCR in rodent blood samples (dark bars) and Ny. whitmani sand flies post exposure to rodents during xenodiagnosis (light bars).

https://doi.org/10.1371/journal.pntd.0010996.s006

(DOCX)

S1 Table. Trapping schedule showing the number of trap-nights (total single traps set over trapping period).

https://doi.org/10.1371/journal.pntd.0010996.s007

(DOCX)

S2 Table. Median recapture interval (days) and number of times individual rodents were recaptured.

https://doi.org/10.1371/journal.pntd.0010996.s008

(DOCX)

S3 Table. Infection prevalences based on conventional PCR compared to quantitative PCR (qPCR) of rodent blood samples using PCR target primers specific to parasites of the genus Leishmania (Viannia) braziliensis.

https://doi.org/10.1371/journal.pntd.0010996.s009

(DOCX)

S4 Table. Infection and infectiousness of individual rodents at the time of xenodiagnosis exposing Ny. whitmani sand flies, and the associated L. (V.) braziliensis parasite loads in the rodent blood and xenopositive sand flies.

https://doi.org/10.1371/journal.pntd.0010996.s010

(DOCX)

S5 Table. Infection and infectiousness of individual rodents at the time of xenodiagnosis exposing Lu. longipalpis sand flies, and the associated L. (V.) braziliensis parasite loads in the rodent blood and xenopositive sand flies.

https://doi.org/10.1371/journal.pntd.0010996.s011

(DOCX)

Acknowledgments

We are grateful to biology students Mariana Verissimo, Priscila Carla da Silva, Luana Patrícia Brito for their laboratory support. Hélio Valença, Amilton Silva, (technicians of FUNASA/MS—National Health Foundation/Ministry of Health) and Fernando Silva (technician of IAM/FIOCRUZ) for their support in the field work. Also to Maria Edileuza Felinto Brito, Éricka Almeida and Andréa Sales (IAM/FIOCRUZ) for their support in the laboratory, and Romário Melo and José Egberto Silva for supporting work with rodents and sandflies. We would also like to thank Professor James L. Patton (Museum of Vertebrate Zoology, University of California, Berkeley, USA) for his help with the rodent taxonomy. We are grateful also for institutional and logistical support from IAM/Fiocruz to vehicle availability for field work.

References

  1. 1. Espinosa OA, Serrano MG, Camargo EP, Teixeira MMG, Shaw JJ. An appraisal of the taxonomy and nomenclature of trypanosomatids presently classified as Leishmania and Endotrypanum. Parasitology. 2018;145(4):430–42 (Epub 2016 Dec 15). Epub 2016/12/16. pmid:27976601.
  2. 2. Kawa H, Sabroza PC. Espacialização da leishmaniose tegumentar na cidade do Rio de Janeiro. Cadernos de Saúde Pública. 2002;18:853–65.
  3. 3. Monteiro WM, Neitzke HC, Silveira TGV, Lonardoni MVC, Teodoro U, Ferreira MEMC. Pólos de produção de leishmaniose tegumentar americana no norte do Estado do Paraná, Brasil. Cadernos de Saúde Pública. 2009;25:1083–92.
  4. 4. Shaw J. The importance of understanding enzootic cycles in the epidemiology of zoonotic diseases with special reference to the American leishmaniases. Trans R Soc Trop Med Hyg. 2019;113(3):108–9. pmid:30358870
  5. 5. Lainson R, Shaw JJ. 17 New World Leishmaniasis. In: Cox FEG, Kreier JP, Wakelin D, editors. Topley & Wilson’s Microbiology and Microbial Infections. 10th ed. London, Sydney, Auckland: Hodder Arnold; 2005. p. 313–49.
  6. 6. de Brito ME, Andrade MS, Dantas-Torres F, Rodrigues EH, Cavalcanti Mde P, de Almeida AM, et al. Cutaneous leishmaniasis in northeastern Brazil: a critical appraisal of studies conducted in State of Pernambuco. Rev Soc Bras Med Trop. 2012;45(4):425–9. pmid:22836662.
  7. 7. Brandão-Filho SP, Brito ME, Carvalho FG, Ishikawa EA, Cupolillo E, Floeter-Winter L, et al. Wild and synanthropic hosts of Leishmania (Viannia) braziliensis in the endemic cutaneous leishmaniasis locality of Amaraji, Pernambuco State, Brazil. Trans R Soc Trop Med Hyg. 2003;97(3):291–6. pmid:15228244.
  8. 8. Courtenay O, Quinnell RJ, Garcez LM, Shaw JJ, Dye C. Infectiousness in a cohort of brazilian dogs: why culling fails to control visceral leishmaniasis in areas of high transmission. J Infect Dis. 2002;186(9):1314–20. pmid:12402201.
  9. 9. Quinnell RJ, Courtenay O. Transmission, reservoir hosts and control of zoonotic visceral leishmaniasis. Parasitology. 2009;136(14):1915–34. pmid:19835643.
  10. 10. Alexandre J, Sadlova J, Lestinova T, Vojtkova B, Jancarova M, Podesvova L, et al. Experimental infections and co-infections with Leishmania braziliensis and Leishmania infantum in two sand fly species, Lutzomyia migonei and Lutzomyia longipalpis. Sci Rep. 2020;10(1):3566. Epub 2020/02/29. pmid:32108151; PubMed Central PMCID: PMC7046706.
  11. 11. Walton BC, Shaw JJ, Lainson R. Observations on the in vitro cultivation of Leishmania braziliensis. J Parasitol. 1977;63(6):1118–9.
  12. 12. Andrade MS, Courtenay O, Brito ME, Carvalho FG, Carvalho AW, Soares F, et al. Infectiousness of Sylvatic and Synanthropic Small Rodents Implicates a Multi-host Reservoir of Leishmania (Viannia) braziliensis. PLoS Negl Trop Dis. 2015;9(10):e0004137. pmid:26448187; PubMed Central PMCID: PMC4598029.
  13. 13. Queiroz A, Sousa R, Heine C, Cardoso M, Guimaraes LH, Machado PR, et al. Association between an emerging disseminated form of leishmaniasis and Leishmania (Viannia) braziliensis strain polymorphisms. J Clin Microbiol. 2012;50(12):4028–34. pmid:23035200.
  14. 14. Rego FD, da Rocha Lima A, Pereira AAS, Quaresma PF, Pascoal-Xavier MA, Shaw JJ, et al. Genetic variant strains of Leishmania (Viannia) braziliensis exhibit distinct biological behaviors. Parasitol Res. 2018;117(10):3157–68. pmid:30022292.
  15. 15. CityPopulation. Amaraji Municipal Disrict, Amaraji (Pernambuco). https://wwwcitypopulationde/en/brazil/pernambuco/amaraji/260090605__amaraji/. 2021.
  16. 16. IBGE. Amaraji População. https://wwwibgegovbr/cidades-e-estados/pe/amarajihtml. 2021.
  17. 17. de Carvalho MR, Lima BS, Marinho-Junior JF, da Silva FJ, Valenca HF, Almeida Fde A, et al. Phlebotomine sandfly species from an American visceral leishmaniasis area in the Northern Rainforest region of Pernambuco State, Brazil. Cad Saude Publica. 2007;23(5):1227–32. pmid:17486244.
  18. 18. Rangel EF, Lainson R, Carvalho BM, Costa SM, Shaw JJ. Sand Fly Vectors of American Cutaneous Leishmaniasis in Brazil. In: Rangel EF, Shaw JJ, editors. Brazilian Sand Flies. Basel: Springer International Publishing; 2018. p. 341–80.
  19. 19. Brandao-Filho SP, de Carvalho FG, de Brito ME, Almeida Fde A, Nascimento LA. American cutaneous leishmaniasis in Pernambuco, Brazil: eco-epidemiological aspects in ’Zona da Mata’ region. Mem Inst Oswaldo Cruz. 1994;89(3):445–9. Epub 1994/07/01. pmid:7476230.
  20. 20. Brandão-Filho SP, Campbell-Lendrum D, BME F., Shaw JJ, Davies CR. Epidemiological Survey based on Montenegro skin test confirm increasing burden of cutaneous leishmaniasis in Pernambuco, Brazil. mio. 1999;94 (Suppl. II):82.
  21. 21. Patton JL, Pardiñas UFJ, D’Elía G. Mammals of South America, Volume 2 Rodents. Chicago: University of Chcago Press; 2015 March 2015. 1384 p.
  22. 22. Brandão-Filho SP, Donalisio MR, da Silva FJ, Valença HF, Costa PL, Shaw JJ, et al. Spatial and temporal patterns of occurrence of Lutzomyia sand fly species in an endemic area for cutaneous leishmaniasis in the Atlantic Forest region of northeast Brazil. J Vect Ecol. 2011;36(s1):S71–S6. pmid:21366783
  23. 23. Volf P, Volfova V. Establishment and maintenance of sand fly colonies. J Vect Ecol. 2011;36 Suppl 1:S1–9. Epub 2011/03/05. pmid:21366760.
  24. 24. Lima BS, Dantas-Torres F, de Carvalho MR, Marinho-Junior JF, de Almeida EL, Brito ME, et al. Small mammals as hosts of Leishmania spp. in a highly endemic area for zoonotic leishmaniasis in North-Eastern Brazil. Trans R Soc Trop Med Hyg. 2013;107(9):592–7. Epub 2013/07/23. pmid:23868744.
  25. 25. de Bruijn MH, Barker DC. Diagnosis of New World leishmaniasis: specific detection of species of the Leishmania braziliensis complex by amplification of kinetoplast DNA. Acta Trop. 1992;52(1):45–58.
  26. 26. de Paiva Cavalcanti M, Dantas-Torres F, de Albuquerque SCG, Silva de Morais RC, de Brito ME, Otranto D, et al. Quantitative real time PCR assays for the detection of Leishmania (Viannia) braziliensis in animals and humans. Mol Cell Probes. 2013;27(3–4):122–8. Epub 2013/02/14. pmid:23402826.
  27. 27. Biosystems A. Real time PCR.Absolute quantification getting started guide. Atlanta. 2005.
  28. 28. Michalsky EM, Rocha MF, da Rocha Lima AC, Franca-Silva JC, Pires MQ, Oliveira FS, et al. Infectivity of seropositive dogs, showing different clinical forms of leishmaniasis, to Lutzomyia longipalpis phlebotomine sand flies. Vet Parasitol. 2007;147(1–2):67–76. Epub 2007/04/24. pmid:17449184.
  29. 29. Volf P, Myskova J. Sand flies and Leishmania: specific versus permissive vectors. Trends Parasitol. 2007;23(3):91–2. pmid:17207663.
  30. 30. Sadlova J, Vojtkova B, Hrncirova K, Lestinova T, Spitzova T, Becvar T, et al. Host competence of African rodents Arvicanthis neumanni, A. niloticus and Mastomys natalensis for Leishmania major. Int J Parasitol Parasites Wildl. 2019;8:118–26. Epub 2019/02/12. pmid:30740304; PubMed Central PMCID: PMC6356118.
  31. 31. Svobodova M, Votypka J, Nicolas L, Volf P. Leishmania tropica in the black rat (Rattus rattus): persistence and transmission from asymptomatic host to sand fly vector Phlebotomus sergenti. Microbes Infect. 2003;5(5):361–4. Epub 2003/05/10. pmid:12737990.
  32. 32. Travi BL, Ferro C, Cadena H, Montoya-Lerma J, Adler GH. Canine visceral leishmaniasis: dog infectivity to sand flies from non-endemic areas. Res Vet Sci. 2002;72(1):83–6. pmid:12002643.
  33. 33. Molina R, Jiménez MI, Cruz I, Iriso A, Martín-Martín I, Sevillano O, et al. The hare (Lepus granatensis) as potential sylvatic reservoir of Leishmania infantum in Spain. Veterinary Parasitology. 2012;190(1):268–71. pmid:22677135
  34. 34. Courtenay O, Marinho-Junior JF, de Carvalho AW, de Carvalho F, de Paiva Cavalcante M, Brandão Filho SP, et al. Force of infection estimates in a cohort of free-ranging wild rodents naturally infected with Leishmania (Viannia) braziliensis, aetiological agent of human Tegumentary Leishmaniasis in the Americas. Submitted.
  35. 35. de Oliveira JA, Gonçalves PR. Genus Oxymycterus Waterhouse, 1837. In: Patton JL, Pardiñas UFJ, D’Elía G, editors. Mammals of South America, Volume 2 Rodents. Chicago: University of Chcago Press; 2015. p. 247–68.
  36. 36. Peçanha WT, Althoff SL, Galiano D, Quintela FM, Maestri R, Gonçalves GL, et al. Pleistocene climatic oscillations in Neotropical open areas: Refuge isolation in the rodent Oxymycterus nasutus endemic to grasslands. PLoS One. 2017;12(11):e0187329. pmid:29176839
  37. 37. Tonelli GB, Tanure A, Rego FD, Carvalho GML, Stumpp R, Ássimos GR, et al. Leishmania (Viannia) braziliensis infection in wild small mammals in ecotourism area of Brazil. PLOS One. 2017;12(12):e0190315. pmid:29284049
  38. 38. Forattini OP, Pattoli DBG, Rabello EX, Ferreira OA. Infecções naturais de mamíferos silvestres em área endêmica de leishmaniose tegumentar do Estado de São Paulo, Brasil. Revista de Saúde Pública. 1972;6(3):255–61.
  39. 39. Fernandez MS, Fraschina J, Acardi S, Liotta DJ, Lestani E, Giuliani M, et al. Assessment of the role of small mammals in the transmission cycle of tegumentary leishmaniasis and first report of natural infection with Leishmania braziliensis in two sigmodontines in northeastern Argentina. Parasitol Res. 2018;117(2):405–12. Epub 2017/12/22. pmid:29264720.
  40. 40. Vasconcelos IA, Vasconcelos AW, Fe Filho NM, Queiroz RG, Santana EW, Bozza M, et al. The identity of Leishmania isolated from sand flies and vertebrate hosts in a major focus of cutaneous leishmaniasis in Baturite, northeastern Brazil. Am J Trop Med Hyg. 1994;50(2):158–64. Epub 1994/02/01. pmid:8116807.
  41. 41. 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. Veterinary Parasitology. 2005;129(3–4):219–27. Epub 2005/04/23. pmid:15845276.
  42. 42. Pereira LOR, Moreira RB, de Oliveira MP, Reis SO, de Oliveira Neto MP, Pirmez C. Is Leishmania (Viannia) braziliensis parasite load associated with disease pathogenesis? Int J Infect Dis. 2017;57:132–7. Epub 2017/02/09. pmid:28167253.
  43. 43. Ferreira E, C., Cruz I, Canavate C, de Melo LA, Pereira AA, Madeira FA, et al. Mixed infection of Leishmania infantum and Leishmania braziliensis in rodents from endemic urban area of the New World. BMC Vet Res. 2015;11:71. Epub 2015/04/19. pmid:25890323; PubMed Central PMCID: PMC4374209.
  44. 44. Lima DO, Pinho GM, Fernandez FAS. Spatial patterns of the semi-aquatic rodent Nectomys squamipes in Atlantic forest streams. Journal of Natural History. 2016;50(7–8):497–511.
  45. 45. Ernest KA. Nectomys squamipes. Mammalian Species. 1986;(265):1–5.
  46. 46. Bonvicino CR, Weksler M. Genus Nectomys, Peters, 1861. In: Patton JL, Pardiñas UFJ, D’Elía G, editors. Mammals of South America, Volume 2 Rodents. 2. Chicago: University of Chcago Press; 2015. p. 369–77.
  47. 47. Mares MA, Ernest KA. Population and Community Ecology of Small Mammals in a Gallery Forest of Central Brazil. J Mammolgy. 1995;76(3):750–68.
  48. 48. D’Andrea PS, Gentile R, Maroja LS, Fernandes FA, Coura R, Cerqueira R. Small mammal populations of an agroecosystem in the Atlantic Forest domain, southeastern Brazil. Braz J Biol. 2007;67(1):179–86. Epub 2007/05/17. pmid:17505767.
  49. 49. Courtenay O, Carson C, Calvo-Bado L, Garcez LM, Quinnell RJ. Heterogeneities in Leishmania infantum infection: using skin parasite burdens to identify highly infectious dogs. PLoS Negl Trop Dis. 2014;8(1):e2583. pmid:24416460.
  50. 50. Marcelino AP, Ferreira EC, Avendanha JS, Costa CF, Chiarelli D, Almeida G, et al. Molecular detection of Leishmania braziliensis in Rattus norvegicus in an area endemic for cutaneous leishmaniasis in Brazil. Vet Parasitol. 2011;183(1–2):54–8. Epub 2011/07/20. pmid:21767914.
  51. 51. Sadlova J, Seblova V, Votypka J, Warburg A, Volf P. Xenodiagnosis of Leishmania donovani in BALB/c mice using Phlebotomus orientalis: a new laboratory model. Parasit Vectors. 2015;8:158. Epub 2015/04/17. pmid:25881258; PubMed Central PMCID: PMC4364506.
  52. 52. Guevara P, Rojas E, Gonzalez N, Scorza JV, Anez N, Valera M, et al. Presence of Leishmania braziliensis in blood samples from cured patients or at different stages of immunotherapy. Clin Diagn Lab Immunol. 1994;1(4):385–9. Epub 1994/07/01. pmid:8556473; PubMed Central PMCID: PMC368272.
  53. 53. de Oliveira CI, Brodskyn CI. The immunobiology of Leishmania braziliensis infection. Front Immunol. 2012;3:145. Epub 2012/06/16. pmid:22701117; PubMed Central PMCID: PMC3370302.
  54. 54. Canario A, Queiroz M, Cunha G, Cavalcante T, Riesz V, Sharma R, et al. Presence of parasite DNA in clinically unaffected nasal mucosa during cutaneous leishmaniasis caused by Leishmania (Viannia) braziliensis. Clin Microbiol Infect. 2019;25(4):515 e5- e7. Epub 2019/01/08. pmid:30616010.
  55. 55. Doehl JSP, Bright Z, Dey S, Davies H, Magson J, Brown N, et al. Skin parasite landscape determines host infectiousness in visceral leishmaniasis. Nat Commun. 2017;8(1):57. Epub 2017/07/07. pmid:28680146; PubMed Central PMCID: PMC5498584.
  56. 56. Lainson R, Shaw JJ. A brief history of the genus Leishmania (Protozoa: Kinetoplastida) in the Americas with particular reference to Amazonian Brazil. Ciência e Cultura. 1992;44:94–106.
  57. 57. Lainson R, Shaw JJ, Povoa M. The importance of edentates (sloths and anteaters) as primary reservoirs of Leishmania braziliensis guyanensis, causative agent of "pianbois" in north Brazil. Transactions of the Royal Society of Tropical Medicine and Hygiene. 1981;75(4):611–2. pmid:7324143.
  58. 58. Sadlova J, Vojtkova B, Becvar T, Lestinova T, Spitzova T, Bates P, et al. Host competence of the African rodents Arvicanthis neumanni, A. niloticus and Mastomys natalensis for Leishmania donovani from Ethiopia and L. (Mundinia) sp. from Ghana. Int J Parasitol Parasites Wildl. 2020;11:40–5. Epub 2019/12/28. pmid:31879594; PubMed Central PMCID: PMC6920283.
  59. 59. Roque AL, Jansen AM. Wild and synanthropic reservoirs of Leishmania species in the Americas. Int J Parasitol Parasites Wildl. 2014;3(3):251–62. pmid:25426421; PubMed Central PMCID: PMC4241529.
  60. 60. Oddone R, Schweynoch C, Schonian G, Dos Santos de Sousa C, Cupolillo E, Espinosa D, et al. Development of a multilocus microsatellite typing approach for discriminating strains of the Leishmania subgenus L. (Viannia). J Clin Microbiol. 2009;47(9):2818–25. pmid:19587302.
  61. 61. Figueiredo de Sá BSL, Rezende AM, Melo Neto OP, Brito MEF, Brandão Filho SP. Identification of divergent Leishmania (Viannia) braziliensis ecotypes derived from a geographically restricted area through whole genome analysis. PloS Negl Trop Dis. 2019;13(6):e0007382. Epub 2019/06/07. pmid:31170148; PubMed Central PMCID: PMC6581274.