Resolving the enigma of Iquitos and Manaus: A modeling analysis of multiple COVID-19 epidemic waves in two Amazonian cities.

The two nearby Amazonian cities of Iquitos and Manaus endured explosive COVID-19 epidemics and may well have suffered the world's highest infection and death rates over 2020, the first year of the pandemic. State-of-the-art epidemiological and modeling studies estimated that the populations of both cities came close to attaining herd immunity (>70% infected) at the termination of the first wave and were thus protected. This makes it difficult to explain the more deadly second wave of COVID-19 that struck again in Manaus just months later, simultaneous with the appearance of a new P.1 variant of concern, creating a catastrophe for the unprepared population. It was suggested that the second wave was driven by reinfections, but the episode has become controversial and an enigma in the history of the pandemic. We present a data-driven model of epidemic dynamics in Iquitos, which we also use to explain and model events in Manaus. By reverse engineering the multiple epidemic waves over 2 y in these two cities, the partially observed Markov process model inferred that the first wave left Manaus with a highly susceptible and vulnerable population (≈40% infected) open to invasion by P.1, in contrast to Iquitos (≈72% infected). The model reconstructed the full epidemic outbreak dynamics from mortality data by fitting a flexible time-varying reproductive number [Formula: see text] while estimating reinfection and impulsive immune evasion. The approach is currently highly relevant given the lack of tools available to assess these factors as new SARS-CoV-2 virus variants appear with different degrees of immune evasion.

Exploring the seasonality and optimal control strategy of HIV/AIDS epidemic in China: The impact of seasonal testing.

In this work, we investigate how the seasonal variation in the number of individuals who are tested for an HIV antibody in outpatient clinics affects the HIV transmission patterns in China, which has not been well studied. Based on the characteristics of outpatient testing data and reported cases, we establish a periodic infectious disease model to study the impact of seasonal testing on HIV transmission. The results indicate that the seasonal testing is a driving factor for the seasonality of new cases. We demonstrate the feasibility of ending the HIV/AIDS epidemic. We find that the diagnostic rates related to testing play a crucial role in controlling the size of the epidemic. Specifically, when considering minimizing both infected individuals and diagnostic rates, the level of attention paid to undiagnosed infected individuals is always positively correlated with the optimal diagnostic rates, while the optimal diagnostic rates are negatively correlated with the size of the epidemic at the terminal time.

Estimation of the serial interval of monkeypox during the early outbreak in 2022.

With increased transmissibility and novel transmission mode, monkeypox poses new threats to public health globally in the background of the ongoing COVID-19 pandemic. Estimates of the serial interval, a key epidemiological parameter of infectious disease transmission, could provide insights into the virus transmission risks. As of October 2022, little was known about the serial interval of monkeypox due to the lack of contact tracing data. In this study, public-available contact tracing data of global monkeypox cases were collected and 21 infector-infectee transmission pairs were identified. We proposed a statistical method applied to real-world observations to estimate the serial interval of the monkeypox. We estimated a mean serial interval of 5.6 days with the right truncation and sampling bias adjusted and calculated the reproduction number of 1.33 for the early monkeypox outbreaks at a global scale. Our findings provided a preliminary understanding of the transmission potentials of the current situation of monkeypox outbreaks. We highlighted the need for continuous surveillance of monkeypox for transmission risk assessment.

Effectiveness of the booster dose of inactivated COVID-19 vaccine against Omicron BA.5 infection: a matched cohort study of adult close contacts.

