Rift Valley fever virus dynamics in a transhumant cattle system in The Gambia

  1. Essa Jarra  Is a corresponding author
  2. Divine Ekwem
  3. Sarah Cleaveland
  4. Daniel T Haydon  Is a corresponding author
  1. Department of Livestock Services, Ministry of Agriculture, The Gambia
  2. School of Biodiversity, One Health, and Veterinary Medicine, College of Medical, Veterinary and Life Sciences, University of Glasgow, United Kingdom
8 figures, 4 tables and 1 additional file

Figures

Figure 1 with 3 supplements
Distribution of (A) degree, (B) normalised betweenness, and (C) eigenvector centrality values based on cross-sectional household survey data on livestock movement.

Only a few households showed high centrality values, highlighting their importance in the network.

Figure 1—figure supplement 1
A topological network depicting the overall connectivity among study households based on shared resource areas in The Gambia.

The network is undirected and unweighted, with grey edges. Coloured circles denote individual households within specific villages, with varying hues representing different study villages. White circles indicate shared water and grazing locations used by households. Squares represent transhumant household herds, which seasonally migrate over long distances (up to 80 km), establishing connections with household herds across villages. The transhumant movements are directed towards the Gambia river valley from within the Sahelian eco-region.

Figure 1—figure supplement 2
Spearman’s rank correlation analysis illustrating the relationship between household (HH) Rift Valley fever virus (RVFV) seropositivity in cattle and unweighted, normalized (A) degree, (B) betweenness, and (C) eigenvector centrality (computed only on the giant component) values derived from this study.
Figure 1—figure supplement 3
Spearman’s rank correlation analysis illustrating the relationship between the geographic distance separating connected households (HH) and their mean Rift Valley fever virus (RVFV) seropositivity of cattle in this study.
Time series of the proportion of infectious cattle in the three subpopulations across the two eco-regions.

(A) The simulated Rift Valley fever (RVF) infection dynamics from the deterministic model (solid lines) together with 20 realisations of the stochastic model (dashed lines). (B) A magnified view highlighting the seasonal dynamics of Rift Valley fever virus (RVFV) transmission focused on the last 2 years of the 20-year deterministic simulation in the cattle subpopulations, highlighting finer seasonal variations (wet season = beige; dry season = cyan). The proportion of infectious cattle peaked at the latter part of the wet season, but infections quickly disappeared in the dry season in the Sahelian eco-region. Once transhumant herds begin to arrive in the Gambia river region, infections are predicted to rise. I=infectious cattle in each subpopulation. A square root transformation was applied to the y-axis for visualisation purposes using coord_trans(y = ‘sqrt’), while the data remained in its original scale.

Figure 3 with 1 supplement
Extinction rate of Rift Valley fever virus (RVFV) over time (red line with 95% CI – grey ribbon), in the transhumant subpopulation (A) and the entire cattle population within the system (B) based on 1000 stochastic realisations, illustrating differences in the timing of local and system-wide extinctions, respectively.

Most of the local extinction occurred shortly after RVFV introduction into a fully susceptible population in the T subpopulation. Note: the timing of the local extinction in the T subpopulation depicted here represents the first extinction event within this subpopulation; re-infection occurs when the subpopulation returned to the river.

Figure 3—figure supplement 1
Distribution of the number of local extinction events in the T subpopulation over an average transmission period of 13.7 years, based on stochastic simulations.

The x-axis shows the number of extinction events per simulation, and the y-axis shows the count of simulations exhibiting each number of extinctions.

The full 10-year simulation of the weekly Rst (green) and force of infection (λi,season) in each eco-region (blue = Sahelian, red = Gambia river eco-region).

Shaded areas correspond to the seasons (wet season = beige; dry season = cyan) in the last 2 years of the simulation. A square root transformation was applied to the y-axis for visualisation purposes using coord_trans(y = ‘sqrt’), while the data remained in its original scale.

Predicted Rift Valley fever virus (RVFV) seroprevalence in each population subject to force of infection (FOI) during the wet and dry seasons.

(A) Cumulative growth of RVFV immunity over a 10-year period, representing the hypothetical lifespan of cattle, with no consideration for decay of RVFV seropositivity (π=0). (B) Proportion of immune cattle after introducing a seropositivity decay parameter (π=0.00217), aligning the predicted seroprevalence with observed data. Observed seroprevalence across different age classes in the three structured cattle populations in The Gambia is shown as dots, with 95% confidence intervals (CI) as error bars. P: cattle exposed to RVFV infection and recovered, assumed to be seropositive.

