- Open Access
An Italian individual-level data study investigating on the association between air pollution exposure and Covid-19 severity in primary-care setting
BMC Public Health volume 21, Article number: 902 (2021)
Several studies have been focusing on the potential role of atmospheric pollutants in the diffusion and impact on health of Covid-19. This study’s objective was to estimate the association between ≤10 μm diameter particulate matter (PM10) exposure and the likelihood of experiencing pneumonia due to Covid-19 using individual-level data in Italy.
Information on Covid-19 patients was retrieved from the Italian IQVIA® Longitudinal Patient Database (LPD), a computerized network of general practitioners (GPs) including anonymous data on patients’ consultations and treatments. All patients with a Covid-19 diagnosis during March 18th, 2020 – June 30th, 2020 were included in the study. The date of first Covid-19 registration was the starting point of the 3-month follow-up (Index Date). Patients were classified based on Covid-19-related pneumonia registrations on the Index date and/or during follow-up presence/absence. Each patient was assigned individual exposure by calculating average PM10 during the 30-day period preceding the Index Date, and according to GP’s office province. A multiple generalized linear mixed model, mixed-effects logistic regression, was used to assess the association between PM10 exposure tertiles and the likelihood of experiencing pneumonia.
Among 6483 Covid-19 patients included, 1079 (16.6%) had a diagnosis of pneumonia. Pneumonia patients were older, more frequently men, more health-impaired, and had a higher individual-level exposure to PM10 during the month preceding Covid-19 diagnosis. The mixed-effects model showed that patients whose PM10 exposure level fell in the second tertile had a 30% higher likelihood of having pneumonia than that of first tertile patients, and the risk for those who were in the third tertile was almost doubled.
The consistent findings toward a positive association between PM10 levels and the likelihood of experiencing pneumonia due to Covid-19 make the implementation of new strategies to reduce air pollution more and more urgent.
Covid-19 infection, whose pathogenic agent is Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2), was firstly reported in Wuhan, Hubei Province, China, in December 2019 . The spectrum of Covid-19 clinical manifestation is very wide, as it includes asymptomatic infection, mild upper respiratory tract symptoms, mild and severe pneumonia, acute respiratory distress syndrome, sepsis and septic shock,  but also myalgia, fatigue, sputum production, headache, hemoptysis, diarrhea, anosmia, and ageusia. Furthermore, as clinical data became available, potential cardiovascular, gastrointestinal, neurological, and cutaneous manifestations of the disease were described by scientific literature . Despite the course of the disease is often mild and undistinguishable from a common flu, in a considerable number of cases Covid-19 infection may require hospitalization, and can lead to an acute respiratory distress syndrome (ARDS) and death  . Covid-19 outbreak was officially confirmed as a pandemic on 11 February 2020,  and since the end of February 2020, contagious has rapidly spread in Italy, particularly in the North (Lombardy, Veneto and Emilia-Romagna), and in many other European countries. As of the end of March, Italy represented the third country worldwide in terms of total number of cases and the first one in terms of total number of deaths . As of March 2021, it is estimated that more than 3 million people have contracted Covid-19, and more than 100,000 persons died as a consequence of this infection in Italy. Regions where most of the new cases and deaths are reported are still those of Northern Italy, and particularly Lombardy . A large part of Northern Italy territory includes the Po basin area (Padan Plain), which is the most industrialized area of the country, and one of the most polluted region of Europe . The World Health Organization (WHO), reported that air pollution is responsible for 7 million deaths worldwide every year, and represented one of the main concern regarding public health . Airborne particulate matter (PM), which is a heterogeneous mixture of solid and liquid, organic and inorganic material suspended in air, is considered as the most relevant component of air pollution .
A growing body of evidences is showing that some countries were subject to a greater spread of Covid-19 and suffered higher lethality, and this undoubtably captured researchers attention . As a consequence, several studies have been focusing on the potential role of atmospheric pollutants, particularly PM, in the diffusion of Covid-19 both in the short- and the long-term, as well as in the impact of the virus on human health [5, 11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46]. One of the ideas underneath the potential relationship between airborne PM and Covid-19 diffusion is that the atmospheric PM might exercise a carrier action along with the virus [5, 26]. Setti and colleagues recently demonstrated the presence of the SARS-CoV-2 RNA on PM . Furthermore, it has been hypothesised that the presence of air-related pollutants can put pressure on the health conditions of the populations at risk, thus offering preconditions for the development of Covid-19 and its complications, including life-threatening ones . Almost all of the studies so far performed succeeded in finding a positive association between air pollution and Covid-19 diffusion and its health-related outcomes. However, it is worth mentioning that most of them were ecological studies, thus relied on aggregated data only. Differently, the analysis here reported used patient-level data. In particular, data available from the Italian Regional Environmental Protection Agencies (ARPA) on ≤10 μm diameter PM (PM10) daily concentrations were collected to assign an individual exposure to a sample of patients experiencing Covid-19 infection. Patients’ data on Covid-19 was retrieved from IQVIA® Longitudinal Patient Database (LPD), a large repository of secondary data fed by approximately 900 Italian General Practitioners (GPs).
To authors’ knowledge, this was the first study ever using patient-level data aimed at estimating the association between short-term exposure to PM10 and the likelihood of experiencing pneumonia due to Covid-19 infection as a proxy of disease severity on a large sample of patients in Italy.
Information on Covid-19 patients was retrieved from the Italian IQVIA® LPD database. Italian IQVIA® LPD is part of a computerized network of GPs from different European countries feeding a centralized database with extensive and anonymous data on patients’ consultations and treatments. This database reflects the clinical practice of a national sample of GPs since it allows the collection and longitudinal analysis of data taken from patients’ records related to prescriptions and healthcare resource utilisations in everyday clinical practice. Drug prescriptions and medical diagnoses are both coded directly by GPs and comply with the Anatomical Therapeutic and Chemical (ATC) classification system, and with the 9th edition of International Classification of Disease (ICD-9-CM), respectively. Currently, about 900 Italian GPs contribute to the IQVIA® LPD, providing data from routinely collected records of ~ 1.2 million patients. The Italian IQVIA® LPD database, established in 1998 by the Italian College of General Practitioners (Società Italiana di Medicina Generale - SIMG), was found to be representative of the Italian general population [48,49,50] and a reliable source of information in numerous previous studies for several disease areas [51,52,53,54,55,56,57].
PM10 daily concentration data detected by the official air quality monitoring stations located on the entire Italian territory were retrieved from ARPA Regional websites for the period January 2020 – June 2020 and are publicly available.
Study population, exposure, and outcomes definitions
On March 17th, 2020, all the GPs collaborating with the Italian IQVIA® LPD were asked by SIMG/HealthSearch to use two different codes to distinguish between Covid-19 patients with or without pneumonia. The present analysis firstly selected all patients who had at least one Covid-19 registration during the period March 18th, 2020 – June 30th, 2020 (selection period). The date of the first Covid-19 registration was considered as the Index Date and patients were grouped based on the presence/absence of a registration of Covid-19 with pneumonia within 90 days since the Index Date (follow-up) as a proxy of disease severity. The final study cohort was composed of Covid-19 patients who had data availability for the study period and no missing information on age, sex, and PM10 exposure. For each patient included in the analysis, diagnoses recorded during the 12-month period preceding the Index Date (baseline) were collected, and the following conditions, defined through ICD-9 codes, were considered as comorbidities of interest: hypertension (ICD-9 code 401.xx), diabetes (ICD-9 code 250.xx), asthma (ICD-9 code 493.xx), chronic obstructive pulmonary disease (COPD) (ICD-9 codes 490.xx, 491.xx, 492.xx, 494.xx, 495.xx, 496.xx), obesity (ICD-9 codes 278.00, 278.01 and/or a body mass index value of at least 30 kg/m2), coronary artery disease (CAD) (ICD-9 codes 410.xx, 411.xx, 412.xx, 413.xx, 414.xx, V45.81, V45.82), and cerebrovascular disease (ICD-9 codes from 430.xx to 438.xx). Information on age and sex on the Index Date, smoking habits, and GP’s office province were also retrieved from the IQVIA® LPD database. The daily average value of PM10 (μg/m3) was calculated at province level based on data from ARPA air quality monitoring stations. Individual PM10 exposure was then assigned to each patient by calculating the mean of PM10 daily values observed during the month preceding the Index Date and according to GP’s office province.
A description of patients’ individual PM10 exposure was given in terms of mean and standard deviation by Region. Descriptive statistics on patients’ characteristics were reported for the two groups of patients with or without pneumonia. A stratification of patients by tertiles of PM10 calculated on the overall cohort and presence or absence of pneumonia was also provided. Because the two groups of Covid-19 patients with or without pneumonia were independent, when dealing with categorical variables, Chi-square tests were performed to investigate on between-groups differences. Because of the large sample size, Kolmogorov-Smirnov test was performed to test normality of PM10 exposure assessed as a continuous variable. P-value resulting from the test was very close to 0 (i.e., < 0.01), thus a non-parametric Wilcoxon signed-rank test was performed to understand whether PM10 exposure differed between the two groups. P-values resulting from Chi-square and Wilcoxon tests were reported as continuous variables. As p-values can range from 0 to 1, the closer the values were to zero, the higher the diffidence in the null hypothesis of no association . Given the multilevel structure of the data, a multiple mixed-effects logistic regression model was built by inserting a random intercept for GPs’ office province to assess the association between PM10 exposure and the likelihood of experiencing pneumonia as a consequence of Covid-19 infection. Indeed, no further details in addition to province were available regarding GPs’ offices location, and then on patients’ residence. Being so, PM10 values were used only at province level, with this possibly introducing within-province correlation despite an individual level of exposure was assigned to each patient depending on his/her Index Date. Furthermore, due to the pandemic, the disease status of subjects in the same province may be correlated, violating the independence assumption of the ordinary regression model . Fixed effects included in the model were PM10 exposure in terms of tertiles, sex, age class, and presence or absence of each comorbidity of interest during baseline period. Given the large number of patients involved, we used the maximum likelihood as an estimation method. Covariance structure was assumed to be variance component and the ratio between the Pearson chi-square statistic and its degrees of freedom was calculated to verify that variability has been properly modeled. Associations were expressed as Odds Ratios (ORs), together with their 95% confidence intervals (CI). Point estimates should be regarded as the most compatible values, while the other values included in the CIs should be regarded as progressively less compatible (but nevertheless still compatible) the greater their distance from the point estimate . Multiple logistic mixed models have been previously used in literature to deal with data having an analogous structure to investigate on the association between air pollution and health outcomes [59,60,61]. Sensitivity analyses were performed to assess the robustness of results. Firstly, we run the above described model on the subgroup of subjects with smoking habits information available to account for the potential confounding effect of cigarettes smoking on the association between exposure to PM10 and risk of pneumonia. A dummy variable indicating whether a patient have ever smoked or not was included among model’s covariates. Secondly, we decided to run the multiple mixed-effects logistic regression model using PM10 tertiles calculated based on patients’ individual exposure during the 14- and 7-day periods preceding the Index Date as exposure variables. The latter analyses were intended to test the robustness of results with respect to the time-window chosen to assess PM10 exposure. Thirdly, stratified analyses were performed by running the multiple mixed-effects logistic regression model separately on men and women, on < 65 and ≥ 65 years old subjects, and on never smoking patients. C-statistic was calculated to assess models’ discrimination performance (i.e., the extent to which patients who were predicted to be high risk exhibited higher pneumonia rates compared with those who were predicted to be low risk); C-statistic can range from 0.5 (poor discrimination) to 1.0 (perfect discrimination) and was used in previous studies as a measure of model’s discrimination ability . All the analyses were performed using SAS software, version 9.4.
A final cohort of 6483 patients was defined according to eligibility criteria. Overall, Covid-19 patients without pneumonia numbered 5404 (83.4%), while those with pneumonia were 1079 (16.6%). Figure 1, which describes individual PM10 exposure by Region, shows that on average, Veneto, Lombardy, Marche, and Campania were the Regions where patients have been exposed to the highest level of PM10 during the month preceding the Index Date. In particular, PM10 mean values ranged from 25.3 to 28.2 μg/m3, thus being higher than the mean value calculated on the total cohort (24.1 μg/m3). For all the remaining Regions, PM10 mean values fell below the Italian mean value and ranged from 14.0 μg/m3 for Valdaosta to 23.1 μg/m3 for Emilia-Romagna (Fig. 1). Molise Region was not included in the analysis due to IQVIA® LPD data unavailability.
Table 1 reports Covid-19 patients’ characteristics stratified by pneumonia presence or absence. Among Covid-19 patients with pneumonia we found a higher proportion of male than the one observed among patients without pneumonia (57.1% versus 45.9%). Overall, mean age of Covid-19 patients was around 55 years, and subjects with pneumonia were older than those without pneumonia, with mean age being around 64 and 53 years respectively for the two groups (data not shown). Consistently, the proportion of Covid-19 patients aged 65 years or older was 50.4 and 26.9% for the group with and without pneumonia, respectively. Ex- and current smokers accounted for a higher proportion among Covid-19 patients without pneumonia. Overall, the most frequently reported comorbidities among those of interest were hypertension, obesity, and diabetes mellitus. The comparison between groups showed that among pneumonia patients the proportions of subjects with each condition were higher than those observed among patients without pneumonia. PM10 exposure during the month preceding the Index Date measured as a continuous variable was higher for pneumonia patients (Table 1).
The stratification of Covid-19 patients by presence/absence of pneumonia and tertiles of PM10 exposure showed that a much higher proportion of subjects whose mean level of exposure fell in the third tertile was found for pneumonia group. Indeed, almost half of patients with pneumonia had a PM10 mean value higher than 28.7 μg/m3 during the month preceding the Index Date, i.e., Covid-19 diagnosis (Fig. 2).
Results from the multiple mixed-effect logistic regression model showed that PM10 exposure level during the month preceding the Index Date was associated with the risk of experiencing pneumonia due to Covid-19 infection. In particular, patients whose level of PM10 exposure fell in the second tertile had a 30% higher likelihood of having pneumonia than patients whose level of exposure fell into the first tertile; the risk for patients in the third PM10 exposure tertile was almost doubled. Among the covariates included in the model, gender, age class, and presence/absence of baseline asthma and or obesity had an influence on the likelihood of experiencing pneumonia. In particular, male sex, older age, as well as the presence of asthma and/or obesity exposed patients to a higher risk of pneumonia as a consequence of Covid-19 infection. The ratio of the Pearson chi-square statistic and its degrees of freedom was 0.96. The C-statistic was 0.8, thus discrimination performance of the model was good (Table 2).
Results from the sensitivity analysis run on the subgroup of patients with available information on smoking habits confirmed those from the model run on the total cohort and did not show any association between smoking habits and risk of pneumonia (data not shown). Furthermore, results from the multiple mixed-effect logistic regression models which accounted for PM10 exposure assessed during the 14- and 7-day periods preceding the Index Date (Tables 1S and 2S), as well as those from stratified analyses (Tables 3S, 4S, 5S, 6S, 7S) confirmed the robustness of our findings.
The main objective of the present study was to understand whether short-term exposure to PM10 may increase Covid-19 patients’ likelihood of experiencing pneumonia as a proxy of disease severity. Findings from the multiple mixed-effect logistic model suggest that short-term exposure to PM10 may represent a risk factor for the development of pneumonia in patients with Covid-19 infection. Furthermore, an increasing trend in the likelihood of experiencing pneumonia was observed corresponding to increasing levels of PM10.
Several studies previously investigated on the association between air pollution and Covid-19 spread and adverse outcomes in Italy [4, 8, 10, 11, 13,14,15, 29, 30, 35, 36, 40, 44, 46]. A study by Conticini and colleagues focused on two Northern Italy Regions, Lombardy and Emilia Romagna, which are part of the Padan plain. According to data from Italian Civil Protection, these Regions had the highest level of virus lethality in the world at the time of the first epidemic wave . Being so, Conticini and colleagues speculated that the high level of pollution should be considered as an additional co-factor of the high level of lethality recorded in that area . Dettori and colleagues examined the role of air pollutants in relation to the number of deaths per each Italian province affected by Covid-19. PM10 was found to be an independent predictor for Covid-19-related mortality . Similarly, Bianconi et al. investigated on the association between PM exposure and Covid-19 cases and related deaths both at Regional and province level in Italy. Study results seemed to suggest that the greater diffusion and lethality of Covid-19 might be at least partially related to the past and cumulative PM exposure . Accarino and colleagues, who investigated on short-term exposure to atmospheric pollutants and spatio-temporal distribution of Covid-19 cases and deaths in Italy, also suggested a potential correlation, particularly with PMs . Coker and colleagues found that a 1 μg/m3 increase in long-term exposure to PM2.5 was associated with a 9% increase in COVID-19 related mortality in Northern Italy . Similar studies have been conducted also in other countries [17, 20, 23, 25, 27, 30, 33, 34, 36, 38, 39, 41,42,43, 45]. A study by Wu and colleagues investigating on the association between PM2.5 exposure and risk of COVID-19 death in the United States found that an increase of 1 μg/m3 in PM2.5 was associated with an 8% increase in the COVID-19 death rate . Cole and colleagues found a positive relationship between air pollution and Covid-19 cases, hospital admissions and deaths using data from 355 Dutch municipalities .
Findings from all the above-mentioned analyses came from ecological studies, which used aggregated data. Undoubtably, ecological studies, whose approach is extremely cost-effective, are crucial in rapidly evolving areas of research. Indeed, they allow drawing area level conclusions, which can be useful for policy-making . However, ecological regression analyses are unable to adjust for individual-level risk factors, which, instead, are known to affect Covid-19 health outcomes. To authors’ knowledge, very few studies have been performed using individual-level data, and none of them was conducted in Italy. Travaglio and colleagues investigated on the associations between several air pollutants and the risk of COVID-19 infection using patient-level data obtained from a cohort of 1450 subjects in the UK. Results from the analysis showed that levels of PM pollutants and nitrogen oxides were associated with an increase in SARS-CoV-2 infections. No investigations were done accounting for Covid-19 severity or mortality . Another study conducted in the UK and using individual-level data found a positive association between exposure to NO2 and Covid-19 mortality, while the association with PM2.5 was uncertain . Finally, a study conducted in Mexico City used patient-level data to estimate the effects of both long- and short-term exposure to PM2.5 on Covid-19 mortality: evidences toward a positive relationship between PM2.5 air pollution and the likelihood for an individual to die following Covid-19 infection did emerge; this relationship increased with age, and, although findings suggested that the association was mainly driven by long-term exposure, authors did not exclude that short-term exposure might also have an effect .
In light of the few studies using patient-level data, authors of the present study do believe that findings here reported should be regarded as extremely significant and add an important contribution in the understanding of the relationship between air pollution and Covid-19 severity. Furthermore, it is worth mentioning that individual-level data in the present study included information on exposure in addition to that on potential confounders. In fact, PM10 exposure was calculated for each patient considering a specific time-window defined based on Covid-19 registration date. Moreover, the inclusion of data from the entire Italian territory and the extension of the study period up to the end of June, allowed to account for a very high variability in terms of exposure. Besides this, important strengths of the present study are the representativeness of the data source used to identify Covid-19 cases , as well as the setting where data was collected. Indeed, GPs are on the front-line in the management of this pandemic as they are the first point of contact for people affected by Covid-19, except for those patients who develop extremely severe forms of the disease since its onset. As such, we were able to retrieve information on a representative sample of patients in terms of disease severity, and this prevented our study from the risk of selection bias which may have occurred using inpatient-setting data. Risk factors for Covid-19-related pneumonia other than PM10 exposure identified by the present study were older age, male sex, asthma, and obesity. These findings agree with those from previous studies, with this further confirming the robustness of our data. In particular, Polverino et al. found that 65 years older age and male sex were among predictors of death in a sample of Covid-19 inpatients . Similarly, increasing age was one of the independent risk factors for all-cause in-hospital mortality in a study conducted on 317 hospitalized adult patients with a diagnosis of Covid-19 . Baronio and colleagues found that admission to intensive care unit (ICU) and poor survival were associated with advanced age and higher body mass index . Furthermore, obesity was found to be a strong, independent risk factor for respiratory failure, admission to the ICU and death in a sample of 482 Covid-19 hospitalized patients . Finally, a study conducted on behalf of the National Health System (NHS) in England linked primary-care electronic medical records of 17,278,392 adults to 10,926 Covid-19-related deaths; among factors associated with Covid-19 death there were male sex, greater age, and severe asthma .
The present study also presents some limits. Firstly, we do not know whether Covid-19 diagnoses were confirmed by a nasopharyngeal swab. However, demographic characteristics of the overall cohort are in line with Covid-19 cases description provided by the Italian Istituto Superiore di Sanità (ISS). Indeed, according to ISS data, the proportion of men among subjects affected by Covid-19 as of the end of June 2020 was 45.8%, compared to the 47.7% observed in the present study. Also, age class distributions were in line, even if patients in the present study were just slightly younger . However, it should be considered that the selection period for the present study was delayed with respect to Covid-19 outbreak, and, differently from sex distribution which remained constant, mean age of Covid-19 patients has progressively decreased . In light of the above consideration, we do believe that IQVIA® LPD is a reliable data source for the identification of Covid-19 cases. Secondly, the main analysis investigating on the association between PM10 exposure and risk of Covid-19-related pneumonia performed on the overall cohort did not account for smoking habits due to the limited availability of such information. However, results from the sensitivity analysis run on the subgroup of patients confirmed all the associations found by the analysis performed on the overall cohort. Furthermore, it is worth mentioning that there were previous studies that did not find a correlation between smoking and adverse outcomes in Covid-19 patients . Thirdly, it is possible that the date of the Covid-19 diagnosis registration did not exactly correspond to the date of infection’s onset. However, authors of the present study are confident that the application of a 30-day time-window to estimate average PM10 exposure should have mitigated the potential effect of this limitation. Finally, PM10 comprises PM2.5 and PM5, and may also contain a heterogeneous group of gaseous compounds such as sulfur dioxide (SO2) and nitrogen dioxide (NO2) . As such, it cannot be excluded that the adverse effect of PM10 on Covid-19 outcomes we observed might be attributed to one or more of PM10 components, instead of PM10 itself. This is the reason why authors of the present study would like to claim the importance of reducing air pollution as a whole.
This is the first study aimed to investigate on the association between PM10 exposure and the risk of developing Covid-19-related pneumonia as a proxy of disease severity using individual-level data in Italy. Consistent findings toward a positive association between PM10 levels and the likelihood of experiencing pneumonia were found. Authors of the present study would like to claim the urgency of implementing new strategies to reduce air pollution.
Availability of data and materials
IQVIA® LPD data that support the findings of this study are available from IQVIA, but restrictions apply to the availability of these data, which were used under license for the current study, and so are not publicly available. Data are, however, available from the corresponding author upon reasonable request and with permission of IQVIA. PM10 daily concentration data detected by the official air quality monitoring stations located on the entire Italian territory were retrieved from ARPA Regional websites for the period January 2020 – June 2020 and are publicly available.
Riou J, Althaus CL. Pattern of early human-to-human transmission of Wuhan 2019 novel coronavirus (2019-nCoV), December 2019 to January 2020. Eurosurveillance. 2020;25(4) [cited 2021 Apr 1]. Available from: https://www.eurosurveillance.org/content/10.2807/1560-7917.ES.2020.25.4.2000058.
Menzella F, Matucci A, Vultaggio A, Barbieri C, Biava M, Scelfo C, et al. COVID-19: general overview, pharmacological options and ventilatory support strategies. Multidiscip Respir Med. 2020 Jan 28;15(1):708. https://doi.org/10.4081/mrm.2020.708.
Mizrahi B, Shilo S, Rossman H, Kalkstein N, Marcus K, Barer Y, et al. Longitudinal symptom dynamics of COVID-19 infection. Nat Commun. 2020 Dec;11(1):6208. https://doi.org/10.1038/s41467-020-20053-y.
Conticini E, Frediani B, Caro D. Can atmospheric pollution be considered a co-factor in extremely high level of SARS-CoV-2 lethality in Northern Italy? Environ Pollut Barking Essex. 2020;261:114465.
Martelletti L, Martelletti P. Air pollution and the novel Covid-19 disease: a putative disease risk factor. SN Compr Clin Med. 2020;15:1–5.
Ortenzi F, Albanese E, Fadda M. A Transdisciplinary analysis of COVID-19 in Italy: the Most affected country in Europe. Int J Environ Res Public Health. 2020 Dec 18;17(24):9488. https://doi.org/10.3390/ijerph17249488.
Protezione Civile Italiana. Map available at: https://opendatadpc.maps.arcgis.com/apps/opsdashboard/index.html#/b0c68bce2cce478eaac82fe38d4138b1. Last accessed on 31 March 2021.
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 Mar 15;760:143355. https://doi.org/10.1016/j.scitotenv.2020.143355.
World Health Organization. Air pollution. Available at: https://www.who.int/health -topics/air-pollution. Last accessed on 14 Feb 2021.
Bianconi V, Bronzo P, Banach M, Sahebkar A, Mannarino MR, Pirro M. Particulate matter pollution and the COVID-19 outbreak: results from Italian regions and provinces. Arch Med Sci AMS. 2020;16(5):985–92. https://doi.org/10.5114/aoms.2020.95336.
Dettori M, Deiana G, Balletto G, Borruso G, Murgante B, Arghittu A, et al. Air pollutants and risk of death due to COVID-19 in Italy. Environ Res. 2021 Jan;192:110459. https://doi.org/10.1016/j.envres.2020.110459.
Li H, Xu X-L, Dai D-W, Huang Z-Y, Ma Z, Guan Y-J. Air pollution and temperature are associated with increased COVID-19 incidence: a time series study. Int J Infect Dis. 2020 Aug;97:278–82. https://doi.org/10.1016/j.ijid.2020.05.076.
Bontempi E. First data analysis about possible COVID-19 virus airborne diffusion due to air particulate matter (PM): the case of Lombardy (Italy). Environ Res. 2020 Jul;186:109639. https://doi.org/10.1016/j.envres.2020.109639.
Coccia M. Factors determining the diffusion of COVID-19 and suggested strategy to prevent future accelerated viral infectivity similar to COVID. Sci Total Environ. 2020 Aug 10;729:138474. https://doi.org/10.1016/j.scitotenv.2020.138474.
Fattorini D, Regoli F. Role of the chronic air pollution levels in the Covid-19 outbreak risk in Italy. Environ Pollut Barking Essex. 2020;264:114732.
Frontera A, Cianfanelli L, Vlachos K, Landoni G, Cremona G. Severe air pollution links to higher mortality in COVID-19 patients: the “double-hit” hypothesis. J Inf Secur. 2020;81(2):255–9.
Liang D, Shi L, Zhao J, Liu P, Schwartz J, Gao S, Sarnat J, Liu Y, Ebelt S, Scovronick N, Chang HH. Urban Air Pollution May Enhance COVID-19 Case-Fatality and Mortality Rates in the United States. medRxiv [Preprint]. 2020:2020.05.04.20090746. https://doi.org/10.1101/2020.05.04.20090746. Update in: Innovation (N Y). 2020 Sep 21;:100047.
Onder G, Rezza G, Brusaferro S. Case-Fatality Rate and Characteristics of Patients Dying in Relation to COVID-19 in Italy. JAMA. 2020; [cited 2021 Jan 19]; Available from: https://jamanetwork.com/journals/jama/fullarticle/2763667.
Riccò M, Ranzieri S, Balzarini F, Bragazzi NL, Corradi M. SARS-CoV-2 infection and air pollutants: correlation or causation? Sci Total Environ. 2020 Sep 10;734:139489. https://doi.org/10.1016/j.scitotenv.2020.139489.
Wu X, Nethery RC, Sabath MB, Braun D, Dominici F. Exposure to air pollution and COVID-19 mortality in the United States: A nationwide cross-sectional study. Epidemiology. 2020; [cited 2021 Jan 21]. Available from: http://medrxiv.org/lookup/doi/10.1101/2020.04.05.20054502.
Zoran MA, Savastru RS, Savastru DM, Tautan MN. Assessing the relationship between surface levels of PM2.5 and PM10 particulate matter impact on COVID-19 in Milan, Italy. Sci Total Environ. 2020 Oct 10;738:139825. https://doi.org/10.1016/j.scitotenv.2020.139825.
Zoran MA, Savastru RS, Savastru DM, Tautan MN. Assessing the relationship between ground levels of ozone (O3) and nitrogen dioxide (NO2) with coronavirus (COVID-19) in Milan, Italy. Sci Total Environ. 2020;740:140005. https://doi.org/10.1016/j.scitotenv.2020.140005.
Cole MA, Ozgen C, Strobl E. Air pollution exposure and Covid-19 in Dutch municipalities. Environ Resour Econ. 2020 Aug;76(4):581–610. https://doi.org/10.1007/s10640-020-00491-4.
Magazzino C, Mele M, Schneider N. The relationship between air pollution and COVID-19-related deaths: an application to three French cities. Appl Energy. 2020 Dec 1;279:115835. https://doi.org/10.1016/j.apenergy.2020.115835.
Zhu Y, Xie J, Huang F, Cao L. Association between short-term exposure to air pollution and COVID-19 infection: evidence from China. Sci Total Environ. 2020 Jul 20;727:138704. https://doi.org/10.1016/j.scitotenv.2020.138704.
Setti L, Passarini F, de Gennaro G, Di Gilio A, Palmisani J, Buono P, et al. SIMA. Società Italiana di Medicina Ambientale. Position Paper. Relazione circa l’effetto dell’inquinamento da particolato atmosferico e la diffusione di virus nella popolazione. Available at: https://www.simaonlus.it/wpsima/wp-content/uploads/2020/03/COVID19_Position-Paper_Relazione-circa-l%E2%80%99effetto-dell%E2%80%99inquinamento-da-particolato-atmosferico-e-la-diffusione-di-virus-nella-popolazione.pdf. Last accessed on 14 Feb 2021.
Adhikari A, Yin J. Short-Term Effects of Ambient Ozone, PM2.5, and Meteorological Factors on COVID-19 Confirmed Cases and Deaths in Queens, New York. Int J Environ Res Public Health. 2020;17(11):4047. https://doi.org/10.3390/ijerph17114047.
Andrée BPJ. Incidence of COVID-19 and Connections with Air Pollution Exposure: Evidence from the Netherlands. Epidemiology; 2020 May [cited 2021 Jan 20]. Available from: http://medrxiv.org/lookup/doi/10.1101/2020.04.27.20081562
Coker ES, Cavalli L, Fabrizi E, Guastella G, Lippo E, Parisi ML, et al. The effects of air pollution on COVID-19 related mortality in northern Italy. Environ Resour Econ. 2020;4:1–24.
Pansini R, Fornacca D. Early evidence of a higher incidence of COVID-19 in the air-polluted regions of eight severely affected countries. Epidemiology; 2020 [cited 2021 Jan 20]. Available from: http://medrxiv.org/lookup/doi/10.1101/2020.04.30.20086496
Konstantinoudis G, Padellini T, Bennett J, Davies B, Ezzati M, Blangiardo M. Long-term exposure to air-pollution and COVID-19 mortality in England: a hierarchical spatial analysis. Environ Int. 2021 Jan;146:106316. https://doi.org/10.1016/j.envint.2020.106316.
Travaglio M, Yu Y, Popovic R, Selley L, Leal NS, Martins LM. Links between air pollution and COVID-19 in England. Environ Pollut. 2021 Jan;268(Pt A):115859. https://doi.org/10.1016/j.envpol.2020.115859.
Vasquez-Apestegui V, Parras-Garrido E, Tapia V, Paz-Aparicio VM, Rojas JP, Sánchez-Ccoyllo OR, Gonzales GF. Association Between Air Pollution in Lima and the High Incidence of COVID-19: Findings from a Post Hoc Analysis. Res Sq [Preprint]. 2020:rs.3.rs-39404. https://doi.org/10.21203/rs.3.rs-39404/v1.
Yao Y, Pan J, Liu Z, Meng X, Wang W, Kan H, et al. Temporal association between particulate matter pollution and case fatality rate of COVID-19 in Wuhan. Environ Res. 2020 Oct;189:109941. https://doi.org/10.1016/j.envres.2020.109941.
Accarino G, Lorenzetti S, Aloisio G. Assessing correlations between short-term exposure to atmospheric pollutants and COVID-19 spread in all Italian territorial areas. Environ Pollut. 2021 Jan;268(Pt A):115714. https://doi.org/10.1016/j.envpol.2020.115714.
Ogen Y. Assessing nitrogen dioxide (NO2) levels as a contributing factor to coronavirus (COVID-19) fatality. Sci Total Environ. 2020 Jul 15;726:138605. https://doi.org/10.1016/j.scitotenv.2020.138605.
López-Feldman A, Heres D, Marquez-Padilla F. Air pollution exposure and COVID-19: a look at mortality in Mexico City using individual-level data. Sci Total Environ. 2021 Feb 20;756:143929. https://doi.org/10.1016/j.scitotenv.2020.143929.
Hou CK, Qin YF, Wang G, Liu QL, Yang XY, Wang H. Impact of a long-term air pollution exposure on the case fatality rate of COVID-19 patients-A multicity study. J Med Virol. 2021;93(5):2938–46. https://doi.org/10.1002/jmv.26807. Epub 2021 Feb 15.
Bilal null, Bashir MF, Komal B, Benghoul M, Bashir MA, Tan D. Nexus Between the COVID-19 Dynamics and Environmental Pollution Indicators in South America. Risk Manag Healthc Policy. 2021;14:67–74.
Cascetta E, Henke I, Di Francesco L. The Effects of Air Pollution, Sea Exposure and Altitude on COVID-19 Hospitalization Rates in Italy. Int J Environ Res Public Health. 2021;18(2):452. https://doi.org/10.3390/ijerph18020452.
Fernández D, Giné-Vázquez I, Liu I, Yucel R, Nai Ruscone M, Morena M, et al. Are environmental pollution and biodiversity levels associated to the spread and mortality of COVID-19? A four-month global analysis. Environ Pollut Barking Essex. 2020;271:116326.
Hassan MS, Bhuiyan MAH, Tareq F, Bodrud-Doza M, Tanu SM, Rabbani KA. Relationship between COVID-19 infection rates and air pollution, geo-meteorological, and social parameters. Environ Monit Assess. 2021 Jan 4;193(1):29. https://doi.org/10.1007/s10661-020-08810-4.
Hutter HP, Poteser M, Moshammer H, Lemmerer K, Mayer M, Weitensfelder L, Wallner P, Kundi M. Air Pollution Is Associated with COVID-19 Incidence and Mortality in Vienna, Austria. Int J Environ Res Public Health. 2020;17(24):9275. https://doi.org/10.3390/ijerph17249275.
Fiasca F, Minelli M, Maio D, Minelli M, Vergallo I, Necozione S, et al. Associations between COVID-19 incidence rates and the exposure to PM2.5 and NO2: a Nationwide observational study in Italy. Int J Environ Res Public Health. 2020;13:17(24).
Hoang T, Nguyen TQ, Tran TTA. Short-term exposure to ambient air pollution in association with COVID-19 of two clusters in South Korea. Trop Med IntHealth. 2021;26(4):478–91.
Rovetta A, Castaldo L. Relationships between demographic, geographic, and environmental statistics and the spread of novel coronavirus disease (COVID-19) in Italy. Cureus. 2020 Nov 9;12(11):e11397. https://doi.org/10.7759/cureus.11397.
Setti L, Passarini F, De Gennaro G, Barbieri P, Perrone MG, Borelli M, et al. SARS-Cov-2RNA found on particulate matter of Bergamo in northern Italy: first evidence. Environ Res. 2020;188:109754. https://doi.org/10.1016/j.envres.2020.109754.
Filippi A, Vanuzzo D, Bignamini AA, Sessa E, Brignoli O, Mazzaglia G. Computerized general practice databases provide quick and cost-effective information on the prevalence of angina pectoris. Ital Heart J Off J Ital Fed Cardiol. 2005;6(1):49–51.
XIII Report Health Search (edizione 2020). Istituto di Ricerca della SIMG. Società Italiana di Medicina Generale e delle Cure Primarie [Internet] Available from: https://www.healthsearch.it/report/
Perera G, Rijnbeek PR, Alexander M, Ansell D, Avillach P, Duarte-Salles T, et al. Vascular and metabolic risk factor differences prior to dementia diagnosis: a multidatabase case–control study using European electronic health records. BMJ Open. 2020 Nov;10(11):e038753. https://doi.org/10.1136/bmjopen-2020-038753.
Ravera M, Cannavò R, Noberasco G, Guasconi A, Cabib U, Pieracci L, et al. High performance of a risk calculator that includes renal function in predicting mortality of hypertensive patients in clinical application. J Hypertens. 2014 Jun;32(6):1245–54. https://doi.org/10.1097/HJH.0000000000000177.
Levi M, Rosselli M, Simonetti M, Brignoli O, Cancian M, Masotti A, et al. Epidemiology of iron deficiency anaemia in four European countries: a population-based study in primary care. Eur J Haematol. 2016 Dec;97(6):583–93. https://doi.org/10.1111/ejh.12776.
Katz P, Pegoraro V, Liedgens H. Characteristics, resource utilization and safety profile of patients prescribed with neuropathic pain treatments: a real-world evidence study on general practices in Europe - the role of the lidocaine 5% medicated plaster. Curr Med Res Opin. 2017 Aug;33(8):1481–9. https://doi.org/10.1080/03007995.2017.1335191.
Di Marco F, Santus P, Terraneo S, Peruzzi E, Muscianisi E, Ripellino C, et al. Characteristics of newly diagnosed COPD patients treated with triple inhaled therapy by general practitioners: a real world Italian study. NPJ Prim Care Respir Med. 2017;27(1):51. https://doi.org/10.1038/s41533-017-0051-9.
Agostoni E, Barbanti P, Frediani F, Trifirò G, Burgio L, di Nola L, et al. Real-world insights on the management of migraine patients: an Italian nationwide study. Curr Med Res Opin. 2019 Sep;35(9):1545–54. https://doi.org/10.1080/03007995.2019.1602032.
Levi M, Simonetti M, Marconi E, Brignoli O, Cancian M, Masotti A, et al. Gender differences in determinants of iron-deficiency anemia: a population-based study conducted in four European countries. Ann Hematol. 2019 Jul;98(7):1573–82. https://doi.org/10.1007/s00277-019-03707-w.
Lavorini F, Bianco A, Blasi F, Braido F, Corsico AG, Di Marco F, et al. What drives inhaler prescription for asthma patients? Results from a real-life retrospective analysis. Respir Med. 2020;166:105937.
Andrade C, The P. Value and statistical significance: misunderstandings, explanations, challenges, and alternatives. Indian J Psychol Med. 2019;41(3):210–5. https://doi.org/10.4103/IJPSYM.IJPSYM_193_19.
Zhang L, Yang Y, Li Y, Qian Z, Xiao W, Wang X, et al. Short-term and long-term effects of PM2.5 on acute nasopharyngitis in 10 communities of Guangdong, China. Sci Total Environ. 2019;688:136–42. https://doi.org/10.1016/j.scitotenv.2019.05.470.
Norbäck D, Lu C, Zhang Y, Li B, Zhao Z, Huang C, et al. Sources of indoor particulate matter (PM) and outdoor air pollution in China in relation to asthma, wheeze, rhinitis and eczema among pre-school children: synergistic effects between antibiotics use and PM10 and second hand smoke. Environ Int. 2019 Apr;125:252–60. https://doi.org/10.1016/j.envint.2019.01.036.
Zhang Z, Dong B, Li S, Chen G, Yang Z, Dong Y, et al. Exposure to ambient particulate matter air pollution, blood pressure and hypertension in children and adolescents: a national cross-sectional study in China. Environ Int. 2019 Jul;128:103–8. https://doi.org/10.1016/j.envint.2019.04.036.
Nikolian VC, Kamdar NS, Regenbogen SE, Morris AM, Byrn JC, Suwanabol PA, et al. Anastomotic leak after colorectal resection: a population-based study of risk factors and hospital variation. Surgery. 2017 Jun;161(6):1619–27. https://doi.org/10.1016/j.surg.2016.12.033.
31 Marzo 2021. Aggiornamento casi Covid-19. Dati aggregati quotidiani Regioni/PPA. Ministero della Salute - Istituto Superiore di Sanità. Available at: https://github.com/pcm-dpc/COVID-19/blob/master/schede-riepilogative/regioni/dpc-covid19-ita-scheda-regioni-20210331.pdf. Last accessed on 31/03/2021.
Wu X, Nethery RC, Sabath MB, Braun D, Dominici F. Air pollution and COVID-19 mortality in the United States: Strengths and limitations of an ecological regression analysis. Sci Adv. 2020;6(45):eabd4049. https://doi.org/10.1126/sciadv.abd4049.
Polverino F, Stern DA, Ruocco G, Balestro E, Bassetti M, Candelli M, et al. Comorbidities, cardiovascular therapies, and COVID-19 mortality: a Nationwide, Italian observational study (ItaliCO). Front Cardiovasc Med. 2020 Oct 9;7:585866. https://doi.org/10.3389/fcvm.2020.585866.
Vena A, Giacobbe DR, Di Biagio A, Mikulska M, Taramasso L, De Maria A, et al. Clinical characteristics, management and in-hospital mortality of patients with coronavirus disease 2019 in Genoa. Italy Clin Microbiol Infect. 2020;26(11):1537–44. https://doi.org/10.1016/j.cmi.2020.07.049.
Baronio M, Freni-Sterrantino A, Pinelli M, Natalini G, Tonini G, Marri M, et al. Italian SARS-CoV-2 patients in intensive care: towards an identikit for subjects at risk? Eur Rev Med Pharmacol Sci. 2020 Sep;24(18):9698–704. https://doi.org/10.26355/eurrev_202009_23061.
Rottoli M, Bernante P, Belvedere A, Balsamo F, Garelli S, Giannella M, et al. How important is obesity as a risk factor for respiratory failure, intensive care admission and death in hospitalised COVID-19 patients? Results from a single Italian Centre. Eur J Endocrinol. 2020 Oct;183(4):389–97. https://doi.org/10.1530/EJE-20-0541.
Williamson EJ, Walker AJ, Bhaskaran K, Bacon S, Bates C, Morton CE, et al. Factors associated with COVID-19-related death using OpenSAFELY. Nature. 2020 Aug;584(7821):430–6. https://doi.org/10.1038/s41586-020-2521-4.
Sorveglianza integrata Covid-19: archivio. Epicentro. Istituto Superiore di Sanità. Avaliable at https://www.epicentro.iss.it/coronavirus/sars-cov-2-sorveglianza-dati-archivio.
Lippi G, Henry BM. Active smoking is not associated with severity of coronavirus disease 2019 (COVID-19). Eur J Intern Med. 2020 May;75:107–8. https://doi.org/10.1016/j.ejim.2020.03.014.
Wang K, Hao Y, Au W, Christiani DC, Xia Z-L. A systematic review and meta-analysis on short-term particulate matter exposure and chronic obstructive pulmonary disease hospitalizations in China. J Occup Environ Med. 2019 Apr;61(4):e112–24. https://doi.org/10.1097/JOM.0000000000001539.
No fundings were received for the conduction of the present study.
Ethics approval and consent to participate
Studies conducted on IQVIA® LPD data do not require an approval by an Ethical Committee. Indeed, IQVIA® LPD relies on anonymous data which are not originated by any clinical trial. Before being stored into IQVIA databases, data is anonymized through non-identifiable encryption process directly on GPs’ computer and patients’ names are not collected into IQVIA® LPD. IQVIA® LPD complies with the European Regulation 679/2016 and the ex-Legislative Decree 196/03 and subsequent modifications.
Consent for publication
All the authors are employees of IQVIA Solutions Italy S.r.l. The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1.
Table 1S – Estimates from the multiple mixed-effect logistic regression model evaluating the likelihood of experiencing pneumonia. Sensitivity analysis with PM10 exposure assessed during the 14-day period preceding the Index Date (i.e., date of first Covid-19 registration). Table 2S – Estimates from the multiple mixed-effect logistic regression model evaluating the likelihood of experiencing pneumonia. Sensitivity analysis with PM10 exposure assessed during the 7-day period preceding the Index Date (i.e., date of first Covid-19 registration). Table 3S – Estimates from the multiple mixed-effect logistic regression model evaluating the likelihood of experiencing pneumonia in men. Table 4S – Estimates from the multiple mixed-effect logistic regression model evaluating the likelihood of experiencing pneumonia in women. Table 5S – Estimates from the multiple mixed-effect logistic regression model evaluating the likelihood of experiencing pneumonia in 65 years younger subjects. Table 6S – Estimates from the multiple mixed-effect logistic regression model evaluating the likelihood of experiencing pneumonia in subjects aged ≥65 years. Table 7S – Estimates from the multiple mixed-effect logistic regression model evaluating the likelihood of experiencing pneumonia in never smoking subjects.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Pegoraro, V., Heiman, F., Levante, A. et al. An Italian individual-level data study investigating on the association between air pollution exposure and Covid-19 severity in primary-care setting. BMC Public Health 21, 902 (2021). https://doi.org/10.1186/s12889-021-10949-9
- Air pollution
- Particulate matter
- Individual-level data