Skip to main content

Rare damaging CCR2 variants are associated with lower lifetime cardiovascular risk

Abstract

Background

Previous work has shown a role of CCL2, a key chemokine governing monocyte trafficking, in atherosclerosis. However, it remains unknown whether targeting CCR2, the cognate receptor of CCL2, provides protection against human atherosclerotic cardiovascular disease.

Methods

Computationally predicted damaging or loss-of-function (REVEL > 0.5) variants within CCR2 were detected in whole-exome-sequencing data from 454,775 UK Biobank participants and tested for association with cardiovascular endpoints in gene-burden tests. Given the key role of CCR2 in monocyte mobilization, variants associated with lower monocyte count were prioritized for experimental validation. The response to CCL2 of human cells transfected with these variants was tested in migration and cAMP assays. Validated damaging variants were tested for association with cardiovascular endpoints, atherosclerosis burden, and vascular risk factors. Significant associations were replicated in six independent datasets (n = 1,062,595).

Results

Carriers of 45 predicted damaging or loss-of-function CCR2 variants (n = 787 individuals) were at lower risk of myocardial infarction and coronary artery disease. One of these variants (M249K, n = 585, 0.15% of European ancestry individuals) was associated with lower monocyte count and with both decreased downstream signaling and chemoattraction in response to CCL2. While M249K showed no association with conventional vascular risk factors, it was consistently associated with a lower risk of myocardial infarction (odds ratio [OR]: 0.66, 95% confidence interval [CI]: 0.54–0.81, p = 6.1 × 10−5) and coronary artery disease (OR: 0.74, 95%CI: 0.63–0.87, p = 2.9 × 10−4) in the UK Biobank and in six replication cohorts. In a phenome-wide association study, there was no evidence of a higher risk of infections among M249K carriers.

Conclusions

Carriers of an experimentally confirmed damaging CCR2 variant are at a lower lifetime risk of myocardial infarction and coronary artery disease without carrying a higher risk of infections. Our findings provide genetic support for the translational potential of CCR2-targeting as an atheroprotective approach.

Background

Atherosclerotic cardiovascular disease (CVD) is the leading cause of morbidity and mortality worldwide [1,2,3]. Over 20 years of preclinical research have provided overwhelming evidence for a causal role of inflammation in atherogenesis [4, 5] and recent trials provided proof-of-concept that targeting inflammation can lead to reductions in adverse cardiovascular events [6,7,8]. The canakinumab anti-inflammatory thrombosis outcome study (CANTOS) demonstrated that treatment with a monoclonal antibody against IL-1β lowers risk of recurrent vascular events among individuals with recent myocardial infarction [6]. The colchicine cardiovascular outcomes trial (COLCOT) [7] and the low-dose colchicine-2 (LoDoCo2) trial [8] further showed that colchicine, an established drug with widespread inhibitory effects on inflammatory pathways, [9, 10] lowers the risk of recurrent vascular events in patients with coronary artery disease (CAD). Targeting inflammation for atheroprotection must be balanced against the impact on any host defense responses. For example, both canakinumab [6] and colchicine [7] were associated with adverse effects including fatal infections in the CANTOS and the colchicine trials. While translational efforts have mostly focused on the inflammasome-IL-1β/IL-6 axis, [11] evidence from preclinical studies and early-phase clinical trials highlights the promise of alternative cytokines [5] for the development of a second generation of atherosclerosis-centered anti-inflammatory treatments [4].

CC-motif chemokine ligand 2 (CCL2) is a pivotal inflammatory chemokine regulating monocyte trafficking [12] that has been studied as a potential target in atherosclerosis. Preclinical data suggest that pharmacological targeting of CCL2 or its receptor CCR2 might lower atherosclerosis burden in experimental models [13]. Furthermore, deleting Ccr2 [14, 15] or Ccl2 [16] in atheroprone apolipoprotein E-deficient (Apoe−/−) or LDL receptor-deficient (Ldlr−/−) mice has been found to lead to reduction of atherosclerotic lesions and decreased intraplaque macrophage accumulation. Conversely, overexpression of Ccl2 accelerated atherosclerosis in Apoe−/− mice [17]. However, only recently large-scale genetic and epidemiological studies have highlighted the relevance of the CCL2/CCR2 pathway in human CVD, calling for clinical translation of strategies targeting this pathway [18]. Prospective observational studies [19, 20] support associations of higher circulating CCL2 levels with ischemic stroke and cardiovascular death, whereas Mendelian randomization analyses from population genetic studies [21, 22] show associations with higher risk of ischemic stroke and coronary artery disease. Furthermore, CCL2 levels are higher in atherosclerotic lesions from patients with symptomatic carotid stenosis, as compared to asymptomatic disease and are associated with features of plaque vulnerability [23].

Although these studies support a role of the CCL2/CCR2 axis in human atherosclerosis, it remains unclear whether pharmacological intervention in this pathway could lead to atheroprotection in humans. Several molecules targeting CCR2 are currently under development for autoimmune disease, liver disease, and cancer and could be repurposed for prevention of atherosclerotic cardiovascular disease [18]. Studies examining the phenotypic effects of rare genetic variants in population-based studies have been instrumental in predicting the consequences of pharmacological interventions [24,25,26,27,28] and might thus serve as a validation step for drug targets under development.

Here, we leveraged whole-exome sequencing (WES) data from 454,775 participants of the population-based UK Biobank (UKB) study to explore whether rare damaging variants in the CCR2 gene are associated with cardiovascular disease risk, burden of atherosclerosis, and traditional vascular risk biomarkers. We experimentally confirmed the damaging effects of one variant (M249K), which has a frequency of 0.15% among UKB participants of European ancestry and replicated its effects on cardiovascular risk in external population- and hospital-based biobanks. Finally, we performed a phenome-wide association study across two biobanks to explore associations of this damaging variant with the risk of infections or other potential safety signals of any CCR2-targeting treatments.

Methods

Study population

We used data from the UKB, a population-based prospective cohort study of UK residents aged 40–69 years recruited between 2006 and 2010 from 22 assessment centers across the UK [29]. This analysis was restricted to 454,775 out of 502,419 participants with available whole-exome sequencing data (UKB Exome 450k release from October 2021). Primary and secondary analyses were performed with an updated Functional Equivalence (FE) protocol that retains original quality scores in the CRAM files (referred to as the OQFE protocol). We included only variants that met published criteria [30]: individual and variant missingness < 10%, Hardy–Weinberg equilibrium p value > 10−15, minimum read coverage depth of 7 for SNPs and 10 for indels, at least one sample per site passed the allele balance threshold > 0.15 for SNPs and 0.20 for indels. We used genotype array data released by the UKB study to assign individuals to continental ancestry super-groups (African (AFR), Hispanic or Latin American (HLA, originally referred to as “AMR” by the 1000 Genomes Project), East Asian (EAS), European (EUR), and South Asian (SAS)) by projecting each sample onto reference principal components (PCs) calculated from the 1000 Genomes reference panel. We included only individuals without evidence of relatedness within the UKB samples, as defined by a KING cut-off of < 0.084.

Detection of damaging genetic variants in CCR2 and experimental validation

Restricting our analyses to the CCR2 gene, we detected predicted LoF or damaging missense variants with a MAF < 1% in the CCR2 exonic region and then explored which of those are associated with monocyte count as a functional readout. Variants from WES were annotated as previously described [31] using VEP v101 [32, 33]. We used the LofTee plugin for predicting LoF variants [34] and a REVEL cutoff of > 0.5 from dbNSFP version 4.0a for predicted damaging missense mutations.

We aimed to further restrict our analyses to damaging variants with proven functional effects. One of the key functions of the CCL2 axis is the recruitment of classical monocytes from the bone marrow to the circulation in a CCR2-dependent way [35,36,37]. Consequently, damaging variants in the CCR2 gene would be expected to be associated with a lower monocyte count in the circulation. To test associations with the circulating counts of monocytes and other white blood cells (WBC), absolute counts were extracted from the UKB fields 30,160 (basophil counts), 30,150 (eosinophil counts), 30,130 (monocyte counts), 30,120 (lymphocyte counts), and 30,140(neutrophil counts). Distributions were visually checked for normal distribution and log-normalized when needed. Associations were tested using regenie v2.2.4 [38]. For WBC analyses, we used sex, age at blood draw, and the first 5 ancestral PCs as covariates. The mixed model parameters were estimated using 200,000 genotyped common variants. Saddle point approximation regression was applied. Genetic variants associated with monocyte count at a false discovery rate (FDR)-corrected p value < 0.05 were prioritized for further experimental validation. To test whether the associations are specific for monocyte count, we also tested associations with other WBC counts (basophils, eosinophils, lymphocytes, neutrophils) for each variant.

Two predicted damaging CCR2 variants showing a significant association with lower monocyte count were brought forward to experimental testing. The two variants were M249K (3:46,358,273:T:A) and W165S (3:46,358,021:G:C, both hg38). First, we tested how human CCR2-knockout monocytic THP-1 cells transfected with these variants or wild-type CCR2 respond to chemoattraction by CCL2 using a Transwell migration assay (details in Additional file 1: Supplementary Methods). Next, we tested whether M249K-mutant CCR2 (which showed an effect in the migration assay) influences cAMP production in response to CCL2 in transfected human HEK293T cells (details in Additional file 1: Supplementary Methods).

Additionally, a widely used and versatile flexible protein–protein docking tool, HADDOCK 2.4 [39] was utilized to re-dock CCR2 in its wild-type and mutant (M249K) form with the CCL2 ligand in both the monomeric and homodimeric form. The input PDB files were 1dok and 5t1a. Following residues were suggested as active residues after literature research: 1dok: 2, 3, 4, 5, 6, 7, 8, 9, 10, 17, 18, 19, 20, 23, 24, 25, 34, 35, 36, 37, 47, 48, 49, 50. 5t1a: 37, 38, 39, 40, 41, 42, 101, 102, 103, 104, 105, 106, 107, 108, 109, 110, 111, 112, 113, 114, 179, 180, 181, 182, 183, 184, 185, 186, 187, 188, 189, 190, 191, 192, 193, 194, 195, 196, 197, 198, 199, 200, 201, 202, 203, 204, 205, 206, 269, 270, 271, 272, 273, 274, 275, 276, 277, 278, 279, 280, 281, 282, 283, 284, 285. We used standard parameters and the resulting HADDOCK score as an outcome measure of CCL2-CCR2 interaction.

Associations with cardiovascular risk, atherosclerosis burden, and vascular risk factors

Focusing on two sets of variants, we then explored associations of CCR2 variants with clinical cardiovascular endpoints including myocardial infarction, coronary artery disease, acute ischemic stroke, peripheral artery disease, and abdominal aortic aneurysm using ICD-10- and ICD-9-coded diagnoses, OPCS4-coded procedures, self-report, and algorithmically defined phenotypes provided by the UKB as detailed in Additional file 2: Table S1. As our hypothesis was that these variants would lead to reductions in risk by lowering the lifetime burden of atherosclerosis, we also constructed a combined phenotype of atherosclerosis manifestations in four vascular beds (coronary arteries, cerebrovascular system, peripheral arteries of the extremities, and aortic atherosclerosis). Severe atherosclerosis was defined as presence of clinical manifestations in at least two vascular beds (Additional file 2: Table S1).

We followed two approaches to test associations with these phenotypes: (i) we performed a burden test combining all predicted damaging variants (LoF or REVEL > 0.5) within CCR2; (ii) focusing on M249K, the only damaging monocyte-lowering CCR2 variant that was experimentally validated in the migration and cAMP assays, we performed logistic regression analysis adjusted for age, sex, and the first 5 ancestral PCs testing. We included both prevalent and incident endpoints as outcomes in all analyses and used Firth’s correction for unbalanced case/control ratios in our logistic regression analysis. To explore whether the associations with cardiovascular risk are mediated through known vascular risk factors, we further tested associations of the damaging genetic variants with the following phenotypes: systolic and diastolic blood pressure, circulating LDL and HDL cholesterol, circulating apolipoprotein B levels, circulating glycated hemoglobin A1c (HbA1c) concentration, body mass index, and C-reactive protein. Individuals under antihypertensive medications were excluded in the analyses for systolic and diastolic blood pressure, individuals under lipid-lowering medications were excluded from the analyses for LDL and HDL cholesterol and apolipoprotein B levels, and individuals under glucose-lowering treatments were excluded from the analyses for HbA1c.

Validation of M249K CCR2 variant in external datasets and meta-analysis

For external validation, we obtained summary statistics from six different data sources based on a pre-defined protocol: the TransOmics and Precision Medicine Program supported by NHLBI (TOPMed) Program including multiple studies and a trans-ancestry population in the USA (n = 51,732), [40] the population-based deCODE dataset in Iceland (n = 345,992), [41] the hospital-based Penn Medicine Biobank (PMBB, n = 43,721), [42] the Million Veteran Program (MVP, n = 438,905), [43] the Geisinger DiscovEHR-MyCode cohort in central and northeastern Pennsylvania (n = 145,826), [44] as well as the hospital-based Mass General Brigham Biobank (MGBB, n = 36,419) [45]. Specifically, we requested data from logistic regression models adjusted for age, sex, race, study-specific variables (e.g., sequencing center), and the first 10 PCs (with the Firth’s correction) for the M249K variant. The variant was either assessed by sequencing or directly genotyped in all cohorts. Details about the individual cohorts are provided in Additional file 1: Supplementary Methods. The derived odds ratios (OR) from the six datasets were meta-analyzed using fixed-effects and random-effects meta-analyses and were subsequently also meta-analyzed with the results from the UKB. We set a significance threshold of p < 0.05 in the replication meta-analysis. Heterogeneity was assessed with the I [2] and the Cochran Q statistic.

Phenome-wide association study

To explore potential adverse effects associated with damaging CCR2 genetic variants, we tested associations of M249K with the full range of clinical phenotypes encoded in the UKB and replicated in an external cohort (Geisinger DiscovEHR-MyCode) [44]. We used the Phecode Map 1.2 to map UKB ICD10-codes to phecodes [46] using all ICD10 codes (main position, secondary position, death records) from the UKB. We excluded phecodes with < 100 cases and phecodes that are male- or female-specific. Individuals were assigned a case status if > 1 ICD10 code was mapped to the respective phecode. To approximate effect size in a logistic regression framework, we used minor allele carrier status as an independent variable and age at baseline, sex, and 5 ancestry PCs as covariates. We used Firth’s correction for unbalanced case/control ratios in our logistic regression analysis for all results with p < 0.05. The results from the two cohorts were meta-analyzed using fixed- and random-effects models. Heterogeneity was assessed with the I [2] and the Cochran Q statistic.

Results

Computationally predicted damaging CCR2 variants

Among 428,191 unrelated (out of 454,775 with whole-exome sequencing data) UK Biobank participants from the whole-exome sequencing data release, we found a total of 45 predicted LoF or damaging (REVEL > 0.5) variants in the exonic region of the CCR2 gene distributed across 787 heterozygous carriers (frequency 0.18%, Fig. 1 and Additional file 2: Table S2). There was no homozygous carrier of any of these variants and variants were predominantly prevalent in individuals of European ancestry (779 carriers, frequency 0.20%), as compared to individuals of African, Hispanic or Latin American, East Asian, and South Asian ancestry (total of 8 carriers, pooled frequency 0.02%). Due to the very low frequency in other ancestries, we restricted our analyses to individuals of European ancestry (n = 393,416), Baseline characteristics of the study population are presented in Additional file 2: Table S3.

Fig. 1
figure 1

Damaging CCR2 variants. Domain structure of the CCR2 protein and position of the predicted loss-of-function (LoF) or missense damaging (REVEL > 0.5) variants present in > 2 UK Biobank participants in the CCR2 exonic regions

Associations of CCR2 variants with cardiovascular endpoints in gene-burden tests

In a discovery analysis including all 45 computationally predicted LoF or damaging variants, we next explored whether genetic variation in CCR2 is associated with atherosclerotic disease using a gene-based burden test. While far from genome-wide statistical significance for gene-based testing (p < 2.7 × 10−6), we found nominally significant associations of the predicted damaging CCR2 variants with myocardial infarction (OR: 0.60, 95%CI: 0.40–0.90, p = 0.008) and coronary artery disease (OR: 0.76, 95%CI: 0.59–0.99, p = 0.03), as well as directionally consistent associations with the odds of all other examined outcomes (ischemic stroke, peripheral artery disease, abdominal aortic aneurysm, Fig. 2A). Variant burden was associated with lower risk of a combined atherosclerotic endpoint (OR: 0.80, 95%CI: 0.65–0.98, p = 0.03), as well as with severe atherosclerotic disease, defined by clinical manifestations in at least two vascular beds (OR: 0.24, 95%CI: 0.07–0.82, p = 0.003; Fig. 2B). We found a trend for a dose–response pattern with lower frequency of clinically manifest atherosclerotic disease across the number of vascular beds involved among carriers for these CCR2 variants (frequency of atherosclerosis presence among carriers vs. non-carriers: in no vascular bed 88.0% vs. 90.7%; in 1 vascular bed 10.5% vs. 9.1%; in 2 vascular beds 1.2% vs. 0.3%; in 3 vascular beds 0.3% vs. 0.1%, in 4 vascular beds 0.03% vs. 0%; OR from ordinal regression: 0.74, 95%CI: 0.54–0.97, p = 0.01).

Fig. 2
figure 2

Associations of rare computationally predicted damaging CCR2 variants with risk of atherosclerotic cardiovascular disease in the UK Biobank. Associations of 45 predicted loss-of-function (LoF) or damaging (REVEL > 0.5) CCR2 variants (minimum allele frequency < 0.01) in burden tests with A risk of cardiovascular disease endpoints and B a combined atherosclerosis endpoint as well as severe clinical atherosclerotic disease defined by manifestations in at least 2 vascular beds (coronary arteries, cerebrovascular system, peripheral arteries of extremities, aorta) among UK Biobank participants of European ancestry. C Prevalence of clinically manifest atherosclerosis across 4 vascular beds (coronary arteries, cerebrovascular system, peripheral arteries of extremities, aorta) among carriers and non-carriers of the 45 predicted LoF or damaging CCR2 variants. The odds ratio (OR) is derived from ordinal regression adjusted for age, sex, and the first 5 ancestral principal components

Functional consequences of CCR2 damaging variants

To move beyond computationally predicted functional effects, we followed a two-step approach to select CCR2 variants with functional relevance. Since one of the key functions of CCL2 is the CCR2-dependent recruitment of classical monocytes from the bone marrow to the circulation, we first tested the association of the predicted damaging variants in CCR2 with circulating monocyte counts. We found two variants (M249K and W165S) to be associated with lower monocyte counts (FDR-corrected p < 0.05, Additional file 2: Table S2) and prioritized them for experimental validation.

We tested whether transfection of human CCR2-knockout THP-1 monocytes with either of these two variants influences directed migration towards CCL2 in a trans-migration cell assay (Fig. 3A). Cells transfected with M249K showed a clearly reduced response to CCL2, when compared to cells transfected with wild-type CCR2, resembling the response of CCR2-knockout cells (Fig. 3B and Additional file 1: Fig. S1). Cells transfected with W165S showed a profile similar to cells transfected with wild-type CCR2 (Additional file 1: Fig. S1). HEK293T cells transfected with M249K further showed a profound reduction of downstream activity (increase in cAMP) in response to CCL2, when compared to HEK293T cells transfected with wild-type CCR2 (Fig. 3C). We found no difference in protein levels of CCR2 in lysates of HEK293T cells transfected with the M249K mutant vs. wild-type CCR2, suggesting no impact of the M249K variant on expression or stability of CCR2 (Additional file 1: Fig. S2). Flow cytometry revealed no difference of CCR2 expression in the surface of the HEK293T cells transfected with the M249K mutant or wild-type CCR2, thus implying no impact on localization of the protein (Additional file 1: Fig. S2). M249K-transfected HEK293T cells further showed decreased ERK1/2 phosphorylation as a result of exposure to increasing concentration of CCL2, suggesting reduced MAP kinase signaling activity that is downstream to beta-arrestin signaling [47, 48] (Additional file 1: Fig. S2). Experiments in both THP-1 monocytes and HEK293T cells showed a bell-shaped curve response to increasing concentrations of CCL2, a phenomenon that is well-characterized for chemokine receptor signaling in the literature [49,50,51].

Fig. 3
figure 3

Validation of the damaging effect of M249K (3:46,358,273:T:A) on CCL2-driven monocyte chemotaxis. A Experimental outline of chemotaxis assay. Transfected CCR2 − / − THP-1 cells (marked in green) and non-transfected cells (marked in red) were used in a Transwell chemotaxis assay. The cells that moved to the lower chamber were collected, stained, and quantified using flow cytometry. B Chemotaxis in response to increasing concentrations of CCL2 for human CCR2-knockout monocytic THP-1 cells transfected with the mutant M249K vs. wild-type CCR2, as determined in a trans-well migration assay (comparisons derived from two-way ANOVA, ns: p > 0.05, ***p < 0.001, ****p < 0.0001). C Results of cyclic AMP (cAMP) assay. Shown is the cAMP activity in HEK293T cells transfected with either empty vector or wild-type or M249K CCR2 in response to different concentrations of CCL2. Results are presented as “Ratio 665/620 × 10,000” (ratio of fluorescence at 665 nm and 620 nm × 10,000). D Associations of M249K with counts of different leukocyte populations in the subset of European ancestry participants of the UK Biobank (N = 393,838), as derived from linear regression models adjusted for age, sex, and the first 5 ancestral principal components

M249K has a frequency of 0.15% among UKB participants of European ancestry and leads to the replacement of methionine by lysine in the sixth transmembrane domain of the CCR2 receptor (Fig. 1). Confirming the specificity of the consequences of M249K on monocyte recruitment, we found no evidence of associations with any of the other leukocyte type counts in the UKB (Fig. 3D). Exploring whether mutant (M249K) CCR2 would impact the interaction between CCR2 and CCL2 with HADDOCK 2.4, docking complexes between WT CCR2 and CCL2 showed a more favorable HADDOCK score than between M249K CCR2 and CCL2 for both monomeric and homodimeric CCL2 (dimer: WT: − 28.0 ± 13.1, M249K: − 9.4 ± 11.3; monomer: WT − 83.8 ± 8.2; M249K − 61.5 ± 12.5), suggesting impaired interaction and thus possibly impaired downstream signaling of the complex.

Associations of M249K with cardiovascular risk and replication in external cohorts

Carriers of the M249K variant were at lower risk of myocardial infarction (OR: 0.62, 95%CI: 0.40–0.97, p = 0.03) and coronary artery disease (OR: 0.73, 95%CI: 0.56–0.97, p = 0.02) in the UKB (Fig. 4A). As expected, given that M249K was by far the most frequent among computationally predicted damaging CCR2 variants, we found associations with a lower risk of the combined atherosclerosis endpoint, as well as significantly lower risk of severe atherosclerosis (manifestations in ≥ 2 vascular beds, Fig. 4B), similar to the burden test. To explore whether the effects of M249K are mediated through effects on risk factors targeted by current preventive approaches, we next tested associations with established biomarkers of vascular risk. We found no associations with any of the tested vascular risk factors including blood pressure, hyperglycemia, or circulating lipids. Furthermore, there was no evidence of association with the levels of the known inflammatory biomarker C-reactive protein, thus suggesting that the effects might be independent of factors targeted by current atheroprotective anti-inflammatory treatments (Fig. 4C).

Fig. 4
figure 4

Association of the damaging M249K CCR2 variant with lifetime cardiovascular risk. Associations of the M249K CCR2 variant with risk of A cardiovascular disease endpoints, B a combined atherosclerosis endpoint and severe clinical atherosclerotic disease defined by manifestations in at least 2 vascular beds (coronary arteries, cerebrovascular system, peripheral arteries of extremities, aorta), and C conventional vascular risk biomarkers in models adjusted for age, sex, and the first five principal components among UK Biobank participants of European ancestry. HbA1c and C-reactive protein levels are log-transformed for normalization. Fixed-effects and random-effects meta-analysis of the effects of M249K on risk of D myocardial infarction and E coronary artery disease across 7 population- and hospital-based biobanks. The effects correspond to odds ratios (OR) derived from logistic regression models adjusted for age, sex, and the first 5 ancestral principal components in each biobank

As these associations were only nominally significant and suggestive of an effect of this variant, we aimed to externally replicate our findings to provide robust evidence for an association with a lower risk of myocardial infarction and coronary artery disease. To externally replicate the effects of M249K, we meta-analyzed data from six external datasets (TOPMed, deCODE, PMBB, MGB, MVP, Geisinger). The frequency of the variant varied between the six datasets (range: 0.0006% in MVP to 0.228% in Geisinger), totaling 767 heterozygous carriers among 1,062,595 individuals (91,420 myocardial infarction and 142,361 coronary artery disease cases). When meta-analyzing data from the six cohorts for the association between carrying M249K and risk of myocardial infarction or coronary artery disease, the ORs were comparable to those of the UKB (OR for myocardial infarction: 0.67, 95%CI: 0.53–0.85, p = 0.0006; OR for coronary artery disease: 0.74, 95%CI: 0.53–0.85, p = 0.003, Additional file 1: Fig. S3 and Additional file 2: Table S4). When meta-analyzing these data with UKB reaching a sample size of 1,456,011 individuals including 1314 M249K carriers (111,394 and 207,019 cases of myocardial infarction and coronary artery disease, respectively), we found M249K carriers to have 34% lower odds of suffering a myocardial infarction (OR: 0.66, 95%CI: 0.54–0.81, p = 6.1 × 10−5) and 26% lower odds of coronary artery disease (OR: 0.74, 95%CI: 0.63–0.87, p = 2.9 × 10−4, Fig. 4D–E and Additional file 2: Table S4).

Phenome-wide association study

As a last step, we performed a PheWAS to explore whether M249K is associated with potential adverse effects that would raise signals for possible side-effects of any CCR2-targeting treatments. Because we lacked statistical power for most of the outcomes, we restricted our analyses to endpoints with ≥ 10 cases in the carrier group of the UKB. To increase statistical power, we meta-analyzed our results with the Geisinger cohort. Across all phenotypes, only myocardial infarction showed an association with M249K that passed the FDR-corrected p < 0.05 threshold. Coronary atherosclerosis and angina pectoris were also less frequent (p < 0.05) among carriers although only at a nominal non-corrected p < 0.05 (Fig. 5A). No infectious disease phenotypes were found to be more common (p < 0.05) among M249K carriers (Additional file 2: Table S5). Senile cataract was more common among carriers (p < 0.05). Among UKB participants, we found no evidence for an association between carrier status for a CCR2 damaging variants and survival over a 15-year follow-up period (Fig. 5B).

Fig. 5
figure 5

Phenome-wide association study of the damaging M249K CCR2 variant and associations with overall survival. A Results from a phenome-wide association study in UK Biobank and the Geisinger DiscovEHR-MyCode cohort for the M249K CCR2 variant among participants of European ancestry. Only phecodes with ≥ 10 cases in UK Biobank were analyzed. We present the names of phenotypes associated with rare CCR2 variants at p < 0.05. Only myocardial infarction reached an FDR-corrected p value of < 0.05. The results are presented as log-odds ratios (log-OR). B Kaplan–Meier curves for overall survival across 15 years of follow-up among carriers and non-carriers of the M249K CCR2 variant in the European subset of UK Biobank participants

Discussion

Using data from 1.5M individuals, we found that heterozygous carriers of a rare experimentally confirmed damaging CCR2 variant (M249K) are at lower lifetime risk of myocardial infarction and coronary artery disease. Carriers of this variant showed no differences in LDL cholesterol, blood pressure, BMI, HbA1c, and C-reactive protein levels, indicating that damaging CCR2 variants exert their effects independently of conventional vascular risk factors targeted by available atheroprotective treatments. On the other hand, we found no evidence of associations with higher risk of infectious diseases or overall mortality among carriers. Collectively, our results provide genetic support for CCR2 as a potential therapeutic target for atherosclerotic cardiovascular disease.

Our findings extend our previous results from a Mendelian randomization study suggesting a higher atherosclerosis risk among individuals with genetic predisposition to elevated CCL2 levels [21]. That study focused on common genetic variants distributed across the genome and shown to influence CCL2 levels in trans. In contrast, the current study assessed damaging rare genetic variants within the CCR2 gene. The functional impact of M249K, the most frequent variant among those associated with atherosclerotic disease, was experimentally confirmed. As such, our results are unlikely to be confounded by pleiotropic effects of the variants. Despite the lack of data from clinical trials targeting CCR2 in patients with cardiovascular disease, studies in experimental models of atherosclerosis have provided supportive evidence for an atheroprotective effect of pharmacologically targeting CCR2 [13].

Our work consistently connects these findings to humans and provides a rationale for conducting a clinical trial with pharmacological targeting of CCR2 in patients with atherosclerosis. Ultimately, only clinical testing can provide evidence in favor of or against pharmacological targeting of CCR2 for patients with atherosclerosis. A major bottleneck for initiating a clinical trial would be the need for proof of concept in a phase 2 trial before a trial with hard endpoint is conducted. Given the evidence that CRP levels are not directly affected by the CCL2/CCR2 axis, [19] an intermediate endpoint is needed that captures plaque inflammation, and in particular, macrophage accumulation [13]. PET imaging, particularly with an FDG tracer, has shown promise for detecting macrophage content in large plaques and has already been used in phase 2 trials targeting plaque inflammation [52,53,54]. In addition, recent advances in tracer development have enabled the development of a tracer that specifically binds to CCR2 [55]. Initial data from patients have shown promise in imaging monocyte recruitment to the heart following myocardial infarction [55].

Interestingly, we found no evidence supporting a connection of M249K with established pharmacological targets of atherosclerosis, such as LDL cholesterol, hyperglycemia, or elevated blood pressure. This implies that any benefit from CCR2-targeting approaches would be expected to be independent of available preventive approaches against atherosclerotic disease. Also, we found no evidence of associations with circulating C-reactive protein levels, a biomarker of the activity of the IL-6 signaling pathway that is also the target of atheroprotective anti-inflammatory treatments currently under development [56]. Based on this, our findings suggest that CCR2 inhibition might have an atheroprotective effect on top of such approaches, which in turn might be an argument to pursue approaches embracing combination therapies.

Our PheWAS analysis was limited by low statistical power preventing us from providing sufficient evidence for associations between M249K and a higher risk of key disease endpoints. Future studies should aim to assess potential safety concerns related to CCR2 inhibition by exploring M249K effects on key safety outcomes with in larger meta-analyses of human biobanks. Nevertheless, it is noteworthy that we found no evidence for an association of M249K with risk of infections, which is a potential barrier to the use of atheroprotective anti-inflammatory treatments. This is in line with previous early phase clinical trials testing inhibitors of CCR2 for other indications, which also found no important safety concerns [57,58,59]. Furthermore, we found no evidence of associations with overall survival, minimizing concerns about a significant impact on unknown fatal adverse effects. The reason why we do not see a protective effect on survival despite the effect on risk of myocardial infarction could be related to statistical power and length of follow-up. It is a common observation that medications with strong effects on lowering cardiovascular risk, such as statins, have failed to show survival benefits.

Our study has limitations. First, despite leveraging the largest available whole-exome sequencing studies totaling almost 1.5 million individuals, our results for the PheWAS analyses are still limited by low statistical power and should be interpreted cautiously. Given the rarity of the damaging CCR2 variants, any association with rare potential side-effects would ultimately be undetectable in the context of this study. While formal statistical replication is precluded by this power issue, the consistent effect directions and the overall significance of the meta-analysis raise confidence that the identified association is valid. Further, our main results for associations between CCR2 variants and cardiovascular disease endpoints failed to reach the formal threshold of statistical significance for gene-based testing (p < 2.7 × 10−6). As such, the gene-based results should be interpreted as suggestive. While the associations for M249K also do not reach significance for variant-based testing (p < 5 × 10−8), the replication across studies provides stronger evidence for a real causal association. Genetic analyses with rare variants are still underpowered, as also highlighted in recent burden test studies from the UK Biobank, where no gene showed significant associations with either myocardial infarction or ischemic stroke [60]. Future research should aim to pool additional data as they become available. Second, UK Biobank consists primarily of European individuals, and consequently, we detected damaging CCR2 variants that are predominantly detected in European populations. As such, our results should not be extended to non-European individuals. Third, all analyses are based on individuals heterozygous for the damaging CCR2 variants, as we found no homozygotes for damaging CCR2 variants. While this lack of homozygotes can be fully explained by the frequency of damaging CCR2 variants in the general population, it remains unknown whether homozygous status for damaging CCR2 variants would lead to potentially fatal complications. Fourth, the data from the UKB do not allow a distinction to be made between classical and non-classical monocytes. The available data suggest that CCL2 only acts on classical CCR2 + monocytes in humans or Ly6Chi monocytes in mice. Experimental studies in mice have revealed that Ccr2-deficient Ly6Chi monocytes are trapped in the bone marrow and are not mobilized after infections [35]. Future studies in humans should further explore the impact of the examined mutations on monocyte subtypes. However, it should be noted that inhibition of monocyte recruitment from the circulation into atherosclerotic lesions is the main suspected mechanism by which CCR2 damaging variants act, and the effects of lowering monocyte counts remain to be tested.

Conclusions

In conclusion, we found that heterozygous carriers of a rare damaging CCR2 variant are at lower lifetime risk of myocardial infarction and coronary artery disease without carrying a higher risk of infections. This provides further genetic support for the concept that pharmacological targeting of CCR2 might be an efficient and viable immunotherapeutic strategy to prevent atherosclerotic cardiovascular disease.

Data availability

Individual-level data from the UK Biobank are available to researchers upon approval of a submitted research proposal through https://ams.ukbiobank.ac.uk/ams/. We accessed UK Biobank data following approval of an application by the UKB Ethics and Governance Council (Application No. 36993 and 2532). Summary results from the external cohorts were obtained following a common pre-defined research protocol to principal investigators of the study and are available in the figures and supplementary tables. The code used for the analysis is available in zenodo (link).

References

  1. Diseases GBD, Injuries C. Global burden of 369 diseases and injuries in 204 countries and territories, 1990–2019: a systematic analysis for the Global Burden of Disease Study 2019. Lancet. 2020;396:1204–22. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/S0140-6736(20)30925-9.

    Article  Google Scholar 

  2. Guedeney P, Claessen BE, Kalkman DN, Aquino M, Sorrentino S, Giustino G, Farhan S, Vogel B, Sartori S, Montalescot G, et al. Residual inflammatory risk in patients with low LDL cholesterol levels undergoing percutaneous coronary intervention. J Am Coll Cardiol. 2019;73:2401–9. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.jacc.2019.01.077.

    Article  PubMed  Google Scholar 

  3. Ridker PM. How common is residual inflammatory risk? Circ Res. 2017;120:617–9. https://doiorg.publicaciones.saludcastillayleon.es/10.1161/CIRCRESAHA.116.310527.

    Article  CAS  PubMed  Google Scholar 

  4. Lutgens E, Atzler D, Doring Y, Duchene J, Steffens S, Weber C. Immunotherapy for cardiovascular disease. Eur Heart J. 2019;40:3937–46. https://doiorg.publicaciones.saludcastillayleon.es/10.1093/eurheartj/ehz283.

    Article  CAS  PubMed  Google Scholar 

  5. Weber C, Noels H. Atherosclerosis: current pathogenesis and therapeutic options. Nat Med. 2011;17:1410–22. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/nm.2538.

    Article  CAS  PubMed  Google Scholar 

  6. Ridker PM, Everett BM, Thuren T, MacFadyen JG, Chang WH, Ballantyne C, Fonseca F, Nicolau J, Koenig W, Anker SD, et al. Antiinflammatory therapy with canakinumab for atherosclerotic disease. N Engl J Med. 2017;377:1119–31. https://doiorg.publicaciones.saludcastillayleon.es/10.1056/NEJMoa1707914.

    Article  CAS  PubMed  Google Scholar 

  7. Tardif JC, Kouz S, Waters DD, Bertrand OF, Diaz R, Maggioni AP, Pinto FJ, Ibrahim R, Gamra H, Kiwan GS, et al. Efficacy and safety of low-dose colchicine after myocardial infarction. N Engl J Med. 2019;381:2497–505. https://doiorg.publicaciones.saludcastillayleon.es/10.1056/NEJMoa1912388.

    Article  CAS  PubMed  Google Scholar 

  8. Nidorf SM, Fiolet ATL, Mosterd A, Eikelboom JW, Schut A, Opstal TSJ, The SHK, Xu XF, Ireland MA, Lenderink T, et al. Colchicine in patients with chronic coronary disease. N Engl J Med. 2020;383:1838–47. https://doiorg.publicaciones.saludcastillayleon.es/10.1056/NEJMoa2021372.

    Article  CAS  PubMed  Google Scholar 

  9. Leung YY, Yao Hui LL, Kraus VB. Colchicine–update on mechanisms of action and therapeutic uses. Semin Arthritis Rheum. 2015;45:341–50. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.semarthrit.2015.06.013.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Opstal TSJ, Hoogeveen RM, Fiolet ATL, Silvis MJM, The SHK, Bax WA, de Kleijn DPV, Mosterd A, Stroes ESG, Cornel JH. Colchicine attenuates inflammation beyond the inflammasome in chronic coronary artery disease: a LoDoCo2 proteomic substudy. Circulation. 2020;142:1996–8. https://doiorg.publicaciones.saludcastillayleon.es/10.1161/CIRCULATIONAHA.120.050560.

    Article  CAS  PubMed  Google Scholar 

  11. Ridker PM. From C-reactive protein to interleukin-6 to interleukin-1: moving upstream to identify novel targets for atheroprotection. Circ Res. 2016;118:145–56. https://doiorg.publicaciones.saludcastillayleon.es/10.1161/CIRCRESAHA.115.306656.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Mackay CR. Chemokines: immunology’s high impact factors. Nat Immunol. 2001;2:95–101.

    Article  CAS  PubMed  Google Scholar 

  13. Zivkovic L, Asare Y, Bernhagen J, Dichgans M, Georgakis MK. Pharmacological targeting of the CCL2/CCR2 axis for atheroprotection: a meta-analysis of preclinical studies. Arterioscler Thromb Vasc Biol. 2022;42:e131–44. https://doiorg.publicaciones.saludcastillayleon.es/10.1161/ATVBAHA.122.317492.

    Article  CAS  PubMed  Google Scholar 

  14. Boring L, Gosling J, Cleary M, Charo IF. Decreased lesion formation in CCR2-/- mice reveals a role for chemokines in the initiation of atherosclerosis. Nature. 1998;394:894–7. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/29788.

    Article  CAS  PubMed  Google Scholar 

  15. Dawson TC, Kuziel WA, Osahar TA, Maeda N. Absence of CC chemokine receptor-2 reduces atherosclerosis in apolipoprotein E-deficient mice. Atherosclerosis. 1999;143:205–11. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/s0021-9150(98)00318-9.

    Article  CAS  PubMed  Google Scholar 

  16. Gu L, Okada Y, Clinton SK, Gerard C, Sukhova GK, Libby P, Rollins BJ. Absence of monocyte chemoattractant protein-1 reduces atherosclerosis in low density lipoprotein receptor-deficient mice. Mol Cell. 1998;2:275–81. https://doiorg.publicaciones.saludcastillayleon.es/S1097-2765(00)80139-2 [pii].

    Article  CAS  PubMed  Google Scholar 

  17. Aiello RJ, Bourassa PA, Lindsey S, Weng W, Natoli E, Rollins BJ, Milos PM. Monocyte chemoattractant protein-1 accelerates atherosclerosis in apolipoprotein E-deficient mice. Arterioscler Thromb Vasc Biol. 1999;19:1518–25.

    Article  CAS  PubMed  Google Scholar 

  18. Georgakis MK, Bernhagen J, Heitman LH, Weber C, Dichgans M. Targeting the CCL2-CCR2 axis for atheroprotection. Eur Heart J. 2022;43:1799–808. https://doiorg.publicaciones.saludcastillayleon.es/10.1093/eurheartj/ehac094.

    Article  CAS  PubMed  Google Scholar 

  19. Georgakis MK, Malik R, Bjorkbacka H, Pana TA, Demissie S, Ayers C, Elhadad MA, Fornage M, Beiser AS, Benjamin EJ, et al. Circulating monocyte chemoattractant protein-1 and risk of stroke: meta-analysis of population-based studies involving 17 180 individuals. Circ Res. 2019;125:773–82. https://doiorg.publicaciones.saludcastillayleon.es/10.1161/CIRCRESAHA.119.315380.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Georgakis MK, de Lemos JA, Ayers C, Wang B, Bjorkbacka H, Pana TA, Thorand B, Sun C, Fani L, Malik R, et al. Association of circulating monocyte chemoattractant protein-1 levels with cardiovascular mortality: a meta-analysis of population-based studies. JAMA Cardiol. 2021;6:587–92. https://doiorg.publicaciones.saludcastillayleon.es/10.1001/jamacardio.2020.5392.

    Article  PubMed  Google Scholar 

  21. Georgakis MK, Gill D, Rannikmae K, Traylor M, Anderson CD, Lee JM, Kamatani Y, Hopewell JC, Worrall BB, Bernhagen J, et al. Genetically determined levels of circulating cytokines and risk of stroke. Circulation. 2019;139:256–68. https://doiorg.publicaciones.saludcastillayleon.es/10.1161/CIRCULATIONAHA.118.035905.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Georgakis MK, Parodi L, Frerich S, Mayerhofer E, Tsivgoulis G, Pirruccello JP, Slowik A, Rundek T, Malik R, Dichgans M, et al. Genetic architecture of stroke of undetermined source: overlap with known stroke etiologies and associations with modifiable risk factors. Ann Neurol. 2022;91:640–51. https://doiorg.publicaciones.saludcastillayleon.es/10.1002/ana.26332.

    Article  CAS  PubMed  Google Scholar 

  23. Georgakis MK, van der Laan SW, Asare Y, Mekke JM, Haitjema S, Schoneveld AH, de Jager SCA, Nurmohamed NS, Kroon J, Stroes ESG, et al. Monocyte-chemoattractant protein-1 levels in human atherosclerotic lesions associate with plaque vulnerability. Arterioscler Thromb Vasc Biol. 2021;41:2038–48. https://doiorg.publicaciones.saludcastillayleon.es/10.1161/ATVBAHA.121.316091.

    Article  CAS  PubMed  Google Scholar 

  24. Interleukin-6 Receptor Mendelian Randomisation Analysis C, Swerdlow DI, Holmes MV, Kuchenbaecker KB, Engmann JE, Shah T, Sofat R, Guo Y, Chung C, Peasey A, et al. The interleukin-6 receptor as a target for prevention of coronary heart disease: a Mendelian randomisation analysis. Lancet. 2012;379:1214–1224. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/S0140-6736(12)60110-X.

  25. Cohen JC, Boerwinkle E, Mosley TH Jr, Hobbs HH. Sequence variations in PCSK9, low LDL, and protection against coronary heart disease. N Engl J Med. 2006;354:1264–72. https://doiorg.publicaciones.saludcastillayleon.es/10.1056/NEJMoa054013.

    Article  CAS  PubMed  Google Scholar 

  26. Clarke R, Peden JF, Hopewell JC, Kyriakou T, Goel A, Heath SC, Parish S, Barlera S, Franzosi MG, Rust S, et al. Genetic variants associated with Lp(a) lipoprotein level and coronary disease. N Engl J Med. 2009;361:2518–28. https://doiorg.publicaciones.saludcastillayleon.es/10.1056/NEJMoa0902604.

    Article  CAS  PubMed  Google Scholar 

  27. Nioi P, Sigurdsson A, Thorleifsson G, Helgason H, Agustsdottir AB, Norddahl GL, Helgadottir A, Magnusdottir A, Jonasdottir A, Gretarsdottir S, et al. Variant ASGR1 associated with a reduced risk of coronary artery disease. N Engl J Med. 2016;374:2131–41. https://doiorg.publicaciones.saludcastillayleon.es/10.1056/NEJMoa1508419.

    Article  CAS  PubMed  Google Scholar 

  28. Stitziel NO, Khera AV, Wang X, Bierhals AJ, Vourakis AC, Sperry AE, Natarajan P, Klarin D, Emdin CA, Zekavat SM, et al. ANGPTL3 deficiency and protection against coronary artery disease. J Am Coll Cardiol. 2017;69:2054–63. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.jacc.2017.02.030.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Sudlow C, Gallacher J, Allen N, Beral V, Burton P, Danesh J, Downey P, Elliott P, Green J, Landray M, et al. UK Biobank: an open access resource for identifying the causes of a wide range of complex diseases of middle and old age. PLoS Med. 2015;12: e1001779. https://doiorg.publicaciones.saludcastillayleon.es/10.1371/journal.pmed.1001779.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Backman JD, Li AH, Marcketta A, Sun D, Mbatchou J, Kessler MD, Benner C, Liu D, Locke AE, Balasubramanian S, et al. Exome sequencing and analysis of 454,787 UK Biobank participants. Nature. 2021;599:628–34. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/s41586-021-04103-z.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Malik R, Beaufort N, Frerich S, Gesierich B, Georgakis MK, Rannikmae K, Ferguson AC, Haffner C, Traylor M, Ehrmann M, et al. Whole-exome sequencing reveals a role of HTRA1 and EGFL8 in brain white matter hyperintensities. Brain. 2021;144:2670–82. https://doiorg.publicaciones.saludcastillayleon.es/10.1093/brain/awab253.

    Article  PubMed  PubMed Central  Google Scholar 

  32. McLaren W, Gil L, Hunt SE, Riat HS, Ritchie GR, Thormann A, Flicek P, Cunningham F. The Ensembl variant effect predictor. Genome Biol. 2016;17:122. https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s13059-016-0974-4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Ioannidis NM, Rothstein JH, Pejaver V, Middha S, McDonnell SK, Baheti S, Musolf A, Li Q, Holzinger E, Karyadi D, et al. REVEL: an ensemble method for predicting the pathogenicity of rare missense variants. Am J Hum Genet. 2016;99:877–85. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.ajhg.2016.08.016.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Karczewski KJ, Francioli LC, Tiao G, Cummings BB, Alfoldi J, Wang Q, Collins RL, Laricchia KM, Ganna A, Birnbaum DP, et al. The mutational constraint spectrum quantified from variation in 141,456 humans. Nature. 2020;581:434–43. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/s41586-020-2308-7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Serbina NV, Pamer EG. Monocyte emigration from bone marrow during bacterial infection requires signals mediated by chemokine receptor CCR2. Nat Immunol. 2006;7:311–7. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/ni1309.

    Article  CAS  PubMed  Google Scholar 

  36. Tsou CL, Peters W, Si Y, Slaymaker S, Aslanian AM, Weisberg SP, Mack M, Charo IF. Critical roles for CCR2 and MCP-3 in monocyte mobilization from bone marrow and recruitment to inflammatory sites. J Clin Invest. 2007;117:902–9. https://doiorg.publicaciones.saludcastillayleon.es/10.1172/JCI29919.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Shi C, Jia T, Mendez-Ferrer S, Hohl TM, Serbina NV, Lipuma L, Leiner I, Li MO, Frenette PS, Pamer EG. Bone marrow mesenchymal stem and progenitor cells induce monocyte emigration in response to circulating toll-like receptor ligands. Immunity. 2011;34:590–601. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.immuni.2011.02.016.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Mbatchou J, Barnard L, Backman J, Marcketta A, Kosmicki JA, Ziyatdinov A, Benner C, O’Dushlaine C, Barber M, Boutkov B, et al. Computationally efficient whole-genome regression for quantitative and binary traits. Nat Genet. 2021;53:1097–103. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/s41588-021-00870-7.

    Article  CAS  PubMed  Google Scholar 

  39. Honorato RV, Trellet ME, Jimenez-Garcia B, Schaarschmidt JJ, Giulini M, Reys V, Koukos PI, Rodrigues J, Karaca E, van Zundert GCP, et al. The HADDOCK2.4 web server for integrative modeling of biomolecular complexes. Nat Protoc. 2024;19:3219–3241. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/s41596-024-01011-0.

  40. Taliun D, Harris DN, Kessler MD, Carlson J, Szpiech ZA, Torres R, Taliun SAG, Corvelo A, Gogarten SM, Kang HM, et al. Sequencing of 53,831 diverse genomes from the NHLBI TOPMed Program. Nature. 2021;590:290–9. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/s41586-021-03205-y.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Aegisdottir HM, Thorolfsdottir RB, Sveinbjornsson G, Stefansson OA, Gunnarsson B, Tragante V, Thorleifsson G, Stefansdottir L, Thorgeirsson TE, Ferkingstad E, et al. Genetic variants associated with syncope implicate neural and autonomic processes. Eur Heart J. 2023;44:1070–80. https://doiorg.publicaciones.saludcastillayleon.es/10.1093/eurheartj/ehad016.

    Article  CAS  PubMed  Google Scholar 

  42. Kember RL, Merikangas AK, Verma SS, Verma A, Judy R, Regeneron Genetics C, Damrauer SM, Ritchie MD, Rader DJ, Bucan M. Polygenic risk of psychiatric disorders exhibits cross-trait associations in electronic health record data from European ancestry individuals. Biol Psychiatry. 2021;89:236–45. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.biopsych.2020.06.026.

    Article  CAS  PubMed  Google Scholar 

  43. Verma A, Huffman JE, Rodriguez A, Conery M, Liu M, Ho Y-L, Kim Y, Heise DA, Guare L, Panickan VA, et al. Diversity and scale: genetic architecture of 2,068 traits in the VA Million Veteran Program. medRxiv. 2023:2023.2006.2028.23291975. https://doiorg.publicaciones.saludcastillayleon.es/10.1101/2023.06.28.23291975.

  44. Kelly MA, Leader JB, Wain KE, Bodian D, Oetjens MT, Ledbetter DH, Martin CL, Strande NT. Leveraging population-based exome screening to impact clinical care: the evolution of variant assessment in the Geisinger MyCode research project. Am J Med Genet C Semin Med Genet. 2021;187:83–94. https://doiorg.publicaciones.saludcastillayleon.es/10.1002/ajmg.c.31887.

    Article  PubMed  Google Scholar 

  45. Koyama S, Wang Y, Paruchuri K, Uddin MM, Cho SMJ, Urbut SM, Haidermota S, Hornsby WE, Green RC, Daly MJ, et al. Decoding genetics, ancestry, and geospatial context for precision health. medRxiv. 2023:2023.2010.2024.23297096. https://doiorg.publicaciones.saludcastillayleon.es/10.1101/2023.10.24.23297096.

  46. Wu P, Gifford A, Meng X, Li X, Campbell H, Varley T, Zhao J, Carroll R, Bastarache L, Denny JC, et al. Mapping ICD-10 and ICD-10-CM codes to phecodes: workflow development and initial evaluation. JMIR Med Inform. 2019;7: e14325. https://doiorg.publicaciones.saludcastillayleon.es/10.2196/14325.

    Article  PubMed  PubMed Central  Google Scholar 

  47. Grundmann M, Merten N, Malfacini D, Inoue A, Preis P, Simon K, Ruttiger N, Ziegler N, Benkel T, Schmitt NK, et al. Lack of beta-arrestin signaling in the absence of active G proteins. Nat Commun. 2018;9:341. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/s41467-017-02661-3.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Rajagopal S, Rajagopal K, Lefkowitz RJ. Teaching old receptors new tricks: biasing seven-transmembrane receptors. Nat Rev Drug Discov. 2010;9:373–86. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/nrd3024.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Islam LN, Wilkinson PC. Chemotactic factor-induced polarization, receptor redistribution, and locomotion of human blood monocytes. Immunology. 1988;64:501–7.

    CAS  PubMed  PubMed Central  Google Scholar 

  50. Schroder JM, Persoon NL, Christophers E. Lipopolysaccharide-stimulated human monocytes secrete, apart from neutrophil-activating peptide 1/interleukin 8, a second neutrophil-activating protein. NH2-terminal amino acid sequence identity with melanoma growth stimulatory activity. J Exp Med. 1990;171:1091–100. https://doiorg.publicaciones.saludcastillayleon.es/10.1084/jem.171.4.1091.

    Article  CAS  PubMed  Google Scholar 

  51. Uguccioni M, D’Apuzzo M, Loetscher M, Dewald B, Baggiolini M. Actions of the chemotactic cytokines MCP-1, MCP-2, MCP-3, RANTES, MIP-1 alpha and MIP-1 beta on human monocytes. Eur J Immunol. 1995;25:64–8. https://doiorg.publicaciones.saludcastillayleon.es/10.1002/eji.1830250113.

    Article  CAS  PubMed  Google Scholar 

  52. Fayad ZA, Mani V, Woodward M, Kallend D, Abt M, Burgess T, Fuster V, Ballantyne CM, Stein EA, Tardif JC, et al. Safety and efficacy of dalcetrapib on atherosclerotic disease using novel non-invasive multimodality imaging (dal-PLAQUE): a randomised clinical trial. Lancet. 2011;378:1547–59. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/s0140-6736(11)61383-4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Gaztanaga J, Farkouh M, Rudd JH, Brotz TM, Rosenbaum D, Mani V, Kerwin TC, Taub R, Tardif JC, Tawakol A, et al. A phase 2 randomized, double-blind, placebo-controlled study of the effect of VIA-2291, a 5-lipoxygenase inhibitor, on vascular inflammation in patients after an acute coronary syndrome. Atherosclerosis. 2015;240:53–60. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.atherosclerosis.2015.02.027.

    Article  CAS  PubMed  Google Scholar 

  54. van Wijk DF, Sjouke B, Figueroa A, Emami H, van der Valk FM, MacNabb MH, Hemphill LC, Schulte DM, Koopman MG, Lobatto ME, et al. Nonpharmacological lipoprotein apheresis reduces arterial inflammation in familial hypercholesterolemia. J Am Coll Cardiol. 2014;64:1418–26. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.jacc.2014.01.088.

    Article  CAS  PubMed  Google Scholar 

  55. Liu Y, Imanishi T, Ikejima H, Tsujioka H, Ozaki Y, Kuroi A, Okochi K, Ishibashi K, Tanimoto T, Ino Y, et al. Association between circulating monocyte subsets and in-stent restenosis after coronary stent implantation in patients with ST-elevation myocardial infarction. Circ J. 2010;74:2585–91. https://doiorg.publicaciones.saludcastillayleon.es/10.1253/circj.cj-10-0544.

    Article  PubMed  Google Scholar 

  56. Ridker PM, Devalaraja M, Baeres FMM, Engelmann MDM, Hovingh GK, Ivkovic M, Lo L, Kling D, Pergola P, Raj D, et al. IL-6 inhibition with ziltivekimab in patients at high atherosclerotic risk (RESCUE): a double-blind, randomised, placebo-controlled, phase 2 trial. Lancet. 2021;397:2060–9. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/S0140-6736(21)00520-1.

    Article  CAS  PubMed  Google Scholar 

  57. Friedman SL, Ratziu V, Harrison SA, Abdelmalek MF, Aithal GP, Caballeria J, Francque S, Farrell G, Kowdley KV, Craxi A, et al. A randomized, placebo-controlled trial of cenicriviroc for treatment of nonalcoholic steatohepatitis with fibrosis. Hepatology. 2018;67:1754–67. https://doiorg.publicaciones.saludcastillayleon.es/10.1002/hep.29477.

    Article  CAS  PubMed  Google Scholar 

  58. de Zeeuw D, Bekker P, Henkel E, Hasslacher C, Gouni-Berthold I, Mehling H, Potarca A, Tesar V, Heerspink HJ, Schall TJ, et al. The effect of CCR2 inhibitor CCX140-B on residual albuminuria in patients with type 2 diabetes and nephropathy: a randomised trial. Lancet Diabetes Endocrinol. 2015;3:687–96. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/S2213-8587(15)00261-2.

    Article  CAS  PubMed  Google Scholar 

  59. Vergunst CE, Gerlag DM, Lopatinskaya L, Klareskog L, Smith MD, van den Bosch F, Dinant HJ, Lee Y, Wyant T, Jacobson EW, et al. Modulation of CCR2 in rheumatoid arthritis: a double-blind, randomized, placebo-controlled clinical trial. Arthritis Rheum. 2008;58:1931–9. https://doiorg.publicaciones.saludcastillayleon.es/10.1002/art.23591.

    Article  CAS  PubMed  Google Scholar 

  60. Karczewski KJ, Solomonson M, Chao KR, Goodrich JK, Tiao G, Lu W, Riley-Gillis BM, Tsai EA, Kim HI, Zheng X, et al. Systematic single-variant and gene-based association testing of thousands of phenotypes in 394,841 UK Biobank exomes. Cell Genom. 2022;2:100168. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.xgen.2022.100168.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

We thank Markus Brandhofer (Institute for Stroke and Dementia Research) for assistance and helpful discussions regarding the HADDOCK docking experiment. We acknowledge the Penn Medicine BioBank (PMBB) for providing data and thank the patient-participants of Penn Medicine who consented to participate in this research program. We would also like to thank the Penn Medicine BioBank team and Regeneron Genetics Center for providing genetic variant data for analysis. We thank Geisinger MyCode Community Health Initiative and Geisinger-Regeneron DiscovEHR collaboration contributors who have been critical in the generation of the original data used for this study. Cohort-specific acknowledgements for studies included in the TOPMed program are provided in Additional file 2: Table S6.

Funding

Open Access funding enabled and organized by Projekt DEAL. Marios K. Georgakis is supported by the German Research Foundation (Deutsche Forschungsgemeinschaft, DFG) with an Emmy Noether grant (GZ: GE 3461/2–1, ID 512461526) and a Munich Cluster for Systems Neurology Early Excellence Academy fellowship (EXC 2145 SyNergy, ID 390857198) as well as by the Fritz-Thyssen Foundation (Ref. 10.22.2.024MN), the Hertie Foundation (Hertie Network of Excellence in Clinical Neuroscience, ID P1230035), and LMU Munich (LMU-TAU Research Cooperation Program). Scott M. Damrauer is supported by IK2-CX001780. The PMBB is supported by Perelman School of Medicine at University of Pennsylvania, a gift from the Smilow family, and the National Center for Advancing Translational Sciences of the National Institutes of Health under CTSA award number UL1TR001878. This publication does not represent the views of the Department of Veterans Affairs or the US Government. Paul S. de Vries, Natalie R. Hasbani, and Alanna C. Morrison were supported by R01 HL146860. Joshua C. Bis and Bruce M. Psaty were funded in part by R01HL105756. Reading of the carotid IMT measures was supported by N01HC85085 and N01HC45133; Coronary Calcium scans were supported by R01HL64587. Kendra A. Young is supported by NHLBIn U01 HL089897 and U01 HL089856. Sharon M. Lutz was supported by R01 MH129337. Nicholette D. Palmer, Barry I. Freedman, and Donald W. Bowden were supported by R01 AG058921. Patricia A. Peyser was supported by R01 HL146860. Lawrence F. Bielak was supported by R01 HL146860. LMR was supported by T32 HL129982 and by the National Center for Advancing Translational Sciences, National Institutes of Health, through Grant KL2TR002490. Christopher D. Anderson receives sponsored research support from the US National Institutes of Health (R01NS103924, U01NS069673) and the American Heart Association (AHA 21SFRN812095). Whole genome sequencing (WGS) for the Trans-Omics in Precision Medicine (TOPMed) program was supported by the National Heart, Lung and Blood Institute (NHLBI). Jürgen Bernhagen is supported by grants from the Deutsche Forschungsgemeinschaft (DFG; CRC 1123 [A3], BE 1977/14–1, and Munich Cluster for Systems Neurology [SyNergy, EXC 2145; ID390857198]). Martin Dichgans is supported by grants from the Deutsche Forschungsgemeinschaft (DFG; CRC 1123 [B3], DI-722/16–1 [ID: 428668490], DI-722/21–1, and Munich Cluster for Systems Neurology [SyNergy, EXC 2145; ID390857198]), the Leducq Foundation, the flagship P4-medicine project DigiMed Bayern, and the Vascular Dementia Research Foundation.

