Urinary microbiome-based metagenomic signature for the noninvasive diagnosis of hepatocellular carcinoma.

BACKGROUND: Gut microbial dysbiosis is implicated in chronic liver disease and hepatocellular carcinoma (HCC), but the role of microbiomes from various body sites remains unexplored. We assessed disease-specific alterations in the urinary microbiome in HCC patients, investigating their potential as diagnostic biomarkers. METHODS: We performed cross-sectional analyses of urine samples from 471 HCC patients and 397 healthy controls and validated the results in an independent cohort of 164 HCC patients and 164 healthy controls. Urinary microbiomes were analyzed by 16S rRNA gene sequencing. A microbial marker-based model distinguishing HCC from controls was built based on logistic regression, and its performance was tested. RESULTS: Microbial diversity was significantly reduced in the HCC patients compared with the controls. There were significant differences in the abundances of various bacteria correlated with HCC, thus defining a urinary microbiome-derived signature of HCC. We developed nine HCC-associated genera-based models with robust diagnostic accuracy (area under the curve [AUC], 0.89; balanced accuracy, 81.2%). In the validation, this model detected HCC with an AUC of 0.94 and an accuracy of 88.4%. CONCLUSIONS: The urinary microbiome might be a potential biomarker for the detection of HCC. Further clinical testing and validation of these results are needed in prospective studies.

Features of T1 pancreatic cancer and validation of the eighth edition AJCC staging system definition using a Korean-Japanese joint cohort and the SEER database.

BACKGROUND/PURPOSE: Little is known about the features of T1 pancreatic ductal adenocarcinoma (PDAC) and its definition in the eighth edition of the American Joint Committee on Cancer (AJCC) staging system needs validation. The aims were to analyze the clinicopathologic features of T1 PDAC and investigate the validity of its definition. METHOD: Data from 1506 patients with confirmed T1 PDAC between 2000 and 2019 were collected and analyzed. The results were validated using 3092 T1 PDAC patients from the Surveillance, Epidemiology, and End Results (SEER) database. RESULTS: The median survival duration of patients was 50 months, and the 5-year survival rate was 45.1%. R0 resection was unachievable in 10.0% of patients, the nodal metastasis rate was 40.0%, and recurrence occurred in 55.2%. The current T1 subcategorization was not feasible for PDAC, tumors with extrapancreatic extension (72.8%) had worse outcomes than those without extrapancreatic extension (median survival 107 vs. 39 months, p < .001). Extrapancreatic extension was an independent prognostic factor whereas the current T1 subcategorization was not. The results of this study were reproducible with data from the SEER database. CONCLUSION: Despite its small size, T1 PDAC displayed aggressive behavior warranting active local and systemic treatment. The subcategorization by the eighth edition of the AJCC staging system was not adequate for PDAC, and better subcategorization methods need to be explored. In addition, the role of extrapancreatic extension in the staging system should be reconsidered.

Proteogenomic landscape of human pancreatic ductal adenocarcinoma in an Asian population reveals tumor cell-enriched and immune-rich subtypes.

We report a proteogenomic analysis of pancreatic ductal adenocarcinoma (PDAC). Mutation-phosphorylation correlations identified signaling pathways associated with somatic mutations in significantly mutated genes. Messenger RNA-protein abundance correlations revealed potential prognostic biomarkers correlated with patient survival. Integrated clustering of mRNA, protein and phosphorylation data identified six PDAC subtypes. Cellular pathways represented by mRNA and protein signatures, defining the subtypes and compositions of cell types in the subtypes, characterized them as classical progenitor (TS1), squamous (TS2-4), immunogenic progenitor (IS1) and exocrine-like (IS2) subtypes. Compared with the mRNA data, protein and phosphorylation data further classified the squamous subtypes into activated stroma-enriched (TS2), invasive (TS3) and invasive-proliferative (TS4) squamous subtypes. Orthotopic mouse PDAC models revealed a higher number of pro-tumorigenic immune cells in TS4, inhibiting T cell proliferation. Our proteogenomic analysis provides significantly mutated genes/biomarkers, cellular pathways and cell types as potential therapeutic targets to improve stratification of patients with PDAC.

