Bayesian workflow for time-varying transmission in stratified compartmental infectious disease transmission models.

Compartmental models that describe infectious disease transmission across subpopulations are central for assessing the impact of non-pharmaceutical interventions, behavioral changes and seasonal effects on the spread of respiratory infections. We present a Bayesian workflow for such models, including four features: (1) an adjustment for incomplete case ascertainment, (2) an adequate sampling distribution of laboratory-confirmed cases, (3) a flexible, time-varying transmission rate, and (4) a stratification by age group. Within the workflow, we benchmarked the performance of various implementations of two of these features (2 and 3). For the second feature, we used SARS-CoV-2 data from the canton of Geneva (Switzerland) and found that a quasi-Poisson distribution is the most suitable sampling distribution for describing the overdispersion in the observed laboratory-confirmed cases. For the third feature, we implemented three methods: Brownian motion, B-splines, and approximate Gaussian processes (aGP). We compared their performance in terms of the number of effective samples per second, and the error and sharpness in estimating the time-varying transmission rate over a selection of ordinary differential equation solvers and tuning parameters, using simulated seroprevalence and laboratory-confirmed case data. Even though all methods could recover the time-varying dynamics in the transmission rate accurately, we found that B-splines perform up to four and ten times faster than Brownian motion and aGPs, respectively. We validated the B-spline model with simulated age-stratified data. We applied this model to 2020 laboratory-confirmed SARS-CoV-2 cases and two seroprevalence studies from the canton of Geneva. This resulted in detailed estimates of the transmission rate over time and the case ascertainment. Our results illustrate the potential of the presented workflow including stratified transmission to estimate age-specific epidemiological parameters. The workflow is freely available in the R package HETTMO, and can be easily adapted and applied to other infectious diseases.

Surveillance of SARS-CoV-2 prevalence from repeated pooled testing: application to Swiss routine data.

Surveillance of SARS-CoV-2 through reported positive RT-PCR tests is biased due to non-random testing. Prevalence estimation in population-based samples corrects for this bias. Within this context, the pooled testing design offers many advantages, but several challenges remain with regards to the analysis of such data. We developed a Bayesian model aimed at estimating the prevalence of infection from repeated pooled testing data while (i) correcting for test sensitivity; (ii) propagating the uncertainty in test sensitivity; and (iii) including correlation over time and space. We validated the model in simulated scenarios, showing that the model is reliable when the sample size is at least 500, the pool size below 20, and the true prevalence below 5%. We applied the model to 1.49 million pooled tests collected in Switzerland in 2021-2022 in schools, care centres, and workplaces. We identified similar dynamics in all three settings, with prevalence peaking at 4-5% during winter 2022. We also identified differences across regions. Prevalence estimates in schools were correlated with reported cases, hospitalizations, and deaths (coefficient 0.84 to 0.90). We conclude that in many practical situations, the pooled test design is a reliable and affordable alternative for the surveillance of SARS-CoV-2 and other viruses.

Health of singleton neonates in Switzerland through time and crises: a cross-sectional study at the population level, 2007-2022.

BACKGROUND: Being exposed to crises during pregnancy can affect maternal health through stress exposure, which can in return impact neonatal health. We investigated temporal trends in neonatal outcomes in Switzerland between 2007 and 2022 and their variations depending on exposure to the economic crisis of 2008, the flu pandemic of 2009, heatwaves (2015 and 2018) and the COVID-19 pandemic. METHODS: Using individual cross-sectional data encompassing all births occurring in Switzerland at the monthly level (2007-2022), we analysed changes in birth weight and in the rates of preterm birth (PTB) and stillbirth through time with generalized additive models. We assessed whether the intensity or length of crisis exposure was associated with variations in these outcomes. Furthermore, we explored effects of exposure depending on trimesters of pregnancy. RESULTS: Over 1.2 million singleton births were included in our analyses. While birth weight and the rate of stillbirth have remained stable since 2007, the rate of PTB has declined by one percentage point. Exposure to the crises led to different results, but effect sizes were overall small. Exposure to COVID-19, irrespective of the pregnancy trimester, was associated with a higher birth weight (+12 grams [95% confidence interval (CI) 5.5 to 17.9 grams]). Being exposed to COVID-19 during the last trimester was associated with an increased risk of stillbirth (odds ratio 1.24 [95%CI 1.02 to 1.50]). Exposure to the 2008 economic crisis during pregnancy was not associated with any changes in neonatal health outcomes, while heatwave effect was difficult to interpret. CONCLUSION: Overall, maternal and neonatal health demonstrated resilience to the economic crisis and to the COVID-19 pandemic in a high-income country like Switzerland. However, the effect of exposure to the COVID-19 pandemic is dual, and the negative impact of maternal infection on pregnancy is well-documented. Stress exposure and economic constraint may also have had adverse effects among the most vulnerable subgroups of Switzerland. To investigate better the impact of heatwave exposure on neonatal health, weekly or daily-level data is needed, instead of monthly-level data.

