Spatial spread of COVID-19 during the early pandemic phase in Italy.

Quantifying the potential spatial spread of an infectious pathogen is key to defining effective containment and control strategies. The aim of this study is to estimate the risk of SARS-CoV-2 transmission at different distances in Italy before the first regional lockdown was imposed, identifying important sources of national spreading. To do this, we leverage on a probabilistic model applied to daily symptomatic cases retrospectively ascertained in each Italian municipality with symptom onset between January 28 and March 7, 2020. Results are validated using a multi-patch dynamic transmission model reproducing the spatiotemporal distribution of identified cases. Our results show that the contribution of short-distance ( ≤ 10 k m ) transmission increased from less than 40% in the last week of January to more than 80% in the first week of March 2020. On March 7, 2020, that is the day before the first regional lockdown was imposed, more than 200 local transmission foci were contributing to the spread of SARS-CoV-2 in Italy. At the time, isolation measures imposed only on municipalities with at least ten ascertained cases would have left uncontrolled more than 75% of spillover transmission from the already affected municipalities. In early March, national-wide restrictions were required to curb short-distance transmission of SARS-CoV-2 in Italy.

Retrospective analysis of the Italian exit strategy from COVID-19 lockdown.

After the national lockdown imposed on March 11, 2020, the Italian government has gradually resumed the suspended economic and social activities since May 4, while maintaining the closure of schools until September 14. We use a model of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) transmission to estimate the health impact of different exit strategies. The strategy adopted in Italy kept the reproduction number Rt at values close to one until the end of September, with marginal regional differences. Based on the estimated postlockdown transmissibility, reopening of workplaces in selected industrial activities might have had a minor impact on the transmissibility. Reopening educational levels in May up to secondary schools might have influenced SARS-CoV-2 transmissibility only marginally; however, including high schools might have resulted in a marked increase of the disease burden. Earlier reopening would have resulted in disproportionately higher hospitalization incidence. Given community contacts in September, we project a large second wave associated with school reopening in the fall.

Estimates of Serial Interval and Reproduction Number of Sudan Virus, Uganda, August-November 2022.

We estimated the mean serial interval for Sudan virus in Uganda to be 11.7 days (95 CI% 8.2-15.8 days). Estimates for the 2022 outbreak indicate a mean basic reproduction number of 2.4-2.7 (95% CI 1.7-3.5). Estimated net reproduction numbers across districts suggest a marked spatial heterogeneity.

Risk of Symptomatic Infection During a Second Coronavirus Disease 2019 Wave in Severe Acute Respiratory Syndrome Coronavirus 2-Seropositive Individuals.

We analyzed 221 coronavirus disease 2019 cases identified between June 2020 and January 2021 in 6074 individuals screened for immunoglobulin G antibodies in May 2020, representing 77% of residents of 5 Italian municipalities. The relative risk of developing symptomatic infection in seropositive participants was 0.055 (95% confidence interval, .014-.220).

Early Estimates of Monkeypox Incubation Period, Generation Time, and Reproduction Number, Italy, May-June 2022.

We analyzed the first 255 PCR-confirmed cases of monkeypox in Italy in 2022. Preliminary estimates indicate mean incubation period of 9.1 (95% CI 6.5-10.9) days, mean generation time of 12.5 (95% CI 7.5-17.3) days, and reproduction number among men who have sex with men of 2.43 (95% CI 1.82-3.26).

Pressure on the Health-Care System and Intensive Care Utilization During the COVID-19 Outbreak in the Lombardy Region of Italy: A Retrospective Observational Study in 43,538 Hospitalized Patients.

