Next Article in Journal
Disparities in Methods Used to Determine Microplastics in the Aquatic Environment: A Review of Legislation, Sampling Process and Instrumental Analysis
Previous Article in Journal
Italians Are the Fastest 3000 m Open-Water Master Swimmers in the World
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Multivariate Age-Structured Stochastic Model with Immunization Strategies to Describe Bronchiolitis Dynamics

by
Mónica López-Lacort
1,
Ana Corberán-Vallet
2,* and
Francisco J. Santonja Gómez
2
1
Vaccine Research Department, Fisabio-Public Health, Avda. Cataluña 21, 46020 Valencia, Spain
2
Department of Statistics and Operations Research, Universitat de València, Dr. Moliner 50, 46100 Burjassot, Spain
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2021, 18(14), 7607; https://doi.org/10.3390/ijerph18147607
Submission received: 14 May 2021 / Revised: 12 July 2021 / Accepted: 14 July 2021 / Published: 17 July 2021
(This article belongs to the Section Infectious Disease Epidemiology)

Abstract

:
Bronchiolitis has a high morbidity in children under 2 years old. Respiratory syncytial virus (RSV) is the most common pathogen causing the disease. At present, there is only a costly humanized monoclonal RSV-specific antibody to prevent RSV. However, different immunization strategies are being developed. Hence, evaluation and comparison of their impact is important for policymakers. The analysis of the disease with a Bayesian stochastic compartmental model provided an improved and more natural description of its dynamics. However, the consideration of different age groups is still needed, since disease transmission greatly varies with age. In this work, we propose a multivariate age-structured stochastic model to understand bronchiolitis dynamics in children younger than 2 years of age considering high-quality data from the Valencia health system integrated database. Our modeling approach combines ideas from compartmental models and Bayesian hierarchical Poisson models in a novel way. Finally, we develop an extension of the model that simulates the effect of potential newborn immunization scenarios on the burden of disease. We provide an app tool that estimates the expected reduction in bronchiolitis episodes for a range of different values of uptake and effectiveness.

1. Background

Bronchiolitis is a common lower respiratory tract infection (LRTI) that mainly affects children under 2 years old, with the greatest burden occurring in infants younger than 6 months [1]. Respiratory syncytial virus (RSV) is by far the most common cause of viral bronchiolitis [2]. Previous data from the Valencia region of Spain showed that approximately 2 out of 10 children younger than 2 years of age are diagnosed of bronchiolitis, 3 out of 100 are hospitalized and 1.6 out of 100 are hospitalized with RSV bronchiolitis [3,4].
At present, palivizumab is the only prophylactic therapy available for RSV. This is a costly humanized monoclonal RSV-specific antibody that is given monthly to infants at increased risk of severe RSV infection: premature babies, infants with chronic lung disease and infants with congenital heart disease [5,6]. However, different RSV immunization strategies, including maternal immunization, immunoprophylaxis with new monoclonal antibodies (mAbs) and pediatric immunization, are currently under development and could be available soon [7,8].
Immaturity of the infant’s immune system does not allow active immunization until 2 months of age. However, this population subgroup is more susceptible to severe RSV infection. Therefore, passive immunization practices (maternal vaccination and mAbs administered at birth) have the potential to be the most effective strategies to protect individuals at an early age. ResVax is the most advanced maternal vaccine in clinical development. In a phase III multi-country, randomized, placebo-controlled trial evaluating its efficacy against RSV-LRTI in infants from birth to 90–180 days of life, ResVax showed 44% efficacy in reducing RSV-LRTI hospitalization. However, it failed to meet the primary outcome of prevention of medically significant LRTI [9]. A phase III trial showed the potential of mAb nirsevimab to protect infants for an entire RSV season with just one injection. A single dose of this mAb with an extended half-life resulted in a lower incidence of medically attended RSV-associated LRTI and hospitalizations than with a placebo in healthy preterm infants entering their first RSV season. They were approximately 70% and 80% lower, respectively [10].
Before the licensure of these products for RSV and the implementation of a particular preventive strategy, policymakers will need to assess their potential health and economic benefits. This requires an exhaustive evaluation of the real burden of disease and the simulation of disease incidence under different immunization scenarios.
Different studies have been proposed to analyze either bronchiolitis or RSV dynamics. Most of these studies consider classical regression models to describe causal relationships between risk factors (baseline characteristics or clinical features) and disease counts [4,11,12,13,14,15]. Deterministic compartmental models based on differential equations have also been commonly implemented to model RSV dynamics [16,17,18,19]. Deterministic models assume that every individual has an equal probability of contacting every other individual in the population, so they may fail if the sizes of the compartments are not large enough to ensure homogeneous mixing. The incorporation of stochasticity into epidemic models is necessary for considering heterogeneous mixing of individuals in the population and also the inherent stochasticity in the transmission of disease [20,21,22,23]. A Bayesian stochastic compartmental model in discrete time has improved the description of RSV disease [24].
Alternative stochastic approaches for the analysis of time series of infectious disease counts have been described by Held et al. [25]. They explain disease counts by means of a Poisson (or negative binomial) model with two components: a parameter-driven component that describes endemic seasonal patterns and an observation-driven component, which is an autoregression on past counts that explains localized epidemics. This model can be extended in a multivariate setting to analyze time series that correspond to disease counts in different spatial units or different age groups jointly [26].
Systematic reviews of model-based evaluations of immunization strategies against RSV have been presented in [27,28]. Many of these studies use mathematical models based on ordinary differential equations (see, for instance, [17,29,30]). Bayesian stochastic compartmental models have also been extended to simulate and assess the effect of a vaccination strategy that consists of vaccinating a proportion of newborns [31]. Since clinical data are not yet available, models exploring the epidemiological and economic outcomes of potential RSV interventions adopt different assumptions regarding the target population, uptake rates, effectiveness and duration of immunity. Therefore, the estimates that are obtained vary considerably.
In this paper, we present a Bayesian multivariate age-structured model in discrete time to describe bronchiolitis dynamics in the region of Valencia in a population of children under 2 years of age. In particular, the population is divided into four age groups, and new infections are described considering interaction among these groups. We base our formulation on a Poisson model (or a negative binomial model if the data exhibit overdispersion). The main innovation is that the mean of the distribution at each time point can be seen as an observation-driven component where the autoregressive parameter varies stochastically over time. A parameter allowing for heterogeneous mixing of individuals in the population is also introduced to capture different contact patterns. This modeling approach does not require information on susceptibles, and it also provides a simplified framework to describe disease counts, since it avoids the modeling of complex transitions between the different compartments. To estimate the impact of potential prevention strategies for RSV, an extension of the model including simulation of bronchiolitis episodes under different passive immunization scenarios is considered.

