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.

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).

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.

Spatial modes for transmission of chikungunya virus during a large chikungunya outbreak in Italy: a modeling analysis.

BACKGROUND: The spatial spread of many mosquito-borne diseases occurs by focal spread at the scale of a few hundred meters and over longer distances due to human mobility. The relative contributions of different spatial scales for transmission of chikungunya virus require definition to improve outbreak vector control recommendations. METHODS: We analyzed data from a large chikungunya outbreak mediated by the mosquito Aedes albopictus in the Lazio region, Italy, consisting of 414 reported human cases between June and November 2017. Using dates of symptom onset, geographic coordinates of residence, and information from epidemiological questionnaires, we reconstructed transmission chains related to that outbreak. RESULTS: Focal spread (within 1 km) accounted for 54.9% of all cases, 15.8% were transmitted at a local scale (1-15 km) and the remaining 29.3% were exported from the main areas of chikungunya circulation in Lazio to longer distances such as Rome and other geographical areas. Seventy percent of focal infections (corresponding to 38% of the total 414 cases) were transmitted within a distance of 200 m (the buffer distance adopted by the national guidelines for insecticide spraying). Two main epidemic clusters were identified, with a radius expanding at a rate of 300-600 m per month. The majority of exported cases resulted in either sporadic or no further transmission in the region. CONCLUSIONS: Evidence suggest that human mobility contributes to seeding a relevant number of secondary cases and new foci of transmission over several kilometers. Reactive vector control based on current guidelines might allow a significant number of secondary clusters in untreated areas, especially if the outbreak is not detected early. Existing policies and guidelines for control during outbreaks should recommend the prioritization of preventive measures in neighboring territories with known mobility flows to the main areas of transmission.

Transmission dynamics of the ongoing chikungunya outbreak in Central Italy: from coastal areas to the metropolitan city of Rome, 2017.

A large chikungunya outbreak is ongoing in Italy, with a main cluster in the Anzio coastal municipality. With preliminary epidemiological data, and a transmission model using mosquito abundance and biting rates, we estimated the basic reproduction number R0 at 2.07 (95% credible interval: 1.47-2.59) and the first case importation between 21 May and 18 June 2017. Outbreak risk was higher in coastal/rural sites than urban ones. Novel transmission foci could occur up to mid-November.

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.

An easier life to come for mosquito researchers: field-testing across Italy supports VECTRACK system for automatic counting, identification and absolute density estimation of Aedes albopictus and Culex pipiens adults.

BACKGROUND: Disease-vector mosquito monitoring is an essential prerequisite to optimize control interventions and evidence-based risk predictions. However, conventional entomological monitoring methods are labor- and time-consuming and do not allow high temporal/spatial resolution. In 2022, a novel system coupling an optical sensor with machine learning technologies (VECTRACK) proved effective in counting and identifying Aedes albopictus and Culex pipiens adult females and males. Here, we carried out the first extensive field evaluation of the VECTRACK system to assess: (i) whether the catching capacity of a commercial BG-Mosquitaire trap (BGM) for adult mosquito equipped with VECTRACK (BGM + VECT) was affected by the sensor; (ii) the accuracy of the VECTRACK algorithm in correctly classifying the target mosquito species genus and sex; (iii) Ae. albopictus capture rate of BGM with or without VECTRACK. METHODS: The same experimental design was implemented in four areas in northern (Bergamo and Padua districts), central (Rome) and southern (Procida Island, Naples) Italy. In each area, three types of traps-one BGM, one BGM + VECT and the combination of four sticky traps (STs)-were rotated each 48 h in three different sites. Each sampling scheme was replicated three times/area. Collected mosquitoes were counted and identified by both the VECTRACK algorithm and operator-mediated morphological examination. The performance of the VECTRACK system was assessed by generalized linear mixed and linear regression models. Aedes albopictus capture rates of BGMs were calculated based on the known capture rate of ST. RESULTS: A total of 3829 mosquitoes (90.2% Ae. albopictus) were captured in 18 collection-days/trap/site. BGM and BGM + VECT showed a similar performance in collecting target mosquitoes. Results show high correlation between visual and automatic identification methods (Spearman Ae. albopictus: females = 0.97; males = 0.89; P < 0.0001) and low count errors. Moreover, the results allowed quantifying the heterogeneous effectiveness associated with different trap types in collecting Ae. albopictus and predicting estimates of its absolute density. CONCLUSIONS: Obtained results strongly support the VECTRACK system as a powerful tool for mosquito monitoring and research, and its applicability over a range of ecological conditions, accounting for its high potential for continuous monitoring with minimal human effort.

Increasing situational awareness through nowcasting of the reproduction number.

Background: The time-varying reproduction number R is a critical variable for situational awareness during infectious disease outbreaks; however, delays between infection and reporting of cases hinder its accurate estimation in real-time. A number of nowcasting methods, leveraging available information on data consolidation delays, have been proposed to mitigate this problem. Methods: In this work, we retrospectively validate the use of a nowcasting algorithm during 18 months of the COVID-19 pandemic in Italy by quantitatively assessing its performance against standard methods for the estimation of R. Results: Nowcasting significantly reduced the median lag in the estimation of R from 13 to 8 days, while concurrently enhancing accuracy. Furthermore, it allowed the detection of periods of epidemic growth with a lead of between 6 and 23 days. Conclusions: Nowcasting augments epidemic awareness, empowering better informed public health responses.