Neurocognitive and mental health outcomes in children with tungiasis: a cross-sectional study in rural Kenya and Uganda

Abstract

Background

Tungiasis, a neglected tropical parasitosis, disproportionately affects children. Few empirical studies have reported neurocognitive and mental health outcomes of children with ectoparasitic skin diseases like tungiasis. Pathophysiology of tungiasis suggests it could detrimentally affect cognition and behaviour. This study pioneered the investigation of neurocognitive and mental health outcomes in children with tungiasis.

Methods

This was a multi-site cross-sectional study including 454 quasi-randomly sampled school-children aged 8–14 from 48 randomly selected schools in two counties in Kenya and a district in Uganda. The participants were stratified into infected and uninfected based on the presence of tungiasis. The infected were further classified into mild and severe infection groups based on the intensity of the infection. Adapted, validated, and standardized measures of cognition and mental health such as Raven Matrices and Child Behaviour Checklist were used to collect data. Statistical tests including a multilevel, generalized mixed-effects linear models with family link set to identity were used to compare the scores of uninfected and infected children and to identify other potential risk factors for neurocognitive and behavioural outcomes.

Results

When adjusted for covariates, mild infection was associated with lower scores in literacy [adjusted β(aβ) = − 8.9; 95% confidence interval (CI) − 17.2, − 0.6], language (aβ = − 1.7; 95% CI − 3.2, − 0.3), cognitive flexibility (aβ = − 6.1; 95% CI − 10.4, − 1.7) and working memory (aβ = − 0.3; 95% CI − 0.6, − 0.1). Severe infection was associated with lower scores in literacy (aβ = − 11.0; 95% CI − 19.3, − 2.8), response inhibition, (aβ = − 2.2; 95% CI − 4.2, − 0.2), fine motor control (aβ = − 0.7; 95% CI − 1.1, − 0.4) and numeracy (aβ = − 3; 95% CI − 5.5, − 0.4).

Conclusions

This study provides first evidence that tungiasis is associated with poor neurocognitive functioning in children. Since tungiasis is a chronic disease with frequent reinfections, such negative effects may potentially impair their development and life achievements.

Graphical abstract

Background

Tungiasis is a common parasitic infection in the tropics. The infection occurs when a female sand flea (Tunga penetrans) embeds in the skin. Tungiasis is endemic in South America, the Caribbean, and sub-Saharan Africa and disproportionately affects low-income communities [1]. Particularly children, but also the elderly and individuals with disabilities are the most susceptible to tungiasis [2]. Children are also at critical stage of brain development and are at risk of impaired neurocognitive and mental problems.

Neurocognitive and behavioural development are complex processes that involve distinct yet interrelated bioecological and psychosocial factors. Viral infections like HIV [3] and parasitic infections like schistosomiasis are associated with poor cognitive abilities [4]. Psychosocial factors such as peer, family, and school interactions also exert a significant influence on human cognitive development [5]. Schooling, parenting behaviour, family socioeconomic status, and orphanhood also influence cognitive development in children [6].

Pathogenesis of tungiasis makes it plausible that it can lead to adverse neurocognitive and mental health outcomes. Inflammation at the entry site and bacterial superinfection [7] cause pain and impair mobility and social interactions [8]. Since learning occurs through reciprocal interaction of the individual and the environment [9], impaired interaction puts children with tungiasis at risk for poor neurocognitive development. Although the pathogenesis of tungiasis may contribute to adverse neurocognitive and mental health outcomes among affected persons, there is no implicating empirical evidence. This study investigated the impact of tungiasis on neurocognitive and mental health in school-going children.

Methods

Study design and setting

Community-based cross-sectional surveys were implemented as part of a larger study investigating the disease ecology of tungiasis in Matuga and Msambweni sub-counties in Kwale, Kenya; Ugenya sub-county in Siaya; and Bugiri in Uganda. The regions have various cultures and ethnicities, including livestock-keeping practices, soil features, and closeness to animal habitats, yet their climate conditions are comparable. The recruitment of participants and data collection was done between February 2020 and April 2021.

Study size

The study aimed to test if the mean of outcomes were significantly different in the two group. Therefore, a sample size for a two-sample means test was computed. The study used category fluency as reference outcome. Previous study in a similar setting but with younger population (mean age = 5.2 years) reported a mean category fluency of 15.97 [10]. Assuming a common standard deviation of 2 the study required at least 506 (253 infected and 253 uninfected) participants to detect at least 0.05 difference in means at α = 0.05 and power of 0.8. The sample size was calculated using Stata [11]. However, the actual study size was 454 (220 infected and 234 uninfected), mean difference of category fluency between the two group was 2.6 and a common standard deviation of 7, giving the study a power of 0.97.

Study population and sampling procedure

The study targeted eight to fourteen-year-old children, the most susceptible to tungiasis infection [2]. The study established inclusion criteria that specified residency in a household with a natural soil floor, given its known association with increased risk of tungiasis, and the availability of an adult caregiver for informed consent and interviews. In addition to those specified for infected participants, eligibility criteria for uninfected participants also required the absence of infected family members. In stage I, sixteen public primary schools were to be randomly chosen within each region from a list of all existing public primary schools. However due to the low prevalence of tungiasis in some regions, additional 12 schools in Siaya and four schools in Kwale were randomly selected. Moreover, due to the exceptionally low prevalence of tungiasis in the Bugiri region, the decision was made to conclude data collection after surveying only eight schools. As a result, the selection outcome for Stage I comprised 28 schools from Siaya, 20 schools from Kwale, and eight schools from Bugiri. In stage II, up to a maximum of 102 school-going children were quasi-randomly selected in each school. This process resulted in a total of 5331 pupils. Hands, and feet of the 5331 pupils were visually inspected for tungiasis. Out of the 5331, 589 pupils were infected while 4742 pupils were uninfected. In each school, up to 10 infected and 20 uninfected pupils were then quasi-randomly selected as index pupils from those with tungiasis and those without, respectively. This process resulted in selection of 361 from 589 infected pupils and 729 pupils out of 4742 uninfected pupils as index pupils as shown in Fig. 1. These index pupils were to participate in the larger study. In stage III, six infected and six uninfected children were quasi-randomly selected in each school from the pool of index pupils. The section of infected pupils at this stage was based on severity of infection-aiming at three with severe (> 10 fleas) and three with a mild infection (< 10 fleas) where possible. Conversely, the selection of uninfected children was done through a simple random process. Overall, stage III resulted in selection of 253 of the 361 infected pupils and 523 of the 729 uninfected pupils. This total of 506 pupils (253 infected and 253 uninfected pupils) formed the the final study group for neurocognitive and mental health assessments (Fig. 1).

Fig. 1

Participant selection flow diagram. Orange represents infected children; green represents uninfected children. Yes, indicates pupils selected for the next stage. No, indicates pupils excluded from the study

Study variables

Neurocognitive function and mental health problems were main outcome variables. Five domains of neurocognitive functioning were: attention, memory, language, perceptual-motor, and executive function. These domains are detailed in S1 of the Additional file 1. The explanatory variable of interest was tungiasis status. Other explanatory variables included in each model as potential confounders were nutritional status (underweight, stunting, wasting), disability, perinatal complications, residence (Siaya, Kwale, Bugiri), socioeconomic status (SES), school absenteeism (school days missed in the week preceding data collection), orphanhood, household size, ill family member, and both household head-related factors (sex, relation to participant, age) and caregiver-related factors (sex, relation to participant, age, spending time with the participant, exposure to hugging or cuddling, correction method, caregiver depression, and caregiver stress).

Data sources and methods of measurement

Diagnosis of tungiasis and classification of infection

Trained community health workers washed children’s feet to expose embedded fleas. The feet and fingers of the children were visually inspected embedded flea. The selected participants were categorised infected if they had at least an imbedded flea and uninfected if they did not present with a flea. The embedded fleas were manually counted and the infected further categorised into infection status as mild infection if they had less than 10 imbedded fleas and severe infection if they presented with ten or more imbedded fleas.

Neurocognitive and mental health measures