BACKGROUND: Although COVID-19 vaccines and their booster regimens protect against symptomatic infections and severe outcomes, there is limited evidence about their protection against asymptomatic and symptomatic infections in real-world settings, particularly when considering that the majority of SARS-CoV-2 Omicron infections were asymptomatic. We aimed to assess the effectiveness of the booster dose of inactivated vaccines in mainland China, i.e., Sinopharm (BBIBP-CorV) and Sinovac (CoronaVac), against Omicron infection in an Omicron BA.5 seeded epidemic. METHODS: Based on an infection-naive but highly vaccinated population in Urumqi, China, the study cohort comprised all 37,628 adults who had a contact history with individuals having SARS-CoV-2 infections, i.e., close contacts, between August 1 and September 7, 2022. To actively detect SARS-CoV-2 infections, RT-PCR tests were performed by local authorities on a daily basis for all close contacts, and a testing-positive status was considered a laboratory-confirmed outcome. The cohort of close contacts was matched at a ratio of 1:5 with the fully vaccinated (i.e., 2 doses) and booster vaccinated groups (i.e., 3 doses) according to sex, age strata, calendar date, and contact settings. Multivariate conditional logistic regression models were adopted to estimate the marginal effectiveness of the booster dose against Omicron BA.5 infection after adjusting for confounding variables. Subgroup analyses were performed to assess vaccine effectiveness (VE) in different strata of sex, age, the time lag from the last vaccine dose to exposure, and the vaccination status of the source case. Kaplan-Meier curves were employed to visualize the follow-up process and testing outcomes among different subgroups of the matched cohort. FINDINGS: Before matching, 37,099 adult close contacts were eligible for cohort enrolment. After matching, the 2-dose and 3-dose groups included 3317 and 16,051 contacts, and the proportions with Omicron infections were 1.03% and 0.62% among contacts in the 2-dose and 3-dose groups, respectively. We estimated that the adjusted effectiveness of the inactivated booster vaccine versus 2 doses against Omicron infection was 35.5% (95% CI 2.0, 57.5). The booster dose provided a higher level of protection, with an effectiveness of 60.2% (95% CI 22.8, 79.5) for 15-180 days after vaccination, but this VE decreased to 35.0% (95% CI 2.8, 56.5) after 180 days. Evidence for the protection of the booster dose was detected among young adults aged 18-39 years, but was not detected for those aged 40 years or older. INTERPRETATION: The receipt of the inactivated vaccine booster dose was associated with a significantly lower Omicron infection risk, and our findings confirmed the vaccine effectiveness (VE) of booster doses against Omicron BA.5 variants. Given the rapid evolution of SARS-CoV-2, we highlight the importance of continuously monitoring the protective performance of vaccines against the genetic variants of SARS-CoV-2, regardless of existing vaccine coverage.

Transmission risk assessment of HIV/AIDS epidemic resulting from sexual transmission in China, 2013-2017.

To estimate the risk of human immunodeficiency virus (HIV) infection through sexual transmission in China from 2013 to 2017 accurately, we divide the total population into three groups, namely, men who have sex with men (MSM) group; non-marital and commercial sex group: female sex workers (FSW) and their clients (FSWC); non-marital and non-commercial sex group: general women (GW) and general men (GM). First, the risks of HIV infection among men who have contacts with infected men or infected women decrease annually. Second, the number of contacts between susceptible FSWC and infected FSW per unit time is greater than that between susceptible MSM and infected MSM, and also greater than that between susceptible FSW and infected FSWC, which suggests that the intervention for commercial sex of heterosexual men should be strengthened. Third, the effective reproduction numbers of the MSM group and non-commercial sex group decrease annually, while the effective reproduction number of the commercial sex group decreases first, then increases, because the risk of women being infected by men decreases first, then increases. Additionally, the effective reproduction number of the commercial sex group exceeds that of the MSM group after 2015, which indicates that commercial sex contributes more and more to the HIV/AIDS epidemic.

Characterizing superspreading potential of infectious disease: Decomposition of individual transmissibility.

In the context of infectious disease transmission, high heterogeneity in individual infectiousness indicates that a few index cases can generate large numbers of secondary cases, a phenomenon commonly known as superspreading. The potential of disease superspreading can be characterized by describing the distribution of secondary cases (of each seed case) as a negative binomial (NB) distribution with the dispersion parameter, k. Based on the feature of NB distribution, there must be a proportion of individuals with individual reproduction number of almost 0, which appears restricted and unrealistic. To overcome this limitation, we generalized the compound structure of a Poisson rate and included an additional parameter, and divided the reproduction number into independent and additive fixed and variable components. Then, the secondary cases followed a Delaporte distribution. We demonstrated that the Delaporte distribution was important for understanding the characteristics of disease transmission, which generated new insights distinct from the NB model. By using real-world dataset, the Delaporte distribution provides improvements in describing the distributions of COVID-19 and SARS cases compared to the NB distribution. The model selection yielded increasing statistical power with larger sample sizes as well as conservative type I error in detecting the improvement in fitting with the likelihood ratio (LR) test. Numerical simulation revealed that the control strategy-making process may benefit from monitoring the transmission characteristics under the Delaporte framework. Our findings highlighted that for the COVID-19 pandemic, population-wide interventions may control disease transmission on a general scale before recommending the high-risk-specific control strategies.

Propensity score analysis with missing data using a multi-task neural network.

BACKGROUND: Propensity score analysis is increasingly used to control for confounding factors in observational studies. Unfortunately, unavoidable missing values make estimating propensity scores extremely challenging. We propose a new method for estimating propensity scores in data with missing values. MATERIALS AND METHODS: Both simulated and real-world datasets are used in our experiments. The simulated datasets were constructed under 2 scenarios, the presence (T = 1) and the absence (T = 0) of the true effect. The real-world dataset comes from LaLonde's employment training program. We construct missing data with varying degrees of missing rates under three missing mechanisms: MAR, MCAR, and MNAR. Then we compare MTNN with 2 other traditional methods in different scenarios. The experiments in each scenario were repeated 20,000 times. Our code is publicly available at https://github.com/ljwa2323/MTNN . RESULTS: Under the three missing mechanisms of MAR, MCAR and MNAR, the RMSE between the effect and the true effect estimated by our proposed method is the smallest in simulations and in real-world data. Furthermore, the standard deviation of the effect estimated by our method is the smallest. In situations where the missing rate is low, the estimation of our method is more accurate. CONCLUSIONS: MTNN can perform propensity score estimation and missing value filling at the same time through shared hidden layers and joint learning, which solves the dilemma of traditional methods and is very suitable for estimating true effects in samples with missing values. The method is expected to be broadly generalized and applied to real-world observational studies.

Twice evasions of Omicron variants explain the temporal patterns in six Asian and Oceanic countries.

BACKGROUND: The ongoing coronavirus 2019 (COVID-19) pandemic has emerged and caused multiple pandemic waves in the following six countries: India, Indonesia, Nepal, Malaysia, Bangladesh and Myanmar. Some of the countries have been much less studied in this devastating pandemic. This study aims to assess the impact of the Omicron variant in these six countries and estimate the infection fatality rate (IFR) and the reproduction number [Formula: see text] in these six South Asia, Southeast Asia and Oceania countries. METHODS: We propose a Susceptible-Vaccinated-Exposed-Infectious-Hospitalized-Death-Recovered model with a time-varying transmission rate [Formula: see text] to fit the multiple waves of the COVID-19 pandemic and to estimate the IFR and [Formula: see text] in the aforementioned six countries. The level of immune evasion and the intrinsic transmissibility advantage of the Omicron variant are also considered in this model. RESULTS: We fit our model to the reported deaths well. We estimate the IFR (in the range of 0.016 to 0.136%) and the reproduction number [Formula: see text] (in the range of 0 to 9) in the six countries. Multiple pandemic waves in each country were observed in our simulation results. CONCLUSIONS: The invasion of the Omicron variant caused the new pandemic waves in the six countries. The higher [Formula: see text] suggests the intrinsic transmissibility advantage of the Omicron variant. Our model simulation forecast implies that the Omicron pandemic wave may be mitigated due to the increasing immunized population and vaccine coverage.

Getting Jab or Regular Test: Observations from an Impulsive Epidemic COVID-19 Model.

Several safe and effective vaccines are available to prevent individuals from experiencing severe illness or death as a result of COVID-19. Widespread vaccination is widely regarded as a critical tool in the fight against the disease. However, some individuals may choose not to vaccinate due to vaccine hesitancy or other medical conditions. In some sectors, regular compulsory testing is required for such unvaccinated individuals. Interestingly, different sectors require testing at various frequencies, such as weekly or biweekly. As a result, it is essential to determine the optimal testing frequency and identify underlying factors. This study proposes a population-based model that can accommodate different personal decision choices, such as getting vaccinated or undergoing regular tests, as well as vaccine efficacies and uncertainties in epidemic transmission. The model, formulated as impulsive differential equations, uses time instants to represent the reporting date for the test result of an unvaccinated individual. By employing well-accepted indices to measure transmission risk, including the basic reproduction number, the peak time, the final size, and the number of severe infections, the study shows that an optimal testing frequency is highly sensitive to parameters involved in the transmission process, such as vaccine efficacy, disease transmission rate, test accuracy, and existing vaccination coverage. The testing frequency should be appropriately designed with the consideration of all these factors, as well as the control objectives measured by epidemiological quantities of great concern.

The immune evasion ability of Delta variant is comparable to that of Beta variant in South Africa.

BACKGROUND: The high immune evasion ability of SARS-COV-2 Omicron variant surprised the world and appears to be far stronger than any previous variant. Previous to Omicron it has been difficult to assess and compare immune evasion ability of different variants, including the Beta and Delta variants, because of the relatively small numbers of reinfections and because of the problems in correctly identifying reinfections in the population. This has led to different claims appearing in the literature. Thus we find claims of both high and low immune evasion for the Beta variant. Some findings have suggested that the Beta variant has a higher immune evasion ability than the Delta variant in South Africa, and others that it has a lower ability. METHOD: In this brief report, we re-analyse a unique dataset of variant-specific reinfection data and a simple model to correct for the infection attack rates of different variants. RESULT: We find that a model with the Delta variant having  an equal or higher immune evasion ability than Beta variant is compatible with the data. CONCLUSION: We conclude that the immune evasion ability of Beta variant is not stronger than Delta variant, and indeed, the immune evasion abilities of both variants are weak in South Africa.

Evaluation of Effectiveness of Global COVID-19 Vaccination Campaign.

To model estimated deaths averted by COVID-19 vaccines, we used state-of-the-art mathematical modeling, likelihood-based inference, and reported COVID-19 death and vaccination data. We estimated that >1.5 million deaths were averted in 12 countries. Our model can help assess effectiveness of the vaccination program, which is crucial for curbing the COVID-19 pandemic.

The non-pharmaceutical interventions may affect the advantage in transmission of mutated variants during epidemics: A conceptual model for COVID-19.

As the COVID-19 pandemic continues, genetic mutations in SARS-CoV-2 emerge, and some of them are found more contagious than the previously identified strains, acting as the major mechanism for many large-scale epidemics. The transmission advantage of mutated variants is widely believed as an innate biological feature that is difficult to be altered by artificial factors. In this study, we explore how non-pharmaceutical interventions (NPI) may affect transmission advantage. A two-strain compartmental epidemic model is proposed and simulated to investigate the biological mechanism of the relationships among different NPIs, the changes in transmissibility of each strain and transmission advantage. Although the NPIs are effective in flattening the epidemic curve, we demonstrate that NPIs probably lead to a decline in transmission advantage, which is likely to occur if the NPIs become intensive. Our findings uncover the mechanistic relationship between NPIs and transmission advantage dynamically, and highlight the important role of NPIs not only in controlling the intensity of epidemics but also in slowing or even containing the growth of the proportion of variants.

Differences in the superspreading potentials of COVID-19 across contact settings.

BACKGROUND: Superspreading events (SSEs) played a critical role in fueling the COVID-19 outbreaks. Although it is well-known that COVID-19 epidemics exhibited substantial superspreading potential, little is known about the risk of observing SSEs in different contact settings. In this study, we aimed to assess the potential of superspreading in different contact settings in Japan. METHOD: Transmission cluster data from Japan was collected between January and July 2020. Infector-infectee transmission pairs were constructed based on the contact tracing history. We fitted the data to negative binomial models to estimate the effective reproduction number (R) and dispersion parameter (k). Other epidemiological issues relating to the superspreading potential were also calculated. RESULTS: The overall estimated R and k are 0.561 (95% CrI: 0.496, 0.640) and 0.221 (95% CrI: 0.186, 0.262), respectively. The transmission in community, healthcare facilities and school manifest relatively higher superspreading potentials, compared to other contact settings. We inferred that 13.14% (95% CrI: 11.55%, 14.87%) of the most infectious cases generated 80% of the total transmission events. The probabilities of observing superspreading events for entire population and community, household, health care facilities, school, workplace contact settings are 1.75% (95% CrI: 1.57%, 1.99%), 0.49% (95% CrI: 0.22%, 1.18%), 0.07% (95% CrI: 0.06%, 0.08%), 0.67% (95% CrI: 0.31%, 1.21%), 0.33% (95% CrI: 0.13%, 0.94%), 0.32% (95% CrI: 0.21%, 0.60%), respectively. CONCLUSION: The different potentials of superspreading in contact settings highlighted the need to continuously monitoring the transmissibility accompanied with the dispersion parameter, to timely identify high risk settings favoring the occurrence of SSEs.

The Heterogeneous Severity of COVID-19 in African Countries: A Modeling Approach.

The COVID-19 pandemic has had a considerable impact on global health and economics. The impact in African countries has not been investigated thoroughly via fitting epidemic models to the reported COVID-19 deaths. We downloaded the data for the 12 most-affected countries with the highest cumulative COVID-19 deaths to estimate the time-varying basic reproductive number ([Formula: see text]) and infection attack rate. We develop a simple epidemic model and fitted it to reported COVID-19 deaths in 12 African countries using iterated filtering and allowing a flexible transmission rate. We observe high heterogeneity in the case-fatality rate across the countries, which may be due to different reporting or testing efforts. South Africa, Tunisia, and Libya were most affected, exhibiting a relatively higher [Formula: see text] and infection attack rate. Thus, to effectively control the spread of COVID-19 epidemics in Africa, there is a need to consider other mitigation strategies (such as improvements in socioeconomic well-being, healthcare systems, the water supply, and awareness campaigns).

Modelling the effects of the contaminated environments on tuberculosis in Jiangsu, China.

Tuberculosis (TB) is still an important public health issue in Jiangsu province, China. In this study, based on the TB transmission routes and the statistical data of TB cases, we formulate a novel TB epidemic model accounting for the effects of the contaminated environments on TB transmission dynamics. The value of this study lies in two aspects. Mathematically, we define the basic reproduction number, R0, and prove that R0 can be used to govern the threshold dynamics of the model. Epidemiologically, we find that the annual average R0 is 1.13,>1 and TB in Jiangsu is an endemic disease. Therefore, in order to control the TB in Jiangsu efficiently, we must decrease the virus shedding rate or increase the recovery rates, and increase the environmental clearance rate.

Shrinkage in serial intervals across transmission generations of COVID-19.

One of the key epidemiological characteristics that shape the transmission of coronavirus disease 2019 (COVID-19) is the serial interval (SI). Although SI is commonly considered following a probability distribution at a population scale, recent studies reported a slight shrinkage (or contraction) of the mean of effective SI across transmission generations or over time. Here, we develop a likelihood-based statistical inference framework with truncation to explore the change in SI across transmission generations after adjusting the impacts of case isolation. The COVID-19 contact tracing surveillance data in Hong Kong are used for exemplification. We find that for COVID-19, the mean of individual SI is likely to shrink with a factor at 0.72 per generation (95%CI: 0.54, 0.96) as the transmission generation increases, where a threshold may exist as the lower boundary of this shrinking process. We speculate that one of the probable explanations for the shrinkage in SI might be an outcome due to the competition among multiple candidate infectors within the same case cluster. Thus, the nonpharmaceutical interventive strategies are crucially important to block the transmission chains, and mitigate the COVID-19 epidemic.

Inferencing superspreading potential using zero-truncated negative binomial model: exemplification with COVID-19.

BACKGROUND: In infectious disease transmission dynamics, the high heterogeneity in individual infectiousness indicates that few index cases generate large numbers of secondary cases, which is commonly known as superspreading events. The heterogeneity in transmission can be measured by describing the distribution of the number of secondary cases as a negative binomial (NB) distribution with dispersion parameter, k. However, such inference framework usually neglects the under-ascertainment of sporadic cases, which are those without known epidemiological link and considered as independent clusters of size one, and this may potentially bias the estimates. METHODS: In this study, we adopt a zero-truncated likelihood-based framework to estimate k. We evaluate the estimation performance by using stochastic simulations, and compare it with the baseline non-truncated version. We exemplify the analytical framework with three contact tracing datasets of COVID-19. RESULTS: We demonstrate that the estimation bias exists when the under-ascertainment of index cases with 0 secondary case occurs, and the zero-truncated inference overcomes this problem and yields a less biased estimator of k. We find that the k of COVID-19 is inferred at 0.32 (95%CI: 0.15, 0.64), which appears slightly smaller than many previous estimates. We provide the simulation codes applying the inference framework in this study. CONCLUSIONS: The zero-truncated framework is recommended for less biased transmission heterogeneity estimates. These findings highlight the importance of individual-specific case management strategies to mitigate COVID-19 pandemic by lowering the transmission risks of potential super-spreaders with priority.