Multi-biomarker panel prediction model for diagnosis of pancreatic cancer.

BACKGROUND/PURPOSE: The current study aimed to develop a prediction model using a multi-marker panel as a diagnostic screening tool for pancreatic ductal adenocarcinoma. METHODS: Multi-center cohort of 1991 blood samples were collected from January 2011 to September 2019, of which 609 were normal, 145 were other cancer (colorectal, thyroid, and breast cancer), 314 were pancreatic benign disease, and 923 were pancreatic ductal adenocarcinoma. The automated multi-biomarker Enzyme-Linked Immunosorbent Assay kit was developed using three potential biomarkers: LRG1, TTR, and CA 19-9. Using a logistic regression model on a training data set, the predicted values for pancreatic ductal adenocarcinoma were obtained, and the result was classification into one of the three risk groups: low, intermediate, and high. The five covariates used to create the model were sex, age, and three biomarkers. RESULTS: Participants were categorized into four groups as normal (n = 609), other cancer (n = 145), pancreatic benign disease (n = 314), and pancreatic ductal adenocarcinoma (n = 923). The normal, other cancer, and pancreatic benign disease groups were clubbed into the non-pancreatic ductal adenocarcinoma group (n = 1068). The positive and negative predictive value, sensitivity, and specificity were 94.12, 90.40, 93.81, and 90.86, respectively. CONCLUSIONS: This study demonstrates a significant diagnostic performance of the multi-marker panel in distinguishing pancreatic ductal adenocarcinoma from normal and benign pancreatic disease states, as well as patients with other cancers.

Combining p-values from various statistical methods for microbiome data.

Motivation: In the field of microbiome analysis, there exist various statistical methods that have been developed for identifying differentially expressed features, that account for the overdispersion and the high sparsity of microbiome data. However, due to the differences in statistical models or test formulations, it is quite often to have inconsistent significance results across statistical methods, that makes it difficult to determine the importance of microbiome taxa. Thus, it is practically important to have the integration of the result from all statistical methods to determine the importance of microbiome taxa. A standard meta-analysis is a powerful tool for integrative analysis and it provides a summary measure by combining p-values from various statistical methods. While there are many meta-analyses available, it is not easy to choose the best meta-analysis that is the most suitable for microbiome data. Results: In this study, we investigated which meta-analysis method most adequately represents the importance of microbiome taxa. We considered Fisher's method, minimum value of p method, Simes method, Stouffer's method, Kost method, and Cauchy combination test. Through simulation studies, we showed that Cauchy combination test provides the best combined value of p in the sense that it performed the best among the examined methods while controlling the type 1 error rates. Furthermore, it produced high rank similarity with the true ranks. Through the real data application of colorectal cancer microbiome data, we demonstrated that the most highly ranked microbiome taxa by Cauchy combination test have been reported to be associated with colorectal cancer.

Comparison of survival prediction models for pancreatic cancer: Cox model versus machine learning models.

A survival prediction model has recently been developed to evaluate the prognosis of resected nonmetastatic pancreatic ductal adenocarcinoma based on a Cox model using two nationwide databases: Surveillance, Epidemiology and End Results (SEER) and Korea Tumor Registry System-Biliary Pancreas (KOTUS-BP). In this study, we applied two machine learning methods-random survival forests (RSF) and support vector machines (SVM)-for survival analysis and compared their prediction performance using the SEER and KOTUS-BP datasets. Three schemes were used for model development and evaluation. First, we utilized data from SEER for model development and used data from KOTUS-BP for external evaluation. Second, these two datasets were swapped by taking data from KOTUS-BP for model development and data from SEER for external evaluation. Finally, we mixed these two datasets half and half and utilized the mixed datasets for model development and validation. We used 9,624 patients from SEER and 3,281 patients from KOTUS-BP to construct a prediction model with seven covariates: age, sex, histologic differentiation, adjuvant treatment, resection margin status, and the American Joint Committee on Cancer 8th edition T-stage and N-stage. Comparing the three schemes, the performance of the Cox model, RSF, and SVM was better when using the mixed datasets than when using the unmixed datasets. When using the mixed datasets, the C-index, 1-year, 2-year, and 3-year time-dependent areas under the curve for the Cox model were 0.644, 0.698, 0.680, and 0.687, respectively. The Cox model performed slightly better than RSF and SVM.

DeepHisCoM: deep learning pathway analysis using hierarchical structural component models.

Many statistical methods for pathway analysis have been used to identify pathways associated with the disease along with biological factors such as genes and proteins. However, most pathway analysis methods neglect the complex nonlinear relationship between biological factors and pathways. In this study, we propose a Deep-learning pathway analysis using Hierarchical structured CoMponent models (DeepHisCoM) that utilize deep learning to consider a nonlinear complex contribution of biological factors to pathways by constructing a multilayered model which accounts for hierarchical biological structure. Through simulation studies, DeepHisCoM was shown to have a higher power in the nonlinear pathway effect and comparable power for the linear pathway effect when compared to the conventional pathway methods. Application to hepatocellular carcinoma (HCC) omics datasets, including metabolomic, transcriptomic and metagenomic datasets, demonstrated that DeepHisCoM successfully identified three well-known pathways that are highly associated with HCC, such as lysine degradation, valine, leucine and isoleucine biosynthesis and phenylalanine, tyrosine and tryptophan. Application to the coronavirus disease-2019 (COVID-19) single-nucleotide polymorphism (SNP) dataset also showed that DeepHisCoM identified four pathways that are highly associated with the severity of COVID-19, such as mitogen-activated protein kinase (MAPK) signaling pathway, gonadotropin-releasing hormone (GnRH) signaling pathway, hypertrophic cardiomyopathy and dilated cardiomyopathy. Codes are available at https://github.com/chanwoo-park-official/DeepHisCoM.

Clinicopathologic and protein markers distinguishing the "polymerase epsilon exonuclease" from the "copy number low" subtype of endometrial cancer.

OBJECTIVE: The need to perform genetic sequencing to diagnose the polymerase epsilon exonuclease (POLE) subtype of endometrial cancer (EC) hinders the adoption of molecular classification. We investigated clinicopathologic and protein markers that distinguish the POLE from the copy number (CN)-low subtype in EC. METHODS: Ninety-one samples (15 POLE, 76 CN-low) were selected from The Cancer Genome Atlas EC dataset. Clinicopathologic and normalized reverse phase protein array expression data were analyzed for associations with the subtypes. A logistic model including selected markers was constructed by stepwise selection using area under the curve (AUC) from 5-fold cross-validation (CV). The selected markers were validated using immunohistochemistry (IHC) in a separate cohort. RESULTS: Body mass index (BMI) and tumor grade were significantly associated with the POLE subtype. With BMI and tumor grade as covariates, 5 proteins were associated with the EC subtypes. The stepwise selection method identified BMI, cyclin B1, caspase 8, and X-box binding protein 1 (XBP1) as markers distinguishing the POLE from the CN-low subtype. The mean of CV AUC, sensitivity, specificity, and balanced accuracy of the selected model were 0.97, 0.91, 0.87, and 0.89, respectively. IHC validation showed that cyclin B1 expression was significantly higher in the POLE than in the CN-low subtype and receiver operating characteristic curve of cyclin B1 expression in IHC revealed AUC of 0.683. CONCLUSION: BMI and expression of cyclin B1, caspase 8, and XBP1 are candidate markers distinguishing the POLE from the CN-low subtype. Cyclin B1 IHC may replace POLE sequencing in molecular classification of EC.

PATHOME-Drug: a subpathway-based polypharmacology drug-repositioning method.

MOTIVATION: Drug repositioning reveals novel indications for existing drugs and in particular, diseases with no available drugs. Diverse computational drug repositioning methods have been proposed by measuring either drug-treated gene expression signatures or the proximity of drug targets and disease proteins found in prior networks. However, these methods do not explain which signaling subparts allow potential drugs to be selected, and do not consider polypharmacology, i.e. multiple targets of a known drug, in specific subparts. RESULTS: Here, to address the limitations, we developed a subpathway-based polypharmacology drug repositioning method, PATHOME-Drug, based on drug-associated transcriptomes. Specifically, this tool locates subparts of signaling cascading related to phenotype changes (e.g. disease status changes), and identifies existing approved drugs such that their multiple targets are enriched in the subparts. We show that our method demonstrated better performance for detecting signaling context and specific drugs/compounds, compared to WebGestalt and clusterProfiler, for both real biological and simulated datasets. We believe that our tool can successfully address the current shortage of targeted therapy agents. AVAILABILITY AND IMPLEMENTATION: The web-service is available at http://statgen.snu.ac.kr/software/pathome. The source codes and data are available at https://github.com/labnams/pathome-drug. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Inclusive Quantification Assay of Serum Des-γ-Carboxyprothrombin Proteoforms for Hepatocellular Carcinoma Surveillance by Targeted Mass Spectrometry.

Hepatocellular carcinoma (HCC) is a malignant cancer with one of the highest mortality rates. Des-γ-carboxyprothrombin (DCP) is an HCC serologic surveillance marker that can complement the low sensitivity of alpha-fetoprotein (AFP). DCP exists in the blood as a mixture of proteoforms from an impaired carboxylation process at glutamic acid (Glu) residues within the N-terminal domain. The heterogeneity of DCP may affect the accuracy of measurements because DCP levels are commonly determined using an immunoassay that relies on antibody reactivity to an epitope in the DCP molecule. In this study, we aimed to improve the DCP measurement assay by applying a mass spectrometry (MS)-based approach for a more inclusive quantification of various DCP proteoforms. We developed a multiple-reaction monitoring-MS (MRM-MS) assay to quantify multiple noncarboxylated peptides included in the various des-carboxylation states of DCP. We performed the MRM-MS assay in 300 patients and constructed a robust diagnostic model that simultaneously monitored three noncarboxylated peptides. The MS-based quantitative assay for DCP had reliable surveillance power, which was evident from the area under the receiver operating characteristic curve (AUROC) values of 0.874 and 0.844 for the training and test sets, respectively. It was equivalent to conventional antibody-based quantification, which had AUROC values at the optimal cutoff (40 mAU/mL) of 0.743 and 0.704 for the training and test sets, respectively. The surveillance performance of the MS-based DCP assay was validated using an independent validation set consisting of 318 patients from an external cohort, resulting in an AUROC value of 0.793. Conclusion: Due to cost effectiveness and high reproducibility, the quantitative DCP assay using the MRM-MS method is superior to antibody-based quantification and has equivalent performance.

Identifying miRNA-mRNA Integration Set Associated With Survival Time.

In the "personalized medicine" era, one of the most difficult problems is identification of combined markers from different omics platforms. Many methods have been developed to identify candidate markers for each type of omics data, but few methods facilitate the identification of multiple markers on multi-omics platforms. microRNAs (miRNAs) is well known to affect only indirectly phenotypes by regulating mRNA expression and/or protein translation. To take into account this knowledge into practice, we suggest a miRNA-mRNA integration model for survival time analysis, called mimi-surv, which accounts for the biological relationship, to identify such integrated markers more efficiently. Through simulation studies, we found that the statistical power of mimi-surv be better than other models. Application to real datasets from Seoul National University Hospital and The Cancer Genome Atlas demonstrated that mimi-surv successfully identified miRNA-mRNA integrations sets associated with progression-free survival of pancreatic ductal adenocarcinoma (PDAC) patients. Only mimi-surv found miR-96, a previously unidentified PDAC-related miRNA in these two real datasets. Furthermore, mimi-surv was shown to identify more PDAC related miRNAs than other methods because it used the known structure for miRNA-mRNA regularization. An implementation of mimi-surv is available at http://statgen.snu.ac.kr/software/mimi-surv.

Development and External Validation of Survival Prediction Model for Pancreatic Cancer Using Two Nationwide Databases: Surveillance, Epidemiology and End Results (SEER) and Korea Tumor Registry System-Biliary Pancreas (KOTUS-BP).

Background/Aims: Several prediction models for evaluating the prognosis of nonmetastatic resected pancreatic ductal adenocarcinoma (PDAC) have been developed, and their performances were reported to be superior to that of the 8th edition of the American Joint Committee on Cancer (AJCC) staging system. We developed a prediction model to evaluate the prognosis of resected PDAC and externally validated it with data from a nationwide Korean database. Methods: Data from the Surveillance, Epidemiology and End Results (SEER) database were utilized for model development, and data from the Korea Tumor Registry System-Biliary Pancreas (KOTUS-BP) database were used for external validation. Potential candidate variables for model development were age, sex, histologic differentiation, tumor location, adjuvant chemotherapy, and the AJCC 8th staging system T and N stages. For external validation, the concordance index (C-index) and time-dependent area under the receiver operating characteristic curve (AUC) were evaluated. Results: Between 2004 and 2016, data from 9,624 patients were utilized for model development, and data from 3,282 patients were used for external validation. In the multivariate Cox proportional hazard model, age, sex, tumor location, T and N stages, histologic differentiation, and adjuvant chemotherapy were independent prognostic factors for resected PDAC. After an exhaustive search and 10-fold cross validation, the best model was finally developed, which included all prognostic variables. The C-index, 1-year, 2-year, 3-year, and 5-year time-dependent AUCs were 0.628, 0.650, 0.665, 0.675, and 0.686, respectively. Conclusions: The survival prediction model for resected PDAC could provide quantitative survival probabilities with reliable performance. External validation studies with other nationwide databases are needed to evaluate the performance of this model.

CORRIGENDUM: Addition of data source: Diagnostic model for pancreatic cancer using a multi-biomarker panel.

[This corrects the article on p. 144 in vol. 100, PMID: 33748028.].

Development of Machine Learning Models to Predict Platinum Sensitivity of High-Grade Serous Ovarian Carcinoma.

To support the implementation of individualized disease management, we aimed to develop machine learning models predicting platinum sensitivity in patients with high-grade serous ovarian carcinoma (HGSOC). We reviewed the medical records of 1002 eligible patients. Patients' clinicopathologic characteristics, surgical findings, details of chemotherapy, treatment response, and survival outcomes were collected. Using the stepwise selection method, based on the area under the receiver operating characteristic curve (AUC) values, six variables associated with platinum sensitivity were selected: age, initial serum CA-125 levels, neoadjuvant chemotherapy, pelvic lymph node status, involvement of pelvic tissue other than the uterus and tubes, and involvement of the small bowel and mesentery. Based on these variables, predictive models were constructed using four machine learning algorithms, logistic regression (LR), random forest, support vector machine, and deep neural network; the model performance was evaluated with the five-fold cross-validation method. The LR-based model performed best at identifying platinum-resistant cases with an AUC of 0.741. Adding the FIGO stage and residual tumor size after debulking surgery did not improve model performance. Based on the six-variable LR model, we also developed a web-based nomogram. The presented models may be useful in clinical practice and research.

Microbiome Markers of Pancreatic Cancer Based on Bacteria-Derived Extracellular Vesicles Acquired from Blood Samples: A Retrospective Propensity Score Matching Analysis.

