Comparing COVID-19 and Influenza Presentation and Trajectory

Background: COVID-19 is a newly recognized illness with a predominantly respiratory presentation. It is important to characterize the differences in disease presentation and trajectory between COVID-19 patients and other patients with common respiratory illnesses. These differences can enhance knowledge of pathogenesis and help in guiding treatment. Methods: Data from electronic medical records were obtained from individuals admitted with respiratory illnesses to Rambam Health Care Campus, Haifa, Israel, between October 1st, 2014 and October 1st, 2020. Four groups of patients were defined: COVID-19 (693), influenza (1,612), severe acute respiratory infection (SARI) (2,292), and Others (4,054). The variable analyzed include demographics (7), vital signs (8), lab tests (38), and comorbidities (15) from a total of 8,651 hospitalized adult patients. Statistical analysis was performed on biomarkers measured at admission and for their disease trajectory in the first 48 h of hospitalization, and on comorobidity prevalence. Results: COVID-19 patients were overall younger in age and had higher body mass index, compared to influenza and SARI. Comorbidity burden was lower in the COVID-19 group compared to influenza and SARI. Severely- and moderately-ill COVID-19 patients older than 65 years of age suffered higher rate of in-hospital mortality compared to hospitalized influenza patients. At admission, white blood cells and neutrophils were lower among COVID-19 patients compared to influenza and SARI patients, while pulse rate and lymphoctye percentage were higher. Trajectories of variables during the first 2 days of hospitalization revealed that white blood count, neutrophils percentage and glucose in blood increased among COVID-19 patients, while decreasing among other patients. Conclusions: The intrinsic virulence of COVID-19 appeared higher than influenza. In addition, several critical functions, such as immune response, coagulation, heart and respiratory function, and metabolism were uniquely affected by COVID-19.


INTRODUCTION
Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), is the virus underlying COVID-19, a newly recognized illness that initially spread throughout Wuhan (Hubei province), and from there, to other provinces in China and then across the globe. As of October 15th 2020, over 39,000,000 infections and over 1,100,000 casualties have been linked to SARS-CoV-2. The clinical spectrum of SARS-CoV-2-associated pneumonia ranges from mild to life-threatening (1,2). Several studies have described general epidemiological findings, clinical presentation and clinical outcomes of SARS-CoV-2 pneumonia, and identified mortality risk factors (3)(4)(5)(6)(7)(8)(9). The need for detailed information on the clinical characteristics of hospitalized patients with COVID- 19 and their clinical course is essential to achieve a thorough understanding of the disease development and progression. Moreover, the exact differences between the clinical presentation and illness trajectory of COVID-19 vs. other respiratory (viral) infections remain illusive. Investigating the clinical features of influenza like illness (ILI) is of paramount importance to identify COVID-19 specificities. This work questioned the virulence of COVID-19 as compared to seasonal influenza and SARI. Recognizing the characteristics discriminating COVID-19 from influenza, will be critical to support the management of the current pandemic.
There are limited works that have pursued the pathophysiological differences between ILIs and COVID-19, and have focused solely on the H1N1 influenza strain (10)(11)(12)(13). They investigated symptoms, comorbidities, laboratory examinations, treatments and scans in relatively small cohorts of patients and consequently had limited statistical power. Hence, the relative virulence of COVID-19 vs. influenza has been and remains under debate. Discriminating biomarkers at clinical presentation and characteristics of the trajectory of COVID-19 vs. influenza are poorly characterized.
To address this knowledge gap, this retrospective analysis reviewed the electronic medical records (EMR) from the Rambam Health Care Campus, located in Haifa, Israel. Demographics, comorbidities, vital signs and laboratory tests were analyzed to identify features that can potentially discriminate between COVID-19, influenza and SARI at the time of admission to the hospital, as well as their trajectory during the first 48 h after admission.

Data Source
A unique cloud-based database, named COV19, was created based on the model of MIMIC III (14). The database contained detailed de-identified clinical information. Specific views (tables) that contained multiple variables related to a given type of medical data, were created. Ethical approval for this research was provided by the local institutional review board (IRB; #0141-20). The de-identified datasets were uploaded to a Microsoft Azure cloud server, which also offers data analysis, visualization, and querying tools. The description of the tables and variables included in COV19 is available on the online resource site (https://cov19-resource.com/). Future access to the cloud can be given to interested researchers, subject to hospital IRB approval.

Study Population
This single center retrospective observational cohort study uses EMR data from Rambam Health Care Campus, a 1,000-bed tertiary academic hospital in Northern Israel, during which the pandemic opened five dedicated COVID-19 departments. The hospital EMR database was queried for hospitalized adult (age 18 and above) cases admitted for COVID-19, influenza or SARI, or tested for COVID-19 during hospitalization (with either a positive or negative result), from October 1st, 2014 until October 1st, 2020.

Disease Groups
The disease groups were defined as follows: COVID-19: At least one positive reverse transcription polymerase chain reaction (RT-PCR) test for SARS-CoV-2 in nasopharyngeal swab. Most COVID-19 cases were positive within a week before admission to the hospital or at admission, and very few of them were diagnosed a few days after admission. COVID-19 were also tested for Influenza. Influenza: tested positive for influenza A or B virus by RT-PCR test and tested negative for COVID-19, or tested positive for influenza A or B virus by RT-PCR test and admission date prior to COVID-19 emergence in Israel, on February 23, 2020. SARI only: Physician report in the EMR that matches World Health Organization (WHO) SARI case definition (15), i.e., an acute respiratory infection with history of fever or measured fever of ≥38 • C, and cough with onset within the last 10 days, hospitalization, and no positive test for COVID-19 or influenza (either negative or not tested). Others: Tested negative for COVID-19 and not classified into influenza or SARI groups.
Forty percent of the SARI cases were tested for Influenza, when indicated by the physician for differential diagnosis or for the purpose of surveillance during the influenza season. As the cohort is defined by inclusion criteria as all patients tested for COVID-19, or either testing positive for influenza or having a SARI diagnosis, the "others" group is a control group inevitably generated by those with negative result for COVID-19, that are also negative for Influenza and SARI. Thus, this group contains zero cases prior to February 23, 2020, the date of COVID-19 emergence in Israel. After this date, which approaches the end of the Influenza season, only 5% of the patients in "others" group were tested for Influenza. Tests for COVID-19 after its emergence were performed for 80% of the SARI patients, and by inclusion criteria, for all patients in "others."

COVID-19 Severity
Following existing guidelines within the context of COVID-19 (16)(17)(18), patients were defined as moderately ill if they were diagnosed for COVID-19 pneumonia clinically or by X-ray. Patients were defined as severely ill if either their breaths number per minute was larger than 30, their unsupported oxygen saturation was of 93% or lower, their PF ratio was below 300, or were critically ill. Critically ill patients were defined as those who either went through mechanical ventilation support (invasive or non-invasive), were hospitalized in an intensive care unit, or were administered vasopressor medications (noradrenaline and vasopressin) or inotropic medications (dopamine, dobutamine, milrinone, and adrenaline).

Collected Data
Demographic information, such as age, sex, ethnic group (Jewish or Arabs, which included Druze and Muslims), weight, body mass index (BMI), length of hospitalization and mortality rates were collected. In addition, comorbidities, vital signs including fever, respiratory variables (breaths count per minutes, oxygen saturation) blood pressure and tests results, such as metabolic profiles, complete blood count, and coagulation tests were collected. Comorbidities were defined by ICD-9 codes as detailed in Supplementary Table 1. A thorough analysis (see "Statistical Analysis" section) of comorbidities, demographics and mortality rate was performed to characterize predispositions and the severity and case fatality of each disease. Moreover, lab tests and vital signs were screened at admission and for the first 2 days of hospitalization in order to identify distinct signature or putative biomarkers of COVID-19, influenza, and SARI.

Statistical Analysis
Demographic variables, comorbidity rate, vital signs, and lab tests were compared between disease groups at admission using the Chi-squared test or Fisher's exact test for categorical variables, and analysis of variance or Kruskal-Wallis test for continuous variables. The p-values across all tests were corrected to control the false discovery rate (FDR) criterion (19). Mortality rates were compared before and after excluding patients at mild severity level. Medians and inter-quartile range (IQR) were used to describe the continuous variables. In addition, standardized scores were calculated for scale unification across variables, using the median for centralization and the median absolute deviation (MAD) for rescaling, thereby allowing their representation in a comparative heatmap. Adjustments for confounders were performed using generalized linear models. Age-adjusted COVID-19 odds ratios for each comorbidity were calculated using multivariate logistic regression, excluding mild cases, to eliminate severity bias.
The trajectory over time was compared for each numeric measure, using a non-parametric repeated measure model for a factorial design (20). In the first stage, for each measure, we used the interaction effect in the model to test for difference in time trend between the disease groups. For variables that showed a significant effect, we conducted post hoc pairwise tests between the groups. The p-values over all tests across the two stages were corrected using a hierarchical FDR controlling procedure (21). Three time intervals were defined to follow trends during hospitalization: 0-6, 6-24, and 24-48 h from admission. Effect size was defined by the difference in slopes across time between each pair of compared disease groups. The larger slope difference among the slopes obtained for (0-6, 6-24) and (6)(7)(8)(9)(10)(11)(12)(13)(14)(15)(16)(17)(18)(19)(20)(21)(22)(23)(24) time gaps was selected as the effect size for each pairwise comparison. The R software (22) was used for statistical analysis, including the R package nparLD (23) for applying the non-parametric model. A 0.05 threshold was used to determine significance.

Database
A total of 9,670 hospitalizations at Rambam Health Care Campus met the initial inclusion criteria between October 1st, 2014 and October 1st, 2020 (Figure 1). Excluded were 33 noncritical COVID-19 cases admitted before April 1, 2020, a period during which all positive cases were systematically hospitalized, including very mild and asymptomatic cases, and 23 COVID-19 cases admitted for reasons unrelated to COVID-19 (e.g., women at labor, traffic accident injury). In order to remove bias and correlations due to repeated per-patient hospitalizations 963 cases, mainly non-COVID-19 hospitalization, with fewer than 30 days between consecutive hospitalization, were excluded. The remaining 8,651 cases were classified into the four disease groups. This included 693 COVID-19 patients, of whom 127 (18.3%) were classified as critical. A total of 68 variables were evaluated: demographics (7), vital signs (8), lab tests (38), and comorbidities (15).

Demographics
As shown in Table 1, COVID-19 patients (median age 59.8 ± 29.7 IQR years) were younger compared to influenza patients (median age 70.5 ± 21.8 IQR years) and SARI patients (median age 70.5 ± 22.9 IQR years). Patients between 18 and 44 years of age made up 25% of the COVID-19, compared to 13-14% among influenza and SARI patients. In contrast 24% of the COVID-19 patients were 75 years or older, compared to 40% among influenza and SARI patients. The sex distribution among COVID-19 patients reflected a small preference toward males (52%), compared to a more significant preference toward males among SARI patients (58.6%), and no preference among influenza patients (50%). The proportion of Arabs among the COVID-19 patients (40%) was larger compared to influenza and SARI patients (25 and 23.5%, respectively), independently of disease severity. Duration of hospitalization was shorter for COVID-19 patients (median 4 ± 7 IQR days) compared to influenza and SARI patients (respectively median 5 ± 5 IQR and median 5 ± 6 IQR days). However, for moderate to severe COVID-19 patients, a longer duration of hospitalization was observed (median 6 ± 7.75 IQR days) as compared to influenza and SARI.
BMI was higher for COVID-19 patients (median 28.7 ± IQR 7.2 kg/m 2 ) compared to influenza (median 27.8 ± IQR 7.3 kg/m 2 ) and SARI (median 26.6 ± 7.4 IQR kg/m 2 ), with 35% of the COVID-19 patients being obese (BMI > 30 kg/m 2 ), compared to 30% of the influenza patients and 24% of the SARI patients. Moreover, 48% of moderate to severe COVID-19 patients under 65 years of age with no comorbidities were obese, compared to 26% of influenza patients and 14% of SARI patients. After adjustment for age, BMI among moderate to severe COVID-19 patients was found larger by 1.7 kg/m 2 compared to influenza patients, and by 2.9 kg/m 2 compared to SARI patients. Overall in-hospital mortality was similar between COVID-19 (8%), influenza (8.1%), and other (8.6%) group, and markedly higher for the SARI group (13.5%). To limit the selection bias introduced by the current pandemic situation that would result in a more permissive hospitalization policy for the COVID-19 patients, we excluded mild COVID-19 cases. This analysis resulted in an in-hospital death rate of 12.6% for moderate to severe COVID-19 patients. When focusing on age groups 65-74 and 75+ years, mortality rates of COVID-19 patients were 15.5 and 21.4%, respectively. Moreover, in-hospital mortality stratified by BMI showed that a lower proportion of patient with BMI over 30 kg/m 2 died compared to the group with BMI between 20 and 25 kg/m 2 (5.2 and 9.2%, respectively).

Comorbidities
As shown in Table 2 and depicted in Figure 2, the incidence rates for almost all analyzed comorbidities were lower among moderate to severe COVID-19 patients compared to other patients. A total of 9.5% of COVID-19 patients had cancer, compared to 17.1% of influenza patients and 23.4% of SARI patients. Accordingly, the age-adjusted odd ratio (AOR) for COVID-19 were 0.4 (95%CI: 0.27; 0.59). Similar results were obtained when separating to solid-type cancer and hematologic cancer. In total, 3.5% of moderate to severe COVID-19 patients had chronic obstructive pulmonary disease (COPD), compared to 10-11% of influenza and SARI patients and the AOR for COVID-19 was 0.29 (95%CI: 0. 16 . However, an interaction effect between dementia and arrival from nursing home was observed (p-value = 0.0047). Stratified AOR led to nearly significant levels due to smaller sample size. COVID-19 patients hospitalized from nursing homes depicted an AOR of 1.82 (p-value = 0.058), while the AOR for COVID-19 patients not coming from nursing homes was 0.54 (p-value = 0.079).
Lastly, diabetes was an important risk factor for COVID-19 severity, with a prevalence of nearly 30% of COVID-19 cases. Influenza and SARI depicted a similar prevalence suggesting diabetes as a common comorbidity of ILIs.

Vital Signs at Admission
The median standardized score (Figure 3A) of the heart rate (pulse) among COVID-19 patients (median 88 ± 23 IQR bpm) was substantially lower compared to influenza (median 93 ± 28 IQR bpm) and SARI (median 97 ± 29 IQR bpm) patients. Systolic blood pressure was lower among COVID-19 patients (median 128 ± 27 mmHg) compared to influenza patients (median 135 ± 37 mmHg), while diastolic blood pressure was higher among COVID-19 patients (median 77 ± 14 mmHg) compared to influenza (median 73 ± 19 mmHg) and SARI patients (median 72 ± 17 mmHg). The standardized scores (in Figure 3) showed that among the vital signs, blood pressure parameters (systolic, diastolic, and average) were regrouped as a single cluster. Similarly, saturation parameters (oxygen and room) formed a single group. A third group included pulse, fever, and breath number. However, neither respiratory measures nor saturation or temperature showed clinically relevant median differences between ILIs.

Lab Tests at Admission
Laboratory examinations are presented in Table 4 and in Figure 3B. Lymphocytes percentage among COVID-19 patients (median 19.9 ± 15.5 IQR %) was substantially higher compared to influenza patients (median 11.5 ± 10.9 IQR %) and SARI patients (median 10.4 ± 11 IQR %). Similarly, albumin levels were higher among COVID-19 patients (median 3.8 ± 0.6 IQR g/dL) compared to influenza patients (median 3.4 ± 0.8 IQR g/dL) and SARI patients (median 3.1 ± 0.8 IQR g/dL). The  (abs), White blood count (WBC), brain natriuretic peptide (BNP) and neutrophils percentage formed a single cluster that indicated substantially lower levels in COVID-19 patients as compared to influenza and SARI. A similar effect with reduced levels in COVID-19 patients as compared to influenza and SARI was shown for the cluster containing met-hemoglobin and carboxyhemoglobin, the cluster containing procalcitonin, bicarbonate, alanine aminotransterrase (ALT), potassium and phosphorous, and the cluster containing creatinine, lactate, blood urea nitrogen (BUN), and glucose.

Disease Trajectories
Trends over time of parameters that showed differences between patient groups, are presented in Table 5. The results for of all pairwise group comparisons are provided as Supplementary Table 2. For each variable and compared pair, the table provides the time gap for which the stronger difference in rank means slopes was found, either (0-6, 6-24) or (6)(7)(8)(9)(10)(11)(12)(13)(14)(15)(16)(17)(18)(19)(20)(21)(22)(23)(24), and the difference itself. Figure 4 represents the trends for each of the variables showing a significant effect on diseases trajectories. As shown by the trends in Figure 4, both systolic and diastolic blood pressures, as well as pulse, decreased more slowly during the 6-48 h after admission among COVID-19 patients compared to influenza and SARI patients. The difference in slopes with respect to its significance level (p-value) are represented in the volcano plot in Figure 5.

DISCUSSION
The relative virulence of COVID-19 vs. influenza has been and is still debated. Therefore, our focus was to assess COVID-19 mortality and virulence compared to other ILIs.
In-hospital death rates were investigated. While COVID-19 depicted a lower overall mortality, stratification by age showed an increased fraction of casualty among elderly persons consistent with other reports (24). Age has been associated with COVID-19 severity and mortality in plethora of studies (1,2,(4)(5)(6)(7)9). The presented comorbidity analysis exposed striking differences between disease groups. COVID-19 patients depicted substantially fewer "overall comorbidities" than influenza or SARI patients (66.3 vs. 81.6 and 81.9%). The difference remained after age adjustment and excluding mild severity, suggesting that COVID-19 patients required hospitalization more frequently even if they were healthier before their infection. Hypertension, cardiovascular disease and diabetes depicted a large prevalence in COVID-19, similarly to previous findings (3,6,25), and represented common risk   (26). A recent study suggested that the differential immune cell profiles of cancer patients treated with immunomodulatory agents, may impact the host response to the SARS-COV2 diminish disease severity (27). The lower prevalence can also be explained by higher self COVID-19 risk perception among cancer patients that leads to higher adherence to isolation methods. Taken together, the present investigation demonstrates that COVID-19 hospitalized patients present a significantly lower amount of comorbidities and are younger. COVID-19 patients hospitalized with moderate to severe disease had a significantly higher mortality rate than hospitalized influenza patients. Overall, this suggests that the intrinsic virulence of COVID-19 is higher than influenza. Our analysis of vital signs revealed that pulse, blood pressure and temperature (fever) were significantly different between ILIs at admission and showed different patterns in the 2 days postadmission despite the small effect size on fever. An independent study on acute respiratory syndrome (10) comparing Influenza and COVID-19 patients (severe) showed that partial pressure of oxygen in arterial blood (PaO 2 ) was remarkably low for COVID-19 patients compared to influenza. Hence, oxygen saturation, heartbeat, temperature, and others continuous measures (i.e., blood glucose, activity) may serve as informative parameters to be continuously monitored in COVID-19 patients. For example, such measurements could be used to anticipate hospitalization, at the point of care, or to predict potential readmission (28)(29)(30)(31)(32)(33)(34)(35). Monitoring vital signs at the hospital is critical for following patient disease trajectory (36) and for guiding physicians on clinical interventions, such as initiation of ventilation assistance or pharmaceutical treatment.
Obesity has been reported as a major risk factor for COVID-19 and type 2 diabetes. Furthermore, impaired metabolic health (i.e., dyslipidemia and insulin resistance) is linked with a higher risk of COVID-19-associated pneumonia (37)(38)(39). The presented analysis showed that BMI was indeed significantly higher in COVID-19 patients with respect to SARI and Influenza, while diabetes prevalence was high but similar to other ILIs (∼30%). The lower BMI group depicted a higher mortality as compared to obese patients. However, due to the high proportion of missing data on BMI among the patients that died and the relatively small sample size, the observed differences warrant further investigations. This obesity paradox has been reported elsewhere (40) and seems to be common between ILIs (41). Several  hypothesis have been proposed to explain why obese patients are more affected by COVID-19 without an increased mortality (42). For instance, Adipocytes contains the ACE2 receptor that enable the entrance of the COVID-19 virus, which turns adipose tissue into a potential target and viral reservoir (43). On the other side, Influenza infection relies on a different mechanism, through hemagglutinin onto sialic acid sugars on the surfaces of epithelial cells (44) which could explain a higher severity of COVID-19 disease in overweight patients. Conversely, in obese patients, an increased circulating levels of adipokines and inflammatory cytokines, such as TNFα, IL-6, or C-Reactive Protein (CRP) was observed (45). This chronic low-grade inflammation may impair the adaptive immune responses to viral infections (46) and consequently reduce the probability of a lethal cytokine storm. However, it remains controversial for the COVID-19 disease and the underlying mechanisms remains unclear (47). At the laboratory examination level, several metabolic intermediates and enzymes were significantly different between COVID-19, SARI, and influenza patients, particularly, lactate, LDH and glucose. Intriguingly, glucose levels in COVID-19 patients were lower at admission but took an opposite trajectory in the first 2 days of admission, compared to SARI and influenza. Impaired glucose homeostasis is associated with poor COVID-19 prognosis and has been hypothesized as an underlying trigger of the cytokine storm in COVID-19 patients (48)(49)(50)(51). In severe cases of COVID-19, a hyperinflammatory response (cytokine storm) is correlated with poor outcome (52)(53)(54). Several studies reported an association between the neutrophils-to-lymphocyte ratio (NLR) and the severity of COVID-19 (55)(56)(57). Accordingly, our investigation of lab examination revealed that COVID-19 have a prominent effect on blood cell-types counts and proportions, such as lymphocytes, neutrophils, WBCs at admission and along the first 2 days of hospitalization. The present analysis reveal that  . The most highly significant result was for average blood pressure, as seen at the left window, which refers to all influenza-COVID-19 comparisons, and at the right window, which refers to all SARI-COVID-19 comparisons. Trends for WBC, glucose, and diastolic blood pressure for COVID-19 patients were found to be different from both influenza and SARI.
the trajectories of glucose and immune cells are affected by the clinical management of COVID-19 patients, probably through corticosteroid treatments (58).
Ferritin is a key player of immune dysregulation through its immune-suppressive and pro-inflammatory effects, supporting the possibility of cytokine storm (59,60). We observed that ferritin levels were more elevated among COVID-19 patients as compared to SARI or influenza patients at admission. However, here the ferritin test was performed only for a small number (∼1.5%) of severe COVID-19 patients.
A non-neglectable risk of thrombosis or disseminated intravascular clot formation has been described for COVID-19 patients (53,61). Histopathological studies in post-mortem lung tissue revealed pulmonary microthrombi in 57% of COVID-19 and 58% of SARS as compared to 24% of H1N1 influenza patients (62). However, the underlying mechanisms remain to be fully characterized. It has been hypothesized that the dysregulated immune responses orchestrated by inflammatory cytokines is involved. Interestingly, our analysis of laboratory tests at admission confirmed higher level of PT and PTT in COVID-19 patients while platelet levels did not show extreme values for COVID patients with respect to influenza or SARI patients. In addition, Ddimers depicted a much lower value in influenza compared to COVID-19 or SARI patients. However, this observation might be biased due to the low number of tests made for influenza patients (∼5%). As the risk due to thrombosis has been described early in the pandemic (63), D-dimer quantification is currently performed routinely for COVID-19 patients.
The diagnosis of dementia has been reported as an important risk factor for mortality in COVID-19 patients (64)(65)(66). The prevalence of demented patient found in the COVID-19 cohort (∼10%) was lower than previous estimates, which vary from 13 to 42% (67). According to our analysis, dementia was found to be more prevalent among COVID-19 patients as compared to patients of the other cohorts. In order to determine if this finding can be explained by the higher proportion of COVID-19 patients arriving from nursing homes, the variable of nursing home residence was added to the model. A significant interaction between the nursing home residence variable and dementia in COVID-19 patients was observed. The results can be potentially explained by a higher exposure and infection risk of demented nursing home residents caused by a lower ability to adhere to the needed isolation behaviors (68). COVID-19 exposure in nursing home has been reported by others and is a common issue in several countries (24,69). More than 20% of all reported COVID-19-associated deaths occurring in nursing homes in countries, such as Canada, Sweden, and the UK. For instance, nursing home population depicted a 1.70fold higher infection attack rate than the general population in France (24,69). Thus, COVID-19 exposure in nursing home has been reported by others and is a common issue in several countries.
The above findings point to vital signs and lab results, as well as to comorbidities, that can be helpful in discriminating between COVID-19 and other ILI's at presentation. We plan to pursue a study to develop data-driven tools to efficiently identify COVID-19 cases and predict their disease trajectory at the point of care.
Haifa and the Northern district of Israel are a multicultural region that encompasses several ethnic groups, who have distinct cultural agendas and socio-economic characteristics. The data collected from the largest regional hospital provided were of value for assessment of differential COVID-19 manifestations and risk factors in various ethnic groups. A striking enrichment of all COVID-19 cases in Arab communities accounting for about 40% of COVID-19 cases was observed, as opposed to roughly 25% for the influenza and SARI cases, similar to the Arab proportion in the Haifa area. Higher exposure risk and infection rates of COVID-19 at Arab settlements may account for this effect as no differences in terms of severity between Arabs and Jews were observed. Further stratification of the Jewish population could not be performed for a more detailed analysis of the orthodox or secular fraction. However, it was previously reported that a similar effect is likely in the orthodox Jewish communities (70)(71)(72).
Other studies that compared COVID-19 to Influenza did not refer to associations between variables, which may be attributed to systematic activity, and to their clinical trajectory, and used limited cohorts of patients or focused only on acute respiratory distress syndrome, for instance (10,11). Furthermore, they did not account for multiple testing of variables, which may lead to increased overall type I error (19). In this study, a comprehensive approach was applied to accommodate simultaneous testing of multiple biomarkers. Heterogeneity in distribution between variables was addressed by unifying their scale through robust standardization, thereby allowing comparison of their effects and detection of common patterns through clustering. Nonparametric repeated measures model was performed to allow for time-course analysis under various forms of non-symmetric, long-tailed distributions.

Study Limitations
The present study has few important limitations. Firstly, this study is based on a single medical center. Secondly this analysis did not include important information, such as symptoms, images (scans, X-rays), waveforms (e.g., electrocardiogram, oxygen measurements in blood, respiration, ventilation), omics data (e.g., genomics, epi-genomics, transcriptomics, lipidomics, and metabolomics) and pharmacological interventions (apart from the severe cases definition). We decided not to include oxygen measurement because of the uncertainly as per whether the measurement was taken on an arterial or venous line. Thirdly, potential selection bias, information bias, or nondifferential misclassification associated with the use of EMR data are existing in the presented dataset. Variables were collected from a unique EMR sources and using a single query tools. These EMR were curated similarly for each groups which allowed to observe clinically relevant significant differences. Differences in the prevalence of personal characteristics between the diseases among hospitalized patients can be related to differences in exposure risk, infection susceptibility, disease severity and admission to hospital policies and biases. To address the effect of selection bias related to differences in hospitalization admission policy of COVID-19 patients compared to the other diseases, asymptomatic and mild cases were excluded in part of the analysis. Asymptomatic and mild cases are usually not hospitalized in influenza and SARI groups unless combined with other risk factors (Berkson's bias) (73).
Furthermore, some of the lab tests are not part of the routine work, hence selection bias by indication is likely. Not all patients within the SARI group were tested for Influenza, COVID-19 or other viruses, thus the SARI group potentially contains undetected Influenza cases, and less likely, COVID-19 cases, as the vast majority patients with Influenza or SARI in our cohort were diagnosed prior to COVID-19 emergence. Nevertheless, the SARI group is heterogeneous in diagnoses, while the Influenza group, which is the focus in this paper for comparison to the COVID-19 group, remains homogeneous and clearly defined. Potential misclassification of influenza within the "others" group is expected to be negligible. This group contains cases with a negative COVID-19 test result, at a period starting at the emergence of COVID-19 in Israel, which was the end of the 2020 influenza season, and ending before the beginning of the next influenza season (74).

Conclusions
The intrinsic virulence of COVID-19 appears higher than influenza. Several critical functions, such as immune response, coagulation, heart and respiratory function and metabolism were markedly impacted by the COVID-19 disease, despite some similarities observed with influenza and SARI. Moreover, COVID-19 seems to differently affect specific segments of the population, potentially due to increased exposure in localized communities or in nursing homes.

DATA AVAILABILITY STATEMENT
The datasets presented in this article are not readily available because of privacy restrictions. Requests to access the datasets should be directed to Dr. Ronit Almog, r_almog@rambam.health.gov.il.

ETHICS STATEMENT
The studies involving human participants were reviewed and approved by the Rambam Health Care Campus Institutional Review Board. Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.

AUTHOR CONTRIBUTIONS
AR: statistical approach and methodology, data analysis, visualization and interpretation, data curation, manuscript writing, and revision. JS: data analysis, visualization, and interpretation, manuscript writing, and revision. RA: study conception, medical and epidemiological interpretation, data curation, and manuscript revision. SL, TB, and AJ: software/infrastructure development. DE: study conception, medical interpretation, and manuscript revision. JB: study conception, supervision, funding, and manuscript revision.