- 1Department of Cardiology, Liuzhou People’s Hospital, Affiliated of Guangxi Medical University, Liuzhou, Guangxi, China
- 2The Key Laboratory of Coronary Atherosclerotic Disease Prevention and Treatment of Liuzhou, Liuzhou, Guangxi, China
Purpose: To investigate the predictive value of leukocyte subsets and C-reactive protein (CRP) in coronary artery disease (CAD).
Methods: We conducted a Mendelian randomization analysis (MR) on leukocyte subsets, C-reactive protein (CRP) and CAD, incorporating data from 68,624 patients who underwent coronary angiography from 2010 to 2022. After initial screening, clinical data from 46,664 patients were analyzed. Techniques employed included propensity score matching (PSM), logistic regression, lasso regression, and random forest algorithms (RF). Risk factors were assessed, and the sensitivity and specificity of the models were evaluated using receiver operating characteristic (ROC) curves. Additionally, survival analysis was conducted based on a 36-month follow-up period.
Results: The inverse variance weight (IVW) analysis showed that basophil count (OR 0.92, 95% CI: 0.84–1.00, P = 0.048), CRP levels (OR 0.87, 95% CI: 0.73–1.00, P = 0.040), and lymphocyte count (OR 1.10, 95% CI: 1.04–1.16, P = 0.001) are significant risk factors for CAD. Using LASSO regression, logistic regression, and RF analysis, both CRP and lymphocyte counts were consistently identified as risk factors for CAD, prior to and following PSM. The ROC curve analysis indicated that the combination of lymphocyte and CRP levels after PSM achieves a higher diagnostic value (0.85). Survival analysis revealed that high lymphocyte counts and low CRP levels are associated with a decreased risk of Major Adverse Cardiovascular Events (MACE) (P < 0.001). Conversely, a higher CRP level combined with lymphocyte counts correlates with a poorer prognosis.
Conclusion: There is a causal relationship between lymphocytes, CRP and CAD. The combined assessment of CRP and lymphocytes offers diagnostic value for CAD. Furthermore, high CRP levels coupled with low lymphocyte counts are associated with a poor prognosis.
1 Introduction
As the global population ages, the incidence of coronary artery disease (CAD) continues to increase. According to the World Health Organization's Global Burden of Ischemic Heart Disease Report from 1990 to 2019, approximately 9.14 million people died from CAD in 2019, while an estimated 197 million people worldwide were affected by CAD (1). Identifying risk factors closely associated with CAD is crucial for its prevention and treatment, as well as for reducing the social burden of the disease.
The development of CAD begins with endothelial damage, leading to lipid accumulation, atherosclerotic plaque formation, and ultimately, the progression of these plaques, causing narrowing or blockage of the coronary arteries (2). Low-density lipoprotein (LDL) plays a pivotal role in the formation of lipid plaques, and reducing LDL levels has been widely accepted in clinical guidelines (3). However, even among populations that have achieved LDL control targets, a significant residual risk of CAD persists, primarily due to other CAD risk factors, including systemic inflammation (4–6). Various inflammatory cells and mediators contribute to the formation and progression of lipid plaques. Studies assessing the role of inflammation in CAD risk have yielded inconsistent results (7–10). Epidemiological and multiple prospective cohort studies have demonstrated associations between C-reactive protein (CRP), leukocytes and their subgroups, and coronary artery disease (CAD) risk (11–20). Some research has even suggested a protective effect on the heart by reducing CRP levels in rats (21). Nonetheless, Mendelian randomization studies have not confirmed a causal relationship between CRP levels and CAD risk (22, 23). These conflicting findings necessitate a thorough investigation of the relationship between CRP, leukocytes, their subgroups, and CAD risk.
Two-Sample Mendelian Randomization (TSMR) is a frequently used method for analyzing disease risk factors (24, 25). Unlike traditional observational studies, TSMR is largely unaffected by confounding factors (26, 27). However, TSMR studies often lack subsequent clinical validation, potentially leading to conclusions that may contradict those of traditional randomized controlled trials and extensive observational studies. Initially, this study employed the TSMR approach to explore the association between leukocyte subsets, CRP, and CAD. It then utilized propensity score matching (PSM) on the study population, followed by logistic regression analysis, Lasso regression, and random forest algorithms to validate the TSMR findings. A clinical model was developed to assess the diagnostic efficiency of the TSMR method. To enhance the accuracy of the survival analysis, the study population was further stratified using ROC curves, thereby improving the reliability of the study conclusions.
2 Method
2.1 Mendelian randomization analysis
2.1.1 Study design
The exposure factors were defined as leukocyte subsets and C-reactive protein (CRP) levels in the blood. Leukocyte subsets included eosinophil count, eosinophil percentage, monocyte count, neutrophil count, and lymphocyte count. The outcome variable was coronary artery disease (CAD).
2.1.2 Data collection
Genome-wide association study (GWAS) data for eosinophil count (bbj-a-12), eosinophil percentage (ebi-a-GCST004600), monocyte count (ebi-a-GCST90002339), neutrophil count (ebi-a-GCST9000235), lymphocyte count (ukb-d-30120-irnt), CRP level (ieu-b-4764), and CAD (ebi-a-GCST005194) were obtained from the open GWAS project website (gwas.mrcieu.ac.uk). The characteristics of the leukocyte subsets and coronary heart disease datasets are shown in Table 1.
2.1.3 Statistical analysis
In the two-sample Mendelian randomization (TSMR) analysis, instrumental variables (IVs) were employed for data analysis (26). The inclusion of IV genetic variants was required to satisfy the following three assumptions (27): (1) the IVs must be associated with leukocyte subset levels and CRP levels; (2) the IVs must not be influenced by confounding factors; and (3) the IVs must affect CAD exclusively through their association with leukocyte subset levels and CRP levels. To ensure a robust correlation between the IVs and leukocyte subset levels and CRP levels, a significance threshold of P < 5 × 10−8 was set, and single nucleotide polymorphisms (SNPs) meeting this criterion were selected as preliminary IVs. Additionally, the linkage disequilibrium coefficient r2r^2r2 was set to 0.001, and the region width was 10,000 kb to mitigate the influence of gene pleiotropy on the results. SNP information related to exposure was extracted from the CAD GWAS data. SNPs with high linkage were used to replace missing SNPs, and SNPs without replacement sites were excluded. The final IV dataset was obtained by combining data on leukocyte subset levels, CRP levels, and CAD. Four regression models were employed to assess the causal relationship between leukocyte subset levels, CRP levels, and CAD: MR-Egger regression, random effects inverse variance weighting (IVW), weighted median method, and weighted mode. Sensitivity analyses were conducted using heterogeneity tests, pleiotropy tests, and leave-one-out sensitivity tests. All analyses were performed using the TwoSampleMR package in R version 4.3.1, with a significance level (α) set at 0.05.
2.2 Case-control study
2.2.1 Study design
This part of the study aimed to leverage large clinical sample data to validate the robustness of Mendelian randomization results through multiple analytical methods and to explore the predictive value of lymphocytes and C-reactive protein (CRP) for major adverse cardiovascular events (MACE). MACE was defined as cardiovascular death (death due to cardiovascular causes such as myocardial infarction, heart failure, or arrhythmia), nonfatal myocardial infarction (a heart attack that did not result in death), and nonfatal stroke (a cerebrovascular event that did not result in death). Figure 1 provides a flowchart of the study design.
Figure 1. Process of recruiting patients for the study. CAD, coronary artery disease. CRP: C-reactive protein. Flowchart of patient recruitment and study design illustrating the selection and matching process to ensure balanced analysis groups.
2.2.2 Data collection
Baseline clinical data were collected from patient medical records using spreadsheets. This study included clinical data from 68,624 patients who underwent coronary angiography in the Department of Cardiology at Liuzhou People's Hospital between 2010 and 2022. Inclusion criteria were as follows: patients who underwent coronary angiography and were diagnosed with coronary artery disease (CAD) with >50% stenosis of the main coronary artery branches. Baseline characteristics, including demographic information, clinical characteristics, and laboratory measurements, were collected. Key variables included left ventricular ejection fraction (LVEF) as an indicator of cardiac function and medication use [including angiotensin-converting enzyme inhibitors (ACEI), β-blocker, and antiplatelet drugs]. These factors were included to better account for variables affecting prognosis. Exclusion criteria were as follows: (1) life expectancy <2 years, (2) missing key clinical data, and (3) conditions such as renal failure, bone marrow hematopoietic dysfunction, or malignant tumors. After screening, 5,467 patients with coronary heart disease and 41,197 patients without coronary heart disease were included in the study. Patients were followed up regularly at 3, 6, 9, and 12 months after enrollment through telephone follow-up, readmission follow-up, and outpatient follow-up. The total follow-up period was 36 months. Propensity score matching (PSM) was performed using a caliper value of 0.02. This study involved a retrospective data analysis, and patient identity information was anonymized during the experimental design phase, ensuring that patient privacy was not compromised.
2.2.3 Statistical analysis
For the clinical data of the study population, continuous data were expressed as mean ± standard deviation, and analysis of variance was used for intergroup comparisons. Categorical data were expressed as counts and ratios, with intergroup comparisons made using the chi-square test or Fisher's exact test. Propensity score matching (PSM) was employed to control potential confounders and create comparable patient groups for the CAD and non-CAD groups. LASSO regression was used to identify significant predictors of MACE, and the results were further validated using the random forest method. Multivariate logistic regression analysis was conducted to evaluate the predictive value of lymphocytes and CRP for MACE. For individuals meeting the diagnostic criteria for CAD, the receiver operating characteristic (ROC) curve was used to calculate the cutoff values of lymphocyte count and CRP level for grouping criteria, dividing the population into four groups: low lymphocyte count, high lymphocyte count, low CRP level, and high CRP level (28, 29). To evaluate the diagnostic performance of the models, we calculated the area under the ROC curve (AUC) for each model. The DeLong test was used to assess whether the differences in AUC between the models were statistically significant. Survival analysis was performed to assess the outcome of major adverse cardiovascular events (MACE).
3 Results
3.1 Mendelian randomization
The results are presented in Table 2. The MR-Egger method was employed to assess horizontal pleiotropy, with all groups demonstrating P-values greater than 0.05, indicating no evidence of horizontal pleiotropy. Heterogeneity among instrumental variables (IVs) was evaluated using Cochran's Q test (23). No heterogeneity was detected in monocyte count (P > 0.05). However, significant heterogeneity was observed among IVs for eosinophil count, eosinophil percentage, neutrophil count, CRP levels, and lymphocyte count, as indicated by Q statistics with P-values less than 0.05. Subsequent analyses were conducted using MR-Egger regression, random effects inverse variance weighting (IVW), the weighted median method, and the weighted mode, with particular emphasis on the IVW method. The IVW analysis revealed significant associations between coronary artery disease (CAD) and basophil count (IVW, OR 0.92, 95% CI: 0.84–1.00, P = 0.048), lymphocyte count (IVW, OR 1.10, 95% CI: 1.04–1.16, P = 0.001), and CRP levels (IVW, OR 0.87, 95% CI: 0.73–1.00, P = 0.040). These findings suggest causal relationships between CAD and eosinophil count, lymphocyte count, and CRP levels, as indicated by the IVW method (P < 0.05).
3.2 population characteristics
Propensity score matching (PSM) was employed to mitigate the influence of confounding variables and enhance the comparability between the experimental and control groups (30). A total of 46,664 participants were included in the study, comprising 5,467 patients with coronary artery disease (CAD) and 41,197 non-CAD individuals. Following PSM, each group consisted of 5,430 individuals. Before PSM, the CAD group exhibited significantly higher levels of age, systolic blood pressure (SBP), eosinophils, monocytes, homocysteine (Hcy), C-reactive protein (CRP), triglycerides (TG), and left ventricular ejection fraction (LVEF) compared to the control group (P < 0.05). Conversely, the CAD group had significantly lower levels of pulse pressure (PP), diastolic blood pressure (DBP), white blood cells (WBC), red blood cells (RBC), platelets (PLT), hematocrit (Hct), hemoglobin (Hb), CD4 + and CD8+ T cells, high-density lipoprotein cholesterol (HDL-C), low-density lipoprotein cholesterol (LDL-C), total cholesterol (TC), apolipoprotein B (ApoB), glycated hemoglobin (HbA1c), estimated glomerular filtration rate (eGFR), use of angiotensin-converting enzyme inhibitors (ACEI) or angiotensin II receptor blockers (ARBs), and beta-blockers (P < 0.001). After PSM, elevated WBC counts were observed only in the CAD group (P < 0.05), indicating a successful reduction of baseline differences between groups (see Table 3 for details).
3.3 Lasso regression
Lasso regression was employed to address multicollinearity and facilitate feature selection (31). In the analysis of pre-matched data, lasso regression identified the most significant variables at the optimal lambda value (−7.61), ranking them in descending order as follows: C-reactive protein (CRP), platelets (PLT), eosinophil count, percentage of eosinophils, and lymphocyte count. Similarly, the analysis of post-matched data at the optimal lambda value (−7.32) ranked the variables as follows: white blood cell count (WBC), CRP, lymphocyte count, eosinophil count, and percentage of eosinophils. Notably, CRP, eosinophil count, percentage of eosinophils, and lymphocyte count were consistently identified as significant variables regardless of matching status. This consistency underscores the robustness of these variables in predicting coronary artery disease (CAD). Figure 2 illustrates the results of the lasso regression analysis.
Figure 2. Lasso regression. (A,B) Demonstrate that before PSM, at the optimal lambda value (lambda = −7.61), 23 variables were retained, with variable importance ranking from highest to lowest as CRP, PLT, eosinophil count, percentage of eosinophils, and lymphocyte count, respectively. (C,D) Illustrate that after PSM, at the optimal lambda value (lambda = −7.32), 21 variables were retained, with variable importance ranking from highest to lowest as WBC, CRP, lymphocyte count, eosinophil count, and percentage of eosinophils, respectively. LASSO regression analysis identified CRP and lymphocyte count as significant predictors of CAD, demonstrating their importance in risk prediction.
3.4 Logistic regression
The dependent variable, coronary artery disease (CAD) diagnosis, was dichotomized using a binary categorical variable assignment method, with CAD coded as 1 and non-CAD as 0. A stepwise approach was subsequently applied to incorporate these variables into the logistic regression model. Before matching, the risk factors considered for CAD included age, gender, triglycerides (TG), monocyte count, glomerular filtration rate, total cholesterol (TC), body mass index (BMI), C-reactive protein (CRP), eosinophil percentage, platelets, hemoglobin (Hb), homocysteine (Hcy), high-density lipoprotein cholesterol (HDL-C), low-density lipoprotein cholesterol (LDL-C), pulse pressure (PP), apolipoprotein B (ApoB), and lymphocyte count. After matching, the risk factors for CAD were identified as Hcy, HDL-C, white blood cell count (WBC), age, Hb, BMI, glomerular filtration rate, CRP, lymphocyte count, and red blood cell count (RBC). CRP and lymphocyte count were consistently identified as significant risk factors regardless of matching status. Detailed results are provided in Tables 4, 5.
3.5 Random forest
The Random Forest (RF) algorithm determines the relative importance of each indicator by evaluating the average reduction in prediction accuracy (32). Figure 3 illustrates the importance rankings of risk factors for coronary artery disease (CAD), where higher values denote greater significance. RF analysis was performed on both pre-matched and post-matched data, with variables arranged in descending order of importance. Before matching, the top three variables were lymphocyte count, C-reactive protein (CRP), and monocyte count. After matching, these same top three variables—lymphocyte count, CRP, and monocyte count—remained the most significant. This consistency across matching statuses further emphasizes the critical role of these variables in predicting CAD.
Figure 3. Random forest. Before matching, (A) show that the top three ranked variables were lymphocyte count, CRP, and monocyte count, (B) show that the error rate of the model stabilizes when the number of classification trees exceeds 500; After matching, (C) show that the top three ranked variables remained lymphocyte count, CRP, and monocyte count, (D) show that the error rate of the model stabilizes when the number of classification trees exceeds 1,000. Random forest analysis highlights key predictors of CAD, emphasizing the importance of CRP and lymphocyte count in the model.
3.6 Diagnostic efficacy
Lymphocyte count and C-reactive protein (CRP) were identified as key variables associated with the onset and progression of coronary artery disease (CAD), warranting further evaluation of their diagnostic and prognostic capabilities. Initially, the cutoff value for lymphocyte count was set at 1.6 × 109/L, yielding a sensitivity of 46.5%, specificity of 65.6%, and an area under the curve (AUC) of 56.7% (P < 0.05). The cutoff for CRP was set at 0.2 mg/L, resulting in a sensitivity of 31.0%, specificity of 76.2%, and an AUC of 53.1% (P < 0.05). After propensity score matching (PSM), the optimal cutoff for lymphocyte count was adjusted to 0.6 × 109/L, showing a sensitivity of 83.2%, specificity of 18.6%, and an AUC of 49.5% (P < 0.05). The adjusted CRP cutoff was 0.4 mg/L, yielding a sensitivity of 20.5%, specificity of 83.1%, and an AUC of 51.9% (P < 0.05) (see Figure 4 for details). Model 1 (CRP + lymphocyte count) after matching achieved an AUC of 85% (P < 0.05), while Model 2 (CRP + lymphocyte count) before matching recorded an AUC of 63.8% (P < 0.05) (see Figure 5 for details). The analysis confirms that both lymphocyte count and CRP have substantial diagnostic significance for CAD, particularly when lymphocyte count exceeds 1.6 × 109/L and CRP exceeds 0.2 mg/L.
Figure 4. Clinical model. (A,B) Show that before matching the ROC curves for Lymphocyte (AUC = 0.567), and CRP (AUC = 0.531), (C) indicates that before matching, the diagnostic efficacy of lymphocyte count is higher than CRP (P < 0.001). (D,E) Show that after matching the ROC curves for Lymphocyte (AUC = 0.495), and CRP (AUC = 0.519), (F) indicates that after matching, the diagnostic efficacy of lymphocyte count is lower than CRP (P < 0.001). Comparison of ROC curves before and after PSM for lymphocyte count and CRP showed that the diagnostic performance was improved after matching.
Figure 5. Clinical model (combine). The ROC curves for Model 1 (Lymphocyte + CRP, after PSM); AUC = 0.85 (95% CI 0.819–0.881). The ROC curves for Model 2 (Lymphocyte + CRP, Before PSM); AUC = 0.638 (95% CI 0.6–0.679). PSM: Propensity Score Matching. The ROC curve of the combination of CRP and lymphocyte count showed a moderate diagnostic value for CAD, with improved accuracy after PSM.
3.7 Survival analysis
In univariate survival analyses, lymphocyte counts ≤1.6 × 109/L were identified as a risk factor for coronary artery disease (CAD) before propensity score matching (PSM). After PSM, lymphocyte counts ≤0.6 × 109/L continued to be a risk factor. Before PSM, C-reactive protein (CRP) levels >0.2 mg/L were consistently associated with an increased risk of CAD, and this association persisted with CRP levels >0.4 mg/L after PSM. Following a three-year follow-up period, higher lymphocyte counts and lower CRP levels were associated with a reduced risk of major adverse cardiovascular events (MACE) both before and after matching (P < 0.001) (Figures 6A,B,D,E). In the combined survival analysis, groups with low lymphocyte counts and low CRP levels served as the reference. Initially, the odds ratio (OR) for participants with high lymphocyte counts and low CRP levels was 2.033 (P < 0.001). The ORs were 2.772 for those with low lymphocyte counts and high CRP levels, and 1.53 for those with high levels of both lymphocytes and CRP (P < 0.001). Post-matching, the OR adjusted to 1.53 for the first group, 1.312 for those with high lymphocyte counts and low CRP levels, 1.47 for those with low lymphocyte counts and high CRP levels, and 1.303 for those with high levels of both indicators (P < 0.001) (Figures 6C,F). Notably, high CRP levels were associated with a more significant impact on the occurrence of MACE.
Figure 6. Survival analysis. (A) Before PSM matching, The follow-up of patients with different Lymphocyte group for 36 months. (B) Before PSM, The follow-up of patients with different CRP group for 36 months. (C) Before PSM matching, the follow-up of patients with different Lymphocyte + CRP group for 36 months (D) after PSM matching, The follow-up of patients with different Lymphocyte group for 36 months. (E) After PSM matching, The follow-up of patients with different CRP group for 36 months. (F) After PSM matching, The follow-up of patients with different Lymphocyte + CRP group for 36 months. Kaplan-Meier survival curves showed that patients with high CRP levels and low lymphocyte counts had a higher incidence of MACE, emphasizing the prognostic significance of these biomarkers.
4 Discussion
In this study, the causal relationship between leukocyte subtypes, CRP and CAD was first analyzed using TSMR. The study included 5,467 patients with confirmed CAD and 41,197 controls. Post PSM analysis, logistic regression, Lasso regression, and RF algorithms demonstrated strong associations between lymphocyte counts, CRP levels, and the development of CAD. ROC curves determined the optimal cutoff values for lymphocyte count (>1.6 × 109/L) and CRP levels (<0.2 mg/L) in the study population, while also evaluating their diagnostic sensitivity and specificity for CAD, both individually and combined. Survival analyses, considering various lymphocyte counts and CRP levels, indicated that lower lymphocyte counts and higher CRP levels were linked to an increased incidence of MACE, with a significant contribution from CRP levels.
CRP is an acute-phase reactant protein that reflects the body's inflammatory state (33). Recent evidence increasingly suggests that inflammation significantly contributes to the development of atherosclerosis and CAD (34–36). Beyond its role in forming atherosclerotic plaques, inflammation also contributes to their instability and rupture (5). CRP influences endothelial activation and dysfunction by modifying endothelial vasoreactivity, primarily through a reduction in endothelial nitric oxide synthase (eNOS) activity, which may heighten the risk of atherothrombosis, hypertension, and CAD (37). Elevated CRP levels (below 10 mg/L) are generally associated with an increased cardiovascular risk (5, 38, 39). Our findings align with this perspective, suggesting that high CRP levels may signal ongoing inflammatory activity within the arteries, potentially leading to plaque instability and an elevated risk of MACE. An eight-year follow-up study by Ridker et al. (40). involving 27,939 female patients showed a correlation between CRP levels and the incidence of cardiovascular events, indicating CRP's association with CAD onset in women. Cushman et al. (41). measured baseline CRP levels in 3,971 older adults without prior cardiovascular disease and conducted a ten-year follow-up. After adjusting for confounding factors, they found that individuals with CRP levels above 3 mg/L had a relative risk (RR) of 1.45 for CAD compared to those with levels above 1 mg/L. In a 1:1 case-control study of 420 CAD patients by Liu et al. (42), individuals in the acute coronary syndrome (ACS) group exhibited significantly higher CRP levels compared to controls. A notably higher incidence of MACE was observed among patients with elevated CRP levels. Even after adjusting for baseline confounders, CRP levels remained an independent predictor of MACE, with higher values linked to significantly greater all-cause mortality and myocardial infarction rates during the follow-up. In a five-year study of 3,802 participants, Zhuang Qian and colleagues (43) found that those with hs-CRP levels ≥1.08 mg/L faced a higher risk of developing CAD than those with hs-CRP levels <1.08 mg/L. However, subsequent MR analysis did not establish a significant causal relationship between hs-CRP and CAD. The difference in data sources from GWAS in our study might explain the divergent conclusions, attributable to variances in sample sources.
Lymphocytes play a crucial role in the human immune system, and alterations in their number and functionality are associated with various pathological conditions. Reduced lymphocyte levels can suggest chronic inflammation and immunosuppression, both involved in the pathophysiological processes of CAD (44). Atherosclerosis, a chronic inflammatory disease, increasingly recognizes the significance of lymphocytes. A decrease in lymphocyte count may reflect a reduced ability of the immune system to manage inflammatory responses, thus expediting the development and progression of coronary artery lesions. The pathological mechanisms involving lymphocytes in CAD are intricate, with different lymphocyte subtypes playing specific roles (44, 45). Zhang et al. (46). demonstrated that T lymphocytes facilitate atherosclerosis development by producing pro-inflammatory cytokines and adhesion molecules, prompting monocyte migration to the subendothelial layer. Additionally, numerous studies (47–49) have underscored the protective roles of natural and induced regulatory T cells (nTreg and iTreg) in atherosclerosis, possibly by diminishing the antigen-presenting activity of dendritic cells, thereby moderating both innate and adaptive immune cell activation. In mouse-based research, Kayw (50) and colleagues showed that B lymphocyte subtypes B1 and B2 jointly affect atherosclerosis development. B1 cells produce IgM antibodies that primarily recognize oxidized LDL (ox-LDL), conferring a protective effect against atherosclerosis, whereas B2 cells produce IgG and IgE, linked to the promotion of the disease. Acanfora et al. (51). demonstrated that lymphocyte levels are significantly correlated with the prognosis of patients with acute coronary syndrome, with low levels associated with an elevated risk of adverse cardiovascular events. Anoop Dinesh Shah et al. (52). conducted a cohort study revealing that low lymphocyte levels (≤ 1.45 × 109/L) were significantly associated with CAD mortality within the first six months (OR = 2.25, 95% CI 1.90–2.67), with a weaker correlation noted thereafter.
The established role of CRP as an inflammatory marker in CAD is well recognized, and the significance of lymphocytes in CAD is increasingly acknowledged. Despite this, the precise mechanisms and practical applications of this relationship warrant further investigation. The current study enhances the existing evidence by demonstrating a correlation between CRP levels and lymphocyte counts in CAD patients. Both peripheral blood lymphocyte counts and CRP levels may serve as indicators for assessing immune status and inflammation in these patients. Regular monitoring of CRP and lymphocyte levels in high-risk populations could enable early detection of individuals susceptible to CAD and facilitate timely preventative or therapeutic interventions. However, additional research is necessary to confirm the prognostic value of CRP and lymphocyte counts in CAD and to guide clinical decision-making. Future studies should explore the underlying mechanisms of CRP and lymphocyte fluctuations in CAD, evaluate their potential as therapeutic targets, and develop tailored intervention strategies to enhance the prevention and management of cardiovascular diseases.
Our study innovatively integrated MR with a comprehensive real-world analysis, which validated prior research conclusions and enhanced the reliability of these findings. However, it is critical to acknowledge the inherent limitations of this study. Being a single-center observational study, the generalizability of the results may be limited. Moreover, as a retrospective study, our research could be influenced by certain confounding factors and selection bias. The study also did not control for all possible confounders such as other chronic diseases, medication use, and lifestyle factors, which might affect CRP and lymphocyte levels, as well as the risk of CAD. In future research, we aim to expand the sample size and strive to include multi-center data. Additionally, more data and extended follow-up are necessary to predict risk factors and prognosis for CAD patients accurately. Therefore, it is essential to conduct further randomized controlled trials to validate our findings.
5 Conclusion
There is a causal relationship between lymphocytes, CRP and CAD. The combined assessment of CRP and lymphocytes offers diagnostic value for CAD. Furthermore, high CRP levels coupled with low lymphocyte counts are associated with a poor prognosis.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics statement
The studies involving humans were approved by Declaration of Helsinki (as revised in 2013). The Ethics Committee (No. LS-2010-008; March 11, 2010). The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent for participation was not required from the participants or the participants’ legal guardians/next of kin in accordance with the national legislation and institutional requirements.
Author contributions
YL: Conceptualization, Methodology, Writing – original draft. XY: Investigation, Writing – original draft. Y-CH: Investigation, Writing – original draft. Z-HB: Data curation, Writing – original draft. S-YL: Data curation, Writing – original draft. YL: Visualization, Writing – original draft. Y-LL: Visualization, Writing – review & editing. LM: Conceptualization, Methodology, Writing – original draft.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. The authors acknowledge the essential role of the funding of National Natural Science Foundation of China (NSFC: 82060072), National Natural Science Foundation of Guangxi (2020GXNSFAA297003), Project of Liuzhou Science and Technology (2022CAC0202), the project of Liuzhou people’s Hospital (LYRGCC202107 and LYRGCC202203), Guangxi Self-Financing Research Projects (Z20210459), Guangxi Medical and health key discipline construction project and The Key Laboratory of Coronary Atherosclerotic Disease Prevention and Treatment of Liuzhou.
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.
References
1. Safiri S, Karamzad N, Singh K, Carson-Chahhoud K, Adams C, Nejadghaderi SA, et al. Burden of ischemic heart disease and its attributable risk factors in 204 countries and territories, 1990–2019. Eur J Prev Cardiol. (2022) 29(2):420–31. doi: 10.1093/eurjpc/zwab213
2. Hansson GK. Inflammation, atherosclerosis, and coronary artery disease. N Engl J Med. (2005) 352(16):1685–95. doi: 10.1056/NEJMra043430
3. Virani SS, Newby LK, Arnold SV, Bittner V, Brewer LC, Demeter SH, et al. 2023 AHA/ACC/ACCP/ASPC/NLA/PCNA Guideline for the management of patients with chronic coronary disease: a report of the American Heart Association/American College of Cardiology joint committee on clinical practice guidelines. Circulation. (2023) 148(9):e9–e119. doi: 10.1161/CIR.0000000000001168. Erratum in: Circulation. (2023) 148(13):e148. doi: 10.1161/CIR.0000000000001183. Erratum in: Circulation. 148(23):e186. doi: 10.1161/CIR.0000000000001195
4. Shaya GE, Leucker TM, Jones SR, Martin SS, Toth PP. Coronary heart disease risk: low-density lipoprotein and beyond. Trends Cardiovasc Med. (2022) 32(4):181–94. doi: 10.1016/j.tcm.2021.04.002
5. Attiq A, Afzal S, Ahmad W, Kandeel M. Hegemony of inflammation in atherosclerosis and coronary artery disease. Eur J Pharmacol. (2024) 966:176338. doi: 10.1016/j.ejphar.2024.176338
6. Libby P, Buring JE, Badimon L, Hansson GK, Deanfield J, Bittencourt MS, et al. Atherosclerosis. Nat Rev Dis Primers. (2019) 5(1):56. doi: 10.1038/s41572-019-0106-z
7. Ridker PM, Everett BM, Thuren T, MacFadyen JG, Chang WH, Ballantyne C, et al. Antiinflammatory therapy with canakinumab for atherosclerotic disease. N Engl J Med. (2017) 377(12):1119–31. doi: 10.1056/NEJMoa1707914
8. Aimo A, Pascual-Figal DA, Barison A, Cediel G, Vicente ÁH, Saccaro LF, et al. Colchicine for the treatment of coronary artery disease. Trends Cardiovasc Med. (2021) 31(8):497–504. doi: 10.1016/j.tcm.2020.10.007
9. Tardif JC, Kouz S, Waters DD, Bertrand OF, Diaz R, Maggioni AP, et al. Efficacy and safety of low-dose colchicine after myocardial infarction. N Engl J Med. (2019) 381(26):2497–505. doi: 10.1056/NEJMoa1912388
10. Ridker PM, Everett BM, Pradhan A, MacFadyen JG, Solomon DH, Zaharris E, et al. Low-dose methotrexate for the prevention of atherosclerotic events. N Engl J Med. (2019) 380(8):752–62. doi: 10.1056/NEJMoa1809798
11. Wheeler J. Associations between differential leucocyte count and incident coronary heart disease: 1764 incident cases from seven prospective studies of 30 374 individuals. Eur Heart J. (2004) 25(15):1287–92. doi: 10.1016/j.ehj.2004.05.002
12. Strang F, Schunkert H. C-reactive protein and coronary heart disease: all said–is not it? Mediators Inflamm. (2014) 2014:757123. doi: 10.1155/2014/757123
13. Emerging Risk Factors Collaboration, Kaptoge S, Di Angelantonio E, Lowe G, Pepys MB, Thompson SG, Collins R, et al. C-reactive protein concentration and risk of coronary heart disease, stroke, and mortality: an individual participant meta-analysis. Lancet. (2010) 375(9709):132–40. doi: 10.1016/S0140-6736(09)61717-7
14. Shrivastava AK, Singh HV, Raizada A, Singh SK. C-reactive protein, inflammation and coronary heart disease. The Egyptian Heart Journal. (2015) 67(2):89–97. doi: 10.1016/j.ehj.2014.11.005
15. Prins BP, Abbasi A, Wong A, Vaez A, Nolte I, Franceschini N, et al. Investigating the causal relationship of C-reactive protein with 32 complex somatic and psychiatric outcomes: a large-scale cross-consortium Mendelian randomization study. PLoS Med. (2016) 13(6):e1001976. doi: 10.1371/journal.pmed.1001976
16. Wang J, Song J, Wu J, He C, Xu C, Liu Y. Leukocyte and leukocyte subset counts reveal compensatory mechanisms in coronary heart disease. Clin Chim Acta. (2013) 418:79–85. doi: 10.1016/j.cca.2012.12.028
17. Goeman JJ. L1 penalized estimation in the cox proportional hazards model. Biom J. (2010) 52(1):70–84. doi: 10.1002/bimj.200900028
18. Galimzhanov A, Sabitov Y, Guclu E, Tenekecioglu E, Mamas MA. Phenotyping for percutaneous coronary intervention and long-term recurrent weighted outcomes. Int J Cardiol. (2023) 374:12–9. doi: 10.1016/j.ijcard.2022.12.035
19. Özmen M. The role of C-reactive protein:albumin ratio and neutrophil:lymphocyte ratio in predicting coronary artery disease. Cardiovasc J Afr. (2024) 34:1–4. doi: 10.5830/cvja-2024-004
20. Kelesoglu S, Yilmaz Y, Elcık D. Relationship between C-reactive protein to albumin ratio and coronary collateral circulation in patients with stable coronary artery disease. Angiology. (2021) 72(9):829–35. doi: 10.1177/00033197211004392
21. Pepys MB, Hirschfield GM, Tennent GA, Gallimore JR, Kahan MC, Bellotti V, et al. Targeting C-reactive protein for the treatment of cardiovascular disease. Nature. (2006) 440(7088):1217–21. doi: 10.1038/nature04672
22. Kuppa A, Tripathi H, Al-Darraji A, Tarhuni WM, Abdel-Latif A. C-reactive protein levels and risk of cardiovascular diseases: a two-sample bidirectional Mendelian randomization study. Int J Mol Sci. (2023) 24(11):9129. doi: 10.3390/ijms24119129
23. Zhao Q, Liu R, Chen H, Yang X, Dong J, Bai M. White blood cells and coronary heart disease: a Mendelian randomization study. Front Genet. (2023) 14:1127820. doi: 10.3389/fgene.2023.1127820
24. Sanderson E, Glymour MM, Holmes MV, Kang H, Morrison J, Munafò MR, et al. Mendelian randomization. Nat Rev Methods Primers. (2022) 2:6. doi: 10.1038/s43586-021-00092-5
25. Miao L, Deng GX, Yin RX, Nie RJ, Yang S, Wang Y, et al. No causal effects of plasma homocysteine levels on the risk of coronary heart disease or acute myocardial infarction: a Mendelian randomization study. Eur J Prev Cardiol. (2021) 28(2):227–34. doi: 10.1177/2047487319894679
26. Larsson SC, Butterworth AS, Burgess S. Mendelian randomization for cardiovascular diseases: principles and applications. Eur Heart J. (2023) 44(47):4913–24. doi: 10.1093/eurheartj/ehad736
27. Jansen H, Samani NJ, Schunkert H. Mendelian randomization studies in coronary artery disease. Eur Heart J. (2014) 35(29):1917–24. doi: 10.1093/eurheartj/ehu208
28. Tiwari A, Chugh A, Sharma A. Ensemble framework for cardiovascular disease prediction. Comput Biol Med. (2022) 146:105624. doi: 10.1016/j.compbiomed.2022.105624
29. Wang Z, Sun Z, Yu L, Wang Z, Li L, Lu X. Machine learning-based prediction of composite risk of cardiovascular events in patients with stable angina pectoris combined with coronary heart disease: development and validation of a clinical prediction model for Chinese patients. Front Pharmacol. (2023) 14:1334439. doi: 10.3389/fphar.2023.1334439
30. Wang Z, Wu J, Zhang D. Hysterectomy and ischemic heart disease: an observational study using propensity score methods in NHANES 2007-2018. Atherosclerosis. (2021) 327:5–12. doi: 10.1016/j.atherosclerosis.2021.04.009
31. Ninomiya K, Kageyama S, Shiomi H, Kotoku N, Masuda S, Revaiah PC, et al. Can machine learning aid the selection of percutaneous vs surgical revascularization? J Am Coll Cardiol. (2023) 82(22):2113–24. doi: 10.1016/j.jacc.2023.09.818
32. Hu J, Szymczak S. A review on longitudinal data analysis with random forest. Brief Bioinform. (2023) 24(2):bbad002. doi: 10.1093/bib/bbad002
33. Calabró P, Willerson JT, Yeh ETH. Inflammatory cytokines stimulated C-reactive protein production by human coronary artery smooth muscle cells. Circulation. (2003) 108(16):1930–2. doi: 10.1161/01.CIR.0000096055.62724.C5
34. Libby P, Ridker PM. Inflammation and atherosclerosis: role of C-reactive protein in risk assessment. Am J Med. (2004) 116(Suppl 6A):9S–16. doi: 10.1016/j.amjmed.2004.02.006
35. Badimon L, Storey RF, Vilahur G. Update on lipids, inflammation and atherothrombosis. Thromb Haemost. (2011) 105(Suppl 1):S34–42. doi: 10.1160/ths10-11-0717
36. Mohammadnia N, Opstal TSJ, El Messaoudi S, Bax WA, Cornel JH. An update on inflammation in atherosclerosis: how to effectively treat residual risk. Clin Ther. (2023) 45(11):1055–9. doi: 10.1016/j.clinthera.2023.08.016
37. Devaraj S, Kumaresan PR, Jialal I. C-reactive protein induces release of both endothelial microparticles and circulating endothelial cells in vitro and in vivo: further evidence of endothelial dysfunction. Clin Chem. (2011) 57(12):1757–61. doi: 10.1373/clinchem.2011.169839
38. Kushner I, Antonelli MJ. What should we regard as an “elevated” C-reactive protein level? Ann Intern Med. (2015) 163(4):326. doi: 10.7326/L15-5126
39. Ridker PM. Clinical application of C-reactive protein for cardiovascular disease detection and prevention. Circulation. (2003) 107(3):363–9. doi: 10.1161/01.CIR.0000053730.47739.3C
40. Ridker PM, Rifai N, Rose L, Buring JE, Cook NR. Comparison of C-reactive protein and low-density lipoprotein cholesterol levels in the prediction of first cardiovascular events. N Engl J Med. (2002) 347(20):1557–65. doi: 10.1056/NEJMoa021993
41. Cushman M, Arnold AM, Psaty BM, Manolio TA, Kuller LH, Burke GL, et al. C-reactive protein and the 10-year incidence of coronary heart disease in older men and women: the cardiovascular health study. Circulation. (2005) 112(1):25–31. doi: 10.1161/CIRCULATIONAHA.104.504159
42. Liu R, Xu F, Ma Q, Zhou Y, Liu T. C-reactive protein level predicts cardiovascular risk in Chinese young female population. Oxid Med Cell Longev. (2021) 2021:6538079. doi: 10.1155/2021/6538079
43. Zhuang Q, Shen C, Chen Y, Zhao X, Wei P, Sun J, et al. Association of high sensitive C-reactive protein with coronary heart disease: a Mendelian randomization study. BMC Med Genet. (2019) 20(1):170. doi: 10.1186/s12881-019-0910-z
44. Ammirati E, Moroni F, Magnoni M, Camici PG. The role of T and B cells in human atherosclerosis and atherothrombosis. Clin Exp Immunol. (2015) 179(2):173–87. doi: 10.1111/cei.12477
45. Saigusa R, Winkels H, Ley K. T cell subsets and functions in atherosclerosis. Nat Rev Cardiol. (2020) 17(7):387–401. doi: 10.1038/s41569-020-0352-5
46. Zhang Y, Liu H, Tang W, Qiu Q, Peng J. Resveratrol prevents TNF-α-induced VCAM-1 and ICAM-1 upregulation in endothelial progenitor cells via reduction of NF-κB activation. J Int Med Res. (2020) 48(9):300060520945131. doi: 10.1177/0300060520945131
47. Kushwah R, Hu J. Role of dendritic cells in the induction of regulatory T cells. Cell Biosci. (2011) 1(1):20. doi: 10.1186/2045-3701-1-20
48. Lan Q, Zhou X, Fan H, Chen M, Wang J, Ryffel B, et al. Polyclonal CD4+Foxp3+ treg cells induce TGFβ-dependent tolerogenic dendritic cells that suppress the murine lupus-like syndrome. J Mol Cell Biol. (2012) 4(6):409–19. doi: 10.1093/jmcb/mjs040
49. Yang J, Yang Y, Fan H, Zou H. Tolerogenic splenic IDO (+) dendritic cells from the mice treated with induced-treg cells suppress collagen-induced arthritis. J Immunol Res. (2014) 2014:831054. doi: 10.1155/2014/831054
50. Kyaw T, Tipping P, Bobik A, Toh BH. Opposing roles of B lymphocyte subsets in atherosclerosis. Autoimmunity. (2017) 50(1):52–6. doi: 10.1080/08916934.2017.1280669
51. Acanfora D, Gheorghiade M, Trojano L, Furgi G, Pasini E, Picone C, et al. Relative lymphocyte count: a prognostic indicator of mortality in elderly patients with congestive heart failure. Am Heart J. (2001) 142(1):167–73. doi: 10.1067/mhj.2001.115792
Keywords: coronary artery disease, lymphocytes, CRP, Mendelian randomization, retrospective analysis
Citation: Liu Y, Yuan X, He Y-C, Bi Z-H, Li S-Y, Li Y, Liu Y-L and Miao L (2024) Exploring the predictive values of CRP and lymphocytes in coronary artery disease based on a machine learning and Mendelian randomization. Front. Cardiovasc. Med. 11:1442275. doi: 10.3389/fcvm.2024.1442275
Received: 1 June 2024; Accepted: 26 August 2024;
Published: 9 September 2024.
Edited by:
Hiroki Teragawa, JR Hiroshima Hospital, JapanReviewed by:
Wenming He, Ningbo First Hospital, ChinaDi Wang, Shanghai Jiao Tong University School of Medicine, China
Hongwei Li, Guangzhou First People’s Hospital, China
Copyright: © 2024 Liu, Yuan, He, Bi, Li, Li, Liu and Miao. 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: Yan-Li Liu, Z3hseWxAMTI2LmNvbQ==; Liu Miao, ZHIubWlhb2xpdUBxcS5jb20=
†These authors have contributed equally to this work