Evolutionary consequences of delaying intervention for monkeypox

Philip L F Johnson, Carl T Bergstrom, Roland R Regoes, Ira M Longini, M Elizabeth Halloran, Rustom Antia

The Lancet

September 21, 2022

ABSTRACT

Since May, 2022, clusters of monkeypox infections have caused global concern. At present, this concern has been tempered by the fact that, even when uncontrolled, the number of infections is growing slowly, indicating a reproductive number (R) not much larger than unity. However, the effect of R on the probability of evolution might not be obvious. We suggest that, compared with zoonotic pathogens with large R values, those pathogens with R values just above 1, such as monkeypox virus, have a higher probability of evolution during the timeframe in which the number of cases remains low. Waiting until the number of cases is high would give monkeypox virus—or any emerging pathogen—the opportunity to adapt substantially to humans.

A platform trial design for preventive vaccines against Marburg virus and other emerging infectious disease threats

Ira M Longini, Yang Yang, Thomas R Fleming, César Muñoz-Fontela, Rui Wang, Susan S Ellenberg, George Qian, M Elizabeth Halloran, Martha Nason, Victor De Gruttola, Sabue Mulangu, Yunda Huang, Christl A Donnelly, Ana-Maria Henao Restrepo

Clinical Trials

July 22, 2022

ABSTRACT

Background:

The threat of a possible Marburg virus disease outbreak in Central and Western Africa is growing. While no Marburg virus vaccines are currently available for use, several candidates are in the pipeline. Building on knowledge and experiences in the designs of vaccine efficacy trials against other pathogens, including SARS-CoV-2, we develop designs of randomized Phase 3 vaccine efficacy trials for Marburg virus vaccines.

Methods:

A core protocol approach will be used, allowing multiple vaccine candidates to be tested against controls. The primary objective of the trial will be to evaluate the effect of each vaccine on the rate of virologically confirmed Marburg virus disease, although Marburg infection assessed via seroconversion could be the primary objective in some cases. The overall trial design will be a mixture of individually and cluster-randomized designs, with individual randomization done whenever possible. Clusters will consist of either contacts and contacts of contacts of index cases, that is, ring vaccination, or other transmission units.

Results:

The primary efficacy endpoint will be analysed as a time-to-event outcome. A vaccine will be considered successful if its estimated efficacy is greater than 50% and has sufficient precision to rule out that true efficacy is less than 30%. This will require approximately 150 total endpoints, that is, cases of confirmed Marburg virus disease, per vaccine/comparator combination. Interim analyses will be conducted after 50 and after 100 events. Statistical analysis of the trial will be blended across the different types of designs. Under the assumption of a 6-month attack rate of 1% of the participants in the placebo arm for both the individually and cluster-randomized populations, the most likely sample size is about 20,000 participants per arm.

Conclusion:

This event-driven design takes into the account the potentially sporadic spread of Marburg virus. The proposed trial design may be applicable for other pathogens against which effective vaccines are not yet available.

Quantifying the importance and location of SARS-CoV-2 transmission events in large metropolitan areas

Alberto Aleta, David Mart´ın-Corral, Michiel A. Bakker, Ana Pastore y Piontti, Marco Ajelli, Maria Litvinova, Matteo Chinazzi,Natalie E. Dean, M. Elizabeth Halloran, Ira M. Longini, Jr., Alex Pentland, Alessandro Vespignani, Yamir Moreno, Esteban Moro

PNAS

June 28, 2022

ABSTRACT

Detailed characterization of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) transmission across different settings can help design less disruptive interventions. We used real-time, privacy-enhanced mobility data in the New York City, NY and Seattle, WA metropolitan areas to build a detailed agent-based model of SARS-CoV-2 infection to estimate the where, when, and magnitude of transmission events during the pandemic’s first wave. We estimate that only 18% of individuals produce most infections (80%), with about 10% of events that can be considered superspreading events (SSEs). Although mass gatherings present an important risk for SSEs, we estimate that the bulk of transmission occurred in smaller events in settings like workplaces, grocery stores, or food venues. The places most important for transmission change during the pandemic and are different across cities, signaling the large underlying behavioral component underneath them. Our modeling complements case studies and epidemiological data and indicates that real-time tracking of transmission events could help evaluate and define targeted mitigation policies.

