Novel Blood Pressure Locus and Gene Discovery Using Genome-Wide Association Study and Expression Data Sets From Blood and the Kidney

The Genetic Variants Influencing Hypertension Prevalence Based on the Risk of Insulin Resistance as Assessed Using the Metabolic Score for Insulin Resistance (METS-IR)

Insulin resistance is a major indicator of cardiovascular diseases, including hypertension. The Metabolic Score for Insulin Resistance (METS-IR) offers a simplified and cost-effective way to evaluate insulin resistance. This study aimed to identify genetic variants associated with the prevalence of hypertension stratified by METS-IR score levels. Data from the Korean Genome and Epidemiology Study (KoGES) were analyzed. The METS-IR was calculated using the following formula: ln [(2 × fasting blood glucose (FBG) + triglycerides (TG)) × body mass index (BMI)]/ ln [high-density lipoprotein cholesterol (HDL-C)]. The participants were divided into tertiles 1 (T1) and 3 (T3) based on their METS-IR scores. Genome-wide association studies (GWAS) were performed for hypertensive cases and non-hypertensive controls within these tertile groups using logistic regression adjusted for age, sex, and lifestyle factors. Among the METS-IR tertile groups, 3517 of the 19,774 participants (17.8%) at T1 had hypertension, whereas 8653 of the 20,374 participants (42.5%) at T3 had hypertension. A total of 113 single-nucleotide polymorphisms (SNPs) reached the GWAS significance threshold (p < 5 × 10-8) in at least one tertile group, mapping to six distinct genetic loci. Notably, four loci, rs11899121 (chr2p24), rs7556898 (chr2q24.3), rs17249754 (ATP2B1), and rs1980854 (chr20p12.2), were significantly associated with hypertension in the high-METS-score group (T3). rs10857147 (FGF5) was significant in both the T1 and T3 groups, whereas rs671 (ALDH2) was significant only in the T1 group. The GWASs identified six genetic loci significantly associated with hypertension, with distinct patterns across METS-IR tertiles, highlighting the role of metabolic context in genetic susceptibility. These findings underscore critical genetic factors influencing hypertension prevalence and provide insights into the metabolic-genetic interplay underlying this condition.

Is the Relationship Between Cardiovascular Disease and Alzheimer's Disease Genetic? A Scoping Review

BACKGROUND/OBJECTIVES

Cardiovascular disease (CVD) and Alzheimer's disease (AD) are two diseases highly prevalent in the aging population and often co-occur. The exact relationship between the two diseases is uncertain, though epidemiological studies have demonstrated that CVDs appear to increase the risk of AD and vice versa. This scoping review aims to examine the current identified overlapping genetics between CVDs and AD at the individual gene level and at the shared pathway level.

METHODS

Following PRISMA-ScR guidelines for a scoping review, we searched the PubMed and Scopus databases from 1990 to October 2024 for articles that involved (1) CVDs, (2) AD, and (3) used statistical methods to parse genetic relationships.

RESULTS

Our search yielded 2918 articles, of which 274 articles passed screening and were organized into two main sections: (1) evidence of shared genetic risk; and (2) shared mechanisms. The genes APOE, PSEN1, and PSEN2 reportedly have wide effects across the AD and CVD spectrum, affecting both cardiac and brain tissues. Mechanistically, changes in three main pathways (lipid metabolism, blood pressure regulation, and the breakdown of the blood-brain barrier (BBB)) contribute to subclinical and etiological changes that promote both AD and CVD progression. However, genetic studies continue to be limited by the availability of longitudinal data and lack of cohorts that are representative of diverse populations.

CONCLUSIONS

Highly penetrant familial genes simultaneously increase the risk of CVDs and AD. However, in most cases, sets of dysregulated genes within larger-scale mechanisms, like changes in lipid metabolism, blood pressure regulation, and BBB breakdown, increase the risk of both AD and CVDs and contribute to disease progression.

The Promoter Polymorphism rs3918226 of the Endothelial Nitric Oxide Synthase Gene as a Novel Susceptibility Marker for Peripheral Artery Disease

BACKGROUND

This pilot study aimed to investigate the association between the single nucleotide polymorphism (SNP) rs3918226 in the promoter of the nitric oxide synthase (NOS3) gene and the risk of peripheral artery disease (PAD).

METHODS

DNA samples from 1,263 unrelated subjects of Slavic origin, including 620 patients with PAD and 643 controls, were genotyped for the SNP rs3918226 using the MassArray-4 system.

RESULTS

The rs3918226 polymorphism was found to be strongly associated with an increased risk of PAD regardless of coronary artery disease, hypertension, or cigarette smoking (odds ratio [OR] = 2.86; 95% confidence interval [CI] 1.89-4.32; Pperm < 0.0001). The SNP-PAD association was almost 3 times stronger in females (OR = 8.31; 95% CI 3.07-22.48) than in males (OR = 1.79; 95% CI 1.10-2.93). SNP rs3918226 was correlated with ankle-brachial index and total plasma cholesterol in patients with PAD (Рperm < 0.05). The NOS3 polymorphism was closely associated with SNPs rs7692387 and rs13139571 in guanylate cyclase soluble subunit alpha-3 (GUCY1A3) to determine the risk of PAD, suggesting that the rs3918226 polymorphism may disrupt signaling in the NO-soluble guanylyl cyclase pathway. Diplotypes with wild-type alleles, such as NOS3 rs3918226-C/C×GUCY1A1 rs7692387G/G and NOS3 rs3918226-C/C×GUCY1A1 rs13139571C/C, showed strong protection against disease risk (false discovery rate ≤ 0.001). Functional SNP annotation revealed that the allele rs3918226-T was associated with decreased expression of NOS3, most strongly in the tibial arteries than in the coronary artery or aorta.

CONCLUSIONS

The present study is the first to show that the rs3918226 polymorphism of NOS3 is a novel susceptibility marker for PAD. Further research in independent populations is necessary to reproduce the association between polymorphism rs3918226 and disease risk.

Mitochondrial related variants associated with cardiovascular traits

Introduction

Cardiovascular disease (CVD) is responsible for over 30% of mortality worldwide. CVD arises from the complex influence of molecular, clinical, social, and environmental factors. Despite the growing number of autosomal genetic variants contributing to CVD, the cause of most CVDs is still unclear. Mitochondria are crucial in the pathophysiology, development and progression of CVDs; the impact of mitochondrial DNA (mtDNA) variants and mitochondrial haplogroups in the context of CVD has recently been highlighted.

Aims

We investigated the role of genetic variants in both mtDNA and nuclear-encoded mitochondrial genes (NEMG) in CVD, including coronary artery disease (CAD), hypertension, and serum lipids in the UK Biobank, with sub-group analysis for diabetes.

Methods

We investigated 371,542 variants in 2,527 NEMG, along with 192 variants in 32 mitochondrial genes in 381,994 participants of the UK Biobank, stratifying by presence of diabetes.

Results

Mitochondrial variants showed associations with CVD, hypertension, and serum lipids. Mitochondrial haplogroup J was associated with CAD and serum lipids, whereas mitochondrial haplogroups T and U were associated with CVD. Among NEMG, variants within Nitric Oxide Synthase 3 (NOS3) showed associations with CVD, CAD, hypertension, as well as diastolic and systolic blood pressure. We also identified Translocase Of Outer Mitochondrial Membrane 40 (TOMM40) variants associated with CAD; Solute carrier family 22 member 2 (SLC22A2) variants associated with CAD and CVD; and HLA-DQA1 variants associated with hypertension. Variants within these three genes were also associated with serum lipids.

Conclusion

Our study demonstrates the relevance of mitochondrial related variants in the context of CVD. We have linked mitochondrial haplogroup U to CVD, confirmed association of mitochondrial haplogroups J and T with CVD and proposed new markers of hypertension and serum lipids in the context of diabetes. We have also evidenced connections between the etiological pathways underlying CVDs, blood pressure and serum lipids, placing NOS3, SLC22A2, TOMM40 and HLA-DQA1 genes as common nexuses.

Genetic and Epigenetic Mechanisms Regulating Blood Pressure and Kidney Dysfunction

The pioneering work of Dr Lewis K. Dahl established a relationship between kidney, salt, and high blood pressure (BP), which led to the major genetic-based experimental model of hypertension. BP, a heritable quantitative trait affected by numerous biological and environmental stimuli, is a major cause of morbidity and mortality worldwide and is considered to be a primary modifiable factor in renal, cardiovascular, and cerebrovascular diseases. Genome-wide association studies have identified monogenic and polygenic variants affecting BP in humans. Single nucleotide polymorphisms identified in genome-wide association studies have quantified the heritability of BP and the effect of genetics on hypertensive phenotype. Changes in the transcriptional program of genes may represent consequential determinants of BP, so understanding the mechanisms of the disease process has become a priority in the field. At the molecular level, the onset of hypertension is associated with reprogramming of gene expression influenced by epigenomics. This review highlights the specific genetic variants, mutations, and epigenetic factors associated with high BP and how these mechanisms affect the regulation of hypertension and kidney dysfunction.

Comparing the effects of CETP in East Asian and European ancestries: a Mendelian randomization study

CETP inhibitors are a class of lipid-lowering drugs in development for treatment of coronary heart disease (CHD). Genetic studies in East Asian ancestry have interpreted the lack of CETP signal with low-density lipoprotein cholesterol (LDL-C) and lack of drug target Mendelian randomization (MR) effect on CHD as evidence that CETP inhibitors might not be effective in East Asian participants. Capitalizing on recent increases in sample size of East Asian genetic studies, we conducted a drug target MR analysis, scaled to a standard deviation increase in high-density lipoprotein cholesterol. Despite finding evidence for possible neutral effects of lower CETP levels on LDL-C, systolic blood pressure and pulse pressure in East Asians (interaction p-values < 1.6 × 10-3), effects on cardiovascular outcomes were similarly protective in both ancestry groups. In conclusion, on-target inhibition of CETP is anticipated to decrease cardiovascular disease in individuals of both European and East Asian ancestries.

An endothelial regulatory module links blood pressure regulation with elite athletic performance

The control of transcription is crucial for homeostasis in mammals. A previous selective sweep analysis of horse racing performance revealed a 19.6 kb candidate regulatory region 50 kb downstream of the Endothelin3 (EDN3) gene. Here, the region was narrowed to a 5.5 kb span of 14 SNVs, with elite and sub-elite haplotypes analyzed for association to racing performance, blood pressure and plasma levels of EDN3 in Coldblooded trotters and Standardbreds. Comparative analysis of human HiCap data identified the span as an enhancer cluster active in endothelial cells, interacting with genes relevant to blood pressure regulation. Coldblooded trotters with the sub-elite haplotype had significantly higher blood pressure compared to horses with the elite performing haplotype during exercise. Alleles within the elite haplotype were part of the standing variation in pre-domestication horses, and have risen in frequency during the era of breed development and selection. These results advance our understanding of the molecular genetics of athletic performance and vascular traits in both horses and humans.

Unveiling blood pressure-associated genes in aortic cells through integrative analysis of GWAS and RNA modification-associated variants

Background

Genome-wide association studies (GWAS) have identified more than a thousand loci for blood pressure (BP). Functional genes in these loci are cell-type specific. The aim of this study was to elucidate potentially functional genes associated with BP in the aorta through the utilization of RNA modification-associated single-nucleotide polymorphisms (RNAm-SNPs).

Methods

Utilizing large-scale genetic data of 757,601 individuals from the UK Biobank and International Consortium of Blood Pressure consortium, we identified associations between RNAm-SNPs and BP. The association between RNAm-SNPs, gene expression, and BP were examined.

Results

A total of 355 RNAm-SNPs related to m6A, m1A, m5C, m7G, and A-to-I modification were associated with BP. The related genes were enriched in the pancreatic secretion pathway and renin secretion pathway. The BP GWAS signals were significantly enriched with m6A-SNPs, highlighting the potential functional relevance of m6A in physiological processes influencing BP. Notably, m6A-SNPs in CYP11B1, PDE3B, HDAC7, ACE, SLC4A7, PDE1A, FRK, MTHFR, NPPA, CACNA1D, and HDAC9 were identified. Differential methylation and differential expression of the BP genes in FTO-overexpression and METTL14-knockdown vascular smooth muscle cells were detected. RNAm-SNPs were associated with ascending and descending aorta diameter and the genes showed differential methylation between aortic dissection (AD) cases and controls. In scRNA-seq study, we identified ARID5A, HLA-DPB1, HLA-DRA, IRF1, LINC01091, MCL1, MLF1, MLXIPL, NAA16, NADK, RERG, SRM, and USP53 as differential expression genes for AD in aortic cells.

Conclusion

The present study identified RNAm-SNPs in BP loci and elucidated the associations between the RNAm-SNPs, gene expression, and BP. The identified BP-associated genes in aortic cells were associated with AD.

Genome-wide analysis in over 1 million individuals of European ancestry yields improved polygenic risk scores for blood pressure traits

Hypertension affects more than one billion people worldwide. Here we identify 113 novel loci, reporting a total of 2,103 independent genetic signals (P < 5 × 10-8) from the largest single-stage blood pressure (BP) genome-wide association study to date (n = 1,028,980 European individuals). These associations explain more than 60% of single nucleotide polymorphism-based BP heritability. Comparing top versus bottom deciles of polygenic risk scores (PRSs) reveals clinically meaningful differences in BP (16.9 mmHg systolic BP, 95% CI, 15.5-18.2 mmHg, P = 2.22 × 10-126) and more than a sevenfold higher odds of hypertension risk (odds ratio, 7.33; 95% CI, 5.54-9.70; P = 4.13 × 10-44) in an independent dataset. Adding PRS into hypertension-prediction models increased the area under the receiver operating characteristic curve (AUROC) from 0.791 (95% CI, 0.781-0.801) to 0.826 (95% CI, 0.817-0.836, ∆AUROC, 0.035, P = 1.98 × 10-34). We compare the 2,103 loci results in non-European ancestries and show significant PRS associations in a large African-American sample. Secondary analyses implicate 500 genes previously unreported for BP. Our study highlights the role of increasingly large genomic studies for precision health research.

Impact of RRP1B Variants on the Phenotype, Progression, and Metastasis of Cervical Cancer

Metastasis is a key determinant of cancer progression, influenced significantly by genetic mechanisms. RRP1B, primarily a nucleolar protein, emerges as a suppressor of metastasis, forming alliances with various cellular components and modulating gene expression. This study investigates the involvement of the ribosomal RNA processing 1 homolog B (RRP1B) gene in metastasis regulation in cervical cancer. Through a comprehensive analysis of 172 cervical cancer patients, we evaluated five RRP1B single nucleotide polymorphisms (SNPs) (rs2838342, rs7276633, rs2051407, rs9306160, and rs762400) for their associations with clinicopathological features and survival outcomes. Significant associations were observed between specific genetic variants and clinicopathological parameters. Notably, the A allele of rs2838342 was associated with reduced odds of advanced tumor size, worse prognosis, and, preliminarily, distant metastasis, while the T allele of rs7276633 correlated with a decreased risk of higher tumor size and worse prognosis. Additionally, the C allele of rs2051407 demonstrated protective effects against larger tumors, metastasis, and adverse prognosis. The rs9306160 C allele exhibited a protective effect against metastasis. The rs762400 G allele was significant for reduced tumor size and metastasis risk. Furthermore, the rs2838342 A allele, rs7276633 T allele, rs2051407 C allele, and rs762400 G allele were associated with improved overall survival, demonstrating their potential significance in predicting prognoses in cervical cancer. Linkage disequilibrium and haplotypes analysis enabled us to evaluate the collective effect of the analyzed SNPs, which was in line with the results of allelic models. Our findings underscore the clinical relevance of RRP1B SNPs as prognostic markers in cervical cancer, shedding light on the intricate interplay between genetic factors and disease-progression dynamics. This research provides critical insights for future investigations and underscores the importance of incorporating RRP1B SNP detection into prognostic-assessment tools for accurate prediction of disease outcomes in cervical cancer.

Genetic imputation of kidney transcriptome, proteome and multi-omics illuminates new blood pressure and hypertension targets

Genetic mechanisms of blood pressure (BP) regulation remain poorly defined. Using kidney-specific epigenomic annotations and 3D genome information we generated and validated gene expression prediction models for the purpose of transcriptome-wide association studies in 700 human kidneys. We identified 889 kidney genes associated with BP of which 399 were prioritised as contributors to BP regulation. Imputation of kidney proteome and microRNAome uncovered 97 renal proteins and 11 miRNAs associated with BP. Integration with plasma proteomics and metabolomics illuminated circulating levels of myo-inositol, 4-guanidinobutanoate and angiotensinogen as downstream effectors of several kidney BP genes (SLC5A11, AGMAT, AGT, respectively). We showed that genetically determined reduction in renal expression may mimic the effects of rare loss-of-function variants on kidney mRNA/protein and lead to an increase in BP (e.g., ENPEP). We demonstrated a strong correlation (r = 0.81) in expression of protein-coding genes between cells harvested from urine and the kidney highlighting a diagnostic potential of urinary cell transcriptomics. We uncovered adenylyl cyclase activators as a repurposing opportunity for hypertension and illustrated examples of BP-elevating effects of anticancer drugs (e.g. tubulin polymerisation inhibitors). Collectively, our studies provide new biological insights into genetic regulation of BP with potential to drive clinical translation in hypertension.

Genetic architecture of cardiac dynamic flow volumes

