Search

Studying the early events that occur after viral infection in humans is difficult unless one intentionally infects volunteers in a human challenge study. Here, we use data about severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in such a study in combination with mathematical modeling to gain insights into the relationship between the amount of virus in the upper respiratory tract and the immune response it generates. We propose a set of dynamic models of increasing complexity to dissect the roles of target cell limitation, innate immunity, and adaptive immunity in determining the observed viral kinetics. We introduce an approach for modeling the effect of humoral immunity that describes a decline in infectious virus after immune activation. We fit our models to viral load and infectious titer data from all the untreated infected participants in the study simultaneously. We found that a power-law with a power h < 1 describes the relationship between infectious virus and viral load. Viral replication at the early stage of infection is rapid, with a doubling time of ~2 h for viral RNA and ~3 h for infectious virus. We estimate that adaptive immunity is initiated ~7 to 10 d postinfection and appears to contribute to a multiphasic viral decline experienced by some participants; the viral rebound experienced by other participants is consistent with a decline in the interferon response. Altogether, we quantified the kinetics of SARS-CoV-2 infection, shedding light on the early dynamics of the virus and the potential role of innate and adaptive immunity in promoting viral decline during infection., Competing Interests: Competing interests statement:The authors declare no competing interest.

Chronic hepatitis B virus (HBV) infection is strongly associated with increased risk of liver cancer and cirrhosis. While existing treatments effectively inhibit the HBV life cycle, viral rebound occurs rapidly following treatment interruption. Consequently, functional cure rates of chronic HBV infection remain low and there is increased interest in a novel treatment modality, capsid assembly modulators (CAMs). Here, we develop a multiscale mathematical model of CAM treatment in chronic HBV infection. By fitting the model to participant data from a phase I trial of the first-generation CAM vebicorvir, we estimate the drug's dose-dependent effectiveness and identify the physiological mechanisms that drive the observed biphasic decline in HBV DNA and RNA, and mechanistic differences between HBeAg-positive and negative infection. Finally, we demonstrate analytically and numerically that HBV RNA is more sensitive than HBV DNA to increases in CAM effectiveness.

In 2022, there was a global resurgence of mpox, with different clinical-epidemiological features compared with previous outbreaks. Sexual contact was hypothesized as the primary transmission route, and the community of men having sex with men (MSM) was disproportionately affected. Because of the stigma associated with sexually transmitted infections, the real burden of mpox could be masked. We quantified the basic reproduction number ( R 0 ) and the underestimated fraction of mpox cases in 16 countries, from the onset of the outbreak until early September 2022, using Bayesian inference and a compartmentalized, risk-structured (high-/low-risk populations) and two-route (sexual/non-sexual transmission) mathematical model. Machine learning (ML) was harnessed to identify underestimation determinants. Estimated R 0 ranged between 1.37 (Canada) and 3.68 (Germany). The underestimation rates for the high- and low-risk populations varied between 25-93% and 65-85%, respectively. The estimated total number of mpox cases, relative to the reported cases, is highest in Colombia (3.60) and lowest in Canada (1.08). In the ML analysis, two clusters of countries could be identified, differing in terms of attitudes towards the 2SLGBTQIAP+ community and the importance of religion. Given the substantial mpox underestimation, surveillance should be enhanced, and country-specific campaigns against the stigmatization of MSM should be organized, leveraging community-based interventions.

