Identification of the BCAR1-CFDP1-TMEM170A locus as a determinant of carotid intima-media thickness and coronary artery disease risk

Genetic Associations of Circulating Cardiovascular Proteins With Gestational Hypertension and Preeclampsia

Importance

Hypertensive disorders of pregnancy (HDPs), including gestational hypertension and preeclampsia, are important contributors to maternal morbidity and mortality worldwide. In addition, women with HDPs face an elevated long-term risk of cardiovascular disease.

Objective

To identify proteins in the circulation associated with HDPs.

Design, Setting, and Participants

Two-sample mendelian randomization (MR) tested the associations of genetic instruments for cardiovascular disease-related proteins with gestational hypertension and preeclampsia. In downstream analyses, a systematic review of observational data was conducted to evaluate the identified proteins' dynamics across gestation in hypertensive vs normotensive pregnancies, and phenome-wide MR analyses were performed to identify potential non-HDP-related effects associated with the prioritized proteins. Genetic association data for cardiovascular disease-related proteins were obtained from the Systematic and Combined Analysis of Olink Proteins (SCALLOP) consortium. Genetic association data for the HDPs were obtained from recent European-ancestry genome-wide association study meta-analyses for gestational hypertension and preeclampsia. Study data were analyzed October 2022 to October 2023.

Exposures

Genetic instruments for 90 candidate proteins implicated in cardiovascular diseases, constructed using cis-protein quantitative trait loci (cis-pQTLs).

Main Outcomes and Measures

Gestational hypertension and preeclampsia.

Results

Genetic association data for cardiovascular disease-related proteins were obtained from 21 758 participants from the SCALLOP consortium. Genetic association data for the HDPs were obtained from 393 238 female individuals (8636 cases and 384 602 controls) for gestational hypertension and 606 903 female individuals (16 032 cases and 590 871 controls) for preeclampsia. Seventy-five of 90 proteins (83.3%) had at least 1 valid cis-pQTL. Of those, 10 proteins (13.3%) were significantly associated with HDPs. Four were robust to sensitivity analyses for gestational hypertension (cluster of differentiation 40, eosinophil cationic protein [ECP], galectin 3, N-terminal pro-brain natriuretic peptide [NT-proBNP]), and 2 were robust for preeclampsia (cystatin B, heat shock protein 27 [HSP27]). Consistent with the MR findings, observational data revealed that lower NT-proBNP (0.76- to 0.88-fold difference vs no HDPs) and higher HSP27 (2.40-fold difference vs no HDPs) levels during the first trimester of pregnancy were associated with increased risk of HDPs, as were higher levels of ECP (1.60-fold difference vs no HDPs). Phenome-wide MR analyses identified 37 unique non-HDP-related protein-disease associations, suggesting potential on-target effects associated with interventions lowering HDP risk through the identified proteins.

Conclusions and Relevance

Study findings suggest genetic associations of 4 cardiovascular disease-related proteins with gestational hypertension and 2 associated with preeclampsia. Future studies are required to test the efficacy of targeting the corresponding pathways to reduce HDP risk.

Cfdp1 Is Essential for Cardiac Development and Function

Cardiovascular diseases (CVDs) are the prevalent cause of mortality worldwide. A combination of environmental and genetic effectors modulates the risk of developing them. Thus, it is vital to identify candidate genes and elucidate their role in the manifestation of the disease. Large-scale human studies have revealed the implication of Craniofacial Development Protein 1 (CFDP1) in Coronary Artery Disease (CAD). CFDP1 belongs to the evolutionary conserved Bucentaur (BCNT) family, and to date, its function and mechanism of action in Cardiovascular Development are still unclear. We utilized zebrafish to investigate the role of cfdp1 in the developing heart due to the high genomic homology, similarity in heart physiology, and ease of experimental manipulations. We showed that cfdp1 was expressed during development, and we tested two morpholinos and generated a cfdp1 mutant line. The cfdp1-/- embryos developed arrhythmic hearts and exhibited defective cardiac performance, which led to a lethal phenotype. Findings from both knockdown and knockout experiments showed that abrogation of cfdp1 leads to downregulation of Wnt signaling in embryonic hearts during valve development but without affecting Notch activation in this process. The cfdp1 zebrafish mutant line provides a valuable tool for unveiling the novel mechanism of regulating cardiac physiology and function. cfdp1 is essential for cardiac development, a previously unreported phenotype most likely due to early lethality in mice. The detected phenotype of bradycardia and arrhythmias is an observation with potential clinical relevance for humans carrying heterozygous CFDP1 mutations and their risk of developing CAD.

Causal analysis of plasma IL-8 on carotid intima media thickness, a measure of subclinical atherosclerosis

BACKGROUND

We investigated the causality of IL-8 on carotid intima-media thickness (c-IMT), a measure of sub-clinical atherosclerosis.

METHODS

The IMPROVE is a multicenter European study (n = 3,711). The association of plasma IL-8 with c-IMT (mm) was estimated by quantile regression. Genotyping was performed using the Illumina CardioMetabo and Immuno chips. Replication was attempted in three independent studies and a meta-analysis was performed using a random model.

RESULTS

In IMPROVE, each unit increase in plasma IL-8 was associated with an increase in median c-IMT measures (all p<0·03) in multivariable analyses. Linear regression identified rs117518778 and rs8057084 as associated with IL-8 levels and with measures of c-IMT. The two SNPs were combined in an IL-8-increasing genetic risk that showed causality of IL-8 on c-IMT in IMPROVE and in the UK Biobank (n = 22,179). The effect of IL-8 on c-IMT measures was confirmed in PIVUS (n = 1,016) and MDCCC (n = 6,103). The association of rs8057084 with c-IMT was confirmed in PIVUS and UK Biobank with a pooled estimate effect (β) of -0·006 with 95%CI (-0·008- -0·003).

CONCLUSION

Our results indicate that genetic variants associated with plasma IL-8 also associate with c-IMT. However, we cannot infer causality of this association, as these variants lie outside of the IL8 locus.

A Review of Vascular Traits and Assessment Techniques, and Their Heritability

Various tools are available to assess atherosclerosis, arterial stiffening, and endothelial function. They offer utility in the assessment of hypertensive phenotypes, in cardiovascular risk prediction, and as surrogate endpoints in clinical trials. We explore the relative influence of participant genetics, with reference to large-scale genomic studies, population-based cohorts, and candidate gene studies. We find heritability estimates highest for carotid intima-media thickness (CIMT 35–65%), followed by pulse wave velocity as a measure of arterial stiffness (26–43%), and flow mediated dilatation as a surrogate for endothelial function (14–39%); data were lacking for peripheral artery tonometry. We furthermore examine genes and polymorphisms relevant to each technique. We conclude that CIMT and pulse wave velocity dominate the existing evidence base, with fewer published genomic linkages for measures of endothelial function. We finally make recommendations regarding planning and reporting of data relating to vascular assessment techniques, particularly when genomic data are also available, to facilitate integration of these tools into cardiovascular disease research.

Twenty-Five Novel Loci for Carotid Intima-Media Thickness: A Genome-Wide Association Study in >45 000 Individuals and Meta-Analysis of >100 000 Individuals

OBJECTIVE

Carotid artery intima-media thickness (cIMT) is a widely accepted marker of subclinical atherosclerosis. Twenty susceptibility loci for cIMT were previously identified and the identification of additional susceptibility loci furthers our knowledge on the genetic architecture underlying atherosclerosis.

APPROACH AND RESULTS

We performed 3 genome-wide association studies in 45 185 participants from the UK Biobank study who underwent cIMT measurements and had data on minimum, mean, and maximum thickness. We replicated 15 known loci and identified 20 novel loci associated with cIMT at P<5×10-8. Seven novel loci (ZNF385D, ADAMTS9, EDNRA, HAND2, MYOCD, ITCH/EDEM2/MMP24, and MRTFA) were identified in all 3 phenotypes. An additional new locus (LOXL1) was identified in the meta-analysis of the 3 phenotypes. Sex interaction analysis revealed sex differences in 7 loci including a novel locus (SYNE3) in males. Meta-analysis of UK Biobank data with a previous meta-analysis led to identification of three novel loci (APOB, FIP1L1, and LOXL4). Transcriptome-wide association analyses implicated additional genes ARHGAP42, NDRG4, and KANK2. Gene set analysis showed an enrichment in extracellular organization and the PDGF (platelet-derived growth factor) signaling pathway. We found positive genetic correlations of cIMT with coronary artery disease rg=0.21 (P=1.4×10-7), peripheral artery disease rg=0.45 (P=5.3×10-5), and systolic blood pressure rg=0.30 (P=4.0×10-18). A negative genetic correlation between average of maximum cIMT and high-density lipoprotein was found rg=-0.12 (P=7.0×10-4).

CONCLUSIONS

Genome-wide association meta-analyses in >100 000 individuals identified 25 novel loci associated with cIMT providing insights into genes and tissue-specific regulatory mechanisms of proatherosclerotic processes. We found evidence for shared biological mechanisms with cardiovascular diseases.

Associations of carotid intima media thickness with gene expression in whole blood and genetically predicted gene expression across 48 tissues

Carotid intima media thickness (cIMT) is a biomarker of subclinical atherosclerosis and a predictor of future cardiovascular events. Identifying associations between gene expression levels and cIMT may provide insight to atherosclerosis etiology. Here, we use two approaches to identify associations between mRNA levels and cIMT: differential gene expression analysis in whole blood and S-PrediXcan. We used microarrays to measure genome-wide whole blood mRNA levels of 5647 European individuals from four studies. We examined the association of mRNA levels with cIMT adjusted for various potential confounders. Significant associations were tested for replication in three studies totaling 3943 participants. Next, we applied S-PrediXcan to summary statistics from a cIMT genome-wide association study (GWAS) of 71 128 individuals to estimate the association between genetically determined mRNA levels and cIMT and replicated these analyses using S-PrediXcan on an independent GWAS on cIMT that included 22 179 individuals from the UK Biobank. mRNA levels of TNFAIP3, CEBPD and METRNL were inversely associated with cIMT, but these associations were not significant in the replication analysis. S-PrediXcan identified associations between cIMT and genetically determined mRNA levels for 36 genes, of which six were significant in the replication analysis, including TLN2, which had not been previously reported for cIMT. There was weak correlation between our results using differential gene expression analysis and S-PrediXcan. Differential expression analysis and S-PrediXcan represent complementary approaches for the discovery of associations between phenotypes and gene expression. Using these approaches, we prioritize TNFAIP3, CEBPD, METRNL and TLN2 as new candidate genes whose differential expression might modulate cIMT.

The Association between HDL-C and Subclinical Atherosclerosis Depends on CETP Plasma Concentration: Insights from the IMPROVE Study

The impact of cholesteryl ester transfer protein (CETP) on atherosclerosis is highly debated. This study aimed to investigate the associations between plasma CETP or CETP genotypes and carotid intima-media thickness (cIMT) and the influence of high-density lipoprotein cholesterol (HDL-C) on these associations. Plasma CETP and HDL-C concentrations were measured in 552 subjects free of any pharmacological treatment from the IMPROVE cohort, which includes 3711 European subjects at high cardiovascular risk. CETP single-nucleotide polymorphisms (SNPs) and cIMT measures (cIMT_max; cIMT_mean-max of bifurcations, common and internal carotids; plaque-free common carotid [PF CC]-IMT_mean) were available for the full cohort. In drug-free subjects, plasma CETP correlated with HDL-C levels (r = 0.19, p < 0.0001), but not with cIMT variables. When stratified according to HDL-C quartiles, CETP positively correlated with cIMT_max and cIMT_mean-max, but not with PF CC-IMT_mean, in the top HDL-C quartile only. Positive associations between the CETP concentration and cIMT_max or cIMT_mean-max were found in the top HDL-C quartile, whereas HDL-C levels were negatively correlated with cIMT_max and cIMT_mean-max when the CETP concentration was below the median (HDL-C × CETP interaction, p = 0.001 and p = 0.003 for cIMT_max and cIMT_mean-max, respectively). In the full cohort, three CETP SNPs (rs34760410, rs12920974, rs12708968) were positively associated with cIMT_max. rs12444708 exhibited a significant interaction with HDL-C levels in the prediction of cIMT_max. In conclusion, a significant interplay was found between plasma CETP and/or CETP genotype and HDL-C in the prediction of carotid plaque thickness, as indexed by cIMT_max. This suggests that the association of HDL-C with carotid atherosclerosis is CETP-dependent.

A fast and efficient colocalization algorithm for identifying shared genetic risk factors across multiple traits