Ecology and public health burden of Keystone virus in Florida

Christopher J.Henry, Alexander N.Pillai, John A. Lednicky, J. Glenn Morris Jr., Thomas J. Hladish

Epidemics

June 13, 2022

ABSTRACT

Keystone virus (KEYV) is an under-studied orthobunyavirus that is transmitted via both horizontal and vertical cycles involving various mosquito species and vertebrate hosts. Historical evidence indicates that KEYV causes sub-clinical infections in humans, and some case studies draw links between this virus and encephalitis. Given KEYV’s potential to cause human infections, it is plausible that it causes an under-appreciated proportion of both generic viral infections and unidentified viral encephalitis cases. This review summarizes current knowledge of KEYV and its disease dynamics in order to better understand the virus’ medical and economic burden on human populations.

Household secondary attack rates of SARS-CoV-2 by variant and vaccination status: an updated systematic review and meta-analysis

Zachary J. Madewell, Yang Yang, Ira M. Longini Jr, M. Elizabeth Halloran, Natalie E. Dean

medrXiv

January 11, 2022

ABSTRACT

We previously reported a household secondary attack rate (SAR) for SARS-CoV-2 of 18.9% through June 17, 2021. To examine how emerging variants and increased vaccination have affected transmission rates, we searched PubMed from June 18, 2021, through January 7, 2022. Meta-analyses used generalized linear mixed models to obtain SAR estimates and 95%CI, disaggregated by several covariates. SARs were used to estimate vaccine effectiveness based on the transmission probability for susceptibility (VES,p), infectiousness (VEI,p), and total vaccine effectiveness (VET,p). Household SAR for 27 studies with midpoints in 2021 was 35.8% (95%CI, 30.6%-41.3%), compared to 15.7% (95%CI, 13.3%-18.4%) for 62 studies with midpoints through April 2020. Household SARs were 38.0% (95%CI, 36.0%-40.0%), 30.8% (95%CI, 23.5%-39.3%), and 22.5% (95%CI, 18.6%-26.8%) for Alpha, Delta, and Beta, respectively. VEI,p, VES,p, and VET,p were 56.6% (95%CI, 28.7%-73.6%), 70.3% (95%CI, 59.3%-78.4%), and 86.8% (95%CI, 76.7%-92.5%) for full vaccination, and 27.5% (95%CI, -6.4%-50.7%), 43.9% (95%CI, 21.8%-59.7%), and 59.9% (95%CI, 34.4%-75.5%) for partial vaccination, respectively. Household contacts exposed to Alpha or Delta are at increased risk of infection compared to the original wild-type strain. Vaccination reduced susceptibility to infection and transmission to others.

Using social contact data to improve the overall effect estimate of a cluster-randomized influenza vaccination program in Senegal

Gail E Potter, Nicole Bohme Carnegie, Jonathan D Sugimoto, Aldiouma Diallo, John C Victor, Kathleen M Neuzil , M Elizabeth Halloran

Royal Statistical Society C

January 10, 2022

ABSTRACT

This study estimates the overall effect of two influenza vaccination programs consecutively administered in a cluster-randomized trial in western Senegal over the course of two influenza seasons from 2009-2011. We apply cutting-edge methodology combining social contact data with infection data to reduce bias in estimation arising from contamination between clusters. Our time-varying estimates reveal a reduction in seasonal influenza from the intervention and a nonsignificant increase in H1N1 pandemic influenza. We estimate an additive change in overall cumulative incidence (which was 6.13% in the control arm) of -0.68 percentage points during Year 1 of the study (95% CI: -2.53, 1.18). When H1N1 pandemic infections were excluded from analysis, the estimated change was -1.45 percentage points and was significant (95% CI, -2.81, -0.08). Because cross-cluster contamination was low (0-3% of contacts for most villages), an estimator assuming no contamination was only slightly attenuated (-0.65 percentage points). These findings are encouraging for studies carefully designed to minimize spillover. Further work is needed to estimate contamination - and its effect on estimation - in a variety of settings.