During the spring of 2020, the coronavirus disease 2019 (COVID-19) epidemic caused an unprecedented demand for intensive-care resources in the Lombardy region of Italy. Using data on 43,538 hospitalized patients admitted between February 21 and July 12, 2020, we evaluated variations in intensive care unit (ICU) admissions and mortality over the course of 3 periods: the early phase of the pandemic (February 21-March 13), the period of highest pressure on the health-care system (March 14-April 25, when numbers of COVID-19 patients exceeded prepandemic ICU bed capacity), and the declining phase (April 26-July 12). Compared with the early phase, patients aged 70 years or more were less often admitted to an ICU during the period of highest pressure on the health-care system (odds ratio (OR) = 0.47, 95% confidence interval (CI): 0.41, 0.54), with longer ICU delays (incidence rate ratio = 1.82, 95% CI: 1.52, 2.18) and lower chances of dying in the ICU (OR = 0.47, 95% CI: 0.34, 0.64). Patients under 56 years of age had more limited changes in the probability of (OR = 0.65, 95% CI: 0.56, 0.76) and delay to (incidence rate ratio = 1.16, 95% CI: 0.95, 1.42) ICU admission and increased mortality (OR = 1.43, 95% CI: 1.00, 2.07). In the declining phase, all quantities decreased for all age groups. These patterns may suggest that limited health-care resources during the peak phase of the epidemic in Lombardy forced a shift in ICU admission criteria to prioritize patients with higher chances of survival.

Estimation of the incubation period and generation time of SARS-CoV-2 Alpha and Delta variants from contact tracing data.

Quantitative information on epidemiological quantities such as the incubation period and generation time of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants is scarce. We analysed a dataset collected during contact tracing activities in the province of Reggio Emilia, Italy, throughout 2021. We determined the distributions of the incubation period for the Alpha and Delta variants using information on negative polymerase chain reaction tests and the date of last exposure from 282 symptomatic cases. We estimated the distributions of the intrinsic generation time using a Bayesian inference approach applied to 9724 SARS-CoV-2 cases clustered in 3545 households where at least one secondary case was recorded. We estimated a mean incubation period of 4.9 days (95% credible intervals, CrI, 4.4-5.4) for Alpha and 4.5 days (95% CrI 4.0-5.0) for Delta. The intrinsic generation time was estimated to have a mean of 7.12 days (95% CrI 6.27-8.44) for Alpha and of 6.52 days (95% CrI 5.54-8.43) for Delta. The household serial interval was 2.43 days (95% CrI 2.29-2.58) for Alpha and 2.74 days (95% CrI 2.62-2.88) for Delta, and the estimated proportion of pre-symptomatic transmission was 48-51% for both variants. These results indicate limited differences in the incubation period and intrinsic generation time of SARS-CoV-2 variants Alpha and Delta compared to ancestral lineages.

Co-circulation of SARS-CoV-2 Alpha and Gamma variants in Italy, February and March 2021.

BackgroundSeveral SARS-CoV-2 variants of concern (VOC) have emerged through 2020 and 2021. There is need for tools to estimate the relative transmissibility of emerging variants of SARS-CoV-2 with respect to circulating strains.AimWe aimed to assess the prevalence of co-circulating VOC in Italy and estimate their relative transmissibility.MethodsWe conducted two genomic surveillance surveys on 18 February and 18 March 2021 across the whole Italian territory covering 3,243 clinical samples and developed a mathematical model that describes the dynamics of co-circulating strains.ResultsThe Alpha variant was already dominant on 18 February in a majority of regions/autonomous provinces (national prevalence: 54%) and almost completely replaced historical lineages by 18 March (dominant across Italy, national prevalence: 86%). We found a substantial proportion of the Gamma variant on 18 February, almost exclusively in central Italy (prevalence: 19%), which remained similar on 18 March. Nationally, the mean relative transmissibility of Alpha ranged at 1.55-1.57 times the level of historical lineages (95% CrI: 1.45-1.66). The relative transmissibility of Gamma varied according to the assumed degree of cross-protection from infection with other lineages and ranged from 1.12 (95% CrI: 1.03-1.23) with complete immune evasion to 1.39 (95% CrI: 1.26-1.56) for complete cross-protection.ConclusionWe assessed the relative advantage of competing viral strains, using a mathematical model assuming different degrees of cross-protection. We found substantial co-circulation of Alpha and Gamma in Italy. Gamma was not able to outcompete Alpha, probably because of its lower transmissibility.

