Background/Aims: Regorafenib is recommended by guidelines and trials as a sequential second-line therapy following progression on first-line sorafenib or lenvatinib in hepatocellular carcinoma (HCC). However, efficacy is limited, highlighting the urgent need to screen suitable patients and develop sensitization strategies. Methods: Acquired sorafenib- or lenvatinib-resistant (SR or LR) HCC cell lines and organoids were established. Genome-wide CRISPR library screen was performed in SR or LR cell strains to identify synthetic lethal targets of regorafenib. RNA-seq and FITC-regorafenib efflux assay were used to elucidate ERBB3-driven downstream signaling. Preclinical mouse models of cell line- and patient-derived xenografts and clinical cohorts of HCC patients were employed to validate the efficacy of ERBB3-guided patient stratification. Results: Screening with CRISPR library, we showed that inhibition of ERBB3 was synthetic lethal with regorafenib in SR or LR cell strains and organoids. Mechanistically, SR or LR triggered feedback activation of ERBB3 signaling and mediated regorafenib efflux via ERBB3-HIF1A-ABCB1 cascade pathway, limiting sensitivity to regorafenib. Moreover, ERBB3-low tumors following SR or LR exhibited significant sensitivity to regorafenib, suggesting its potential as a predictive biomarker to screen optimal candidates for sequential therapy. Seribantumab, an ERBB3-targeting monoclonal antibody, inhibited ERBB3-HIF1A-ABCB1 cascade, and its combination with regorafenib exerted marked synergistic anti-tumor effects on ERBB3-high tumors resistant to sorafenib or lenvatinib both in vitro and in vivo. Conclusions This study revealed that ERBB3 was a key resistance factor driving limited efficacy to sequential regorafenib, but also an effective therapeutic target whose inhibition enhanced regorafenib sensitivity after SR or LR.

BACKGROUND:The most common complication of traumatic femoral neck fractures after internal fixation is femoral head necrosis.Currently,many studies have reported on the risk factors that affect the occurrence and development of postoperative femoral head necrosis,but there is still a lack of tools to predict the risk of femoral head necrosis after internal fixation of femoral neck fractures.OBJECTIVE:To develop a predictive model that estimates the risk of femoral head necrosis shortly after patients with femoral neck fractures receive cannulated screw internal fixation.METHODS:A retrospective analysis reviewed clinical records of 172 patients who underwent cannulated screw internal fixation for femoral neck fractures at Department of Orthopedics of Mianyang Central Hospital from January 2013 to June 2023.Patients were categorized into two groups based on the presence or absence of femoral head necrosis within one year post-operation:the necrosis group and the non-necrosis group.Univariate analysis,Lasso regression,and multivariate Logistic regression techniques were employed to identify the determinants of femoral head necrosis.A nomogram prediction model was constructed using R language's"rms"package,version 4.0.The receiver operating characteristic curve was used to evaluate the discriminatory ability of the model.The Hosmer-Lemeshow test was used to evaluate the goodness of fit of the model,and the decision curve analysis was used to determine its clinical application benefits.Internal validation of the study was conducted using the Bootstrap method,involving 1 000 repeated samplings.To delve deeper into the primary factors influencing femoral head necrosis post-internal fixation of the femoral neck,this paper employed the SHAP method for data set analysis.RESULTS AND CONCLUSION:(1)The risk factors leading to femoral head necrosis in the short term after cannulated screw fixation of femoral neck fractures include:smoking,diabetes,Garden classification,fracture line location,reduction quality,age,and operation time.(2)The prediction model demonstrated robust performance,evidenced by an area under the curve of 0.940(95％Confidence Interval:0.903 to 0.977),indicating a high level of prediction accuracy.The model achieved a sensitivity of 90.2％and a specificity of 87.6％,indicating that its diagnostic performance was stable.The Hosmer-Lemeshow goodness-of-fit test yielded a chi-square value of 6.593 with a P-value of 0.581,confirming that the model's predictions closely align with the observed outcomes.(3)The calibration curve of the model also performed well,and its overall trend was very close to the ideal curve,further proving the high accuracy of the model.(4)The internal validation was carried out by the Bootstrap method with 1 000 repeated samplings,and the area under the curve of the model internal validation was still as high as 0.939,proving that the model had good stability.(5)Through the decision curve,it is found that within the probability threshold range of 1％to 92％,the model can obtain the maximum net benefit value.(6)The SHAP analysis results show that among the risk factors analyzed in this study,the location of the fracture line serves as the most significant predictor of femoral head necrosis following internal fixation with cannulated screws in femoral neck fractures,and subcapital fractures are extremely prone to femoral head necrosis after surgery.(7)It is concluded that the validated prediction model demonstrates strong discriminative power and reliability,offering practical clinical utility.It serves as a useful reference tool for short-term risk assessment of femoral head necrosis following internal fixation of femoral neck fractures.

BACKGROUND:The most common complication of traumatic femoral neck fractures after internal fixation is femoral head necrosis.Currently,many studies have reported on the risk factors that affect the occurrence and development of postoperative femoral head necrosis,but there is still a lack of tools to predict the risk of femoral head necrosis after internal fixation of femoral neck fractures.OBJECTIVE:To develop a predictive model that estimates the risk of femoral head necrosis shortly after patients with femoral neck fractures receive cannulated screw internal fixation.METHODS:A retrospective analysis reviewed clinical records of 172 patients who underwent cannulated screw internal fixation for femoral neck fractures at Department of Orthopedics of Mianyang Central Hospital from January 2013 to June 2023.Patients were categorized into two groups based on the presence or absence of femoral head necrosis within one year post-operation:the necrosis group and the non-necrosis group.Univariate analysis,Lasso regression,and multivariate Logistic regression techniques were employed to identify the determinants of femoral head necrosis.A nomogram prediction model was constructed using R language's"rms"package,version 4.0.The receiver operating characteristic curve was used to evaluate the discriminatory ability of the model.The Hosmer-Lemeshow test was used to evaluate the goodness of fit of the model,and the decision curve analysis was used to determine its clinical application benefits.Internal validation of the study was conducted using the Bootstrap method,involving 1 000 repeated samplings.To delve deeper into the primary factors influencing femoral head necrosis post-internal fixation of the femoral neck,this paper employed the SHAP method for data set analysis.RESULTS AND CONCLUSION:(1)The risk factors leading to femoral head necrosis in the short term after cannulated screw fixation of femoral neck fractures include:smoking,diabetes,Garden classification,fracture line location,reduction quality,age,and operation time.(2)The prediction model demonstrated robust performance,evidenced by an area under the curve of 0.940(95％Confidence Interval:0.903 to 0.977),indicating a high level of prediction accuracy.The model achieved a sensitivity of 90.2％and a specificity of 87.6％,indicating that its diagnostic performance was stable.The Hosmer-Lemeshow goodness-of-fit test yielded a chi-square value of 6.593 with a P-value of 0.581,confirming that the model's predictions closely align with the observed outcomes.(3)The calibration curve of the model also performed well,and its overall trend was very close to the ideal curve,further proving the high accuracy of the model.(4)The internal validation was carried out by the Bootstrap method with 1 000 repeated samplings,and the area under the curve of the model internal validation was still as high as 0.939,proving that the model had good stability.(5)Through the decision curve,it is found that within the probability threshold range of 1％to 92％,the model can obtain the maximum net benefit value.(6)The SHAP analysis results show that among the risk factors analyzed in this study,the location of the fracture line serves as the most significant predictor of femoral head necrosis following internal fixation with cannulated screws in femoral neck fractures,and subcapital fractures are extremely prone to femoral head necrosis after surgery.(7)It is concluded that the validated prediction model demonstrates strong discriminative power and reliability,offering practical clinical utility.It serves as a useful reference tool for short-term risk assessment of femoral head necrosis following internal fixation of femoral neck fractures.

Objective To understand mosquito density, species composition, and seasonal variations around the Beijing Capital International Airport(BCIA)and provide scientific evidence for mosquito-borne disease prevention and control. Methods Meteorological data were collected from 2016 to 2022, and mosquito density was monitored at seven surveillance sites around the BCIA from May to October each year using CO2-baited mosquito traps to analyze the relationship between mosquito density and meteorological factors. Results In total,68 518 female mosquitoes were captured, with a mosquito density of 10.20 per light · hour. The dominant species was Culex pipiens pallens(73.90%), followed by Aedes albopictus(9.71%)and Ae. vexans(9.29%);the mosquito density was highest in 2016(14.8 per light·hour)and lowest in 2022(4.46 per light·hour). Significant statistical differences were observed in mosquito density among different species(F=18.118, P<0.05). The density of Cx. pipiens pallens and other mosquito species showed a mutually exclusive trend; the peak of mosquito density varied based on years and habitats, with significant statistical differences in mosquito density among different habitats(F=8.504, P<0.05). The highest mosquito density was observed near the Wenyu River(21.50 per light·hour), whereas the lowest was recorded at the BCIA construction site(2.44 per light· hour). The monthly average temperature(r=0.595), monthly average highest temperature(r=0.575), and monthly average lowest temperature(r=0.624)showed moderate positive correlations, whereas the monthly average air pressure(r=-0.484)showed a moderate negative correlation. The monthly average minimum temperature was included in the regression model(F=25.575, P<0.000), and the equation Y=1.029X4-8.181 was used. The monthly average air pressure 2 months prior(b=-2.418, β=-1.619, P<0.05)and the monthly average relative humidity(b=-0.739, β=-1.201, P<0.05)significantly negatively predicted the mosquito density. Notably, the regression equation used was y=2526.170+(-2.418)X9+(-0.739)X13. The exposure-response analysis revealed that the density of mosquitoes was not linearly related to the average monthly air pressure, rainfall, and duration of sunshine. Conclusions The dominant mosquito species around the BCIA was Cx. pipiens pallens, Ae. albopictus, and Ae. vexans. During the control of Cx. pipiens pallens, the impact on other mosquito species must also be considered. Environmental management, breeding sites, and scientific use of pesticides should be prioritized based on the activity periods of mosquitoes. Further studies on mosquitoes and meteorological factors should be conducted to provide references and novel avenues for mosquito control.

Objective To investigate the capacity of municipal and county-level Centers for Disease Control (CDC) in Shandong Province for γ radionuclides in food as well as gross α and gross β in drinking water, and to provide evidence for strengthening radiological laboratories. Method Data regarding radioactive detection capacity were collected from 16 prefecture-level cities and 136 county-level CDCs. Chi-square test and Fisher's exact tests were used to compare capacity differences across regions and between cities with and without nuclear power plants. Correlation analysis was used to evaluate the spatial equity of gross α and gross β detection capacity for drinking water. Result The detection capacity of gross α and gross β in drinking water varied significantly among counties (χ²=81.026, P<0.001），no significant difference was found between coastal and inland cities (χ²=2.021，P>0.155); Cities with nuclear power plants showed better detection capacity (χ²=20.253，P<0.01). The detection capacity for γ radionuclides in food was primarily concentrated in the prefecture-level cities that host nuclear power plants. Conclusion The overall detection capacity for radioactivity in drinking water was acceptable but spatially uneven. The capacity for monitoring radioactive contamination in food was generally weak. It is necessary to strengthen capacity building at the grassroots level.

Peripheral nerve injury (PNI) severs peripheral connections and induces profound reorganization of primary somatosensory (S1) and motor (M1) cortical maps, a dynamic process that critically shapes the extent and quality of functional recovery. Brain-derived neurotrophic factor (BDNF) has emerged as a central molecular hub linking the initial peripheral insult to subsequent central plastic changes. This review systematically examines BDNF-mediated cortical remodeling following PNI, with a particular emphasis on its spatiotemporally specific regulatory mechanisms and the therapeutic opportunities they present. Immediately after PNI, deafferented cortical territories undergo rapid functional silencing and, over time, become progressively encroached upon by adjacent intact representations. This maladaptive reorganization—characterized by sensory map fusion and motor compensatory invasion—actively impedes successful reinnervation and optimal recovery. Such phenomena have been consistently documented in both non-human primate models and human fMRI studies, spanning sensory and motor modalities. Motor remapping is also prominently observed in facial nerve injury and chronic nerve compression models, underscoring the universality of these plastic changes. Notably, a range of interventions including regional local anesthesia, auditory-tactile substitution, and guided tactile imagery can effectively mitigate this aberrant remodeling, revealing promising avenues for the therapeutic modulation of cortical plasticity during critical post-injury windows. BDNF coordinates post-PNI cortical remodeling through three mechanistically interconnected pathways. First, BDNF drives dendritic and axonal growth via TrkB receptor signaling: retrogradely transported axonal signaling endosomes integrate CREB-dependent transcriptional programs with mTOR-mediated local protein synthesis to bolster dendritic arborization and axonal elongation, while BDNF-induced Limk1 translation precisely fine-tunes actin cytoskeletal dynamics, thereby stabilizing dendritic spine morphology. Second, BDNF governs GABAergic interneuron maturation and perineuronal net formation through the JNK signaling cascade, thereby acting as a “plasticity brake” that preserves the delicate excitatory/inhibitory balance and constrains indiscriminate, maladaptive rewiring. Third, BDNF robustly potentiates NMDAR-dependent long-term potentiation (LTP) by elevating NMDA channel open probability and increasing synaptic receptor density, thus tightly coupling structural remodeling with functional synaptic strengthening. Concurrently, long-term depression pathways are subject to modulation, with NMDAR antagonism exerting time-dependent, bidirectional effects that vary with the post-injury phase. Superimposed on these mechanisms, the functionally antagonistic actions of proBDNF/p75NTR (favoring pruning and growth cone collapse) and mature BDNF/TrkB (promoting survival and stabilization) add further layers of regulatory complexity, highlighting the absolute necessity for precise spatiotemporal control in any therapeutic strategy. Extending our previous experimental findings, we introduce a “cortical reclamation threshold” model. In this framework, early upregulation of BDNF-TrkB signaling enhances inhibitory tone to restrict cortical encroachment by neighboring intact regions, whereas delayed but precisely timed BDNF delivery following surgical nerve repair actively facilitates the reclamation of the original, deafferented cortical maps by regenerating afferent fibers. Importantly, the common BDNF Val66Met polymorphism may significantly modulate individual thresholds, providing a biologically grounded rationale for patient stratification and personalized timing of intervention. Despite this mechanistic promise, clinical translation faces considerable hurdles. These include inherently poor pharmacokinetics of BDNF, severely limited penetration across the blood-brain barrier, dose-dependent neurotoxicity mediated through the p75NTR receptor, and the overarching dual-edged nature of BDNF signaling, which can simultaneously promote both adaptive and maladaptive plasticity. Emerging delivery platforms——including engineered exosomes, polymeric nanoparticles, cryogel microcarriers, and adeno-associated virus (AAV)-based gene therapy——offer potential solutions, but all require rigorous pharmacokinetic, toxicological, and functional validation. Future research priorities should encompass: (1) comprehensive single-cell and spatial transcriptomic mapping to resolve BDNF cellular sources and real-time signaling dynamics across distinct cortical layers and cell types; (2) genotype-stratified dose-response studies to establish safe, effective, and personalized delivery protocols; and (3) large-scale multicenter clinical trials that seamlessly integrate BDNF-targeted interventions with multimodal rehabilitation strategies, such as repetitive transcranial magnetic stimulation and mirror therapy, with longitudinal neuroimaging serving as a biomarker for cortical reclamation. By systematically addressing these priorities, we aim to transform BDNF from a critical endogenous regulator into a precisely controllable and titratable therapeutic target, ultimately enabling a paradigm shift from mere structural repair to true functional neural reconstruction following PNI.

Objective:To investigate the clinical value of laparoscopic evaluation of anato-mical position of inferior mesenteric artery (IMA), inferior mesenteric vein (IMV) and left colic artery (LCA).Methods:The prospective one-arm study was conducted. The clinical data of 229 pati-ents who underwent laparoscopic left hemicolectomy for left colon or laparoscopic radical resection of rectal cancer in The Second Affiliated Hospital of Air Force Medical University from December 2022 to December 2023 were selected. The distance between the origin point of IMA and the origin point of the first branch (L1) as well as the distance from the origin point of LCA root to the junction of LCA and IMV (L2) were measured during the operation. IMA classification, the location relation-ship of LCA and IMV junction were recorded. Observation indicators: (1) situations of enrolled patients; (2) difference analysis between L1, L2 and clinical features; (3) distribution characteristics of the location relationship between LCA and IMV in different types of IMA. Mann-Whitney U test was used for comparison of measurement data with skewed distribution between groups, Kruskal-Wallis H test was used for comparison between multiple groups, and Dunn-Bonferroni test was used for pairwise comparison. Comparison of count data between groups was performed by chi-square test. Pearson or Spearman correlation analysis was conducted for correlation of continuous variables. Results:(1) Situations of enrolled patients. A total of 229 eligible patients were screened out, including 146 males and 83 females, aged 64(range, 55-71)years. The height of 229 patients was 168(range, 160-172)cm, the weight was 65.0(55.5,71.5)kg, the body surface area was (1.68±0.17)m 2, the tumor maximum diameter was 3.0(2.5,4.0)cm. The total number of lymph nodes dissected was 19(17,21), and the number of No.253 lymph node dissected was 4(3,5). The L1 was 3.50(1.20,8.00)cm, and the L2 was 2.20(0.50,7.30)cm. There were 58, 31, 32, 71, 22, 90, 26 and 212 patients with smoking, alcohol drinking, diabetes, hypertension, coronary heart disease, neoadjuvant chemo-therapy, neoadjuvant radiotherapy and preservation of the LCA, respectively. Among 229 patients, cases with BMI <18.5 kg/m 2, 18.5-23.9 kg/m 2 and >23.9 kg/m 2 were 11, 133 and 85, respectively. There were 153 cases in pathological stage Ⅰ-Ⅱ and 76 cases in stage Ⅲ. There were 168 cases of Dixon operation, 6 cases of Miles operation and 55 cases of sigmoid colon resection. There were 135 cases of IMA type 1, 44 cases of IMA type 2, 23 cases of IMA type 3, 2 cases of IMA type 4, and 25 cases of IMA type unable to judge. (2) Difference analysis between L1, L2 and clinical features.Correlation analysis showed negative correlation between the height, body surface area and L1 ( r=-0.17, -0.15, P<0.05). The L1 was 3.20(2.68,4.00)cm for male patients and 3.60(3.00,4.20)cm for female patients, respectively, showing a significant difference between the two groups ( Z=-2.37, P<0.05). The L1 of patients with IMA type 1, 2, and 3 was 3.20(2.80,4.00)cm, 3.85(3.00,4.48)cm, and 3.20(2.50,4.30)cm, respectively, showing a significant difference among them ( H=7.54, P<0.05). Further pairwise com-parison showed that there was a significant difference in L1 between patients with IMA type 2 and those with IMA type 1 ( P<0.05). The L2 of smokers and non-smokers were 2.50(1.95,3.20)cm and 2.20(1.60,2.80)cm, respectively, showing a significant difference between the two groups ( Z=-2.24, P<0.05). (3)Distribution characteristics of the location relationship between LCA and IMV in different types of IMA. There was no significant difference in LCA distribution between the anterior and posterior positions of IMV among the three IMA types (type 1, 2, 3) ( χ2=1.63, P>0.05). Conclusions:Patients with greater height have larger body surface area and shorter L1. L1 is significantly longer in female patients than in male patients. L1 is significantly longer in patients with IMA type 2 than in those with type 1. L2 is significantly longer in smokers than in non-smokers. There was no significant difference in the distribution location between LCA and IMV among patients of IMA type 1, 2 and 3.

Body weight abnormalities, including overweight, obesity, and underweight, have become a dual public health challenge in Chinese adults: overweight and obesity lead to a variety of chronic complications, while underweight increases the risks of malnutrition, sarcopenia, and organ dysfunction. To systematically address these issues, multidisciplinary experts in endocrinology, sports science, nutrition, and psychiatry from various regions have held multiple weight management seminars. Based on the latest epidemiological data and clinical evidence, they expanded the guideline to include assessment and intervention strategies for underweight, in addition to the core content of obesity management. This guideline outlines the etiological mechanisms, evaluation methods, and multidimensional management strategies for overweight and obesity, covering key areas such as diagnosis and assessment, medical nutrition therapy, exercise prescription, pharmacological intervention, and psychological support. It is intended to provide a scientific and standardized approach to weight management across the adult population, aiming to curb the rising prevalence of obesity, mitigate complications associated with abnormal body weight, and improve nutritional status and overall quality of life.

Objective To evaluate the efficacy and safety of transcatheter aortic valve implantation(TAVI)for the treatment of primary aortic valve regurgitation(NAVR)and to compare the difference in the choice of prosthetic valve size and the difference in complications with aortic stenosis(AS).Methods According to the definition of Valve Academic Research Consortium(VARC-3),143 patients with NAVR/AS treated with TAVI and patients with NAVR treated with surgical aortic valve replacement(SAVR)at Tongji Hospital,Tongji Medical College,Huazhong University of Science and Technology,China,from March 2019 to September 2024 were selected,and clinical data on baseline,perioperative,and primary endpoint events were were retrospectively collected and compared.Results Forty-three patients with NAVR were treated with TAVI,with a device success rate of 86.0％and a surgical success rate of 95.3％.Subgroup comparisons:(1)NAVR-TAVI group than NAVR-SAVR group:patients in the TAVI group had a significantly shorter operative time than those in the SAVR group(P＜0.001);complete left bundle branch block was more likely to occur after TAVI(P=0.042),and complete right bundle branch block was more likely to occur after SAVR(P=0.044).SAVR postoperatively The incidence of congestive heart failure was higher(P=0.013),and the mortality rate was significantly higher in the SAVR group than in the TAVI group(P=0.019).(2)NAVR-TAVI group than AS-TAVI group:the differences in access selection,THV size[28(22,34)mm vs.24(22,32)mm,P=0.044]and proportion of THV overdiameter[14％(7％,20％)vs.7％(3％,11％),P＜0.001]were statistically significant.patients in AS and NAVR groups had 1 case of permanent pacing after TAVI treatment.In the AS and NAVR groups,there was 1 case of permanent pacemaker implantation after TAVI.2 patients in the AS group were converted to surgical treatment,and 6 patients died.Conclusions The use of"off-label"(transfemoral)and"on-label"(transapical)TAVI devices(both from domestic sources)is safer than SAVR for the treatment of NAVR,especially in elderly and high-risk patients.Compared with patients with AS treated with TAVI,larger diameter annulas are usually selected for NAVR,with higher rates of valve migration,but overall safety and efficacy are comparable to AS.

Objectives:To explore the value of non-invasive myocardial work combined with myocardial contrast echocardiography(MCE)in the early diagnosis of coronary artery disease and its efficacy in stratifying the severity of coronary vessel lesions.Methods:A total of 130 patients with suspected coronary artery disease admitted to the First Affiliated Hospital of Dalian Medical University from June 2024 to January 2025 were enrolled in this study.All patients underwent echocardiography and MCE after admission,and coronary angiography(CAG).Based on CAG results,patients were divided into non-CAD group(n=45,coronary artery stenosis<50%)and CAD group(n=85,coronary artery stenosis≥50%).Patients in CAD group were further divided into low-score CAD group(≤49 points,n=43)and high-score CAD group(>49 points,n=42)according to the median of Gensini score(49 points).Non-invasive MW indices and quantitative MCE parameters were assessed.A binary logistic regression model was used to construct a combined diagnostic model,and the value of each parameter in diagnosing CAD and evaluating the severity of coronary lesions was calculated.The receiver operating characteristic(ROC)curve of subjects was estimated,and the area under the curve(AUC)was calculated to evaluate its sensitivity and specificity for early diagnosis of coronary heart disease.Results:Compared with the non-CAD group,the global longitudinal strain,global work index(GWI),and global constructive work(GCW)in both low-score and high-score CAD groups were significantly lower(all P<0.05),and the global work efficiency in the high-score CAD group was significantly reduced(P<0.05).MCE indices in both low-score and high-score CAD groups were significantly lower than those in the non-CAD group(all P<0.05).The multivariate logistic stepwise regression analysis and ROC curve showed that GWI(OR=0.997,95%CI:0.995-0.999,P=0.003)and A value(representing the peak intensity of the curve,reflecting myocardial blood volume(OR=0.415,95%CI:0.246-0.698,P=0.001)were independent predictors of low-score CAD.The combined diagnostic sensitivity and specificity for low-score coronary artery disease were 72.1%and 88.9%respectively,with an AUC of 0.851.GCW(OR=0.997,95%CI:0.995-1.000,P=0.019)and β-value(OR=0.000,95%CI:0.000-0.003,P<0.001)were independent predictors of high-score CAD.The combined diagnostic sensitivity and specificity for high-score coronary artery disease were 88.1%and 88.9%respectively,with an AUC of 0.934.Conclusions:Both non-invasive myocardial work parameters and MCE parameters have high diagnostic efficacy for coronary artery lesions of various degrees.The combined application of the two methods significantly improves the accuracy of coronary artery disease diagnosis,with improved sensitivity and specificity than single technique.Our results provide a new non-invasive comprehensive diagnostic model for clinical early diagnosis and risk stratification of coronary artery disease.