Influenza A Gradual and Epochal Evolution: Insights from Simple Models

The recurrence of influenza A epidemics has originally been explained by a “continuous antigenic drift” scenario. Recently, it has been shown that if genetic drift is gradual, the evolution of influenza A main antigen, the haemagglutinin, is punctuated. As a consequence, it has been suggested that influenza A dynamics at the population level should be approximated by a serial model. Here, simple models are used to test whether a serial model requires gradual antigenic drift within groups of strains with the same antigenic properties (antigenic clusters). We compare the effect of status based and history based frameworks and the influence of reduced susceptibility and infectivity assumptions on the transient dynamics of antigenic clusters. Our results reveal that the replacement of a resident antigenic cluster by a mutant cluster, as observed in data, is reproduced only by the status based model integrating the reduced infectivity assumption. This combination of assumptions is useful to overcome the otherwise extremely high model dimensionality of models incorporating many strains, but relies on a biological hypothesis not obviously satisfied. Our findings finally suggest the dynamical importance of gradual antigenic drift even in the presence of punctuated immune escape. A more regular renewal of susceptible pool than the one implemented in a serial model should be part of a minimal theory for influenza at the population level.     

Preparation and in vivo toxicity study of solid lipid microparticles as carrier for pulmonary administration

The purpose of this research was to investigate the effects of processing conditions on the characteristics of solid lipid microparticles (SLM) with a potential application as carriers for pulmonary administration. Compritol (5.0% wt/wt) SLM dispersions were prepared by rotor-stator homogenization, at different surfactant concentrations and emulsification times. The SLM were characterized, in terms of morphology and size, after lyophilization and sterilization by autoclaving process. In vivo assessment was carried out in rats by intratracheal instillation of either placebo or SLM dispersion, and by bronchoalveolar lavage for cytological analysis. Mean particle size of 4 to 5 μm was achieved using 0.3% and 0.4% (wt/wt) of emulsifier (Poloxamer 188) and emulsification times of 2 and 5 minutes. The particles showed spherical shape and smooth surface. The morphology of microparticles, the size, and the size distribution were not substantially modified after lyophilization and sterilization. Total cell counts showed no significant differences between placebo and SLM 0.5% or 2.5% groups. Regarding cytology, percentage of polymorphonuclear neutrophils and macrophages did not significantly differ between groups. These results suggest that a single intratracheal administration of the SLMs does not induce a significant inflammatory airway response in rats and that the SLMs might be a potential carrier for encapsulated drug via the pulmonary route.     

Does the Effectiveness of Control Measures Depend on the Influenza Pandemic Profile?

Background

Although strategies to contain influenza pandemics are well studied, the characterization and the implications of different geographical and temporal diffusion patterns of the pandemic have been given less attention.

Methodology/Main Findings

Using a well-documented metapopulation model incorporating air travel between 52 major world cities, we identified potential influenza pandemic diffusion profiles and examined how the impact of interventions might be affected by this heterogeneity. Clustering methods applied to a set of pandemic simulations, characterized by seven parameters related to the conditions of emergence that were varied following Latin hypercube sampling, were used to identify six pandemic profiles exhibiting different characteristics notably in terms of global burden (from 415 to >160 million of cases) and duration (from 26 to 360 days). A multivariate sensitivity analysis showed that the transmission rate and proportion of susceptibles have a strong impact on the pandemic diffusion. The correlation between interventions and pandemic outcomes were analyzed for two specific profiles: a fast, massive pandemic and a slow building, long-lasting one. In both cases, the date of introduction for five control measures (masks, isolation, prophylactic or therapeutic use of antivirals, vaccination) correlated strongly with pandemic outcomes. Conversely, the coverage and efficacy of these interventions only moderately correlated with pandemic outcomes in the case of a massive pandemic. Pre-pandemic vaccination influenced pandemic outcomes in both profiles, while travel restriction was the only measure without any measurable effect in either.

Conclusions

Our study highlights: (i) the great heterogeneity in possible profiles of a future influenza pandemic; (ii) the value of being well prepared in every country since a pandemic may have heavy consequences wherever and whenever it starts; (iii) the need to quickly implement control measures and even to anticipate pandemic emergence through pre-pandemic vaccination; and (iv) the value of combining all available control measures except perhaps travel restrictions.     

Using singular perturbations to reduce an epidemiological model: application to bovine viral diarrhoea virus within-herd spread

Studying the spread of a pathogen in a managed metapopulation such as cattle herds in a geographical region often requires to take into account both the within- and between-herd transmission dynamics. This can lead to high-dimensional metapopulation systems resulting from the coupling of several within-herd transmission models. To tackle this problem, we aim in this paper at reducing the dimension of a tractable but realistic dynamical system reproducing the within-herd spread. The context chosen to illustrate our purpose is bovine viral diarrhoea virus (BVDV) transmission in a cattle herd structured in two age classes and several epidemiological states, including two infectious states (transiently and persistently infected). Different time scales, corresponding to the epidemiological and demographic processes, are identified which allow to build a reduced model. Singular perturbation technique is used to prove that, under some non-restrictive conditions on parameter values, the behaviour of the original system is quite accurately approximated by that of the reduced system. Simulations are also performed to corroborate the approximation quality. Our study illustrates the methodological interest of using singular perturbations to reduce model complexity. It also rigorously proves the biologically intuitive assumption that transiently infected individuals can be neglected in a homogeneous population, when capturing the global dynamics of BVDV spread.     

Impact of the Infection Period Distribution on the Epidemic Spread in a Metapopulation Model

Epidemic models usually rely on the assumption of exponentially distributed sojourn times in infectious states. This is sometimes an acceptable approximation, but it is generally not realistic and it may influence the epidemic dynamics as it has already been shown in one population. Here, we explore the consequences of choosing constant or gamma-distributed infectious periods in a metapopulation context. For two coupled populations, we show that the probability of generating no secondary infections is the largest for most parameter values if the infectious period follows an exponential distribution, and we identify special cases where, inversely, the infection is more prone to extinction in early phases for constant infection durations. The impact of the infection duration distribution on the epidemic dynamics of many connected populations is studied by simulation and sensitivity analysis, taking into account the potential interactions with other factors. The analysis based on the average nonextinct epidemic trajectories shows that their sensitivity to the assumption on the infectious period distribution mostly depends on , the mean infection duration and the network structure. This study shows that the effect of assuming exponential distribution for infection periods instead of more realistic distributions varies with respect to the output of interest and to other factors. Ultimately it highlights the risk of misleading recommendations based on modelling results when models including exponential infection durations are used for practical purposes.     

Modelling the effect of heterogeneity of shedding on the within herd Coxiella burnetii spread and identification of key parameters by sensitivity analysis

Coxiella burnetii is the bacterium responsible for Q fever, a worldwide zoonosis. Ruminants, especially cattle, are recognized as the most important source of human infections. Although a great heterogeneity between shedder cows has been described, no previous studies have determined which features such as shedding route and duration or the quantity of bacteria shed have the strongest impact on the environmental contamination and thus on the zoonotic risk. Our objective was to identify key parameters whose variation highly influences C. burnetii spread within a dairy cattle herd, especially those related to the heterogeneity of shedding.To compare the impact of epidemiological parameters on different dynamical aspects of C. burnetii infection, we performed a sensitivity analysis on an original stochastic model describing the bacterium spread and representing the individual variability of the shedding duration, routes and intensity as well as herd demography. This sensitivity analysis consisted of a principal component analysis followed by an ANOVA. Our findings show that the most influential parameters are the probability distribution governing the levels of shedding, especially in vaginal mucus and faeces, the characteristics of the bacterium in the environment (i.e. its survival and the fraction of bacteria shed reaching the environment), and some physiological parameters related to the intermittency of shedding (transition probability from a non-shedding infected state to a shedding state) or to the transition from one type of shedder to another one (transition probability from a seronegative shedding state to a seropositive shedding state).Our study is crucial for the understanding of the dynamics of C. burnetii infection and optimization of control measures. Indeed, as control measures should impact the parameters influencing the bacterium spread most, our model can now be used to assess the effectiveness of different control strategies of Q fever within dairy cattle herds.     

Inference in Gaussian state-space models with mixed effects for multiple epidemic dynamics

