Hostname: page-component-7c8c6479df-xxrs7 Total loading time: 0 Render date: 2024-03-28T11:15:19.772Z Has data issue: false hasContentIssue false

Individual and contextual risk factors for chikungunya virus infection: the SEROCHIK cross-sectional population-based study

Published online by Cambridge University Press:  03 May 2018

A. Fred
Affiliation:
INSERM CIC1410, CHU Réunion, Saint Pierre, Reunion, France Department of Social and Preventive Medicine, School of Public Health, Montreal University, Montreal, Canada
A. Fianu
Affiliation:
INSERM CIC1410, CHU Réunion, Saint Pierre, Reunion, France
M. Béral
Affiliation:
UMR 1309 CMAEE ‘Contrôle des maladies animales, exotiques et émergentes’, CIRAD, Sainte Clotilde, Reunion, France Ministère de l'Agriculture, Direction Régionale de l'Agriculture et de la Forêt, Dijon, France
V. Guernier
Affiliation:
Université de La Réunion, CRVOI ‘Centre de Recherche et de Veille de l'océan Indien’, CYROI, Sainte Clotilde, Reunion, France Australian Institute of Tropical Health and Medicine, James Cook University, Townsville, Queensland, Australia
D. Sissoko
Affiliation:
Department of Tropical Medicine and Clinical International Health, CHU Bordeaux, Bordeaux, France Infectious Diseases in Low Income Countries (IDLIC), Bordeaux Population Health Research Centre (INSERM U1219, Université de Bordeaux, ISPED), Bordeaux, France
M. Méchain
Affiliation:
Department of Tropical Medicine and Clinical International Health, CHU Bordeaux, Bordeaux, France Infectious Diseases in Low Income Countries (IDLIC), Bordeaux Population Health Research Centre (INSERM U1219, Université de Bordeaux, ISPED), Bordeaux, France
A. Michault
Affiliation:
Bacteriology, Virology and Parasitology lab, CHU Réunion, Saint Pierre, Reunion, France
V. Boisson
Affiliation:
Polyvalent Intensive Care Unit, CHU Réunion, Saint Pierre, Reunion, France
B.-A. Gaüzère
Affiliation:
Polyvalent Intensive Care Unit, CHU Réunion, Saint Denis, Reunion, France
F. Favier
Affiliation:
INSERM CIC1410, CHU Réunion, Saint Pierre, Reunion, France
D. Malvy
Affiliation:
Department of Tropical Medicine and Clinical International Health, CHU Bordeaux, Bordeaux, France Infectious Diseases in Low Income Countries (IDLIC), Bordeaux Population Health Research Centre (INSERM U1219, Université de Bordeaux, ISPED), Bordeaux, France
P. Gérardin*
Affiliation:
INSERM CIC1410, CHU Réunion, Saint Pierre, Reunion, France Université de La Réunion, CNRS 9192, INSERM U1187, IRD 249, CHU Reunion, UM 134 PIMIT ‘Processus Infectieux en Milieu Tropical’, CYROI, Sainte Clotilde, Reunion, France
*
Author for correspondence: P. Gérardin, E-mail: patrick.gerardin@chu-reunion.fr
Rights & Permissions [Opens in a new window]

Abstract

The purpose of the study was to weigh the community burden of chikungunya determinants on Reunion island. Risk factors were investigated within a subset of 2101 adult persons from a population-based cross-sectional serosurvey, using Poisson regression models for dichotomous outcomes. Design-based risk ratios and population attributable fractions (PAF) were generated distinguishing individual and contextual (i.e. that affect individuals collectively) determinants. The disease burden attributable to contextual determinants was twice that of individual determinants (overall PAF value 89.5% vs. 44.1%). In a model regrouping both categories of determinants, the independent risk factors were by decreasing PAF values: an interaction term between the reporting of a chikungunya history in the neighbourhood and individual house (PAF 45.9%), a maximal temperature of the month preceding the infection higher than 28.5 °C (PAF 25.7%), a socio-economically disadvantaged neighbourhood (PAF 19.0%), altitude of dwelling (PAF 13.1%), cumulated rainfalls of the month preceding the infection higher than 65 mm (PAF 12.6%), occupational inactivity (PAF 11.6%), poor knowledge on chikungunya transmission (PAF 7.3%) and obesity/overweight (PAF 5.2%). Taken together, these covariates and their underlying causative factors uncovered 80.8% of chikungunya at population level. Our findings lend support to a major role of contextual risk factors in chikungunya virus outbreaks.

Type
Original Paper
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted re-use, distribution, and reproduction in any medium, provided the original work is properly cited.
Copyright
Copyright © Cambridge University Press 2018

Introduction

Chikungunya virus (CHIKV) is responsible for chikungunya fever (CHIKF), a dengue-like illness characterised by persistent incapacitating polyarthralgia [Reference Staples, Breiman and Powers1,Reference Gérardin2]. Between its first description in Tanzania in the early 50s [Reference Lumsden3], and its re-emergence in Lamu island (Kenya) in 2004 and subsequent spread in the Indian ocean in 2005–2006 [Reference Renault4], CHIKF was regarded as a rather harmless neglected tropical disease limited to developing countries. Since the occurrence of atypical and severe forms of CHIKF [Reference Gérardin5Reference Josseran8], and its more recent spread to Europe [Reference Rezza9] and to the Americas [Reference Leparc-Goffart10], CHIKV has increasingly been recognised as a leading threat for public health, through the global expansion of its mosquito vectors, Aedes (Ae) albopictus and Ae aegypti.

Between March 2005 and August 2006, Reunion island experienced a major outbreak of CHIKF [Reference Gérardin11] due to the circulation of a new genotype of the African variant of the virus, renamed Indian ocean lineage, together with the permissive adaption of the local competent vector, Ae albopictus [Reference Schuffenecker12]. At the end of the epidemic, a population-based serosurvey highlighted a huge penetration of the virus in the community with a burden of ~266 000 clinical cases (attack rate: ~35%) and 300 000 infected persons (prevalence rate: 38.2%) [Reference Renault4,Reference Gérardin11]. So far, the reasons why transmission was halted at a prevalence of ~40% remain to be elicited, as such a figure contradicts prediction of herd immunity that estimated the protective level at above 70% [Reference Boëlle13].

