# Integrated Analysis of Environment, Cattle and Human Serological Data: Risks and Mechanisms of Transmission of Rift Valley Fever in Madagascar

Marie-Marie Olive, Véronique Chevalier, Vladimir Grosbois, Annelise Tran, Soa-Fy Andriamandimby, Benoit Durand, Jean-Pierre Ravalohery, Seta Andriamamonjy, Fanjasoa Rakotomanana, Christophe Rogier, et al.

<figure><img src="image_1.png" /></figure>

## ► To cite this version:

Marie-Marie Olive, Véronique Chevalier, Vladimir Grosbois, Annelise Tran, Soa-Fy Andriamandimby, et al.. Integrated Analysis of Environment, Cattle and Human Serological Data: Risks and Mechanisms of Transmission of Rift Valley Fever in Madagascar. PLoS Neglected Tropical Diseases, 2016, 10 (7), pp.e0004827. 10.1371/journal.pntd.0004827. pasteur-01675206

HAL Id: pasteur-01675206

https://riip.hal.science/pasteur-01675206

Submitted on 4 Jan 2018

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

<figure><img src="image_4.png" /></figure>

CrossMark
click for updates

RESEARCH ARTICLE

# Integrated Analysis of Environment, Cattle and Human Serological Data: Risks and Mechanisms of Transmission of Rift Valley Fever in Madagascar

Marie-Marie Olive¹,²*, Véronique Chevalier¹, Vladimir Grosbois¹, Annelise Tran¹, Soa-Fy Andriamandimby², Benoit Durand³, Jean-Pierre Ravalohery², Seta Andriamamonjy², Fanjasoa Rakotomanana⁴, Christophe Rogier⁵, Jean-Michel Heraud²

1 CIRAD, Animal and Integrated Risk Management (AGIRs) Unit, Montpellier, France, 2 Institut Pasteur de Madagascar, Virology Unit, Antananarivo, Madagascar, 3 Paris-Est University, ANSES, Laboratory for Animal Health, Epidemiology Unit, Maisons-Alfort, France, 4 Institut Pasteur de Madagascar, Epidemiology Unit, Antananarivo, Madagascar, 5 Institut Pasteur de Madagascar, Direction, Antananarivo, Madagascar

* marie-marie.olive@cirad.fr

<figure><img src="image_5.png" /></figure>

OPEN ACCESS

**Citation:** Olive M-M, Chevalier V, Grosbois V, Tran A, Andriamandimby S-F, Durand B, et al. (2016) Integrated Analysis of Environment, Cattle and Human Serological Data: Risks and Mechanisms of Transmission of Rift Valley Fever in Madagascar. PLoS Negl Trop Dis 10(7): e0004827. doi:10.1371/journal.pntd.0004827

**Editor:** Christopher M. Barker, University of California, Davis, UNITED STATES

**Received:** January 7, 2016

**Accepted:** June 16, 2016

**Published:** July 14, 2016

**Copyright:** © 2016 Olive et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** All relevant data are within the paper and its Supporting Information files.

**Funding:** This study was supported in part by funds raised by (i) CIRAD (ii) the Centre National d'Etudes Spatiales (CNES; grant # 131163/00) (iii) the Institut Pasteur de Madagascar. Human field work was funded by the Institut Pasteur of Madagascar (Internal Project through the ZORA (ZOonoses, Rodent and Arboviruses) project and the Centers for Disease Control and Prevention (CDC; Cooperative Agreement #U51/IP000327). Moreover, this material

## Abstract

### Background

Rift Valley fever (RVF) is a vector-borne disease affecting ruminants and humans. Madagascar was heavily affected by RVF in 2008–2009, with evidence of a large and heterogeneous spread of the disease. The identification of at-risk environments is essential to optimize the available resources by targeting RVF surveillance in Madagascar. Herein, the objectives of our study were: (i) to identify the environmental factors and areas favorable to RVF transmission to both cattle and human and (ii) to identify human behaviors favoring human infections in Malagasy contexts.

### Methodology/Principal Findings

First, we characterized the environments of Malagasy communes using a Multiple Factor Analysis (MFA). Then, we analyzed cattle and human serological data collected at national level using Generalized Linear Mixed Models, with the individual serological status (cattle or human) as the response, and MFA factors, as well as other potential risk factors (cattle density, human behavior) as explanatory variables. Cattle and human seroprevalence rates were positively associated to humid environments (p<0.001). Areas with high cattle density were at risk (p<0.01; OR = 2.6). Furthermore, our analysis showed that frequent contact with raw milk contributed to explain human infection (OR = 1.6). Finally, our study highlighted the eastern-coast, western and north-western parts as high-risk areas for RVF transmission in cattle.

### Conclusions/Significance

Our integrated approach analyzing environmental, cattle and human datasets allow us to bring new insight on RVF transmission patterns in Madagascar. The association between
is based upon work supported by the U. S.
Department of Homeland Security under the Center
of Excellence of Emerging and Zoonotic Animal
Diseases Grant Award Number 2010-ST061-
AG0001. Cattle study was supported in part by funds
raised by World Health Organization (WHO) and the
Food and Agriculture Organization of the United
Nations (FAO) through the Central Emergency
Response Fund of the United Nations. The views and
conclusions contained in this paper are those of the
authors and should not be interpreted as necessarily
representing the official policies, either explicit or
implicit, of the U. S. Department of Homeland
Security. The funders had no role in study design,
data collection and analysis, decision to publish, or
preparation of the manuscript.

**Competing Interests:** The authors have declared
that no competing interests exist.

cattle seroprevalence, humid environments and high cattle density suggests that concomi-
tant vectorial and direct transmissions are critical to maintain RVF enzootic transmission.
Additionally, in the at-risk humid environment of the western, north-western and the east-
ern-coast areas, suitable to Culex and Anopheles mosquitoes, vectorial transmission prob-
ably occurs in both cattle and human. The relative contribution of vectorial or direct
transmissions could be further assessed by mathematic modelling.

Author Summary

Rift Valley fever virus (RVFV) is a pathogen that causes a vector-borne tropical disease. The disease affects ruminants and humans and severely impacts the health and economy of affected countries. Madagascar was heavily affected by Rift Valley fever (RVF) in 2008–2009, with evidence of a large and heterogeneous spread of the disease. Our study aims at identifying environmental and human-related risk factors for RVFV transmission. First, we characterized Malagasy environments according to their putative influence on RVFV mosquito density and population dynamics. Then we statistically analyzed cattle and human serological data collected at a national level with the individual serological status as response, and Malagasy environments previously characterized by climatic and landscape variables as well as other potential risk factors as explanatory variables. Our results identified humid environments of the western, north-western and eastern parts of the island as risky areas. The identification of at-risk environments is essential to focus veterinary surveillance and control of RVFV.

**Introduction**

Rift Valley fever virus (RVFV) is an arthropod-borne zoonotic virus belonging to the *Bunyaviridae* family and affecting ruminants and humans. Infection causes abortion in pregnant ruminants and acute deaths in newborns [1,2]. In the majority of human cases, infection is asymptomatic or causes mild symptoms such as fever, headaches and muscle pains [2]. However severe cases occur, characterized by retinitis, encephalitis, or hemorrhagic forms that may lead to death [2]. Ruminants are infected through vector bites and probably also by direct contact with infected tissues or fluids, such as blood or abortion products [2,3]. Humans are mainly infected through direct contact with infectious tissues or fluids of ruminants but vectorial transmission has been suspected in Central African Republic (RCA) and Gabon [4,5]. Virus circulation has been reported in several eco-climatic areas: arid in Western Africa and Arabic Peninsula [1,6], sub-humid in Eastern Africa [7,8], wet forests in central Africa [5], dam and irrigated agricultural land under hot climatic conditions in Egypt, Mauritania and Sudan [9–11] and recently humid highlands in Madagascar [3,12]. The respective roles of direct and vectorial transmissions remain unevaluated in both human and cattle and probably vary among these eco-climatic areas.

Madagascar experienced two major Rift Valley fever (RVF) outbreaks: 1990–91 in the eastern-coast and central highlands and 2008–09 in the south, the north and the highlands [13–15]. The last outbreaks occurred in two epidemic waves during the two successive rainy seasons of 2007–08 and 2008–09. Following the first wave, passive surveillance and emergency response were developed. Sentinel surveillance in herds were set up with field veterinarians [16]. This sentinel surveillance allowed the early detection of the second wave of outbreak in
cattle at the end of 2008 and thus the implementation of local control measures to prevent the
spreading of RVF outside the region [15,16]. At the end of the epidemic, about 700 suspected
human cases were recorded from which 26 were fatal. About 400 human and cattle samples
were received for laboratory analyses and RVF infection was confirmed or considered as proba-
ble in 86 human and 46 ruminant samples [15]. Following the 2008–09 epidemics, studies
showed a wide and heterogeneous spread of RVFV infection both in human and cattle [15,17]
suggesting that some areas were more favorable than others to transmission [17]. Madagascar
has a large variety of eco-climatic patterns, including semi-arid in the south, tropical in the
west and on the eastern-coast, and temperate in the central highlands [18]. Apart from the
highlands [3,12,19], RVF epidemiology is poorly understood in this country [15,17]. Since
2007, a human syndromic-based surveillance system has been developed which has allowed
the detection of the first case of RVF in humans in 2008 [15]. Besides, retrospective investiga-
tions suggested that RVFV circulated among livestock since December 2007 [15], revealing a
dearth in veterinary surveillance. The main difficulty to implement veterinary surveillance in
Madagascar is the lack of basic means to collect and communicate veterinary information [20].
Thus, the identification of at-risk environments is essential to optimize the available resources
by targeting RVF surveillance. In addition, there is a need to provide insight into the role of the
two transmission routes and better adapt available control measures.

Herein, the objectives of our study were: (i) to identify the environmental factors and areas favorable to RVFV transmission to both cattle and human and (ii) to identify human behaviors favoring human infections in Malagasy contexts.

**Materials and Methods**

To achieve these goals, we characterized the environments of Malagasy communes using a
Multiple Factor Analysis (MFA). Then we analyzed cattle and human serological data using a
Generalized Linear Mixed Models (GLMMs), with the individual serological status (cattle or
human) as the response, and MFA factors, as well as potential other risk factors (covariates), as
explanatory variables.

Cattle and human datasets

The cattle dataset contained results of a national cross-sectional serological survey performed in August 2008 on 3,450 ruminants [17]. Only cattle with known breeding location were included in the study (n = 1,432; Fig 1; [17]). The human dataset contained data from a national cross-sectional serological survey conducted from November 2011 to April 2012 and from October 2012 to May 2013 in 56 sites (cities or villages). Six percent of these sera were used in a previous study [21]. In each of these sites, 30 adults were randomly chosen and sampled on a voluntary basis (Fig 1). Potential contacts with ruminants or fresh ruminant fluids (secretion, blood, milk) and socio-professional categories—butcher, farmer, health worker, worker in contact with environment (water, forest) and others (teacher, student, administrative worker, retired)—were documented through a dedicated questionnaire.

Ethics statement

The cattle study was performed in collaboration with the Malagasy Veterinary Services and animals were sampled by qualified veterinarians [17]. The human study protocol was approved by the Malagasy competent authorities, the Malagasy Ethic National Committee (authorization N°066/MSAMP/CE, 26th July 2011). After reading of the informed consent letter, written and oral consent was obtained from volunteering individuals. Participants were sampled by qualified investigators and the data were analyzed anonymously.
<figure><img src="image_6.png" /><figcaption>**Fig 1. Cattle and human sampling sites [17].** Animal and human sera were analyzed using commercial ELISA kits (Biological Diagnostic Supplies Ltd., BDSL) to detect anti-RVFV immunoglobulin (Ig) G [17,22,23]. Cattle and human data were aggregated at the commune level (n = 1,578).</figcaption></figure>

doi:10.1371/journal.pntd.0004827.g001

## Covariates

The following covariates were selected according to their putative influence on mosquito density and population dynamics or on the risk of contact with ruminants:

* Cattle density. This variable has previously been identified as a risk factor for RVF transmission [7,24].

* Surface covered by water bodies and landscape classes (such as forest, shrub, and agricultural areas). Density and population dynamics of vectors are influenced by environmental factors such as climate, the presence of water bodies and other landscape features [1,25]. The presence of temporary water bodies and floodplains are known as risk factors for RVF in semi-arid areas in eastern Africa, the Arabian Peninsula and Western Africa [1]. Artificial water bodies such as dam and irrigated rice fields are also known to be associated with high abundance of RVFV vectors in western Africa [1]. Furthermore, RVFV transmission occurred in forested or shrubby areas [5,8,26,27]. A recent study details the mosquito species and their habitat in Madagascar [28]. Briefly, RVF potential mosquito vectors belong to the genera *Aedes*, *Anopheles*, *Culex*, *Eretmapodites* and *Mansonia* [28]. The breeding areas of the *Aedes* genus are mostly associated with temporary water bodies such as flooded area, temporary pond, puddles, rice field [28]. *Culex* and *Anopheles* mosquitoes breeding areas are diversified and could be temporary (rice field, swamps) or permanent (lakes, pond). Stagnant and permanent water bodies are the habitat of *Eretmapodites* and *Mansonia* respectively [28].

* Rainfall, Normalized Difference Vegetation Index (NDVI) and temperatures. The risk of RVFV infection in Eastern and Southern Africa has been shown to vary as a function of rainfall, NDVI and temperatures [29,30].

* Human related factors: habitat, gender, profession, contact with ruminant and ruminant products [31].