Participants underwent approximately two-hour battery of tests administered by trained research assistants to assess their neurocognitive abilities across multiple domains. Language function was evaluated with the Early Grade Reading Assessment (EGRA) [12] and Category Fluency Test (CFT) [13], while attention was assessed with the Comprehensive Trail-Making Test (CTMT) [14] and Stroop Color and Word Test (SCWT) [15]. Working memory was evaluated using the backward digit span task [16], and fine motor control was assessed with the bead threading test [17]. The battery also included the Early Grade Maths Assessments (EGMA) [18], Standard and Coloured Raven Progressive Matrices (RPM) [19] to evaluate numeracy and nonverbal intelligence, respectively. The lower scores in these cognitive tests indicate poor neurocognitive function. These tests are valid and reliable for assessing their respective domains and have been adapted for use in Kenya [20] and Uganda [10]. The Child Behavior Checklist [CBCL] [21] was used to assess mental health outcomes. In this study, the total score was used to assess mental health problems, with higher scores indicating more problems. The neurocognitive and mental health measures are detailed in S1 in supplementary material.

Covariates

Anthropometric measurements including height, weight, and Mid-arm circumference (MUAC) were used to assess nutritional status. Height was measured using a stadiometer, weight was measured using a calibrated scale, and MUAC was measured using a flexible tape measure. Height-for-age (HAZ) and weight-for-age (WAZ) z-scores were calculated according to the Centers for Disease Control and Prevention (CDC) [22], with HAZ z-scores < − 2 and WAZ z-scores < − 2 indicating stunting and underweight, respectively, while MUAC was used to evaluate wasting. Structured questionnaires were used to collect data on factors associated with poor neurocognitive and mental health outcomes, including disability child perinatal complications, region of residence, child age, child sex, and school grade level. Household socioeconomic status (SES) was assessed using tetrachoric principal component analysis (PCA) and the resulting wealth index was created based on eigenvalue and scree plot as detailed in S2, S3 and S4 in the Additional file 1. Psychosocial covariates including orphanhood, school absenteeism, household size, and caregiver information were also collected. Respondents’ relation to participant, age, spending time with the participant, exposure to hugging or cuddling, correction method, and caregiver mental health (depression, and stress).

Psychosocial covariates covered various topics such as orphanhood (also assessed using binary response options), school absenteeism (measured as the number of days absent from school in the week preceding data collection), household size (classified as either more than 2 adults or less than 2 adults), having a chronically ill family member (also assessed using binary response options), sex of household head, relation of the household head to the child (classified as either "child" or "not child"), age of household head, sex of caregiver, age of caregiver, and the relation of the caregiver to the child. Other factors assessed in the questionnaires included the amount of time the caregiver spent with the child (reported as "a lot of time" or "not a lot of time"), whether the child was hugged or cuddled (assessed using a binary response option of "yes" or "no"), and the caregiver's methods of correcting the child (reported as "beating" or "other methods").

Caregiver mental health was measured using Patient Health Questionnaire-9 (PHQ-9) and Parental Stress Scale (PSS). The PHQ-9 is a self-report questionnaire that measures depression by asking the respondent to rate the frequency of their symptoms over the past 2 weeks on a scale of 0–3 [23]. The questionnaire covers various areas related to depression, and scores range from 0 to 27. In this study depression cutoff was scores of 10 and above. The Parental Stress Scale (PSS), a questionnaire used to assess parental stress, includes 18 items covering various aspects of parenting, and respondents rate how often they experience stress related to each item on a 5-point scale [24]. The total score is calculated by summing the scores, with higher scores indicating higher levels of parental stress. In this study CBCL, PHQ-9 and PSS had acceptable alphas of 0.94, 0.87 and 0.68, respectively an indication of internal consistency.

Measures to address bias and errors

The outcome assessors were distinct from the infection assessors and were kept unaware of the participant's status to minimize potential biases, however, the status could be known in participants with visible signs of infection. Questionnaires and test score sheets were adapted to the Research Electronic Data Capture (REDCap) database [25], hosted at the International Centre of Insect Physiology and Ecology (ICIPE). Responses were recorded on tablets conditionally formatted to ensure the validity, accuracy, and completeness of the data. The study's data collection involved in-person interviews and assessments performed by assessors trained in neurocognitive evaluations. To accommodate diverse local languages, the questionnaires were initially written in English and subsequently translated into Kiswahili, Dholuo, and Kisoga. Information bias was reduced through pretesting, adjustments, and clarifications of the questions before commencing data collection. To promote consistency in the interview process, a workshop was held involving all assessors and investigators to establish a common interpretation of the responses. Moreover, to minimize misclassification bias, the assessors verified the presence of other infection indicators (like pain and itching) in infected.

Statistical analysis

Data were analysed utilizing STATA software version 15 [11]. The difference in distributions between the uninfected group and each infected group, as well as among all groups together, were compared using binomial and multinomial tests, respectively. Continuous data were presented as means and standard errors (SE) if normally distributed or medians and interquartile range (IQR) if skewed. For bivariate analyses of continuous data, Mann Whitney U Test, or the Kruskal–Wallis rank test was used to compare skewed distributions while Two Sample Students t test or analysis of variance (ANOVA) was used to compare means of normally distributed data. Categorical data were presented as frequencies with respective percentages and their proportions compared using the Pearson Chi-Square Test or Fisher's Exact Test. Analysis of covariates (ANCOVA), followed by the Scheffe Test was used for pairwise comparison of adjusted means. Covariates adjusted for included age, sex, grade, nutritional status, care giver education, disability, school absenteeism and SES.

For regression analyses, multilevel mixed effects generalized linear models with gaussian family and identity link were used for bivariable and multivariable analyses with the unique school identifier as a random effect to identify factors associated with neurocognitive and mental health outcome scores in children adjusting for age, sex, and grade as priori confounders. Multivariate analyses were conducted separately for each neurocognitive and mental health outcome. Backward stepwise selection was used to identify the most significant variables for the model. An exhaustive model containing all predictor variables was initially established. The variable with the highest p-value was subsequently eliminated from the model. This iterative process was continued until the stopping criterion (P < 0.05) was reached and the model selected as final.

Results

Characteristics of participants

This study involved 506 participants, however 52 (approximately 10%) had crucial variables like age, sex or grade missing at random and were excluded from analysis. The exclusion resulted in imbalanced infected and uninfected groups. Analysis was done on data from 454 participants, including 234 uninfected children, 109 children with mild infection, and 111 with severe infection (Table 1).

Table 1 Socio-demographic, biological, and psychosocial characteristics of participants—comparison between uninfected, mild disease, and severe disease groups

The three groups; uninfected, mild infection and severe infection were differently disturbed by sex, age, grade and MUAC (P < 0.001) as shown in Table 1. A higher proportion of participants with severe disease were underweight than the uninfected (18.9% vs 9.8%; P = 0.018). Similarly, higher proportion severe infection group had male caregivers [13/111(11.8%)] than infected [12/109 (5.1%)]. Compared to the uninfected, the severe infection group had a higher proportion of caregivers with depression (47.3% vs 35.5%; P = 0.029). Other characteristics of participants are presented in Table 1.

Neurocognitive and behavioural outcomes

In the current study, the infected group had significantly lower scores than the uninfected group in literacy, language, response inhibition, working memory, fine motor control, non-verbal intelligence, and numeracy, as detailed in Table 2. The lower neurocognitive tests scores suggest impaired neurocognitive ability. Infected group had higher mean total problem scores than uninfected group (34.4 vs 32.2) indicating the presence of behavioural and emotional issues. The uninfected, mild infection and severe infections groups had significantly different adjusted means in literacy, language, cognitive flexibility, working memory, fine motor, and behavioral problems after adjusting for included age, sex, grade, nutritional status, care giver education, disability, school absenteeism and SES. Compared to uninfected, mild infection group had significantly lower scores in literacy, language, cognitive flexibility, working memory, and fine motor while severe infection had significantly lower scores in literacy, fine motor, and behavioral problems as shown in Table 3.

Table 2 Comparison of neurocognitive and behavioral outcomes between uninfected and infected groups—summary of main domains and sub-domains scores

Table 3 Comparison of neurocognitive and behavioral outcomes among uninfected, mild tungiasis, and severe tungiasis groups—adjusted means and group differences corrected for effects of covariates

Covariates of neurocognitive and mental health outcomes

