Machine Learning Approaches to Identify Patient Comorbidities and Symptoms That Increased Risk of Mortality in COVID-19

Personalized Assessment of Mortality Risk and Hospital Stay Duration in Hospitalized Patients with COVID-19 Treated with Remdesivir: A Machine Learning Approach

Background: Despite advancements in vaccination, early treatments, and understanding of SARS-CoV-2, its impact remains significant worldwide. Many patients require intensive care due to severe COVID-19. Remdesivir, a key treatment option among viral RNA polymerase inhibitors, lacks comprehensive studies on factors associated with its effectiveness. Methods: We conducted a retrospective study in 2022, analyzing data from 252 hospitalized COVID-19 patients treated with remdesivir. Six machine learning algorithms were compared to predict factors influencing remdesivir's clinical benefits regarding mortality and hospital stay. Results: The extreme gradient boost (XGB) method showed the highest accuracy for both mortality (95.45%) and hospital stay (94.24%). Factors associated with worse outcomes in terms of mortality included limitations in life support, ventilatory support needs, lymphopenia, low albumin and hemoglobin levels, flu and/or coinfection, and cough. For hospital stay, factors included vaccine doses, lung density, pulmonary radiological status, comorbidities, oxygen therapy, troponin, lactate dehydrogenase levels, and asthenia. Conclusions: These findings underscore XGB's effectiveness in accurately categorizing COVID-19 patients undergoing remdesivir treatment.

Predictive Model for Mortality in Severe COVID-19 Patients across the Six Pandemic Waves

The impact of SARS-CoV-2 infection remains substantial on a global scale, despite widespread vaccination efforts, early therapeutic interventions, and an enhanced understanding of the disease's underlying mechanisms. At the same time, a significant number of patients continue to develop severe COVID-19, necessitating admission to intensive care units (ICUs). This study aimed to provide evidence concerning the most influential predictors of mortality among critically ill patients with severe COVID-19, employing machine learning (ML) techniques. To accomplish this, we conducted a retrospective multicenter investigation involving 684 patients with severe COVID-19, spanning from 1 June 2020 to 31 March 2023, wherein we scrutinized sociodemographic, clinical, and analytical data. These data were extracted from electronic health records. Out of the six supervised ML methods scrutinized, the extreme gradient boosting (XGB) method exhibited the highest balanced accuracy at 96.61%. The variables that exerted the greatest influence on mortality prediction encompassed ferritin, fibrinogen, D-dimer, platelet count, C-reactive protein (CRP), prothrombin time (PT), invasive mechanical ventilation (IMV), PaFi (PaO2/FiO2), lactate dehydrogenase (LDH), lymphocyte levels, activated partial thromboplastin time (aPTT), body mass index (BMI), creatinine, and age. These findings underscore XGB as a robust candidate for accurately classifying patients with COVID-19.

Predicting Multimorbidity Using Saudi Health Indicators (Sharik) Nationwide Data: Statistical and Machine Learning Approach

The Saudi population is at high risk of multimorbidity. The risk of these morbidities can be reduced by identifying common modifiable behavioural risk factors. This study uses statistical and machine learning methods to predict factors for multimorbidity in the Saudi population. Data from 23,098 Saudi residents were extracted from the "Sharik" Health Indicators Surveillance System 2021. Participants were asked about their demographics and health indicators. Binary logistic models were used to determine predictors of multimorbidity. A backpropagation neural network model was further run using the predictors from the logistic regression model. Accuracy measures were checked using training, validation, and testing data. Females and smokers had the highest likelihood of experiencing multimorbidity. Age and fruit consumption also played a significant role in predicting multimorbidity. Regarding model accuracy, both logistic regression and backpropagation algorithms yielded comparable outcomes. The backpropagation method (accuracy 80.7%) was more accurate than the logistic regression model (77%). Machine learning algorithms can be used to predict multimorbidity among adults, particularly in the Middle East region. Different testing methods later validated the common predicting factors identified in this study. These factors are helpful and can be translated by policymakers to consider improvements in the public health domain.

Effect of Cytomegalovirus Reactivation on Inflammatory Status and Mortality of Older COVID-19 Patients

Herpesviridae reactivation such as cytomegalovirus (CMV) has been described in severe COVID-19 (COronaVIrusDisease-2019). This study aimed to understand if CMV reactivation in older COVID-19 patients is associated with increased inflammation and in-hospital mortality. In an observational single-center cohort study, 156 geriatric COVID-19 patients were screened for CMV reactivation by RT-PCR. Participants underwent a comprehensive clinical investigation that included medical history, functional evaluation, laboratory tests and cytokine assays (TNF-α, IFN-α, IL-6, IL-10) at hospital admission. In 19 (12.2%) of 156 COVID-19 patients, CMV reactivation was detected. Multivariate Cox regression models showed that in-hospital mortality significantly increased among CMV positive patients younger than 87 years (HR: 9.94, 95% CI: 1.66–59.50). Other factors associated with in-hospital mortality were C-reactive protein (HR: 1.17, 95% CI: 1.05–1.30), neutrophil count (HR: 1.20, 95% CI: 1.01–1.42) and clinical frailty scale (HR:1.54, 95% CI: 1.04–2.28). In patients older than 87 years, neutrophil count (HR: 1.13, 95% CI: 1.05–1.21) and age (HR: 1.15, 95% CI: 1.01–1.31) were independently associated with in-hospital mortality. CMV reactivation was also correlated with increased IFN-α and TNF-α serum levels, but not with IL-6 and IL-10 serum changes. In conclusion, CMV reactivation was an independent risk factor for in-hospital mortality in COVID-19 patients younger than 87 years old, but not in nonagenarians.

Identification of Comorbidities, Genomic Associations, and Molecular Mechanisms for COVID-19 Using Bioinformatics Approaches

Several studies have been done to identify comorbidities of COVID-19. In this work, we developed an analytical bioinformatics framework to reveal COVID-19 comorbidities, their genomic associations, and molecular mechanisms accomplishing transcriptomic analyses of the RNA-seq datasets provided by the Gene Expression Omnibus (GEO) database, where normal and infected tissues were evaluated. Using the framework, we identified 27 COVID-19 correlated diseases out of 7,092 collected diseases. Analyzing clinical and epidemiological research, we noticed that our identified 27 diseases are associated with COVID-19, where hypertension, diabetes, obesity, and lung cancer are observed several times in COVID-19 patients. Therefore, we selected the above four diseases and performed assorted analyses to demonstrate the association between COVID-19 and hypertension, diabetes, obesity, and lung cancer as comorbidities. We investigated genomic associations with the cross-comparative analysis and Jaccard's similarity index, identifying shared differentially expressed genes (DEGs) and linking DEGs of COVID-19 and the comorbidities, in which we identified hypertension as the most associated illness. We also revealed molecular mechanisms by identifying statistically significant ten pathways and ten ontologies. Moreover, to understand cellular physiology, we did protein-protein interaction (PPI) analyses among the comorbidities and COVID-19. We also used the degree centrality method and identified ten biomarker hub proteins (IL1B, CXCL8, FN1, MMP9, CXCL10, IL1A, IRF7, VWF, CXCL9, and ISG15) that associate COVID-19 with the comorbidities. Finally, we validated our findings by searching the published literature. Thus, our analytical approach elicited interconnections between COVID-19 and the aforementioned comorbidities in terms of remarkable DEGs, pathways, ontologies, PPI, and biomarker hub proteins.

Lethality risk markers by sex and age-group for COVID-19 in Mexico: a cross-sectional study based on machine learning approach

BACKGROUND

Mexico ranks fifth worldwide in the number of deaths due to COVID-19. Identifying risk markers through easily accessible clinical data could help in the initial triage of COVID-19 patients and anticipate a fatal outcome, especially in the most socioeconomically disadvantaged regions. This study aims to identify markers that increase lethality risk in patients diagnosed with COVID-19, based on machine learning (ML) methods. Markers were differentiated by sex and age-group.

METHODS

A total of 11,564 cases of COVID-19 in Mexico were extracted from the Epidemiological Surveillance System for Viral Respiratory Disease. Four ML classification methods were trained to predict lethality, and an interpretability approach was used to identify those markers.

RESULTS

Models based on Extreme Gradient Boosting (XGBoost) yielded the best performance in a test set. This model achieved a sensitivity of 0.91, a specificity of 0.69, a positive predictive value of 0.344, and a negative predictive value of 0.965. For female patients, the leading markers are diabetes and arthralgia. For males, the main markers are chronic kidney disease (CKD) and chest pain. Dyspnea, hypertension, and polypnea increased the risk of death in both sexes.

CONCLUSIONS

ML-based models using an interpretability approach successfully identified risk markers for lethality by sex and age. Our results indicate that age is the strongest demographic factor for a fatal outcome, while all other markers were consistent with previous clinical trials conducted in a Mexican population. The markers identified here could be used as an initial triage, especially in geographic areas with limited resources.

Prediction of mortality risk and duration of hospitalization of COVID-19 patients with chronic comorbidities based on machine learning algorithms

Background

The severity of coronavirus (COVID-19) in patients with chronic comorbidities is much higher than in other patients, which can lead to their death. Machine learning (ML) algorithms as a potential solution for rapid and early clinical evaluation of the severity of the disease can help in allocating and prioritizing resources to reduce mortality.

Objective

The objective of this study was to predict the mortality risk and length of stay (LoS) of patients with COVID-19 and history of chronic comorbidities using ML algorithms.

Methods

This retrospective study was conducted by reviewing the medical records of COVID-19 patients with a history of chronic comorbidities from March 2020 to January 2021 in Afzalipour Hospital in Kerman, Iran. The outcome of patients, hospitalization was recorded as discharge or death. The filtering technique used to score the features and well-known ML algorithms were applied to predict the risk of mortality and LoS of patients. Ensemble Learning methods is also used. To evaluate the performance of the models, different measures including F1, precision, recall, and accuracy were calculated. The TRIPOD guideline assessed transparent reporting.

Results

This study was performed on 1291 patients, including 900 alive and 391 dead patients. Shortness of breath (53.6%), fever (30.1%), and cough (25.3%) were the three most common symptoms in patients. Diabetes mellitus(DM) (31.3%), hypertension (HTN) (27.3%), and ischemic heart disease (IHD) (14.2%) were the three most common chronic comorbidities of patients. Twenty-six important factors were extracted from each patient's record. Gradient boosting model with 84.15% accuracy was the best model for predicting mortality risk and multilayer perceptron (MLP) with rectified linear unit function (MSE = 38.96) was the best model for predicting the LoS. The most common chronic comorbidities among these patients were DM (31.3%), HTN (27.3%), and IHD (14.2%). The most important factors in predicting the risk of mortality were hyperlipidemia, diabetes, asthma, and cancer, and in predicting LoS was shortness of breath.

Conclusion

The results of this study showed that the use of ML algorithms can be a good tool to predict the risk of mortality and LoS of patients with COVID-19 and chronic comorbidities based on physiological conditions, symptoms, and demographic information of patients. The Gradient boosting and MLP algorithms can quickly identify patients at risk of death or long-term hospitalization and notify physicians to do appropriate interventions.

Network-Based Data Analysis Reveals Ion Channel-Related Gene Features in COVID-19: A Bioinformatic Approach

Coronavirus disease 2019 (COVID-19) seriously threatens human health and has been disseminated worldwide. Although there are several treatments for COVID-19, its control is currently suboptimal. Therefore, the development of novel strategies to treat COVID-19 is necessary. Ion channels are located on the membranes of all excitable cells and many intracellular organelles and are key components involved in various biological processes. They are a target of interest when searching for drug targets. This study aimed to reveal the relevant molecular features of ion channel genes in COVID-19 based on bioinformatic analyses. The RNA-sequencing data of patients with COVID-19 and healthy subjects (GSE152418 and GSE171110 datasets) were obtained from the Gene Expression Omnibus (GEO) database. Ion channel genes were selected from the Hugo Gene Nomenclature Committee (HGNC) database. The RStudio software was used to process the data based on the corresponding R language package to identify ion channel-associated differentially expressed genes (DEGs). Based on the DEGs, Gene Ontology (GO) functional and pathway enrichment analyses were performed using the Enrichr web tool. The STRING database was used to generate a protein-protein interaction (PPI) network, and the Cytoscape software was used to screen for hub genes in the PPI network based on the cytoHubba plug-in. Transcription factors (TF)-DEG, DEG-microRNA (miRNA) and DEG-disease association networks were constructed using the NetworkAnalyst web tool. Finally, the screened hub genes as drug targets were subjected to enrichment analysis based on the DSigDB using the Enrichr web tool to identify potential therapeutic agents for COVID-19. A total of 29 ion channel-associated DEGs were identified. GO functional analysis showed that the DEGs were integral components of the plasma membrane and were mainly involved in inorganic cation transmembrane transport and ion channel activity functions. Pathway analysis showed that the DEGs were mainly involved in nicotine addiction, calcium regulation in the cardiac cell and neuronal system pathways. The top 10 hub genes screened based on the PPI network included KCNA2, KCNJ4, CACNA1A, CACNA1E, NALCN, KCNA5, CACNA2D1, TRPC1, TRPM3 and KCNN3. The TF-DEG and DEG-miRNA networks revealed significant TFs (FOXC1, GATA2, HINFP, USF2, JUN and NFKB1) and miRNAs (hsa-mir-146a-5p, hsa-mir-27a-3p, hsa-mir-335-5p, hsa-let-7b-5p and hsa-mir-129-2-3p). Gene-disease association network analysis revealed that the DEGs were closely associated with intellectual disability and cerebellar ataxia. Drug-target enrichment analysis showed that the relevant drugs targeting the hub genes CACNA2D1, CACNA1A, CACNA1E, KCNA2 and KCNA5 were gabapentin, gabapentin enacarbil, pregabalin, guanidine hydrochloride and 4-aminopyridine. The results of this study provide a valuable basis for exploring the mechanisms of ion channel genes in COVID-19 and clues for developing therapeutic strategies for COVID-19.

A Machine Learning Model for Predicting Hospitalization in Patients with Respiratory Symptoms during the COVID-19 Pandemic

A machine learning approach is a useful tool for risk-stratifying patients with respiratory symptoms during the COVID-19 pandemic, as it is still evolving. We aimed to verify the predictive capacity of a gradient boosting decision trees (XGboost) algorithm to select the most important predictors including clinical and demographic parameters in patients who sought medical support due to respiratory signs and symptoms (RAPID RISK COVID-19). A total of 7336 patients were enrolled in the study, including 6596 patients that did not require hospitalization and 740 that required hospitalization. We identified that patients with respiratory signs and symptoms, in particular, lower oxyhemoglobin saturation by pulse oximetry (SpO₂) and higher respiratory rate, fever, higher heart rate, and lower levels of blood pressure, associated with age, male sex, and the underlying conditions of diabetes mellitus and hypertension, required hospitalization more often. The predictive model yielded a ROC curve with an area under the curve (AUC) of 0.9181 (95% CI, 0.9001 to 0.9361). In conclusion, our model had a high discriminatory value which enabled the identification of a clinical and demographic profile predictive, preventive, and personalized of COVID-19 severity symptoms.

Deep-Risk: Deep Learning-Based Mortality Risk Predictive Models for COVID-19

The SARS-CoV-2 virus has proliferated around the world and caused panic to all people as it claimed many lives. Since COVID-19 is highly contagious and spreads quickly, an early diagnosis is essential. Identifying the COVID-19 patients' mortality risk factors is essential for reducing this risk among infected individuals. For the timely examination of large datasets, new computing approaches must be created. Many machine learning (ML) techniques have been developed to predict the mortality risk factors and severity for COVID-19 patients. Contrary to expectations, deep learning approaches as well as ML algorithms have not been widely applied in predicting the mortality and severity from COVID-19. Furthermore, the accuracy achieved by ML algorithms is less than the anticipated values. In this work, three supervised deep learning predictive models are utilized to predict the mortality risk and severity for COVID-19 patients. The first one, which we refer to as CV-CNN, is built using a convolutional neural network (CNN); it is trained using a clinical dataset of 12,020 patients and is based on the 10-fold cross-validation (CV) approach for training and validation. The second predictive model, which we refer to as CV-LSTM + CNN, is developed by combining the long short-term memory (LSTM) approach with a CNN model. It is also trained using the clinical dataset based on the 10-fold CV approach for training and validation. The first two predictive models use the clinical dataset in its original CSV form. The last one, which we refer to as IMG-CNN, is a CNN model and is trained alternatively using the converted images of the clinical dataset, where each image corresponds to a data row from the original clinical dataset. The experimental results revealed that the IMG-CNN predictive model outperforms the other two with an average accuracy of 94.14%, a precision of 100%, a recall of 91.0%, a specificity of 100%, an F1-score of 95.3%, an AUC of 93.6%, and a loss of 0.22.

Mortality Prediction Analysis among COVID-19 Inpatients Using Clinical Variables and Deep Learning Chest Radiography Imaging Features

The emergence of the COVID-19 pandemic over a relatively brief interval illustrates the need for rapid data-driven approaches to facilitate clinical decision making. We examined a machine learning process to predict inpatient mortality among COVID-19 patients using clinical and chest radiographic data. Modeling was performed with a de-identified dataset of encounters prior to widespread vaccine availability. Non-imaging predictors included demographics, pre-admission clinical history, and past medical history variables. Imaging features were extracted from chest radiographs by applying a deep convolutional neural network with transfer learning. A multi-layer perceptron combining 64 deep learning features from chest radiographs with 98 patient clinical features was trained to predict mortality. The Local Interpretable Model-Agnostic Explanations (LIME) method was used to explain model predictions. Non-imaging data alone predicted mortality with an ROC-AUC of 0.87 ± 0.03 (mean ± SD), while the addition of imaging data improved prediction slightly (ROC-AUC: 0.91 ± 0.02). The application of LIME to the combined imaging and clinical model found HbA1c values to contribute the most to model prediction (17.1 ± 1.7%), while imaging contributed 8.8 ± 2.8%. Age, gender, and BMI contributed 8.7%, 8.2%, and 7.1%, respectively. Our findings demonstrate a viable explainable AI approach to quantify the contributions of imaging and clinical data to COVID mortality predictions.

ISW-LM: An intensive symptom weight learning mechanism for early COVID-19 diagnosis

The novel coronavirus disease 2019 (COVID-19) pandemic has severely impacted the world. The early diagnosis of COVID-19 and self-isolation can help curb the spread of the virus. Besides, a simple and accurate diagnostic method can help in making rapid decisions for the treatment and isolation of patients. The analysis of patient characteristics, case trajectory, comorbidities, symptoms, diagnosis, and outcomes will be performed in the model. In this paper, a symptom-based machine learning (ML) model with a new learning mechanism called Intensive Symptom Weight Learning Mechanism (ISW-LM) is proposed. The proposed model designs three new symptoms' weight functions to identify the most relevant symptoms used to diagnose and classify COVID-19. To verify the efficiency of the proposed model, multiple laboratory and clinical datasets containing epidemiological symptoms and blood tests are used. Experiments indicate that the importance of COVID-19 infection symptoms varies between countries and regions. In most datasets, the most frequent and significant predictive symptoms for diagnosing COVID-19 are fever, sore throat, and cough. The experiment also compares the state-of-the-art methods with the proposed method, which shows that the proposed model has a high accuracy rate of up to 97.1711%. The positive results indicate that the proposed learning mechanism can help clinicians quickly diagnose and screen patients for COVID-19 at an early stage.

The Renin-Angiotensin System as a Component of Biotrauma in Acute Respiratory Distress Syndrome

Acute respiratory distress syndrome (ARDS) is a major concern in critical care medicine with a high mortality of over 30%. Injury to the lungs is caused not only by underlying pathological conditions such as pneumonia, sepsis, or trauma, but also by ventilator-induced lung injury (VILI) resulting from high positive pressure levels and a high inspiratory oxygen fraction. Apart from mechanical factors that stress the lungs with a specific physical power and cause volutrauma and barotrauma, it is increasingly recognized that lung injury is further aggravated by biological mediators. The COVID-19 pandemic has led to increased interest in the role of the renin-angiotensin system (RAS) in the context of ARDS, as the RAS enzyme angiotensin-converting enzyme 2 serves as the primary cell entry receptor for severe acute respiratory syndrome (SARS) coronavirus (CoV)-2. Even before this pandemic, studies have documented the involvement of the RAS in VILI and its dysregulation in clinical ARDS. In recent years, analytical tools for RAS investigation have made major advances based on the optimized precision and detail of mass spectrometry. Given that many clinical trials with pharmacological interventions in ARDS were negative, RAS-modifying drugs may represent an interesting starting point for novel therapeutic approaches. Results from animal models have highlighted the potential of RAS-modifying drugs to prevent VILI or treat ARDS. While these drugs have beneficial pulmonary effects, the best targets and application forms for intervention still have to be determined to avoid negative effects on the circulation in clinical settings.

Using Machine Learning to Predict Mortality for COVID-19 Patients on Day 0 in the ICU

Rationale: Given the expanding number of COVID-19 cases and the potential for new waves of infection, there is an urgent need for early prediction of the severity of the disease in intensive care unit (ICU) patients to optimize treatment strategies. Objectives: Early prediction of mortality using machine learning based on typical laboratory results and clinical data registered on the day of ICU admission. Methods: We retrospectively studied 797 patients diagnosed with COVID-19 in Iran and the United Kingdom (U.K.). To find parameters with the highest predictive values, Kolmogorov-Smirnov and Pearson chi-squared tests were used. Several machine learning algorithms, including Random Forest (RF), logistic regression, gradient boosting classifier, support vector machine classifier, and artificial neural network algorithms were utilized to build classification models. The impact of each marker on the RF model predictions was studied by implementing the local interpretable model-agnostic explanation technique (LIME-SP). Results: Among 66 documented parameters, 15 factors with the highest predictive values were identified as follows: gender, age, blood urea nitrogen (BUN), creatinine, international normalized ratio (INR), albumin, mean corpuscular volume (MCV), white blood cell count, segmented neutrophil count, lymphocyte count, red cell distribution width (RDW), and mean cell hemoglobin (MCH) along with a history of neurological, cardiovascular, and respiratory disorders. Our RF model can predict patient outcomes with a sensitivity of 70% and a specificity of 75%. The performance of the models was confirmed by blindly testing the models in an external dataset. Conclusions: Using two independent patient datasets, we designed a machine-learning-based model that could predict the risk of mortality from severe COVID-19 with high accuracy. The most decisive variables in our model were increased levels of BUN, lowered albumin levels, increased creatinine, INR, and RDW, along with gender and age. Considering the importance of early triage decisions, this model can be a useful tool in COVID-19 ICU decision-making.

A deep learning approach using effective preprocessing techniques to detect COVID-19 from chest CT-scan and X-ray images

Coronavirus disease-19 (COVID-19) is a severe respiratory viral disease first reported in late 2019 that has spread worldwide. Although some wealthy countries have made significant progress in detecting and containing this disease, most underdeveloped countries are still struggling to identify COVID-19 cases in large populations. With the rising number of COVID-19 cases, there are often insufficient COVID-19 diagnostic kits and related resources in such countries. However, other basic diagnostic resources often do exist, which motivated us to develop Deep Learning models to assist clinicians and radiologists to provide prompt diagnostic support to the patients. In this study, we have developed a deep learning-based COVID-19 case detection model trained with a dataset consisting of chest CT scans and X-ray images. A modified ResNet50V2 architecture was employed as deep learning architecture in the proposed model. The dataset utilized to train the model was collected from various publicly available sources and included four class labels: confirmed COVID-19, normal controls and confirmed viral and bacterial pneumonia cases. The aggregated dataset was preprocessed through a sharpening filter before feeding the dataset into the proposed model. This model attained an accuracy of 96.452% for four-class cases (COVID-19/Normal/Bacterial pneumonia/Viral pneumonia), 97.242% for three-class cases (COVID-19/Normal/Bacterial pneumonia) and 98.954% for two-class cases (COVID-19/Viral pneumonia) using chest X-ray images. The model acquired a comprehensive accuracy of 99.012% for three-class cases (COVID-19/Normal/Community-acquired pneumonia) and 99.99% for two-class cases (Normal/COVID-19) using CT-scan images of the chest. This high accuracy presents a new and potentially important resource to enable radiologists to identify and rapidly diagnose COVID-19 cases with only basic but widely available equipment.

Effects of Bacille Calmette Guerin (BCG) vaccination during COVID-19 infection

The coronavirus disease 2019 (COVID-19) is caused by the infection of highly contagious severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), also known as the novel coronavirus. In most countries, the containment of this virus spread is not controlled, which is driving the pandemic towards a more difficult phase. In this study, we investigated the impact of the Bacille Calmette Guerin (BCG) vaccination on the severity and mortality of COVID-19 by performing transcriptomic analyses of SARS-CoV-2 infected and BCG vaccinated samples in peripheral blood mononuclear cells (PBMC). A set of common differentially expressed genes (DEGs) were identified and seeded into their functional enrichment analyses via Gene Ontology (GO)-based functional terms and pre-annotated molecular pathways databases, and their Protein-Protein Interaction (PPI) network analysis. We further analysed the regulatory elements, possible comorbidities and putative drug candidates for COVID-19 patients who have not been BCG-vaccinated. Differential expression analyses of both BCG-vaccinated and COVID-19 infected samples identified 62 shared DEGs indicating their discordant expression pattern in their respected conditions compared to control. Next, PPI analysis of those DEGs revealed 10 hub genes, namely ITGB2, CXCL8, CXCL1, CCR2, IFNG, CCL4, PTGS2, ADORA3, TLR5 and CD33. Functional enrichment analyses found significantly enriched pathways/GO terms including cytokine activities, lysosome, IL-17 signalling pathway, TNF-signalling pathways. Moreover, a set of identified TFs, miRNAs and potential drug molecules were further investigated to assess their biological involvements in COVID-19 and their therapeutic possibilities. Findings showed significant genetic interactions between BCG vaccination and SARS-CoV-2 infection, suggesting an interesting prospect of the BCG vaccine in relation to the COVID-19 pandemic. We hope it may potentially trigger further research on this critical phenomenon to combat COVID-19 spread.