MERS-Coronavirus across Kenya: a spatial examination of social and environmental drivers

1. Introduction

A combination of interacting social and environmental factors often drives the emergence and spread of infectious diseases (Heffernan, 2018). In particular, climate and agricultural land use intersect in ways that influence environmental conditions (Brodie, 2016), which can trigger the emergence and spread of infectious diseases. Specifically, ambient temperature affects infection rates, reproduction, and incubation time of pathogens, with higher temperatures accelerating pathogen maturation (Baker et al., 2022)(Semenza et al., 2022). Further, agriculture for food production has been associated with more than 25% of all infectious diseases and more than 50% of all zoonotic infectious diseases that have emerged in humans (Rohr et al., 2019). Lastly, the loss of biodiversity and the ensuing loss of host heterogeneity due in part to climate and agricultural land-use change has been linked to disease susceptibility and transfer (Heffernan, 2018).

In East Africa, changes in climate (i.e., hotter and drier trends) and agricultural land-use change have increased the likelihood of infectious disease emergence and transmission, such as Ebola virus, Flaviviruses, Usutu viruses, Chikungunya and O’nyong-nyong viruses, Bunyaviruses, and Rift Valley Fever and Crimean-Congo haemorrhagic viruses (Duygu et al., 2018; Fanelli & Buonavoglia, 2021; Fenollar & Mediannikov, 2018; Muturi et al., 2023; Pandit et al., 2022). Further, a combination of climate, agricultural and economic changes are supporting the spread of emerging pathogens from East Africa into the Middle East and Europe (Ryan et al., 2019; Xiao et al., 2015). Still, East Africa is considered one of the most at-risk regions in Africa to the impacts of climate change as the livelihoods of a large proportion of the region’s population depends on rain-fed agriculture (Serdeczny et al., 2017). Agro-pastoralists who depend on both livestock keeping and rain-fed crop production are considered the most vulnerable groups to climate change (Hughes & Anderson, 2020). Concurrently, population growth rates in East Africa are among the highest in the world, which increases the pressure for land conversion, and specifically for the expansion of cropland that encroaches on wildlife habitat (Bullock et al., 2021). In all, hotter and drier conditions, population growth, cropland expansion, and encroachment on wildlife habitats exacerbates infectious disease emergence and transmission (Lee-Cruz et al., 2021).

Despite the synergistic influence that climate and agricultural land use has on emerging infectious diseases, these drivers are often empirically examined separately and the potential synergistic effects are often missed (Brodie, 2016). Seldom, the influence of climate and agricultural land use have been examined related to emerging infectious diseases in a spatially explicit way at regional scales. This is especially true for Middle East Respiratory Syndrome coronavirus (MERS-CoV) that is prevalent in dromedary camels across Kenya where climate and agricultural land-use change are conspicuous (Lawrence et al., 2023a,b).

MERS-CoV is an infectious zoonotic disease that in humans targets the lower respiratory tract and can lead to multi-organ failure, resulting in death. Dromedary camels have been shown to be a natural reservoir of the virus from where spill-over to humans can occur (Adney et al, 2014). The infection is spread when a person comes into close contact with an infected dromedary camel, or possibly when a person consumes contaminated camel products such as milk and meat. Humans can also spread the virus to each other through very-close contact with infected individuals, similar to the current SARS-CoV-2 transmission (Aguanno et al., 2018).

The disease was first detected in humans in 2012 in Saudi Arabia (WHO, 2019). MERS-CoV has spread globally with more than 2,000 human infections resulting in nearly 850 identified deaths in 27 countries across North America, Europe, Asia, and Africa as of December 2019; while no new morbidity and mortality have been reported in the past several years (WHO, 2018). While no human infections have been documented in Kenya previously, the biophysical environment provides an opportunity to predict the conditions in which this may occur. Infection rates among and between camels and humans have been investigated throughout the Middle East (Reeves, 2015), Africa (Gikonyo et al., 2018; Miguel et al., 2017), and Asia (Saqib et al., 2017), which has provided preliminary mapping of infection risk. Since the discovery of MERS-CoV, serological and molecular evidence have demonstrated that the virus in dromedary camels is genetically similar to the one occurring in humans confirming the hypothesis that dromedary camels are the primary transmission reservoirs, which sheds the virus in high numbers and likely serve as reservoirs for human infections (Adney et al, 2014).

