Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Public Health, 15 September 2025

Sec. Infectious Diseases: Epidemiology and Prevention

Volume 13 - 2025 | https://doi.org/10.3389/fpubh.2025.1641156

Association between COVID-19 vaccination and coronary heart disease: based on 2023 national health interview survey data

Juncong MaoJuncong Mao1Yunfei Hong
Yunfei Hong1*Rui Bao
Rui Bao2*
  • 1Department of Cardiology, The Third People’s Hospital of Yunnan, Kunming, Yunnan, China
  • 2Department of SICU, Fuwai Yunnan Hospital, Chinese Academy of Medical Sciences, Affiliated Cardiovascular Hospital of Kunming Medical University, Kunming, China

Background: Coronary heart disease (CHD) represents a critical cardiovascular ailment necessitating thorough investigation. This research endeavors to explore the potential link between COVID-19 vaccination and CHD, using data from the National Health Interview Survey (NHIS).

Methods: The study encompasses 20,906 participants from the 2023 NHIS cohort, and these participants were stratified into two groups: CHD patients and non-CHD individuals (controls). To determine the protective factors for CHD, both univariate and multivariate logistic regression analyses were carried out. Furthermore, Receiver Operating Characteristic (ROC) curves were plotted to assess the predictive performance of models that consider COVID-19 vaccination as a potential protective factor against CHD.

Results: In this study, a number of potential risk factors were investigated, including age (AGEP_A), sex (SEX_A), and race (RACEALLP_A) et al. Among them, the number of COVID-19 vaccinations was confirmed to be an effective protective factor for preventing coronary heart disease. Univariate logistic regression analysis showed that the risk of coronary heart disease was reduced in people who received 2 doses (OR = 0.68, 95% CI 0.49–0.92, p = 0.016), 3 doses (OR = 0.47, 95% CI 0.34–0.63, p < 0.001), 4 doses (OR = 0.39, 95% CI 0.28–0.53, p < 0.001), 5 doses (OR = 0.31, 95% CI 0.22–0.43, p < 0.001), and 6 or more doses (OR = 0.21, 95% CI 0.14–0.32, p < 0.001) of the COVID-19 vaccine. In the multivariate logistic regression analysis, Model 3 (after adjusting for multiple covariates) showed that the OR for those who received 6 or more doses of the vaccine was 0.459 (95% CI 0.289–0.726, p < 0.001). Receiver Operating Characteristic (ROC) curve analysis showed that the AUC for predicting coronary heart disease based on COVID-19 vaccination status was 0.845 (95% CI 0.8357–0.8539).

Conclusion: Based on NHIS database, a predictive model for CHD has been developed, and COVID-19 vaccinations were identified as a protective factor against CHD. This model holds potential clinical value.

1 Introduction

As the leading global cause of mortality, the disease burden of Coronary heart disease (CHD) continues to escalate. World Health Organization (WHO) data indicate that approximately 9 million deaths annually are attributed to CHD, accounting for 16% of global mortality (1). The economic burden is substantial, with annual healthcare costs exceeding $200 billion in the United States, while developing countries face treatment disparities due to resource limitations (2, 3). Despite advances in treatment, the five-year incidence of major adverse cardiovascular events remains high (20%) (4).

In China, the rising incidence among older adults populations is particularly pronounced, projecting a severe public health crisis by 2030 (2). CHD exhibits distinct sex-specific disparities: males face higher lifetime risks due to physiological hormonal differences and behavioral factors (e.g., smoking, alcohol use), while postmenopausal females experience a sharp decline in estrogen-mediated cardioprotection, leading to significantly increased morbidity (5, 6). These findings have highlighted the need to develop novel prevention strategies, including an evaluation of the role of COVID-19 vaccination.

The pandemic has profoundly disrupted CHD care delivery. Dual challenges emerged: (1) delayed elective procedures (e.g., a 38% reduction in global coronary intervention volume) and suspended cardiac rehabilitation programs escalated disease progression risks (7); (2) although remote monitoring technologies (e.g., wearable ECG devices) enabled home-based care, the “digital divide” (device utilization <30% among older adults/low-income groups) limited equitable access (8). Mechanistically, SARS-CoV-2 invades cardiomyocytes via ACE2 receptors, triggering endothelial inflammatory storms (5-8-fold increases in IL-6/TNF-α levels) and plaque activation, thereby augmenting 1-year acute myocardial infarction risk by 63% in infected individuals (9). Prior studies have demonstrated associations between COVID-19 vaccination and various cardiovascular conditions (10, 11), with additional research examining vaccination rates and predictive factors in adults with CHD (4), no conclusive evidence has established a direct causal relationship between COVID-19 vaccination and CHD incidence or progression. This knowledge gap underscores the critical need for rigorous epidemiological investigations to clarify potential biological interactions.

Despite advances facilitated by NHIS and other large-scale databases, critical gaps persist: (1) mechanisms underlying racial/sex heterogeneity in vaccine efficacy remain elusive; (2) dose-dependent obesity-mediated attenuation of vaccine effectiveness lacks therapeutic targets; (3) longitudinal validation of temporal associations between long COVID and CHD progression is warranted. This study aims to investigate the association between COVID-19 vaccination and the incidence of coronary heart disease (CHD) based on data from the 2023 National Health Interview Survey (NHIS), to develop a CHD prediction model incorporating vaccination status and other covariates, and to provide a novel theoretical basis for the early diagnosis and prevention of CHD.

2 Materials and methods

2.1 Source of data

As the largest health monitoring system in the United States, the NHIS used a multi-stage stratified probability sampling method to track population-level health behaviors and chronic disease patterns, providing nationally representative data (6).

The study contained adults with a history of coronary heart disease (CHD; n = 29,522) from the 2023 NHIS.1 The criteria for exclusion were as outlined below: (1) participants under the age of 18 years; (2) participants with uncertain CHD diagnoses or incomplete information pertaining to CHD; (3) those with missing data for any variables of interest. Following the application of these criteria, the number of participants, which summed up to 20,906, was selected for inclusion in the analysis (Supplementary Table 1). Subsequently, sample size estimation was performed using the R package “pwr” (v 1.3–0) (12). The expected effect size was set at 0.1, the significance level (α) at 0.05, and the statistical power (1-β) at 0.8. The results revealed that the minimum required sample size was 1,570, while the actual sample size included in this study was much larger than this value, ensuring the statistical power of the research. To determine the presence of notable disparities in 16 baseline characteristics (or covariates) across the two groups, the enrolled participants were divided into a CHD group and a control group for comparison. The statistical software package “tableone,” version 0.13.2, from the R programming environment (13) was utilized.

2.2 Dependent variable (outcome)

The definition of CHD: Data from the 2023 NHIS were utilized to classify patients with and without CHD into the CHD and control groups. Specifically, the “Adult Sample” file was selected from the “Data Release” section, and the relevant variable was identified by accessing the “Variables” layout. CHD status (self-reported) was determined based on responses to the variable ID: CHDEV_A, which asks, “Have you ever been told by a doctor or other health professional that you have coronary heart disease?” Participants giving an affirmative response were assigned to the CHD group, whereas those providing a negative reply were allocated to the control group.

2.3 Definition of the exposure factor

The determination of the exposure factor (number of COVID-19 vaccinations) was based on responses to the variable ID: SHTCVD19NM1_A, which inquired about the “Number of COVID-19 vaccinations.” The answers were as follows: those who received one vaccination, those who received two vaccinations, those who received three vaccinations, those who received four vaccinations, those who received five vaccinations, and those who received six or more vaccinations.

2.4 Covariates

The covariates were classified into 4 categories: sociodemographic, health status, related diseases, and healthcare data. Specifically, these include age, sex, race, smoking, number of COVID-19 vaccinations, marital status, health insurance, household income poverty rate, educational level, diabetes, hypertension, high cholesterol, food security, depression, pulmonary disease and body mass index (BMI).

2.5 Model establishment

Firstly, the R package “tableone” (v 0.13.2) (10) was employed to construct the baseline table. Next, to determine whether each variable served as a risk factor for CHD, univariate logistic regression analyses were carried out to elucidate the relationship between every single variable and the incidence of CHD. Additionally, three logistic regression models were developed to assess the link between exposure factors and CHD, with the computation of odds ratio values and their corresponding 95% confidence intervals. Model 1 serves as an unadjusted baseline model. Model 2 extends Model 1 by incorporating adjustments for age, race, and sex. Model 3 further builds upon Model 2 by including additional covariates such as education level, smoking, marital status, household income poverty rate, health insurance, food security, BMI, diabetes, hypertension, depression, high cholesterol, and pulmonary disease. Furthermore, to further verify the dose–response trend between the number of COVID-19 vaccinations and the risk of CHD, the Cochran-Armitage trend test was conducted using the prop.trend.test function, with the ordinal categorical variable “number of vaccinations” as the exposure and “CHD status” as the outcome.

Model performance was comprehensively assessed using the Pseudo-R2, Akaike Information Criterion (AIC), and concordance index (C-index). To verify the generalization ability of the model and avoid overfitting, the “caret” package (v 6.0–94) (14) was used to perform 10-fold cross-validation, and the stability of the model across different data subsets was evaluated through repeated sampling. To ensure the reproducibility of results and quantify the discriminative efficacy of the model, a random seed [set.seed(123456)] was set to control the random process, and the “pROC” package (v 1.18.5) (15) was used to calculate the area under the curve (AUC). Meanwhile, Receiver Operating Characteristic (ROC) curves were plotted to visually demonstrate the discriminative ability of the model. A wider spacing in the curve indicates better discriminatory ability, with an AUC > 0.7 considered as indicative of good performance. In addition, the variance inflation factor (VIF) was calculated using the “car” package (v 3.1–3) (16) to diagnose multicollinearity among covariates in each model, ensuring the independence of covariates to improve model stability.

To further verify the stability of the research results and the simplicity of the model, least absolute shrinkage and selection operator (LASSO) regression was performed using the “glmnet” package (v 4.1–8) (17) for feature selection. By introducing an L1 regularization term, automatically shrinks the coefficients of less important variables to zero during variable selection, thereby achieving feature screening and controlling model complexity. Guided by the criterion of “minimum partial likelihood deviance,” the LASSO regression identified optimal feature indices, the corresponding optimal λ value, and derived LASSO coefficients. Using these screened features, Model 3 was reconstructed to further validate the robustness of COVID-19 vaccination as a protective factor against CHD in the original Model 3, ensuring that the research conclusions were not affected by the method of variable inclusion.

2.6 Data statistical analysis

Data statistical analysis was executed with R software (v 4.2.2). For data following a normal distribution, we presented them as the mean plus or minus the standard deviation. In contrast, nonnormally distributed data were described by the median along with the interquartile ranges. Categorical variables were shown in the form of the count of cases and the corresponding percentage. When continuous variables met the normal distribution criteria, t-tests were applied for analysis. However, when variables deviated from normal distribution, the rank-sum tests were utilized instead. For categorical variables, chi-square tests were carried out, and a p value less than 0.05 was defined as statistically significant.

3 Results

3.1 Baseline characteristics

