key: cord-0940898-m8p9xnfu authors: Ronderos Botero, Diana Maria; Omar, Alaa Mabrouk Salem; Sun, Haozhe Keith; Mantri, Nikhitha; Fortuzi, Ked; Choi, Yongsub; Adrish, Muhammad; Nicu, Marin; Bella, Jonathan N.; Chilimuri, Sridhar title: COVID-19 in the Healthy Patient Population: Demographic and Clinical Phenotypic Characterization and Predictors of In-Hospital Outcomes date: 2020-09-10 journal: Arterioscler Thromb Vasc Biol DOI: 10.1161/atvbaha.120.314845 sha: 428570f58bcb46fd7f6bf2a9b23f5a4184f74999 doc_id: 940898 cord_uid: m8p9xnfu Coronavirus disease 2019 (COVID-19) can infect patients in any age group including those with no comorbid conditions. Understanding the demographic, clinical, and laboratory characteristics of these patients is important toward developing successful treatment strategies. APPROACH AND RESULTS: In a retrospective study design, consecutive patients without baseline comorbidities hospitalized with confirmed COVID-19 were included. Patients were subdivided into ≤55 and >55 years of age. Predictors of in-hospital mortality or mechanical ventilation were analyzed in this patient population, as well as subgroups. Stable parameters in overall and subgroup models were used to construct a cluster model for phenotyping of patients. Of 1207 COVID-19–positive patients, 157 met the study criteria (80≤55 and 77>55 years of age). Most reliable predictors of outcomes overall and in subgroups were age, initial and follow-up d-dimer, and LDH (lactate dehydrogenase) levels. Their predictive cutoff values were used to construct a cluster model that produced 3 main clusters. Cluster 1 was a low-risk cluster and was characterized by younger patients who had low thrombotic and inflammatory features. Cluster 2 was intermediate risk that also consisted of younger population that had moderate level of thrombosis, higher inflammatory cells, and inflammatory markers. Cluster 3 was a high-risk cluster that had the most aggressive thrombotic and inflammatory feature. CONCLUSIONS: In healthy patient population, COVID-19 remains significantly associated with morbidity and mortality. While age remains the most important predictor of in-hospital outcomes, thromboinflammatory interactions are also associated with worse clinical outcomes regardless of age in healthy patients. GRAPHIC ABSTRACT: A graphic abstract is available for this article. C oronavirus disease 2019 is caused by the recently identified coronavirus also known as severe acute respiratory syndrome coronavirus 2 (or SARS-CoV-2). 1 Since December 2019, COVID-19 has caused a global outbreak and is still spreading quickly in >100 countries. 2 Our understanding of this disease has been increasing with time. A disease entity that was initially identified as a primarily respiratory illness has slowly emerged as a systemic syndrome that causes endothelial dysfunction leading to microthrombosis and severe inflammatory response leading to a cytokine storm. 3 However, the pathophysiological and clinical characterization of the disease is still evolving. The disease still remains significantly heterogeneous both in its demographic characteristics and clinical features. One example of such heterogeneity is that this disease was thought to vastly affect the older population with significant comorbidities such as hypertension and diabetes mellitus. 4 However, in our clinical experience, we have started to see significant morbidity and mortality in younger and otherwise healthy patient population. With Patients with symptoms suggestive of COVID-19 infection were identified in the emergency department. After initial assessment, basic clinical data including time from symptom onset and admission laboratory investigations were obtained. In addition, all patients had pulmonary imaging (chest radiograph or computed tomography of the chest) done. Patients subsequently had a nasal swab for COVID-19 RNA to confirm the infection with SARS-Cov-2. Electronic medical charts were reviewed for the presence of baseline comorbidities, initial and subsequent laboratory values, treatment and therapeutic modalities, need for mechanical ventilation, and in-hospital mortality. All laboratory data were obtained on the same day of admission, while follow-up laboratory values were obtained for LDH (lactate dehydrogenase) on the second day of admission and for d-dimer if a change in the clinical condition (change in the requirement of oxygen or vital signs or clinical suspension of thromboembolic episode). Patients were then classified based on their age as younger (≤55 years of age) or older (>55 [older] years of age). We selected a 55-year-old cutoff value to differentiate between younger and older patients based in the US Census Bureau report, which defines older adults as ≥55 years of age and elderly as ≥65 years of age. 5 Subgroups were compared for their demographic, clinical, and laboratory variables and outcomes. The primary end point was in-hospital mortality and need for mechanical ventilation. Continuous variables were expressed as mean±SD, and nominal and categorical variables were expressed as numbers (%). The independent samples Student t test and 1-way ANOVA were used to compare the mean values of different groups, and χ 2 test was used for comparison of nominal and categorical variables. Predictors of mortality and mechanical ventilation were checked using univariate and multivariate Cox regression models. The most stable predictors in all models were used to construct an SPSS unsupervised 2-step cluster model to test the presence of the natural subgroups and were subsequently compared for characterization of each cluster of patients. Kaplan-Meier survival curves were used to test the difference in cumulative in-hospital outcomes in different clusters. For all statistical tests, P<0.05 was considered statistically significant. All analyses were performed with commercially available software (SPSS, version 23.0; SPSS, Inc). During the study period, 1207 patients with confirmed COVID-19 were identified. The baseline demographic, clinical, laboratory, and imaging criteria for all patients are summarized in Table I One hundred fifty-seven patients met the study criteria and were included for final analysis. The mean age of the study population was 52.6±17 years. Eighty (50.9%) patients were ≤55 years of age, and 46 (29%) were women. Of the study population, 105 (67%) patients received hydroxychloroquine/azithromycin combination, and 13 (8%) patients received When patients were classified based on their age, it was found that patients >55 years of age had lower lymphocyte count, and lower albumin levels. Older patients were also noted to have higher d-dimer, LDH (lactate dehydrogenase), lactate, ferritin, and creatinine levels ( Table 1 ). There were no differences between the two groups with regard to sex, body mass index, and laboratory values including high-sensitivity C-reactive protein, hemoglobin levels, white cell count, platelet count, prothrombin time, partial thromboplastin time, bilirubin, serum protein, and procalcitonin levels ( Table 1 ). The two groups were similar in terms of medication use. Older patient subgroup had higher mortality and higher need for mechanical ventilation (Table 1) . In a subgroup analysis of patients who died, we noted that these patients were older and had higher d-dimer, LDH, lactate, and procalcitonin levels ( Table 1) . The overall predictors of composite of in-hospital mortality or need for mechanical ventilation in patients without comorbidities are summarized in Table 2 . Univariate Cox regression model showed that the independent predictors of in-hospital mortality were age, d-dimer, and LDH levels. Multivariate Cox regression models, however, showed that d-dimer and LDH levels were the only predictors of in-hospital mortality, while age lost statistical significance. In patients ≤55 years of age with no comorbidities (younger healthier patients), the predictors of in-hospital mortality in the univariate Cox regression model were lower hemoglobin, higher neutrophil count, d-dimer, LDH, and total bilirubin level. In the multivariate regression analysis, neutrophil count, d-dimer, and LDH levels were the only independent predictors of in-hospital mortality ( Table 2) . In older patients' subgroup, higher d-dimers, LDH, and partial thromboplastin time levels were the predictors of in-hospital mortality by univariate Cox regression model. Multivariate regression analysis, however, showed that only higher d-dimer predicted in-hospital mortality ( Table 2) . Most consistent predictors of the composite outcomes regardless of age in all models were found to be d-dimer and LDH. Follow-up values of d-dimer and LDH were found to be significantly higher in patients >55 years of age and in patients who died, compared with patients ≤55 years of age and those who survived (P=0.02 and 0.005, respectively; Table 1 ). We also analyzed Δd-dimer and ΔLDH and noted that they were not significantly different between younger and older patient subgroups (P=0.390 and 0.356, respectively). However, these values were significantly higher in patients who died compared with patients who survived (P=0.002 and 0.026, respectively; Table 1 ). Receiver operating characteristic was done for initial and follow-up (second) d-dimer and LDH, as well as their delta values (Figure 1 ), and we found that the best predictors for composite outcomes for initial, follow-up, and Δd-dimer were when their values were >461, 491, and 798 ng/mL, respectively. Similarly, the best predictors for initial, follow-up, and ΔLDH were when their values were >467, 505, and 128 units/L, respectively. Based on these cutoff values, Cox regression models were repeated with adjustment for both age and sex for all patients (Table 3) , and it was noted that all adjusted variables can significantly predict outcomes based on univariate models, whereas in multivariate models, only adjusted Δd-dimer >798 ng/mL and repeat LDH level of >505 units/L predicted mortality. Outcome Based on Age, D-Dimer, and LDH Next, a 2-step cluster model was initiated for all patients based on age, initial d-dimer >461 ng/mL, follow-up d-dimer >491 ng/mL, initial LDH >467 unit/L, and follow-up LDH >505 unit/L ( Figure 2 ). The output revealed that patients were classified into 3 different clusters. Based on Kaplan-Meier curve for survival free of the composite outcomes, it was noted that cluster 1 was a lowrisk cluster, cluster 2 was an intermediate-risk cluster, and cluster 3 was a high-risk cluster (mortality and mechanical ventilation: 0 [0%], 16 [27%], and 11 [55%], respectively; P<0.001; Table 4 ). Post hoc examination of the clusters revealed the following characteristics for each cluster: Age Clusters 1 and 2 were significantly younger than cluster 3 (51.4±14.7, 54.3±15.9, and 64.8±8.7 years of age, respectively; P=0.005; Table 4 ). Thrombosis burden based on d-dimer level was found to progressively increase from cluster 1 to clusters 2 and 3 (260±126, 2133±5053.3, and 13332±19692 ng/mL, respectively; P<0.001; Table 4 ). While hemoglobin and platelet levels were not different between clusters, white cell count was the lowest in cluster 1 compared with the other two clusters (6.9±3.6, 8.7±3.7, and 10.7±4.4 10/ 3 mL, respectively; P=0.002; Table 4 ). Deferential white count showed that this difference was mainly due to neutrophil count, which progressively increased from cluster 1 to cluster 3 (5.4±3.6, 7.4±3.6, and 9.1±4.1 10 3 /mL, respectively; P=0.002; Table 4 ), while lymphocyte count was not different between clusters. It was noted that for acute-phase reactants such as highsensitivity C-reactive protein and ferritin significantly and progressively increased from cluster 1 to cluster 3 (highsensitivity C-reactive protein: 104.7±105, 166±101, and 158±80 mg/d, respectively, P=0.024; ferritin: 756±669, 1408±1097, and 1880±1938 mg/mL, respectively, P=0.006; Table 4 ). Parameters suggestive of target organ damage such as renal function tests and liver function tests and parameters of tissue hypoperfusion such as lactate levels were checked for all clusters. Target organ damage and tissue hypoperfusion were the hallmark of cluster 3 as suggested by difference in creatinine (cluster 1, 1±0. 2 Table 4 ). In addition, tissue damage as suggested by LDH levels progressively increased from cluster 1 to clusters 2 and 3 (cluster 1, 315±87 unit/L; cluster 2, 545±187 unit/L; cluster 3, 822±395 unit/L; P<0.001; Table 4 ). To the best of our knowledge, this is the first report describing the demographic, clinical, and laboratory parameters of healthier patient population hospitalized with COVID-19. The main findings of the current study can be summarized as follows: the main predictors of in-hospital outcomes (mortality or need for mechanical ventilation) are older age and high d-dimer and LDH levels. d-dimer and LDH seemed to be the most consistent predictors of adverse outcomes regardless of age. Moreover, the predictive ability of both variables increased when these were followed up during hospitalization with the best sensitivity acquired by second high d-dimer measurement and the best specificity observed for significant positive change of d-dimer (Δd-dimer) on subsequent testing. Finally, when patients were classified with an unsupervised cluster using their age, initial, and follow-up d-dimer, and LDH, the model output created 3 different phenotypes with distinct in-hospital outcomes. Importantly, two of these clusters were relatively younger patients who were differentiated based on their d-dimer, inflammatory markers, and parameters of end-organ damage, and the third cluster consisted of older patients with the significantly worst clinical and laboratory profile. Despite the great experience gained in the clinical characterization and diagnosis, COVID-19 remains a dynamic disease process for which our understanding continues to evolve from what was initially identified as a primarily respiratory illness toward a systemic syndrome that involves multiorgan endothelial dysfunction strongly linked to microthrombosis and severe inflammatory response. The heterogeneity of the disease is demonstrated by the fact that it continues to significantly affect all age groups and not just the older population or those with comorbidities, particularly diabetes mellitus and hypertension. 6 Recent studies showed that children and young adults who died had also had baseline comorbidities and that those who did not have baseline comorbidities could be managed on general wards and were discharged successfully. 7, 8 Demographic, Laboratory, and Clinical Characteristics of the Disease In our study, of 1207 patients admitted with COVID-19 during the earlier period of the pandemic, 157 of them (13%) had no baseline comorbidities, which was slightly more than that reported by the Centers for Disease Control (7.9%). 9 Twenty-seven (17%) of all patients reported here had in-hospital outcomes. Importantly, 80 of these patients were <55 years of age, of whom 7 (9%) had in-hospital outcomes, suggesting that rates of serious adverse outcomes are not as low as initially expected at the beginning of the pandemic. While age remained an important predictor of mortality, sex failed to predict mortality overall and in subgroups. When patients were classified based on their age, it was found that older patients (>55 years of age) had lower lymphocyte count, serum albumin, and platelet count and higher LDH, creatinine, and lactic acid, and the most remarkable difference between both groups was the significantly elevated d-dimer levels in the older group. This was confirmed in the Cox regression models when only d-dimer and LDH were the stable predictors of mortality regardless of the age. This suggests that the primary pathophysiology in this group of patients is thrombotic in nature while all other changes can be secondary. The fact that the follow-up d-dimer and Δd-dimer levels in hospital were the parameters that could predict mortality with high sensitivity and specificity, respectively, suggests that this thrombotic burden is dynamic and continues to rebound during the course of the disease. When unsupervised cluster model was applied to all patients, an important link between the inflammatory status and thrombosis was uncovered pointing toward age-related progressive worsening of inflammatory state coupled with the inflammatory cells and the inflammatory markers that seems linked to higher mortality. It is important to note that such natural classification of patients revealed that younger patients may not be similar in their presentation in the disease and only those who have higher thromboinflammatory burden are at an increased risk for in-hospital outcomes, suggesting that there may be underlying physiological or genetic factors that enhance such process in some individuals versus the others. On the other hand, the cluster model suggested that in older age, this thromboinflammatory relationship seems to be the most aggressive, leading to higher in-hospital outcomes. SARS-Cov-2 entry into the cells is mediated by ACE-2 (angiotensin-converting enzyme 2) receptor on cellular membranes, 10 abolishing ACE-2 enzyme activity. 11 The expression of ACE-2 receptors may decrease with advanced age. 12 This may lead to a decreased metabolism of angiotensin II and a systemic imbalance between angiotensin II (proinflammatory/procoagulant) and angiotensin 1-7 (anti-inflammatory/anticoagulant), causing a worse baseline proinflammatory/procoagulant state in older patients that is aggravated by COVID-19 infection, which induces further downregulation of ACE-2 in older patients (Figure 3 ). The aging process can also impact the immune system leading to a process of a less-competent immune system referred to as immune senescence, 13 which incorporates dysfunction in cellular signaling that directly affects the ability to contain infections and in case of COVID-19 can lead to dysfunctional (decreased or overexerted) response to viral load and inability to control the viral spread. Another possible explanation is inflammation-induced thrombosis (IIT)-a phenomenon reported in several inflammatory and autoimmune diseases. 14 The pathogenesis of IIT is complicated and involves bidirectional interactions between the inflammatory state and the coagulation system ( Figure 3 ). Reports have shown that IIT can be initiated either on the proinflammatory end or the procoagulant end. An overlap of both ends of IIT may occur, causing a vicious cycle that leads to rebound inflammation and thrombosis with a baseline increased proinflammatory state as explained previously. This may explain why in our study the predictive ability of d-dimer elevation may shed light on that a re-elevation of d-dimer may signify multiple episodes of thrombosis that is linked to high mortality in those patients. Importantly, such IIT is also reported to develop in the absence of vascular endothelial damage, which may explain the pathology in younger patients observed in our study. Recent pathological reports have shown that viral elements exist within endothelial cells of patients with COVID-19 associated with accumulation of inflammatory cells, causing endotheliitis. 15 Endothelial dysfunction occurring in COVID-19 can shift the vascular equilibrium toward more vasoconstriction, subsequent organ hypoperfusion, and inflammation with associated procoagulant state. This can explain the elevated lactate, LDH, and inflammatory markers that seem to be related to the severity of the disease in our study. In patients with COVID-19, endotheliitis can lead to activation of the coagulation cascade to form fibrin meshes in an attempt to contain the viral spread by forming microthrombi in situ, 16 which can potentially dislodge systemically leading to arterial and venous thrombosis. As such, the integrity of the endothelial function and immune system in the younger people is crucial for a regulated immune response good enough to control the virus spread by forming microthrombi for viral entrapment, and this in turn may partially lead to a relative quiescent systemic inflammatory state in the younger healthier patients. On the other hand, direct endothelial infection can also aggravate the baseline endothelial dysfunction that can exist because of comorbidities or due to aging, 17 which may partially explain the high levels of d-dimers and inflammatory response in the older patients in our study. An alternate explanation for the change in d-dimer with time can be the multiple injuries, or multiple hit to the endothelium by different successive mechanisms while the disease progresses. The successive hits can be composed of endothelial infection with the associated thrombo-entrapment of the virus, inflammatory response to the virus mediated by IIT, and hypoxia-or hypotensioninduced thrombosis. 18 In our study, we have attempted to uncover the natural distribution of parameters between patients to understand phenotyping, and for that, we have used a machine learning-based unsupervised cluster model. Computational clustering is an exploratory statistical method designed to uncover natural groups within data that would otherwise be invisible using traditional classification methods. Unsupervised clustering, particularly, separates patients into groups without a priori classification. In our study, unsupervised 2-step cluster analysis was used. In the first step, small preclusters were created, which are then merged according to the greatest change in the distance measure in the second step into the most meaningful clusters. The parameters used for clustering in our study were those that had the best prediction of mortality in terms of Cox regression models and receiver operating characteristic curves, namely age, initial and follow-up d-dimer, and LDH levels. The model spit out 3 distinct clusters in terms of in-hospital outcomes that could be identified as low, intermediate, and high risk. The most interesting finding was that the lowand intermediate-risk clusters were identical in terms of age (both mostly represent younger patients) and other demographic variables; however, the low-risk cluster had no in-hospital outcomes and was characterized by the better profile across all predefined disease domains, namely thrombosis, inflammation, inflammatory cells, and end-organ damage. On the other hand, the intermediaterisk cluster showed an intermediate level of thrombosis, elevation of inflammatory cells, and target organ damage and a higher level of inflammatory markers. The highrisk cluster had the worst profile across all domains and had older age. Such clusters seem reasonably explained based on the previously proposed mechanisms. It seemed that the presumed more preserved endothelial function in the intermediate-risk group leads to a less thrombosis burden compared with the high-risk group in which the mechanisms of thromboinflammation were the most aggressive probably because of a worse state of endothelial dysfunction created by more pronounced aging process as explained previously. We acknowledge the following limitations: first, the sample size is small, which can affect the results of all models. However, the incidence and prevalence of COVID-19 in patients with no comorbidities remain small, making it difficult to gather information from an appropriate number of patients. As such, the findings of our study should be confirmed by further larger multicenter studies. Second, the data and variables collected for our patients depended on the laboratory analyses that were done in the initial stages of the disease. Other relevant laboratory findings that can confirm our hypothesis such as troponin levels, and longer trends of d-dimers and LDH, as well as echocardiography, computed tomography chest with contrast, or other imaging modalities, were not available for the vast majority of the patients. Moreover, the time between the first and second d-dimer evaluation was based on clinical suspicion only, as our study was not powered to detect the effect of time to change of d-dimer clinical outcomes; further studies need to address this limitation. Further studies should focus on recruiting patients after measuring these variables to get more insight into the pathophysiological mechanisms. Third, most of the presumed mechanisms in the current study are just assumptions, and despite that pathological reports exist, they are still rarely done, and to confirm the presumed mechanisms, more detailed pathological reports are needed. Fourth, the 2-step cluster model that was used in our study functions best in the presence of large scale data, and as such, results from this model should be treated with care. Traditionally, sample size selection is based on the number of variables fed into the cluster model. The best approach to assess the appropriateness of the sample size is to determine whether the dimensionality of the data is not too high for the number of cases to be grouped. The minimal sample size to include should be no <2 k cases (k=number of variables), preferably 5×2 k . 19, 20 In our analysis, the number of variables fed into the model was 5, and accordingly, the appropriate sample size for a relatively stable model would be anywhere between 32 and 160. Thus, our sample size of 157 patients would be sufficient. However, validation of all the models described in this report needs to be done on larger scale studies before any of these findings can be adopted into clinical practice. Fifth, the differences between the laboratory values in the conventional analysis (Table 1) , despite statistically significant, were only subtle, and the mean values of parameters such as those of neutrophils, lymphocytes, creatinine, and liver function tests appear to be normal. However, considerable number of patients had abnormal values of these parameters when categorically classified (Table II in the Data Supplement) . It is to be noted that when patients were classified based on the unsupervised cluster model using d-dimer, LDH, and age, the differences between these variables changed, demonstrating how a cluster analysis without a priori assumptions can change our look into laboratory results based on the natural distribution of the classifying attributes, compared with conventional methods of comparison, which depends on artificial separation of data based on arbitrary cutoff values such as age of 55 years used in the current study. Finally, while d-dimer elevates as the result of thrombolysis, which represents an indirect evidence of thrombosis and thus used clinically to exclude thromboembolic episodes, d-dimer can also be elevated in a variety of other situations including, pregnancy, inflammation, malignancy, trauma, postsurgical treatment, liver disease (decreased clearance), and heart disease. [21] [22] [23] It is also frequently high in hospitalized patients. Given that the population studied is healthy at baseline excludes most of the nonthrombotic causes of d-dimer elevation and, moreover, the high levels for patients on presentation and follow-up during their hospital stay, is somewhat not familiar to those noted for hospitalized patients. In patients without baseline comorbidities, COVID-19 can be fatal and is associated with significant in-hospital outcomes defined as in-hospital mortality or need for mechanical ventilation. While age remains the most important predictor, d-dimer and LDH elevation at baseline and follow-up during the hospital stay seems to be linked to worse in-hospital outcomes regardless of age in these patients. The underlying mechanisms seem to be related to an aging process, IIT, endothelial infection, or multiple hit from different successive pathological insults. Understanding of such pathological processes and their impact on patients with and without comorbidities among all age categories is crucial toward developing successful treatment strategies. Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and coronavirus disease-2019 (COVID-19): the epidemic and the challenges Global trends in epidemiology of coronavirus disease 2019 (COVID-19) Weathering the COVID-19 storm: lessons from hematologic cytokine syndromes Prevalence and impact of cardiovascular metabolic diseases on COVID-19 in China Classifying age groups D-dimer levels on admission to predict in-hospital mortality in patients with COVID-19 Coronavirus disease 2019 in children -United States Chinese Pediatric Novel Coronavirus Study Team. SARS-CoV-2 infection in children COVIDView: a weekly surveillance summary of U.S. COVID-19 activity SARS-CoV-2 cell entry depends on ACE2 and TMPRSS2 and is blocked by a clinically proven protease inhibitor Angiotensin-converting enzyme 2 (ACE2) as a SARS-CoV-2 receptor: molecular mechanisms and potential therapeutic target Ronderos Botero et al COVID-19 in the Healthy Patient Population The dilemma of coronavirus disease 2019, aging, and cardiovascular disease: insights from cardiovascular aging science Oxidative-inflammatory stress in immune cells from adult mice with premature aging Inflammation-induced thrombosis: mechanisms, disease associations and management Endothelial cell infection and endotheliitis in COVID-19 Neutrophil extracellular traps: a new source of tissue factor in atherothrombosis Endothelial dysfunction and aging: an update Pulmonary vascular endothelialitis, thrombosis, and angiogenesis in COVID-19 Die latent-class-analyse: Einführung in theorie und anwendung A review of unquestioned standards in using cluster analysis for data-driven market segmentation Towards evidence based emergency medicine: best BETs from the manchester royal infirmary. BET 2: D-dimer levels during normal menstrual cycle Relation of D-dimer and troponin I in patients with new-onset atrial fibrillation The influence of various patient characteristics on D-dimer concentration in critically ill patients and its role as a prognostic indicator in the intensive care unit setting None.