Area-level excess mortality in times of COVID-19 in Switzerland: geographical, socioeconomic and political determinants.

The coronavirus disease 2019 (COVID-19)-related excess mortality in Switzerland is well documented, but no study examined mortality at the small-area level. We analysed excess mortality in 2020 for 2141 Swiss municipalities using a Bayesian spatiotemporal model fitted to 2011-19 data. Areas most affected included the Ticino, the Romandie and the Northeast. Rural areas, municipalities within cross-border labour markets, of lower socioeconomic position and with less support for control measures in the popular vote on the COVID-19 Act had greater excess mortality. Particularly vulnerable municipalities require special efforts to mitigate the impact of pandemics.

Projecting the development of antimicrobial resistance in Neisseria gonorrhoeae from antimicrobial surveillance data: a mathematical modelling study.

BACKGROUND: The World Health Organization recommends changing the first-line antimicrobial treatment for gonorrhoea when ≥ 5% of Neisseria gonorrhoeae cases fail treatment or are resistant. Susceptibility to ceftriaxone, the last remaining treatment option has been decreasing in many countries. We used antimicrobial resistance surveillance data and developed mathematical models to project the time to reach the 5% threshold for resistance to first-line antimicrobials used for N. gonorrhoeae. METHODS: We used data from the Gonococcal Resistance to Antimicrobials Surveillance Programme (GRASP) in England and Wales from 2000-2018 about minimum inhibitory concentrations (MIC) for ciprofloxacin, azithromycin, cefixime and ceftriaxone and antimicrobial treatment in two groups, heterosexual men and women (HMW) and men who have sex with men (MSM). We developed two susceptible-infected-susceptible models to fit these data and produce projections of the proportion of resistance until 2030. The single-step model represents the situation in which a single mutation results in antimicrobial resistance. In the multi-step model, the sequential accumulation of resistance mutations is reflected by changes in the MIC distribution. RESULTS: The single-step model described resistance to ciprofloxacin well. Both single-step and multi-step models could describe azithromycin and cefixime resistance, with projected resistance levels higher with the multi-step than the single step model. For ceftriaxone, with very few observed cases of full resistance, the multi-step model was needed to describe long-term dynamics of resistance. Extrapolating from the observed upward drift in MIC values, the multi-step model projected ≥ 5% resistance to ceftriaxone could be reached by 2030, based on treatment pressure alone. Ceftriaxone resistance was projected to rise to 13.2% (95% credible interval [CrI]: 0.7-44.8%) among HMW and 19.6% (95%CrI: 2.6-54.4%) among MSM by 2030. CONCLUSIONS: New first-line antimicrobials for gonorrhoea treatment are needed. In the meantime, public health authorities should strengthen surveillance for AMR in N. gonorrhoeae and implement strategies for continued antimicrobial stewardship. Our models show the utility of long-term representative surveillance of gonococcal antimicrobial susceptibility data and can be adapted for use in, and for comparison with, other countries.

Socio-demographic characteristics associated with COVID-19 vaccination uptake in Switzerland: longitudinal analysis of the CoMix study.