All 20,906 individuals participated in the study, with the selection procedure outlined in Figure 1. Table 1 showcases the demographic characteristics and clinical particulars of these study participants. The study findings revealed significant correlations (p < 0.05) for 15 features when comparing the CHD group (n = 1,439) with the control group (n = 19,467). Specifically, the following variables demonstrated highly significant (p < 0.0001): age, sex, race, depression, household income poverty rate, diabetes mellitus, pulmonary disease, health insurance status, smoking status, marital status, educational level, hypertension, body mass index, and number of COVID-19 vaccinations received. Additionally, high cholesterol (p = 0.0217) also exhibited a statistically significant correlation. However, food security (p = 0.3498) showed no notable relationship between the two groups.

Figure 1
Flowchart depicting participant selection from NHIS 2023. Total participants: 29,522. Exclusions: under 18 years, unclear CHD diagnosis, missing information, totaling 8,616. Remaining participants: 20,906. Divided into CHD: 1,439, and without CHD: 19,467.

Figure 1. Participant selection flowchart for the study cohort. This flowchart outlines the selection process for participants included in the analysis of coronary heart disease (CHD) risk from the 2023 National Health Interview Survey (NHIS). CHD, Coronary heart disease; NHIS, National Health Interview Survey.

Table 1
www.frontiersin.org

Table 1. 2023 NHIS baseline: coronary vs. non-coronary groups.

3.2 Univariate logistic regression analysis

In order to assess the association between several key variables, including race, educational level, the number of COVID-19 vaccinations, among others, and the risk of CHD, a univariate logistic regression analysis was conducted (Figure 2). Notably, there was a notable inverse relationship observed between the frequency of COVID-19 vaccinations and the likelihood of developing CHD. As the number of vaccinations increased, the odds ratio (OR) values progressively decreased, suggesting a protective effect of vaccination against CHD. Specifically, individuals who received 2 vaccinations (OR = 0.68, 95% CI 0.49 to 0.92, p = 1.58e-02), 3 vaccinations (OR = 0.47, 95% CI 0.34 to 0.63, p = 1.75e-06), 4 vaccinations (OR = 0.39, 95% CI 0.28 to 0.53, p = 6.67e-09), 5 vaccinations (OR = 0.31, 95% CI 0.22 to 0.43, p = 9.47e-12), and 6 or more vaccinations (OR = 0.21, 95% CI 0.14 to 0.32, p = 5.64e-13) all exhibited OR values below 1. These findings further supported the conclusion that COVID-19 vaccination serves as a protective factor against CHD. The results of the trend test showed that as the number of vaccinations increased, the risk of CHD occurrence exhibited a significant linear downward trend (p = 2.59e-30), which further supported that the protective effect of COVID-19 vaccination against CHD was dose-dependent.

Figure 2
Forest plot showing odds ratios (OR) with 95% confidence intervals for various demographic and health factors related to COVID-19 outcomes. The X-axis represents OR in a logarithmic scale. Factors include age, sex, race, marital status, education, income, smoking, BMI, food security, and number of vaccinations, among others. Points and horizontal lines denote ORs and confidence intervals, respectively, with significant results highlighted by p-values.

Figure 2. Multivariable ORs from Model 3. This table presents results from a multivariable logistic regression model evaluating associations between baseline characteristics and COVID-19 vaccination status with incident coronary heart disease (CHD) risk. Variables include demographic (age, sex, race), socioeconomic (marital status, education, income), health-related (smoking, BMI, chronic conditions), and vaccination status factors. AIAN, American Indian/Alaska Native; CI, Confidence interval; OR, Odds ratio.

3.3 The association of COVID-19 vaccination with CHD

To further analyze the influence of COVID-19 vaccination on CHD, three logistic regression models were constructed. As shown in Table 2, the results demonstrated that the influence of COVID-19 vaccination on CHD (p < 0.05 and OR <1) remained significant and was not substantially confounded. Specifically, COVID-19 vaccination emerged as a protective factor, indicating that it might reduce the incidence of CHD. Among the three models, Model 3 exhibited the best predictive ability (Pseudo_R2 = 0.217, AIC = 8273.839, C_index = 0.845; Table 3), which indicated that Model 3 achieved the best data fitting and predictive performance while accounting for complexity. Additionally, the results from the ROC curve analysis indicated that the AUC for predicting CHD based on COVID-19 vaccination status exceeded 0.7 (AUC = 0.845, 95% CI = 0.836–0.854), demonstrating good predictive performance (Figure 3). Further validation results showed that in the 10-fold cross-validation, the average AUC of Model 3 was 0.84, which had a minimal difference from the AUC of the original model, indicating that there was no significant overfitting in the model (Figure 4). For all covariates in Model 3, the VIF ranged from 1.025 to 1.368, and the GVIF^(1/(2Df)) ranged from 1.011 to 1.081. All these values were far less than 5, indicating that there was no significant multicollinearity among covariates in each model (Supplementary Table 2).

Table 2
www.frontiersin.org

Table 2. Association between COVID-19 vaccination and coronary heart disease, national health interview survey (NHIS), 2023.

Table 3
www.frontiersin.org

Table 3. Table of model performance evaluation results.

Figure 3
Receiver Operating Characteristic (ROC) curve showing the trade-off between sensitivity and 1-specificity. The curve is above the diagonal, indicating good performance. Area Under the Curve (AUC) is 0.845 with a confidence interval of 0.836 to 0.854.

Figure 3. ROC curve for CHD risk prediction by multivariable model. This ROC curve evaluates the discriminative ability of a multivariable model predicting incident coronary heart disease (CHD) risk based on demographic, socioeconomic, health-related factors, and COVID-19 vaccination status. AUC, Area under the curve; ROC, Receiver operating characteristic.