Genome-wide association studies (GWAS) have identified thousands of genomic regions affecting complex diseases. The next challenge is to elucidate the causal genes and mechanisms involved. One approach is to use statistical colocalization to assess shared genetic aetiology across multiple related traits (e.g. molecular traits, metabolic pathways and complex diseases) to identify causal pathways, prioritize causal variants and evaluate pleiotropy. We propose HyPrColoc (Hypothesis Prioritisation for multi-trait Colocalization), an efficient deterministic Bayesian algorithm using GWAS summary statistics that can detect colocalization across vast numbers of traits simultaneously (e.g. 100 traits can be jointly analysed in around 1 s). We perform a genome-wide multi-trait colocalization analysis of coronary heart disease (CHD) and fourteen related traits, identifying 43 regions in which CHD colocalized with ≥1 trait, including 5 previously unknown CHD loci. Across the 43 loci, we further integrate gene and protein expression quantitative trait loci to identify candidate causal genes.

Sex-dimorphic genetic effects and novel loci for fasting glucose and insulin variability

Differences between sexes contribute to variation in the levels of fasting glucose and insulin. Epidemiological studies established a higher prevalence of impaired fasting glucose in men and impaired glucose tolerance in women, however, the genetic component underlying this phenomenon is not established. We assess sex-dimorphic (73,089/50,404 women and 67,506/47,806 men) and sex-combined (151,188/105,056 individuals) fasting glucose/fasting insulin genetic effects via genome-wide association study meta-analyses in individuals of European descent without diabetes. Here we report sex dimorphism in allelic effects on fasting insulin at IRS1 and ZNF12 loci, the latter showing higher RNA expression in whole blood in women compared to men. We also observe sex-homogeneous effects on fasting glucose at seven novel loci. Fasting insulin in women shows stronger genetic correlations than in men with waist-to-hip ratio and anorexia nervosa. Furthermore, waist-to-hip ratio is causally related to insulin resistance in women, but not in men. These results position dissection of metabolic and glycemic health sex dimorphism as a steppingstone for understanding differences in genetic effects between women and men in related phenotypes.

Genetic Variants Associated with Non-Alcoholic Fatty Liver Disease Do Not Associate with Measures of Sub-Clinical Atherosclerosis: Results from the IMPROVE Study

Non-alcoholic fatty liver disease (NAFLD) and atherosclerosis-related cardiovascular diseases (CVD) share common metabolic pathways. We explored the association between three NAFLD-associated single nucleotide polymorphisms (SNPs) rs738409, rs10401969, and rs1260326 with sub-clinical atherosclerosis estimated by the carotid intima-media thickness (c-IMT) and the inter-adventitia common carotid artery diameter (ICCAD) in patients free from clinically overt NAFLD and CVD. The study population is the IMPROVE, a multicenter European study (n = 3711). C-IMT measures and ICCAD were recorded using a standardized protocol. Linear regression with an additive genetic model was used to test for association of the three SNPs with c-IMT and ICCAD. In secondary analyses, the association of the three SNPs with c-IMT and ICCAD was tested after stratification by alanine aminotransferase levels (ALT). No associations were found between rs738409, rs1260326, rs10401969, and c-IMT or ICCAD. Rs738409-G and rs10401969-C were associated with ALT levels (p < 0.001). In patients with ALT levels above 28 U/L (highest quartile), we observed an association between rs10401969-C and c-IMT measures of c-IMT_max and c-IMT_mean-max (p = 0.018 and 0.021, respectively). In conclusion, NAFLD-associated SNPs do not associate with sub-clinical atherosclerosis measures. However, our results suggest a possible mediating function of impaired liver function on atherosclerosis development.

Carotid Intima-Media Thickness: Novel Loci, Sex-Specific Effects, and Genetic Correlations With Obesity and Glucometabolic Traits in UK Biobank

OBJECTIVE

Atherosclerosis is the underlying cause of most cardiovascular disease, but mechanisms underlying atherosclerosis are incompletely understood. Ultrasound measurement of the carotid intima-media thickness (cIMT) can be used to measure vascular remodeling, which is indicative of atherosclerosis. Genome-wide association studies have identified many genetic loci associated with cIMT, but heterogeneity of measurements collected by many small cohorts have been a major limitation in these efforts. Here, we conducted genome-wide association analyses in UKB (UK Biobank; N=22 179), the largest single study with consistent cIMT measurements. Approach and Results: We used BOLT-LMM software to run linear regression of cIMT in UKB, adjusted for age, sex, and genotyping chip. In white British participants, we identified 5 novel loci associated with cIMT and replicated most previously reported loci. In the first sex-specific analyses of cIMT, we identified a locus on chromosome 5, associated with cIMT in women only and highlight VCAN as a good candidate gene at this locus. Genetic correlations with body mass index and glucometabolic traits were also observed. Two loci influenced risk of ischemic heart disease.

CONCLUSIONS

These findings replicate previously reported associations, highlight novel biology, and provide new directions for investigating the sex differences observed in cardiovascular disease presentation and progression.

Subclinical atherosclerosis and its progression are modulated by PLIN2 through a feed-forward loop between LXR and autophagy

BACKGROUND

Hyperlipidaemia is a major risk factor for cardiovascular disease, and atherosclerosis is the underlying cause of both myocardial infarction and stroke. We have previously shown that the Pro251 variant of perilipin-2 reduces plasma triglycerides and may therefore be beneficial to reduce atherosclerosis development.

OBJECTIVE

We sought to delineate putative beneficial effects of the Pro251 variant of perlipin-2 on subclinical atherosclerosis and the mechanism by which it acts.

METHODS

A pan-European cohort of high-risk individuals where carotid intima-media thickness has been assessed was adopted. Human primary monocyte-derived macrophages were prepared from whole blood from individuals recruited by perilipin-2 genotype or from buffy coats from the Karolinska University hospital blood central.

RESULTS

The Pro251 variant of perilipin-2 is associated with decreased intima-media thickness at baseline and over 30 months of follow-up. Using human primary monocyte-derived macrophages from carriers of the beneficial Pro251 variant, we show that this variant increases autophagy activity, cholesterol efflux and a controlled inflammatory response. Through extensive mechanistic studies, we demonstrate that increase in autophagy activity is accompanied with an increase in liver-X-receptor (LXR) activity and that LXR and autophagy reciprocally activate each other in a feed-forward loop, regulated by CYP27A1 and 27OH-cholesterol.

CONCLUSIONS

