Print ISSN: 2476-535X, Online ISSN: 2717-2910

Document Type : Original Article

Authors

1 Central Veterinary Research Institute, P.O. box 33980, Balmoral, Lusaka Zambia Chreso University, Faculty of Health Sciences, Makeni, Lusaka, Zambia

2 The University of Zambia, School of Veterinary Medicine, Lusaka Zambia

3 The University of Zambia, College of Veterinary Sciences, Lusaka Zambia

4 Sokoine University of Agriculture, College of Veterinary and Biomedical Sciences, Morogoro Tanzania

5 Chreso University, Faculty of Health Sciences, Makeni, Lusaka, Zambia

Abstract

Fleas (Siphonaptera) are important vectors of several animal and human disease pathogens, while rodents are considered as reservoirs of most pathogens, including Yersinia pestis Factors that influence the parasitism rate of fleas, ecological aspects that modulate their distribution, and host-flea relationship in Eastern Zambia remain unknown.  Furthermore, there is little information on the biodiversity and abundance of rodents and fleas in the study area. A total of 1212 mammals were sampled and examined. These included rodents (n=329), Crocidura (n=113), domestic pigs (n=254), small ruminants (n=346) and carnivores (n=168), and 1578 fleas, where five species were identified. There were nine genera and species of rodents with one genus of Crocidura captured. The results showed that 27(8.2%) and 19(5.8%) rodents and 8(7.0%) and 2(1.8%) Crocidura were positive for antibodies and pla gene for Y. pestis, respectively. Echidnophaga larina were the most mean abundant (MA=8.58), while Xenopsylla cheopis had the least mean abundant (MA=0.14), nevertheless it was the most infected with Y. pestis. Mastomys. natalensis was highest in plague positivity 31/56, followed by Crocidura spp 10/56 and Rattus rattus 6/56. The results indicated that three flea species were infected with Yersinia pestis. Shannon-Weiner (H) and dominance (D) indices of rodents were 1.5 and 0.2789, while the flea indices were 0.5310 and 0.8389, respectively. There was a strong association between richness of fleas and plague disease (p=0.01; x2=65.3). It’s established that rodents were more biodiversity than fleas while both were unevenly distributed. It’s recommended that control measures of fleas be intensified and sustained to lessen the spread of their associated diseases.

Keywords

Summary

Fleas (Siphonaptera) are important vectors of several animal and human disease pathogens, while rodents are considered as reservoirs of most pathogens, including Yersinia pestis Factors that influence the parasitism rate of fleas, ecological aspects that modulate their distribution, and host-flea relationship in Eastern Zambia remain unknown.  Furthermore, there is little information on the biodiversity and abundance of rodents and fleas in the study area. A total of 1212 mammals were sampled and examined. These included rodents (n=329), Crocidura (n=113), domestic pigs (n=254), small ruminants (n=346) and carnivores (n=168), and 1578 fleas, where five species were identified. There were nine genera and species of rodents with one genus of Crocidura captured. The results showed that 27(8.2%) and 19(5.8%) rodents and 8(7.0%) and 2(1.8%) Crocidura were positive for antibodies and pla gene for Y. pestis, respectively. Echidnophaga larina were the most mean abundant (MA=8.58), while Xenopsylla cheopis had the least mean abundant (MA=0.14), nevertheless it was the most infected with Y. pestis. Mastomys. natalensis was highest in plague positivity 31/56, followed by Crocidura spp 10/56 and Rattus rattus 6/56. The results indicated that three flea species were infected with Yersinia pestis. Shannon-Weiner (H) and dominance (D) indices of rodents were 1.5 and 0.2789, while the flea indices were 0.5310 and 0.8389, respectively. There was a strong association between richness of fleas and plague disease (p=0.01; x2=65.3). It’s established that rodents were more biodiversity than fleas while both were unevenly distributed. It’s recommended that control measures of fleas be intensified and sustained to lessen the spread of their associated diseases.

