The emergence of COVID-19 in Indonesia: analysis of predictors of infection and mortality using independent and clustered data approaches ========================================================================================================================================= * Erlina Burhan * Ari Fahrial Syam * Ahmad Jabir Rahyussalim * Prasenohadi * Navy G Lolong Wulung * Agus Dwi Susanto * I Gede Ketut Sajinadiyasa * Dewi Puspitorini * Dewi Lestari * Indah Suci Widyahening * Vivi Setiawaty * Dwiana Ocviyanti * Kartika Qonita Putri * Aswin Guntara * Davrina Rianda * Anuraj H Shankar * Rina Agustina ## Abstract **Background** Analyses of correlates of SARS-CoV-2 infection or mortality have usually assessed individual predictors. This study aimed to determine if patterns of combined predictors may better identify risk of infection and mortality **Methods** For the period of March 2nd to 10th 2020, the first 9 days of the COVID-19 pandemic in Indonesia, we selected all 18 confirmed cases, of which 6 died, and all 60 suspected cases, of which 1 died; and 28 putatively negative patients with pneumonia and no travel history. We recorded data for travel, contact history, symptoms, haematology, comorbidities, and chest x-ray. Hierarchical cluster analyses (HCA) and principal component analyses (PCA) identified cluster and covariance patterns for symptoms or haematology which were analysed with other predictors of infection or mortality using logistic regression. **Results** For univariate analyses, no significant association with infection was seen for fever, cough, dyspnoea, headache, runny nose, sore throat, gastrointestinal complaints (GIC), or haematology. A PCA symptom component for fever, cough, and GIC tended to increase risk of infection (OR 3.41; 95% CI 1.06-14; p=0.06), and a haematology component with elevated monocytes decreased risk (OR 0.26; 0.07-0.79; 0.027). Multivariate analysis revealed that an HCA cluster of 3-5 symptoms, typically fever, cough, headache, runny nose, sore throat but little dyspnoea and no GIC tended to reduce risk (aOR 0.048; <0.001–0.52; 0.056). In univariate analyses for death, an HCA cluster of cough, fever and dyspnoea had increased risk (OR 5.75; 1.06 − 31.3, 0.043), but no other individual predictor, cluster or component was associated. Other significant predictors of infection were age ≥ 45, international travel, contact with COVID-19 patient, and pneumonia. Diabetes and history of contact were associated with higher mortality. **Conclusions** Cluster groups and co-variance patterns may be stronger correlates of SARS-CoV-2 infection than individual predictors. Comorbidities may warrant careful attention as would COVID-19 exposure levels. Keywords * Coronavirus * SARS-CoV-2 * COVID-19 * hierarchical cluster analysis * Indonesia * principal component analysis ## Introduction The coronavirus SARS-CoV-2 is a deadly respiratory pathogen first reported in Wuhan City, Hubei Province of China in early 2020, and rapidly spread globally. On 11 March 2020, the World Health Organization declared a pandemic(1), and by 2 May 2020 the number of confirmed cases worldwide reached nearly 3.3 millions with more than 230,000 deaths, and a global case fatality rate of 7.04%.(2) Indonesia, a large archipelago country in South-east Asia with a population of 260 million, did not report a confirmed COVID-19 case until 2 March 2020, several weeks after a total of 217 cases had been reported from neighboring countries of Singapore, Malaysia, Thailand and Australia. This raised concerns from groups both inside and outside the country regarding diagnostic and reporting practice. However, by 2 May 2020 the number of confirmed cases were 10,551, with 800 deaths, yielding a case fatality rate of 7.58%, amongst the highest in the world.(2) Again, this raised concerns about both the diagnostic capacity and care of cases. The availability of routine data for case history, symptoms, haematology and radiology presents an opportunity to define presentation patterns predictive for COVID-19. Knowledge of anomalous risk patterns may be especially useful at the early phase of the epidemic in new locations when high throughput molecular testing is not yet available. Analyses of syndromic and haematology data in some studies showed the most common symptoms of COVID-19 were fever, cough, and lymphocytopenia. Patients may also be asymptomatic yet with abnormal findings from lung imaging.(3) Previous analyses typically examined each symptom, exposure or hematological characteristic individually or as separate predictors of SARS-CoV-2 infection, pathogenicity, or mortality. Few studies have examined clusters or co-variance patterns of symptoms, exposures, haematological data, or combinations thereof, and their associations with outcomes. We therefore analysed routine clinical data and patient history from all hospital admissions during the first nine days of the epidemic in Indonesia who were confirmed by real time polymerase chain reaction (RT-PCR) as COVID-19 cases, and all suspect but virus negative cases, and other lung infection cases presumed negative for the virus. We applied hierarchical cluster analysis (HCA) and principal components analysis (PCA) to discern specific clusters or co-variance patterns of syndromic and haemotological data, and analysed these along with radiological features and patient history for associations with SARS-CoV-2 infection or death using univariate and multivariate logistic regression. These analytic approaches may prove useful for more rapid recognition of emergent COVID-19 cases, thereby hastening treatment, and enhancing survival. ## Materials and methods ### Study setting This was a retrospective cohort study of patients from four collaborating hospitals in Jakarta, Bali and Jambi provinces which had reported suspected and confirmed case of COVID-19 to the Ministry of Health. We collected data from confirmed, suspected and lung infection cases admitted from March 2nd to 10th 2020, the first nine days of cases in Indonesia. Ethical approval was obtained from the Research Ethics Committee of the Faculty of Medicine, Universitas Indonesia – Dr. Cipto Mangunkusumo General Hospital, with approval number 20030331. After consideration of ethical issues, logistics and urgency of this work, the Committee waived the requirement for written individual informed consent and approved the sharing of anonymized data. ### Patients identification As per hospital protocols, patients admitted with respiratory infections had been registered as potential COVID-19 cases and been classified as under-supervision or monitored. Under-supervision had included patients with clinical symptoms of fever (>38 °C) or history of fever, runny nose, sore throat, and pneumonia-like symptoms such as cough, dyspnoea, sweating, sharp chest pain worsened by breathing deeply, or cough with a history of travel to a COVID-19 affected country in the last 14 days, or having contact with a person with acute upper respiratory tract infection. Monitored had been defined as upper respiratory symptoms without pneumonia-like symptoms, plus travel history to an affected country within 14 days. A confirmed COVID-19 case was an under-supervision or monitored patient with a positive SARS-CoV-2 test by RT-PCR for naso-oropharyngeal swab, sputum, or bronchoalveolar lavage (18 patients); a suspected case was a patient twice-confirmed negative by RT-PCR test (47 patients) or not assessed (13 patients); and formal classification required the case had been reported to the Ministry of Health of the Republic Indonesia. For comparison, we selected patients with community acquired pneumonia who were admitted between January 1st to March 10th 2020 and reported no history of international travel. A total of 18 confirmed COVID-19 patients, 60 suspected patients, and 28 lung patients were identified. ### Data collection Study physicians reviewed medical records of hospitalized patients, and in accord with World Health Organization diagnostic guidelines, and hospital procedures, data were extracted for demographics, nationality, travel history, clinical signs, radiological features, and haemotology, including haemoglobin and hematocrit, and counts for thrombocytes, total leucocytes, monocytes, lymphocytes, basophils, neutrophils, and eosinophils. Patients with symptoms of acute respiratory tract infection or travel history from COVID-19-affected countries within the last 14 days had been assessed in the emergency room or out-patient clinic(4), and case history, physical examination and chest x-rays had been done within the hospital. Venous blood samples had been taken for complete white blood cell count. Some samples were assessed for electrolytes, lactic acid products, C-reactive protein, procalcitonin, and liver as well as renal function tests; these tests were selectively assessed and gaps in data precluded inclusion in this analysis. Chest computed tomography (CT) scan was not reported for any patient.(5) Naso-oropharyngeal swab, sputum or bronchoalveolar lavage specimens had been sent to the National Institutes for Health Research and Development (NIHRD) laboratory in Jakarta for RT-PCR to assess presence of the virus. Tests were performed after virus inactivation, and RNA was extracted using the QiAmp RNA viral kit based on the instructions and standard laboratory techniques, followed by RT-PCR using SuperScript™ III Platinum® One-Step Quantitative RT-PCR for the N gene in accord with US CDC protocols.(6) ### Statistical analyses Descriptive statistics for continuous data are presented as the mean and standard deviation (SD) for normally distributed data, and median and interquartile range (IQR) for non-normal data. Normality of haematology data was assessed by the Shapiro-Wilk test and QQ plots and, if needed, data were normalized by log, square root, or cube root transforms, or reflection followed by transform. Normalized data were standardized to a mean of zero and standard deviation of one. Comparisons for differences between groups was by one-way ANOVA or the Kruskal-Wallis test for non-parametric data. Categorical and binary data were expressed as the number and proportion (%) with comparison using the χ2 test. Univariate logistic regression was used to identify the association of risk factors, or independent variables, such as demography, symptoms, comorbidities, chest x-ray abnormalities, and laboratory assessments, with SARS-CoV-2 infection or death as dependent variables. Multivariate logistic regression with backward selection was used with a cutoff of p<0.15 from univariate regression, and retained variables previously reported to affect the outcomes. Firth’s logistic regression was used to permit calculation of regression coefficients when standard logistic regression failed due to the binary outcome prevalence being very high or very low for a particular risk factor. A two-sided α of <0.05 was considered statistically significant. SPSS V.26.0 and SAS 9.4 were used for all analyses. PCA was performed to identify composite symptom and laboratory components based on their co-variance constructs. Categorical PCA was used for binary symptom data, and conventional PCA for continuous, normalized, and standardized laboratory data. For laboratory data, missing values were imputed as the mean of each variable. A component was retained following cross validation if eigenvalues were above the cutoffs defined by Horn’s parallel analysis. Principal component scores were recoded to tertiles and used as predictors of SARS-CoV-2 infection or death by logistic regression with the lowest tertile as the reference category. HCA was used to define clusters of association of symptoms and other characteristics. The number of retained clusters were determined as 5 based on the highest adjusted Rand index with a high number of clusters and a high pseudo T2 statistic. Cluster memberships were used as predictors of SARS-CoV-2 infection or death by logistic regression. Symptoms and hematological data were analysed as separate variables as risk factors for associations with SARS-CoV-2 infection or death using univariate and multivariate logistic regression, and after classification into patterns or clusters of risk factors defined by PCA or HCA, followed by logistic regression. ## Results In Table 1 we summarize patient characteristics, including age, gender, history of international exposure, history of contact, symptoms, haematology, comorbidities, chest x-ray abnormalities, and diagnosis at discharge. A total of 106 subjects comprised of 18 confirmed COVID-19, 60 suspected, and 28 lung infection cases. The mean age of all subjects was 44.8 ± 18.7 years with 42% being 45-65 years old. The mean age of confirmed cases was 52.6 ±10.4 years, older than suspected patients who were 35.6 ± 17.3 years, but younger than for lung infection cases who were 59.3 ± 14. The gender distribution was nearly equal, except in the confirmed case group wherein 66.7% were male. International travel history was found in 50% of confirmed, and 85.0% of suspected cases. Moreover, 50% of confirmed and 16.7% of suspected cases had contact history with COVID-19 patients. View this table: [Table 1.](http://medrxiv.org/content/early/2020/07/11/2020.07.10.20147942/T1) Table 1. General characteristics of the study patients by category Fever and cough, followed by dyspnoea, were the most frequent symptoms followed by runny nose, sore throat, headache, and gastrointestinal complaints (GIC). Comorbidities were present in nearly half of patients, with combined severe diseases as the most common, followed by pre-existing respiratory disorders and diabetes mellitus. However, there were no pre-existing respiratory disorders in confirmed COVID-19 cases. The two most common disorders of confirmed cases were diabetes (27.8%) and hypertension (27.8%). Chest x-ray showed pneumonia in 94.4% of confirmed and 53.3% of suspected cases. Among the suspected cases, 25% were within normal limits. Consolidation was less frequently reported in confirmed cases. Haematology did not show lymphopenia or thrombocytopenia in confirmed or suspected cases. However, a lower level of leucocytes was seen in suspected and confirmed cases, whereas leucocytosis was observed more frequently in 85% of the lung infection cases. In univariate logistic regression, age ≥ 45 years, positive history of contact, and pneumonia by chest x-ray, were significantly associated with higher risk of SARS-CoV-2 infection. Among all symptoms, runny nose tended toward lower risk for infection. Hypertension and malignancy were comorbidities significantly associated with confirmed cases (Table 2). Multivariate logistic regression found that age ≥ 45 years (aOR 7.81; 95% CI 1.61-38.0), travel history to an affected area (aOR 63.2; 95% CI 4.65-858.9), history of contact with a COVID-19 positive patient (aOR 232; 95% CI 15.1-3569) and pneumonia by chest x-ray (aOR 31.5; 95% CI 1.74-571.3) were independently associated with SARS-CoV-2 infection. View this table: [Table 2.](http://medrxiv.org/content/early/2020/07/11/2020.07.10.20147942/T2) Table 2. Univariate and multivariate logistic regression analyses for predictors of SARS-CoV-2 infection (n=106) Among all patients, including 7 deaths and 99 survivors, univariate analysis (Table 3) showed that age ≥ 45 years, international travel, contact with confirmed COVID-19 case, and diabetes mellitus were associated with death. Dyspnoea tended toward increased risk for death. In the multivariate logistic regression, travel history, presence of history of contact with confirmed COVID-19 cases, and diabetes mellitus retained significant association with mortality. No significant associations were observed between laboratory markers and death. Gender was not associated with SARS-CoV-2 infection or death. View this table: [Table 3.](http://medrxiv.org/content/early/2020/07/11/2020.07.10.20147942/T3) Table 3. Univariate and multivariate regression analyses for predictors of COVID-19 death (n=106) Results for the HCA of symptoms are seen in the dendrogram (Figure 1). HCA cluster 1 (HCA-C1) comprised 7 (38.9%) confirmed, 17 (28.3%) suspected, and 1 (3.6%) lung infection cases. Patients reported predominantly fever with cough, with some also reporting sore throat and GIC, but none with dyspnoea. HCA-C2 comprised 2 (11.1%) confirmed, 11 (18.3%) suspected, and 5 (17.9%) lung infection cases. Patients had predominantly fever, with few combined with dyspnoea, and runny nose. However, cough was not a primary symptom, and none had GIC. Most importantly, patients had fewer symptoms compared to other HCA clusters, and with some reporting no symptoms. HCA-C3 comprised of 17 (28.3%) suspected, but no confirmed or lung infection cases. Patients tended to present with more symptoms, typically 3-5, which were mainly fever, cough, headache, runny nose, and sore throat. A few experienced dyspnoeas and none had GIC. HCA-C4 comprised 1 (5.6%) confirmed, 4 (6.7%) suspected, and 6 (21.4%) lung infection cases. All reported only cough and dyspnoea. HCA-C5 comprised 8 (44.4%) confirmed, 11 (18.3%) suspected and 16 (57%) lung infection cases. All exhibited cough, dyspnoea, and fever. ![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/07/11/2020.07.10.20147942/F1.medium.gif) [Figure 1.](http://medrxiv.org/content/early/2020/07/11/2020.07.10.20147942/F1) Figure 1. Dendrogram of symptoms cluster for confirmed and suspected COVID-19 cases, and non-COVID-19 pneumonia patients The PCA for symptoms revealed two components. PCA symptom 1 (PCA-S1) represented the covariance pattern of patients with headache, runny nose, sore throat, and less dyspnoea with factor loadings of 0.47, 0.43, 0.45, and −0.5, respectively, and accounted for 30% of the total variance across symptoms. PCA-S2 represented a pattern of fever, cough, and GIC with factor loadings of 0.65, 0.51 and 0.63, respectively, and accounted for 18% of variance. The PCA for haematology revealed three components. PCA-H1 represented haemoglobin and haematocrit with loadings of 0.54 and 0.52 and accounted for 42% of total variance. PCA-H2 represented lymphocyte and neutrophil counts with loadings of 0.76 and − 0.83 and accounted for 18% of variance. PCA-H3 had a high factor loading of 1.15 for monocytes and tended toward lower factor loadings of −0.05, −0.05, −0.02, 0.09, −0.12, −0.05, 0.13 and −0.05 for haemoglobin, haematocrit, thrombocytes, leucocytes, lymphocytes, basophils, neutrophils, and eosinophils, respectively, and accounted for 14% of total variance. Univariate logistic regression of symptom data on SARS-CoV-2 infection revealed that runny nose tended toward decreased risk (OR 0.1; 95%CI < 0.001 – 0.83; p = 0.08). No other symptoms exhibited associations in either univariate or multivariate analysis (Table 4). Multivariate analysis of HCA groups revealed that HCA-C3 tended toward reduced risk of COVID-19 cases (aOR 0.048; 95%CI < 0.01 – 0.52; p = 0.056). HCA-C3 comprised patients with symptoms of fever, cough, headache, runny nose, or sore throat, in which only a few experienced dyspnoea, and none had GIC. No other clusters showed any significant associations with COVID-19 cases. However, in all multivariate models with HCA predictors, age ≥ 45 years, presence of international travel history, previous contact with a COVID-19 patient, and pneumonia by x-ray were significant associations, consistent with previous results. For mortality, none of the symptoms nor HCA clusters showed any association in univariate or multivariate analyses (Table 5). However, the presence of diabetes and history of contact with COVID-19 patients were consistently associated with higher mortality in all HCA models. View this table: [Table 4.](http://medrxiv.org/content/early/2020/07/11/2020.07.10.20147942/T4) Table 4. Univariate and multivariate logistic regression analyses of HCA clusters and SARS-CoV-2 infection (n=106) View this table: [Table 5.](http://medrxiv.org/content/early/2020/07/11/2020.07.10.20147942/T5) Table 5. Univariate and multivariate logistic regression analyses of HCA clusters and COVID-19 death (n=106) Analysis of PCA revealed a symptom component for fever, cough, and GIC that tended to increase risk of infection (OR 3.41; 95% CI 1.06-14; p=0.06), and a haematology component with elevated monocytes decreased risk (OR 0.26; 95%CI 0.07-0.79; p=0.027) (Table 6). Multivariate models that included either PCA clusters or the individual predictors were consistent in that presence of travel history and contact with COVID-19 patients were significant predictors of COVID-19 cases. View this table: [Table 6.](http://medrxiv.org/content/early/2020/07/11/2020.07.10.20147942/T6) Table 6. Univariate and multivariate logistic regression analyses of PCA components and SARS-CoV-2 infection (n=106) In univariate analyses for death, HCA-C5 of cough, fever and dyspnoea had increased risk (OR 5.75; P5%CI 1.06 − 31.3, p=0.043). No other individual symptoms or haematology measures, nor HCA clusters or PCA components were associated with death in univariate or multivariate analyses (Table 7). However, we note that diabetes tended toward increased risk for death (aOR 8.86; 95%CI 0.72 – 17.0; p = 0.07), and that history of contact with COVID-19 cases was consistently associated with death (aOR 28.6; 95%CI 1.63->999.9; p = 0.01). View this table: [Table 7.](http://medrxiv.org/content/early/2020/07/11/2020.07.10.20147942/T7) Table 7. Univariate and multivariate logistic regression analyses of PCA components and COVID-19 death (n=106) ## Discussion To our knowledge, this is the first report of confirmed and suspected cases of COVID-19 in Indonesia. It is unique in examining patterns, clusters, or covariance constructs of symptoms and haematological data, or combinations thereof, with other predictors and their associations with key outcomes of infection and death. The present study found consistent results in multivariate analyses for age ≥ 45 years, presence of travel history, contact with a COVID-19 patient, and pneumonia by x-ray as significant predictors of COVID-19 cases. The associations persisted when adjusted with data for symptoms and haematology either individually or as clusters or components. In contrast, analysis of individual symptoms in either univariate or multivariate analysis did not reveal any associations, but HCA yielded a cluster that tended toward reduced risk of COVID-19 cases, and included patients that typically exhibited 3-5 symptoms, which were mainly fever, cough, headache, runny nose, or sore throat, few GIC, and lacking dyspnoea. Other PCA components tended toward associations in univariate analysis, although none was significant, nor were any significant in multivariate analysis. Nevertheless, the stronger tendency for associations of cluster and covariance constructs underscore the importance of interpreting infection as a systemic event. Our results are in line with previous studies that identified older age, history of exposure, chronic underlying disease, symptoms, laboratory findings and chest x-ray examination. Specifically, we observed that persons ≥ 45 years had increased risk of infection, with 75% of confirmed cases being 45-65 years of age. This echoes previous reports wherein mean age of patients was over 40 years, and where critically ill patients were older than non-critical patients.(7) Liu W, et al., showed that age was significantly associated with the COVID-19 cases (OR, 10.6; 95% CI: 2.10– 53.4; p value = 0.004) and disease progression (OR, 8.55; 95% CI: 1.63– 44.8; P = 0.011) (8), similar to our study with aORs ranging from 5.32 – 9.01. We also showed travel history increased the risk of COVID-19 cases with aORs 2.8 − 19.5, and with travel to Singapore or Malaysia being the most common endemic countries recently visited. This is supported by findings from Wells CR, et al. wherein the daily risk of exporting a confirmed case via international travel exceeded 95% by January 13, 2020 (9). They estimated that 64.3% of exported cases were pre-symptomatic while travelling. As expected, we observed that reported history of contact with a COVID-19 case was a very strong predictor of infection (OR 232; 95% CI 15.1-3569), similar to Liang W-hua, et al., where 1,334 of 1,590 (83.9%) had exposure history (10), and Wang R, et al. where 57.6% of patients had a history of exposure to persons from Wuhan within the past 2 weeks.(11) However, the association of recall of exposure with increased mortality was unanticipated, and warrants further exploration. Our observation that pneumonia by chest x-ray strongly associated with infection with aORs 3.15 – 19.8 would be expected, and seen in previous report from Fuyang, China where 80% patients showed bilateral pneumonia by chest x-ray, and all of the critical patients showed pneumonia by x-ray.(11) Similar to other findings, fever, cough and dyspnoea occurred in most patients, with cough and fever associated with infection, although not significant. Univariate regression with cluster analyses confirmed that patients with multiple symptoms, but without dyspnoea, were less likely to be infected (OR 0.04; 95% CI: <0.00 –0.93; P= 0.037). This is consistent with reports that respiratory failure is associated with infection (OR, 8.021; 95% CI: 2.022– 31.821; P = 0.003).(8) The results on GIC warrant further study as angiotensin converting enzyme 2 (ACE-2), a virus receptor, is expressed in gastrointestinal epithelial cells. Viral nucleocapsid protein has also been identified in the gastric, duodenal, and rectum glandular epithelial cell cytoplasm. GI infection by SARS-CoV-2 may be important, and GI symptoms should not be overlooked among suspected cases. We did not observe any significant associations of any laboratory results with infection or death. We note that a high neutrophil-to-lymphocyte ratio (NLR) with cut-off ≥ 3.13 was associated with a non-significant increased risk of infection. This cut-off value was based on previous studies which highlighted the importance of the NLR which reflects systemic inflammation.(7) Other than mortality, we did not explore the association of markers on with clinical severity of illness as evaluated in previous studies. We found diabetes was a consistent and strong predictor for COVID-19 death. Diabetes was previously reported as a risk factor for severe disease and mortality in SARS and MERS (Middle East respiratory syndrome) coronavirus infections, and H1N1 influenza. A recent review summarized that diabetes increased the risk for SARS-CoV-2 infection and hospitalization, with a two-fold risk of severe disease or admission to the intensive care unit (ICU), and increased risk of death.(12) In our study 28% of COVID-19 patients had diabetes or hypertension. Others report increased comorbidity of cardiovascular, endocrine, and digestive disease, of which hypertension (30%) and diabetes (12.1%) were the most common.(13-16) For death, we observed age ≥ 45 years tended toward an association but was not consistent across analyses. A study by Zhou, et al., showed that age is a predictor for mortality. (17) We also report no significant difference in the risk for infection or death based on sex. As mentioned, reported contact history was a predictor for mortality with aORs 51 – 144. One study of 3,019 suspected cases among health workers documented 1,716 confirmed cases with nearly 15% being severe and with 5 deaths. It was noted that transmission occurred mostly among close contacts.(18) The underlying mechanism for how diabetes and older age are associated with adverse outcomes remains unclear, a current study of case series of COVID-19 patients demonstrated the potential contribution of pre-existing endothelial dysfunction. Post-mortem analysis of a renal transplant recipient aged 71 years whose condition deteriorated after COVID-19 diagnosis showed structures of viral inclusion in endothelial cells and accumulation of endothelium-associated inflammatory cells. And post-mortem analysis of a 58-year diabetic woman revealed lymphocytic endotheliosis in lung, heart, kidney, and liver, in addition to evidence of myocardial infarction, but without signs of lymphocytic myocarditis. Endothelial cells express the ACE-2 virus receptor(19), and there may be increased expression of ACE-2 in the presence of comorbidities.(20-22) Hence, COVID-19 could cause widespread endothelial dysfunction from direct infection or immune-mediated recruitment of inflammatory cells, and provoke multiple organ failure. However, the relatively small size of our study precluded meaningful analysis of comorbidities. Indonesia hosts many popular tourist destinations which created logistic, economic and political challenges to limiting foreign visitors. As part of the outbreak response, the Ministry of Health ordered the use of thermal scanners at ports of entry on 5 February 2020, prepared approximately 100 hospitals as COVID-19 referral hospitals, and assigned the National Institute for Health Research and Development to establish a coronavirus diagnostic laboratory supported by WHO and the USA Centers for Disease Control. The first two positive cases in Indonesia were announced by the President on 2 March 2020, but the actual number of cases remains unknown due to lack of active surveillance, lags in testing, and a testing strategy that limits assessment of patients with mild symptoms or persons asymptomatic for respiratory infection. Indonesia now focuses on monitoring and preventing further spread by isolation, self-quarantine, and social distancing especially in large cities, such as Jakarta and Surabaya.(1, 23) The predictors in our study for increased risk of COVID-19 cases (age, international travel history, contact with a COVID-19 patient, pneumonia diagnosis by x-ray) and higher mortality (presence of diabetes and contact history), can be used as essential data for screening and surveillance, and help guide action to reduce transmission and mitigate effects from a possible second wave of COVID-19. There are limitations to our study. First, the sample size was small as this is a preliminary study on COVID-19 during its emergence in Indonesia, and was limited to data that were accessible from that period. We also note the lack of data on smoking and body mass index. This is may be important because ACE-2 expression has been reported to increase in smokers and persons with complications of high body mass index.(24-26) Further investigation of correlates of expression, activation, and interaction of ACE-2 with SARS-CoV-2 in the Indonesian population would be important. The likelihood of underreported cases needs to be considered given the variety of clinical manifestations and suboptimal availability of tools for diagnosis. Actions need to be taken such as a simulation to assess and verify pandemic preparedness for influx of cases at referral hospitals for proper triage, and establishment of high throughput diagnostics laboratories for COVID-19. ## Conclusions In independent and clustered variables analysis, older age, international travel history, contact with a COVID-19 patient, and pneumonia by x-ray were the predictors for increased risk of SARS-CoV-2 infection. The presence of diabetes and history of contact with COVID-19 cases were consistently associated with higher mortality. Association of individual symptoms with infection or mortality were not observed, whereas clustered symptoms or components tended to show stronger associations with infection. Use of such clusters in patient screening and assessment will be useful at the clinic and community level. The Ministry of Health needs to focus on preparedness measures, surveillance and public information to prevent movement of SARS-CoV-2 infection, and create messaging to alleviate public anxiety to enhance action and responsible behaviours. ## Data Availability The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request. ## Declarations ### Ethical approval and consent to participate Ethical approval was obtained from the Research Ethics Committee of the Faculty of Medicine, Universitas Indonesia – Dr. Cipto Mangunkusumo General Hospital, with approval number 20030331. After consideration of ethical issues, logistics and urgency of this work, the Committee waived the requirement for written individual informed consent and approved the sharing of anonymized data. ## Consent for publication Not applicable. ## Availability of data and materials The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request. ## Competing interests The authors declare that they have no competing interests ## Funding This work was supported by grant award from the Directorate Research and Community Services Universitas Indonesia (PUTI Q1) to AJR. ## Authors’ contributions RA, AJR, AFS, DO, EB and AHS designed and coordinated this study. EB, RA, AG, KQP, IGS, P, FI, DP, DL, NGL, VS and ADS collected the data. The data was analysed by RA and interpreted by EB, RA, ISW, and AHS. EB, RA, AHS, ISW, AG, KQP and DR drafted the manuscript together with the help of AJR. AHS provided feedback and helped streamline the content and design analytical procedures. The corresponding author RA coordinated the final submission of the paper. All authors contributed to the final version of the manuscript. ## Acknowledgements We thank Diah Handayani, Andriansjah Rukmana, Muhammad Luqman Labib Zufar, Prasetio Adinugroho, and Annisa Feby Canintika for their insights and outstanding support in drafting this study. We acknowledge the contributions from COVID-19 team of Persahabatan Hospital, Sanglah Hospital, Gatot Subroto Central Army Hospital, and Raden Mattaher Hospital to several authors and contributors to this work. This study was supported by the in kind contributions of the Dean’s Office and the Indonesian Medical Education and Research Institute of the Faculty of Medicine Universitas Indonesia. ## List of abbreviations ACE-2 : angiotensin converting enzyme 2 COVID-19 : Coronavirus disease 2019 CT : computed tomography GIC : gastrointestinal complaints HCA : hierarchical cluster analyses IQR : interquartile range NIHRD : National Institutes for Health Research and Development NLR : neutrophil-to-lymphocyte ratio PCA : principal component analyses RT-PCR : real time polymerase chain reaction SARS-CoV-2 : Severe acute respiratory syndrome coronavirus 2 SD : standard deviation TB : lung tuberculosis * Received July 10, 2020. * Revision received July 10, 2020. * Accepted July 11, 2020. * © 2020, Posted by Cold Spring Harbor Laboratory The copyright holder for this pre-print is the author. All rights reserved. The material may not be redistributed, re-used or adapted without the author's permission. ## References 1. 1.WHO. Update on coronavirus disease in Indonesia 2020 [Available from: <[https://www.who.int/indonesia/news/novel-coronavirus](https://www.who.int/indonesia/news/novel-coronavirus)> 2. 2.WHO. Coronavirus disease (COVID-19) Situation Dashboard 2020 [Available from: [https://covid19.who.int/](https://covid19.who.int/). 3. 3.Liu S LH, Wang Y, Wang D, Ju S, Yang Y. Characteristics and Associations with Severity in COVID-19 Patients: A Multicentre Cohort Study from Jiangsu Province, China. SSRN Electron J. 2020. 4. 4.Indonesia KKR. Pedoman Kesiapsiagaan Menghadapi Coronavirus Disease (COVID-19). 2020. 5. 5.Indonesia KKR. Kesiapsiagaan COVID-19 di Indonesia 2020 [Available from: [http://pusatkrisis.kemkes.go.id/kesiapsiagaan-covid-19-di-indonesia-per-tanggal-15-februari-2020](http://pusatkrisis.kemkes.go.id/kesiapsiagaan-covid-19-di-indonesia-per-tanggal-15-februari-2020). 6. 6.Prevention CfDCa. CDC 2019-Novel Coronavirus (2019-nCoV) Real-Time RT-PCR Diagnostic Panel. 2020. 7. 7.Liu J, Liu Y, Xiang P, Pu L, Xiong H, Li C, et al. Neutrophil-to-Lymphocyte Ratio Predicts Severe Illness Patients with 2019 Novel Coronavirus in the Early Stage 2020. 8. 8.Liu W, Tao ZW, Lei W, Ming-Li Y, Kui L, Ling Z, et al. Analysis of factors associated with disease outcomes in hospitalized patients with 2019 novel coronavirus disease. Chin Med J (Engl). 2020. 9. 9.Wells CR, Sah P, Moghadas SM, Pandey A, Shoukat A, Wang Y, et al. Impact of international travel and border control measures on the global spread of the novel 2019 coronavirus outbreak. Proc Natl Acad Sci U S A. 2020;117(13):7504–9. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMToiMTE3LzEzLzc1MDQiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMC8wNy8xMS8yMDIwLjA3LjEwLjIwMTQ3OTQyLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 10. 10.Liang WH, Guan WJ, Li CC, Li YM, Liang HR, Zhao Y, et al. Clinical characteristics and outcomes of hospitalised patients with COVID-19 treated in Hubei (epicenter) and outside Hubei (non-epicenter): A Nationwide Analysis of China. Eur Respir J. 2020. 11. 11.Wang R, Pan M, Zhang X, Fan X, Han M, Zhao F, et al. Epidemiological and clinical features of 125 Hospitalized Patients with COVID-19 in Fuyang, Anhui, China. Int J Infect Dis. 2020. 12. 12.Madsbad S. COVID-19 INFECTION IN PEOPLE WITH DIABETES 2020. Available from: <[https://www.touchendocrinology.com/insight/covid-19-infection-in-people-with-diabetes/](https://www.touchendocrinology.com/insight/covid-19-infection-in-people-with-diabetes/)>. 13. 13.Zhang JJ, Dong X, Cao YY, Yuan YD, Yang YB, Yan YQ, et al. Clinical characteristics of 140 patients infected with SARS-CoV-2 in Wuhan, China. Allergy. 2020. 14. 14.Guan WJ, Ni ZY, Hu Y, Liang WH, Ou CQ, He JX, et al. Clinical Characteristics of Coronavirus Disease 2019 in China. N Engl J Med. 2020. 15. 15.Kim ES, Chin BS, Kang CK, Kim NJ, Kang YM, Choi JP, et al. Clinical Course and Outcomes of Patients with Severe Acute Respiratory Syndrome Coronavirus 2 Infection: a Preliminary Report of the First 28 Patients from the Korean Cohort Study on COVID-19. J Korean Med Sci. 2020;35(13):e142. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3346/jkms.2020.35.e142&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F11%2F2020.07.10.20147942.atom) 16. 16.Hussain A, Bhowmik B, do Vale Moreira NC. COVID-19 and diabetes: Knowledge in progress. Diabetes Res Clin Pract. 2020;162:108142. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.diabres.2020.108142&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32278764&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F11%2F2020.07.10.20147942.atom) 17. 17.Zhou F, Yu T, Du R, Fan G, Liu Y, Liu Z, et al. Clinical course and risk factors for mortality of adult inpatients with COVID-19 in Wuhan, China: a retrospective cohort study. The Lancet. 2020;395(10229):1054–62. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/SO140-6736(20)30566-3&link_type=DOI) 18. 18. M ZWJ, McGoogan. Characteristics of and Important Lessons From the Coronavirus Disease 2019 (COVID-19) Outbreak in China Summary of a Report of 72.314 Cases From the Chinese Center for Disease Control and Prevention. 2020. 19. 19.Lau SK, Woo PC, Yip CC, Tse H, Tsoi HW, Cheng VC, et al. Coronavirus HKU1 and other coronavirus infections in Hong Kong. J Clin Microbiol. 2006;44(6):2063–71. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamNtIjtzOjU6InJlc2lkIjtzOjk6IjQ0LzYvMjA2MyI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzA3LzExLzIwMjAuMDcuMTAuMjAxNDc5NDIuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 20. 20.Oudit GY, Imai Y, Kuba K, Scholey JW, Penninger JM. The role of ACE2 in pulmonary diseases--relevance for the nephrologist. Nephrol Dial Transplant. 2009;24(5):1362–5. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/ndt/gfp065&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19228756&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F11%2F2020.07.10.20147942.atom) 21. 21.Riviere G, Michaud A, Breton C, VanCamp G, Laborie C, Enache M, et al. Angiotensin-converting enzyme 2 (ACE2) and ACE activities display tissue-specific sensitivity to undernutrition-programmed hypertension in the adult rat. Hypertension. 2005;46(5):1169–74. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1161/01.HYP.0000185148.27901.fe&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16203874&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F11%2F2020.07.10.20147942.atom) 22. 22.Roca-Ho H, Riera M, Palau V, Pascual J, Soler MJ. Characterization of ACE and ACE2 Expression within Different Organs of the NOD Mouse. Int J Mol Sci. 2017;18(3). 23. 23.Jakarta GD. Instruksi Gubernur No. 23 Tahun 2020. 2020. 24. 24.Cai G. 2020. Tobacco-use disparity in gene expression of ACE2, the receptor of 2019-nCov. 25. 25.Wu C ZM. Single-cell RNA expression profiling shows that ACE2, the putative receptor of Wuhan has significant expression in the nasal, mouth,lung and colon tissues, and tends to be co-expressed with HLA-DRB1 in the four tissues. 2020. 26. 26.Zhao Y ZZ, Wang Y, Zhou Y, Ma Y, Zuo W. Single-cell RNA expression profiling of ACE2, the putative receptor of Wuhan 2019-nCov. bioRxiv. 2020.