Multi-Trait Genome-Wide Association Study of Atherosclerosis Detects Novel Pleiotropic Loci

Although affecting different arterial territories, the related atherosclerotic vascular diseases coronary artery disease (CAD) and peripheral artery disease (PAD) share similar risk factors and have shared pathobiology. To identify novel pleiotropic loci associated with atherosclerosis, we performed a joint analysis of their shared genetic architecture, along with that of common risk factors. Using summary statistics from genome-wide association studies of nine known atherosclerotic (CAD, PAD) and atherosclerosis risk factors (body mass index, smoking initiation, type 2 diabetes, low density lipoprotein, high density lipoprotein, total cholesterol, and triglycerides), we perform 15 separate multi-trait genetic association scans which resulted in 25 novel pleiotropic loci not yet reported as genome-wide significant for their respective traits. Colocalization with single-tissue eQTLs identified candidate causal genes at 14 of the detected signals. Notably, the signal between PAD and LDL-C at the PCSK6 locus affects PCSK6 splicing in human liver tissue and induced pluripotent derived hepatocyte-like cells. These results show that joint analysis of related atherosclerotic disease traits and their risk factors allowed identification of unified biology that may offer the opportunity for therapeutic manipulation. The signal at PCSK6 represent possible shared causal biology where existing inhibitors may be able to be leveraged for novel therapies.

Introduction

Atherosclerotic vascular disease is a leading cause of death worldwide (Lozano et al., 2012; Kobiyama and Ley, 2018; Virani et al., 2020) and can affect multiple arterial territories. Although clear differences in disease pathobiology exist (Lin et al., 2013), epidemiological analyses have shown both coronary artery disease (CAD) and peripheral artery disease (PAD) share similar risk factors and frequently co-occur (Ozkaramanli Gur et al., 2018; Klarin et al., 2019; Sundaram et al., 2020). These risk factors include dyslipidemia, obesity, hypertension, diabetes, and tobacco use (Criqui and Aboyans, 2015). PAD patients with concomitant CAD are known to experience more extensive and aggressive disease (Hussein et al., 2011).

The genetics of CAD have been well characterized and a number of genome-wide association studies (GWAS) have identified over 200 genetic risk loci with robust connections to CAD (Khera and Kathiresan, 2017; Van Der Harst and Verweij, 2018; Koyama et al., 2020). For most loci, however, underlying mechanisms by which these loci influence CAD risk remains unclear. Although PAD has been less intensively studied, recent work has identified 21 total risk loci associated with PAD risk (Matsukura et al., 2015; Klarin et al., 2019). Genetic correlation studies have demonstrated a high degree of shared genetic architecture between CAD and PAD (LD-score regression-based genetic correlation r_g = 0.67) (Purcell et al., 2007). This genetic correlation, based on shared pathobiology, can be leveraged to identify novel pleiotropic genetic architecture common to both disease traits (Zhao et al., 2017; Baselmans et al., 2019).

The development of statistical approaches for multi-trait GWAS meta-analysis has facilitated joint analyses of traits with substantial evidence for a common pathophysiological basis to elucidate shared genetic etiology (Klarin et al., 2019). Furthermore, correlated causal risk factors can also be included in these multi-trait GWAS analyses to provide insight on their shared genetic pathways (Holmes et al., 2015; Zhao et al., 2017; Riaz et al., 2018; Siewert and Voight, 2018; Larsson et al., 2020). Our previous work has analyzed CAD pairwise with secondary traits to understand shared genetic etiology to successfully identify new risk loci (Zhao et al., 2017; Siewert and Voight, 2018). Yet, there have been no studies which evaluate atherosclerosis endpoints jointly with multiple cardiometabolic causal risk factors for discovery and quantitative interpretation.

In this study, we performed a series of N-weighted multivariate genome-wide-association meta-analyses (N-GWAMA) (Baselmans et al., 2019) using different combinations of nine atherosclerotic or atherosclerosis risk factor traits, and identified 31 unique pleiotropic loci not previously associated with any analyzed trait combination. We subsequently used single-tissue expression quantitative trait loci (eQTL) colocalization analysis at these loci to identify candidate causal genes and their tissue site of action. Some of these causal gene candidates have potential opportunities for drug target repurposing to treat atherosclerotic vascular disease, including PCSK6. Ultimately, this study provides a better understanding of biological pathways underlying atherosclerosis to inform future therapeutic development.

Methods

This study was approved by the U.S. Department of Veterans Affairs Central Institutional Review Board. All participants gave written informed consent for study participation.

Genetic Association Data

We collected the summary statistics from the largest published GWAS to maximize our power for novel discovery. PAD summary statistics were the European ancestry subjects from the recent VA Million Veteran Program analysis which consisted of 24,009 PAD cases and 150,983 PAD controls (Klarin et al., 2019). These data can be accessed from dbGAP (phs001672). CAD data were taken from CARDIoGRAMplusC4D combined with the United Kingdom BioBank (UKBB) (Van Der Harst and Verweij, 2018) and consisted of 122,733 CAD cases and 424,528 CAD controls. Data for body mass index (BMI) (meta-analysis of GIANT and UKBB; 806,834 individuals; (Yengo et al., 2018)), type 2 diabetes (T2D) (meta-analysis of consortia; 228,499 cases and 1,178,783 controls; (Vujkovic et al., 2020)), smoking initiation (smoking) (UKBB; 462,690 individuals; (Wootton et al., 2020)), and 4 lipid traits (meta-analysis of MVP and GLGC data; 723,000 participants; (Klarin et al., 2018)). Access urls for all data obtained from the public domain are provided in Supplementary Table S1.

N-GWAMA Multi-Trait GWAS

Using the summary statistics from publicly available single-trait GWAS (Supplementary Table S1), we performed 15 N-GWAMA (Baselmans et al., 2019) multi-trait GWAS centered around PAD, CAD, and the following atherosclerotic risk factor traits: BMI, smoking, T2D, LDL-C, HDL-C, TC, and TG. Full details are provided in Supplementary Methods. Briefly, we first performed a bivariate GWAS for PAD and CAD followed by a series of trivariate GWAS combining PAD, CAD, and one of seven correlated traits that represented traditional atherosclerotic risk factors. We also performed a series of bivariate GWAS between PAD and these seven traits individually, given that a series of bivariate GWAS between CAD and most of these seven traits has already been performed (Siewert and Voight, 2018).

Each N-GWAMA multi-trait GWAS resulted in a set of independent loci represented by a sentinel SNP. We defined an independent locus as the genomic region that includes all variants within 1 megabase (Mb) of the sentinel SNP and any other variants that were in linkage diseqilibrium (LD) of r² > 0.2 with the sentinel SNP using the 1,000 Genomes European ancestry cohort (1 kG EUR) (Purcell et al., 2007). We then applied a series of filters to remove loci that were not plausibly pleiotropic, or did not represent novel associations. To ensure there was evidence a locus was pleiotropic and that a single trait was not driving the association, we required that the sentinel SNP was at least nominally associated (p < 5 × 10⁻³) with all the individual traits involved in the multi-trait analysis. We also required that none of the variants at an independent locus were previously associated with any of the traits used in the N-GWAMA multi-trait GWAS by applying two filters. First, we required each sentinel SNP was not genome-wide significant for any of the individual traits (p > 5 × 10⁻⁸). It was also necessary that none of the SNPs at the independent locus were previously reported to be genome-wide significant for any of the individual traits involved in the multi-trait analysis in the GWAS Catalog (Buniello et al., 2019) (Figure 1). Finally, we excluded loci in the HLA region from these experiments due to the difficulty of interpreting the independent signals of these loci. Code for the pipeline is available at: https://github.com/Bellomot/Athero_NGWAMA_Multitrait_GWAS.