Our objective in this study was to spatially examine the climate, agriculture, and socio-demographic factors related to agro-pastoralism that can influence the prevalence of MERS-CoV in dromedary camels across northern Kenya. Our research questions were: 1) how has MERS-CoV in dromedary camels varied across geographic regions of northern Kenya, and 2) what climate, agriculture, and socio-demographic factors of agro-pastoralism were spatially related to the geographic variation in MERS-CoV cases? To answer our questions, we analyzed the spatial distribution of historical cases of serological evidence of MERS-CoV at the county level and applied spatial statistical analysis to examine the spatial relationships of the MERS-CoV cases between 2016 and 2018 to climate, agriculture, and socio-demographic factors of agro-pastoralism.

2. Study Site, Data Description, and Methods

2.3 Data analysis, model evaluation and selection

We analyzed the data via a combination of non-spatial and spatial analyses using R (R Core Team, 2013), as summarized in Figure 1. Initially, we mapped MERS-CoV cases at the county level in northern Kenya using ArcGIS Pro 2.9.2 (ESRI 2022) and examined seropositivity relative to the total samples across the region, upper Rift Valley, Eastern, and Northeastern regions during the study period. In the upper Rift Valley region, we included the counties of Turkana, Laikipia, Baringo, and W. Pokot. In the upper Eastern region, we included the counties of Marsabit, Isiolo, and Samburu. In the upper Northeastern region, we included the counties of Mandera and Wajir. We then formally investigated and quantified spatial correlation of and between variables using variographic analysis, which decomposes the spatial variability of observed variables among distance classes. In this process, first, we examined the spatial autocorrelation of MERS-CoV cases. Next, we examined spatial correlation between MERS-CoV cases and a) socio-demographic variables related to agro-pastoralism, b) agriculture, and c) climate variables (Appendix C). Given that much of the socio-demographic data were from after 2010 and that the MERS-CoV data in the upper Northeastern region were only from 2008, we focused the spatial correlation analysis of MERS-CoV on 2016 through 2018 in the upper Rift Valley and upper Eastern Regions. Further, we distributed the data evenly between the northern and southern parts of each county to satisfy the practical rule that a variogram should only be applied over a specified distance for which the number of pairs is greater than 30, and because the data were not further spatially identified (Journel and Huijbregts, 1978). We fitted the data using spherical and exponential models and used a distance of between 250 to 350 km, which was deemed conservative, based on the maximum distance of 707 km for the combined upper Rift Valley and Eastern regions of northern Kenya (Journel & Huijbregts, 1978; Crawley, 2013).

• Download figure
• Open in new tab

Figure 1.

Summary process of spatial modeling and analysis of MERS-CoV and climate, agriculture, and socio-demographic factors related to agro-pastoralism that can influence the prevalence. The initial step involved mapping MERS-CoV cases and examining relative seropositivity across northern parts of regions across Kenya. The next step involved quantifying spatial correlation between the variables. Finally, simultaneous autoregressive models were tested.

After confirming spatial correlation among the MERS-CoV cases and socio-environmental variables, we specified simultaneous autoregressive (SAR) models, a statistical method that augments linear regression models with an additional term to account for the spatial correlation structure in a dataset (Kissling & Carl, 2008). To include the spatial correlation structure of our dataset into the SAR models, we defined neighbors among the northern and southern parts of and between each county based on shared borders, and created a spatially weighted matrix. Using shared borders to define neighbors, rather than including counties beyond those with shared borders, allowed us to account for spatial correlation if it diminished over an increasing distance. We weighted each county’s neighbor equally, such that the weights of all neighbors of a sub-county summed to one. Equation 1 shows the general SAR model in matrix form that includes the spatial structure of our dataset.