Cardiac blood flow is a critical determinant of human health. However, the definition of its genetic architecture is limited by the technical challenge of capturing dynamic flow volumes from cardiac imaging at scale. We present DeepFlow, a deep-learning system to extract cardiac flow and volumes from phase-contrast cardiac magnetic resonance imaging. A mixed-linear model applied to 37,653 individuals from the UK Biobank reveals genome-wide significant associations across cardiac dynamic flow volumes spanning from aortic forward velocity to aortic regurgitation fraction. Mendelian randomization reveals a causal role for aortic root size in aortic valve regurgitation. Among the most significant contributing variants, localizing genes (near ELN, PRDM6 and ADAMTS7) are implicated in connective tissue and blood pressure pathways. Here we show that DeepFlow cardiac flow phenotyping at scale, combined with genotyping data, reinforces the contribution of connective tissue genes, blood pressure and root size to aortic valve function.

Genetic Variants of Adrenoceptors

Adrenoceptors are class A G-protein-coupled receptors grouped into three families (α1-, α2-, and β-adrenoceptors), each one including three members. All nine corresponding adrenoceptor genes display genetic variation in their coding and adjacent non-coding genomic region. Coding variants, i.e., nucleotide exchanges within the transcribed and translated receptor sequence, may result in a difference in amino acid sequence thus altering receptor function and signaling. Such variants have been intensely studied in vitro in overexpression systems and addressed in candidate-gene studies for distinct clinical parameters. In recent years, large cohorts were analyzed in genome-wide association studies (GWAS), where variants are detected as significant in context with specific traits. These studies identified two of the in-depth characterized 18 coding variants in adrenoceptors as repeatedly statistically significant genetic risk factors - p.Arg389Gly in the β1- and p.Thr164Ile in the β2-adrenoceptor, along with 56 variants in the non-coding regions adjacent to the adrenoceptor gene loci, the functional role of which is largely unknown at present. This chapter summarizes current knowledge on the two coding variants in adrenoceptors that have been consistently validated in GWAS and provides a prospective overview on the numerous non-coding variants more recently attributed to adrenoceptor gene loci.

The role of depression and antidepressant treatment in antihypertensive medication adherence and persistence: Utilising electronic health record data

Higher blood pressure levels in patients with depression may be associated with lower adherence to antihypertensive medications (AHMs). Here, we use electronic health record (EHR) data from the Estonian Biobank (EstBB) to investigate the role of lifetime depression in AHM adherence and persistence. We also explore the relationship between antidepressant initiation and intraindividual change in AHM adherence among hypertension (HTN) patients with newly diagnosed depression. Diagnosis and pharmacy refill data were obtained from the National Health Insurance database. Adherence and persistence to AHMs were determined for hypertension (HTN) patients initiating treatment between 2009 and 2017 with a three-year follow-up period. Multivariable regression was used to explore the associations between depression and AHM adherence or persistence, adjusting for sociodemographic, genetic, and health-related factors. A linear mixed-effects model was used to estimate the effect of antidepressant treatment initiation on antihypertensive medication adherence, adjusting for age and sex. We identified 20,724 individuals with newly diagnosed HTN (6294 depression cases and 14,430 controls). Depression was associated with 6% lower probability of AHM adherence (OR = 0.943, 95%CI = 0.909-0.979) and 12% lower odds of AHM persistence (OR = 0.876, 95%CI = 0.821-0.936). Adjusting for sociodemographic, genetic, and health-related factors did not significantly influence these associations. AHM adherence increased 8% six months after initiating antidepressant therapy (N = 132; β = 0.078; 95%CI = 0.025-0.131). Based on the EHR data on EstBB participants, depression is associated with lower AHM adherence and persistence. Additionally, antidepressant therapy may help improve AHM adherence in patients with depression.

Cardiac Molecular Analysis Reveals Aging-Associated Metabolic Alterations Promoting Glycosaminoglycans Accumulation Via Hexosamine Biosynthetic Pathway

Age is a prominent risk factor for cardiometabolic disease, and often leads to heart structural and functional changes. However, precise molecular mechanisms underlying cardiac remodeling and dysfunction resulting from physiological aging per se remain elusive. Understanding these mechanisms requires biological models with optimal translation to humans. Previous research demonstrated that baboons undergo age-related reduction in ejection fraction and increased heart sphericity, mirroring changes observed in humans. The goal of this study was to identify early cardiac molecular alterations that precede functional adaptations, shedding light on the regulation of age-associated changes. We performed unbiased transcriptomics of left ventricle (LV) samples from female baboons aged 7.5-22.1 years (human equivalent ~30-88 years). Weighted-gene correlation network and pathway enrichment analyses were performed to identify potential age-associated mechanisms in LV, with histological validation. Myocardial modules of transcripts negatively associated with age were primarily enriched for cardiac metabolism, including oxidative phosphorylation, tricarboxylic acid cycle, glycolysis, and fatty-acid β-oxidation. Transcripts positively correlated with age suggest upregulation of glucose uptake, pentose phosphate pathway, and hexosamine biosynthetic pathway (HBP), indicating a metabolic shift towards glucose-dependent anabolic pathways. Upregulation of HBP commonly results in increased glycosaminoglycan precursor synthesis. Transcripts involved in glycosaminoglycan synthesis, modification, and intermediate metabolism were also upregulated in older animals, while glycosaminoglycan degradation transcripts were downregulated with age. These alterations would promote glycosaminoglycan accumulation, which was verified histologically. Upregulation of extracellular matrix (ECM)-induced signaling pathways temporally coincided with glycosaminoglycan accumulation. We found a subsequent upregulation of cardiac hypertrophy-related pathways and an increase in cardiomyocyte width. Overall, our findings revealed a transcriptional shift in metabolism from catabolic to anabolic pathways that leads to ECM glycosaminoglycan accumulation through HBP prior to upregulation of transcripts of cardiac hypertrophy-related pathways. This study illuminates cellular mechanisms that precede development of cardiac hypertrophy, providing novel potential targets to remediate age-related cardiac diseases.

Genetically Predicted Vegetable Intake and Cardiovascular Diseases and Risk Factors: An Investigation with Mendelian Randomization

BACKGROUND

The associations between vegetable intake and cardiovascular diseases have been demonstrated in observational studies, but less sufficiently in randomized trials. Mendelian randomization has been considered a promising alternative in causal inference. The separate effects of cooked and raw vegetable intake remain unclear. This study aimed to investigate the associations between cooked and raw vegetable intake with cardiovascular outcomes using MR.

METHODS

We identified 15 and 28 genetic variants statistically and biologically associated with cooked and raw vegetable intake, respectively, from previous genome-wide association studies, which were used as instrumental variables to estimate associations with coronary heart disease (CHD), stroke, heart failure (HF), and atrial fibrillation (AF). The independent effects of genetically predicted cooked and raw vegetable intake were examined using multivariable MR analysis. We performed one-sample and two-sample MR analyses and combined their results using meta-analysis. Bonferroni correction was applied for multiple comparisons. We performed two-sample MR analysis for cardiometabolic risk factors (serum lipids, blood pressure, body mass index, and glycemic traits) to explore the potential mechanisms.

RESULTS

In the MR meta-analysis of 1.2 million participants, we found null evidence for associations between genetically predicted cooked and raw vegetable intake with CHD, HF, or AF. Raw vegetable intake was nominally associated with stroke (odds ratio [95% confidence interval] 0.82 [0.69-0.98] per 1 daily serving increase, p = 0.03), but this association did not pass the corrected significance level. We found consistently null evidence for associations with serum lipids, blood pressure, body mass index, or glycemic traits.

CONCLUSIONS

We found null evidence for associations between genetically predicted vegetable intake with CHD, AF, HF, or cardiometabolic risk factors in this MR study. Raw vegetable intake may reduce risk of stroke, but this warrants more research. True associations between vegetable intake and CVDs cannot be completely ruled out, and future investigations are required for causal inference in nutritional research.

Systematic characterization of regulatory variants of blood pressure genes

High blood pressure (BP) is the major risk factor for cardiovascular disease. Genome-wide association studies have identified genetic variants for BP, but functional insights into causality and related molecular mechanisms lag behind. We functionally characterize 4,608 genetic variants in linkage with 135 BP loci in vascular smooth muscle cells and cardiomyocytes by massively parallel reporter assays. High densities of regulatory variants at BP loci (i.e., ULK4, MAP4, CFDP1, PDE5A) indicate that multiple variants drive genetic association. Regulatory variants are enriched in repeats, alter cardiovascular-related transcription factor motifs, and spatially converge with genes controlling specific cardiovascular pathways. Using heuristic scoring, we define likely causal variants, and CRISPR prime editing finally determines causal variants for KCNK9, SFXN2, and PCGF6, which are candidates for developing high BP. Our systems-level approach provides a catalog of functionally relevant variants and their genomic architecture in two trait-relevant cell lines for a better understanding of BP gene regulation.

Evaluating the use of blood pressure polygenic risk scores across race/ethnic background groups

We assess performance and limitations of polygenic risk scores (PRSs) for multiple blood pressure (BP) phenotypes in diverse population groups. We compare "clumping-and-thresholding" (PRSice2) and LD-based (LDPred2) methods to construct PRSs from each of multiple GWAS, as well as multi-PRS approaches that sum PRSs with and without weights, including PRS-CSx. We use datasets from the MGB Biobank, TOPMed study, UK biobank, and from All of Us to train, assess, and validate PRSs in groups defined by self-reported race/ethnic background (Asian, Black, Hispanic/Latino, and White). For both SBP and DBP, the PRS-CSx based PRS, constructed as a weighted sum of PRSs developed from multiple independent GWAS, perform best across all race/ethnic backgrounds. Stratified analysis in All of Us shows that PRSs are better predictive of BP in females compared to males, individuals without obesity, and middle-aged (40-60 years) compared to older and younger individuals.

Development of a genetic risk score to predict the risk of hypertension in European adolescents from the HELENA study

Introduction

From genome wide association study (GWAS) a large number of single nucleotide polymorphisms (SNPs) have previously been associated with blood pressure (BP) levels. A combination of SNPs, forming a genetic risk score (GRS) could be considered as a useful genetic tool to identify individuals at risk of developing hypertension from early stages in life. Therefore, the aim of our study was to build a GRS being able to predict the genetic predisposition to hypertension (HTN) in European adolescents.

Methods

Data were extracted from the Healthy Lifestyle in Europe by Nutrition in Adolescence (HELENA) cross-sectional study. A total of 869 adolescents (53% female), aged 12.5-17.5, with complete genetic and BP information were included. The sample was divided into altered (≥130 mmHg for systolic and/or ≥80 mmHg for diastolic) or normal BP. Based on the literature, a total of 1.534 SNPs from 57 candidate genes related with BP were selected from the HELENA GWAS database.

Results

From 1,534 SNPs available, An initial screening of SNPs univariately associated with HTN (p < 0.10) was established, to finally obtain a number of 16 SNPs significantly associated with HTN (p < 0.05) in the multivariate model. The unweighted GRS (uGRS) and weighted GRS (wGRS) were estimated. To validate the GRSs, the area under the curve (AUC) was explored using ten-fold internal cross-validation for uGRS (0.802) and wGRS (0.777). Further covariates of interest were added to the analyses, obtaining a higher predictive ability (AUC values of uGRS: 0.879; wGRS: 0.881 for BMI z-score). Furthermore, the differences between AUCs obtained with and without the addition of covariates were statistically significant (p < 0.05).

Conclusions

Both GRSs, the uGRS and wGRS, could be useful to evaluate the predisposition to hypertension in European adolescents.

Risk Effects of rs1799945 Polymorphism of the HFE Gene and Intergenic Interactions of GWAS-Significant Loci for Arterial Hypertension in the Caucasian Population of Central Russia

The aim of this case-control replicative study was to investigate the link between GWAS-impact for arterial hypertension (AH) and/or blood pressure (BP) gene polymorphisms and AH risk in Russian subjects (Caucasian population of Central Russia). AH (n = 939) and control (n = 466) cohorts were examined for ten GWAS AH/BP risk loci. The genotypes/alleles of these SNP and their combinations (SNP-SNP interactions) were tested for their association with the AH development using a logistic regression statistical procedure. The genotype GG of the SNP rs1799945 (C/G) HFE was strongly linked with an increased AH risk (ORrecGG = 2.53; 95%CIrecGG1.03-6.23; ppermGG = 0.045). The seven SNPs such as rs1173771 (G/A) AC026703.1, rs1799945 (C/G) HFE, rs805303 (G/A) BAG6, rs932764 (A/G) PLCE1, rs4387287 (C/A) OBFC1, rs7302981 (G/A) CERS5, rs167479 (T/G) RGL3, out of ten regarded loci, were related with AH within eight SNP-SNP interaction models (<0.001 ≤ pperm-interaction ≤ 0.047). Three polymorphisms such as rs8068318 (T/C) TBX2, rs633185 (C/G) ARHGAP42, and rs2681472 (A/G) ATP2B1 were not linked with AH. The pairwise rs805303 (G/A) BAG6-rs7302981 (G/A) CERS5 combination was a priority in determining the susceptibility to AH (included in six out of eight SNP-SNP interaction models [75%] and described 0.82% AH entropy). AH-associated variants are conjecturally functional for 101 genes involved in processes related to the immune system (major histocompatibility complex protein, processing/presentation of antigens, immune system process regulation, etc.). In conclusion, the rs1799945 polymorphism of the HFE gene and intergenic interactions of BAG6, CERS5, AC026703.1, HFE, PLCE1, OBFC1, RGL3 have been linked with AH risky in the Caucasian population of Central Russia.

Sex-Specific Features of the Correlation between GWAS-Noticeable Polymorphisms and Hypertension in Europeans of Russia

The aim of the study was directed at studying the sex-specific features of the correlation between genome-wide association studies (GWAS)-noticeable polymorphisms and hypertension (HTN). In two groups of European subjects of Russia (n = 1405 in total), such as men (n = 821 in total: n = 564 HTN, n = 257 control) and women (n = 584 in total: n = 375 HTN, n = 209 control), the distribution of ten specially selected polymorphisms (they have confirmed associations of GWAS level with blood pressure (BP) parameters and/or HTN in Europeans) has been considered. The list of studied loci was as follows: (PLCE1) rs932764 A > G, (AC026703.1) rs1173771 G > A, (CERS5) rs7302981 G > A, (HFE) rs1799945 C > G, (OBFC1) rs4387287 C > A, (BAG6) rs805303 G > A, (RGL3) rs167479 T > G, (ARHGAP42) rs633185 C > G, (TBX2) rs8068318 T > C, and (ATP2B1) rs2681472 A > G. The contribution of individual loci and their inter-locus interactions to the HTN susceptibility with bioinformatic interpretation of associative links was evaluated separately in men's and women's cohorts. The men-women differences in involvement in the disease of the BP/HTN-associated GWAS SNPs were detected. Among women, the HTN risk has been associated with HFE rs1799945 C > G (genotype GG was risky; OR_GG = 11.15 p_permGG = 0.014) and inter-locus interactions of all 10 examined SNPs as part of 26 intergenic interactions models. In men, the polymorphism BAG6 rs805303 G > A (genotype AA was protective; OR_AA = 0.30 p_permAA = 0.0008) and inter-SNPs interactions of eight loci in only seven models have been founded as HTN-correlated. HTN-linked loci and strongly linked SNPs were characterized by pronounced polyvector functionality in both men and women, but at the same time, signaling pathways of HTN-linked genes/SNPs in women and men were similar and were represented mainly by immune mechanisms. As a result, the present study has demonstrated a more pronounced contribution of BP/HTN-associated GWAS SNPs to the HTN susceptibility (due to weightier intergenic interactions) in European women than in men.

Polygenic Susceptibility to Hypertension and Blood Pressure Control in Stroke Survivors

BACKGROUND AND OBJECTIVES

Blood pressure (BP) is often not at goal in stroke survivors, leaving individuals vulnerable to additional vascular events. Given that BP is a highly heritable trait, we hypothesize that a higher polygenic susceptibility to hypertension (PSH) leads to worse BP control in stroke survivors.

METHODS

We conducted a study within the UK Biobank evaluating persons of European ancestry who survived an ischemic or hemorrhagic stroke. To model the PSH, we created polygenic risk scores (PRSs) for systolic and diastolic BP using 732 genetic variants. We divided the PRSs into quintiles and used linear/logistic regression to test whether higher PSH led to higher observed BP, uncontrolled BP (systolic BP > 140 mm Hg or diastolic BP > 90 mm Hg), and resistant BP (uncontrolled BP despite being on ≥3 antihypertensive drugs). We conducted an independent replication using data from the Vitamin Intervention for Stroke Prevention (VISP) trial.

RESULTS

We analyzed 5,940 stroke survivors. When comparing stroke survivors with very low vs very high PSH, the mean systolic BP was 137 (SD 18) vs 143 (SD 20, p < 0.001), the mean diastolic BP was 81 (SD 10) vs 84 (SD 11, p < 0.001), the prevalence of uncontrolled BP was 42.8% vs 57.2% (p < 0.001), and the prevalence of resistant hypertension was 3.9% vs 11% (p < 0.001). Results remained significant using multivariable models (p < 0.001) and were replicated in the VISP study (all tests with p < 0.05).

DISCUSSION

A higher PSH is associated with worse BP control in stroke survivors. These findings point to genetic predisposition as an important determinant of poorly controlled BP in this population.

A Genomic Risk Score Identifies Individuals at High Risk for Intracerebral Hemorrhage

BACKGROUND

Intracerebral hemorrhage (ICH) has an estimated heritability of 29%. We developed a genomic risk score for ICH and determined its predictive power in comparison to standard clinical risk factors.

METHODS

We combined genome-wide association data from individuals of European ancestry for ICH and related traits in a meta-genomic risk score ([metaGRS]; 2.6 million variants). We tested associations with ICH and its predictive performance in addition to clinical risk factors in a held-out validation dataset (842 cases and 796 controls). We tested associations with risk of incident ICH in the population-based UK Biobank cohort (486 784 individuals, 1526 events, median follow-up 11.3 years).

RESULTS

One SD increment in the metaGRS was significantly associated with 31% higher odds for ICH (95% CI, 1.16-1.48) in age-, sex- and clinical risk factor-adjusted models. The metaGRS identified individuals with almost 5-fold higher odds for ICH in the top score percentile (odds ratio, 4.83 [95% CI, 1.56-21.2]). Predictive models for ICH incorporating the metaGRS in addition to clinical predictors showed superior performance compared to the clinical risk factors alone (c-index, 0.695 versus 0.686). The metaGRS showed similar associations for lobar and nonlobar ICH, independent of the known APOE risk locus for lobar ICH. In the UK Biobank, the metaGRS was associated with higher risk of incident ICH (hazard ratio, 1.15 [95% CI, 1.09-1.21]). The associations were significant within both a relatively high-risk population of antithrombotic medications users, as well as among a relatively low-risk population with a good control of vascular risk factors and no use of anticoagulants.

CONCLUSIONS

We developed and validated a genomic risk score that predicts lifetime risk of ICH beyond established clinical risk factors among individuals of European ancestry. Whether implementation of the score in risk prognostication models for high-risk populations, such as patients under antithrombotic treatment, could improve clinical decision making should be explored in future studies.

Genome-wide association study of early ischaemic stroke risk in Brazilian individuals with sickle cell disease implicates ADAMTS2 and CDK18 and uncovers novel loci

Ischaemic stroke is a common complication of sickle cell disease (SCD) and without intervention can affect 11% of children with SCD before the age of 20. Within the Trans-Omics for Precision Medicine (TOPMed), a genome-wide association study (GWAS) of ischaemic stroke was performed on 1333 individuals with SCD from Brazil (178 cases, 1155 controls). Via a novel Cox proportional-hazards analysis, we searched for variants associated with ischaemic stroke occurring at younger ages. Variants at genome-wide significance (p < 5 × 10-8 ) include two near genes previously linked to non-SCD early-onset stroke (<65 years): ADAMTS2 (rs147625068, p = 3.70 × 10-9 ) and CDK18 (rs12144136, p = 2.38 × 10-9 ). Meta-analysis, which included the independent SCD cohorts Walk-PHaSST and PUSH, exhibited consistent association for variants rs1209987 near gene TBC1D32 (p = 3.36 × 10-10 ), rs188599171 near CUX1 (p = 5.89 × 10-11 ), rs77900855 near BTG1 (p = 4.66 × 10-8 ), and rs141674494 near VPS13C (1.68 × 10-9 ). Findings from this study support a multivariant model of early ischaemic stroke risk and possibly a shared genetic architecture between SCD individuals and non-SCD individuals younger than 65 years.

A systems biology approach identifies the role of dysregulated PRDM6 in the development of hypertension

Genetic variants in the third intron of the PRDM6 gene have been associated with BP traits in multiple GWAS. By combining fine mapping, massively parallel reporter assays, and gene editing, we identified super enhancers that drive the expression of PRDM6 and are partly regulated by STAT1 as the causal variants for hypertension. The heterozygous disruption of Prdm6 in mice expressing Cre recombinase under the control of mouse smooth muscle cell protein 22-α promoter (Prdm6fl/+ SM22-Cre) exhibited a markedly higher number of renin-producing cells in the kidneys at E18.5 compared with WT littermates and developed salt-induced systemic hypertension that was completely responsive to the renin inhibitor aliskiren. Strikingly, RNA-Seq analysis of the mouse aortas identified a network of PRDM6-regulated genes that are located in GWAS-associated loci for blood pressure, most notably Sox6, which modulates renin expression in the kidney. Accordingly, the smooth muscle cell-specific disruption of Sox6 in Prdm6fl/+ SM22-Cre mice resulted in a dramatic reduction of renin. Fate mapping and histological studies also showed increased numbers of neural crest-derived cells accompanied by increased collagen deposition in the kidneys of Prdm6fl/+ Wnt1Cre-ZsGreen1Cre mice compared with WT mice. These findings establish the role of PRDM6 as a regulator of renin-producing cell differentiation into smooth muscle cells and as an attractive target for the development of antihypertensive drugs.

Haplotypic variants of COVID-19 related genes are associated with blood pressure and metabolites levels

The genetic association of coronavirus disease 2019 (COVID-19) with its complications has not been fully understood. This study aimed to identify variants and haplotypes of candidate genes implicated in COVID-19 related traits by combining the literature review and pathway analysis. To explore such genes, the protein-protein interactions and relevant pathways of COVID-19-associated genes were assessed. A number of variants on candidate genes were identified from Genome-wide association studies (GWASs) which were associated with COVID-19 related traits (p ˂ 10^-6 ). Haplotypic blocks were assessed using haplotypic structures among the 1000 Genomes Project (r²  ≥ 0.8, D' ≥ 0.8). Further functional analyses were performed on the selected variants. The results demonstrated that a group of variants in ACE and AGT genes were significantly correlated with COVID-19 related traits. Three haplotypes were identified to be involved in the blood metabolites levels and the development of blood pressure. Functional analyses revealed that most GWAS index variants were expression quantitative trait loci and had transcription factor binding sites, exonic splicing enhancers or silencer activities. Furthermore, the proxy haplotype variants, rs4316, rs4353, rs4359, and three variants, namely rs2493133, rs2478543, and rs5051, were associated with blood metabolite and systolic blood pressure, respectively. These variants exerted more regulatory effects compared with other GWAS variants. The present study indicates that the genetic variants and candidate haplotypes of COVID-19 related genes are associated with blood pressure and blood metabolites. However, further observational studies are warranted to confirm these results.

Effects of Pre-Pregnancy Overweight/Obesity on the Pattern of Association of Hypertension Susceptibility Genes with Preeclampsia

The aim of this study was to explore the effects of pre-pregnancy overweight/obesity on the pattern of association of hypertension susceptibility genes with preeclampsia (PE). Ten single-nucleotide polymorphisms (SNPs) of the 10 genome-wide association studies (GWAS)-significant hypertension/blood pressure (BP) candidate genes were genotyped in 950 pregnant women divided into two cohorts according to their pre-pregnancy body mass index (preBMI): preBMI ≥ 25 (162 with PE and 159 control) and preBMI < 25 (290 with PE and 339 control). The PLINK software package was utilized to study the association (analyzed four genetic models using logistic regression). The functionality of PE-correlated loci was analyzed by performing an in silico database analysis. Two SNP hypertension/BP genes, rs805303 BAG6 (OR: 0.36−0.66) and rs167479 RGL3 (OR: 1.86), in subjects with preBMI ≥ 25 were associated with PE. No association between the studied SNPs and PE in the preBMI < 25 group was determined. Further analysis showed that two PE-associated SNPs are functional (have weighty eQTL, sQTL, regulatory, and missense values) and could be potentially implicated in PE development. In conclusion, this study was the first to discover the modifying influence of overweight/obesity on the pattern of association of GWAS-significant hypertension/BP susceptibility genes with PE: these genes are linked with PE in preBMI ≥ 25 pregnant women and are not PE-involved in the preBMI < 25 group.

Polymorphisms of hypertension susceptibility genes as a risk factors of preeclampsia in the Caucasian population of central Russia

INTRODUCTION

The study was designed to assess the effects of hypertension (HT) susceptibility genes polymorphisms in the development of preeclampsia (PE) in Caucasians from Central Russia.

METHODS

PE patients (n = 452) and women control group (n = 498) were genotyped for 10 polymorphisms of HT/blood pressure (BP) susceptibility genes (according to the previously published GWAS in Caucasian populations) including AC026703.1 (rs1173771), HFE (rs1799945), BAG6 (rs805303), PLCE1 (rs932764), OBFC1 (rs4387287), ARHGAP42 (rs633185), CERS5 (rs7302981), ATP2B1 (rs2681472), TBX2 (rs8068318) and RGL3 (rs167479). A logistic regression method was applied to search for associations between SNPs and PE. The relationship between SNP-SNP interactions and PE risk was analyzed by performing MB-MDR.

RESULTS

The rs1799945 gene in HFE significantly independently increased the risk of developing PE (OR = 2.24) and rs805303 in BAG6 was associated with a reduced risk in the occurrence of PE (OR = 0.55-0.78). Among the 10 SNPs examined, nine SNPs were associated with PEs within the 10 most significant SNP-SNP interaction models. Loci rs7302981 CERS5, rs805303 BAG6 and rs932764 PLCE1 contributed to the largest number of epistatic models (50% or more).

DISCUSSION

The present study is the first to report an association between polymorphisms of HT/BP susceptibility genes important for GWAS and the risk of PE in Caucasians from Central Russia. Our pathway-based functional annotation of the PE risk variants highlights the potential regulatory function (epigenetic/eQTL/sQTL/non-synonymous) that nine genetic risk markers and their 115 highly correlated variants exert on 155 genes. The study shows that these genes may function cooperatively in key signaling pathways in PE biology.

Genomic basis of insularity and ecological divergence in barn owls (Tyto alba) of the Canary Islands

Islands, and the particular organisms that populate them, have long fascinated biologists. Due to their isolation, islands offer unique opportunities to study the effect of neutral and adaptive mechanisms in determining genomic and phenotypical divergence. In the Canary Islands, an archipelago rich in endemics, the barn owl (Tyto alba), present in all the islands, is thought to have diverged into a subspecies (T. a. gracilirostris) on the eastern ones, Fuerteventura and Lanzarote. Taking advantage of 40 whole-genomes and modern population genomics tools, we provide the first look at the origin and genetic makeup of barn owls of this archipelago. We show that the Canaries hold diverse, long-standing and monophyletic populations with a neat distinction of gene pools from the different islands. Using a new method, less sensitive to structure than classical F_ST, to detect regions involved in local adaptation to insular environments, we identified a haplotype-like region likely under selection in all Canaries individuals and genes in this region suggest morphological adaptations to insularity. In the eastern islands, where the subspecies is present, genomic traces of selection pinpoint signs of adapted body proportions and blood pressure, consistent with the smaller size of this population living in a hot arid climate. In turn, genomic regions under selection in the western barn owls from Tenerife showed an enrichment in genes linked to hypoxia, a potential response to inhabiting a small island with a marked altitudinal gradient. Our results illustrate the interplay of neutral and adaptive forces in shaping divergence and early onset speciation.

Kidney omics in hypertension: from statistical associations to biological mechanisms and clinical applications

Hypertension is a major cardiovascular disease risk factor and contributor to premature death globally. Family-based investigations confirmed a significant heritable component of blood pressure (BP), whereas genome-wide association studies revealed >1000 common and rare genetic variants associated with BP and/or hypertension. The kidney is not only an organ of key relevance to BP regulation and the development of hypertension, but it also acts as the tissue mediator of genetic predisposition to hypertension. The identity of kidney genes, pathways, and related mechanisms underlying the genetic associations with BP has started to emerge through integration of genomics with kidney transcriptomics, epigenomics, and other omics as well as through applications of causal inference, such as Mendelian randomization. Single-cell methods further enabled mapping of BP-associated kidney genes to cell types, and in conjunction with other omics, started to illuminate the biological mechanisms underpinning associations of BP-associated genetic variants and kidney genes. Polygenic risk scores derived from genome-wide association studies and refined on kidney omics hold the promise of enhanced diagnostic prediction, whereas kidney omics-informed drug discovery is likely to contribute new therapeutic opportunities for hypertension and hypertension-mediated kidney damage.

Identification and single-base gene-editing functional validation of a cis-EPO variant as a genetic predictor for EPO-increasing therapies

Hypoxia-inducible factor prolyl hydroxylase inhibitors (HIF-PHIs) are currently under clinical development for treating anemia in chronic kidney disease (CKD), but it is important to monitor their cardiovascular safety. Genetic variants can be used as predictors to help inform the potential risk of adverse effects associated with drug treatments. We therefore aimed to use human genetics to help assess the risk of adverse cardiovascular events associated with therapeutically altered EPO levels to help inform clinical trials studying the safety of HIF-PHIs. By performing a genome-wide association meta-analysis of EPO (n = 6,127), we identified a cis-EPO variant (rs1617640) lying in the EPO promoter region. We validated this variant as most likely causal in controlling EPO levels by using genetic and functional approaches, including single-base gene editing. Using this variant as a partial predictor for therapeutic modulation of EPO and large genome-wide association data in Mendelian randomization tests, we found no evidence (at p < 0.05) that genetically predicted long-term rises in endogenous EPO, equivalent to a 2.2-unit increase, increased risk of coronary artery disease (CAD, OR [95% CI] = 1.01 [0.93, 1.07]), myocardial infarction (MI, OR [95% CI] = 0.99 [0.87, 1.15]), or stroke (OR [95% CI] = 0.97 [0.87, 1.07]). We could exclude increased odds of 1.15 for cardiovascular disease for a 2.2-unit EPO increase. A combination of genetic and functional studies provides a powerful approach to investigate the potential therapeutic profile of EPO-increasing therapies for treating anemia in CKD.

Appraising the Causal Association between Systemic Iron Status and Heart Failure Risk: A Mendelian Randomisation Study

Although observational studies have shown that abnormal systemic iron status is associated with an increased risk of heart failure (HF), it remains unclear whether this relationship represents true causality. We aimed to explore the causal relationship between iron status and HF risk. Two-sample Mendelian randomisation (MR) was applied to obtain a causal estimate. Genetic summary statistical data for the associations (p < 5 × 10−8) between single nucleotide polymorphisms (SNPs) and four iron status parameters were obtained from the Genetics of Iron Status Consortium in genome-wide association studies involving 48,972 subjects. Statistical data on the association of SNPs with HF were extracted from the UK biobank consortium (including 1088 HF cases and 360,106 controls). The results were further tested using MR based on the Bayesian model averaging (MR-BMA) and multivariate MR (MVMR). Of the twelve SNPs considered to be valid instrumental variables, three SNPs (rs1800562, rs855791, and rs1799945) were associated with all four iron biomarkers. Genetically predicted iron status biomarkers were not causally associated with HF risk (all p > 0.05). Sensitivity analysis did not show evidence of potential heterogeneity and horizontal pleiotropy. Convincing evidence to support a causal relationship between iron status and HF risk was not found. The strong relationship between abnormal iron status and HF risk may be explained by an indirect mechanism.

Large-Scale Multi-Omics Studies Provide New Insights into Blood Pressure Regulation

Recent genome-wide association studies uncovered part of blood pressure’s heritability. However, there is still a vast gap between genetics and biology that needs to be bridged. Here, we followed up blood pressure genome-wide summary statistics of over 750,000 individuals, leveraging comprehensive epigenomic and transcriptomic data from blood with a follow-up in cardiovascular tissues to prioritise likely causal genes and underlying blood pressure mechanisms. We first prioritised genes based on coding consequences, multilayer molecular associations, blood pressure-associated expression levels, and coregulation evidence. Next, we followed up the prioritised genes in multilayer studies of genomics, epigenomics, and transcriptomics, functional enrichment, and their potential suitability as drug targets. Our analyses yielded 1880 likely causal genes for blood pressure, tens of which are targets of the available licensed drugs. We identified 34 novel genes for blood pressure, supported by more than one source of biological evidence. Twenty-eight (82%) of these new genes were successfully replicated by transcriptome-wide association analyses in a large independent cohort (n = ~220,000). We also found a substantial mediating role for epigenetic regulation of the prioritised genes. Our results provide new insights into genetic regulation of blood pressure in terms of likely causal genes and involved biological pathways offering opportunities for future translation into clinical practice.

Genetic analyses of circulating PUFA-derived mediators identifies heritable dihydroxyeicosatrienoic acid species

Estimates of heritability are the first step in identifying a trait with substantial variation due to genetic factors. Large-scale genetic analyses can identify the DNA variants that influence the levels of circulating lipid species and the statistical technique Mendelian randomisation can use these DNA variants to address potential causality of these lipids in disease. We estimated the heritability of plasma eicosanoids, octadecanoids and docosanoids to identify those lipid species with substantial heritability. We analysed plasma lipid mediators in 31 White British families (196 participants) ascertained for high blood pressure and deeply clinically and biochemically phenotyped over a 25-year period. We found that the dihydroxyeicosatrienoic acid (DHET) species, 11,12-DHET and 14,15-DHET, products of arachidonic acid metabolism by cytochrome P450 (CYP) monooxygenase and soluble epoxide hydrolase (sEH), exhibited substantial heritability (h2 = 33%-37%; Padj<0.05). Identification of these two heritable bioactive lipid species allows for future large-scale, targeted, lipidomics-genomics analyses to address causality in cardiovascular and other diseases.

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.

Roles of Cardiometabolic Factors in Mediating the Causal Effect of Type 2 Diabetes on Cardiovascular Diseases: A Two-Step, Two-Sample Multivariable Mendelian Randomization Study

Objective

The objective of this study is to investigate the roles of cardiometabolic factors (including blood pressure, blood lipids, thyroid function, body mass, and insulin sensitivity) in mediating the causal effect of type 2 diabetes (T2DM) on cardiovascular disease (CVD) outcomes.

Design

Two-step, two-sample multivariable Mendelian randomization (MVMR) study.

Setting

International genome-wide association study (GWAS) consortia data.

Exposure

Type 2 diabetes, blood pressure: systolic blood pressure (SBP), diastolic blood pressure (DBP); blood lipids: low-density lipoprotein (LDL), high-density lipoprotein (HDL), total cholesterol (TC), triglycerides (TG); thyroid function: hyperthyroidism, hypothyroidism; body mass index (BMI), waist-hip-ratio (WHR), and insulin sensitivity.

Main Outcomes

Cardiovascular disease includes coronary heart disease (CHD), myocardial infarction (MI), and stroke.

Methods

Summary-level data for exposures and main outcomes were extracted from GWAS consortia. We used two-sample MR to illustrate the causal effect of T2DM on CVD subtypes and regression-based MVMR to quantify the possible mediation effects of cardiometabolic factors on CVD.

Results

Each additional unit of log odds of T2DM increased 16% risk of CHD [odds ratio (OR): 1.16, 95% CI: 1.12–1.21], 15% risk of myocardial infarction (MI) (OR: 1.15, 95% CI: 1.10–1.20), and 10% risk of stroke (OR: 1.10, 95% CI: 1.06–1.13). In mediation analysis, SBP, DBP, and TG were found as main mediators, while the mediation effects of other cardiometabolic factors were not significant. The proportion of total effect of T2DM on CHD mediated by SBP, DBP, and TG was 16% (95% CI: 8–24%), 7% (95% CI: 1–13%) and 10% (95% CI: 2–18%), respectively. Mediation effect of SBP and DBP on MI and stroke, TG on MI was also prominent, while mediation effect of TG on stroke was not significant. The combined mediation effect of all the three mediators accounted for 29%, 26%, and 13% of the total effect of T2DM on CHD, MI, and stroke, respectively.

Conclusion

Systolic blood pressure, DBP, and TG mediate a substantial proportion of the causal effect of T2DM on CVD and thus interventions on these factors might reduce the considerable excess risk of CVD among patients with T2DM.

Promoter DNA Methylation in GWAS-Identified Genes as Potential Functional Elements for Blood Pressure: An Observational and Mendelian Randomization Study

Genome-wide association studies have identified numerous genetic loci for blood pressure (BP). However, the relationships of functional elements inside these loci with BP are not fully understood. This study represented an effort to determine if promoter DNA methylations inside BP-associated loci were associated with BP.We conducted a cross-sectional study investigating the association between promoter DNA methylations of 10 candidate genes and BP in 1,241 Chinese individuals. Twenty-one genomic fragments in the CpG Islands were sequenced. The associations of methylation levels with BP and hypertension were assessed in regression models. Mendelian randomization (MR) analysis was then applied to find supporting evidence for the identified associations.A total of 413 DNA methylation sites were examined in an observational study. Methylation levels of 24 sites in PRDM6, IGFBP3, SYT7, PDE3A, TBX2 and C17orf82 were significantly associated with BP. Methylation levels of PRDM6 and SYT7 were significantly associated with hypertension. Methylation levels of five sites (including cg06713098) in IGFBP3 were significantly associated with DBP. MR analysis found associations between the methylation levels of six CpG sites (cg06713098, cg14228300, cg23193639, cg21268650, cg10677697 and cg04812164) around the IGFBP3 promoter and DBP. Methylation levels of cg14228300 and cg04812164 were associated with SBP. By further applying several MR methods we showed that the associations may not be due to pleiotropy. Association between IGFBP3 mRNA levels in blood cells and BP was also found in MR analysis. This study identified promoter methylation as potential functional element for BP. The identified methylations may be involved in the regulatory pathway linking genetic variants to BP.

Deep learning enables genetic analysis of the human thoracic aorta

Enlargement or aneurysm of the aorta predisposes to dissection, an important cause of sudden death. We trained a deep learning model to evaluate the dimensions of the ascending and descending thoracic aorta in 4.6 million cardiac magnetic resonance images from the UK Biobank. We then conducted genome-wide association studies in 39,688 individuals, identifying 82 loci associated with ascending and 47 with descending thoracic aortic diameter, of which 14 loci overlapped. Transcriptome-wide analyses, rare-variant burden tests and human aortic single nucleus RNA sequencing prioritized genes including SVIL, which was strongly associated with descending aortic diameter. A polygenic score for ascending aortic diameter was associated with thoracic aortic aneurysm in 385,621 UK Biobank participants (hazard ratio = 1.43 per s.d., confidence interval 1.32-1.54, P = 3.3 × 10-20). Our results illustrate the potential for rapidly defining quantitative traits with deep learning, an approach that can be broadly applied to biomedical images.

Genetic Support of A Causal Relationship Between Iron Status and Type 2 Diabetes: A Mendelian Randomization Study

CONTEXT

Iron overload is a known risk factor for type 2 diabetes (T2D); however, iron overload and iron deficiency have both been associated with metabolic disorders in observational studies.

OBJECTIVE

Using mendelian randomization (MR), we assessed how genetically predicted systemic iron status affected T2D risk.

METHODS

A 2-sample MR analysis was used to obtain a causal estimate. We selected genetic variants strongly associated (P < 5 × 10-8) with 4 biomarkers of systemic iron status from a study involving 48 972 individuals performed by the Genetics of Iron Status consortium and applied these biomarkers to the T2D case-control study (74 124 cases and 824 006 controls) performed by the Diabetes Genetics Replication and Meta-analysis consortium. The simple median, weighted median, MR-Egger, MR analysis using mixture-model, weighted allele scores, and MR based on a Bayesian model averaging approaches were used for the sensitivity analysis.

RESULTS

Genetically instrumented serum iron (odds ratio [OR]: 1.07; 95% CI, 1.02-1.12), ferritin (OR: 1.19; 95% CI, 1.08-1.32), and transferrin saturation (OR: 1.06; 95% CI, 1.02-1.09) were positively associated with T2D. In contrast, genetically instrumented transferrin, a marker of reduced iron status, was inversely associated with T2D (OR: 0.91; 95% CI, 0.87-0.96).

CONCLUSION

Genetic evidence supports a causal link between increased systemic iron status and increased T2D risk. Further studies involving various ethnic backgrounds based on individual-level data and studies regarding the underlying mechanism are warranted for reducing the risk of T2D.

Search for a Functional Genetic Variant Mimicking the Effect of SGLT2 Inhibitor Treatment

SGLT2 inhibitors (SGLT2i) block renal glucose reabsorption. Due to the unexpected beneficial observations in type 2 diabetic patients potentially related to increased natriuresis, SGLT2i are also studied for heart failure treatment. This study aimed to identify genetic variants mimicking SGLT2i to further our understanding of the potential underlying biological mechanisms. Using the UK Biobank resource, we identified 264 SNPs located in the SLC5A2 gene or within 25kb of the 5′ and 3′ flanking regions, of which 91 had minor allele frequencies >1%. Twenty-seven SNPs were associated with glycated hemoglobin (HbA1c) after Bonferroni correction in participants without diabetes, while none of the SNPs were associated with sodium excretion. We investigated whether these variants had a directionally consistent effect on sodium excretion, HbA1c levels, and SLC5A2 expression. None of the variants met these criteria. Likewise, we identified no common missense variants, and although four SNPs could be defined as 5′ or 3′ prime untranslated region variants of which rs45612043 was predicted to be deleterious, these SNPs were not annotated to SLC5A2. In conclusion, no genetic variant was found mimicking SGLT2i based on their location near SLC5A2 and their association with sodium excretion or HbA1c and SLC5A2 expression or function.

Genetic correlation and causal relationships between cardio-metabolic traits and lung function impairment

BACKGROUND

Associations of low lung function with features of poor cardio-metabolic health have been reported. It is, however, unclear whether these co-morbidities reflect causal associations, shared genetic heritability or are confounded by environmental factors.

METHODS

We performed three analyses: (1) cardio-metabolic health to lung function association tests in Northern Finland Birth cohort 1966, (2) cross-trait linkage disequilibrium score regression (LDSC) to compare genetic backgrounds and (3) Mendelian randomisation (MR) analysis to assess the causal effect of cardio-metabolic traits and disease on lung function, and vice versa (bidirectional MR). Genetic associations were obtained from the UK Biobank data or published large-scale genome-wide association studies (N > 82,000).

RESULTS

We observed a negative genetic correlation between lung function and cardio-metabolic traits and diseases. In Mendelian Randomisation analysis (MR), we found associations between type 2 diabetes (T2D) instruments and forced vital capacity (FVC) as well as FEV1/FVC. Body mass index (BMI) instruments were associated to all lung function traits and C-reactive protein (CRP) instruments to FVC. These genetic associations provide evidence for a causal effect of cardio-metabolic traits on lung function. Multivariable MR suggested independence of these causal effects from other tested cardio-metabolic traits and diseases. Analysis of lung function specific SNPs revealed a potential causal effect of FEV1/FVC on blood pressure.

CONCLUSIONS

The present study overcomes many limitations of observational studies by using Mendelian Randomisation. We provide evidence for an independent causal effect of T2D, CRP and BMI on lung function with some of the T2D effect on lung function being attributed to inflammatory mechanisms. Furthermore, this analysis suggests a potential causal effect of FEV1/FVC on blood pressure. Our detailed analysis of the interplay between cardio-metabolic traits and impaired lung function provides the opportunity to improve the quality of existing intervention strategies.

Genome-wide association study on blood pressure traits in the Iranian population suggests ZBED9 as a new locus for hypertension

High blood pressure is the heritable risk factor for cardiovascular and kidney diseases. Genome-wide association studies(GWAS) on blood pressure traits increase our understanding of its underlying genetic basis. However, a large proportion of GWAS was conducted in Europeans, and some roadblocks deprive other populations to benefit from their results. Iranians population with a high degree of genomic specificity has not been represented in international databases to date, so to fill the gap, we explored the effects of 652,919 genomic variants on Systolic Blood Pressure (SBP), Diastolic Blood Pressure (DBP), and Hypertension (HTN) in 7694 Iranian adults aged 18 and over from Tehran Cardiometabolic Genetic Study (TCGS). We identified consistent signals on ZBED9 associated with HTN in the genome-wide borderline threshold after adjusting for different sets of environmental predictors. Moreover, strong signals on ABHD17C and suggestive signals on FBN1 were detected for DBP and SBP, respectively, while these signals were not consistent in different GWA analysis. Our finding on ZBED9 was confirmed for all BP traits by linkage analysis in an independent sample. We found significant associations with similar direction of effects and allele frequency of genetic variants on ZBED9 with DBP (genome-wide threshold) and HTN (nominal threshold) in GWAS summary data of UK Biobank. Although there is no strong evidence to support the function of ZBED9 in blood pressure regulation, it provides new insight into the pleiotropic effects of hypertension and other cardiovascular diseases.

Circulating ceramides as biomarkers of cardiovascular disease: Evidence from phenotypic and genomic studies

There is a need for new biomarkers of atherosclerotic cardiovascular disease (ACVD), the main cause of death globally. Ceramides, a class of potent bioactive lipid mediators, have signalling roles in apoptosis, cellular stress and inflammation. Recent studies have highlighted circulating ceramides as novel biomarkers of coronary artery disease, type-2 diabetes and insulin resistance. Ceramides are highly regulated by enzymatic reactions throughout the body in terms of their activity and metabolism, including production, degradation and transport. The genetic studies that have been completed to date on the main ceramide species found in circulation are described, highlighting the importance of DNA variants in genes involved in ceramide biosynthesis as key influencers of heritable, circulating ceramide levels. We also review studies of disease associations with ceramides and discuss mechanistic insights deriving from recent genomic studies. The signalling activities of ceramides in vascular inflammation and apoptosis, associations between circulating ceramides and coronary artery disease risk, type-2 diabetes and insulin resistance, and the potential importance of ceramides with regard to ACVD risk factors, such as blood pressure, lipoproteins and lifestyle factors, are also discussed.

Precision Medicine Approaches to Vascular Disease: JACC Focus Seminar 2/5

In this second of a 5-part Focus Seminar series, we focus on precision medicine in the context of vascular disease. The most common vascular disease worldwide is atherosclerosis, which is the primary cause of coronary artery disease, peripheral vascular disease, and a large proportion of strokes and other disorders. Atherosclerosis is a complex genetic disease that likely involves many hundreds to thousands of single nucleotide polymorphisms, each with a relatively modest effect for causing disease. Conversely, although less prevalent, there are many vascular disorders that typically involve only a single genetic change, but these changes can often have a profound effect that is sufficient to cause disease. These are termed "Mendelian vascular diseases," which include Marfan and Loeys-Dietz syndromes. Given the very different genetic basis of atherosclerosis versus Mendelian vascular diseases, this article was divided into 2 parts to cover the most promising precision medicine approaches for these disease types.

Pharmacoepidemiology: A time for a new multidisciplinary approach to precision medicine

The advent of the genomic age has created a rapid increase in complexity for the development and selection of drug treatments. A key component of precision medicine is the use of genetic information to improve therapeutic effectiveness of drugs and prevent potential adverse drug reactions. Pharmacoepidemiology, as a field, uses observational methods to evaluate the safety and effectiveness of drug treatments in populations. Pharmacoepidemiology by virtue of its focus, tradition, and research orientation can provide appropriate study designs and analysis methods for precision medicine. The objective of this manuscript is to demonstrate how pharmacoepidemiology can impact and shape precision medicine and serve as a reference for pharmacoepidemiologists interested in contributing to the science of precision medicine. This paper depicts the state of the science with respect to the need for pharmacoepidemiology and pharmacoepidemiological methods, tools and approaches for precision medicine; the need for and how pharmacoepidemiologists use their skills to engage with the precision medicine community; and recommendations for moving the science of precision medicine pharmacoepidemiology forward. We propose a new integrated multidisciplinary approach dedicated to the emerging science of precision medicine pharmacoepidemiology.

Familial genetic and environmental risk profile and high blood pressure event: a prospective cohort of cardio-metabolic and genetic study

BACKGROUND AND AIMS

High blood pressure is the heritable risk factor for cardiovascular diseases. We investigated whether the presence of familial genetic and environmental risk factors are associated with increased risk of high blood pressure.

METHODS

A total of 4,559 individuals from 401 families were included in this study. Familial aggregation analysis was carried out on systolic blood pressure (SBP), diastolic blood pressure (DBP), body mass index (BMI) and waist circumference (WC), and heritability was estimated for SBP and DBP. The association between familial risk factors and blood pressure traits including, incidence of hypertension, SBP and DBP was estimated separately using regression-based two-level Haseman-Elston (HE) method, with individual and familial BMI and WC as environmental exposures and familial genetic profile of known variants as genetic risk factors in 210 index families (≥2 hypertensive cases). Models were adjusted for the two nested sets of covariates.

RESULTS

During a follow-up of 15 years, the SBP, DBP, BMI and WC were highly correlated in inter class of mother-offspring and intraclass of sister-sister with heritability of 30 and 25% for DBP and SBP, respectively. Among index families, those whose members with higher familial BMI or WC had significantly increased risk of hypertension and consistent, strong signals of rs2493134 (AGT) linked with SBP and DBP, rs976683 (NLGN1) linked with SBP and HTN, and epistasis of rs2021783 (TNXB) and known genetic variants linked with all blood pressure traits.

CONCLUSIONS

Findings from this study show that familial genetic and environmental risk profile increase risk for high blood pressure beyond the effect of the individuals' own risk factors.

The relationship between blood pressure and risk of atrial fibrillation: a Mendelian randomization study

AIMS

Observational studies suggest elevated blood pressure (BP) as the leading risk factor for incident atrial fibrillation (AF), but whether this relationship is causal remains unknown. In this study, we used Mendelian randomization (MR) to investigate the potential causal association of BP levels with the risk of developing AF.

METHODS AND RESULTS

Genetic variants associated with the BP traits were retrieved from the International Consortium of Blood Pressure-Genome Wide Association Studies (N = 299 024). From 901 reported variants, 894 were assessed in a dedicated Genome-Wide Association Study of AF genetics, including >1 000 000 subjects of European ancestry. We used two-sample MR analyses to examine the potential causal association of systolic BP (SBP) and diastolic BP (DBP) as well as of pulse pressure (PP) with AF. MR analysis identified a potentially causal association between AF and SBP [odds ratio (OR): 1.018 per 1 mmHg increase, 95% confidence interval (CI): 1.012-1.024, P < 0.001], DBP (OR: 1.026, 95% CI: 1.016-1.035, P < 0.001), and PP (OR: 1.014, 95% CI: 1.001-1.028, P = 0.033). These findings were robust in sensitivity analyses, including the MR-Egger method and the MR pleiotropy residual sum and outlier test (MR-PRESSO). The causal relationship of BP and AF did not change when single-nucleotide polymorphisms associated with possible confounders (i.e. coronary artery disease and obesity) of the causal relationship were excluded.

CONCLUSIONS

The association between increased BP levels and the risk of AF is likely causal and applies for different BP indices. Independently from other risk factors, optimal BP control might represent an important therapeutic target for AF prevention in the general population.

Multi-trait transcriptome-wide association studies with probabilistic Mendelian randomization

A transcriptome-wide association study (TWAS) integrates data from genome-wide association studies and gene expression mapping studies for investigating the gene regulatory mechanisms underlying diseases. Existing TWAS methods are primarily univariate in nature, focusing on analyzing one outcome trait at a time. However, many complex traits are correlated with each other and share a common genetic basis. Consequently, analyzing multiple traits jointly through multivariate analysis can potentially improve the power of TWASs. Here, we develop a method, moPMR-Egger (multiple outcome probabilistic Mendelian randomization with Egger assumption), for analyzing multiple outcome traits in TWAS applications. moPMR-Egger examines one gene at a time, relies on its cis-SNPs that are in potential linkage disequilibrium with each other to serve as instrumental variables, and tests its causal effects on multiple traits jointly. A key feature of moPMR-Egger is its ability to test and control for potential horizontal pleiotropic effects from instruments, thus maximizing power while minimizing false associations for TWASs. In simulations, moPMR-Egger provides calibrated type I error control for both causal effects testing and horizontal pleiotropic effects testing and is more powerful than existing univariate TWAS approaches in detecting causal associations. We apply moPMR-Egger to analyze 11 traits from 5 trait categories in the UK Biobank. In the analysis, moPMR-Egger identified 13.15% more gene associations than univariate approaches across trait categories and revealed distinct regulatory mechanisms underlying systolic and diastolic blood pressures.