The general linear model (GLM) (solid line) and Loess (dotted line) smoothing plots of the relationship between percentage change in predicted seroprevalence of each structured subpopulation (M, L, and T) and percentage change in parameter values constrained within ±20% of the mean posterior value.
Map of study location – The Gambia.

(A) Illustrates the location of villages selected for household survey and ruminant sampling. Map created in QGIS 3.28.3; land cover data derived from Global Land Cover 2000 (GLC2000) Project. (B) Spatial range of transhumant cattle movements in The Gambia identified from the household survey. Movements originated from study villages (represented as red squares) and extended to destination villages, which are either other study villages or villages not selected for this study (represented as black triangles). This study was conducted during the dry season, and all transhumant herds were sampled at their destination villages. Distinct directional movements between homestead villages and their respective destinations are illustrated with blue arrows. The Gambia river is depicted as the meandering white line running the length of the country. The map was generated using the igraph package in R.

Epidemiological model of Rift Valley fever virus (RVFV) transmission and infection dynamics among the cattle subpopulations in The Gambia.

(A) Schematic representation of ecoclimatic region and seasonal combinations that influence RVFV transmission between the Sahelian areas and the Gambia river eco-region. (B) The transmission matrices that determined possible RVFV transmissions between the subpopulations during each season. (C) Transmission framework of the within eco-region RVFV transmission; parameters are defined within the main text.

Tables

Table 1
Mean posterior values and 95% credible intervals (CrI) of the estimated parameters.
Parameter descriptionNotationUnitMean value95% CrI
Per capita birth ratebday–10.00220.0012–0.0033
Per capita natural death rateμday–10.00170.0010–0.0023
RVF specific mortality rateγday–10.07580.0358–0.1868
RVF recovery rateδday–10.11610.0500–0.1917
Scaling factorψ2.05080.9224–3.1282
Wet season transmission parameter in the Sahelian eco-regionβs,wet*1.4×10–66.0×10–7 – 2.2×10–6
Wet season transmission parameter in the Gambia river eco-regionβr,wet*3.7×10–61.9×10–6 – 5.3×10–6
Dry season transmission parameter in the Gambia river eco-regionβr,dry*3.1×10–61.7×10–6 – 4.4×10–6
  1. *

    (Infected·susceptible·day)⁻¹.

Table 2
Summary of estimated seasonal reproduction number (Rst) and regio-specific force of infection λi.
SeasonsMean RstMaximum RstMinimum RstAvg. weeks Rst>1
Wet1.151.400.919/17
Dry1.021.110.8723/35
SeasonsRegionMean λiMaximum λiMinimum λi
WetSahelian0.0080.050.0001
river0.0030.0050.0002
DrySahelian
river0.0020.0180.0001
Table 3
The elasticity for all eight parameters in the general linear model (GLM).

The parameters with the largest coefficient are shown in bold.

Subpopulation MSubpopulation LSubpopulation T
ParameterInterceptCoefficientInterceptCoefficientInterceptCoefficient
βs,wet–0.0040.140–0.005–0.0270.008–0.032
βr,wet–0.0140.0180.064–1.0230.059–1.485
βr,dry–0.0280.1800.0741.0100.0561.604
μ–0.104–3.9650.151–4.7310.091–4.251
b–0.1543.6030.1084.7840.0594.230
γ–0.001–1.0160.055–1.7350.057–1.875
δ–0.0190.8740.0591.5000.0731.571
τ–0.0040.140–0.005–0.0270.008–0.032
Table 4
Preliminary parameter values for the susceptible-infectious-recovered (SIR) model.
Parameter descriptionNotationValueUnitReference
Per capita birth rateb0.00215day–1Gachohi et al., 2016
Per capita natural death rateμ0.0016day–1Nicolas et al., 2014
RVF specific mortality rateγ0.075day–1Gachohi et al., 2016
RVF recovery rateδ1/8day–1Durand et al., 2020
Scaling factorψ2Assumed
Decay rate of RVFV seropositivityπ0.005week–1Assumed
Total M populationNM177,0042016 Livestock Census
Total L populationNL80,4332016 Livestock Census
Total T populationNT35,4002016 Livestock Census

Additional files

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Essa Jarra
  2. Divine Ekwem
  3. Sarah Cleaveland
  4. Daniel T Haydon
(2026)
Rift Valley fever virus dynamics in a transhumant cattle system in The Gambia
eLife 14:RP107346.
https://doi.org/10.7554/eLife.107346.3