For the first time, we show that perilipin-2 affects susceptibility to human atherosclerosis through activation of autophagy and stimulation of cholesterol efflux. We demonstrate that perilipin-2 modulates levels of the LXR ligand 27OH-cholesterol and initiates a feed-forward loop where LXR and autophagy reciprocally activate each other; the mechanism by which perilipin-2 exerts its beneficial effects on subclinical atherosclerosis.

Targeted sequencing of linkage region in Dominican families implicates PRIMA1 and the SPATA7-PTPN21-ZC3H14-EML5-TTC8 locus in carotid-intima media thickness and atherosclerotic events

Carotid intima-media thickness (cIMT) is a subclinical marker for atherosclerosis. Previously, we reported a quantitative trait locus (QTL) for total cIMT on chromosome 14q and identified PRiMA1, FOXN3 and CCDC88C as candidate genes using a common variants (CVs)-based approach. Herein, we further evaluated the genetic contribution of the QTL to cIMT by resequencing. We sequenced all exons within the QTL and genomic regions of PRiMA1, FOXN3 and CCDC88C in Dominican families with evidence for linkage to the QTL. Unrelated Dominicans from the Northern Manhattan Study (NOMAS) were used for validation. Single-variant-based and gene-based analyses were performed for CVs and rare variants (RVs). The strongest evidence for association with CVs was found in PRiMA1 (p = 8.2 × 10⁻⁵ in families, p = 0.01 in NOMAS at rs12587586), and in the five-gene cluster SPATA7-PTPN21-ZC3H14-EML5-TTC8 locus (p = 1.3 × 10⁻⁴ in families, p = 0.01 in NOMAS at rs2274736). No evidence for association with RVs was found in PRiMA1. The top marker from previous study in PRiMA1 (rs7152362) was associated with fewer atherosclerotic events (OR = 0.67; p = 0.02 in NOMAS) and smaller cIMT (β = −0.58, p = 2.8 × 10⁻⁴ in Family). Within the five-gene cluster, evidence for association was found for exonic RVs (p = 0.02 in families, p = 0.28 in NOMAS), which was enriched among RVs with higher functional potentials (p = 0.05 in NOMAS for RVs in the top functional tertile). In summary, targeted resequencing provided validation and novel insights into the genetic architecture of cIMT, suggesting stronger effects for RVs with higher functional potentials. Furthermore, our data support the clinical relevance of CVs associated with subclinical atherosclerosis.

Genetic overlap between birthweight and adult cardiometabolic diseases has implications for genomic medicine

Before implementing therapeutic genomic interventions for optimizing health in early life, comprehensive understanding of their effect on several traits across the life course is warranted. Abnorml  birthweight is associated with cardiometabolic disease risk in adulthood; however, the extent of genetic pleiotropy in the association has not been comprehensively investigated. We tested for pleiotropy and enrichment of functional loci between birthweight and 15 cardiometabolic disease traits (CMD). We found significantly abundant genetic pleiotropy (P < 3.3 × 10⁻³) and enrichment of functional annotations (P < 3.3 × 10⁻³) in loci influencing both birthweight and CMD. We did not observe consistent effect directions of pleiotropic loci on the traits. A total of 67 genetic loci, of which 65 loci have been reported in previous genome-wide association studies, were associated with both birthweight and CMD at a false discovery rate of 5%. Two novel loci were associated with birthweight and adult coronary artery disease (rs2870463 in CTRB1) and with birthweight and adult waist circumference (rs12704673 in CALCR). Both loci are known to have regulatory effects on expression of nearby genes. In all, our findings revealed pervasive genetic pleiotropy in early growth and adulthood cardiometabolic diseases, implying the need for caution when considering genetic loci as therapeutic targets.

GWAS and colocalization analyses implicate carotid intima-media thickness and carotid plaque loci in cardiovascular outcomes

Carotid artery intima media thickness (cIMT) and carotid plaque are measures of subclinical atherosclerosis associated with ischemic stroke and coronary heart disease (CHD). Here, we undertake meta-analyses of genome-wide association studies (GWAS) in 71,128 individuals for cIMT, and 48,434 individuals for carotid plaque traits. We identify eight novel susceptibility loci for cIMT, one independent association at the previously-identified PINX1 locus, and one novel locus for carotid plaque. Colocalization analysis with nearby vascular expression quantitative loci (cis-eQTLs) derived from arterial wall and metabolic tissues obtained from patients with CHD identifies candidate genes at two potentially additional loci, ADAMTS9 and LOXL4. LD score regression reveals significant genetic correlations between cIMT and plaque traits, and both cIMT and plaque with CHD, any stroke subtype and ischemic stroke. Our study provides insights into genes and tissue-specific regulatory mechanisms linking atherosclerosis both to its functional genomic origins and its clinical consequences in humans.

Antidepressant drug-specific prediction of depression treatment outcomes from genetic and clinical variables

Individuals with depression differ substantially in their response to treatment with antidepressants. Specific predictors explain only a small proportion of these differences. To meaningfully predict who will respond to which antidepressant, it may be necessary to combine multiple biomarkers and clinical variables. Using statistical learning on common genetic variants and clinical information in a training sample of 280 individuals randomly allocated to 12-week treatment with antidepressants escitalopram or nortriptyline, we derived models to predict remission with each antidepressant drug. We tested the reproducibility of each prediction in a validation set of 150 participants not used in model derivation. An elastic net logistic model based on eleven genetic and six clinical variables predicted remission with escitalopram in the validation dataset with area under the curve 0.77 (95%CI; 0.66-0.88; p = 0.004), explaining approximately 30% of variance in who achieves remission. A model derived from 20 genetic variables predicted remission with nortriptyline in the validation dataset with an area under the curve 0.77 (95%CI; 0.65-0.90; p < 0.001), explaining approximately 36% of variance in who achieves remission. The predictive models were antidepressant drug-specific. Validated drug-specific predictions suggest that a relatively small number of genetic and clinical variables can help select treatment between escitalopram and nortriptyline.

Genome-wide analysis yields new loci associating with aortic valve stenosis

