key: cord-103995-vok1dni9 authors: Chin, Taylor; Buckee, Caroline O.; Mahmud, Ayesha S. title: Quantifying the success of measles vaccination campaigns in the Rohingya refugee camps date: 2019-10-05 journal: nan DOI: 10.1101/19008417 sha: doc_id: 103995 cord_uid: vok1dni9 In the wake of the Rohingya population's mass migration from Myanmar, one of the world's largest refugee settlements was constructed in Cox's Bazar, Bangladesh to accommodate nearly 900,000 new refugees. Refugee populations are particularly vulnerable to infectious disease outbreaks due to many population and environmental factors. A large measles outbreak, with over 2,500 cases, occurred among the Rohingya population between September and December 2017. Here, we estimate key epidemiological parameters and use a dynamic mathematical model of measles transmission to evaluate the effectiveness of the reactive vaccination campaigns in the refugee camps. We also estimate the potential for subsequent outbreaks under different vaccination coverage scenarios. Our modeling results highlight the success of the vaccination campaigns in rapidly curbing transmission and emphasize the public health importance of maintaining high levels of vaccination in this population, where high birth rates and historically low vaccination coverage rates create suitable conditions for future measles outbreaks. 3 . CC-BY-NC-ND 4.0 International license It is made available under a author/funder, who has granted medRxiv a license to display the preprint in perpetuity. is the (which was not peer-reviewed) The copyright holder for this preprint . https://doi.org/10.1101/19008417 doi: medRxiv preprint that occurs in other settings [28] . Figure 1 : Schematic representation of disease states in a model of a measles outbreak in Cox's Bazar. Individuals are either susceptible to measles (S), infectious with measles (I), or recovered from and immune to measles infection (R), either due natural immunity from infection or prior vaccination. β represents the transmission coefficient, γ is the recovery rate, and µ is the birth and death rate. The model starts with one infected (and infectious) measles case (i.e., I(0) = 1). The number of recovered individuals at time = 0 was calculated as R(0) = N − S(0) − I(0). The model tracks the number of individuals that move between the compartments each day and is run for one year using the ordinary differential equations 1-3. The transmission coefficient, β, was estimated according to the relationship β = R0/N D, where N is the total 93 population (i.e., S + I + R) and D is the average duration of infectiousness ( The model uses one infected case as the initial value for the number of individuals in the infected compartment, I(0). The initial value for the recovered compartment was calculated as (Table 1) . The model runs for 365 days to reflect one year following the start of the outbreak. LHS was used to take into account uncertainty in parameter estimates; 1,000 values for each parameter were sampled 105 from the assumed parameter probability distributions (Table 1) where VC is ≥ 1 dose coverage of the Measles, Mumps, and Rubella (MMR) vaccine, VE is the median effectiveness . CC-BY-NC-ND 4.0 International license It is made available under a author/funder, who has granted medRxiv a license to display the preprint in perpetuity. is the (which was not peer-reviewed) The copyright holder for this preprint . https://doi.org/10.1101/19008417 doi: medRxiv preprint now. Table 1 reports all parameters and calculations used as inputs in the mathematical model. LHS was again used 113 to take into account uncertainty in parameter estimates. 114 For this analysis, the SIR model was run for two years to estimate the impact of vaccination on the next outbreak. The cumulative number of cases over two years was calculated for vaccination coverage rates of 10%, 50%, and 90% of 116 the at-risk population. Cumulative cases were estimated as the cumulative proportion of the population that moves 117 into the infectious compartment in the SIR framework. Distributions of our estimates for the cumulative number of 118 cases at two years under the different vaccination rate scenarios are presented. To contrast these results with the 119 scenario of the next measles outbreak in the absence of vaccination, the same model used to estimate the size of the 120 2017 outbreak in the absence of vaccination was used with two changes: 1) the total population size was increased to 121 870,534 [1] , and 2) the proportion of the population that was assumed to be susceptible at the start of the outbreak 122 was estimated as su instead of s0 (Table 1) Figure 3 ). Vaccination therefore has the potential to avert between approximately 4,000 and 19,000 cases over two 160 years depending on the coverage rate. Figure 4B shows estimates in terms of boxplots to reflect the uncertainty in 161 our model's results due to parameter uncertainty. . 167 We find that nearly 77,000 cases were averted over the course of one year due to vaccination in the 2017 outbreak. The daily effective reproductive number declined from its maximum value of 4.3 at the end of week 5 to below 1 by 169 the middle of week 11. When we apply our model to estimate the impact of vaccination on subsequent outbreaks, 170 6 . CC-BY-NC-ND 4.0 International license It is made available under a author/funder, who has granted medRxiv a license to display the preprint in perpetuity. is the (which was not peer-reviewed) The copyright holder for this preprint . https://doi.org/10.1101/19008417 doi: medRxiv preprint CC-BY-NC-ND 4.0 International license It is made available under a author/funder, who has granted medRxiv a license to display the preprint in perpetuity. is the (which was not peer-reviewed) The copyright holder for this preprint . https://doi.org/10.1101/19008417 doi: medRxiv preprint isolation of infectious cases were implemented during the outbreak, which would have the effect of shortening the it may be reasonable to assume that these interventions were not widely used given their infeasibility in this setting. 179 Secondly, the contraction of the serial interval described by Kenah, Lipsitch, and Robins [35] in high transmission 180 settings is not considered in this analysis. Instead, a constant serial interval distribution is assumed. For these reasons, Second, reporting rates and how these rates changed over time are also not considered in this analysis. One 187 possible scenario is that identification of measles was worse at the start of the epidemic relative to later in the 188 epidemic, and that fewer measles cases were therefore reported in the early phase of the epidemic. If this were the 189 case, the initial case counts may underestimate the true extent of measles transmission at the start of the epidemic, CC-BY-NC-ND 4.0 International license It is made available under a author/funder, who has granted medRxiv a license to display the preprint in perpetuity. is the (which was not peer-reviewed) The copyright holder for this preprint . https://doi.org/10.1101/19008417 doi: medRxiv preprint Figure 1 : Histogram of estimated R 0 values using LHS with 1,000 parameter combinations, where the maximum has been restricted to 20. Summary statistics: 13.67 (mean), 13.69 (median), 7.76 (2.5%) and 19.49 (97.5%) 13 . CC-BY-NC-ND 4.0 International license It is made available under a author/funder, who has granted medRxiv a license to display the preprint in perpetuity. is the (which was not peer-reviewed) The copyright holder for this preprint . https://doi.org/10.1101/19008417 doi: medRxiv preprint CC-BY-NC-ND 4.0 International license It is made available under a author/funder, who has granted medRxiv a license to display the preprint in perpetuity. is the (which was not peer-reviewed) The copyright holder for this preprint . https://doi.org/10.1101/19008417 doi: medRxiv preprint Situation Report: Rohingya Refugee Crisis, Cox's Bazar | URL more-60-rohingya-babies-born-bangladesh-refugee-camps-every-day-unicef The Theory of Measles Elimination: Implications for the Design of Elimination Strategies. The 268 Different epidemic curves for severe acute respiratory syndrome reveal similar 270 impacts of control measures Serial Intervals of Respiratory Infectious Diseases: A 272 Systematic Review and Analysis Refugee Settlements Situation Report: Rohingya Refugee Crisis 000 Rohingya Babies Set to be Born in Overcrowded Refugee Camps in 281 Republic of Korea pledges US$ 500,000 to scale up UNFPA's response to the Rohingya refugee crisis in 284 99s-response-rohingya-refugee-crisis-bangladesh Prevention of Measles, Rubella Large Sample Properties of Simulations Using Latin Hypercube Sampling An age-structured model of pre-and post-vaccination measles transmission Mortality and Morbidity Weekly Bulletin (MMWB) Waning antibody levels and avidity: 298 implications for MMR vaccine-induced protection Impact of Public Health Responses During a Measles Outbreak in an Amish 300 Community in Ohio: Modeling the Dynamics of Transmission Health Sector Bulletin #9 Rohingya Crisis in Cox's Bazar The basic reproduction number (R0) of measles: a systematic review Generation interval contraction and epidemic data analysis Estimation of the reproductive number and the serial interval in early phase of the 2009 312 influenza A/H1n1 pandemic in the USA Weekly Situation Report #7, Rohingya Refugees Weekly Situation Report #57, Rohingya Refugees Situation Report #10, Rohingya Refugees Epidemiological Highlights