Population based cohort
The Atherosclerosis Risk in Communities Study (ARIC) study is an observational, population-based cohort recruiting adults from 4 US communities: Forsyth County, North Carolina; Washington County, Maryland; Jackson, Mississippi; and Minneapolis, Minnesota. Study participants have been prospectively followed since the baseline visit (1987 to 1989)(12). From 2005 through 2014, medical records for cohort residents ≥55 years of age hospitalized with heart failure were abstracted and formed a HF community surveillance study (13, 14).
Demographics (age, sex, smoking status, BMI, systolic and diastolic blood pressure, heart rate), medical histories (history of hypertension, diabetes), and laboratory values (BNP, hemoglobin) were collected from the hospital record by trained abstractors. For the purposes of this analysis, the last laboratory values were analyzed as suggested in previous study (13). Left ventricular ejection fraction (EF) was prioritized first for inpatient transthoracic echocardiography reports and replaced with EF obtained from other imaging methods if even the within 2-year transthoracic echocardiography reports was not available. Patients were defined as statin users if they answered “yes” for the question of “Lipid Lowering Statins At hospital discharge” and non-users if they answered “no”. 28-day and 1-year mortality were considered as the outcome of the population study. All cause deaths within 28 days and within 1 year were ascertained by linking patients records within the national Death Index. NICM herein was referred to idiopathic or dilated cardiomyopathy.
Initially, 496 patients were identified as having idiopathic or dilated cardiomyopathy. Patients were excluded if comorbiding with ischemic cardiomyopathy (n=66), with missing information on statins use (n=17) or without outcome data for 28-day or 1-year mortality (n=41).
All research activities are approved by local institutional review boards from the 4 ARIC communities and all participants were received the written informed consent from the original study. Informed consent was not required as all data used in this study were anonymized.
Statistical analysis
Baseline characteristics were compared between statin users and non-users. Data were presented as mean± SD (standard deviation) or median ± IQR (interquartile range) for continuous variables and number (percentage) for categorical variables. The Kruskal-Wallis or student t-test were used for continuous variables and Chi-square test for categorical variables.
Kaplan-Meier estimates were used to evaluate the 28-day and 1-year mortality risk between statin users and non-users. In addition, logistic regression analysis was applied to calculate the odds ratios (ORs) and 95% confidence intervals (CI) of the mortality risk with statin use. We constructed two models. One was not adjusted (model 1) and the other was adjusted for age, sex, history of hypertension, diabetes, body mass index, systolic blood pressure, diastolic blood pressure and smoking status (model 2).
Mendelian Randomization Analysis
MR uses publicly available summary-level GWAS data. Data for statins SNPs were obtained through UK Biobank, with 462933 individuals’ information of statin use (atorvastatin, rosuvastatin and simvastatin, separately). Data for NICM were obtained from FinnGen GWAS, including 11400 cases and 175752 controls of predominantly European ancestry. All data sources used (exposure from UK Biobank; outcome from FinnGen) have existing approvals from their respective Institutional Review Boards.
We conducted five MR methods, including inverse variance weighted (IVW), MR-Egger, weighted median, simple mode weighted mode to investigate whether potential causality between statin therapy and NICM exist. Effect size estimates (beta) and standard errors (SEs) for SNP-exposure and SNP-outcome associations from independent samples are used to estimate the causal exposure-outcome association. All MR studies rely on 3 basic assumptions as described in previous studies (11). Each instrumental variable (IV) was constructed from SNPs showing GWAS significant association (P<5×10–8) with the respective trait. All included SNPs were independent with the linkage disequilibrium (r2=0.001, kb=10,000). Palindromic SNPs with intermediate allele frequencies were removed during analyses. Leave-one-out analysis was applied to estimate the effect of one single SNP. Heterogeneity between IVs was tested using the Cochran’s Q statistic. MR-Egger intercept test was used to examine any possible horizontal pleiotropic effects.
All statistical analyses were performed using R version 4.1.2 or Stata 15.1 (StataCorp/SE, College Station, TX). “TwoSampleMR” R package was used for MR analyses (15).