Tracking the progressive spread of the SARS-CoV-2 Omicron variant in Italy, December 2021 to January 2022.

BackgroundThe SARS-CoV-2 variant of concern Omicron was first detected in Italy in November 2021.AimTo comprehensively describe Omicron spread in Italy in the 2 subsequent months and its impact on the overall SARS-CoV-2 circulation at population level.MethodsWe analyse data from four genomic surveys conducted across the country between December 2021 and January 2022. Combining genomic sequencing results with epidemiological records collated by the National Integrated Surveillance System, the Omicron reproductive number and exponential growth rate are estimated, as well as SARS-CoV-2 transmissibility.ResultsOmicron became dominant in Italy less than 1 month after its first detection, representing on 3 January 76.9-80.2% of notified SARS-CoV-2 infections, with a doubling time of 2.7-3.3 days. As of 17 January 2022, Delta variant represented < 6% of cases. During the Omicron expansion in December 2021, the estimated mean net reproduction numbers respectively rose from 1.15 to a maximum of 1.83 for symptomatic cases and from 1.14 to 1.36 for hospitalised cases, while remaining relatively stable, between 0.93 and 1.21, for cases needing intensive care. Despite a reduction in relative proportion, Delta infections increased in absolute terms throughout December contributing to an increase in hospitalisations. A significant reproduction numbers' decline was found after mid-January, with average estimates dropping below 1 between 10 and 16 January 2022.ConclusionEstimates suggest a marked growth advantage of Omicron compared with Delta variant, but lower disease severity at population level possibly due to residual immunity against severe outcomes acquired from vaccination and prior infection.

Impact of a Nationwide Lockdown on SARS-CoV-2 Transmissibility, Italy.

On March 11, 2020, Italy imposed a national lockdown to curtail the spread of severe acute respiratory syndrome coronavirus 2. We estimate that, 14 days after lockdown, the net reproduction number had dropped below 1 and remained stable at ¼0.76 (95% CI 0.67-0.85) in all regions for >3 of the following weeks.

Modeling the interplay between demography, social contact patterns, and SARS-CoV-2 transmission in the South West Shewa Zone of Oromia Region, Ethiopia.

BACKGROUND: COVID-19 spread may have a dramatic impact in countries with vulnerable economies and limited availability of, and access to, healthcare resources and infrastructures. However, in sub-Saharan Africa, a low prevalence and mortality have been observed so far. METHODS: We collected data on individuals' social contacts in the South West Shewa Zone (SWSZ) of Ethiopia across geographical contexts characterized by heterogeneous population density, work and travel opportunities, and access to primary care. We assessed how socio-demographic factors and observed mixing patterns can influence the COVID-19 disease burden, by simulating SARS-CoV-2 transmission in remote settlements, rural villages, and urban neighborhoods, under school closure mandate. RESULTS: From national surveillance data, we estimated a net reproduction number of 1.62 (95% CI 1.55-1.70). We found that, at the end of an epidemic mitigated by school closure alone, 10-15% of the population residing in the SWSZ would have been symptomatic and 0.3-0.4% of the population would require mechanical ventilation and/or possibly result in a fatal outcome. Higher infection attack rates are expected in more urbanized areas, but the highest incidence of critical disease is expected in remote subsistence farming settlements. School closure contributed to reduce the reproduction number by 49% and the attack rate of infections by 28-34%. CONCLUSIONS: Our results suggest that the relatively low burden of COVID-19 in Ethiopia observed so far may depend on social mixing patterns, underlying demography, and the enacted school closures. Our findings highlight that socio-demographic factors can also determine marked heterogeneities across different geographical contexts within the same region, and they contribute to understand why sub-Saharan Africa is experiencing a relatively lower attack rate of severe cases compared to high-income countries.

Potential short-term outcome of an uncontrolled COVID-19 epidemic in Lombardy, Italy, February to March 2020.

