ORIGINAL RESEARCH article

Front. Psychiatry, 26 September 2023

Sec. Behavioral and Psychiatric Genetics

Volume 14 - 2023 | https://doi.org/10.3389/fpsyt.2023.1279962

Determination of genetic correlation between tobacco smoking and coronary artery disease

  • 1. Joint Institute of Tobacco and Health, Kunming, Yunnan, China

  • 2. State Key Laboratory for Diagnosis and Treatment of Infectious Diseases, The First Affiliated Hospital, Collaborative Innovation Center for Diagnosis and Treatment of Infectious Diseases, Zhejiang University School of Medicine, Hangzhou, Zhejiang, China

Article metrics

View details

5

Citations

3k

Views

1,4k

Downloads

Abstract

Backgrounds:

Tobacco smoking is an important risk factor for coronary artery disease (CAD), but the genetic mechanisms linking smoking to CAD remain largely unknown.

Methods:

We analyzed summary data from the genome-wide association study (GWAS) of the UK Biobank for CAD, plasma lipid concentrations (n = 184,305), and smoking (n = 337,030) using different biostatistical methods, which included LD score regression and Mendelian randomization (MR).

Results:

We identified SNPs shared by CAD and at least one smoking behavior, the genes where these SNPs are located were found to be significantly enriched in the processes related to lipoprotein metabolic, chylomicron-mediated lipid transport, lipid digestion, mobilization, and transport. The MR analysis revealed a positive correlation between smoking cessation and decreased risk for CAD when smoking cessation was considered as exposure (p = 0.001), and a negative correlation between the increased risk for CAD and smoking cessation when CAD was considered as exposure (p = 2.95E-08). This analysis further indicated that genetic liability for smoking cessation increased the risk of CAD.

Conclusion:

These findings inform the concomitant conditions of CAD and smoking and support the idea that genetic liabilities for smoking behaviors are strongly associated with the risk of CAD.

Introduction

Tobacco smoking is one of the most important public health problems world-wide, accounting for 9% of deaths (1). Several epidemiologic studies have proved that tobacco smoking is a major risk factor for many diseases such as lung cancer and cardiovascular and respiratory diseases (1). Lung cancer is the most prevalent smoking-associated cause of death, followed by ischemic heart disease and chronic airway obstruction. Tobacco smoking and its impact on the respiratory system caused an estimated 8 million deaths per year, with more than 10% of these deaths related to second-hand smoke (2). In the United States, tobacco smoking is associated with 30% of all CAD-related deaths each year (3) and also doubles the risk of premature cardiovascular deaths (4).

Several epidemiologic studies have revealed that tobacco smoking increases the incidence of fatal CAD and associates with various cardiovascular diseases (5–8). Extensive clinical evidence has supported the idea that tobacco smoking causes multiple genetic and epigenetic abnormalities in the respiratory epithelium (9, 10). In a previous study, Sabater-Ileal and colleagues identified a genetic locus that influences both lung function and CAD (11), although the findings were not genome-wide in scale and were underpowered due to a small sample size.

Tobacco contains more than 4,000 chemicals (12), and the exact toxic components and the mechanisms involved in tobacco-related CAD and cardiovascular dysfunction are still unknown. Recently, a genetic predisposition to the development of atherogenesis in individuals exposed to cigarette smoke has been reported. The commonly documented examples are CYP1A1 MSP polymorphism and certain endothelial NO synthase intron 4 polymorphisms. Both increase the susceptibility to cigarette smoke exposure-related atherosclerotic diseases including multi-vessel CAD and myocardial infarction (MI) (13, 14). Given that much of the available data were derived from observational studies, which are unable to account fully for confounding and reverse causation, the genetic correlation and causal relations between smoking and cardiovascular diseases remain to be determined.

The principle of Mendelian randomization (MR) relies on the basic laws of Mendelian genetics, segregation, and independent association. When these principles hold at a population level, the influence of confounding factors can be reduced because religion, growth, environment, and other confounding factors are considered to be random (15). Given that alleles are randomly allocated and become fixed at conception, MR studies are less susceptible to reverse causality than are observational studies.

In this study, we examined the pleiotropic effect of tobacco smoking and CAD using publicly available GWAS summary statistics. Then, we used bidirectional MR method to reveal the nature of the causal relations between CAD and tobacco smoking. Finally, we determined the biological processes or pathways involved in the comorbidity of these two diseases.

Materials and methods

GWAS summary data sets

