- Open Access
The association between first and second wave COVID-19 mortality in Italy
BMC Public Health volume 21, Article number: 2069 (2021)
The relation between the magnitude of successive waves of the COVID-19 outbreak within the same communities could be useful in predicting the scope of new outbreaks.
We investigated the extent to which COVID-19 mortality in Italy during the second wave was related to first wave mortality within the same provinces. We compared data on province-specific COVID-19 2020 mortality in two time periods, corresponding to the first wave (February 24–June 30, 2020) and to the second wave (September 1–December 31, 2020), using cubic spline regression.
For provinces with the lowest crude mortality rate in the first wave (February–June), i.e. < 22 cases/100,000/month, mortality in the second wave (September–December) was positively associated with mortality during the first wave. In provinces with mortality greater than 22/100,000/month during the first wave, higher mortality in the first wave was associated with a lower second wave mortality. Results were similar when the analysis was censored at October 2020, before the implementation of region-specific measures against the outbreak. Neither vaccination nor variant spread had any role during the study period.
These findings indicate that provinces with the most severe initial COVID-19 outbreaks, as assessed through mortality data, faced milder second waves.
The COVID-19 pandemic in many countries has been characterized by waves of infection. These recurrent local exacerbations of the pandemic present an opportunity to study the spread of the virus within a population. Key features of the pandemic are still not well understood, such as the susceptibility of the population to subsequent waves after the first outbreak, the threshold for herd immunity, the role of superspreaders [1,2,3,4,5,6,7] as well as of meteorological and environmental factors [8,9,10,11].
We compared the two COVID-19 waves within Italy, where the geographical distribution of the SARS-CoV-2 infection spread and the COVID-19 incidence was uneven during the first wave of the pandemic . From an analysis within provinces of the official case counts during the first and second wave of SARS-CoV-2 infection, we found evidence of an association in the occurrence between the two periods. When incidence in the first wave was low (< 500 cases/100,000/day), the second wave incidence tended to be positively correlated, whereas a high first-wave incidence was strongly and inversely correlated with the second wave incidence. This observation raised the possibility that greater spread at the beginning of the pandemic could have induced some protection at the population level, resulting in a milder second wave, despite low levels (10–20%) of anti-SARS-CoV-2 antibody seroprevalence even in the areas most severely hit.
One problem with our earlier finding is that the incidence of SARS-CoV-2 infection is difficult to measure. It depends on the severity of the clinical course, the availability of molecular tests and the policy underlying their use at the population level. In Italy, for instance, the implementation of nasopharyngeal swabs and molecular tests to detect SARS-CoV-2 infection in the first wave was restricted primarily to symptomatic patients with suspected COVID-19. In the summer of 2020, however, the public-health policy was modified to extend testing to contacts of cases and to individuals without clinical symptoms, nearly tripling the number of daily molecular tests from the first to the second wave time periods . An additional change of policy, starting on November 6, 2020, during the second wave, was the mandatory adoption of differential area-specific public-health measures [14, 15].
The area-specific relation between the two waves may be better assessed by a more objective measures of COVID-19 spread such as mortality rate, instead of incidence rates that are subject to swab testing policy and availability. Recently, COVID-19 death data for Italy in 2020 have become available with greater geographic detail, with figures available for the 107 provinces rather than the 20 larger regions within Italy. We used these mortality data to revisit the relation of the dynamics of the first and second waves of the outbreak.
Study population and outcome
We used data made freely available from public sources. We accessed 2020 monthly province-specific COVID-19 deaths as reported by the Italian National Institute of Statistics – ISTAT . Provinces are the administrative entities that are intermediate between smaller municipalities and larger regions within Italy. Within 2020, we examined two time frames: we considered the first wave to be the period from the pandemic onset (February 24) to June 30, when the number of new infection cases had substantially dropped after the March–April peak. From July to the end of August the number of new cases remained low, mostly related to screening of subjects returning from vacation with swab testing. Cases increased again beginning in early September. Therefore, we considered the second wave to be the period from September 1 to December 31, 2020. During these periods, neither circulation of virus variants nor implementation of SARS-CoV-2 vaccination at the population level had yet started in Italy. Only a handful of Italians had been vaccinated by the end of December, but large scale administration started in January 2021 [17, 18].
Based on these data and the Italian population at January 1, 2020, and January 1, 2021 available at the ISTAT website , we computed wave-specific COVID-19 mortality during the first and second waves. For the former period, we used as reference population data the provincial population size as of January 1, 2020. For the second period and for the overall 2020 mortality we averaged the population size at January 1, 2020 and January 1, 2021 for each province. We retrieved the seroprevalence data for each Italian region and province determined through a national survey made available by ISTAT and the National Ministry of Health with reference to the period May 25, 2020 through July 15, 2020 .
We examined the relation between first and second wave province-specific crude mortality rates from COVID-19 by modeling the mean mortality rate using a restricted cubic spline regression , a model that fits a curvilinear pattern to the data as previously reported [22,23,24]. In particular, we fitted a restricted cubic spline model that weighted provincial data by population size, performing both a crude analysis and adjusting for potential confounders. Variables included in the model were aging index (i.e. a ratio between resident population aged ≥65 years and those aged ≤14 years), percentage commuting outside the municipality of residence on a daily basis, and percentage of dwellings occupied by only one resident (available at a provincial level from the National Institute of Statistics ). In the regression model, we used three knots at fixed percentiles (10th, 50th and 90th) of the first wave distribution, and computed a pointwise 95% confidence interval (CI) [21, 25]. We used ‘mkspline’, ‘regress’, and ‘xbrcsplinei’ routines of Stata software (v17.0, College Station-TX, 2021) for all analyses.
Table 1 and Fig. 1 lists province-specific mortality rates in the investigated periods. National cumulative mortality (cumulative incidence of death) during the first wave averaged 58.2/100,000 persons, ranging from 6.5/100,000 persons in the Southern Italy, Basilicata region, to 165/100,000 persons in Northern Italy, Lombardy region. Corresponding figures for the second wave period were 67.6/100,000 persons for the national average, ranging from 19.4/100,000 in the Calabria region to 192/100,000 persons in the Aosta Valley region. In particular, seven provinces, five of the in Lombardy region (Bergamo, Brescia, Cremona, Lodi, and Pavia) and two in Emilia-Romagna region (Parma and Piacenza) experienced the highest mortality rates during the first wave. Second wave mortality was lower than first wave, always ≤100/100,000 persons and generally lower than the regional or national averages. The provinces with the highest first wave mortality, namely Piacenza (with the overall highest mortality in Italy with 333/100,000 persons), Cremona, Lodi and Bergamo, experienced the highest absolute decrease in mortality rates.
In the spline analysis, the relation between mortality in the first and second wave was U-shaped (Fig. 2), with a direct association between the estimates of the two waves up to 88 deaths /100,000, and above that an inverse pattern, with low second wave figures for those provinces with the highest rates during the first wave (88 through 192 deaths/100,000). The latter corresponded to a mortality rate of 22/100,000/month, as only 5 deaths were recorded in February 2020 following the first disease diagnosis in Italian residents on February 20. Crude analysis, ignoring control of potential confounders such as proportion elderly, living alone, and degree of mobility, gave similar results. When we repeated the main analysis by limiting the second wave to September–October 2020, i.e. by removing the November–December period when region-specific public health and social distancing policies were first allowed and implemented, results were not substantially changed, though mortality in this truncated period was much lower (Fig. 2).
Italy was the first country where the SARS-CoV-2 infection swept out of control, following the initial outbreak in China, where it was kept under control and geographically confined. The early arrival of the pandemic in Italy allowed a longer time frame to monitor and study the behavior of the epidemic. Italy did not allow the adoption of area-specific policies to curtail the outbreak until November, using instead a nationwide approach imposed by the Italian Government. It has been suggested recently that severity of COVID-19 pneumonia in hospitalized subjects was lower during the second wave , possibly due to SARS-CoV-2 genomic variation . However, during the time frame of this study there was no evidence of spread of SARS-CoV-2 variants , thus avoiding any confounding effect from differential geographic spread of viral lineages with different transmission or virulence features. In addition, recent assessment of fatality rates between first and second waves in Italy showed that they were comparable . Nonetheless, we cannot rule out that improvement in COVID-19 treatment may have decreased mortality during the pandemic, especially in the second wave, despite the lack of fully effective treatments [29,30,31]. Finally, vaccination had no role in curbing either of the two waves analyzed here, since almost no vaccine doses were administered in Italy in 2020. Indeed, the vaccination campaign had not begun yet in the investigated period, apart from a handful of ‘demonstration’ vaccinations administered to health professionals on December 27, 2020 at the national hub of the Spallanzani Hospital in Rome. Large scale administration started in January 2021 .
In contrast to our earlier analysis that relied on SARS-CoV-2 infection incidence, when area-specific mortality was still unavailable, we used mortality from COVID-19 to gauge the severity of the waves.
Mortality is a far more reliable indicator than incidence for documenting the actual trends of the pandemic, because the detection and notification of newly-infected cases from the general population depends on the availability and implementation of molecular testing, which in turn depends on health authority guidelines and policies . Therefore, the number and type of individuals to be subject to SARS-CoV-2 testing may differ considerably across geographic areas or time periods independently of real differences in incidence, due to changing referral guidelines (such as the presence of COVID-19 symptoms) and availability of molecular tests. In Italy, first wave health policy allowed COVID-19 molecular testing only for COVID-19 symptomatic individuals, while in the second wave much more testing occurred, with systematic extension to asymptomatic or weakly symptomatic individuals and therefore a great increase in number of notifications. As a consequence, the ‘official’ SARS-CoV-2 infection incidence strongly increased in the second wave. This increase, however, could have been an artifact due to the aforementioned changes in testing policy and availability, especially in light of the roughly comparable numbers of deaths in the two waves. However, incidence is less reliable compared to mortality in documenting the real trend of the pandemic, also due to its inherent differences compared with SARS-CoV-2 infection incidence and with the hybrid pattern of the related assessments at the population level, due to changes in molecular testing availability and policy.
We were able to confirm the previously reported incidence patterns for SARS-CoV-2 in Italy during 2020 . The key findings were a positive association between the two waves when the epidemic was not particularly strong in the first part of the year, and a negative association when the first wave was severe. This finding was similar when using a different time frame for the second wave to account for the possible influence of local policies. Local policies were permitted beginning as of November 6, 2020. In addition, the pattern we investigated is based on ‘real’ data, compared to studies presenting predictive modeling [34, 35], thus strengthening our findings.
Provinces that were relatively lightly affected by the first wave had a similar experience in the second wave, with rates closely mirroring those experienced during the first phase of the outbreak. In contrast, the Italian provinces severely hit by the first wave, such as Bergamo, Lodi, Cremona and Piacenza, experienced a marked reduction of mortality in the second wave, and presumably in the spread of COVID-19. The reasons for such an inverse pattern, however, were not directly identifiable through this study and may only be hypothesized. One possibility is a higher likelihood of immunity for superspreaders and individuals with high mobility and propensity to transmission during the first wave . Another possible explanation is the development of an additional cellular immunity against SARS-CoV-2 in the most severely hit population during the first wave [37, 38]. The provinces most severely affected by the first wave showed a population-based seroprevalence in the 10–25% range, far below the estimated threshold for herd immunity of 50–70% . It is possible that the prevalence of immunity, despite being below the herd immunity threshold, was high enough in provinces with a severe first wave to interact with behavioral factors and lessen the intensity of the second wave.
Availability of data and materials
We used data made freely available from public sources. All data generated or analysed during this study are included in this published article.
Kault D. Superspreaders, asymptomatics and COVID-19 elimination. Med J Aust. 2020;213(10):447–8 e1. https://doi.org/10.5694/mja2.50835.
Xia Y, Zhong L, Tan J, Zhang Z, Lyu J, Chen Y, et al. How to understand "herd immunity" in COVID-19 pandemic. Front Cell Dev Biol. 2020;8:547314. https://doi.org/10.3389/fcell.2020.547314.
Okell LC, Verity R, Watson OJ, Mishra S, Walker P, Whittaker C, et al. Have deaths from COVID-19 in Europe plateaued due to herd immunity? Lancet. 2020;395(10241):e110–e1. https://doi.org/10.1016/S0140-6736(20)31357-X.
Randolph HE, Barreiro LB. Herd immunity: understanding COVID-19. Immunity. 2020;52(5):737–41. https://doi.org/10.1016/j.immuni.2020.04.012.
Kwok KO, Lai F, Wei WI, Wong SYS, Tang JWT. Herd immunity - estimating the level required to halt the COVID-19 epidemics in affected countries. J Inf Secur. 2020;80(6):e32–e3. https://doi.org/10.1016/j.jinf.2020.03.027.
Omer SB, Yildirim I, Forman HP. Herd immunity and implications for SARS-CoV-2 control. JAMA. 2020;324(20):2095–6. https://doi.org/10.1001/jama.2020.20892.
Papachristodoulou E, Kakoullis L, Parperis K, Panos G. Long-term and herd immunity against SARS-CoV-2: implications from current and past knowledge. Pathog Dis. 2020;78:ftaa025. https://doi.org/10.1093/femspd/ftaa025.
Marques M, Rovira J, Nadal M, Domingo JL. Effects of air pollution on the potential transmission and mortality of COVID-19: a preliminary case-study in Tarragona Province (Catalonia, Spain). Environ Res. 2021;192:110315. https://doi.org/10.1016/j.envres.2020.110315.
Filippini T, Rothman KJ, Cocchio S, Narne E, Mantoan D, Saia M, et al. Associations between mortality from COVID-19 in two Italian regions and outdoor air pollution as assessed through tropospheric nitrogen dioxide. Sci Total Environ. 2021;760:143355. https://doi.org/10.1016/j.scitotenv.2020.143355.
Cherrie M, Clemens T, Colandrea C, Feng Z, Webb DJ, Weller RB, et al. Ultraviolet a radiation and COVID-19 deaths in the USA with replication studies in England and Italy. Br J Dermatol. 2021;185(2):363–70. https://doi.org/10.1111/bjd.20093.
Filippini T, Rothman KJ, Goffi A, Ferrari F, Maffeis G, Orsini N, et al. Satellite-detected tropospheric nitrogen dioxide and spread of SARS-CoV-2 infection in northern Italy. Sci Total Environ. 2020;739:140278. https://doi.org/10.1016/j.scitotenv.2020.140278.
Vinceti M, Filippini T, Rothman KJ, Di Federico S, Orsini N. SARS-CoV-2 infection incidence during the first and second COVID-19 waves in Italy. Environ Res. 2021;197:111097. https://doi.org/10.1016/j.envres.2021.111097.
Romagnani P, Gnone G, Guzzi F, Negrini S, Guastalla A, Annunziato F, et al. The COVID-19 infection: lessons from the Italian experience. J Public Health Policy. 2020;41(3):238–44. https://doi.org/10.1057/s41271-020-00229-y.
Presidente del Consiglio dei Ministri. DPCM 6 Novembre 2020. 2020 http://www.governo.it/sites/new.governo.it/files/dpcm_20200322.pdf [Accessed August 30, 2021].
Vinceti M, Filippini T, Rothman KJ, Ferrari F, Goffi A, Maffeis G, et al. Lockdown timing and efficacy in controlling COVID-19 using mobile phone tracking. EClinicalMedicine. 2020;25:100457. https://doi.org/10.1016/j.eclinm.2020.100457.
ISTAT. COVID-19 deaths in Italy. Report March 5, 2021. https://www.istat.it/it/files//2020/03/tabella-decessi-provinciali_fonte_istat_decessi_provinciali_per_covid_fonte_ISS_5marzo.xlsx [Accessed August 30, 2021].
ELIXIR-IT. SARS-CoV-2 lineages of potential epidemiological interest in Italy. 2021 https://www.covid19dataportal.it/ [Accessed August 30, 2021].
Ministry of Health. Anti-COVID vaccination campaign. 2021 www.salute.gov.it [Accessed August 30, 2021].
ISTAT. Populations and households. Italian National Institute of Statistics. 2020 https://www.istat.it/en/population-and-households [Accessed August 31, 2021].
ISTAT. First results of the seroprevalence investigation on SARS-COV-2. Italian National Institute of Statistics. 2020 https://www.istat.it/it/files//2020/08/ReportPrimiRisultatiIndagineSiero.pdf [Accessed August 30, 2021].
Harrell FE. Regression modeling strategies: with applications to linear models, logistic regression, and survival analysis. New York: Springer; 2001. https://doi.org/10.1007/978-1-4757-3462-1.
Filippini T, Adani G, Malavolti M, Garuti C, Cilloni S, Vinceti G, et al. Dietary habits and risk of early-onset dementia in an Italian case-control study. Nutrients. 2020;12(12):3682. https://doi.org/10.3390/nu12123682.
Filippini T, Mandrioli J, Malagoli C, Costanzini S, Cherubini A, Maffeis G, et al. Risk of amyotrophic lateral sclerosis and exposure to particulate matter from vehicular traffic: a case-control study. Int J Environ Res Public Health. 2021;18(3):973. https://doi.org/10.3390/ijerph18030973.
Vinceti M, Rothman KJ, Crespi CM, Sterni A, Cherubini A, Guerra L, et al. Leukemia risk in children exposed to benzene and PM10 from vehicular traffic: a case-control study in an Italian population. Eur J Epidemiol. 2012;27(10):781–90. https://doi.org/10.1007/s10654-012-9727-1.
Orsini N. XBRCSPLINE: Stata module to tabulate differences in predicted responses after restricted cubic spline models: Boston College Department of economics; 2009.
Borghesi A, Golemi S, Carapella N, Zigliani A, Farina D, Maroldi R. Lombardy, northern Italy: COVID-19 second wave less severe and deadly than the first? A preliminary investigation. Infect Dis (Lond). 2021;53(5):370–5. https://doi.org/10.1080/23744235.2021.1884745.
Toyoshima Y, Nemoto K, Matsumoto S, Nakamura Y, Kiyotani K. SARS-CoV-2 genomic variations associated with mortality rate of COVID-19. J Hum Genet. 2020;65(12):1075–82. https://doi.org/10.1038/s10038-020-0808-9.
Filippini T, Zagnoli F, Bosi M, Giannone ME, Marchesi C, Vinceti M. An assessment of case-fatality and infection-fatality rates of first and second COVID-19 waves in Italy: COVID-19 fatality rate in Italy. Acta Biomed. 2021;92:e2021420. https://doi.org/10.23750/abm.v92iS6.12241.
COVID-19 RISK and Treatments (CORIST) Collaboration. Use of hydroxychloroquine in hospitalised COVID-19 patients is associated with reduced mortality: findings from the observational multicentre Italian CORIST study. Eur J Intern Med. 2020;82:38–47. https://doi.org/10.1016/j.ejim.2020.08.019.
Docea AO, Tsatsakis A, Albulescu D, Cristea O, Zlatian O, Vinceti M, et al. A new threat from an old enemy: reemergence of coronavirus (review). Int J Mol Med. 2020;45:1631–43. https://doi.org/10.3892/ijmm.2020.4555.
Di Castelnuovo A, Costanzo S, Antinori A, Berselli N, Blandi L, Bonaccio M, et al. Lopinavir/ritonavir and darunavir/cobicistat in hospitalized COVID-19 patients: findings from the multicenter Italian CORIST study. Front Med (Lausanne). 2021;8:639970. https://doi.org/10.3389/fmed.2021.639970.
CPD - Italian Civil Protection Department. Anti COVID-19 Vaccine Report. 2021 https://italia.github.io/covid19-dashboard-vaccini/ [Accessed August 30, 2021].
De Matteis A, Turkmen Ceylan FB, Adrianopoli M, Ertek S. The contribution of testing in the fight against COVID-19: evidence from Italy. Eur J Pub Health. 2021;31(4):913–7. https://doi.org/10.1093/eurpub/ckab064.
Qian Y, Xie W, Zhao J, Xue M, Liu S, Wang L, et al. Investigating the effectiveness of re-opening policies before vaccination during a pandemic: SD modelling research based on COVID-19 in Wuhan. BMC Public Health. 2021;21(1):1638. https://doi.org/10.1186/s12889-021-11631-w.
Ahmed DA, Ansari AR, Imran M, Dingle K, Bonsall MB. Mechanistic modelling of COVID-19 and the impact of lockdowns on a short-time scale. PLoS One. 2021;16(10):e0258084. https://doi.org/10.1371/journal.pone.0258084.
Signorelli C, Zucchi A, Tersalvi CA, Ciampichini R, Beato E, Balzarini F, et al. High seroprevalence of SARS-CoV-2 in Bergamo: evidence for herd immunity or reason to be cautious? Int J Public Health. 2020;65(9):1815–7. https://doi.org/10.1007/s00038-020-01524-x.
Braun J, Loyal L, Frentsch M, Wendisch D, Georg P, Kurth F, et al. SARS-CoV-2-reactive T cells in healthy donors and patients with COVID-19. Nature. 2020;587(7833):270–4. https://doi.org/10.1038/s41586-020-2598-9.
Mateus J, Grifoni A, Tarke A, Sidney J, Ramirez SI, Dan JM, et al. Selective and cross-reactive SARS-CoV-2 T cell epitopes in unexposed humans. Science. 2020;370(6512):89–94. https://doi.org/10.1126/science.abd3871.
Fontanet A, Cauchemez S. COVID-19 herd immunity: where are we? Nat Rev Immunol. 2020;20(10):583–4. https://doi.org/10.1038/s41577-020-00451-5.
This study was supported by the grants ‘UNIMORE FAR 2019 and 2020 Interdisciplinare’ of the Fondazione di Modena’ to Dr. Vinceti and Dr. Filippini and by the FISR 2020-COVID19 grant by the Italian Ministry of the University and the Research to Dr. Vinceti.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Vinceti, M., Filippini, T., Rothman, K.J. et al. The association between first and second wave COVID-19 mortality in Italy. BMC Public Health 21, 2069 (2021). https://doi.org/10.1186/s12889-021-12126-4
- Public health