Article Text

Download PDFPDF

Original research
History of coronary heart disease increased the mortality rate of patients with COVID-19: a nested case–control study
  1. Tian Gu1,
  2. Qiao Chu2,
  3. Zhangsheng Yu3,
  4. Botao Fa3,
  5. Anqi Li4,
  6. Lei Xu5,
  7. Ruijun Wu4,
  8. Yaping He2,6
  1. 1Department of Biostatistics, University of Michigan, Ann Arbor, Michigan, USA
  2. 2School of Public Health, Shanghai Jiao Tong University School of Medicine, Shanghai, China
  3. 3Department of Bioinformatics and Biostatistics, Shanghai Jiao Tong University, Shanghai, China
  4. 4School of Social Development, East China Normal University, Shanghai, China
  5. 5Department of Cardiology, Shanghai Chest Hospital, Shanghai, China
  6. 6Center for Health Technology Assessment, Shanghai Jiao Tong University, Shanghai, China
  1. Correspondence to Dr. Yaping He; hypcyr{at}sina.com; Dr. Ruijun Wu; rjwu{at}re.ecnu.edu.cn

Abstract

Objective Evaluate the risk of pre-existing comorbidities on COVID-19 mortality, and provide clinical suggestions accordingly.

Setting A nested case–control design using confirmed case reports released from the news or the national/provincial/municipal health commissions of China between 18 December 2019 and 8 March 2020.

Participants Patients with confirmed SARS-CoV-2 infection, excluding asymptomatic patients, in mainland China outside of Hubei Province.

Outcome measures Patient demographics, survival time and status, and history of comorbidities.

Method A total of 94 publicly reported deaths in locations outside of Hubei Province, mainland China, were included as cases. Each case was matched with up to three controls, based on gender and age ±1 year old (94 cases and 181 controls). The inverse probability-weighted Cox proportional hazard model was performed, controlling for age, gender and the early period of the outbreak.

Results Of the 94 cases, the median age was 72.5 years old (IQR=16), and 59.6% were men, while in the control group the median age was 67 years old (IQR=22), and 64.6% were men. Adjusting for age, gender and the early period of the outbreak, poor health conditions were associated with a higher risk of COVID-19 mortality (HR of comorbidity score, 1.31 [95% CI 1.11 to 1.54]; p=0.001). The estimated mortality risk in patients with pre-existing coronary heart disease (CHD) was three times that of those without CHD (p<0.001). The estimated 30-day survival probability for a profile patient with pre-existing CHD (65-year-old woman with no other comorbidities) was 0.53 (95% CI 0.34 to 0.82), while it was 0.85 (95% CI 0.79 to 0.91) for those without CHD. Older age was also associated with increased mortality risk: every 1-year increase in age was associated with a 4% increased risk of mortality (p<0.001).

Conclusion Extra care and early medical interventions are needed for patients with pre-existing comorbidities, especially CHD.

  • epidemiology
  • coronary heart disease
  • survival analysis
http://creativecommons.org/licenses/by-nc/4.0/

This is an open access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited, appropriate credit is given, any changes made indicated, and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/.

Statistics from Altmetric.com

Request Permissions

If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.

Strengths and limitations of this study

  • Since we used data outside of the epicentre of the outbreak in mainland China, we avoided the competing mortality risk caused by insufficient healthcare resources.

  • The study controlled for the confounding effects of age, gender and the early period of the pandemic, which were, in general, not considered in the existing literatures.

  • The survival-time-related statistical results can help guide the early clinical intervention to reduce mortality.

  • The lack of medical test results in the publicly reported data restricted further investigation on the association between comorbidity-related clinical index and mortality.

  • The missing at random assumption could not be verified.

Introduction

Since the first report of COVID-19 in December 2019 in Wuhan, Hubei Province, China, the novel virus infection has rapidly spread to other cities in China and has now been detected in 215 countries and locations internationally.1 On 11 March 2020, the WHO declared COVID-19 a pandemic and has called for aggressive actions from all countries to fight the disease. Current research has indicated that COVID-19 is caused by SARS-CoV-2, a beta-coronavirus similar to the SARS-CoV in its genetic sequence.2 Epidemiological evidence suggests that the initial reported cases in China had a history of exposure to the Huanan seafood market.3 4 With the escalated spread of the infection, there has been clear evidence of human-to-human transmission.5 6 The most common symptoms include fever, dry cough and fatigue,5–8 and there are asymptomatic yet contagious cases.9

According to the COVID-19 situation reports of World Health Organization (WHO)1, as of 31 July 2020, there were 87 956 confirmed cases in mainland China, including 4666 deaths. Internationally, a total of 17 million confirmed cases and 664 244 deaths have been reported outside of China. Considering the global public health threat posed by COVID-19, unravelling the prognostic factors for patients, especially the risk factors of mortality associated with COVID-19, has important implications for clinical practice and is urgently warranted.

Studies have indicated that in severe cases, patients tend to be older in age6 8 and are more likely to have had pre-existing medical conditions, including but not limited to hypertension,3 6 8 diabetes,6 8 cardiovascular diseases,3 6 8 10–12 cerebrovascular diseases,6 chronic obstructive pulmonary disease (COPD),3 8 cancer13 and digestive diseases,14 in comparison to those with non-severe cases.3 6 8–15

Recently, scholarly attention has focused on identifying the risk factors for death from COVID-19. Some evidence suggests that pre-existing medical conditions are likely risk factors for death from COVID-19. For example, a study based on 72 314 cases in China indicated that the case-fatality rate (CFR) tends to be higher among those who are older in age and who have pre-existing cardiovascular disease, diabetes and/or hypertension.9 Similarly, by conducting logistic regression on odds of in-hospital deaths among 54 diseased patients and 137 recovered patients in Wuhan City, Hubei Province, Zhou et al16 found that older age, higher Sequential Organ Failure Assessment score and D-dimer greater than 1 μg/mL at hospital admission were associated with increased odds of in-hospital death. Chen et al4 found that pre-existing hypertension and other cardiovascular complications were more common among diseased patients than recovered patients.

However, several gaps remain in the understanding of risk factors for mortality of COVID-19. First, most current research on pre-existing comorbidities of COVID-19 was based on univariate comparison, which did not account for important confounders such as age and gender.17–20 Second, no studies have investigated the hazard of the identified risk factors over time or the probability of survival at a given time. Under the rapidly changing pandemic situation, it is crucial to provide timely survival-time guidance for implementing the targeted treatment to the high-risk patients in clinical practice. Third, most existing studies on mortality risk factors were focused on patients diagnosed in Wuhan, Hubei Province, with little understanding about the mortality risk factors outside of Hubei Province. The risk factors are likely different inside and outside of Hubei Province since current research has found the clinical symptom severity5 and the CFR9 21 to be higher in Hubei Province (the centre of outbreak) than in cities outside of Hubei Province in China. Fourth, no studies thus far have taken into account the pandemic stage when evaluating mortality risk factors. It has been found that the average daily infection rate in China was different before and after 22 January 2020, since non-pharmaceutical interventions were taken by the government before this date.22 The change of time in the pandemic stage may also influence the risk factors for fatality associated with COVID-19.

To fill the above research gaps in the existing literature about the mortality risk factors for COVID-19, the present study was conducted using a nested case–control (NCC) study, which aimed to evaluate the risk of the common pre-existing comorbidities (hypertension, CHD, diabetes, etc) for mortality associated with COVID-19 in mainland China outside of Hubei Province. NCC, also called risk set sampling, has been widely used in studying fatal disease risk effect in large pharmacoepidemiological studies23–28 and risk prediction in pandemic influenza A (H1N1) 2009 (pH1N1).29 NCC is cost-effective in data collection and is especially suitable for research on the mortality risk of diseases such as COVID-19, where the number of event-free people largely exceeds those who are symptomatic.30 To attain this goal, we employed survival analysis on 275 publicly reported confirmed cases, adjusting for age, gender and the change in the pandemic stage in China (ie, before and after 22 January 2020).

Method

Study design and rationale

This study performed survival analysis under an NCC design to assess the roles of common comorbidities (cardiocerebrovascular, endocrine and respiratory disease, etc) in predicting mortality for COVID-19 among patients in mainland China outside of Hubei Province. The study period was from 18 December 2019, when the first laboratory-confirmed case was announced in China, to 8 March 2020.

