bioRxiv preprint doi: https://doi.org/10.1101/665646; this version posted June 22, 2019. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
1 Vaccination strategies to control Ebola epidemics in the context of variable
2 household inaccessibility levels
3
4 G. Chowell*1, A. Tariq1, M. Kiskowski2 5 6 1Department of Population Health Sciences, School of Public Health, Georgia State
7 University, 30303, Atlanta, GA, USA 8 2Department of Mathematics and Statistics, University South Alabama, 36688, Mobile, 9 AL, USA 10 11 *Corresponding author 12 13 14 15
1
bioRxiv preprint doi: https://doi.org/10.1101/665646; this version posted June 22, 2019. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
16 Abstract 17 Despite a very effective vaccine, active conflict and community distrust during the 18 ongoing DRC Ebola epidemic are undermining control efforts, including a ring 19 vaccination strategy that require the prompt immunization of close contacts of 20 infected individuals. However, 20% or more of close contacts cannot be reached or 21 refuse vaccination [1], and it is predicted that the ring vaccination strategy would not 22 be effective with such a high level of inaccessibility. In May 2019, the DRC expanded 23 their control efforts to include a community vaccine strategy. To assess the impact of 24 vaccination strategies for controlling Ebola epidemics in the context of variable levels 25 of community resistance, we employed an individual-level stochastic transmission 26 model that incorporates four sources of heterogeneity: a proportion of the population 27 is inaccessible for contact tracing and vaccination due to lack of confidence in 28 interventions or geographic inaccessibility, two levels of population mixing resembling 29 household and community transmission, two types of vaccine doses with different 30 time periods until immunity, and transmission rates that depend on spatial distance. 31 Our results indicate that a ring vaccination strategy alone would not be effective for 32 containing the epidemic in the context of significant delays to vaccinating contacts 33 even for low levels of household inaccessibility and affirm the positive impact of a 34 supplemental community vaccination strategy. Our key results are that as levels of 35 inaccessibility increase, there is a qualitative change in the effectiveness of the 36 vaccination strategy. For higher levels of vaccine access, the probability that the 37 epidemic will end steadily increases over time, even if probabilities are lower than they 38 would be otherwise with full community participation. For levels of vaccine access that 39 are too low, however, the vaccination strategies are not expected to be successful in 40 ending the epidemic even though they help lower incidence levels, which saves lives, 41 makes the epidemic easier to contain and reduces spread to other communities. This 42 qualitative change occurs for both types of vaccination strategies: ring vaccination is 43 effective for containing an outbreak until the levels of inaccessibility exceeds 44 approximately 10% in the context of significant delays to vaccinating contacts, a 45 combined ring and community vaccination strategy is effective until the levels of 46 inaccessibility exceeds approximately 50%. More broadly, our results underscore the 47 need to enhance community engagement to public health interventions in order to 48 enhance the effectiveness of control interventions to ensure outbreak containment. 49 50 51
2
bioRxiv preprint doi: https://doi.org/10.1101/665646; this version posted June 22, 2019. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
52 Keywords: Mathematical modeling, individual-based model, Ebola, vaccination 53 strategies, ring vaccination, community vaccination, community accessibility, 54 community resistance, epidemic containment, endemic state, epidemic wave. 55 56 57 Author summary 58 59 Containing Ebola outbreaks hinges on educating local communities. In the context of 60 the ongoing Ebola epidemic in DRC, active conflict and community distrust are 61 undermining control efforts, including vaccination strategies. In this paper, we 62 employed an individual-level stochastic structured transmission model to assess the 63 impact of vaccination strategies on epidemic control in the context of variable levels of 64 household inaccessibility. We found that a ring vaccination strategy would not be 65 effective for containing the epidemic in the context of significant delays to vaccinating 66 contacts even for low levels of household inaccessibility and evaluate the impact of a 67 supplemental community vaccination strategy. For lower levels of inaccessibility, the 68 probability of epidemic containment increases over time. For higher levels of 69 inaccessibility, vaccination strategies are not expected to contain the epidemic even 70 though they help lower incidence levels, which saves lives, makes the epidemic easier 71 to contain and reduces spread to other communities. We found that ring vaccination is 72 effective for containing an outbreak until the levels of inaccessibility exceeds 73 approximately 10%, a combined ring and community vaccination strategy is effective 74 until the levels of inaccessibility exceeds approximately 50%. Our findings underscore 75 the need to enhance community engagement to public health interventions. 76 77 78
3
bioRxiv preprint doi: https://doi.org/10.1101/665646; this version posted June 22, 2019. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
79 Introduction
80 A key lesson from past Ebola outbreaks indicates that their containment hinges on
81 building community trust of response teams and educating local communities in order
82 to fight misinformation relating to the transmission of the virus in the population and
83 the role of control interventions [2, 3]. These educational activities are critical for the
84 successful deployment of contact tracing and surveillance efforts, treatment and
85 management of Ebola cases, and promoting behaviors that mitigate transmission rates
86 [4]. However, the ongoing Ebola outbreak in the Democratic Republic of Congo (DRC),
87 which has been active in the region for about 10 months (August 2018 to May 2019),
88 has become the most complex Ebola epidemic to date, threatening to spread to
89 neighboring countries. In particular, this unprecedented epidemic is unfolding in a
90 climate of community distrust in the interventions [5, 6] in an active conflict zone [7],
91 where over one hundred attacks have deliberately targeted healthcare workers and
92 treatment centers involved in the Ebola response efforts, allowing sustained
93 transmission in the region. Consequently, the DRC Ebola epidemic has accumulated
94 2084 cases including 1405 deaths as of June 11, 2019 and has spilled over to Uganda
95 where 3 cases have been confirmed as of June 12, 2019 [8, 9]. A recent surge in Ebola
96 incidence coincided with an increasing trend in the number of violent attacks on health
97 centers and health teams fighting the epidemic in the region [10].
98
99 Although control efforts now employ a highly effective emergency vaccine, the
100 ongoing Ebola epidemic in the DRC is the first to occur in an active conflict zone.
101 Deliberate violent attacks and threats to health workers on a scale not seen in previous
102 Ebola outbreaks are directly undermining active control efforts in the region [7]. A ring
4
bioRxiv preprint doi: https://doi.org/10.1101/665646; this version posted June 22, 2019. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
103 vaccination strategy [11, 12] that was highly effective in the capital [13] was
104 introduced rapidly on August 8th, 2018 within 8 days of the outbreak declaration in the
105 provinces of North Kivu and Ituri [14]. However, it has been challenging to implement
106 an effective ring vaccination strategy in a climate of violence targeting healthcare
107 workers and decreasing accessibility to a mobile, fractionated population as it relies on
108 the identification of contacts and contacts of contacts [15]. This issue is compounded
109 by a large fraction of Ebola cases unconnected to known chains of transmission [16]. In
110 order to boost the impact of the vaccine in a challenging transmission setting with a
111 large fraction of hidden or underground transmission, the vaccination strategy is now
112 incorporating community-level vaccination. The community-level vaccination strategy
113 will immunize members of communities even if they are not known contacts of Ebola
114 cases using a lower vaccine dose, which will take longer to elicit protective antibody
115 levels in vaccinated individuals [17, 18].
116
117 Here we sought to evaluate the impact of ongoing Ebola vaccination campaigns by
118 extending an individual-based model with household-community mixing [19, 20],
119 which has been previously used to analyze the effects of lower versus higher rates of
120 population mixing on Ebola transmission dynamics. This model has been able to fit the
121 growth patterns of the 2014-15 Western Africa Ebola epidemic (e.g., the growth
122 pattern of the epidemic in Guinea, versus Sierra Leone or Liberia) by calibrating the
123 extent of community mixing [19]. This structured transmission model predicts
124 outbreaks that propagate through the population as spatial waves with an endemic
125 state [20] and has been useful to gain insight on the level of control that would be
126 required to contain Ebola epidemics [20]. One of the less intuitive results of the
5
bioRxiv preprint doi: https://doi.org/10.1101/665646; this version posted June 22, 2019. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
127 community model was that even a low daily incidence can indicate an epidemic that is
128 difficult to extinguish if saturation effects decrease incidence and are masking a more
129 robust transmission rate. Saturation effects occur when contacts of infectious
130 individuals are already infected by other members of the community, such as family
131 members or members of other close groups, and decrease the incidence when
132 community mixing is low.
133
134 Assessing the effect of targeted vaccination efforts requires mathematical models that
135 capture the contact structure of the community network [21-25]. In our paper, we
136 aimed to investigate the effect of ring [21, 23, 26-28] and community vaccination
137 strategies on outbreak control. For this purpose, our baseline stochastic model, which
138 has been calibrated using data for the Western African Ebola epidemic [19, 20], was
139 adapted to incorporate key features of the Ebola epidemic in the DRC. Specifically, we
140 evaluate the effectiveness of control strategies in the context of a fraction of the
141 population that is accessible to vaccination teams, tied to the success of contact
142 tracing and vaccination efforts. Our baseline network transmission model also adapts
143 for the ring vaccination strategy by integrating heterogeneity in community
144 transmission rates that scale with the distance between an infectious individual and
145 each member of that individual’s community. As the distance between an infectious
146 individual and a contact increases, the transmission rate decreases exponentially (or by
147 any other function). Ring and community vaccination strategies are compared by
148 accounting for the two different spatial scales for the radius of contacts that are
149 vaccinated (either as a ring or community-wide) and the two different vaccine doses (a
6
bioRxiv preprint doi: https://doi.org/10.1101/665646; this version posted June 22, 2019. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
150 full or half dose), that result in different time periods for the immunization of
151 vaccinated contacts.
152
153 Methods
154
155 Household-Community Model for Ebola Spread and Two Types of Vaccination with
156 Accessible and Inaccessible Households
157 We previously described our model for disease spread [19, 20] and vaccination [29] in
158 a contact network with household-community structure. A new component of our
159 model is that we add a second type of household that is inaccessible to vaccination
160 teams. These households may not participate for a variety of reasons including Ebola
161 cases that are not identified, lack of confidence in the vaccination program and
162 geographic inaccessibility. The inaccessible households do not participate in
163 vaccination and do not provide contacts lists. We also modify the model to include two
164 vaccine types. A larger vaccine dose is used for ring vaccination with faster
165 immunization occurring 10 days later, and a smaller vaccine dose is used for
166 community vaccination with immunization occurring 28 days later. This corresponds to
167 the dosing and predicted timing of the ring and community vaccination regimens of 0.5
168 ml and 0.2 ml doses, respectively, recently described for the DRC Ebola epidemic [15,
169 17].
170
171 Disease and Vaccination States
172 The progression of Ebola disease, vaccination and immunization are modeled with
173 seven epidemiological states (S, E, I, R, Svr, Svc, M); including four SEIR states for Ebola
7
bioRxiv preprint doi: https://doi.org/10.1101/665646; this version posted June 22, 2019. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
174 disease progression (susceptible S, exposed E, infectious I and refractory R) and three
175 states for vaccination progression (two susceptible but vaccinated states SVR and SVC
176 and an immune state M) (see Figure 1). The states Svc and Svr are assigned to
177 individuals vaccinated by community or ring vaccination. Individuals with these states
178 remain susceptible to Ebola exposure. After a delay of 10 days for Svr states or 28 days
179 for Svc states, the individuals become immune by transitioning to state M.
180 Individuals in three epidemiological states (S, Svc and Svr) are at risk of contracting
181 Ebola from infectious individuals (I), and transitioning to the latent period (exposed
182 state E). Transition rates from any of these three susceptible states to the exposed
183 state depend on the contact network and increase with the number of infectious
184 contacts. In our model, individuals in the exposed state E (latency) are no longer
185 available for effective vaccination and immunization. Thus, if an individual is indicated
186 to be vaccinated by either the ring or community vaccination programs, they may be
187 approached by a vaccination team and given the vaccine since their exposure status is
188 unknown, but once in the exposed class, they may not change states to the susceptible
189 states Svc or Svr. Once exposed, nodes transition from the exposed (E) to infectious
190 class (I) with probability per day where 9 days is the average incubation period
191 [30-34] and once infectious, nodes transition from the infectious (I) to refractory class
192 (R) with probability per day, where 5.6 days is the average infectious period [30-
193 33]. Parameter definitions and values used in simulations and their sources are given
194 in Table 1.
195
196
8
bioRxiv preprint doi: https://doi.org/10.1101/665646; this version posted June 22, 2019. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
197 Contact Network with Household and Community Structure
198 As in refs. [19, 20, 29], we model the spread of Ebola on a contact network that
199 consists of households of size H that are organized in communities of size C
200 households. Households are arranged as a continuous line of households within an
201 LxH grid where L is the number of households in the total contact network N.
202 Households are indexed by their position in the line , ,… and a network
203 distance η between two households hi and hj is defined as the number of households
204 separating them (η=|i-j|). Each household has its own “community” that overlaps with
205 the communities of other nearby households. In particular, the ith community is
206 centered at the ith household and includes all of the households within a radius of the
207 ith household (for a community containing C households, the community radius is
208 Rc=(C-1)/2). This contact network is highly mutually connected, the extent of overlap
209 between the ith and jth communities depend on the network distance of the ith and jth
210 households in that communities of nearby households overlap and share most
211 households. Due to this high mutual connectivity and linear arrangement of the
212 contact network, saturation effects build up quickly and disease progresses spatially as
213 1-dimensional wave through the lattice network. The size of the contact network and
214 the size of the simulated population is unlimited since the lattice size LxH is dynamic
215 and extends as needed as the disease and immunity extends through the population
216 by adding households to the ends of the array.
217
218 In this version of the model, each household is designated as accessible or inaccessible
219 as illustrated in Figure 2. This assignment is random when the household is created
220 and once a household is labeled, this assignment does not change. During the initial
9
bioRxiv preprint doi: https://doi.org/10.1101/665646; this version posted June 22, 2019. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
221 construction of the contact network, and as the array is dynamically extended,
222 households are labeled inaccessible with probability β, otherwise the household is
223 accessible. Here we compare results as a function of the fraction of inaccessible
224 households between 0 and 50%. (This range includes the estimate of 30% non-
225 participating individuals that was estimated for a vaccination program when it is
226 expanded to younger age groups in [35]; and indeed the program was approved for
227 expansion to younger age groups in the first half of 2019 [36] This range also includes
228 estimates of 25% for the number of cases that are missed for the ongoing DRC
229 epidemic[37].
230
231 Distance Dependent Transmission Rates
232 At the onset of an epidemic, before the accumulation of saturation effects when all
233 nodes of the network are susceptible except for a single infectious case, the household
234 reproductive number R0H is the average number of secondary infections within the
235 household of the infected individual and the community reproductive number R0C is
236 the average number of secondary infections within the community. As in ref. [20], we
237 set R0H =2.0 infections within the household and R0C =0.7 infections within the
238 community, based on total R0 and how they were observed to distribute within
239 household and communities in historical Ebola epidemics [4, 19, 30, 38].
240
241 Transmission rates for homogenous network transmission: In the simplest case of
242 homogeneous transmission rates, all members of a household have an equal
243 probability of infection from an infectious contact, and all members of the broader
244 community have an equal probability of infection from an infectious contact.
10
bioRxiv preprint doi: https://doi.org/10.1101/665646; this version posted June 22, 2019. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
245 Transmission rates within the household and within the community are normalized so
246 that the reproductive number of household and community infections are always R0H
247 and R0C, respectively. Since there are H-1 household contacts within the household of
248 a single infectious individual, and community contacts within the
249 community of a single infectious individual, for a fixed infectious period λ, household
250 and community transmission rates for each infectious-susceptible contact are given by:
251 tā , tü . ·
252 However, as described in refs. [39, 40], for exponentially distributed infectious periods
253 λ, household and community transmission rates for each infectious-susceptible contact
254 would be: