Multi-ancestry genome-wide study identifies effector genes and druggable pathways for coronary artery calcification

The role of coffee and potential mediators in subclinical atherosclerosis: insights from Mendelian randomization study

Background and aims

Coffee contains many bioactive compounds, and its inconsistent association with subclinical atherosclerosis has been reported in observational studies. In this Mendelian randomization study, we investigated whether genetically predicted coffee consumption is associated with subclinical atherosclerosis, as well as the role of potential mediators.

Methods

We first conducted a two-sample Mendelian randomization analysis to examine the causal effect of coffee and its subtypes on subclinical atherosclerosis inferred from coronary artery calcification (CAC). Next, the significant results were validated using another independent dataset. Two-step Mendelian randomization analyses were utilized to evaluate the causal pathway from coffee to subclinical atherosclerosis through potential mediators, including blood pressure, blood lipids, body mass index, and glycated hemoglobin. Mendelian randomization analyses were performed using the multiplicative random effects inverse-variance weighted method as the main approach, followed by a series of complementary methods and sensitivity analyses.

Results

Coffee, filtered coffee, and instant coffee were associated with the risk of CAC (β = 0.79, 95% CI: 0.12 to 1.47, p = 0.022; β = 0.66, 95% CI: 0.17 to 1.15, p = 0.008; β = 0.66, 95% CI: 0.20 to 1.13, p = 0.005; respectively). While no significant causal relationship was found between decaffeinated coffee and CAC (β = -1.32, 95% CI: -2.67 to 0.04, p = 0.056). The association between coffee and CAC was validated in the replication analysis (β = 0.27, 95% CI: 0.07 to 0.48, p = 0.009). Body mass index mediated 39.98% of the effect of coffee on CAC (95% CI: 9.78 to 70.19%, p = 0.009), and 5.79% of the effect of instant coffee on CAC (95% CI: 0.54 to 11.04%, p = 0.030).

Conclusion

Our study suggests that coffee other than decaffeinated coffee increases the risk of subclinical atherosclerosis inferred from CAC. Body mass index mediated 39.98 and 5.79% of the causal effects of coffee and instant coffee on CAC, respectively. Coffee should be consumed with caution, especially in individuals with established cardiovascular risk factors, and decaffeinated coffee appears to be a safer choice.

A Large-Scale Genome-Wide Study of Gene-Sleep Duration Interactions for Blood Pressure in 811,405 Individuals from Diverse Populations

Although both short and long sleep duration are associated with elevated hypertension risk, our understanding of their interplay with biological pathways governing blood pressure remains limited. To address this, we carried out genome-wide cross-population gene-by-short-sleep and long-sleep duration interaction analyses for three blood pressure traits (systolic, diastolic, and pulse pressure) in 811,405 individuals from diverse population groups. We discover 22 novel gene-sleep duration interaction loci for blood pressure, mapped to 23 genes. Investigating these genes' functional implications shed light on neurological, thyroidal, bone metabolism, and hematopoietic pathways that necessitate future investigation for blood pressure management that caters to sleep health lifestyle. Non-overlap between short sleep (12) and long sleep (10) interactions underscores the plausible nature of distinct influences of both sleep duration extremes in cardiovascular health. Several of our loci are specific towards a particular population background or sex, emphasizing the importance of addressing heterogeneity entangled in gene-environment interactions, when considering precision medicine design approaches for blood pressure management.

Association of PHACTR1 with Coronary Artery Calcium Differs by Sex and Cigarette Smoking

Background: Coronary artery calcium (CAC) is a marker of subclinical atherosclerosis and is a complex heritable trait with both genetic and environmental risk factors, including sex and smoking. Methods: We performed genome-wide association (GWA) analyses for CAC among all participants and stratified by sex in the COPDGene study (n = 6144 participants of European ancestry and n = 2589 participants of African ancestry) with replication in the Diabetes Heart Study (DHS). We adjusted for age, sex, current smoking status, BMI, diabetes, self-reported high blood pressure, self-reported high cholesterol, and genetic ancestry (as summarized by principal components computed within each racial group). For the significant signals from the GWA analyses, we examined the single nucleotide polymorphism (SNP) by sex interactions, stratified by smoking status (current vs. former), and tested for a SNP by smoking status interaction on CAC. Results: We identified genome-wide significant associations for CAC in the chromosome 9p21 region [CDKN2B-AS1] among all COPDGene participants (p = 7.1 × 10-14) and among males (p = 1.0 × 10-9), but the signal was not genome-wide significant among females (p = 6.4 × 10-6). For the sex stratified GWA analyses among females, the chromosome 6p24 region [PHACTR1] had a genome-wide significant association (p = 4.4 × 10-8) with CAC, but this signal was not genome-wide significant among all COPDGene participants (p = 1.7 × 10-7) or males (p = 0.03). There was a significant interaction for the SNP rs9349379 in PHACTR1 with sex (p = 0.02), but the interaction was not significant for the SNP rs10757272 in CDKN2B-AS1 with sex (p = 0.21). In addition, PHACTR1 had a stronger association with CAC among current smokers (p = 6.2 × 10-7) than former smokers (p = 7.5 × 10-3) and the SNP by smoking status interaction was marginally significant (p = 0.03). CDKN2B-AS1 had a strong association with CAC among both former (p = 7.7 × 10-8) and current smokers (p = 1.7 × 10-7) and the SNP by smoking status interaction was not significant (p = 0.40). Conclusions: Among current and former smokers of European ancestry in the COPDGene study, we identified a genome-wide significant association in the chromosome 6p24 region [PHACTR1] with CAC among females, but not among males. This region had a significant SNP by sex and SNP by smoking interaction on CAC.

Associations of genetically predicted vitamin D status and deficiency with the risk of carotid artery plaque: a Mendelian randomization study

Low concentrations of circulating 25-hydroxy-vitamin D are observationally associated with an increased risk of subclinical atherosclerosis and cardiovascular disease. However, randomized controlled trials have not reported the beneficial effects of vitamin D supplementation on atherosclerotic cardiovascular disease (ASCVD) outcomes. Whether genetically predicted vitamin D status confers protection against the development of carotid artery plaque, a powerful predictor of subclinical atherosclerosis, remains unknown. We conducted a two-sample Mendelian randomization (MR) study to explore the association of genetically predicted vitamin D status and deficiency with the risk of developing carotid artery plaque. We leveraged three genome-wide association studies (GWAS) of vitamin D status and one GWAS of vitamin D deficiency. We used the inverse-variance weighted (IVW) approach as our main method, and MR-Egger, weighted-median, and radialMR as MR sensitivity analyses. We also conducted sensitivity analyses using biologically plausible genetic instruments located within genes encoding for vitamin D metabolism (GC, CYP2R1, DHCR7, CYP24A1). We did not find significant associations between genetically predicted vitamin D status (Odds ratio (OR) = 0.99, P = 0.91) and deficiency (OR = 1.00, P = 0.97) with the risk of carotid artery plaque. We additionally explored the potential causal effect of vitamin D status on coronary artery calcification (CAC) and carotid intima-media thickness (cIMT), two additional markers of subclinical atherosclerosis, and we did not find any significant association (βCAC = - 0.14, P = 0.23; βcIMT = 0.005, P = 0.19). These findings did not support the causal effects of vitamin D status and deficiency on the risk of developing subclinical atherosclerosis.

Exploring the mediating role of calcium homeostasis in the association between diabetes mellitus, glycemic traits, and vascular and valvular calcifications: a comprehensive Mendelian randomization analysis

BACKGROUND

The interplay between diabetes mellitus (DM), glycemic traits, and vascular and valvular calcifications is intricate and multifactorial. Exploring potential mediators may illuminate underlying pathways and identify novel therapeutic targets.

METHODS

We utilized univariable and multivariable Mendelian randomization (MR) analyses to investigate associations and mediation effects. Additionally, the multivariable MR analyses incorporated cardiometabolic risk factors, allowing us to account for potential confounders.

RESULTS

Type 2 diabetes mellitus (T2DM) and glycated hemoglobin (HbA1c) were positively associated with both coronary artery calcification (CAC) and calcific aortic valvular stenosis (CAVS). However, fasting glucose (FG) was only linked to CAVS and showed no association with CAC. Additionally, CAVS demonstrated a causal effect on FG. Calcium levels partially mediated the impact of T2DM on both types of calcifications. Specifically, serum calcium was positively associated with both CAC and CAVS. The mediation effects of calcium levels on the impact of T2DM on CAC and CAVS were 6.063% and 3.939%, respectively. The associations between T2DM and HbA1c with calcifications were influenced by body mass index (BMI) and smoking status. However, these associations were generally reduced after adjusting for hypertension.

CONCLUSION

Our findings suggest a genetically supported causal relationship between DM, glycemic traits, and vascular and valvular calcifications, with serum calcium playing a critical mediating role.

Unraveling interindividual variation of trimethylamine N-oxide and its precursors at the population level

Trimethylamine N-oxide (TMAO) is a circulating microbiome-derived metabolite implicated in the development of atherosclerosis and cardiovascular disease (CVD). We investigated whether plasma levels of TMAO, its precursors (betaine, carnitine, deoxycarnitine, choline), and TMAO-to-precursor ratios are associated with clinical outcomes, including CVD and mortality. This was followed by an in-depth analysis of their genetic, gut microbial, and dietary determinants. The analyses were conducted in five Dutch prospective cohort studies including 7834 individuals. To further investigate association results, Mendelian Randomization (MR) was also explored. We found only plasma choline levels (hazard ratio [HR] 1.17, [95% CI 1.07; 1.28]) and not TMAO to be associated with CVD risk. Our association analyses uncovered 10 genome-wide significant loci, including novel genomic regions for betaine (6p21.1, 6q25.3), choline (2q34, 5q31.1), and deoxycarnitine (10q21.2, 11p14.2) comprising several metabolic gene associations, for example, CPS1 or PEMT. Furthermore, our analyses uncovered 68 gut microbiota associations, mainly related to TMAO-to-precursors ratios and the Ruminococcaceae family, and 16 associations of food groups and metabolites including fish-TMAO, meat-carnitine, and plant-based food-betaine associations. No significant association was identified by the MR approach. Our analyses provide novel insights into the TMAO pathway, its determinants, and pathophysiological impact on the general population.

Association between BMD and coronary artery calcification: an observational and Mendelian randomization study

Observational studies have reported inconsistent associations between bone mineral density (BMD) and coronary artery calcification (CAC). We examined the observational association of BMD with CAC in 2 large population-based studies and evaluated the evidence for a potential causal relation between BMD and CAC using polygenic risk scores (PRS), 1- and 2-sample Mendelian randomization (MR) approaches. Our study populations comprised 1414 individuals (mean age 69.9 yr, 52.0% women) from the Rotterdam Study and 2233 individuals (mean age 56.5 yr, 50.9% women) from the Framingham Heart Study with complete information on CAC and BMD measurements at the total body (TB-), lumbar spine (LS-), and femoral neck (FN-). We used linear regression models to evaluate the observational association between BMD and CAC. Subsequently, we compared the mean CAC across PRSBMD quintile groups at different skeletal sites. In addition, we used the 2-stage least squares regression and the inverse variance weighted (IVW) model as primary methods for 1- and 2-sample MR to test evidence for a potentially causal association. We did not observe robust associations between measured BMD levels and CAC. These results were consistent with a uniform random distribution of mean CAC across PRSBMD quintile groups (P-value > .05). Moreover, neither 1- nor 2-sample MR supported the possible causal association between BMD and CAC. Our results do not support the contention that lower BMD is (causally) associated with an increased CAC risk. These findings suggest that previously reported epidemiological associations of BMD with CAC are likely explained by unmeasured confounders or shared etiology, rather than by causal pathways underlying both osteoporosis and vascular calcification processes.

A Large-Scale Genome-Wide Study of Gene-Sleep Duration Interactions for Blood Pressure in 811,405 Individuals from Diverse Populations

Although both short and long sleep duration are associated with elevated hypertension risk, our understanding of their interplay with biological pathways governing blood pressure remains limited. To address this, we carried out genome-wide cross-population gene-by-short-sleep and long-sleep duration interaction analyses for three blood pressure traits (systolic, diastolic, and pulse pressure) in 811,405 individuals from diverse population groups. We discover 22 novel gene-sleep duration interaction loci for blood pressure, mapped to genes involved in neurological, thyroidal, bone metabolism, and hematopoietic pathways. Non-overlap between short sleep (12) and long sleep (10) interactions underscores the plausibility of distinct influences of both sleep duration extremes in cardiovascular health. With several of our loci reflecting specificity towards population background or sex, our discovery sheds light on the importance of embracing granularity when addressing heterogeneity entangled in gene-environment interactions, and in therapeutic design approaches for blood pressure management.

Multi-ancestry genetic analysis of gene regulation in coronary arteries prioritizes disease risk loci

Genome-wide association studies (GWASs) have identified hundreds of risk loci for coronary artery disease (CAD). However, non-European populations are underrepresented in GWASs, and the causal gene-regulatory mechanisms of these risk loci during atherosclerosis remain unclear. We incorporated local ancestry and haplotypes to identify quantitative trait loci for expression (eQTLs) and splicing (sQTLs) in coronary arteries from 138 ancestrally diverse Americans. Of 2,132 eQTL-associated genes (eGenes), 47% were previously unreported in coronary artery; 19% exhibited cell-type-specific expression. Colocalization revealed subgroups of eGenes unique to CAD and blood pressure GWAS. Fine-mapping highlighted additional eGenes, including TBX20 and IL5. We also identified sQTLs for 1,690 genes, among which TOR1AIP1 and ULK3 sQTLs demonstrated the importance of evaluating splicing to accurately identify disease-relevant isoform expression. Our work provides a patient-derived coronary artery eQTL resource and exemplifies the need for diverse study populations and multifaceted approaches to characterize gene regulation in disease processes.

Unravelling novel and pleiotropic genes for cannon bone circumference and bone mineral density in Yorkshire pigs

Leg weakness is a prevalent health condition in pig farms. The augmentation of cannon bone circumference and bone mineral density can effectively improve limb strength in pigs and alleviate leg weakness. This study measured forelimb cannon bone circumference (fCBC) and rear limb cannon bone circumference (rCBC) using an inelastic tapeline and rear limb metatarsal area bone mineral density (raBMD) using a dual-energy X-ray absorptiometry bone density scanner. The samples of Yorkshire castrated boars were genotyped using a 50K single-nucleotide polymorphism (SNP) array. The SNP-chip data were imputed to the level of whole-genome sequencing data (iWGS). This study used iWGS data to perform genome-wide association studies and identified novel significant SNPs associated with fCBC on SSC6, SSC12, and SSC13, rCBC on SSC12 and SSC14, and raBMD on SSC7. Based on the high phenotypic and genetic correlations between CBC and raBMD, multi-trait meta-analysis was performed to identify pleiotropic SNPs. A significant potential pleiotropic quantitative trait locus (QTL) regulating both CBC and raBMD was identified on SSC15. Bayes fine mapping was used to establish the confidence intervals for these novel QTLs with the most refined confidence interval narrowed down to 56 kb (15.11 to 15.17 Mb on SSC12 for fCBC). Furthermore, the confidence interval for the potential pleiotropic QTL on SSC15 in the meta-analysis was narrowed down to 7.45 kb (137.55 to137.56 Mb on SSC15). Based on the biological functions of genes, the following genes were identified as novel regulatory candidates for different phenotypes: DDX42, MYSM1, FTSJ3, and MECOM for fCBC; SMURF2, and STC1 for rCBC; RGMA for raBMD. Additionally, RAMP1, which was determined to be located 23.68 kb upstream of the confidence interval of the QTL on SSC15 in the meta-analysis, was identified as a potential pleiotropic candidate gene regulating both CBC and raBMD. These findings offered valuable insights for identifying pathogenic genes and elucidating the genetic mechanisms underlying CBC and BMD.

Genetic Variants at the 9p21.3 Locus Are Associated with Risk for Non-Compressible Artery Disease: Results from the ARTPER Study

Peripheral artery disease (PAD) and non-compressible artery disease (NCAD) constitute predictors of subclinical atherosclerosis easily assessed through the ankle brachial index (ABI). Although both diseases show substantial genetic influences, few genetic association studies have focused on the ABI and PAD, and none have focused on NCAD. To overcome these limitations, we assessed the role of several candidate genes on the ABI, both in its continuous distribution and in the clinical manifestations associated to its extreme values: PAD and NCAD. We examined 13 candidate genomic regions in 1606 participants from the ARTPER study, a prospective population-based cohort, with the ABI assessed through ultrasonography. Association analyses were conducted independently for individuals with PAD (ABI < 0.9) or with NCAD (ABI > 1.4) vs. healthy participants. After including potential covariates and correction for multiple testing, minor alleles in the genetic markers rs10757278 and rs1333049, both in the 9p21.3 region, were significantly associated with a decreased risk of NCAD. Associations with the ABI showed limited support to these results. No significant associations were detected for PAD. The locus 9p21.3 constitutes the first genetic locus associated with NCAD, an assessment of subclinical atherosclerosis feasible for implementation in primary healthcare settings that has been systematically neglected from genetic studies.

Beyond the Basics: Unraveling the Complexity of Coronary Artery Calcification

Coronary artery calcification (CAC) is mainly associated with coronary atherosclerosis, which is an indicator of coronary artery disease (CAD). CAC refers to the accumulation of calcium phosphate deposits, classified as micro- or macrocalcifications, that lead to the hardening and narrowing of the coronary arteries. CAC is a strong predictor of future cardiovascular events, such as myocardial infarction and sudden death. Our narrative review focuses on the pathophysiology of CAC, exploring its link to plaque vulnerability, genetic factors, and how race and sex can affect the condition. We also examined the connection between the gut microbiome and CAC, and the impact of genetic variants on the cellular processes involved in vascular calcification and atherogenesis. We aimed to thoroughly analyze the existing literature to improve our understanding of CAC and its potential clinical and therapeutic implications.

Type 2 Diabetes Modifies the Association of CAD Genomic Risk Variants With Subclinical Atherosclerosis

BACKGROUND

Individuals with type 2 diabetes (T2D) have an increased risk of coronary artery disease (CAD), but questions remain about the underlying pathology. Identifying which CAD loci are modified by T2D in the development of subclinical atherosclerosis (coronary artery calcification [CAC], carotid intima-media thickness, or carotid plaque) may improve our understanding of the mechanisms leading to the increased CAD in T2D.

METHODS

We compared the common and rare variant associations of known CAD loci from the literature on CAC, carotid intima-media thickness, and carotid plaque in up to 29 670 participants, including up to 24 157 normoglycemic controls and 5513 T2D cases leveraging whole-genome sequencing data from the Trans-Omics for Precision Medicine program. We included first-order T2D interaction terms in each model to determine whether CAD loci were modified by T2D. The genetic main and interaction effects were assessed using a joint test to determine whether a CAD variant, or gene-based rare variant set, was associated with the respective subclinical atherosclerosis measures and then further determined whether these loci had a significant interaction test.

RESULTS

Using a Bonferroni-corrected significance threshold of P<1.6×10-4, we identified 3 genes (ATP1B1, ARVCF, and LIPG) associated with CAC and 2 genes (ABCG8 and EIF2B2) associated with carotid intima-media thickness and carotid plaque, respectively, through gene-based rare variant set analysis. Both ATP1B1 and ARVCF also had significantly different associations for CAC in T2D cases versus controls. No significant interaction tests were identified through the candidate single-variant analysis.

CONCLUSIONS

These results highlight T2D as an important modifier of rare variant associations in CAD loci with CAC.

Integrative single-cell meta-analysis reveals disease-relevant vascular cell states and markers in human atherosclerosis

Coronary artery disease (CAD) is characterized by atherosclerotic plaque formation in the arterial wall. CAD progression involves complex interactions and phenotypic plasticity among vascular and immune cell lineages. Single-cell RNA-seq (scRNA-seq) studies have highlighted lineage-specific transcriptomic signatures, but human cell phenotypes remain controversial. Here, we perform an integrated meta-analysis of 22 scRNA-seq libraries to generate a comprehensive map of human atherosclerosis with 118,578 cells. Besides characterizing granular cell-type diversity and communication, we leverage this atlas to provide insights into smooth muscle cell (SMC) modulation. We integrate genome-wide association study data and uncover a critical role for modulated SMC phenotypes in CAD, myocardial infarction, and coronary calcification. Finally, we identify fibromyocyte/fibrochondrogenic SMC markers (LTBP1 and CRTAC1) as proxies of atherosclerosis progression and validate these through omics and spatial imaging analyses. Altogether, we create a unified atlas of human atherosclerosis informing cell state-specific mechanistic and translational studies of cardiovascular diseases.