key: cord-0766404-fg6hod8k authors: Hu, Jiao; Han, zhengyuan; Heidari, Ali Asghar; Shou, Yeqi; Ye, Hua; Wang, Liangxing; Huang, Xiaoying; Chen, Huiling; Chen, Yanfan; Wu, Peiliang title: Detection of COVID-19 severity using blood gas analysis parameters and Harris hawks optimized extreme learning machine date: 2021-12-24 journal: Comput Biol Med DOI: 10.1016/j.compbiomed.2021.105166 sha: 0a7fc77db37ed7313ca5351e76be7bb47001f39c doc_id: 766404 cord_uid: fg6hod8k Coronavirus disease-2019 (COVID-19) has made the world more cautious about widespread viruses, and a tragic pandemic that was caused by a novel coronavirus has harmed human beings in recent years. The new coronavirus pneumonia outbreak is spreading rapidly worldwide. We collect arterial blood samples from 51 patients with a COVID-19 diagnosis. Blood gas analysis is performed using a Siemens RAPID Point 500 blood gas analyzer. To accurately determine the factors that play a decisive role in the early recognition and discrimination of COVID-19 severity, a prediction framework that is based on an improved binary Harris hawk optimization (HHO) algorithm in combination with a kernel extreme learning machine is proposed in this paper. This method uses specular reflection learning to improve the original HHO algorithm and is referred to as HHOSRL. The experimental results show that the selected indicators, such as age, PaO(2), SaO(2)%, Na(+), and LAC, are essential for the early accurate assessment of COVID-19 severity by the feature selection method that is proposed in this paper. The simulation results show that the established methodlogy can achieve promising performance. We believe that our proposed model provides an effective strategy for accurate early assessment of COVID-19 and distinguishing disease severity. The International Committee on the Taxonomy of Viruses (ICTV) described a new virus, namely, severe acute respiratory syndrome coronavirus (SARS-CoV-2) [1] . SARS-CoV-2 is believed to be the pathogen that causes viral pneumonia, which has caused a worldwide pandemic. On February 11, 2020, the viral pneumonia that is caused by SARS-CoV-2 was named coronavirus disease 2019 by the World Health Organization (WHO) [2] . The case fatality ratio (CFR) of COVID-19 is significantly lower than those of severe acute respiratory syndrome (SARS) and Middle East respiratory syndrome (MERS) [3, 4] . However, COVID-19's long incubation period of approximately two weeks and the presence of asymptomatic infections increase the risk of COVID-19 virus infection and promote its spread [5] . To date, according to WHO data, COVID-19 has spread globally, with more than 70,000,000 cases of COVID-19 diagnosed in various countries, which have resulted in more than 1,590,000 related deaths (https://www.who.int/). The most common clinical manifestation of 2019-nCoV infection includes fever, dry cough, dyspnea, chest pain, myalgia, and fatigue [6] . According to a previous clinical study, 20 .7%-31.4% of COVID-19 patients developed a severe form of this disease, such as adult respiratory distress syndrome (ARDS). Furthermore, 4.9% to 11.5% of patients with COVID-19 needed advanced life support techniques in the intensive care unit (ICU) [7] . The main characteristics of severe COVID-19 are rapid progression, ARDS, multiple organ dysfunction (MOF), and high fatality rate [8] . According to previous reports, elderly patients with comorbidities are more susceptible to being infected by COVID-19 [9, 10] . Moreover, the mortality rate of elderly patients (> 65 years old) with comorbidities and ARDS is significantly increased [11] . The rapid progression of COVID-19 places substantial strain on health care systems and hospital critical care resources. Severe COVID-19 patients may experience rapid deterioration if not treated in the ICU timely [12] . Therefore, it is vital to conduct accurate and frequent clinical assessments. However, due to resources being stretched thin and lack of experience and prior knowledge of physicians, accurate and frequent assessment is not easy. Effective novel prognostic model systems that can help monitor changes in the condition of COVID-19 patients can guide the effective use of hospital resources. Recently, machine learning methods have been used by many medical workers to help solve medical problems. According to many studies, machine learning (ML) techniques and artificial intelligence (AI) methods have been widely implemented in the diagnosis of COVID-19 [13] [14] [15] [16] , and many studies have applied swarm intelligence methods to COVID-19 image segmentation [17, 18] . Pham et al. [19] trained convolutional neural networks (CNNs) to fine-tune COVID-19 detection in chest slices. Canayaz [20] used deep learning models such as AlexNet, VGG19, GoogleNet, and ResNet to conduct feature extraction and select optimal potential features; two meta-heuristic algorithms, namely, binary particle swarm optimization and binary gray wolf optimization, were used. Al-Falluji [21] et al. utilized X-ray images for deep learning and retrieved essential biomarkers that were related to COVID-19 disease detection. Shaban et al. [22] combined a fuzzy inference engine with deep neural networks (DNNs) to propose a new hybrid diagnostic strategy (HDS) for classifying newly infected individuals to determine whether they are infected. Sun et al. [23] proposed an adaptive feature selection guided deep forest (AFS-DF) approach that is based on chest CT images for COVID-19 classification. Shaban et al. [24] introduced a novel COVID-19 patient detection strategy (CPDS) that J o u r n a l P r e -p r o o f is based on hybrid feature selection and enhanced the KNN classifier, which significantly improved the diagnostic accuracy. Dey et al. [25] built various machine learning models for predicting protein-protein interactions (PPIs) between COVID-19 and human proteins and further validated them by biological experiments. Abraham et al. [26] combined multiple CNN-extracted features with a correlation-based feature selection (CFS) technique and a Bayesian classifier for COVID-19 prediction. Liu et al. [27] developed and validated a complete machine learning framework for chest CT images to distinguish COVID-19 from global pneumonia (GP). Tuncer et al. [28] proposed a novel intelligent computer vision method for automatic detection of the COVID-19 virus and conducted 10-fold cross-validation based on the SVM classifier, which showed a classification accuracy of 100.0%. Casiraghi et al. [29] proposed an interpretable machine learning system that provides simple decision criteria for use by clinicians in assessing patient risk. Novitasari et al. [30] used convolutional neural network methods such as feature extraction and support vector machines (SVM) as classification methods to detect whether a patient being examined was healthy, coronavirus-positive, or only had pneumonia. In this work, we proposed a new machine learning framework, namely, the swarm intelligence augmented kernel extreme learning machine (KELM) [25] , for predicting the severity of COVID-19. Notably, for the first time, the binary Harris hawk optimization algorithm (HHO), which is improved using specular reflection learning, is used in combination with the KELM classifier for feature selection. The experimental results and simulation results demonstrate the superior performance of the method. The results show that bHHOSRL_KELM performs well in determining which factors play a decisive role in the outcome in COVID-19 diagnosis. The proposed bHHOSRL_KELM was compared with other classifiers, such as FKNN, KNN, MLP, and SVM, in terms of four classification metrics, namely, classification accuracy, sensitivity, specificity, and Mathews correlation coefficient (MCC). In addition, nine other feature selection methods that are based on swarm intelligence algorithms were used to evaluate the performance of the developed bHHOSRL_KELM based on the fitness values throughout the iterative process. The experimental results show that the developed bHHOSRL_KELM has high predictive performance. The main contributions of this study are as follows: (1) an efficient diagnostic aid for COVID-19 in blood specimens was developed; (2) a promising hybrid model was proposed, and the potential of KELM was enhanced using an improved HHO method; and (3) the bHHOSRL_KELM feature selection method was used to identify the most critical features effectively. The remainder of this paper is organized as bellow. The first section describes the current state of research on novel coronaviruses and the application of machine learning methods to novel coronavirus research. The second section presents the dataset that is used in this paper. The third section presents the improved HHO method and its integration with KELM. The fourth section analyze the proposed method and presents the experimental results on COVID-19. The fifth section discusses the proposed method and analyzes the final medical data results. We separated the COVID-19 patient dataset into severe COVID-19, which corresponded to 21 samples, and nonsevere COVID-19, which corresponded to 30 samples. For a diagnosis of severe COVID-19, the following requirements should be satisfied: (1) respiratory rate greater than 30/min (respiratory distress); (2) resting oxygen saturation lower than 93%; and (3) oxygenation index (OI) lower than 300 mmHg. We collected arterial blood samples from all 51 patients with a COVID-19 diagnosis. Blood gas analysis was performed by a SIEMENS RAPID Point 500 blood gas analyzer. The basic clinical information and 22 blood gas analysis parameters (features) are listed in Table 1 . All continuity variables are presented as the mean ± standard deviation (SD) and were analyzed with SPSS Statistics 24.0. An independent sample t test was used to analyze the continuous variables (age and blood gas analysis parameters). p < 0.05 was considered statistically significant. The statistical analysis results are presented in Table 2 . In this study, the improved HHO with specular reflection learning (HHOSRL) mechanism is used as a learning algorithm in the parcel-based feature selection method, and the binary HHOSRL method is used as a feature selection tool to identify the critical features and evaluate the feature subsets using the KELM model. HHO [32] is a new meta-heuristic algorithm that was proposed by Heidari et al. in 2019. Since its introduction, it has been used in many situations, such as parameter estimation of photovoltaic models. It is a superior method compared to all other algorithms in recent years [33] [34] [35] [36] , such as the colony predation algorithm (CPA) [37] . HHO has the unique feature that Harris hawks can cooperate in groups to chase prey and adjust the chase pattern according to the dynamics of the situation and the escape pattern of the prey. Since its introduction, HHO has been applied to solve many problems such as parameters identification of photovoltaic cells and modules [38] [39] [40] [41] , feature selection [42, 43] , optimizing the machine learning models [16, 44] , engineering design problems [41, [45] [46] [47] , web service composition [46] , bankruptcy prediction [48] , and multi-objective problems [49] . Harris hawks perch at random in a specified location and find their prey through 2 strategies, which are expressed mathematically in Eq. (1). where The HHO algorithm shifts between exploration and various exploitation behaviors based on the prey's escape energy, which is defined as follows: where 0 is the initial energy of the prey, which is a random number in [-1,1] that is updated automatically at each iteration; is the number of iterations; and _ is the maximum number of iterations. The search phase is entered when | | ≥ 1, and the exploitation phase is entered when | | < 1. Define r as a random number between [0,1] for choosing the exploitation strategy. When 0.5 ≤ | | < 1 and ≥ 0.5, a soft besiege strategy is used to update the position. When | | < 0.5 and ≥ 0.5, a hard besiege strategy is used to update the position. When 0.5 ≤ | | < 1 and r<0.5, a soft besiege strategy with a progressive rapid dive is used to update the position. Eq. (7) where ( * ) is the fitness function; S is a 2-dimensional random vector with elements that are random numbers in [0,1]; and ( * ) is the mathematical expression for Levy flight. When | | < 0.5 and r<0.5, a soft besiege strategy with a progressive rapid dive is used to update the position. In this paper, we address a binary optimization case, which corresponds to whether a feature is selected or not. In a discrete binary space, a solution is restricted to two numbers (0 and 1). For our case, an upgraded binary HHO is proposed. In this study, we represent a solution as a d To overcome this problem, the method discretizes the updated hawks' position vector to a binary vector. The updating equation is as follows: where is a random number in [0,1] and ( + 1) is the updated location at the t-th iteration. The expression for ( ) is presented below. Eq. (12) This paper updates the solution using a variation operator based on specular reflection learning (SRL) to further explore more combinations of properties. Eq. (13) where ( , ) is the binary value of the −th dimension of the i-th individual in the t-th iteration after mutation and is a variable that controls the range of the new position after specular reflection, which is calculated via Eq. (14) . where 1 and 2 are two random numbers that are uniformly distributed between 0 and 1 and 0 is the neighborhood radius, which is also called the elasticity factor. From Eq. (14), if 1 > 2 , will appear in the right neighborhood of ; otherwise, will appear in the left neighborhood of . Therefore, it is reasonable to set 0 to [0,1] to balance the search space of the left and right neighbors. Additionally, is also set to [0,1], namely, all possible values within a radius 0 of the neighborhood are obtained. Initialize a set of binary hawks (X , i=1…N) and the probability of mutation ; Calculate the fitness with the selected feature subset for each hawk; Set the individual with the lowest fitness to ; Update each individual's escape ability and jump strength ; After the excitation function and the number of implicit layer nodes are established, the following 3 steps are conducted to train the ELM on a training dataset. Step 1: Randomly generate input weights and , 1 ≤ ≤ ; Step 2: Calculate the output matrix of the hidden layer ; Step 3: Calculate the output weight matrix = + ; where + is the Moore-Penrose generalized inverse of the implicit layer output matrix . When is nonsingular, + = ( ) −1 . To eliminate errors in the results of the "sick matrix," according to the strategy of ridge regression, a regularization factor is introduced, and the least-squares solution of the output weights of the network is Therefore, the corresponding ELM output function is Eq. (21) where γ is the nuclear parameter of the RBF kernel function. The two key parameters have been shown to have a large impact on the performance of KELM on many problems, such as second major prediction [50] , medical diagnosis [51] [52] [53] [54] [55] [56] , financial stress prediction [57, 58] , and recognition of foreign fibers in cotton [59] . To find more representative attributes in the dataset to help us make a medical diagnosis, this paper uses HHOSRL-KELM as a feature selection method to choose the optimal feature subset. First, we use HHOSRL as the optimization algorithm to find the optimal subset, and after finding the optimal features, we use KELM as the classifier for the classification task. Figure 1 shows a HHOSRL-KELM flowchart. It illustrates the process of finding key factors for a new coronary pneumonia blood sample. The main steps of the HHOSRL-KEML method are described below. J o u r n a l P r e -p r o o f Step 1: Initialize the parameters of the HHOSRL method, such as population size, search space boundary, variance probability, maximum number of iterations, and initial escape energy. Step 2: Randomly initialize the binary population of Harris hawks. Step 3: Use the agent's binary value in each dimension to represent the subset selection of the dataset (1 indicates that the feature is selected, and 0 indicates that it is not selected). Step 4: Calculate the fitness value of the selected feature subset for each hawk as follows: where is the classification error rate of KELM, | | is the number of features in the data sample, | | represents the number of features in the selected feature subset, and are two weights that reflect the importance of the classification error rate and the length of the selected features, respectively In this paper, we set = 0.99 and = 0.01, which are values that are commonly used in many works [42, 43] . Step 5: Update the population of agents according to the HHOSRL algorithm. Step 6: Select the individual with the smallest fitness value as the optimal solution. Step 7: Determine whether the termination condition has been reached, namely, whether the maximum number of iterations has been reached. If yes, proceed to the next step; otherwise, repeat step 3 until the termination condition is met. Step 8: Return the final optimal solution as the selected feature subset. Step 9: Use the final feature subset as the input parameter of KELM to obtain the final classification result. Step 10: Using the classification results that were obtained in Step 9, calculate the classification error accuracy, number of selected feature subsets, sensitivity, specificity, and other evaluation criteria. J o u r n a l P r e -p r o o f We utilized four mutual rules based on the confusion matrix to validate the efficacy of the classifier. Full definitions of these metrics are provided in [64, 65] . Here, we present their formulas to avoid discussions that are outside the scope of this study: where all variables are as defined in [65] . In addition, the MCC was used to carefully evaluate the classifier's performance because it provides an objective predictive valuation [66] . It is well known that classification methods play an important role in the performance of wrapped feature selection approches [15, [67] [68] [69] [70] [71] [72] . Therefore, to evaluate the performance of the proposed HHOSRL method, it was combined with five classification methods, namely, SVM, k-nearest neighbor J o u r n a l P r e -p r o o f (KNN), fuzzy k-nearest neighbor (FKNN), multilayer perceptron (MLP), and KELM, and applied to the feature selection of blood samples. To evaluate the reliability of the results, we analyzed these five feature selection methods on blood data in terms of four aspects: MCC, ACC, sensitivity, and specificity. The parameter settings of these five classification methods are presented in Table 3 . Figure 2 shows error bar graphs for the five methods on these four criteria. The error line is the magnitude of the mean confidence space, which represents the standard error and is used to visualize the magnitude of the standard deviation. It can be easily observed from the graph that all five classification methods achieve a result of 1 on this criterion and that the error results are all 0. This is sufficient to prove that HHOSRL has a strong learning performance and can accurately find the globally optimal solution. It also shows that HHOSRL has strong stability. In addition to these four evaluation criteria, to compare the advantages and disadvantages of HHOSRL using these five classifiers, the five feature selection algorithms were also analyzed in terms of time consumption, fitness value, size of the selected feature subset and classification error rate. iterations in Figure 4 shows that bHHOSRL_KELM has the fastest convergence speed and the smallest convergence value. However, according to the error line, these four algorithms have many errors. bHHO does not reach the optimum on any of the four criteria. Finally, the error bar graph of bHHOSRL_KELM shows that the overall accuracy is the highest and the error is the smallest; hence, bHHOSRL_KELM has the best results on these ten feature selection methods. J o u r n a l P r e -p r o o f Figure 8 . According to the graph, age, PaO2, SaO2%, Na +, and LAC were selected the most times. PaO2 was selected 57 times; hence, PaO2 was the most important factor that influenced the determination of whether a patient was infected with the new coronavirus. The next most important attribute was SaO2%, which was selected ten times. This also suggests that SaO2% plays a role in the final outcome diagnosis. Age and Na + were checked nine times, thereby indicating that age and Na + can facilitate final diagnosis. Finally, LAC was selected eight times. This shows that LAC is also a factor that cannot be ignored in diagnosis. Figure 9 shows the convergence graphs of these ten algorithms, which were used to analyze the convergence accuracy and convergence speed of bHHOSRL_KELM. As shown in the graph, in terms of both convergence speed and convergence accuracy, bHHOSRL_KELM performs the best, and BGSA and bGWO perform the next best. In this study, a KELM-based feature selection method was used to screen a dataset of COVID-19 patients. Then, we screened several key features, including age, PaO2, SaO2%, Na +, and LAC. Subsequently, the bHHOSRL_KELM model was constructed to accurately assess the severity of COVID-19 and monitor its early progression. Therefore, we believe that the use of the bHHOSRL_KELM model can make a more accurate clinical decision-making. Due to its great optimization capability, the proposed HHOSRL algorithm can also be applied to solve other problems, such as problems in video deblurring [73] , microgrid planning [74] , information retrieval services [75] [76] [77] , image dehazing [78] , kayak cycle phase segmentation [79] , human motion capture [80] , fault detection [81] , video coding optimization [82] , outlier detection [83] , location-based services [84, 85] , image retrieval [86] , multivariate time series analysis [87] and multi-objective problems [88] . Previous studies have shown that age is an important independent indicator for predicting the prognosis of SARS and MERS [89, 90] . In the study, Smits et al. found that elderly macaques that were vaccinated with SARS coronavirus showed stronger host innate responses than younger macaques, which may be related to the increased differential expression of inflammation-related genes and reduced expression of type I interferon β [91] . Similar to SARS and MERS, Zhou et al. found that the mortality rate of COVID-19 patients increases with age [92] . Through multivariate analysis, other studies show that in COVID-19 patients, advanced age is an independent risk factor for pneumonia. Moreover, patients with ARDS are older than those without ARDS [93] . Furthermore, another study found that elderly patients have a higher probability of developing sepsis, which may be related to age-dependent defects in lymphocyte function and excessive production of type 2 cytokines [94] . In summary, age is an essential predictor of the prognosis of COVID-19 patients. SaO2 is the percentage of the volume of HbO2 that is bound by oxygen in the blood relative to the total volume of bound hemoglobin. PO2 is the tension that is produced by oxygen that is physically dissolved in the blood [95] . Xie et al. found that hypoxemia in COVID-19 patients is associated with mortality [96] . Hypoxia has been reported to increase angiotensin I-converting enzyme 2 (ACE-2) expression at the transcription and protein levels in human cells [97, 98] . According to reports, ACE2 is the target receptor that SARS-CoV-2 enters. Moreover, the binding affinity of the SARS-CoV-2 spike protein to the ACE2 receptor is 10-20 times higher than that of SARS-CoV [99] . The high affinity of ACE2 and the SARS-CoV-2 spike protein may explain some complications of COVID-19 patients, including acute renal failure, cardiovascular and cerebrovascular diseases. Diffuse endotheliitis and microthrombi formation may be an important pathogenesis mechanism for COVID-19 [100, 101] . Therefore, hypoxia may aggravate the severity of COVID-19 by upregulating the target receptor for virus entry. Additionally, in the presence of microthrombi, hypoxemia will cause a higher degree of hypoxia and damage to peripheral tissues [102] . Based on the above considerations, SaO2 and PO2 may be powerful predictors of the severity of COVID-19. Sodium ions (Na + ), which are the most abundant cations in the extracellular fluid, are essential for maintaining the volume of the extracellular fluid, regulating the acid-base balance, maintaining the normal osmotic pressure and cell physiological functions, and participating in the normal physiological activities of the nerve and muscle system [103] . In early COVID-19 studies, the sodium ion concentration in the severe COVID-19 group was significantly lower than the concentration in the nonsevere COVID-19 group [104, 105] . In a study that was based on 59 COVID-19 patients, Huang J o u r n a l P r e -p r o o f and colleagues found that the sodium ion concentration in ICU care groups was lower than that in no ICU care groups [9] . Our study also found that the average sodium ion concentration in severe patients, namely, 131.00 mmol/L, was lower than that in nonsevere patients, namely, 136.23 mmol/L (p<0.01). Therefore, serum sodium ion concentration may also be an essential indicator for predicting severe COVID-19. Lactic acid (LAC) is an intermediate product of sugar metabolism. Increased glucose metabolism or decreased pyruvate metabolism will increase the production of lactic acid. In sepsis, the inflammation response is associated with increased glycolysis and impaired pyruvate dehydrogenase. At this time, the metabolism of pyruvate is limited; hence, the concentration of pyruvate increases. To maintain a normal ratio of pyruvate to lactic acid, the lactic acid concentration will increase [106] . Therefore, LAC can be used as an indicator of the severity of inflammation in patients. A study of 1461 patients showed that the mortality rate of COVID-19 patients with high lactate levels was approximately twice that of patients with low lactate levels [107] . Although cell breakdown under conditions of critical illness may also play a role, the increase in lactic acid in COVID-19 patients may be caused by increased glycolytic activity in multiple cell types [108] . Moreover, an inflammatory factor, namely, interleukin 6, which induces lactic acid production, is present at high concentration in patients with COVID-19. In conclusion, LAC can be used as a predictor of the severity of COVID-19 disease. Up to now, few studies have utilized blood gas analysis parameters or clinical information to distinguish the severity of the COVID-19. To the best of our knowledge, this is the first study to predict the prognosis of COVID-19 using machine learning methods based on age, PaO2, SaO2%, Na + , and LAC. However, this study has several limitations. First, our data originated from a single third-level grade-A hospital in the east of China, and the sample size is not large enough. This limits the accuracy of the prediction. We will expand the sample size in future studies. Second, multicenter, large independent/external datasets and prospective studies need to be conducted to validate the results. This article begins with a description of the current status of COVID-19 and the tremendous strain it places on health care systems and hospital critical care resources. It also describes in detail the sources of the data that are used in this paper. Later, the methodology that is proposed is described in detail. First, the basic HHO algorithm is described. Then, the improvement strategy that is proposed in this paper is presented. Finally, its fusion with the KELM classifier is used to filter out essential features from a blood sample dataset. To evaluate the performance of the proposed method, HHOSRL is fused with various classifiers, and it is demonstrated that HHOSRL produces satisfactory results using many of these classifiers. An accuracy of almost 100.00% is achieved. Then, it is determined that fusion with the KELM classifier enables HHOSRL to perform best on blood samples. After determining the best classifier, the proposed method is compared with various swarm intelligence feature selection methods. It is found that the proposed bHHOSRL_KELM achieves almost 100.00% specificity, accuracy, sensitivity, and MCC. The time consumption of bHHOSRL_KELM is much less than those of other feature selection methods. Finally, the five features that were selected most frequently in the experiment were selected by this method. Moreover, the roles of these five features in medical diagnosis are discussed. For future work, several issues merit further consideration. Additional influencing factors and coefficients can be included in the investigation, and in our paper, only the available data are presented. Moreover, parallel computing can be used to decrease the computational load during various applications. Additionally, more data samples can be collected to construct a more efficient and reliable framework. Finally, bHHOSRL_KELM can be used for the diagnosis of other diseases, and the algorithm's application scope can be expanded, e.g., to clustering and CT image segmentation. A whale optimization algorithm with chaos mechanism based on quasi-opposition for global optimization problems Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and coronavirus disease-2019 (COVID-19): The epidemic and the challenges SARS-CoV-2, SARS-CoV, and MERS-COV: A comparative overview COVID-19, SARS and MERS: are they closely related? 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 COVID-19): A literature review Vibrational responses of a MHC viscoelastic thick annular plate in thermal environment using GDQ method, Mechanics Based Design of Structures and Machines Epidemiological and clinical characteristics of 99 cases of 2019 novel coronavirus pneumonia in Wuhan, China: a descriptive study Clinical features of patients infected with 2019 novel coronavirus in Diabetic patients with COVID-19 infection are at higher risk of ICU admission and poor short-term outcome Clinical course and outcomes of critically ill patients with SARS-CoV-2 pneumonia in Wuhan, China: a single-centered, retrospective, observational study COVID-19 and Italy: what next? Analysis of COVID-19 severity from the perspective of coagulation index using evolutionary machine learning with enhanced brain storm optimization Early Recognition and Discrimination of COVID-19 Severity Using Slime Mould Support Vector Machine for Medical Decision-Making Evolutionary warning system for COVID-19 severity: Colony predation algorithm enhanced extreme learning machine Diagnosing Coronavirus Disease 2019 (COVID-19): Efficient Harris Hawks-Inspired Fuzzy K-Nearest Neighbor Prediction Methods Gaussian Barebone Salp Swarm Algorithm with Stochastic Fractal Search for medical image segmentation: A COVID-19 case study Ant colony optimization with Cauchy and greedy Levy mutations for multilevel COVID 19 X-ray image segmentation Classification of COVID-19 chest X-rays with deep learning: new models or fine tuning? Diagnosis of COVID-19 using deep neural networks and meta-heuristic-based feature selection on X-ray images Automatic Detection of COVID-19 Using Chest X-Ray Images and Modified ResNet18-Based Convolution Neural Networks Detecting COVID-19 patients based on fuzzy inference engine and Deep Neural Network Adaptive Feature Selection Guided Deep Forest for COVID-19 Classification With Chest CT Patients Detection Strategy (CPDS) based on hybrid feature selection and enhanced KNN classifier, Knowledge-Based Systems Machine learning techniques for sequence-based prediction of viral-host interactions between SARS-CoV-2 and human proteins Computer-aided detection of COVID-19 from X-ray images using multi-CNN and Bayesnet classifier Differentiating novel coronavirus pneumonia from general pneumonia based on machine learning An automated Residual Exemplar Local Binary Pattern and iterative ReliefF based COVID-19 detection method using chest X-ray image Explainable Machine Learning for Early Assessment of COVID-19 Risk Prediction in Emergency Departments DETECTION OF COVID-19 CHEST X-RAY USING SUPPORT VECTOR MACHINE AND CONVOLUTIONAL NEURAL NETWORK An effective machine learning approach for identifying non-severe and severe coronavirus disease 2019 patients in a rural Chinese population: the wenzhou retrospective study Harris hawks optimization: Algorithm and applications, Future Generation Computer Systems-the Synthesizing Multi-Layer Perceptron Network with Ant Lion Biogeography-Based Dragonfly Algorithm Evolutionary Strategy Invasive Weed and League Champion Optimization Hybrid Algorithms in Predicting Heating Load in J o u r n a l P r e -p r o o f Residential Buildings Suggesting a Stochastic Fractal Search Paradigm in Combination with Artificial Neural Network for Early Prediction of Cooling Load in Residential Buildings An Innovative Metaheuristic Strategy for Solar Energy Management through a Neural Networks Framework, Energies Electrical Power Prediction through a Combination of Multilayer Perceptron with Water Cycle Ant Lion and Satin Bowerbird Searching Optimizers The Colony Predation Algorithm Parameters identification of photovoltaic cells and modules using diversification-enriched Harris hawks optimization with chaotic drifts Boosted mutation-based Harris hawks optimizer for parameters identification of single-diode solar cell models, Energy Conversion and Management Orthogonally adapted Harris hawks optimization for parameter estimation of photovoltaic models Horizontal and vertical crossover of Harris hawk optimizer with Nelder-Mead simplex for parameter estimation of photovoltaic models Boosted binary Harris hawks optimizer and feature selection Memory-based Harris hawk optimization with learning agents: a feature selection approach Predicting Entrepreneurial Intention of Students: An Extreme Learning Machine With Gaussian Barebone Harris Hawks Optimizer Multipopulation differential evolution-assisted Harris hawks optimization: Framework and case studies Survival exploration strategies for Harris Hawks Optimizer Multi-strategy Gaussian Harris hawks optimization for fatigue life of tapered roller bearings Enhanced Harris hawks optimization with multi-strategy for global optimization tasks, Expert Systems with Applications Elitist non-dominated sorting Harris hawks optimization: Framework and developments for multi-objective problems Evolving an optimal kernel extreme learning machine by using an enhanced grey wolf optimization strategy An enhanced Bacterial Foraging Optimization and its application for training kernel extreme learning machine A New Kernel Extreme Learning Machine Framework for Somatization Disorder Diagnosis A New Effective Machine Learning Framework for Sepsis Diagnosis Toward an optimal kernel extreme learning machine using a chaotic moth-flame optimization strategy with applications in medical diagnoses An enhanced grey wolf optimization based feature selection wrapped kernel extreme learning machine for medical diagnosis, Computational and mathematical methods in medicine An efficient hybrid kernel extreme learning machine approach for early diagnosis of Parkinson‫׳‬s disease An improved grasshopper optimization algorithm with application to financial stress prediction An Effective Computational Model for Bankruptcy Prediction Using Kernel Extreme Learning Machine Approach An Efficient and Effective Automatic Recognition System for Online Recognition of Foreign Fibers in Cotton Moth-flame optimization algorithm: A novel nature-inspired heuristic paradigm, Knowledge-Based Systems The Whale Optimization Algorithm Extreme Learning Machine for Regression and Multiclass Classification Efficient multi-population outpost fruit fly-driven optimizers: Framework and advances in support vector machines An effective model for predicting serum albumin levels in hemodialysis patients Extreme learning machine for regression and multiclass classification Dispersed foraging slime mould algorithm: Continuous and binary variants for global optimization and wrapper-based feature selection, Knowledge-Based Systems Chaotic diffusionlimited aggregation enhanced grey wolf optimizer: Insights, analysis, binarization, and feature selection An effective model for predicting serum albumin level in hemodialysis patients MFeature: Towards high performance evolutionary tools for feature selection Evolving fuzzy k-nearest neighbors using an enhanced sine cosine algorithm: Case study of lupus nephritis Orthogonal learning covariance matrix for defects of grey wolf optimizer: Insights, balance, diversity, and feature selection, Knowledge-Based Systems Recursive Neural Network for Video Deblurring, IEEE Transactions on Circuits and Systems for Video Technology A Risk-Averse Conic Model for Networked Microgrids Planning With Reconfiguration and Reorganizations A user sensitive subject protection approach for book search service A dummy-based user privacy protection approach for text information retrieval, Knowledge-Based Systems An effective approach for the protection of user commodity viewing privacy in e-commerce website Pyramid Channel-based Feature Attention Network for image dehazing, Computer Vision and Image Understanding Sensor Combination Selection Strategy for Kayak Cycle Phase Segmentation Based on Body Sensor Networks Sensor network oriented human motion capture via wearable intelligent system A Novel Advancing Signal Processing Method Based on Coupled Multi-Stable Stochastic Resonance for Fault Detection Video coding optimization in AVS2, Information Processing & Management Anomaly Detection With Representative Neighbors Constructing dummy query sequences to protect location privacy and query privacy in location-based services A location privacy-preserving system based on query range cover-up for location-based services Adaptive Hashing With Sparse Matrix Factorization Evaluation of constraint in photovoltaic cells using ensemble multi-strategy shuffled frog leading algorithms An enhanced fast non-dominated solution sorting genetic algorithm for multi-objective problems Outcomes and prognostic factors in 267 patients with severe acute respiratory syndrome in Hong Kong Predictors of mortality in Middle East respiratory syndrome (MERS) Exacerbated innate host response to SARS-CoV in aged non-human primates Clinical course and risk factors for mortality of adult inpatients with COVID-19 in Wuhan, China: a retrospective cohort study Association Between Clinical Manifestations and Prognosis in Patients with COVID-19 The immunopathogenesis of sepsis in elderly patients Relating oxygen partial pressure, saturation and content: the haemoglobin-oxygen dissociation curve Association Between Hypoxemia and Mortality in Patients With COVID-19 Role of HIF-1alpha in the regulation ACE and ACE2 expression in hypoxic human pulmonary artery smooth muscle cells Transcriptome-based network analysis reveals renal cell type-specific dysregulation of hypoxia-associated transcripts Cryo-EM structure of the 2019-nCoV spike in the prefusion conformation Endothelial cell infection and endotheliitis in COVID-19 COVID-19-Related Severe Hypercoagulability in Patients Admitted to Intensive Care Unit for Acute Respiratory Failure Considerations for target oxygen saturation in COVID-19 patients: are we under-shooting? Integration of acid-base and electrolyte disorders Clinical Characteristics of Coronavirus Disease 2019 in China Epidemiologic and clinical characteristics of 91 hospitalized patients with COVID-19 in Zhejiang, China: a retrospective, multi-centre case series The ten pitfalls of lactate clearance in sepsis Early predictors of in-hospital mortality in patients with COVID-19 in a large American cohort Characterization of the Inflammatory Response to Severe COVID-19 Illness This work was supported in part by the Natural Science Foundation of Zhejiang Province (LZ22F020005), National Natural Science Foundation of China (62076185 and U1809209). We acknowledge the comments of the reviewers.