No Arabic abstract
The unexpected Ebola virus outbreak in West Africa in 2014 involving the Zaire ebolavirus made clear that other regions outside Central Africa, its previously documented niche, were at risk of future epidemics. The complex transmission cycle and a lack of epidemiological data make mapping areas at risk of the disease challenging. We used a Geographic Information System-based multicriteria evaluation (GIS-MCE), a knowledge-based approach, to identify areas suitable for Ebola virus spillover to humans in regions of Guinea, Congo and Gabon where Ebola viruses already emerged. We identified environmental, climatic and anthropogenic risk factors and potential hosts from a literature review. Geographical data layers, representing risk factors, were combined to produce suitability maps of Ebola virus spillover at the landscape scale. Our maps show high spatial and temporal variability in the suitability for Ebola virus spillover at a fine regional scale. Reported spillover events fell in areas of intermediate to high suitability in our maps, and a sensitivity analysis showed that the maps produced were robust. There are still important gaps in our knowledge about what factors are associated with the risk of Ebola virus spillover. As more information becomes available, maps produced using the GIS-MCE approach can be easily updated to improve surveillance and the prevention of future outbreaks.
Identifying directed interactions between species from time series of their population densities has many uses in ecology. This key statistical task is equivalent to causal time series inference, which connects to the Granger causality (GC) concept: $x$ causes $y$ if $x$ improves the prediction of $y$ in a dynamic model. However, the entangled nature of nonlinear ecological systems has led to question the appropriateness of Granger causality, especially in its classical linear Multivariate AutoRegressive (MAR) model form. Convergent-cross mapping (CCM), a nonparametric method developed for deterministic dynamical systems, has been suggested as an alternative. Here, we show that linear GC and CCM are able to uncover interactions with surprisingly similar performance, for predator-prey cycles, 2-species deterministic (chaotic) or stochastic competition, as well as 10- and 20-species interaction networks. There is no correspondence between the degree of nonlinearity of the dynamics and which method performs best. Our results therefore imply that Granger causality, even in its linear MAR($p$) formulation, is a valid method for inferring interactions in nonlinear ecological networks; using GC or CCM (or both) can instead be decided based on the aims and specifics of the analysis.
Mathematical modelling has successfully been used to provide quantitative descriptions of many viral infections, but for the Ebola virus, which requires biosafety level 4 facilities for experimentation, modelling can play a crucial role. Ebola modelling efforts have primarily focused on in vivo virus kinetics, e.g., in animal models, to aid the development of antivirals and vaccines. But, thus far, these studies have not yielded a detailed specification of the infection cycle, which could provide a foundational description of the virus kinetics and thus a deeper understanding of their clinical manifestation. Here, we obtain a diverse experimental data set of the Ebola infection in vitro, and then make use of Bayesian inference methods to fully identify parameters in a mathematical model of the infection. Our results provide insights into the distribution of time an infected cell spends in the eclipse phase (the period between infection and the start of virus production), as well as the rate at which infectious virions lose infectivity. We suggest how these results can be used in future models to describe co-infection with defective interfering particles, which are an emerging alternative therapeutic.
Population dynamics and evolutionary genetics underly the structure of ecosystems, changing on the same timescale for interacting species with rapid turnover, such as virus (e.g. HIV) and immune response. Thus, an important problem in mathematical modeling is to connect ecology, evolution and genetics, which often have been treated separately. Here, extending analysis of multiple virus and immune response populations in a resource - prey (consumer) - predator model from Browne and Smith cite{browne2018dynamics}, we show that long term dynamics of viral mutants evolving resistance at distinct epitopes (viral proteins targeted by immune responses) are governed by epistasis in the virus fitness landscape. In particular, the stability of persistent equilibrium virus-immune (prey-predator) network structures, such as nested and one-to-one, and bifurcations are determined by a collection of circuits defined by combinations of viral fitnesses that are minimally additive within a hypercube of binary sequences representing all possible viral epitope sequences ordered according to immunodominance hierarchy. Numerical solutions of our ordinary differential equation system, along with an extended stochastic version including random mutation, demonstrate how pairwise or multiplicative epistatic interactions shape viral evolution against concurrent immune responses and convergence to the multi-variant steady state predicted by theoretical results. Furthermore, simulations illustrate how periodic infusions of subdominant immune responses can induce a bifurcation in the persistent viral strains, offering superior host outcome over an alternative strategy of immunotherapy with strongest immune response.
Persistent Organic Pollutants represent a global ecological concern due to their ability to accumulate in organisms and to spread species-by-species via feeding connections. In this work we focus on the estimation and simulation of the bioaccumulation dynamics of persistent pollutants in the marine ecosystem, and we apply the approach for reconstructing a model of PCBs bioaccumulation in the Adriatic sea, estimated after an extensive review of trophic and PCBs concentration data on Adriatic species. Our estimations evidence the occurrence of PCBs biomagnification in the Adriatic food web, together with a strong dependence of bioaccumulation on trophic dynamics and external factors like fishing activity.
We recently described a dynamic causal model of a COVID-19 outbreak within a single region. Here, we combine several of these (epidemic) models to create a (pandemic) model of viral spread among regions. Our focus is on a second wave of new cases that may result from loss of immunity--and the exchange of people between regions--and how mortality rates can be ameliorated under different strategic responses. In particular, we consider hard or soft social distancing strategies predicated on national (Federal) or regional (State) estimates of the prevalence of infection in the population. The modelling is demonstrated using timeseries of new cases and deaths from the United States to estimate the parameters of a factorial (compartmental) epidemiological model of each State and, crucially, coupling between States. Using Bayesian model reduction, we identify the effective connectivity between States that best explains the initial phases of the outbreak in the United States. Using the ensuing posterior parameter estimates, we then evaluate the likely outcomes of different policies in terms of mortality, working days lost due to lockdown and demands upon critical care. The provisional results of this modelling suggest that social distancing and loss of immunity are the two key factors that underwrite a return to endemic equilibrium.