We have witnessed the 2-year-long global rampage of COVID-19 caused by the wide spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). However, knowledge about biomarkers of the entire COVID-19 process is limited. Identification of the systemic features of COVID-19 will lead to critical biomarkers and therapeutic targets for early intervention and clinical disease course prediction. Here, we performed a comprehensive analysis of clinical measurements and serum metabolomics in 199 patients with different stages of COVID-19. In particular, our study is the first serum metabolomic analysis of critical rehabilitation patients and critical death patients. We found many differential metabolites in the comparison of metabolomic results between ordinary, severe, and critical patients and uninfected patients. Through the metabolomic results of COVID-19 patients in various stages, and critical rehabilitation patients and critical death patients, we identified a series of differential metabolites as biomarkers, a separate queue and precise distinction, and predicted COVID-19 verification. These differentially expressed metabolites, included 1,2-di-(9Z,12Z-octadecadienoyl)-sn-glycero-3-phosphate, propylparaben, 20-hydroxyeicosatetraenoic acid, triethanolamine, chavicol, disialosyl galactosyl globoside, 1-arachidonoylglycerophosphoinositol, and alpha-methylstyrene, all of which have been identified for the first time as biomarkers in COVID-19 progression. These biomarkers are involved in many pathological and physiological pathways of COVID-19, for example, immune responses, platelet degranulation, and metabolism which might result in pathogenesis. Our results showed valuable information about metabolites obviously altered in COVID-19 patients with different stages, which could shed light on the pathogenesis as well as serve as potential therapeutic agents of COVID-19.