BACKGROUND: Vaccination is an effective strategy to reduce morbidity and mortality from coronavirus disease 2019 (COVID-19). However, the uptake of COVID-19 vaccination has varied across and within countries. Switzerland has had lower levels of COVID-19 vaccination uptake in the general population than many other high-income countries. Understanding the socio-demographic factors associated with vaccination uptake can help to inform future vaccination strategies to increase uptake. METHODS: We conducted a longitudinal online survey in the Swiss population, consisting of six survey waves from June to September 2021. Participants provided information on socio-demographic characteristics, history of testing for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), social contacts, willingness to be vaccinated, and vaccination status. We used a multivariable Poisson regression model to estimate the adjusted rate ratio (aRR) and 95% confidence intervals (CI) of COVID-19 vaccine uptake. RESULTS: We recorded 6,758 observations from 1,884 adults. For the regression analysis, we included 3,513 observations from 1,883 participants. By September 2021, 600 (75%) of 806 study participants had received at least one vaccine dose. Participants who were older, male, and students, had a higher educational level, household income, and number of social contacts, and lived in a household with a medically vulnerable person were more likely to have received at least one vaccine dose. Female participants, those who lived in rural areas and smaller households, and people who perceived COVID-19 measures as being too strict were less likely to be vaccinated. We found no significant association between previous SARS-CoV-2 infections and vaccination uptake. CONCLUSIONS: Our results suggest that socio-demographic factors as well as individual behaviours and attitudes played an important role in COVID-19 vaccination uptake in Switzerland. Therefore, appropriate communication with the public is needed to ensure that public health interventions are accepted and implemented by the population. Tailored COVID-19 vaccination strategies in Switzerland that aim to improve uptake should target specific subgroups such as women, people from rural areas or people with lower socio-demographic status.

Historically High Excess Mortality During the COVID-19 Pandemic in Switzerland, Sweden, and Spain.

BACKGROUND: Excess mortality quantifies the overall mortality impact of a pandemic. Mortality data have been accessible for many countries in recent decades, but few continuous data have been available for longer periods. OBJECTIVE: To assess the historical dimension of the COVID-19 pandemic in 2020 for 3 countries with reliable death count data over an uninterrupted span of more than 100 years. DESIGN: Observational study. SETTING: Switzerland, Sweden, and Spain, which were militarily neutral and not involved in combat during either world war and have not been affected by significant changes in their territory since the end of the 19th century. PARTICIPANTS: Complete populations of these 3 countries. MEASUREMENTS: Continuous series of recorded deaths (from all causes) by month from the earliest available year (1877 for Switzerland, 1851 for Sweden, and 1908 for Spain) were jointly modeled with annual age group-specific death and total population counts using negative binomial and multinomial models, which accounted for temporal trends and seasonal variability of prepandemic years. The aim was to estimate the expected number of deaths in a pandemic year for a nonpandemic scenario and the difference in observed and expected deaths aggregated over the year. RESULTS: In 2020, the number of excess deaths recorded per 100 000 persons was 100 (95% credible interval [CrI], 60 to 135) for Switzerland, 75 (CrI, 40 to 105) for Sweden, and 155 (CrI, 110 to 195) for Spain. In 1918, excess mortality was 6 to 7 times higher. In all 3 countries, the peaks of monthly excess mortality in 2020 were greater than most monthly excess mortality since 1918, including many peaks due to seasonal influenza and heat waves during that period. LIMITATION: Historical vital statistics might be affected by minor completeness issues before the beginning of the 20th century. CONCLUSION: In 2020, the COVID-19 pandemic led to the second-largest infection-related mortality disaster in Switzerland, Sweden, and Spain since the beginning of the 20th century. PRIMARY FUNDING SOURCE: Foundation for Research in Science and the Humanities at the University of Zurich, Swiss National Science Foundation, and National Institute of Allergy and Infectious Diseases.

Estimating Tuberculosis Transmission Risks in a Primary Care Clinic in South Africa: Modeling of Environmental and Clinical Data.

BACKGROUND: Congregate settings, such as healthcare clinics, may play an essential role in Mycobacterium tuberculosis (Mtb) transmission. Using patient and environmental data, we studied transmission at a primary care clinic in South Africa. METHODS: We collected patient movements, cough frequency, and clinical data, and measured indoor carbon dioxide (CO2) levels, relative humidity, and Mtb genomes in the air. We used negative binomial regression model to investigate associations. RESULTS: We analyzed 978 unique patients who contributed 14 795 data points. The median patient age was 33 (interquartile range [IQR], 26-41) years, and 757 (77.4%) were female. Overall, median CO2 levels were 564 (IQR 495-646) parts per million and were highest in the morning. Median number of coughs per day was 466 (IQR, 368-503), and overall median Mtb DNA copies/µL/day was 4.2 (IQR, 1.2-9.5). We found an increased presence of Mtb DNA in the air of 32% (95% credible interval, 7%-63%) per 100 additional young adults (aged 15-29 years) and 1% (0-2%) more Mtb DNA per 10% increase of relative humidity. Estimated cumulative transmission risks for patients attending the clinic monthly for at least 1 hour range between 9% and 29%. CONCLUSIONS: We identified young adults and relative humidity as potentially important factors for transmission risks in healthcare clinics. Our approach should be used to detect transmission and evaluate infection control interventions.

Survival among people hospitalized with COVID-19 in Switzerland: a nationwide population-based analysis.

BACKGROUND: Increasing age, male sex, and pre-existing comorbidities are associated with lower survival from SARS-CoV-2 infection. The interplay between different comorbidities, age, and sex is not fully understood, and it remains unclear if survival decreases linearly with higher ICU occupancy or if there is a threshold beyond which survival falls. METHOD: This national population-based study included 22,648 people who tested positive for SARS-CoV-2 infection and were hospitalized in Switzerland between February 24, 2020, and March 01, 2021. Bayesian survival models were used to estimate survival after positive SARS-CoV-2 test among people hospitalized with COVID-19 by epidemic wave, age, sex, comorbidities, and ICU occupancy. Two-way interactions between age, sex, and comorbidities were included to assess the differential risk of death across strata. ICU occupancy was modeled using restricted cubic splines to allow for a non-linear association with survival. RESULTS: Of 22,648 people hospitalized with COVID-19, 4785 (21.1%) died. The survival was lower during the first epidemic wave than in the second (predicted survival at 40 days after positive test 76.1 versus 80.5%). During the second epidemic wave, occupancy among all available ICU beds in Switzerland varied between 51.7 and 78.8%. The estimated survival was stable at approximately 81.5% when ICU occupancy was below 70%, but worse when ICU occupancy exceeded this threshold (survival at 80% ICU occupancy: 78.2%; 95% credible interval [CrI] 76.1 to 80.1%). Periods with higher ICU occupancy (>70 vs 70%) were associated with an estimated number of 137 (95% CrI 27 to 242) excess deaths. Comorbid conditions reduced survival more in younger people than in older people. Among comorbid conditions, hypertension and obesity were not associated with poorer survival. Hypertension appeared to decrease survival in combination with cardiovascular disease. CONCLUSIONS: Survival after hospitalization with COVID-19 has improved over time, consistent with improved management of severe COVID-19. The decreased survival above 70% national ICU occupancy supports the need to introduce measures for prevention and control of SARS-CoV-2 transmission in the population well before ICUs are full.

Estimating the cumulative incidence of SARS-CoV-2 with imperfect serological tests: Exploiting cutoff-free approaches.

Large-scale serological testing in the population is essential to determine the true extent of the current SARS-CoV-2 pandemic. Serological tests measure antibody responses against pathogens and use predefined cutoff levels that dichotomize the quantitative test measures into sero-positives and negatives and use this as a proxy for past infection. With the imperfect assays that are currently available to test for past SARS-CoV-2 infection, the fraction of seropositive individuals in serosurveys is a biased estimator of the cumulative incidence and is usually corrected to account for the sensitivity and specificity. Here we use an inference method-referred to as mixture-model approach-for the estimation of the cumulative incidence that does not require to define cutoffs by integrating the quantitative test measures directly into the statistical inference procedure. We confirm that the mixture model outperforms the methods based on cutoffs, leading to less bias and error in estimates of the cumulative incidence. We illustrate how the mixture model can be used to optimize the design of serosurveys with imperfect serological tests. We also provide guidance on the number of control and case sera that are required to quantify the test's ambiguity sufficiently to enable the reliable estimation of the cumulative incidence. Lastly, we show how this approach can be used to estimate the cumulative incidence of classes of infections with an unknown distribution of quantitative test measures. This is a very promising application of the mixture-model approach that could identify the elusive fraction of asymptomatic SARS-CoV-2 infections. An R-package implementing the inference methods used in this paper is provided. Our study advocates using serological tests without cutoffs, especially if they are used to determine parameters characterizing populations rather than individuals. This approach circumvents some of the shortcomings of cutoff-based methods at exactly the low cumulative incidence levels and test accuracies that we are currently facing in SARS-CoV-2 serosurveys.

