Substantial transmission of SARS-CoV-2 through casual contact in retail stores: Evidence from matched administrative microdata on card payments and testing.

This paper presents quasiexperimental evidence of Covid-19 transmission through casual contact between customers in retail stores. For a large sample of individuals in Denmark, we match card payment data, indicating exactly where and when each individual made purchases, with Covid-19 test data, indicating when each individual was tested and whether the test was positive. The resulting dataset identifies more than 100,000 instances where an infected individual made a purchase in a store and, in each instance, allows us to track the infection dynamics of other individuals who made purchases in the same store around the same time. We estimate transmissions by comparing the infection rate of exposed customers, who made a purchase within 5 min of an infected individual, and nonexposed customers, who made a purchase in the same store 16 to 30 min before. We find that exposure to an infected individual in a store increases the infection rate by around 0.12 percentage points (P < 0.001) between day 3 and day 7 after exposure. The estimates imply that transmissions in stores contributed around 0.04 to the reproduction number for the average infected individual and significantly more in the period where Omicron was the dominant variant.

Estimating the reproduction number and transmission heterogeneity from the size distribution of clusters of identical pathogen sequences.

Quantifying transmission intensity and heterogeneity is crucial to ascertain the threat posed by infectious diseases and inform the design of interventions. Methods that jointly estimate the reproduction number R and the dispersion parameter k have however mainly remained limited to the analysis of epidemiological clusters or contact tracing data, whose collection often proves difficult. Here, we show that clusters of identical sequences are imprinted by the pathogen offspring distribution, and we derive an analytical formula for the distribution of the size of these clusters. We develop and evaluate an inference framework to jointly estimate the reproduction number and the dispersion parameter from the size distribution of clusters of identical sequences. We then illustrate its application across a range of epidemiological situations. Finally, we develop a hypothesis testing framework relying on clusters of identical sequences to determine whether a given pathogen genetic subpopulation is associated with increased or reduced transmissibility. Our work provides tools to estimate the reproduction number and transmission heterogeneity from pathogen sequences without building a phylogenetic tree, thus making it easily scalable to large pathogen genome datasets.

Genotypic-phenotypic landscape computation based on first principle and deep learning.

The relationship between genotype and fitness is fundamental to evolution, but quantitatively mapping genotypes to fitness has remained challenging. We propose the Phenotypic-Embedding theorem (P-E theorem) that bridges genotype-phenotype through an encoder-decoder deep learning framework. Inspired by this, we proposed a more general first principle for correlating genotype-phenotype, and the P-E theorem provides a computable basis for the application of first principle. As an application example of the P-E theorem, we developed the Co-attention based Transformer model to bridge Genotype and Fitness model, a Transformer-based pre-train foundation model with downstream supervised fine-tuning that can accurately simulate the neutral evolution of viruses and predict immune escape mutations. Accordingly, following the calculation path of the P-E theorem, we accurately obtained the basic reproduction number (${R}_0$) of SARS-CoV-2 from first principles, quantitatively linked immune escape to viral fitness and plotted the genotype-fitness landscape. The theoretical system we established provides a general and interpretable method to construct genotype-phenotype landscapes, providing a new paradigm for studying theoretical and computational biology.

Phylodynamic signatures in the emergence of community-associated MRSA.

Community-associated, methicillin-resistant Staphylococcus aureus (MRSA) lineages have emerged in many geographically distinct regions around the world during the past 30 y. Here, we apply consistent phylodynamic methods across multiple community-associated MRSA lineages to describe and contrast their patterns of emergence and dissemination. We generated whole-genome sequencing data for the Australian sequence type (ST) ST93-MRSA-IV from remote communities in Far North Queensland and Papua New Guinea, and the Bengal Bay ST772-MRSA-V clone from metropolitan communities in Pakistan. Increases in the effective reproduction number (Re) and sustained transmission (Re > 1) coincided with spread of progenitor methicillin-susceptible S. aureus (MSSA) in remote northern Australian populations, dissemination of the ST93-MRSA-IV genotype into population centers on the Australian East Coast, and subsequent importation into the highlands of Papua New Guinea and Far North Queensland. Applying the same phylodynamic methods to existing lineage datasets, we identified common signatures of epidemic growth in the emergence and epidemiological trajectory of community-associated S. aureus lineages from America, Asia, Australasia, and Europe. Surges in Re were observed at the divergence of antibiotic-resistant strains, coinciding with their establishment in regional population centers. Epidemic growth was also observed among drug-resistant MSSA clades in Africa and northern Australia. Our data suggest that the emergence of community-associated MRSA in the late 20th century was driven by a combination of antibiotic-resistant genotypes and host epidemiology, leading to abrupt changes in lineage-wide transmission dynamics and sustained transmission in regional population centers.

