- 1Servicio de Cardiología, Complejo Hospitalario Universitario Insular Materno-Infantil, Las Palmas de Gran Canaria, Spain
- 2IMDEA Research Institute of Food and Health Sciences, Madrid, Spain
- 3Departamento de Matemáticas, Universidad de Las Palmas de Gran Canaria, Las Palmas de Gran Canaria, Spain
- 4SECUGEN SL, Madrid, Spain
- 5Servicio de Oftalmología, Hospital Universitario Gran Canaria Doctor Negrín, Las Palmas de Gran Canaria, Spain
- 6Servicio de Bioquímica Clínica y Análisis Clínicos, Complejo Hospitalario Universitario Insular Materno-Infantil, Las Palmas de Gran Canaria, Spain
- 7Unidad de Investigación, Complejo Hospitalario Universitario Insular Materno-Infantil, Las Palmas de Gran Canaria, Spain
Monocytes participate in the development of atherosclerosis through the action of cytokines and other inflammatory mediators. Among them, CCR2 and its ligands, CCL2 and CCL7 play an important role, so the main objective of this work was to determine whether genetic variants affecting their activity were associated with cardiovascular disease. A cohort of 519 patients that have suffered coronary events was analyzed under a propensity score-matching protocol selecting a homogeneous set of cases and controls, according to age, sex, smoking status, dyslipidemia, arterial hypertension and type 2 diabetes as risk factors. While dyslipidemia and arterial hypertension were more prevalent among patients with angina pectoris, current smoking status and elevated inflammatory markers, including total leukocyte and monocyte counts, were more likely associated with acute coronary events. Propensity score matching analysis, performed to eliminate the influence of these risk factors and highlight genetic modifiers, revealed that a single nucleotide variant, rs17735770 at the 3'untranslated region of the CCL7 gene transcript, was associated with decreased cardiovascular risk in a group represented mostly by men, with an average age of 57, and without significant differences in traditional risk factors. Furthermore, the presence of this variant altered the local mRNA structure encompassing a binding site for miR-23ab, resulting in increased translation of a reporter gene in a miR23 independent fashion. The rs17735770 genetic variant led to increased expression of CCL7, a potential antagonist of CCR2 at inflammatory sites, where it could play a meaningful role during the evolution of atherosclerosis.
Introduction
Beyond classical cardiovascular risk factors, inflammation appears as a major player in the formation, evolution and maturation of atheroma plaques in the endothelium (1, 2). Monocytes play a major role from start in this process, adhering to the endothelium and infiltrating into inflammatory sites in response to chemoatractive chemokines, becoming macrophages and, later, foam cells which are the origin of the intimal fat strike (3). After death, foam cells contribute to the formation of the necrotic core of the plaque and, by releasing metalloproteinases and cathepsins, destabilize the plaque by thinning the fibrous outer capsular layer which, upon breakage and tissue factor release, favors coronary thrombosis and the occurrence of an acute coronary event (1–3).
Chemokines and their receptors are important mediators of inflammation and tissue morphogenesis, including hematopoiesis and angiogenesis, through the induction of chemotaxis, proliferation and apoptosis of many cell types (4). In particular the C-C type chemokine receptor type 2 (CCR2) and its ligand CCL2, and later CX3CR1 and CCR5 have been associated with the recruitment of monocytes to inflammatory sites (5–8).
Mice deficient in the CCR2 gene show defects in monocyte migration, reduced leukocyte adhesion, and monocyte recruitment at the plaques (9–11) and rescue the development of atherosclerotic lesions under a high fat diet, when crossed to hyperlipidemic mice, deficient in ApoE (12, 13). Accordingly, mice deficient in CCL2 are more resistant to the development of atherosclerosis in response to a high cholesterol diet than their wild type littermates (14, 15).
While CCL2 is an exclusive CCR2 ligand, CCL7 may also bind to other chemokine receptors present in a variety of cell types (4). At inflammatory sites, where there is a presence of metalloproteinases such as gelatinases (MMP-2 y MMP-9) and MMP-7, the CCL7 molecule may be proteolyzed to form a competitive CCR2 antagonist that may modulate its response (16).
Given the importance of the CCR2 receptor in the development of atherosclerosis in mice, some efforts have been directed toward the identification of genetic variants that may modulate its activity in humans. The frequent polymorphic variant p.Val64Ile (SNP Id: rs1799864) has been associated with increased risk of myocardial infarction and heart failure (17), but not always (18–21), and is also associated with decreased coronary calcification (22). Likewise, the association of variations in the ligand CCL2 with cardiovascular disease is also controversial (23–25). Therefore, the possibility that genetic variants in this receptor-ligand system could influence the development of cardiovascular disease in our patients was explored.
Materials and Methods
Subjects
All participants gave written informed consent before being included in the study, which was also approved by the Ethics Committee of our institution. The study included 519 patients, between 25 and 90 years old, that presented with a coronary event and required hospital admission, and a set of 760 controls (26). Excluded from the study were patients without a clear coronary event diagnosis, those with cancer, fever, or any known inflammatory disease, and those that did not wish to participate. In-hospital data at admission included: anthropomorphic and demographic data, analytical data and the type of acute coronary event: stable angina, unstable angina, non-ST (NSTEMI) and ST-elevation myocardial infarction (STEMI) (27).
Genetic Analysis
Genomic DNA was extracted from white blood cells with the Puregene kit (QUIAGEN, Hilden, Germany). An initial screen for the detection of genetic variants in the CCR2, CCL2 and CCL7 genes was performed in 70 cases and 30 healthy controls by heterozygous double pass Sanger sequencing of amplified DNA fragments containing coding exons and flanking sequences. All primers used for polymerase chain reaction (PCR) are shown in a Supplementary File.
Sequences were analyzed with the 4Peaks software (Nucleobytes Inc., The Netherlands), and aligned against reference sequences using BLAST (28) to identify variants, that were later searched in the dbSNP database (29). Genotyping of the rs17735770 variant involved the design of an RFLP (restriction fragment length polymorphism) assay, as the wild type to variant transition generates a novel SnaBI restriction site. Primers used for PCR (see Supplementary Methods for sequences) were designed with the PrimerBLAST tool (30). Genotyping was performed on 771 controls and 519 cases by using PCR followed by fragment electrophoretic analysis essentially as described (31).
Bioinformatic Analysis
Species-specific searches to detect alignment of the CCL7 gene 3' UTR region in mammals was performed with BLAST (28). In order to detect predicted miRNA binding motifs in the region containing the rs17735770 variant, the TargetScanHuman® v. 5.1 application (32) was used as implemented on the UCSC Genome Browser (33). The algorithms RNAfold and RNAcofold were used on the Vienna RNA websuite (34).
3'UTR Luciferase Reporter Assays
cDNA fragments corresponding to the entire 3'UTR of the human CCL7 gene were amplified from human genomic DNA by PCR with Phusion Hot Start II DNA Polymerase (Thermo Fisher Scientific, Waltham, MA) using primers flanked by XhoI (5' forward primer) and NotI (3' reverse primer) overhangs (sequences in Supplementary Methods). The resulting PCR products were directionally cloned into the XhoI and NotI sites of the psiCHECK-2 vector (Promega Biotech, Madison, WI), downstream of the Renilla luciferase (hRluc) ORF. The vector also contains a constitutively expressed humanized Photynus luciferase (hluc) gene, which is used to normalize transfections. The rs17735770 variant and a disrupting point mutation in the seed region of the predicted miR-23a-3p site (AATGTGAA to AATCAGGA) within the CCL7 3'UTR were generated by using QuikChange Multi Site Directed-Mutagenesis Kit (Agilent Technologies Inc, Wood Dale, IL) according to the manufacturer's protocol (primers sequences in Supplementary Methods). All constructs were checked by Sanger sequencing of individual clones.
HEK and COS-7 cells were plated into 12-well plates and co-transfected with 1 μg of the indicated 3'UTR luciferase reporter vectors and 40nM of the miR-23a-3p mimic or nontargeting control mimic (CM) sequence (Horizon Discovery, formerly Dharmacon, Lafayette, CO) utilizing Lipofectamine 2000 (Invitrogen-Thermo Fisher Scientific). In the case of THP-1, cells were treated with 100 nM phorbol myristate acetate (61 mg/ml PMA) for 24 h before transfection. Luciferase activity was measured 24 h after transfection using the Dual-Glo Luciferase Assay System (Promega). Renilla luciferase activity was normalized to the corresponding Photynus luciferase activity and plotted as a percentage of the control (cells co-transfected with the corresponding concentration of control mimic). At least three independent experiments were performed in triplicates.
Statistical Analysis
Categorical variables were expressed as frequencies and percentages and continuous as means and standard deviations (SD). The percentages were compared, as appropriate, using the Chi-square (χ2) test or the exact Fisher test and the means by the t-test. Multiple comparisons were carried out, as appropriate, using Scheffe's method or a non-parametric approach (35).
In order to determine the association between the presence of coronary heart disease (CHD) and the rs17735770 variant (T → C) at the CCL7 gene, a similar control was selected for each case (matching), based on a propensity score, defined as the conditional probability of suffering CHD given a set of covariates obtained for each individual using logistic regression, in which the outcome variable was CHD. The covariates included in the model were selected using the best subset regression procedure together with and Akaike Information Criterion (AIC) (36).
A 1-to-1 matched analysis without replacement on the basis of the estimated propensity score of each patient was performed. The caliper chosen was 0.3. After propensity score matching, baseline characteristics were compared with the McNemar tests for binary variables and the t-tests or Wilcoxon test, as appropriate, for continuous variables, both for paired data. In addition, we assessed the success of propensity score matching to balance covariates in the two groups using standardized differences. Standardized differences of < 10% support the assumption of balance between the 2 groups (37). Statistical significance was set at p < 0.05. Data were analyzed using the R package, version 3.3.1 (38).
Results
Association of Cardiovascular Risk Factors With Disease
A total of 519 cases presenting with a cardiovascular event were compared with a group of 760 controls. Both groups showed strong differences in age (61.3 in cases vs. 48.9 years in controls, p < 0.001), sex distribution (72.5% were male in the cases group vs. 45.4% in controls, p < 0.001), smoking (61.9 vs. 39.0%, p < 0.001), arterial hypertension (65.7 vs. 22.9%; p < 0.001), diabetes mellitus type 2 (T2DM) (40.3 vs. 27.1%; p < 0.001) and dyslipidemia (53.2 vs. 33.4%; p < 0.001) (Table 1).
Patients were also classified according to the type of event. Current smokers were more frequent among the STEMI patients. While body mass index differences were not significant among all four groups, diabetes and dyslipidemia appeared associated with angina rather than with acute events. Total leukocyte and platelet counts, along with levels of C reactive protein (CRP) and fibrinogen, increased in the NSTEMI group, showing maximal values in the STEMI patients (Table 2).
Association of Genetic Variants in the CCR2, CCL2, and CCL7 Genes With Coronary Heart Disease
The significant elevation of leukocytes and inflammatory markers in acute events emphasized the association of inflammation with coronary heart disease in this cohort. Given the potential role of the CCR2 chemokine receptor in progression of atherosclerosis, the hypothesis that genetic variants that could affect function of the receptor and its ligands could represent a risk factor in this setting was explored. Seventy cases and 30 healthy controls were selected to sequence all coding and flanking sequences of the CCR2, CCL2 and CCL7 genes. Assuming a dominant model, this approach would examine 70 “disease” and 130 “control” alleles.
A total of 16 genetic variants were found in CCR2, 2 in CCL2 and 4 in CCL7. Out of these, there was a single common non synonymous variant (rs1799864) affecting CCR2. Two novel variants were also found on this gene (3:46354176; 3:46354213), and a third one, an indel (rs372263390), had no frequency data available. Among all other variants found, only one appeared to be underrepresented in cases (2/70) vs. controls (6/30), rs17735770, located at the 3'untranslated region (UTR) of the CCL7 gene (see Supplementary Table).
Before attempting an analysis of the association with the rs17735770 variant with disease, and because of significant differences observed between affected and controls, a matching based on the propensity score was carried out. First, all related variants associated with disease were identified through a multivariant logistic regression analysis, using the complete enumeration algorithm and Akaike Information Criterion (AIC). Table 3 summarizes the logistic model in which was based the propensity score. The covariates that showed independent association with the coronary heart disease (CHD) variable were age (p < 0.001), sex (male) (p < 0.001), smoker status (p < 0.001), dyslipidemia (p = 0.053), arterial hypertension (p < 0.001) and type 2 diabetes mellitus (p = 0.025). From the 519 patients with CHD, a subset of 305 was obtained, which resulted in a set of homogeneous pairs of cases and controls in terms of the variables considered.
After propensity score matching of where the related variables were introduced, 305 patients were selected for which a homogenous pair was found for the referred variables using a caliper of 0.5. As shown (Table 4) a set of homogeneous pairs of cases and controls in terms of the variables considered was selected, with similar characteristics in patients with CHD vs. controls: mostly by men, with an average age of 57, and without significant differences in traditional risk factors such as smoking, dyslipidemia, arterial hypertension or type 2 diabetes, as all standardized differences diverged < 10%, supporting the assumption of a balance between both groups. After pairing, a conditional logistic regression analysis revealed that carrying the C variant at the rs17735770 CCL7 gene variant appeared as a protection factor for CHD (OR = 0.333, 95% CI = 0.157–0.709) (Table 5).
Predicted Functional Consequence of the rs17735770 Variant
The TargetScanHuman® v. 5.1 application (32), which uses all known miRNA heptanucleotide binding sequences present in the 3'UTR sequences from all known mammalian genes, was used to detect that the rs17735770 variant was in close proximity (6 base pairs away) of a high score (98/100) 8-mer binding site, for the microRNAs miR-23a and b, generically known as miR-23ab (Figure 1A). The CCL7 gene, as revealed by similarity BLAST (28) searches, was present, exclusively, in mammals and, despite being present in the least conserved 3'UTR of the gene, the immediate vicinity of the variant is highly conserved among different mammalian species (Figure 1B). Indeed, this consensus octamer binding site, at the 3'UTR of the CCL7 gene has been suggested to be recognized by miR-23 and regulate CCL7 expression in mice (39).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/908070/fcvm-09-908070-HTML/image_m/fcvm-09-908070-g001.jpg)
Figure 1. Predicted functional significance of the of the rs1773577 variant. (A) Context sequence of the rs1773577 variant aligned with mammalian species. Numbers indicate the position on the mRNA. The position affected by the rs1773577 variant is red, while the miR28ab binding site predicted by the TargetScan software.is shown in green. Non-conserved positions are indicated in gray. Reference sequences shown are: NM_006273.3 (Homo sapiens); XM_003278361.3 (Nomascus leucogenys); XM_002748371.3 (Callithrix jacchus); XM_008695392.1 (Ursus maritimus); XM_004747090.2 (Mustela putorius furo); XM_005597581.1 (Equus caballus); NM_013654.3 (Mus musculus). (B) Alignment of the predicted miR23 binding sequence with the corresponding mature miR23a and b. (C–E) RNAfold algorithm analysis of the effect of the rs1773577 variant in its RNA structural context. Panel (C) shows a minimum free energy (MFE) drawing of a local hairpin for the wt (T) allele pairing prediction, while panel (D) shows the same prediction for the variant (C). Color coding indicates base-pair probabilities as shown on the figure, at the bottom of panel (C). The green bar on the side of the predicted structures indicates the miR23ab predicted binding site. (E) Graphical representation of the level of entropy for each position along both structures, as indicated.
As the rs17735770 variant lied within the 3ÚTR of the CCL7 gene, in close proximity to a miR-23ab binding site, the possibility that this variant could affect the local RNA structure was investigated in silico using the RNAfold algorithm (34), which calculates the minimum free energy (MFE) structure for every particular sequence. Therefore, the effect of a single nucleotide variant on the secondary structure of an RNA molecule can be evaluated by comparing the resulting structures for each variant. As shown (Figures 1C,D), the RNA structure-analysis algorithm predicted the formation of a hairpin encompassing the miR-23ab binding site, shown with an estimation of base-pairing probability. The minimum free energy of the thermodynamic ensemble was −4.80 Kcal/mol for the most common T variant vs. −7.30 Kcal/mol for the C variant, indicating greater stability for the latter. This observation agrees with the calculated positional-entropy for each structure (Figure 1E). The minimum free energies for miR-23-mRNA dimer formation under both circumstances were calculated with RNAcofold (34), revealing predicted values of −13.95 Kcal/mol for the T allele and −16.13 Kcal/mol for the C variant allele, meaning that heterodimer formation with miR-23 was also favored by the rs17735770 variant.
Functional Evaluation of miR23 Binding
The functional relevance of the putative miR-23 binding site at the human CCL7 3'UTR was evaluated in the context of heterologous cellular assays. Firstly, synthetic DNA plasmid constructs containing the human CCL7 3'UTR downstream of the coding sequences of the Renilla luciferase gene were transfected into human embryonic kidney HEK-293 cells along with either miR-23a-3p mimic or nontargeting control mimic (CM) sequence As shown, overexpression of miR-23 significantly repressed CCL7 3'UTR activity (Figure 2A), while the mutation of the miR-23 binding site relieved the repression, which is consistent with its direct interaction with the human CCL7 3'UTR (Figure 2B) (40).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/908070/fcvm-09-908070-HTML/image_m/fcvm-09-908070-g002.jpg)
Figure 2. Functional evaluation of the miR-23 binding site at the human CCL7 3'UTR. (A) The human CCL7 gene 3'UTR was cloned downstream of a Renilla luciferase reporter and transfected into HEK-293 cells along with a miR-23a-3p mimic or a nontargeting control mimic (CM) as indicated. After 24 h, Renilla luciferase activity was measured and normalized according to Firefly luciferase activity, present in the same construct under a housekeeping promoter and a SV40 polyadenylation signal. The relative luciferase activity of the control was assigned to 100% activity as shown. (B) Alternatively, the miR-23 consensus site on the CCL7 gene 3'UTR was mutated from AATGTGAA to AATCAGGA and tested as above. (C) The rs17735770 variant was introduced into the CCL7 gene 3'UTR and cloned downstream of a Renilla luciferase reporter as the wt counterpart. Both the wt and rs17735770 variant versions were transfected into HEK-293 cells and tested as above in the presence of CM or miR-23a-3p mimics as indicated. A similar experimental approach was followed in COS-7 (D) and THP-1 cells treated with 100 nM PMA 24 h before transfections (E). Data are expressed as relative luciferase activities compared to the activity in control samples assigned to 100% activity co-transfected with an equal concentration of the CM and correspond to the means ± SEM from three experiments performed in triplicate. * P < 0.05 (significantly different from cells co-transfected with CM and the WT or mut 3′ UTR) ***p ≤ 0.001.
Once the miR-23 binding site was validated, similar analysis were performed to compare constructs containing either the common (wild type) or the human CCL7 3'UTR containing the rs17735770 variant. As shown, the construct containing the variant showed higher activity than the wt counterpart while inhibitory effect of miR-23a on the 3'UTR remained intact (Figure 2C).
To further explore this effect in other cell lines, COS-7 transformed African green monkey kidney fibroblasts were subjected to transfection with either common or variant 3'UTR along with CM or miR-23a mimic as indicated (Figure 2D). These analyses revealed that the levels of expression of the Renilla luciferase reporter under either the wild type or the variant CCL7 3'UTR were similar, and both were inhibited upon co-transfection with miR-23a.
Finally, we further studied CCL7 3'UTR activity in a native CCL7-expressing background, the human monocyte THP-1 cell line differentiated toward macrophage-like cells by addition of 100 nM PMA (41), leading to increased endogenous CCL7 expression (42). In this setting, the 3'UTR activity containing either the common or the variant were both significantly repressed by miR-23a and, in agreement with the results obtained in HEK-293 cells, a higher basal activity was detected for the construct containing the variant CCL7 3'UTR downstream of the reporter gene (Figure 2E).
Discussion
Atherosclerosis is the result of the chronification of an inflammatory process that fails to resolve and where monocytes play a key role from start. While much is known about the mechanisms that participate in the initiation and evolution of the disease, much less is known about the signals that may lead to its resolution (43). In this study, we have approached the study of classical risk factors associated with coronary heart disease in combination with genetic variants in genes that regulate monocyte function, in particular the CCR2 receptor and its ligands.
The characteristics of the population of patients presented in this study are similar to other populations and, accordingly, associated with known risk factors (44, 45). Elderly males, smokers, and patients presenting with arterial hypertension, diabetes or dyslipidemia were significantly associated with cardiovascular disease when compared to the control set. It was interesting to observe how specific risk factors segregated differentially according to the type of event suffered. While dyslipidemia and arterial hypertension were more prevalent among patients with angina pectoris, current smoking status and elevated inflammatory markers, including total leukocyte and monocyte counts, were more likely associated with acute coronary events in younger patients. In fact, cigarette smoking is the most common and modifiable risk factor in young adult patients with myocardial infarction (46) whereas diabetes and hypertension are less frequent (47). Also, single vessel coronary disease is more common than three vessel diseases in younger patients, as they are exposed to fewer cardiovascular risk factors, leading to reduced collateral circulation, a greater frequency of STEMI (48) and a higher myocardial area at risk ultimately leading to myocyte necrosis and to an increase in inflammatory markers.
A multivariant analysis further supported strong independent associations of age, smoking status, arterial hypertension, and dyslipidemia with cardiovascular risk in this population. It is well established that there is a genetic heritage behind some of these factors and, possibly, other physiological mechanisms that may also influence the susceptibility to suffer chronic diseases (49). It could be expected that specific genetic traits would become evident before age-dependent or other environmental factors exert their effect (50). Therefore, there might be a need to exclude general risk factors in order to detect genetic differences that could be risk factors per se. So, to approach this analysis, cases and controls were matched considering all the variables that were identified to be associated with the disease phenotype through a multivariant logistic regression analysis, generating a set of homogeneous pairs of cases and controls with similar characteristics.
The CCR2 receptor and its ligands were selected as candidate genes for a genetic study because of their importance in macrophage biology and because both agonist and antagonist functions are known for this setting, providing a window of opportunity for a negative feedback control mechanism that could lead to inflammation resolution (16). CCR2 is necessary for the development of atherosclerosis, as it has been demonstrated in murine models of atherosclerosis where either the CCR2 gene or the CCL2 gene, encoding for one of its ligands, are inactivated (11–15). Although later work reveals a key role of both CCR1 and CCR5 in monocyte recruitment leading to atherogenesis (51), it is also evident that, beyond facilitating the egress of monocytes from the bone marrow to increase their circulating numbers (9, 10, 52), CCR2 and its ligands are also autonomously involved in remodeling the endothelium in inflammatory processes and in the mobilization of monocytes associated with it (53–55).
From all CCR2, CCL2 and CCL7 variants identified, only rs17735770, located at the 3'end of the CCL7 gene, was found to segregate preferentially with the absence of cardiovascular disease, an observation that was further substantiated when analyzed in a greater case-control-matched population, where all pairs shared identical risk factors. This group was represented mostly by men with an average age of 57, indicating that this genetic variable influenced or modified disease risk before age-related factors take place.
The association of this variant, at the 3'end of the CCL7 gene, with the healthy phenotype is of interest and, accordingly, since the propensity score approach used in this study eliminates differences in other risk factors, setting a background to detect genetic modifiers, it has not been previously linked with cardiovascular disease, as suggested by GWAS exploring the association of genetic variants with a variety of cardiovascular disease-associated traits (56). This variant was predicted to alter the local RNA structure (Figure 1), lied close to a binding site for miR-23ab, and its presence resulted in increased stability of both an internal hairpin and miR-23ab dimer formation, and also in enhanced expression of a reporter in a heterologous assay, regardless of exogenous miR-23a addition.
Single nucleotide polymorphisms in regulatory non-coding RNAs such as miRNAs, as well as their direct binding sites or their vicinity at the 3‘UTR regions have been previously associated with human disease (57). Indeed, similar to our observations in the CCL7 3'UTR, a variant closely located to a miR-199a binding site at the 3'UTR of the HIF1A gene, is associated with pancreatic ductal adenocarcinoma risk and worse clinical outcomes (58). In this case, however, the rs17735770 variant did not hamper miR-23 binding, but rather resulted in increased CCL7 3'UTR activity independently of the participation of miR-23 in a heterologous assay in two different cell lines. These results suggest that variation resulting from rs17735770 might interfere with additional factors that regulate CCL7 expression at the posttranscriptional level. Given the function of miRNAs in the control of mRNA translation and stability (59), it is tempting to speculate that the presence of the rs17735770 variant results in increased expression of CCL7 at the sites of inflammation, where it could become an antagonist and ultimately contribute to the resolution of inflammation.
Little is known about the association of these miRNAs with cardiovascular disease and/or inflammation. Increased miR-23 levels are associated with the maintenance of the integrity of endothelial structures (60), and has been proposed to evaluate the presence and severity of coronary lesions in CHD patients (61) while, in an experimental autoinmmune murine model of encephalomyelitis, miR-23b suppresses leukocyte migration and pathogenesis by targeting CCL7 (39).
Our findings invite to revisit the experimentation on therapeutic approaches based on CCR2 antagonists for the prevention of cardiovascular disease. Because their influence in inflammation and monocyte biology, chemokines and their receptors have been the subject of pharmaceutical targeting for atherosclerosis and related conditions in recent years (62). In murine preclinical models of atherosclerosis, both directed antagonism to either CCR2 or a subset of chemokine receptors, including CCR2, result in a reduction in plaque formation (63, 64). However, clinical trials in humans have not been very informative in this area. A beneficial reduction of serum C Reactive Protein (CRP) has been recorded in patients treated with MLN1202, a humanized monoclonal antibody disrupting the interaction of CCR2 with CCL2 and, likely, with other ligands as well (65). These effects, accompanied by an overall reduction in circulating monocytes, are seen after a systemic exposure to MLN1202 in patients that have a promoter genetic variant at the CCL2 gene which results in increased CCL2 expression. In this setting, antagonism at CCR2 not only blocks monocyte migration into inflammatory lesions but also their release into circulation, also mediated by the CCR2-CCL2 interaction (6), and could lead to unwanted side effects derived from monocyte retrieval from circulation.
The impact of increased expression of CCL7 at inflammatory sites, and its function as an antagonist, as we speculate from our results, supports the idea that antagonism directed at the plaque should cause local benefit without compromising innate immunity elsewhere. Indeed, selective targeting of anti-inflammatory microRNAs to inflamed endothelium in a mouse model of atherosclerosis, lead to a local reduction in the expression of CCL2, among other chemokines, accompanied by reduced inflammation, macrophage numbers and plaque size (66). Recent observations revealing a circadian nature for the CCR2-CCL2 interaction at the plaque open the possibility for specific pharmaceutical intervention when inflammatory activity at the plaque is higher, a timely fashion to achieve selective targeting (67).
We acknowledge that further experimental evidence will be necessary to give a solid functional significance to the rs17735770 genetic variant as cardioprotective. The benefit of increased local CCL7 expression at the plaque and a differential expression between carriers and non-carriers of the rs17735770 variant is an issue that is difficult to approach in humans, but may inspire the design of novel therapeutic approaches.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Materials, further inquiries can be directed to the corresponding author.
Ethics Statement
The studies involving human participants were reviewed and approved by CEIm PROVINCIAL LAS PALMAS Barranco de la Ballena, s/n Edificio Anexo al Hospital Universitario de Gran Canaria Dr. Negrín 35019-Las Palmas de Gran Canaria, Spain, email: ceimprovlpa.scs@gobiernodecanarias.org. The patients/participants provided their written informed consent to participate in this study.
Author Contributions
JM-G, EM-Q, and FR-G collected, processed, and analyzed patient information. AP-G and CR designed and performed miRNA functional studies. PS-S was performed the statistical analysis. AD-C, PG-S, and AT performed sequencing and genotyping. MR supervised biochemical laboratory work. AT conceived the work, participated in data collection and analysis, prepared the figures and wrote the manuscript, which was finally supervised, modified, and finally approved by all authors.
Funding
This project has been initially funded by FUNCIS Project (PI 16/06) to EM-Q, by the Servicio Canario de Salud, and by grants RTI2018-095061-B-I00 (Convocatoria 2018 de proyectos de I+D+i ≪RETOS INVESTIGACIÓN≫) from the Ministerio de Ciencia, Investigación y Universidades, and TALENTO Grant 2017-T1/BMD-5333 from Comunidad de Madrid both to CR. AP-G is employed through PEDJ-2018-POST/BDM-8900 (Ayudas para la contratación de investigadores postdoctorales) from Consejería de Educación e Investigación from the Madrid Government, Spain, awarded to CR and AP-G.
Conflict of Interest
AD-C was employed by SECUGEN SL, Madrid, Spain and designed all the primers and performed all the sequencing associated with this project at SECUGEN SL, a service paid for by the authors. However, SECUGEN SL was not involved in the study design, collection, analysis, interpretation of data, the writing of this article or the decision to submit it for publication.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher's Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
We are indebted to patients and their families for their enthusiastic participation in this study. Dr. Cristina Balagué kindly reviewed the manuscript.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcvm.2022.908070/full#supplementary-material
References
1. Libby P. Mechanisms of acute coronary syndromes and their implications for therapy. N Engl J Med. (2013) 368:2004–13. doi: 10.1056/NEJMra1216063
2. Hansson GK, Libby P, Tabas I. Inflammation and plaque vulnerability. J Intern Med. (2015) 278:483–93. doi: 10.1111/joim.12406
3. Moore KJ, Tabas I. Macrophages in the pathogenesis of atherosclerosis. Cell. (2011) 145:341–55. doi: 10.1016/j.cell.2011.04.005
4. Hughes CE, Nibbs R. A guide to chemokines and their receptors. FEBS J. (2018) 285:2944–71. doi: 10.1111/febs.14466
5. Ingersoll MA, Platt AM, Potteaux S, Randolph GJ. Monocyte trafficking in acute and chronic inflammation. Trends Immunol. (2011) 32:470–7. doi: 10.1016/j.it.2011.05.001
6. Gautier EL, Jakubzick C, Randolph GJ. Regulation of the migration and survival of monocyte subsets by chemokine receptors and its relevance to atherosclerosis. Arterioscler Thromb Vasc Biol. (2009) 29:1412–8. doi: 10.1161/ATVBAHA.108.180505
7. Tacke F, Alvarez D, Kaplan TJ, Jakubzick C, Spanbroek R, Llodra J, et al. Monocyte subsets differentially employ CCR2, CCR5, and CX3CR1 to accumulate within atherosclerotic plaques. J Clin Invest. (2007) 117:185–94. doi: 10.1172/JCI28549
8. Combadière C, Potteaux S, Rodero M, Simon T, Pezard A, Esposito B, et al. Combined inhibition of CCL2, CX3CR1, and CCR5 abrogates Ly6C(hi) and Ly6C(lo) monocytosis and almost abolishes atherosclerosis in hypercholesterolemic mice. Circulation. (2008) 117:1649–57. doi: 10.1161/CIRCULATIONAHA.107.745091
9. Boring L, Gosling J, Chensue SW, Kunkel SL, Farese RV Jr, Broxmeyer HE, et al. Impaired monocyte migration and reduced type 1 (Th1) cytokine responses in C-C chemokine receptor 2 knockout mice. J Clin Invest. (1997) 100:2552–61. doi: 10.1172/JCI119798
10. Kuziel WA, Morgan SJ, Dawson TC, Griffin S, Smithies O, Ley K, et al. Severe reduction in leukocyte adhesion and monocyte extravasation in mice deficient in CC chemokine receptor 2. Proc Natl Acad Sci U S A. (1997) 94:12053–8. doi: 10.1073/pnas.94.22.12053
11. 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. doi: 10.1038/29788
12. 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. doi: 10.1016/S0021-9150(98)00318-9
13. Guo J, Van Eck M, Twisk J, Maeda N, Benson GM, Groot PH, et al. Transplantation of monocyte CC-chemokine receptor 2-deficient bone marrow into ApoE3-Leiden mice inhibits atherogenesis. Arterioscler Thromb Vasc Biol. (2003) 23:447–53. doi: 10.1161/01.ATV.0000058431.78833.F5
14. Gu L, Okada Y, Clinton SK, Gerard C, Sukhova GK, Libby P, et al. Absence of monocyte chemoattractant protein-1 reduces atherosclerosis in low density lipoprotein receptor-deficient mice. Mol Cell. (1998) 2:275–81. doi: 10.1016/S1097-2765(00)80139-2
15. Liu XL, Zhang PF, Ding SF, Wang Y, Zhang M, Zhao YX, et al. Local gene silencing of monocyte chemoattractant protein-1 prevents vulnerable plaque disruption in apolipoprotein E-knockout mice. PLoS ONE. (2012) 7:e33497. doi: 10.1371/journal.pone.0033497
16. McQuibban GA, Gong JH, Tam EM, McCulloch CA, Clark-Lewis I, Overall CM. Inflammation dampened by gelatinase A cleavage of monocyte chemoattractant protein-3. Science. (2000) 289:1202–6. doi: 10.1126/science.289.5482.1202
17. Ortlepp JR, Vesper K, Mevissen V, Schmitz F, Janssens U, Franke A, et al. Chemokine receptor (CCR2) genotype is associated with myocardial infarction and heart failure in patients under 65 years of age. J Mol Med. (2003) 81:363–7. doi: 10.1007/s00109-003-0435-x
18. González P, Alvarez R, Batalla A, Reguero JR, Alvarez V, Astudillo A, et al. Genetic variation at the chemokine receptors CCR5/CCR2 in myocardial infarction. Genes Immun. (2001) 2:191–5. doi: 10.1038/sj.gene.6363760
19. Szalai C, Duba J, Prohászka Z, Kalina A, Szabó T, Nagy B, et al. Involvement of polymorphisms in the chemokine system in the susceptibility for coronary artery disease (CAD). Coincidence of elevated Lp(a) and MCP-1−2518 G/G genotype in CAD patients. Atherosclerosis. (2001) 158:233–9. doi: 10.1016/S0021-9150(01)00423-3
20. Dow DJ, McMahon AD, Gray IC, Packard CJ, Groot PH. CCR2 and coronary artery disease: a woscops substudy. BMC Res Notes. (2010) 3:31. doi: 10.1186/1756-0500-3-31
21. Golbus JR, Stitziel NO, Zhao W, Xue C, Farrall M, McPherson R, et al. Common and rare genetic variation in CCR2, CCR5, or CX3CR1 and risk of atherosclerotic coronary heart disease and glucometabolic traits. Circ Cardiovasc Genet. (2016) 9:250–8. doi: 10.1161/CIRCGENETICS.115.001374
22. Valdes AM, Wolfe ML, O'Brien EJ, Spurr NK, Gefter W, Rut A, et al. Val64Ile polymorphism in the C-C chemokine receptor 2 is associated with reduced coronary artery calcification. Arterioscler Thromb Vasc Biol. (2002) 22:1924–8. doi: 10.1161/01.ATV.0000038486.48400.E7
23. McDermott DH, Yang Q, Kathiresan S, Cupples LA, Massaro JM, Keaney JF Jr, et al. CCL2 polymorphisms are associated with serum monocyte chemoattractant protein-1 levels and myocardial infarction in the Framingham Heart Study. Circulation. (2005) 112:1113–20. doi: 10.1161/CIRCULATIONAHA.105.543579
24. Iwai N, Kajimoto K, Kokubo Y, Okayama A, Miyazaki S, Nonogi H, et al. Assessment of genetic effects of polymorphisms in the MCP-1 gene on serum MCP-1 levels and myocardial infarction in Japanese. Circ J. (2006) 70:805–9. doi: 10.1253/circj.70.805
25. Jemaa R, Rojbani H, Kallel A, Ben Ali S, Feki M, Chabrak S, et al. Association between the−2518G/A polymorphism in the monocyte chemoattractant protein-1 (MCP-1) gene and myocardial infarction in Tunisian patients. Clin Chim Acta. (2008) 390:122–5. doi: 10.1016/j.cca.2008.01.004
26. Wägner AM, Wiebe JC, Boronat M, Saavedra P, Marrero D, Varillas F, et al. Insulin resistance explains the relationship between novel cardiovascular risk factors and hypertension. The telde study. J Endocrinol Invest. (2011) 34:e409–12. doi: 10.3275/7325
27. Rodríguez-González F, Martínez-Quintana E, Saavedra P, Medina-Gil JM, Riaño M, Garay-Sánchez P, et al. CYP2C19 or CYP3A5 genotyping does not predict clinical response to clopidogrel. J Clin Pharmacol. (2018) 58:1274–83. doi: 10.1002/jcph.1144
28. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. (1990) 215:403–10. doi: 10.1016/S0022-2836(05)80360-2
29. Available, online at: https://www.ncbi.nlm.nih.gov/snp/
30. Ye J, Coulouris G, Zaretskaya I, Cutcutache I, Rozen S, Madden TL. Primer-BLAST: a tool to design target-specific primers for polymerase chain reaction. BMC Bioinformatics. (2012) 13:134. doi: 10.1186/1471-2105-13-134
31. Martínez-Quintana E, Medina-Gil JM, Rodríguez-González F, Garay-Sánchez P, Limiñana JM, Saavedra P, et al. Positive clinical response to clopidogrel is independent of paraoxonase 1 Q192R and CYP2C19 genetic variants. J Clin Pharmacol. (2014) 54:843–9. doi: 10.1002/jcph.275
32. Friedman RC, Farh KK, Burge CB, Bartel DP. Most mammalian mRNAs are conserved targets of microRNAs. Genome Res. (2009) 19:92–105. doi: 10.1101/gr.082701.108
33. Available, online at: https://genome.ucsc.edu/
34. Gruber AR, Lorenz R, Bernhart SH, Neuböck R, Hofacker IL. The Vienna RNA websuite. Nucleic Acids Res. (2008) 36:W70–4. doi: 10.1093/nar/gkn188
35. Konietschke F, Placzek M, Schaarsachmidt F, Hothorn LA. nparcomp: An R software package for nonparametric multiple comparisons and simultaneous confidence intervals. J Stat Softw. (2015) 64:1–19. doi: 10.18637/jss.v064.i09
36. Morgan JA, Tatar JF. Calculation of the residual sum of squares for all possible regressions. Technometrics. (1972) 14:317–25. doi: 10.1080/00401706.1972.10488918
37. Sekhon JS. Multivariate and propensity score matching software with automated balance optimization: the matching package for R. J Stat Softw. (2011) 42:1–52. doi: 10.18637/jss.v042.i,07
38. R Development Core Team. R: A language environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing (2016). Available online at: https://www.Rproject.org/
39. Zhang Y, Han JJ, Liang XY, Zhao L, Zhang F, Rasouli J, et al. miR-23b Suppresses Leukocyte Migration and Pathogenesis of Experimental Autoimmune Encephalomyelitis by Targeting CCL7. Mol Ther. (2018) 26:582–92. doi: 10.1016/j.ymthe.2017.11.013
40. Fernández-de Frutos M, Galán-Chilet I, Goedeke L, Kim B, Pardo-Marqués V, Pérez-García A, et al. MicroRNA 7 impairs insulin signaling and regulates Aβ levels through posttranscriptional regulation of the insulin receptor substrate 2, insulin receptor, insulin-degrading enzyme, and liver X receptor pathway. Mol Cell Biol. (2019) 39: e00170–19. doi: 10.1128/MCB.00170-19
41. Park EK, Jung HS, Yang HI, Yoo MC, Kim C, Kim KS. Optimized THP-1 differentiation is required for the detection of responses to weak stimuli. Inflamm Res. (2007) 56:45–50. doi: 10.1007/s00011-007-6115-5
42. Sozzani S, Allavena P, Proost P, Van Damme J, Mantovani A. Chemokines as targets for pharmacological intervention. Prog Drug Res. (1996) 47:53–80. doi: 10.1007/978-3-0348-8998-8_2
43. Libby P, Tabas I, Fredman G, Fisher EA. Inflammation and its resolution as determinants of acute coronary syndromes. Circ Res. (2014) 114:1867–79. doi: 10.1161/CIRCRESAHA.114.302699
44. Kannel WB, Dawber TR, Kagan A, Revotskie N, Stokes J 3rd. Factors of risk in the development of coronary heart disease–six year follow-up experience. The Framingham Study. Ann Intern Med. (1961) 55:33–50. doi: 10.7326/0003-4819-55-1-33
45. Wilson PW, D'Agostino RB, Levy D, Belanger AM, Silbershatz H, Kannel WB. Prediction of coronary heart disease using risk factor categories. Circulation. (1998) 97:1837–47. doi: 10.1161/01.CIR.97.18.1837
46. Zimmerman FH, Cameron A, Fisher LD, Ng G. Myocardial infarction in young adults: angiographic characterization, risk factors and prognosis (Coronary Artery Surgery Study Registry). J Am Coll Cardiol. (1995) 26:654–61. doi: 10.1016/0735-1097(95)00254-2
47. Dolder MA, Oliver MF. Myocardial infarction in young men. Study of risk factors in nine countries. Br Heart J. (1975) 37:493–503. doi: 10.1136/hrt.37.5.493
48. Avezum A, Makdisse M, Spencer F, Gore JM, Fox KA, Montalescot G, et al. Impact of age on management and outcome of acute coronary syndrome: observations from the Global Registry of Acute Coronary Events (GRACE). Am Heart J. (2005) 149:67–73. doi: 10.1016/j.ahj.2004.06.003
49. Kathiresan S, Srivastava D. Genetics of human cardiovascular disease. Cell. (2012) 148:1242–57. doi: 10.1016/j.cell.2012.03.001
50. Lloyd-Jones DM, Nam BH, D'Agostino RBSr, Levy D, Murabito JM, Wang TJ, et al. Parental cardiovascular disease as a risk factor for cardiovascular disease in middle-aged adults: a prospective study of parents and offspring. JAMA. (2004) 291:2204–11. doi: 10.1001/jama.291.18.2204
51. Soehnlein O, Drechsler M, Döring Y, Lievens D, Hartwig H, Kemmerich K, et al. Distinct functions of chemokine receptor axes in the atherogenic mobilization and recruitment of classical monocytes. EMBO Mol Med. (2013) 5:471–81. doi: 10.1002/emmm.201201717
52. Tsou CL, Peters W, Si Y, Slaymaker S, Aslanian AM, Weisberg SP, et al. 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. doi: 10.1172/JCI29919
53. Schepers A, Eefting D, Bonta PI, Grimbergen JM, de Vries MR, van Weel V, et al. Anti-MCP-1 gene therapy inhibits vascular smooth muscle cells proliferation and attenuates vein graft thickening both in vitro and in vivo. Arterioscler Thromb Vasc Biol. (2006) 26:2063–9. doi: 10.1161/01.ATV.0000235694.69719.e2
54. Inoue S, Egashira K, Ni W, Kitamoto S, Usui M, Otani K, et al. Anti-monocyte chemoattractant protein-1 gene therapy limits progression and destabilization of established atherosclerosis in apolipoprotein E-knockout mice. Circulation. (2002) 106:2700–6. doi: 10.1161/01.CIR.0000038140.80105.AD
55. Schober A, Zernecke A, Liehn EA, von Hundelshausen P, Knarren S, Kuziel WA, et al. Crucial role of the CCL2/CCR2 axis in neointimal hyperplasia after arterial injury in hyperlipidemic mice involves early monocyte recruitment and CCL2 presentation on platelets. Circ Res. (2004) 95:1125–33. doi: 10.1161/01.RES.0000149518.86865.3e
56. Available, online at: https://cvd.hugeamp.org/variant.html?variant=rs17735770
57. Moszyńska A, Gebert M, Collawn JF, Bartoszewski R. SNPs in microRNA target sites and their potential role in human disease. Open Biol. (2017) 7:170019. doi: 10.1098/rsob.170019
58. Wang X, Ren H, Zhao T, Ma W, Dong J, Zhang S, et al. Single nucleotide polymorphism in the microRNA-199a binding site of HIF1A gene is associated with pancreatic ductal adenocarcinoma risk and worse clinical outcomes. Oncotarget. (2016) 7:13717–29. doi: 10.18632/oncotarget.7263
59. Iwakawa HO, Tomari Y. The Functions of MicroRNAs: mRNA Decay and Translational Repression. Trends Cell Biol. (2015) 25:651–65. doi: 10.1016/j.tcb.2015.07.011
60. Wang KC, Garmire LX, Young A, Nguyen P, Trinh A, Subramaniam S, et al. Role of microRNA-23b in flow-regulation of Rb phosphorylation and endothelial cell growth. Proc Natl Acad Sci U S A. (2010) 107:3234–9. doi: 10.1073/pnas.0914825107
61. Di Y, Zhang D, Hu T, Li D. miR-23 regulate the pathogenesis of patients with coronary artery disease. Int J Clin Exp Med. (2015) 8:11759–69.
62. Márquez AB, van der Vorst E, Maas SL. Key Chemokine Pathways in Atherosclerosis and Their Therapeutic Potential. J Clin Med. (2021) 10:3825. doi: 10.3390/jcm10173825
63. Ni W, Egashira K, Kitamoto S, Kataoka C, Koyanagi M, Inoue S, et al. New anti-monocyte chemoattractant protein-1 gene therapy attenuates atherosclerosis in apolipoprotein E-knockout mice. Circulation. (2001) 103:2096–101. doi: 10.1161/01.CIR.103.16.2096
64. Ravindran D, Ridiandries A, Vanags LZ, Henriquez R, Cartland S, Tan JT, et al. Chemokine binding protein 'M3' limits atherosclerosis in apolipoprotein E-/- mice. PLoS ONE. (2017) 12:e0173224. doi: 10.1371/journal.pone.0173224
65. Gilbert J, Lekstrom-Himes J, Donaldson D, Lee Y, Hu M, Xu J, et al. Effect of CC chemokine receptor 2 CCR2 blockade on serum C-reactive protein in individuals at atherosclerotic risk and with a single nucleotide polymorphism of the monocyte chemoattractant protein-1 promoter region. Am J Cardiol. (2011) 107:906–11. doi: 10.1016/j.amjcard.2010.11.005
66. Ma S, Tian XY, Zhang Y, Mu C, Shen H, Bismuth J, et al. E-selectin-targeting delivery of microRNAs by microparticles ameliorates endothelial inflammation and atherosclerosis. Sci Rep. (2016) 6:22910. doi: 10.1038/srep22910
Keywords: chemokine, CCL7, genetics, microRNA (miRNA), cardiovascular disease
Citation: Medina-Gil JM, Pérez-García A, Saavedra-Santana P, Díaz-Carrasco A, Martínez-Quintana E, Rodríguez-González F, Ramírez CM, Riaño M, Garay-Sánchez P and Tugores A (2022) A Common Variant at the 3'untranslated Region of the CCL7 Gene (rs17735770) Is Associated With Decreased Susceptibility to Coronary Heart Disease. Front. Cardiovasc. Med. 9:908070. doi: 10.3389/fcvm.2022.908070
Received: 30 March 2022; Accepted: 12 May 2022;
Published: 31 May 2022.
Edited by:
Shizuka Uchida, Aalborg University Copenhagen, DenmarkReviewed by:
Subhi Marwari, University of Pennsylvania, United StatesGiuseppe Militello, Mirimus Inc., United States
Copyright © 2022 Medina-Gil, Pérez-García, Saavedra-Santana, Díaz-Carrasco, Martínez-Quintana, Rodríguez-González, Ramírez, Riaño, Garay-Sánchez and Tugores. 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: Antonio Tugores, atugores@yahoo.com
†These authors have contributed equally to this work and share first authorship