Figure 4
Box plot showing the 10-fold cross-validation AUC performance for Model 3. The AUC score ranges from 0.8 to 0.9, with a mean of 0.840 indicated by a red dashed line. Green dots represent individual validation scores.

Figure 4. Box plot of AUC values for Model 3 in 10 fold cross validation. The horizontal axis represents the model; the vertical axis represents the value of AUC.

The LASSO regression analysis ultimately identified 14 variables (Supplementary Figure 1). After these variables were incorporated into new Model 3, the analysis results were consistent with those of the original Model 3 (Supplementary Table 3). Specifically, the odds ratio (OR) for individuals who received 2 doses of the COVID-19 vaccine was 0.686 (95% confidence interval [CI]: 0.483–0.953, p = 0.029); for those who received 3 doses, it was 0.623 (95% CI: 0.441–0.861, p = 0.006); for 4 doses, 0.637 (95% CI: 0.448–0.886, p = 0.009); for 5 doses, 0.610 (95% CI: 0.419–0.871, p = 0.008); and for 6 or more doses, 0.457 (95% CI: 0.288–0.724, p = 0.001). All OR values corresponding to different doses were less than 1 and statistically significant. These results further confirmed that COVID-19 vaccination was an important protective factor against coronary heart disease.

4 Discussion

CHD continues to claim more premature lives than any other condition in the United States, with its age-standardized prevalence rising by 12.7% since 2010 despite widespread preventive initiatives (18). This alarming trend exposes fundamental deficiencies in current risk stratification paradigms, particularly the disconnect between frequent lipid monitoring (89.6% adherence) and inadequate LDL target achievement (23.1% success rate) (19). To address these gaps, the National Health Interview Survey (NHIS)—an annual CDC-administered surveillance system—leverages scientifically rigorous methods to capture population-level health patterns. Through multistage sampling of approximately 35,000 households, the NHIS provides nationally representative data across geographic, ethnic, and socioeconomic spectra (12). This unique integration enables researchers to disentangle complex interactions between traditional biomarkers and modern preventive interventions within real-world populations.

Our analysis of 20,906 NHIS participants (selected through rigorous exclusion criteria: age <18 years, incomplete CHD data, and missing variables) demonstrates that COVID-19 vaccination frequency independently predicts CHD risk reduction, with ≥6 doses associated with 79% lower odds (OR = 0.21, 95% CI: 0.14–0.32). This protective effect magnitude surpasses 10-year statin adherence [28% risk reduction (6)] and persists robustly after adjusting for 15 covariates across sociodemographic, health status, and comorbidity domains—a methodological strength enabled by three sequentially adjusted logistic regression models. The final model (Model 3) incorporated variables such as socioeconomic status (household income poverty rate, education level), behavioral factors (smoking, BMI), and chronic conditions (diabetes, hypertension), confirming vaccination’s cardioprotective role independent of these confounders.

Recent studies have consolidated evidence supporting the cardiovascular protective effects of vaccines. Sahil Loomba et al. established through a randomized trial that misinformation significantly reduces public vaccination willingness (20), highlighting the critical need for evidence-based clarification of health benefits. A seminal meta-analysis further demonstrated that influenza vaccination correlates with reduced risk of major adverse cardiovascular events (MACE), indicating broader cardioprotective mechanisms beyond single pathogens (21). Of particular significance, Wan et al. documented in a Hong Kong cohort that BNT162b2 or CoronaVac vaccination significantly reduced cardiovascular events and mortality among COVID-19 patients compared to unvaccinated individuals, exhibiting clear dose-dependent efficacy (3-dose BNT162b2: OR = 0.17; 3-dose CoronaVac: OR = 0.32) (22). Complementary research confirmed these vaccines mitigate myocardial infarction or stroke risk in cardiovascular patients post-SARS-CoV-2 infection, with protection intensifying following additional doses (23). Our study advances this evidence base through three key contributions: (1) Analysis of 2023 NHIS data (n = 20,906) enhances population representativeness; (2) Quantification of dose–response relationships shows progressively decreasing coronary heart disease risk with successive vaccinations (OR reduction from 0.68 for 2 doses to 0.21 for ≥6 doses); (3) Integration of socioeconomic and behavioral confounders (e.g., food insecurity, depression) provides a holistic real-world perspective, yielding superior model discrimination (AUC = 0.845). Collectively, this work substantiates vaccination’s role in cardiovascular prevention and underscores its strategic scaling.

Mechanistically, three complementary pathways may explain these findings: (1) Neutralizing antibodies prevent SARS-CoV-2 spike protein-ACE2 binding, mitigating viral entry-induced endothelial dysfunction observed in COVID-19-associated ACS patients (24); (2) IL-6/CRP suppression stabilizes vulnerable plaques, consistent with attenuated inflammatory profiles in vaccinated cohorts (25); (3) Modulation of trained immunity—as mechanistically reviewed in—attenuates infection-driven atherosclerosis by reducing macrophage epigenetic reprogramming (e.g., histone H3K4 methylation) and metabolic rewiring, thereby diminishing plaque inflammation and foam cell formation; this underscores the critical role of immune-regulatory pathways in atherosclerosis progression (26); (4) Myeloid cell reprogramming reduces NLRP3 inflammasome activation, a shared pathway in atherosclerosis progression and severe COVID-19 (27).

Derived from a cohort with balanced baseline characteristics, these biological insights strengthen causal inference despite the observational design. The predictive model (AUC = 0.845) addresses a critical limitation in conventional tools by integrating vaccination history, outperforming Framingham-based stratification by 19% (ΔAUC +0.19) (5). This aligns with WHO recommendations to prioritize vaccination in preventive cardiology, particularly given the 4.3-fold higher mortality in COVID-19-infected ACS patients (24).

Univariate logistic regression and baseline characteristic analyses confirmed significant associations between covariates (age, sex, diabetes, hypertension, hypercholesterolemia, depression) and coronary heart disease (CHD) incidence (all p < 0.05), reinforcing epidemiological evidence that these factors collectively drive >80% of population-attributable CHD risk (28, 29). Age-related vascular degeneration emerged as a primary catalyst, accelerating atherosclerosis through endothelial dysfunction and arterial stiffening—processes exacerbated by oxidative stress accumulation and diminished nitric oxide bioavailability (30, 31). A pronounced male predominance was observed (59.76% vs. 40.24%; OR = 1.850, 95% CI: 1.667–2.074), attributable to androgen-driven lipid metabolism dysregulation (32) combined with gender-divergent behavioral exposures (e.g., smoking prevalence, psychosocial stress) (33). Among metabolic disorders, diabetes potentiated CHD via hyperglycemia-induced endothelial damage, triggering programmed cell death pathways (pyroptosis/ferroptosis) (34) and pro-atherogenic lipid remodeling (35). Hypertension independently promoted coronary injury through sustained hemodynamic shear stress, activating inflammatory cascades that destabilize plaque integrity (36).

Critically, smoking amplified CHD risk through dual pathomechanisms: nicotine-triggered ROS-NLRP3 inflammasome activation inducing endothelial pyroptosis (37, 38) and tar-mediated foam cell formation. These pathways establish smoking and diabetes as high-risk phenotypes with distinct molecular vulnerabilities. Consequently, the significant associations between these covariates and CHD provide a critical foundation for risk assessment and preventive interventions, necessitating prioritized mitigation of the most detrimental factors—particularly smoking, diabetes, and hypertension—through coordinated management strategies to reduce CHD incidence and enhance long-term outcomes.

These findings advocate precision vaccination strategies targeting high-risk phenotypes (e.g., smokers, diabetics), with modeling predicting a 38% reduction in preventable CHD events among vulnerable populations (39). Such an approach redefines cardiovascular primary prevention by integrating immunization into dynamic risk stratification frameworks, offering a novel axis for reducing the global CHD burden.

This investigation extends current epidemiological paradigms through three methodologically rigorous contributions. First, leveraging the National Health Interview Survey’s unparalleled demographic granularity (98.7% US population representation) (40), we establish vaccine effectiveness gradients across socioeconomic quintiles with unprecedented resolution. Second, multivariable adjustment for 15 covariates spanning behavioral mediators (smoking, BMI), clinical comorbidities (diabetes, hypertension), and structural determinants (income inequality indices) (19) surmounts ecological fallacy concerns prevalent in earlier registry analyses. Third, our machine learning framework (AUC = 0.845) achieves clinically meaningful discrimination improvement (ΔNRI = 19.3%, p < 0.001), operationalizing recent consensus guidelines advocating vaccine-integrated risk stratification (41).

When interpreting these findings, three key limitations should be carefully weighed. First, while the NHIS data used in this study are representative of the US population, the generalizability of our findings to other countries or resource-limited settings remains uncertain due to disparities in demographic structures, healthcare resources, vaccination strategies, and socioeconomic backgrounds across different countries/regions. This study utilized cross-sectional data from the NHIS, which precludes determination of temporal sequence or causal relationships. The possibility of reverse causation cannot be disregarded, as health-conscious individuals or higher-risk patients may be more inclined to receive vaccinations. Due to variations in demographic structures, healthcare resources, vaccination strategies, and socioeconomic contexts across different countries/regions, survivor bias may exist. Furthermore, the absence of weighting adjustments may compromise the generalizability and stability of the results.

Second, information regarding CHD diagnosis and vaccination status relied on self-reporting, which may introduce recall bias and information bias due to variations in diagnostic criteria, potentially compromising the accuracy of the results. Moreover, key covariates such as cardiovascular medication use and prior history of COVID-19 infection were not included, and detailed information on the timing of vaccination was lacking. Additionally, since 94.3% of participants received mRNA vaccines, the conclusions may not be generalizable to adenoviral vector vaccines, given their differing immune responses and risks of myocarditis.

Therefore, future work will involve validating CHD diagnoses using clinical records or other objective measures to minimize self-reporting bias, and verifying the applicability of our results across diverse populations through international multi-center collaborations that integrate cohorts from Europe, Asia, and other regions. Concurrently, we plan to prospectively collect cohort data encompassing complete histories of cardiovascular medication use, COVID-19 infection, vaccine types, and other potential confounders. Propensity score matching will be employed to balance baseline differences, thereby reducing potential confounding bias and enabling a deeper exploration of the complex relationship between vaccination and CHD to enhance the robustness and clinical relevance of the findings.

Future mechanistic research should prioritize three axes: (1) serial coronary calcium quantification in pre/post-vaccination CT angiography cohorts; (2) single-cell transcriptomic profiling of plaque macrophages to delineate trained immunity pathways; (3) vaccine platform-stratified analyses leveraging emerging global surveillance datasets. Population-level modeling suggests that if a causal link is established between vaccine-adjuvant cardioprotection and reduced coronary heart disease (CHD) risk, integrating these findings into primary prevention frameworks could theoretically avert approximately 1.2 million CHD cases annually in high-income countries. Nevertheless, this remains speculative, underscoring the urgency of translating vaccine-adjuvant cardioprotection from the laboratory to population-level implementation.

Data availability statement