Time Scales of Human Mpox Transmission in The Netherlands.

Mpox has spread rapidly to many countries in nonendemic regions. After reviewing detailed exposure histories of 109 pairs of mpox cases in the Netherlands, we identified 34 pairs where transmission was likely and the infectee reported a single potential infector with a mean serial interval of 10.1 days (95% credible interval, 6.6-14.7 days). Further investigation into pairs from 1 regional public health service revealed that presymptomatic transmission may have occurred in 5 of 18 pairs. These findings emphasize that precaution remains key, regardless of the presence of recognizable symptoms of mpox.

Characteristics of the Sexual Networks of Men Who Have Sex With Men in Montréal, Toronto, and Vancouver: Insights from Canada's 2022 Mpox Outbreak.

BACKGROUND: The 2022-2023 global mpox outbreak disproportionately affected gay, bisexual, and other men who have sex with men (GBM). We investigated differences in GBM's sexual partner distributions across Canada's 3 largest cities and over time, and how they shaped transmission. METHODS: The Engage Cohort Study (2017-2023) recruited GBM via respondent-driven sampling in Montréal, Toronto, and Vancouver (n = 2449). We compared reported sexual partner distributions across cities and periods: before COVID-19 (2017-2019), pandemic (2020-2021), and after lifting of restrictions (2021-2023). We used Bayesian regression and poststratification to model partner distributions. We estimated mpox's basic reproduction number (R0) using a risk-stratified compartmental model. RESULTS: Pre-COVID-19 pandemic distributions were comparable: fitted average partners (past 6 months) were 10.4 (95% credible interval: 9.4-11.5) in Montréal, 13.1 (11.3-15.1) in Toronto, and 10.7 (9.5-12.1) in Vancouver. Sexual activity decreased during the pandemic and increased after lifting of restrictions, but remained below prepandemic levels. Based on reported cases, we estimated R0 of 2.4 to 2.7 and similar cumulative incidences (0.7%-0.9%) across cities. CONCLUSIONS: Similar sexual partner distributions may explain comparable R0 and cumulative incidence across cities. With potential for further recovery in sexual activity, mpox vaccination and surveillance strategies should be maintained.

SARS-CoV-2 Transmission in Alberta, British Columbia, and Ontario, Canada, January 2020-January 2022.

We estimated COVID-19 transmission potential and case burden by variant type in Alberta, British Columbia, and Ontario, Canada, during January 23, 2020-January 27, 2022; we also estimated the effectiveness of public health interventions to reduce transmission. We estimated time-varying reproduction number (Rt) over 7-day sliding windows and nonoverlapping time-windows determined by timing of policy changes. We calculated incidence rate ratios (IRRs) for each variant and compared rates to determine differences in burden among provinces. Rt corresponding with emergence of the Delta variant increased in all 3 provinces; British Columbia had the largest increase, 43.85% (95% credible interval [CrI] 40.71%-46.84%). Across the study period, IRR was highest for Omicron (8.74 [95% CrI 8.71-8.77]) and burden highest in Alberta (IRR 1.80 [95% CrI 1.79-1.81]). Initiating public health interventions was associated with lower Rt and relaxing restrictions and emergence of new variants associated with increases in Rt.

Insights into COVID-19 epidemiology and control from temporal changes in serial interval distributions in Hong Kong.

The serial interval distribution is used to approximate the generation time distribution, an essential parameter to infer the transmissibility (${R}_t$) of an epidemic. However, serial interval distributions may change as an epidemic progresses. We examined detailed contact tracing data on laboratory-confirmed cases of COVID-19 in Hong Kong during the five waves from January 2020 to July 2022. We reconstructed the transmission pairs and estimated time-varying effective serial interval distributions and factors associated with longer or shorter intervals. Finally, we assessed the biases in estimating transmissibility using constant serial interval distributions. We found clear temporal changes in mean serial interval estimates within each epidemic wave studied and across waves, with mean serial intervals ranged from 5.5 days (95% CrI: 4.4, 6.6) to 2.7 (95% CrI: 2.2, 3.2) days. The mean serial intervals shortened or lengthened over time, which were found to be closely associated with the temporal variation in COVID-19 case profiles and public health and social measures and could lead to the biases in predicting ${R}_t$. Accounting for the impact of these factors, the time-varying quantification of serial interval distributions could lead to improved estimation of ${R}_t$, and provide additional insights into the impact of public health measures on transmission.

Season of death, pathogen persistence and wildlife behaviour alter number of anthrax secondary infections from environmental reservoirs.

An important part of infectious disease management is predicting factors that influence disease outbreaks, such as R, the number of secondary infections arising from an infected individual. Estimating R is particularly challenging for environmentally transmitted pathogens given time lags between cases and subsequent infections. Here, we calculated R for Bacillus anthracis infections arising from anthrax carcass sites in Etosha National Park, Namibia. Combining host behavioural data, pathogen concentrations and simulation models, we show that R is spatially and temporally variable, driven by spore concentrations at death, host visitation rates and early preference for foraging at infectious sites. While spores were detected up to a decade after death, most secondary infections occurred within 2 years. Transmission simulations under scenarios combining site infectiousness and host exposure risk under different environmental conditions led to dramatically different outbreak dynamics, from pathogen extinction (R < 1) to explosive outbreaks (R > 10). These transmission heterogeneities may explain variation in anthrax outbreak dynamics observed globally, and more generally, the critical importance of environmental variation underlying host-pathogen interactions. Notably, our approach allowed us to estimate the lethal dose of a highly virulent pathogen non-invasively from observational studies and epidemiological data, useful when experiments on wildlife are undesirable or impractical.

Supercritical and homogenous transmission of monkeypox in the capital of China.

Starting from May 31, 2023, the local transmission of monkeypox (Mpox) in mainland China began in Beijing. Till now, the transmission characteristics have not been explored. Based on the daily Mpox incidence data in the first 3 weeks of Beijing (from May 31 to June 21, 2023), we employed the instant-individual heterogeneity transmission model to simultaneously calculate the effective reproduction number (Re ) and the degree of heterogeneity (k) of the Beijing epidemic. We additionally simulated the monthly infection size in Beijing from July to November and compared with the reported data to project subsequent transmission dynamics. We estimated Re to be 1.68 (95% highest posterior density [HPD]: 1.12-2.41), and k to be 2.57 [95% HPD: 0.54-83.88], suggesting the transmission of Mpox in Beijing was supercritical and didn't have considerable transmission heterogeneity. We projected that Re fell in the range of 0.95-1.0 from July to November, highlighting more efforts needed to further reduce the Mpox transmissibility. Our findings revealed supercritical and homogeneous transmission of the Mpox epidemic in Beijing. Our results could serve as a reference for understanding and predicting the ongoing Mpox transmission in other regions of China and evaluating the effect of control measures.

Who acquires infection from whom? A sensitivity analysis of transmission dynamics during the early phase of the COVID-19 pandemic in Belgium.

Age-related heterogeneity in a host population, whether due to how individuals mix and contact each other, the nature of host-pathogen interactions defining epidemiological parameters, or demographics, is crucial in studying infectious disease dynamics. Compartmental models represent a popular approach to address the problem, dividing the population of interest into a discrete and finite number of states depending on, for example, individuals' age and stage of infection. We study the corresponding linearised system whose operator, in the context of a discrete-time model, equates to a square matrix known as the next generation matrix. Performing formal perturbation analysis of the entries of the aforementioned matrix, we derive indices to quantify the age-specific variation of its dominant eigenvalue (i.e., the reproduction number) and explore the relevant epidemiological information we can derive from the eigenstructure of the matrix. The resulting method enables the assessment of the impact of age-related population heterogeneity on virus transmission. In particular, starting from an age-structured SEIR model, we demonstrate the use of this approach for COVID-19 dynamics in Belgium. We analyse the early stages of the SARS-CoV-2 spread, with particular attention to the pre-pandemic framework and the lockdown lifting phase initiated as of May 2020. Our results, influenced by our assumption on age-specific susceptibility and infectiousness, support the hypothesis that transmission was only influenced to a small extent by children in the age group [0,18) and adults over 60 years of age during the early phases of the pandemic and up to the end of July 2020.

Modelling the impact of human behavior using a two-layer Watts-Strogatz network for transmission and control of Mpox.

PURPOSE: This study aims to evaluate the effectiveness of mitigation strategies and analyze the impact of human behavior on the transmission of Mpox. The results can provide guidance to public health authorities on comprehensive prevention and control for the new Mpox virus strain in the Democratic Republic of Congo as of December 2023. METHODS: We develop a two-layer Watts-Strogatz network model. The basic reproduction number is calculated using the next-generation matrix approach. Markov chain Monte Carlo (MCMC) optimization algorithm is used to fit Mpox cases in Canada into the network model. Numerical simulations are used to assess the impact of mitigation strategies and human behavior on the final epidemic size. RESULTS: Our results show that the contact transmission rate of low-risk groups and susceptible humans increases when the contact transmission rate of high-risk groups and susceptible humans is controlled as the Mpox epidemic spreads. The contact transmission rate of high-risk groups after May 18, 2022, is approximately 20% lower than that before May 18, 2022. Our findings indicate a positive correlation between the basic reproduction number and the level of heterogeneity in human contacts, with the basic reproduction number estimated at 2.3475 (95% CI: 0.0749-6.9084). Reducing the average number of sexual contacts to two per week effectively reduces the reproduction number to below one. CONCLUSION: We need to pay attention to the re-emergence of the epidemics caused by low-risk groups when an outbreak dominated by high-risk groups is under control. Numerical simulations show that reducing the average number of sexual contacts to two per week is effective in slowing down the rapid spread of the epidemic. Our findings offer guidance for the public health authorities of the Democratic Republic of Congo in developing effective mitigation strategies.

Mathematical Assessment of the Role of Human Behavior Changes on SARS-CoV-2 Transmission Dynamics in the United States.

The COVID-19 pandemic has not only presented a major global public health and socio-economic crisis, but has also significantly impacted human behavior towards adherence (or lack thereof) to public health intervention and mitigation measures implemented in communities worldwide. This study is based on the use of mathematical modeling approaches to assess the extent to which SARS-CoV-2 transmission dynamics is impacted by population-level changes of human behavior due to factors such as (a) the severity of transmission (such as disease-induced mortality and level of symptomatic transmission), (b) fatigue due to the implementation of mitigation interventions measures (e.g., lockdowns) over a long (extended) period of time, (c) social peer-pressure, among others. A novel behavior-epidemiology model, which takes the form of a deterministic system of nonlinear differential equations, is developed and fitted using observed cumulative SARS-CoV-2 mortality data during the first wave in the United States. The model fits the observed data, as well as makes a more accurate prediction of the observed daily SARS-CoV-2 mortality during the first wave (March 2020-June 2020), in comparison to the equivalent model which does not explicitly account for changes in human behavior. This study suggests that, as more newly-infected individuals become asymptomatically-infectious, the overall level of positive behavior change can be expected to significantly decrease (while new cases may rise, particularly if asymptomatic individuals have higher contact rate, in comparison to symptomatic individuals).

Prevention and control of Ebola virus transmission: mathematical modelling and data fitting.

The Ebola virus disease (EVD) has been endemic since 1976, and the case fatality rate is extremely high. EVD is spread by infected animals, symptomatic individuals, dead bodies, and contaminated environment. In this paper, we formulate an EVD model with four transmission modes and a time delay describing the incubation period. Through dynamical analysis, we verify the importance of blocking the infection source of infected animals. We get the basic reproduction number without considering the infection source of infected animals. And, it is proven that the model has a globally attractive disease-free equilibrium when the basic reproduction number is less than unity; the disease eventually becomes endemic when the basic reproduction number is greater than unity. Taking the EVD epidemic in Sierra Leone in 2014-2016 as an example, we complete the data fitting by combining the effect of the media to obtain the unknown parameters, the basic reproduction number and its time-varying reproduction number. It is shown by parameter sensitivity analysis that the contact rate and the removal rate of infected group have the greatest influence on the prevalence of the disease. And, the disease-controlling thresholds of these two parameters are obtained. In addition, according to the existing vaccination strategy, only the inoculation ratio in high-risk areas is greater than 0.4, the effective reproduction number can be less than unity. And, the earlier the vaccination time, the greater the inoculation ratio, and the faster the disease can be controlled.

Modeling insect growth regulators for pest management.

Insect growth regulators (IGRs) have been developed as effective control measures against harmful insect pests to disrupt their normal development. This study is to propose a mathematical model to evaluate the cost-effectiveness of IGRs for pest management. The key features of the model include the temperature-dependent growth of insects and realistic impulsive IGRs releasing strategies. The impulsive releases are carefully modeled by counting the number of implements during an insect's temperature-dependent development duration, which introduces a surviving probability determined by a product of terms corresponding to each release. Dynamical behavior of the model is illustrated through dynamical system analysis and a threshold-type result is established in terms of the net reproduction number. Further numerical simulations are performed to quantitatively evaluate the effectiveness of IGRs to control populations of harmful insect pests. It is interesting to observe that the time-changing environment plays an important role in determining an optimal pest control scheme with appropriate release frequencies and time instants.

Transmission dynamics of a reaction-advection-diffusion dengue fever model with seasonal developmental durations and intrinsic incubation periods.