FIGURE 1

FIGURE 1. Flowchart of multi-trait analysis and candidate gene results. 9 traits were analyzed in 15 different bivariate and trivariate scans that best represented atherosclerosis. The summary statistics from all scans were filtered by single trait p-values and loci within 500 kb or in LD (EUR r² > 0.21 kG EUR) with the known trait being tested according to the GWAS Catalog, resulting in 150 unique loci. Trait-to-trait colocalization with a threshold of a conditional posterior probability of colocalization >0.8 was performed to ensure evidence of a shared causal SNP between each trait. The resulting 31 unique loci were run through single tissue eQTL analysis using GTEx v8 to identify candidate causal genes and tissues for each locus. 34 unique genes were identified among 14 loci.

Given that we performed 15 multi-trait GWAS for combinations of related traits, we next implemented a multiple testing correction procedure to assess significance. Due to the high correlation between each of the multi-trait GWAS we performed, a Bonferroni correction (p < 3.3 × 10⁻⁹) for each trait combination test is conservative. Thus, we constructed a null distribution Z-score sampling strategy to estimate an α = 0.05 p-value threshold given the set of N-GWAMA multi-trait GWAS that we performed. Under the assumption that the correlation of the Z-scores across the entire genome that resulted from the N-GWAMA multi-trait GWAS are a reasonable estimate of the correlation in the multivariate null distribution of Z-scores, we can use the correlation matrix of the Z-scores from the 15 multi-trait GWAS to model the 15-dimentional multivariate standard normal distribution that is the theoretical null distribution of these results.

To get this estimate of the appropriate α = 0.05 p-value threshold, we first drew 10,000 sets of 1 million samples from a 15-dimentional multivariate normal distribution centered at the origin, and used the correlation between the Z-scores of all the SNPs that were tested across all 15 N-GWAMA multi-trait GWAS as the correlation matrix (Supplementary Table S2). We kept the most extreme Z-score from each of the 10,000 sets and then identified the 95th percentile of the most extreme Z-scores as our α = 0.05 threshold. We defined experiment-wide significance as the 95th percentile Z-score of 5.87, which corresponds to a p-value of 4.3 × 10⁻⁹.

Trait-Trait Colocalization

For each multi-trait GWAS, we assessed the evidence of a shared causal variant at each significant locus by performing colocalization analysis between the trait signals using COLOC for bivariate GWAS and MOLOC for trivariate GWAS (Giambartolomei et al., 2014; Giambartolomei et al., 2018) (Supplementary Table S3). For this analysis, we applied a 500 kilobases (Kb) window (+/− 250 Kb) around the sentinel SNP. A conditional probability of colocalization is defined as the posterior probability of colocalization conditioned on the presence of a signal for each trait. For example, when using COLOC, this would be represented as the posterior probability of hypothesis four (PP4) over the sum of the posterior probabilities of hypothesis three (PP3) and PP4 (PP4/(PP3 + PP4)) (Giambartolomei et al., 2014). A probability of ≥0.8 was considered of interest. Loci that had a conditional probability of colocalization >0.5 and <0.8 were visually inspected using LocusZoom plots (Supplementary Table S4). If the LD structure suggested additional associations unlinked to the sentinel SNP, approximate conditional analysis was performed (see details on Approximate Conditional Analysis below).

Single-Tissue Gene Expression Colocalization

We performed single-tissue colocalization analysis to prioritize candidate causal genes implicated in each N-GWAMA multi-trait GWAS using RNA-seq data obtained from the Genotype-Tissue Expression (GTEx) project (Giambartolomei et al., 2014). We identified the list of genes and tissues for which each N-GWAMA sentinel SNP was a significant single-tissue eQTL in any GTEx v8 tissue as reported in the “.signif_variant_gene_pairs.txt.gz” files available from the GTEx Portal (Supplementary Table S5). We then performed colocalization between either CAD or PAD, as determined by which trait had the most significant sentinel SNP at each locus, and each single-tissue eQTL signal (Giambartolomei et al., 2014). We selected the window of colocalization to be 500 Kb spanning the sentinel SNP. Similar to trait-trait colocalization analysis, our threshold to classify the traits as colocalized was a conditional probability of colocalization (PP4/(PP3 + PP4)) ≥ 0.8. We visually inspected LocusZoom plots for loci where colocalization analysis resulted in a conditional probability of colocalization <0.8 but >0.5 and performed approximate conditional analysis when the LD structure suggested possible allelic heterogeneity (see details on Approximate Conditional Analysis below).

Approximate Conditional Analysis

For each locus that showed evidence of multiple independent signals, we performed approximate conditional analysis on variants that appeared to be associated with the trait of interest independently of the sentinel SNP (Supplementary Table S6). This analysis was necessary given that the presence of multiple associated variants in a region violates the assumptions of COLOC (Giambartolomei et al., 2014). Potential nearby association signals were identified using LocusZoom plots and the LDlink LDassoc tool between sentinel and putative distinctive signal variant quantified using LDlink (Pruim et al., 2010; Machiela and Chanock, 2018). We performed approximate conditional analysis using GCTA-COJO with 1,000 Genome Project data (European samples, n = 503) as a reference panel (Yang et al., 2011; 1000 Genomes Project Consortium et al., 2012). We conditioned the sentinel SNP on the most associated variant for each potential confounding signal identified at the locus. We then repeated the colocalization experiment on the locus using the conditional variant p-values. A full list of traits, the sentinel SNPs, and the conditioned variants for each conditional analysis are provided in the supplement (Supplementary Table S6).

Splicing Quantitative Trait Locus Colocalization

We performed a colocalization analysis between the PAD signal at the PCSK6 locus and the GTEx v8 liver tissue splicing quantitative trait locus (sQTL) signal with the intron ID: 101365044:101366196:clu_14775. We also identified this intron signal in the Phenotyping Lipid traits in iPS derived hepatocytes Study (PhLiPS) hepatocyte-like cell (HLC) sQTL data by lifting over the start and stop of this intron to hg19 (101905249:101906401) (Gawronski et al., 2019). We then performed colocalization analysis between the HLC sQTL signal and the PAD signal as well as the HLC sQTL signal and the GTEx v8 liver tissue sQTL signal (Figure 2).

FIGURE 2

FIGURE 2. PCSK6 locus with a sentinel SNP of rs1531817. Pleiotropic signal between PAD and LDL-C with an sQTL for PCSK6 in liver tissue. This locus also colocalized with hepatocyte-like cells (HLCs) in vitro.

Results

Multi-Trait GWAS Analysis Results

We first calculated the genetic correlation between PAD and CAD and the seven atherosclerosis risk factors (BMI, smoking, T2D, LDL-C, HDL-C, TC, and TG) using the summary statistics files for these GWAS using LD score regression (Supplementary Figure S1; Supplementary Table S7). We then performed 15 N-GWAMA scans centered around PAD and CAD to detect novel loci not previously reported as genome-wide significant for any of their respective traits (Supplementary Table S8). A total of 150 sentinel SNPs were multivariate genome-wide significant with all single trait p-values between <5 × 10⁻³ and >5 × 10⁻⁸ (Supplementary Table S3). Of these sentinel SNPs, 31 were nominal genome-wide significant (multi-trait p < 5 × 10⁻⁸) and met our trait-to-trait colocalization criteria, and thirteen were experiment-wide significant (multi-trait p < 4.3 × 10⁻⁹) and met our trait-to-trait colocalization criteria. If we had used the Bonferroni correction threshold, only one locus, NUP85 (p = 3.55 × 10⁻⁹), would change from experiment-wide significant to genome-wide significant. The nominal genome-wide significant sentinel SNPs represent 25 independent loci and 11 experiment-wide significant independent loci (>1 Mb from any of our other reported loci, Table 1, Supplementary Figures S2–S34). Fourteen of the nominal genome-wide significant sentinel SNPs and eight experiment-wide significant sentinel SNPs colocalized with one or more single-tissue eQTLs (Supplementary Table S5). Finally, five loci colocalized with eQTLs for genes that have been implicated in atherosclerosis by previous studies.

TABLE 1

TABLE 1. Atherosclerosis trait N-GWAMA analysis and results. Trait 3 p value will have a value of NA if there were only 2 traits analyzed. Conditional posterior probability represents the probability of the trait-to-trait colocalization analysis (e.g., PP4/(PP3 + PP4)).

We noted that two of our signals that exceeded multi-trait experiment-wide significance mapped to a previously established locus for several cardiometabolic traits. That signal was tagged by rs2076211, and associated with PAD, CAD, and BMI (bivariate p = 7.6 × 10⁻¹⁰) or PAD, CAD, and LDL (bivariate p = 4 × 10⁻⁹). This variant mapped to the nearby gene PNPLA3, a well-established locus associated for non-alcoholic fatty liver disease (Speliotes et al., 2011), multiple liver enzymes measures (Yuan et al., 2008), hemotological traits (Kamatani et al., 2010), sex-hormone binding globulin levels (Ruth et al., 2020), and T2D (Mahajan et al., 2018). Our sentinel SNP that tagged both multi-trait signals was in strong LD with rs738409 (r² = 0.73 1 kG EUR), the previously reported sentinel variant associated with these additional traits. A previous multi-trait scan for T2D and LDL reported strong association for both traits at this locus (Klimentidis et al., 2020). However, association with this locus and atherosclerotic disease (PAD, CAD) or to BMI to our knowledge has not been previously reported, but are compelling given the extensive pleiotropy for atherosclerotic causal risk factors here.

PCSK6 Locus

We detected a signal that exceeded multi-trait experiment-wide significance with PAD and LDL-C (bivariate p = 3.2 × 10⁻¹⁰) at the PCSK6 locus. A rare coding variant in this region has been reported to associate with LDL-C (Klimentidis et al., 2020; Sinnott-Armstrong et al., 2021), however, the coding variant (NP_002561.1:p.Thr964Met, rs34631529) and our sentinel SNP (rs1531817) are not in linkage disequilibrium (r² = 0.0086 1 kG EUR) based on data from the 1,000 Genomes Project (Purcell et al., 2007), indicating that we detected a different signal at this locus. To further differentiate whether our signal was novel, we performed an additional conditional analysis on the coding variant rs34631529 in PAD data without any notable changes in the PCSK6 locus signal (Supplementary Figure S36). We further note that previous GWAS have found that variants at this locus are associated with inflammatory markers (Iyengar et al., 2015; Hackinger et al., 2018; Nath et al., 2019; Folkersen et al., 2020; Richardson et al., 2020; Ruotsalainen et al., 2021).

To better understand how genetic variation at the PCSK6 locus influences circulating lipid levels, we investigated the association of the bivariate sentinel SNP at this locus in the publicly available GWAS of NMR lipid subfractions: extra-small subfrations (XS), extra-large subfractions (XL), HDL, intermediate density lipoprotein (IDL), LDL, and very-low density lipoprotein (VLDL) (Kettunen et al., 2016). We found our sentinel SNP (rs1531817) had a nominal association with medium VLDL particles (β = 0.03, SE = 0.01, p = 9 × 10⁻³), total lipids in medium VLDL (β = 0.03, SE = 0.01, p = 0.02), TG in medium VLDL (β = 0.02, SE = 0.01, p = 0.03), serum TG (β = 0.02, SE = 0.01, p = 0.03), and TG in large VLDL (β = 0.02, SE = 0.01, p = 0.03).

Our sentinel SNP was a sQTL for PCSK6 in GTEx v8 liver tissue (Figure 2). To identify a potential experimental model of this splicing change, we searched for this sQTL in PhLiPS HLC summary data (PhLiPS Study) (Pashos et al., 2017). The signal at PCSK6 colocalized with an sQTL in these data as well (Figure 2), which suggests that derived hepatocyte-like cells could be a good model for further studied of the effect of this locus on LDL-C and PAD risk.

SORCS3 Locus

We detected a nominal genome-wide significant signal with PAD and TC (bivariate p = 2.4 × 10⁻⁸) at the SORCS3 locus rs11599236 (Table 1). This locus was previously observed to be genome-wide significant in GWAS studies of mood disorders (Howard et al., 2018; Ward et al., 2020) (Supplementary Table S4). This signal colocalized with SORCS3 mRNA expression in pituitary tissue (Supplementary Figure S34). The opposite direction of effect was noted for both traits and the gene-tissue pair: decreased SORCS3 associated with increased PAD and TC (β = 0.31, SE = 0.05, p = 5.5 × 10⁻⁸).

Other Candidate Genes Indicated With Known Atherosclerotic Biology

We detected a trivariate GWAS signal with PAD, CAD, and smoking (trivariate p = 1.9 × 10⁻⁹) at the HMBS locus rs1006195 (Table 1). This variant was genome-wide significant in previous GWAS studies for several cardiometabolic traits including Apolipoprotein A1 levels, waist-hip ratio, BMI, fat mass percentage, HDL-C, and T2D (Turcot et al., 2018; Van Der Harst and Verweij, 2018; Pulit et al., 2019; Richardson et al., 2020; Zhu et al., 2020). This pleiotropic signal colocalized with HMBS and VPS11 mRNA expression in several tissues (Supplementary Figure S14; Supplementary Table S5). HMBS (β = 0.28, SE = 0.03, p = 4.1 × 10⁻²³) demonstrated the same direction of effect with PAD, CAD, and smoking, indicating that increased HMBS expression is associated with increased PAD, CAD, and smoking risk. VPS11 (β = −0.25, SE = 0.03, p = 6.0 × 10⁻¹⁵) demonstrated the opposite direction of effect in all tissue except skeletal muscle and the left ventricle of the heart, meaning increased VPS11 expression in skeletal muscle and left ventricle is associated with increased PAD, CAD, and smoking risk.