Bayesian workflow for disease transmission modeling in Stan.

This tutorial shows how to build, fit, and criticize disease transmission models in Stan, and should be useful to researchers interested in modeling the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic and other infectious diseases in a Bayesian framework. Bayesian modeling provides a principled way to quantify uncertainty and incorporate both data and prior knowledge into the model estimates. Stan is an expressive probabilistic programming language that abstracts the inference and allows users to focus on the modeling. As a result, Stan code is readable and easily extensible, which makes the modeler's work more transparent. Furthermore, Stan's main inference engine, Hamiltonian Monte Carlo sampling, is amiable to diagnostics, which means the user can verify whether the obtained inference is reliable. In this tutorial, we demonstrate how to formulate, fit, and diagnose a compartmental transmission model in Stan, first with a simple susceptible-infected-recovered model, then with a more elaborate transmission model used during the SARS-CoV-2 pandemic. We also cover advanced topics which can further help practitioners fit sophisticated models; notably, how to use simulations to probe the model and priors, and computational techniques to scale-up models based on ordinary differential equations.

Drivers of HIV-1 drug resistance to non-nucleoside reverse-transcriptase inhibitors (NNRTIs) in nine southern African countries: a modelling study.

INTRODUCTION: The rise of HIV-1 drug resistance to non-nucleoside reverse-transcriptase inhibitors (NNRTI) threatens antiretroviral therapy's long-term success (ART). NNRTIs will remain an essential drug for the management of HIV-1 due to safety concerns associated with integrase inhibitors. We fitted a dynamic transmission model to historical data from 2000 to 2018 in nine countries of southern Africa to understand the mechanisms that have shaped the HIV-1 epidemic and the rise of pretreatment NNRTI resistance. METHODS: We included data on HIV-1 prevalence, ART coverage, HIV-related mortality, and survey data on pretreatment NNRTI resistance from nine southern Africa countries from a systematic review, UNAIDS and World Bank. Using a Bayesian hierarchical framework, we developed a dynamic transmission model linking data on the HIV-1 epidemic to survey data on NNRTI drug resistance in each country. We estimated the proportion of resistance attributable to unregulated, off-programme use of ART. We examined each national ART programme's vulnerability to NNRTI resistance by defining a fragility index: the ratio of the rate of NNRTI resistance emergence during first-line ART over the rate of switching to second-line ART. We explored associations between fragility and characteristics of the health system of each country. RESULTS: The model reliably described the dynamics of the HIV-1 epidemic and NNRTI resistance in each country. Predicted levels of resistance in 2018 ranged between 3.3% (95% credible interval 1.9-7.1) in Mozambique and 25.3% (17.9-33.8) in Eswatini. The proportion of pretreatment NNRTI resistance attributable to unregulated antiretroviral use ranged from 6% (2-14) in Eswatini to 64% (26-85) in Mozambique. The fragility index was low in Botswana (0.01; 0.0-0.11) but high in Namibia (0.48; 0.16-10.17), Eswatini (0.64; 0.23-11.8) and South Africa (1.21; 0.83-9.84). The combination of high fragility of ART programmes and high ART coverage levels was associated with a sharp increase in pretreatment NNRTI resistance. CONCLUSIONS: This comparison of nine countries shows that pretreatment NNRTI resistance can be controlled despite high ART coverage levels. This was the case in Botswana, Mozambique, and Zambia, most likely because of better HIV care delivery, including rapid switching to second-line ART of patients failing first-line ART.

Estimation of SARS-CoV-2 mortality during the early stages of an epidemic: A modeling study in Hubei, China, and six regions in Europe.

BACKGROUND: As of 16 May 2020, more than 4.5 million cases and more than 300,000 deaths from disease caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) have been reported. Reliable estimates of mortality from SARS-CoV-2 infection are essential for understanding clinical prognosis, planning healthcare capacity, and epidemic forecasting. The case-fatality ratio (CFR), calculated from total numbers of reported cases and reported deaths, is the most commonly reported metric, but it can be a misleading measure of overall mortality. The objectives of this study were to (1) simulate the transmission dynamics of SARS-CoV-2 using publicly available surveillance data and (2) infer estimates of SARS-CoV-2 mortality adjusted for biases and examine the CFR, the symptomatic case-fatality ratio (sCFR), and the infection-fatality ratio (IFR) in different geographic locations. METHOD AND FINDINGS: We developed an age-stratified susceptible-exposed-infected-removed (SEIR) compartmental model describing the dynamics of transmission and mortality during the SARS-CoV-2 epidemic. Our model accounts for two biases: preferential ascertainment of severe cases and right-censoring of mortality. We fitted the transmission model to surveillance data from Hubei Province, China, and applied the same model to six regions in Europe: Austria, Bavaria (Germany), Baden-Württemberg (Germany), Lombardy (Italy), Spain, and Switzerland. In Hubei, the baseline estimates were as follows: CFR 2.4% (95% credible interval [CrI] 2.1%-2.8%), sCFR 3.7% (3.2%-4.2%), and IFR 2.9% (2.4%-3.5%). Estimated measures of mortality changed over time. Across the six locations in Europe, estimates of CFR varied widely. Estimates of sCFR and IFR, adjusted for bias, were more similar to each other but still showed some degree of heterogeneity. Estimates of IFR ranged from 0.5% (95% CrI 0.4%-0.6%) in Switzerland to 1.4% (1.1%-1.6%) in Lombardy, Italy. In all locations, mortality increased with age. Among individuals 80 years or older, estimates of the IFR suggest that the proportion of all those infected with SARS-CoV-2 who will die ranges from 20% (95% CrI 16%-26%) in Switzerland to 34% (95% CrI 28%-40%) in Spain. A limitation of the model is that count data by date of onset are required, and these are not available in all countries. CONCLUSIONS: We propose a comprehensive solution to the estimation of SARS-Cov-2 mortality from surveillance data during outbreaks. The CFR is not a good predictor of overall mortality from SARS-CoV-2 and should not be used for evaluation of policy or comparison across settings. Geographic differences in IFR suggest that a single IFR should not be applied to all settings to estimate the total size of the SARS-CoV-2 epidemic in different countries. The sCFR and IFR, adjusted for right-censoring and preferential ascertainment of severe cases, are measures that can be used to improve and monitor clinical and public health strategies to reduce the deaths from SARS-CoV-2 infection.

Impact of scaling up dolutegravir on antiretroviral resistance in South Africa: A modeling study.

BACKGROUND: Rising resistance of HIV-1 to non-nucleoside reverse transcriptase inhibitors (NNRTIs) threatens the success of the global scale-up of antiretroviral therapy (ART). The switch to WHO-recommended dolutegravir (DTG)-based regimens could reduce this threat due to DTG's high genetic barrier to resistance. We used mathematical modeling to predict the impact of the scale-up of DTG-based ART on NNRTI pretreatment drug resistance (PDR) in South Africa, 2020 to 2040. METHODS AND FINDINGS: We adapted the Modeling Antiretroviral drug Resistance In South Africa (MARISA) model, an epidemiological model of the transmission of NNRTI resistance in South Africa. We modeled the introduction of DTG in 2020 under 2 scenarios: DTG as first-line regimen for ART initiators, or DTG for all patients, including patients on suppressive NNRTI-based ART. Given the safety concerns related to DTG during pregnancy, we assessed the impact of prescribing DTG to all men and in addition to (1) women beyond reproductive age; (2) women beyond reproductive age or using contraception; and (3) all women. The model projections show that, compared to the continuation of NNRTI-based ART, introducing DTG would lead to a reduction in NNRTI PDR in all scenarios if ART initiators are started on a DTG-based regimen, and those on NNRTI-based regimens are rapidly switched to DTG. NNRTI PDR would continue to increase if DTG-based ART was restricted to men. When given to all men and women, DTG-based ART could reduce the level of NNRTI PDR from 52.4% (without DTG) to 10.4% (with universal DTG) in 2040. If only men and women beyond reproductive age or on contraception are started on or switched to DTG-based ART, NNRTI PDR would reach 25.9% in 2040. Limitations include substantial uncertainty due to the long-term predictions and the current scarcity of knowledge about DTG efficacy in South Africa. CONCLUSIONS: Our model shows the potential benefit of scaling up DTG-based regimens for halting the rise of NNRTI resistance. Starting or switching all men and women to DTG would lead to a sustained decline in resistance levels, whereas using DTG-based ART in all men, or in men and women beyond childbearing age, would only slow down the increase in levels of NNRTI PDR.