Our analyses involved testing SAR models in a stepwise process. First, we tested MERS-CoV cases against individual predictor variables. The predictor variables that were statistically significant and spatial correlated with MERS-CoV cases were retained for further testing in the next set of models. The next set of models tested began with the predictor variable with the lowest AIC in the previous set of models tested. Each of the other statistically significant predictor variables from the previously tested set of models were then individually tested in the new best model at that point in the process. The stepwise testing continued until all possible variables were tested while retaining only the predictor variables that were statistically significant, exhibited spatial correlation, and had the lowest AIC. We evaluated and compared the SAR models relative to each other using the p-value (with at least a 90% level) of the likelihood ratio test where a model with no spatial correlation (i.e., λ = 0) is compared to the fitted model with a non-zero spatial correlation parameter (Kissling & Carl, 2008). Ultimately, we used the Akaike information criterion (AIC) to choose the best performing model.

3. Results

The total number of samples and seropositive cases between 1992 through 2018 in the upper Rift Valley (Turkana, Baringo, W. Pokot, Laikipia) (n=2229, pos.= 914) had a 41% positivity, in the upper Eastern region (Marsabit, Isiolo, Samburu) (n=2486, pos.=1876) had a 75% positivity, and in the upper North Eastern region (Mandera, Wajir) (n=3506, pos.=1983) had a 57% positivity (Figure 2). The total number of samples and seropositive cases between 2016 through 2018 in the upper Rift Valley (n=935, pos.= 414) had a 51% positivity, and the upper Eastern region (n=1964, pos.=1510) had a 77% positivity. Spatial correlation of MERS-CoV extended a distance of roughly 400 kilometers across the upper Rift Valley and Eastern regions between 2016 and 2018.

• Download figure
• Open in new tab

Figure 2.

Geographic distribution and spatial correlation of MERS-CoV across northern regions of Kenya. The upper Rift Valley (counties of Turkana, Baringo, W. Pokot, Laikipia) had an average seropositivity of 41% between 1992 and 2018, and an average seropositivity of 51% between 2016 and 2018. The upper Eastern region (counties Marsabit, Isiolo, Samburu) had an average seropositivity of 75% between 1992 and 2018, and an average seropositivity of 77% between 2016 and 2018. The upper Northeastern region (counties of Mandera and Wajir) had a seropositivity of 57% in 2008. The spatial autocorrelation of MERS-CoV extends a distance of approximately 400 km in 2016 to 2018 across the upper Rift Valley and Eastern regions of northern Kenya.

The climate, agriculture, and socio-demographic variables (Appendix C) that were individually spatially correlated with MERS-CoV cases between 2016 and 2018 and statistically significant according to the p-value were included in the SAR models. While all of the six spatial covariates included in the SAR models did not achieve p < 0.10, they were deemed close to the cutoff and important theoretically to assess the spatial relationship (Table 2). According to the AIC, the spatial relationship of ethno-religious camel practices with MERS-CoV cases was significantly better than the other independent variables. Specifically, the p-value (p=0.015) was below the significance threshold of p < 0.05 and the AIC was statistically different from all the other individually tested independent variables with a -24.92 change in the AIC from the weakest performing SAR model, which was the spatial relationship of human population with MERS-CoV cases.

The multi-variate SAR models that were statistically significant and spatially correlated with MERS-CoV cases are shown in Table 3, while all multi-variate SAR models that were tested are shown in Appendix D through Appendix H. The first synergistic effect appeared with the inclusion of agricultural land (model 9), which was statistically different from ethno-religious camel practices alone (model 6) based on the AIC, as well as spatially correlated with MERS-CoV cases. The SAR model improved, according to the AIC, and a broader synergistic effect demonstrated when ethnic diversity and agro-pastoral activities were included (model 19).

However, the best performing SAR model according to the AIC, and the most complex synergistic effect was when all of the six predictors (model 23) were included, and were spatially correlated with MERS-CoV cases and statistically significant. Thus, ethno-religious camel practices and agriculture land were the two initially interacting variables of importance. Next, ethnic diversity and agro-pastoral activities combined with the previous two independent variables were the next set of interacting variables of importance. Finally, average air temperature and human population combined with the previous four independent variables were the final set of interacting variables of importance.

4. Discussion