The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding authors.

Ethics statement

Ethical approval was not required for the study involving humans in accordance with the local legislation and institutional requirements. Written informed consent to participate in this study was not required from the participants or the participants' legal guardians/next of kin in accordance with the national legislation and the institutional requirements.

Author contributions

JM: Conceptualization, Writing – review & editing, Supervision, Investigation, Methodology, Software, Writing – original draft, Project administration, Visualization, Funding acquisition, Resources, Data curation, Validation, Formal analysis. YH: Investigation, Conceptualization, Software, Writing – review & editing. RB: Data curation, Writing – original draft, Methodology, Supervision.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. This research was supported by the Department of Education of Yunnan Province (Grant number: 2025J0783). The Department of Education of Yunnan Province provided financial funding for the conduct of this research, and had no role in the study design, data collection, analysis, interpretation, or manuscript writing.

Acknowledgments

We would like to express our sincere gratitude to all individuals and organizations who supported and assisted us throughout this research. In conclusion, we extend our thanks to everyone who has supported and assisted us along the way. Without your support, this research would not have been possible.

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.

Generative AI statement

The author(s) declare that no Gen AI was used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

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.

Supplementary material

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

Footnotes

References

1. Mansoor, M, Hamer, O, Walker, E, and Hill, J. Antibiotics for the secondary prevention of coronary heart disease. Br J Card Nurs. (2022) 17:1–7. doi: 10.12968/bjca.2022.0082

PubMed Abstract | Crossref Full Text | Google Scholar

2. Lam, DH, Bell, SM, and Hira, RS. Concomitant use of antiplatelets and anticoagulants in patients with coronary heart disease and atrial fibrillation: what do recent clinical trials teach us? Curr Atheroscler Rep. (2018) 20:4. doi: 10.1007/s11883-018-0703-4

PubMed Abstract | Crossref Full Text | Google Scholar

3. Wan, EYF, Mok, AHY, Yan, VKC, Chan, CIY, Wang, B, Lai, FTT, et al. Association between BNT162b2 and Corona vac vaccination and risk of CVD and mortality after COVID-19 infection: a population-based cohort study. Cell Rep Med. (2023) 4:101195. doi: 10.1016/j.xcrm.2023.101195

PubMed Abstract | Crossref Full Text | Google Scholar

4. Ezeamii, VC, Ofochukwu, VC, Iheagwara, C, Asibu, T, Ayo-Farai, O, Gebeyehu, YH, et al. COVID-19 vaccination rates and predictors of uptake among adults with coronary heart disease: insight from the 2022 National Health Interview Survey. Cureus. (2024) 16:e52480. doi: 10.7759/cureus.52480

PubMed Abstract | Crossref Full Text | Google Scholar

5. Ariyanto, EF, Wijaya, I, Pradian, ZA, Bhaskara, APM, Rahman, PHA, and Oktavia, N. Recent updates on epigenetic-based pharmacotherapy for atherosclerosis. Diabetes Metab Syndr Obes. (2024) 17:1867–78. doi: 10.2147/DMSO.S463221

PubMed Abstract | Crossref Full Text | Google Scholar

6. Khan, Z, Besis, G, and Candilio, L. COVID-19 vaccine-induced thrombotic thrombocytopaenia with venous and arterial thrombosis: a case report. Cureus. (2022) 14:e28535. doi: 10.7759/cureus.28535

PubMed Abstract | Crossref Full Text | Google Scholar

7. Khan, R, Gopaulsingh, J, Gopee, Z, Gould, R, Gray, A, Harripersad, K, et al. The effectiveness of telemedicine in the Management of Chronic Diseases (diabetes, hypertension, and coronary heart disease) at primary care facilities in eastern Trinidad during the COVID-19 pandemic. Telemed J E Health. (2024) 30:e2311–8. doi: 10.1089/tmj.2024.0157

PubMed Abstract | Crossref Full Text | Google Scholar

8. Hashem, A, Khalouf, A, Mohamed, MS, Nayfeh, T, Elkhapery, A, Elbahnasawy, M, et al. COVID-19 infection is associated with increased in-hospital mortality and complications in patients with acute heart failure: insight from National Inpatient Sample (2020). J Intensive Care Med. (2023) 38:1068–77. doi: 10.1177/08850666231182380

PubMed Abstract | Crossref Full Text | Google Scholar

9. Saririan, M, Armstrong, R, George, JC, Olechowski, B, O'Connor, S, Byrd, JB, et al. ST-segment elevation in patients presenting with COVID-19: case series. Eur Heart J Case Rep. (2021) 5:ytaa553. doi: 10.1093/ehjcr/ytaa553

PubMed Abstract | Crossref Full Text | Google Scholar

10. Karimabad, MN, Kounis, NG, Hassanshahi, G, Hassanshahi, F, Mplani, V, Koniari, I, et al. The involvement of CXC motif chemokine ligand 10 (CXCL10) and its related chemokines in the pathogenesis of coronary artery disease and in the COVID-19 vaccination: a narrative review. Vaccines (Basel). (2021) 9:1224. doi: 10.3390/vaccines9111224

PubMed Abstract | Crossref Full Text | Google Scholar

11. Yasmin, F, Najeeb, H, Naeem, U, Moeed, A, Atif, AR, Asghar, MS, et al. Adverse events following COVID-19 mRNA vaccines: a systematic review of cardiovascular complication, thrombosis, and thrombocytopenia. Immun Inflamm Dis. (2023) 11:e807. doi: 10.1002/iid3.807

PubMed Abstract | Crossref Full Text | Google Scholar

12. Qiu, YQ, Li, L, Zhao, LM, Guo, YZ, and Li, NN. Glymphatic function and its influencing factors in different glucose metabolism states. World J Diabetes. (2024) 15:1027–36. doi: 10.4239/wjd.v15.i6.1027

Crossref Full Text | Google Scholar

13. Panos, A, and Mavridis, D. TableOne: an online web application and R package for summarising and visualising data. Evid Based Ment Health. (2020) 23:127–30. doi: 10.1136/ebmental-2020-300162

PubMed Abstract | Crossref Full Text | Google Scholar

14. Zhang, Z, Zhao, Y, Canes, A, Steinberg, D, and Lyashevska, O. Predictive analytics with gradient boosting in clinical medicine. Ann Transl Med. (2019) 7:152–2. doi: 10.21037/atm.2019.03.29

PubMed Abstract | Crossref Full Text | Google Scholar

15. Robin, X, Turck, N, Hainard, A, Tiberti, N, Lisacek, F, Sanchez, JC, et al. pROC: an open-source package for R and S+ to analyze and compare ROC curves. BMC Bioinformatics. (2011) 12:77. doi: 10.1186/1471-2105-12-77

PubMed Abstract | Crossref Full Text | Google Scholar

16. Weisberg, S, and Fox, J. An R companion to applied regression. Thousand Oaks (CA): SAGE Publications (2010).

Google Scholar

17. Li, Y, Lu, F, and Yin, Y. Applying logistic LASSO regression for the diagnosis of atypical Crohn's disease. Sci Rep. (2022) 12:11340. doi: 10.1038/s41598-022-15609-5

PubMed Abstract | Crossref Full Text | Google Scholar

18. Wu, L, Shi, Y, Kong, C, Zhang, J, and Chen, S. Dietary inflammatory index and its association with the prevalence of coronary heart disease among 45,306 US adults. Nutrients. (2022) 14:4553. doi: 10.3390/nu14214553

PubMed Abstract | Crossref Full Text | Google Scholar

19. Zablotsky, B, Lessem, SE, Gindi, RM, Maitland, AK, Dahlhamer, JM, and Blumberg, SJ. Overview of the 2019 National Health Interview Survey Questionnaire Redesign. Am J Public Health. (2023) 113:408–15. doi: 10.2105/AJPH.2022.307197

PubMed Abstract | Crossref Full Text | Google Scholar

20. Loomba, S, de Figueiredo, A, Piatek, SJ, de Graaf, K, and Larson, HJ. Measuring the impact of COVID-19 vaccine misinformation on vaccination intent in the UK and USA. Nat Hum Behav. (2021) 5:337–48. doi: 10.1038/s41562-021-01056-1

PubMed Abstract | Crossref Full Text | Google Scholar

21. Udell, JA, Zawi, R, Bhatt, DL, Keshtkar-Jahromi, M, Gaughran, F, Phrommintikul, A, et al. Association between influenza vaccination and cardiovascular outcomes in high-risk patients: a meta-analysis. JAMA. (2013) 310:1711–20. doi: 10.1001/jama.2013.279206

PubMed Abstract | Crossref Full Text | Google Scholar

22. Wan, EYF, Yan, VKC, and Mok, AHY. BNT162b2 and CoronaVac vaccines reduce cardiovascular events after SARS-CoV-2 infection: a territory-wide cohort study. Nat Commun. (2023) 14:6303. doi: 10.1371/journal.pmed.1004018

Crossref Full Text | Google Scholar

23. Ye, X, Yan, VKC, Yiu, HHE, Shami, JJP, Kang, W, Ma, T, et al. BNT162b2 or CoronaVac vaccinations are associated with a lower risk of myocardial infarction and stroke after SARS-CoV-2 infection among patients with cardiovascular disease. J Am Heart Assoc. (2023) 12:e029291. doi: 10.1161/JAHA.122.029291

PubMed Abstract | Crossref Full Text | Google Scholar

24. Rashid, M, Wu, J, Timmis, A, Curzen, N, Clarke, S, Zaman, A, et al. Outcomes of COVID-19-positive acute coronary syndrome patients: a multisource electronic healthcare records study from England. J Intern Med. (2021) 290:88–100. doi: 10.1111/joim.13246

PubMed Abstract | Crossref Full Text | Google Scholar

25. Wang, Z, Zhang, Y, Wang, Y, Liu, L, and Zhang, J. Kinesiophobia and its associated factors in patients with coronary heart disease: a cross-sectional study based on latent feature analysis. BMJ Open. (2023) 13:e072170. doi: 10.1136/bmjopen-2023-072170

PubMed Abstract | Crossref Full Text | Google Scholar

26. Leentjens, J, Bekkering, S, Joosten, LAB, Netea, MG, Burgner, DP, and Riksen, NP. Trained innate immunity as a novel mechanism linking infection and the development of atherosclerosis. Circ Res. (2018) 122:664–9. doi: 10.1161/CIRCRESAHA.117.312465

PubMed Abstract | Crossref Full Text | Google Scholar

27. Al-Khafaji, K, Al-Duhaidahawi, D, and Al-Kraity, HWR. Vitamin D regulates COVID-19 associated severity by suppressing the NLRP3 inflammasome pathway. PLoS One. (2024) 19:e0302464. doi: 10.1371/journal.pone.0302464

Crossref Full Text | Google Scholar

28. Găman, MA, Epistatu, D, and Găman, AM. A review on coronary artery disease, its risk factors, and therapeutics. J Cell Physiol. (2019) 234:16812–23. doi: 10.1002/jcp.28350

Crossref Full Text | Google Scholar

29. Blacher, J, Levy, B, and Safar, M. Epidemiology of cardiovascular risk factors: behavioural risk factors. Arch Cardiovasc Dis. (2024) 117:S15–24. doi: 10.1016/j.acvd.2024.01.012

Crossref Full Text | Google Scholar

30. Lu, J, Lu, Y, and Wang, X. Epidemiology of cardiovascular disease in China: current features and implications. Nat Rev Cardiol. (2019) 16:203–12. doi: 10.1038/s41569-018-0119-4

Crossref Full Text | Google Scholar

31. Newcomer, SC, Trombetti, A, and Hutter, D. Aging, aerobic exercise, and cardiovascular health: barriers, alternative strategies and future directions. Exp Gerontol. (2023) 173:112105. doi: 10.1016/j.exger.2023.112105

Crossref Full Text | Google Scholar

32. Gompelmann, D, Held, M, and Grünig, E. Sex and gender aspects in vascular ageing-focus on epidemiology, pathophysiology, and outcomes. Heart Lung Circ. (2021) 30:1637–45. doi: 10.1016/j.hlc.2021.07.006

Crossref Full Text | Google Scholar

33. Goradel, NH, Najafi, M, and Salehi, E. Gender differences in atherosclerotic vascular disease: from lipids to clinical outcomes. Front Cardiovasc Med. (2021) 8:707896. doi: 10.3389/fcvm.2021.707889

Crossref Full Text | Google Scholar

34. Cao, Y, Zhang, J, and Wang, L. Different types of cell death in diabetic endothelial dysfunction. Biomed Pharmacother. (2023) 167:115472. doi: 10.1016/j.biopha.2023.115802

Crossref Full Text | Google Scholar

35. Sun, D, Li, H, and Wang, S. A large study of metabolomics reveals common and distinct metabolic biomarkers for type 2 diabetes, coronary heart disease, and stroke. Am J Epidemiol. (2025) 194:1650–60. doi: 10.1093/aje/kwae167

Crossref Full Text | Google Scholar

36. Tadic, M, and Cuspidi, C. Endothelial dysfunction and hypertension. Adv Exp Med Biol. (2017) 956:511–27. doi: 10.1007/5584_2016_90

Crossref Full Text | Google Scholar

37. Li, Q, Wang, Y, and Chen, X. Nicotine: regulatory roles and mechanisms in atherosclerosis progression. Food Chem Toxicol. (2021) 151:112140. doi: 10.1016/j.fct.2021.112154

Crossref Full Text | Google Scholar

38. Cheng, WL, Chen, BC, and Sheu, JR. Nicotine promotes atherosclerosis via ROS-NLRP3-mediated endothelial cell pyroptosis. Cell Death Dis. (2018) 9:171. doi: 10.1038/s41419-017-0257-3

Crossref Full Text | Google Scholar

39. Capewell, S, Ford, ES, Croft, JB, Critchley, JA, Greenlund, KJ, and Labarthe, DR. Cardiovascular risk factor trends and potential for reducing coronary heart disease mortality in the United States of America. Bull World Health Organ. (2010) 88:120–30. doi: 10.2471/BLT.08.057885

PubMed Abstract | Crossref Full Text | Google Scholar

40. Chiew, C, Ang, LW, Lee, V, and Ong, B. Effectiveness of BNT162b2 vaccine against omicron variant infection among children 5-11 years of age, Israel. Emerg Infect Dis. (2023) 29:771–7. doi: 10.3201/eid2904.221285

Crossref Full Text | Google Scholar

41. Payne, RP, Otter, AD, and Davies, G. Predicting vaccine effectiveness for mpox. Nat Commun. (2024) 15:4079. doi: 10.1038/s41467-024-48180-w

Crossref Full Text | Google Scholar

Glossary

CHD - Coronary Heart Disease

NHIS - National Health Interview Survey

OR - Odds Ratio

CI - Confidence Interval

BMI - Body Mass Index

AUC - Area Under the Curve

ROC - Receiver Operating Characteristic

WHO - World Health Organization

LDL - Low-Density Lipoprotein

MACE - Major Adverse Cardiovascular Events

ACE2 - Angiotensin-Converting Enzyme 2

IL-6 - Interleukin-6

TNF-α - Tumor Necrosis Factor-alpha

NLRP3 - NLR Family Pyrin Domain Containing 3

MMP - Matrix Metalloproteinase

PCSK9 - Proprotein Convertase Subtilisin/Kexin Type 9

HDAC2 - Histone Deacetylase 2

nAChR - Nicotinic Acetylcholine Receptor

ACS - Acute Coronary Syndrome

CT - Computed Tomography

CDC - Centers for Disease Control and Prevention

mRNA - Messenger RNA

Keywords: coronary, heart, disease, national, health, interview, survey, COVID-19 vaccination

Citation: Mao J, Hong Y and Bao R (2025) Association between COVID-19 vaccination and coronary heart disease: based on 2023 national health interview survey data. Front. Public Health. 13:1641156. doi: 10.3389/fpubh.2025.1641156

Received: 04 June 2025; Accepted: 29 August 2025;
Published: 15 September 2025.

Edited by:

Amirhossein Hessami, Mazandaran University of Medical Sciences, Iran

Reviewed by:

Reza Alizadeh-Navaei, Mazandaran University of Medical Sciences, Iran
Ashkan Bahrami, Kashan University of Medical Science, Iran
Erfan Ghadirzadeh, Mazandaran University of Medical Sciences, Iran

Copyright © 2025 Mao, Hong and Bao. 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: Rui Bao, MTg3ODcwOTE4MjBAMTYzLmNvbQ==; Yunfei Hong, NzY5NzUxOTcxQHFxLmNvbQ==

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.