Aortic valve stenosis (AS) is the most common valvular heart disease, and valve replacement is the only definitive treatment. Here we report a large genome-wide association (GWA) study of 2,457 Icelandic AS cases and 349,342 controls with a follow-up in up to 4,850 cases and 451,731 controls of European ancestry. We identify two new AS loci, on chromosome 1p21 near PALMD (rs7543130; odds ratio (OR) = 1.20, P = 1.2 × 10⁻²²) and on chromosome 2q22 in TEX41 (rs1830321; OR = 1.15, P = 1.8 × 10⁻¹³). Rs7543130 also associates with bicuspid aortic valve (BAV) (OR = 1.28, P = 6.6 × 10⁻¹⁰) and aortic root diameter (P = 1.30 × 10⁻⁸), and rs1830321 associates with BAV (OR = 1.12, P = 5.3 × 10⁻³) and coronary artery disease (OR = 1.05, P = 9.3 × 10⁻⁵). The results implicate both cardiac developmental abnormalities and atherosclerosis-like processes in the pathogenesis of AS. We show that several pathways are shared by CAD and AS. Causal analysis suggests that the shared risk factors of Lp(a) and non-high-density lipoprotein cholesterol contribute substantially to the frequent co-occurence of these diseases.

Aortic valve stenosis (AS) is the most common valvular heart disease. Here the authors identify two new AS loci that also associate with bicuspid aortic valve, aortic root diameter and/or coronary artery disease implicating both developmental abnormalities and atherosclerosis-like processes in AS.

A Large-Scale Multi-ancestry Genome-wide Study Accounting for Smoking Behavior Identifies Multiple Significant Loci for Blood Pressure

Genome-wide association analysis advanced understanding of blood pressure (BP), a major risk factor for vascular conditions such as coronary heart disease and stroke. Accounting for smoking behavior may help identify BP loci and extend our knowledge of its genetic architecture. We performed genome-wide association meta-analyses of systolic and diastolic BP incorporating gene-smoking interactions in 610,091 individuals. Stage 1 analysis examined ∼18.8 million SNPs and small insertion/deletion variants in 129,913 individuals from four ancestries (European, African, Asian, and Hispanic) with follow-up analysis of promising variants in 480,178 additional individuals from five ancestries. We identified 15 loci that were genome-wide significant (p < 5 × 10-8) in stage 1 and formally replicated in stage 2. A combined stage 1 and 2 meta-analysis identified 66 additional genome-wide significant loci (13, 35, and 18 loci in European, African, and trans-ancestry, respectively). A total of 56 known BP loci were also identified by our results (p < 5 × 10-8). Of the newly identified loci, ten showed significant interaction with smoking status, but none of them were replicated in stage 2. Several loci were identified in African ancestry, highlighting the importance of genetic studies in diverse populations. The identified loci show strong evidence for regulatory features and support shared pathophysiology with cardiometabolic and addiction traits. They also highlight a role in BP regulation for biological candidates such as modulators of vascular structure and function (CDKN1B, BCAR1-CFDP1, PXDN, EEA1), ciliopathies (SDCCAG8, RPGRIP1L), telomere maintenance (TNKS, PINX1, AKTIP), and central dopaminergic signaling (MSRA, EBF2).

Integrative studies implicate matrix metalloproteinase-12 as a culprit gene for large-artery atherosclerotic stroke

BACKGROUND

Ischaemic stroke and coronary heart disease are important contributors to the global disease burden and share atherosclerosis as the main underlying cause. Recent evidence from a genome-wide association study (GWAS) suggested that single nucleotide polymorphisms (SNP) near the MMP12 gene at chromosome 11q22.3 were associated with large-vessel ischaemic stroke. Here, we evaluated and extended these results by examining the relationship between MMP12 and atherosclerosis in clinical and experimental studies.

METHODS AND RESULTS

Plasma concentrations of MMP12 were measured at baseline in 3394 subjects with high-risk for cardiovascular disease (CVD) using the Olink ProSeek CVD I array. The plasma MMP12 concentration showed association with incident cardiovascular and cerebrovascular events (130 and 67 events, respectively, over 36 months) and carotid intima-media thickness progression (P = 3.6 × 10^-5 ). A GWAS of plasma MMP12 concentrations revealed that SNPs rs499459, rs613084 and rs1892971 at chr11q22.3 were independently associated with plasma MMP12 (P < 5 × 10^-8 ). The lead SNPs showed associations with mRNA levels of MMP12 and adjacent MMPs in atherosclerotic plaques. MMP12 transcriptomic and proteomic levels were strongly significantly increased in carotid plaques compared with control arterial tissue and in plaques from symptomatic versus asymptomatic patients. By combining immunohistochemistry and proximity ligation assay, we demonstrated that MMP12 localizes to CD68 + macrophages and interacts with elastin in plaques. MMP12 silencing in human THP-1-derived macrophages resulted in reduced macrophage migration.

CONCLUSIONS

Our study supports the notion that MMP12 is implicated in large-artery atherosclerotic stroke, functionally by enhancing elastin degradation and macrophage invasion in plaques.

Large-scale genome-wide analysis identifies genetic variants associated with cardiac structure and function

BACKGROUND

Understanding the genetic architecture of cardiac structure and function may help to prevent and treat heart disease. This investigation sought to identify common genetic variations associated with inter-individual variability in cardiac structure and function.

METHODS

A GWAS meta-analysis of echocardiographic traits was performed, including 46,533 individuals from 30 studies (EchoGen consortium). The analysis included 16 traits of left ventricular (LV) structure, and systolic and diastolic function.

RESULTS

The discovery analysis included 21 cohorts for structural and systolic function traits (n = 32,212) and 17 cohorts for diastolic function traits (n = 21,852). Replication was performed in 5 cohorts (n = 14,321) and 6 cohorts (n = 16,308), respectively. Besides 5 previously reported loci, the combined meta-analysis identified 10 additional genome-wide significant SNPs: rs12541595 near MTSS1 and rs10774625 in ATXN2 for LV end-diastolic internal dimension; rs806322 near KCNRG, rs4765663 in CACNA1C, rs6702619 near PALMD, rs7127129 in TMEM16A, rs11207426 near FGGY, rs17608766 in GOSR2, and rs17696696 in CFDP1 for aortic root diameter; and rs12440869 in IQCH for Doppler transmitral A-wave peak velocity. Findings were in part validated in other cohorts and in GWAS of related disease traits. The genetic loci showed associations with putative signaling pathways, and with gene expression in whole blood, monocytes, and myocardial tissue.

CONCLUSION

The additional genetic loci identified in this large meta-analysis of cardiac structure and function provide insights into the underlying genetic architecture of cardiac structure and warrant follow-up in future functional studies.

FUNDING

For detailed information per study, see Acknowledgments.

Rare Variants in NOD1 Associated with Carotid Bifurcation Intima-Media Thickness in Dominican Republic Families

