1. Introduction
Atherosclerosis is a complex disease, which pathogenesis involving numerous pathological processes increased plasma cholesterol, its deposition in the arterial wall, endothelial dysfunction and vascular remodeling. These processes are influenced both by environmental and genetic factors [
1,
2]. It is widely agreed that dyslipidemia plays a crucial role in the development of atherosclerosis and coronary artery disease (CAD) [
2]. Lipid metabolism abnormalities and coronary heart disease are determined by substantial effects of genetic factors [
3,
4,
5]. According to the GWAS Catalog (
https://www.ebi.ac.uk/gwas/home), approximately 1700 genome-wide association studies (GWAS) on over 450 cardiovascular phenotypes have been conducted, yielding substantial insights into the genetic etiology of cardiovascular disease and related traits [
6]. Numerous GWAS including large-scale studies of the Global Lipids Genetics Consortium have discovered loci associated with plasma lipids [
7,
8,
9,
10,
11,
12]. GWAS were also conducted to find genetic determinants of statin pharmacogenetics, including JUPITER, the largest study of lipid-lowering effects of rosuvastatin [
13,
14,
15]. Statins are the first line of lipid-lowering drugs widely used in both primary and secondary prevention of CAD. For one of the most popular statins, rosuvastatin, it was shown that it could reduce both blood lipid levels, including total cholesterol (TC), LDL-C (low-density lipoprotein cholesterol), triglycerides (TG), and carotid intima-media thickness (CIMT) [
16].
Carotid intima-media thickness is a surrogate marker for the presence and progression of atherosclerosis, useful for evaluating risk and incidence of cardiovascular disease [
17]. A causal relationship between LDL-C and CIMT has been observed [
18], and CIMT is being considered as a marker of preclinical atherosclerosis [
18,
19]. A growing body of evidence has been provided to consider CIMT as a predictive marker for coronary artery atherosclerosis, its severity, and the extent of plaque burden [
20,
21,
22,
23]. Importantly, rosuvastatin may reduce not only CIMT but also atherosclerotic plaque growth [
24,
25]. It is impossible to claim definitely whether hypolipidemic effects of rosuvastatin are solely responsible for the regression of atherosclerosis, because pleiotropic effects of the drug may contribute to in arterial wall changes through inflammation reduction [
25,
26]. The mechanisms by which rosuvastatin is responsible for vascular wall changes may be unraveled by pharmacogenetic studies; however, a majority of such studies have investigated only lipid-lowering effects of the drug [
13,
14,
27,
28]. Our research team has conducted several studies on the effect of gene polymorphisms on plasma lipids and CIMT, as well as the dynamics of their changes during rosuvastatin lipid-lowering therapy [
29,
30,
31]. It has been discovered that some lipid-associated GWAS loci contribute to lipid-lowering effects associated with rosuvastatin therapy and determine CIMT regression. Without a doubt, CIMT seems to be a surrogate marker of atherosclerosis in CAD patients, despite it correlates with disease severity [
20,
21,
22,
23]. However, in terms of clinical application, CIMT regression during rosuvastatin therapy may mirror vascular wall changes and be indicative of drug efficacy. The aim of our pharmacogenetic study was to identify the effects of fifteen lipid-associated GWAS loci such as rs4846914 (
GALNT2), rs11220463 (
ST3GAL4), rs881844 (
STARD3), rs1689800 (
ZNF648), rs12328675 (
COBLL1), rs9987289 (
PPP1R3B), rs55730499 (
LPA), rs3136441 (
F2), rs6065906 (
PLTP), rs838880 (
SCARB1), rs386000 (
LILRA3), rs1883025 (
ABCA1), rs3764261(
CETP), rs217406 (
NPC1L1), and rs16942887 (
PSKH1) on plasma lipids and CIMT in CAD patients taking rosuvastatin. These loci have recently been observed to be associated with CAD susceptibility and CIMT [
32].
2. Materials and Methods
Study design. The present research was designed as a prospective study, in which the participants received rosuvastatin at an initial dose of 5 mg. Then, in 4 weeks, we controlled the levels of LDL-C. In cases where the levels were below the target level (1.8 mmol/l), patients continued taking the dose, which let them attain the target level. If it was not attained, the dose of rosuvastatin was gradually increased to 10, 20, and finally the maximum dose of 40 mg. The control of LDL-C levels was performed every 4 weeks after the dose increase. For the association analysis, we used the change (delta, Δ) in TC, LDL-C, TG, and CIMT values over the first 6 and 12 months of the study.
Study participants. The patient cohort comprised 116 Russian patients with the diagnosis of CAD, stable angina pectoris grade II–III. The study participants were described in detail in our previous work [
31]. In brief, the cohort included 85 men (73%) and 31 postmenopausal women (27%), with the mean age of 61.0 ± 7.25 years. Approximately 58% of CAD patients experienced a myocardial infarction in the past, and 97.5% of study participants had arterial hypertension. Diagnosis of CAD was confirmed by qualified cardiologists according to the Canadian Cardiovascular Society, ECG stress tests (treadmill test), and 24-h Holter’s ECG monitoring. TC levels higher than 4.0 mmol/l and LDL-C levels higher than 1.8 mmol/l were considered as criteria for dyslipidemia.
Biochemical and ultrasound investigations were described in our previous paper [
31]. Plasma lipid levels, except LDL-C, were detected using an automatic laboratory analyzer. LDL-C were calculated using Friedewald’s equation. CIMT measurement was performed in B-mode using the standard method at the distal third of the common carotid artery at a distance of 1–1.5 cm proximal to the bifurcation along the posterior wall [
33]. The following two parameters of carotid intima-media thickness were used in the study: mean CIMT (mean value of all measurements on the right and left sides) and maximal CIMT (the maximal value obtained from measurements on both sides with the further assessment of its change on the side where it was the highest).
Genetic analysis. Single nucleotide polymorphisms (SNPs) selected for this study have previously been found to be associated with plasma lipids in several genome-wide association studies [
7,
8,
9,
10,
11,
12] (GWAS Catalog (
https://www.ebi.ac.uk/gwas, accessed date 13 March 2023). As mentioned above, some of the SNPs have been recently found to be associated with CAD susceptibility and CIMT [
32], and we are interested in investigating whether these variants determine the lipid- and CIMT-lowering effects of rosuvastatin in CAD patients. The investigated SNPs included rs4846914 (
GALNT2), rs11220463 (
ST3GAL4), rs881844 (
STARD3), rs1689800 (
ZNF648), rs12328675 (
COBLL1), rs9987289 (
PPP1R3B), rs55730499 (
LPA), rs3136441 (
F2), rs6065906 (
PLTP), rs838880 (
SCARB1), rs386000 (
LILRA3), rs1883025 (
ABCA1), rs3764261 (
CETP), rs217406 (
NPC1L1), and rs16942887 (
PSKH1).
Extraction of DNA was performed from venous blood samples using phenol-chloroform method and precipitation with ethanol. Genotyping was performed using iPLEX technology on the MassARRAY 4 system (Agena Bioscience, San-Diego, CA, USA). Software MassARRAY Assay Design Suite was used to select a primer set and to design a multiplex panel for SNP genotyping. The sequences of primers are available on request.
Statistical analysis. The normality of the distribution of lipid and CIMT values was determined using the Kolmogorov-Smirnov and Shapiro-Wilk tests. As the trait distributions deviated from the normal one, the values were expressed as median and interquartile range (Me, Q1; Q3). Statistical package STATISTICA v13.0 (Statsoft, Tulsa, OK, USA) was utilized for descriptive statistics, distribution analysis, and determining the significance of lipid and CIMT changes. The significance of lipid and CIMT change during rosuvastatin therapy was tested by the Wilcoxon’s matched pairs test. The distribution of genotype frequencies according to Hardy–Weinberg equilibrium was assessed with Fisher’s exact test. For the association analysis, we used linear regression with adjustments for sex, age, body mass index (calculated as body mass in kilograms divided by height in meters in a square), and rosuvastatin dose, which allowed the participants to attain the target LDL-C level. Change in lipid levels was calculated as the difference between natural log-transformed on- and off-treatment levels divided by the natural log-transformed off-treatment level as described by Postmus et al. [
13]. We applied the same approach for the CIMT change, but before the logarithmic transformation, the CIMT data were multiplied by 10. Empirical P-values (Pperm) were calculated through the adaptive permutation procedure. We tested three genetic models for SNP-phenotype associations, such as additive, dominant, and recessive. The PLINK v1.92 software [
34] was used for all genetic calculations, including estimation of minor allele frequencies (MAF), tests for Hardy-Weinberg equilibrium (HWE), multiple regression analysis, and permutation procedures. When the p-value was less than 0.05, all of the results were declared statistically significant.
4. Discussion
The present study demonstrated, for the first time, the impact of lipid-associated GWAS loci on CIMT (
Figure 1) and plasma levels of TC, LDL-C, and TG in CAD patients during rosuvastatin therapy (
Figure 2).
In particular, we found associations of SNPs located at
ZNF648, LPA, ST3GAL4, PSKH1, GALNT2, COBLL1, PPP1R3B, and
STARD3 genes with reduction of CIMT. Pharmacogenetic associations of the variants in
ZNF648, LPA, ST3GAL4, PSKH1, PLTP, SCARB1, and
ABCA1 with the change in lipid levels were also found for the first time. Thus, lipid-associated GWAS loci are associated not only with lipids and the risk of coronary artery disease [
32], but are also responsible for the antiatherogenic effects of rosuvastatin after 6 and 12 months of lipid-lowering therapy.
Potential mechanisms by which polymorphic variants of the investigated lipid metabolism genes can impact the antiatherogenic effects of rosuvastatin therapy are of interest. The
ZNF648 gene encodes zinc finger protein 648, which may be involved in the DNA-templated transcription [
35]. In the present study rs1689800 polymorphism was associated with CIMT change. The possible mechanism of the effect of the variant on vascular wall may include the influence on oxidative stress (GLUL catabolizes glutamate, an amino acid required for glutathione synthesis), because such a mechanism is described for the nearest SNP rs10911021 located in the same genomic region [
36].
The
GALNT2 gene encodes N-acetylgalactosaminyltransferase 2, an enzyme involved in O-linked glycosylation of the substrates, including those regulating lipid metabolism such as apolipoprotein C-III (APOC-III), angiopoietin-related protein 3 (ANGPTL3) and phospholipid transfer protein (PLTP); in these ways it influences HDL-C and TG metabolism [
37,
38]. The rs4846914G allele has been linked to atherogenic changes in plasma lipid metabolism, such as an increase in TG and a decrease in HDL-C levels [
9,
10]. We found no associations of this variant with changes in any plasma lipid, but this SNP was significantly associated with a better CIMT regression during a 6-month period of rosuvastatin therapy. The presence of the link with CIMT change and absence of any association with lipid level change suggest that the pharmacogenetic effect on CIMT change might be mediated through non-lipid-related mechanisms, taking into account that SNP rs4846914 is known to be associated with endothelial function, serum levels of insulin and glucose [
39], as well as hypertension [
40]. In particular, the latter study demonstrated higher promoter methylation of the
GALNT2 gene, higher levels of ApoB, and lower levels of ApoA1 in hypertensives with the GG genotype [
40].
The
COBLL1 gene encodes a cordon-bleu WH2 repeat protein like 1, having actin monomer and cadherin binding activity [
41]. In the present study, the variant rs12328675 was associated with worse CIMT regression (for the carriers of the C-allele, dominant effect) during 12 months of observation. This association was the most significant among all studied variants over a 12-month period. There were no associations with lipid levels change. The link between this polymorphism and CIMT change can be explained by the ability of the risk allele C, associated with the higher CAD risk, to form transcription factor binding sites for the factors involved in the processes of vascular inflammation regulation and angiogenesis (AP-1 (syn. JUN), SMAD2, SMAD3, SMAD4, and E2F8) [
42].
The
LPA gene encodes lipoprotein (a), which is a well-known and independent risk factor for coronary artery disease [
11]. The studied rs55730499 variant is associated with Lp(a) concentrations [
11], coronary artery disease risk [
8], myocardial infarction risk [
43], and stroke risk [
44]. The T-allele is associated with higher TC levels out of therapy [
43], and in our study, the TT genotype (homozygous for the minor allele) was associated with reduced drug response in terms of TC, LDL-C, and CIMT change. The association with TC change had the highest significance among all associations with TC change, and moreover, rs55730499 was the only variant studied that was associated with both TC, LDL-C, and CIMT change in plasma at the same period of observation (12 months). Based on these findings, the studied
LPA variant can be considered the first of all studied variants to be used as a predictor for rosuvastatin therapy personalization.
The
PPP1R3B gene encodes protein phosphatase 1 regulatory subunit 3B, a key protein in hepatic glycogen metabolism [
12]. The rs9987289 variant is associated with plasma TC and LDL-C levels [
9,
12], but was not tested for carotid atherosclerosis traits or the pharmacogenetics of rosuvastatin. In our study, the presence of the minor A-allele was associated with a worse CIMT response to rosuvastatin. Such an effect could be explained by the association of the studied SNP with inflammatory markers, such as C-reactive protein [
45], known for its influence on cardiovascular risk and atherosclerosis progression [
46], and also by the link of the variant with metabolic syndrome [
47] and glucose levels [
48], taking into account the known role of hyperglycemia in the dysfunction of the endothelium [
49,
50].
The
ABCA1 gene encodes phospholipid-transporting ATPase ABCA1, which provides the efflux of intracellular cholesterol to apolipoproteins and the formation of nascent high-density lipoproteins [
51]. The minor T-allele of the rs1883025 SNP is associated in GWAS with an anti-atherogenic phenotype: lower TC [
9], LDL-C [
52], TG [
53], and HDL-C levels [
52], and with a better response to statin therapy ([
54], the particular drug and lipid are not specified). In the present study, we confirmed a better response to statin therapy in terms of triglyceride reduction in carriers of the minor T-allele (the dominant effect of the SNP).
The
ST3GAL4 gene encodes CMP-N-acetylneuraminate-beta-galactosamide-alpha-2,3-sialyltransferase 4, an enzyme involved in the terminal sialylation of glycoproteins and glycolipids. A beta-galactoside alpha2-3 sialyltransferase takes part in hemostasis (sialylation of plasma von Willebrand factor) and the inflammatory process (selectin-mediated rolling and adhesion of leukocytes during extravasation) [
55]. The studied SNP in this gene, rs11220463, is associated with both blood lipids (TC and LDL-C) [
53] and inflammation (CRP levels) [
56]. In the present study, we have found pharmacogenetic associations of the studied variant with TC and CIMT reduction during therapy. The association with CIMT can be possibly explained by the involvement of the product of
ST3GAL4 in inflammatory processes, taking into account the known association of systemic inflammation (assessed by the systemic immune-inflammatory index) with CIMT [
57]. For rs11220463 the association with inflammation was found in terms of CRP levels [
56], and CRP is known to be associated with CIMT [
58,
59], but this association is not proven to be causal [
58].
The
PSKH1 gene encodes serine/threonine-protein kinase H1, involved in intracellular protein trafficking and pre-mRNA processing [
60]. For the rs16942887 variant, we found associations with LDL-C and CIMT changes on rosuvastatin therapy. There are no reported associations in the literature between this variation and changes in LDL-C and CIMT when taking rosuvastatin. The possible mechanism of the effect of SNP is difficult to predict because of the lack of information on the mechanisms of lipid and vascular influence of
PSKH1.
The
STARD3 gene encodes StAR-related lipid transfer protein 3, which mediates cholesterol transport from the endoplasmic reticulum to endosomes [
61]. In the literature, the rs881844 polymorphism in
STARD3 was known to be associated with plasma lipids, including TC and HDL-C [
43,
53]. In the present study we didn’t find any associations with changes in plasma lipid levels on the rosuvastatin therapy, however, this SNP showed the strongest effects on CIMT changes in both 6- and 12-month periods. The effect on the vascular wall can be explained by the involvement of
STARD3 in the regulation of cholesterol-dependent inflammation and sensitivity to proinflammatory cytokines [
62].
The
PLTP gene encodes phospholipid transfer protein, which is involved in the transfer of phospholipids and free cholesterol from LDLs and VLDLs into HDLs [
63], as well as the uptake of cholesterol from peripheral cells [
64]. Taking into account the function of the gene product, it was rather logical to find the influence of rs6065906 in
PLTP on the lipid-lowering effect of rosuvastatin in terms of TC reduction in the present study, which was not reported before.
Thus, the molecular mechanisms underlying the effects of the studied loci in terms of lipid and CIMT reduction remains unknown, but taking into account the biological functions of the proteins, these genes possess pleiotropic effects on biological and pathological processes including lipid metabolism, inflammation, leukocyte adhesion to the endothelium, endothelial dysfunction, and glucose metabolism. Further experimental studies are needed to explain the functional effects of polymorphisms. Many of the studied loci have been associated with HDL-C levels in GWAS. However, because of non-significant HDL-C change on the therapy, observed in the present study, we didn’t test the SNPs for associations with HDL-C change. This is the limitation of our study, but we can’t say it’s a big omission, because one of the main phenotypes studied here was CIMT, which is known to have a causal relationship with LDL-C and with other factors but not with HDL-C itself [
18].