Next Article in Journal
Can Machine Learning Assist in Diagnosis of Primary Immune Thrombocytopenia? A Feasibility Study
Next Article in Special Issue
Applications of ChatGPT in Heart Failure Prevention, Diagnosis, Management, and Research: A Narrative Review
Previous Article in Journal
High-Resolution Imaging in Macular Telangiectasia Type 2: Case Series and Literature Review
Previous Article in Special Issue
Enhancing Clinical Validation for Early Cardiovascular Disease Prediction through Simulation, AI, and Web Technology
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Combining the Strengths of the Explainable Boosting Machine and Metabolomics Approaches for Biomarker Discovery in Acute Myocardial Infarction

by
Ahmet Kadir Arslan
1,
Fatma Hilal Yagin
1,*,
Abdulmohsen Algarni
2,
Fahaid AL-Hashem
3 and
Luca Paolo Ardigò
4,*
1
Department of Biostatistics and Medical Informatics, Faculty of Medicine, Inonu University, Malatya 44280, Türkiye
2
Department of Computer Science, King Khalid University, Abha 61421, Saudi Arabia
3
Department of Physiology, College of Medicine, King Khalid University, Abha 61421, Saudi Arabia
4
Department of Teacher Education, NLA University College, 0166 Oslo, Norway
*
Authors to whom correspondence should be addressed.
Diagnostics 2024, 14(13), 1353; https://doi.org/10.3390/diagnostics14131353
Submission received: 29 May 2024 / Revised: 20 June 2024 / Accepted: 24 June 2024 / Published: 26 June 2024
(This article belongs to the Special Issue Artificial Intelligence in Cardiology Diagnosis )

Abstract

:
Acute Myocardial Infarction (AMI), a common disease that can have serious consequences, occurs when myocardial blood flow stops due to occlusion of the coronary artery. Early and accurate prediction of AMI is critical for rapid prognosis and improved patient outcomes. Metabolomics, the study of small molecules within biological systems, is an effective tool used to discover biomarkers associated with many diseases. This study intended to construct a predictive model for AMI utilizing metabolomics data and an explainable machine learning approach called Explainable Boosting Machines (EBM). The EBM model was trained on a dataset of 102 prognostic metabolites gathered from 99 individuals, including 34 healthy controls and 65 AMI patients. After a comprehensive data preprocessing, 21 metabolites were determined as the candidate predictors to predict AMI. The EBM model displayed satisfactory performance in predicting AMI, with various classification performance metrics. The model’s predictions were based on the combined effects of individual metabolites and their interactions. In this context, the results obtained in two different EBM modeling, including both only individual metabolite features and their interaction effects, were discussed. The most important predictors included creatinine, nicotinamide, and isocitrate. These metabolites are involved in different biological activities, such as energy metabolism, DNA repair, and cellular signaling. The results demonstrate the potential of the combination of metabolomics and the EBM model in constructing reliable and interpretable prediction outputs for AMI. The discussed metabolite biomarkers may assist in early diagnosis, risk assessment, and personalized treatment methods for AMI patients. This study successfully developed a pipeline incorporating extensive data preprocessing and the EBM model to identify potential metabolite biomarkers for predicting AMI. The EBM model, with its ability to incorporate interaction terms, demonstrated satisfactory classification performance and revealed significant metabolite interactions that could be valuable in assessing AMI risk. However, the results obtained from this study should be validated with studies to be carried out in larger and well-defined samples.

1. Introduction

Acute myocardial infarction (AMI) stands as a critical juncture in cardiovascular health, characterized by the abrupt interruption of blood flow within the heart muscle [1]. This condition demands swift and precise intervention to avert severe consequences [2]. Comprising both ST-segment elevation myocardial infarction (STEMI) and non-ST-segment elevation myocardial infarction (NSTEMI), AMI subjects individuals to a heightened susceptibility to recurrent ischemic events [3]. The urgency of preventing AMI requires a comprehensive investigation of the complex pathways and contributing factors that lead to the risk.
In recent years, scientific studies have been carried out with “omics” approaches to identify biomarkers for various diseases. Metabolomics, one of these important approaches, involves the comprehensive analysis of small molecules called metabolites in cells, biological fluids, tissues, or organisms. Possible biomarkers to be determined by metabolomic analyses can be used in early diagnosis and follow-up of diseases and guidance of treatment strategies. Metabolomics is an area of omics sciences that focuses on the exhaustive study of small compounds or metabolites inside biological systems. The metabolites comprise diverse components such as carbohydrates, amino acids, lipids, and other chemical molecules. Metabolomics strives to comprehend the complex metabolic processes that occur within organisms and provides insights into their health, physiology, and responses to environmental factors [4,5].
Explainable Boosting Machines (EBM) is a sophisticated machine learning method that puts together the utilities of both boosting and interpretability. EBM is an interpretable machine learning model aimed to ensure both accurate predictions and insights into the model’s decision-making process and is particularly beneficial in cases where model interpretability is critical, such as credit scoring and healthcare. EBM is built on the concept of boosting, where numerous weak models are joined to form a powerful and useful predictive model. It consecutively fits a sequence of decision trees/weak learners, each one dedicated to correcting the errors committed by the former models [6,7].
The main objective of this study is to discover prognostic metabolite features that are candidates for predicting AMI using EBM. EBM was chosen in this research due to its speed and interpretability/explainability at local and global levels, as well as its advanced variable engineering features such as automatic binary interaction term detection. The secondary aim of the study is to investigate the contribution of binary interactions of metabolite biomarkers to the AMI classification performance and explainability of the model through automatic binary interaction term detection. For these purposes, a combined prognostic prediction model with the biomarker discovery process for AMI was developed. This framework enhances the precision of prognostic predictions for AMI and provides intuitive insights into these predictions through the consideration of metabolomics risk factors.