Cardiovascular disorders including ischemic stroke (IS) and myocardial infarction (MI) are heritable; however, few replicated loci have been identified. One strategy to identify loci influencing these complex disorders is to study subclinical phenotypes, such as carotid bifurcation intima-media thickness (bIMT). We have previously shown bIMT to be heritable and found evidence for linkage and association with common variants on chromosome 7p for bIMT. In this study, we aimed to characterize contributions of rare variants (RVs) in 7p to bIMT. To achieve this aim, we sequenced the 1 LOD unit down region on 7p in nine extended families from the Dominican Republic (DR) with strong evidence for linkage to bIMT. We then performed the family-based sequence kernel association test (famSKAT) on genes within the 7p region. Analyses were restricted to single nucleotide variants (SNVs) with population based minor allele frequency (MAF) <5%. We first analyzed all exonic RVs and then the subset of only non-synonymous RVs. There were 68 genes in our analyses. Nucleotide-binding oligomerization domain (NOD1) was the most significantly associated gene when analyzing exonic RVs (famSKAT p = 9.2x10^-4; number of SNVs = 14). We achieved suggestive replication of NOD1 in an independent sample of twelve extended families from the DR (p = 0.055). Our study provides suggestive statistical evidence for a role of rare variants in NOD1 in bIMT. Studies in mice have shown Nod1 to play a role in heart function and atherosclerosis, providing biologic plausibility for a role in bIMT thus making NOD1 an excellent bIMT candidate.

Genetic Research and Women's Heart Disease: a Primer

PURPOSE OF REVIEW

This review provides a brief synopsis of sexual dimorphism in atherosclerosis with an emphasis on genetic studies aimed to better understand the atherosclerotic process and clinical outcomes in women. Such studies are warranted because development of atherosclerosis, impact of several traditional risk factors, and burden of coronary heart disease (CHD) differ between women and men.

RECENT FINDINGS

While most candidate gene studies pool women and men and adjust for sex, some sex-specific studies provide evidence of association between candidate genes and prevalent and incident CHD in women. So far, most genome-wide association studies (GWAS) also failed to consider sex-specific associations. The few GWAS focused on women tended to have small sample sizes and insufficient power to reject the null hypothesis of no association even if associations exist. Few studies consider that sex can modify the effect of gene variants on CHD. Sufficiently large-scale genetic studies in women of different race/ethnic groups, taking into account possible gene-gene and gene-environment interactions as well as hormone-mediated epigenetic mechanisms, are needed. Using the same disease definition for women and men might not be appropriate. Accurate phenotyping and inclusion of relevant outcomes in women, together with targeting the entire spectrum of atherosclerosis, could help address the contribution of genes to sexual dimorphism in atherosclerosis. Discovered genetic loci should be taken forward for replication and functional studies to elucidate the plausible underlying biological mechanisms. A better understanding of the etiology of atherosclerosis in women would facilitate future prevention efforts and interventions.

Genetic loci on chromosome 5 are associated with circulating levels of interleukin-5 and eosinophil count in a European population with high risk for cardiovascular disease

IL-5 is a Th2 cytokine which activates eosinophils and is suggested to have an atheroprotective role. Genetic variants in the IL5 locus have been associated with increased risk of CAD and ischemic stroke. In this study we aimed to identify genetic variants associated with IL-5 concentrations and apply a Mendelian randomisation approach to assess IL-5 levels for causal effect on intima-media thickness in a European population at high risk of coronary artery disease. We analysed SNPs within robustly associated candidate loci for immune, inflammatory, metabolic and cardiovascular traits. We identified 2 genetic loci for IL-5 levels (chromosome 5, rs56183820, BETA=0.11, P=6.73E(-5) and chromosome 14, rs4902762, BETA=0.12, P=5.76E(-6)) and one for eosinophil count (rs72797327, BETA=-0.10, P=1.41E(-6)). Both chromosome 5 loci were in the vicinity of the IL5 gene, however the association with IL-5 levels failed to replicate in a meta-analysis of 2 independent cohorts (rs56183820, BETA=0.04, P=0.2763, I(2)=24, I(2)-P=0.2516). No significant associations were observed between SNPs associated with IL-5 levels or eosinophil count and IMT measures. Expression quantitative trait analyses indicate effects of the IL-5 and eosinophil-associated SNPs on RAD50 mRNA expression levels (rs12652920 (r2=0.93 with rs56183820) BETA=-0.10, P=8.64E(-6) and rs11739623 (r2=0.96 with rs72797327) BETA=-0.23, P=1.74E(-29), respectively). Our data do not support a role for IL-5 levels and eosinophil count in intima-media thickness, however SNPs associated with IL-5 and eosinophils might influence stability of the atherosclerotic plaque via modulation of RAD50 levels.

Intima-media thickness of the descending aorta in patients with bicuspid aortic valve

Objective

A bicuspid aortic valve (BAV) is associated with accelerated aortic valve disease (AVD) and abnormalities in aortic elasticity. We investigated the intima-media thickness of the descending aorta (AoIMT) in patients with AVD with or without an ascending aortic aneurysm (AscAA), in relation to BAV versus tricuspid aortic valve (TAV) phenotype, type of valve disease, cardiovascular risk factors, and single-nucleotide polymorphisms (SNPs) with a known association with carotid IMT.

Methods and results

368 patients (210 with BAV, 158 with TAV,); mean age 64 ± 13 years) were examined using transesophageal echocardiography (TEE) before valvular and/or aortic surgery. No patient had a coronary disease (CAD). The AoIMT was measured on short-axis TEE images of the descending aorta using a semi-automated edge-detection technique. AoIMT was univariately (P < 0.05) related to age, blood pressure, smoking, creatinine, highly sensitive C-reactive protein, HDL, valve hemodynamics and BAV. In the TAV subgroup it was also associated with the rs200991 SNP. Using multivariate regression analysis, age was the main determinant for AoIMT (P < 0.001), followed by male gender (P = 0.02), BAV was no longer a significant predictor of AoIMT. AoIMT was still related to the rs200991 SNP in TAV (P = 0.034), and to creatinine in BAV (P = 0.019), when other variables were accounted for.

Conclusions

Intima-media thickness of the descending aorta is not affected by aortic valve morphology (BAV/TAV); age is the main determinant of AoIMT. Genetic markers (SNPs) known to influence IMT in the carotid artery seem to correlate to IMT in the descending aorta only in patients with TAV.

Impact of carotid atherosclerosis loci on cardiovascular events

BACKGROUND

Genome-wide association studies (GWAS) have identified six single-nucleotide polymorphisms (SNPs) related to carotid intima media thickness (cIMT) or plaque. However, whether these loci relate to other vascular diseases and subsequent vascular events is unclear.

METHODS AND RESULTS

We tested six SNPs (rs4888378, rs11781551, rs445925, rs6601530, rs17398575 and rs1878406) for association with subclinical atherosclerotic measures (cIMT, plaque presence and ankle-brachial index), as well as ischemic stroke, abdominal aortic aneurysm, peripheral or coronary artery disease (CAD) in the Second Manifestations of ARTerial disease (SMART) cohort. Four SNPs were associated with cIMT and two with plaque (p < 0.05). One SNP was also significantly associated to CAD (rs1878406, OR = 1.24, 95% CI = 1.08-1.42, p = 2 × 10(-3)). A genetic risk score (GRS) based on the cIMT-related SNPs was associated to increased risk of cIMT itself (p = 1 × 10(-3)), but not to other secondary outcomes or vascular events during follow-up (p = 0.86).

CONCLUSIONS

In addition to replicating previously published associations for cIMT, we confirmed a nominally significant effect between the GRS and cIMT.

Sequencing of candidate genes in Dominican families implicates both rare exonic and common non-exonic variants for carotid intima-media thickness at bifurcation

Through linkage and tagSNP-based association studies in 100 Dominican Republic (DR) families, we previously identified ANLN and AOAH (7p14.3) as candidate genes for carotid intima-media thickness at bifurcation (bIMT). Introns, exons, and flanking regions of ANLN and AOAH were re-sequenced in 151 individuals from nine families with evidence for linkage at 7p14.3. For common variants [CV, minor allele frequency (MAF) ≥5 %], single variant-based analysis was performed. For rare variants (RV, MAF < 5 %), gene-based analysis aggregating all RVs within a gene was performed. CV analysis revealed the strongest signal at rs3815483 (P = 0.0003) in ANLN and rs60023210 (P = 0.00005) in AOAH. In ANLN, RV analysis found suggestive evidence for association with exonic RVs (P = 0.08), and in particular non-synonymous RVs (P = 0.04) but not with all RVs (P = 0.15). The variant alleles of all non-synonymous RVs segregated with the major allele of rs3815483 and were associated with lower bIMT while a novel synonymous RV segregated with the minor allele of rs3815483 and was associated with greater bIMT. Additional analysis in 561 DR individuals found suggestive evidence for association with all ANLN non-synonymous RVs (P = 0.08). In AOAH, no evidence for association with RVs was detected. Instead, conditional analysis revealed that multiple independent intronic CVs are associated with bIMT in addition to rs60023210. We demonstrate the utility of using family-based studies to evaluate the contribution of RVs. Our data suggest two modes of genetic architecture underlying the linkage and association at ANLN (multiple exonic RVs) and AOAH (multiple intronic CVs with uncharacterized functions).

Sex-Specific Effects of Adiponectin on Carotid Intima-Media Thickness and Incident Cardiovascular Disease

BACKGROUND

Plasma adiponectin levels have previously been inversely associated with carotid intima-media thickness (IMT), a marker of subclinical atherosclerosis. In this study, we used a sex-stratified Mendelian randomization approach to investigate whether adiponectin has a causal protective influence on IMT.

METHODS AND RESULTS

Baseline plasma adiponectin concentration was tested for association with baseline IMT, IMT progression over 30 months, and occurrence of cardiovascular events within 3 years in 3430 participants (women, n=1777; men, n=1653) with high cardiovascular risk but no prevalent disease. Plasma adiponectin levels were inversely associated with baseline mean bifurcation IMT after adjustment for established risk factors (β=-0.018, P<0.001) in men but not in women (β=-0.006, P=0.185; P for interaction=0.061). Adiponectin levels were inversely associated with progression of mean common carotid IMT in men (β=-0.0022, P=0.047), whereas no association was seen in women (0.0007, P=0.475; P for interaction=0.018). Moreover, we observed that adiponectin levels were inversely associated with coronary events in women (hazard ratio 0.57, 95% CI 0.37 to 0.87) but not in men (hazard ratio 0.82, 95% CI 0.54 to 1.25). A gene score of adiponectin-raising alleles in 6 loci, reported recently in a large multi-ethnic meta-analysis, was inversely associated with baseline mean bifurcation IMT in men (β=-0.0008, P=0.004) but not in women (β=-0.0003, P=0.522; P for interaction=0.007).

CONCLUSIONS

This report provides some evidence for adiponectin protecting against atherosclerosis, with effects being confined to men; however, compared with established cardiovascular risk factors, the effect of plasma adiponectin was modest. Further investigation involving mechanistic studies is warranted.

Functional Analysis of a Carotid Intima-Media Thickness Locus Implicates BCAR1 and Suggests a Causal Variant

BACKGROUND

Carotid intima-media thickness (IMT) is a marker of subclinical atherosclerosis that can predict cardiovascular disease events over traditional risk factors. This study examined the BCAR1-CFDP1-TMEM170A locus on chromosome 16, associated with carotid IMT and coronary artery disease in the IMT and IMT-Progression as Predictors of Vascular Events (IMPROVE) cohort, to identify the functional variant.

METHODS AND RESULTS

In analysis of the locus lead single nucleotide polymorphism (SNP; rs4888378, intronic in CFDP1) in Progressione della Lesione Intimale Carotidea (PLIC), the protective AA genotype was associated with slower IMT progression in women (P=0.04) but not in men. Meta-analysis of 5 cohort studies also supported a protective effect of the A allele on common carotid IMT in women only (women: β=-0.0047, P=1.63 × 10(-4); men: β=-0.0029, P=0.0678). Two hundred fourteen noncoding variants in strong linkage disequilibrium (r(2) ≥ 0.8) with rs4888378 were identified from 1000 Genome Project. ENCODE regulatory chromatin marks were used to create a shortlist of 6 possible regulatory variants. Electrophoretic mobility shift assays on the shortlist detected allele-specific protein binding to the lead SNP rs4888378; multiplexed competitor electrophoretic mobility shift assays implicated FOXA as the protein. Luciferase reporter assays on rs4888378 showed a significant 35% to 92% (P=0.0057; P=4.0 × 10(-22)) decrease in gene expression with the A allele. Expression quantitative trait loci analysis confirmed previously reported associations of rs4888378 with BCAR1 in vascular tissues.