This study was conducted to collate data from disparate studies and datasets to provide greater insight to the prediction of MERS-CoV infections that may be facing Kenya, and regions with similar geographies across the world. Most often, the occurrence of infectious diseases has been studied independently, overlooking the potential combined effects of climate change, agricultural land use, and human-environment interaction. Here, we comprehensively and spatially evaluated the influence of various components that may contribute to the emergence of MERS-CoV throughout the northern regions of Kenya. This study highlighted the benefit of bringing together diverse areas of data and research to spatially inform infectious disease risk.

The geographical variations of serological evidence of MERS-CoV in dromedary camels and spatial correlations with the social and environmental variables may be partly attributed to the dominant ethnic groups in the different regions and their particular agro-pastoral management practices. Our results showed low MERS-CoV seropositivity in the Rift Valley region where tribal communities practice a more diversified form of livestock management with camels reared alongside other livestock, such as cattle, sheep, and goats (Iiyama et al., 2008).

Particularly, the Turkana communities and the Pokot have adapted livestock diversification as part of a long-term adaptation strategy to manage drought and diseases (Opiyo et al., 2015), and tend to acquire cattle and camels through cultural practices, such as dowry rather than directly from markets (De Vries et al., 2006). Thus, tribal communities, such as Pokot, Maasai or Turkana have a limited dependence on camels for their economic activity that can result in lower seropositivity of MER-CoV compared to other regions (Deem et al., 2015).

Tribal communities, such as the Somali, Gabra, and Borana in the eastern and northeastern regions of Kenya primarily focus on camels as a pastoral livelihood strategy due to the camel’s adaptability to arid environments and ability to withstand droughts (Ngere et al., 2020). Also, the Somali primarily inhabit the Northeastern region with high regard for camels for their annual religious practices and ritual migrations (Watson, 2010), and a greater number of cultural practices and traditions associated with livestock rearing, trade, and consumption in tribal communities (Dan et al., 2021; Kagunyu et al., 2018). Further, the tribal communities in both the Eastern and Northeastern regions tend to migrate to different areas in search of fresh water and pastures, and restock their camels through direct markets, which can lead to higher prevalence rates of MERS-CoV (De Vries et al., 2006; (Hughes & Anderson, 2020).

Recent findings have shown that the expanding arid regions of Kenya have reduced agricultural land use and increased the reliance on pastoralism, and particularly camel rearing (Lawrence et al., 2023a). Also, the increasing trend in air temperatures has resulted in decreased rainfed agricultural production, and further leading to increased pastoralism (Lawrence et al., 2023b). As a consequence, these effects of climate change may continue to pose a major challenge in creating favorable conditions for the emergence and transmission of viral zoonotic diseases resulting in the observed high seropositivity rates of MERS-CoV. Similar results of spatial clustering of zoonotic diseases have been studied in different agro-ecological climate zones and sociodemographics of Kenya to identify potential disease hotspots of anthrax (Nderitu et al., 2021). Overall, these findings are significant as shifting climatic zones continue to impact agricultural systems while zoonotic diseases like MERS-CoV pose a major risk to communities living in arid regions of Kenya warranting immediate interventions and continued surveillance.

Due to these climatic and agricultural shifts, there is an urgent need to predict where and under what conditions would viral transmission of MERS-CoV occur. These findings suggest that continued monitoring of camels, humans, and their spatial patterns will be an integral component of informing prediction models. This study relied on previously collected data that differed in methods and specificity of data, which also assumes equal probability of detecting seropositivity. Further, the county-level finding in this study can vary at smaller more local spatial scales. Thus, results are limited in power and yet, the broader spatial patterns observed and insights that this type of study can provide are important in the planning and implementation of prediction modeling, as well as location-based risk management.

5. Conclusion

Climate and agricultural land-use change have increased the likelihood of infectious disease emergence and transmissions, but these drivers are often examined separately as synergistic effects are ignored. Further, seldom are the influence of climate and agricultural land use on emerging infectious diseases examined in a spatially explicit way at regional scales. Our objective in this study was to spatially examine the climate, agriculture, and socio-demographic factors related to agro-pastoralism that can influence the prevalence of MERS-CoV in dromedary camels across northern Kenya. Overall, this study can provide important insights in the planning and implementation of prediction modeling, as well as location-based risk management.