The study cohort was defined as all the publicly reported confirmed patients with COVID-19 outside of Hubei Province in mainland China during the study period. During this period, 112 deaths outside of Hubei Province were reported by the National Health Committee of China, and 18 were excluded from the present study due to unavailablility of case reports. A total of 448 publicly reported laboratory-confirmed COVID-19 cases (94 deaths and 354 survivors) were initially collected. To avoid selection bias due to intentionally collecting patients with certain pre-existing comorbidities, two authors independently collected, compared and reviewed the full text of each case report. Following the typical NCC design setting where all deaths were included as cases, and each case was matched with up to three controls on gender and age ±1 year old (94 cases and 181 controls). The sample distribution across all 32 province-level regions in mainland China is presented in online supplemental table S1.

Data collection procedure

We routinely searched for daily news and public health reports on confirmed COVID-19 cases in all areas in mainland China outside of Hubei Province. Patients’ clinical and comorbidity characteristics were recorded and doubly confirmed by national/provincial/municipal health commission websites, the official COVID-19 data reporting websites in China. Follow-up time was defined as the duration from the date of disease onset until the end of observation on 8 March 2020, or when the participant died, whichever came first. For each eligible patient, we followed local reports to update their survival status until the end of follow-up time. Pre-existing comorbidities were recorded based on the descriptions in case reports.

As illustrated in figure 1, the inclusion criteria were publicly reported patients with COVID-19 who had complete information on basic demographics (age, gender and region), disease onset date (the first time they became symptomatic) and the history of comorbidities (including hypertension, CHD, cardiac failure, cerebral infarction, diabetes, chronic bronchitis, COPD, renal failure and history of surgery). Asymptomatic patients were not included in this study. In addition, we defined ‘comorbidity-free patients’ as those who were specifically described as ‘no pre-existing medical condition/comorbidity’ on the national/provincial/municipal health commission websites.

Figure 1

Patient flow diagram detailing included subjects and exclusion criteria.

In the following three steps, we used patient no. 213 in the sample as an example to introduce the dynamic tracking method we used to identify any missing dates.

Step 1. We conducted an internet search on confirmed cases on baidu.com, the largest search engine in China, using the keywords “confirmed COVID-19 cases report” and “pre-existing comorbidities”. A search result pertained to one confirmed case reported on the website of the Municipal Health Commission of Binzhou (Shandong Province) on 17 February, described as ‘the 15th confirmed case: 30-year-old man without pre-existing morbidities, who lives in the neighbourhood of Xincun Village. This patient was diagnosed positive on 16 February and is being treated with precaution in Bincheng hospital’. We recorded the age, gender and region of this patient and that he was comorbidity-free.

Step 2. We then determined the onset date of COVID-19 for this patient based on another announcement on the same website. In the announcement titled ‘Possible exposure locations and times of the 15th confirmed case’, it says, ‘the patient was symptomatic on 14 February’.

Step 3. Finally, following the updates on the website, we confirmed the event status of this patient as discharged on 3 March 2020.

Statistical analysis

Analyses were performed in R V.3.6.2 (R Foundation for Statistical Computing) through RStudio V.1.2.5042. Data and code are available online at GitHub (https://github.com/GuTian-TianGu/COVID-19_NCCstudy). Baseline clinical characteristics were shown as mean (SD), median (range) or number (%), with a comparison of characteristics in subjects stratified by case and control via the non-parametric Mann-Whitney U test for continuous variables and χ2 or Fisher’s exact test for binary variables.

In order to use the time-to-event information under the NCC design, the inverse probability-weighted Cox proportional hazard regression model was employed.31 The matching between cases and controls, as well as the relative weights were simultaneously obtained via KMprob function in multiple NCC R package32 by specifying the Kaplan-Meier type weights with additional matching on gender and age ±1 year old. Only survivors were assigned weights since all cases (deaths) were included as designed with a weight of one. A total of 113 survivors (mean age 46.5) with sampling probabilities of zero were considered ‘fail to match’ and excluded from the study, mainly due to younger age than cases. A majority of the excluded patients were from Shandong Province (38.1%) due to the relatively high representation of the sample (detailed information of excluded survivors is available online at GitHub). In a sensitivity analysis adjusting for Shandong Province (results not shown here), we observed the consistent results as the main analysis.

The comorbidity score, ranging from zero to nine, was defined as the summation of nine comorbidities that have been specifically mentioned in relation to COVID-19 outcomes (CHD, hypertension, cardiac failure, cerebral infarction, chronic bronchitis, COPD, diabetes, renal failure and history of surgery). The Kaplan-Meier curve was plotted to check the proportional hazard assumption, and the Pearson correlation test was used to rule out the multicollinearity concern before fitting any model. Univariate weighted Cox models were performed for each comorbidity. The multivariate weighted Cox model was used to determine if pre-existing comorbidity yielded prognostic hazard information. We included those comorbidities that were marginally significant (p<0.1) in the univariate analysis to the multivariate model. Other than the common risk factors (age and gender), the multivariate model also adjusted for the early period of the pandemic (after vs before 22 January 2020, when no intervention was taken by the government).22 Although matching was based on age and gender, we adjusted for the matching covariates since the matching was broken with inverse probability weighting.31 A separate multivariate model was built by using the comorbidity score as a continuous predictor, adjusting for the same covariates. Hazard ratios (HRs) from the weighted Cox model were reported along with 95% CIs and p values. Sensitivity analysis was performed using multivariate logistic regression to provide estimated odds ratio (ORs), which included the same covariates as the multivariate weighted Cox model.

Weighted Cox model-based survival estimates were plotted for a sample patient profile (65-year-old woman with no other comorbidities) to compare the survival probability over time with and without CHD. The log-rank test was used to compare the median survival difference.

Results

Sample description

Table 1 summarises patient demographics and pre-existing medical conditions. Results are presented for all patients in the study (n=275) as well as for cases (n=94) and controls (n=181), respectively. Patients were 24 to 94 years old (Meanage=66.4, SDage=14.5). The average age tended to be older in the case group (70.7 years old) than in the control group (64.2 years old). Median ages were similar to mean ages in both groups. A majority (62.9%) of the patients were men. Overall, 25.5% of the total patients had clinical symptoms associated with COVID-19 before 22 January 2020. A relatively small proportion of the total sample had COPD, renal failure, history of surgery and hepatic failure (4.4%, 4.4%, 3.6% and 1.1%, respectively). Among all pre-existing comorbidities with over 5% of the total sample, hypertension was the most common (39.6%), followed by diabetes (26.2%), CHD (14.5%), cardiac failure (8%), cerebral infarction (6.9%) and chronic bronchitis (6.9%). Patients in the case group had more CHD (p<0.001) and more cerebral infarction (p=0.05) than those in the control group. The mean comorbidity score was 1.22 (SD=1.21) in the overall sample and 1.6 (SD=1.32) in the case and 1.02 (SD=1.10) in the control group (p<0.001).

Table 1

Patientcharacteristics, stratified by survival status*

Model results

Results of the Pearson correlation test showed no significant correlations with the presence of comorbidities of interest, and the assumption of the proportional hazard was not violated.

Table 2 presents the results of univariate and multivariate weighted Cox models. Older age was associated with significantly higher risk of death with similar magnitude in univariate and multivariate models. In the adjusted model, every 1-year increase in age was associated with an estimated 4% higher risk of death (p<0.001). No significant hazard difference was found between male and female patients. Disease infection during the early no-intervention period was associated with a higher risk of death but was not statistically significant.

Table 2

Univariate and multivariate model result from weighted Cox proportional hazard regression

In a separate model using the comorbidity score as a predictor, we observed that a higher comorbidity score was associated with a higher mortality risk in both unadjusted and adjusted models (p=0.009 and p=0.03, respectively). All pre-existing comorbidities had an HR over 1 in the univariate model, of which CHD had the largest HR of 4.2 (p<0.001), followed by cerebral infarction (HR=2.9, p=0.004), COPD (HR=2.6, p=0.01), renal failure (HR=2.3, p=0.09), cardiac failure (HR=1.9, p=0.1), history of surgery (HR=1.7, p=0.34), hypertension (HR=1.4, p=0.17), diabetes (HR=1.1, p=0.61) and chronic bronchitis (HR=1.1, p=0.55), but not all were statistically significant. After adjusting for age, gender and the early period of the pandemic in China, CHD was the only comorbidity that yielded a significant mortality risk: patients who had COVID-19 with pre-existing CHD had an estimated 2.9 times higher risk of death than those without CHD. In addition, cerebral infarction, COPD and renal failure all had an estimated HR of around 2.0. Similar results were observed by using unweighted logistic regression in a sensitivity analysis (online supplemental table S2).

The overall median follow-up was 40 days, during which 94 deaths were observed. Figure 2 shows the estimated survival probability over 70 days for a sample patient profile with and without CHD (65-year-old woman with no other comorbidities). For such a patient profile, having pre-existing CHD led to a significantly shorter survival probability over time compared with those without CHD (p<0.001). For those with CHD, the estimated 30-day survival probability was 0.53 (95% CI 0.34 to 0.82); on the other hand, this number was 0.85 (95% CI 0.79 to 0.91) for those without CHD.

Figure 2

Estimated survival probability over time from the adjusted Cox proportional hazard model for an example patient profile (65-year-old woman without CHD (dashed line) or with CHD (solid line), who had no other comorbidities). The estimated 30-day survival probability was 0.53 (95% CI 0.34 to 0.82) for patients with pre-existing CHD, while 0.85 (95% CI 0.79 to 0.91) for those without (p<0.001). CHD, coronary heart disease.

Discussion

In our research, based on publicly reported confirmed cases and adjusted for the confounding effect of age, gender and the early period of the pandemic when no intervention was taken, we used survival analysis to estimate the fatality risk of pre-existing comorbidities in COVID-19. There were three major findings: first, poor health condition was associated with higher mortality risk of COVID-19. Second, after adjusting for confounders, CHD was the only significant risk factor for COVID-19 mortality. Patients with pre-existing CHD were 3.11 times more likely to die than those without CHD (p<0.001). For one patient profile (65-year-old woman without other comorbidities), we saw an estimated 30-day survival probability of 0.85 (p<0.001). However, for those with CHD, the 30-day survival probability was 0.53. Third, older age was associated with an increased risk of death. Specifically, every 1-year increase in age was associated with a 4% increased risk of mortality (p<0.001).

To the best of our knowledge, the present study is the first to provide substantial statistical evidence showing the effect of CHD in predicting mortality for COVID-19. This result is consistent with previous studies that found higher CFR among patients with cardiovascular disease.9 15 It is worth pointing out that the existing studies9 16 that investigate prognostic factors for death used χ2 tests or univariate logistic regression that does not control for potential confounders. In contrast, by conducting weighted Cox proportional hazard regression models, our study used time-to-event outcomes, which offer more survival-time-related information to help guide the clinical intervention and more statistical power to detect risk factors.33 34

Previous studies have indicated that cardiovascular events following pneumonia may increase the risk of mortality,4 35–40 which explained our findings from the viewpoint of pathophysiological mechanisms. One potential mechanism underlying the association between pneumonia and cardiovascular events is inflammation.37 Specifically, the inflammatory reaction following pneumonia can result in plaque instability and damage in the blood vessels. Evidence of elevated local inflammation in the atherosclerotic coronary arteries following acute systemic infections has been shown in many studies.37 39 Thus, infections may result in heightened loading that is imposed on cardiomyocytes and lead to sympathetic hyperactivity or ischaemia, which may increase the risk of arrhythmia and heart failure in patients who had COVID-19 with pre-existing CHD.36

Given the limited understanding of the prognostic factors for COVID-19, more research, including potentially prospective studies, is needed to investigate the mechanism by which pre-existing CHD may influence the survival probability of patients with COVID-19. From the clinical point of view, early evaluation of a patient’s medical history is necessary to implement early medical interventions and decrease the mortality risk. We suggest monitoring the dynamic heart rate for patients with pre-existing CHD. For those severe symptomatic patients who had pre-existing heart ischaemia and abnormal heart function, early medical intervention may be needed.40

Furthermore, our results indicated that the risk of death from COVID-19 was significantly higher in older patients. Adjusting for others, every 1-year increase in age was associated with a 4% increased risk of death, which is similar to what was found in previous studies.8 9 16 Alongside the evidence of prognostic risk in CHD, we suggest that extra care is needed for those with CHD, especially for elderly patients.

Previous studies yielded mixed results regarding gender differences in mortality risk of COVID-19. Some studies found male sex was associated with a higher risk of death from COVID-19,41 whereas other studies did not find gender to be a significant factor predicting the mortality risk of COVID-19.10 16 In the current study, gender was not a significant mortality risk factor for COVID-19. More research is needed to further our understanding of gender differences in the outcome of COVID-19 and underlying mechanisms.

The design of excluding patients from Hubei Province was based on the concerns of unknown confounders caused by insufficient medical resources in the epicentre. Centers for Disease Control and Precention (CDC)’s report has pointed out that the rapidly increasing number of infections could easily crash the healthcare system by exceeding its maximum capacity.42 Therefore, analysing patient data outside of Hubei Province avoided the competing mortality risk caused by insufficient healthcare resources and revealed the true underlying impact of pre-existing comorbidities on COVID-19 mortality.43

This study has several limitations. It is worth noting that the data were collected when COVID-19 was spreading rapidly in China and the health authorities and researchers had limited understanding of the incubation period, modes of viral transmission and effective treatment. Whether our findings can be generalised to later epidemic phases warrants future research. One limitation of the present study lies in the nature of publicly reported data. Researchers have pointed out that severe cases may be over-represented in publicly reported data.44 Nevertheless, we have managed to reduce the potential bias caused by severe case over-representation by appropriately matching the cases and controls in the NCC design. Following the NCC design, we allowed the controls to be matched with cases on age ±1 year old instead of the exact matching, which caused the significant age difference between two groups (table 1). Thus, we adjusted for the matching covariates in all the models to address this.31 The auto-matching procedure via the statistical programme also prevented the possibility of tendentiously selecting survivors with comorbidity-free history during data collection. In addition, NCC design is favoured in our situation where the risk factor data and event of interest can be identified opportunistically from publicly reported confirmed cases.30 Therefore, NCC was the optimal choice, given the restricted availability of public data. Moreover, due to the lack of information of treatment in the health reports published by the local health commission websites, we did not include treatment information into analysis, which may produce confounding effects. It calls for future research to investigate the mortality risk effect of pre-existing comorbidities by adding treatment as a covariate in the model. Finally, we were not able to verify the missing at random assumption of the 18 missing deaths (four from Heilongjiang, five from Henan, five from Beijing, two from Hunan, one from Liaoning and one from Yunnan) as well as the survivors.

In conclusion, our findings provided preliminary yet strong evidence supporting the association between pre-existing CHD and mortality risk for patients with COVID-19. Based on our findings, close monitoring, extra care and early medical intervention are needed for patients with pre-existing CHD to reduce the mortality risk associated with COVID-19.

References

Footnotes

  • Twitter @GuTian_TianGu

  • Contributors TG conducted the analysis, interpreted the data, drafted the article and rechecked the transcribed manuscript. QC conducted the literature review and drafted the article. ZY guided and supervised the statistical analysis. BF helped with data analysis. AL helped with literature review, data collection, management and the data quality check, and was responsible for reference organization. LX provided clinical support and interpretation of the results. RW collected, preprocessed and managed data, and conducted preliminary data analysis. YH guided and supervised the research process. RW and YH provided funding support. All authors critically revised the manuscript for intellectual content, approved the final draft and agreed to accountability for all aspects of the work.

  • Funding This project was funded by the National Natural Science Foundation of China (No. 71874111), China Ministry of Education Key Research Institute of Humanities and Social Sciences at Universities (No. 17JJD790008), Shanghai Municipal Health Bureau Foundation (No. 201740116), and Shanghai Jiao Tong University Scientific and Technological Innovation Funds (YG2020YQ01, YG2020YQ06).

  • Competing interests None declared.

  • Patient and public involvement Patients and/or the public were not involved in the design, conduct, reporting, or dissemination plans of this research.

  • Patient consent for publication Exempt.

  • Ethics approval The data published in news reports and websites were open to the public and free of identifiers.The study was approved by Shanghai Jiao Tong University Public Health and Nursing Medical Research Ethics Committee (SJUPN-202001).

  • Provenance and peer review Not commissioned; externally peer reviewed.

  • Data availability statement Data are available in a public, open access repository. Raw data are available at https://github.com/GuTian-TianGu/COVID-19_NCCstudy.git.