CONCLUSIONS

Molecular studies suggest the lead SNP as a potentially causal SNP at the BCAR1-CFDP1-TMEM170A locus, and expression quantitative trait loci studies implicate BCAR1 as the causal gene. This variant showed stronger effects on common carotid IMT in women, raising questions about the mechanism of the causal SNP on atherosclerosis.

p130Cas/BCAR1 scaffold protein in tissue homeostasis and pathogenesis

BCAR1 (also known as p130Cas/BCAR1) is an adaptor protein that belongs to the CAS family of scaffold proteins. In the past years, increasing evidence has demonstrated the ability of p130Cas/BCAR1 to activate signaling originating from mechanical stimuli, cell-extracellular matrix (ECM) adhesion and growth factor stimulation cascades during normal development and disease in various biological models. In this review we will specifically discuss the more recent data on the contribution of p130Cas/BCAR1 in the regulation of tissue homeostasis and its potential implications in pathological conditions.

Human genetic evidence for involvement of CD137 in atherosclerosis

Atherosclerosis is an inflammatory disease and the main cause of cardiovascular disease. Inflammation promotes plaque instability and clinical disease, such as myocardial infarction, stroke and peripheral vascular disease. Subclinical atherosclerosis begins with thickening of the arterial intimal layer, and increased intima-media thickness (IMT) in the carotid artery is a widely used measurement of subclinical atherosclerosis. Activation of CD137 (tumor necrosis factor receptor super family 9) promotes inflammation and disease development in murine atherosclerosis. CD137 is expressed in human atherosclerosis, but its role is largely unknown. This study uses a genetic approach to investigate CD137 in human atherosclerotic disease. In publicly available data on genotype and gene expression from the HapMap project, the minor T allele of rs2453021, a single nucleotide polymorphism in CD137, was significantly associated with CD137 gene expression. In the PROCARDIS and Wellcome Trust Case Control Consortium (WTCCC) cohorts of 13,029 cases and controls, no significant association was detected between the minor T allele of rs2453021 and risk for coronary artery disease or myocardial infarction. However, in the IMPROVE multicenter study of 3,418 individuals, the minor T allele of rs2453021 was associated with increased IMT of the common carotid artery (CCA), as measured by ultrasonography, with presence of plaque in CCA and with increased incidence of adverse noncardiac vascular events. Taken together, this study shows that the minor T allele of rs2453021 is associated with increased IMT in the CCA and increased risk of incident noncardiac vascular events, thus providing the first human genetic evidence for involvement of CD137 in atherosclerosis.

Common genetic determinants of lung function, subclinical atherosclerosis and risk of coronary artery disease

Chronic obstructive pulmonary disease (COPD) independently associates with an increased risk of coronary artery disease (CAD), but it has not been fully investigated whether this co-morbidity involves shared pathophysiological mechanisms. To identify potential common pathways across the two diseases, we tested all recently published single nucleotide polymorphisms (SNPs) associated with human lung function (spirometry) for association with carotid intima-media thickness (cIMT) in 3,378 subjects with multiple CAD risk factors, and for association with CAD in a case-control study of 5,775 CAD cases and 7,265 controls. SNPs rs2865531, located in the CFDP1 gene, and rs9978142, located in the KCNE2 gene, were significantly associated with CAD. In addition, SNP rs9978142 and SNP rs3995090 located in the HTR4 gene, were associated with average and maximal cIMT measures. Genetic risk scores combining the most robustly spirometry–associated SNPs from the literature were modestly associated with CAD, (odds ratio (OR) (95% confidence interval (CI₉₅) = 1.06 (1.03, 1.09); P-value = 1.5×10⁻⁴, per allele). In conclusion, our study suggests that some genetic loci implicated in determining human lung function also influence cIMT and susceptibility to CAD. The present results should help elucidate the molecular underpinnings of the co-morbidity observed across COPD and CAD.

CAS proteins in health and disease: an update

The CAS family of scaffolding proteins has increasingly attracted scrutiny as important for regulation of cancer-associated signaling. BCAR1 (also known as p130Cas), NEDD9 (HEF1, Cas-L), EFS (Sin), and CASS4 (HEPL) are regulated by and mediate cell attachment, growth factor, and chemokine signaling. Altered expression and activity of CAS proteins are now known to promote metastasis and drug resistance in cancer, influence normal development, and contribute to the pathogenesis of heart and pulmonary disease. In this article, we provide an update on recently published studies describing signals regulating and regulated by CAS proteins, and evidence for biological activity of CAS proteins in normal development, cancer, and other pathological conditions.

Genome-wide interaction study identifies RCBTB1 as a modifier for smoking effect on carotid intima-media thickness

OBJECTIVE

Carotid intima-media thickness (cIMT), a marker for atherosclerosis, is affected by smoking and has substantial interindividual variation. We sought to identify the genetic moderators influencing the effect of smoking on cIMT.

APPROACH AND RESULTS

With a multistage design using 722 379 single nucleotide polymorphisms (SNP), a genome-wide interaction study was performed in a discovery sample of 669 Hispanics, followed by replication in 589 subjects (264 Hispanics, 172 non-Hispanic blacks, 153 non-Hispanic whites). Assuming an additive genetic model, regression analysis was performed to test for smoking-SNP interaction on cIMT while controlling for age, sex, and the top 3 principal components of ancestry. The strongest interaction in Hispanics was found with a synonymous splicing SNP (rs3751383) in exon 9 of RCBTB1 (P=2.5e(-6) in discovery sample; P=0.01 in the Hispanic replication sample; P<8.8e(-9) in the combined Hispanic sample). Stratification analysis in the combined Hispanic sample showed that smoking had no effect on cIMT among rs3751383 G homozygote (P=0.15), a moderate effect among rs3751383 heterozygote (P=0.01), and a strong effect among rs3751383 A homozygote (P=2.1e(-7)). A consistent trend was observed in the non-Hispanic white and black data sets, leading to an interaction effect of P<2.9e(-9) in the meta-analysis of all 1258 subjects.

CONCLUSIONS

Our study represents the first genome-wide smoking-SNP interaction study of cIMT and identifies RCBTB1 as a modifier of the smoking effect on cIMT. Testing for gene-environment interactions can help uncover genetic factors that contribute to the interindividual variation in response to the same environmental exposure.