2. Materials and Methods

In this study, by following the MINIMAR (MINimum Information for Medical AI Reporting) [8], minimum reporting standards were tried to be ensured, especially in the model architecture and interpretation/evaluation of the results. The summary of both the coverage status of the MINIMAR and the key points of the Material and Methods section was provided in Supplementary Table S1.

2.1. The Data Set

This data is available at the NIH Common Fund’s National Metabolomics Data Repository (NMDR) website, the Metabolomics Workbench, https://www.metabolomicsworkbench.org (access date: 30 February 2024) where it has been assigned Project ID (ST001908) [9]. The dataset contains intensity values of 102 prognostic metabolites as continuous numeric type obtained by LC–MS (Liquid chromatography–mass spectrometry) technique and a binary outcome variable (Control/AMI). The list of 102 metabolites analyzed in the study and 21 metabolites included in the model training after the feature selection phase was provided in Supplementary Table S2. The current dataset consisted of 99 individuals, 34 healthy controls, and 65 individuals diagnosed with AMI. The gender distribution of the dataset is male (74, 74.7%) and female (25, 25.3%). According to the metadata in the Metabolomics Workbench database from which the dataset was downloaded, the study consisted of 100 individuals, but the downloaded dataset contained mass intensity values for 99 individuals.

2.2. Superficial Data Set Quality Check

The metabolite features were examined in terms of missingness and measures of dispersion by standard deviation, min–max scaling, and near-zero variance. No problems were encountered in the checks that could affect the next stages of data analysis.

2.2.1. Outlier Analysis Phase

In this study, the OutlierTree [10] Python library (also the name of the method) was applied to carry out the outlier detection process. This method that leverages decision trees offers researchers an explainable perspective on outlier detection. To discover outliers in a dataset by constructing a tree-based model specifically optimized for outlier detection is aimed. This method identifies outliers without assigning scores to individual observations. For this, it builds decision trees to understand how different columns in the data relate to each other. It then examines these relationships to determine observations that stand out as significantly different from the rest [10,11].

2.2.2. Missing Value Imputation Phase

To handle missing values, miceforest [12] Python library was employed. This tool aims to impute missing values using the LightGBM [13] model-based Multiple Imputation by Chained Equations (MICE) [14] approach. Using the MICE technique combined with the power of LightGBM, this approach provides computational speed, memory efficiency, and data type flexibility for missing value imputations [12].

2.2.3. Feature Selection (FS) Phase

In this study, the Boruta method [15] was carried out for the FS task. Boruta is a robust and data-driven FS technique that can help researchers identify the most relevant predictors in their datasets. It is particularly beneficial in circumstances where it is wished to automate the FS process and avoid manually establishing feature importance thresholds. This feature enabled the Boruta technique to be used in this study.

2.2.4. Model Training Phase

  • The employed model in a nutshell
EBM, founded in InterpretML [16] Python package and from the interpretable machine learning model family, is a model based on generalized additive models [17] with binary interaction terms.
EBM uses a set of interpretable predictors for each instance to generate predictions. Making the model more interpretable compared to sophisticated models like deep neural networks, predictors are split into bins. The primary strength of EBM is its transparency. It ensures explicit insights into how each feature affects the prediction of the model. To achieve this mission, EBM generates global and local explanations (overall predictor importance and how features affect individual predictions). EBM produces predictor importance scores, which express the impact of each predictor on the model’s predictions. These scores are used to rank and prioritize predictors based on their contributions. The complexity of the EBM can be regulated by adjusting the number of weak models (trees) in the ensemble. This can help trade-offs among interpretability and prediction strength [16].
For any ith data point in the dataset (xi, yi), the form of EBM is characterized as [16]:
g ( E [ y ] ) = β 0 + f i ( x i ) + f i , j ( x i , x j ) ,
where g(.), β0, fi, and fi,j represent link function, intercept, shape/smooth functions, and the interaction effect between features x i ,   and   x j , respectively. The functions are learned by EBM employing enhanced ensemble learning tools. EBM provides a rapid implementation of the GA2M (Generalized Additive Model with Pairwise Interactions) algorithm. The binary interaction terms in the formula are automatically added to the model using the FAST (Feature-based Algorithm for Subgroup Discovery) technique [18].
b.
The rest of the model training details
Before model training, the data set was randomly split into two parts for validating the model: training (75%) and test (25%) sets. To determine the value ranges of the optimization parameters of the EBM model, a comprehensive literature review was made, and the ranges specified in Table 1 were used in the model training process. In an attempt to find optimal hyperparameters, the grid search technique was used together with the 5-fold cross-validation method. The detailed information about the hyperparameters experimented during the training process is presented in Table 1.
After determining the optimal values for the hyperparameters, the model was retrained using the relevant values, and the prediction process was completed on the test data set.

2.2.5. Model Performance Evaluations

The binary classification performance of the EBM model was performed using accuracy, sensitivity, specificity, F1 measure, and area under the ROC (Receiver Operating Characteristics) curve. A graphical representation called the ROC curve is used to assess how well a binary classification model performs. At different threshold values, it compares the True Positive Rate (TPR) against the False Positive Rate (FPR). Additionally, 95% bootstrapped confidence intervals with 1000 repetitions were calculated.

3. Results

Table 2 presents the distribution statistics and p-values of AMI and Control groups in terms of gender, age, Body Mass Index (BMI), and smoking.
The superficial checks of the dataset showed there were no problems that could affect the data preprocessing and modeling steps. In the analyses conducted for outlier detection, a total of 13 suspected outliers were detected in different features. There were no missing values in the original data set. However, after the outlier analysis, 13 suspected outliers were removed from the dataset and replaced with “NA” (Not available). The missing values were imputed by the LightGBM-powered MICE method. In the feature selection analysis conducted after missing value imputation, the Boruta method selected 21 features out of 102 metabolite features, gender, age, body mass index (BMI), and smoking, and proceeded to the EBM modeling stage with these selected features. After parameter optimization, two EBM models were trained with and without the addition of interaction terms. The results of the confusion matrices and ROC plots of both models related to test classification performances are presented in Figure 1. Figure 1a shows that the EBM model with the interaction term misclassified only 2 instances and the area under the ROC curve of the model was 0.95. Figure 1b shows the performance outputs of the EBM model without the interaction term, which misclassified 4 instances and had an area under the ROC curve of 0.93. In addition, the measurements for both training and test classification performances of the related models by various metrics are given in Table 3 with 95% confidence intervals.
Upon evaluation of classification performances in Table 3, a notable sensitivity score was achieved in predicting AMI. A heightened sensitivity implies a reduced occurrence of false negatives (FN). In comparative biological studies, false positive and false negative errors are prevalent. Hence, establishing the probability of a genuine effect being significant holds paramount importance. A diminished FN value signifies a promising outcome for AMI instances. This outcome holds significant importance as the primary objective of this research is to minimize overlooked AMI cases (false negatives).
In Figure 2, the weighted mean absolute score-based feature importances for the two EBM models were demonstrated. The top three most important features determined by the EBM model without interaction terms were “Nicotinamide”, “Creatinine”, and “Isocitrate” (Figure 2).
Figure 3 shows the feature importance scores of the interaction terms added EBM model. When the feature importance outputs of the EBM model with interaction terms are examined, it is observed that the first 5 metabolites (Creatinine and nicotinamide, Creatinine and isocitrate, nicotinamide and UDP-D-glucose, Creatinine and Citraconic acid, nicotinamide) are mostly composed of the 3 metabolites mentioned above.
Figure 4 is shown divided into 4 rows. Row (a) shows the heat plots for the first two interaction terms with the highest weighted average absolute score according to the output of the EBM model with interaction terms. In this row, Nicotinamide x Creatinine and Isocitrate x Creatinine interactions are demonstrated, respectively. The pink-colored rectangles represent areas of increased risk of AMI. Row (b) shows the distribution of single-term metabolites for the same model and the effect of intensity variation on classification prediction and AMI risk. In the related graphs in Row (b), it can be said that the increase in creatinine and nicotinamide concentrations increases the risk of AMI. On the other hand, an increase in isocitrate concentration decreases the risk of AMI. In rows (c) and (d), boxplots and ROC plots from training data for creatinine, nicotinamide, and isocitrate metabolites are presented. Row (c) indicates that creatinine and nicotinamide concentration medians were higher and the isocitrate concentration median was lower in the AMI group. Row (d) shows that the creatinine metabolite feature has the highest ROC AUC value.
In Figure 5, the local explanations of two samples having the highest-class probabilities in terms of the true predicted Control (a) and AMI (b) groups. The orange bars indicate the terms, contributing to the AMI classification and the blue bars indicate the terms contributing to the control group classification. For the sample labeled “Control” by the model, nicotinamide is the term that contributes most to the prediction. Also, for the sample labeled as “AMI” by the model, the feature that contributes the most to the prediction is the creatinine and isocitrate interaction term. For the relevant sample, the interaction terms, both involving the creatinine metabolite, are contributing to the model’s prediction of AMI.

4. Discussion