Indeed, the understanding of the community burden of CHIKF is a puzzling problem and constitutes a challenge for epidemiologists and public health stakeholders aiming to implement control measures in the event of future outbreaks.

So far, a few population-based cross-sectional serosurveys have proposed the investigation of CHIKF determinants in distinguishing individual and contextual risk factors [Reference Sissoko14Reference Azami17]. None of these has been able to guide decision-making based on appropriate public health impact measures such as population attributable fractions (PAF). With the exception of one study [Reference Feldstein18], all previous studies have based identification of CHIKF determinants on logistic regression approach taking prevalence odds ratio, best known as the odds ratio (OR), as the ‘effect measure’ of interest [Reference Skov19]. These however represent only proxies of the log-binomial model and prevalence proportion ratio (PPR), the gold standard method and estimator in cross-sectional studies, because the OR often tends to overstate the PPR when the outcome is common [Reference Skov19], and hence the OR is more difficult to interpret (i.e. the former can be interpreted only in terms of strength of association, while the latter is multiplicative and can be interpreted in terms of multiplying/dividing the prevalence estimates) [Reference Barros and Hirakata20]. In this context, given the propensity of the OR to overestimate the PPR [Reference Coutinho, Scazufca and Menezes21], we chose the Poisson regression model and the incidence rate ratio (IRR), as two proxies of the log-binomial model and PPR, respectively. We also used PAF, free of hypotheses on causation, with the aim of uncovering risk patterns hosting the determinants of increased CHIKV susceptibility to be investigated in further study of causal pathways.

The objective of the present study was to weigh the community burden of CHIKF risk factors using IRR and related PAF. Determinants were identified to be from primarily individual or contextual origin, and investigated separately before being pooled. This purpose was chosen to guide further prioritisation of the most effective public health interventions for mitigating future outbreaks, based on amendable determinants.

Methods

The SEROCHIK cross-sectional serosurvey was conducted between 17 August and 20 October 2006, shortly after the end of the chikungunya epidemic on Reunion island.

A brief description of the study setting is displayed in Supplementary Appendix 1.

Survey design and procedures

The number of subjects needed to obtain a representative random sample of the general population was estimated to be 2640, assuming 35 ± 2% expected prevalence [Reference Renault4], an α risk of 5%, 20% proportion of refusal or absenteeism and a cluster effect.

The sample was built using a two-level probability sampling procedure. At level 1, a random draw for households (primary sampling unit) was carried out. Selection was conducted in six putative strata defined by the crossed combinations of the dwelling type (individual/collective 2–20 units/collective >20 units) and size of municipality (⩽10 000/>10 000 inhabitants). Overall, five strata were obtained after removing an empty stratum. At level 2, following a predefined Kish method [Reference Kish, Kalton and Heeringa22], the field investigator randomly assigned one ‘index person’ (the Kish individual) per household to be interviewed within the selected dwellings. In accordance, a sampling fraction ranging from 0.071 to 1 was generated.

To ensure that the sample was representative of the underlying population, we corrected the sample for age, gender, dwelling and residence area based on socio-demographic information provided by the 2006 census data. Thus, a set of weights was used to increase or decrease the influence of under-/over-represented groups selected from the source population.

Data collection

Structured questionnaires, administered by 25 field investigators, aimed at collecting data on demographics, health, knowledge on transmission and prevention of CHIKF, dwelling and close environment of residence. These questionnaires are listed in Supplementary Appendix 2.

Individual characteristics included gender, age, occupation, chronic disease, body mass index (BMI), four items dedicated to the knowledge of CHIKV transmission (scored 0–4) and 11 preventive behaviours against Ae albopictus (scored 0–11).

Contextual household variables included the type and altitude of dwelling, household size, area of residence, recent history of CHIKF in the neighbourhood and a social deprivation proxy-index characterising the socio-economic status of the municipality at ecological level.

The construction of this deprivation index is detailed in Supplementary Appendix 3.

Outcome measure

Each participant consented to a fingertip prick and the outcome measured was the result (positive or negative) of CHIKV-specific IgG ELISA serology detected on filter paper-absorbed blood (Whatman 903® Protein Saver™ Card, Schleicher & Schuell) [Reference Gérardin11]. The full procedure of blood testing is detailed in Supplementary Appendix 4.

Additional climate variables

After the geo-referencing of each participant's address at IRIS (‘Ilots Regroupés pour l'Information Statistique’) level, the smallest geographical unit for acquiring aggregated socio-environmental information in France, monthly meteorological data gathered from the 21 closest weather stations were obtained to complete the questionnaire data, using website (http://www.margouilla.net).

The climatic parameters considered were ‘cumulated rainfall’, average, maximal and minimal temperatures, and average solar radiation. The rationale for choosing these variables is explained in Supplementary Appendix 5.

For CHIKV-infected (CHIK+) symptomatic subjects, we inputted the average values of the closest station from the data recorded in the 2 months preceding the first CHIKF case included in the IRIS. For CHIKV-naïve (CHIK−) individuals, we used the average values of the closest station between the first and the last cases reported in the same IRIS.

Statistical analysis

Given the subjective nature and need for intelligibility of some individual variables, comparison of outcome was restricted to people aged 15 years or more to ensure reliable data input based on personal information. First, we examined associations between CHIKF (design-based seroprevalence weighted by sampling fraction) and both individual and contextual characteristics. Design-based crude IRR and 95% confidence intervals (CI) were assessed using χ 2 likelihood ratio tests weighted by sampling fraction.

Second, we investigated independent risk factors for CHIKF in a four-step multivariate Poisson regression for dichotomous outcome, modelling the fixed effects (no random intercept to allow contextual effects). In the first step, we fitted two distinct full models with respect to the intrinsic nature of the variables, individual or contextual (grouping within households) with all the relevant covariates found in bivariate analysis. In the ‘individual-variable’ model, we forced gender and age, potentially associated with unadjustable determinants to minimise residual confounding. In the second step, we fitted one unique parsimonious, ‘explicative’ model with all the putative risk factors previously identified using a manual stepwise backward elimination procedure (P value < 0.05 to be retained in the model). All interaction terms between the variables included in this model were tested using the Mantel–Haënszel method [Reference Greenland, Rothmann, Lash and Greenland23]. In the third step, we studied the contribution of climate variables. Finally, in the fourth step, we added the significant climate variables to the precedent ‘explicative’ model. All the covariates within this final ‘decision-making’ model were set as binary. At each step, we determined the pseudo-likelihood of the Poisson regression model using a survey-adjusted Wald test.

For all these analyses, our purpose was to weigh the ‘explicative’ part of each model using the combined adjusted PAF for cross-sectional data, which were produced from each specific PAF value (for overall PAF calculations, see Supplementary Appendix 6). PAF indicates the percentage of cases that would not occur in a population if the factor was eliminated. Subsequently, we gauged each key determinant of the final ‘decision-making’ model according to its amendable potential, and classified the risk factors into amendable and non-amendable, if they appeared to be causal.

All the analyses were performed using Stata14® (StataCorp. 2015, College Station, TX, USA) excluding observations with missing data.

Ethics and funding

The SEROCHIK serosurvey was approved by the ethical committee for studies with human subjects of Bordeaux (No 2006/47) and the National Commission for Informatics and Liberty, the French Data Protection Authority. All participants provided their informed consent to answer the questionnaire and for blood collection. Parent or guardian of all child participants provided informed consent on their behalf.

The study was funded by the National Institute of Health and Medical Research. The funding source had no role in study design, data collection, analysis and interpretation. The study respected the STROBE statement (Checklist in Supplementary Appendix 7).

Results

Study population

The selection of the study population is presented in Figure 1. Of the 3032 subjects sampled at level 1, 2442 Kish individuals were surveyed by field investigators and 2101 participants aged 15 years or more were analysed. This sample was representative of the habitat in terms of areas of residence and altitude of dwelling place. Beyond these contextual features, it was also representative of individual characteristics such as obesity, hypertension, diabetes and asthma. Of note, our population was shifted towards the over-representation of females (59.8% vs. 49.3%, P < 0.01), elder people (⩾ 50 years, 40.1% vs 33.4%, P ⩽ 0.01) and residents of individual house (78.5% vs 56.4%, P < 0.01) (Table 1).

Fig. 1. Study population. Flow chart of the study population.

Table 1. Characteristics of the 2101 subjects (⩾15 years) analysed for chikungunya risk factors related to population structure and missing data, SEROCHIK study, August–October 2006, Reunion island

a Institut National de la Statistique et des Etudes Economiques.

b 2006 census.

c La Réunion island is divided into four administrative residential areas (microregions): North, West, South and East.

Data are given as numbers and percentages in parentheses.

Individual determinants of chikungunya

The subjects aged 60–69 or 70 years or more, overweight (25 ⩽ BMI < 30 kg/m2) or obese (BMI ⩾ 30 kg/m2), reporting a chronic disease or no occupation were more likely to be infected (Supplementary Table S1). Among active people, farm workers (<10% of the population) were more likely to be infected (P = 0.01). Self-use of repellent sprays or creams (used within more than half of the population) was protective against CHIKF (P < 0.01). Indoor insecticide users (~40% of the population) were slightly protected compared with non-users (P = 0.04). Paradoxically, infection was more common in subjects reporting one or more of the following altruistic behaviours aimed at reducing the amount of breeding sites: covering tanks and water supplies (P = 0.01), putting sand in containers (P = 0.02), pruning shrubs and cleaning wastes (P < 0.01), removing clutter in the courtyard (P < 0.01). Other protective behaviours against day-biting Ae albopictus adult females were not associated with CHIKF (data not shown). However, when adding 11 of the 16 preventive behaviours in a score, we observed a negative correlation between the reported behaviours and exposure level, the subjects cumulating 6–11 protective behaviours being the more likely to be infected (Supplementary Table S1).

Importantly, scoring knowledge of CHIKV transmission by summing good answers to four key questions revealed higher risk for the subjects denying Ae albopictus as the main vector. Interestingly, this score was driven by literacy (P < 0.01), place of birth (P < 0.01) and marital status (P = 0.02) (Supplementary Table S2). Knowledge features showed a protective effect restricted to individuals aged <50 years (data not shown).

Finally, when adjusting eight individual covariates together (forcing gender but excluding the behaviour score), only the absence of occupation was associated with CHIKF (adjusted IRR 1.42, 95% CI 1.20–1.66), even though obesity (adjusted IRR 1.29, 95% CI 1.05–1.58) and poor knowledge on CHIKV transmission (score = 1: adjusted IRR 1.27; 95% CI 1.01–1.59; or score = 0: adjusted IRR 1.33; 95% CI 1.05–1.69) remained linked to CHIKF (Supplementary Table S3).

Taken together, the ‘individual-covariate full model’ and its underlying causative factors uncovered 44.1% of CHIKF (95% CI 41.5–46.6%). In other words, if CHIKF had been fully preventable, control measures aimed at protecting the specific groups of individuals presenting the eight abovementioned risk factors would have diminished the disease burden by more than 40%.

Contextual determinants of chikungunya

Individual house or dwelling located under 750 m height (peak of prevalence observed in the 250–500 m stratum), or in the eastern, southern or western microregions (peak of prevalence in the eastern rainy microregion), were associated with infection (Supplementary Table S4).

Higher risks were observed in people living alone or in households of at least five persons, or reporting a CHIKF history in the neighbourhood. In addition, makeshift constructions (wooden creole house under sheet) and absence of windows screen increased exposure (P < 0.01 for both factors, respectively). At least, risk was increased in municipalities of low-to-intermediate socio-economic level (both at P < 0.001). When integrating climate variables in a model controlling altitude and area of residence (Supplementary Table S5 and Fig. S6), maximal temperature in the month preceding introduction of CHIKV (Tmax-1) in the IRIS was predictive of infection (adjusted IRR 1.14; 95% CI 1.10–1.19) while rainfall (Pluv-1) was very marginal (adjusted IRR 1.1001; 95% CI 1.1000–1.002).

However, when controlling eight of the abovementioned covariates, while individual house (adjusted IRR 1.67; 95% CI 1.34–2.08), CHIKF history in the neighbourhood (adjusted IRR 1.92; 95% CI 1.41–2.60), low-to-intermediate socio-economic level (intermediate: adjusted IRR 1.31; 95% CI 1.05–1.64; deprived: IRR 1.33; 95% CI 1.12–1.56), Tmax-1 (>28.5 °C: adjusted IRR 1.59; 95% CI 1.37–1.83) and Pluv-1 (>65 mm: adjusted IRR 1.31; 95% CI 1.15–1.49) remained independently associated with infection, altitude of dwelling, area of residence and household size were no longer linked to CHIKF, except for categories of people living alone (adjusted IRR 1.15; 95% CI 1.01–1.30) or below 250 m (adjusted IRR 1.94; 95% CI 1.10–3.43) (Supplementary Table S7). Adjusting on climate variables unravelled a weak negative interaction between Tmax-1 and altitude (Mantel–Haënszel OR 0.99; P < 0.01) that we chose to ignore in this intermediate analysis.

Taken together, the ‘contextual-covariate full model’ and its underlying causative factors uncover 89.5% of CHIKF (95% CI 83.4–93.3%). In other words, control measures aimed at protecting the specific groups of individuals living under the eight abovementioned contextual conditions would have diminished the disease burden by ~90%.

Herein, and given PAF values, contextual characteristics uncovered a greater proportion of infections (about double) than individual characteristics.

Influence on modelling on key determinants of chikungunya

For explicative purposes, we decided to weigh individual and contextual risk factors in a full model (data not shown). Overall, adjustment for contextual variables altered marginally the strength of individual characteristics (Table 2). Similarly, adjustment for individual characteristics ruled out the ‘household effect’ while the strength of other contextual variables was slightly modified, with the exception of altitude. Indeed, controlling both individual and contextual characteristics, without considering climate variables, exhibited a negative dose–response effect between altitude and CHIKF, with a maximal exposure below 250 m. Of note, this was observed despite a strong synergistic interaction between altitude and dwelling (Mantel–Haënszel OR 2.58; P = 0.01). Thus, individual house increased the risk substantially for those living in the 250–500 m stratum. As expected, combining all relevant risk factors in the ‘explicative’ model allowed the more comprehensive overview of CHIKF at community level. This was evidenced by increasing overall PAF values (Table 2), the seven covariates uncovering 92.2% (95% CI 85.7–95.7%) of infections. Importantly, this model allowed prioritisation of risk factors using PAF values, the main contributors being altitude of dwelling, then by decreasing order, history of CHIKF in the neighbourhood (hereafter designated as clustering), dwelling, socio-economic level of municipality (indicative of social deprivation), occupation, knowledge of transmission and the BMI.

Table 2. Multivariate explicative model of individual and contextual risk factors for chikungunya among 2101 subjects (⩾15 years), SEROCHIK survey, August–October 2006, Reunion island

Bold characters highlight significant risk factors.

a Poisson regression model. Pei: proportion of exposed individuals in the infected population are given as percentages. aIRR: adjusted incidence rate ratios and 95% CI: 95% confidence intervals. Population attributable fractions (PAF) and 95% CI are given as percentages.

b Individual variables are defined for personal (individual) exposures.

c Regular working or studying occupation.

d Score based on four questions (agree/disagree/1 point): Is chikungunya a mosquito-borne virus? Can chikungunya be transmitted by all species of mosquito? Can the mosquito transmit chikungunya to human? Can the human transmit the virus to the mosquito?

e Collective (grouping) variables are defined for contextual (household or area-level) exposures.

f Previous history of chikungunya fever in the neighbourhood indicative of clustering.

g Derived from a homemade socio-economic index categorising the 24 municipalities of the island into tree levels based on three indices: socio-economic composition (three variables), spatial segregation of ethnic minorities (one variable), existence of measures promoting social cohesion (one variable).

Given the possibility that altitude be confounded by temperature, we tested the robustness of our model by adding Tmax-1 and Pluv-1 in a final ‘decision-making’ model. Prior to this model, by pooling the participants reporting a CHIKF history in their neighbours with those for whom this information was missing, we identified the dwelling as an effect modifier of clustering (Mantel–Haënszel OR 2.33; P < 0.01), this latter being relevant in individual house. Overall, this model retained nine significant covariates plus the interaction term (Table 3). Interestingly, the hierarchy of each specific PAF value was drastically modified, and the final ‘decision-making’ model highlighted the prominence of contextual variables, while the overall PAF of the model was slightly diminished, its covariates and underlying causative factors uncovering 80.8% of CHIKF (95% CI 77.7–83.4%). Importantly, only one individual key determinant, the knowledge score, was found to be amendable, whereas three of the six contextual variables seemed relevant to guide the prioritisation of real-time control measures. Thus, we deemed that the community-based interventions around CHIKF clusters, focused on people living in individual houses within intermediate or deprived areas should be the most appropriate control measures for public health action in our context.

Table 3. Multivariate decision-making model of individual and contextual binary risk factors for chikungunya among 2101 subjects (⩾15 years), SEROCHIK survey, August–October 2006, Reunion island

a Poisson regression model. Pei: proportion of exposed individuals in the infected population are given as percentages. aIRR: adjusted incidence rate ratios and 95% CI: 95% confidence intervals. Population attributable fractions (PAF) and 95% CI are given as percentages.

b Individual variables are defined for personal (individual) exposures.

c Regular working or studying occupation.

d Score based on four questions (agree/disagree/1 point): Is chikungunya a mosquito-borne virus? Can chikungunya be transmitted by all species of mosquito? Can the mosquito transmit chikungunya to human? Can the human transmit the virus to the mosquito?

e Collective (grouping) variables are defined for contextual (household or area-level) exposures.

f Previous history of chikungunya fever in the neighbourhood indicative of clustering. Chikungunya in the neighbours ^ dwelling is an interaction term between the two abovementioned covariates.

g Derived from a homemade socio-economic index categorising the 24 municipalities of the island into tree levels based on three indices: socio-economic composition (three variables), spatial segregation of ethnic minorities (one variable) and existence of measures promoting social cohesion (one variable).

Discussion

Here, we report the findings of a large serosurvey conducted on Reunion island, aimed to guide public health interventions for CHIKF control. By combining survey data and timely acquired climate variables, and using a multi-step multivariate analysis with Poisson regression models for dichotomous outcomes, our analysis has identified and weighed the community burden of individual and contextual determinants of CHIKF.

Interestingly, our results, which fulfil the standards of adequate public health impact measures for decision-making purposes, suggest a high disease burden of risk factors of primarily contextual origin, in comparison with that of individual determinants. In addition, we identified the most susceptible contextual variables to public health interventions.

Thus, in a nine-covariate Poisson model with an additional interaction term, the contributors for CHIKF, ranked by decreasing PAF values, were a CHIKF history in the neighbourhood and individual house, high maximal temperature the month preceding the infection, a socio-economically disadvantaged neighbourhood, low altitude of dwelling, high precipitations the month preceding the infection, occupational inactivity, poor knowledge on CHIKV transmission and obesity/overweight. Together, these covariates and their underlying causative factors uncovered over 80% of infections in the community.

Consistently, the addition of the maximal temperature and cumulated rainfall the month before the introduction of CHIKV in the area provided supportive information showing both an interaction and confounding between altitude and temperature. In accordance, the adjustment for climate variables decreased drastically the risk associated with altitude.

Of utmost importance, the overall PAF value of eight relevant contextual covariates and their underlying causative factors uncovered at least 89.5% of the CHIKF observed on the island, while the overall PAF value of eight relevant individual covariates and their underlying causative factors uncovered up to 44.1%.

This substantial role of contextual variables found in Reunion contrasts the results of the serosurvey conducted in Mayotte [Reference Sissoko14]. Herein, the overall PAF value of three independent household features (construction type, household size and Asset Index) was estimated to be 42% and combining the four independent individual characteristics (gender, birthplace, length of schooling, occupation) was estimated to be 89% in the same age population [Reference Sissoko14]. This discrepancy may underscore both the influence of different methods and distinctive contextual and individual characteristics of the disease burden. Even though epidemiologists have long recognised the ‘neighbourhood effects’ of social deprivation and climate as risk factors for vector-borne diseases such as dengue [Reference Brunkard24], health sociologists do emphasise the relevance of both categories of determinants, environmental changes and landscapes being primarily shaped by the habitat characteristics. Dwellings, in turn, harbour distinct populations with different mental (cognitive) representations and different behaviours against vector-borne diseases owed to different cultural or educational backgrounds [Reference Raude and Setbon15].

The major burden of the altitude may underlie several explicative factors closely related, including climate variables, vegetation, host and dwelling densities. Herein, we show that the maximum risk was observed in the 250–500 m range, which coincides with the traditional Reunionese dwelling, the individual house with garden, an ideal niche for Aedes vectors breeding under tropical climate conditions [Reference Raude and Setbon15], given a synergistic interaction between altitude and dwelling [Reference Serpa25]. As expected, the effect of mid altitudes (250–750 m) was confounded by Tmax-1 and Pluv-1 in the nine-covariate final model, which unravels the critical roles of maximum temperature under 250 m and that of rainfall [Reference Serpa25]. Importantly, ~80% of the Reunionese population lives under the height of 250 m along the coastal area of the island. Aedes albopictus was found as high as 1200 m in Reunion [Reference Delatte26] and even as 2100 m in central Nepal [Reference Dhimal27]. Temperature and precipitations play a key role in CHIKV transmission [Reference Perkins28]. Higher temperatures shorten the extrinsic incubation of CHIKV, lead to an increase in Ae albopictus vectors biting frequency and to an extension in their lifespan [Reference Delatte29]. The effect of precipitations is dual [Reference Anyamba30]. Rainfalls create ideal conditions for mosquitoes to spawn [Reference Ditsuwan31], but may decrease host-seeking female abundance [Reference Roiz32]. Drought periods often lead to inadequate water storage practices [Reference Chretien33], whereas wet climates hasten the hatching of Ae albopictus eggs [Reference Dhimal27]. We believe these pathways likely underlie the effects of temperature and rainfall also in our context.

We found that a history of CHIKF in the neighbourhood, a proxy of disease clustering, uncovered a significant proportion of infections at population level. This is especially obvious, as spatiotemporal clustering is a feature of mosquito-borne diseases amenable to vector control interventions [Reference Dhimal27,Reference Anyamba30,Reference Ditsuwan31,Reference Eisen and Lozano-Fuantes34]. It has been well documented, both in urban and rural environments. This neighbourhood effect on CHIKV transmission dynamic was reported in three clusters investigated during the first epidemic wave in Reunion island [Reference Taglioni and Dehecq35], and is supported by mapping of participant's residencies revealing spatial heterogeneity of risk (data not shown). This influence of neighbourhood was later confirmed in Singapore [Reference Tan36]. It is also in agreement with more recent findings from a CHIKF outbreak in a small rural area of Bangladesh [Reference Salje37]. In this study, using Bayesian modelling of transmission, the authors reported that 58% of infections occurred at neighbourhood level.

The impact of individual house was found consistent with the findings of two other studies conducted on Reunion island [Reference Taglioni and Dehecq35,Reference Setbon and Raude38]. Thus, in the three abovementioned clusters involved in the emergence, age, density and surface of the dwelling increased risk [Reference Taglioni and Dehecq35]. Moreover, makeshift constructions and the absence of window screens were associated with infection, consistent with other observations [Reference Sissoko14Reference Moro16].

In our study, people living in socio-economically deprived areas were more likely to be infected, as previously shown for several infectious pathogens including CHIKV [Reference Setbon and Raude38,Reference Farmer39]. Indeed, neighbourhood environments (here assessed at municipality level) may contribute to CHIKV exposure due to residential segregation (i.e. disparities in dwelling places by race/ethnicity or socio-economic position) [Reference Diez-Roux and Mair40], or due to differences in mosquito control effectiveness, this latter being modified by household expenditures for chemically based protective measures, spatial coverage of interventions or community mobilisation (i.e. acceptability or reluctance behaviours) [Reference Thuilliez, Bellia, Dehecq and Reilhes41Reference Anderson43]. This theoretical framework was confirmed in Mayotte where social disparities in infection rate were primarily structured by housing conditions as well as cognitive representations of the disease, ‘legitimate images’ being found in rich urban settings, ‘folk theories’ being found in poor suburban settings [Reference Raude and Setbon15]. Interestingly, the neighbourhood environment was found different between dengue (or Zika) and CHIKF-affected populations in French Guiana and Rio Janeiro, Brazil, where CHIKF tended to impact the most impoverished communities living in overcrowded areas, which may explain the higher basic reproductive number (R0) for CHIKV than for flaviviruses [Reference Bonifay44,Reference Fuller45].

With respect to individual variables, occupational inactivity was a key determinant of infection, which is likely coherent with increasing exposure to mosquito bites in the home outdoor environment. In agreement, outdoor activities (e.g. being employed as farm worker) have been increasingly reported as a risk factor for CHIKF [Reference Azami17]. In our study, farm workers were the most at-risk profession. We believe that the effect of occupational inactivity may thus reflect the seasonal and undeclared character of farm work on Reunion island, or the importance of other on-off outdoor jobs in a community largely affected by unemployment.

We found knowledge of the disease having an impact on infection. Knowledge of a disease is believed to drive attitudes, beliefs and practices towards better protection [Reference Fritzell46]. Interestingly, appropriate behaviours were not shown to be protective when knowledge was adjusted. Poor knowledge on CHIKV transmission was relevant among elder individuals (⩾80 years), indigenous, illiterate or living alone, which sheds light on the vulnerability of Creole populations. These data are coherent with the irrational ideas of catastrophe thinking reported as ‘folk theories’ among Indian ocean communities along CHIKF epidemics, where the absence of perceived controllability was the most common belief [Reference Raude and Setbon15,Reference Setbon and Raude38]. They are also in line with the findings of the KAP study in French Guiana, which revealed that adoption of protective behaviours against vector-borne diseases is a multi-factorial process that depends on both socio-economic and cognitive factors [Reference Fritzell46].

Last, obesity and to a lesser extent being overweight were slightly associated with infection, which may likely be the conjunction of several possible factors: lack of education and illiteracy, lesser physical activity, larger body surface, higher Aedes mosquito response to CO2, lactic acid or sweat compounds [Reference Torres-Estrada and Rodriguez47].

Our study has some strengths and limitations. First, it was population-based, and our estimates were weighted on the sampling fraction to ensure the best possible representativeness. Notwithstanding, our sample was a little shifted towards women and elderly people and we had to correct the weighting of the sample upon analysis to minimise the selection bias. Second, we used a Poisson regression model for dichotomous outcomes to overcome convergence problems. In the framework of cross-sectional studies, this user-friendly model has gained credit over the last decade to be increasingly used as a good proxy of the log binomial model, including in seroepidemiologic studies [Reference Feldstein18,Reference Coutinho, Scazufca and Menezes21]. In our study, the necessity to compensate the selection bias using survey-readjusted estimates (reweighted on the sampling plan) precluded us to use with Stata the robust variance option, conditional to this model [Reference Coutinho, Scazufca and Menezes21]. However, the estimates of our final model did not differ from unweighted estimates with robust variance (data not shown), so that it is unlikely that this statistical constraint might have changed the overall sense of our results. Third, the use of a Kish sampling method did not allow assessing intra-household random effects, which likely biased the neighbourhood effect estimates. Given the clustering of CHIKF cases [Reference Setbon and Raude38], we strongly recommend cluster sampling for investigating future outbreaks [Reference Sissoko14]. Fourth, the KAP component of the survey was likely sensitive to social desirability bias [Reference Alobuia48], which may have limited the interpretation of the role of protective behaviours. This was suggested in our study by the absence of positive correlation between the knowledge score and the behaviour score. Alternatively, identification of clusters may have also promoted protective behaviours within the infested areas, which may have skewed the relationship between protective behaviours and infections towards a negative correlation.

In conclusion, our findings lend support to a major role of contextual risk factors in CHIKF outbreaks, which in turn highlights the appropriateness of community-based interventions in socio-economically disadvantaged and vector-friendly neighbourhoods [Reference Zellweger49]. Investigating pathways linked to socio-environmental determinants might thus be useful to unravel CHIKF causative factors and to understanding the drivers of future arboviral emergences.

Supplementary material

The supplementary material for this article can be found at https://doi.org/10.1017/S0950268818000341.

Acknowledgements

The authors are grateful to the population of Reunion island for being so co-operative and to the SEROCHIK full investigation team, especially Christine Catteau Joëlle Perrau, Corinne Mussard, Sylvaine Porcherat and Nadège Naty. The authors thank Coralie Foray and Météo France at the CIRAD for providing climate variables. The authors also thank all their collaborators, especially the French Research coordinating Task Force for Chikungunya and Dengue which commissioned the SEROCHIK survey. Full list of collaborators can be found on Supplementary Appendix 8. The authors dedicate this work to the memory of Dr Alain Michault (1950–2015). An obituary has been published in the year 2016.

Footnotes

*

Deceased.

These authors contributed equally.

References

1.Staples, JE, Breiman, RF and Powers, AM (2009) Chikungunya fever: an epidemiological review of a reemerging infectious disease. Clinical Infectious Diseases 49, 942948.Google Scholar
2.Gérardin, P, et al. (2013) Predictors of chikungunya rheumatism: a prognostic survey ancillary to the TELECHIK cohort study. Arthritis Research and Therapy 15, R9.Google Scholar
3.Lumsden, WHR (1955) An epidemic of virus disease in southern province, Tanganyika. II. General description and epidemiology. Transactions of the Royal Society of Tropical Medicine and Hygiene 49, 3357.CrossRefGoogle ScholarPubMed
4.Renault, P, et al. (2007) A major epidemic of chikungunya virus infection in Reunion Island, France, 2005–2006. American Journal of Tropical Medicine and Hygiene 77, 727731.CrossRefGoogle Scholar
5.Gérardin, P, et al. (2008) Multidisciplinary prospective study of mother-to-child chikungunya virus infections on the island of La Réunion. PLoS Medicine 5, e60.CrossRefGoogle ScholarPubMed
6.Economopoulou, A, et al. (2008) Atypical chikungunya virus infections: clinical manifestations, mortality and risk factors for severe disease during the 2005–2006 outbreak on Réunion. Epidemiology and Infection 137, 534541.CrossRefGoogle ScholarPubMed
7.Lemant, J, et al. (2008) Serious acute chikungunya virus infection requiring intensive care during the Reunion Island outbreak in 2005–2006. Critical Care Medicine 36, 25362541.Google Scholar
8.Josseran, L, et al. (2006) Chikungunya disease outbreak, Reunion Island. Emerging Infectious Diseases 12, 19941995.Google Scholar
9.Rezza, G, et al. (2007) Infection with chikungunya virus in Italy: an outbreak in a temperate region. The Lancet 370, 18401846.Google Scholar
10.Leparc-Goffart, I, et al. (2014) Chikungunya in the Americas. The Lancet 383, 514.CrossRefGoogle ScholarPubMed
11.Gérardin, P, et al. (2008) Estimating chikungunya prevalence in La Réunion Island outbreak by serosurveys: two methods for two critical times of the epidemic. BMC Infectious Diseases 8, 99.CrossRefGoogle Scholar
12.Schuffenecker, I, et al. (2006) Genome microevolution of chikungunya viruses causing the Indian Ocean outbreak. PLoS Medicine 3, e263.Google Scholar
13.Boëlle, PY, et al. Investigating transmission in a two-wave epidemic of chikungunya fever, Réunion Island. Vector Borne Zoonotic Diseases 2008; 8, 207218.Google Scholar
14.Sissoko, D, et al. (2008) Seroprevalence and risk factors of chikungunya virus infection in Mayotte, Indian Ocean, 2005–2006: a population-based survey. PLoS ONE 3, e3066.CrossRefGoogle ScholarPubMed
15.Raude, J and Setbon, M (2009) The role of environmental and individual factors in the social epidemiology of chikungunya disease. Health & Place 15, 689699.Google Scholar
16.Moro, ML, et al. (2010) Chikungunya virus in North-Eastern Italy: a seroprevalence survey. American Journal of Tropical Medicine and Hygiene 82, 508511.CrossRefGoogle ScholarPubMed
17.Azami, NA, et al. (2013) Emergence of chikungunya seroprevalence in healthy Malaysian adults residing in outbreak-free locations: chikungunya seroprevalence results from the Malaysian cohort. BMC Infectious Diseases 13, 67.CrossRefGoogle ScholarPubMed
18.Feldstein, LR, et al. (2016) The first reported outbreak of chikungunya in the US Virgin Islands, 2014–2015. American Journal of Tropical Medicine and Hygiene 95, 885889.Google Scholar
19.Skov, T, et al. (1998) Prevalence for proportion ratios: estimation and hypothesis testing. International Journal of Epidemiology 27, 9195.Google Scholar
20.Barros, AD and Hirakata, VN (2003) Alternatives for logistic regression in cross-sectional studies: an empirical comparison of models that directly estimate the prevalence ratio. BMC Medical Research Methodology 3, 21.Google Scholar
21.Coutinho, LMS, Scazufca, M and Menezes, PR (2008) Methods for estimating prevalence ratios in cross-sectional studies. Revista de Saude Publica 42, 16.Google Scholar
22.Kish, L (1965) Some selection techniques. In Kalton, G and Heeringa, S (eds). Survey Sampling. New York, USA: Wiley and Sons, pp. 440501.Google Scholar
23.Greenland, S (2008) Application of stratified analysis methods. In Rothmann, KL, Lash, TL and Greenland, S (eds). Modern Epidemiology, 3rd edn. Philadelphia, USA: Lippincott Williams & Wilkins, pp. 283302.Google Scholar
24.Brunkard, JM, et al. (2007) Dengue fever seroprevalence and risk factors, Texas–Mexico border, 2004. Emerging Infectious Diseases 13, 14771483.CrossRefGoogle ScholarPubMed
25.Serpa, LL, et al. (2006) Seasonal variation of Aedes aegypti and Aedes albopictus in a city of Southeastern Brazil. Revista de Saude Publica 40, 11011105.Google Scholar
26.Delatte, H, et al. (2008) Geographic distribution and developmental sites of Aede albopictus (Diptera: Culicidae) during a chikungunya epidemic event. Vector Borne Zoonotic Diseases 8, 2534.Google Scholar
27.Dhimal, M, et al. (2015) Risk factors for the presence of chikungunya and dengue vectors (Aedes aegypti and Aedes albopictus), their altitudinal distribution and climatic determinants of their abundance in Central Nepal. PLoS Neglected Tropical Diseases 9, e0003545.CrossRefGoogle ScholarPubMed
28.Perkins, TA, et al. (2015) Estimating drivers of autochthonous transmission of chikungunya virus in its invasion of the Americas. PLoS Currents 7, pii: ecurrents.outbreaks.a4c7b6ac10e0420b1788c9767946d1fc.Google Scholar
29.Delatte, H, et al. (2009) Influence of temperature on immature development, survival, longevity, fecundity, and gonotrophic cycles of Aedes albopictus, vector of chikungunya and dengue in the Indian Ocean. Journal of Medical Entomology 46, 3341.Google Scholar
30.Anyamba, A, et al. (2012) Climate teleconnections and recent patterns of human and animal disease outbreaks. PLoS Neglected Tropical Diseases 6, e1465.CrossRefGoogle ScholarPubMed
31.Ditsuwan, T, et al. (2011) Assessing the spreading patterns of dengue infection and chikungunya fever outbreaks in lower southern Thailand using a geographic information system. Annals of Epidemiology 21, 253261.Google Scholar
32.Roiz, D, et al. (2010) Effects of temperature and rainfall on the activity and dynamics of host-seeking Aedes albopictus females in northern Italy. Vector Borne Zoonotic Diseases 10, 811816.CrossRefGoogle ScholarPubMed
33.Chretien, JP, et al. (2007) Drought-associated chikungunya emergence along coastal east Africa. American Journal of Tropical Medicine and Hygiene 76, 405407.CrossRefGoogle ScholarPubMed
34.Eisen, L and Lozano-Fuantes, S (2009) Use of mapping and spatial and space-time modeling approaches in operational control of Aedes aegypti and Dengue. PLoS Neglected Tropical Diseases 3, e411.Google Scholar
35.Taglioni, F and Dehecq, JS (2009) L'environnement socio-spatial comme facteur d’émergence des maladies infectieuses: Le chikungunya dans l'océan Indien. EchoGéo 9, 135.Google Scholar
36.Tan, CH, et al. (2011) Entomological investigation and control of a chikungunya cluster in Singapore. Vector Borne Zoonotic Diseases 11, 383390.Google Scholar
37.Salje, H, et al. (2016) How social structures, space, and behaviours shape the spread of infectious diseases using chikungunya as a case study. Proceedings of the National Academy of Sciences USA 113, 1342013425.Google Scholar
38.Setbon, M and Raude, J (2008) Chikungunya on Réunion Island: social, environmental and behavioural factors in an epidemic context. Population-E 63, 555583.CrossRefGoogle Scholar
39.Farmer, P (1996) Social inequalities and emerging infectious diseases. Emerging Infectious Diseases 2, 259269.Google Scholar
40.Diez-Roux, AV, Mair, C (2010). Neighborhoods and health. Annals of New-York Academy of Sciences 1186, 125145.CrossRefGoogle ScholarPubMed
41.Thuilliez, J, Bellia, C, Dehecq, JS and Reilhes, O (2014) Household-level expenditure on protective measures against mosquitoes on the island of La Réunion, France. PLoS Neglected Tropical Diseases 8, e2609.Google Scholar
42.Lutambi, AM, et al. (2014) Clustering of vector control interventions has important consequences for their effectiveness: a modelling study. PLoS ONE 9, e101625.Google Scholar
43.Anderson, N, et al. (2015) Evidence based community mobilization for dengue prevention in Nicaragua and Mexico (Camino Verde, the Green Way): cluster randomized trial. British Medical Journal 351, h3267.Google Scholar
44.Bonifay, T, et al. (2017) Poverty and arbovirus outbreaks: when chikungunya virus hits more precarious populations than dengue virus in French Guiana. Open Forum Infectious Diseases 4, https://doi.org/10.1093/ofid/ofx247.Google Scholar
45.Fuller, TL, et al. (2017) Behavioral, climatic, and environmental risk factors for Zika and Chikungunya virus infections in Rio de Janeiro, Brazil, 2015–16. PLoS ONE 12, e0188002.CrossRefGoogle Scholar
46.Fritzell, C, et al. (2016) Knowledge, attitude and practices of vector-borne disease prevention during the emergence of a new arbovirus: implications for the control of chikungunya virus in French Guiana. PLoS Neglected Tropical Diseases 10, e0005081.Google Scholar
47.Torres-Estrada, JL and Rodriguez, MH (2003) Physico-chemical signals involved in host localization and the induction of mosquito bites. Salude Publica de Mexico 45, 497505.Google Scholar
48.Alobuia, WM, et al. (2015) Knowledge, attitude and practices regarding vector-borne diseases in Western Jamaica. Annals of Global Health 81, 654663.CrossRefGoogle ScholarPubMed
49.Zellweger, RM, et al. (2017) Socioeconomic and environmental determinants of dengue transmission in an urban setting: an ecological study in Noumea. PLoS Neglected Tropical Diseases 11, e0005471.CrossRefGoogle Scholar
Figure 0

Fig. 1. Study population. Flow chart of the study population.

Figure 1

Table 1. Characteristics of the 2101 subjects (⩾15 years) analysed for chikungunya risk factors related to population structure and missing data, SEROCHIK study, August–October 2006, Reunion island

Figure 2

Table 2. Multivariate explicative model of individual and contextual risk factors for chikungunya among 2101 subjects (⩾15 years), SEROCHIK survey, August–October 2006, Reunion island

Figure 3

Table 3. Multivariate decision-making model of individual and contextual binary risk factors for chikungunya among 2101 subjects (⩾15 years), SEROCHIK survey, August–October 2006, Reunion island

Supplementary material: File

Fred et al. supplementary material

Fred et al. supplementary material 1

Download Fred et al. supplementary material(File)
File 1.8 MB