Keywords: Biodiversity Crocidura, Fleas, Plague, Rodents

 

 

Introduction

Fleas (Insecta: Siphonaptera) encompasses at least 2500 species with 15 subspecies globally (Whiting et al., 2008). They are small, laterally flattened and wingless, and highly specialized arthropods. They are of great importance as vectors of many viral, bacterial, and parasitic pathogens. Being hematophagous arthropods, fleas may transmit pathogens through soiled mouthparts (e.g., Ypestis, and viral pathogens), regurgitation of gut contents (e.g., Ypestis), and infectious saliva (e.g., R. felis in salivary glands) (Eisen et al., 2012). Also, transmission can occur through contaminated feces (e.g., Rickettsia typhi, Bartonella,hanselae) (Amatre et al., 2009) (Bitam et al., 2010). Fleas have complex life cycles, which include the production of non-parasitic omnivorous larvae and obligate-hematophagous imagoes (Sanchez and Lareschi, 2019). Once the flea emerges from the cocoon, it immediately seeks a potential host to get a blood meal. Most fleas prefer the host habitat (nest fleas), while others stay more or less permanently on the body of the host itself (body fleas) (Bitam et al., 2010).

Fleas have tremendous medical and economic importance because they are vectors of several causative agents of diseases in animals and humans.  Such diseases as  the cat scratch fever  caused by B. henselae, Q fever (Coxiella burnetii), murine typhus (R. typhii), flea-borne spotted fever (R. felis),and bubonic plague (Y. pestis) (Bitam et al., 2010). Due to the potential vector capacity of fleas, their distribution and biodiversity may have essential consequences for host survival and disease dynamics (Foottit and Galloway, 2018). Arthropod-transmitted diseases represented approximately 17% of all infectious diseases globally and are influenced by a complex and dynamic ecosystem that involves vectors, hosts, and infectious agents, as well as environmental factors(Young et al., 2015). Fleas are able  to bite animals, suck blood and mechanically or biologically transmit disease pathogens such as Yersinia pestis to susceptible hosts, including humans and rodents (Ago et al., 2015). The abundance of fleas in an area may suggest an increased number or frequency of flea-borne diseases. This happens when fleas are numerous and consequently look for animals and humans to suck blood and subsequently transmit disease-causing pathogens (Nziza et al., 2019; Miarinjara and Boyer, 2016).

Rodents are small mammals of the order Rodentia, which comprises more than 2000 species and approximately 30 families. There are diverse biological and ecological differences among rodents in shape, size, weight, and habitat. The smallest rodents (Mus minutoide) may weigh 5 grams, while the largest (American cabybard) weighs more than 70 kg. Rodent's ecological requirements vary widely, which include domestic and semi-domestic field and forest habitats. These animals are potential reservoirs of infectious diseases of humans and other mammals and are also suitable hosts for flea and other ectoparasites (Esfandiari et al., 2017).

Shrews are insectivores and able to harbor disease pathogens, such as bacteria, virus, and protozoa. These disease pathogens could be transmitted to a suitable host and cause disease through the vector such as fleas (Moore et al., 2015).

The abundance of rodents and fleas depends on climatic and environmental conditions. Climate changes caused by global warming and human intervention have contributed to vicissitudes in the biological parameters and distribution ranges of vector/fleas and their hosts (Kausrud et al., 2007). These conditions are likely to affect the biodiversity and abundance or dominance of both rodents and their flea ectoparasites. At times of heavy rainfall and availability of food crops, there are corresponding large numbers of rodents and their fleas, contrary to times of adverse climate when the rodents and their fleas diminish in number. Rodents have the potential to transmit sylvatic plague and other zoonotic diseases directly to humans and other animals. Such diseases can be transmitted by contamination of foods with rodent fecal materials or urine, direct contacts with infected animals, droplets, or bites by appropriate insect vectors such as fleas in case of plague disease (Backhans and Fellström, 2012). These vectors then transmit the pathogens from one reservoir host to another susceptible animal, including humans, small carnivores, and other animals.  In view of the aforegoing description of the relationship between rodents and fleas regarding plague epidemiology, it was felt desirable to investigate and establish the current status of biodiversity and abundance of these creatures in relation to plague epidemiology; therefore the aim of the study was to accomplish this desirability in the eastern region of Zambia.

 

Materials and methods

Area and time of study

This study took place in two districts, Nyimba and Sinda, both in the eastern province of Zambia from 2013 to 2017. The choice of these areas was based on their recent history of plague outbreaks, whereby the last human plague was detected in Nyimba and Sinda districts in 2015 and 2008, respectively. The eastern Zambia districts experience high rainfall between January and March of each year. The communities grow crops such as maize, groundnuts, cassava, which are all suitable for rodents and other animals.

Study design

A cross-sectional study design was conducted, where the abundance and biodiversity of the rodents, shrews, and flea ectoparasites were examined.

Sample collection

Domestic animals

Villages in the study area were randomly selected in which 20 or more domestic animals of the selected species, comprising pigs, dogs, cats, goats, and sheep were sampled. From each selected client, a verbal consent was acquired (as per approved Plague Research policy provided in Assurance No. FWA0000338).  Each selected animal was assigned a unique identification number. The first animal was randomly selected, followed by a systematic technique (Systematic Random Sampling). The selected animal was immobilized and rested on the white plastic sheet for blood collection and inspection of fleas and other ectoparasites. The animal was groomed with cotton wool containing 90% diethyl ether to sedate the ectoparasites. The appropriate brush was used to scrub the animal to detach fleas and other ectoparasites from its fur and skin. A pair of forcepswas used to gently removed fleas that fell onto the white sheet and those which remained attached on the animal’s fur/skin into vials with 70% ethanol. The study did not include animals that came from other villages in the past six months (Nyirenda et al., 2017).

Rodents and shrews trapping

The study villages were selected randomly and were alienated into six arbitrary sectors from where three zones were randomly nominated for trapping rodents and shrews. Sherman’s live traps (50×65×157 mm) were baited with peanut butter mixed with soya beans flour and were placed at a distance of 10m apart in nearby bushes overnight. Wire cage traps (145×100×230 mm) (Hoga-lab, Kyoto, Japan), baited with fish, Stolothrissa tanganicae (Kapenta), and ripe tomatoes, were set in chosen houses in the zones. Inspection of traps was conducted the following morning. The captured animals were taken for flea ectoparasite inspection and collection in a mobile laboratory installed at a distant place away from the community to collect sera, organs, and fleas. Other processes included identification and enumeration of rodents and fleas were carried out in such a laboratory. The ensnaring process was continuous for three consecutive days in the same area as previously described (Kilonzo, 1976). Sera and organs were kept at -20oC until needed for ELISA and PCR, respectively (Nyirenda et al., 2018).

Collection of flea ectoparasites from rodents and shrews

Fleas and other ectoparasites from the rodents and shrews were collected by introducing the latter into a plastic bag with a soaked cotton wool, in 90% diethyl ether to sedate both the animal and fleas. Once the animals were sedated, they were placed into a silver basin and brushed thoroughly with a hard toothbrush to disengage the ectoparasites. Fleas that dropped into the basin were gently picked using either a pair of forceps or fine camel brush into vials with 70% ethanol as a preservative. Other ectoparasites were also collected into separate vials with the same concentration of ethanol (Nyirenda et al., 2018).

Rodent identification

All the captured rodents were identified to genus and/or specie level using the key features described by Kingdon (Kingdon, 1974). In addition, the body and tail lengths were measured and recorded. The carcasses were then counted, and numbers of each species were established.

Flea processing and identification

The fleas were preliminary identified and pooled (1–5) according to their species, location of the host, and species of its host. One to two (1–2) fleas were removed and processed from each pool to confirm the identity. The fleas were processed and identified to specie level using the method described by Kilonzo (Kilonzo, 1999), which basically involve drying the insects on the No. 1 filter paper, boiling them for 10 minutes in heat resistant tubes containing 10% Sodium Hydroxide. The fleas were then put in cold water for one hour and into acidified water (water + glacial acetic acid of equal volume) for 30 minutes to neutralize Sodium Hydroxide. After this, fleas were dehydrated through increasing strengths of ethanol (from 50%, 70%, 95%, and absolute alcohol) for one hour at each stage. After this process, the fleas were transferred into the vial  with  clove oil and incubated at room temperature overnight or until they were clear and transparent, after which they were mounted with Dibutylphthalate Polystyrene Xylene (DPX) on the glass slide with a coverslip for examination under the light microscope using dry objectives (x4, x10 or x40) magnification for detection of common taxonomic features (Kilonzo, 1999).

Identification of the fleas was based on  main flea features such as pronotal combs, genal combs, and shape of the head and reproductive organs (spermathecae in females and penis plates in males (Nyirenda et al., 2018; Kilonzo, 1999; Pratt and Stojanovich, 1966; Traub, 1962).

DNA extraction from organs and fleas

The DNA from tissues of rodents was extracted using the DNA extraction kit (ZR genomic DNA-Tissue Mini-Prep Catalog No. D3051) following the manufacturer’s instructions (Zymo Research Irvine, CA, USA). In the case of fleas, these were identified before preparing them for DNA extraction using the heat treatment method as described elsewhere (Nyirenda et al., 2017).

 

PCR Technique

PCR was performed using Phusion flash high fidelity master mix (Finnzymes Oy, Finland) in a highly PCR specialized laboratory. The primers Yp pla1 (5'TGC TTT ATG ACG CAG AAA CAG G3') as the forward primer and Yp pla2 (5'CTG TAG CTG TCC AAC TGA AAC G3') as the reverse primer, that amplifies a 344-bp region spanning residues 425 to 769 of the plasminogen activator gene, were used (Sodeinde and Goguen, 1989). The PCR reactions were performed as described elsewhere (Nyirenda et al., 2018) (Nyirenda et al., 2017). The positive and negative controls, which came together with the kit, were also included in the test.

ELISA testing

The stored sera were removed from the freezer (−20 °C) and left at room temperature to thaw and processed using the protocol as previously described (Chu, 2000).

Data analysis

Data were entered in Microsoft Excel software and analyzed using Epi infoTM 7.0.8.0, a computer statistical package from the Centre for Disease Control and Prevention (CDC), where confidence intervals (CI) and positive percentages were generated. The p-value and the chi-square (x2) were also calculated.

Flea diversity was estimated by calculating the Specific richness (S=number of flea species) and Shannon diversity species index (H = -Σpi ln pi), where pi =proportion of each flea species in the sample). Flea evenness was estimated by calculating the Simpson dominance index (D = Σp2i), where pi is the proportion of individuals found in species ‘i' as calculated elsewhere (Sanchez and Lareschi, 2019; Ralaizafisoloarivony et al., 2014).

For each host species, the following indices and parameters were calculated: flea species richness (S=number of flea species), Simpson index (Mean Abundance (MA) = total number of individuals’ of a parasite species on a host/total number of host species, including infested and non-infested species), and prevalence ([P = number of infested animals with one or more individuals of a parasite species/total number of examined mammals for that parasite species] multiplied by 100).

 

Results

In this study, 1212 mammals were examined for flea vectors, which comprised rodents (n=329), shrews (n=113), domestic pigs (n=254), small domestic ruminants (n=364), and domestic carnivores (n=168) were examined for fleas. A total of 1578 fleas were collected from such animals. The rodent population composed of nine genera and species, while all the shrews belonged to the genera of Crocidura spp. (n=113) The rodent species included: Mastomys natalensis (n=189), Rattus rattus (n=60), Saccostomus spp (pouched mouse (n=43), Tatera spp (n=2), Graphiurus spp (n=5), Mus spp (n=2), Acomys spp (n=1), Gerbillurus (Gerbil) spp (n=22) and Steastomys parvus (fat mouse) (n=5). On ELISA and PCR, 27 (8.2%) and 19 (5.8%) were positive for Yersinia pestis antibodies and pla gene, respectively, while 8 (7.0%) and 2 (1.8%) shrews were also positive for Yersinia pestis antibodies and pla genes, respectively. Five species of fleas collected from the animals in the study area were; X. cheopis (n=59), C. canis (n=64), E. gallinacea (n=389), E. larina (n=1064),and C. felis (n=2).The results also showed E. gallinacea had the highest mean abundance (MA) of the (8.58), while X. cheopis had the lowest (MA=0.14)in both districts (Table 1). It was revealed that Mastomys natalensis (31/56) was the highest-ranking rodent species carrying Y. pestis. Mus spp, Acomy spp, and Steatomys parvus were negative for Y. pestis on both tests. Statistical analysis presented no significant statistical difference between the rodents and plague disease (p = 0.347; x2=10) (Table 2).

The study further revealed that three species, namely X. cheopis, E. gallinacean, and C. canis, were infected with the Y. pestis bacterium. The plague positive fleas were collected from rodents, pigs, and goats, respectively. Statistical analysis demonstrated a significant statistical difference between richness of fleas and plague disease (p = 0.01; x2=65) (Table 3). Our results also showed that the biodiversity Shannon-Weiner (H) and dominance (D) of rodents were 1.5 and 0.2789 respectively in the study area (Table 4), while the biodiversity Shannon-Weiner (H) of the fleas showed 0.5310 and the dominance (D) showed 0.8389 (Table 5).

Discussion

Rodents

The current observations that M. natalensis is the most abundant rodent species in the study area, followed by R. rattus, and that the two species had the highest infection rates with Yersinia pestis (Table 2).  This suggests that the plague organisms are endemically in the area and that the rodent species in question are the most suitable reservoirs of the plague disease in the area. These observations are consistent with earlier findings described by Kilonzo (Kilonzo, 1999). In view of their ecological nature and close association with humans in houses and the fields, these rodents could easily transmit Y. pestis to humans and other animals as they were diverse and abundant in the study area. The observation that the two species are abundant and the fact that they are domestic and/or semi-domestic suggest that they can destroy substantial food and cash crops, thus causing hunger and poverty in the area. Furthermore, the current results indicate that two rodent species, as well as shrews were plague positive in both Sinda and Nyimba districts. Detection of Yersinia pestis antibodies and the pla gene in the Crocidura spp suggest that these small mammals were in contact with the bacterium and could serve as the reservoir host for the disease in the area(Moore et al., 2015). Despite the detection of plague in these animals in both ELISA and PCR, there was no outbreak during our study. This is at least partially consistent with observation elsewhere that rodents are merely reservoirs of Yersinia pestis, and the presence of suitable vectors alone could not stimulate the outbreak of the disease.

 The general picture is that rodents were biodiverse in the area, and they were unevenly distributed, thus resulting in sporadic numbers in some areas. This scenario is probably attributable to climate and environmental changes(Morand, 2011)(Eisen et al., 2015). In view of the current observations, establishment and adherence to regular surveillance services for plague and other flea/rodent zoonoses is desirable.

Fleas

The observation in the current study that Echidnophaga larina collected from domestic pigs were the most abundant fleas species in the study area, followed by Echidnophaga gallinacea collected from the same animals, probably suggest that the animals (pigs) were not regularly dipped in or spray/dusted with appropriate pesticides and their houses are unhygienic. The presence of Yersinia pestis antigen in 1% E. gallinacea collected from pigs and 31.3% of C. canis collected from goats in Nyimba districts suggests that the two flea species are suitable reservoirs/carriers of pathogen. However, the two flea species are potential, but not efficient vectors of the disease. On the other hand, the presence of Y. pestis in 5.1% X. cheopis collected from rodents trapped in Nyimba district suggests the potential plague outbreak of human plague in the area if and when conditions become favorable, since the flea is known to be the efficient vector of the disease globally (Miarinjara and Boyer, 2016).

Furthermore, the current observation showed that the fleas were neither biodiverse nor evenly distributed, thus suggesting that they are seasonal as reported elsewhere in Argentina, Tanzania, and Uganda (Sanchez and Lareschi, 2019;Young et al., 2015; Eisen et al., 2012; Kilonzo et al. 2006; Njunwa et al., 1989). During the hot season, fleas can breed more efficiently, and their survival rate is higher than in other seasons. Flea density dynamic is high and is influenced by several factors, including the host, altitude, temperature, and relative humidity (RH) (Eisen et al., 2015b). The latter is the most appropriate condition for vectors to breed and transmit Y. pestis and thus facilitate disease outbreaks(Eisen et al., 2015a).

Besides transmitting disease pathogens, fleas also inflict and injure the skin of their host, thus  allowing opportunistic bacteria to infect the host open wounds (Krasnov et al., 2019). The current observation of the large number of fleas on domestic pigs, therefore, suggests that such animals are at risk of getting severe skin problems or infection caused by skin-penetrating opportunistic organisms. The low flea index of X. cheopis (efficient vectors) on rodents (suitable vectors) suggests a low plague outbreak at the moment. However, monitoring of these ectoparasites in the disease quiescent period is essential to assess seasonal increase and subsequently predict the plague outbreak in the area (Ralaizafisoloarivony et al., 2014). The significant association (p=0.01) between flea numbers and plague is a further indication that plague depends on the bites from infected fleas.

Conclusion and recommendations

It can be justifiably concluded from the observations in this study that rodents are more bio-diverse than fleas while both of them are unevenly distributed in the study area. It can also be concluded that Echidnophaga larina was the most abundant domestic animal flea species in the area and it is most hosted by domestic pigs. Likewise, it is conclusive from the study a well-known efficient vector of plague was the only rodent flea species in the area and its population density (flea per rodent) was too low to maintain transmission of plague at the time of the study. It can further be concluded from the study that a thorough understanding of the flea-host association and rodent and flea distribution and diversity may be useful in epidemiological studies of plague disease and could provide a reliable base for prediction and surveillance of plague and other related re-emerging zoonotic diseases. It is, therefore, recommended that management of flea ectoparasites should be strengthened and surveillance services are maintained to forecast and consequently prevent or minimize the outbreaks of such diseases.

 

Acknowledgment

We would like to thank the Department of veterinary services in the Ministry of Fisheries and Livestock for allowing us to sample the livestock in the study area. Special thanks go to Mr. Siwale G of Nyanje mission hospital, Veterinary officers at Sinda and Nyimba districts to organize the local community and livestock farmers during the time of sampling.

Funding

The authors received no specific funding for this work

Ethical Approval

Ethical approval to conduct this study was sought the Biomedical Research Ethics Committee (BREC), Zambia (Assurance No. is FWA00000338). Verbal consent to sample the animals was granted by the livestock farmers from the study area.

 

Conflict of Interest Statement

The authors have declared that no competing interests exist.

Ago Y., Rasmussen S., Allentoft M.E., Nielsen K., Nielsen R., Kristiansen K., Willerslev E. and et al. (2015). Early Divergent Strains of Yersinia Pestis in Eurasia Article Early Divergent Strains of Yersinia Pestis in Eurasia 5000 Years Ago. Cell, 163, pp. 571–82.
Amatre G., Babi N., Enscore R.E., Ogen-Odoi A., Atiku L. Akol A., Gage K.L. and Eisen R.J. (2009). Flea Diversity and Infestation Prevalence on Rodents in a Plague-Endemic Region of Uganda. The American Journal of Tropical Medicine and Hygiene, 81 (4), pp. 718–24.
Backhans A. and Fellström C. (2012). Rodents on Pig and Chicken Farms - a Potential Threat to Human and Animal Health. Infection Ecology & Epidemiology, 2, 17093.
Bitam I., Dittmar K., Parola P., Whiting M.F. and Raoult D. (2010). Fleas and Flea-Borne Diseases. International Journal of Infectious Diseases (Official Publication of the International Society for Infectious Diseases), 14 (8), pp. e667-e676
Chu M.C. (2000). Laboratory Manual of Plague Diagnostics. Geneva: . US Centers for Disease Control and Prevention and World Health Organization.
Eisen, R.J, Borchert J.N, Mpanga J.T, Atiku L.A, MacMillan K, Boegler K.A, Montenieri J.A, Monaghan A and Gage K.L. (2012). Flea Diversity as an Element for Persistence of Plague Bacteria in an East African Plague Focus. PLoS ONE, 7 (4), pp. 1–8.
Eisen, R.J, Dennis D.T and Gage K.L. (2015a). The Role of Early-Phase Transmission in the Spread of Yersinia Pestis. Journal of Medical Entomology, 52 (6), pp. 1183–92.
Esfandiari B., Nahrevanian H., Khaki P., Esfandiari B., Nahrevanian M.R.P., Gouya M.M., Hanifi H., Khaki P., Mostafavi E., Darvish J. (2017). “Epidemiological Distribution of Rodents as Potent Reservoirs for Infectious Diseases in the Provinces of Mazandaran, Gilan and Golestan, Northern Iran.” Infectious Disease Reports, 9 (6900), pp. 62–65.
Foottit R.G, Adler P.H, and Galloway T.D. (2018). Biodiversity of Ectoparasites: Lice (Phthiraptera) and Fleas (Siphonaptera) In: Insect Biodiversity. Science and Society II, John Wiley & Sons.
Kausrud K.L., Viljugrein H., Frigessi A., Begon M., Davis S., Leirs H., Vladimir D. and Stenseth N.C. (2007). Climatically Driven Synchrony of Gerbil Populations Allows Large-Scale Plague Outbreaks. Proceedings. Biological Sciences / The Royal Society, 274, pp. 1963–69.
Kilonzo, B.S. (1999). Basic Procedures for Collecting, Processing and Identification of Common Fleas. Morogoro, Tanzania.
Kilonzo B.S, Mbise T.J, Mwalimu D.C. and Kindamba L. (2006). Observations on the Endemicity of Plague in Karatu and Ngorongoro, Northern Tanzania. Tanzania Health Research Bulletin 8 (1), pp. 1–6.
Kilonzo B.S. (1976). A Survey of Rodents and Their Flea Ectoparasites in North-Eastern Tanzania. East African Journal of Medical Research, 3, pp. 117–126.
Kilonzo B.S. (1999). Plague Epidemiology and Control in Eastern and Southern Africa during the Period 1978 to 1997. The Central African Journal of Medicine, 35, pp. 70–76.
Kingdon, Jonathan. (1974). East African Mammals: An Atlas of Evolution in Africa. (Hares and Rodents). London: Academic press.
Krasnov B.R, Shenbrot G.I, Warburton E.M, Van Der Mescht L., Surkova E.N., Medvedev S.G., Pechnikova N., Ermolova N., Kotti B.K. and Khokhlova I.S. (2019). Species and Site Contributions to β-Diversity in Fleas Parasitic on the Palearctic Small Mammals: Ecology, Geography and Host Species Composition Matter the Most. Parasitology, 146 (5), pp. 653–61.
Miarinjara A and Boyer S. (2016). Current Perspectives on Plague Vector Control in Madagascar: Susceptibility Status of Xenopsylla Cheopis to 12 Insecticides. PLoS Neglected Tropical Diseases, 10 (2), e0004414.
Moore S.M, Monaghan A, Borchert J.N, Mpanga J.T, Atiku L.A, Boegler K.A, Montenieri J, Macmillan K, Gage K.L and Eisen R.J. (2015). Seasonal Fluctuations of Small Mammal and Flea Communities in a Ugandan Plague Focus: Evidence to Implicate Arvicanthis Niloticus and Crocidura Spp. as Key Hosts in Yersinia Pestis Transmission. Parasites and Vectors, 8 (11), pp. 1–15.
Morand S. (2011). Infectious Diseases, Biodiversity and Global Changes: How the Biodiversity Sciences May Help, In: The Importance of Biological Interactions in the Study of Biodiversity. pp. 231–54.
Njunwa K.J., Mwaiko G.L., Kilonzo B.S. and Mhina J.I.K. (1999). Seasonal Patterns of Rodents, Fleas and Plague Status in the Western Usambara Mountains, Tanzania. Medical and Veterinary Entomology, 3 (1), pp. 17–22.
Nyirenda S.S., Hang’ombe B.M., Kilonzo B.S., Kangwa H.L., Mulenga E. and Moonga L. (2017). Potential roles of pigs, small ruminants, rodents, and their flea vectors in plague epidemiology in Sinda district, eastern Zambia. Journal of Medical entomology, 54(3), pp. 719-725
Nyirenda S.S., Hang’ombe B.M., Mulenga E. and Kilonzo B.S. (2017). Serological and PCR Investigation of Yersinia Pestis in Potential Reservoir Hosts from a Plague Outbreak Focus in Zambia. BMC Research Notes, 10 (1), 345.  
Nyirenda, S.S., Hang’ombe B.M., Simulundu E., Mulenga E., Moonga L., Machang’u R.S., Misinzo G. and Kilonzo B.S. (2018). Molecular Epidemiological Investigations of Plague in Eastern Province of Zambia. BMC Microbiology, 18, 2.
Nziza J., Tumushime J.C., Cranfield M., Ntwari A.E., Modrý D., Mudakikwa A., Šlapeta J. (2019). Fleas from Domestic Dogs and Rodents in Rwanda Carry Rickettsia Asembonensis and Bartonella Tribocorum. Medical and Veterinary Entomology, 33 (1), pp. 177–84.
Pratt H.D. and Stojanovich C.H.J. (1966). Illustrated Key to Species Found during Plague Investigations, In: US Department of Health, Education and Welfare Ed. Atlanta, USA: Communicable Disease Center. pp. 171–74.
Ralaizafisoloarivony N.A., Kimaro D.N., Kihupi N.I., Mulungu L.S., Leirs H., Msanya B.M., Jozef A., Deckers J.A. and Gulinck H. (2014). Small Mammals Distribution and Diversity in a Plague Endemic Area in West Usambara Mountains, Tanzania. Tanzania Journal of Health Research, 16 (3), pp. 1–9.
Sanchez J. and Lareschi M. (2019). Diversity, Distribution and Parasitism Rates of Fleas (Insecta: Siphonaptera) on Sigmodontine Rodents (Cricetidae) from Argentinian Patagonia. Bulletin of Entomological Research, 109 (1), pp. 72–83.
Sodeinde O.A. and Goguen J.D. (1989). Nucleotide Sequence of the Plasminogen Activator Gene of Yersinia Pestis: Relationship to OmpT of Escherichia Coli and Gene E of Salmonella Typhimurium. Infection and Immunity, 57 (5), pp. 1517–23.
Traub R. (1962). Rothschild Collection of Fleas. Nature, 196 (4852), 304.
Whiting M.F., Whiting A.S., Hastriter M.W. and Dittmar K. (2008). A Molecular Phylogeny of Fleas (Insecta: Siphonaptera): Origins and Host Associations. Cladistics, 24 (5), pp. 677–707.
Young H.S., Dirzo R., McCauley D.J., Agwanda B., Cattaneo L., Dittmar K., Eckerlin R.P., et al. (2015). Drivers of Intensity and Prevalence of Flea Parasitism on Small Mammals in East African Savanna Ecosystems. Journal of Parasitology, 101 (3), pp. 327–35.