Blog

KARACHI – A recent breakdown of commercial bank accounts’ balances in Pakistan reveals a significant concentration of low-value deposits, highlighting financial disparities across the banking population.

According to data compiled by the State Bank of Pakistan (SBP) and analyzed by Topline Securities, a staggering 51 percent of bank accounts in the country have less than Rs5,000 in deposits. The report sheds light on the limited financial capacity of a majority of account holders, indicating widespread financial under-inclusion.

The data further revealed that 20pc of accounts maintain balances between Rs5,000 and Rs50,000.

It adds that 26% hold deposits ranging from Rs50,000 to Rs1 million while only **3%** of accounts contain balances above Rs1 million.

This concentration of low-balance accounts suggests that while banking access may be growing, a large portion of the population is either unable or unwilling to maintain substantial deposits.

The findings also emphasize the role of microfinance, digital wallets, and low-cost banking solutions to uplift the underbanked segments of the population.

On the other hand, SBP’s foreign exchange reserves closed at US$ 14.51 billion as on 30 June 2025.

During FY25, SBP’s Fx reserves has recorded an increase of US$ 5.12 billion to reach US$ 14.51 billion as on 30 June 2025 compared to US$ 9.39 billion as on 30 June 2024. This reflects a noticeable improvement in the country’s current account balance and realization of planned inflows during the year.

The federal government is considering lifting the ban on new gas connections as a potential solution to the ongoing surplus in liquefied natural gas (LNG) supply, which is causing strain on Pakistan’s gas infrastructure and its international commitments. 

According to media reports, the Ministry of Finance is evaluating this demand-side intervention, which would allow new connections for industrial and residential consumers, as more than 3.5 million applications for new gas connections are pending with Sui gas companies. 

The moratorium, imposed in 2009, was partially lifted but reimposed in 2022 due to rising gas shortages.

The government is also facing a dilemma, as high electricity tariffs and the economic slowdown have forced industries to shut down their captive gas plants, while LNG infrastructure, valued at over $7 billion, remains underused. 

The original agreements for LNG supply mandated that power plants use a significant portion of imported LNG, but due to underutilisation, circular debt has shifted from the power sector to the petroleum sector, affecting upstream and downstream stakeholders.

To manage the LNG surplus, Pakistan has curtailed domestic gas production by 300 million cubic feet per day (mmcfd), a move that has resulted in significant losses for local producers and put pressure on future gas exploration. Additionally, several LNG cargoes from Qatar are surplus, and five cargoes from last winter still need to be rescheduled.

The government is proposing a new strategy to allocate domestic gas exclusively to residential consumers and use imported LNG for industrial, power, and fertilizer sectors. 

This would be combined with a revision of the weighted average cost of gas (WACOG) to introduce a uniform pricing model for local and imported gas.

As the government faces increasing pressure to address the gas supply issues, a 50% hike in fixed gas charges and up to a 17% increase in per-unit rates for industrial, power, and bulk users has been announced. This adjustment is aimed at addressing the revenue shortfalls of gas utilities, with a projected financial burden of Rs85 billion for FY2025-26. 

Additionally, the imposition of a captive levy on industrial gas users has led to a significant drop in demand, further distorting the market.

Introduction

Coronary atherosclerotic heart disease refers to the accumulation of fatty plaques in the epicardial arteries, coronary atherosclerosis leads to stenosis or occlusion of the lumen and leads to myocardial ischemia, hypoxia or necrosis, abbreviated as coronary heart disease(CHD),which is the most common type of organ lesions caused by atherosclerosis. It is a disease with high morbidity and mortality worldwide, with a high incidence of recurrent events in patients and insufficient prediction tools available.^1,2 Therefore, enhancing the prediction and prevention of MACE in CHD is essential to improve clinical prognosis and patient quality of life.

Studies have shown^3–5 that myocardial injury markers: creatine kinase isoenzyme (CK-MB), creatine kinase (CK), lactate dehydrogenase (LDH), aspartate aminotransferase (AST), troponin I (cTnI), troponin T (cTnT), ultrasensitive troponin (hs-cTnI), myoglobin (MYO), cardiac-type fatty acid binding protein (H -FABP), brain natriuretic peptide (BNP), and N-terminal B-type natriuretic peptide precursor (NT-proBNP) can be used as the core indexes for the diagnosis and prognosis of cardiovascular diseases, which need to be further studied in depth. Currently, some studies have shown that the prediction accuracy of hs-cTnI in predicting short-term MACE in patients with non-ST-segment elevation myocardial infarction reached 92.3% by logistic regression analysis;⁶ the use of randomised survival forest model for predicting cardiovascular outcomes has revealed that cTnT is an important predictor of the occurrence of heart failure (HF), and the model has a better prediction accuracy.⁷

The 4 core diagnostic methods of traditional Chinese medicine(TCM) – inspection, auscultation and olfaction, inquiry and palpation – are collectively known as the “Four Diagnostic Methods”, which form the basis of syndrome differentiation and treatment in TCM. Each method collects patient information from different perspectives and ultimately makes a comprehensive judgment on the nature of the disease. Among the four diagnostic methods in TCM, inspection is the most important. It is an important means for doctors to understand a patient’s physical condition and disease changes by observing their complexion, tongue body, tongue coating, eyes and other external manifestations. Inspection not only enables a preliminary diagnosis of the disease but also provides an important basis for the other three diagnostic methods. At present, the research on inspection mainly focuses on tongue diagnosis and facial diagnosis. Auscultation and olfaction is a diagnostic method that acquires patient information by smelling odors and listening to sounds. It is of great significance for judging internal organ diseases, the severity of the condition, and the prognosis of the disease. The position of inquiry in the four diagnostic methods has always been highly valued by all medical practitioners. The process of inquiry not only enables a judgment on the patient’s condition, reducing the interference of subjective consciousness in the other three diagnostic methods, but also allows for an understanding of the patient’s mental state, emotions, etc. It serves as a “bridge” between doctors and patients. However, unlike the other three diagnostic methods, inquiry does not have specific images, sounds, tastes, and pulses that can be referred to and transformed. The content of the medical consultation is not easy to quantify, and it also varies according to each doctor’s academic thoughts and clinical experience. At the same time, it is an important basis for the final determination of the disease. Palpation is a process where doctors use their fingers or palms to perform various forms of palpation on certain parts of the patient’s body to assess the condition and understand the illness. At present, palpation in TCM mainly refers to pulse diagnosis, which has been highly valued by medical practitioners throughout history. Different schools of thought have different interpretations of pulse diagnosis, which depends on the clinical experience and theoretical knowledge of medical practitioners, making it difficult to standardize pulse diagnosis in TCM.⁸

Tongue diagnosis is an important method of diagnosis in TCM, reflecting the state of systemic diseases through tongue images. However, traditional tongue diagnosis relies on the subjective experience of doctors, and the diagnostic results lack consistency and repeatability.⁹ With the application of artificial intelligence and image processing technology in the study of objectivation of the four diagnostic methods of TCM, the research of modern diagnostic techniques based on objectivated tongue diagnostic parameters has been initially developed in recent years. For example, the non-invasive diagnosis and identification model of CHD established by feature fusion algorithm, both of which have an accuracy rate of more than 80%, provide a new idea for large-scale screening and convenient diagnosis of CHD, and at the same time indicate that the tongue information may become an effective marker for the diagnosis of CHD;^10,11 in the preliminary stage of this project, we analysed the relationship between the tongue colour of the patients with coronary arteriography and the severity of the stenosis in coronary arteries through the retrospective cohort study, and the results showed that tongue colour parameters could provide a reference for predicting the degree of coronary artery stenosis, and it was also suggested that future studies could extend the tongue features, optimise the prediction model, and explore the application in other cardiovascular diseases.¹²

In recent years, the rapid development of artificial intelligence has provided new ideas for exploring non-invasive prognostic methods for CHD. Machine learning (ML) has advantages over traditional statistical methods when clinical data exhibit complexity and multidimensionality, eg, ML involves the selection and integration of multiple models, and also automatically learns to deal with non-linear relationships and select useful predictive features. Therefore ML can reveal hidden relationships in data more effectively and is increasingly used for diagnosis and risk prediction of clinical diseases.⁷ A study has developed an ML model that predicts mortality in patients with HF caused by CHD and evaluated the performance of the ML model to provide an explicit interpretation of individualised risk prediction through the combination of ML and SHAP, and to intuitively understand the impact of the key features in the model, leading to a better treatment plan and optimal allocation of resources for patients. In addition, the interpretable framework increases the transparency of the model and facilitates physicians to understand the reliability of the predictive model.^13,14

The innovativeness of this study lies in exploring the joint analysis of multimodal data (myocardial injury markers and tongue images) with the help of modern technologies such as ML, providing important theoretical support and methodological basis for the construction of the MACE prediction model for patients with CHD with the combination of TCM and Western medicines, promoting the individualised secondary prevention, and providing a certain research basis for the future development of a more efficient and accurate prognostic prediction tool for MACE in patients with CHD. The flow chart of this study is shown in Figure 1.

Figure 1 Flowchart.

Material

Research Subjects

All patients with CHD originated from the cardiology outpatient clinics and wards of Shuguang Hospital affiliated to Shanghai University of Traditional Chinese Medicine, Yueyang Combined Hospital of Traditional Chinese Medicine and Western Medicine, Jiading Central Hospital, AnTu Division of Shanghai YangPu Central Hospital and Shanghai YangPu Hospital of Traditional Chinese Medicine in the year 2023–2024. A total of 1,398 patients were collected, follow-up observation was performed at 6 months, 9 months and 12 months after the first collection, and clinical information and MACE were recorded, and the final total valid sample was 1,293 patients according to the inclusion criteria and exclusion criteria nadir, among which the number of cases without MACE was 1014; the number of MACE cases was 279. Retrospective analyses were performed based on the participants’ original clinical data and follow-up data.

Diagnostic Criteria for CHD

Coronary angiography (CAG) examination reveals that the diameter of the subepicardial coronary artery is narrowed by more than 50%, and the patient has typical angina pectoris symptoms, evidence of ischemic changes is shown through electrocardiogram (such as ST-segment elevation or reduction, T-wave inversion, etc)., which can be diagnosed as CHD.¹⁵

Inclusion Criteria

① Meet the western medical diagnostic criteria of CHD; ② Men and women are not limited, age 35–85 years old; ③ Patients’ informed consent; ④ Clinical literature and information are complete and accurate.

Exclusion Criteria

① Acute and chronic nephritis, urinary tract infections, hypertrophic cardiomyopathy, trauma, surgery, acute infectious diseases and severe hepatic and renal insufficiency; ② pregnancy, lactating women, mental disorders; ③Subjects who voluntarily propose to withdraw; ④ Serious complications or deterioration of the condition during the study, requiring emergency measures; ⑤Incomplete clinical information or failure to comply with the follow-up and data collection process.

Clinical Information

General Clinical Data

Basic information such as participants’ age, gender, BMI, blood pressure, heart rate, disease duration, smoking history, alcohol consumption history, and past medical history; laboratory test results (LVEF,Gensini and myocardial injury marker indicators) and follow-up results (whether recurrence of MACE), and tongue parameters were recorded.

MACEs in Patients with CHD

MACE is a composite indicator event of clinical events,^16,17 and is usually measured in clinical trials for cardiovascular patients. It may include multiple endpoint indicators, such as ① recurrent unstable angina pectoris: Patients report that the sudden onset lasts for a short period of time, with pain centered on the posterior part of the sternum, or spreading to the upper arm, left chest, and neck; ② Revascularization of the target vessels again: It refers to the revascularization surgery (PCI or coronary artery bypass grafting) performed on the infarct-related vessels treated by the previous PCI, including target vessel revascularization (TVR) and target lesion revascularization (TLR). The former refers to the revascularization of any part of the target vessel, including Target Lesion (TLR) and Non-Target Lesion (non-TLR); The latter specifically refers to the repeated intervention of the stent implantation site or within a 5mm range of its proximal/distal end; ③ Non-fatal recurrent myocardial infarction: The results of echocardiography, myocardial enzymes, etc. indicate myocardial ischemia and myocardial injury, myocardial cell death, myocardial necrosis, and local circulatory disorders and other pathological conditions; ④ Heart failure: According to the guidelines of the European Society of Cardiology, the diagnosis of heart failure refers to ventricular filling and/or ejection dysfunction caused by various structural or functional diseases of the heart, where the cardiac output cannot meet the metabolic needs of body tissues, resulting in clinical manifestations such as breathing difficulties, limited physical activity, and fluid retention. ⑤ Cardiac death. The fundamental and direct cause of death lies in the structural or functional abnormalities of the heart itself, mainly caused by sudden cardiac death, acute congestive heart failure, acute myocardial infarction, severe arrhythmia and other structural/functional heart diseases. ⑥ Stroke. The definition of stroke is based on imaging findings or typical symptoms.

Tongue Image Collection and Parameter Extraction

① Tongue image acquisition: Canon PowerShot SX720 HS digital camera and Aisle X-rite ColorChecker Classic mini professional 24-colour card were used to acquire tongue images of the CHD patients and extract the tongue parameters; the subjects avoided drinking water for half an hour, avoided eating for one hour, rested quietly for more than 5 minutes, took the orthopedic sitting position or the supine position, opened their mouths and stretched their tongues, relax the facial muscles, hold the position for 3 seconds, and take the tongue picture; the collection of tongue diagnostic parameters was carried out by 2 experienced TCM practitioners (M-Z, JY-L).

② Tongue parameters extraction: tongue body and tongue coat were separated and segmented using “SMX System 2.0 tongue image analysis software” (registration number: 2008SR12316). The tongue image parameters used in this study mainly take the RGB model and HSV model as measurement indicators to describe the color parameters of tongue tissue (TC_) and tongue coating (CC_), including the tongue colour parameters RGB, HSV, and the indices of tongue fatness and thinness, pitting, dentition, petechiae, etc.; and the coat colour parameters RGB, HSV, and the indices of coat thickthin, curdygeasy and peeledcoating.

As shown in Figure 2¹² the RGB model is an additive colour model that represents various colours through different combinations of intensities of the primary colours red, green and blue. The intensity of each colour is usually represented by an integer between 0 and 255, where 0 means no colour and 255 means maximum intensity. On the other hand, the HSV model is a colour model based on human visual perception. In this particular model, the H-value (hue) denotes the basic attributes such as red, green, and blue; the S-value (saturation) denotes purity or brightness; and the V–value denotes lightness.

Figure 2 RGB and HSV colour models.

In addition, since the tongue is a complex organ whose surface is covered with papillae and taste buds, different regions of the tongue may exhibit specific tongue features under the influence of different diseases or pathological processes. As shown in Figure 3, in order to observe and analyse these features in a more systematic way, the system parameters divide the tongue into root(ROOT), middle(MID), tip(TIP), left (LEFT)and right(RIGHT).

Figure 3 Methods for collecting and extracting parameters from tongue images.

Data Collection and Quality Control

In this study, all data were obtained from the original clinical literature and clinical data collection was carried out by postgraduate students. It was checked by the researchers and imported into the EpiData database. To ensure the accuracy of the data, each case was checked for consistency and revised by at least two data managers for any inconsistencies until the information was accurate. Two experienced postgraduate students were designated as quality control supervisors to oversee the data collection process, data checking and data quality control. All patient-related information in this study was kept intact for patient privacy considerations.

Statistical Analysis

Statistical Methods

All statistical tests were conducted using two-sided hypothesis tests at a level of P= 0.05 (P≤ value 0.05 indicates statistical significance). The specific principles were as follows: (1) Quantitative data were described by means and standard deviations, and hypothesis tests were performed using the t-test (normal distribution) or Wilcoxon rank sum test. (2) Qualitative data were described by frequencies and percentages, and hypothesis testing was done using chi-square test or Fisher’s test. Rank-level data were analysed using Wilcoxon rank sum test. After statistical analysis, radar charts were used to visualise the changes in the characteristic parameters of tongue body and tongue coat colour in different MACE groups.

Analysing Risk Factors for MACE Using Non-Parametric Tests

Non-parametric test was performed with MACE as the primary endpoint and myocardial injury markers and tongue parameters as covariates. Non-parametric test analysis was used to screen out significant variables (P= 0.05), which were estimated using maximum likelihood algorithm, and odd ratios (OR) and 95% CIs were calculated. Software SPSS 27 was used for data analysis. Feature selection is crucial in the process of building clinical prediction models, as some features may have no effect on the prediction results or even interfere with the validity of the prediction model. Although ML-based feature selection methods have the ability to automatically select features, their results are convoluted and difficult to interpret.By comparing the differences in myocardial injury markers and tongue parameters between different groups, we can directly understand which features are associated with MACE.

Feature Screening Is Carried Out After LASSO Dimension Reduction

Lasso regression compresses the variable coefficients in the regression model by generating a penalty function to prevent overfitting and solve the problem of severe collinearity. Lasso regression was first proposed by the British Robert Tibshirani and is currently widely used in prediction models. In the New Gran literature, it is recommended that for model fitting with too many variables and a small sample size, Lasso regression should be considered first for variable screening. Default cross-validation is 10% off.

ML Methods to Build Prediction Models and Verification

ML algorithms can extract useful information from massive data and make objective and accurate predictions through algorithmic models. We hope to analyse myocardial injury markers and tongue parameters through ML to more comprehensively assess the application value of myocardial injury markers and tongue parameters in predicting MACE. In this study, we explored and validated the hypothesis that the analysis of myocardial injury markers and tongue colour parameters can be used as a clinical auxiliary diagnostic method by building a MACE prediction model based on myocardial injury markers and tongue colour parameters. To compare the predictive ability of 5 different models (XGBoost, logistic regression, KNN, LightGBM, AdaBoost), a cross-validation approach was used with a number of times of 2, a random seed of 44, automated grid-seeking of parameters, calculation of accuracy, sensitivity, and specificity, and plotting of decision curve analysis (DCA). These are the classical ML models used to build predictive models that have been shown to be effective in previous heart disease studies.^18–20 In each modelling process, the dataset was randomly divided into a training set and a test set in a ratio of 8:2.The training set is used for model training and the test set is used to validate the model performance. All feature values are normalised by a standard scaler, which works by subtracting the mean and dividing by the variance so that all values are centred on zero with a variance of 1. During model training, parameters are optimised using hierarchical 2-fold cross-validation and grid search. The categorical metrics used to assess the performance of model training and validation are accuracy, sensitivity, specificity, and F1 score: ①Accuracy is the most commonly used categorical metric. It indicates the proportion of samples correctly predicted by the model to the total number of samples; ② Sensitivity indicates the proportion of all samples predicted as positive by the model that are actually positive. This helps us to understand how accurately the model predicts the positive class; ③ Specificity is the proportion of the number of samples correctly predicted as negative by the model to the number of actual negative samples; ④Recall indicates the proportion of samples correctly predicted by the model to be positive out of all samples that were actually positive;⑤The F1 score is the reconciled mean of precision and recall, which takes into account the performance of precision and recall; the higher the F1 score, the better the model performs in terms of precision and recall.The ROC curve is a graphical tool used to evaluate the performance of a classification model, which shows the relationship between the True Positive Rate (TPR) and the False Positive Rate (FPR) at different classification thresholds. The AUC value (area under the curve) represents the area under the ROC curve and reflects the overall performance of the model at different classification thresholds. The AUC values range from 0.5 to 1, where 0.5 indicates that the model performance is comparable to random guessing and 1 indicates perfect model performance. By comparing the AUC values of the five models, it is possible to see the performance ability of myocardial injury markers and tongue feature parameters in predicting MACE.

In order to improve the credibility of the model and facilitate its use by clinicians, it is important not only to report the prediction results but also to interpret the model. Traditional feature importance can only indicate the importance of the feature values, but cannot clearly describe their impact on the prediction results. Therefore, in our study, we used the SHapley Additive Explanatory Values (SHAP) for visualisation and analysis.The SHAP values not only reveal the importance of the features, but also provide a quantitative metric for assessing the contribution of each feature to the predicted outcome. This allows us to accurately compare the impact of different markers of myocardial injury, and different tongue features, on model output results. In addition, the SHAP values reveal the interactions between features, allowing us to gain a deeper understanding of the interactions between these features and a comprehensive understanding of the relationship between myocardial injury markers and MACE. In this study, the best ML algorithm for MACE will be used to plot the SHAP values for each sample using the Extreme Intelligence Analytics platform.

Results

Baseline Information

A total of 1398 patients with CHD from 1st January 2023 to 31st December 2024, finally 1293 participants were included in the final study according to the inclusion criteria and exclusion criteria, with a mean age of 70 years (the youngest was 34 years and the oldest was 85 years), 688 males (53.2%) and 605 females (46.8%);913 had hypertension (70.61%), 427 had diabetes mellitus (30.30%), 115 had hyperlipidaemia (8.90%), 475 had a history of smoking (36.73%),416 had a history of alcohol consumption (32.17%). According to the follow-up results, of the 1293 patients included in the study, MACE occurred in 279 (21.6%) participants.Since the sample sizes of the MACE group and the non-MACE group differ significantly, the SMOTE method was used to balance the samples of 279 MACE patients according to the characteristics of specific MACE events. Eventually, there were 1116 cases in the MACE group and 1014 cases in the non-MACE group.There was no statistically significant difference between the MACE group and the non-MACE group in terms of age, gender, systolic blood pressure, diastolic blood pressure, heart rate, BMI, history of smoking, history of alcohol consumption, duration of the disease, and previous history (history of hypertension, diabetes mellitus, and history of hyperlipidaemia) (P>0.05), and in terms of whether or not the group had PCI (P< 0.05) was statistically significant (Table 1).

Table 1 Comparison of Baseline Data Between the Non-MACE Group and the MACE Group [n(%) or M(Q1, Q3)]

Non-Parametric Test of Risk Factors for MACE in Patients with CHD

Non-parametric test was used to screen the risk factors for the occurrence of MACE in patients with coronary heart disease. Myocardial injury markers, tongue parameters and other indicators were included in the data. The results showed that:Gensini, LVEF,LDH,MYO,CKMB,cTnT,BNP,NTproBNP,TC_ROOT_R,TC_ROOT_G,TC_ROOT_B,TC_ROOT_S,TC_ROOT_V,TC_MID_S,TC_TIP_B,TC_TIP_S,TC_RIGHT_R,TC_RIGHT_V,TC_R,TC_S,TC_V,CC_ROOT_H,CC_ROOT_S,CC_MID_G,CC_MID_B,CC_MID_H,CC_MID_S,CC_TIP_R,CC_TIP_G,CC_TIP_B,CC_TIP_H,CC_TIP_V,CC_LEFT_R,CC_LEFT_G,CC_LEFT_B,CC_LEFT_V,CC_RIGHT_G,CC_RIGHT_S,CC_G,CC_B,CC_H,CC_S,CURDYGREASY_VAL,THICKTHIN_VAL,PEELEDCOATING_VAL (p <0.05) by Mann–Whitney-U test showed a statistically significant difference between the groups, indicating that these indicators led to the occurrence of MACE in this study (p=0.05) (Table 2, note: Table 2 only includes statistically significant results; see Table 1S for full table). Radar charts visualised the changes in tongue colour,coat colour and tongue shape characteristic parameters in different MACE groups (Figures 4–5).

Table 2 Non-Parametric Tests for the Presence or Absence of Risk Factors for the Development of MACE in Patients with CHD

Figure 4 Radar map of tongue colour distribution.

Figure 5 Radar map of coat colour distribution.

Multivariate Logistic Regression Analysis of Risk Factors for MACE in Patients with CHD

Multivariate logistic regression analysis was conducted on imaging indicators and myocardial injury markers in patients with CHD based on whether MACE occurred. Multivariate regression analysis was carried out when the single-factor P value was <0.1, and two-way stepwise regression analysis was also conducted. The results showed that the differences in Gensini, AST, LDH and NT-proBNP in the multivariate regression analysis were statistically significant (P<0.05). Due to the large collinearity of tongue image parameters, it was not involved in the multivariate regression analysis. See Table 3 for details.

Table 3 Multivariate Logistic Regression Analysis of Risk Factors for MACE in Patients with CHD

Feature Dimension Reduction

Feature dimension reduction was performed using LASSO, and the obtained feature factors were: Gensini score, AST, LDH, NT-proBNP, TC_ROOT_R, TC_MID_S, TC_TIP_B, CC_MID_B, CC_TIP_G, and PEELEDCOATING_VAL. Figure 6 is the cross-validation curve of LASSO regression, and Figure 7 is the cross-sectional view of LASSO coefficients. The LASSO coefficient Table can be found in Table 2S.

Figure 6 LASSO regression cross-validation curve.

Figure 7 LASSO coefficient profile.

ML Modelling to Compare the Joint Predictive Ability of Different Indicators

The XGBoost, logistic, KNN, LightGBM and AdaBoost models were used to predict MACE in order to train and verify the myocardial injury markers and tongue color parameters of the two groups of patients. The inclusion factor is the feature screening factor after LASSO dimension reduction. As shown in Table 4, except for the logistic model, the other 4 algorithms all exhibit relatively high verification results. In terms of the training set, XGBoost, KNN, and LightGBM have superior prediction performance, with accuracy, precision, sensitivity, specificity, and F1 scores all exceeding 70%. As shown in Table 5, in terms of the test set, XGBoost and LightGBM have the highest prediction performance. We integrated the training process of each model and 2x cross-validation using the training dataset. Based on the evaluation results of the two datasets (training and validation), XGBoost has been proven to be a highly stable model based on tongue image parameters and also the best predictive model for determining MACE. Figure 8 shows the ROC curves of different models using myocardial injury markers and tongue image parameter data. In this ROC curve graph, lines of different colors represent whether MACE occurs. The results of the validation set show that the predictive abilities of XGBoost and LightGBM for the occurrence of MACE are greater than 0.9 (AUC > 0.97). The ROC curve of the validation set is shown in Figure 4S, the calibration curve is shown in Figure 1S, the decision curve validation is shown in Figure 2S, the forest plot is shown in Figure 3S, and the PR curve of the training set is shown in Figure 4S.

Table 4 Whether MACE Has Occurred Training Set

Table 5 Whether MACE Has Occurred Validation Set

Figure 8 ROC curves for different models using myocardial injury markers, tongue parameter data.

SHAP Feature Values

According to the XGBoost model visualisation report, the magnitude of the output contribution of the 10 parametric features is shown in Figure 9, they are Gensini, AST, LDH, NT-proBNP, TC_ROOT_R, TC_MID_S, TC_TIP_B, CC_MID_B, CC_TIP_G, PEELEDCOATING_VAL,These features are associated with MACE.

Figure 9 Summary of the characteristics of MACE-related impact parameters.

In order to show more visually the output contribution of different groups of myocardial injury markers and tongue features, the SHAP value plot further shows the effect of myocardial injury markers and tongue colour parameter values on MACE. Each row represents a parameter and each dot represents a theme. The colour represents the size of the feature parameter, the larger the value, the redder the colour. The results are shown in Figure 10: High values of Gensini, CC_TIP_G, and low values of NT-proBNP, LDH, TC_ROOT_R, AST and TC_TIP_B positively predict MACE.

Figure 10 SHAP value plot of XGBoost for MACE features.

Conduct Subgroup Analysis on Specific MACE Events

A total of 1,116 MACE events were included in this study. According to the above classification and statistics, there were 636 cases of unstable angina pectoris, 308 cases of emergency revascularization, 52 cases of myocardial infarction, 60 cases of heart failure, 8 cases of cardiogenic death, and 52 cases of stroke.

In previous studies, MACE could be directly subdivided into different types of events, such as hospitalization and death.²¹ Some studies also grouped MACE according to severity (no, mild, moderate, severe) to analyze the association between different severity and clinical outcomes.²² There are also studies that classify MACE into in-hospital MACE and non-in-hospital MACE based on whether MACE occurs in the hospital,²³ and some studies adopt stratified reporting (major MACE vs Extended MACE) has been widely accepted and applied in cardiovascular clinical research and is supported by a large number of authoritative literatures and guidelines. The core purpose of this stratified approach is to provide more comprehensive clinical event information while ensuring the comparability of core endpoints.^24–26 For example, in its numerous ACS and secondary prevention studies (such as PROVE IT-TIMI 22, TRITON-TIMI 38, PEGASUS-TIMI 54, etc)., the standard primary endpoint is usually a triple composite endpoint of cardiovascular death, myocardial infarction, and stroke. Secondary endpoints or exploratory endpoints usually include emergency revascularization, hospitalization due to unstable angina pectoris, hospitalization due to heart failure, etc. These events are often combined into extended composite endpoints.^27,28 Therefore, in this study, only 112 cases of the primary endpoint events were discussed (8 cases of cardiovascular death, 52 cases of myocardial infarction, and 52 cases of stroke).

Baseline Analysis of Patients with Different Primary Endpoint Events of MACE

As shown in Table 6: There were statistically significant differences among the groups of patients in terms of heart rate, BMI, disease duration, smoking history, history of hypertension, history of diabetes, and history of hyperlipidemia (P<0.05).

Table 6 Baseline Analysis of Patients with Different Primary Endpoint Events of MACE

Comparison of Myocardial Injury Markers in Patients with Different Primary Endpoint Events of MACE

As shown in Table 7, AST,LDH,MYO and BNP of the three groups of patients were (P<0.05) statistically significant by Kruskal–Wallis test. The AST of patients with myocardial infarction was significantly higher than that of patients with stroke and death, the MYO was significantly lower than that of patients with stroke and death (P<0.01), the LDH of patients with stroke was significantly higher than that of patients with myocardial infarction and death (P<0.05), and the BNP of patients with death was significantly higher than that of patients with stroke and myocardial infarction (P<0.01).

Table 7 Analysis of Myocardial Injury Markers in Patients with Different Primary Endpoint Events of MACE

Construction of Prediction Models for Different Primary Endpoint Events of MACE and Comparison of the Combined Predictive Capabilities of Different Indicators

Stroke Patients and Non-Stroke Patients

In this study, a total of 1,116 patients with MACE were mainly included, among whom 52 were stroke patients and 1,064 were non-stroke patients.When using LASSO for feature dimension reduction, the variable selection of the corresponding model is:Gensini+LVEF+LDH+MYO+CKMB+cTnT+BNP+NTproBNP+TC_ROOT_R+TC_ROOT_B+TC_ROOT_V+TC_MID_R+TC_MID_G+TC_MID_B+TC_MID_H+TC_MID_V+TC_TIP_R+TC_TIP_G+TC_TIP_B+TC_TIP_H+TC_TIP_S+TC_LEFT_R+TC_LEFT_B+TC_LEFT_H+TC_LEFT_V+TC_RIGHT_R+TC_RIGHT_B+TC_RIGHT_H+TC_RIGHT_V+TC_R+TC_G+TC_B+TC_H+TC_S+CC_ROOT_G+CC_ROOT_B+CC_ROOT_H+CC_ROOT_S+CC_MID_R+CC_MID_B+CC_MID_S+CC_MID_V+CC_TIP_B+CC_TIP_H+CC_TIP_S+CC_LEFT_R+CC_LEFT_G+CC_LEFT_B+CC_LEFT_S+CC_LEFT_V+CC_RIGHT_R+CC_RIGHT_B+CC_RIGHT_H+CC_RIGHT_S+CC_RIGHT_V+CC_R+CC_G+CC_B+CC_H+CC_S+CRACK_VAL+SPOTSTHORN_VAL+CURDYGREASY_VAL+MOISTDRY_VAL+THICKTHIN_VAL+PEELEDCOATING_VAL.

The XGBoost, logistic, KNN, LightGBM and AdaBoost models were used to predict the occurrence of stroke events. Five-fold cross-validation was employed to train and validate the myocardial injury markers and tongue color parameters of the two groups of patients. The inclusion factor is the feature screening factor after LASSO dimension reduction. Figure 11 shows the ROC curves of different models for verifying the occurrence of stroke.The results show that the three models of XGBoost, AdaBoost, and LightGBM in the training set and validation set have overfitting (AUC=1), which may be caused by the excessive difference in sample size. It is recommended to expand the sample size in the future. Tables 3S and 4S show the specific results of the training set and feature set respectively.

Figure 11 The ROC curves of stroke occurrence were verified by different models.

Patients with Myocardial Infarction and Patients Without Myocardial Infarction

In this study, a total of 1,116 patients with MACE were mainly included, among whom 52 were patients with myocardial infarction and 1,064 were non- myocardial infarction patients.

When using LASSO for feature dimension reduction, the variable selection of the corresponding model is:Gensini+LVEF+AST+LDH+MYO+CKMB+cTnT+BNP+NTproBNP+TC_ROOT_R+TC_ROOT_H+TC_ROOT_S+TC_ROOT_V+TC_MID_R+TC_MID_H+TC_MID_S+TC_TIP_R+TC_TIP_G+TC_TIP_H+TC_TIP_V+TC_LEFT_R+TC_LEFT_G+TC_LEFT_H+TC_LEFT_S+TC_LEFT_V+TC_RIGHT_R+TC_RIGHT_G+TC_RIGHT_H+TC_RIGHT_S+TC_R+TC_G+TC_B+TC_H+TC_V+CC_ROOT_R+CC_ROOT_G+CC_ROOT_B+CC_ROOT_H+CC_ROOT_V+CC_MID_R+CC_MID_G+CC_MID_S+CC_TIP_H+CC_TIP_S+CC_LEFT_R+CC_LEFT_G+CC_LEFT_V+CC_RIGHT_R+CC_RIGHT_H+CC_RIGHT_S+CC_RIGHT_V+CC_R+CC_V+CRACK_VAL+SPOTSTHORN_VAL+ENLARGEDTHIN_VAL+MOISTDRY_VAL+PEELEDCOATING_VAL.

The XGBoost, logistic, KNN, LightGBM and AdaBoost models were used to predict the occurrence of myocardial infarction. Five-fold cross-validation was employed to train and validate the myocardial injury markers and tongue color parameters of the two groups of patients. The inclusion factor is the feature screening factor after LASSO dimension reduction. Figure 12 shows the ROC curves of different models for verifying the occurrence of myocardial infarction. The results show that the three models of XGBoost, AdaBoost, and LightGBM in the training set and validation set have overfitting (AUC=1), which may be caused by the excessive difference in sample size. It is recommended to expand the sample size in the future. Tables 5S and 6S show the specific results of the training set and feature set.

Figure 12 The ROC curves of different models for verifying the occurrence of myocardial infarction.

Deceased Patients and Non-Deceased Patients

In this study, a total of 1,116 patients with MACE were mainly included, among whom 8 died and 1,108 were non-deceased.

When using LASSO for feature dimension reduction, the variable selection of the corresponding model is:Gensini+AST+BNP+TC_ROOT_B+TC_ROOT_H+TC_MID_S+TC_TIP_H+TC_TIP_S+TC_LEFT_H+TC_RIGHT_R+TC_RIGHT_V+CC_MID_R+CC_MID_V+CC_TIP_S+CC_LEFT_S+CC_RIGHT_S+CRACK_VAL+SPOTSTHORN_VAL+ENLARGEDTHIN_VAL+CURDYGREASY_VAL+PEELEDCOATING_VAL.

The XGBoost, logistic, KNN, LightGBM and AdaBoost models were used to predict the occurrence of myocardial infarction. Five-fold cross-validation was employed to train and validate the myocardial injury markers and tongue color parameters of the two groups of patients. The inclusion factor is the feature screening factor after LASSO dimension reduction. Figure 13 shows the ROC curves of different models for verifying the occurrence of death. Tables 7S and 8S show the specific results of the training set and the feature set.

Figure 13 The ROC curves of death were verified by different models.

Discussion

Relationship Between Tongue and MACE

The diagnostic and therapeutic method of inferring the heart’s vitality by observing the tongue image is called “Si-wai-chuai-nei”. The tongue image can objectively and accurately reflect the deficiency of the internal organs, the vitality of the qi and blood, the surplus and deficit of the fluid, the nature of the disease and the evil, the severity of the disease, its regression and the prognosis, etc., and it is a very important basis for the clinical diagnosis and treatment. “As the heart opens its orifices through the tongue”, tongue image plays an important role in the diagnosis of CHD, and modern medicine also believes that nutritional status, microcirculation status, and hormone levels can affect the formation of tongue image.

In this study, the following tongue factors were found to be related to MACE:TC_ROOT_R, TC_MID_S, TC_TIP_B, CC_MID_B, CC_TIP_G, PEELEDCOATING_VAL.These tongue factors are broadly focused on the tongue colour at the root of the tongue, the tip of the tongue coat colour and peeling coat. According to the division of the tongue and the relationship between the organs, the root of the tongue corresponds to the kidneys, the tip of the tongue corresponds to the heart, lungs, CHD in the heart, but with the lungs, liver, spleen, kidneys, and other organs are related to the heart, the heart is the master of the blood, the lungs are the master of the joints, the two are coordinated with each other, the operation of qi and blood from the smooth. Heart disease yang qi insufficiency, can not promote the blood flow on the tongue, lung qi governance section failure, then blood stasis, tongue loss of moistening, so the tip of the tongue to see the light greenish purple; kidney is the foundation of the innate, kidney yang for the root of the body’s yang qi, kidney yang insufficiency will lead to the heart yang is not vibrant, the heart veins lose warmth, which triggered thoracic paralysis cardiac pain, renal yin insufficiency is the deficiency of the heart yin, the heart fire is exuberant, burns the Jin for phlegm, phlegm stasis paralysis of the heart veins, resulting in thoracic paralysis cardiac pain. The findings are consistent with the pathogenesis of chest paralysis in TCM.²⁹ The Shanghanlunbenzhi – identify the tongue coat’said: ‘tongue coat by the stomach in the anger of the appearance of the stomach qi from the heart and spleen, so no disease often have thin coat, is the stomach in the anger of the grass, such as the ground of the micro-herb also’. If the evil invasion of the body or the body of positive Qi weakness leads to the imbalance of yin and yang, qi and blood, fluid, tongue coating also occurs a variety of changes.

Previous studies have suggested that a sudden retreat of the tongue coating or the appearance of floral peeling of the coating suggests that the gastric qi is about to be extinguished and the disease is in danger, and the analysis of the present study also concludes that the peeling of the coating is correlated with MACE. All of the above risk factors reflect that “tongue is a candidate outside the heart”, and thus tongue image can provide an important basis for predicting the judgement of MACE in CHD.

Currently, domestic studies on tongue and MACE still focus on small-sample, single-centre, correlation studies, and there are relatively few large-sample, multi-centre, prospective studies. The results of previous studies have shown^30–35 that old tongue, young tongue, tooth-marked tongue, green or blue-purple tongue, pale white tongue, peeling moss (excluding type of peeling moss), and purplish or reddish-purple sublingual veins are all relevant tongue factors for MACE in CHD; the scientific research team led by academician Chen KX has found that the indicators of purplish-purple sublingual veins, old tongue, and greenish-purple or greenish-purple tongue colour are closely correlated to MACE; a study has adopted the methodology of retrospective cohort study to establish and validate that stable CHD occurs in a single-centre, multi-centre, and prospective study. A study used a retrospective cohort study to establish and validate a combined Chinese and Western medicine prediction model for the occurrence of MACE in stable CHD, in which the dark red tongue and blood stasis were risk factors for the occurrence of cardiovascular events, with a validation set AUC value of 0.714; a study explored the correlation between the greenish purple tongue and the degree of risk of patients with acute myocardial infarction by observing the tongue of patients with acute myocardial infarction, and the results found that patients in the greenish purple tongue group had a higher degree of risk, more serious myocardial cell injury. The results found that the patients in the group with blue and purple tongue had a higher degree of danger, worse cardiomyocyte damage, worse cardiac function, and more serious inflammatory reaction, indicating that the appearance of blue and purple tongue in patients with acute myocardial infarction had a certain correlation with auxiliary examination indexes of the degree of danger of the disease, suggesting that tongue diagnosis of TCM, with blue and purple tongue has an important clinical value in the judgement of the degree of danger; it has been found that the level of BNP of patients with purple tongue is significantly higher than those with non-purple and purple tongue (P<0.05); the level of heart rate shock (HRT) is was significantly attenuated (P<0.05), suggesting that purple-darkened tongue may have some predictive value for the cardiovascular prognosis of patients with myocardial infarction; mirrored tongue is an independent risk factor for the occurrence of MACE in patients with acute heart failure (AHF) during the first year, and influences the occurrence of MACE in conjunction with age, length of hospitalisation, ejection fraction and NT-proBNP.

Relationship Between NT-proBNP and MACE

NT-proBNP is the most commonly used natriuretic peptide and plays a diagnostic role in the assessment of HF. It may be increased due to cardiac systolic or diastolic dysfunction, left ventricular hypertrophy, valvular heart disease, ischaemia or a combination of these factors. In addition to its use in the diagnosis of HF, NT-proBNP levels are very useful for risk stratification and management of patients with suspected HF: a trend towards decreasing levels of natriuretic peptide in the management of HF suggests an effective management strategy. In addition, several studies have demonstrated that NT-proBNP is also strongly associated with MACE in a variety of primary prevention and general populations;³⁶ in prospective cohort studies, it was found that patients undergoing multiple stenting were at increased risk of MACE with significantly elevated levels of NT-proBNP post-PCI, whilst other studies have reported a significant increase in the risk of MACE in patients who underwent multiple stents and had an increase in NT- proBNP levels increased more than threefold;³⁷ a review has shown that the predictive performance of the Revised Cardiac Risk Index (RCRI) in predicting MACE is improved with the addition of NT-proBNP, troponin, or their combination, and other studies have shown that NT-proBNP alone is more effective in predicting MACE. NT-proBNP may even have a higher discriminatory performance than RCRI when used alone;³⁸

In a related study, NT-proBNP and hsTnT levels were found to be associated with the risk of future MACE in patients with primary and secondary prevention of T2DM, and these two cardiac biomarkers can help to identify patients at very high risk of atherosclerotic events;³⁹ and in a related cohort study, NYHA class III–IV, NT-proBNP, and β-blocker medication were found, five variables, including the presence of LGE and overall longitudinal strain (GLS) of the left ventricle, were significantly associated with MACE and were used to construct column line plots, which obtained good C-index differentiation in the development and validation cohorts, respectively;⁴⁰ and a study found that traditional cardiovascular risk factors (NT-proBNP) were associated with a higher risk of MACE.⁴¹

The Relationship Between LDH and MACE

LDH is widely present in various tissues of the human body (such as the heart, liver, skeletal muscle, red blood cells, kidneys, lungs, brain, etc). When these tissues are damaged due to ischemia (such as myocardial infarction), hypoxia, inflammation, necrosis, malignant tumors or mechanical injury (such as hemolysis), LDH within the cells will be released into the blood, leading to an increase in serum LDH levels. Elevated LDH levels may be associated with dyslipidemia, elevated low-density lipoprotein cholesterol (LDL-C), and the formation of atherosclerotic plaques, further promoting the occurrence of MACE.⁴² The core events of MACE (such as myocardial infarction and heart failure) themselves can cause damage to myocardial cells and lead to an increase in LDH. Therefore, after acute myocardial infarction (AMI), LDH and its isoenzyme (LDH-1) were once regarded as important diagnostic and monitoring indicators (although troponin now has higher specificity and sensitivity). The fundamental causes of MACE often involve ischemia and hypoxia of the myocardium or brain tissue caused by atherosclerosis and thrombosis. Under hypoxic conditions, cell metabolism shifts from aerobic respiration to anaerobic glycolysis, generating a large amount of lactic acid. At the same time, the activity of LDH (which catalyzes the interconversion of pyruvate and lactic acid) increases. The continuous hypoxia eventually leads to cell damage and the release of LDH. The LDH level can be used as the basis for patient stratification. For example, in acute aortic dissection, 204 U/L is the optimal cutoff value of LDH, which is used to distinguish high-risk patients.⁴³ LDH combined with other biomarkers (such as AST, CK-MB, cTnI, α-HBDH) can significantly improve the accuracy of MACE prediction.⁴⁴

The Relationship Between AST and MACE

AST is abundant in cardiomyocytes. When cardiomyocytes are damaged due to ischemia (such as coronary artery occlusion), hypoxia, inflammation or necrosis, the cell membrane ruptures, and the AST inside the cells is released into the blood, resulting in an increase in serum AST levels. In diagnosed AMI patients, the peak level of AST or the area under the curve is often used to roughly estimate the extent of myocardial infarction. The larger the infarction area, the more extensive the myocardial necrosis, the more AST is released, the higher the peak value, and the poorer the prognosis is usually (the risk of MACE such as heart failure, cardiogenic shock, and death is higher). Among AMI patients, the degree of elevated AST is independently associated with short-term and long-term poor prognosis (death, heart failure, reinfarction). A high peak or continuous increase often indicates a large infarction area and severe impairment of cardiac function. Studies have found that an elevated AST/ALT ratio is significantly associated with the risk of MACCE (Major adverse cardiovascular and cerebrovascular events) in patients with unstable angina pectoris, with an optimal threshold of 1.29. This suggests that the AST/ALT ratio can be used as an indicator to assess cardiovascular risk. However, AST itself is not a direct indicator for independently predicting MACE. The judgment needs to be made in combination with the clinical background (such as coronary artery lesions, metabolic syndrome) and comprehensive assessment tools (such as CT-FFR, F4I).⁴⁵ It was also found in this study that the AST of patients with myocardial infarction was significantly higher than that of patients with stroke and death (P<0.01).

The Relationship Between Gensini and MACE

The Gensini Score is an angiographic scoring system used to quantify the severity of coronary atherosclerotic lesions. It has a clear and close positive correlation with the risk of MACE (Major Adverse Cardiovascular Events). In short, the higher the Gensini score, the more severe the coronary artery lesion, and the higher the risk of developing MACE in the future. A high score indicates a more extensive and severe coronary artery stenosis, leading to a larger range and more severe degree of myocardial insufficiency (ischemia), which is the direct pathological basis for events such as myocardial infarction, angina pectoris, and ischemic heart failure. Although Gensini mainly assesses the degree of stenosis, extensive severe lesions are often accompanied by the presence of more vulnerable (unstable) plaques, and the rupture of these plaques is the main cause of acute coronary syndrome (ACS, such as STEMI/NSTEMI). Diffuse severe stenosis causes the myocardium to remain on the edge of ischemia for a long time, resulting in poor cardiac functional reserve and making it more prone to decompensation (heart failure) under stress (such as infection, surgery, arrhythmia). Patients with high scores often require more complex interventional therapy (PCI), and their risk of postoperative restenosis or stent thrombosis is relatively high. Even without intervention, severely stenotic areas are more prone to thrombosis and occlusion. Previous studies have mentioned that the Gensini score is an independent risk factor for predicting MACE in ACS patients after PCI. This indicates that the Gensini score has a certain role in predicting MACE. Although the sensitivity and specificity of the Gensini score in predicting MACE when used alone are limited, as a standardized assessment tool for the degree of coronary artery lesions, It can still serve as an important reference for comprehensive risk stratification, especially having significant advantages when combined with other biomarkers or clinical indicators.⁴⁶ In this study, the SHAP value graph also revealed that a high value of Gensini could positively predict MACE.

Application of ML in Predicting MACE

ML has a wide range of applications in predicting MACE; one study used cardiac magnetic resonance (CMR) parameters and clinical features for ML-empowered risk stratification, which was able to surpass the classical Risk-SCD model of sudden cardiac death (SCD) in hypertrophic cardiomyopathy (HCM) and accurately predict MACE; in addition, the non-linear correlation between the CMR features (LGE and left ventricular pressure gradient) and MACE found in this study provides valuable insights into the clinical assessment and management of HCM;⁴⁷ a study evaluated the ability of different ML and traditional statistical methods in predicting MACE and found that ML methods improved the identification and prediction of patients at risk of MACE compared to traditional statistical methods even with a small sample size;⁴⁸ a study found that ML-based noninvasive fractional flow reserve (FFRCT) has prognostic value in predicting MACE after PCI in low and intermediate-risk patients;⁴⁹ and it was found that ML can accurately predict the risk of MACE in patients with suspected CHD who underwent SPECT MPI and CCTA, and that the ML feature ordering also shows which features are key to making such predictions, both at the sample and patient level.⁵⁰

Potential Causes of Insufficient Model Performance

Data limitations: myocardial injury markers (such as troponin and CK-MB) have high specificity in predicting coronary heart disease events, but insufficient sample size or poor data quality (such as improper marker detection time window) may lead to model failure; (2) Feature fusion problem: The physiological correlation between tongue image parameters and myocardial markers should be clear (such as the pathological relationship between TCM tongue diagnosis and myocardial ischemia). If the feature selection does not reflect the synergistic effect, it may weaken the effect of the model; (3) Method suitability: Traditional models (such as Logistic regression) have limited ability to capture nonlinear relationships, while complex models (such as random forests) are easy to overfit in small samples.

Strengths and Limitations

Although this study adopted strict quality control measures throughout the process of data collection, processing and analysis to ensure the accuracy and reliability of the results, and systematically investigated the relationship between markers of myocardial injury, tongue parameters and MACE, there are still some limitations: firstly, the sample size is relatively small, which may not fully reflect the real situation, and future studies will expand the sample size and external validation to better verify the persuasiveness of this prediction model; secondly, the number of groups with and without MACE varied greatly. Although the SMOTE processing method was adopted, it still had a certain bias effect on the results; in addition, as a retrospective study, it may bring information bias and selection bias, and prospective design in future studies can mitigate the impact of these biases; finally, the current study lacked in-depth exploration of the pathogenesis of MACE and further studies are needed in the future to investigate the molecular association between tongue microvascular morphology and coronary microcirculatory disorders (eg, ET-1, VEGE expression), and combined with the dynamic changes of myocardial markers (such as the early increase of myoglobin and the long window period of troponin), the insufficiency of time series feature extraction was discussed.

Summary

In this study, through a retrospective cohort analysis, the characteristic changes of myocardial injury markers and tongue image parameters in patients with or without MACE were explored. LASSO was used for feature dimension reduction to verify the effectiveness of myocardial injury markers and tongue image parameters in predicting MACE. Finally, subgroup analyses were conducted for different MACE events. These findings not only provide a modern scientific basis for TCM tongue diagnosis, but also provide new ideas and methods for early diagnosis and condition assessment of MACE. Future studies can further expand the scope of combined Chinese and Western medicine research, optimise the prediction model and explore its application value in other cardiovascular diseases.

Data Sharing Statement

The authors will provide the raw data supporting the conclusions of this paper without reservation.

Ethical Review Statement

This research has been reviewed and approved by the Ethics Committee of Shanghai University of Traditional Chinese Medicine, with the approval number: 2023-3-10-08-07. At the same time, it is confirmed that informed consent has been obtained from the research participants and that the guidelines outlined in the Declaration of Helsinki have been adhered to and confirmed.

Author Contributions

All authors made a significant contribution to the work reported, whether that is in the conception, study design, execution, acquisition of data, analysis and interpretation, or in all these areas; took part in drafting, revising or critically reviewing the article; gave final approval of the version to be published; have agreed on the journal to which the article has been submitted; and agree to be accountable for all aspects of the work.

Funding

The authors declare that the research, writing and/or publication of this article was supported by funding. This study was supported by grants from the National Natural Science Foundation of China (82374336 and 82074333), the Shanghai Key Laboratory of Health Identification and Assessment (21DZ2271000) and the Shanghai Pudong New Area Famous Chinese Medicine Studio Construction Project (PDZY-2025-0727, PDZY-2025-0951).

Disclosure

The authors declare that this study was conducted without any business or financial relationship that could be considered a potential conflict of interest.

References

1. ESC Scientific Document Group. ESC guidelines for the management of chronic coronary syndromes. Eur Heart J. 2024;45(36):3415–3537.

2. National Health and Family Planning Commission Expert Committee of Rational Drug Use, Chinese Pharmacists Association. Coronary heart disease and reasonable medication guide (2nd edition). Chin J Med Front. 2018;10(6):1–130.

3. Xintong L, Peng L. Research status of the predictive value of TCM syndrome differentiation in coronary heart disease and heart failure. Xinjiang J Tradl Chin Med. 2024;42(01):88–91.

4. Lihong Z, Yali L, Jiuyang W. Analysis of changes of myocardial injury markers in patients with myocardial infarction and their relationship with prognosis. J Med Forum. 2024;45(16):1686–1691.

5. Haodong W. Development of Metal Composite Nano-Enzyme-Labeled Side-Stream Immunochromatography Test Strips for Rapid Detection of Myocardial Injury Markers. Jilin Univ. 2024.

6. Zebin G, Yan L, Yahui L, et al. The predictive value of high-sensitivity troponin I for the short-term prognosis of patients suspected of acute coronary syndrome. Chin J Geriatric Cardiovascu Cerebrovascu Dis. 2024;26(07):751–754.

7. Ambale-Venkatesh B, Yang X, Wu CO. Cardiovascular event prediction by machine learning: the multi-ethnic study of atherosclerosis. Circ Res. 2017;121(9):1092–1101. doi:10.1161/CIRCRESAHA.117.311312

8. Zezheng K, Linhua Z, Chongxiang X, et al. Research progress on the modernization of the four diagnostic methods in traditional Chinese medicine. World Chin Med. 2025. 25,20;02:354–359.

9. Wang L, Nan J, Bin C, et al. Visual analysis of the development of multi-dimensional tongue diagnosis technology based on citespace. J Shaanxi Univ Chin Med. 2024;25:1–7.

10. Duan M. Research on machine learning-assisted diagnosis of coronary heart disease in traditional Chinese and western medicine based on tongue image characteristics. Beijing Univ Chin Med. 2023.

11. Duan M, Zhang Y, Liu Y. Machine learning aided non-invasive diagnosis of coronary heart disease based on tongue features fusion. Technol Health Care. 2024;32(1):441–457. doi:10.3233/THC-230590

12. Li J, Xiong D, Xu Z. Tongue color parameters in predicting the degree of coronary stenosis: a retrospective cohort study of 282 patients with coronary angiography. Front Cardiovasc Med. 2024;11:1436278. doi:10.3389/fcvm.2024.1436278

13. Wang K, Tian J, Zheng C, et al. Interpretable prediction of 3-year all-cause mortality in patients with heart failure caused by coronary heart disease based on machine learning and SHAP. Comput Biol Med. 2021;137:104813. doi:10.1016/j.compbiomed.2021.104813

14. Li J, Liu S, Hu Y, Zhu L, Mao Y, Liu J. Predicting mortality in intensive care unit patients with heart failure using an interpretable machine learning model: retrospective cohort study. J Med Internet Res. 2022;24(8):e38082. doi:10.2196/38082

15. Interventional Cardiology Group. Cardiovascular branch of the Chinese medical association, etc. guidelines for the diagnosis and treatment of stable coronary heart disease. Chin J Cardiovascu Dis. 2018;46(9):680–694.

16. Gallone G, Bellettini M, d’Ascenzo F. Coronary plaque characteristics associated with major adverse cardiovascular events in atherosclerotic patients and lesions: a systematic review and meta-analysis. JACC Cardiovasc Imaging. 2023;16(12):1584–1604. doi:10.1016/j.jcmg.2023.08.006

17. Wang J-X, Liang M-M, Gao J. Using machine learning to predict MACEs risk in patients with premature myocardial infarction. Rev Cardiovasc Med. 2025;26(5):31298. doi:10.31083/RCM31298

18. Li X, Zhao Y, Hu D. Development of an interpretable machine learning model associated with heavy metals’ exposure to identify coronary heart disease among US adults via SHAP: findings of the US NHANES from 2003 to 2018. Chemosphere. 2023;311(Pt 1):137039. doi:10.1016/j.chemosphere.2022.137039

19. Hou XZ, Wu Q, Wang SH. Development and external validation of a risk prediction model for depression in patients with coronary heart disease. J Affect Disord. 2024;367:137–147. doi:10.1016/j.jad.2024.08.218

20. Chen B, Ruan L, Li T. Machine learning improves risk stratification of coronary heart disease and stroke. Ann Transl Med. 2022;10(21):1156. doi:10.21037/atm-22-1916

21. Yudistira KA, Yasmin, Wiryawan N. Comparison between 6 minutes walking distance (6MWD), changes in rate pressure product (RPP) and heart rate walking speed index (HRWSI) of 6 minutes walk test (6MWT) pre-discharge as a predictor major adverse cardiac event (MACE) in patients with acute heart failure. Inter J Sci Adv. 2023;4(5):676–687.

22. Chen C, Wang N, Yu Y, et al. Association of echocardiographic cardiac valve calcification with major adverse cardiovascular events in patients on acute myocardial infarction. medRxiv. 2024;2024.10.15.24315566.

23. Rong Y, Zheng H, Zhichao G, et al. Study on the predictive value of serum omentin-1 for major in-hospital adverse cardiovascular events after percutaneous coronary intervention in patients with ST-segment elevation myocardial infarction. J Prac Cardiovascu Cerebrovasc Pulm Dis. 2021;29(12):14–18.

24. Knuuti J, Wijns W, Bax JJ, ESC Scientific Document Group. 2019 ESC guidelines for the diagnosis and management of chronic coronary syndromes. Eur Heart J. 2020;41(3):407–477. doi:10.1093/eurheartj/ehz425

25. Collet JP, Thiele H, Siontis GCM. ESC scientific document group. 2020. ESC guidelines for the management of acute coronary syndromes in patients presenting without persistent ST-segment elevation. Eur Heart J. 2021;42(14):1289–1367. doi:10.1093/eurheartj/ehaa575

26. McDonagh TA, Metra M, Francesco Piepoli M, et al. ESC Scientific Document Group. 2021 ESC guidelines for the diagnosis and treatment of acute and chronic heart failure. Eur Heart J. 2021;42(36):3599–3726. doi:10.1093/eurheartj/ehab368

27. Nicholls SJ, Puri R, Nissen SE. Effect of evolocumab on progression of coronary disease in statin-treated patients: the GLAGOV randomized clinical trial. JAMA. 2016;316(22):2373–2384. doi:10.1001/jama.2016.16951

28. Sandesara PB, Virani SS, Shapiro MD, Shapiro MD. The forgotten lipids: triglycerides, remnant cholesterol, and atherosclerotic cardiovascular disease risk. Endocr Rev. 2019;40(2):537–557. doi:10.1210/er.2018-00184

29. Wenjie X, Liangliang H, Yuanyuan L, et al. Analysis of tongue color characteristics in Framingham populations with different cardiovascular risk levels and patients with coronary atherosclerotic heart disease. Shanghai J Traditional Chin Med. 2020;54(09):23–26.

30. Yang J. Research on the correlation between the characteristics of TCM syndromes in the stable stage of coronary heart disease and recurrent cardiovascular events. Beijing Univ Chin Med. 2014;13:806300.

31. Qian L. Establishment and validation of a prognostic evaluation model of integrated traditional Chinese and western medicine for elderly patients with stable coronary heart disease. China Acad Chin Med Sci. 2020.

32. Gongyu L. Research on the integrated traditional Chinese and western medicine prediction model for stable coronary heart disease cardiovascular events based on machine learning. Beijing Univ Chin Med. 2023.

33. Chen J. Study on the correlation between cyanosis tongue and the risk level of patients with acute myocardial infarction. Chongqing Med Univ. 2022.

34. Ming W, Lifu M, Ruilin Z, et al. Study on the relationship between purple-dark tongue and prognosis of patients with myocardial infarction. J Beijing Univ Chin Med. 2012;35(07):494–496+504.

35. Yunhu C, Moqing Y, Hongfeng XU. Mirror-like tongue is an important predictor of acute heart failure: a cohort study of acute heart failure in Chinese patients. J Tradit Chin Med. 2023;43(6):1243–1251. doi:10.19852/j.cnki.jtcm.20230904.004

36. Wang XY, Zhang F, Yang J, Zheng L-R, Yang J. The biomarkers for acute myocardial infarction and heart failure. Biomed Res Int. 2020;2020(1):2018035. doi:10.1155/2020/2018035

37. Subkhan M, Heriansyah T, Dimiati H, Purnawarman A, Dimiati H. Association between NT-proBNP level and the number of stents with major advanced cardiovascular events (MACE) in patients with multivessel coronary artery disease treated with percutaneous coronary intervention: a prospective cohort study. Narra J. 2024;4(1):e710. doi:10.52225/narra.v4i1.710

38. Vernooij LM, van Klei WA, M Damen JA, Takada T, van Waes J, Damen JA. The comparative and added prognostic value of biomarkers to the revised cardiac risk index for preoperative prediction of major adverse cardiac events and all-cause mortality in patients who undergo noncardiac surgery. Cochrane Database Syst Rev. 2021;12(12):CD013139. doi:10.1002/14651858.CD013139.pub2

39. Zelniker TA, Wiviott SD, Morrow DA. Association of cardiac biomarkers with major adverse cardiovascular events in high-risk patients with diabetes: a secondary analysis of the DECLARE-TIMI 58 trial. JAMA Cardiol. 2023;8(5):503–509. doi:10.1001/jamacardio.2023.0019

40. Xiang X, Zhao K, Zhao S. Prediction of adverse outcomes in nonischemic dilated cardiomyopathy: a CMR-based nomogram. Int J Cardiol. 2023;390:131136. doi:10.1016/j.ijcard.2023.131136

41. Kraler S, Liberale L, Nopp S, Moik F. Biomarker-enhanced cardiovascular risk prediction in patients with cancer: a prospective cohort study. J Thromb Haemost. 2024;22(11):3125–3136. doi:10.1016/j.jtha.2024.07.019

42. Sijie Q, Song G, Ben H. Analysis of the prognostic predictive value of different health-related indicators for elderly patients with diabetes mellitus complicated with coronary heart disease. Public Health Prev Med. 2023;34(05):103–106.

43. Zhang YJ, Sun Y, Zhao YB, Ma D. α-HBDH is a superior to LDH in predicting major adverse cardiovascular events in patients with acute aortic dissection. Heliyon. 2024;10(8):e29155. doi:10.1016/j.heliyon.2024.e29155

44. Xiaomin Q, Sirong L, Hui L, et al. Relationship between apoB/apoA-1, LDH and ALP levels and major adverse cardiovascular events in patients with coronary heart disease. J Molecular Diagnosis Therap. 2021;13(10):1635–1638+1643.

45. Lv D, Guo Y, Li X, Zhang L. Increased transferase ratio is associated with adverse cardio-cerebral events in patients with unstable angina: a retrospective cohort study. Medicine. 2023;102(31):e34563. doi:10.1097/MD.0000000000034563

46. Xinyi H, Xudong Y, Shikui S. The predictive value of FAI and CT-FFR combined with Gensini score for the occurrence of MACE after percutaneous coronary intervention for ACS. Chin J CT and MRI. 2025;23(04):90–92+135.

47. Zhao K, Zhu Y, Zheng H. Machine learning in hypertrophic cardiomyopathy: nonlinear model from clinical and CMR features predicting cardiovascular events. JACC Cardiovasc Imaging. 2024;17(8):880–893. doi:10.1016/j.jcmg.2024.04.013

48. Juan-Salvadores P, Veiga C, Romo AI. Using machine learning techniques to predict MACE in very young acute coronary syndrome patients. Diagnostics. 2022;12(2):422. doi:10.3390/diagnostics12020422

49. Tang CX, Guo BJ, Zhang LJ. Feasibility and prognostic role of machine learning-based FFR_CT in patients with stent implantation. Eur Radiol. 2021;31(9):6592–6604. doi:10.1007/s00330-021-07922-w

50. Alahdab F, El Shawi R, Al-Mallah M, Han Y, Al-Mallah M. Patient-level explainable machine learning to predict major adverse cardiovascular events from SPECT MPI and CCTA imaging. PLoS One. 2023;18(11):e0291451. doi:10.1371/journal.pone.0291451