2. Data

2.1. Population of Interest

The region of Valencia, one of the 17 autonomous regions of Spain, has approximately 4,900,000 inhabitants. Around 2% of the Valencia population is younger than 2 years old (approximately 100,000 children). The Regional Health System (RHS) is divided into 241 health care districts structured into 24 health departments. It includes 34 public hospitals, 24 of them attending acute pediatric patients. This study includes children under the age of 2 years who were born in the region of Valencia between January 2009 and December 2012 and are covered by the RHS.

2.2. Data Sources

The region of Valencia has a health system integrated database (VID) that gathers health and sociodemographic data from 98% of the population [32]. In particular, we used the population information system (an administrative database that collects and updates sociodemographic data from both residents and non-residents with access to public health services) to determine the population of interest. Primary care electronic medical notes (SIA) were implemented in 2006, and all medical visits are registered and ICD-coded. Hospitalizations were obtained using discharge reports from the Spanish Minimum Basic Data Set (MBDS).

2.3. Age-Structured Bronchiolitis Cases

Because around 90% of bronchiolitis cases in children less than 2 years of age are treated in primary care offices and very few cases have a RSV microbiological confirmation, the analysis of laboratory-confirmed cases of RSV bronchiolitis is less useful for incidence estimation due to underreporting. Hence, it is more practical to analyze bronchiolitis-associated outpatient visits and hospitalizations to evaluate the potential impact of RSV immunization strategies.
In particular, we analyze bronchiolitis episodes identified from hospitalization and primary care attendance through a search of the first appearance of the following ICD-9-codes 466.1, 466.11 and 466.19, in MBDS and SIA. Based on the McConnochie criterion, only the first health care encounter (either outpatient visit or hospitalization) with a bronchiolitis ICD-9-MC code was counted as a case of bronchiolitis.
Considering that children aged less than 6 months have an elevated morbidity, we consider four age groups: 0–5, 6–11, 12–17 and 18–23 months of age. This grouping allows us to analyze in greater detail both transmission dynamics and the impact of potential immunization measures. Weekly aggregated cases from July 2010 (week 28 June–4 July) to December 2012 (week 24 December–30 December) are considered (T = 131 weeks in total). Since this is a population-based cohort study, data for all the age groups were not available until July 2010, so data from January 2009 to June 2010 are excluded in the analysis. In total, the study includes 30,555 cases of bronchiolitis: 14,635 (47.90%) of them are in children aged less than 6 months, 10,600 (34.69%) cases in the age group of 6–11 months, 3462 (11.33%) in the age group of 12–17 months and 1858 (6.08%) in the age group of 18–23 months. Figure 1 shows the time plot of the series for the four age groups. As can be seen, there is a higher incidence for younger children and a clear seasonal pattern, with more cases during the weeks of autumn and winter. These data can be downloaded from https://rotapp.shinyapps.io/APP_IMPACTO_RSV, (accessed on 14 May 2021).

3. Model Description

3.1. Background

The most commonly used models in the analysis of infectious disease counts are compartmental models, which divide the population being studied into different compartments according to disease status and describe the evolution of infection through changes in the number of individuals in each compartment. Corberán et al. [24] proposed a Bayesian stochastic compartmental model in discrete time to describe RSV dynamics in the Region of Valencia. Let y t , t = 1 , 2 , , T , be the number of newly infected individuals at week t. The model assumes that:
y t B i ( S t 1 , p t ) ,
where S t 1 is the susceptible population at time t 1 , which is updated weekly using some recursion equations, and p t is the probability of becoming infected at time t. To take into account the transmissible nature of the infection, p t is modeled as:
p t = y t 1 h · exp ( r t ) 1 + y t 1 h · exp ( r t ) ,
where the mixing parameter h allows for heterogeneous mixing of individuals in the population and r t = α 0 + w = 1 W β w · sin 2 · π · w · t 52 +   γ w · cos 2 · π · w · t 52 + ε t . By allowing the transmission rate exp r t to vary over time, the stochastic model provided an improved and accurate description of the pattern of disease.
When the focus is on describing counts of new infections, the use of a Bayesian hierarchical model provides an alternative framework. Held et al. [25] proposed a stochastic model for the analysis of disease counts based on a Poisson (or negative binomial) model with two components, which describe endemic seasonal patterns and localized epidemics. The starting point was a simple branching process model, which was later extended to include seasonal terms in the endemic rate or to adjust for overdispersion. In particular, weekly new counts of disease y t , t = 1 , 2 , , T , are modeled as:
y t P o i s s o n ν t + λ t · y t 1 ,
where ν t is the parameter-driven (or endemic) component and λ t · y t 1 the observation-driven (or epidemic) component, which allows for occasional outbreaks. A main feature of the proposed model is that the autoregressive parameter λ is allowed to vary over time. To do this, a Bayesian changepoint model with an unknown number of changepoints is used to capture sudden changes in infectiousness. The endemic component is described as follows:
ν t = exp α 0 + β · sin 2 · π · t 52 + γ · cos 2 · π · t 52 .
Multivariate extensions of that model can be found in [26] for the joint analysis of multiple time series of counts, where each component corresponds to a geographical region or a certain age group. Let y t j denote the count of disease in age group j at week t. In the multivariate scenario, counts of disease can be described as:
y t j P o i s s o n ν t j + λ · y t 1 j , ν t j = exp α 0 j + α 1 · t + w = 1 W β w · sin 2 · π · w · t 52 + γ w · cos 2 · π · w · t 52 ,
where parameters α 0 j allow for different incidence levels in the different age groups. The epidemic component here depends only on previous counts in the corresponding age group.
A more general model considering previous counts in other age groups as potential explanatory variables is formulated as:
y t j P o i s s o n ν t j + λ · y t 1 j + ϕ · i j y t 1 i ,
where the additional parameter ϕ captures the autoregressive effect of the other age groups. Paul et al. [33] extended this multivariate model by allowing the autoregressive parameters to depend on the age-group; that is, λ j and ϕ j for each time series. However, in these multivariate models, the autoregressive parameters λ and ϕ are not allowed to vary over time.

3.2. Our Proposal

Let y t j denote the count of disease in age group j , j = 1 , 2 , 3 , 4 , observed at week t, t = 1 , 2 , ,   T = 131 . We use a Bayesian hierarchical model, which provides a straightforward framework to describe counts of new infections at each week. Our model assumes that the probability governing the counts is a Poisson distribution or a negative binomial distribution if the data exhibit overdispersion. In those situations in which the population size is not large enough in comparison with the observed counts of disease, the binomial distribution could be used in a similar way.
Seasonal variation in disease transmission that is persistent with a stable pattern is modeled through sine–cosine waves, which are common for all groups. This is a sensible assumption since annual epidemic peaks occur simultaneously in the different age groups. On the other hand, it is reasonable to assume that contacts among the different age groups considered here are equally likely to occur. Hence, counts of disease at time t 1 are summed up to describe new infections. In addition, heterogeneous mixing of individuals in the population is allowed. Finally, due to the epidemiology of bronchiolitis, it is important to consider a different transmission rate for each age group.
Taking into account these considerations, the proposed model is given by the following equations:
(5) y t j P o i s s o n μ t j , (6) μ t j = exp ( λ t j ) · i = 1 4 y t 1 i h , (7) λ t j = α 0 j + β · sin 2 · π · t 52 +   γ · cos 2 · π · t 52 + ε t ,
where the autoregressive parameter exp ( λ t j ) is allowed to vary stochastically over time by means of random effects ε t N ( 0 , σ ε 2 ) that represent unspecified features of week t. To avoid overfitting, these random effects are common for all the age groups; that is, if counts of disease at a particular week are higher (or lower) than expected for some age groups, that increase (or decrease) should also be observed in the remaining age groups. h [ 0 , 1 ] is the mixing parameter ( h = 1 would correspond to the assumption of mass action). Non-informative flat prior distributions are considered for parameters α 0 j , β and γ . The uniform distribution in the interval [0,1] is considered as a prior for h, and the uniform distribution in the interval (0,2) is considered for the standard deviation σ ε .
To adjust for overdispersion, the negative binomial distribution can be used instead. In that case, the model can be formulated as
(8) y t j N e g B i n ( p t j , k ) , (9) p t j = k k + μ t j ,
where μ t j is defined as in Equation (6). Hence, the mean of the distribution is equal to μ t j and the variance is given by μ t j + ( μ t j ) 2 k . The second parameter k > 0 incorporates the extra-Poisson variation (the limiting case k = corresponds to the Poisson distribution). As a prior distribution for parameter k, we consider the Gamma distribution G a ( 1 , 0.01 ) .
It is important to emphasize that our formulation can be seen as a Poisson model with only one component, the observation-driven component [25,26], modified in a novel way to allow for contact heterogeneity. By considering a time-varying autoregressive parameter, the model properly describes both endemic an epidemic periods. As shown in the next section, this formulation will allow us to simulate disease incidence in a scenario with immunization measures. On the other hand, the mean of the distribution is equivalent to that proposed in [24], where the process intensity mean was dependent on the product of S t 1 , y t 1 h and the transmission rate exp r t . In our modeling framework, where information on susceptibles is not incorporated, the autoregressive parameter exp ( λ t j ) describes how infected individuals in the previous week produce new infections.

3.3. Our Extension with Immunization Strategies

We show here how the proposed model can be used to simulate the impact of passive immunization through maternal vaccination or mAb administered at birth. The main parameters to consider when estimating the effects of a particular immunization strategy are the uptake, the efficacy and the duration of protection.
Based on vaccination coverage data provided by the Ministry of Health, Consumer Affairs and Social Welfare of Spain (https://www.mscbs.gob.es/profesionales/saludPublica/prevPromocion/vacunaciones/calendario-y-coberturas/coberturas/home.htm, accesed on 14 May 2021) and the clinical trial results presented in [9,10], the following assumptions are incorporated into the simulation process:
In the case of maternal vaccination:
  • A range of vaccine uptake from 50 to 85%;
  • 40% efficacy in reducing RSV bronchiolitis episodes;
  • Duration of protection of 6 months.
With respect to the use of mAb:
  • The uptake parameter can be considered to be 95%;
  • 70% efficacy in reducing RSV bronchiolitis episodes;
  • 6 months of induced immunity.
Because the duration of protection is assumed to be 6 months, children are immune from birth to 6 months, when they move to the next age group. Hence, in our modeling scheme, the evaluation of these two strategies come down to a decrease in the number of cases in the first age group (0–5 months) throughout the simulation study. Since we analyze total counts of bronchiolitis, it is also important to take into account that approximately 70% of bronchiolitis cases are due to RSV.
Table 1 shows the percentage of children in the age group 0–5 months that should be removed according to the specified values for the uptake and the efficacy parameters. These values are also in accordance with the assumptions made in some of the studies discussed in [27].
Let π rsv represent the proportion of bronchiolitis cases due to RSV and π im the proportion of children in the age group 0–5 months that are immune as a result of the implementation of an immunization strategy. For age group 1 (0–5 months), the number of infected children at week t can be estimated from our model as:
y ^ t 1 = ( 1 π im ) · π rsv · exp ( λ t 1 ) · i = 1 4 y ^ t 1 i h + ( 1 π rsv ) · exp ( λ t 1 ) · i = 1 4 y ^ t 1 i h ;
that is, a proportion π im is removed from the number of RSV bronchiolitis cases that would be expected from contacts with infected children at week t 1 . Note that, under the immunization scenario, contacts of immune children will not end in contagion.
For the age groups j = 2 , 3 , 4 , infected children at week t can be estimated as:
y ^ t j = exp ( λ t j ) · i = 1 4 y ^ t 1 i h .
In the previous equations, exp ( λ t j ) · i = 1 4 y ^ t 1 i h represents the mean μ t j of the distribution governing the counts, which is either the Poisson or the negative binomial distribution, based on the simulated counts under the immunization scenario. Parameters { exp ( λ t j ) } and h are, respectively, the autoregressive parameters and the mixing parameter of the model without a immunization strategy (see Equation (6)). It is important to emphasize that these parameters represent features of bronchiolitis dynamics (how the disease spreads at each time point for the different age groups) and contact patterns that do not depend on the number of infected children in previous weeks. Hence, they can be used to simulate the number of infections after a particular preventive intervention has been implemented.

4. Results

R Statistical Software (Foundation for Statistical Computing, Vienna, Austria) and WinBUGS program [34] were used to perform the analysis using MCMC simulation methods. A total of 25,000 iterations were used as the burn-in period of the MCMC. After that, 75,000 iterations were run, and only 1 in every 150 of them was kept to reduce correlation. Two chains were simulated, so M = 1000 values were simulated in total from the posterior distribution. MCMC convergence was assessed by visual inspection of the trace plots, the Brooks–Gelman–Rubin scale reduction factor (Rhat equal to 1 means good convergence) and the effective sample size (n.eff above 100 means good convergence). All statistical analyses of the study are completely reproducible. The BUGS code used can be found as supplementary material to the paper.

4.1. Results without Immunization Strategies

We fitted both the Poisson and the binomial negative model to the bronchiolitis data. The corresponding DIC [35] values are, respectively, 3811.77 (pD 105.89) and 3710.97 (pD 83.96), so adjusting for overdispersion provides a better fit. The posterior mean of parameter k is 44.45, which is the 95% credible interval [32.65,59.41]. Figure 2 displays the posterior mean bronchiolitis temporal profile (red line) together with real counts of disease (black line). As can be seen, the proposed model is able to accurately recover disease dynamics in all the age groups.
Figure 3 shows the estimated autoregressive parameter exp λ t j in each age group and time point together with its seasonal component, which is given by:
exp α 0 j + β · sin 2 · π · t 52 +   γ · cos 2 · π · t 52 .
As expected, previous counts of disease have a higher impact on the youngest children. In fact, the transmission of disease decreases as the age increases. This figure also demonstrates that seasonality has a strong influence on transmission dynamics. However, because the seasonal pattern varies slightly from year to year, the incorporation of random effects accounting for stochasticity in the transmission is fundamental to providing a more accurate description of the data. These results agree with those obtained in [24]. In that paper, the proposed model with a stochastic transmission rate led to an improved goodness of fit in comparison with its counterpart where the seasonal pattern repeated over time; that is, there were no weekly random effects in the model assumed for the transmission rate.

Model Comparison

In order to assess the performance of our model, we also analyzed the bronchiolitis data with a Poisson model whose mean is given by the sum of an endemic and an epidemic components. Taking into account the assumptions made in this particular case study (common seasonality for all ages and contacts among the different age groups that are equally likely to occur), the Poisson model with two components is formulated as:
(12) y t j P o i s s o n ν t j + λ j · i = 1 4 y t 1 i (13) ν t j = exp α 0 j + β · sin 2 · π · t 52 +   γ · cos 2 · π · t 52 .
To ensure that the mean of the Poisson is non-negative, parameter λ j is modeled as λ j = exp ( δ j ) , with δ j N 0 , σ δ .
Figure 4 displays the posterior mean bronchiolitis temporal profile (red line) corresponding to this Poisson model with two components together with real counts of disease (black line). The DIC value for this model is 4506.18 (pD 9.88).
This two-component model is able to describe the overall dynamics of bronchiolitis. However, the Poisson model proposed here provides a better fit, as judged by a lower DIC (3811.77 against 4506.18). This may be due to the fact that the epidemic component, which depends on a time-constant autoregressive parameter, cannot explain the stochasticity inherent in transmission dynamics.

4.2. Results with Immunization Strategies

Once we have obtained a sample from the posterior distribution of the parameters of our model without a immunization program, we can evaluate the effect of the immunization strategies previously described. We assume here that the immunization strategy started six months before the first week of July 2010 (week 28 June to 4 July), so that the specified percentage (see Table 1) of children in the age group of 0–5 months can be assumed to be immune at time t = 1 of our analysis. We simulate the evolution of bronchiolitis counts using Equations (10) and (11), assuming π rsv = 0.7 . Note that for each iteration of the MCMC simulation, we have one simulated value of the regression parameters and the mixing parameter { exp ( λ t j ( m ) ) } and h ( m ) , m = 1 , 2 , , M = 1000 , so we have M simulated values for the number of infections for each age group and week. Using these values, we can get point estimates (for instance, the mean of the simulated values) as well as credible intervals.
Figure 5 shows real counts of bronchiolitis and point estimates (mean values) for three simulation scenarios, corresponding to π im = 0.2 (green line), 0.34 (blue line) and 0.67 (red line). As expected, the number of bronchiolitis cases decreases as the proportion of immune children in the age group of 0–5 month increases. Because of herd protection, the decrease can be observed in all the age groups. The decrease expected in each age group and the total decrease are shown in Table 2.
We also developed an R Shiny application (https://rotapp.shinyapps.io/APP_IMPACTO_RSV, accessed on 14 May 2021) that simulates the evolution of bronchiolitis that would be obtained by varying the uptake and the efficacy parameters, as well as the percentage of bronchiolitis cases due to RSV.

5. Conclusions

We have developed an age-structured stochastic model that allows us to accurately explain bronchiolitis dynamics in the region of Valencia. Our modeling approach combines ideas from compartmental models and Bayesian hierarchical Poisson models in a novel way. By using a hierarchical Poisson (or a binomial negative) model, we do not require information on susceptibles and we avoid the modeling of complex transitions between the multiple compartments. Unlike standard formulations, the main innovation of our model is that the mean of the distribution at each time point depends only on an observation-driven component where the autoregressive parameter is allowed to vary stochastically over time. A mixing parameter is also introduced to capture heterogeneous mixing of individuals in the population. By modeling the dynamics of disease as a function of previous counts, we can simulate disease incidence under different immunization scenarios. Note that the incorporation of a parameter-driven component would hinder simulation of disease evolution when a decrease in the number of cases at previous weeks has been obtained as a result of the assumed immunity.
Even though we have developed the model for the analysis of bronchiolitis, it can be adapted for other infectious diseases with (or without) a seasonal pattern, replacing the transmission rate according to the nature of disease.
The extension proposed in this paper provides a useful framework to address one of the important needs in RSV bronchiolitis incidence control: the implementation of an immunization strategy. We have assessed here the effects of passive immunization through maternal vaccination or mAb administered at birth. In our modelling scheme, the evaluation of these two strategies come down to a decrease in the number of cases in the first age group (0–5 months) throughout the simulation study. We have simulated the evolution of bronchiolitis counts for different values of the uptake and efficacy. The duration of immunity has been assumed to be equal to 6 months. If immunity lasts more than 6 months, the model can be adapted so that a percentage of children in the age group of 6–11 months is also removed, since children would move to this age group being immune for a certain number of weeks. On the other hand, if immunity lasts less than 6 months, the percentage of children removed in the age group 0–5 months should be reduced accordingly.
The results obtained show that newborn immunization contributes to a substantial decrease in the number of bronchiolitis infections. Because of herd immunity, this decrease is also observed in all the age groups considered. As expected, the decrease is greater when the uptake and/or the efficacy increase. As the profiles of the products under development become better defined, further alignment will be possible.
We have developed an R Shiny application that simulates the evolution of bronchiolitis episodes for different values of the uptake and efficacy, as well as the percentage of bronchiolitis cases due to RSV. This application is an easy tool that allows users and decision makers to interact with the simulation analysis and improves visualization of the results.
We have not assessed here the impact of immunization of high-risk children, such as preterm infants or children with high-risk comorbid conditions. It may be important to consider a separate analysis for these groups most at risk.
Another fruitful area for further research would be the consideration of smaller age groups so that different immunization strategies such as vaccination of infants older than 3 months could also be tested. Stratification by gestational age may also be important when evaluating maternal vaccination strategies, since transfer of maternal antibodies for preterm infants may be incomplete.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/ijerph18147607/s1.

Author Contributions

Data Curation: M.L.-L. Investigation: M.L.-L., A.C.-V. and F.J.S.G. Methodology: M.L.-L., A.C.-V. and F.J.S.G. Formal Analysis: M.L.-L. and A.C.-V. Software: M.L.-L. and A.C.-V. Writing—Original Draft Preparation: M.L.-L., A.C.-V. and F.J.S.G. All authors have read and agreed to the published version of the manuscript.

Funding

This work has been supported by Grant Number MTM2017-83850-P from the Spanish Ministry of Economy and Competitiveness.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors are grateful to Alejandro Orrico Sánchez, Cintia Muñoz Quiles and Javier Díez Domingo, from the Vaccine Research Department, Fisabio-Public Health, for providing practical guidance and insightful comments on a previous version of this manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Korsten, K.; Bont, L. Seasonal immunisation against respiratory syncytial virus disease. Lancet Public Health 2017, 2, e344–e345. [Google Scholar] [CrossRef]
  2. Ralston, S.L.; Lieberthal, A.S.; Meissner, H.C.; Alverson, B.K.; Baley, J.E.; Gadomski, A.M.; Johnson, D.W.; Light, M.J.; Maraqa, N.F.; Mendonca, E.A.; et al. Clinical Practice Guideline: The Diagnosis, Management, and Prevention of Bronchiolitis. Pediatria 2014, 134, e1474–e1502. [Google Scholar] [CrossRef] [Green Version]
  3. Díez-Domingo, J.; Ridao-López, M.; Úbeda-Sansano, I.; Ballester-Sanz, A. Incidence and Cost of Hospitalizations for Bronchiolitis and Respiratory Syncytial Virus Infections in the Autonomous Community of Valencia in Spain (2001 and 2002). Anal. Pediatr. 2006, 65, 325–330. [Google Scholar] [CrossRef] [PubMed]
  4. Muñoz-Quiles, C.; López-Lacort, M.; Úbeda-Sansano, I.; Alemán-Sánchez, S.; Perez-Vilar, S.; Puig-Barberà, J.; Diez-Domingo, J. Population-based Analysis of Bronchiolitis Epidemiology in Valencia, Spain. Pediatr. Infect. Dis. J. 2016, 35, 275–280. [Google Scholar] [CrossRef] [PubMed]
  5. Resch, B. Respiratory Syncytial Virus Infection in High-risk Infants—An Update on Palivizumab Prophylaxis. Open Microbiol. J. 2014, 8, 71–77. [Google Scholar] [CrossRef] [Green Version]
  6. Shook, B.C.; Lin, K. Recent Advances in Developing Antiviral Therapies for Respiratory Syncytial Virus. Top. Curr. Chem. 2017, 375, 40. [Google Scholar] [CrossRef] [PubMed]
  7. Rainisch, G.; Adhikari, B.; Meltzer, M.I.; Langley, G. Estimating the impact of multiple immunization products on medically-attended respiratory syncytial virus (RSV) infections in infants. Vaccine 2020, 38, 251–257. [Google Scholar] [CrossRef] [PubMed]
  8. Domachowske, J.B.; Anderson, E.J.; Goldstein, M. The Future of Respiratory Syncytial Virus Disease Prevention and Treatment. Infect. Dis. Ther. 2021, 10, 47–60. [Google Scholar] [CrossRef] [PubMed]
  9. Madhi, S.A.; Polack, F.P.; Piedra, P.A.; Munoz, F.M.; Trenholme, A.A.; Simões, E.A.; Swamy, G.K.; Agrawal, S.; Ahmed, K.; August, A.; et al. Respiratory Syncytial Virus Vaccination during Pregnancy and Effects in Infants. N. Engl. J. Med. 2020, 383, 426–439. [Google Scholar] [CrossRef]
  10. Griffin, M.P.; Yuan, Y.; Takas, T.; Domachowske, J.B.; Madhi, S.A.; Manzoni, P.; Simões, E.A.; Esser, M.T.; Khan, A.A.; Dubovsky, F.; et al. Single-Dose Nirsevimab for Prevention of RSV in Preterm Infants. N. Engl. J. Med. 2020, 383, 415–425. [Google Scholar] [CrossRef] [PubMed]
  11. Lanari, M.; Giovannini, M.; Giuffré, L.; Marini, A.; Rondini, G.; Rossi, G.A.; Merolla, R.; Zuccotti, G.V.; Salvioli, G. Prevalence of respiratory syncytial virus infection in Italian infants hospitalized for acute lower respiratory tract infections, and association between respiratory syncytial virus infection risk factors and disease severity. Pediatr. Pulmonol. 2002, 33, 458–465. [Google Scholar] [CrossRef]
  12. Koehoorn, M.; Karr, C.J.; Demers, P.A.; Lencar, C.; Tamburic, L.; Brauer, M. Descriptive Epidemiological Features of Bronchiolitis in a Population-Based Cohort. Pediatria 2008, 122, 1196–1203. [Google Scholar] [CrossRef] [PubMed]
  13. García, C.G.; Bhore, R.; Soriano-Fallas, A.; Trost, M.; Chason, R.; Ramilo, O.; Mejias, A. Risk Factors in Children Hospitalized with RSV Bronchiolitis Versus Non-RSV Bronchiolitis. Pediatria 2010, 126, e1453–e1460. [Google Scholar] [CrossRef] [Green Version]
  14. Flaherman, V.J.; I Ragins, A.; Li, S.X.; Kipnis, P.; Masaquel, A.; Escobar, G.J. Frequency, duration and predictors of bronchiolitis episodes of care among infants ≥32 weeks gestation in a large integrated healthcare system: A retrospective cohort study. BMC Health Serv. Res. 2012, 12, 144. [Google Scholar] [CrossRef] [Green Version]
  15. Murray, J.; Bottle, A.; Sharland, M.; Modi, N.; Aylin, P.; Majeed, A.; Saxena, S.; on behalf of the Medicines for Neonates Investigator Group. Risk Factors for Hospital Admission with RSV Bronchiolitis in England: A Population-Based Birth Cohort Study. PLoS ONE 2014, 9, e89186. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Weber, A.; Weber, M.; Milligan, P. Modeling Epidemics Caused by Respiratory Sincytial Virus (RSV). Math. Biosci. 2001, 172, 95–113. [Google Scholar] [CrossRef]
  17. Acedo, L.; Diez-Domingo, J.; Moraño, J.-A.; Villanueva, R.-J. Mathematical modelling of respiratory syncytial virus (RSV): Vaccination strategies and budget applications. Epidemiol. Infect. 2009, 138, 853–860. [Google Scholar] [CrossRef]
  18. Leecaster, M.; Gesteland, P.; Greene, T.; Walton, N.; Gundlapalli, A.; Rolfs, R.; Byington, C.; Samore, M. Modeling the variations in pediatric respiratory syncytial virus seasonal epidemics. BMC Infect. Dis. 2011, 11, 105. [Google Scholar] [CrossRef] [Green Version]
  19. Moore, H.C.; Jacoby, P.; Hogan, A.; Blyth, C.; Mercer, G.N. Modelling the Seasonal Epidemics of Respiratory Syncytial Virus in Young Children. PLoS ONE 2014, 9, e100422. [Google Scholar] [CrossRef] [Green Version]
  20. Isham, V. Stochastic models for epidemics. In Celebrating Statistics: Papers in Honour of Sir David Cox on the Occasion of his 80th Birthday; Cox, D.R., Davison, A.C., Yadolah, D., Wermuth, N., Eds.; Oxford University Press: Oxford, UK, 2005; pp. 27–53. [Google Scholar]
  21. Allen, L.J.S. An Introduction to Stochastic Epidemic Models. In Mathematical Epidemiology. Lecture Notes in Mathematics; Brauer, F., van den Driessche, P., Wu, J., Eds.; Springer: Berlin, Germany, 2008; pp. 81–130. [Google Scholar] [CrossRef]
  22. Brauer, F. Compartmental Models in Epidemiology. In Mathematical Epidemiology. Lecture Notes in Mathematics; Brauer, F., van den Driessche, P., Wu, J., Eds.; Springer: Berlin, Germany, 2008; pp. 19–79. [Google Scholar] [CrossRef]
  23. Greenwood, P.E.; Gordillo, L.F. Stochastic Epidemic Modeling. In Mathematical and Statistical Estimation Approaches in Epidemiology; Chowell, G., Hyman, J.M., Bettencourt, L.M.A., Castillo-Chavez, C., Eds.; Springer: Berlin, Germany, 2009; pp. 31–52. [Google Scholar] [CrossRef]
  24. Corberán-Vallet, A.; Santonja, F.J. A Bayesian SIRS model for the analysis of respiratory syncytial virus in the region of Valencia, Spain. Biom. J. 2014, 56, 808–818. [Google Scholar] [CrossRef]
  25. Held, L.; Hofmann, M.; Höhle, M.; Schmid, V. A Two-Component Model for Counts of Infectious Diseases. Biostatistics 2006, 7, 422–437. [Google Scholar] [CrossRef] [Green Version]
  26. Held, L.; Höhle, M.; Hofmann, M. A statistical framework for the analysis of multivariate infectious disease surveillance counts. Stat. Model. Int. J. 2005, 5, 187–199. [Google Scholar] [CrossRef] [Green Version]
  27. Treskova, M.; Pozo-Martin, F.; Scholz, S.; Schönfeld, V.; Wichmann, O.; Harder, T. Assessment of the Effects of Active Immunisation against Respiratory Syncytial Virus (RSV) using Decision-Analytic Models: A Systematic Review with a Focus on Vaccination Strategies, Modelling Methods and Input Data. PharmacoEconomics 2021, 39, 287–315. [Google Scholar] [CrossRef] [PubMed]
  28. Mezei, A.; Cohen, J.; Renwick, M.J.; Atwell, J.; Portnoy, A. Mathematical modelling of respiratory syncytial virus (RSV) in low- and middle-income countries: A systematic review. Epidemics 2021, 35, 100444. [Google Scholar] [CrossRef] [PubMed]
  29. Hogan, A.B.; Campbell, P.; Blyth, C.; Lim, F.; Fathima, P.; Davis, S.; Moore, H.C.; Glass, K. Potential impact of a maternal vaccine for RSV: A mathematical modelling study. Vaccine 2017, 35, 6172–6179. [Google Scholar] [CrossRef]
  30. Hodgson, D.; Pebody, R.; Panovska-Griffiths, J.; Baguelin, M.; Atkins, K.E. Evaluating the next generation of RSV intervention strategies: A mathematical modelling study and cost-effectiveness analysis. BMC Med. 2020, 18, 348. [Google Scholar] [CrossRef] [PubMed]
  31. Jornet-Sanz, M.; Corberán-Vallet, A.; Santonja, F.J.; Villanueva, R.J. A Bayesian Stochastic SIRS Model with a Vaccination Strategy for the Analysis of Respiratory Syncytial Virus. Sort Stat. Oper. Res. Trans. 2017, 41, 159–175. [Google Scholar] [CrossRef]
  32. García-Sempere, A.; Sanchez, A.O.; Muñoz-Quiles, C.; Hurtado, I.; Peiró, S.; Sanfélix-Gimeno, G.; Diez-Domingo, J. Data Resource Profile: The Valencia Health System Integrated Database (VID). Int. J. Epidemiol. 2020, 49, 740–741e. [Google Scholar] [CrossRef] [Green Version]
  33. Paul, M.; Held, L.; Toschke, A.M. Multivariate modelling of infectious disease surveillance data. Stat. Med. 2008, 27, 6250–6267. [Google Scholar] [CrossRef]
  34. Lunn, D.J.; Thomas, A.; Best, N.; Spiegelhalter, D. WinBUGS—A Bayesian modelling framework: Concepts, structure, and extensibility. Stat. Comput. 2000, 10, 325–337. [Google Scholar] [CrossRef]
  35. Spiegelhalter, D.J.; Best, N.G.; Carlin, B.P.; Van Der Linde, A. Bayesian measures of model complexity and fit. J. R. Stat. Soc. Ser. B Stat. Methodol. 2002, 64, 583–639. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Weekly counts of bronchiolitis for the four age groups from July 2010 (week 28 June to 4 July) to December 2012 (week 24 December to 30 December).
Figure 1. Weekly counts of bronchiolitis for the four age groups from July 2010 (week 28 June to 4 July) to December 2012 (week 24 December to 30 December).
Ijerph 18 07607 g001
Figure 2. Model accuracy. Bronchiolitis counts (black line) together with posterior mean estimates (red line) from July 2010 to December 2012 (T = 131 weeks) obtained with the proposed model based on the negative binomial distribution.
Figure 2. Model accuracy. Bronchiolitis counts (black line) together with posterior mean estimates (red line) from July 2010 to December 2012 (T = 131 weeks) obtained with the proposed model based on the negative binomial distribution.
Ijerph 18 07607 g002
Figure 3. Estimated autoregressive parameter (red line) together with its seasonal component (black line).
Figure 3. Estimated autoregressive parameter (red line) together with its seasonal component (black line).
Ijerph 18 07607 g003
Figure 4. Bronchiolitis counts (black line) together with posterior mean estimates (red line) obtained with a Poisson model with two components.
Figure 4. Bronchiolitis counts (black line) together with posterior mean estimates (red line) obtained with a Poisson model with two components.
Ijerph 18 07607 g004
Figure 5. Bronchiolitis counts (black line) together with simulated counts (mean values) for three different immunization proportions: π im = 0.2 (green line), 0.34 (blue line) and 0.67 (red line).
Figure 5. Bronchiolitis counts (black line) together with simulated counts (mean values) for three different immunization proportions: π im = 0.2 (green line), 0.34 (blue line) and 0.67 (red line).
Ijerph 18 07607 g005
Table 1. Possible scenarios in the evaluation of immunization strategies.
Table 1. Possible scenarios in the evaluation of immunization strategies.
Maternal VaccinationmAb at Birth
UptakeEfficacy% RemovedUptakeEfficacy% Removed
50%40%20%95%70%66.50%
85%40%34%
Table 2. Expected decrease in the number of bronchiolitis counts for three different immunization proportions: π im = 0.2, 0.34 and 0.67.
Table 2. Expected decrease in the number of bronchiolitis counts for three different immunization proportions: π im = 0.2, 0.34 and 0.67.
π im 0–5 Months6–11 Months12–17 Months18–23 MonthsOverall
0.228.94%14.63%15.48%16.25%21.68%
0.3447.48%28.75%29.78%30.03%37.92%
0.6776.32%54.17%54.56%54.84%64.87%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

López-Lacort, M.; Corberán-Vallet, A.; Santonja Gómez, F.J. A Multivariate Age-Structured Stochastic Model with Immunization Strategies to Describe Bronchiolitis Dynamics. Int. J. Environ. Res. Public Health 2021, 18, 7607. https://doi.org/10.3390/ijerph18147607

AMA Style

López-Lacort M, Corberán-Vallet A, Santonja Gómez FJ. A Multivariate Age-Structured Stochastic Model with Immunization Strategies to Describe Bronchiolitis Dynamics. International Journal of Environmental Research and Public Health. 2021; 18(14):7607. https://doi.org/10.3390/ijerph18147607

Chicago/Turabian Style

López-Lacort, Mónica, Ana Corberán-Vallet, and Francisco J. Santonja Gómez. 2021. "A Multivariate Age-Structured Stochastic Model with Immunization Strategies to Describe Bronchiolitis Dynamics" International Journal of Environmental Research and Public Health 18, no. 14: 7607. https://doi.org/10.3390/ijerph18147607

APA Style

López-Lacort, M., Corberán-Vallet, A., & Santonja Gómez, F. J. (2021). A Multivariate Age-Structured Stochastic Model with Immunization Strategies to Describe Bronchiolitis Dynamics. International Journal of Environmental Research and Public Health, 18(14), 7607. https://doi.org/10.3390/ijerph18147607

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop