Predicting COVID-19 severity using major risk factors and received vaccines =========================================================================== * Ariel Israel * Alejandro A. Schäffer * Eugene Merzon * Ilan Green * Eli Magen * Avivit Golan-Cohen * Shlomo Vinker * Eytan Ruppin ## Abstract **Background** Vaccines are highly effective in preventing severe disease and death from COVID-19, and new medications that can reduce severity of disease have been approved. However, many countries are facing limited supply of vaccine doses and medications. A model estimating the probabilities for hospitalization and mortality according to individual risk factors and vaccine doses received could help prioritize vaccination and yet scarce medications to maximize lives saved and reduce the burden on hospitalization facilities. **Methods** Electronic health records from 101,039 individuals infected with SARS-CoV-2, since the beginning of the pandemic and until November 30, 2021 were extracted from a national healthcare organization in Israel. Logistic regression models were built to estimate the risk for subsequent hospitalization and death based on the number of BNT162b2 mRNA vaccine doses received and few major risk factors (age, sex, body mass index, hemoglobin A1C, kidney function, and presence of hypertension, pulmonary disease and malignancy). **Results** The models built predict the outcome of newly infected individuals with remarkable accuracy: area under the curve was 0.889 for predicting hospitalization, and 0.967 for predicting mortality. Even when a breakthrough infection occurs, having received three vaccination doses significantly reduces the risk of hospitalization by 66% (OR=0.339) and of death by 78% (OR=0.223). **Conclusions** The models enable rapid identification of individuals at high risk for hospitalization and death when infected. These patients can be prioritized to receive booster vaccination and the yet scarce medications. A calculator based on these models is made publicly available on [http://covidest.web.app](http://covidest.web.app) Keywords * COVID-19 * disease severity * calculator * diabetes * obesity * kidney disease ## Introduction Since the start of the COVID-19 pandemic, over 250 million individuals have been infected and over 5 million individuals have died ([https://coronavirus.jhu.edu/map.html](https://coronavirus.jhu.edu/map.html)). Since late 2020, multiple vaccines have been developed1,2, and more recently, various treatments including monoclonal antibodies, and anti-viral molecules have been FDA approved3. Unfortunately, there are still not enough vaccine doses and medications available worldwide. Thus, it is an urgent unmet need to identify objectively which patients would most benefit from booster vaccination doses or early access to the yet scarce medicines4–6. To address these questions, we conducted a retrospective study among members of Leumit Health Services (LHS), one of the four main health maintenance organizations (HMOs) in Israel, which has over 700,000 members. Our study starts at the beginning of the pandemic and hence covers periods before any vaccines were available and after. Israel was one of the first countries to implement a large-scale vaccination plan (using Pfizer/BioNTech BNT162b2 vaccine)7 and to deploy a third vaccine booster dose. Nevertheless, there exists a large variation in vaccination uptake8, which provides an opportunity to assess the beneficial effects of vaccination in conjunction with other patient risk factors. In addition to vaccinations, we considered other traits associated with COVID-19 severity, including type II diabetes 9–14, kidney disease 11,14–16, chronic obstructive pulmonary disease (COPD) 16–20, obesity10–12,21,22, hypertension 23–25, malignancy26 and advanced age 11,27–29. We constructed predictive models that estimate the risks that patients newly infected with SARS-CoV-2 (as reflected by positive PCR tests) would require hospitalization during the disease course and would die from COVID-19. Predictions are based on patient’s age, sex, the clinical factors mentioned above and vaccination status at the time of infection (0,1,2, or 3 doses). Importantly, all these factors are typically part of a patient’s medical records and measured in routine laboratory testing. Models were built based on the outcomes of LHS COVID-19 patients during the study period. To keep the models simple and interpretable, and to allow deployment of the models in any health provider information system, we used multivariable logistic regression models based on the most essential risk factors. Regression coefficients and odds ratio (OR) for each factor are provided, and these can be calculated using a simple formula to obtain risk estimates for any individual. The resulting models additionally allow to estimate the decrease in hospitalization and mortality risk that can be achieved for any individual by providing additional vaccine doses. These estimates can help identify individuals and segments of population for which additional vaccination could help most. Notably, most treatments available should be given early in the course of the disease to those at high risk, to be most effective. A web-based calculator is provided, and the approach to run or adapt the models is fully described. ## Methods ### Study subjects and study design This is a population study performed in Leumit Health Services (LHS), a nationwide healthcare provider in Israel, which provides services to around 700,000 members throughout the country. LHS uses centrally managed electronic health records (EHRs), continuously updated regarding subjects’ demographics, medical diagnoses, medical encounters, hospitalizations and laboratory tests. All LHS members have similar health insurance and similar access to healthcare services. The study is based on members of LHS, of age ≥ 5 (eligible for vaccination), who had at least one positive PCR test for SARS-CoV-2 between April 2020 and November 30, 2021. Patients’ data were extracted from LHS central data warehouse on January 3, 2022. For each COVID-19 episode, the date of the first positive PCR test was taken as the index date. Number of vaccine doses received were calculated at the index date. Diagnosis and laboratory data were queried as known 15 days before the index date. The following factors were included in the analysis: sex, age, the Body Mass Index (BMI) as a categorical variable (<18.5; 18.5-25; 25-30; 30-35; ≥35), hemoglobin A1C range (<6.5; 6.5-8; 8-10; ≥10), last glomerular filtration rate (GFR) as an estimate of kidney function as a categorical variable (categories: ≥90; 60-89; 45-59; 30-44; <30). Presence of comorbidity conditions was assessed by presence of an active chronic diagnosis at this date. Chronic diagnoses, coded according to the International Classification of Diseases 9th revision (ICD-9), are regularly added, updated or ended, by the treating physician, at each encounter. The validity of chronic diagnoses in the registry has been previously examined and confirmed as high30–32. To keep the models as simple as possible, we deliberately limited the models to the conditions that we identified as having the most significant effect on disease severity: hypertension, chronic obstructive pulmonary disease (COPD) and malignancy (solid or hematologic). Individuals who had a pregnancy diagnosis up to 210 days before the PCR test were excluded, as hospitalization would often pertain to pregnancy surveillance or delivery, and not reflect disease severity. ### Ethics statement The study protocol was approved by the statutory clinical research committee of Leumit Health Services and the Shamir Medical Center Institutional Review Board (129-2-LEU). Informed consent was waived because this is a large-scale retrospective study and all data were deidentified. ### SARS-CoV-2 testing by real-time RT-PCR Nasopharyngeal swabs were taken and examined for SARS-CoV-2 by real-time RT-PCR performed with internal positive and negative controls, according to World Health Organization guidelines, using TaqPath™ COVID-19 Combo Kit (Thermo Fisher Scientific) and COBAS SARS-Cov-2 6800/8800 (Roche Diagnostics) assays. ### Statistical analyses Standard descriptive statistics were used to present the demographic characteristics of individuals included in the study cohort. Statistical analyses were done with R version 4.0.4 (R Foundation for Statistical Computing). Multivariable logistic regression models were fitted using the “glm” procedure, using age as a continuous variable, sex as a binary variable, number of vaccine doses, BMI category, hemoglobin A1C range, and GFR estimate33 as categorical variables; and presence of hypertension, pulmonary disease or malignancy as binary variables. Receiver operating characteristic (ROC) curves were used to assess model performance34, using 10-fold cross-validation. A two-sided P<0.05 was considered statistically significant. Missing values, which appeared only in BMI, kidney function and hemoglobin A1C variables, were treated by two different approaches. We used k-nearest-neighbors imputation to replace missing values, and also performed an alternative analysis in which missing values were treated as separate “missing” categories. We display here the regression coefficients obtained after imputation. Both methods resulted in very similar classifier performance. ## Results ### Factors associated with hospitalization of SARS-CoV-2 positive individuals 101,039 COVID-19 episodes were included, based on a positive test for SARS-CoV-2 obtained between April 1, 2020 and November 30, 2021. 393 (0.4%) resulted in patient death during hospitalization or within 30 days of the disease, and 1,752 (1.7%) required patient hospitalization for COVID-19 that did not end in patient death. **Table 1A** shows the baseline characteristics of individuals included in the cohort, according to their outcomes. **Table 1B** displays the distribution of categorical variables after imputation of missing values. Generally, patients hospitalized and who died of the disease were older, had a greater proportion of males, had higher BMIs and hemoglobin A1C values, and were more likely to be affected with hypertension, pulmonary disease, malignancy, and impaired kidney function. View this table: [Table 1A:](http://medrxiv.org/content/early/2022/01/03/2021.12.31.21268575/T1) Table 1A: Demographic and clinical characteristics of the study population View this table: [Table 1B:](http://medrxiv.org/content/early/2022/01/03/2021.12.31.21268575/T2) Table 1B: Clinical characteristics of the study population after missing variables imputation We built multivariable logistic regression models to predict both the hospitalization and mortality outcomes. The odds ratios from multivariable regression models reflect the extent with which each risk factor affects the outcome, after adjustment for the others. **Table 2** displays the model for hospitalization risk based on the comparison of 2,145 episodes that resulted either in hospitalization or death vs. 98,894 infections that did not necessitate hospitalization. For each variable, the regression coefficient, with the corresponding odds ratio and 95% confidence interval and p-values are displayed. The footnote explains how to calculate the outcome probability for any given patient data. View this table: [Table 2:](http://medrxiv.org/content/early/2022/01/03/2021.12.31.21268575/T3) Table 2: Logistic regression model for hospitalization risk We emphasize a few key findings arising from the multivariable regression analysis. First, increased age is significantly associated with the risk of hospitalization: each year of age increases the odds for hospitalization by a multiplicative factor of 1.061, which means that compared to an individual aged twenty with similar other risk factors, a patient aged sixty is 10 times more likely to get hospitalized, and a patient aged eighty, 34 times more likely to get hospitalized. Female sex reduces odds for hospitalization by 34%. Obesity increases its risk in a gradual manner (OR=1.324 for BMI between 25 and 30, OR=1.664 for BMI between 30-35, and OR=2.932 for BMI over 35, compared to the reference category of normal BMI, p<0.001 for all). Diabetes mellitus, as reflected by most recent hemoglobin A1C values is independently associated with increased risk in a gradual manner (OR=1.454 for A1C between 6.5 and 8%, OR=1.908 for A1C between 8 and 10%, and OR=3.048 for A1C above 10% compared to the reference category of A1C below 6.5%, p<0.001 for all). Impaired kidney function is also associated with increased risk in a gradual manner (OR=1.568 for GFR between 45 and 59, OR=2.774 for GFR between 30 and 44, and OR=4.000 for GFR below 30 compared to the reference category of GFR above 90, p<0.001 for all). Hospitalization risks significantly increase with hypertension (OR=1.270, 95% CI [1.130-1.428]), pulmonary disease (OR=1.331, 95% CI [1.134-1.563]) and cancer (OR=1.197, 95% CI [1.030-1.390]). As expected, compared to unvaccinated patients, vaccination significantly decreases the hospitalization risk with (OR= 0.602, 95% CI [0.521-0.697] for two vaccine doses and OR=0.339, 95% CI [0.241-0.476] for three vaccine doses, p<0.001 for both categories), and even for the relatively small group of single vaccinated individuals, with OR=0.823, 95% CI [0.694-0.976] (P=0.025). The ROC curve shows the diagnostic ability of a classifier as its discrimination threshold is varied. We performed a 10-fold cross validation and calculated the ROC curve to estimate the performance of our model (**Figure 1**). The performance of the hospitalization risk model was remarkably accurate, with an area under the curve (AUC) of 0.889. The model was able to predict 50%, 80% and 90% of hospitalizations, with respective specificities of 95.3%, 82.2% and 70.2%. ![Figure 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/01/03/2021.12.31.21268575/F1.medium.gif) [Figure 1:](http://medrxiv.org/content/early/2022/01/03/2021.12.31.21268575/F1) Figure 1: Receiver Operating Curve for hospitalization risk model The ROC curve shows the sensitivity and the specificity of the hospitalization model as its discrimination threshold is varied. With a threshold of 8.71% for risk, 50% of the COVID-19 episodes necessitating hospitalization can be identified (sensitivity=50%), and specificity is 95.3% (false positive rate=4.7%); with a risk threshold of 2.39%, sensitivity is 80% and specificity is 82.2% (false positive rate=7.8%); with a risk threshold of 1.14%, sensitivity is 90% and specificity is 70.2% (false positive rate=30.8%). ### Factors associated with mortality for SARS-CoV-2 positive individuals **Table 3** displays the model for mortality risk. It is based on the comparison of 393 fatal cases of COVID-19 compared to 101,039 disease episodes that did not end in patient death. The smaller size of the outcome group limits the power of the model, nevertheless a few factors are associated with a large statistically significant effect: Advanced age is even more strikingly associated with increased mortality risk. Each year of age increases the odds for death by a factor of 1.105: compared to an individual aged twenty with similar other risk factors, a patient aged sixty is 54 times more likely to die of the disease, and a patient aged eighty is 393 times more likely to die. Female sex is also associated with reduced risk of death, reducing this risk by about a half. View this table: [Table 3:](http://medrxiv.org/content/early/2022/01/03/2021.12.31.21268575/T4) Table 3: Logistic regression model for mortality risk Both extreme obesity and underweight increased the risk of death (OR=1.963 for BMI above 35, and OR=2.179 for BMI below 18.5, compared to normal BMI), while other BMI categories were not significantly associated with death. Impaired kidney function was also associated with augmented mortality risk in a gradual manner (OR=2.000 for GFR between 45 and 59, OR=3.097 for GFR between 30 and 44, and OR=6.888 for GFR below 30 compared to the reference category of GFR above 90, p<0.001 for all). Diabetes mellitus, as reflected by last hemoglobin A1C also increases mortality risk in a gradual manner, although more moderately than for hospitalization. The other comorbidities coefficients are overall similar to those for hospitalization, although the smaller outcome group allowed to detect statistical significance only for hypertension and pulmonary disease. Vaccination with a booster dose significantly decreased the mortality risk by 78% (OR= 0.223, 95% CI [0.091-0.551], p=0.001). We performed a 10-fold cross validation and plotted the ROC curve to estimate the performance of the mortality risk model (**Figure 2**). The model was very accurate, with an AUC of 0.967, and was able to predict 50%, 80% and 90% of death events, with respective specificities of 98.6%, 95.2% and 91.2%. ![Figure 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/01/03/2021.12.31.21268575/F2.medium.gif) [Figure 2:](http://medrxiv.org/content/early/2022/01/03/2021.12.31.21268575/F2) Figure 2: Receiver Operating Curve for mortality risk model The ROC curve shows the sensitivity and the specificity of the mortality model as its discrimination threshold is varied. With a threshold of 5.42% for risk, 50% of the COVID-19 episodes ending in patient death can be identified (sensitivity=50%), and specificity is 98.6% (false positive rate=1.4%); with a risk threshold of 1.32%, sensitivity is 80% and specificity is 95.2% (false positive rate =4.8%); with a risk threshold of 0.57%, sensitivity is 90% and specificity is 91.2% (false positive rate=8.8%). ### Risk calculators **Tables 2 and 3** provide all the coefficients, as well as the formula that is used to calculate the absolute risk of a given individual, using the regression models described above. Basically, the coefficients are to be multiplied by the corresponding variables and summed to obtain the natural logarithm of the odds ratio. After exponentiation, the odds ratio can be converted to a probability by dividing it by its value plus one. This calculator is made available online and can be used to calculate hospitalization and mortality probabilities for any given individual. ## Discussion This study developed models that can estimate the risks of subsequent hospitalization and death for any individual newly infected with SARS-CoV-2, using health records from a large healthcare provider in Israel. It also provides further large-scale confirmation of recently published studies from other Israeli health care providers that showed that a third, booster vaccine provides a sharp and almost immediate increase in protection35,36. Considering that vaccination has also been shown to substantially decrease the risk for symptomatic infection, its overall cumulative effect on hospitalization and death are even greater than the odds ratio reported here. Our study capitalizes on the centrally managed, comprehensive electronic health records maintained in our health organization, that are continuously updated with hospitalization and death from COVID-19; the early adoption of a homogenous vaccination program in Israel, and of the booster dose; and the relatively large number of individuals that have been tested positive for COVID-19. We deliberately opted for model simplicity by limiting the number of input variables to recognized clinical factors associated with disease severity that are documented in most health organizations. We intentionally did not include country specific demographic variables, to generate a model that is generalizable to other populations. Even with these limited inputs, the resulting AUCs are highly accurate, achieving 0.889 for hospitalization and 0.967 for mortality. The better AUC performance of the model for mortality risk likely reflects that death is mostly determined by the patient health status, while the decision to hospitalize a patient is additionally impacted by availability of family or social support at home and the patient’s own preferences, which are not accounted for by our model. Importantly, these mortality and hospitalization risk estimates can be given as soon as the disease is diagnosed, allowing to identify which patients are most at risk for a severe course, so that they could be given treatment early in the infection. The models can also prioritize which populations to vaccinate, or urged to receive booster doses, to maximize lives saved and to reduce the load on hospitalization facilities. Several attempts have already been made to build predictive models for COVID-19 severity, notably 4,16,27,37,38. Nevertheless, there is still an unmet need for simple models, which like CURB-6539, could allow instant triage of new patients, using few essential parameters, including vaccination status. This approach allowed us to produce models with remarkably high accuracy. Importantly, in contrast to what we and others found for the infection risk30, we did not observe that time elapsed since vaccination significantly increased hospitalization and mortality risks. For these outcomes, the protective effect of vaccination was largely determined by the cumulative number of vaccine doses received. This may indicate that the immune system response elicited by mRNA vaccine injection has a more lasting effect on hospitalization and mortality risk than on the risk of symptomatic infection following exposure to the virus. Our study has several limitations. First, it is based on a population which was vaccinated almost exclusively with the Pfizer/Biotech BNT162b2 vaccine, and with the first two doses spaced by 21 days. It is uncertain how the estimated effect of vaccination under these conditions would apply to populations vaccinated using different vaccines or using a different vaccination schedule. Moreover, factors specific to our health organization may have affected the results, such as criteria for hospital admission, and treatment decisions that influence mortality. Evolving patient management policies could have a confounding effect on the number of vaccine doses at different times. Last, we have yet to assess the model ability to predict disease severity with new viral variants, such as the recently spreading Omicron. Additional studies in different populations would help to ascertain the validity of our models in different settings. To enable such validation studies, we provide the full model formulas and encourage their use. In conclusion, the models described here, and available online as a calculator, allow to identify individuals most at risk for severe disease or death if infected, using very few essential parameters and vaccination status. This approach can guide public health decisions to optimally allocate vaccines and scarce medicines to maximize lives saved5. Calculator address: [https://covidest.web.app/](https://covidest.web.app/) ## Data Availability Data were obtained from patients' electronic health records and IRB approval restrains its use to researchers inside Leumit Health Services. [https://covidest.web.app/calculator.html](https://covidest.web.app/calculator.html) ## Acknowledgments This research was supported in part by the Intramural Research Program, National Institutes of Health, National Cancer Institute, Center for Cancer Research. The content of this publication does not necessarily reflect the views or policies of the Department of Health and Human Services, nor does mention of trade names, commercial products, or organizations imply endorsement by the U.S. Government. Thanks to Dr. Doug Lowy for helpful comments. The authors had no conflict of interest to report ## Footnotes * Added k-nearest-neighbors imputation for missing data; Minor english corrections ## Abbreviations AUC : Area Under the Curve BMI : Body mass Index COPD : chronic obstructive pulmonary disease COVID-19 : Coronavirus disease of 2019 GFR : glomerular filtration rate LHS : Leumit Health Services OR : odds ratio ROC : Receiver Operating Characteristic SARS-CoV-2 : Severe acute respiratory syndrome coronavirus 2 * Received December 31, 2021. * Revision received January 3, 2022. * Accepted January 3, 2022. * © 2022, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/) ## References 1. 1.Fan Y, Chan KH, Hung IFN. Safety and efficacy of COVID-19 vaccines: A systematic review and meta-analysis of different vaccines at phase 3. Vaccines [Internet] 2021;9(9). Available from: [http://www.ncbi.nlm.nih.gov/pubmed/34579226](http://www.ncbi.nlm.nih.gov/pubmed/34579226) 2. 2.Samaranayake LP, Seneviratne CJ, Fakhruddin KS. Coronavirus disease 2019 (COVID-19) vaccines: A concise review. Oral Dis [Internet] 2021;Available from: [http://www.ncbi.nlm.nih.gov/pubmed/33991381](http://www.ncbi.nlm.nih.gov/pubmed/33991381) 3. 3.Dyer O. Covid-19: FDA expert panel recommends authorising molnupiravir but also voices concerns. Bmj [Internet] 2021 [cited 2021 Dec 28];375:n2984. Available from: [https://pubmed.ncbi.nlm.nih.gov/34857644/](https://pubmed.ncbi.nlm.nih.gov/34857644/) 4. 4.Ioannou GN, Green P, Fan VS, et al. Development of COVIDVax Model to Estimate the Risk of SARS-CoV-2-Related Death among 7.6 Million US Veterans for Use in Vaccination Prioritization. JAMA Netw Open [Internet] 2021;4(4):e214347. Available from: [http://www.ncbi.nlm.nih.gov/pubmed/33822066](http://www.ncbi.nlm.nih.gov/pubmed/33822066) 5. 5.Patel MK. Booster Doses and Prioritizing Lives Saved. N Engl J Med [Internet] 2021;Available from: [http://www.ncbi.nlm.nih.gov/pubmed/34879192](http://www.ncbi.nlm.nih.gov/pubmed/34879192) 6. 6.Persad G, Peek ME, Shah SK. Fair Allocation of Scarce Therapies for COVID-19. Clin Infect Dis [Internet] 2021;Available from: [http://www.ncbi.nlm.nih.gov/pubmed/34922352](http://www.ncbi.nlm.nih.gov/pubmed/34922352) 7. 7.Rosen B, Waitzberg R, Israeli A. Israel’s rapid rollout of vaccinations for COVID-19. Isr J Health Policy Res [Internet] 2021;10(1):6. Available from: [http://www.ncbi.nlm.nih.gov/pubmed/33499905](http://www.ncbi.nlm.nih.gov/pubmed/33499905) 8. 8.Raz A, Keshet Y, Popper-Giveon A, Karkabi MS. One size does not fit all: Lessons from Israel’s Covid-19 vaccination drive and hesitancy. Vaccine 2021;39(30):4027–8. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.vaccine.2021.06.004&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=34120767&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F03%2F2021.12.31.21268575.atom) 9. 9.Dennis JM, Mateen BA, Sonabend R, et al. Type 2 diabetes and covid-19– related mortality in the critical care setting: A national cohort study in england, march–july 2020. Diabetes Care [Internet] 2021;44(1):50–7. Available from: [http://www.ncbi.nlm.nih.gov/pubmed/33097559](http://www.ncbi.nlm.nih.gov/pubmed/33097559) 10. 10.Drucker DJ. Diabetes, obesity, metabolism, and SARS-CoV-2 infection: the end of the beginning. Cell Metab [Internet] 2021;33(3):479–98. Available from: [http://www.ncbi.nlm.nih.gov/pubmed/33529600](http://www.ncbi.nlm.nih.gov/pubmed/33529600) 11. 11.Gimeno-Miguel A, Bliek-Bueno K, Poblador-Plou B, et al. Chronic diseases associated with increased likelihood of hospitalization and mortality in 68,913 COVID-19 confirmed cases in Spain: A population-based cohort study. PLoS One [Internet] 2021;16(11 November):e0259822. Available from: [http://www.ncbi.nlm.nih.gov/pubmed/34767594](http://www.ncbi.nlm.nih.gov/pubmed/34767594) 12. 12.Hamer M, Gale CR, Kivimäki M, Batty GD. Overweight, obesity, and risk of hospitalization for COVID-19: A community-based cohort study of adults in the United Kingdom. Proc Natl Acad Sci U S A [Internet] 2020;117(35):21011–3. Available from: [http://www.ncbi.nlm.nih.gov/pubmed/32788355](http://www.ncbi.nlm.nih.gov/pubmed/32788355) 13. 13.Longmore DK, Miller JE, Bekkering S, et al. Diabetes and overweight/obesity are independent, nonadditive risk factors for in-hospital severity of covid-19: An international, multicenter retrospective meta-analysis. Diabetes Care [Internet] 2021;44(6):1281–90. Available from: [http://www.ncbi.nlm.nih.gov/pubmed/33858854](http://www.ncbi.nlm.nih.gov/pubmed/33858854) 14. 14.Zaki N, Alashwal H, Ibrahim S. Association of hypertension, diabetes, stroke, cancer, kidney disease, and high-cholesterol with COVID-19 disease severity and fatality: A systematic review. Diabetes Metab Syndr Clin Res Rev [Internet] 2020;14(5):1133–42. Available from: [http://www.ncbi.nlm.nih.gov/pubmed/32663789](http://www.ncbi.nlm.nih.gov/pubmed/32663789) 15. 15.Oetjens MT, Luo JZ, Chang A, et al. Electronic health record analysis identifies kidney disease as the leading risk factor for hospitalization in confirmed COVID-19 patients. PLoS One [Internet] 2020;15(11 November):e0242182. Available from: [http://www.ncbi.nlm.nih.gov/pubmed/33180868](http://www.ncbi.nlm.nih.gov/pubmed/33180868) 16. 16.Lian Z, Li Y, Wang W, et al. The Prediction Model of Risk Factors for COVID-19 Developing into Severe Illness Based on 1046 Patients with COVID-19. Emerg Med Int 2021;2021:1–11. 17. 17.Gülsen A, König IR, Jappe U, Drömann D. Effect of comorbid pulmonary disease on the severity of COVID-19: A systematic review and meta-analysis. Respirology [Internet] 2021;26(6):552–65. Available from: [http://www.ncbi.nlm.nih.gov/pubmed/33955623](http://www.ncbi.nlm.nih.gov/pubmed/33955623) 18. 18.Hosseinzadeh R, Goharrizi MASB, Bahardoust M, et al. Should all patients with hypertension be worried about developing severe coronavirus disease 2019 (COVID-19)? Clin Hypertens [Internet] 2021;27(1):3. Available from: [http://www.ncbi.nlm.nih.gov/pubmed/33451360](http://www.ncbi.nlm.nih.gov/pubmed/33451360) 19. 19.Jung Y, Wee JH, Kim JH, Choi HG. The effects of previous asthma and copd on the susceptibility to and severity of covid-19: A nationwide cohort study in south korea. J Clin Med [Internet] 2021;10(20). Available from: [http://www.ncbi.nlm.nih.gov/pubmed/34682749](http://www.ncbi.nlm.nih.gov/pubmed/34682749) 20. 20.Song J, Zeng M, Wang H, et al. Distinct effects of asthma and COPD comorbidity on disease expression and outcome in patients with COVID-19. Allergy Eur J Allergy Clin Immunol [Internet] 2021;76(2):483–96. Available from: [http://www.ncbi.nlm.nih.gov/pubmed/32716553](http://www.ncbi.nlm.nih.gov/pubmed/32716553) 21. 21.Czernichow S, Beeker N, Rives-Lange C, et al. Obesity Doubles Mortality in Patients Hospitalized for Severe Acute Respiratory Syndrome Coronavirus 2 in Paris Hospitals, France: A Cohort Study on 5,795 Patients. Obesity (Silver Spring) [Internet] 2020;28(12):2282–9. Available from: [http://www.ncbi.nlm.nih.gov/pubmed/32815621](http://www.ncbi.nlm.nih.gov/pubmed/32815621) 22. 22.Recalde M, Roel E, Pistillo A, et al. Characteristics and outcomes of 627 044 COVID-19 patients living with and without obesity in the United States, Spain, and the United Kingdom. Int J Obes [Internet] 2021;45(11):2347–57. Available from: [http://www.ncbi.nlm.nih.gov/pubmed/34267326](http://www.ncbi.nlm.nih.gov/pubmed/34267326) 23. 23.Cook TM. The importance of hypertension as a risk factor for severe illness and mortality in COVID-19. Anaesthesia [Internet] 2020;75(7):976–7. Available from: [http://www.ncbi.nlm.nih.gov/pubmed/32339251](http://www.ncbi.nlm.nih.gov/pubmed/32339251) 24. 24.Du Y, Zhou N, Zha W, Lv Y. Hypertension is a clinically important risk factor for critical illness and mortality in COVID-19: A meta-analysis. Nutr Metab Cardiovasc Dis [Internet] 2021 [cited 2021 Dec 21];31(3):745–55. Available from: /pmc/articles/PMC7831720/ 25. 25.Liang X, Shi L, Wang Y, et al. The association of hypertension with the severity and mortality of COVID-19 patients: Evidence based on adjusted effect estimates. J Infect [Internet] 2020 [cited 2021 Dec 21];81(3):e44–7. Available from: /pmc/articles/PMC7315979/ 26. 26.Meng Y, Meng Y, Lu W, et al. Cancer history is an independent risk factor for mortality in hospitalized COVID-19 patients: A propensity score-matched analysis. J Hematol Oncol [Internet] 2020;13(1):75. Available from: [http://www.ncbi.nlm.nih.gov/pubmed/32522278](http://www.ncbi.nlm.nih.gov/pubmed/32522278) 27. 27.Experton B, Tetteh HA, Lurie N, et al. A predictive model for severe covid-19 in the medicare population: A tool for prioritizing primary and booster covid-19 vaccination. Biology (Basel) [Internet] 2021;10(11). Available from: [http://www.ncbi.nlm.nih.gov/pubmed/34827181](http://www.ncbi.nlm.nih.gov/pubmed/34827181) 28. 28.Machado-Alba JE, Valladales-Restrepo LF, Machado-Duque ME, et al. Factors associated with admission to the intensive care unit and mortality in patients with COVID-19, Colombia. PLoS One [Internet] 2021;16(11 November):e0260169. Available from: [http://www.ncbi.nlm.nih.gov/pubmed/34797857](http://www.ncbi.nlm.nih.gov/pubmed/34797857) 29. 29.Vahey GM, McDonald E, Marshall K, et al. Risk factors for hospitalization among persons with COVID-19-Colorado. PLoS One [Internet] 2021;16(9 September):e0256917. Available from: [http://www.ncbi.nlm.nih.gov/pubmed/34473791](http://www.ncbi.nlm.nih.gov/pubmed/34473791) 30. 30.Israel A, Merzon E, Schäffer AA, et al. Elapsed time since BNT162b2 vaccine and risk of SARS-CoV-2 infection: test negative design study. Bmj [Internet] 2021 [cited 2021 Nov 25];375:e067873. Available from: [https://www.bmj.com/content/375/bmj-2021-067873](https://www.bmj.com/content/375/bmj-2021-067873) 31. 31.Rennert G, Peterburg Y. Prevalence of selected chronic diseases in Israel. Isr Med Assoc J [Internet] 2001;3(6):404–8. Available from: [http://www.ncbi.nlm.nih.gov/pubmed/11433630](http://www.ncbi.nlm.nih.gov/pubmed/11433630) 32. 32.Hamood R, Hamood H, Merhasin I, Keinan-Boker L. A feasibility study to assess the validity of administrative data sources and self-reported information of breast cancer survivors. Isr J Health Policy Res 2016;5(1):50. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s13584-016-0111-6&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27980719&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F03%2F2021.12.31.21268575.atom) 33. 33.Levey AS, Coresh J, Greene T, et al. Using standardized serum creatinine values in the modification of diet in renal disease study equation for estimating glomerular filtration rate. Ann Intern Med 2006;145(4):247–54. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.7326/0003-4819-145-4-200608150-00004&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16908915&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F03%2F2021.12.31.21268575.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000239858800002&link_type=ISI) 34. 34.Lasko TA, Bhagwat JG, Zou KH, Ohno-Machado L. The use of receiver operating characteristic curves in biomedical informatics. J Biomed Inform [Internet] 2005;38(5):404–15. Available from: [http://www.ncbi.nlm.nih.gov/pubmed/16198999](http://www.ncbi.nlm.nih.gov/pubmed/16198999) 35. 35.Arbel R, Hammerman A, Sergienko R, et al. BNT162b2 Vaccine Booster and Mortality Due to Covid-19. N Engl J Med [Internet] 2021;Available from: [http://www.ncbi.nlm.nih.gov/pubmed/34879190](http://www.ncbi.nlm.nih.gov/pubmed/34879190) 36. 36.Bar-On YM, Goldberg Y, Mandel M, et al. Protection of BNT162b2 Vaccine Booster against Covid-19 in Israel. N Engl J Med [Internet] 2021;385(15):1393–400. Available from: [http://www.ncbi.nlm.nih.gov/pubmed/34525275](http://www.ncbi.nlm.nih.gov/pubmed/34525275) 37. 37.Mamidi TKK, Tran-Nguyen TK, Melvin RL, Worthey EA. Development of An Individualized Risk Prediction Model for COVID-19 Using Electronic Health Record Data. Front Big Data [Internet] 2021;4:675882. Available from: [http://www.ncbi.nlm.nih.gov/pubmed/34151259](http://www.ncbi.nlm.nih.gov/pubmed/34151259) 38. 38.Hippisley-Cox J, Coupland CA, Mehta N, et al. Risk prediction of covid-19 related death and hospital admission in adults after covid-19 vaccination: national prospective cohort study. BMJ [Internet] 2021;374:n2244. Available from: [http://www.ncbi.nlm.nih.gov/pubmed/34535466](http://www.ncbi.nlm.nih.gov/pubmed/34535466) 39. 39.Lim WS, Van Der Eerden MM, Laing R, et al. Defining community acquired pneumonia severity on presentation to hospital: an international derivation and validation study. Thorax [Internet] 2003 [cited 2021 Dec 31];58(5):377. Available from: /pmc/articles/PMC1746657/?report=abstract