We also detected a nominal genome-wide significant signal with PAD and CAD (bivariate p = 7.8 × 10⁻⁹) at the VDAC2 locus rs7088974 (Table 1). Variants at this locus have been found to be associated with BMI in previous GWAS(Pulit et al., 2019) (Supplementary Table S4). Although variants near this locus have also been associated with smoking behavior, our results suggest that the locus we detected is independent of smoking behavior (Liu et al., 2019) (Supplementary Figure S35). This signal colocalized with VDAC2 mRNA expression in multiple vascular tissues relevant to atherosclerosis, including aorta and tibial artery (Supplementary Table S5). The direction of effect in all tissue was opposite to the direction of effect for PAD and CAD: the allele associated with increased VDAC2 expression (EA = C, EAF = 0.57) is associated with decreased PAD and CAD (β = −0.14, SE = 0.02, p = 7.7 × 10⁻⁹).

Finally, we detected a trivariate signal with PAD, CAD, and HDL-C (trivariate p = 3.4 × 10⁻¹¹) at the LRCH1 locus rs9526214 (Table 1). This locus had evidence of allelic heterogeneity when we reviewed the regional association plots, which led us to perform approximate conditional analyses on the pleiotropic signal sentinel SNP rs9316223 and the resulting conditional probability of colocalization met our criteria (Supplementary Table S5). This locus has been found to be genome-wide significant in previous GWAS studies for platelets, systolic blood pressure, and stroke (Evangelou et al., 2018; Malik et al., 2018) (Supplementary Table S4). This signal colocalized with LRCH1 mRNA expression in tibial artery, whole blood, and other tissues (Supplementary Figure S12; Supplementary Table S5). The opposite direction of effect was noted for all three traits and the gene tissue pair: the allele associated with increased LRCH1 was also associated with decreased PAD, CAD, and HDL-C (β = −0.09, SE = 0.02, p = 1.3 × 10⁻⁸).

Discussion

To advance our understanding of the genetic etiology of atherosclerosis, different combinations of nine known atherosclerotic or atherosclerosis risk factor traits were used to perform 15 N-GWAMA scans which resulted in 25 unique novel pleiotropic loci (Figure 1). Colocalization with single-tissue eQTLs identified 34 candidate causal genes across 14 of the detected signals. Five of these loci had candidate causal genes previously associated with atherosclerosis through other studies. While candidate causal genes remain elusive for the remaining loci, the patterns of association represent physiology that appears compelling. For example, our top association at the SATB1 locus was modestly associated with CAD, PAD, and T2D in the same direction of effect, implying that a perturbation informed by the human genetics data might be expected to be ameliorative for all three traits. Functional work to elucidate causal variant, genes, and mechanism at these loci thus may provide new insights into the etiological pathways for this collection of disease endpoints.

PCSK6 Activity Effects Lipid Levels, Plaque Formation and Stability

We identified a signal at the PCSK6 locus that has a bivariate association with PAD and LDL-C and provide strong evidence in support of PCSK6 as the causal gene at the PCSK6 locus. PCSK6 is a calcium-dependent serine endoprotease that cleaves proteins to active and inactive forms depending on the target protein (Kiefer et al., 1991).

There is convincing experimental evidence to suggest that PCSK6 directly influences plaque development and plaque stability. Smooth muscle cell migration in injured arteries is facilitated by cytokine induced PCSK6 expression that activates matrix metalloproteinases (MMP14/MMP2) (Perisic et al., 2013). This smooth muscle cell mechanism may explain the association of PCSK6 with carotid intima-media thickness in a candidate gene study (Rykaczewska et al., 2020). The sentinel SNP associated with maximum progression of carotid intima-media thickness was the same as the sentinel SNP identified in our bivariate scan between LDL-C and PAD. There is also evidence that PCSK6 activates MMP9, which enhances degredation of the extra cellular matrix and thus promotes plaque instability (Li et al., 2020; Testa et al., 2021).

Our data suggest PCSK6 also influences lipid metabolism, a known upstream cause of atherosclerotic progression. This is in agreement with the known role of PCSK6 in lipid metabolism, where it cleaves and inactivates endothelial lipase (EL) and lipoprotein lipase (LPL) (Jin et al., 2005), which can lead to hyperlipidemia (Choi and Korstanje, 2013). This finding is further supported by the recently reported rare coding variant in PCSK6 associated with decreased LDL-C (Klimentidis et al., 2020; Sinnott-Armstrong et al., 2021).

It remains to be determined how changes in PCSK6 activity result in altered lipid metabolism and risk of PAD. It is possible that changes in PCSK6 activity alter lipid metabolism and therefore have downstream effects on PAD or altered PCSK6 activity may effect LDL-C and PAD separately through independent mechanisms. In all likelihood, it is a combination of LDL-C dependent and independent mechanisms that link PCSK6 to PAD.

From a drug targeting standpoint, we would expect that inhibiting PCSK6 reduces LDL-C levels and PAD risk. There are several non-FDA approved, non-specific PCSK6 inhibitors that influence lipid metabolism: alpha1-antitrypsin Portland (alpha1-PDX) (Gauster et al., 2005), profurin (Jin et al., 2007), Pf-pep (Byun et al., 2010), and dicoumarols (DC), specifically DC2 (Komiyama et al., 2009).

SORCS3 Influences Energy Metabolism

We identified a novel pleiotropic association at the SORCS3 locus between PAD and TC. Colocalization experiments support SORCS3 as the causal gene at this locus. SORCS3 is a type I transmembrane protein that is a member of the Vps10p receptor family (Hermey, 2009). The post-synaptic sorting receptor SORCS3 is highly expressed in the hippocampus and binds NGF and PDGF-BB to modulate several signal transduction pathways involved in neuronal activity (Hermey, 2009; Christiansen et al., 2017). A recent knockout mouse study has shown that even loss of one Sorcs3 allele decreases lipid metabolism as a source of energy and increases adiposity (Subkhangulova et al., 2018). The proposed mechanism for this occurrence is that loss of SORCS3 expressed in agouti-related peptide (AGRP) neurons leads to unchecked production of circulating AGRP (Henry et al., 2015). ARGP is well described to cause enhanced food intake, decreased locomotor activity, decreased use of lipids as metabolic fuel, and overall increased adiposity (Maier et al., 2018; Beutler et al., 2020). One limitation of GTEx data is a very small sample size of hypothalamus data, our data show that reduced expression of SORCS3 in the pituitary gland is related to high levels of TC and an increased overall risk for PAD, potentially through increased circulating AGRP diverting energy metabolism away from lipid sources.

Limitations

We acknowledge there are several limitations to this study. First, there is sample overlap between several of our single trait summary statistics files. The N-GWAMA method attempts to account for this, but if the correction for the overlap was insufficient this could inflate our false discovery rate. Second, in using the null Z-score resampling strategy to estimate the appropriate multiple testing correction, we assume that the correlation of the Z-scores across the entire genome that resulted from the N-GWAMA multi-trait GWAS are a reasonable approximation of the correlation between the multivariate null distribution of Z-scores. This assumption could be broken due to highly associated SNPs leading to an altered correlation matrix compared to the true null distribution and the Z-scores having an inflated variance from confounders that may not have been corrected for by LDSC. Finally, we selected nine atherosclerotic and cardiometabolic traits based on conventional relationships with atherosclerosis; however, there are likely multifactorial and multidirectional relationships within this group of traits. It is possible that some of the novel loci represent the interaction between traits instead of the intended representation of atherosclerosis as we have interpreted it. It is also possible that including other sets of cardiometabolic risk factors may identify additional novel loci.

Conclusion

We have shown that publicly available GWAS data can be leveraged to perform multi-trait scans with N-GWAMA methods to identify novel loci that unify atherosclerosis. In this study, 25 nominal genome-wide significant loci were associated jointly with PAD and other atherosclerotic traits. These loci may represent novel genetic etiologies of atherosclerosis. A total of 34 candidate causal genes were identified across 14 novel pleiotropic loci and among those, PCSK6 represents possible causal biology with known inhibitors that has large potential to be a therapeutic target for atherosclerosis. These results highlight the biological underpinnings of atherosclerosis and the potential to develop non-invasive medical treatments for atherosclerosis.

Data Availability Statement

The data analyzed in this study is subject to the following licenses/restrictions: The GWAS and GTEx eQTL data are publicly available. GWAS Summary Statistics also in Supplementary Table S1: PAD: dbGAP accession code phs001672.v2.p1; CAD: Mendeley doi: 10.17632/gbbsrpx6bs.1; Lipids: MVP, dbGAP accession codes phs001672.v4.p1, pha004828.1, pha004831.1, pha004837.1, pha004834.1; GLGC, http://csg.sph.umich.edu/willer/public/lipids/; T2D: dbGAP under accession number phs001672.v3.p1; BMI: https://portals.broadinstitute.org/collaboration/giant/index.php/GIANT_consortium_data_files#20; Smoking: https://data.bris.ac.uk/data/dataset/10i96zb8gm0j81yz0q6ztei23d; GTEx v8 data is available at: https://gtexportal.org/home/datasets Additional data that support the findings of this study are available on request from the coauthor (SD); these data are not publicly available due to U.S. Government and Department of Veteran’s Affairs restrictions relating to participant privacy and consent. Requests to access these datasets should be directed to SD c2NvdHQuZGFtcmF1ZXJAcGVubm1lZGljaW5lLnVwZW5uLmVkdQ==. Code to perform analyses in this article are available from the authors upon request (BV, SD), or from github at https://github.com/Bellomot/Athero_NGWAMA_Multitrait_GWAS.

Author Contributions

TB and WB ran experiments, interpreted results, and wrote the manuscript BC, KG, DZ, JP ran experiments for follow up analyses on loci BC, ML, NT wrote code for performing analyses DK, JL, TA, JG, PW, KC, MV, CO, K-MC, PT collected data to make these experiments possible DR and MR helped interpret results, provided data for follow up experiments, and helped write the manuscript SD and BV oversaw the project, helped interpret results, and helped write the manuscript.

Funding

This research is based on data from the Million Veteran Program, Office of Research and Development, Veterans Health Administration, and was supported by awardno. MVP000. This publication does not represent the views of the Department of Veteran Affairs or the United States Government. This research was also supported by funding from: the Department of Veterans Affairs awards nos. I01-BX03340 (KC and PW), I01-BX003362 (PT and KC), and IK2-CX001780 (SD), the National Institutes of Health (DK101478 and DK126194 to BV), the American Heart Association (20PRE35120109 to WB).

Conflict of Interest

SD receives research support to his institution from RenalytixAI and personal consulting frees from Calico Labs, both outside the current work. CO is employed by Novartis Institutes of Biomedical Research.

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

Publisher’s Note

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

Supplementary Material

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

References

1000 Genomes Project Consortium et al., 2012 1000 Genomes Project Consortium Abecasis, G. R., Auton, A., Brooks, L. D., Depristo, M. A., Durbin, R. M., et al. (2012). An Integrated Map of Genetic Variation from 1,092 Human Genomes. Nature 491, 56–65. doi:10.1038/nature11632

PubMed Abstract | CrossRef Full Text | Google Scholar

Baselmans, B. M. L., Jansen, R., Ip, H. F., van Dongen, J., Abdellaoui, A., van de Weijer, M. P., et al. (2019). Multivariate Genome-wide Analyses of the Well-Being Spectrum. Nat. Genet. 51, 445–451. doi:10.1038/s41588-018-0320-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Beutler, L. R., Corpuz, T. V., Ahn, J. S., Kosar, S., Song, W., Chen, Y., et al. (2020). Obesity Causes Selective and Long-Lasting Desensitization of AgRP Neurons to Dietary Fat. Elife 9, e55909. doi:10.7554/eLife.55909

PubMed Abstract | CrossRef Full Text | Google Scholar

Buniello, A., Macarthur, J. A. L., Cerezo, M., Harris, L. W., Hayhurst, J., Malangone, C., et al. (2019). The NHGRI-EBI GWAS Catalog of Published Genome-wide Association Studies, Targeted Arrays and Summary Statistics 2019. Nucleic Acids Res. 47, D1005–D1012. doi:10.1093/nar/gky1120

PubMed Abstract | CrossRef Full Text | Google Scholar

Byun, S., Tortorella, M. D., Malfait, A.-M., Fok, K., Frank, E. H., and Grodzinsky, A. J. (2010). Transport and Equilibrium Uptake of a Peptide Inhibitor of PACE4 into Articular Cartilage Is Dominated by Electrostatic Interactions. Arch. Biochem. Biophys. 499, 32–39. doi:10.1016/j.abb.2010.04.019

PubMed Abstract | CrossRef Full Text | Google Scholar

Choi, S., and Korstanje, R. (2013). Proprotein Convertases in High-Density Lipoprotein Metabolism. Biomark Res. 1, 27. doi:10.1186/2050-7771-1-27

PubMed Abstract | CrossRef Full Text | Google Scholar

Christiansen, G. B., Andersen, K. H., Riis, S., Nykjaer, A., Bolcho, U., Jensen, M. S., et al. (2017). The Sorting Receptor SorCS3 Is a Stronger Regulator of Glutamate Receptor Functions Compared to GABAergic Mechanisms in the hippocampus. Hippocampus 27, 235–248. doi:10.1002/hipo.22689

PubMed Abstract | CrossRef Full Text | Google Scholar

Criqui, M. H., and Aboyans, V. (2015). Epidemiology of Peripheral Artery Disease. Circ. Res. 116, 1509–1526. doi:10.1161/circresaha.116.303849

PubMed Abstract | CrossRef Full Text | Google Scholar

Evangelou, E., Warren, H. R., Mosen-Ansorena, D., Mifsud, B., Pazoki, R., Gao, H., et al. (2018). Genetic Analysis of over 1 Million People Identifies 535 New Loci Associated with Blood Pressure Traits. Nat. Genet. 50, 1412–1425. doi:10.1038/s41588-018-0205-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Folkersen, L., Gustafsson, S., Wang, Q., Hansen, D. H., Hedman, Å. K., Schork, A., et al. (2020). Genomic and Drug Target Evaluation of 90 Cardiovascular Proteins in 30,931 Individuals. Nat. Metab. 2, 1135–1148. doi:10.1038/s42255-020-00287-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Gauster, M., Hrzenjak, A., Schick, K., and Frank, S. (2005). Endothelial Lipase Is Inactivated upon Cleavage by the Members of the Proprotein Convertase Family. J. Lipid Res. 46, 977–987. doi:10.1194/jlr.m400500-jlr200

CrossRef Full Text | Google Scholar

Gawronski, K. A. B., Bone, W. B., Park, Y., Pashos, E., Wang, X., Yang, W., et al. (2019). Evaluating the Contribution of Cell-type Specific Alternative Splicing to Variation in Lipid Levels. biorxiv. doi:10.1101/659326

CrossRef Full Text | Google Scholar

Giambartolomei, C., Vukcevic, D., Schadt, E. E., Franke, L., Hingorani, A. D., Wallace, C., et al. (2014). Bayesian Test for Colocalisation between Pairs of Genetic Association Studies Using Summary Statistics. Plos Genet. 10, e1004383. doi:10.1371/journal.pgen.1004383

PubMed Abstract | CrossRef Full Text | Google Scholar

Giambartolomei, C., Zhenli Liu, J., Zhang, W., Hauberg, M., Shi, H., Boocock, J., et al. (2018). A Bayesian Framework for Multiple Trait Colocalization from Summary Association Statistics. Bioinformatics 34, 2538–2545. doi:10.1093/bioinformatics/bty147

PubMed Abstract | CrossRef Full Text | Google Scholar

Hackinger, S., Prins, B., Mamakou, V., Zengini, E., Marouli, E., Brčić, L., et al. (2018). Evidence for Genetic Contribution to the Increased Risk of Type 2 Diabetes in Schizophrenia. Transl Psychiatry 8, 252. doi:10.1038/s41398-018-0304-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Henry, F. E., Sugino, K., Tozer, A., Branco, T., and Sternson, S. M. (2015). Cell Type-specific Transcriptomics of Hypothalamic Energy-Sensing Neuron Responses to Weight-Loss. Elife 4, e09800. doi:10.7554/eLife.09800

PubMed Abstract | CrossRef Full Text | Google Scholar

Hermey, G. (2009). The Vps10p-Domain Receptor Family. Cell. Mol. Life Sci. 66, 2677–2689. doi:10.1007/s00018-009-0043-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Holmes, M. V., Asselbergs, F. W., Palmer, T. M., Drenos, F., Lanktree, M. B., Nelson, C. P., et al. (2015). Mendelian Randomization of Blood Lipids for Coronary Heart Disease. Eur. Heart J. 36, 539–550. doi:10.1093/eurheartj/eht571

PubMed Abstract | CrossRef Full Text | Google Scholar

Howard, D. M., Adams, M. J., Shirali, M., Clarke, T.-K., Marioni, R. E., Davies, G., et al. (2018). Genome-wide Association Study of Depression Phenotypes in UK Biobank Identifies Variants in Excitatory Synaptic Pathways. Nat. Commun. 9, 1470. doi:10.1038/s41467-018-03819-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Hussein, A. A., Uno, K., Wolski, K., Kapadia, S., Schoenhagen, P., Tuzcu, E. M., et al. (2011). Peripheral Arterial Disease and Progression of Coronary Atherosclerosis. J. Am. Coll. Cardiol. 57, 1220–1225. doi:10.1016/j.jacc.2010.10.034

CrossRef Full Text | Google Scholar

Iyengar, S. K., Sedor, J. R., Freedman, B. I., Kao, W. H., Kretzler, M., Keller, B. J., et al. (2015). Genome-Wide Association and Trans-ethnic Meta-Analysis for Advanced Diabetic Kidney Disease: Family Investigation of Nephropathy and Diabetes (FIND). Plos Genet. 11, e1005352. doi:10.1371/journal.pgen.1005352

PubMed Abstract | CrossRef Full Text | Google Scholar

Jin, W., Fuki, I. V., Seidah, N. G., Benjannet, S., Glick, J. M., and Rader, D. J. (2005). Proprotein Covertases Are Responsible for Proteolysis and Inactivation of Endothelial Lipase. J. Biol. Chem. 280, 36551–36559. doi:10.1074/jbc.m502264200

CrossRef Full Text | Google Scholar

Jin, W., Wang, X., Millar, J. S., Quertermous, T., Rothblat, G. H., Glick, J. M., et al. (2007). Hepatic Proprotein Convertases Modulate HDL Metabolism. Cel. Metab. 6, 129–136. doi:10.1016/j.cmet.2007.07.009

CrossRef Full Text | Google Scholar

Kamatani, Y., Matsuda, K., Okada, Y., Kubo, M., Hosono, N., Daigo, Y., et al. (2010). Genome-wide Association Study of Hematological and Biochemical Traits in a Japanese Population. Nat. Genet. 42, 210–215. doi:10.1038/ng.531

PubMed Abstract | CrossRef Full Text | Google Scholar

Kettunen, J., Demirkan, A., Würtz, P., Draisma, H. H. M., Haller, T., Rawal, R., et al. (2016). Genome-wide Study for Circulating Metabolites Identifies 62 Loci and Reveals Novel Systemic Effects of LPA. Nat. Commun. 7, 11122. doi:10.1038/ncomms11122

PubMed Abstract | CrossRef Full Text | Google Scholar

Khera, A. V., and Kathiresan, S. (2017). Genetics of Coronary Artery Disease: Discovery, Biology and Clinical Translation. Nat. Rev. Genet. 18, 331–344. doi:10.1038/nrg.2016.160

PubMed Abstract | CrossRef Full Text | Google Scholar

Kiefer, M. C., Tucker, J. E., Joh, R., Landsberg, K. E., Saltman, D., and Barr, P. J. (1991). Identification of a Second Human Subtilisin-like Protease Gene in thefes/fpsRegion of Chromosome 15. DNA Cel. Biol. 10, 757–769. doi:10.1089/dna.1991.10.757

CrossRef Full Text | Google Scholar

Klarin, D., Damrauer, S. M., Cho, K., Sun, Y. V., Teslovich, T. M., Honerlaw, J., et al. (2018). Genetics of Blood Lipids Among ∼300,000 Multi-Ethnic Participants of the Million Veteran Program. Nat. Genet. 50, 1514–1523. doi:10.1038/s41588-018-0222-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Klarin, D., Lynch, J., Aragam, K., Chaffin, M., Assimes, T. L., Huang, J., et al. (2019). Genome-wide Association Study of Peripheral Artery Disease in the Million Veteran Program. Nat. Med. 25, 1274–1279. doi:10.1038/s41591-019-0492-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Klimentidis, Y. C., Arora, A., Newell, M., Zhou, J., Ordovas, J. M., Renquist, B. J., et al. (2020). Phenotypic and Genetic Characterization of Lower LDL Cholesterol and Increased Type 2 Diabetes Risk in the UK Biobank. Diabetes 69, 2194–2205. doi:10.2337/db19-1134

PubMed Abstract | CrossRef Full Text | Google Scholar

Kobiyama, K., and Ley, K. (2018). Atherosclerosis. Circ. Res. 123 (10), 1118–1120. doi:10.1161/circresaha.118.313816

PubMed Abstract | CrossRef Full Text | Google Scholar

Komiyama, T., Coppola, J. M., Larsen, M. J., Van Dort, M. E., Ross, B. D., Day, R., et al. (2009). Inhibition of Furin/proprotein Convertase-Catalyzed Surface and Intracellular Processing by Small Molecules. J. Biol. Chem. 284, 15729–15738. doi:10.1074/jbc.m901540200

CrossRef Full Text | Google Scholar

