key: cord-0976535-qij4d3a2 authors: Vignatelli, Luca; Zenesini, Corrado; Belotti, Laura M.B.; Baldin, Elisa; Bonavina, Giuseppe; Calandra‐Buonaura, Giovanna; Cortelli, Pietro; Descovich, Carlo; Fabbri, Giovanni; Giannini, Giulia; Guarino, Maria; Pantieri, Roberta; Samoggia, Giuseppe; Scaglione, Cesa; Trombetti, Susanna; D'Alessandro, Roberto; Nonino, Francesco title: Risk of Hospitalization and Death for COVID‐19 in People With Parkinson's Disease or Parkinsonism date: 2020-12-02 journal: Mov Disord DOI: 10.1002/mds.28408 sha: 60dfa9255a51a6cccd8d7e958ea421b9917fb525 doc_id: 976535 cord_uid: qij4d3a2 BACKGROUND: The risk of COVID‐19 and related death in people with Parkinson's disease or parkinsonism is uncertain. The aim of the study was to assess the risk of hospitalization for COVID‐19 and death in a cohort of patients with Parkinson's disease or parkinsonism compared with a control population cohort, during the epidemic bout (March–May 2020) in Bologna, northern Italy. METHODS: Participants of the ParkLink study with the clinical diagnosis of Parkinson's disease or parkinsonism and people anonymously matched (ratio 1:10) for sex, age, district, and Charlson Index were included. The hospital admission rate for COVID‐19 (February 26–May 31, 2020) and the death rate for any cause were the outcomes of interest. RESULTS: The ParkLink cohort included 696 subjects with Parkinson's disease and 184 with parkinsonism, and the control cohort had 8590 subjects. The 3‐month hospitalization rate for COVID‐19 was 0.6% in Parkinson's disease, 3.3% in parkinsonism, and 0.7% in controls. The adjusted hazard ratio (age, sex, district, Charlson Index) was 0.8 (95% CI, 0.3–2.3, P = 0.74) in Parkinson's disease and 3.3 (1.4–7.6, P = 0.006) in parkinsonism compared with controls. Twenty‐nine of the infected subjects died; 30‐day fatality rate was 35.1%, without difference among the 3 groups. Six of 10 Parkinson's disease/parkinsonism patients had the infection during hospitalization or in a nursing home. CONCLUSIONS: Parkinson's disease per se probably is not a risk factor for COVID‐19 hospitalization. Conversely, parkinsonism is an independent risk factor probably because of a more severe health status, entailing higher care dependence and placement in high‐infection‐risk accommodations. © 2020 The Authors. Movement Disorders published by Wiley Periodicals LLC. on behalf of International Parkinson and Movement Disorder Society. elderly affected by chronic conditions. [6] [7] [8] Persons with Parkinson's disease (PD) and parkinsonisms (PS) show a paradigmatic pattern of vulnerability, being elderly and often with multiple age-related comorbidities, including cognitive impairment and higher polypharmacy risk. 9 Whether people with PD or PS have an increased risk of COVID-19 and an associated higher risk of death is still debatable. [10] [11] [12] [13] [14] A telephone survey, including an outpatient clinic cohort of people with PD and their family members as controls, reported a similar prevalence of "probable" or "confirmed" COVID-19. 11 Conversely, another telephone survey suggested a higher prevalence of COVID-19 compared with that expected for the general population. 15 Trying to overcome the limitations of these studies, we designed a population-based study investigating the risk for hospital admission of COVID-19 and the risk of death in people with PD or PS. ParkLink Bologna is a population-based record linkage project started in 2016 in the province of the city of Bologna, Italy, including a cohort of people with a diagnosis of PD or PS confirmed by a neurologist according to international criteria. In ParkLink Bologna clinical data of individual patients are linked to administrative databases. 16 The aim of the present study was to assess the risk of hospital admission for COVID-19 and related death in the ParkLink Bologna cohort of patients with PD or PS compared with a control population cohort during the COVID-19 epidemic bout in the local health trust of Bologna (March to May 2020). The Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) guidelines 17 were followed. A historical cohort design was applied. The local health trust (LHT) of Bologna, in northern Italy, had a population of 886,000 (752,104 adults), on December 31, 2019. At the beginning of the outbreak, the area included 9 general hospitals, 3 academic hospitals (1 general, 1 neurologic, 1 orthopedic), 11 private clinics operating within the Italian National Health Service, for a total of about 4400 beds, of which about 120 were intensive care. In response to the COVID-19 emergency, hospital facilities underwent major reorganization; since mid-March, 8 COVID-19 hospitals/ wards have been created, with a total of about 400 beds for symptomatic COVID patients and about 150 COVID-19 intensive care unit beds. 18 Dedicated diagnostic protocols for suspected COVID-19 patients to be applied at every emergency department access were explicitly implemented by the LHT of Bologna. This study is based on the ParkLink Bologna project (https://www.ausl.bologna.it/isnb/ricerca/parklink), an ongoing record linkage system for PD and PS approved in 2015 by the local Ethics Committee (see Statements of Ethics section). The project population includes the whole population living in the area covered by the LHT of Bologna during the period of observation, July 1, 2015-June 30, 2030. For the purpose of the present study we included people with PD or PS (ParkLink cohort) and people anonymously matched from the general population (control cohort), alive on December 31, 2019. People included in the ParkLink cohort are both prevalent and incident cases of PD or PS recruited consecutively after signing a form giving consent to link their diagnosis to their stored administrative data for any possible research purpose. Since 2015, all neurologists working in the area of the LHT of Bolognaincluding 3 hospital outpatient services dedicated to movement disorders and several public and private outpatient serviceshave been invited to enroll people with PD or PS, during clinical practice. Up to December 2019, 24 neurologists have voluntarily joined the project. The following data have been recorded in an e-CRF and stored in a secure database linked to the other administrative databases: unique anonymous identification code, date of birth, diagnosis, year of onset, motor symptoms at onset (tremor or bradykinesia), side of onset (unilateral or bilateral), and Hoehn and Yahr score. Diagnosis of PD is defined by the recruiting neurologist according to Gelb criteria. 19 Specified PS, [20] [21] [22] [23] other rare neurological causes affecting the basal ganglia, or unspecified PS are recorded in the system and labeled as parkinsonism tout-court (PS). People with drug-induced PS are excluded. The control cohort consisted of people matched with the ParkLink cohort, with a ratio of 1:10 for age, sex, district of residence, and comorbidity according to the Charlson Index. 24 Subjects who used drugs for PD (levodopa, dopaminergics, monoamine oxidase-B inhibitors) for at least 180 consecutive days during 2019 were excluded. Being that the Italian health system is universal, all resident citizens are anonymously recorded in a homogenous way when accessing any public or private health facility and personal data stored by the qualified LHT in secure databases. The ParkLink Bologna project allows, for each individual, the linkage of the clinical diagnosis with administrative databases and the extraction of population-matched controls. The Unità Operativa Epidemiologia e Statistica of the IRCCS Istituto delle Scienze Neurologiche di Bologna (part of the LHT of Bologna) is in charge of the ParkLink Bologna project and has access to several administrative databases (drug prescription, hospital discharge, copayment exemption, medical home care, mortality). The Charlson Index was based on data from the hospital discharge database, and the index was calculated according to the adapted method for International Classification of Disease-9-CM (ICD-9-CM). 25 All hospital facilities, both public and private, of the LHT of Bologna provide data on hospital discharge. Hospital admission rate for COVID-19 from February 26 to May 31, 2020, was the primary outcome. Diagnosis was obtained from hospital discharge data. The ICD-9-CM codes for the diagnosis of COVID-19 were 079.82 and 480.3, according to the Emilia-Romagna region regulation (March 10, 2020). Hospital admission rate for any reason was assessed to estimate the general trend of hospitalization in the same time window. We assessed also the death rate for any cause from inception until June 30, 2020, as reported by the demographic database, and measured as time to event from COVID-19 hospitalization (index date) to death or the end of follow-up, whichever came first. When possible, qualitative assessment of hospital charts was carried out for patients with PD/PS and COVID-19 to clarify in which circumstances COVID-19 infection occurred. In the descriptive analysis, the characteristics of the 2 cohorts are presented as mean ± standard deviation (SD) for the continuous variables and with absolute (n) and relative (%) frequency for categorical variables. The chi-square test and Kruskal-Wallis test were used to evaluate the univariable association between conditions (PD/PS/controls) and categorical or continuous variables, respectively. To evaluate the outcomes (hospitalization for COVID-19 and death), 2 analyses were performed. In the first one, the time to entry in the analysis was February 26, 2020, and the time to end point was the date of hospital admission or May 31, 2020. In the second analysis, the time to entry was the date of hospital admission, and the time to end point was the date of death or June 30, 2020. We used univariable and multivariable Cox regression models to estimate the hazard ratio (HR) with 95% confidence interval (95% CI) associating COVID-19 hospitalization and death in the presence of PD/PS. Results of the univariable analysis are presented as Kaplan-Meier curves. In the multivariable models we included sex, age, Charlson Index, and district of residence as covariates to take into account the impact of the comorbidity. The proportional hazards assumption was tested (P > 0.05) using Schoenfeld residuals. Data linkage and statistical analyses were conducted using Stata SE version 14.2. 26 The ParkLink cohort included 696 subjects with PD (mean age, 75.0 years) and 184 with PS (mean age, 80.5 years), followed up for 2149 and 561 personmonths, respectively. The control cohort included 8590 subjects (mean age, 76.0 years; 26,506 person-months of follow-up), after the exclusion of 27 subjects on treatment with drugs for PD. Demographic and clinical characteristic data are summarized by group in Table 1 . People with PS were older, with higher Hoehn and Yahr scale score and with more frequent comorbidities than people with PD. The epidemic bout is depicted in Figure 1 as crude number of hospital admissions for COVID-19 in the whole population of the LHT of Bologna. In our cohorts, the first COVID-19 patient was admitted on February 27, 2020. Hospital admissions for COVID-19 were 74: 4 among PD subjects (3 men), 6 among PS subjects (3 men), and 64 in the control cohort (44 men). The 3-month admission rates were 0.6% in PD (95% CI, 0.2%-1.5%), 3.3% (95% CI, 1.5%-7.1%) in PS, and 0.7% (95% CI, 0.6%-0.9%) in controls. Kaplan-Meier curves showed a progressively higher risk of hospitalization for PS compared with the PD group and the control group during the epidemic bout (Fig. 2) . The crude HR compared with the control cohort was 0.8 (95% CI, 0.3-2.1; P = 0.61) for PD and 4.4 (95% CI, 1.9-10.2; P = 0.01) for PS. Univariable analysis also found the following factors associated with hospital admission ( Table 2) : age, Charlson Index, district of residence, myocardial infarction, congestive heart failure, cerebrovascular disease, dementia, diabetes, renal disease, and any malignancy. After adjustment (Table 2) , HR was 0.8 (95% CI, 0.3-2.3; P = 0.74) in PD and 3.3 (95% CI, 1.4-7.6; P = 0.006) in PS. In the same period, the adjusted risk of hospital admission for any reason was slightly increased in PD subjects (HR, 1.51; 95% CI, 1.1-2.1; P = 0.011) but not in PS subjects (HR, 1.2; 95% CI, 0.7-2.2; P = 0.42). Considering only the ParkLink cohort, the univariable analysis with the Hoehn and Yahr scale score as factor found an HR for hospital admission of 1.63 (95% CI, 0.97-2.74; P = 0.064). The overall mean length of hospital stay was 15.0 ± 11.5 days (range, 1-52 days) in controls, 29.5 ± Movement Disorders, 2020 Overall, 29 of the infected subjects died; 30-day fatality rate was 35.1% (95% CI, 25.5%-47.2%). Twenty-five controls (39%) died a mean of 16 ± 14.8 days after hospital admission (range, 0-60 days): The 4 people with PD and COVID-19 infection had a mean Hoehn and Yahr scale score of 2.37 ± 0.63 and showed the following comorbidities (Table 3) : 2 with congestive heart failure, 1 with cerebrovascular disease, 1 with dementia, and 1 with peripheral vascular disease. One female subject had 3 comorbidities and died after 11 days of hospitalization. The 6 people with PS had a mean Hoehn and Yahr scale score of 3.67 ± 0.82 and had the following diagnoses (and comorbidities; Table 3 ): 3 vascular parkinsonism (3 cerebrovascular disease, 2 dementia, 1 renal disease), 1 progressive supranuclear palsy (none), 1 multiple system atrophy (cerebrovascular disease, dementia, diabetes), 1 undetermined parkinsonism (perinatal hydrocephalus, dementia). The latter 3 subjects died after hospitalization for 6 days (2 subjects) and 21 days (1 subject). In 6 patients the COVID-19 infection occurred during hospitalization or as a guest in a nursing home. In the other 2 cases it occurred probably when attending senior citizens/day care centers. During the SARS-CoV-2 epidemic bout in Emilia-Romagna region, Italy, we observed a 3-month (March-May 2020) hospital admission rate for COVID-19 of 0.6% in people with PD, not different from a control population. Conversely, people with PS showed a 3-fold risk (admission rate of 3.3%) after adjustment for comorbidities. Both the ParkLink cohort and the control cohort included mainly elderly people, and we found that COVID-19 risk is associated with age and number of comorbidities. Finally, the risk of death after hospitalization was high in both cohorts (around 35%). No association with PD or PS was detected, possibly because of the small number of events in these groups. It has been hypothesized that patients with PD or PS could have a heightened risk of COVID-19 13, 27 because of the indirect association with old age, multiple comorbidities, and polypharmacy frequency. 9 Respiratory dysfunction, common among people with PD or PS, may represent a direct risk factor for worse respiratory complications of COVID-19 infection. 12, 27 Our data suggest that PD per se is not associated with a higher risk for COVID-19 hospitalization. A case-controlled telephonic survey by Fasano et al 11 found that patients with PD (selected from a movement disorder tertiary center in the Lombardy region of Italy) did not have a COVID-19 risk (7.1%) different from that of family members (7.6%), defined as "probable" or "confirmed". Conversely, another Italian study reported a higher risk of infection in a movement disorder tertiary center cohort of people with PD, making an indirect comparison of the prevalence found in the cohort with that expected in the general population. Recently, a telephone survey 28 on a population of PD patients in the Piedmont region of Italy reported a prevalence of 0.57% COVID-19 positivity in a population of PD patients that was indirectly compared with the 0.63% prevalence of the general population. All studies were limited by the method of ascertainment of the infection (self-reported in most cases), 1 study 11 by the controls (family members) probably sharing the same environmental risk of infection, therefore increasing the risk of recall bias, and 2 studies ,28 by the indirect comparison of the prevalence in the general population without adjusting for other factors. Our novel observation of a higher risk of COVID-19 hospitalization for people with PS could be explained by the following hypotheses. First, the category of PS, even if including a heterogeneous spectrum of disorders (some specific such as progressive supranuclear palsy or multiple system atrophy, other less specific "vascular parkinsonism" or "unspecified" PS), can be considered more vulnerable than PD because of a relatively more rapid disease progression, [29] [30] [31] [32] [33] [34] reaching milestones such as dysphagia, falls, cognitive impairment, and psychiatric symptoms earlier. Second, people with PS, showing on average higher disability and being cognitively impaired more often than PD patients, have a higher probability of being admitted to nursing homes. 35 Risk of SARS-CoV-2 infection in Italy was higher among nursing home guests 36, 37 ; therefore, people with PS could have had this risk heightened by being part of this population. In our study, indeed at least 6 of 10 patients were infected as guests in the hospital or in the nursing home. Moreover, patients with atypical parkinsonism may more frequently require tracheostomy, percutaneous endoscopic gastrostomy, urinary catheterization, or supplementary nutrition and/or rehydration. These procedures need regular (eg, weekly or monthly) assistance for their revision, checkup and replacement of consumables (eg, cannula, indwelling catheter, drainage bag) from specialized professional nurses and palliative care specialists. Finally, we can hypothesize that during the lockdown period people with PS had a higher risk of hospitalization for any reason; however, we did not find any difference. Risk factors for severe outcomes in hospitalized COVID-19 patients are reported to be chronic kidney disease, chronic obstructive pulmonary disease, cerebrovascular disease, and coronary heart disease. 38 However, risk factors for hospitalization for COVID-19 in the general population are less extensively studied and are reported to be obesity, chronic kidney disease, diabetes, hypertension, and asthma after adjustment for age, sex, and ethnicity. 39 In our 2 cohorts of elderly people (with or without PD/PS), neurologic (cerebrovascular disease, dementia) or nonneurologic comorbidities (myocardial infarction, congestive heart failure, diabetes, renal disease, malignancies) were found to be strongly associated with worse outcomes among hospitalized patients with COVID-19. Case fatality after hospital admission for COVID-19 is extremely variable and affected by many factors, 40 such as case definition, case mix of patients, geographical context, and capacity of the health care system. 41 In Italy, namely, in Lombardy during the epidemic bout, 30-day fatality rate was 19.7% in patients admitted to a large general hospital (median age, 61 years). 42 A large population-based study network of hospitalized patients (median age, 69 years) of northern Italy found a 30-day case fatality rate of 27.6%. 40 In other European countries (Germany, Spain, England), 43-45 the 30day case fatality rate ranged between 14% and 27% in patients with varying mean age (63-73 years). We found an impressive 30-day fatality rate of 35% that could be explained by the criteria of selection of our cohorts, which included people matched to the target population of PD patients and thus aged people with a relatively higher burden of comorbidities. No difference was found between the ParkLink cohort and the control cohort, possibly because of the insufficient power to detect differences. Two studies 28,46 suggested a possibly higher risk of death in PD patients after COVID-19 infection. However, serious methodological limitations in both studies (eg, PD diagnosis not on a clinical basis, 46 no adjustment for comorbidities, 28, 46 telephonic ascertainment of COVID-19 infection 28 ) did not allow drawing conclusions. Evidence suggests that subjects with chronic and serious neurodegenerative diseases experienced worsening of their clinical status during the COVID-19 epidemic. A telemedicine assessment survey on carers of 32 patients with frontotemporal dementia during the COVID-19 emergency in Italy showed worsening in behavior, language, and cognitive function in most cases. 47 Among COVID-19-infected patients such impact can be explained as a possible direct cause of the infection itself, but in uninfected patients it may be indirectly because of a generalized restriction of care provision (see Ferini-Strambi and Salsone 48 ). Indeed, people with PD/PS or with dementia require follow-up and supportive care provided in specialized health care settings, and in many countries, including Italy, the cancellation of nonurgent outpatient appointments could have caused a drastic reduction of such dedicated care. Population-based studies on the direct and indirect impact of the COVID-19 epidemic on neurodegenerative diseases are needed. The main strength of this study is having the clinical diagnosis of PD and PS performed by a neurologist according to international criteria. In a retrospective study, diagnostic assessment by a qualified neurologist is highly preferable than relying on other methods based on administrative claims data only. Second, it is a population-based cohort study including patients with PD and PS representative of the entire spectrum of patients (not selected from a movement disorder tertiary centre), compared with a matched general population cohort. Third, the diagnosis of COVID-19 infection is obtained from hospital discharge-coded administrative data, therefore excluding the risk of recall bias. Our study has some limitations. The use of ICD codes to label comorbidities could have underestimated the presence of other diseases in patients with a burdensome condition such as PD or PS. In particular, people with PS could have a residual presence of comorbidities not recorded by the SDO system. However, the different risk found for PD and PS suggests that this factor could have biased the results only in a small part. Moreover, the ICD code to diagnose COVID-19 (079.82, SARS-associated coronavirus; and 480.3, "Pneumonia due to SARS-associated coronavirus) has been recommended with little delay compared with the epidemic start in our region and has not been validated for COVID-19. Nevertheless, coding errors would have the same random probability in the 2 cohorts. In addition, because of the difficulty of the health system in dealing with the epidemic bout, we cannot exclude that some COVID-19 symptomatic patients were not hospitalized. However, because of the prompt strengthening of the local health facilities dedicated to COVID-19, we can assume that only mild symptomatic patients have been missed. In any case, there is no reason that the hospitalization risk could be biased in favor/disfavor of one of the groups of patients included. Finally, an updated evaluation of the Hoehn and Yahr scale score was suboptimal (in 36% of the ParkLink cohort patients the last available Hoehn and Yahr evaluation was before 2018); thus, the association of this factor to hospitalization could not be reliably tested. In conclusion, our study found that PD is not per se a risk factor for COVID-19 hospitalization, whereas, conversely, being affected by PS is an independent risk factor. A more severe health status of these patients, entailing higher care dependence and placement in accommodations (eg, hospitals, nursing homes) with high infection risk may explain such association. These aspects should be confirmed in larger population-based cohorts and studied for a longer period. The Global epidemiology of coronavirus disease 2019 (COVID-19): disease incidence, daily cumulative index, mortality, and their association with country healthcare resources and economic status Critical care utilization for the COVID-19 outbreak in Lombardy, Italy: early experience and forecast during an emergency response Covid-19 outbreak progression in Italian regions: approaching the peak by the end of march in northern Italy and first week of April in southern Italy Case-Fatality Rate BS. Characteristics of patients dying in relation to COVID-19 in Italy Clinical characteristics of hospitalized individuals dying with COVID-19 by age Group in Italy Global, regional, and national estimates of the population at increased risk of severe COVID-19 due to underlying health conditions in 2020: a modelling study Co-morbidity and polypharmacy in Parkinson's disease: insights from a large Scottish primary care database Outcome of Parkinson's disease patients affected by COVID-19 COVID-19 in Parkinson's disease patients living in Lombardy, Italy The impact of the COVID-19 pandemic on Parkinson's disease: hidden sorrows and emerging opportunities Impact of the COVID-19 pandemic on Parkinson's disease and movement disorders May Parkinson's disease be a protective factor against CNS involvement by SARS-CoV2? Prevalence and impact Q12 of COVID-19 in Parkinson's disease: evidence from a multicenter survey in Tuscany region Use of antidepressants and the risk of Parkinson's disease in the local health Trust of Bologna: a historical cohort study The strengthening the reporting of observational studies in epidemiology (STROBE) statement: guidelines for reporting observational studies Emergenza Covid-19. Piano per la gestione dei letti ospedalieri Diagnostic criteria for Parkinson disease Second consensus statement on the diagnosis of multiple system atrophy Clinical research criteria for the diagnosis of progressive supranuclear palsy (Steele-Richardson-Olszewski syndrome): report of the NINDS-SPSP international workshop Criteria for the diagnosis of corticobasal degeneration Clinicopathological investigation of vascular parkinsonism, including clinical criteria for diagnosis A new method of classifying prognostic comorbidity in longitudinal studies: development and validation Adapting a clinical comorbidity index for use with ICD-9-CM administrative databases StataCorp. Stata Statistical Software COVID-19: an early review of its global impact and considerations for Parkinson's disease patient care COVID-19 in Parkinson's disease: report on prevalence and outcome Survival and causes of death among people with clinically diagnosed Synucleinopathies with parkinsonism: a population-based study Predictors of survival in progressive supranuclear palsy and multiple system atrophy: a systematic review and meta-analysis The clinical features of pathologically confirmed vascular parkinsonism Palliative care and its emerging role in multiple system atrophy and progressive Supranuclear palsy Clinical outcomes of progressive supranuclear palsy and multiple system atrophy Clinical outcomes of two main variants of progressive supranuclear palsy and multiple system atrophy: a prospective natural history study Care of Late Stage Parkinsonism (CLaSP) Consortium. Characteristics of patients with latestage parkinsonism who are nursing home residents compared with those living at home Nursing homes or besieged castles: COVID-19 in northern Italy High impact of COVID-19 in long-term care facilities, suggestion for monitoring in the EU/EEA Epidemiological, comorbidity factors with severity and prognosis of COVID-19: a systematic review and metaanalysis Risk factors for COVID-19-associated hospitalization: COVID-19-associated hospitalization surveillance network and behavioral risk factor surveillance system Survival of hospitalized COVID-19 patients in Northern Italy: a population-based cohort study by the ITA-COVID19 Network2020 The many estimates of the COVID-19 case fatality rate 30-day mortality in patients hospitalized with COVID-19 during the first wave of the Italian epidemic: a prospective cohort study Clinical course and factors associated with outcomes among 1904 patients hospitalized with COVID-19 in Germany: an observational study Epidemiological characteristics of the COVID-19 outbreak in a secondary hospital in Spain Clinical characteristics and 28-day mortality of medical patients admitted with COVID-19 to a Central London teaching hospital Coronavirus disease 2019 Case Fatality and Parkinson's disease Telemedicine for delivery of Care in Frontotemporal Lobar Degeneration during COVID-19 pandemic: results from southern Italy COVID-19 and neurological disorders: are neurodegenerative or neuroimmunological diseases more vulnerable? SGML and CITI Use Only DO NOT PRINT Author Roles Luca Vignatelli: research project conception, organization and execution research project conception; statistical analysis design and execution; writing the first draft of the manuscript statistical analysis execution; review and critique of the manuscript. Elisa Baldin: research project conception; statistical analysis design; review and critique of the manuscript. Giuseppe Bonavina: research project execution research project organization and execution; review and critique of the manuscript. Pietro Cortelli: research project organization and execution; review and critique of the manuscript. Carlo Descovich: research project organization; review and critique of the manuscript. Giovanni Fabbri: research project execution; review and critique of the manuscript. Giulia Giannini: research project execution research project organization and execution; review and critique of the manuscript. Roberta Pantieri: research project execution; review and critique of the manuscript. Giuseppe Samoggia: research project execution; review and critique of the manuscript. Cesa Scaglione: research project organization; review and critique of the manuscript. Susanna Trombetti: research project organization Alessandro: research project conception, organization and execution; review and critique of the manuscript. Francesco Nonino: research project conception None.