Bridging the gap between HIV epidemiology and antiretroviral resistance evolution: Modelling the spread of resistance in South Africa.

The scale-up of antiretroviral therapy (ART) in South Africa substantially reduced AIDS-related deaths and new HIV infections. However, its success is threatened by the emergence of resistance to non-nucleoside reverse-transcriptase inhibitors (NNRTI). The MARISA (Modelling Antiretroviral drug Resistance In South Africa) model presented here aims at investigating the time trends and factors driving NNRTI resistance in South Africa. MARISA is a compartmental model that includes the key aspects of the local HIV epidemic: continuum of care, disease progression, and gender. The dynamics of NNRTI resistance emergence and transmission are then added to this framework. Model parameters are informed using data from HIV cohorts participating in the International epidemiology Databases to Evaluate AIDS (IeDEA) and literature estimates, or fitted to UNAIDS estimates. Using this novel approach of triangulating clinical and resistance data from various sources, MARISA reproduces the time trends of HIV in South Africa in 2005-2016, with a decrease in new infections, undiagnosed individuals, and AIDS-related deaths. MARISA captures the dynamics of the spread of NNRTI resistance: high levels of acquired drug resistance (ADR, in 83% of first-line treatment failures in 2016), and increasing transmitted drug resistance (TDR, in 8.1% of ART initiators in 2016). Simulation of counter-factual scenarios reflecting alternative public health policies shows that increasing treatment coverage would have resulted in fewer new infections and deaths, at the cost of higher TDR (11.6% in 2016 for doubling the treatment rate). Conversely, improving switching to second-line treatment would have led to lower TDR (6.5% in 2016 for doubling the switching rate) and fewer new infections and deaths. Implementing drug resistance testing would have had little impact. The rapid ART scale-up and inadequate switching to second-line treatment were the key drivers of the spread of NNRTI resistance in South Africa. However, even though some interventions could have substantially reduced the level of NNRTI resistance, no policy including NNRTI-based first line regimens could have prevented this spread. Thus, by combining epidemiological data on HIV in South Africa with biological data on resistance evolution, our modelling approach identified key factors driving NNRTI resistance, highlighting the need of alternative first-line regimens.

Pattern of early human-to-human transmission of Wuhan 2019 novel coronavirus (2019-nCoV), December 2019 to January 2020.

Since December 2019, China has been experiencing a large outbreak of a novel coronavirus (2019-nCoV) which can cause respiratory disease and severe pneumonia. We estimated the basic reproduction number R0 of 2019-nCoV to be around 2.2 (90% high density interval: 1.4-3.8), indicating the potential for sustained human-to-human transmission. Transmission characteristics appear to be of similar magnitude to severe acute respiratory syndrome-related coronavirus (SARS-CoV) and pandemic influenza, indicating a risk of global spread.

Interhuman transmissibility of Middle East respiratory syndrome coronavirus: estimation of pandemic risk.