Attach importance of the bootstrap t test against Student's t test in clinical epidemiology: a demonstrative comparison using COVID-19 as an example.

Student's t test is valid for statistical inference under the normality assumption or asymptotically. By contrast, although the bootstrap t test was proposed in 1993, it is seldom adopted in medical research. We aim to demonstrate that the bootstrap t test outperforms Student's t test under normality in data. Using random data samples from normal distributions, we evaluated the testing performance, in terms of true-positive rate (TPR) and false-positive rate and diagnostic abilities, in terms of the area under the curve (AUC), of the bootstrap t test and Student's t test. We explore the AUC of both tests with varying sample size and coefficient of variation. We compare the testing outcomes using the COVID-19 serial interval (SI) data in Shenzhen and Hong Kong, China, for demonstration. With fixed TPR, the bootstrap t test maintained the equivalent accuracy in TPR, but significantly improved the true-negative rate from the Student's t test. With varying TPR, the diagnostic ability of bootstrap t test outperformed or equivalently performed as Student's t test in terms of the AUC. The equivalent performances are possible but rarely occur in practice. We find that the bootstrap t test outperforms by successfully detecting the difference in COVID-19 SI, which is defined as the time interval between consecutive transmission generations, due to sex and non-pharmaceutical interventions against the Student's t test. We demonstrated that the bootstrap t test outperforms Student's t test, and it is recommended to replace Student's t test in medical data analysis regardless of sample size.

A Zika Endemic Model for the Contribution of Multiple Transmission Routes.

Zika virus disease is a viral disease primarily transmitted to humans through the bite of infected female mosquitoes. Recent evidence indicates that the virus can also be sexually transmitted in hosts and vertically transmitted in vectors. In this paper, we propose a Zika model with three transmission routes, that is, vector-borne transmission between humans and mosquitoes, sexual transmission within humans and vertical transmission within mosquitoes. The basic reproduction number [Formula: see text] is computed and shown to be a sharp threshold quantity. Namely, the disease-free equilibrium is globally asymptotically stable as [Formula: see text], whereas there exists a unique endemic equilibrium which is globally asymptotically stable as [Formula: see text]. The relative contributions of each transmission route on the reproduction number, and the short- and long-term host infections are analyzed. Numerical simulations confirm that vectorial transmission contributes the most to the initial and subsequent transmission. The role of sexual transmission in the early phase of a Zika outbreak is greater than the long term, while vertical transmission is the opposite. Reducing mosquito bites is the most effective measure in lowering the risk of Zika virus infection.

Synchronized nonpharmaceutical interventions for the control of COVID-19.

The world is experiencing an ongoing pandemic of coronavirus disease-2019 (COVID-19), which is caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). In attempts to control the pandemic, a range of nonpharmaceutical interventions (NPIs) has been implemented worldwide. However, the effect of synchronized NPIs for the control of COVID-19 at temporal and spatial scales has not been well studied. Therefore, a meta-population model that incorporates essential nonlinear processes was constructed to uncover the transmission characteristics of SARS-CoV-2 and then assess the effectiveness of synchronized NPIs on COVID-19 dynamics in China. Regional synchronization of NPIs was observed in China, and it was found that a combination of synchronized NPIs (the travel restrictions, the social distancing and the infection isolation) prevented 93.7% of SARS-CoV-2 infections. The use of synchronized NPIs at the time of the Wuhan lockdown may have prevented as much as 38% of SARS-CoV-2 infections, compared with the unsynchronized scenario. The interconnectivity of the epicenter, the implementation time of synchronized NPIs, and the number of regions considered all affected the performance of synchronized NPIs. The results highlight the importance of using synchronized NPIs in high-risk regions for the control of COVID-19 and shed light on effective strategies for future pandemic responses. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1007/s11071-021-06505-0.