Sustained coronavirus disease (COVID-19) transmission is ongoing in Italy, with 7,375 reported cases and 366 deaths by 8 March 2020. We provide a model-based evaluation of patient records from Lombardy, predicting the impact of an uncontrolled epidemic on the healthcare system. It has the potential to cause more than 250,039 (95% credible interval (CrI): 147,717-459,890) cases within 3 weeks, including 37,194 (95% CrI: 22,250-67,632) patients requiring intensive care. Aggressive containment strategies are required.

Age-specific SARS-CoV-2 infection fatality ratio and associated risk factors, Italy, February to April 2020.

We analysed 5,484 close contacts of coronavirus disease (COVID-19) cases in Italy, all tested for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Infection fatality ratio was 0.43% (95% confidence interval (CI): 0.21-0.79) for individuals younger than 70 years and 10.5% (95% CI: 8.0-13.6) for older individuals. Risk of death after infection was 62% lower (95% CI: 31-80) in clusters identified after 16 March 2020 and 1.8-fold higher for males (95% CI: 1.03-3.16).

Epidemiological characteristics of COVID-19 cases and estimates of the reproductive numbers 1 month into the epidemic, Italy, 28 January to 31 March 2020.

BackgroundOn 20 February 2020, a locally acquired coronavirus disease (COVID-19) case was detected in Lombardy, Italy. This was the first signal of ongoing transmission of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in the country. The number of cases in Italy increased rapidly and the country became the first in Europe to experience a SARS-CoV-2 outbreak.AimOur aim was to describe the epidemiology and transmission dynamics of the first COVID-19 cases in Italy amid ongoing control measures.MethodsWe analysed all RT-PCR-confirmed COVID-19 cases reported to the national integrated surveillance system until 31 March 2020. We provide a descriptive epidemiological summary and estimate the basic and net reproductive numbers by region.ResultsOf the 98,716 cases of COVID-19 analysed, 9,512 were healthcare workers. Of the 10,943 reported COVID-19-associated deaths (crude case fatality ratio: 11.1%) 49.5% occurred in cases older than 80 years. Male sex and age were independent risk factors for COVID-19 death. Estimates of R0 varied between 2.50 (95% confidence interval (CI): 2.18-2.83) in Tuscany and 3.00 (95% CI: 2.68-3.33) in Lazio. The net reproduction number Rt in northern regions started decreasing immediately after the first detection.ConclusionThe COVID-19 outbreak in Italy showed a clustering onset similar to the one in Wuhan, China. R0 at 2.96 in Lombardy combined with delayed detection explains the high case load and rapid geographical spread. Overall, Rt in Italian regions showed early signs of decrease, with large diversity in incidence, supporting the importance of combined non-pharmacological control measures.

Modeling the impact of changes in day-care contact patterns on the dynamics of varicella transmission in France between 1991 and 2015.

Annual incidence rates of varicella infection in the general population in France have been rather stable since 1991 when clinical surveillance started. Rates however show a statistically significant increase over time in children aged 0-3 years, and a decline in older individuals. A significant increase in day-care enrolment and structures' capacity in France was also observed in the last decade. In this work we investigate the potential interplay between an increase of contacts of young children possibly caused by earlier socialization in the community and varicella transmission dynamics. To this aim, we develop an age-structured mathematical model, informed with historical demographic data and contact matrix estimates in the country, accounting for longitudinal linear increase of early childhood contacts. While the reported overall varicella incidence is well reproduced independently of mixing variations, age-specific empirical trends are better captured by accounting for an increase in contacts among pre-school children in the last decades. We found that the varicella data are consistent with a 30% increase in the number of contacts at day-care facilities, which would imply a 50% growth in the contribution of 0-3y old children to overall yearly infections in 1991-2015. Our findings suggest that an earlier exposure to pathogens due to changes in day-care contact patterns, represents a plausible explanation for the epidemiological patterns observed in France. Obtained results suggest that considering temporal changes in social factors in addition to demographic ones is critical to correctly interpret varicella transmission dynamics.

The impact of demographic changes, exogenous boosting and new vaccination policies on varicella and herpes zoster in Italy: a modelling and cost-effectiveness study.

BACKGROUND: The present study aims to evaluate the cost-effectiveness of the newly introduced varicella and herpes zoster (HZ) vaccination programmes in Italy. The appropriateness of the introduction of the varicella vaccine is highly debated because of concerns about the consequences on HZ epidemiology and the expected increase in the number of severe cases in case of suboptimal coverage levels. METHODS: We performed a cost-utility analysis based on a stochastic individual-based model that considers realistic demographic processes and two different underlying mechanisms of exogenous boosting (temporary and progressive immunity). Routine varicella vaccination is given with a two-dose schedule (15 months, 5-6 years). The HZ vaccine is offered to the elderly (65 years), either alone or in combination with an initial catch-up campaign (66-75 years). The main outcome measures are averted cases and deaths, costs per quality-adjusted life years gained, incremental cost-effectiveness ratios, and net monetary benefits associated with the different vaccination policies. RESULTS: Demographic processes have contributed to shaping varicella and HZ epidemiology over the years, decreasing varicella circulation and increasing the incidence of HZ. The recent introduction of varicella vaccination in Italy is expected to produce an enduring reduction in varicella incidence and, indirectly, a further increase of HZ incidence in the first decades, followed by a significant reduction in the long term. However, the concurrent introduction of routine HZ vaccination at 65 years of age is expected to mitigate this increase and, in the longer run, to reduce HZ burden to its minimum. From an economic perspective, all the considered policies are cost-effective, with the exception of varicella vaccination alone when considering a time horizon of 50 years. These results are robust to parameter uncertainties, to the two different hypotheses on the mechanism driving exogenous boosting, and to different demographic projection scenarios. CONCLUSIONS: The recent introduction of a combined varicella and HZ vaccination programme in Italy will produce significant reductions in the burden of both diseases and is found to be a cost-effective policy. This programme will counterbalance the increasing trend of zoster incidence purely due to demographic processes.

The impact of demographic changes on the epidemiology of herpes zoster: Spain as a case study.

Varicella zoster virus (VZV) causes varicella upon first exposure and may reactivate later in life into herpes zoster (HZ), with a risk that is thought to be reduced by re-exposures to VZV. Given the decades-long time scales of reactivation and its dependence on the accumulation of re-exposure episodes, adopting a long-term perspective may be useful to correctly interpret current epidemiological trends of VZV. In this study, we investigate the possible impact of demographic changes on varicella and HZ in Spain, using an age-structured mathematical model informed with historical demographic data and calibrated against age-specific profiles of varicella seroprevalence and HZ incidence data. The model qualitatively reproduces the remarkable growth of HZ incidence observed in Spain between 1997 and 2004, before the introduction of varicella vaccination programmes. We demonstrate that this growth may be partially ascribed to the reduction of varicella circulation that followed the overall decline of the birth rate in the twentieth century. Model predictions further suggest that, even under the most optimistic projections, HZ incidence will continue its rise until at least 2040. Considering the effect of demographic changes can help interpreting variations in epidemiological trends of HZ, contributing to a more accurate evaluation of vaccination programmes against VZV.

The impact of underreported infections on vaccine effectiveness estimates derived from retrospective cohort studies.

BACKGROUND: Surveillance data and vaccination registries are widely used to provide real-time vaccine effectiveness (VE) estimates, which can be biased due to underreported (i.e. under-ascertained and under-notified) infections. Here, we investigate how the magnitude and direction of this source of bias in retrospective cohort studies vary under different circumstances, including different levels of underreporting, heterogeneities in underreporting across vaccinated and unvaccinated, and different levels of pathogen circulation. METHODS: We developed a stochastic individual-based model simulating the transmission dynamics of a respiratory virus and a large-scale vaccination campaign. Considering a baseline scenario with 22.5% yearly attack rate and 30% reporting ratio, we explored fourteen alternative scenarios, each modifying one or more baseline assumptions. Using synthetic individual-level surveillance data and vaccination registries produced by the model, we estimated the VE against documented infection taking as reference either unvaccinated or recently vaccinated individuals (within 14 days post-administration). Bias was quantified by comparing estimates to the known VE assumed in the model. RESULTS: VE estimates were accurate when assuming homogeneous reporting ratios, even at low levels (10%), and moderate attack rates (<50%). A substantial downward bias in the estimation arose with homogeneous reporting and attack rates exceeding 50%. Mild heterogeneities in reporting ratios between vaccinated and unvaccinated strongly biased VE estimates, downward if cases in vaccinated were more likely to be reported and upward otherwise, particularly when taking as reference unvaccinated individuals. CONCLUSIONS: In observational studies, high attack rates or differences in underreporting between vaccinated and unvaccinated may result in biased VE estimates. This study underscores the critical importance of monitoring data quality and understanding biases in observational studies, to more adequately inform public health decisions.

The decline of the 2022 Italian mpox epidemic: Role of behavior changes and control strategies.

In 2022, a global outbreak of mpox occurred, predominantly impacting men who have sex with men (MSM). The rapid decline of this epidemic is yet to be fully understood. We investigated the Italian outbreak by means of an individual-based mathematical model calibrated to surveillance data. The model accounts for transmission within the MSM sexual contact network, in recreational and sex clubs attended by MSM, and in households. We indicate a strong spontaneous reduction in sexual transmission (61-87%) in affected MSM communities as the possible driving factor for the rapid decline in cases. The MSM sexual contact network was the main responsible for transmission (about 80%), with clubs and households contributing residually. Contact tracing prevented about half of the potential cases, and a higher success rate in tracing contacts could significantly amplify its effectiveness. Notably, immunizing the 23% of MSM with the highest sexual activity (10 or more partners per year) could completely prevent new mpox resurgences. This research underscores the importance of augmenting contact tracing, targeted immunization campaigns of high-risk groups, and fostering reactive behavioral changes as key strategies to manage and prevent the spread of emerging sexually transmitted pathogens like mpox within the MSM community.

Estimating the potential risk of transmission of arboviruses in the Americas and Europe: a modelling study.

BACKGROUND: Estimates of the spatiotemporal distribution of different mosquito vector species and the associated risk of transmission of arboviruses are key to design adequate policies for preventing local outbreaks and reducing the number of human infections in endemic areas. In this study, we quantified the abundance of Aedes albopictus and Aedes aegypti and the local transmission potential for three arboviral infections at an unprecedented spatiotemporal resolution in areas where no entomological surveillance is available. METHODS: We developed a computational model to quantify the daily abundance of Aedes mosquitoes, leveraging temperature and precipitation records. The model was calibrated on mosquito surveillance data collected in 115 locations in Europe and the Americas between 2007 and 2018. Model estimates were used to quantify the reproduction number of dengue virus, Zika virus, and chikungunya in Europe and the Americas, at a high spatial resolution. FINDINGS: In areas colonised by both Aedes species, A aegypti was estimated to be the main vector for the transmission of dengue virus, Zika virus, and chikungunya, being associated with a higher estimate of R0 when compared with A albopictus. Our estimates highlighted that these arboviruses were endemic in tropical and subtropical countries, with the highest risks of transmission found in central America, Venezuela, Colombia, and central-east Brazil. A non-negligible potential risk of transmission was also estimated for Florida, Texas, and Arizona (USA). The broader ecological niche of A albopictus could contribute to the emergence of chikungunya outbreaks and clusters of dengue autochthonous cases in temperate areas of the Americas, as well as in mediterranean Europe (in particular, in Italy, southern France, and Spain). INTERPRETATION: Our results provide a comprehensive overview of the transmission potential of arboviral diseases in Europe and the Americas, highlighting areas where surveillance and mosquito control capacities should be prioritised. FUNDING: EU and Ministero dell'Università e della Ricerca, Italy (Piano Nazionale di Ripresa e Resilienza Extended Partnership initiative on Emerging Infectious Diseases); EU (Horizon 2020); Ministero dell'Università e della Ricerca, Italy (Progetti di ricerca di Rilevante Interesse Nazionale programme); Brazilian National Council of Science, Technology and Innovation; Ministry of Health, Brazil; and Foundation of Research for Minas Gerais, Brazil.