In the bivariable regression analyses, mild infection was associated with lower scores in literacy, language, cognitive flexibility, and working memory, while severe disease was associated with lower scores in fine motor control and higher scores in behavioural problems (Table 4) Other factors associated with lower neurocognitive and behavioural problem are shown in Table 4. After controlling for covariates in multivariable analyses, the mild infection was associated with significant lower scores in literacy [adjusted beta co-efficient (aβ) = − 8.9; 95% confidence interval (CI): − 17.2, − 0.6] where on average, there was 8.9 unit decrease in literacy score in mild infection compared to uninfected as shown in Table 5. Similarly, mild infection was associated with 1.7 (aβ = − 1.7; 95% CI − 3.2, − 0.3), 6.1 (aβ = − 6.1; 95% CI − 10.4, − 1.7) and 0.3 (aβ = − 0.3; 95% CI − 0.6, − 0.1) unit decrease in language cognitive flexibility and working memory, respectively (Table 5). Severe infection was associated with significant lower scores in response inhibition, fine motor control, and numeracy. Averagely, severe infection was associated with a 2.2 unit decrease in response inhibition (aβ = − 2.2; 95% CI − 4.2, − 0.2) a 0.7 unit decrease in fine motor control (aβ = − 0.7; 95% CI − 1.1, − 0.4) and a 3 unit decrease in numeracy (aβ = − 3; 95% CI − 5.5, − 0.4) (Table 5). Stunting was associated with lower scores in language (aβ = − 2.0; 95% CI − 3.9, − 0.2) and numeracy (aβ = − 5.7; 95% CI − 9, − 2.4). Unexpectedly wasting and perinatal complications were associated with better response inhibition scores. Other factors independently associated with lower scores in various domains included residing in a specific geographic region, belonging to households of larger size, having a chronically ill family member, relation to the household head, caregiver depression and caregiver stress (Table 5). After adjusting for covariates, mental health outcomes were neither associated with mild nor severe infection. However, residing in Kwale (aβ = 9.1; 95% CI 3.5, 14.7), having an ill family member (aβ = 4.5; 95% CI 0.7, 8.2), caregiver depression (aβ = 11.4; 95% CI 7.7, 15.2), and caregiver stress (aβ = 0.4; 95% CI 0.1, 0.6) remained significantly associated with mental health outcomes.

Table 4 Beta coefficients (β) and 95% confidence interval (CIs) for the association of characteristics of participants with neurocognitive and behavioral outcomes (N = 454)

Table 5 Adjusted beta coefficients (aβ) and 95% confidence intervals (CIs) for the association of characteristics of participants with neurocognitive and behavioral outcomes (n = 454)

Discussion

The primary objective of this study was to evaluate the association between tungiasis and neurocognitive and mental health outcomes in school-aged children. No previous studies have reported the effect of T. penetrans infections on neurocognitive and mental health outcomes using validated assessment tools as employed in the current investigation. This study analyses effects of any ectoparasitosis on neurocognitive and behavioural functions in children. This gap of knowledge, which the present study pioneers to close, is particularly problematic since it leads to underestimation of effects parasitic diseases that are highly abundant in poor communities, particularly in the tropics. Similar studies targeting other skin disease such as scabies, pediculosis, and even cutaneous larva migrans caused by hookworms, preferably based on a common set of tools to access neurocognitive and behavioural outcomes, would help recognise the impact of ectoparasitosis on health in children.

A significant association between tungiasis and various neurocognitive domains in school-aged children was observed. Even after adjusting for potential confounders, tungiasis remained significantly associated with poor literacy, language, cognitive flexibility, response inhibition, working memory, fine motor skills, and numeracy scores. In the present study, various other factors were associated with poor neurocognitive outcomes. These included stunting, residing in a specific geographic region, belonging to households of larger size, having a chronically ill family member, relation to the household head, and poor mental health among caregivers. These are well documented risk factors for poor neurocognitive outcomes and for children with tungiasis these factors exacerbate the negative impact of tungiasis on outcomes.