Author information

Authors and Affiliations

Authors

Contributions

MKG conceived the study, analyzed the data of the UK Biobank, performed the meta-analyses for the M249K variant, supervised this work, drafted the first version of the manuscript and reviewed the manuscript for intellectual content. RM performed the burden test analyses and phenome-wide association study in the UK Biobank, drafted the first version of the manuscript and reviewed the manuscript for intellectual content. OEB performed the variant transfection experiments and the chemotaxis assays and reviewed the manuscript for intellectual content. NRB analyzed the TOPMed data and reviewed the manuscript for intellectual content. JL analyzed the original Geisinger data and reviewed the manuscript for intellectual content. JEH analyzed the original MVP data and reviewed the manuscript for intellectual content. GS analyzed the original PMBB data and reviewed the manuscript for intellectual content. RWPT and TNK analyzed the original MGB data and reviewed the manuscript for intellectual content. YA assisted in the chemotaxis assays and reviewed the manuscript for intellectual content. ACM, MLT, RJ, BDM, HX, MEM, RD, EEK, RJFL, JGT, JJC, JCB, BMP, WTL, LAY, SML, MHC, JB, ATK, FFW, NHC, SS, RSV, NDP, BIF, DWB, LRY, BGK, LCB, PAP, LFB, FA, APC, MEH, LMR, SSR, WSP, RPT, KDT, XG, MCM, JEC, JB, SLC, JWH, YH, TLA, CK contributed data to the TOPMed analyses and reviewed the manuscript for intellectual content. JB designed and supervised the variant transfection experiments and the chemotaxis assays and reviewed the manuscript for intellectual content. CDA designed and supervised the analyses in MGB and reviewed the manuscript for intellectual content. SMD designed and supervised the analyses in PMBB and reviewed the manuscript for intellectual content. RZ designed and supervised the analyses in Geisinger and reviewed the manuscript for intellectual content. JIR and PSdV designed and supervised the analyses in TOPMed and reviewed the manuscript for intellectual content. MD conceived and supervised this work, drafted the first version of the manuscript and reviewed the manuscript for intellectual content. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Marios K. Georgakis or Martin Dichgans.

Ethics declarations

Ethics approval and consent to participate

All studies have received ethical approvals from the respective ethical authorities and participants in all studies have provided written informed consent. Data from the UKB are available upon submission and approval of a research proposal. The UKB has institutional review board approval from the Northwest Multi‐Center Research Ethics Committee (Manchester, UK). deCODE has been approved by the National Bioethics Committee of Iceland, and the study was conducted in agreement with conditions issued by the Data Protection Authority of Iceland (VSN_14-015). All common research protocols for the TOPMed Program have been approved by the institutional review board at the University of Maryland Baltimore, whereas individual participating studies have obtained ethical approval from their local ethical authorities, as described previously [60]. The Penn Medicine BioBank is approved by the University of Pennsylvania under IRB protocol# 813913, MVP by the Veterans Affairs Central Institutional Review Board, the Geisinger DiscovEHR-MyCode cohort by Geisinger Institutional Review Boards, and MGBB by the Ethics Committee of Mass General Brigham.

Consent for publication

Not applicable.

Competing interests

MKG has received consulting fees from Tourmaline Bio Inc., and Gerson Lehrman Group Inc. (unrelated to this work) and is a member of the Editorial Board of Neurology. SMD receives research support from RenalytixAI and personal fees from Calico Labs, both outside the current work. MEM receives funding from Regeneron Pharmaceutical Inc. unrelated to this work. RD has received research support from AstraZeneca and Goldfinch Bio, not related to this work. BMP serves on the Steering Committee of the Yale Open Data Access Program funded by Johnson & Johnson. CDA has received sponsored research support from Bayer AG, has consulted for ApoPharma, and is a member of the Editorial Board of Neurology. LMR is a consultant for the TOPMed Administrative Coordinating Center (through WeStat). The remaining authors have no conflicts of interest to disclose.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Georgakis, M.K., Malik, R., Bounkari, O.E. et al. Rare damaging CCR2 variants are associated with lower lifetime cardiovascular risk. Genome Med 17, 27 (2025). https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s13073-025-01456-2

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s13073-025-01456-2

Keywords