Quantifying the importance and location of SARS-CoV-2 transmission events in large metropolitan areas

Alberto Aleta, David Martín-Corral, Michiel A. Bakker, Ana Pastore y Piontti, Marco Ajelli, Maria Litvinova, Matteo Chinazzi, Natalie E. Dean, M. Elizabeth Halloran, Ira M. Longini Jr., Alex Pentland, Alessandro Vespignani, Yamir Moreno, Esteban Moro

medRxiv

December 26, 2021

ABSTRACT

Detailed characterization of SARS-CoV-2 transmission across different settings can help design less disruptive interventions. We used real-time, privacy-enhanced mobility data in the New York City and Seattle metropolitan areas to build a detailed agent-based model of SARS-CoV-2 infection to estimate the where, when, and magnitude of transmission events during the pandemic’s first wave. We estimate that only 18% of individuals produce most infections (80%), with about 10% of events that can be considered super-spreading events (SSEs). Although mass-gatherings present an important risk for SSEs, we estimate that the bulk of transmission occurred in smaller events in settings like workplaces, grocery stores, or food venues. The places most important for transmission change during the pandemic and are different across cities, signaling the large underlying behavioral component underneath them. Our modeling complements case studies and epidemiological data and indicates that real-time tracking of transmission events could help evaluate and define targeted mitigation policies.

Markov Genealogy Processes

AARON A. KING, QIANYING LIN, EDWARD L. IONIDES

Theoretical Population Biology

December 9, 2021

Abstract

We construct a family of genealogy-valued Markov processes that are induced by a continuous-time Markov population process. We derive exact expressions for the likelihood of a given genealogy conditional on the history of the underlying population process. These lead to a nonlinear filtering equation which can be used to design efficient Monte Carlo inference algorithms. We demonstrate these calculations with several examples. Existing full-information approaches for phylodynamic inference are special cases of the theory.

Cryptic transmission of SARS-CoV-2 and the first COVID-19 wave

Jessica T. Davis, Matteo Chinazzi, Nicola Perra, Kunpeng Mu, Ana Pastore y Piontti, Marco Ajelli, Natalie E. Dean, Corrado Gioannini, Maria Litvinova, Stefano Merler, Luca Rossi, Kaiyuan Sun, Xinyue Xiong, Ira M. Longini Jr, M. Elizabeth Halloran, Cécile Viboud & Alessandro Vespignani

Nature

October 25, 2021

ABSTRACT

Considerable uncertainty surrounds the timeline of introductions and onsets of local transmission of SARS-CoV-2 globally1–7. Although a limited number of SARS-CoV-2 introductions were reported in January and February 20208,9, the narrowness of the initial testing criteria, combined with a slow growth in testing capacity and porous travel screening10, left many countries vulnerable to unmitigated, cryptic transmission. Here we use a global metapopulation epidemic model to provide a mechanistic understanding of the early dispersal of infections, and the temporal windows of the introduction and onset of SARS-CoV-2 local transmission in Europe and the United States. We find that community transmission of SARS-CoV-2 was likely in several areas of Europe and the United States by January 2020, and estimate that by early March, only 1 to 3 in 100 SARS-CoV-2 infections were detected by surveillance systems. The modelling results highlight international travel as the key driver of the introduction of SARS-CoV-2 with possible introductions and transmission events as early as December 2019–January 2020. We find a heterogeneous, geographic distribution of cumulative infection attack rates by 4 July 2020, ranging from 0.78%–15.2% across US states and 0.19%–13.2% in European countries. Our approach complements phylogenetic analyses and other surveillance approaches and provides insights that can be used to design innovative, model-driven surveillance systems that guide enhanced testing and response strategies.

Challenges of evaluating and modelling vaccination in emerging infectious diseases

Zachary J Madewell, Natalie E Dean, Jesse A Berlin, Paul M Coplan, Kourtney J Davis, Claudio J Struchiner, M Elizabeth Halloran

Epidemics

October 5, 2021

ABSTRACT

Outbreaks of emerging pathogens pose unique methodological and practical challenges for the design, implementation, and evaluation of vaccine efficacy trials. Lessons learned from COVID-19 highlight the need for innovative and flexible study design and application to quickly identify promising candidate vaccines. Trial design strategies should be tailored to the dynamics of the specific pathogen, location of the outbreak, and vaccine prototypes, within the regional socioeconomic constraints. Mathematical and statistical models can assist investigators in designing infectious disease clinical trials. We introduce key challenges for planning, evaluating, and modelling vaccine efficacy trials for emerging pathogens.

Transition to endemicity: Understanding COVID-19

Rustom Antia, M. Elizabeth Halloran

Immunity

September 23, 2021

ABSTRACT

The emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and its associated disease, coronavirus disease 2019 (COVID-19), has caused a devastating pandemic worldwide. Here, we explain basic concepts underlying the transition from an epidemic to an endemic state, where a pathogen is stably maintained in a population. We discuss how the number of infections and the severity of disease change in the transition from the epidemic to the endemic phase and consider the implications of this transition in the context of COVID-19.

Durability of protection after 5 doses of acellular pertussis vaccine among 5–9 year old children in King County, Washington

Madhura S. Rane, Pejman Rohani, M. Elizabeth Halloran

Vaccine

September 4, 2021

ABSTRACT

Purpose
Waning of immunity after vaccination with the acellular Pertussis (aP) vaccine has been proposed as one of the main reasons for pertussis resurgence in the US. In this study, we estimated time-varying vaccine effectiveness after 5 doses of aP vaccine.

Methods
We conducted a retrospective cohort study among children 5–9 years old (born between 2008 and 2012) living in King County, Washington, USA, who participated in the Washington State Immunization Information System. We estimated time-varying vaccine effectiveness after 5 doses of aP using smoothed scaled Schoenfeld residuals obtained from fitting Cox proportional hazards models to the data as well as piecewise constant Poisson regression.

Results
There were 55 pertussis cases in this cohort, of whom 22 (40%) were fully-vaccinated and 33 (60%) were under-vaccinated. Vaccine effectiveness (VE) remained high for up to 42 months after the fifth dose (VE(t) = 89%; 95% CI: 64%, 97%) as estimated using survival analysis methods and up to 4 years (VE(t) = 93%; 95% CI: 67%, 98%) as estimated using Poisson regression.

Conclusion
We did not find evidence for waning of vaccine effectiveness for up to four years after 5 doses of aP among 5 –9 years old children in King County, WA.

Association of Diphtheria-Tetanus–Acellular Pertussis Vaccine Timeliness and Number of Doses With Age-Specific Pertussis Risk in Infants and Young Children

Madhura S. Rane, Pejman Rohani, M. Elizabeth Halloran

JAMA
August 10, 2021

ABSTRACT

Importance In most countries, the diphtheria-tetanus–acellular pertussis (DTaP) vaccine is administered as a 3-dose infant series followed by additional booster doses in the first 5 years of life. Short-term immunity from the DTaP vaccine can depend on the number, timing, and interval between doses. Not receiving doses in a timely manner might be associated with a higher pertussis risk.

Objective To examine the association between number and timeliness of vaccine doses and age-specific pertussis risk.

Design, Setting, and Participants This population-based, retrospective cohort study used Washington State Immunization Information System data and pertussis surveillance data from Public Health Seattle and King County, Washington. Included participants were children aged 3 months to 9 years born or living in King County, Washington, between January 1, 2008, and December 31, 2017. Data were analyzed from June 30 to December 1, 2019.

Exposures Being undervaccinated (receiving fewer than recommended doses at a given age) or delayed vaccination (not receiving doses within time frames recommended by Centers for Disease Control and Prevention).

Main Outcomes and Measures Suspected, probable, and confirmed pertussis diagnosis.

Results A total of 316 404 children (median age, 65.2 months [interquartile range, 35.3-94.1 months]; 162 025 boys [51.2%]) as of December 31, 2017, with 17.4 million person-months of follow-up were included in the analysis. A total of 19 943 children (6.3%) had no vaccines recorded in the Immunization Information System, 116 193 (36.7%) received a vaccine with a delay, and 180 268 (56.9%) were fully vaccinated with no delay. Delayed vaccination and undervaccination rates were higher for older children (17.6% delayed or undervaccinated at age 2 months for dose 1 at 3 months vs 41.6% at age 5 years for dose 5) but improved for successive birth cohorts (52.2% for 2008 birth cohort vs 32.3% for 2017 birth cohort). Undervaccination was significantly associated with higher risk of pertussis for the 3-dose primary series (adjusted relative risk [aRR], 4.8; 95% CI, 3.1-7.6), the first booster (aRR, 3.2; 95% CI, 2.3-4.5), and the second booster (aRR, 4.6; 95% CI, 2.6-8.2). However, delay in vaccination among children who received the recommended number of vaccine doses was not associated with pertussis risk.

Conclusions and Relevance The results of this cohort study suggest that undervaccination is associated with higher pertussis risk. Short delays in vaccine receipt may be less important if the age-appropriate number of doses is administered, but delaying doses is not recommended. Ensuring that children receive all doses of pertussis vaccine, even if there is some delay, is important.

Using simulated infectious disease outbreaks to inform site selection and sample size for individually randomized vaccine trials during an ongoing epidemic

Zachary J Madewell, Ana Pastore Y Piontti, Qian Zhang, Nathan Burton, Yang Yang, Ira M Longini, M Elizabeth Halloran, Alessandro Vespignani, Natalie E Dean

Clinical Trials

July 3, 2021

ABSTRACT

Background:
Novel strategies are needed to make vaccine efficacy trials more robust given uncertain epidemiology of infectious disease outbreaks, such as arboviruses like Zika. Spatially resolved mathematical and statistical models can help investigators identify sites at highest risk of future transmission and prioritize these for inclusion in trials. Models can also characterize uncertainty in whether transmission will occur at a site, and how nearby or connected sites may have correlated outcomes. A structure is needed for how trials can use models to address key design questions, including how to prioritize sites, the optimal number of sites, and how to allocate participants across sites.

Methods:
We illustrate the added value of models using the motivating example of Zika vaccine trial planning during the 2015–2017 Zika epidemic. We used a stochastic, spatially resolved, transmission model (the Global Epidemic and Mobility model) to simulate epidemics and site-level incidence at 100 high-risk sites in the Americas. We considered several strategies for prioritizing sites (average site-level incidence of infection across epidemics, median incidence, probability of exceeding 1% incidence), selecting the number of sites, and allocating sample size across sites (equal enrollment, proportional to average incidence, proportional to rank). To evaluate each design, we stochastically simulated trials in each hypothetical epidemic by drawing observed cases from site-level incidence data.

Results:
When constraining overall trial size, the optimal number of sites represents a balance between prioritizing highest-risk sites and having enough sites to reduce the chance of observing too few endpoints. The optimal number of sites remained roughly constant regardless of the targeted number of events, although it is necessary to increase the sample size to achieve the desired power. Though different ranking strategies returned different site orders, they performed similarly with respect to trial power. Instead of enrolling participants equally from each site, investigators can allocate participants proportional to projected incidence, though this did not provide an advantage in our example because the top sites had similar risk profiles. Sites from the same geographic region may have similar outcomes, so optimal combinations of sites may be geographically dispersed, even when these are not the highest ranked sites.


Conclusion:
Mathematical and statistical models may assist in designing successful vaccination trials by capturing uncertainty and correlation in future transmission. Although many factors affect site selection, such as logistical feasibility, models can help investigators optimize site selection and the number and size of participating sites. Although our study focused on trial design for an emerging arbovirus, a similar approach can be made for any infectious disease with the appropriate model for the particular disease.

Modeling the impacts of clinical influenza testing on influenza vaccine effectiveness estimates

Leora R Feldstein, Jill M Ferdinands, Wesley H Self, Adrienne G Randolph, Michael Aboodi, Adrienne H Baughman, Samuel M Brown, Matthew C Exline, D Clark Files, Kevin Gibbs, Adit A Ginde, Michelle N Gong, Carlos G Grijalva, Natasha Halasa, Akram Khan, Christopher J Lindsell, Margaret Newhams, Ithan D Peltan, Matthew E Prekker, Todd W Rice, Nathan I Shapiro, Jay Steingrub, H Keipp Talbot, M Elizabeth Halloran, Manish Patel

Journal of Infectious Diseases

ABSTRACT

Background

Test-negative design studies for evaluating influenza vaccine effectiveness (VE) enroll patients with acute respiratory infection. Enrollment typically occurs before influenza status is determined, resulting in over-enrollment of influenza-negative patients. With availability of rapid and accurate molecular clinical testing, influenza status could be ascertained prior to enrollment, thus improving study efficiency. We estimate potential biases in VE when using clinical testing.

Methods

We simulate data assuming 60% vaccinated, 25% of those vaccinated are influenza positive, and VE of 50%. We show the effect on VE in five scenarios.

Results

VE is affected only when clinical testing preferentially targets patients based on both vaccination and influenza status. VE is overestimated by 10% if non-testing occurs in 39% of vaccinated influenza-positive patients and 24% of others; and if non-testing occurs in 8% of unvaccinated influenza-positive patients and 27% of others. VE is underestimated by 10% if non-testing occurs in 32% of unvaccinated influenza-negative patients and 18% of others.

Conclusions

Although differential clinical testing by vaccine receipt and influenza positivity may produce errors in estimated VE, bias in testing would have to be substantial and overall proportion of patients tested would have to be small to result in a meaningful difference in VE.

Estimating population-level effects of the acellular pertussis vaccine using routinely collected immunization data

Madhura S Rane, M Elizabeth Halloran

Clinical Infectious Diseases

April 21, 2021

ABSTRACT

Background
Measuring and reporting the different population-level effects of the acellular pertussis vaccine on pertussis disease in addition to direct effects can increase the cost-effectiveness of a vaccine.

Methods
We conducted a retrospective cohort study of children born between January 1, 2008, and December 31, 2017, in King County, Washington, who were enrolled in the Washington State Immunization Information System. Diphtheria-Tetanus-acellular-Pertussis (DTaP) vaccination data from WA-IIS was linked with pertussis case data from Public Health Seattle and King County. Census-level vaccination coverage was estimated as proportion of age-appropriately vaccinated children residing in it. Direct vaccine effectiveness was estimated by comparing pertussis risk in fully-vaccinated and under-vaccinated children. Population-level vaccine effects were estimated by comparing pertussis risk in census tracts in the highest vaccination coverage quartile to that in the lowest vaccination coverage quartile.

Results
For direct protection, estimated vaccine effectiveness was 76% (95% CI: 63% - 84%) in low vaccination coverage clusters and it decreased to 47% (95% CI: 13% - 68%) in high vaccination coverage clusters, after adjusting for potential confounders. The estimated indirect effect was 45.0% (95% CI: 1%, 70%), total effect was 93.9% (95% CI: 91%, 96%), and overall effect was 42.2% (95% CI: 19%, 60%).

Conclusion
Our findings suggest that DTaP vaccination provided direct as well as indirect protection in the highly immunized King County, WA. Routine DTaP vaccination programs may have the potential to provide not only protection for vaccinated individuals but also for the under-vaccinated individuals living in the same area.

Inverse probability weighted estimators of vaccine effects accommodating partial interference and censoring

Sujatro Chakladar, Samuel Rosin, Michael G. Hudgens, M. Elizabeth Halloran, John D. Clemens, Mohammad Ali, Michael E. Emch

Biometrics

March 25, 2021

ABSTRACT

Estimating population-level effects of a vaccine is challenging because there may be interference, that is, the outcome of one individual may depend on the vaccination status of another individual. Partial interference occurs when individuals can be partitioned into groups such that interference occurs only within groups. In the absence of interference, inverse probability weighted (IPW) estimators are commonly used to draw inference about causal effects of an exposure or treatment. Tchetgen Tchetgen and VanderWeele proposed a modified IPW estimator for causal effects in the presence of partial interference. Motivated by a cholera vaccine study in Bangladesh, this paper considers an extension of the Tchetgen Tchetgen and VanderWeele IPW estimator to the setting where the outcome is subject to right censoring using inverse probability of censoring weights (IPCW). Censoring weights are estimated using proportional hazards frailty models. The large sample properties of the IPCW estimators are derived, and simulation studies are presented demonstrating the estimators' performance in finite samples. The methods are then used to analyze data from the cholera vaccine study.

Cryptic transmission of SARS-CoV-2 and the first COVID-19 wave in 2 Europe and the United States

Jessica T Davis, Matteo Chinazzi, Nicola Perra, Kunpeng Mu, Ana Pastore Y Piontti, Marco Ajelli, Natalie E Dean, Corrado Gioannini, Maria Litvinova, Stefano Merler, Luca Rossi, Kaiyuan Sun, Xinyue Xiong, M Elizabeth Halloran, Ira M Longini, Cécile Viboud, Alessandro Vespignani

medRxiv

March 25, 2021

ABSTRACT

Given the narrowness of the initial testing criteria, the SARS-CoV-2 virus spread through cryptic transmission in January and February, setting the stage for the epidemic wave experienced in March and April, 2020. We use a global metapopulation epidemic model to provide a mechanistic understanding of the global dynamic underlying the establishment of the COVID-19 pandemic in Europe and the United States (US). The model is calibrated on international case introductions at the early stage of the pandemic. We find that widespread community transmission of SARS-CoV-2 was likely in several areas of Europe and the US by January 2020, and estimate that by early March, only 1 - 3 in 100 SARS-CoV-2 infections were detected by surveillance systems. Modeling results indicate international travel as the key driver of the introduction of SARS-CoV-2 with possible importation and transmission events as early as December, 2019. We characterize the resulting heterogeneous spatio-temporal spread of SARS-CoV-2 and the burden of the first COVID-19 wave (February-July 2020). We estimate infection attack rates ranging from 0.78%-15.2% in the US and 0.19%-13.2% in Europe. The spatial modeling of SARS-CoV-2 introductions and spreading provides insights into the design of innovative, model-driven surveillance systems and preparedness plans that have a broader initial capacity and indication for testing.

Estimating and interpreting secondary attack risk: Binomial considered biased

Yushuf Sharker, Eben Kenah

PLOS Computational Biology

January 20, 2021

ABSTRACT

The household secondary attack risk (SAR), often called the secondary attack rate or secondary infection risk, is the probability of infectious contact from an infectious household member A to a given household member B, where we define infectious contact to be a contact sufficient to infect B if he or she is susceptible. Estimation of the SAR is an important part of understanding and controlling the transmission of infectious diseases. In practice, it is most often estimated using binomial models such as logistic regression, which implicitly attribute all secondary infections in a household to the primary case. In the simplest case, the number of secondary infections in a household with m susceptibles and a single primary case is modeled as a binomial(m, p) random variable where p is the SAR. Although it has long been understood that transmission within households is not binomial, it is thought that multiple generations of transmission can be neglected safely when p is small. We use probability generating functions and simulations to show that this is a mistake. The proportion of susceptible household members infected can be substantially larger than the SAR even when p is small. As a result, binomial estimates of the SAR are biased upward and their confidence intervals have poor coverage probabilities even if adjusted for clustering. Accurate point and interval estimates of the SAR can be obtained using longitudinal chain binomial models or pairwise survival analysis, which account for multiple generations of transmission within households, the ongoing risk of infection from outside the household, and incomplete follow-up. We illustrate the practical implications of these results in an analysis of household surveillance data collected by the Los Angeles County Department of Public Health during the 2009 influenza A (H1N1) pandemic.