Koyama, S., Ito, K., Terao, C., Akiyama, M., Horikoshi, M., Momozawa, Y., et al. (2020). Population-specific and Trans-ancestry Genome-wide Analyses Identify Distinct and Shared Genetic Risk Loci for Coronary Artery Disease. Nat. Genet. 52, 1169–1177. doi:10.1038/s41588-020-0705-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Larsson, S. C., Mason, A. M., Bäck, M., Klarin, D., Damrauer, S. M., Michaëlsson, K., et al. (2020). Genetic Predisposition to Smoking in Relation to 14 Cardiovascular Diseases. Eur. Heart J. 41, 3304–3310. doi:10.1093/eurheartj/ehaa193

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, T., Li, X., Feng, Y., Dong, G., Wang, Y., and Yang, J. (2020). The Role of Matrix Metalloproteinase-9 in Atherosclerotic Plaque Instability. Mediators Inflamm. 2020, 3872367. doi:10.1155/2020/3872367

PubMed Abstract | CrossRef Full Text | Google Scholar

Lin, J. S., Olson, C. M., Johnson, E. S., and Whitlock, E. P. (2013). The Ankle-Brachial index for Peripheral Artery Disease Screening and Cardiovascular Disease Prediction Among Asymptomatic Adults: a Systematic Evidence Review for the U.S. Preventive Services Task Force. Ann. Intern. Med. 159, 333–341. doi:10.7326/0003-4819-159-5-201309030-00007

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, M., Jiang, Y., Wedow, R., Li, Y., Brazel, D. M., Chen, F., et al. (2019). Association Studies of up to 1.2 Million Individuals Yield New Insights into the Genetic Etiology of Tobacco and Alcohol Use. Nat. Genet. 51, 237–244. doi:10.1038/s41588-018-0307-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Lozano, R., Naghavi, M., Foreman, K., Lim, S., Shibuya, K., Aboyans, V., et al. (2012). Global and Regional Mortality from 235 Causes of Death for 20 Age Groups in 1990 and 2010: a Systematic Analysis for the Global Burden of Disease Study 2010. Lancet 380 (9859), 2095–2128. doi:10.1016/S0140-6736(12)61728-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Machiela, M. J., and Chanock, S. J. (2018). LDassoc: an Online Tool for Interactively Exploring Genome-wide Association Study Results and Prioritizing Variants for Functional Investigation. Bioinformatics 34, 887–889. doi:10.1093/bioinformatics/btx561

PubMed Abstract | CrossRef Full Text | Google Scholar

Mahajan, A., Wessel, J., Willems, S. M., Zhao, W., Robertson, N. R., Chu, A. Y., et al. (2018). Refining the Accuracy of Validated Target Identification through Coding Variant fine-mapping in Type 2 Diabetes. Nat. Genet. 50, 559–571. doi:10.1038/s41588-018-0084-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Maier, M. T., Vilhelmsson, A., Louie, S. M., Vagena, E., Nomura, D. K., Koliwad, S. K., et al. (2018). Regulation of Hepatic Lipid Accumulation and Distribution by Agouti-Related Protein in Male Mice. Endocrinology 159, 2408–2420. doi:10.1210/en.2018-00040

PubMed Abstract | CrossRef Full Text | Google Scholar

Malik, R., Chauhan, G., Traylor, M., Sargurupremraj, M., Okada, Y., Mishra, A., et al. (2018). Multiancestry Genome-wide Association Study of 520,000 Subjects Identifies 32 Loci Associated with Stroke and Stroke Subtypes. Nat. Genet. 50, 524–537. doi:10.1038/s41588-018-0058-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Matsukura, M., Ozaki, K., Takahashi, A., Onouchi, Y., Morizono, T., Komai, H., et al. (2015). Genome-Wide Association Study of Peripheral Arterial Disease in a Japanese Population. PLoS One 10, e0139262. doi:10.1371/journal.pone.0139262

PubMed Abstract | CrossRef Full Text | Google Scholar

Nath, A. P., Ritchie, S. C., Grinberg, N. F., Tang, H. H.-F., Huang, Q. Q., Teo, S. M., et al. (2019). Multivariate Genome-wide Association Analysis of a Cytokine Network Reveals Variants with Widespread Immune, Haematological, and Cardiometabolic Pleiotropy. Am. J. Hum. Genet. 105, 1076–1090. doi:10.1016/j.ajhg.2019.10.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Ozkaramanli Gur, D., Guzel, S., Akyuz, A., Alpsoy, S., and Guler, N. (2018). The Role of Novel Cytokines in Inflammation: Defining Peripheral Artery Disease Among Patients with Coronary Artery Disease. Vasc. Med. 23, 428–436. doi:10.1177/1358863x18763096

PubMed Abstract | CrossRef Full Text | Google Scholar

Pashos, E. E., Park, Y., Wang, X., Raghavan, A., Yang, W., Abbey, D., et al. (2017). Large, Diverse Population Cohorts of hiPSCs and Derived Hepatocyte-like Cells Reveal Functional Genetic Variation at Blood Lipid-Associated Loci. Cell Stem Cell 20, 558–570. doi:10.1016/j.stem.2017.03.017

PubMed Abstract | CrossRef Full Text | Google Scholar

Perisic, L., Hedin, E., Razuvaev, A., Lengquist, M., Osterholm, C., Folkersen, L., et al. (2013). Profiling of Atherosclerotic Lesions by Gene and Tissue Microarrays Reveals PCSK6 as a Novel Protease in Unstable Carotid Atherosclerosis. Arterioscler Thromb. Vasc. Biol. 33, 2432–2443. doi:10.1161/atvbaha.113.301743

PubMed Abstract | CrossRef Full Text | Google Scholar

Pruim, R. J., Welch, R. P., Sanna, S., Teslovich, T. M., Chines, P. S., Gliedt, T. P., et al. (2010). LocusZoom: Regional Visualization of Genome-wide Association Scan Results. Bioinformatics 26, 2336–2337. doi:10.1093/bioinformatics/btq419

PubMed Abstract | CrossRef Full Text | Google Scholar

Pulit, S. L., Stoneman, C., Morris, A. P., Wood, A. R., Glastonbury, C. A., Tyrrell, J., et al. (2019). Meta-analysis of Genome-wide Association Studies for Body Fat Distribution in 694 649 Individuals of European Ancestry. Hum. Mol. Genet. 28, 166–174. doi:10.1093/hmg/ddy327

PubMed Abstract | CrossRef Full Text | Google Scholar

Purcell, S., Neale, B., Todd-Brown, K., Thomas, L., Ferreira, M. A. R., Bender, D., et al. (2007). PLINK: a Tool Set for Whole-Genome Association and Population-Based Linkage Analyses. Am. J. Hum. Genet. 81, 559–575. doi:10.1086/519795

CrossRef Full Text | Google Scholar

Riaz, H., Khan, M. S., Siddiqi, T. J., Usman, M. S., Shah, N., Goyal, A., et al. (2018). Association between Obesity and Cardiovascular Outcomes. JAMA Netw. Open 1, e183788. doi:10.1001/jamanetworkopen.2018.3788

PubMed Abstract | CrossRef Full Text | Google Scholar

Richardson, T. G., Sanderson, E., Palmer, T. M., Ala-Korpela, M., Ference, B. A., Davey Smith, G., et al. (2020). Evaluating the Relationship between Circulating Lipoprotein Lipids and Apolipoproteins with Risk of Coronary Heart Disease: A Multivariable Mendelian Randomisation Analysis. Plos Med. 17, e1003062. doi:10.1371/journal.pmed.1003062

PubMed Abstract | CrossRef Full Text | Google Scholar

Ruotsalainen, S. E., Partanen, J. J., Cichonska, A., Lin, J., Benner, C., Surakka, I., et al. (2021). An Expanded Analysis Framework for Multivariate GWAS Connects Inflammatory Biomarkers to Functional Variants and Disease. Eur. J. Hum. Genet. 29, 309–324. doi:10.1038/s41431-020-00730-8

CrossRef Full Text | Google Scholar

Ruth, K. S., Day, F. R., Tyrrell, J., Thompson, D. J., Wood, A. R., Mahajan, A., et al. (2020). Using Human Genetics to Understand the Disease Impacts of Testosterone in Men and Women. Nat. Med. 26, 252–258. doi:10.1038/s41591-020-0751-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Rykaczewska, U., Suur, B. E., Röhl, S., Razuvaev, A., Lengquist, M., Sabater-Lleal, M., et al. (2020). PCSK6 Is a Key Protease in the Control of Smooth Muscle Cell Function in Vascular Remodeling. Circ. Res. 126, 571–585. doi:10.1161/circresaha.119.316063

PubMed Abstract | CrossRef Full Text | Google Scholar

Siewert, K. M., and Voight, B. F. (2018). Bivariate Genome-wide Association Scan Identifies 6 Novel Loci Associated with Lipid Levels and Coronary Artery Disease. Circ. Genom. Precis Med. 11, e002239. doi:10.1161/CIRCGEN.118.002239

PubMed Abstract | CrossRef Full Text | Google Scholar

Sinnott-Armstrong, N., Tanigawa, Y., Amar, D., Mars, N., Benner, C., Aguirre, M., et al. (2021). Genetics of 35 Blood and Urine Biomarkers in the UK Biobank. Nat. Genet. 53, 185–194. doi:10.1038/s41588-020-00757-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Speliotes, E. K., Yerges-Armstrong, L. M., Wu, J., Hernaez, R., Kim, L. J., Palmer, C. D., et al. (2011). Genome-wide Association Analysis Identifies Variants Associated with Nonalcoholic Fatty Liver Disease that Have Distinct Effects on Metabolic Traits. Plos Genet. 7, e1001324. doi:10.1371/journal.pgen.1001324

PubMed Abstract | CrossRef Full Text | Google Scholar

Subkhangulova, A., Malik, A. R., Hermey, G., Popp, O., Dittmar, G., Rathjen, T., et al. (2018). SORCS1 and SORCS3 Control Energy Balance and Orexigenic Peptide Production. EMBO Rep. 19, e44810. doi:10.15252/embr.201744810

PubMed Abstract | CrossRef Full Text | Google Scholar

Sundaram, V., Bloom, C., Zakeri, R., Halcox, J., Cohen, A., Bowrin, K., et al. (2020). Temporal Trends in the Incidence, Treatment Patterns, and Outcomes of Coronary Artery Disease and Peripheral Artery Disease in the UK, 2006-2015. Eur. Heart J. 41, 1636–1649. doi:10.1093/eurheartj/ehz880

PubMed Abstract | CrossRef Full Text | Google Scholar

Testa, G., Staurenghi, E., Giannelli, S., Sottero, B., Gargiulo, S., Poli, G., et al. (2021). Up-regulation of PCSK6 by Lipid Oxidation Products: A Possible Role in Atherosclerosis. Biochimie 181, 191–203. doi:10.1016/j.biochi.2020.12.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Turcot, V., Lu, Y., Highland, H. M., Schurmann, C., Justice, A. E., Fine, R. S., et al. (2018). Protein-altering Variants Associated with Body Mass index Implicate Pathways that Control Energy Intake and Expenditure in Obesity. Nat. Genet. 50, 26–41. doi:10.1038/s41588-017-0011-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Van Der Harst, P., and Verweij, N. (2018). Identification of 64 Novel Genetic Loci Provides an Expanded View on the Genetic Architecture of Coronary Artery Disease. Circ. Res. 122, 433–443. doi:10.1161/circresaha.117.312086

PubMed Abstract | CrossRef Full Text | Google Scholar

Virani, S. S., Alonso, A., Benjamin, E. J., Bittencourt, M. S., Callaway, C. W., Carson, A. P., et al. (2020). Heart Disease and Stroke Statistics-2020 Update: A Report from the American Heart Association. Circulation 141 (9), e139–e596. doi:10.1161/CIR.0000000000000757

PubMed Abstract | CrossRef Full Text | Google Scholar

Vujkovic, M., Keaton, J. M., Lynch, J. A., Miller, D. R., Zhou, J., Tcheandjieu, C., et al. (2020). Discovery of 318 New Risk Loci for Type 2 Diabetes and Related Vascular Outcomes Among 1.4 Million Participants in a Multi-Ancestry Meta-Analysis. Nat. Genet. 52, 680–691. doi:10.1038/s41588-020-0637-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Ward, J., Tunbridge, E. M., Sandor, C., Lyall, L. M., Ferguson, A., Strawbridge, R. J., et al. (2020). The Genomic Basis of Mood Instability: Identification of 46 Loci in 363,705 UK Biobank Participants, Genetic Correlation with Psychiatric Disorders, and Association with Gene Expression and Function. Mol. Psychiatry 25, 3091–3099. doi:10.1038/s41380-019-0439-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Wootton, R. E., Richmond, R. C., Stuijfzand, B. G., Lawn, R. B., Sallis, H. M., Taylor, G. M. J., et al. (2020). Evidence for Causal Effects of Lifetime Smoking on Risk for Depression and Schizophrenia: a Mendelian Randomisation Study. Psychol. Med. 50, 2435–2443. doi:10.1017/s0033291719002678

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, J., Lee, S. H., Goddard, M. E., and Visscher, P. M. (2011). GCTA: a Tool for Genome-wide Complex Trait Analysis. Am. J. Hum. Genet. 88, 76–82. doi:10.1016/j.ajhg.2010.11.011

CrossRef Full Text | Google Scholar

Yengo, L., Sidorenko, J., Kemper, K. E., Zheng, Z., Wood, A. R., Weedon, M. N., et al. (2018). Meta-analysis of Genome-wide Association Studies for Height and Body Mass index in ∼700000 Individuals of European Ancestry. Hum. Mol. Genet. 27, 3641–3649. doi:10.1093/hmg/ddy271

PubMed Abstract | CrossRef Full Text | Google Scholar

Yuan, X., Waterworth, D., Perry, J. R. B., Lim, N., Song, K., Chambers, J. C., et al. (2008). Population-based Genome-wide Association Studies Reveal Six Loci Influencing Plasma Levels of Liver Enzymes. Am. J. Hum. Genet. 83, 520–528. doi:10.1016/j.ajhg.2008.09.012

CrossRef Full Text | Google Scholar

Zhao, W., Rasheed, A., Tikkanen, E., Lee, J.-J., Butterworth, A. S., Howson, J. M. M., et al. (2017). Identification of New Susceptibility Loci for Type 2 Diabetes and Shared Etiological Pathways with Coronary Heart Disease. Nat. Genet. 49, 1450–1457. doi:10.1038/ng.3943

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhu, Z., Guo, Y., Shi, H., Liu, C.-L., Panganiban, R. A., Chung, W., et al. (2020). Shared Genetic and Experimental Links between Obesity-Related Traits and Asthma Subtypes in UK Biobank. J. Allergy Clin. Immunol. 145, 537–549. doi:10.1016/j.jaci.2019.09.035

CrossRef Full Text | Google Scholar