The GWAS summary data for CAD and plasma lipid concentrations have been described in a previous report (16). Briefly, the summary statistics of a large GWAS meta-analysis comprising more than 120,000 CAD cases and 339,115 controls were obtained from CARDIoGRAMplusC4D Consortium website (http://www.cardiogramplusc4d.org/data-downloads/). A total of 9,149,595 variants were included either in the CARDIoGRAMplusC4D 1,000 Genomes–imputed GWAS or the MIGen/CARDIoGRAM Exome chip study. The smoking data were obtained from the large GSCAN summary statistics (17). The GSCAN investigated four smoking-related phenotypes, including age at initiation of regular smoking (AgeSmk; n = 341,427), whether an individual had ever smoked regularly (SmkInit; n = 1,232,091), cigarettes smoked per day (CPD; n = 337,334), and smoking cessation (SmkCes; n = 547,219). The GWAS summary statistics for different smoking phenotypes can be found at https://conservancy.umn.edu/handle/11299/201564.

Further, we obtained published GWAS meta-analysis association data for lipid concentrations from the Center for Statistical Genetics, which was a joint analysis that examined 188,577 individuals whose genomic DNA samples were genotyped with two platforms from multiple studies (18). Complete GWAS summary statistics were downloaded from webpage http://csg.sph.umich.edu/willer/public/lipids2013/.

Estimation of genetic correlation by LD score regression (LDSC)

The genetic correlations (rg) between CAD and smoking behaviors were estimated by LDSC (19). Pairwise LD r2 among SNPs was conducted using pre-computed LD scores with the 1,000 Genomes Project reference panel of subjects of European ancestry. Quality control steps were adopted from LD scores default procedures, including imputation quality >0.9 and minor allele frequency > 0.1. Moreover, all SNPs retained for further analysis were merged with SNPs in the HapMap 3 reference panel. Correlation was considered significant at a corrected p value of <0.05 by Bonferroni correction.

Mendelian randomization (MR)

We extracted the effect estimates and standard errors (Ses) from relevant GWAS and employed TwoSampleMR (v. 0.4.22) R package to clarify the potential causal effect for both smoking and CAD (20). The following strategies were used to identify genetic instruments. First, we filtered GWAS summary datasets to require shared susceptible loci in both smoking and CAD. The variants showing genome-wide significance (p < 5 × 10−8) in GWAS for CAD were considered to be candidate variants, and then we checked the significance of these genetic loci separately in four other smoking behaviors: smoking initiation age, smoking cessation, CPD, and age of initiation in GSCAN studies. The common SNPs were harmonized using default parameters within the built-in “harmonize data” function and then trimmed by PLINK (v. 1.07) to obtain independent risk variants for each disease (21).

To start, the MR analysis was performed by generating instrumental variable estimates for each SNP. The averaged causal estimate of each SNP was calculated using the inverse–variance-weighted (IVW) method, i.e., specifically defined as the beta coefficient associated with SNP-CAD divided by the beta coefficient associated with SNP-smoking behaviors (22). In addition, we used a series of sensitivity analyses, which included weighted median and MR Egger, to evaluate the reliability of our results.

Gene and pathway analysis

The gene-based analysis that links SNPs to genes was conducted using MAGMA with default settings. To gain biological insights into shared genes, we used the WebGestalt tool (23) to assess enrichment of the identified shared gene set in the Gene Ontology (GO) biological processes with redundant GO terms been removed. Both analyses were based on shared genes that were identified from cross-trait meta-analysis. Pathways with a false discovery rate (FDR) < 0.05 were considered significant.

Results

Susceptible loci shared by CAD and smoking behaviors

To investigate the genetic overlap between CAD and smoking behaviors, we used GWAS summary data from large-scale genome-wide studies (Supplementary Table S1). We detected a great number of significantly associated SNPs overlapped between CAD risk loci and at least one smoking phenotype (Supplementary Table S2). Of the CAD GWAS loci, 2091 SNPs (35.01%) showed nominal significance in CPD, 526 SNPs (25.16%) showed nominal significance for smoking initiation, 317 SNPs (15.16%) for smoking cession, and 85 SNPs (4.07%) for age at smoking initiation. Notably, 24 SNPs reached genome-wide significance (Figure 1).

Figure 1

Genetic correlations between CAD and smoking behaviors

We used the LDSC method to test for the genome-wide correlations between CAD and smoking behaviors. Significant genetic correlations were found between CAD and all the smoking-related traits with the smallest p values <1 × 10−18 (Figure 2). We observed significant positive genetic correlations between CAD and CPD, smoking initiation, and smoking cessation (rg > 0.2), but a negative genetic correlation between CAD and age at initiation (rg < −0.2).

Figure 2

Mendelian randomization analysis

Considering the presence of potential LD relations among those significant SNPs of interest, we performed p value-informed LD pruning with the goal of obtaining independent GWAS SNPs. This led to the identification of 15, 63, 18, and 152 independent SNPs for Age of Initiation, CPD, Smoking Cessation and Smoking Initiation, respectively. Bi-directional MR analysis provided strong evidence that smoking initiation increased the risk of CAD (IVW: β = 0.191; p = 2.59 × 10−6) with a consistent direction of effect in all three MR methods (Table 1). There also was evidence for a consistent but weaker genetic liability for smoking cessation on CAD (IVW: β = 0.234; p = 0.001). The same findings were observed for age at initiation as the instrument on CAD (IVW: β = 0.295; p = 0.039).

Table 1

ExposureOutcomeNo. SNPsBetap-valueMethod
Age of initiationCAD150.2950.039Inverse-variance weighted
150.2950.026Weighted median
15−0.4430.609MR Egger
CPDCAD63−0.0390.406Inverse-variance weighted
63−0.0720.097Weighted median
63−0.1750.025MR Egger
Smoking cessationCAD180.2340.001Inverse-variance weighted
180.2340.001Weighted median
180.0910.559MR Egger
Smoking initiationCAD1520.1912.59E-06Inverse-variance weighted
1520.1933.35E-06Weighted median
1520.4190.052MR Egger

Effect of smoking behaviors on CAD using two-sample MR analysis.

On the other hand, we obtained negative correlations between CPD and CAD, although the statistic was less significant. Only the result of MR Egger was significant. When treating CAD as an instrument, strong evidence of decreased risk of smoking cessation by CAD was observed (IVW: β = −0.046; p = 2.95 × 10−8; Table 2). The effects of CAD on other smoking behaviors became even weaker and inconsistent across the different methods.

Table 2

ExposureOutcomeNo. SNPsBetap-valueMethod
Age of initiation680.0120.013Inverse-variance weighted
CAD680.0120.03Weighted median
68−0.0010.932MR Egger
68−0.0470.075Inverse-variance weighted
CADCPD68−0.0630.018Weighted median
68−0.0550.366MR Egger
Smoking cessation75−0.0462.95E-08Inverse-variance weighted
CAD75−0.0464.29E-06Weighted median
75−0.0796.10E-04MR Egger
Smoking initiation680.0040.584Inverse-variance weighted
CAD680.0030.637Weighted median
68−0.0190.299MR Egger

Effect of CAD on smoking behaviors using two-sample MR analysis.

Biological pathway and enrichment analysis

We performed pathway analyses to identify biological pathways enriched for shared genetic loci related to smoking and CAD based on significant cross-trait meta-analysis results. For a detailed list of the overlapped genes and SNPs between CAD and smoking-related phenotypes, please refer to Table 3. Pathway analysis showed that the SNP-related genes were significantly enriched in lipoprotein metabolic, chylomicron-mediated lipid transport, lipid digestion, mobilization, and transport (Table 4). The GO analysis suggested that shared genes in CAD and smoking behaviors were enriched in triglyceride-rich lipoprotein particle clearance, blood vessel development, and very-low-density lipoprotein particle clearance (Table 5).

Table 3

GeneGene functionSNPs
ABI2ABL Interactor 2rs116426890
ADAMTS7ADAM metallopeptidase with Thrombospondin type 1 motif 7rs7182642; rs4887096; rs55834964; rs7174367; rs12916326; rs12916648; rs1809419; rs1809409; rs11635870; rs11635931; rs1807006; rs12286; rs1807007; rs8043119; rs11633351; rs6495267; rs1809420; rs4886591; rs2904223; rs4887102; rs35934157
ALDH2Aldehyde dehydrogenase 2 family memberrs11513729
APOBApolipoprotein Brs503105
APOC1Apolipoprotein C1rs12721051; rs56131196; rs4420638; rs12721046; rs111789331; rs66626994; rs10414043; rs7256200
APOC1P1Apolipoprotein C1 pseudogene 1rs111789331; rs66626994
APOEApolipoprotein Ers429358; rs10414043; rs7256200; rs769449
ATXN2Ataxin 2rs4766578;rs10774625;rs653178;rs597808;rs7137828;rs11065979; rs11065987
BRAPBRCA1 associated proteinrs11065979; rs11065987; rs11065991
C19orf38Chromosome 19 open reading frame 38rs11881156; rs12973042; rs12979495; rs67987899; rs11881438; rs7260254; rs35443547
CARFCalcium responsive transcription factorrs140244541; rs72928613; rs72936866; rs72932559; rs114702158; rs72928620; rs72932566; rs72928610; rs114123510; rs72928609; rs72932560; rs72926798; rs141002954; rs72932554;rs72932561; rs9749722; rs74421437; rs72936872; rs72936873
CARM1Coactivator associated arginine methyltransferase 1rs12971616; rs4804547; rs17616105; rs12976693; rs35991287; rs74179956; rs12977506; rs73009538; rs4804142; rs11880628; rs35734575; rs11879571; rs36045835; rs7256879; rs34857893; rs2053064; rs1971038; rs8105092; rs2053065; rs1541595; rs1971039; rs2304089
CDKN2B-AS1CDKN2B antisense RNA 1rs568447; rs545226; rs1333045; rs10217586; rs10738604; rs10811641; rs3217992; rs3218020; rs2891168; rs4977574; rs10757274; rs504318; rs10757264; rs7028268; rs567453; rs496892; rs490005; rs10757272; rs10738606; rs10120688; rs10738607; rs10811645; rs10811643; rs10738610; rs7035484; rs7049105; rs1591136; rs10811647; rs3218012; rs10115049; rs10811656; rs10757266; rs10757279; rs10811650
CELSR2Cadherin EGF LAG seven-pass G-type receptor 2rs11102967; rs611917
CHRNB4Cholinergic receptor nicotinic beta 4 subunitrs56354501; rs8031513; rs56195905; rs1809412; rs1814880; rs11072799; rs1809414; rs11856441; rs12901228; rs2219939; rs1825086; rs1383635; rs1383634; rs58717592; rs59671175; rs1809415
CNNM2Cyclin and CBS domain divalent metal cation transport mediator 2rs12414777
COL4A1Collagen type IV alpha 1 chainrs12866570
COL4A2Collagen type IV alpha 2 chainrs55940034; rs9515201
CYP20A1Cytochrome P450 family 20 subfamily A member 1rs72938351; rs72938315; rs115600411; rs114863726; rs114407963; rs116443099
DMRTA1DMRT like family A1rs10738610; rs10811656; rs10757279; rs10757277; rs1333046; rs1333049; rs10757278; rs1333048; rs4977575; rs7857118; rs1333047
DNM2Dynamin 2rs117159625; rs11880613; rs117786851; rs8111962; rs1109375; rs2287029
EDNRAEndothelin receptor type Ars7668383
FAM117BFamily with sequence similarity 117 member Brs146289328
FESFES proto-oncogene, tyrosine kinasers2521501
GGCXGamma-glutamyl carboxylasers7568458; rs35215812; rs10198569; rs10187424; rs12473819; rs6547624; rs6757263
HECTD4HECT domain E3 ubiquitin protein ligase 4rs11066188; rs7953257
ICA1LIslet cell autoantigen 1 likers11675462; rs4675290; rs2036927; rs12693975; rs10207567; rs1971819; rs934287; rs72932720; rs72932716; rs72932711; rs72932709; rs72932707; rs76122535; rs72932722; rs72932723; rs72932731; rs72932737; rs72934734; rs72932725; rs143911965; rs6435168; rs6732078; rs72932746; rs80087860; rs6705330; rs72932752; rs78128841
IL6RInterleukin 6 receptorrs6694258; rs6689393; rs4845625; rs12129500; rs6686750; rs4845619; rs7549338; rs6667434; rs7553796; rs4553185; rs59632925; rs7549250; rs4845371; rs11265612; rs6687726; rs4845618; rs6694817; rs11265611; rs6427658; rs12118721; rs10908838; rs12118770; rs12117832; rs6689306
KCNE2Potassium voltage-gated channel subfamily E regulatory subunit 2rs1018757; rs2211693; rs2211695; rs2211694; rs8128536; rs7279974; rs8129119; rs4817636; rs2096469; rs4817637; rs2096467; rs2834439; rs7280276; rs9982672; rs1023354; rs28451064; rs762158; rs9980618; rs743339; rs7277800; rs9982601; rs4817639; rs9974878; rs7281592; rs60687229; rs2834431; rs28591415; rs16991453; rs9978142; rs7283231; rs8132042; rs8131284
KCNK5Potassium two pore domain channel subfamily K member 5rs1155347; rs10456468; rs56015508; rs55902013; rs9394577; rs12211281; rs55856036; rs1544935
LDLRLow density lipoprotein receptorrs8103309; rs11666925; rs876794; rs10409001; rs10415811; rs11085757; rs8102273; rs11879026; rs73013202; rs73013198
LINC00310Long intergenic non-protein coding RNA 310rs1018757; rs2211693; rs2211695; rs2211694; rs8128536; rs7279974; rs8129119; rs4817636; rs2096469; rs4817637; rs2096467; rs2834439; rs7280276; rs9982672; rs1023354; rs28451064
LINC00841Long intergenic non-protein coding RNA 841rs2051120; rs6593388; rs2624688; rs6593393; rs4948800; rs2085797; rs7478408; rs898549; rs10899956; rs12359058; rs10899955; rs7924201; rs7088951; rs10899963; rs10899954; rs898551; rs10793515; rs66887775
LOC646938TBC1 domain family member 2B pseudogeners56354501; rs8031513; rs55834964; rs56195905; rs1809412; rs1814880; rs11072799; rs1809414; rs11856441; rs12901228; rs2219939; rs1825086
MAP3K4Mitogen-activated protein kinase kinase kinase 4rs117340856
MAPKAPK5-AS1MAPKAPK5 antisense RNA 1rs11513729
MC4RMelanocortin 4 receptorrs523288; rs538656; rs571312; rs663129; rs6567160; rs34633411; rs11152213; rs66723169; rs1942872; rs12958167; rs12955983; rs12954782; rs11663816
MIA3MIA SH3 domain ER export factor 3rs1909196; rs34679168; rs34767248; rs2378584; rs17163301; rs4618978; rs4846769; rs4846384; rs4846770; rs17163345; rs17163313
MIR4422Microrna 4,422rs55694910; rs72664304; rs72664303
MORF4L1Mortality factor 4 like 1rs11634042; rs8034274; rs7168915; rs11852830; rs11857877; rs4439728; rs12232282; rs4420501; rs4438276; rs8037171; rs11637783
MRASMuscle RAS oncogene homologrs13324341; rs1199337; rs1199338; rs185244
MTAPMethylthioadenosine phosphorylasers7041637
NAA25N-alpha-acetyltransferase 25, NatB auxiliary subunitrs17696736
NBEAL1Neurobeachin like 1rs72932566; rs76461893; rs145299755; rs72932573; rs148707292; rs72932572; rs72932574; rs72932583; rs6728861; rs72932575; rs4675310; rs148812085; rs2351524; rs140750546; rs72934573; rs140168762; rs72934505; rs145538381; rs115654617
NOS3Nitric oxide synthase 3rs3918226
PHACTR1Phosphatase and actin regulator 1rs10807323; rs13209002; rs9381401; rs12202891; rs9395224; rs9472790; rs1571997; rs9381462; rs6916397; rs2026458; rs6916421; rs7776079
PLGPlasminogenrs117340856
PMAIP1Phorbol-12-myristate-13-acetate-induced protein 1rs523288; rs538656; rs571312; rs663129; rs6567160; rs34633411; rs11152213; rs66723169; rs1942872; rs12958167; rs12955983; rs12954782; rs11663816; rs11664883; rs2045438; rs35476226; rs17175602; rs2045439; rs11660069; rs12957325; rs12970134
RAPH1RAS association (RalGDS/AF-6) and pleckstrin homology domains 1rs139644567; rs140274075
SAYSD1SAYSVFN motif domain containing 1rs1155347; rs10456468; rs56015508; rs55902013; rs9394577; rs12211281; rs55856036; rs1544935
SH2B3SH2B adaptor protein 3rs7310615; rs3184504; rs10774624
SH2D6SH2 domain containing 6rs6750847; rs7591175; rs6722691; rs11895399; rs11895401; rs1446669; rs2166529; rs6547620; rs6721924; rs6739015; rs2044474; rs13394343; rs6719046; rs6733913; rs17026396
SMARCA4SWI/SNF related, matrix associated, actin dependent regulator of chromatin, subfamily a, member 4rs8103309; rs60314748; rs60448955; rs73013159; rs7275; rs55677033; rs3786727; rs11666925; rs11670205; rs3786725; rs10417578; rs68010235; rs6511718; rs876794
SPECC1L-ADORA2ASPECC1L-ADORA2A readthrough (NMD candidate)rs5760347; rs62233136; rs5760359; rs5751841; rs5760368; rs62233133; rs5760350; rs2298379
SVILSupervillinrs1774241; rs1774240; rs1832864; rs12259037; rs12779954; rs1418276; rs4749520; rs11007851; rs10826749
TAF1A-AS1TAF1A antisense RNA 1rs1909196; rs34679168; rs4618978
TMED1Transmembrane p24 trafficking protein 1rs73007593; rs11881156
TRAFD1TRAF-type zinc finger domain containing 1rs17630235
TTC29Tetratricopeptide repeat domain 29rs7668383
VAMP8Vesicle associated membrane protein 8rs35215812; rs10198569; rs10187424; rs6547624; rs6757263
WDR12WD repeat domain 12rs72934734; rs35212307; rs114395475; rs7582720; rs72936856; rs72934751; rs72934763; rs72934749; rs6738618; rs150788469; rs6725887; rs6435169; rs7560547; rs77931721; rs142603618; rs72936846; rs72936852; rs72936830; rs143035655; rs77268589
YIPF2Yip1 domain family member 2rs17850995
ZC3HC1Zinc finger C3HC-type Containing 1rs11556924; rs56179563

The overlapped genes and SNPs between CAD and smoking-related phenotypes.

Table 4

Pathway IDPathway namep-valueq-value (FDR B&H)Genes included
PW:0000482Lipoprotein metabolic5.15E-069.66E-04APOB, APOC1, APOE
1,270,005Lipoprotein metabolism4.57E-053.86E-03APOB, APOC1, APOE, LDLR
1,270,006Chylomicron-mediated lipid transport5.74E-054.22E-03APOB, APOE, LDLR
1,270,002Lipid digestion, mobilization, and transport2.54E-041.18E-02APOB, APOC1, APOE, LDLR
M12950Angiotensin-converting enzyme 2 regulates heart function5.49E-041.62E-02COL4A1, COL4A2
1,270,115Metabolism of nitric oxide1.60E-033.41E-02NOS3, DNM2
1,427,851VLDL interactions1.75E-033.41E-02APOB, APOC1
172,847Protein digestion and absorption1.83E-033.41E-02KCNK5, COL4A1, COL4A2
1,270,008LDL-mediated lipid transport1.91E-033.41E-02APOB, LDLR

Detected shared pathway between smoking and CAD based on pathway analysis.

Table 5

Gene Ontology (GO) IDGene Ontology (GO) termp-valueq-value
(FDR B&H)
Genes included
GO:0071830Triglyceride-rich lipoprotein particle clearance5.75E-097.17E-06APOB, APOC1, APOE, LDLR
GO:0001568Blood vessel development1.94E-061.34E-03APOB, APOE, NOS3, PLG, COL4A1, COL4A2, DNM2, EDNRA, MIA3, SMARCA4, IL6R, LDLR
GO:0034447Very-low-density lipoprotein particle clearance2.16E-061.34E-03APOB, APOC1, APOE
GO:0001944Vasculature development2.95E-061.43E-03APOB, APOE, NOS3, PLG, COL4A1, COL4A2, DNM2, EDNRA, MIA3, SMARCA4, IL6R, LDLR
GO:0072358Cardiovascular system development3.43E-061.43E-03APOB, APOE, NOS3, PLG, COL4A1, COL4A2, DNM2, EDNRA, MIA3, SMARCA4, IL6R, LDLR
GO:0042159Lipoprotein catabolic process8.10E-062.53E-03APOB, APOE, LDLR
GO:0030195Negative regulation of blood coagulation1.53E-053.72E-03APOE, NOS3, PLG, SH2B3
GO:0048514Blood vessel morphogenesis2.76E-054.31E-03APOB, APOE, NOS3, PLG, COL4A1, COL4A2, EDNRA, MIA3, SMARCA4, LDLR
GO:1990777Lipoprotein particle2.51E-062.30E-04APOB, APOC1, APOE, LDLR
GO:0034358Plasma lipoprotein particle2.51E-062.30E-04APOB, APOC1, APOE, LDLR
GO:0032994Protein-lipid complex3.09E-062.30E-04APOB, APOC1, APOE, LDLR
GO:0034362Low-density lipoprotein particle6.77E-063.00E-04APOB, APOE, LDLR
GO:0034364High-density lipoprotein particle4.26E-051.13E-03APOB, APOC1, APOE

Detected shared pathways between CAD and smoking based on GO analysis.

Discussion

In this study, we revealed the genetic correlation and causal relations between smoking and CAD, providing a comprehensive evaluation of the shared genetic etiology of tobacco smoking and cardiovascular diseases. Our findings have highlighted the discovery that different smoking behaviors have strong associations with CAD, specifically, the correlation between smoking initiation, smoking cessation, and CAD.

The approach to MR is based on the assumptions that: (1) the genetic marker is associated with the exposure; (2) the genetic marker is independent of any confounding factors; and (3) there is no association between the genetic marker and outcome except through confounding factors. However, it should be acknowledged that these assumptions generally are not easy to evaluate. Results from the present MR study were based on data from the GWAS, which has corroborated the results obtained from conventional prospective observational studies that confirmed that tobacco smoking is a risk factor for CAD (5–8).

To our knowledge, this study represents one of a few large-scale genome-wide analysis to investigate the genetic overlap between smoking and CAD (24–26). Similar to the findings from these reports, our analyses also revealed strong associations between smoking initiation, smoking cessation, and CAD. Further, we found a significant positive association between smoking initiation and CAD when smoking initiation was considered as exposure (inverse–variance-weighted: β = 0.191; p = 2.59E-06; weighted median: β = 0.193; p = 3.35E-06), suggesting that smokers are more susceptible to CAD. In addition, we found a negative correlation between CAD and smoking cessation when CAD was considered as exposure (inverse–variance-weighted: β = −0.046; p = 2.95E-08; weighted median: β = 0.193; p = 3.35E-06). This indicates that patients with CAD are less likely to quit smoking, possibly because of tobacco addiction. Together, these findings demonstrate the presence of shared genetic etiologies between tobacco smoking and CAD.

Moreover, we found strong evidence of a genetic correlation between CAD and serum high-density lipoprotein (HDL) and low-density lipoprotein (LDL) as well, which is consistent with the results reported previously (27–29). We observed a positive correlation between LDL and CAD when LDL was considered as the exposure (inverse–variance-weighted: β = 0.387; p = 4.96E-43; weighted median: β = 0.406; p = 7.41E-39), and a negative correlation between HDL and CAD when HDL was considered as the exposure (inverse–variance-weighted: β = −0.245; p = 1.19E-14; weighted median: β = −0.188; p = 2.66E-10). As is well documented in the literature, LDL is a strong risk factor for CAD (27–29) whereas HDL is an anti-atherosclerotic plasma lipoprotein and a protective factor against CAD (30–32).

Further, we found that smoking has significant associations with HDL and LDL. A positive correlation between CPD and LDL was observed when CPD was considered as exposure (inverse–variance-weighted: β =0.06; p = 0.01; weighted median: β = 0.06; p = 0.02), and a negative correlation between CPD and HDL when CPD was considered as exposure (inverse–variance-weighted: β = −0.06; p = 0.005; weighted median: β = −0.06; p = 0.02). There also is a positive correlation between LDL and smoking cessation when LDL was considered as exposure (inverse–variance-weighted: β = −0.03; p = 0.002; weighted median: β = −0.04; p = 2.43E-5). Taken together, these findings provided a clear indication that smoking increases the risk of CAD by affecting the regulation of LDL and HDL, which needs to be further investigated.

We also performed GO enrichment and KEGG pathway analyses based on the genes where the SNPs overlapped between CAD and smoking-related phenotypes are located. We found several functions and pathways to be related to the lipoprotein metabolic and blood vessel development, which are all closely associated with CAD. It has been reported that the APOE–APOC1–APOC2–APOC4 cluster was significantly related to lipoprotein-associated phospholipase A2 mass and activity and CAD (33). Interestingly, the SNP rs4420638 located downstream of the APOC1 gene was found to be significantly related to smoking cessation (p = 7.4E-6) (34). Moreover, as a brain eQTL based on the information from BRAINEAC, rs4420638, this SNP has been linked to Alzheimer’s disease (35–40) and cognitive decline (41).

We used non-overlapping data sources in the context of summary-level MR analysis of exposure and outcome, which greatly improved the confidence in the causal effect estimates. In addition, through a range of sensitivity analysis methods, similar causal estimates and consistent causal inferences could be drawn. However, this study had limitations as well. First, there was a stark difference in sample sizes among different phenotypes, which might contribute to discrepancies in statistical power. Second, the information available on the summary-level GWAS data had limited us to divide samples into subgroups, which prevented us from studying the age-related heterogeneity.

In conclusion, this was a systematic analysis of the shared etiology and possible causal relations of smoking and CAD by employing large-scale GWASs. Genetic methods represent another option for assessing causality when there are too many confounding factors in randomized controlled traits, and our findings strongly support the hypothesis that smoking behavior is causally related to CAD risk. We found significant genetic overlap and correlations between CAD and smoking at the SNP level. Taken together, the data from this study enhance the understanding of the genetic etiology of the relations between CAD and smoking and might help to dissect smoking behaviors and develop preventive strategies to reduce the burden of cardiovascular disease in public health.

Funding

The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This study was supported by the Joint Institute of Tobacco and Health Open Project Fund, the China Precision Medicine Initiative (2016YFC0906300), Research Center for Air Pollution and Health of Zhejiang University, and the State Key Laboratory for Diagnosis and Treatment of Infectious Diseases of the First Affiliated Hospital of Zhejiang University.

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.

Statements

Data availability statement

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

Author contributions

ZZ: Writing – original draft, Writing – review & editing, Data curation, Formal analysis. QL: Data curation, Formal analysis, Writing – original draft, Methodology. MeL: Data curation, Formal analysis, Methodology, Writing – review & editing. YY: Data curation, Formal analysis, Methodology, Software, Writing – original draft. FQ: Data curation, Formal analysis, Software, Writing – original draft. YX: Data curation, Formal analysis, Methodology, Software, Writing – original draft. SL: Data curation, Formal analysis, Methodology, Writing – review & editing. ZY: Data curation, Project administration, Writing – review & editing. YG: Project administration, Writing – review & editing, Investigation, Resources, Supervision. MiL: Funding acquisition, Investigation, Resources, Writing – review & editing, Conceptualization, Supervision, Writing – original draft. JY: Funding acquisition, Investigation, Project administration, Resources, Writing – review & editing.

Acknowledgments

This research was conducted using GWAS summary data from PGC, GSCAN, and UK Biobank resources. We are deeply grateful to the investigators who comprised the PGC, GSCAN, and UK Biobank research teams and to all the subjects recruited into these projects.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision.

Supplementary material

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

References

  • 1.

    CDC. Current cigarette smoking among adults–United States. MMWR Morb Mortal Wkly Rep. (2012) 61:88994.

  • 2.

    World Health Organization. WHO report on the global tobacco epidemic, 2011: warning about the dangers of tobacco. Geneva, Switzerland: World Health Organization (2011).

  • 3.

    National Center for Chronic Disease Prevention and Health Promotion (US) Office on Smoking and Health In: The health consequences of Smoking-50 years of progress: a report of the surgeon general. Atlanta, GA: Centers for Disease Control and Prevention (US) (2014).

  • 4.

    MonsUMuezzinlerAGellertCSchottkerBAbnetCCBobakMet al. Impact of smoking and smoking cessation on cardiovascular events and mortality among older adults: meta-analysis of individual participant data from prospective cohort studies of the CHANCES consortium. BMJ. (2015) 350:h1551. doi: 10.1136/bmj.h1551

  • 5.

    DollRHillAB. Lung cancer and other causes of death in relation to smoking; a second report on the mortality of British doctors. Br Med J. (1956) 2:107181. doi: 10.1136/bmj.2.5001.1071

  • 6.

    JonasMAOatesJAOckeneJKHennekensCH. Statement on smoking and cardiovascular disease for health care professionals. Am Heart Assoc Circul. (1992) 86:16649. doi: 10.1161/01.cir.86.5.1664

  • 7.

    PriceJFMowbrayPILeeAJRumleyALoweGDFowkesFG. Relationship between smoking and cardiovascular risk factors in the development of peripheral arterial disease and coronary artery disease: Edinburgh artery study. Eur Heart J. (1999) 20:34453. doi: 10.1053/euhj.1998.1194

  • 8.

    WillettWCGreenAStampferMJSpeizerFEColditzGARosnerBet al. Relative and absolute excess risks of coronary heart disease among women who smoke cigarettes. N Engl J Med. (1987) 317:13039. doi: 10.1056/NEJM198711193172102

  • 9.

    FranklinWAGazdarAFHaneyJWistubaIILa RosaFGKennedyTet al. Widely dispersed p53 mutation in respiratory epithelium. A novel mechanism for field carcinogenesis. J Clin Invest. (1997) 100:21337. doi: 10.1172/JCI119748

  • 10.

    WistubaIILamS., BehrensC., Virmani, A.K., Fong, K.M., LeRicheJ., et al. (1997). Molecular damage in the bronchial epithelium of current and former smokers. J Natl Cancer Inst89, 13661373. doi: 10.1093/jnci/89.18.1366

  • 11.

    Sabater-LlealMMalarstigAFolkersenLSoler ArtigasMBaldassarreDKavousiMet al. Common genetic determinants of lung function, subclinical atherosclerosis and risk of coronary artery disease. PLoS One. (2014) 9:e104082. doi: 10.1371/journal.pone.0104082

  • 12.

    BurnsDM. Cigarettes and cigarette smoking. Clin Chest Med. (1991) 12:63142. doi: 10.1016/S0272-5231(21)00813-3

  • 13.

    WangXLGrecoMSimASDuarteNWangJWilckenDE. Effect of CYP1A1 MspI polymorphism on cigarette smoking related coronary artery disease and diabetes. Atherosclerosis. (2002) 162:3917. doi: 10.1016/s0021-9150(01)00723-7

  • 14.

    WangXLSimASWangMXMurrellGATrudingerBWangJ. Genotype dependent and cigarette specific effects on endothelial nitric oxide synthase gene expression and enzyme activity. FEBS Lett. (2000) 471:4550. doi: 10.1016/s0014-5793(00)01356-9

  • 15.

    SmithGDEbrahimS. 'Mendelian randomization': can genetic epidemiology contribute to understanding environmental determinants of disease?Int J Epidemiol. (2003) 32:122. doi: 10.1093/ije/dyg070

  • 16.

    ChenLYaoYJinCWuSLiuQLiJet al. Integrative genomic analysis identified common regulatory networks underlying the correlation between coronary artery disease and plasma lipid levels. BMC Cardiovasc Disord. (2019) 19:310. doi: 10.1186/s12872-019-01271-9

  • 17.

    LiuMJiangYWedowRLiYBrazelDMChenFet al. Association studies of up to 1.2 million individuals yield new insights into the genetic etiology of tobacco and alcohol use. Nat Genet. (2019) 51:23744. doi: 10.1038/s41588-018-0307-5

  • 18.

    WillerCJSchmidtEMSenguptaSPelosoGMGustafssonSKanoniSet al. Discovery and refinement of loci associated with lipid levels. Nat Genet. (2013) 45:127483. doi: 10.1038/ng.2797

  • 19.

    Bulik-SullivanBFinucaneHKAnttilaVGusevADayFRLohPRet al. An atlas of genetic correlations across human diseases and traits. Nat Genet. (2015) 47:123641. doi: 10.1038/ng.3406

  • 20.

    HemaniGZhengJElsworthBWadeKHHaberlandVBairdDet al. The MR-base platform supports systematic causal inference across the human phenome. elife. (2018) 7. doi: 10.7554/eLife.34408

  • 21.

    PurcellSNealeBTodd-BrownKThomasLFerreiraMABenderDet al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. (2007) 81:55975. doi: 10.1086/519795

  • 22.

    BurgessSButterworthAThompsonSG. Mendelian randomization analysis with multiple genetic variants using summarized data. Genet Epidemiol. (2013) 37:65865. doi: 10.1002/gepi.21758

  • 23.

    ZhangBKirovSSnoddyJ. WebGestalt: an integrated system for exploring gene sets in various biological contexts. Nucleic Acids Res. (2005) 33:W7418. doi: 10.1093/nar/gki475

  • 24.

    ChenSYangFXuTWangYZhangKFuGet al. Smoking and coronary artery disease risk in patients with diabetes: a mendelian randomization study. Front Immunol. (2023) 14:891947. doi: 10.3389/fimmu.2023.891947

  • 25.

    JiaYWangRGuoDSunLShiMZhangKet al. Contribution of metabolic risk factors and lifestyle behaviors to cardiovascular disease: a mendelian randomization study. Nutr Metab Cardiovasc Dis. (2022) 32:197281. doi: 10.1016/j.numecd.2022.04.019

  • 26.

    LevinMGKlarinDAssimesTLFreibergMSIngelssonELynchJet al. Genetics of smoking and risk of atherosclerotic cardiovascular diseases: a mendelian randomization study. JAMA Netw Open. (2021) 4:e2034461. doi: 10.1001/jamanetworkopen.2020.34461

  • 27.

    MoonJYKwonHMKwonSWYoonSJKimJSLeeSJet al. Lipoprotein(a) and LDL particle size are related to the severity of coronary artery disease. Cardiology. (2007) 108:2829. doi: 10.1159/000099097

  • 28.

    van den BergVJVroegindeweyMMKardysIBoersmaEHaskardDHartleyAet al. Anti-oxidized LDL antibodies and coronary artery disease: a systematic review. Antioxidants. (2019) 8. doi: 10.3390/antiox8100484

  • 29.

    ZhaoXZhangHWXuRXGuoYLZhuCGWuNQet al. Oxidized-LDL is a useful marker for predicting the very early coronary artery disease and cardiovascular outcomes. Per Med. (2018) 15:5219. doi: 10.2217/pme-2018-0046

  • 30.

    HeineckeJW. The HDL proteome: a marker--and perhaps mediator--of coronary artery disease. J Lipid Res. (2009) 50:S16771. doi: 10.1194/jlr.R800097-JLR200

  • 31.

    NatarajanPCollierTSJinZLyassALiYIbrahimNEet al. Association of an HDL Apolipoproteomic score with coronary atherosclerosis and cardiovascular death. J Am Coll Cardiol. (2019) 73:213545. doi: 10.1016/j.jacc.2019.01.073

  • 32.

    Weissglas-VolkovDPajukantaP. Genetic causes of high and low serum HDL-cholesterol. J Lipid Res. (2010) 51:203257. doi: 10.1194/jlr.R004739

  • 33.

    GrallertHDupuisJBisJCDehghanABarbalicMBaumertJet al. Eight genetic loci associated with variation in lipoprotein-associated phospholipase A2 mass and activity and coronary heart disease: meta-analysis of genome-wide association studies from five community-based studies. Eur Heart J. (2012) 33:23851. doi: 10.1093/eurheartj/ehr372

  • 34.

    HorwitzTLamKChenYXiaYLiuC. A decade in psychiatric GWAS research. Mol Psychiatry. (2019) 24:37889. doi: 10.1038/s41380-018-0055-z

  • 35.

    CoonKDDunckleyTLStephanDA. A generic research paradigm for identification and validation of early molecular diagnostics and new therapeutics in common disorders. Mol Diagn Ther. (2007) 11:114. doi: 10.1007/BF03256218

  • 36.

    HeroldCHooliBVMullinKLiuTRoehrJTMattheisenMet al. Family-based association analyses of imputed genotypes reveal genome-wide significant association of Alzheimer's disease with OSBPL6, PTPRG, and PDCL3. Mol Psychiatry. (2016) 21:160812. doi: 10.1038/mp.2015.218

  • 37.

    KambohMIBarmadaMMDemirciFYMinsterRLCarrasquilloMMPankratzVSet al. Genome-wide association analysis of age-at-onset in Alzheimer's disease. Mol Psychiatry. (2012) 17:13406. doi: 10.1038/mp.2011.135

  • 38.

    KambohMIDemirciFYWangXMinsterRLCarrasquilloMMPankratzVSet al. Genome-wide association study of Alzheimer's disease. Transl Psychiatry. (2012) 2:e117. doi: 10.1038/tp.2012.45

  • 39.

    LiHWettenSLiLSt JeanPLUpmanyuRSurhLet al. Candidate single-nucleotide polymorphisms from a genomewide association study of Alzheimer disease. Arch Neurol. (2008) 65:4553. doi: 10.1001/archneurol.2007.3

  • 40.

    WebsterJAMyersAJPearsonJVCraigDWHu-LinceDCoonKDet al. Sorl1 as an Alzheimer's disease predisposition gene?Neurodegener Dis. (2008) 5:604. doi: 10.1159/000110789

  • 41.

    De JagerPLShulmanJMChibnikLBKeenanBTRajTWilsonRSet al. A genome-wide scan for common variants affecting the rate of age-related cognitive decline. Neurobiol Aging. (2012) 33:1017 e10115. doi: 10.1016/j.neurobiolaging.2011.09.033

Summary

Keywords

smoking, CAD, coronary artery disease, genome-wide association study, mendelian randomization analysis

Citation

Zhu Z, Liu Q, Li M, Yao Y, Qi F, Xu Y, Lu S, Yang Z, Guan Y, Li MD and Yao J (2023) Determination of genetic correlation between tobacco smoking and coronary artery disease. Front. Psychiatry 14:1279962. doi: 10.3389/fpsyt.2023.1279962

Received

19 August 2023

Accepted

13 September 2023

Published

26 September 2023

Volume

14 - 2023

Edited by

Huiping Zhang, Boston University, United States

Reviewed by

Yunlong Ma, Wenzhou Medical University, China; Xiangning Chen, University of Texas Health Science Center at Houston, United States

Updates

Copyright

*Correspondence: Ming D. Li, ; Jianhua Yao,

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics