Skip to main content

ORIGINAL RESEARCH article

Front. Cardiovasc. Med., 08 February 2022
Sec. Atherosclerosis and Vascular Medicine

Thrombosis and Major Bleeding Risk After Primary PCI Among Patients With Multivessel Coronary Artery Disease

\nXiaoxiao ZhaoXiaoxiao Zhao1Chen LiuChen Liu1Peng ZhouPeng Zhou1Zhaoxue ShengZhaoxue Sheng1Jiannan LiJiannan Li1Jinying ZhouJinying Zhou1Runzhen ChenRunzhen Chen1Ying WangYing Wang1Yi ChenYi Chen1Li SongLi Song1Hanjun Zhao
Hanjun Zhao1*Hongbing Yan
Hongbing Yan2*
  • 1Department of Cardiology, National Center for Cardiovascular Diseases, Peking Union Medical College, Fuwai Hospital, Chinese Academy of Medical Sciences, Beijing, China
  • 2Fuwai Hospital, Chinese Academy of Medical Sciences, Shenzhen, China

Background and Aim: This study aimed to develop and validate separate risk prediction models for thrombosis events (TEs) and major bleeding (MB) in patients with multivessel coronary artery lesions who had undergone primary percutaneous coronary intervention (PCI).

Methods and Results: Thrombosis events (TEs) were defined as the composite of myocardial infarction recurrence or ischemic cerebrovascular events, whereas MB was defined as the occurrence of bleeding academic research consortium (BARC) three or five bleeding. The derivation and validation cohorts comprised 2,976 patients who underwent primary PCI between January 2010 and June 2017. At a median follow-up of 3.07 years (1,122 days), TEs and MB occurred in 167 and 98 patients, respectively. Independent predictors of TEs were older age, prior PCI, non-ST elevated MI (NSTEMI), and stent thrombosis (ST). Independent predictors of MB were triple therapy at discharge, coronary artery bifurcation lesions, lesion restenosis, target lesion of the left main coronary artery, stent thrombosis, non-use of IABP during primary PCI, type A/B according to the American College of Cardiology classification of the coronary lesion, and PTCA. In the derivation and validation cohorts, the areas under the curve were 0.817 and 0.82 for thrombosis and 0.886 and 0.976 for bleeding, respectively. In the derivation cohort, high thrombotic risk (n = 755) was associated with higher 3-year incidence of TEs, major adverse cardiovascular events (MACEs), and all-cause death compared to low risk (n = 1,275) (p = 0.0022, 0.019, and 0.012, respectively). High bleeding risk (n = 1,675) was associated with higher incidence of bleeding, MACEs, and cardiac death compared to low risk (n = 355) (p < 0.0001).

Conclusion: Simple risk scores can be useful in predicting risks of ischemic and bleeding events after primary PCI, thereby stratifying thrombotic or MB risks and facilitating clinical decisions.

Introduction

The risk of thrombotic events (TEs), such as myocardial infarction (MI) and stent thrombosis (ST), is lower in patients with the acute coronary syndrome who receive dual antiplatelet therapy (DAPT) with aspirin or who have undergone primary percutaneous coronary intervention (PCI) (1). The marginally higher mortality associated with bleeding events has been reported to be comparable to the risk associated with MI (2, 3). Hence, avoiding bleeding events is becoming increasingly important. These results suggest that clinical decision-making concerning the optimal duration of DAPT for individual patients following PCI must be predicated on balancing the long-term risks of intensive antithrombotic therapy and avoiding major bleeding (MB) (4, 5). In this context, it is essential to develop stratification tools for distinguishing high-risk ischemic patients from high-risk bleeding patients. To date, most post-PCI algorithms are single scoring systems or focused on in-hospital events or short-term risks (611). Although it is practicable to inform clinical decisions with respect to the short-term provision of DAPT, respective weights of underlying risk factors vary from early as opposed to later, and predicting risks of thrombotic and bleeding events is more important in the long term. Accordingly, this study aimed to explore prediction rules for long-term outcomes of TEs and MB events separately in a large Chinese observational database of patients with multivessel coronary artery lesions who had undergone primary PCI. We prepared and presented the current article in accordance with the TRIPOD reporting checklist (Appendix File)1.

Materials and Methods

Study Population: Enrollment and Randomization

A prospective observational study was conducted on patients who had undergone primary PCI in Fuwai Hospital (National Center for Cardiovascular Diseases, Peking Union Medical College and Chinese Academy of Medical Sciences) in Beijing, China, between January 2010 and June 2017. This study was designed to investigate the validity of separate risk prediction models for subsequent clinical adverse events. With respect to eligibility criteria, adult patients (1) who had undergone primary PCI, including stent implantation, thrombus aspiration, and balloon dilation in the coronary artery, and (2) who provided written informed consent were included in the study. Patients (1) who refused participation, (2) who were lost to follow-up when contacted, and (3) who did not have coronary angiography parameters (coronary angiography, used to diagnose ischemic heart disease after chest pain, is a procedure that uses contrast dye and x-ray pictures to detect blockages in coronary arteries; the coronary angiography parameters mentioned above including whether patients used thrombus aspiration, stent implantation, use of IABP, percutaneous transluminal coronary angioplasty, and complex procedure) or multivessel coronary artery lesions were, however, excluded from the analysis. Following the application of the inclusion and exclusion criteria, a total of 3,976 subjects with acute myocardial infarction with multivessel lesions remained. The patients were administered aspirin 300 mg, clopidogrel 600 mg or ticagrelor 180 mg, and heparin 100 IU/kg before the procedure of intervention. The access of primary PCI was performed via radial or femoral artery. Thrombus aspiration was performed to reduce the burden of the thrombus. Duration of dual antiplatelet therapy (DAPT) consisted of oral aspirin and a P2Y12 inhibitor for at least 12 months following primary PCI.

The included patients were randomly and proportionally (70:30%) divided into the derivation cohort (n = 2,084) and the validation cohort (n = 892) (divide the training dataset into training and validation sets, ideally 7,030, and model on 70% of the training dataset; then, use the 30% validation data set for cross-validation and performance evaluation using evaluation metrics).

This study was conducted according to the principles outlined in the Declaration of Helsinki and was approved by the Ethics Committee of Fuwai Hospital. All the study subjects gave informed consent.

Study Definitions

Thrombosis events (TEs) were defined as the occurrence of coronary thrombotic complications such as ischemic cerebrovascular events or MI recurrence. Consistent with the universal definition, MI recurrence was defined as the recurrence of chest pain accompanied by either re-ST-segment elevation or ST-segment depression attributed to myocardial ischemia and re-elevation of cardiac troponin I >25% (12). MB was defined as the occurrence of Bleeding Academic Research Consortium type three or five bleeding (13) and was adjusted by a blinded committee. Multivessel coronary artery lesions, characterized by significant stenosis (diameter ≥1.5 mm and significant stenosis >50%) in all three major coronary arteries determined by a cardiologist. Stroke was defined as rapidly developing focal or general brain dysfunction that lasted for more than 24 h or caused death, excluding non-vascular causes (e.g., trauma, metabolic disorders, tumors, and any neurological abnormalities due to central nervous system infection). Additionally, ischemic stroke included cerebral thrombosis and cerebral embolism. Imaging data are as follows: computed tomography (CT) showed insular signs, namely, disappearance and blurring of the gray matter interface in the conduction zone, consistency of the density of the insular cortex with the outer capsule, and disappearance or narrowing of the cerebral sulci in the cortex; the abnormal high signal shadow was found in the responsible lesion area under magnetic resonance imaging (MRI) detection; major adverse cardiovascular events (MACEs) were identified as the composite of overall mortality, MI recurrence, and ischemic cerebrovascular events; events were identified using physician-reported diagnoses extracted from cardiac catheterization laboratory report, hospital discharge records, or clinical notes in the event of death. Anemia was defined as hemoglobin level <12 g/dl in men and <11 g/dl in women (14). Complex PCI procedures included bifurcation, total occlusion, thrombus, or >2 stent. Chronic kidney disease (CKD) was defined as chronic renal structural and functional impairments due to a variety of causes such as eGFR lesions (<60 ml/min·1.73 m2), imaging abnormalities, and abnormal blood/urine composition for at least the past 3 months. The estimated glomerular filtration rate (eGFR) was calculated by the formula of Modification of Diet in Renal Disease (MDRD). No-reflow phenomenon was defined as Thrombolysis in Myocardial Infarction flow grade <3 after primary PCI. The patients were subsequently divided into the low-risk and high-risk groups according to the best cutoff value of the prognostic index. Triple therapy on discharge was defined as a combination of DAPT (aspirin plus thienopyridine) and oral anticoagulation therapy.

Follow-Up

The patients were followed-up at least 1 year after discharge by physicians. The health status of the enrolled patients was confirmed via telephone calls and review of health records, and this method, which was approved by the Review Board of Fuwai Hospital. The follow-up primary endpoints including TEs, MB, MACEs, cardiac death, and all cause death, were identified and extracted primary endpoints from recordings of hospital records, laboratory reports, and clinical notes in the event of death by physicians who in charge of follow-up.

Statistical Analysis

Normal distribution of outcome variables was conducted using the method of Kolmogorov–Smirnov test. Baseline clinical and procedural characteristics of primary PCI according to the presence or absence of TEs and MB were compared between patients with and without multivessel lesions, and between the derivation and validation cohorts by chi-square test and Student's t-test for categorical variables and continuous variables, respectively. Cox proportional hazards regression and least absolute shrinkage and selection operator (LASSO) regression generated separate prediction models to achieve reduction and simplification of the models and to prevent the occurrence of overfitting, with time to the first occurrence of TEs or MB serving as the dependent variable in each respective model. Event-free patients were censored at the end of the study, or at the time of death or last contact, whichever came first.

To account for missing data on serum low-density lipoprotein cholesterol at baseline (n = 9, 0.4%), history of smoking (n = 192, 9.2%), body mass index (n = 114, 5.5%), and prior hyperlipidemia (n = 210, 10.1%), covariates for each model were identified through an iterative process involving multiple imputation with automated variable selection. Five imputations were generated from the original dataset, and multivariate normal regression was used to substitute the missing data on serum low-density lipoprotein cholesterol, history of smoking, body mass index, and prior hyperlipidemia within each impute in the first step.

Candidate covariates for each model were age, sex, triple therapy at discharge, body mass index (<25, 25–35 [reference], >35 kg/m2), smoking status, prior hyperlipidemia, hypertension, diabetes mellitus, prior coronary artery bypass grafting, prior PCI, CKD, abnormal liver function, malignancy, Killip classification, length of the lesion, complex procedure (bifurcation, total occlusion, thrombus, or >2 stents), anemia, eGFR <60 ml/min, platelet count <1,00,000/μl, percutaneous transluminal coronary angioplasty, the target of bifurcation, use of intra-aortic balloon pump (IABP), stent implantation, no-reflow phenomenon, revascularization after discharge, lesion restenosis, and target lesion of the left main, right, left circumflex, or left anterior descending coronary artery. Finally, one imputes with regression coefficients combined across all imputed datasets, as described by Rubin, was generated using the covariates (15). Variables that remained significant at a threshold of min Se were retained as final predictors by LASSO regression. Model discrimination was quantified using Harrell's c-statistic and calibration chart for the derivation and validation cohorts. At the beginning of the model establishment of LASSO regression, all identified independent variables were selected to minimize model deviation caused by the non-inclusion of significant independent variables, which were selected by univariable regression. Therefore, multivariable regression was not conducted. To improve prediction accuracy, the established model needs to find the set of independent variables with the strongest explanatory power for the dependent variables. A more refined model is obtained by constructing a function that compresses some coefficients and sets some coefficients to 0, 0.5, or minimization. LASSO regression is a biased estimation of data with complex collinearity retaining the advantage of contraction. The Lars algorithm software package of R language provided LASSO programming. Hence, variable selection and dimensionality reduction can be achieved accurately by LASSO regression.

Integer risk scores for the outcomes of TEs and MB were generated using fully adjusted regression coefficients, as described by Sullivan et al. (16). The survival receiver operating characteristic (ROC) curve for prognostic index (PI) was fitted using the Kaplan–Meier method, and the best cutoff value was obtained. The patients were subsequently divided into the low-risk and high-risk groups according to the best cutoff value of the prognostic index (PI). Observed event rates were calculated as Kaplan–Meier (K-M) estimates of time to the first event. Predicted event rates were estimated using fully adjusted Cox regression models. The main software used for statistical analysis in this study used survival and rms package in R language version I 386 3.6.3. Other analyses were performed using SPSS version 20.0 (IBM Corp., Armonk, NY, United States). All p-values were two-tailed, and statistical significance was set at p < 0.05.

Performance and Internal Validation Cohort

External validation of each score was performed. Each subject in the validation cohort was assigned with a TE risk score and MB event risk score in the same manner as in the derivation cohort. The patients were subsequently categorized into subjects with low and high thrombotic and bleeding risks using the same thresholds as in the derivation cohort. The 3-year adverse event rates were counted for each risk division using the K-M method, and the difference was determined by the method of the log-rank test. Discriminations of both the derivation and validation cohorts were assessed by calculating the area under the ROC curve (AUC) and expressed as c-statistic with the use of the MedCalc software for Windows, version 18.2.1.0 (MedCalc Software, Mariakerke, Belgium). The accuracy of the new model and acuity (bleeding model)/autar (thrombosis model) risk score model predicting MACEs among patients with MI who underwent PPCI was compared according to the area under the ROC (AUC) curve by a non-parametric test developed. MedCalc for Windows version 18.2.1 (MedCalc Software, Mariakerke, Belgium) was used for comparison.

Results

Baseline Characteristics

Among 4,151 enrolled patients who had undergone PCI in Fuwai Hospital (Beijing, China) between January 2010 and June 2017, those who were lost to follow-up (n = 97), had no coronary angiography parameters (n = 78), and had no multivessel coronary artery lesions (n = 1,000) were excluded. Hence, the study population comprised 2,976 patients in total. At a median follow-up of 3.07 years (1,122 days), 167 patients sustained TEs, whereas 98 patients experienced MB.

Of the patients suffering MB, the incidence of cerebral hemorrhage, fundus bleeding, gastrointestinal bleeding, urogenital bleeding, nasal mucosa bleeding was 11.22 (11), 19.39 (19), 53.06 (52), 7.14 (7), and 9.18% (9), respectively. The baseline characteristics are summarized in Tables 1.1, 1.2. Patients with TEs more frequently presented with MACEs (p < 0.001) and showed a higher incidence of hyperlipidemia (p = 0.017) than those without TEs. Patients with MB had lower degree of lesion stenosis (p = 0.034) and showed higher incidence of all-cause death (p = 0.033), cardiac death (p = 0.041), and cerebral hemorrhage (p < 0.001) than their counterparts without MB.

TABLE 1.1
www.frontiersin.org

Table 1.1. Baseline clinical characteristics of patients with vs. without thrombotic events or major bleeding events.

TABLE 1.2
www.frontiersin.org

Table 1.2. Baseline clinical characteristics of subgroup patients with vs. without thrombotic events or major bleeding. events.

There was no significant difference in baseline characteristics between the derivation cohort and the validation cohort with respect to several aspects, such as clinical variables, laboratory examination, and procedural characteristics (Table 1.2). For the purpose of assessing the strength of association between 39 potential predictors and TEs or MB events in the derivation cohort, we constructed univariate Cox regression models and presented the results as hazard ratios (95% CI) with p-values. Variables identified to show association (with p < 0.5) in the univariate logistic regression models were included in the multivariate models (Table 2).

TABLE 2
www.frontiersin.org

Table 2. Univariate Cox analysis for major bleeding events and thrombotic events in the derivation cohort.

Predictors of Thrombotic Risk Scores and MB Risk Scores

Point estimates for each predictive covariate in the final prediction models are shown in Tables 3.1, 3.2. The strongest predictors of TEs, quantified and ranked using the values in the nomogram, were older age, history of prior PCI, non-ST-elevation MI, and stent thrombosis status. Correspondingly, the strongest contributors to MB were percutaneous transluminal coronary angioplasty, non-use of IABP during primary PCI, target lesion of the left main coronary artery, triple therapy at discharge, lesion restenosis, and type A/B according to the American College of Cardiology classification of coronary lesions, without stent thrombosis and coronary artery bifurcation lesions. The coronary TE prediction model had a moderate level of discrimination, with a c-index of 0.616 and adequate calibration for the entire population. Analogous parameters of model performance for MB in the entire cohort had a c-statistic of 0.676. Using the fully adjusted regression coefficients and nomogram graph, we developed integer-based risk projects for both MB and TEs (Tables 3.1, 3.2, respectively).

TABLE 3.1
www.frontiersin.org

Table 3.1. Integer risk score for major bleeding.

TABLE 3.2
www.frontiersin.org

Table 3.2. Integer risk score for thrombotic events.

Furthermore, 627 (83%) out of 755 subjects with high TEs risk scores in the derivation cohort and 258 (51.7%) out of 499 subjects with high TEs risk scores in the validation cohort also had high MB risk scores; mortality and MB rates for these subjects were very high (Figure 1). Among those with high thrombotic risk scores, only 128 patients (17.0%) in the derivation cohort and 241 patients (48.3%) in the validation cohort had low bleeding risk scores (Figure 1). The majority of patients with low thrombotic risk scores had high bleeding risk scores (Figure 1). Within each thrombotic risk level, the frequency of high MB risk increased in the majority of patients as the thrombotic risk increased.

FIGURE 1
www.frontiersin.org

Figure 1. Distribution of bleeding risk score categories according to thrombotic risk score categories in the (A) derivation cohort (N = 2,030, p = 0.336) and (B) validation cohort (N = 892, p = 0.275). Blue represents high bleeding risk, and red represents low bleeding risk.

Clinical Outcomes of Thrombosis and Clinically Relevant Bleeding in the Derivation and Validation Cohorts

The 3-year cumulative incidences of TEs/MB, MACEs, cardiac death, and all-cause death according to the thrombotic risk score categories and MB risk score categories in the derivation and validation cohorts are presented in Figures 2A–D. In the derivation cohort, TEs (p = 0.0022), MACEs (p = 0.019), and all-cause death (p = 0.012) were significantly different between the low-risk and high-risk groups divided by the PI of thrombotic risk score categories (Figure 2A). Analogously, in the validation cohort, TEs (p = 0.023), MACEs (p = 0.00057), cardiac death (p = 0.024), and all-cause death (p = 0.00057) were significantly different between the low-risk and high-risk groups divided by the PI of thrombotic risk score categories (Figure 2B). Similar results are presented in Figures 2C,D.

FIGURE 2
www.frontiersin.org

Figure 2. (A) Cumulative 3-year incidence of thrombotic events (p = 0.0022), MACEs (p = 0.019), cardiac-caused death (p = 0.076), and all-caused death (p = 0.012) according to the thrombotic risk score categories in the derivation cohort (N = 2,084). (B) Cumulative 3-year incidence of thrombotic events (p = 0.023), MACEs (p = 0.00057), cardiac-caused death (p = 0.0024), and all-caused death (p = 0.00057) according to the thrombotic risk score categories in the validation cohort (N = 892). (C) Cumulative 3-year incidence of major bleeding events (p = 0.0003), MACEs (p < 0.0001), cardiac-caused death (p < 0.0001), and all-caused death (p < 0.0001) according to the thrombotic risk score categories in the derivation cohort (N = 2,084). (D) Cumulative 3-year incidence of major bleeding events (p = 0.00039), MACEs (p < 0.0001), cardiac-caused death (p = 0.002), and all-caused death (p = 0.00041) according to the thrombotic risk score categories in the validation cohort (N = 892). MACEs, major adverse cardiovascular events; PI group, prognostic index group; PI group = 0 in Figures 3A,B, low risk group classified by the PI of thrombotic risk score categories; PI group = 1 in Figures 3A,B, high-risk group classified by the PI of thrombotic risk score categories; PI group = 0 in Figures 3C,D, low-risk group classified by the PI of major bleeding risk score categories; PI group = 1 in Figures 3C,D, high risk group classified by the PI of major bleeding risk score categories.

Figure 3 shows the ROC curves for the discriminatory value of the 3-year evaluation performance of the risk prediction model in the derivation and validation cohorts. The AUCs by the TE prediction project were 0.817, 0.771, 0.927, and 0.893 for TEs, MACEs, cardiac death, and all-cause death, respectively, in the derivation cohort (Figure 3A). The AUCs by the TE prediction project were 0.820, 0.782,0.973, and 0.906 for TEs, MACEs, cardiac death, and all-cause death, respectively, in the validation cohort (Figure 3B). The AUCs by the MB prediction project were 0.886, 0.791, 0.939, and 0.906 for bleeding, MACEs, cardiac death, and all-cause death, respectively, in the derivation cohort (Figure 3C). The AUCs by the MB prediction project were 0.976, 0.86, 0.863, and 0.961 for bleeding, MACEs, cardiac death, and all-cause death, respectively, in the validation cohort (Figure 3D).

FIGURE 3
www.frontiersin.org

Figure 3. (A) ROC curve of thrombotic events (p < 0.001), MACEs (p < 0.001), cardiac-caused death (p < 0.001), and all-caused death (p < 0.001) according to the thrombotic risk score categories in the derivation cohort (N = 2,084). (B) Cumulative 3-year incidence of thrombotic events (p < 0.001), MACEs (p < 0.001), cardiac-caused death (p < 0.001), and all-caused death (p < 0.001) according to the thrombotic risk score categories in the validation cohort (N = 892). (C) Cumulative 3-year incidence of major bleeding events (p < 0.001), MACEs (p < 0.001), cardiac-caused death (p < 0.001), and all-caused death (p < 0.001) according to the major bleeding risk score categories in the derivation cohort (N = 2,084). (D) Cumulative 3-year incidence of major bleeding events (p < 0.001), MACEs (p < 0.001), cardiac-caused death (p < 0.001), and all-caused death (p < 0.001) according to the major bleeding risk score categories in the validation cohort (N = 892). MACEs, major adverse cardiovascular events.

Evaluation of the Risk Prediction Model

Appendix Figure 1 illustrates the evaluation performed, including the calibration curve (Appendix Figures 1A,B) and decision curve analysis (DCA) curve (Appendix Figure 1C) for TEs and MB, respectively, in the derivation and validation cohorts. It is comparable between the observed and predicted risks projects, illustrating that the model calibration was excellent for both predicting scores. Figure 4A shows a pairwise comparison of ROC curves between the new bleeding model and acuity risk score model. The AUC of the new bleeding model is 0.743, and the AUC of acuity score is 0.721. Figure 4B shows a pairwise comparison of ROC curves between the new thrombotic model and autar risk score model. The AUC of the new thrombotic model is 0.818, and the AUC of acuity score is 0.829.

FIGURE 4
www.frontiersin.org

Figure 4. (A). Pairwise comparison of ROC curves between the new bleeding model and the acuity risk score model (N = 3,976). The area under the ROC curve (AUC) of the new bleeding model is 0.743 [confidence interval (CI)0.723–0.763], and the AUC of acuity score is 0.721 (CI.706–0.747). P = 0.4047, z statistic, 0.833. (B). Pairwise comparison of ROC curves between the new thrombotic model and the autar risk score model (N = 3,976). The AUC of new the thrombotic model is 0.818 (CI.803–0.833), and the AUC of acuity score is 0.829 (0.814–0.843). p = 0.1079, z statistic, 1.608.

Discussion

This study, which involved 2,976 real-world patients with multivessel coronary artery disease who had undergone primary PCI in China, yielded the following main findings: first, we reported the development of separate models for predicting the risks of TEs and MB, which demonstrated moderate accuracy in discrimination concordant and stratified the risk in the derivation and validation cohorts; second, we showed that subjects with high thrombotic risk also had high bleeding risk in a large proportion of the study population.

The objective of this study was to identify readily available characteristics that were independently correlated with TEs and MB in patients with multivessel lesions who had undergone primary PCI. The study described risk indicators (e.g., clinical characteristics, angiography). While various tools have been developed to stratify risk after undergoing PCI, the majority of them are focused on peri-procedural short-term outcomes (69, 11, 1719). We modeled events occurring after discharge and found that clinical risk factors and procedural parameters could predict the risk of TEs. This result is consistent with the findings of previous studies that emphasized the importance of the complexity of coronary artery lesions with respect to the risk of TEs (6, 20). Nevertheless, this result differs from the finding of a previous study that highlighted the importance of clinical risk factors alone (21). These discrepancies may be because thrombotic risk factors are not static but dynamic over time. A thrombotic risk score was proposed in the Thrombin Receptor Antagonist in Secondary Prevention of Atherothrombotic Ischemic Events–Thrombolysis in Myocardial Infarction 50 trial (22), which included ischemic stroke as a TE; age was the common independent predictor of TEs in that previous study. According to the most recent ESC guidelines for the management of patients presenting with NSTEMI, based on the result of the ISAR-REACT 5 trial, prasugrel is the recommended drug in patients who proceed to PCI. No patient in this study assumed prasugrel, because this was not the antithrombotic drug of choice during the period in which the patients were enrolled. Analogously, we identified older age as an independent determinant of long-term TEs, which is consistent with the previously observed association between parameters and thrombosis. It is plausible to include ischemic stroke as one of the components of the TE composite endpoint. Among subjects who underwent PCI, ischemic stroke, which demands intensive antithrombotic therapy, is as clinically important as MI and stent thrombosis.

The risk factors constituting the bleeding risk score established in our study were consistent with those in previous studies generally (2, 7, 23). Not surprisingly, we identified triple therapy [defined as a combination of oral anticoagulant therapy and DAPT (aspirin plus thienopyridine)] at discharge as an independent determinant of long-term bleeding, which is consistent with previous observations (21). Nonetheless, female sex, older age, and renal dysfunction were absent in our final bleeding model. It is possible that earlier studies that focused on shorter-term events accounted for these discrepancies and that underlying risk factors for bleeding were not constant but variable over time resulting in differences. Analogously, the hypothesis is similar to the findings of Genereux et al. (2) and Ducrocq et al. (24), who reported that bifurcation lesions were associated with post-discharge bleeding. However, no prospective study has indicated clinical utility to guide treatment decisions. The use of proton pump inhibitors has been shown to reduce the incidence of gastrointestinal bleeding in patients undergoing PCI, and liberal prophylactic use is essential for the prognosis of these patients. The present prediction project assessing bleeding events was generally consistent with previous studies, and the difference in the risk factors identified might be attributable to the selected population, race, and study design. Therefore, it is preferable to practice the prediction rule derived from the cohort with homologous characteristics of the race. The present prediction project evaluating MB risks showed modest accuracy in both the derivation and validation cohorts. AUCs ranged from 0.791 to 0.976, indicating that the risk score was helpful for discrimination in the clinical prediction models. We performed a pairwise comparison of the new bleeding model with the acuity risk score model, and the new thrombotic model with the autar risk score model. The AUC of the new bleeding model is larger than the acuity score, which showed excellent performance. Although the AUC of the new thrombotic model is similar to the autar risk score, it still simplified the model of predicting thrombosis events and is of benefit for application in clinical assessment.

It is recommended and well-validated by guidelines that risk stratification tools assist with therapeutic decision-making for patients with acute coronary syndrome (2528). This study included more complex patients than previous trials, that is, this study included patients with two-vessel and triple-vessel diseases. In the context of the growing trend toward individualization and evidence-based therapy, risk stratification could meet patient preferences and enhance patient compliance while balancing against the adverse effects of some therapies (e.g., thrombosis and bleeding) in patients with multivessel disease. The risk stratification strategy outlined in this study provides clinic doctors with an opportunity to select potential candidates with the greatest absolute gains, and it is important to offer therapeutic interventions for secondary prevention of acute MI. We did not include the type of stent used, because previous trials had reported no significant discrepancy in the incidence of TEs between bare-metal stents and first-generation drug-eluting stents (29, 30).

Out-of-hospital stroke results in substantial mortality and morbidity. Considering that substantial mortality and morbidity are correlated with post-PCI ischemic stroke, more studies evaluating risk factors are required to prevent post-PCI ischemic stroke. The Organization to Assess Strategies for Ischemic Syndromes I and II studies (31) reported that stroke in subjects with the coronary disease was correlated with a 6-month mortality rate of 27%. Therefore, we assessed the incidence of stroke in patients who underwent primary PCI. Ischemic stroke and hemorrhagic stroke were categorized as a TE and an MB event, respectively. IABP use and cerebral hemodynamic impairment, especially ischemic stroke, have potential associations. This study assessed risk factors for TEs after primary PCI; however, no association was identified. According to the most recent ESC guideline for the management of patients presenting with NSTEMI, triple antithrombotic therapy is suggested only for 1 week; this recent recommendation has probably reduced the incidence of bleeding.

In this study, we have found that a lower degree of lesion stenosis contributed to MB. After reviewing the literature, we did not find a reasonable explanation for this result. Therefore, basic research on this aspect should be carried out accordingly in the future to carry out relevant research from the perspective of the mechanism. The study of Marco Zimarino et al. (32) has made the conclusion that PCI of bifurcation lesions is associated with increased risk of thrombotic events and investigated the theme of the duration of DAPT after PCI of bifurcation lesion. The literature has exposed the state of art concerning of the choice of antithrombotic drugs, timing of initiation, the DAPT duration, risk stratification and overall the identification of patients at high bleeding risk, with a decision-making algorithm for DAPT duration in PCI in coronary bifurcation.

The PRECISE-DAPT score, which has been validated in two large independent patient populations with ACS, is a five-item bleeding risk prediction model developed to estimate the bleeding risk in patients who receive DAPT after stent implantation (4). Based on the PRECISE-DAPT score, categorization of patients has been proved to be valuable to inform decision-making for the duration of DAPT in stented patients (33, 34). There are two retrospective analyses (35, 36) that showed an absolute bigger reduction of ischemic risk in patients who are receiving long-term DAPT after the complex intervention. Furthermore, patients who have undergone complex intervention simultaneously carry features that greatly increase their bleeding risk, including renal disease, multiple comorbidities, and previous bleeding.

Study Limitations

This study has several important limitations that should not be ignored. First of all, it has an observational and prospective design, which precludes causal inference and carries inherent limitations. Second, information on the history of previous bleeding events, which could be significantly correlated with exceedingly great risk for MB events, has not been collected. Third, potential reporting bias might have been introduced, accounting for the fact that bleeding events were not adjudicated by a blinded clinical event committee independently. Fourth, the study is derived between 2010 and 2017, and early-generation drug-eluting stents are not distinguished from new-generation drug-eluting stents. Fifth, follow-up information is prospectively recorded according to prespecified definitions, which may limit the power to identify other predictors of stroke. Additionally, genetic characteristics and unidentified biochemical parameters should be considered to provide additional optimization of antithrombotic benefit and reduce bleeding risk. The study population was mainly male, which induced an obvious gender bias. Therefore, a larger prospective study should be needed to explore the separate models for predicting the risks of TEs and MB. This study did not find any correlation between CKD, which is considered as a dichotomous variable in this study, and incidence of bleeding. However, different scores developed in other studies (7, 23) examined not only the presence of renal impairment but also the grade of the dysfunction. Therefore, we might discriminate the grades of chronic kidney disease (stages 3, 4, and 5) in further larger studies. Finally, the cohort in this study was derived from patients who agreed to participate, which might have resulted in unaccounted selection pressures that affect the generalizability of the total cohort.

Conclusion

This study reported the development of separate models for predicting the risks of TEs and MB in subjects with multivessel coronary artery disease who had undergone primary PCI, which demonstrated moderate accuracy in discrimination concordant and stratified the risk in the derivation and validation cohorts. Furthermore, this study showed that subjects with a high incidence of thrombotic risk had greater bleeding risk in a large proportion of the study population.

Data Availability Statement

The datasets used and/or analyzed during this study are available from the corresponding author on reasonable request. Requests to access these datasets should be directed to hbyanfuwai2018@163.com.

Ethics Statement

The studies involving human participants were reviewed and approved by Ethics Committee of Fuwai Hospital. The patients/participants provided their written informed consent to participate in this study.

Author Contributions

HY, XZ, CL, PZ, ZS, JL, JZ, RC, YW, YC, LS, and HZ: substantial contributions to conception and design, data acquisition, or data analysis and interpretation, drafting the article or critically revising it for important intellectual content, final approval of the version to be published, and agreement to be accountable for all aspects of the study in ensuring that questions related to the accuracy or integrity of the study are appropriately investigated and resolved. All authors contributed to the article and approved the submitted version.

Funding

This study was supported by the Chinese Academy of Medical Sciences Innovation Fund for Medical Sciences (2016-I2M-1-009), National Natural Science Funds (Number: 81970308), the Fund of Sanming Project of Medicine in Shenzhen (Number: SZSM201911017) and Shenzhen Key Medical Discipline Construction Fund (number: SZXK001).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher's Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

The authors gratefully acknowledge all the individuals who participated in this study. We would like to thank Editage (www.editage.cn) for the English language editing.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcvm.2021.729432/full#supplementary-material

Footnote

1. ^The authors have completed the TRIPOD reporting checklist.

References

1. Capodanno D, Alfonso F, Levine NG, Valgimigli M, Angiolillo DJ. Acc/aha versus esc guidelines on dual antiplatelet therapy. J Am Coll Cardiol. (2018). 72:2915–31. doi: 10.1016/j.jacc.2018.09.057

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Généreux P, Giustino G, Witzenbichler B, Weisz G, Stuckey TD, Rinaldi MJ, et al. Incidence, predictors, and impact of post-dischargebleeding after percutaneous coronary intervention. J Am Coll Cardiol. (2015) 66:1036–45. doi: 10.1016/j.jacc.2015.06.1323

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Valgimigli M, Costa F, Lokhnygina Y, Clare MR, Wallentin L, Moliterno DJ, et al. Trade-off ofmyocardial infarction vs bleeding typeson mortality after acute coronary syndrome:lessons from the Thrombin Receptor Antagonist forClinical Event Reduction in Acute CoronarySyndrome (TRACER) randomized trial. Eur Heart J. (2017) 38:804–10. doi: 10.1093/eurheartj/ehw525

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Costa F, Klaveren DV, James S, Heg K, Räber L, Feres F, et al. Derivation and validation of the predictingbleeding complications in patients undergoingstent implantation and subsequentdual antiplatelet therapy (PRECISE-DAPT). score: a pooled analysis of individualpatientdatasets from clinical trials. Lancet. (2017) 389:1025–34. doi: 10.1016/S0140-6736(17)30397-5

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Robert WY, Eric AS, Dean JK, Sharon-Lise TN, Anthony HG, David JC, et al. Development and validation of a predictionrule for benefit and harm of dualantiplatelet therapy beyond 1 year afterpercutaneous coronaryintervention. JAMA. (2016) 315:1735–49. doi: 10.1001/jama.2016.3775

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Palmerini T, Genereux P, Caixeta A, Cristea E, Lansky A, Mehran R, et al. A newscore for risk stratification of patients with acutecoronary syndromes undergoing percutaneouscoronary intervention: the ACUITY-PCI (AcuteCatheterization and Urgent Intervention TriageStrategy-Percutaneous CoronaryIntervention) riskscore. JACC Cardiovasc Interv. (2012) 5:1108–16. doi: 10.1016/j.jcin.2012.07.011

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Subherwal S, Bach RG, Chen AY, Gage BF, Rao SV, Newby LK, et al. Baselinerisk of major bleeding innon-ST-segmentelevationmyocardial infarction: the CRUSADE (Can Rapid Risk Stratification of Unstable AnginaPatients Suppress Adverse Outcomes With EarlyImplementation of the ACC/AHA Guidelines) bleeding score. Circulation. (2009) 119:1873–82. doi: 10.1161/CIRCULATIONAHA.108.828541

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Peterson ED, Dai D, DeLong ER, Brennan JM, Singh M, Rao SV, et al. Contemporary mortality risk prediction forpercutaneouscoronary intervention: results from 588, 398 procedures in the nationalcardiovasculardata registry. J Am CollCardiol. (2010). 55:1923–32. doi: 10.1016/j.jacc.2010.02.005

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Boersma E, Pieper KS, Steyerberg EW, Wilcox RG, Chang WC, Lee KL, et al. Predictors ofoutcome in patients with acute coronary syndromeswithout persistent ST-segment elevation. Results from an international trial of 9461patients. Circulation. (2000) 101:2557–67. doi: 10.1161/01.CIR.101.22.2557

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Morrow DA, Antman EM, Charlesworth A, Cairns R, Murphy SA, de Lemos JA, et al. TIMI risk score for ST-elevation myocardialinfarction: a convenient, bedside, clinical score forrisk assessment at presentation: an intravenousnPA for treatment of infarcting myocardium earlyII trial substudy. Circulation. (2000) 102:2031–7. doi: 10.1161/01.CIR.102.17.2031

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Rao SV, McCoy LA, Spertus JA, Krone RJ, Singh M, Fitzgerald S, et al. Anupdated bleeding model to predict the risk ofpost-procedure bleeding among patients undergoingpercutaneous coronary intervention: areport using an expanded bleeding definition from the National Cardiovascular Data Registry CathPCIRegistry. JACC Cardiovasc Interv. (2013) 6:897–904. doi: 10.1016/j.jcin.2013.04.016

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Thygesen K, Alpert JS, Jaffe AS, Simoons ML, Chaitman BR, White HD. Third universal definition of myocardial infarction. J Am Coll Cardiol. (2012) 60:1581–98. doi: 10.1016/j.jacc.2012.08.001

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Mehran R, Rao SV, Bhatt DL, Gibson CM, Caixeta A, Eikelboom J, et al. Standardized bleeding definitions for cardiovascular clinical trials: a consensus report from the bleeding academic research consortium. Circulation. (2011). 123:2736–47. doi: 10.1161/CIRCULATIONAHA.110.009449

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Astor BC, Muntner P, Levin A, Eustace JA, Coresh J. Association of kidney function with anemia: the third national health and nutrition examination survey (1988–1994). Arch Intern Med. (2002) 162:1401–8. doi: 10.1001/archinte.162.12.1401

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Ines R, Gray AM, Crispin J, Murray DW, Oliver RA. Multiple imputation for patient reported outcome measures in randomised controlled trials: advantages and disadvantages of imputing at the item, subscale or composite score level. BMC Med Res Methodol. (2018). 18:87. doi: 10.1186/s12874-018-0542-6

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Sullivan LM, Massaro JM, D'Agostino RB Sr. Presentation of multivariate data for clinical use: the framingham study risk score functions. StatMed. (2004) 23:1631–60. doi: 10.1002/sim.1742

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Dangas GD, Claessen BE, Mehran R, Xu K, Fahy M, Parise H, et al. Development and validation of a stent thrombosisrisk score in patients with acute coronary syndromes. JACC Cardiovasc Interv. (2012) 5:1097–105. doi: 10.1016/j.jcin.2012.07.012

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Bohula EA, Aylward PE, Bonaca MP, Corbalan RL, Kiss RG, Murphy SA, et al. Efficacy and safety of vorapaxar with and without a thienopyridine for secondary prevention in patients with previous myocardial infarction and no history of stroke or transient ischemic attack: results from TRA 2°P-TIMI 50. Circulation. (2015) 132:1871–9. doi: 10.1161/CIRCULATIONAHA.114.015042

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Morrow DA. Cardiovascular risk prediction in patients with stableand unstable coronary heart disease. Circulation. (2010) 121:2681–91. doi: 10.1161/CIRCULATIONAHA.109.852749

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Généreux P, Madhavan MV, Mintz GS, Maehara A, Palmerini T, Lasalle L, et al. Ischemic outcomes after coronary intervention of calcified vessels in acute coronary syndromes. Pooled analysis from the HORIZONS-AMI (Harmonizing Outcomes With Revascularization and Stents in Acute Myocardial Infarction) and ACUITY (Acute Catheterization and Urgent Intervention Triage Strategy). TRIALS J Am CollCardiol. (2014) 63:1845–54. doi: 10.1016/j.jacc.2014.01.034

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Baber U, Mehran R, Giustino G, Cohen DJ, Henry TD, Sartori S, et al. Coronary thrombosis and major bleeding after PCI with drug-eluting stentsrisk scores from PARIS. J Am CollCardiol. (2016) 67:2224–34. doi: 10.1016/j.jacc.2016.02.064

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Bohula EA, Bonaca MP, Braunwald E, Aylward PE, Corbalan R, De Ferrari GM, et al. Atherothrombotic risk stratification and the efficacy and safety of vorapaxar inpatients with stable ischemic heart disease and previous myocardialinfarction. Circulation. (2016) 134:304–13. doi: 10.1161/CIRCULATIONAHA.115.019861

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Mehran R, Pocock SJ, Nikolsky E, Clayton T, Dangas GD, Kirtane AJ, et al. A risk score to predict bleeding in patients with acutecoronary syndromes. J Am CollCardiol. (2010) 55:2556–66. doi: 10.1016/j.jacc.2009.09.076

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Ducrocq G, Wallace JS, Baron G, Ravaud P, Alberts MK, Wilson PW, et al. Riskscore to predict serious bleeding in stable outpatientswith or at risk of atherothrombosis. EurHeart J. (2010) 31:1257–65. doi: 10.1093/eurheartj/ehq021

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Granger CB, Goldberg RJ, Dabbous O, Pieper KS, Eagle KA, Cannon CP, et al. Global registry of acute coronary events investigators. Predictors of hospital mortality in the global registry of acute coronary events. Arch Intern Med. (2003) 163:2345–53. doi: 10.1001/archinte.163.19.2345

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Amsterdam EA, Wenger NK, Brindis RG, Casey DE, Ganiats TG, Holmes DR Jr, et al. 2014 AHA/ACC guideline for the management of patients with non-ST-elevation acute coronary syndromes: executive summary: a report of the American College of Cardiology/ American Heart Association Task Force on Practice Guidelines. Circulation. (2014) 130:2354–94. doi: 10.1161/CIR.0000000000000133

PubMed Abstract | CrossRef Full Text | Google Scholar

27. O'Gara PT, Kushner FG, Ascheim DD, Casey DE, Chung MK, Lemos JA, et al. CF/AHA Task Force. 2013 ACCF/ AHA guideline for the management of ST-elevation myocardial infarction: executive summary: a report of the American College of Cardiology Foundation/American Heart Association Task Force on Practice Guidelines. Circulation. (2013) 127:529–55. doi: 10.1161/CIR.0b013e3182742c84

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Roffi M, Patrono C, Collet JP, Mueller C, Valgimigli M, Andreotti F, et al. 2015 ESC guidelines for the management of acutecoronary syndromes in patients presenting without persistent ST-segment elevation: task force for the management of acute coronary syndromes in patients presenting without persistent ST-segment elevation of the European Society of Cardiology (ESC). Eur Heart. (2016) J37:267–315. doi: 10.1093/eurheartj/ehv320

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Aoki J, Lansky AJ, Mehran R, Moses J, Bertrand ME, McLaurin BT, et al. Early stent thrombosis in patientswith acute coronary syndromes treated with drug-eluting and baremetal stents: the acute catheterization and urgent interventiontriage strategy trial. Circulation. (2009) 119:687–98. doi: 10.1161/CIRCULATIONAHA.108.804203

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Kukreja N, Onuma Y, Garcia-Garcia HM, Daemen J, Domburg RV, Serruys PW. The risk of stent thrombosis in patients with acute coronary syndromes treated with bare-metal and drug-eluting stents. JACC Cardiovasc Interv. (2009) 2:534–41. doi: 10.1016/j.jcin.2009.04.003

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Cronin L, Mehta SR, Zhao F, Pogue J, Budaj A, Hunt D, et al. Stroke in relation to cardiac proceduresin patients with non–ST-elevation acute coronary syndrome: astudy involving 18, 000 patients. Circulation. (2001) 104:269–74. doi: 10.1161/01.CIR.104.3.269

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Zimarino M, Angiolillo DJ, Dangas G, Capodanno D, Barbato E, Hahn JY, et al. Antithrombotic therapy after percutaneous coronary intervention of bifurcation lesions. EuroIntervention. (2021) 17:59–66. doi: 10.4244/EIJ-D-20-00885

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Costa F, Van Klaveren D, Feres F, James S, Räber L, Pilgrim T, et al. Dual antiplatelet therapy duration basedon ischemic and bleeding risks after coronarystenting. J Am Coll Cardiol. (2019) 73 741–54. doi: 10.1016/j.jacc.2018.11.048

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Mehta SR, Bainey KR, Cantor WJ, Lordkipanidzé M, Marquis-Gravel G, Robinson SD, et al. 2018 Canadian Cardiovascular Society (CCS) Canadian Association of Interventional Cardiology (CAIC) focused update of the guidelines for the use of antiplatelet therapy. Can J Cardiol. (2018) 34:214–33. doi: 10.1016/j.cjca.2017.12.012

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Giustino G, Chieffo A, Palmerini T, Valgimigli M, Feres F, Abizaid A, et al. Efficacy and safety of dual antiplatelet therapy after complex PCI. J Am Coll Cardiol. (2016) 68:1851–64. doi: 10.1016/j.jacc.2016.07.760

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Yeh RW, Kereiakes DJ, Steg PG, Cutlip DE, Croce KJ, Massaro JM, et al. Lesion complexity and outcomes of extended dual anti platelet therapy after percutaneous coronary intervention. J Am Coll Cardiol. (2017) 70:2213–23. doi: 10.1016/j.jacc.2017.09.011

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Antman EM, Cohen M, Bernink PJ, McCabe CH, Horacek T, Papuchis G, et al. TheTIMI risk score for unstable angina/non-ST elevationMI: a method for prognostication and therapeuticdecision making. JAMA. (2000) 284:835–42. doi: 10.1001/jama.284.7.835

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: thrombosis, major bleeding, primary PCI, multi-vessels coronary artery disease, inflammation

Citation: Zhao X, Liu C, Zhou P, Sheng Z, Li J, Zhou J, Chen R, Wang Y, Chen Y, Song L, Zhao H and Yan H (2022) Thrombosis and Major Bleeding Risk After Primary PCI Among Patients With Multivessel Coronary Artery Disease. Front. Cardiovasc. Med. 8:729432. doi: 10.3389/fcvm.2021.729432

Received: 22 September 2021; Accepted: 23 December 2021;
Published: 08 February 2022.

Edited by:

Per Morten Sandset, University of Oslo, Norway

Reviewed by:

Yasuhiro Nakano, Kyushu University, Japan
Marco Zimarino, Asl Lanciano Vasto Chieti, Italy

Copyright © 2022 Zhao, Liu, Zhou, Sheng, Li, Zhou, Chen, Wang, Chen, Song, Zhao and Yan. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Hongbing Yan, hbyanfuwai2018@163.com; Hanjun Zhao, 15210020808@163.com

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.