Values of covariates were computed at the commune level (except for human behaviors which were at the individual level).
**Cattle density.** For each of the 1,578 communes considered, cattle density was estimated using the new global distribution maps for cattle produced by the Food and Agriculture Organization of the United Nations (FAO; http://www.fao.org/Ag/againinfo/resources/en/glw/GLW_dens.html; [32]).

**Water bodies and landscape classes.** A landscape map of Madagascar was obtained from Globcover project [33]. The GlobCover 2009 landscape product is a 300-m global landscape map produced from an automated classification of Medium Resolution Imaging Spectrometer (MERIS) time series. The global landscape map included 22 landscape classes defined with the United Nations (UN) Land Cover Classification System (LCCS). Among these 22 classes, we identified 5 relevant LCCS categories: “Cultivated Terrestrial Areas and Managed Lands” (so-called Crops), “Woody/ Trees”, “Shrubs”, “Herbaceous”, “Artificial Surfaces (so-called Urbanization)”. To reflect the availability of potential breeding habitats of RVF vectors in Madagascar such as artificial, irrigated, permanent and temporary water bodies, we needed to combine different data sources extracted from several GIS databases. The first one described inland permanent water point, such as lake, and was available from DIVA-GIS (http://www.diva-gis.org/). Marshland data representing temporary water bodies were obtained from Geographical Information Systems at the Royal Botanic Gardens, Kew [34]. Wetland locations representing temporary water bodies were extracted from the International Panel on Climate Change (IPCC; [35]). Irrigated area locations came from Global Map of Irrigation Areas (GMIA) from AQUASTAT-FAO [36].

**Climatic variables: Precipitation, temperatures and NDVI.** To depict the climatic conditions at each commune of Madagascar, day and night Land Surface Temperature (LST) and NDVI were retrieved from Moderate Resolution Imaging Spectroradiometer (MODIS; http://iridl.ldeo.columbia.edu/). For the period 2001 to 2010, day and night LST were extracted from MODIS data produced every 8 days at 1 km spatial resolution (MODIS MOD11A2 product: Land Surface Temperature and Emissivity). For the same period, NDVI data were obtained from MODIS data produced every 16 days at 250 m spatial resolution (MODIS MOD13A1: Vegetation Indices). Rainfall data were retrieved from the Tropical Rainfall Measuring Mission (TRMM; http://pmm.nasa.gov/trmm/mission-end). These data were produced at 25 km spatial resolution. Finally, for each commune and the same period, we computed the annual mean of day and night LST, NDVI and precipitation. Seasonality of NDVI and precipitation was also considered by computing the difference between the cumulated value over 3 months of the rainy season (November, December and January) and the cumulated value over 3 months of the dry season (June, July and August).

**Human related factors.** Human density was computed for each commune using data generated from Landscan 2007 Global Population Grid from Oak Ridge National Laboratory & the US Department of Defense (OCHA, 2007). Based on our field knowledge, the communes with more than 5,000 persons per square kilometer were considered as “urban”. Other communes were classified as “rural”. Human behaviors were documented through a dedicated questionnaire.

For each of the 1,578 communes considered, the percentage of surface covered by each landscape class (vegetation and water bodies), as well as the values of climatic, NDVI and cattle density covariates were computed with the Quantum GIS software [37]. Malagasy commune administrative boundaries and data come from the layers data merged by the Office for the Coordination of Humanitarian Affairs (OCHA) and based on data obtained from the Malagasy National Disaster Management Office in 2011.

## Multiple Factor Analysis

Synthetic variables characterizing the environment of communes were computed using a MFA combining the previously mentioned climatic and landscape variables [38,39]. By performing a
factor analysis inside each variable category and then between categories, MFA produces a quantitative summary of the initial set of variables taking the form of a set of linear combination of variables, referred to as factors [39]. The climatic category included the annual means of day and night LST, the annual mean and seasonality of precipitation. The landscape category included the percentage of the surface of the commune covered by each landscape category and the annual mean and seasonality of NDVI. The value of each factor was computed for each of the 1,578 Malagasy communes. Correlation between MFA factor values and cattle density distribution was assessed using Pearson product-moment correlation coefficient test.

## Statistical analysis

As a first step univariate analyses of association between suspected risk factors and cattle or human RVFV serological status were undertaken using Chi square tests for categorical factors and generalized linear models for quantitative factors. Risk factors with significance level ≤0.20 were then included as explanatory variables in GLMMs, with cattle or human individual serological status as the binomial response. In these models, it was assumed that the relationships between serological prevalence and quantitative factors were linear on the logit scale. To account for interdependency of serological status of individuals sampled in the same locality, the smallest administrative unit—the commune for the cattle model and the city/village for human model—were included in the models as a random effect. Multicollinearity among explanatory variables was assessed using Variance Inflation Factors (VIF) and correlation tests. Collinear factors were not included in a same model. The selection of the best models was based on the Akaike Information Criterion (AIC). When needed, a multi-model inference approach was used to estimate model-averaged fixed effects (mafe) and the relative importance (RI) of each explanatory variable [40]. Within the set of models tested, only those with an AIC within 2 units difference from the best model were considered [40].

Internal validity of sets of models was evaluated using the Receiver Operating Characteristic (ROC) curve method [41].

In addition, we calculated the 10-fold cross-validation prediction. Because, it is not possible to perform 10-fold cross-validation on GLMM, this procedure was applied to Generalized Linear Models that were similar to the selected GLMM except that did not include the site of sampling as random effect. Firstly, the cattle seroprevalence dataset was split randomly into 10 parts. Then, the model was fitted to 90% of the data and used to predict the serological status of the remaining 10% individuals as validation step. The procedure was performed 10 times, each time with 1 of the 10 parts as validation step. [42].

Finally, parameter estimations derived from the best cattle model were used to predict and map cattle seroprevalence at the commune scale for the whole island.

Data analyses were performed using R software version 3.0.1 [43–49].

## Results

### Environmental characterization of Malagasy communes

Four MFA factors contributing to 60% of the total variance were selected. Table 1 shows the correlation between each quantitative covariate included in the MFA and each of these four factors:

* Factor 1 separated areas based on seasonality in primary productivity (photosynthetic activity measured by NDVI), vegetation, land use and temperature. Large positive values described ecosystems with high seasonal primary productivity dominated by herbaceous vegetation and with low surfaces of crops under dry and hot climatic conditions (Fig 2A in
<table><caption>Table 1. Correlation between each quantitative covariate included in the MFA and each factor (Factor 1, Factor 2, Factor 3 and Factor 4).</caption><thead><tr><td>Covariate</td><td>Factor 1</td><td>Factor 2</td><td>Factor 3</td><td>Factor 4</td></tr></thead><tbody><tr><td>Mean LST-day</td><td>0.92</td><td>-0.19</td><td>0.11</td><td>/</td></tr><tr><td>Mean LST-night</td><td>0.50</td><td>-0.66</td><td>0.14</td><td>0.26</td></tr><tr><td>Mean precipitation</td><td>-0.70</td><td>/</td><td>0.32</td><td>0.31</td></tr><tr><td>Seasonality of precipitation</td><td>0.17</td><td>-0.15</td><td>0.82</td><td>0.09</td></tr><tr><td>Mean NDVI</td><td>-0.83</td><td>-0.34</td><td>/</td><td>/</td></tr><tr><td>NDVI seasonality</td><td>0.63</td><td>0.45</td><td>0.08</td><td>0.08</td></tr><tr><td>Herbaceous</td><td>0.84</td><td>-0.12</td><td>-0.24</td><td>0.11</td></tr><tr><td>Shrubs</td><td>0.11</td><td>0.40</td><td>0.30</td><td>-0.17</td></tr><tr><td>Wood-Trees</td><td>-0.33</td><td>0.56</td><td>0.37</td><td>-0.19</td></tr><tr><td>Urbanization</td><td>/</td><td>0.14</td><td>-0.30</td><td>0.27</td></tr><tr><td>Crops</td><td>-0.62</td><td>-0.61</td><td>-0.24</td><td>0.10</td></tr><tr><td>Irrigated area</td><td>/</td><td>0.66</td><td>-0.08</td><td>0.37</td></tr><tr><td>Wetlands</td><td>/</td><td>0.24</td><td>-0.39</td><td>0.46</td></tr><tr><td>Water bodies</td><td>/</td><td>/</td><td>0.07</td><td>0.22</td></tr><tr><td>Marshlands</td><td>/</td><td>0.07</td><td>0.18</td><td>0.71</td></tr></tbody></table>

/: The correlation coefficients were not significantly different from zero and so not included in the results

doi:10.1371/journal.pntd.0004827.t001

green). Large negative values described ecosystems with low seasonal primary productivity including crops under wet and less hot climatic conditions (Fig 2A in brown). The communes with the largest positive values for Factor 1 are located in the south-western part of Madagascar (Fig 2A in green) while the communes with the largest negative values for Factor 1 are located on the north-eastern part (Fig 2A in brown);

* Factor 2 separated areas based on seasonality in primary productivity, vegetation, land use and temperature. Large positive values described ecosystems with high seasonal primary

<figure><img src="image_7.png" /><figcaption>Fig 2. Geographical representation of the MFA factor values and cattle density of the 1,578 Malagasy communes. (A) Factor 1, (B) Factor 2, (C) Factor 3, (D) Factor 4, (E) cattle density categories. For each factor, green colors represent positive values and brown negative values. The darkest colors represent the highest values. Cattle were sampled in communes surrounded in black and human were enrolled in communes surrounded in purple.</figcaption></figure>

doi:10.1371/journal.pntd.0004827.g002
productivity including ligneous vegetation and irrigated areas (rice fields) under climatic conditions characterized by low night temperatures (Fig 2B in green). Large negative values described ecosystems with low seasonal primary productivity including crops under climatic conditions characterized by warm night temperatures (Fig 2B in brown). The communes with the largest positive values of Factor 2 are located in the central highlands (Fig 2B in green). The communes with the largest negative values are mostly located in the eastern part of the island (Fig 2B in brown);

* Factor 3 was a rainfall seasonality index. The highest values of Factors3 (highly seasonal rainfall) are observed in the north-western part of the island (Fig 2C in green);

* Factor 4 represented a humid areas (marshlands, wetland and irrigated lands) index. The highest values are mostly located on the eastern-coast and the north-western part of the island (Fig 2D in green).

Considering each of the 1,578 communes of Madagascar, MFA factors values ranged from -1.9 to 3.3 (Factor 1), -1.9 to 2.8 (Factor 2), -5.1 to 2.7 (Factor 3) and -1.1 to 7.3 (Factor 4).

## Description of data and univariate statistical analysis (Table 2)

A total of 1,432 individuals from the initial cattle dataset, sampled in 26 communes belonging to 22 Malagasy districts were included in the study (Fig 1). The number of animals sampled per commune ranged from 1 to 110. Cattle ages ranged from 1 to 12 years (mean age 4.5 years). The overall seropositivity rate was 19.3% (CI 95% [17.3–21.8]). Age was categorized in 4 groups: 1–2, 3–4, 5–6 and more than 7 years old. Cattle density was classified as follows according to quartiles: below 9.7; 9.7 to 14.3; 14.3 to 19.1 and more than 19.1 per square kilometer. MFA factor values of the 26 sampling communes ranged from -1.7 to 2.6 (Factor 1), -0.9 to 1.5 (Factor 2), -1.5 to 2.3 (Factor 3) and -1.1 to 0.6 (Factor 4). Age category, cattle density category, Factor 1, Factor 3 and Factor 4 were statistically associated with cattle seroprevalence ($p \le 0.20$). A total of 1,679 people were sampled, 91% (n = 1,529) living in rural areas and 9% (n = 150) living in urban areas (Fig 1). Age of volunteers ranged from 18 to 99 years (mean age 37.6 years) with a ratio of 1.03 (male/female). The overall seropositivity rate was 9.5% (95% CI [8.1–11.0]). Age was categorized in 4 groups: 18 to 26, 26 to 37, 37 to 46 and more than 46 years old. Cattle density of the related communes was classified as following: below 6.3; 6.3 to 11.7; 11.7 to 22.0 and more than 22.0 per square kilometer. A total of 267 individuals declared no contact with live animals or animal product and were categorized as "other profession". Among them, 24 individuals were seropositive (9.0% 95% CI [5.8–13.1]). MFA factor values of the 48 communes ranged from -1.86 to 3.29 (Factor 1), -1.87 to 2.77 (Factor 2), -5.08 to 1.75 (Factor 3) and -0.75 to 4.51 (Factor 4). Habitat, gender, contact with ruminants, contact with raw milk, profession, age category, cattle density category, Factor 2, Factor 3 and Factor 4 were statistically associated with human seroprevalence ($p \le 0.20$).

## Multivariate analysis

Both cattle and human seroprevalences increased gradually with age categories (Table 2). It was thus assumed that the relationship between age and seroprevalence was linear: age was thus included as a quantitative variable in multivariate models. Since the variation in cattle or human seroprevalence along cattle density categories was not clearly gradual, cattle density was included as a categorical variable in both cattle and human multivariate models. Cattle density was correlated with Factor 1, Factor 2 and Factor 3 and thus was included separately
<table><caption>Table 2. Descriptive and univariate analyses for cattle and human seroprevalences.</caption><thead><tr><th colspan="2">Characteristics</th><th></th><th>Positive</th><th>Total</th><th>Seroprevalence [95% CI]</th><th>Chi2</th></tr></thead><tbody><tr><td rowspan="10">Cattle</td><td rowspan="4">Age</td><td>1 to 2</td><td>46</td><td>353</td><td>13.0 [9.7–17.0]</td><td>p&lt;0.001</td></tr><tr><td>3 to 4</td><td>69</td><td>422</td><td>16.4 [13.0–20.2]</td><td></td></tr><tr><td>5 to 6</td><td>72</td><td>361</td><td>19.9 [15.9–24.4]</td><td></td></tr><tr><td>&gt; 7</td><td>90</td><td>296</td><td>30.4 [25.2–36.0]</td><td></td></tr><tr><td rowspan="4">Cattle density per sq. km</td><td>&lt; 9.7</td><td>69</td><td>359</td><td>19.2 [15.3–23.7]</td><td>p &lt; 0.001</td></tr><tr><td>9.7–14.3</td><td>55</td><td>357</td><td>15.4 [11.8–19.6]</td><td></td></tr><tr><td>14.3–19.1</td><td>37</td><td>362</td><td>10.2 [7.3–13.8]</td><td></td></tr><tr><td>&gt; 19.1</td><td>116</td><td>354</td><td>32.8 [27.9–37.9]</td><td></td></tr><tr><td colspan="2">Factor 1</td><td>/</td><td>/</td><td>/</td><td>p &lt; 0.01</td></tr><tr><td colspan="2">Factor 2</td><td>/</td><td>/</td><td>/</td><td>p &gt; 0.20</td></tr><tr><td colspan="2">Factor 3</td><td>/</td><td>/</td><td>/</td><td>p &lt; 0.10</td></tr><tr><td colspan="2">Factor 4</td><td>/</td><td>/</td><td>/</td><td>p &lt; 0.10</td></tr><tr><td colspan="2">Total cattle</td><td></td><td>277</td><td>1432</td><td>15.9 [14.0–17.8]</td><td>/</td></tr><tr><td rowspan="10">Human</td><td rowspan="2">Habitat</td><td>Urban</td><td>9</td><td>150</td><td>6.0 [2.8–11.1]</td><td>p &lt; 0.20</td></tr><tr><td>Rural</td><td>150</td><td>1,529</td><td>9.8 [8.4–11.4]</td><td></td></tr><tr><td rowspan="2">Gender</td><td>F</td><td>50</td><td>851</td><td>5.9 [4.4–7.7]</td><td>p &lt; 0.001</td></tr><tr><td>M</td><td>109</td><td>828</td><td>13.2 [10.9–15.7]</td><td></td></tr><tr><td rowspan="2">Contact with ruminant</td><td>No</td><td>103</td><td>1,209</td><td>8.5 [7.0–10.2]</td><td>p &lt; 0.05</td></tr><tr><td>Yes</td><td>56</td><td>470</td><td>11.9 [9.1–15.2]</td><td></td></tr><tr><td rowspan="2">Contact with raw milk</td><td>No</td><td>140</td><td>1,576</td><td>8.9 [7.5–10.4]</td><td>p &lt; 0.005</td></tr><tr><td>Yes</td><td>19</td><td>103</td><td>18.4 [11.5–27.3]</td><td></td></tr><tr><td rowspan="2">Contact with fresh ruminant fluids</td><td>No</td><td>158</td><td>1,675</td><td>9.4 [8.1–10.9]</td><td>NA</td></tr><tr><td>Yes</td><td>1</td><td>5</td><td>20 [0.1–71.6]</td><td></td></tr><tr><td rowspan="5">Profession</td><td>Butcher</td><td>1</td><td>6</td><td>16.7 [0.0–64.4]</td><td>p &lt; 0.005</td></tr><tr><td>Farmers</td><td>95</td><td>755</td><td>12.6 [10.3–15.2]</td><td></td></tr><tr><td>Health</td><td>1</td><td>19</td><td>5.3 [0.0–26.0]</td><td></td></tr><tr><td>Contact with environment</td><td>9</td><td>52</td><td>17.3 [8.2–30.3]</td><td></td></tr><tr><td>Others</td><td>53</td><td>847</td><td>6.3 [4.7–8.1]</td><td></td></tr><tr><td rowspan="4">Age</td><td>18 to 26</td><td>30</td><td>455</td><td>6.6 [4.5–9.3]</td><td>p &lt; 0.05</td></tr><tr><td>26 to 37</td><td>35</td><td>423</td><td>8.3 [5.8–11.3]</td><td></td></tr><tr><td>37 to 46</td><td>40</td><td>361</td><td>11.1 [8.0–14.8]</td><td></td></tr><tr><td>&gt; 46</td><td>54</td><td>440</td><td>12.3 [9.4–15.7]</td><td></td></tr><tr><td rowspan="4">Cattle density per sq. km</td><td>&lt; 6.3</td><td>51</td><td>450</td><td>11.3 [8.6–14.6]</td><td>p &lt; 0.05</td></tr><tr><td>6.3–11.7</td><td>51</td><td>420</td><td>12.1 [9.2–15.7]</td><td></td></tr><tr><td>11.7–22.0</td><td>28</td><td>389</td><td>7.2 [4.8–10.2]</td><td></td></tr><tr><td>&gt; 22.0</td><td>29</td><td>420</td><td>6.9 [4.7–9.8]</td><td></td></tr><tr><td colspan="2">Factor 1</td><td>/</td><td>/</td><td>/</td><td>p &gt; 0.20</td></tr><tr><td colspan="2">Factor 2</td><td>/</td><td>/</td><td>/</td><td>p &lt; 0.01</td></tr><tr><td colspan="2">Factor 3</td><td>/</td><td>/</td><td>/</td><td>p &lt; 0.2</td></tr><tr><td colspan="2">Factor 4</td><td>/</td><td>/</td><td>/</td><td>p &lt; 0.10</td></tr><tr><td colspan="2">Total human</td><td></td><td>159</td><td>1,679</td><td>9.5 [8.1–11.0]</td><td>/</td></tr></tbody></table>

doi:10.1371/journal.pntd.0004827.t002

from Factor 1, Factor 2 and Factor 3 in both cattle and human multivariate models. The multi-collinearity test did not detect any correlation between human related factors (VIF < 1.5).

For cattle, the single selected model (weight 0.99; S1 Table) included age, cattle density and Factor 4 as explanatory variables (S1 Table and Table 3). Factor 4 and age had a significant positive effect on seroprevalence (estimation of fixed effect at 1.73 and 0.17 respectively; p<0.001
<table><caption>Table 3. Results from the best cattle model.</caption><thead><tr><th>Variable</th><th></th><th>Estimate</th><th>95% CI</th><th>p-value</th></tr></thead><tbody><tr><td>Intercept</td><td>/</td><td>-2.34</td><td>[-3.02–-1.72]</td><td>/</td></tr><tr><td>Age</td><td>/</td><td>0.17</td><td>[0.10–0.23]</td><td>p < 0.001</td></tr><tr><td rowspan="3">Cattle density per sq. km</td><td>&lt; 6.3</td><td>Reference</td><td>/</td><td>/</td></tr><tr><td>6.3–11.7</td><td>-0.24</td><td>[-1.01–-0.54]</td><td>NS</td></tr><tr><td>11.7–22.0</td><td>-0.66</td><td>[-1.61–-0.24]</td><td>NS</td></tr><tr><td rowspan="2">Factor 4</td><td>&gt; 22.0</td><td>0.97</td><td>[0.30–1.69]</td><td>p < 0.01</td></tr><tr><td>/</td><td>1.73</td><td>[0.96–2.55]</td><td>p < 0.001</td></tr></tbody></table>

NS = not significant

doi:10.1371/journal.pntd.0004827.t003

for both explanatory variables; Table 3). Areas with high cattle density (> 19.1 per sq. km) were at risk (p<0.01; OR = 2.6 95% CI [1.3–5.4]; Table 3). According to AIC, seven models were considered as suitable for describing seroprevalence in humans and thus were analyzed using a multi-model inference approach (S1 Table). These models included age, Factor 2, Factor 3, Factor 4, gender, habitat, contact with raw milk, contact with fresh ruminant product, with live ruminant as explanatory variables (S1 Table and Table 4). Age, gender (male; OR = 2.3 95% CI [1.6–3.3]) and Factor 4 had a significant positive effect on seroprevalence (p<0.001, p<0.05 and p<0.05 respectively; Table 4). Factor 2 had a significant negative effect on seroprevalence (p<0.05) whereas Factor 3 had a minor importance in this set of models (RI = 0.27; Table 4). Contact with raw milk had a moderate effect on individual seroprevalence (OR = 1.6 95% CI [1.0–3.5]) whereas direct contacts ruminants and/or with fresh ruminant fluids, and habitat had a low impact on seroprevalence (RI = 0.12 or less; Table 4).

Internal validity of both cattle and human sets of models were satisfactory with an Area Under the Curve (AUC) of 0.82 (95% CI [0.79–0.84]) and 0.80 (95% CI [0.77–0.84]) for cattle and human models respectively. The 10-fold cross-validation estimated an individual prediction error of about 14%.

Cattle seroprevalence was predicted according to Factor 4, cattle density categories and for a fixed cattle age of 5 years. To avoid biased estimations resulting from extrapolations, the prediction of seroprevalence was restricted to communes included in the range of the Factor 4 values corresponding to communes where cattle were sampled (i.e [-1.1–-0.6]; n = 1,368). The

<table><caption>Table 4. Results from the multi-model inference approach for human dataset analysis.</caption><thead><tr><th>Variables</th><th>model-averaged fixed effects (mafe)</th><th>95% CI</th><th>p-value</th><th>Relative importance (RI)</th><th>Number of models</th></tr></thead><tbody><tr><td>Age</td><td>0.02</td><td>[0.01–0.03]</td><td>0.001</td><td>1</td><td>7</td></tr><tr><td>Factor 2</td><td>-0.41</td><td>[-0.74–-0.09]</td><td>0.05</td><td>1</td><td>7</td></tr><tr><td>Factor 3</td><td>0.17</td><td>[-0.08–0.41]</td><td>NS</td><td>0.27</td><td>2</td></tr><tr><td>Factor 4</td><td>0.34</td><td>[0.08–0.61]</td><td>0.05</td><td>1</td><td>7</td></tr><tr><td>Gender</td><td>0.83</td><td>[0.52–1.14]</td><td>0.001</td><td>1</td><td>7</td></tr><tr><td>Contact with raw milk</td><td>0.60</td><td>[0.05–1.15]</td><td>NS</td><td>0.75</td><td>5</td></tr><tr><td>Contact with fresh ruminant fluids</td><td>1.04</td><td>[-1.26–-3.36]</td><td>NS</td><td>0.12</td><td>1</td></tr><tr><td>Cattle density categories</td><td>/</td><td>/</td><td>NS</td><td>/</td><td>0</td></tr><tr><td>Profession</td><td>/</td><td>/</td><td>NS</td><td>/</td><td>0</td></tr><tr><td>Contact with ruminant</td><td>-0.07</td><td>[-0.44–-0.29]</td><td>NS</td><td>0.10</td><td>1</td></tr><tr><td>Habitat</td><td>-0.42</td><td>[-1.42–-0.57]</td><td>NS</td><td>0.12</td><td>1</td></tr></tbody></table>

NS = not significant

doi:10.1371/journal.pntd.0004827.t004
<figure><img src="image_8.png" /><figcaption>**Fig 3.** Predicted cattle seroprevalence in Madagascar and areas affected by RVF outbreaks in ruminant during 1990–1991 and 2008–2009. The cattle seroprevalence (SP) was predicted per commune and according to the best cattle model (Factor 4, cattle density and fixed age 5 years old).</figcaption></figure>

doi:10.1371/journal.pntd.0004827.g003

prediction map highlights the western, north-western part and eastern-coast of Madagascar as high-risk areas for RVF transmission (Fig 3). Nineteen percent of the communes affected by outbreaks in ruminants during the 1990–91 and 2008–09 epizootics are located in areas with a predicted seroprevalence higher than 25%. Yet, 24% of the communes affected by these epizootics are located in low risk areas (predicted seroprevalence lower than 10%). Observed and predicted seroprevalence at the district level are compared in the S1 Appendix.

Discussion

Following the 2008–09 epidemics, studies showed that RVFV spread widely but heterogeneously over Madagascar in both cattle and human populations [15,17]. This could be explained by the presence of ecosystems that are more or less suitable to the RVF candidate vector genera in Madagascar, including mosquitoes in the *Aedes*, *Anopheles*, *Culex*, *Eretmapo-dites* and *Mansonia* genera [25, 28]. Indeed, vector density and population dynamics are influenced by environmental factors such as climate and landscape features [1,25]. However, to date, environmental factors linked to the transmission of RVFV have never been investigated
in Madagascar. To characterize Malagasy environments, we used MFA methods to generate environmental indicators that combined climatic, NDVI and landscape variables selected according to their putative influence on mosquito vector populations. Our MFA classification is in accordance with the known Malagasy ecosystems [18].

The risk of transmission and respective roles of direct and vectorial transmission are probably different among eco-climatic areas. In the case of direct transmission the force of infection is expected to depend on the number of potentially infectious contacts that a susceptible individual experiences over a time unit. This contact rate is expected to depend, among others, on cattle density. A positive association between cattle density and IgG seroprevalence rate in cattle and/or humans would thus suggest a direct transmission of RVFV. In the case of direct cattle to human transmission, the force of infection in the human population is also expected to depend on the frequency of human behaviors resulting in exposure to ruminant fluids or products. By contrast, in the case of vectorial transmission, due to the so-called "dilution effect" and for a fixed vector density, increased cattle density would decrease the probability for a susceptible individual to be bitten by an infectious vector over a time unit [50,51]. Therefore, a negative association between cattle density and IgG seroprevalence rates in cattle and/or in humans would rather suggest a vectorial transmission. The force of infection is also expected to increase with vector density. In Madagascar the density of vectors mostly depends on climatic and landscape factors [1, 25, 27]. High cattle densities are not systematically associated with high vector densities, as the main RVF vectors breed rather in large water bodies [25] than in artificial containers created from livestock-related activities. Thus, a positive association between cattle and/or human seroprevalence, local environmental and climatic conditions favorable to mosquitoes is expected under the hypothesis of vectorial transmission.

According to our analysis, cattle seroprevalence increased with age suggesting an enzootic circulation. Cattle seropositivity was positively associated to humid environment (large surface of permanent wetlands, marshlands and irrigated lands) each of these factors being favorable to *Culex* and *Anopheles* mosquitoes [25]. Actually, during the 2008–09 epidemics, 3 mosquito species were found to be naturally infected by RVFV: *Anopheles coustani*, *An. squamosus* and *Culex antennatus* [52]. *Cx. antennatus* is considered a RVFV vector and both *Anopheles* species as candidate vectors [25]. Cattle seroprevalence was also positively associated with cattle density suggesting the existence of a direct transmission between cattle, as suggested by Nicolas et al [3,19]. However, in our study cattle density and environmental factors were not independent (correlation with Factor 1, Factor 2 and Factor 3). Because of such associations it was impossible to disentangle the influence of cattle density from the influence of environmental conditions and thus to thoroughly assess the relative importance of vectorial and direct transmission. The prediction map of cattle seroprevalence highlighted the eastern-coast, western and north-western parts as high-risk areas. Surprisingly, some districts affected by RVFV outbreaks are located in the predicted low risk area [13–15]. The last outbreaks were mostly reported in the highlands, which are highly connected by road to the capital city, Antananarivo. Outbreaks occurring in isolated areas may have been missed explaining why a low proportion of outbreaks were located in predicted at-risk areas. On the other hand, enzootic transmission could have maintained a sufficient level of immunity in cattle in the high risk area restraining the outbreak magnitude in these regions. RVF could have been introduced in low risk areas through cattle trade and because of the low level of cattle immunity in these zones, trigger outbreaks. Nevertheless, as RVF cases were suspected to be under-reported, it was not possible to assess the relationship between the prediction of the herd immunity and the case notifications. Using satellite measurements (sea surface temperatures, rainfall and NDVI) and human cases as model output, Anyamba et al. [29] identified mainly the east-coast and some small areas of northern and north-western parts as at-risk for 2008–2009 RVF outbreaks in Madagascar.
Considering that, during the 2008–09 outbreaks, several human cases occurred from the con-
tact with infected fresh meat from traded ruminants [15], all the human infections could not
be attributed to local infection [15]. Moreover, the detection of human cases depends on the
intensity of the local circulation between ruminants and vectors, the probability of human
exposure, the presence of clinical signs and the declaration to health services. Then, the human
case data were probably not an optimum indicator of spatial distribution of RVF cases as sug-
gested by Anyamba et al. [29]. Our prediction map is based on cattle for which the infection
could be attributed to a local infection and identifies larger at-risk areas on western part of
Madagascar than Anyamba et al. [29]. The discrepancy between results of Anyamba et al.
study [29] and our study may be due the methodological differences: environmental variables
included in both models and human clinical cases as model output on one side, bovine serolog-
ical results on the other hand.

The estimated overall human seroprevalence was 9.5% (IC95% [8.2–11.0]). This seropreva-
lence is higher than adult seroprevalence observed in the island of Mayotte (2011) and Tanza-
nia (2007–08) [24,53] but lower than adult seroprevalence in Kenya or Saudi Arabia [54,55].
Additionally, this seroprevalence is higher than the seroprevalence estimated for Madagascar
in Gray et al [21]. The difference in the sampling area could explain this difference. Indeed,
sera from the study of Gray et al. [21] were mainly sampled in south where RVF seroprevalence
in human is low. Because of the different eco-epidemiological contexts and survey settings it is
difficult to compare our results with the studies performed in Mayotte, Tanzanian, Kenya and
Saudi Arabia [24, 53–55]. Human RVF seropositivity increased with age, suggesting an
endemic transmission in human populations. As observed in cattle, human seropositivity was
positively associated with the presence of temporary and artificial water points. In addition, 24
seropositive individuals declared no contact with ruminant or ruminant products, and the 3
mosquito species considered as potential vectors in Madagascar are zoo-anthropophilic feeders
[25,52]: these results strongly suggest the existence of a vectorial transmission from ruminant
to humans. Our analysis showed that frequent contact with raw milk contributed to explain
human infection as previously suspected in Kenya [31]. Direct contact with fresh blood was
not identified as human related risk factor whereas this way is suspected to be the main route
of human infection in other studies [31]. In our sample, the number of people in contact with
fresh blood was very low resulting in a low statistical power. However, this way of transmission
has still to be considered, especially in the areas unfavorable to mosquitoes where direct contact
could explain human infections [15].

Our integrated approach analyzing environmental, cattle and human datasets allow us to
bring new insight on RVF transmission patterns in Madagascar. The association between cattle
seroprevalence, humid environments and high cattle density suggests that concomitant vecto-
rial and direct transmissions are critical to maintain RVFV enzootic transmission.

Even if the 2008–09 outbreaks are suspected to be associated with infected domestic animals
imported from east Africa [56], our study confirms that enzootic and endemic circulations
occur in Madagascar as suggested before [3,12,21].

The identification of at-risk environments is essential to focus veterinary surveillance and
control of RVFV. Because of the variety of ecosystems and socio-cultural practices in Madagas-
car, it is likely that some areas are more favorable to direct transmission [3,19], while others are
more favorable to vectorial transmission or to both transmission pathways. In the at-risk humid
environment of the western, north-western and the eastern-coast areas, suitable for *Culex* and
*Anopheles* mosquitoes, vectorial transmission probably occur in both cattle and human. In the
future, mathematical modeling may be used to decipher the relative contribution of each trans-
mission pathway in both human and ruminants, integrate the role of animal trade in disease
spread in the Malagasy context, and thus propose adapted surveillance and control measures.
**Supporting Information**

S1 Table. Comparison of the values and weight of AIC for the cattle and human models.
(DOCX)

**S1 Appendix. Scatterplot of observed versus predicted seroprevalences at the district level.**
Seroprevalence has been predicted for each age category in each communes sampled. For each district the sampling has been reconstructed taking into account the communes sampled and the number of animals sampled in each commune. Grey points correspond to districts where less than 5 animals were sampled.
(DOCX)

**Acknowledgments**

We especially thank the population of Madagascar who participated to the studies. We thank those who facilitated the survey, i.e., heads of fokontany, local administration authorities and health authorities from Ministry of Health. We also thank the Plague Unit at the Institut Pasteur de Madagascar for data collection and supporting (S. Telfer, C. Rahaingosoamamitiana, F. M. Andriamiarimanana, S. Rahelinirina, M. Rajerison), S. Andrimasinoro for the management of data, B.S. Rahoilijaona H.A. Rakotoarison, H. Raharimampianina and A.M. Rakotohaingomahefa for their field supports. We are grateful to the authors of the cattle survey and especially E. Jeanmaire, J.M. Reynes and S. de la Rocque for providing the data of cattle survey. We thank G. Gray from the Division of Infectious Diseases of Duke University for its support. Finally, we thank three anonymous reviewers for their careful reading of our manuscript and their comments and suggestions.

**Author Contributions**

Conceived and designed the experiments: MMO VC VG AT SFA BD FR CR JMH. Performed the experiments: MMO SFA JPR SA FR. Analyzed the data: MMO VC VG AT. Wrote the paper: MMO VC VG AT BD JMH.

**References**

1. Chevalier V. Relevance of Rift Valley fever to public health in the European Union. Clin Microbiol Infect. 2013 August; 19(8):705-708. doi: 10.1111/1469-0691.12163 PMID: 23517372

2. Pepin M, Bouloy M, Bird BH, Kemp A, Paweska J. Rift Valley fever virus(Bunyaviridae: Phlebovirus): an update on pathogenesis, molecular epidemiology, vectors, diagnostics and prevention. Vet Res. 2010 December; 41(6):61. doi: 10.1051/vetres/2010033 PMID: 21188836

3. Nicolas G, Durand B, Rakotoarimanana TT, Lacote S, Chevalier V, Marianneau P. A 3-year serological and virological cattle follow-up in Madagascar highlands suggests a non-classical transmission route of Rift Valley fever virus. Am J Trop Med Hyg. 2014 February; 90(2):265-266. doi: 10.4269/ajtmh.13-0538 PMID: 24366500

4. Nakounné E, Selekon B, Morvan J. [Microbiological surveillance: viral hemorrhagic fever in Central African Republic: current serological data in man]. Bull Societe Pathol Exot. 2000 January; 93(5):340-347.

5. Pourrut X, Nkoghé D, Souris M, Paupy C, Paweska J, Padilla C, et al. Rift Valley fever virus seroprevalence in human rural populations of Gabon. PLoS Negl Trop Dis. 2010 July 27; 4(7):e763. doi: 10.1371/journal.pntd.0000763 PMID: 20668541

6. Jupp PG, Kemp A, Grobbelaar A, Lema P, Burt FJ, Alahmed AM, et al. The 2000 epidemic of Rift Valley fever in Saudi Arabia: mosquito vector studies. Med Vet Entomol. 2002 September; 16(3):245-252. doi: 10.1046/j.1365-2915.2002.00371.x PMID: 12243225

7. Hightower A, Kinkade C, Nguku PM, Anyangu A, Mutonga D, Omolo J, et al. Relationship of climate, geography, and geology to the incidence of Rift Valley fever in Kenya during the 2006–2007 outbreak. Am J Trop Med Hyg. 2012 February; 86(2):373-380. doi: 10.4269/ajtmh.2012.11-0450 PMID: 22302875
8. Sindato C, Karimuribo ED, Pfeiffer DU, Mboera LEG, Kivaria F, Dautu G, et al. Spatial and temporal pattern of Rift Valley fever outbreaks in Tanzania; 1930 to 2007. PLoS One. 2014 February 25; 9(2): e88897. doi: 10.1371/journal.pone.0088897 PMID: 24586433

9. Meegan JM. The Rift Valley fever epizootic in Egypt 1977–78. 1. Description of the epizzotic and virological studies. Trans R Soc Trop Med Hyg. 1979; 73(6):618-623. doi: 10.1016/0035-9203(79)90004-X

10. Digoutte JP, Peters CJ. General aspects of the 1987 Rift Valley fever epidemic in Mauritania. Res Virol. 1989 February; 140(1):27-30. doi: 10.1016/S0923-2516(89)80081-0 PMID: 2711042

11. Hassan OA, Ahlm C, Sang R, Evander M. The 2007 Rift Valley fever outbreak in Sudan. PLoS Negl Trop Dis. 2011 September 27; 5(9):e1229. doi: 10.1371/journal.pntd.0001229 PMID: 21980543

12. Chevalier V, Rakotondrafara T, Jourdan M, Heraud JM, Andriamanivo HR, Durand B, et al. An unexpected recurrent transmission of Rift Valley fever virus in cattle in a temperate and mountainous area of Madagascar. PLoS Negl Trop Dis. 2011 December 20; 5(12):e1423. doi: 10.1371/journal.pntd.0001423 PMID: 22206026

13. Morvan J, Saluzzo JF, Fontenille D, Rollin PE, Coulanges P. Rift Valley fever on the east coast of Madagascar. Res Virol. 1991 December; 142(6):475-482. doi: 10.1016/0923-2516(91)90070-J PMID: 1687082

14. Morvan J, Rollin PE, Laventure S, Rakotoarivony I, Roux J. Rift Valley fever epizootic in the central highlands of Madagascar. Res Virol. 1992 December; 143(6):407-415. doi: 10.1016/S0923-2516(06) 80134-2 PMID: 1297176

15. Andriamandimby SF, Randrianarivo-Solofoniaina AE, Jeanmaire EM, Ravololomanana L, Razafimantsoa LT, Rakotojoelinandrasana T, et al. Rift Valley fever during rainy seasons, Madagascar, 2008 and 2009. Emerg Infect Dis. 2010 June; 16(6):963-970. doi: 10.3201/eid1606.091266 PMID: 20507747

16. Jeanmaire E, de La Rocque S. Rift Valley fever in Madagascar: an updated map of the distribution of the disease in 2008. Rome: Food and Agriculture Organization of the United Nations. EMPRESS. 2010; 35:23-30. Available: http://www.fao.org/3/a-i1648e/i1648e05.pdf. Accessed 22 March 2016.

17. Jeanmaire EM, Rabenarivahiny R, Biarmann M, Rabibisoa L, Ravaomanana F, Randriamparany T, et al. Prevalence of Rift Valley fever infection in ruminants in Madagascar after the 2008 outbreak. Vector Borne Zoonotic Dis. 2011 April 12; 11(4):395-402. doi: 10.1089/vbz.2009.0249 PMID: 21395414

18. Cornet A. Essai de cartographie bioclimatique à Madagascar. (Notice Explicative; 55). Paris: ORS-TOM. 1974; Available: http://www/documentation.ird.fr/hor/fdi:06946. Accessed 10 October 2015.

19. Nicolas G, Chevalier V, Tantely LM, Fontenille D, Durand B. A spatially explicit metapopulation model and cattle trade analysis suggests key determinants for the recurrent circulation of rift valley Fever virus in a pilot area of Madagascar highlands. PLoS Negl Trop Dis. 2014 December 4; 8(12):e3346. doi: 10.1371/journal.pntd.0003346 PMID: 25474116

20. Goutard FL, Binot A, Duboz R, Rasamoelina-Andriamanivo H, Pedrono M, Holl D, et al. How to reach the poor? Surveillance in low-income countries, lessons from experiences in Cambodia and Madagascar. Prev Vet Med. 2015 June 1; 120(1):12-26. doi: 10.1016/j.prevetmed.2015.02.014 PMID: 25842000

21. Gray GC, Anderson BD, LaBeaud AD, Heraud J-M, Fèvre EM, Andriamandimby SF, et al. Seroepidemiological Study of Interepidemic Rift Valley Fever Virus Infection Among Persons with Intense Ruminant Exposure in Madagascar and Kenya. Am J Trop Med Hyg. 2015 December 9; 93(6):1364-1370. doi: 10.4269/ajtmh.15-0383 PMID: 26458775

22. Paweska JT, Burt FJ, Anthony F, Smith SJ, Grobbelaar AA, Croft JE, et al. IgG-sandwich and IgM-capture enzyme-linked immunosorbent assay for the detection of antibody to Rift Valley fever virus in domestic ruminants. J Virol Methods. 2003 November; 113(2):103-112. doi: 10.1016/S0166-0934(03) 00228-3 PMID: 14553896

23. Paweska JT, Burt FJ, Swanepoel R. Validation of IgG-sandwich and IgM-capture ELISA for the detection of antibody to Rift Valley fever virus in humans. J Virol Methods. 2005 March; 124(1-2):173-181. doi: 10.1016/j.jviromet.2004.11.020 PMID: 15664066

24. Heinrich N, Saathoff E, Weller N, Clowes P, Kroidl I, Ntinginya E, et al. High seroprevalence of Rift Valley fever and evidence for endemic circulation in Mbeya region, Tanzania, in a cross-sectional study. PLoS Negl Trop Dis. 2012 March; 6(3):e1557. doi: 10.1371/journal.pntd.0001557 PMID: 22479657

25. Tantely LM, Boyer S, Fontenille D. A review of mosquitoes associated with Rift Valley fever virus in Madagascar. Am J Trop Med Hyg. 2015 April; 92(4):722-729. doi: 10.4269/ajtmh.14-0421 PMID: 25732680

26. Muiruri S, Kabiru EW, Muchiri EM, Hussein H, Kagondu F, LaBeaud AD, et al. Cross-sectional survey of Rift Valley fever virus exposure in Bodhei village located in a transitional coastal forest habitat in Lamu county, Kenya. Am J Trop Med Hyg. 2015 February; 92(2):394-400. doi: 10.4269/ajtmh.14-0440 PMID: 25535309
27. Owange NO, Ogara WO, Affognon H, Peter GB, Kasiiti J, Okuthe S, et al. Occurrence of rift valley fever in cattle in Ijara district, Kenya. Prev Vet Med. 2014 November; 117(1):121–8. doi: 10.1016/j. prevetmed.2014.08.008 PMID: 25217406

28. Tantely ML, Le Goff G, Boyer S, Fontenille D. An updated checklist of mosquito species (Diptera: Culicidae) from Madagascar. Parasite 2016 April 21; 23:20. doi: 10.1051/parasite/2016018 PMID: 27101839

29. Anyamba A, Linthicum KJ, Small J, Britch SC, Pak E, de La Rocque S, et al. Prediction, assessment of the Rift Valley fever activity in East and Southern Africa 2006–2008 and possible vector control strategies. Am J Trop Med Hyg. 2010 August; 83(2 Suppl):43-51. doi: 10.4269/ajtmh.2010.09-0289 PMID: 20682905

30. Metras R, Jewell C, Porphyre T, Thompson PN, Pfeiffer DU, Collins LM, et al. Risk factors associated with Rift Valley fever epidemics in South Africa in 2008–11. Sci Rep. 2015 March 25; 5:9492. doi: 10.1038/srep09492 PMID: 25804974

31. Nicholas DE, Jacobsen KH, Waters NM. Risk factors associated with human Rift Valley fever infection: systematic review and meta-analysis. Trop Med Int Health. 2014 December; 19(12):1420-1429. doi: 10.1111/tmi.12385 PMID: 25252137

32. Robinson TP, Wint GRW, Conchedda G, Van Boeckel TP, Ercoli V, Palamara E, et al. Mapping the global distribution of livestock. PloS One. 2014 May 29; 9(5):e96084. doi: 10.1371/journal.pone.0096084 PMID: 24875496

33. Bontemps S, Defourny P, Van Bogaert E, Arino O, Kalogirou V, Perez JR. GLOBCOVER 2009: Products Description and Validation Report. UCLouvain & ESA Team; 2011. Available: http://epic.awi.de/31014/16/GLOBCOVER2009_Validation_Report_2-2.pdf. Accessed 10 June 2015.

34. Moat J, Du Puy D. Madagascar, Remaining Primary Vegetation [Internet]. Royal Botanic Gardens, Kew (Eds.), London; 1997. Available: www.kew.org/gis/projects/madagascar/veg_mapping.html. Accessed 14 November 2014.

35. Batjes NH. IPCC default soil classes derived from the Harmonized World Soil Data Base (Ver. 1.1). Carbon Benefits Project (CBP) and ISRIC–World Soil Information; 2010. Available: http://www.isric.org/isric/webdocs/docs/ISRIC_Report_2009_02.pdf. Accessed 12 November 2014.

36. Siebert S, Frenken K. Irrigated areas. In: Atlas of African agriculture research & development; 2014. Available: http://www.fao.org/nr/water/aquastat/irrigationmap/index.stm. Accessed 12 November 2014.

37. QGIS Development Team. QGIS Geographic Information System [Internet]. Open Source Geospatial Foundation; 2009. Available: http://qgis.osgeo.org. Accessed 10 June 2014.

38. Escofier B, Pagès J. Multiple factor analysis (AFMULT package). Comput Stat Data Anal. 1994 August; 18(1):121-140. doi: 10.1016/0167-9473(94)90135-X

39. Costard S, Porphyre V, Messad S, Rakotondrahanta S, Vidon H, Roger F, et al. Multivariate analysis of management and biosecurity practices in smallholder pig farms in Madagascar. Prev Vet Med. 2009 November 15; 92(3):199-209. doi: 10.1016/j.prevetmed.2009.08.010 PMID: 19781801

40. Grueber CE, Nakagawa S, Laws RJ, Jamieson IG. Multimodel inference in ecology and evolution: challenges and solutions. J Evol Biol. 2011 April; 24(4):699-711. doi: 10.1111/j.1420-9101.2010.02210.x PMID: 21272107

41. Park SH, Goo JM, Jo C-H. Receiver operating characteristic (ROC) curve: practical review for radiologists. Korean J Radiol. 2004 March; 5(1):11–18. doi: 10.3348/kjr.2004.5.1.11 PMID: 15064554

42. Franco AO, Davies CR, Mylne A, Dedet J-P, GáLlego M, Ballart C, et al. Predicting the distribution of canine leishmaniasis in western Europe based on environmental variables. Parasitology. 2011 December; 138(14):1878–91. doi: 10.1017/S003118201100148X PMID: 21914251

43. R Development Core Team. R: A Language and Environment for Statistical Computing. Vienna, Austria: the R Foundation for Statistical Computing. 2011.

44. Lê S, Josse J, Husson F. FactoMineR: An R Package for Multivariate Analysis. J Stat Softw. 2008 March 18; 25(1):1. doi: 10.18637/jss.v025.i01

45. Fox J, Weisberg S. An {R} Companion to Applied Regression. Second. Sage; 2011. Available: http://socserv.socsci.mcmaster.ca/jfox/Books/Companion. Accessed 10 October 2015.

46. Bates D, Maechler M, Bolker B, Walker S. Fitting Linear Mixed-Effects Models using lme4. Journal of Statistical Software. 2015;51 pages. Available: http://arxiv.org/abs/1406.5823. Accessed 10 October 2015.

47. Bartoń K. MuMln: Multi-Model Inference; 2015. Available: https://cran.r-project.org/web/packages/MuMln/index.html. Accessed 26 August 2015.

48. Robin X, Turck N, Hainard A, Tiberti N, Lisacek F, Sanchez J-C, et al. pROC: an open-source package for R and S+ to analyze and compare ROC curves. BMC Bioinformatics. 2011 March 17; 12(1):77. doi: 10.1186/1471-2105-12-77
49. Canty A, Ripley B. Package "boot": Bootstrap Functions. Cran.2016. Available: https://cran.r-project.org/web/packages/boot/boot.pdf. Accessed 24 March 2016.

50. Keeling MJ, Rohani P. Modeling infectious diseases in humans and animals. Princeton: Princeton University Press; 2008.

51. Durand B, Zanella G, Biteau-Coroller F, Locatelli C, Baurier F, Simon C, et al. Anatomy of bluetongue virus serotype 8 epizootic wave, France, 2007–2008. Emerg Infect Dis. 2010 December; 16 (12):1861-1868. doi: 10.3201/eid1612.100412 PMID: 21122214

52. Ratovonjato J, Olive M-M, Tantely LM, Andrianavolambo L, Tata E, Razainirina J, et al. Detection, isolation, and genetic characterization of Rift Valley fever virus from Anopheles (Anopheles) coustani, Anopheles (Anopheles) squamosus, and Culex (Culex) antennatus of the Haute Matsiatra region, Madagascar. Vector Borne Zoonotic Dis. 2011 June 16; 11(6):753-759. doi: 10.1089/vbz.2010.0031 PMID: 21028960

53. Lernout T, Cardinale E, Jego M, Desprès P, Collet L, Zumbo B, et al. Rift valley fever in humans and animals in Mayotte, an endemic situation? PloS One. September 2013; 8(9):e74192. doi: 10.1371/journal.pone.0074192 PMID: 24098637

54. LaBeaud AD, Ochiai Y, Peters CJ, Muchiri EM, King CH. Spectrum of Rift Valley fever virus transmission in Kenya: insights from three distinct regions. Am J Trop Med Hyg. 2007 May; 76(5):795-800. PMID: 17488893

55. Memish ZA, Masri MA, Anderson BD, Heil GL, Merrill HR, Khan SU, et al. Elevated Antibodies Against Rift Valley Fever Virus Among Humans with Exposure to Ruminants in Saudi Arabia. Am J Trop Med Hyg. 2015 January4; 92(4):739-43. doi: 10.4269/ajtmh.14-0575 PMID: 25646253

56. Carroll SA, Reynes J-M, Khristova ML, Andriamandimby SF, Rollin PE, Nichol ST. Genetic evidence for Rift Valley fever outbreaks in Madagascar resulting from virus introductions from the East African mainland rather than enzootic maintenance. J Virol. 2011 July; 85(13):6162-6167. doi: 10.1128/JVI.00335-11 PMID: 21507967