In this study, a pipeline including extensive data preprocessing stages and EBM, an explainable machine learning approach, was constructed to identify potential metabolite biomarkers that can be used to predict AMI. Additionally, the classification performance of the EBM model on the preprocessed training data was compared with and without adding interaction terms.
A study [19], which produced the dataset used in this study, aimed to find even more in-depth and interesting patterns for predicting AMI using various omics approaches than single omics studies. In this study, Random Forest was used as the classification model, 27 metabolite features were included in the modeling and the classification performance of the model was found to be 0.836 in terms of AUC ROC. In this study, the ROC AUC values of the EBM models trained in two different scenarios with 21 metabolite features were 0.93 and 0.95. Although the data sets considered are the same, it would not be fair to compare the two models as they are not subjected to the same data preprocessing process. However, it can be argued that EBM is more convenient than Random Forest for relevant tasks in terms of being boosting-based, detecting binary interaction terms automatically, and producing locally/globally interpretable results.
In another LC-MS/MS-based metabolomics research [20], unstable angina and AMI, subtypes of acute coronary syndrome (ACS), and healthy control groups were extensively compared in terms of serum and urine metabolites. In this multi-classification task performing Multinominal Adaptive LASSO and Random Forest models, 2-ketobutyric acid, LysoPC (18:2(9Z,12Z)), argininosuccinic acid, and cyclic GMP metabolites were identified as possible biomarkers for ACS prediction.
It is observed in Table 3 that the interaction term added EBM model outperforms the single term EBM model in terms of all classification performance metrics. The inclusion of automatically detected binary interaction terms in the EBM model is intended to improve model performance at the same time without compromising explainability, as stated by the developers of the model [14]. This is important because, in general, model prediction performance and explainability/interpretability are assumed as negatively correlated [21]. It can be stated that this is also confirmed in the current study when the classification metrics and unique explainable items produced by the models are considered.
When the feature importance plots obtained from both EBM models are examined, it is pointed out that the single and interaction effects of creatinine, nicotinamide, and isocitrate metabolites stand out in terms of contribution to the AMI classification performance. The discussion of the results will be based on these 3 metabolites.
Considering Figure 2 and Figure 3 together, it is observed that nicotinamide is the metabolite with the highest feature importance in EBM modeling without adding the interaction terms. In the EBM model with added interaction terms, it was observed that the metabolite with the fourth-highest feature importance (highest among single effects). Furthermore, the AUC value (0.83) obtained from the model-independent ROC analysis for nicotinamide has a good degree of discrimination (Figure 4c). Nicotinamide, also known as niacinamide, is a form of vitamin B3 and a versatile compound with many important functions such as energy metabolism, DNA repair, cellular signaling, and gene expression in the body [22,23,24,25]. Nicotinamide is a precursor to nicotinamide adenine dinucleotide (NAD+), a crucial coenzyme involved in various cellular processes, including energy production, DNA repair, and cell signaling. During AMI, the heart muscle experiences significant energy depletion. Nicotinamide’s potential to boost NAD+ levels could play a role in restoring energy balance and mitigating cellular damage. Nicotinamide exhibits antioxidant properties that might protect cardiac cells from oxidative stress, a major contributor to tissue damage during AMI [26,27,28].
In addition, the relationship between nicotinamide metabolite and the smoking status of the participants in this study was also examined. A total of 54 individuals in the dataset were current smokers, while the remaining 45 were either never smokers or former smokers. When these two groups were compared in terms of normalized nicotinamide concentration, a statistically significant difference was found between the groups (p = 0.002, t = 3.2, df = 97). Furthermore, the mean normalized nicotinamide concentrations were lower in the smoker group (Figure 4c). Also, when the percentages were analyzed, it was observed that the percentage of smokers in the AMI group (40.7%) was lower than in the Control group (59.3%) (Table 1). This is confirmed by the information that nicotinamide concentration is lower in smokers [29]. Perhaps this may suggest whether smoking status has a confounding effect on the determination of nicotinamide as a biomarker.
According to the results, the creatinine metabolite has the second highest feature importance in the single-effect EBM model. If the model with interaction terms is considered, the creatinine effect is present in 3 of the first 4 interaction terms that contribute the most to the classification performance. Similar to nicotinamide, it is observed that creatinine levels above certain normalized intensity values may increase the risk of AMI (Figure 4a,b). In addition, when the ROC analysis outputs were analyzed, the metabolite with the highest AUC value among the 3 metabolites was creatinine Figure 4b. The boxplots show that the median of nicotinamide and creatinine mass density values were higher in the AMI group (Figure 4c).
Creatinine is a waste product from the breakdown of creatinine phosphate, an amino acid that muscles use to produce energy [30]. Elevated levels of AMI may indicate damage to the heart muscle as a direct consequence of the infarction (heart attack). Creatinine is also an important indicator of kidney function. During AMI, decreased blood flow to the kidneys can impair their function, leading to reduced creatinine clearance and subsequently elevated serum levels. Creatine, the precursor to creatinine, plays a vital role in cellular energy metabolism. Disrupted energy production during AMI could contribute to alterations in creatinine levels [31,32,33]. The studies investigating the utility of creatinine as a risk factor for heart disease have shown that people with high levels of creatinine metabolites have a higher risk of developing heart disease [34,35]. In the case of AMI, an observational study in 2019 investigated the metabolic features of CHD patients using a targeted metabolomics approach. Blood samples were collected from 302 patients with CHD and 59 normal coronary artery (NCA) subjects and analyzed using the LC-MS technique. According to the results, it was determined that creatinine concentration was significantly higher in the CHD group than in the NCA group, including patients with AMI [36].
Furthermore, in the related study [37] that produced the dataset used in this study, both creatinine and nicotinamide metabolite levels were found to be up-regulated. This indicates that the analyte concentrations of these two biomarkers are higher in patients with AMI than in the control group.
In this study, isocitrate was observed as the 3rd most important metabolite feature. Isocitrate is a molecule involved in a series of chemical reactions known as the Krebs cycle (TCA cycle or citric acid cycle). The Krebs cycle is a series of reactions that cells use to produce energy [38,39]. Isocitrate is an intermediate in the tricarboxylic acid (TCA) cycle, the central pathway for cellular energy production. Disruptions to the TCA cycle during AMI can impact isocitrate levels. Elevated isocitrate might indicate impaired energy metabolism in the heart muscle. The TCA cycle occurs within mitochondria, the cellular powerhouses. Alterations in isocitrate levels could reflect mitochondrial dysfunction, a hallmark of AMI [40,41].
However, the study [37] that provided the data set used in the current study found that the relevant metabolite was down-regulated in AMI patients compared to the control group, indicating a lower intensity of the analyte. This output is in line with the results of the current study.

5. Conclusions

In this study, to determine the effect of potential biomarkers that can be evaluated to predict AMI and their interactions on AMI prediction was aimed. For this purpose, two different variants of the EBM model were trained on the training dataset and these two cases were compared with each other using various classification metrics and feature importance values. The EBM model with interaction terms provided satisfactory classification performance and identified metabolite interactions that can be considered in assessing the risk of AMI. Furthermore, the fact that it contributes to the examination of the metabolite contribution of AMI risk on a sample-based and that it generates these findings on its own and not using external model agnostic methods such as SHAP (Shapley Additive explanations), etc., makes the EBM model attractive for use as an explainable/interpretable model in metabolomics research. In conclusion, this study successfully developed a pipeline incorporating extensive data preprocessing and the EBM model to identify potential metabolite biomarkers for predicting AMI. The EBM model, with its ability to incorporate interaction terms, demonstrated satisfactory classification performance and revealed significant metabolite interactions that could be valuable in assessing AMI risk.

6. Limitations and Future Works

The two main limitations of this study are the relatively small sample size and the lack of an external validation process that would enable a more robust assessment of the classification performance and other outputs of the model. The application of the model in a multi-omics study in combination with other omics approaches to identify AMI risk factors and determine their contribution more robustly is suggested as future research.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/diagnostics14131353/s1, Table S1: Highlights for Material and Methods and coverage status of the MINIMAR; Table S2: The list of 102 metabolites analyzed in the study and 21 metabolites included in the model training after the feature selection phase.

Author Contributions

Conceptualization, A.K.A. and F.H.Y.; Data curation, A.K.A. and F.H.Y.; Formal analysis, A.K.A.; Investigation, A.K.A. and F.H.Y.; Methodology, A.K.A. and F.H.Y.; Resources, A.K.A. and F.H.Y.; Software, A.K.A.; Validation, A.K.A.; Writing—original draft, A.K.A., F.H.Y., A.A., F.A.-H. and L.P.A.; Writing—review and editing, A.K.A., F.H.Y., A.A., F.A.-H. and L.P.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research was financially supported by the Deanship of Scientific Research at King Khalid University under research grant number (R.G.P.2/93/45).

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki and was approved by the Inonu University Health Sciences Non-Interventional Clinical Research Ethics Committee (protocol code = 2024/5692, approval date 19 March 2024).

Informed Consent Statement

Not applicable.

Data Availability Statement

The data is open access and can be requested from the corresponding author upon appropriate request.

Acknowledgments

This research was financially supported by the Deanship of Scientific Research at King Khalid University under research grant number (R.G.P.2/93/45).

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Reed, G.W.; Rossi, J.E.; Cannon, C.P. Acute myocardial infarction. Lancet 2017, 389, 197–210. [Google Scholar] [CrossRef]
  2. Hajar, R. Evolution of myocardial infarction and its biomarkers: A historical perspective. Heart Views Off. J. Gulf Heart Assoc. 2016, 17, 167. [Google Scholar] [CrossRef] [PubMed]
  3. Kim, Y.H.; Her, A.-Y.; Rha, S.-W.; Choi, C.U.; Choi, B.G.; Kim, J.B.; Park, S.; Kang, D.O.; Park, J.Y.; Park, S.-H. Comparison of clinical outcomes after non-ST-segment and ST-segment elevation myocardial infarction in diabetic and nondiabetic populations. J. Clin. Med. 2022, 11, 5079. [Google Scholar] [CrossRef]
  4. Clish, C.B. Metabolomics: An emerging but powerful tool for precision medicine. Cold Spring Harb. Mol. Case Stud. 2015, 1, a000588. [Google Scholar] [CrossRef] [PubMed]
  5. Johnson, C.H.; Ivanisevic, J.; Siuzdak, G. Metabolomics: Beyond biomarkers and towards mechanisms. Nat. Rev. Mol. Cell Biol. 2016, 17, 451–459. [Google Scholar] [CrossRef]
  6. Guldogan, E.; Yagin, F.H.; Pinar, A.; Colak, C.; Kadry, S.; Kim, J. A proposed tree-based explainable artificial intelligence approach for the prediction of angina pectoris. Sci. Rep. 2023, 13, 22189. [Google Scholar] [CrossRef]
  7. Konstantinov, A.V.; Utkin, L.V. Interpretable machine learning with an ensemble of gradient boosting machines. Knowl.-Based Syst. 2021, 222, 1–28. [Google Scholar] [CrossRef]
  8. Hernandez-Boussard, T.; Bozkurt, S.; Ioannidis, J.P.A.; Shah, N.H. MINIMAR (MINimum Information for Medical AI Reporting): Developing reporting standards for artificial intelligence in health care. J. Am. Med. Inform. Assoc. 2020, 27, 2011–2015. [Google Scholar] [CrossRef]
  9. The Metabolomics Workbench. Available online: https://www.metabolomicsworkbench.org/ (accessed on 30 February 2024).
  10. Cortes, D. Explainable outlier detection through decision tree conditioning. arXiv 2020, arXiv:2001.00636. [Google Scholar]
  11. Zenkl-Galaz, M.A.; Loyola-González, O.; Medina-Pérez, M.A. IOGOD: An interpretable outlier generation-based outlier detector for categorical databases. Expert Syst. Appl. 2022, 195, 116570. [Google Scholar] [CrossRef]
  12. Wilson, S. miceforest: Fast, Memory Efficient Imputation with LightGBM. 2022. Available online: https://github.com/AnotherSamWilson/miceforest (accessed on 4 March 2024).
  13. Ke, G.; Meng, Q.; Finley, T.; Wang, T.; Chen, W.; Ma, W.; Ye, Q.; Liu, T.-Y. Lightgbm: A highly efficient gradient boosting decision tree. Adv. Neural Inf. Process. Syst. 2017, 30, 1–9. [Google Scholar]
  14. Van Buuren, S.; Groothuis-Oudshoorn, K. mice: Multivariate imputation by chained equations in R. J. Stat. Softw. 2011, 45, 1–67. [Google Scholar] [CrossRef]
  15. Kursa, M.B.; Jankowski, A.; Rudnicki, W.R. Boruta–a system for feature selection. Fundam. Informaticae 2010, 101, 271–285. [Google Scholar] [CrossRef]
  16. Nori, H.; Jenkins, S.; Koch, P.; Caruana, R. Interpretml: A unified framework for machine learning interpretability. arXiv 2019, arXiv:1909.09223. [Google Scholar]
  17. Hastie, T.; Tibshirani, R. Generalized Additive Models; Some Applications. In Generalized Linear Models. Lecture Notes in Statistics; Gilchrist, R., Francis, B., Whittaker, J., Eds.; Springer: New York, NY, USA, 1985; Volume 32. [Google Scholar]
  18. Lou, Y.; Caruana, R.; Gehrke, J.; Hooker, G. Accurate intelligible models with pairwise interactions. In Proceedings of the 19th ACM SIGKDD International Conference on Knowledge Discovery and Data Mining, Chicago, IL, USA, 11–14 August 2013; pp. 623–631. [Google Scholar]
  19. Lim, S.Y.; Lim, F.L.S.; Criado-Navarro, I.; Yeo, X.H.; Dayal, H.; Vemulapalli, S.D.; Seah, S.J.; Laserna, A.K.C.; Yang, X.; Tan, S.H. Multi-omics investigation into acute myocardial infarction: An integrative method revealing interconnections amongst the metabolome, lipidome, glycome, and metallome. Metabolites 2022, 12, 1080. [Google Scholar] [CrossRef] [PubMed]
  20. Fu, M.; He, R.; Zhang, Z.; Ma, F.; Shen, L.; Zhang, Y.; Duan, M.; Zhang, Y.; Wang, Y.; Zhu, L. Multinomial machine learning identifies independent biomarkers by integrated metabolic analysis of acute coronary syndrome. Sci. Rep. 2023, 13, 20535. [Google Scholar]
  21. Herm, L.-V.; Wanner, J.; Seubert, F.; Janiesch, C. I Don’t Get IT, but IT seems Valid! The Connection between Explainability and Comprehensibility in (X) AI Research. 2021. Available online: https://www.researchgate.net/profile/Christian_Janiesch/publication/351157118_I_Don't_Get_It_but_It_Seems_Valid_The_Connection_Between_Explainability_and_Comprehensibility_in_XAI_Research/links/608acb9ea6fdccaebdf89b4d/I-Dont-Get-It-but-It-Seems-Valid-The-Connection-Between-Explainability-and-Comprehensibility-in-XAI-Research.pdf (accessed on 25 March 2024).
  22. Makarov, M.V.; Trammell, S.A.; Migaud, M.E. The chemistry of the vitamin B3 metabolome. Biochem. Soc. Trans. 2019, 47, 131–147. [Google Scholar] [PubMed]
  23. Surjana, D.; Halliday, G.M.; Damian, D.L. Role of nicotinamide in DNA damage, mutagenesis, and DNA repair. J. Nucleic Acids 2010, 2010, 157591. [Google Scholar] [CrossRef]
  24. Anderson, G.D.; Peterson, T.C.; Farin, F.M.; Beyer, R.P.; Hoane, M. The effect of nicotinamide on gene expression in a traumatic brain injury model. Front. Neurosci. 2013, 7, 40698. [Google Scholar]
  25. Xie, N.; Zhang, L.; Gao, W.; Huang, C.; Huber, P.E.; Zhou, X.; Li, C.; Shen, G.; Zou, B. NAD+ metabolism: Pathophysiologic mechanisms and therapeutic potential. Signal Transduct. Target. Ther. 2020, 5, 227. [Google Scholar]
  26. Yao, W.; Pei, Z.; Zhang, X. NAD(+): A key metabolic regulator with great therapeutic potential for myocardial infarction via Sirtuins family. Heliyon 2023, 9, e21890. [Google Scholar] [CrossRef]
  27. Fontecha-Barriuso, M.; Lopez-Diaz, A.M.; Carriazo, S.; Ortiz, A.; Sanz, A.B. Nicotinamide and acute kidney injury. Clin. Kidney J. 2021, 14, 2453–2462. [Google Scholar] [CrossRef]
  28. Abdellatif, M.; Trummer-Herbst, V.; Koser, F.; Durand, S.; Adão, R.; Vasques-Nóvoa, F.; Freundt, J.K.; Voglhuber, J.; Pricolo, M.-R.; Kasa, M.; et al. Nicotinamide for the treatment of heart failure with preserved ejection fraction. Sci. Transl. Med. 2021, 13, eabd7064. [Google Scholar] [CrossRef]
  29. Liou, S. Nicotinamide. Available online: http://161.35.229.200/index.php/2010/06/29/nicotinamide/ (accessed on 9 March 2024).
  30. Williamson, L.; New, D. How the use of creatine supplements can elevate serum creatinine in the absence of underlying kidney pathology. Case Rep. 2014, 2014, bcr2014204754. [Google Scholar] [CrossRef]
  31. Babu, P.J.; Tirkey, A.; Rao, T.J.M.; Chanu, N.B.; Lalchhandama, K.; Singh, Y.D. Conventional and nanotechnology based sensors for creatinine (A kidney biomarker) detection: A consolidated review. Anal. Biochem. 2022, 645, 114622. [Google Scholar]
  32. Jose, P.; Skali, H.; Anavekar, N.; Tomson, C.; Krumholz, H.M.; Rouleau, J.L.; Moye, L.; Pfeffer, M.A.; Solomon, S.D. Increase in creatinine and cardiovascular risk in patients with systolic dysfunction after myocardial infarction. J. Am. Soc. Nephrol. JASN 2006, 17, 2886–2891. [Google Scholar] [CrossRef]
  33. Gibson, C.M.; Pinto, D.S.; Murphy, S.A.; Morrow, D.A.; Hobbach, H.P.; Wiviott, S.D.; Giugliano, R.P.; Cannon, C.P.; Antman, E.M.; Braunwald, E. Association of creatinine and creatinine clearance on presentation in acute myocardial infarction with subsequent mortality. J. Am. Coll. Cardiol. 2003, 42, 1535–1543. [Google Scholar] [CrossRef]
  34. Chen, X.; Jin, H.; Wang, D.; Liu, J.; Qin, Y.; Zhang, Y.; Zhang, Y.; Xiang, Q. Serum creatinine levels, traditional cardiovascular risk factors and 10-year cardiovascular risk in Chinese patients with hypertension. Front. Endocrinol. 2023, 14, 1140093. [Google Scholar]
  35. Vikse, B.; Vollset, S.; Tell, G.; Refsum, H.; Iversen, B. Distribution and determinants of serum creatinine in the general population: The Hordaland Health Study. Scand. J. Clin. Lab. Investig. 2004, 64, 709–722. [Google Scholar] [CrossRef]
  36. Zhong, Z.; Liu, J.; Zhang, Q.; Zhong, W.; Li, B.; Li, C.; Liu, Z.; Yang, M.; Zhao, P. Targeted metabolomic analysis of plasma metabolites in patients with coronary heart disease in southern China. Medicine 2019, 98, e14309. [Google Scholar] [CrossRef]
  37. Lim, S.Y.; Ng, B.H.; Vermulapalli, D.; Lau, H.; Carrasco Laserna, A.K.; Yang, X.; Tan, S.H.; Chan, M.Y.; Li, S.F.Y. Simultaneous Polar Metabolite and N-Glycan Extraction Workflow for Joint-Omics Analysis: A Synergistic Approach for Novel Insights into Diseases. J. Proteome Res. 2022, 21, 643–653. [Google Scholar] [CrossRef]
  38. Gasmi, A.; Peana, M.; Arshad, M.; Butnariu, M.; Menzel, A.; Bjørklund, G. Krebs cycle: Activators, inhibitors and their roles in the modulation of carcinogenesis. Arch. Toxicol. 2021, 95, 1161–1178. [Google Scholar] [CrossRef]
  39. Williams, N.C.; O’Neill, L.A. A role for the Krebs cycle intermediate citrate in metabolic reprogramming in innate immunity and inflammation. Front. Immunol. 2018, 9, 317864. [Google Scholar] [CrossRef]
  40. Jiang, M.; Xie, X.; Cao, F.; Wang, Y. Mitochondrial metabolism in myocardial remodeling and mechanical unloading: Implications for ischemic heart disease. Front. Cardiovasc. Med. 2021, 8, 789267. [Google Scholar] [CrossRef]
  41. Watany, M.M.; Abd-Ellatif, R.N.; Abdeldayem, M.E.; El-Horany, H.E.-s. Association between genetic variations of mitochondrial isocitrate dehydrogenase (IDH2) and acute myocardial infarction. Gene 2022, 829, 146497. [Google Scholar] [CrossRef]
Figure 1. Confusion matrices and ROC curves for trained Explainable Boosting Machine Model Test Set predictions: There are two types of performance: (a) without interaction terms and (b) with interaction terms.
Figure 1. Confusion matrices and ROC curves for trained Explainable Boosting Machine Model Test Set predictions: There are two types of performance: (a) without interaction terms and (b) with interaction terms.
Diagnostics 14 01353 g001
Figure 2. The feature importances of metabolites generated by the EBM model without adding interaction terms.
Figure 2. The feature importances of metabolites generated by the EBM model without adding interaction terms.
Diagnostics 14 01353 g002
Figure 3. The feature importances of metabolites generated by the EBM model were trained by adding interaction terms.
Figure 3. The feature importances of metabolites generated by the EBM model were trained by adding interaction terms.
Diagnostics 14 01353 g003
Figure 4. Interaction (a), distribution line (b), boxplot (c), and ROC (d) plots for the top 3 metabolites contributing the most to the model according to both the EBM models.
Figure 4. Interaction (a), distribution line (b), boxplot (c), and ROC (d) plots for the top 3 metabolites contributing the most to the model according to both the EBM models.
Diagnostics 14 01353 g004
Figure 5. The local explanations of two samples that have the highest-class probabilities in terms of the Control (a) and AMI (b) groups.
Figure 5. The local explanations of two samples that have the highest-class probabilities in terms of the Control (a) and AMI (b) groups.
Diagnostics 14 01353 g005
Table 1. Detailed information about the optimal hyperparameters.
Table 1. Detailed information about the optimal hyperparameters.
Hyperparameters Candidate Value Determined Value
“Outer bags”1 to 11; step = 110
“Learning rate”[0.001, 0.005, 0.01]0.01
“Early stopping rounds”35 to 41; step = 137
“Max rounds”9000 to 10,000; step = 10010,000
“Max leaves”5 to 11; step = 110
Table 2. The descriptive and inferential statistics of AMI and Control groups in terms of some factors.
Table 2. The descriptive and inferential statistics of AMI and Control groups in terms of some factors.
VariablesGroupsp
AMI Control
(n = 65)(n = 34)
Gender<0.001
Female4 (16%)21 (84%)
Male61 (82.4%)13 (17.6%)
Age55.72 ± 9.0156.35 ± 8.860.74
BMI25.26 ± 3.7425.07 ± 3.860.82
Smoking
Yes22 (40.7%)32 (59.3%)<0.001
No43 (95.6%)2 (4.4%)
Table 3. The outputs show the binary classification performance of the EBM model with and without the addition of interaction terms.
Table 3. The outputs show the binary classification performance of the EBM model with and without the addition of interaction terms.
MetricInteraction Terms Added?Data SourceValueBCI * (95%)
AccuracyYesTrain1.00(0.99–1.00)
Test0.92(0.80–1.00)
NoTrain1.00(0.99–1.00)
Test0.84(0.68–0.96)
SensitivityYesTrain1.00(0.99–1.00)
Test0.89(0.67–1.00)
NoTrain1.00(0.99–1.00)
Test0.83(0.65–1.00)
SpecificityYesTrain1.00(0.99–1.00)
Test0.94(0.80–1.00)
NoTrain1.00(0.99–1.00)
Test0.86(0.50–1.00)
F1 scoreYesTrain1.00(0.99–1.00)
Test0.94(0.83–1.00)
NoTrain1.00(0.99–1.00)
Test0.88(0.73–0.98)
AUCYesTrain1.00(0.99–1.00)
Test0.95(0.83–1.00)
NoTrain1.00(0.99–1.00)
Test0.93(0.79–1.00)
*: Bootstrapped confidence interval with 1000 repetitions.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Arslan, A.K.; Yagin, F.H.; Algarni, A.; AL-Hashem, F.; Ardigò, L.P. Combining the Strengths of the Explainable Boosting Machine and Metabolomics Approaches for Biomarker Discovery in Acute Myocardial Infarction. Diagnostics 2024, 14, 1353. https://doi.org/10.3390/diagnostics14131353

AMA Style

Arslan AK, Yagin FH, Algarni A, AL-Hashem F, Ardigò LP. Combining the Strengths of the Explainable Boosting Machine and Metabolomics Approaches for Biomarker Discovery in Acute Myocardial Infarction. Diagnostics. 2024; 14(13):1353. https://doi.org/10.3390/diagnostics14131353

Chicago/Turabian Style

Arslan, Ahmet Kadir, Fatma Hilal Yagin, Abdulmohsen Algarni, Fahaid AL-Hashem, and Luca Paolo Ardigò. 2024. "Combining the Strengths of the Explainable Boosting Machine and Metabolomics Approaches for Biomarker Discovery in Acute Myocardial Infarction" Diagnostics 14, no. 13: 1353. https://doi.org/10.3390/diagnostics14131353

APA Style

Arslan, A. K., Yagin, F. H., Algarni, A., AL-Hashem, F., & Ardigò, L. P. (2024). Combining the Strengths of the Explainable Boosting Machine and Metabolomics Approaches for Biomarker Discovery in Acute Myocardial Infarction. Diagnostics, 14(13), 1353. https://doi.org/10.3390/diagnostics14131353

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop