Next Article in Journal
Depression, Anxiety and Associated Factors among Frontline Hospital Healthcare Workers in the Fourth Wave of COVID-19: Empirical Findings from Vietnam
Next Article in Special Issue
The Geographical Distribution and Influencing Factors of COVID-19 in China
Previous Article in Journal
Pilot Malacology Surveys for the Intermediate Hosts of Schistosomiasis in Rural and Semi-Urban Areas of the Moyen-Ogooué Province, Gabon
Previous Article in Special Issue
Larvicidal Activities against Aedes aegypti of Supernatant and Pellet Fractions from Cultured Bacillus spp. Isolated from Amazonian Microenvironments
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Environmental and Household-Based Spatial Risks for Tungiasis in an Endemic Area of Coastal Kenya

1
Graduate School of Biomedical Sciences, Nagasaki University, 1-12-4 Sakamoto, Nagasaki-shi 852-8523, Nagasaki, Japan
2
Department of Eco-Epidemiology, Institute of Tropical Medicine, Nagasaki University, 1-12-4 Sakamoto, Nagasaki-shi 852-8523, Nagasaki, Japan
3
Nagasaki University Institute of Tropical Medicine-Kenya Medical Research Institute (NUITM-KEMRI) Project, C/O Centre for Microbiology Research, KEMRI, Nairobi P.O. Box 19993-00202, Kenya
4
Social Environment and Health, Survey Research Center, Institute for Social Research, University of Michigan, Ann Arbor, MI 48109, USA
5
Department of Epidemiology, School of Public Health, University of Michigan, Ann Arbor, MI 48109, USA
6
Centre of Public Health Research, Kenya Medical Research Institute (KEMRI), Off Mbagathi Road, Nairobi P.O. Box 54840-00200, Kenya
7
Division of Health and Safety Sciences Education, Department of Educational Collaboration, Osaka Kyoiku University, 4-698-1 Asahigaoka, Kashiwara-shi 582-8582, Osaka, Japan
*
Author to whom correspondence should be addressed.
Trop. Med. Infect. Dis. 2022, 7(1), 2; https://doi.org/10.3390/tropicalmed7010002
Submission received: 11 November 2021 / Revised: 16 December 2021 / Accepted: 17 December 2021 / Published: 23 December 2021
(This article belongs to the Special Issue Spatial Epidemiology of Infectious Diseases)

Abstract

:
Tungiasis is a cutaneous parasitosis caused by an embedded female sand flea. The distribution of cases can be spatially heterogeneous even in areas with similar risk profiles. This study assesses household and remotely sensed environmental factors that contribute to the geographic distribution of tungiasis cases in a rural area along the Southern Kenyan Coast. Data on household tungiasis case status, demographic and socioeconomic information, and geographic locations were recorded during regular survey activities of the Health and Demographic Surveillance System, mainly during 2011. Data were joined with other spatial data sources using latitude/longitude coordinates. Generalized additive models were used to predict and visualize spatial risks for tungiasis. The household-level prevalence of tungiasis was 3.4% (272/7925). There was a 1.1% (461/41,135) prevalence of infection among all participants. A significant spatial variability was observed in the unadjusted model (p-value < 0.001). The number of children per household, earthen floor, organic roof, elevation, aluminum content in the soil, and distance to the nearest animal reserve attenuated the odds ratios and partially explained the spatial variation of tungiasis. Spatial heterogeneity in tungiasis risk remained even after a factor adjustment. This suggests that there are possible unmeasured factors associated with the complex ecology of sand fleas that may contribute to the disease’s uneven distribution.

1. Introduction

Tungiasis is a zoonotic, cutaneous parasitosis caused by the embedding of female sand fleas (Tunga spp.) into the upper strata of the skin [1]. Tungiasis is common to poor communities in sub-Saharan Africa, South America, and the Caribbean [2]. The female flea burrows headfirst into the upper strata of the skin. After copulation, the female produces eggs, which cause the flea to expand to several times its original size, compressing and causing stress to the surrounding skin tissues. Eggs are expelled through the exposed posterior, fall to the ground, and mature in sands and soils [3,4].
Tungiasis represents the classic neglected tropical disease: it afflicts the poorest people mainly in very rural communities in tropical countries, and there is little research and almost no investment in efforts to prevent, control, and/or treat the disease [1]. Few studies have reported on affected populations globally, though over 20 million in the Americas and 668 million people in sub-Saharan Africa are estimated to be at risk of tungiasis [5,6].
To prevent and control tungiasis, it is essential to elucidate specific risk factors that raise or mitigate infection. Several environmental, demographic, and behavioral factors have been associated with disease incidence. These include age groups, specifically children and the elderly [7,8,9,10,11]; low socioeconomic status associated with inadequate housing [7,9,11,12,13,14,15,16]; low affordability, access to, and knowledge of hygiene and health behaviors [11,12,14,15]; and the lack of clean water and sanitation at the household level [7,16]. Relevant to the zoonotic aspects of tungiasis, the presence of animals of various species in and around household compounds was also reported to raise the risk of human tungiasis [9,11,12,13,14,17]. These risk factors are also common to poor communities in tropical developing countries. More detailed information on how risk for tungiasis is distributed within communities is needed for the efficient development and implementation of preventive measures.
The geographical distribution of tungiasis can vary widely, even between neighboring areas that share similar risk profiles [1,7,16]. The environmental suitability of tungiasis varies across the sub-Saharan African continent, and comprises a complex range of topographic, environmental, and human-influenced factors (e.g., livestock density) [6]. In particular, soil suitability for sand fleas may influence their distribution, as the soil is the site of their off-host development [3,18]. Therefore, we should simultaneously consider environmental and human risk factors in the analysis for tungiasis infection and the recorded geographic distribution of tungiasis infection when assessing areas of risk. To date, however, there have been few epidemiological studies of tungiasis that consider how potential environmental factors and risk factors related to individuals and households determine the spatial distributions of the disease using on the ground, house to house prevalence surveys in sub-Saharan Africa.
This research tested three main hypotheses. First, we suspect that tungiasis risk will be heterogeneously distributed, i.e., households at a high risk for tungiasis will be proximal to other households at risk. Second, we suspect that environmental factors, such as soil type, wetness, and proximity to a wildlife reserve, will explain the spatial variation of tungiasis cases in this area. Third, we speculate that a combination of environmental and household factors will work in concert to raise or lower risk.
This study aims to address spatial risks that cause a spatially heterogeneous distribution of tungiasis by considering risk factors related to demographics and socio-economy, as well as spatial ecological factors, in a tungiasis endemic area using a comprehensive, household-level tungiasis survey in a resource-poor, rural area of East Africa.

2. Materials and Methods

2.1. Study Population and Area

Tungiasis infection data and household information were collected by the Health and Demographic Surveillance System (HDSS) managed by the Nagasaki University Institute of Tropical Medicine and Kenya Medical Research Institute (NUITM-KEMRI) [19]. The HDSS conducts survey activities in Kwale County, along the southern coastal region of Kenya on the Indian Ocean. Kwale is an arid area that belongs to the savanna climatic zone. The HDSS in the study area was set in 2010 to monitor the population dynamics of around 40,000 residents (10,000 households), such as births, deaths, and migrations, in the coverage area of approximately 385 km2, where the geographical location ranges between 5°42′52″ S, 39°14′35″ E and 5°55′31″ S, 39°31′6″ E along the northern side of the Shimba Hills National Reserve and surrounding the Mwaluganje Elephant Sanctuary (Figure 1) [19]. All residents enrolled in the HDSS who were present at home at the time of data collection and who consented to participate were included in the analysis.
  • Inclusion and exclusion criteria
All permanent residents of the HDSS area over the age of nine months were considered eligible for inclusion in the study. Non-permanent residents, defined as people not intending to stay in the HDSS area for more than three months, are not normally included in routine data collections and thus were not included in the study.

2.2. Tungiasis Infection and Household Information

We used survey information on tungiasis infection collected between March and December 2011 during the HDSS follow-up survey. Survey teams were trained by clinical officers to identify signs and symptoms of tungiasis. At the time of data collection, cases were identified by first asking individuals if they were experiencing a tungiasis infection. If the individual reported that he or she was infected, either through self-report or through a household representative, the case status was confirmed by visual inspection. Household members who were not reported as being infected were considered negative and recorded as such. No attempt was made to grade lesions using an assessment tool, such as the Fortaleza classification [4], and to identify the species of Tunga present in participants’ lesions. Our ethical approval did not allow for the invasive extraction of gravid female sand fleas from the skin, which is necessary for species identification [20].
In the case of data duplication as a result of visiting the household more than once, information from the first time the individual was surveyed was retained for analysis. This study considered household-level infection, rather than individual infection, when assessing spatial risks. We recorded households as “tungiasis present” when at least one individual in the home was positive for tungiasis. This binary measure was used as an outcome in all analyses for this research.
Household-level explanatory variables included latitude/longitude coordinates, wall, floor, and roof materials, and toilets. These were available from the HDSS baseline survey conducted in July 2010 and from follow-up surveys through 2012. Households with missing information were excluded from the analysis.
Household demographic information composition was aggregated to the household level. We calculated the numbers of children under 15 years of age and elderly over 60 years of age per household to assign the age structure of household members. For those who belonged to more than one household (e.g., polygamy or second house owner), we assigned a value divided by the number of households they belonged to (e.g., 0.5 people if a person belongs to two households) as the number of residents per household (number of people per household). After conversion, households with more than 19 people (upper quartile plus three quartile ranges) were excluded as outliers.

2.3. Ecological Factors

We included several variables related to the ecology of Tunga spp.; specifically, we chose variables relevant to off-host stage development based on a review of relevant literature [3,6,12,21]. The data comprised remote sensing and freely available data sets. See Table S1 for details. Pointwise values were extracted using latitude/longitude coordinates of each household. All raster to point data extraction was conducted on QGIS 3.14 [22]. The following is a brief description of the environmental data used in this study.

2.3.1. Normalized Difference Vegetation Index (NDVI) and Yearly Land Cover in 2011

NDVI data were obtained from the United States National Aeronautics and Space Administration Land Processes Distributed Active Archive Center (NASA LP DAAC) through the Earth Engine Data Catalog [23]. The median NDVI of each pixel was calculated from merged layers of 16-day NDVI with masking low-quality pixels. Then, NDVI was integerized using the following equation: NDVI   ( 8   bit ) = ( NDVI + 1 ) × 100 . The generated value ranges from 0–200 instead of the original range, which is from −1 to 1.

2.3.2. Topographic Wetness Index (TWI)

TWI is a proxy for soil wetness index, representing where water will accumulate, taking elevation differences into account [24,25]. Source data to compute TWI, the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) Global Digital Elevation Model (GDEM), was obtained from the Ministry of Economy, Trade, and Industry (METI) of Japan and NASA through Earthdata Search [26]. Required components (catchment area, flow width, and slope) calculating the TWI were computed using recommended algorithms [25]. Initially, terrain sinks were removed from the Digital Elevation Model data using Sink Removal with Fill Sinks in SAGA GIS 2.3.2. Then, the total catchment area was calculated by the Multiple Flow Direction algorithm [27] with the convergence at 1.0. Secondly, the slope was calculated in the algorithm of 10 parameters and third-order polynomial [28]. TWI was lastly computed using the pseudo-specific catchment area.

2.3.3. Elevation Data

Elevation data were obtained from Advanced Land Observing Satellite (ALOS) World 3D–30 m provided by the Japan Aerospace Exploration Agency (JAXA) [29].

2.3.4. Soil Properties

Data on soil properties were obtained from the International Soil Reference and Information Centre (ISRIC) [30]: Soil pH in H2O at a depth of 0 cm [31], soil texture at a depth of 0 cm [31], soil organic carbon content [31], the extractable aluminum content of the soil at a depth of 0–30 cm [32], and the extractable iron content of the soil at a depth of 0–30 cm [32]. Layer data closest to a depth of 5 cm, where sand flea larvae exist, were chosen [3]. Tunga spp. are found in sandy, clay, or laterite soils [18,33]. Many descriptions can be found for soil texture, such as sandy and clay soils, but there is little for laterite [2,9,11,12,34,35]. Therefore, the latter two mineral contents characterizing laterite were included in the exploratory analysis.

2.3.5. Distance from Households to the National Reserve or the Elephant Sanctuary

As wild animals in the Shimba Hills National Reserve and the Mwaluganje Elephant Sanctuary might have played a potential role in transmitting tungiasis around the buffer zone, Euclidean distance from households to the nearest reserve was calculated in QGIS 3.14.

2.4. Statistical Analysis

Generalized additive models (GAMs) [36] were used to create a predictive map of odds ratios (ORs) of household tungiasis infection over the study area. They were also used to estimate ORs for other risk factors. We used the MapGAM [37] and the gam [38] packages in R version 4.0.3 [39]. GAMs are extensions of typical linear regressions models (e.g., logistic regression model for binary data) that include smoothed terms that allow for the relaxing of typical assumptions of linearity between predictors and outcomes [40]. In the field of spatial epidemiology, mixed-mode GAMs, in which only location coordinates are non-parametrically smoothed, whereas other covariates keep a form of linear parametric prediction [36], have been used to investigate the geographic variation of diseases [36,41,42,43,44,45,46,47]. The below model was used to calculate ORs in this study:
log ( p i 1 p i ) = S ( x i ,   y i ) + Z i β + α
where pi denotes the probability that the ith household has one or more tungiasis cases; xi and yi are projected coordinates for the ith household; S(xi, yi) is a 2-dimensional smoothing of location; Zi is a vector of covariates for the ith household (household-level variables); β is a vector of regression parameters; α is an intercept [36].
A locally weighted scatterplot smoother (loess) was used for two-dimensional smoothing. A span size, specifying a percentage of neighborhood data points used for smoothing, was determined by minimizing the Akaike information criterion [36]. Regular grid points used for predictions were created using a shapefile of the HDSS area in the same projection as the household’s coordinates. The numbers of rows and columns of the grid were adjusted to make square grids rather than rectangle ones by reflecting the east-west and north-south distances. A function, modgam in the MapGAM package, fit a GAM model and estimated the log-odds at each of the grid points. The OR was calculated by the odds at each grid point and the median log odds over the study area [48]. To test the global null hypothesis that geolocation of household is unassociated with tungiasis occurrence, we conducted a permutation test with 1000 conditional permutations for a global p-value. As the conditional permutation test with a fixed span size induces an inflated type I error rate, the alpha risk was set to be 0.025 so that the type I error rate was within 0.05 [49]. At the same time, a pointwise permutation test was conducted to determine areas with significantly increased or decreased odds (hot spots and cold spots) [36]. The permutation distribution was calculated at each grid point using the same permutation data set for the global test (p = 0.05). Significantly increased or decreased OR area was defined as an area with grid points that rank in the higher 97.5% or the lower 2.5% of the pointwise permutation distribution, and contour lines of 0.975 and 0.025 indicate increased or decreased OR area, respectively [36]. The distribution pattern of the disease was recognized by significant hot and cold spots.
Spatial risk factors were examined by comparing a risk surface of a map between the unadjusted and an adjusted model. First, each variable was added to the unadjusted model to determine whether it was a spatial risk or not, as previously mentioned [42,46]. Then, variables that changed the surfaces by 10% were added to the unadjusted model together [42,46]. If the odds ratio increased or decreased by 10% from the unadjusted model at any grid point in the adjusted model, the adjusted variable was considered a spatial risk. As a result, all variables showed a risk surface change of more than 10%, and they were included in models as follows. In an adjusted model (model 1), we included the following: the number of children per household, the number of elderly per household, the number of males per household, type of floor (earthen/non-earthen), type of wall (earthen/non-earthen), type of roof (organic/non-organic), type of toilet (flush or pit latrine/none), NDVI, land cover (grasslands/croplands/savannas), TWI, elevation, soil pH, soil texture (sandy clay loam/sandy loam), soil organic carbon content, aluminum content in the soil, and iron content in the soil. In another adjusted model (model 2), we added the distance to the nearest animal reserve area to explore the effect of the distance on tungiasis risk. Models were assessed for multicollinearity by the variance inflation factor using a reduced model (i.e., without smoothing of location).

3. Results

3.1. Demographic Characteristics

Among 41,618 participants where tungiasis status was known, 483 individuals were excluded due to missing information, and 41,135 individuals comprising 7925 households were retained and used for analysis. A total of 461 people (1.1%) from 272 households (3.4%) had tungiasis at the time of the data collection. See Figure 1 for the scatter of the tungiasis positive and negative households. The age of participants ranged from one to 113 years (mean = 22 years), and the gender was almost equal: 50.7% female and 49.3% male. The data collection extended over both the dry and rainy seasons: 12.0% in March, 35.2% in April, 27.4% in May, 17.2% in June, 2.4% in July, 2.7% in August, 2.8% in September, 0.2% in November, and 0.1% in December. More than half of the data were collected in the rainy season in April, May, November, and December.

3.2. Spatial Analysis

In the unadjusted analysis, the variation of ORs displayed a pattern that corresponded with the geographical distribution of households with positive cases (Figure 2a). The optimal span size that minimized the AIC was 0.1, meaning that 10% of neighboring data were used for smoothing. Most areas had significantly increased or decreased ORs. A p-value for the global permutation test was less than 0.001 (<0.025: α), indicating that the occurrence of tungiasis was significantly associated with the geolocation of the household. Significantly increased ORs were observed among households in the east and west sides of the Mwaluganje Elephant Sanctuary and the southern part centering on the region convex southeastward (Figure 2a). In contrast, households in the northwest and northeast showed significantly decreased ORs.
Figure 2b shows the adjusted map (model 1) for the spatial risks. While the southern part still showed significantly increased ORs, there was a shrinking of areas at increased ORs on the east side of the Mwaluganje Elephant Sanctuary (See Figure 2b). The heterogeneous distribution of the selected variables (e.g., spatial aggregation of risks) partially explained the geographic distribution of hot spots of tungiasis. However, even after adjustment, a significant association between location and tungiasis remained (global p-value < 0.001), and the adjusted map did not become fully flat (Figure 2b). The optimal span size was 0.1 in model 1.
Figure 2c shows model 2, which was additionally adjusted for the distance to the nearest animal reserve area to see the additional impact. The optimal span size was 0.15 in model 2. The centering of the hot spot moved from the convex (Figure 2b) to the west nook (Figure 2c), reflecting the adjustment for the distance to the nearest animal reserve area. However, the southern region still had increased ORs, and the unexplained geographical risk variance remained (global p-value < 0.001). The variance inflation factors of the models were less than ten, indicating that multicollinearity might not have existed. No variables were excluded from the analysis.
Table 1 shows the used factors in the multivariate GAMs. In model 1, the number of children per household was associated with tungiasis (OR: 1.4, 95%CI: 1.3–1.5) at a host-related level. Earthen floors (OR: 3.2, 95%CI: 1.4–7.7) and organic roofs (OR: 1.7, 95%CI: 1.1–2.8) significantly affected the tungiasis status in socioeconomic factors. A positive association was observed among ecological factors at the household point in elevation (OR: 1.2, 95%CI: 1.1–1.3), aluminum content (OR: 1.1, 95%CI: 1.03–1.1), and iron content (OR: 1.2, 95%CI: 1.01–1.5). Although insignificant, indexes related to greenness and wetness, NDVI, and TWI showed positive associations with tungiasis. In model 2, the distance to the nearest animal reserve area was also significantly associated with tungiasis (OR: 0.6, 95%CI: 0.5–0.7). One kilometer away from the nearest animal reserve area, the OR was reduced to almost half. On the other hand, the iron content was not a significant spatial risk through an additional adjustment for the distance to the nearest animal reserve.

4. Discussion

This study revealed both the variability in the spatial distribution of tungiasis cases in a tungiasis-endemic area and the risk factors for tungiasis by adjusting for the spatial distribution of the potential risk factors. The pointwise permutation test found that the distribution of the disease was heterogeneous, with some areas having significantly increased or decreased ORs. The significance remained even after adjusting for other covariates: the number of children per household, the number of elderly per household, the number of males per household, type of floor, type of wall, type of roof, type of toilet, NDVI, land cover, TWI, elevation, soil pH, soil texture, soil organic carbon content, aluminum content in the soil, iron content in the soil, and distance to the nearest animal reserve area. Risk factors identified in this study were the number of children in the home, earthen floor, organic roof, elevation, aluminum content in the soil, and distance to the nearest animal reserve.
Our results support other research and have many implications for the complex network of factors that determine the household risk for tungiasis. First, an increased number of children in the home, earthen floor, and organic roof were associated with an increased spatial risk for household tungiasis risk, which agrees with the results of numerous studies on tungiasis [7,8,9,10,11,12,13,14,15,16]. Children spend much of their time on bare feet, especially in poor rural areas, which increases the chance of infection and transmission of tungiasis. Thus, it is more likely for tungiasis to be maintained in a household as a result of having more susceptible hosts. Organic material from grasses and leaves falling onto floors from thatched roofs might provide suitable sustenance for sand fleas living in the soil.
Next, elevation was shown to be positively associated with tungiasis. Tunga spp. tolerates high altitudes, even surviving and thriving at altitudes of over several thousand meters above sea level [50,51]. We found only one study on the relationship between elevation and tungiasis infection. This study reported no apparent spatial deviation of cases over elevation by visual comparison [12]. Taken together, there is a possibility that elevation was not a spatial driver, and other factors with a topographic distribution similar to elevation would have played a role instead in this study. Such factors would include water access, poverty, soil types, vegetation, and wildlife distributed according to elevation, some of which may be ecologically associated. As for elevation, more research is needed to account for the association of tungiasis with elevation.
Laterite, reddish-brown soil containing rich iron and aluminum oxide, can be found widely in tropical Africa [52,53], and its distribution in the African continent seems to overlap with the environmental suitability of tungiasis shown by Deka [6] (See Figure S1). A higher aluminum content in the soil was associated with a higher risk of the disease in this study. On the other hand, another characteristic of laterite, iron content in the soil, was not a significant spatial risk in model 2, which was additionally controlled for the distance to the nearest animal reserve area. This suggested that the iron content had no association with the distribution of tungiasis in the study area. One possible reason for this is the diverse composition of iron and aluminum in laterite [54]. Another possible reason is that we used the result of the model rather than actual measurements for the iron content [32].
As for aluminum, it is generally toxic to living things. Exposure to aluminum in food or medium leads to behavior abnormality, decreased fertility, growth inhibition, and a shortened lifespan in insects [55,56,57,58,59]. Considering that sand fleas’ off-host stage development occurs in soil and that larvae feed on organic matter in the environment [3,18,60], it is not impossible that a higher soil aluminum content negatively affects sand fleas’ propagation. In addition, aluminum content data might not reflect actual measurement values because they are based on a prediction by a model [32]. On the other hand, aluminum ions in the soil are strongly toxic to the plant and cause a low crop yield. A lower soil pH than 5.5, in which aluminum ions can be soluble and toxic [61,62], was found in the study area (pH 5.3–6.9), while the aluminum content in the study area (509 to 845 mg/kg) took around the median of that over the African continent (203 to 1784 mg/kg). It is thus not impossible that the environment and vegetation according to such soil may have affected sand fleas. As for the biological effect of aluminum on sand fleas, an aluminum-controlled artificial culture would be necessary for precise insight, although it is not available for sand fleas yet. Given the scarce evidence available, we suggest future studies on the relationships between soil metal content and tungiasis. Future studies might test hypotheses of the aluminum soil content and tungiasis risk on a larger scale using actual measurements of the aluminum soil content.
We found that households located closer to animal reserves had increased odds of tungiasis cases in the home. Since dogs, cats, pigs, and other livestock have been suggested to serve as animal reservoirs of sand fleas [9,11,12,13,14,17,63,64,65], wildlife in the reserve areas possibly contributed to the spreading of sand flea eggs and adults around households. Previous research in the same study area has suggested that small mammals and dogs, presumably moving in and out of the wildlife park, might bring parasites in proximity to the home, furthering the household transmission risk [17]. Wildlife themselves, or intermediary dogs and small animals that roam in and out of the bounds of the wildlife reserve, might explain why the infection risk is high in households near the park. Anecdotal reports also suggest that humans encroach on the park to forage for food and medicines and to graze animals, possibly raising the infection risk. It is unknown, however, exactly how animal habitats and human residences overlap and how likely it is that the development and infection of sand fleas occur outside living areas, where the majority of transmission is believed to occur [3,18,66]. Further studies are needed to address the dynamic ecology of sand fleas and their association with human, livestock, and wildlife interactions.
The seasonality of tungiasis [67] and the importance of moisture in the development of sand flea eggs [3] suggest the existence of an optimal range of moisture for sand fleas. However, soil moisture-related variables (NDVI and TWI) were not associated spatially with tungiasis in the study. A possible reason is that the soil moisture-related variables would have been within the tolerable ranges of sand fleas’ development, considering that the Kwale HDSS area is not so large as to have climatic differences. Another reason would be that the median value of NDVI in 2011 might have summarized the variances meaninglessly over dry and rainy seasons, although there was no choice but to use the median due to data collection over ten months in 2011. Furthermore, artificial differences in the sand flea’s habitat (i.e., ground water or drainage structure of a house) would have been more significant than natural differences.
As per our hypotheses, the combination of environmental and household factors accounted for the spatial pattern of tungiasis; however, the unexplained variance of risk remained even after adjustment. It was suggested that unmeasured spatial factors would have contributed to the spatial disparity of the disease. As noted, as a limitation, data of some known risk factors, such as socioeconomic status and animal existence at the household, were not included in the analysis. Thus, these factors might explain the residual differences in the risk of tungiasis.
There were several limitations to this study. First, as the data used in this study were collected in 2011, our results may not reflect current ecological and public health conditions. However, the results may reflect the spatial variation of the tungiasis distribution and spatial risk factors. These results can inform future data collection and research efforts to disentangle the complex ecology of tungiasis. Second, data collection on the tungiasis status took ten months across the dry and rainy seasons in 2011. Since more than half of the data were collected in rainy seasons, when transmission of sand fleas tends to be low [67], the prevalence and distribution of tungiasis would not have reflected the most infected state in the high season. Although the potentially missed cases may have affected the spatial analysis, this study identified spatial risk factors successfully. At the same time, the extended window of data collection would explain why the prevalence rate of 1.1% at the individual level in the study was remarkably lower than other studies conducted in the dry season in Kenya, which showed a 20% or more prevalence rate [7,12,16,68]. Third, the data for the socioeconomic status and livestock possession were not included in the analysis because of a five-year gap in data collection. There is a possibility that these variables, known risk factors of tungiasis, could have explained the residual geographic variation if included in the analysis. Moreover, a possibility of information bias was noticed as a potential cause of the heterogeneous distribution. In the HDSS area, of nine data collectors assigned to each area for data collection, only one covered the convex area where the center of the positive cases was located. Thus, the bias of the data collector could not be ruled out, although the distribution of tungiasis was consistent with the local awareness. A significant challenge with this study is within the method of data collection for infection status. Infection was based on self-reports or from third-party reports from other household members. In addition, the stigmatization of tungiasis in the community might have enhanced underreporting [69,70,71]. While we are unlikely to have seen false positives, there is a high possibility that some cases were missed, inflating the number of negatives. While we recognize this as a problem, we believe that the measure of household infection is sufficient and might provide a good proxy of community transmission. Despite these limitations, the present study highlights the need for further spatial studies to reveal underlying factors on tungiasis distribution.

5. Conclusions

In conclusion, the study showed the spatial analysis for heterogeneously distributed tungiasis using the geolocation of the household, household-based variables, and point-based ecological factors. As a result, the number of children per household, earthen floor, organic roof, elevation, aluminum content in the soil, and distance to the nearest animal reserve area were spatial risk factors at the household level in this study. Three established risk factors (children, earthen floor, and organic roof) could be promising predictors of the geographical distribution of tungiasis, whereas the roles of aluminum, elevation, and wildlife need further investigation. Considering that the variables did not fully explain the geographical disparity of tungiasis, a more complementary inclusion and consideration of variables are recommended for future work to elucidate spatial risk factors. Such spatial work will contribute to identifying high-risk populations and areas, leading to adequate control approaches.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/tropicalmed7010002/s1, Table S1: Ecological data used in the analysis and its sources and details, Figure S1: Distribution of environmental suitability for tungiasis, laterite, aluminum content in soil, and iron content in soil in the African continent. R code used for analysis in this study is available online at https://github.com/hayalpaca/Spatial_risk_of_tungiasis (accessed on 10 November 2021).

Author Contributions

Conceptualization, A.H., P.S.L. and S.K.; data curation, M.N., S.W.M. and M.C.; formal analysis, A.H.; investigation, M.N., S.W.M., M.C. and K.G.; methodology, A.H.; project administration, M.K. and S.K.; resources, M.K.; software, A.H.; supervision, P.S.L., M.K., K.G. and S.K.; validation, P.S.L. and S.K.; visualization, A.H.; writing—original draft, A.H.; writing—review and editing, P.S.L. and S.K. All authors have read and agreed to the published version of the manuscript.

Funding

The study was supported by the Basic Expense for the Kenya Research Station (Management Expenses Grants from the Government of Japan).

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the Scientific and Ethics Review Unit of Kenya Medical Research Institute (SSC number 1088) and the Institutional Review Board of the Institute of Tropical Medicine, Nagasaki University (IRB # 06060604).

Informed Consent Statement

At the first visit of a household, a data collector obtained written consent from the head of the household on behalf of his/her family.

Data Availability Statement

The datasets generated and/or analyzed during the current study are not publicly available to secure the anonymity of the geolocation of households who participated in this study, but are available from the corresponding author on reasonable request.

Acknowledgments

We are grateful to all the participants in this study. We would also like to thank the Kenyan Ministry of Health, the Director of KEMRI, NUITM-KEMRI, the Kwale County Government for their support, and NASA LP DAAC, METI/NASA, JAXA, and ISRIC for their data. Chemi Chemi community-based organization contributed to data collection in the fieldwork. We humbly acknowledge the help of and the Umoja Women’s Group (Kwale). We also thank the Institute for Global Change Biology at the University of Michigan. A.H. is a Research Fellow of Japan Society for the Promotion of Science and grateful for academic and financial support to the Program for Nurturing Global Leaders in Tropical and Emerging Communicable Diseases, Graduate School of Biomedical Sciences, Nagasaki University and Nagasaki University “Doctoral Program for World-leading Innovative and Smart Education” for Global Health, “Global Health Elite Programme for Building a Healthier World.” Finally, we would like to express our gratitude to research assistants in the HDSS, Kazuchiyo Miyamichi and Tarek Numair in the Department of Eco-epidemiology, NUITM for their great support, and Brittni Delmaine for English language editing and reviewing.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Feldmeier, H.; Heukelbach, J.; Ugbomoiko, U.S.; Sentongo, E.; Mbabazi, P.; von Samson-Himmelstjerna, G.; Krantz, I.; International Expert Group for Tungiasis. Tungiasis—A neglected disease with many challenges for global public health. PLoS Negl. Trop. Dis. 2014, 8, e3133. [Google Scholar] [CrossRef] [Green Version]
  2. Heukelbach, J.; De Oliveira, F.A.S.; Hesse, G.; Feldmeier, H. Tungiasis: A neglected health problem of poor communities. Trop. Med. Int. Health 2001, 6, 267–272. [Google Scholar] [CrossRef]
  3. Nagy, N.; Abari, E.; D’Haese, J.; Calheiros, C.; Heukelbach, J.; Mencke, N.; Feldmeier, H.; Mehlhorn, H. Investigations on the life cycle and morphology of Tunga penetrans in Brazil. Parasitol. Res. 2007, 101 (Suppl. S2), 233–242. [Google Scholar] [CrossRef]
  4. Eisele, M.; Heukelbach, J.; Van Marck, E.; Mehlhorn, H.; Meckes, O.; Franck, S.; Feldmeier, H. Investigations on the biology, epidemiology, pathology and control of Tunga penetrans in Brazil: I. Natural history of tungiasis in man. Parasitol. Res. 2003, 90, 87–99. [Google Scholar] [CrossRef] [PubMed]
  5. World Health Organization. Tungiasis. Available online: https://www.who.int/news-room/fact-sheets/detail/tungiasis (accessed on 1 August 2021).
  6. Deka, M.A. Mapping the Geographic Distribution of Tungiasis in Sub-Saharan Africa. Trop. Med. Infect. Dis. 2020, 5, 122. [Google Scholar] [CrossRef]
  7. Wiese, S.; Elson, L.; Reichert, F.; Mambo, B.; Feldmeier, H. Prevalence, intensity and risk factors of tungiasis in Kilifi County, Kenya: I. Results from a community-based study. PLoS Negl. Trop. Dis. 2017, 11, e0005925. [Google Scholar] [CrossRef] [Green Version]
  8. Muehlen, M.; Heukelbach, J.; Wilcke, T.; Winter, B.; Mehlhorn, H.; Feldmeier, H. Investigations on the biology, epidemiology, pathology and control of Tunga penetrans in Brazil. II. Prevalence, parasite load and topographic distribution of lesions in the population of a traditional fishing village. Parasitol. Res. 2003, 90, 449–455. [Google Scholar] [CrossRef]
  9. Muehlen, M.; Feldmeier, H.; Wilcke, T.; Winter, B.; Heukelbach, J. Identifying risk factors for tungiasis and heavy infestation in a resource-poor community in northeast Brazil. Trans. R. Soc. Trop. Med. Hyg. 2006, 100, 371–380. [Google Scholar] [CrossRef]
  10. Wilcke, T.; Heukelbach, J.; Moura, R.C.S.; Kerr-Pontes, L.R.S.; Feldmeier, H. High prevalence of tungiasis in a poor neighbourhood in Fortaleza, Northeast Brazil. Acta Trop. 2002, 83, 255–258. [Google Scholar] [CrossRef]
  11. Ugbomoiko, U.S.; Ariza, L.; Ofoezie, I.E.; Heukelbach, J. Risk factors for tungiasis in Nigeria: Identification of targets for effective intervention. PLoS Negl. Trop. Dis. 2007, 1, e87. [Google Scholar] [CrossRef] [Green Version]
  12. Nyangacha, R.M.; Odongo, D.; Oyieke, F.; Bii, C.; Muniu, E.; Chasia, S.; Ochwoto, M. Spatial distribution, prevalence and potential risk factors of Tungiasis in Vihiga County, Kenya. PLoS Negl. Trop. Dis. 2019, 13, e0007244. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Wafula, S.T.; Ssemugabo, C.; Namuhani, N.; Musoke, D.; Ssempebwa, J.; Halage, A.A. Prevalence and risk factors associated with tungiasis in Mayuge district, Eastern Uganda. Pan Afr. Med. J. 2016, 24, 77. [Google Scholar] [CrossRef]
  14. Girma, M.; Astatkie, A.; Asnake, S. Prevalence and risk factors of tungiasis among children of Wensho district, southern Ethiopia. BMC Infect. Dis. 2018, 18, 456. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Nsanzimana, J.; Karanja, S.; Kayongo, M.; Nyirimanzi, N.; Umuhoza, H.; Murangwa, A.; Muganga, R.; Musafili, A. Factors associated with tungiasis among primary school children: A cross-sectional study in a rural district in Rwanda. BMC Public Health 2019, 19, 1192. [Google Scholar] [CrossRef]
  16. Elson, L.; Wiese, S.; Feldmeier, H.; Fillinger, U. Prevalence, intensity and risk factors of tungiasis in Kilifi County, Kenya II: Results from a school-based observational study. PLoS Negl. Trop. Dis. 2019, 13, e0007326. [Google Scholar] [CrossRef]
  17. Larson, P.S.; Ono, M.; Changoma, M.; Goto, K.; Kaneko, S.; Moji, K.; Minakawa, N. Presence of dogs and proximity to a wildlife reserve increase household level risk of tungiasis in Kwale, Kenya. Trop. Med. Health 2021, 49, 54. [Google Scholar] [CrossRef]
  18. Linardi, P.M.; Calheiros, C.M.; Campelo-Junior, E.B.; Duarte, E.M.; Heukelbach, J.; Feldmeier, H. Occurrence of the off-host life stages of Tunga penetrans (Siphonaptera) in various environments in Brazil. Ann. Trop. Med. Parasitol. 2010, 104, 337–345. [Google Scholar] [CrossRef] [PubMed]
  19. Kaneko, S.; K’Opiyo, J.; Kiche, I.; Wanyua, S.; Goto, K.; Tanaka, J.; Changoma, M.; Ndemwa, M.; Komazawa, O.; Karama, M.; et al. Health and Demographic Surveillance System in the Western and coastal areas of Kenya: An infrastructure for epidemiologic studies in Africa. J. Epidemiol. 2012, 22, 276–285. [Google Scholar] [CrossRef] [Green Version]
  20. Pampiglione, S.; Trentini, M.; Fioravanti, M.L.; Onore, G.; Rivasi, F. Additional description of a new species of Tunga (Siphonaptera) from Ecuador. Parasite 2003, 10, 9–15. [Google Scholar] [CrossRef] [Green Version]
  21. Ngunjiri, J.; Keiyoro, P. Soil Factors Influencing the Occurrence of Tunga penetrans in Kenya: Investigation of Soil Factors Influencing the Prevalence of the Tunga penetrans (Jigger Flea) in Muranga County—Kenya; LAMBERT Academic Publishing: Chisinau, Moldova, 2011. [Google Scholar]
  22. QGIS Organization. QGIS Geographic Information System; QGIS Association: Zurich, Switzerland, 2021. [Google Scholar]
  23. Earth Engine Data Catalog. Available online: https://developers.google.com/earth-engine/datasets/ (accessed on 17 February 2021).
  24. Gruber, S.; Peckham, S. Chapter 7 Land-Surface Parameters and Objects in Hydrology. In Geomorphometry—Concepts, Software, Applications; Developments in Soil Science; Hengl, T., Reuter, H.I., Eds.; Elsevier: Amsterdam, The Netherlands, 2009; Volume 33, pp. 171–194. [Google Scholar]
  25. Kopecky, M.; Macek, M.; Wild, J. Topographic Wetness Index calculation guidelines based on measured soil moisture and plant species composition. Sci. Total Environ. 2021, 757, 143785. [Google Scholar] [CrossRef]
  26. Earthdata Search. Available online: https://search.earthdata.nasa.gov/search (accessed on 3 March 2021).
  27. Freeman, T.G. Calculating catchment area with divergent flow based on a regular grid. Comput. Geosci. 1991, 17, 413–422. [Google Scholar] [CrossRef]
  28. Haralick, R.M. Ridges and valleys on digital images. Comput. Vis. Graph. Image Process. 1983, 22, 28–38. [Google Scholar] [CrossRef]
  29. ALOS Global Digital Surface Model. Available online: https://www.eorc.jaxa.jp/ALOS/en/aw3d30/index.htm (accessed on 18 March 2021).
  30. ISRIC Data Hub. Available online: https://data.isric.org/geonetwork/srv/eng/catalog.search#/home (accessed on 20 March 2021).
  31. Hengl, T.; Mendes de Jesus, J.; Heuvelink, G.B.; Ruiperez Gonzalez, M.; Kilibarda, M.; Blagotic, A.; Shangguan, W.; Wright, M.N.; Geng, X.; Bauer-Marschallinger, B.; et al. SoilGrids250m: Global gridded soil information based on machine learning. PLoS ONE 2017, 12, e0169748. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Hengl, T.; Leenaars, J.G.B.; Shepherd, K.D.; Walsh, M.G.; Heuvelink, G.B.M.; Mamo, T.; Tilahun, H.; Berkhout, E.; Cooper, M.; Fegraus, E.; et al. Soil nutrient maps of Sub-Saharan Africa: Assessment of soil nutrient content at 250 m spatial resolution using machine learning. Nutr. Cycl. Agroecosyst. 2017, 109, 77–102. [Google Scholar] [CrossRef] [Green Version]
  33. Heukelbach, J. Tungiasis. Available online: http://www.orpha.net/data/patho/GB/uk-Tungiasis.pdf (accessed on 26 October 2021).
  34. Winter, B.; Oliveira, F.A.; Wilcke, T.; Heukelbach, J.; Feldmeier, H. Tungiasis-related knowledge and treatment practices in two endemic communities in northeast Brazil. J. Infect. Dev. Ctries. 2009, 3, 458–466. [Google Scholar] [CrossRef] [Green Version]
  35. Chadee, D.D. Tungiasis among five communities in south-western Trinidad, West Indies. Ann. Trop. Med. Parasitol. 1998, 92, 107–113. [Google Scholar] [CrossRef]
  36. Webster, T.; Vieira, V.; Weinberg, J.; Aschengrau, A. Method for mapping population-based case-control studies: An application using generalized additive models. Int. J. Health Geogr. 2006, 5, 26. [Google Scholar] [CrossRef] [Green Version]
  37. Bai, L.; Bartell, S.; Bliss, R.; Vieira, V. MapGAM: Mapping Smoothed Effect Estimates from Individual-Level Data, R package version 1.2-5.; R Foundation for Statistical Computing: Vienna, Austria, 2019. [Google Scholar]
  38. Hastie, T. gam: Generalized Additive Models, R package version 1.20; R Foundation for Statistical Computing: Vienna, Austria, 2020. [Google Scholar]
  39. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2020. [Google Scholar]
  40. Hastie, T.; Tibshirani, R. Generalized Additive Models. Statist. Sci. 1986, 1, 297–310. [Google Scholar] [CrossRef]
  41. Hoffman, K.; Aschengrau, A.; Webster, T.F.; Bartell, S.M.; Vieira, V.M. Associations between residence at birth and mental health disorders: A spatial analysis of retrospective cohort data. BMC Public Health 2015, 15, 688. [Google Scholar] [CrossRef] [Green Version]
  42. Hoffman, K.; Weisskopf, M.G.; Roberts, A.L.; Raz, R.; Hart, J.E.; Lyall, K.; Hoffman, E.M.; Laden, F.; Vieira, V.M. Geographic Patterns of Autism Spectrum Disorder among Children of Participants in Nurses’ Health Study II. Am. J. Epidemiol. 2017, 186, 834–842. [Google Scholar] [CrossRef]
  43. Bristow, R.E.; Chang, J.; Ziogas, A.; Anton-Culver, H.; Vieira, V.M. Spatial analysis of adherence to treatment guidelines for advanced-stage ovarian cancer and the impact of race and socioeconomic status. Gynecol. Oncol. 2014, 134, 60–67. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Ribeiro, A.I.; Santos, A.C.; Vieira, V.M.; Barros, H. Hotspots of childhood obesity in a large metropolitan area: Does neighbourhood social and built environment play a part? Int. J. Epidemiol. 2020, 49, 934–943. [Google Scholar] [CrossRef]
  45. Vieira, V.M.; VoPham, T.; Bertrand, K.A.; James, P.; DuPre, N.; Tamimi, R.M.; Laden, F.; Hart, J.E. Contribution of socioeconomic and environmental factors to geographic disparities in breast cancer risk in the Nurses’ Health Study II. Environ. Epidemiol. 2020, 4, e080. [Google Scholar] [CrossRef]
  46. Vieira, V.M.; Fabian, M.P.; Webster, T.F.; Levy, J.I.; Korrick, S.A. Spatial Variability in ADHD-Related Behaviors among Children Born to Mothers Residing Near the New Bedford Harbor Superfund Site. Am. J. Epidemiol. 2017, 185, 924–932. [Google Scholar] [CrossRef] [Green Version]
  47. Vieira, V.; Webster, T.; Weinberg, J.; Aschengrau, A.; Ozonoff, D. Spatial analysis of lung, colorectal, and breast cancer on Cape Cod: An application of generalized additive models to case-control data. Environ. Health 2005, 4, 11. [Google Scholar] [CrossRef] [Green Version]
  48. Bai, L.; Gillen, D.L.; Bartell, S.M.; Vieira, V.M. Mapping Smoothed Spatial Effect Estimates from Individual-Level Data: MapGAM. R J. 2020, 12, 32. [Google Scholar] [CrossRef]
  49. Young, R.L.; Weinberg, J.; Vieira, V.; Ozonoff, A.; Webster, T.F. Generalized Additive Models and Inflated Type I Error Rates of Smoother Significance Tests. Comput. Stat. Data Anal. 2011, 55, 366–374. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Gadisa, E.; Jote, K. Prevalence and factors associated with intestinal parasitic infection among under-five children in and around Haro Dumal Town, Bale Zone, Ethiopia. BMC Pediatrics 2019, 19, 385. [Google Scholar] [CrossRef]
  51. Pampiglione, S.; Fioravanti, M.L.; Gustinelli, A.; Onore, G.; Mantovani, B.; Luchetti, A.; Trentini, M. Sand flea (Tunga spp.) infections in humans and domestic animals: State of the art. Med. Vet. Entomol. 2009, 23, 172–186. [Google Scholar] [CrossRef]
  52. Areola, O. Soils. In The Physical Geography of Africa; Adams, W., Goudie, A.S., Orme, A.R., Eds.; Oxford University Press: New York, NY, USA, 1996; pp. 134–147. [Google Scholar]
  53. Fookes, P.G. Tropical Residual Soils, a Geological Society Engineering Group Working Party Revised Report; The Geological Society: London, UK, 1997. [Google Scholar]
  54. Ghosh, S.; Guchhait, S.K. Introduction to Laterite Study. In Laterites of the Bengal Basin: Characterization, Geochronology and Evolution; Ghosh, S., Guchhait, S.K., Eds.; Springer International Publishing: Cham, Switzerland, 2020; pp. 1–16. [Google Scholar]
  55. Bayliak, M.M.; Lylyk, M.P.; Gospodaryov, D.V.; Kotsyubynsky, V.O.; Butenko, N.V.; Storey, K.B.; Lushchak, V.I. Protective effects of alpha-ketoglutarate against aluminum toxicity in Drosophila melanogaster. Comp. Biochem. Physiol. C Toxicol. Pharmacol. 2019, 217, 41–53. [Google Scholar] [CrossRef] [PubMed]
  56. Wu, Z.; Du, Y.; Xue, H.; Wu, Y.; Zhou, B. Aluminum induces neurodegeneration and its toxicity arises from increased iron accumulation and reactive oxygen species (ROS) production. Neurobiol. Aging 2012, 33, 199.e1–199.e12. [Google Scholar] [CrossRef] [PubMed]
  57. Kijak, E.; Rosato, E.; Knapczyk, K.; Pyza, E. Drosophila melanogaster as a model system of aluminum toxicity and aging. Insect Sci. 2014, 21, 189–202. [Google Scholar] [CrossRef]
  58. Anand, A.S.; Gahlot, U.; Prasad, D.N.; Kohli, E. Aluminum oxide nanoparticles mediated toxicity, loss of appendages in progeny of Drosophila melanogaster on chronic exposure. Nanotoxicology 2019, 13, 977–989. [Google Scholar] [CrossRef] [PubMed]
  59. Liu, L.; Qian, X.; Chao, M.; Zhao, Y.; Huang, J.; Wang, T.; Sun, F.; Ling, E.; Song, H. Aluminum toxicity related to SOD and expression of presenilin and CREB in Bombyx mori. Arch. Insect Biochem. Physiol. 2018, 99, e21480. [Google Scholar] [CrossRef] [PubMed]
  60. Harvey, T.V.; Linardi, P.M.; Carlos, R.S.A.; Heukelbach, J. Tungiasis in domestic, wild, and synanthropic animals in Brazil. Acta Trop. 2021, 222, 106068. [Google Scholar] [CrossRef]
  61. Krstic, D.; Djalovj, I.; Nikezic, D.; Bjelj, D. Aluminium in Acid Soils: Chemistry, Toxicity and Impact on Maize Plants. In Food Production Approaches, Challenges and Tasks; Intech: Rijeka, Croatia, 2012. [Google Scholar] [CrossRef] [Green Version]
  62. Yamamoto, Y. Aluminum toxicity in plant cells: Mechanisms of cell death and inhibition of cell elongation. Soil Sci. Plant Nutr. 2019, 65, 41–55. [Google Scholar] [CrossRef] [Green Version]
  63. Heukelbach, J.; Costa, A.M.L.; Wilcke, T.; Mencke, N.; Feldmeier, H. The animal reservoir of Tunga penetrans in severely affected communities of north-east Brazil. Med. Vet. Entomol. 2004, 18, 329–335. [Google Scholar] [CrossRef]
  64. Mutebi, F.; Krucken, J.; Feldmeier, H.; Waiswa, C.; Mencke, N.; Sentongo, E.; von Samson-Himmelstjerna, G. Animal Reservoirs of Zoonotic Tungiasis in Endemic Rural Villages of Uganda. PLoS Negl. Trop. Dis. 2015, 9, e0004126. [Google Scholar] [CrossRef] [Green Version]
  65. Mutebi, F.; Krucken, J.; Feldmeier, H.; Waiswa, C.; Mencke, N.; von Samson-Himmelstjerna, G. Tungiasis-associated morbidity in pigs and dogs in endemic villages of Uganda. Parasites Vectors 2016, 9, 44. [Google Scholar] [CrossRef] [Green Version]
  66. Calheiros, C.M.L. Aspectos Biológicos e Ecológicos de Tunga penetrans (L., 1758) (Siphonaptera Tungidae) em Áreas Endêmicas Brasileiras. Ph.D. Thesis, Universidade Federal de Minas Gerais, Belo Horizonte, Minas Gerais, Brazil, 2007. [Google Scholar]
  67. Heukelbach, J.; Wilcke, T.; Harms, G.; Feldmeier, H. Seasonal variation of tungiasis in an endemic community. Am. J. Trop. Med. Hyg. 2005, 72, 145–149. [Google Scholar] [CrossRef]
  68. Mwangi, J.N.; Ozwara, H.S.; Gicheru, M.M. Epidemiology of tunga penetrans infestation in selected areas in Kiharu constituency, Murang’a County, Kenya. Trop. Dis. Travel Med. Vaccines 2015, 1, 13. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  69. Feldmeier, H.; Sentongo, E.; Krantz, I. Tungiasis (sand flea disease): A parasitic disease with particular challenges for public health. Eur. J. Clin. Microbiol. Infect. Dis. 2013, 32, 19–26. [Google Scholar] [CrossRef]
  70. Kimani, B.; Nyagero, J.; Ikamari, L. Knowledge, attitude and practices on jigger infestation among household members aged 18 to 60 years: Case study of a rural location in Kenya. Pan Afr. Med. J. 2012, 13 (Suppl. S1), 7. [Google Scholar] [PubMed]
  71. Kimotho, S.; Miller, A.N.; Ngure, P. Managing communication surrounding tungiasis stigma in Kenya. Communicatio 2015, 41, 523–542. [Google Scholar] [CrossRef]
Figure 1. Distribution of the households with or without tungiasis in the Kwale-HDSS area.
Figure 1. Distribution of the households with or without tungiasis in the Kwale-HDSS area.
Tropicalmed 07 00002 g001
Figure 2. Geographical distribution of estimated odds ratios of tungiasis in the household. For easy comparison, a color band is fixed in an arbitrary range. Regions at the edge of the study area and no data or regions with odds ratios (ORs) outside the fixed range are shown in white. Based on the permutation test, contour lines indicate areas with significantly increased or decreased ORs. (a) Unadjusted; (b) Model 1 adjusted for the number of children per household, the number of elderly per household, the number of males per household, type of floor, type of wall, type of roof, type of toilet, Normalized Difference Vegetation Index (NDVI), land cover, Topographic Wetness Index (TWI), elevation, soil pH, soil texture, soil organic carbon content, aluminum content in the soil, and iron content in the soil; (c) Model 2 adjusted for the number of children per household, the number of elderly per household, the number of males per household, type of floor, type of wall, type of roof, type of toilet, NDVI, land cover, TWI, elevation, soil pH, soil texture, soil organic carbon content, aluminum content in the soil, iron content in the soil, and distance to the nearest animal reserve area.
Figure 2. Geographical distribution of estimated odds ratios of tungiasis in the household. For easy comparison, a color band is fixed in an arbitrary range. Regions at the edge of the study area and no data or regions with odds ratios (ORs) outside the fixed range are shown in white. Based on the permutation test, contour lines indicate areas with significantly increased or decreased ORs. (a) Unadjusted; (b) Model 1 adjusted for the number of children per household, the number of elderly per household, the number of males per household, type of floor, type of wall, type of roof, type of toilet, Normalized Difference Vegetation Index (NDVI), land cover, Topographic Wetness Index (TWI), elevation, soil pH, soil texture, soil organic carbon content, aluminum content in the soil, and iron content in the soil; (c) Model 2 adjusted for the number of children per household, the number of elderly per household, the number of males per household, type of floor, type of wall, type of roof, type of toilet, NDVI, land cover, TWI, elevation, soil pH, soil texture, soil organic carbon content, aluminum content in the soil, iron content in the soil, and distance to the nearest animal reserve area.
Tropicalmed 07 00002 g002
Table 1. Spatial risks used for adjustment and odds ratios in the multivariate analysis.
Table 1. Spatial risks used for adjustment and odds ratios in the multivariate analysis.
Household Based VariablesTungiasis n = 7925Mean (SD)Model 1Model 2
NegativePositiveOR95% CIp-ValueOR95% CIp-Value
n = 7653n = 272
Number of children per household 2.3 (1.9)1.41.3-1.5<0.011.41.3-1.5<0.01
Number of elderly per household 0.28 (0.6)1.10.8-1.4>0.051.10.8-1.4>0.05
Number of males per household 2.6 (1.8)1.10.97-1.2>0.051.10.97-1.2>0.05
FloorNon-earthen11957 1.0 1.0
Earthen6458265 3.21.4-7.7<0.013.11.3-7.6<0.05
WallNon-earthen8916 1.0 1.0
Earthen6762266 0.50.2-1.5>0.050.50.2-1.5>0.05
RoofNon-organic219327 1.0 1.0
Organic5460245 1.71.1-2.8<0.051.81.1-2.9<0.05
ToiletNone4235142 1.0 1.0
Flush or pit latrine3418130 0.90.7-1.2>0.050.90.7-1.2>0.05
Normalized Difference Vegetation Index (NDVI) 157 (7.1)1.010.97-1.04>0.051.0040.97-1.04>0.05
Land coverGrasslands6722234 1.0 1.0
Croplands73510 1.10.5-2.4>0.051.20.5-2.8>0.05
Savannas19628 0.70.4-1.4>0.050.80.4-1.5>0.05
Topographic Wetness Index (TWI) 6.0 (1.7)1.020.9-1.1>0.051.020.9-1.1>0.05
Elevation (10 m) 24 (8.0) *1.21.1-1.3<0.011.21.1-1.3<0.01
Soil pH 6.0 (0.3)0.90.3-2.2>0.050.90.3-2.2>0.05
Soil textureSandy clay loam5988267 1.0 1.0
Sandy loam16655 0.90.3-2.5>0.050.90.3-2.5>0.05
Soil organic carbon content (g/kg) 35 (30)0.9990.99-1.003>0.051.0010.995-1.01>0.05
Aluminum content of soil (10 mg/kg) 69 (7.2) *1.11.03-1.1<0.011.11.02-1.1<0.01
Iron content of soil (10 mg/kg) 15 (1.1) *1.21.01-1.5<0.051.20.97-1.5>0.05
Distance to the nearest animal reserve (km) 3.6 (3.4) 0.60.5-0.7<0.01
Model 1: a generalized additive model adjusted for the number of children per household, the number of elderly per household, the number of males per household, type of floor, type of wall, type of roof, type of toilet, NDVI, land cover, TWI, elevation, soil pH, soil texture, soil organic carbon content, aluminum content in the soil, and iron content in the soil; Model 2: a generalized additive model additionally adjusted for the distance to the nearest animal reserve area; *: Note that the units of these variables are ten times larger than usual. Actual mean values are 240 m for elevation, 690 mg/kg for aluminum of soil, and 150 mg/kg for iron content of soil, as well as standard deviations; CI: confidence interval; OR: odds ratio; SD: standard deviation.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Hyuga, A.; Larson, P.S.; Ndemwa, M.; Muuo, S.W.; Changoma, M.; Karama, M.; Goto, K.; Kaneko, S. Environmental and Household-Based Spatial Risks for Tungiasis in an Endemic Area of Coastal Kenya. Trop. Med. Infect. Dis. 2022, 7, 2. https://doi.org/10.3390/tropicalmed7010002

AMA Style

Hyuga A, Larson PS, Ndemwa M, Muuo SW, Changoma M, Karama M, Goto K, Kaneko S. Environmental and Household-Based Spatial Risks for Tungiasis in an Endemic Area of Coastal Kenya. Tropical Medicine and Infectious Disease. 2022; 7(1):2. https://doi.org/10.3390/tropicalmed7010002

Chicago/Turabian Style

Hyuga, Ayako, Peter S. Larson, Morris Ndemwa, Sheru W. Muuo, Mwatasa Changoma, Mohamed Karama, Kensuke Goto, and Satoshi Kaneko. 2022. "Environmental and Household-Based Spatial Risks for Tungiasis in an Endemic Area of Coastal Kenya" Tropical Medicine and Infectious Disease 7, no. 1: 2. https://doi.org/10.3390/tropicalmed7010002

Article Metrics

Back to TopTop