Relationships between demographic, geographic, and environmental statistics and the spread of novel coronavirus disease (COVID-19) in Italy =========================================================================================================================================== * Alessandro Rovetta * Lucia Castaldo ## Abstract **Background** Since January 2020, the COVID-19 pandemic has raged around the world, causing nearly a million deaths and hundreds of severe economic crises. In this scenario, Italy has been one of the most affected countries. **Objective** This study investigated significant correlations between COVID-19 cases and demographic, geographical, and environmental statistics of each Italian region from February 26 to August 12, 2020. We further investigated the link between the spread of SARS-CoV-2 and particulate matter (PM) 2.5 and 10 concentrations before the lockdown in Lombardy. **Methods** All demographic data were obtained from the *AdminStat Italia* website, and geographic data were from the *Il Meteo* website. The collection frequency was one week. Data on PM2.5 and PM10 average daily concentrations were collected from previously published articles. We used Pearson’s coefficients to correlate the quantities that followed a normal distribution, and Spearman’s coefficient to correlate quantities that did not follow a normal distribution. **Results** We found significant strong correlations between COVID-19 cases and population number in 60.0% of the regions. We also found a significant strong correlation between the spread of SARS-CoV-2 in the various regions and their latitude, and with the historical averages (last 30 years) of their minimum temperatures. We identified a significant strong correlation between the number of COVID-19 cases until August 12 and the average daily concentrations of PM2.5 in Lombardy until February 29, 2020. No significant correlation with PM10 was found in the same periods. However, we found that 40 *μg*/*m*3 for PM2.5 and 50 *μg*/*m*3 PM10 are plausible thresholds beyond which particulate pollution clearly favors the spread of SARS-CoV-2. **Conclusion** Since SARS-CoV-2 is correlated with historical minimum temperatures and PM10 and 2.5, health authorities are urged to monitor pollution levels and to invest in precautions for the arrival of autumn. Furthermore, we suggest creating awareness campaigns for the recirculation of air in enclosed places and to avoid exposure to the cold. Keywords * COVID-19 * Italy * Pandemic * Epidemiology * Coronavirus-2019 ## Introduction The Chief of the World Health Organization (WHO) has declared the COVID-19 pandemic the most severe pandemic in recent human history [1]. To date, over 200 countries have been involved, with over 29 million cases and 900,000 deaths [2]. Between February and April 2020, Italy was the most affected nation in terms of the number of new cases and deaths [3]. The severity of the emergency was so important that, despite a drastic drop in infections during the summer months, it is still among the top 20 nations afflicted by the novel coronavirus [4]. On January 23, two Chinese tourists tested positive for COVID-19 near Rome [5]. However, the patients appeared to have been readily isolated, averting extensive contagion. Toward the end of February, the situation began to accelerate and fall outside the control of the authorities. From February 21, to counter the spread of the virus, the Italian government declared various lockdowns that extended around the outbreak of Lodi, in the Lombardy region. On March 10, the lockdown went into effect nationwide [6]. For these reasons, we believe Italy is one of the main sources of information for fully understanding the behavior of SARS-CoV-2. This is the first study to provide a complete and detailed history of the correlations between the spread of SARS-CoV-2 and demographic, geographic, and environmental statistics in Italy. From the analysis of our results, it was possible to highlight anomalous and/or local properties and behaviors of the virus as well as test the statistical significance of the hypotheses and scenarios proposed by other research. ## Methods We collected data on Italian demographic statistics and pollution from the *AdminStat Italia* website and a previously published article [7, 8]. We looked for significant Pearson (*R*) and Spearman (*r*) correlations between COVID-19 cases per province in each region from February 26 until August 12, 2020, and birth rate, median age, population density, death rate, old age index, population number, percentage of unmarried people, family members, growth rate, percentage of divorcees, percentage of foreigners, foreign growth rate, and percentage of widowers. The frequency of data collection was one week. For each correlation identified, we calculated the angular coefficient of the interpolating line and correlated the latter with geographical characteristics such as regions’ latitude and minimum temperatures in February and March (last 30 years of historical data). We collected geographic data from the *Il Meteo* website [9]. The survey was conducted in the following regions: Abruzzo, Calabria, Campania, Emilia-Romagna, Friuli Venezia Giulia, Lazio, Liguria, Lombardia, Marche, Piemonte, Sardegna, Sicilia, Toscana, Veneto (number of provinces > 3). In the *Results* section, we report only significant results. The week in which the correlation approached the threshold of statistical significance is reported in parentheses, for example, Abruzzo (4). Finally, we have deepened the results of a previous study on the search for correlations between the daily average concentrations of particulate matter (PM) 10 and 2.5 until February 29 and COVID-19 cases for every province in the Lombardy region from February 26 to August 12, 2020 [7]. ### Statistical analysis Each result was reported together with its standard deviation (*SD*) and p-value (*p*); we chose a statistical significance threshold of *α* ≤ 05. For each sample, kurtosis (*k*) and skewness (*s*) were calculated using Microsoft Excel 2020 software. We used the formulas ![Graphic][1] and ![Graphic][2], with n sample size, to obtain their respective standard deviations [10]. We considered data compatible with a normal distribution only when *t**k* = *k* · (24/*n*)-1/2 ≤ 1.5 and *t**s* =*s* · (6/*n*)-1/2 ≤ 1.5; here, the Pearson correlation index was deemed more reliable. When *t**k* ∈ [1.5, 3], *t**s* ≤ 3 or *t**k* ≤ 3, *t**s* ∈ [1.5, 3], we considered it appropriate to evaluate both correlations. Finally, when *t**k*, *t**s* > 3, we judged the Spearman correlation index to be more appropriate. When we highlighted linear correlations, we used the *Igor Pro 6*.*37* software to interpolate the data through the equation *y* = *a* + *bx*. In the results section, we report the mean values *R**best*, *r**best* of the correlations identified, with their 95% confidence intervals. When the coefficients exceeded the value. 700 with p ≫ .05, they were reported to specify the absence of statistical significance. To verify the importance of the correlation, we constructed a suitable correlation matrix. We have defined *pure* correlations independent from other quantities correlated with COVID-19. ## Results ### Correlation between COVID-19 cases and population number We found significant strong correlations between COVID-19 cases and population number in 60.0% of regions such as Calabria (5), Campania (1), Lazio (1), Liguria (2), Lombardy (4), Piedmont (2), Sardinia (3), Sicily (1), and Veneto (4) (*R**best* = .935, 95% *CI*: .830 – 1.000, *p**best* = .046, 95% *CI*: .006 - .040). We identified suspicious correlations in Emilia-Romagna, Marche, and Puglia (Table 1). Since correlations have occurred in many regions from the first weeks, it is likely that the virus had spread as early as January. On the contrary, Veneto and the most affected area, Lombardy, seem to have experienced a gradual contagion. The average of the angular coefficients resulting from the linear interpolations of the pairs (COVID-19 cases, population number) is *b* = .0037 (95% *CI*. 0009, .0065). View this table: [Table 1.](http://medrxiv.org/content/early/2020/09/24/2020.09.18.20196980/T1) Table 1. Pearson and Spearman correlations between COVID-19 cases and population density and between COVID-19 cases and population numbers. ### Speed of spread of COVID-19: Correlation with latitude and minimum temperatures We found a significant strong correlation between the angular coefficients *b* of the various regions and their latitude (Table 2). These data indicate the dependence of COVID-19 on geographical and/or climatic factors (*R* = .895, *p* < .0001, *r* = .874, *p* < .0001). In particular, we found a significant correlation with the historical averages (last 30 years) of the minimum temperatures in the regions (*R* = −. 576, *p* = .039; *r* = −. 629, *p* = .021 for March, *R* = −. 685, *p* = .010, *r* = −. 615, *p* = .025 for February). By narrowing the analysis to the regions that showed a net correlation, we obtained much stronger correlations (*R* = −. 849, *p* = .008, *r* = −. 940, *p* = .005 for March, *R* = −. 923, *p* = .001, *r* = −. 872, *p* = .005 for February). The same is true for the correlation with latitude (*R* = .926, *p* = .001, *r* = .886, *p* = .003). View this table: [Table 2.](http://medrxiv.org/content/early/2020/09/24/2020.09.18.20196980/T2) Table 2. Geographical and environmental data of the regions in which a correlation was found (including suspicious ones). ### Correlation between COVID-19 cases and PM10 and PM2.5 daily averages (Lombardy only) We identified a strong significant correlation between the number of COVID-19 cases until August 12 and the average daily concentrations of PM2.5 in Lombardy until February 29, 2020 (r=.76, p=.004). We found no significant correlation with PM10 in the same periods. Therefore, in the long run, the correlation with PM2.5 was more statistically incisive than that with PM10. However, in the early stages of the outbreak (until February 26, 2020), we found a correlation with PM2.5 (*r*= .63, *p* = .029) and PM10 (*r* = .72, *p* = .009). In the second week of March, the correlation with PM10 disappeared, while that with PM2.5 continued to exist until now. We identified a drastic surge in COVID-19 cases near 40 *μg*/*m*3 PM2.5 and 50*μg*/*m*3 PM10 (Figure 1); therefore, these may be thresholds beyond which particulate pollution clearly favors the spread of SARS-CoV-2. All the correlations found are statistically valid as they are not related to the other quantities analyzed. ![Figure 1.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/09/24/2020.09.18.20196980/F1.medium.gif) [Figure 1.](http://medrxiv.org/content/early/2020/09/24/2020.09.18.20196980/F1) Figure 1. COVID-19 cases – PM10 and PM2.5 scatterplot evolution from February 26 to March 18, 2020. ### Correlation between COVID-19 cases and population density We found significant strong correlations between COVID-19 cases and population density in % of regions such as Abruzzo (4), Campania (1), Lazio (1), Sicily (2), and Veneto (1) (*R**best* = .880, 95% *CI*: .757 – 1.000, *p**best* = .023, 95% *CI*: .006 - .040). Suspicious correlations have been found in Piedmont and Liguria (Table 1). The only pure correlation was that in Abruzzo since there was a correlation between the population density and number in all other cases (*R**best* = .939, 95% *CI*: .880 – 997, *p**best* = .003, 95% *CI*: .001 - .006). However, we highlighted discrepancies between the onset of the above in Sicily and Veneto. ### Plausible scenarios The results are compatible with the following scenarios: 1. The almost immediate correlations between COVID-19 cases and the number of inhabitants per province in Campania, Lazio, Sicily, Liguria, and Piedmont strongly indicate that SARS-CoV-2 was in circulation for a long time before the first confirmed case. Observing the Lombardy trend, we deduced that the virus seems to have taken four weeks to correlate with the demographic dimension. It is therefore plausible that COVID-19 spread in Italy from January 2020 (or earlier). 2. Low temperatures can weaken the immune defenses, favoring the contagion from SARS-CoV-2. 3. Low temperatures can push people to have gatherings indoors and without air circulation, favoring the spread of the novel coronavirus. 4. Low temperatures could promote the survival of the virus. 5. PM2.5 can weaken the immune defenses, favoring the contagion from novel coronavirus. 6. PM2.5 can serve as an excellent carrier for the novel coronavirus. 7. High concentrations of PM10 may have contributed to the spread of the virus by acting as a carrier. This is supported not only by the initial correlation with PM10 but also by the extremely high PM10 concentration found in the first outbreak in Lodi (67 /). 8. It is possible that SARS-CoV-2 has undergone evolutionary mutations in northern Italy (particularly in Lombardy). This agrees with the possibility that COVID-19 spread nationwide, mistaken for flu or common colds until the aforementioned mutation occurred. ## Discussion The first two confirmed cases of COVID-19 in Italy were two Chinese tourists who landed on January 19, 2020 [5]. They could travel freely around the town, making a brief stop in Parma and staying in a hotel in Rome [11]. We believe that it is statistically unlikely and morally incorrect to associate them with the arrival of the novel coronavirus in Italy. In particular, our results suggest that SARS-CoV-2 has been circulating in Italy since early January, probably mistaken for flu or the common cold. In fact, having ascertained that: * Those suffering from COVID-19 were largely asymptomatic [12] * Symptoms associated with COVID-19 are milder in children compared with adults [13] * The incubation period ranges from 2 to 14 days with an average of 5–6 days [14] * 80% of people with COVID-19 have mild symptoms [15] * Average mortality is about 1% [16, 17] * The fatality rate is proportional to age and extremely high for patients over 65 years [18] the most likely hypothesis is that the virus had been present in Italy before their arrival. All these factors drastically increase the time needed to identify and isolate a person infected with COVID-19 since: a) the asymptomatic individuals infected without their knowledge and the people with whom they have been in contact, b) children (besides the problem of asymptomaticity) showed milder symptoms than adults, inducing parents and relatives to associate them with normal colds or flu, c) the symptoms were identified up to 14 days after contagion, allowing infected subjects to infect other people unknowingly, d) in the overwhelming majority of cases, symptomatic patients showed mild symptoms not causing concern in work colleagues and relatives, and e) the fact that the death risk was extremely concentrated in the high age groups prevented an easy assessment of the extent of the epidemic since these groups were naturally more exposed to fatal phenomena, that is, until high numbers have been reached, the collective psychological impact was low. The marked correlation between the speed of the spread of the virus among the population and the minimum temperatures of each region suggests that the late autumn and winter seasons can strongly favor the pandemic. In fact, due to low external temperatures, people more frequently have indoor gatherings without recirculation of air, which creates favorable conditions for the proliferation of viruses [19]. Just as rhinoviruses, adenoviruses, and influenza viruses, the novel coronavirus may also survive better in colder and drier climates [20-23]. Furthermore, sudden changes in temperature and cold can lower immune defenses [22]. The strong and prolonged correlations we found with fine PM2.5 in the Lombardy region indicate that this type of pollution played an important role in the epidemic. This may be linked to the fact that fine particles substantially reduce immune defenses as well as increase the severity of symptoms due to damage induced in the respiratory system [24, 25]. Moreover, PM2.5 (and PM10) could act as a virus carrier [26]. In fact, from Figure 1, it can be seen that there is a drastic increase in the number of COVID-19 cases when PM2.5 and PM10 exceed approximately 40 *μg*/*m*3 and 50 *μg*/*m*3, respectively. These could be the thresholds beyond which particulate matter significantly favors the spread of SARS-CoV-2. In addition to empirical data and other studies, this hypothesis is consistent with the first outbreak in Lodi, where the PM10 average daily concentration was the highest of the month (67 *μg*/*m*3) [7, 27]. However, our results show a greater incisiveness of PM2.5 compared to PM10 in the spread of the virus in Lombardy (Figure 1). Finally, considering that in Lombardy, the correlation between COVID-19 cases and the number of inhabitants per region became significant after four weeks, the severity of symptoms was more severe than in other regions, and the basic reproduction number (*R*0) seems to have been the highest, we suggest that an evolutionary genetic mutation may have occurred in Lombardy [7]. In fact, although the genome does not appear to have changed substantially, Zhang et al. showed that even small mutations can cause significant changes in SARS-CoV-2 behavior [28]. ### Limitations Statistical correlations can provide valid support for hypotheses and theories as well as fundamental indicators of phenomena to be explored; however, they cannot demonstrate the causal nature of a phenomenon. ## Conclusion We found significant strong and lasting correlations between the spread of SARS-CoV-2 and the number of inhabitants of each region, between the spread-speed of COVID-19 and the historical minimum temperatures in the months of February and March, and between the number of COVID-19 cases and the average daily concentrations of fine PM2.5. Correlations with the average daily concentrations of PM10 were found up to the first week of March, indicating that this type of pollution also played a role in the spread of the virus, linked to exceeding specific daily peaks. Therefore, we suggest that the health authorities pay particular attention to the arrival of the winter months, not only by investing in adequate precautions but also by launching various awareness campaigns on air recirculation indoors and in avoiding exposure to the cold. In addition, it will also be necessary to carefully monitor levels of PM10 and 2.5, even by limiting the use of cars if necessary. ## Data Availability All the data necessary for carrying out this study are presented in the paper or in the articles reported in the references. * Received September 18, 2020. * Revision received September 24, 2020. * Accepted September 24, 2020. * © 2020, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NonCommercial 4.0 International), CC BY-NC 4.0, as described at [http://creativecommons.org/licenses/by-nc/4.0/](http://creativecommons.org/licenses/by-nc/4.0/) ## References 1. [1].Kaiser Family Foundation (KFF). WHO says COVID-19 pandemic worst ever global health emergency, calls on nations to do more to prevent spread. NYT, WSJ examine impacts on developing countries’ economies. July 27, 2020. Accessed: September 16, 2020. URL: [https://www.kff.org/news-summary/who-says-covid-19-pandemic-worst-ever-global-health-emergency-calls-on-nations-to-do-more-to-prevent-spread-nyt-wsj-examine-impacts-on-developing-countries-economies/](https://www.kff.org/news-summary/who-says-covid-19-pandemic-worst-ever-global-health-emergency-calls-on-nations-to-do-more-to-prevent-spread-nyt-wsj-examine-impacts-on-developing-countries-economies/) 2. [2].World Health Organization (WHO). Coronavirus disease (COVID-19) pandemic. Accessed: September 16, 2020. URL: [https://www.who.int/emergencies/diseases/novel-coronavirus-2019](https://www.who.int/emergencies/diseases/novel-coronavirus-2019) 3. [3].Coronavirus: Italy death toll tops 30,000, highest in EU. BBC News, May 8, 2020. Accessed: September 16, 2020. URL: [https://www.bbc.com/news/world-europe-52594570](https://www.bbc.com/news/world-europe-52594570) 4. [4].Worldometer. COVID-19 coronavirus pandemic. Accessed: September 16, 2020. URL: [https://www.worldometers.info/coronavirus/](https://www.worldometers.info/coronavirus/) 5. [5].Giovanetti M, Benvenuto D, Angeletti S, Ciccozzi M. The first two cases of 2019LnCoV in Italy: Where they come from? J Med Virol. 2020 May; 92(5): 518–521. DOI: 10.1002/jmv.25699 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/jmv.25699&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32022275&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F24%2F2020.09.18.20196980.atom) 6. [6].Coronavirus: Italy extends emergency measures nationwide. BBC News. Accessed: September 16, 2020. URL: [https://www.bbc.com/news/world-europe-51810673](https://www.bbc.com/news/world-europe-51810673) 7. [7].Rovetta A, Bhagavathula AS, Castaldo L. Modeling the epidemiological trend and behavior of COVID-19 in Italy. Cureus 12(8): e9884. DOI: 10.7759/cureus.9884 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.7759/cureus.9884&link_type=DOI) 8. [8].Mappe tematiche, curiosità, confronti e classifiche per i comuni, le province e le regioni sulla base di 20 indicatori socio-demografici. AdminStat Italia (Italian website). Accessed: August 12, 2020. URL: [https://ugeo.urbistat.com/AdminStat/it/it/classifiche/dati-sintesi/province/italia/380/1](https://ugeo.urbistat.com/AdminStat/it/it/classifiche/dati-sintesi/province/italia/380/1) 9. [9].Medie Climatiche Italia. Il Meteo (Italian website). Accessed: September 12, 2020. URL: [https://www.ilmeteo.it/portale/medie-climatiche?refresh_ce](https://www.ilmeteo.it/portale/medie-climatiche?refresh_ce) 10. [10].Wright BD, Herrington AJ. Problematic standard errors and confidence intervals for skewness and kurtosis. Behav Res 43, 8–17 (2011). DOI: 10.3758/s13428-010-0044-x [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3758/s13428-010-0044-x&link_type=DOI) 11. [11].Turisti cinesi a Roma contagiati dal coronavirus. L’allarme all’hotel Palatino: “Mio marito ha la febbre” (article in Italian). La Repubblica (Italian newspaper), January 31, 2020. Accessed: September 16, 2020. URL: [https://www.repubblica.it/cronaca/2020/01/31/news/coronavirus\_italia\_turisti\_cinesi-247205434/](https://www.repubblica.it/cronaca/2020/01/31/news/coronavirus_italia_turisti_cinesi-247205434/) 12. [12].Gandhi M, Yokoe DS, Havlir DV. Asymptomatic transmission, the Achilles’ heel of current strategies to control Covid-19, N Engl J Med 2020; 382:2158–2160, DOI: 10.1056/NEJMe2009758 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMe2009758&link_type=DOI) 13. [13].Bixler D, Miller AD, Mattison CP, et al. SARS-CoV-2–Associated deaths among persons aged <21 Years — United States, February 12–July 31, 2020. MMWR Morb Mortal Wkly Rep. ePub: 15 September 2020. DOI: 10.15585/mmwr.mm6937e4 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.15585/mmwr.mm6937e4&link_type=DOI) 14. [14].World Health Organization (WHO). Coronavirus disease 2019 (COVID-19) Situation Report – 73. URL: [https://www.who.int/docs/default-source/coronaviruse/situation-reports/20200402-sitrep-73-covid-19.pdf](https://www.who.int/docs/default-source/coronaviruse/situation-reports/20200402-sitrep-73-covid-19.pdf) 15. [15].I’ve heard 80% of people with coronavirus (COVID-19) have mild symptoms. What are mild symptoms? WebMD, August 11, 2020. Accessed: September 17, 2020. URL: [https://www.webmd.com/lung/qa/ive-heard-80-of-people-with-coronavirus-covid19-have-mild-symptoms-what-are-mild-symptoms](https://www.webmd.com/lung/qa/ive-heard-80-of-people-with-coronavirus-covid19-have-mild-symptoms-what-are-mild-symptoms) 16. [16].Fauci AS, Lane CH, Redfield RR. COVID-19 — Navigating the uncharted, N Engl J Med 2020; 382:1268–1269, DOI: 10.1056/NEJMe2002387 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMe2002387&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32109011&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F24%2F2020.09.18.20196980.atom) 17. [17].Mallapaty S. How deadly is the coronavirus? Scientists are close to an answer, Nature 582, 467–468 (2020), DOI: 10.1038/d41586-020-01738-2 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/d41586-020-01738-2&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32546810&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F24%2F2020.09.18.20196980.atom) 18. [18].COVID-19 hospitalization and death by age, Centers for Disease Control and Prevention (CDC), August 18, 2020. Accessed: September 17, 2020. URL: [https://www.cdc.gov/coronavirus/2019-ncov/covid-data/investigations-discovery/hospitalization-death-by-age.html](https://www.cdc.gov/coronavirus/2019-ncov/covid-data/investigations-discovery/hospitalization-death-by-age.html) 19. [19].Why is the flu more common in the winter? WebMD, July 20, 2020. Accessed: September 17, 2020. URL: [https://www.webmd.com/cold-and-flu/qa/why-is-the-flu-more-common-in-the-winter](https://www.webmd.com/cold-and-flu/qa/why-is-the-flu-more-common-in-the-winter) 20. [20].Elert M. Why is there a winter flu season? Popular Science, January 17, 2013. URL: [https://www.popsci.com/science/article/2013-01/fyi-why-winter-flu-season/](https://www.popsci.com/science/article/2013-01/fyi-why-winter-flu-season/) 21. [21].Feltman R. Why winter is cold and flu season—and what you can do about it. Popular Science, January 4, 2018. URL: [https://www.popsci.com/cold-flu-season-sick-winter/](https://www.popsci.com/cold-flu-season-sick-winter/) 22. [22].Study shows why the flu likes winter. New York Times, December 5, 2007. Accessed: September 17, 2020. URL: [https://www.nytimes.com/2007/12/05/health/research/05flu.html](https://www.nytimes.com/2007/12/05/health/research/05flu.html) 23. [23].Coronavirus: How bad will winter really be? BBC News, July 27, 2020. Accessed: September 17, 2020. URL: [https://www.bbc.com/news/health-53515077](https://www.bbc.com/news/health-53515077) 24. [24].Xing YF, Xu YH, Lian YX. The impact of PM2.5 on the human respiratory system, J Thorac Dis. 2016 Jan; 8(1): E69–E74. DOI: 10.3978/j.issn.2072-1439.2016.01.19 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3978/j.issn.2072-1439.2016.01.19&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26904255&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F24%2F2020.09.18.20196980.atom) 25. [25].Lee A, Kinney P, Chillrud S, Jack D. A systematic review of innate immunomodulatory effects of household air pollution secondary to the burning of biomass fuels, Ann Glob Health. 2015 May-Jun; 81(3): 368–374. DOI: 10.1016/j.aogh.2015.08.006 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.aogh.2015.08.006&link_type=DOI) 26. [26].Comunian S, Dongo D, Milani C, Palestini P, Air Pollution and COVID-19: The Role of particulate matter in the spread and increase of COVID-19’s morbidity and mortality, Int J Environ Res Public Health. 2020 Jun; 17(12): 4487. DOI: 10.3390/ijerph17124487 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3390/ijerph17124487&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32580440&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F24%2F2020.09.18.20196980.atom) 27. [27].Comunian S, Dongo D, Milani C, Palestini P. Air Pollution and COVID-19: The Role of particulate matter in the spread and increase of COVID-19’s morbidity and mortality, Int J Environ Res Public Health. 2020 Jun; 17(12): 4487. DOI: 10.3390/ijerph17124487 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3390/ijerph17124487&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32580440&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F24%2F2020.09.18.20196980.atom) 28. [28].Zhang L, Jackson CB, Mou H, et al. The D614G mutation in the SARS-CoV-2 spike protein reduces S1 shedding and increases infectivity, Version 1. bioRxiv. Preprint. 2020 Jun 12. DOI: 10.1101/2020.06.12.148726 [PREPRINT] [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiYmlvcnhpdiI7czo1OiJyZXNpZCI7czoxOToiMjAyMC4wNi4xMi4xNDg3MjZ2MSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzA5LzI0LzIwMjAuMDkuMTguMjAxOTY5ODAuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) [1]: /embed/inline-graphic-1.gif [2]: /embed/inline-graphic-2.gif