Effects of tungiasis might be either direct (e.g., because pain and itching hamper concentration in school) or indirect (e.g., due to stigmatisation). Potential direct effects due to tungiasis-related pain and itching are obvious but have not been investigated directly. It is well established that pain can affect multiple neurocognitive and behavioural functions [26]. The tungiasis-related stigma has been suggested to negatively impact social interaction and participation in educational activities among affected children. A study in Kenya reported that children with tungiasis experience difficulty in borrowing books from their peers and catching up on missed schoolwork [27]. These findings may partially explain the poor performance in language, attention, memory, perceptual and motor control, and executive functions observed among children with tungiasis in the present study.

We are only beginning to understand the underlying pathophysiology of the impact of tungiasis on neurocognitive ability and mental health in children. However, several hypotheses have been formulated. One possibility is that tungiasis-induced persistent pain may lead to deficits in cognitive flexibility and working memory in children [28]. Furthermore, intense itching and pain caused by tungiasis have been postulated to be linked to poor concentration and sleep disturbances [8]. The relationship between sleep disturbance and cognitive impairment is well-documented [29]. Significant poor performance even after accounting for confounders suggests that other mechanisms, such as upregulation of some host immune responses especially during acute phases and mild stages, are at play. Chronicity and severity of infection, in this case, may desensitize [30] or exhaust [31] these immune mediators, which may partially explain why some neurocognitive scores such as cognitive flexibility and working memory were negatively associated with mild infections and insignificantly associated with severe infections (inferring chronicity). However, none of the published studies directly investigated effects of tungiasis on cognitive and behavioural aspects. The present study found no statistically significant association between tungiasis and mental health outcomes after controlling for other covariates. These results contrast with previous research on other chronic skin disease (atopic dermatitis, psoriasis and vitiligo) which showed a relationship between disease and mental health outcomes [32].

Overall, the findings support the hypothesis that tungiasis has a considerable negative impact on multiple neurocognitive functioning in children and may contribute to neurocognitive impairment. The potential long-term effects of such impairment may include difficulty in learning and performing academic tasks, reduced productivity and earning potential in adulthood, and lower overall quality of life. Of particular interest is the bidirectional relationship between tungiasis neurocognitive and mental health outcomes. The study not only highlights the potential for tungiasis to detrimentally affect cognition and behaviour but also invites consideration of the reverse scenario, where cognitive and behavioural aspects might exert an influence on the occurrence and severity of tungiasis. While this bidirectional perspective adds complexity to the interpretation of the current findings, the study design and analytical approach were primarily geared towards exploring the impact of tungiasis on neurocognitive and mental health outcomes. There is need for further research to understand the mechanisms underlying the observed associations, and to develop effective interventions to mitigate the impact of tungiasis on neurocognitive functioning and overall health in affected children.

One of the limitations of the present study is the fact that presence of other pathogens was not accessed, including parasites such as schistosomes, which are known to be prevalent in the study population and have been previously shown to have cognitive effects [4]. Due to the cross-sectional study design, causality cannot be inferred from the results obtained. Bidirectional causality between tungiasis and neurocognitive function was not thoroughly investigated, with our study primarily focussing on assessing the impact of tungiasis on neurocognitive mental health outcomes. Therefore, future longitudinal research including also potential effects of other infectious diseases is recommended to confirm the effects observed here. Utilizing a combination of neuropsychological and neurophysiological measures would also provide a more comprehensive understanding of the short- and long-term effects of tungiasis on cognitive outcomes in school-aged children. It would also be beneficial to investigate potential mechanisms, such as sleep and attention, through which tungiasis may directly impact cognitive ability.

Conclusions

This study uncovers the profound detrimental association between tungiasis and the neurocognitive abilities of school-aged children, highlighting a previously overlooked link between this parasitic disease and poor performance across multiple domains. While emphasizing the urgent requirement for comprehensive interventions and targeted support systems, it is important to acknowledge the limitation of not thoroughly investigating the potential bidirectional influence between tungiasis and neurocognitive function. This calls for further research to understand the underlying mechanisms and develop effective interventions aimed at mitigating the impact of tungiasis on children's neurocognitive functioning. The findings strongly advocate for heightened awareness, improved healthcare measures, and resource allocation to effectively address the far-reaching consequences of tungiasis. Failure to address this critical issue perpetuates the cycle of underestimating the devastating effects of parasitic diseases on the cognitive development and overall well-being of vulnerable communities.