Novel biomarkers for early diagnosis of pancreatic cancer (PC) are necessary to improve prognosis. We aimed to discover candidate biomarkers by identifying compositional differences of microbiome between patients with PC (n = 38) and healthy controls (n = 52), using microbial extracellular vesicles (EVs) acquired from blood samples. Composition analysis was performed using 16S rRNA gene analysis and bacteria-derived EVs. Statistically significant differences in microbial compositions were used to construct PC prediction models after propensity score matching analysis to reduce other possible biases. Between-group differences in microbial compositions were identified at the phylum and genus levels. At the phylum level, three species (Verrucomicrobia, Deferribacteres, and Bacteroidetes) were more abundant and one species (Actinobacteria) was less abundant in PC patients. At the genus level, four species (Stenotrophomonas, Sphingomonas, Propionibacterium, and Corynebacterium) were less abundant and six species (Ruminococcaceae UCG-014, Lachnospiraceae NK4A136 group, Akkermansia, Turicibacter, Ruminiclostridium, and Lachnospiraceae UCG-001) were more abundant in PC patients. Using the best combination of these microbiome markers, we constructed a PC prediction model that yielded a high area under the receiver operating characteristic curve (0.966 and 1.000, at the phylum and genus level, respectively). These microbiome markers, which altered microbial compositions, are therefore candidate biomarkers for early diagnosis of PC.

Powerful p-value combination methods to detect incomplete association.

Meta-analyses increase statistical power by combining statistics from multiple studies. Meta-analysis methods have mostly been evaluated under the condition that all the data in each study have an association with the given phenotype. However, specific experimental conditions in each study or genetic heterogeneity can result in "unassociated statistics" that are derived from the null distribution. Here, we show that power of conventional meta-analysis methods rapidly decreases as an increasing number of unassociated statistics are included, whereas the classical Fisher's method and its weighted variant (wFisher) exhibit relatively high power that is robust to addition of unassociated statistics. We also propose another robust method based on joint distribution of ordered p-values (ordmeta). Simulation analyses for t-test, RNA-seq, and microarray data demonstrated that wFisher and ordmeta, when only a small number of studies have an association, outperformed existing meta-analysis methods. We performed meta-analyses of nine microarray datasets (prostate cancer) and four association summary datasets (body mass index), where our methods exhibited high biological relevance and were able to detect genes that the-state-of-the-art methods missed. The metapro R package that implements the proposed methods is available from both CRAN and GitHub ( http://github.com/unistbig/metapro ).

Identification of genetic loci affecting body mass index through interaction with multiple environmental factors using structured linear mixed model.

Multiple environmental factors could interact with a single genetic factor to affect disease phenotypes. We used Struct-LMM to identify genetic variants that interacted with environmental factors related to body mass index (BMI) using data from the Korea Association Resource. The following factors were investigated: alcohol consumption, education, physical activity metabolic equivalent of task (PAMET), income, total calorie intake, protein intake, carbohydrate intake, and smoking status. Initial analysis identified 7 potential single nucleotide polymorphisms (SNPs) that interacted with the environmental factors (P value < 5.00 × 10-6). Of the 8 environmental factors, PAMET score was excluded for further analysis since it had an average Bayes Factor (BF) value < 1 (BF = 0.88). Interaction analysis using 7 environmental factors identified 11 SNPs (P value < 5.00 × 10-6). Of these, rs2391331 had the most significant interaction (P value = 7.27 × 10-9) and was located within the intron of EFNB2 (Chr 13). In addition, the gene-based genome-wide association study verified EFNB2 gene significantly interacting with 7 environmental factors (P value = 5.03 × 10-10). BF analysis indicated that most environmental factors, except carbohydrate intake, contributed to the interaction of rs2391331 on BMI. Although the replication of the results in other cohorts is warranted, these findings proved the usefulness of Struct-LMM to identify the gene-environment interaction affecting disease.

Diagnostic model for pancreatic cancer using a multi-biomarker panel.

PURPOSE: Diagnostic biomarkers of pancreatic ductal adenocarcinoma (PDAC) have been used for early detection to reduce its dismal survival rate. However, clinically feasible biomarkers are still rare. Therefore, in this study, we developed an automated multi-marker enzyme-linked immunosorbent assay (ELISA) kit using 3 biomarkers (leucine-rich alpha-2-glycoprotein [LRG1], transthyretin [TTR], and CA 19-9) that were previously discovered and proposed a diagnostic model for PDAC based on this kit for clinical usage. METHODS: Individual LRG1, TTR, and CA 19-9 panels were combined into a single automated ELISA panel and tested on 728 plasma samples, including PDAC (n = 381) and normal samples (n = 347). The consistency between individual panels of 3 biomarkers and the automated multi-panel ELISA kit were accessed by correlation. The diagnostic model was developed using logistic regression according to the automated ELISA kit to predict the risk of pancreatic cancer (high-, intermediate-, and low-risk groups). RESULTS: The Pearson correlation coefficient of predicted values between the triple-marker automated ELISA panel and the former individual ELISA was 0.865. The proposed model provided reliable prediction results with a positive predictive value of 92.05%, negative predictive value of 90.69%, specificity of 90.69%, and sensitivity of 92.05%, which all simultaneously exceed 90% cutoff value. CONCLUSION: This diagnostic model based on the triple ELISA kit showed better diagnostic performance than previous markers for PDAC. In the future, it needs external validation to be used in the clinic.

Clinical Characteristics and Outcomes of COVID-19 Cohort Patients in Daegu Metropolitan City Outbreak in 2020.

BACKGROUND: A coronavirus disease 2019 (COVID-19) outbreak started in February 2020 and was controlled at the end of March 2020 in Daegu, the epicenter of the coronavirus outbreak in Korea. The aim of this study was to describe the clinical course and outcomes of patients with COVID-19 in Daegu. METHODS: In collaboration with Daegu Metropolitan City and Korean Center for Diseases Control, we conducted a retrospective, multicenter cohort study. Demographic, clinical, treatment, and laboratory data, including viral RNA detection, were obtained from the electronic medical records and cohort database and compared between survivors and non-survivors. We used univariate and multi-variable logistic regression methods and Cox regression model and performed Kaplan-Meier analysis to determine the risk factors associated with the 28-day mortality and release from isolation among the patients. RESULTS: In this study, 7,057 laboratory-confirmed patients with COVID-19 (total cohort) who had been diagnosed from February 18 to July 10, 2020 were included. Of the total cohort, 5,467 were asymptomatic to mild patients (77.4%) (asymptomatic 30.6% and mild 46.8%), 985 moderate (14.0%), 380 severe (5.4%), and 225 critical (3.2%). The mortality of the patients was 2.5% (179/7,057). The Cox regression hazard model for the patients with available clinical information (core cohort) (n = 2,254) showed the risk factors for 28-day mortality: age > 70 (hazard ratio [HR], 4.219, P = 0.002), need for O2 supply at admission (HR, 2.995; P = 0.001), fever (> 37.5°C) (HR, 2.808; P = 0.001), diabetes (HR, 2.119; P = 0.008), cancer (HR, 3.043; P = 0.011), dementia (HR, 5.252; P = 0.008), neurological disease (HR, 2.084; P = 0.039), heart failure (HR, 3.234; P = 0.012), and hypertension (HR, 2.160; P = 0.017). The median duration for release from isolation was 33 days (interquartile range, 24.0-46.0) in survivors. The Cox proportional hazard model for the long duration of isolation included severity, age > 70, and dementia. CONCLUSION: Overall, asymptomatic to mild patients were approximately 77% of the total cohort (asymptomatic, 30.6%). The case fatality rate was 2.5%. Risk factors, including older age, need for O2 supply, dementia, and neurological disorder at admission, could help clinicians to identify COVID-19 patients with poor prognosis at an early stage.