https://orcid.org/0000-0002-7556-6928
Figures
Abstract
Background
COVID-19 rapidly escalated into a pandemic, threatening 213 countries, areas, and territories the world over. We aimed to identify potential province-level socioeconomic determinants of the virus’s dissemination, and explain between-province differences in the speed of its spread, based on data from 36 provinces of Northern Italy.
Methods
This is an ecological study. We included all confirmed cases of SARS-CoV-2 reported between February 24th and March 30th, 2020. For each province, we calculated the trend of contagion as the relative increase in the number of individuals infected between two time endpoints, assuming an exponential growth. Pearson’s test was used to correlate the trend of contagion with a set of healthcare-associated, economic, and demographic parameters by province. The virus’s spread was input as a dependent variable in a stepwise OLS regression model to test the association between rate of spread and province-level indicators.
Results
Multivariate analysis showed that the spread of COVID-19 was correlated negatively with aging index (p-value = 0.003), and positively with public transportation per capita (p-value = 0.012), the % of private long-term care hospital beds and, to a lesser extent (p-value = 0.070), the % of private acute care hospital beds (p-value = 0.006).
Conclusion
Demographic and socioeconomic factors, and healthcare organization variables were found associated with a significant difference in the rate of COVID-19 spread in 36 provinces of Northern Italy. An aging population seemed to naturally contain social contacts. The availability of healthcare resources and their coordination could play an important part in spreading infection.
Citation: Buja A, Paganini M, Cocchio S, Scioni M, Rebba V, Baldo V (2020) Demographic and socio-economic factors, and healthcare resource indicators associated with the rapid spread of COVID-19 in Northern Italy: An ecological study. PLoS ONE 15(12): e0244535. https://doi.org/10.1371/journal.pone.0244535
Editor: William Joe, Institute of Economic Growth, INDIA
Received: April 27, 2020; Accepted: December 11, 2020; Published: December 28, 2020
Copyright: © 2020 Buja et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: Data analysis was performed on ecological and not on individual data, derived from the publicly available Italian National Institute of Statistics databases (http://dati.istat.it/Index.aspx, https://www.istat.it/it/archivio/14562) and from the Italian National Olympic Committee database (https://www.coni.it/it/registro-societa-sportive.html). Variables are listed in Tables 1 and 2. The authors of the present study have no special access privileges in accessing these datasets.
Funding: The author(s) received no specific funding for this work.
Competing interests: The authors have declared that no competing interests exist.
Background
The disease caused by the COVID-19 virus first detected in Wuhan, China, in December 2019 [1] escalated rapidly into a pandemic [2], coming to threaten 213 countries, areas and territories across the globe [3].
Italy recorded its first imported cases in January 2020 [4], and local cases emerged towards the end of February in Northern Italy, in the Lombardia and Veneto regions [4, 5]. Restrictions on international travel, domestic mobility, mass gatherings, and sporting events were gradually enacted by the Italian Government in an effort to contain the contagion [4], until the country’s complete lockdown on March 9th, 2020 [6]. Once implemented, these mitigation strategies (which sounded draconian to the population) proved effective in limiting the diffusion of COVID-19 to other parts of the country, but the infection continued to spread with different trends in Northern Italy.
Many demographic and socioeconomic factors, and healthcare organization variables can be associated with the variability of contagious disease propagation rates. It is important for public health systems to identify these factors and variables to produce solid evidence that can drive mitigation strategies. As well as having a direct influence on health outcomes in a given area, the findings of public health research also have important implications for ongoing global efforts to contain COVID-19. For example, a recent study found that even a small increase in long-term exposure to PM2.5 leads to a large increase in the virus-related death rate, with the magnitude of increase 20 times that observed for PM2.5 and all cause mortality [7].
Besides environmental conditions, factors such as employment rates, population density, and healthcare resources may help to account for the observed differences in rapidity of COVID-19 propagation.
With this ecological study we aimed to identify potential province-level demographic, socioeconomic and healthcare determinants of the COVID-19 pandemic’s dissemination, and to explain between-province differences in the rate of the disease’s spread, based on data from 36 provinces in four regions of Northern Italy.
Materials and methods
Context
The Italian National Health Service (NHS) was set up in 1978 with universal coverage, solidarity, human dignity, and human health as its guiding principles. It is regionally based, and organized at national, regional and local levels. Under the Italian Constitution, central government controls the distribution of tax revenue for publicly-financed healthcare and decides on a national statutory benefits package of “essential levels of care”—offered to all residents in every region. The country’s 19 regions and two autonomous provinces are responsible for organizing and delivering healthcare services through local health units. Depending on the region, public funds are allocated by local health units to public hospitals and accredited private clinics. In 2017, there were approximately 213,700 hospital beds, including day hospital and day surgery beds: 170,000 in public hospitals (2.8 beds per 1,000 population, including 2.5/1,000 for acute care and 0.3/1,000 for long-stay patients), and 43,700 in accredited private clinics (0.7 beds per 1,000 population, including 0.4/1,000 for acute care and 0.3/1,000 for long-stay patients) [8]. Public hospitals are either managed directly by local health units and coordinated with other local health services such as preventive medicine departments and primary care districts, or they operate as semi-independent public enterprises. A prospective payment system based on diagnosis-related groups (DRG) operates across the country and accounts for most hospital revenue, but is generally not applied to hospitals run directly by local health units, where global budgets are common.
Study design and data sources
In the present ecological study we included all confirmed cases of SARS-CoV-2 reported between February 24th and March 30th, 2020, drawn from the Italian Civil Protection Department’s official data [9]. This agency operates certified surveillance systems for the Italian government, with a two-stage testing procedure for SARS-CoV-2, one at local health unit level, and one at national laboratory level for confirmation purposes. All reported cases of COVID-19 infection met the laboratory criteria for confirmation established in the EU case definition procedure, which includes using RT-PCR.
We collected data on 36 provinces in four regions of Northern Italy–Piemonte, Lombardia, Veneto, and Emilia-Romagna–which have had the highest numbers of infections during the COVID-19 pandemic. In order to standardize as much as possible the period analyzed, we chose to evaluate the onset phase of the disease and we considered as starting point for each province i (T0 i) the first day when that province reached 10 cases. For each province, we calculated the trend of contagion was computed as the relative increase in the number of individuals infected between two time endpoints: T0 i and T20 i (20 days after T0). Demographics were extracted from the most recently updated databases on the website of the Italian Institute of Statistics (ISTAT) [10], and the ISTAT’s Health for All database [11]. The number of sporting associations for each province was obtained from the registry of the Italian National Olympic Committee (CONI) [12].
In particular, the variables introduced in the last model are defined as subsequent:
- Aging Index: Number of elderly in a population for 100 youngster, source ISTAT
- Population Density: Number of people on the province territory divided by the territorial extension, source ISTAT
- Density Inside Housing: Number of family members per square meter, source ISTAT
- % Sport Associations / Resident Population: Number of sport associations divided by the resident population, source CONI–Registro Società Sportive
- Unemployment Rate: Ratio of people searching for a job and the work forces population, source ISTAT
- Employment Rate: Ratio of employed people and the work forces population, source ISTAT
- % of Active Companies / Resident Population: Number of active companies divided by the resident population, source ISTAT
- Local Public Transport Use Per Capita: Local public transport demand in provinces capital towns, source ISTAT
- % of Private Long-term Care Beds: Percentage of hospital beds for long-term care in private clinics on the total number of hospital beds for long-term care, source Health For All–ISTAT
- % of Private Acute Care Hospital Beds: Percentage of hospital beds for acute care in private clinics on the total number of hospital beds for long-term care, source Health For All–ISTAT
- % of Anesthesiologists / Resident Population: Number of anesthesiologist divided by the resident population, source Health For All–ISTAT
- General Practitioners Rate: Number of general practitioners for 10.000 inhabitants, source Health For All–ISTAT
- Dentists Rate: Number of dentists for 10.000 inhabitants, source Health For All–ISTAT
- Nursing Staff Rate: Number of nurse for 10.000 inhabitants, source Health For All–ISTAT
- Acute Care Hospital Beds Rate: Acute care hospital beds for 10.000 inhabitants, source Health For All–ISTAT
- Long-term Care Hospital Beds Rate: Long-term care hospital beds for 10.000 inhabitants, source Health For All–ISTAT
- Emergency Medical Service Rate: Number of emergency medical services for 10.000 inhabitants, source Health For All–ISTAT
Data analysis
We modeled the disease’s spread assuming an exponential growth with the number of infections multiplying by a certain factor each day. Assuming no recoveries during the three-week period considered, the total number of people infected at the time t, y_t can be written as:
where y_1 is the number of cases on day 1 and β is the transmission factor.
Since we know the number of cases on day 1 and day 20, we can compute β_i, the transmission factor, for the i-th province as:
A descriptive analysis was conducted, calculating minimum, maximum and mean values, to illustrate the spread of the disease. Pearson’s test was applied to correlate β_i with a set of healthcare-associated, economic and demographic parameters at province level (see Table 1).
β_i was input as a dependent variable in an OLS regression model to study the association between the rate of spread of COVID-19 infection and the variables that correlated with the dependent variable on a bivariate level with a p-value of less than 0.20. Density inside housing and the unemployment rate were not included in the multivariate model due to multicollinearity issues. A stepwise selection procedure was used. A residuals analysis was performed to check the assumptions of the model. This analysis included testing for the residuals’ normal distribution (Shapiro-Wilk normality test) and independence (Durbin-Watson test).
A sensitivity analysis was conducted, considering different dependent variables as rates (t1/t0) or incremental rates ([t1-t0]/t0), which produced the same results (data not shown). All data obtained were encoded in a master sheet using a Microsoft Office Excel spreadsheet (Version 2016, Microsoft Corporation, Redmond, WA). Data analysis was done with Stata 14 statistical software (Stata Corp LP, College Station, TX) and R Studio Integrated Development for R (RStudio, Inc., Boston, MA).
Results
Table 2 shows that the increase in the numbers of cases in 20 days varied considerably among the various provinces from a minimum of 115 cases in the province of Rovigo to a maximum of 3398 cases in Bergamo, with an average growth of 852 cases per province in 20 days. The rate of this increase was lowest in the province of Asti (4.5%) and highest in Brescia (330%), with an average over all provinces of 59.9%.
The scatter plot (Fig 1) shows significant correlations between the increase in transmission and several demographic, socioeconomic and healthcare variables.
The spread of COVID-19 was correlated negatively with aging index and positively with the other indicators showed.
A significant negative correlation emerged between increment β and population aging rates, while a significant positive correlation was found with employment rates, public transportation rates, in-house density, population density, and the proportions of private acute and long-term care beds in clinics and nursing homes.
Table 3 shows the bivariate and multivariate analyses, demonstrating that the spread of COVID-19 in the regional setting considered was negatively associated (p-value = 0.003) with the population aging index, i.e. the ratio of the number of people aged 65 or more to the number of people aged 0–14 years multiplied by 100. It was positively associated with: public transportation per capita, i.e. the ratio of the yearly number of public transport passengers to the total population(p-value = 0.012). It was also positively associated, albeit to a lesser extent (p-value = 0.070), with the % of private long-term care hospital beds, and more strongly with the % of private acute care hospital beds (p-value = 0.006).
The Shapiro-Wilk test showed that the residuals had a normal distribution (p-value = 0.870) and the Rubin-Watson test showed that the residuals were independent (p-value = 0.48). The high value of the adjusted R squared indicates that the model had a good fit.
Discussion
With this ecological study, we found correlations between higher SARS-CoV-2 contagion rates and certain demographic, socioeconomic and healthcare factors in the Italian provinces considered. Unlike the first cases imported from outbreak regions, the local spread of communicable respiratory diseases seems to favor certain patterns and conditions, especially relating to social contact patterns [13].
We found a negative association between COVID-19 contagion rates and aging. This seems somewhat at odds with the currently-known epidemiology of the virus, that in Italy is more likely to infect and kill older people. The median age of people infected with COVID-19 nationwide is 62 years old, and the median age of the related deaths is 80 years old [14]. On the other hand, provinces with a higher aging index had a more gentle contagion curve. A possible explanation is that older people tend to move less outside their home and travel less far from where they live [15]. They are also less likely to take part in mass gatherings or social events, whereas younger people tend to spend more time with others outside their family [16]. Taken together, these behavioral factors could predispose the elderly to social distancing and self-isolation, naturally containing the propagation of communicable diseases in the provinces with the oldest populations.
In our analysis, higher levels of employment, public transportation usage, in-house density, and population density correlated positively with the spread of infection. What these socioeconomic factors have in common is the mobility of individuals and their exposure to close social contacts, both of which facilitate the propagation of SARS-CoV-2. This is consistent with recent analyses that associated a higher risk of COVID-19 transmission in mainland China with the use of trains [17], buses and flights from Wuhan, and the risk was higher the longer a journey lasted [18]. In-house density has also already been identified as a factor predisposing to higher transmission rates of pandemic influenza in Great Britain [19].
A greater spread of the virus was seen in provinces with larger proportions of private acute care beds in clinics, and long-term care beds in nursing homes. These data should not be interpreted as a lack of involvement of the private sector in providing hospital care for COVID patients. These facilities are an important part of the services provided by the Italian health system and they responded promptly to the emergency situation. They have ordinary and intensive care beds for COVID patients just like public hospitals. The explanation should be sought by looking at the issue more in depth and speculating on the outcome variable analyzed in the present study, i.e. the rate of spread of the infection. Hospitals, clinics and care homes are only responsible for containing it within their walls, while the preventive medicine departments of local health units are in charge of wider containment efforts, based on epidemiological intelligence. It may be that provinces with higher proportions of private healthcare facilities experienced more communication difficulties with the reporting of cases to preventive medicine departments to gather epidemic intelligence because private facilities are usually less integrated with preventive medicine and community healthcare services than public hospitals. This could hamper the early identification of patients through the sort of epidemic intelligence activities strongly recommended by the WHO in the case of COVID-19 [20]. Another explanation could relate to how private clinics are funded, i.e. based on DRG-based tariffs for each patient admitted. Owing to incentives linked to this particular funding mechanism, provinces with a relatively high concentration of private hospital beds might have paid less attention to avoiding the hospitalization of patients with COVID-19 symptoms who could have been successfully treated at home with the involvement of community and primary care services. In fact, primary healthcare services have a crucial role in ensuring that people caring for a family member suffering from COVID 19 manage their contact with the patient appropriately, and follow national or local policies regarding home quarantine. Caregivers should wear medical masks or the best available protection against respiratory droplets when in close contact with the patient, and observe hand hygiene recommendations. The WHO emphasized the importance of notifying healthcare providers of the diagnosis in order to receive instructions on where to seek care, when and where to enter healthcare facilities, and what precautions to take to prevent and control the infection [20]. Our findings would suggest the importance of improving communications between private clinics and the Italian NHS responsible for their accreditation (i.e. the issue of the public license to provide hospital care on the NHS’s behalf) to ensure a real consistency of their activities with regional health plans (including their integration with the preventive medicine and primary care services of local health units), and an effective assessment of their activities and results.
This study has at least two important limitations. First, we collected data from available health indicators, which do not measure all the phenomena relevant to understanding and explaining the findings, even though our models performed well in explaining their variability (MR2 = ~60%). Second, the real number of SARS-CoV-2 contagions is known to be underestimated because not all individuals in the population considered were screened, and swabs were handled differently by the various regional healthcare systems. Even with such well-known limitations [21], ecological studies can still help healthcare workers and stakeholders to contain infections and fight pandemics, especially in the early stages of emerging diseases when clinical data are still limited [22].
Conclusions
Several demographic, socioeconomic and healthcare factors were found associated with significant differences in the rate of COVID-19 spread in 36 provinces of Northern Italy. An older population seemed to naturally slow the contagion due to fewer social contacts. Socioeconomic factors (especially the rates of employment and public transportation usage) and organizational features of local healthcare systems (particularly the proportion of private healthcare facilities) also seemed to play an important part in the infection’s spread. These correlations could be helpful for the purpose of designing measures to reduce the propagation of COVID-19 in other countries now dealing with the first phase of this pandemic.
References
- 1.
World Health Organization. Disease Outbreak News Update–January 12, 2020: Novel Coronavirus–China. [cited 2020 April 13]. Available from: https://www.who.int/csr/don/12-january-2020-novel-coronavirus-china/en/
- 2.
World Health Organization. WHO Director-General's opening remarks at the media briefing on COVID-19–11 March 2020. [cited 2020 April 13]. Available from: https://www.who.int/dg/speeches/detail/who-director-general-s-opening-remarks-at-the-media-briefing-on-covid-19—11-march-2020
- 3.
World Health Organization. Coronavirus disease (COVID-19) Pandemic–April 12 Update. [cited 2020 April 13]. Available from: https://www.who.int/emergencies/diseases/novel-coronavirus-2019
- 4.
Italian Ministry of Health. COVID-19 –Italian Situation. [cited 2020 April 13]. Available from: http://www.salute.gov.it/portale/nuovocoronavirus/dettaglioContenutiNuovoCoronavirus.jsp?area=nuovoCoronavirus&id=5351&lingua=italiano&menu=vuoto
- 5.
Italian Civil Protection. Coronavirus disease 2019 (COVID-19) Italian Situation Map. [cited 2020 April 13]. Available from: http://opendatadpc.maps.arcgis.com/apps/opsdashboard/index.html#/b0c68bce2cce478eaac82fe38d4138b1
- 6.
The Prime Minister of the Italian Republic–Decree March 9th, 2020. Urgent Measures of Containment and Management of Epidemiological Emergency COVID-19, further dispositions. [cited 2020 April 13]. Available from: https://www.gazzettaufficiale.it/eli/id/2020/03/09/20A01558/sg
- 7.
Wu X, Nethery RC, Sabath MB, Braun D, Dominici F. Exposure to air pollution and COVID-19 mortality in the United States. MedRxiv 2020.04.05.20054502 [Preprint]. doi: https://doi.org/https://doi.org/10.1101/2020.04.05.20054502
- 8.
Italian Ministry of Health. Annuario statistico del Servizio Sanitario Nazionale, Assetto organizzativo, attività e fattori produttivi del SSN—Anno 2017. Published on June 11th, 2018. [cited 2020 April 13]. Available from: http://www.salute.gov.it/portale/documentazione/p6_2_2_1.jsp?lingua=italiano&id=2370
- 9.
Italian Civil Protection–Contagion Rates. [cited 2020 April 13]. Available from: https://github.com/pcm-dpc/COVID-19/tree/master/dati-province
- 10.
National Institute of Statistics (ISTAT). Main database. [cited 2020 April 13]. Available from: http://dati.istat.it/Index.aspx
- 11.
National Institute of Statistics (ISTAT). Health for All database. [cited 2020 April 13]. Available from: https://www.istat.it/it/archivio/14562.
- 12.
Italian National Olympic Committee (CONI). Registry of Sports Associations. [cited 2020 April 13]. Available from: https://www.coni.it/it/registro-societa-sportive.html
- 13. Leung K, Jit M, Lau EHY, Wu JT. Social contact patterns relevant to the spread of respiratory infectious diseases in Hong Kong. Sci Rep. 2017;7(1):7974. Published 2017 Aug 11. pmid:28801623
- 14.
Istituto Superiore di Sanità. Characteristics of SARS-CoV-2 patients dying in Italy. Report based on available data on April 13th, 2020. [cited 2020 April 13]. Available from: https://www.epicentro.iss.it/en/coronavirus/bollettino/Report-COVID-2019_13_april_2020.pdf
- 15.
Mollenkopf H, Hieber A, Wahl H-W. Continuity and change in older adults’ out-of-home mobility over ten years: a qualitative-quantitative approach. In: Meusburger P, Werlen B, Suarsana L, editors. Knowledge and Action. Knowledge and Space. Cham: Springer International Publishing; 2017:267–289. https://doi.org/10.1007/978-3-319-44588-5_15
- 16. Marcum CS. Age differences in daily social activities. Res Aging. 2013;35(5):612–640. pmid:25190898
- 17. Zhao S, Zhuang Z, Ran J, Lin J, Yang G, Yang L. The association between domestic train transportation and novel coronavirus (2019-nCoV) outbreak in China from 2019 to 2020: a data-driven correlational report. Trav Med Infect Dis. 2020;30:101568 pmid:32006656
- 18. Zheng R, Xu Y, Wang W, Ning G, Bi Y. Spatial transmission of COVID-19 via public and private transportation in China [published online ahead of print, 2020 Mar 14]. Travel Med Infect Dis. 2020;101626. pmid:32184132
- 19. House T, Keeling MJ. Household structure and infectious disease transmission. Epidemiol Infect. 2009;137(5):654–661. pmid:18840319
- 20.
World Health Organization. (2007). Infection prevention and control of epidemic- and pandemic-prone acute respiratory diseases in health care: WHO interim guidelines. Geneva: World Health Organization. [cited 2020 April 13]. Available from: https://apps.who.int/iris/handle/10665/69707
- 21. Cohen HW. Limitations of an ecological study: a review. Am J Hypertens. 2005;18(6):750. pmid:15925730
- 22. Fauci AS. Pandemic influenza threat and preparedness. Emerg Infect Dis. 2006;12:73–7 pmid:16494721