BACKGROUND: The new Middle East respiratory syndrome coronavirus (MERS-CoV) infection shares many clinical, epidemiological, and virological similarities with that of severe acute respiratory syndrome (SARS)-CoV. We aimed to estimate virus transmissibility and the epidemic potential of MERS-CoV, and to compare the results with similar findings obtained for prepandemic SARS. METHODS: We retrieved data for MERS-CoV clusters from the WHO summary and subsequent reports, and published descriptions of cases, and took into account 55 of the 64 laboratory-confirmed cases of MERS-CoV reported as of June 21, 2013, excluding cases notified in the previous 2 weeks. To assess the interhuman transmissibility of MERS-CoV, we used Bayesian analysis to estimate the basic reproduction number (R0) and compared it to that of prepandemic SARS. We considered two scenarios, depending on the interpretation of the MERS-CoV cluster-size data. RESULTS: With our most pessimistic scenario (scenario 2), we estimated MERS-CoV R0 to be 0·69 (95% CI 0·50-0·92); by contrast, the R0 for prepandemic SARS-CoV was 0·80 (0·54-1·13). Our optimistic scenario (scenario 1) yielded a MERS-CoV R0 of 0·60 (0·42-0·80). Because of recent implementation of effective contact tracing and isolation procedures, further MERS-CoV transmission data might no longer describe an entire cluster, but only secondary infections directly caused by the index patient. Hence, we calculated that, under scenario 2, eight or more secondary infections caused by the next index patient would translate into a 5% or higher chance that the revised MERS-CoV R0 would exceed 1--ie, that MERS-CoV might have pandemic potential. INTERPRETATION: Our analysis suggests that MERS-CoV does not yet have pandemic potential. We recommend enhanced surveillance, active contact tracing, and vigorous searches for the MERS-CoV animal hosts and transmission routes to human beings. FUNDING: Agence Nationale de la Recherche (Labex Integrative Biology of Emerging Infectious Diseases), and the European Community's Seventh Framework Programme project PREDEMICS.

Gender/Sex Disparities in the COVID-19 Cascade From Testing to Mortality: An Intersectional Analysis of Swiss Surveillance Data.

Objectives: This study investigates gender and sex disparities in COVID-19 epidemiology in the Canton of Vaud, Switzerland, focusing on the interplay with socioeconomic position (SEP) and age. Methods: We analyzed COVID-19 surveillance data from March 2020 to June 2021, using an intersectional approach. Negative binomial regression models assessed disparities between women and men, across SEP quintiles and age groups, in testing, positivity, hospitalizations, ICU admissions, and mortality (Incidence Rate Ratios [IRR], with 95% Confidence Intervals [CI]). Results: Women had higher testing and positivity rates than men, while men experienced more hospitalizations, ICU admissions, and deaths. The higher positivity in women under 50 was mitigated when accounting for their higher testing rates. Within SEP quintiles, gender/sex differences in testing and positivity were not significant. In the lowest quintile, women's mortality risk was 68% lower (Q1: IRR 0.32, CI 0.20-0.52), with decreasing disparities with increasing SEP quintiles (Q5: IRR 0.66, CI 0.41-1.06). Conclusion: Our findings underscore the complex epidemiological patterns of COVID-19, shaped by the interactions of gender/sex, SEP, and age, highlighting the need for intersectional perspectives in both epidemiological research and public health strategy development.

Direct and indirect effects of the COVID-19 pandemic on mortality in Switzerland.

The direct and indirect impact of the COVID-19 pandemic on population-level mortality is of concern to public health but challenging to quantify. Using data for 2011-2019, we applied Bayesian models to predict the expected number of deaths in Switzerland and compared them with laboratory-confirmed COVID-19 deaths from February 2020 to April 2022 (study period). We estimated that COVID-19-related mortality was underestimated by a factor of 0.72 (95% credible interval [CrI]: 0.46-0.78). After accounting for COVID-19 deaths, the observed mortality was -4% (95% CrI: -8 to 0) lower than expected. The deficit in mortality was concentrated in age groups 40-59 (-12%, 95%CrI: -19 to -5) and 60-69 (-8%, 95%CrI: -15 to -2). Although COVID-19 control measures may have negative effects, after subtracting COVID-19 deaths, there were fewer deaths in Switzerland during the pandemic than expected, suggesting that any negative effects of control measures were offset by the positive effects. These results have important implications for the ongoing debate about the appropriateness of COVID-19 control measures.