Climate, Cattle Rearing Systems and African Animal Trypanosomosis Risk in Burkina Faso

Background In sub-Saharan countries infested by tsetse flies, African Animal Trypanosomosis (AAT) is considered as the main pathological constraint to cattle breeding. Africa has known a strong climatic change and its population was multiplied by four during the last half-century. The aim of this study was to characterize the impact of production practices and climate on tsetse occurrence and abundance, and the associated prevalence of AAT in Burkina Faso. Methodology/Principal Findings Four sites were selected along a South-north transect of increasing aridity. The study combines parasitological and entomological surveys. For the parasitological aspect, blood samples were collected from 1,041 cattle selected through a stratified sampling procedure including location and livestock management system (long transhumance, short transhumance, sedentary). Parasitological and serological prevalence specific to livestock management systems show a gradual increase from the Sahelian to the Sudano-Guinean area (P<0.05). Livestock management system had also a significant impact on parasitological prevalence (P<0.05). Tsetse diversity, apparent densities and their infection rates overall decreased with aridity, from four species, an apparent density of 53.1 flies/trap/day and an infection rate of 13.7% to an absence at the northern edge of the transect, where the density and diversity of other biting flies were on the contrary highest (p<0.001). Conclusions/Significance The climatic pressure clearly had a negative impact on tsetse abundance and AAT risk. However, the persistency of tsetse habitats along the Mouhoun river loop maintains a high risk of cyclical transmission of T. vivax. Moreover, an “epidemic mechanical livestock trypanosomosis” cycle is likely to occur in the northern site, where trypanosomes are brought in by cattle transhuming from the tsetse infested area and are locally transmitted by mechanical vectors. In Burkina Faso, the impact of tsetse thus extends to a buffer area around their distribution belt, corresponding to the herd transhumance radius.


Introduction
The climatic change and its impact on environments vary extensively in space and time [1], due to local socio-cultural and biophysical conditions. However, the increase of the frequency and amplitude of extreme climatic events (extreme temperatures, droughts, flooding, etc.) is generally fast, like the extreme vulnerability of West African small producers to cope with these changes. These events are not caused by natural variations of the climate only: they are exemplified by an increased exploitation of natural resources, particularly land use [2]. For example, from some 30 million inhabitants in 1900, West African population increased to 306 million in 2010, and could reach between 550 and 700 million in 2050 [3]. The development of extensive agriculture and the increase in human density are associated to landscape fragmentation in the Mouhoun river basin [4].
Associated to extreme climatic events, they have huge consequences on the environment and sustainable development, particularly for African producers who depend strongly on ecosystem services like natural grazing [5]. These changes have also an impact on vectors' environment, and thus on their geographical distribution and density. This situation causes a modification of the hosts-vector contacts and hence, of the epidemiology of vector diseases, particularly African Animal Trypanosomoses (AAT) [6].
Indeed, these diseases represent the main constraint to the development of more intensive livestock production systems in sub-Saharan Africa [7], notably in areas with high agricultural potential. Cattle breeding has an important place in the economy and socio-cultural activities of populations. Thus, controlling trypanosomoses and their vectors becomes a major step to increase the productivity of livestock production systems and contribute to small producers' food security. Several methods are proposed for the control and the reduction of trypanosomosis negative impacts on animal production. However, the impact of the present instability of tsetse flies natural habitats caused by climatic changes and anthropic pressure must be taken into consideration [6]. The global objective of this survey was to analyze the relative impact of production practices, particularly transhumance, and climatic conditions on the epidemiology of AAT and their main vectors (tsetse) in Burkina Faso.

Study area
The survey was achieved in four sites along a South-North transect in Burkina Faso (Fig. 1). These sites were selected according to cattle production systems and the agro-ecological area to which they belong: Folonzo, in the Sudano-Guinean area (16u609W and 9u879N), which is a welcome area for transhumant herds; Koumbia, in the Sudanese area (15u509W and 11u079N), which is a welcome, transit and attachment area for transhumant herds; Dédougou, in the Sudano-Sahelian area (15u559W and 12u499N), a welcome, transit and attachment area for transhumant herds; Djibo, in the Sahelian area (13u609W and 14u109N), an attachment area for transhumant herds.

Parasitological surveys
Blood samples were taken at the jugular vein from a total of 1,041 cattle (Table 1). In each livestock rearing system, cattle were selected randomly without criteria of age, sex or breed. However, the sex and the age of the animals, the date of the last treatment against AAT and the type of trypanocid used were recorded. The sample was composed of Fulani zebus and half-bred between Fulani zebus and ''Baoulé'' trypanotolerant cattle and in majority by females (73.63%). Animals belonged mainly to Fulani cattle farmers.
Two diagnosis methods were used for blood examination of the samples: the Buffy-coat (BCM) [8] to detect active infections by trypanosomes and the indirect ELISA-method to detect antitrypanosomes antibodies in plasma, attesting past or present infections of the animal with trypanosomes [9]. The BCM allowed  the diagnosis of trypanosome species based on morphological, mobility and size criteria [8].

Entomological surveys
All biting flies were trapped, immediately followed by tsetse flies dissection. In each of the 4 sites ( Fig. 1), ten standard biconical traps [10] were set at intervals of 100 meters along cattle watering points. The trapping lasted 72 hours with harvests every 24 hours. Tsetse and other biting flies were numbered by species, by trap, and ADT (apparent density per trap per day) was calculated.
Non teneral tsetse flies were dissected using a binocular microscope. The dissection started with the proboscis, then the salivary glands and finally the mid-gut. After dissection, these organs were placed between a slide and a cover slip in a drop of Ringer's solution, and then observed using a microscope (640) for parasite detection. Flies with an infected organ were collected in eppendorf tubes of 0.5 ml containing 30 ml of sterile distilled water (each organ independently) and stored at 220uC for ulterior PCR (Polymerase Chain Reaction) analysis, to determine the species of trypanosome involved, with primers for T. vivax and T. congolense savannah type [11]. In the same way, females' physiological age was determined by the dissection of the reproductive system [12].

Statistical analyses
Rates of infection of cattle and tsetse in different sites were compared using the Khi square test and binomial mixed effects models, where the cattle rearing system and the site represented the stationary effects. The herd was considered as a random effect.
The ADP were compared overall using a Kruskal-Wallis rank sum test [13] and then by pairs using the Steel-type none parametric multiple comparisons test (npmc) [14]. Tests were performed using the R-2.9.2-win32 software.
The mean number of infectious flies by trap and by day was calculated as an indicator of cyclical transmission risk [15]. The relative risks and their confidence intervals were obtained from bootstrapping in the ADT and infectious rate distributions from each site, assuming spatial homogeneity within a given site (5 000 Monte Carlo simulations, @risk software).

Parasitological surveys
Parasitological and serological prevalences were different between sites (p,0.05). They gradually increased along the climatic gradient of the Sahelian area toward the Sudano-Guinean area (Table 2), from 0 to 20% and 18 to 83% for parasitological and serological prevalence respectively (p,10 23 ).
Parasitological investigations allowed detecting active infections by Trypanosoma vivax and T. congolense with predominance of T. vivax in Dédougou and Koumbia and of T. congolense in Folonzo. No case of active infection by T. brucei brucei was identified in any of the sites (Table 2). In Djibo, no case of active infection by any trypanosome species was found. The cattle rearing system had a significant impact on the parasitological prevalence (p = 0.02) but not on the serological prevalence (p = 0.77).
In Dédougou and Koumbia, parasitological prevalences were significantly higher in the sedentary rearing system that in those practicing long transhumance toward the subhumid and humid areas of the south (p,0.05). In Djibo, animals carriers of antitrypanosomes antibodies were identified in the three cattle rearing systems (Table 2). T. congolense was predominant in the transhuming herds in all the sites (p,10 23 ).

Entomological surveys
Tsetse flies densities varied significantly along the climatic gradient (p,10 23 ). No tsetse was captured in Djibo, contrary to Folonzo where the ADT was 53.10670. 15; 8.1665.55 in Koumbia and 27.24623.72 in Dédougou (Fig. 2). The difference between Koumbia and Dédougou was not significant (p = 0.09) whereas all the other pair-comparisons were highly significant (p,0.02).
The specific density analysis shows that Glossina tachinoides has the largest spectrum of distribution. The highest ADT for this species, 39.63656.69, found in Folonzo, was significantly higher than in Dédougou (21.62620.18, p,0.05) and Koumbia sites (7.0365.34, p,0.001). The ADT of G. tachinoides in Dédougou was also significantly higher than in Koumbia (p,0.001). G. palpalis gambiensis was present in 3 sites, but with lower ADT (1.8762.71 in Koumbia, 0.3660.76 in Dédougou and 0.7961.28 in Folonzo) comparing to those of G. tachinoides. G. morsitans submorsitans and G. medicorum were found only in Folonzo with ADT of 10.5615.23 and 0.661.30 respectively. The tsetse flies diversity decreased with the aridity gradient, from 4 to 0. At the opposite, the diversity of other biting flies increased from 1 to 4 with the aridity gradient (Table 3), and the apparent densities were higher (p,0.001) in the drier site (Djibo) than in the three other sites, where they were similar (p.0.05).
Overall, the risk of cyclical transmission, measured as the apparent density of infectious fly per trap per day, was 1 29] times more important in Folonzo than it the two other sites (p,0.05), between which it was not significantly different (p.0.05).
The physiological age of 59 tsetse in Dédougou, 23 in Koumbia and 43 in Folonzo were measured. The mean age was significantly lower in the site of Koumbia (28618 days) than the two others (p = 0.007). There was no difference between the mean age of the flies in Dédougou (39620 days) and Folonzo (38613 days) (p = 0.74).

Discussion
The importance of AAT prevalence decreased with the aridity degree in our study area, with the exception of Dédougou where the Mouhoun river loop allows the persistency of high tsetse densities [15,16]. The fact that some sites are more infected than others demonstrates the spatial heterogeneity of the AAT along of climatic gradient.
The predominance of T. vivax over T. congolense in Koumbia and in Dédougou could be related to a better transmission by riverine tsetse [17]. On the contrary, the high prevalence of T. congolense in Folonzo reveals the importance of contacts between animals and G. morsitans submorsitans, known as efficient vectors for this trypanosome [18,19]. This hypothesis is reinforced by the fact that, at the time of the entomological survey, G. morsitans submorsitans was found abundant in the area (10.5615.23 tsetse/ trap/day). The absence of T. brucei brucei in cattle and its weak seroprevalence, confirms tendencies observed in recent parasitological studies conducted in Burkina [20,21]. In the study of Dayo et al. (2009), T. congolense was also found much predominant in cattle, in a site close to Folonzo.
Three typical cycles have been described to characterize the epidemiological settings of African Animal Trypanosomosis, including sylvatic trypanosomosis, interface trypanosomoses and endemic livestock trypanosomosis. Koumbia and Dédougou can be considered as within the endemic livestock trypanosomosis cycle, where cattle are the main hosts and wild fauna tends to disappear. On the contrary, Folonzo represents a typical interface trypanosomoses situation where the wild fauna is still abundant and the pressure of cyclical vectors very high. It would be interesting to investigate whether the strains of T. vivax transmitted in the letter cycle are more virulent than in the two other sites, as observed for T. congolense in Southern Africa [22]. The higher parasitological prevalence in sedentary animals in Dédougou and  Koumbia could be explained by higher trypanocid treatment frequencies in transhumant cattle. Even if this was not evidenced from interviews with livestock owners, previous studies have found that such information can be very unreliable and often many more treatments can be given than would be apparent. The resident cattle are concentrated along the last water points (Mouhoun and Bougouriba respectively), and experience high contact intensity with riverine tsetse that are good vectors of T. vivax. This transmission might be seasonally relayed by mechanical vectors, particularly abundant in late rainy season and cold dry season in our study area (September-November) [23]. Similar observations were done in Gambia using parasitological follow-up of transhumant and sedentary animals [24]. The absence of parasitological prevalence in all cattle rearing systems in Djibo might be explained by the combined effect of the absence of cyclical vectors and the frequency of trypanocid treatments by the breeders (before, during and after the return from transhumance in the trypanosomian endemic area). This site is located far upon the northern limit of tsetse, which was very well characterized during former surveys in the study area [4,16,25]. The majority of antibodies detected in animals in Djibo probably result from past infections. This hypothesis is supported by the absence of significant difference between the hematocrit of seropositive and seronegative animals. Sedentary and low transhumance herds are probably contaminated locally by the mechanical vectors from animals infected during their transhumance in the endemic AAT area. Indeed, the letter could maintain an intermittent epidemic transmission of AAT [26], explaining the serological predominance of T. vivax in this site. All these animals (sedentary animals infected mechanically and animals infected during transhumance) are treated several times during the year and were exempt of parasites during the survey, but remain seropositive for a long time (,3 months) [9]. This migratory situation, that always existed, could be exemplified by the climatic accidents and anthropic pressure effects that increase herd movements [6]. The sedentary herds in Djibo cannot be considered as within the three epidemiological cycles formerly described [6] and we thus hereby make the hypothesis of a fourth cycle named ''epidemic mechanical livestock trypanosomosis'' where local transmission is seasonal and ensured by mechanical vectors, where cattle are the main hosts, and where trypanosomes are brought in by cattle coming back from transhumance in the tsetse infested area.
With a mean apparent density of infectious tsetse of 0.74 [95% CI 0.04-2.01], the cyclical transmission risk can be considered as high in Dédougou and Koumbia, according to former studies [15,27], and very high in Folonzo 7.27 [3.66-12.59]. This observation, as well as the growing diversity of tsetse species from the Sahelian to the Sudano-Guinean area is probably related to climatic conditions (rainfall and temperature) [6], the fragmentation degree of vegetal formations caused by the demographic pressure [4] and the presence of wild fauna in Folonzo.
Nowadays, tsetse flies have regressed in the sahelian area where they have been captured until 1935 [25]. The presence of permanent hydrographic network, as the Mouhoun river loop in Dédougou [15], or the Niayes in Senegal, however allow riverine tsetse to persist in high densities in some sahelian sites [16,28], as observed in this survey in Dédougou. Tsetse flies absence in Djibo is attributable to the modification of the hygrometric conditions and temperatures (drought episodes of 1970-1990) associated to anthropic pressure on the plant and animal resources leading to the disappearance of tsetse forest habitats and wild hosts. The same factors lead to the disappearance of G. morsitans submorsitans in Dédougou and Koumbia [29]. Riverine flies, however, seem more resilient thanks to an opportunistic feeding behaviour associated to learning capacities [30] and linear dispersal along rivers [31].
The reduction of the mean age of the tsetse population in Koumbia indicates a higher adult mortality, which is probably a result of the perturbation/deterioration of their habitat. The reduction of the lifespan also favors trypanosome species having a short extrinsic cycle life as T. vivax (10 days), in comparison to those with longer cycle as T. congolense (14 days) and T. brucei brucei (30 days) [32].
Thus, this study confirmed that aridity and landscape fragmentation are associated to a reduction of the biodiversity of hosts, cyclical vectors and transmitted parasites in Burkina Faso [6]. Moreover, the impact of tsetse appeared not limited to their distribution area, but extended to a buffer area corresponding to the transhumance radius of herds, where transmission is probably relayed by mechanical vectors. The situation documented here only partially contributes to characterizing the impact of production systems and climate on tsetse populations and prevalence of AAT in Burkina Faso. It would be interesting to investigate whether similar situations are encountered elsewhere and using other methods (namely monitoring of AAT incidence).