# Genome-Wide Association of Proprotein Convertase Subtilisin/Kexin Type 9 Plasma Levels in the ELSA-Brasil Study

Isabela Bensenor, Kallyandra Padilha, Isabella Ramos Lima, Raul Dias

Santos, Gilles Lambert, Stéphane Ramin-Mangata, Marcio Bittencourt,

Alessandra Goulart, Itamar Santos, Jose Mill, et al.

► To cite this version:

Isabela Bensenor, Kallyandra Padilha, Isabella Ramos Lima, Raul Dias Santos, Gilles Lambert, et al.. Genome-Wide Association of Proprotein Convertase Subtilisin/Kexin Type 9 Plasma Levels in the ELSA-Brasil Study. Frontiers in Genetics, 2021, 12, 10.3389/fgene.2021.728526. hal-04465439

HAL Id: hal-04465439

https://hal.science/hal-04465439v1

Submitted on 10 Oct 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Check for
updates

# Genome-Wide Association of Proprotein Convertase Subtilisin/Kexin Type 9 Plasma Levels in the ELSA-Brasil Study

Isabela Bensenor<fnref n="1" />, Kallyandra Padilha², Isabella Ramos Lima², Raul Dias Santos³, Gilles Lambert⁴, Stéphane Ramin-Mangata⁴, Marcio S Bittencourt³, Alessandra C Goulart<fnref n="1" />, Itamar S. Santos<fnref n="1" />, Jose G Mill⁵, Jose E Krieger², Paulo A. Lotufo<fnref n="1" /> and Alexandre C. Pereira²*

## OPEN ACCESS

**Edited by:**
Muhammad Abu-Elmagd,
King Abdulaziz University, Saudi Arabia

**Reviewed by:**
Massimiliano Ruscica,
University of Milan, Italy
Jianxiang Shi,
Zhengzhou University, China

*Correspondence:
Alexandre C. Pereira
alexandre.pereira@incor.usp.br

**Specialty section:**
This article was submitted to Human and Medical Genomics, a section of the journal Frontiers in Genetics

**Received:** 21 June 2021
**Accepted:** 16 August 2021
**Published:** 29 September 2021

**Citation:**
Bensenor I, Padilha K, Lima IR, Santos RD, Lambert G, Ramin-Mangata S, Bittencourt MS, Goulart AC, Santos IS, Mill JG, Krieger JE, Lotufo PA and Pereira AC (2021) Genome-Wide Association of Proprotein Convertase Subtilisin/Kexin Type 9 Plasma Levels in the ELSA-Brasil Study. *Front. Genet.* 12:728526.
doi: 10.3389/fgene.2021.728526

Pharmacological inhibition of PCSK9 (proprotein convertase subtilisin/kexin type 9) is an established therapeutic option to treat hypercholesterolemia, and plasma PCSK9 levels have been implicated in cardiovascular disease incidence. A number of genetic variants within the PCSK9 gene locus have been shown to modulate PCSK9 levels, but these only explain a very small percentage of the overall PCSK9 interindividual variation. Here we present data on the genetic association structure between PCSK9 levels and genome-wide genetic variation in a healthy sample from the general population. We performed a genome-wide association study of plasma PCSK9 levels in a sample of Brazilian individuals enrolled in the Estudo Longitudinal de Saude do Adulto cohort (n=810). Enrolled individuals were free from cardiovascular disease, diabetes and were not under lipid-lowering medication. Genome-wide genotyping was conducted using the Axiom_PMRA.r3 array, and imputation was performed using the TOPMED multi-ancestry sample panel as reference. Total PCSK9 plasma concentrations were determined using the Quantikine SPC900 ELISA kit. We observed two genome-wide significant loci and seven loci that reached the pre-defined value of *p* threshold of $1 \times 10^{-6}$. Significant variants were near KCNA5 and KCNA1, and LINC00353. Genetic variation at the PCSK9 locus was able to explain approximately 4% of the overall interindividual variations in PCSK9 levels. Colocalization analysis using eQTL data suggested RWDD3, ATXN7L1, KCNA1, and FAM177A1 to be potential mediators of some of the observed associations. Our results suggest that PCSK9 levels may be modulated by *trans* genetic variation outside of the PCSK9 gene and this may have clinical implications. Understanding both environmental and genetic predictors of PCSK9 levels may help identify new targets for cardiovascular disease treatment and contribute to a better assessment of the benefits of long-term PCSK9 inhibition.

<footnote n="1">¹Center for Clinical and Epidemiologic Research, University of São Paulo, São Paulo, Brazil, ²Laboratory of Genetics and Molecular Cardiology, Heart Institute (InCor), University of São Paulo Medical School Hospital, São Paulo, Brazil, ³Lipid Clinic, Heart Institute (InCor), University of São Paulo Medical School Hospital, São Paulo, Brazil, ⁴Inserm UMR 1188 DéTROI, Université La Réunion, Sainte Clotilde, France, ⁵Department of Physiological Sciences, Federal University of Espírito Santo, Vitória, Brazil</footnote>
Keywords: genetic variation, genome-wide association study, atherosclerosis, proprotein convertase subtilisin/kexin type 9, colocalization analysis

## INTRODUCTION

Proprotein convertase subtilisin/kexin type 9 (PCSK9) is a key modulator of LDL receptor (LDLR) degradation and, consequently, LDL-cholesterol (LDL-C) serum levels. Gain-of-function mutations in PCSK9 have been shown to cause familial hypercholesterolemia and increased cardiovascular risk (Hopkins et al., 2015). On the other hand, loss-of-function variants have been shown to associate with low LDL-C levels and reduced cardiovascular risk (Kent et al., 2017). Furthermore, plasma PCSK9 has been independently associated with other components of the lipid profile (Leander et al., 2016; Brumpton et al., 2019). As a result, pharmacological inhibition of PCSK9 became mainstream as a lipid reduction strategy (Reyes-Soffer et al., 2017).

Understanding the factors that modulate interindividual variability of PCSK9 plasma levels is important for the better understanding of individual responses to treatment as well as the identification of new targets for cardiovascular disease treatment. The use of unbiased genetic approaches has the potential to contribute to increase our understanding of the two.

We conducted a genome-wide association study (GWAS) in healthy individuals from the general population aiming at the identification of genetic variation associated to plasma PCSK9 levels.

## MATERIALS AND METHODS

### Study Population

The study sample belongs to the Estudo Longitudinal de Saude do Adulto, NCT02320461 (ELSA-Brasil). For the present analysis, we included 810 participants that have both PCSK9 plasma levels and genome-wide genotype information.

The ELSA-Brasil study design and cohort profile have been published elsewhere (Aquino et al., 2012). Briefly, ELSA-Brasil enrolled 15,105 civil servants living in six large Brazilian urban areas (Belo Horizonte, Porto Alegre, Rio de Janeiro, Salvador, Sao Paulo, and Vitoria), aged between 35 and 74 years at baseline. Information on sociodemographic, clinical history, family history of diseases, lifestyle factors, mental health, cognitive status, and occupational exposure was assessed from August 2008 to December 2010. Anthropometric, laboratory and imaging measurements were also obtained. In addition to baseline measurements, samples of plasma and DNA were collected and stored for further analysis at -80°C (Pereira et al., 2013). All participants signed an informed consent before enrollment. The study was conducted in accordance with the Declaration of Helsinki and was approved by the Research Ethics Committees and by the National Research Ethics Committee (CONEP).

Participants enrolled in the Sao Paulo site (5,061 people in total) without diabetes (exclusion criteria: fasting plasma glucose-FPG >126 mg/dl and/or 2-h post-load glucose >200 mg/dl

and/or history of treatment with oral anti-diabetic agents or insulin), without cardiovascular, renal or hepatic diseases (exclusion criteria: self-reported history of medical diagnosis of these pathologies), and who did not report prescription of lipid-lowering agents, were eligible for a PCSK9 ancillary and exploratory study. From the 1,751 randomly selected participants fulfilling the inclusion criteria for subsequent PCSK9 plasma concentration measurements (Ramin-Mangata et al., 2020), we used the 810 who had genome-wide genotype information for the present analysis.

### Biochemical Analyses

A 12-h fasting blood sample was drawn in the morning soon after arrival at the research clinic, following standardized procedures for sample collection and processing. A standardized 75 g oral glucose tolerance test was performed in all participants without known diabetes utilizing an anhydrous glucose solution. For measurement of fasting and post-load glucose, we used the hexokinase method (ADVIA 1200, Siemens); for fasting and post-load insulin, an immunoenzymatic assay, and for HbA1c, high-pressure liquid chromatography. Total cholesterol (TC), high-density lipoprotein-cholesterol (HDL-C), and triglycerides (TG) were measured with enzymatic colorimetric assays (ADVIA Chemistry). LDL-C was calculated using the Friedewald equation. When TG were ≥400 mg/dl, LDL-C was measured directly with an enzymatic colorimetric assay (ADVIA Chemistry).

Total PCSK9 plasma concentrations were determined using the Quantikine SPC900 ELISA kit (R&D Systems, Lille, France; Ramin-Mangata et al., 2020). Briefly, plasma samples were diluted 1: 20 in the calibrator diluent onto ELISA plates and incubated for 2 h on a plate shaker at 450 rpm. Wells were rinsed with wash buffer using an automated Hydroflex TECAN microplate washer. The detection HRP-conjugated antibody was added to each well and plates were incubated for 2 h at 450 rpm. Wells were rinsed. The TMB substrate solution was added to each well and plates were further incubated in the dark for 30 min at 450 rpm. Reactions were stopped by the addition of 0.2 N acid sulphuric solution. Absorbance was read at 450 nm with reference at 540 nm on an Infinite 200 pro TECAN plate-reader. The same experimenter (SR-M) performed the PCSK9 measurements on the same site and at the same time. The reported intra-assay precision coefficient of variation was 5.4%, and the minimum detectable dose of human PCSK9 ranged from 0.030-0.219 ng/ml.

### SNP Genotyping and Imputation

Genomic DNA extraction has been previously described (Chor et al., 2019). ELSA-Brasil DNA samples were genotyped using Axiom_PMRA.r3 array (ThermoFisher) and genotypes annotated using the Axiom_PMRA.na35.annot.db provided at the ThermoFisher site. Genotype calling was performed using
Affymetrix Power Tools. Initial VCF file containing 850,483 variants fulfilled all quality criteria.

Imputation was performed using the Haplotype Reference Consortium Michigan Imputation Server using the TOPMED reference haplotype panel as reference. After imputation data were exported in the standard PLINK format, downstream QC procedures and statistical analysis were conducted using the latest PLINK<fnref n="1" /> and R software packages,<fnref n="2" /> installed on a Linux-based computation resource. Imputation markers were kept if R2>0.3, and minor allele frequency (MAF)>0.01. A HWE value of *p* <1×10⁻²⁰ was used to control for potential genotyping clustering problems. Genetic population structure was studied through PCA analysis after LD-pruning of associated markers (see also Statistical Analysis section). A total of 11,524,071 SNPs were used for genome-wide analysis, 11,289,274 for autosomal, and 234,797 for X-chromosomal analyses.

## Colocalization Analysis

For colocalization analysis, we defined a window spanning 500 Kb center at the most associated variant in all regions classified as having a suggestive association signal. Information on all variants within this region was used for colocalization testing. We used the LocusFocus<fnref n="3" /> analytical approach for colocalization testing. Briefly, all genes residing in each selected region with their expression quantitative trait loci (eQTL) summary statistics available in GTEx were sequentially tested for colocalization with the results obtained for PCSK9 association. As reference LD structure, we used 1,000 genomes 2012 European LD matrix (our sample has approximately 80% European ancestry). Colocalization was tested against all 48 tissues available in GTEx and the most significant signal was selected.

## Statistical Analysis

PCSK9 levels were log-transformed for all analyses. Baseline categorical parameters are presented using frequencies (proportions), and continuous parameters are presented using mean±SD. Before GWAS, we adjusted a linear model for log (PCSK9) adjusting for age. The residuals of this model were used for GWAS as a continuous variable. Confounding effects for age, sex, smoking and BMI were later tested for all genome-wide and suggestive GWA hits.

Genome-wide association analyses were conducted using plink. We conducted two analyses – one without any further adjustment and one adjusting for the first four principal components. The threshold for genome-wide significance was set to *p* <5×10⁻⁸. Associations with *p* <1×10⁻⁶ were considered as suggestive and presented as a list of top SNPs.

Due to the high level of admixture and complex genetic population structure present in the Brazilian population, we conducted two different sensitivity analyses taking into consideration self-referred race and a particular individual

position in a PCA plot generated using the 2 first principal components. Briefly, for the self-referred race sensitivity analysis, association summary statistics were generated in each of three established subgroups: whites, blacks and browns ("pardos" in Portuguese). For the PCA-defined subgroup analysis, we used k-means clustering with *k*=3 and defined three different subgroups with higher European, African and Native-American ancestries. Meta-analysis used a fixed-effect model and was calculated using plink -meta-analysis routine.

Local association plots were created using LocusZoom (Pruim et al., 2010). Local linkage disequilibrium structure was determined using Haploview (Barrett, 2009).

Mediation analysis was conducted for selected loci. To select mark for a genetic risk score for plasma PCSK9 levels, we determined independently associated variants at the PCSK9 genomic locus (cis-pQTL) through fitting a multiple linear regression model using 20 nominally associated markers at this locus and a stepwise variable selection procedure. Genetic risk score was derived as the sum of weighted genotypes by their final regression coefficients.

# RESULTS

## Relationship Between Cardiovascular Risk Factors and Plasma Pcsk9

Clinical and laboratory characteristics of the ELSA-Brasil sample used in the present analysis are summarized in Table 1. Plasma PCSK9 levels were associated with TC (*p*=0.0006), TG (*p*=0.003), and LDL-C (*p*=0.003; Table 1).

## Genome-Wide Association Analysis of PCSK9 Plasma Levels

We performed a GWA of the age-adjusted residuals of the log transformed values of plasma PCSK9. In the primary analysis adjusted for the four first principal components, we identified two loci that reached the pre-defined genome-wide significant level of 5×10⁻⁸ (Figure 1). Notably, no significant genomic inflation was observed (lambda = 1.06). In addition, no significant difference was observed in the effect sizes and values of *p* in genome-wide significant and suggestive loci when running an unadjusted analysis (Supplementary Figure S1).

In addition, we explored the effect that genetic population structure might have on observed results. For this, we conducted two different sensitivity analyses. The first, involved a trans-ethnic meta-analysis after separating individuals into the three most commonly self-referred races in Brazil (i.e., Whites, Blacks and Browns; Supplementary Figure S2). In the second, we first derived clusters using a k-means cluster algorithm with *k*=3 and data from the 2 first PCs (Supplementary Figure S3). Samples were divided into 3 subgroups and analysis proceeded as described for the sensitivity analysis using self-referred race. Of note, the overall results were very similar to the overall analysis.

<footnote n="1">¹http://pngu.mgh.harvard.edu/_purcell/plink</footnote>

<footnote n="2">²http://www.r-project.org/</footnote>

<footnote n="3">³https://locusfocus.research.sickkids.ca/</footnote>
<table><caption>**TABLE 1** | Clinical and laboratory characteristics of studied subjects according to tertiles of plasma PCSK9 concentrations.</caption><thead><tr><th></th><th>1st<br>(n=270)</th><th>2nd<br>(n=270)</th><th>3rd<br>(n=270)</th><th>Overall<br>(n=810)</th></tr></thead><tbody><tr><td><strong>log(PCSK9; ng/ml)</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>Mean(SD)</td><td>5.4(±0.16)</td><td>5.7(±0.057)</td><td>6.0(±0.14)</td><td>5.7(±0.26)</td></tr><tr><td><strong>Sex</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>Male</td><td>123(46%)</td><td>129(48%)</td><td>124(46%)</td><td>376(46%)</td></tr><tr><td>Female</td><td>147(54%)</td><td>141(52%)</td><td>146(54%)</td><td>434(54%)</td></tr><tr><td><strong>Age(years)</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>Mean(SD)</td><td>50(±7.9)</td><td>50(±8.8)</td><td>52(±7.9)</td><td>51(±8.2)</td></tr><tr><td><strong>Race</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>Black</td><td>38(14%)</td><td>24(9%)</td><td>34(13%)</td><td>96(12%)</td></tr><tr><td>Mixed</td><td>54(20%)</td><td>63(23%)</td><td>61(23%)</td><td>178(22%)</td></tr><tr><td>White</td><td>159(59%)</td><td>163(60%)</td><td>158(59%)</td><td>480(59%)</td></tr><tr><td>Asian</td><td>12(4%)</td><td>15(6%)</td><td>12(4%)</td><td>39(5%)</td></tr><tr><td>Indigenous</td><td>3(1%)</td><td>2(1%)</td><td>4(1%)</td><td>9(1%)</td></tr><tr><td>Missing</td><td>4(1.5%)</td><td>3(1.1%)</td><td>1(0.4%)</td><td>8(1.0%)</td></tr><tr><td><strong>Smoking</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>Never smoker</td><td>147(54%)</td><td>162(60%)</td><td>133(49%)</td><td>442(55%)</td></tr><tr><td>Former smoker</td><td>65(24%)</td><td>74(27%)</td><td>94(35%)</td><td>233(29%)</td></tr><tr><td>Smoker</td><td>58(21%)</td><td>34(13%)</td><td>43(16%)</td><td>135(17%)</td></tr><tr><td><strong>BMI(kg/m²)</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>Mean(SD)</td><td>26(±4.6)</td><td>27(±4.8)</td><td>27(±4.9)</td><td>27(±4.8)</td></tr><tr><td><strong>Hypertension</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>No</td><td>210(78%)</td><td>201(74%)</td><td>197(73%)</td><td>608(75%)</td></tr><tr><td>Yes</td><td>60(22%)</td><td>69(26%)</td><td>73(27%)</td><td>202(25%)</td></tr><tr><td><strong>Glucose(mg/dl)</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>Mean(SD)</td><td>110(±6.5)</td><td>110(±6.3)</td><td>110(±6.7)</td><td>110(±6.5)</td></tr><tr><td><strong>HbA1c(%)</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>Mean(SD)</td><td>5.2(±0.46)</td><td>5.2(±0.47)</td><td>5.2(±0.46)</td><td>5.2(±0.46)</td></tr><tr><td><strong>Uric acid(mg/dl)</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>Mean(SD)</td><td>5.6(±1.4)</td><td>5.7(±1.5)</td><td>5.7(±1.6)</td><td>5.6(±1.5)</td></tr><tr><td><strong>Total cholesterol(mg/dl)</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>Mean(SD)</td><td>210(±37)</td><td>210(±38)</td><td>220(±38)</td><td><strong>220(±38)***</strong></td></tr><tr><td><strong>Triglycerides(mg/dl)</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>Mean(SD)</td><td>120(±65)</td><td>130(±76)</td><td>140(±88)</td><td>130(±77)**</td></tr><tr><td><strong>HDL-C(mg/dl)</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>Mean(SD)</td><td>55(±14)</td><td>56(±14)</td><td>56(±13)</td><td>56(±14)</td></tr><tr><td><strong>LDL-C(mg/dl)</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>Mean(SD)</td><td>130(±31)</td><td>130(±32)</td><td>140(±33)</td><td><strong>130(±32)**</strong></td></tr></tbody></table>

## Genome-Wide Significant Loci

In our main analysis, we observed two genome-wide significant loci and seven loci that reached the pre-defined value of *p* threshold of $1 \times 10^{-6}$ (Table 2).

The strongest associations with PCSK9 plasma levels were observed on chromosome 12p13.32, top lead SNP rs116367042 (value of *p* 5.97e-09). A regional association plot of the locus is shown in Figure 2. The closest gene is KCNA5 and significant eQTLs have been observed in the region for AKAP3, DYRK4, KCNA5, KCNA1, NDUFA9, and GALTN8. The region has been

described as associated with serum uric acid levels in a previous GWAS. All summary statistics from the main analysis can be found at Supplementary File 1.

The second genome-wide significant hit was observed in 13q31.2, in the region coding for LINC00353. A regional association plot of the locus is shown in Supplementary Figure S5. Notably, only a single marker was associated with PCSK9 levels at this locus, reducing, thus, its credibility.

## Suggestive Loci

Using a pre-defined suggestive significance threshold of $1 \times 10^{-6}$ we identified additional 7 loci (Supplementary Table S2; Supplementary Figure S5). Those loci were on chr1p32.3 (nearest gene TXNDC12.), chr1p21.3 (nearest gene RWDD3), chr1q32.1 (nearest gene MAPKAPK2), chr7q22.3 (nearest gene ATXN7L1), chr9q31.3 (nearest gene MUSK), chr14q13.2 (nearest gene KIAA0391), and chr18q12.2 (nearest gene KIAA1328).

## PCSK9 Locus Association Structure

Previous GWAS and candidate-gene association studies have observed significant associations between PCSK9, LDL-C, and TC levels and genetic variants at the PCSK9 locus. Here we extend these observations using a multi-ethnic sample (Supplementary Figure S6). Of note, stronger associations are located at the 3' region of PCSK9 and within the nearby USP24 gene. Interestingly, previous studies in Europeans, African and other admixed samples have also described stronger associations for total cholesterol and LDL levels at this same region.

Linkage disequilibrium of the PCSK9 locus was resolved in four main haplotype blocks (Supplementary Table S1). Fifty-seven markers were nominally associated with PCSK9 levels being the most associated rs505151, rs662145, rs487230, and rs555687. Tagging associated SNPs in the PCSK9 locus, we were able to reduce the number of associated variants from 57 to 20, capturing 100% of the initial variation.

Using information from all 20 tagged markers and a stepwise regression approach, we were able to derive a PCSK9 instrumental variable made of four independently associated markers at the PCSK9 locus (cis-pQTLs; Supplementary Table S1). The R-squared for the multiple regression model containing all 4 markers was 0.036. Of particular importance, a model containing independently associated markers, BMI, age and smoking status, although highly significant ($p=5.186e-08$) was only able to explain 5.2% of the overall variation in PCSK9 plasma levels in our sample, genetic information being the variable with the highest effect size in our model.

## Colocalization Analysis of Associated Loci

Finally, we studied the colocalization pattern between the identified loci and expression traits of the genes in the vicinity of the association signal. For this, we used data from all the available tissues in the GTEx database. Colocalization analysis suggested that RWDD3, ATXN7L1, KCNA1, and FAM177A1 are potential mediators of the observed associations on chromosomes 1, 7, 12, and 14, respectively (Supplementary Table S2).

<footnote>**value of *p* 0.001–0.01; ***$p<0.001$.</footnote>
<figure><img src="image_3.png" /><figcaption>**TABLE 2** | Genome-wide and suggestive associated loci.</figcaption></figure>

<table><thead><tr><th>Marker</th><th>Chromosome</th><th>Position (Hg19)</th><th>Ref allele</th><th>Alt allele</th><th>value of p</th></tr></thead><tbody><tr><td>rs35120342</td><td>1</td><td>52,497,077</td><td>G</td><td>T</td><td>1.03E-07</td></tr><tr><td>rs12125618</td><td>1</td><td>95,781,754</td><td>T</td><td>G</td><td>4.28E-07</td></tr><tr><td>rs17015194</td><td>1</td><td>206,911,203</td><td>C</td><td>G</td><td>8.80E-07</td></tr><tr><td>rs112386665</td><td>7</td><td>105,335,698</td><td>C</td><td>A</td><td>7.77E-08</td></tr><tr><td>rs4574919</td><td>9</td><td>113,501,372</td><td>T</td><td>G</td><td>7.18E-08</td></tr><tr><td>rs116367042</td><td>12</td><td>5,109,095</td><td>T</td><td>C</td><td>5.97E-09</td></tr><tr><td>rs9555910</td><td>13</td><td>90,188,022</td><td>T</td><td>G</td><td>2.00E-08</td></tr><tr><td>rs10444669</td><td>14</td><td>35,659,171</td><td>T</td><td>A</td><td>4.91E-07</td></tr><tr><td>rs12457651</td><td>18</td><td>34,790,416</td><td>T</td><td>C</td><td>5.24E-07</td></tr></tbody></table>

## DISCUSSION

PCSK9 is a serine protease with protein-protein interaction with the LDL receptor that has both genetic and clinical validation (Petrilli et al., 2020). PCSK9 binds to the LDLR and is thought to reduce the recycling of these proteins from the cell surface (sending them to lysosomes instead), inhibiting LDL-particle removal from the extracellular fluid (Deng et al., 2019). Blocking PCSK9 can lower blood LDL-C concentrations, and low PCSK9 levels are associated with lower LDL-C levels and reduced incidence of atherosclerotic cardiovascular disease. Despite the

elusive importance of PCSK9 in lipoprotein homeostasis, few studies have analyzed PCSK9 plasma levels as a function of global genetic variation (Paquette and Baass, 2018; Pott et al., 2018). Understanding the genetic architecture that modulates PCSK9 levels may help dissect the mechanisms by which PCSK9 inhibition improves vascular function and overall cardiovascular morbidity and mortality.

It is assumed that PCSK9 modulates cardiovascular risk through cholesterol levels, more specifically LDL-C levels. Indeed, pharmacological inhibition of PCSK9 leads to significant decreases in LDL-C and reduction in the incidence of cardiovascular
<figure><img src="image_4.png" /><figcaption>**FIGURE 2** | Local association plot for rs116367042.</figcaption></figure>

events. However, it is not known whether PCSK9 has other actions independent of plasmatic LDL-C levels (Cesaro et al., 2020), and the association between PCSK9 inhibition and inflammatory markers has not been as consistent as its association to lipid levels (Ruscica et al., 2019). For instance, PCSK9 is substantially expressed in arterial walls and macrophages (Cariou et al., 2016), and it is possible that its actions in these cells are not directly linked to LDL-C metabolism. PCSK9 has also been shown to be associated with metabolic factors other than lipoproteins. It is positively associated with albumin, liver enzymes (ALT, ALP, AST, GGT) and with hepatic steatosis, although whether this association is confounded by or mediated by LDL-C is still unclear (Paquette et al., 2020). Finally, PCSK9 is also expressed in the intestine, endocrine pancreas and brain, and non-lipid-lowering effects of PCSK9 inhibition could also be linked to platelet activation (Qi et al., 2021), cell proliferation and apoptosis (Macchi et al., 2021).

Importantly, there is great interindividual variation in both PCSK9 levels and response to PCSK9 inhibitors (Chasman

et al., 2012; Qamar et al., 2019; Ramin-Mangata et al., 2020). In addition, only about 20% of circulating PCSK9 variance can be explained by clinical variables. Previously identified genetic variation only adds less than 5% to this figure, almost all of it from eQTL and pQTL within the PCSK9 locus itself (Anderson et al., 2014; Lim, 2017).

We conducted a GWAS study aiming at identifying genetic determinants of PCSK9 plasma levels. To our knowledge, this is the second GWAS conducted for PCSK9 levels and the first using a sample from a multi-ethnic population. Despite the relatively small sample size, we were able to observe two genome-wide significant association loci and a number of loci with suggestive association signals. In addition, we have confirmed the previously described association between PCSK9 levels and common genetic variation at the *PCSK9* locus (Pott et al., 2018). The most interesting observed genome-wide significant locus was at Chr12 within the *KCNA* gene cluster (Figure 2). Colocalization analysis was able to detect a significant colocalization signal with the *KCNA1* gene expression profile
in adipose subcutaneous tissue. Mutations in KCNA1 have been shown to cause episodic ataxia/myokymia syndrome type 1. The gene is lowly expressed in the adipose tissue and liver; and no metabolic phenotype has been associated with manipulations in KCNA1. Despite not being able to replicate this finding in an independent GWAS, further work characterizing the role of genetic variants nearby KCNA1 in PCSK9 levels is warranted.

We are not the first to describe genome-wide significant variants associated with plasma PCSK9 levels outside of the PCSK9 gene. Pott et al. in a GWAS conducted in 3290 individuals from the LIFE-Heart cohort identified variations within the FBXL18 gene to be associated with PCSK9 levels (Pott et al., 2018). We did not identify any association in this region and together with the low imputation quality the authors of this previous GWAS described, we suggest the association between FBXL18 and PCSK9 to be targeted in future studies aiming at clarifying the role of this locus on potentially regulating PCSK9 serum levels.

Plasma PCSK9 levels have been associated with several cardiovascular and metabolic risk factors (Caselli et al., 2019). Notwithstanding the understanding of PCSK9 mechanism at the molecular level, completely understanding the directionality of the associations between PCSK9 levels and other metabolic traits has been ill explored. In fact, most studies assume that PCSK9 is associated with lipid levels because of the interaction between PCSK9 and LDLR at the molecular level. However, it is unknown whether different predictors of PCSK9 levels are indeed associated with the same degree of increased cardiovascular risk. Indeed, recent experimental and clinical studies have also reported that higher circulating PCSK9 levels contributed to coronary atherosclerosis by enhancing the expression of pro-inflammatory genes, promoting apoptosis of human endothelial cells and activating platelet reactivity (Ricci et al., 2018; Yurtseven et al., 2020). The causal directionality of these associations, however, has not been fully explored and neither has the relationship of predictors of PCSK9 interindividual variation and clinically actionable management strategies.

By identifying potentially *trans* associations with PCSK9 levels, our data give rise to the possibility of a more complex mechanism, where different genetic factors may modulate PCSK9. It remains to be determined if PCSK9 levels driven by *trans* genetic factors carry the same increased risk of cardiovascular disease as PCSK9 levels determined by genetic variation at the PCSK9 locus. It is important to note that similar *trans* mechanisms have been identified for other important genes regulating lipid metabolism, such as LPA (Li et al., 2015). In summary, our data suggest that PCSK9 levels may be modulated by upstream targets other than genetic variation in the PCSK9 gene, which are well-known proxies for PCSK9 levels. The lack of large-effect size loci that modulate PCSK9 serum levels also points to the possibility that interindividual variation in PCSK9 is mostly a function of epigenetic modulation and not of a polygenic component. Indeed, several studies have shown that the PCSK9 promoter is dynamically methylated following several different exposures.

This study has some potential limitations. First and foremost, we have not been able to find a suitable replication sample for our GWAS results. The observed genome-wide significant loci still need to be replicated in an independent sample to be, in effect, taken as drivers of PCSK9 plasma levels (Trinder and Brunham, 2021). In addition, the reduced sample size of our study may have prevented us to identify other genome-wide significant loci with decreased effect size. In fact, post hoc calculation of our statistical power assuming a MAF of 0.2 was 0.73 to detect a difference of 0.5 standard deviation in the mean values of genotype groups, but only 0.11 to detect the same, different for alleles at a MAF of 0.1. Mendelian randomization analysis using our sample lacked the necessary statistical power to derive robust conclusions regarding the causality of *trans* PCSK9 variants and coronary artery disease or even LDL-C levels. Finally, lack of data on hsCRP and lipoprotein (a) precluded efforts in trying to understand the relevance of the described genetic associations in these variables. These aspects should be better defined in further studies.

## CONCLUSION

In conclusion, we describe new genome-wide significant loci associated with PCSK9 plasma levels in a sample from a healthy population. Our results suggest that PCSK9 levels may be modulated by *trans* genetic variation outside of the PCSK9 gene. Understanding both environmental and genetic predictors of PCSK9 levels may help identify new targets for cardiovascular disease treatment and contribute to better assessment of the benefits of long-term PCSK9 inhibition.

## DATA AVAILABILITY STATEMENT

All summary statistics from the main analysis are publicly available as Supplementary Material to this manuscript (Supplementary File Data Sheet 1). The data that support secondary findings of this study are available from ELSA-Brasil study on reasonable request. Requests for access to more detailed summary statistics, replication results, and analytic methods will be considered by the authors.

## ETHICS STATEMENT

The studies involving human participants were reviewed and approved by Research Ethics Committees and by the National Research Ethics Committee (CONEP). The patients/participants provided their written informed consent to participate in this study.

## AUTHOR CONTRIBUTIONS

IB conceptualized and designed the study, researched data, and wrote the manuscript. KP, IL, RS, GL, MB, AG, IS, JM,
and JK acquired data, assisted with the data analysis, and edited the manuscript. SR-M performed experiments and analysed and interpreted the results. PL and AP supervised the work, performed statistical analyses, and wrote the manuscript. All authors reviewed the manuscript.

## FUNDING

RS is recipient of a scholarship from Conselho Nacional de Pesquisa e Desenvolvimento Tecnológico, Brazil (CNPq) #303734/2018-3. MB is recipient of a scholarship from Conselho Nacional de Pesquisa e Desenvolvimento Tecnológico, Brazil (CNPq) #310255/2018-0. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. The ELSA-Brasil baseline study

## REFERENCES

Anderson, J. M., Cerda, A., Hirata, M. H., Rodrigues, A. C., Dorea, E. L., Bernik, M. M. S., et al. (2014). Influence of PCSK9 polymorphisms on plasma lipids and response to atorvastatin treatment in Brazilian subjects. *J. Clin. Lipidol.* 8, 256-264. doi: 10.1016/j.jacl.2014.02.008

Aquino, E. M. L., Barreto, S. M., Bensenor, I. M., Carvalho, M. S., Chor, D., Duncan, B. B., et al. (2012). Brazilian longitudinal study of adult health (ELSA-Brasil): objectives and design. *Am. J. Epidemiol.* 175, 315-324. doi: 10.1093/aje/kwr294

Barrett, J. C. (2009). Haploview: visualization and analysis of SNP genotype data. *Cold Spring Harb. Protoc.* 2009:pdb.ip71. doi: 10.1101/pdb.ip71

Brumpton, B. M., Fritsche, L. G., Zheng, J., Nielsen, J. B., Mannila, M., Surakka, I., et al. (2019). Variation in serum PCSK9 (Proprotein convertase Subtilisin/Kexin type 9), cardiovascular disease risk, and an investigation of potential unanticipated effects of PCSK9 inhibition. *Circ. Genomic Precis. Med.* 12:e002235. doi: 10.1161/CIRCGEN.118.002335

Cariou, B., Ding, Z., and Mehta, J. L. (2016). PCSK9 and atherosclerosis: Beyond LDL-cholesterol lowering. *Atherosclerosis* 253, 275-277. doi: 10.1016/j. atherosclerosis.2016.08.007

Caselli, C., Del Turco, S., Ragusa, R., Lorenzoni, V., De Graaf, M., Basta, G., et al. (2019). Association of PCSK9 plasma levels with metabolic patterns and coronary atherosclerosis in patients with stable angina. *Cardiovasc. Diabetol.* 18:144. doi: 10.1186/s12933-019-0949-3

Cesaro, A., Bianconi, V., Gragnano, F., Moscarella, E., Fimiani, F., Monda, E., et al. (2020). Beyond cholesterol metabolism: The pleiotropic effects of proprotein convertase subtilisin/kexin type 9 (PCSK9). Genetics, mutations, expression, and perspective for long-term inhibition. *Biofactors* 46, 367-380. doi: 10.1002/biof.1619

Chasman, D. I., Giulianini, F., MacFadyen, J., Barratt, B. J., Nyberg, F., and Ridker, P. M. (2012). Genetic determinants of statin-induced low-density lipoprotein cholesterol reduction. *Circ. Cardiovasc. Genet.* 5, 257-264. doi: 10.1161/CIRCGENETICS.111.961144

Chor, D., Pereira, A., Pacheco, A. G., Santos, R. V., Fonseca, M. J. M., Schmidt, M. I., et al. (2019). Context-dependence of race self-classification: results from a highly mixed and unequal middle-income country. *PLoS One* 14:e0216653. doi: 10.1371/journal.pone.0216653

Deng, S., Alabi, A., Gu, H., Adijiang, A., Qin, S., and Zhang, D. (2019). Identification of amino acid residues in the ligand binding repeats of LDL receptor important for PCSK9 binding. *J. Lipid Res.* 60, 516-527. doi: 10.1194/jlr.M089193

Hopkins, P. N., Defesche, J., Fouchier, S. W., Bruckert, E., Luc, G., Cariou, B., et al. (2015). Characterization of autosomal dominant hypercholesterolemia caused by PCSK9 gain of function mutations and its specific treatment With Alirocumab, a PCSK9 monoclonal antibody. *Circ. Cardiovasc. Genet.* 8, 823-831. doi: 10.1161/CIRCGENETICS.115.001129

was supported by the Brazilian Ministry of Health (Science and Technology Department) and the Brazilian Ministry of Science and Technology (Financiadora de Estudos e Projetos and CNPq National Research Council).

## ACKNOWLEDGMENTS

The authors thank the staff and participants of the ELSA-Brasil for their important contributions.

## SUPPLEMENTARY MATERIAL

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

Kent, S. T., Rosenson, R. S., Avery, C. L., Chen, Y.-D. I., Correa, A., Cummings, S. R., et al. (2017). PCSK9 loss-of-function variants, low-density lipoprotein cholesterol, and risk of coronary heart disease and stroke. *Circ. Cardiovasc. Genet.* 10:e001632. doi: 10.1161/CIRCGENETICS.116.001632

Leander, K., Målarstig, A., van't Hooft, F. M., Hyde, C., Hellénius, M.-L., Troutt, J. S., et al. (2016). Circulating Proprotein convertase Subtilisin/Kexin type 9 (PCSK9) predicts future risk of cardiovascular events independently of established risk FactorsCLINICAL PERSPECTIVE. *Circulation* 133, 1230-1239. doi: 10.1161/CIRCULATIONAHA.115.018531

Li, J., Lange, L. A., Sabourin, J., Duan, Q., Valdar, W., and Willis, M. S. (2015). Genome- and exome-wide association study of serum lipoprotein (a) in the Jackson heart study. *J. Hum. Genet.* 60, 755-761. doi: 10.1038/jhg.2015.107

Lim, G. B. (2017). HMGCR and PCSK9 variants and cardiovascular risk. *Nat. Rev. Cardiol.* 14:66. doi: 10.1038/nrcardio.2016.204

Macchi, C., Ferri, N., Sirtori, C. R., Corsini, A., Banach, M., and Ruscica, M. (2021). Proprotein convertase Subtilisin/Kexin type 9: A view beyond the canonical cholesterol-lowering impact. *Am. J. Pathol.* 191, 1385-1397. doi: 10.1016/j.ajpath.2021.04.016

Paquette, M., and Baass, A. (2018). Genetic regulation of PCSK9 (Proprotein convertase Subtilisin/Kexin type 9) plasma levels. *Circ. Genomic Precis. Med.* 11:e002165. doi: 10.1161/CIRCGEN.118.002165

Paquette, M., Gauthier, D., Chamberland, A., Prat, A., De Lucia Rolfe, E., Rasmussen, I. J., et al. (2020). Circulating PCSK9 is associated with liver biomarkers and hepatic steatosis. *Clin. Biochem.* 77, 20-25. doi: 10.1016/j.clinbiochem.2020.01.003

Pereira, A. C., Bensenor, I. M., Fedeli, L. M., Castilhos, C., Vidigal, P. G., Maniero, V., et al. (2013). Delineamento e implementacao do biobanco do ELSA-Brasil: estudo prospectivo na populacao brasileira. *Rev. Saude Publica* 47, 72-78. doi: 10.1590/S0034-8910.2013047003822

Petrilli, W. L., Adam, G. C., Erdmann, R. S., Abeywickrema, P., Agnani, V., Ai, X., et al. (2020). From screening to targeted degradation: strategies for the discovery and optimization of small molecule ligands for PCSK9. *Cell Chem. Biol.* 27, 32.e3-40.e3. doi: 10.1016/j.chembiol.2019.10.002

Pott, J., Schlegel, V., Teren, A., Horn, K., Kirsten, H., Bluecher, C., et al. (2018). Genetic regulation of PCSK9 (Proprotein convertase Subtilisin/Kexin type 9) plasma levels and its impact on atherosclerotic vascular disease phenotypes. *Circ. Genomic Precis. Med.* 11:e001992. doi: 10.1161/CIRCGEN.117.001992

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

Qamar, A., Giugliano, R. P., Keech, A. C., Kuder, J. F., Murphy, S. A., Kurtz, C. E., et al. (2019). Interindividual variation in low-density lipoprotein cholesterol level reduction With Evolocumab. *JAMA Cardiol.* 4:59. doi: 10.1001/jamacardio.2018.4178

Qi, Z., Hu, L., Zhang, J., Yang, W., Liu, X., Jia, D., et al. (2021). PCSK9 (Proprotein convertase Subtilisin/Kexin 9) enhances platelet activation,
thrombosis, and myocardial infarct expansion by binding to platelet CD36. *Circulation* 143, 45-61. doi: 10.1161/CIRCULATIONAHA.120.046290

Ramin-Mangata, S., Wargny, M., Pichelin, M., Le May, C., Thédrez, A., Blanchard, V., et al. (2020). Circulating PCSK9 levels are not associated with the conversion to type 2 diabetes. *Atherosclerosis* 293, 49-56. doi: 10.1016/j.atherosclerosis.2019.11.027

Reyes-Soffer, G., Pavlyha, M., Ngai, C., Thomas, T., Holleran, S., Ramakrishnan, R., et al. (2017). Effects of PCSK9 inhibition With Alirocumab on lipoprotein metabolism in healthy humans. *Circulation* 135, 352-362. doi: 10.1161/CIRCULATIONAHA.116.025253

Ricci, C., Ruscica, M., Camera, M., Rossetti, L., Macchi, C., Colciago, A., et al. (2018). PCSK9 induces a pro-inflammatory response in macrophages. *Sci. Rep.* 8:2267. doi: 10.1038/s41598-018-20425-x

Ruscica, M., Tokgözoğlu, L., Corsini, A., and Sirtori, C. R. (2019). PCSK9 inhibition and inflammation: A narrative review. *Atherosclerosis* 288, 146-155. doi: 10.1016/j.atherosclerosis.2019.07.015

Trinder, M., and Brunham, L. R. (2021). Polygenic scores for dyslipidemia: the emerging genomic model of plasma lipoprotein trait inheritance. *Curr. Opin. Lipidol.* 32, 103-111. doi: 10.1097/MOL.0000000000000737

Yurtseven, E., Ural, D., Baysal, K., and Tokgözoğlu, L. (2020). An update on the role of PCSK9 in atherosclerosis. *J. Atheroscler. Thromb.* 27, 909-918. doi: 10.5551/jat.55400

**Conflict of Interest:** RS has received honoraria related to consulting, research and/or speaker activities from: Abbott, Ache, Amgen, Astra Zeneca, Esperion,

EMS, Kowa, Libbs, Novo-Nordisk, Merck, MSD, Pfizer, PTC pharmaceuticals and Sanofi/Regeneron. MB has received honoraria related to consulting, research and/or speaker activities from: Boston Scientific, EMS, GE HealthCare, Novo-Nordisk, and Sanofi/Regeneron.

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.

The reviewer MR declared a past co-authorship with one of the authors RS to the handling editor.

**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.

Copyright © 2021 Bensenor, Padilha, Lima, Santos, Lambert, Ramin-Mangata, Bittencourt, Goulart, Santos, Mill, Krieger, Lotufo and Pereira. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.