Journal of Mathematical Biology - Respiratory syncytial virus (RSV) is a leading cause of acute lower respiratory tract infection worldwide, resulting in approximately sixty thousand annual...     

Spread of Q fever within dairy cattle herds: key parameters inferred using a Bayesian approach

Q fever is a worldwide zoonosis caused by Coxiella burnetii. Although ruminants are recognized as the most important source of human infection, no previous studies have focused on assessing the characteristics of the bacterial spread within a cattle herd and no epidemic model has been proposed in this context. We assess the key epidemiological parameters from field data in a Bayesian framework that takes into account the available knowledge, missing data and the uncertainty of the observation process owing to the imperfection of diagnostic tests. We propose an original individual-based Markovian model in discrete time describing the evolution of the infection for each animal. Markov chain Monte Carlo methodology is used to estimate parameters of interest from data consisting of individual health states of 217 cows of five chronically infected dairy herds sampled every week for a four-week period. Outputs are the posterior distributions of the probabilities of transition between health states and of the environmental bacterial load. Our findings show that some herds are characterized by a very low infection risk while others have a mild infection risk and a non-negligible intermittent shedding probability. Moreover, the antibody status seems to be a key point in the bacterial spread (shedders with antibodies shed for a longer period of time than shedders without antibodies). In addition to the biological insights, these estimates also provide information for calibrating simulation models to assess control strategies for C. burnetii infection.     

A mechanistic and data-driven reconstruction of the time-varying reproduction number: Application to the COVID-19 epidemic

The effective reproduction number R_eff is a critical epidemiological parameter that characterizes the transmissibility of a pathogen. However, this parameter is difficult to estimate in the presence of silent transmission and/or significant temporal variation in case reporting. This variation can occur due to the lack of timely or appropriate testing, public health interventions and/or changes in human behavior during an epidemic. This is exactly the situation we are confronted with during this COVID-19 pandemic. In this work, we propose to estimate R_eff for the SARS-CoV-2 (the etiological agent of the COVID-19), based on a model of its propagation considering a time-varying transmission rate. This rate is modeled by a Brownian diffusion process embedded in a stochastic model. The model is then fitted by Bayesian inference (particle Markov Chain Monte Carlo method) using multiple well-documented hospital datasets from several regions in France and in Ireland. This mechanistic modeling framework enables us to reconstruct the temporal evolution of the transmission rate of the COVID-19 based only on the available data. Except for the specific model structure, it is non-specifically assumed that the transmission rate follows a basic stochastic process constrained by the observations. This approach allows us to follow both the course of the COVID-19 epidemic and the temporal evolution of its R_eff(t). Besides, it allows to assess and to interpret the evolution of transmission with respect to the mitigation strategies implemented to control the epidemic waves in France and in Ireland. We can thus estimate a reduction of more than 80% for the first wave in all the studied regions but a smaller reduction for the second wave when the epidemic was less active, around 45% in France but just 20% in Ireland. For the third wave in Ireland the reduction was again significant (>70%).     

Apoptosis and survival of osteoblast-like cells are regulated by surface attachment

We tested the hypothesis that RGDS peptides regulate osteoblast survival in culture. Osteoblast-like MC3T3-E1 cells were allowed to attach to RGDS peptides that had been tethered to a silicone surface utilizing a previously described grafting technique. The RGDS-modified surface caused up-regulation of alpha(v)beta(3) integrin. We noted that there was an increase in expression of activated focal adhesion kinase and activated Akt. There was no change in the expression level of the anti-apoptotic protein Bcl-2, the pro-apoptotic protein Bad, or the inactivated form of Bad, pBad. Attachment to the RGDS-treated membrane completely abolished apoptosis induced by staurosporine, the Ca(2+).P(i) ion pair, and sodium nitroprusside. However, the surface modification did not interfere with apoptosis mediated by the free RGDS peptide or serum-free medium. When the activity of the phosphatidylinositol 3-kinase pathway was inhibited, RGDS-dependent resistance to apoptosis was eliminated. These results indicated that the binding of cells to RGDS abrogated apoptosis via the mitochondrial pathway and that the suppression of apoptosis was dependent on the activity of phosphatidylinositol 3-kinase.