The role of time-varying viral shedding in modelling environmental surveillance for public health: revisiting the 2013 poliovirus outbreak in Israel.

Environmental pathogen surveillance is a sensitive tool that can detect early-stage outbreaks, and it is being used to track poliovirus and other pathogens. However, interpretation of longitudinal environmental surveillance signals is difficult because the relationship between infection incidence and viral load in wastewater depends on time-varying shedding intensity. We developed a mathematical model of time-varying poliovirus shedding intensity consistent with expert opinion across a range of immunization states. Incorporating this shedding model into an infectious disease transmission model, we analysed quantitative, polymerase chain reaction data from seven sites during the 2013 Israeli poliovirus outbreak. Compared to a constant shedding model, our time-varying shedding model estimated a slower peak (four weeks later), with more of the population reached by a vaccination campaign before infection and a lower cumulative incidence. We also estimated the population shed virus for an average of 29 days (95% CI 28-31), longer than expert opinion had suggested for a population that was purported to have received three or more inactivated polio vaccine (IPV) doses. One explanation is that IPV may not substantially affect shedding duration. Using realistic models of time-varying shedding coupled with longitudinal environmental surveillance may improve our understanding of outbreak dynamics of poliovirus, SARS-CoV-2, or other pathogens.

Assessing the sensitivity of the polio environmental surveillance system.

BACKGROUND: The polio environmental surveillance (ES) system has been an incredible tool for advancing polio eradication efforts because of its ability to highlight the spatial and temporal extent of poliovirus circulation. While ES often outperforms, or is more sensitive than AFP surveillance, the sensitivity of the ES system has not been well characterized. Fundamental uncertainty of ES site sensitivity makes it difficult to interpret results from ES, particularly negative results. METHODS AND FINDINGS: To study ES sensitivity, we used data from Afghanistan and Pakistan to examine the probability that each ES site detected the Sabin 1, 2, or 3 components of the oral polio vaccine (OPV) as a function of virus prevalence within the same district (estimated from AFP data). Accounting for virus prevalence is essential for estimating site sensitivity because Sabin detection rates should vary with prevalence-high immediately after supplemental immunization activities (SIAs), but low in subsequent months. We found that most ES sites in Pakistan and Afghanistan are highly sensitive for detecting poliovirus relative to AFP surveillance in the same districts. For example, even when Sabin poliovirus is at low prevalence of ~0.5-3% in AFP surveillance, most ES sites have ~34-50% probability of detecting Sabin. However, there was considerable variation in ES site sensitivity and we flagged several sites for re-evaluation based on low sensitivity rankings and low wild polio virus detection rates. In these areas, adding new sites or modifying collection methods in current sites could improve sensitivity of environmental surveillance. CONCLUSIONS: Relating ES detections to virus prevalence significantly improved our ability to evaluate site sensitivity compared to evaluations based solely on ES detection rates. To extend our approach to new sites and regions, we provide a preliminary framework for relating ES and AFP detection rates, and descriptions of how detection rates might relate to SIAs and natural seasonality.

Spatial model for risk prediction and sub-national prioritization to aid poliovirus eradication in Pakistan.

BACKGROUND: Pakistan is one of only three countries where poliovirus circulation remains endemic. For the Pakistan Polio Eradication Program, identifying high risk districts is essential to target interventions and allocate limited resources. METHODS: Using a hierarchical Bayesian framework we developed a spatial Poisson hurdle model to jointly model the probability of one or more paralytic polio cases, and the number of cases that would be detected in the event of an outbreak. Rates of underimmunization, routine immunization, and population immunity, as well as seasonality and a history of cases were used to project future risk of cases. RESULTS: The expected number of cases in each district in a 6-month period was predicted using indicators from the previous 6-months and the estimated coefficients from the model. The model achieves an average of 90% predictive accuracy as measured by area under the receiver operating characteristic (ROC) curve, for the past 3 years of cases. CONCLUSIONS: The risk of poliovirus has decreased dramatically in many of the key reservoir areas in Pakistan. The results of this model have been used to prioritize sub-national areas in Pakistan to receive additional immunization activities, additional monitoring, or other special interventions.

Evaluating cessation of the type 2 oral polio vaccine by modeling pre- and post-cessation detection rates.

The globally synchronized removal of the attenuated Sabin type 2 strain from the oral polio vaccine (OPV) in April 2016 marked a major change in polio vaccination policy. This change will provide a significant reduction in the burden of vaccine-associated paralytic polio (VAPP), but may increase the risk of circulating vaccine-derived poliovirus (cVDPV2) outbreaks during the transition period. This risk can be monitored by tracking the disappearance of Sabin-like type 2 (SL2) using data from the polio surveillance system. We studied SL2 prevalence in 17 countries in Africa and Asia, from 2010 to 2016 using acute flaccid paralysis surveillance data. We modeled the peak and decay of SL2 prevalence following mass vaccination events using a beta-binomial model for the detection rate, and a Ricker function for the temporal dependence. We found type 2 circulated the longest of all serotypes after a vaccination campaign, but that SL2 prevalence returned to baseline levels in approximately 50days. Post-cessation model predictions identified 19 anomalous SL2 detections outside of model predictions in Afghanistan, India, Nigeria, Pakistan, and western Africa. Our models established benchmarks for the duration of SL2 detection after OPV2 cessation. As predicted, SL2 detection rates have plummeted, except in Nigeria where OPV2 use continued for some time in response to recent cVDPV2 detections. However, the anomalous SL2 detections suggest specific areas that merit enhanced monitoring for signs of cVDPV2 outbreaks.

Recruitment limitation of long-lived conifers: implications for climate change responses.

Seed availability and suitable microsites for germination are likely to severely constrain the responses of plant species to climate change, especially at and beyond range edges. For example, range shifts may be slow if seed availability is low at range edges due to low parent-tree abundance or reduced fecundity. Even when seeds are available, climatic and biotic factors may further limit the availability of suitable microsites for recruitment. Unfortunately, the importance of seed and microsite limitation during range shifts remains unknown, since few studies have examined both factors simultaneously, particularly across species' ranges. To address this issue, we assessed seed availability and the factors influencing germination for six conifer species across a large environmental gradient encompassing their elevational ranges. Specifically, we assessed (1) how parent-tree abundance influences annual seed availability; (2) how seed limitation varies across species' ranges; (3) how climatic and biotic factors affect germination; and (4) how seed and suitable microsite availability covary annually within and among species. We found that seed availability declined toward species' upper range edges for most species, primarily due to low parent-tree abundance rather than declining fecundity. Range expansions are thus likely to be lagged with respect to climate change, as long generation times preclude rapid increases in tree density. Negative impacts of canopy cover on germination rates suggest range shifts will further be slowed by competition with existing vegetation. Moreover, years of high seed production were generally correlated among species, but not correlated with the availability of suitable microsites, implying that seedling competition and the interaction between seed and microsite limitation will further constrain recruitment. However, the nature of microsite limitation varied strongly between treeline and low-elevation species due to differing responses to snowpack duration and. competition, suggesting that treeline species may be quicker to shift their ranges in response to warming than low-elevation species. In all, our results demonstrate that seed and microsite limitation will likely result in lagged responses to climate change but with differences among species leading to complex range shift dynamics.