LitCovid-sentences

PMC:7730445 / 13888-13892 JSON TXT

Death by SARS-CoV 2: a Romanian COVID-19 multi-centre comorbidity study Abstract Evidence regarding the relation between SARS-CoV-2 mortality and the underlying medical condition is scarce. We conducted an observational, retrospective study based on Romanian official data about location, age, gender and comorbidities for COVID-19 fatalities. Our findings indicate that males, hypertension, diabetes, obesity and chronic kidney disease were most frequent in the COVID-19 fatalities, that the burden of disease was low, and that the prognosis for 1-year survival probability was high in the sample. Evidence shows that age-dependent pairs of comorbidities could be a negative prognosis factor for the severity of disease for the SARS-CoV 2 infection. Introduction The rapid spread of the new type of coronavirus named novel Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-COV-2) from Wuhan (China) in the world determined the World Health Organization (WHO) to declare novel coronavirus disease (2019 COVID-19 or nCOVID-19) a pandemic on March 11, 20201,2. The first confirmed case in Romania was on February 26, 2020, and, by April 20 2020, 8936 confirmed case, 2017 healed cases and 520 deaths had been recorded, 88 cases with missing data and 98,491 tests were done3. The case definitions, as well as the therapies applied by the Romanian physicians, were also used in Italy, Spain and France4–7. The diagnostic was established based on the clinical symptomatology like cough, temperature, breathing difficulty without a prescribed etiology, on patient's direct exposure to SARS-COV-2, and the corroboration by testing with the RT-PCR method8,9. The Romanian government took strict measures to limit the outbreak, and a Coronavirus task force was set up to coordinate this effort under the direct supervision of the Ministry of Health (MS). Even though some cases of intra-hospital or community infection arose, they were isolated and contained with the help of special hygiene measures, social distancing and a very tight lockdown imposed by the authorities. However, a large number of deaths were confirmed among the adult population with no patient under 20 years of age dying because of COVID-19. Hence, issues regarding the COVID-19 mortality, especially the relationship with the comorbidities and the burden of disease are of great importance for the medical community8. Several studies about the comorbidities of the patients infected with SARS-COV-2 indicate hypertension, diabetes, obesity, neoplasms, chronic kidney disease and chronic obstructive pulmonary disease (COPD) as leading risk factors for a lethal evolution of the disease10–12. Romania holds a top place in Europe for deaths caused by cardiovascular diseases with hypertension affecting over 45% of the adult population. Almost half of the population aged over 65 has at least one chronic disease, and 1 in 5 adults smokes daily (32% men as compared to 18% women)13. Moreover, Romania has a record number of deaths due to preventable causes like ischemic cardiac disease, pulmonary cancer or alcoholism14. The prevalence of diabetes is approximately 11.6% and doubles for prediabetes as confirmed by Predatorr study15 and Mentor study16. Therefore, it is expected that a high number of infections with SARS-CoV-2 in Romania to have a severe course of the disease or to end deadly. To investigate the comorbidity profiles for SARS-CoV2 deaths, an observational retrospective study was conducted starting from the official public communications of Romanian Ministry of Internal Affairs (MAI) regarding the fatalities declared as COVID-19 deaths by the National Institute of Public Health (INSP)4,8. Additionally, data from the National Center for Statistics and Informatics in Public Health were utilized to run a case-case study between the COVID-19 death population and the fatalities due to hospital pneumonia between March 22 and April 20 in the years 2016–2018 in Romania. As far as we know, this is the first study in Eastern Europe that investigates patterns of comorbidities for SARS-COV-2 fatalities. Our principal research question is to generate hypothesis about comorbidity and burden of disease in the population deceased because of the COVID-19 virus and to investigate their association with gender and age. Additionally, the comparison between the fatalities due to COVID-19 and to pneumonia in terms of gender, age and comorbidity profiles allows us to characterize the specific medical pre-conditions of COVID-19 population. Results Population characteristics Figure 1 illustrates the history of the SARS-CoV-2 pandemic in Romania from February 26 until April 20, including dates, number of fatalities and number of the missing subjects due to incomplete information. The Romanian COVID-19 mortality study has a sample size of 432 patients with complete data regarding gender, age and comorbidities. It consists of a multi-case, multi-center approach focused on the mortality frequency and of a case-case evaluation with the aim of comparing hospital-pneumonia and COVID-19 mortality (see the section on the study design for the details). Our goal is to investigate the association between the mortality and comorbidities, gender, age and hospital-pneumonia based on the statistical methods presented in detail in the last section of the present article. Since the Romanian population mainly consists of 84.1% Ethnic Romanians and of a small Hungarian minority that makes up 5.4% of the population demographics that both belong to the white ethnicity, an effect of the ethnicity was not considered in this article17. Figure 1 History of SARS-COV-2 pandemic in Romania up to April 20. First, the basic descriptive statistics about our sample are presented. There are 282 (65.3%) men and 150 (34.7%) women with a statistically significant difference between the size of these two groups (Table S1). The general mean age is 67 years, with an SD of 13.1, and the median value is 68 (Table S2). Sources of infection with deaths due to COVID-19 virus were reported in 36 counties in Romania with significant case concentration in Suceava-county n. = 67 (15.51%), Arges-county n. = 44 (10.19%) and Hunedoara-county n. = 41 (9.49%). Details regarding the spatial distribution of the distinct diseases and aggregated groups of diseases are displayed in Table S6 and Table S7 in the supplementary material. A total of 170 different conditions and 20 groups of diseases were identified in the study, with a broader diversity of 149 conditions for men over 107 conditions for women. The diseases with the highest prevalence in the sample are hypertension, obesity, diabetes and chronic kidney disease as well as diseases of the circulatory system and nutritional or metabolic disorders. The sample size of the pneumonia study is equal to 874 persons with complete data as to gender, age and comorbidities. The sample includes 492 (56.3%) men and 382 (43.7%) women, the average age is 73.3 years, with an SD of 13.5, and the median age is 76. Association between gender or age and comorbidities One-sample proportion test indicates a statistically significant difference in the relative frequency of man and women in the sample and in the subpopulations with hypertension, diabetes, diabetes mellitus type 2, chronic kidney diseases, dependence on renal dialysis, chronic obstructive pulmonary disease, acute ischaemic heart disease or respiratory failure as displayed in Table S1. This is also valid for the subpopulations without comorbidities and for those corresponding to the aggregated group of diseases with a relative frequency of more than 5% except for diseases of the nervous system. On the other hand, the variable gender is not significantly associated with single conditions or group of illnesses except for heart failure (see Table 1). Table 1 Association between gender and distinct diseases or aggregated group of diseases, as absolute numbers and as percentage. Characteristicsn,% Menn = 282 Womenn = 150 p-value Distinct diseases Hypertension 113 40.07% 49 32.67% 0.16 Obesity 32 11.35% 21 14.00% 0.52 Diabetes Total 105 37.23% 48 32% 0.33 Diabetes mellitus, unspecified 40 14.18% 26 17% 0.47 Diabetes mellitus type 2 56 19.86% 19 13% 0.081 Diabetes mellitus type 1 9 3.19% 3 2% 0.68 Diseases of the circulatory system, unspecified 41 14.54% 27 18.00% 0.42 Supraventricular tachyarrhythmia 16 5.67% 9 6.00% 1 Heart failure, unspecified 14 4.96% 16 10.67% 0.043 Cerebral ischaemic stroke 16 5.67% 12 8% 0.47 Chronic kidney disease 31 10.99% 13 8.67% 0.55 Kidney failure, unspecified 17 6.03% 9 6.00% 1 Dependence on renal dialysis 29 10.28% 14 9.33% 0.88 Chronic obstructive pulmonary disease 18 6.38% 7 4.67% 0.61 Acute ischaemic heart disease, unspecified 15 5.32% 4 2.67% 0.3 Respiratory failure 15 6.32% 4 2.67% 0.3 Aggregated group of diseases No comorbidities 11 3.90% 2 1.33% 0.23 Diseases of the circulatory system 168 59.57% 88 58.67% 0.94 Endocrine, nutritional or metabolic diseases 122 43.26% 63 42.00% 0.88 Diseases of the genitourinary system 54 19.15% 25 16.67% 0.61 Diseases of the respiratory system 47 16.67% 20 13.33% 0.44 Diseases of the nervous system 35 12.41% 26 17.33% 0.21 Diseases of the digestive system 30 10.64% 14 9.33% 0.79 Neoplasms 26 9.22% 14 9.33% 1 Mental, behavioural or neurodevelopmental disorders 25 8.87% 11 7.33% 0.71 Factors influencing health status or contact with health services 37 13.12% 17 11.33% 0.7 Symptoms, signs or clinical findings, not elsewhere classified 13 4.61% 11 7.33% 0.34 Certain infectious or parasitic diseases 15 5.32% 5 3.33% 0.49 In bold the differences that reached the statistical significance. Because age is an established risk factor for mortality, we also performed an age-stratified analysis with the lifespan split in 10-year intervals starting with age 50. The chi-squared test of association between the gender and age factors is statistically significant, with the largest male relative frequency (81.16%) between 50 and 59 years. Moreover, the factor age is also significantly associated with most of the comorbidities, with hypertension and dependence on renal dialysis following also this trend but at a significance level of 10% (see Table 2). Table 2 Association between age and distinct diseases or group of diseases, as absolute numbers and as percentage. Characteristicsn, % < 504410.19% 50–596915.97% 60–6912027.78% 70–7912829.93% ≥ 807116.44% p-value Gender Male 29 66% 56 81.16% 73 60.83% 80 62.5% 44 61.97% 0.05 Female 15 34.09% 13 19.84% 47 39.17% 48 37.5% 27 38.03% Distinct diseases Hypertension 12 27.28% 21 30.43% 49 40.83% 45 35.16% 35 49.30% 0.074 Obesity 13 29.55% 9 13.04% 20 16.67% 8 6.25% 3 4.23% 0.0001 Diabetes Total 14 31.82% 23 33.33% 53 44.17% 48 37.5% 16 22.53% 0.045 Diabetes mellitus, unspecified 7 15.91% 8 11.59% 18 15.00% 24 18.75% 10 14.08% 0.74 Diabetes mellitus type 2 7 15.91% 12 17.39% 31 25.83% 20 15.63% 5 7.04% 0.021 Diabetes mellitus type 1 – 3 4.34% 4 3.33% 4 3.13% 1 1.40% – Diseases of the circulatory system, unspecified 6 13.64% 7 10.14% 18 15.00% 23 17.97% 14 19.72% 0.53 Supraventricular tachyarrhythmia 1 2.27% – – 12 9.38% 10 14.08% – Heart failure, unspecified – 5 7.25% 10 8.33% ( 7 5.47% 8 11.27% – Cerebral ischaemic stroke 2 4.55% 1 1.49% 9 7.50% 9 7.03% 7 9.86% 0.31 Acute ischaemic heart disease, unspecified – 2 2.90% 7 5.83% 4 3.13% 6 8.45% – Chronic kidney disease 4 9.09% 3 4.35% 14 11.67% 14 10.94% 9 12.67% 0.48 Dependence on renal dialysis 3 6.82% 13 13.04% 13 10.83% 15 11.72% 3 4.23% 0.072 Kidney failure, unspecified 3 6.82% 3 4.35% 9 7.50% 8 6.25% 3 4.23% 0.87 Respiratory failure 2 4.55% 4 5.80% ) 3 2.50% 4 3.13% 6 8.45% 0.33 Chronic obstructive pulmonary disease 1 2.27% – 12 10.00% 10 7.81% 2 2.82% – Aggregated group of diseases Diseases of the circulatory system 14 31.82% 33 47.83% 74 61.67% 82 64.06% 53 74.65% 0.0001 Endocrine, nutritional or metabolic diseases 21 47.73% 30 43.48% 63 52.50% 53 41.41% 18 25.35% 0.007 Diseases of the genitourinary system 7 15.91% 7 10.14% 24 20.00% 29 22.66% 12 16.90% 0.27 Diseases of the respiratory system 6 13.64% 10 14.49% 18 15.00% 20 15.63% 13 18.31% 0.96 Diseases of the nervous system 3 6.82% 2 2.90% 18 15.00% 19 14.84% 19 26.76% 0.001 Diseases of the digestive system 4 9.09% 8 11.59% 18 15.00% 10 7.81% 4 5.63% 0.23 Neoplasms 1 2.27% 7 10.14% 12 10.00% 16 12.50% 4 5.63% 0.25 Mental, behavioural or neurodevelopmental disorders 3 6.82% 2 2.90% 6 5.00% 15 11.72% 10 14.08% 0.049 Factors influencing health status or contact with health services 4 9.09% 10 14.49% 15 12.50% 19 14.84% 6 8.45% 0.66 Symptoms, signs or clinical findings, not elsewhere classified 2 4.55% 3 4.35% 6 5.00% 11 8.59% 2 2.82% 0.46 Certain infectious or parasitic diseases 3 6.82% 2 2.90% 6 5.00% 4 3.13% 5 7.04% 0.63 In bold the p-values that reached the statistical significance. Comparison between COVID-19 and pneumonia fatalities The comparison with pneumonia fatalities yielded that the proportion of males in the COVID-19 death sample is also significantly higher than in the pneumonia death sample as shown in Table 3. Age and the relative frequency of essential hypertension, diseases of the circulatory system, obesity, dependence of renal dialysis, and chronic kidney failure are also different between the COVID-19 and pneumonia groups. Additionally, the odd to die for COVID-19 versus hospital pneumonia is 150% greater for men than women with similar comorbidities (see Table 4). Moreover, the adjusted odds to die for COVID-19 compared to pneumonia are 22 times higher for persons with the chronic kidney failure, doubled for persons with hypertension, and almost equal for patients with type 2 diabetes. Table 3 Comparison between the characteristics of the COVID-19 and pneumonia groups. Characteristics COVID-19 Pneumonia p-value Sample size n, % 432 100% 711 100% Gender n, % 0.0001 Female 150 35% 320 45% Male 282 65% 391 54.99% Age Mean (sd) Median 66.97 (13.07) 68 73.14 (13.68) 78 Age group n, % 0.0001 < 50 44 10.19% 49 6.89% 50–59 69 15.97% 51 7.17% 60–69 120 27.78% 141 19.83% 70–79 128 29.63% 199 27.99% ≥ 80 71 16.44% 271 38.12% Comorbidities Hypertension 162 37.50% 181 25.46% 0.0001 Obesity 56 12.96% 63 8.86% 0.03 Diabetes mellitus type I 12 2.78% 5 0.70% 0.01 Diabetes mellitus type II 75 17.36% 145 20.39% 0.21 Dialysis 43 9.95% 4 0.56% 0.0001 Chronic kidney disease 43 9.95% 88 12.38% 0.21 Supraventricular tachyarrhythmia 25 5.79% 201 28.27% 0.0001 Congestive heart failure 12 2.78% 204 28.69% 0.0001 COPD 24 5.56% 72 10.13% 0.01 Heart failure 29 6.71% 27 3.80% 0.03 Kidney failure 26 6.02% 4 0.56% 0.0001 Respiratory failure 19 4.40% 88 12.38% 0.0001 Cerebral ischaemic stroke 28 6.48% – 0.0001 In bold the p-values that reached the statistical significance. Table 4 Logistic regression analysis for the comparison between COVID-19 and pneumonia death groups. Covariates CoefficientΒ OR (COVID-19 vs. pneumonia) 95% CI p-value Age − 0.022 0.978 [0.968–0.988] 0.0001 Gender (0 = female) 0.414 1.513 [1.145–2.000] 0.01 Hypertension 0.740 2.096 [1.560–2.818] 0.0001 Obesity 0.266 1.305 [0.843–2.019] 0.23 Diabetes type 2 − 0.355 0.701 [0.494–0.996] 0.049 Dialysis 3.104 22.290 [7.244–68.590] 0.0001 Chronic kidney disease − 0.522 0.593 [0.369–0.955] 0.03 Heart failure 1.095 2.990 [1.612–5.546] 0.0001 COPD − 0.730 0.482 [0.284–0.817] 0.01 Supraventricular tachyarrhythmia − 1.687 0.185 [0.116–0.295] 0.0001 Respiratory failure − 1.045 0.352 [0.203–0.609] 0.0001 Intercept 0.884 2.420 0.02 In bold the p-values that reached the statistical significance. The burden of disease and survival probability We also extended the analysis to the coexistence of multiple health conditions, and we focused on the multimorbidity. This has a mean value of 2.109, an SD of 1.66 and a median value of 2 in our study. The multimorbidity factor with levels 0, 1, 2, 3, 4 or more than 5 for single or aggregated comorbidities is significantly associated with the factor age but not with gender, confirming the pattern displayed at the level of single conditions and group of diseases (Tables S3 and S4). The Charlson comorbidity index (CCI) had a mean of 1.324, SD of 0.95 and median of one indicating that the severity of the underlying medical condition in the study was mainly mild. This was confirmed by the high relative frequency of 91.20% for the mild level for the factor CCI. There were no significant differences between males and females or among age groups with respect to CCI (Tables S3 and S4). Moreover, 1-year survival probability averaged out at 81%, and 50% of the study individuals had a prognosis for it of 85%. Co-occurrence patterns The co-occurrence of the diseases is of great importance since the mean multimorbidity of 2.109 indicates that not single but pair of diseases are pre-conditions in our mortality study. Pearson's phi coefficient is a robust measure of the comorbidity association, and its values for diseases in general and separated for men and women are displayed in the Table S5. Different co-occurrence patterns emerge between men and women with the striking lack of the correlation between pneumonia or respiratory failure, respectively stroke and dissociative neurological symptom disorder in the female profile. The largest value of the parameter phi of 0.47 respectively 0.40 for female and 0.59 for male corresponds to kidney failure and dependence to dialysis, while the latter is also correlated with chronic kidney disease as the correlation lies by 0.30 in general and 0.38 in men. They are followed by pneumonia or respiratory failure with phi equal to 0.45 in general and to 0.55 in male and by stroke and dissociative neurological symptom disorder with a correlation coefficient of 0.41 in the whole study and of 0.49 in men. Diabetes mellitus type 2 is correlated with hypertension not only for the whole study (phi = 0.11) but also in the men (phi = 0.10) and women subsamples (phi = 0.11). Hypertension appears mainly together with stroke with phi equal to 0.15, 0.16 respectively 0.14 for the whole study, men and women. Finally, a comorbidity network analysis rounded off our exploratory investigation. This identified 11 clusters in the general study, 9 clusters in male, and 6 clusters in female subgroups, each including two or three diseases reinforcing the previous results about pair of diseases as pre-condition for COVID-19 fatalities (see Fig. 2). Figure 2 Co-occurrence clusters for All Subjects (A), Women (B) and Men (C). Discussion The high relative frequency of men in the whole study (65.3% males vs. 34.7% females, p-value < 0.001) and in all subpopulations with distinct diseases and aggregated group of diseases shows that an important risk factor for death by COVID-19 in Romania is being male. This is consistent with other studies about predictors for severity and mortality due to infection with SARS-CoV 2 in other populations10,18. In addition, the comorbidities with the highest prevalence in our study, such as hypertension (37.5%), diabetes (35.4%), obesity (12.27%) and chronic kidney disease (10.19%) as well as diseases of the circulatory system (59.26%) and nutritional or metabolic disorders (42.82%), had similar frequencies with those reported in Chinese or American studies but contradicts the findings of the Italian CORIST study10,19–21. Indeed, the percentages that we found for the distinct comorbidities are consistent to those reported for China; for instance. Li et al., found that the majority of patients deceased in the hospital (N = 128) were men (73%), and that the most prevalent comorbidities associated with COVID-19 mortality were hypertension (47%) and diabetes (13%), while Zhou et al. reported that, among 54 fatalities for COVID-19, 62% were men, and the most prevalent comorbidities were hypertension (30%), diabetes (19%) and coronary heart disease (8%)22,23. These diseases also display significant differences in the relative frequency of men vs. women in our study in Romanian population: 74.67% men vs. 25.33% women for diabetes mellitus type 2, 69.75% men vs. 30.25% women for hypertension, 70.45% men vs. 29.55% women for chronic kidney disease, as well as 65.62% vs. 34.38% respectively 65.95% vs. 34.95% for diseases of the circulatory system and for nutritional or metabolic disorders. Almost all the comorbidities in our study are not associated with gender but with age, as also previously reported in other populations24,25. We found in the present study the highest number of deaths in the age group between 70–79 years (N = 128, 29.63%), followed by the age group 60–69 years (N = 120, 27.78%) and the majority were men (62.5% men vs. 37.5% women). The Romanian age-specific mortality statistics are similar to the estimates of crude-case fatality from Italy, Spain or Germany at the beginning of April, and in contrast to the United Kingdom, where the death count was higher in 65+ years population cluster with a relative frequency of 44% for persons aged 65–79 years and 46% for the 80+ age group4–8,26. The Romanian age structure of the COVID-10 fatalities is similar to the Chinese and lower than the Italian group9,25. Subjects younger than 50 years are less susceptible to die by COVID-19 with a relative frequency of 10% as reported in other articles, and those over 70 years are almost half of the study participants (46%)27. In the survey by Du et al. it was found that subjects who died for SARS-CoV-2 infection were older than 65 years, and that the main comorbidities related to death were hypertension and CVD28. On the other side, in ten European countries and Canada, fatalities due to COVID-19 for individuals younger than 65 years represented 4.8–9.3% of all. Only 13.0% of the COVID-19 fatalities in the UK and 7.8–23.9% in the USA were younger than 65 years. On the other hand, the prevalence for subjects older than 80 years was 54–69% in Europe, 67% in Canada and 36–63% in US27. Another important meta-analysis made by Verity et al. has shown that 13.4% persons older than 80 with more than one comorbidity died for SARS-CoV-2 and that the estimated case fatality rate was 1.38% for 38 countries29. Also, in the Spanish study by Perez-Tanoira et al., men older than 85 years represented 27% of their cohort (N = 108) and the most prevalent comorbidity was coronary heart disease, with a crude mortality rate of 19%30. Statistics displayed in Table S1 and Table 3 show that gender and age are decisive factors for COVID-19 mortality in our Romanian population. This is of great importance for tailored personalized medicine, where treatment management needs to be adapted to patients’ characteristics, and in the epidemiology, where population stratification is the key factor for the identification of the highest risk subgroups. In the whole world, the incidence of COVID-19 and its mortality rate is higher for people with diabetes and/or obesity. The precise mechanisms of SARS-CoV-2 infection are still not fully elucidated, but it has been suggested that this may be due to the chronic inflammatory state and/or insulin-resistance31,32. In our study, diabetes mellitus was classified as type-1, type-2 and unspecified. The unspecified diabetes is most probably the form of type-2 diabetes, since most of the deceased people are between 50 and 79 years and the incidence of the diabetes of type 2 in Romania is extremely higher than type-1, as also shown in the Predatorr study15. Hence, we can assume that the relative frequency of diabetes type-2 was higher than hypertension for the three clusters of age 50–59, 60–69 and 70–79 years. The association between age and different comorbidities suggest the existence of age-stratified positive and negative markers for COVID-19 mortality. Generally, obesity correlates positively with diabetes30,31, although this is confirmed in our study only in the group of youngest age. with a prevalence of 29.55%. Hence, we can interpret obesity as a negative marker for the severity of COVID-19 infection in combination with age below 50 years. In addition, we found the following negative markers: for the age group 50–59 years diabetes (33.33%), for the age group 60–69 years diabetes (44.17%) and hypertension (35.16%), for the age group 70–79 diabetes (44.17%), hypertension (35.16%) and chronic kidney diseases (17.87%), and for the age group over 80 years hypertension (49%), diabetes (22.53%) and the diseases of the circulatory system (19.72%). Further, the main risk factors for death regarding aggregated groups of diseases were the following: the endocrine diseases (diabetes, obesity) with a prevalence of 47.73% for subjects under 50 years, and diseases of the circulatory system, with a prevalence between 47.83 and 74.65% for subjects in all the other age groups. The prevalence of respiratory failure and chronic pulmonary diseases in our study was almost the lowest compared to the rest of the comorbidities but it was similar with the data reported elsewhere22,23. Diseases of the respiratory system were more frequent in the age group over 80 years (18.31%) and had higher prevalence in men (16.67%) versus women (13.33%). This is somewhat surprising since pneumonia had a higher prevalence in the Romanian population in the last three years. Usually, patients with comorbidities from this group of diseases are more vulnerable to viral respiratory infections33. Moreover, it is known that a symptom for the infection with SARS-CoV 2 is an atypical pneumonia with a variable degree of severity34. Since pneumonia could be a negative prognostic factor for COVID-19 patients35, we investigated the relation between the comorbidities of the two diseases. The results show significant differences between these two diseases as shown in Table 4. The hypothesis about the diseases of the respiratory system and the neoplasms as prognosis factors for the mortality through COVID-19 is not confirmed in our study. Therefore, it is likely that patients died because of the impact and virulence of the SARS-CoV-2 infection over comorbidities like hypertension, diabetes mellitus or obesity, in combination with age and gender36. Surprisingly, the average multimorbidity for the COVID-19 fatalities is low but confirms previous findings about infected populations or groups with a severe form of the disease19. The high relative frequency of the mild CCI corroborates this conclusion regarding subjects with a low burden of disease. Our prognosis regarding the survival probability confirms other observations about the estimated year's life lost by the COVID-19 fatalities being "considerably more than the "1–2 years" as presented in37. This is the first article where the co-occurrence profiles are investigated with the help of Pearson's co-occurrence coefficient and where co-occurrence patterns are explicitly described through clustering. Our findings on clustering data are not surprising since diabetes appears mainly with cardiovascular diseases, and kidney diseases with dialysis. Also, the sparsity of respiratory diseases, particularly in the women subgroup, should be further investigated. In conclusion, our study characterizes comorbidities of COVID-19 fatalities in Romania and their association with gender and age. Most of our results are consistent with previous studies from other populations, and the main differences in our Romanian population lie in the lower age of the deceased patients, on their small multimorbidity count and their high 1-year survival probability. In addition, the comparison with the pneumonia group offers a new perspective on the medical pre-conditions of SARS-CoV 2. Fortunately, we registered a low incidence of COVID-19 deaths in Romania until 20th of April 2020, and future studies are needed to determine precisely the pathogenic mechanism that affected these patients, especially males with hypertension, diabetes, chronic kidney disease or obesity and potential protective conditions like those suggested in38–41. Methods Study design The cohort The Romanian COVID-19 mortality sub-study was conducted between March 22 and April 20, 2020, with the aim to explore the comorbidity patterns in the patients who died because of this new virus. A multiple-case, multiple-center design was chosen to answer this research question with data collected from the official public notifications of MAI concerning the COVID-19 deaths registered by the INSP3,4,8. Initially, 451 COVID-19 deaths were reported by INSP from the beginning of the epidemic until April 20 (Fig. 1). The list was completed afterwards with 69 additional cases due to the delay in the reporting process from some centers of COVID-19 infection. The official report included information about gender, age, county of residence, history of hospitalization and diagnostics, and/ or comorbidities. Only cases with information about gender, age and comorbidities were included in the study due to the main focus of the research42. Eventually, 432 patients who fulfilled these criteria were kept in the Romanian COVID-19 mortality sub-study. Since we aim at an analytic and not statistical generalization of the results regarding the comorbidity patterns, statistical testing was just serving exploratory purposes42. Hence, the complete-case analysis is justified in this case both by purpose and by the state of things. Case-case study Additionally, we run a case-case study utilizing data regarding the deaths due to hospital pneumonia between March 22 and April 20 in the years 2016–2018 extracted from the administrative database of the INSP. The study was carried out with the approval of the Local Bioethics Council of the INSP (number 6342/04 MAY 2020) and with the permission of the Local Bioethics Council of the “Carol Davila” University of Medicine, Bucharest Romania (9984/11 MAY 2020). Statistical analysis An exploratory data analysis was applied to identify general, spatial-, gender—or age-specific comorbidity patterns in our study. In a first step, the various conditions were sorted according to the standard diagnosis tool ICD-11 (The International Classification of the Diseases) as it was recommended for the mortality and morbidity statistics by the WHO43. The ICD-11 taxonomy includes categories of single diseases or the aggregated groups of diseases. We investigated the comorbidities at both levels of aggregation because the study included a large number of distinct conditions with low prevalence. The statistical analysis was done with the help of the statistical software R version 3.6.2. Continuous variables were expressed as mean, standard deviation (SD) and median. We focused in this paper on comparing the individual crude relative frequencies of the single and aggregated comorbidities between different gender and age groups. Since age is a significant risk factor for mortality, we investigated the frequency of the comorbidities stratified on ten-year intervals starting with the 50+ age groups. Categorical variables were tabulated in frequency (n) and relative frequency (%). A binomial and a chi-squared test for one proportion (only the last was reported) and a chi-squared test for equal proportions were applied when it was appropriate at the 5% level of significance. The post-hoc tests were not corrected for multiple comparisons since the focus of our study was on generating hypotheses for single diseases or group of diseases and not on an experiment-wise approach. The spatial distribution of the COVID-19 mortality corresponded to the country-specific administrative units– counties and was illustrated by single or aggregated disease frequency and relative frequency (%). The association between different comorbidities and the odds ratio between COVID-19 and pneumonia deaths was investigated with the help of the logistic regression, and the results were reported as coefficients β, as point estimates for odds ratio and their 95% confidence intervals. Multimorbidity was computed for each person in the sample since this is a good indicator of the burden of disease44,45 and transformed in a multimorbidity factor with levels 0, 1, 2, 3, 4, or more the 5 for single or aggregated comorbidities. A controversial issue regarding the epidemic of COVID-19 is the burden of disease in the deceased population and the hypothetical survival probability in the case that these people would not have been infected with this virus. A simplified version of the Charlson Comorbidity Index (CCI) was calculated on a scale from 0 to 8, as indicated in24 for an assessment of the individual disease severity. We introduced the CCI factor with two levels (mild for CCI < 3 and severe for CCI ≥ 3), and we tested its differences between gender and age groups by a chi-squared test. Following the same paper24, we adapted their statistical model for the individual prognosis of the death risk based on the reported hazard ratio for the age levels and for the CCI in multiple European studies as explained in the section 8 of the supplementary material. This gave us a prognosis for the mean, the standard deviation and the median of the crude 1-year death risk and survival probability. The next step in the analysis of the COVID-19 comorbidities concerns the co-occurrence of the diseases. This was explored with the help of the Pearson's phi coefficient as defined and investigated, e.g. in46 and computed for each pair of single or aggregated diseases with an incidence of at least 4. A value of 1 means a "perfect match", 0 means no co-occurrence and negative values up to − 1 means "inverse occurrence". A comorbidity network approach was considered to investigate the relationships between diseases. We computed the dissimilarity between the conditions studied in the previous step as one minus the Pearson's phi. Hence, a value of 1 corresponds to "no association" and values larger than 1 to "inverse association". The hierarchical classification was performed through an agglomerative clustering algorithm with the Ward method since this proved to have the largest agglomeration coefficient. Furthermore, the dendrogram was cut at nodes corresponding to branch length of 1 based both on a visual inspection and on predefined length of the branch. This study has several limitations. We considered only the deaths up to April 20 when the peak of the first pandemic wave was reached in Romania. Moreover, total case design was used because of the missing data in official public communications. Disclaimer This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors. Ethical compliance This research was carried out in accordance with the relevant guidelines and regulations. Informed consent We have obtained informed consent from all participants and/or their legal guardians. Supplementary Information Supplementary Information. Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. These authors contributed equally: Anca Pantea Stoian, Mihaela Pricop-Jeckstadt, Adrian Pana, Ruxandra Schitea, Antonio Ceriello and Manfredi Rizzo. Supplementary Information The online version contains supplementary material available at 10.1038/s41598-020-78575-w. Author contributions A.P.S., M.P.J., A.P. conceived the research question. A.P.S., M.P.J., A.P., B.V.I., M.R., A.C., R.S. and V.J. designed the study and analysis plan. A.P.S., M.P.J., A.P., R.S., M.G., M.R., A.C., S.P., C.N., S.P., R.H., A.I.S., D.C., L.M. and C.P. prepared the data. M.P.J. did the statistical analysis. A.P.S., M.P.J., A.P., B.T., A.T., C.S., A.A., A.C., M.R., and V.J. drafted the initial and final versions of the manuscript. All authors critically reviewed early and approved final versions of the manuscript. Competing interests APS, MR, AC have given talks, attended conferences, and participated in advisory boards and clinical trials sponsored by various pharmaceutical companies; yet, no financial or professional help was received for the preparation of this manuscript. APS is currently Vice-President, National Diabetes Commission, Ministry of Health, Romania. MR is currently Chief Medical and Scientific Advisor, Diabetes, Novo Nordisk South East Europe, Middle East and Africa (SEEMEA). The rest of the authors do not have a conflict of interests to declare for this article.

Document structure show

article-title	Death by SARS-CoV 2: a Romanian COVID-19 multi-centre comorbidity study
abstract	Evidence regarding the relation between SARS-CoV-2 mortality and the underlying medical condition is scarce. We conducted an observational, retrospective study based on Romanian official data about location, age, gender and comorbidities for COVID-19 fatalities. Our findings indicate that males, hypertension, diabetes, obesity and chronic kidney disease were most frequent in the COVID-19 fatalities, that the burden of disease was low, and that the prognosis for 1-year survival probability was high in the sample. Evidence shows that age-dependent pairs of comorbidities could be a negative prognosis factor for the severity of disease for the SARS-CoV 2 infection.
p	Evidence regarding the relation between SARS-CoV-2 mortality and the underlying medical condition is scarce. We conducted an observational, retrospective study based on Romanian official data about location, age, gender and comorbidities for COVID-19 fatalities. Our findings indicate that males, hypertension, diabetes, obesity and chronic kidney disease were most frequent in the COVID-19 fatalities, that the burden of disease was low, and that the prognosis for 1-year survival probability was high in the sample. Evidence shows that age-dependent pairs of comorbidities could be a negative prognosis factor for the severity of disease for the SARS-CoV 2 infection.
body	Introduction The rapid spread of the new type of coronavirus named novel Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-COV-2) from Wuhan (China) in the world determined the World Health Organization (WHO) to declare novel coronavirus disease (2019 COVID-19 or nCOVID-19) a pandemic on March 11, 20201,2. The first confirmed case in Romania was on February 26, 2020, and, by April 20 2020, 8936 confirmed case, 2017 healed cases and 520 deaths had been recorded, 88 cases with missing data and 98,491 tests were done3. The case definitions, as well as the therapies applied by the Romanian physicians, were also used in Italy, Spain and France4–7. The diagnostic was established based on the clinical symptomatology like cough, temperature, breathing difficulty without a prescribed etiology, on patient's direct exposure to SARS-COV-2, and the corroboration by testing with the RT-PCR method8,9. The Romanian government took strict measures to limit the outbreak, and a Coronavirus task force was set up to coordinate this effort under the direct supervision of the Ministry of Health (MS). Even though some cases of intra-hospital or community infection arose, they were isolated and contained with the help of special hygiene measures, social distancing and a very tight lockdown imposed by the authorities. However, a large number of deaths were confirmed among the adult population with no patient under 20 years of age dying because of COVID-19. Hence, issues regarding the COVID-19 mortality, especially the relationship with the comorbidities and the burden of disease are of great importance for the medical community8. Several studies about the comorbidities of the patients infected with SARS-COV-2 indicate hypertension, diabetes, obesity, neoplasms, chronic kidney disease and chronic obstructive pulmonary disease (COPD) as leading risk factors for a lethal evolution of the disease10–12. Romania holds a top place in Europe for deaths caused by cardiovascular diseases with hypertension affecting over 45% of the adult population. Almost half of the population aged over 65 has at least one chronic disease, and 1 in 5 adults smokes daily (32% men as compared to 18% women)13. Moreover, Romania has a record number of deaths due to preventable causes like ischemic cardiac disease, pulmonary cancer or alcoholism14. The prevalence of diabetes is approximately 11.6% and doubles for prediabetes as confirmed by Predatorr study15 and Mentor study16. Therefore, it is expected that a high number of infections with SARS-CoV-2 in Romania to have a severe course of the disease or to end deadly. To investigate the comorbidity profiles for SARS-CoV2 deaths, an observational retrospective study was conducted starting from the official public communications of Romanian Ministry of Internal Affairs (MAI) regarding the fatalities declared as COVID-19 deaths by the National Institute of Public Health (INSP)4,8. Additionally, data from the National Center for Statistics and Informatics in Public Health were utilized to run a case-case study between the COVID-19 death population and the fatalities due to hospital pneumonia between March 22 and April 20 in the years 2016–2018 in Romania. As far as we know, this is the first study in Eastern Europe that investigates patterns of comorbidities for SARS-COV-2 fatalities. Our principal research question is to generate hypothesis about comorbidity and burden of disease in the population deceased because of the COVID-19 virus and to investigate their association with gender and age. Additionally, the comparison between the fatalities due to COVID-19 and to pneumonia in terms of gender, age and comorbidity profiles allows us to characterize the specific medical pre-conditions of COVID-19 population. Results Population characteristics Figure 1 illustrates the history of the SARS-CoV-2 pandemic in Romania from February 26 until April 20, including dates, number of fatalities and number of the missing subjects due to incomplete information. The Romanian COVID-19 mortality study has a sample size of 432 patients with complete data regarding gender, age and comorbidities. It consists of a multi-case, multi-center approach focused on the mortality frequency and of a case-case evaluation with the aim of comparing hospital-pneumonia and COVID-19 mortality (see the section on the study design for the details). Our goal is to investigate the association between the mortality and comorbidities, gender, age and hospital-pneumonia based on the statistical methods presented in detail in the last section of the present article. Since the Romanian population mainly consists of 84.1% Ethnic Romanians and of a small Hungarian minority that makes up 5.4% of the population demographics that both belong to the white ethnicity, an effect of the ethnicity was not considered in this article17. Figure 1 History of SARS-COV-2 pandemic in Romania up to April 20. First, the basic descriptive statistics about our sample are presented. There are 282 (65.3%) men and 150 (34.7%) women with a statistically significant difference between the size of these two groups (Table S1). The general mean age is 67 years, with an SD of 13.1, and the median value is 68 (Table S2). Sources of infection with deaths due to COVID-19 virus were reported in 36 counties in Romania with significant case concentration in Suceava-county n. = 67 (15.51%), Arges-county n. = 44 (10.19%) and Hunedoara-county n. = 41 (9.49%). Details regarding the spatial distribution of the distinct diseases and aggregated groups of diseases are displayed in Table S6 and Table S7 in the supplementary material. A total of 170 different conditions and 20 groups of diseases were identified in the study, with a broader diversity of 149 conditions for men over 107 conditions for women. The diseases with the highest prevalence in the sample are hypertension, obesity, diabetes and chronic kidney disease as well as diseases of the circulatory system and nutritional or metabolic disorders. The sample size of the pneumonia study is equal to 874 persons with complete data as to gender, age and comorbidities. The sample includes 492 (56.3%) men and 382 (43.7%) women, the average age is 73.3 years, with an SD of 13.5, and the median age is 76. Association between gender or age and comorbidities One-sample proportion test indicates a statistically significant difference in the relative frequency of man and women in the sample and in the subpopulations with hypertension, diabetes, diabetes mellitus type 2, chronic kidney diseases, dependence on renal dialysis, chronic obstructive pulmonary disease, acute ischaemic heart disease or respiratory failure as displayed in Table S1. This is also valid for the subpopulations without comorbidities and for those corresponding to the aggregated group of diseases with a relative frequency of more than 5% except for diseases of the nervous system. On the other hand, the variable gender is not significantly associated with single conditions or group of illnesses except for heart failure (see Table 1). Table 1 Association between gender and distinct diseases or aggregated group of diseases, as absolute numbers and as percentage. Characteristicsn,% Menn = 282 Womenn = 150 p-value Distinct diseases Hypertension 113 40.07% 49 32.67% 0.16 Obesity 32 11.35% 21 14.00% 0.52 Diabetes Total 105 37.23% 48 32% 0.33 Diabetes mellitus, unspecified 40 14.18% 26 17% 0.47 Diabetes mellitus type 2 56 19.86% 19 13% 0.081 Diabetes mellitus type 1 9 3.19% 3 2% 0.68 Diseases of the circulatory system, unspecified 41 14.54% 27 18.00% 0.42 Supraventricular tachyarrhythmia 16 5.67% 9 6.00% 1 Heart failure, unspecified 14 4.96% 16 10.67% 0.043 Cerebral ischaemic stroke 16 5.67% 12 8% 0.47 Chronic kidney disease 31 10.99% 13 8.67% 0.55 Kidney failure, unspecified 17 6.03% 9 6.00% 1 Dependence on renal dialysis 29 10.28% 14 9.33% 0.88 Chronic obstructive pulmonary disease 18 6.38% 7 4.67% 0.61 Acute ischaemic heart disease, unspecified 15 5.32% 4 2.67% 0.3 Respiratory failure 15 6.32% 4 2.67% 0.3 Aggregated group of diseases No comorbidities 11 3.90% 2 1.33% 0.23 Diseases of the circulatory system 168 59.57% 88 58.67% 0.94 Endocrine, nutritional or metabolic diseases 122 43.26% 63 42.00% 0.88 Diseases of the genitourinary system 54 19.15% 25 16.67% 0.61 Diseases of the respiratory system 47 16.67% 20 13.33% 0.44 Diseases of the nervous system 35 12.41% 26 17.33% 0.21 Diseases of the digestive system 30 10.64% 14 9.33% 0.79 Neoplasms 26 9.22% 14 9.33% 1 Mental, behavioural or neurodevelopmental disorders 25 8.87% 11 7.33% 0.71 Factors influencing health status or contact with health services 37 13.12% 17 11.33% 0.7 Symptoms, signs or clinical findings, not elsewhere classified 13 4.61% 11 7.33% 0.34 Certain infectious or parasitic diseases 15 5.32% 5 3.33% 0.49 In bold the differences that reached the statistical significance. Because age is an established risk factor for mortality, we also performed an age-stratified analysis with the lifespan split in 10-year intervals starting with age 50. The chi-squared test of association between the gender and age factors is statistically significant, with the largest male relative frequency (81.16%) between 50 and 59 years. Moreover, the factor age is also significantly associated with most of the comorbidities, with hypertension and dependence on renal dialysis following also this trend but at a significance level of 10% (see Table 2). Table 2 Association between age and distinct diseases or group of diseases, as absolute numbers and as percentage. Characteristicsn, % < 504410.19% 50–596915.97% 60–6912027.78% 70–7912829.93% ≥ 807116.44% p-value Gender Male 29 66% 56 81.16% 73 60.83% 80 62.5% 44 61.97% 0.05 Female 15 34.09% 13 19.84% 47 39.17% 48 37.5% 27 38.03% Distinct diseases Hypertension 12 27.28% 21 30.43% 49 40.83% 45 35.16% 35 49.30% 0.074 Obesity 13 29.55% 9 13.04% 20 16.67% 8 6.25% 3 4.23% 0.0001 Diabetes Total 14 31.82% 23 33.33% 53 44.17% 48 37.5% 16 22.53% 0.045 Diabetes mellitus, unspecified 7 15.91% 8 11.59% 18 15.00% 24 18.75% 10 14.08% 0.74 Diabetes mellitus type 2 7 15.91% 12 17.39% 31 25.83% 20 15.63% 5 7.04% 0.021 Diabetes mellitus type 1 – 3 4.34% 4 3.33% 4 3.13% 1 1.40% – Diseases of the circulatory system, unspecified 6 13.64% 7 10.14% 18 15.00% 23 17.97% 14 19.72% 0.53 Supraventricular tachyarrhythmia 1 2.27% – – 12 9.38% 10 14.08% – Heart failure, unspecified – 5 7.25% 10 8.33% ( 7 5.47% 8 11.27% – Cerebral ischaemic stroke 2 4.55% 1 1.49% 9 7.50% 9 7.03% 7 9.86% 0.31 Acute ischaemic heart disease, unspecified – 2 2.90% 7 5.83% 4 3.13% 6 8.45% – Chronic kidney disease 4 9.09% 3 4.35% 14 11.67% 14 10.94% 9 12.67% 0.48 Dependence on renal dialysis 3 6.82% 13 13.04% 13 10.83% 15 11.72% 3 4.23% 0.072 Kidney failure, unspecified 3 6.82% 3 4.35% 9 7.50% 8 6.25% 3 4.23% 0.87 Respiratory failure 2 4.55% 4 5.80% ) 3 2.50% 4 3.13% 6 8.45% 0.33 Chronic obstructive pulmonary disease 1 2.27% – 12 10.00% 10 7.81% 2 2.82% – Aggregated group of diseases Diseases of the circulatory system 14 31.82% 33 47.83% 74 61.67% 82 64.06% 53 74.65% 0.0001 Endocrine, nutritional or metabolic diseases 21 47.73% 30 43.48% 63 52.50% 53 41.41% 18 25.35% 0.007 Diseases of the genitourinary system 7 15.91% 7 10.14% 24 20.00% 29 22.66% 12 16.90% 0.27 Diseases of the respiratory system 6 13.64% 10 14.49% 18 15.00% 20 15.63% 13 18.31% 0.96 Diseases of the nervous system 3 6.82% 2 2.90% 18 15.00% 19 14.84% 19 26.76% 0.001 Diseases of the digestive system 4 9.09% 8 11.59% 18 15.00% 10 7.81% 4 5.63% 0.23 Neoplasms 1 2.27% 7 10.14% 12 10.00% 16 12.50% 4 5.63% 0.25 Mental, behavioural or neurodevelopmental disorders 3 6.82% 2 2.90% 6 5.00% 15 11.72% 10 14.08% 0.049 Factors influencing health status or contact with health services 4 9.09% 10 14.49% 15 12.50% 19 14.84% 6 8.45% 0.66 Symptoms, signs or clinical findings, not elsewhere classified 2 4.55% 3 4.35% 6 5.00% 11 8.59% 2 2.82% 0.46 Certain infectious or parasitic diseases 3 6.82% 2 2.90% 6 5.00% 4 3.13% 5 7.04% 0.63 In bold the p-values that reached the statistical significance. Comparison between COVID-19 and pneumonia fatalities The comparison with pneumonia fatalities yielded that the proportion of males in the COVID-19 death sample is also significantly higher than in the pneumonia death sample as shown in Table 3. Age and the relative frequency of essential hypertension, diseases of the circulatory system, obesity, dependence of renal dialysis, and chronic kidney failure are also different between the COVID-19 and pneumonia groups. Additionally, the odd to die for COVID-19 versus hospital pneumonia is 150% greater for men than women with similar comorbidities (see Table 4). Moreover, the adjusted odds to die for COVID-19 compared to pneumonia are 22 times higher for persons with the chronic kidney failure, doubled for persons with hypertension, and almost equal for patients with type 2 diabetes. Table 3 Comparison between the characteristics of the COVID-19 and pneumonia groups. Characteristics COVID-19 Pneumonia p-value Sample size n, % 432 100% 711 100% Gender n, % 0.0001 Female 150 35% 320 45% Male 282 65% 391 54.99% Age Mean (sd) Median 66.97 (13.07) 68 73.14 (13.68) 78 Age group n, % 0.0001 < 50 44 10.19% 49 6.89% 50–59 69 15.97% 51 7.17% 60–69 120 27.78% 141 19.83% 70–79 128 29.63% 199 27.99% ≥ 80 71 16.44% 271 38.12% Comorbidities Hypertension 162 37.50% 181 25.46% 0.0001 Obesity 56 12.96% 63 8.86% 0.03 Diabetes mellitus type I 12 2.78% 5 0.70% 0.01 Diabetes mellitus type II 75 17.36% 145 20.39% 0.21 Dialysis 43 9.95% 4 0.56% 0.0001 Chronic kidney disease 43 9.95% 88 12.38% 0.21 Supraventricular tachyarrhythmia 25 5.79% 201 28.27% 0.0001 Congestive heart failure 12 2.78% 204 28.69% 0.0001 COPD 24 5.56% 72 10.13% 0.01 Heart failure 29 6.71% 27 3.80% 0.03 Kidney failure 26 6.02% 4 0.56% 0.0001 Respiratory failure 19 4.40% 88 12.38% 0.0001 Cerebral ischaemic stroke 28 6.48% – 0.0001 In bold the p-values that reached the statistical significance. Table 4 Logistic regression analysis for the comparison between COVID-19 and pneumonia death groups. Covariates CoefficientΒ OR (COVID-19 vs. pneumonia) 95% CI p-value Age − 0.022 0.978 [0.968–0.988] 0.0001 Gender (0 = female) 0.414 1.513 [1.145–2.000] 0.01 Hypertension 0.740 2.096 [1.560–2.818] 0.0001 Obesity 0.266 1.305 [0.843–2.019] 0.23 Diabetes type 2 − 0.355 0.701 [0.494–0.996] 0.049 Dialysis 3.104 22.290 [7.244–68.590] 0.0001 Chronic kidney disease − 0.522 0.593 [0.369–0.955] 0.03 Heart failure 1.095 2.990 [1.612–5.546] 0.0001 COPD − 0.730 0.482 [0.284–0.817] 0.01 Supraventricular tachyarrhythmia − 1.687 0.185 [0.116–0.295] 0.0001 Respiratory failure − 1.045 0.352 [0.203–0.609] 0.0001 Intercept 0.884 2.420 0.02 In bold the p-values that reached the statistical significance. The burden of disease and survival probability We also extended the analysis to the coexistence of multiple health conditions, and we focused on the multimorbidity. This has a mean value of 2.109, an SD of 1.66 and a median value of 2 in our study. The multimorbidity factor with levels 0, 1, 2, 3, 4 or more than 5 for single or aggregated comorbidities is significantly associated with the factor age but not with gender, confirming the pattern displayed at the level of single conditions and group of diseases (Tables S3 and S4). The Charlson comorbidity index (CCI) had a mean of 1.324, SD of 0.95 and median of one indicating that the severity of the underlying medical condition in the study was mainly mild. This was confirmed by the high relative frequency of 91.20% for the mild level for the factor CCI. There were no significant differences between males and females or among age groups with respect to CCI (Tables S3 and S4). Moreover, 1-year survival probability averaged out at 81%, and 50% of the study individuals had a prognosis for it of 85%. Co-occurrence patterns The co-occurrence of the diseases is of great importance since the mean multimorbidity of 2.109 indicates that not single but pair of diseases are pre-conditions in our mortality study. Pearson's phi coefficient is a robust measure of the comorbidity association, and its values for diseases in general and separated for men and women are displayed in the Table S5. Different co-occurrence patterns emerge between men and women with the striking lack of the correlation between pneumonia or respiratory failure, respectively stroke and dissociative neurological symptom disorder in the female profile. The largest value of the parameter phi of 0.47 respectively 0.40 for female and 0.59 for male corresponds to kidney failure and dependence to dialysis, while the latter is also correlated with chronic kidney disease as the correlation lies by 0.30 in general and 0.38 in men. They are followed by pneumonia or respiratory failure with phi equal to 0.45 in general and to 0.55 in male and by stroke and dissociative neurological symptom disorder with a correlation coefficient of 0.41 in the whole study and of 0.49 in men. Diabetes mellitus type 2 is correlated with hypertension not only for the whole study (phi = 0.11) but also in the men (phi = 0.10) and women subsamples (phi = 0.11). Hypertension appears mainly together with stroke with phi equal to 0.15, 0.16 respectively 0.14 for the whole study, men and women. Finally, a comorbidity network analysis rounded off our exploratory investigation. This identified 11 clusters in the general study, 9 clusters in male, and 6 clusters in female subgroups, each including two or three diseases reinforcing the previous results about pair of diseases as pre-condition for COVID-19 fatalities (see Fig. 2). Figure 2 Co-occurrence clusters for All Subjects (A), Women (B) and Men (C). Discussion The high relative frequency of men in the whole study (65.3% males vs. 34.7% females, p-value < 0.001) and in all subpopulations with distinct diseases and aggregated group of diseases shows that an important risk factor for death by COVID-19 in Romania is being male. This is consistent with other studies about predictors for severity and mortality due to infection with SARS-CoV 2 in other populations10,18. In addition, the comorbidities with the highest prevalence in our study, such as hypertension (37.5%), diabetes (35.4%), obesity (12.27%) and chronic kidney disease (10.19%) as well as diseases of the circulatory system (59.26%) and nutritional or metabolic disorders (42.82%), had similar frequencies with those reported in Chinese or American studies but contradicts the findings of the Italian CORIST study10,19–21. Indeed, the percentages that we found for the distinct comorbidities are consistent to those reported for China; for instance. Li et al., found that the majority of patients deceased in the hospital (N = 128) were men (73%), and that the most prevalent comorbidities associated with COVID-19 mortality were hypertension (47%) and diabetes (13%), while Zhou et al. reported that, among 54 fatalities for COVID-19, 62% were men, and the most prevalent comorbidities were hypertension (30%), diabetes (19%) and coronary heart disease (8%)22,23. These diseases also display significant differences in the relative frequency of men vs. women in our study in Romanian population: 74.67% men vs. 25.33% women for diabetes mellitus type 2, 69.75% men vs. 30.25% women for hypertension, 70.45% men vs. 29.55% women for chronic kidney disease, as well as 65.62% vs. 34.38% respectively 65.95% vs. 34.95% for diseases of the circulatory system and for nutritional or metabolic disorders. Almost all the comorbidities in our study are not associated with gender but with age, as also previously reported in other populations24,25. We found in the present study the highest number of deaths in the age group between 70–79 years (N = 128, 29.63%), followed by the age group 60–69 years (N = 120, 27.78%) and the majority were men (62.5% men vs. 37.5% women). The Romanian age-specific mortality statistics are similar to the estimates of crude-case fatality from Italy, Spain or Germany at the beginning of April, and in contrast to the United Kingdom, where the death count was higher in 65+ years population cluster with a relative frequency of 44% for persons aged 65–79 years and 46% for the 80+ age group4–8,26. The Romanian age structure of the COVID-10 fatalities is similar to the Chinese and lower than the Italian group9,25. Subjects younger than 50 years are less susceptible to die by COVID-19 with a relative frequency of 10% as reported in other articles, and those over 70 years are almost half of the study participants (46%)27. In the survey by Du et al. it was found that subjects who died for SARS-CoV-2 infection were older than 65 years, and that the main comorbidities related to death were hypertension and CVD28. On the other side, in ten European countries and Canada, fatalities due to COVID-19 for individuals younger than 65 years represented 4.8–9.3% of all. Only 13.0% of the COVID-19 fatalities in the UK and 7.8–23.9% in the USA were younger than 65 years. On the other hand, the prevalence for subjects older than 80 years was 54–69% in Europe, 67% in Canada and 36–63% in US27. Another important meta-analysis made by Verity et al. has shown that 13.4% persons older than 80 with more than one comorbidity died for SARS-CoV-2 and that the estimated case fatality rate was 1.38% for 38 countries29. Also, in the Spanish study by Perez-Tanoira et al., men older than 85 years represented 27% of their cohort (N = 108) and the most prevalent comorbidity was coronary heart disease, with a crude mortality rate of 19%30. Statistics displayed in Table S1 and Table 3 show that gender and age are decisive factors for COVID-19 mortality in our Romanian population. This is of great importance for tailored personalized medicine, where treatment management needs to be adapted to patients’ characteristics, and in the epidemiology, where population stratification is the key factor for the identification of the highest risk subgroups. In the whole world, the incidence of COVID-19 and its mortality rate is higher for people with diabetes and/or obesity. The precise mechanisms of SARS-CoV-2 infection are still not fully elucidated, but it has been suggested that this may be due to the chronic inflammatory state and/or insulin-resistance31,32. In our study, diabetes mellitus was classified as type-1, type-2 and unspecified. The unspecified diabetes is most probably the form of type-2 diabetes, since most of the deceased people are between 50 and 79 years and the incidence of the diabetes of type 2 in Romania is extremely higher than type-1, as also shown in the Predatorr study15. Hence, we can assume that the relative frequency of diabetes type-2 was higher than hypertension for the three clusters of age 50–59, 60–69 and 70–79 years. The association between age and different comorbidities suggest the existence of age-stratified positive and negative markers for COVID-19 mortality. Generally, obesity correlates positively with diabetes30,31, although this is confirmed in our study only in the group of youngest age. with a prevalence of 29.55%. Hence, we can interpret obesity as a negative marker for the severity of COVID-19 infection in combination with age below 50 years. In addition, we found the following negative markers: for the age group 50–59 years diabetes (33.33%), for the age group 60–69 years diabetes (44.17%) and hypertension (35.16%), for the age group 70–79 diabetes (44.17%), hypertension (35.16%) and chronic kidney diseases (17.87%), and for the age group over 80 years hypertension (49%), diabetes (22.53%) and the diseases of the circulatory system (19.72%). Further, the main risk factors for death regarding aggregated groups of diseases were the following: the endocrine diseases (diabetes, obesity) with a prevalence of 47.73% for subjects under 50 years, and diseases of the circulatory system, with a prevalence between 47.83 and 74.65% for subjects in all the other age groups. The prevalence of respiratory failure and chronic pulmonary diseases in our study was almost the lowest compared to the rest of the comorbidities but it was similar with the data reported elsewhere22,23. Diseases of the respiratory system were more frequent in the age group over 80 years (18.31%) and had higher prevalence in men (16.67%) versus women (13.33%). This is somewhat surprising since pneumonia had a higher prevalence in the Romanian population in the last three years. Usually, patients with comorbidities from this group of diseases are more vulnerable to viral respiratory infections33. Moreover, it is known that a symptom for the infection with SARS-CoV 2 is an atypical pneumonia with a variable degree of severity34. Since pneumonia could be a negative prognostic factor for COVID-19 patients35, we investigated the relation between the comorbidities of the two diseases. The results show significant differences between these two diseases as shown in Table 4. The hypothesis about the diseases of the respiratory system and the neoplasms as prognosis factors for the mortality through COVID-19 is not confirmed in our study. Therefore, it is likely that patients died because of the impact and virulence of the SARS-CoV-2 infection over comorbidities like hypertension, diabetes mellitus or obesity, in combination with age and gender36. Surprisingly, the average multimorbidity for the COVID-19 fatalities is low but confirms previous findings about infected populations or groups with a severe form of the disease19. The high relative frequency of the mild CCI corroborates this conclusion regarding subjects with a low burden of disease. Our prognosis regarding the survival probability confirms other observations about the estimated year's life lost by the COVID-19 fatalities being "considerably more than the "1–2 years" as presented in37. This is the first article where the co-occurrence profiles are investigated with the help of Pearson's co-occurrence coefficient and where co-occurrence patterns are explicitly described through clustering. Our findings on clustering data are not surprising since diabetes appears mainly with cardiovascular diseases, and kidney diseases with dialysis. Also, the sparsity of respiratory diseases, particularly in the women subgroup, should be further investigated. In conclusion, our study characterizes comorbidities of COVID-19 fatalities in Romania and their association with gender and age. Most of our results are consistent with previous studies from other populations, and the main differences in our Romanian population lie in the lower age of the deceased patients, on their small multimorbidity count and their high 1-year survival probability. In addition, the comparison with the pneumonia group offers a new perspective on the medical pre-conditions of SARS-CoV 2. Fortunately, we registered a low incidence of COVID-19 deaths in Romania until 20th of April 2020, and future studies are needed to determine precisely the pathogenic mechanism that affected these patients, especially males with hypertension, diabetes, chronic kidney disease or obesity and potential protective conditions like those suggested in38–41. Methods Study design The cohort The Romanian COVID-19 mortality sub-study was conducted between March 22 and April 20, 2020, with the aim to explore the comorbidity patterns in the patients who died because of this new virus. A multiple-case, multiple-center design was chosen to answer this research question with data collected from the official public notifications of MAI concerning the COVID-19 deaths registered by the INSP3,4,8. Initially, 451 COVID-19 deaths were reported by INSP from the beginning of the epidemic until April 20 (Fig. 1). The list was completed afterwards with 69 additional cases due to the delay in the reporting process from some centers of COVID-19 infection. The official report included information about gender, age, county of residence, history of hospitalization and diagnostics, and/ or comorbidities. Only cases with information about gender, age and comorbidities were included in the study due to the main focus of the research42. Eventually, 432 patients who fulfilled these criteria were kept in the Romanian COVID-19 mortality sub-study. Since we aim at an analytic and not statistical generalization of the results regarding the comorbidity patterns, statistical testing was just serving exploratory purposes42. Hence, the complete-case analysis is justified in this case both by purpose and by the state of things. Case-case study Additionally, we run a case-case study utilizing data regarding the deaths due to hospital pneumonia between March 22 and April 20 in the years 2016–2018 extracted from the administrative database of the INSP. The study was carried out with the approval of the Local Bioethics Council of the INSP (number 6342/04 MAY 2020) and with the permission of the Local Bioethics Council of the “Carol Davila” University of Medicine, Bucharest Romania (9984/11 MAY 2020). Statistical analysis An exploratory data analysis was applied to identify general, spatial-, gender—or age-specific comorbidity patterns in our study. In a first step, the various conditions were sorted according to the standard diagnosis tool ICD-11 (The International Classification of the Diseases) as it was recommended for the mortality and morbidity statistics by the WHO43. The ICD-11 taxonomy includes categories of single diseases or the aggregated groups of diseases. We investigated the comorbidities at both levels of aggregation because the study included a large number of distinct conditions with low prevalence. The statistical analysis was done with the help of the statistical software R version 3.6.2. Continuous variables were expressed as mean, standard deviation (SD) and median. We focused in this paper on comparing the individual crude relative frequencies of the single and aggregated comorbidities between different gender and age groups. Since age is a significant risk factor for mortality, we investigated the frequency of the comorbidities stratified on ten-year intervals starting with the 50+ age groups. Categorical variables were tabulated in frequency (n) and relative frequency (%). A binomial and a chi-squared test for one proportion (only the last was reported) and a chi-squared test for equal proportions were applied when it was appropriate at the 5% level of significance. The post-hoc tests were not corrected for multiple comparisons since the focus of our study was on generating hypotheses for single diseases or group of diseases and not on an experiment-wise approach. The spatial distribution of the COVID-19 mortality corresponded to the country-specific administrative units– counties and was illustrated by single or aggregated disease frequency and relative frequency (%). The association between different comorbidities and the odds ratio between COVID-19 and pneumonia deaths was investigated with the help of the logistic regression, and the results were reported as coefficients β, as point estimates for odds ratio and their 95% confidence intervals. Multimorbidity was computed for each person in the sample since this is a good indicator of the burden of disease44,45 and transformed in a multimorbidity factor with levels 0, 1, 2, 3, 4, or more the 5 for single or aggregated comorbidities. A controversial issue regarding the epidemic of COVID-19 is the burden of disease in the deceased population and the hypothetical survival probability in the case that these people would not have been infected with this virus. A simplified version of the Charlson Comorbidity Index (CCI) was calculated on a scale from 0 to 8, as indicated in24 for an assessment of the individual disease severity. We introduced the CCI factor with two levels (mild for CCI < 3 and severe for CCI ≥ 3), and we tested its differences between gender and age groups by a chi-squared test. Following the same paper24, we adapted their statistical model for the individual prognosis of the death risk based on the reported hazard ratio for the age levels and for the CCI in multiple European studies as explained in the section 8 of the supplementary material. This gave us a prognosis for the mean, the standard deviation and the median of the crude 1-year death risk and survival probability. The next step in the analysis of the COVID-19 comorbidities concerns the co-occurrence of the diseases. This was explored with the help of the Pearson's phi coefficient as defined and investigated, e.g. in46 and computed for each pair of single or aggregated diseases with an incidence of at least 4. A value of 1 means a "perfect match", 0 means no co-occurrence and negative values up to − 1 means "inverse occurrence". A comorbidity network approach was considered to investigate the relationships between diseases. We computed the dissimilarity between the conditions studied in the previous step as one minus the Pearson's phi. Hence, a value of 1 corresponds to "no association" and values larger than 1 to "inverse association". The hierarchical classification was performed through an agglomerative clustering algorithm with the Ward method since this proved to have the largest agglomeration coefficient. Furthermore, the dendrogram was cut at nodes corresponding to branch length of 1 based both on a visual inspection and on predefined length of the branch. This study has several limitations. We considered only the deaths up to April 20 when the peak of the first pandemic wave was reached in Romania. Moreover, total case design was used because of the missing data in official public communications. Disclaimer This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors. Ethical compliance This research was carried out in accordance with the relevant guidelines and regulations. Informed consent We have obtained informed consent from all participants and/or their legal guardians. Supplementary Information Supplementary Information.
sec	Introduction The rapid spread of the new type of coronavirus named novel Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-COV-2) from Wuhan (China) in the world determined the World Health Organization (WHO) to declare novel coronavirus disease (2019 COVID-19 or nCOVID-19) a pandemic on March 11, 20201,2. The first confirmed case in Romania was on February 26, 2020, and, by April 20 2020, 8936 confirmed case, 2017 healed cases and 520 deaths had been recorded, 88 cases with missing data and 98,491 tests were done3. The case definitions, as well as the therapies applied by the Romanian physicians, were also used in Italy, Spain and France4–7. The diagnostic was established based on the clinical symptomatology like cough, temperature, breathing difficulty without a prescribed etiology, on patient's direct exposure to SARS-COV-2, and the corroboration by testing with the RT-PCR method8,9. The Romanian government took strict measures to limit the outbreak, and a Coronavirus task force was set up to coordinate this effort under the direct supervision of the Ministry of Health (MS). Even though some cases of intra-hospital or community infection arose, they were isolated and contained with the help of special hygiene measures, social distancing and a very tight lockdown imposed by the authorities. However, a large number of deaths were confirmed among the adult population with no patient under 20 years of age dying because of COVID-19. Hence, issues regarding the COVID-19 mortality, especially the relationship with the comorbidities and the burden of disease are of great importance for the medical community8. Several studies about the comorbidities of the patients infected with SARS-COV-2 indicate hypertension, diabetes, obesity, neoplasms, chronic kidney disease and chronic obstructive pulmonary disease (COPD) as leading risk factors for a lethal evolution of the disease10–12. Romania holds a top place in Europe for deaths caused by cardiovascular diseases with hypertension affecting over 45% of the adult population. Almost half of the population aged over 65 has at least one chronic disease, and 1 in 5 adults smokes daily (32% men as compared to 18% women)13. Moreover, Romania has a record number of deaths due to preventable causes like ischemic cardiac disease, pulmonary cancer or alcoholism14. The prevalence of diabetes is approximately 11.6% and doubles for prediabetes as confirmed by Predatorr study15 and Mentor study16. Therefore, it is expected that a high number of infections with SARS-CoV-2 in Romania to have a severe course of the disease or to end deadly. To investigate the comorbidity profiles for SARS-CoV2 deaths, an observational retrospective study was conducted starting from the official public communications of Romanian Ministry of Internal Affairs (MAI) regarding the fatalities declared as COVID-19 deaths by the National Institute of Public Health (INSP)4,8. Additionally, data from the National Center for Statistics and Informatics in Public Health were utilized to run a case-case study between the COVID-19 death population and the fatalities due to hospital pneumonia between March 22 and April 20 in the years 2016–2018 in Romania. As far as we know, this is the first study in Eastern Europe that investigates patterns of comorbidities for SARS-COV-2 fatalities. Our principal research question is to generate hypothesis about comorbidity and burden of disease in the population deceased because of the COVID-19 virus and to investigate their association with gender and age. Additionally, the comparison between the fatalities due to COVID-19 and to pneumonia in terms of gender, age and comorbidity profiles allows us to characterize the specific medical pre-conditions of COVID-19 population.
title	Introduction
p	The rapid spread of the new type of coronavirus named novel Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-COV-2) from Wuhan (China) in the world determined the World Health Organization (WHO) to declare novel coronavirus disease (2019 COVID-19 or nCOVID-19) a pandemic on March 11, 20201,2.
p	The first confirmed case in Romania was on February 26, 2020, and, by April 20 2020, 8936 confirmed case, 2017 healed cases and 520 deaths had been recorded, 88 cases with missing data and 98,491 tests were done3. The case definitions, as well as the therapies applied by the Romanian physicians, were also used in Italy, Spain and France4–7. The diagnostic was established based on the clinical symptomatology like cough, temperature, breathing difficulty without a prescribed etiology, on patient's direct exposure to SARS-COV-2, and the corroboration by testing with the RT-PCR method8,9.
p	The Romanian government took strict measures to limit the outbreak, and a Coronavirus task force was set up to coordinate this effort under the direct supervision of the Ministry of Health (MS). Even though some cases of intra-hospital or community infection arose, they were isolated and contained with the help of special hygiene measures, social distancing and a very tight lockdown imposed by the authorities. However, a large number of deaths were confirmed among the adult population with no patient under 20 years of age dying because of COVID-19. Hence, issues regarding the COVID-19 mortality, especially the relationship with the comorbidities and the burden of disease are of great importance for the medical community8.
p	Several studies about the comorbidities of the patients infected with SARS-COV-2 indicate hypertension, diabetes, obesity, neoplasms, chronic kidney disease and chronic obstructive pulmonary disease (COPD) as leading risk factors for a lethal evolution of the disease10–12. Romania holds a top place in Europe for deaths caused by cardiovascular diseases with hypertension affecting over 45% of the adult population. Almost half of the population aged over 65 has at least one chronic disease, and 1 in 5 adults smokes daily (32% men as compared to 18% women)13. Moreover, Romania has a record number of deaths due to preventable causes like ischemic cardiac disease, pulmonary cancer or alcoholism14. The prevalence of diabetes is approximately 11.6% and doubles for prediabetes as confirmed by Predatorr study15 and Mentor study16. Therefore, it is expected that a high number of infections with SARS-CoV-2 in Romania to have a severe course of the disease or to end deadly.
p	To investigate the comorbidity profiles for SARS-CoV2 deaths, an observational retrospective study was conducted starting from the official public communications of Romanian Ministry of Internal Affairs (MAI) regarding the fatalities declared as COVID-19 deaths by the National Institute of Public Health (INSP)4,8. Additionally, data from the National Center for Statistics and Informatics in Public Health were utilized to run a case-case study between the COVID-19 death population and the fatalities due to hospital pneumonia between March 22 and April 20 in the years 2016–2018 in Romania. As far as we know, this is the first study in Eastern Europe that investigates patterns of comorbidities for SARS-COV-2 fatalities.
p	Our principal research question is to generate hypothesis about comorbidity and burden of disease in the population deceased because of the COVID-19 virus and to investigate their association with gender and age. Additionally, the comparison between the fatalities due to COVID-19 and to pneumonia in terms of gender, age and comorbidity profiles allows us to characterize the specific medical pre-conditions of COVID-19 population.
sec	Results Population characteristics Figure 1 illustrates the history of the SARS-CoV-2 pandemic in Romania from February 26 until April 20, including dates, number of fatalities and number of the missing subjects due to incomplete information. The Romanian COVID-19 mortality study has a sample size of 432 patients with complete data regarding gender, age and comorbidities. It consists of a multi-case, multi-center approach focused on the mortality frequency and of a case-case evaluation with the aim of comparing hospital-pneumonia and COVID-19 mortality (see the section on the study design for the details). Our goal is to investigate the association between the mortality and comorbidities, gender, age and hospital-pneumonia based on the statistical methods presented in detail in the last section of the present article. Since the Romanian population mainly consists of 84.1% Ethnic Romanians and of a small Hungarian minority that makes up 5.4% of the population demographics that both belong to the white ethnicity, an effect of the ethnicity was not considered in this article17. Figure 1 History of SARS-COV-2 pandemic in Romania up to April 20. First, the basic descriptive statistics about our sample are presented. There are 282 (65.3%) men and 150 (34.7%) women with a statistically significant difference between the size of these two groups (Table S1). The general mean age is 67 years, with an SD of 13.1, and the median value is 68 (Table S2). Sources of infection with deaths due to COVID-19 virus were reported in 36 counties in Romania with significant case concentration in Suceava-county n. = 67 (15.51%), Arges-county n. = 44 (10.19%) and Hunedoara-county n. = 41 (9.49%). Details regarding the spatial distribution of the distinct diseases and aggregated groups of diseases are displayed in Table S6 and Table S7 in the supplementary material. A total of 170 different conditions and 20 groups of diseases were identified in the study, with a broader diversity of 149 conditions for men over 107 conditions for women. The diseases with the highest prevalence in the sample are hypertension, obesity, diabetes and chronic kidney disease as well as diseases of the circulatory system and nutritional or metabolic disorders. The sample size of the pneumonia study is equal to 874 persons with complete data as to gender, age and comorbidities. The sample includes 492 (56.3%) men and 382 (43.7%) women, the average age is 73.3 years, with an SD of 13.5, and the median age is 76. Association between gender or age and comorbidities One-sample proportion test indicates a statistically significant difference in the relative frequency of man and women in the sample and in the subpopulations with hypertension, diabetes, diabetes mellitus type 2, chronic kidney diseases, dependence on renal dialysis, chronic obstructive pulmonary disease, acute ischaemic heart disease or respiratory failure as displayed in Table S1. This is also valid for the subpopulations without comorbidities and for those corresponding to the aggregated group of diseases with a relative frequency of more than 5% except for diseases of the nervous system. On the other hand, the variable gender is not significantly associated with single conditions or group of illnesses except for heart failure (see Table 1). Table 1 Association between gender and distinct diseases or aggregated group of diseases, as absolute numbers and as percentage. Characteristicsn,% Menn = 282 Womenn = 150 p-value Distinct diseases Hypertension 113 40.07% 49 32.67% 0.16 Obesity 32 11.35% 21 14.00% 0.52 Diabetes Total 105 37.23% 48 32% 0.33 Diabetes mellitus, unspecified 40 14.18% 26 17% 0.47 Diabetes mellitus type 2 56 19.86% 19 13% 0.081 Diabetes mellitus type 1 9 3.19% 3 2% 0.68 Diseases of the circulatory system, unspecified 41 14.54% 27 18.00% 0.42 Supraventricular tachyarrhythmia 16 5.67% 9 6.00% 1 Heart failure, unspecified 14 4.96% 16 10.67% 0.043 Cerebral ischaemic stroke 16 5.67% 12 8% 0.47 Chronic kidney disease 31 10.99% 13 8.67% 0.55 Kidney failure, unspecified 17 6.03% 9 6.00% 1 Dependence on renal dialysis 29 10.28% 14 9.33% 0.88 Chronic obstructive pulmonary disease 18 6.38% 7 4.67% 0.61 Acute ischaemic heart disease, unspecified 15 5.32% 4 2.67% 0.3 Respiratory failure 15 6.32% 4 2.67% 0.3 Aggregated group of diseases No comorbidities 11 3.90% 2 1.33% 0.23 Diseases of the circulatory system 168 59.57% 88 58.67% 0.94 Endocrine, nutritional or metabolic diseases 122 43.26% 63 42.00% 0.88 Diseases of the genitourinary system 54 19.15% 25 16.67% 0.61 Diseases of the respiratory system 47 16.67% 20 13.33% 0.44 Diseases of the nervous system 35 12.41% 26 17.33% 0.21 Diseases of the digestive system 30 10.64% 14 9.33% 0.79 Neoplasms 26 9.22% 14 9.33% 1 Mental, behavioural or neurodevelopmental disorders 25 8.87% 11 7.33% 0.71 Factors influencing health status or contact with health services 37 13.12% 17 11.33% 0.7 Symptoms, signs or clinical findings, not elsewhere classified 13 4.61% 11 7.33% 0.34 Certain infectious or parasitic diseases 15 5.32% 5 3.33% 0.49 In bold the differences that reached the statistical significance. Because age is an established risk factor for mortality, we also performed an age-stratified analysis with the lifespan split in 10-year intervals starting with age 50. The chi-squared test of association between the gender and age factors is statistically significant, with the largest male relative frequency (81.16%) between 50 and 59 years. Moreover, the factor age is also significantly associated with most of the comorbidities, with hypertension and dependence on renal dialysis following also this trend but at a significance level of 10% (see Table 2). Table 2 Association between age and distinct diseases or group of diseases, as absolute numbers and as percentage. Characteristicsn, % < 504410.19% 50–596915.97% 60–6912027.78% 70–7912829.93% ≥ 807116.44% p-value Gender Male 29 66% 56 81.16% 73 60.83% 80 62.5% 44 61.97% 0.05 Female 15 34.09% 13 19.84% 47 39.17% 48 37.5% 27 38.03% Distinct diseases Hypertension 12 27.28% 21 30.43% 49 40.83% 45 35.16% 35 49.30% 0.074 Obesity 13 29.55% 9 13.04% 20 16.67% 8 6.25% 3 4.23% 0.0001 Diabetes Total 14 31.82% 23 33.33% 53 44.17% 48 37.5% 16 22.53% 0.045 Diabetes mellitus, unspecified 7 15.91% 8 11.59% 18 15.00% 24 18.75% 10 14.08% 0.74 Diabetes mellitus type 2 7 15.91% 12 17.39% 31 25.83% 20 15.63% 5 7.04% 0.021 Diabetes mellitus type 1 – 3 4.34% 4 3.33% 4 3.13% 1 1.40% – Diseases of the circulatory system, unspecified 6 13.64% 7 10.14% 18 15.00% 23 17.97% 14 19.72% 0.53 Supraventricular tachyarrhythmia 1 2.27% – – 12 9.38% 10 14.08% – Heart failure, unspecified – 5 7.25% 10 8.33% ( 7 5.47% 8 11.27% – Cerebral ischaemic stroke 2 4.55% 1 1.49% 9 7.50% 9 7.03% 7 9.86% 0.31 Acute ischaemic heart disease, unspecified – 2 2.90% 7 5.83% 4 3.13% 6 8.45% – Chronic kidney disease 4 9.09% 3 4.35% 14 11.67% 14 10.94% 9 12.67% 0.48 Dependence on renal dialysis 3 6.82% 13 13.04% 13 10.83% 15 11.72% 3 4.23% 0.072 Kidney failure, unspecified 3 6.82% 3 4.35% 9 7.50% 8 6.25% 3 4.23% 0.87 Respiratory failure 2 4.55% 4 5.80% ) 3 2.50% 4 3.13% 6 8.45% 0.33 Chronic obstructive pulmonary disease 1 2.27% – 12 10.00% 10 7.81% 2 2.82% – Aggregated group of diseases Diseases of the circulatory system 14 31.82% 33 47.83% 74 61.67% 82 64.06% 53 74.65% 0.0001 Endocrine, nutritional or metabolic diseases 21 47.73% 30 43.48% 63 52.50% 53 41.41% 18 25.35% 0.007 Diseases of the genitourinary system 7 15.91% 7 10.14% 24 20.00% 29 22.66% 12 16.90% 0.27 Diseases of the respiratory system 6 13.64% 10 14.49% 18 15.00% 20 15.63% 13 18.31% 0.96 Diseases of the nervous system 3 6.82% 2 2.90% 18 15.00% 19 14.84% 19 26.76% 0.001 Diseases of the digestive system 4 9.09% 8 11.59% 18 15.00% 10 7.81% 4 5.63% 0.23 Neoplasms 1 2.27% 7 10.14% 12 10.00% 16 12.50% 4 5.63% 0.25 Mental, behavioural or neurodevelopmental disorders 3 6.82% 2 2.90% 6 5.00% 15 11.72% 10 14.08% 0.049 Factors influencing health status or contact with health services 4 9.09% 10 14.49% 15 12.50% 19 14.84% 6 8.45% 0.66 Symptoms, signs or clinical findings, not elsewhere classified 2 4.55% 3 4.35% 6 5.00% 11 8.59% 2 2.82% 0.46 Certain infectious or parasitic diseases 3 6.82% 2 2.90% 6 5.00% 4 3.13% 5 7.04% 0.63 In bold the p-values that reached the statistical significance. Comparison between COVID-19 and pneumonia fatalities The comparison with pneumonia fatalities yielded that the proportion of males in the COVID-19 death sample is also significantly higher than in the pneumonia death sample as shown in Table 3. Age and the relative frequency of essential hypertension, diseases of the circulatory system, obesity, dependence of renal dialysis, and chronic kidney failure are also different between the COVID-19 and pneumonia groups. Additionally, the odd to die for COVID-19 versus hospital pneumonia is 150% greater for men than women with similar comorbidities (see Table 4). Moreover, the adjusted odds to die for COVID-19 compared to pneumonia are 22 times higher for persons with the chronic kidney failure, doubled for persons with hypertension, and almost equal for patients with type 2 diabetes. Table 3 Comparison between the characteristics of the COVID-19 and pneumonia groups. Characteristics COVID-19 Pneumonia p-value Sample size n, % 432 100% 711 100% Gender n, % 0.0001 Female 150 35% 320 45% Male 282 65% 391 54.99% Age Mean (sd) Median 66.97 (13.07) 68 73.14 (13.68) 78 Age group n, % 0.0001 < 50 44 10.19% 49 6.89% 50–59 69 15.97% 51 7.17% 60–69 120 27.78% 141 19.83% 70–79 128 29.63% 199 27.99% ≥ 80 71 16.44% 271 38.12% Comorbidities Hypertension 162 37.50% 181 25.46% 0.0001 Obesity 56 12.96% 63 8.86% 0.03 Diabetes mellitus type I 12 2.78% 5 0.70% 0.01 Diabetes mellitus type II 75 17.36% 145 20.39% 0.21 Dialysis 43 9.95% 4 0.56% 0.0001 Chronic kidney disease 43 9.95% 88 12.38% 0.21 Supraventricular tachyarrhythmia 25 5.79% 201 28.27% 0.0001 Congestive heart failure 12 2.78% 204 28.69% 0.0001 COPD 24 5.56% 72 10.13% 0.01 Heart failure 29 6.71% 27 3.80% 0.03 Kidney failure 26 6.02% 4 0.56% 0.0001 Respiratory failure 19 4.40% 88 12.38% 0.0001 Cerebral ischaemic stroke 28 6.48% – 0.0001 In bold the p-values that reached the statistical significance. Table 4 Logistic regression analysis for the comparison between COVID-19 and pneumonia death groups. Covariates CoefficientΒ OR (COVID-19 vs. pneumonia) 95% CI p-value Age − 0.022 0.978 [0.968–0.988] 0.0001 Gender (0 = female) 0.414 1.513 [1.145–2.000] 0.01 Hypertension 0.740 2.096 [1.560–2.818] 0.0001 Obesity 0.266 1.305 [0.843–2.019] 0.23 Diabetes type 2 − 0.355 0.701 [0.494–0.996] 0.049 Dialysis 3.104 22.290 [7.244–68.590] 0.0001 Chronic kidney disease − 0.522 0.593 [0.369–0.955] 0.03 Heart failure 1.095 2.990 [1.612–5.546] 0.0001 COPD − 0.730 0.482 [0.284–0.817] 0.01 Supraventricular tachyarrhythmia − 1.687 0.185 [0.116–0.295] 0.0001 Respiratory failure − 1.045 0.352 [0.203–0.609] 0.0001 Intercept 0.884 2.420 0.02 In bold the p-values that reached the statistical significance. The burden of disease and survival probability We also extended the analysis to the coexistence of multiple health conditions, and we focused on the multimorbidity. This has a mean value of 2.109, an SD of 1.66 and a median value of 2 in our study. The multimorbidity factor with levels 0, 1, 2, 3, 4 or more than 5 for single or aggregated comorbidities is significantly associated with the factor age but not with gender, confirming the pattern displayed at the level of single conditions and group of diseases (Tables S3 and S4). The Charlson comorbidity index (CCI) had a mean of 1.324, SD of 0.95 and median of one indicating that the severity of the underlying medical condition in the study was mainly mild. This was confirmed by the high relative frequency of 91.20% for the mild level for the factor CCI. There were no significant differences between males and females or among age groups with respect to CCI (Tables S3 and S4). Moreover, 1-year survival probability averaged out at 81%, and 50% of the study individuals had a prognosis for it of 85%. Co-occurrence patterns The co-occurrence of the diseases is of great importance since the mean multimorbidity of 2.109 indicates that not single but pair of diseases are pre-conditions in our mortality study. Pearson's phi coefficient is a robust measure of the comorbidity association, and its values for diseases in general and separated for men and women are displayed in the Table S5. Different co-occurrence patterns emerge between men and women with the striking lack of the correlation between pneumonia or respiratory failure, respectively stroke and dissociative neurological symptom disorder in the female profile. The largest value of the parameter phi of 0.47 respectively 0.40 for female and 0.59 for male corresponds to kidney failure and dependence to dialysis, while the latter is also correlated with chronic kidney disease as the correlation lies by 0.30 in general and 0.38 in men. They are followed by pneumonia or respiratory failure with phi equal to 0.45 in general and to 0.55 in male and by stroke and dissociative neurological symptom disorder with a correlation coefficient of 0.41 in the whole study and of 0.49 in men. Diabetes mellitus type 2 is correlated with hypertension not only for the whole study (phi = 0.11) but also in the men (phi = 0.10) and women subsamples (phi = 0.11). Hypertension appears mainly together with stroke with phi equal to 0.15, 0.16 respectively 0.14 for the whole study, men and women. Finally, a comorbidity network analysis rounded off our exploratory investigation. This identified 11 clusters in the general study, 9 clusters in male, and 6 clusters in female subgroups, each including two or three diseases reinforcing the previous results about pair of diseases as pre-condition for COVID-19 fatalities (see Fig. 2). Figure 2 Co-occurrence clusters for All Subjects (A), Women (B) and Men (C).
title	Results
sec	Population characteristics Figure 1 illustrates the history of the SARS-CoV-2 pandemic in Romania from February 26 until April 20, including dates, number of fatalities and number of the missing subjects due to incomplete information. The Romanian COVID-19 mortality study has a sample size of 432 patients with complete data regarding gender, age and comorbidities. It consists of a multi-case, multi-center approach focused on the mortality frequency and of a case-case evaluation with the aim of comparing hospital-pneumonia and COVID-19 mortality (see the section on the study design for the details). Our goal is to investigate the association between the mortality and comorbidities, gender, age and hospital-pneumonia based on the statistical methods presented in detail in the last section of the present article. Since the Romanian population mainly consists of 84.1% Ethnic Romanians and of a small Hungarian minority that makes up 5.4% of the population demographics that both belong to the white ethnicity, an effect of the ethnicity was not considered in this article17. Figure 1 History of SARS-COV-2 pandemic in Romania up to April 20. First, the basic descriptive statistics about our sample are presented. There are 282 (65.3%) men and 150 (34.7%) women with a statistically significant difference between the size of these two groups (Table S1). The general mean age is 67 years, with an SD of 13.1, and the median value is 68 (Table S2). Sources of infection with deaths due to COVID-19 virus were reported in 36 counties in Romania with significant case concentration in Suceava-county n. = 67 (15.51%), Arges-county n. = 44 (10.19%) and Hunedoara-county n. = 41 (9.49%). Details regarding the spatial distribution of the distinct diseases and aggregated groups of diseases are displayed in Table S6 and Table S7 in the supplementary material. A total of 170 different conditions and 20 groups of diseases were identified in the study, with a broader diversity of 149 conditions for men over 107 conditions for women. The diseases with the highest prevalence in the sample are hypertension, obesity, diabetes and chronic kidney disease as well as diseases of the circulatory system and nutritional or metabolic disorders. The sample size of the pneumonia study is equal to 874 persons with complete data as to gender, age and comorbidities. The sample includes 492 (56.3%) men and 382 (43.7%) women, the average age is 73.3 years, with an SD of 13.5, and the median age is 76.
title	Population characteristics
p	Figure 1 illustrates the history of the SARS-CoV-2 pandemic in Romania from February 26 until April 20, including dates, number of fatalities and number of the missing subjects due to incomplete information. The Romanian COVID-19 mortality study has a sample size of 432 patients with complete data regarding gender, age and comorbidities. It consists of a multi-case, multi-center approach focused on the mortality frequency and of a case-case evaluation with the aim of comparing hospital-pneumonia and COVID-19 mortality (see the section on the study design for the details). Our goal is to investigate the association between the mortality and comorbidities, gender, age and hospital-pneumonia based on the statistical methods presented in detail in the last section of the present article. Since the Romanian population mainly consists of 84.1% Ethnic Romanians and of a small Hungarian minority that makes up 5.4% of the population demographics that both belong to the white ethnicity, an effect of the ethnicity was not considered in this article17. Figure 1 History of SARS-COV-2 pandemic in Romania up to April 20.
figure	Figure 1 History of SARS-COV-2 pandemic in Romania up to April 20.
label	Figure 1
caption	History of SARS-COV-2 pandemic in Romania up to April 20.
p	History of SARS-COV-2 pandemic in Romania up to April 20.
p	First, the basic descriptive statistics about our sample are presented. There are 282 (65.3%) men and 150 (34.7%) women with a statistically significant difference between the size of these two groups (Table S1). The general mean age is 67 years, with an SD of 13.1, and the median value is 68 (Table S2).
p	Sources of infection with deaths due to COVID-19 virus were reported in 36 counties in Romania with significant case concentration in Suceava-county n. = 67 (15.51%), Arges-county n. = 44 (10.19%) and Hunedoara-county n. = 41 (9.49%). Details regarding the spatial distribution of the distinct diseases and aggregated groups of diseases are displayed in Table S6 and Table S7 in the supplementary material.
p	A total of 170 different conditions and 20 groups of diseases were identified in the study, with a broader diversity of 149 conditions for men over 107 conditions for women. The diseases with the highest prevalence in the sample are hypertension, obesity, diabetes and chronic kidney disease as well as diseases of the circulatory system and nutritional or metabolic disorders.
p	The sample size of the pneumonia study is equal to 874 persons with complete data as to gender, age and comorbidities. The sample includes 492 (56.3%) men and 382 (43.7%) women, the average age is 73.3 years, with an SD of 13.5, and the median age is 76.
sec	Association between gender or age and comorbidities One-sample proportion test indicates a statistically significant difference in the relative frequency of man and women in the sample and in the subpopulations with hypertension, diabetes, diabetes mellitus type 2, chronic kidney diseases, dependence on renal dialysis, chronic obstructive pulmonary disease, acute ischaemic heart disease or respiratory failure as displayed in Table S1. This is also valid for the subpopulations without comorbidities and for those corresponding to the aggregated group of diseases with a relative frequency of more than 5% except for diseases of the nervous system. On the other hand, the variable gender is not significantly associated with single conditions or group of illnesses except for heart failure (see Table 1). Table 1 Association between gender and distinct diseases or aggregated group of diseases, as absolute numbers and as percentage. Characteristicsn,% Menn = 282 Womenn = 150 p-value Distinct diseases Hypertension 113 40.07% 49 32.67% 0.16 Obesity 32 11.35% 21 14.00% 0.52 Diabetes Total 105 37.23% 48 32% 0.33 Diabetes mellitus, unspecified 40 14.18% 26 17% 0.47 Diabetes mellitus type 2 56 19.86% 19 13% 0.081 Diabetes mellitus type 1 9 3.19% 3 2% 0.68 Diseases of the circulatory system, unspecified 41 14.54% 27 18.00% 0.42 Supraventricular tachyarrhythmia 16 5.67% 9 6.00% 1 Heart failure, unspecified 14 4.96% 16 10.67% 0.043 Cerebral ischaemic stroke 16 5.67% 12 8% 0.47 Chronic kidney disease 31 10.99% 13 8.67% 0.55 Kidney failure, unspecified 17 6.03% 9 6.00% 1 Dependence on renal dialysis 29 10.28% 14 9.33% 0.88 Chronic obstructive pulmonary disease 18 6.38% 7 4.67% 0.61 Acute ischaemic heart disease, unspecified 15 5.32% 4 2.67% 0.3 Respiratory failure 15 6.32% 4 2.67% 0.3 Aggregated group of diseases No comorbidities 11 3.90% 2 1.33% 0.23 Diseases of the circulatory system 168 59.57% 88 58.67% 0.94 Endocrine, nutritional or metabolic diseases 122 43.26% 63 42.00% 0.88 Diseases of the genitourinary system 54 19.15% 25 16.67% 0.61 Diseases of the respiratory system 47 16.67% 20 13.33% 0.44 Diseases of the nervous system 35 12.41% 26 17.33% 0.21 Diseases of the digestive system 30 10.64% 14 9.33% 0.79 Neoplasms 26 9.22% 14 9.33% 1 Mental, behavioural or neurodevelopmental disorders 25 8.87% 11 7.33% 0.71 Factors influencing health status or contact with health services 37 13.12% 17 11.33% 0.7 Symptoms, signs or clinical findings, not elsewhere classified 13 4.61% 11 7.33% 0.34 Certain infectious or parasitic diseases 15 5.32% 5 3.33% 0.49 In bold the differences that reached the statistical significance. Because age is an established risk factor for mortality, we also performed an age-stratified analysis with the lifespan split in 10-year intervals starting with age 50. The chi-squared test of association between the gender and age factors is statistically significant, with the largest male relative frequency (81.16%) between 50 and 59 years. Moreover, the factor age is also significantly associated with most of the comorbidities, with hypertension and dependence on renal dialysis following also this trend but at a significance level of 10% (see Table 2). Table 2 Association between age and distinct diseases or group of diseases, as absolute numbers and as percentage. Characteristicsn, % < 504410.19% 50–596915.97% 60–6912027.78% 70–7912829.93% ≥ 807116.44% p-value Gender Male 29 66% 56 81.16% 73 60.83% 80 62.5% 44 61.97% 0.05 Female 15 34.09% 13 19.84% 47 39.17% 48 37.5% 27 38.03% Distinct diseases Hypertension 12 27.28% 21 30.43% 49 40.83% 45 35.16% 35 49.30% 0.074 Obesity 13 29.55% 9 13.04% 20 16.67% 8 6.25% 3 4.23% 0.0001 Diabetes Total 14 31.82% 23 33.33% 53 44.17% 48 37.5% 16 22.53% 0.045 Diabetes mellitus, unspecified 7 15.91% 8 11.59% 18 15.00% 24 18.75% 10 14.08% 0.74 Diabetes mellitus type 2 7 15.91% 12 17.39% 31 25.83% 20 15.63% 5 7.04% 0.021 Diabetes mellitus type 1 – 3 4.34% 4 3.33% 4 3.13% 1 1.40% – Diseases of the circulatory system, unspecified 6 13.64% 7 10.14% 18 15.00% 23 17.97% 14 19.72% 0.53 Supraventricular tachyarrhythmia 1 2.27% – – 12 9.38% 10 14.08% – Heart failure, unspecified – 5 7.25% 10 8.33% ( 7 5.47% 8 11.27% – Cerebral ischaemic stroke 2 4.55% 1 1.49% 9 7.50% 9 7.03% 7 9.86% 0.31 Acute ischaemic heart disease, unspecified – 2 2.90% 7 5.83% 4 3.13% 6 8.45% – Chronic kidney disease 4 9.09% 3 4.35% 14 11.67% 14 10.94% 9 12.67% 0.48 Dependence on renal dialysis 3 6.82% 13 13.04% 13 10.83% 15 11.72% 3 4.23% 0.072 Kidney failure, unspecified 3 6.82% 3 4.35% 9 7.50% 8 6.25% 3 4.23% 0.87 Respiratory failure 2 4.55% 4 5.80% ) 3 2.50% 4 3.13% 6 8.45% 0.33 Chronic obstructive pulmonary disease 1 2.27% – 12 10.00% 10 7.81% 2 2.82% – Aggregated group of diseases Diseases of the circulatory system 14 31.82% 33 47.83% 74 61.67% 82 64.06% 53 74.65% 0.0001 Endocrine, nutritional or metabolic diseases 21 47.73% 30 43.48% 63 52.50% 53 41.41% 18 25.35% 0.007 Diseases of the genitourinary system 7 15.91% 7 10.14% 24 20.00% 29 22.66% 12 16.90% 0.27 Diseases of the respiratory system 6 13.64% 10 14.49% 18 15.00% 20 15.63% 13 18.31% 0.96 Diseases of the nervous system 3 6.82% 2 2.90% 18 15.00% 19 14.84% 19 26.76% 0.001 Diseases of the digestive system 4 9.09% 8 11.59% 18 15.00% 10 7.81% 4 5.63% 0.23 Neoplasms 1 2.27% 7 10.14% 12 10.00% 16 12.50% 4 5.63% 0.25 Mental, behavioural or neurodevelopmental disorders 3 6.82% 2 2.90% 6 5.00% 15 11.72% 10 14.08% 0.049 Factors influencing health status or contact with health services 4 9.09% 10 14.49% 15 12.50% 19 14.84% 6 8.45% 0.66 Symptoms, signs or clinical findings, not elsewhere classified 2 4.55% 3 4.35% 6 5.00% 11 8.59% 2 2.82% 0.46 Certain infectious or parasitic diseases 3 6.82% 2 2.90% 6 5.00% 4 3.13% 5 7.04% 0.63 In bold the p-values that reached the statistical significance.
title	Association between gender or age and comorbidities
p	One-sample proportion test indicates a statistically significant difference in the relative frequency of man and women in the sample and in the subpopulations with hypertension, diabetes, diabetes mellitus type 2, chronic kidney diseases, dependence on renal dialysis, chronic obstructive pulmonary disease, acute ischaemic heart disease or respiratory failure as displayed in Table S1. This is also valid for the subpopulations without comorbidities and for those corresponding to the aggregated group of diseases with a relative frequency of more than 5% except for diseases of the nervous system. On the other hand, the variable gender is not significantly associated with single conditions or group of illnesses except for heart failure (see Table 1). Table 1 Association between gender and distinct diseases or aggregated group of diseases, as absolute numbers and as percentage. Characteristicsn,% Menn = 282 Womenn = 150 p-value Distinct diseases Hypertension 113 40.07% 49 32.67% 0.16 Obesity 32 11.35% 21 14.00% 0.52 Diabetes Total 105 37.23% 48 32% 0.33 Diabetes mellitus, unspecified 40 14.18% 26 17% 0.47 Diabetes mellitus type 2 56 19.86% 19 13% 0.081 Diabetes mellitus type 1 9 3.19% 3 2% 0.68 Diseases of the circulatory system, unspecified 41 14.54% 27 18.00% 0.42 Supraventricular tachyarrhythmia 16 5.67% 9 6.00% 1 Heart failure, unspecified 14 4.96% 16 10.67% 0.043 Cerebral ischaemic stroke 16 5.67% 12 8% 0.47 Chronic kidney disease 31 10.99% 13 8.67% 0.55 Kidney failure, unspecified 17 6.03% 9 6.00% 1 Dependence on renal dialysis 29 10.28% 14 9.33% 0.88 Chronic obstructive pulmonary disease 18 6.38% 7 4.67% 0.61 Acute ischaemic heart disease, unspecified 15 5.32% 4 2.67% 0.3 Respiratory failure 15 6.32% 4 2.67% 0.3 Aggregated group of diseases No comorbidities 11 3.90% 2 1.33% 0.23 Diseases of the circulatory system 168 59.57% 88 58.67% 0.94 Endocrine, nutritional or metabolic diseases 122 43.26% 63 42.00% 0.88 Diseases of the genitourinary system 54 19.15% 25 16.67% 0.61 Diseases of the respiratory system 47 16.67% 20 13.33% 0.44 Diseases of the nervous system 35 12.41% 26 17.33% 0.21 Diseases of the digestive system 30 10.64% 14 9.33% 0.79 Neoplasms 26 9.22% 14 9.33% 1 Mental, behavioural or neurodevelopmental disorders 25 8.87% 11 7.33% 0.71 Factors influencing health status or contact with health services 37 13.12% 17 11.33% 0.7 Symptoms, signs or clinical findings, not elsewhere classified 13 4.61% 11 7.33% 0.34 Certain infectious or parasitic diseases 15 5.32% 5 3.33% 0.49 In bold the differences that reached the statistical significance.
table-wrap	Table 1 Association between gender and distinct diseases or aggregated group of diseases, as absolute numbers and as percentage. Characteristicsn,% Menn = 282 Womenn = 150 p-value Distinct diseases Hypertension 113 40.07% 49 32.67% 0.16 Obesity 32 11.35% 21 14.00% 0.52 Diabetes Total 105 37.23% 48 32% 0.33 Diabetes mellitus, unspecified 40 14.18% 26 17% 0.47 Diabetes mellitus type 2 56 19.86% 19 13% 0.081 Diabetes mellitus type 1 9 3.19% 3 2% 0.68 Diseases of the circulatory system, unspecified 41 14.54% 27 18.00% 0.42 Supraventricular tachyarrhythmia 16 5.67% 9 6.00% 1 Heart failure, unspecified 14 4.96% 16 10.67% 0.043 Cerebral ischaemic stroke 16 5.67% 12 8% 0.47 Chronic kidney disease 31 10.99% 13 8.67% 0.55 Kidney failure, unspecified 17 6.03% 9 6.00% 1 Dependence on renal dialysis 29 10.28% 14 9.33% 0.88 Chronic obstructive pulmonary disease 18 6.38% 7 4.67% 0.61 Acute ischaemic heart disease, unspecified 15 5.32% 4 2.67% 0.3 Respiratory failure 15 6.32% 4 2.67% 0.3 Aggregated group of diseases No comorbidities 11 3.90% 2 1.33% 0.23 Diseases of the circulatory system 168 59.57% 88 58.67% 0.94 Endocrine, nutritional or metabolic diseases 122 43.26% 63 42.00% 0.88 Diseases of the genitourinary system 54 19.15% 25 16.67% 0.61 Diseases of the respiratory system 47 16.67% 20 13.33% 0.44 Diseases of the nervous system 35 12.41% 26 17.33% 0.21 Diseases of the digestive system 30 10.64% 14 9.33% 0.79 Neoplasms 26 9.22% 14 9.33% 1 Mental, behavioural or neurodevelopmental disorders 25 8.87% 11 7.33% 0.71 Factors influencing health status or contact with health services 37 13.12% 17 11.33% 0.7 Symptoms, signs or clinical findings, not elsewhere classified 13 4.61% 11 7.33% 0.34 Certain infectious or parasitic diseases 15 5.32% 5 3.33% 0.49 In bold the differences that reached the statistical significance.
label	Table 1
caption	Association between gender and distinct diseases or aggregated group of diseases, as absolute numbers and as percentage.
p	Association between gender and distinct diseases or aggregated group of diseases, as absolute numbers and as percentage.
table	Characteristicsn,% Menn = 282 Womenn = 150 p-value Distinct diseases Hypertension 113 40.07% 49 32.67% 0.16 Obesity 32 11.35% 21 14.00% 0.52 Diabetes Total 105 37.23% 48 32% 0.33 Diabetes mellitus, unspecified 40 14.18% 26 17% 0.47 Diabetes mellitus type 2 56 19.86% 19 13% 0.081 Diabetes mellitus type 1 9 3.19% 3 2% 0.68 Diseases of the circulatory system, unspecified 41 14.54% 27 18.00% 0.42 Supraventricular tachyarrhythmia 16 5.67% 9 6.00% 1 Heart failure, unspecified 14 4.96% 16 10.67% 0.043 Cerebral ischaemic stroke 16 5.67% 12 8% 0.47 Chronic kidney disease 31 10.99% 13 8.67% 0.55 Kidney failure, unspecified 17 6.03% 9 6.00% 1 Dependence on renal dialysis 29 10.28% 14 9.33% 0.88 Chronic obstructive pulmonary disease 18 6.38% 7 4.67% 0.61 Acute ischaemic heart disease, unspecified 15 5.32% 4 2.67% 0.3 Respiratory failure 15 6.32% 4 2.67% 0.3 Aggregated group of diseases No comorbidities 11 3.90% 2 1.33% 0.23 Diseases of the circulatory system 168 59.57% 88 58.67% 0.94 Endocrine, nutritional or metabolic diseases 122 43.26% 63 42.00% 0.88 Diseases of the genitourinary system 54 19.15% 25 16.67% 0.61 Diseases of the respiratory system 47 16.67% 20 13.33% 0.44 Diseases of the nervous system 35 12.41% 26 17.33% 0.21 Diseases of the digestive system 30 10.64% 14 9.33% 0.79 Neoplasms 26 9.22% 14 9.33% 1 Mental, behavioural or neurodevelopmental disorders 25 8.87% 11 7.33% 0.71 Factors influencing health status or contact with health services 37 13.12% 17 11.33% 0.7 Symptoms, signs or clinical findings, not elsewhere classified 13 4.61% 11 7.33% 0.34 Certain infectious or parasitic diseases 15 5.32% 5 3.33% 0.49
tr	Characteristicsn,% Menn = 282 Womenn = 150 p-value
th	Characteristicsn,%
th	Menn = 282
th	Womenn = 150
th	p-value
tr	Distinct diseases
td	Distinct diseases
tr	Hypertension 113 40.07% 49 32.67% 0.16
td	Hypertension
td	113 40.07%
p	113
p	40.07%
td	49 32.67%
p	49
p	32.67%
td	0.16
tr	Obesity 32 11.35% 21 14.00% 0.52
td	Obesity
td	32 11.35%
p	32
p	11.35%
td	21 14.00%
p	21
p	14.00%
td	0.52
tr	Diabetes
td	Diabetes
tr	Total 105 37.23% 48 32% 0.33
td	Total
td	105 37.23%
p	105
p	37.23%
td	48 32%
p	48
p	32%
td	0.33
tr	Diabetes mellitus, unspecified 40 14.18% 26 17% 0.47
td	Diabetes mellitus, unspecified
td	40 14.18%
p	40
p	14.18%
td	26 17%
p	26
p	17%
td	0.47
tr	Diabetes mellitus type 2 56 19.86% 19 13% 0.081
td	Diabetes mellitus type 2
td	56 19.86%
p	56
p	19.86%
td	19 13%
p	19
p	13%
td	0.081
tr	Diabetes mellitus type 1 9 3.19% 3 2% 0.68
td	Diabetes mellitus type 1
td	9 3.19%
p	9
p	3.19%
td	3 2%
p	3
p	2%
td	0.68
tr	Diseases of the circulatory system, unspecified 41 14.54% 27 18.00% 0.42
td	Diseases of the circulatory system, unspecified
td	41 14.54%
p	41
p	14.54%
td	27 18.00%
p	27
p	18.00%
td	0.42
tr	Supraventricular tachyarrhythmia 16 5.67% 9 6.00% 1
td	Supraventricular tachyarrhythmia
td	16 5.67%
p	16
p	5.67%
td	9 6.00%
p	9
p	6.00%
td	1
tr	Heart failure, unspecified 14 4.96% 16 10.67% 0.043
td	Heart failure, unspecified
td	14 4.96%
p	14
p	4.96%
td	16 10.67%
p	16
p	10.67%
td	0.043
tr	Cerebral ischaemic stroke 16 5.67% 12 8% 0.47
td	Cerebral ischaemic stroke
td	16 5.67%
p	16
p	5.67%
td	12 8%
p	12
p	8%
td	0.47
tr	Chronic kidney disease 31 10.99% 13 8.67% 0.55
td	Chronic kidney disease
td	31 10.99%
p	31
p	10.99%
td	13 8.67%
p	13
p	8.67%
td	0.55
tr	Kidney failure, unspecified 17 6.03% 9 6.00% 1
td	Kidney failure, unspecified
td	17 6.03%
p	17
p	6.03%
td	9 6.00%
p	9
p	6.00%
td	1
tr	Dependence on renal dialysis 29 10.28% 14 9.33% 0.88
td	Dependence on renal dialysis
td	29 10.28%
p	29
p	10.28%
td	14 9.33%
p	14
p	9.33%
td	0.88
tr	Chronic obstructive pulmonary disease 18 6.38% 7 4.67% 0.61
td	Chronic obstructive pulmonary disease
td	18 6.38%
p	18
p	6.38%
td	7 4.67%
p	7
p	4.67%
td	0.61
tr	Acute ischaemic heart disease, unspecified 15 5.32% 4 2.67% 0.3
td	Acute ischaemic heart disease, unspecified
td	15 5.32%
p	15
p	5.32%
td	4 2.67%
p	4
p	2.67%
td	0.3
tr	Respiratory failure 15 6.32% 4 2.67% 0.3
td	Respiratory failure
td	15 6.32%
p	15
p	6.32%
td	4 2.67%
p	4
p	2.67%
td	0.3
tr	Aggregated group of diseases
td	Aggregated group of diseases
tr	No comorbidities 11 3.90% 2 1.33% 0.23
td	No comorbidities
td	11 3.90%
p	11
p	3.90%
td	2 1.33%
p	2
p	1.33%
td	0.23
tr	Diseases of the circulatory system 168 59.57% 88 58.67% 0.94
td	Diseases of the circulatory system
td	168 59.57%
p	168
p	59.57%
td	88 58.67%
p	88
p	58.67%
td	0.94
tr	Endocrine, nutritional or metabolic diseases 122 43.26% 63 42.00% 0.88
td	Endocrine, nutritional or metabolic diseases
td	122 43.26%
p	122
p	43.26%
td	63 42.00%
p	63
p	42.00%
td	0.88
tr	Diseases of the genitourinary system 54 19.15% 25 16.67% 0.61
td	Diseases of the genitourinary system
td	54 19.15%
p	54
p	19.15%
td	25 16.67%
p	25
p	16.67%
td	0.61
tr	Diseases of the respiratory system 47 16.67% 20 13.33% 0.44
td	Diseases of the respiratory system
td	47 16.67%
p	47
p	16.67%
td	20 13.33%
p	20
p	13.33%
td	0.44
tr	Diseases of the nervous system 35 12.41% 26 17.33% 0.21
td	Diseases of the nervous system
td	35 12.41%
p	35
p	12.41%
td	26 17.33%
p	26
p	17.33%
td	0.21
tr	Diseases of the digestive system 30 10.64% 14 9.33% 0.79
td	Diseases of the digestive system
td	30 10.64%
p	30
p	10.64%
td	14 9.33%
p	14
p	9.33%
td	0.79
tr	Neoplasms 26 9.22% 14 9.33% 1
td	Neoplasms
td	26 9.22%
p	26
p	9.22%
td	14 9.33%
p	14
p	9.33%
td	1
tr	Mental, behavioural or neurodevelopmental disorders 25 8.87% 11 7.33% 0.71
td	Mental, behavioural or neurodevelopmental disorders
td	25 8.87%
p	25
p	8.87%
td	11 7.33%
p	11
p	7.33%
td	0.71
tr	Factors influencing health status or contact with health services 37 13.12% 17 11.33% 0.7
td	Factors influencing health status or contact with health services
td	37 13.12%
p	37
p	13.12%
td	17 11.33%
p	17
p	11.33%
td	0.7
tr	Symptoms, signs or clinical findings, not elsewhere classified 13 4.61% 11 7.33% 0.34
td	Symptoms, signs or clinical findings, not elsewhere classified
td	13 4.61%
p	13
p	4.61%
td	11 7.33%
p	11
p	7.33%
td	0.34
tr	Certain infectious or parasitic diseases 15 5.32% 5 3.33% 0.49
td	Certain infectious or parasitic diseases
td	15 5.32%
p	15
p	5.32%
td	5 3.33%
p	5
p	3.33%
td	0.49
table-wrap-foot	In bold the differences that reached the statistical significance.
p	In bold the differences that reached the statistical significance.
p	Because age is an established risk factor for mortality, we also performed an age-stratified analysis with the lifespan split in 10-year intervals starting with age 50. The chi-squared test of association between the gender and age factors is statistically significant, with the largest male relative frequency (81.16%) between 50 and 59 years. Moreover, the factor age is also significantly associated with most of the comorbidities, with hypertension and dependence on renal dialysis following also this trend but at a significance level of 10% (see Table 2). Table 2 Association between age and distinct diseases or group of diseases, as absolute numbers and as percentage. Characteristicsn, % < 504410.19% 50–596915.97% 60–6912027.78% 70–7912829.93% ≥ 807116.44% p-value Gender Male 29 66% 56 81.16% 73 60.83% 80 62.5% 44 61.97% 0.05 Female 15 34.09% 13 19.84% 47 39.17% 48 37.5% 27 38.03% Distinct diseases Hypertension 12 27.28% 21 30.43% 49 40.83% 45 35.16% 35 49.30% 0.074 Obesity 13 29.55% 9 13.04% 20 16.67% 8 6.25% 3 4.23% 0.0001 Diabetes Total 14 31.82% 23 33.33% 53 44.17% 48 37.5% 16 22.53% 0.045 Diabetes mellitus, unspecified 7 15.91% 8 11.59% 18 15.00% 24 18.75% 10 14.08% 0.74 Diabetes mellitus type 2 7 15.91% 12 17.39% 31 25.83% 20 15.63% 5 7.04% 0.021 Diabetes mellitus type 1 – 3 4.34% 4 3.33% 4 3.13% 1 1.40% – Diseases of the circulatory system, unspecified 6 13.64% 7 10.14% 18 15.00% 23 17.97% 14 19.72% 0.53 Supraventricular tachyarrhythmia 1 2.27% – – 12 9.38% 10 14.08% – Heart failure, unspecified – 5 7.25% 10 8.33% ( 7 5.47% 8 11.27% – Cerebral ischaemic stroke 2 4.55% 1 1.49% 9 7.50% 9 7.03% 7 9.86% 0.31 Acute ischaemic heart disease, unspecified – 2 2.90% 7 5.83% 4 3.13% 6 8.45% – Chronic kidney disease 4 9.09% 3 4.35% 14 11.67% 14 10.94% 9 12.67% 0.48 Dependence on renal dialysis 3 6.82% 13 13.04% 13 10.83% 15 11.72% 3 4.23% 0.072 Kidney failure, unspecified 3 6.82% 3 4.35% 9 7.50% 8 6.25% 3 4.23% 0.87 Respiratory failure 2 4.55% 4 5.80% ) 3 2.50% 4 3.13% 6 8.45% 0.33 Chronic obstructive pulmonary disease 1 2.27% – 12 10.00% 10 7.81% 2 2.82% – Aggregated group of diseases Diseases of the circulatory system 14 31.82% 33 47.83% 74 61.67% 82 64.06% 53 74.65% 0.0001 Endocrine, nutritional or metabolic diseases 21 47.73% 30 43.48% 63 52.50% 53 41.41% 18 25.35% 0.007 Diseases of the genitourinary system 7 15.91% 7 10.14% 24 20.00% 29 22.66% 12 16.90% 0.27 Diseases of the respiratory system 6 13.64% 10 14.49% 18 15.00% 20 15.63% 13 18.31% 0.96 Diseases of the nervous system 3 6.82% 2 2.90% 18 15.00% 19 14.84% 19 26.76% 0.001 Diseases of the digestive system 4 9.09% 8 11.59% 18 15.00% 10 7.81% 4 5.63% 0.23 Neoplasms 1 2.27% 7 10.14% 12 10.00% 16 12.50% 4 5.63% 0.25 Mental, behavioural or neurodevelopmental disorders 3 6.82% 2 2.90% 6 5.00% 15 11.72% 10 14.08% 0.049 Factors influencing health status or contact with health services 4 9.09% 10 14.49% 15 12.50% 19 14.84% 6 8.45% 0.66 Symptoms, signs or clinical findings, not elsewhere classified 2 4.55% 3 4.35% 6 5.00% 11 8.59% 2 2.82% 0.46 Certain infectious or parasitic diseases 3 6.82% 2 2.90% 6 5.00% 4 3.13% 5 7.04% 0.63 In bold the p-values that reached the statistical significance.
table-wrap	Table 2 Association between age and distinct diseases or group of diseases, as absolute numbers and as percentage. Characteristicsn, % < 504410.19% 50–596915.97% 60–6912027.78% 70–7912829.93% ≥ 807116.44% p-value Gender Male 29 66% 56 81.16% 73 60.83% 80 62.5% 44 61.97% 0.05 Female 15 34.09% 13 19.84% 47 39.17% 48 37.5% 27 38.03% Distinct diseases Hypertension 12 27.28% 21 30.43% 49 40.83% 45 35.16% 35 49.30% 0.074 Obesity 13 29.55% 9 13.04% 20 16.67% 8 6.25% 3 4.23% 0.0001 Diabetes Total 14 31.82% 23 33.33% 53 44.17% 48 37.5% 16 22.53% 0.045 Diabetes mellitus, unspecified 7 15.91% 8 11.59% 18 15.00% 24 18.75% 10 14.08% 0.74 Diabetes mellitus type 2 7 15.91% 12 17.39% 31 25.83% 20 15.63% 5 7.04% 0.021 Diabetes mellitus type 1 – 3 4.34% 4 3.33% 4 3.13% 1 1.40% – Diseases of the circulatory system, unspecified 6 13.64% 7 10.14% 18 15.00% 23 17.97% 14 19.72% 0.53 Supraventricular tachyarrhythmia 1 2.27% – – 12 9.38% 10 14.08% – Heart failure, unspecified – 5 7.25% 10 8.33% ( 7 5.47% 8 11.27% – Cerebral ischaemic stroke 2 4.55% 1 1.49% 9 7.50% 9 7.03% 7 9.86% 0.31 Acute ischaemic heart disease, unspecified – 2 2.90% 7 5.83% 4 3.13% 6 8.45% – Chronic kidney disease 4 9.09% 3 4.35% 14 11.67% 14 10.94% 9 12.67% 0.48 Dependence on renal dialysis 3 6.82% 13 13.04% 13 10.83% 15 11.72% 3 4.23% 0.072 Kidney failure, unspecified 3 6.82% 3 4.35% 9 7.50% 8 6.25% 3 4.23% 0.87 Respiratory failure 2 4.55% 4 5.80% ) 3 2.50% 4 3.13% 6 8.45% 0.33 Chronic obstructive pulmonary disease 1 2.27% – 12 10.00% 10 7.81% 2 2.82% – Aggregated group of diseases Diseases of the circulatory system 14 31.82% 33 47.83% 74 61.67% 82 64.06% 53 74.65% 0.0001 Endocrine, nutritional or metabolic diseases 21 47.73% 30 43.48% 63 52.50% 53 41.41% 18 25.35% 0.007 Diseases of the genitourinary system 7 15.91% 7 10.14% 24 20.00% 29 22.66% 12 16.90% 0.27 Diseases of the respiratory system 6 13.64% 10 14.49% 18 15.00% 20 15.63% 13 18.31% 0.96 Diseases of the nervous system 3 6.82% 2 2.90% 18 15.00% 19 14.84% 19 26.76% 0.001 Diseases of the digestive system 4 9.09% 8 11.59% 18 15.00% 10 7.81% 4 5.63% 0.23 Neoplasms 1 2.27% 7 10.14% 12 10.00% 16 12.50% 4 5.63% 0.25 Mental, behavioural or neurodevelopmental disorders 3 6.82% 2 2.90% 6 5.00% 15 11.72% 10 14.08% 0.049 Factors influencing health status or contact with health services 4 9.09% 10 14.49% 15 12.50% 19 14.84% 6 8.45% 0.66 Symptoms, signs or clinical findings, not elsewhere classified 2 4.55% 3 4.35% 6 5.00% 11 8.59% 2 2.82% 0.46 Certain infectious or parasitic diseases 3 6.82% 2 2.90% 6 5.00% 4 3.13% 5 7.04% 0.63 In bold the p-values that reached the statistical significance.
label	Table 2
caption	Association between age and distinct diseases or group of diseases, as absolute numbers and as percentage.
p	Association between age and distinct diseases or group of diseases, as absolute numbers and as percentage.
table	Characteristicsn, % < 504410.19% 50–596915.97% 60–6912027.78% 70–7912829.93% ≥ 807116.44% p-value Gender Male 29 66% 56 81.16% 73 60.83% 80 62.5% 44 61.97% 0.05 Female 15 34.09% 13 19.84% 47 39.17% 48 37.5% 27 38.03% Distinct diseases Hypertension 12 27.28% 21 30.43% 49 40.83% 45 35.16% 35 49.30% 0.074 Obesity 13 29.55% 9 13.04% 20 16.67% 8 6.25% 3 4.23% 0.0001 Diabetes Total 14 31.82% 23 33.33% 53 44.17% 48 37.5% 16 22.53% 0.045 Diabetes mellitus, unspecified 7 15.91% 8 11.59% 18 15.00% 24 18.75% 10 14.08% 0.74 Diabetes mellitus type 2 7 15.91% 12 17.39% 31 25.83% 20 15.63% 5 7.04% 0.021 Diabetes mellitus type 1 – 3 4.34% 4 3.33% 4 3.13% 1 1.40% – Diseases of the circulatory system, unspecified 6 13.64% 7 10.14% 18 15.00% 23 17.97% 14 19.72% 0.53 Supraventricular tachyarrhythmia 1 2.27% – – 12 9.38% 10 14.08% – Heart failure, unspecified – 5 7.25% 10 8.33% ( 7 5.47% 8 11.27% – Cerebral ischaemic stroke 2 4.55% 1 1.49% 9 7.50% 9 7.03% 7 9.86% 0.31 Acute ischaemic heart disease, unspecified – 2 2.90% 7 5.83% 4 3.13% 6 8.45% – Chronic kidney disease 4 9.09% 3 4.35% 14 11.67% 14 10.94% 9 12.67% 0.48 Dependence on renal dialysis 3 6.82% 13 13.04% 13 10.83% 15 11.72% 3 4.23% 0.072 Kidney failure, unspecified 3 6.82% 3 4.35% 9 7.50% 8 6.25% 3 4.23% 0.87 Respiratory failure 2 4.55% 4 5.80% ) 3 2.50% 4 3.13% 6 8.45% 0.33 Chronic obstructive pulmonary disease 1 2.27% – 12 10.00% 10 7.81% 2 2.82% – Aggregated group of diseases Diseases of the circulatory system 14 31.82% 33 47.83% 74 61.67% 82 64.06% 53 74.65% 0.0001 Endocrine, nutritional or metabolic diseases 21 47.73% 30 43.48% 63 52.50% 53 41.41% 18 25.35% 0.007 Diseases of the genitourinary system 7 15.91% 7 10.14% 24 20.00% 29 22.66% 12 16.90% 0.27 Diseases of the respiratory system 6 13.64% 10 14.49% 18 15.00% 20 15.63% 13 18.31% 0.96 Diseases of the nervous system 3 6.82% 2 2.90% 18 15.00% 19 14.84% 19 26.76% 0.001 Diseases of the digestive system 4 9.09% 8 11.59% 18 15.00% 10 7.81% 4 5.63% 0.23 Neoplasms 1 2.27% 7 10.14% 12 10.00% 16 12.50% 4 5.63% 0.25 Mental, behavioural or neurodevelopmental disorders 3 6.82% 2 2.90% 6 5.00% 15 11.72% 10 14.08% 0.049 Factors influencing health status or contact with health services 4 9.09% 10 14.49% 15 12.50% 19 14.84% 6 8.45% 0.66 Symptoms, signs or clinical findings, not elsewhere classified 2 4.55% 3 4.35% 6 5.00% 11 8.59% 2 2.82% 0.46 Certain infectious or parasitic diseases 3 6.82% 2 2.90% 6 5.00% 4 3.13% 5 7.04% 0.63
tr	Characteristicsn, % < 504410.19% 50–596915.97% 60–6912027.78% 70–7912829.93% ≥ 807116.44% p-value
th	Characteristicsn, %
th	< 504410.19%
th	50–596915.97%
th	60–6912027.78%
th	70–7912829.93%
th	≥ 807116.44%
th	p-value
tr	Gender
td	Gender
tr	Male 29 66% 56 81.16% 73 60.83% 80 62.5% 44 61.97% 0.05
td	Male
td	29 66%
p	29
p	66%
td	56 81.16%
p	56
p	81.16%
td	73 60.83%
p	73
p	60.83%
td	80 62.5%
p	80
p	62.5%
td	44 61.97%
p	44
p	61.97%
td	0.05
tr	Female 15 34.09% 13 19.84% 47 39.17% 48 37.5% 27 38.03%
td	Female
td	15 34.09%
p	15
p	34.09%
td	13 19.84%
p	13
p	19.84%
td	47 39.17%
p	47
p	39.17%
td	48 37.5%
p	48
p	37.5%
td	27 38.03%
p	27
p	38.03%
tr	Distinct diseases
td	Distinct diseases
tr	Hypertension 12 27.28% 21 30.43% 49 40.83% 45 35.16% 35 49.30% 0.074
td	Hypertension
td	12 27.28%
p	12
p	27.28%
td	21 30.43%
p	21
p	30.43%
td	49 40.83%
p	49
p	40.83%
td	45 35.16%
p	45
p	35.16%
td	35 49.30%
p	35
p	49.30%
td	0.074
tr	Obesity 13 29.55% 9 13.04% 20 16.67% 8 6.25% 3 4.23% 0.0001
td	Obesity
td	13 29.55%
p	13
p	29.55%
td	9 13.04%
p	9
p	13.04%
td	20 16.67%
p	20
p	16.67%
td	8 6.25%
p	8
p	6.25%
td	3 4.23%
p	3
p	4.23%
td	0.0001
tr	Diabetes
td	Diabetes
tr	Total 14 31.82% 23 33.33% 53 44.17% 48 37.5% 16 22.53% 0.045
td	Total
td	14 31.82%
p	14
p	31.82%
td	23 33.33%
p	23
p	33.33%
td	53 44.17%
p	53
p	44.17%
td	48 37.5%
p	48
p	37.5%
td	16 22.53%
p	16
p	22.53%
td	0.045
tr	Diabetes mellitus, unspecified 7 15.91% 8 11.59% 18 15.00% 24 18.75% 10 14.08% 0.74
td	Diabetes mellitus, unspecified
td	7 15.91%
p	7
p	15.91%
td	8 11.59%
p	8
p	11.59%
td	18 15.00%
p	18
p	15.00%
td	24 18.75%
p	24
p	18.75%
td	10 14.08%
p	10
p	14.08%
td	0.74
tr	Diabetes mellitus type 2 7 15.91% 12 17.39% 31 25.83% 20 15.63% 5 7.04% 0.021
td	Diabetes mellitus type 2
td	7 15.91%
p	7
p	15.91%
td	12 17.39%
p	12
p	17.39%
td	31 25.83%
p	31
p	25.83%
td	20 15.63%
p	20
p	15.63%
td	5 7.04%
p	5
p	7.04%
td	0.021
tr	Diabetes mellitus type 1 – 3 4.34% 4 3.33% 4 3.13% 1 1.40% –
td	Diabetes mellitus type 1
td	–
td	3 4.34%
p	3
p	4.34%
td	4 3.33%
p	4
p	3.33%
td	4 3.13%
p	4
p	3.13%
td	1 1.40%
p	1
p	1.40%
td	–
tr	Diseases of the circulatory system, unspecified 6 13.64% 7 10.14% 18 15.00% 23 17.97% 14 19.72% 0.53
td	Diseases of the circulatory system, unspecified
td	6 13.64%
p	6
p	13.64%
td	7 10.14%
p	7
p	10.14%
td	18 15.00%
p	18
p	15.00%
td	23 17.97%
p	23
p	17.97%
td	14 19.72%
p	14
p	19.72%
td	0.53
tr	Supraventricular tachyarrhythmia 1 2.27% – – 12 9.38% 10 14.08% –
td	Supraventricular tachyarrhythmia
td	1 2.27%
p	1
p	2.27%
td	–
td	–
td	12 9.38%
p	12
p	9.38%
td	10 14.08%
p	10
p	14.08%
td	–
tr	Heart failure, unspecified – 5 7.25% 10 8.33% ( 7 5.47% 8 11.27% –
td	Heart failure, unspecified
td	–
td	5 7.25%
p	5
p	7.25%
td	10 8.33% (
p	10
p	8.33%
p	(
td	7 5.47%
p	7
p	5.47%
td	8 11.27%
p	8
p	11.27%
td	–
tr	Cerebral ischaemic stroke 2 4.55% 1 1.49% 9 7.50% 9 7.03% 7 9.86% 0.31
td	Cerebral ischaemic stroke
td	2 4.55%
p	2
p	4.55%
td	1 1.49%
p	1
p	1.49%
td	9 7.50%
p	9
p	7.50%
td	9 7.03%
p	9
p	7.03%
td	7 9.86%
p	7
p	9.86%
td	0.31
tr	Acute ischaemic heart disease, unspecified – 2 2.90% 7 5.83% 4 3.13% 6 8.45% –
td	Acute ischaemic heart disease, unspecified
td	–
td	2 2.90%
p	2
p	2.90%
td	7 5.83%
p	7
p	5.83%
td	4 3.13%
p	4
p	3.13%
td	6 8.45%
p	6
p	8.45%
td	–
tr	Chronic kidney disease 4 9.09% 3 4.35% 14 11.67% 14 10.94% 9 12.67% 0.48
td	Chronic kidney disease
td	4 9.09%
p	4
p	9.09%
td	3 4.35%
p	3
p	4.35%
td	14 11.67%
p	14
p	11.67%
td	14 10.94%
p	14
p	10.94%
td	9 12.67%
p	9
p	12.67%
td	0.48
tr	Dependence on renal dialysis 3 6.82% 13 13.04% 13 10.83% 15 11.72% 3 4.23% 0.072
td	Dependence on renal dialysis
td	3 6.82%
p	3
p	6.82%
td	13 13.04%
p	13
p	13.04%
td	13 10.83%
p	13
p	10.83%
td	15 11.72%
p	15
p	11.72%
td	3 4.23%
p	3
p	4.23%
td	0.072
tr	Kidney failure, unspecified 3 6.82% 3 4.35% 9 7.50% 8 6.25% 3 4.23% 0.87
td	Kidney failure, unspecified
td	3 6.82%
p	3
p	6.82%
td	3 4.35%
p	3
p	4.35%
td	9 7.50%
p	9
p	7.50%
td	8 6.25%
p	8
p	6.25%
td	3 4.23%
p	3
p	4.23%
td	0.87
tr	Respiratory failure 2 4.55% 4 5.80% ) 3 2.50% 4 3.13% 6 8.45% 0.33
td	Respiratory failure
td	2 4.55%
p	2
p	4.55%
td	4 5.80% )
p	4
p	5.80%
p	)
td	3 2.50%
p	3
p	2.50%
td	4 3.13%
p	4
p	3.13%
td	6 8.45%
p	6
p	8.45%
td	0.33
tr	Chronic obstructive pulmonary disease 1 2.27% – 12 10.00% 10 7.81% 2 2.82% –
td	Chronic obstructive pulmonary disease
td	1 2.27%
p	1
p	2.27%
td	–
td	12 10.00%
p	12
p	10.00%
td	10 7.81%
p	10
p	7.81%
td	2 2.82%
p	2
p	2.82%
td	–
tr	Aggregated group of diseases
td	Aggregated group of diseases
tr	Diseases of the circulatory system 14 31.82% 33 47.83% 74 61.67% 82 64.06% 53 74.65% 0.0001
td	Diseases of the circulatory system
td	14 31.82%
p	14
p	31.82%
td	33 47.83%
p	33
p	47.83%
td	74 61.67%
p	74
p	61.67%
td	82 64.06%
p	82
p	64.06%
td	53 74.65%
p	53
p	74.65%
td	0.0001
tr	Endocrine, nutritional or metabolic diseases 21 47.73% 30 43.48% 63 52.50% 53 41.41% 18 25.35% 0.007
td	Endocrine, nutritional or metabolic diseases
td	21 47.73%
p	21
p	47.73%
td	30 43.48%
p	30
p	43.48%
td	63 52.50%
p	63
p	52.50%
td	53 41.41%
p	53
p	41.41%
td	18 25.35%
p	18
p	25.35%
td	0.007
tr	Diseases of the genitourinary system 7 15.91% 7 10.14% 24 20.00% 29 22.66% 12 16.90% 0.27
td	Diseases of the genitourinary system
td	7 15.91%
p	7
p	15.91%
td	7 10.14%
p	7
p	10.14%
td	24 20.00%
p	24
p	20.00%
td	29 22.66%
p	29
p	22.66%
td	12 16.90%
p	12
p	16.90%
td	0.27
tr	Diseases of the respiratory system 6 13.64% 10 14.49% 18 15.00% 20 15.63% 13 18.31% 0.96
td	Diseases of the respiratory system
td	6 13.64%
p	6
p	13.64%
td	10 14.49%
p	10
p	14.49%
td	18 15.00%
p	18
p	15.00%
td	20 15.63%
p	20
p	15.63%
td	13 18.31%
p	13
p	18.31%
td	0.96
tr	Diseases of the nervous system 3 6.82% 2 2.90% 18 15.00% 19 14.84% 19 26.76% 0.001
td	Diseases of the nervous system
td	3 6.82%
p	3
p	6.82%
td	2 2.90%
p	2
p	2.90%
td	18 15.00%
p	18
p	15.00%
td	19 14.84%
p	19
p	14.84%
td	19 26.76%
p	19
p	26.76%
td	0.001
tr	Diseases of the digestive system 4 9.09% 8 11.59% 18 15.00% 10 7.81% 4 5.63% 0.23
td	Diseases of the digestive system
td	4 9.09%
p	4
p	9.09%
td	8 11.59%
p	8
p	11.59%
td	18 15.00%
p	18
p	15.00%
td	10 7.81%
p	10
p	7.81%
td	4 5.63%
p	4
p	5.63%
td	0.23
tr	Neoplasms 1 2.27% 7 10.14% 12 10.00% 16 12.50% 4 5.63% 0.25
td	Neoplasms
td	1 2.27%
p	1
p	2.27%
td	7 10.14%
p	7
p	10.14%
td	12 10.00%
p	12
p	10.00%
td	16 12.50%
p	16
p	12.50%
td	4 5.63%
p	4
p	5.63%
td	0.25
tr	Mental, behavioural or neurodevelopmental disorders 3 6.82% 2 2.90% 6 5.00% 15 11.72% 10 14.08% 0.049
td	Mental, behavioural or neurodevelopmental disorders
td	3 6.82%
p	3
p	6.82%
td	2 2.90%
p	2
p	2.90%
td	6 5.00%
p	6
p	5.00%
td	15 11.72%
p	15
p	11.72%
td	10 14.08%
p	10
p	14.08%
td	0.049
tr	Factors influencing health status or contact with health services 4 9.09% 10 14.49% 15 12.50% 19 14.84% 6 8.45% 0.66
td	Factors influencing health status or contact with health services
td	4 9.09%
p	4
p	9.09%
td	10 14.49%
p	10
p	14.49%
td	15 12.50%
p	15
p	12.50%
td	19 14.84%
p	19
p	14.84%
td	6 8.45%
p	6
p	8.45%
td	0.66
tr	Symptoms, signs or clinical findings, not elsewhere classified 2 4.55% 3 4.35% 6 5.00% 11 8.59% 2 2.82% 0.46
td	Symptoms, signs or clinical findings, not elsewhere classified
td	2 4.55%
p	2
p	4.55%
td	3 4.35%
p	3
p	4.35%
td	6 5.00%
p	6
p	5.00%
td	11 8.59%
p	11
p	8.59%
td	2 2.82%
p	2
p	2.82%
td	0.46
tr	Certain infectious or parasitic diseases 3 6.82% 2 2.90% 6 5.00% 4 3.13% 5 7.04% 0.63
td	Certain infectious or parasitic diseases
td	3 6.82%
p	3
p	6.82%
td	2 2.90%
p	2
p	2.90%
td	6 5.00%
p	6
p	5.00%
td	4 3.13%
p	4
p	3.13%
td	5 7.04%
p	5
p	7.04%
td	0.63
table-wrap-foot	In bold the p-values that reached the statistical significance.
p	In bold the p-values that reached the statistical significance.
sec	Comparison between COVID-19 and pneumonia fatalities The comparison with pneumonia fatalities yielded that the proportion of males in the COVID-19 death sample is also significantly higher than in the pneumonia death sample as shown in Table 3. Age and the relative frequency of essential hypertension, diseases of the circulatory system, obesity, dependence of renal dialysis, and chronic kidney failure are also different between the COVID-19 and pneumonia groups. Additionally, the odd to die for COVID-19 versus hospital pneumonia is 150% greater for men than women with similar comorbidities (see Table 4). Moreover, the adjusted odds to die for COVID-19 compared to pneumonia are 22 times higher for persons with the chronic kidney failure, doubled for persons with hypertension, and almost equal for patients with type 2 diabetes. Table 3 Comparison between the characteristics of the COVID-19 and pneumonia groups. Characteristics COVID-19 Pneumonia p-value Sample size n, % 432 100% 711 100% Gender n, % 0.0001 Female 150 35% 320 45% Male 282 65% 391 54.99% Age Mean (sd) Median 66.97 (13.07) 68 73.14 (13.68) 78 Age group n, % 0.0001 < 50 44 10.19% 49 6.89% 50–59 69 15.97% 51 7.17% 60–69 120 27.78% 141 19.83% 70–79 128 29.63% 199 27.99% ≥ 80 71 16.44% 271 38.12% Comorbidities Hypertension 162 37.50% 181 25.46% 0.0001 Obesity 56 12.96% 63 8.86% 0.03 Diabetes mellitus type I 12 2.78% 5 0.70% 0.01 Diabetes mellitus type II 75 17.36% 145 20.39% 0.21 Dialysis 43 9.95% 4 0.56% 0.0001 Chronic kidney disease 43 9.95% 88 12.38% 0.21 Supraventricular tachyarrhythmia 25 5.79% 201 28.27% 0.0001 Congestive heart failure 12 2.78% 204 28.69% 0.0001 COPD 24 5.56% 72 10.13% 0.01 Heart failure 29 6.71% 27 3.80% 0.03 Kidney failure 26 6.02% 4 0.56% 0.0001 Respiratory failure 19 4.40% 88 12.38% 0.0001 Cerebral ischaemic stroke 28 6.48% – 0.0001 In bold the p-values that reached the statistical significance. Table 4 Logistic regression analysis for the comparison between COVID-19 and pneumonia death groups. Covariates CoefficientΒ OR (COVID-19 vs. pneumonia) 95% CI p-value Age − 0.022 0.978 [0.968–0.988] 0.0001 Gender (0 = female) 0.414 1.513 [1.145–2.000] 0.01 Hypertension 0.740 2.096 [1.560–2.818] 0.0001 Obesity 0.266 1.305 [0.843–2.019] 0.23 Diabetes type 2 − 0.355 0.701 [0.494–0.996] 0.049 Dialysis 3.104 22.290 [7.244–68.590] 0.0001 Chronic kidney disease − 0.522 0.593 [0.369–0.955] 0.03 Heart failure 1.095 2.990 [1.612–5.546] 0.0001 COPD − 0.730 0.482 [0.284–0.817] 0.01 Supraventricular tachyarrhythmia − 1.687 0.185 [0.116–0.295] 0.0001 Respiratory failure − 1.045 0.352 [0.203–0.609] 0.0001 Intercept 0.884 2.420 0.02 In bold the p-values that reached the statistical significance.
title	Comparison between COVID-19 and pneumonia fatalities
p	The comparison with pneumonia fatalities yielded that the proportion of males in the COVID-19 death sample is also significantly higher than in the pneumonia death sample as shown in Table 3. Age and the relative frequency of essential hypertension, diseases of the circulatory system, obesity, dependence of renal dialysis, and chronic kidney failure are also different between the COVID-19 and pneumonia groups. Additionally, the odd to die for COVID-19 versus hospital pneumonia is 150% greater for men than women with similar comorbidities (see Table 4). Moreover, the adjusted odds to die for COVID-19 compared to pneumonia are 22 times higher for persons with the chronic kidney failure, doubled for persons with hypertension, and almost equal for patients with type 2 diabetes. Table 3 Comparison between the characteristics of the COVID-19 and pneumonia groups. Characteristics COVID-19 Pneumonia p-value Sample size n, % 432 100% 711 100% Gender n, % 0.0001 Female 150 35% 320 45% Male 282 65% 391 54.99% Age Mean (sd) Median 66.97 (13.07) 68 73.14 (13.68) 78 Age group n, % 0.0001 < 50 44 10.19% 49 6.89% 50–59 69 15.97% 51 7.17% 60–69 120 27.78% 141 19.83% 70–79 128 29.63% 199 27.99% ≥ 80 71 16.44% 271 38.12% Comorbidities Hypertension 162 37.50% 181 25.46% 0.0001 Obesity 56 12.96% 63 8.86% 0.03 Diabetes mellitus type I 12 2.78% 5 0.70% 0.01 Diabetes mellitus type II 75 17.36% 145 20.39% 0.21 Dialysis 43 9.95% 4 0.56% 0.0001 Chronic kidney disease 43 9.95% 88 12.38% 0.21 Supraventricular tachyarrhythmia 25 5.79% 201 28.27% 0.0001 Congestive heart failure 12 2.78% 204 28.69% 0.0001 COPD 24 5.56% 72 10.13% 0.01 Heart failure 29 6.71% 27 3.80% 0.03 Kidney failure 26 6.02% 4 0.56% 0.0001 Respiratory failure 19 4.40% 88 12.38% 0.0001 Cerebral ischaemic stroke 28 6.48% – 0.0001 In bold the p-values that reached the statistical significance. Table 4 Logistic regression analysis for the comparison between COVID-19 and pneumonia death groups. Covariates CoefficientΒ OR (COVID-19 vs. pneumonia) 95% CI p-value Age − 0.022 0.978 [0.968–0.988] 0.0001 Gender (0 = female) 0.414 1.513 [1.145–2.000] 0.01 Hypertension 0.740 2.096 [1.560–2.818] 0.0001 Obesity 0.266 1.305 [0.843–2.019] 0.23 Diabetes type 2 − 0.355 0.701 [0.494–0.996] 0.049 Dialysis 3.104 22.290 [7.244–68.590] 0.0001 Chronic kidney disease − 0.522 0.593 [0.369–0.955] 0.03 Heart failure 1.095 2.990 [1.612–5.546] 0.0001 COPD − 0.730 0.482 [0.284–0.817] 0.01 Supraventricular tachyarrhythmia − 1.687 0.185 [0.116–0.295] 0.0001 Respiratory failure − 1.045 0.352 [0.203–0.609] 0.0001 Intercept 0.884 2.420 0.02 In bold the p-values that reached the statistical significance.
table-wrap	Table 3 Comparison between the characteristics of the COVID-19 and pneumonia groups. Characteristics COVID-19 Pneumonia p-value Sample size n, % 432 100% 711 100% Gender n, % 0.0001 Female 150 35% 320 45% Male 282 65% 391 54.99% Age Mean (sd) Median 66.97 (13.07) 68 73.14 (13.68) 78 Age group n, % 0.0001 < 50 44 10.19% 49 6.89% 50–59 69 15.97% 51 7.17% 60–69 120 27.78% 141 19.83% 70–79 128 29.63% 199 27.99% ≥ 80 71 16.44% 271 38.12% Comorbidities Hypertension 162 37.50% 181 25.46% 0.0001 Obesity 56 12.96% 63 8.86% 0.03 Diabetes mellitus type I 12 2.78% 5 0.70% 0.01 Diabetes mellitus type II 75 17.36% 145 20.39% 0.21 Dialysis 43 9.95% 4 0.56% 0.0001 Chronic kidney disease 43 9.95% 88 12.38% 0.21 Supraventricular tachyarrhythmia 25 5.79% 201 28.27% 0.0001 Congestive heart failure 12 2.78% 204 28.69% 0.0001 COPD 24 5.56% 72 10.13% 0.01 Heart failure 29 6.71% 27 3.80% 0.03 Kidney failure 26 6.02% 4 0.56% 0.0001 Respiratory failure 19 4.40% 88 12.38% 0.0001 Cerebral ischaemic stroke 28 6.48% – 0.0001 In bold the p-values that reached the statistical significance.
label	Table 3
caption	Comparison between the characteristics of the COVID-19 and pneumonia groups.
p	Comparison between the characteristics of the COVID-19 and pneumonia groups.
table	Characteristics COVID-19 Pneumonia p-value Sample size n, % 432 100% 711 100% Gender n, % 0.0001 Female 150 35% 320 45% Male 282 65% 391 54.99% Age Mean (sd) Median 66.97 (13.07) 68 73.14 (13.68) 78 Age group n, % 0.0001 < 50 44 10.19% 49 6.89% 50–59 69 15.97% 51 7.17% 60–69 120 27.78% 141 19.83% 70–79 128 29.63% 199 27.99% ≥ 80 71 16.44% 271 38.12% Comorbidities Hypertension 162 37.50% 181 25.46% 0.0001 Obesity 56 12.96% 63 8.86% 0.03 Diabetes mellitus type I 12 2.78% 5 0.70% 0.01 Diabetes mellitus type II 75 17.36% 145 20.39% 0.21 Dialysis 43 9.95% 4 0.56% 0.0001 Chronic kidney disease 43 9.95% 88 12.38% 0.21 Supraventricular tachyarrhythmia 25 5.79% 201 28.27% 0.0001 Congestive heart failure 12 2.78% 204 28.69% 0.0001 COPD 24 5.56% 72 10.13% 0.01 Heart failure 29 6.71% 27 3.80% 0.03 Kidney failure 26 6.02% 4 0.56% 0.0001 Respiratory failure 19 4.40% 88 12.38% 0.0001 Cerebral ischaemic stroke 28 6.48% – 0.0001
tr	Characteristics COVID-19 Pneumonia p-value
th	Characteristics
th	COVID-19
th	Pneumonia
th	p-value
tr	Sample size n, % 432 100% 711 100%
td	Sample size n, %
p	Sample size
p	n, %
td	432 100%
p	432
p	100%
td	711 100%
p	711
p	100%
tr	Gender n, % 0.0001
td	Gender n, %
p	Gender
p	n, %
td	0.0001
tr	Female 150 35% 320 45%
td	Female
td	150 35%
p	150
p	35%
td	320 45%
p	320
p	45%
tr	Male 282 65% 391 54.99%
td	Male
td	282 65%
p	282
p	65%
td	391 54.99%
p	391
p	54.99%
tr	Age Mean (sd) Median 66.97 (13.07) 68 73.14 (13.68) 78
td	Age Mean (sd) Median
p	Age
p	Mean (sd)
p	Median
td	66.97 (13.07) 68
p	66.97 (13.07)
p	68
td	73.14 (13.68) 78
p	73.14 (13.68)
p	78
tr	Age group n, % 0.0001
td	Age group n, %
p	Age group
p	n, %
td	0.0001
tr	< 50 44 10.19% 49 6.89%
td	< 50
td	44 10.19%
p	44
p	10.19%
td	49 6.89%
p	49
p	6.89%
tr	50–59 69 15.97% 51 7.17%
td	50–59
td	69 15.97%
p	69
p	15.97%
td	51 7.17%
p	51
p	7.17%
tr	60–69 120 27.78% 141 19.83%
td	60–69
td	120 27.78%
p	120
p	27.78%
td	141 19.83%
p	141
p	19.83%
tr	70–79 128 29.63% 199 27.99%
td	70–79
td	128 29.63%
p	128
p	29.63%
td	199 27.99%
p	199
p	27.99%
tr	≥ 80 71 16.44% 271 38.12%
td	≥ 80
td	71 16.44%
p	71
p	16.44%
td	271 38.12%
p	271
p	38.12%
tr	Comorbidities
td	Comorbidities
tr	Hypertension 162 37.50% 181 25.46% 0.0001
td	Hypertension
td	162 37.50%
p	162
p	37.50%
td	181 25.46%
p	181
p	25.46%
td	0.0001
tr	Obesity 56 12.96% 63 8.86% 0.03
td	Obesity
td	56 12.96%
p	56
p	12.96%
td	63 8.86%
p	63
p	8.86%
td	0.03
tr	Diabetes mellitus type I 12 2.78% 5 0.70% 0.01
td	Diabetes mellitus type I
td	12 2.78%
p	12
p	2.78%
td	5 0.70%
p	5
p	0.70%
td	0.01
tr	Diabetes mellitus type II 75 17.36% 145 20.39% 0.21
td	Diabetes mellitus type II
td	75 17.36%
p	75
p	17.36%
td	145 20.39%
p	145
p	20.39%
td	0.21
tr	Dialysis 43 9.95% 4 0.56% 0.0001
td	Dialysis
td	43 9.95%
p	43
p	9.95%
td	4 0.56%
p	4
p	0.56%
td	0.0001
tr	Chronic kidney disease 43 9.95% 88 12.38% 0.21
td	Chronic kidney disease
td	43 9.95%
p	43
p	9.95%
td	88 12.38%
p	88
p	12.38%
td	0.21
tr	Supraventricular tachyarrhythmia 25 5.79% 201 28.27% 0.0001
td	Supraventricular tachyarrhythmia
td	25 5.79%
p	25
p	5.79%
td	201 28.27%
p	201
p	28.27%
td	0.0001
tr	Congestive heart failure 12 2.78% 204 28.69% 0.0001
td	Congestive heart failure
td	12 2.78%
p	12
p	2.78%
td	204 28.69%
p	204
p	28.69%
td	0.0001
tr	COPD 24 5.56% 72 10.13% 0.01
td	COPD
td	24 5.56%
p	24
p	5.56%
td	72 10.13%
p	72
p	10.13%
td	0.01
tr	Heart failure 29 6.71% 27 3.80% 0.03
td	Heart failure
td	29 6.71%
p	29
p	6.71%
td	27 3.80%
p	27
p	3.80%
td	0.03
tr	Kidney failure 26 6.02% 4 0.56% 0.0001
td	Kidney failure
td	26 6.02%
p	26
p	6.02%
td	4 0.56%
p	4
p	0.56%
td	0.0001
tr	Respiratory failure 19 4.40% 88 12.38% 0.0001
td	Respiratory failure
td	19 4.40%
p	19
p	4.40%
td	88 12.38%
p	88
p	12.38%
td	0.0001
tr	Cerebral ischaemic stroke 28 6.48% – 0.0001
td	Cerebral ischaemic stroke
td	28 6.48%
p	28
p	6.48%
td	–
td	0.0001
table-wrap-foot	In bold the p-values that reached the statistical significance.
p	In bold the p-values that reached the statistical significance.
table-wrap	Table 4 Logistic regression analysis for the comparison between COVID-19 and pneumonia death groups. Covariates CoefficientΒ OR (COVID-19 vs. pneumonia) 95% CI p-value Age − 0.022 0.978 [0.968–0.988] 0.0001 Gender (0 = female) 0.414 1.513 [1.145–2.000] 0.01 Hypertension 0.740 2.096 [1.560–2.818] 0.0001 Obesity 0.266 1.305 [0.843–2.019] 0.23 Diabetes type 2 − 0.355 0.701 [0.494–0.996] 0.049 Dialysis 3.104 22.290 [7.244–68.590] 0.0001 Chronic kidney disease − 0.522 0.593 [0.369–0.955] 0.03 Heart failure 1.095 2.990 [1.612–5.546] 0.0001 COPD − 0.730 0.482 [0.284–0.817] 0.01 Supraventricular tachyarrhythmia − 1.687 0.185 [0.116–0.295] 0.0001 Respiratory failure − 1.045 0.352 [0.203–0.609] 0.0001 Intercept 0.884 2.420 0.02 In bold the p-values that reached the statistical significance.
label	Table 4
caption	Logistic regression analysis for the comparison between COVID-19 and pneumonia death groups.
p	Logistic regression analysis for the comparison between COVID-19 and pneumonia death groups.
table	Covariates CoefficientΒ OR (COVID-19 vs. pneumonia) 95% CI p-value Age − 0.022 0.978 [0.968–0.988] 0.0001 Gender (0 = female) 0.414 1.513 [1.145–2.000] 0.01 Hypertension 0.740 2.096 [1.560–2.818] 0.0001 Obesity 0.266 1.305 [0.843–2.019] 0.23 Diabetes type 2 − 0.355 0.701 [0.494–0.996] 0.049 Dialysis 3.104 22.290 [7.244–68.590] 0.0001 Chronic kidney disease − 0.522 0.593 [0.369–0.955] 0.03 Heart failure 1.095 2.990 [1.612–5.546] 0.0001 COPD − 0.730 0.482 [0.284–0.817] 0.01 Supraventricular tachyarrhythmia − 1.687 0.185 [0.116–0.295] 0.0001 Respiratory failure − 1.045 0.352 [0.203–0.609] 0.0001 Intercept 0.884 2.420 0.02
tr	Covariates CoefficientΒ OR (COVID-19 vs. pneumonia) 95% CI p-value
th	Covariates
th	CoefficientΒ
th	OR (COVID-19 vs. pneumonia)
th	95% CI
th	p-value
tr	Age − 0.022 0.978 [0.968–0.988] 0.0001
td	Age
td	− 0.022
td	0.978
td	[0.968–0.988]
td	0.0001
tr	Gender (0 = female) 0.414 1.513 [1.145–2.000] 0.01
td	Gender (0 = female)
td	0.414
td	1.513
td	[1.145–2.000]
td	0.01
tr	Hypertension 0.740 2.096 [1.560–2.818] 0.0001
td	Hypertension
td	0.740
td	2.096
td	[1.560–2.818]
td	0.0001
tr	Obesity 0.266 1.305 [0.843–2.019] 0.23
td	Obesity
td	0.266
td	1.305
td	[0.843–2.019]
td	0.23
tr	Diabetes type 2 − 0.355 0.701 [0.494–0.996] 0.049
td	Diabetes type 2
td	− 0.355
td	0.701
td	[0.494–0.996]
td	0.049
tr	Dialysis 3.104 22.290 [7.244–68.590] 0.0001
td	Dialysis
td	3.104
td	22.290
td	[7.244–68.590]
td	0.0001
tr	Chronic kidney disease − 0.522 0.593 [0.369–0.955] 0.03
td	Chronic kidney disease
td	− 0.522
td	0.593
td	[0.369–0.955]
td	0.03
tr	Heart failure 1.095 2.990 [1.612–5.546] 0.0001
td	Heart failure
td	1.095
td	2.990
td	[1.612–5.546]
td	0.0001
tr	COPD − 0.730 0.482 [0.284–0.817] 0.01
td	COPD
td	− 0.730
td	0.482
td	[0.284–0.817]
td	0.01
tr	Supraventricular tachyarrhythmia − 1.687 0.185 [0.116–0.295] 0.0001
td	Supraventricular tachyarrhythmia
td	− 1.687
td	0.185
td	[0.116–0.295]
td	0.0001
tr	Respiratory failure − 1.045 0.352 [0.203–0.609] 0.0001
td	Respiratory failure
td	− 1.045
td	0.352
td	[0.203–0.609]
td	0.0001
tr	Intercept 0.884 2.420 0.02
td	Intercept
td	0.884
td	2.420
td	0.02
table-wrap-foot	In bold the p-values that reached the statistical significance.
p	In bold the p-values that reached the statistical significance.
sec	The burden of disease and survival probability We also extended the analysis to the coexistence of multiple health conditions, and we focused on the multimorbidity. This has a mean value of 2.109, an SD of 1.66 and a median value of 2 in our study. The multimorbidity factor with levels 0, 1, 2, 3, 4 or more than 5 for single or aggregated comorbidities is significantly associated with the factor age but not with gender, confirming the pattern displayed at the level of single conditions and group of diseases (Tables S3 and S4). The Charlson comorbidity index (CCI) had a mean of 1.324, SD of 0.95 and median of one indicating that the severity of the underlying medical condition in the study was mainly mild. This was confirmed by the high relative frequency of 91.20% for the mild level for the factor CCI. There were no significant differences between males and females or among age groups with respect to CCI (Tables S3 and S4). Moreover, 1-year survival probability averaged out at 81%, and 50% of the study individuals had a prognosis for it of 85%.
title	The burden of disease and survival probability
p	We also extended the analysis to the coexistence of multiple health conditions, and we focused on the multimorbidity. This has a mean value of 2.109, an SD of 1.66 and a median value of 2 in our study. The multimorbidity factor with levels 0, 1, 2, 3, 4 or more than 5 for single or aggregated comorbidities is significantly associated with the factor age but not with gender, confirming the pattern displayed at the level of single conditions and group of diseases (Tables S3 and S4).
p	The Charlson comorbidity index (CCI) had a mean of 1.324, SD of 0.95 and median of one indicating that the severity of the underlying medical condition in the study was mainly mild. This was confirmed by the high relative frequency of 91.20% for the mild level for the factor CCI. There were no significant differences between males and females or among age groups with respect to CCI (Tables S3 and S4). Moreover, 1-year survival probability averaged out at 81%, and 50% of the study individuals had a prognosis for it of 85%.
sec	Co-occurrence patterns The co-occurrence of the diseases is of great importance since the mean multimorbidity of 2.109 indicates that not single but pair of diseases are pre-conditions in our mortality study. Pearson's phi coefficient is a robust measure of the comorbidity association, and its values for diseases in general and separated for men and women are displayed in the Table S5. Different co-occurrence patterns emerge between men and women with the striking lack of the correlation between pneumonia or respiratory failure, respectively stroke and dissociative neurological symptom disorder in the female profile. The largest value of the parameter phi of 0.47 respectively 0.40 for female and 0.59 for male corresponds to kidney failure and dependence to dialysis, while the latter is also correlated with chronic kidney disease as the correlation lies by 0.30 in general and 0.38 in men. They are followed by pneumonia or respiratory failure with phi equal to 0.45 in general and to 0.55 in male and by stroke and dissociative neurological symptom disorder with a correlation coefficient of 0.41 in the whole study and of 0.49 in men. Diabetes mellitus type 2 is correlated with hypertension not only for the whole study (phi = 0.11) but also in the men (phi = 0.10) and women subsamples (phi = 0.11). Hypertension appears mainly together with stroke with phi equal to 0.15, 0.16 respectively 0.14 for the whole study, men and women. Finally, a comorbidity network analysis rounded off our exploratory investigation. This identified 11 clusters in the general study, 9 clusters in male, and 6 clusters in female subgroups, each including two or three diseases reinforcing the previous results about pair of diseases as pre-condition for COVID-19 fatalities (see Fig. 2). Figure 2 Co-occurrence clusters for All Subjects (A), Women (B) and Men (C).
title	Co-occurrence patterns
p	The co-occurrence of the diseases is of great importance since the mean multimorbidity of 2.109 indicates that not single but pair of diseases are pre-conditions in our mortality study. Pearson's phi coefficient is a robust measure of the comorbidity association, and its values for diseases in general and separated for men and women are displayed in the Table S5. Different co-occurrence patterns emerge between men and women with the striking lack of the correlation between pneumonia or respiratory failure, respectively stroke and dissociative neurological symptom disorder in the female profile. The largest value of the parameter phi of 0.47 respectively 0.40 for female and 0.59 for male corresponds to kidney failure and dependence to dialysis, while the latter is also correlated with chronic kidney disease as the correlation lies by 0.30 in general and 0.38 in men. They are followed by pneumonia or respiratory failure with phi equal to 0.45 in general and to 0.55 in male and by stroke and dissociative neurological symptom disorder with a correlation coefficient of 0.41 in the whole study and of 0.49 in men. Diabetes mellitus type 2 is correlated with hypertension not only for the whole study (phi = 0.11) but also in the men (phi = 0.10) and women subsamples (phi = 0.11). Hypertension appears mainly together with stroke with phi equal to 0.15, 0.16 respectively 0.14 for the whole study, men and women.
p	Finally, a comorbidity network analysis rounded off our exploratory investigation. This identified 11 clusters in the general study, 9 clusters in male, and 6 clusters in female subgroups, each including two or three diseases reinforcing the previous results about pair of diseases as pre-condition for COVID-19 fatalities (see Fig. 2). Figure 2 Co-occurrence clusters for All Subjects (A), Women (B) and Men (C).
figure	Figure 2 Co-occurrence clusters for All Subjects (A), Women (B) and Men (C).
label	Figure 2
caption	Co-occurrence clusters for All Subjects (A), Women (B) and Men (C).
p	Co-occurrence clusters for All Subjects (A), Women (B) and Men (C).
sec	Discussion The high relative frequency of men in the whole study (65.3% males vs. 34.7% females, p-value < 0.001) and in all subpopulations with distinct diseases and aggregated group of diseases shows that an important risk factor for death by COVID-19 in Romania is being male. This is consistent with other studies about predictors for severity and mortality due to infection with SARS-CoV 2 in other populations10,18. In addition, the comorbidities with the highest prevalence in our study, such as hypertension (37.5%), diabetes (35.4%), obesity (12.27%) and chronic kidney disease (10.19%) as well as diseases of the circulatory system (59.26%) and nutritional or metabolic disorders (42.82%), had similar frequencies with those reported in Chinese or American studies but contradicts the findings of the Italian CORIST study10,19–21. Indeed, the percentages that we found for the distinct comorbidities are consistent to those reported for China; for instance. Li et al., found that the majority of patients deceased in the hospital (N = 128) were men (73%), and that the most prevalent comorbidities associated with COVID-19 mortality were hypertension (47%) and diabetes (13%), while Zhou et al. reported that, among 54 fatalities for COVID-19, 62% were men, and the most prevalent comorbidities were hypertension (30%), diabetes (19%) and coronary heart disease (8%)22,23. These diseases also display significant differences in the relative frequency of men vs. women in our study in Romanian population: 74.67% men vs. 25.33% women for diabetes mellitus type 2, 69.75% men vs. 30.25% women for hypertension, 70.45% men vs. 29.55% women for chronic kidney disease, as well as 65.62% vs. 34.38% respectively 65.95% vs. 34.95% for diseases of the circulatory system and for nutritional or metabolic disorders. Almost all the comorbidities in our study are not associated with gender but with age, as also previously reported in other populations24,25. We found in the present study the highest number of deaths in the age group between 70–79 years (N = 128, 29.63%), followed by the age group 60–69 years (N = 120, 27.78%) and the majority were men (62.5% men vs. 37.5% women). The Romanian age-specific mortality statistics are similar to the estimates of crude-case fatality from Italy, Spain or Germany at the beginning of April, and in contrast to the United Kingdom, where the death count was higher in 65+ years population cluster with a relative frequency of 44% for persons aged 65–79 years and 46% for the 80+ age group4–8,26. The Romanian age structure of the COVID-10 fatalities is similar to the Chinese and lower than the Italian group9,25. Subjects younger than 50 years are less susceptible to die by COVID-19 with a relative frequency of 10% as reported in other articles, and those over 70 years are almost half of the study participants (46%)27. In the survey by Du et al. it was found that subjects who died for SARS-CoV-2 infection were older than 65 years, and that the main comorbidities related to death were hypertension and CVD28. On the other side, in ten European countries and Canada, fatalities due to COVID-19 for individuals younger than 65 years represented 4.8–9.3% of all. Only 13.0% of the COVID-19 fatalities in the UK and 7.8–23.9% in the USA were younger than 65 years. On the other hand, the prevalence for subjects older than 80 years was 54–69% in Europe, 67% in Canada and 36–63% in US27. Another important meta-analysis made by Verity et al. has shown that 13.4% persons older than 80 with more than one comorbidity died for SARS-CoV-2 and that the estimated case fatality rate was 1.38% for 38 countries29. Also, in the Spanish study by Perez-Tanoira et al., men older than 85 years represented 27% of their cohort (N = 108) and the most prevalent comorbidity was coronary heart disease, with a crude mortality rate of 19%30. Statistics displayed in Table S1 and Table 3 show that gender and age are decisive factors for COVID-19 mortality in our Romanian population. This is of great importance for tailored personalized medicine, where treatment management needs to be adapted to patients’ characteristics, and in the epidemiology, where population stratification is the key factor for the identification of the highest risk subgroups. In the whole world, the incidence of COVID-19 and its mortality rate is higher for people with diabetes and/or obesity. The precise mechanisms of SARS-CoV-2 infection are still not fully elucidated, but it has been suggested that this may be due to the chronic inflammatory state and/or insulin-resistance31,32. In our study, diabetes mellitus was classified as type-1, type-2 and unspecified. The unspecified diabetes is most probably the form of type-2 diabetes, since most of the deceased people are between 50 and 79 years and the incidence of the diabetes of type 2 in Romania is extremely higher than type-1, as also shown in the Predatorr study15. Hence, we can assume that the relative frequency of diabetes type-2 was higher than hypertension for the three clusters of age 50–59, 60–69 and 70–79 years. The association between age and different comorbidities suggest the existence of age-stratified positive and negative markers for COVID-19 mortality. Generally, obesity correlates positively with diabetes30,31, although this is confirmed in our study only in the group of youngest age. with a prevalence of 29.55%. Hence, we can interpret obesity as a negative marker for the severity of COVID-19 infection in combination with age below 50 years. In addition, we found the following negative markers: for the age group 50–59 years diabetes (33.33%), for the age group 60–69 years diabetes (44.17%) and hypertension (35.16%), for the age group 70–79 diabetes (44.17%), hypertension (35.16%) and chronic kidney diseases (17.87%), and for the age group over 80 years hypertension (49%), diabetes (22.53%) and the diseases of the circulatory system (19.72%). Further, the main risk factors for death regarding aggregated groups of diseases were the following: the endocrine diseases (diabetes, obesity) with a prevalence of 47.73% for subjects under 50 years, and diseases of the circulatory system, with a prevalence between 47.83 and 74.65% for subjects in all the other age groups. The prevalence of respiratory failure and chronic pulmonary diseases in our study was almost the lowest compared to the rest of the comorbidities but it was similar with the data reported elsewhere22,23. Diseases of the respiratory system were more frequent in the age group over 80 years (18.31%) and had higher prevalence in men (16.67%) versus women (13.33%). This is somewhat surprising since pneumonia had a higher prevalence in the Romanian population in the last three years. Usually, patients with comorbidities from this group of diseases are more vulnerable to viral respiratory infections33. Moreover, it is known that a symptom for the infection with SARS-CoV 2 is an atypical pneumonia with a variable degree of severity34. Since pneumonia could be a negative prognostic factor for COVID-19 patients35, we investigated the relation between the comorbidities of the two diseases. The results show significant differences between these two diseases as shown in Table 4. The hypothesis about the diseases of the respiratory system and the neoplasms as prognosis factors for the mortality through COVID-19 is not confirmed in our study. Therefore, it is likely that patients died because of the impact and virulence of the SARS-CoV-2 infection over comorbidities like hypertension, diabetes mellitus or obesity, in combination with age and gender36. Surprisingly, the average multimorbidity for the COVID-19 fatalities is low but confirms previous findings about infected populations or groups with a severe form of the disease19. The high relative frequency of the mild CCI corroborates this conclusion regarding subjects with a low burden of disease. Our prognosis regarding the survival probability confirms other observations about the estimated year's life lost by the COVID-19 fatalities being "considerably more than the "1–2 years" as presented in37. This is the first article where the co-occurrence profiles are investigated with the help of Pearson's co-occurrence coefficient and where co-occurrence patterns are explicitly described through clustering. Our findings on clustering data are not surprising since diabetes appears mainly with cardiovascular diseases, and kidney diseases with dialysis. Also, the sparsity of respiratory diseases, particularly in the women subgroup, should be further investigated. In conclusion, our study characterizes comorbidities of COVID-19 fatalities in Romania and their association with gender and age. Most of our results are consistent with previous studies from other populations, and the main differences in our Romanian population lie in the lower age of the deceased patients, on their small multimorbidity count and their high 1-year survival probability. In addition, the comparison with the pneumonia group offers a new perspective on the medical pre-conditions of SARS-CoV 2. Fortunately, we registered a low incidence of COVID-19 deaths in Romania until 20th of April 2020, and future studies are needed to determine precisely the pathogenic mechanism that affected these patients, especially males with hypertension, diabetes, chronic kidney disease or obesity and potential protective conditions like those suggested in38–41.
title	Discussion
p	The high relative frequency of men in the whole study (65.3% males vs. 34.7% females, p-value < 0.001) and in all subpopulations with distinct diseases and aggregated group of diseases shows that an important risk factor for death by COVID-19 in Romania is being male. This is consistent with other studies about predictors for severity and mortality due to infection with SARS-CoV 2 in other populations10,18.
p	In addition, the comorbidities with the highest prevalence in our study, such as hypertension (37.5%), diabetes (35.4%), obesity (12.27%) and chronic kidney disease (10.19%) as well as diseases of the circulatory system (59.26%) and nutritional or metabolic disorders (42.82%), had similar frequencies with those reported in Chinese or American studies but contradicts the findings of the Italian CORIST study10,19–21. Indeed, the percentages that we found for the distinct comorbidities are consistent to those reported for China; for instance. Li et al., found that the majority of patients deceased in the hospital (N = 128) were men (73%), and that the most prevalent comorbidities associated with COVID-19 mortality were hypertension (47%) and diabetes (13%), while Zhou et al. reported that, among 54 fatalities for COVID-19, 62% were men, and the most prevalent comorbidities were hypertension (30%), diabetes (19%) and coronary heart disease (8%)22,23. These diseases also display significant differences in the relative frequency of men vs. women in our study in Romanian population: 74.67% men vs. 25.33% women for diabetes mellitus type 2, 69.75% men vs. 30.25% women for hypertension, 70.45% men vs. 29.55% women for chronic kidney disease, as well as 65.62% vs. 34.38% respectively 65.95% vs. 34.95% for diseases of the circulatory system and for nutritional or metabolic disorders. Almost all the comorbidities in our study are not associated with gender but with age, as also previously reported in other populations24,25.
p	We found in the present study the highest number of deaths in the age group between 70–79 years (N = 128, 29.63%), followed by the age group 60–69 years (N = 120, 27.78%) and the majority were men (62.5% men vs. 37.5% women). The Romanian age-specific mortality statistics are similar to the estimates of crude-case fatality from Italy, Spain or Germany at the beginning of April, and in contrast to the United Kingdom, where the death count was higher in 65+ years population cluster with a relative frequency of 44% for persons aged 65–79 years and 46% for the 80+ age group4–8,26. The Romanian age structure of the COVID-10 fatalities is similar to the Chinese and lower than the Italian group9,25. Subjects younger than 50 years are less susceptible to die by COVID-19 with a relative frequency of 10% as reported in other articles, and those over 70 years are almost half of the study participants (46%)27. In the survey by Du et al. it was found that subjects who died for SARS-CoV-2 infection were older than 65 years, and that the main comorbidities related to death were hypertension and CVD28. On the other side, in ten European countries and Canada, fatalities due to COVID-19 for individuals younger than 65 years represented 4.8–9.3% of all. Only 13.0% of the COVID-19 fatalities in the UK and 7.8–23.9% in the USA were younger than 65 years. On the other hand, the prevalence for subjects older than 80 years was 54–69% in Europe, 67% in Canada and 36–63% in US27. Another important meta-analysis made by Verity et al. has shown that 13.4% persons older than 80 with more than one comorbidity died for SARS-CoV-2 and that the estimated case fatality rate was 1.38% for 38 countries29. Also, in the Spanish study by Perez-Tanoira et al., men older than 85 years represented 27% of their cohort (N = 108) and the most prevalent comorbidity was coronary heart disease, with a crude mortality rate of 19%30.
p	Statistics displayed in Table S1 and Table 3 show that gender and age are decisive factors for COVID-19 mortality in our Romanian population. This is of great importance for tailored personalized medicine, where treatment management needs to be adapted to patients’ characteristics, and in the epidemiology, where population stratification is the key factor for the identification of the highest risk subgroups.
p	In the whole world, the incidence of COVID-19 and its mortality rate is higher for people with diabetes and/or obesity. The precise mechanisms of SARS-CoV-2 infection are still not fully elucidated, but it has been suggested that this may be due to the chronic inflammatory state and/or insulin-resistance31,32. In our study, diabetes mellitus was classified as type-1, type-2 and unspecified. The unspecified diabetes is most probably the form of type-2 diabetes, since most of the deceased people are between 50 and 79 years and the incidence of the diabetes of type 2 in Romania is extremely higher than type-1, as also shown in the Predatorr study15. Hence, we can assume that the relative frequency of diabetes type-2 was higher than hypertension for the three clusters of age 50–59, 60–69 and 70–79 years.
p	The association between age and different comorbidities suggest the existence of age-stratified positive and negative markers for COVID-19 mortality. Generally, obesity correlates positively with diabetes30,31, although this is confirmed in our study only in the group of youngest age. with a prevalence of 29.55%. Hence, we can interpret obesity as a negative marker for the severity of COVID-19 infection in combination with age below 50 years. In addition, we found the following negative markers: for the age group 50–59 years diabetes (33.33%), for the age group 60–69 years diabetes (44.17%) and hypertension (35.16%), for the age group 70–79 diabetes (44.17%), hypertension (35.16%) and chronic kidney diseases (17.87%), and for the age group over 80 years hypertension (49%), diabetes (22.53%) and the diseases of the circulatory system (19.72%). Further, the main risk factors for death regarding aggregated groups of diseases were the following: the endocrine diseases (diabetes, obesity) with a prevalence of 47.73% for subjects under 50 years, and diseases of the circulatory system, with a prevalence between 47.83 and 74.65% for subjects in all the other age groups.
p	The prevalence of respiratory failure and chronic pulmonary diseases in our study was almost the lowest compared to the rest of the comorbidities but it was similar with the data reported elsewhere22,23. Diseases of the respiratory system were more frequent in the age group over 80 years (18.31%) and had higher prevalence in men (16.67%) versus women (13.33%). This is somewhat surprising since pneumonia had a higher prevalence in the Romanian population in the last three years. Usually, patients with comorbidities from this group of diseases are more vulnerable to viral respiratory infections33. Moreover, it is known that a symptom for the infection with SARS-CoV 2 is an atypical pneumonia with a variable degree of severity34. Since pneumonia could be a negative prognostic factor for COVID-19 patients35, we investigated the relation between the comorbidities of the two diseases. The results show significant differences between these two diseases as shown in Table 4. The hypothesis about the diseases of the respiratory system and the neoplasms as prognosis factors for the mortality through COVID-19 is not confirmed in our study. Therefore, it is likely that patients died because of the impact and virulence of the SARS-CoV-2 infection over comorbidities like hypertension, diabetes mellitus or obesity, in combination with age and gender36.
p	Surprisingly, the average multimorbidity for the COVID-19 fatalities is low but confirms previous findings about infected populations or groups with a severe form of the disease19.
p	The high relative frequency of the mild CCI corroborates this conclusion regarding subjects with a low burden of disease. Our prognosis regarding the survival probability confirms other observations about the estimated year's life lost by the COVID-19 fatalities being "considerably more than the "1–2 years" as presented in37.
p	This is the first article where the co-occurrence profiles are investigated with the help of Pearson's co-occurrence coefficient and where co-occurrence patterns are explicitly described through clustering. Our findings on clustering data are not surprising since diabetes appears mainly with cardiovascular diseases, and kidney diseases with dialysis. Also, the sparsity of respiratory diseases, particularly in the women subgroup, should be further investigated.
p	In conclusion, our study characterizes comorbidities of COVID-19 fatalities in Romania and their association with gender and age. Most of our results are consistent with previous studies from other populations, and the main differences in our Romanian population lie in the lower age of the deceased patients, on their small multimorbidity count and their high 1-year survival probability. In addition, the comparison with the pneumonia group offers a new perspective on the medical pre-conditions of SARS-CoV 2.
p	Fortunately, we registered a low incidence of COVID-19 deaths in Romania until 20th of April 2020, and future studies are needed to determine precisely the pathogenic mechanism that affected these patients, especially males with hypertension, diabetes, chronic kidney disease or obesity and potential protective conditions like those suggested in38–41.
sec	Methods Study design The cohort The Romanian COVID-19 mortality sub-study was conducted between March 22 and April 20, 2020, with the aim to explore the comorbidity patterns in the patients who died because of this new virus. A multiple-case, multiple-center design was chosen to answer this research question with data collected from the official public notifications of MAI concerning the COVID-19 deaths registered by the INSP3,4,8. Initially, 451 COVID-19 deaths were reported by INSP from the beginning of the epidemic until April 20 (Fig. 1). The list was completed afterwards with 69 additional cases due to the delay in the reporting process from some centers of COVID-19 infection. The official report included information about gender, age, county of residence, history of hospitalization and diagnostics, and/ or comorbidities. Only cases with information about gender, age and comorbidities were included in the study due to the main focus of the research42. Eventually, 432 patients who fulfilled these criteria were kept in the Romanian COVID-19 mortality sub-study. Since we aim at an analytic and not statistical generalization of the results regarding the comorbidity patterns, statistical testing was just serving exploratory purposes42. Hence, the complete-case analysis is justified in this case both by purpose and by the state of things. Case-case study Additionally, we run a case-case study utilizing data regarding the deaths due to hospital pneumonia between March 22 and April 20 in the years 2016–2018 extracted from the administrative database of the INSP. The study was carried out with the approval of the Local Bioethics Council of the INSP (number 6342/04 MAY 2020) and with the permission of the Local Bioethics Council of the “Carol Davila” University of Medicine, Bucharest Romania (9984/11 MAY 2020). Statistical analysis An exploratory data analysis was applied to identify general, spatial-, gender—or age-specific comorbidity patterns in our study. In a first step, the various conditions were sorted according to the standard diagnosis tool ICD-11 (The International Classification of the Diseases) as it was recommended for the mortality and morbidity statistics by the WHO43. The ICD-11 taxonomy includes categories of single diseases or the aggregated groups of diseases. We investigated the comorbidities at both levels of aggregation because the study included a large number of distinct conditions with low prevalence. The statistical analysis was done with the help of the statistical software R version 3.6.2. Continuous variables were expressed as mean, standard deviation (SD) and median. We focused in this paper on comparing the individual crude relative frequencies of the single and aggregated comorbidities between different gender and age groups. Since age is a significant risk factor for mortality, we investigated the frequency of the comorbidities stratified on ten-year intervals starting with the 50+ age groups. Categorical variables were tabulated in frequency (n) and relative frequency (%). A binomial and a chi-squared test for one proportion (only the last was reported) and a chi-squared test for equal proportions were applied when it was appropriate at the 5% level of significance. The post-hoc tests were not corrected for multiple comparisons since the focus of our study was on generating hypotheses for single diseases or group of diseases and not on an experiment-wise approach. The spatial distribution of the COVID-19 mortality corresponded to the country-specific administrative units– counties and was illustrated by single or aggregated disease frequency and relative frequency (%). The association between different comorbidities and the odds ratio between COVID-19 and pneumonia deaths was investigated with the help of the logistic regression, and the results were reported as coefficients β, as point estimates for odds ratio and their 95% confidence intervals. Multimorbidity was computed for each person in the sample since this is a good indicator of the burden of disease44,45 and transformed in a multimorbidity factor with levels 0, 1, 2, 3, 4, or more the 5 for single or aggregated comorbidities. A controversial issue regarding the epidemic of COVID-19 is the burden of disease in the deceased population and the hypothetical survival probability in the case that these people would not have been infected with this virus. A simplified version of the Charlson Comorbidity Index (CCI) was calculated on a scale from 0 to 8, as indicated in24 for an assessment of the individual disease severity. We introduced the CCI factor with two levels (mild for CCI < 3 and severe for CCI ≥ 3), and we tested its differences between gender and age groups by a chi-squared test. Following the same paper24, we adapted their statistical model for the individual prognosis of the death risk based on the reported hazard ratio for the age levels and for the CCI in multiple European studies as explained in the section 8 of the supplementary material. This gave us a prognosis for the mean, the standard deviation and the median of the crude 1-year death risk and survival probability. The next step in the analysis of the COVID-19 comorbidities concerns the co-occurrence of the diseases. This was explored with the help of the Pearson's phi coefficient as defined and investigated, e.g. in46 and computed for each pair of single or aggregated diseases with an incidence of at least 4. A value of 1 means a "perfect match", 0 means no co-occurrence and negative values up to − 1 means "inverse occurrence". A comorbidity network approach was considered to investigate the relationships between diseases. We computed the dissimilarity between the conditions studied in the previous step as one minus the Pearson's phi. Hence, a value of 1 corresponds to "no association" and values larger than 1 to "inverse association". The hierarchical classification was performed through an agglomerative clustering algorithm with the Ward method since this proved to have the largest agglomeration coefficient. Furthermore, the dendrogram was cut at nodes corresponding to branch length of 1 based both on a visual inspection and on predefined length of the branch. This study has several limitations. We considered only the deaths up to April 20 when the peak of the first pandemic wave was reached in Romania. Moreover, total case design was used because of the missing data in official public communications. Disclaimer This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors. Ethical compliance This research was carried out in accordance with the relevant guidelines and regulations. Informed consent We have obtained informed consent from all participants and/or their legal guardians.
title	Methods
sec	Study design The cohort The Romanian COVID-19 mortality sub-study was conducted between March 22 and April 20, 2020, with the aim to explore the comorbidity patterns in the patients who died because of this new virus. A multiple-case, multiple-center design was chosen to answer this research question with data collected from the official public notifications of MAI concerning the COVID-19 deaths registered by the INSP3,4,8. Initially, 451 COVID-19 deaths were reported by INSP from the beginning of the epidemic until April 20 (Fig. 1). The list was completed afterwards with 69 additional cases due to the delay in the reporting process from some centers of COVID-19 infection. The official report included information about gender, age, county of residence, history of hospitalization and diagnostics, and/ or comorbidities. Only cases with information about gender, age and comorbidities were included in the study due to the main focus of the research42. Eventually, 432 patients who fulfilled these criteria were kept in the Romanian COVID-19 mortality sub-study. Since we aim at an analytic and not statistical generalization of the results regarding the comorbidity patterns, statistical testing was just serving exploratory purposes42. Hence, the complete-case analysis is justified in this case both by purpose and by the state of things.
title	Study design
sec	The cohort The Romanian COVID-19 mortality sub-study was conducted between March 22 and April 20, 2020, with the aim to explore the comorbidity patterns in the patients who died because of this new virus. A multiple-case, multiple-center design was chosen to answer this research question with data collected from the official public notifications of MAI concerning the COVID-19 deaths registered by the INSP3,4,8. Initially, 451 COVID-19 deaths were reported by INSP from the beginning of the epidemic until April 20 (Fig. 1). The list was completed afterwards with 69 additional cases due to the delay in the reporting process from some centers of COVID-19 infection. The official report included information about gender, age, county of residence, history of hospitalization and diagnostics, and/ or comorbidities. Only cases with information about gender, age and comorbidities were included in the study due to the main focus of the research42. Eventually, 432 patients who fulfilled these criteria were kept in the Romanian COVID-19 mortality sub-study. Since we aim at an analytic and not statistical generalization of the results regarding the comorbidity patterns, statistical testing was just serving exploratory purposes42. Hence, the complete-case analysis is justified in this case both by purpose and by the state of things.
title	The cohort
p	The Romanian COVID-19 mortality sub-study was conducted between March 22 and April 20, 2020, with the aim to explore the comorbidity patterns in the patients who died because of this new virus. A multiple-case, multiple-center design was chosen to answer this research question with data collected from the official public notifications of MAI concerning the COVID-19 deaths registered by the INSP3,4,8.
p	Initially, 451 COVID-19 deaths were reported by INSP from the beginning of the epidemic until April 20 (Fig. 1).
p	The list was completed afterwards with 69 additional cases due to the delay in the reporting process from some centers of COVID-19 infection. The official report included information about gender, age, county of residence, history of hospitalization and diagnostics, and/ or comorbidities. Only cases with information about gender, age and comorbidities were included in the study due to the main focus of the research42. Eventually, 432 patients who fulfilled these criteria were kept in the Romanian COVID-19 mortality sub-study. Since we aim at an analytic and not statistical generalization of the results regarding the comorbidity patterns, statistical testing was just serving exploratory purposes42. Hence, the complete-case analysis is justified in this case both by purpose and by the state of things.
sec	Case-case study Additionally, we run a case-case study utilizing data regarding the deaths due to hospital pneumonia between March 22 and April 20 in the years 2016–2018 extracted from the administrative database of the INSP. The study was carried out with the approval of the Local Bioethics Council of the INSP (number 6342/04 MAY 2020) and with the permission of the Local Bioethics Council of the “Carol Davila” University of Medicine, Bucharest Romania (9984/11 MAY 2020).
title	Case-case study
p	Additionally, we run a case-case study utilizing data regarding the deaths due to hospital pneumonia between March 22 and April 20 in the years 2016–2018 extracted from the administrative database of the INSP.
p	The study was carried out with the approval of the Local Bioethics Council of the INSP (number 6342/04 MAY 2020) and with the permission of the Local Bioethics Council of the “Carol Davila” University of Medicine, Bucharest Romania (9984/11 MAY 2020).
sec	Statistical analysis An exploratory data analysis was applied to identify general, spatial-, gender—or age-specific comorbidity patterns in our study. In a first step, the various conditions were sorted according to the standard diagnosis tool ICD-11 (The International Classification of the Diseases) as it was recommended for the mortality and morbidity statistics by the WHO43. The ICD-11 taxonomy includes categories of single diseases or the aggregated groups of diseases. We investigated the comorbidities at both levels of aggregation because the study included a large number of distinct conditions with low prevalence. The statistical analysis was done with the help of the statistical software R version 3.6.2. Continuous variables were expressed as mean, standard deviation (SD) and median. We focused in this paper on comparing the individual crude relative frequencies of the single and aggregated comorbidities between different gender and age groups. Since age is a significant risk factor for mortality, we investigated the frequency of the comorbidities stratified on ten-year intervals starting with the 50+ age groups. Categorical variables were tabulated in frequency (n) and relative frequency (%). A binomial and a chi-squared test for one proportion (only the last was reported) and a chi-squared test for equal proportions were applied when it was appropriate at the 5% level of significance. The post-hoc tests were not corrected for multiple comparisons since the focus of our study was on generating hypotheses for single diseases or group of diseases and not on an experiment-wise approach. The spatial distribution of the COVID-19 mortality corresponded to the country-specific administrative units– counties and was illustrated by single or aggregated disease frequency and relative frequency (%). The association between different comorbidities and the odds ratio between COVID-19 and pneumonia deaths was investigated with the help of the logistic regression, and the results were reported as coefficients β, as point estimates for odds ratio and their 95% confidence intervals. Multimorbidity was computed for each person in the sample since this is a good indicator of the burden of disease44,45 and transformed in a multimorbidity factor with levels 0, 1, 2, 3, 4, or more the 5 for single or aggregated comorbidities. A controversial issue regarding the epidemic of COVID-19 is the burden of disease in the deceased population and the hypothetical survival probability in the case that these people would not have been infected with this virus. A simplified version of the Charlson Comorbidity Index (CCI) was calculated on a scale from 0 to 8, as indicated in24 for an assessment of the individual disease severity. We introduced the CCI factor with two levels (mild for CCI < 3 and severe for CCI ≥ 3), and we tested its differences between gender and age groups by a chi-squared test. Following the same paper24, we adapted their statistical model for the individual prognosis of the death risk based on the reported hazard ratio for the age levels and for the CCI in multiple European studies as explained in the section 8 of the supplementary material. This gave us a prognosis for the mean, the standard deviation and the median of the crude 1-year death risk and survival probability. The next step in the analysis of the COVID-19 comorbidities concerns the co-occurrence of the diseases. This was explored with the help of the Pearson's phi coefficient as defined and investigated, e.g. in46 and computed for each pair of single or aggregated diseases with an incidence of at least 4. A value of 1 means a "perfect match", 0 means no co-occurrence and negative values up to − 1 means "inverse occurrence". A comorbidity network approach was considered to investigate the relationships between diseases. We computed the dissimilarity between the conditions studied in the previous step as one minus the Pearson's phi. Hence, a value of 1 corresponds to "no association" and values larger than 1 to "inverse association". The hierarchical classification was performed through an agglomerative clustering algorithm with the Ward method since this proved to have the largest agglomeration coefficient. Furthermore, the dendrogram was cut at nodes corresponding to branch length of 1 based both on a visual inspection and on predefined length of the branch. This study has several limitations. We considered only the deaths up to April 20 when the peak of the first pandemic wave was reached in Romania. Moreover, total case design was used because of the missing data in official public communications.
title	Statistical analysis
p	An exploratory data analysis was applied to identify general, spatial-, gender—or age-specific comorbidity patterns in our study. In a first step, the various conditions were sorted according to the standard diagnosis tool ICD-11 (The International Classification of the Diseases) as it was recommended for the mortality and morbidity statistics by the WHO43. The ICD-11 taxonomy includes categories of single diseases or the aggregated groups of diseases. We investigated the comorbidities at both levels of aggregation because the study included a large number of distinct conditions with low prevalence. The statistical analysis was done with the help of the statistical software R version 3.6.2.
p	Continuous variables were expressed as mean, standard deviation (SD) and median. We focused in this paper on comparing the individual crude relative frequencies of the single and aggregated comorbidities between different gender and age groups. Since age is a significant risk factor for mortality, we investigated the frequency of the comorbidities stratified on ten-year intervals starting with the 50+ age groups. Categorical variables were tabulated in frequency (n) and relative frequency (%). A binomial and a chi-squared test for one proportion (only the last was reported) and a chi-squared test for equal proportions were applied when it was appropriate at the 5% level of significance.
p	The post-hoc tests were not corrected for multiple comparisons since the focus of our study was on generating hypotheses for single diseases or group of diseases and not on an experiment-wise approach. The spatial distribution of the COVID-19 mortality corresponded to the country-specific administrative units– counties and was illustrated by single or aggregated disease frequency and relative frequency (%).
p	The association between different comorbidities and the odds ratio between COVID-19 and pneumonia deaths was investigated with the help of the logistic regression, and the results were reported as coefficients β, as point estimates for odds ratio and their 95% confidence intervals.
p	Multimorbidity was computed for each person in the sample since this is a good indicator of the burden of disease44,45 and transformed in a multimorbidity factor with levels 0, 1, 2, 3, 4, or more the 5 for single or aggregated comorbidities.
p	A controversial issue regarding the epidemic of COVID-19 is the burden of disease in the deceased population and the hypothetical survival probability in the case that these people would not have been infected with this virus. A simplified version of the Charlson Comorbidity Index (CCI) was calculated on a scale from 0 to 8, as indicated in24 for an assessment of the individual disease severity. We introduced the CCI factor with two levels (mild for CCI < 3 and severe for CCI ≥ 3), and we tested its differences between gender and age groups by a chi-squared test.
p	Following the same paper24, we adapted their statistical model for the individual prognosis of the death risk based on the reported hazard ratio for the age levels and for the CCI in multiple European studies as explained in the section 8 of the supplementary material. This gave us a prognosis for the mean, the standard deviation and the median of the crude 1-year death risk and survival probability.
p	The next step in the analysis of the COVID-19 comorbidities concerns the co-occurrence of the diseases. This was explored with the help of the Pearson's phi coefficient as defined and investigated, e.g. in46 and computed for each pair of single or aggregated diseases with an incidence of at least 4. A value of 1 means a "perfect match", 0 means no co-occurrence and negative values up to − 1 means "inverse occurrence".
p	A comorbidity network approach was considered to investigate the relationships between diseases. We computed the dissimilarity between the conditions studied in the previous step as one minus the Pearson's phi. Hence, a value of 1 corresponds to "no association" and values larger than 1 to "inverse association". The hierarchical classification was performed through an agglomerative clustering algorithm with the Ward method since this proved to have the largest agglomeration coefficient. Furthermore, the dendrogram was cut at nodes corresponding to branch length of 1 based both on a visual inspection and on predefined length of the branch.
p	This study has several limitations. We considered only the deaths up to April 20 when the peak of the first pandemic wave was reached in Romania. Moreover, total case design was used because of the missing data in official public communications.
sec	Disclaimer This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors.
title	Disclaimer
p	This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors.
sec	Ethical compliance This research was carried out in accordance with the relevant guidelines and regulations.
title	Ethical compliance
p	This research was carried out in accordance with the relevant guidelines and regulations.
sec	Informed consent We have obtained informed consent from all participants and/or their legal guardians.
title	Informed consent
p	We have obtained informed consent from all participants and/or their legal guardians.
sec	Supplementary Information Supplementary Information.
title	Supplementary Information
sec	Supplementary Information.
p	Supplementary Information.
caption	Supplementary Information.
p	Supplementary Information.
back	Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. These authors contributed equally: Anca Pantea Stoian, Mihaela Pricop-Jeckstadt, Adrian Pana, Ruxandra Schitea, Antonio Ceriello and Manfredi Rizzo. Supplementary Information The online version contains supplementary material available at 10.1038/s41598-020-78575-w. Author contributions A.P.S., M.P.J., A.P. conceived the research question. A.P.S., M.P.J., A.P., B.V.I., M.R., A.C., R.S. and V.J. designed the study and analysis plan. A.P.S., M.P.J., A.P., R.S., M.G., M.R., A.C., S.P., C.N., S.P., R.H., A.I.S., D.C., L.M. and C.P. prepared the data. M.P.J. did the statistical analysis. A.P.S., M.P.J., A.P., B.T., A.T., C.S., A.A., A.C., M.R., and V.J. drafted the initial and final versions of the manuscript. All authors critically reviewed early and approved final versions of the manuscript. Competing interests APS, MR, AC have given talks, attended conferences, and participated in advisory boards and clinical trials sponsored by various pharmaceutical companies; yet, no financial or professional help was received for the preparation of this manuscript. APS is currently Vice-President, National Diabetes Commission, Ministry of Health, Romania. MR is currently Chief Medical and Scientific Advisor, Diabetes, Novo Nordisk South East Europe, Middle East and Africa (SEEMEA). The rest of the authors do not have a conflict of interests to declare for this article.
footnote	Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
p	Publisher's note
p	Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
footnote	These authors contributed equally: Anca Pantea Stoian, Mihaela Pricop-Jeckstadt, Adrian Pana, Ruxandra Schitea, Antonio Ceriello and Manfredi Rizzo.
p	These authors contributed equally: Anca Pantea Stoian, Mihaela Pricop-Jeckstadt, Adrian Pana, Ruxandra Schitea, Antonio Ceriello and Manfredi Rizzo.
sec	Supplementary Information The online version contains supplementary material available at 10.1038/s41598-020-78575-w.
title	Supplementary Information
p	The online version contains supplementary material available at 10.1038/s41598-020-78575-w.
notes	Author contributions A.P.S., M.P.J., A.P. conceived the research question. A.P.S., M.P.J., A.P., B.V.I., M.R., A.C., R.S. and V.J. designed the study and analysis plan. A.P.S., M.P.J., A.P., R.S., M.G., M.R., A.C., S.P., C.N., S.P., R.H., A.I.S., D.C., L.M. and C.P. prepared the data. M.P.J. did the statistical analysis. A.P.S., M.P.J., A.P., B.T., A.T., C.S., A.A., A.C., M.R., and V.J. drafted the initial and final versions of the manuscript. All authors critically reviewed early and approved final versions of the manuscript.
title	Author contributions
p	A.P.S., M.P.J., A.P. conceived the research question. A.P.S., M.P.J., A.P., B.V.I., M.R., A.C., R.S. and V.J. designed the study and analysis plan. A.P.S., M.P.J., A.P., R.S., M.G., M.R., A.C., S.P., C.N., S.P., R.H., A.I.S., D.C., L.M. and C.P. prepared the data. M.P.J. did the statistical analysis. A.P.S., M.P.J., A.P., B.T., A.T., C.S., A.A., A.C., M.R., and V.J. drafted the initial and final versions of the manuscript. All authors critically reviewed early and approved final versions of the manuscript.
notes	Competing interests APS, MR, AC have given talks, attended conferences, and participated in advisory boards and clinical trials sponsored by various pharmaceutical companies; yet, no financial or professional help was received for the preparation of this manuscript. APS is currently Vice-President, National Diabetes Commission, Ministry of Health, Romania. MR is currently Chief Medical and Scientific Advisor, Diabetes, Novo Nordisk South East Europe, Middle East and Africa (SEEMEA). The rest of the authors do not have a conflict of interests to declare for this article.
title	Competing interests
p	APS, MR, AC have given talks, attended conferences, and participated in advisory boards and clinical trials sponsored by various pharmaceutical companies; yet, no financial or professional help was received for the preparation of this manuscript. APS is currently Vice-President, National Diabetes Commission, Ministry of Health, Romania. MR is currently Chief Medical and Scientific Advisor, Diabetes, Novo Nordisk South East Europe, Middle East and Africa (SEEMEA). The rest of the authors do not have a conflict of interests to declare for this article.

Annnotations TAB TSV DIC JSON TextAE

last updated at 2021-10-18 14:17:05 UTC

Denotations: 0
Blocks: 0
Relations: 0

PMC:7730445 / 13888-13892 JSONTXT

Document structure show

Annnotations TAB TSV DIC JSON TextAE

PMC:7730445 / 13888-13892 JSON TXT