The evolution of Covid-19 in Italy during the Summer 2020: analysis and interpretation of an unpredicted rest period ==================================================================================================================== * Giuseppe De Natale * Lorenzo De Natale * Claudia Troise * Renato Somma * Vito Marchitelli * Antonio Coviello * Karen G. Holmberg ## Abstract The coronavirus (COVID-19) pandemic was particularly invasive in Italy during the period of March to the end of April 2020 then displayed a significant decrease both in the number of infections and in the seriousness of illness throughout the summer of 2020. In this discussion, we measure the seriousness of the disease by the ratio of Intensive Care Units (ICU) spaces occupied by COVID-19 patients and the number of still Active Cases (AC) each month from April to August 2020. We also use the ratio between the number of Deaths (D) and the number of Active Cases. What clearly emerges, from rigorous statistical analysis, is a progressive decrease of both the ratios, indicating progressive mitigation of the disease. This is particularly evident when comparing March-April with July-August; during the summer period the two ratios have become roughly 18 times lower. We test such sharp decreases against possible bias in counting active cases, and we confirm their statistical significance. We then interpret such evidence in terms of the well-known seasonality of the human immune system and the virus-inactivating effect of stronger UV rays in the summer. ## Introduction COVID-19 had devastating effects in the months of March-May 2020 in Europe. The CFR (Case Fatality Ratio) in European countries (updated on May 26, 2020) reached peaks close to 19% in France, about 16% in Belgium, and around 14% in Italy, UK and Hungary (Oke and Henegan, 2020). De Natale et al. (2020) discuss the reasons for the high CFR in Italy and these results could potentially be applied to other European countries with very high CFR rates. The main cause, recently confirmed by widespread randomized serological tests in Italy, was in fact the large underestimation of the true number of infections during the peak of pandemic. Recent studies indicate that the true number of infected people in Italy was around 1.5 million people, i.e. about six times the tested confirmed cases (ISTAT, 2020). With such a correction, the Infection Fatality Rate (IFR), which represents the true lethality of this infection, drops to about 2.3. As occurred in Italy, it is likely that other European countries experienced a higher infection rate than was understood at the time. Another problem, experienced in the most severely impacted Italian region Lombardy was the near-collapse of the health infrastructure accompanied by crisis management errors during the infection peak (De Natale et al., 2020). Lombardy had the highest CFR in Italy, at close to 20%. Other countries with very high CFR likely had similar contexts. It is clearly evidenced that countries with a very strong health system, like Germany, are characterised by much lower CFR (Oke and Henegan, 2020). Starting in May 2020, COVID-19 seemed to lose much of its severity. This was very clearly evidenced to medical staff with direct experience with patients in the main hospitals and raised lively debates in Italy that were reported (besides Italian media) also by International press agencies (Reuters, 2020). Such limited, clinical observations raised political-social discussions in Italy about the need to continue strict containment measures. The evolution of the infection transfers and consequent illnesses during the summer and after the relaxation of the lockdown and other containment measures was far milder than expected by epidemiological forecasts (Vespignani et al., 2020; Vollmer et al., 2020). In this paper, we statistically analyse data of ICU occupation and deaths due to COVID-19 as related to the number of active cases from the end of March to August. When rigorously tested, the ratios of ICU occupation to active cases and the ratios of deaths to active cases show significant changes, thus indicating a change in the evolution of the illness from spring to summer 2020. The likely implications of such changes are then interpreted and discussed, taking into account the possible factors affecting the illness: weakening of the virus, counting bias, and changes in the immune response. The results and interpretation are then discussed in light of a possible forecast of what kind of evolution we could expect in the coming months (Autumn-Winter 2020-2021). Finally, we discuss the implications of these observations, made in Italy, for the more general worldwide situation. ## Data analysis We study the global, clinical evolution of COVID-19 in Italy using the data of ICU occupation numbers, deaths, and active cases in different periods. As De Natale et al. (2020) point out, ICU numbers and deaths are rather robust data whereas recording of active cases can be strongly biased by heterogeneous testing, loss of asymptomatic cases, etc. Recently, the first results of a massive testing campaign to randomly check the percentage of people expressing antibodies to SARS-CoV-2 were released (ISTAT report, 2020). Such tests confirm, as first hypothesized by De Natale et al. (2020), that the number of infected people was about six times larger than indicated by official tests, reaching about 1.5 million people (instead of about 250.000 officially tested positive). In the following analyses, we first assume that such high underestimation of active cases has been almost constant during the analysed period; then, we test our results with respect to the maximum bias implied by such an underestimation. Data on active cases, ICU occupation numbers, and deaths in Italy are from the Department of Italian Civil Protection Repository (2020). Here, we report the time evolution of the ratio between the number of people in ICU and the total number of ‘active’ cases (i.e., total less recovered and deaths, at that time), indicated by ICU/AC. In order to make our estimates more robust, we choose to consider another important and robust indicator: the number of deaths. We therefore also use the ratio between the number of deaths and the number of active cases, indicated as D/AC. In order to get more accurate estimates, we also considered the average time lag from COVID-19 confirmation (the actual data on active infection we have) and the ICU admission as from COVID-19 confirmation to death. According to Wilson et al. (2020), the average time lag between infection confirmation and ICU admission is 6 days whereas the average time lag between infection confirmation and death is 13 days. For this reason, we shift the median day in which we compute the active cases to 6 days before the median day of the ICU number computation; we accordingly shift the median day to compute the number of deaths to 7 days after the ICU median day (so that there are 13 days between the days of active cases computation and the corresponding days of deaths computation). As April 3, 2020 was the day of maximum ICU occupation in Italy for COVID-19 cases we used the third day of each month, from April to August, as the median day to consider ICU occupation numbers. The median days for considering active cases and deaths are chosen accordingly with the described shifts: so that for active cases we take the 27th or 28th of the month before (depending if it has 30 or 31 days), and for deaths we take the 10th of each month considered. In order to obtain more robust estimates of the various data, we choose to average the data during 7 days around each median day (considering also 3 days before and 3 days after the median day). We then computed the two quantities - ICU/AC and D/AC - as the respective averages for each month. The results, for months from April to August, are reported in Table1 and Figure1. The data in Figure 1 show the statistical uncertainties at 95% (2 standard deviations) and indicate that both the ratios sharply decrease, with decreasing values clearly separated, well above the statistical uncertainty. Another way to look at such marked change is to consider the ratios of ICU/AC and D/AC computed for April to each one of the following months. View this table: [Table 1.](http://medrxiv.org/content/early/2020/09/02/2020.08.29.20184127/T1) Table 1. Ratios ICU/AC and D/AC. Indicated uncertainties are statistical errors computed at 95% probability level (2σ). ![Fig.1a](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/09/02/2020.08.29.20184127/F1.medium.gif) [Fig.1a](http://medrxiv.org/content/early/2020/09/02/2020.08.29.20184127/F1) Fig.1a Ratio between ICU occupation and active cases in each month, from April to August. Error bars are as in table 1 (statistical errors at 95% probability level, i.e. two standard deviations). ![Fig.1b](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/09/02/2020.08.29.20184127/F2.medium.gif) [Fig.1b](http://medrxiv.org/content/early/2020/09/02/2020.08.29.20184127/F2) Fig.1b Ratio between deaths and active cases in each month, from April to August. Error bars are as in table 1 (statistical errors at 95% probability level, i.e. two standard deviations). The computed values of such ratios are reported in Table 2, and shown in Fig.2. ![Fig.2a](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/09/02/2020.08.29.20184127/F3.medium.gif) [Fig.2a](http://medrxiv.org/content/early/2020/09/02/2020.08.29.20184127/F3) Fig.2a Relative ratios of ICU/AC in April, with respect to the following months May to August. Error bars are as in table 1 (statistical errors at 95% probability level, i.e. two standard deviations). ![Fig.2b](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/09/02/2020.08.29.20184127/F4.medium.gif) [Fig.2b](http://medrxiv.org/content/early/2020/09/02/2020.08.29.20184127/F4) Fig.2b Relative ratios of D/AC in April, with respect to the following months May to August. Error bars are as in table 1 (statistical errors at 95% probability level, i.e. two standard deviations). View this table: [Table 2.](http://medrxiv.org/content/early/2020/09/02/2020.08.29.20184127/T2) Table 2. Relative ratios ICU/AC and D/AC between April and following months. Indicated uncertainties are statistical errors computed at 95% probability level (2σ). As is clear, both the ratios progressively and rapidly increase: from April to June, the relative ratio (ICU/CA)April/(ICU/CA)June increases of a factor 7.9, and the relative ratios (D/CA)April/(D/CA)June of a factor 6.2; from April to August, the relative ratio (ICU/CA)April/(ICU/CA)August increases of a factor 17.6 and the relative ratio (D/CA)April/(D/CA)August of a factor 16.7. These increases are exceptionally marked. Assuming the problems of estimating the true number of active cases are constant in time they unequivocally indicate that the illness became progressively milder during the summer months. The next step, however, is to consider the maximum possible bias due to heterogeneity in the counting procedures of infected people. As we said, we are now able to consider the massive underestimation of infection rates that occurred in Italy in the first months of the pandemic; this error indicates a factor of six more people infected than officially tested positive (ISTAT report, 2020). Then, in order to compute the maximum bias such a fluctuating underestimation could imply in our computations we assume that all of the underestimation occurred in the months of March-April and that in the last summer months we succeeded in testing all COVID-19 positive cases. Such an assumption divides the March-April ratios (ICU/AC and D/AC) by a factor 6 to test the obtained values against the values computed for August (with data from July-August). We then performed a rigorous test of hypothesis to determine if the observed increase of the ratios in August, with respect to the ratios in April though decreased of a factor 6, is significant. We applied the well-known Student’s test (Spiegel, 1992). The Student’s test starts by assuming the following formula for the t-variable (Spiegel, 1992): ![Formula][1] where: ![Graphic][2] and ![Graphic][3] are the average values of ICU/AC (or D/AC) for April and August respectively, and ![Graphic][4] and ![Graphic][5] are the respective variances as computed from the samples. In the assumption of the null hypothesis, *μ*1 and *μ*2, the true average values of the two samples are equal, and the value of the t-variable can be tested with respect to the Student’s distribution. For our data, the values result t=5.00 for the ratios ICU/AC, and t=3.7 for the ratios D/AC, which are both largely out of 99% probability limit for the Student’s distribution (respectively, out of 99.9% and of 99.5%). So, we can state that the decrease of the ratios ICU/AC and D/AC, from April to August, is significant for both the variables, even if the maximum bias due to infection counting procedures is assumed, with less than 0.5% probability of being wrong. ## Discussion and Conclusions Interpreting the very clear evidence that the disease has become progressively less severe, particularly in July-August, is not as simple as it would appear. The simplest explanation could be that the virus itself, which is continuously mutating and adapting to the host (Geoghegan et al., 2018), lost much of its lethality. With COVID-19, however, we should consider that the severity of the disease is mostly determined by the response of the human immune system (Song et al., 2020). The scientific literature currently provides no clear evidence for the virus becoming substantially mutated and less aggressive. There is evidence for some mutations which actually made the disease less severe (Barnaby et al., 2020; Benedetti et al., 2020); however, there is no indication that such mutations currently have a large diffusion (or have increased compared to past months) in order to significantly mitigate the disease. Hence, the hypothesis that the observed mitigation is due to the large diffusion of a significant virus mutation lacks scientific basis. It is well known that infectious diseases (Fishman, 2007) as well as auto-immune diseases are characterized by an evident seasonality (Iikuni et al., 2007; Moltchanova et al., 2009). The main reason for this appears to be the seasonal variation of the human immune system response, which has been assessed even in terms of gene expression (Dopico et al., 2015). In particular, during the summer the immune system response is more effective and less inflammatory. It is proven that during European winters the human immune system has a marked pro-inflammatory character, with increased levels of soluble Il-6 receptor and C reactive protein (Dopico et al., 2015). An inflammatory response by the immune system, with cytokine storms, has been recognised as the main factor leading to lung and/or other organs failure and death (e.g. Song et al., 2020). Hence, besides widespread evidence that almost all the flu-like epidemics are strongly dampened during the summer, there is a very fundamental reason because this is even more so in the case of COVID-19: in the acute phase, COVID-19 behaves like an autoimmune syndrome and so is particularly sensitive to the seasonality of the immune response (Dopico et al., 2015; Iikuni et al., 2007). Other researchers previously noted some seasonal/climate effects (e.g. Scafetta, 2020; Carleton and Meng, 2020). In addition, it has been proven the summer sunlight rapidly inactivates Sars-Cov-2 (Ratnesar-Shumate et al., 2020). It is important to highlight that other explanations, based on the possible bias due to inhomogeneous counting and/or mean age of the infected people in the different periods, as often claimed also on the basis of a relatively younger age inferred for recently recorded infections (ISS, 2020), do not appear to be sufficient to explain the significant decrease of the ICU/AC and D/AC in the summer months. In addition to our demonstration here relative to the possible bias of inhomogeneous counting, there is compelling evidence (Davies et al., 2020) that susceptibility to the virus of population older than 20 years is more than double with respect to younger population. Such evidence makes obviously very unlikely to imagine that only very young people in Italy are being infected, and such infection does not spread also into the older population. The most likely explanation of the inferred lower mean age of recent infections is that, actually, due to the predominance of rather asymptomatic cases, the testing is mostly operated to control young people coming back from summer vacations. Then, a very strong decrease of the seriousness of the disease in the summer months appears the only realistic explanation of data. What we have described in the Italian case seems to also be able to explain the comparatively low lethality of the virus in recent months observed even in countries experiencing a very large number of infections. Countries like the USA, India, Brazil, and almost all the South American and North African countries, in which epidemics are active, are experiencing much lower CFR (Case Fatality Rates) compared to what European countries experienced in the spring of 2020 (Oke and Henegan, 2020). Among them, the rate (CFR=0.8) in Israel is extremely low. Except in few cases (like Israel, for instance) it is not probable that all of these countries are recording the real number of infection cases with much higher precision than the European countries. The largely mitigated effect of COVID-19 epidemics in Italy, despite the acute situation in the months of March-April (De Natale et al., 2020), strongly validates its management of the epidemics and the Italian lockdown strategy. We should note, in fact, that just after the first reopening, in the middle of May 2020, epidemiological studies predicted a very large increase of infections, severe cases, and deaths (Vespignani et al., 2020; Vollmer et al., 2020). Instead, despite the complete reopenings of June 3, 2020, such dire forecasts have been completely discounted, confirming the appropriateness of the response and progression from complete lockdown to gradual and then complete re-opening. As we conclude here, the recent mitigation of the illness is mostly due to the seasonality of immune response, which is more effective and less inflammatory in the summer, combined with the strongly germicidal power of the summer sun’s UV rays. It remains to be seen what the practical implications are for future procedures in order to avoid new epidemic waves. Obviously, despite several media (and some doctors with direct clinical experience) claims that the virus has lost power, if the illness mitigation is only attributable to the summer season we should be very careful to continue the ‘mild’ actions of containment (social distancing, face masks in close environments, etc.) in order to avoid, in autumn and winter, the possible escalation of epidemics with renewed need for partial or total lockdowns. It is also possible that, superimposed on the seasonal effects, there is also a weakening of the virus due to natural mutations and adaptations to the host; however, we have no clear evidence for this in scientific literature. An eventual weakening of the SARS-CoV-2 virus may be evident in autumn 2020, when the concomitant seasonal effect will no longer apply. Before then, it will be good practice to maintain all the precautions to contain the epidemics. ## Data Availability All data used are in public repository, both Italian and International * Received August 29, 2020. * Revision received August 29, 2020. * Accepted September 2, 2020. * © 2020, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/) ## References 1. Young, B.E.; Fong, S.W; Chan; Y-H; Mak, T-M.; Ang, L.W.; Anderson, D.E.; Lee, C.Y-P.; Amrun, S.N.; Lee, B.; Goh, Y.S.; Su, Y.C.F.; E Wei, W.E.; Kalimuddin, S., Chai, L.Y.A.; Pada, S.; Tan, S.Y.; Sun, L.; Parthasarathy, P.; Chen, Y.Y.C.; Barkham, T.; Lin, R.T.P.; Maurer-Stroh, S.; Leo, Y-S.; Wang, L-F.; Renia, L.; Lee, V.J.; Smith, G.J.D.; Lye, D.C.; Ng, L.F.N.. Effects of a major deletion in the SARS-CoV-2 genome on the severity of infection and the inflammatory response: an observational cohort study.The Lancet 2020.DOI: 10.1016/S0140-6736(20)31757-8. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0140-6736(20)31757-8&link_type=DOI) 2. Benedetti, F.; Snyder, G.; Giovanetti, M.; Angeletti, S.; Gallo, R.; Ciccozzi, M.; Zella, D. Emerging of a SARS-CoV-2 viral strain with a deletion in nsp1. 2020 doi:10.21203/rs.3.rs-62592/v1 (preprint). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.21203/rs.3.rs-62592/v1&link_type=DOI) 3. Carleton, T.; Meng, K.C. Causal empirical estimates suggest COVID-19 transmission rates are highly seasonal, medRxiv, 2020 doi: [https://doi.org/10.1101/2020.03.26.20044420](https://doi.org/10.1101/2020.03.26.20044420) 4. Davies, N.G.; Klepac, P.; Liu, Y.; Prem, K.; Jit, M.; CMMID COVID-19 working group, Eggo, R.M. Age-dependent effects in the transmission and control of COVID-19 epidemics. Nat Med 2020, 26,1205-1211 [https://doi.org/10.1038/s41591-020-0962-9](https://doi.org/10.1038/s41591-020-0962-9) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41591-020-0962-9&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32546824&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F02%2F2020.08.29.20184127.atom) 5. De Natale, G.; Ricciardi, V.; De Luca, G.; De Natale, D.; Di Meglio, G.; Ferragamo, A.; Marchitelli, V.; Piccolo, A.; Scala, A.; Somma, R.; Spina, E.; Troise, C. The COVID-19 Infection in Italy: A Statistical Study of an Abnormally Severe Disease.J. Clin. Med. 2020, 9, 1564. 6. Department of Italian Civil Protection Repository. Available online: [http://www.protezionecivile.gov.it/attivita-rischi/rischio-sanitario/emergenze/coronavirus](http://www.protezionecivile.gov.it/attivita-rischi/rischio-sanitario/emergenze/coronavirus) (accessed on 4 August 2020). 7. Dopico, X.C.; Evangelou, M.; Ferreira, R.C.; Guo, H.; Pekalski, M.L.; Smyth, D.J.; Cooper, N.; Burren, O.S.; Fulford, A.J.; Hennig, B.J.; Prentice, A.M.; Ziegler, A-G.; Bonifacio, E.; Wallace, C.; Todd, J.A.., Widespread seasonal gene expression reveals annual differences in human immunity and physiology. Nat Commun 2015, 6, 7000. [https://doi.org/10.1038/ncomms8000](https://doi.org/10.1038/ncomms8000). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ncomms8000&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25965853&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F02%2F2020.08.29.20184127.atom) 8. Fisman, D. N. Seasonality of infectious diseases. Annu. Rev. Public Health 2007, 28, 127-143 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1146/annurev.publhealth.28.021406.144128&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17222079&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F02%2F2020.08.29.20184127.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000246436700010&link_type=ISI) 9. Geoghegan, J.L.; Holmes, E.C. The phylogenomics of evolving virus virulence. Nat Rev Genet 2018, 19, 756-769. [https://doi.org/10.1038/s41576-018-0055-5](https://doi.org/10.1038/s41576-018-0055-5) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41576-018-0055-5&link_type=DOI) 10. Iikuni, N.; Nakajima, A.; Inoue, E.; Tanaka, E.; Okamoto, H.; Hara, M.; Tomatsu, T.; Kamatani, N.; Yamanaka, H. What's in season for rheumatoid arthritis patients? Seasonal fluctuations in disease activity. Rheumatology 2007, 46, 846-848. doi: 10.1093/rheumatology/kel414 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/rheumatology/kel414&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17264092&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F02%2F2020.08.29.20184127.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000246121100024&link_type=ISI) 11. ISS Epidemia COVID-19 published August 18, 2020 pp.22 Available online: [https://www.epicentro.iss.it/coronavirus/bollettino/Bollettino-sorveglianza-integrata-COVID-1918-agosto-2020.pdf](https://www.epicentro.iss.it/coronavirus/bollettino/Bollettino-sorveglianza-integrata-COVID-1918-agosto-2020.pdf) (In Italian, accessed on August 28, 2020) 12. ISTAT Report, Primi risultati sull’indagine di sieroprevalenza sul Sars-Cov-2. 2020. Available online:[https://www.istat.it/it/files/2020/08/ReportPrimiRisultatiIndagineSiero.pdf](https://www.istat.it/it/files/2020/08/ReportPrimiRisultatiIndagineSiero.pdf) (In Italian, accessed on August 3, 2020) 13. Moltchanova, E. V.; Schreier, N.; Lammi, N.; Karvonen, M. Seasonal variation of diagnosis of Type 1 diabetes mellitus in children worldwide. Diabet. Med. 2009, 26, 673-678. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1464-5491.2009.02743.x&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19573115&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F02%2F2020.08.29.20184127.atom) 14. Oke, J.; Heneghan, C. Global COVID-19 Case Fatality Rates; Nuffield Department of Primary Care Health Sciences: Oxford, UK, 2020; Available online: [https://www.cebm.net/global-COVID-19-case-fatality-rates/](https://www.cebm.net/global-COVID-19-case-fatality-rates/) (accessed on 25 August 2020). 15. Reuters. New coronavirus losing potency, top Italian doctor says, Thomson Reuters, World News May 31, 2020 Available online: [https://www.reuters.com/article/us-health-coronavirus-italy-virus/new-coronavirus-losing-potency-top-italian-doctor-says-idUSKBN237000](https://www.reuters.com/article/us-health-coronavirus-italy-virus/new-coronavirus-losing-potency-top-italian-doctor-says-idUSKBN237000) (accessed on August 28, 2020). 16. Ratnesar-Shumate, S.; Williams, G.; Green, B.; Krause, M.; Holland, B.; Wood, S.; Bohannon, J.; Boydston, J.; Freeburger, D.; Hooper, I.; Beck, K.; Yeager, J.; Altamura, L.A.; Biryukov, J.; Yolitz, J.; Schuit, M.; Wahl, V.; Hevey, M.; Dabisch, P. Simulated Sunlight Rapidly Inactivates SARS-CoV-2 on Surfaces, The Journal of Infectious Diseases, 2020, 222, 2, 15, 214-222, [https://doi.org/10.1093/infdis/jiaa274](https://doi.org/10.1093/infdis/jiaa274) 17. Scafetta, N. Distribution of the SARS-CoV-2 Pandemic and Its Monthly Forecast Based on Seasonal Climate Patterns. Int. J. Environ. Res. Public Health, 2020, 17, 3493. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F02%2F2020.08.29.20184127.atom) 18. Song, P.; Li, W.; Xie, J.; Hou, Y.; You, C. Cytokine storm induced by sars-cov-2. Clin. Chim. Acta 2020, 509, 280-287. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cca.2020.06.017&link_type=DOI) 19. Spiegel, M. R. Theory and Problems of Probability and Statistics, New York: McGraw-Hill, 1992, pp. 116-117 20. Vespignani et al., 2020. Available online: Il Sole 24 Ore Infodata: [https://www.infodata.ilsole24ore.com/2020/05/24/covid-19-modeling-italy-quattro-scenari-prevedere-contagio/](https://www.infodata.ilsole24ore.com/2020/05/24/covid-19-modeling-italy-quattro-scenari-prevedere-contagio/) (accessed on June 1, 2020). 21. Vollmer, M.A.C et al. Report 20: Using mobility to estimate the transmission intensity of COVID-19 in Italy: A subnational analysis with future scenarios Technical Report; Imperial College COVID-19 Response Team: London, UK, 2020. Available online: [https://www.imperial.ac.uk/mrc-global-infectious-disease-analysis/covid-19/report-20-italy/](https://www.imperial.ac.uk/mrc-global-infectious-disease-analysis/covid-19/report-20-italy/) (accessed on August 28, 2020) 22. WHO (World Health Organization).Rolling Updates on Coronavirus Disease (COVID-19); WHO: Geneve, Switzerland, 2020. Available online: [https://www.who.int/emergencies/diseases/novel-coronavirus-2019?gclid=Cj0KCQjw-af6BRC5ARIsAALPIlXu7eg96lFx3xa-l\_DLkvSmCgRFGfB9FKIE6J8oLurUuiX80JiS\_5kaAqZGEALw\_wcB](https://www.who.int/emergencies/diseases/novel-coronavirus-2019?gclid=Cj0KCQjw-af6BRC5ARIsAALPIlXu7eg96lFx3xa-l_DLkvSmCgRFGfB9FKIE6J8oLurUuiX80JiS_5kaAqZGEALw_wcB) (accessed on August 28, 2020) 23. Wilson, N.; Kvalsvig, A.; Barnard, L.; Baker, M.G. Case-Fatality Risk Estimates for COVID-19 Calculated by Using a Lag Time for Fatality. Emerg. Infect. Dis. 2020, 26(6): 1339–1441. [https://dx.doi.org/10.3201/eid2606.200320](https://dx.doi.org/10.3201/eid2606.200320) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3201/eid2606.200320&link_type=DOI) [1]: /embed/graphic-7.gif [2]: /embed/inline-graphic-1.gif [3]: /embed/inline-graphic-2.gif [4]: /embed/inline-graphic-3.gif [5]: /embed/inline-graphic-4.gif