In this paper, we propose a reaction-advection-diffusion dengue fever model with seasonal developmental durations and intrinsic incubation periods. Firstly, we establish the well-posedness of the model. Secondly, we define the basic reproduction number ℜ 0 for this model and show that ℜ 0 is a threshold parameter: if ℜ 0 < 1 , then the disease-free periodic solution is globally attractive; if ℜ 0 > 1 , the system is uniformly persistent. Thirdly, we study the global attractivity of the positive steady state when the spatial environment is homogeneous and the advection of mosquitoes is ignored. As an example, we use the model to investigate the dengue fever transmission case in Guangdong Province, China, and explore the impact of model parameters on ℜ 0 . Our findings indicate that ignoring seasonality may underestimate ℜ 0 . Additionally, the spatial heterogeneity of transmission may increase the risk of disease transmission, while the increase of seasonal developmental durations, intrinsic incubation periods and advection rates can all reduce the risk of disease transmission.

Threshold dynamics of a reaction-advection-diffusion schistosomiasis epidemic model with seasonality and spatial heterogeneity.

Most water-borne disease models ignore the advection of water flows in order to simplify the mathematical analysis and numerical computation. However, advection can play an important role in determining the disease transmission dynamics. In this paper, we investigate the long-term dynamics of a periodic reaction-advection-diffusion schistosomiasis model and explore the joint impact of advection, seasonality and spatial heterogeneity on the transmission of the disease. We derive the basic reproduction number R 0 and show that the disease-free periodic solution is globally attractive when R 0 < 1 whereas there is a positive endemic periodic solution and the system is uniformly persistent in a special case when R 0 > 1 . Moreover, we find that R 0 is a decreasing function of the advection coefficients which offers insights into why schistosomiasis is more serious in regions with slow water flows.

Modelling phytoplankton-virus interactions: phytoplankton blooms and lytic virus transmission.

A dynamic reaction-diffusion model of four variables is proposed to describe the spread of lytic viruses among phytoplankton in a poorly mixed aquatic environment. The basic ecological reproductive index for phytoplankton invasion and the basic reproduction number for virus transmission are derived to characterize the phytoplankton growth and virus transmission dynamics. The theoretical and numerical results from the model show that the spread of lytic viruses effectively controls phytoplankton blooms. This validates the observations and experimental results of Emiliana huxleyi-lytic virus interactions. The studies also indicate that the lytic virus transmission cannot occur in a low-light or oligotrophic aquatic environment.

Evaluation of age-structured vaccination strategies for curbing the disease spread.

Age structure is one of the crucial factors in characterizing the heterogeneous epidemic transmission. Vaccination is regarded as an effective control measure for prevention and control epidemics. Due to the shortage of vaccine capacity during the outbreak of epidemics, how to design vaccination policy has become an urgent issue in suppressing the disease transmission. In this paper, we make an effort to propose an age-structured SVEIHR model with the disease-caused death to take account of dynamics of age-related vaccination policy for better understanding disease spread and control. We present an explicit expression of the basic reproduction number R 0 , which determines whether or not the disease persists, and then establish the existence and stability of endemic equilibria under certain conditions. Numerical simulations are illustrated to show that the age-related vaccination policy has a tremendous influence on curbing the disease transmission. Especially, vaccination of people over 65 is better than for people aged 21-65 in terms of rapid eradication of the disease in Italy.

Final epidemic size of a two-community SIR model with asymmetric coupling.

Communities are commonly not isolated but interact asymmetrically with each other, allowing the propagation of infectious diseases within the same community and between different communities. To reveal the impact of asymmetrical interactions and contact heterogeneity on disease transmission, we formulate a two-community SIR epidemic model, in which each community has its contact structure while communication between communities occurs through temporary commuters. We derive an explicit formula for the basic reproduction number R 0 , give an implicit equation for the final epidemic size z, and analyze the relationship between them. Unlike the typical positive correlation between R 0 and z in the classic SIR model, we find a negatively correlated relationship between counterparts of our model deviating from homogeneous populations. Moreover, we investigate the impact of asymmetric coupling mechanisms on R 0 . The results suggest that, in scenarios with restricted movement of susceptible individuals within a community, R 0 does not follow a simple monotonous relationship, indicating that an unbending decrease in the movement of susceptible individuals may increase R 0 . We further demonstrate that network contacts within communities have a greater effect on R 0 than casual contacts between communities. Finally, we develop an epidemic model without restriction on the movement of susceptible individuals, and the numerical simulations suggest that the increase in human flow between communities leads to a larger R 0 .