The recent mpox outbreak (in 2022-2023) has different clinical and epidemiological features compared with previous outbreaks of the disease. During this outbreak, sexual contact was believed to be the primary transmission route of the disease. In addition, the community of men having sex with men (MSM) was disproportionately affected by the outbreak. This population is also disproportionately affected by HIV infection. Given that both diseases can be transmitted sexually, the endemicity of HIV, and the high sexual behavior associated with the MSM community, it is essential to understand the effect of the two diseases spreading simultaneously in an MSM population. Particularly, we aim to understand the potential effects of HIV on an mpox outbreak in the MSM population. We develop a mechanistic mathematical model of HIV and mpox co-infection. Our model incorporates the dynamics of both diseases and considers HIV treatment with anti-retroviral therapy (ART). In addition, we consider a potential scenario where HIV infection increases susceptibility to mpox, and investigate the potential impact of this mechanism on mpox dynamics. Our analysis shows that HIV can facilitate the spread of mpox in an MSM population, and that HIV treatment with ART may not be sufficient to control the spread of mpox in the population. However, we showed that a moderate use of condoms or reduction in sexual contact in the population combined with ART is beneficial in controlling mpox transmission. Based on our analysis, it is evident that effective control of HIV, specifically through substantial ART use, moderate condom compliance, and reduction in sexual contact, is imperative for curtailing the transmission of mpox in an MSM population and mitigating the compounding impact of these intertwined epidemics., Competing Interests: The authors declare that there is no competing interest with this publication titled “Understanding the impact of HIV on mpox transmission in an MSM population: a mathematical modeling study’’., (© 2024 The Authors.)

The outbreak of the severe acute respiratory syndrome coronavirus 2 started in Wuhan, China, towards the end of 2019 and spread worldwide. The rapid spread of the disease can be attributed to many factors including its high infectiousness and the high rate of human mobility around the world. Although travel/movement restrictions and other non-pharmaceutical interventions aimed at controlling the disease spread were put in place during the early stages of the pandemic, these interventions did not stop COVID-19 spread. To better understand the impact of human mobility on the spread of COVID-19 between regions, we propose a hybrid gravity-metapopulation model of COVID-19. Our modeling framework has the flexibility of determining mobility between regions based on the distances between the regions or using data from mobile devices. In addition, our model explicitly incorporates time-dependent human mobility into the disease transmission rate, and has the potential to incorporate other factors that affect disease transmission such as facemasks, physical distancing, contact rates, etc. An important feature of this modeling framework is its ability to independently assess the contribution of each factor to disease transmission. Using a Bayesian hierarchical modeling framework, we calibrate our model to the weekly reported cases of COVID-19 in thirteen local health areas in Metro Vancouver, British Columbia (BC), Canada, from July 2020 to January 2021. We consider two main scenarios in our model calibration: using a fixed distance matrix and time-dependent weekly mobility matrices. We found that the distance matrix provides a better fit to the data, whilst the mobility matrices have the ability to explain the variance in transmission between regions. This result shows that the mobility data provides more information in terms of disease transmission than the distances between the regions., Competing Interests: The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright: © 2023 Iyaniwura et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

Monkeypox, a zoonotic disease, is emerging as a potential sexually transmitted infection/disease, with underlying transmission mechanisms still unclear. We devised a risk-structured, compartmental model, incorporating sexual behavior dynamics. We compared different strategies targeting the high-risk population: a scenario of control policies geared toward the use of condoms and/or sexual abstinence (robust control strategy) with risk compensation behavior change, and a scenario of control strategies with behavior change in response to the doubling rate (adaptive control strategy). Monkeypox's basic reproduction number is 1.464, 0.0066, and 1.461 in the high-risk, low-risk, and total populations, respectively, with the high-risk group being the major driver of monkeypox spread. Policies imposing condom use or sexual abstinence need to achieve a 35% minimum compliance rate to stop further transmission, while a combination of both can curb the spread with 10% compliance to abstinence and 25% to condom use. With risk compensation, the only option is to impose sexual abstinence by at least 35%. Adaptive control is more effective than robust control where the daily sexual contact number is reduced proportionally and remains constant thereafter, shortening the time to epidemic peak, lowering its size, facilitating disease attenuation, and playing a key role in controlling the current outbreak., (© 2023 The Authors. Journal of Medical Virology published by Wiley Periodicals LLC.)

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has been spreading worldwide for over two years, with millions of reported cases and deaths. The deployment of mathematical modeling in the fight against COVID-19 has recorded tremendous success. However, most of these models target the epidemic phase of the disease. The development of safe and effective vaccines against SARS-CoV-2 brought hope of safe reopening of schools and businesses and return to pre-COVID normalcy, until mutant strains like the Delta and Omicron variants, which are more infectious, emerged. A few months into the pandemic, reports of the possibility of both vaccine- and infection-induced immunity waning emerged, thereby indicating that COVID-19 may be with us for longer than earlier thought. As a result, to better understand the dynamics of COVID-19, it is essential to study the disease with an endemic model. In this regard, we developed and analyzed an endemic model of COVID-19 that incorporates the waning of both vaccine- and infection-induced immunities using distributed delay equations. Our modeling framework assumes that the waning of both immunities occurs gradually over time at the population level. We derived a nonlinear ODE system from the distributed delay model and showed that the model could exhibit either a forward or backward bifurcation depending on the immunity waning rates. Having a backward bifurcation implies that $ R&#95;c < 1 $ is not sufficient to guarantee disease eradication, and that the immunity waning rates are critical factors in eradicating COVID-19. Our numerical simulations show that vaccinating a high percentage of the population with a safe and moderately effective vaccine could help in eradicating COVID-19.

Background: Vaccine hesitancy threatens efforts to bring the coronavirus disease 2019 (COVID-19) pandemic to an end. Given that social or interpersonal contact is an important driver for COVID-19 transmission, understanding the relationship between contact rates and vaccine hesitancy may help identify appropriate targets for strategic intervention. The purpose of this study was to assess the association between interpersonal contact and COVID-19 vaccine hesitancy among a sample of unvaccinated adults in the Canadian province of British Columbia (BC)., Methods: Unvaccinated individuals participating in the BC COVID-19 Population Mixing Patterns Survey (BC-Mix) were asked to indicate their level of agreement to the statement, " I plan to get the COVID-19 vaccine ." Multivariable multinomial logistic regression was used to assess the association between self-reported interpersonal contact and vaccine hesitancy, adjusting for age, sex, ethnicity, educational attainment, occupation, household size and region of residence. All analyses incorporated survey sampling weights based on age, sex, geography, and ethnicity., Results: Results were based on survey responses collected between March 8, 2021 and December 6, 2021, by a total of 4,515 adults aged 18 years and older. Overall, 56.7% of respondents reported that they were willing to get the COVID-19 vaccine, 27.0% were unwilling and 16.3% were undecided. We found a dose-response association between interpersonal contact and vaccine hesitancy. Compared to individuals in the lowest quartile (least contact), those in the fourth quartile (highest contact), third quartile and second quartile groups were more likely to be vaccine hesitant, with adjusted odd ratios (aORs) of 2.85 (95% CI: 2.02, 4.00), 1.91(95% CI: 1.38, 2.64), 1.78 (95% CI: 1.13, 2.82), respectively., Conclusion: Study findings show that among unvaccinated people in BC, vaccine hesitancy is greater among those who have high contact rates, and hence potentially at higher risk of acquiring and transmitting infection. This may also impact future uptake of booster doses., Competing Interests: Author NJ has advised and spoken for AbbVie, not related to this project or COVID-19. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Adu, Iyaniwura, Mahmood, Jeong, Makuza, Cua, Binka, García, Ringa, Wong, Yu, Irvine, Otterstatter and Janjua.)

Monkeypox is an emerging zoonotic disease caused by the monkeypox virus, which is an infectious agent belonging to the genus Orthopoxvirus . Currently, commencing from the end of April 2022, an outbreak of monkeypox is ongoing, with more than 43,000 cases reported as of 23 August 2022, involving 99 countries and territories across all the six World Health Organization (WHO) regions. On 23 July 2022, the Director-General of the WHO declared monkeypox a global public health emergency of international concern (PHEIC), since the outbreak represents an extraordinary, unusual, and unexpected event that poses a significant risk for international spread, requiring an immediate, coordinated international response. However, the real magnitude of the burden of disease could be masked by failures in ascertainment and under-detection. As such, underestimation affects the efficiency and reliability of surveillance and notification systems and compromises the possibility of making informed and evidence-based policy decisions in terms of the adoption and implementation of ad hoc adequate preventive measures. In this review, synthesizing 53 papers, we summarize the determinants of the underestimation of sexually transmitted diseases, in general, and, in particular, monkeypox, in terms of all their various components and dimensions (under-ascertainment, underreporting, under-detection, under-diagnosis, misdiagnosis/misclassification, and under-notification)., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Bragazzi, Woldegerima, Iyaniwura, Han, Wang, Shausan, Badu, Okwen, Prescod, Westin, Omame, Converti, Mellado, Wu and Kong.)

Following the emergence of COVID-19 at the end of 2019, several mathematical models have been developed to study the transmission dynamics of this disease. Many of these models assume homogeneous mixing in the underlying population. However, contact rates and mixing patterns can vary dramatically among individuals depending on their age and activity level. Variation in contact rates among age groups and over time can significantly impact how well a model captures observed trends. To properly model the age-dependent dynamics of COVID-19 and understand the impacts of interventions, it is essential to consider heterogeneity arising from contact rates and mixing patterns. We developed an age-structured model that incorporates time-varying contact rates and population mixing computed from the ongoing BC Mix COVID-19 survey to study transmission dynamics of COVID-19 in British Columbia (BC), Canada. Using a Bayesian inference framework, we fit four versions of our model to weekly reported cases of COVID-19 in BC, with each version allowing different assumptions of contact rates. We show that in addition to incorporating age-specific contact rates and mixing patterns, time-dependent (weekly) contact rates are needed to adequately capture the observed transmission dynamics of COVID-19. Our approach provides a framework for explicitly including empirical contact rates in a transmission model, which removes the need to otherwise model the impact of many non-pharmaceutical interventions. Further, this approach allows projection of future cases based on clear assumptions of age-specific contact rates, as opposed to less tractable assumptions regarding transmission rates., (Copyright © 2022 The Author(s). Published by Elsevier B.V. All rights reserved.)

We extended a class of coupled PDE-ODE models for studying the spatial spread of airborne diseases by incorporating human mobility. Human populations are modeled with patches, and a Lagrangian perspective is used to keep track of individuals' places of residence. The movement of pathogens in the air is modeled with linear diffusion and coupled to the SIR dynamics of each human population through an integral of the density of pathogens around the population patches. In the limit of fast diffusion pathogens, the method of matched asymptotic analysis is used to reduce the coupled PDE-ODE model to a nonlinear system of ODEs for the average density of pathogens in the air. The reduced system of ODEs is used to derive the basic reproduction number and the final size relation for the model. Numerical simulations of the full PDE-ODE model and the reduced system of ODEs are used to assess the impact of human mobility, together with the diffusion of pathogens on the dynamics of the disease. Results from the two models are consistent and show that human mobility significantly affects disease dynamics. In addition, we show that an increase in the diffusion rate of pathogen leads to a lower epidemic., (© 2022. The Author(s), under exclusive licence to Society for Mathematical Biology.)

Background: Close-contact rates are thought to be a driving force behind the transmission of many infectious respiratory diseases. Yet, contact rates and their relation to transmission and the impact of control measures, are seldom quantified. We quantify the response of contact rates, reported cases and transmission of COVID-19, to public health contact-restriction orders, and examine the associations among these three variables in the province of British Columbia, Canada., Methods: We derived time series data for contact rates, daily cases and transmission of COVID-19 from a social contacts survey, reported case counts and by fitting a transmission model to reported cases, respectively. We used segmented regression to investigate impacts of public health orders; Pearson correlation to determine associations between contact rates and transmission; and vector autoregressive modeling to quantify lagged associations between contacts rates, daily cases, and transmission., Results: Declines in contact rates and transmission occurred concurrently with the announcement of public health orders, whereas declines in cases showed a reporting delay of about 2 weeks. Contact rates were a significant driver of COVID-19 and explained roughly 19 and 20% of the variation in new cases and transmission, respectively. Interestingly, increases in COVID-19 transmission and cases were followed by reduced contact rates: overall, daily cases explained about 10% of the variation in subsequent contact rates., Conclusion: We showed that close-contact rates were a significant time-series driver of transmission and ultimately of reported cases of COVID-19 in British Columbia, Canada and that they varied in response to public health orders. Our results also suggest possible behavioral feedback, by which increased reported cases lead to reduced subsequent contact rates. Our findings help to explain and validate the commonly assumed, but rarely measured, response of close contact rates to public health guidelines and their impact on the dynamics of infectious diseases., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Ringa, Iyaniwura, David, Irvine, Adu, Spencer, Janjua and Otterstatter.)

The pandemic of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) took the world by surprise. Following the first outbreak of COVID-19 in December 2019, several models have been developed to study and understand its transmission dynamics. Although the spread of COVID-19 is being slowed down by vaccination and other interventions, there is still a need to have a clear understanding of the evolution of the pandemic across countries, states and communities. To this end, there is a need to have a clearer picture of the initial spread of the disease in different regions. In this project, we used a simple SEIR model and a Bayesian inference framework to estimate the basic reproduction number of COVID-19 across Africa. Our estimates vary between 1.98 (Sudan) and 9.66 (Mauritius), with a median of 3.67 (90% CrI: 3.31-4.12). The estimates provided in this paper will help to inform COVID-19 modeling in the respective countries/regions., Competing Interests: The authors have declared that no competing interests exist.

Estimates of the basic reproduction number ( R 0 ) for COVID-19 are particularly variable in the context of transmission within locations such as long-term healthcare (LTHC) facilities. We sought to characterize the heterogeneity of R 0 across known outbreaks within these facilities. We used a unique comprehensive dataset of all outbreaks that occurred within LTHC facilities in British Columbia, Canada as of 21 September 2020. We estimated R 0 in 18 LTHC outbreaks with a novel Bayesian hierarchical dynamic model of susceptible, exposed, infected and recovered individuals, incorporating heterogeneity of R 0 between facilities. We further compared these estimates to those obtained with standard methods that use the exponential growth rate and maximum likelihood. The total size of outbreaks varied dramatically, with range of attack rates 2%-86%. The Bayesian analysis provided an overall estimate of R 0 = 2.51 (90% credible interval 0.47-9.0), with individual facility estimates ranging between 0.56 and 9.17. Uncertainty in these estimates was more constrained than standard methods, particularly for smaller outbreaks informed by the population-level model. We further estimated that intervention led to 61% (52%-69%) of all potential cases being averted within the LTHC facilities, or 75% (68%-79%) when using a model with multi-level intervention effect. Understanding of transmission risks and impact of intervention are essential in planning during the ongoing global pandemic, particularly in high-risk environments such as LTHC facilities., (© 2022 The Authors.)

We developed an endemic model of COVID-19 to assess the impact of vaccination and immunity waning on the dynamics of the disease. Our model exhibits the phenomenon of backward bifurcation and bi-stability, where a stable disease-free equilibrium coexists with a stable endemic equilibrium. The epidemiological implication of this is that the control reproduction number being less than unity is no longer sufficient to guarantee disease eradication. We showed that this phenomenon could be eliminated by either increasing the vaccine efficacy or by reducing the disease transmission rate (adhering to non-pharmaceutical interventions). Furthermore, we numerically investigated the impacts of vaccination and waning of both vaccine-induced immunity and post-recovery immunity on the disease dynamics. Our simulation results show that the waning of vaccine-induced immunity has more effect on the disease dynamics relative to post-recovery immunity waning and suggests that more emphasis should be on reducing the waning of vaccine-induced immunity to eradicate COVID-19., Competing Interests: Competing interestThe authors declare that there is no competing interest., (© The Author(s), under exclusive licence to Springer Nature B.V. 2022.)

Adherence to public health policies such as the non-pharmaceutical interventions implemented against COVID-19 plays a major role in reducing infections and controlling the spread of the diseases. In addition, understanding the transmission dynamics of the disease is also important in order to make and implement efficient public health policies. In this paper, we developed an SEIR-type compartmental model to assess the impact of adherence to COVID-19 non-pharmaceutical interventions and indirect transmission on the dynamics of the disease. Our model considers both direct and indirect transmission routes and stratifies the population into two groups: those that adhere to COVID-19 non-pharmaceutical interventions (NPIs) and those that do not adhere to the NPIs. We compute the control reproduction number and the final epidemic size relation for our model and study the effect of different parameters of the model on these quantities. Our results show that there is a significant benefit in adhering to the COVID-19 NPIs.

Following successful non-pharmaceutical interventions (NPI) aiming to control COVID-19, many jurisdictions reopened their economies and borders. As little immunity had developed in most populations, re-establishing higher contact carried substantial risks, and therefore many locations began to see resurgence in COVID-19 cases. We present a Bayesian method to estimate the leeway to reopen, or alternatively the strength of change required to re-establish COVID-19 control, in a range of jurisdictions experiencing different COVID-19 epidemics. We estimated the timing and strength of initial control measures such as widespread distancing and compared the leeway jurisdictions had to reopen immediately after NPI measures to later estimates of leeway. Finally, we quantified risks associated with reopening and the likely burden of new cases due to introductions from other jurisdictions. We found widely varying leeway to reopen. After initial NPI measures took effect, some jurisdictions had substantial leeway (e.g., Japan, New Zealand, Germany) with > 0.99 probability that contact rates were below 80% of the threshold for epidemic growth. Others had little leeway (e.g., the United Kingdom, Washington State) and some had none (e.g., Sweden, California). For most such regions, increases in contact rate of 1.5-2 fold would have had high (> 0.7) probability of exceeding past peak sizes. Most jurisdictions experienced June-August trajectories consistent with our projections of contact rate increases of 1-2-fold. Under such relaxation scenarios for some regions, we projected up to ∼100 additional cases if just one case were imported per week over six weeks, even between jurisdictions with comparable COVID-19 risk. We provide an R package covidseir to enable jurisdictions to estimate leeway and forecast cases under different future contact patterns. Estimates of leeway can establish a quantitative basis for decisions about reopening. We recommend a cautious approach to reopening economies and borders, coupled with strong monitoring for changes in transmission., (Crown Copyright © 2021. Published by Elsevier B.V. All rights reserved.)

Importance: An effective vaccine against SARS-CoV-2 will reduce morbidity and mortality and allow substantial relaxation of physical distancing policies. However, the ability of a vaccine to prevent infection or disease depends critically on protecting older individuals, who are at highest risk of severe disease., Objective: We quantitatively estimated the relative benefits of COVID-19 vaccines, in terms of preventing infection and death, with a particular focus on effectiveness in elderly people., Design: We applied compartmental mathematical modelling to determine the relative effects of vaccines that block infection and onward transmission, and those that prevent severe disease. We assumed that vaccines showing high efficacy in adults would be deployed, and examined the effects of lower vaccine efficacy among the elderly population., Setting and Participants: Our mathematical model was calibrated to simulate the course of an epidemic among the entire population of British Columbia, Canada. Within our model, the population was structured by age and levels of contact., Main Outcome(s) and Measure(s): We assessed the effectiveness of possible vaccines in terms of the predicted number of infections within the entire population, and deaths among people aged 65 years and over., Results: In order to reduce the overall rate of infections in the population, high rates of deployment to all age groups will be critical. However, to substantially reduce mortality among people aged 65 years and over, a vaccine must directly protect a high proportion of people in that group., Conclusions and Relevance: Effective vaccines deployed to a large fraction of the population are projected to substantially reduce infection in an otherwise susceptible population. However, even if transmission were blocked highly effectively by vaccination of children and younger adults, overall mortality would not be substantially reduced unless the vaccine is also directly protective in elderly people. We strongly recommend: (i) the inclusion of people aged 65 years and over in future trials of COVID-19 vaccine candidates; (ii) careful monitoring of vaccine efficacy in older age groups following vaccination., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2021 Elsevier Ltd. All rights reserved.)

Extensive non-pharmaceutical and physical distancing measures are currently the primary interventions against coronavirus disease 2019 (COVID-19) worldwide. It is therefore urgent to estimate the impact such measures are having. We introduce a Bayesian epidemiological model in which a proportion of individuals are willing and able to participate in distancing, with the timing of distancing measures informed by survey data on attitudes to distancing and COVID-19. We fit our model to reported COVID-19 cases in British Columbia (BC), Canada, and five other jurisdictions, using an observation model that accounts for both underestimation and the delay between symptom onset and reporting. We estimated the impact that physical distancing (social distancing) has had on the contact rate and examined the projected impact of relaxing distancing measures. We found that, as of April 11 2020, distancing had a strong impact in BC, consistent with declines in reported cases and in hospitalization and intensive care unit numbers; individuals practising physical distancing experienced approximately 0.22 (0.11-0.34 90% CI [credible interval]) of their normal contact rate. The threshold above which prevalence was expected to grow was 0.55. We define the "contact ratio" to be the ratio of the estimated contact rate to the threshold rate at which cases are expected to grow; we estimated this contact ratio to be 0.40 (0.19-0.60) in BC. We developed an R package 'covidseir' to make our model available, and used it to quantify the impact of distancing in five additional jurisdictions. As of May 7, 2020, we estimated that New Zealand was well below its threshold value (contact ratio of 0.22 [0.11-0.34]), New York (0.60 [0.43-0.74]), Washington (0.84 [0.79-0.90]) and Florida (0.86 [0.76-0.96]) were progressively closer to theirs yet still below, but California (1.15 [1.07-1.23]) was above its threshold overall, with cases still rising. Accordingly, we found that BC, New Zealand, and New York may have had more room to relax distancing measures than the other jurisdictions, though this would need to be done cautiously and with total case volumes in mind. Our projections indicate that intermittent distancing measures-if sufficiently strong and robustly followed-could control COVID-19 transmission. This approach provides a useful tool for jurisdictions to monitor and assess current levels of distancing relative to their threshold, which will continue to be essential through subsequent waves of this pandemic., Competing Interests: The authors have declared that no competing interests exist.

We formulated and analyzed a class of coupled partial and ordinary differential equation (PDE-ODE) model to study the spread of airborne diseases. Our model describes human populations with patches and the movement of pathogens in the air with linear diffusion. The diffusing pathogens are coupled to the SIR dynamics of each population patch using an integro-differential equation. Susceptible individuals become infected at some rate whenever they are in contact with pathogens (indirect transmission), and the spread of infection in each patch depends on the density of pathogens around the patch. In the limit where the pathogens are diffusing fast, a matched asymptotic analysis is used to reduce the coupled PDE-ODE model into a nonlinear system of ODEs, which is then used to compute the basic reproduction number and final size relation for different scenarios. Numerical simulations of the reduced system of ODEs and the full PDE-ODE model are consistent, and they predict a decrease in the spread of infection as the diffusion rate of pathogens increases. Furthermore, we studied the effect of patch location on the spread of infections for the case of two population patches. Our model predicts higher infections when the patches are closer to each other.